梵净山典型植被类型的根系真菌群落结构及功能类群特征研究

刘敏　李龙　李海波等

关键词：植被类型；FUNGuild;根系真菌；功能类群；细根性状

中图分类号：S718.81 文献标志码：A 文章编号：1001-1498（2023）01-0047-12

梵净山地处热带、亚热带生物区系向温带生物区系的生态交错区，由于独特的地理位置、复杂的生物和生态演化过程、多样化的植物类群特征、众多地方特有和濒危动植物种类等特点，引起国内外不同领域学者的广泛关注，相关研究涉及到土壤理化性质、物种种类及多样性、珍稀濒危物种种群结构特征等多个方面。关于地下微生物群落的研究则集中在丛枝菌根真菌，对真菌群落结构及功能类群特征尚无系统报道。真菌特别是某些功能类群在促进有机质分解和养分循环、增强植物对矿物质营养元素的吸收、提高植物对生物和非生物胁迫的耐受性、维持生态系统的稳定性和生产力等方面发挥着重要作用。此外，真菌群落组成不仅受土壤理化性质、气候等非生物因素影响，其对植物种类、群落组成和多样性等生物因素也较为敏感，而植被类型可通过影响非生物和生物因素作用于真菌群落。尤为特殊的是，不同植被类型还会塑造特定的真菌功能类群：Tedersoo等通过对全世界范围内不同植被类型的土壤真菌群落进行分析后发现，外生菌根真菌在温带落叶林中较为丰富，热带潮湿雨林和热带干旱森林中植物病原菌更为多样；Liu等发现在中国东部森林地区，温带森林以外生菌根真菌为主，热带和亚热带森林中具有丰富的丛枝菌根真菌和植物病原菌；Sheng等研究了神农架山地森林生态系统中的土壤真菌功能类群，发现外生菌根真菌在落叶阔叶林、针阔混交林和针叶林中尤为丰富，植物病原菌一木腐生菌主要出现在常绿阔叶林，植物病原菌和寄生性真菌在山顶灌丛中具有较高相对丰度。菌在山顶灌丛中具有较高相对丰度。

梵净山拥有的典型植物区系过渡性及地球上同纬度保存完好的森林生态系统，为丰富且独特真菌种类的形成创造良好生存条件。此外，由于水热条件等差异造就了梵净山多样且分异明显的植被类型，为探讨植被类型与真菌群落间的关系提供了天然理想实验室。本研究以梵净山国家级自然保护区典型且具有代表性的常绿阔叶林、常绿落叶阔叶混交林及位于山顶附近的亚高山矮林3种森林植被类型为研究对象，分析真菌群落结构和功能类群特征与植被类型间的关系，探讨土壤性质、细根性状和叶性状对真菌群落的影响。研究结果可为后续探讨亚热带地区不同植被类型树种的环境适应机制及地上和地下部分的养分交换策略提供来自于地下真菌群落方面的数据支撑。

1材料与方法

1.1研究区概况及样本采集

研究区域位于梵净山国家级自然保护区（27°46′50″N～28°1′30″N，108°35′55″E～108°48′30″E）内，地处贵州铜仁市江口县、印江自治县、松桃自治县交界处，面积约为419 km²。亚热带季风湿润气候，年平均气温16℃，年平均降雨量1300 mm，平均相对湿度80%以上，土壤以山地黄壤和黄棕壤为主。植被类型丰富，具有代表性和面积较大的是常绿阔叶林和常绿落叶阔叶混交林，前者主要植物种类为栲（Castanopsisfargesii Franch.）、米槠（C.carlesii（Hemsl.）Hayata.）、巴东栎（Quercus engleriana Seem.）、山茶（Camellia japonica L.）等，后者植物种类包括水青冈（Fagus longipetiolata Seem.）、贵州青冈（Cyclobalanopsis argyrotricha （A.Camus）Chun et Y.T.Chang ex Y.C.Hsu et H.W. Jen）和扇叶槭（Acer flabellatum Rehd.）等。亚高山矮林是位于山顶附近的特殊植被类型，以吊钟花（Enkianthus quinqueflorus Lour.）、耳叶杜鹃（Rhododendron auriculatum Hemsl.）和长萼杜鹃（Rhododendron longicalyx Fang f.）等杜鹃花科植物为主。

在常绿阔叶林海拔600m处、常绿落叶阔叶混交林海拔1500 m处和亚高山矮林海拔2200 m处分别垂直于海拔高度设置1条100 m的样方线，每隔10 m设置一个采样点，中心象限法采集样本，即在每个采样点设置4个方向，选择距离采样点最近的2个树种（起测胸径为2cm），去除表层石块和落叶等杂物，从树干开始挖掘完整根系，抖落法收集根际土壤。由于部分采样点树木生长在坡度较大区域或岩石缝隙中，采集根系和土壤极为困难，仅48个样本用于后续分析（每种植被类型16个）。其中，常绿阔叶林中，栲5个，巴东栎4个，山茶3个，黄丹木姜子（Litsea elongata （Wall.ex Nees） Benth.et Hook.f.）2个，四照花（Cornus kousa subsp.chinensis（Osborn）Q.Y.Xiang）1个，柃木（Eurya japonica Thunb.）1个；常绿落叶阔叶混交林中，贵州青冈5个，水青冈4个，多脉青冈（Cyclobalanopsis multinervis W.C.Cheng＆T.Hong）3个，云山青冈（C.sessilifolia （Blume） Schottky）3个，青冈（C.glauca （Thunberg） Oersted）1个；亚高山矮林中，耳叶杜鹃8个，吊钟花6个，长萼杜鹃2个。同时，在每个植株上随机摘取20片向阳、完全展开、没有病虫害的叶片，用于叶性状指标分析。根系样本装在放有干冰的泡沫箱运回实验室，一部分保存在4℃冰箱用于分析根性状指标，另一部分保存在-70℃冰箱用于真菌测定。

1.2细根性状、叶性状和土壤理化性质分析

根系用自来水清洗掉附着土壤后，依据Pregitzer等根序分级法，按照序级进行分级，前两个序级归为吸收根并进行分离。EPSON Perfection V700（EPSON America Inc.）扫描细根样本，根系图像分析软件WinRhizo（Regent Instruments Inc.，Quebec，Canada）分析扫描的图片获得平均根直径/mm、根长/cm、根表面积／cm²和根体积／cm³。将扫描过的细根称量质量后，放入信封袋在70℃下烘48 h，再次称量质量得到根生物量／g。计算比根长（根长／生物量），（cm·g^-1）、根比表面积（根表面积／生物量），（cm²cm·g^-1）、根组织密度（生物量，根体积），（g·g^-3）和根干物质含量（根干质量，根鲜质量）／（g·g^-1）。

叶片放平整后利用Nikon Coolpi×5000数码相机（Nikon Corp.，Tokyo，Japan）拍照，ImageJ（Systat Software Inc.，Richmond，CA）软件计算叶面积／cm²。称重后，在70℃下烘72 h，再次称量质量得到叶生物量，g。计算比叶面积（叶面积／叶生物量），（cm²·g^-1）和叶干物质含量（叶干质量／叶鲜质量）／（g·g^-1）。

所有土壤、细根和叶样本均重复测定3次。土壤pH测定采用pH计法（土水比例1：2.5）；土壤、细根和叶碳、氮含量测定利用CN元素分析仪（Vario EL Ⅲ，Hanau.Germany）;土壤全磷测定采用HClO₄-H₂SO₄消煮一钼锑抗比色法，细根和叶磷含量测定采用H2S04-H202消煮一钼锑抗比色法。

1.3根系DNA提取、PCR扩增和测序

根系用蒸馏水冲洗10 min后，70%乙醇浸泡1min，3.3%次氯酸钠浸泡3min，75%乙醇浸泡30 s，无菌水反复冲洗后验证消毒效果。灭菌后根系用液氮冷冻研磨，DNeasy Plant Mini Kit（Qiagen Gmbh， Hilden， Germany）提取DNA，0.8%琼脂糖凝胶电泳检测DNA提取质量，紫外分光光度计（Thermo Fisher Scientific，MA，USA）准确定量DNA浓度。真菌rDNA-ITS保守区引物ITSIF12（5-GTGARTCATCGAATCTTTG-3）和ITS2 （5-TCCTCCGCTTATTGATATGC-3）进行PCR扩增。扩增体系包括12.5μL Phusion（R）Hot Start Flex 2X Master Mix（New England Biolabs，Ipswich，USA）、10 μmoKl·L^-1正向引物和反向引物各2.5 pL、Genomic DNA 50 ng，补充灭菌超纯水到25 pL。扩增条件：98℃预变性30 s；98℃变性10 s，54℃退火30 s、72℃延伸45 s，共35个循环；72℃延伸10 min。扩增产物利用2%琼脂糖凝胶电泳检测，AxyPrepPCR Cleanup Kit（Axygen Biosciences， CA， USA）纯化，Quant-iTTM PicoGreen^TMdsDNA Assay Kit（Thermo Fisher Scientific，MA，USA）定量浓度，Illumina MiSeq平台测序。

1.4生物信息学分析和真菌功能预测

DADA2（divisive amplicon denoisingalgorithm）进行序列拼接、去除错误碱基序列、去掉barcode和引物序列、过滤嵌合体后，Ribosomal Database Project （RDP）Naive Bayesian Classifier得到扩增序列变异（amplicon sequence variants，ASVs）。ASVs相当于以100%相似度聚类OTUs（operational taxonomic units），但与OTUs比，具有较高数据精确度与物种分辨率。在UNITE数据库中比对，得到序列注释信息，FUNGuildc比对和解析真菌营养型和功能类群，置信度为“极可能”（highly probable）和“很可能”（probable）等级的功能分组被用于后续分析。

1.5数据统计分析

混合效应模型探讨植被类型对土壤理化性质、细根性状、叶性状、和真菌α多样性指数的影响，植被类型为固定因子，树种为随机因子，并检验模型残差的正态性、齐性和独立性，若不满足则对数据进行转换。LEfSe（http：//huttenhower.sph.harv-ard.edu/lefse/）判别不同植被类型间真菌的显著差异类群。非度量多维标度（nonmetric multidimensional scaling，NMDS）排序和置换多元方差分析（permutational multivariate analysis of variance，PERMANOVA）对不同植被类型间真菌群落组成进行比较，Kruskal-Wallis test和Wilcoxon signed-rank tests判别植被类型间真菌Bray-Curtis相异性指数的高低。植被类型、树种、细根性状、叶性状和土壤性质对真菌α多样性指数影响的相对重要性基于ABT（aggregated boosted tree）分析，并计算其与真菌α多样性指数的Pearson相关性。降趋势对应分析（detrended correspondence analysis，DCA）检验真菌数据为单峰分布还是线性分布，DCA排序轴梯度长度最大值为2.657，小于3.0，选择线性排序模型中的冗余分析（redundancy analysis，RDA）探讨影响真菌群落的因子；执行RDA前，对所涉及解释变量进行前向选择。方差分解分析评估植被类型、土壤性质、细根性状和叶性状对真菌群落影响的大小。除特别说明外，统计分析均在R中完成。

2结果

2.1不同植被类型根系真菌群落特征概况

48个根系样本共得到5 103个ASVs，属于9门、193科和385个属（不包括unidentified fungi）。子囊菌门（Ascomycota）和担子菌门（Basidiomycota）为优势门，占全部ASVs数的89.93%。小蔓毛壳科（Herpotrichiellaceae）、球囊霉科（Glomeraceae）和柔膜菌科（Helotiaceae）为优势科，晶杯菌科（Hyaloscyphaceae）、Sebacinales_Group_A、被孢霉科（Mortierellaceae）、皮盘菌科（Dermateaceae）、红菇科（Russulaceae）、小菇科（Mycenaceae）、丛赤壳科（Nectriaceae）、水盘菌科（Vibrisseaceae）、锤舌菌科（Leotiaceae）、肉座菌科（Hypocreaceae）和革菌科（Thelephoraceae）相对丰度超过1%。常绿阔叶林具有最多ASVs数（2287个），子囊菌门为优势门，球囊霉科为优势科；亚高山矮林共得到ASVs数为1636个，子囊菌门为优势门，柔膜菌科和晶杯菌科为优势科；常绿落叶阔叶混交林ASVs数最少，为1608个，子囊菌门和担子菌门为优势门，球囊霉科和小蔓毛壳科为优势科。LEfSe分析（图1A）显示，常绿阔叶林显著差异类群为球囊菌门（Glomeromycota）、球囊菌纲（Glomeromycetes）、粪壳菌纲（ Sordariomycetes）、Pezizomycotina_Incertae_sedis纲、Pezizomycotina_lncertae_sedis目、粪壳菌目（Sordariales）、鸡油菌目（Cantharellales）、球囊霉目（Glomerales）、球囊霉科、丛赤壳科、Pezizomycotina_lncertae_sedis科和发菌科（Trichocomaceae），常绿落叶阔叶混交林中显著差异类群为革菌目（Thelephorales）、红菇目（Russulales）、革菌科和红菇科，锤舌菌纲（Leotiomycetes）、柔膜菌目（Helotiales）、蜡壳耳目（Sebacinales）、Sebacinales_Group_B科、柔膜菌科和锤舌菌科是亚高山矮林的显著差异类群。

2.2不同植被类型根系真菌功能类群分析

3种植被类型树种根系真菌在“极可能”和“很可能”水平上进行功能分组的ASVs数为2231个，共被分为病理营养型（pathotroph）、共生营养型（symbiotroph）、腐生营养型（saprotroph）、病理-共生营养型（pathotroph-symbiotroph）和病理-腐生营养型（pathotroph-saprotroph）5个营养型。其中，41.86%的ASVs被归为腐生营养型，其次是共生营养型（28.73%）和病理营养型（15.37%），其余2个营养型所占比例较低。在功能亚类水平上，undefined saprotroph相对丰度最高，其次是丛枝菌根（arbuscular mycorrhizal）真菌、外生菌根（ectomycorrhizal）真菌、植物病原菌（plant pathogen）、内生真菌（endophyte）、寄生性真菌（fungal parasite）和杜鹃花类菌根（ericoid mycorrhizal）真菌，其余功能亚类相对丰度较低。此外，丛枝菌根真菌仅出现在常绿阔叶林和常绿落叶阔叶混交林，且在常绿阔叶林（7.77%）的相对丰度高于常绿落叶阔叶混交林（7.30%）；外生菌根真菌在常绿落叶阔叶混交林（10.20%）的相对丰度明显高于常绿阔叶林（3.28%）和亚高山矮林（0.98%），内生真菌和杜鹃花类菌根真菌只出现在亚高山矮林。LEfSe分析（图1B）还发现，丛枝菌根真菌是常绿阔叶林的显著差异功能类群，外生菌根真菌是常绿落叶阔叶混交林的显著差异功能类群，亚高山矮林的显著差异功能类群为杜鹃花类菌根真菌。

2.3不同植被类型根系真菌群落结构分析

混合效应模型结果表明，植被类型能显著影响真菌丰富度和Shannon指数（P<0.05），对Simpson指数和Pielou指数无显著影响（P>0.05），常绿落叶林中真菌丰富度显著高于常绿落叶阔叶混交林和亚高山矮林，Shannon指数显著高于亚高山矮林（图2）。

NMDS分析结果显示真菌群落组成在植被类型间存在差异（图3A），PERMANOVA揭示了相同的结果（F=2.856，P=0.001）。此外，Bray-Curtis相异性指数（图3B）在3种植被类型间存在极显著差异（P<0.001），且在常绿落叶阔叶混交林中最高，常绿阔叶林次之，亚高山矮林最低，表明真菌群落组成在亚高山矮林样本间具有最高相似性，其次是常绿阔叶林，常绿落叶阔叶混交林的相似性最低，NMDS分析也印证了该结果。

2.4植被类型、树种、土壤性质、细根性状和叶性状与根系真菌群落的关系分析

以植被类型为固定因子，树种为随机因子构建的混合效应模型发现，植被类型可显著影响土壤pH、土壤全碳、土壤全氮、土壤全磷、比根长、根组织密度、平均根直径、根比表面积、根碳含量、根氮含量、叶面积和叶碳含量（P<0.05）。此外，常绿阔叶林和常绿落叶阔叶混交林土壤pH显著高于亚高山矮林，亚高山矮林土壤具有最高的全碳和全氮含量，常绿阔叶林具有最高的土壤全磷含量；比根长、根比表面积和根碳含量在亚高山矮林显著高于常绿阔叶林和常绿落叶阔叶混交林，常绿阔叶林具有最高的根组织密度、平均根直径和根氮含量；叶面积在常绿阔叶林和常绿落叶阔叶混交林显著高于亚高山矮林，叶碳含量在亚高山矮林显著高于常绿阔叶林和常绿落叶阔叶混交林（表1）。

ABT分析用来揭示植被类型、树种、土壤理化性质、细根性状和叶性状对真菌丰富度和多样性指数影响的相对重要性（图4）。结果表明：根碳含量是影响真菌丰富度的首要因子，相对重要性可达到24.47%（图4A）；影响Shannon指数、Simpson指数和Pielou指数的首要因子均为叶干物质含量，相对重要性分别为14.63%、18.08%和15.55%（图4B、图4C、图4D）。Pearson相关性分析结果指出，根碳含量与真菌丰富度（r=-0.626，P<0.001）、Shannon指数（r=-0.407，P=0.004）呈显著负相关，叶干物质含量与Shannon指数（r=-0.293，P=0.043）、Simpson指数（r=-0.591，P=0.002）、Pielou指数（r=-0.531，P=0.0259）呈显著负相关。

对包括植被类型和树种在内的22个解释变量进行前向选择后剩余7个，分别是植被类型、土壤pH、土壤全磷、根干物质含量、根碳含量、叶干物质含量和叶碳含量，与真菌群落构建RDA模型（图5A）。结果发现，植被类型、土壤pH、根碳含量和叶碳含量是影响真菌群落的显著因子（P<0.05）。方差分解分析进一步揭示植被类型、土壤性质（土壤pH、土壤全磷）、细根性状（根干物质含量、根碳含量）和叶性状（叶干物质含量、叶碳含量）对真菌群落的影响（图5B），发现植被类型（26.96%）是影响真菌群落的首要因子，其次是细根性状（20.10%）和土壤性质（13.53%），叶性状（3.92%）影响最小。

3讨论

3.1不同植被类型根系真菌和功能类群

本研究基于高通量测序技术探讨世界自然遗产地梵净山地区3种植被类型根系真菌群落，发现子囊菌门和担子菌门为优势真菌门，森林生态系统中枯枝落叶和枯死木是植物和微生物主要养分来源，子囊菌门和担子菌门的大多数真菌能促进纤维素和木质素等的分解和矿化从而调控生物地球化学循环。此外，丛枝菌根真菌为常绿阔叶林的显著差异功能类群，外生菌根真菌是常绿落叶阔叶混交林的显著差异功能类群，可能因为常绿阔叶林养分周转速率一般高于落叶阔叶林，不同菌根真菌采取的养分循环策略存在差异，丛枝菌根真菌通常在养分周转速率快的地区占优势，外生菌根真菌普遍存在于养分循环慢的区域。同时，本研究还发现，杜鹃花类菌根真菌主要出现在亚高山矮林，一方面梵净山地区亚高山矮林的优势植物主要以杜鹃花科为主，这为杜鹃花类菌根真菌的生存提供大量可供选择的宿主；另一方面亚高山矮林生长在山顶土层浅薄、土壤贫瘠、昼夜温差大、温度低的恶劣生境中，丰富和多样化的杜鹃花类菌根真菌可通过促进杜鹃花科植物对矿物质营养元素的吸收及产生具有抗氧化和水解作用的酶来帮助宿主植物获得复合有机物中的养分，从而增强杜鹃花科植物对山顶恶劣环境的适应能力。

3.2不同植被类型真菌群落结构的差异

常绿阔叶林中植物根系真菌具有最高的α多样性，可能原因如下：第一，常绿阔叶林分布区域的温度高于常绿落叶阔叶混交林和亚高山矮林，而较高的温度能够增强初级生产力及枯枝落叶物等的分解速率，为异养真菌生存提供更多资源；第二，亚热带森林中大部分树种为菌根依赖型树种，菌根真菌的存在能帮助宿主植物获取更多可利用土壤养分。虽然常绿阔叶林土壤磷含量高于其它两种植被类型，但在整个亚热带森林，磷是植物生长的重要限制元素，再加上常绿阔叶林具有最低的土壤氮含量。一方面某些菌根真菌的存在有助于溶解土壤中不能被植物吸收利用的磷酸盐和固定生物氮，另一方面宿主植物会倾向于选择更多的菌根真菌搭档来获得生长所必须的氮和磷；使得常绿阔叶林具有丰富且多样的真菌种类。此外，真菌群落组成在3种植被类型间存在差异，RDA结果也发现植被类型是影响真菌群落组成的显著因子。同在中国亚热带地区，Han等的研究表明真菌群落组成在常绿阔叶林、落叶阔叶林、针阔混交林、针叶林和次生灌木林存在显著差异。植被类型间不仅植物种类及群落组成、植物功能性状等生物因素存在差异，凋落物厚度及成分、土壤理化性质等非生物因素也有较大差别，从而可能直接或间接影响真菌群落组成。亚高山矮林真菌具有最低的Bray-Curtis相异性指数，是由于亚高山矮林所处区域植物种类（以杜鹃花科植物为主）及土壤、气候条件的独特性而孕育出高度相似的真菌群落。

3.3驱动根系真菌群落组成的主要影响因素

植物功能性状是植物长期适应外部环境变化的结果，近年来越来越多的研究关注根和叶功能性状与真菌群落的关系：Prada-Salcedo等指出根氮含量和根组织密度与真菌群落组成密切相关；Sweeney等发现根直径、根氮含量和比根长在驱动真菌群落组成中发挥着重要作用；而本研究发现根碳含量和叶碳含量能显著影响真菌群落组成，出现这些不同结果可能与宿主植物种类及采样部位的差异有关。当土壤养分受到限制时，宿主植物通常会将自身大量碳交易给菌根共生体，让更多菌根真菌搭档去竞争碳水化合物从而提高养分吸收效益，可能是根碳含量和叶碳含量显著影响真菌群落组成且根碳含量与真菌丰富度和Shannon指数呈显著负相关的原因。本研究还发现土壤pH是影响真菌群落组成的显著因子，先前同在森林生态系统中的研究也指出，土壤pH可以影响腐生真菌、丛枝菌根真菌以及外生菌根真菌，这是因为土壤pH会改变养分可利用性。此外，也有研究进一步探讨了土壤性质、根性状和叶性状对真菌群落的影响大小：Sweeney等指出与叶性状相比，根性状能决定丛枝菌根真菌和腐生真菌的丰富度和多样性；Dong等发现根系真菌多样性和群落组成主要受土壤性质（土壤pH、土壤碳磷比和氮磷比）的显著影响，而与叶性状无显著相关性。在本研究中，方差分解分析结果表明，细根性状对真菌群落的影响大于土壤性质和叶性状，可能原因：一方面植物细根（吸收根）是大多数真菌的栖息场所，较高比根长和根比表面积能为其生存提供更多空间；另一方面细根是某些真菌与宿主植物进行物质交换区域，植物可能会更加依赖于共生真菌来应对亚热带地区缺磷的土壤环境，使某些真菌功能类群（如丛枝菌根真菌、外生菌根真菌和杜鹃花类菌根真菌等）对细根性状变化较为敏感。对上述原因分析只是基于当前生态学知识的推测，后续研究会从植物种类与真菌群落结构的关系、不同菌根类型与宿主植物的养分获取途径等方面深入解析梵净山地区地上植物群落和地下真菌群落的关联机制。

4结论

世界自然遗产地梵净山3种典型植被类型的植物根系中，子囊菌门和担子菌门为优势真菌门，小蔓毛壳科、球囊霉科和柔膜菌科为优势真菌科，腐生营养型真菌所占比例较高。丛枝菌根真菌是常绿阔叶林的显著差异功能类群，外生菌根真菌是常绿落叶阔叶混交林的显著差异功能类群，而亚高山矮林的显著差异功能类群为杜鹃花类菌根真菌。植被类型、土壤pH、根碳含量和叶碳含量显著影响根系真菌群落组成，且细根性状对真菌群落的影响大于土壤性质和叶性状。研究结果不仅为丰富世界自然遗产地的物种多样性提供了数据支撑，还能够为深入探讨亚热带地区土壤·真菌·植物养分交易策略奠定基础。