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

    黃淮麥區(qū)小麥粒重基因等位變異的分子鑒定及育種應(yīng)用

    2021-08-31 06:46:16張福彥程仲杰陳曉杰王嘉歡范家霖張建偉楊保安
    作物學(xué)報 2021年11期
    關(guān)鍵詞:等位粒重克隆

    張福彥 程仲杰 陳曉杰 王嘉歡 陳 鋒 范家霖 張建偉,* 楊保安,*

    黃淮麥區(qū)小麥粒重基因等位變異的分子鑒定及育種應(yīng)用

    張福彥1程仲杰1陳曉杰1王嘉歡1陳 鋒2范家霖1張建偉1,*楊保安1,*

    1河南省科學(xué)院同位素研究所有限責(zé)任公司 / 河南省核農(nóng)學(xué)重點實驗室, 河南鄭州 450015;2河南農(nóng)業(yè)大學(xué)農(nóng)學(xué)院 / 省部共建小麥玉米作物學(xué)國家重點實驗室, 河南鄭州 450002

    采用特異性引物PCR擴增方法對183份黃淮海麥區(qū)的小麥品種(系)的粒重基因、和的等位變異進行分子鑒定, 并結(jié)合2016—2017和2017—2018年度的千粒重表型數(shù)據(jù), 分析不同等位變異類型對小麥粒重的影響, 從而找出優(yōu)勢基因型組合。結(jié)果表明, 不同年份間參試品種(系)的千粒重差異達到極顯著水平(<0.01);位點上發(fā)現(xiàn)和兩種等位變異, 其分布頻率分別為66.7%和33.3%;位點上等位變異分布頻率較高, 為94.5%, 而等位變異分布頻率極低, 僅為5.5%;位點上發(fā)現(xiàn)和兩種等位變異, 其分布頻率分別為79.8%和20.2%。不同等位變異組合的品種千粒重存在顯著差異(<0.05), 其中具有3個高千粒重等位變異組合品種的平均千粒重最高, 與具有品種的平均千粒重差異不顯著, 但是顯著高于其他組合(<0.05)。基因型組合小麥品種(系)的平均千粒重最低。、和位點上的不同等位變異均會導(dǎo)致小麥千粒重的顯著變化, 其中和位點上的等位變異對小麥粒重的影響更為重要。在參試材料中沒有發(fā)現(xiàn)具有3個低千粒重等位變異組合的品種, 在7種不同等位變異組合中, 具有3個高千粒重等位變異組合品種的平均千粒重最高, 是優(yōu)勢基因型組合。

    小麥; 粒重基因; 功能標記; 等位變異; 千粒重

    小麥是全球最為重要的糧食作物之一, 世界上約40%的人口以小麥為主要食糧, 我國小麥的豐歉直接影響到國家乃至世界糧食安全, 因此高產(chǎn)、穩(wěn)產(chǎn)、抗病、優(yōu)質(zhì)、廣適是現(xiàn)階段我國小麥育種中最主要的育種目標之一[1]。小麥產(chǎn)量由單位面積穗數(shù)、穗粒數(shù)以及粒重三因素構(gòu)成, 粒重是三要素中遺傳相對穩(wěn)定的性狀, 主要受遺傳因子影響且主要受加性效應(yīng)的基因控制[2]。粒重是影響小麥產(chǎn)量最為重要因素之一, 粒重在我國不同的小麥生態(tài)區(qū)均被不同程度的正向選擇和提高, 是現(xiàn)代小麥育種中遺傳改良最為顯著的產(chǎn)量性狀, 對我國小麥單產(chǎn)水平的提高做出了較大貢獻[3-4], 但目前關(guān)于小麥粒重基因的遺傳基礎(chǔ)解析和粒重相關(guān)的功能標記在育種中的應(yīng)用仍不夠深入。

    近年來, 隨著小麥全基因組序列信息的公布及現(xiàn)代分子克隆技術(shù)的快速發(fā)展, 關(guān)于小麥粒重基因的克隆和功能標記開發(fā)也取得較快發(fā)展。Ma等[5]在普通小麥的2AL染色體上克隆得到細胞壁轉(zhuǎn)化酶基因, 發(fā)現(xiàn)該基因的2種等位變異, 分別命名為和, 并根據(jù)二者序列差異開發(fā)了一對互補的顯性功能標記CWI21和CWI22, 進一步研究發(fā)現(xiàn)變異類型與高千粒重的品種相關(guān), 而類型則與低千粒重的品種相關(guān)。之后, Jiang等[6]在普通小麥中克隆出基因, 其中和基因也都發(fā)現(xiàn)了相應(yīng)的單倍型, 進一步研究發(fā)現(xiàn)具有和單倍型的品種千粒重顯著高于具有和單倍型的品種。Zhang等[7]利用同源克隆從普通小麥7DS染色體上克隆得到基因, 發(fā)現(xiàn)其含有和兩種等位變異, 根據(jù)基因序列多態(tài)性開發(fā)的共顯性標記GS7D, 能夠準確檢測和類型, 發(fā)現(xiàn)擁有不同等位變異類型品種的千粒重間存在顯著性差異。Wang等[8]利用同源克隆技術(shù)從普通小麥中克隆了一個控制籽粒大小的基因, 該基因被定位在3AS染色體上, 并發(fā)現(xiàn)其至少存在2種等位變異, 分別命名為和。之后, 又在該基因啟動子上游順式作用元件Sp1內(nèi)發(fā)現(xiàn)單堿基G的插入, 這種變化對小麥的粒寬和粒重具有顯著影響[9]。Hanif等[10]和Hu等[11]利用水稻中的基因在普通小麥進行同源克隆, Hanif等[10]克隆得到位于普通小麥3AL染色體上的基因, 并發(fā)現(xiàn)該基因在不同品種中存在單核苷酸多態(tài)性, 分別命名為和, 而Hu等[11]克隆得到基因, 并將其定位在普通小麥4AL染色體上, 且在位點發(fā)現(xiàn)3種等位變異類型, 分別命名為:、和, 并開發(fā)了共顯性標記TG23, 進一步研究發(fā)現(xiàn)和變異類型的品種在粒徑和粒重方面顯著高于類型, 但和變異類型在當(dāng)前品種中分布概率較低。最近, Yang等[12]同源克隆發(fā)現(xiàn)了控制小麥粒長和粒重的基因, 被定位在5AL染色體上, 根據(jù)多單核苷酸多態(tài)性發(fā)現(xiàn)存在2種等位變異, 分別命名為和, 后者是小麥粒長和粒重的優(yōu)勢等位變異, 但在我國小麥中分布頻率較低, 認為該類型在我國小麥育種利用上具有較大潛力。Ma等[13]利用同源克隆的方法獲得小麥基因, 并利用瞬時沉默和過表達技術(shù)揭示表達水平與小麥粒重呈正相關(guān)。之后, 司文潔等[14]又根據(jù)不同小麥品種啟動子區(qū)序列內(nèi)SNP位點差異開發(fā)啟動子區(qū)功能標記CAPS-5Ap, 并對其進行了驗證。Yan等[15]利用同源克隆的方法獲得控制小麥粒重的基因, 命名為, 位于7BS上。不同籽粒大小的小麥品種基因多態(tài)性分析, 發(fā)現(xiàn)基因存在2種等位變異, 命名為和, 并開發(fā)了功能標記TaGW8-7B, 并發(fā)現(xiàn)不同等位變異對小麥粒重有顯著影響。高產(chǎn)和穩(wěn)產(chǎn)是我國乃至世界小麥育種家永恒的育種目標, 粒重的增加為推動我國乃至世界小麥單產(chǎn)水平的整體提升做出重要貢獻, 而關(guān)于小麥粒重遺傳基礎(chǔ)解析已成為當(dāng)前研究熱點之一。

    截至目前, 盡管已有關(guān)于小麥粒重基因等位變異檢測的相關(guān)報道[16-17], 但同時對多個粒重基因進行分子鑒定及粒重功能標記育種應(yīng)用的研究少之甚少。本研究以183份黃淮海麥區(qū)的小麥品種(系)為材料, 利用功能標記CWI21、CWI22、TaGW8-7B以及GS7D分別對2AL、7BS和7DS染色體上的粒重基因等位變異進行檢測, 同時結(jié)合2016—2017和2017—2018年度的千粒重表型數(shù)據(jù), 分析不同等位變異類型對小麥粒重的影響, 并找出優(yōu)勢的粒重基因型組合, 旨在為培育出高產(chǎn)、穩(wěn)產(chǎn)的小麥新品種提供一些優(yōu)異親本資源, 也為粒重基因功能標記輔助育種提供參考依據(jù)。

    1 材料與方法

    1.1 試驗材料

    豫豐11、鄭品麥22、百農(nóng)207、濟麥22、石4185、西農(nóng)511、西農(nóng)979、周麥27、周麥32、周麥36、豫同194、豫農(nóng)202、魯原502、良星99、豫同843、豐德存麥1號、豐德存麥5號等183份小麥品種(系)來自河南省科學(xué)院同位素研究所小麥育種室收集的種質(zhì)資源圃, 這些材料主要涵蓋了黃淮海麥區(qū)近20年來生產(chǎn)上種植的小麥品種(系)。

    1.2 表型測定

    2016—2017和2017—2018年度在河南省科學(xué)院新鄭試驗基地種植, 每個品種(系)種植2行, 行長2.5 m, 行距0.25 m, 每行播50粒種子, 待出苗后每行定苗至30株左右。田間管理同大田生產(chǎn), 正常成熟后人工收獲, 及時晾曬, 儲藏備用。

    千粒重測定, 采用德國福弗(pfeuffer)公司生產(chǎn)的康達得型(contador)全自動數(shù)粒儀進行種子計數(shù), 每份材料均勻取1000粒, 重復(fù)3次, 同時采用常熟市天量儀器有限責(zé)任公司生產(chǎn)的LT1002E型電子天平進行稱重。

    1.3 基因型檢測

    每個材料選取3粒有代表性的種子, 粉碎后分別放入3個2 mL離心管中, 根據(jù)Chen等[18]提取籽粒DNA的方法, 采用SLS (十二酰肌氨酸鈉)快速提取小麥籽粒DNA法提取參試材料DNA, 并根據(jù)每個品種(系) 3種功能標記檢測判斷該品種(系)的基因型。

    利用已開發(fā)功能標記CWI21、CWI22、TaGW8- 7B和GS7D對所有參試材料中粒重基因、和的等位變異進行分子檢測(表1)。PCR反應(yīng)體系為20 μL, 主要包括buffer緩沖液2.5 μL,DNA 聚合酶(北京天根生化科技有限公司) (2.5 U μL–1) 0.4 μL, dNTP (200 μmol L–1) 1.0 μL, 上下游引物(10 pmol L–1)各0.5 μL, 模板DNA 2.0 μL, ddH2O 13.1 μL。在Tprofessional standard型PCR儀(Biometra)中進行擴增, PCR反應(yīng)程序為94℃預(yù)變性5 min, 94℃變性30 s, 退火溫度具體參照(表1)進行, 72℃延伸30 s, 30個循環(huán); 72℃延伸10 min, 反應(yīng)程序結(jié)束后取5 μL PCR擴增產(chǎn)物用1.5%瓊脂糖凝膠電泳分離、利用GeneGreen Nucleic Acid DyeGeneGreen核酸染料進行染色, Fire Reader凝膠成像系統(tǒng)下進行掃描成像, 并保存至計算機以備統(tǒng)計分析。

    1.4 統(tǒng)計分析

    采用SPSS 18.0軟件進行數(shù)據(jù)分析,對具有不同粒重基因型的小麥品種(系)間千粒重表型進行差異顯著性分析和方差分析(ANOVA),并利用LSD法進行多重比較。

    2 結(jié)果與分析

    2.1 不同年份間粒重的差異分析

    從表2可知,不同年份間供試品種(系)的千粒重差異達到極顯著水平。兩個年度間供試品種(系)的千粒重平均值分別為44.11 g和46.29 g,變幅分別為34.50~53.35 g和34.20~58.45 g,其變異系數(shù)均較高,分別為8.48%和8.75%??梢姡S淮海地區(qū)小麥粒重的遺傳變異非常豐富,具有較大的遺傳改良潛力。

    2.2 粒重基因不同等位變異分子檢測

    利用顯性互補功能標記CWI21和CWI22檢測基因在供試材料中的不同等位變異, 發(fā)現(xiàn)CWI21標記在石4185、鄭品麥22、鄭育麥16、豫豐11、農(nóng)大2011、輪選166、豫農(nóng)186和淮麥19等61份材料中可擴增出404 bp的片段, 屬于類型(圖1-B), 而CWI22標記在商麥167、中育1220、渦麥66、泛育麥17、鄭麥103、鄭麥369和周麥36等122個材料中擴增出402 bp的片段, 屬于類型(圖1-A)。

    表2 不同年份中供試材料的千粒重方差分析

    **表示差異達到< 0.01顯著水平。**Significant at< 0.01.

    圖1 標記CWI21 (A)和CWI22 (B)在部分小麥品種中擴增結(jié)果

    M: DL2000; 1: 石4185; 2: 鄭品麥22; 3: 商麥167; 4: 鄭育麥16; 5: 豫豐11; 6: 中育1220; 7: 渦麥66; 8: 泛育麥17; 9: 豫農(nóng)186; 10: 賽德麥1號; 11: 鄭麥369; 12: 淮麥19。

    M: DL2000; 1: Shi 4185; 2: Zhengpinmai 22: 3: Shangmai 167; 4: Zhengyumai 16; 5: Yufeng 11; 6: Zhongyu 1220; 7: Guomai 66; 8: Fanyumai 17; 9: Yunong 186; 10: Saidemai 1; 11: Zhengmai 369; 12: Huaimai 19.

    利用共顯性功能標記TaGW8-7B對所有參試材料進行分子檢測(圖2), 發(fā)現(xiàn)在商麥167、鑫農(nóng)518、中新16、豫豐11、皖科06725、輪選166、泉麥29、農(nóng)大2011、鄭麥1860等173份材料中擴增出1097 bp的片段, 屬于類型, 而在溫麥6號、鄭麥518、豫農(nóng)202、石84-7111、淮麥19等10份材料中擴增出1373 bp的片段, 屬于類型。

    利用共顯性標記GS7D對183份小麥品種(系)進行檢測(圖3), 發(fā)現(xiàn)在魯原502、濮興0369、輪選987、周8425B、蘭考906、周麥27、西農(nóng)511、中育1220、豫豐11等146份材料中擴增出562 bp的片段, 屬于類型, 而在豫農(nóng)186、賽德麥1號、龍科1221、中育9302等37份材料中可擴增出522 bp的片段, 屬于類型。

    圖2 TaGW8-7B標記在部分小麥品種中擴增結(jié)果

    1: 鄭品麥22; 2: 商麥167; 3: 溫麥6號; 4: 鄭麥518; 5: 豫農(nóng)202; 6: 中新16; 7: 豫豐11; 8: 皖科06725; 9: 豫麥2號; 10: 冀麥38; 11: 石84-7111; 12: 豐德存麥5號; 13: 良星99; 14: 洛麥28; 15: 淮麥19; 16: 百農(nóng)3217; M: DL2000。

    1: Zhengpinmai 22; 2: Shangmai 167; 3: Wenmai 6; 4: Zhengmai 518; 5: Yunong 202; 6: Zhongxin 16; 7: Yufeng 11; 8: Wanke 06725; 9: Yumai 2; 10: Jimai 38; 11: Shi 84-7111; 12: Fengdecunmai 5; 13: Liangxing 99; 14: Luomai 28; 15: Huaimai 19; 16: Bainong 3217; M: DL2000.

    圖3 GS7D標記檢測部分小麥品種

    1: 豫農(nóng)186; 2: 賽德麥1號; 3: 魯原502; 4: 鄭麥366; 5: 中育9302; 6: 龍科1221; 7: 濮興0369; 8: 豫豐11; 9: 農(nóng)大2011; 10: 鄭品麥8號; 11: 高麥6號; 12: 豐德存麥5號; M: DL2000。

    1: Yunong 186; 2: Saidemai 1; 3: Luyuan 502; 4: Zhengmai 366; 5: Zhongyu 9302; 6: Longke 1221; 7: Puxing 0369; 8: Yufeng 11; 9: Nongda 2011; 10: Zhengpinmai 8; 11: Gaomai 6; 12: Fengdecunmai 5; M: DL2000.

    2.3 粒重基因不同等位變異的分布及其對小麥粒重的影響

    在所有參試材料中,、和位點具上均發(fā)現(xiàn)2種等位變異。從表3可知,位點上鑒定發(fā)現(xiàn)183份參試材料中等位變異的分布頻率為66.67%,為33.33%, 同時利用SPSS 18.0軟件對不同等位變異對應(yīng)的千粒重進行方差分析, 發(fā)現(xiàn)具有等位變異類型品種的千粒重顯著高于具有等位變異類型品種的千粒重(<0.05),為優(yōu)異等位基因。位點上鑒定所有參試材料中等位變異分布頻率較高,為94.54%, 而等位變異分布頻率極低, 僅為5.46%, 方差分析發(fā)現(xiàn)具有等位變異類型品種的千粒重與具有等位變異類型品種的千粒重, 二者差異達到極顯著水平(<0.01),為優(yōu)異等位基因。位點上鑒定所有參試材料中等位變異分布頻率為79.78%,為20.22%, 方差分析發(fā)現(xiàn)具有等位變異類型品種的千粒重與具有等位變異類型品種的千粒重, 二者差異也達到極顯著水平(<0.01),為優(yōu)異等位基因。說明、和位點上不同等位變異均會導(dǎo)致小麥千粒重的顯著變化, 而和位點上的等位變異對粒重具有更為重要的影響。

    表3 粒重基因的不同等位變異對小麥千粒重影響

    *和**分別表示差異達到< 0.05和< 0.01顯著水平。

    *and**indicate significant difference at< 0.05 and< 0.01, respectively.

    2.4 不同基因型組合與小麥粒重的關(guān)系

    利用功能標記CWI21、CWI22、TaGW8-7B和GS7D對所有參試材料進行檢測。在、和位點上共發(fā)現(xiàn)、、、、、和共7種不同的等位變異組合, 分別占參試材料的52.5%、23.5%、7.1%、11.5%、2.7%、1.6%和1.1%。從表4可知, 具有不同等位變異組合品種的千粒重間存在顯著差異。具有3個高千粒重等位變異組合品種的平均千粒重最高, 為44.69 g, 與具有品種的平均千粒重差異不顯著, 但是顯著高于其他組合; 在參試材料中沒有發(fā)現(xiàn)具有3個低千粒重等位變異組合的品種;具有、、和等位變異組合品種(系)間的平均千粒重差異也不顯著; 而基因型組合小麥品種(系)的平均千粒重最低, 僅為41.39 g。可見,、和位點上不同等位變異組合會引起小麥粒重的顯著變化。

    表4 TaCwi-A1、TaGS-D1和TaGw8-B1位點上不同等位變異組合對小麥千粒重的影響

    ?H: 高千粒重; M: 中等千粒重; ML: 中低千粒重; L: 低千粒重。?不同小寫字母表示在< 0.05水平顯著差異。

    ?H: high thousand-grain weight; M: medium thousand-grain weight; ML: medium low thousand-grain weight; L:low thousand-grain weight.?Different lowercase letters indicate significant difference at< 0.05.

    3 討論

    小麥產(chǎn)量由單位面積穗數(shù)、穗粒數(shù)以及千粒重三因素構(gòu)成, 其中穗粒數(shù)的增加是建立在單位面積穗數(shù)減少的基礎(chǔ)上, 而千粒重的增加則是相對獨立的, 且受遺傳特性影響最大, 其遺傳力高達89%[19]。Wu等[20]研究發(fā)現(xiàn)1945年至2010年以來我國不同的小麥主產(chǎn)區(qū)小麥產(chǎn)量平均每年增產(chǎn)幅度都在0.33%~1.42%之間, 而粒重和單穗重的顯著增加對小麥產(chǎn)量的持續(xù)增長起到了舉足輕重的作用。Rasheed等[21]研究表明, 在普通小麥中目前已知的有、、、、和等控制粒長、粒寬、粒重的基因均具有正向遺傳效應(yīng), 且對小麥粒重高低具有極其重要的影響。同時, 除了基因型影響外, 脫落酸、赤霉素等內(nèi)源激素以及光照、溫度、水分、養(yǎng)分等外界環(huán)境等多種因素對小麥粒重也產(chǎn)生重要影響。本文僅對粒重控制基因?qū)π←溋V氐挠绊懸约袄昧V毓δ軜擞涍M行育種應(yīng)用進行了研究和討論。

    小麥粒重是一個受多基因控制的復(fù)雜數(shù)量性狀,與其相關(guān)的大多數(shù)QTL位點的表型貢獻率較小, 重復(fù)性較差, 且較為分散。Gupta等[22]對小麥粒長、粒寬和千粒重的QTL遺傳定位的研究結(jié)果進行了系統(tǒng)總結(jié), 發(fā)現(xiàn)受親本遺傳背景差異、QTL作圖方法、標記類型等因素的影響, 在不同的作圖群體中小麥千粒重的QTL定位的結(jié)果存在較大的差異, 且遠遠不能滿足分子標記輔助選擇育種和粒重相關(guān)基因克隆的需要。隨著圖位克隆、同源克隆等現(xiàn)代生物學(xué)技術(shù)和標記開發(fā)技術(shù)的快速發(fā)展, 目前已先后克隆出、、、、、、、和等小麥粒重形成的相關(guān)基因, 并開發(fā)出GS7D、TG23、TaGW8-7B、TaSus2-2Btgw、CWI22、CWI21和TKX3D等相應(yīng)的功能標記, 可快速準確鑒定不同小麥品種中控制粒重的相關(guān)基因及其變異類型的分布情況。相吉山等[23]和劉永偉等[24]利用CWI22、CWI21標記分別對新疆小麥品種資源和黃淮麥區(qū)的小麥品種資源的等位變異類型進行了分子檢測, 均發(fā)現(xiàn)具有基因型材料的千粒重顯著高于類型, 進一步驗證了是粒重的優(yōu)異等位變異。Lu等[25]在普通小麥3DS染色體上電子克隆獲得了細胞分裂素氧化酶基因, 并根據(jù)其序列開發(fā)共顯性功能標記TKX3D, 在不同小麥品種能夠分子檢測出和變異類型, 發(fā)現(xiàn)利用該標記能區(qū)分小麥粒重的高低, 可滿足小麥粒重高低的分子標記輔助選擇需求。之后, 簡大為等[26]利用上述功能標記揭示了新疆小麥改良品種與地方品種的遺傳變異, 發(fā)現(xiàn)高千粒重等位變異在新疆改良品種中分布頻率明顯高于地方品種, 而且大部分優(yōu)異等位變異分布頻率隨著育種時期的推進呈現(xiàn)不連續(xù)性上升的趨勢。本研究從基因型的角度來分析黃淮地區(qū)小麥品種(系)粒重的變化情況, 結(jié)果在所有參試材料中共發(fā)現(xiàn)7種不同的等位變異組合, 且具有不同等位變異組合的品種千粒重存在顯著差異。此外, 沒有發(fā)現(xiàn)具有3個低千粒重等位變異組合的品種, 這可能在長期的人工選擇過程中低千粒重的劣勢等位變異類型逐漸被淘汰, 也可能是由于本研究所選擇試驗材料的代表性存在一定的局限性。

    本研究利用小麥粒重功能標記進行分子檢測, 并結(jié)合其粒重表型測定進行分析, 結(jié)果發(fā)現(xiàn)少數(shù)品種基因型與千粒重表型不一致的情況, 如具有2個高千粒重和1個低千粒重等位變異組合品種的千粒重也較低, 這可能是由于粒重除受上述檢測位點控制外, 還存在其他控制小麥粒重基因位點及其相關(guān)基因, 也可能是由于本試驗材料該類型品種數(shù)量較少(僅2份), 不具有代表性。從本研究的粒重表型結(jié)果來看, 黃淮麥區(qū)品種(系)的千粒重變異范圍較大, 且變異系數(shù)高, 遺傳變異非常豐富, 具有較大的遺傳改良潛力, 也為篩選較高粒重的優(yōu)異種質(zhì)資源創(chuàng)造了有利條件。

    4 結(jié)論

    利用功能標記分子鑒定的方法對來自于黃淮麥區(qū)的183份小麥品種(系)檢測, 發(fā)現(xiàn)在、和位點上不同等位變異均會導(dǎo)致小麥千粒重的顯著變化, 而和位點上的等位變異則對小麥粒重的影響更為重要。參試材料中發(fā)現(xiàn)共7種不同的等位變異組合, 不同等位變異組合的品種千粒重存在顯著性差異, 沒有發(fā)現(xiàn)具有3個低千粒重等位變異組合的品種, 同時發(fā)現(xiàn)具有3個高千粒重等位變異組合品種的平均千粒重最高, 是優(yōu)勢基因型組合, 所占比例為52.5%, 具有這些優(yōu)勢基因型組合的優(yōu)異材料將為我國小麥高產(chǎn)育種過程中親本的選配提供有效資源。

    [1] 劉志勇, 王道文, 張愛民, 梁翰文, 呂慧穎, 鄧向東, 葛毅強, 魏珣, 楊維才. 小麥育種行業(yè)創(chuàng)新現(xiàn)狀與發(fā)展趨勢. 植物遺傳資源學(xué)報, 2018, 19: 430–434.

    Liu Z Y, Wang D W, Zhang A M, Liang H W, Lyu H Y, Deng X D, Ge Y Q, Wei X, Yang W C. Current status and perspective of wheat genomics, genetics, and breeding., 2018, 19: 430–434 (in Chinese with English abstract).

    [2] Zanke C D, Ling J, Plieske J, Kollers S, Ebmeyer E, Korzun V, Argillier O, Stiewe G, Hinze M, Neumann F, Eichhorn A, Polley A, Jaenecke C, Ganal M W, R?der M S. Analysis of main effect QTL for thousand grain weight in European winter wheat (L.) by genome-wide association mapping., 2015, 6: 644.

    [3] 宋健民, 戴雙, 李豪圣, 程敦公, 劉愛峰, 曹新有, 劉建軍, 趙振東. 山東省近年來審定小麥品種農(nóng)藝和品質(zhì)性狀演變分析. 中國農(nóng)業(yè)科學(xué), 2013, 46: 1114–1126.

    Song J M, Dai S, Li H S, Cheng D G, Liu A F, Cao X Y, Liu J J, Zhao Z D. Evolution of agronomic and quality traits of wheat cultivars released in Shandong province recently., 2013, 46: 1114–1126 (in Chinese with English abstract).

    [4] 曹廷杰, 趙虹, 王西成, 崔黨群, 詹克慧. 河南省半冬性小麥品種主要農(nóng)藝性狀的演變規(guī)律. 麥類作物學(xué)報, 2010, 30: 439–442.

    Cao T J, Zhao H, Wang X C, Cui D Q, Zhan K H. Evolution of main agronomic traits for semi winter wheat varieties in Henan province., 2010, 30: 439–442 (in Chinese with English abstract).

    [5] Ma D Y, Yan J, He Z H, Wu L, Xia X C. Characterization of a cell wall invertase gene-on common wheat chromosome 2A and development of functional markers., 2012, 29: 43–52.

    [6] Jiang Y M, Jiang Q Y, Hao C Y, Hou J, Wang L F, Zhang H N, Zhang S N, Chen X H, Zhang X Y. A yield-associated gene, in wheat: its function, selection and evolution in global breeding revealed by haplotype analysis., 2015, 128: 131–143.

    [7] Zhang Y J, Liu J D, Xia X C, He Z H.-, an ortholog of rice, is associated with grain weight and grain length in common wheat., 2014, 34: 1097–1107.

    [8] Wang S S, Zhang X F, Chen F, Cui D Q. A single-nucleotide polymorphism ofgene revealed its association with kernel weight in Chinese bread wheat., 2015, 6: 1166.

    [9] Wang S S, Yan X F, Wang Y Y, Liu H M, Cui D Q, Chen F. Haplotypes of thegene are its association with high thousand-kernel weight in Chinese bread wheat., 2016, 7: 783.

    [10] Hanif M, Gao F M, Liu J D, Wen W E, Zhang Y J, Rasheed A, Xia X C, He Z H, Cao S H., an ortholog of rice, is associated with grain weight and yield in bread wheat., 2016, 36: 1?8.

    [11] Hu M J, Zhang H P, Cao J J, Zhu X P, Wang S X, Jiang H, Wu Z Y, Lu J, Chang C, Sun G L, Ma C X. Characterization of an IAA-glucose hydrolase geneassociated with grain weight in common wheat (L.)., 2016, 36: 1–11.

    [12] Yang J, Zhou Y J, Wu Q H, Chen Y X, Zhang P P, Zhang Y, Hu W G, Wang X C, Zhao H, Dong L L, Han J, Liu Z Y, Cao T J. Molecular characterization of a novelallele and its association with grain length in wheat (L.)., 2019, 132: 1799–1814.

    [13] Ma M, Zhao H X, Li Z J, Hu S W, Song W N, Liu X L.regulates seed size in wheat ()., 2015, 67: 1397–1410.

    [14] 司文潔, 吳林楠, 郭利建, 周夢蝶, 劉香利, 馬猛, 趙惠賢. 小麥粒重相關(guān)基因功能標記開發(fā)及驗證. 作物學(xué)報, 2019, 45: 1905–1911.

    Si W J, Wu L N, Guo L J, Zhou M D, Liu X L, Ma M, Zhao H X. Development and validation of the functional marker of grain weight-related genein wheat (L.)., 2019, 45: 1905–1911 (in Chinese with English abstract).

    [15] Yan X F, Zhao L, Ren Y, Dong Z D, Cui D Q, Chen F. Genome-wide association study revealed that thegene was associated with kernel size in Chinese bread wheat., 2019, 9: 2702.

    [16] 寇程, 高欣, 李立群, 李揚, 王中華, 李學(xué)軍. 小麥粒重基因等位變異的組成分析及育種選擇. 作物學(xué)報, 2015, 41: 1640–1647.

    Kou C, Gao X, Li L Q, Li Y, Wang Z H, Li X J. Composition and selection ofalleles for wheat kernel weight., 2015, 41: 1640–1647 (in Chinese with English abstract).

    [17] 時佳, 白璐, 任毅, 穆培源, 梁曉東, 瑪依拉, 耿洪偉. 新疆小麥、、等位變異對粒重的影響及應(yīng)用. 分子植物育種, 2018, 16: 848–858.

    Shi J, Bai L, Ren Y, Mu P Y, Liang X D, Ma Y L, Geng H W. Effects and application of allelic variation of,andon grain weight of Xinjiang wheat., 2018, 16: 848–858 (in Chinese with English abstract).

    [18] Chen F, Zhang F Y, Xia X C, Dong Z D, Cui D Q. Distribution of puroindoline alleles in bread wheat cultivars of the Yellow and Huai valley of China and discovery of a novel puroindoline a allele without PINA protein., 2012, 29: 371–378.

    [19] Cui F, Ding A M, Li J, Zhao C H, Li X F, Feng D S, Wang X Q, Wang L, Gao J R, Wang H G. Wheat kernel dimensions: how do they contribute to kernel weight at an individual QTL level., 2011, 90: 409–425.

    [20] Wu W, Li C J, Ma B L, Shah F, Liu Y, Liao Y C. Genetic progress in wheat yield and associated traits in China since 1945 and future prospects., 2014, 196: 155–168.

    [21] Rasheed A, Xia X C, Ogbonnaya F, Mahmood T, Zhang Z W, Mujeeb-Kazi A, He Z H. Genome-wide association for grain morphology in synthetic hexaploid wheats using digital imaging analysis., 2014, 14: 128.

    [22] Gupta P K, Rustgi S, Kumar N. Genetic and molecular basis of grain size and grain number and its relevance to grain productivity in higher plants., 2006, 49: 565–571

    [23] 相吉山, 穆培源, 桑偉, 聶迎彬, 徐紅軍, 莊麗, 崔鳳娟, 韓新年, 鄒波. 小麥粒重基因功能標記CWI22、CWI21的驗證及應(yīng)用. 中國農(nóng)業(yè)科學(xué), 2014, 47: 2671–2679.

    Xiang J S, Mu P Y, Sang W, Nie Y B, Xu H J, Zhuang L, Cui F J, Han X N, Zou B. Validation and application of function markers CWI22 and CWI21 ofgene related to wheat kernel weight., 2014, 47: 2671–2679 (in Chinese with English abstract).

    [24] 劉永偉, 周碩, 王雪征, 孫果忠, 朱金永, 韓秋芬, 李春杰, 趙和, 王海波. 粒重基因等位變異在黃淮麥區(qū)小麥品種(系)中的分布及功能分析. 華北農(nóng)學(xué)報, 2017, 32: 131–137.

    Liu Y W, Zhou S, Wang X Z, Sun G Z, Zhu J Y, Han Q F, Li C J, Zhao H, Wang H B. Functional analysis and distribution of allelic variations ofgene related to kernel weight in Yellow and Huai River Valleys facultative wheat zone., 2017, 32: 131–137 (in Chinese with English abstract).

    [25] Lu J, Chang C, Zhang H P, Wang S X, Sun G L, Xiao S H, Ma C X. Identification of a novel allele ofassociated with grain size, filling rate and weight of common wheat., 2015, 10: e0144765.

    [26] 簡大為, 周陽, 劉宏偉, 楊麗, 買春艷, 于立強, 韓新年, 張宏軍, 李洪杰. 利用功能標記揭示新疆小麥改良品種與地方品種的遺傳變異. 作物學(xué)報, 2018, 44: 657–671.

    Jian D W, Zhou Y, Liu H W, Yang L, Mai C Y, Yu L Q, Han X N, Zhang H J, Li H J. Functional markers reveal genetic variations in wheat improved cultivars and landraces from Xinjiang., 2018, 44: 657–671 (in Chinese with English abstract).

    Molecular identification and breeding application of allelic variation of grain weight gene in wheat from the Yellow-Huai-River Valley

    ZHANG Fu-Yan1, CHENG Zhong-Jie1, CHEN Xiao-Jie1, WANG Jia-Huan1, CHEN Feng2, FAN Jia-Lin1, ZHANG Jian-Wei1,*, and YANG Bao-An1,*

    1Isotope Institute Co., Ltd, Henan Academy of Sciences / Henan Key Laboratory of Nuclear Agricultural Sciences, Zhengzhou, 450015, Henan, China;2Agronomy College, Henan Agricultural University / National Key Laboratory of Wheat and Maize Crop Science in Henan Province, Zhengzhou, 450002, Henan, China

    Grain weight is one of the most important yield traits in wheat. To accelerate the application process of dominant allelic variation of grain weight gene in wheat breeding, the allelic variation of grain weight gene was identified by using functional markers and the combinations of dominant grain weight genotypes were investigated the allelic variations of grain weight genes,, andin 183 wheat varieties (lines) from the Yellow-Huai-River Valley were identified by PCR amplification with specific primers. To identify the dominant genotype combinations, the effects of different allelic variation genotypes on wheat grain weight were studied by combining the phenotypic data of the 1000-grain weight (TGW) from 2016–2017 and 2017–2018. The results showed that the difference of TGW between different years was highly significant at< 0.01. Two alleles,and, were detected atlocus, with the frequencies of 66.7% and 33.3%, res-pectively. The frequency ofallele onlocus was up to 94.5%, while the frequency ofallele was only 5.5%. In addition, two alleles,and, were found inlocus, and their frequencies were 79.8% and 20.2%, respectively. Further results indicated that there were significant differences in TGW among different allelic variation combinations at< 0.05. Among them, the average TGW of varieties withgenotype was the highest. There was not significant difference in TGW betweenandgenotype, but it was significantly higher than other genotypes at< 0.05. The average TGW ofgenotype was the lowest. The allelic variations at the loci of,,andall led to significant changes in TGW, and the allelic variations ofandloci were more important to TGW in wheat. There were no varieties with three low TGW allelic variation combinationsin the tested materials. Among the seven different allele combinations, the average TGW of three high TGW allelic variation combinationswas the highest, which was the dominant genotype combination.

    wheat; grain weight gene; functional marker; allelic variation; TGW

    10.3724/SP.J.1006.2021.01083

    本研究由河南省小麥產(chǎn)業(yè)技術(shù)體系建設(shè)專項(Z2010-01-04), 河南省重點研發(fā)與推廣專項(科技攻關(guān))(202102110028), 河南省科學(xué)院基本科研業(yè)務(wù)專項(190604015)和省部共建小麥玉米作物學(xué)國家重點實驗室開放課題(30500884)資助。

    This study was supported by the Henan Wheat Research System (Z2010-01-04), the Key Research and Development and Promotion Program of Henan Province (202102110028), the Basic Scientific Research Project of Henan Academy of Sciences (190604015), and the National Key Laboratory of Wheat and Maize Crop Science (30500884).

    張建偉, E-mail: zjw10308@163.com; 楊保安, E-mail: yangcorn@163.com

    E-mail: zhangfuyan704@163.com

    2020-11-03;

    2021-03-19;

    (網(wǎng)絡(luò)出版日期): 2021-04-14.

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

    猜你喜歡
    等位粒重克隆
    基于GBS測序和連鎖分析的藜麥單株粒重QTL定位
    種子(2023年9期)2023-11-22 13:10:56
    克隆狼
    芥菜種子顏色調(diào)控基因TT8的等位變異及其地理分布分析
    ·術(shù)語解析·
    離體穗培養(yǎng)條件下C、N供給對小麥穗粒數(shù)、粒重及蛋白質(zhì)含量的影響
    浙江:誕生首批體細胞克隆豬
    玉米自交系京92改良后代單穗粒重的雜種優(yōu)勢研究
    外引大麥農(nóng)藝性狀SSR關(guān)聯(lián)位點及等位變異表型效應(yīng)分析
    花時間在餐廳門口等位值嗎
    抗BP5-KLH多克隆抗體的制備及鑒定
    后天国语完整版免费观看| 国产99白浆流出| 亚洲国产精品999在线| 国产三级黄色录像| 在线观看午夜福利视频| 99热这里只有精品一区 | 中文字幕人成人乱码亚洲影| 搡老妇女老女人老熟妇| 国产单亲对白刺激| 亚洲成av人片在线播放无| 国产一级毛片七仙女欲春2| av欧美777| 偷拍熟女少妇极品色| 18禁黄网站禁片午夜丰满| 免费在线观看日本一区| 男人舔女人的私密视频| 欧美性猛交╳xxx乱大交人| 两性夫妻黄色片| 国产精品久久电影中文字幕| av天堂中文字幕网| 国产一区二区在线观看日韩 | 18美女黄网站色大片免费观看| 精品久久蜜臀av无| 国产毛片a区久久久久| 午夜影院日韩av| 90打野战视频偷拍视频| 三级毛片av免费| 免费在线观看日本一区| 国产精品久久久久久亚洲av鲁大| av在线天堂中文字幕| 非洲黑人性xxxx精品又粗又长| 日日夜夜操网爽| 不卡一级毛片| 色综合站精品国产| x7x7x7水蜜桃| 国产亚洲精品一区二区www| www.自偷自拍.com| 偷拍熟女少妇极品色| aaaaa片日本免费| 国产成人精品久久二区二区免费| 男人舔奶头视频| 国产亚洲欧美在线一区二区| 精品无人区乱码1区二区| 国产av麻豆久久久久久久| 亚洲一区二区三区色噜噜| 国产私拍福利视频在线观看| 亚洲色图 男人天堂 中文字幕| 两个人看的免费小视频| 午夜福利高清视频| 俄罗斯特黄特色一大片| 成人亚洲精品av一区二区| 成人高潮视频无遮挡免费网站| 少妇人妻一区二区三区视频| 日本免费a在线| 两人在一起打扑克的视频| 9191精品国产免费久久| 国产又色又爽无遮挡免费看| cao死你这个sao货| 久久午夜综合久久蜜桃| 99久久久亚洲精品蜜臀av| 波多野结衣高清作品| 性色avwww在线观看| 国产又黄又爽又无遮挡在线| 999久久久国产精品视频| 亚洲中文av在线| 欧美日韩一级在线毛片| 好看av亚洲va欧美ⅴa在| 欧美日本亚洲视频在线播放| 舔av片在线| 亚洲欧美激情综合另类| 色老头精品视频在线观看| 国产一区二区在线观看日韩 | av天堂在线播放| 九色成人免费人妻av| 亚洲国产色片| 搡老熟女国产l中国老女人| 午夜影院日韩av| 午夜亚洲福利在线播放| 特大巨黑吊av在线直播| 国产精品98久久久久久宅男小说| 亚洲精品色激情综合| 真人一进一出gif抽搐免费| 中文字幕久久专区| 99国产极品粉嫩在线观看| 亚洲九九香蕉| 色av中文字幕| 午夜福利在线观看免费完整高清在 | 国产一区二区三区视频了| 岛国视频午夜一区免费看| 给我免费播放毛片高清在线观看| 精品久久久久久久人妻蜜臀av| 欧美+亚洲+日韩+国产| 亚洲av电影在线进入| 久久久久性生活片| 岛国在线观看网站| 久久中文看片网| 真人一进一出gif抽搐免费| 99久久精品一区二区三区| 伊人久久大香线蕉亚洲五| 久久久精品欧美日韩精品| 最新在线观看一区二区三区| 亚洲欧美精品综合一区二区三区| 亚洲av成人av| 欧美3d第一页| 日韩中文字幕欧美一区二区| 久久久久久国产a免费观看| 国产成人av教育| 黑人操中国人逼视频| 欧美日本亚洲视频在线播放| 中亚洲国语对白在线视频| 国产aⅴ精品一区二区三区波| 成年女人看的毛片在线观看| a级毛片a级免费在线| 日本撒尿小便嘘嘘汇集6| 超碰成人久久| 午夜免费观看网址| 无人区码免费观看不卡| 色老头精品视频在线观看| 日本 av在线| 欧美日韩中文字幕国产精品一区二区三区| 91av网站免费观看| 久久人妻av系列| 欧美乱码精品一区二区三区| 免费看光身美女| 久久国产乱子伦精品免费另类| 深夜精品福利| 一个人观看的视频www高清免费观看 | 黄频高清免费视频| 日韩三级视频一区二区三区| 国产精品一及| 一个人观看的视频www高清免费观看 | 高清毛片免费观看视频网站| 男人的好看免费观看在线视频| 午夜精品在线福利| 后天国语完整版免费观看| 国产成人av教育| 久久久久久久久免费视频了| 亚洲午夜理论影院| 国产精品久久久久久久电影 | 午夜福利欧美成人| 亚洲av成人av| 亚洲精华国产精华精| 日本撒尿小便嘘嘘汇集6| xxx96com| 特级一级黄色大片| 久久久久久九九精品二区国产| 日韩欧美一区二区三区在线观看| 听说在线观看完整版免费高清| 母亲3免费完整高清在线观看| 天堂影院成人在线观看| 日本五十路高清| 亚洲欧美日韩卡通动漫| 免费在线观看日本一区| 99re在线观看精品视频| 中国美女看黄片| 国产av麻豆久久久久久久| 亚洲乱码一区二区免费版| 啪啪无遮挡十八禁网站| 小蜜桃在线观看免费完整版高清| 丰满的人妻完整版| xxxwww97欧美| 欧美三级亚洲精品| 黄色视频,在线免费观看| 精品99又大又爽又粗少妇毛片 | 美女高潮的动态| 变态另类成人亚洲欧美熟女| 黄色视频,在线免费观看| 国产成人啪精品午夜网站| 精品国产亚洲在线| 一夜夜www| 一本精品99久久精品77| 亚洲美女黄片视频| 99在线视频只有这里精品首页| 久久久久久久久中文| 亚洲国产看品久久| 热99re8久久精品国产| 欧美精品啪啪一区二区三区| 一二三四社区在线视频社区8| 首页视频小说图片口味搜索| 色哟哟哟哟哟哟| 免费大片18禁| 日本一本二区三区精品| 日本五十路高清| 久久天躁狠狠躁夜夜2o2o| 欧美最黄视频在线播放免费| 97人妻精品一区二区三区麻豆| 国产成人影院久久av| 色综合婷婷激情| 中文字幕熟女人妻在线| 成人鲁丝片一二三区免费| 黄色成人免费大全| 最近视频中文字幕2019在线8| 亚洲欧美精品综合久久99| 日韩中文字幕欧美一区二区| 国产成人精品无人区| 国产成人av激情在线播放| 午夜影院日韩av| 亚洲精品久久国产高清桃花| 无限看片的www在线观看| 久久天堂一区二区三区四区| 精品久久蜜臀av无| 国产精品久久久久久亚洲av鲁大| 香蕉国产在线看| 老司机深夜福利视频在线观看| 丰满人妻一区二区三区视频av | 亚洲国产精品久久男人天堂| 国产精品影院久久| 久久久色成人| 亚洲欧美精品综合一区二区三区| 我要搜黄色片| 脱女人内裤的视频| 丰满的人妻完整版| 精品午夜福利视频在线观看一区| 国产精品综合久久久久久久免费| 性色avwww在线观看| 91在线观看av| 波多野结衣高清作品| 亚洲精品456在线播放app | 91麻豆av在线| 日本黄大片高清| 日本a在线网址| 国产探花在线观看一区二区| 国产精品久久久久久精品电影| 精品欧美国产一区二区三| 看免费av毛片| 又粗又爽又猛毛片免费看| 免费大片18禁| 精品久久久久久久久久久久久| 国产视频一区二区在线看| 久久天堂一区二区三区四区| 午夜免费激情av| 真实男女啪啪啪动态图| 亚洲人成网站在线播放欧美日韩| 91麻豆精品激情在线观看国产| 亚洲国产精品成人综合色| 国产单亲对白刺激| 婷婷亚洲欧美| 国产精品香港三级国产av潘金莲| 国产黄色小视频在线观看| av片东京热男人的天堂| 国产精品99久久久久久久久| 此物有八面人人有两片| 精品久久久久久久毛片微露脸| 久久久久性生活片| 亚洲午夜理论影院| 日韩欧美 国产精品| 美女高潮喷水抽搐中文字幕| 色精品久久人妻99蜜桃| 欧美日韩亚洲国产一区二区在线观看| 网址你懂的国产日韩在线| 欧美一区二区国产精品久久精品| 日本一本二区三区精品| 免费av毛片视频| 夜夜夜夜夜久久久久| 午夜福利在线观看免费完整高清在 | 亚洲乱码一区二区免费版| 一本一本综合久久| 一个人看的www免费观看视频| 久久久久久大精品| 啦啦啦观看免费观看视频高清| 午夜a级毛片| 在线看三级毛片| 国产99白浆流出| 啪啪无遮挡十八禁网站| 国内精品一区二区在线观看| 精品国产乱码久久久久久男人| 中文字幕久久专区| 亚洲精品色激情综合| 亚洲avbb在线观看| 一级a爱片免费观看的视频| 一级作爱视频免费观看| 一区二区三区激情视频| 亚洲精华国产精华精| 国产精品亚洲美女久久久| 国产成人福利小说| 91老司机精品| 国产精品免费一区二区三区在线| www国产在线视频色| 欧美+亚洲+日韩+国产| 免费观看的影片在线观看| 757午夜福利合集在线观看| a级毛片a级免费在线| 精品一区二区三区av网在线观看| 国产私拍福利视频在线观看| 国产av一区在线观看免费| 久久国产乱子伦精品免费另类| 熟女人妻精品中文字幕| 黄色 视频免费看| 母亲3免费完整高清在线观看| 国产亚洲精品一区二区www| 亚洲国产欧美人成| 怎么达到女性高潮| 波多野结衣高清作品| 精品久久久久久久久久久久久| 美女免费视频网站| 国产单亲对白刺激| 天天躁狠狠躁夜夜躁狠狠躁| 欧美丝袜亚洲另类 | 最新在线观看一区二区三区| 精品人妻1区二区| 黄片大片在线免费观看| 婷婷亚洲欧美| 免费人成视频x8x8入口观看| 欧美激情久久久久久爽电影| 黄频高清免费视频| 99国产精品一区二区三区| 久久久久久久午夜电影| 99精品在免费线老司机午夜| www.999成人在线观看| 久久久久精品国产欧美久久久| 欧美在线黄色| 狂野欧美激情性xxxx| 久久久久久人人人人人| 最近在线观看免费完整版| 在线观看免费视频日本深夜| 性色avwww在线观看| 国语自产精品视频在线第100页| 国内毛片毛片毛片毛片毛片| 国产亚洲精品一区二区www| 视频区欧美日本亚洲| 日韩大尺度精品在线看网址| 免费高清视频大片| 99热精品在线国产| 19禁男女啪啪无遮挡网站| 黑人欧美特级aaaaaa片| 亚洲av免费在线观看| av黄色大香蕉| 国产精品久久久久久精品电影| 熟女少妇亚洲综合色aaa.| 无遮挡黄片免费观看| 很黄的视频免费| 精品电影一区二区在线| 一进一出抽搐动态| 国产三级中文精品| 欧美色欧美亚洲另类二区| 97碰自拍视频| 99久久精品国产亚洲精品| 91字幕亚洲| 别揉我奶头~嗯~啊~动态视频| 不卡一级毛片| 久久久国产成人免费| 身体一侧抽搐| 麻豆国产av国片精品| 欧美日韩黄片免| 男插女下体视频免费在线播放| 精品福利观看| www.www免费av| 黄色成人免费大全| 久久久久九九精品影院| 国产精品久久久久久亚洲av鲁大| 午夜精品久久久久久毛片777| 又紧又爽又黄一区二区| 午夜免费激情av| 亚洲成人免费电影在线观看| 免费观看的影片在线观看| 日韩中文字幕欧美一区二区| tocl精华| 国产极品精品免费视频能看的| 禁无遮挡网站| 色哟哟哟哟哟哟| 后天国语完整版免费观看| 不卡一级毛片| 亚洲成人久久爱视频| 午夜福利高清视频| 免费在线观看成人毛片| 色在线成人网| 成熟少妇高潮喷水视频| 热99在线观看视频| 午夜久久久久精精品| 丰满的人妻完整版| 99热这里只有是精品50| 久久精品影院6| 一级毛片女人18水好多| 亚洲精品国产精品久久久不卡| 精品国内亚洲2022精品成人| 久久久久久久午夜电影| 好看av亚洲va欧美ⅴa在| 免费看美女性在线毛片视频| 欧美激情久久久久久爽电影| 国模一区二区三区四区视频 | 国产成人av教育| 欧美+亚洲+日韩+国产| 久久这里只有精品19| 国产精品久久久久久亚洲av鲁大| 女人被狂操c到高潮| 97人妻精品一区二区三区麻豆| 精品久久蜜臀av无| 国产伦精品一区二区三区视频9 | 中文资源天堂在线| 我要搜黄色片| 欧美日韩精品网址| 高清毛片免费观看视频网站| 我的老师免费观看完整版| 欧美日韩一级在线毛片| 精品人妻1区二区| 欧美3d第一页| 热99在线观看视频| 国产在线精品亚洲第一网站| 女生性感内裤真人,穿戴方法视频| 成人精品一区二区免费| 观看美女的网站| 亚洲最大成人中文| 一边摸一边抽搐一进一小说| 亚洲电影在线观看av| 久久久国产成人免费| 18禁观看日本| 女人被狂操c到高潮| 啦啦啦免费观看视频1| www.精华液| 他把我摸到了高潮在线观看| 一级毛片女人18水好多| 嫩草影院精品99| 一本综合久久免费| 一个人观看的视频www高清免费观看 | 宅男免费午夜| 窝窝影院91人妻| 亚洲av日韩精品久久久久久密| 在线观看舔阴道视频| 亚洲av第一区精品v没综合| 成人国产综合亚洲| 桃红色精品国产亚洲av| 一区二区三区高清视频在线| 99久久精品国产亚洲精品| 亚洲成人中文字幕在线播放| av黄色大香蕉| 国产精品九九99| 中文资源天堂在线| 啪啪无遮挡十八禁网站| 亚洲 欧美 日韩 在线 免费| 无人区码免费观看不卡| 精华霜和精华液先用哪个| 成人18禁在线播放| АⅤ资源中文在线天堂| 91九色精品人成在线观看| 性色av乱码一区二区三区2| 日韩精品中文字幕看吧| 午夜免费激情av| 一进一出抽搐gif免费好疼| 99热6这里只有精品| 在线视频色国产色| 国产高清videossex| 长腿黑丝高跟| 国产午夜精品论理片| 两个人视频免费观看高清| 网址你懂的国产日韩在线| 青草久久国产| 日本五十路高清| 亚洲av电影不卡..在线观看| 在线免费观看的www视频| 成人一区二区视频在线观看| 欧美性猛交╳xxx乱大交人| 岛国在线观看网站| 国产亚洲av嫩草精品影院| 大型黄色视频在线免费观看| 91麻豆av在线| 亚洲av免费在线观看| 麻豆一二三区av精品| 欧美日韩乱码在线| 亚洲人成网站在线播放欧美日韩| 欧美+亚洲+日韩+国产| 女警被强在线播放| 欧美av亚洲av综合av国产av| 中文字幕最新亚洲高清| 女生性感内裤真人,穿戴方法视频| 欧美日韩国产亚洲二区| 欧美xxxx黑人xx丫x性爽| 悠悠久久av| 美女高潮的动态| 久久精品国产亚洲av香蕉五月| 精品国产乱子伦一区二区三区| 久久久久久大精品| 激情在线观看视频在线高清| 久久九九热精品免费| 国产精品影院久久| 国产精品久久电影中文字幕| 老汉色av国产亚洲站长工具| 搞女人的毛片| 免费观看的影片在线观看| 亚洲av成人不卡在线观看播放网| 色综合欧美亚洲国产小说| 18禁国产床啪视频网站| 波多野结衣巨乳人妻| 欧美成人性av电影在线观看| 欧美日韩黄片免| 夜夜夜夜夜久久久久| 九九久久精品国产亚洲av麻豆 | 免费电影在线观看免费观看| 一级毛片精品| 精品免费久久久久久久清纯| 成年女人看的毛片在线观看| 夜夜夜夜夜久久久久| 久久久成人免费电影| 免费一级毛片在线播放高清视频| 亚洲欧美一区二区三区黑人| 精品国产超薄肉色丝袜足j| 欧美色欧美亚洲另类二区| 久久精品aⅴ一区二区三区四区| 一级毛片女人18水好多| 亚洲精品中文字幕一二三四区| 一区二区三区激情视频| 成年免费大片在线观看| 国产精品,欧美在线| 麻豆一二三区av精品| 亚洲欧美日韩卡通动漫| 两个人的视频大全免费| 色综合欧美亚洲国产小说| 五月玫瑰六月丁香| 国产午夜精品论理片| 99riav亚洲国产免费| 此物有八面人人有两片| 国产精品永久免费网站| a级毛片a级免费在线| 国语自产精品视频在线第100页| 99精品久久久久人妻精品| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久久久久久久久免费视频| 免费看美女性在线毛片视频| 日本黄色视频三级网站网址| 亚洲一区高清亚洲精品| 亚洲av美国av| 欧美日韩国产亚洲二区| 免费观看的影片在线观看| 国产野战对白在线观看| 男女那种视频在线观看| 精品国产乱码久久久久久男人| 全区人妻精品视频| 国产69精品久久久久777片 | 亚洲一区二区三区不卡视频| 亚洲五月婷婷丁香| 最好的美女福利视频网| 国产高清三级在线| 欧美日韩乱码在线| 三级国产精品欧美在线观看 | 黄片小视频在线播放| 91九色精品人成在线观看| 18禁裸乳无遮挡免费网站照片| 日本五十路高清| 国产综合懂色| 观看美女的网站| 身体一侧抽搐| 欧美日韩亚洲国产一区二区在线观看| 国产成人福利小说| 国产精品一区二区精品视频观看| 非洲黑人性xxxx精品又粗又长| 国产精品久久久人人做人人爽| 看黄色毛片网站| 丝袜人妻中文字幕| 2021天堂中文幕一二区在线观| av黄色大香蕉| 少妇裸体淫交视频免费看高清| 精品一区二区三区视频在线观看免费| 免费看日本二区| 五月玫瑰六月丁香| 岛国在线观看网站| 女警被强在线播放| 色综合亚洲欧美另类图片| 国产亚洲精品av在线| 桃红色精品国产亚洲av| 国产成人系列免费观看| 99精品在免费线老司机午夜| 亚洲国产欧美一区二区综合| 国产av在哪里看| 黄色 视频免费看| 人妻久久中文字幕网| 国产一区二区激情短视频| 手机成人av网站| 国产淫片久久久久久久久 | 欧美不卡视频在线免费观看| 一进一出抽搐动态| 欧美午夜高清在线| 欧美绝顶高潮抽搐喷水| 在线观看日韩欧美| 国产精品久久久人人做人人爽| 国产熟女xx| 给我免费播放毛片高清在线观看| 国产真人三级小视频在线观看| 天堂av国产一区二区熟女人妻| 欧美成人一区二区免费高清观看 | 亚洲欧洲精品一区二区精品久久久| or卡值多少钱| 夜夜看夜夜爽夜夜摸| 免费在线观看成人毛片| 18禁国产床啪视频网站| 超碰成人久久| 操出白浆在线播放| 热99在线观看视频| 亚洲无线观看免费| 精品一区二区三区视频在线 | 国产成人影院久久av| 亚洲国产精品sss在线观看| 国产欧美日韩精品一区二区| 成人无遮挡网站| 久久人人精品亚洲av| 脱女人内裤的视频| 久久天躁狠狠躁夜夜2o2o| 亚洲自拍偷在线| 亚洲 欧美一区二区三区| 中国美女看黄片| 成人三级黄色视频| 天堂av国产一区二区熟女人妻| a在线观看视频网站| 欧美在线黄色| 99久久99久久久精品蜜桃| 香蕉国产在线看| 一级a爱片免费观看的视频| 波多野结衣巨乳人妻| 青草久久国产| 久久亚洲精品不卡| 国产伦一二天堂av在线观看| 欧美色视频一区免费| 丰满的人妻完整版| 不卡av一区二区三区| 久久这里只有精品19| 波多野结衣高清作品| 变态另类成人亚洲欧美熟女|