张 玮 洪艳云 刘登望 张博文 易图永,* 李 林,* 湖南农业大学植物保护学院 / 植物病虫害生物学与防控湖南省重点实验室, 湖南长沙 408; 湖南农业大学旱地作物研究所,湖南长沙 408

花生(Arachis hypogaeaL.)是一种重要的油料作物, 在我国境内广泛种植[1]。然而我国花生主产区受海洋大气环流影响, 降雨量高, 钙素极易流失, 从而导致土壤酸化缺钙[2]。花生是嗜钙作物, 缺钙会影响果实发育, 形成空壳, 造成减产甚至绝收[3], 严重影响经济效益[4-6]。同时缺钙容易诱发籽仁生理病变使之更易感染真菌病害。在植物与病原菌互作过程中, 钙参与植物的防御反应和组织结构的形成, 增强植物细胞稳定性和抗病能力[7]。钙能直接抑制病原菌菌丝生长、孢子萌发、附着胞形成等,甚至能参与病原真菌的致病性调控[8-11]。

近期研究表明, 微生物菌群可以影响植物和病原物之间的相互作用, 有益微生物能够促进植物的生长, 同时增加植物对非生物胁迫的耐受性[12]。许多根际促生菌能够增强宿主抵抗病原物的能力, 通过诱导系统抗性(induced systemic resistance, ISR)使地上部分能抵抗病原微生物的侵袭和害虫的取食[13-14]。同时植物根系分泌物也能通过改变土壤理化性质吸引附近土壤环境微生物向根部聚集, 进而影响根际微生物的群落结构[15]。Yuan 等[16]证实, 拟南芥在受到叶部病原菌丁香假单胞菌番茄致病变种(Pseudomonas syringaepv.tomato,Pst)侵染后, 可通过增加氨基酸和长链有机酸分泌量等途径改变根系分泌物成分, 从而调控土壤微生物群落来加强自身抵御病害的能力。Molina 等[17]通过将恶臭假单胞菌(Pseudomonas putida)、鞘氨醇单胞菌(Sphingomonassp.)、巴西固氮螺菌(Azospirillum brasilense)以及不动杆菌(Acinetobactersp.)的混合菌群接种至玉米根部,从而提高了玉米的抗旱性。Berendsen 等[14]用活体营养型卵菌拟南芥无色霜霉菌(Hyaloperonospora arabidopsidis)侵染拟南芥叶片后, 通过促进根际3种特异诱导ISR 的细菌增殖, 起到抗霜霉病的作用。这些研究进一步证实了接种微生物混合种比单一菌株在抵御病害与应对胁迫上具有更好的效果。

目前, 植物与微生物互作的研究大多数基于模式生物拟南芥或者水稻、玉米等主要粮食作物, 而对花生根系内生菌群落结构鲜有报道。本研究通过对不同钙处理下花生根系内生菌的群落基因组进行16SrRNA 基因V3–V4 区测序, 试图探析对钙敏感花生植株在酸性低钙条件下, 施钙处理对其根系微生物群落结构的影响, 为合理施用矿质营养元素、改良酸性土壤品质、降低花生病害发生提供理论依据。

1 材料与方法

1.1 供试材料

试验地点为湖南农业大学耘园科研基地。花生采用国家主推的大籽品种湘花2008, 由湖南农业大学旱地作物研究所选育, 对钙敏感; 酸性低钙土壤取自湖南省浏阳市普迹镇书院村月光坪的第四纪红壤表层土, pH 4.5, 含有机质27.3 g kg-1、全氮(N)0.82 g kg-1、全磷(P) 0.33 g kg-1、全钾(K) 17 g kg-1、碱解氮(N) 64 mg kg-1、速效钾(K) 82 mg kg-1、阳离子交换量13 cmol (+) kg-1、交换性钙(Ca) 0.74 cmol(1/2 Ca2+) kg-1、交换性镁(Mg) 0.13 cmol (1/2 Mg2+)kg-1、有效锌(Zn) 0.63 mg kg-1、有效硫(S) 143.7 mg kg-1、有效硼(B) 0.12 mg kg-1、有效钼(Mo) 0.08 mg kg-1, 有效磷(P)未检出, 为典型酸性缺钙红壤(依发明专利“花生种质钙敏感性的土壤鉴定方法”进行土壤鉴定, 具体鉴定方法详见专利附录)。

选取内径36 cm、高75 cm、壁厚度0.5 cm 的PVC 管组装成塑料大桶, 一端置于地膜之上以防止供试土壤与试验地原有土壤接触。每桶内装供试土壤90 kg。施钙组每桶施加生石灰(CaO) 12.5 g (125 kg hm-2), 与桶中表层供试土壤充分混匀, 对照组不施加生石灰(CaO)。为保证除钙以外的其他基本养分平衡一致的供给, 本试验除进行施钙及对照处理外,各样品组其他养分均依临界值进行养分平衡: 于土柱表层0～20 cm 处, 参照大田标准(45%氮磷钾等比例复合肥750 kg hm-2)每桶施尿素4.02 g、磷酸二氢钾8.23 g 及六水氯化镁8.50 g, 并混匀。施钙组与对照组各相同处理10 桶作为重复。依十字将大桶划分为4 个区域, 每桶定苗4 株(作为不同生育期采集检测样本), 始花期每桶喷施等量硼肥。花生生长发育阶段进行正常的水分和病虫草害的管理。

1.2 试验方法

1.2.1 样品采集及处理 分别于苗期(2019-06-04)、花针期(2019-06-24)与饱果期(2019-08-14)采集施钙组与对照组植株根系样品。每组设置3 个重复。将各组样品植株全株拔出, 抖动根系以去除疏松的大块土壤。用2 把无菌剪刀分别在施钙组与对照组样品根茎交界处正下方切下约5 cm 的根。抖落大块土壤, 并用灭菌毛刷清除表面小颗粒状土壤, 再用无菌水多次反复冲洗样本表面, 并摘去小侧根。将处理后的根放入5%次氯酸钠溶液浸泡5 min, 再转移至装有适量无菌玻璃珠及无菌水的50 mL 灭菌离心管中, 剧烈摇晃以便物理去除根表土壤微生物, 用无菌水重复振荡冲洗3 次。称取2 g 根系组织, 用灭菌剪刀尽可能剪碎, 放入灭菌三角瓶, 加入20 mL 无菌的1×PBS 缓冲液, 摇床振荡4 h, 40 KHz 超声波超声振荡15 min, 5 μm 滤膜真空抽滤。将滤液15,000×g离心15 min, 弃上清收集菌体沉淀。使用EasyPure Genomic DNA Kit 试剂盒提取菌体沉淀中的根系内生细菌基因组DNA, 使用NanoDrop2000 检测所提样品的核酸质量, 合格后送天津诺禾致源生物信息科技有限公司进行16S 扩增子测序, 样品编号如表1 所示。

表1 16S 扩增子测序样品分组编号Table 1 Group number of 16S amplicon sequencing sample

1.2.2 16S 扩增子测序 取适量样本DNA 于离心管中, 使用无菌水稀释样本至1 ng μL-1, 以稀释后的基因组DNA 为模板, 根据测序区域16S V3-V4 选择 特 异 引 物 341F (Forward primer 5′-CCTAYG GGRBGCASCAG-3′)与806R (Reverse primer 5′-GG ACTACNNGGGTATCTAAT-3′)进行PCR。PCR 产物使用 2%浓度的琼脂糖凝胶进行电泳检测后使用GeneJET 试剂盒纯化回收目标条带。使用Ion Plus Fragment Library Kit 试剂盒进行文库构建。构建好的文库经质检合格后使用Thermofisher 的Ion S5 XL进行上机测序, 对测序后下机的原始数据进行拼接和过滤, 得到Clean Reads。

1.2.3 OTUs 聚类分析 默认以 97%的一致性(Identity)对 Clean Reads 进行 OTUs (operational Taxonomic units)聚类, 然后通过与数据库Silva132的比对进行物种注释。分析微生物群落结构差异并绘制差异OTUs 韦恩图、UPGMA 聚类树及属水平下物种相对丰度柱形图。选取属分类水平下Top40优势物种, 使用SPSS 软件对施钙处理与对照下不同生育期各物种相对丰度差异进行独立样品T 检验统计分析, 显著差异的分析结果使用OriginPro 2019软件绘制柱状图。

1.2.4 LEfSe 分析 LEfSe (LDA effect size)是一种用于发现和解释高维度生物标识(基因、通路和分类单元)的分析工具, 它强调统计意义和生物相关性,能够在组与组之间寻找具有统计学差异的Biomarker, 让研究人员能够识别不同丰度的特征以及相关联的类别, 以便找出哪些核心菌群引起组间群落的显著差异。本研究以3 作为LDA Score 的筛选值对各样品组测序数据进行LEfSe 分析并绘制进化分支图。在进化分支图中, 由内至外辐射的圆圈代表了由门至属(或种)的分类级别。在不同分类级别上的每一个小圆圈代表该水平下的一个分类, 小圆圈直径大小与相对丰度大小呈正比。着色原则为无显著差异的物种统一着色为黄色, 差异物种Biomarker 跟随组进行着色。

1.2.5 共发生网络构建 分别对施钙组与对照组样品中相对丰度排名Top100 的微生物属进行Pearson相关系数计算, 以优势微生物属作为核心节点, 选取Pearson 相关系数大于0.6 或小于-0.6 且显著差异(P<0.05)的连接作为有效连接, 使用Cytoscape 软件进行共发生网络图绘制, 以便更加清晰的看出不同样品组中某个微生物属和其他微生物的关联状况。

2 结果与分析

2.1 根系微生物基因组16S V3-V4 的PCR 扩增结果

以稀释后基因组DNA 为模板, 选择16S V3-V4的特异引物341F (5′-CCTAYGGGRBGCASCAG-3′)与806R (5′-GGACTACNNGGGTATCTAAT-3′)进行PCR 扩增。部分样品基因组特异引物扩增后电泳检测结果如图1 所示。

2.2 16S 扩增子测序数据分析

2.2.1 OTU 聚类及物种注释概况 经过物种注释及对不同分类层级统计分析发现, 共有2378 个OTUs, 其中能够注释到数据库的OTUs 数目为2372(99.75%), 在门水平上, 发现占据主导地位的物种包括厚壁菌门 (Firmicutes) 、变形菌门(Proteobacteria)、蓝藻门(Cyanobacteria)、拟杆菌门(Bacteroidetes)及放线菌门(Actinobacteria); 在纲水平上主要有芽孢杆菌纲(Bacilli)、γ-变形菌纲(Gammaproteobacteria)、梭菌纲(Clostridia)、α-变形菌纲(Alphaproteobacteria)及拟杆菌纲(Bacteroidia);在目水平上的优势物种为芽孢杆菌目(Bacillales)、肠杆菌目(Enterobacteriales)、根瘤菌目(Rhizobiales)、梭 菌 目 (Clostridiales) 以及黄单胞菌目(Xanthomonadales); 在科水平上的主要物种为芽孢杆菌科(Bacillaceae)、肠杆菌科(Enterobacteriaceae)、伯克氏菌科(Burkholderiaceae) 、 根瘤菌科(Rhizobiaceae)及黄单胞菌科(Xanthomonadaceae)。

2.2.2 花生根系内生菌群落结构 为分析不同生育期花生根系内生细菌的群落结构, 研究不同样本组之间的相似性, 本研究通过对样本以Binary Jaccard 距离矩阵进行UPGMA 等级聚类, 构建样本的聚类树。同时根据聚类得到的OTUs 分析不同样本组之间共有及特有的OTUs 并绘制成韦恩图, 如图2 所示。根据差异OTUs 韦恩图分析可以发现, 除43 个共有的核心OTUs 以外, 随着花生植株的生长发育, 各生育期特有的OTUs 逐渐减少: 苗期施钙与不施钙样品组特有的OTUs 数目为177、117; 花针期施钙组特有24 个差异OTUs, 而不施钙组特有47 个差异OTUs; 饱果期施钙组与对照组的差异OTUs 数分别为5、8。

选取不同生育期不同处理下各样品组在属(Genus)分类水平上相对丰度排名前20 的物种绘制物种相对丰度柱形图, 以便直观查看各样品组中相对丰度较高的优势物种及其所占比例(图3)。根据测序结果, 统计了各生育期不同处理样品组的优势物种(相对丰度排名前5 的优势属), 如表2 所示。苗期不施钙样品组优势菌属有不动杆菌属(Acinetobacter:16.9%)、克雷伯氏菌属(Klebsiella: 16.8%)、肠杆菌属(Enterobacter: 7.5%) 、厌氧芽孢杆菌属(Anoxybacillus: 5.0%)、草螺菌属(Herbaspirillum:4.9%)等; 苗期施钙样品组优势菌属包括克雷伯氏菌属(Klebsiella: 30.1%)、芽孢杆菌属(Bacillus: 26.3%)、肠杆菌属(Enterobacter: 6.1%) 、水栖菌属(Enhydrobacter: 5.1%) 、慢生根瘤菌属(Bradyrhizobium: 1.8%)等; 花针期不施钙样品组优势菌属有克雷伯氏菌属(Klebsiella: 16.9%)、肠杆菌属(Enterobacter: 8.9%)、不动杆菌属(Acinetobacter:6.7%)、泛菌属(Pantoea: 5.0%)、戴氏菌属(Dyella:3.7%)等; 花针期施钙样品组优势菌属有克雷伯氏菌属 (Klebsiella: 37.2%) 、寡养单胞菌属(Stenotrophomonas: 13.1%) 、假单胞菌属(Pseudomonas: 8.7%)、从毛单胞菌属(Comamonas:8.3%)、肠杆菌属(Enterobacter: 7.4%)等; 饱果期不施钙样品组优势菌属为克雷伯氏菌属(Klebsiella:33.4%)、肠杆菌属(Enterobacter: 11.3%)、类芽孢杆菌属(Paenibacillus: 6.3%)、假单胞菌属(Pseudomonas:4.9%)、芽孢杆菌属(Bacillus: 1.7%)等; 饱果期施钙样品组优势菌属为克雷伯氏菌属(Klebsiella: 30.9%)假单胞菌属(Pseudomonas: 8.9%)、类芽孢杆菌属(Paenibacillus: 8.7%) 、 赖氨酸芽孢杆菌属(Lysinibacillus: 8.3%)、芽孢杆菌属(Bacillus: 3.2%)等; 部分主要优势物种未鉴定到属分类水平, 如伯克氏菌科 (Burkholderiaceae) 、 根瘤菌科(Rhizobiaceae)、 梭菌目(Clostridiales)、蓝藻门(Cyanobacteria)等故而未列入上述不同样品组中优势物种排名。综合各样品组数据发现, 克雷伯氏菌属(Klebsiella)、肠杆菌属(Enterobacter)在所有样品组中均为优势物种且相对丰度较高。

表2 属分类水平下不同生育期不同处理相对丰度百分比Top 10 的根系内生菌Table 2 Percentage of the Top 10 relative abundance under different growth stage and different treatments at the genus level

2.2.3 优势菌属差异 为了寻找不同处理下不同生育期样品组中具有显著差异的微生物属, 选取属分类水平下Top40 优势物种, 对施钙处理与对照下的不同生育期各物种相对丰度差异进行独立样品T检验(T-test)统计分析(图4)发现, 施钙组与对照组样品在苗期与饱果期时, 虽然部分微生物属在相对丰度上存在差异, 但并不显著(P＞0.05), 而花针期施钙处理与对照之间, 假单胞菌属(Pseudomonas:P=0.025)、赖氨酸芽孢杆菌属(Lysinibacillus:P=0.037)及劳尔氏菌属(Ralstonia:P=0.038)的相对丰度存在显著差异(P≤0.05)。其中假单胞菌属(Pseudomonas)与赖氨酸芽孢杆菌属(Lysinibacillus)施钙组中的相对丰度较不施钙组明显增加而劳尔氏菌属(Ralstonia)的相对丰度则在施钙组中显著降低。前人研究发现假单胞菌属(Pseudomonas)与赖氨酸芽孢杆菌属(Lysinibacillus)的部分菌株具有抑菌活性, 降解有害物质甚至抵御病害侵袭的作用, 而劳尔氏菌属(Ralstonia)与花生青枯病病原菌在分类水平上为同一属。施钙对这些优势微生物属相对丰度的显著影响极有可能与这些菌属的生物学功能息息相关。施钙可能在一定程度上能抑制某些病原菌生长, 从而有利于花生抵御相关病害。

LEfSe 统计分析结果如图5 (LEfSe 进化分支图)所示, 红色节点表示在苗期不施钙样品组中起到重要作用的微生物类群, 绿色节点表示在苗期施钙样品组中起到重要作用的微生物类群, 蓝色节点表示在花针期不施钙样品组中起到重要作用的微生物类群, 紫色节点表示在花针期施钙样品组中起到重要作用的微生物类群, 浅蓝色节点表示在饱果期不施钙样品组中起到重要作用的微生物类群, 橙色节点表示在饱果期施钙样品组中起到重要作用的微生物类群。进化分支图结果显示, 紫色扇形区域与橙色扇形区域较大且相关节点较多, 表明花针期与饱果期施钙样品组中能起到关键生物学作用的微生物类群较多。属分类水平上: 施钙样品组中苗期的慢生根瘤菌属(Bradyrhizobium), 花针期的希瓦氏菌属(Shewanella)、假单胞菌属(Pseudomonas)及寡养单胞菌属(Stenotrophomonas), 饱果期的类芽孢杆菌属(Paenibacillus)、赖氨酸芽孢杆菌属(Lysinibacillus)与短波单胞菌属(Brevundimonas)等优势物种分别在不同的生育期起关键性生物学作用。

2.2.4 Network 共发生网络分析 网络分析被广

泛应用于复杂微生物群落的探究中。为进一步探究不同处理下花生根系内生细菌间的相互作用, 本研究计算差异菌属两两之间pearson 相关系数并构建了共发生网络(图6)。在施钙组微生物菌属共发生网络中, 核心菌属类诺卡式菌属(Nocardioides)与18 个菌属显著正相关, 溶杆菌属(Lysobacter)和中生根瘤菌属(Mesorhizobium)均与17 个菌属显著正相关, 细链孢菌属(Catenulispora)与16 个菌属显著正相关。而相应的在不施钙处理下样品组根系内生菌关联互作网络中, 发现核心菌属潘多拉菌属(Pandoraea)与其他菌属之间的总关联系数最高, 而戴氏菌属(Dyella)与 16 个菌属显著正相关。丛毛单胞菌属(Comamonas)、罗河杆菌属(Rhodanobacter)、芽孢杆菌属(Bacillus)、劳尔氏菌属(Ralstonia)等优势菌属也与其他菌属之间存在不同程度的关联互作。推测在施钙处理的花生根系内生细菌群落中, 类诺卡式菌属(Nocardioides)、溶杆菌属(Lysobacter)、中生根瘤菌属(Mesorhizobium)以及细链孢菌属(Catenulispora)作为核心菌属具有重要调控作用, 且对网络内其他微生物种群影响较大。而在不施钙处理的花生根系内生细菌群落中, 潘多拉菌属(Pandoraea)、戴氏菌属(Dyella)、丛毛单胞菌属(Comamonas)、罗河杆菌属(Rhodanobacter)、芽孢杆菌属(Bacillus)、劳尔氏菌属(Ralstonia)等可能为该微生物群落的核心调控菌属, 对其他微生物种群具有较大的影响。同时可以明显看出, 以类诺卡式菌属(Nocardioides)为核心的施钙组菌群之间互作联系更为紧密, 有效关联更多。表明施钙能在一定程度上增强根系内生菌群之间的互作联系, 使相互关联紧密。

3 讨论

植物根系是植物-微生物互作的重要生境场所。在此微生态区域, 微生物更易受到植物种类、土壤类型、植物生长阶段等多因素的影响。本研究采用高通量测序技术, 摒弃了繁琐复杂的分离纯化菌株,以新的视角全面探究, 比较分析了施钙处理与对照之间花生根系内生菌群落结构的显著差异: 差异OTUs 分析发现, 苗期与饱果期施钙样品组中所含根系内生细菌差异OTUs 数显著高于不施钙样品,即这2 个不同生育期施钙花生样品中根系内生细菌多样性较对照组更为丰富。优势菌属差异分析发现,假单胞菌属(Pseudomonas)与赖氨酸芽孢杆菌属(Lysinibacillus)的相对丰度在施钙组中显著增加, 而劳尔氏菌属(Ralstonia)相对丰度显著降低。克雷伯氏菌属(Klebsiella)、肠杆菌属(Enterobacter)、类芽孢杆菌属 (Paenibacillus) 、寡养单胞菌属(Stenotrophomonas)等优势物种作为施钙与对照组间显著差异菌群分别在不同的生育期施钙组样品中起关键性生物学作用。同时, 共发生网络的构建展示出施钙能在一定程度上使根系内生细菌种群之间联系更为紧密, 增强优势菌群之间的关联互作; 同时分析结果还表明, 类诺卡式菌属(Nocardioides)、溶杆菌属(Lysobacter)、细链孢菌属(Catenulispora)等是施钙组根系内生菌群落中的核心菌属。

Ibanez 等[18]的研究发现, 克雷伯氏菌属(Klebsiella)能在花生根系结瘤固氮和促进植株生长,Sharm 等[19]发现, 克雷伯氏菌 MBE02 (KlebsiellaMBE02)能诱导花生产生ISR 同时抵御真菌病害。Gong 等[20]的研究论述了阿氏肠杆菌Vt-7 (Enterobacter asburiaeVt-7)对黄曲霉的抑菌活性及对黄曲霉毒素的减毒作用, 与此同时Yang 等[21]分离的荧光假单胞菌3JW1 (Pseudomonas fluorescens3JW1)、Wang 等[22]发现的青岛假单胞菌新种(Pseudomonas qingdaonensissp. nov.)、Wang 等[23]首次报道的解木糖赖氨酸芽孢杆菌(Lysinibacillus xylanilyticus)等菌株同样具有抑制黄曲霉菌活性, 降解黄曲霉毒素或者产生黄曲霉毒素抑制剂的作用。Haggag 等[24]研究证实, 多粘类芽孢杆菌(Paenibacillus polymyxa)能有效生物防治花生冠腐病。魏兰芳等[25]发现, 抗生素溶杆菌13-1 (Lysobacter antibioticus13-1)菌株产生的抗菌物质能对水稻白叶枯病具有防治效果。类诺卡式菌属(Nocardioides)大多能降解有毒物质, 对石油污染土壤进行生物修复, 雒晓芳等[26]却分离鉴定并探究了白色类诺卡氏菌(Nocardioides albus)对部分病原菌的抑菌活性。云天艳等[27]还从木薯根际分离到能拮抗尖孢镰刀菌的专化型放线菌MS13, 鉴定后属于细链孢菌属(Catenulispora)。

结合这些以往的研究成果可以合理推测, 克雷伯氏菌属(Klebsiella)作为所有样品组中的优势物种在某种程度上具有结瘤固氮促进植株生长以及提高花生产量的作用。施钙能显著增加假单胞菌属(Pseudomonas)与赖氨酸芽孢杆菌属(Lysinibacillus)在根系内生菌群落中的相对丰度, 联合肠杆菌属(Enterobacter)、类芽孢杆菌属(Paenibacillus)、寡养单胞菌属(Stenotrophomonas)等优势物种能在一定条件下抑制黄曲霉活性同时缓解黄曲霉毒素的毒性甚至降低花生冠腐病的发生概率。与对照组相比, 相对丰度显著降低的劳尔氏菌属(Ralstonia)则间接说明施钙极有可能同时抑制花生青枯病的发生。而类诺卡式菌属(Nocardioides)、溶杆菌属(Lysobacter)、细链孢菌属(Catenulispora)等作为施钙组中的核心菌属能对部分病原菌产生抑菌活性。

4 结论

施钙能丰富根系内生细菌群落的多样性, 使群落间更为紧密的关联互作; 也能使根系微生物组具有更好的促生作用, 使植株能更好地抵御病害侵袭;施钙能增强特定物种的相对丰度从而促进花生植株的生长发育, 提高产量与品质同时可能抑制常见病害如花生冠腐病、花生青枯病的发生, 抑制黄曲霉菌活性甚至降解黄曲霉毒素; 这些研究结果在一定程度上为南方酸性红壤条件下种植花生提供了施肥管理及病害防控等方面的合理模式, 同时为如何改良酸性土壤品质及提高花生产量与品质奠定了理论基础。