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

    普通小麥粒重相關基因克隆的研究進展

    2017-02-05 02:29:24呂廣德孫憲印亓曉蕾王瑞霞孫盈盈牟秋煥錢兆國
    麥類作物學報 2017年10期
    關鍵詞:關聯(lián)產(chǎn)量

    呂廣德,孫憲印,亓曉蕾,王 超,王瑞霞,孫盈盈,牟秋煥,米 勇,錢兆國,吳 科

    (泰安市農(nóng)業(yè)科學研究院,山東泰安 271000)

    普通小麥粒重相關基因克隆的研究進展

    呂廣德,孫憲印,亓曉蕾,王 超,王瑞霞,孫盈盈,牟秋煥,米 勇,錢兆國,吳 科

    (泰安市農(nóng)業(yè)科學研究院,山東泰安 271000)

    隨著分子生物學特別是生物信息學和基因組學的發(fā)展,小麥粒重相關基因的克隆也取得了很大進展。本文主要介紹了與小麥粒重有關的 TaGS5、 TaGASR7、 TaSus、 TaWSC、 Ta-6-SFT、 TaGW2、 (TaCKX、) TaTGW6、 TaCWI、 TaGS1基因以及TEF、SAP等家族基因 TaTEF、 TaSAP、 TaSnRK2s的克隆及其功能研究的進展,以期為小麥高產(chǎn)育種提供參考。

    普通小麥;粒重相關基因;克隆

    小麥作為重要的糧食作物,和水稻、玉米一起提供了人類50%的能量。隨著全球人口增長,預計到2020年對小麥的需求量將增加40%[1]。小麥產(chǎn)量要素包括穗粒數(shù)、穗數(shù)和千粒重,其中,千粒重的遺傳力較高,證明其是最穩(wěn)定的產(chǎn)量構(gòu)成要素[2]。小麥粒重與籽粒大小、形狀呈正相關關系,而籽粒形狀是由粒長、粒寬、籽粒長寬比和粒厚共同組成[3]。小麥粒重受很多帶有加性效應和上位性效應的基因所控制,也受環(huán)境因素的影響[3]。控制籽粒大小和重量的復雜遺傳因素促使了小麥染色體上粒重相關基因/QTLs的發(fā)現(xiàn)。為了更好地將產(chǎn)量相關基因應用于小麥分子育種中,小麥粒重基因克隆的研究一直是人們研究的熱點和重點。近年來,隨著分子生物學特別是生物信息學和基因組學的發(fā)展,小麥粒重相關基因的克隆取得了很大進展。在本文中,我們對當前已克隆的與粒重有關的基因加以闡述,以期為小麥高產(chǎn)育種提供參考。

    1 小麥籽粒大小和重量的遺傳研究進展

    小麥的馴化是在10 000年以前開始的,其和其他農(nóng)作物最重要的區(qū)別之一是籽粒大小和重量的改變[3]。籽粒大小(粒長、粒寬、粒厚、籽粒密度等)是小麥產(chǎn)量因素中的一個重要組成部分,較大的籽粒不僅和產(chǎn)量有直接關系,也對籽?;盍驮缙谏L有非常大的影響,從而促進和穩(wěn)定小麥本身的生產(chǎn)能力。幾何模型顯示,籽粒形狀、大小的改變能使面粉產(chǎn)量增加5%左右[4]。籽粒的長度、寬度、厚度和籽粒密度等性狀也是小麥形態(tài)品質(zhì)檢測的重要指標[5]。在小麥育種過程中,較大的籽粒被千粒重所代表,千粒重的遺傳力較高證明其是最穩(wěn)定的產(chǎn)量構(gòu)成要素[2]。在小麥產(chǎn)量的構(gòu)成要素中,千粒重的增加是相對獨立的。Breseghello等[6-7]對S×O、R×G兩個作圖群體研究發(fā)現(xiàn),小麥產(chǎn)量與籽粒寬/長之比存在極顯著相關關系,而與籽粒厚/長之比無顯著相關關系;而粒型的各指標間緊密相關,其中,千粒重與籽粒體積的相關性最高(r=0.95和0.91),粒長與粒寬的相關性卻不大(r=0.30和0.22),表明兩者在遺傳上是相對獨立的性狀。王瑞霞等[8]對粒型各性狀間的相關性進行研究表明,不同性狀間的相關關系都達到了(極)顯著水平,其中,籽粒體積和粒重的相關系數(shù)最高,而粒長和粒寬、粒厚的相關系數(shù)相對較低。Sun等[9]把密切相關的粒型和粒重性狀一起進行分析表明,粒長與籽寬、粒重呈顯著正相關,粒長與容重呈為負相關,其中,籽粒寬度和粒重的相關度最高,但與籽粒容重關系不大,而粒重與籽粒容重呈正相關關系。Gegas等[3]對6個DH群體中粒型各性狀進行研究發(fā)現(xiàn),千粒重與籽粒體積、粒寬和籽粒密度都存在極顯著正相關關系(r≥0.75,P<0.001),但與粒長的相關性較弱(r≥0.23,P<0.001)。

    2 粒重相關基因克隆的研究進展

    2.1 TaGS5基因的克隆

    籽粒大小是粒重的重要組成部分。在水稻中研究表明, OsGS5通過促進細胞分裂來調(diào)控籽粒大小和重量。Wang等[10]利用電子克隆技術在六倍體小麥中克隆了TaGS5基因并將其定位在3AS和3DS染色體上,通過對該基因的序列分析發(fā)現(xiàn) TaGS5-A1基因的第六外顯子有一個SNP位點。與農(nóng)藝性狀的關聯(lián)分析結(jié)果顯示,等位基因 TaGS5-A1a在株高、穗長、穗下節(jié)間長、籽粒長寬比等性狀上比等位基因 TaGS5-A1b的表達水平高,而在千粒重、粒寬方面, TaGS5-A1b的貢獻更高[10]。另外, TaGS5-A1的表達分析表明,在籽粒發(fā)育的不同階段, TaGS5-A1b比 TaGS5-A1a的表達水平高[2]。Ma等[11]在3A/3B/3D染色體上也克隆到了 TaGS5基因,單倍型分析以及轉(zhuǎn)基因試驗分析結(jié)果顯示, TaGS5-3A對籽粒大小是正向調(diào)控的,優(yōu)異等位基因 TaGS5-3A-T對小麥高產(chǎn)育種起到一個積極地推動作用。

    2.2 TaGASR7基因的克隆

    在高等植物中, GASR7(gibberelinacid stimulated related 7)基因是赤霉素調(diào)控基因家族中的一員。在水稻和小麥中,該基因與籽粒長度密切關聯(lián)。Ling等[12]對烏拉爾圖小麥的全基因組草圖繪制過程中發(fā)現(xiàn)一個赤霉素調(diào)控基因,并命名為 TaGASR7,該基因與水稻中的赤霉素調(diào)控基因 OsGASR7同源。Dong等[13]在小麥中克隆了三個 GASR7同源基因,分別命名為 TaGASR7-A1、 TaGASR7-B1、 TaGASR7-D1。在 TaGASR7-A1中,單倍型 H1C比 H1G對粒長的貢獻大。Kumar等[14]研究發(fā)現(xiàn),小麥粒長和粒寬在遺傳上是相對獨立的,在其構(gòu)建的ND705×PI414566重組自交系群體中共檢測到29個與粒型和籽粒大小有關的QTL,其中 QTL-27與 TaGASR7基因的位置一致,均位于7AL染色體上。 TaGASR7基因可能對小麥的籽粒發(fā)育和籽粒大小等性狀具有重要作用[14]。

    2.3 TaSus基因的克隆

    小麥開花期和成熟期的籽粒灌漿機制決定著籽粒最后的重量,籽粒灌漿整個過程包含三個生理進程,分別是干物質(zhì)的積累、水分積累和干燥、籽粒形態(tài)膨脹。其中,干物質(zhì)包括淀粉、蛋白質(zhì)以及其他的營養(yǎng)成分,而淀粉占籽粒沉淀物重量的70%[4]。Hou等[15]克隆得到小麥蔗糖合成酶(sucrose synthesis, Sus)基因 TaSus1和 TaSus2(分別位于7A/7B/7D和2A/2B/2D染色體),通過對基因單倍型和千粒重的關聯(lián)分析共檢測到4個與千粒重關聯(lián)密切的單倍型,分別是 TaSus2-2A(Hap-A)、 TaSus2-2B( TaSus2-2B-Hap-H)、 TaSus1-7A(Hap-1和Hap-2)和 TaSus1-7B(Hap-T)。 TaSAP1是一個與非生物脅迫有關的蛋白。Chang等[16]發(fā)現(xiàn)編碼該蛋白的基因位于7A染色體上,與千粒重、穗粒數(shù)、穗下節(jié)長、穗長、小穗數(shù)顯著關聯(lián)。Mohler等[17]研究發(fā)現(xiàn), TaSus1-7A與 TaSAP-A1緊密連鎖,并與位于該染色體上控制粒寬和粒長的QTL緊密連鎖。

    2.4 TaWSC基因的克隆

    小麥營養(yǎng)器官中暫貯存性可溶性碳水化合物(water soluble carbohydrates,WSC)是維系其生存和產(chǎn)量形成的重要代謝物質(zhì),WSC積累轉(zhuǎn)運是被微效多基因所控制[18-21],不僅通過滲透調(diào)節(jié)來緩沖逆境脅迫對小麥的傷害,而且也是小麥籽粒灌漿所需的重要碳源。小麥生育后期營養(yǎng)器官WSC向籽粒的高效轉(zhuǎn)運和再分配對小麥產(chǎn)量形成和粒重增加至關重要[19]。研究者利用不同遺傳背景材料和遺傳圖譜,對小麥WSC積累轉(zhuǎn)運相關性狀進行了QTL定位和遺傳分析。Zhang等[22]研究表明,與WSC關聯(lián)的QTL位于1A、2D、4A、4B、5D、6B、7B和7D染色體上。Yang等[23]發(fā)現(xiàn)了20個與莖稈WSC積累轉(zhuǎn)運相關的QTL,其中,9個QTL具有增效效應,位于1A、1D、2D、7A、7B和7D染色體上;11個QTL具有減效效應,位于1A、2A、2D、3B、4A、5A、6B和7B染色體上。Fu等[23]定位了13個控制籽粒4種可溶性糖含量的加性QTL,位于1B、1D、2B、2D、3B、4A、4B、5A、5D、6B和7B等11條染色體上,可解釋4.97%~53.57%的表型變異。Zhang等[22]鑒定到了對千粒重正向調(diào)控的WSC等位基因位點,分別是 Xcfd17-2D223、Xcfd53-2D263、Xgwm181-3B140、Xgwm181-3B161、Xgwm389-3B116、Xbarc125-3D147、Xgwm358-5D162和 Xgwm537-7B205。Li等[24]在16個WSC優(yōu)異等位基因中檢測出5個單獨對千粒重有貢獻的位點,分別是 Xgwm181-3B187、Xgwm148-2B165、Xgwm261-2D203、Xgwm149-4B153和 Xgwm358-5D162。通過對WSC等位基因位點的遺傳研究,莖稈中的WSC的含量已經(jīng)是小麥育種的一個選擇標準,WSC優(yōu)異等位基因的標記輔助選擇對小麥育種將起到很大的促進作用。

    2.5 Ta-6-SFT基因的克隆

    在高等植物中,果聚糖是一種重要的儲藏性可溶性碳水化合物,6-SFT (sucrose:fructan 6-fructosyltransferase)是果聚糖合成過程中的關鍵酶之一。Kawakami等[25]在普通小麥中克隆到了 6-SFT基因(又名 wft1)。Yue等[26]在 6-SFT-A2基因中檢測到13個SNP/InDel位點,該基因分成三個單倍型( HapⅠ、 HapⅡ和 HapⅢ),通過對 6-SFT基因的核苷酸序列差異性分析和SNP與小麥農(nóng)藝性狀的關聯(lián)分析, HapⅢ的千粒重比 HapⅠ、 HapⅡ的千粒重要高,單倍型 HapⅢ能顯著提高小麥粒重,含有該位點的等位基因在提高粒重方面應該是一個優(yōu)異等位基因。

    2.6 TaGW2基因的克隆

    宿振起等[27]根據(jù)水稻中的粒寬基因 OsGW2,在普通小麥中克隆到 TaGW2基因,該基因位于小麥第6同源染色體上,編碼具有E3泛素連接酶功能的蛋白質(zhì)。根據(jù)大粒型和小粒型 TaGW2-6A等位基因在啟動子區(qū)的序列差異,開發(fā)了CAPS標記,關聯(lián)分析結(jié)果顯示, TaGW2-6A主要影響小麥粒寬和粒重,大粒型等位變異 Hap-6A-A使小麥粒重增加3.1 g[27]。Huang等[28-29]在小麥6AS的 Xgwm786和 Xbarc146位點檢測到影響小麥千粒重的QTL,該QTL與 TaGW2-6A在染色體上處同一位點區(qū)域。Sun等[9]利用重組自交系群體在小麥6AS的 TaGW2-6A位點附近檢測到了影響小麥粒寬和千粒重的QTL。Sun等[30]在小麥的6AS與 TaGW2-6A共定位的位點區(qū)同樣檢測到了影響小麥籽粒粒寬和千粒重的QTL,通過遺傳圖譜定位顯示, TaGW2-6A等位基因與以往檢測到的小麥染色體上的控制粒寬及千粒重的位點相吻合。通過比較定位結(jié)果進一步證明了 TaGW2-6A對小麥粒寬及千粒重的調(diào)控功能。Jaiswal等[31]對 TaGW2-6A的啟動子區(qū)域進行序列差異性分析,在其研究鑒定的5個單倍型中, Hap5和 Hap2比其他的單倍型對粒重有更顯著地影響。Qin等[32]對 TaGW2-6B的單倍型與農(nóng)藝性狀的關聯(lián)分析發(fā)現(xiàn), TaGW2-6B對千粒重的影響比 TaGW2-6A更顯著,在增加粒寬和粒重方面, Hap-6B-1是一個優(yōu)異的等位基因。另外,通過對 TaGW2-6A和 TaGW2-6B的單倍型相互作用分析發(fā)現(xiàn),優(yōu)異單倍型組合 Hap-6A-A/Hap-6B-1對粒重的增加更顯著[32]。最近研究發(fā)現(xiàn), TaGW2-6A與粒寬是負相關的關系,但是在同一個小麥品種中, TaGW2-6B和 TaGW2-6D的轉(zhuǎn)錄水平對粒寬存在正向效應, TaGW2的三個同源基因在籽粒發(fā)育過程中可能有不同的功能,可能是三個基因的平衡作用最終影響籽粒的大小[33]。

    2.7 TaCKX基因的克隆

    細胞分裂素(cytokinin,CK)是一個通過調(diào)控胚乳細胞數(shù)量以及調(diào)節(jié)籽粒灌漿模式來控制小麥籽粒大小和籽粒重量的關鍵激素。細胞分裂素氧化/脫氫酶(cytokinin oxidase/dehydrogenase, CKX)是降解細胞分裂素的主要酶類,因此,深入了解CKX在植物生長發(fā)育及形態(tài)建成等過程中的作用機理,對作物產(chǎn)量、抗性等農(nóng)藝性狀的改良有重要的現(xiàn)實意義。Zhang等[34]在普通小麥中利用同源克隆方法克隆了 TaCKX5基因,并將其定位在3A/3B/3D染色體上。Zhang等[35]根據(jù)水稻 OsCKX2基因在小麥中克隆到 TaCKX6-D1,根據(jù)遺傳連鎖圖譜分析、關聯(lián)分析以及表達分析發(fā)現(xiàn),該基因單倍型 TaCKX6-D1a與粒重有顯著關聯(lián)。Chang等[36]在普通小麥中也克隆得到CKX基因,其中, TaCKX1和 TaCKX4定位在3A染色體上, TaCKX2定位在7A和7B染色體上,而 TaCKX2.1、 TaCKX2.2、 TaCKX3、 TaCKX5和 TaCKX6定位在3DS染色體上。Lu等[37]對RIL群體研究顯示,在不同的環(huán)境條件下,等位變異 TaCKX6a02與籽粒大小、粒重和籽粒灌漿速率有顯著的關聯(lián),解釋了17.1%~38.2%的表型變異;同時,開發(fā)了一個特異標記TKX3D,用來對籽粒大小、粒重和籽粒灌漿速率等性狀進行標記輔助選擇。Chang等[36]還對灌漿期旗葉葉綠素含量和粒重等性狀與CKX基因進行關聯(lián)分析,結(jié)果顯示, TaCKX4與葉綠素含量和粒重顯著關聯(lián),通過遺傳圖譜分析, TaCKX4基因與3AL染色體上的分子標記Xwmc169緊密連鎖,該標記與很多控制粒重和授粉后5~15 d旗葉中葉綠素含量的QTL共分離,這些QTL能解釋這兩種性狀8.9%~22.3%的表型變異。另外,該研究也發(fā)現(xiàn) TaCKX4基因的拷貝數(shù)變化能顯著影響小麥粒重和葉綠素含量[36]。

    2.8 TaTGW6基因的克隆

    TGW6基因編碼一個全新的IAA-葡萄糖水解酶,可顯著提高水稻粒重和產(chǎn)量。Hanif等[38]在小麥3AL染色體上克隆到 TaTGW6-A1基因,并利用開發(fā)出的功能標記 TaTGW6-A1-CAPS鑒定該基因與千粒重和產(chǎn)量的關系,對產(chǎn)量性狀QTL分析發(fā)現(xiàn),該基因位點能解釋17.4%的產(chǎn)量性狀變異;兩個單倍型( TaTGW6-A1a和 TaTGW6-A1b)中, TaTGW6-A1a比 TaTGW6-A1b有更高的粒重和產(chǎn)量,說明 TaTGW6-A1a是一個優(yōu)異的等位基因,在小麥育種過程中可對產(chǎn)量性狀進行標記輔助選擇。

    2.9 TaCWI基因的克隆

    CWI基因主要參與蔗糖的分配、控制細胞的分化、庫器官的建立及籽粒的發(fā)育等過程,是庫器官形成中決定庫強的一個最關鍵酶,并最終影響作物的產(chǎn)量。Ma等[39]克隆了普通小麥 2A 染色體上的細胞壁轉(zhuǎn)化酶基因 TaCwi-A1的全長編碼序列,并根據(jù)其等位變異 TaCwi-A1a和 TaCwi-A1b在第4內(nèi)含子2 727 bp處的SNP差異開發(fā)了一對互補顯性功能標記CWI22和CWI21。韓利明等[40]利用標記CWI22和CWI21對來自21個國家的745份小麥品種進行檢測,結(jié)果發(fā)現(xiàn),CWI22和CWI21能很好地區(qū)分等位變異 TaCwi-A1a和 TaCwi-A1b。Ma等[39]通過遺傳連鎖分析把 TaCwi-A1定位在2AL上,與其共分離的主效QTL能夠解釋4.8%的千粒重變異。通過對千粒重進一步檢測發(fā)現(xiàn),CWI22擴增的402 bp條帶與高粒重關聯(lián),CWI21所擴增的404 bp條帶與低粒重關聯(lián)。Jiang等[41]在4A/5B/5D染色體上也能克隆到 TaCWI基因,其中, TaCWI-4A的啟動子區(qū)域有2個SNP位點, TaCWI-5D基因中有4個SNP位點和兩個Indels位點,最終根據(jù)這兩個基因的差異區(qū)域開發(fā)了兩個CAPS標記CAPS4A和CAPS5D;利用這兩個CAPS標記對378份小麥品種分析發(fā)現(xiàn), TaCWI-5D中的單倍型 Hap-5D-C與高粒重有顯著關聯(lián)。Jiang等[41]研究還發(fā)現(xiàn),在干旱區(qū)域,Hap-4A-C可以產(chǎn)生更多的種子,但是在水澆地區(qū)域,Hap-4A-T可以產(chǎn)生更高的粒重,是更優(yōu)異的單倍型。

    2.10 TaGS1基因的克隆

    谷氨酰胺合成酶(glutamine synthetase, GS)是參與將氨催化合成有機氮第一步的關鍵酶,在小麥中,三個GS基因被克隆,包括 GS1、 GSe和 GSr[42]。 GS1基因在籽粒灌漿過程中能夠?qū)再動員,在水稻中,GS1蛋白質(zhì)含量與籽粒數(shù)量和籽粒大小顯著相關[43-44];在玉米[45-47]和小麥[48-49]中,GS1蛋白質(zhì)和GS1活性的QTL與籽粒產(chǎn)量的QTL存在共位性關系。Guo等[50]克隆得到 TaGS1a基因的gDNA全長序列,并根據(jù)差異位點開發(fā)了一個CAPS標記,關聯(lián)分析結(jié)果顯示,該基因與粒重、粒寬、粒長等農(nóng)藝性狀顯著關聯(lián)。

    2.11 家族基因的克隆

    Zheng等[51]對小麥 TEF(transcript elongation factors)基因家族進行了研究, TaTEF-7A位于7A染色體上,與穗粒數(shù)顯著關聯(lián)。進一步對近等基因系的表型分析得知, TaTEF-7A可以增加籽粒產(chǎn)量,并增加產(chǎn)量相關性狀水平,其中,優(yōu)異單倍型 Hap-7A-3與產(chǎn)量顯著關聯(lián)[51]。以上結(jié)果顯示, TaTEF-7A對穗粒數(shù)來說是一個功能性調(diào)控因子。

    小麥 SAP(stress-associated protein)基因家族也與其千粒重有關。在各種環(huán)境脅迫(如干旱、鹽堿以及凍害)下,Chang等[16]將 TaSAP1-A1基因定位在7A染色體上,基于 TaSAP1-A1啟動子區(qū)的3個多態(tài)性位點 InDel5-1810、 SNP-2606和 InDel39-1637開發(fā)了3個分子標記,分別命名為T7AM5、T7AM2606和T7AM39。關聯(lián)分析表明,以上3個分子標記與穗長、穗下節(jié)間長度、每穗總小穗數(shù)、穗粒數(shù)和千粒重5個農(nóng)藝性狀顯著關聯(lián)[16]。

    蔗糖非發(fā)酵相關蛋白激酶家族2(sucrose non-fermenting1-related protein kinase2, SnRK2)參與逆境條件下植物體內(nèi)逆境信號傳遞和能量代謝,在抵御非生物脅迫和調(diào)節(jié)植物的生長發(fā)育等方面發(fā)揮重要作用。Zhang等[52]克隆到 TaSnRK2.10基因,根據(jù)序列多態(tài)性開發(fā)功能標記并與農(nóng)藝性狀進行關聯(lián)分析,發(fā)現(xiàn) TaSnRK2.10-4A兩種單倍型的千粒重在所有的4個環(huán)境條件下均達到了極顯著差異(P≤0.01),單倍型 Hap-4A-H的千粒重較 Hap-4A-L平均高2.17 g。王 倩等[53]對 TaSnRK2.10基因的啟動子區(qū)域進行序列分析并開發(fā)了4個分子標記,將自然群體的262份材料分為4種單倍型,分別與千粒重、單株穗數(shù)和每穗小穗數(shù)顯著相關或極顯著相關,并檢測出 HapⅡ和 HapⅢ是提高千粒重的優(yōu)異單倍型。

    3 討 論

    3.1 基因發(fā)掘方法有所創(chuàng)新

    現(xiàn)代分子生物學特別是基因組學和生物信息學的迅猛發(fā)展推動了基因發(fā)掘策略和方法的改進[54]。目前,基因發(fā)掘方法主要包含正向遺傳學方法和反向遺傳學方法[55]?;谀繕嘶蛟谌旧w位置的圖位克隆方法和基于連鎖不平衡將標記或候選基因的遺傳變異(等位基因變異)與目標性狀表型聯(lián)系起來的關聯(lián)分析方法屬于從表型到基因的正向遺傳學的范疇[56-57];通過遺傳轉(zhuǎn)化研究基因功能以及所引起的可能表型變異進行的基因克隆屬于從基因到表型的反向遺傳學范疇。隨著擬南芥、水稻等模式植物基因組序列測序完成,推動了比較基因組學的發(fā)展,在小麥中通過同源克隆、插入或刪除、mRNA差異顯示、轉(zhuǎn)座子標簽或RNA干涉[58]等方法發(fā)掘到的基因已經(jīng)有很多。另外,將基因關聯(lián)分析與表達分析[59]或轉(zhuǎn)基因[60]相結(jié)合,同源克隆與利用重組自交系群體定位[61]相結(jié)合,以及整合連鎖圖譜、表達譜和功能互補分析來克隆微效抗性基因/QTL[62]等各種方法進行整合以研究基因功能已成為當前基因發(fā)掘的主要方法。

    作物的很多重要性狀都受微效多基因控制遺傳,因此,提高這些基因/QTL的檢測效率對于基因克隆至關重要?;谏衔恍缘母黜椫餍盎蛐汀镰h(huán)境互作效應而創(chuàng)建的發(fā)育性狀條件QTL定位分析法[63],QTL加性、顯性、上位性主效應及其與環(huán)境互作的基因定位遺傳模型[64-65],數(shù)量性狀基因定位及效應分析的基于混合模型的復合區(qū)間作圖(MCIM)方法[66-67],應用混合分布模型對多環(huán)境試驗數(shù)據(jù)進行聯(lián)合分析以及測驗主基因型與環(huán)境互作及環(huán)境效應的統(tǒng)計方法[68],檢測作物數(shù)量性狀基因與遺傳標記連鎖關系的方差分析法[69],四向雜交設計QTL分析的極大似然方法[70]以及數(shù)量性狀基因的完備區(qū)間作圖方法[71]等極大地提高了QTL的檢測效率。同時,QTLIci-Mapping等計算機軟件的發(fā)明也為QTL定位研究提供了有效的計算工具。

    3.2 建立基因功能研究平臺, 加強特異基因克隆

    在小麥的近緣屬植物中,以水稻為研究對象定位或克隆到的基因最多,而其他植物則較少。從基因的功能驗證研究來看,通過水稻轉(zhuǎn)基因進行功能驗證的數(shù)量較多,而其他作物特別是小麥進行轉(zhuǎn)基因功能驗證的基因屈指可數(shù)[72]。作物的基因可通過轉(zhuǎn)化模式植物如擬南芥、煙草等進行功能驗證,如大豆特有的結(jié)瘤特性可以轉(zhuǎn)結(jié)瘤作物的模式植物蒺藜苜蓿[73]和百脈根[74]進行驗證。但是需要指出的是,在模式植物中的功能驗證并不都能完全代表其在來源作物中的直接驗證,如水稻的候選基因與擬南芥基因的直系同源性在很多情況下都是可行的[75],但水稻與擬南芥的共線性也是有限的。構(gòu)建不同作物的高效轉(zhuǎn)化系統(tǒng)[76]是檢驗基因功能不容忽視的技術平臺,水稻基因發(fā)掘研究所取得的成就是與技術平臺的建立分不開的。因此,小麥應借鑒新理論、新技術和新方法加強其特有基因的發(fā)掘,促進小麥自身特異基因發(fā)掘的快速發(fā)展。

    3.3 小麥粒重基因在育種中的應用

    當前作物基因發(fā)掘研究進展比較快,但這些基因的發(fā)掘多處在理論研究的層面。邱麗娟等[72]對在主要農(nóng)作物中克隆得到的300多個基因進行研究發(fā)現(xiàn),只有約18%的基因轉(zhuǎn)入來源作物并進行了功能檢測,表現(xiàn)出應用潛力,但在育種中進行利用價值評估的基因數(shù)目卻寥寥無幾[2]。小麥粒重性狀是由多基因控制的數(shù)量性狀,到目前為止還沒有粒重基因在小麥育種中應用情況的報道。在現(xiàn)代小麥育種過程中,單獨一個基因可能對提高粒重起到很微小的作用(甚至不起作用),況且產(chǎn)量性狀也受環(huán)境因素以及小麥植株本身的株型葉相等性狀的影響,對于我們闡述的已經(jīng)克隆到的與粒重相關的基因,其單獨在小麥中的轉(zhuǎn)入來增加小麥產(chǎn)量是行不通的,必須協(xié)調(diào)好基因、環(huán)境以及植株本身的優(yōu)良性狀才能較大幅度的提高小麥產(chǎn)量。另外,根據(jù)這些基因序列開發(fā)的特異性分子標記進行分子標記輔助選擇也將推動小麥產(chǎn)量的提高和高產(chǎn)育種的快速發(fā)展。

    [1] RAJARAM S.Role of conventional plant breeding and biotechnology in future wheat production [J].TurkishJournalofAgricultureandForestry,2005,29:105.

    [2] GIYRA A,SAULESCU N.Chromosomal location of genes controlling grain size in a large grained selection of wheat (TriticumaestivumL.) [J].Euphytica,1996,89:77.

    [3] GEGAS V C,NAZARI A,GRIFFITHS S,etal.A genetic framework for grain size and shape variation in wheat [J].ThePlantCell,2010,22:1046.

    [4] CRISTINA D,CIUCA M,CORNEA P C.Genetic contral of grain size and weight in wheat-where are we now? [J].ScientificBulletin.SeriesF.Biotechnologies,2016,2285:27.

    [5] DHOLAKIA B,AMMIRAJU J,SINGH H,etal.Molecular marker analysis of kernel size and shape in bread wheat [J].PlantBreed,2003,122:392.

    [6] BRESEGHELLO F,SORRELLS M E.QTL analysis of kernel size and shape in two hexaploid wheat mapping populations [J].FieldCropResearch,2007,101:172.

    [7] 馬 崗,徐兆飛,張哲夫,等.小麥粒形與產(chǎn)量關系的初步研究 [J].麥類作物學報,1997,17(5):18.

    MA G,XU Z F,ZHANG Z F,etal.The preliminary research on the relationship between wheat grain shape and yield [J].JournalofTriticeaeCrops,1997,17(5):18.

    [8] 王瑞霞,張秀英,伍 玲,等.不同生態(tài)環(huán)境下冬小麥籽粒大小相關性狀的QTL分析[J].中國農(nóng)業(yè)科學,2009,42(2):398.

    WANG R X,ZHANG X Y,WU L,etal.QTL analysis of grain size and related traits in winter wheat under different ecological environments [J].ScientiaAgricultureSinica,2009,42(2):398.

    [9] SUN X Y,WU K,ZHAO Y,etal.QTL analysis of kernel shape and weight using recombinant inbred lines in wheat [J].Euphytica,2009,165:615.

    [10] WANG S,ZHANG X,CHEN F,etal.A single nucleotide polymorphism of TaGS5 gene revealed its association with kernel weight in Chinese bread wheat [J].FrontiersinPlantScience,2015,6:1166.

    [11] MA L,LI T,HAO C,etal.TaGS5-3A,a grain size gene selected during wheat improvement for larger kernel and yield [J].PlantBiotechnologyJournal,2015,5(14):1269.

    [12] LING H Q,ZHAO S,LIU D,etal.Draft genome of the wheat A-genome progenitorTriticumurartu[J].Nature,2013,496:87.

    [13] DONG L,WANG F,LIU T,etal.Natural variation of TaGASR7-A1 affects grain length in common wheat under multiple cultivation conditions [J].MolecularBreeding,2014,34:937.

    [14] KUMAR A,MANTOVANI E E,SEETAN R,etal.Dissection of genetic factors underlying wheat kernel shape and size in an Elite×Non adapted cross using hight density SNP linkage map [J].ThePlantGenome,2016,9:1.

    [15] HOU J,JIANG Q,HAO C,etal.Global selection on sucrose synthase haplotypes during a century of wheat breeding [J].PlantPhysiology,2014,164:1918.

    [16] CHANG J,ZHANG J,MAO X,etal.Polymorphism of TaSAP1-A1 and its association with agronomic traits in wheat [J].Planta,2013,237:1495.

    [17] MOHLER V,ALBRECHT T,CASTELL A,etal.Considering causal genes in the genetic dissection of kernel traits in common wheat [J].JournalofAppliedGenetics,2016,57:467.

    [18] PATIL R M,TAMHANKAR S A,OAK M D,etal.Mapping of QTL for agronomic traits and kernel characters in durum wheat (TriticumdurumDesf.) [J].Euphytica,2013,190(1):117.

    [19] YANG D L,JING R L,CHANG X P,etal.Identification of quantitative trait loci and environmental interactions for accumulation and remobilization of water-soluble carbohydrates in wheat (TriticumaestivumL.) stems [J].Genetics,2007,176:571.

    [20] HUYNH B L,WALLWORK H,PALMER L,etal.Quantitative trait loci for grain fructan concentration in wheat (TriticumaestivumL.) [J].TheoreticalandAppliedGenetics,2008,117:701.

    [21] REBETZKE G J,HERWAARDEN A F,JENKINS C,etal.Quantitative trait loci for soluble stem carbohydrate production in wheat [J].AustralianJournalofAgriculturalResearch,2008,59:891.

    [22] ZHANG B,LI W,CHANG X,etal.Effects of favorable alleles for water-soluble carbohydrates at grain filling on grain weight under drought and heat stresses in wheat [J].PLoSOne,2014,9(7):e102917.

    [23] FU X Y,QI Z L,LI S S,etal.QTL detection for water-soluble oligosaccharide content of grain in common wheat(TriticumaestivumL.) [J].MolecularPlantBreeding,2011,2(10):68.

    [24] LI W Y,ZHANG B,LI R Z,etal.Favorable alleles for stem water-soluble carbohydrates identified by association analysis contribute to grain weight under drought stress conditions in wheat [J].PLoSOne,2015,10(3):e0119438.

    [25] KAWAKAMI A,YOSHIDA M.Molecular characterization of sucrose:sucrose1-fructosyltrans ferase and sucrose:fructan 6-fructosy ltrans ferase associated with fructan accumulation in winter wheat during cold hardening [J].Bioscience,Biotechnology,andBiochemistry,2002,66(11):2297.

    [26] YUE A Q,LI A,MAO X G,etal.Identification and development of a functional marker from 6-SFT-A2 associated with grain weight in wheat [J].MolecularBreeding,2015,35(2):63.

    [27] SU Z Q,HAO C,WANG L,etal.Identification and development of a functional marker of TaGW2 associated with grain weight in bread wheat (TriticumaestivumL.) [J].TheoreticalandAppliedGenetics,2011,122:211.

    [28] HUANG X Q,KEMPF H,GANAL M W,etal.Advanced back-cross QTL analysis in progenies derived from a cross between a German elite winter wheat variety and a synthetic wheat (TriticumaestivumL.) [J].TheoreticalandAppliedGenetics,2004,109:933.

    [29] HUANG X Q,CLOUTIER S,LYCAR L,etal.Molecular detection of QTL for agronomic and quality traits in a doubled haploid population derived from two Canadian wheats (TriticumaestivumL.) [J].TheoreticalandAppliedGenetics,2006,113:753.

    [30] SUN X C,MARZA F,MA H X,etal.Mapping quantitative trait loci for quality factors in an inter-class of USA and Chinese wheat [J].TheoreticalandAppliedGenetics,2010,120:1041.

    [31] JAISWAL V,GAHLAUT V,MATHUR S,etal.Identification of novel SNP in promoter sequence of TaGW2-6A associated with grain weight and other agronomic traits in wheat (TriticumaestivumL.) [J].PLoSOne,2015,10:e0129400.

    [32] QIN L,HAO C,HOU J,etal.Homologous haplotypes,expression,genetic effects and geographic distribution of thewheat yield gene TaGW2 [J].PlantBiology,2014,14:107.

    [33] HONG Y,CHEN L,DU L,etal.Transcript suppression of TaGW2 increased grain width and weight in bread wheat [J].FunctionalIntegrativeGenomics,2014,14:341.

    [34] 張 磊,張寶石,周榮華,等.小麥細胞分裂素氧化/脫氫酶基因(TaCKX5)的克隆及其染色體定位[J].中國農(nóng)業(yè)科學,2008,41(3):636.

    ZHANG L,ZHANG B S,ZHOU R H,etal.Isolation and chromosomal localization of cytokine in oxidase/dehydrogenase gene (TaCKX5) in wheat [J].ScientiaAgricultureSinica,2008,41(3):636.

    [35] ZHANG L,ZHAO Y L,GAO L F,etal. TaCKX6-D1,the ortholog of rice OsCKX2,is associated with grain weight in hexaploid wheat [J].NewPhytologist,2012,195:574.

    [36] CHANG C,LU J,ZHANG H P,etal.Copy number variation of cytokinin oxidase gene Tackx4 associated with grain weight and chlorophyll content of flag leaf in common wheat [J].PLoSOne,2015,10:e0145970.

    [37] LU J,CHANG C,ZHANG H P,etal.Identification of a novel allele of TaCKX6a02 associated with grain size,filling rate and weight of common wheat [J].PLoSOne,2015,10:e0144765.

    [38] HANIF M,GAO F,LIU J,etal. TaTGW6-A1,anortholog of rice TGW6,is associated with grain weight and yield in bread wheat [J].MolecularBreeding,2016,36:1.

    [39] MA D Y,YAN J,HE Z H,etal.Characterization of a cell wall invertase gene TaCwiA1 on common wheat chromosome 2A and development of functional markers [J].MolecularBreeding,2012,29:43.

    [40] 韓利明,楊芳萍,夏先春,等. 株高、粒重及抗病相關基因在不同國家小麥品種中的分布 [J].麥類作物學報,2011,31(5):824.

    HAN L M,YANG F P,XIA X C,etal.Distribution of genes related to plant height,gernel weight and disease resistance among wheat cultivars from major countries [J].JournalofTriticeaeCrops,2011,31(5):824.

    [41] JIANG Y,JIANG Q,HAO C,etal.A yield-associated gene TaCWI in wheat:its function,selection and evolution in global breeding revealed by haplotype analysis [J].TheoreticalandAppliedGenetics,2015,128:131.

    [42] BERMARD S M,MOLLER A L B,DIONISIO G,etal.Gene expression,cellular localisation and function of glutamine synthetase isozymes in wheat (TriticumaestivumL.) [J].PlantMolecularBiology,2008,67,89.

    [43] YAMAYA T,OBARA M,NAKAJIMA H,etal.Genetic manipulation and quantitative-trait loci mapping for nitrogen recycling in rice [J].JournalofExperimentalBotany,2002,53:917.

    [44] OBARA M,SATO T,SASAKI S,etal.Identification and characterization of a QTL on chromosome 2 for cytosolic glutamine synthetase content and panicle number in rice [J].TheoreticalandAppliedGenetics,2004,110:1.

    [45] HIRAL B,BERTIN P,QUILLERE I,etal.Towards a better understanding of the genetic and physiological basis for nitrogen use efficiency in maize [J].PlantPhysiology,2001,125:1258.

    [46] HIREL B,GOUIS J,NEY B,etal.The challenge of improving nitrogen use efficiency in crop plants:towards a more central role for genetic variability and quantitative genetics within integrated approaches [J].JournalofExperimentalBotany,2007,58:2369.

    [47] GALLAIS A,HIREL B.An approach to the genetics of nitrogen use efficiency in maize [J].JournalofExperimentalBotany,2004,55:295.

    [48] HABASH D Z,BEMARD S,SCHONDELMAIER J,etal.The genetics of nitrogen use in hexaploid wheat:N utilisation,development and yield [J].TheoreticalandAppliedGenetics,2007,114:403.

    [49] FEDOROVA L,FEDOROV A.Introns in gene evolution [J].Genetica,2003,118:123.

    [50] GUO Y,SUN J J,ZHANG G Z,etal.Haplotype,molecular marker and phenotype effects associated with mineral nutrient and grain size traits of TaGS1a in wheat [J].FieldCropsResearch,2003,154:119.

    [51] ZHENG J,LIU H,WANG Y,etal. TEF-7A,a transcript elongation factor gene,influences yield-related traits in bread wheat (TriticumaestivumL.) [J].JournalofExperimentalBotany,2014,65(18):5351.

    [52] ZHANG Z G,Lü G D,LI B,etal.Isolation and characterization of the TaSnRK2.10 gene and its association with agronomic traits in wheat (TriticumaestivumL.) [J].PLoSOne,2017,12(3):e0174425.

    [53] 王 倩,毛新國,昌小平,等.小麥 TaSnRK2.10的多態(tài)性及與農(nóng)藝性狀的關聯(lián)[J].中國農(nóng)業(yè)科學,2014,47(10):1865.

    WANG Q,MAO X G,CHANG X P,etal.Polymorphism of TaSnRK2.10 and its association with yield-related traits in wheat [J].ScientiaAgricultureSinica,2014,47(10):1865.

    [54] VARSHNEY R K,GRANER A,SORRELLS M E.Genomics-assisted breeding for crop improvement [J].TrendsinPlantScience,2005,10:621.

    [55] TAKAJASHI J S,PINTO L H,VITATEMA M H.Forward and reverse genetic approaches to behavior in the mouse [J].Science,1994,264:1724.

    [56] COULSON A,SULSTON J,BRENNER S,etal.Toward a physical map of the genome of the nematodeCaenorhabditiselegans[J].ProceedingoftheNationalAcademyofSciencesUSA,1986,83:7821.

    [57] GUPTA P K,RUSTGI S,KULWAL P L.Linkage disequilibrium and association studies in higher plants:present status and future prospects [J].PlantMolecularBiology,2005,57:461.

    [58] MATZKE M,MATZKE A J M,KOOTER J M.RNA:guiding gene silencing [J].Science,2001,293:1080.

    [59] TIAN Z,WANG X,LEE R,etal.Artificial selection for determinate growth habit in soybean [J].ProceedingoftheNationalAcademyofSciencesUSA,2010,107:8563.

    [60] FAN L J,BAO J D,WANG Y,etal.Post domestication selection in the maize starch pathway [J].PLoSOne,2009,4:e7612.

    [61] YIN Z T,MENG F F,SONG H N,etal.Expression quantitative trait loci analysis of two genes encoding rubisco activase in soybean [J].PlantPhysiology,2010,152:1625.

    [62] HU K M,QIU D Y,SHEN X L,etal.Isolation and manipulation of quantitative trait loci for disease resistance in rice using a candidate gene approach [J].MolecularPlant,2008,1:786.

    [63] ZHU J.Mixed model approaches of mapping genes for complex quantitative traits [J].JournalofZhejiangUniversity,1999,33(3):327.

    [64] WANG D L,ZHU J,LI Z K,etal.Mapping QTLs with epistatic effects and QTL-environment interactions by mixed linear model approaches [J].TheoreticalandAppliedGenetics,1999,99:1255.

    [65] WU C C,SUN S K,NIMMAKAYALA P,etal.A BAC- and BIBAC-based physical map of the soybean genome [J].GenomeResearch,2004,14:319.

    [66] YANG J,ZHU J,WILLIAMS R W.Mapping the genetic architecture of complex traits in experimental populations [J].Bioinformatics,2007,23:1527.

    [67] YANG J,HU C,HU H,etal.QTL Network:mapping and visualizing genetic architecture of complex traits in experimental populations [J].Bioinformatics,2008,24:721.

    [68] WU W R.Genetic analysis of qualitative-quantitative trait in multiple-environment experiments [J].ActaAgronomicaSinica,1999,25(6):723.

    [69] LI W M,WU W R,LU H R.A method of detecting linkage between quantitative trait loci and genetic marker and its application in wheat [J].ActaAgronomicaSinica,1993,19(2):97.

    [70] AO Y,XU C W.Maximum likelihood method for mapping QTL in four-way cross design [J].ActaAgronomicaSinica,2006,32(1):51 .

    [71] WANG J K.Inclusive composite interval mapping of quantitative trait genes [J].ActaAgronomicaSinica,2009,35(2):239.

    [72] 邱麗娟,郭 勇,黎 裕,等.中國作物新基因發(fā)掘:現(xiàn)狀、挑戰(zhàn)與展望[J].作物學報,2011,37(1):1.

    QIU L J,GUO Y,LI Y,etal.Novel gene discovery of crops in China:status,challenging,and perspective [J].ActaAgronomicaSinica,2011,37(1):1.

    [73] BARKER D G,BIANCHI S,BLONDON F,etal.Medicagotruncatula,a model plant for studying the molecular genetics of theRhizobium-legume symbiosis [J].PlantMolecularBiologyReporter,1990,8:40.

    [74] HANDBERG K,STOUGAARD J.Lotusjaponicus,an autogamous,diploid legume species for classical and molecular genetics [J].ThePlantJournal,1992,2:487.

    [75] YU J,HU S N,WANG J,etal.A draft sequence of the rice genome (OryzasativaL.ssp.indica) [J].Science,2002,296:79.

    [76] YU J,WANG J,LIN W,etal.The genomes ofOryzasativa:a history of duplications [J].PLoSBiology,2005,3:e38.

    ResearchProgressofGeneCloningRelatedtoGrainWeightinCommonWheat

    LüGuangde,SUNXianyin,QIXiaolei,WANGChao,WANGRuixia,SUNYingying,MUQiuhuan,MIYong,QIANZhaoguo,WUKe

    (Tai’an Academy of Agricultural Science, Tai’an, Shandong 271000,China)

    With the rapid development of molecular biology, especially bioinformatics and genomics,gene cloning related to grain weight in common wheat has also been further studied. This paper mainly introduces the progress of genetic and functional researches on grain weight related genes in wheat, including TaGS5, TaGASR7, TaSus, TaWSC, Ta-6-SFT, TaGW2, TaCKX, TaTGW6, TaCWI, and TaGS1 and gene family of TaTEF, TaSAP, and TaSnRK2s, which provides the corresponding theory and technology support in molecular breeding of wheat.

    Common wheat; Genes related to grain weight; Cloning

    2017-01-11

    2017-06-18

    國家小麥現(xiàn)代產(chǎn)業(yè)技術體系項目(CARS-3-2-22);山東省現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術體系項目(SDAIT-04-021-12);山東省重點研發(fā)計劃項目(科技攻關部分)(2016GNC113004)

    E-mail:2007guangd@163.com

    吳 科(E-mail:sdtawuke1964@126.com)

    時間:2017-10-11

    網(wǎng)絡出版地址:http://kns.cnki.net/kcms/detail/61.1359.S.20171011.1601.008.html

    S512.1;S330

    A

    1009-1041(2017)10-1301-08

    猜你喜歡
    關聯(lián)產(chǎn)量
    2022年11月份我國鋅產(chǎn)量同比增長2.9% 鉛產(chǎn)量同比增長5.6%
    提高玉米產(chǎn)量 膜下滴灌有效
    不懼于新,不困于形——一道函數(shù)“關聯(lián)”題的剖析與拓展
    “苦”的關聯(lián)
    當代陜西(2021年17期)2021-11-06 03:21:36
    世界致密油產(chǎn)量發(fā)展趨勢
    海水稻產(chǎn)量測評平均產(chǎn)量逐年遞增
    2018年我國主要水果產(chǎn)量按省(區(qū)、市)分布
    “一帶一路”遞進,關聯(lián)民生更緊
    當代陜西(2019年15期)2019-09-02 01:52:00
    2018上半年我國PVC產(chǎn)量數(shù)據(jù)
    聚氯乙烯(2018年9期)2018-02-18 01:11:34
    奇趣搭配
    欧美成人午夜免费资源| 久久国产亚洲av麻豆专区| 国产午夜精品久久久久久一区二区三区| 日韩一本色道免费dvd| 成人毛片60女人毛片免费| 人成视频在线观看免费观看| 91久久精品电影网| 日韩一本色道免费dvd| 精品国产一区二区久久| 国产精品秋霞免费鲁丝片| 国产欧美另类精品又又久久亚洲欧美| 精品久久蜜臀av无| 亚洲色图 男人天堂 中文字幕 | 美女视频免费永久观看网站| 国产视频内射| 男人爽女人下面视频在线观看| 国产高清有码在线观看视频| 各种免费的搞黄视频| 亚洲精品乱久久久久久| 人人妻人人澡人人看| 制服诱惑二区| 3wmmmm亚洲av在线观看| 大话2 男鬼变身卡| 永久免费av网站大全| 天堂中文最新版在线下载| 满18在线观看网站| 久久久午夜欧美精品| 伊人亚洲综合成人网| 一级,二级,三级黄色视频| 亚洲av成人精品一二三区| 免费不卡的大黄色大毛片视频在线观看| 97超视频在线观看视频| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久伊人网av| 国产精品人妻久久久影院| 午夜激情久久久久久久| 日韩中文字幕视频在线看片| 中文字幕最新亚洲高清| 下体分泌物呈黄色| 国产精品一二三区在线看| 亚洲精品一区蜜桃| 91精品一卡2卡3卡4卡| 午夜福利,免费看| 久久女婷五月综合色啪小说| 国产成人精品福利久久| 丰满迷人的少妇在线观看| 人人妻人人澡人人看| 搡老乐熟女国产| 精品久久久精品久久久| 国产精品久久久久久久久免| 午夜影院在线不卡| 国产又色又爽无遮挡免| 亚洲人成77777在线视频| 汤姆久久久久久久影院中文字幕| 国产精品人妻久久久影院| 久久青草综合色| videossex国产| 欧美日韩av久久| 纯流量卡能插随身wifi吗| 秋霞伦理黄片| 精品人妻熟女毛片av久久网站| 91精品三级在线观看| 一级毛片黄色毛片免费观看视频| 中文天堂在线官网| 最新的欧美精品一区二区| 少妇被粗大猛烈的视频| 99热国产这里只有精品6| 又粗又硬又长又爽又黄的视频| av国产久精品久网站免费入址| 下体分泌物呈黄色| 亚洲国产欧美日韩在线播放| 欧美97在线视频| 国产淫语在线视频| 久久女婷五月综合色啪小说| 中文天堂在线官网| 久久精品国产自在天天线| 搡女人真爽免费视频火全软件| 国产男人的电影天堂91| 狂野欧美激情性xxxx在线观看| 亚洲高清免费不卡视频| 亚洲欧美清纯卡通| 国产一级毛片在线| 免费少妇av软件| 欧美+日韩+精品| 免费高清在线观看日韩| 日韩一本色道免费dvd| 一级毛片 在线播放| 亚洲欧洲精品一区二区精品久久久 | 好男人视频免费观看在线| 美女国产高潮福利片在线看| 另类精品久久| 乱人伦中国视频| 一级毛片aaaaaa免费看小| 九草在线视频观看| 波野结衣二区三区在线| 亚洲av日韩在线播放| 超色免费av| 极品人妻少妇av视频| av免费在线看不卡| 亚洲av福利一区| 久久午夜综合久久蜜桃| 最后的刺客免费高清国语| av国产久精品久网站免费入址| 久久精品夜色国产| 99久久人妻综合| 亚洲性久久影院| 成人亚洲欧美一区二区av| 日韩av在线免费看完整版不卡| 成人黄色视频免费在线看| 亚洲综合色惰| 日韩av免费高清视频| 久久97久久精品| 亚洲精品久久午夜乱码| 日本欧美国产在线视频| 国产 精品1| 久久精品国产鲁丝片午夜精品| 久热这里只有精品99| 日韩在线高清观看一区二区三区| 亚洲人成网站在线观看播放| 在线看a的网站| 超色免费av| 多毛熟女@视频| 亚洲精品国产av成人精品| 蜜桃国产av成人99| 在线观看免费视频网站a站| 我的老师免费观看完整版| 精品一区在线观看国产| 成年av动漫网址| 色网站视频免费| 亚州av有码| 熟女人妻精品中文字幕| 国产成人精品福利久久| 高清av免费在线| 欧美日韩av久久| 国产色婷婷99| 高清视频免费观看一区二区| 国产精品一区二区三区四区免费观看| 欧美精品一区二区大全| 最黄视频免费看| 日本午夜av视频| 色哟哟·www| 美女cb高潮喷水在线观看| 三上悠亚av全集在线观看| 两个人免费观看高清视频| 国产极品天堂在线| 午夜91福利影院| 夜夜爽夜夜爽视频| 69精品国产乱码久久久| 下体分泌物呈黄色| 黄色配什么色好看| 一边亲一边摸免费视频| 韩国高清视频一区二区三区| 国产成人精品在线电影| 国产亚洲最大av| 秋霞在线观看毛片| 一本一本综合久久| 午夜av观看不卡| 极品人妻少妇av视频| 国产视频首页在线观看| videossex国产| 免费黄网站久久成人精品| 久久av网站| 汤姆久久久久久久影院中文字幕| 免费av中文字幕在线| 一本一本综合久久| 九草在线视频观看| 精品午夜福利在线看| 国产日韩欧美视频二区| 国产av码专区亚洲av| 精品久久久久久久久av| 久久精品国产a三级三级三级| 日本黄大片高清| 亚洲情色 制服丝袜| 国产乱人偷精品视频| 成人毛片a级毛片在线播放| 9色porny在线观看| 新久久久久国产一级毛片| 国产女主播在线喷水免费视频网站| 狂野欧美激情性xxxx在线观看| 性色av一级| 精品亚洲成a人片在线观看| 亚洲精华国产精华液的使用体验| 一个人看视频在线观看www免费| 国产男女超爽视频在线观看| 亚洲精品日韩av片在线观看| 亚洲欧洲日产国产| 插阴视频在线观看视频| 久久久久久久国产电影| 天天影视国产精品| 三级国产精品欧美在线观看| 自线自在国产av| 少妇精品久久久久久久| 亚洲欧美成人精品一区二区| 成年美女黄网站色视频大全免费 | 我要看黄色一级片免费的| 国产精品久久久久久精品古装| 成人综合一区亚洲| 欧美xxⅹ黑人| 人人妻人人澡人人看| 毛片一级片免费看久久久久| 亚洲精品一区蜜桃| 成人毛片60女人毛片免费| 国产高清国产精品国产三级| 99热这里只有是精品在线观看| 99国产精品免费福利视频| 免费播放大片免费观看视频在线观看| 一区二区三区乱码不卡18| 大又大粗又爽又黄少妇毛片口| 大陆偷拍与自拍| 男女免费视频国产| 精品久久久久久久久av| 最近2019中文字幕mv第一页| 18+在线观看网站| 亚洲精品视频女| 亚洲情色 制服丝袜| 美女视频免费永久观看网站| 色5月婷婷丁香| 成人黄色视频免费在线看| 久久精品久久精品一区二区三区| 天堂俺去俺来也www色官网| 黄色一级大片看看| 亚洲综合色惰| 国产精品嫩草影院av在线观看| 精品久久久噜噜| 日韩一本色道免费dvd| 夫妻午夜视频| 欧美xxxx性猛交bbbb| 久久久精品免费免费高清| av有码第一页| 欧美亚洲 丝袜 人妻 在线| 大片电影免费在线观看免费| 免费观看无遮挡的男女| 久久99热6这里只有精品| 午夜视频国产福利| 美女脱内裤让男人舔精品视频| 18禁裸乳无遮挡动漫免费视频| 色婷婷久久久亚洲欧美| 蜜桃久久精品国产亚洲av| 日韩熟女老妇一区二区性免费视频| 国产成人免费观看mmmm| 人人妻人人添人人爽欧美一区卜| 日本欧美国产在线视频| 91成人精品电影| 中文字幕人妻丝袜制服| 97在线视频观看| 晚上一个人看的免费电影| 亚洲中文av在线| 亚洲怡红院男人天堂| 美女福利国产在线| 午夜福利视频在线观看免费| 美女脱内裤让男人舔精品视频| 国产成人免费观看mmmm| 成人无遮挡网站| 在线观看人妻少妇| 久久久亚洲精品成人影院| 黄色怎么调成土黄色| 成年av动漫网址| 成人免费观看视频高清| 久久久亚洲精品成人影院| 欧美少妇被猛烈插入视频| 国产av国产精品国产| 免费看av在线观看网站| 成年人免费黄色播放视频| 精品国产乱码久久久久久小说| 美女福利国产在线| 欧美成人午夜免费资源| 全区人妻精品视频| 这个男人来自地球电影免费观看 | 国产一区二区在线观看日韩| 午夜久久久在线观看| 成人二区视频| 国产成人精品福利久久| 观看av在线不卡| 毛片一级片免费看久久久久| 纯流量卡能插随身wifi吗| 免费大片黄手机在线观看| 草草在线视频免费看| 2021少妇久久久久久久久久久| 国产色爽女视频免费观看| 寂寞人妻少妇视频99o| 少妇的逼好多水| 国产精品99久久99久久久不卡 | 五月开心婷婷网| 日本wwww免费看| 少妇猛男粗大的猛烈进出视频| 亚洲精品中文字幕在线视频| 少妇被粗大的猛进出69影院 | h视频一区二区三区| 最近中文字幕2019免费版| 18在线观看网站| 综合色丁香网| 亚洲欧美日韩另类电影网站| 国产成人午夜福利电影在线观看| 汤姆久久久久久久影院中文字幕| 在线观看美女被高潮喷水网站| 国产精品久久久久久精品古装| 国产精品一区二区三区四区免费观看| 男人操女人黄网站| 黄片播放在线免费| 韩国高清视频一区二区三区| 伦精品一区二区三区| 国产黄频视频在线观看| 9色porny在线观看| 如何舔出高潮| 欧美最新免费一区二区三区| av又黄又爽大尺度在线免费看| 下体分泌物呈黄色| 高清视频免费观看一区二区| 日本vs欧美在线观看视频| 亚洲精品日本国产第一区| 99精国产麻豆久久婷婷| av国产久精品久网站免费入址| 看免费成人av毛片| videosex国产| 午夜福利在线观看免费完整高清在| 久久久久精品久久久久真实原创| 精品少妇内射三级| 亚洲av成人精品一区久久| freevideosex欧美| 精品少妇内射三级| kizo精华| 亚洲av男天堂| 日韩 亚洲 欧美在线| 成年美女黄网站色视频大全免费 | 欧美+日韩+精品| 国产伦精品一区二区三区视频9| 午夜91福利影院| 涩涩av久久男人的天堂| 久久久久国产网址| 视频区图区小说| 亚洲国产色片| 中国美白少妇内射xxxbb| 精品一区二区免费观看| 最近的中文字幕免费完整| 色网站视频免费| 99九九线精品视频在线观看视频| 又黄又爽又刺激的免费视频.| 成人国产麻豆网| 中文字幕制服av| 亚洲色图 男人天堂 中文字幕 | 国产成人a∨麻豆精品| 蜜桃在线观看..| 777米奇影视久久| 高清在线视频一区二区三区| 交换朋友夫妻互换小说| 免费大片18禁| 午夜影院在线不卡| 观看av在线不卡| 99久国产av精品国产电影| 日韩成人伦理影院| 亚洲精品色激情综合| 国产69精品久久久久777片| 国产av精品麻豆| 成人国语在线视频| 久久久精品区二区三区| 黑丝袜美女国产一区| 久久久久人妻精品一区果冻| 热re99久久精品国产66热6| 精品酒店卫生间| 国产精品三级大全| 熟妇人妻不卡中文字幕| 91午夜精品亚洲一区二区三区| 卡戴珊不雅视频在线播放| av在线app专区| 在线观看免费日韩欧美大片 | 人人妻人人爽人人添夜夜欢视频| www.色视频.com| 在线免费观看不下载黄p国产| 伦精品一区二区三区| 久久久久久伊人网av| 日本黄色片子视频| 久久久久久久久久成人| 又大又黄又爽视频免费| 亚洲国产最新在线播放| 久久精品国产亚洲av天美| 欧美精品一区二区大全| 国产视频首页在线观看| 国产日韩欧美亚洲二区| 2022亚洲国产成人精品| 99九九在线精品视频| 考比视频在线观看| 精品久久久久久电影网| 韩国av在线不卡| 两个人的视频大全免费| 久久午夜福利片| 伦精品一区二区三区| 汤姆久久久久久久影院中文字幕| av网站免费在线观看视频| 日本与韩国留学比较| 亚洲五月色婷婷综合| 母亲3免费完整高清在线观看 | 美女国产高潮福利片在线看| 免费av不卡在线播放| 中文字幕av电影在线播放| 国产精品不卡视频一区二区| 一级毛片 在线播放| 青春草国产在线视频| 永久网站在线| 国语对白做爰xxxⅹ性视频网站| 人妻制服诱惑在线中文字幕| 国产亚洲最大av| 国产成人aa在线观看| 久久女婷五月综合色啪小说| 亚洲精品国产色婷婷电影| 国产片特级美女逼逼视频| 少妇高潮的动态图| 午夜激情久久久久久久| 久久久久久人妻| 蜜臀久久99精品久久宅男| 欧美 亚洲 国产 日韩一| a级片在线免费高清观看视频| 激情五月婷婷亚洲| av专区在线播放| 爱豆传媒免费全集在线观看| 晚上一个人看的免费电影| 午夜福利网站1000一区二区三区| 精品国产国语对白av| 伦精品一区二区三区| 国产亚洲最大av| 亚洲国产毛片av蜜桃av| 亚洲怡红院男人天堂| 黑人巨大精品欧美一区二区蜜桃 | 亚洲一区二区三区欧美精品| 国内精品宾馆在线| 大香蕉久久网| 18禁在线播放成人免费| 日韩精品免费视频一区二区三区 | 国产极品天堂在线| 欧美日韩视频高清一区二区三区二| 亚洲精品国产av蜜桃| 免费人妻精品一区二区三区视频| 一区二区三区乱码不卡18| 国产精品一区二区三区四区免费观看| 日日摸夜夜添夜夜爱| 赤兔流量卡办理| 日日撸夜夜添| 在线观看www视频免费| 国产探花极品一区二区| 免费大片18禁| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 热re99久久精品国产66热6| 国产成人精品无人区| 国产黄频视频在线观看| 制服丝袜香蕉在线| 亚洲国产av新网站| 国产精品免费大片| 日韩免费高清中文字幕av| 黄片播放在线免费| 成人手机av| 久久影院123| 国产精品久久久久久精品古装| 一个人看视频在线观看www免费| 成人手机av| av在线播放精品| 亚洲美女搞黄在线观看| 久久精品人人爽人人爽视色| 最黄视频免费看| 国模一区二区三区四区视频| 三上悠亚av全集在线观看| 色94色欧美一区二区| 国产伦精品一区二区三区视频9| 日韩亚洲欧美综合| 亚洲欧美成人精品一区二区| 国产一区二区在线观看日韩| 能在线免费看毛片的网站| 精品一品国产午夜福利视频| 精品久久久久久电影网| 日本黄色片子视频| 黄色配什么色好看| 亚洲av男天堂| 欧美日韩亚洲高清精品| 日本免费在线观看一区| 欧美少妇被猛烈插入视频| 中文字幕制服av| 免费人妻精品一区二区三区视频| 高清不卡的av网站| av在线老鸭窝| 蜜臀久久99精品久久宅男| 国产成人免费无遮挡视频| 乱码一卡2卡4卡精品| 2022亚洲国产成人精品| 在线观看免费高清a一片| 秋霞伦理黄片| 观看av在线不卡| 女性被躁到高潮视频| 久久这里有精品视频免费| 国产伦精品一区二区三区视频9| av女优亚洲男人天堂| 五月开心婷婷网| 久久午夜综合久久蜜桃| 亚洲av在线观看美女高潮| 亚洲av日韩在线播放| 99九九线精品视频在线观看视频| 午夜影院在线不卡| 一区二区av电影网| 国产精品久久久久久精品古装| 国产精品一区二区在线观看99| 欧美另类一区| 亚洲精品久久成人aⅴ小说 | 欧美bdsm另类| 曰老女人黄片| 色94色欧美一区二区| 一区二区三区四区激情视频| 亚洲精品一二三| 亚洲精品色激情综合| 韩国高清视频一区二区三区| 久久女婷五月综合色啪小说| 国产精品欧美亚洲77777| 久久久久久久大尺度免费视频| 欧美日韩综合久久久久久| 免费久久久久久久精品成人欧美视频 | 国产黄色视频一区二区在线观看| 免费看不卡的av| 99热国产这里只有精品6| av不卡在线播放| 日韩成人伦理影院| 日日撸夜夜添| 如何舔出高潮| 亚洲精品日韩在线中文字幕| 热99国产精品久久久久久7| 免费大片黄手机在线观看| 大片电影免费在线观看免费| 18禁动态无遮挡网站| 日韩人妻高清精品专区| 国产免费一区二区三区四区乱码| 人妻一区二区av| 高清av免费在线| 丝瓜视频免费看黄片| 制服诱惑二区| 亚洲五月色婷婷综合| 大话2 男鬼变身卡| 在线观看美女被高潮喷水网站| 黄色欧美视频在线观看| 最新的欧美精品一区二区| 99久久精品国产国产毛片| 亚洲无线观看免费| 国产欧美另类精品又又久久亚洲欧美| 妹子高潮喷水视频| 午夜日本视频在线| 人妻少妇偷人精品九色| 亚洲经典国产精华液单| 国产高清国产精品国产三级| 精品人妻一区二区三区麻豆| 亚洲av中文av极速乱| 久久人人爽av亚洲精品天堂| videos熟女内射| 国产一区二区三区av在线| av电影中文网址| 久久av网站| 欧美精品国产亚洲| 蜜桃久久精品国产亚洲av| 91午夜精品亚洲一区二区三区| 人妻人人澡人人爽人人| 伦理电影大哥的女人| 亚洲欧洲精品一区二区精品久久久 | 亚洲欧洲精品一区二区精品久久久 | 国产熟女午夜一区二区三区 | 国产在线视频一区二区| 成人漫画全彩无遮挡| 免费少妇av软件| 成人国产av品久久久| 久久久精品免费免费高清| av福利片在线| 九色亚洲精品在线播放| 波野结衣二区三区在线| 国产日韩一区二区三区精品不卡 | 午夜福利,免费看| 日本爱情动作片www.在线观看| 麻豆成人av视频| 18禁在线播放成人免费| 午夜免费男女啪啪视频观看| xxx大片免费视频| 国产在线一区二区三区精| 五月伊人婷婷丁香| 青春草国产在线视频| 美女大奶头黄色视频| 最后的刺客免费高清国语| 性高湖久久久久久久久免费观看| 午夜久久久在线观看| 久久精品国产a三级三级三级| 26uuu在线亚洲综合色| 亚洲不卡免费看| 少妇人妻久久综合中文| 成年av动漫网址| 亚洲精品乱久久久久久| 免费人成在线观看视频色| 99热网站在线观看| av在线app专区| 久久久国产一区二区| 国产精品不卡视频一区二区| 日韩成人av中文字幕在线观看| 免费看光身美女| 精品久久久精品久久久| 中文天堂在线官网| 另类精品久久| 精品一区二区免费观看| 在线观看一区二区三区激情| √禁漫天堂资源中文www| 在线播放无遮挡| 观看美女的网站| 免费大片黄手机在线观看| 国产精品久久久久久av不卡| 亚洲经典国产精华液单| 久久午夜福利片| 久久精品人人爽人人爽视色| 少妇猛男粗大的猛烈进出视频| a级片在线免费高清观看视频| 亚洲内射少妇av| 新久久久久国产一级毛片| 精品国产乱码久久久久久小说| 啦啦啦啦在线视频资源| 这个男人来自地球电影免费观看 | 国产亚洲午夜精品一区二区久久| 午夜影院在线不卡| 国产一区亚洲一区在线观看|