• <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)控機理的研究進展
    上位?
    av线在线观看网站| 亚洲av成人精品一区久久| 日本wwww免费看| 在线精品无人区一区二区三| 男人爽女人下面视频在线观看| 最新中文字幕久久久久| 高清视频免费观看一区二区| 狂野欧美白嫩少妇大欣赏| 亚洲综合色网址| 亚洲欧美一区二区三区黑人 | 亚洲成人av在线免费| 在线亚洲精品国产二区图片欧美 | 美女内射精品一级片tv| 亚洲精品美女久久av网站| 黄色视频在线播放观看不卡| 五月玫瑰六月丁香| 你懂的网址亚洲精品在线观看| 久久精品国产a三级三级三级| 一区二区三区四区激情视频| 国产 一区精品| 欧美激情国产日韩精品一区| 韩国av在线不卡| 久久久久人妻精品一区果冻| 午夜日本视频在线| 久久这里有精品视频免费| 美女国产高潮福利片在线看| 777米奇影视久久| 男女免费视频国产| 天天影视国产精品| 国内精品宾馆在线| 欧美精品国产亚洲| 91久久精品电影网| 考比视频在线观看| 亚洲五月色婷婷综合| 狂野欧美激情性xxxx在线观看| 午夜日本视频在线| 日韩成人伦理影院| 热re99久久精品国产66热6| 久久影院123| 国产精品一国产av| 日本黄色日本黄色录像| 三上悠亚av全集在线观看| 一二三四中文在线观看免费高清| 精品国产乱码久久久久久小说| 国产白丝娇喘喷水9色精品| 内地一区二区视频在线| 少妇猛男粗大的猛烈进出视频| 搡女人真爽免费视频火全软件| 人人妻人人澡人人爽人人夜夜| 十八禁高潮呻吟视频| 国产亚洲精品久久久com| 成人国产av品久久久| 超色免费av| 在线精品无人区一区二区三| 国产欧美日韩一区二区三区在线 | 蜜桃在线观看..| 99热这里只有是精品在线观看| 日韩成人av中文字幕在线观看| 亚洲天堂av无毛| 成人亚洲精品一区在线观看| videos熟女内射| 五月天丁香电影| 免费av中文字幕在线| 91国产中文字幕| 国产在线免费精品| 久久精品国产鲁丝片午夜精品| 一级毛片我不卡| 精品久久久久久久久av| 伊人久久国产一区二区| 精品久久久久久久久av| 汤姆久久久久久久影院中文字幕| 伦精品一区二区三区| 内地一区二区视频在线| 日韩免费高清中文字幕av| 亚洲av福利一区| 欧美最新免费一区二区三区| 国国产精品蜜臀av免费| 卡戴珊不雅视频在线播放| 91精品三级在线观看| 久久国产亚洲av麻豆专区| 国产有黄有色有爽视频| 国内精品宾馆在线| 国产永久视频网站| 国产乱来视频区| 最黄视频免费看| 免费高清在线观看日韩| 成年女人在线观看亚洲视频| 欧美3d第一页| 欧美日韩国产mv在线观看视频| 国产精品三级大全| 免费观看a级毛片全部| 插逼视频在线观看| 赤兔流量卡办理| 中文字幕免费在线视频6| 热99国产精品久久久久久7| 国产精品女同一区二区软件| 国产伦精品一区二区三区视频9| 亚洲精品aⅴ在线观看| 伊人久久精品亚洲午夜| 欧美精品高潮呻吟av久久| 三级国产精品片| 久久午夜福利片| 黑人猛操日本美女一级片| 久久久久精品久久久久真实原创| 久久99一区二区三区| 精品人妻熟女毛片av久久网站| 午夜福利视频在线观看免费| 少妇熟女欧美另类| 国产精品久久久久久av不卡| 色哟哟·www| 日本wwww免费看| 这个男人来自地球电影免费观看 | 国产在线一区二区三区精| 丝袜脚勾引网站| 久久 成人 亚洲| 一级毛片电影观看| 草草在线视频免费看| 亚洲精品亚洲一区二区| 女人久久www免费人成看片| 一区二区日韩欧美中文字幕 | 黄片播放在线免费| 国产黄色免费在线视频| 欧美老熟妇乱子伦牲交| 亚洲综合精品二区| 亚洲精品久久午夜乱码| 少妇的逼好多水| 亚洲怡红院男人天堂| 国产极品天堂在线| 精品国产露脸久久av麻豆| 麻豆成人av视频| 青春草视频在线免费观看| 亚洲国产成人一精品久久久| 亚洲欧美中文字幕日韩二区| 91精品伊人久久大香线蕉| 亚洲精品国产色婷婷电影| 最近最新中文字幕免费大全7| 中文欧美无线码| 高清黄色对白视频在线免费看| 久久久久国产网址| 国产成人午夜福利电影在线观看| 亚洲少妇的诱惑av| 美女视频免费永久观看网站| 黄色视频在线播放观看不卡| 日韩欧美精品免费久久| 国产综合精华液| 人人妻人人澡人人爽人人夜夜| 考比视频在线观看| 涩涩av久久男人的天堂| 美女视频免费永久观看网站| 黄色视频在线播放观看不卡| 国语对白做爰xxxⅹ性视频网站| 黄色怎么调成土黄色| 欧美国产精品一级二级三级| 一级毛片 在线播放| 亚洲熟女精品中文字幕| 日韩av在线免费看完整版不卡| 九九在线视频观看精品| 99re6热这里在线精品视频| 午夜免费观看性视频| av在线老鸭窝| 高清毛片免费看| 亚洲精品456在线播放app| 九色亚洲精品在线播放| 99久久精品国产国产毛片| 91精品国产国语对白视频| 青青草视频在线视频观看| 精品久久久久久久久av| 桃花免费在线播放| 看非洲黑人一级黄片| 免费高清在线观看视频在线观看| 亚洲av日韩在线播放| 最近中文字幕2019免费版| 久久精品国产亚洲av涩爱| 免费看不卡的av| 亚洲精品久久成人aⅴ小说 | 成年人免费黄色播放视频| 亚洲精华国产精华液的使用体验| 午夜影院在线不卡| 新久久久久国产一级毛片| 蜜桃久久精品国产亚洲av| 91精品国产九色| 欧美丝袜亚洲另类| 亚洲五月色婷婷综合| 国产精品一区www在线观看| 婷婷色综合www| 丰满少妇做爰视频| 毛片一级片免费看久久久久| 视频中文字幕在线观看| 国产精品一区二区在线不卡| 欧美国产精品一级二级三级| 日韩av免费高清视频| av又黄又爽大尺度在线免费看| 桃花免费在线播放| 亚洲伊人久久精品综合| 一区二区av电影网| 欧美bdsm另类| 欧美日韩视频精品一区| 黄片无遮挡物在线观看| 91国产中文字幕| 大香蕉久久成人网| 亚洲精品久久久久久婷婷小说| av网站免费在线观看视频| 两个人的视频大全免费| 亚洲欧美精品自产自拍| 国产av精品麻豆| 亚洲成人av在线免费| 成人亚洲欧美一区二区av| 国产免费现黄频在线看| 岛国毛片在线播放| 国产欧美亚洲国产| 久久久久国产网址| 国产在线免费精品| 性高湖久久久久久久久免费观看| 欧美激情极品国产一区二区三区 | 亚洲丝袜综合中文字幕| 中文字幕久久专区| 亚洲三级黄色毛片| 亚洲三级黄色毛片| h视频一区二区三区| 考比视频在线观看| 亚洲欧美成人精品一区二区| 国产精品国产av在线观看| 久久久欧美国产精品| 久久亚洲国产成人精品v| 亚洲无线观看免费| 久久av网站| 国产成人一区二区在线| 精品人妻熟女毛片av久久网站| 久久久久国产网址| 成人免费观看视频高清| 久久国产精品大桥未久av| 在线天堂最新版资源| 国产免费福利视频在线观看| 天堂中文最新版在线下载| 制服丝袜香蕉在线| 久久99一区二区三区| 国产永久视频网站| 欧美亚洲 丝袜 人妻 在线| 黄片无遮挡物在线观看| 日本午夜av视频| 草草在线视频免费看| 国产男女内射视频| 精品一区在线观看国产| 亚洲国产精品一区三区| 亚洲精品久久午夜乱码| 嫩草影院入口| 午夜福利影视在线免费观看| 老司机影院成人| 少妇 在线观看| 婷婷成人精品国产| 特大巨黑吊av在线直播| 简卡轻食公司| 在线观看美女被高潮喷水网站| 日韩伦理黄色片| 你懂的网址亚洲精品在线观看| 久久狼人影院| 亚洲av综合色区一区| 亚洲成人手机| 高清av免费在线| 黄片无遮挡物在线观看| 亚洲成色77777| 日韩电影二区| 国产亚洲一区二区精品| 毛片一级片免费看久久久久| 日韩av免费高清视频| 看非洲黑人一级黄片| 亚洲国产精品一区二区三区在线| 欧美激情 高清一区二区三区| av线在线观看网站| 性色av一级| 亚洲av在线观看美女高潮| 欧美97在线视频| 伦精品一区二区三区| 精品国产一区二区三区久久久樱花| 边亲边吃奶的免费视频| 亚洲熟女精品中文字幕| 亚洲一级一片aⅴ在线观看| 欧美另类一区| 在线观看免费高清a一片| 国产片特级美女逼逼视频| 国产日韩欧美视频二区| 内地一区二区视频在线| 日韩精品免费视频一区二区三区 | 99热全是精品| 乱人伦中国视频| 午夜免费男女啪啪视频观看| 看十八女毛片水多多多| 永久免费av网站大全| 五月伊人婷婷丁香| 大片电影免费在线观看免费| 日本av手机在线免费观看| 春色校园在线视频观看| 国产视频内射| 国产亚洲av片在线观看秒播厂| 国产又色又爽无遮挡免| .国产精品久久| 2022亚洲国产成人精品| 日韩制服骚丝袜av| 中文字幕最新亚洲高清| 国语对白做爰xxxⅹ性视频网站| 国产日韩欧美视频二区| 日日啪夜夜爽| 国产成人一区二区在线| 亚洲第一区二区三区不卡| 91国产中文字幕| 777米奇影视久久| 国产高清有码在线观看视频| 国产老妇伦熟女老妇高清| 亚洲国产成人一精品久久久| 日韩一区二区视频免费看| 国产白丝娇喘喷水9色精品| 交换朋友夫妻互换小说| 一区二区av电影网| 成人18禁高潮啪啪吃奶动态图 | 亚洲国产成人一精品久久久| 久久婷婷青草| 亚洲,欧美,日韩| 在线播放无遮挡| 岛国毛片在线播放| 国产在线免费精品| 视频在线观看一区二区三区| 十分钟在线观看高清视频www| 综合色丁香网| 欧美日韩av久久| 久久精品国产a三级三级三级| 91精品一卡2卡3卡4卡| 一级毛片黄色毛片免费观看视频| 制服诱惑二区| 国国产精品蜜臀av免费| 日韩免费高清中文字幕av| 欧美 日韩 精品 国产| 成人午夜精彩视频在线观看| 美女脱内裤让男人舔精品视频| 日本av免费视频播放| 亚洲欧洲日产国产| 亚洲国产最新在线播放| 丝袜在线中文字幕| 80岁老熟妇乱子伦牲交| 国产熟女欧美一区二区| 五月伊人婷婷丁香| 国产精品一区二区在线不卡| 国内精品宾馆在线| 亚洲不卡免费看| 人妻一区二区av| 久久热精品热| 黑人高潮一二区| 久久午夜福利片| 国产精品久久久久久久久免| 日本欧美视频一区| 亚洲av二区三区四区| 欧美日韩在线观看h| 激情五月婷婷亚洲| 在现免费观看毛片| 中文字幕人妻熟人妻熟丝袜美| 日韩,欧美,国产一区二区三区| 午夜日本视频在线| 高清av免费在线| 亚洲欧美清纯卡通| 久久精品国产亚洲av天美| 日韩在线高清观看一区二区三区| 久久人妻熟女aⅴ| 99热网站在线观看| 91精品国产国语对白视频| 国产69精品久久久久777片| 考比视频在线观看| 国产女主播在线喷水免费视频网站| 亚洲av.av天堂| 国产精品久久久久久av不卡| 丝袜在线中文字幕| 日韩成人伦理影院| 最新的欧美精品一区二区| 午夜激情久久久久久久| 最近手机中文字幕大全| 日日啪夜夜爽| 亚洲成色77777| 国产极品粉嫩免费观看在线 | 久久人人爽人人片av| 午夜免费观看性视频| 制服丝袜香蕉在线| 一本色道久久久久久精品综合| 欧美日韩国产mv在线观看视频| 欧美成人午夜免费资源| 精品酒店卫生间| 少妇被粗大猛烈的视频| 校园人妻丝袜中文字幕| 欧美丝袜亚洲另类| 亚洲国产精品999| 只有这里有精品99| 亚洲国产精品一区三区| 国产日韩欧美视频二区| 免费看不卡的av| 国产黄色免费在线视频| 亚洲欧洲精品一区二区精品久久久 | 国产乱人偷精品视频| 欧美最新免费一区二区三区| 日本wwww免费看| 99九九在线精品视频| 极品少妇高潮喷水抽搐| 精品国产露脸久久av麻豆| 久久久久久久久大av| av国产久精品久网站免费入址| av线在线观看网站| 免费高清在线观看日韩| 亚洲综合色惰| 国产在视频线精品| 新久久久久国产一级毛片| 成人免费观看视频高清| 日本黄大片高清| 亚洲av欧美aⅴ国产| 如日韩欧美国产精品一区二区三区 | 国产深夜福利视频在线观看| 精品久久久噜噜| 国产成人一区二区在线| 国产精品一区二区在线观看99| 日本-黄色视频高清免费观看| 亚洲精品视频女| 国产国语露脸激情在线看| 爱豆传媒免费全集在线观看| 亚洲国产av影院在线观看| 日本wwww免费看| 一本色道久久久久久精品综合| 搡老乐熟女国产| 能在线免费看毛片的网站| 久久精品国产亚洲网站| 新久久久久国产一级毛片| 亚洲精品一二三| 久久99精品国语久久久| 国产黄片视频在线免费观看| 九九爱精品视频在线观看| 夜夜看夜夜爽夜夜摸| 亚洲欧洲日产国产| 国产69精品久久久久777片| 岛国毛片在线播放| 99热网站在线观看| 精品一区二区三区视频在线| 国产精品国产av在线观看| 亚洲天堂av无毛| 亚洲国产最新在线播放| 大又大粗又爽又黄少妇毛片口| 亚洲av中文av极速乱| 精品少妇久久久久久888优播| 亚洲图色成人| 亚洲av日韩在线播放| 久久久久人妻精品一区果冻| 18禁动态无遮挡网站| 国产av码专区亚洲av| 有码 亚洲区| 婷婷成人精品国产| 最近中文字幕2019免费版| 亚洲内射少妇av| 婷婷色av中文字幕| 亚洲成人av在线免费| 2018国产大陆天天弄谢| 天堂俺去俺来也www色官网| 伦理电影大哥的女人| 日韩中文字幕视频在线看片| 亚洲欧美日韩另类电影网站| 久久鲁丝午夜福利片| 亚州av有码| a级毛色黄片| 久久女婷五月综合色啪小说| 少妇的逼好多水| 91精品国产九色| 亚洲第一区二区三区不卡| 在线观看国产h片| 熟妇人妻不卡中文字幕| 日本黄大片高清| 青青草视频在线视频观看| 亚洲少妇的诱惑av| 亚洲激情五月婷婷啪啪| 高清不卡的av网站| 黄色欧美视频在线观看| 欧美亚洲 丝袜 人妻 在线| 嫩草影院入口| av视频免费观看在线观看| 一区二区三区精品91| 国产成人freesex在线| 精品国产乱码久久久久久小说| 久久久精品区二区三区| a级毛片在线看网站| 国产老妇伦熟女老妇高清| 免费观看av网站的网址| 最后的刺客免费高清国语| 欧美最新免费一区二区三区| 亚洲成人手机| 人人妻人人澡人人爽人人夜夜| 99久久精品国产国产毛片| 久久热精品热| 久久久久久久大尺度免费视频| 一本色道久久久久久精品综合| 国产精品人妻久久久久久| 91精品国产国语对白视频| 亚洲av日韩在线播放| 在线观看www视频免费| 美女国产视频在线观看| 精品人妻在线不人妻| av在线播放精品| 久久人妻熟女aⅴ| a级毛色黄片| 国产一区二区在线观看av| h视频一区二区三区| 91在线精品国自产拍蜜月| 日本-黄色视频高清免费观看| 三级国产精品片| 精品一区二区免费观看| 国产精品国产av在线观看| 午夜福利视频精品| 久久热精品热| 蜜桃久久精品国产亚洲av| 九色成人免费人妻av| 日本黄色日本黄色录像| 国产精品嫩草影院av在线观看| 91国产中文字幕| 丁香六月天网| 欧美性感艳星| 狠狠婷婷综合久久久久久88av| 99热这里只有是精品在线观看| 亚洲国产精品专区欧美| kizo精华| 国产免费一区二区三区四区乱码| 母亲3免费完整高清在线观看 | 日日爽夜夜爽网站| 国产女主播在线喷水免费视频网站| 国产一级毛片在线| 中文精品一卡2卡3卡4更新| 久久久久久久久久成人| 免费观看的影片在线观看| 91国产中文字幕| 一个人免费看片子| 国产成人av激情在线播放 | 国产成人a∨麻豆精品| 伊人亚洲综合成人网| 狠狠精品人妻久久久久久综合| 亚洲内射少妇av| av有码第一页| 最近2019中文字幕mv第一页| 国产一区二区在线观看av| 国产又色又爽无遮挡免| 精品久久久精品久久久| 国产精品成人在线| 亚洲av成人精品一区久久| 天堂8中文在线网| 肉色欧美久久久久久久蜜桃| 亚洲成人一二三区av| 大码成人一级视频| 欧美日本中文国产一区发布| 曰老女人黄片| 26uuu在线亚洲综合色| 久热这里只有精品99| 嫩草影院入口| 欧美亚洲 丝袜 人妻 在线| 国产亚洲av片在线观看秒播厂| 亚洲情色 制服丝袜| 永久网站在线| 在线亚洲精品国产二区图片欧美 | 九草在线视频观看| 欧美精品亚洲一区二区| 国产伦精品一区二区三区视频9| 国产成人午夜福利电影在线观看| 人妻人人澡人人爽人人| 免费人妻精品一区二区三区视频| 又黄又爽又刺激的免费视频.| 精品少妇久久久久久888优播| 高清黄色对白视频在线免费看| 午夜精品国产一区二区电影| 精品一区在线观看国产| 国产熟女欧美一区二区| 国产一级毛片在线| 精品久久久久久久久av| 亚洲人成77777在线视频| 麻豆乱淫一区二区| 久久久久视频综合| 国产欧美亚洲国产| 97精品久久久久久久久久精品| 女人精品久久久久毛片| 亚洲在久久综合| 国产精品久久久久久精品电影小说| 搡女人真爽免费视频火全软件| 97超碰精品成人国产| 国产亚洲精品久久久com| av福利片在线| 午夜免费男女啪啪视频观看| 最近手机中文字幕大全| 美女主播在线视频| 国产精品国产av在线观看| 亚洲精品国产av成人精品| 观看av在线不卡| 午夜视频国产福利| 亚洲精品国产av成人精品| 美女主播在线视频| 美女中出高潮动态图| 国产精品女同一区二区软件| 夫妻性生交免费视频一级片| 久久久久国产精品人妻一区二区| 最近最新中文字幕免费大全7| 成年人免费黄色播放视频| 狂野欧美白嫩少妇大欣赏| 亚洲一级一片aⅴ在线观看| 岛国毛片在线播放| 三级国产精品片| 大片电影免费在线观看免费| 老司机影院毛片| 婷婷成人精品国产| 久久久精品94久久精品| 人体艺术视频欧美日本| 特大巨黑吊av在线直播| 亚洲欧洲精品一区二区精品久久久 | 在线观看三级黄色| 亚洲经典国产精华液单| 性色av一级| 色哟哟·www| 日本免费在线观看一区| 成年av动漫网址| 欧美日韩视频高清一区二区三区二| 母亲3免费完整高清在线观看 |