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

    黍稷高基元EST-SSR標(biāo)記開發(fā)及200份核心種質(zhì)資源遺傳多樣性分析

    2023-06-09 09:23:44張立媛劉少雄薛亞鵬王瑞云劉敏軒
    作物學(xué)報(bào) 2023年8期
    關(guān)鍵詞:糜子基元種質(zhì)

    王 倩 張立媛 許 月 李 海 劉少雄 薛亞鵬 陸 平 王瑞云 劉敏軒,*

    研究簡(jiǎn)報(bào)

    黍稷高基元EST-SSR標(biāo)記開發(fā)及200份核心種質(zhì)資源遺傳多樣性分析

    王 倩1,2張立媛3許 月4李 海5劉少雄2薛亞鵬1,2陸 平2王瑞云1,*劉敏軒2,*

    1山西農(nóng)業(yè)大學(xué)農(nóng)學(xué)院, 山西太谷 030801;2中國(guó)農(nóng)業(yè)科學(xué)院作物科學(xué)研究所, 北京 100081;3赤峰市農(nóng)牧科學(xué)研究所, 內(nèi)蒙古赤峰 024031;4吉林大學(xué)生命科學(xué)學(xué)院, 吉林長(zhǎng)春 130012;5山西農(nóng)業(yè)大學(xué)高寒區(qū)作物研究所, 山西大同 037008

    為搞清黍稷核心種質(zhì)資源的遺傳背景, 開發(fā)新的高基元EST-SSR標(biāo)記用以評(píng)估國(guó)內(nèi)不同生態(tài)區(qū)資源的遺傳多樣性, 為加快黍稷育種進(jìn)程和挖掘優(yōu)異種質(zhì)提供依據(jù)。本研究基于轉(zhuǎn)錄組測(cè)序開發(fā)了200個(gè)高基元EST-SSR標(biāo)記并進(jìn)行多態(tài)性篩選, 利用篩選到的多態(tài)性標(biāo)記評(píng)價(jià)200份黍稷核心種質(zhì)資源的遺傳多樣性。結(jié)果表明, 200個(gè)標(biāo)記在6份地理來源差異顯著的黍稷材料中有52個(gè)呈多態(tài)性, 開發(fā)效率為26%, 其中四、五和六堿基重復(fù)多態(tài)性標(biāo)記分別為17個(gè)(32.7%)、17個(gè)(32.7%)和18個(gè)(34.6%)。利用52個(gè)高基元EST-SSR多態(tài)性標(biāo)記對(duì)200份核心資源進(jìn)行檢測(cè), 共得到129個(gè)觀測(cè)等位變異, 每個(gè)位點(diǎn)檢測(cè)到2~3個(gè); Shannon多樣性指數(shù)為0.3093~1.0910, 平均為0.7277; 多態(tài)性信息含量為0.1948~0.8211, 平均為0.5104?;赨PGMA將200份資源劃分為6個(gè)群組?;赟tructure的遺傳結(jié)構(gòu)(=7)將資源劃分為7個(gè)類群, 黍稷資源主要來自4個(gè)(東北地區(qū)、華北地區(qū)、黃土高原和北方地區(qū))基因庫(kù)?;谥鞒煞址治鰧⒉牧暇蹫?個(gè)類群, 與其地理來源一致。構(gòu)建了52個(gè)四–六堿基重復(fù)EST-SSR標(biāo)記, 開發(fā)率為26%; 用其分析材料的遺傳差異發(fā)現(xiàn)黍稷地方核心種質(zhì)遺傳多樣性較為豐富。

    黍稷; 核心種質(zhì); 高基元EST-SSR; 遺傳多樣性; 群體結(jié)構(gòu)

    黍稷(L.), 又叫糜子、糜黍和粢, 是起源于中國(guó)的一年生暖季型作物, 栽培歷史超過10,000年[1-3]。黍稷生長(zhǎng)周期短, 具有耐旱、耐瘠、耐鹽堿和水資源利用效率高的特性, 適應(yīng)不良環(huán)境的能力強(qiáng), 在世界范圍內(nèi)廣泛種植[4-6]。黍稷富含蛋白質(zhì)、礦物質(zhì)、維生素和微量元素, 包括銅、鋅、鐵、錳[7], 其籽粒中的必需氨基酸(亮氨酸、異亮氨酸、蛋氨酸)含量高于小麥[8]。黍稷脫殼后稱黃米, 可制作粘糕、釀制黃酒。目前中國(guó)保存黍稷資源9800余份[9], 對(duì)其進(jìn)行遺傳多樣性研究, 有助于探究其起源與進(jìn)化關(guān)系, 同時(shí)為良種選育、資源保護(hù)等提供材料基礎(chǔ)。

    評(píng)估種質(zhì)資源遺傳多樣性可以通過表型性狀或者分子標(biāo)記等方法[10]。表型性狀具有不穩(wěn)定性, 容易受到外界環(huán)境的影響, 從而導(dǎo)致評(píng)價(jià)結(jié)果不準(zhǔn)確[11]; 而分子標(biāo)記技術(shù)是從DNA水平上揭示遺傳變異, 受外界環(huán)境影響小, 結(jié)果穩(wěn)定可靠[12]。目前廣泛應(yīng)用于黍稷遺傳差異分析的分子標(biāo)記有SSR。根據(jù)SSR的來源可將其分為基因組SSR和EST-SSR (表達(dá)序列標(biāo)簽, expressed sequence tag, EST) 2個(gè)類型[13]。基因組SSR已廣泛用于黍稷遺傳多樣性分析、遺傳圖譜構(gòu)建等領(lǐng)域[1,6,14-15]。連帥等[16]用63個(gè)基因組SSR分析國(guó)內(nèi)外192份黍稷地方品種和野生種質(zhì)遺傳多樣性, 檢測(cè)出161個(gè)等位變異, 發(fā)現(xiàn)不同地理來源群體間表現(xiàn)出顯著的遺傳多樣性差異, 其中內(nèi)蒙地區(qū)、東北地區(qū)和黃土高原地區(qū)種質(zhì)資源遺傳多樣性最豐富。Liu等[1]用67個(gè)基因組SSR分析88份黍稷種質(zhì), 檢測(cè)到179個(gè)等位基因, 多態(tài)性信息含量平均值為0.376, 不同生態(tài)區(qū)內(nèi)遺傳變異豐富, 而且不同群體間遺傳關(guān)系復(fù)雜。薛延桃等[17]用103個(gè)基因組SSR評(píng)估216份甘肅和新疆黍稷種質(zhì), 共檢測(cè)到299個(gè)等位基因, 發(fā)現(xiàn)甘肅資源遺傳多樣性較高。EST-SSR (表達(dá)序列標(biāo)簽, expressed sequence tag, EST)標(biāo)記是基于轉(zhuǎn)錄組測(cè)序技術(shù)開發(fā)得出的, 具有重復(fù)性好、多態(tài)性高、穩(wěn)定性強(qiáng)、信息量大等特點(diǎn)[18]。目前, 基于轉(zhuǎn)錄組測(cè)序技術(shù)挖掘了大量的黍稷EST-SSR標(biāo)記。劉笑瑜等[19]從一個(gè)抗旱和抗鹽品種轉(zhuǎn)錄組測(cè)序結(jié)果中開發(fā)了151個(gè)高基元SSR標(biāo)記, 檢測(cè)發(fā)現(xiàn)85個(gè)標(biāo)記具有多態(tài)性(四、五和六堿基重復(fù)標(biāo)記分別為71個(gè)、10個(gè)和4個(gè)), 開發(fā)效率為56.29%。王璐琳等[20]根據(jù)實(shí)驗(yàn)室前期黍稷轉(zhuǎn)錄組測(cè)序結(jié)果設(shè)計(jì)了70個(gè)三堿基重復(fù)SSR標(biāo)記, 篩選出17個(gè)多態(tài)性標(biāo)記, 開發(fā)效率為25.37%。Jiang等[21]從雁黍5號(hào)中發(fā)掘了8139個(gè)EST-SSR, 隨機(jī)驗(yàn)證24個(gè)標(biāo)記發(fā)現(xiàn), 3個(gè)具有多態(tài)性, 開發(fā)效率為12.5%。陳小紅等[22]基于黃糜子和鎮(zhèn)原大糜子轉(zhuǎn)錄組測(cè)序結(jié)果構(gòu)建了200個(gè)低基元SSR, 驗(yàn)證發(fā)現(xiàn)80個(gè)標(biāo)記具有多態(tài)性(一、二和三堿基重復(fù)SSR標(biāo)記分別為10個(gè)、15個(gè)和55個(gè)), 開發(fā)效率為40%。高基元SSR引物分辨率高, 對(duì)不同基因型分辨力強(qiáng), 可以檢測(cè)到低基元SSR擴(kuò)增不出的變異[19]。何杰麗等[23]用80個(gè)低基元EST-SSR標(biāo)記對(duì)來自國(guó)內(nèi)外的144份黍稷種質(zhì)進(jìn)行遺傳差異評(píng)估, 檢測(cè)到206個(gè)等位變異, PIC平均值為0.4573, 基于主成分分析將材料分為7個(gè)類群, 該結(jié)果與材料的地理來源一致。王舒婷等[24]用27個(gè)高基元SSR擴(kuò)增57份山西黍稷種質(zhì), 檢測(cè)到71個(gè)等位變異, 發(fā)現(xiàn)山西省黍稷資源的遺傳多樣性豐富。王瑞云等[25]用85個(gè)高基元SSR評(píng)估國(guó)內(nèi)6個(gè)生態(tài)栽培區(qū)的96份黍稷資源, 檢測(cè)到232個(gè)等位變異, 基于UPGMA和遺傳結(jié)構(gòu)均將材料劃分為4個(gè)類群。陳小紅等[26]用30個(gè)高基元SSR標(biāo)記評(píng)估國(guó)內(nèi)4個(gè)生態(tài)栽培區(qū)的133份種質(zhì), 檢測(cè)到90個(gè)等位變異, 基于UPGMA將試材聚為3個(gè)類群, 主成分分析將材料劃歸4類。

    核心種質(zhì)是以較小的種質(zhì)資源最大限度地代表整個(gè)種質(zhì)資源的遺傳多樣性, 因此對(duì)于核心種質(zhì)遺傳多樣性的評(píng)估極為重要[27]。隨著分子生物學(xué)的發(fā)展, 分子標(biāo)記在核心種質(zhì)遺傳多樣性檢測(cè)方面應(yīng)用也較為廣泛。姜俊燁等[27]用24個(gè)SSR分析了國(guó)內(nèi)外1075份蠶豆核心種質(zhì)的遺傳多樣性。車卓等[28]用13個(gè)SRAP標(biāo)記分析了100份黑芝麻核心種質(zhì)遺傳多樣性。徐彥剛等[29]用45個(gè)SSR對(duì)不同來源的78份西瓜核心種質(zhì)進(jìn)行了遺傳多樣性分析。種質(zhì)資源是創(chuàng)造新性狀、培育新品種的材料基礎(chǔ), 課題組前期構(gòu)建了516份國(guó)內(nèi)外黍稷核心種質(zhì), 含有國(guó)內(nèi)外地方品種、育成品種和野生材料, 從中篩選出來自18個(gè)省份的200份具有代表性的優(yōu)良性狀的黍稷地方核心資源, 以促進(jìn)黍稷種質(zhì)資源的進(jìn)一步鑒定和有效利用。

    黍稷雖然是一種小宗作物, 但其栽培歷史悠久、種質(zhì)資源豐富, 且目前已完成了全基因組測(cè)序工作[30-31], 并廣泛開發(fā)了SSR等分子標(biāo)記, 從分子水平開展了較多遺傳多樣性和群體遺傳結(jié)構(gòu)研究。但所用遺傳標(biāo)記多為低基元SSR, 對(duì)研究材料的基因型分辨能力不足, 而開發(fā)適用于黍稷研究的高基元多態(tài)性SSR標(biāo)記并利用其進(jìn)行遺傳多樣性分析的相關(guān)研究較少。課題組前期對(duì)抗旱和干旱敏感品種進(jìn)行了轉(zhuǎn)錄組測(cè)序, 發(fā)現(xiàn)存在大量高基元EST-SSR標(biāo)記, 其中五、六堿基重復(fù)標(biāo)記較多, 因此, 本研究開發(fā)了200個(gè)新的高基元SSR標(biāo)記。并且前期構(gòu)建了一批黍稷核心種質(zhì)資源, 可以大大提高黍稷種質(zhì)資源利用效率?;诖? 本研究擬首先利用6份地理來源差異顯著的材料對(duì)本試驗(yàn)開發(fā)的200個(gè)高基元SSR標(biāo)記進(jìn)行多態(tài)性初篩, 隨后利用篩選到的高基元SSR標(biāo)記對(duì)構(gòu)建的200份黍稷核心種質(zhì)進(jìn)行遺傳多樣性分析, 為鑒定黍稷優(yōu)異種質(zhì)、育種親本選配提供理論依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)材料

    參試材料為中國(guó)農(nóng)業(yè)科學(xué)院作物科學(xué)研究所資源中心小宗作物創(chuàng)新組所構(gòu)建的黍稷地方品種核心種質(zhì)200份, 來源于國(guó)內(nèi)6個(gè)生態(tài)栽培區(qū)(表1和附表1)。

    表1 參試黍稷資源的生態(tài)區(qū)分布

    其他糜子區(qū): 南方秋冬糜子區(qū)和青藏高原春糜子區(qū)。

    Other sowing ecotope: the Southern autumn & winter-sowing ecotope and the Qinghai-Tibet plateau spring-sowing ecotope.

    1.2 試驗(yàn)方法

    1.2.1 DNA提取與檢測(cè) 選取15 d左右的幼葉, 用DNA提取試劑盒(北京鼎國(guó)昌盛生物有限公司)提取全基因組DNA。用1%的瓊脂糖凝膠電泳檢測(cè)DNA的質(zhì)量, 用超微量分光光度計(jì)測(cè)定其濃度。調(diào)試其終濃度至30~60 ng μL–1, 保存于–20℃?zhèn)溆谩?/p>

    1.2.2 EST-SSR標(biāo)記分析 用地理來源差異較大的6份材料(表2)對(duì)前期通過轉(zhuǎn)錄組測(cè)序獲得的200對(duì)高基元EST-SSR引物進(jìn)行多態(tài)性篩選, 其中52對(duì)引物(附表2, 由英濰捷基(上海)貿(mào)易有限公司合成)的擴(kuò)增條帶清晰且多態(tài)性穩(wěn)定, 可用于下一步評(píng)估200份黍稷核心種質(zhì)資源的遺傳多樣性。

    表2 用于篩選EST-SSR引物的黍稷材料

    在MY-CYCLER PCR儀上進(jìn)行PCR擴(kuò)增反應(yīng)。PCR反應(yīng)體系(20 μL)包括PCR Master Mix (北京擎科生物技術(shù)有限公司) 10 μL、正向引物和反向引物各0.8 μL、基因組DNA 1 μL和ddH2O 7.4 μL。PCR反應(yīng)程序?yàn)?8℃ 5 min; 98℃ 10 s,m值10 s, 72℃ 15 s, 36個(gè)循環(huán); 72℃ 2 min, 4℃保存。用8%聚丙烯酰胺凝膠電泳分離PCR產(chǎn)物, 硝酸銀染色檢測(cè)擴(kuò)增產(chǎn)物。

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

    讀取電泳結(jié)果, 同一位置有擴(kuò)增條帶的記為1, 否則記為0。用PopGen 1.32計(jì)算每對(duì)引物的遺傳多樣性參數(shù), 包括觀測(cè)等位基因()、有效等位基因()、Shannon多樣性指數(shù)()等。用PowerMarker 3.25計(jì)算多態(tài)性信息含量(PIC)和遺傳距離。用MEGA 6.06構(gòu)建聚類圖。分別用NTSYSpc 2.11和Structure 2.3.4進(jìn)行主成分和群體遺傳結(jié)構(gòu)分析。

    2 結(jié)果與分析

    2.1 高基元EST-SSR引物篩選

    利用6份地理來源差異顯著的黍稷對(duì)200個(gè)高基元EST-SSR標(biāo)記進(jìn)行多態(tài)性擴(kuò)增發(fā)現(xiàn), 14對(duì)引物無擴(kuò)增條帶, 186對(duì)引物有擴(kuò)增條帶, 其中52對(duì)在6份材料間呈多態(tài)性,開發(fā)效率為26%。52對(duì)多態(tài)性引物中四、五和六堿基重復(fù)引物分別有17對(duì)(32.7%)、17對(duì)(32.7%)和18對(duì)(34.6%)。

    2.2 基于52對(duì)高基元EST-SSR標(biāo)記的200份黍稷資源遺傳多樣性分析

    利用52對(duì)高基元EST-SSR引物評(píng)估200份黍稷核心種質(zhì)資源的遺傳多樣性(附表3)發(fā)現(xiàn), 52個(gè)位點(diǎn)共檢測(cè)出129個(gè)觀測(cè)等位變異, 每個(gè)位點(diǎn)檢測(cè)到2~3個(gè), 平均為2.4804個(gè); 有效等位變異變化范圍為1.1893 (BMGJY80)~ 2.9565 (BMGJY141), 平均為1.9729; Shannon多樣性指數(shù)為0.3093 (BMGJY155)~1.0910 (BMGJY141), 平均為0.7277; 觀測(cè)雜合度為0 (BMGJY5、BMGJY12和BMGJY19等)~0.9282 (BMGJY177), 平均為0.2695; 期望雜合度為0.1597 (BMGJY80)~0.6636 (BMGJY141), 平均為0.4609; Nei’s期望雜合度為0.1592 (BMGJY80)~0.6618 (BMGJY141), 平均為0.4593; 多態(tài)性信息含量為0.1948 (BMGJY86)~0.8211 (BMGJY46), 平均為0.5104。

    2.3 不同生態(tài)區(qū)黍稷資源的遺傳多樣性和相似性分析

    分析不同生態(tài)區(qū)黍稷資源的遺傳多樣性參數(shù)(表3)發(fā)現(xiàn), 各生態(tài)區(qū)觀測(cè)等位基因變化范圍為2.1346~2.4808, 最高的是北方春糜子區(qū)和黃土高原春夏糜子區(qū), 最低的是其他糜子區(qū)。有效等位基因變化范圍為1.7636~1.9570, 最高的是東北春糜子區(qū), 最低的是華北夏糜子區(qū)。從Nei’s基因雜合度分析, 觀測(cè)雜合度最高的是黃土高原春夏糜子區(qū)(0.2803), 最低的是華北夏糜子區(qū)(0.2272)。期望雜合度最高的是東北春糜子區(qū)(0.4558), 最低的是華北夏糜子區(qū)(0.3715)。就Shannon多樣性指數(shù)和PIC值而言, 黃土高原春夏糜子區(qū)材料最高, 其他糜子區(qū)材料最低, 說明前者遺傳多樣性比較豐富。

    分析不同生態(tài)區(qū)黍稷資源的遺傳距離與遺傳一致度參數(shù)(表4)發(fā)現(xiàn), 遺傳距離變化范圍為0.0177~0.1238, 平均值為0.0670; 遺傳一致度變化范圍為0.8835~0.9824, 平均值為0.9357。北方春糜子區(qū)和黃土高原春夏糜子區(qū)間的遺傳距離最小(0.0177), 遺傳一致度最高(0.9824), 東北春糜子區(qū)和其他糜子區(qū)間的遺傳距離最大(0.1238), 遺傳一致度最低(0.8835)。由此表明地理位置越近, 遺傳距離越小, 遺傳相似性越高, 親緣關(guān)系越近; 地理位置越遠(yuǎn), 遺傳距離越大, 遺傳相似性越低, 親緣關(guān)系越遠(yuǎn)。

    2.4 基于UPGMA的黍稷資源聚類分析

    基于UPGMA對(duì)200份黍稷核心種質(zhì)進(jìn)行聚類分析(圖1)發(fā)現(xiàn), 200份黍稷核心種質(zhì)分為6個(gè)組群(表5)。組群I有6份材料, 主要來自青海, 屬于黃土高原春夏糜子區(qū)。組群II有8份材料, 主要來自新疆, 屬于西北春夏糜子區(qū)。組群III有18份, 主要來自內(nèi)蒙古, 屬于北方春糜子區(qū)。組群IV有33份, 主要來自甘肅, 屬于黃土高原春夏糜子區(qū)。組群V有112份, 主要來自山西, 屬于北方春糜子區(qū)。組群VI有23份, 主要來自山東, 屬于華北夏糜子區(qū)。

    表3 6個(gè)生態(tài)區(qū)黍稷的遺傳多樣性參數(shù)

    其他糜子區(qū): 南方秋冬糜子區(qū)和青藏高原春糜子區(qū)。

    Other sowing ecotope: the Southern autumn & winter-sowing ecotope and the Qinghai-Tibet plateau spring-sowing ecotope. NES: Northeast spring-sowing ecotope; NSU: Northern summer-sowing ecotope; NSP: Northern spring-sowing ecotope; LPSS: Loess plateau spring & summer-sowing ecotope; NWSS: Northwest spring & summer-sowing ecotope; OSP: other sowing ecotope.

    表4 不同黍稷群體間的遺傳距離與遺傳一致度

    Nei氏遺傳一致度(對(duì)角線以上)和遺傳距離(對(duì)角線以下)。其他糜子區(qū): 南方秋冬糜子區(qū)和青藏高原春糜子區(qū)。

    Nei’s genetic identity (above diagonal) and genetic distance (below diagonal). NES: Northeast spring-sowing ecotope; NSU: Northern summer-sowing ecotope; NSP: Northern spring-sowing ecotope; LPSS: Loess plateau spring & summer-sowing ecotope; NWSS: Northwest spring & summer-sowing ecotope; OSP: other sowing ecotope. Other sowing ecotope: the Southern autumn & winter-sowing ecotope and the Qinghai-Tibet plateau spring-sowing ecotope.

    表5 基于UPGMA的200份黍稷資源組群分布

    圖1 基于UPGMA的200份黍稷資源聚類分析

    NES: 東北春糜子區(qū); NSU: 華北夏糜子區(qū); NSP: 北方春糜子區(qū); LPSS: 黃土高原春夏糜子區(qū); NWSS: 西北春夏糜子區(qū); OSP: 其他糜子區(qū)。其他糜子區(qū): 南方秋冬糜子區(qū)和青藏高原春糜子區(qū)。

    NES: Northeast spring-sowing ecotope; NSU: Northern summer-sowing ecotope; NSP: Northern spring-sowing ecotope; LPSS: Loess plateau spring & summer-sowing ecotope; NWSS: Northwest spring & summer-sowing ecotope; OSP: other sowing ecotope. Other sowing ecotope:the Southern autumn & winter-sowing ecotope and the Qinghai-Tibet plateau spring-sowing ecotope.

    2.5 基于Structure的群體遺傳結(jié)構(gòu)分析

    用Structure軟件對(duì)200份供試材料基因庫(kù)數(shù)目進(jìn)行建模(圖2)發(fā)現(xiàn), ?在=2和=7處存在峰值。由圖3可知,=2時(shí), 將試材分為2個(gè)類群, 紅色類群(121份), 主要為山西和陜西種質(zhì), 屬于北方春和黃土高原春夏糜子區(qū); 綠色類群(79份), 主要為甘肅和河北種質(zhì), 屬于北方春和黃土高原春夏糜子區(qū)。=7時(shí), 劃分為7個(gè)類群, 紅色類群(32份), 主要來自山西和陜西種質(zhì), 代表黃土高原基因庫(kù); 綠色類群(26份), 主要來自青海種質(zhì), 代表黃土高原基因庫(kù); 深藍(lán)色類群(31份), 主要來自黑龍江和吉林種質(zhì), 代表東北基因庫(kù); 黃色類群(19份), 主要來自山東種質(zhì), 代表華北基因庫(kù); 紫色類群(25份), 主要來自山西種質(zhì), 代表北方基因庫(kù); 淺藍(lán)色類群(34份), 主要來自山西種質(zhì), 代表北方基因庫(kù); 橙色類群(33份), 主要來自甘肅種質(zhì), 代表黃土高原基因庫(kù)。

    圖2 基于Structure對(duì)200份黍稷資源群體建模

    圖3 基于Structure的黍稷核心種質(zhì)資源群體遺傳結(jié)構(gòu)圖

    顏色代表類群。橫坐標(biāo)數(shù)字分別代表黍稷材料序號(hào)和生態(tài)區(qū)類型。NES: 東北春糜子區(qū); NSU: 華北夏糜子區(qū); NSP: 北方春糜子區(qū); LPSS: 黃土高原春夏糜子區(qū); NWSS: 西北春夏糜子區(qū); OSP: 其他糜子區(qū)。其他糜子區(qū): 南方秋冬糜子區(qū)和青藏高原春糜子區(qū)。

    Different color represents the different group. Numbers in the horizontal axis represent accession number and type of ecotope area, respectively.NES: Northeast spring-sowing ecotope; NSU: Northern summer-sowing ecotope; NSP: Northern spring-sowing ecotope; LPSS: Loess plateau spring & summer-sowing ecotope; NWSS: Northwest spring & summer-sowing ecotope; OSP: other sowing ecotope. Other sowing ecotope:the Southern autumn & winter-sowing ecotope and the Qinghai-Tibet plateau spring-sowing ecotope.

    2.6 200份黍稷資源遺傳差異的主成分分析

    對(duì)200份供試材料進(jìn)行三維PCA分析(圖4)發(fā)現(xiàn), 前3個(gè)主成分PC1、PC2和PC3分別解釋總變異的49.7%、5.4%和3.5%, 累積解釋變異58.6%。200份黍稷資源可以劃分為6個(gè)類群(東北春糜子區(qū)28份、華北夏糜子區(qū)23份、北方春糜子區(qū)63份、黃土高原春夏糜子區(qū)67份、西北春夏糜子區(qū)12份、其他糜子區(qū)7份), 與材料的地理來源一致。

    圖4 200份黍稷資源遺傳多樣性的主成分分析

    NES: 東北春糜子區(qū); NSU: 華北夏糜子區(qū); NSP: 北方春糜子區(qū); LPSS: 黃土高原春夏糜子區(qū); NWSS: 西北春夏糜子區(qū); OSP: 其他糜子區(qū)。其他糜子區(qū): 南方秋冬糜子區(qū)和青藏高原春糜子區(qū)。

    NES: Northeast spring-sowing ecotope; NSU: Northern summer-sowing ecotope; NSP: Northern spring-sowing ecotope; LPSS:Loess plateau spring & summer-sowing ecotope; NWSS: Northwest spring & summer-sowing ecotope; OSP: other sowing ecotope. Other sowing ecotope: the Southern autumn & winter-sowing ecotope and the Qinghai-Tibet plateau spring-sowing ecotope.

    3 討論

    3.1 黍稷高基元SSR標(biāo)記的開發(fā)利用

    黍稷是世界上最早馴化的谷類作物之一, 具有重要的農(nóng)業(yè)、歷史和進(jìn)化意義。黍稷種質(zhì)資源豐富, 因此開發(fā)新的分子標(biāo)記對(duì)于促進(jìn)其遺傳多樣性的研究意義重大。高基元SSR標(biāo)記的引物分辨率高, 對(duì)不同基因型的分辨能力強(qiáng), 2017年, 劉笑瑜[19]用151個(gè)高基元SSR擴(kuò)增6份黍稷材料, 發(fā)現(xiàn)85個(gè)標(biāo)記呈多態(tài)性, 多態(tài)率為56.3%。本研究從200對(duì)通過轉(zhuǎn)錄組測(cè)序得到的高基元SSR引物中篩選得到52對(duì)多態(tài)性引物, 多態(tài)率為26%, 低于劉笑瑜的結(jié)果, 這可能與引物篩選標(biāo)準(zhǔn)不同以及獲得引物渠道有關(guān)。本研究開發(fā)的四、五和六堿基重復(fù)多態(tài)性引物分別為17對(duì)、17對(duì)和18對(duì), 而劉笑瑜的研究結(jié)果分別為71對(duì)、10對(duì)和4對(duì), 表明轉(zhuǎn)錄組測(cè)序材料不同, 獲得的SSR標(biāo)記存在差異, 本研究開發(fā)的五、六堿基SSR標(biāo)記數(shù)量占比顯著高于劉笑瑜的結(jié)果。一般來說, PIC<0.25, 為低度多態(tài)性引物, 0.25≤PIC<0.5為中度多態(tài)性引物, PIC≥0.5, 為高度多態(tài)性引物[32]。王舒婷等[24]利用27對(duì)高基元SSR檢測(cè)57份黍稷種質(zhì)的PIC平均值為0.5667, 陳小紅等[26]用30對(duì)高基元SSR擴(kuò)增130份黍稷材料發(fā)現(xiàn)PIC平均值為0.6966。本研究與前人結(jié)果一致, PIC均大于0.5, 所選引物為高度多態(tài)性引物, 說明都可用于黍稷遺傳多樣性的研究。

    3.2 黍稷種質(zhì)資源的遺傳差異

    本研究參試200份黍稷地方核心種質(zhì)來源于國(guó)內(nèi)6個(gè)生態(tài)區(qū)的18個(gè)省、市、自治區(qū), 具有很強(qiáng)的代表性。因此, 利用新開發(fā)的高基元SSR分子標(biāo)記對(duì)200份黍稷核心種質(zhì)進(jìn)行遺傳多樣性分析, 將為黍稷種質(zhì)資源的收集保存以及育種研究提供理論依據(jù)。

    本研究評(píng)估200份黍稷核心種質(zhì)資源遺傳差異發(fā)現(xiàn)Shannon多樣性指數(shù)()變化范圍為0.3093~1.0910, 平均為0.7277, 低于王舒婷等[24](0.7686)和王瑞云等[25](0.7708)結(jié)果, 這與選取試驗(yàn)材料差異有關(guān), 前人所選試材中含有農(nóng)家種和育成品種, 材料來源更廣, 而本研究材料均為地方品種; 高于連帥等[16](0.6275)結(jié)果, 這可能與我們所用高基元EST-SSR標(biāo)記有關(guān), 該標(biāo)記可以檢測(cè)到的變異更多。

    有效等位基因數(shù)、Shannon多樣性指數(shù)()和多態(tài)信息含量(PIC)可以綜合的反映群體遺傳多樣性水平[33]。從生態(tài)區(qū)角度對(duì)黍稷資源進(jìn)行遺傳多樣性分析, 結(jié)果表明黃土高原春夏糜子區(qū)材料的有效等位變異、基因多樣性和PIC指標(biāo)均高于其他6個(gè)生態(tài)區(qū)材料, 顯現(xiàn)出較豐富的遺傳多樣性, 這與王瑞云等[25]的結(jié)果一致。何杰麗等[23]、石甜甜等[34]和寇淑君等[35]研究均發(fā)現(xiàn)北方春糜子區(qū)與黃土高原春夏糜子區(qū)黍稷資源間的遺傳距離最小, 遺傳一致度最高, 親緣關(guān)系較近, 本研究也發(fā)現(xiàn)類似結(jié)果, 表明不同生態(tài)區(qū)黍稷種質(zhì)資源間遺傳距離的遠(yuǎn)近與其栽培環(huán)境及地理分布密切相關(guān)。

    3.3 黍稷種質(zhì)資源的聚類分析與遺傳結(jié)構(gòu)

    聚類分析與群體結(jié)構(gòu)分析是研究作物種質(zhì)資源遺傳多樣性和遺傳背景的有效手段, 是作物種質(zhì)資源有效保護(hù)和利用的前提[36-37]。董俊麗等[38]對(duì)96份黍稷資源的表型性狀進(jìn)行遺傳結(jié)構(gòu)和聚類分析結(jié)果表明, 類群劃分與地理來源有一定的相關(guān)性。王瑞云等[25]利用高基元微衛(wèi)星標(biāo)記對(duì)96份中國(guó)糜子遺傳多樣性研究發(fā)現(xiàn), 聚類分析和遺傳結(jié)構(gòu)分析結(jié)果均與地理起源相關(guān)。本研究基于遺傳距離的聚類分析和遺傳結(jié)構(gòu)分析發(fā)現(xiàn), 地理來源相近的資源基本聚在同一類群, 與前人結(jié)果一致。因此, 本研究采用高基元EST-SSR標(biāo)記方法對(duì)200份黍稷核心種質(zhì)資源的遺傳關(guān)系和群體結(jié)構(gòu)進(jìn)行分析, 明確了黍稷核心種質(zhì)資源間遺傳背景和群體間的遺傳多樣性情況, 對(duì)今后育種材料的選用具有指導(dǎo)意義。

    3.4 核心種質(zhì)的構(gòu)建

    作物種質(zhì)資源是作物新品種的選育、特異種質(zhì)材料的利用以及種質(zhì)創(chuàng)新的物質(zhì)基礎(chǔ)[39]。隨著科學(xué)技術(shù)的發(fā)展和育種水平的不斷提高, 國(guó)內(nèi)外已構(gòu)建了包括四大糧食作物水稻[40]、小麥[41]、大麥[42]、玉米[43]等多種作物的核心種質(zhì)。中國(guó)是黍稷的起源中心, 之前也有報(bào)道構(gòu)建了黍稷核心種質(zhì), 包括胡興雨等[44]以國(guó)家作物種質(zhì)庫(kù)中保存的8016份國(guó)內(nèi)外黍稷資源為材料, 在11個(gè)表型性狀聚類的基礎(chǔ)上, 按比例法取樣, 并依各組的遺傳多樣性指數(shù)進(jìn)行適當(dāng)調(diào)整, 構(gòu)建了780份國(guó)內(nèi)外黍稷核心種質(zhì); 董孔軍等[45]以國(guó)內(nèi)黍稷地方種質(zhì)資源566份為基礎(chǔ), 按照20%取樣比例確定相應(yīng)數(shù)量類群, 根據(jù)各基因型在類群中的平均離差度確定取樣樣本, 構(gòu)建了甘肅省黍稷地方核心種質(zhì)122份。課題組前期構(gòu)建了516份國(guó)內(nèi)外黍稷核心種質(zhì)資源, 其中含有國(guó)內(nèi)外地方品種、育成品種及野生材料, 本試驗(yàn)材料200份是從中選出的具有優(yōu)良性狀、代表性強(qiáng)的國(guó)內(nèi)地方品種。因此, 本研究利用新開發(fā)的高基元SSR標(biāo)記技術(shù)評(píng)估200份黍稷核心種質(zhì)資源的遺傳多樣性, 對(duì)于了解黍稷親緣關(guān)系、挖掘優(yōu)異種質(zhì)、提高育種效率和種質(zhì)資源利用率意義重大。

    附表 請(qǐng)見網(wǎng)絡(luò)版: 1) 本刊網(wǎng)站http://zwxb. chinacrops.org/; 2) 中國(guó)知網(wǎng)http://www.cnki.net/; 3) 萬(wàn)方數(shù)據(jù)http://c.wanfangdata.com.cn/Periodical-zuowxb.aspx。

    [1] Liu M X, Xu Y, He J H, Zhang S, Wang Y Y, Lu P. Genetic diversity and population structure of broomcorn millet (L.) cultivars and landraces in China based on microsatellite markers., 2016, 17: 370.

    [2] Lu H, Zhang J, Liu K B, Wu N, LI Y, Zhou K, Ye M, Zhang T, Zhang H, Yang X, Shen L, Xu D, Li Q. Earliest domestication of common millet (L.) in East Asia extended to 10,000 years ago., 2009, 106: 7367–7372.

    [3] Wang R Y, Hunt H V, Qiao Z J, Wang L, Han Y H. Diversity and cultivation of broomcorn millet (L.) in China: a review., 2016, 70: 332–342.

    [4] 董明, 降彥苗, 李海權(quán), 耿玲玲, 劉建燁, 喬志紅, 劉國(guó)慶. 光周期變化對(duì)糜子形態(tài)建成及幼穗發(fā)育進(jìn)程的影響. 中國(guó)農(nóng)業(yè)科學(xué), 2020, 53: 1118–1125. Dong M, Jiang Y M, Li H Q, Geng L L, Liu J Y, Qiao Z H, Liu G Q. Effects of photoperiod changes on morphological characters and young panicle development in proso millet (L.)., 2020, 53: 1118–1125 (in Chinese with English abstract).

    [5] Amadoubr I, Le M. Millets: nutritional composition, some health benefits and processing: a review., 2013, 25: 501.

    [6] 薛延桃, 陸平, 喬治軍, 劉敏軒, 王瑞云. 基于SSR標(biāo)記的黍稷種質(zhì)資源遺傳多樣性及親緣關(guān)系研究. 中國(guó)農(nóng)業(yè)科學(xué), 2018, 51: 2846–2859. Xue Y T, Lu P, Qiao Z J, Liu M X, Wang R Y. Genetic diversity and genetic relationship of broomcorn millet (L.) germplasm based on SSR markers., 2018, 51: 2846–2859 (in Chinese with English abstract).

    [7] Boukail S, Macharia M, Miculan M, Masoni A, Calamai A, Palchetti E, Dell’Acqua M. Genome wide association study of agronomic and seed traits in a world collection of proso millet (L.)., 2021, 21: 330.

    [8] Kalinova J, Moudry J. Content and quality of protein in proso millet (L.) varieties., 2006, 61: 43–47.

    [9] 王瑞云. 糜子遺傳多樣性及進(jìn)化研究進(jìn)展. 北京: 中國(guó)農(nóng)業(yè)出版社, 2017. pp 2–3. Wang R Y. Genetic Diversity and Evolution Advancement in Common Millet (L.). Beijing: China Agriculture Press, 2017. pp 2–3 (in Chinese with English abstract).

    [10] 宋蕓, 張?chǎng)稳? 李政, 孫哲, 李澳旋, 杜曉蓉, 喬永剛. 基于SSR分子標(biāo)記的柴胡遺傳多樣性與遺傳結(jié)構(gòu)分析. 藥學(xué)學(xué)報(bào), 2022, 57: 1193–1202. Song Y, Zhang X R, Li Z, Sun Z, Li A X, Du X R, Qiao Y G. Genetic diversity and genetic structure analysis ofDC. based on SSR molecular marker., 2022, 57: 1193–1202 (in Chinese with English abstract).

    [11] 葛淑俊, 孟義江, 李廣敏, 馬峙英. 我國(guó)藥用植物遺傳多樣性研究進(jìn)展. 中草藥, 2006, (10): 1584–1589. Ge S J, Meng Y J, Li G M, Ma S Y. Research progress on genetic diversity in Chinese medicinal plants., 2006, (10): 1584–1589 (in Chinese with English abstract).

    [12] Grover A, Sharma P C. Development and use of molecular markers: past and present., 2016, 36: 290–302.

    [13] 程小毛, 黃曉霞. SSR標(biāo)記開發(fā)及其在植物中的應(yīng)用. 中國(guó)農(nóng)學(xué)通報(bào), 2011, 27(5): 304–307. Cheng X M, Huang X X. Development and application of SSR markers in plants., 2011, 27(5): 304–307 (in Chinese with English abstract).

    [14] 王銀月, 劉敏軒, 陸平, 喬治軍, 楊天育, 李海, 崔喜艷. 構(gòu)建黍稷分子遺傳圖譜SSR引物的篩選. 作物雜志, 2014, (4): 32–38. Wang Y Y, Liu M X, Lu P, Qiao Z J, Yang T Y, Li H, Cui X Y. The SSR marker selection of broomcorn millet (L.) for construction of genetic linkage map., 2014, (4): 32–38 (in Chinese with English abstract).

    [15] 王銀月. 黍稷種質(zhì)遺傳資源多樣性及遺傳圖譜構(gòu)建研究. 吉林農(nóng)業(yè)大學(xué)碩士學(xué)位論文, 吉林長(zhǎng)春, 2015. Wang Y Y. Assessment of Genetic Diversity and Construction of Genetic Linkage Map of Broomcorn Millet (L.). MS Thesis of Jilin Agricultural University, Changchun, Jilin, China, 2015 (in Chinese with English abstract).

    [16] 連帥, 陸平, 喬治軍, 張琦, 張茜, 劉敏軒, 王瑞云. 利用SSR分子標(biāo)記研究國(guó)內(nèi)外黍稷地方品種和野生資源的遺傳多樣性. 中國(guó)農(nóng)業(yè)科學(xué), 2016, 49: 3264–3275. Lian S, Lu P, Qiao Z J, Zhang Q, Zhang Q, Liu M X, Wang R Y. Genetic diversity in broomcorn millet (L.) from China and abroad by using SSR markers., 2016, 49: 3264–3275 (in Chinese with English abstract).

    [17] 薛延桃, 陸平, 史夢(mèng)莎, 孫昊月, 劉敏軒, 王瑞云. 新疆、甘肅黍稷資源的遺傳多樣性與群體遺傳結(jié)構(gòu)研究. 作物學(xué)報(bào), 2019, 45: 1511–1521. Xue Y T, Lu P, Shi M S, Sun H Y, Liu M X, Wang R Y. Genetic diversity and population genetic structure of broomcorn millet accessions in Xinjiang and Gansu., 2019, 45: 1511–1521 (in Chinese with English abstract).

    [18] 馬杰, 屈雯, 陳春艷, 王磊, 馬俊, 劉針杉, 馬維, 周平, 何遠(yuǎn)寬, 孫勃. 基于轉(zhuǎn)錄組序列的葉用芥菜奶奶青菜EST-SSR標(biāo)記開發(fā)與遺傳多樣性分析. 浙江農(nóng)業(yè)學(xué)報(bào), 2021, 33: 1640–1649. Ma J, Qu W, Chen C Y, Wang L, Ma J, Liu Z S, Ma W, Zhou P, He Y K, Sun B. Development of SSR markers based on transcriptome sequencing and genetic diversity analysis of Nainaiqingcai leaf mustard., 2021, 33: 1640–1649 (in Chinese with English abstract).

    [19] 劉笑瑜. 利用高基元SSR分析中國(guó)糜子資源的遺傳多樣性. 山西農(nóng)業(yè)大學(xué)碩士學(xué)位論文, 山西太谷, 2017. Liu X Y. Analysis of Genetic Diversity of Chinese Proso Millet Resources by High Motif SSR. MS Thesis of Shanxi Agricultural University, Taigu, Shanxi, China, 2017 (in Chinese with English abstract).

    [20] 王璐琳, 王瑞云, 何杰麗, 薛延桃, 陳凌, 王海崗, 喬治軍. 糜子特異性 SSR 標(biāo)記的開發(fā). 山西農(nóng)業(yè)科學(xué), 2018, 46(1): 1–4. Wang L L, Wang R Y, He J L, Xue Y T, Chen L, Wang H G, Qiao Z J. Development of specific SSR marker in common millet., 2018, 46(1): 1–4 (in Chinese with English abstract).

    [21] Jiang Y M, Li H Q, Zhang J N, Xiang J Y, Cheng R H, Liu G Q. Whole genomic EST-SSR development based on high-throughput transcript sequencing in proso millet (L.)., 2018, 20: 617–620.

    [22] 陳小紅, 何杰麗, 石甜甜, 邵歡歡, 王海崗, 陳凌, 高志軍, 王瑞云, 喬治軍. 基于轉(zhuǎn)錄組測(cè)序開發(fā)糜子SSR標(biāo)記. 中國(guó)農(nóng)業(yè)科學(xué), 2020, 53: 1940–1949. Chen X H, He J L, Shi T T, Shao H H, Wang H G, Chen L, Gao Z J, Wang R Y, Qiao Z J. Developing SSR markers of proso millet based on transcriptome sequencing., 2020, 53: 1940–1949 (in Chinese with English abstract).

    [23] 何杰麗, 石甜甜, 陳凌, 王海崗, 高志軍, 楊美紅, 王瑞云, 喬治軍. 糜子EST-SSR的開發(fā)及種質(zhì)資源遺傳多樣性分析. 植物學(xué)報(bào), 2019, 54: 723–732. He J L, Shi T T, Chen L, Wang H G, Gao Z J, Yang M H, Wang R Y, Qiao Z J. The genetic diversity of common millet () germplasm resources based on the EST-SSR markers., 2019, 54: 723–732 (in Chinese with English abstract).

    [24] 王舒婷, 何杰麗, 石甜甜, 陳凌, 王海崗, 王瑞云, 喬治軍. 利用微衛(wèi)星標(biāo)記分析山西糜子的遺傳多樣性. 植物遺傳資源學(xué)報(bào), 2019, 20: 69–78. Wang S T, He J L, Shi T T, Chen L, Wang H G, Wang R Y, Qiao Z J. Genetic diversity analysis of broomcorn millet (L.) of Shanxi province using microsatellite markers., 2019, 20: 69–78 (in Chinese with English abstract).

    [25] 王瑞云, 劉笑瑜, 王海崗, 陸平, 劉敏軒, 陳凌, 喬治軍. 用高基元微衛(wèi)星標(biāo)記分析中國(guó)糜子遺傳多樣性. 中國(guó)農(nóng)業(yè)科學(xué), 2017, 50: 3848–3859. Wang R Y, Liu X Y, Wang H G, Lu P, Liu M X, Chen L, Qiao Z J. Evaluation of genetic diversity of common millet () germplasm available in China using high motif nucleotide repeat SSR markers., 2017, 50: 3848–3859 (in Chinese with English abstract).

    [26] 陳小紅, 林元香, 王倩, 丁敏, 王海崗, 陳凌, 高志軍, 王瑞云,喬治軍. 基于高基元SSR構(gòu)建黍稷種質(zhì)資源的分子身份證. 作物學(xué)報(bào), 2022, 48: 908–919. Chen X H, Lin Y X, Wang Q, Ding M, Wang H G, Chen L, Gao Z J, Wang R Y, Qiao Z J. Development of DNA molecular ID card in hog millet germplasm based on high motif SSR., 2022, 48: 908–919 (in Chinese with English abstract).

    [27] 姜俊燁, 楊濤, 王芳, 方俐, 仲偉文, 關(guān)建平, 宗緒曉. 國(guó)內(nèi)外蠶豆核心種質(zhì)SSR遺傳多樣性對(duì)比及微核心種質(zhì)構(gòu)建. 作物學(xué)報(bào), 2014, 40: 1311–1319. Jiang J Y, Yang T, Wang F, Fang L, Zhong W W, Guan J P, Zong X X. Genetic diversity analysis of germplasm resources and construction of mini-core collections forL. at home and abroad., 2014, 40: 1311–1319 (in Chinese with English abstract).

    [28] 車卓, 張艷欣, 孫建, 張秀榮, 尚勛武, 王化俊. 應(yīng)用SRAP標(biāo)記分析黑芝麻核心種質(zhì)遺傳多樣性. 作物學(xué)報(bào), 2009, 35: 1936–1941. Che Z, Zhang Y X, Sun J, Zhang X R, Shang X W, Wang H J. Genetic diversity analysis of black sesame (DC) core collection of China using SRAP markers., 2009, 35: 1936–1941 (in Chinese with English abstract).

    [29] 徐彥剛, 姚協(xié)豐, 李良俊, 徐錦華, 任潤(rùn)生, 張曼, 劉廣, 婁麗娜, 朱凌麗, 侯茜, 徐建, 羊杏平. 西瓜核心種質(zhì)遺傳多樣性分析及指紋圖譜構(gòu)建. 分子植物育種, 2021, 19: 305–317. Xu Y G, Yao X F, Li L J, Xu J H, Ren R S, Zhang M, Liu G, Lou L N, Zhu L L, Hou Q, Xu J,Yang X P. Analysis of genetic diversity and construction of DNA fingerprints for watermelon core germplasms., 2021, 19: 305–317 (in Chinese with English abstract).

    [30] Shi J, Ma X, Zhang J, Zhou Y, Liu M, Huang L, Sun S, Zhang X, Gao X, Zhan W, Li P, Wang L, Lu P, Zhao H, Song W, Lai J. Chromosome conformation capture resolved near complete genome assembly of broomcorn millet., 2019, 10: 464.

    [31] Zou C, Li L, Miki D, Li D, Tang Q, Xiao L, Rajput S, Deng P, Peng L, Jia W, Huang R, Zhang M, Sun Y, Hu J, Fu X, Schnable P S, Chang Y, Li F, Zhang H, Feng B, Zhu X, Liu R, Schnable J C, Zhu J K, Zhang H. The genome of broomcorn millet., 2019, 10: 436.

    [32] Sharma P, Mehta G, Shefail, Muthusamy S K, Singh S K, Singh G P. Development and validation of heat-responsive candidate gene and miRNA gene based SSR markers to analysis genetic diversity in wheat for heat tolerance breeding., 2021, 48: 381–393.

    [33] 何瑞超. 綠豆遺傳多樣性研究及種子萌發(fā)期耐鹽性評(píng)價(jià). 內(nèi)蒙古農(nóng)業(yè)大學(xué)碩士學(xué)位論文, 內(nèi)蒙古呼和浩特, 2021.He R C. Study on Genetic Diversity of Mung Bean and Evaluation of Salt Tolerance During Seed Germination. MS Thesis of Inner Mongolia Agricultural University, Hohhot, Inner Mongolia, China, 2021 (in Chinese with English abstract).

    [34] 石甜甜, 何杰麗, 高志軍, 陳凌, 王海崗, 喬治軍, 王瑞云. 利用EST-SSR評(píng)估糜子資源遺傳差異. 中國(guó)農(nóng)業(yè)科學(xué), 2019, 52: 4100–4109. Shi T T, He J L, Gao Z J, Chen L, Wang H G, Qiao Z J, Wang R Y. Genetic diversity of common millet resources assessed with EST-SSR markers., 2019, 52: 4100–4109 (in Chinese with English abstract).

    [35] 寇淑君, 霍阿紅, 付國(guó)慶, 紀(jì)軍建, 王瑤, 左振興, 劉敏軒, 陸平. 利用熒光SSR分析中國(guó)糜子的遺傳多樣性和群體遺傳結(jié)構(gòu). 中國(guó)農(nóng)業(yè)科學(xué), 2019, 52: 1475–1487. Kou S J, Huo A H, Fu G Q, Ji J J, Wang Y, Zuo Z X, Liu M X, Lu P. Genetic diversity and population structure of broomcorn millet in China based on fluorescently labeled SSR., 2019, 52: 1475–1487 (in Chinese with English abstract).

    [36] 陳志祥, 羅小燕, 李拴林, 吳如月, 王文強(qiáng), 丁西朋. 基于SSR標(biāo)記的木豆種質(zhì)資源遺傳多樣性與群體結(jié)構(gòu)分析. 草地學(xué)報(bào), 2021, 29: 904–911. Chen Z X, Luo X Y, Li S L, Wu R Y, Wang W Q, Ding X P. Genetic diversity and population structure analysis of the pigeonpea (L.) germplasm resources based on SSR markers., 2021, 29: 904–911 (in Chinese with English abstract).

    [37] 謝立峰, 李寧, 李燁, 姚明華. 茄子種質(zhì)遺傳多樣性及群體結(jié)構(gòu)的SRAP分析. 植物學(xué)報(bào), 2019, 54: 58–63. Xie L F, Li N, Li Y, Yao M H. Genetic diversity and population structure of eggplant () germplasm resources based on SRAP method., 2019, 54: 58–63 (in Chinese with English abstract).

    [38] 董俊麗, 王海崗, 陳凌, 王君杰, 曹曉寧, 王綸, 喬治軍. 糜子骨干種質(zhì)遺傳多樣性和遺傳結(jié)構(gòu)分析. 中國(guó)農(nóng)業(yè)科學(xué), 2015, 48: 3121–3131. Dong J L, Wang H G, Chen L, Wang J J, Cao X N, Wang L, Qiao Z J. Analysis of genetic diversity and genetic structure of broomcorn millet of the Chinese., 2015, 48: 3121–3131 (in Chinese with English abstract).

    [39] 韋祖生, 付海天, 田益農(nóng). 作物核心種質(zhì)評(píng)價(jià)與應(yīng)用研究進(jìn)展.現(xiàn)代農(nóng)業(yè)科技, 2015, (21): 9–11. Wei Z S, Fu H T, Tian Y N. Overview of evaluation and application for crops core collection., 2015, (21): 9–11 (in Chinese with English abstract).

    [40] 李自超, 張洪亮, 曾亞文, 楊忠義, 申時(shí)全, 孫傳清, 王象坤. 云南地方稻種資源核心種質(zhì)取樣方案研究. 中國(guó)農(nóng)業(yè)科學(xué), 2000, 33(5): 1–7. Li Z C, Zhang H L, Zeng Y W, Yang Z Y, Shen S Q, Sun C Q, Wang X K. Study on sampling schemes of core collection of local varieties of rice in Yunnan, China., 2000, 33(5): 1–7 (in Chinese with English abstract).

    [41] 董玉琛, 曹永生, 張學(xué)勇, 劉三才, 王蘭芬, 游光霞, 龐斌雙, 李立會(huì), 賈繼增. 中國(guó)普通小麥初選核心種質(zhì)的產(chǎn)生. 植物遺傳資源學(xué)報(bào), 2003, 4: 1–8. Dong Y C, Cao Y S, Zhang X Y, Liu S C, Wang L F, You G X, Pang B S, Li L H, Jia J Z. Establishment of candidate core collections in Chinese common wheat germplasm., 2003, 4: 1–8 (in Chinese with English abstract).

    [42] 胥婷婷, 林峰, 華為, 汪軍妹, 朱靖環(huán), 賈巧君, 尚毅, 楊建明.我國(guó)大麥選育品種核心種質(zhì)的初步構(gòu)建. 浙江農(nóng)業(yè)學(xué)報(bào), 2011, 23: 483–488. Xu T T, Lin F, Hua W, Wang J M, Zhu J H, Jia Q J, Shang Y, Yang J M. Establishment of core collection on bred cultivars of Chinese barley., 2011, 23: 483–488 (in Chinese with English abstract).

    [43] 李健英, 解睿, 李正莉, 曹宇賓, 郭俊宏. 玉米核心種質(zhì)及構(gòu)建方法研究進(jìn)展. 山西農(nóng)業(yè)科學(xué), 2010, 38(9): 81–85. Li J Y, Xie R, Li Z L, Cao Y B, Guo J H. Maize core collection and its establishment., 2010, 38(9): 81–85 (in Chinese with English abstract).

    [44] 胡興雨, 王綸, 張宗文, 陸平, 張紅生. 中國(guó)黍稷核心種質(zhì)的構(gòu)建. 中國(guó)農(nóng)業(yè)科學(xué), 2008, 41: 3489–3502. Hu X Y, Wang L, Zhang Z W, Lu P, Zhang H S. Establishment of broomcorn millet core collection in China., 2008, 41: 3489–3502 (in Chinese with English abstract).

    [45] 董孔軍, 楊天育, 何繼紅, 任瑞玉, 張磊. 甘肅省糜子地方品種資源核心種質(zhì)的構(gòu)建. 甘肅農(nóng)業(yè)科技, 2012, (7): 7–11. Dong K J, Yang T Y, He J H, Ren R Y, Zhang L. Construction of core germplasm of local variety resources of millet in Gansu province., 2012, (7): 7–11 (in Chinese with English abstract).

    High motif EST-SSR markers development and genetic diversity evaluation for 200 core germplasms in proso millet

    WANG Qian1,2, ZHANG Li-Yuan3, XU Yue4, LI Hai5, LIU Shao-Xiong2, XUE Ya-Peng1,2,LU Ping2, WANG Rui-Yun1,*, and LIU Min-Xuan2,*

    1College of Agronomy, Shanxi Agricultural University, Taigu 030801, Shanxi, China;2Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China;3Chifeng Academy of Agricultural and Animal Husbandry Sciences, Chifeng 024031, Inner Mongolia, China;4College of Life Sciences, Jilin University, Changchun 130012, Jilin, China;5Institute of High Latitude Crops, Shanxi Agricultural University, Datong 037008, Shanxi, China

    The objective of this study is to clarify the genetic background for the core germplasm resources of proso millet and to evaluate the genetic diversity from different ecological regions in China. To provide a basis for accelerating proso millet breeding process and mining excellent germplasm, a new high motif EST-SSR markers were developed. In this study, 200 high motif EST-SSR markers were developed based on transcriptome sequencing results, and used to evaluate the genetic diversity for 200 core germplasm accessions of proso millet. Results showed that 52 of the 200 markers were polymorphic in 6 proso millet accessions with significant geographical differences, with the development efficiency of 26%, including 17 tetra-nucleotide repeat SSRs (32.7%), 17 penta-nucleotide repeat ones (32.7%), and 18 hexa-nucleotide repeat ones (34.6%). Among the examined 200 accessions, a total of 129 observed alleles were detected in 52 loci, changed 2–3 alleles at each locus. The Shannon diversity index was 0.3093–1.0910, with an average of 0.7277. The polymorphism information content was 0.1948–0.8211, with an average of 0.5104. Based on UPGMA, the 200 proso millet accessions were divided into six groups. Based on Structure (=7), proso millet were mainly from four gene pools in Northeast, Northern, Loess Plateau, and North China. Based on principal component analysis, the accessions were classified into six clusters, consistent with the geographic region. This study constructed a set of 52 tetra-, penta- and hexa-nucleotide repeat EST-SSR markers, with the efficiency of 26%. Based on the above markers, genetic diversity of proso millet for landrace core germplasm was abundant.

    proso millet; core germplasm; high motif EST-SSR; genetic diversity; population structure

    2023-02-13;

    2023-02-28.

    10.3724/SP.J.1006.2023.24201

    通信作者(Corresponding authors):王瑞云, E-mail: wry925@126.com; 劉敏軒, E-mail: liuminxuan@caas.cn

    E-mail: 457836483@qq.com

    2022-08-31;

    本研究由財(cái)政部和農(nóng)業(yè)農(nóng)村部國(guó)家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng)(CARS-06-14.5-A2), 農(nóng)作物種質(zhì)資源保護(hù)與利用專項(xiàng)(2022NWB036-06), 中國(guó)農(nóng)業(yè)科學(xué)院科技創(chuàng)新工程, 吉林省自然科學(xué)基金(20200201153JC)和山西省現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng)資金(2023CYJSTX03-12)資助。

    This study was supported by the China Agriculture Research System of MOF and MARA (CARS-06-14.5-A2), the Protection and Utilization for Crop Germplasm Resources (2022NWB036-06), the Agricultural Science and Technology Innovation Program of the Chinese Academy of Agricultural Sciences, the Natural Science Foundation of Jilin Province (20200201153JC), and the Earmarked Fund for Modern Agro-industry Technology Research System (2023CYJSTX03-12).

    URL: https://kns.cnki.net/kcms/detail//11.1809.S.20230227.1702.014.html

    This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

    猜你喜歡
    糜子基元種質(zhì)
    關(guān)注基元反應(yīng)的考查
    華南地區(qū)最大農(nóng)作物種質(zhì)資源保護(hù)庫(kù)建成
    糜子品種理化特性與體外抗氧化性研究
    山西構(gòu)建糜子DNA分子身份證
    基于網(wǎng)絡(luò)調(diào)研的我國(guó)糜子消費(fèi)現(xiàn)狀分析
    人體細(xì)胞內(nèi)存在全新DNA結(jié)構(gòu)
    亞麻抗白粉病種質(zhì)資源的鑒定與篩選
    貴州玉米種質(zhì)資源遺傳多樣性及核心種質(zhì)庫(kù)構(gòu)建
    紅錐種質(zhì)早期生長(zhǎng)表現(xiàn)
    小小糜子富了一鎮(zhèn)百姓
    国产精品人妻久久久久久| 亚洲av成人精品一区久久| 插阴视频在线观看视频| 国产精品一及| 日韩人妻高清精品专区| 看免费成人av毛片| 日本欧美视频一区| 亚洲成人一二三区av| 午夜老司机福利剧场| 多毛熟女@视频| 久久精品国产自在天天线| 我的老师免费观看完整版| 日日摸夜夜添夜夜添av毛片| 观看免费一级毛片| 免费在线观看成人毛片| 国产精品久久久久久av不卡| 有码 亚洲区| 亚洲欧美成人综合另类久久久| 天堂中文最新版在线下载| 中文资源天堂在线| 乱码一卡2卡4卡精品| 性色av一级| 日本一二三区视频观看| 免费大片18禁| 亚洲国产欧美在线一区| 七月丁香在线播放| 精品人妻视频免费看| 国产在线免费精品| 亚洲三级黄色毛片| 国产美女午夜福利| 综合色丁香网| 精品一区二区免费观看| 久久人妻熟女aⅴ| 精品久久久噜噜| 毛片女人毛片| xxx大片免费视频| 男人舔奶头视频| 亚洲国产日韩一区二区| 夜夜爽夜夜爽视频| 少妇熟女欧美另类| 天美传媒精品一区二区| 一区在线观看完整版| 人人妻人人看人人澡| 少妇熟女欧美另类| 欧美xxxx黑人xx丫x性爽| 搡女人真爽免费视频火全软件| 色5月婷婷丁香| 啦啦啦啦在线视频资源| 少妇裸体淫交视频免费看高清| 黄片无遮挡物在线观看| 亚洲欧美日韩东京热| 精品人妻一区二区三区麻豆| 老师上课跳d突然被开到最大视频| 最近2019中文字幕mv第一页| 18禁在线无遮挡免费观看视频| 国产视频首页在线观看| 十分钟在线观看高清视频www | 最近的中文字幕免费完整| 中文字幕精品免费在线观看视频 | 视频中文字幕在线观看| 最近最新中文字幕免费大全7| av免费观看日本| 熟女电影av网| 日本黄色日本黄色录像| 大又大粗又爽又黄少妇毛片口| 特大巨黑吊av在线直播| 熟女人妻精品中文字幕| 国产成人a∨麻豆精品| 日韩强制内射视频| 自拍欧美九色日韩亚洲蝌蚪91 | 国产精品偷伦视频观看了| 男女啪啪激烈高潮av片| 精品熟女少妇av免费看| 三级经典国产精品| 亚洲久久久国产精品| 91精品一卡2卡3卡4卡| 免费少妇av软件| 国产 一区 欧美 日韩| 成年av动漫网址| 免费看av在线观看网站| 麻豆乱淫一区二区| 亚洲欧美中文字幕日韩二区| 我要看黄色一级片免费的| videossex国产| 久久97久久精品| 久久人人爽av亚洲精品天堂 | 男女无遮挡免费网站观看| 一个人免费看片子| 内地一区二区视频在线| 晚上一个人看的免费电影| 美女高潮的动态| 秋霞在线观看毛片| 五月天丁香电影| 我要看黄色一级片免费的| 亚洲va在线va天堂va国产| 国产成人a∨麻豆精品| 女人十人毛片免费观看3o分钟| 欧美一区二区亚洲| 国产视频首页在线观看| 国国产精品蜜臀av免费| 国产精品人妻久久久影院| 美女国产视频在线观看| 99久久精品热视频| 在线观看一区二区三区| 亚洲天堂av无毛| 亚洲欧美精品专区久久| 这个男人来自地球电影免费观看 | 九色成人免费人妻av| 国产亚洲91精品色在线| 免费不卡的大黄色大毛片视频在线观看| 亚洲美女视频黄频| 搡女人真爽免费视频火全软件| 亚洲欧洲国产日韩| 精品国产乱码久久久久久小说| 国产 一区精品| 亚洲美女黄色视频免费看| av在线观看视频网站免费| 在线观看免费日韩欧美大片 | 蜜桃亚洲精品一区二区三区| 女的被弄到高潮叫床怎么办| 男人狂女人下面高潮的视频| 99热6这里只有精品| 少妇 在线观看| 日韩一区二区三区影片| 亚洲精品日本国产第一区| 国产精品久久久久久精品电影小说 | 一二三四中文在线观看免费高清| av在线播放精品| 久久青草综合色| 久久 成人 亚洲| 国产成人一区二区在线| 亚洲国产精品国产精品| 丝袜脚勾引网站| 欧美精品亚洲一区二区| 狂野欧美激情性bbbbbb| 青春草国产在线视频| 超碰97精品在线观看| 国产精品久久久久久精品古装| 国产精品久久久久久精品古装| 晚上一个人看的免费电影| 亚洲精品乱码久久久v下载方式| 亚洲av在线观看美女高潮| 观看免费一级毛片| 精品一区二区三卡| 国产欧美日韩精品一区二区| 精品一品国产午夜福利视频| 人人妻人人看人人澡| 麻豆成人av视频| 女人久久www免费人成看片| 黑丝袜美女国产一区| 天堂中文最新版在线下载| av网站免费在线观看视频| 成人毛片60女人毛片免费| 亚洲av中文av极速乱| 只有这里有精品99| 最黄视频免费看| 亚洲怡红院男人天堂| 黑人高潮一二区| 国产精品免费大片| 舔av片在线| 高清不卡的av网站| 性色avwww在线观看| 高清欧美精品videossex| 国产美女午夜福利| 热re99久久精品国产66热6| 如何舔出高潮| 国产在线一区二区三区精| 欧美一级a爱片免费观看看| 亚洲丝袜综合中文字幕| 亚洲无线观看免费| 国产精品一区二区三区四区免费观看| 最近中文字幕2019免费版| 欧美xxⅹ黑人| 久久精品人妻少妇| 亚洲欧美精品专区久久| 国产精品av视频在线免费观看| 国产一区亚洲一区在线观看| 亚洲国产欧美人成| 久久久久久久精品精品| 丰满少妇做爰视频| 一本久久精品| 国产人妻一区二区三区在| 亚洲成人一二三区av| 尤物成人国产欧美一区二区三区| 精品视频人人做人人爽| 人人妻人人爽人人添夜夜欢视频 | 中文字幕制服av| 国产日韩欧美亚洲二区| 久久婷婷青草| 国产一区二区三区综合在线观看 | 菩萨蛮人人尽说江南好唐韦庄| 成人漫画全彩无遮挡| 男女国产视频网站| 在线观看国产h片| 国产黄片视频在线免费观看| 卡戴珊不雅视频在线播放| 久久韩国三级中文字幕| 国产在视频线精品| 插阴视频在线观看视频| 亚洲va在线va天堂va国产| 亚洲精品国产av成人精品| 亚洲国产最新在线播放| 18禁裸乳无遮挡免费网站照片| 欧美xxxx黑人xx丫x性爽| 国产老妇伦熟女老妇高清| 三级国产精品片| 国产无遮挡羞羞视频在线观看| 麻豆乱淫一区二区| 下体分泌物呈黄色| 99热6这里只有精品| 最近的中文字幕免费完整| 成年人午夜在线观看视频| 日本欧美国产在线视频| 一本久久精品| 一本久久精品| 一本一本综合久久| 男人和女人高潮做爰伦理| 热99国产精品久久久久久7| 免费观看的影片在线观看| 一区二区三区乱码不卡18| 中文乱码字字幕精品一区二区三区| 久久人人爽人人片av| 国产乱人偷精品视频| 亚洲欧美成人精品一区二区| 亚洲色图av天堂| 在线看a的网站| 黄片wwwwww| 国产国拍精品亚洲av在线观看| 在线观看美女被高潮喷水网站| 色视频在线一区二区三区| 91久久精品国产一区二区三区| 国产一级毛片在线| 久久久久久久久久久丰满| kizo精华| 熟女电影av网| 欧美日韩视频精品一区| 日韩三级伦理在线观看| 嫩草影院新地址| 国产精品一区www在线观看| 国产久久久一区二区三区| 蜜桃在线观看..| 全区人妻精品视频| 欧美日韩在线观看h| 日韩三级伦理在线观看| 久久97久久精品| 久久精品夜色国产| 中文在线观看免费www的网站| 久久亚洲国产成人精品v| 亚洲第一区二区三区不卡| 日本av免费视频播放| 免费久久久久久久精品成人欧美视频 | 久久国产精品大桥未久av | 只有这里有精品99| 精品午夜福利在线看| 亚洲,一卡二卡三卡| 最近的中文字幕免费完整| 丝瓜视频免费看黄片| av在线老鸭窝| 少妇高潮的动态图| 国产亚洲5aaaaa淫片| 18禁裸乳无遮挡动漫免费视频| 久久97久久精品| 国产精品国产三级国产专区5o| 精品亚洲成a人片在线观看 | 在线观看国产h片| 久久久色成人| 久久久久久久久久久丰满| 在线精品无人区一区二区三 | 色婷婷久久久亚洲欧美| 欧美成人午夜免费资源| 人人妻人人看人人澡| 1000部很黄的大片| 亚洲综合色惰| 高清不卡的av网站| 成年免费大片在线观看| 成人一区二区视频在线观看| 国产男女超爽视频在线观看| 精华霜和精华液先用哪个| 国产视频首页在线观看| 国产老妇伦熟女老妇高清| 新久久久久国产一级毛片| 人妻制服诱惑在线中文字幕| 国产极品天堂在线| 777米奇影视久久| 午夜福利影视在线免费观看| 99热这里只有精品一区| 久久亚洲国产成人精品v| 久久久久久九九精品二区国产| 日韩大片免费观看网站| 美女cb高潮喷水在线观看| 亚洲av在线观看美女高潮| 99精国产麻豆久久婷婷| 日韩伦理黄色片| 高清不卡的av网站| 国产91av在线免费观看| 久久鲁丝午夜福利片| 亚洲精品亚洲一区二区| www.色视频.com| 大香蕉97超碰在线| 中文字幕久久专区| 少妇被粗大猛烈的视频| 午夜福利视频精品| 涩涩av久久男人的天堂| 成人毛片a级毛片在线播放| 2018国产大陆天天弄谢| 国产 精品1| 亚洲欧美成人精品一区二区| 久久精品国产自在天天线| 大片免费播放器 马上看| 中文字幕亚洲精品专区| 国产男女超爽视频在线观看| 97超视频在线观看视频| 久久国产亚洲av麻豆专区| 国产日韩欧美亚洲二区| 亚洲激情五月婷婷啪啪| 舔av片在线| 51国产日韩欧美| 欧美日本视频| 国内揄拍国产精品人妻在线| 国产精品久久久久久精品电影小说 | 亚洲欧美日韩东京热| 各种免费的搞黄视频| 中文乱码字字幕精品一区二区三区| 一个人免费看片子| 亚洲成人av在线免费| 又大又黄又爽视频免费| 欧美一区二区亚洲| 欧美丝袜亚洲另类| av视频免费观看在线观看| 天天躁夜夜躁狠狠久久av| 日本猛色少妇xxxxx猛交久久| 老司机影院成人| 国产一区有黄有色的免费视频| 欧美日韩综合久久久久久| 亚洲国产av新网站| 一二三四中文在线观看免费高清| 久久久久久久大尺度免费视频| 免费大片黄手机在线观看| 最黄视频免费看| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品国产av蜜桃| 性高湖久久久久久久久免费观看| 精华霜和精华液先用哪个| 久久97久久精品| 免费黄网站久久成人精品| 大香蕉久久网| 久久精品国产自在天天线| 国产成人免费无遮挡视频| 欧美高清性xxxxhd video| 亚洲中文av在线| 国内揄拍国产精品人妻在线| 日日摸夜夜添夜夜爱| 高清毛片免费看| 成人亚洲欧美一区二区av| 日日撸夜夜添| 性色avwww在线观看| tube8黄色片| 高清黄色对白视频在线免费看 | 七月丁香在线播放| 婷婷色综合www| 成人综合一区亚洲| 久久热精品热| 青春草亚洲视频在线观看| 婷婷色av中文字幕| 黄色一级大片看看| 久久精品国产自在天天线| 毛片一级片免费看久久久久| 丝瓜视频免费看黄片| 九九在线视频观看精品| 最近最新中文字幕免费大全7| 国内精品宾馆在线| 亚洲丝袜综合中文字幕| 午夜日本视频在线| 成人影院久久| 蜜桃亚洲精品一区二区三区| 国产白丝娇喘喷水9色精品| 高清日韩中文字幕在线| 成年人午夜在线观看视频| 国产欧美日韩精品一区二区| 汤姆久久久久久久影院中文字幕| 人人妻人人添人人爽欧美一区卜 | 国产成人精品福利久久| 久久久久久伊人网av| 老司机影院成人| 国产午夜精品一二区理论片| 日本av免费视频播放| 久久韩国三级中文字幕| 免费黄频网站在线观看国产| 欧美精品人与动牲交sv欧美| 一级二级三级毛片免费看| 18禁动态无遮挡网站| 国产av精品麻豆| 免费久久久久久久精品成人欧美视频 | 看非洲黑人一级黄片| 国产精品久久久久久av不卡| 99热这里只有是精品50| 亚洲一区二区三区欧美精品| 亚洲精品久久午夜乱码| av黄色大香蕉| 国产永久视频网站| 午夜免费鲁丝| 中文字幕人妻熟人妻熟丝袜美| 国产高潮美女av| 天堂8中文在线网| 亚洲国产精品国产精品| 交换朋友夫妻互换小说| 免费大片18禁| 国产乱人偷精品视频| 欧美精品一区二区大全| 中国国产av一级| 丝瓜视频免费看黄片| 美女脱内裤让男人舔精品视频| 一区在线观看完整版| 性色avwww在线观看| 久久久精品免费免费高清| 欧美日韩在线观看h| 能在线免费看毛片的网站| a级一级毛片免费在线观看| 一区在线观看完整版| 亚洲欧美清纯卡通| 国产又色又爽无遮挡免| 视频区图区小说| 成年人午夜在线观看视频| 欧美人与善性xxx| 我的女老师完整版在线观看| 日韩 亚洲 欧美在线| 少妇高潮的动态图| 中文精品一卡2卡3卡4更新| 热re99久久精品国产66热6| 欧美日韩综合久久久久久| 激情五月婷婷亚洲| 欧美成人午夜免费资源| 久久久久久九九精品二区国产| 啦啦啦中文免费视频观看日本| 精品亚洲成a人片在线观看 | 美女cb高潮喷水在线观看| 91精品一卡2卡3卡4卡| 少妇人妻久久综合中文| 妹子高潮喷水视频| 22中文网久久字幕| 在线观看免费高清a一片| 欧美一级a爱片免费观看看| 亚洲精品日韩av片在线观看| 视频中文字幕在线观看| 乱码一卡2卡4卡精品| 亚洲av中文av极速乱| 毛片一级片免费看久久久久| 国产免费又黄又爽又色| 日韩av不卡免费在线播放| 黄色怎么调成土黄色| 国产日韩欧美亚洲二区| 日本av免费视频播放| 久久久色成人| 国产精品一区二区在线观看99| 久久综合国产亚洲精品| 精品一区二区三卡| 亚洲va在线va天堂va国产| 国产在线免费精品| 天堂8中文在线网| 亚洲av中文字字幕乱码综合| 丰满人妻一区二区三区视频av| 欧美 日韩 精品 国产| 这个男人来自地球电影免费观看 | 精品人妻一区二区三区麻豆| 97精品久久久久久久久久精品| 日韩在线高清观看一区二区三区| www.色视频.com| 国产视频内射| 亚洲电影在线观看av| 最近中文字幕高清免费大全6| 99热这里只有是精品50| av免费观看日本| 国产精品一区www在线观看| 成人亚洲精品一区在线观看 | 99九九线精品视频在线观看视频| 极品少妇高潮喷水抽搐| 亚洲欧美精品专区久久| 日韩亚洲欧美综合| 亚洲中文av在线| 亚洲精品国产色婷婷电影| 丝袜喷水一区| 久久精品国产自在天天线| 久久精品夜色国产| 高清午夜精品一区二区三区| 日韩成人av中文字幕在线观看| 欧美97在线视频| 亚洲美女视频黄频| 国产乱来视频区| av在线蜜桃| 99久久综合免费| 伦精品一区二区三区| 精品久久久久久久久av| 精品人妻一区二区三区麻豆| 成人黄色视频免费在线看| 黄色视频在线播放观看不卡| 久久久久久久久久人人人人人人| 久久久久人妻精品一区果冻| 黄色欧美视频在线观看| 亚洲怡红院男人天堂| 亚洲av综合色区一区| 国产永久视频网站| 2021少妇久久久久久久久久久| 大香蕉久久网| 欧美一区二区亚洲| 99久久精品一区二区三区| 国产国拍精品亚洲av在线观看| 免费黄频网站在线观看国产| 一区在线观看完整版| 国产亚洲一区二区精品| 国产永久视频网站| 亚洲欧美中文字幕日韩二区| 插逼视频在线观看| 亚洲国产最新在线播放| 免费看光身美女| 久久ye,这里只有精品| 青春草视频在线免费观看| 婷婷色综合大香蕉| av不卡在线播放| 国产一区二区三区综合在线观看 | 男人舔奶头视频| 性色avwww在线观看| 高清视频免费观看一区二区| 各种免费的搞黄视频| 欧美日韩亚洲高清精品| 22中文网久久字幕| 欧美xxxx黑人xx丫x性爽| 国产精品不卡视频一区二区| 最近最新中文字幕免费大全7| 欧美3d第一页| 永久免费av网站大全| 六月丁香七月| 国产成人一区二区在线| 欧美精品国产亚洲| 国产乱人偷精品视频| 久久精品国产自在天天线| 日本黄色片子视频| 国产极品天堂在线| 欧美xxxx性猛交bbbb| 成人毛片a级毛片在线播放| 国内精品宾馆在线| 亚洲国产高清在线一区二区三| 国产免费一区二区三区四区乱码| 黄色配什么色好看| 亚洲精品久久午夜乱码| 亚洲,欧美,日韩| 99九九线精品视频在线观看视频| 欧美高清成人免费视频www| 夫妻午夜视频| 久久这里有精品视频免费| 少妇的逼好多水| 美女中出高潮动态图| 啦啦啦在线观看免费高清www| 国产一区有黄有色的免费视频| 最近最新中文字幕免费大全7| 国产日韩欧美亚洲二区| 色婷婷av一区二区三区视频| 亚洲精品一区蜜桃| 久久99热这里只有精品18| av线在线观看网站| 亚洲国产av新网站| 亚洲图色成人| 一区二区av电影网| 香蕉精品网在线| 国产一区二区在线观看日韩| 欧美极品一区二区三区四区| 亚洲精品日韩av片在线观看| 嘟嘟电影网在线观看| 亚洲天堂av无毛| 晚上一个人看的免费电影| 欧美日韩在线观看h| 一级av片app| 黄片wwwwww| 欧美高清性xxxxhd video| 一区二区av电影网| 国产乱人偷精品视频| av免费在线看不卡| 亚洲人成网站高清观看| 亚洲av福利一区| 国产欧美日韩精品一区二区| 国产精品三级大全| av视频免费观看在线观看| 国产视频首页在线观看| 中文资源天堂在线| 天美传媒精品一区二区| 国产欧美亚洲国产| 美女福利国产在线 | 国产 一区精品| 嘟嘟电影网在线观看| 久久精品久久久久久噜噜老黄| 亚洲美女搞黄在线观看| 菩萨蛮人人尽说江南好唐韦庄| 国产日韩欧美亚洲二区| 日韩三级伦理在线观看| 美女中出高潮动态图| 男女边吃奶边做爰视频| 91精品一卡2卡3卡4卡| 91在线精品国自产拍蜜月| 女人久久www免费人成看片| 欧美成人一区二区免费高清观看| 成人亚洲欧美一区二区av| 卡戴珊不雅视频在线播放| 伦理电影大哥的女人| 亚洲成人手机| av福利片在线观看| 99久久精品国产国产毛片| 人妻夜夜爽99麻豆av| 水蜜桃什么品种好| 亚洲av中文av极速乱| 亚洲成人一二三区av| 搡女人真爽免费视频火全软件| 久久久a久久爽久久v久久| 中文精品一卡2卡3卡4更新| 国产在线男女| 精品久久久久久久久av| 午夜免费男女啪啪视频观看|