枯草芽孢桿菌M-23對柑橘砂皮病防效及柑橘葉際細菌群落多樣性的影響

李審微 洪艷云 李新文 何可佳 戴良英 盧曉鵬 宋娜 易圖永

摘要：【目的】明確柑橘砂皮病不同發病等級及噴施生防菌枯草芽孢桿菌M-23后柑橘葉際細菌群落結構的變化，揭示生防菌防治柑橘砂皮病的葉際微生物學機制，為進一步研究柑橘砂皮病的生物防治提供理論基礎?！痉椒ā吭囼炘O對照區和噴施生防菌區2個區域，噴施生防菌區分別在2019年3月中旬柑橘春梢長2～3 cm、4月中旬、謝花2/3、果實蠶豆大小和果實乒乓球大小時噴施M-23菌株發酵液，噴藥液量3 L/株;同時對照區噴施等量的滅菌NB液體培養基，在最后一次噴施生防菌15 d 后進行防效調查。采集噴施生防菌區柑橘葉片（R4處理，病情等級0級）和未噴施生防菌區不同發病等級的柑橘葉片（R1、R2和R3處理，病情等級分別為0、5和9級），使用FastDNA[?] Spin Kit for Soil試劑盒提取柑橘葉片葉際微生物總DNA，通過Illumina MiseqTM高通量平臺測序分析微生物群落多樣性?！窘Y果】生防菌M-23對柑橘砂皮病的防治效果達74.02%。高通量測序共獲得571091條序列，單一樣品序列數17833～96147條，在97%的相似水平下進行生物信息統計分析聚類后獲得2932個OTUs。4個處理的OTUs數量分別為1944、1545、1205和607個，表現為R1>R2>R3>R4。在門水平的組成上主要由變形菌門（Proteobacteria）、厚壁菌門（Firmicutes）、擬桿菌門（Bacteroidetes）、放線菌門（Actinobacteria）、浮霉菌門（Planctomycetes）、梭桿菌門（Fusobacteria）、裝甲菌門（Armatimonadetes）、酸桿菌門（Acidobacteria）和疣微菌門（Verrucomicrobia）組成。變形菌門是優勢菌門，在R1、R2、R3和R4處理中的相對豐度分別為81.84%、87.93%、86.80%和68.82%。對照區的葉際微生物多樣性與發病等級成反比;噴施生防菌后柑橘葉片葉際細菌群落多樣性降低，在屬水平上芽孢桿菌屬（Bacillus）和假單胞菌屬（Pseudomonas）的相對豐度得到明顯增加，芽孢桿菌屬的相對豐度達33.41%?！窘Y論】噴施生防菌M-23能有效防治柑橘砂皮病。M-23菌株能在柑橘葉片上定殖和增殖，改善葉際細菌群落結構，提升有益菌的比例，使病情指數降低，具有一定的應用潛力。

關鍵詞： 柑橘砂皮病;葉際細菌;生防菌;高通量測序;多樣性

中圖分類號： S436.66? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼： A 文章編號：2095-1191（2020）07-1699-07

Abstract：【Objective】Focusing on different severities of Diaporthe citri and the shifts of bacterial community diversity in citrus phyllosphere after spraying with Bacillus subtilis M-23， this study was to reveal the phyllospheric microbiological mechanism of controlling D. citri with biocontrol bacteria and provide theoretical basis for further studies on biological control of D. citri. 【Method】The experiment was conducted in two areas—the control area and the treatment area. In the treatment area， biocontrol bacteria M-23 was sprayed in the middle of March， 2019， when the citrus spring tip grew 2-3 cm， in the middle of April，? 2/3 of flowers fell， the citrus size was like broad bean and the citrus size was like ping-pong ball， evenly onto citrus leaves with an inoculum of 3 L/plant. At the same time， the control area was sprayed with the same amount of sterilized NB medium. The control effect was investigated 15 d after the last application of biocontrol agents. Collecting leaves with different incidence grades in the sprayed areas（R4 treatment， incidence grade 0） and unsprayed areas（R1， R2 and R3 treatments，incidence grades 0， 5 and 9 respectively）. FastDNA[?] Spin Kit for Soil kit was used to extract the total DNA of phyllosphere microorganisms， and the diversity of? microbial community was analyzed by Illumina MiseqTM high throughput platform sequencing method. 【Result】The result showed that the control effect of biocontrol bacterium M-23 on D. citri was 74.02%. A total of 571091 high quality sequences were obtained by high throughput sequencing， the number of single sample sequences ranged from 17833 to 96147， and 2932 OTUs were obtained by bioinformatics statistical analysis clustering at a similar level of 97%. The number of OTUs of the four treatments were 1944， 1545， 1205 and 607， respectively，in the following order：R1>R2>R3>R4. The horizontal composition of the phylum was mainly composed of Proteobacteria， Firmicutes， Bacteroidetes， Actinomycetes， Fusobacteria， Armatimonadetes， Acidobacteria and Verruciformis. Proteobacteria were dominant， and the relative abundances in R1， R2， R3 and R4 treatments were 81.84%， 87.93%， 86.80% and 68.82%， respectively. In the control area，phyllosphere microbial diversity was inversely proportio-nal to the incidence level. After spraying the biocontrol agent， the diversity of the bacterial community in the phyllosphere decreased， and the relative abundance of Bacillus and Pseudomonas increased obviously at the genus level， the relative abundance of Bacillus reached 33.41%. 【Conclusion】To sum up， spraying M-23 can effectively prevent the occurrence of D. citri. M-23 can colonize and proliferate on citrus leaves， which improves the community structure of phyllosphere bacteria， increases the proportion of beneficial bacteria， and reduces the disease index. Therefore， B. subtilis M-23 has certain application potential.

Key words： Diaporthe citri; phyllosphere bacteria; biocontrol bacteria; high-throughput sequencing; diversity

Foundation item： Post Scientist Project of Modern Agriculture（Citrus） Industrial Technology System（CARS-26）; Citrus Plant Protection Post Expert of Hunan Fruit Industry Technology System（Xiangnongfa〔2015〕137）; Hunan Agricultural Universitys “Double First-Class” Construction Sub-project（SYL2019027）; Hunan Agricultural University Youth Science Foundation Project（19QN30）

0 引言

【研究意義】近年來，柑橘砂皮病在湖南省大面積發生，且日益嚴重，給湖南省柑橘產業帶來巨大經濟損失（馬靖艷等，2019）。柑橘砂皮病菌（Diaporthe citri）為弱寄生菌，在樹體受傷或樹勢衰弱時進行侵染，特別是侵染葉片和果實，常在表面形成黑色麻點，嚴重影響果實外觀和品質，影響其經濟價值（方天露，2017）。目前，防治柑橘砂皮病主要依靠化學藥劑防治，每年需進行6次以上噴施才有較好的控制效果。大量化學藥劑的使用極易導致病原菌產生抗性和農藥殘留，影響環境微生物的平衡（周娜等，2015）。生物防治措施因對環境友好、安全程度高、不會造成藥物殘留而受到研究者們的重視。因此，研究柑橘砂皮病不同發病等級柑橘葉片葉際微生物多樣性，探究噴施生防菌對葉際微生物多樣性的影響，有助于從微生態分析中發現柑橘砂皮病病害與微生態之間的關系，明確芽孢桿菌對柑橘砂皮病的抗病機制，為防治柑橘砂皮病提供新思路?！厩叭搜芯窟M展】高爽等（2016）綜述了植物葉際微生物群落的組成、特點及與外界環境的互作，發現不同因素均會對葉際微生物群落組成造成一定影響，其變化有規律可循，微生物群落組成的改變會影響寄主植物的生長?；瘜W農藥的使用對植物葉際微生物群落的影響較大，謝蘭芬等（2018）研究發現，使用百菌清處理玉米葉片對其葉際微生物的影響較大，而解淀粉芽孢桿菌與噴施環境兼容性良好，是一種生態、環境友好的防治手段。近年來，隨著高通量測序技術的發展，擴增子測序分析技術逐漸在植物葉際微生物多樣性的研究中得到廣泛使用（Bodenhausen et al.，2013）。羅路云等（2017）比較了不同發病等級下植物葉際微生物種群的多樣性，發現由于病害的發生，植物葉際細菌種群會發生顯著變化，發病葉片葉際細菌多樣性明顯低于健康葉片葉際細菌多樣性。Zhang等（2008）研究蘇云金芽孢桿菌對微生物群落的影響，發現接種蘇云金芽孢桿菌后辣椒葉際微生物結構受到影響，微生物多樣性降低。岑?。?016）將多粘類芽孢桿菌噴施在茶葉上，發現多粘類芽孢桿菌可在葉際定殖，施用生防菌后在實現良好促生和病害抑制作用的基礎上減少了葉際細菌群落的相對豐度，但對于群落的組成和體系的穩定影響較小?！颈狙芯壳腥朦c】目前，針對柑橘砂皮病的研究主要集中在其發生規律和化學防治上，生物防治也僅限于基礎的生防菌篩選，沒有對柑橘砂皮病不同發病等級葉際細菌群落多樣性的差異性比較和施用生防菌發酵液后葉際細菌群落結構變化方面的研究。生防菌M-23由湖南農業大學植物保護學院園藝作物病害研究室篩選獲得，通過前期研究發現其能抑制柑橘砂皮病病原菌菌絲生長和孢子萌發，但未對其田間防治效果及對柑橘葉際微生物的影響進行研究?！緮M解決的關鍵問題】通過生防菌M-23對柑橘砂皮病防治的田間試驗，利用高通量測序技術比較不同處理下柑橘葉際細菌群落多樣性的差異，探究生防菌防治柑橘砂皮病的機制，為柑橘砂皮病的生物防治提供理論參考。

1 材料與方法

1. 1 試驗材料

供試生防菌株為枯草芽孢桿菌（Bacillus subti-lis）M-23，由湖南農業大學植物保護學院園藝作物病害研究室分離獲得。

1. 2 試驗方法

1. 2. 1 試驗地概況 試驗地位于湖南省懷化市洪江市雙溪鎮柑橘園藝場，該橘園于1986年建園，品種為冰糖橙，柑橘砂皮病發生較重。

1. 2. 2 試驗設計 試驗設對照區和噴施生防菌區2個區域。將生防菌M-23活化后以0.8%的接種量接種于滅菌的NB液體培養基中，28 ℃、180 r/min培養5 d后備用。噴施生防菌區分別在2019年3月中旬柑橘春梢長2～3 cm、4月中旬、謝花2/3、果實蠶豆大小和果實乒乓球大小時噴施生防菌發酵液，使用人工背負式電動噴霧器對葉片進行噴灑，噴藥液量3 L/株;同時對照區以相同條件噴施等量的滅菌NB液體培養基。葉片采集設4個處理，其中對照區葉片采集根據柑橘砂皮病的病情等級（R1：0級，無病斑;R2：5級，每葉病斑51～75個;R3：9級，每葉病斑100個以上）對應設R1、R2和R3處理，噴施生防菌區設R4處理（葉片無病斑）。

1. 2. 3 樣品采集方法 最后一次噴施生防菌15 d后，從對照區選擇3株柑橘樹分別采集R1、R2和R3處理的柑橘葉片。葉片采集采用五點取樣法，每株柑橘樹于東西南北中各選擇1個枝條，每枝條采集R1處理同等大小柑橘葉片1片，每處理共15片葉混合放入無菌采樣袋中，進行3次重復，共獲得3個樣品，分別命名為SP01、SP02和SP03;以同樣的方式采集R2和R3處理柑橘葉片樣品，R2處理的樣品命名為SP51、SP52和SP53，R3處理的樣品命名為SP91、SP92和SP93。從噴施生防菌區3棵不同柑橘樹上采集無病的柑橘葉片，方法同上，樣品命名為BS01、BS02和BS03。

1. 2. 4 防治效果調查 參照張岳等（2018）的方法對柑橘砂皮病進行分級：0級，無病;1級，葉片或果實上有1～25個病斑;3級：葉片或果實上有26～50個病斑;5級：葉片或果實上有51～75個病斑;7級：葉片或果實上有76～100個病斑;9級：葉片或果實上有100個病斑以上。收獲期間分別調查對照區和噴施生防菌區柑橘砂皮病病情指數并計算防治效果，每區調查3株柑橘樹，每株5個枝條，共調查15個枝條上葉片的發病情況。

病情指數=∑（各級病葉數×各級代表值）/（調查

總葉數×最高級代表值）×100

防治效果（%）=（對照區病情指數-防治區病情

指數）/對照區病情指數×100

1. 2. 5 葉際微生物總DNA提取 取葉片6 g加入滅菌的錐形瓶中，向錐形瓶加入pH 7的0.1 mol/L PBS緩沖液100 mL、吐溫-80 100 ?L。錐形瓶置于超聲波振蕩器中振蕩10 min后，在恒溫搖床中以28 ℃、180 r/min的條件培養30 min，然后再次用超聲波振蕩器振蕩5 min，在無菌條件下用鑷子取出葉片。采用真空泵對錐形瓶中的液體進行抽濾，用0.45 ?m微孔濾膜收集微生物。使用FastDNA[?] Spin Kit for Soil試劑盒提取葉際微生物DNA，要求DNA濃度和純度的A260/A280為1.8～2.0。

1. 2. 6 16S rDNA的V3～V4區PCR擴增及測序 第一輪擴增，PCR引物為341F[5'-CCCTACACGACGCT CTTCCGATCTG（barcode）CCTACGGGNGGCWGCA G-3']和805R（5'-GACTGGAGTTCCTTGGCACCCG AGAATTCCAGACTACHVGGGTATCTAATCC-3'）。反應體系30 ?L：2×Taq Master Mix 15 ?L，上、下游引物（10 μmol/L）各1 ?L，DNA模板10～20 ng，ddH2O補足至30 ?L。擴增程序：94 ℃預變性3 min;94 ℃ 30 s，45 ℃ 20 s，65 ℃ 30 s，進行5個循環;94 ℃預變性20 s;94 ℃ 30 s，55 ℃ 20 s，72 ℃ 30 s，進行20個循環;72 ℃延伸5 min，4 ℃保存。第二輪擴增，引物Illumina橋式PCR兼容引物。反應體系30 ?L：2×Taq Master Mix 15 ?L，上、下游引物（10 μmol/L）各1 ?L，PCR產物20 ng（第一輪），ddH2O補足至30 ?L。擴增程序：95 ℃預變性3 min;94 ℃ 20 s，55 ℃ 20 s，72 ℃ 30 s，進行5個循環;72 ℃延伸5 min。PCR產物經1%瓊脂糖凝膠電泳檢測后送至生工生物工程（上海）股份有限公司進行高通量測序。

1. 3 數據處理與分析

Illumina MiseqTM平臺獲得的原始數據，通過質控過濾后剔除嵌合體和非特異性擴增序列，得到的有效數據使用Usearc軟件以97%的默認相似性將序列聚類成為一個OTU。以聚類similarity值（0.86～0.99）為橫坐標、OTU數目變化（0～25000）為縱坐標制圖，本研究OTU分析和分析學分析選擇的similarity值為0.97。使用97%相似度的OTU，用Mothur進行Rarefaction分析（Asnicar et al.，2015），利用R制作稀釋曲線用來評估各樣品在測序強度間的差異。根據每個樣品文庫的OTU豐度信息，完成Alpha和Beta多樣性分析。計算Chao1、ACE、Shannon和Simpson指數，評估序列文庫的多樣性。通過NMDS（Non-metric multidimensional scaling）分析4組不同處理微生物群落結構間的差異。

2 結果與分析

2. 1 生防菌M-23對柑橘砂皮病的防治效果

由表1可知，對照區發病較嚴重，病情指數高達36.45;噴施生防菌區施菌后病情指數降低至9.47，防治效果明顯，防效達74.02%。說明生防菌M-23對柑橘砂皮病具有較強的拮抗作用。

2. 2 高通量測序數據預處理結果

利用高通量測序共得到571091條有效序列，單一樣品序列數17833～96147條，平均為47591條，對序列在97%相似水平下進行生物信息統計分析，獲得2932個OTUs。采用稀釋曲線來比較樣品間物種的豐度和試驗數據的合理性，如圖1所示，Shannon指數的稀釋曲線趨于平緩，可認為測序深度已基本覆蓋樣品中絕大多數物種，更多的數據量只會產生少量新的OTU，測序數據量足以反映樣品中的物種多樣性。

2. 3 OTU聚類分析結果

由圖2可看出，R1、R2、R3和R4處理分別有1944、1545、1205和607個OTUs，4個處理間的OTUs數量表現為R1>R2>R3>R4，說明未噴施生防菌且未發病的R1處理微生物種類最豐富，隨著病情指數的增加微生物種類豐富度隨之降低，噴施生防菌的R4處理微生物種類最少。

2. 4 Alpha多樣性指數分析結果

目前，普遍認可種（Species）的分類水平為97%序列相似度（Gu et al.，2019），在該相似度水平下使用獲取的OTU分析各樣品的Alpha多樣性，通過Chao1、ACE、Shannon和Simpson指數對柑橘葉際微生物多樣性進行評估，所得的參數如表2所示。由表2可知，Shannon、Simpson、Chao1和ACE指數的變化范圍分別為1.41～3.77、0.084～0.372、365.6～1256.5和413.7～1302.2。對噴施滅菌培養基的R1、R2和R3處理的Shannon指數進行比較，表現為R1>R2>R3，表明微生物多樣性隨著病情指數的升高而降低;噴施生防菌的R4處理與噴施滅菌培養基的R1、R2和R3處理相比，噴施生防菌后Simpson指數顯著降低（P<0.05），與OTU聚類分析結果吻合。

2. 5 NMDS非度量多維尺度分析結果

根據非度量多維尺度分析結果（圖3），R2和R3處理與R1處理在坐標軸上有一定距離，表明柑橘砂皮病的發生導致葉際細菌群落結構改變;R4處理與R1、R2和R3處理距離較遠，說明生防菌的加入改變了葉際細菌的群落結構。如圖3所示，相同處理的重復能緊密靠在一起，不同處理間可明顯區分，各處理間有明顯差異性，相同處理間群落相似度較高，未噴施生防菌不同病情等級條件下的葉際細菌群落與生防菌處理的細菌群落可明顯區分開來。使用非參數檢驗方法Anosim做進一步分析，根據分析結果（表3）可知，R接近1，說明4個處理的組間差異大于組內差異，P=0.011，小于0.05，說明分析結果具有統計學上的顯著性差異。

2. 6 高通量測序分析柑橘葉際細菌群落結構

柑橘葉際細菌群落門水平和屬水平在各分組樣品中相對豐度≥1%的組成情況見圖4和圖5。在門水平上主要由變形菌門（Proteobacteria）、厚壁菌門（Firmicutes）、擬桿菌門（Bacteroidetes）、放線菌門（Actinobacteria）、浮霉菌門（Planctomycetes）、梭桿菌門（Fusobacteria）、裝甲菌門（Armatimonadetes）、酸桿菌門（Acidobacteria）和疣微菌門（Verrucomicrobia）組成（圖4）。其中變形菌門占絕對優勢，在R1、R2、R3和R4處理中的相對豐度分別為81.84%、87.93%、86.80%和68.82%，噴施滅菌培養基處理的R1處理與R2和R3處理相比，變形菌門的相對豐度隨著病情等級升高而升高，R2處理與R3處理間無明顯差異，噴施生防菌后R4處理中變形菌門的相對豐度有所下降;發病后降低了浮霉菌門的相對豐度;噴施滅菌培養基處理組的厚壁菌門的相對豐度均在10.00%以下，而噴施生防菌后R4處理組厚壁菌門的相對豐度為34.89%，相對豐度明顯增加。

進一步對各分組樣品在屬水平（OTU數相對豐度≥1%）的相對豐度進行研究，結果（圖5）表明，不同樣品的主要葉際細菌群落結構存在明顯差異，其中，假單胞菌屬（Pseudomonas）、鞘氨醇單胞菌屬（Sphingomonas）、甲基桿菌屬（Methylobacterium）和馬賽菌屬（Massilia）為R1、R2和R3處理的共有優勢菌屬，黃單胞菌屬（Xanthomonas）為R2和R3處理的共有優勢菌屬。橘園發病后假單胞菌屬的豐度降低，噴施生防菌后明顯上升，在R4處理組中達最高（42.79%）。鞘氨醇單胞菌屬和黃單胞菌屬在發病后相對豐度先明顯升高，然后略微下降，可能與致病菌存在協同作用。噴施生防菌后R4處理的細菌主要由假單胞菌屬和芽孢桿菌屬（Bacillus）組成，芽孢桿菌屬在R4處理中的相對豐度達33.41%，說明M-23菌株在葉際得到較好的定殖。

3 討論

近年來，湖南省懷化地區由于降水增加、溫度上升、果農對柑橘園管理不當、化學農藥濫用等導致柑橘砂皮病發生越來越嚴重，給柑橘的品質帶來嚴重影響（馬靖艷等，2019），造成重大經濟損失。生物防治是防治柑橘砂皮病最有效的綠色防治措施，通常是在植物生長的自然環境中獲得對病原菌有拮抗作用的生防菌株來防治病害，這種生物防治的方法可維持葉際微生物生態系統的平衡（Wang et al.，2019）。

植物病害的發生會影響其環境微生物結構變化，多樣性和均勻性隨著病害的發生而降低，化學農藥的噴施會影響環境微生物的生態平衡，而生防菌可優化環境微生物的組成（李春宏等，2019;李忠奎等，2019）。岑?。?016）研究多粘類芽孢桿菌在茶葉上的定殖及其對葉際細菌群落的影響時發現，生防菌在葉片上定殖作為一種生物刺激，可改良葉際細菌的生態環境，并不會產生生物入侵的作用打破葉際微生物的生態平衡造成生態環境的紊亂。謝蘭芬等（2018）在生防菌對玉米葉際微生物影響的研究中證明殺菌劑嚴重降低了玉米葉際微生物的代謝活動，而解淀粉芽孢桿菌與噴施環境兼容性良好。Zhang等（2008）探究黃瓜和辣椒的葉際使用蘇云金芽孢桿菌對其葉際微生物結構的影響，結果顯示蘇云芽孢桿菌會降低微生物的多樣性。本研究中R1處理（發病等級為0級）的Chao1、Shannon和ACE指數均大于其他處理，Simpson指數小于其他處理，在病害入侵后其葉際微生物多樣性和豐富度較健康葉際微生物有所降低;噴施生防菌的R4處理Chao1、Shannon和ACE指數的降低可能是因為柑橘葉片噴施生防菌后，突然有大量生防菌加入，其他葉際微生物的定殖受到一定影響，從而提升了有協同作用菌群的相對豐度，抑制了部分菌群的定殖，因此導致葉際微生物的多樣性有所降低，與Zhang等（2008）的研究結果一致。

在葉際細菌菌群結構的組成中，變形菌門、厚壁菌門、擬桿菌門和放線菌門是4個處理共有的菌門，變形菌門和厚壁菌門是其優勢菌門，厚壁菌門在R4處理中所占的比例最高，為34.89%，明顯高于未噴施生防菌的R1、R2和R3處理，其原因是噴施生防菌后枯草芽孢桿菌M-23在葉際得到良好定殖，增加了生防菌株所屬門的相對豐度;變形菌門主要包括α-變形菌綱、β-變形菌綱和γ-變形菌綱，盡管變形菌門在不同物種葉際的相對豐度不一致，但均為葉際細菌中的優勢門類，與Laforest-Lapointe和Whitaker（2019）對葉際細菌群落研究得出的結果一致。變形菌門中很多細菌對周邊環境的平衡起到重要作用，其中γ-變形菌綱中的假單胞菌屬在4個處理中占據優勢地位，與劉暢等（2020）測定的煙草葉際微生物多樣性的研究結果基本一致。隨著病原菌的侵染假單胞菌屬的相對豐度下降，噴施生防菌后抑制了病害的發生，其相對豐度明顯上升，遠高于其他處理，推測是假單胞菌屬與噴施的生防菌M-23在葉際的互作對植物起到一定的保護作用。沙月霞和沈瑞清（2019）研究發現，3種芽胞桿菌浸種處理顯著改變水稻根和莖部內生細菌群落結構，提高了假單胞菌屬的相對豐度，本研究結果與之一致。假單胞菌屬作為葉際微生物的組成，隨著病害的發生或生防菌的噴施其相對豐度發生變化，可能與植物防御病原物入侵有關。有研究表明，丁香假單胞菌（Pseudomonas syringae）產生的生物活性分子（冠菌素Coronatine和一種蛋白酶體抑制劑Sy- ringolin A）能引起氣孔關閉，從而影響病原體進入質外體（Hunter et al.，2010）。Ko等（2014）分離得到的一株假單胞菌菌株THJ609-3對柑橘砂皮病有較好的防效，可抑制菌絲生長和孢子萌發。此前也有學者發現多種生防菌聯合使用可提高對病害的防治效果，譙天敏等（2015）將銅綠假單胞菌（P. aeruginosa）與長枝木霉菌（Trichoderma longibrachiatum）2種生防菌聯合使用，對雜交竹梢枯病的協同防效得到顯著提高;Wang等（2019）使用枯草芽孢桿菌和哈茨木霉組成的生物產物能有效防治馬鈴薯瘡痂病。因此，芽孢桿菌屬與假單胞菌屬聯合使用以有效防治柑橘砂皮病有待進一步探究。

本研究運用高通量測序技術檢測柑橘葉際細菌群落結構，比較柑橘砂皮病不同發病等級下葉際細菌群落組成的區別和噴施生防菌后葉際菌群的變化，有助于深層次了解柑橘健康葉片與砂皮病葉的微生態差異，為研究柑橘砂皮病的防治提供理論基礎。

4 結論

噴施生防菌M-23能有效防治柑橘砂皮病。在柑橘砂皮病發生較嚴重的柑橘園中噴施生防菌M-23后，M-23菌株能在柑橘葉片上定殖和增殖，改善葉際細菌群落結構，提升有益菌的比例，使病情指數降低，具有一定的應用潛力。

參考文獻：

岑浴. 2016. 多粘類芽孢桿菌在茶葉上的定殖及其對葉際細菌群落的影響[D]. 石家莊：河北科技大學. [Cen Y. 2016. Colonization of Paenibacillus polymyxa in the phyllosphere of Camellia sinensis and its effect on the bacterial community[D]. Shijiazhuang：Hebei University of? Science and Technology.]

方天露，李怡嵐，易圖永，戴良英，何可佳，謝深喜. 2017. 柑橘砂皮病病原菌的分離鑒定和藥劑篩選試驗[J]. 中國果樹，（2）：53-55. [Fang T L，Li Y L，Yi T Y，Dai L Y，He K J，Xie S X. 2017. Isolation，identification and fungicide screening of the pathogen of Diaporthe citri[J]. China Fruits，（2）：53-55.]

高爽，劉笑塵，董錚，劉茂炎，戴良英. 2016. 葉際微生物及其與外界互作的研究進展[J]. 植物科學學報，34（4）：654-661. [Gao S，Liu X C，Dong Z，Liu M Y，Dai L Y. 2016. Advance of phyllosphere microorganisms and their intera-ction with the outside environment[J]. Plant Science Journal，34（4）：654-661.]

李春宏，殷劍美，王立，蔣璐，郭文琦，韓曉勇，張培通. 2019. 連作對芋頭根際土壤理化性狀和微生物特性的影響[J]. 江蘇農業學報，35（4）：825-833. [Li C H，Yin J M，Wang L，Jiang L，Guo W Q，Han X Y，Zhang P T. 2019. Effects of continuous cropping on soil physicochemical and microbial properties of taro[J]. Jiangsu Journal of Agricultural Sciences，35（4）：825-833.]

劉暢，汪漢成，謝紅煉，向立剛，黃宇，陳乾麗，余知和，鄒光進. 2020. 感染赤星病煙草葉際細菌的多樣性分析[J]. 煙草科技，53（2）：8-14. [Liu C，Wang H C，Xie H L，Xiang L G，Huang Y，Chen Q L，Yu Z H，Zou G J. 2020. Biodiversity analysis of phyllosphere bacterial genus from tobacco leaves infected by brown spot disease[J].Tobacco Science and Technology，53（2）：8-14.]

李忠奎，凌愛芬，李紅麗，陳娟，朱顯俊，王勇，陳玉藍，王巖. 2019. 基于多樣性測序對健康與易感病煙田根際土壤微生物群落分析[J]. 河南農業大學學報 ，53（6）：918-925. [Li Z K，Ling A F，Li H L，Chen J，Zhu X J，Wang Y，Chen Y L，Wang Y. 2019. Analysis of rhizosphere soil microbial communities in healthy and susceptible tobacco fields based on diversity sequencing[J]. Journal of Henan Agricultural University，53（6）：918-925.]

羅路云，張卓，金德才，譚新球，張德詠，田云，盧向陽，劉勇. 2017. 南瓜白粉病不同病情等級下葉際細菌群落結構和多樣性[J]. 植物病理學報，47（5）：688-695. [Luo L Y，Zhang Z，Jin D C，Tan X Q，Zhang D Y，Tian Y，Lu X Y，Liu Y. 2017. Phyllosphere bacterial diversities and community structures on pumkin with different severities of powdery mildew[J]. Acta Phytopathologica Sinica，47（5）：688-695.]

馬靖艷，楊付前，劉自君，劉斌，普建才，譚娟. 2019. 懷化地區柑橘砂皮病發生動態及綜合防控措施[J]. 中國植保導刊，39（3）：55-58. [Ma J Y，Yang F Q，Liu Z J，Liu B，Pu J C，Tan J. 2019. Occurrence trends and comprehensive control measures of Diaporthe citri in Huaihua area[J]. China Plant Protection，39（3）：55-58.]

譙天敏，張靜，朱天輝. 2015. 銅綠假單胞菌與長枝木霉對雜交竹梢枯病的協同增效生防研究[J]. 北京林業大學學報，37（2）：113-120. [Qiao T M，Zhang J，Zhu T H. 2015. Synergistic biocontrol mechanism of Pseudomonas aeruginosa and Trichoderma longibrachiatum on Arthri-nium phaeospermum[J]. Journal of Beijing Forestry University，37（2）：113-120.]

沙月霞，沈瑞清. 2019. 芽胞桿菌浸種對水稻內生細菌群落結構的影響[J]. 生態學報，39（22）：8442-8451. [Sha Y X，Shen R Q. 2019. Impact of seed soaked by Bacillus on endophytic bacterial community structure of rice[J]. Acta Ecologica Sinica，39（22）：8442-8451.]

謝蘭芬，何鵬飛，吳毅歆，劉春明，肖春，何月秋. 2018. 解淀粉芽孢桿菌B9601-Y2對玉米葉際微生物群落的影響[J]. 玉米科學，26（2）：149-155. [Xie L F，He P F，Wu Y X，Liu C M，Xiao C，He Y Q. 2018 .Effect of Bacillus amyloliquefaciens 9601-Y2 on maize phyllosphere microbial communities[J]. Journal of Maize Sciences，26（2）：149-155.]

張岳，易圖永，何可佳. 2018. 柑橘砂皮病發生規律研究[J]. 湖南農業科學，（9）：78-80. [Zhang Y，Yi T Y，He K J. 2018. Research on the occurrence regularity of Diaporthe citri[J]. Hunan Agricultural Sciences，（9）：78-80.]

周娜，胡軍華，姚廷山，王雪蓮，王娟，彭鳳格，洪棋斌，江東. 2015. 柑橘種質抗柑橘蒂腐病菌擴展能力的評價[J]. 園藝學報，42（1）：1889-1898. [Zhou N，Hu J H，Yao T S，Wang X L，Wang J，Peng F G，Hong Q B，Jiang D. 2015. Evaluation of anti-expansion capacity of different citrus germplasm against Diaporthe citri[J]. Acta Horticulturae Sinica，42（1）：1889-1898.]

Asnicar F，Weingart G，Tickle T L，Huttenhower C，Segata N. 2015. Compact graphical representation of phylogenetic data and metadata with GraPhlAn[J]. PeerJ，3：e1029.

Bodenhausen N，Horton M W，Bergelson J. 2013. Bacterial communities associated with the leaves and the roots of Arabidopsis thaliana[J]. PLoS One，8（2）：e56329.

Gu S S，Hu Q L，Cheng Y Q，Bai L Y，Liu Z H，Xiao W J，Gong Z H，Wu Y N，Feng K，Deng Y，Tan L. 2019. App-lication of organic fertilizer improves microbial community diversity and alters microbial network structure in tea（Camellia sinensis） plantation soils[J]. Soil and Tillage Research，195：1-10.

Hunter P J，Hand P，Pink D，Whipps J M，Bending G D. 2010. Both leaf properties and microbe-microbe interactions influence within-species variation in bacterial population diversity and structure in the lettuce（Lactuca species） phyllosphere[J]. Apple Environ mental Microbio-logy，76（24）：8117-8125.

Ko Y J，Kim J S，Kim K D，Jeun Y C. 2014. Microscopical observation of inhibition behaviors against Diaporthe citri by pretreated with Pseudomonas putida strain THJ609-3 on the leaves of citrus plants[J]. Journal of Microbiology，52（10）：879-883.

Laforest-Lapointe I，Whitaker B K. 2019. Decrypting the phyllosphere microbiota：Progress and challenges[J]. American Journal of Botany，106（2）： 171-173.

Wang Z S，Li Y，Zhuang L B，Yu Y，Liu J，Zhang L X，Gao Z J，Wu Y F，Gao W，Ding G H，Wang Q. 2019. A rhizosphere-derived consortium of Bacillus subtilis and Trichoderma harzianum suppresses common scab of potato and increases yield[J]. Computational and Structural Biotechnology Journal，17：645-653.

Zhang B G，Bai Z H，Hoefel D，Tang L，Yang Z G，Zhuang G Q，Yang J Z，Zhang H X. 2008. Assessing the impact of the biological control agent Bacillus thuringiensis on the indigenous microbial community within the pepper plant phyllosphere[J]. FEMS Microbiology Letters，284：102-108.

（責任編輯 麻小燕）