999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

高溫脅迫下金釵石斛花的生理及轉錄組響應分析

2025-02-09 00:00:00蒙紅霞李佳蔚龔建英王華新
熱帶作物學報 2025年2期
關鍵詞:植物差異分析

摘""要:金釵石斛具有較高的觀賞價值,是春石斛雜交品種群的重要親本之一,為探索金釵石斛花響應高溫脅迫的生理和分子機制,本研究結合水勢、細胞膜通透性監測、乙烯相關化合物含量測定,以及轉錄組學分析,以36"℃高溫脅迫5、10、24"h為處理組,常溫培養為對照組,研究36"℃下金釵石斛花的生理特性,分析其對高溫的響應轉錄組特征。結果表明:36"℃高溫脅迫在24"h內對金釵石斛花水分狀況和細胞膜通透性幾乎沒有影響,3個高溫處理組的金釵石斛花乙烯前體(ACC)含量與對照組(Dno-CK)相比均顯著增加;高溫處理5"h與常溫對照組(Dno-CK"vs"Dno-5h)、高溫處理10"h與常溫對照組(Dno-CK"vs"Dno-10h)、高溫處理24"h與常溫對照組(Dno-CK"vs"Dno-24h)3個比較組的差異基因分別為25"204、25"528、26"878個,上調表達與下調表達的基因量約為1∶1,隨著高溫脅迫時間的延長,金釵石斛花的差異表達的基因總數逐漸上升。KEGG富集分析結果表明,3個比較組差異基因均富集到單萜生物合成通路,說明高溫促進金釵石斛花揮發性萜類化合物的合成。GO功能富集分析結果表明,金釵石斛花受高溫脅迫大部分差異基因富集到遺傳信息處理途徑、代謝途徑相關的通路,在核苷酸切除修復、淀粉和蔗糖代謝、內質網中的蛋白質加工3條通路中差異基因數量最多,說明高溫激發金釵石斛花提高各類代謝活動響應熱脅迫。對差異表達基因中熱激蛋白基因(HSP)篩選發現,金釵石斛花的HSP基因大量上調表達響應高溫脅迫,金釵石斛花在24"h內響應高溫脅迫的HSP基因主要是HSP20基因,其次是HSP90基因。金釵石斛花的熱響應基因的進一步研究可針對HSP20基因和HSP90基因。本研究結果可為石斛耐熱育種研究提供參考,為耐高溫春石斛品系的培育提供理論支持。

關鍵詞:金釵石斛;高溫脅迫;生理特性;轉錄組分析;熱激蛋白中圖分類號:S567.239""""""文獻標志碼:A

Physiological"and"Transcriptome"Response"Analysis"of"Dendrobium"nobile"Flowers"under"High"Temperature"Stress

MENG"Hongxia1,2,"LI"Jiawei1,2*,"GONG"Jianying3,"WANG"Huaxin3

1."College"of"Forestry,"Guangxi"University"/"Guangxi"Key"Laboratory"of"Forest"Ecology"and"Conservation,"Nanning,"Guangxi""530004,"China;"2."State"Key"Laboratory"for"Conservation"and"Utilization"of"Subtropical"Agro-bioresources,"College"of"Forestry,"Guangxi"University,"Nanning,"Guangxi"530004,"China;"3."Guangxi"Forestry"Research"Institute,"Nanning,"Guangxi"530002,"China

Abstract:"Dendrobium"nobile"has"a"high"ornamental"value"and"is"one"of"the"important"parental"plants"in"the"hybridization"group"of"Dendrobium."To"explore"the"physiological"and"molecular"mechanisms"of"D."nobile"flowers"in"response"to"high-temperature"stress,"this"study"combined"water"potential,"cell"membrane"permeability"monitoring,"ethylene-related"compound"content"determination,"and"transcriptome"analysis."Treatment"groups"were"subjected"to"36"℃"high-tempe rat ure"stress"for"5"h,"10"h"and"24"h,"with"a"normal"temperature"culture"as"the"control."The"physiological"characteristics"of""D."nobile"flowers"under"36"℃"and"the"transcriptome"response"to"high"temperature"were"studied."The"results"shoed"that"36"℃"high-temperature"stress"had"almost"no"effect"on"the"water"status"and"cell"membrane"permeability"of"D."nobile"flowers"within"24h."The"content"of"ethylene"precursor"(ACC)"in"the"flowers"significantly"increased"in"all"three"high-temperature"treatment"groups"compared"to"the"control"(Dno-CK)."The"number"of"differentially"expressed"genes"between"the"5"h"high-temperature"treatment"and"the"control"group"(Dno-CK"vs"Dno-5h),"the"10"h"high-temperature"treatment"and"the"control"(Dno-CK"vs"Dno-10h),"and"the"24"h"high-temperature"treatment"and"the"control"(Dno-CK"vs"Dno-24h)"was"25"204,"25"528"and"26"878,"respectively,"with"the"number"of"genes"up-and"downregulated"was"approximately"1∶1."As"the"duration"of"high-temperature"stress"increased,"the"total"number"of"differentially"expressed"genes"in"D."nobile"flowers"gradually"rose."KEGG"enrichment"analysis"showed"that"the"differential"genes"of"the"three"comparison"groups"were"enriched"tonbsp;the"monoterpene"biosynthesis"pathway,"indicating"that"high"temperature"promoted"the"synthesis"of"volatile"terpenoids"in"D."nobile."Functional"enrichment"analysis"indicated"that"most"of"the"differentially"ex-pressed"genes"in"response"to"high-temperature"stress"were"enriched"in"pathways"related"to"genetic"information"processing"and"metabolic"pathways,"with"the"highest"number"of"differential"genes"in"the"nucleotide"excision"repair,"starch"and"sucrose"metabolism,"and"protein"processing"in"the"endoplasmic"reticulum"pathways."This"suggests"that"high"temperature"stimulates"D."nobile"flowers"to"enhance"various"metabolic"activities"in"response"to"heat"stress."Screening"of"differentially"expressed"genes"revealed"that"heat"shock"protein"(HSP)"genes"in"D."nobile"flowers"were"significantly"up-regulated"in"response"to"high-temperature"stress."The"HSP"genes"responding"to"high-temperature"stress"within"24"h"in"D."nobile"flowers"were"primarily"the"HSP20"genes,"followed"by"the"HSP90"genes."Further"research"on"the"heat"re-sponse"genes"of"D."nobile"flowers"can"focus"on"the"HSP20"and"HSP90"genes."This"study"on"the"physiological"and"transcriptome"response"of"""""D."nobile"flowers"to"high-temperature"stress"would"provide"a"reference"for"heat-tolerant"breeding"research"and"theoretical"support"for"the"cultivation"of"heat-resistant"D."varieties.

Keywords:"Dendrobium"nobile;"high"temperature"stress;"physiological"property;"transcriptome"analysis;"heat"shock"protein

DOI:"10.3969/j.issn.1000-2561.2025.02.003

蘭科石斛屬(Dendrobium)園藝種是盆花與切花中的名貴類群,可分為春石斛和秋石斛兩類雜交品種群。春石斛是以原產于東亞低緯度高海拔地區的落葉種類為親本,主要親本包括原產于中國的金釵石斛(D."nobile)。金釵石斛性喜溫暖濕潤,生于山地林中樹干上或山谷巖石上,目前在中國廣西、湖北、貴州、云南、西藏等地區均有分布[1]。花香可以增強觀賞植物的美學特性,在植物的繁殖過程中起著重要作用,金釵石斛花香輕清,具有較高觀賞價值的同時還可以用于鎮靜解郁、抗炎等[2]。金釵石斛生境的特殊性導致其生長發育容易受到各種環境因子的影響,高溫便是其中之一。目前針對植物高溫脅迫的研究主要集中在植物營養器官的形態、生理生化特性、光合作用與呼吸作用、分子調節機制等方面[3],對于花器官的直接關注較少。

植物水分生理狀態可以通過植物器官的水勢直接反映,植物在面臨脅迫時水勢變化較為敏感,大部分植物通過降低水勢應對脅迫,如景天植物在低鹽與高鹽脅迫下葉片水勢均呈下降趨勢[4],新疆大葉苜蓿在輕度與重度水分脅迫下葉片水勢均呈下降趨勢[5]。也有部分植物如豬毛菜采用升高葉水勢的方式適應脅迫環境[6]。高溫脅迫會導致細胞膜受損,內穩態失衡,電導率可以在一定程度上反映其受損程度。一般情況下,耐熱性強的植株細胞膜穩定性高于耐熱性差的植株,植物遭受溫度脅迫的程度與電導率大小呈正比關系[7-8]。氣態植物激素乙烯參與植物對脅迫環境的響應,如40"℃高溫脅迫下黑油冬白菜在48"h內乙烯含量呈直線上升狀態[9]。乙烯的合成主要由2種特定的酶促過程完成,ACC合酶(ACS)將底物轉化為ACC后,氧化酶(ACO)將ACC轉化為乙烯,研究表明,ACC(乙烯前體)預處理可以提高擬南芥在40"℃下的存活率,減少氧化損傷,并誘導其耐熱性[10-11]。

高溫脅迫啟動各種代謝途徑,目前已有大量與響應高溫脅迫高度相關的途徑被鑒別,如谷胱甘肽代謝途徑(glutathione"metabolism"pathways)、光合作用-天線蛋白代謝通路(photosyn thetic-an t enna"protein"metabolism"pathway)、葉綠素代謝通路(the"chlorophyll"metabolism"pathway)、碳代謝途徑(carbon"metabolic"pathways)、單萜類生物合成(monoterpene"class"biosynthesis)等[12-15]。處于高溫脅迫的植物通過各種熱響應基因的表達適應脅迫環境,近年研究發現,熱激蛋白(heat"shock"proteins,"HSP)在植物遭受高溫脅迫過程中大量表達,如火龍果在40"℃高溫脅迫下處理5"h后HSP21基因開始大量上調表達[16]。月季RcHSP18.1基因常溫下不表達,受到46"℃高溫脅迫2.5"h內大量表達[17]。三棱蝦脊蘭在高溫脅迫后的差異基因中鑒定到XTH、HSP等基因家族,且脅迫后HSP20、HSP40和HSP70基因全部上調表達[18]。在42"℃高溫脅迫處理6"h內,黃瓜HSP20、HSP70和HSP90基因也全部上調表達[19]。高溫脅迫會導致生物體中的蛋白質變性,HSP作為分子伴侶,協助蛋白質復性、穩定、細胞內轉運和降解,防止受損蛋白質的積累,維持細胞內環境的穩定性[20-21]。

前人針對石斛高溫脅迫分子機制的相關研究表明,翅萼石斛花在40"℃高溫脅迫3"h后有82個HSP基因明顯上調表達[22];鐵皮石斛經歷35"℃高溫脅迫后,與跨膜轉運功能相關的基因均下調表達,推測高溫脅迫會降低石斛的質膜轉運能力[23];40"℃高溫脅迫1"h后,鼓槌石斛花中表達量最高的10個基因均為HSP基因[24]。目前有關金釵石斛高溫脅迫轉錄組響應分析的相關研究較少,基于前人對石斛屬植物高溫脅迫的研究,本研究以金釵石斛作為試驗材料,常溫下的植株作為對照組,利用人工氣候箱進行36"℃高溫處理5、10、24"h,利用RNA-seq技術對處理后的金釵石斛花被進行轉錄組測序,利用液相色譜串聯質譜(LC-MS/MS)法測定花被乙烯類激素含量,同時測定水勢、電導率參數。通過解讀金釵石斛高溫脅迫生理響應及轉錄組響應,以期為金釵石斛的耐高溫品種選育提供科學依據。

1""材料與方法

1.1""材料

供試材料為金釵石斛(圖1),采自廣西壯族自治區林業科學研究院,于花期(2024年4月)選取長勢一致且健康的12株石斛進行試驗。高溫處理:在花朵完全展開時將植株放入人工氣候箱,設置36"℃高溫5"h(Dno-5h)、10"h(Dno-10h)、24"h(Dno-24h)3組處理,光照12"h/d。常溫培養(24~26"℃)的金釵石斛植株作為對照(Dno-CK),每個處理3株金釵石斛。各處理結束后采集花朵的花瓣放入50"mL離心管中,液氮速凍6~8"min,于?80"℃冰箱中保存,用于后續轉錄組測序及乙烯類激素含量測定。

1.2""方法

1.2.1""花朵水勢及電導率測定""對照組(Dno-"CK)采集常溫下的花朵,高溫處理組(Dno-5h、Dno-10h、Dno-24h)在處理結束后將植株拿出人工氣候箱后,迅速采集測定。花朵水勢采用壓力勢水勢儀(PMS,"Corvallis,Oregon,USA)測定,花朵電導率采用電導率儀(DDSJ-319L型)測定,每個處理在3株金釵石斛上各自采集1朵花朵測定。

1.2.2""乙烯類激素含量測定""采用液相色譜串聯質譜法(LC-MS/MS)測定,研磨Dno-CK、Dno-5h、Dno-10h、Dno-24h四個處理的樣品后加入內標混合溶液、提取劑,離心后80%甲醇溶液復溶,過濾,置于樣瓶里。設置液相條件和質譜條件提取代謝物采集質譜數據,采用Analyst"1.6.3軟件和MultiQuant"3.0.3軟件處理質譜數據,繪制激素的標準曲線并計算其含量。

1.2.3""轉錄組測序""(1)RNA提取、cDNA文庫構建及測序。采用乙醇沉淀法和CTAB-"PBIOZOL提取試驗樣本總RNA,使用Qubit熒光定量儀和Qsep400高通量生物片段分析儀對總RNA進行鑒定和定量。cDNA文庫構建步驟主要包括:富集mRNA,以小片段mRNA為模板反轉錄生成第一鏈cDNA并合成第二鏈cDNA,對雙鏈cDNA進行測序接頭連接后獲取250~350"bp的插入片段的文庫,經PCR擴增后構建完成文庫。庫檢合格后進行Illumina測序。

(2)轉錄組數據處理。使用fastp對原始數據進行過濾得到clean"reads,轉錄本組裝使用Trinity軟件完成,聚類去冗余由Corset軟件完成。利用DIAMOND軟件在KEGG、NR、Swiss-Prot、GO、COG/KOG、Trembl數據庫中對去冗余后的轉錄本序列進行注釋分析。

(3)差異表達基因分析。使用RSEM軟件計算轉錄本的FPKM值,以此定量基因表達水平。使用DESeq"2分析Dno-CK"vs"Dno-5h、Dno-CK"vs"Dno-10h、Dno-CK"vs"Dno-24h三組之間的差異表達基因(DEGs),差異基因的篩選條件為|log2"(Fold"Change)|≥1,且FDRlt;0.05。基于超幾何檢驗進行富集分析,KEGG富集分析以Pathway為單位進行超幾何分布檢驗,GO富集分析則基于GO-term進行。

1.2.4""揮發性成分分析""采用HS-SPME-GC-MS進行揮發性成分分析。選擇3株生長良好的金釵石斛于盛花期采集花朵,每株采集2朵花,保存于?80"℃冰箱中。樣品進行液氮研磨,每個樣本稱取約500"mg于頂空瓶中,加入飽和NaCl溶液和20"μL內標溶液,設置HS-SPME萃取條件、色譜條件和質譜條件,進行樣本萃取及GC-MS分析。通過MassHunter軟件處理質譜分析后的原始數據,對樣本的代謝物進行質譜定性定量分析,根據離子流峰面積歸一化法計算各組分在總揮發物中的相對含量。

1.3""數據處理

使用Excel"2016軟件處理數據,使用SPSS"26.0軟件評估差異的顯著性,使用Origin"2021軟

件制圖。

2""結果與分析

2.1""高溫脅迫下金釵石斛花朵生理響應

在電導率方面,3個高溫組處理的花朵電導率與對照組基本一致,保持在4.29~4.95"μs/cm,無顯著性差異(圖2A),說明36"℃高溫脅迫在24"h內對金釵石斛花朵細胞膜幾乎無影響。在水勢方面,3個高溫組處理的花水勢基本一致,保持在?1.36~?1.34"MPa,與對照組相比略有降低但無顯著性差異(圖2B),說明36"℃高溫脅迫在24"h內對金釵石斛花朵水勢的影響不大。乙烯前體(ACC)含量,與Dno-CK相比,Dno-5h增加了76.9%,差異顯著;Dno-10h增加了17.9%,差異顯著;Dno-24h增加了92.9%,差異顯著。金釵石斛在36"℃高溫脅迫24"h時ACC變化量達到最高,增長量最低點則在脅迫10"h時(圖2C)。

2.2""測序數據質量評估

本研究對4個處理結束后的金釵石斛花朵進行轉錄組測序,得到12組金釵石斛花朵的原始數據和質控數據(表1)。其中,過濾后獲得的clean"reads量共102.17"Gb,占原始數據(raw"reads)的98.01%以上。12組金釵石斛花朵的GC含量均大于45%,Q30堿基百分比保持在95%以上,Q20堿基百分比則保持在98%以上。表明本試驗的轉錄組測序質量較好,能滿足基因功能注釋和分析要求。

2.3""差異表達基因篩選和分析

3個高溫處理組與對照組Dno-CK差異表達基因的篩選結果顯示(圖3),金釵石斛經過5"h的36"℃高溫脅迫后共差異表達25"204個基因,包括11"703個上調表達基因和13"501個下調表達基因;經過10"h的36"℃高溫脅迫后共差異表達25"528個基因,表達上調和下調的差異基因分別為12"532個和12"996個;在24"h的36"℃高溫脅迫后,金釵石斛差異表達的基因總數有26"878個,包括13"876個上調表達基因和13"002個下調表達基因。隨著高溫脅迫時間的延長,金釵石斛差異表達的基因總數逐漸上升。與5"h相比,10"h高溫脅迫多了324個差異表達基因,而經過24"h高溫脅迫的金釵石斛比Dno-5h、Dno-10h分別多了1674個、1350個差異表達基因。說明金釵石斛遭受高溫脅迫5"h后差異基因數量有較大程度上升,且在熱處理24"h的差異基因數量保持持續增長的狀態。

2.4""差異表達基因KEGG分析

對高溫處理組與對照組的差異表達基因進行KEGG信號通路富集分析(圖4),結果表明,Dno-"CK"vs"Dno-5h組富集最顯著的20條通路中差異基因數量排名前5的通路主要包括剪接體(ko03040)、內質網中的蛋白質加工(ko04141)、RNA聚合酶(ko03020)、核苷酸切除修復(ko03420)、淀粉和蔗糖代謝(ko00500);Dno-CK"vs"Dno-10h組富集的通路主要包括內質網中的蛋白質加工(ko04141)、淀粉和蔗糖代謝(ko00500)、植物激素信號轉導(ko04075)、剪接體(ko03040)、次生代謝物的生物合成(ko01110);Dno-CK"vs"Dno-24h組富集的通路主要包括次生代謝物的生物合成(ko01110)、淀粉和蔗糖代謝(ko00500)、植物激素信號轉導(ko04075)、內質網中的蛋白質加工(ko04141)、核苷酸切除修復(ko03420)。

此外,Dno-CK"vs"Dno-5h、Dno-CK"vs"Dno-"10h、Dno-CK"vs"Dno-24h三組中差異表達基因顯著富集的相同通路有8個,分別是歸屬于遺傳信息處理過程的內質網中的蛋白質加工(ko04141)、核苷酸切除修復(ko03420)通路;歸屬于代謝過程的角質、木栓質和蠟生物合成(ko00073)、牛磺酸和低牛磺酸的代謝(ko00430)、氰氨基酸代紅色方框代表3個比較組合差異表達基因顯著富集的相同通路。

謝(ko00460)、淀粉和蔗糖代謝(ko00500)、脂肪酸降解(ko00071)、單萜生物合成(ko00902)通路。說明金釵石斛在經高溫脅迫5~24"h的過程中,遺傳信息處理途徑和代謝途徑發生了較大變化,其中內質網中的蛋白質加工(ko04141)、核苷酸切除修復(ko03420)、淀粉和蔗糖代謝(ko00500)3條通路富集程度較顯著,且差異基因數量較多,說明這3條信號通路可能在金釵石斛響應高溫脅迫上具有重要作用。

2.5""差異表達基因GO富集分析

對高溫處理組與對照組的差異表達基因進行GO富集分析(圖5)發現,Dno-CK"vs"Dno-5h組的差異表達基因GO富集有14"696個基因,Dno-CK"vs"Dno-10h組有14"865個基因,Dno-CK"vs"Dno-24h組有15"654個基因。分別取3組富集最顯著的9個GO"term進行分析,結果表明,Dno-CK"vs"Dno-5h組的差異表達基因主要富集于2-羥戊二酸依賴性雙加氧酶活性、伴侶蛋白結合、綠葉揮發性物質生物合成過程、sm樣蛋白家族復合物、二十碳烷酸代謝過程、解剖結構形態發生的調控、剪接體snRNP復合物、細胞對熱的響應和生長素穩態。Dno-CK"vs"Dno-10h組的差異表達基因主要富集于果膠分解代謝過程、酶抑制劑活性、伴侶蛋白結合、綠葉揮發性物質生物合成過程、二十碳烷酸代謝過程、天冬氨酸酯酶活性、L-甲硫氨酸從甲基硫代腺苷的回收、氨基酸回收和L-甲硫氨酸回收。Dno-CK"vs"Dno-24h組的差異表達基因主要富集于氧化還原酶活性、綠葉揮發性物質生物合成過程、黃酮類生物合成過程、外肽酶活性、絲氨酸型外肽酶活性、絲氨酸型羧基肽酶活性、羧肽酶活性、二十碳烷酸代謝過程和脂氧酶途徑。此外,Dno-CK"vs"Dno-5h、Dno-CK"vs"Dno-10h、Dno-CK"vs"Dno-24h三組中差異表達基因顯著富集的相同通路為綠葉揮發性物質生物合成過程和二十碳烷酸代謝過程。

2.6""熱激蛋白基因分析

熱激蛋白(heat"shock"proteins,"HSP)是在高溫下合成的重要蛋白質,已被證明參與對環境壓力的各種反應并調節許多發育過程。在高溫處理組與對照組的差異表達基因功能注釋列表中篩選紅色方框代表3個比較組合差異表達基因顯著富集的相同通路。

熱激蛋白表達相關基因,結果表明(附表1),Dno-CK"vs"Dno-5h組差異表達基因中有71個上調表達的熱激蛋白基因(48個HSP20基因,23個HSP90基因)和8個下調表達的熱激蛋白基因(4個HSP20基因,4個HSP90基因);Dno-CK"vs"Dno-10h組差異表達基因中有78個上調表達的熱激蛋白基因(49個HSP20基因,1個HSP70基因,28個HSP90基因)和9個下調表達的熱激蛋白基因(7個HSP20基因,2個HSP90基因);Dno-CK"vs"Dno-24h組差異表達基因中有69個上調表達的熱激蛋白基因(41個HSP20基因,28個HSP90基因)和6個下調表達的熱激蛋白基因(5個HSP20基因,1個HSP90基因)。說明金釵石斛熱激蛋白基因大量上調表達響應高溫脅迫且表達的熱激蛋白基因以HSP20為主。

對Dno-CK"vs"Dno-5h、Dno-CK"vs"Dno-10h、Dno-CK"vs"Dno-24h三組中共同上調表達的49個熱激蛋白基因進行分析,結果如圖6所示。有28個基因在5~24"h的高溫脅迫處理下表達量呈下降趨勢(圖6A);有12個基因在5~10"h的表達量呈上升趨勢,10~24"h的表達量呈下降趨勢,10"h的表達量最高(圖6B),而余下9個基因趨勢與其相反,在10"h的表達量最低(圖6C)。這49個熱激蛋白基因由33個HSP20基因和16個HSP90基因組成。此外,Cluster-61924.7(HSP90基因)的表達量在3個高溫脅迫處理中明顯高于其余48個基因,在Dno-CK"vs"Dno-5h組中的表達量為其余基因的1.17~8.63倍,在Dno-CK"vs"Dno-10h組中的表達量為1.82~15.28倍,在Dno-CK"vs"Dno-24h組中的表達量為1.81~14.19倍。推測金釵石斛的Cluster-61924.7在響應高溫脅迫過程中發揮重要作用。

2.7""金釵石斛花朵揮發性成分分析

對金釵石斛盛花期的花朵揮發性成分進行鑒定,共檢測出1480種揮發性物質,其中,萜類化合物的種類最多(296種),占總量的20%,其次為酯類(266種),占總量的18%,其余主要成分包括酮類、醇類、雜環化合物、醛類等(圖7A)。揮發性成分相對含量較高的化合物為雜環化合物、脂類、萜類、酮類,分別達到271.74、271.67、267.39、266.06"μg/g(圖7B)。對金釵石斛花中的揮發性萜類化合物進行分析,共有21種成分對花香具有貢獻(圖7C、圖7D),其中2,6-二甲基-7-辛烯-2-醇和3,7-二甲基-1,3,7-辛三烯的貢獻最大,"其相對含量分別為0.91、0.45"μg/g。

3""討論

短期高溫脅迫會對植物基因表達、生理活動、植株形態等造成一定的影響,為了生存和適應高溫的逆境,植物在進化過程中發展出了復雜的自我防御機制。一般情況下,植物組織通過升高或降低水勢響應脅迫環境,其水勢與吸水能力成反比關系,金釵石斛在面臨干旱脅迫時,葉片通過升高水勢延遲脫水來應對干旱環境[25]。本研究發現高溫處理下的花朵水勢與常溫對照組相比略有降低但無顯著性差異,說明36"℃高溫脅迫在24"h內對金釵石斛花朵水分狀況影響不高。同一高溫脅迫下不同種金線蓮葉片電導率趨勢不同,福建金線蓮于脅迫的第2天電導率顯著下降,第6天顯著上升,而臺灣金線蓮葉片在高溫脅迫6"d的電導率與對照組相比無顯著性差異,這與大部分研究高溫脅迫與電導率關系的結果不完全一致[26]。本研究中4個處理下的金釵石斛花朵電導率均無顯著性差異,與臺灣金線蓮葉片電導率的變化情況相似,說明金釵石斛對36"℃高溫脅迫具有較強的適應能力,此溫度處理對其細胞膜通透性幾乎沒有影響。趙暉等[27]研究表明在高溫脅迫下杜鵑花葉片乙烯釋放速率增高且乙烯合成相關酶活性升高,本研究中3個高溫處理組的金釵石斛花乙烯前體(ACC)含量與對照組Dno-CK相比均顯著增加,這與趙暉等的研究結果相似。

本研究利用轉錄組測序篩選出Dno-CK"vs"Dno-5h、Dno-CK"vs"Dno-10h、Dno-CK"vs"Dno-24h三個比較組的差異基因分別為25"204、25"528、26"878個,上調表達與下調表達的基因量約為1∶1,這些基因可以作為后續研究與金釵石斛響應高溫脅迫密切相關基因的候選基因。隨著高溫脅迫時間的延長,金釵石斛差異表達的基因總數逐漸上升。這與程秋如等[28]低溫脅迫處理油棕的研究結果一致,說明植物通過激活大量基因表達響應溫度脅迫。本研究前期對常溫狀態下金釵石斛進行揮發性成分分析發現,萜類化合物的種類及相對含量均較高,賦花香的主要成分為2,6-二甲基-7-辛烯-2-醇、3,7-二甲基-1,3,7-辛三烯。對差異基因進行KEGG功能富集分析發現,Dno-CK"vs"Dno-5h、Dno-CK"vs"Dno-10h、Dno-CK"vs"Dno-24h三個比較組差異基因均富集到單萜生物合成通路,說明高溫促進金釵石斛花揮發性萜類化合物的合成。對差異基因進行GO功能富集分析發現,金釵石斛受高溫脅迫大部分差異基因富集到遺傳信息處理途徑、代謝途徑相關的通路,在核苷酸切除修復、淀粉和蔗糖代謝、內質網中的蛋白質加工3條通路中差異基因數量最多。說明高溫激發金釵石斛提高各類代謝活動響應熱脅迫,這與菜心、茶樹等植物受高溫脅迫后差異基因富集結果一致[29-30]。

熱激蛋白是植物在溫度脅迫下產生的一類保護性蛋白質,高度保守并存在于幾乎所有生物體中,參與許多生長發育過程并對逆境脅迫做出反應[31]。對差異表達基因中HSP基因的篩選發現,金釵石斛熱激蛋白基因大量上調表達,響應高溫脅迫,Dno-CK"vs"Dno-5h、Dno-CK"vs"Dno-10h、Dno-CK"vs"Dno-24h三組中共同上調表達的49個熱激蛋白基因由33個HSP20基因和16個HSP90基因組成。這說明金釵石斛在24"h內響應高溫脅迫的HSP基因主要是HSP20基因,其次是HSP90基因。研究表明,HSP20是植物中最大的HSP家族,HSP20蛋白是許多高等植物中高溫脅迫誘導的最豐富的蛋白質,HSP20基因在植物的非生物脅迫中起著關鍵作用[32];HSP90基因與植物脅迫信號轉導、受體折疊、轉錄因子和激酶以及輔助細胞在脅迫下的存活有關[33-35]。共同上調表達的HSP基因大部分在5"h的表達量最高,5~24"h的表達量呈下降趨勢。這與朱白婢等[36]研究高溫處理后冬瓜耐熱基因在0~4"h的表達模式相同,熱響應基因在高溫脅迫短時間內迅速表達,植物適應高溫后開始逐漸下降。金釵石斛熱響應基因的進一步研究可針對HSP20基因和HSP90基因進行,本研究對金釵石斛高溫脅迫轉錄組響應的分析可為石斛耐熱育種研究提供參考。

參考文獻

[1]"明興加,"李博然,"趙紀峰,"張植瑋,"袁良琛,"王剛,"羅毅波."金釵、金釵石斛的名實考證[J]."中國中藥雜志,"2016,"41(10):"1956-1964.MING"X"J,"LI"B"R,"ZHAO"J"F,"ZHANG"Z"W,"YUAN"L"C,"WANG"G,"LUO"Y"B."Investigation"of"original"materials"under"traditonal"Chinese"medicine"names"of"“Jinchai”"and"“Jinchai"Shihu”[J]."China"Journal"of"Chinese"Materia"Medica,"2016,"41(10):"1956-1964."(in"Chinese)

[2]"顏沛沛,"葉煒,"江金蘭,"王培育."金釵石斛花香氣日變化規律[J]."亞熱帶農業研究,"2021,"17(3):"179-183.YAN"P"P,"YE"W,"JIANG"J"L,"WANG"P"Y."Diurnal"variation"of"floral"scents"emitted"from"Dendrobium"nobile"Lindl."[J]."Subtropical"Agriculture"Research,"2021,"17(3):"179-183."(in"Chinese)

[3]"江珊,"吳龍英,"趙寶生,"黃佳惠,"蔣宇喆,"焦元,"黃進."植物耐受高溫脅迫的分子機制研究進展[J]."中國農學通報,"2024,"40(9):"132-138.JIANG"S,"WU"Lnbsp;Y,"ZHAO"B"S,"HUANG"J"H,"JIANG"Y"J,"JIAO"Y,"HUANG"J."Molecular"mechanism"of"heat"stress"tolerance"in"plants:"a"review[J]."Chinese"Agricultural"Science"Bulletin,"2024,"40(9):"132-138."(in"Chinese)

[4]"張國軍."鹽脅迫對4種景天幼苗水勢、熒光效率、丙二醛的影響[J]."農業科技與裝備,"2012(7):"3-6.ZHANG"G"J."Effect"of"salt"stress"on"the"water"potential,"fluorescence"efficiency,"malonaldehyde"of"four"stonecrop"seedlings[J]."Agricultural"Science"amp;"Technology"and"Equipment,"2012(7):"3-6."(in"Chinese)

[5]"南思睿,"羅永忠,"于思敏,"仝慧鑫,"何鈺."水分脅迫對紫花苜蓿葉片水勢及生物量分配的影響[J]."草原與草坪,"2022,"42(6):"63-70.NAN"S"R,"LUO"Y"Z,"YU"S"M,"TONG"H"X,"HE"Y."Effects"of"water"stress"on"leaf"water"potential"and"biomass"allocation"of"Medicago"sativa[J]."Grassland"and"Turf,"2022,"42(6):"63-70."(in"Chinese)

[6]"付曉玥,"閆建成,"梁存柱,"賈成朕,"史斌,"張景慧,"齊小軍."干旱與半干旱區一年生植物水勢對模擬降水變化的響應[J]."內蒙古大學學報(自然科學版),"2012,"43(2):"160-167.FU"X"Y,"YAN"J"C,"LIANG"C"Z,"JIA"C"Z,"SHI"B,"ZHANG"J"H,"QI"X"J."Water"potentials"of"annual"plants"responseto"simulated"rainfall"in"aridand"semiarid"regions[J]."Journal"of"Inner"Mongolia"University"(Natural"Science"Edition),"2012,"43(2):"160-167."(in"Chinese)

[7]"胡娜,"孫玉紅,"張國新,"陳思."6個品種櫻花在北京地區的適應性研究[J]."亞熱帶植物科學,"2021,"50(1):"33-38.HU"N,"SUN"Y"H,"ZHANG"G"X,"CHEN"S."Adaptability"of"six"cultivars"of"cerasus"in"Beijing[J]."Subtropical"Plant"Science,"2021,"50(1):"33-38."(in"Chinese)

[8]"王宏輝,"顧俊杰,"房偉民,"陳發棣,"張棟梁."不同品種紅掌葉片細胞膜穩定性和滲透調節物質對高溫脅迫的響應[J]."上海農業學報,"2020,"36(6):"24-27.WANG"H"H,"GU"J"J,"FANG"W"M,"CHEN"F"D,"ZHANG"D"L."Response"of"cell"membrane"stability"and"osmotic"adjustment"substances"in"leaves"of"different"cultivars"of"Anthurium"andraeanum"under"high"temperature"stress[J]."Acta"Agriculturae"Shanghai,"2020,"36(6):"24-27."(in"Chinese)

[9]"高璐,"李曉鋒,"奚丹丹,"朱玉英,"朱紅芳."不同溫度條件下淹水脅迫對不結球白菜植物激素的影響[J]."上海農業學報,"2022,"38(1):"33-40.GAO"L,"LI"X"F,"XI"D"D,"ZHU"Y"Y,"ZHU"H"F."Effects"of"waterlogging"stress"on"plant"hormones"in"non-heading"Chinese"cabbage"under"different"temperature"conditions[J]."Acta"Agriculturae"Shanghai,"2022,"38(1):"33-40."(in"Chinese)

[10]"BOLLER"T,"HERNER"R"C,"KENDE"H."Assay"for"and"enzymatic"formation"of"an"ethylenenbsp;precursor,"1-aminocyclo pro p ane-1-carboxylic"acid[J]."Planta,"1979,"145:"293-303.

[11]"LARKINDALE"J,"KNIGHT"M"R."Protection"against"heat"stress-induced"oxidative"damage"in"Arabidopsis"involves"calcium,"abscisic"acid,"ethylene,"and"salicylic"acid[J]."Plant"Physiol,"2002,"128:"682-695.

[12]"朱丹,"秦亞東,"李林華,"高煥."高溫脅迫誘導薄荷差異表達基因的分析[J]."安徽科技學院學報,"2022,"36(6):"54-62.ZHU"D,"QIN"Y"D,"LI"L"H,"GAO"H."Analysis"of"differentially"expressed"genesin"Mentha"haplocalyx"Briq."inducedby"high"temperature"stress[J]."Journal"of"Anhui"Science"and"Technology"University,"2022,"36(6):"54-62."(in"Chinese)

[13]"周娜,"鄭陽,"陸景偉,"胡燕,"陶偉林,"雷開榮,"潘曉雪."高溫脅迫下蘿卜苗期的轉錄組分析[J]."南方農業學報,2022,"53(3):"665-675.ZHOU"N,"ZHENG"Y,"LU"J"W,"HU"Y,"TAO"W"L,"LEI"K"R,"PAN"X"X."Transcriptome"analysis"of"Raphanus"sativus"seedlings"in"response"to"high"temperature"stress[J]."Journal"of"Southern"Agriculture,"2022,"53(3):"665-675."(in"Chinese)

[14]"郝力慧,"董彬,"朱紹華,"馬進."牡丹響應高溫脅迫的轉錄組分析及PsHSP基因表達[J]."浙江農林大學學報,"2021,"38(4):"802-811.HAO"L"H,"DONG"B,"ZHU"S"H,"MA"J."Transcriptome"analysis"and"PsHSP"gene"expression"of"Paeonia"suffruticosa"in"response"to"high"temperature"stress[J]."Journal"of"Zhejiang"A"amp;"F"University,"2021,"38(4):"802-811."(in"Chinese)

[15]"宋姣敏,"徐小萍,"陳裕坤,"陳桂信,"賴鐘雄."熱脅迫下水仙花葉片的轉錄組分析[J]."熱帶作物學報,2019,"40(10):"1991-2000.SONG"J"M,"XU"X"P,"CHEN"Y"K,"CHEN"G"X,"LAI"Z"X."An a l yses"of"leaf"transcriptome"under"heat"stress"in"narcissus[J]."Chinese"Journal"of"Tropical"Crops,"2019,"40(10):"1991-2000."(in"Chinese)

[16]"李俊成,"黨蕓芝,"孫清明."高溫脅迫下火龍果轉錄組及熱激蛋白響應分析[J]."浙江農業學報,"2024,"36(5):"1067-"1075.LI"J"C,"DANG"Y"Z,"SUN"Q"M."Transcriptome"analysis"and"HSP"response"of"pitaya"(Selenicereus"spp.)"under"heat"stress[J]."Acta"Agriculturae"Zhejiangensis,"2024,"36(5):"1067-1075."(in"Chinese)

[17]"楊碧楠,"李博文,"楊振宇,"徐奕芃,"閻韻清,"婁玉霞,"奉樹成,"明鳳."月季‘安吉拉’熱脅迫反應機制及功能基因的挖掘[J]."園藝學報,"2024,"51(6):"1284-1296.YANG"B"N,"LI"B"W,"YANG"Z"N,"XU"Y"P,"YAN"Y"Q,"LOU"Y"X,"FENG"S"C,"MING"F."The"mechanism"of"heat"stress"response"and"the"exploration"of"functional"genes"in"Rosa"chinensis"'Angela'[J]."Acta"Horticulturae"Sinica,"2024,"51(6):"1284-1296."(in"Chinese)

[18]"陳倫妍,"朱雅婷,"陸艾鮮,"吳悠,"冀宏宇,"陳志清,"吳沙沙,"翟俊文."2種蝦脊蘭屬(Calanthe"R."Br.)植物熱脅迫相關基因表達模式分析[J]."植物遺傳資源學報,"2023,"24(5):"1413-1425.CHEN"L"Y,"ZHU"Y"T,"LU"A"X,"WU"Y,"JI"H"Y,"CHEN"Z"Q,"WU"S"S,"ZHAI"J"W."The"analysis"of"expression"patterns"of"heat"stress-related"genes"in"two"species"of"Calanthe"R."Br.[J]."Journal"of"Plant"Genetic"Resources,"2023,"24(5):"1413-1425."(in"Chinese)

[19]"陳雪倩,"韓靜,"任仲海."黃瓜響應高溫脅迫轉錄組分析[J]."分子植物育種,"2021,"19(12):"3905-3914.CHEN"X"Q,"HAN"J,"REN"Z"H."Transcriptome"analysis"of"cucumber"in"response"to"high"temperature"stress[J]."Molecular"Plant"Breeding,"2021,"19(12):"3905-3914."(in"Chinese)

[20]"SONG"N"H,"AHN"Y"J."DcHsp17."7,"a"small"heat"shock"protein"in"carrot,"is"tissue-specifically"expressed"under"salt"stress"and"confers"tolerance"to"salinity[J]."New"Biotechnol,"2011,"28:"698-704.

[21]"QIN"F,"YU"B,"LI"W."Heat"shock"protein"101"(HSP101)"promotes"flowering"under"nonstress"conditions[J]."Plant"Physiol,"2021,"186:"407-419".

[22]"陳宏宇,"于瑩."翅萼石斛花在溫度脅迫下的轉錄組分析[J]."安徽農學通報,"2023,"29(15):"37-42.CHEN"H"Y,"YU"Y."Transcriptome"analysis"of"Dendrobium"cariniferum"under"temperature"stress[J]."Anhui"Agricultural"Science"Bulletin,"2023,"29(15):"37-"42."(in"Chinese)

[23]"徐載萍,"王孟穎,"孫福英,"趙雯,"楊莘鍇,"萬林春,"曾黎明."高溫脅迫下鐵皮石斛的轉錄組分析[J]."中國中藥雜志,"2024,"49(20):"5433-5440.XU"Z"P,"WANG"M"Y,"SUN"F"Y,"ZHAO"W,"YANG"X"K,"WAN"L"C,"ZENG"L"M."Transcriptome"analysis"of"Dendr obi um"officinale"exposed"to"high"temperature"stress[J]."China"Journal"of"Chinese"Materia"Medica,"2024,"49(20):"5433-"5440."(in"Chinese)

[24]"陳宏宇,"于瑩,"李菲."基于轉錄組鑒定鼓槌石斛花溫度脅迫響應基因[J]."高師理科學刊,"2023,"43(11):"44-49.CHEN"H"Y,"YU"Y,"LI"F."Identifying"genes"responsive"to"temperature"stress"in"Dendrobium"chrysotoxum"flower"based"on"transcriptome"analysis[J]."Journal"of"Science"of"Teachers'"College"and"University,"2023,"43(11):"44-49."(in"Chinese)

[25]"呂朝燕,"高智席,"嚴羽,"王利霞,"魏英,"葉麗,"盧珊珊,"張琴."干旱-復水對兩種石斛屬植物葉水勢的影響[J]."廣西植物,"2021,"41(2):"177-182.LYU"C"Y,"GAO"Z"X,"YAN"Y,"WANG"L"X,"WEI"Y,"YE"L,"LU"S"S,"ZHANG"Q."Effects"of"drought-rewatering"on"leaf"water"potential"of"two"Dendrobium"plants[J]."Guihaia,"2021,"41(2):"177-182."(in"Chinese)

[26]"林曉紅,"施木田,"林三睦."高溫脅迫對金線蓮葉綠素熒光參數及SOD活性與電導率的影響[J]."熱帶作物學報,"2014,"35(6):"1137-1142.LIN"X"H,"SHI"M"T,"LIN"S"M."Effects"of"high-temperature"stress"on"chlorophyll"fluorescence"parameters,"SOD"activity"and"electrolyte"leakage"of"Anoectochilus"roxburghii"(Wall.)"Lindl."and"Anoectochilus"formosanus"Hayata[J]."Chinese"Jour n al"of"Tropical"Crops,"2014,"35(6):"1137-1142."(in"Chinese)

[27]"趙暉,"耿興敏,"王露露,"許世達."乙烯在杜鵑花耐熱機制中的作用研究[J]."園藝學報,"2022,"49(3):"561-570.ZHAO"H,"GENG"X"M,"WANG"L"L,"XU"S"D."Research"on"the"effect"of"ethylene"in"heat"resistance"mechanism"of"Rhododendron[J]."Acta"Horticulturae"Sinica,"2022,"49(3):"561-570."(in"Chinese)

[28]"程秋如,"劉加能,"劉子凡,"潘登浪,"曾憲海,"曾精."油棕響應低溫脅迫的轉錄組分析[J/OL]."分子植物育種,"(2024-"06-25)[2024-08-27]."http://kns.cnki.net/kcms/d etail/46.1068."S.20240624.1745.008.html.CHENG"Q"R,"LIU"J"N,"LIU"Z"F,"PAN"D"L,"ZENG"X"H,"ZENG"J."Transcriptome"analysis"of"oil"palm"in"response"to"low"temperature"stress[J]."Molecular"Plant"Breeding,"(2024-06-25)[2024-08-27]."http://kns.cnki.net/kcms/detail/46."1068.S.20240624.1745.008.html."(in"Chinese)

[29]"江定,"李光光,"袁凡崇,"雷世康,"張華,"戴修純,"鄭巖松."菜心響應高溫脅迫的轉錄組分析與基因挖掘[J]."南方農業學報,"2024,"55(3):"766-783.JIANG"D,"LI"G"G,"YUAN"F"C,"LEI"S"K,"ZHANG"H,"DAI"X"C,"ZHENG"Y"S."Transcriptome"analysis"and"gene"mining"of"flowering"Chinese"cabbage"in"response"to"high"temperature"stress[J]."Journal"of"Southern"Agriculture,"2024,"55(3):"766-"783."(in"Chinese)

[30]"翟秀明,"唐敏,"李解,"陳世春,"胡方潔,"侯渝嘉."基于RNA-Seq技術的茶樹響應高溫脅迫轉錄組差異性分析[J]."分子植物育種,"2020,"18(17):"5629-5637.ZHAI"X"M,"TANG"M,"LI"J,"CHEN"S"C,"HU"F"J,"HOU"Y"J."Difference"analysis"of"heat"stress-responsive"transcriptome"of"Camellia"sinensis"based"on"RNA-Seq"technology[J]."Molecular"Plant"Breeding,"2020,"18(17):"5629-5637."(in"Chinese)

[31]"FEDER"M"E."Organismal,"ecological,"and"evolutionary"aspects"of"heat-shock"proteins"and"the"stress"response:"established"conclusions"and"unresolved"issues[J]."American"Zoologist,"1999,"39:"857-864.nbsp;(in"Chinese)

[32]"OUYANG"Y,"CHEN"J,"XIE"W,"WANG"L,"ZHANG"Q."Comp rehensive"sequence"and"expression"profile"analysis"of"Hsp20"gene"family"in"rice[J]."Plant"Molecular"Biology,"2009,"70(3):"341-57.

[33]"HAN"Y,"CHEN"Y,"YIN"S,"ZHANG"M,"WANG"W."Over-ex pre ssion"of"TaEXPB23,"a"wheat"expansin"gene,"improves"oxidative"stress"tolerance"in"transgenic"tobacco"plants[J]."Plant"Physiol,"2015,"173:"62-71.

[34]"CHA"J,"AHN"G,"KIM"J"Y,"KANG"S"B,"KIM"M"R,"SU'UDI"M,"KIM"W"Y,"SON"D."Structural"and"functional"differences"of"cytosolic"90-kDa"heat-shock"proteins"(Hsp90s)"in"Arabidopsis"thaliana[J]."Plant"Physiology"and"Biochemistry,"2013,"70:"368-373.

[35]"JACKSON"S"E,"QUEITSCH"C,"TOFT"D."Hsp90:"from"structure"to"phenotype[J]."Nature"Structural"amp;"Molecular"Biology,"2004,"11(12):"1152-1155.

[36]"朱白婢,"王靜,"張文,"胡艷平,"李雪嶠."冬瓜響應高溫脅迫的轉錄組測序及分析[J/OL]."分子植物育種,"(2024-04-"30)[2024-08-27]."http://kns.cnki.net/kcms/detail/46.1068.S."20 240430.0955.004.html.ZHU"B"B,"WANG"J,"ZHANG"W,"HU"Y"P,"LI"X"Q."Tran sc ri pt ome"sequencing"and"analysis"of"wax"gourd"in"resp on se"to"high"temperature"stress[J]."Molecular"Plant"Breeding,"(2024-"04-30)[2024-08-27]."http://kns.cnki.net/kcms/detail/46.1068."S.20240430.0955.004.html."(in"Chinese)

猜你喜歡
植物差異分析
相似與差異
音樂探索(2022年2期)2022-05-30 21:01:37
隱蔽失效適航要求符合性驗證分析
找句子差異
電力系統不平衡分析
電子制作(2018年18期)2018-11-14 01:48:24
生物為什么會有差異?
哦,不怕,不怕
將植物穿身上
電力系統及其自動化發展趨勢分析
植物罷工啦?
植物也瘋狂
主站蜘蛛池模板: 国产一区二区丝袜高跟鞋| 亚洲天堂色色人体| 国产在线观看99| 这里只有精品在线| 免费可以看的无遮挡av无码| 免费在线视频a| 在线播放精品一区二区啪视频| 无码日韩视频| 一本大道东京热无码av| 欧美日本一区二区三区免费| 国产在线八区| 嫩草国产在线| 国产网站免费| 无码日韩人妻精品久久蜜桃| 欧美日韩中文字幕二区三区| 亚洲无码37.| 综合色区亚洲熟妇在线| 最新亚洲av女人的天堂| 精品国产成人av免费| 国产人成乱码视频免费观看| 国产成人1024精品| 国产91全国探花系列在线播放 | 99视频全部免费| 亚洲精品无码av中文字幕| 亚洲视频在线观看免费视频| 伊人无码视屏| 欧美成一级| 亚洲视频免| 国产在线观看一区精品| 热思思久久免费视频| 亚洲天堂成人在线观看| 成人在线亚洲| 国产精品免费p区| 99999久久久久久亚洲| 谁有在线观看日韩亚洲最新视频 | 精品亚洲麻豆1区2区3区 | 欧美日韩精品综合在线一区| 国产精品青青| 国产免费怡红院视频| 国产v欧美v日韩v综合精品| 波多野结衣爽到高潮漏水大喷| 精品国产欧美精品v| 亚洲日韩Av中文字幕无码| 永久免费无码日韩视频| 亚洲精品麻豆| 久久亚洲日本不卡一区二区| 国产精品尤物铁牛tv | 波多野结衣一区二区三视频 | 精品国产电影久久九九| 亚洲三级a| 日韩在线网址| 国产精品女人呻吟在线观看| 综合久久五月天| 98超碰在线观看| 国产99视频精品免费观看9e| 亚洲中文字幕手机在线第一页| 日韩在线永久免费播放| 色综合网址| 亚洲国语自产一区第二页| 欧美视频二区| 91人人妻人人做人人爽男同| 青草91视频免费观看| 性色在线视频精品| 91精品最新国内在线播放| 美女无遮挡被啪啪到高潮免费| 国产精品99久久久久久董美香| 国产精品永久在线| 欧美成人免费午夜全| 国产迷奸在线看| 国产成人亚洲日韩欧美电影| 国产一区二区三区在线无码| 国产高潮流白浆视频| 久久久精品国产亚洲AV日韩| 国产精品网址你懂的| 国产亚洲男人的天堂在线观看| 精品撒尿视频一区二区三区| 天天综合网色中文字幕| 欧美精品aⅴ在线视频| 久久久久免费看成人影片| 日本不卡视频在线| 亚洲色图综合在线| 成人日韩精品|