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

    草莓應(yīng)答炭疽菌侵染的轉(zhuǎn)錄組分析

    2020-06-08 10:55:23陳哲黃靜趙佳梁宏
    植物保護(hù) 2020年3期
    關(guān)鍵詞:差異表達(dá)基因草莓

    陳哲 黃靜 趙佳 梁宏

    摘要 本研究以栽培草莓為材料,用炭疽菌對植株進(jìn)行侵染,對健康草莓和患病草莓進(jìn)行轉(zhuǎn)錄組分析(RNA-seq),共得到86 988個(gè)序列以及2 127個(gè)差異表達(dá)基因。利用GO數(shù)據(jù)庫對差異表達(dá)基因進(jìn)行了聚類分析,結(jié)果發(fā)現(xiàn)差異表達(dá)基因主要聚類在光合作用、氧化還原過程、氧化還原酶活性、碳固定過程以及糖代謝過程等方面。同時(shí),利用KEGG數(shù)據(jù)庫對差異表達(dá)基因影響的通路進(jìn)行注釋后發(fā)現(xiàn),炭疽菌侵染主要影響了光合作用、倍半萜和三萜生物途徑、類黃酮生物合成、谷胱甘肽代謝和植物激素信號(hào)轉(zhuǎn)導(dǎo)等途徑中關(guān)鍵基因的表達(dá)水平。挑選了20個(gè)差異表達(dá)基因進(jìn)行qRT-PCR,結(jié)果有19個(gè)基因的表達(dá)趨勢與轉(zhuǎn)錄測序結(jié)果一致。本研究獲得的草莓應(yīng)答炭疽菌侵染轉(zhuǎn)錄組數(shù)據(jù)信息,將有助于豐富草莓抗炭疽病的調(diào)控機(jī)制。

    關(guān)鍵詞 草莓; 炭疽菌; 轉(zhuǎn)錄組分析; 差異表達(dá)基因

    中圖分類號(hào): S 436.63

    文獻(xiàn)標(biāo)識(shí)碼: A

    DOI: 10.16688/j.zwbh.2019096

    Transcriptomics analysis of strawberry response to

    Colletotrichum theobromicola infection

    CHEN Zhe, HUANG Jing, ZHAO Jia, LIANG Hong

    (Biotechnology Research Center, Shanxi Academy of Agriculture Sciences, Taiyuan 030031, China)

    Abstract

    In this study, we analyzed the response of strawberry plants(Fragaria×ananassa) to infection with Colletotrichum theobromicola using the next-generation sequencing (RNA-seq) and investigated the alterations in gene expression between the healthy and infected plants. We have identified 86 988 unigenes, of which 2 127 were differentially expressed in the fungi-infected plants (DEGs). The unigenes were annotated and classified with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. The DEGs annotated to the GO database were distributed in many functional types, including photosynthesis, oxidation-reduction process, oxidoreductase activity and growth, development and metabolism. The KEGG pathway analysis of these genes suggested that strawberry disease caused by C.theobromicola may affect multiple processes including photosynthesis, sesquiterpenoid and triterpenoid biosynthesis, flavonoid biosynthesis, glutathione metabolism and plant hormone signal transduction. In order to validate our DGE data, 20 DEGs with annotations were selected for qRT-PCR analysis. Among them,19 out of 20 DGEs exhibited a consistent expression pattern between RNA-Seq and qRT-PCR.The comprehensive transcriptome information obtained in this study will help to understand the regulatory mechanisms involved in strawberry response to C.theobromicola infection.

    Key words

    strawberry; Colletotrichum theobromicola; transcriptomics analysis; differentially expressed genes

    炭疽病的病原菌為炭疽菌屬Colletotrichum的多種真菌[1],其中可以侵染草莓的病原菌多是C.acutatum復(fù)合種[2](C.nymphaeae、C.cuscutae、C.godetiae、C.salicis和C.simmondsii)和C.gloeosporioides復(fù)合種[3](C.fructicola和C.theobromicola)。炭疽病可在草莓、蘋果、梨、棗、柑橘和虎耳草等多種植物中發(fā)生,給這些植物的生產(chǎn)造成嚴(yán)重?fù)p失。近年來,炭疽病在草莓生產(chǎn)中發(fā)生日趨嚴(yán)重,已經(jīng)成為草莓苗期甚至整個(gè)生長期的主要病害,是設(shè)施草莓移栽后成活率低的主要原因之一,給草莓生產(chǎn)帶來了嚴(yán)重影響。目前, 生產(chǎn)上針對炭疽病的防治主要有化學(xué)防治及一些微生物防治的方法,前者易造成藥物殘留、環(huán)境污染等問題,后者則往往以預(yù)防為主。種植抗性品種是防治草莓炭疽病最根本的方法,這樣才能從源頭上解決炭疽病侵染[4]。

    植物在逆境脅迫下,轉(zhuǎn)錄水平上的基因表達(dá)會(huì)發(fā)生很大的改變。近年來轉(zhuǎn)錄組測序技術(shù)(RNA-seq)發(fā)展成為研究逆境應(yīng)答、挖掘新基因及代謝途徑確定的重要技術(shù)之一[5]。張明偉等研究發(fā)現(xiàn)野生稻基因?qū)胂礧6023在受到白葉枯菌侵染后,差異表達(dá)基因主要富集在次生代謝產(chǎn)物的生物合成、植物激素信號(hào)轉(zhuǎn)導(dǎo)及糖代謝途徑等方面,并從中挖掘出了多個(gè)抗病基因[6]。劉之慧等發(fā)現(xiàn),不同配比的紅藍(lán)光對草莓組培苗有機(jī)物的積累與代謝、光合作用、葉綠素的產(chǎn)生與代謝、植物激素的產(chǎn)生均有明顯作用[7]。張小芳等發(fā)現(xiàn)耐旱型野生大豆受到干旱脅迫后,其差異表達(dá)基因廣泛涉及糖、脂類、蛋白質(zhì)和核酸等生物大分子代謝、能量代謝以及次生代謝過程,代謝途徑則涉及植物-病原體互作、植物激素信號(hào)轉(zhuǎn)導(dǎo)、RNA降解和ABC轉(zhuǎn)運(yùn)蛋白等[8]。

    本研究用膠孢炭疽菌侵染草莓根部,待草莓發(fā)病后,對健康草莓和患病草莓進(jìn)行轉(zhuǎn)錄組分析,找到二者在轉(zhuǎn)錄水平基因表達(dá)的差異,分析參與植物抗病過程的相關(guān)基因和途徑,為草莓抗病品種的選育提供一些理論依據(jù)。

    1 材料與方法

    1.1 供試草莓和炭疽菌菌株供試草莓品種為‘紅顏。

    供試真菌草莓炭疽菌C.theobromicola TAN-9為中國農(nóng)業(yè)大學(xué)農(nóng)學(xué)與生物技術(shù)學(xué)院張國珍老師贈(zèng)送。

    1.2 培養(yǎng)基

    PDA培養(yǎng)基:馬鈴薯200 g,葡萄糖20 g,瓊脂15~20 g,蒸餾水1 000 mL,自然pH。

    1.3 病原菌侵染試驗(yàn)

    將健康、生長狀態(tài)一致的2葉1心草莓幼苗移植到滅菌基質(zhì)中,處理組(D)和對照組(CK)各15株。待移苗成功10 d后對處理組草莓進(jìn)行傷根處理,并澆灌濃度為1×106 個(gè)/mL草莓炭疽菌TAN-9的孢子懸浮液;對照組澆灌無菌水。此后,結(jié)合草莓生長狀態(tài)和基質(zhì)干濕情況,每隔5~7 d進(jìn)行補(bǔ)施病原菌,一共澆灌3次。

    在處理組接種病原菌50 d后進(jìn)行樣品采集?;疾∪~柄的選擇:截取包含病斑的10 mm長的葉柄段,病斑長度為5.1~10.0 mm。每株草莓只收集1個(gè)葉柄,5株草莓即5個(gè)葉柄段為1組,處理組和對照組各3組重復(fù)。采集到的樣品在液氮迅速冷凍,-80℃冷藏備用。

    1.4 轉(zhuǎn)錄組分析

    1.4.1 草莓葉柄RNA提取

    將采取的樣品用干冰盒包裝好,送上海派森諾生物科技股份有限公司提取 RNA和測序。利用瓊脂糖凝膠電泳和Agilent 2100 Bioanalyzer檢測RNA完整性。

    1.4.2 草莓葉柄cDNA文庫構(gòu)建和質(zhì)控

    用Oligo(dT)磁珠富集總RNA中帶有polyA結(jié)構(gòu)的mRNA,采用離子打斷的方式,將RNA打斷成200~300 bp片段。以RNA為模板,用6堿基隨機(jī)引物和逆轉(zhuǎn)錄酶合成cDNA第一鏈,并以第一鏈cDNA為模板合成第二鏈cDNA,第二鏈cDNA合成時(shí),其中的堿基T,被替換成U,從而達(dá)到鏈特異性文庫的目的。

    文庫構(gòu)建完成后,采用PCR擴(kuò)增進(jìn)行文庫片段富集,之后根據(jù)片段大小進(jìn)行文庫選擇,文庫大小在300~400 bp。接著,通過Agilent 2100 Bioanalyzer對文庫進(jìn)行質(zhì)檢,再對文庫總濃度及文庫有效濃度進(jìn)行檢測。然后根據(jù)文庫的有效濃度以及文庫所需數(shù)據(jù)量,將含有不同Index序列的文庫按比例進(jìn)行混合。混合文庫統(tǒng)一稀釋到2 nmol/L,通過堿變性,形成單鏈文庫。

    1.4.3 高通量測序

    草莓葉柄的轉(zhuǎn)錄組分析: 利用第二代測序技術(shù),基于Illumina HiSeq測序平臺(tái),對庫進(jìn)行雙末端 (paired-end,PE)測序。將測序得到的原始讀序(raw reads)經(jīng)過截除reads中的測序接頭以及引物序列,過濾低質(zhì)量值數(shù)據(jù)等處理,獲得高質(zhì)量的干凈讀序(clean reads),確保數(shù)據(jù)質(zhì)量。進(jìn)行堿基組成和堿基質(zhì)量分析,控制數(shù)據(jù)質(zhì)量。

    提取每個(gè)基因下最長的轉(zhuǎn)錄本作為該基因的代表序列,稱為Unigene。對 Unigene 進(jìn)行基因功能注釋?;蚬δ茏⑨屗玫降臄?shù)據(jù)庫包括NR(NCBI non-redundant protein sequences)、GO (gene ontology)、KEGG(Kyoto encyclopedia of genes and genome)、eggNOG(evolutionary genealogy of genes: non-supervised orthologous groups)、Swiss-Prot。

    1.4.4 差異基因的篩選

    采用DESeq對基因表達(dá)進(jìn)行差異分析,篩選差異表達(dá)基因條件為:表達(dá)差異倍數(shù)|log2FoldChange|>1,顯著性 P-value<0.05。

    1.4.5 差異表達(dá)基因的GO功能顯著性富集分析

    采用BLAST2 GO軟件將差異基因映射到GO數(shù)據(jù)庫的各個(gè)分類,計(jì)算每個(gè)分類中的差異基因數(shù)目,以整個(gè)基因組為背景,采用超幾何分布計(jì)算差異基因顯著富集的GO分類信息。

    1.4.6 差異表達(dá)基因的 KEGG 通路富集顯著性分析

    將差異基因與KEGG數(shù)據(jù)庫進(jìn)行比對,獲得差異表達(dá)基因主要參與的代謝途徑和信號(hào)通路。以整個(gè)基因組為背景,采用超幾何分布計(jì)算差異基因顯著富集的通路。

    1.4.7 差異表達(dá)基因的熒光定量PCR(qRT-PCR)驗(yàn)證

    綜合分析NCBI數(shù)據(jù)庫,并結(jié)合本研究轉(zhuǎn)錄組數(shù)據(jù)提供的信息,對轉(zhuǎn)錄組中20個(gè)可能與植物抗性或者病原菌侵染過程等機(jī)制相關(guān)的差異表達(dá)基因進(jìn)行聚類分析,并對其進(jìn)行qRT-PCR驗(yàn)證。

    利用 Premier 5.0軟件設(shè)計(jì)qPCR引物,以草莓β-actin作為內(nèi)參基因(表1)。用PrimeScriptTM 1st stand cDNA Synthesis Kit將檢測合格并定量的草莓葉柄總RNA逆轉(zhuǎn)錄成 cDNA。熒光定量PCR反應(yīng)體系(20 μL): cDNA 1 μL、上下游引物各 0.8 μL、2×SYBR real-time PCR premixture 10 μL,加超純水補(bǔ)足至20 μL。熒光定量PCR的擴(kuò)增條件為 95℃預(yù)變性5 min,95℃變性 15 s,60℃延伸 30 s,循環(huán) 40 次。每份樣品進(jìn)行3次生物學(xué)重復(fù)。整個(gè)流程由上海派森諾生物科技股份有限公司完成。

    2 結(jié)果與分析

    2.1 病原菌侵染草莓植株

    ‘紅顏草莓在接種病原菌45 d后,葉柄開始發(fā)病,在接種50 d后處理組中超過半數(shù)的草莓葉柄出現(xiàn)了明顯的紅褐色病斑(圖1),此時(shí)采集樣品進(jìn)行分析。

    2.2 測序數(shù)據(jù)結(jié)果統(tǒng)計(jì)

    通過RNA-seq技術(shù)對健康草莓葉柄(CK1、CK2、CK3)和患病草莓葉柄(D1、D2、D3)的6個(gè)樣品進(jìn)行分析,獲得的序列數(shù)據(jù)見表2。

    將上述序列進(jìn)行拼接和分析,共獲得86 988個(gè)Unigene序列。 分別在NR數(shù)據(jù)庫、GO數(shù)據(jù)庫、KEGG數(shù)據(jù)庫、eggNOG數(shù)據(jù)庫和Swiss-Prot數(shù)據(jù)庫中對Unigene序列進(jìn)行比對分析,發(fā)現(xiàn)有4 364個(gè)序列全部獲得注釋(表3)。

    2.3 差異表達(dá)基因的篩選

    將健康草莓葉柄(CK)與患病草莓葉柄(D)的差異基因進(jìn)行雙向聚類分析(圖2),結(jié)果顯示聚類結(jié)果良好。對二者差異基因的表達(dá)統(tǒng)計(jì)分析發(fā)現(xiàn),顯著差異表達(dá)基因的數(shù)目是2 127個(gè),其中上調(diào)表達(dá)數(shù)量是1 021個(gè),下調(diào)表達(dá)基因數(shù)量是1 106個(gè),上調(diào)表達(dá)基因數(shù)量略低于下調(diào)表達(dá)基因。大量的基因在草莓響應(yīng)炭疽病時(shí)呈現(xiàn)差異表達(dá)模式,表明草莓響應(yīng)炭疽病是多基因共同作用的復(fù)雜過程。

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

    對差異表達(dá)基因進(jìn)行GO數(shù)據(jù)庫分析,一共比對到了1 122個(gè)GO功能條目,篩選出p值小于0.01的68個(gè)條目進(jìn)行后續(xù)分析。 這68個(gè)條目分為3大類功能,生物學(xué)過程(biological process)、細(xì)胞組分(cellular component)和分子功能(molecular function),每類功能分別包括26個(gè)、14個(gè)和10個(gè)功能亞組,圖3展示了每類功能之下的10個(gè)功能亞分組。

    分析發(fā)現(xiàn),感染炭疽菌的草莓葉柄的光合作用過程受到了極大的影響,這一結(jié)果和病癥相對應(yīng)。炭疽菌侵染到草莓葉柄處,逐漸積累到一定程度時(shí)會(huì)在葉柄外部出現(xiàn)紅褐色病斑,嚴(yán)重影響正常的光合作用,進(jìn)而草莓的很多生長代謝過程也會(huì)出現(xiàn)異常。分子功能中的氧化還原酶活性功能,還有生物學(xué)過程中的氧化還原過程都富集了大量的差異表達(dá)基因等,這二者都是植物抗逆反應(yīng)中重要的環(huán)節(jié)[9],說明在草莓應(yīng)答炭疽病的過程中,氧化還原反應(yīng)很可能發(fā)揮了關(guān)鍵性的作用。

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

    利用KEGG數(shù)據(jù)庫可以分析差異表達(dá)基因參與的最主要生化代謝途徑和信號(hào)轉(zhuǎn)導(dǎo)途徑。本研究中有556個(gè)DEGs映射到154個(gè)pathway,圖4中展示差異顯著的前20個(gè)通路,包括光合作用、碳代謝、光合生物中的碳固定、光合作用觸角蛋白質(zhì)、甘氨酸、絲氨酸和蘇氨酸代謝、類黃酮生物合成、谷胱甘肽代謝、植物激素信號(hào)轉(zhuǎn)導(dǎo)、牛磺酸和次?;撬岽x、氨基酸的生物合成、苯丙醇生物合成、細(xì)胞色素P450對外源物的代謝、藥物代謝細(xì)胞色素P450、乙醛和二羧酸代謝、甲烷代謝、苯丙氨酸、酪氨酸和色氨酸的生物合成、糖酵解/糖異生、戊糖磷酸途徑以及倍半萜和三萜生物合成。

    從差異表達(dá)基因的富集程度可以看出,其主要還是集中在植物正常生長的代謝途徑中,包括光合作用、牛磺酸和次?;撬岽x、倍半萜和三萜生物途徑等,這表明病原菌的入侵嚴(yán)重干擾了植物的正常生長,尤其是光合作用。同時(shí),與植物抵御非生物及生物脅迫有關(guān)的途徑,包括類黃酮生物合成、谷胱甘肽代謝以及植物激素信號(hào)轉(zhuǎn)導(dǎo)都有差異表達(dá)基因的富集,而在類黃酮途徑中更是有超過三分之一的基因出現(xiàn)了變化,充分說明炭疽菌侵染后,觸發(fā)了草莓體內(nèi)的抗性反應(yīng)[1011],而類黃酮途徑可能起到了關(guān)鍵性的作用。

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

    利用實(shí)時(shí)定量PCR(qRT-PCR)技術(shù)對轉(zhuǎn)錄組數(shù)據(jù)進(jìn)行驗(yàn)證,選取了20個(gè)基因進(jìn)行驗(yàn)證。分析結(jié)果顯示(圖5),有19個(gè)差異表達(dá)基因的表達(dá)趨勢和轉(zhuǎn)錄組測序結(jié)果一致,表明 RNA-Seq 分析結(jié)果可靠。雖然差異基因的表達(dá)倍數(shù)略有不同,但這可能與兩種檢測方法的檢測范圍和靈敏度不同有關(guān)。

    19個(gè)差異表達(dá)基因中包含10個(gè)上調(diào)基因和9個(gè)下調(diào)基因。上調(diào)基因包括:1-氨基環(huán)丙烷-1-羧酸合成酶(ACS)、2-氨基乙硫醇雙加氧酶(ADO)、過氧化物酶21(POD21)、無色花色苷雙加氧酶(LDOX)、查爾酮合酶(FrCHS5)、查爾酮黃酮異構(gòu)酶3(CHI3)、谷胱甘肽S-轉(zhuǎn)移酶F11(GST F11)、MYB、苯丙氨酸解氨酶1(PAL1)和呼吸爆發(fā)氧化酶同源蛋白A(RhoA);下調(diào)基因有:NF-X1型鋅指蛋白(NF-XL1)、色氨酸合酶α鏈(trpA)、生長素誘導(dǎo)蛋白15A(AX15A)、細(xì)胞分裂素脫氫酶5(CKX5)、ADP-葡萄糖焦磷酸化酶大亞基(gpL2)、類胡蘿卜素雙加氧酶4(CCD4)、9-順式環(huán)氧類胡蘿卜素雙加氧酶1(NCED1)、核酮糖1,5-二磷酸羧化酶/加氧酶(rbcL)和葉綠體果糖-1,6-二磷酸酶Ⅰ(cfbI)。

    從表達(dá)量上分析,在上調(diào)和下調(diào)基因中各有3個(gè)表達(dá)量倍數(shù)超過40倍的基因。ACS、POD21和ADO基因,在處理組中表達(dá)倍數(shù)分別是對照組的927、586和326倍;而trpA、AX15A和CKX5在對照組的表達(dá)量則是處理組的49、252和271倍。

    3 討論

    草莓炭疽病是草莓在生長期內(nèi)面臨的主要威脅之一[1213],直接影響其生長發(fā)育,給草莓產(chǎn)量和品質(zhì)帶來極大危害,培育抗炭疽病的草莓品種是防控病害的最直接有效的方法,能從源頭上預(yù)防炭疽病的發(fā)生[14]。病原菌在侵入寄主植物后,二者在共同進(jìn)化及相互識(shí)別過程中形成了錯(cuò)綜復(fù)雜的關(guān)系。病原菌從植物獲得營養(yǎng)的同時(shí)也激發(fā)了植物體內(nèi)一系列生理生化反應(yīng),增強(qiáng)其抗病性,從而達(dá)到抵制病原菌入侵的目的。近年來,越來越多的研究者從轉(zhuǎn)錄組水平上揭示植物抗病系統(tǒng)中的代謝通路以及分子調(diào)控網(wǎng)絡(luò),不僅有助于篩選與鑒定關(guān)鍵作用基因,也為抗病性育種提供了有力支持[1517]。

    本研究利用高通量測序手段,在轉(zhuǎn)錄組水平上分析健康草莓葉柄和患炭疽病草莓葉柄的差異表達(dá)基因,結(jié)果發(fā)現(xiàn)草莓應(yīng)答病原菌的過程中,共有2 127個(gè)基因表達(dá)出現(xiàn)顯著差異,其中有1 021個(gè)差異基因的表達(dá)水平顯著升高,1 106個(gè)差異基因的表達(dá)水平顯著下降。GO數(shù)據(jù)庫和KEGG數(shù)據(jù)庫分析顯示,當(dāng)炭疽病菌在進(jìn)入草莓體內(nèi)后,草莓植株的正常生長受到了極大的影響,多個(gè)代謝通路的基因表達(dá)發(fā)生了改變。根據(jù)轉(zhuǎn)錄組差異基因分析結(jié)果,挑選了20個(gè)基因進(jìn)行qRT-PCR檢測,有19個(gè)基因的表達(dá)趨勢與轉(zhuǎn)錄組結(jié)果一致,說明轉(zhuǎn)錄組數(shù)據(jù)真實(shí)可靠,可以進(jìn)行后續(xù)分析。通過差異基因和代謝通路進(jìn)行分析,發(fā)現(xiàn)本研究中‘紅顏草莓應(yīng)答炭疽菌主要涉及活性氧的平衡調(diào)節(jié)、乙烯的生物合成以及植物正常發(fā)育等生理過程。

    3.1 活性氧平衡與草莓應(yīng)答炭疽病密切相關(guān)

    正常條件下,植物體內(nèi)活性氧 ROS 的生成和清除處于動(dòng)態(tài)平衡狀態(tài),當(dāng)植物遭遇脅迫時(shí),穩(wěn)態(tài)被破壞,大量的 ROS (主要由NADPH氧化酶催化產(chǎn)生,NADPH氧化酶由Rbohs編碼合成[18])迅速累積,使植物細(xì)胞受到氧化損傷[19]。為了維持平衡,植物會(huì)通過兩種方式清除活性氧,第一種是酶促清除系統(tǒng),通過抗氧化酶(SOD、GST、CAT和POD等)將活性氧分解;第二種是非酶系統(tǒng)的抗氧化次級(jí)代謝物,他們可以與活性氧發(fā)生反應(yīng),直接清除ROS,例如多酚、類黃酮、類胡蘿卜素等[20]。

    根據(jù)本研究結(jié)果,我們認(rèn)為草莓接種炭疽菌后,上述兩種方式都會(huì)開始工作,進(jìn)行活性氧的清除。在酶促清除系統(tǒng)中,草莓主要通過體內(nèi)過氧化物氧化酶POD21來進(jìn)行。POD是活性氧酶促清除系統(tǒng)中常見的酶,可以分解體內(nèi)產(chǎn)生的H2O2,它與植物的抗逆境能力有著密切的聯(lián)系。大豆在重度干旱脅迫下,葉片內(nèi)POD活性最高[21];馬鈴薯[22]和番茄[23]在受到病原菌侵染后,葉片POD活性較未接種對照均有顯著提高。在本試驗(yàn)中其表達(dá)量在炭疽菌侵染后升高到對照的586倍,足以證明POD21在活性氧清除過程中發(fā)揮了關(guān)鍵性的作用。

    在第二種方式中,草莓體內(nèi)主要起作用的則是類黃酮物質(zhì)。類黃酮代謝途徑是植物抵抗病原菌的防衛(wèi)反應(yīng)之一,有實(shí)驗(yàn)證實(shí),海島棉在接種枯萎病菌之后,類黃酮代謝通路中的多個(gè)相關(guān)基因表達(dá)量發(fā)生了顯著升高,同時(shí)認(rèn)為此通路與海島棉的抗性相關(guān)[24]。在本試驗(yàn)中,類黃酮代謝途徑(包括上游的苯丙烷類代謝途徑)的相關(guān)基因PAL1、LDOX、FrCHS5與CHI3均在處理組中發(fā)生了上調(diào),因此推測當(dāng)草莓受到病原菌侵染,上述基因發(fā)生上調(diào),加速了類黃酮物質(zhì)[10]的積累,從而清除ROS,維持ROS的平衡,緩解了炭疽菌給植物帶來的氧化傷害,增強(qiáng)了草莓對病菌的抗性。

    值得注意的是,草莓被炭疽菌侵染后,2-氨基乙硫醇雙加氧酶ADO的表達(dá)量上升到對照的326倍。ADO可以將半胱氨酸氧化成次牛磺酸(hypotaurine, HT),后者在動(dòng)物體內(nèi)是一種強(qiáng)抗氧化劑[25];盡管次?;撬嵩谥参镏械难芯亢苌?,但是也有類似的作用。有研究證實(shí),水稻幼苗經(jīng)過HT的處理能夠降低堿脅迫下葉片中活性氧的積累,緩解其對幼苗生長的抑制作用[26]。ADO表達(dá)量的大量增加,很可能預(yù)示著其在應(yīng)答過程中會(huì)起到一定的關(guān)鍵作用,所以說其作用機(jī)制是值得關(guān)注的一個(gè)新的研究點(diǎn),它能否對活性氧平衡發(fā)揮作用則需要進(jìn)一步的研究探索。

    3.2 ACS在草莓應(yīng)答炭疽病中發(fā)揮作用

    乙烯是植物體內(nèi)一種重要的氣態(tài)激素,參與調(diào)節(jié)植物生長發(fā)育的多個(gè)過程,也在植物響應(yīng)生物和非生物脅迫過程中起重要的調(diào)控作用。ACS是乙烯生物合成過程的關(guān)鍵調(diào)控酶,它催化SAM(S-腺苷甲硫氨酸)合成乙烯前體ACC(1-氨基環(huán)丙烷-1-羧酸),ACC在 ACC氧化酶的作用下氧化生成乙烯。耿夙的研究發(fā)現(xiàn),在輪紋病菌侵染蘋果果實(shí)后,ACS5a、ACS5b、ACS4被大量誘導(dǎo)表達(dá),乙烯也被大量誘導(dǎo)合成[27]。這和本研究結(jié)果類似,當(dāng)炭疽菌侵染草莓后,ACS基因的表達(dá)量驟然上升到對照的927倍,是所有差異表達(dá)基因中表達(dá)倍數(shù)最高的基因。推測ACS基因的顯著上調(diào)會(huì)導(dǎo)致ACC的過表達(dá),進(jìn)而促進(jìn)乙烯的合成,而乙烯的累積會(huì)促進(jìn)草莓體內(nèi)逆境脅迫的信號(hào)轉(zhuǎn)導(dǎo)過程,起到防御炭疽菌侵害的作用[28]。

    3.3 植物正常生長受到抑制

    色氨酸合酶由兩個(gè)不同的亞基組成,α-亞基和β-亞基依次催化色氨酸生物合成的最后兩個(gè)步驟[29],并且色氨酸合酶還參與了甘氨酸、絲氨酸和蘇氨酸代謝以及苯丙氨酸和酪氨酸的生物合成;生長素誘導(dǎo)蛋白基因編碼合成植物特有的生長素響應(yīng)蛋白 SAUR (small auxin up RNA),而且此基因在生長旺盛的部位表達(dá)量要高于其他部位[30];細(xì)胞分裂素脫氫酶是分解細(xì)胞分裂素的唯一降解酶[3132]。因此,我們推測,當(dāng)炭疽菌侵染草莓后,草莓體內(nèi)涉及正常生長的一些基因發(fā)生了重大改變。色氨酸合酶 α鏈基因trpA的下調(diào),表明色氨酸合酶α-亞基會(huì)受到影響,進(jìn)而影響植物體內(nèi)包括色氨酸在內(nèi)的多個(gè)氨基酸的合成和代謝;AX15A的顯著下調(diào)會(huì)導(dǎo)致SAUR的缺乏,無法響應(yīng)生長素信號(hào);而CKX 5表達(dá)量的下降,導(dǎo)致細(xì)胞分裂素在草莓體內(nèi)失去代謝平衡。所以,這3個(gè)基因的顯著下調(diào),影響了植物激素和氨基酸代謝途徑,最終導(dǎo)致草莓的生長受到阻礙,出現(xiàn)植株矮小、根系不發(fā)達(dá)等現(xiàn)象。

    本文通過RNA-seq分析和qRT-PCR驗(yàn)證發(fā)現(xiàn)了一些與草莓應(yīng)答炭疽病相關(guān)的差異表達(dá)基因,其中表達(dá)量有顯著差異的基因有6個(gè),分析推測這些差異基因可能在活性氧平衡、乙烯合成以及植物正常生長過程發(fā)揮了重要作用。那么,這些基因在草莓應(yīng)答病原菌過程中是怎樣發(fā)揮作用的,其調(diào)控機(jī)制又是什么,它們之間又會(huì)形成怎樣的關(guān)聯(lián)網(wǎng)絡(luò),尤其是次?;撬嵩谥参锘钚匝跚宄^程中所起到的作用是十分值得研究的,所以我們需要進(jìn)一步開展分子試驗(yàn)來解決這些深層次的問題。轉(zhuǎn)錄組分析為深入研究植物病原菌互作提供了良好的平臺(tái)[33],可以從整體上了解患病草莓與健康草莓在轉(zhuǎn)錄組水平上的差異表達(dá)基因,有助于挖掘草莓防御炭疽病菌的關(guān)鍵基因,為研究草莓應(yīng)答炭疽病的分子機(jī)理提供重要的基因資源,從而為抗性草莓品種的選育提供有力的科學(xué)理論。

    參考文獻(xiàn)

    [1] 陳哲, 黃靜, 趙佳, 等. 防治草莓炭疽病的芽胞桿菌組合的篩選[J]. 中國生物防治學(xué)報(bào), 2018, 34(4): 582588.

    [2] DAMM U, CANNON P F, WOUDENBERG J H C, et al. The Colletotrichum acutatum species complex [J]. Studies in Mycology, 2012, 73(1):37113.

    [3] WEIR B S, JOHNSTON P R, DAMM U. The Colletotrichum gloeosporioides species complex [J]. Studies in Mycology, 2012, 73(1):115180.

    [4] CHEN Jing, ZHANG Hanping, FENG Mingfeng, et al. Transcriptome analysis of woodland strawberry (Fragaria vesca) response to the infection by strawberry vein banding virus (SVBV) [J]. Virology Journal, 2016, 13(1): 128.

    [5] 張健, 唐露, 張雅潔, 等. 轉(zhuǎn)錄組測序技術(shù)在植物水淹脅迫研究中的應(yīng)用[J]. 分子植物育種: 2019(2): 113.

    [6] 張明偉, 徐飛飛, 郝巍, 等. 野生稻基因?qū)胂礧6023對白葉枯菌的抗譜及轉(zhuǎn)錄組差異表達(dá)基因分析[J]. 植物遺傳資源學(xué)報(bào), 2017, 18(2): 298309.

    [7] 劉之慧, 蔡雯婷, 董清華. 紅藍(lán)光處理草莓組培苗的轉(zhuǎn)錄組測序與分析[J]. 北京農(nóng)學(xué)院學(xué)報(bào), 2017, 32(4): 2835.

    [8] 張小芳, 王冰冰, 徐燕, 等. PEG模擬干旱脅迫下野生大豆轉(zhuǎn)錄組分析[J]. 大豆科學(xué), 2018, 37(5): 681689.

    [9] 劉洪博, 劉新龍, 蘇火生, 等. 干旱脅迫下割手密根系轉(zhuǎn)錄組差異表達(dá)分析[J]. 中國農(nóng)業(yè)科學(xué), 2017, 50(6): 11671178.

    [10]潘媛, 陳大霞, 宋旭紅, 等. 基于轉(zhuǎn)錄組測序挖掘國槐花蕾黃酮類物質(zhì)生物合成相關(guān)基因[J]. 中國中藥雜志, 2018, 43(13): 26822689.

    [11]張亞真, 張芬, 王麗鴛, 等. 植物谷胱甘肽轉(zhuǎn)移酶在類黃酮累積中的作用[J]. 植物生理學(xué)報(bào), 2015, 51(11): 18151820.

    [12]谷春艷, 蘇賢巖, 楊雪, 等. 解淀粉芽胞桿菌WH1G與咪鮮胺協(xié)同防治草莓炭疽病[J]. 植物保護(hù), 2018, 44(2): 184189.

    [13]馬世龍, 王國良, 趙君, 等. 大棚‘紅頰草莓病害控制關(guān)鍵技術(shù)[J]. 現(xiàn)代園藝, 2017(6): 6566.

    [14]NAMAI K, MATSUSHIMA Y, MORISHIMA M. Resistance to anthracnose is decreased by tissue culture but increased with longer acclimation in the resistant strawberry cultivar [J]. Journal of General Plant Pathology, 2013, 79(6): 402411.

    [15]梁國平, 李文芳, 馬宗桓, 等. 基于轉(zhuǎn)錄組研究葡萄糖對葡萄試管苗碳代謝的影響[J]. 園藝學(xué)報(bào), 2019, 46(1): 121134.

    [16]王冕. 花生應(yīng)答果腐病病原菌侵染的轉(zhuǎn)錄組和miRNAs研究及抗性相關(guān)基因克隆[D]. 泰安: 山東農(nóng)業(yè)大學(xué), 2018.

    [17]楊仕美, 喬光, 毛永亞, 等. 基于火龍果轉(zhuǎn)錄組測序的SSR標(biāo)記開發(fā)及種質(zhì)親緣關(guān)系分析[J]. 分子植物育種, 2018, 16(24): 80968110.

    [18]ZHANG Yunting, LI Yali, HE Yuwei, et al. Identification of NADPH oxidase family members associated with cold stress in strawberry [J]. FEBS Open Bio, 2018, 8(4): 593605.

    [19]郭鵬飛, 雷健, 羅佳佳, 等. 柱花草苯丙氨酸解氨酶(SgPALs)對生物脅迫與非生物脅迫的響應(yīng)[J].熱帶作物學(xué)報(bào),2019,40(9):17421751.

    [20]王福祥, 肖開轉(zhuǎn), 姜身飛, 等. 干旱脅迫下植物體內(nèi)ROS的作用機(jī)制[J].科學(xué)通報(bào),2019,64(17):17651779.

    [21]王文佳,李爽,馬澤眾,等.水分脅迫對春大豆葉片保護(hù)酶活性及相對電導(dǎo)率的影響[J].中國農(nóng)學(xué)通報(bào),2019,35(11):1418.

    [22]胡耀華, 李清宇, 唐翊. 基于高光譜的晚疫病脅迫下馬鈴薯葉片的過氧化物酶活性檢測[J]. 江蘇大學(xué)學(xué)報(bào)(自然科學(xué)版), 2018, 39(6): 683688.

    [23]張菁, 連清貴, 陳婧, 等. 淡紫擬青霉QLP12對感染灰霉病后番茄植株的促生作用及抗病相關(guān)酶活性變化[J]. 植物保護(hù)學(xué)報(bào), 2018, 45(5): 10881095.

    [24]黃啟秀,曲延英,姚正培,等. 海島棉枯萎病抗性與類黃酮代謝途徑基因表達(dá)量的相關(guān)性[J].作物學(xué)報(bào),2017,43(12):17911801.

    [25]TERRIENTE-PALACIONS C, DIAZ I, CASTELLARI M. A validated ultra-performance liquid chromatography with diode array detection coupled to electrospray ionization and triple quadrupole mass spectrometry method to simultaneously quantify taurine, homotaurine, hypotaurine and amino acids in macro-and microalgae [J].Journal of Chromatography A, 2019,1589:8392.

    [26]祝一文,車永梅,趙方貴,等.堿脅迫下H2S參與活性氧代謝和水稻幼苗生長的調(diào)控[J].農(nóng)業(yè)生物技術(shù)學(xué)報(bào),2018,26(7):11241131.

    [27]耿夙. 輪紋病致病菌Botryosphaeria dothidea侵染對蘋果抗性反應(yīng)、乙烯合成和相關(guān)基因表達(dá)的影響[D].泰安:山東農(nóng)業(yè)大學(xué),2016.

    [28]張弢, 董春海. 乙烯信號(hào)轉(zhuǎn)導(dǎo)及其在植物逆境響應(yīng)中的作用[J]. 生物技術(shù)通報(bào), 2016, 32(10): 1117.

    [29]FLEMING J R, SCHUPFNER M, BUSCH F, et al. Evolutionary morphing of tryptophan synthase: functional mechanisms for the enzymatic channeling of indole [J]. Journal of Molecular Biology, 2018:50675079.

    [30]張歡,楊英杰,李鼎立,等.杜梨根莖葉特異表達(dá)基因的RNA-Seq分析[J].園藝學(xué)報(bào),2018,45(10):18811894.

    [31]張晨, 鄭偉, 周濤, 等. 太子參細(xì)胞分裂素氧化/脫氫酶基因PhCKX的克隆及表達(dá)分析[J]. 分子植物育種, 2018, 16(19): 62426249.

    [32]李志康,嚴(yán)冬,薛張逸,等.細(xì)胞分裂素對植物生長發(fā)育的調(diào)控機(jī)理研究進(jìn)展及其在水稻生產(chǎn)中的應(yīng)用探討[J].中國水稻科學(xué),2018,32(4):311324.

    [33]ZANARDO L G, DE SOUZA G B, ALVES M S. Transcriptomics of plant-virus interactions: a review [J]. Theoretical and Experimental Plant Physiology. 2019, 31:103125.

    (責(zé)任編輯:田 喆)

    收稿日期: 20190302?? 修訂日期: 20190511

    基金項(xiàng)目:山西省農(nóng)業(yè)科學(xué)院特色農(nóng)業(yè)技術(shù)攻關(guān)項(xiàng)目(YGG17067);山西省應(yīng)用基礎(chǔ)研究項(xiàng)目(201601D202062)

    致? 謝: 參加本試驗(yàn)部分工作的還有江代禮、譚翰杰、張能和紀(jì)燁斌等同學(xué),特此一并致謝。

    通信作者E-mail:lh1964@126.com

    #為并列第一作者

    猜你喜歡
    差異表達(dá)基因草莓
    愛提問的小草莓
    鹽脅迫對花生硝酸鹽積累及信號(hào)轉(zhuǎn)導(dǎo)的影響
    草莓
    生物信息學(xué)分析患有乳腺癌的乳腺球樣本中與自我更新相關(guān)的關(guān)鍵基因
    草莓
    空腹吃草莓
    條斑紫菜優(yōu)良品系的基因芯片表達(dá)譜分析
    高溫脅迫下草坪草高羊茅差異表達(dá)基因的分子研究
    采草莓
    抑制消減雜交技術(shù)在瓜菜作物研究中的應(yīng)用
    中國瓜菜(2009年4期)2009-09-23 08:46:18
    av天堂在线播放| 亚洲熟妇熟女久久| 亚洲av成人不卡在线观看播放网| 麻豆成人av在线观看| 岛国在线观看网站| 亚洲成av人片在线播放无| 国产高清三级在线| av黄色大香蕉| 亚洲欧美精品综合久久99| 成人性生交大片免费视频hd| 欧美日韩乱码在线| 欧美黑人欧美精品刺激| 亚洲激情在线av| avwww免费| 此物有八面人人有两片| 午夜两性在线视频| 国产视频一区二区在线看| 国产高清视频在线播放一区| 在线a可以看的网站| 免费av不卡在线播放| 国产伦一二天堂av在线观看| 午夜免费男女啪啪视频观看 | 欧美另类亚洲清纯唯美| 亚洲一区高清亚洲精品| 国产美女午夜福利| 一进一出抽搐动态| 国产美女午夜福利| aaaaa片日本免费| 欧美成狂野欧美在线观看| 午夜福利在线观看吧| 亚洲国产精品成人综合色| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品影视一区二区三区av| 久久精品91无色码中文字幕| 午夜亚洲福利在线播放| av视频在线观看入口| 午夜福利在线观看吧| 99精品久久久久人妻精品| 51国产日韩欧美| 老司机午夜十八禁免费视频| 亚洲成a人片在线一区二区| 在线观看日韩欧美| 免费在线观看亚洲国产| 亚洲av不卡在线观看| 精品久久久久久成人av| 国内久久婷婷六月综合欲色啪| 亚洲avbb在线观看| 国产久久久一区二区三区| 一本精品99久久精品77| 九色国产91popny在线| 日本一本二区三区精品| 看黄色毛片网站| 中文资源天堂在线| 久久精品亚洲精品国产色婷小说| 一本综合久久免费| 亚洲五月婷婷丁香| 日日干狠狠操夜夜爽| 一个人观看的视频www高清免费观看| 丝袜美腿在线中文| 日本 欧美在线| 男人和女人高潮做爰伦理| 免费在线观看成人毛片| 欧美中文日本在线观看视频| 日韩欧美在线二视频| 中文字幕人妻熟人妻熟丝袜美 | 成人三级黄色视频| 久久99热这里只有精品18| 久久久久久久久久黄片| 国产国拍精品亚洲av在线观看 | 两人在一起打扑克的视频| 色综合欧美亚洲国产小说| 最新美女视频免费是黄的| 天堂网av新在线| av专区在线播放| 天天一区二区日本电影三级| 婷婷精品国产亚洲av| 国产成人av激情在线播放| 日韩大尺度精品在线看网址| 国产伦在线观看视频一区| 波多野结衣巨乳人妻| 香蕉av资源在线| or卡值多少钱| 在线播放无遮挡| 亚洲欧美日韩东京热| av天堂中文字幕网| 国产精品影院久久| 99久久九九国产精品国产免费| 身体一侧抽搐| 国产真人三级小视频在线观看| 国产午夜精品久久久久久一区二区三区 | 亚洲无线观看免费| 国产色婷婷99| 韩国av一区二区三区四区| 夜夜看夜夜爽夜夜摸| 亚洲美女视频黄频| 岛国视频午夜一区免费看| 舔av片在线| 免费搜索国产男女视频| 免费观看的影片在线观看| 成年免费大片在线观看| 日韩欧美国产在线观看| 欧美日韩国产亚洲二区| 国产乱人伦免费视频| 精品一区二区三区人妻视频| 97人妻精品一区二区三区麻豆| 男人舔女人下体高潮全视频| 嫩草影院精品99| 给我免费播放毛片高清在线观看| www国产在线视频色| 久久人妻av系列| 中文资源天堂在线| 中文亚洲av片在线观看爽| 国产成人系列免费观看| 国产精品一区二区三区四区久久| 男女视频在线观看网站免费| 成年免费大片在线观看| 国产探花在线观看一区二区| 日韩国内少妇激情av| 九九在线视频观看精品| 男女下面进入的视频免费午夜| 国产一区二区亚洲精品在线观看| 国产精品久久视频播放| aaaaa片日本免费| 久久伊人香网站| 免费观看精品视频网站| 在线播放无遮挡| 亚洲人成网站在线播放欧美日韩| 日韩精品青青久久久久久| 草草在线视频免费看| 19禁男女啪啪无遮挡网站| 久久久久久久精品吃奶| 日韩亚洲欧美综合| 老司机午夜福利在线观看视频| 天堂网av新在线| 精品福利观看| 在线观看av片永久免费下载| 国产不卡一卡二| 国产欧美日韩一区二区三| 亚洲狠狠婷婷综合久久图片| 免费搜索国产男女视频| x7x7x7水蜜桃| 久久久久久久亚洲中文字幕 | 91在线观看av| 在线观看免费视频日本深夜| 国产精品综合久久久久久久免费| 无遮挡黄片免费观看| 长腿黑丝高跟| 午夜福利在线观看吧| xxxwww97欧美| 日韩欧美在线二视频| 脱女人内裤的视频| 老熟妇乱子伦视频在线观看| 亚洲 国产 在线| 亚洲美女视频黄频| 偷拍熟女少妇极品色| 成年人黄色毛片网站| 亚洲欧美精品综合久久99| 午夜精品一区二区三区免费看| 亚洲av中文字字幕乱码综合| 欧美激情久久久久久爽电影| 大型黄色视频在线免费观看| 男女做爰动态图高潮gif福利片| 19禁男女啪啪无遮挡网站| 亚洲黑人精品在线| 99在线视频只有这里精品首页| 亚洲国产中文字幕在线视频| 丁香六月欧美| 在线观看66精品国产| 国产午夜精品久久久久久一区二区三区 | 免费看光身美女| 亚洲国产高清在线一区二区三| 日韩欧美精品v在线| 欧美午夜高清在线| 亚洲av第一区精品v没综合| 在线观看日韩欧美| 国产在视频线在精品| 欧美日韩国产亚洲二区| 我的老师免费观看完整版| 国产免费av片在线观看野外av| 欧美一区二区国产精品久久精品| 少妇裸体淫交视频免费看高清| 国产亚洲精品久久久com| 国产精品亚洲av一区麻豆| 国产在线精品亚洲第一网站| 成人国产综合亚洲| 97碰自拍视频| 日韩欧美三级三区| 久久亚洲真实| а√天堂www在线а√下载| 男人和女人高潮做爰伦理| 欧美另类亚洲清纯唯美| netflix在线观看网站| 九九在线视频观看精品| 日本 av在线| 我要搜黄色片| 在线观看日韩欧美| 国产色爽女视频免费观看| 成年版毛片免费区| 午夜影院日韩av| 床上黄色一级片| 免费人成视频x8x8入口观看| 午夜精品在线福利| 日本在线视频免费播放| 老司机深夜福利视频在线观看| 日韩精品青青久久久久久| 一二三四社区在线视频社区8| 99久久无色码亚洲精品果冻| 久久久国产成人免费| 国产精品永久免费网站| 久久久久久久久中文| 日本黄色视频三级网站网址| 搡老岳熟女国产| 国产亚洲欧美在线一区二区| 欧美xxxx黑人xx丫x性爽| 欧美3d第一页| 少妇的逼水好多| av在线蜜桃| 天堂动漫精品| 国产精品久久视频播放| 成熟少妇高潮喷水视频| 亚洲中文日韩欧美视频| 午夜福利在线观看免费完整高清在 | 免费在线观看亚洲国产| 男女那种视频在线观看| 伊人久久精品亚洲午夜| 日韩欧美精品v在线| 一个人免费在线观看的高清视频| 亚洲欧美日韩高清专用| 欧美日韩亚洲国产一区二区在线观看| 欧美日韩国产亚洲二区| 人人妻,人人澡人人爽秒播| 99在线人妻在线中文字幕| 中国美女看黄片| 高清在线国产一区| 激情在线观看视频在线高清| 亚洲乱码一区二区免费版| 99热只有精品国产| 嫩草影院入口| 久久久久性生活片| 一本综合久久免费| 婷婷精品国产亚洲av在线| 国产一区二区亚洲精品在线观看| 亚洲av免费高清在线观看| 两个人看的免费小视频| 成人特级黄色片久久久久久久| 少妇丰满av| 成年女人永久免费观看视频| 国产欧美日韩一区二区精品| 日日夜夜操网爽| 身体一侧抽搐| 国产黄a三级三级三级人| 麻豆成人午夜福利视频| 免费观看精品视频网站| 身体一侧抽搐| 欧美日韩中文字幕国产精品一区二区三区| 51国产日韩欧美| 一级毛片高清免费大全| 国产99白浆流出| 深夜精品福利| 国产精品嫩草影院av在线观看 | 国产午夜福利久久久久久| 深夜精品福利| 免费观看的影片在线观看| xxx96com| 成年女人看的毛片在线观看| 99久久精品热视频| 日韩高清综合在线| 两个人看的免费小视频| 成人特级黄色片久久久久久久| 亚洲美女视频黄频| 老鸭窝网址在线观看| 熟女少妇亚洲综合色aaa.| 欧美日韩瑟瑟在线播放| 中文字幕人成人乱码亚洲影| 露出奶头的视频| 日本撒尿小便嘘嘘汇集6| 搞女人的毛片| av欧美777| 91在线精品国自产拍蜜月 | 亚洲成人久久性| 国产成年人精品一区二区| tocl精华| 人妻丰满熟妇av一区二区三区| 国产极品精品免费视频能看的| 91麻豆精品激情在线观看国产| 日本 欧美在线| 在线播放国产精品三级| 少妇高潮的动态图| 国产一区在线观看成人免费| 久久久久久久精品吃奶| 国产精品久久电影中文字幕| 麻豆成人av在线观看| 亚洲av美国av| 性色avwww在线观看| 日韩欧美国产在线观看| 老司机在亚洲福利影院| 99久久精品热视频| 深爱激情五月婷婷| 精品午夜福利视频在线观看一区| 国产精品乱码一区二三区的特点| av中文乱码字幕在线| 亚洲18禁久久av| 18禁裸乳无遮挡免费网站照片| 久久人妻av系列| 女同久久另类99精品国产91| 亚洲av美国av| 午夜精品在线福利| 波多野结衣巨乳人妻| 国产午夜精品久久久久久一区二区三区 | 99热这里只有是精品50| 激情在线观看视频在线高清| 亚洲精品日韩av片在线观看 | 免费看日本二区| 亚洲最大成人中文| 国产精品久久久人人做人人爽| 青草久久国产| 窝窝影院91人妻| 国产精品久久久久久亚洲av鲁大| 一个人免费在线观看的高清视频| 成年女人毛片免费观看观看9| 内射极品少妇av片p| 欧美bdsm另类| 18美女黄网站色大片免费观看| 级片在线观看| 日本黄色视频三级网站网址| 中文字幕人妻丝袜一区二区| 国产成人aa在线观看| 欧美日韩国产亚洲二区| 国产一区在线观看成人免费| 俺也久久电影网| 亚洲人成网站高清观看| 精品电影一区二区在线| 国产精品自产拍在线观看55亚洲| 老司机福利观看| 精品不卡国产一区二区三区| 亚洲,欧美精品.| 国产精品美女特级片免费视频播放器| 欧美乱码精品一区二区三区| 国产真实伦视频高清在线观看 | 欧美三级亚洲精品| 一二三四社区在线视频社区8| 在线观看日韩欧美| 欧美日本视频| 亚洲国产精品sss在线观看| 一个人观看的视频www高清免费观看| 蜜桃久久精品国产亚洲av| 全区人妻精品视频| 久久精品影院6| 美女高潮的动态| 国产乱人伦免费视频| 88av欧美| 夜夜夜夜夜久久久久| 又黄又爽又免费观看的视频| 国产精品久久久人人做人人爽| 精品熟女少妇八av免费久了| 欧美+日韩+精品| 少妇丰满av| 成人性生交大片免费视频hd| 免费人成在线观看视频色| 成人国产一区最新在线观看| 色播亚洲综合网| 日本在线视频免费播放| 国内精品美女久久久久久| 丰满乱子伦码专区| 国产单亲对白刺激| av天堂中文字幕网| 天堂动漫精品| 色视频www国产| 一级黄色大片毛片| 国产私拍福利视频在线观看| 午夜老司机福利剧场| 看免费av毛片| 国产午夜精品久久久久久一区二区三区 | 亚洲中文字幕一区二区三区有码在线看| 国产精品一区二区三区四区久久| 亚洲成人久久爱视频| 亚洲欧美日韩东京热| 黄色女人牲交| 女警被强在线播放| 亚洲国产精品合色在线| 亚洲精品美女久久久久99蜜臀| 精品日产1卡2卡| 国产免费一级a男人的天堂| 极品教师在线免费播放| 精品久久久久久久久久久久久| 国产99白浆流出| 亚洲av熟女| 99精品在免费线老司机午夜| 波多野结衣高清无吗| 精品免费久久久久久久清纯| 两个人看的免费小视频| 日韩欧美国产一区二区入口| 免费搜索国产男女视频| 日韩大尺度精品在线看网址| 男人舔女人下体高潮全视频| 麻豆国产av国片精品| 一本综合久久免费| 久久中文看片网| 国产亚洲欧美98| 女同久久另类99精品国产91| 噜噜噜噜噜久久久久久91| 欧美一级毛片孕妇| 熟女电影av网| 90打野战视频偷拍视频| www国产在线视频色| 又紧又爽又黄一区二区| 国产精品三级大全| e午夜精品久久久久久久| www日本在线高清视频| 美女高潮喷水抽搐中文字幕| 午夜精品久久久久久毛片777| 国产伦一二天堂av在线观看| 黄片小视频在线播放| 日韩有码中文字幕| 国产毛片a区久久久久| 亚洲18禁久久av| 欧美一级毛片孕妇| 久久精品91蜜桃| 女人十人毛片免费观看3o分钟| 国产在线精品亚洲第一网站| 国产真实伦视频高清在线观看 | 午夜两性在线视频| 国产黄色小视频在线观看| 成人亚洲精品av一区二区| 99久久久亚洲精品蜜臀av| 国产精品 国内视频| 在线播放无遮挡| 97超视频在线观看视频| 女人被狂操c到高潮| 精品一区二区三区视频在线观看免费| 国产黄色小视频在线观看| 天堂网av新在线| 亚洲片人在线观看| 国产成人aa在线观看| 成人国产综合亚洲| 国产91精品成人一区二区三区| 久99久视频精品免费| 无限看片的www在线观看| 亚洲精品美女久久久久99蜜臀| www国产在线视频色| 国产午夜福利久久久久久| 国产三级黄色录像| 高清在线国产一区| 2021天堂中文幕一二区在线观| 俺也久久电影网| 久久久久国产精品人妻aⅴ院| 亚洲熟妇中文字幕五十中出| 日本与韩国留学比较| 国产欧美日韩一区二区精品| 亚洲最大成人中文| 国产一区二区在线av高清观看| 特大巨黑吊av在线直播| 黄色日韩在线| 免费看美女性在线毛片视频| 99热精品在线国产| 国产不卡一卡二| 亚洲成a人片在线一区二区| 九九在线视频观看精品| 一夜夜www| 看片在线看免费视频| 免费av毛片视频| 免费看日本二区| 99久久久亚洲精品蜜臀av| 国产午夜精品久久久久久一区二区三区 | 国产免费一级a男人的天堂| 国产欧美日韩精品亚洲av| 日本在线视频免费播放| 2021天堂中文幕一二区在线观| 精品熟女少妇八av免费久了| 国产精品1区2区在线观看.| 性色av乱码一区二区三区2| a级一级毛片免费在线观看| 波多野结衣高清作品| 国产成人啪精品午夜网站| 国产成人系列免费观看| 精品熟女少妇八av免费久了| 中文亚洲av片在线观看爽| 国产精品久久久久久人妻精品电影| 欧美日韩国产亚洲二区| 哪里可以看免费的av片| 国产三级黄色录像| 免费电影在线观看免费观看| 一个人免费在线观看电影| 老司机福利观看| 日韩人妻高清精品专区| 国内毛片毛片毛片毛片毛片| 一级毛片高清免费大全| 亚洲精品粉嫩美女一区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 99精品在免费线老司机午夜| 国产精品久久久久久精品电影| 欧美丝袜亚洲另类 | 精品欧美国产一区二区三| 中文亚洲av片在线观看爽| www.www免费av| 亚洲专区国产一区二区| 午夜亚洲福利在线播放| 国产精品电影一区二区三区| 精品国产美女av久久久久小说| 蜜桃亚洲精品一区二区三区| 国产爱豆传媒在线观看| 嫩草影院入口| 成人高潮视频无遮挡免费网站| 偷拍熟女少妇极品色| 成人精品一区二区免费| 免费高清视频大片| 两个人看的免费小视频| 国产精品av视频在线免费观看| 午夜久久久久精精品| 国产三级在线视频| 精品一区二区三区视频在线 | 精品久久久久久久末码| 欧美日韩国产亚洲二区| 综合色av麻豆| 麻豆一二三区av精品| 亚洲精品美女久久久久99蜜臀| 精品人妻一区二区三区麻豆 | 国产极品精品免费视频能看的| 国产伦人伦偷精品视频| 一级毛片高清免费大全| 国产黄a三级三级三级人| 天堂网av新在线| 岛国视频午夜一区免费看| 丰满人妻熟妇乱又伦精品不卡| 美女大奶头视频| 色视频www国产| 欧美日本亚洲视频在线播放| 午夜激情福利司机影院| 精品不卡国产一区二区三区| 69av精品久久久久久| 国产精品99久久久久久久久| 一个人免费在线观看电影| 国产精品一区二区三区四区免费观看 | 88av欧美| svipshipincom国产片| 亚洲片人在线观看| 精品午夜福利视频在线观看一区| 久久人人精品亚洲av| x7x7x7水蜜桃| 欧美+亚洲+日韩+国产| 中文亚洲av片在线观看爽| 日韩欧美三级三区| 欧美大码av| 国产免费av片在线观看野外av| 成年女人毛片免费观看观看9| 最好的美女福利视频网| 在线观看免费午夜福利视频| 亚洲 欧美 日韩 在线 免费| 久久天躁狠狠躁夜夜2o2o| 黄色日韩在线| 欧美成人a在线观看| 婷婷精品国产亚洲av| 国产精品久久视频播放| 亚洲av二区三区四区| 国产精品久久视频播放| 黄色成人免费大全| 日本 欧美在线| av福利片在线观看| 男女视频在线观看网站免费| 又黄又爽又免费观看的视频| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产精品合色在线| 丝袜美腿在线中文| 亚洲精品一卡2卡三卡4卡5卡| 两个人视频免费观看高清| 99久久九九国产精品国产免费| 一个人观看的视频www高清免费观看| 亚洲欧美日韩东京热| 亚洲av第一区精品v没综合| 欧美精品啪啪一区二区三区| 99久久无色码亚洲精品果冻| 黄色日韩在线| www.www免费av| 动漫黄色视频在线观看| 一本一本综合久久| 一区二区三区高清视频在线| 一级作爱视频免费观看| 日韩免费av在线播放| 中文字幕精品亚洲无线码一区| 日韩欧美国产一区二区入口| 国产亚洲欧美98| 国产真实伦视频高清在线观看 | 色综合亚洲欧美另类图片| 欧美国产日韩亚洲一区| 国产精品免费一区二区三区在线| 亚洲在线自拍视频| 中文字幕人妻丝袜一区二区| 午夜精品一区二区三区免费看| 床上黄色一级片| 亚洲国产欧美网| 性色av乱码一区二区三区2| 亚洲av五月六月丁香网| 日韩高清综合在线| 免费在线观看影片大全网站| 99精品久久久久人妻精品| 日韩精品青青久久久久久| 97超级碰碰碰精品色视频在线观看| 国内少妇人妻偷人精品xxx网站| 亚洲人与动物交配视频| 亚洲精品在线观看二区| www国产在线视频色| 久久精品夜夜夜夜夜久久蜜豆| 久久久久精品国产欧美久久久| 免费电影在线观看免费观看| 国产在线精品亚洲第一网站| 久久精品综合一区二区三区| 国产熟女xx| 超碰av人人做人人爽久久 | 欧美av亚洲av综合av国产av| 熟女少妇亚洲综合色aaa.| 久久精品国产亚洲av香蕉五月| 国产综合懂色| 国产黄a三级三级三级人| 在线播放无遮挡| 99国产精品一区二区蜜桃av|