基于高通量测序的黑柄乳菇菌塘细菌多样性分析

秦婉琳 马永远 李双燕 伍燕

摘 要 乳菇系贵州食用菌“黔菌”品牌重点打造的项目之一。黑柄乳菇（Lactarius gerardii）属于外生菌根菌，是一种美味的野生菌，市场价格不菲。为了解野生黑柄乳菇的菌根生境，给下一步研发黑柄乳菇的人工栽培技术奠定一些基础，以黑柄乳菇菌塘土壤为研究对象，采用高通量测序技术，结合生物信息学对来自于黔西南州不同地区的5个菌塘土壤的细菌组成及丰度进行研究。结果表明：黑柄乳菇菌塘菌群分属于9个门27个属，菌塘中主要门有变形菌门、酸杆菌门、放线菌门和浮霉菌门;主要属有伯克氏菌属、缓生根瘤菌属、Gp1、Gp2、Gp3、假单胞菌属和WPS-1 genera incertae sedis。五个采集点菌塘中的菌群丰富，多样性指数水准高，但因各个菌塘生境不同使得每个菌塘的优势菌属丰度有差异。

关键词 黑柄乳菇;菌塘;高通量测序;细菌;菌落多样性

中图分类号：S646.1;S182 文献标志码：A DOI：10.19415/j.cnki.1673-890x.2022.11.001

黑柄乳菇（Lactarius gerardii）又名茅草菌，隶属伞菌目（Agaricales）红菇科（Russulaceae）乳菇属（Lactarius），菌盖和菌柄呈黑褐色，菌褶白色延伸至菌柄，菌肉白色、质脆，受伤后会溢出白色乳浆，因此又叫黑盖奶浆菌;夏季多生长在针阔混交林地上，是一种美味的野生菌[1]。图1（参见封三）为野生黑柄乳菇子实体，采集地点为兴义市万屯镇云屯生态公园松林下。乳菇属在贵州比较常见，例如松乳菇（紫花菌）、多汁乳菇、红汁乳菇、靓丽乳菇等[2-3]。近年研究表明，通过人工合成松乳菇菌体并接种到天然马尾松林下，能得到松乳菇子实体[4-9]。松乳菇营养丰富，含多种氨基酸和微量元素，具有增强机体免疫力、抗氧化和抗肿瘤等生物活性[10-15]，松乳菇的栽培成功对开发利用黔西南地区的黑柄乳菇具有参考价值。

菌塘（Shiro）是由子实体的根部菌丝、宿主根系及周围的腐殖质土壤扭结形成的颗粒状结构，是真菌菌丝体萌发和形成子实体的重要条件[16-17]。深入研究菌塘，对于开发人工驯化技术，了解菇类栽培土壤细菌群的组成及丰度非常重要。本文以黔西南州的野生黑柄乳菇菌塘土壤为研究对象，高通量测序菌塘菌群后采用统计聚类的方法确定不同样本优势菌群结构的分型，旨在探究黑柄乳菇生境土壤微生物群落特征，找出优势菌群、分析共性并得出规律性认识，为下一步研发黑柄乳菇的人工栽培技术奠定一些基础。

1  材料与方法

1.1  材料

1.1.1  样品采集

选取黑柄乳菇出产丰富的五个采集点：兴义市万屯镇云屯生态公园（WT）、安龙县新桥镇（XQ）、安龙县海子镇（HZ）、安龙县平乐镇（PL）、安龙县龙广镇（LG）。每个地点选取黑柄乳菇成簇生长的菌塘，在距菌根2 cm半径周围处取土样，混合后放入4 ℃冰箱保存。

1.1.2  药品及试剂

Qubit3.0 DNA检测试剂盒，美国Life公司;E.Z.N.ATM Mag-Bind Soil DNA Kit试剂盒，美国OMEGA公司;2×Taq Master Mix，南京诺唯赞生物技术公司;MagicPure Size Selection DNA Beads，北京全式金生物技术公司。

1.2  仪器和设备

Pico-21台式离心机，美国赛默飞公司;Bio-Rad PCR扩增仪、凝胶成像仪，美国伯乐公司;Q32866 Qubit® 3.0荧光计，美国英杰生物技术公司;DYCZ-21电泳槽、DYY-6C型电泳仪电源，北京市六一仪器厂。

1.3  方法

1.3.1  样品预处理

称取200 mg土样，放入灭菌的2 mL离心管中，加入1 mL PBS缓冲液混匀，10 000 rpm室温离心3 min，弃上清，样品干燥待用。

1.3.2  总DNA提取和检测

按照OMEGA试剂盒E.Z.N.ATM Mag-Bind Soil DNA Kit使用说明操作提取总DNA，1%琼脂糖凝胶检测DNA完整性。

1.3.3  PCR扩增及测序

宏基因组测序类群为细菌;扩增区域为细菌V3～V4变区，引物1为341F：CCTACGGGNGGCWGCAG，引物2为805R：GACTACHVGGGTATCTAATCC，测序平台为Miseq，由生工生物工程（上海）公司完成。

1.3.4  数据分析

高通量测序序列（reads）通过拼接、过滤、去除嵌合体及非特异性扩增序列，进行OTU（Operational Taxonomic Unit）聚类分析和物种分类学分析。基于OTU聚类分析结果，对OTU进行多样性指数分析、多样性的群落组成分析，对系统发育信息进行Beta多样性分析、分组检验分析、差异显著性检验、环境因子分析等深入的统计学和可视化拟合。

2  结果与分析

2.1  菌塘Alpha多样性指数分析

表1为各个菌塘土壤高通量测序数据Alpha多样性指数分析结果，可以反映所测微生物群落的丰度和多样性。由表可知，五个菌塘土样测序共得到140 903条有效序列，过滤后得到4 385个有效OTU，样品的Coverage指数≥0.99，覆盖率高，能真实反映样品的测序情况。

样品的ACE指数和Chao指数反映群落分布丰度（Communit richness），数值越大说明样本的菌群分布丰度越好;Shannoneven指数越高或Simpson指数越低，说明采样点的菌群多樣性较丰富。从表1可以看出，龙广（LG）菌落较丰富、多样性较好。

图2为五个菌塘菌群样品的稀释曲线，Shannon指数先呈直线上升，后趋向平稳，表明测序数据能较好地覆盖菌塘中菌群的信息序列，曲线在平滑稳定后产出的OTU数趋于稳定，其中来自海子（HZ）的样品Shannon指数稍低，说明该采集地的细菌多样性低于其余四个菌塘。

2.2  各样品中细菌类群的相关性分析

图3（参见封三）为五个菌塘菌群基于OTU序列所作的聚类树，可以看出HZ和PL聚成一类，并和LG、XQ位于同一簇上，WT位于较远分支，直观反映了各个菌塘菌群的相似性和差异关系。

图4（参见封三）的韦恩图为各个菌塘菌群间OTU关联分析，反映了各自生境中土壤菌群的复杂性和共性。五个菌塘共有的OTU数是471个，各个菌塘菌群互有重叠，反映的是菌塘间的共性;不重叠部分除了万屯（WT）较多，有104个OTU外，其余菌塘均较少。除了万屯（WT）位于生态公园内，人类活动影响较小之外，其余四个采集地都在安龙县农业生产区。数据分析结果直观反映出各菌塘土壤基质微生物种类相似性和差异性。

2.3  菌塘中细菌组成分析

五个菌塘细菌菌群归属于9个门，即变形菌门（Proteobacteria，68.88%～43.97%）、酸杆菌门（Acidobacteria，30.24%～10.33%）、放线菌门（Actinobacteria，15.70%～4.78%）、浮霉菌门（Planctomycetes）、拟杆菌门（Bacteroidetes）、疣微菌门（Verrucomicrobia）、candidate_division_WPS-1、绿弯菌门（Chloroflexi）和candidate_division_WPS-2，从图5（参见封三）中可知，以变形菌门、酸杆菌门和放线菌门为主。

由图6（参见封三）可见，各菌塘在属级水平上比较丰度，有27个属，平乐（PL）菌塘主要的属有缓生根瘤菌属（Bradyrhizobium，7.52%）、伯克氏菌属（Burkholderia，5.89%）、Gp1（7.02%）;万屯（WT）菌塘主要的属有缓生根瘤菌属（Bradyrhizobium，6.35%）、根瘤菌属（Rhizobium，5.6%）、假单胞菌属（Pseudomonas，5.55%）;新桥（XQ）菌塘主要的属有伯克氏菌属（Burkholderia，9.53%）、假单胞菌属（Pseudomonas，5.67%）、Gp1（4.39%）;海子（HZ）菌塘主要的属有Gp1（10.24%）、Gp2（9.5%）、慢生根瘤菌属（Bradyrhizobium，8.86%）、伯克氏菌属（Burkholderia，6.54%）;龙广（LG）菌塘主要的属有缓生根瘤菌属（Bradyrhizobium，6.51%）、Gp2（4.72%）、Gp1（4.16%）。

结合各个菌塘门、属水平统计及表2所列前几种优势属可知，各个菌塘中的优势菌门和属均是变形菌门及其下属，这些菌群能起到固氮及吸收转化磷、钾等无机盐的作用，并与黑柄乳菇菌丝、土壤腐殖质等扭结形成菌根菌共生体，外部表现为根瘤形态。同时注意到各个菌塘的酸杆菌门及酸杆菌属Gp1、Gp2和Gp3含量较高，通常认为这些菌增多可对被污染土壤进行修复。

3  结论

采用高通量测序技术结合生物信息学对黔西南州的野生黑柄乳菇五个菌塘进行分析，对菌群的组成和丰度形成了一个规律性的认识。各个菌塘中以变形菌门、酸杆菌门、放线菌门为主，菌群的组成类型变化不大但丰度上有区别，例如万屯（WT）菌塘变形菌门丰度达到68.88%，龙广（LG）菌塘变形菌门丰度为43.97%;海子（HZ）菌塘酸杆菌门丰度为30.24%，万屯（WT）菌塘酸桿菌门丰度为10.33%。采样点万屯镇云屯生态公园环境保护相对较好，土壤有机质丰富，黑柄乳菇生活的菌塘中以固氮和转化磷、钾元素为主的变形菌门居多[18-19];海子镇是安龙县粮食主产区，自然环境受人类生产活动影响较大，与万屯相比酸杆菌门丰度较高（通常认为在人类生产活动较多的环境，例如过度施肥、喷施农药等会引起土壤中酸杆菌属含量升高）。酸杆菌门的Gp1、Gp2、Gp3能降低土壤中的重金属和农药污染[20]，有助于恢复土壤生态。

乳菇作为贵州食用菌“黔菌”品牌重点打造的项目之一，其中的松乳菇已有少量人工栽培，但远不能满足市场需要。兴义野生黑柄乳菇在当地出产丰富，味道鲜美，深受消费者喜爱，市场价格也在逐年攀升。对黑柄乳菇的菌塘菌群组成和丰度开展研究，进行深入了解，对下一步做“精”、做“特”食用菌产业具有理论指导意义。通过本次分析，发现生境中的酸杆菌门菌群较丰富，说明人类在野生菌产出区域活动过于频繁，对土壤的理化性质已造成影响，提醒我们注意保护当地生态环境，保护野生黑柄乳菇及其他野生食用菌的可持续生长和出产品质。

参考文献：

[1]  戴玉成，周丽伟，杨祝良，等.中国食用菌名录[J].菌物学报，2010，29（1）：1-21.

[2]  易筑刚，陈春旭，陈华，等.贵州两种乳菇的鉴定及同源性分析[J].福建农业学报，2021，36（7）：766-770.

[3]  陈新，康恒，边银丙.乳菇类外生菌根食用菌研究进展[J].菌物学报，2018，37（12）：1562-1571.

[4]  郑华英，徐丽丽，徐明.松乳菇松林下栽培及产菇情况初步研究[J].江苏林业科技，2020，47（5）：28-31.

[5]  王迪，何晓兰，熊卫萍，等.人工合成松乳菇菌根及杂菌侵染宿主的形态学观察和分子鉴定[J].食用菌学报，2021，28（3）：129-134.

[6]  易筑刚，陈春旭，陈华，等.松乳菇目的菌宿主辨识方法及非目的菌清除实验效果分析[J].贵州林业科技，2021，49（4）：46-50.

[7]  谭云，申爱荣，谭著明，等.红汁乳菇培养基和菌株优选[J].湖南林业科技，2021，48（6）：1-8.

[8]  徐彦军.3种松乳菇菌种在不同培养基中生长表现[J].种子，2020，39（11）：107-109.

[9]  郑旋，康超，杨玲，王晶，等.不同分离培养条件对松乳菇菌丝的影响研究[J].种子科技，2020，38（15）：16-17.

[10] 王秋艳，丁慧敏，朱亚男，等.多汁乳菇多糖对小鼠急性酒精性肝损伤的保护作用[J].食品工业科技，2021，42（24）：313-319.

[11] 宋盛英，王勇，苏上，等.橙黄色、紫褐色松乳菇的营养成分比较[J].食用菌，2020，42（1）：60-62.

[12] 白瑞，鮑中英，段淑红，等.松乳菇多糖治疗原发性肝癌的细胞实验及机制[J].临床与病理杂志，2019，39（10）：2110-2117.

[13] 陈璇，高慧，冯丽娜，等.燕山山脉一种香料型乳菇的香味成分分析[J].菌物学报，2019，38（10）：

1670-1680.

[14] 龚燕思，林菲菲，黄仁鸿，等.多汁乳菇对糖尿病小鼠血糖的影响[J].中国临床药理学杂志，2019，35（18）：2068-2070.

[15] 江波，胡芬，何标.松乳菇多糖通过PI3K/Akt信号通路对人骨肉瘤细胞MG-63裸鼠移植瘤生长的抑制作用[J].中国药师，2019，22（9）：1591-1595.

[16] 徐漫，傅婉秋，戴传超，等.外生菌根真菌促生微生物生态功能研究进展[J].生态学杂志，2018，37（4）：1246-1256.

[17] 于文海，左园园，刘平稳.外生菌根食用菌研究概况[J].现代食品，2018（15）：23-26.

[18] 刘平静，肖杰，孙本华，等.长期不同施肥措施下土细菌群落结构变化及其主要影响因素[J].植物营养与肥料学报，2020，26（2）：307-315.

[19] 吴庆梅，黄宝灵，吕成群，等.不同松林根际促生菌的数量研究[J].广西农业科学，2009，40（5）：521-526.

[20] 王继玥，刘燕，刘勇，等.基于高通量测序检测Pb污染对三叶草根际土壤细菌多样性的影响[J].环境科学研究，2018，31（1）：102-110.

收稿日期：2022-01-25

基金项目：2021年贵州省大学生创新创业计划项目（S202110666097）。

作者简介：秦婉琳（2001—），女，贵州毕节人，在读本科生。

*为通信作者：伍燕（1974—），女，硕士，教授，主要从事天然产物开发应用研究。E-mail：767692523@qq.com。