• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    不同水分環(huán)境下玉米葉形相關(guān)性狀與SSR標記的關(guān)聯(lián)分析*

    2020-11-10 06:46:26趙小強李文麗
    關(guān)鍵詞:葉形葉長自交系

    鐘 源, 趙小強, 李文麗

    不同水分環(huán)境下玉米葉形相關(guān)性狀與SSR標記的關(guān)聯(lián)分析*

    鐘 源, 趙小強**, 李文麗

    (甘肅農(nóng)業(yè)大學甘肅省干旱生境作物學重點實驗室 蘭州 730070)

    玉米的葉形結(jié)構(gòu)與其抗旱性緊密相關(guān), 挖掘不同水分環(huán)境下調(diào)控玉米葉形相關(guān)性狀的顯著關(guān)聯(lián)分子標記位點, 可為揭示玉米葉形結(jié)構(gòu)的分子遺傳機理、克隆相關(guān)調(diào)控基因并進行抗旱理想株型育種提供參考。本研究在不同水分環(huán)境下分析187份玉米自交系葉長、葉寬、葉夾角、葉向值、葉面積、葉形系數(shù)及葉片卷曲度等7個葉形相關(guān)性狀的變化, 采用SSR標記對這些材料進行全基因組掃描并分析其遺傳多樣性, 采用一般線性模型(GLM)尋找不同水分環(huán)境下與玉米7個葉形相關(guān)性狀關(guān)聯(lián)的分子標記。結(jié)果表明: 1)干旱環(huán)境下187份玉米自交系的葉長、葉寬、葉夾角、葉面積均顯著降低, 而葉向值、葉形系數(shù)及葉片卷曲度均明顯升高, 且干旱環(huán)境下這7個葉形相關(guān)性狀的變異比率為30.53%~198.31%。2)145對SSR標記共檢測出652個等位變異, 變異范圍在2~13個; 多態(tài)性信息量(PIC)變異范圍為0.201~0.966, 平均0.478。采用UPGMA聚類及群體遺傳結(jié)構(gòu)分析均將供試材料分為旅大紅骨(LRC)、唐四平頭(TSPT)、蘭卡斯特(Lan)、P及瑞德(Reid)等5大類群。3)采用GLM模型, 在不同水分環(huán)境下共檢測到15個SSR標記與玉米的7個葉形相關(guān)性狀在<0.01水平下關(guān)聯(lián), 各標記對表型的解釋率為2.25%~27.30%, 72.97%的標記可在干旱環(huán)境下被檢測到; 其中umc1124、umc2363、umc2214、umc1742、phi331888、umc1378、bnlg1863、umc2134和umc1345標記同時與不同水分環(huán)境下的多個葉形相關(guān)性狀連鎖, 表現(xiàn)出明顯的“一因多效”現(xiàn)象。這些研究結(jié)果將為玉米葉形結(jié)構(gòu)的遺傳改良及抗旱理想株型分子標記輔助選擇育種提供參考。

    玉米; 葉形結(jié)構(gòu); 干旱; 遺傳結(jié)構(gòu); 關(guān)聯(lián)分析

    玉米()是集糧食、飼料及生物質(zhì)能源為一體的多用途C4作物, 其在穩(wěn)定國家糧食安全、促進畜牧業(yè)發(fā)展及改善人民生活水平等方面扮演著重要角色。干旱是限制玉米增產(chǎn)、穩(wěn)產(chǎn)最為重要的非生物逆境脅迫因素之一。Zhao等[1]報道, 近30 a來干旱導致全球玉米減產(chǎn)39.3%。因此, 提高玉米的抗旱性并培育優(yōu)良抗旱新品種是解決玉米安全生產(chǎn)的有效育種策略。葉片不僅是玉米進行光合作用的重要組織源, 而且是玉米感知干旱、溫度等外界壓力信號及與環(huán)境進行水分、氣體交換的重要器官[2-5], 與玉米的產(chǎn)量和品質(zhì)緊密相關(guān)。葉長、葉寬、葉夾角、葉向值、葉面積、葉形系數(shù)及葉片卷曲度等玉米葉形結(jié)構(gòu)是屬于多個微效基因共同調(diào)控的典型數(shù)量遺傳性狀, 受環(huán)境效應的影響較大[6]。Zhao等[1]、Zhao等[2]、Ribaut等[7]、Bennetzen等[8]和趙小強等[9]研究揭示, 干旱缺水時不同株型玉米的葉形結(jié)構(gòu)、氣孔形態(tài)、葉片的各種生理代謝及光合特性均會發(fā)生適應性的改變, 其葉長、葉寬、葉夾角及葉面積均顯著減小, 葉形系數(shù)和葉向值均會顯著增大; 葉片上/下表皮氣孔個數(shù)會增多, 氣孔長度和寬度均會變小; 葉片的葉綠素a、b、a+b含量及葉綠素a/b會顯著降低, 葉片的氣孔導度、凈光合速率、葉片溫度、蒸騰速率和RuBPCase活性均會明顯減小, 胞間CO2濃度會明顯升高, 玉米的這些葉形結(jié)構(gòu)的適應性變化共同構(gòu)成了玉米的株型抗旱機制。

    近年來, 數(shù)量遺傳學家先后開發(fā)的加性-顯性-母體遺傳模型(additive-dominant-maternal effect genetic model, ADM)[3]、混合線性加顯模型(mixed additive-dominance model, AD)[10]、主基因+多基因遺傳模型(major genes plus polygenes genetic model, MG+PG)[11]、ABC法[10]、Geoffrey分析法[12]及世代分析法[13]等遺傳模型與分析方法在玉米葉形相關(guān)性狀遺傳機理研究中得到了廣泛應用, 為玉米葉形結(jié)構(gòu)的遺傳改良及理想株型育種奠定了堅實的基礎(chǔ)。加之借助分子標記技術(shù)鑒定不同環(huán)境下調(diào)控玉米葉形相關(guān)性狀的基因/QTL (quantitative trait locus)位點, 為深入解析其葉形結(jié)構(gòu)的遺傳機制提供了可行的手段。Mickelson等[14]利用B73×Mo17發(fā)展了1套RILs (recombinant inbred lines)群體, 通過CIM (composite interval mapping)法在不同生態(tài)環(huán)境下定位到了10個QTLs與葉夾角相關(guān)。Pan等[15]以20份葉形結(jié)構(gòu)差異較大的自交系構(gòu)建1套含有1 887份家系的ROAM (random-open-parent association mapping)群體, 對玉米的葉長、葉寬、葉夾角、穗上葉片數(shù)、穗下葉片數(shù)及總?cè)~數(shù)等進行了QTL、QTL×E (QTL by environment interaction)及GWAS (genome-wide association study)分析。Ku等[5]通過Yu82×Shen137構(gòu)建了1套F2∶3群體, 采用CIM法在兩種生態(tài)環(huán)境下檢測到了15個調(diào)控玉米葉長、葉寬、葉夾角及葉向值的QTLs位點。Zhao等[2]利用Langhuang×TS141和Chang7-2×TS141構(gòu)建了2套F2和2套F2∶3群體, 采用CIM法在9種不同水分環(huán)境下共識別到了29個調(diào)控葉夾角、葉向值、葉長、葉寬、葉面積及葉形系數(shù)的sQTLs (stable QTLs), 并通過meta分析, 檢測到了21個meta-QTLs (mQTLs)和24個調(diào)控玉米葉形結(jié)構(gòu)發(fā)育的候選基因。甚至, Ku等[4]、郭書磊等[6]和Wang等[16]先后整合了多套不同遺傳背景下檢測到的玉米葉長、葉寬、葉夾角及葉向值QTLs信息, 通過meta分析分別識別到了17、73和12個相關(guān)mQTLs, 并預測到了15、44和7個與玉米葉片發(fā)育緊密關(guān)聯(lián)的候選基因。從這些研究不難看出, 通過構(gòu)建定位群體雖然可以找到一些重要QTLs/基因信息, 但是該方法不僅需要花費大量人力和時間, 且由于構(gòu)建的定位群體所采用材料的遺傳背景、群體結(jié)構(gòu)與大小及QTL定位方法不同, 導致QTL定位結(jié)果不盡相同, 嚴重影響了相關(guān)性狀的基因克隆及分子標記輔助選擇(marker-assisted selection, MAS)育種的應用。

    基于連鎖不平衡作圖(linkage disequilibrium, LD)的關(guān)聯(lián)分析法(association analysis), 是通過對供試群體候選基因的分子標記進行掃描或檢測, 進而確定相應群體不同遺傳材料間所攜帶的等位基因位點信息并分析等位基因?qū)δ繕诵誀钬暙I大小的一種方法[17]。該方法作為數(shù)量遺傳性狀研究中的重要補充方法, 已被廣泛應用到玉米[18]、小麥()[19]及水稻()[20]等作物重要數(shù)量遺傳性狀等位基因位點挖掘中?;谝陨峡紤], 本研究采用SSR (simple sequence repeats)標記技術(shù)揭示了187份不同玉米自交系的遺傳多樣性, 在不同水分環(huán)境處理下通過對這些自交系的葉長、葉寬、葉夾角、葉向值、葉面積、葉形系數(shù)及葉片卷曲度等葉形相關(guān)性狀進行了田間鑒定, 并采用關(guān)聯(lián)分析法定位不同水分環(huán)境下與玉米葉形相關(guān)性狀緊密連鎖的SSR標記, 以期為玉米葉形相關(guān)抗旱優(yōu)良新品種的MAS育種提供參考。

    1 材料與方法

    1.1 試驗材料

    以黃早四(來源于四平頭雜株)、B73(來源于BSSSC5)、丹340(來源于旅9×有稃玉米的雜交)、齊319(選自美國雜交種P78799)、Mo17(來源于C103×C1187-2)、掖478(來源于沈5003×U8112)等6份玉米雜種優(yōu)勢劃分標準測驗種在內(nèi)的84份國內(nèi)核心玉米種質(zhì)及本試驗自育的103份玉米自交系為試材。這些材料均由甘肅省干旱生境作物學重點實驗室-分子標記輔助選育實驗室提供。

    1.2 試驗方法

    1.2.1 田間玉米自交系葉形結(jié)構(gòu)鑒定

    試驗于2017年4月中旬至9月下旬在甘肅省農(nóng)業(yè)科學院黃羊麥類作物育種試驗站(37.97oN, 102.63oE, 海拔1 508 m)進行。試驗均設(shè)正常供水(CK, 玉米全生育期內(nèi)只要降水不足就及時供水)及干旱脅迫處理(DS, 玉米大喇叭口前期至花期結(jié)束不供水, 其他生育時期每隔20 d供水1次)[21-23]。每一處理均采用完全隨機區(qū)組設(shè)計, 3次重復。每一自交系設(shè)單行區(qū), 行長5 m, 株距25 cm, 行距40 cm。播前于2017年3月統(tǒng)一整地、施肥并采用平膜覆蓋試驗地, 其他管理同一般大田。2017年該試驗地玉米生育期(4—9月)內(nèi)月平均氣溫為17.8 ℃, 總降水量為169.6 mm?;ê蟮?0 d從每一處理相應自交系中選擇長勢整體一致的10株, 并參照Zhao等[2]和趙小強等[24]方法測定其棒三葉的葉長(leaf length, LL; 葉基至葉頂端的長度)、葉寬(leaf width, LW; 葉片中間展開寬度)、葉夾角(leaf angle, LA; 葉片與莖稈的銳角角度), 并計算其葉面積(leaf size, LS)、葉向值(leaf orientation value, LOV)、葉形系數(shù)(leaf shape coefficient, LSC)及葉片卷曲度(leaf rolling index, LRI), 計算公式為:

    [26](2)

    [2](3)

    式中: LF (length from beginning of ligule to flagging point of leaf, 葉基至葉最高點的長度)為高點長; LCW (width of a leaf in the middle of natural curl)為葉片中間自然卷曲狀態(tài)下的寬度;為3, 代表測定葉片數(shù)。

    1.2.2 SSR標記分析

    采用CTAB法[28]提取從田間剪取的五葉期玉米幼苗的全基因組DNA, 并用1%的瓊脂糖凝膠電泳檢測DNA質(zhì)量, 紫外分光光度計檢測DNA濃度并稀釋至50 ng?μL?1, ?20 ℃低溫保存?zhèn)湓囼炈琛腗aizeGDB (Maize Genome Database)數(shù)據(jù)庫網(wǎng)站(http://www.maizegdb.org/)中選取均勻分布于玉米10個連鎖群多態(tài)性較好的SSR標記145對, 并參照鐘源等[28]的方法進行PCR擴增、聚丙烯酰胺凝膠電泳及銀染。

    1.3 數(shù)據(jù)統(tǒng)計分析

    田間不同水分環(huán)境處理下測定的187份玉米自交系的葉形相關(guān)性狀數(shù)據(jù)均采用統(tǒng)計軟件IBM SPSS 19.0 (SPSS Inc., Chicago, IL, USA)進行描述性統(tǒng)計分析、聯(lián)合方差分析及Pearson相關(guān)分析。SSR擴增產(chǎn)物構(gòu)建0、1矩陣, 并以A、B、C等不同大寫字母記錄基因型。SSR位點的多態(tài)性信息量(polymorphism information content value, PIC)按以下公式計算:

    2 結(jié)果與分析

    2.1 不同水分環(huán)境下參試玉米自交系葉形相關(guān)性狀的表型鑒定

    在正常供水下, 187份玉米自交系的7個葉形相關(guān)性狀的變異系數(shù)介于13.21%(葉寬)~59.23%(葉片卷曲度); 而干旱脅迫處理下, 這些玉米自交系的7個葉形相關(guān)性狀的變異系數(shù)介于9.62%(葉長)~ 28.22%(葉夾角)(表1)。說明玉米自交系的7個葉形相關(guān)性狀在不同水分處理下均存在豐富的遺傳變異。方差分析表明, 187份玉米自交系的7個葉形相關(guān)性狀在自交系間均達到<0.01或<0.05水平差異顯著(表2)。說明玉米的這些葉形結(jié)構(gòu)受自交系自身遺傳特性的調(diào)控, 因此, 育種中可以重點改良玉米的這些葉形結(jié)構(gòu), 進而達到改良玉米株型的目的。這7個葉形相關(guān)性狀在水分處理間及自交系與水分處理互作間差異均極顯著(<0.01)或顯著(<0.05), 且在干旱脅迫處理下, 這些玉米自交系的葉長、葉寬、葉夾角、葉面積均顯著降低, 而葉向值、葉形系數(shù)及葉片卷曲度均顯著增大, 且這7個葉形相關(guān)性狀的變異比率介于30.53%(葉片卷曲度)~198.31%(葉夾角)(表1和表2)。表明玉米為了避免水分環(huán)境對其造成的傷害, 植株的葉長、葉寬、葉夾角、葉向值、葉面積、葉形系數(shù)及葉片卷曲度等葉形結(jié)構(gòu)明顯地發(fā)生適應性的改變, 且這些性狀還顯著受基因型與環(huán)境互作效應的共同調(diào)控。此外, 在兩種水分環(huán)境下, 187份玉米自交系7個葉形相關(guān)性狀的偏度和峰度均基本介于?1.0~1.0(表1)。表明玉米的這些葉形相關(guān)性狀呈典型的數(shù)量遺傳特性, 因此可以對其進行QTL定位或關(guān)聯(lián)分析。

    2.2 參試玉米自交系葉形相關(guān)性狀的相關(guān)性分析

    187份玉米自交系的7個葉形相關(guān)性狀間Pearson相關(guān)分析結(jié)果(表3)表明, 葉長與葉面積、葉形系數(shù)及葉寬間均極顯著正相關(guān)(<0.01), 與葉夾角、葉向值及葉片卷曲度間均顯著負相關(guān)(<0.05); 葉寬與葉面積間極顯著正相關(guān), 與葉形系數(shù)、葉片卷曲度、葉夾角及葉向值間極顯著(<0.01)或顯著(<0.05)負相關(guān); 葉夾角與葉向值、葉面積及葉形系數(shù)間均極顯著(<0.01)或顯著(<0.05)負相關(guān); 葉向值與葉面積間極顯著負相關(guān)(<0.01); 葉面積與葉形系數(shù)和葉片卷曲度間顯著負相關(guān); 葉形系數(shù)與葉片卷曲度間顯著(<0.05)正相關(guān)(表3)。說明玉米的這些葉形相關(guān)性狀間彼此關(guān)聯(lián), 其相互協(xié)同或抑制作用共同建成了不同玉米自交系特定的葉形結(jié)構(gòu)。

    2.3 參試玉米自交系的SSR標記遺傳多樣性分析

    選擇均勻分布于玉米10個連鎖群的145對SSR標記對187份玉米自交系進行了全基因組掃描, 共檢測出了652個等位基因, 每對SSR標記檢測出了2~13個, 平均為4.50個, 其中位于Bin 3.05處的標記phi073的等位基因最多檢測出了13個, 而等位基因變異最少的是位于Bin 1.05處的umc2025等13個標記, 均為2個等位基因(表4)。每對SSR標記的多態(tài)性信息量(PIC)的變異范圍介于0.200~0.966, 平均為0.478。其中標記umc1278 (Bin 1.07)的PIC最大, 而標記bnlg1877 (Bin 10.07)的PIC最小(表4)。進一步利用遺傳相似系數(shù)(GS)和UPGMA分析, 參試所有玉米自交系的GS變異范圍為0.610~1.000, 并在GS在0.612時可將187份玉米自交系明顯劃分為5大類群, 即包含丹340等25份材料的旅大紅骨類群(LRC)、黃早四為代表的唐四平頭類群(TSPT)為34份材料、Mo17在內(nèi)的28份蘭卡斯特類群(Lan)材料、包括齊319等38份材料的P群(P)及B73與掖478等62份材料聚成的瑞德類群(Reid)材料(表5)。

    表1 不同水分環(huán)境下187份玉米自交系葉形相關(guān)性狀的表型值

    LL: 葉長; LW: 葉寬; LA: 葉夾角; LOV: 葉向值; LS: 葉面積; LSC: 葉形系數(shù); LRI: 葉片卷曲度; CK: 正常供水處理; DS: 干旱脅迫處理。LL: leaf length; LW: leaf width; LA: leaf angle; LOV: leaf orientation value; LS: leaf size; LSC: leaf shape coefficient; LRI: leaf rolling index; CK: normal water treatment; DS: drought stress treatment.

    表2 不同水分環(huán)境下187份玉米自交系葉形相關(guān)性狀的方差分析(F值)

    *或**分別代表<0.05或<0.01水平差異顯著。LL: 葉長; LW: 葉寬; LA: 葉夾角; LOV: 葉向值; LS: 葉面積; LSC: 葉形系數(shù); LRI: 葉片卷曲度。* or ** indicate significant difference at<0.05 or<0.01 levels, respectively. LL: leaf length; LW: leaf width; LA: leaf angle; LOV: leaf orientation value; LS: leaf size; LSC: leaf shape coefficient; LRI: leaf rolling index.

    表3 187份玉米自交系葉形相關(guān)性狀的相關(guān)分析

    *或**分別代表<0.05或<0.01水平顯著相關(guān)。LL: 葉長; LW: 葉寬; LA: 葉夾角; LOV: 葉向值; LS: 葉面積; LSC: 葉形系數(shù); LRI: 葉片卷曲度。* or ** indicate significant correlation at<0.05 or<0.01 levels, respectively. LL: leaf length; LW: leaf width; LA: leaf angle; LOV: leaf orientation value; LS: leaf size; LSC: leaf shape coefficient; LRI: leaf rolling index.

    表4 145對SSR標記在187份玉米自交系中檢測到的等位基因數(shù)及多態(tài)性信息量(PIC)

    續(xù)表4

    ***為不確定堿基; —為缺失。*** means uncertain base; — means missing.

    2.4 參試玉米自交系群體的遺傳結(jié)構(gòu)分析

    利用Structure 2.3.1軟件對187份玉米自交系進行群體遺傳結(jié)構(gòu)和模型聚類分析。結(jié)果表明當=1~10時, Δ在=5時出現(xiàn)明顯峰值(圖1), 因此可以判定187份玉米自交系可劃分為5大類群。5大類群的遺傳固定指數(shù)均>0.05。說明這5大類群均能獨立成群。對照代表性自交系可知, 這5大類群分別為LRC(紅色)、TSPT(綠色)、Lan(粉色)、P(黃色)及Reid(藍色)類群等(圖2)。這與UPGMA聚類分析大部分材料的分群結(jié)果基本一致。

    2.5 不同水分環(huán)境下玉米葉形相關(guān)性狀與SSR標記的關(guān)聯(lián)分析

    采用一般線性模型(GLM)對SSR標記位點與不同水分環(huán)境下玉米7個葉形相關(guān)性狀進行了關(guān)聯(lián)分析。結(jié)果(表6)表明, 共檢測到15對SSR標記與不同水分環(huán)境下玉米的葉長、葉寬、葉夾角、葉向值、葉面積、葉形系數(shù)及葉片卷曲度等葉形相關(guān)性狀在<0.01水平顯著關(guān)聯(lián)。表型變異的解釋率介于2.25%(位于Bin 10.03處的標記umc1345與干旱脅迫下的葉形系數(shù)相關(guān))~27.30%(位于Bin 1.01處的bnlg149位點與正常供水環(huán)境下的葉夾角連鎖), 其主要位于第1、2、3、4、5、7、8、9及10號染色體上。除bnlg149 (Bin 1.01)與葉夾角、umc2376 (Bin 3.01)與葉片卷曲度、bnlg1621a (Bin 4.06)與葉寬、umc2356 (Bin 8.05)與葉面積、umc2218 (Bin 8.08)與葉寬及bnlg1451 (Bin 10.01)與葉向值等性狀相關(guān)外, 其余9對SSR標記均與2~7個葉形相關(guān)性狀間同時關(guān)聯(lián)。

    3 討論

    3.1 干旱脅迫對玉米葉形結(jié)構(gòu)的影響

    葉片是玉米進行生命活動的基礎(chǔ), 其主要參與玉米的光合作用、蒸騰作用及玉米對各種逆境脅迫的應答響應中, 最終影響玉米的產(chǎn)量和品質(zhì)。Zhao等[2]以抗旱性差異較強的3份玉米自交系及其構(gòu)建的2個F1雜交種、2套F2及2套F2∶3群體在9種不同水分環(huán)境下分析發(fā)現(xiàn), 干旱脅迫時這些材料的葉長(2.8%~13.4%)、葉寬(11.0%~26.2%)、葉夾角(13.4%~26.3%)、葉面積(18.3%~36.3%)均顯著降低, 而葉向值(4.2%~29.1%)和葉形系數(shù)(3.2%~25.1%)均顯著升高。張振平等[29]通過對7份抗旱能力不同的玉米品種在中度水分脅迫下研究表明, 玉米的葉片卷曲度和葉片衰老指數(shù)與其抗旱性間均呈顯著負相關(guān), 葉向值和葉形系數(shù)與其抗旱性間達到了顯著正相關(guān); 進一步灰色關(guān)聯(lián)度分析揭示, 與葉片卷曲度和葉片衰老指數(shù)相比, 葉形系數(shù)和葉向值在玉米抗旱性中發(fā)揮更為重要的作用。趙小強等[9]研究闡明, 干旱脅迫下玉米通過葉形結(jié)構(gòu)發(fā)生的適應性改變, 調(diào)節(jié)影響葉片的氣孔形態(tài)、各種生理代謝及光合特性, 以減少植株對太陽輻射的接收、降低植株的葉片溫度及蒸騰散失, 提高植株的水分利用效率、改善植株的生理代謝, 最終達到改善玉米抗旱性的目的。本研究也得到了高度一致的結(jié)果, 即在花期干旱脅迫處理下, 187份玉米自交系的7個葉形相關(guān)性狀在水分處理間差異均達極顯著水平。說明不同水分環(huán)境均會顯著影響玉米葉形結(jié)構(gòu)的塑造, 進而影響后期玉米籽粒的發(fā)育及產(chǎn)量的形成。進一步分析發(fā)現(xiàn), 與對照正常供水環(huán)境相比, 干旱脅迫后187份玉米自交系的葉長、葉寬、葉夾角、葉面積平均降低了6.56%、7.16%、7.32%和7.97%, 而葉向值、葉形系數(shù)及葉片卷曲度平均升高了4.34%、6.16%和32.16%, 且干旱脅迫下這7個葉形相關(guān)性狀的變異比率為30.53%~198.31%。這與趙小強等[9]在玉米、Yang等[27]和Isidro等[30]在小麥、Sohrabi等[31]在水稻、張飛等[32]在高粱()、朱宗河等[33]在油菜()及Masarirambi等[34]在馬鈴薯()上的研究相類似。因此, 植物的葉形結(jié)構(gòu)與其抗旱性緊密相關(guān), 特別是在玉米上, 以上7個葉形相關(guān)性狀均可作為玉米抗旱性鑒定的重要參考指標。

    表5 187份玉米自交系UPGMA聚類分析

    字體加粗的材料為相應雜種優(yōu)勢類群的標準測驗種。The bolded materials are the standard test species for the corresponding heterosis groups.

    圖1 187份玉米自交系的群體遺傳結(jié)構(gòu)分析的ΔK值隨K值變化曲線

    : 群體個數(shù)The number of population.

    圖2 基于SSR標記的187份玉米自交系群體遺傳結(jié)構(gòu)

    LRC: 旅大紅骨類群; TSPT: 唐四平頭類群; Lan: 蘭卡斯特類群; P: P群; Reid: 瑞德類群; 不同顏色代表不同類群, 顏色多少表示該類群的概率。LRC: luda red cob group; TSPT: tangsipingtou group; Lan: lancaster group; P: P group; Reid: Reid group; Different colors indicate the corresponding groups, areas of colors indicate probability of the corresponding groups.

    表6 不同水分環(huán)境下與玉米葉形性狀相關(guān)的SSR標記及其對表型變異的解釋率(P<0.01)

    LL: 葉長; LW: 葉寬; LA: 葉夾角; LOV: 葉向值; LS: 葉面積; LSC: 葉形系數(shù); LRI: 葉片卷曲度。CK: 正常供水處理; DS: 干旱脅迫處理。LL: leaf length; LW: leaf width; LA: leaf angle; LOV: leaf orientation value; LS: leaf size; LSC: leaf shape coefficient; LRI: leaf rolling index. CK: normal water treatment; DS: drought stress treatment.

    3.2 玉米自交系遺傳多樣性及其葉形相關(guān)性狀間的變異分析

    玉米自交系彼此間的遺傳基礎(chǔ)狹窄一直是困擾重大突破性玉米優(yōu)良新品種選育的關(guān)鍵因素之一。因此, 從分子水平及表型上全面剖析玉米自交系間的遺傳多樣性, 并比較親本自交系間的親緣遠近, 對于有目的地組配優(yōu)良雜交組合及選育出優(yōu)良新品種具有十分重要的指導意義。目前, 國內(nèi)外學者已運用SSR標記技術(shù)對不同來源玉米自交系的遺傳多樣性及雜種優(yōu)勢類群進行了廣泛研究, 認為SSR標記的遺傳距離與雜種優(yōu)勢的表現(xiàn)呈顯著正相關(guān), 但由于自交系的類型及來源不同導致其雜種優(yōu)勢類群劃分結(jié)果不盡一致[28,35]。本研究采用均勻分布于玉米10個連鎖群的145對SSR標記對國內(nèi)84份核心自交系及本課題組近年來自育的103份自交系總計187份材料進行了遺傳多樣性和遺傳結(jié)構(gòu)分析, 結(jié)果表明, 187份玉米自交系的GS介于0.610~1.000, 并采用UPGMA聚類法當GS在0.612時可將187份玉米自交系明顯劃分為LRC、TSPT、Lan、P及Reid等5大類群, 且這與本研究采用Structure 2.3.1軟件群體遺傳結(jié)構(gòu)分析結(jié)果基本相一致。Wu等[36]采用MaizeSNP50 BeadChip芯片(包含56 110個SNP標記)也將367份中國玉米育種自交系劃分成了LRC、TSPT、Lan、P及Reid等類群。此外, 本研究分析還表明, 187份玉米自交系在2種水分環(huán)境下7個葉形相關(guān)性狀的變異系數(shù)介于9.62%~59.23%。說明玉米這187份自交系的遺傳多樣性較大, 在7個葉形相關(guān)性狀間均存在豐富的遺傳變異。

    3.3 玉米自交系葉形相關(guān)性狀與SSR標記的關(guān)聯(lián)分析

    Ingvarsson等[37]指出關(guān)聯(lián)分析可根據(jù)作物自然群體(自交系)長期演化過程中積累的重組信息, 進而檢測出群體內(nèi)處于非隨機狀態(tài)的標記或候選基因的變異與目標性狀顯著性關(guān)聯(lián)的頻率, 并能對目標性狀基因進行精細定位。與QTL定位相比, 該方法可有效避免作圖群體構(gòu)建過程中費時費力等缺點。因此, 本研究采用GLM模型對2種水分環(huán)境下187份玉米自交系的7個葉形相關(guān)性狀進行了關(guān)聯(lián)分析, 在葉長(5對)、葉寬(6對)、葉夾角(5對)、葉向值(5對)、葉面積(7對)、葉形系數(shù)(5對)及葉片卷曲度(4對)間檢測到了15對關(guān)聯(lián)的SSR標記, 其主要分布于除第6染色體以外的其余9號染色體上, 單個關(guān)聯(lián)SSR標記的表型變異貢獻率為2.25%~27.30%。在所有找到的關(guān)聯(lián)SSR標記中, 大約72.97%的關(guān)聯(lián)SSR標記可在干旱脅迫環(huán)境下被檢測到。因此, 干旱脅迫下檢測到的與玉米7個葉形相關(guān)性狀相關(guān)的SSR標記為今后玉米葉形相關(guān)抗旱分子遺傳機理的解析及抗旱MAS育種奠定基礎(chǔ)。此外, Tuberosa等[38]認為由于不同性狀間存在共同的遺傳機制, 因此“一因多效”或是控制不同性狀的基因普遍呈緊密連鎖的現(xiàn)象, 而在本研究中總共檢測到了9個SSR標記(umc1124、umc2363、umc2214、umc1742、phi331888、umc1378、bnlg1863、umc2134和umc1345)同時與多個玉米葉形相關(guān)性狀連鎖, 且這與本研究7個葉形相關(guān)性狀間的Pearson相關(guān)分析的結(jié)果高度一致, 即單個葉形相關(guān)性狀與其余4~6個葉形相關(guān)性狀呈顯著或極顯著正/負相關(guān)。說明玉米染色體的這些SSR標記位點共同調(diào)控著玉米的7個葉形相關(guān)性狀的遺傳。

    其中, 位于Bin 2.02處的umc2363同時與2種水分環(huán)境下的葉長、葉夾角、葉片卷曲度及正常供水環(huán)境下的葉面積間相關(guān), 解釋的表型變異為2.42%~4.11%。趙小強[10]采用2套F2∶3群體也在umc2363附近找到了多個干旱脅迫環(huán)境下調(diào)控葉夾角的主效QTL, 其受基因加性遺傳效應的控制。于多個正常供水和干旱脅迫下, Zhao等[2]也在Bin 2.02的umc2363-umc2403區(qū)間也檢測到了1個調(diào)控葉夾角的sQTL, 其可解釋6.81%~19.53%的葉夾角變異。進一步采用mQTL分析, Zhao等[2]還在該區(qū)域的umc2363-umc1227區(qū)域內(nèi)檢測到了1個同時調(diào)控葉夾角、葉向值及葉寬的mQTL2-1和1個調(diào)控玉米葉耳和葉舌發(fā)育及葉片直立生長的()候選基因。標記umc2214(Bin 2.09)同時與2種水分環(huán)境下的葉長、葉寬、葉面積及葉形系數(shù)間相連鎖, 表型變異的解釋率為2.42%~3.41%。Zhao等[2]也在umc2184-umc2214區(qū)間內(nèi)挖掘到了1個調(diào)控玉米葉長的mQTL2-3。Li等[39]還在此區(qū)域的233.55~236.61 Mb處預測到了1個調(diào)控葉長的()候選基因; 甚至, 采用圖位克隆的方法, Bolduc等[40]證實了調(diào)控玉米葉片的形成。phi331888 (Bin 5.04)位點同時與2種水分環(huán)境下的葉長、葉寬、葉夾角、葉向值及正常供水環(huán)境下的葉面積、葉形系數(shù)、葉片卷曲度相連鎖, 解釋的表型變異為2.39%~4.52%。Ku等[5]也在Bin 5.04的bnlg1287附近同時定位到了1個同時調(diào)控葉夾角和葉向值的QTL, 受基因部分顯性遺傳效應的調(diào)控, 表型貢獻率分別為9.7%和9.8%。同時劉鵬飛等[41]在Bin 5.04處的umc1171附近挖掘到1個調(diào)控葉夾角的QTL, 其受基因加性遺傳效應的調(diào)控, 表型貢獻率為7.4%。在Bin 7.00 (umc1378)處同時與2種水分環(huán)境下的葉長、葉夾角、葉向值、葉面積及正常供水環(huán)境下的葉寬相關(guān), 解釋的表型變異為2.53%~4.20%。劉鵬飛等[41]在umc1378-umc1545區(qū)間檢測到了1個調(diào)控葉夾角的QTL, 其受基因部分顯性遺傳效應的調(diào)控, 可解釋3.3%的葉夾角變異。在umc2177-umc378-umc1426 (Bin 7.00)區(qū)間Zhao等[2]在4套F2和F2∶3群體的6種不同水分環(huán)境下同時檢測到1個調(diào)控葉夾角和葉向值的sQTL。甚至, 趙小強[10]對umc1378-umc2177區(qū)間同時調(diào)控玉米葉夾角和葉向值的主效QTL進行精細定位和親本間Illumina HiSeq 2500重測序后, 最終在該主效QTL區(qū)間挖掘到了1個調(diào)控玉米葉片發(fā)育的()基因; 有趣的是, 郭書磊等[6]也在此區(qū)域檢測到了1個調(diào)控葉寬的mQTL, 并找到了和候選基因, 其中是1個的基因并對葉片基部-末端軸向的形成中起調(diào)控作用[42], 而是1個基因, 該基因能夠上調(diào)分生組織中家族基因的表達, 進而保持葉片細胞分裂活性和抑制葉片中和家族基因, 進而促進葉片近軸組織的分化發(fā)育[43]。bnlg1863 (Bin 8.03)位點同時與2種水分環(huán)境下的葉面積及干旱脅迫下的葉片卷曲度相關(guān)聯(lián), 解釋的表型變異為9.44%~22.17%。在Bin 8.03~8.04的umc1904-umc1460區(qū)域劉鵬飛等[41]在1套F2群體上識別到了1個調(diào)控葉夾角的QTL, 其受基因加性遺傳效應的調(diào)控, 解釋的表型貢獻率為7.1%。趙小強[10]在Bin 8.03處的umc1530-bnlg1863、bnlg1863-umc2075、umc2075-bnlg162區(qū)間分別檢測到多個水分環(huán)境下同時調(diào)控玉米葉長、葉面積及葉向值的QTL/sQTL。甚至, Zhao等[2]還在bnlg1863-umc1302區(qū)間上檢測到了1個同時調(diào)控葉夾角、葉長、葉面積及葉形系數(shù)的mQTL, 并挖掘到了1個負向調(diào)控葉原基發(fā)育的轉(zhuǎn)錄因子()[44]。umc2134 (Bin 9.05)標記同時與2種水分環(huán)境下的葉形系數(shù)及正常供水下的葉向值相關(guān)聯(lián), 解釋的表型變異為3.18%~7.86%。Zhao等[2]在Bin 9.05 (bnlg1270a-umc1231)處識別到了1個同時調(diào)控葉寬、葉向值及葉面積的mQTL, 并挖掘出了()基因, 該基因作為一種ADP/ATP載體, 在玉米葉基部分生組織發(fā)育中積累。表明Bin 2.02、Bin 2.09、Bin 5.04、Bin 7.00、Bin 8.03及Bin 9.05等附近是調(diào)控玉米多個葉形相關(guān)性狀的重要“一因多效”或富集位點, 且這些區(qū)域存在一些重要調(diào)控基因, 為不同水分環(huán)境下玉米葉片結(jié)構(gòu)的塑造發(fā)揮重要的積極效應。此外, 本研究還找到了3個新的SSR標記位點與不同水分環(huán)境下的玉米的葉形相關(guān)性狀在<0.01水平下關(guān)聯(lián), 分別是位于Bin 1.05 (umc1124)處同時與2種水分環(huán)境下的葉夾角及正常供水環(huán)境下的葉向值, 位于umc1742 (Bin 3.04)處同時與正常供水環(huán)境下的葉長、葉寬及葉形系數(shù), 及位于umc1345 (Bin 10.03)位點與2種水分環(huán)境下的葉面積和葉形系數(shù)間相連鎖SSR標記, 這3個連鎖標記可為不同水分環(huán)境下調(diào)控玉米葉形相關(guān)性狀QTL/基因挖掘提供新信息。

    4 結(jié)論

    干旱脅迫環(huán)境下不同玉米自交系的葉長、葉寬、葉夾角和葉面積顯著降低, 而葉向值、葉形系數(shù)和葉片卷曲度明顯升高, 且干旱脅迫后玉米自交系的葉長(52.60%)、葉寬(76.46%)、葉夾角(198.31%)、葉向值(118.62%)、葉面積(65.27%)、葉形系數(shù)(66.57%)及葉片卷曲度(30.53%)間的變異比率間差異顯著, 這7個葉形性狀與玉米的抗旱性緊密相關(guān), 可作為玉米抗旱性鑒定的重要參考指標。采用均勻覆蓋玉米全基因組的145對SSR標記可將供試玉米自交系劃分為LRC、TSPT、Lan、P及Reid等5大類群, 其遺傳多樣性較大, PIC介于0.200~0.966。進一步通過GLM分析在2種水分環(huán)境下總共檢測到了15個SSR標記與玉米的7個葉形性狀相關(guān)聯(lián), 其中9個標記同時與2種水分環(huán)境下的多個葉形相關(guān)性狀連鎖, 表現(xiàn)出明顯的“一因多效”現(xiàn)象。

    [1] ZHAO X Q, ZHANG J W, FANG P, et al. Comparative QTL analysis for yield components and morphological traits in maize (L.) under water-stressed and well-watered conditions[J]. Breeding Science, 2019, 69(4): 621–632

    [2] ZHAO X Q, FANG P, ZHANG J W, et al. QTL mapping for six ear leaf architecture traits under water-stressed and well-watered conditions in maize (L.)[J]. Plant Breeding, 2018, 137(1): 60–72

    [3] 趙小強, 姬祥卓, 彭云玲. 花期玉米()穗三葉葉形結(jié)構(gòu)的遺傳效應及配合力分析[J]. 分子植物育種, 2020, 18(9): 3024–3031 ZHAO X Q, JI X Z, PENG Y L. Analysis of combining ability and genetic effect for leaf shape about the three-ear-leaves at flowering stage in[J]. Molecular Plant Breeding, 2020, 18(9): 3024–3031

    [4] KU L X, ZHANG J, GUO S L, et al. Integrated multiple population analysis of leaf architecture traits in maize (L.)[J]. Journal of Experimental Botany, 2012, 63(1): 261–274

    [5] KU L X, ZHAO W M, ZHANG J, et al. Quantitative trait loci mapping of leaf angle and leaf orientation value in maize (L.)[J]. Theoretical and Applied Genetics, 2010, 121(5): 951–959

    [6] 郭書磊, 張君, 齊建雙, 等. 玉米葉形相關(guān)性狀的Meta-QTL及候選基因分析[J]. 植物學報, 2018, 53(4): 487–501 GUO S L, ZHANG J, QI J S, et al. Analysis of Meta- quantitative trait loci and their candidate genes related to leaf shape in maize[J]. Chinese Bulletin of Botany, 2018, 53(4): 487–501

    [7] RIBAUT J M, BETRAN J, MONNEVEUX P, et al. Drought tolerance in maize[M]//BENNETZEN J L, HAKE S C. Handbook of Maize: Its Biology. New York: Springer, 2009: 311–344

    [8] BENNETZEN J L, HAKE S C. Handbook of Maize: Its Biology[M]. New York: Springer, 2010: 324

    [9] 趙小強, 陸晏天, 白明興, 等. 不同株型玉米基因型對干旱脅迫的響應分析[J]. 草業(yè)學報, 2020, 29(2): 149–162 ZHAO X Q, LU Y T, BAI M X, et al. Response of maize genotypes with different plant architecture to drought stress[J]. Acta Prataculturae Sinica, 2020, 29(2): 149–162

    [10] 趙小強. 玉米株型相關(guān)耐旱遺傳機理研究[D]. 蘭州: 甘肅農(nóng)業(yè)大學, 2018 ZHAO X Q. Genetic mechanisms study of drought tolerance related to plant architecture in maize (L.)[D]. Lanzhou: Gansu Agricultural University, 2018

    [11] 馬娟, 王鐵固, 張懷勝, 等. 玉米葉夾角、葉向值主基因+多基因遺傳模型分析[J]. 河南農(nóng)業(yè)科學, 2012, 41(5): 15–19 MA J, WANG T G, ZHANG H S, et al. The major genes plus polygenes genetic analysis on leaf angle and leaf orientation value in maize[J]. Journal of Henan Agricultural Sciences, 2012, 41(5): 15–19

    [12] GRIFFING B. A generalised treatment of the use of diallel crosses in quantitative inheritance[J]. Heredity, 1956, 10(1): 31–50

    [13] 陳宛秋, 張彪, 康繼偉, 等. 幾個常用玉米自交系株型性狀的遺傳及其在株型育種中的利用[J]. 四川農(nóng)業(yè)大學學報, 1993, 11(4): 563–567 CHEN W Q, ZHANG B, KANG J W, et al. Inheritance of plant characters and yield of some corn inbred lines often used in breeding and its potential ultilization[J]. Journal of Sichuan Agricultural University, 1993, 11(4): 563–567

    [14] MICKELSON S M, STUBER C S, SENIOR L, et al. Quantitative trait loci controlling leaf and tassel traits in a B73×Mo17 population of maize[J]. Crop Science, 2002, 42(6): 1902–1909

    [15] PAN Q C, XU Y C, LI K, et al. The genetic basis of plant architecture in 10 maize recombinant inbred line populations[J]. Plant Physiology, 2017, 175(2): 858–873

    [16] WANG Y J, XU J, DENG D X, et al. A comprehensive meta-analysis of plant morphology, yield, stay-green, and virus disease resistance QTL in maize (L.)[J]. Planta, 2016, 243(2): 459–471

    [17] ZONDERVAN K T, CARDON L R. The complex interplay among factors that influence allelic association[J]. Nature Reviews Genetics, 2004, 5(2): 89–100

    [18] DUCROCQ S, MADUR D, VEYRIERAS J B, et al. Key impact ofon flowering time adaptation in maize: Evidence from association mapping and ecogeographical information[J]. Genetics, 2008, 178(4): 2433–2437

    [19] BRESEGHELLO F, SORRELLS M E. Association mapping of kernel size and milling quality in wheat (L.) cultivars[J]. Genetics, 2006, 172(2): 1165–1177

    [20] EIZENGA G C, AGRAMA H A, LEE F N, et al. Identifying novel resistance genes in newly introduced blast resistant rice germplasm[J]. Crop Science, 2006, 46(5): 1870–1878

    [21] ZHAO X Q, PENG Y L, ZHANG J W, et al. Identification of QTLs and meta-QTLs for seven agronomic traits in multiple maize populations under well-watered and water-stressed conditions[J]. Crop Science, 2018, 58(2): 507–520

    [22] ZHAO X Q, PENG Y L, ZHANG J W, et al. Mapping QTLs and meta-QTLs for two inflorescence architecture traits in multiple maize populations under different watering environments[J]. Molecular Breeding, 2017, 37(7): 91

    [23] 趙小強, 任斌, 彭云玲, 等. 8種水旱環(huán)境下2個玉米群體穗部性狀QTL間的上位性及環(huán)境互作效應分析[J]. 作物學報, 2019, 45(6): 856–871 ZHAO X Q, REN B, PENG Y L, et al. Epistatic and QTL×environment interaction effects for ear related traits in two maize () populations under eight watering environments[J]. Acta Agronomica Sinica, 2019, 45(6): 856–871

    [24] 趙小強, 方鵬, 彭云玲, 等. 基于兩個相關(guān)群體的玉米7個主要農(nóng)藝性狀遺傳分析和QTL定位[J]. 草業(yè)學報, 2018, 27(9): 152–165 ZHAO X Q, FANG P, PENG Y L, et al. Genetic analysis and QTL mapping for seven agronomic traits in maize () using two connected populations[J]. Acta Prataculturae Sinica, 2018, 27(9): 152–165

    [25] PEPPER G E, PEARCE R B, MOCK J J. Leaf orientation and yield of maize[J]. Crop Science, 1977, 17(6): 883–886

    [26] YANG G H, DONG Y B, LI Y L, et al. Integrative detection and verification of QTL for plant traits in two connected RIL populations of high-oil maize[J]. Euphytica, 2015, 206: 203–223

    [27] YANG D L, LIU Y, CHANG H B, et al. Genetic dissection of flag leaf morphology in wheat (L.) under diverse water regimes[J]. BMC Genetics, 2016, 17(1): 94

    [28] 趙小強, 方鵬, 彭云玲, 等. 基于兩個相關(guān)群體的玉米6個穗部性狀QTL定位[J]. 農(nóng)業(yè)生物技術(shù)學報, 2018, 26(5): 729-742 ZHAO X Q, FANG P, PENG Y L, et al. QTL mapping for six ear-related traits based on two maize () related populations[J]. Journal of Agricultural Biotechnology, 2018, 26(5): 729-742

    [29] 張振平, 孫世賢, 張悅, 等. 玉米葉部形態(tài)指標與抗旱性的關(guān)系研究[J]. 玉米科學, 2009, 17(3): 68–70 ZHANG Z P, SUN S X, ZHANG Y, et al. Study on the relations of leaf morphology indexes and drought resistance of maize[J]. Journal of Maize Sciences, 2009, 17(3): 68–70

    [30] ISIDRO J, KNOX R, CLARKE F, et al. Quantitative genetic analysis and mapping of leaf angle in durum wheat[J]. Planta, 2012, 236(6): 1713–1723

    [31] SOHRABI M, RAFII M Y, HANAFI M M, et al. Genetic diversity of upland rice germplasm in Malaysia based on quantitative traits[J]. The Scientific World Journal, 2012, 2012: 416291

    [32] 張飛, 王艷秋, 朱凱, 等. 干旱脅迫下適宜機械化生產(chǎn)高粱品種株型變化及生理響應[J]. 江蘇農(nóng)業(yè)科學, 2018, 46(1): 49–51 ZHANG F, WANG Y Q, ZHU K, et al. Plant type changes and physiological characteristics of sorghum suitable to mechanization under drought stress[J]. Jiangsu Agricultural Sciences, 2018, 46(1): 49–51

    [33] 朱宗河, 鄭文寅, 張學昆. 甘藍型油菜耐旱相關(guān)性狀的主成分分析及綜合評價[J]. 中國農(nóng)業(yè)科學, 2011, 44(9): 1775–1787ZHU Z H, ZHENG W Y, ZHANG X K. Principal component analysis and comprehensive evaluation on morphological and agronomic traits of drought tolerance in rapeseed (L.)[J]. Scientia Agricultura Sinica, 2011, 44(9): 1775–1787

    [34] MASARIRAMBI M T, MANDISODZA F C, MASHINGAIDZE A B, et al. Influence of plant population and seed tuber size on growth and yield components of potato ()[J]. International Journal of Agriculture and Biology, 2012, 14(4): 545–549

    [35] SENIOR M L, MURPHY J P, GOODMAN M M, et al. Utility of SSRs for determining genetic similarities an relationship in maize using an agarose gel system[J]. Crop Science, 1998, 38(4): 1088–1098

    [36] WU X, LI Y X, SHI Y S, et al. Fine genetic characterization of elite maize germplasm using high-throughput SNP genotyping[J]. Theoretical and Applied Genetics, 2014, 127(3): 621–631

    [37] INGVARSSON P K, STREET N R. Association genetics of complex traits in plants[J]. New Phytologist, 2011, 189(4): 909–922

    [38] TUBEROSA R, SALVI S, SANGUINETI M C, et al. Mapping QTLs regulating morpho-physiological traits and yield: Case studies, shortcomings and perspectives in drought-stressed maize[J]. Annals of Botany, 2002, 89(7): 941–963

    [39] LI C H, LI Y X, SHI Y S, et al. Genetic control of the leaf angle and leaf orientation value as revealed by ultra-high density maps in three connected maize populations[J]. PLoS One, 2015, 10(3): e0121624

    [40] BOLDUC N, HAKE S. The maize transcription factor KNOTTED1 directly regulates the gibberellin catabolism gene[J]. The Plant Cell, 2009, 21(6): 1647–1658

    [41] 劉鵬飛, 蔣鋒, 王漢寧, 等. 玉米葉夾角和葉向值的QTL定位[J]. 核農(nóng)學報, 2012, 26(2): 231–237 LIU P F, JIANG F, WANG H N, et al. QTL mapping for leaf angle and leaf orientation in corn[J]. Journal of Nuclear Agricultural Sciences, 2012, 26(2): 231–237

    [42] NAGASAKI H, SAKAMOTO T, SATO Y, et al. Functional analysis of the conserved domains of a rice KNOX homeodomain protein, OSH15[J]. The Plant Cell, 2001, 13(9): 2085–2098

    [43] JUN S E, CHO K H, HWANG J Y, et al. Comparative analysis of the conserved functions of Arabidopsis DRL1 and yeast KTI12[J]. Molecular and Cells, 2015, 38(3): 243–250

    [44] EFRONI I, HAN S K, KIM H J, et al. Regulation of leaf maturation by chromatin-mediated modulation of cytokinin responses[J]. Developmental Cell, 2013, 24(4): 438–445

    Association analysis of SSR markers with leaf morphology in maize () under diverse watering regimes*

    ZHONG Yuan, ZHAO Xiaoqiang**, LI Wenli

    (Gansu Provincial Key Laboratory of Aridland Crop Science, Gansu Agricultural University, Lanzhou 730070, China)

    The leaf morphology of maize () is closely related to its drought resistance. The molecular markers correlated with leaf morphology traits under different watering regimes may provide insights into the genetic mechanisms and identify the breeding varieties. In this study, changes in the leaf length (LL), leaf width (LW), leaf angle (LA), leaf orientation value (LOV), leaf size (LS), leaf shape coefficient (LSC), and leaf rolling index (LRI) of 187 maize inbred lines were analyzed under different watering regimes. A genome-wide scan and the genetic diversity were analyzed with simple-sequence repeat (SSR) markers, and the associations between the markers and the leaf morphological traits were assessed via a general linear model (GLM). The results showed that: 1) the coefficients of variation for the leaf morphological traits of the 187 inbred lines were 13.21%?59.23% and 9.62%?28.22% under normal water and dry conditions, respectively, indicating abundant genetic variation. There were significant differences in the seven leaf morphological traits, so the leaf morphology was adjusted by the hereditary character of the inbred line. LL, LW, LA, and LS decreased significantly under drought stress, whereas LOV, LSC, and LRI increased, and the variance rate of the seven traits was 30.53%?198.31%. 2) A total of 652 alleles were detected using 145 SSR markers. The polymorphism information content (PIC) ranged from 0.201 to 0.966, with an average of 0.478. The tested materials were divided into five groups based on unweighted pair group method with arithmetic mean clustering and population structure; the Luda red cob group, Tangsipingtou group, Lancaster group, P group, and Reid group. 3) A total of 15 SSR markers were associated with seven leaf morphologies under different watering regimes (GLM,<0.01), and the amount of explained phenotypic variance was 2.25%?27.30%. Approximately 72.97% of the SSR markers were detected under drought conditions. The umc1124, umc2363, umc1742, phi331888, umc1378, bnlg1863, umc2134, and umc1345 markers were simultaneously associated with multiple leaf morphologies under different watering regimes, indicating pleiotropy. These results provide useful information for the genetic improvement of leaf morphology and marker-assisted selection breeding programs for maize drought resistance and optimization.

    Maize; Leaf morphology; Drought; Genetic structure; Association analysis

    , E-mail: zhaoxq3324@163.com

    Jun. 27, 2020;

    Aug. 19, 2020

    S513.035.2

    10.13930/j.cnki.cjea.200507

    鐘源, 趙小強, 李文麗. 不同水分環(huán)境下玉米葉形相關(guān)性狀與SSR標記的關(guān)聯(lián)分析[J]. 中國生態(tài)農(nóng)業(yè)學報(中英文), 2020, 28(11): 1726-1738

    ZHONG Y, ZHAO X Q, LI W L. Association analysis of SSR markers with leaf morphology in maize () under diverse watering regimes[J]. Chinese Journal of Eco-Agriculture, 2020, 28(11): 1726-1738

    * 甘肅農(nóng)業(yè)大學甘肅省干旱生境作物學重點實驗室開放基金(GSCS-2019-8, GSCS-2020-5)、甘肅農(nóng)業(yè)大學科技創(chuàng)新基金-公招博士科研啟動基金(GAU-KYQD-2018-19, GAU-KYQD-2018-12)、甘肅省高等學校創(chuàng)新能力提升項目(2019A-052, 2019A-054)、國家自然科學基金項目(32060486)和蘭州青綠儀器技術(shù)有限公司橫向項目(WT20191025)資助

    趙小強, 主要研究方向為現(xiàn)代生物技術(shù)及其在作物遺傳育種中的應用。E-mail: zhaoxq3324@163.com

    鐘源, 主要從事生物化學與分子生物學研究。E-mail: zhongy@gsau.edu.cn

    2020-06-27

    2020-08-19

    * This study was supported by the Open Fund of Gansu Provincial Key Laboratory of Aridland Crop Science, Gansu Agricultural University, China (GSCS-2019-8, GSCS-2020-5), the Scientific Research Start-up Funds for Openly-recruited Doctors of Science and Technology Innovation Funds of Gansu Agricultural University, China (GAU-KYQD-2018-19, GAU-KYQD-2018-12), the Developmental Funds of Innovation Capacity in Higher Education of Gansu Province, China (2019A-052, 2019A-054), National Natural Science Foundation of China (32060486) and the Transverse Project of Lanzhou Qinglü Instrument and Technology Company, China (WT20191025).

    猜你喜歡
    葉形葉長自交系
    馬奶子葡萄葉面積評估模型的建立
    回歸方程法測定兩個龍眼品種葉面積研究
    玉米自交系京92遺傳改良研究
    不同來源玉米自交系穗粒性狀的多樣性分析與改良
    洛陽市13種綠化樹種葉面積的回歸測算
    SRAP結(jié)合SSR標記分析油菜自交系的遺傳多樣性
    干旱脅迫對4份玉米自交系生理與光合特性的影響
    楓葉
    水稻葉形遺傳調(diào)控機理的研究進展
    上位?
    久久人人97超碰香蕉20202| 欧美变态另类bdsm刘玥| 汤姆久久久久久久影院中文字幕| 人人妻人人爽人人添夜夜欢视频| 久久精品国产鲁丝片午夜精品| 水蜜桃什么品种好| 久久 成人 亚洲| 久久久久久伊人网av| 久久女婷五月综合色啪小说| 一区二区日韩欧美中文字幕| 黑丝袜美女国产一区| 中文天堂在线官网| 制服诱惑二区| 成年美女黄网站色视频大全免费| 黄片播放在线免费| 国产成人精品无人区| 美女午夜性视频免费| 在线免费观看不下载黄p国产| 一区二区三区乱码不卡18| 成人二区视频| 岛国毛片在线播放| 少妇被粗大猛烈的视频| 91精品国产国语对白视频| 一级黄片播放器| 自线自在国产av| 免费黄色在线免费观看| 日本午夜av视频| 老女人水多毛片| 免费在线观看视频国产中文字幕亚洲 | 国产片特级美女逼逼视频| 女的被弄到高潮叫床怎么办| 国产精品女同一区二区软件| 日本色播在线视频| 两性夫妻黄色片| 国产精品三级大全| 大片免费播放器 马上看| 2021少妇久久久久久久久久久| 综合色丁香网| 春色校园在线视频观看| 国产女主播在线喷水免费视频网站| 国产成人精品在线电影| 这个男人来自地球电影免费观看 | 热re99久久国产66热| 黄网站色视频无遮挡免费观看| 黑人巨大精品欧美一区二区蜜桃| 又黄又粗又硬又大视频| 各种免费的搞黄视频| 少妇被粗大的猛进出69影院| 一二三四在线观看免费中文在| 国产精品 欧美亚洲| 青青草视频在线视频观看| 国产成人免费无遮挡视频| 亚洲久久久国产精品| 精品国产一区二区久久| 大话2 男鬼变身卡| 国产女主播在线喷水免费视频网站| 菩萨蛮人人尽说江南好唐韦庄| √禁漫天堂资源中文www| 亚洲美女视频黄频| 1024视频免费在线观看| av有码第一页| 亚洲人成网站在线观看播放| videos熟女内射| 老汉色∧v一级毛片| 久久精品夜色国产| 国产在线一区二区三区精| 九草在线视频观看| 亚洲综合色惰| 9191精品国产免费久久| 亚洲精品第二区| 国产精品久久久久久精品古装| 欧美精品亚洲一区二区| 精品少妇内射三级| 熟女少妇亚洲综合色aaa.| 精品一区二区三区四区五区乱码 | 麻豆乱淫一区二区| 国产精品蜜桃在线观看| 99久久人妻综合| 成人影院久久| av福利片在线| 日韩伦理黄色片| 国产福利在线免费观看视频| 欧美最新免费一区二区三区| av福利片在线| 国产成人av激情在线播放| 咕卡用的链子| 99久久精品国产国产毛片| 久久久久久伊人网av| 2022亚洲国产成人精品| 天天操日日干夜夜撸| 亚洲精品aⅴ在线观看| 久久精品国产亚洲av高清一级| 久久久精品国产亚洲av高清涩受| 日韩欧美精品免费久久| 成年动漫av网址| 男的添女的下面高潮视频| 国产精品免费大片| 中文乱码字字幕精品一区二区三区| 亚洲av综合色区一区| 久久久久久人妻| 最近2019中文字幕mv第一页| 国产在视频线精品| 亚洲婷婷狠狠爱综合网| 日韩成人av中文字幕在线观看| 中文字幕av电影在线播放| 国产精品久久久久久久久免| 99精国产麻豆久久婷婷| 可以免费在线观看a视频的电影网站 | 下体分泌物呈黄色| 在线观看国产h片| 久久久久视频综合| 黑人猛操日本美女一级片| 国产精品一区二区在线不卡| 亚洲成人手机| 国产爽快片一区二区三区| 人人妻人人爽人人添夜夜欢视频| 熟女电影av网| 成人毛片60女人毛片免费| 久久久久久久亚洲中文字幕| 麻豆乱淫一区二区| 香蕉精品网在线| 精品人妻一区二区三区麻豆| 男人添女人高潮全过程视频| 久久99一区二区三区| 国产亚洲欧美精品永久| 免费看av在线观看网站| 免费观看在线日韩| 少妇的丰满在线观看| 国产一区二区 视频在线| 校园人妻丝袜中文字幕| 哪个播放器可以免费观看大片| 天天躁夜夜躁狠狠久久av| 久久精品国产亚洲av高清一级| 伊人亚洲综合成人网| 日本色播在线视频| 国产精品三级大全| 久久久久久人人人人人| 色94色欧美一区二区| 免费高清在线观看视频在线观看| 国产成人91sexporn| 国产97色在线日韩免费| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产欧美在线一区| 国产男女内射视频| 黄色怎么调成土黄色| 亚洲国产最新在线播放| 亚洲四区av| 男人舔女人的私密视频| 午夜福利网站1000一区二区三区| 国产男女内射视频| 欧美 亚洲 国产 日韩一| 啦啦啦在线免费观看视频4| av在线播放精品| 亚洲精品久久午夜乱码| av天堂久久9| 女人精品久久久久毛片| 狠狠婷婷综合久久久久久88av| 一个人免费看片子| 国产黄色视频一区二区在线观看| 高清欧美精品videossex| 狂野欧美激情性bbbbbb| 曰老女人黄片| 丝袜美足系列| 日韩中文字幕视频在线看片| 成人亚洲精品一区在线观看| 久久精品国产鲁丝片午夜精品| 麻豆乱淫一区二区| 在线观看免费高清a一片| 国精品久久久久久国模美| 国产精品免费视频内射| 日韩电影二区| 永久免费av网站大全| 亚洲av福利一区| 精品少妇久久久久久888优播| 日本欧美国产在线视频| 亚洲综合精品二区| 美国免费a级毛片| 九色亚洲精品在线播放| 欧美 日韩 精品 国产| 成年女人在线观看亚洲视频| 国产精品av久久久久免费| 欧美日韩亚洲国产一区二区在线观看 | 久久精品久久久久久久性| 国产精品99久久99久久久不卡 | 亚洲欧美日韩另类电影网站| 久久精品国产亚洲av天美| 成年美女黄网站色视频大全免费| 有码 亚洲区| 水蜜桃什么品种好| 一区福利在线观看| 国产精品二区激情视频| 亚洲,欧美,日韩| 99re6热这里在线精品视频| 国产av精品麻豆| 91国产中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 精品视频人人做人人爽| 国产女主播在线喷水免费视频网站| 性少妇av在线| 人成视频在线观看免费观看| 天堂中文最新版在线下载| 亚洲国产精品一区二区三区在线| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产欧美在线一区| 精品国产国语对白av| 国产日韩欧美在线精品| 国产一区二区三区av在线| 欧美人与性动交α欧美软件| 看十八女毛片水多多多| 精品第一国产精品| 国产乱人偷精品视频| av.在线天堂| 26uuu在线亚洲综合色| 在线观看三级黄色| 色哟哟·www| 有码 亚洲区| 亚洲,一卡二卡三卡| 午夜91福利影院| 黄色 视频免费看| 国产精品不卡视频一区二区| 18在线观看网站| 2018国产大陆天天弄谢| 伦理电影免费视频| 国产熟女午夜一区二区三区| 69精品国产乱码久久久| 欧美精品av麻豆av| 午夜老司机福利剧场| 老熟女久久久| 大香蕉久久成人网| 免费在线观看视频国产中文字幕亚洲 | 成人亚洲欧美一区二区av| 美女脱内裤让男人舔精品视频| 波多野结衣一区麻豆| 丰满乱子伦码专区| 久久国产亚洲av麻豆专区| 成人二区视频| 国产成人精品久久久久久| 日韩成人av中文字幕在线观看| 国产色婷婷99| 少妇人妻久久综合中文| 亚洲国产色片| 天天影视国产精品| 人人妻人人澡人人爽人人夜夜| 欧美日韩一区二区视频在线观看视频在线| 最新的欧美精品一区二区| 国产极品粉嫩免费观看在线| 人人澡人人妻人| 777久久人妻少妇嫩草av网站| 男女无遮挡免费网站观看| 国产精品国产三级专区第一集| av电影中文网址| 国产探花极品一区二区| www.精华液| 日日啪夜夜爽| 成人国产麻豆网| 一边摸一边做爽爽视频免费| 母亲3免费完整高清在线观看 | 国产乱人偷精品视频| 亚洲第一青青草原| 青草久久国产| 天美传媒精品一区二区| 免费观看av网站的网址| 大话2 男鬼变身卡| 日韩大片免费观看网站| 毛片一级片免费看久久久久| 成人二区视频| 中文字幕人妻丝袜一区二区 | 久久久a久久爽久久v久久| 欧美日韩综合久久久久久| 少妇人妻 视频| 咕卡用的链子| 美女午夜性视频免费| av在线观看视频网站免费| 99久国产av精品国产电影| 美女高潮到喷水免费观看| 丰满乱子伦码专区| 午夜免费观看性视频| 99久久中文字幕三级久久日本| 欧美成人精品欧美一级黄| 视频在线观看一区二区三区| 亚洲欧美一区二区三区国产| 中文天堂在线官网| 精品少妇久久久久久888优播| 韩国精品一区二区三区| 免费观看av网站的网址| 亚洲av男天堂| 国产综合精华液| 国产乱来视频区| 老汉色∧v一级毛片| 欧美+日韩+精品| 午夜福利网站1000一区二区三区| xxx大片免费视频| 高清欧美精品videossex| 男人添女人高潮全过程视频| 美女国产视频在线观看| 熟妇人妻不卡中文字幕| 日韩精品有码人妻一区| 18禁观看日本| 国产熟女午夜一区二区三区| 丝袜在线中文字幕| 国产极品天堂在线| 国产日韩欧美视频二区| 1024视频免费在线观看| 在线免费观看不下载黄p国产| 寂寞人妻少妇视频99o| 欧美精品av麻豆av| 三上悠亚av全集在线观看| 中文字幕另类日韩欧美亚洲嫩草| 性少妇av在线| 毛片一级片免费看久久久久| 国产一级毛片在线| 男女无遮挡免费网站观看| 亚洲av欧美aⅴ国产| 色吧在线观看| 日韩免费高清中文字幕av| 亚洲国产毛片av蜜桃av| 亚洲精品日韩在线中文字幕| 亚洲国产精品成人久久小说| 国产成人精品久久二区二区91 | 性色avwww在线观看| 精品酒店卫生间| 两个人看的免费小视频| 18+在线观看网站| 熟女少妇亚洲综合色aaa.| 欧美av亚洲av综合av国产av | 人人妻人人添人人爽欧美一区卜| 亚洲精品一区蜜桃| 中文天堂在线官网| 欧美另类一区| 欧美人与性动交α欧美软件| 韩国av在线不卡| 伊人久久国产一区二区| 免费女性裸体啪啪无遮挡网站| 国产亚洲欧美精品永久| 日韩 亚洲 欧美在线| 久久久久精品人妻al黑| 日韩av在线免费看完整版不卡| 欧美另类一区| 国产深夜福利视频在线观看| 国产极品粉嫩免费观看在线| 日韩制服丝袜自拍偷拍| 如日韩欧美国产精品一区二区三区| av电影中文网址| 亚洲视频免费观看视频| 欧美人与性动交α欧美精品济南到 | 日韩,欧美,国产一区二区三区| 男女边摸边吃奶| 色婷婷av一区二区三区视频| av在线播放精品| 纯流量卡能插随身wifi吗| 香蕉丝袜av| 制服人妻中文乱码| 热99久久久久精品小说推荐| 精品第一国产精品| 亚洲在久久综合| 在线观看美女被高潮喷水网站| 欧美日韩视频精品一区| 欧美日韩国产mv在线观看视频| 免费黄频网站在线观看国产| 国产成人精品无人区| 日韩成人av中文字幕在线观看| 日韩视频在线欧美| 亚洲国产欧美日韩在线播放| 人人妻人人澡人人爽人人夜夜| 国产亚洲最大av| 涩涩av久久男人的天堂| 国产黄色免费在线视频| 亚洲欧美一区二区三区久久| 久久久欧美国产精品| 国产一区有黄有色的免费视频| 老鸭窝网址在线观看| av在线播放精品| 在线观看免费视频网站a站| 日本vs欧美在线观看视频| 婷婷色av中文字幕| 国产黄色免费在线视频| 国产精品久久久久久精品电影小说| 成人免费观看视频高清| 丰满迷人的少妇在线观看| 在线观看www视频免费| 亚洲三区欧美一区| 亚洲国产日韩一区二区| 成人国产麻豆网| 国语对白做爰xxxⅹ性视频网站| 秋霞伦理黄片| 男女边摸边吃奶| xxxhd国产人妻xxx| 国产野战对白在线观看| 亚洲激情五月婷婷啪啪| 亚洲,一卡二卡三卡| 国产色婷婷99| 欧美人与性动交α欧美精品济南到 | 边亲边吃奶的免费视频| 国产男女超爽视频在线观看| 午夜免费鲁丝| 精品一品国产午夜福利视频| 久久久精品国产亚洲av高清涩受| 亚洲三区欧美一区| 男女下面插进去视频免费观看| 国产免费一区二区三区四区乱码| 久久精品国产亚洲av涩爱| www.熟女人妻精品国产| 91成人精品电影| 国产熟女午夜一区二区三区| 人妻少妇偷人精品九色| 成年人免费黄色播放视频| 免费看不卡的av| 免费高清在线观看视频在线观看| 午夜91福利影院| 欧美变态另类bdsm刘玥| 两性夫妻黄色片| 好男人视频免费观看在线| 亚洲人成网站在线观看播放| 最黄视频免费看| 久久国内精品自在自线图片| 成人午夜精彩视频在线观看| 久久精品亚洲av国产电影网| 在线观看免费高清a一片| 日韩欧美一区视频在线观看| 亚洲美女视频黄频| 天堂俺去俺来也www色官网| 电影成人av| 一本久久精品| 亚洲,欧美,日韩| av免费在线看不卡| 最近2019中文字幕mv第一页| 亚洲,欧美,日韩| 久久久久久伊人网av| 亚洲精品第二区| 大陆偷拍与自拍| 欧美另类一区| 91成人精品电影| 最新的欧美精品一区二区| 天天躁日日躁夜夜躁夜夜| 午夜免费观看性视频| 久久这里只有精品19| 欧美成人午夜精品| 制服丝袜香蕉在线| 综合色丁香网| 国产精品久久久av美女十八| 久久久亚洲精品成人影院| 国产一区二区 视频在线| 日韩制服丝袜自拍偷拍| 永久网站在线| 亚洲 欧美一区二区三区| 亚洲国产精品国产精品| 精品福利永久在线观看| 大片免费播放器 马上看| 久久精品人人爽人人爽视色| 女性被躁到高潮视频| 亚洲精品日韩在线中文字幕| 国产欧美日韩一区二区三区在线| 99热网站在线观看| 久久婷婷青草| 看免费av毛片| 精品国产露脸久久av麻豆| 美女国产高潮福利片在线看| 国产激情久久老熟女| 久久av网站| 久久韩国三级中文字幕| 欧美激情高清一区二区三区 | 婷婷色av中文字幕| 色婷婷av一区二区三区视频| 秋霞伦理黄片| 免费久久久久久久精品成人欧美视频| 国产极品粉嫩免费观看在线| 人妻系列 视频| 最近最新中文字幕大全免费视频 | 少妇的逼水好多| 亚洲美女视频黄频| 欧美日韩精品网址| 夫妻性生交免费视频一级片| 一级毛片电影观看| 午夜福利视频在线观看免费| 成人国语在线视频| 欧美变态另类bdsm刘玥| 亚洲国产精品一区三区| 日韩欧美一区视频在线观看| 大码成人一级视频| 久久久精品区二区三区| 色网站视频免费| 成人18禁高潮啪啪吃奶动态图| 999久久久国产精品视频| 最近最新中文字幕大全免费视频 | 伦精品一区二区三区| 亚洲综合色惰| 成年av动漫网址| 99香蕉大伊视频| 国产精品久久久久久精品古装| 亚洲天堂av无毛| 亚洲国产色片| 色网站视频免费| 中国国产av一级| 亚洲人成网站在线观看播放| 中文精品一卡2卡3卡4更新| 亚洲少妇的诱惑av| 国产国语露脸激情在线看| 国产精品免费视频内射| 老司机影院毛片| 国产精品 欧美亚洲| 国产日韩欧美亚洲二区| 如日韩欧美国产精品一区二区三区| av女优亚洲男人天堂| 国产白丝娇喘喷水9色精品| 校园人妻丝袜中文字幕| 亚洲精华国产精华液的使用体验| 亚洲国产av影院在线观看| 涩涩av久久男人的天堂| av一本久久久久| 亚洲成av片中文字幕在线观看 | 大片免费播放器 马上看| 在线观看免费视频网站a站| 成人国产av品久久久| 久久久久国产网址| 精品亚洲乱码少妇综合久久| 亚洲精品自拍成人| 欧美精品亚洲一区二区| 欧美精品一区二区大全| 日韩免费高清中文字幕av| 夫妻性生交免费视频一级片| av一本久久久久| 国产xxxxx性猛交| 日本午夜av视频| 人妻少妇偷人精品九色| 天堂中文最新版在线下载| 王馨瑶露胸无遮挡在线观看| 成人免费观看视频高清| 大片免费播放器 马上看| 久久女婷五月综合色啪小说| 午夜福利视频精品| 九草在线视频观看| 伦理电影免费视频| 热99久久久久精品小说推荐| 国产精品人妻久久久影院| 黄色 视频免费看| 久久99蜜桃精品久久| 午夜激情av网站| 国产成人精品久久久久久| av国产久精品久网站免费入址| 女人久久www免费人成看片| 啦啦啦在线免费观看视频4| 不卡视频在线观看欧美| 天美传媒精品一区二区| 国产探花极品一区二区| 亚洲av成人精品一二三区| 日韩一区二区三区影片| 丝袜美腿诱惑在线| 毛片一级片免费看久久久久| 亚洲成色77777| 熟女少妇亚洲综合色aaa.| 久久韩国三级中文字幕| 久久久久国产网址| 国产成人免费无遮挡视频| 国产精品.久久久| 王馨瑶露胸无遮挡在线观看| 亚洲欧美一区二区三区黑人 | 伊人久久大香线蕉亚洲五| 国产精品久久久av美女十八| 亚洲人成电影观看| 亚洲精品在线美女| 国产精品一国产av| 亚洲精品中文字幕在线视频| 免费高清在线观看视频在线观看| 久久精品久久久久久久性| 成人黄色视频免费在线看| 久热久热在线精品观看| 人人妻人人澡人人看| 乱人伦中国视频| 亚洲国产av影院在线观看| 国产一区有黄有色的免费视频| 亚洲情色 制服丝袜| 久久国产精品大桥未久av| 久久久精品区二区三区| 久久久久视频综合| 亚洲成av片中文字幕在线观看 | 啦啦啦啦在线视频资源| 日韩伦理黄色片| 精品亚洲乱码少妇综合久久| 亚洲国产av影院在线观看| 亚洲伊人色综图| 日韩大片免费观看网站| 少妇精品久久久久久久| 日韩中文字幕视频在线看片| 欧美精品高潮呻吟av久久| 亚洲av福利一区| 国产又爽黄色视频| 欧美国产精品一级二级三级| 1024香蕉在线观看| 亚洲色图综合在线观看| 精品国产一区二区三区久久久樱花| 午夜福利影视在线免费观看| 午夜福利网站1000一区二区三区| 精品少妇一区二区三区视频日本电影 | 亚洲av免费高清在线观看| 91午夜精品亚洲一区二区三区| 精品国产超薄肉色丝袜足j| 亚洲色图 男人天堂 中文字幕| 亚洲伊人色综图| 高清av免费在线| 18禁国产床啪视频网站| 丰满迷人的少妇在线观看| 亚洲精品美女久久av网站| 在线观看免费高清a一片| 成人免费观看视频高清| 精品卡一卡二卡四卡免费| 国产精品.久久久| 亚洲精品av麻豆狂野| av有码第一页| 美国免费a级毛片| av网站免费在线观看视频| 一区二区三区激情视频| 波多野结衣一区麻豆| 亚洲少妇的诱惑av| 国产精品免费大片| 久久精品国产亚洲av涩爱|