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

    分蘗相關(guān)基因在再生稻腋芽萌發(fā)期的表達(dá)分析

    2021-02-09 00:32:08朱容劉歆吳建偉劉念梁大成劉章勇
    關(guān)鍵詞:表達(dá)分析腋芽再生稻

    朱容 劉歆 吳建偉 劉念 梁大成 劉章勇

    摘要:【目的】明確分蘗相關(guān)基因在再生稻腋芽萌發(fā)期的表達(dá)動(dòng)態(tài)變化,為闡明水稻再生季腋芽發(fā)育的調(diào)控機(jī)理提供理論依據(jù)?!痉椒ā恳栽偕痔Y能力強(qiáng)的材料i21和再生分蘗能力弱的材料i89為研究對(duì)象,運(yùn)用實(shí)時(shí)熒光定量PCR(qRT-PCR)檢測(cè)8個(gè)分蘗相關(guān)基因在水稻頭季(分蘗始期和拔節(jié)孕穗期)和再生季不同時(shí)期(頭季去穗后1、24、48和72 h)的表達(dá)動(dòng)態(tài)變化,并分析分蘗相關(guān)基因間及其與最大再生分蘗的相關(guān)性,篩選出與腋芽生長(zhǎng)相關(guān)性高的基因?!窘Y(jié)果】自頭季稻稻穗收割后,材料i21和i89的再生分蘗數(shù)的變化趨勢(shì)基本一致,均隨著腋芽的生長(zhǎng)再生分蘗數(shù)逐漸增多,收割后20 d達(dá)最多并趨于穩(wěn)定,但材料i21再生季發(fā)苗多且快,再生分蘗發(fā)生速度明顯高于材料i89;材料i21和i89的最大再生分蘗數(shù)為33和10個(gè)。LAX1、LAX2和MOC1基因的相對(duì)表達(dá)量在不同材料不同時(shí)期間均存在明顯差異,其中自頭季去穗后,除MOC1基因在頭季去穗后72 h的相對(duì)表達(dá)量外,LAX1、LAX2和MOC1基因在水稻再生分蘗形成期均表現(xiàn)為材料i21中的相對(duì)表達(dá)量極顯著高于材料i89(P<0.01,下同),MOC1基因在頭季去穗后1 h達(dá)峰值,LAX1和LAX2基因在頭季去穗后24 h達(dá)峰值,隨后顯著降低(P<0.05,下同)并趨于穩(wěn)定。頭季去穗后,除TAD1基因在頭季去穗后24 h時(shí)外,材料i89中D10、HTD1、OsTB1和TAD1基因的相對(duì)表達(dá)量均極顯著高于材料i21中的相對(duì)表達(dá)量,D10和HTD1基因在頭季去穗后1 h達(dá)峰值,OsTB1和TAD1基因在頭季去穗后24和48 h達(dá)峰值,D27基因在材料i21和i89中的相對(duì)表達(dá)量均先升高后降低,且頭季去穗后24~72 h,D27基因在材料i21中的表達(dá)量極顯著高于材料i89。最大再生分蘗數(shù)分別與LAX2和MOC1基因的表達(dá)水平呈極顯著和顯著正相關(guān),與D10、HTD1、OsTB1和TAD1基因的表達(dá)量呈極顯著負(fù)相關(guān)。此外,不同基因的表達(dá)水平也存在不同相關(guān)性?!窘Y(jié)論】LAX1、LAX2和MOC1基因高表達(dá)有利于水稻再生季分蘗的形成,D10、HTD1、OsTB1和TAD1基因高表達(dá)會(huì)抑制水稻再生季分蘗的形成。不同分蘗相關(guān)基因間存在相互作用,共同調(diào)控水稻再生季分蘗的形成。

    關(guān)鍵詞: 分蘗基因;再生稻;腋芽;表達(dá)分析

    中圖分類號(hào):S511.035.3? ?; ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼: A 文章編號(hào):2095-1191(2021)10-2632-09

    Abstract:【Objective】 The expression of tiller-related genes during axillary bud outgrowth in ratoon rice were profiled to provide insight into the molecular mechanisms underlying axillary bud development in the ratoon season. 【Method】 The rice varieties i21 and i89 with high and low ratooning abilities, respectively, were used as the study materials. Real-time fluorescence quantitatire PCR (qRFRCR) was used to characterize the expression dynamics of eight tiller-related genes at the stages of tiller initiation, elongation and booting of the first season and at 1, 24, 48 and 72 h after heading removal in the ratoon season. The correlation between tiller-related genes and maximum ratoon tillers was analyzed and the genes with high correlation with axillary bud growth were selected.【Result】The number of ratoon tillers in i21 and i89 increased gradually with the growth of axillary buds and reached a stable maximum 20 d after the first-season harvest. However, i21 ratoon tillers grew much faster than those of i89 and the maximum number of ratoon tillers of i21 and i89 were 33 and 10, respectively. The relative expression levels of LAX1, LAX2 and MOC1 genes differed significantly in the two varieties. The relative expression levels of LAX1, LAX2 and MOC1 genes in i21 were significantly higher than those in i89 during the ratooning period (P<0.01, the same below). MOC1 reached peak levels at 1 h after beheading, whereas the expressions of LAX1 and LAX2? peaked at 24 h after beheading then decreased significantly (P<0.05, the same below) to stable levels thereafter. During ratooning, the relative expression levels of D10, HTD1, OsTB1 and TAD1 genes in i89 were significantly higher than those in i21, except for TAD1 at 24 h after beheading. D10 and HTD1 genes peaked at 1 h after beheading, while OsTB1 and TAD1 genes reached their peak values later at 24 and 48 h after beheading. The relative expression levels of D27 in i21 and i89 increased initially but decreased at 24-72 h after heading removal, although this decrease in i21 was significantly less than that in i89. The maximum tiller number was positively and significantly correlated with the expression levels of LAX2 and MOC1 and negatively correlated with those of D10, HTD1, OsTB1 and TAD1. In addition, different gene expression levels also had different correlations.【Conclusion】The high expression of LAX1, LAX2 and MOC1 genes was beneficial to the formation of ratoon tillers, while the high expression of D10, HTD1, OsTB1 and TAD1 genes inhibited ratoon till outgrowth. This work indicates that the different tiller-related genes interacted to jointly regulate ratoon tiller formation.

    Key words:tiller gene; ratoon rice; axillary bud; expression analysis

    Foundation item:Key Scientific and Technological Project of Hubei (2020BBA044)

    0 引言

    【研究意義】水稻(Oryza sativa L.)是我國(guó)乃至全世界重要的糧食作物,其產(chǎn)量一直備受關(guān)注(劉貴富等,2018)。目前增加水稻播種面積和單產(chǎn)日益困難,增加復(fù)種指數(shù)成為提高稻谷產(chǎn)量的關(guān)鍵舉措(Ray Foley,2013;Dong et al.,2017)。再生稻作為一種土地高效利用、資源節(jié)約型的栽培模式,具有米質(zhì)優(yōu)、低污染、低能耗的特點(diǎn),是提高復(fù)種指數(shù)、增加稻谷產(chǎn)量的有效途徑(朱永川等,2013;徐富賢,2015;Jiang et al.,2016;Tao et al.,2016;劉歆等,2020)。由于分蘗力是水稻的重要農(nóng)藝性狀,直接影響稻谷產(chǎn)量(李學(xué)勇等,2003;李經(jīng)勇等,2009;孫秀紅和胡波,2020),因此,腋芽萌發(fā)能力是再生稻產(chǎn)量的關(guān)鍵因素(李貴勇等,2015;習(xí)敏等,2017;胡香玉等,2019;汪浩等,2019)。對(duì)再生稻分蘗力進(jìn)行分子生物學(xué)方面的研究,以期解析再生稻分蘗產(chǎn)生和調(diào)控的分子機(jī)理,對(duì)選育強(qiáng)再生力的水稻品種具有重要的指導(dǎo)意義。【前人研究進(jìn)展】目前,已研究證實(shí)水稻分蘗相關(guān)基因MONOCULM1(MOC1)、LAX1、LAX2、D10、HTD1、D27、OsTB1和TAD1的表達(dá)影響水稻分蘗形成,其中MOC1基因是第一個(gè)發(fā)現(xiàn)的控制水稻分蘗的單蘗基因,其等位基因?yàn)镚NP6(Li et al.,2003;Shao et al.,2019;Zhang et al.,2021),編碼植物特有的GRAS家族蛋白(Bolle,2004),在腋生分生組織的發(fā)生、分蘗芽的形成過(guò)程中發(fā)揮促進(jìn)作用;LAX1基因在植株的側(cè)生分生組織的起始發(fā)育發(fā)揮重要作用,編碼非典型的bHLH(basic-Helix-Loop-Helix)型轉(zhuǎn)錄因子,在頂端分生組織與新形成的分生組織之間的邊界區(qū)域表達(dá),主要調(diào)節(jié)水稻腋芽原基的形成(Maekawa et al.,2001;Komatsu et al.,2003;Oikawa and Kyozuka,2009);LAX2基因編碼一個(gè)核蛋白,該蛋白含有植物特有的保守結(jié)構(gòu)域,能與LAX1發(fā)生蛋白互作,Gnp4為其等位基因,LAX2/Gnp4基因編碼一個(gè)功能未知的核蛋白,調(diào)控水稻腋生分生組織的形成(Tabuchi et al.,2011;Zhang et al.,2011);D10、HTD1和D27基因參與獨(dú)腳金內(nèi)酯(Strigolactones,SLs)的生物合成,SLs是一種控制植物分枝的新激素,其可抑制植物分枝和側(cè)芽生長(zhǎng)(Gomez-Roldan et al.,2008;Umehara et al.,2008);D10基因與擬南芥的MAX4、豌豆的RMS1和矮牽牛的DAD1是直系同源基因,編碼胡蘿卜素裂解雙加氧酶OsCCD8,該酶抑制側(cè)芽的向外生長(zhǎng)從而使水稻分蘗數(shù)減少(Arite et al.,2007);HTD1基因編碼胡蘿卜素裂解雙加氧酶OsCCD7,該酶抑制水稻分枝的生長(zhǎng)從而減少水稻分蘗數(shù)(Zou et al.,2005,2006);D27基因編碼定位于葉綠體上的α-β類胡蘿卜素異構(gòu)酶含鐵蛋白,主要在莖頂端和根部表達(dá),在SLs生物合成中發(fā)揮關(guān)鍵作用,使水稻分蘗數(shù)減少(Lin et al.,2009;Adrian et al.,2012);OsTB1/FC1位于同一基因位點(diǎn),作用于SLs的下游,編碼含TCP結(jié)構(gòu)域的bHLH轉(zhuǎn)錄因子,對(duì)分蘗發(fā)揮負(fù)調(diào)控作用;OsTB1基因?qū)Ψ痔Y數(shù)的負(fù)調(diào)控作用是通過(guò)控制腋芽的生長(zhǎng)過(guò)程來(lái)實(shí)現(xiàn),腋芽中OsTB1基因過(guò)量表達(dá)會(huì)阻礙其形成,最終導(dǎo)致水稻分蘗數(shù)顯著減少,而OsTB1基因功能喪失型突變體fc1分蘗數(shù)顯著增加(Takeda et al.,2003;Minakuchi et al.,2010);TAD1基因編碼細(xì)胞分裂后期啟動(dòng)復(fù)合物(Anaphase-promoting complex,APC/C)的共激活子,與MOC1基因作用于同一途徑,通過(guò)細(xì)胞周期依賴的方式降低靶標(biāo)MOC1基因的表達(dá)水平,抑制腋生分生組織的起始和形成,從而抑制其分蘗(Lin et al.,2012;Xu et al.,2012)。以上研究表明MOC1、LAX1和LAX2基因促進(jìn)水稻分蘗的形成,而D10、HTD1、D27、OsTB1和TAD1基因抑制水稻分蘗的發(fā)生?!颈狙芯壳腥朦c(diǎn)】雖然分蘗相關(guān)基因的功能及其對(duì)水稻分蘗形成影響已有大量研究報(bào)道,但目前鮮見(jiàn)有關(guān)分蘗相關(guān)基因在水稻再生季的表達(dá)特性的研究報(bào)道?!緮M解決的關(guān)鍵問(wèn)題】以本課題組篩選出的再生力極端表型水稻品種為材料,利用實(shí)時(shí)熒光定量PCR(qRT-PCR)檢測(cè)水稻分蘗相關(guān)基因在極端表型材料中再生季腋芽發(fā)生不同時(shí)期的表達(dá)差異,以期為水稻再生季穗數(shù)形成的分子機(jī)理提供理論依據(jù)。

    1 材料與方法

    1. 1 試驗(yàn)材料

    192份核心水稻種質(zhì)材料,由中國(guó)農(nóng)業(yè)科學(xué)院作物科學(xué)研究所提供,用于篩選再生能力極端表型的材料。植物總RNA提取試劑盒E.Z.N.A.? Plant RNA Kit購(gòu)自上海北諾科技生物科技有限公司;2×EasyTaq PCR SuperMix(AS111-01)購(gòu)自北京全式金生物技術(shù)有限公司;cDNA反轉(zhuǎn)錄試劑盒Fast Quant RT Kit(with gDNA)和熒光定量試劑盒SuperReal PreMix Plus(SYBR Green)購(gòu)自天根生化科技(北京)有限公司;100 bp DNA Ladder、1 kb DNA Ladder、4S Green Plus染料和6×Loading Buffer均購(gòu)自生工生物工程(上海)股份有限公司;Tris堿、硼酸和EDTA購(gòu)自北京索萊寶科技有限公司。5×TBE貯存液是由54 g/L Tris堿、27.5 g/L硼酸、20 mL 0.5 mol/L EDTA(pH 8.0)配制而成,將其稀釋10倍成0.5×TBE;50×TAE貯存液用242 g/L Tris堿、37.2 g/L EDTA和57.1 mL醋酸加水定容至1 L配制。主要設(shè)備儀器:超速離心機(jī)(Eppendorf,德國(guó))、PCR儀(Bio-Rad,美國(guó))、實(shí)時(shí)熒光定量PCR儀(Bio-Rad CFX96TM Real-Time System,美國(guó))、凝膠成像系統(tǒng)(Bio-Rad,美國(guó))、Nanodrop ND-1000分光光度計(jì)(美國(guó))和水平瓊脂糖凝膠電泳儀(北京君意東方電泳設(shè)備有限公司)。

    1. 2 試驗(yàn)方法

    1. 2. 1 材料篩選 192份核心水稻種質(zhì)材料于2018年分別在荊州、海南陵水和廣東英德共3個(gè)生態(tài)點(diǎn)種植,調(diào)查頭季和再生季最大分蘗數(shù),篩選出再生能力極端表型的材料,編號(hào)為i21和i89,其中,i21為再生分蘗能力強(qiáng)的材料,i89為再生分蘗能力弱的材料。

    1. 2. 2 田間種植 篩選出的材料于2019年夏季在長(zhǎng)江大學(xué)農(nóng)學(xué)院試驗(yàn)基地采用盆栽方式種植,塑料盆上直徑36 cm、高31 cm,每盆裝18±0.3 kg干土,土表距離盆口約4 cm。3月31日播種,4月25日移栽,每盆種4穴,單穴單苗,每個(gè)品種3盆,3次重復(fù),正常水肥管理。

    1. 2. 3 再生季分蘗動(dòng)態(tài)調(diào)查 自頭季稻去穗后,每隔5 d調(diào)查5株水稻,掛牌記錄再生分蘗數(shù),直至達(dá)到最大分蘗為止。

    1. 2. 4 取樣方法 設(shè)6次取樣時(shí)期,分別為頭季分蘗始期(CK)、頭季拔節(jié)孕穗期(R1)、頭季去穗后1 h (R2)、頭季去穗后24 h(R3)、頭季去穗后48 h(R4)和頭季去穗后72 h(R5),取來(lái)自5個(gè)單株節(jié)間倒2節(jié)腋芽的混合樣品(CK為5個(gè)單株水稻基部2 cm左右的混合樣品,統(tǒng)一在水稻達(dá)到4個(gè)分蘗時(shí)取樣),取樣后立即經(jīng)液氮冷凍30 min并于-80 ℃冰箱保存。每個(gè)時(shí)期設(shè)3個(gè)生物學(xué)重復(fù)。

    1. 2. 5 總RNA提取及反轉(zhuǎn)錄合成cDNA 利用E.Z.N.A.? Plant RNA Kit植物總RNA提取試劑盒提取水稻樣品總RNA,用1%瓊脂糖凝膠電泳檢測(cè)總RNA的完整性,取1 μL用Nanodrop ND-1000分光光度計(jì)測(cè)定其純度和濃度,結(jié)果顯示RNA提取效果好,且RNA無(wú)降解;OD260/OD280為1.8~2.0,純度較好,可滿足后續(xù)的試驗(yàn)要求。

    利用Fast Quant RT Kit(with gDNA)試劑盒反轉(zhuǎn)錄合成cDNA。以Actin為內(nèi)參基因,采用普通PCR檢測(cè)各材料cDNA質(zhì)量,結(jié)果證實(shí)cDNA無(wú)gDNA污染,質(zhì)量較好,可用于后續(xù)試驗(yàn)。

    1. 2. 6 引物設(shè)計(jì) 通過(guò)國(guó)家水稻數(shù)據(jù)中心(http://www.ricedata.cn/)和已報(bào)道的分蘗基因的相關(guān)文獻(xiàn)資料(Komatsu et al.,2003;Li et al.,2003;Takeda et al.,2003;Arite et al.,2007;Hong et al.,2009),查找分蘗相關(guān)基因信息,然后通過(guò)NCBI數(shù)據(jù)庫(kù)(https://www.ncbi.nlm.nih.gov/)獲取分蘗相關(guān)基因LAX1、LAX2、MOC1、D10、HTD1、D27、OsTB1和TAD1基因序列,最后通過(guò)BLAST比對(duì)確定其編碼蛋白的保守區(qū)域,然后利用Primer Premier 5.0設(shè)計(jì)長(zhǎng)度合適、GC含量合理、無(wú)發(fā)卡結(jié)構(gòu)、無(wú)錯(cuò)配及無(wú)雜帶的引物(表1)。

    1. 2. 7 實(shí)時(shí)熒光定量PCR檢測(cè) 以Actin作為內(nèi)參基因,采用實(shí)時(shí)熒光定量PCR檢測(cè)8個(gè)水稻分蘗相關(guān)基因的表達(dá)情況。按照SuperReal PreMix Plus說(shuō)明配制反應(yīng)體系20.0 μL:10 μL 2×SuperReal PreMix Plus,cDNA模板2.0 μL,10 μmol/L上、下游引物各0.6 μL,ddH20補(bǔ)足至20.0 μL。反應(yīng)結(jié)束后,通過(guò)檢測(cè)熔解曲線確定反應(yīng)產(chǎn)物特異性,采用比較CT值法(2-△△Ct)計(jì)算基因的相對(duì)表達(dá)量(Gomez-Roldan et al.,2008)。

    1. 3 統(tǒng)計(jì)分析

    使用Excel 2007進(jìn)行數(shù)據(jù)整理,用SPSS 21.0對(duì)數(shù)據(jù)進(jìn)行顯著性分析。

    2 結(jié)果與分析

    2. 1 再生分蘗特性比較結(jié)果

    由圖1-A可知,自頭季稻稻穗收割后,材料i21和i89再生分蘗數(shù)的變化趨勢(shì)基本一致,均隨著腋芽的生長(zhǎng)再生分蘗數(shù)逐漸增多,收割后20 d達(dá)最多并趨于穩(wěn)定,但材料i21再生季發(fā)苗多且快,再生分蘗發(fā)生速度明顯高于材料i89。由圖1-B可知,這2個(gè)材料的最大再生分蘗數(shù)存在極顯著差異(P<0.01,下同),其中材料i21的最大再生分蘗數(shù)為33個(gè);材料i89的最大再生分蘗數(shù)為10個(gè),且不同年限間2個(gè)材料的再生能力表現(xiàn)穩(wěn)定(圖1)。

    2. 2 分蘗正調(diào)控基因在不同材料不同時(shí)期間的表達(dá)分析結(jié)果

    由圖2可知,自頭季去穗后,LAX1基因在材料i21和i89的相對(duì)表達(dá)量顯著升高(P<0.05,下同),至頭季去穗后24 h達(dá)峰值,隨后顯著降低并趨于穩(wěn)定;LAX2基因在材料i21的相對(duì)表達(dá)量變化趨勢(shì)與LAX1基因基本一致,均呈先升高再降低的變化趨勢(shì),在頭季去穗后24 h達(dá)峰值,但LAX2基因在材料i89中的表達(dá)量自頭季去穗后呈降低─升高─降低的變化趨勢(shì),最后逐漸趨于穩(wěn)定;MOC1基因在材料i21和i89的相對(duì)表達(dá)量的變化趨勢(shì)與LAX2基因基本一致。

    由圖2還可知,自頭季去穗后,除MOC1基因在頭季去穗后72 h外,LAX1、LAX2和MOC1基因在水稻再生分蘗形成期(R2~R5)均表現(xiàn)為再生分蘗能力強(qiáng)的材料i21中的相對(duì)表達(dá)量極顯著高于分蘗能力弱的材料i89??梢?jiàn),LAX1、LAX2和MOC1基因高表達(dá)可促進(jìn)水稻再生季分蘗形成,表明這3個(gè)基因是參與水稻分蘗芽形成的關(guān)鍵基因,對(duì)分蘗芽的形成發(fā)揮促進(jìn)作用,可提高水稻分蘗數(shù)。

    2. 3 分蘗負(fù)調(diào)控基因在不同材料不同時(shí)期的表達(dá)分析結(jié)果

    2. 3. 1 SLs合成基因的表達(dá)差異 由圖3可知,自頭季去穗后,D10和HTD1基因在材料i21中的相對(duì)表達(dá)量顯著降低并逐漸趨于穩(wěn)定;D10基因在材料i89的相對(duì)表達(dá)量顯著升高,頭季去穗后1 h達(dá)峰值,隨后顯著降低并趨于穩(wěn)定;材料i89中HTD1基因的相對(duì)表達(dá)量在去頭季穗后1 h顯著升高并達(dá)峰值,隨后先降低再升高并逐漸趨于穩(wěn)定,且其在材料i89的相對(duì)表達(dá)量均高于材料i21;D27基因在材料i21和i89中的相對(duì)表達(dá)量先升高后降低,且頭季去穗24 h后其在材料i21中的相對(duì)表達(dá)量極顯著高于材料i89,說(shuō)明D10和HTD1基因高表達(dá)會(huì)抑制水稻再生季分蘗形成。

    2. 3. 2 OsTB1和TAD1基因的表達(dá)差異 由圖4可知,自頭季去穗后,OsTB1基因在材料i21中的相對(duì)表達(dá)量出現(xiàn)2次先升高后降低,頭季去穗后48 h達(dá)最大值,OsTB1基因在材料i89中的相對(duì)表達(dá)量先升高后降低,頭季去穗后24 h達(dá)峰值,隨后顯著降低,且OsTB1基因在材料i21中的相對(duì)表達(dá)量均顯著低于材料i89;TAD1基因在材料i21和i89中的相對(duì)表達(dá)量均呈先升高后降低趨勢(shì),但在材料i21中的相對(duì)表達(dá)量在頭季去穗后24 h達(dá)峰值,在材料i89中的相對(duì)表達(dá)量在頭季去穗后48 h達(dá)峰值,且其在材料i89中的相對(duì)表達(dá)量顯著高于材料i21,說(shuō)明OsTB1和TAD1基因高表達(dá)可抑制水稻再生季分蘗形成。

    2. 4 分蘗相關(guān)基因間及其與最大再生分蘗數(shù)間的相關(guān)分析結(jié)果

    由表2可知,最大再生分蘗數(shù)分別與LAX2和MOC1基因的表達(dá)水平呈極顯著和顯著正相關(guān),與TAD1、OsTB1、D10和HTD1基因的表達(dá)水平呈極顯著負(fù)相關(guān)。此外,不同基因的表達(dá)水平也存在不同相關(guān)性,D27基因與LAX1、LAX2和MOC1基因呈極顯著正相關(guān);LAX1基因與LAX2基因、LAX2基因與MOC1基因、D10基因與HTD1基因均呈極顯著正相關(guān);LAX2基因與HTD1基因、D27基因與OsTB1基因均呈極顯著負(fù)相關(guān);OsTB1基因與MOC1和LAX2基因均呈顯著負(fù)相關(guān);D27基因與TAD1基因也呈顯著負(fù)相關(guān)??梢?jiàn),分蘗相關(guān)基因的表達(dá)水平與水稻再生季分蘗形成密切相關(guān),而且不同分蘗相關(guān)基因也存在相互作用,共同調(diào)控水稻再生季分蘗的形成。同時(shí),進(jìn)一步證實(shí)LAX1、LAX2和MOC1基因的表達(dá)水平升高有利于水稻再生季分蘗形成,D10、HTD1、OsTB1和TAD1基因的表達(dá)水平升高會(huì)抑制水稻再生季分蘗形成。

    3 討論

    水稻分蘗是一個(gè)受遺傳基因和環(huán)境雙重影響的復(fù)雜數(shù)量性狀(譚震波等,1997;楊川航等,2012;李興星等,2016)。已有研究表明,LAX1、LAX2和MOC1基因均能促進(jìn)水稻腋芽的生長(zhǎng)萌發(fā),正向調(diào)控水稻分蘗的形成(Maekawa et al.,2001;Li et al.,2003;Tabuchi et al.,2011)。本研究發(fā)現(xiàn),除MOC1基因在頭季去穗后72 h的相對(duì)表達(dá)量外,LAX1、LAX2和MOC1基因在水稻再生分蘗形成期(R2~R5)均表現(xiàn)為再生分蘗能力強(qiáng)的材料i21中的相對(duì)表達(dá)量極顯著高于分蘗能力弱的材料i89,說(shuō)明LAX1、LAX2和MOC1基因?qū)τ谒驹偕痔Y的形成發(fā)揮促進(jìn)作用,與其在水稻頭季的研究結(jié)果一致,均能正向調(diào)節(jié)水稻分蘗的形成。MOC1基因在頭季去穗后72 h的相對(duì)表達(dá)量表現(xiàn)為再生分蘗能力強(qiáng)的材料i21顯著低于分蘗能力弱的材料i89,與頭季去穗后其他時(shí)期(R2~R4)剛好相反,其原因可能是頭季去穗后72 h水稻再生腋芽已萌發(fā),處于伸長(zhǎng)生長(zhǎng)階段,MOC1基因的相對(duì)表達(dá)量升高,與其他通路的基因存在關(guān)聯(lián),抑制了分蘗,具體的機(jī)制尚不清楚,有待進(jìn)一步研究探索。

    D10、HTD1、D27、OsTB1和TAD1基因是從一系列水稻突變體材料定位克隆出的分蘗相關(guān)基因,均可抑制水稻側(cè)芽的生長(zhǎng),從而負(fù)調(diào)節(jié)水稻分蘗的形成(Takeda et al.,2003;Zou et al.,2005;Arite et al.,2007;Lin et al.,2009;Xu et al.,2012)。本研究中,D10、HTD1、OsTB1和TAD1基因均在再生分蘗能力弱的材料i89中相對(duì)表達(dá)量較高,說(shuō)明這些基因可抑制水稻再生分蘗的形成,與其在水稻頭季的研究結(jié)果一致,均能負(fù)向調(diào)節(jié)水稻分蘗的形成。D27基因在頭季拔節(jié)孕穗期和頭季去穗后1 h的相對(duì)表達(dá)量呈現(xiàn)出i89高于i21,在頭季去穗后24~72 h的相對(duì)表達(dá)量為i21高于i89,與前人研究結(jié)果不一致,導(dǎo)致無(wú)法判斷其對(duì)再生分蘗的形成是否存在影響,其原因可能是D27基因在無(wú)效分蘗的控制中起作用。造成再生季中基因表達(dá)差異的原因很多,其一,再生季的生長(zhǎng)過(guò)程本身是一個(gè)復(fù)雜過(guò)程,其生長(zhǎng)環(huán)境不同,基因表達(dá)調(diào)控也會(huì)存在很大差異;其二,材料本身的遺傳背景差異很大,造成最終分蘗結(jié)果存在差異;其三,頭季去穗后會(huì)產(chǎn)生基因表達(dá)量無(wú)變化或產(chǎn)生與基因功能不一致的結(jié)果,這可能是植株本身對(duì)外界變化的適應(yīng)能力體現(xiàn)(Cline and Oh,2006)或是再生季的調(diào)控過(guò)程本身的不同所導(dǎo)致,推測(cè)水稻再生季通過(guò)其他途徑調(diào)控再生苗的生長(zhǎng)。

    水稻分蘗相關(guān)基因形成一個(gè)復(fù)雜的調(diào)控網(wǎng)絡(luò),且存在組織特異性、生長(zhǎng)時(shí)期特異性及特異性表達(dá)。Tabuchi等(2011)通過(guò)對(duì)稀穗突變體lax1與lax2的研究發(fā)現(xiàn),LAX2是通過(guò)與LAX1相互作用從而調(diào)控水稻分蘗的形成;Xu等(2012)研究表明TAD1基因可通過(guò)細(xì)胞周期依賴的方式降低MOC1基因的表達(dá)水平從而抑制水稻分蘗的發(fā)生。本研究相關(guān)分析結(jié)果顯示,LAX1基因與LAX2基因的呈極顯著正相關(guān)(相關(guān)系數(shù)為0.780**),MOC1基因與TAD1基因呈顯著負(fù)相關(guān)(相關(guān)系數(shù)為-0.490*),與前人的研究結(jié)果(Tabuchi et al.,2011;Xu et al.,2012)一致。此外,本研究還發(fā)現(xiàn),不同分蘗相關(guān)基因間存在復(fù)雜的相關(guān)關(guān)系,說(shuō)明水稻再生分蘗的形成可能與水稻頭季分蘗一樣,存在著復(fù)雜的分子調(diào)控機(jī)制,且不同基因之間存在相互作用,共同調(diào)控水稻再生季分蘗的形成,但其具體的分子遺傳機(jī)理還有待進(jìn)一步研究。

    本研究分析了分蘗相關(guān)基因在再生稻腋芽萌發(fā)期的表達(dá)差異,后續(xù)的研究工作可從以下兩方面開(kāi)展:一方面可利用本課題已有的遺傳多樣性豐富的材料群體進(jìn)行重測(cè)序,鑒定其基因型,并結(jié)合再生稻分蘗能力的表型進(jìn)行全基因組關(guān)聯(lián)分析,獲得顯著性SNP位點(diǎn);另一方面,構(gòu)建F2群體,對(duì)其進(jìn)行重測(cè)序,鑒定基因型,結(jié)合表型進(jìn)行QTL分析,從而鑒定再生季分蘗發(fā)生的關(guān)鍵調(diào)控基因。借助分子手段深入研究水稻再生季分蘗的基因調(diào)控網(wǎng)絡(luò),對(duì)于揭示再生稻的遺傳機(jī)理具有重要意義。

    4 結(jié)論

    LAX1、LAX2和MOC1基因高表達(dá)有利于水稻再生季分蘗的形成,D10、HTD1、OsTB1和TAD1基因高表達(dá)會(huì)抑制水稻再生季分蘗的形成。不同分蘗相關(guān)基因間存在相互作用,共同調(diào)控水稻再生季分蘗的形成。

    參考文獻(xiàn):

    胡香玉,鐘旭華,彭碧琳,劉彥卓,黃農(nóng)榮,梁開(kāi)明,潘俊峰,傅友強(qiáng). 2019. 不同氮肥運(yùn)籌下低樁機(jī)收再生稻的產(chǎn)量和經(jīng)濟(jì)效益[J]. 中國(guó)稻米,(4):1-7. doi:10.3969/j.issn.1006-8082. 2019.04.004. [Hu X Y,Zhong X H,Peng B L,Liu Y Z,Huang N R,Liang K M,Pan J F,F(xiàn)u Y Q. 2019. Grain yield and profit of machine-harvested low stubble ratoon rice under different nitrogen management[J]. China Rice,(4):1-7.]

    李貴勇,譚麗明,劉玉文,黃潔,陳雙云,夏瓊梅,龍瑞平,楊從黨. 2015. 再生稻的產(chǎn)量形成及其高產(chǎn)群體分析[J]. 江西農(nóng)業(yè)大學(xué)學(xué)報(bào),37(4):577-583. [Li G Y,Tan L M,Liu Y W,Huang J,Chen S Y,Xia Q M,Long R P,Yang C D. 2015. Analysis on yield formation and its higher yielding population quality of ratooning rice[J]. Acta Agriculturae Universitatis Jiangxiensis,37(4):577-583.] doi: 10. 13836/j.jjau. 2015089.

    李經(jīng)勇,張洪松,唐永群. 2009. 中國(guó)再生稻研究與應(yīng)用[J]. 南方農(nóng)業(yè),3(3):88-92. [Li J Y,Zhang H S,Tang Y Q. 2009. Research and application of regenerated rice in China[J]. South China Agriculture,3(3):88-92.] doi:10. 19415/j.cnki.1673-890x.2009.03.030.

    李興星,鄭劍,周軍杰,楊永清,秦小健,張漢馬,梁永書(shū). 2016. 利用秈粳交“熱粳35/協(xié)B”F2群體對(duì)水稻再生力的QTL分析[J]. 分子植物育種,14(9):2383-2391. [Li X X,Zheng J,Zhou J J,Yang Y Q,Qin X J,Zhang H M,Liang Y S. 2016. QTL analysis of ratooning ability by using Indica-japonica “Rejing 35/Xie B” F2 populations in rice[J]. Molecular Plant Breeding,14(9):2383-2391.] doi:10.13271/j.mpb.014.002383.

    李學(xué)勇,錢前,李家洋. 2003. 水稻分蘗的分子機(jī)理研究[J]. 中國(guó)科學(xué)院院刊,18(4):274-276. [Li X Y,Qian Q,Li J Y. 2003. Progress in elucidating the molecular mechanism of rice tillering[J]. Bulletin of the Chinese Academy of Sciences,18(4):274-276.] doi:10.16418/j.issn.1000-3045. 2003.04.009.

    劉貴富,陳明江,李明,呂慧穎,葛毅強(qiáng),魏珣,楊維才. 2018.水稻育種行業(yè)創(chuàng)新進(jìn)展[J]. 植物遺傳資源學(xué)報(bào),19(3):416-429.[Liu G F,Chen M J,Li M,Lü H Y,Ge Y Q,Wei X,Yang W C. 2018. Advances and innovation of rice breeding[J]. Journal of Plant Genetic Resources,19(3):416-429.] doi:10.13430/j.cnki.jpgr.2018.03.006.

    劉歆,朱容,楊梅,劉章勇. 2020. 再生稻高產(chǎn)種質(zhì)資源篩選及高產(chǎn)群體分析[J]. 作物雜志,(2):28-33. [Liu X,Zhu R,Yang M,Liu Z Y. 2020. Screening of high-yield germplasms for ratoon rice and analysis of high yield composition[J]. Crops,(2):28-33.] doi:10.16035/j.issn.1001-7283.2020.02.005.

    孫秀紅,胡波. 2020. 影響水稻分蘗發(fā)生的主要栽培措施和環(huán)境因素研究進(jìn)展[J]. 現(xiàn)代農(nóng)業(yè)科學(xué),(7):6-8. [Sun X H,Hu B. 2020. Research progress on main cultivation measures and environmental factors affecting tillering of rice[J]. Modern Agricultural Science,(7):6-8.] doi:10.3969/ j.issn.1007-5739.2020.07.002.

    譚震波,沈利爽,袁祚廉,陸朝福,陳英,周開(kāi)達(dá),朱立煌. 1997. 水稻再生能力和頭季稻產(chǎn)量性狀的QTL定位及其遺傳效應(yīng)分析[J]. 作物學(xué)報(bào),23(3):289-295. [Tan Z B,Shen L S,Yuan Z L,Lu C F,Chen Y,Zhou K D,Zhu L H. 1997. Identification of QTLs for ratooning ability and grain yield traits of rice and analysis of their genetic effects[J]. Acta Agronomica Sinica,23(3):289-295.] doi:10.3321/j.issn:0496-3490.1997.03.005.

    汪浩,張強(qiáng),張文地,李思宇,黃健,朱安,劉立軍. 2020. 腋芽萌發(fā)能力對(duì)再生稻產(chǎn)量影響的研究進(jìn)展[J]. 中國(guó)水稻科學(xué),34(3):205-216. [Wang H,Zhang Q,Zhang W D,Li S Y,Huang J ,Zhu A,Liu L J. 2020. Advances in the effects of the ability of axillary bud germination on grain yield in ratoon rice[J]. Chinese Journal Rice Science,34(3):205-216.] doi:10.16819/j.1001-7216.2020.9097.

    習(xí)敏,吳文革,汪靖桂,汪宏偉,陳剛,許有尊. 2017. 再生稻種植產(chǎn)量差形成的研究[J]. 華北農(nóng)學(xué)報(bào),32(1):104-110. [Xi M,Wu W G,Wang J G,Wang H W,Chen G,Xu Y Z. 2017. Study on formation of grain yield differences in ratooning rice cultivation[J]. Acta Agriculturae Boreali-Sinica,32(1):104-110.] doi:10.7668/hbnxb.2017.01.017.

    徐富賢,熊洪,張林,朱永川,蔣鵬,郭曉藝,劉茂. 2015. 再生稻產(chǎn)量形成特點(diǎn)與關(guān)鍵調(diào)控技術(shù)研究進(jìn)展[J]. 中國(guó)農(nóng)業(yè)科學(xué),48(9):1702-1717. [Xu F X,Xiong H,Zhang L,Zhu Y C,Jiang P,Guo X Y,Liu M. 2015. Progress in research of yield formation of ratooning rice andIts high-yielding key regulation technologies[J]. Scientia Agricultura Sinica,48(9):1702-1717.] doi:10.3864/j.issn.0578-1752.2015.09.04.

    楊川航,王玉平,涂斌,李婷,胡亮,李仕貴. 2012. 利用RIL群體對(duì)水稻再生力及相關(guān)農(nóng)藝性狀的QTL分析[J]. 作物學(xué)報(bào),38(7):1240-1246. [Yang C H,Wang Y P,Tu B,Li T,Hu L,Li S G. 2012. QTL analysis of rice ratooning ability and related agronomic traits by using RIL populations[J]. Acta Agronomica Sinica,38(7):1240-1246.] doi:10.3724/SP.J.1006.2012.01240.

    朱永川,熊洪,徐富賢,郭曉藝,張林,劉茂,周興兵. 2013. 再生稻栽培技術(shù)的研究進(jìn)展[J]. 中國(guó)農(nóng)學(xué)通報(bào),29(36):1-8. [Zhu Y C,Xiong H,Xu F X,Guo X Y,Zhang L,Liu M,Zhou X B. 2013. Progress on research of ratoon rice cultivation technology[J].Chinese Agricultural Science Bulletin,29(36):1-8.] doi:10.11924/j.issn.1000-6850. 2013-0363.

    Adrian A,Muhammad J,Mattia M,Mark B,Martina V,Peter B,Sandro G,Harro B,Peter B,Salim A B. 2012. The path from β-carotene to carlactone,a strigolactone-like plant hormone[J]. Science,335(6074):1348-1351. doi:10.1126/ science.1218094.

    Arite T,lwata H,Ohshima K,Maekawa M,NakajimaM,Kojima M,Sakakibara H,Kyozuka J. 2007.DWARF10,an RMS1/MAX4/DAD1 ortholog,controls lateral bud outgrowth in rice[J]. The Plant Journal,51(6):1019-1029. doi:10.1111/j.1365-313X.2007.03210.x.

    Bolle C. 2004. The role of GRAS proteins in plant signal transduction and development[J]. Planta,2(18):683-692. doi:10.1007/s00425-004-1203-z.

    Cline M G,Oh C. 2006. A reappraisal of the role of abscisic acid and its interaction with auxin in apical dominance[J]. Annals of Botany,98(4):891-897. doi:10.1093/aob/mcl173.

    Dong H L,Chen Q,Wang W Q,Peng S B,Huang J L,Cui K H,Nie L X. 2017. The growth and yield of a wet-seeded rice-ratoon rice system in central China[J]. Field Crops Research,208:55-59. doi:10.1016/j.fcr.2017.04.003.

    Gomez-Roldan V,F(xiàn)ermas S,Brewer P B,Puech-Pages V,Dun E A,Pillot J P,Letisse F,Matusova R,Danoun S,Portais J C,Bouwmeester H,Becard G,Beveridge C A,Rameau C,Rochange S F. 2008. Strigolactone inhibition of shoot branching[J]. Nature,455(7210):189-194. doi:10.1038/nature07271.

    Hong C Y,Chao Y Y,Yang M Y,Cheng S Y,Cho S C,Kao C H. 2009. NaCl-induced expression of glutathione reductase in roots of rice(Oryza sativa L.) seedlings is mediated through hydrogen peroxide but not abscisic acid[J]. Plant and Soil,320:103-115. doi:10.1007/s11104-008-9874-z.

    Jiang Q W,Wang W Q,Chen Q,Peng S B,Huang J L,Cui K H,Nie L X. 2016. Response of first flood irrigation timing after rice dry-direct-seeding:Productivity and greenhouse gas emissions in Central China[J]. Agricultural Water Management,177:241-247. doi:10.1016/j.agwat.2016. 08. 006.

    Komatsu K,Maekawa M,Ujiie S,Satake Y,F(xiàn)urutani I,Okamoto H,Shimamoto K,Kyozuka J. 2003. LAX and SPA:Major regulators of shoot branching in rice[J]. Procee-dings of the National Academy of Sciences of the United States of America,2003,100(20):11765-11770. doi:10. 1073/pnas.1932414100.

    Li X Y,Qian Q,F(xiàn)u Z M,Wang Y H,Xiong G S,Zeng D L,Wang X Q,Liu X F,Teng S,Yuan M,Luo D,Han B,Li J Y. 2003. Control of tillering in rice[J]. Nature,422(6932):618-621. doi:10.1038/nature01518.

    Lin H,Wang R X,Qian Q,Yan M X,Meng X B,F(xiàn)u Z M,Yan C Y,Jiang B,Su Z,Li J Y,Wang Y H. 2009. DWARF27,an iron-containing protein required for the biosynthesisof strigolactones,regulates rice tiller bud outgrowth[J]. The Plant Cell,21(5):1512-1525. doi:10.11 05/tpc.109.065987.

    Lin Q B,Wang D,Dong H,Gu S H,Cheng Z J,Gong J,Qin R Z,Jiang L,Li G,Wang J L,Wu F Q,Guo X P,Zhang X,Lei C L,WangHY,WanJM. 2012. Rice APC/CTE controls tillering by mediating the degradation of MONOCULM 1[J]. Nature Communications,(3):644-656. doi:10.1038/ncomms1716.

    Maekawa M,Komatsu M,Shimamoto K,Kyozuka J. 2001. The LAX1 and FRIZZY PANICLE 2 genes determine the inflorescence architecture of rice by controlling rachis-branch and spikeletdevelopment[J]. Developmental Biology,231(2):364-373. doi:10.1006/dbio.2000. 9988.

    Minakuchi K,Kameoka H,Yasuno N,Umehara M,Luo L,Kobayashi K,Hanada A,Ueno K,Asami T,Yamaguchi S,Kyozuka J. 2010. FINE CULM1(FC1) works downstream of strigolactones to inhibit the outgrowth of axillary buds in rice[J]. Plant and Cell Physiology,51(7):1127-1135. doi:10. 1093/pcp/pcq083.

    Oikawa T,Kyozuka J. 2009. Two-step regulation of LAX PANICLE1 protein accumulation in axillary meristem formation in rice[J]. The Plant Cell,21(4):1095-1108. doi:10.1105/tpc.108.065425.

    Ray D K,F(xiàn)oley J A. 2013. Increasing global crop harvest frequency:Recent trends and future directions[J]. Environmental Research Letters,8:44041-44050. doi:10.1088/1748-9326/8/4 /044041.

    Shao G N,Lu Z F,Xiong J S,Wang B,Jing Y H,Meng X B,Liu G F,Ma H Y,Yan L,Yu H. 2019. Tiller bud formation regulators MOC1 and MOC3 cooperatively promote tiller bud outgrowth by activating FON1 expression in rice[J].Molecular Plant,12(8):1090-1102. doi:10.1016/j.molp. 2019.04.008.

    Tabuchi H,Zhang Y,Hattori S,Omae M,Shimizu-Sato S,Oikawa T,Qian Q,Nishimura M,Kitano H,Xie H,F(xiàn)ang X H,Yoshida H,Kyozuka J,Chen F,Sato Y. 2011. LAX PANICLE2 of rice encodes a novel nuclear protein and regulates the formation of axillary meristems[J].The Plant Cell,23(9):3276-3287. doi:10.1105/tpc.111.088 765.

    Takeda T,Suwa Y,Suzuki M,Kitano H,Ueguchi-Tanaka M,Ashikari M,Matsuoka M,Ueguchi C. 2003. The OsTB1 gene negatively regulates lateral branching in rice[J].The Plant Journal,33(3):513-520. doi:10.1046/j.1365-313X. 2003.01648.x.

    Tao Y,Chen Q,Peng S B,Wang W Q,Nie L X. 2016. Lower global warming potential and higher yield of wet direct-seeded rice in Central China[J].Agronomy for Sustaina-ble Development,36(2):24. doi:10.1007/s13593-016-0361-2.

    Umehara M,Hanada A,Yoshida S,Akiyama K,Arite T,Takeda-Kamiya N,Magome H,Kamiya Y,Shirasu K,Yoneyama K,Kyozuka J,Yamaguchi S. 2008. Inhibition of shoot branching by new terpenoid plant hormones[J]. Nature,455(7210):195-200. doi:10.1038/nature07272.

    Xu C,Wang Y H,Yu Y C,Duan J,Duan J B,Liao Z G,Xiong G S,Meng X B,Liu G F,Qian Q,Li J Y. 2012. Degradation of MONOCULM1 by APC/C(TAD1) regulates rice tillering[J]. Nature Communications,(3)743-750. doi:10.1038/ncomms1743.

    Zhang Z Y,Li J J,Yao G X,Zhang H L,Dou H J,Shi H L,Sun X M,Li Z C. 2011. Fine mapping and cloning of the grain number per-panicle gene(Gnp4) on chromosome 4 in rice(Oryza sativa L.)[J]. Agricultural Scien-ces in China,10(12):1825-1833. doi:10.1016/S1671- 2927 (11)60182-X.

    Zhang Z Y,Sun X M,Ma X Q,Xu B X,Zhao Y,Ma Z Q,Li G L,Khan N U,Pan Y H,Liang Y T,Zhang H L,Li J J,Li Z C. 2021. GNP6,a novel allele of MOC1,regulates panicle and tiller development in rice[J]. The Crop Journal,9(1):57-67. doi:10.1016/j.cj.2020.04.011.

    Zou J H,Chen Z X,Zhang S Y,Zhang W P,Jiang G G,Zhao X F,Zhai W X,Pan X B,Zhu L H. 2005. Characterizations and fine mapping of a mutant gene for high tillering and dwarf in rice(Oryza sativa L.)[J]. Planta,222(4):604-612. doi:10.1007/s00425-005-0007-0.

    Zou J H,Zhang S Y,Zhang W P,Li G,Chen Z X,Zhai W X,Zhao X F,Pan X P,Xie Q,Zhu L H. 2006. The rice HIGH-TILLERING DWARF1 encoding an ortholog of Arabidopsis MAX3 is required for negative regulation of the outgrowth of axillary buds[J]. The Plant Journal,48(5):687-698. doi:10.1111/j.1365-313X.2006.02916.x.

    (責(zé)任編輯 陳 燕)

    猜你喜歡
    表達(dá)分析腋芽再生稻
    外源獨(dú)腳金內(nèi)酯對(duì)煙草腋芽伸長(zhǎng)及獨(dú)腳金內(nèi)酯代謝途徑相關(guān)基因表達(dá)的影響
    “稻+再生稻+馬鈴薯”種植模式適用栽培技術(shù)
    再生稻促芽追肥機(jī)研制
    生長(zhǎng)素調(diào)控植物腋芽發(fā)育的研究進(jìn)展
    茶樹(shù)帶腋芽莖段組織培養(yǎng)研究
    糧食 超級(jí)稻再生稻平均畝產(chǎn)創(chuàng)世界紀(jì)錄
    仲丁靈對(duì)西瓜腋芽?jī)?nèi)源激素的影響
    紅花生育酚環(huán)化酶基因的克隆及表達(dá)分析
    湖北省再生稻面積達(dá)140萬(wàn)畝
    玉米紋枯病病菌y—谷氨酰轉(zhuǎn)肽酶基因克隆與表達(dá)分析
    日韩熟女老妇一区二区性免费视频| 亚洲综合色网址| 国产成人精品福利久久| 亚洲av电影在线观看一区二区三区| 丰满乱子伦码专区| 哪个播放器可以免费观看大片| 五月伊人婷婷丁香| 国产福利在线免费观看视频| 久久久欧美国产精品| 麻豆精品久久久久久蜜桃| 精品国产超薄肉色丝袜足j| 亚洲精品在线美女| 国产片特级美女逼逼视频| 一级a爱视频在线免费观看| 99热网站在线观看| 亚洲综合色网址| 国产亚洲最大av| 国产国语露脸激情在线看| 亚洲精品日韩在线中文字幕| 亚洲天堂av无毛| 欧美亚洲 丝袜 人妻 在线| 欧美人与性动交α欧美精品济南到 | 下体分泌物呈黄色| 免费黄色在线免费观看| 热re99久久精品国产66热6| 男人舔女人的私密视频| 高清在线视频一区二区三区| 亚洲三级黄色毛片| 亚洲视频免费观看视频| av视频免费观看在线观看| 韩国高清视频一区二区三区| 黄网站色视频无遮挡免费观看| 777久久人妻少妇嫩草av网站| 色婷婷久久久亚洲欧美| 十分钟在线观看高清视频www| 中国三级夫妇交换| 一边亲一边摸免费视频| 国产亚洲午夜精品一区二区久久| 国产亚洲欧美精品永久| 大码成人一级视频| 免费av中文字幕在线| 精品少妇久久久久久888优播| 亚洲国产最新在线播放| 亚洲欧美精品综合一区二区三区 | 夫妻午夜视频| 国产激情久久老熟女| 亚洲欧洲精品一区二区精品久久久 | 人人妻人人爽人人添夜夜欢视频| 中文字幕人妻丝袜一区二区 | 美女中出高潮动态图| 男人添女人高潮全过程视频| 少妇的逼水好多| 热99久久久久精品小说推荐| 国产在线视频一区二区| 日韩熟女老妇一区二区性免费视频| 国产高清不卡午夜福利| 纵有疾风起免费观看全集完整版| 欧美日韩成人在线一区二区| 国产精品秋霞免费鲁丝片| 高清av免费在线| 国产视频首页在线观看| 国产精品一区二区在线不卡| 色播在线永久视频| 热99国产精品久久久久久7| 国产精品一二三区在线看| 天天躁夜夜躁狠狠躁躁| 一区福利在线观看| 在线精品无人区一区二区三| 伦理电影大哥的女人| 久久久久久久久久久久大奶| 国产黄频视频在线观看| 成人国语在线视频| 一级毛片电影观看| 午夜日韩欧美国产| 少妇熟女欧美另类| 女人高潮潮喷娇喘18禁视频| 日韩视频在线欧美| 久久久久国产网址| 精品人妻偷拍中文字幕| 成人黄色视频免费在线看| 久久久久国产一级毛片高清牌| 丰满迷人的少妇在线观看| 日日爽夜夜爽网站| 欧美日韩精品网址| 九九爱精品视频在线观看| 精品一品国产午夜福利视频| av有码第一页| 老鸭窝网址在线观看| 国产精品偷伦视频观看了| 免费在线观看黄色视频的| 国产精品免费大片| 久久午夜福利片| 亚洲精品第二区| videos熟女内射| 亚洲精品一二三| www.av在线官网国产| 人成视频在线观看免费观看| 亚洲美女黄色视频免费看| 桃花免费在线播放| 成人漫画全彩无遮挡| 9191精品国产免费久久| 欧美精品亚洲一区二区| 丝瓜视频免费看黄片| 久热这里只有精品99| 五月开心婷婷网| 男女边摸边吃奶| 国产一区二区在线观看av| 亚洲精品国产色婷婷电影| videossex国产| 韩国高清视频一区二区三区| 国产免费又黄又爽又色| av有码第一页| 免费高清在线观看视频在线观看| 亚洲成国产人片在线观看| 看免费av毛片| 久久久久网色| 国产亚洲一区二区精品| 伊人久久大香线蕉亚洲五| 国产精品香港三级国产av潘金莲 | 午夜福利,免费看| 午夜免费男女啪啪视频观看| 亚洲av电影在线进入| 亚洲av中文av极速乱| 国产精品无大码| 大片电影免费在线观看免费| 精品少妇久久久久久888优播| 丝袜脚勾引网站| 午夜精品国产一区二区电影| 国产一区有黄有色的免费视频| 亚洲第一青青草原| 成年美女黄网站色视频大全免费| 免费人妻精品一区二区三区视频| 桃花免费在线播放| 国产综合精华液| 中文字幕精品免费在线观看视频| 欧美日韩av久久| 女的被弄到高潮叫床怎么办| 男人添女人高潮全过程视频| 欧美激情极品国产一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 国产黄色免费在线视频| 亚洲av男天堂| 伦理电影免费视频| 午夜福利视频精品| 韩国高清视频一区二区三区| 国产野战对白在线观看| 日韩在线高清观看一区二区三区| 午夜福利在线免费观看网站| 亚洲人成电影观看| 亚洲国产欧美在线一区| 桃花免费在线播放| 国产精品三级大全| 毛片一级片免费看久久久久| av国产久精品久网站免费入址| 亚洲国产精品成人久久小说| 亚洲图色成人| 人妻一区二区av| 成人毛片60女人毛片免费| 欧美人与性动交α欧美精品济南到 | 久久99精品国语久久久| 一级毛片 在线播放| 青青草视频在线视频观看| 亚洲四区av| 欧美成人精品欧美一级黄| 赤兔流量卡办理| 在线天堂最新版资源| 日本av免费视频播放| 大片电影免费在线观看免费| 高清不卡的av网站| 久久久久网色| 欧美亚洲日本最大视频资源| 老女人水多毛片| 午夜福利,免费看| 午夜福利网站1000一区二区三区| 亚洲天堂av无毛| 精品久久蜜臀av无| 成人18禁高潮啪啪吃奶动态图| 亚洲成人手机| 777久久人妻少妇嫩草av网站| 韩国精品一区二区三区| 美女午夜性视频免费| 精品酒店卫生间| 1024视频免费在线观看| 成人国产av品久久久| 赤兔流量卡办理| 黄色配什么色好看| 麻豆av在线久日| 桃花免费在线播放| 校园人妻丝袜中文字幕| 欧美国产精品va在线观看不卡| 九色亚洲精品在线播放| 乱人伦中国视频| av视频免费观看在线观看| 欧美日韩综合久久久久久| 黄片无遮挡物在线观看| 亚洲av福利一区| av国产精品久久久久影院| 日韩 亚洲 欧美在线| 熟妇人妻不卡中文字幕| 国产精品不卡视频一区二区| 中文天堂在线官网| 久久久久国产一级毛片高清牌| 中文字幕最新亚洲高清| 99九九在线精品视频| 亚洲国产欧美日韩在线播放| 另类亚洲欧美激情| 看十八女毛片水多多多| 电影成人av| 午夜久久久在线观看| 美女福利国产在线| 婷婷色av中文字幕| 久久久精品免费免费高清| 女性生殖器流出的白浆| 国产精品无大码| 男女边吃奶边做爰视频| 一级毛片我不卡| 丝袜在线中文字幕| 午夜免费男女啪啪视频观看| 国产av码专区亚洲av| 又大又黄又爽视频免费| 欧美少妇被猛烈插入视频| 免费日韩欧美在线观看| 男人舔女人的私密视频| 欧美激情 高清一区二区三区| 久久人妻熟女aⅴ| 亚洲欧美成人综合另类久久久| 欧美另类一区| 99久久精品国产国产毛片| 自线自在国产av| 亚洲男人天堂网一区| 国产福利在线免费观看视频| 最近最新中文字幕大全免费视频 | av不卡在线播放| 成人二区视频| 精品视频人人做人人爽| 老司机影院毛片| 国产精品成人在线| 亚洲av福利一区| 人妻人人澡人人爽人人| 婷婷成人精品国产| 国产又色又爽无遮挡免| 中文字幕另类日韩欧美亚洲嫩草| 最新的欧美精品一区二区| 亚洲国产av新网站| 久久人人爽av亚洲精品天堂| 97人妻天天添夜夜摸| 乱人伦中国视频| 黄色一级大片看看| 两性夫妻黄色片| 欧美少妇被猛烈插入视频| 天天躁日日躁夜夜躁夜夜| 男人舔女人的私密视频| 在线精品无人区一区二区三| freevideosex欧美| 日本猛色少妇xxxxx猛交久久| 欧美日韩精品网址| 欧美av亚洲av综合av国产av | 精品酒店卫生间| 成人黄色视频免费在线看| 欧美少妇被猛烈插入视频| 精品卡一卡二卡四卡免费| 热re99久久精品国产66热6| 久久久欧美国产精品| 成年美女黄网站色视频大全免费| 五月伊人婷婷丁香| 亚洲欧洲精品一区二区精品久久久 | 欧美日韩综合久久久久久| 五月伊人婷婷丁香| 啦啦啦中文免费视频观看日本| 最近2019中文字幕mv第一页| 天天影视国产精品| 亚洲,欧美精品.| 这个男人来自地球电影免费观看 | 美女xxoo啪啪120秒动态图| 亚洲国产色片| 少妇精品久久久久久久| 国产欧美亚洲国产| 最黄视频免费看| 国产激情久久老熟女| 亚洲欧美中文字幕日韩二区| 丝袜喷水一区| 精品国产露脸久久av麻豆| 99香蕉大伊视频| 久久精品久久精品一区二区三区| 狂野欧美激情性bbbbbb| 一级爰片在线观看| 日本免费在线观看一区| 免费播放大片免费观看视频在线观看| 综合色丁香网| 欧美亚洲日本最大视频资源| 伦精品一区二区三区| 美女大奶头黄色视频| 欧美中文综合在线视频| 成年人午夜在线观看视频| 国产精品久久久av美女十八| 午夜福利在线观看免费完整高清在| 一区二区av电影网| 久久午夜福利片| 天天躁夜夜躁狠狠躁躁| 性高湖久久久久久久久免费观看| 久久精品亚洲av国产电影网| av天堂久久9| 少妇人妻精品综合一区二区| 深夜精品福利| 男女啪啪激烈高潮av片| 又粗又硬又长又爽又黄的视频| 最近2019中文字幕mv第一页| 亚洲视频免费观看视频| 亚洲精品,欧美精品| 日韩av免费高清视频| 亚洲欧美日韩另类电影网站| 久久久欧美国产精品| 2022亚洲国产成人精品| 成人漫画全彩无遮挡| 亚洲成人一二三区av| 在线天堂中文资源库| av免费观看日本| 亚洲精品av麻豆狂野| 韩国精品一区二区三区| 亚洲成人av在线免费| 久久久久国产一级毛片高清牌| 2021少妇久久久久久久久久久| 一级爰片在线观看| 中文字幕人妻丝袜制服| 日韩大片免费观看网站| 国产精品嫩草影院av在线观看| 久热这里只有精品99| 国产精品国产三级专区第一集| 成人影院久久| 久热久热在线精品观看| 色94色欧美一区二区| 看非洲黑人一级黄片| 亚洲欧美一区二区三区久久| 欧美成人午夜免费资源| 在线天堂最新版资源| 一级片'在线观看视频| 欧美精品亚洲一区二区| 精品人妻在线不人妻| 麻豆av在线久日| 91在线精品国自产拍蜜月| 九九爱精品视频在线观看| 欧美激情高清一区二区三区 | 又粗又硬又长又爽又黄的视频| 久久精品夜色国产| 亚洲欧美成人精品一区二区| 欧美人与性动交α欧美精品济南到 | 国产亚洲午夜精品一区二区久久| 久久久久国产精品人妻一区二区| 国产精品二区激情视频| 免费播放大片免费观看视频在线观看| av卡一久久| 国产又爽黄色视频| 电影成人av| 日韩av免费高清视频| 一边亲一边摸免费视频| 2021少妇久久久久久久久久久| 精品人妻一区二区三区麻豆| 成人毛片60女人毛片免费| 久久影院123| 女人被躁到高潮嗷嗷叫费观| av有码第一页| 精品久久久久久电影网| 在现免费观看毛片| kizo精华| 如日韩欧美国产精品一区二区三区| 精品少妇久久久久久888优播| 亚洲av在线观看美女高潮| 日本-黄色视频高清免费观看| 成人国语在线视频| 欧美精品亚洲一区二区| 99re6热这里在线精品视频| 国产熟女午夜一区二区三区| 色播在线永久视频| 国产精品久久久av美女十八| 纵有疾风起免费观看全集完整版| 国产精品 欧美亚洲| 如何舔出高潮| 亚洲人成77777在线视频| 午夜久久久在线观看| 欧美亚洲日本最大视频资源| 亚洲成色77777| 日韩av免费高清视频| 黄网站色视频无遮挡免费观看| 久久久久精品人妻al黑| 亚洲一级一片aⅴ在线观看| 国产精品久久久久久精品电影小说| 免费看不卡的av| 日韩中字成人| 亚洲欧美清纯卡通| 青春草视频在线免费观看| 麻豆av在线久日| 久久国产亚洲av麻豆专区| 成年av动漫网址| 久久亚洲国产成人精品v| 如日韩欧美国产精品一区二区三区| 国产1区2区3区精品| 亚洲国产看品久久| 两性夫妻黄色片| 春色校园在线视频观看| 国产精品不卡视频一区二区| 搡女人真爽免费视频火全软件| 国产精品熟女久久久久浪| 久久久欧美国产精品| 日韩不卡一区二区三区视频在线| av不卡在线播放| 久久毛片免费看一区二区三区| 精品国产乱码久久久久久小说| 国产精品久久久久成人av| 久久久欧美国产精品| 制服诱惑二区| 婷婷色综合www| 熟女av电影| 精品一区二区免费观看| 99久国产av精品国产电影| 亚洲精华国产精华液的使用体验| 久久久精品国产亚洲av高清涩受| 侵犯人妻中文字幕一二三四区| 国产成人精品在线电影| √禁漫天堂资源中文www| 如何舔出高潮| 国产福利在线免费观看视频| av国产精品久久久久影院| 女人高潮潮喷娇喘18禁视频| freevideosex欧美| 午夜91福利影院| h视频一区二区三区| 国产精品偷伦视频观看了| 在现免费观看毛片| 日韩,欧美,国产一区二区三区| 精品久久蜜臀av无| 国产精品久久久久久精品古装| 久久人妻熟女aⅴ| 久久久久久久久免费视频了| 女人高潮潮喷娇喘18禁视频| 成人毛片a级毛片在线播放| 黄色毛片三级朝国网站| 久久精品aⅴ一区二区三区四区 | 国产精品女同一区二区软件| 久久精品国产a三级三级三级| 日韩人妻精品一区2区三区| 9191精品国产免费久久| 一区二区日韩欧美中文字幕| 国产一区二区在线观看av| 亚洲欧美精品综合一区二区三区 | 国产精品不卡视频一区二区| 日韩av不卡免费在线播放| av卡一久久| 午夜免费男女啪啪视频观看| 日日撸夜夜添| 国产精品一二三区在线看| 日本-黄色视频高清免费观看| 久久久久久久精品精品| 午夜福利视频在线观看免费| 制服诱惑二区| 2021少妇久久久久久久久久久| 欧美日本中文国产一区发布| 熟女少妇亚洲综合色aaa.| 免费黄网站久久成人精品| 精品人妻偷拍中文字幕| 久久 成人 亚洲| 亚洲人成网站在线观看播放| 久久这里只有精品19| 亚洲成人一二三区av| 亚洲国产精品一区二区三区在线| 久久 成人 亚洲| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 天堂8中文在线网| 免费少妇av软件| 亚洲国产av影院在线观看| 最近的中文字幕免费完整| 欧美日韩一区二区视频在线观看视频在线| 黄色毛片三级朝国网站| 亚洲精品美女久久久久99蜜臀 | 日本黄色日本黄色录像| 国产片内射在线| 国产淫语在线视频| 国产97色在线日韩免费| 老汉色∧v一级毛片| 亚洲精品视频女| 春色校园在线视频观看| 欧美日韩综合久久久久久| 精品第一国产精品| 精品人妻在线不人妻| 久久久久视频综合| 在线亚洲精品国产二区图片欧美| 免费高清在线观看视频在线观看| 五月天丁香电影| 最近中文字幕2019免费版| 少妇猛男粗大的猛烈进出视频| 少妇被粗大猛烈的视频| 亚洲av免费高清在线观看| 国产精品.久久久| 亚洲av电影在线观看一区二区三区| 亚洲国产看品久久| 十分钟在线观看高清视频www| 亚洲成人一二三区av| 日韩av不卡免费在线播放| 欧美国产精品va在线观看不卡| 人体艺术视频欧美日本| 中文字幕av电影在线播放| 深夜精品福利| 国产高清国产精品国产三级| 日本欧美视频一区| 桃花免费在线播放| 乱人伦中国视频| 免费播放大片免费观看视频在线观看| 激情视频va一区二区三区| 麻豆av在线久日| 五月天丁香电影| 亚洲欧美精品自产自拍| 国语对白做爰xxxⅹ性视频网站| 久久久国产精品麻豆| 中文字幕av电影在线播放| 国产一区二区三区综合在线观看| 国产视频首页在线观看| kizo精华| 美女xxoo啪啪120秒动态图| kizo精华| 国产综合精华液| 国产又色又爽无遮挡免| 97在线视频观看| 日韩一区二区视频免费看| 中文精品一卡2卡3卡4更新| 男女午夜视频在线观看| av免费观看日本| 老汉色∧v一级毛片| 青草久久国产| 国产亚洲一区二区精品| 亚洲精品,欧美精品| 欧美成人午夜精品| 日产精品乱码卡一卡2卡三| 久久久久国产精品人妻一区二区| 一级片免费观看大全| a级毛片黄视频| av.在线天堂| 亚洲精品国产av成人精品| 成人亚洲欧美一区二区av| 丝袜美足系列| 色94色欧美一区二区| 久久久久人妻精品一区果冻| 婷婷成人精品国产| 伊人亚洲综合成人网| 亚洲精品国产av成人精品| 国产免费现黄频在线看| 亚洲精品久久午夜乱码| 欧美激情 高清一区二区三区| 亚洲精品国产色婷婷电影| 亚洲人成网站在线观看播放| 国产成人精品婷婷| 男女无遮挡免费网站观看| 亚洲精品乱久久久久久| 97在线人人人人妻| 国产麻豆69| 伊人久久国产一区二区| 免费在线观看黄色视频的| 飞空精品影院首页| 一区二区三区精品91| 亚洲欧美精品综合一区二区三区 | 婷婷色综合大香蕉| 建设人人有责人人尽责人人享有的| 久久久久久久久久久久大奶| 久久久久网色| 国产女主播在线喷水免费视频网站| 亚洲一区中文字幕在线| 亚洲av欧美aⅴ国产| 在现免费观看毛片| 亚洲av日韩在线播放| 欧美精品亚洲一区二区| 久久99精品国语久久久| 超碰97精品在线观看| 国产国语露脸激情在线看| 精品亚洲成a人片在线观看| 欧美精品国产亚洲| 亚洲视频免费观看视频| 午夜免费鲁丝| 女人久久www免费人成看片| 国产免费一区二区三区四区乱码| 香蕉丝袜av| a 毛片基地| 一区二区三区精品91| 亚洲国产欧美日韩在线播放| 亚洲美女视频黄频| 免费少妇av软件| 丰满乱子伦码专区| 欧美+日韩+精品| 好男人视频免费观看在线| 国产有黄有色有爽视频| 曰老女人黄片| 欧美日韩成人在线一区二区| 成年av动漫网址| 久久久久国产一级毛片高清牌| 两个人免费观看高清视频| 丰满迷人的少妇在线观看| 国产一区二区激情短视频 | 中文字幕制服av| 桃花免费在线播放| 国产日韩欧美在线精品| www.熟女人妻精品国产| 十八禁网站网址无遮挡| 9191精品国产免费久久| 国产福利在线免费观看视频| 日韩一区二区三区影片| 日韩伦理黄色片| 一个人免费看片子| 91久久精品国产一区二区三区| 久久久久国产网址| 高清欧美精品videossex| 欧美精品人与动牲交sv欧美| 婷婷色av中文字幕| 久久久久久伊人网av| 国产又爽黄色视频| 亚洲国产精品一区三区| 美女高潮到喷水免费观看| 高清欧美精品videossex| 亚洲av中文av极速乱|