陳哲 黃靜 趙佳 梁宏

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

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

中圖分類號： S 436.63

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

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

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

1 材料與方法

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

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

1.2 培養(yǎng)基

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

1.3 病原菌侵染試驗

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

在處理組接種病原菌50 d后進行樣品采集?；疾∪~柄的選擇：截取包含病斑的10 mm長的葉柄段，病斑長度為5.1～10.0 mm。每株草莓只收集1個葉柄，5株草莓即5個葉柄段為1組，處理組和對照組各3組重復。采集到的樣品在液氮迅速冷凍，-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堿基隨機引物和逆轉(zhuǎn)錄酶合成cDNA第一鏈，并以第一鏈cDNA為模板合成第二鏈cDNA，第二鏈cDNA合成時，其中的堿基T，被替換成U，從而達到鏈特異性文庫的目的。

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

1.4.3 高通量測序

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

提取每個基因下最長的轉(zhuǎn)錄本作為該基因的代表序列，稱為Unigene。對 Unigene 進行基因功能注釋?；蚬δ茏⑨屗玫降臄?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對基因表達進行差異分析，篩選差異表達基因條件為：表達差異倍數(shù)|log2FoldChange|>1，顯著性 P-value<0.05。

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

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

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

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

1.4.7 差異表達基因的熒光定量PCR（qRT-PCR）驗證

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

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

2 結(jié)果與分析

2.1 病原菌侵染草莓植株

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

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

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

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

2.3 差異表達基因的篩選

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

2.4 差異表達基因的GO富集分析

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

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

2.5 差異表達基因的 Pathway 富集分析

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

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

2.6 差異表達基因的qRT-PCR驗證

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

19個差異表達基因中包含10個上調(diào)基因和9個下調(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）、生長素誘導蛋白15A（AX15A）、細胞分裂素脫氫酶5（CKX5）、ADP-葡萄糖焦磷酸化酶大亞基（gpL2）、類胡蘿卜素雙加氧酶4（CCD4）、9-順式環(huán)氧類胡蘿卜素雙加氧酶1（NCED1）、核酮糖1，5-二磷酸羧化酶/加氧酶（rbcL）和葉綠體果糖-1，6-二磷酸酶Ⅰ（cfbI）。

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

3 討論

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

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

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

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

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

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

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

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

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

3.3 植物正常生長受到抑制

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

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

參考文獻

[1] 陳哲， 黃靜， 趙佳， 等. 防治草莓炭疽病的芽胞桿菌組合的篩選[J]. 中國生物防治學報， 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ù)在植物水淹脅迫研究中的應用[J]. 分子植物育種： 2019（2）： 113.

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

[7] 劉之慧， 蔡雯婷， 董清華. 紅藍光處理草莓組培苗的轉(zhuǎn)錄組測序與分析[J]. 北京農(nóng)學院學報， 2017， 32（4）： 2835.

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

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

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

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

[12]谷春艷， 蘇賢巖， 楊雪， 等. 解淀粉芽胞桿菌WH1G與咪鮮胺協(xié)同防治草莓炭疽病[J]. 植物保護， 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]. 園藝學報， 2019， 46（1）： 121134.

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

[17]楊仕美， 喬光， 毛永亞， 等. 基于火龍果轉(zhuǎn)錄組測序的SSR標記開發(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）對生物脅迫與非生物脅迫的響應[J].熱帶作物學報，2019，40（9）：17421751.

[20]王福祥， 肖開轉(zhuǎn)， 姜身飛， 等. 干旱脅迫下植物體內(nèi)ROS的作用機制[J].科學通報，2019，64（17）：17651779.

[21]王文佳，李爽，馬澤眾，等.水分脅迫對春大豆葉片保護酶活性及相對電導率的影響[J].中國農(nóng)學通報，2019，35（11）：1418.

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

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

[24]黃啟秀，曲延英，姚正培，等. 海島棉枯萎病抗性與類黃酮代謝途徑基因表達量的相關(guān)性[J].作物學報，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ù)學報，2018，26（7）：11241131.

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

[28]張弢， 董春海. 乙烯信號轉(zhuǎn)導及其在植物逆境響應中的作用[J]. 生物技術(shù)通報， 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]張歡，楊英杰，李鼎立，等.杜梨根莖葉特異表達基因的RNA-Seq分析[J].園藝學報，2018，45（10）：18811894.

[31]張晨， 鄭偉， 周濤， 等. 太子參細胞分裂素氧化/脫氫酶基因PhCKX的克隆及表達分析[J]. 分子植物育種， 2018， 16（19）： 62426249.

[32]李志康，嚴冬，薛張逸，等.細胞分裂素對植物生長發(fā)育的調(diào)控機理研究進展及其在水稻生產(chǎn)中的應用探討[J].中國水稻科學，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.

（責任編輯：田 喆）

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

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

致? 謝： 參加本試驗部分工作的還有江代禮、譚翰杰、張能和紀燁斌等同學，特此一并致謝。

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

#為并列第一作者