沿高山沼泽湿地退化梯度丰富与稀有细菌及真菌亚群落的差异化生态响应——青海湖北岸土壤宏基因组证据

《Frontiers in Microbiology》：Contrasting ecological responses of abundant and rare bacterial and fungal subcommunities along an alpine marsh wetland degradation gradient

　　

引言：丰富与稀有微生物类群可能对湿地退化响应不同，但其对群落重组的相对贡献在高山沼泽湿地中仍未被充分解析。方法：研究人员在青海湖（Qinghai Lake）北岸利用15个土壤宏基因组（soil metagenomes）和9种土壤理化性质，调查轻度退化（LD）、中度退化（MD）和重度退化（SD）高山沼泽湿地。根据域特异性流行度（prevalence）和平均相对丰度，将细菌和真菌属划分为丰富、中间、稀有或低流行类群。沿退化梯度比较了样本水平丰富度和原始丰度贡献、α和β多样性、群落组成、分布广度、候选指示属和环境关联。结果：总检测细菌丰富度从LD的每样本346.2个属降至SD的287.0个属。稀有细菌在检测丰富度中的占比从45.3%降至40.4%，其原始丰度贡献从6.79%降至3.70%。丰富真菌占真菌总丰度的97.1–98.5%。稀有细菌丰富度随退化单调下降，而丰富真菌的Shannon多样性、Simpson多样性和Pielou's evenness下降。所有四个焦点亚群落在退化水平间组成差异显著，丰富真菌的PERMANOVA（置换多元方差分析）效应最大，稀有真菌最小。稀有部分表现出更高的基于发生率的β多样性和更窄的观测分布。土壤PC1和PC2分别与丰富细菌、稀有细菌和丰富真菌群落变异的51.6%、41.4%和55.1%相关，但仅与稀有真菌群落变异的0.6%相关。讨论：将微生物群落划分为丰富和稀有部分，揭示了沿湿地退化梯度显著的丰度依赖和域依赖性重组。稀有细菌表现出丰富度和丰度贡献的明显损失，而丰富真菌表现出更强的群落内多样性下降和更强的环境结构化。这些对比性响应在全群落汇总中被掩盖，凸显了在评估微生物对高山沼泽湿地退化的响应时明确分离丰富和稀有类群的价值。

论文解读：沿高山沼泽湿地退化梯度丰富与稀有细菌及真菌亚群落的差异化生态响应

高山沼泽湿地是青藏高原高海拔生态系统的重要组成部分，提供水源调节、碳固存、养分循环和生物多样性保护等生态系统服务。由于低温、高土壤湿度和有机质分解缓慢，高山湿地积累大量土壤有机碳，是山地重要碳库。但青藏高原高山沼泽湿地对气候变化和人为干扰高度脆弱，许多湿地已显著退化。退化常伴随植被退化、水文改变、土壤湿度下降、养分失衡和土壤理化条件变化，从而改变地下生境和生态系统过程。土壤微生物调节有机质分解、养分转化和碳循环，因此退化诱导的植被和土壤变化会通过资源有效性和环境过滤影响微生物群落结构。理解土壤微生物如何响应高山沼泽湿地退化，对认识环境恶化的生态后果和预测高海拔湿地生态系统稳定性至关重要。土壤微生物群落丰度分布高度不均，少数丰富类群（abundant taxa）与大量低丰度稀有类群（rare taxa）共存，构成微生物稀有生物圈。丰富类群通常占群落生物量和代谢活动的主要部分，稀有类群则贡献大量分类多样性，并可能作为系统发育和生态变异储库。丰富与稀有类群代表对比性生态策略：丰富类群局部丰度高、出现频繁、环境分布较广；稀有类群丰度低、分布较窄、生境特化较强，部分可通过休眠、扩散或功能冗余在不利条件下持续存在。尽管微生物对生态系统退化的响应已有大量研究，但多数集中于总微生物群落，较少同时解析不同丰度类群，尤其缺乏高山沼泽湿地中丰富与稀有细菌和真菌亚群落的宏基因组证据。因此，研究人员利用青海湖流域轻度退化（LD）、中度退化（MD）和重度退化（SD）空间代替时间梯度，结合15个土壤宏基因组和土壤理化测量，检验丰富、稀有、中间和低流行类群的丰富度与原始丰度贡献、α多样性与组成、基于发生率的β多样性组分、观测分布广度、退化相关指示属以及组成变异与环境梯度的关联。

关键技术方法：研究人员在中国青海湖（Qinghai Lake）北岸刚察县建立LD、MD、SD高山沼泽湿地空间梯度，2025年6月采集15个0–20 cm复合土壤样本，每退化级别5个样地，每样地3个土芯混合，并测定9种土壤理化性质。对土壤DNA进行宏基因组测序（DNBSEQ-T7，2×150 bp），经fastp质控、MEGAHIT组装、Prodigal基因预测、MMseqs2去冗余、CoverM以TPM（transcripts per million，每百万转录本）定量，DIAMOND比对细菌GTDB和真菌NR数据库，CAT_pack分类注释。属水平相对丰度按域分别定义丰富类群（平均相对丰度≥0.10%）、稀有类群（<0.05%）、中间类群（0.05–0.10%）和低流行类群（流行度<20%）。采用α多样性、PCoA、PERMANOVA、PERMDISP、Jaccard β多样性、Levins生态位宽度、IndVal指示属、PCA、dbRDA和Mantel分析。

研究结果：
3.1 丰富与稀有微生物类群的组成和丰度贡献：通过比较所有丰度类别，检测细菌总丰富度从LD的346.2属/样本降至MD的297.6和SD的287.0，真菌从53.0降至50.2和47.8。细菌中，稀有属在LD占检测丰富度45.3%和原始丰度6.79%，在SD分别为40.4%和3.70%；丰富细菌丰度贡献从84.81%升至91.07%，丰富度贡献从31.8%升至38.9%；低流行细菌丰富度贡献从5.3%降至1.7%。真菌中，丰富属占检测丰富度68.5–72.4%和原始丰度97.06–98.46%，稀有属约贡献8%丰富度和低于0.7%丰度。结果表明，退化相关损失集中于低丰度细菌丰富度，而群落丰度仍由丰富属主导。
3.2 丰富与稀有亚群落的α多样性响应：通过Kruskal–Wallis、Spearman和Dunn检验，丰富细菌丰富度在退化水平间有差异（q=0.037），但无单调趋势，LD低于MD，SD居中；稀有细菌丰富度从LD的157.6属/样本单调降至SD的116.2（趋势q=0.028），但三组Kruskal–Wallis校正后不显著。丰富真菌的Shannon多样性、Simpson多样性和Pielou's均匀度沿梯度下降（趋势q=0.028–0.029；Kruskal–Wallis q=0.037），Dunn检验显示LD与MD和SD分离，MD与SD无差异。稀有真菌α多样性指标经FDR校正后不显著。支持的模式是稀有细菌丰富度单调下降，以及丰富真菌多样性和均匀度协同下降。
3.3 丰富与稀有亚群落的组成重组：通过PCoA和PERMANOVA，丰富细菌组成在LD、MD、SD间差异显著（R²=0.730，q<0.001），所有两两比较显著。稀有细菌总体组成显著重组（R²=0.637，q<0.001），但主要体现为LD与MD和SD分离，MD与SD不显著。丰富真菌组成差异最强（R²=0.823，q<0.001），LD与MD和SD显著分离，MD与SD较相似。稀有真菌总体组成差异显著但解释比例较低（R²=0.235，q=0.020）。PERMDISP经FDR校正后四个亚群落均不显著，表明组成差异主要反映组中心偏移而非组内离散度差异。总体上，丰富细菌和真菌组成响应最强，稀有真菌响应较弱。
3.4 β多样性组分和生态位宽度：基于存在–缺席矩阵的Jaccard相异度，丰富细菌总Jaccard相异度在LD、MD、SD接近零（0.022、0.004、0.016），嵌套组分绝对量仅0.019、0.004、0.012，说明组内组成一致性较强而非嵌套驱动重组。稀有细菌总相异度较高且相似（0.282、0.274、0.282），周转组分从LD的0.075增至MD的0.129和SD的0.111，嵌套相关组分分别为0.208、0.145、0.171。丰富真菌总Jaccard相异度为0.278、0.192、0.276。稀有真菌样本间相异度很高（0.849、0.869、0.886），绝对周转组分为0.698、0.688、0.775。丰富属的Levins生态位宽度在细菌和真菌中均大于稀有属（均q<0.001），对应其在15个样本中观测分布更广。由于检测概率随丰度下降，分布广度和基于发生率的β多样性需在此限制下解释。
3.5 跨退化水平的指示属：IndVal分析识别出70个经FDR校正后仍显著的退化水平指示属。LD含最多指示属，尤其稀有细菌中有37个LD相关指示属；LD还含10个丰富细菌和10个丰富真菌指示属。丰富细菌LD指示属代表包括Rubrivivax、Candidatus Deferrimicrobium、Hydrogenophaga、Thiobacillus、Anaerolinea和Gallionella，丰富真菌包括Cryphonectria、Diaporthe、Xylogone、Bimuria、Periconia和Pyrenophora。MD含较少显著指示属，包括3个稀有细菌属和丰富真菌Glomus、Ambispora。SD含4个丰富细菌和4个稀有细菌指示属，代表丰富细菌包括Sphingomonas、Candidatus Acidiferrum、Sphingomicrobium和Gillisia；未检测到FDR显著的真菌SD指示属，也未鉴定到稀有真菌指示属。指示属在LD、尤其稀有细菌中的集中表明，退化不仅伴随丰富类群替换，也伴随许多与轻退化条件相关的稀有类群下降或再分布。
3.6 丰富与稀有亚群落组成的环境解释：土壤理化性质PCA显示PC1和PC2分别解释53.5%和25.1%，合计78.6%。LD样本与较高土壤湿度和较高N、P、K有效性相关，SD样本趋向较高pH，MD居中。Mantel分析显示，丰富细菌、稀有细菌和丰富真菌与多个土壤性质显著关联，尤其pH、有效磷、有效钾、土壤湿度和土壤温度；稀有真菌与所测土壤变量无显著关联。dbRDA显示，土壤PC1和PC2显著解释丰富细菌51.6%、稀有细菌41.4%和丰富真菌55.1%的组成变异（FDR校正后均q<0.001），稀有真菌仅0.6%且不显著（q=0.438）。稀有细菌不可解释比例大于丰富细菌。总体表明，丰富细菌和真菌更紧密关联土壤退化梯度，稀有细菌关联中等，稀有真菌的大部分组成变异不能由所测土壤性质解释。

讨论部分总结：4.1 丰富与稀有亚群落的不对称响应：全群落评估常掩盖退化响应差异。稀有类群占细菌分类多样性较大比例，但丰富类群占细菌和真菌总丰度多数，这种分类丰富度与相对丰度不对称符合微生物稀有生物圈概念。退化未导致统一生物多样性下降：稀有细菌丰富度显著下降，丰富真菌多样性和均匀度下降，丰富细菌与稀有真菌相对稳定。稀有类群不应普遍视为更脆弱，特化不等于低持久性，部分稀有微生物可通过休眠、时间波动或功能冗余维持长期稳定。丰富类群生态位较宽、扩散机会较多，稀有类群分布较窄、生境关联较强，因此退化通过丰度类别间选择性再分布而非简单总多样性降低来重塑微生物多样性。
4.2 群落周转和生态特化：稀有细菌和真菌生态位宽度较窄，稀有真菌组内相异度高且由周转主导。丰富类群常占据更广生态范围，稀有类群可能与特定环境背景关联更强。但丰度类别不应解释为固定生态策略。稀有真菌高周转可能反映异质退化条件下的组成变异，而非退化单独驱动稀有类群替换。土壤环境梯度解释了丰富细菌、稀有细菌和丰富真菌的大量变异，但稀有真菌组成与所测环境变量关联有限，可能响应细尺度生境异质性、植物源资源或未测量因素。指示属集中在稀有细菌群落，尤其轻退化湿地，表明稀有类群可能提供生境转换的敏感生态信号，但应解释为生态关联而非退化功能驱动。丰富类群可能通过较广环境分布促进群落持久性，稀有类群通过局域发生和较高替换率贡献群落变异性。
4.3 细菌和真菌的不同环境过滤：土壤pH、湿度、有效养分和温度共同刻画退化相关环境梯度，丰富细菌、稀有细菌和丰富真菌与这些环境变化关联强于稀有真菌。pH是微生物组成的重要决定因素，细菌常比真菌表现出更强pH响应。湿地退化中土壤湿度和养分有效性下降可改变底物扩散、氧条件和资源可及性，对不同生态策略类群施加选择压力。丰富类群可能跟踪广尺度退化生境变化，稀有真菌更依赖局域微生境或未测生态因素。细菌和真菌资源获取策略与生态位不同，细菌因世代短和依赖易利用底物而对无机养分和理化条件响应较快，真菌可通过腐生和共生策略占据较广资源生态位。研究中土壤有机碳、总碳和无机碳未随退化显著变化，但碳可利用性仍可能通过底物质量、植物输入和碳可及性影响微生物组织，整体碳库测量未完全反映。因此，退化相关微生物重组涉及亚群落间差异环境过滤，而非单一主导环境驱动。
4.4 退化和恢复监测的生态意义：稀有细菌群落含大量退化相关指示属，尤其在轻退化条件，表明稀有类群可捕捉优势微生物类群未完全反映的细微生境转换。微生物群落可作为生态系统状态指标，低丰度类群可能对时间和环境变化贡献不成比例，是生态响应储库。但稀有类群应解释为敏感生态指标而非退化直接驱动，因为其分布受生境条件、扩散过程和种群动态共同影响。将稀有微生物部分与环境测量结合，可能提高退化监测敏感性，尤其在微生物重组先于生态系统水平性质可检测变化时。仅基于丰富类群评估恢复可能低估群落重组复杂性。丰富类群可能代表主导生态功能，稀有类群提供群落灵活性和微生物多样性维持信息。恢复评估高山沼泽湿地时应同时考虑丰富和稀有微生物库，以更完整理解微生物稳定性和恢复潜力。

研究结论翻译：高山沼泽湿地退化以亚群落特异性方式重组了土壤微生物群落。丰富类群与土壤环境梯度的关联更强，而稀有类群表现出更大的空间周转，并且其组成变异中由已测量土壤变量解释的比例更小。退化相关指示属在稀有细菌群落中集中，以及稀有真菌群落的显著周转，进一步表明不同稀有微生物部分可提供湿地退化的不同生态信号。这些对比性响应表明，仅进行全群落分析可能掩盖丰富与稀有类群在退化响应、生态位特征和群落维持机制方面的重要差异。