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

    孕穗期低溫脅迫下寒地水稻轉(zhuǎn)錄組分析

    2024-12-31 00:00:00郭震華馬文東潘國(guó)君杜曉東周雪松陸文靜蔡麗君劉傳雪馮延江
    江蘇農(nóng)業(yè)科學(xué) 2024年19期
    關(guān)鍵詞:轉(zhuǎn)錄組寒地水稻

    摘要:為了明確寒地水稻孕穗期低溫應(yīng)答機(jī)制,以孕穗期耐冷性差異較大的龍粳25、龍粳11水稻品種為材料,對(duì)幼穗穎花開展轉(zhuǎn)錄組測(cè)序研究。結(jié)果表明,在龍粳11中共檢測(cè)到9 468個(gè)差異表達(dá)基因(DEG),在龍粳25中共檢測(cè)到6 784個(gè)DEG。其中3 508個(gè)DEG為二者共同表達(dá),3 276個(gè)DEG在龍粳25中特異表達(dá)。GO富集分析表明,龍粳25、龍粳11的共有DEG主要富集在生物進(jìn)程調(diào)控、氧化還原反應(yīng)、非生物刺激應(yīng)答等類別中,KEGG通路主要富集在植物激素信號(hào)轉(zhuǎn)導(dǎo)、碳代謝、淀粉蔗糖代謝等途徑。而龍粳25的特有DEG多富集在生物進(jìn)程調(diào)控、非生物刺激應(yīng)答、氧化還原反應(yīng)進(jìn)程等類別,KEGG通路主要富集在植物激素信號(hào)轉(zhuǎn)導(dǎo)、內(nèi)質(zhì)網(wǎng)蛋白加工、苯丙烷合成等途徑。綜上表明,寒地水稻對(duì)孕穗期的低溫應(yīng)答是一個(gè)復(fù)雜過程,包含多個(gè)代謝調(diào)控途徑及多個(gè)相應(yīng)基因。

    關(guān)鍵詞:寒地;水稻;轉(zhuǎn)錄組;孕穗期低溫

    中圖分類號(hào):S511.01" 文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1002-1302(2024)19-0034-07

    收稿日期:2023-10-31

    基金項(xiàng)目:黑龍江省博士后面上資助項(xiàng)目(編號(hào):LBH-Z21027);黑龍江省省屬科研業(yè)務(wù)費(fèi)項(xiàng)目(編號(hào):CZKYF2023-1-C007、CZKYF2021-2-C002);黑龍江省農(nóng)業(yè)科技創(chuàng)新跨越工程重大需求科技創(chuàng)新攻關(guān)項(xiàng)目(編號(hào):CX23ZD01);黑龍江省自然科學(xué)基金優(yōu)秀青年項(xiàng)目(編號(hào):YQ2021C033)。

    作者簡(jiǎn)介:郭震華(1985—),男,內(nèi)蒙古集寧人,博士,副研究員,主要從事水稻分子育種研究。E-mail:hljsdsgzh@163.com。

    通信作者:劉傳雪,碩士,研究員,主要從事水稻育種研究,E-mail:liuchuanxue2007@163.com;馮延江,博士,研究員,主要從事水稻育種栽培研究,E-mail:zixuanfeng2008@163.com。

    水稻是全世界最重要的糧食作物之一[1,全球一半以上的人口以水稻為主食。水稻起源地多為熱帶或亞熱帶地區(qū),因此對(duì)低溫十分敏感。通常認(rèn)為,秈稻低于18 ℃,粳稻低于15 ℃,就會(huì)對(duì)其正常生長(zhǎng)造成不利影響[2。目前,水稻冷害已然成為世界性問題,其中尤以日本、韓國(guó)及我國(guó)東北稻區(qū)最為突出[3。黑龍江省地處我國(guó)最北部,屬于寒地稻作區(qū),是世界水稻栽培的北限,年平均氣溫為全國(guó)最低;低溫冷害是限制黑龍江水稻安全生產(chǎn)的重要因素。

    水稻在全生育期內(nèi)都可能受到低溫冷害的威脅。孕穗期是花粉形成的關(guān)鍵時(shí)期,若遭受低溫會(huì)直接影響花粉發(fā)育,造成花粉不育或花粉管無法萌發(fā),導(dǎo)致結(jié)實(shí)率下降,最終造成水稻減產(chǎn)[4。具體來說,以減數(shù)分裂期對(duì)低溫最為敏感;若此時(shí)遭遇低溫,水稻花粉壁上的絨氈層細(xì)胞不能正常降解,影響淀粉在花粉中的正常轉(zhuǎn)運(yùn),花粉發(fā)育所需的能量供給不足,最終造成花粉敗育5-6

    低溫脅迫除對(duì)結(jié)實(shí)率等直觀表型有影響外,對(duì)作物各方面生理功能都會(huì)造成一系列的影響,包括光合作用、能量代謝、滲透調(diào)節(jié)、過氧化物酶系統(tǒng)保護(hù)等[7-8。其中涉及到的不同代謝途徑中多個(gè)調(diào)控基因,均被報(bào)道在低溫脅迫應(yīng)答中起到不同的作用。Ctb1則是第1個(gè)通過圖位克隆得到的水稻孕穗期耐冷調(diào)控基因,該基因編碼1個(gè)F-box蛋白,通過與E3泛素連接酶亞基Skp1間的互作,來啟動(dòng)耐冷應(yīng)答反應(yīng)[9。Osmyb4被稱作水稻耐冷應(yīng)答的總開關(guān),低溫脅迫下正向調(diào)控下游多個(gè)基因的表達(dá)[10。過氧化物酶基因OsPOX1、抗壞血酸過氧化物酶基因OsAPXa、過氧化氫酶基因OsCATC等均受低溫誘導(dǎo)[11-13。

    轉(zhuǎn)錄組測(cè)序技術(shù)是近年來在基因芯片、基因表達(dá)譜等技術(shù)上發(fā)展起來的新一代測(cè)序技術(shù)。與之前的測(cè)序方法相比,轉(zhuǎn)錄組測(cè)序通量高,能夠快速全面地對(duì)某一物種的特定組織或器官在某一特定狀態(tài)下進(jìn)行全轉(zhuǎn)錄本測(cè)序分析。目前此項(xiàng)技術(shù)已經(jīng)成功應(yīng)用在多種作物中[14-16。在水稻低溫脅迫應(yīng)答方面,亦有眾多相關(guān)報(bào)道17-18。本研究針對(duì)黑龍江的生態(tài)類型和氣候條件,以當(dāng)?shù)卦兴肫谀屠湫圆町愝^大的粳稻品種龍粳25(耐冷)、龍粳11(不耐冷)作為研究對(duì)象,通過在孕穗期低溫脅迫下的轉(zhuǎn)錄組測(cè)序分析,進(jìn)一步闡明黑龍江水稻孕穗期耐冷性的遺傳調(diào)控分子機(jī)制,探索發(fā)掘寒地水稻孕穗期耐冷調(diào)控基因,以期為后續(xù)耐冷分子育種提供理論參考。

    1 材料與方法

    1.1 試驗(yàn)材料

    試驗(yàn)材料為黑龍江省農(nóng)業(yè)科學(xué)院水稻研究所育成并提供的龍粳25、龍粳11,二者均為黑龍江省第三積溫區(qū)適宜種植品種,孕穗期耐冷性差異明顯。

    1.2 材料種植

    本研究采用盆栽試驗(yàn)。試驗(yàn)地點(diǎn)為黑龍江省農(nóng)業(yè)科學(xué)院水稻研究所(下簡(jiǎn)稱水稻所)田間試驗(yàn)場(chǎng)。盆栽所需土壤均采自水稻所田間。土壤進(jìn)行干燥粉碎并過篩混勻后,裝入直徑34 cm、高35 cm的塑料桶中,每桶裝干土量約為12 kg。試驗(yàn)材料于2018年4月20日播種,待生長(zhǎng)到3葉1心時(shí)(5月15日),挑選長(zhǎng)勢(shì)一致的幼苗移栽至塑料桶中,每桶均勻栽植15株。每品種分別種植12桶,各取6桶作為對(duì)照。生長(zhǎng)過程中,為保證各單株生育進(jìn)程一致,所有材料均不留分蘗,只留主莖。孕穗期判定依據(jù)Sato等的方法[12,略作改動(dòng)。當(dāng)劍葉與倒2葉的葉枕距達(dá)-4~0 cm時(shí),即認(rèn)為其進(jìn)入孕穗期,并對(duì)該單株掛牌標(biāo)記。將同時(shí)進(jìn)入孕穗期的試驗(yàn)材料搬至光照培養(yǎng)箱(HPG-280,哈爾濱市東聯(lián)電子技術(shù)開發(fā)有限公司生產(chǎn))中,12 ℃低溫持續(xù)處理 2 d,相對(duì)濕度為75%,光—暗周期設(shè)定為8 h—16 h。低溫處理結(jié)束后,將試驗(yàn)材料移至正常環(huán)境生長(zhǎng)至成熟。收獲時(shí)計(jì)算結(jié)實(shí)率。

    1.3 材料取樣及RNA抽提

    取水稻穎花作為試驗(yàn)材料進(jìn)行轉(zhuǎn)錄組測(cè)序分析。分別在未處理和低溫處理2 d的龍粳25、龍粳11中選取生育進(jìn)程一致的水稻主莖,剝出幼穗,取幼穗上1/3的穎花,從中挑取長(zhǎng)度為3~5 mm 的穎花,稱重約0.5 g,錫箔紙包裹,過液氮迅速冷凍,于-80 ℃保存,每處理重復(fù)3次。RNA抽提采用Trizol試劑盒(Thermo Fisher Scientific公司生產(chǎn))的,方法參照試劑盒提取說明,RNA的純凈和完整程度通過瓊脂糖凝膠電泳檢測(cè)。

    1.4 文庫(kù)構(gòu)建及測(cè)序

    文庫(kù)構(gòu)建及轉(zhuǎn)錄組測(cè)序委托北京百邁克生物技術(shù)有限公司(www.biomarker.com.cn)進(jìn)行。文庫(kù)構(gòu)建完成后,對(duì)其質(zhì)量進(jìn)行檢測(cè),檢測(cè)合格后的文庫(kù)按照目標(biāo)下機(jī)數(shù)據(jù)量進(jìn)行pooling,用Illumina HiSeq 2500平臺(tái)根據(jù)操作標(biāo)準(zhǔn)進(jìn)行測(cè)序,產(chǎn)生 150 bp 的雙端reads。

    1.5 測(cè)序數(shù)據(jù)過濾及組裝

    測(cè)序后下機(jī)產(chǎn)生的原始數(shù)據(jù)(Raw Data)通常以FASTQ格式提供。過濾掉包含接頭的低質(zhì)量數(shù)據(jù)后即為可用數(shù)據(jù)(clean reads),通過Q20、Q30、GC 含量以及序列重復(fù)水平等參數(shù)來評(píng)估所得可用數(shù)據(jù)能否用于進(jìn)一步分析。以O(shè)ryza sativa L. ssp. japonica (Oryza sativa IRGSP-1.0)作為本試驗(yàn)水稻的參考基因組,通過TopHat2將可用數(shù)據(jù)比對(duì)到參考基因組上。以FPKM 值判定轉(zhuǎn)錄本或基因表達(dá)水平。

    1.6 差異表達(dá)基因的篩選及生物信息學(xué)分析

    通過DEGseq2軟件分析基因表達(dá)水平。參數(shù)設(shè)置具體為|log2 Fold Change (FC)|≥1,1 discovery rate (FDR)≤0.01。GO功能分類中基因以及基因數(shù)目的篩選是基于GO數(shù)據(jù)庫(kù)(http://www.geneontology.org/)的映射完成。GO功能富集是通過GOSeq R軟件包,以Wallenius 非中心超幾何分布檢驗(yàn)進(jìn)行分析。KEGG Pathway 的富集分析是通過KOBAS軟件(http://kobas.cbi.pku.edu.cn/home.do)進(jìn)行。GO及KEGG的富集標(biāo)準(zhǔn)皆為P≤0.01。

    1.7 差異表達(dá)基因的qRT-PCR分析

    選用Thermo Fisher Scientific公司的Trizol試劑盒,分別對(duì)龍粳25、龍粳11在處理前(0 d)和處理后(2 d)的穎花進(jìn)行抽提RNA,穎花選取部位及方法同“1.3”節(jié)。純化后的RNA通過 HiScript Q RT SuperMix for qPCR試劑盒反轉(zhuǎn)錄合成 cDNA。選取表達(dá)差異大的基因,通過primer3軟件設(shè)計(jì)qRT-PCR引物(表1),以actin為內(nèi)參。qRT-PCR反應(yīng)在ABI 7500 qPCR thermal cycler (Applied Biosystems Inc.,Carlsbad,CA,USA)儀器上運(yùn)行,重復(fù)3次,通過2-ΔΔCT法計(jì)算基因的相對(duì)表達(dá)量變化。

    1.8 數(shù)據(jù)處理

    用WPS Office分析數(shù)據(jù),用SPSS 22.0進(jìn)行顯著性檢驗(yàn)(α=0.05),用GraphPad Prism 9.0作圖。

    2 結(jié)果與分析

    2.1 孕穗期低溫下結(jié)實(shí)率變化

    孕穗期低溫直接影響水稻的結(jié)實(shí)率。通過對(duì)龍粳25、龍粳11在孕穗期進(jìn)行12 ℃低溫處理2 d,由圖1可以看出,未經(jīng)過低溫處理時(shí),龍粳25、龍粳11的結(jié)實(shí)率分別為94.5%、93.05%,二者間未有明顯差異。當(dāng)?shù)蜏靥幚? d后,龍粳25的結(jié)實(shí)率為89.27% 較對(duì)照未有明顯降低。而龍粳11在低溫處理2 d后,結(jié)實(shí)率降為37.85%,較對(duì)照呈極顯著降低??梢钥闯?,龍粳25、龍粳11在正常環(huán)境中結(jié)實(shí)率沒有明顯差異,低溫是造成二者結(jié)實(shí)率差異的原因,并且龍粳25的孕穗期耐冷性要明顯強(qiáng)于龍粳11。

    2.2 測(cè)序結(jié)果組裝分析

    本研究通過對(duì)龍粳25、龍粳11低溫處理2 d的穎花進(jìn)行轉(zhuǎn)錄組測(cè)序分析,共構(gòu)建12個(gè)文庫(kù),去掉raw reads,共得到689.23 Mb的clean reads,GC含量為54.52%~56.58%,Q30均高于91.10%,平均81.70%的clean reads比對(duì)到參考基因組上(表2)。以上結(jié)果表明,龍粳25、龍粳11的轉(zhuǎn)錄組測(cè)序質(zhì)量較高,可以繼續(xù)進(jìn)行后續(xù)分析。

    2.3 差異表達(dá)基因識(shí)別與分析

    低溫處理下,基因的表達(dá)量會(huì)根據(jù)自身對(duì)低溫的敏感程度發(fā)生相應(yīng)變化,其差異表達(dá)程度可反映對(duì)低溫的敏感度。因此,在數(shù)據(jù)質(zhì)控分析之后,進(jìn)一步對(duì)差異表達(dá)基因(differentially expressed gene,DEG)進(jìn)行統(tǒng)計(jì)分析。由圖2-a可知,以|log2 FC|≥1,F(xiàn)DR≤0.01為篩選閾值,在龍粳11中,低溫處理2 d后共有9 468個(gè)DEG,其中包括上調(diào)表達(dá)基因 4 915 個(gè)、下調(diào)表達(dá)基因4 553個(gè)。而在龍粳25中,低溫處理2 d后共有6 784個(gè)基因的表達(dá)量發(fā)生變化,其中包括上調(diào)表達(dá)基因3 499個(gè)、下調(diào)表達(dá)基因 3 285 個(gè)(圖2-b)。進(jìn)一步通過韋恩圖(圖2-c)分析可以看出,龍粳25、龍粳11共有的DEG數(shù)目為3 508個(gè),而龍粳25、龍粳11中特有的DEG分別為3 276、5 960個(gè)。結(jié)果表明,低溫處理2 d后,龍粳11的DEG數(shù)目明顯多于龍粳25,表明龍粳11對(duì)低溫更為敏感,低溫反應(yīng)更加明顯。

    2.4 差異表達(dá)基因的GO富集分析

    通過分析DEG的韋恩圖,發(fā)現(xiàn)有3 508個(gè)DEG是龍粳25、龍粳11在低溫處理后共有的,而在龍粳25中有3 276個(gè)DEG是其低溫下特有的。為近一步明確其功能分類,對(duì)這些基因分別展開GO富集分析。通過對(duì)二者共有的DEG進(jìn)行GO分析,以KS≤0.01為篩選標(biāo)準(zhǔn),生物進(jìn)程、細(xì)胞組分、分子功能的GO功能類別明顯富集到的類別分為95、20、41條。圖3-a分別列出每個(gè)類別中含有基因數(shù)最多的前20條terms。生物進(jìn)程方面,主要富集在生物進(jìn)程調(diào)控(GO:0050789,631個(gè)DEG)、氧化還原反應(yīng)(GO:0055114,301個(gè)DEG)、非生物刺激應(yīng)答(GO:0009628,302個(gè)DEG)、蛋白磷酸化(GO:0006468,279個(gè)DEG)等類別中。細(xì)胞組分方面,主要富集在膜(GO:0016020,967個(gè)DEG)、細(xì)胞質(zhì)膜包圍的囊泡(GO:0016023,773個(gè)DEG)、原生質(zhì)膜(GO:0005886,427個(gè)DEG)、細(xì)胞核(GO:0005634,382個(gè)DEG)等類別中。在分子功能方面,主要富集在ATP結(jié)合功能(GO:0005524,439個(gè)DEG)、絲氨酸/蘇氨酸蛋白激酶活性(GO:0004674,211個(gè)DEG)、轉(zhuǎn)運(yùn)活性(GO:0005215,189個(gè)DEG)、序列特異性DNA結(jié)合轉(zhuǎn)錄(GO:0003700,156個(gè)DEG)等類別中。

    由于龍粳25具有孕穗期強(qiáng)耐冷性,為了發(fā)現(xiàn)其低溫下特有的功能類別,對(duì)其特有的DEG進(jìn)行GO富集分析。同樣以KS≤0.01為篩選標(biāo)準(zhǔn),明顯富集到的生物進(jìn)程、細(xì)胞組分、分子功能的GO功能類別分別為67、24、32條。由圖3-b可見,在生物進(jìn)程方面,龍粳25特有基因多富集在生物進(jìn)程調(diào)控(GO:0019222,226個(gè)DEG)、非生物刺激應(yīng)答(GO:0009628,222個(gè)DEG)、氧化還原反應(yīng)(GO:0055114,215個(gè)DEG)、碳水化合物代謝過程(GO:0005975,176個(gè)DEG)等類別中。在細(xì)胞組分方面,龍粳25特有基因富集最多的前4條類別與龍粳25、龍粳11共有基因富集的類別相同,而所包含的基因數(shù)則分別為723、613、324、289個(gè)。在分子功能方面,龍粳25特有基因主要富集在絲氨酸/蘇氨酸蛋白激酶(GO:0004674,158個(gè)DEG)、序列特異性DNA結(jié)合轉(zhuǎn)錄(GO:0003700,90個(gè)DEG)、轉(zhuǎn)移酶活性/轉(zhuǎn)移糖基(GO:0016757,81個(gè)DEG)、血紅素結(jié)合部位(GO:0020037,69個(gè)DEG)等類別中。

    2.5 差異表達(dá)基因的KEGG富集分析

    為近一步探明低溫對(duì)寒地水稻代謝途徑的影響,更加明確地闡釋DEG的生物學(xué)功能,對(duì)龍粳25、龍粳11低溫處理下共有的DEG以及龍粳25特有的DEG分別進(jìn)行KEGG富集分析。結(jié)果表明,二者共有的3 508個(gè)DEG明顯富集到112個(gè)KEGG通路中。通過列出富集程度最高的前20個(gè)通路(圖4-a)可以看出,DEG主要富集在植物激素信號(hào)轉(zhuǎn)導(dǎo)(ko04075,46個(gè)DEG,富集因子1.81)、碳代謝(ko01200,44個(gè)DEG,富集因子1.61)、淀粉蔗糖代謝(ko00500,43個(gè)DEG,富集因子1.86)、苯丙烷合成(ko00940,35個(gè)DEG,富集因子1.61)、糖酵解/糖原異生通路(ko00010,33個(gè)DEG,富集因子2.67)等途徑。而在龍粳25特有的DEG中,共富集到110個(gè)KEGG通路中。同樣列出富集程度最高的前20個(gè)通路,DEG多富集在植物激素與信號(hào)轉(zhuǎn)導(dǎo)(ko04075,41個(gè)DEG,富集因子1.94)、內(nèi)質(zhì)網(wǎng)蛋白加工(ko04141,34個(gè)DEG,富集因子1.61)、苯丙烷合成(ko00940,32個(gè)DEG,富集因子1.76)、氨基糖和核酸糖代謝(ko00520,21個(gè)DEG,富集因子1.87)、谷胱甘肽代謝(ko00480,17個(gè)DEG,富集因子2.10)等通路中(圖4-b)。綜合比較2個(gè)部分KEGG通路,發(fā)現(xiàn)植物激素信號(hào)轉(zhuǎn)導(dǎo)(ko04075)在龍粳25、龍粳11共有DEG及龍粳25特有DEG中都被明顯富集到含DEG均為最多的通路,表明這條通路在寒地水稻低溫脅迫應(yīng)答中起著重要的作用。

    2.6 差異表達(dá)基因的qRT-PCR驗(yàn)證

    為驗(yàn)證轉(zhuǎn)錄組測(cè)序結(jié)果的準(zhǔn)確性,從測(cè)序結(jié)果中隨機(jī)選取6個(gè)表達(dá)差異大的基因進(jìn)行qRT-PCR分析。分別以龍粳11未做低溫處理時(shí)的各個(gè)基因表達(dá)量作為對(duì)照(相對(duì)表達(dá)量為1),qRT-PCR結(jié)果表明,經(jīng)過低溫處理后2個(gè)品種的6個(gè)基因表達(dá)均產(chǎn)生差異變化;與轉(zhuǎn)錄組測(cè)序結(jié)果相比,雖然基因表達(dá)量的變幅與轉(zhuǎn)錄組結(jié)果不盡一致,但變化趨勢(shì)與轉(zhuǎn)錄測(cè)序結(jié)果均一致,證實(shí)了轉(zhuǎn)錄組測(cè)序結(jié)果的準(zhǔn)確性及可靠性(圖5)。

    3 結(jié)論與討論

    本研究通過RNA-seq技術(shù)對(duì)寒地水稻低溫下展開轉(zhuǎn)錄組測(cè)序,比較2個(gè)品種低溫脅迫后的差異表達(dá)基因。研究結(jié)果表明,低溫脅迫后,龍粳11的DEG為9 468個(gè),龍粳25為6 784個(gè);低溫處理后,龍粳11的DEG明顯多于龍粳25,表明龍粳11對(duì)低溫更加敏感,而龍粳25對(duì)低溫則表現(xiàn)出更強(qiáng)的耐受性,這也可以從二者在低溫下的結(jié)實(shí)率得到驗(yàn)證。GO富集分析發(fā)現(xiàn),低溫脅迫下,生物進(jìn)程調(diào)控、氧化還原反應(yīng)、非生物刺激應(yīng)答、蛋白磷酸化等GO類別的反應(yīng)比較大。KEGG富集分析得出,檸檬烯和蒎烯的降解、糖酵解/糖異生、半乳糖代謝、(角質(zhì)、木質(zhì)、蠟質(zhì))生物合成、植物激素信號(hào)轉(zhuǎn)導(dǎo)等途徑在低溫脅迫后富集的程度最高。以上分析表明,低溫對(duì)寒地水稻的影響涉及到生長(zhǎng)發(fā)育的各個(gè)方面,包括各類復(fù)雜的新陳代謝變化。

    本研究所用的材料為孕穗期耐冷的龍粳25和不耐冷的龍粳11,在探究寒地水稻特有的低溫應(yīng)答防御機(jī)制及相關(guān)基因時(shí),龍粳25特有的DEG及其GO類別值得進(jìn)一步深入研究。本研究中,龍粳25富集到67條生物進(jìn)程、24條細(xì)胞組分、32條分子功能類別。由于生物進(jìn)程能夠直接反映生物體生長(zhǎng)發(fā)育情況,對(duì)此類別展開分析。氧化還原反應(yīng)進(jìn)程不但在二者共有的DEG中得到明顯富集,并且在龍粳25特有的DEG中也被明顯富集,表明此GO類別在寒地水稻低溫應(yīng)答中起重要作用。

    植物體內(nèi)源激素雖然含量很低,但是在生物脅迫、非生物脅迫等眾多脅迫因子方面起到積極的防衛(wèi)作用[19。目前研究較多的激素包括脫落酸(ABA)、赤霉素(GA)、生長(zhǎng)素(IAA)、乙烯(ETH)、水楊酸(SA)、茉莉酸(JA)、細(xì)胞分裂素(CK)、油菜素內(nèi)酯(BR)等[20。而IAA、ABA、ETH等激素在水稻低溫應(yīng)答中起積極作用[21-22。本研究中,植物激素信號(hào)轉(zhuǎn)導(dǎo)通路在龍粳11、龍粳25共有的DEG以及龍粳25特有的DEG都被明顯富集,且含有基因數(shù)也是最多的,表明此通路在寒地水稻低溫應(yīng)答中同樣起關(guān)鍵調(diào)控作用,這一結(jié)果也與前人研究結(jié)果[23一致。

    低溫脅迫會(huì)對(duì)水稻植株的正常生長(zhǎng)產(chǎn)生各種生理及分子調(diào)控機(jī)制方面的影響;水稻通過自身調(diào)整,盡可能適應(yīng)溫度脅迫,以保持其正常生長(zhǎng)[24-25。黑龍江地處寒稻作區(qū),由于其特殊的地理環(huán)境及生態(tài)類型,低溫是寒地水稻正常生長(zhǎng)的重要限制因子[26。本研究以寒地水稻孕穗期耐冷差異明顯的品種龍粳25、龍粳11為研究材料,在前人對(duì)其結(jié)實(shí)率研究的基礎(chǔ)上,對(duì)低溫脅迫下穎花進(jìn)行轉(zhuǎn)錄組測(cè)序分析。通過差異表達(dá)基因分析,發(fā)現(xiàn)龍粳11的DEG數(shù)量明顯多于龍粳25,表明龍粳11對(duì)低溫反應(yīng)更為敏感。GO、KEGG分析表明,氧化還原反應(yīng)進(jìn)程、低溫應(yīng)答、植物激素信號(hào)轉(zhuǎn)導(dǎo)等眾多GO類別及KEGG通路均參與到低溫脅迫應(yīng)答中。此外,通過qRT-PCR驗(yàn)證,驗(yàn)證了RNA-seq結(jié)果的準(zhǔn)確性,期待本研究結(jié)果對(duì)寒地水稻孕穗期耐冷的分子遺傳調(diào)控機(jī)理的解析有一定的推動(dòng)作用。

    參考文獻(xiàn):

    [1]程式華,李 建. 現(xiàn)代中國(guó)水稻[M]. 北京:金盾出版社,2007.

    [2]Howarth C J,Ougham H J. Gene expression under temperature stress[J]. The New Phytologist,1993,125(1):1-26.

    [3]Guo Z H,Cai L J,Chen Z Q,et al. Identification of candidate genes controlling chilling tolerance of rice in the cold region at the booting stage by BSA-Seq and RNA-Seq[J]. Royal Society Open Science,2020,7(11):201081.

    [4]Guo Z H,Liu C X,Xiao W M,et al. Comparative transcriptome profile analysis of anther development in reproductive stage of rice in cold region under cold stress[J]. Plant Molecular Biology Reporter,2019,37(3):129-145.

    [5]Mamun E A,Alfred S,Cantrill L C,et al. Effects of chilling on male gametophyte development in rice[J]. Cell Biology International,2006,30(7):583-591.

    [6]Oda S,Kaneko F,Yano K,et al. Morphological and gene expression analysis under cool temperature conditions in rice anther development[J]. Genes amp; Genetic Systems,2010,85(2):107-120.

    [7]Qiu C P,Ethier G,Pepin S,et al. Persistent negative temperature response of mesophyll conductance in red raspberry (Rubus idaeus L.) leaves under both high and low vapour pressure deficits:a role for abscisic acid?[J]. Plant,Cell amp; Environment,2018,41(4):876.

    [8]李健陵,霍治國(guó),吳麗姬,等. 孕穗期低溫對(duì)水稻產(chǎn)量的影響及其生理機(jī)制[J]. 中國(guó)水稻科學(xué),2014,28(3):277-288.

    [9]Saito K,Hayano-Saito Y,Kuroki M,et al. Map-based cloning of the rice cold tolerance gene Ctb1[J]. Plant Science,2010,179(1/2):97-102.

    [10]Vannini C,Campa M,Iriti M,et al. Evaluation of transgenic tomato plants ectopically expressing the rice Osmyb4 gene[J]. Plant Science,2007,173(2):231-239.

    [11]Kim S H,Choi H S,Cho Y C,et al. Cold-responsive regulation of a flower-preferential class Ⅲ peroxidase gene,OsPOX1,in rice (Oryza sativa L.)[J]. Journal of Plant Biology,2012,55(2):123-131.

    [12]Sato Y,Masuta Y,Saito K,et al. Enhanced chilling tolerance at the booting stage in rice by transgenic overexpression of the ascorbate peroxidase gene,OsAPXa[J]. Plant Cell Reports,2011,30(3):399-406.

    [13]Guo Z H,Cai L J,Liu C X,et al. Global analysis of differentially expressed genes between two Japonica rice varieties induced by low temperature during the booting stage by RNA-Seq[J]. Royal Society Open Science,2020,7(6):192243.

    [14]韋秀葉,趙信林,郭 媛,等. 麻類作物轉(zhuǎn)錄組測(cè)序分析研究進(jìn)展[J]. 中國(guó)麻業(yè)科學(xué),2020,42(3):128-134.

    [15]吳 倩,張 磊,黃志平,等. 轉(zhuǎn)錄組測(cè)序及其在野生大豆基因資源發(fā)掘中的應(yīng)用[J]. 大豆科學(xué),2013,32(6):845-851.

    [16]許 波,張偉強(qiáng),馮曉曦,等. 轉(zhuǎn)錄組測(cè)序技術(shù)在玉米中應(yīng)用研究進(jìn)展[J]. 玉米科學(xué),2014,22(1):67-72,78.

    [17]Bai B,Wu J,Sheng W T,et al. Comparative analysis of anther transcriptome profiles of two different rice male sterile lines genotypes under cold stress[J]. International Journal of Molecular Sciences,2015,16(5):11398-11416.

    [18]Zhang L,Li X Y,Li D K,et al. CARK1 mediates ABA signaling by phosphorylation of ABA receptors[J]. Cell Discovery,2018,4:30.

    [19]叢 慶,張 琪,宋麗莉,等. 激素在植物冷脅迫應(yīng)答中的角色[J]. 核農(nóng)學(xué)報(bào),2016,30(3):614-619.

    [20]張明菊,朱 莉,夏啟中. 植物激素對(duì)脅迫反應(yīng)調(diào)控的研究進(jìn)展[J]. 湖北大學(xué)學(xué)報(bào)(自然科學(xué)版),2021,43(3):242-253,263.

    [21]劉 芬. 低溫脅迫對(duì)細(xì)枝木麻黃無性系生理指標(biāo)和轉(zhuǎn)錄組的影響[D]. 長(zhǎng)沙:中南林業(yè)科技大學(xué),2015:5-7.

    [22]Jain M,Khurana J P. Transcript profiling reveals diverse roles of auxin-responsive genes during reproductive development and abiotic stress in rice[J]. The FEBS Journal,2009,276(11):3148-3162.

    [23]高紅秀,朱 琳,劉天奇,等. 水稻植物激素響應(yīng)低溫脅迫反應(yīng)的轉(zhuǎn)錄組分析[J]. 分子植物育種,2021,19(13):4188-4197.

    [24]胡瀟婕,毛東海. 基于RNA-Seq技術(shù)分析植物激素信號(hào)途徑在水稻幼苗中對(duì)低溫脅迫的應(yīng)答規(guī)律[J]. 農(nóng)業(yè)現(xiàn)代化研究,2019,40(5):878-890.

    [25]劉琳帥,卞景陽,孫興榮,等. 水稻低溫冷害的研究進(jìn)展[J]. 江蘇農(nóng)業(yè)科學(xué),2022,50(24):9-15.

    [26]賀 梅,宋冬明. 寒地稻作區(qū)冷害防御措施[J]. 北方水稻,2015,45(1):43-44.

    猜你喜歡
    轉(zhuǎn)錄組寒地水稻
    什么是海水稻
    基于寒地公園適老化需求的研究——以長(zhǎng)春為例
    園林科技(2021年2期)2022-01-19 02:22:22
    水稻種植60天就能收獲啦
    軍事文摘(2021年22期)2021-11-26 00:43:51
    寒地黑土無公害水產(chǎn)品健康養(yǎng)殖發(fā)展思路
    寒地鴨場(chǎng)飼養(yǎng)管理與疾病防治
    一季水稻
    文苑(2020年6期)2020-06-22 08:41:52
    水稻花
    文苑(2019年22期)2019-12-07 05:29:00
    寒地多年生毛櫻桃改接李優(yōu)良品種技術(shù)
    基于轉(zhuǎn)錄組測(cè)序的山茱萸次生代謝生物合成相關(guān)基因的挖掘
    金釵石斛轉(zhuǎn)錄組SSR位點(diǎn)信息分析
    国产欧美日韩一区二区三区在线| av一本久久久久| 18在线观看网站| 欧美日韩视频精品一区| 一级爰片在线观看| 免费av中文字幕在线| 欧美日韩一区二区视频在线观看视频在线| 黑人巨大精品欧美一区二区蜜桃 | 亚洲国产精品成人久久小说| 亚洲av福利一区| 天天操日日干夜夜撸| av在线老鸭窝| 啦啦啦中文免费视频观看日本| 亚洲精品乱码久久久久久按摩| 赤兔流量卡办理| 国产精品99久久99久久久不卡 | 美女视频免费永久观看网站| 成人18禁高潮啪啪吃奶动态图| 中文天堂在线官网| 国产精品麻豆人妻色哟哟久久| 国产精品熟女久久久久浪| 草草在线视频免费看| 国产伦理片在线播放av一区| 国产精品国产三级国产专区5o| 成人18禁高潮啪啪吃奶动态图| 波野结衣二区三区在线| 飞空精品影院首页| 久久国产亚洲av麻豆专区| 在线观看国产h片| 多毛熟女@视频| 中文字幕另类日韩欧美亚洲嫩草| 免费观看av网站的网址| 免费高清在线观看日韩| 精品酒店卫生间| 在线精品无人区一区二区三| 国产一区二区激情短视频 | 久久精品国产自在天天线| 22中文网久久字幕| 一级毛片黄色毛片免费观看视频| 90打野战视频偷拍视频| 男女无遮挡免费网站观看| 在线亚洲精品国产二区图片欧美| 丝袜脚勾引网站| 国产免费一区二区三区四区乱码| 亚洲国产精品一区三区| 777米奇影视久久| 在线观看免费日韩欧美大片| 久久99一区二区三区| 色94色欧美一区二区| 熟女av电影| 国产一区有黄有色的免费视频| 女的被弄到高潮叫床怎么办| 国产黄色视频一区二区在线观看| 国产精品久久久av美女十八| 日本爱情动作片www.在线观看| 大话2 男鬼变身卡| 日本欧美视频一区| 韩国av在线不卡| 日韩免费高清中文字幕av| 免费观看性生交大片5| 亚洲欧洲日产国产| 国产麻豆69| 国产一区二区三区av在线| a级毛片在线看网站| 99视频精品全部免费 在线| 精品人妻熟女毛片av久久网站| 午夜免费观看性视频| 精品亚洲成a人片在线观看| 婷婷成人精品国产| 亚洲五月色婷婷综合| 老司机影院成人| 精品少妇内射三级| 黄片播放在线免费| 国产激情久久老熟女| 一本色道久久久久久精品综合| 最黄视频免费看| 丝袜人妻中文字幕| 久久久久国产精品人妻一区二区| 成人漫画全彩无遮挡| 高清av免费在线| 亚洲精品色激情综合| 日韩不卡一区二区三区视频在线| 制服诱惑二区| 超色免费av| 国产精品 国内视频| 日本免费在线观看一区| 亚洲国产日韩一区二区| 欧美激情 高清一区二区三区| 国产成人精品一,二区| 高清在线视频一区二区三区| 秋霞在线观看毛片| a级毛片黄视频| 欧美国产精品va在线观看不卡| 精品少妇内射三级| 日本爱情动作片www.在线观看| 成人18禁高潮啪啪吃奶动态图| 欧美国产精品一级二级三级| xxx大片免费视频| 菩萨蛮人人尽说江南好唐韦庄| 日韩伦理黄色片| av播播在线观看一区| 婷婷色麻豆天堂久久| 久久免费观看电影| 永久网站在线| 香蕉丝袜av| 色网站视频免费| 国产免费现黄频在线看| 国产在线视频一区二区| 午夜福利,免费看| 国产精品不卡视频一区二区| 国产av码专区亚洲av| 美女脱内裤让男人舔精品视频| 99re6热这里在线精品视频| 欧美激情 高清一区二区三区| 18+在线观看网站| 久热这里只有精品99| 日韩欧美一区视频在线观看| 日韩电影二区| 午夜影院在线不卡| av在线老鸭窝| 夜夜爽夜夜爽视频| 久久国产精品男人的天堂亚洲 | 亚洲av国产av综合av卡| 亚洲av国产av综合av卡| 美女大奶头黄色视频| 一级片免费观看大全| 日韩精品有码人妻一区| 大话2 男鬼变身卡| 夫妻性生交免费视频一级片| 欧美人与善性xxx| 在线观看一区二区三区激情| 黄色 视频免费看| 国产黄色视频一区二区在线观看| 婷婷色综合www| 久久国产亚洲av麻豆专区| 天堂8中文在线网| 免费av中文字幕在线| 各种免费的搞黄视频| 美女内射精品一级片tv| 国产亚洲最大av| av又黄又爽大尺度在线免费看| 最黄视频免费看| 国产精品无大码| 男女边摸边吃奶| 国产免费福利视频在线观看| 免费高清在线观看日韩| 国产成人精品福利久久| 久热这里只有精品99| 高清在线视频一区二区三区| 如何舔出高潮| 国产69精品久久久久777片| 成年人午夜在线观看视频| 最后的刺客免费高清国语| 久久97久久精品| 日韩成人av中文字幕在线观看| 国产视频首页在线观看| 18禁裸乳无遮挡动漫免费视频| 26uuu在线亚洲综合色| 街头女战士在线观看网站| 美女国产视频在线观看| 18禁裸乳无遮挡动漫免费视频| 亚洲一级一片aⅴ在线观看| 91久久精品国产一区二区三区| 一本久久精品| 久久精品夜色国产| 自线自在国产av| 亚洲三级黄色毛片| 久久毛片免费看一区二区三区| 九草在线视频观看| 日本免费在线观看一区| 两性夫妻黄色片 | 国产激情久久老熟女| 国产片内射在线| 亚洲精品久久午夜乱码| 久久精品国产综合久久久 | 亚洲av电影在线进入| 日本-黄色视频高清免费观看| 纵有疾风起免费观看全集完整版| 蜜臀久久99精品久久宅男| 日韩伦理黄色片| 精品久久蜜臀av无| 伦理电影免费视频| 国精品久久久久久国模美| 日韩免费高清中文字幕av| 亚洲人与动物交配视频| 妹子高潮喷水视频| 午夜精品国产一区二区电影| 亚洲国产精品国产精品| 欧美97在线视频| 久久av网站| 90打野战视频偷拍视频| 男女免费视频国产| 亚洲激情五月婷婷啪啪| 国产 一区精品| 国产一区二区在线观看av| 国产欧美另类精品又又久久亚洲欧美| 国产激情久久老熟女| 综合色丁香网| 亚洲成人手机| 曰老女人黄片| 久久精品久久久久久久性| 欧美成人午夜免费资源| 日本欧美视频一区| 国产毛片在线视频| 看十八女毛片水多多多| 亚洲欧美一区二区三区黑人 | 99视频精品全部免费 在线| 寂寞人妻少妇视频99o| 精品国产乱码久久久久久小说| 成人18禁高潮啪啪吃奶动态图| 久久人人爽av亚洲精品天堂| 一二三四中文在线观看免费高清| 99久国产av精品国产电影| 99热全是精品| 乱人伦中国视频| 午夜福利乱码中文字幕| 国产一区二区三区综合在线观看 | 天美传媒精品一区二区| 亚洲精品国产av成人精品| 久久久久国产精品人妻一区二区| 欧美 日韩 精品 国产| 亚洲欧美色中文字幕在线| 色5月婷婷丁香| 熟女电影av网| 一区二区三区精品91| 日本vs欧美在线观看视频| 免费av中文字幕在线| 国产精品 国内视频| 如日韩欧美国产精品一区二区三区| 国产成人91sexporn| 欧美精品av麻豆av| 精品国产国语对白av| 97精品久久久久久久久久精品| 欧美成人午夜免费资源| 天天影视国产精品| 日本黄色日本黄色录像| 久久99一区二区三区| 观看av在线不卡| 国产激情久久老熟女| av卡一久久| 久久久久久人人人人人| 18禁观看日本| 午夜免费观看性视频| 久久久国产精品麻豆| 久久久亚洲精品成人影院| 99香蕉大伊视频| 大码成人一级视频| av天堂久久9| 中文字幕精品免费在线观看视频 | 高清黄色对白视频在线免费看| 日本午夜av视频| 自拍欧美九色日韩亚洲蝌蚪91| 综合色丁香网| 中文乱码字字幕精品一区二区三区| 久久久久久人人人人人| 一区二区日韩欧美中文字幕 | av免费观看日本| 日韩av免费高清视频| 男男h啪啪无遮挡| 99热全是精品| av在线老鸭窝| 精品人妻熟女毛片av久久网站| 永久免费av网站大全| 欧美日韩视频精品一区| 亚洲综合色网址| 亚洲精品av麻豆狂野| 亚洲精品成人av观看孕妇| 97在线人人人人妻| 久久久久久久久久成人| 国产福利在线免费观看视频| 99香蕉大伊视频| 狂野欧美激情性xxxx在线观看| 美国免费a级毛片| 蜜桃国产av成人99| 久热这里只有精品99| 最近中文字幕2019免费版| 亚洲熟女精品中文字幕| 欧美另类一区| 国产亚洲精品久久久com| 亚洲精品国产av蜜桃| 亚洲欧美成人精品一区二区| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品色激情综合| 国产成人a∨麻豆精品| a级毛片黄视频| 精品99又大又爽又粗少妇毛片| 成年女人在线观看亚洲视频| 看免费成人av毛片| 国产成人精品久久久久久| 有码 亚洲区| 久久人人爽av亚洲精品天堂| 最后的刺客免费高清国语| 51国产日韩欧美| 99久久人妻综合| 宅男免费午夜| 你懂的网址亚洲精品在线观看| 91国产中文字幕| 国产精品人妻久久久影院| 国产成人av激情在线播放| 女性被躁到高潮视频| 欧美3d第一页| 久久久久久久亚洲中文字幕| 日韩视频在线欧美| 尾随美女入室| 一级片免费观看大全| 看免费av毛片| 麻豆精品久久久久久蜜桃| 十八禁网站网址无遮挡| 成人手机av| 王馨瑶露胸无遮挡在线观看| 亚洲国产精品国产精品| 国产免费现黄频在线看| 免费播放大片免费观看视频在线观看| 国产免费视频播放在线视频| 伊人久久国产一区二区| 久久精品国产亚洲av涩爱| 国产一区二区在线观看日韩| 亚洲,一卡二卡三卡| 亚洲国产欧美日韩在线播放| 亚洲,欧美,日韩| 插逼视频在线观看| av片东京热男人的天堂| 国产在视频线精品| 国产日韩欧美亚洲二区| 成人亚洲精品一区在线观看| 欧美日韩精品成人综合77777| 人人澡人人妻人| 丰满饥渴人妻一区二区三| 宅男免费午夜| 丝袜人妻中文字幕| 亚洲精品国产色婷婷电影| 国产成人aa在线观看| 在线观看免费视频网站a站| 午夜福利网站1000一区二区三区| 国产精品人妻久久久影院| av不卡在线播放| 在线观看一区二区三区激情| 欧美日韩一区二区视频在线观看视频在线| 99热网站在线观看| 性色avwww在线观看| 日本91视频免费播放| 免费观看av网站的网址| 日韩大片免费观看网站| 国语对白做爰xxxⅹ性视频网站| 午夜免费鲁丝| 欧美丝袜亚洲另类| 高清视频免费观看一区二区| 亚洲精品,欧美精品| 亚洲av电影在线观看一区二区三区| 免费黄色在线免费观看| 国产精品国产av在线观看| 亚洲激情五月婷婷啪啪| 男女高潮啪啪啪动态图| 亚洲精品第二区| 日韩在线高清观看一区二区三区| 亚洲一区二区三区欧美精品| 激情五月婷婷亚洲| 国产乱来视频区| 国产成人一区二区在线| 欧美精品一区二区免费开放| 亚洲国产精品一区三区| 丰满饥渴人妻一区二区三| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一区二区三区四区激情视频| 纯流量卡能插随身wifi吗| 欧美日本中文国产一区发布| 国产欧美日韩综合在线一区二区| 亚洲 欧美一区二区三区| 午夜精品国产一区二区电影| 国产一区二区激情短视频 | 新久久久久国产一级毛片| 肉色欧美久久久久久久蜜桃| 亚洲av福利一区| 水蜜桃什么品种好| 亚洲精华国产精华液的使用体验| 国产精品成人在线| 亚洲欧美日韩卡通动漫| 亚洲av综合色区一区| 欧美激情极品国产一区二区三区 | 午夜福利视频精品| 男人添女人高潮全过程视频| 在线天堂最新版资源| av又黄又爽大尺度在线免费看| 久久久久网色| 亚洲精品一二三| 观看av在线不卡| videossex国产| 青春草视频在线免费观看| 纵有疾风起免费观看全集完整版| 一级毛片电影观看| 99热这里只有是精品在线观看| 亚洲国产日韩一区二区| 我的女老师完整版在线观看| 国产精品嫩草影院av在线观看| 午夜福利乱码中文字幕| av在线观看视频网站免费| 啦啦啦中文免费视频观看日本| 内地一区二区视频在线| 欧美日韩综合久久久久久| 一级片'在线观看视频| 久久久精品区二区三区| 精品酒店卫生间| 久久影院123| 国产免费又黄又爽又色| 欧美日韩视频精品一区| 亚洲国产精品成人久久小说| 国产av一区二区精品久久| 国产精品偷伦视频观看了| 卡戴珊不雅视频在线播放| av在线老鸭窝| 视频在线观看一区二区三区| 免费观看性生交大片5| 国产乱人偷精品视频| 男的添女的下面高潮视频| 人人澡人人妻人| 亚洲av在线观看美女高潮| 天天操日日干夜夜撸| 中国国产av一级| 久久久久久人妻| 天天操日日干夜夜撸| 综合色丁香网| 国产亚洲av片在线观看秒播厂| 日韩一本色道免费dvd| 校园人妻丝袜中文字幕| 久久av网站| 国产精品免费大片| 黄网站色视频无遮挡免费观看| 亚洲成av片中文字幕在线观看 | 国产精品久久久久久精品电影小说| 午夜免费观看性视频| 亚洲精品久久久久久婷婷小说| 久久久久人妻精品一区果冻| 黄色 视频免费看| 国产白丝娇喘喷水9色精品| 日韩制服骚丝袜av| 国内精品宾馆在线| 少妇精品久久久久久久| 欧美精品人与动牲交sv欧美| 国产亚洲最大av| 中文字幕精品免费在线观看视频 | 最新中文字幕久久久久| 最近手机中文字幕大全| 久久国产精品男人的天堂亚洲 | 欧美日韩国产mv在线观看视频| 丁香六月天网| a级毛片在线看网站| 视频中文字幕在线观看| 韩国精品一区二区三区 | 日日摸夜夜添夜夜爱| 91成人精品电影| 日韩人妻精品一区2区三区| 青春草亚洲视频在线观看| 欧美 日韩 精品 国产| 九九在线视频观看精品| 中文精品一卡2卡3卡4更新| 免费黄色在线免费观看| 在线天堂中文资源库| 中国三级夫妇交换| 97精品久久久久久久久久精品| 捣出白浆h1v1| 99久久中文字幕三级久久日本| 国产精品一区www在线观看| 亚洲成av片中文字幕在线观看 | 免费av中文字幕在线| 一边亲一边摸免费视频| 丝袜在线中文字幕| 纵有疾风起免费观看全集完整版| 亚洲精品久久午夜乱码| 精品视频人人做人人爽| 久久久久精品久久久久真实原创| 老司机亚洲免费影院| 成年人免费黄色播放视频| 久久久国产一区二区| 晚上一个人看的免费电影| 国产女主播在线喷水免费视频网站| av国产精品久久久久影院| 秋霞伦理黄片| 亚洲国产av影院在线观看| 99热全是精品| 精品熟女少妇av免费看| 亚洲精品久久午夜乱码| 亚洲三级黄色毛片| 日韩精品有码人妻一区| 在线免费观看不下载黄p国产| 色94色欧美一区二区| 久久影院123| 国产精品麻豆人妻色哟哟久久| 亚洲天堂av无毛| 97人妻天天添夜夜摸| 性高湖久久久久久久久免费观看| www.熟女人妻精品国产 | 亚洲熟女精品中文字幕| 久久久久精品性色| 久久久国产精品麻豆| 久久久久人妻精品一区果冻| 免费黄频网站在线观看国产| 欧美精品一区二区免费开放| 久久久久网色| 大香蕉久久网| 久久久久人妻精品一区果冻| tube8黄色片| 69精品国产乱码久久久| 国产成人aa在线观看| 99视频精品全部免费 在线| 美女xxoo啪啪120秒动态图| 久久国产精品男人的天堂亚洲 | 天堂俺去俺来也www色官网| 在线看a的网站| 嫩草影院入口| 免费在线观看完整版高清| 啦啦啦在线观看免费高清www| 中文字幕免费在线视频6| 汤姆久久久久久久影院中文字幕| 日本欧美国产在线视频| 女性生殖器流出的白浆| 欧美性感艳星| 人妻 亚洲 视频| 亚洲精品乱久久久久久| 久久久久久人妻| 久久99蜜桃精品久久| 黄色视频在线播放观看不卡| 亚洲成人av在线免费| 久久精品久久精品一区二区三区| 午夜福利乱码中文字幕| 成人国语在线视频| 黑人巨大精品欧美一区二区蜜桃 | 欧美精品亚洲一区二区| 欧美日韩一区二区视频在线观看视频在线| 亚洲 欧美一区二区三区| a级片在线免费高清观看视频| 久久人人爽人人爽人人片va| 国产深夜福利视频在线观看| 国产一区二区在线观看av| 久热久热在线精品观看| 日本与韩国留学比较| 男人爽女人下面视频在线观看| 亚洲人成77777在线视频| 亚洲,一卡二卡三卡| 亚洲,欧美,日韩| 成人影院久久| 久热这里只有精品99| 亚洲av福利一区| 男女午夜视频在线观看 | 黄网站色视频无遮挡免费观看| 伊人亚洲综合成人网| 色婷婷久久久亚洲欧美| 大陆偷拍与自拍| 亚洲高清免费不卡视频| www.av在线官网国产| 色视频在线一区二区三区| 欧美老熟妇乱子伦牲交| 九色成人免费人妻av| 免费少妇av软件| 爱豆传媒免费全集在线观看| 精品人妻一区二区三区麻豆| 一区二区av电影网| 狠狠精品人妻久久久久久综合| 波多野结衣一区麻豆| 超碰97精品在线观看| 男的添女的下面高潮视频| 国产在视频线精品| av黄色大香蕉| 欧美bdsm另类| 蜜桃国产av成人99| 国产精品无大码| 另类精品久久| 制服诱惑二区| 国产精品久久久久久精品古装| 大香蕉久久成人网| 黄网站色视频无遮挡免费观看| 亚洲激情五月婷婷啪啪| 熟女人妻精品中文字幕| 最近最新中文字幕大全免费视频 | 亚洲精品av麻豆狂野| 99re6热这里在线精品视频| 美女xxoo啪啪120秒动态图| 尾随美女入室| 久久久久人妻精品一区果冻| 久热这里只有精品99| av一本久久久久| av视频免费观看在线观看| 午夜久久久在线观看| 精品少妇黑人巨大在线播放| 日韩制服骚丝袜av| 亚洲国产成人一精品久久久| 亚洲精品,欧美精品| 热re99久久精品国产66热6| 国产精品无大码| 美女视频免费永久观看网站| 69精品国产乱码久久久| 中文字幕制服av| 亚洲天堂av无毛| 如何舔出高潮| 免费播放大片免费观看视频在线观看| 男人舔女人的私密视频| 亚洲国产av影院在线观看| 新久久久久国产一级毛片| 各种免费的搞黄视频| 精品一区二区三卡| 多毛熟女@视频| 97超碰精品成人国产| 国产综合精华液| www.av在线官网国产| 亚洲高清免费不卡视频| 国产熟女欧美一区二区| 老司机影院成人| 91午夜精品亚洲一区二区三区| 国产69精品久久久久777片| 免费观看a级毛片全部| 中文字幕亚洲精品专区| 精品一区二区三区四区五区乱码 | 国产1区2区3区精品| 久久国产精品大桥未久av| 久久精品国产综合久久久 | 美女内射精品一级片tv|