患擬油壺菌病條斑紫菜表面附生菌群分析*

閻永偉 楊慧超,2 莫照蘭,2,3 李 杰

患擬油壺菌病條斑紫菜表面附生菌群分析*

閻永偉1楊慧超1,2莫照蘭1,2,3李 杰1①

(1. 中國水產(chǎn)科學研究院黃海水產(chǎn)研究所 青島海洋科學與技術試點國家實驗室海洋漁業(yè)科學與食物產(chǎn)出過程功能實驗室 農(nóng)業(yè)農(nóng)村部海水養(yǎng)殖病害防治重點實驗室 山東 青島 266071； 2. 上海海洋大學 水產(chǎn)科學國家級實驗教學示范中心 上海 201306； 3. 中國海洋大學三亞海洋研究院 熱帶海洋生物種質資源開發(fā)與種業(yè)工程實驗室 海南 三亞 572024)

擬油壺菌病(disease)是海上栽培紫菜(sensu lato)的主要病害之一，常引起紫菜大面積病爛并造成嚴重經(jīng)濟損失。本研究利用擬油壺菌感染海區(qū)內不同健康狀態(tài)下的條斑紫菜()[未發(fā)生任何病爛(PyOlpH)、部分紫菜發(fā)生病爛(PyOlpM)和發(fā)生嚴重病爛(PyOlpS)]，分析其附生菌群多樣性、群落結構和主要類群之間的相互作用。結果顯示，3種紫菜附生菌群α多樣性指數(shù)不存在顯著差異，但PyOlpM組指數(shù)高于其他組。3種紫菜附生菌群共有可操作分類單元(operational taxonomic unit, OTU)數(shù)僅占總OTU數(shù)的22.7%，菌群之間存在顯著差異(置換多元方差分析,2=0.405,<0.05)。紫菜感染程度越高，與PyOlpH之間的差異類群數(shù)量越多。共注釋出23門208屬，α-變形菌綱(α-Proteobacteria)、γ-變形菌綱(γ-Proteobacteria)和厚壁菌門(Firmicutes)在所有樣品中均占優(yōu)勢，相對豐度前20個屬中有16個也位于這3個類群中，且隨著感染程度的增加分別出現(xiàn)遞增或遞減的趨勢。其中，貪銅菌屬()和鞘氨醇單胞菌屬()是共現(xiàn)網(wǎng)絡中連接度最高的細菌類群，二者及其依靠正相互作用連接的細菌類群之間存在負相互作用。本研究可為闡明擬油壺菌致病的微生態(tài)機制及尋找生防細菌提供一定的數(shù)據(jù)支持。

條斑紫菜；擬油壺菌病；附生菌群；細菌互作

紫菜(sensu lato，亦稱laver)是一種深受世界人民喜愛的高營養(yǎng)價值海洋紅藻，其廣泛栽培主要集中在中國、韓國和日本等亞洲國家。由于部分海區(qū)的水質貧瘠、潮流不通、水體污染、溫度不適、養(yǎng)殖管理不當?shù)纫蛩?，近年來我國紫菜栽培業(yè)飽受頻發(fā)病害的困擾(黃健等, 2000)，部分栽培海區(qū)的病爛面積可達40%～50% (Yan, 2019; 何禮娟等, 2021)。擬油壺菌病(disease)是海上栽培紫菜的重要病害之一，可在短時間內蔓延并導致大面積病爛(Klochkova, 2016)。已確定該病病原屬于專性細胞內寄生的卵菌綱擬油壺菌屬()，已發(fā)現(xiàn)的、.、.var.、var.、、.以及等7種菌均具有感染紫菜的能力(West, 2006; Sekimoto, 2008、2009; Klochkova, 2016; Kwak, 2017; Badis, 2019)。我國最早于1992年在江蘇省南部栽培的條斑紫菜(，原名[Yang, 2020])中發(fā)現(xiàn)該病(馬家海, 1992)，并先后在栽培條斑紫菜和壇紫菜(, 原名[Yang, 2020])中發(fā)現(xiàn)感染該病(馬家海等, 2007; 劉一萌等, 2012; 何禮娟等, 2021)。擬油壺菌病是對紫菜栽培威脅最大的病害之一，目前對病原尚沒有體外分離培養(yǎng)的方法(Sekimoto, 2008)，因此，對其致病機制還知之甚少，也缺乏對該病害長期有效的防控方法。已經(jīng)施用的措施包括降低栽培密度、網(wǎng)簾干出以及酸洗等，但效果并不理想(Kim, 2014)。因此，探索病害發(fā)生的機理及其生物防控方法十分有必要。

近年來，附生細菌逐漸成為探究藻類病害發(fā)生機理和生物防控的熱點和靶點。這是因為藻類附生細菌除了為藻類生長提供重要的營養(yǎng)物質和生長因子、幫助藻類完成生命周期等重要作用(Goecke, 2010; Singh, 2014)，還與藻類病害密切相關。一方面，環(huán)境脅迫條件下藻類附生微生物的生態(tài)失調是導致病害的重要原因。某些致病微生物類群或其毒力相關的功能基因群在受脅迫個體上大量富集并最終致病(Egan, 2014、2016)。另一方面，藻類附生細菌還能通過分泌抗菌物質或群感效應(quorum sensing, QS)抑制劑等次級代謝產(chǎn)物干擾和抵抗致病菌的附著，保護海藻不受環(huán)境中病原的侵染等(Egan, 2014)。經(jīng)統(tǒng)計，藻類附生細菌中約有12%～50%的菌株具有抗菌活力，是生防菌和抗病物質的重要來源(Singh, 2015)。

此外，對植物和動物卵菌病害的研究也充分體現(xiàn)了附生細菌在促病、防病方面的重要作用。一方面，植物表面附生細菌可產(chǎn)生卵菌毒力必需的細菌毒素或通過共感染以加強卵菌效應基因表達來促進卵菌的致病性(Larousse, 2017)。另一方面，附生細菌能抑制卵菌病原的附著或者利用其基因產(chǎn)物抑制卵菌孢子的萌發(fā)及菌絲的生長。例如，大西洋鮭(L.)的卵受水霉病侵染時，低發(fā)病率的魚卵上附著的高豐度葉居菌屬()細菌能有效抑制病原的附著(Liu, 2014)；而陰溝腸桿菌()則能利用和基因表達的氧化酶和蛋白與終極腐霉()競爭植物產(chǎn)生的不飽和脂肪酸，從而抑制腐霉的萌發(fā)并達到抑病的效果(van Dijk, 2000)。

然而迄今為止，附生細菌與紫菜擬油壺菌病的關系尚未明確，一定程度上阻礙了附生細菌在促病、抗病方面的研究進展。因此，研究紫菜附生細菌與擬油壺菌病之間的關系，對病害的生物防控具有十分重要的意義。本研究分析了受擬油壺菌病感染的不同健康狀態(tài)下的條斑紫菜附生菌群變化，擬為探究病害發(fā)生的微生態(tài)機制及尋找潛在生防細菌類群提供一定的基礎數(shù)據(jù)。

1 材料與方法

1.1 樣品采集

樣品采集地點為條斑紫菜的重要栽培區(qū)域，位于江蘇省鹽城市大豐區(qū)沿岸東沙海區(qū)(121°03′25.88″E, 33°20′11.28″N)，海區(qū)面積約為1333hm2。2018年10月—2019年1月，該海區(qū)條斑紫菜栽培面積共約333 hm2。2019年1月4日開始，約133 hm2區(qū)域條斑紫菜發(fā)生病爛，病爛面積達40%左右。從該區(qū)域選擇3種不同狀況的網(wǎng)簾進行條斑紫菜采樣：網(wǎng)簾上紫菜未發(fā)生任何病爛(PyOlpH)、網(wǎng)簾上部分紫菜已發(fā)生病爛(PyOlpM)以及網(wǎng)簾上紫菜已經(jīng)發(fā)生嚴重病爛(PyOlpS)(圖1)。每種類型網(wǎng)簾均采集3個生物重復，每個生物重復采集同一網(wǎng)簾上相對較為完整的多株(30～40株)條斑紫菜藻體，均勻混合后置于無菌采樣袋中，同時用無菌采水袋采集網(wǎng)簾下部的原位海水以浸沒采樣袋中的紫菜樣品。采集好的紫菜樣品存放于低溫保溫箱中帶回實驗室進行后續(xù)研究。在每個生物重中隨機挑選5～6株條斑紫菜，通過鏡檢確認感染情況。然后，在每個生物重復中再隨機挑選5～6株紫菜藻體，混合后用于附生菌群DNA提取及微生物群落分析。

1.2 病理觀察及病原鑒定

病理觀察及病原鑒定過程見何禮娟等(2021)。將紫菜樣品用無菌海水沖洗3～5次后用滅菌的毛筆去除表面雜質，并用無菌刀片切取面積約為1 cm2的條斑紫菜典型病爛部位，在顯微鏡下觀察后，記錄病爛紫菜的組織病理變化。未發(fā)生病爛網(wǎng)簾的紫菜則隨機切取相同面積的組織進行顯微觀察。

用真菌基因組DNA提取試劑盒(思科捷, 青島)提取PyOlpM和PyOlpS組樣品基因組DNA，并分別用真菌ITS通用引物ITS5和ITS4 (White, 1990)、腐霉()1序列特異引物1-pyth-F和1-pyth-R (Lee, 2015)以及擬油壺菌1特異引物1-olpi-F和1-olpi-R (Kwak, 2017)進行PCR擴增。最后，PCR產(chǎn)物用1%瓊脂糖凝膠電泳進行檢測，并將目的片段膠回收、T載體連接及轉化，將陽性克隆進行測序分析。所得序列在NCBI中進行BLAST同源性比對，并在MEGA-X軟件中用鄰接法(Neighbor-Joining, N-J)構建系統(tǒng)發(fā)育進化樹，模型選擇為Kimura2-paramete (Bootstrap設置為1000)。

圖1 不同健康狀態(tài)下的紫菜栽培網(wǎng)簾比較

1.3 附生菌群DNA提取、PCR擴增及高通量測序

紫菜樣品表面附生細菌群落DNA提取參照Burke等(2009)的方法并做適當修改：每個生物重復的5～6株紫菜藻體先用滅菌海水清洗3次以去除松散結合的浮游微生物；然后，將藻體置于15 mL CMFSW緩沖液(含0.45 mol/L NaCl、10 mmol/L KCl、7 mmol/L Na2SO4和0.5 mmol/L NaHCO3)中，并添加10 mmol/L EDTA和150 μL過濾除菌的快速多酶消化液(3M, 澳大利亞)于室溫下孵育2 h；移除藻體后，清液用等體積的酚∶氯仿∶異戊醇(25∶24∶1)和氯仿∶異戊醇(24∶1)各抽提一次，取上清液并加入3倍體積的無水乙醇和0.1倍體積的醋酸鈉溶液(3 mol/L, pH 5.2)，混合均勻后于?20℃過夜沉淀；4℃20000×離心30 min獲得DNA沉淀后，用70%的乙醇溶液洗滌2次；最后，用200 μL的無酶滅菌水溶解DNA并用瓊脂糖凝膠電泳確定DNA片段大小及質量。

取2 μL DNA，用Q5高保真DNA聚合酶及地球微生物組計劃(Earth Microbiome Project)推薦的16S rRNA基因通用引物515F (5′-GTGYCAGCMGCCGCG GTAA-3′)和926R(5′-CCGYCAATTYMTTTRAGTTT -3′)擴增V4～V5可變區(qū)(Walters, 2016)。每個樣品進行擴增時，515F引物5′端會附加6堿基的barcode以對不同樣品加以區(qū)分。擴增條件：98℃預變性2 min；98℃變性15 s，55℃退火30 s，72℃延伸30 s，共25～27個循環(huán)；72℃后延伸5 min。每個樣品的PCR產(chǎn)物用AxyPrepDNA凝膠回收試劑盒(Axygen, 美國)切膠回收，并用Quant-iT PicoGreen dsDNA Assay Kit (ThermoFisher Scientific, 美國)對PCR產(chǎn)物進行定量。PCR擴增子在上海派森諾公司用Illumina MiSeq平臺進行2×300的雙端測序。

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

每個樣品的雙端下機序列用USEARCH 10.0.240 (Edgar, 2010)進行拼接、質控、去除嵌合體，并按照97%的相似度聚類成可操作分類單元(operational taxonomic unit, OTU)。采用RDP classifier貝葉斯算法在QIIME 1.9.1 (Caporaso, 2010)中基于SILVA132 數(shù)據(jù)庫對OTU代表序列進行分類(cutoff值設為0.8)并生成OTU表，將注釋為Unassignable、Unclassified、Chloroplast以及Mitochondria等的OTU信息從OTU表中去除，所得OTU表用以計算α多樣性、β多樣性和附生菌群組成等。用R (R core team, 2014)和Excel對微生物多樣性及群落進行統(tǒng)計分析，包括各個類群相對豐度統(tǒng)計、α多樣性指數(shù)的單因素方差分析(one-way ANOVA)、韋恩圖分析、主坐標分析(PCoA)及置換多元方差分析(PERMANOVA)。

1.5 差異菌群分析

用LEfSe [line discriminant analysis (LDA) effect size] (Segata, 2011)分析不同感染程度紫菜與未感染紫菜附生菌群中的差異類群。其中，本研究中LEfSe分析的參數(shù)分別設定為LDA score>2.5，=0.05。

1.6 共現(xiàn)網(wǎng)絡分析

用相對豐度前20的細菌屬構建共現(xiàn)網(wǎng)絡(co- occurrence network)。用pysch包中的corr.test ()函數(shù)在R中計算各細菌屬間成對的Spearman相關系數(shù)矩陣，并用Gephi 0.9.2 (Bastian, 2009)對共現(xiàn)網(wǎng)絡進行可視化。網(wǎng)絡圖中的每一個節(jié)點(node)均代表某個細菌屬，節(jié)點之間的連線(也叫作邊, edge)表示2個細菌屬間具有顯著的Spearman相關性(||>0.6, 且FDR校正的<0.01)。

2 結果

2.1 病理觀察及病原鑒定結果

不同狀況網(wǎng)簾上紫菜的顯微觀察結果見圖2。如圖2所示，PyOlpH組紫菜樣品未觀察到任何感染跡象。而PyOlpM和PyOlpS組樣品在顯微鏡下觀察到紫菜細胞受擬油壺菌感染后呈現(xiàn)的典型病理變化特征(紅色箭頭)，且后者病變程度明顯高于前者。詳細的病理觀察結果及分析見何禮娟等(2021)。

利用ITS通用引物、腐霉1特異引物及擬油壺菌1特異引物對患病紫菜DNA進行檢測，發(fā)現(xiàn)僅ITS通用引物和擬油壺菌1特異引物出現(xiàn)擴增條帶(圖3A)，經(jīng)測序及構建N-J進化樹發(fā)現(xiàn)該病原與聚為一支(圖3B)，說明引起大豐地區(qū)紫菜病爛的病原為擬油壺菌。

2.2 各樣本測序結果統(tǒng)計

各樣品16S rRNA基因高通量測序結果見表1。如表1所示，下機序列經(jīng)過拼接、質控等預處理，共獲得有效序列283 911條。經(jīng)過OTU聚類、代表序列注釋以及剔除無效序列等，共獲得247 657條高質量序列，分屬于432個OTU。對測序覆蓋度進行統(tǒng)計，結果顯示，各樣本Good′s coverage指數(shù)均大于0.99 (表1)，表明該測序深度能夠在很大程度上反映各樣品中的附生菌群多樣性。

圖2 不同狀態(tài)網(wǎng)簾上條斑紫菜的顯微觀察

PyOlpS組條斑紫菜顯微觀察結果引自何禮娟等(2021)

Microscopic observation ofin PyOlpS is a quotation from He(2021)

圖3 江蘇省大豐病爛條斑紫菜病原鑒定

A：ITS和1基因擴增的瓊脂糖凝膠電泳圖；B：基于1基因構建的N-J系統(tǒng)進化樹

A: Argarose gel elecgtrophoresis for the amplification of ITS and1; B: N-J phylogenetic tree constructed based on1 gene sequences

M：DL2000 DNA分子量標準；1、2：嚴重和中度病爛紫菜ITS擴增；3、4：嚴重和中度病爛紫菜的腐霉1擴增；5、6：嚴重和中度病爛紫菜的擬油壺菌1擴增；7、8和9分別為ITS、腐霉1和擬油壺菌1陰性對照(無菌水作模板)

M: DL 2000 DNA marker; 1, 2: ITS amplification for PyOlpS and PyOlpM; 3, 4:1 amplification for PyOlpS and PyOlpM; 5, 6:1 amplification for PyOlpS and PyOlpM; 7, 8 and 9: Negative controls using sterilized water as template for ITS,1 and1 amplification, respectively

表1 各樣本測序結果統(tǒng)計

Tab.1 Statistics of the sequencing results of all samples

2.3 α多樣性和β多樣性

用Observed species (Obs)、Chao1、PD whole tree (PD)和Shannon等指數(shù)來表征α多樣性(圖4)。各組之間的多樣性差異均不顯著(>0.05, ANOVA, Tukey′s HSD test)。但是，PyOlpM組樣品的4個α多樣性指數(shù)(分別為133.00、143.00、11.01和5.30)均高于PyOlpH組(分別為122.00、134.00、10.40和5.26)和PyOlpS組(分別為125.00、132.00、9.91和4.41)，而PyOlpS組除了Obs高于PyOlpH組外，其余指數(shù)也均低于PyOlpH組。

根據(jù)各組間附生菌群落結構多樣性的相互關系在OTU水平上繪制了韋恩圖(圖5A)。結果顯示，3組樣品特有的OTU分別為77、86和92個。有98個OTU在3組樣品中均存在，占OTU總數(shù)的22.7%。3組樣品兩兩之間也存在共有OTU，其中，PyOlpH和PyOlpM之間共有31個，PyOlpM和PyOlpS之間共有30個OTU，而PyOlpH和PyOlpS之間僅共有18個OTU。基于Bray-Curtis距離的PCoA顯示，紫菜附生菌群可根據(jù)條斑紫菜被侵染的程度而分開(PCo1解釋度為41.4%)，表明不同狀態(tài)紫菜附生菌群之間存在明顯差異(圖5B)。基于Bray-Curtis距離的PERMANOVA表明，這種組間差異具有顯著性(2= 0.405,=0.042<0.05)。

圖4 不同狀態(tài)紫菜附生菌群的α多樣性指數(shù)比較

圖5 基于OTU的韋恩圖分析(A)和主坐標分析(B)

2.4 微生物群落組成及差異分析

經(jīng)過序列注釋及分類信息統(tǒng)計，所得高質量序列共注釋出23門208屬。其中，變形菌門(Proteobacteria)中的α-變形菌綱(α-Proteobacteria)和γ-變形菌綱(γ-Proteobacteria)以及厚壁菌門(Firmicutes)在所有樣品中均占優(yōu)勢，在PyOlpH、PyOlpM和PyOlpS組中共分別占84.7%、83.4%和88.0% (圖6)。α-變形菌綱(PyOlpH: 26.91%; PyOlpM: 24.59%; PyOlpS: 20.44%)和厚壁菌門(PyOlpH: 17.29%; PyOlpM: 9.56%; PyOlpS: 5.95%)相對豐度隨著感染程度增加呈遞減趨勢，而γ-變形菌綱相對豐度(PyOlpH: 40.5%; PyOlpM: 49.2%; PyOlpS: 61.6%)則隨著感染程度增加呈遞增趨勢。

圖6 不同狀態(tài)條斑紫菜附生菌群在門(或變形菌綱)水平上的相對豐度變化

從細菌屬水平來看，相對豐度前20的細菌屬中共有16個屬分布于α-變形菌綱、γ-變形菌綱和厚壁菌門(圖7)。例如，α-變形菌綱中的甲基桿菌屬()、鞘氨醇單胞菌屬()和柄桿菌屬()，前二者相對豐度隨著感染程度的增加而降低，而后者則相反；γ-變形菌綱中的貪銅菌屬()、假單胞菌屬()和屬的相對豐度均隨感染程度增加而增加；厚壁菌門中的屬的相對豐度隨感染程度增加而減少。

LefSe分析得到的差異菌群結果顯示，PyOlpM與PyOlpH的附生菌群之間共有16個類群存在相對豐度的顯著差異，其中鞘脂單胞菌目(Sphigomonadales)細菌在PyOlpH中顯著富集，而黃單胞菌目(Xanthomonadales)則在PyOlpM中顯著富集(圖8A)。PyOlpS與PyOlpH的附生菌群之間相對豐度存在顯著差異的菌群數(shù)量則達到37個，如PyOlpS中的假單胞菌目(Pseudomonadales)和柄桿菌目(Caulobacterales)，表明有更多的細菌類群參與到擬油壺菌的后期感染中(圖8B)。

2.5 附生細菌的共現(xiàn)網(wǎng)絡分析

利用相對豐度前20的細菌屬構建了一個節(jié)點數(shù)為18、邊數(shù)為27的共現(xiàn)網(wǎng)絡，每個節(jié)點的平均度為3，表示平均每個細菌屬連接3個其他細菌屬，正相互作用的邊數(shù)為14 (51.85%)，負相互作用的邊數(shù)為13 (48.15%)。其中，貪銅菌屬的度為最高(8)，且與柄桿菌屬、假單胞菌屬和屬之間存在正相互作用(圖9)；鞘氨醇單胞菌屬的度次之(7)，該屬細菌與甲基桿菌屬、酸桿菌屬()、蒼白桿菌屬()以及高溫放線菌屬()之間存在正相互作用，且這些細菌屬均與貪銅菌屬細菌或與貪銅菌屬細菌存在正相互作用的細菌類群之間具有負相互作用(圖9)。

3 討論

本研究以江蘇省大豐區(qū)某栽培海區(qū)遭受擬油壺菌感染的條斑紫菜為研究對象，分析了不同感染狀態(tài)下條斑紫菜附生菌群的多樣性和組成差異以及主要主要細菌類群之間的互作關系。

圖7 不同狀態(tài)條斑紫菜附生菌群中前30個細菌屬的相對豐度熱圖

綠色表示相對豐度較高(>1%)的屬，紅色表示相對豐度較低(<1%)的屬

Genus with high relative abundance (>1%) are marked in green, and those with low relative abundance (<1%) are marked in red

圖8 健康紫菜附生菌群與不同感染程度條斑紫菜附生菌之間的差異菌群分析

紅色表示相對豐度在PyOlpH組顯著更高的微生物類群，綠色表示相對豐度在病爛紫菜中顯著更高的類群

Microbial taxa with relative abundance significantly higher in PyOlpH are marked in red, and those with relative abundance significantly higher in diseasedare marked in green

圖9 前20個細菌屬之間的共現(xiàn)網(wǎng)絡

節(jié)點的大小與某一屬連接其他屬數(shù)量的多少成正相關; 邊的粗細與Spearman相關系數(shù)的大小成正比, 其中灰色的邊表示細菌屬之間的正顯著相互作用, 而紅色的邊則表示負顯著相互作用

The size of each node is proportional to the number of connections. The thickness of the edge is proportional to the Spearman’s correlation coefficient. The gray edge indicates significantly positive interaction between bacteria genera, while the red edge indicates significantly negative interaction

研究結果顯示，附生菌群的α多樣性指數(shù)在不同狀態(tài)樣品之間不存在顯著差異。這與赤腐病感染紫菜的菌群多樣性變化較為一致(Yan, 2019)。但兩項研究均發(fā)現(xiàn)，初期感染樣品附生菌群的α多樣性指數(shù)高于未發(fā)生感染樣品，而且β多樣性的結果顯示了不同感染狀態(tài)樣品附生菌群之間存在顯著差異，表明病原的感染與菌群的紊亂密切相關。另外，擬油壺菌病的發(fā)生與栽培環(huán)境的變化緊密相關，如水溫的升高、海水鹽度的降低、水流交換不佳以及干出時間不足等(馬家海, 1992; Klochkova, 2012、2016)。綜合兩方面，與Egan等(2016)提出的觀點較為符合，即環(huán)境脅迫條件下的微生物生態(tài)失調是造成藻類病害發(fā)生的主要原因，因此，藻類病害的發(fā)生發(fā)展是一個復雜的微生態(tài)過程。

研究表明，變形菌門和厚壁菌門通常是海藻表面占主導地位的2個細菌類群，在海藻表面廣泛存在的細菌中它們可被看作是“微生物的雜草種類” (microbial weed species)，并且具有多種多樣的生物功能，例如幫助宿主抵抗多種脅迫環(huán)境因子、具有高效的產(chǎn)能系統(tǒng)以及誘導宿主的形態(tài)發(fā)生和后續(xù)的生長，表明這些類群的細菌對藻類宿主健康的重要性(Singh, 2014)。由于絕大多數(shù)細菌在藻類特別是紫菜葉狀體上的功能還沒有詳盡的解析，因此，只能通過細菌在植物等宿主或者環(huán)境中的相應功能在一定程度上解析其在本研究中的生態(tài)功能。例如，α-變形菌綱中的主要類群(如甲基桿菌屬和鞘氨醇單胞菌屬)以及γ-變形菌綱中的假單胞菌屬是陸生植物的葉際最常檢出的微生物類群。其中，甲基桿菌屬和假單胞菌屬都能分泌AHLs (acyl homoserine lactones)，是植物葉際QS的重要參與者，且后者往往對植物具有致病性，并通過QS分泌胞外聚合物增強其對附著生活的適應性(Vorholt, 2012)。此外，甲基桿菌屬細菌還具有競爭利用植物葉際碳源的功能，因此，它可能通過QS或底物競爭等形式調節(jié)微生物的正常附著和聚集(Vorholt, 2012)。而善于從植物葉際吸收低濃度營養(yǎng)物質的鞘氨醇單胞菌屬細菌則對病原假單胞菌的感染具有強烈的抑制作用，不僅能夠抑制感染癥狀，還能夠減弱病原的生長(Innerebner, 2011)。另外，土壤、水體以及植物等環(huán)境中常見的貪銅菌屬細菌能利用多種碳源進行生長，并且具有抗生素和消毒劑的抗性，甚至可導致免疫力低下的人患上嚴重并發(fā)癥(Balada-Llasat, 2010; Khan, 2016)。這些微生物的生態(tài)功能在一定程度上解釋了這些細菌類群的相對豐度變化規(guī)律，亦即對宿主健康有益的細菌則相對豐度隨著感染程度增加而減少，不利宿主健康的細菌則反之。

除此之外，細菌之間的共現(xiàn)網(wǎng)絡圖也能比較直觀地反映細菌的功能和互作關系。比如，與鞘氨醇單胞菌屬細菌具有正相互作用的細菌類群(如甲基桿菌等，圖9藍色圓圈)往往與假單胞菌和貪銅菌等病害相關的細菌類群(圖9紫色圓圈)具有負相互作用，表明二者之間具有競爭、抑制的作用。這樣的負相互作用為抑病生防菌的查找提供了重要的理論指導(Poudel, 2016)。研究還發(fā)現(xiàn)，PyOlpM比PyOlpH組的附生菌群顯著富集了黃單胞菌目的細菌，這類細菌也是植物葉際表面QS的參與者之一。Berendsen等(2018)研究表明，葉片受霜霉病病原菌侵染的擬南芥在根際招募黃單胞菌目的黃單胞菌屬、寡養(yǎng)單胞菌屬()以及放線菌目(Actinomycetales)的微桿菌屬()等細菌協(xié)同幫助植物抵抗病害。因此可推測，擬油壺菌對紫菜的感染也可使宿主招募某些類群的益生菌來抵抗病害。但是，由于擬油壺菌和霜霉均為內寄生的病原，而由益生菌誘導出的宿主抗病力并不足以治愈已經(jīng)感染的個體(Berendsen, 2018)，因此，紫菜最終仍會被病原及后續(xù)生長起來的其他致病菌侵染。

基于本研究結果可提出以下設想：當環(huán)境條件變得適宜擬油壺菌侵染時，紫菜附生菌群發(fā)生紊亂并朝著利于病原及其他有害微生物附著生長的方向演變(如甲基桿菌、鞘氨醇單胞菌豐度的不斷降低，貪銅菌和假單胞菌豐度的不斷升高)，這使得擬油壺菌病原不斷附著到紫菜上并侵染、殺死宿主細胞；隨著病原的持續(xù)感染，雖然紫菜會招募一些益生菌(如黃單胞菌目細菌)來幫助其抗病，但并不足以使其痊愈；貪銅菌和假單胞菌等有害細菌通過不斷增殖甚至產(chǎn)生抵抗抗生素和消毒劑的能力，導致病害難以控制，紫菜組織持續(xù)病爛并凋亡。

綜上所述，本研究探究了遭受擬油壺菌感染后不同健康狀態(tài)下條斑紫菜的附生菌群多樣性、群落結構和主要類群之間的相互作用，發(fā)現(xiàn)不同健康狀態(tài)下的紫菜附生菌群存在明顯差異，且主要細菌類群之間存在顯著的互作關系，表明附生菌群的變化和互作在病害發(fā)生過程中起到了重要的作用。這些結果為闡明擬油壺菌致病的微生態(tài)機制及尋找生防細菌提供了一定的數(shù)據(jù)支持，今后仍需要大量實驗驗證發(fā)病的生態(tài)機制及附生細菌的生物功能。

BADIS Y, KLOCHKOVA T A, STRITTMATTER M,. Novel species of the oomycetepotentially threaten European red algal cultivation. Journal of Applied Phycology, 2019, 31(2): 1239–1250

BALADA-LLASAT J M, ELKINS C, SWYERS L,. Pseudo-outbreak ofinfection at an outpatient clinic related to rinsing culturette swabs in tap water. Journal of Clinical Microbiology, 2010, 48(7): 2645– 2647

BASTIAN M, HEYMANN S, JACOMY M. Gephi: An open source software for exploring and manipulating networks. Dublin: Third International ICWSM Conference, 2009, 361–362

BERENDSEN R L, VISMANS G, YU K,. Disease-induced assemblage of a plant-beneficial bacterial consortium. ISME Journal, 2018, 12(6): 1496–1507

BURKE C, STAFFAN K, TORSTEN T. Selective extraction of bacterial DNA from the surfaces of macroalgae. Applied and Environmental Microbiology, 2009, 75(1): 252–256

CAPORASO J G, KUCZYNSKI J, STOMBAUGH J,. QIIME allows analysis of high-throughput community sequencing data. Nature Methods, 2010, 7(5): 335–336

EDGAR R C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics, 2010, 26(19): 2460–2461

EGAN S, FERNANDES N D, KUMAR V,. Bacterial pathogens, virulence mechanism and host defence in marine macroalgae. Environmental Microbiology, 2014, 16(4): 925–938

EGAN S, GARDINER M. Microbial dysbiosis: Rethinking disease in marine ecosystems. Frontiers in Microbiology, 2016, 7: 991

GOECKE F, LABES A, WIESE J,. Chemical interactions between marine macroalgae and bacteria. Marine Ecology Progress Series, 2010, 409: 267–299

HE L J, YANG H C, LI J,. Diagnosis ofdisease inProgress in Fishery Sciences, 2021, 42(4): 177–183 [何禮娟, 楊慧超, 李杰, 等. 一起條斑紫菜擬油壺菌病的診斷. 漁業(yè)科學進展, 2021, 42(4): 177– 183]

HUANG J, TANG X X, YAN X J,. Review and forecast of the study on the disease of kelp and laver in China. Marine Science, 2000, 24(9): 26–27 [黃健, 唐學璽, 嚴小軍, 等. 我國海帶和紫菜病害的研究展望. 海洋科學, 2000, 24(9): 26–27]

INNEREBNER G, KNIEF C, VORHOLT J A. Protection ofagainst leaf-pathogenicbystrains in a controlled model system. Applied and Environmental Microbiology, 2011, 77(10): 3202–3210

KHAN S, BEATTIE T K, KNAPP C W,. Relationship between antibiotic- and disinfectant- resistance profiles in bacteria harvested from tap water. Chemosphere, 2016, 152: 132–141

KIM G H, MOON K H, KIM J Y,. A revaluation of algal diseases in Korean() sea farms and their economic impact. Algae, 2014, 29(4): 249–265

KLOCHKOVA T A, SHIM J B, HWANG M S,. Host-parasite interactions and host species susceptibility of the marine oomycete parasite,sp., from Korea that infects red algae. Journal of Applied Phycology, 2012, 24(1): 135–144

KLOCHKOVA T A, SHIN Y J, MOON K H,. New species of unicellular obligate parasite,sp. nov., that plaguessea farms in Korea. Journal of Applied Phycology, 2016, 28(1): 73–83

KWAK M S, KLOCHKOVA T A, JEONG S,.var., an endemic endoparasite infecting cultivatedin Korea. Journal of Applied Phycology, 2017, 29(7): 1–10

LAROUSSE M, GALIANA E. Microbial partnerships of pathogenic Oomycetes. PLoS Pathogens, 2017, 13(1): e1006028

LEE S J, HWANG M S, PARK M A,. Molecular identification of the algal pathogen(Oomycetes) from(Rhodophyta) using ITS and1 markers. Algae, 2015, 30(3): 217–222

LIU Y M, MA J H, WEN Q. Simultaneous infection by red rot and olpidiops disease in laver. Journal of Dalian Ocean University, 2012, 27(6): 546–550 [劉一萌, 馬家海, 文茜. 壇紫菜赤腐病與擬油壺菌病并發(fā)病的初步研究. 大連海洋大學學報, 2012, 27(6): 546–550]

LIU Y, DE BRUIJN I, JACK A L,. Deciphering microbial landscapes of fish eggs to mitigate emerging diseases. ISME Journal, 2014, 8(10): 2002–2014

MA J H, LIN Q S, MIN J,. Preliminary study on the olpidiops-disease of. Journal of Fisheries of China, 2007, 31(6): 860–864 [馬家海, 林秋生, 閔建, 等. 條斑紫菜擬油壺菌病的初步研究. 水產(chǎn)學報, 2007, 31(6): 860–864]

MA J H. An investigation of chytrid-disease inin the coastal water of south Jiangsu. Journal of Shanghai Fisheries University, 1992, 1(3/4): 185–188 [馬家海. 江蘇省南部沿海條斑紫菜壺狀菌病的調查研究. 上海水產(chǎn)大學學報, 1992, 1(3/4): 185–188]

POUDEL R, JUMPPONEN A, SCHLATTER D C,. Microbiome networks: A systems framework for identifying candidate microbial assemblages for disease management. Phytopathology, 2016, 106(10): 1083–1096

R CORE Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing, 2014

SEGATA N, IZARD J, WALDRON L,. Metagenomic biomarker discovery and explanation. Genome Biology, 2011, 12(6): R60

SEKIMOTO S, KLOCHKOVA T A, WEST J A,.sp. nov.: An endoparasitic oomycete that infectsand other red algae (Rhodophyta). Phycologia, 2009, 48(6): 460–472

SEKIMOTO S, YOKOO K, KAWAMURA Y,. Taxonomy, molecular phylogeny, and ultrastructural morphology ofsp. nov. (), a unicellular obligate endoparasite ofandspp. (). Mycological Research, 2008, 112(3): 361–374

SINGH R P, KUMARI P, REDDY C R K. Antimicrobial compounds from seaweeds-associated bacteria and fungi. Applied Microbiology and Biotechnology, 2015, 99(4): 1571–1586

SINGH R P, REDDY C R K. Seaweed-microbial interactions: Key functions of seaweed-associated bacteria. FEMS Microbiology Ecology, 2014, 88(2): 213–230

VAN DIJK K, NELSON E B. Fatty acid competition as a mechanism by whichsuppressessporangium germination and damping-off. Applied and Environmental Microbiology, 2000, 66(12): 5340–5347

VORHOLT J A. Microbial life in the phyllosphere. Nature Reviews Microbiology, 2012, 10(12): 828–840

WALTERS W, HYDE E R, BERG-LYONS D,. Improved bacterial 16S rRNA gene (V4 and V4-5) and fungal internal transcribed spacer marker gene primers for microbial community surveys. mSystems, 2016, 1(1): e00009–15

WEST J A, KLOCHKOVA T A, KIM G H,.sp., an oomycete from Madagascar that infectsand other red algae: Host species susceptibility. Phycological Research, 2006, 54(1): 72–85

WHITE T J, BRUNS T, LEE S,. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: INNIS M A, GELFAND D H, SNINSKY J J,. PCR Protocols–A guide to methods and applications. San Diego: Academic Press, 1990, 315–322

YAN Y W, YANG H C, TANG L,. Compositional shifts of bacterial communities associated withand surrounding seawater co-occurring with red rot disease. Frontiers in Microbiology, 2019, 10: 1666

YANG L E, DENG Y Y, XU G P,. Redefining(Bangiales, Rhodophyta): four new genera, resurrection ofand description ofsp. nov. from China. Journal of Phycology, 2020, 56(4): 862–879

Epiphytic Microbial Communities Associated withwithDisease

YAN Yongwei1, YANG Huichao1,2, MO Zhaolan1,2,3, LI Jie1①

(1. Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Laboratory for Marine Fisheries Science and Food Production Processes, Pilot National Laboratory for Marine Science and Technology (Qingdao), Key Laboratory of Maricultural Organism Disease Control, Ministry of Agriculture and Rural Affairs, Qingdao, Shandong 266071, China; 2. National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai 201306, China; 3. Laboratory of Tropical Marine Germplasm Resources and Breeding Engineering, Sanya Oceanographic Institution, Ocean University of China, Sanya, Hainan 572024, China)

disease is one of the main diseases affecting laver cultivation in the sea, often causing large-scale rot and serious economic losses. Although the oomycetic genushas been identified as the pathogen, little is known regarding its mechanisms, including its interaction with epiphytic microbial communities. In the present study, diversities, structures, and major interactions of epiphytic microbial communities were analyzed usingin sea infected by. Theused in the experiment had different health statuses, including no infection (PyOlpH), partial infection (PyOlpM), and serious infection (PyOlpS). The results showed that no significant difference was observed for the α diversity indices among microbial communities associated with the three groups of, but those of PyOlpM were higher than the other two. Shared operational taxonomic units (OTUs) among the three groups of communities only accounted for 22.7% of the total OTUs, in addition to significant community dissimilarity (permutational multivariate analysis of variance,2=0.405,<0.05). More microbial taxa were obtained when communities associated with more seriously infectedcompared to those associated with PyOlpH. In total, 208 genera belonging to 23 phyla were annotated, and α-Proteobacteria, γ-Proteobacteria, and Firmicutes were dominant in all samples. Sixteen of the top 20 genera were also assigned to the three groups, with their relative abundances increasing or decreasing across the infection status from none to severe. Among them, generaandwere the most connected bacterial taxa in the co-occurrence network, and negative interactions were determined between the two taxa and their positively interacting bacterial taxa. This study provides primary datasets for clarifying the microbial ecological mechanisms and searching for probiotics fordisease.

;disease; Epiphytic microbial communities; Bacterial interaction

S917.1

A

2095-9869(2022)03-0165-11

10.19663/j.issn2095-9869.20210326002

http://www.yykxjz.cn/

閻永偉, 楊慧超, 莫照蘭, 李杰. 患擬油壺菌病條斑紫菜表面附生菌群分析. 漁業(yè)科學進展, 2022, 43(3): 165–175

YAN Y W, YANG H C, MO Z L, LI J. Epiphytic microbial communities associated withwithdisease. Progress in Fishery Sciences, 2022, 43(3): 165–175

LI Jie, E-mail: lijie@ysfri.ac.cn

* 財政部和農(nóng)業(yè)農(nóng)村部: 國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術體系資助 [This work was supported by China Agriculture Research System of MOF and MARA]. 閻永偉，E-mail: yanyw@ysfri.ac.cn

李 杰，副研究員，E-mail: lijie@ysfri.ac.cn

2021-03-26,

2021-04-19

(編輯 馬璀艷)