基于高通量测序的黑柄乳菇菌塘细菌多样性分析
黑柄乳菇（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 ）。每个地点选取黑柄乳菇成簇生长的菌塘，在距菌根2cm 半径周围处取土样，混合后放入4℃冰箱保存。1.1.2
药品及试剂
Qubit3.0DNA 检测试剂盒，美国Life 公司；E.Z.N.A TM Mag-Bind Soil DNA Kit 试剂盒，美国OMEGA 公
司；2×Taq Master Mix ，南京诺唯赞生物技术公司；
收稿日期：2022-01-25
基金项目：2021年贵州省大学生创新创业计划项目（S202110666097）。作者简介：秦婉琳（2001—），女，贵州毕节人，在读本科生。*为通信作者：伍燕（1974—），女，硕士，教授，主要从事天然产物开发应用研究。E-mail ：****************。
秦婉琳，马永远，李双燕，等.基于高通量测序的黑柄乳菇菌塘细菌多样性分析[J ].南方农业，2022，16（11）：1-4.
基于高通量测序的黑柄乳菇菌塘细菌多样性分析
秦婉琳，马永远，李双燕，伍燕*
（兴义民族师范学院，贵州兴义562400）
摘要乳菇系贵州食用菌“黔菌”品牌重点打造的项目之一。黑柄乳菇（Lactarius gerardii ）属于外
生菌根菌，是一种美味的野生菌，市场价格不菲。为了解野生黑柄乳菇的菌根生境，给下一步研发黑柄乳菇的人工栽培技术奠定一些基础，以黑柄乳菇菌塘土壤为研究对象，采用高通量测序技术，结合生物信息学对来自于黔西南州不同地区的5个菌塘土壤的细菌组成及丰度进行研究。结果表明：黑柄乳菇菌塘菌分属于9个门27个属，菌塘中主要门有变形菌门、酸杆菌门、放线菌门和浮霉菌门；主要属有伯克氏菌属、缓生根瘤菌属、Gp 1、Gp 2、Gp 3、假单胞菌属和WPS-1genera incertae sedis 。
五个采集点菌塘中的菌丰富，多样性指数水准高，但因各个菌塘生境不同使得每个菌塘的优势菌属丰度有差异。关键词
美国吃人事件
黑柄乳菇；菌塘；高通量测序；细菌；菌落多样性
中图分类号：S646.1；S182
文献标志码：A
DOI ：10.19415/jki.1673-890x.2022.11.001
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样品预处理
称取200mg土样，放入灭菌的2mL离心管中，加入1mL PBS缓冲液混匀，10000rpm室温离心3min，弃上清，样品干燥待用。
1.3.2总DNA提取和检测
按照OMEGA试剂盒E.Z.N.ATM Mag-Bind Soil DNA Kit使用说明操作提取总DNA，1%琼脂糖凝胶检测DNA完整性。
1.3.3PCR扩增及测序
宏基因组测序类为细菌；扩增区域为细菌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多样性指数分析结果，可以反映所测微生物落的丰度和多样性。由表可知，五个菌塘土样测序共得到140903条有效序列，过滤后得到4385个有效OTU，样品的Coverage指数≥0.99，覆盖率高，能真实反映样品的测序情况。
表1菌塘土壤中细菌Alpha多样性指数分析
样品HZ LG PL WT XQ 有效序列数/条
25040
26755
老师辛勤付出的短句27646
36478
24984
OTU/个
804
956
906
832
887
Coverage
0.99
1.00
1.00
1.00
0.99
ACE
919.06
1019.77
980.94
884.52
980.30
Chao
922.65
1034.01
1007.02
894.19
998.78
Shannoneven
0.77
0.83
0.79
0.73
0.74
Simpson
0.02
8.7e-03
0.01
0.02
0.03
注：HZ-海子，LG-龙广，PL-平乐，WT-万屯，XQ-新桥。
样品的ACE指数和Chao指数反映落分布丰度（Communit richness），数值越大说明样本的菌分布丰度越好；Shannoneven指数越高或Simpson指数越低，说明采样点的菌多样性较丰富。从表1可以看出，龙广（LG）菌落较丰富、多样性较好。
图2为五个菌塘菌样品的稀释曲线，Shannon指数先呈直线上升，后趋向平稳，表明测序数据能较好地覆盖菌塘中菌的信息序列，曲线在平滑稳定后产出的OTU数趋于稳定，其中来自海子（HZ）的样品Shannon指数稍低，说明该采集地的细菌多样性低于其余四个菌塘。
2.2各样品中细菌类的相关性分析
图3（参见封三）为五个菌塘菌基于OTU
序列图2五个菌塘中菌的稀释曲线
所作的聚类树，可以看出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含量较高，通常认为这些菌增多可对被污染土壤进行修复。
表2五个样品中的优势细菌属及其相对丰度
样品PL WT XQ HZ LG
优势细菌属（相对丰度）
Bradyrhizobium（7.52%）、Gp1（7.02%）、Burkholderia（5.89%）、Gp2（3.87%）Bradyrhizobium（6.35%）、Rhizobium（5.60%）、Pseudomonas（5.55%）、WPS-1genera incertae sedis（4.54%）Burkholderia（9.53%）、Pseudomonas（5.67%）、Gp1（4.39%）、Bradyrhizobium（3.59%）
Gp1（10.24%）、Gp2（9.50%）、Bradyrhizobium（8.86%）、Burkholderia（6.54%）、Gp3（3.77%）
Bradyrhizobium（6.51%）、Gp2（4.72%）、Gp1（4.16%）、Gp3（3.58%）、Burkholderia（3.55%）
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.
（责任编辑：丁志祥）

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。