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

    鹽堿脅迫下不同抗性楊樹品種轉(zhuǎn)錄組分析

    2025-03-13 00:00:00李文穎李曉宇楊成超
    中南林業(yè)科技大學學報 2025年2期
    關(guān)鍵詞:差異基因楊樹

    摘 要:【目的】研究不同抗性楊樹品種響應鹽堿脅迫的轉(zhuǎn)錄水平,挖掘楊樹耐鹽堿相關(guān)基因,為后續(xù)深入探討楊樹響應鹽堿脅迫的分子機制及后續(xù)遺傳轉(zhuǎn)化提供理論支撐。【方法】以楊樹耐鹽堿品種‘遼胡1號’和鹽堿敏感品種‘中遼1號’為材料,采用盆栽法,加入 NaHCO3溶液模擬鹽堿脅迫,利用轉(zhuǎn)錄組測序技術(shù)、生物信息學等方法分析鹽堿脅迫下不同抗性楊樹品種的轉(zhuǎn)錄水平,篩選與鹽堿脅迫相關(guān)的差異表達基因,隨機選取6個差異表達基因進行實時熒光定量PCR試驗,驗證轉(zhuǎn)錄組數(shù)據(jù)的準確性?!窘Y(jié)果】轉(zhuǎn)錄組測序結(jié)果表明,鹽堿脅迫下不同抗性楊樹品種之間共獲得5 273個差異表達基因,其中上調(diào)表達2 709個,下調(diào)表達2 564個。GO分析結(jié)果發(fā)現(xiàn),共有4 203個差異基因被注釋到3大類37條通路中。KEGG分析結(jié)果發(fā)現(xiàn),共有3 616個差異基因被注釋,差異基因顯著富集在植物信號轉(zhuǎn)導途徑(植物與病原體相互作用、植物激素信號轉(zhuǎn)導)、代謝途徑(淀粉與蔗糖代謝、碳代謝)和生物合成途徑(苯丙素的生物合成、氨基酸的生物合成)。鹽堿脅迫下MYB、ABC、WRKY、bHLH是差異基因較多的基因家族,篩選出CBL4、SOD2、bHLH35、ABCC10、CYP736A12、P5CS等差異基因表達顯著上調(diào)。qRT-PCR結(jié)果驗證轉(zhuǎn)錄組數(shù)據(jù)結(jié)果可靠?!窘Y(jié)論】本研究通過轉(zhuǎn)錄組測序技術(shù)對鹽堿脅迫下不同抗性楊樹品種進行比較分析,挖掘出與耐鹽堿性相關(guān)差異表達基因,為后續(xù)基因克隆和遺傳轉(zhuǎn)化奠定理論基礎(chǔ)。

    關(guān)鍵詞:楊樹;鹽堿脅迫;轉(zhuǎn)錄組測序;差異基因

    中圖分類號:S792.11 文獻標志碼:A 文章編號:1673-923X(2025)02-0122-09

    基金項目:國家科技創(chuàng)新2030-重大項目(2023ZD0405603);遼寧省農(nóng)業(yè)科學院院長基金項目(2023QN2416)。

    Transcriptome analysis of different resistant poplar varieties under salt alkali stress

    LI Wenying, LI Xiaoyu, YANG Chengchao

    (Poplar Research Institute of Liaoning Province, Gaizhou 115200, Liaoning, China)

    Abstract:【Objective】To study the transcriptional levels of different resistant poplar varieties in response to salt alkali stress, explore the genes related to salt alkali tolerance in poplar, and provide theoretical support for further exploration of the molecular mechanisms of poplar response to salt alkali stress and subsequent genetic transformation.【Method】Using poplar salt alkali tolerant variety P. simonii×P.euphratica and salt alkali sensitive variety P.×canadensiscv.‘Zhongliao1’as materials, using the potted plant method, NaHCO3 was added to simulate salt alkali stress. Transcriptome sequencing technology, bioinformatics, and other techniques were used to analyze the transcription levels of different resistant poplar varieties under salt alkali stress. Differential expression genes related to salt alkali stress were screened, and six differentially expressed genes were randomly selected for real-time fluorescence quantitative PCR testing to verify the accuracy of transcriptome data.【Result】Through transcriptome sequencing, it was found that under salt alkali stress, different resistant poplar varieties obtained a total of 5 273 differentially expressed genes, of which 2 709 were upregulated and 2 564 were downregulated. GO analysis results revealed that a total of 4 203 differentially expressed genes were annotated into 37 pathways across three major categories. The KEGG analysis results showed that a total of 3616 differentially expressed genes were annotated, and the differentially expressed genes were significantly enriched in plant signal transduction pathways (plant pathogen interaction, plant hormone signal transduction), metabolic pathways (starch and sucrose metabolism, carbon metabolism), and biosynthetic pathways (phenylpropanoid biosynthesis, amino acid biosynthesis). MYB, ABC, WRKY, and bHLH are gene families with more differentially expressed genes under salt alkali stress. CBL4, SOD2, bHLH35, ABCC10, CYP736A12, P5CS and other differentially expressed genes were significantly upregulated through screening. The qRT PCR results validate the reliability of transcriptome data.【Conclusion】This study compared and analyzed different resistant poplar varieties under salt alkali stress using transcriptome sequencing technology, and identified differentially expressed genes related to salt alkali tolerance, laying a theoretical foundation for subsequent gene cloning and genetic transformation.

    Keywords: poplar tree; salt alkali stress; transcriptome sequencing; differentially expressed genes

    由于自然因素如氣候干旱,以及人為因素如灌溉方式不合理、排水不暢、大量施用化肥等,使得鹽堿地面積日益擴大[1]。鹽堿脅迫是一種非常復雜且嚴重的非生物脅迫,當土壤含鹽量達到0.3%以上時,土壤中的植物就會遭受傷害,嚴重影響植物的生長和發(fā)育[2]。據(jù)統(tǒng)計,目前全球鹽堿地面積約為9.5×108 hm2,我國鹽堿地面積近1億hm2,約占全球鹽堿地面積的10%[3],且分布較廣,從東北松嫩平原、東部濱海沿岸到西北內(nèi)陸干旱、半干旱地區(qū)均有分布[4],特別是東北松嫩平原是我國北方土壤鹽堿化最廣泛和最嚴重的地區(qū)。其中堿性鹽(NaHCO3)所造成的土壤鹽堿化要比中性鹽(NaCl)造成的土壤鹽堿化問題更加嚴重[5],導致生態(tài)系統(tǒng)的平衡被破壞,森林資源和農(nóng)作物的產(chǎn)量降低,對人類的生產(chǎn)和生活造成了嚴重的不良影響。

    楊樹Populus L.具有生長速度快、易成活等特點,可以適應不同的氣候和土壤環(huán)境,在干旱、貧瘠、鹽漬的環(huán)境下也可以正常生長,是我國北方平原地區(qū)和沙區(qū)營造防護林和用材林的主要樹種,在生態(tài)林的建設(shè)中占有重要地位[4-5]。隨著高通量測序技術(shù)的發(fā)展,轉(zhuǎn)錄組學研究在植物相關(guān)機理研究、優(yōu)異功能基因的挖掘、分子標記開發(fā)等方面已有較廣泛的應用[6]。李菡[7]在對鹽堿脅迫下不同耐鹽性小麥進行轉(zhuǎn)錄組測序時,挖掘出38個表達趨勢相反的差異表達基因,這些基因可能是響應小麥鹽堿脅迫的關(guān)鍵基因,為小麥抗逆育種及鹽堿地小麥栽培管理提供了理論依據(jù)。呂博文等[8]在對蕨麻Potentilla anserine響應鹽堿脅迫的形態(tài)及轉(zhuǎn)錄組分析時發(fā)現(xiàn),ERF、MYB、NAC、WRKY、MYB和bZIP是差異基因最多的轉(zhuǎn)錄因子家族。在鑒定油茶Camellia oleifera中WRKY基因家族在高鹽脅迫下的表達模式發(fā)現(xiàn),CoWRKY11、CoWRKY14、CoWRKY20、CoWRKY29 和CoWRKY56基因在高鹽脅迫下均被誘導表達[9]。目前關(guān)于鹽堿脅迫下楊樹的響應分子機制已有研究報道,如84k楊PagERF114[10]、銀中楊Populus alba × P. Berolinensis的PabCIPK24等基因均被證實參與了植株的耐鹽性[11],這些研究結(jié)果和數(shù)據(jù)都為探究楊樹在響應鹽堿脅迫的分子機制提供了理論基礎(chǔ)。

    本研究對鹽堿脅迫下不同抗性楊樹品種進行轉(zhuǎn)錄組測序,分析差異基因變化差異,挖掘楊樹耐鹽堿相關(guān)基因,為深入研究楊樹耐鹽堿的分子機制提供理論依據(jù)。

    1 材料與方法

    1.1 試驗材料

    試驗地選在遼寧省蓋州市,以楊樹耐鹽堿品系‘遼胡1號楊’P. simonii×P.euphratica(以下用LH表示)和鹽堿敏感品系‘中遼1號楊’P.×canadensis cv. ‘Zhongliao1’(以下用ZL表示)為材料,采用盆栽方法進行試驗。選取當年生扦插苗,截取直徑 1.6 cm、長20 cm左右的穗材,于2023年4月扦插到裝有河砂與珍珠巖比例為5∶1的花盆中,扦插后每個花盆重4.5 kg,加入60 mmol/L NaHCO3溶液模擬鹽脅迫,以加蒸餾水為對照(CK),每4 d澆一次鹽溶液和蒸餾水,37 d后分別采集同一高度位置新鮮葉片樣品,設(shè)3次獨立生物學重復,裝至凍存管于液氮速凍后,置于超低溫冰箱(-80 ℃)進行保存,分別記為ZL60-1,ZL60-2,ZL60-3和LH60-1,LH60-2,LH60-3。

    1.2 試驗方法

    1.2.1 RNA提取和檢測

    委托百邁客(青島)有限公司完成轉(zhuǎn)錄組測序。采用試劑盒方法提取樣品總RNA,使用NanoDrop 2000分光光度計和Agient2100生物分析儀檢測提取RNA的濃度、純度(OD260/OD280)和完整性。

    1.2.2 文庫構(gòu)建和測序

    樣品檢測合格后,用帶有Oligo(dT)的磁珠富集真核生物mRNA,隨后片段化,以mRNA為模板,合成cDNA鏈并進行純化;純化的cDNA再進行末端修復、加A尾并連接測序接頭,最后通過PCR富集得到cDNA文庫,文庫檢測合格后,利用Illumina高通量測序平臺對cDNA文庫進行測序。

    1.2.3 數(shù)據(jù)質(zhì)控和序列比對

    為了保證數(shù)據(jù)分析的質(zhì)量及可靠性,使用fastp軟件對測序所得到的原始數(shù)據(jù)(Raw data)進行過濾質(zhì)控,得到干凈數(shù)據(jù)(Clean data),去除含有接頭和低質(zhì)量的數(shù)據(jù)(Reads),再經(jīng)原始數(shù)據(jù)過濾、測序錯誤率檢測、GC含量分布檢查,最終得到高質(zhì)量的有效數(shù)據(jù)(Clean reads)。使用HISAT2軟件將Clean reads比對到已公開發(fā)表的毛果楊Populus trichocarpa參考基因組上,獲取reads在參考基因組上的定位信息[12]。

    1.2.4 差異表達基因篩選

    新轉(zhuǎn)錄本組裝后,使用TPM(Transcripts per million)法計算基因表達量,利用DESeq軟件對基因進行差異分析,差異表達基因檢測過程中,通過統(tǒng)計學模型進行假設(shè)檢驗概率的計算,對差異顯著性P值進行校正得到錯誤發(fā)現(xiàn)率(False discovery rate,F(xiàn)DR),將差異倍數(shù)Fold change(Fc)≥2且FDR<0.01作為差異基因的篩選標準。為了方便比較,將差異倍數(shù)取對數(shù)值,以log2Fc表示,log2Fc的絕對值越大,F(xiàn)DR值越小的基因在2組樣本中的差異變化越明顯。將篩選得到的差異表達GO、KEGG數(shù)據(jù)庫中進行注釋[21],分析差異表達基因在不同數(shù)據(jù)庫中的分布,以圖形和表格的形式展示不同功能的差異基因數(shù)目及其生物學功能。

    1.2.5 差異表達基因的熒光定量PCR驗證

    轉(zhuǎn)錄組測序分析完成后,隨機挑選6個差異表達基因進行熒光定量驗證。以UBQ基因為內(nèi)參基因,使用Primer5軟件設(shè)計特異性引物(表1),按如下反應體系進行配制,依次加入PowerUp SYBR Green Master Mix 5 μL,上游引物(Primer F)0.2 μL,下游引物(Primer R)0.2 μL,cDNA 1 μL,加ddH2O水至10 μL。PCR反應體系為:95 ℃ 5 min;95 ℃ 15 s,60 ℃ 20 s,72 ℃ 20 s,40個循環(huán);95 ℃ 15 s,65 ℃ 1 min。將各樣品的目的基因和內(nèi)參分別進行 Real Time PCR 反應,每個樣本檢測3個復孔,數(shù)據(jù)采用2-△△CT法進行分析。

    2 結(jié)果與分析

    2.1 提取RNA濃度檢測

    提取的樣本RNA濃度在504~838 ng/μL,OD260/OD280值在2.13~2.14,完整性(RIN)均值為7.51,檢測結(jié)果均為A,綜上結(jié)果表明提取RNA 濃度較高,完整性好,符合質(zhì)量要求(表2)。

    2.2 轉(zhuǎn)錄組測序數(shù)據(jù)質(zhì)控結(jié)果分析

    對原始序列進行過濾和質(zhì)控處理后,各處理樣本有效數(shù)據(jù)均達到5.81 G,質(zhì)控過程中GC含量保持在44%左右,Q30堿基百分比在92.79%及以上。通過將reads數(shù)據(jù)比對到參考基因組上,獲得各個樣本的基因組比對情況,比對比例在85%以上,表明測序質(zhì)量良好(表3)。

    2.3 差異表達基因在不同數(shù)據(jù)庫分布

    以Fold change(Fc)≥2且FDR<0.01作為差異基因的篩選標準,統(tǒng)計‘中遼1號’和‘遼胡1號’在NaHCO3脅迫后的差異表達基因的表達情況,共有5 273個差異表達基因,其中上調(diào)表達2 709個,下調(diào)表達2 564個。差異基因在各數(shù)據(jù)庫中的注釋情況如表4所示,在各數(shù)據(jù)庫中均有注釋,其中在NR(5 134條)、GO(4 203條)、KEGG(3 616條)數(shù)據(jù)庫中差異基因分布較多。

    2.4 差異表達基因GO注釋分析

    通過對差異基因進行GO注釋分析,共有4 203條差異基因分別注釋到生物學過程(Biological process)、細胞組分(Cellular component)和分子功能(Molecular function)3大類的37個功能組中。在生物學過程中涉及的差異基因共分成21個功能組,其中差異基因在細胞過程(Cellular process)、代謝過程(Metabolic process)、生物調(diào)節(jié)(Biological regulation),對刺激的響應(Response to stimulus)等顯著富集。細胞組分中涉及的差異基因分成3個功能組,富集在細胞解剖實體(Cellular anatomical entity)、細胞內(nèi)(Intracellular)和蛋白質(zhì)復合物(Protein-containing complex)。分子功能中涉及的差異基因共分成13個功能組,其中差異基因在結(jié)合物(Binding)和催化活性(Catalytic activity)中顯著富集(圖1)。

    2.5 差異表達基因KEGG分析

    通過對差異基因進行KEGG注釋分析,共有3 616個差異基因定位到131個具體的代謝途徑分支中,其中差異基因在信號轉(zhuǎn)導途徑中富集較多,分別是植物與病原體相互作用(228條)和植物激素信號轉(zhuǎn)導途徑(150條),在代謝途徑和生物合成中也有富集,分別是淀粉與蔗糖代謝(104條)和碳代謝(78條)途徑,生物合成中的苯丙素(85條)和氨基酸的生物合成(72條)也顯著富集(表5)。

    2.6 鹽堿脅迫下基因家族分析

    在楊樹響應鹽堿脅迫時,MYB、bHLH、WRKY、bZIP、NAC、PYL和ABC等基因家族都發(fā)揮了作用,其中注釋較多的是MYB、ABC、WRKY和bHLH基因家族(表6)。

    2.7 鹽堿脅迫相關(guān)基因篩選

    篩選出楊樹鹽堿脅迫相關(guān)的差異表達基因,其中CBL4、SOD2、bHLH35、ABCC10、CYP736A12、WRKY72、WRKY31、P5CS、bZIP11、DIV、ABCB11等基因在楊樹耐鹽堿品種中表達上調(diào),ABCG11、PYL4、NAC29基因在耐鹽堿品種中表達下調(diào)(表7)。這些基因可能是參與楊樹耐鹽堿性的關(guān)鍵基因。

    2.8 差異表達基因qRT-PCR驗證

    以鹽堿脅迫下的楊樹葉片轉(zhuǎn)錄組測序的cDNA為模板,對隨機選取的6個差異表達基因進行qRT-PCR驗證,熒光定量結(jié)果與轉(zhuǎn)錄組測序分析結(jié)果基本一致,進一步證明轉(zhuǎn)錄組測序結(jié)果的可靠性(圖2)。

    3 討 論

    土壤鹽堿化是影響植物生長和產(chǎn)量的一個世界難題,因此,探討植物耐鹽堿的分子機制具有重要的意義[13]。近年來,植物響應鹽堿脅迫的機制已成為不少學者研究的熱點,對大豆Glycine max[14]、駿棗[15]、藜麥Chenopodium quinoa[16]等植物分別從生長特性、生理指標、葉片結(jié)構(gòu)、轉(zhuǎn)錄組、代謝組等方面對植物的耐鹽堿機制進行了研究,加深了人們對不同植物耐鹽機制的了解。

    目前,關(guān)于楊樹耐鹽堿的分子機制研究較少,本研究以楊樹抗鹽堿性品種‘遼胡1號楊’和鹽堿敏感品種‘中遼1號楊’為材料,進行轉(zhuǎn)錄組測序,結(jié)果表明,鹽堿脅迫下2個楊樹品種共得到5 273條差異表達基因,其中表達上調(diào)2 709條,表達下調(diào)2 564條,對這些差異基因進行GO和KEGG富集分析,明確差異基因的生物功能和信號通路。GO功能注釋結(jié)果表明,差異基因在生物過程中的細胞過程和代謝過程顯著富集,在分子功能方面,差異基因在結(jié)合物和催化活性顯著富集。KEGG代謝通路分析發(fā)現(xiàn),差異基因富集到苯丙素生物合成、植物激素轉(zhuǎn)導、淀粉和蔗糖代謝等通路。這些通路可能是影響楊樹耐鹽堿性的重要因素,這與在對大豆耐鹽堿轉(zhuǎn)錄組測序結(jié)果一致[17]。碳代謝途徑在差異基因中也顯著富集,碳代謝包括植物呼吸作用和光合作用,為植物生長發(fā)育提供能量,有研究發(fā)現(xiàn)在50 mmol/L鹽濃度下,胡楊、新疆楊和小葉楊和胡楊的雜交種的光合呼吸作用會受到一定的促進作用[18],以上結(jié)果表明這些功能和途徑在楊樹響應鹽堿脅迫中發(fā)揮重要作用。

    細胞色素P450參與植物體內(nèi)內(nèi)源激素合成以及次生代謝[19] ,有研究發(fā)現(xiàn)人參Panax ginseng中細胞色素P450編碼的PgCYP736A12基因在Nacl脅迫后表達上調(diào),該基因參與了除草劑的代謝[ 20]。本研究中CYP736A12基因表達上調(diào),表明該基因在鹽堿脅迫下發(fā)揮一定的作用。CBL蛋白被鑒定為擬南芥Arabidopsis thaliana和水稻Oryza sativa鈣信號通路中的鈣離子傳感器,由CBL和其受體CIPKs形成的復雜信號網(wǎng)絡系統(tǒng),在鈣離子的介導下能夠使植物遭受逆境脅迫時做出應答[21]。有研究表明,楊樹在遭受寒冷、干旱、鹽脅迫等非生物脅迫時CBL基因家族發(fā)揮重要的調(diào)節(jié)作用,在煙草中過表達MsCBL4能夠提高轉(zhuǎn)基因株系的耐鹽堿能力[22];瞬時過表達ThCBL4基因提高剛毛檉柳Tamarix hispida耐鹽能力[23];本研究中,PsCBL4在耐鹽堿性品種中表達上調(diào),推測該基因在楊樹耐鹽堿性上發(fā)揮積極作用,但其具體功能還需進一步研究。ABCB和ABCC都是ABC轉(zhuǎn)運蛋白家族的一個亞家族蛋白,參與植物體內(nèi)的跨膜轉(zhuǎn)運過程。水稻中的OsABCB11的基因表達受到 IAA、PEG及NaCl的誘導,在生長素及響應脅迫的過程中均發(fā)揮作用[24];有研究發(fā)現(xiàn)在對水稻進行鹽誘導后,OsABCC10基因隨鹽脅迫濃度的升高及脅迫時間的延長,表達量也逐漸增加,并鑒定出OsABCC10是耐鹽基因[25]。本研究中,PsABCC10基因在耐鹽性強品種中表達上調(diào),推測該基因是楊樹響應鹽堿脅迫的關(guān)鍵基因。P5CS基因是植物體內(nèi)脯氨酸生物合成的關(guān)鍵酶,過表達P5CS 基因能有效提高植物的耐鹽性,有研究發(fā)現(xiàn)將大豆中的P5CS基因成功轉(zhuǎn)入水稻,并在NaCl脅迫后,轉(zhuǎn)入P5CS基因的水稻中的脯氨酸含量明顯提高,且耐鹽性得到增強[26]。本研究中PsP5CS基因在耐鹽性強品種中表達上調(diào),推測在耐鹽堿性上發(fā)揮重要作用,可為后續(xù)挖掘利用P5CS基因并深入研究其功能奠定基礎(chǔ)。WRKY、bHLH、bZIP基因家族參與植物生長發(fā)育、組織分化等生物過程,在植物響應鹽堿脅迫的過程中也發(fā)揮重要作用[27]。本研究中PsWRKY72基因在耐鹽性強楊樹品種中表達上調(diào),推測該基因是耐鹽堿的關(guān)鍵基因,有研究證實了這一點,對不同鹽脅迫下甘薯Ipomoea batatas進行轉(zhuǎn)錄組測序,挖掘到耐鹽相關(guān)基因 IbWRKY72,并對其表達模式進行了分析發(fā)現(xiàn)鹽脅迫誘導IbWRKY72 在甘薯根部的表達[28]。本研究中PsWRKY31在耐鹽堿性強楊樹品種中上調(diào)表達,有研究將歐美楊Populus×euramericana中WRKY31基因轉(zhuǎn)入到煙草中能增強植物的耐鹽能力[29],可為后續(xù)遺傳轉(zhuǎn)化奠定基礎(chǔ)。PsbHLH35基因在楊樹鹽堿脅迫下也發(fā)揮關(guān)鍵作用,在本研究中,在鹽脅迫下耐鹽堿強品種中PsbHLH35表達量上調(diào),bHLH35是一種鹽誘導的基因,主要集中在發(fā)芽種子和幼苗。OsbHLH35分別通過ABA依賴性和ABA非依賴性途徑介導發(fā)芽,減輕鹽脅迫下 ABA 對水稻發(fā)芽的抑制作用[30]。由此推測該基因具有調(diào)控植物抗鹽脅迫的功能,可為后續(xù)培育轉(zhuǎn)基因植株提供理論基礎(chǔ)。有研究對轉(zhuǎn)MsbZIP11基因下的擬南芥在抵御鹽堿脅迫方面的功能進行研究,發(fā)現(xiàn)在鹽誘導下轉(zhuǎn)MsbZIP11基因擬南芥相對野生型擬南芥根系更長,存活率更高[31]。本研究中PsbZIP11基因在鹽堿性強的品種中表達量高,推測PsMsbZIP11基因在楊樹響應鹽脅迫的應激反應過程中起到關(guān)鍵性作用。ABCG11[32]、PYL4[33]、NAC29[34]基因在植物響應鹽堿脅迫時也發(fā)揮重要作用。

    本研究缺少對與鹽堿脅迫相關(guān)基因的生物學功能的深入研究,尚不清楚上述基因參與鹽堿脅迫的分子調(diào)控網(wǎng)絡。在后續(xù)研究中,將通過耐鹽堿相關(guān)基因克隆、功能驗證和遺傳轉(zhuǎn)化等試驗,探究上述基因在提升楊樹抗鹽堿能力的分子機制。

    4 結(jié) 論

    本研究通過轉(zhuǎn)錄組測序技術(shù)對鹽堿脅迫下不同抗性楊樹品種進行比較分析,獲得了質(zhì)量較好的轉(zhuǎn)錄組數(shù)據(jù)和差異表達基因,并且篩選出與鹽堿脅迫相關(guān)的差異表達基因,為后續(xù)深入分析楊樹耐鹽堿的分子機理和遺傳轉(zhuǎn)化奠定了基礎(chǔ)。

    參考文獻:

    [1] 朱建峰,崔振榮,吳春紅,等.我國鹽堿地綠化研究進展與展望[J].世界林業(yè)研究,2018,31(4):70-75. ZHU J F, CUI Z R, WU C H, et al. Research advances and prospect of saline and alkali land greening in China[J]. World Forestry Research,2018,31(4):70-75.

    [2] 孫羽,王麒,宋秋來,等.寒地水稻對鹽堿脅迫的響應及耐鹽堿生理機制研究進展[J].黑龍江農(nóng)業(yè)學,2022(8):57-61. SUN Y, WANG Q, SONG Q L, et al. Research progress on response of rice to salin-alkali stress and physiological mechanism of saline-alkali tolerance in cold region[J]. Heilongjiang Agricultural Sciences,2022(8):57-61.

    [3] 王才林,張亞東,趙凌,等.耐鹽堿水稻研究現(xiàn)狀、問題與建議[J].中國稻米,2019,25(1):1-6. WANG C L, ZHANG Y D, ZHAO L, et al. Research status, problems, and suggestions on salt alkali tolerant rice[J]. China Rice,2019,25(1):1-6.

    [4] 費小鈺.鹽堿脅迫對轉(zhuǎn)Lc-CDPK基因水稻生理生化特性的影響[D].長春:吉林農(nóng)業(yè)大學,2017. FEI X Y. Effects of saline-alkali stress on the physiological and biochemical characteristics of transgenic Lc-CDPK genes rice[D]. Changchun:Jilin Agricultural University,2017.

    [5] 崔曉濤.新西伯利亞銀白楊耐干旱和耐鹽堿能力研究[D].哈爾濱:東北林業(yè)大學,2009. CUI X T. Saline-alkali and drought tolerance of Populus bachofenii ×P. pyramidalis ROZ[D]. Harbin:Northeast Forestry University,2009.

    [6] 慧芳,劉秀巖,李宗諭,等.轉(zhuǎn)錄組測序技術(shù)在藥用植物研究中的應用[J].中草藥,2019,50(24):6149-6155. HUI F, LIU X Y, LI Z Y, et al. Application of transcriptome sequencing in study of medicinal plants[J]. Chinese Traditional and Herbal Drugs,2019,50(24):6149-6155.

    [7] 李菡.小麥苗期響應鹽堿脅迫的生理和基因轉(zhuǎn)錄表達分析[D].秦皇島:河北科技師范學院,2023. LI H. Physiological and gene transcriptional expression analysis of wheat seedling stage in response to saline-alkali stress[D]. Qinhuangdao:Hebei Normal University of Science Technology, 2023.

    [8] 呂博文,李軍喬,閆格,等.蕨麻響應鹽堿脅迫的形態(tài)及轉(zhuǎn)錄組分析[J].種子,2024,43(1):14-28. LYU B W, LI J Q, YAN G, et al. Morphological and transcriptome analysis on response of Potentilla anserina to salt-alkali stress[J]. Seed,2024,43(1):14-28.

    [9] 蘇文娟,曹瑞蘭,周增亮,等.油茶WRKY基因家族鑒定及逆境脅迫表達分析[J].中南林業(yè)科技大學學報,2023,43(3): 155-166,174. SU W J, CAO R L, ZHOU Z L, et al. Identification and expression analysis of the WRKY gene family in Camellia oleifera under stress conditions[J]. Journal of Central South University of Forestry Technology,2023,43(3):155-166,174.

    [10] 劉森堯.84K楊PagERF114基因耐鹽功能分析[D].哈爾濱:東北林業(yè)大學,2023. LIU S Y. Analysis of salt tolerance function of PagERF114 gene in Populus alba×P.grandulos[D]. Harbin:Northeast Forestry University,2023.

    [11] 白向東.銀中楊PabCIPK24基因耐鹽研究[D].哈爾濱:東北林業(yè)大學,2023. BAI X D. Salt tolerance of PabCIPK24 in Populus alba × P. Berolinensis [D]. Harbin:Northeast Forestry University,2023.

    [12] KIM D, LANGMEAD B, SALZBERG S. HISAT: A fast spliced aligner with low memory requirements[J]. Nat Methods, 2015,12(4):357-360.

    [13] YE X, WANG H, CAO X, et al. Transcriptome profiling of Puccinellia tenuiflora during seed germination under a long-term saline-alkali stress[J]. BMC Genomics,2019,20(1):1-17.

    [14] 王春雨,朱冠雄,田藝心,等.大豆耐鹽堿生理機制及種質(zhì)篩選研究進展[J].大豆科學,2024,43(1):107-113. WANG C Y, ZHU G X, TIAN Y X, et al. Research progress on physiological mechanism of salt alkali tolerance and germplasm screening in soybean[J]. Soybean Science,2024,43(1):107-113.

    [15] 閆敏,王艷,王程成,等.混合鹽堿脅迫對駿棗葉片結(jié)構(gòu)和光合熒光特性的影響[J].中國農(nóng)業(yè)科技導報,2024,26(3):57-65. YAN M, WANG Y, WANG C C, et al. Effect of mixed salinealkali stress on leaf structure and photosynthetic fluorescence properties of jujube[J]. Journal of Agricultural Science and Technology,2024,26(3):57-65.

    [16] 鮑慶晗.藜麥響應鹽堿脅迫的生理特性及轉(zhuǎn)錄組與代謝組機制解析[D].呼和浩特:內(nèi)蒙古農(nóng)業(yè)大學,2023. BAO Q H. Physiological characteristics and mechanism analysis of transcriptome and metabolome of quinoa under salt and alkali stress[D]. Hohhot: Inner MongoliaAgricultural University,2023.

    [17] 張子呈.耐鹽堿大豆品種鑒定及標記基因篩選[D].哈爾濱:東北農(nóng)業(yè)大學,2024. ZHANG Z C. Identification of saline-alkali tolerant soybean varieties and screening of marker genes[D]. Harbin:Northeast Agricultural University,2024.

    [18] 張震中.幾種楊樹抗旱耐鹽特性研究及其對土壤微生物群落的影響[D].北京:北京林業(yè)大學,2021. ZHANG Z Z. Study on characteristics of drought and salt tolerance of several species of poplars and their effects on soil microbial communities[D]. Beijing:Beijing Forestry University, 2021.

    [19] XU J, WANG X, GUO W. The cytochrome P450 superfamily: Key players in plant development and defense[J]. Journal of Integrative Agriculture,2015,14(9):1673-1686.

    [20] KHANOM S, JANG J, LEE R. Overexpression of ginseng cytochrome P450 CYP736A12 alters plant growth and confers phenylurea herbicide tolerance in Arabidopsis[J]. Journal of Ginseng Research,2019,43(4):645-653.

    [21] KUDLA J, BATISTIC O, HASHIMOTO K. Calcium signals: the lead currency of plant information processing[J]. Plant Cell, 2010,22(3):541-563.

    [22] YI M A, XIAO X Y, LI S Z, et al. Alfalfa MsCBL4 enhances calcium metabolism but not sodium transport in transgenic tobacco under salt and saline-alkali stress[J]. Plant Cell Reports, 2020,39(8):997-1011.

    [23] 鄒全程,唐緋緋,劉中原,等.瞬時過表達ThCBL4基因提高剛毛檉柳耐鹽能力[J].林業(yè)科學研究,2018,31(3):60-67. ZOU Q C, TANG F F, LIU Z Y, et al. A novel Calcineurin B-like proteins gene (ThCBL4) improving the salt tolerance in transient overexpression Tamarix hispida[J]. Forest Research, 2018,31(3):60-67.

    [24] 賀禎媚.水稻OsABCB11基因功能研究[D].杭州:浙江大學, 2020. HE Z M. Research on function of OsABCB11 gene in rice[D]. Hangzhou:Zhejiang University,2020.

    [25] 陸優(yōu)社.水稻OsABCC10基因的克隆及功能分析[D].南寧:廣西大學,2019. LU Y S. Cloning and functional analysis of OsABCC10 in rice[D]. Nanning:Guangxi University,2019.

    [26] 李大紅,劉喜平,甄萍萍. NaCl 對轉(zhuǎn) P5CS 基因水稻種子萌發(fā)的影響[J].種子,2010,29(12):23-26. LI D H, LIU X P, ZHEN P P. The effect of NaCl on seed germination of the transgenic P5CS Rice[J]. Seed,2010,29(12): 23-26.

    [27] 張勇,沈丙權(quán),秦偉捷,等.糖基化蛋白質(zhì)組學:結(jié)構(gòu)、功能和研究方法[J].生命科學,2018,30(4):480-490. ZHANG Y, SHEN B Q, QIN W J, et al. Glycoproteomics: structure, function and methods[J]. Chinese Bulletin of Life Sciences,2018,30(4):480-490.

    [28] 張毅.基于轉(zhuǎn)錄組和共表達網(wǎng)絡分析挖掘甘薯花青素合成和耐鹽相關(guān)基因[D].晉中:山西農(nóng)業(yè)大學,2023. ZHANG Y. Exploring genes related with anthocyanin biosynthesis and salt tolerance in sweet potato based on analysis of transcriptome and co-expression network[D]. Jinzhong: Shanxi Agricultural University,2023.

    [29] 潘鈺.歐美楊轉(zhuǎn)錄因子PeWRKY31的克隆及功能鑒定[D].保定:河北農(nóng)業(yè)大學,2024. PAN Y. Cloning and functional identification of transcription factor PeWRKY31 in Populus×euramericana[D]. Baoding:Hebei Agricultural University,2024.

    [30] CHEN H C, CHENG W H, HONG C Y, et al. The transcription factor OsbHLH035 mediates seed germination and enables seedling recovery from salt stress through ABA-dependent and ABA-independent pathways, respectively[J]. Rice,2018,11(1):50.

    [31] 杜秉昊.紫花苜蓿MsbZIP11基因抗鹽堿的功能研究[D].哈爾濱:哈爾濱師范大學,2021. DU B H. Functional characterization of saline-alkali resistance of MSbZIP11 gene in alfalfa[D]. Harbin:Harbin Normal University, 2021.

    [32] MALIK A W. ABC轉(zhuǎn)運蛋白基因(ABC)和谷胱甘肽蛋白基因(GRX)在陸地棉鹽脅迫響應中的交互作用機制[D].北京:中國農(nóng)業(yè)科學院,2022. MALIK A W. Crosstalk of ABC-transporter and glutaredoxin genes to mediate salt stress in Gossypium hirsutum[D]. Beijing: Chinese Academy of Agricultural Sciences,2022.

    [33] 王宏,藺經(jīng),李曉剛,等.鹽脅迫下的杜梨PbPYL4基因克隆及其與PbNCED2基因表達分析[J].果樹學報,2014,31(6): 1017-1023. WANG H, LIN J, LI X G, et al. Molecular cloning of PbPYL4 gene and expression analysis of PbPYL4 and PbNCED2 in Pyrus betulaefolia under salt stress[J]. Journal of Fruit Science, 2014,31(6):1017-1023.

    [34] 丁海濱.蘋果屬山定子MbNAC25與MbNAC29基因的克隆與功能分析[D].哈爾濱:東北農(nóng)業(yè)大學,2019. DING H B. Isolation and functional analysis of Malus baccata( L.) Borkh NAC transcription factors MbNAC25 and MbNAC29[D]. Harbin: Northeast Agricultural University,2019.

    [本文編校:吳 彬]

    猜你喜歡
    差異基因楊樹
    膠質(zhì)母細胞瘤早期復發(fā)和假性進展的癌癥免疫差異表達基因分析及活性化合物預測
    ICR鼠肝和腎毒性損傷生物標志物的篩選
    基于RNA 測序研究人參二醇對大鼠心血管內(nèi)皮細胞基因表達的影響 (正文見第26 頁)
    基于GEO數(shù)據(jù)庫對骨關(guān)節(jié)炎滑膜差異基因的篩選及生物信息學分析
    跟蹤導練(五)
    口算本失蹤記
    楊樹山漫畫作品欣賞
    喜劇世界(2017年10期)2017-06-01 12:39:29
    楊樹山漫畫作品欣賞
    喜劇世界(2017年9期)2017-06-01 12:39:19
    霧霾來襲
    喜劇世界(2017年4期)2017-03-08 13:00:56
    SSH技術(shù)在絲狀真菌功能基因篩選中的應用
    亚洲美女搞黄在线观看| 99久久精品热视频| 国产高潮美女av| 色播亚洲综合网| 人人妻人人爽人人添夜夜欢视频 | 97精品久久久久久久久久精品| 国产毛片在线视频| 亚洲av福利一区| 又黄又爽又刺激的免费视频.| 久久久久久久久久久免费av| 久久久精品欧美日韩精品| 国产黄色视频一区二区在线观看| 国产亚洲精品久久久com| 中国美白少妇内射xxxbb| 免费高清在线观看视频在线观看| 亚洲国产色片| 成人免费观看视频高清| 极品教师在线视频| 亚洲成人久久爱视频| 一级毛片电影观看| 蜜臀久久99精品久久宅男| 亚洲欧美成人精品一区二区| 国产精品熟女久久久久浪| 永久网站在线| 亚洲人与动物交配视频| 亚洲精品日韩在线中文字幕| 国产乱来视频区| 国产精品无大码| 欧美成人精品欧美一级黄| 亚洲欧美日韩无卡精品| 国产一级毛片在线| 日韩av免费高清视频| av免费在线看不卡| 边亲边吃奶的免费视频| 国产又色又爽无遮挡免| 日韩一区二区三区影片| 欧美一区二区亚洲| 99精国产麻豆久久婷婷| 在线免费观看不下载黄p国产| 亚洲av中文av极速乱| 97精品久久久久久久久久精品| 国产精品成人在线| 亚洲国产成人一精品久久久| 亚洲av成人精品一二三区| 亚洲天堂国产精品一区在线| 国产国拍精品亚洲av在线观看| 99久久精品一区二区三区| 激情 狠狠 欧美| 亚洲国产日韩一区二区| 日韩欧美精品免费久久| 边亲边吃奶的免费视频| 99热这里只有是精品在线观看| 亚洲成人中文字幕在线播放| 熟妇人妻不卡中文字幕| 七月丁香在线播放| 成年女人在线观看亚洲视频 | 免费看日本二区| 亚洲怡红院男人天堂| 国产一区有黄有色的免费视频| 高清av免费在线| 日日摸夜夜添夜夜添av毛片| 日本爱情动作片www.在线观看| av福利片在线观看| 久久久久网色| 亚洲精品,欧美精品| 极品教师在线视频| 日韩三级伦理在线观看| 18+在线观看网站| 视频中文字幕在线观看| 亚洲av二区三区四区| 亚洲国产精品成人久久小说| 亚洲av不卡在线观看| 国产乱人视频| 国产精品久久久久久精品古装| 99久久精品热视频| 午夜免费观看性视频| 日韩成人av中文字幕在线观看| 免费av不卡在线播放| 久久久久精品性色| 亚洲天堂av无毛| 亚洲天堂国产精品一区在线| 97在线视频观看| 日韩,欧美,国产一区二区三区| 插阴视频在线观看视频| 男人舔奶头视频| 蜜桃久久精品国产亚洲av| 国产成年人精品一区二区| 18+在线观看网站| 国模一区二区三区四区视频| 亚洲精品自拍成人| 亚洲av在线观看美女高潮| 精品久久久久久久久av| 99re6热这里在线精品视频| 久久97久久精品| 91久久精品电影网| 干丝袜人妻中文字幕| 丰满人妻一区二区三区视频av| 亚洲人成网站在线观看播放| 欧美zozozo另类| av在线蜜桃| 2021少妇久久久久久久久久久| 亚洲av国产av综合av卡| 久久精品国产亚洲av天美| 国产一区二区亚洲精品在线观看| 国产精品.久久久| 熟女人妻精品中文字幕| 又大又黄又爽视频免费| 国产伦精品一区二区三区四那| 欧美激情在线99| 男的添女的下面高潮视频| 狠狠精品人妻久久久久久综合| 国产黄a三级三级三级人| 大片电影免费在线观看免费| 日韩中字成人| 国内少妇人妻偷人精品xxx网站| 国产精品一区二区三区四区免费观看| 久久这里有精品视频免费| 乱系列少妇在线播放| 亚洲av电影在线观看一区二区三区 | 日韩av在线免费看完整版不卡| 欧美日韩亚洲高清精品| 亚洲精品国产色婷婷电影| 天堂中文最新版在线下载 | 日韩电影二区| 欧美日韩国产mv在线观看视频 | 在线a可以看的网站| 国产69精品久久久久777片| 美女高潮的动态| 精品久久久精品久久久| 精品一区在线观看国产| 亚州av有码| 少妇人妻久久综合中文| 免费观看a级毛片全部| 国产免费又黄又爽又色| 日本一二三区视频观看| 黄片wwwwww| 美女被艹到高潮喷水动态| 丰满少妇做爰视频| 哪个播放器可以免费观看大片| 日日啪夜夜爽| 亚洲精品久久午夜乱码| 亚洲一区二区三区欧美精品 | 亚洲欧美日韩无卡精品| 久久久欧美国产精品| 高清视频免费观看一区二区| 一区二区三区乱码不卡18| 美女主播在线视频| 特级一级黄色大片| 麻豆精品久久久久久蜜桃| 97在线视频观看| 交换朋友夫妻互换小说| 黄色欧美视频在线观看| 日本一二三区视频观看| 一本色道久久久久久精品综合| 夫妻性生交免费视频一级片| 久久久久久久大尺度免费视频| 色视频在线一区二区三区| 午夜日本视频在线| 国产在线男女| 麻豆成人午夜福利视频| 中国三级夫妇交换| 亚洲色图av天堂| 如何舔出高潮| 激情五月婷婷亚洲| 亚洲经典国产精华液单| 内地一区二区视频在线| 欧美性猛交╳xxx乱大交人| 三级经典国产精品| 欧美 日韩 精品 国产| 久久久精品欧美日韩精品| 日韩欧美 国产精品| 观看免费一级毛片| 一区二区三区免费毛片| 日韩不卡一区二区三区视频在线| 伦精品一区二区三区| 久久久久国产精品人妻一区二区| 日本猛色少妇xxxxx猛交久久| 成人毛片60女人毛片免费| 成人美女网站在线观看视频| 国模一区二区三区四区视频| 免费av毛片视频| av.在线天堂| 高清在线视频一区二区三区| 六月丁香七月| 免费电影在线观看免费观看| 99视频精品全部免费 在线| 国产高清三级在线| 99久国产av精品国产电影| 中文精品一卡2卡3卡4更新| 人体艺术视频欧美日本| 中文字幕人妻熟人妻熟丝袜美| 国产免费视频播放在线视频| 日韩精品有码人妻一区| 国产黄片美女视频| 老师上课跳d突然被开到最大视频| videossex国产| 精品久久久久久久久av| 亚洲av国产av综合av卡| 欧美人与善性xxx| 亚洲精品国产av蜜桃| 丰满少妇做爰视频| 97超视频在线观看视频| 久久久久久九九精品二区国产| 亚洲精华国产精华液的使用体验| 热re99久久精品国产66热6| 日本av手机在线免费观看| 日韩欧美精品v在线| 久久久久久久久久久丰满| 午夜免费鲁丝| av天堂中文字幕网| 亚洲在线观看片| 国产亚洲一区二区精品| 亚洲高清免费不卡视频| 少妇人妻久久综合中文| 我的老师免费观看完整版| 亚洲精品一二三| 亚洲精品一区蜜桃| 亚洲成人久久爱视频| 少妇裸体淫交视频免费看高清| 午夜福利视频1000在线观看| 狂野欧美白嫩少妇大欣赏| 日本黄色片子视频| 国产91av在线免费观看| 直男gayav资源| 免费人成在线观看视频色| 国产精品99久久99久久久不卡 | 欧美一级a爱片免费观看看| 亚洲四区av| 国产成人免费观看mmmm| 亚洲精品视频女| 精品视频人人做人人爽| 伊人久久国产一区二区| 十八禁网站网址无遮挡 | 一级毛片电影观看| 国产免费视频播放在线视频| 亚洲国产精品国产精品| 人妻制服诱惑在线中文字幕| 日韩一本色道免费dvd| 噜噜噜噜噜久久久久久91| 18禁裸乳无遮挡免费网站照片| 99热这里只有是精品在线观看| 国产淫语在线视频| 少妇高潮的动态图| 蜜臀久久99精品久久宅男| 免费电影在线观看免费观看| 激情五月婷婷亚洲| 精品一区二区三区视频在线| 日本一本二区三区精品| 国产免费福利视频在线观看| 国产 一区精品| 国产精品精品国产色婷婷| 国产精品嫩草影院av在线观看| 日日啪夜夜爽| 国产一区二区三区综合在线观看 | 午夜福利视频精品| 99热这里只有是精品在线观看| 成人欧美大片| 麻豆久久精品国产亚洲av| 亚洲成人中文字幕在线播放| 观看美女的网站| 国产在线一区二区三区精| 亚洲性久久影院| av在线观看视频网站免费| 一个人观看的视频www高清免费观看| a级毛色黄片| 欧美97在线视频| 国产爱豆传媒在线观看| 亚洲欧美日韩无卡精品| 亚洲国产欧美人成| 久久99热6这里只有精品| tube8黄色片| av天堂中文字幕网| 免费av毛片视频| 天天躁日日操中文字幕| 久久久色成人| 国产欧美日韩一区二区三区在线 | 噜噜噜噜噜久久久久久91| 天堂俺去俺来也www色官网| 干丝袜人妻中文字幕| 国产毛片在线视频| av.在线天堂| 自拍偷自拍亚洲精品老妇| 免费看a级黄色片| 国产一区二区亚洲精品在线观看| 女人十人毛片免费观看3o分钟| 韩国高清视频一区二区三区| 午夜免费观看性视频| 黄色日韩在线| 街头女战士在线观看网站| 欧美国产精品一级二级三级 | 好男人视频免费观看在线| 九九爱精品视频在线观看| 少妇的逼水好多| 你懂的网址亚洲精品在线观看| 国产永久视频网站| 高清午夜精品一区二区三区| 中文精品一卡2卡3卡4更新| 亚洲精品456在线播放app| 三级经典国产精品| 搡女人真爽免费视频火全软件| 嫩草影院新地址| 亚洲精品中文字幕在线视频 | 亚洲一级一片aⅴ在线观看| 在线观看美女被高潮喷水网站| 色视频www国产| 熟女av电影| 亚洲精品视频女| 久久99热这里只频精品6学生| 欧美成人精品欧美一级黄| 亚洲丝袜综合中文字幕| 22中文网久久字幕| 亚洲国产欧美人成| 青春草国产在线视频| 国产探花极品一区二区| 国产精品偷伦视频观看了| 欧美日韩视频高清一区二区三区二| 国产精品蜜桃在线观看| 成年女人看的毛片在线观看| 嫩草影院新地址| 欧美老熟妇乱子伦牲交| 久久国产乱子免费精品| 99热6这里只有精品| 午夜福利高清视频| 精品一区二区三区视频在线| 女人久久www免费人成看片| 日本av手机在线免费观看| 高清毛片免费看| 我的老师免费观看完整版| 日本色播在线视频| 18禁裸乳无遮挡动漫免费视频 | 男人狂女人下面高潮的视频| 欧美三级亚洲精品| 精品久久国产蜜桃| 久久久久久久国产电影| 99久久九九国产精品国产免费| 日本午夜av视频| 欧美日韩国产mv在线观看视频 | 三级国产精品片| 99久久人妻综合| h日本视频在线播放| 国产乱人视频| 国产乱人视频| 激情五月婷婷亚洲| 各种免费的搞黄视频| 国产成人免费观看mmmm| 2018国产大陆天天弄谢| 丰满乱子伦码专区| av卡一久久| 久热久热在线精品观看| 亚洲欧美成人精品一区二区| 国产熟女欧美一区二区| 成人免费观看视频高清| 久久99蜜桃精品久久| 国产亚洲最大av| 国产色爽女视频免费观看| 22中文网久久字幕| 天天躁日日操中文字幕| 美女国产视频在线观看| 爱豆传媒免费全集在线观看| 日本爱情动作片www.在线观看| 91在线精品国自产拍蜜月| 春色校园在线视频观看| 高清在线视频一区二区三区| 麻豆成人av视频| 在线看a的网站| 国产精品秋霞免费鲁丝片| av在线播放精品| 欧美xxxx黑人xx丫x性爽| 又爽又黄无遮挡网站| 男女无遮挡免费网站观看| 国产女主播在线喷水免费视频网站| 中文欧美无线码| 男女国产视频网站| 青青草视频在线视频观看| 麻豆乱淫一区二区| 91午夜精品亚洲一区二区三区| 天天躁日日操中文字幕| 一个人看视频在线观看www免费| 国产极品天堂在线| 日本-黄色视频高清免费观看| 久久国产乱子免费精品| 99久久精品国产国产毛片| a级一级毛片免费在线观看| 青春草视频在线免费观看| 国产精品无大码| 亚洲熟女精品中文字幕| 青青草视频在线视频观看| av在线播放精品| 18禁动态无遮挡网站| 中国美白少妇内射xxxbb| 亚洲精品日韩在线中文字幕| 91久久精品国产一区二区成人| 亚洲精品成人av观看孕妇| 国产伦理片在线播放av一区| 日韩制服骚丝袜av| av国产精品久久久久影院| 国产在视频线精品| 在线观看一区二区三区激情| 午夜日本视频在线| 国产亚洲最大av| 国产乱来视频区| av播播在线观看一区| 亚洲婷婷狠狠爱综合网| 国产av码专区亚洲av| 日本av手机在线免费观看| 国产毛片a区久久久久| 亚洲精品国产色婷婷电影| kizo精华| 丰满乱子伦码专区| 尾随美女入室| 九色成人免费人妻av| 国产色婷婷99| 日韩制服骚丝袜av| 青春草视频在线免费观看| 日本wwww免费看| 97在线人人人人妻| 免费观看av网站的网址| 国产 精品1| 免费av不卡在线播放| 亚洲va在线va天堂va国产| 午夜视频国产福利| 综合色av麻豆| 99久久人妻综合| 黑人高潮一二区| 日本一本二区三区精品| 国产精品国产三级国产av玫瑰| 中文在线观看免费www的网站| 2021天堂中文幕一二区在线观| 成人二区视频| 91aial.com中文字幕在线观看| 嘟嘟电影网在线观看| 亚洲欧美一区二区三区黑人 | 精品熟女少妇av免费看| 偷拍熟女少妇极品色| 不卡视频在线观看欧美| 99九九线精品视频在线观看视频| 国产精品99久久久久久久久| 国产精品女同一区二区软件| 五月开心婷婷网| 国产高清国产精品国产三级 | 国产免费一级a男人的天堂| 精品人妻熟女av久视频| 欧美日韩一区二区视频在线观看视频在线 | 高清av免费在线| 黄色视频在线播放观看不卡| 免费看光身美女| 日韩大片免费观看网站| 亚洲欧美日韩另类电影网站 | 亚洲精品,欧美精品| 熟女人妻精品中文字幕| 亚洲熟女精品中文字幕| 2021少妇久久久久久久久久久| 边亲边吃奶的免费视频| 禁无遮挡网站| 边亲边吃奶的免费视频| 精品一区二区三区视频在线| 中文欧美无线码| 女的被弄到高潮叫床怎么办| 精品国产乱码久久久久久小说| 国产亚洲av片在线观看秒播厂| 日本av手机在线免费观看| 国产成人a区在线观看| 久久精品久久久久久噜噜老黄| 啦啦啦啦在线视频资源| 国产黄色视频一区二区在线观看| 欧美精品一区二区大全| 啦啦啦啦在线视频资源| 国产精品无大码| 高清毛片免费看| 亚洲综合色惰| 免费黄频网站在线观看国产| 大码成人一级视频| 又爽又黄a免费视频| 国产成人福利小说| 久热这里只有精品99| 99九九线精品视频在线观看视频| 免费大片黄手机在线观看| 69人妻影院| 国产有黄有色有爽视频| 性插视频无遮挡在线免费观看| 国产男女超爽视频在线观看| 国产免费视频播放在线视频| 天天一区二区日本电影三级| 一级a做视频免费观看| 又黄又爽又刺激的免费视频.| 午夜激情久久久久久久| 内地一区二区视频在线| 日本猛色少妇xxxxx猛交久久| 久热这里只有精品99| 久久人人爽人人爽人人片va| 内射极品少妇av片p| 搡老乐熟女国产| 网址你懂的国产日韩在线| 天美传媒精品一区二区| 色网站视频免费| 国产精品一区www在线观看| 日本免费在线观看一区| 国产精品熟女久久久久浪| 国产免费一级a男人的天堂| 午夜激情福利司机影院| 人妻夜夜爽99麻豆av| 色婷婷久久久亚洲欧美| 天堂中文最新版在线下载 | 国产伦在线观看视频一区| 欧美zozozo另类| 午夜免费观看性视频| 黄色怎么调成土黄色| 色网站视频免费| 久久久久久久精品精品| 精品少妇黑人巨大在线播放| 久久久久久久久久人人人人人人| 午夜福利在线观看免费完整高清在| 一区二区av电影网| 亚洲国产成人一精品久久久| 日本黄大片高清| 狂野欧美激情性xxxx在线观看| 新久久久久国产一级毛片| 日产精品乱码卡一卡2卡三| 午夜激情福利司机影院| 三级经典国产精品| 日韩精品有码人妻一区| 麻豆乱淫一区二区| 成年版毛片免费区| 又粗又硬又长又爽又黄的视频| 女人被狂操c到高潮| 久久99蜜桃精品久久| 老司机影院成人| 一级毛片黄色毛片免费观看视频| 91久久精品国产一区二区三区| 久久久久久久大尺度免费视频| 亚洲精品视频女| 色播亚洲综合网| 女人久久www免费人成看片| 又粗又硬又长又爽又黄的视频| 亚洲图色成人| 国产亚洲5aaaaa淫片| 婷婷色综合www| 亚洲经典国产精华液单| 丰满乱子伦码专区| 国产探花在线观看一区二区| 少妇丰满av| 下体分泌物呈黄色| 婷婷色麻豆天堂久久| 国产一级毛片在线| 国产精品人妻久久久影院| 免费看av在线观看网站| 在线观看一区二区三区| 天堂网av新在线| 国产男女内射视频| 十八禁网站网址无遮挡 | 制服丝袜香蕉在线| 欧美xxxx性猛交bbbb| eeuss影院久久| 国产精品国产三级国产专区5o| 中文欧美无线码| 新久久久久国产一级毛片| 国产精品人妻久久久影院| 国产毛片在线视频| 91在线精品国自产拍蜜月| 成人毛片a级毛片在线播放| 亚洲国产精品成人久久小说| 日韩伦理黄色片| av专区在线播放| 日韩强制内射视频| 美女脱内裤让男人舔精品视频| 大片免费播放器 马上看| 国内少妇人妻偷人精品xxx网站| 午夜福利视频1000在线观看| 两个人的视频大全免费| 丝袜脚勾引网站| a级一级毛片免费在线观看| 看十八女毛片水多多多| 超碰av人人做人人爽久久| 精品少妇黑人巨大在线播放| 欧美国产精品一级二级三级 | 亚洲成人精品中文字幕电影| 男人和女人高潮做爰伦理| 黄色一级大片看看| 各种免费的搞黄视频| 久久精品国产鲁丝片午夜精品| 欧美日韩视频精品一区| 一个人观看的视频www高清免费观看| 精品酒店卫生间| 超碰av人人做人人爽久久| 久久久a久久爽久久v久久| 三级经典国产精品| 日韩中字成人| 国产亚洲91精品色在线| 久久久久性生活片| 成人欧美大片| 日本wwww免费看| 日韩 亚洲 欧美在线| 22中文网久久字幕| 久久午夜福利片| 亚洲精品中文字幕在线视频 | 久热这里只有精品99| 身体一侧抽搐| 日本爱情动作片www.在线观看| 久久鲁丝午夜福利片| 成人特级av手机在线观看| 黄色欧美视频在线观看| 永久免费av网站大全| 精品久久久久久久人妻蜜臀av| 国产精品av视频在线免费观看| 如何舔出高潮| 免费观看av网站的网址| 国产极品天堂在线| 蜜桃亚洲精品一区二区三区| 国产成人精品久久久久久| 亚洲高清免费不卡视频| 人人妻人人看人人澡| 青春草视频在线免费观看| 成人国产av品久久久| 麻豆精品久久久久久蜜桃| 人体艺术视频欧美日本|