• <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á)分析
    久久午夜综合久久蜜桃| 成人毛片a级毛片在线播放| 波多野结衣一区麻豆| 亚洲欧洲日产国产| 国产 精品1| 精品人妻一区二区三区麻豆| 亚洲成色77777| 日本av免费视频播放| 国产片内射在线| 久久久久网色| 精品人妻在线不人妻| av.在线天堂| 国产一区二区在线观看av| 啦啦啦在线免费观看视频4| 青春草国产在线视频| 一区二区日韩欧美中文字幕| 日本-黄色视频高清免费观看| 国产精品一国产av| 黄片播放在线免费| 免费黄色在线免费观看| 亚洲人成网站在线观看播放| 亚洲激情五月婷婷啪啪| 成人漫画全彩无遮挡| 少妇 在线观看| 久久婷婷青草| 视频在线观看一区二区三区| 男女边吃奶边做爰视频| 97精品久久久久久久久久精品| 欧美97在线视频| 黄频高清免费视频| 2018国产大陆天天弄谢| 久久 成人 亚洲| 五月天丁香电影| 亚洲国产精品成人久久小说| 亚洲国产色片| 精品国产超薄肉色丝袜足j| 三上悠亚av全集在线观看| 国产日韩欧美亚洲二区| 两个人免费观看高清视频| 看免费成人av毛片| 中文字幕人妻丝袜制服| 这个男人来自地球电影免费观看 | 黑丝袜美女国产一区| 纵有疾风起免费观看全集完整版| 免费观看a级毛片全部| 国产亚洲欧美精品永久| 国产精品国产三级国产专区5o| 国产精品国产av在线观看| 边亲边吃奶的免费视频| 国产精品免费大片| a 毛片基地| 久久婷婷青草| 亚洲美女黄色视频免费看| 9色porny在线观看| 纯流量卡能插随身wifi吗| 日韩,欧美,国产一区二区三区| 日本色播在线视频| 一二三四中文在线观看免费高清| 18禁观看日本| 国产免费现黄频在线看| 精品久久蜜臀av无| 一级黄片播放器| 热99久久久久精品小说推荐| 免费人妻精品一区二区三区视频| 久久久久久久精品精品| av国产久精品久网站免费入址| 国产精品久久久久久精品古装| 国产黄色视频一区二区在线观看| 久久久精品国产亚洲av高清涩受| 交换朋友夫妻互换小说| 精品国产乱码久久久久久小说| 久久久久久伊人网av| 观看av在线不卡| av免费观看日本| 成人国产麻豆网| 99热网站在线观看| 亚洲内射少妇av| 国产老妇伦熟女老妇高清| 精品一区二区三卡| 亚洲激情五月婷婷啪啪| 黑人欧美特级aaaaaa片| 国产欧美日韩一区二区三区在线| 日本免费在线观看一区| 91久久精品国产一区二区三区| 欧美成人午夜精品| 国产在线免费精品| 日韩大片免费观看网站| 亚洲五月色婷婷综合| 亚洲精品一二三| 日本色播在线视频| 99久久精品国产国产毛片| 一级爰片在线观看| 日韩人妻精品一区2区三区| 精品一区在线观看国产| 少妇 在线观看| 国产av一区二区精品久久| 超色免费av| 青青草视频在线视频观看| 国产精品蜜桃在线观看| 亚洲中文av在线| av片东京热男人的天堂| 国产片特级美女逼逼视频| 一区福利在线观看| 日本av免费视频播放| 国产1区2区3区精品| 亚洲成人一二三区av| 久久久a久久爽久久v久久| 亚洲,欧美精品.| 99精国产麻豆久久婷婷| freevideosex欧美| 99热网站在线观看| 9色porny在线观看| 国产成人免费无遮挡视频| 九色亚洲精品在线播放| 久久久a久久爽久久v久久| 欧美日本中文国产一区发布| 亚洲av欧美aⅴ国产| 亚洲av.av天堂| 91在线精品国自产拍蜜月| 少妇人妻精品综合一区二区| 精品人妻在线不人妻| 男的添女的下面高潮视频| 搡老乐熟女国产| 国产无遮挡羞羞视频在线观看| 欧美激情极品国产一区二区三区| 777米奇影视久久| 亚洲五月色婷婷综合| 精品人妻熟女毛片av久久网站| 日韩伦理黄色片| 天天躁日日躁夜夜躁夜夜| 日本欧美视频一区| 久久精品亚洲av国产电影网| 欧美日韩国产mv在线观看视频| 老熟女久久久| 久久免费观看电影| 国产免费视频播放在线视频| 欧美精品高潮呻吟av久久| 亚洲激情五月婷婷啪啪| 午夜91福利影院| 日韩精品免费视频一区二区三区| 精品国产乱码久久久久久小说| 一本久久精品| 午夜福利,免费看| 久久久久精品久久久久真实原创| 丝袜脚勾引网站| 久久av网站| 亚洲欧洲国产日韩| 水蜜桃什么品种好| 人人妻人人澡人人爽人人夜夜| 一本—道久久a久久精品蜜桃钙片| 777久久人妻少妇嫩草av网站| 亚洲男人天堂网一区| 日韩制服骚丝袜av| 一级黄片播放器| 五月开心婷婷网| 国产一区亚洲一区在线观看| 亚洲欧美日韩另类电影网站| 免费高清在线观看日韩| 丰满迷人的少妇在线观看| 成年av动漫网址| 午夜福利在线观看免费完整高清在| www.熟女人妻精品国产| 日日啪夜夜爽| 国产精品无大码| 9191精品国产免费久久| 国产不卡av网站在线观看| 精品国产乱码久久久久久小说| 2021少妇久久久久久久久久久| 亚洲人成网站在线观看播放| 桃花免费在线播放| 性色av一级| 90打野战视频偷拍视频| 制服人妻中文乱码| 久久久a久久爽久久v久久| 汤姆久久久久久久影院中文字幕| 亚洲综合精品二区| 十八禁网站网址无遮挡| 午夜久久久在线观看| 国产一级毛片在线| 男女下面插进去视频免费观看| 精品第一国产精品| 久久精品国产鲁丝片午夜精品| 亚洲精品aⅴ在线观看| 精品少妇久久久久久888优播| 一个人免费看片子| 99热国产这里只有精品6| 老司机亚洲免费影院| 欧美老熟妇乱子伦牲交| 精品午夜福利在线看| 国产乱来视频区| 亚洲av电影在线进入| 欧美日韩亚洲国产一区二区在线观看 | 免费黄频网站在线观看国产| 一区在线观看完整版| 欧美 亚洲 国产 日韩一| 18禁国产床啪视频网站| 欧美精品国产亚洲| 免费看av在线观看网站| 交换朋友夫妻互换小说| 午夜福利影视在线免费观看| 亚洲国产日韩一区二区| 午夜老司机福利剧场| 国产精品一二三区在线看| 日本黄色日本黄色录像| 丰满乱子伦码专区| 黑人巨大精品欧美一区二区蜜桃| 免费人妻精品一区二区三区视频| 国产精品久久久久久精品电影小说| 看免费成人av毛片| 亚洲内射少妇av| 纯流量卡能插随身wifi吗| 亚洲一码二码三码区别大吗| 欧美日韩国产mv在线观看视频| kizo精华| xxx大片免费视频| 亚洲国产色片| 久久狼人影院| 久久人妻熟女aⅴ| 99香蕉大伊视频| 亚洲成国产人片在线观看| 国产高清国产精品国产三级| av不卡在线播放| 在线观看www视频免费| 熟女av电影| 亚洲精品aⅴ在线观看| 久久久精品国产亚洲av高清涩受| 国产日韩一区二区三区精品不卡| 曰老女人黄片| 少妇的逼水好多| 男女啪啪激烈高潮av片| 欧美另类一区| 天堂俺去俺来也www色官网| 老汉色av国产亚洲站长工具| 免费播放大片免费观看视频在线观看| 啦啦啦中文免费视频观看日本| 极品人妻少妇av视频| 1024香蕉在线观看| 久久久久精品久久久久真实原创| 曰老女人黄片| 看免费成人av毛片| 99久久中文字幕三级久久日本| 韩国av在线不卡| 国产免费现黄频在线看| 成年女人毛片免费观看观看9 | 美女视频免费永久观看网站| 黄色一级大片看看| 街头女战士在线观看网站| 午夜免费鲁丝| 久久久久久人妻| 永久网站在线| 天美传媒精品一区二区| 97在线人人人人妻| 午夜福利网站1000一区二区三区| av免费在线看不卡| 欧美人与性动交α欧美软件| 精品国产一区二区三区四区第35| 老鸭窝网址在线观看| 国产亚洲最大av| 亚洲精品中文字幕在线视频| 久久久国产精品麻豆| 日韩大片免费观看网站| a级毛片黄视频| 热99久久久久精品小说推荐| 如日韩欧美国产精品一区二区三区| 国精品久久久久久国模美| 999久久久国产精品视频| 国产精品一区二区在线观看99| 亚洲av男天堂| 女性被躁到高潮视频| 成年女人在线观看亚洲视频| 久久99蜜桃精品久久| 久久久久久久精品精品| 女人高潮潮喷娇喘18禁视频| 美女高潮到喷水免费观看| 捣出白浆h1v1| 免费播放大片免费观看视频在线观看| 亚洲精品日本国产第一区| 99国产精品免费福利视频| 日韩大片免费观看网站| 日韩一区二区三区影片| 欧美日韩综合久久久久久| 久久久久久久久免费视频了| 97在线人人人人妻| 国产女主播在线喷水免费视频网站| 国产人伦9x9x在线观看 | 国产免费福利视频在线观看| av在线播放精品| 亚洲三级黄色毛片| 久久精品国产自在天天线| 亚洲欧美清纯卡通| 热99久久久久精品小说推荐| 日韩精品有码人妻一区| 少妇人妻精品综合一区二区| 亚洲经典国产精华液单| 久久影院123| 精品亚洲乱码少妇综合久久| 国产黄色免费在线视频| 69精品国产乱码久久久| 在线观看美女被高潮喷水网站| 精品亚洲乱码少妇综合久久| 免费高清在线观看视频在线观看| 午夜日韩欧美国产| a级片在线免费高清观看视频| 日韩精品有码人妻一区| 成年人午夜在线观看视频| 亚洲国产精品999| 极品人妻少妇av视频| 国产乱来视频区| 狠狠婷婷综合久久久久久88av| 欧美少妇被猛烈插入视频| 国产在线一区二区三区精| 黄频高清免费视频| 国产精品一区二区在线不卡| 欧美成人午夜免费资源| 肉色欧美久久久久久久蜜桃| 啦啦啦中文免费视频观看日本| 尾随美女入室| 天堂8中文在线网| 国产成人午夜福利电影在线观看| 欧美精品av麻豆av| 久久毛片免费看一区二区三区| 777米奇影视久久| 精品亚洲成国产av| 日日撸夜夜添| tube8黄色片| 国产成人免费观看mmmm| 日本91视频免费播放| 精品一区二区免费观看| 自线自在国产av| 日本色播在线视频| 妹子高潮喷水视频| 黄片无遮挡物在线观看| 日韩中文字幕欧美一区二区 | 大香蕉久久网| 色视频在线一区二区三区| 亚洲第一青青草原| 9色porny在线观看| 国产探花极品一区二区| 亚洲情色 制服丝袜| 精品人妻在线不人妻| 亚洲精品国产一区二区精华液| 久久久a久久爽久久v久久| 美女国产视频在线观看| 午夜精品国产一区二区电影| 日韩伦理黄色片| 建设人人有责人人尽责人人享有的| av网站免费在线观看视频| 色婷婷久久久亚洲欧美| 亚洲一码二码三码区别大吗| 777久久人妻少妇嫩草av网站| 永久免费av网站大全| 欧美日韩国产mv在线观看视频| 90打野战视频偷拍视频| 热re99久久国产66热| 精品人妻熟女毛片av久久网站| 丰满迷人的少妇在线观看| 18+在线观看网站| 国产精品欧美亚洲77777| 女人高潮潮喷娇喘18禁视频| 免费不卡的大黄色大毛片视频在线观看| 十八禁高潮呻吟视频| 极品人妻少妇av视频| 999久久久国产精品视频| 看非洲黑人一级黄片| 三级国产精品片| 国产免费现黄频在线看| 亚洲av日韩在线播放| 国产成人精品久久二区二区91 | a级片在线免费高清观看视频| 又大又黄又爽视频免费| 久久久久网色| 久久久欧美国产精品| 边亲边吃奶的免费视频| 在线观看国产h片| 中文字幕制服av| a级毛片黄视频| 国产一区二区激情短视频 | 色播在线永久视频| 另类亚洲欧美激情| 精品一区二区三区四区五区乱码 | 大话2 男鬼变身卡| 综合色丁香网| 成年女人在线观看亚洲视频| 欧美bdsm另类| 国产 精品1| 午夜福利在线免费观看网站| 王馨瑶露胸无遮挡在线观看| 色网站视频免费| 一本大道久久a久久精品| 久久久久久久久久人人人人人人| 中文字幕精品免费在线观看视频| av.在线天堂| 亚洲精品日本国产第一区| 亚洲一码二码三码区别大吗| 9热在线视频观看99| 纵有疾风起免费观看全集完整版| 精品少妇一区二区三区视频日本电影 | 十八禁网站网址无遮挡| 国产成人免费无遮挡视频| 性高湖久久久久久久久免费观看| 日韩欧美精品免费久久| 免费观看在线日韩| 丝袜脚勾引网站| a级毛片在线看网站| 国产日韩欧美视频二区| 街头女战士在线观看网站| 黄色 视频免费看| 日韩中字成人| 97在线人人人人妻| 国产一区亚洲一区在线观看| 考比视频在线观看| 久久99热这里只频精品6学生| 精品人妻熟女毛片av久久网站| 黄色毛片三级朝国网站| 国产av精品麻豆| 国产精品熟女久久久久浪| 亚洲国产精品一区二区三区在线| 国产成人aa在线观看| 亚洲av国产av综合av卡| 久久久久久人人人人人| 亚洲美女视频黄频| 巨乳人妻的诱惑在线观看| 欧美日韩亚洲国产一区二区在线观看 | 中文字幕av电影在线播放| 国产欧美亚洲国产| 国产在视频线精品| 成人毛片a级毛片在线播放| 成年人午夜在线观看视频| 国产精品 国内视频| 天堂俺去俺来也www色官网| 欧美日韩视频精品一区| 九九爱精品视频在线观看| 亚洲国产av影院在线观看| 久久久久久久亚洲中文字幕| 男女边吃奶边做爰视频| 国产探花极品一区二区| 性色avwww在线观看| 黄频高清免费视频| 秋霞在线观看毛片| 亚洲av欧美aⅴ国产| 在线观看国产h片| 日韩av免费高清视频| 亚洲伊人色综图| 日韩精品有码人妻一区| 婷婷成人精品国产| 丝袜美腿诱惑在线| av卡一久久| 精品酒店卫生间| 久久 成人 亚洲| 亚洲av电影在线进入| 国产免费视频播放在线视频| 黄色一级大片看看| 久久97久久精品| 女人久久www免费人成看片| 国精品久久久久久国模美| 免费大片黄手机在线观看| 亚洲精品久久午夜乱码| 高清视频免费观看一区二区| 中文精品一卡2卡3卡4更新| 欧美日韩视频精品一区| 黄片播放在线免费| 涩涩av久久男人的天堂| 亚洲国产av影院在线观看| 999久久久国产精品视频| 26uuu在线亚洲综合色| 性色av一级| 青草久久国产| 亚洲精品第二区| 国产福利在线免费观看视频| av电影中文网址| 女性生殖器流出的白浆| 亚洲,一卡二卡三卡| 免费在线观看完整版高清| 中文字幕最新亚洲高清| 久久狼人影院| 久久精品国产亚洲av天美| av福利片在线| 免费看av在线观看网站| 日韩av在线免费看完整版不卡| 热re99久久国产66热| 伦理电影大哥的女人| 99久久综合免费| 国产精品一区二区在线观看99| 一级毛片我不卡| 黄频高清免费视频| 男人爽女人下面视频在线观看| 寂寞人妻少妇视频99o| 2021少妇久久久久久久久久久| 亚洲第一区二区三区不卡| www日本在线高清视频| 一级毛片黄色毛片免费观看视频| 最近手机中文字幕大全| 亚洲经典国产精华液单| 亚洲国产最新在线播放| 狠狠婷婷综合久久久久久88av| www日本在线高清视频| 中文天堂在线官网| 国产 精品1| 电影成人av| 麻豆精品久久久久久蜜桃| 一区二区三区精品91| 热re99久久精品国产66热6| 亚洲天堂av无毛| 久久精品国产a三级三级三级| 极品人妻少妇av视频| 欧美精品高潮呻吟av久久| 又粗又硬又长又爽又黄的视频| 一区二区日韩欧美中文字幕| 97在线视频观看| 伦精品一区二区三区| 免费观看性生交大片5| 黄色配什么色好看| 成年女人毛片免费观看观看9 | 亚洲精品一二三| 久久综合国产亚洲精品| 欧美+日韩+精品| 日本av手机在线免费观看| 色94色欧美一区二区| av电影中文网址| 人成视频在线观看免费观看| 亚洲国产精品成人久久小说| 寂寞人妻少妇视频99o| 国产精品嫩草影院av在线观看| 国产福利在线免费观看视频| 五月天丁香电影| 欧美精品一区二区大全| 欧美 亚洲 国产 日韩一| 尾随美女入室| 国产高清国产精品国产三级| 91精品三级在线观看| 综合色丁香网| 国产精品.久久久| 午夜91福利影院| 新久久久久国产一级毛片| 午夜福利网站1000一区二区三区| 欧美人与性动交α欧美软件| 国产成人精品久久久久久| www日本在线高清视频| 日韩一本色道免费dvd| 97人妻天天添夜夜摸| 麻豆av在线久日| 秋霞在线观看毛片| 亚洲熟女精品中文字幕| 伊人久久大香线蕉亚洲五| 极品人妻少妇av视频| 国产精品麻豆人妻色哟哟久久| 国产男人的电影天堂91| 美女中出高潮动态图| 亚洲欧洲精品一区二区精品久久久 | 人人妻人人添人人爽欧美一区卜| 午夜福利,免费看| 久久ye,这里只有精品| 在线精品无人区一区二区三| 日本猛色少妇xxxxx猛交久久| 欧美变态另类bdsm刘玥| 国产日韩一区二区三区精品不卡| 1024香蕉在线观看| 青春草亚洲视频在线观看| 国产毛片在线视频| 国产1区2区3区精品| 国产有黄有色有爽视频| 色婷婷久久久亚洲欧美| 伊人亚洲综合成人网| 侵犯人妻中文字幕一二三四区| 中文天堂在线官网| 国产片特级美女逼逼视频| 日韩av不卡免费在线播放| 亚洲欧美精品综合一区二区三区 | 欧美精品av麻豆av| 欧美+日韩+精品| 大片免费播放器 马上看| 你懂的网址亚洲精品在线观看| 国产 精品1| av在线观看视频网站免费| 午夜精品国产一区二区电影| 色婷婷av一区二区三区视频| 精品国产一区二区三区四区第35| 国产成人精品婷婷| 日本欧美视频一区| 日韩av免费高清视频| 80岁老熟妇乱子伦牲交| 色视频在线一区二区三区| 国产一区二区三区综合在线观看| 精品少妇内射三级| 美女高潮到喷水免费观看| 免费在线观看视频国产中文字幕亚洲 | 久久精品人人爽人人爽视色| 国产在视频线精品| a级毛片在线看网站| 中文字幕制服av| 日韩三级伦理在线观看| 蜜桃在线观看..| av国产精品久久久久影院| 夫妻性生交免费视频一级片| 日本欧美国产在线视频| 视频在线观看一区二区三区| 热re99久久精品国产66热6| 精品一区二区免费观看| 国产淫语在线视频| av国产精品久久久久影院| 90打野战视频偷拍视频| 久久久久精品性色| 欧美日韩视频精品一区| 在线观看人妻少妇| 夫妻午夜视频| 捣出白浆h1v1| 一本大道久久a久久精品| 久久精品亚洲av国产电影网| 高清在线视频一区二区三区| 大片电影免费在线观看免费| 精品一区二区三卡| 亚洲精品久久午夜乱码| 亚洲成人av在线免费| 欧美日韩av久久| 99热国产这里只有精品6|