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

    基于WGCNA 的刺葡萄抗白腐病關(guān)鍵基因的發(fā)掘

    2023-06-12 18:35:14譚西北李鵬孫磊樊秀彩劉崇懷姜建福張穎
    果樹學(xué)報(bào) 2023年4期
    關(guān)鍵詞:轉(zhuǎn)錄組抗病葡萄

    譚西北 李鵬 孫磊 樊秀彩 劉崇懷 姜建福 張穎

    摘要:【目的】葡萄白腐病是危害最嚴(yán)重的病害之一,不同的葡萄種質(zhì)對(duì)白腐病的抗性存在差異。在前期的研究中篩選到了1 份抗病種質(zhì)刺葡萄0941,該種質(zhì)為中國野生葡萄,具有耐濕熱、抗病性強(qiáng)的特點(diǎn)。旨在通過轉(zhuǎn)錄組測(cè)序和WGCNA篩選,挖掘刺葡萄抗白腐病關(guān)鍵基因?!痉椒ā恳钥共〉拇唐咸?941(Vd)和感病的歐亞種美人指(Vv)為試材,通過室內(nèi)接菌,設(shè)置0 h、24 h、48 h 三個(gè)時(shí)間點(diǎn),研究病原菌誘導(dǎo)下感病和抗病種質(zhì)的基因表達(dá)差異,篩選出抗病基因,并利用qRT-PCR進(jìn)行驗(yàn)證?!窘Y(jié)果】根據(jù)|Log2Fold Change|≥2 和p≤0.05 的篩選條件,在0 h、24 h、48 h 抗病品種刺葡萄和感病品種美人指之間篩選到差異表達(dá)基因分別有5266、4725、4653 個(gè)。其中上調(diào)基因數(shù)分別為1737、1788、1727 個(gè),下調(diào)基因數(shù)為3529、2937、2926 個(gè);進(jìn)一步對(duì)這些差異基因進(jìn)行GO和KEGG富集分析,發(fā)現(xiàn)了具有防御反應(yīng)、跨膜運(yùn)輸和一些植物激素信號(hào)轉(zhuǎn)導(dǎo)等功能的基因或者通路。為了明確不同品種葡萄果實(shí)抗病之間的調(diào)控網(wǎng)絡(luò),進(jìn)行了WGCNA分析,得到7 個(gè)基因模塊,篩選出2 個(gè)與抗病高度相關(guān)的關(guān)鍵的模塊并且通過hub 基因分析明確模塊內(nèi)基因的互作網(wǎng)絡(luò)關(guān)系,鑒定出2 個(gè)類黃酮代謝相關(guān)基因、2 個(gè)植物激素信號(hào)轉(zhuǎn)導(dǎo)相關(guān)基因以及2 個(gè)hub 基因,并且通過熒光定量驗(yàn)證其表達(dá)情況與轉(zhuǎn)錄組數(shù)據(jù)一致?!窘Y(jié)論】獲得了刺葡萄與美人指響應(yīng)白腐病侵染的差異表達(dá)基因,并顯著富集在防御反應(yīng)、跨膜運(yùn)輸和一些植物激素信號(hào)轉(zhuǎn)導(dǎo)等通路,篩選出6 個(gè)抗病候選基因,包括F35H、EIN和MYB等。

    關(guān)鍵詞:葡萄;抗??;轉(zhuǎn)錄組;基因挖掘;WGCNA 分析

    中圖分類號(hào):S663.1 文獻(xiàn)標(biāo)志碼:A? 文章編號(hào):1009-9980(2023)04-0653-16

    葡萄在我國具有悠久的栽培歷史,因其具有豐富的營養(yǎng)物質(zhì)和較高的經(jīng)濟(jì)價(jià)值,深受廣大民眾的喜愛。近年來,隨著葡萄栽培面積和產(chǎn)量的持續(xù)增長(zhǎng),病害造成的損失也逐年增加,其中葡萄白腐病是危害最嚴(yán)重的病害之一。該病是由白腐墊殼孢[Conielladiplodiella (Speg.) Petrk & Sydow.]引起的真菌性病害[1],19 世紀(jì)在意大利首次報(bào)道,我國最早發(fā)現(xiàn)于1899 年,目前主要發(fā)生在我國東北、華北、西北和華東北部地區(qū),植株感病后,葉片邊緣出現(xiàn)水漬狀病斑,慢慢侵染到葉片中部,呈現(xiàn)出不太明顯的同心輪紋狀的褐色病斑;葡萄果粒逐漸變褐變軟,發(fā)病后期整個(gè)果穗腐爛,造成大面積減產(chǎn)[2-3]。目前在農(nóng)業(yè)生產(chǎn)上主要是化學(xué)藥劑防治葡萄白腐病,化學(xué)防治法雖然效果好、見效快,但是長(zhǎng)期用藥,不僅增加成本、降低果實(shí)品質(zhì),而且還會(huì)存在農(nóng)藥殘留,嚴(yán)重危害環(huán)境及食品安全。

    不同葡萄種質(zhì)對(duì)白腐病的抗性存在差異。歐亞種對(duì)白腐病的抗性較差,東亞種群的葡萄抗病性較強(qiáng)[4]。中國是葡萄屬植物重要的起源地之一,基本包括了所有的東亞種群,野生葡萄種質(zhì)資源豐富[5],利用野生種質(zhì)資源導(dǎo)入異質(zhì)抗性基因是解決栽培品種抗性基因貧乏的有效途徑。我國的刺葡萄對(duì)包括白腐病在內(nèi)的多種真菌性病害都具有很強(qiáng)的抗性[6-7],是葡萄白腐病抗性研究的良好材料。因此,解析葡萄抗病機(jī)制、挖掘抗性基因、培育具有較強(qiáng)抗病性的葡萄品種尤為重要。

    轉(zhuǎn)錄組測(cè)序已被廣泛用于解析擬南芥、番茄、葡萄和其他植物的抗病機(jī)制;研究結(jié)果能夠提供植物在病原菌侵染下應(yīng)答的信號(hào)通路和相關(guān)基因表達(dá)的信息。陳哲等[8]利用轉(zhuǎn)錄組測(cè)序分析,發(fā)現(xiàn)草莓受到炭疽病病原菌侵染后能夠影響類黃酮生物合成和植物激素信號(hào)傳導(dǎo)等通路中關(guān)鍵基因的表達(dá)水平,并鑒定出了19 個(gè)感病相關(guān)的差異表達(dá)基因。暹羅芽孢桿菌提高了貯藏杧果果實(shí)的抗病性,Jiang 等[9]利用暹羅芽孢桿菌的誘導(dǎo)機(jī)制,建立了杧果果實(shí)樣品在貯藏期間的比較轉(zhuǎn)錄組分析,部分基因(JAZ、BAK1 和PR1)被暹羅芽孢桿菌處理上調(diào),引發(fā)應(yīng)激反應(yīng),并誘導(dǎo)酚類抗性物質(zhì)的合成,從而提高杧果果實(shí)的抗病性。其次,一些基因(WRKY22、HSP90、CNGCs、SOD、PAL、4CL、CHS 和HCT)如過氧化物酶體、苯丙素、黃酮類和姜酚生物合成,在杧果果實(shí)中上調(diào),從而增強(qiáng)系統(tǒng)的抗病能力,刺激免疫反應(yīng)。姜長(zhǎng)岳[10]以毛葡萄丹鳳-2 葉片為材料接種白粉菌進(jìn)行轉(zhuǎn)錄組測(cè)序,從MYB家族中篩選出抗病相關(guān)轉(zhuǎn)錄因子VqMYB154,發(fā)現(xiàn)轉(zhuǎn)錄因子VqMYB154 具有正向調(diào)控3 個(gè)芪合酶靶基因表達(dá)以增強(qiáng)植物抗病性的功能。Jiao 等[11]通過轉(zhuǎn)錄組分析中國野生葡萄對(duì)白粉病的抗病機(jī)制,發(fā)現(xiàn)感病品種湖南一號(hào)下調(diào)基因比例明顯高于山-24,而上調(diào)基因主要與防御反應(yīng)和抗病代謝物的生物合成相關(guān),同時(shí)發(fā)現(xiàn)12 號(hào)和13 號(hào)染色體上的NLR基因?qū)Π追劬蠲舾?,其中許多基因在易感種質(zhì)的所有感染階段始終下調(diào)。Liu 等[12]利用轉(zhuǎn)錄組測(cè)序分析,從河岸葡萄中篩選了與抗病性相關(guān)的轉(zhuǎn)錄因子ERF、MYB、WRKY和bHLH,確定了197 個(gè)候選基因,可能參與霜霉病菌引起的防御反應(yīng),其中有5 個(gè)基因富含亮氨酸重復(fù)序列(LRR)。沈才琦等[13]通過對(duì)感炭疽病歐亞種里扎馬特和抗炭疽病刺葡萄黑珍珠的雜交后代中感炭疽病7-2-6 株系和抗炭疽病7-1-8 株系進(jìn)行轉(zhuǎn)錄組測(cè)序分析,發(fā)現(xiàn)葡萄響應(yīng)膠孢炭疽菌侵染的差異表達(dá)基因,顯著富集在植物與病原體相互作用、類黃酮生物合成和MAPK信號(hào)途徑等3 個(gè)通路,并篩選出8 個(gè)抗病候選基因,包括MLO蛋白及WRKY、ERF轉(zhuǎn)錄因子等。轉(zhuǎn)錄組分析已廣泛用于發(fā)掘植物響應(yīng)病病原菌脅迫的調(diào)控因子,這為研究刺葡萄抗病機(jī)制提供了理論依據(jù)。筆者在本研究中對(duì)白腐病病原菌侵染前后的感病歐亞種美人指和抗病中國野生刺葡萄的果實(shí)進(jìn)行轉(zhuǎn)錄組測(cè)序分析,通過轉(zhuǎn)錄組分析已經(jīng)鑒定了幾個(gè)可能在抗葡萄白腐病中起關(guān)鍵作用的基因。

    1 材料和方法

    1.1 試驗(yàn)材料

    以抗白腐病刺葡萄0941 和感白腐病歐亞種美人指轉(zhuǎn)色期的果實(shí)為試驗(yàn)材料,材料采自于國家葡萄種質(zhì)資源圃(鄭州)。在葡萄果實(shí)的果柄處接種直徑1 cm 的白腐病原菌菌塊,每組接種18 粒葡萄,保濕并置于28 ℃恒溫箱中培養(yǎng),同時(shí)在果實(shí)果柄處接種無菌水作為對(duì)照。分別在接菌0 h、24 h、48 h 收集去除感染部分后的果實(shí),每組均設(shè)置3 個(gè)生物學(xué)重復(fù)。采用多糖多酚植物總RNA提取試劑盒(華越洋,貨號(hào)0416-50),參照說明書提取樣品總RNA,置于-80 ℃超低溫冰箱保存?zhèn)溆谩?/p>

    1.2 文庫構(gòu)建及測(cè)序

    總RNA質(zhì)量檢測(cè)合格后,通過深圳華大基因科技服務(wù)有限公司進(jìn)行轉(zhuǎn)錄組測(cè)序。對(duì)得到的原始數(shù)據(jù)進(jìn)行質(zhì)控,去除低質(zhì)量的reads 獲得的高質(zhì)量的質(zhì)控?cái)?shù)據(jù),并儲(chǔ)存于格式為fastq 的數(shù)據(jù)文件中。以葡萄基因組12 ×(https://urgi.versailles.inra.fr/Species/Vitis)為參考基因組進(jìn)行比對(duì)。利用HISAT2 軟件將得到的clean reads 比對(duì)參考基因組進(jìn)行組裝。使用RSEM 軟件計(jì)算基因表達(dá)水平,然后轉(zhuǎn)換成FPKM(每千堿基轉(zhuǎn)錄序列的預(yù)期片段數(shù)/百萬堿基對(duì)測(cè)序)值。使用對(duì)比工具Bowtie2 將clean reads 與參考基因組進(jìn)行比對(duì)以統(tǒng)計(jì)基因比對(duì)率,比對(duì)率都在90%以上。利用組裝的轉(zhuǎn)錄本與GO 數(shù)據(jù)庫、KEGG數(shù)據(jù)庫、NR數(shù)據(jù)庫、Swiss-Prot 數(shù)據(jù)庫、Pfam數(shù)據(jù)庫、COG數(shù)據(jù)庫以及KOG數(shù)據(jù)庫進(jìn)行基因功能注釋。

    1.3 差異表達(dá)基因篩選和富集分析

    比對(duì)品種間基因的表達(dá)情況,篩選差異表達(dá)的部分進(jìn)一步富集分析。使用DESeq2 軟件進(jìn)行差異基因的篩選,將|log2FoldChange|≥2、p≤0.05 的基因作為顯著的差異表達(dá)基因,對(duì)篩選出的差異基因進(jìn)行GO功能富集和KEGG通路分析,以獲得差異表達(dá)基因的主要生物學(xué)功能及主要生化代謝途徑。

    1.4 WGCNA分析通過WGCNA分析進(jìn)一步探究差異表達(dá)基因之間的調(diào)控網(wǎng)絡(luò),旨在鑒定相關(guān)基因模塊及模塊內(nèi)的hub 基因。對(duì)每個(gè)基因模塊進(jìn)行評(píng)估,檢驗(yàn)每個(gè)模塊中基因的相關(guān)性,利用Cytoscape 3.8.0 計(jì)算每個(gè)基因的相關(guān)性,鑒定每個(gè)模塊內(nèi)的hub 基因。通過計(jì)算每個(gè)模塊的特征基因值,檢驗(yàn)各個(gè)基因白腐病之間的相關(guān)性。此外,為了進(jìn)一步探究模塊中基因的功能,利用R包(ClusterProfiler)進(jìn)行GO和KEGG富集分析。

    1.5 實(shí)時(shí)熒光定量PCR(qRT-PCR)分析

    為驗(yàn)證RNA-seq 測(cè)序檢測(cè)到DEGs,本研究基于分析挑取9 個(gè)差異表達(dá)基因進(jìn)行qRT-PCR 分析。使用Primer3plus 在線網(wǎng)站(https://www.primer3plus.com/index.html)設(shè)計(jì)引物,并由上海生工生物公司進(jìn)行合成,引物具體信息見表1。然后用Fasting-King RT Kit(With gDNA)反轉(zhuǎn)錄試劑盒(天根)將RNA 反轉(zhuǎn)錄為cDNA。每個(gè)樣本3 個(gè)重復(fù),使用FastStart Essential DNA Green Master 試劑盒按照其說明在ABI-7300 系統(tǒng)進(jìn)行qRT-PCR 分析,反應(yīng)程序:94 ℃ 30 s;94 ℃ 5 s,60 ℃ 30 s,40 個(gè)循環(huán)。以VvEF1-α 為內(nèi)參基因[14],用2- ΔΔCT 公式計(jì)算相對(duì)表達(dá)量。

    2 結(jié)果與分析

    2.1 轉(zhuǎn)錄組測(cè)序質(zhì)量分析

    轉(zhuǎn)錄組測(cè)序后對(duì)原始數(shù)據(jù)進(jìn)行過濾,一共得到798.77 百萬條高質(zhì)量的過濾序列,共產(chǎn)生119.82 Gb過濾堿基,每個(gè)樣本的Q30 均大于92%,清潔讀數(shù)率都在96%以上。轉(zhuǎn)錄組測(cè)序結(jié)果與葡萄參考基因組的比對(duì)率高,測(cè)序質(zhì)量好,可用于后續(xù)分析(表2)。對(duì)各組重復(fù)之間進(jìn)行相關(guān)性分析,每組樣品之間的相關(guān)系數(shù)均高于0.8(圖1),以確保試驗(yàn)結(jié)果的可靠性。

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

    篩選|log2FoldChange|≥2、p≤0.05 的基因作為顯著差異表達(dá)基因,在0 h、24 h、48 h 抗病品種刺葡萄(Vd)和感病品種美人指(Vv)之間的差異表達(dá)基因分別有5266、4725、4653 個(gè)(圖2-A)。其中上調(diào)基因數(shù)分別為1737、1788、1727 個(gè),下調(diào)基因數(shù)分別為3529、2937、2926 個(gè)。在抗病品種和感病品種之間的差異表達(dá)基因均是下調(diào)基因數(shù)目遠(yuǎn)大于上調(diào)基因數(shù)目,且24 h 和48 h 較0 h 下調(diào)基因數(shù)目大幅減少,3 個(gè)時(shí)期對(duì)比組共有差異表達(dá)基因有2684 個(gè)(圖2-B)。

    2.3 GO、KEGG功能注釋及富集分析

    將篩選出來的所有時(shí)期的差異表達(dá)基因進(jìn)行GO富集分析,結(jié)果分為3 類,分別為生物學(xué)過程、細(xì)胞組分和分子功能。結(jié)果僅展示出富集程度最為顯著的前20 條的代謝途徑(圖3)。在膜、細(xì)胞、細(xì)胞部分等組分中具有較高程度的富集,在生物學(xué)過程中防御反應(yīng)、跨膜運(yùn)輸、代謝過程、細(xì)胞過程等具有較高程度的富集,在分子功能中轉(zhuǎn)錄調(diào)節(jié)活性和催化活性等具有較高程度的富集,結(jié)果表明這些功能在刺葡萄、美人指與白腐病的防御反應(yīng)中發(fā)揮著重要作用。

    對(duì)刺葡萄和美人指之間的差異表達(dá)基因進(jìn)行KEGG代謝通路富集分析,取富集最為顯著的前20個(gè)通路進(jìn)行分析(圖4)。通過KEGG分析,表明上調(diào)差異基因主要富集在類黃酮生物合成、黃酮和黃酮醇生物合成、苯丙烷生物合成、芪類和二芳基庚烷類及姜酚生物合成、α-亞麻酸代謝、亞油酸代謝等代謝途徑中,而下調(diào)差異基因主要富集在植物與病原菌的相互作用、類黃酮生物合成、苯丙烷生物合成、植物激素信號(hào)和絲裂原活化蛋白激酶(MAPK)信號(hào)轉(zhuǎn)導(dǎo)途徑、ABC運(yùn)輸、苯丙氨酸代謝、氨基糖和核苷酸糖代謝、二萜生物合成、DNA復(fù)制等通路中。

    2.4 刺葡萄抗白腐病的相關(guān)差異表達(dá)基因

    刺葡萄和美人指接種白腐病病原菌0 h、24 h、48 h 后進(jìn)行轉(zhuǎn)錄組分析,從KEGG富集的通路中選擇白腐病病原菌誘導(dǎo)后表達(dá)差異較大的基因,作為重要的差異表達(dá)基因。差異表達(dá)基因主要來自次生代謝和植物激素信號(hào)轉(zhuǎn)導(dǎo)兩個(gè)通路。

    葡萄果實(shí)接種白腐病病原菌后,與苯丙烷、類黃酮、黃酮和黃酮醇生物合成等次生代謝途徑的相關(guān)基因在24 h 和48 h 均受到誘導(dǎo)而產(chǎn)生表達(dá)差異。圖5展示了0 h、24 h和48 h中苯丙烷、類黃酮、黃酮和黃酮醇生物合成途徑相關(guān)的差異表達(dá)基因(differentiallyexpressed genes,DEGs)及其轉(zhuǎn)錄水平。其中苯丙氨酸解氨酶PAL(Vitvi00g01431、Vitvi16g01503)、四香豆酸酯CoA 連接酶4CL(Vitvi16g00139)、查耳酮合成酶CHS(Vitvi16g00990、Vitvi16g01480)、莽草酸/奎寧酸羥基肉桂酰轉(zhuǎn)移酶HCT(Vitvi11g00735)、黃酮醇合酶FLS(Vitvi08g01575)這7個(gè)基因在受到白腐病原菌誘導(dǎo)48 h 后,在抗病品種刺葡萄中的表達(dá)量遠(yuǎn)高于感病品種美人指。芪合酶STS(Vitvi16g01465、Vitvi16g01466)的表達(dá)量在受到病原菌誘導(dǎo)后,在刺葡萄和美人指中均呈現(xiàn)上升趨勢(shì),且美人指中的表達(dá)量高于刺葡萄。

    白腐病病原菌侵染葡萄果實(shí)時(shí),植物激素信號(hào)轉(zhuǎn)導(dǎo)途徑中生長(zhǎng)素、細(xì)胞分裂素、赤霉素、脫落酸、油菜素甾醇、茉莉酸和水楊酸等重要防御相關(guān)的信號(hào)通路均發(fā)生富集,82 個(gè)差異表達(dá)基因顯著富集于植物激素信號(hào)轉(zhuǎn)導(dǎo)途徑(圖6),39 個(gè)DEGs(TCH4,4個(gè);AUXIAA,1 個(gè);CH3,2 個(gè);TF,1 個(gè);JAZ,2 個(gè);SAUR,16 個(gè);PR1,2 個(gè);MYC2,1 個(gè);ETR;2 個(gè);TIR,1 個(gè))在刺葡萄中的表達(dá)水平明顯高于美人指。其中,水楊酸信號(hào)途徑下游基因PR1(Vitvi03g00860)和茉莉酸信號(hào)途徑下游基因MYC2(Vitvi02g00231)在刺葡萄接菌處理后表達(dá)上調(diào),且PR1(Vitvi03g00860)在接菌48 h 表達(dá)水平很高,最高上調(diào)1000 倍(0.27<FPKM<1100),MYC2(Vitvi02g00231)在接菌0 h 和24 h 表達(dá)水平無明顯變化,但在48 h 上調(diào)10 倍(0<FPKM<13);茉莉酸信號(hào)途徑下游基因JAZ(Vitvi01g02293)在刺葡萄接菌24 h 下調(diào)表達(dá)3 倍,而后上調(diào)100 倍;乙烯信號(hào)途徑上游基因ETR(Vitvi05g00684)在刺葡萄和美人指接菌后都上調(diào)。

    2.5 接菌處理后轉(zhuǎn)錄因子的表達(dá)模式

    DEGs 中轉(zhuǎn)錄因子有344 個(gè),屬于39 個(gè)轉(zhuǎn)錄因子家族,主要包含40 個(gè)AP2/ERF、66 個(gè)MYB、22 個(gè)WRKY、35 個(gè)bHLH、24 個(gè)NAC、16 個(gè)LBD、16 個(gè)C2H2、11 個(gè)HD-ZIP 等(圖7)。在接菌處理后,其中10 個(gè)轉(zhuǎn)錄因子在美人指中的表達(dá)量明顯高于刺葡萄,例如AP2/ERF(Vitvi05g01073)在接菌24 h 美人指中的表達(dá)量是刺葡萄的60 倍,在接菌48 h 美人指的表達(dá)量是刺葡萄的35 倍。16 個(gè)轉(zhuǎn)錄因子在刺葡萄中的表達(dá)量顯著高于美人指,例如參與植物苯丙烷類次生代謝途徑的MYB(Vitvi02g01019)在接菌24 h 刺葡萄中的表達(dá)量是美人指的100 倍,在接菌48 h 刺葡萄的表達(dá)量是美人指的40 倍。8 個(gè)轉(zhuǎn)錄因子在受到白腐病原菌的誘導(dǎo)后的刺葡萄和美人指中表達(dá)模式相反,例如響應(yīng)激素信號(hào)轉(zhuǎn)導(dǎo)的AP2/ERF(Vitvi16g01438)在接菌24 h 刺葡萄中表達(dá)上調(diào)12倍,而在美人指中下調(diào)表達(dá)10 倍;在接菌48 h 刺葡萄中表達(dá)下調(diào)表達(dá)10 倍,而美人指上調(diào)3 倍。

    2.6 WGCNA分析鑒定刺葡萄抗病的關(guān)鍵基因

    經(jīng)過濾FPKM<1 的基因,4085 個(gè)差異基因用于構(gòu)建加權(quán)基因共表達(dá)網(wǎng)絡(luò)。通過計(jì)算18 個(gè)樣本表達(dá)水平的相關(guān)系數(shù)聚類分析,樣本間聚類良好,未出現(xiàn)離群樣本。根據(jù)差異基因的FPKM 值進(jìn)行相關(guān)度分析并聚類,相關(guān)度較高的基因被分配到同一個(gè)模塊中,圖中聚類樹的不同分支代表不同的基因共表達(dá)模塊,不同顏色代表不同的模塊,共劃分為7 個(gè)模塊,并繪制基因共表達(dá)網(wǎng)絡(luò)熱圖(圖8)。根據(jù)基因表達(dá)模塊的趨勢(shì)與抗病性狀的相關(guān)性分析,以皮爾遜相關(guān)系數(shù)(r>0.5)和顯著性p 值的大小(p<0.05)為篩選條件MEyellow、MEblack、MEgreen模塊與抗病性狀呈顯著正相關(guān),其中MEblack模塊在Vd0h、Vd24h、Vd48h 與抗病性狀呈顯著正相關(guān),MEyellow 模塊在Vd48h 與抗病性狀呈顯著正相關(guān)。對(duì)這兩個(gè)抗病相關(guān)模塊中的差異表達(dá)基因進(jìn)行KEGG富集分析,其富集的通路主要集中在類黃酮生物合成、MAPK信號(hào)途徑、苯丙烷生物合成、植物相互作用(圖9)。

    選取Meyellow、MEblack 模塊內(nèi)連通度最高的前150 個(gè)基因作為模塊的hub 基因,其中候選轉(zhuǎn)錄因子共有13 個(gè),并通過WGCNA 預(yù)測(cè)與它們可能存在互作關(guān)系的DEGs,利用Cytoscape 軟件進(jìn)行可視化分析(圖10)。在MEblcak 模塊中有4 個(gè)轉(zhuǎn)錄因子(MYB,2 個(gè);EIN 和TCP,各1 個(gè))與植物防御反應(yīng)有關(guān),其中與MYB4-Like(Vitvi17g00232)相關(guān)DEGs 的數(shù)量最多,共有104 個(gè);這些基因主要與植物的代謝反應(yīng)有關(guān),如苯丙氨酸生物合成、色氨酸代謝、糖酵解、泛素介導(dǎo)的蛋白水解等。在MEyellow 模塊中有9 個(gè)轉(zhuǎn)錄因子(MYB,3 個(gè);ERF,3 個(gè);bHLH,2 個(gè);bZIP,1 個(gè))與植物防御反應(yīng)有關(guān),其中MYB1R1(Vitvi04g00076)有相關(guān)性的DEGs 數(shù)量為134 個(gè),這些基因的功能注釋發(fā)現(xiàn)與多種植物激素有關(guān),如生長(zhǎng)素、乙烯、水楊酸和茉莉酸,還和氧化還原反應(yīng)、類黃酮生物合成等有關(guān)。

    通過KEGG富集分析、轉(zhuǎn)錄因子家族分析以及WGCNA中與抗病性相關(guān)聯(lián)的模塊,結(jié)合基因注釋,篩選出6 個(gè)與抗病相關(guān)的候選基因,分別是Vitvi01g02263(FAOMT)、Vitvi06g01183(F35H)、Vitvi06g01036(EIN3)、Vitvi04g00076(MYB1R1)、Vitvi17g00232(MYB4)和Vitvi07g02070(ERF098)。

    2.7 DEGs的RT-qPCR

    為了驗(yàn)證RNA-Seq 數(shù)據(jù)的可靠性,從DEGs 中選擇9 個(gè)基因,包括2 個(gè)類黃酮代謝相關(guān)的DEGs、2個(gè)植物激素信號(hào)轉(zhuǎn)導(dǎo)相關(guān)的DEGs、2 個(gè)共表達(dá)模塊篩選到的關(guān)鍵基因和3 個(gè)隨機(jī)挑選的DEGs,通過RT-qPCR方法分析了這些基因在刺葡萄和美人指受白腐病病原菌侵染后的表達(dá)量變化,熒光定量顯示以上基因在接種病菌前后的相對(duì)表達(dá)量與轉(zhuǎn)錄組數(shù)據(jù)一致,這說明轉(zhuǎn)錄組數(shù)據(jù)是可靠的(圖11)。其中FAOMT、F35H、EIN3、MYB4 的表達(dá)量在所有時(shí)期中刺葡萄均明顯高于美人指;MYB1R1、ERF098 在接菌48 h 后,在刺葡萄中的表達(dá)量較其他時(shí)期和美人指中的均有顯著提升。

    3 討論

    本研究中,利用RNA-Seq 技術(shù)對(duì)接種白腐病原菌后不同時(shí)間點(diǎn)的抗病和感病葡萄葉片進(jìn)行轉(zhuǎn)錄組測(cè)序,在0 h、24 h 和48 h 三個(gè)時(shí)期兩者之間分別篩得差異基因5266、4725和4653 個(gè)。當(dāng)白腐病病原菌侵染葡萄果實(shí)后,多個(gè)代謝通路的基因表達(dá)發(fā)生改變,其中次生代謝物可能積極參與葡萄抗病反應(yīng)。苯丙烷代謝物對(duì)植物的發(fā)育和生存至關(guān)重要,苯丙烷代謝物的生物合成和多樣性可通過加氧酶、氧化還原酶、連接酶和轉(zhuǎn)移酶等一系列酶來實(shí)現(xiàn),這些酶通過?;?、甲基化、糖基化和羥基化對(duì)代謝物的基本骨架進(jìn)行化學(xué)修飾[15]。苯丙氨酸由莽草酸途徑產(chǎn)生,苯丙類化合物是類黃酮、異黃酮和孜然芹素等多種酚類化合物的前體,研究發(fā)現(xiàn)苯丙氨酸在植物體內(nèi)積累,能顯著降低植物對(duì)病原菌的敏感性[16]。苯丙氨酸解氨酶(PAL)是苯丙烷代謝途徑中的關(guān)鍵酶,能夠?qū)⒈奖彼徂D(zhuǎn)為肉桂酸,是催化苯丙烷類次級(jí)代謝的第一步反應(yīng),起到樞紐的作用[17]。PAL參與花色苷、植物保護(hù)素、酚類代謝等物質(zhì)的合成,該基因在大多數(shù)植物中是以多基因家族存在,PAL基因家族成員調(diào)控的代謝途徑也各有差異。有研究表明,植物中PAL活性降低時(shí),抵抗外界干擾的能力下降[18]。肉桂酸4-羥基化酶(C4H)為通用的苯丙氨酸途徑限速酶,參與苯丙烷代謝途徑第二步催化反應(yīng),起到承上啟下的作用。該酶能提高自身的活性,來抵抗外界的環(huán)境脅迫。C4H將反式肉桂酸催化為β-香豆酸,通過4-香豆酰輔酶A連接酶轉(zhuǎn)化為香豆酞輔酶A,為下游類黃酮、木質(zhì)素等次級(jí)代謝途徑提供前體物質(zhì)[19]。PAL催化苯丙烷代謝合成途徑中的第一步反應(yīng)產(chǎn)生的物質(zhì)與C4H催化生成4-香豆酸。PAL 與C4H這兩個(gè)酶通常情況下均以協(xié)同的方式進(jìn)行表達(dá),來控制黃酮物質(zhì)的合成,但通常情況下,C4H的活性比PAL的活性低很多[20]。在本試驗(yàn)中PAL(Vitvi16g01503)、4CL(Vitvi16g00139)基因在受到病原菌誘導(dǎo)48 h 后,其在抗病品種刺葡萄中的表達(dá)量遠(yuǎn)高于美人指。在植物體內(nèi)PAL 作為苯丙烷類代謝定速酶,與植物的抗病性有重要關(guān)系。已報(bào)道在金針菇感染異型葡枝霉菌之后PAL活性明顯增強(qiáng)[21],β-1,3-葡寡糖10 倍誘抗劑處理后的葡萄葉片苯丙氨酸解氨酶(PAL)活性增強(qiáng),說明處理之后提高了體內(nèi)防御酶和抗病相關(guān)蛋白活性,進(jìn)而誘導(dǎo)植株產(chǎn)生系統(tǒng)抗病性[22],以上結(jié)果表明PAL與植物的抗病性呈現(xiàn)正相關(guān)關(guān)系。

    類黃酮代謝途徑是苯丙烷代謝途徑通過查爾酮合成酶(CHS)延伸的另一個(gè)重要分支,類黃酮代謝途徑廣泛存在于各種植物體內(nèi),類黃酮化合物是類黃酮代謝途徑中產(chǎn)生的大量次生代謝物質(zhì),大自然中類黃酮化合物一般分為六大類:黃酮、黃酮醇、黃烷酮、黃烷醇、花色素、異黃酮[23]。在植物中,葉片和花等器官的類黃酮化合物含量比較高,類黃酮化合物作為次生代謝物質(zhì)在植物中可以作為植物和微生物之間的關(guān)鍵信號(hào)物質(zhì)發(fā)揮植物防御素的功能,在植物中有很強(qiáng)的抗病抑菌作用。CHS 可以通過催化形成不同的類黃酮化合物來提高植物的抗病能力,也可通過改變植保素、黃酮素、抗氧化酶SOD和POD 等的積累以及防御酶的活性來提高防御能力[24-27]。研究表明,在番茄和蘋果中過表達(dá)CHS 能夠增強(qiáng)對(duì)病原菌等環(huán)境脅迫的抗性[28-29]。抗枯萎病的棉花品種中CHS表達(dá)量明顯上調(diào)后,類黃酮的含量提高從而增強(qiáng)其對(duì)枯萎病的抗性[30]。此外,郭澤西[31]已經(jīng)證明抗性品種中查爾酮合成酶基因CHSI的表達(dá)量比感病品種高,過表達(dá)CHSl能夠提高感病葡萄品種對(duì)灰霉病和霜霉病的抗性。在本研究中發(fā)現(xiàn)接種白腐病病原菌誘導(dǎo)48 h 后,在抗病品種刺葡萄中的表達(dá)量遠(yuǎn)高于感病品種美人指,這與前人的研究一致。同時(shí)筆者通過加權(quán)共表達(dá)網(wǎng)絡(luò)分析發(fā)現(xiàn)的兩個(gè)候選基因:F35H和FAOMT,均是類黃酮代謝途徑中的關(guān)鍵基因,且FAOMT是F35H的下游基因[32],二者在刺葡萄中的表達(dá)量均明顯高于美人指。類黃酮代謝通路由結(jié)構(gòu)基因和調(diào)節(jié)基因共同控制,調(diào)節(jié)基因包括MYB、bHLH 等轉(zhuǎn)錄因子[33],此前發(fā)現(xiàn)MYB4 能夠控制肉桂酸4-羥化酶(C4H)的負(fù)表達(dá)。筆者發(fā)現(xiàn)類黃酮代謝途徑中的關(guān)鍵基因和轉(zhuǎn)錄因子在抗病品種刺葡萄和感病品種美人指中表現(xiàn)出明顯的差異,因此,在刺葡萄響應(yīng)白腐病病原菌侵染時(shí),類黃酮物質(zhì)可能起到重要作用。

    此外,EIN3 和ERF098 是模塊中篩選出刺葡萄響應(yīng)白腐病侵染的關(guān)鍵基因,轉(zhuǎn)錄因子EIN3 是乙烯信號(hào)轉(zhuǎn)導(dǎo)途徑中的調(diào)控因子,可激活乙烯反應(yīng)因子(ERF)[34]。本研究中乙烯信號(hào)途徑下游基因EIN3的表達(dá)水平在處理不同時(shí)間點(diǎn)中刺葡萄均高于美人指,而EIN3 下游基因ERF098 在48 h 抗病品種刺葡萄中的表達(dá)量明顯高于其他時(shí)間和感病品種美人指,而在刺葡萄中48 h 的表達(dá)量明顯上升這可能與反應(yīng)時(shí)間有關(guān)。由此可知,乙烯可能在刺葡萄抗白腐病反應(yīng)的激素信號(hào)轉(zhuǎn)導(dǎo)中發(fā)揮了重要作用。

    參考文獻(xiàn)References:

    [1] AKKURT M,WELTER L,MAUL E,TOPFER R,ZYPRIAN E.Development of scar markers linked to powdery mildew (Uncinulanecator) resistance in grapevine (Vitis vinifera L.and Vitissp.)[J]. Molecular Breeding,2007,19(2):103-111.

    [2] ALI K,MALTESE F,ZYPRIAN E,REX M,CHOI Y H,VERPOORTER. NMR metabolic fingerprinting based identificationof grapevine metabolites associated with downy mildew resistance[J]. Journal of Agricultural and Food Chemistry,2009,57(20):9599-9606.

    [3] 王躍進(jìn),賀普超. 中國葡萄屬野生種抗白腐病機(jī)制研究[J]. 北方園藝,1988(1):5-8.

    WANG Yuejin,HE Puchao. Study of white rot disease resistantmechanism in ChineseVitis wild species[J]. Northern Horticulture,1988(1):5-8.

    [4] 劉崇懷. 中國葡萄屬(Vitis L.)植物分類與地理分布研究[D].鄭州:河南農(nóng)業(yè)大學(xué),2012.

    LIU Chonghuai. Studies on taxonomy and geographical distributionof Chinese wild grape species[D]. Zhengzhou:Henan AgriculturalUniversity,2012.

    [5] 董陽輝,徐佩娟,王藝平,錢劍銳,何鐵海,朱朝磊,陳寶忠,王洪興,張建萍. 葡萄白腐病研究[J]. 江西農(nóng)業(yè)學(xué)報(bào),2011,23(2):107-110.

    DONG Yanghui,XU Peijuan,WANG Yiping,QIAN Jianrui,HETiehai,ZHU Chaolei,CHEN Baozhong,WANG Hongxing,ZHANG Jianping. Study on white rot of grape[J]. Acta AgriculturaeJiangxi,2011,23(2):107-110.

    [6] 張穎,樊秀彩,孫海生,姜建福,劉崇懷. 葡萄種質(zhì)資源對(duì)葉片白腐病的抗性鑒定及評(píng)價(jià)[J]. 果樹學(xué)報(bào),2017,34(9):1095-1105.

    ZHANG Ying,F(xiàn)AN Xiucai,SUN Haisheng,JIANG Jianfu,LIUChonghuai. Identification and evaluation of resistance to whiterot in grape resources[J]. Journal of Fruit Science,2017,34(9):1095-1105.

    [7] 賈云飛. 刺葡萄0943VdWRKY70 轉(zhuǎn)錄因子抗病功能分析[D].洛陽:河南科技大學(xué),2018.

    JIA Yunfei. Function analysis ofVdWRKY70 transcription factordisease resistance fromVitis davidii 0943[D]. Luoyang:HenanUniversity of Science and Technology,2018.

    [8] 陳哲,黃靜,趙佳,梁宏. 草莓應(yīng)答炭疽菌侵染的轉(zhuǎn)錄組分析[J].植物保護(hù),2020,46(3):138-146.

    CHEN Zhe,HUANG Jing,ZHAO Jia,LIANG Hong. Transcriptomicsanalysis of strawberry response toColletotrichum theobromicolainfection[J]. Plant Protection,2020,46(3),138-146.

    [9] JIANG Z C,LI R,TANG Y,CHENG Z Y,QIAN M J,LI W,SHAO Y Z. Transcriptome analysis reveals the inducing effectofBacillus siamensis on disease resistance in postharvest mango[J]. Foods,2022,11(1):107.

    [10] 姜長(zhǎng)岳. 中國野生毛葡萄VqMYB154 轉(zhuǎn)錄因子調(diào)控抗白粉病機(jī)理研究[D]. 楊凌:西北農(nóng)林科技大學(xué),2021.

    JIANG Changyue. Mechanism analysis of transcription factorVqMYB154 regulating the resistance to powdery mildew in ChinesewildVitis quinquangularis[D]. Yangling:Northwest A&FUniversity,2021.

    [11] JIAO C,SUN X P,YAN X X,XU X Z,YAN Q,GAO M,F(xiàn)EI ZJ,WANG X P. Comparative transcriptome profiling of Chinesewild grapes provides insights into powdery mildew resistance[J].Phytopathology,2021,111(11):2041-2051.

    [12] LIU L,ZHANG B,WANG H,YU S Y,GUAN T S,HUANG YF,LIU C Y. Candidate resistance genes selection and transcriptomeanalysis for the early responses toPlasmopara viticola infectionin grape cultivars[J]. Journal of Plant Pathology,2020,102(3):857-869.

    [13] 沈才琦,孫磊,張穎,姜建福,劉崇懷,樊秀彩. 膠孢炭疽菌侵染不同抗性葡萄葉片的轉(zhuǎn)錄組學(xué)分析[J]. 果樹學(xué)報(bào),2022,39(5):730-742.

    SHEN Caiqi,SUN Lei,ZHANGYing,JIANG Jianfu,LIU Chonghuai,F(xiàn)AN Xiucai. Transcriptomic analysis of different resistantleaves infected byColletotrichum gloeosporioides in grape[J].Journal of Fruit Science,2022,39(5):730-742.

    [14] 查倩,奚曉軍,蔣愛麗,田益華,王世平. 葡萄實(shí)時(shí)定量PCR 中穩(wěn)定內(nèi)參基因的篩選[J]. 果樹學(xué)報(bào),2016,33(3):268-274.

    ZHAQian,XI Xiaojun,JIANG Aili,TIANYihua,WANG Shiping.Identification of the appropriate reference gene through using areal- time quantitative PCR in grapes[J]. Journal of Fruit Science,2016,33(3):268-274.

    [15] WANG S C,ALSEEKH S,F(xiàn)ERNIE A R,LUO J. The structureand function of major plant metabolite modifications[J]. MolecularPlant,2019,12(7):899-919.

    [16] KUMAR V,HATAN E,BAR E,DAVIDOVICH-RIKANATI R,DORON- FAIGENBOIM A,SPITZER- RIMON B,ELAD Y,ALKAN N,LEWINSOHN E,OREN-SHAMIR M. Phenylalanineincreases chrysanthemum flower immunity againstBotrytiscinerea attack[J]. The Plant Journal,2020,104(1):226-240.

    [17] 梁冬. 大豆PAL2-3 基因的功能研究[D]. 哈爾濱:東北農(nóng)業(yè)大學(xué),2017.

    LIANG Dong. Functional study ofPAL2-3 gene in soybean[D].Harbin:Northeast Agricultural University,2017.

    [18] 董艷珍. 植物苯丙氨酸解氨酶基因的研究進(jìn)展[J]. 生物技術(shù)通報(bào),2006(S1):31-33.

    DONG Yanzhen. Research advances on plant phenylalanine ammonia-lyase gene[J]. Biotechnology Bulletin,2006(S1):31-33.

    [19] 鄒廣平,李雅麗,馮夢(mèng)薇,楊修,許智偉. 甘草懸浮細(xì)胞肉桂酸-4-羥基化酶(C4H)基因的克隆及表達(dá)分析[J]. 生物技術(shù)通報(bào),2017,33(11):96-100.

    ZOU Guangping,LI Yali,F(xiàn)ENG Mengwei,YANG Xiu,XU Zhiwei.Cloning and expression analysis of cinnamate 4-hydroxylasegene fromGlycyrrhiza uralensis suspension cultured cells[J].Biotechnology Bulletin,2017,33(11):96-100.

    [20] 陳蒙. 山葡萄苯丙烷代謝相關(guān)基因的克隆表達(dá)與功能驗(yàn)證[D].延吉:延邊大學(xué),2019.

    CHEN Meng. Cloning and expression and functional verificationof phenypropane pathway gene inVitis amurensis[D]. Yanji:Yanbian University,2019.

    [21] 黃靜,李長(zhǎng)田,賈傳文,王曉敏,呂瑞娜. 異型葡枝霉感染不同金針菇后防御酶系的活性變化[J]. 中國食用菌,2021,40(1),103-108.

    HUANG Jing,LI Changtian,JIA Chuanwen,WANG Xiaomin,L? Ruina. Changes of defensive enzyme activities in differentFlammulina filiformis varieties after inoculation withCladobotryumvarium[J]. Edible Fungi of China,2021,40(1):103-108.

    [22] 王曉琳,吳琴燕,彭燕瓊,鄔劼,黃潔雪,毛妮妮,吉沐祥. 2 種生物誘抗劑對(duì)葡萄霜霉病的誘導(dǎo)抗病作用[J]. 中國農(nóng)學(xué)通報(bào),2021,37(32):127-131.

    WANG Xiaolin,WU Qinyan,PENG Yanqiong,WU Jie,HUANG Jiexue,MAO Nini,JI Muxiang. Effects of two inducerson induced resistance to grapevine downy mildew[J]. ChineseAgricultural Science Bulletin,2021,37(32):127-131.

    [23] TAYLOR L P,GROTEWOLD E. Flavonoids as developmentalregulators[J]. Current Opinion in Plant Biology,2005,8(3):317-323.

    [24] QIAO Y,MENG X J,JIN X X,DING G H. Genes expression ofkey enzymes in phenylpropanes metabolism pathway in cucumberwith RT-PCR[J]. Advanced Materials Research,2013,746:53-57.[25] CHEN L J,GUO H M,LIN Y,CHENG H M. Chalcone synthaseEaCHS1 fromEupatorium adenophorum functions in saltstress tolerance in tobacco[J]. Plant Cell Reports,2015,34(5):885-894.

    [26] 陳帥. 煙草類黃酮代謝途徑中關(guān)鍵酶CHS 基因與R2R3MYB 類轉(zhuǎn)錄抑制因子功能研究[D]. 雅安:四川農(nóng)業(yè)大學(xué),2017.

    CHEN Shuai. Function analysis ofCHS genes and R2R3 MYBrepressors related to flavonoid biosynthesis pathway inNicotianatabacum[D]. Yaan:Sichuan Agricultural University,2017.

    [27] RAHMAN F U,KHAN I A,ASLAM A,LIU R T,SUN L,WUY D,ASLAM M M,KHAN A U,LI P,JIANG J F,F(xiàn)AN X C,LIU C H,ZHANG Y. Transcriptome analysis reveals pathogenesis-related gene 1 pathway against salicylic acid treatment ingrapevine (Vitis vinifera L.)[J]. Front Genet,2022,13,1033288.

    [28] KROL P,IGIELSKI R,POLLMANN S,KEPCZYNSKA E.Priming of seeds with methyl jasmonate induced resistance tohemi-biotrophFusarium oxysporum f. sp.lycopersici in tomatovia 12-oxo-phytodienoic acid,salicylic acid,and flavonol accumulation[J]. Journal of Plant Physiology,2015,179:122-132.

    [29] JIAO C,SUN X P,YAN X X,XU X Z,YAN Q,GAO M,F(xiàn)EI ZJ,WANG X P. Grape transcriptome response to powdery mildewinfection:comparative transcriptome profiling of Chinese wildgrapes provides insights into powdery mildew resistance[J]. Phytopathology,2022,111(11):2041-2051.

    [30] 黃啟秀. 海島棉類黃酮代謝途徑中與枯萎病抗性相關(guān)基因的克隆及功能驗(yàn)證[D]. 烏魯木齊:新疆農(nóng)業(yè)大學(xué),2017.

    HUANG Qixiu. Cloning and functional verification of genes relatedtoFusarium wilt resistance in flavonoid mtabolism pathwayof island cotton (Gossypium barbadense L.) [D]. Urumqi:Xinjiang Agricultural University,2017.

    [31] 郭澤西,孫大運(yùn),曲俊杰,潘鳳英,劉露露,尹玲. 查爾酮合成酶基因在葡萄抗灰霉病和霜霉病中的作用[J]. 中國農(nóng)業(yè)科學(xué),2022,55(6):1139-1148.

    GUO Zexi,SUN Dayun,QU Junjie,PAN Fengying,LIU Lulu,YIN Ling. The role of chalcone synthase gene in grape resistanceto gray mold and downy mildew[J]. Scientia AgriculturaSinica,2022,55(6):1139-1148.

    [32] LI P Q,RUAN Z,F(xiàn)EI Z X,YAN J J,TANG G H. Integratedtranscriptome and metabolome analysis revealed that flavonoidbiosynthesis may dominate the resistance ofZanthoxylum bungeanumagainst stem canker[J]. Journal of Agricultural andFood Chemistry,2021,69(22):6360-6378.

    [33] SHVACHKO N A,LOSKUTOV I G,SEMILET T V,POPOV VS,KOVALEVA O N,KONAREV A V. Bioactive components inoat and barley grain as a promising breeding trend for functionalfood production[J]. Molecules,2021,26(8):2260.

    [34] GU S Y,WANG L C,CHEUH C M,LO W S.CHITINASELIKE1 regulates root development of dark- grown seedlings bymodulating ethylene biosynthesis inArabidopsis thaliana[J].Frontiers in Plant Science,2019,10:600.

    猜你喜歡
    轉(zhuǎn)錄組抗病葡萄
    我國小麥基因組編輯抗病育種取得突破
    葡萄熟了
    當(dāng)葡萄成熟時(shí)
    女報(bào)(2020年10期)2020-11-23 01:42:42
    bZIP轉(zhuǎn)錄因子在植物激素介導(dǎo)的抗病抗逆途徑中的作用
    基于轉(zhuǎn)錄組測(cè)序的山茱萸次生代謝生物合成相關(guān)基因的挖掘
    金釵石斛轉(zhuǎn)錄組SSR位點(diǎn)信息分析
    人參屬藥用植物轉(zhuǎn)錄組研究進(jìn)展
    葡萄新品種 優(yōu)質(zhì)又抗病
    大豆轉(zhuǎn)錄組測(cè)序研究進(jìn)展綜述
    精品欧美国产一区二区三| 精品不卡国产一区二区三区| 99久国产av精品国产电影| 亚洲av电影在线观看一区二区三区 | 日本黄色片子视频| 成人性生交大片免费视频hd| 99久久精品国产国产毛片| 春色校园在线视频观看| 美女内射精品一级片tv| 国产综合精华液| 女人十人毛片免费观看3o分钟| 深爱激情五月婷婷| 少妇猛男粗大的猛烈进出视频 | 亚洲av一区综合| 69av精品久久久久久| 18禁裸乳无遮挡免费网站照片| 2018国产大陆天天弄谢| 色哟哟·www| 国产精品一区二区三区四区久久| 精品午夜福利在线看| 少妇的逼好多水| 嫩草影院精品99| 日韩成人av中文字幕在线观看| 中国国产av一级| 免费观看无遮挡的男女| 久久国内精品自在自线图片| 国产精品一区www在线观看| 麻豆精品久久久久久蜜桃| 久久草成人影院| 99热这里只有精品一区| 精品亚洲乱码少妇综合久久| 精品久久久噜噜| 少妇的逼水好多| 尾随美女入室| 国产精品一区二区在线观看99 | 草草在线视频免费看| 草草在线视频免费看| 亚洲av成人av| 亚洲最大成人av| 久久6这里有精品| 亚洲欧美清纯卡通| 成人av在线播放网站| 国产亚洲av嫩草精品影院| 赤兔流量卡办理| 韩国av在线不卡| 一个人观看的视频www高清免费观看| 床上黄色一级片| 男的添女的下面高潮视频| 免费av观看视频| 国产乱人偷精品视频| 国产精品嫩草影院av在线观看| 亚洲av电影在线观看一区二区三区 | 尤物成人国产欧美一区二区三区| 熟妇人妻不卡中文字幕| 成人亚洲欧美一区二区av| 久久这里只有精品中国| 嫩草影院入口| 国产亚洲av片在线观看秒播厂 | 亚洲精品国产成人久久av| 春色校园在线视频观看| 黄色欧美视频在线观看| 精品人妻偷拍中文字幕| 一级二级三级毛片免费看| 99热6这里只有精品| 青春草视频在线免费观看| 肉色欧美久久久久久久蜜桃 | 免费大片18禁| 久久久久久久大尺度免费视频| 亚洲人成网站在线观看播放| 欧美日韩在线观看h| 成人漫画全彩无遮挡| 99热6这里只有精品| 午夜免费观看性视频| 麻豆久久精品国产亚洲av| 亚洲国产精品专区欧美| 午夜精品在线福利| 成人午夜精彩视频在线观看| 亚洲乱码一区二区免费版| a级毛片免费高清观看在线播放| 国内精品宾馆在线| 美女被艹到高潮喷水动态| 国产成人精品久久久久久| 97超视频在线观看视频| 人妻少妇偷人精品九色| 亚洲无线观看免费| 久久这里只有精品中国| 亚洲va在线va天堂va国产| 精品亚洲乱码少妇综合久久| 成人毛片60女人毛片免费| 国产男人的电影天堂91| 在线 av 中文字幕| 18+在线观看网站| 亚洲精品,欧美精品| 青青草视频在线视频观看| 成人二区视频| av福利片在线观看| 美女高潮的动态| 看黄色毛片网站| 国产精品国产三级专区第一集| 精品久久久久久久久av| 超碰97精品在线观看| 内射极品少妇av片p| 国产 一区精品| 一夜夜www| 免费看日本二区| 尾随美女入室| 精品国产三级普通话版| 大香蕉97超碰在线| kizo精华| 婷婷六月久久综合丁香| 午夜福利视频1000在线观看| 免费观看精品视频网站| 国产午夜精品久久久久久一区二区三区| 噜噜噜噜噜久久久久久91| 韩国av在线不卡| 欧美人与善性xxx| 少妇人妻一区二区三区视频| 视频中文字幕在线观看| 精品亚洲乱码少妇综合久久| 在线观看av片永久免费下载| 日本爱情动作片www.在线观看| 青春草国产在线视频| 久久精品国产亚洲av涩爱| 精品久久久久久久久亚洲| 色播亚洲综合网| 中文字幕免费在线视频6| 欧美最新免费一区二区三区| 国产精品三级大全| 男人狂女人下面高潮的视频| 欧美日韩综合久久久久久| 亚洲av免费在线观看| 在线 av 中文字幕| 伊人久久精品亚洲午夜| 爱豆传媒免费全集在线观看| 麻豆av噜噜一区二区三区| 亚洲欧美日韩东京热| 精华霜和精华液先用哪个| 超碰av人人做人人爽久久| 中文字幕免费在线视频6| 欧美区成人在线视频| xxx大片免费视频| 人妻少妇偷人精品九色| 少妇高潮的动态图| 久久久国产一区二区| 亚洲精品成人av观看孕妇| 伊人久久精品亚洲午夜| 国产精品人妻久久久影院| 91精品伊人久久大香线蕉| 久久久精品欧美日韩精品| 乱码一卡2卡4卡精品| 国产淫片久久久久久久久| 精品久久国产蜜桃| 国产精品无大码| 日韩精品青青久久久久久| 肉色欧美久久久久久久蜜桃 | 小蜜桃在线观看免费完整版高清| 国产 一区 欧美 日韩| 成人一区二区视频在线观看| 尾随美女入室| h日本视频在线播放| 九九在线视频观看精品| 久久久久久九九精品二区国产| 亚洲乱码一区二区免费版| 精品国产三级普通话版| 有码 亚洲区| 欧美一区二区亚洲| 久久精品夜色国产| 黄色欧美视频在线观看| 精品不卡国产一区二区三区| av国产免费在线观看| av又黄又爽大尺度在线免费看| 精品久久久久久久末码| 亚洲最大成人手机在线| 男人舔奶头视频| 麻豆久久精品国产亚洲av| 亚洲电影在线观看av| 亚洲第一区二区三区不卡| 日韩中字成人| 欧美一区二区亚洲| 最近2019中文字幕mv第一页| 午夜精品一区二区三区免费看| 黄色一级大片看看| 亚洲真实伦在线观看| 久久99热这里只有精品18| 欧美日韩视频高清一区二区三区二| 亚洲四区av| 免费黄频网站在线观看国产| 一个人看的www免费观看视频| 国产三级在线视频| 亚洲国产欧美人成| 亚洲av免费高清在线观看| 国产亚洲av片在线观看秒播厂 | 91久久精品电影网| 免费不卡的大黄色大毛片视频在线观看 | 国产亚洲91精品色在线| 亚洲av一区综合| 亚洲欧美一区二区三区国产| 亚洲人成网站高清观看| www.av在线官网国产| 国产亚洲5aaaaa淫片| 久久久精品94久久精品| 简卡轻食公司| 夜夜爽夜夜爽视频| av.在线天堂| 高清日韩中文字幕在线| 亚洲综合精品二区| 日本三级黄在线观看| 18+在线观看网站| 亚洲激情五月婷婷啪啪| 国产精品国产三级国产av玫瑰| 国产男人的电影天堂91| 卡戴珊不雅视频在线播放| 亚洲欧美中文字幕日韩二区| 狠狠精品人妻久久久久久综合| a级一级毛片免费在线观看| 久久久精品欧美日韩精品| 一级毛片我不卡| 亚洲天堂国产精品一区在线| 黄色一级大片看看| 国产不卡一卡二| 亚洲aⅴ乱码一区二区在线播放| 亚洲一区高清亚洲精品| 国产精品99久久久久久久久| 哪个播放器可以免费观看大片| 我要看日韩黄色一级片| 国产黄色视频一区二区在线观看| 亚洲美女搞黄在线观看| 国产伦一二天堂av在线观看| 又黄又爽又刺激的免费视频.| 80岁老熟妇乱子伦牲交| 精品国产露脸久久av麻豆 | 亚洲av一区综合| 免费看av在线观看网站| 不卡视频在线观看欧美| 婷婷色综合大香蕉| 两个人视频免费观看高清| 三级国产精品片| 国产探花极品一区二区| 在线观看一区二区三区| 日日干狠狠操夜夜爽| 久热久热在线精品观看| 黄色一级大片看看| 亚洲av成人精品一二三区| 亚洲va在线va天堂va国产| 国产成人精品一,二区| 寂寞人妻少妇视频99o| 日韩精品有码人妻一区| 国产欧美另类精品又又久久亚洲欧美| 亚洲内射少妇av| 18禁裸乳无遮挡免费网站照片| 又爽又黄a免费视频| 免费观看av网站的网址| 黑人高潮一二区| 亚洲精品日本国产第一区| 国产探花在线观看一区二区| 三级国产精品欧美在线观看| 国产男人的电影天堂91| 国产精品美女特级片免费视频播放器| 中文字幕久久专区| 国产有黄有色有爽视频| 在线免费观看的www视频| 在线 av 中文字幕| 青春草视频在线免费观看| 日韩av在线免费看完整版不卡| 综合色av麻豆| 午夜激情福利司机影院| 赤兔流量卡办理| 亚洲国产高清在线一区二区三| 日本色播在线视频| 日日干狠狠操夜夜爽| 国产v大片淫在线免费观看| 日本熟妇午夜| 精品人妻视频免费看| 午夜福利视频精品| 国产大屁股一区二区在线视频| 99视频精品全部免费 在线| 男女啪啪激烈高潮av片| 联通29元200g的流量卡| 黄色一级大片看看| 3wmmmm亚洲av在线观看| 搡老乐熟女国产| 色视频www国产| 亚洲在久久综合| 2022亚洲国产成人精品| 嫩草影院入口| 亚洲av中文字字幕乱码综合| 一区二区三区高清视频在线| 嫩草影院精品99| 国产麻豆成人av免费视频| 十八禁国产超污无遮挡网站| 天天躁夜夜躁狠狠久久av| 在线观看美女被高潮喷水网站| 麻豆乱淫一区二区| 久久久久久久久久成人| 国产真实伦视频高清在线观看| 日韩av在线大香蕉| 老女人水多毛片| 自拍偷自拍亚洲精品老妇| 一个人观看的视频www高清免费观看| 国产亚洲av片在线观看秒播厂 | 亚洲精品久久久久久婷婷小说| 亚洲国产欧美在线一区| 免费看日本二区| 中国美白少妇内射xxxbb| 岛国毛片在线播放| 亚洲一级一片aⅴ在线观看| 日韩亚洲欧美综合| 成人美女网站在线观看视频| 国产真实伦视频高清在线观看| 久久综合国产亚洲精品| 久久精品国产鲁丝片午夜精品| 寂寞人妻少妇视频99o| 搡老妇女老女人老熟妇| 中文字幕人妻熟人妻熟丝袜美| 在线观看av片永久免费下载| 男女国产视频网站| 蜜桃亚洲精品一区二区三区| 人人妻人人澡人人爽人人夜夜 | 欧美成人一区二区免费高清观看| 精品亚洲乱码少妇综合久久| 大话2 男鬼变身卡| 久久久久久国产a免费观看| 国内精品宾馆在线| 中文字幕亚洲精品专区| 91精品一卡2卡3卡4卡| 丰满人妻一区二区三区视频av| 久99久视频精品免费| 在线播放无遮挡| 大香蕉久久网| 亚洲av福利一区| 两个人视频免费观看高清| 国产69精品久久久久777片| av.在线天堂| 日日摸夜夜添夜夜爱| 国产淫语在线视频| eeuss影院久久| 搡老妇女老女人老熟妇| 国产亚洲一区二区精品| 国产男人的电影天堂91| 毛片一级片免费看久久久久| 久久这里有精品视频免费| 国产伦精品一区二区三区视频9| 国产黄色小视频在线观看| 高清毛片免费看| 大香蕉久久网| 亚洲国产欧美人成| 又爽又黄无遮挡网站| 三级国产精品欧美在线观看| 国产高清有码在线观看视频| 欧美一级a爱片免费观看看| 亚洲激情五月婷婷啪啪| 国产高清有码在线观看视频| 国产淫语在线视频| 日韩伦理黄色片| 亚洲欧美中文字幕日韩二区| 国产精品一区二区在线观看99 | 春色校园在线视频观看| 亚洲欧美日韩东京热| 丰满人妻一区二区三区视频av| 国产在视频线精品| 99re6热这里在线精品视频| 亚洲自偷自拍三级| 久久久久久久大尺度免费视频| 天天一区二区日本电影三级| 亚洲精品日韩av片在线观看| 毛片一级片免费看久久久久| 国产成人精品一,二区| av国产免费在线观看| 国产亚洲av嫩草精品影院| 尾随美女入室| 精品午夜福利在线看| 亚洲av电影不卡..在线观看| 国产精品久久久久久久电影| 视频中文字幕在线观看| 精品国内亚洲2022精品成人| 亚洲成人中文字幕在线播放| 日韩亚洲欧美综合| 大陆偷拍与自拍| 啦啦啦中文免费视频观看日本| 99久久人妻综合| 午夜福利在线观看免费完整高清在| 一级毛片 在线播放| 国产单亲对白刺激| 日韩精品青青久久久久久| 欧美日韩国产mv在线观看视频 | 国产视频首页在线观看| 亚洲精品一区蜜桃| 哪个播放器可以免费观看大片| 熟女人妻精品中文字幕| 你懂的网址亚洲精品在线观看| 99久久中文字幕三级久久日本| 成人性生交大片免费视频hd| 免费无遮挡裸体视频| 日韩伦理黄色片| 日本黄色片子视频| 22中文网久久字幕| 丝袜美腿在线中文| 欧美成人精品欧美一级黄| 国内精品美女久久久久久| 一个人免费在线观看电影| 国产在线男女| 黑人高潮一二区| 精品人妻一区二区三区麻豆| 91av网一区二区| 午夜爱爱视频在线播放| 岛国毛片在线播放| 亚洲精品国产成人久久av| 日日啪夜夜爽| 日韩伦理黄色片| 亚洲一区高清亚洲精品| 国内少妇人妻偷人精品xxx网站| 在线观看美女被高潮喷水网站| 久久久久久久久久久丰满| 男女边摸边吃奶| 毛片女人毛片| www.av在线官网国产| 青青草视频在线视频观看| 亚洲天堂国产精品一区在线| 中文字幕av成人在线电影| 天天躁日日操中文字幕| 亚洲国产日韩欧美精品在线观看| 日韩大片免费观看网站| 国产在线一区二区三区精| 亚洲精品日韩av片在线观看| 熟女人妻精品中文字幕| 国产伦在线观看视频一区| 丰满人妻一区二区三区视频av| 欧美3d第一页| 欧美精品一区二区大全| 亚洲最大成人手机在线| 一级毛片黄色毛片免费观看视频| 深爱激情五月婷婷| 日韩电影二区| 国产 亚洲一区二区三区 | 好男人视频免费观看在线| kizo精华| a级毛片免费高清观看在线播放| 亚洲天堂国产精品一区在线| 日韩av在线免费看完整版不卡| 久久97久久精品| 一级爰片在线观看| 七月丁香在线播放| 又爽又黄a免费视频| 亚洲欧美成人综合另类久久久| 亚洲一区高清亚洲精品| 亚洲精品亚洲一区二区| 日韩av在线大香蕉| 久久久久久久久中文| 99久久中文字幕三级久久日本| 天堂俺去俺来也www色官网 | 午夜免费观看性视频| 汤姆久久久久久久影院中文字幕 | 中文字幕亚洲精品专区| 免费观看a级毛片全部| 最近的中文字幕免费完整| 三级国产精品欧美在线观看| 精品久久久久久久末码| 国产伦理片在线播放av一区| 青春草亚洲视频在线观看| 美女被艹到高潮喷水动态| 久久人人爽人人爽人人片va| 国产日韩欧美在线精品| av又黄又爽大尺度在线免费看| 内地一区二区视频在线| 91av网一区二区| 免费看日本二区| 国产av在哪里看| 97热精品久久久久久| 亚洲最大成人中文| 国产男人的电影天堂91| 精品人妻一区二区三区麻豆| 菩萨蛮人人尽说江南好唐韦庄| 色视频www国产| 国产免费福利视频在线观看| 婷婷色综合www| 欧美人与善性xxx| 乱人视频在线观看| 男女那种视频在线观看| 男人舔女人下体高潮全视频| 免费观看的影片在线观看| 色网站视频免费| 欧美丝袜亚洲另类| 99久久精品国产国产毛片| 免费无遮挡裸体视频| 国产午夜精品久久久久久一区二区三区| 97热精品久久久久久| 亚洲成色77777| 激情五月婷婷亚洲| 午夜精品一区二区三区免费看| 国产亚洲5aaaaa淫片| 狂野欧美白嫩少妇大欣赏| 赤兔流量卡办理| 精品久久久久久久久av| 尤物成人国产欧美一区二区三区| 街头女战士在线观看网站| 亚洲精品第二区| 日韩一本色道免费dvd| 午夜爱爱视频在线播放| 在线观看av片永久免费下载| 午夜激情福利司机影院| 国产黄色免费在线视频| 国产中年淑女户外野战色| 日韩av不卡免费在线播放| 日韩 亚洲 欧美在线| 亚洲av国产av综合av卡| 水蜜桃什么品种好| 日韩精品青青久久久久久| 精品人妻熟女av久视频| 亚洲av中文av极速乱| 三级毛片av免费| 国产黄色免费在线视频| 亚洲成人一二三区av| av国产免费在线观看| 韩国高清视频一区二区三区| 久久精品人妻少妇| 国产老妇伦熟女老妇高清| 夜夜看夜夜爽夜夜摸| 性色avwww在线观看| 精品国产一区二区三区久久久樱花 | 777米奇影视久久| 韩国高清视频一区二区三区| 婷婷色综合大香蕉| 女人被狂操c到高潮| 看非洲黑人一级黄片| 高清毛片免费看| 国产成人91sexporn| 国产高清三级在线| 日韩制服骚丝袜av| 久久久久久久久久久丰满| 亚洲av.av天堂| 国产一区有黄有色的免费视频 | 丰满少妇做爰视频| 在线a可以看的网站| 国产在视频线精品| 国产成人精品久久久久久| 五月天丁香电影| 国产永久视频网站| 久久久久性生活片| 色综合色国产| 国产精品女同一区二区软件| 看免费成人av毛片| 我要看日韩黄色一级片| 免费看光身美女| videossex国产| 成人亚洲精品一区在线观看 | av卡一久久| 看十八女毛片水多多多| 一区二区三区乱码不卡18| 麻豆久久精品国产亚洲av| 美女内射精品一级片tv| 欧美日韩综合久久久久久| 色播亚洲综合网| 色网站视频免费| 男人爽女人下面视频在线观看| 天天躁夜夜躁狠狠久久av| 久久99热6这里只有精品| 久久精品国产自在天天线| 日韩欧美国产在线观看| 国产伦精品一区二区三区视频9| 国产久久久一区二区三区| 久久久久精品久久久久真实原创| 婷婷色av中文字幕| 搡女人真爽免费视频火全软件| 午夜精品国产一区二区电影 | 成年女人看的毛片在线观看| 99久久精品国产国产毛片| 久久久久久九九精品二区国产| 亚洲精品成人av观看孕妇| 亚洲精品日韩av片在线观看| 精品人妻一区二区三区麻豆| 国产黄a三级三级三级人| 小蜜桃在线观看免费完整版高清| 深爱激情五月婷婷| 久久鲁丝午夜福利片| 成人午夜高清在线视频| 国产伦一二天堂av在线观看| 午夜亚洲福利在线播放| 中文字幕av在线有码专区| a级毛片免费高清观看在线播放| 久久久欧美国产精品| 免费看美女性在线毛片视频| 精品少妇黑人巨大在线播放| 国产精品久久久久久av不卡| 国产av码专区亚洲av| 97精品久久久久久久久久精品| 国产极品天堂在线| 欧美高清性xxxxhd video| 亚洲av成人av| 亚洲性久久影院| 一二三四中文在线观看免费高清| 精品不卡国产一区二区三区| 日韩av在线大香蕉| 黄色一级大片看看| 日本一二三区视频观看| 午夜福利高清视频| 精品人妻熟女av久视频| 少妇被粗大猛烈的视频| 精品亚洲乱码少妇综合久久| 国产精品日韩av在线免费观看| 国产亚洲精品av在线| 久久6这里有精品| 欧美区成人在线视频| 能在线免费观看的黄片| 黄片无遮挡物在线观看| 高清毛片免费看| 亚洲欧洲日产国产| 岛国毛片在线播放| 久久久色成人| 亚洲精品日韩av片在线观看| 精品久久久噜噜| 深爱激情五月婷婷| 亚洲成色77777| 国产精品久久久久久精品电影小说 | 亚洲怡红院男人天堂| 黄色欧美视频在线观看|