紅豆杉健康與感葉枯病針葉內(nèi)生微生物多樣性

王俊凱 劉崢 申東晨 董愛榮

摘 要：采用Illumina Miseq高通量測序技術(shù)來探究人工林與天然林東北紅豆杉（Taxus cuspidata）健康針葉與感葉枯病的病葉之間內(nèi)生真菌和細菌群落結(jié)構(gòu)的差異。結(jié)果表明，紅豆杉針葉感病前后真菌的OTU （Operational Taxonomic Units）變化較小，細菌的OTU變化較大，天然林比人工林具有更多的真菌OTU數(shù)，細菌OTU數(shù)差距較小。4組樣品中共有相同的真菌OTU 9個，占真菌OTU 總數(shù)的9.37%，4組樣品中共有相同的細菌OTU 34個，占細菌OTU總數(shù)的4.19%。在門水平下，感病前后樣品中的真菌群落結(jié)構(gòu)組成基本相同，差異較大的是相對多度，相比于天然林，子囊菌門（ Ascomycota ）在人工林中的相對多度更高，子囊菌門在4組樣品占據(jù)絕對優(yōu)勢。感病前后樣品中的細菌群落結(jié)構(gòu)組成及相對多度都發(fā)生了較大變化，相比于天然林，人工林中細菌種類更多。在屬水平下，Genolevuria屬是人工林健康葉中最主要的優(yōu)勢真菌，穴殼屬（Dothiora）是人工林病葉中最主要的優(yōu)勢真菌，Buckleyzyma屬是天然林健康葉已知屬中最主要的優(yōu)勢真菌，Genolevuria屬是天然林病葉已知屬中最主要的優(yōu)勢真菌。勞爾氏菌屬（Ralstonia）是人工林健康葉中最主要的優(yōu)勢細菌，1174-901-12屬是人工林病葉中最主要的優(yōu)勢細菌，鞘氨醇單胞菌屬（Sphingomonas）是天然林健康葉已知屬中最主要的優(yōu)勢細菌，馬賽菌屬（Massilia）是天然林病葉已知屬中最主要的優(yōu)勢細菌。推測出鏈格孢屬（Alternaria）和盤雙端毛孢屬（Seimatosporium）可能會導(dǎo)致紅豆杉的葉枯病。

關(guān)鍵詞：高通量測序；紅豆杉；葉枯??；真菌多樣性；細菌多樣性

中圖分類號：S763.1 ???文獻標(biāo)識碼：A ??文章編號：1006-8023（2023）04-0010-09

Diversity of Endophytic Microorganisms in Healthy and Susceptible

Leaf Blight Needles of Taxus cuspidata

WANG Junkai， LIU Zheng， SHEN Dongchen， DONG Airong*

（School of Forestry， Northeast Forestry University， Harbin 150040， China）

Abstract：：Illumina Miseq high-throughput sequencing technology was used to investigate the differences in endophytic fungal and bacterial community structure between healthy and diseased needles of Taxus cuspidata in artificial and natural forests. The results showed that the changes of fungal OTU （Operational Taxonomic Units） were small and the changes of bacterial OTU were large. The natural forest had more fungal OTU number than the artificial forest， and the difference of bacterial OTU number was smaller. There were 9 fungal OTUs in the four groups， accounting for 9.37% of the total fungal OTU. There were 34 identical bacterial OTUs in the four groups of samples， accounting for 4.19% of the total bacterial OTU. At the phylum level， the fungal community composition in the samples before and after infection was basically the same， but degree of relative abundance was significantly different. Compared with the natural forest， the relative abundance of Ascomycota was higher in the artificial forest. Ascomycota was dominant in the four groups of samples. The bacterial community structure and relative abundance changed greatly before and after infection. There were more bacterial species in artificial forests than in natural forests. At the genus level， Genolevuria was the most dominant fungus in healthy leaves of plantations. Dothiora was the most dominant fungus in diseased leaves of plantations， Buckleyzyma was the most dominant fungus in known genera of healthy leaves of natural forests. Genolevuria was the most dominant fungus in the known genera of diseased leaves in natural forests. Ralstonia was the most dominant bacteria in healthy leaves of plantations， 1174-901-12 was the most dominant bacteria in diseased leaves of plantations， Sphingomonas was the most dominant bacteria in known genera of healthy leaves of natural forests. Massilia was the most dominant bacterium in the known genus of natural forest leaf disease. It was speculated that Alternaria and Seimatosporium may cause leaf blight of Taxus cuspidata.

Keywords：High-throughput sequencing; Taxus cuspidate; leaf blight; fungal diversity;bacterial diversity

收稿日期：2022-10-08

基金項目：國家自然科學(xué)基金項目（31670494）

第一作者簡介：王俊凱，碩士研究生。研究方向為森林病理學(xué)。E-mail： wangjk19970912@163.com

通信作者：董愛榮，博士，副教授。研究方向為森林病理學(xué)。E-mail： darlmy@tom.com

引文格式：王俊凱，劉崢，申東晨，等. 紅豆杉健康與感葉枯病針葉內(nèi)生微生物多樣性[J].森林工程，2023，39（4）：10-18.

WANG J K， LIU Z， SHEN D C， et al. Diversity of endophytic microorganisms in healthy and susceptible leaf blight needles of Taxus cuspidata[J]. Forest Engineering， 2023，39（4）：10-18.

東北紅豆杉（Taxus cuspidata）為紅豆杉屬（Taxus）植物，是珍貴的第三紀(jì)孑遺樹種[1]，為國家一級重點保護樹種[2， 3]，具有極高的經(jīng)濟價值、觀賞價值和生態(tài)效益[4]。東北紅豆杉終年常綠，是優(yōu)秀的園林綠化樹種，可以吸附有害氣體，釋放植物精氣，驅(qū)蚊防蟲，抗病蟲害能力強[5-6]。在 20 世紀(jì) 70 年代，Wani等[7]的研究團隊在對名為短葉紅豆杉（Taxus brevifolia）的植物成分進行研究后，成功提取出具有較為明顯抗癌效果的紫杉醇，紅豆杉屬植物便成為眾多學(xué)者研究熱點。東北紅豆杉由于其地理環(huán)境以及氣候因素的特殊條件，紫杉醇的平均含量明顯高于其他種的含量[8-9]。東北紅豆杉天然種群本身生長較緩慢，天然結(jié)實率低，且對生長環(huán)境要求較嚴(yán)，種群競爭力較弱[1， 10]，再加上外界環(huán)境影響，東北紅豆杉天然林已瀕臨滅絕，其保護形勢十分嚴(yán)峻[11-13]。

紅豆杉葉枯病會給東北紅豆杉帶來嚴(yán)重的經(jīng)濟損失[14]。2021年在黑龍江省穆棱市穆棱鎮(zhèn)人工林內(nèi)發(fā)現(xiàn)了紅豆杉葉枯病的發(fā)生，經(jīng)調(diào)查發(fā)現(xiàn)，該病害在染病針葉上發(fā)病，并最終導(dǎo)致針葉枯死。近些年來，隨著生物防治的逐漸完善，運用植物內(nèi)生菌拮抗病原菌來防治葉枯病自然是一種較為可行的方法[13， 15]。高通量測序能夠更為準(zhǔn)確快捷了解微生物的性質(zhì)及其在自然界中的生態(tài)功能，通過高通量測序技術(shù)能夠一次性檢測出樣品中可能感染的多種病原物基因組序列， 對病原檢測鑒定具有很好的應(yīng)用[16- 17]。本研究基于高通量測序技術(shù)，分析人工林與天然林感紅豆杉葉枯病與健康針葉內(nèi)生微生物的多樣性以及群落結(jié)構(gòu)的變化，探知優(yōu)勢微生物的分布，從而可用于指導(dǎo)篩選紅豆杉針葉內(nèi)生優(yōu)勢的葉枯病拮抗菌，并探究人工林與天然林內(nèi)生微生物的差異[18-20]。

1 材料與方法

1.1 研究區(qū)概況

采樣地點位于黑龍江省穆棱市穆棱鎮(zhèn)東北紅豆杉人工林及天然林。東北紅豆杉天然林25 hm2樣地，位于穆棱東北紅豆杉國家級自然保護區(qū)內(nèi)（43°49′～44°06′ N，130°00′～130°28′ E）。樣地地形走勢西南高，海拔最高781 m，最低658 m，高差123? m；年平均溫度在-2? ℃，年平均降水量530? mm，土壤為暗棕壤。樣地內(nèi)森林群落優(yōu)勢種明顯，垂直層次復(fù)雜，是典型的中齡復(fù)層異齡林，其主要建群樹種為紫椴（Tilia amurensis）、紅松（Pinus koraiensis）、臭冷杉（Abies nephrolepis）和楓樺（Betula costata）等，樣地內(nèi)共有胸徑（DBH）大于等于1 cm木本植物57種，隸屬于22科38屬。獨立個體數(shù)為63 877株，平均胸徑為7.83? cm，胸高斷面積26.4? m2。樣地內(nèi)現(xiàn)有東北紅豆杉（Taxus cuspidate）成樹169株，天然更新幼苗350～450株。

東北紅豆杉人工林（43°45′ N，131°00′ E）約4 hm2，1995年營建， 土壤為暗棕壤，東北紅豆杉20 a以上開始結(jié)實，40 a以上果實質(zhì)量高，株行距為2.5? m×3? m，主要建群樹種為落葉松（Larix gmelinii），樣地內(nèi)現(xiàn)有東北紅豆杉共計一百多萬株。

1.2 試驗材料

本實驗分別于人工林及天然林采集健康紅豆杉及感葉枯病紅豆杉針葉。2021 年6月18日至20日，采用“S”形十點采樣法，分別剪取 10 棵 15年生健康及感葉枯病紅豆杉針葉。從 10 棵樹上采集的樣品中分別選取 1.0? g，用無菌剪刀剪碎，將其混勻后，使用75%乙醇表面消毒1? min，無菌水沖洗1次，再使用5%（有效氯）次氯酸鈉溶液表面消毒5 min，無菌水沖洗3次，從而得到人工林及天然林各 10.0? g 健康樣品及 10.0? g 感病樣品。穆棱鎮(zhèn)人工林健康樣品、感葉枯病樣品分別用S1、L1表示；穆棱鎮(zhèn)天然林健康樣品、感葉枯病樣品分別用S2、L2表示。將樣品分別置于2? mL的無菌離心管中并于－80? ℃的超低溫冰箱中保存至DNA提取。

1.3 紅豆杉內(nèi)生微生物基因組 DNA 的提取及質(zhì)量檢測

紅豆杉樣品的基因組DNA參照 MP 試劑盒 Fast DNA Spin Kit for Soil 116560-200 的使用說明書進行提取，將所提取到的 DNA 溶于50? μL無菌水中，混合均勻，置于－80? ℃保存?zhèn)溆?。之后使用NanoDrop2000進行基因組DNA濃度及純度檢測，使用瓊脂糖凝膠電泳進行基因組DNA完整性檢測。

1.4 PCR擴增和高通量測序

利用提取的 DNA進行PCR擴增，真菌選擇ITS1 和ITS2R 作為引物，分別為 5′－CTTGGTCATTTAGAGGAAGTAA－3′和 5′－GCTGCGTTCTTCATC-GATGC－3′；細菌選擇 799F 和 1193R 作為引物，引物序列分別為 5′－AACMGGATTAGATACCCKG－3′和 5′－ACGTCATCCCCACCTTCC－3′。將其送至上海美吉生物醫(yī)藥科技有限公司進行高通量測序，利用 Illumina MiSeq 平臺進行測序和數(shù)據(jù)分析。

1.5 數(shù)據(jù)分析

下機數(shù)據(jù)經(jīng)過 QIIME（ v1．8．0） 軟件過濾、拼接和去除嵌合體，去除各樣本中reads 尾部質(zhì)量值在 20 以下的堿基、切除 reads 中含N 部分序列，并去除數(shù)據(jù)中的短序列（ 長度閾值 200? bp） ，隨后再對低復(fù)雜度的序列進行過濾。采用Usearch（ v7.1） 軟件進行數(shù)據(jù)處理，物種比對注釋使用RDP classifier 軟件，保留置信區(qū)間大于0.8 的注釋結(jié)果。用R語言工具繪制稀釋性曲線、韋恩（Venn）圖、群落結(jié)構(gòu)圖及主坐標(biāo)分析圖（Principal Co-ordinates Analysis，PcoA）。利用 Mothur 軟件（version v.1.30.2 http：//www.mothur.org/wiki/Schloss_ SOP#Alpha_diversity）進行 Chao1 （估計物種總數(shù)）指數(shù)、香農(nóng)指數(shù)（ Shannon）、氣旋能量指數(shù)（Accumulated cyclone energy， Ace）、辛普森（Simpson）指數(shù)、Coverage（覆蓋度）指數(shù)計算分析，并在各分類水平上進行群落結(jié)構(gòu)的統(tǒng)計分析，得到微生物群落結(jié)構(gòu)組成。生境內(nèi)的多樣性指數(shù) （within-habitat diversity，Alpha）可用來反映生物群落的豐富度和多樣性。其中，Ace、Chao1指數(shù)可反映生物物種的豐富度，其數(shù)值越大表示該樣品物種豐富度越大，其中Ace指數(shù)用來估計樣本中物種組成的均勻度的指數(shù)，Chao1指數(shù)用來估計樣本中所含OTU（Operational Taxonomic Units）數(shù)目的指數(shù)；Shannon、Simpson指數(shù)可以反映生物群落的多樣性，其中Shannon指數(shù)用來衡量群落的異質(zhì)性，其數(shù)值越大表示群落多樣性越高，Simpson指數(shù)代表隨機取樣的2個體屬于不同種的概率，其數(shù)值越大表示群落物種分布越均勻；Coverage指數(shù)則可以反映生物群落的覆蓋度，其數(shù)值越高，則樣本中序列被測出的概率越高，該指數(shù)反映本次測序結(jié)果是否代表了樣本中微生物的真實情況。生境間多樣性 （between-habitat diversity，Beta）分析通過對不同微生物群落間的物種多樣性進行組間比較分析，探索不同分組樣本間群落組成的相似性或差異性。

2 結(jié)果與分析

2.1 測序質(zhì)量分析

通過對人工林及天然林健康紅豆杉及感葉枯病紅豆杉針葉4種樣品進行高通量測序數(shù)據(jù)統(tǒng)計，S1、L1、S2、L2分別得到102 022、87 750、100 624、65 973條有效真菌序列及16 090、15 336、18 905和15 983條有效細菌序列，采用對測序獲得的序列進行隨機抽樣的方法， 以抽到的序列數(shù)與其所代表的OTU數(shù)構(gòu)建稀釋曲線（圖1）。由圖1可以看出，4種針葉內(nèi)生微生物的稀釋曲線均趨于平坦，表明測序數(shù)據(jù)能很好地反映各樣品中真菌及細菌的群落結(jié)構(gòu)和多樣性，更多的測序數(shù)據(jù)對發(fā)現(xiàn)新的OTU貢獻率較小。

將測序所得的有效序列按97%的相似度進行OTU的劃分后整理繪制的Venn圖（圖2）。由圖2（a）可知，在4組處理中檢測到的真菌總OTU數(shù)為950個，其中S1組有222個、L1組有204個、S2組有306個、L2組有218個。4組樣品中共有相同的OTU89個，占 OTU 總數(shù)的9.37%，4組樣品中特有的 OTU 數(shù)目分別為58、17、74、31個；由圖2（b）可知，在4組處理中檢測到的細菌總OTU數(shù)為811個，其中S1組有274個、L1組有121個、S2組有279個、L2組有137個。4組樣品中共有相同的OTU34個，占 OTU 總數(shù)的4.19%，4組樣品中特有的 OTU 數(shù)目分別為144、17、100、17個，以上結(jié)果表明天然林比人工林具有更多的真菌與細菌OTU數(shù)，細菌OTU數(shù)量差距較??；感病后，針葉的真菌與細菌OTU數(shù)量減少。

2.2 感病前后紅豆杉針葉的微生物群落門類組成

在97％的相似度水平對樣品序列進行OTU聚類，由表1可知，4組紅豆杉針葉微生物多樣性豐富，將不同樣品中相同的真菌進行合并，共鑒定得出真菌6個門、24個綱、71個目、164個科、300個屬、427個種，將不同樣品中相同的細菌進行合并，共鑒定得出細菌23個門、47個綱、104個目、161個科、248個屬、350個種。其中，人工林健康針葉與天然林健康針葉的真菌及細菌種類數(shù)量明顯高于人工林感病針葉與天然林感病針葉（人工林健康針葉的真菌種類數(shù)量稍高于人工林感病針葉）。

在相似水平為97%的條件下，與16S嵌合體數(shù)據(jù)庫（Gold database）及ITS嵌合體數(shù)據(jù)庫：（UNITE）對比，對所得序列進行聚類分析， 圖3為感病前后人工林與天然林內(nèi)生真菌和細菌在門水平上的群落結(jié)構(gòu)特征（物種多度占比大于0.1%的門）。由圖3（a）可知，從真菌分類門水平看，人工林和天然林4種樣品檢測出3個相對多度較高的真菌類群，即： 子囊菌門（ Ascomycota）、擔(dān)子菌門（ Basidiomycota）及未鑒定出的真菌（ unclassified_k__Fungi）。其中S1、L1、S2和L2中子囊菌門的相對多度依次降低，擔(dān)子菌門的相對多度依次升高。S1 和L1中Ascomycota的相對多度分別為83.02%與78.81%，明顯高于S2和L2；S2和L2中unclassified_k_Fungi的相對多度分別為29.95%與33.37%，明顯高于S1和L1。結(jié)果表明，在門水平下，感病前后樣品中的真菌群落結(jié)構(gòu)組成基本相同，差異較大的是相對多度；相比于天然林，Ascomycota在人工林中的相對多度更高。

由圖3（b）可知，從細菌分類門水平看，人工林和天然林4種樣品主要檢測出6個相對多度較高的細菌類群，即： 變形菌門（Proteobacteria）、放線菌門（Actinobacteriota）、厚壁菌門（Firmicutes）、擬桿菌門（Bacteroidota）、粘球菌門 （ Myxococcota）和酸桿菌門（ Acidobacteriota）。其中Proteobacteria與Actinobacteriota為4組針葉共有的優(yōu)勢細菌，S1、L1、S2和L2中Proteobacteria的相對多度分別為67.44%、35.26%、52.24%和86.29%；Actinobacteriota的相對多度分別為8.14%、58.38%、41.39%和11.38%。Myxococcota為L1 及S2共有的優(yōu)勢細菌，分別占比4.74%和4.66%。Firmicutes與Bacteroidota只在S1中為優(yōu)勢細菌，Acidobacteriota只在L1中為優(yōu)勢細菌。結(jié)果表明，在門水平下，感病前后樣品中的細菌群落結(jié)構(gòu)組成及相對多度都發(fā)生了較大變化；相比于天然林，人工林中細菌種類更多。

根據(jù)物種分類結(jié)果，由于樣品中所檢測出的真菌及細菌種類繁多，許多物種含量很少，因此篩選出優(yōu)勢物種，并對優(yōu)勢真菌及優(yōu)勢細菌進行分類統(tǒng)計并將在2組樣品中含量均小于1%的物種合并至其他類別。由圖4（a）可知，從真菌分類屬水平看，Genolevuria屬為S1、L1、S2和L2中共有的優(yōu)勢真菌，占比為4.97%～16.62%。S1中特有的優(yōu)勢真菌為尾孢屬（Cercospora）、副孢霉屬（Paramycosphaerella）、柱隔孢屬（Ramularia）、Mycodiella屬、膜盤菌屬（Pezoloma）、小莢孢腔菌屬（Sporormiella）和Arthrocatena屬；L1中特有的優(yōu)勢真菌為鏈核盤菌屬（Monilinia）、穴殼屬（Dothiora）和白孔座殼菌屬（Leucostoma）；S2中特有的優(yōu)勢真菌為Buckleyzyma屬和Retiarius屬；L2中特有的優(yōu)勢真菌為銀耳屬（Tremella）、酵母菌屬（Mrakia）、彎擔(dān)菌屬（Curvibasidium）、附球菌屬（Epicoccum）、褐銀耳屬（Phaeotremella）、鐮孢屬（Fusarium）和赤霉菌屬（Gibberella）。健康樣品S1和S2中特有的真菌為維希尼克氏酵母屬（Vishniacozyma）、柱隔孢屬（Ramularia）、Dioszegia屬和Arthrocatena屬；感病樣品L1和L2中特有的真菌為鏈格孢屬（Alternaria）和盤雙端毛孢屬（Seimatosporium）；人工林樣品S1和L1中特有的真菌為外囊菌屬（Taphrina）、短梗霉屬（Aureobasidium）、亞隔孢殼屬（Didymella）和Claussenomyces屬；天然林樣品S2和L2中特有的真菌為Buckleyzyma屬、多臂菌屬（Trichomerium）和Retiarius屬。

由圖4（b）可知，4組樣品中的細菌群落結(jié)構(gòu)組成不盡相同，且相對多度的差異也較大。從細菌分類屬水平看，鞘氨醇單胞菌屬（Sphingomonas）為S1、L1、S2和L2中共有的優(yōu)勢細菌，占比分別為1.94%、2.71%、14.23%和6.42%。S1中特有的優(yōu)勢細菌為勞爾氏菌屬（Ralstonia）、克雷伯菌屬（Klebsiella）、地芽孢桿菌屬（Geobacillus）、擬桿菌屬（Bacteroides）、代爾夫特菌屬（Delftia）、噬氫菌屬（Ottowia）、大腸埃希菌-志賀菌屬（Escherichia-Shigella）、湖沉積桿菌屬（Limnobacter）、norank_f__PHOS-HE36屬、黃桿菌屬（Flavobacterium）、Clostridium_sensu_stricto_8屬和亞硝化單胞菌（Nitrosomonas）；L1中特有的優(yōu)勢細菌為內(nèi)桿菌屬（Endobacter）；S2中特有的優(yōu)勢細菌為動球菌屬（Kineococcus）、中村氏菌屬（Nakamurella）、湖沼菌屬（Limnohabitans）、norank_f__norank_o__Gaiellales屬、Klenkia屬和諾卡氏菌屬（Nocardioides）；L2中特有的優(yōu)勢細菌為嗜酸菌屬（Acidiphilium）。健康樣品S1和S2中特有的細菌為勞爾氏菌屬（Ralstonia）、貪噬菌屬（Variovorax）、擬桿菌屬（Bacteroides）、中村氏菌屬（Nakamurella）、湖沼菌屬（Limnohabitans）、norank_f__norank_o__Gaiellales屬和Klenkia屬；感病樣品L1和L2中特有的細菌為嗜酸菌屬（Acidiphilium）；天然林樣品S2和L2中特有的細菌為弗萊德門菌屬（Friedmanniella）和克雷伯菌屬（Klebsiella）。

2.3 群落結(jié)構(gòu)多樣性分析

2.3.1 Alpha多樣性

由表2可知，本研究將相似度為97%的序列劃分為1個 OTU，進行多樣性指數(shù)計算，并以5個衡量指數(shù)詳細列出了4組樣品中微生物的豐富度和多樣性分布情況。所有樣品真菌及細菌的覆蓋度（Coverage）均超過 99.67%，表明真菌及細菌的測序深度已達到較高水平，數(shù)據(jù)可靠。

感病前后人工林與天然林4組紅豆杉針葉樣品中真菌與細菌群落的 Alpha 多樣性指數(shù)部分有差異。 4組樣品的真菌與細菌群落中S2（天然林健康樣品）的 Chao1指數(shù)最高，即測得OTU 數(shù)最高，其次分別為 L2（天然林感病樣品）、S1（人工林健康樣品）和L1（人工林感病樣品），樣品S1及L1的Chao1指數(shù)差異不顯著，Ace多樣性指數(shù)所表現(xiàn)的多樣性規(guī)律與 Chao1多樣性指數(shù)結(jié)果一致。表明天然林樣品的真菌與細菌豐富度指數(shù)高于人工林樣品；在天然林及人工林中，健康樣品的真菌與細菌豐富度指數(shù)高于感病樣品，但人工林健康樣品與感病樣品相差不大。

4組樣品的真菌與細菌群落中S1的 Shannon指數(shù)最高，其次分別為 S2、L2、L1，4組樣品的真菌與細菌群落中Shannon指數(shù)均差異較大，Simpson多樣性指數(shù)所表現(xiàn)的多樣性規(guī)律與Shannon多樣性指數(shù)結(jié)果相反。表明健康樣品的真菌與細菌多樣性指數(shù)高于感病樣品；相比于天然林，在人工林中，健康樣品的真菌與細菌多樣性指數(shù)明顯高于感病樣品。

2.3.2 Beta多樣性

PCoA分析可用來研究樣本群落組成的相似性或差異性，此分析以點的形式將多組數(shù)據(jù)的差異性反映在二維坐標(biāo)圖上，點與點之間距離的遠近代表樣品中群落的差異程度。由圖5可以看出，S1與L1、S2與L2的微生物構(gòu)成差距較大，即同地區(qū)健康與感病樣品的微生物構(gòu)成差距較大；S1與S2、L1與L2的微生物構(gòu)成差距較大，即不同地區(qū)同為健康與感病樣品的微生物構(gòu)成差距也較大。

3 討論

本研究通過Illumina MiSeq測序平臺，分別對人工林與天然林紅豆杉健康葉和感葉枯病針葉內(nèi)生微生物的群落結(jié)構(gòu)進行了分析，用數(shù)據(jù)說明了上述兩者之間共生細菌和真菌群落結(jié)構(gòu)的差異。關(guān)于病害對微生物群落的影響，Wei等[21]發(fā)現(xiàn)病原菌侵染微生物群落會導(dǎo)致群落的多樣性下降，而本研究得出了相同結(jié)論：健康針葉的Ace、Chao1及Shannon多樣性指數(shù)均高于感病針葉，可能是病原菌侵入破壞了原有微生物群落的平衡。

本研究發(fā)現(xiàn)，相比于健康樣品，人工林與天然林感病樣品中特有的真菌為鏈格孢屬（Alternaria）和盤雙端毛孢屬（Seimatosporium）。目前盤雙端毛孢屬（Seimatosporium）已記錄有38種，這一屬中的某些種可以致病，比如大白花杜鵑病害[22]。商娜等[23]確定白三葉草葉斑病病原菌為喬木鏈格孢（Alternaria arborescens），楊志英等[24]確定茶葉斑病病原菌為長柄鏈格孢菌（Alternaria longipes），劉俏等[25]確定櫻桃葉斑病病原菌為鏈格孢（Alternaria alternata），目前鏈格孢屬真菌已記錄有500多種，且95%寄生在不同植物上，鏈格孢屬真菌寄主范圍廣，侵染能力強并且危害程度大。通過分析發(fā)現(xiàn)鏈格孢屬（Alternaria）和盤雙端毛孢屬（Seimatosporium）可能為紅豆杉葉枯病的病原菌。

利用內(nèi)生菌進行植物病害生物防治是近年來的研究熱點之一。研究發(fā)現(xiàn)，貪噬菌屬（Variovorax）在土壤與植物環(huán)境中廣泛分布，且該菌具有固氮作用，目前并未作為致病病原細菌被發(fā)現(xiàn)，陳丹[26]發(fā)現(xiàn)Variovorax可以作為拮抗菌抑制灰葡萄孢霉菌。故接下來的研究方向可以著重研究Variovorax對于病原菌屬的拮抗能力及拮抗原理。

4 結(jié)論

紅豆杉發(fā)生紅豆杉葉枯病后，感葉枯病病葉與健康葉相比，細菌、真菌的多樣性及物種總數(shù)都顯著下降，在門及屬水平上差異都較為顯著；天然林比人工林具有更多的真菌與細菌物種數(shù)，細菌物種數(shù)差距較小。本研究確定了人工林與天然林內(nèi)生微生物的群落結(jié)構(gòu)組成及相對多度，明確了感病前后，病葉與健康葉在門與屬水平上的具體差異，推測出鏈格孢屬（Alternaria）和盤雙端毛孢屬（Seimatosporium）可能會導(dǎo)致東北紅豆杉的葉枯病，并通過查閱文獻，確定了一種可能拮抗病原菌的細菌：貪噬菌屬。研究結(jié)果為接下來拮抗菌拮抗紅豆杉病的研究與應(yīng)用奠定了基礎(chǔ)。

【參 考 文 獻】

[1]徐超，龍婷，吳鑫磊，等.東北紅豆杉種群的回歸成效及影響因素[J].北京林業(yè)大學(xué)學(xué)報，2020，42（8）：34-42.

XU C， LONG T， WU X L， et al. Regression effect and influencing factors of Taxus chinensis population[J]. Journal of Beijing Forestry University， 2020， 42（8）：34-42.

[2]臧潤國，董鳴，李俊清，等.典型極小種群野生植物保護與恢復(fù)技術(shù)研究[J].生態(tài)學(xué)報，2016，36（22）：7130-7135.

ZANG R G， DONG M， LI J Q， et al. Study on conservation and restoration techniques of wild plants with typical minimal population[J]. Acta Ecologica Sinica， 2016， 36（22）：7130-7135.

[3]葛麗麗，黃云峰，薛俊剛，等.東北紅豆杉雌、雄球花發(fā)育期內(nèi)源激素含量動態(tài)變化[J].北華大學(xué)學(xué)報（自然科學(xué)版），2021，22（4）：456-461.

GE L L， HUANG Y F， XUE J G， et al. Dynamic changes of endogenous hormones in female and male Taxus cuspus during development[J]. Journal of Beihua University （Natural Science Edition）， 2021， 22（4）：456-461.

[4]周志強，胡丹，劉彤.天然東北紅豆杉種群生殖力與開花結(jié)實特性[J].林業(yè)科學(xué)，2009，45（5）：80-86.

ZHOU Z Q， HU D， LIU T. Fecundity and characteristics of flowering and fruiting of natural Taxus cuspidata population[J]. Scientia Silvae Sinicae， 2009， 45（5）：80-86.

[5]張朕.東北紅豆杉新品種黃金杉繁殖過程中的關(guān)鍵技術(shù)[D].長春：長春大學(xué)，2021.

ZHANG Z. Key technologies in the propagation process of a new Taxus chinensis variety[D]. Changchun： Changchun University， 2021.

[6]吳磊.園林植物精氣的檢測方法及高效健康植物的篩選[D].長春：長春師范大學(xué)，2019.

WU L. Detection method of garden plant essence and screening of efficient healthy plants[D]. Changchun： Changchun Normal University， 2019.

[7]WANI M C， TAYLOR H L， WALL M E， et al. Plant antitumor agents. Ⅵ. The isolation and structure of taxol， a novel antileukemic and antitumor agent from Taxus brevifolia[J]. Journal of the American Chemical Society， 1971， 93（9）：2325-2327.

[8]崔海超，鄭文藝，張智慧，等.不同紅豆杉屬植物中7種紫杉烷類化合物的HPLC檢測與分析[J].森林工程，2022，38（4）：118-124.

CUI H C， ZHENG W Y， ZHANG Z H， et al. Determination and analysis of seven taxanes from Taxus species by HPLC[J]. Forest Engineering， 2022， 38（4）：118-124.

[9]王書凱，陳凡.東北紅豆杉中紫杉醇含量的研究[J].林業(yè)科技通訊，1999（7）：18-19.

WANG S K， CHEN F. Study on the content of taxol in Taxus cuspidate[J]. Journal of Forestry Science and Technology， 1999（7）：18-19.

[10]吳世雄，劉艷紅，張利民，等.不同產(chǎn)地東北紅豆杉幼苗遷地保護的生長穩(wěn)定性分析[J].北京林業(yè)大學(xué)學(xué)報，2018，40（12）：27-37.

WU S X， LIU Y H， ZHANG L M， et al. Growth stability analysis of Taxus cuspidate seedlings from different habitats under ex situ protection[J]. Journal of Beijing Forestry University， 2018， 40（12）：27-37.

[11]周志強，劉彤，胡林林，等.穆棱東北紅豆杉年輪-氣候關(guān)系及其瀕危機制[J].生態(tài)學(xué)報，2010，30（9）：2304-2310.

ZHOU Z Q， LIU T， HU L L， et al. Tree-ring-climate relationship and endangered mechanism of Taxus cuspidate in Muling[J]. Acta Ecologica Sinica， 2010， 30（9）：2304-2310.

[12]姬少玲.東北紅豆杉繁育技術(shù)研究進展[J].遼寧林業(yè)科技，2015（4）：61-63.

JI S L. Research progress of Taxus chinensis breeding technology[J]. Liaoning Forestry Science and Technology， 2015（4）：61-63.

[13]劉亞苓.兩種拮抗菌對細辛葉枯病的防治效果及機制研究[D].北京：中國農(nóng)業(yè)科學(xué)院，2020.

LIU Y L. Study on the control effect and mechanism of two antagonists against Asarum blight[D]. Beijing： Chinese Academy of Agricultural Sciences， 2020.

[14]李陽.紅豆杉葉枯病病原鑒定及防治藥劑篩選研究[D].福州：福建農(nóng)林大學(xué)，2017.

LI Y. Pathogen identification and control agent screening of Taxus cuspidate leaf blight[D]. Fuzhou： Fujian Agriculture and Forestry University， 2017.

[15]于占晶，侯曉杰，高群英，等.6種細菌對蘋果炭疽葉枯病病菌的抑制作用[J].江蘇農(nóng)業(yè)科學(xué)，2022，50（3）：121-124.

YU Z J， HOU X J， GAO Q Y， et al. Inhibitory effects of 6 types of bacteria on Anthrax leaf blight of apple[J]. Jiangsu Agricultural Science， 2022， 50 （3）：121-124.

[16]劉永杰，王淵，付強，等.高通量測序技術(shù)在病原生物學(xué)方面的研究進展[J].口岸衛(wèi)生控制，2019，24（1）：6-9.

LIU Y J， WANG Y， FU Q， et al. Research progress of high-throughput sequencing technology in pathogenic biology[J]. Port Health Control， 2019， 24（1）：6-9.

[17]吾爾恩·阿合別爾迪，陳雪靜.基于高通量測序技術(shù)馬腸道菌群研究進展[J].微生物學(xué)雜志，2022，42（1）：80-88.

URN A， CHEN X J. Research progress of intestinal microflora of equine based on high-throughput sequencing technology[J]. Chinese Journal of Microbiology， 2022， 42（1）：80-88.

[18]唐濤，袁斌，王帆帆，等.黃連根腐病發(fā)生與根際土壤、根莖內(nèi)生、葉內(nèi)生細菌群落結(jié)構(gòu)變化關(guān)系[J].植物保護，2022，48（1）：73-81.

TANG T， YUAN B， WANG F F， et al. Relationship between occurrence of Rhizoma coptidis root rot and changes of rhizosphere soil， endogenous rhizosphere and endogenous leaf bacterial community structure[J]. Plant Protection， 2002， 48（1）：73-81.

[19]向立剛，郭華，周浩，等.健康與感染青枯病煙株根際土壤與莖稈真菌群落結(jié)構(gòu)與多樣性[J].植物保護，2020，46（1）：189-196.

XIANG L G， GUO H， ZHOU H， et al. Community structure and diversity of rhizosphere soil and stem fungi in healthy and infected tobacco plants[J]. Plant Protection， 2020， 46（1）：189-196.

[20]韓永琴，陳新建，羅路云，等.生防菌劑多黏類芽胞桿菌對辣椒根際土壤細菌群落的影響[J].植物保護，2020，46（2）：135-142.

HAN Y Q， CHEN X J， LUO L Y， et al. Effects of biocontrol agent Bacillus polymyxidae on bacterial community in rhizosphere soil of pepper[J]. Plant Protection， 2020， 46（2）：135-142.

[21]WEI Z， HU J， GU Y， et al. Ralstonia solanacearum pathogen disrupts bacterial rhizosphere microbiome during an invasion[J]. Soil Biology & Biochemistry， 2018， 118：8-17.

[22]馮泉清，付文，陳秀虹，等.盤雙端毛孢屬中國一新紀(jì)錄種[J].大理學(xué)院學(xué)報，2015，14（6）：66-67.

FENG Q Q， FU W， CHEN X H， et al. A New record specie of Seimatosporium from China[J]. Journal of Dali University， 2015， 14（6）：66-67.

[23]商娜，任愛芝，劉冰，等.白三葉草一種葉斑病病原菌的分離鑒定[J].中國農(nóng)學(xué)通報，2022，38（4）：81-85.

SHANG N， REN A Z， LIU B， et al. Isolation and identification of a pathogen of leaf spot disease in white clover[J]. Chinese Agricultural Science Bulletin， 2022， 38（4）：81-85.

[24]楊志英，盧佳欣，李濤，等.貴州開陽地區(qū)茶葉斑病相關(guān)病原菌長柄鏈格孢菌（Alternaria longipes）的研究[J].中國植保導(dǎo)刊，2020，40（11）：5-12.

YANG Z Y， LU J X， LI T， et al. Study on Alternaria longipes， a pathogen related to tea leaf spot disease in Kaiyang， Guizhou Province[J]. China Plant Protection Guide， 2020， 40（11）：5-12.

[25]劉俏，寧楠楠，馬永強，等.青海省櫻桃葉斑病病原菌的分離與鑒定[J].植物保護，2020，46（2）：48-55.

LIU Q， NING N N， MA Y Q， et al. Isolation and identification of the pathogen of cherry leaf spot disease in Qinghai Province[J]. Plant Protection， 2020， 46（2）：48-55.

[26]陳丹.拮抗菌對灰葡萄孢霉菌的抑制作用[D].大連：大連工業(yè)大學(xué)，2019.

CHEN D. Inhibitory effect of antagonistic bacteria on Borealis cinerea[D]. Dalian： Dalian Polytechnic University， 2019.