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

    低溫脅迫下海濱木槿的轉(zhuǎn)錄組分析

    2023-06-30 12:00:49趙文靜鄭旭葛金濤劉興滿(mǎn)繆美華孔祥偉
    江蘇農(nóng)業(yè)科學(xué) 2023年5期
    關(guān)鍵詞:轉(zhuǎn)錄組低溫脅迫差異基因

    趙文靜 鄭旭 葛金濤 劉興滿(mǎn) 繆美華 孔祥偉

    摘要:海濱木槿(Hibiscus hamabo Sieb.et Zucc.)是一種典型的耐鹽植物,甚至能夠在海水中生長(zhǎng),是沿海灘涂綠化的絕佳綠化樹(shù)種,但在栽培試驗(yàn)中發(fā)現(xiàn)海濱木槿不耐寒,很難在北方地區(qū)越冬。為探究海濱木槿在低溫脅迫下的分子應(yīng)答機(jī)制,挖掘海濱木槿的抗寒基因,培育耐寒海濱木槿,通過(guò)轉(zhuǎn)錄組測(cè)序(RNA-seq)的試驗(yàn)方法分析海濱木槿在低溫脅迫下的差異基因表達(dá)情況。試驗(yàn)結(jié)果表明,低溫脅迫下的處理組與常溫對(duì)照相比,共有23 689個(gè)差異表達(dá)基因,其中上調(diào)表達(dá)基因10 256個(gè),下調(diào)表達(dá)基因13 433個(gè)。通過(guò)GO富集分析,發(fā)現(xiàn)差異基因大多富集在代謝過(guò)程、細(xì)胞過(guò)程、膜、催化活性等方面;通過(guò) KEGG的富集分析,發(fā)現(xiàn)這些差異基因大多富集在代謝途徑、次生代謝物的生物合成、光合作用、植物晝夜節(jié)律、植物激素信號(hào)轉(zhuǎn)導(dǎo)等方面,本研究還篩選了一些潛在的耐寒差異基因,如FKF1、 GI、PYL、PP2C、CML、CPK、SNRK2、NCED和BIN2等。本試驗(yàn)結(jié)果可為研究海濱木槿的耐寒機(jī)制以及培育海濱木槿耐寒品種提供技術(shù)參考和基礎(chǔ)。

    關(guān)鍵詞:海濱木槿;低溫脅迫;轉(zhuǎn)錄組;差異基因;耐寒

    中圖分類(lèi)號(hào):S718文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1002-1302(2023)05-0056-09

    海濱木槿(Hibiscus hamabo Sieb.et Zucc.) 是錦葵科木槿屬的小喬木或灌木,在我國(guó)主要分布于浙江舟山、寧波等沿海一帶,在國(guó)外,日本和朝鮮也有發(fā)現(xiàn)報(bào)道[1-2]。海濱木槿極耐鹽堿,能夠在含鹽量1.5%的泥灘生長(zhǎng),即使被海水淹浸,仍能正常生長(zhǎng)和開(kāi)花結(jié)實(shí);同時(shí)對(duì)土壤的要求不高,在pH值為5.5~8.0之間,無(wú)論是泥灘、沙灘、丘陵地都能正常生長(zhǎng)[3]。在以往的研究中發(fā)現(xiàn),海濱木槿是一種典型的鹽生植物,一種優(yōu)良的園林綠化樹(shù)種,能夠應(yīng)用到鹽堿地改良和海岸灘涂綠化中,具有很高的生態(tài)價(jià)值和經(jīng)濟(jì)價(jià)值[4-8]。江蘇地區(qū)灘涂面積較大,海濱木槿的應(yīng)用前景廣闊,但在引種栽培試驗(yàn)中,我們發(fā)現(xiàn)海濱木槿在江蘇北部地區(qū),如連云港地區(qū),地上部分枝條基本無(wú)法安全越冬,只能在第2年從根部另發(fā)新枝。目前,國(guó)內(nèi)外對(duì)海濱木槿的研究主要集中于馴化、形態(tài)及耐鹽機(jī)制等方面[4-9],對(duì)耐寒機(jī)制方面的研究較少。因此,基于轉(zhuǎn)錄組技術(shù)的海濱木槿耐寒機(jī)制的研究將為探索海濱木槿分子標(biāo)記等方面提供新的思路,促進(jìn)海濱木槿耐寒遺傳育種的發(fā)展。

    1 材料與方法

    1.1 試驗(yàn)材料與試驗(yàn)地點(diǎn)

    供試材料為1年生海濱木槿扦插盆栽苗,其母本為從浙江舟山引進(jìn),并在連云港經(jīng)過(guò)5年低溫馴化篩選出的耐寒品種。

    試驗(yàn)地點(diǎn)為連云港市農(nóng)業(yè)科學(xué)院內(nèi)。

    1.2 試驗(yàn)時(shí)間與方法

    于2019年6月20日將長(zhǎng)勢(shì)一致的1年生海濱木槿扦插苗置于人工氣候培養(yǎng)箱內(nèi),模擬低溫脅迫。以20 ℃處理為對(duì)照(CK),脅迫處理(T)溫度設(shè)置為4 ℃,低溫處理48 h采集組織,每個(gè)處理3次重復(fù),每組處理各取10張大小均勻的葉片,用錫箔紙包裹后迅速放入液氮中冷凍,然后放入-80 ℃超低溫冰箱中保存,用于總RNA的提取。

    1.3 總RNA提取

    使用Omega Bio-Tek RNA提取試劑盒(貨號(hào)R6827)提取RNA。利用瓊脂糖凝膠電泳和NanoDrop,檢驗(yàn)提取的RNA的完整性和純度;最后用Agilent 2100生物分析儀確定RNA的完整性。

    1.4 建庫(kù)和轉(zhuǎn)錄組測(cè)序

    建庫(kù)、Unigene基本注釋和序列分析方法參見(jiàn)趙春旭等的試驗(yàn)方法[10]。

    1.5 差異表達(dá)基因的GO分析和KEGG分析

    差異表達(dá)基因(differentially expressed genes,DEGs)分析采用DESeq2[11]進(jìn)行分析,篩選閾值FDR<0.05 且|log2FC|>1。Pathway顯著性富集分析以KEGG Pathway為單位,應(yīng)用超幾何檢驗(yàn),找出與整個(gè)基因組背景相比,在基因中顯著性富集的Pathway。

    2 結(jié)果與分析

    2.1 測(cè)序數(shù)據(jù)及質(zhì)量情況

    通過(guò)對(duì)原始數(shù)據(jù)進(jìn)行數(shù)據(jù)過(guò)濾,以減小無(wú)效數(shù)據(jù)所來(lái)帶的分析干擾。對(duì)下機(jī)的 raw reads利用fastp[12]進(jìn)行質(zhì)控,去除不合格的reads得到clean reads。從表1可以看出,未過(guò)濾前的原始序列數(shù)據(jù)(raw datas)各處理均在3 900萬(wàn)以上,過(guò)濾后的數(shù)據(jù)(clean data)占比均超過(guò)99.8%,不合格的序列數(shù)據(jù)較少。從表2可以看出,Phred數(shù)值大于20的堿基占總體堿基的百分比(Q20)在97%以上,Phred數(shù)值大于30的堿基占總體堿基的百分比(Q30)在93%以上。堿基G和C的數(shù)量總和占總堿基數(shù)量的百分比(GC含量)在46%~47% 之間。說(shuō)明測(cè)序的建庫(kù)質(zhì)量良好。

    2.2 轉(zhuǎn)錄組的組裝

    采用組裝序列的長(zhǎng)度分布來(lái)衡量轉(zhuǎn)錄組組裝質(zhì)量。由表3可知,共組裝了123 276個(gè)Unigene,平均長(zhǎng)度為807 bp,N50長(zhǎng)度為1 225 bp,最大長(zhǎng)度為 17 146 bp,最小長(zhǎng)度為201 bp。由圖1可知,組裝出來(lái)的 Unigene集中在0~1 000 bp之間,數(shù)據(jù)量及質(zhì)量均能保證后續(xù)分析要求。

    2.3 基因功能注釋

    本試驗(yàn)分別在NR、SwissProt、KEGG 和COG/KOG這4個(gè)數(shù)據(jù)庫(kù)繼續(xù)注釋?zhuān)Y(jié)果如圖2所示,4個(gè)數(shù)據(jù)庫(kù)可注釋的Unigene數(shù)量共76 175個(gè),其中有34 723個(gè)Unigene能在所有數(shù)據(jù)庫(kù)中得到注釋。

    通過(guò)NR數(shù)據(jù)庫(kù)比對(duì),結(jié)果(圖3)表明,與海濱木槿同源基因相似度最高的是陸地棉(Gossypium hirsutum),相似率達(dá)到24.36%;其次是亞洲棉(G. arboreum),相似率達(dá)到21.94%;然后是雷蒙德氏棉(G. raimondil),相似率達(dá)到20.68%;與可可樹(shù)(Theobroma cacao)的相似率達(dá)到18.60%。

    2.4 差異基因表達(dá)分析

    差異基因表達(dá)分析采用DESep2[11],篩選條件為FDR<0.05 且 |log2FC|>1的基因?yàn)轱@著差異基因。如圖4所示,低溫脅迫下的處理組與常溫對(duì)照相比,共有23 689個(gè)差異表達(dá)基因。其中上調(diào)表達(dá)基因有10 256個(gè),下調(diào)表達(dá)基因有 13 433 個(gè)。

    2.5 差異表達(dá)基因的GO分析

    將差異表達(dá)基因進(jìn)行GO富集功能分類(lèi),結(jié)果見(jiàn)圖5,主要從生物過(guò)程(biological process,BP)、細(xì)胞組分(cellular component,CC)、分子功能(molecular function,MF)3個(gè)方面進(jìn)行分析。

    在生物學(xué)過(guò)程中,試驗(yàn)材料DEGs 富集最多的依次為代謝過(guò)程(metabolic process),有6 595 個(gè)Unigene得到注釋?zhuān)急?1.3%;細(xì)胞過(guò)程(cellular process),有6 019 個(gè)Unigene得到注釋?zhuān)急?74.2%;單有機(jī)體過(guò)程(single-organism process),有4 653 個(gè)Unigene得到注釋?zhuān)急?57.36%;生物調(diào)節(jié)(biological regulation),有2 493 個(gè)Unigene得到注釋?zhuān)?占比30.73%; 生物過(guò)程的調(diào)節(jié)(regulation of biological process),有2 284個(gè)Unigene得到注釋?zhuān)急?28.16%;對(duì)刺激的反應(yīng)(response to stimulus),有1 954 個(gè)Unigene得到注釋?zhuān)急?4.09%;本土化(localization),有1 556個(gè)Unigene得到注釋?zhuān)急?9.18%;發(fā)展過(guò)程(developmental process),有1 089個(gè)Unigene得到注釋?zhuān)急?3.42%;細(xì)胞成分組織(cellular component organization or biogenesis),有 1 197 個(gè)Unigene得到注釋?zhuān)急?4.76%;信號(hào)(signaling),有712 個(gè)Unigene得到注釋?zhuān)急?.78%; 多細(xì)胞生物的過(guò)程(multicellular organismal process),有834 個(gè)Unigene得到注釋?zhuān)急?0.28%;多有機(jī)體過(guò)程(multi-organism process),有390個(gè)Unigene得到注釋?zhuān)急?.81%;生殖(reproduction),有516個(gè)Unigene得到注釋?zhuān)急?.36%;生殖過(guò)程(reproductive process),有511個(gè)Unigene得到注釋?zhuān)急?.3%;免疫系統(tǒng)過(guò)程(immune system process),有175個(gè)Unigene得到注釋?zhuān)急?.16%;生物過(guò)程的正向調(diào)控(positive regulation of biological process),有203個(gè)Unigene得到注釋?zhuān)急?.5%;生長(zhǎng)(growth),有159個(gè)Unigene得到注釋?zhuān)急?.96%;生物過(guò)程的負(fù)調(diào)控(negative regulation of biological process),有189個(gè)Unigene得到注釋?zhuān)急?.33%;節(jié)奏過(guò)程(rhythmic process),有62 個(gè)Unigene得到注釋?zhuān)急?.76%;生物黏附(biological adhesion),有40個(gè)Unigene得到注釋?zhuān)急?.49%;排毒(detoxification),有10 個(gè)Unigene得到注釋?zhuān)急?.12%;運(yùn)動(dòng)(locomotion),有5個(gè)Unigene得到注釋?zhuān)急?.06%;細(xì)胞死亡(cell killing),有1個(gè)Unigene得到注釋?zhuān)急?.01%。

    在細(xì)胞組分中,試驗(yàn)材料DEGs 富集最多的依次為細(xì)胞(cell和cell part),有3 832 個(gè)Unigene得到注釋?zhuān)急?4.99%;細(xì)胞器(organelle),有3 173個(gè)Unigene得到注釋?zhuān)急?62.09%;膜(membrane),有2 142個(gè)Unigene得到注釋?zhuān)急?41.92%;膜部分(membrane part),有1 536個(gè)Unigene得到注釋?zhuān)急?30.06%;細(xì)胞器部分(organelle part),有1 392 個(gè)Unigene得到注釋?zhuān)急?7.24%;大分子復(fù)合物(macromolecular complex),有742個(gè)Unigene得到注釋?zhuān)急?4.52%;細(xì)胞連接(cell junction),有228 個(gè)Unigene得到注釋?zhuān)急?.46%;細(xì)胞外區(qū)(extracellular region),有134 個(gè)Unigene得到注釋?zhuān)急?.62%;膜封閉腔(membrane-enclosed lumen),有74 個(gè)Unigene得到注釋?zhuān)急?.45%;病毒粒子(virion),有20個(gè)Unigene得到注釋?zhuān)急?.39%;類(lèi)核素(nucleoid),有8 個(gè)Unigene得到注釋?zhuān)急?.16%;細(xì)胞外區(qū)部分(extracellular region part),有4個(gè)Unigene得到注釋?zhuān)急?.08%;細(xì)胞外基質(zhì)(extracellular matrix),有4 個(gè)Unigene得到注釋?zhuān)急?.08%;超分子纖維(supramolecular fiber),有3個(gè)Unigene得到注釋?zhuān)急?.06%。

    在分子功能中,試驗(yàn)材料DEGs 富集最多的依次為催化活性(catalytic activity),有5 259 個(gè)Unigene得到注釋?zhuān)急?8.86%;黏合物(binding),有4 632個(gè)Unigene得到注釋?zhuān)急?0.65%;轉(zhuǎn)運(yùn)蛋白活性(transporter activity),有454個(gè)Unigene得到注釋?zhuān)急?.94%;核酸結(jié)合轉(zhuǎn)錄因子活性(nucleic acid binding transcription factor activity),有382個(gè)Unigene得到注釋?zhuān)急?%;信號(hào)傳感器活動(dòng)(signal transducer activity),有119個(gè)Unigene得到注釋?zhuān)急?.56%;結(jié)構(gòu)分子活性(structural molecule activity),有85個(gè)Unigene得到注釋?zhuān)急?.11%;分子換能器活性(molecular transducer activity),有81 個(gè)Unigene得到注釋?zhuān)急?.06%;分子功能調(diào)節(jié)劑(molecular function regulator),有42 個(gè)Unigene得到注釋?zhuān)急?.55%;抗氧化活性(antioxidant activity),有18個(gè)Unigene得到注釋?zhuān)急?.24%;轉(zhuǎn)錄因子活性,蛋白質(zhì)結(jié)合(transcription factor activity,protein binding),有14個(gè)Unigene得到注釋?zhuān)急?.18%;電子載流子活性(electron carrier activity),有10個(gè)Unigene得到注釋?zhuān)急?.13%;翻譯調(diào)節(jié)活性(translation regulator activity),有1個(gè)Unigene得到注釋?zhuān)急?.01%。

    其中,按富集程度(qvalue<0.05)由高到底排序依次為質(zhì)體(plastid)、質(zhì)體部分(plastid part)、類(lèi)囊體(thylakoid)、質(zhì)體類(lèi)囊體(plastid thylakoid)、葉綠體(chloroplast)、葉綠體部分(chloroplast part)、質(zhì)體基質(zhì)(plastid stroma)、質(zhì)體被膜(plastid envelope)等(圖6)。

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

    生物個(gè)體的功能表達(dá)通常需要大量基因的共同協(xié)作才能完成,因此通過(guò)對(duì)代謝通路的分析,能夠更加清晰地表達(dá) DEGs 的生物學(xué)功能。將差異表達(dá)基因進(jìn)行KEGG分析,并在130個(gè)代謝通路中得到注釋。富集程度(qvalue<0.05)最高的20個(gè)通路分別是:代謝途徑(metabolic pathways)、次生代謝物的生物合成(biosynthesis of secondary metabolites)、光合作用(photosynthesis)、植物的晝夜節(jié)律(circadian rhythm-plant)、卟啉和葉綠素代謝(porphyrin and chlorophyll metabolism)、光合作用-天線蛋白(photosynthesis-antenna proteins)、甘氨酸、絲氨酸和蘇氨酸代謝(glycine,serine and threonine metabolism)、光合生物中的碳固定(carbon fixation in photosynthetic organisms)、二萜生物合成(diterpenoid biosynthesis)、類(lèi)胡蘿卜素生物合成(carotenoid biosynthesis)、碳代謝(carbon metabolism)、植物-病原體相互作用(plant-pathogen interaction)、半胱氨酸和蛋氨酸代謝(cysteine and methionine metabolism)、氮代謝(nitrogen metabolism)、植物激素信號(hào)轉(zhuǎn)導(dǎo)(plant hormone signal transduction)等。表示富集程度的qvalue均<0.001,其中得到注釋的 DEGs數(shù)量最多的代謝通路是代謝途徑,共有1 752個(gè)。其次,次生代謝物的生物合成有995個(gè)DEGs,植物的晝夜節(jié)律有121個(gè)DEGs,甘氨酸、絲氨酸和蘇氨酸代謝有100個(gè)DEGs(圖7)。

    2.7 抗寒差異表達(dá)基因

    本試驗(yàn)在對(duì)差異表達(dá)基因功能分析的基礎(chǔ)上,對(duì)試驗(yàn)材料特異表達(dá)的基因進(jìn)行進(jìn)一步的分析,挑選出一些與植物耐寒相關(guān)的基因,發(fā)現(xiàn)這些基因主要分布在代謝途徑、次生代謝物的生物合成、卟啉和葉綠素代謝、氨基酸代謝、類(lèi)胡蘿卜素生物合成、光合作用與碳代謝、晝夜節(jié)律、植物激素信號(hào)轉(zhuǎn)導(dǎo)等通路中。

    如表4所示,與代謝途徑相關(guān)的表達(dá)基因有E1.6.5.4、NCED、E2.4.1.13;與次生代謝物的生物合成相關(guān)的表達(dá)基因有katE,CAT,catB,srpA、G6PD,zwf;與卟啉和葉綠素代謝相關(guān)的表達(dá)基因有hemD,UROS、hemA;與氨基酸代謝相關(guān)的表達(dá)基因有AGXT2、serC,PSAT1;與類(lèi)胡蘿卜素生物合成相關(guān)的表達(dá)基因有CYP707A、AAO3、ZEP,ABA1;與光合用與碳代謝相關(guān)的表達(dá)基因有petH、psaG、LHCB6、MDH2、IDH3、PGD,gnd,gntZ;與植物晝夜節(jié)律相關(guān)的表達(dá)基因有GI、FKF1;與植物激素信號(hào)轉(zhuǎn)導(dǎo)相關(guān)的表達(dá)基因有CML、CPK、SNRK2、BIN2、IAA、PP2C、PYL。

    3 討論與結(jié)論

    溫度是植物生長(zhǎng)發(fā)育過(guò)程中最受環(huán)境制約的因素之一,限制了植物的生長(zhǎng)發(fā)育和地理分布[13],以至于有了“江南冬景似春華,江北冬嶺只枯丫”的對(duì)比。植物在受到低溫脅迫的時(shí)候,細(xì)胞膜系統(tǒng)首先遭到破壞,膜透性改變,繼而擾亂細(xì)胞內(nèi)正常代謝,如光合作用、呼吸作用等,嚴(yán)重時(shí)導(dǎo)致細(xì)胞凋亡[14-19]。為了研究植物體內(nèi)應(yīng)答低溫脅迫的調(diào)控機(jī)制,隨著近年來(lái)技術(shù)的發(fā)展,RNA-seq、Chip-seq、iTRAQ、TMT、LC-MS/MS 等新一代技術(shù)應(yīng)運(yùn)而生,分別從轉(zhuǎn)錄組、蛋白組、代謝組及轉(zhuǎn)錄翻譯調(diào)控等不同方向?qū)χ参锏目购畽C(jī)制進(jìn)行解析[20]。其中,轉(zhuǎn)錄組測(cè)序(RNA-seq)是對(duì)某一物種的mRNA進(jìn)行高通量測(cè)序的技術(shù),反映在特定條件和特定時(shí)間下該物種的基因表達(dá)情況[21-22],因此,通過(guò)轉(zhuǎn)錄組測(cè)序,有利于揭示作物在冷脅迫下的生理及分子機(jī)制。目前為止,轉(zhuǎn)錄組測(cè)序手段已經(jīng)揭示了小麥、水稻、茶樹(shù)、薰衣草、草地早熟禾、煙草、擬南芥、西葫蘆等經(jīng)濟(jì)作物的耐寒分子機(jī)制,為培育其耐寒品種提供了信息技術(shù)參考[10,13,15,23-24]。

    本研究利用轉(zhuǎn)錄組測(cè)序方法,研究了經(jīng)過(guò)低溫馴化的耐寒海濱木槿在低溫脅迫處理下的基因表達(dá)情況。經(jīng)過(guò)轉(zhuǎn)錄組測(cè)序、數(shù)據(jù)組裝和注釋后,通過(guò)NR數(shù)據(jù)庫(kù)比對(duì)與海濱木槿同源基因相似度最高的是陸地棉(G. hirsutum),相似率達(dá)到24.36%;其次是亞洲棉(G. arboreum),相似率達(dá)到21.94%;然后是雷蒙德氏棉(G. raimondil)的相似率達(dá)到20.68%;可可樹(shù)(Theobroma cacao),相似率達(dá)到18.60%。其中,陸地棉、亞洲棉、雷蒙德氏棉均為錦葵科木槿族棉屬,與海濱木槿同為錦葵科植物,可可樹(shù)也是錦葵目植物,說(shuō)明測(cè)序準(zhǔn)確性是比較高的。采用DESep2進(jìn)行差異基因表達(dá)分析,篩選出FDR<0.05 且 |log2FC|>1 的基因?yàn)轱@著差異基因,共得到23 689個(gè)差異表達(dá)基因,其中上調(diào)表達(dá)基因有10 256個(gè),下調(diào)表達(dá)基因有13 433個(gè)。通過(guò)GO富集分析,發(fā)現(xiàn)DEGs主要富集在代謝過(guò)程、細(xì)胞過(guò)程、膜、催化活性、黏合物等方面。這與相關(guān)報(bào)道較為一致,植物在遭遇低溫脅迫后,首先細(xì)胞膜透性會(huì)改變,傳遞信號(hào),調(diào)控細(xì)胞代謝,以應(yīng)對(duì)低溫帶來(lái)的傷害[25-26]。通過(guò) KEGG的富集分析,發(fā)現(xiàn)差異基因主要富集在代謝途徑、次生代謝物的生物合成、光合作用、植物晝夜節(jié)律、卟啉和葉綠素代謝、甘氨酸、絲氨酸和蘇氨酸代謝、二萜生物合成、類(lèi)胡蘿卜素生物合成、碳代謝、植物-病原體相互作用、半胱氨酸和蛋氨酸代謝、氮代謝、植物激素信號(hào)轉(zhuǎn)導(dǎo)等,說(shuō)明這些途徑可能在海濱木槿低溫脅迫響應(yīng)中發(fā)揮作用。

    通過(guò)進(jìn)一步的發(fā)掘,發(fā)現(xiàn)耐寒相關(guān)的差異基因主要分布在代謝途徑、次生代謝物的生物合成、卟啉和葉綠素代謝、氨基酸代謝、類(lèi)胡蘿卜素生物合成、光合作用與碳代謝、晝夜節(jié)律、植物激素信號(hào)轉(zhuǎn)導(dǎo)等通路中。梁雨晨等認(rèn)為,在薰衣草晝夜節(jié)律通路中,植物受光敏色素(PHY)調(diào)控,誘導(dǎo)產(chǎn)生相關(guān)抗性基因,通過(guò)蛋白網(wǎng)絡(luò)互作分析和篩選,發(fā)現(xiàn)樞紐基因 FKF1和GI,在晝夜節(jié)律通路中上調(diào)表達(dá),參與調(diào)控植物晝夜節(jié)律信號(hào)通路,表明生物鐘調(diào)控在薰衣草應(yīng)答寒冷脅迫中發(fā)揮了重要作用[24]。本研究中,篩選出晝夜節(jié)律中的差異基因GI和FKF1均為上調(diào)表達(dá),與梁雨晨的研究結(jié)果[24]一致。Chen等通過(guò)對(duì)經(jīng)過(guò)低溫脅迫處理后耐寒基因型菠蘿的轉(zhuǎn)錄組測(cè)序分析,DEGs 通過(guò) KEGG 分析,發(fā)現(xiàn)與細(xì)胞壁特性、氣孔關(guān)閉以及脫落酸(ABA)和活性氧(ROS)信號(hào)轉(zhuǎn)導(dǎo)相關(guān)的基因在菠蘿耐寒性中起重要作用,其中,所有與ABA信號(hào)轉(zhuǎn)導(dǎo)相關(guān)的冷脅迫響應(yīng)基因表達(dá)均上調(diào),包括PYR/PYL同源基因和PP2C同源基因[27]。本研究與之相一致,在本研究中,篩選出的差異基因PYL和PP2C均上調(diào)表達(dá)。趙春旭等在探究青海野生草地早熟禾低溫脅迫下的分子應(yīng)答機(jī)制時(shí),發(fā)現(xiàn)在低溫脅迫下,差異基因主要富集在光合作用、氧化還原反應(yīng)過(guò)程、碳水化合物代謝、細(xì)胞膜系統(tǒng)、轉(zhuǎn)運(yùn)蛋白以及次生物代謝中,另外通過(guò)耐寒和不耐寒2種材料轉(zhuǎn)錄組的對(duì)比,篩選出一批在耐寒材料中上調(diào)表達(dá)的潛在的抗寒基因,如CML、CPK、CALM、DHAR、GST、NCED、SNRK2、BSK、CKX、BIN2、ARF和PEK等[11]。本研究的結(jié)果與之一致,在本研究中,篩選出的CML、CPK、SNRK2、NCED、BIN2等與趙春旭研究中重合的基因,也均上調(diào)表達(dá)。

    本研究利用高通量轉(zhuǎn)錄組測(cè)序技術(shù)建立了海濱木槿低溫脅迫前后葉片的轉(zhuǎn)錄組數(shù)據(jù)庫(kù),篩選出差異表達(dá)基因,為研究海濱木槿的耐寒機(jī)制以及培育海濱木槿耐寒品種提供了技術(shù)參考和基礎(chǔ)。

    參考文獻(xiàn):

    [1]張朝芳. 浙江植物志:第一卷 蕨類(lèi)植物·裸子植物[M]. 杭州:浙江科學(xué)技術(shù)出版社,1993.

    [2]張若蕙. 浙江珍稀瀕危植物[M]. 杭州:浙江科學(xué)技術(shù)出版社,1994:330-339.

    [3]唐麗紅,馬明睿,韓 華,等. 上海市景觀水體水生植物現(xiàn)狀及配置評(píng)價(jià)[J]. 生態(tài)學(xué)雜志,2013,32(3):563-570.

    [4]肖祖飛,趙 姣,張 杰,等. 海濱木槿繁殖技術(shù)[J]. 北方園藝,2017(12):74-77.

    [5]於朝廣,殷云龍,蘆治國(guó),等. 海濱木槿研究進(jìn)展及應(yīng)用前景[J]. 現(xiàn)代農(nóng)業(yè)科技,2016(19):140-143,145.

    [6]李翠云,姜彥成,喬桂榮,等. 海濱木槿葉片蛋白質(zhì)雙向電泳體系的建立[J]. 西南林業(yè)大學(xué)學(xué)報(bào),2012,32(4):30-35.

    [7]王宗星. 海濱木槿和蠟楊梅對(duì)海水脅迫的光合生理響應(yīng)[D]. 北京:中國(guó)林業(yè)科學(xué)研究院,2011.

    [8]王連吉. 新優(yōu)樹(shù)種海濱木槿在風(fēng)景園林中的應(yīng)用[J]. 中國(guó)園林,2010,26(4):49-50.

    [9]王秀麗,張 荻,劉紅梅,等. 海濱木槿耐鹽性的初步研究[J]. 上海交通大學(xué)學(xué)報(bào)(農(nóng)業(yè)科學(xué)版),2010,28(3):248-254.

    [10]趙春旭,馬 祥,董文科,等. 低溫脅迫下不同青海野生草地早熟禾的轉(zhuǎn)錄組比較分析[J]. 草地學(xué)報(bào),2020,28(2):305-318.

    [11]Love M I,Huber W,Anders S.Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2[J]. Genome Biology,2014,15(12):550.

    [12]Chen S F,Zhou Y Q,Chen Y R,et al. Fastp:an ultra-fast all-in-one FASTQ preprocessor[J]. Bioinformatics,2018,34(17):1884-1890.

    [13]張 健,唐 露,冉啟凡,等. 植物響應(yīng)低溫脅迫轉(zhuǎn)錄組測(cè)序研究進(jìn)展[J]. 分子植物育種,2020,18(6):1849-1866.

    [14]Wang H S,Yu C,Zhu Z J,et al. Overexpression in tobacco of a tomato GMPase gene improves tolerance to both low and high temperature stress by enhancing antioxidation capacity[J]. Plant Cell Reports,2011,30(6):1029-1040.

    [15]李 楠,鄭勇奇,丁紅梅,等. 低溫脅迫下短枝木麻黃耐寒相關(guān)基因的差異表達(dá)分析[J]. 林業(yè)科學(xué),2017,53(7):62-71.

    [16]劉文盈. 低溫脅迫對(duì)黑果枸杞愈傷組織的影響[J]. 基因組學(xué)與應(yīng)用生物學(xué),2018,37(1):408-412.

    [17]Colinet H,Renault D,Roussel D. Cold acclimation allows Drosophila flies to maintain mitochondrial functioning under cold stress[J]. Insect Biochemistry and Molecular Biology,2017,80:52-60.

    [18]Walker D J,Romero P,Correal E. Cold tolerance,water relations and accumulation of osmolytes in Bituminaria bituminosa[J]. Biologia Plantarum,2010,54(2):293-298.

    [19]宋麗莉. 紫花苜??购D(zhuǎn)錄組測(cè)序分析及MsNCX和MsCML基因的功能研究[D]. 哈爾濱:哈爾濱師范大學(xué),2017.

    [20]張 純,唐承晨,王吉永,等. 轉(zhuǎn)錄組學(xué)在植物應(yīng)答逆境脅迫中的研究進(jìn)展[J]. 生物學(xué)雜志,2017,34(2):86-90.

    [21]慧 芳,劉秀巖,李宗諭,等. 轉(zhuǎn)錄組測(cè)序技術(shù)在藥用植物研究中的應(yīng)用[J]. 中草藥,2019,50(24):6149-6155.

    [22]王麗鴛. 基于EST數(shù)據(jù)庫(kù)和轉(zhuǎn)錄組測(cè)序的茶樹(shù)DNA分子標(biāo)記開(kāi)發(fā)與應(yīng)用研究[D]. 北京:中國(guó)農(nóng)業(yè)科學(xué)院,2011.

    [23]張宏亮. 溫脅迫下西葫蘆轉(zhuǎn)錄組分析與SSR分子標(biāo)記開(kāi)發(fā)[D]. 太谷:山西農(nóng)業(yè)大學(xué),2015.

    [24]梁雨晨,藺海嬌,侯一航,等. 狹葉薰衣草葉片轉(zhuǎn)錄組耐寒相關(guān)差異表達(dá)基因分析[J/OL]. 分子植物育種:1-14[2022-11-10]. http://kns.cnki.net/kcms/detail/46.1068.S.20210604.1045.010.html.

    [25]Orvar B L,Sangwan V,Omann F,et al. Early steps in cold sensing by plant cells:the role of actin cytoskeleton and membrane fluidity[J]. The Plant Journal,2000,23(6):785-794.

    [26]范吉標(biāo). 狗牙根耐寒生理及分子機(jī)制解析[D]. 武漢:中國(guó)科學(xué)院研究生院(武漢植物園),2016.

    [27]Chen C J,Zhang Y F,Xu Z Q,et al. Transcriptome profiling of the pineapple under low temperature to facilitate its breeding for cold tolerance[J]. PLoS One,2016,11(9):e0163315.

    收稿日期:2022-05-16

    基金項(xiàng)目:江蘇省林業(yè)科技創(chuàng)新與推廣項(xiàng)目(編號(hào):LYKJ-連云港[2021]01);連云港市財(cái)政專(zhuān)項(xiàng)(編號(hào):QNJJ1921)。

    作者簡(jiǎn)介:趙文靜(1991—),女,江蘇連云港人,碩士,助理研究員,主要從事植物病理學(xué)和植物育種研究。E-mail:ZWJ6691@163.com。

    通信作者:劉興滿(mǎn),碩士,副研究員,主要從事園藝和植物育種。E-mail:13611551930@126.com。

    猜你喜歡
    轉(zhuǎn)錄組低溫脅迫差異基因
    ICR鼠肝和腎毒性損傷生物標(biāo)志物的篩選
    基于RNA 測(cè)序研究人參二醇對(duì)大鼠心血管內(nèi)皮細(xì)胞基因表達(dá)的影響 (正文見(jiàn)第26 頁(yè))
    基于轉(zhuǎn)錄組測(cè)序的山茱萸次生代謝生物合成相關(guān)基因的挖掘
    金釵石斛轉(zhuǎn)錄組SSR位點(diǎn)信息分析
    根據(jù)萌發(fā)率和出苗率篩選高粱種子萌發(fā)期耐低溫材料
    人參屬藥用植物轉(zhuǎn)錄組研究進(jìn)展
    遮陽(yáng)網(wǎng)覆蓋對(duì)枇杷幼果抗寒性影響機(jī)理初探
    轉(zhuǎn)ICE1基因水稻耐冷性與膜脂過(guò)氧化和抗氧化酶活性的關(guān)系
    江蘇農(nóng)業(yè)科學(xué)(2015年8期)2015-09-10 14:49:14
    大豆轉(zhuǎn)錄組測(cè)序研究進(jìn)展綜述
    日日摸夜夜添夜夜添av毛片| 一级黄片播放器| 大香蕉97超碰在线| 久久午夜福利片| av卡一久久| 国精品久久久久久国模美| 日韩大片免费观看网站| 国产一区二区三区av在线| 午夜福利在线在线| 你懂的网址亚洲精品在线观看| 99re6热这里在线精品视频| 欧美变态另类bdsm刘玥| 真实男女啪啪啪动态图| 内地一区二区视频在线| 内地一区二区视频在线| 最近最新中文字幕大全电影3| 免费无遮挡裸体视频| 自拍偷自拍亚洲精品老妇| 色哟哟·www| 成年人午夜在线观看视频 | 人人妻人人澡欧美一区二区| 国产综合精华液| 97精品久久久久久久久久精品| 国产免费又黄又爽又色| 我的老师免费观看完整版| av免费在线看不卡| 国产在视频线精品| 最近最新中文字幕大全电影3| 色网站视频免费| 亚洲国产av新网站| 国产亚洲av嫩草精品影院| 高清在线视频一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 一个人看视频在线观看www免费| 亚洲av国产av综合av卡| 国产黄色小视频在线观看| 麻豆精品久久久久久蜜桃| 夫妻性生交免费视频一级片| 少妇被粗大猛烈的视频| 一区二区三区免费毛片| 国产精品麻豆人妻色哟哟久久 | 水蜜桃什么品种好| 99久久精品一区二区三区| 欧美xxⅹ黑人| 久久久精品94久久精品| 最近2019中文字幕mv第一页| 中文资源天堂在线| 亚洲欧美精品专区久久| 欧美一区二区亚洲| 美女高潮的动态| 亚洲国产日韩欧美精品在线观看| 亚洲熟女精品中文字幕| 亚洲aⅴ乱码一区二区在线播放| 成人亚洲精品av一区二区| 国产探花极品一区二区| 身体一侧抽搐| 我的女老师完整版在线观看| 国语对白做爰xxxⅹ性视频网站| 高清在线视频一区二区三区| 成年女人看的毛片在线观看| 国国产精品蜜臀av免费| 国产69精品久久久久777片| 91精品一卡2卡3卡4卡| 亚洲精品自拍成人| 日韩成人av中文字幕在线观看| 国产精品伦人一区二区| 久久鲁丝午夜福利片| 国产 一区精品| 少妇裸体淫交视频免费看高清| 啦啦啦啦在线视频资源| 亚洲经典国产精华液单| 五月玫瑰六月丁香| 少妇人妻精品综合一区二区| 国产探花极品一区二区| 国产精品一及| 乱系列少妇在线播放| 精品久久久久久久人妻蜜臀av| 国产黄色免费在线视频| 日日摸夜夜添夜夜爱| 亚洲国产成人一精品久久久| .国产精品久久| 亚洲伊人久久精品综合| 好男人在线观看高清免费视频| 亚洲在线观看片| 中文天堂在线官网| 国产视频首页在线观看| 女的被弄到高潮叫床怎么办| 国产午夜福利久久久久久| 国产成人精品婷婷| 99热这里只有是精品在线观看| 久久久久久久久久成人| 婷婷六月久久综合丁香| 老司机影院成人| 欧美潮喷喷水| 蜜臀久久99精品久久宅男| 国产精品女同一区二区软件| 国产单亲对白刺激| 精品国产露脸久久av麻豆 | 亚洲美女搞黄在线观看| 精品国产露脸久久av麻豆 | 伊人久久国产一区二区| 男的添女的下面高潮视频| 激情五月婷婷亚洲| 一级av片app| 精品酒店卫生间| 亚洲欧洲国产日韩| 中文乱码字字幕精品一区二区三区 | 又大又黄又爽视频免费| 日韩亚洲欧美综合| 久久久精品免费免费高清| 人人妻人人澡欧美一区二区| av在线播放精品| 成年免费大片在线观看| 极品教师在线视频| 中文精品一卡2卡3卡4更新| 亚洲怡红院男人天堂| 免费电影在线观看免费观看| 午夜激情久久久久久久| 精品一区在线观看国产| 午夜福利视频精品| 国产精品一区www在线观看| 91精品国产九色| 欧美97在线视频| 国产毛片a区久久久久| 亚洲精品视频女| 亚洲国产av新网站| 偷拍熟女少妇极品色| 亚洲高清免费不卡视频| 99热6这里只有精品| 中文字幕人妻熟人妻熟丝袜美| 尾随美女入室| 国产精品一区二区性色av| 国产精品日韩av在线免费观看| 色视频www国产| 少妇高潮的动态图| 九色成人免费人妻av| 秋霞在线观看毛片| .国产精品久久| 亚洲av福利一区| 午夜激情欧美在线| 人人妻人人看人人澡| 日韩av在线免费看完整版不卡| 国产黄a三级三级三级人| 免费在线观看成人毛片| 欧美日韩综合久久久久久| 久久人人爽人人爽人人片va| 亚洲国产成人一精品久久久| 午夜福利高清视频| 搡老妇女老女人老熟妇| 青春草视频在线免费观看| 只有这里有精品99| 丰满少妇做爰视频| 欧美变态另类bdsm刘玥| 亚洲av.av天堂| 国产成人aa在线观看| 伦精品一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 婷婷色综合www| 久久久精品免费免费高清| 97人妻精品一区二区三区麻豆| 亚洲国产最新在线播放| 赤兔流量卡办理| 欧美日韩精品成人综合77777| 午夜老司机福利剧场| 欧美潮喷喷水| 少妇丰满av| 观看美女的网站| 国产精品伦人一区二区| 日韩国内少妇激情av| 亚洲国产精品sss在线观看| 最近中文字幕高清免费大全6| 免费观看无遮挡的男女| 干丝袜人妻中文字幕| 高清视频免费观看一区二区 | 久久精品综合一区二区三区| 午夜激情福利司机影院| 免费看不卡的av| 亚洲激情五月婷婷啪啪| 听说在线观看完整版免费高清| 99热这里只有是精品在线观看| 国产伦一二天堂av在线观看| 久久国产乱子免费精品| 日韩强制内射视频| 久久久久性生活片| 国产激情偷乱视频一区二区| 蜜臀久久99精品久久宅男| 国精品久久久久久国模美| 乱码一卡2卡4卡精品| 欧美区成人在线视频| 我的女老师完整版在线观看| 亚洲av免费高清在线观看| 中文字幕久久专区| 精品人妻偷拍中文字幕| 久久久久久国产a免费观看| 免费黄色在线免费观看| 免费观看精品视频网站| 国产有黄有色有爽视频| 精品一区二区免费观看| 久久午夜福利片| videos熟女内射| av福利片在线观看| 国产精品av视频在线免费观看| 少妇裸体淫交视频免费看高清| 国产成人aa在线观看| 国产高清三级在线| 看十八女毛片水多多多| 国产真实伦视频高清在线观看| 爱豆传媒免费全集在线观看| 嫩草影院入口| 久久久久久久大尺度免费视频| 国产伦理片在线播放av一区| 成人毛片a级毛片在线播放| 99久久中文字幕三级久久日本| 免费大片18禁| 女人被狂操c到高潮| 黄色日韩在线| 偷拍熟女少妇极品色| 国产高清三级在线| 婷婷六月久久综合丁香| 麻豆av噜噜一区二区三区| 亚洲国产高清在线一区二区三| 青春草亚洲视频在线观看| 观看免费一级毛片| 少妇丰满av| 最近2019中文字幕mv第一页| 亚洲欧美精品自产自拍| 精品久久久久久久末码| 国内精品一区二区在线观看| 深爱激情五月婷婷| 精品一区二区三卡| 我的女老师完整版在线观看| 蜜桃久久精品国产亚洲av| 男女啪啪激烈高潮av片| videossex国产| 色5月婷婷丁香| 黑人高潮一二区| 久热久热在线精品观看| 免费观看在线日韩| av国产久精品久网站免费入址| 亚洲欧美成人精品一区二区| 国产成人福利小说| 免费黄频网站在线观看国产| 日韩精品有码人妻一区| 免费看av在线观看网站| 久久久午夜欧美精品| 禁无遮挡网站| 国产精品.久久久| 国产亚洲最大av| 午夜久久久久精精品| av播播在线观看一区| 大香蕉97超碰在线| 人体艺术视频欧美日本| 在线观看免费高清a一片| 在线播放无遮挡| av又黄又爽大尺度在线免费看| 久久草成人影院| 成人毛片60女人毛片免费| 亚洲国产精品成人综合色| 小蜜桃在线观看免费完整版高清| 亚洲18禁久久av| 午夜福利成人在线免费观看| 国产 一区精品| 水蜜桃什么品种好| 欧美xxⅹ黑人| 免费观看的影片在线观看| 天天一区二区日本电影三级| 欧美另类一区| 成人欧美大片| av又黄又爽大尺度在线免费看| 国产91av在线免费观看| 街头女战士在线观看网站| 五月天丁香电影| 国产男女超爽视频在线观看| 久久6这里有精品| 少妇的逼水好多| 国产一区亚洲一区在线观看| 中文字幕人妻熟人妻熟丝袜美| 精华霜和精华液先用哪个| 欧美日韩国产mv在线观看视频 | 青青草视频在线视频观看| 99热这里只有是精品50| www.色视频.com| 欧美丝袜亚洲另类| 91午夜精品亚洲一区二区三区| 高清欧美精品videossex| 免费大片黄手机在线观看| a级毛片免费高清观看在线播放| 免费黄频网站在线观看国产| 亚洲av.av天堂| 亚洲av电影不卡..在线观看| 婷婷色综合大香蕉| 久热久热在线精品观看| 91aial.com中文字幕在线观看| 日韩 亚洲 欧美在线| 亚洲国产高清在线一区二区三| 久久97久久精品| 在线播放无遮挡| 大码成人一级视频| 超碰97精品在线观看| 色网站视频免费| 国产精品三级大全| 下体分泌物呈黄色| 午夜免费观看性视频| 国产极品天堂在线| 国产一区二区 视频在线| 婷婷色av中文字幕| 在线天堂最新版资源| 丰满乱子伦码专区| 久久久久久人妻| 日韩熟女老妇一区二区性免费视频| 中文字幕人妻丝袜一区二区 | 国产白丝娇喘喷水9色精品| 黑丝袜美女国产一区| 久久ye,这里只有精品| 18禁动态无遮挡网站| 久久久国产精品麻豆| 在线观看三级黄色| 免费看av在线观看网站| 免费高清在线观看日韩| 欧美少妇被猛烈插入视频| 少妇的丰满在线观看| 欧美日韩视频高清一区二区三区二| 在线观看www视频免费| 老司机影院毛片| 爱豆传媒免费全集在线观看| 三级国产精品片| 亚洲,欧美精品.| 久久这里只有精品19| 免费日韩欧美在线观看| 国产毛片在线视频| 日韩制服骚丝袜av| 国产 精品1| 午夜免费男女啪啪视频观看| 亚洲一区二区三区欧美精品| 亚洲色图综合在线观看| 亚洲伊人色综图| 久久精品国产亚洲av高清一级| 国产精品一区二区在线不卡| 自线自在国产av| 麻豆av在线久日| 色视频在线一区二区三区| 午夜福利乱码中文字幕| 国产成人91sexporn| 亚洲av在线观看美女高潮| 人人妻人人澡人人看| 香蕉精品网在线| 欧美日韩视频精品一区| 中文字幕制服av| 日本午夜av视频| 久久久久久免费高清国产稀缺| 国产极品天堂在线| 日韩中文字幕欧美一区二区 | 国产av码专区亚洲av| 大片免费播放器 马上看| 国产在线免费精品| 久久午夜福利片| 免费高清在线观看日韩| 最近手机中文字幕大全| 亚洲av日韩在线播放| 亚洲精品av麻豆狂野| 午夜免费观看性视频| 久久久精品国产亚洲av高清涩受| 丝瓜视频免费看黄片| 日韩不卡一区二区三区视频在线| 国产乱人偷精品视频| 亚洲伊人久久精品综合| av在线app专区| 日日撸夜夜添| 天天影视国产精品| 日本猛色少妇xxxxx猛交久久| 日本欧美视频一区| 男女边吃奶边做爰视频| 18禁国产床啪视频网站| 久久久久久久精品精品| 久久久国产精品麻豆| 亚洲成人手机| 国产探花极品一区二区| 一级毛片电影观看| 自线自在国产av| 成人手机av| a 毛片基地| 亚洲精品美女久久久久99蜜臀 | 成年动漫av网址| av免费观看日本| 一级黄片播放器| 一级毛片黄色毛片免费观看视频| 国产精品一区二区在线不卡| 99久久人妻综合| 欧美亚洲 丝袜 人妻 在线| 日韩欧美精品免费久久| 久久av网站| 国产男人的电影天堂91| 999精品在线视频| 人成视频在线观看免费观看| 精品国产国语对白av| 亚洲人成77777在线视频| 日韩欧美精品免费久久| 国产福利在线免费观看视频| 亚洲五月色婷婷综合| 亚洲av国产av综合av卡| 美女视频免费永久观看网站| 免费在线观看黄色视频的| 亚洲美女搞黄在线观看| 一本—道久久a久久精品蜜桃钙片| 男女午夜视频在线观看| 国产精品秋霞免费鲁丝片| 999精品在线视频| 亚洲精品,欧美精品| 欧美人与善性xxx| 亚洲人成电影观看| 不卡av一区二区三区| 人人妻人人爽人人添夜夜欢视频| 99香蕉大伊视频| 九色亚洲精品在线播放| 免费女性裸体啪啪无遮挡网站| 日韩成人av中文字幕在线观看| 亚洲av综合色区一区| 肉色欧美久久久久久久蜜桃| 久久久久久久精品精品| 两个人免费观看高清视频| 久久精品国产自在天天线| 韩国av在线不卡| 男女免费视频国产| 午夜91福利影院| 看免费成人av毛片| 国产精品女同一区二区软件| 国产精品二区激情视频| 丝袜美足系列| 国产成人精品久久二区二区91 | 国产精品久久久久久精品电影小说| 日韩av在线免费看完整版不卡| 国产精品一二三区在线看| av免费观看日本| 久久精品国产亚洲av天美| 美女国产视频在线观看| 久久久国产欧美日韩av| 丰满乱子伦码专区| 久久久久国产网址| 午夜老司机福利剧场| 亚洲一码二码三码区别大吗| 一区二区三区激情视频| 精品国产一区二区久久| 亚洲三级黄色毛片| 两个人看的免费小视频| 国产一区二区激情短视频 | 一级,二级,三级黄色视频| 热re99久久精品国产66热6| 最新中文字幕久久久久| 亚洲第一av免费看| 成人漫画全彩无遮挡| 黄色一级大片看看| 超碰成人久久| 只有这里有精品99| 国产深夜福利视频在线观看| 国产亚洲欧美精品永久| tube8黄色片| 国产97色在线日韩免费| 国产精品人妻久久久影院| 成年女人毛片免费观看观看9 | 成年女人在线观看亚洲视频| 美女国产视频在线观看| 久久久a久久爽久久v久久| 国产精品人妻久久久影院| 欧美少妇被猛烈插入视频| 九九爱精品视频在线观看| 一区二区av电影网| 自线自在国产av| 国产精品秋霞免费鲁丝片| 黑人巨大精品欧美一区二区蜜桃| 免费在线观看黄色视频的| 18禁观看日本| 五月伊人婷婷丁香| 成人二区视频| 日产精品乱码卡一卡2卡三| av女优亚洲男人天堂| 天堂中文最新版在线下载| 2021少妇久久久久久久久久久| 国产av精品麻豆| 国精品久久久久久国模美| 欧美日本中文国产一区发布| 日韩一本色道免费dvd| 欧美日韩国产mv在线观看视频| 大香蕉久久网| 啦啦啦中文免费视频观看日本| 高清视频免费观看一区二区| 在线观看国产h片| 久久久精品免费免费高清| 亚洲精品国产色婷婷电影| 日韩在线高清观看一区二区三区| 久久综合国产亚洲精品| 精品人妻熟女毛片av久久网站| 亚洲,一卡二卡三卡| www.av在线官网国产| 啦啦啦在线免费观看视频4| 日韩熟女老妇一区二区性免费视频| 王馨瑶露胸无遮挡在线观看| 国产白丝娇喘喷水9色精品| 亚洲欧美一区二区三区久久| 欧美日韩一区二区视频在线观看视频在线| 精品国产国语对白av| 波野结衣二区三区在线| 亚洲欧美一区二区三区久久| 亚洲经典国产精华液单| 99九九在线精品视频| 女人久久www免费人成看片| 高清不卡的av网站| 老汉色av国产亚洲站长工具| 久久久久久久久免费视频了| 免费看不卡的av| 一级毛片电影观看| 国产午夜精品一二区理论片| 夫妻午夜视频| 少妇被粗大的猛进出69影院| www.熟女人妻精品国产| 天天影视国产精品| 一级毛片电影观看| 亚洲,欧美,日韩| 国产成人a∨麻豆精品| av.在线天堂| 狂野欧美激情性bbbbbb| 制服诱惑二区| 伊人久久国产一区二区| 亚洲av福利一区| 亚洲精华国产精华液的使用体验| 久久青草综合色| 亚洲av电影在线进入| 免费日韩欧美在线观看| 亚洲av国产av综合av卡| 飞空精品影院首页| 女人精品久久久久毛片| 一区福利在线观看| 黄片播放在线免费| 国产精品一区二区在线观看99| 免费在线观看视频国产中文字幕亚洲 | 国产免费福利视频在线观看| 国产精品99久久99久久久不卡 | 久久久国产一区二区| a级片在线免费高清观看视频| 91在线精品国自产拍蜜月| 丝袜美腿诱惑在线| 亚洲欧美精品综合一区二区三区 | 日韩免费高清中文字幕av| 一级片免费观看大全| 麻豆乱淫一区二区| 亚洲色图 男人天堂 中文字幕| 熟妇人妻不卡中文字幕| 日韩一卡2卡3卡4卡2021年| 国产日韩欧美在线精品| 精品一区在线观看国产| 人妻一区二区av| 一二三四中文在线观看免费高清| 90打野战视频偷拍视频| 热99国产精品久久久久久7| 日韩av在线免费看完整版不卡| 婷婷色av中文字幕| 日本wwww免费看| 国产免费又黄又爽又色| 伦理电影免费视频| 九草在线视频观看| 午夜日韩欧美国产| 亚洲精品aⅴ在线观看| 一级黄片播放器| 日韩中字成人| 波多野结衣av一区二区av| 久久久精品94久久精品| 久久久久精品人妻al黑| 一级毛片我不卡| 国产成人精品久久二区二区91 | 国产成人免费无遮挡视频| 美女脱内裤让男人舔精品视频| 免费日韩欧美在线观看| 麻豆av在线久日| 麻豆乱淫一区二区| 日本午夜av视频| 老女人水多毛片| 国产亚洲最大av| 欧美+日韩+精品| av卡一久久| 亚洲欧美精品自产自拍| 午夜激情av网站| 亚洲色图 男人天堂 中文字幕| 久久久久久久国产电影| 亚洲国产av新网站| 欧美精品高潮呻吟av久久| 色哟哟·www| a级片在线免费高清观看视频| 亚洲,一卡二卡三卡| 夫妻性生交免费视频一级片| 日韩免费高清中文字幕av| 国产精品 国内视频| 亚洲国产欧美日韩在线播放| 日本vs欧美在线观看视频| 国产一区有黄有色的免费视频| 亚洲内射少妇av| 两个人免费观看高清视频| 91在线精品国自产拍蜜月| 寂寞人妻少妇视频99o| 麻豆乱淫一区二区| 日韩伦理黄色片| av在线app专区| av网站在线播放免费| 亚洲欧美成人精品一区二区| 高清不卡的av网站| 午夜福利在线观看免费完整高清在| 欧美精品亚洲一区二区| 国产 一区精品| 国产野战对白在线观看| 亚洲欧洲日产国产| 国产一区二区三区综合在线观看| 夜夜骑夜夜射夜夜干| av福利片在线| 国产成人午夜福利电影在线观看| 2021少妇久久久久久久久久久| 视频区图区小说| 高清视频免费观看一区二区|