硫酸富集岩溶含水层中不同水文条件下地下水微生物群落结构

《Frontiers in Microbiology》：Groundwater microbial community structure under contrasting hydrological conditions in a sulfate-rich karst aquifer

　　

岩溶含水层是动态的地下环境，水文波动与氧化还原（redox）变异性可强烈影响地下水微生物组及相关生物地球化学过程。硫酸富集非热液岩溶系统中的微生物硫循环仍知之甚少，尤其缺乏将水化学梯度与微生物代谢潜力联系的基因组解析（genome-resolved）视角。研究人员以巴西东北部Irecê盆地Salitre组硫酸富集岩溶含水层为对象，整合水化学分析、16S rRNA基因测序与基因组解析宏基因组学，表征地下水微生物群落，并评估水化学变异性与硫相关代谢潜力的关系。样品表现出显著水化学异质性，包括氧化–还原电位（ORP）、溶解氧（DO）、硫酸盐浓度及溶解硫化物积累的差异。10个样品的扩增子测序均获得超过10,000条高质量读数，群落组成在环境与对比水文条件间差异显著，显示明显空间与生境水平异质性。从3个样品恢复185个去冗余中高质量宏基因组组装基因组（MAGs），多数代表潜在新谱系。M00596相关硫转化模块在3个宏基因组样品的MAGs中均检出，而完整SOX模块仅限于相对更氧化的地下水。若干MAGs同时含M00596与SOX模块，表明同一基因组内可能存在多条硫转化途径。结果表明，硫循环是该硫酸富集岩溶含水层地下水生物地球化学的核心组分，水化学异质性与微生物群落结构及硫循环代谢潜力变化相关。硫酸富集碳酸盐岩性、异质氧化还原条件与多样硫循环微生物类群共同凸显硫转化对该含水层生态与生物地球化学功能的重要意义。

论文解读：《Frontiers in Microbiology》刊发的硫酸富集岩溶含水层微生物研究

研究背景与立项依据

岩溶含水层是全球重要地下水储层，其水文、地质与生物地球化学过程强烈耦合。地下水穿过碳酸盐岩形成物理化学条件高度异质的环境，常呈寡营养与能量受限特征，适合具有精简基因组和依赖共生代谢的微生物类群生存。候选门辐射（Candidate Phyla Radiation, CPR）与DPANN古菌在地下环境中常见，但其对硫循环的贡献仍不清楚。硫循环是地下环境主要生物地球化学过程之一，氧化与还原条件交替可支持硫酸盐还原与硫氧化微生物共存。已有研究多关注硫化洞穴与硫酸驱动溶蚀过程，而对缺乏开阔洞穴网络的岩溶地下水系统，在水化学变异、氧化还原波动与基因组解析代谢潜力交叉层面仍研究不足。因此，研究人员在巴西东北部Irecê盆地南部Salitre组新元古代灰岩岩溶含水层开展研究，以明确空间水化学异质性、对比水文条件与微生物群落及硫循环代谢潜力的关系。

研究内容与总体结论

研究人员整合水化学、16S rRNA基因扩增子测序与基因组解析宏基因组学，发现该含水层水化学高度异质，微生物群落具显著空间与生境水平差异；从三个代表性样品重建185个中高质量MAGs，其中大量为潜在新谱系。M00596相关硫转化模块广泛存在，完整SOX（sulfur-oxidation）模块仅出现于较氧化水样，部分MAGs共编码还原与氧化相关模块。研究证明硫循环是该含水层生物地球化学核心过程，水化学异质性驱动微生物结构与硫代谢潜力变化，对理解非热液硫酸富集岩溶系统生态功能具有重要意义。

主要技术方法概述

研究人员以巴西Irecê盆地Salitre组岩溶含水层为样本队列，在湿季与干季采集洞穴涌水、复苏带及社区/农用井水；对10个DNA足量样品做16S rRNA V4–V5区扩增子测序，对3个代表对比水化学条件的样品做鸟枪法宏基因组测序。生物信息学上采用DADA2推导扩增子序列变体（ASV），用MEGAHIT独立组装、anvi’o与CONCOCT分箱、CheckM/CheckM2评估质量、dRep去冗余、GTDB-Tk定类，并基于KEGG模块注释重建碳、氮、硫代谢潜力。

研究结果

3.1 样品间水化学变异性

通过现场探头与实验室离子分析，研究人员发现pH近中性至弱碱，DO从2.10%到65.40%，ORP从-233.60到81.60 mV，硫酸盐从2.91到1406.67 mg/L，溶解H₂S仅出现于干季部分样品。湿季更氧化、硫化物低于检测限，干季还原增强、H₂S累积，显示水文期与水化学梯度耦合。

3.2 基于16S rRNA扩增子的微生物群落组成

ASV分析显示样品间丰富度90–874不等。α多样性在湿干季间无显著差异；Bray–Curtis的PCoA显示水文期解释约16.2%群落变异，PERMANOVA显著但全模型不显著，说明空间与生境因素同样重要。dbRDA未达显著，但湿季ASVs更富含Proteobacteria、Verrucomicrobiota、Patescibacteria，干季BI-262被Halothiobacillaceae主导达90.6%，证明群落受环境而非单一季节控制。

3.3 基因组解析微生物多样性

三个宏基因组共恢复226个中高质量MAGs，95%平均核苷酸一致性去冗余后为185个；细菌占83.8%，古菌16.2%，仅约1.1%可归已知物种。Patescibacteriota等CPR类群在BI-04与BI-28占优势，Pseudomonadota、Desulfobacterota、Nitrospirota在多样品出现，显示岩溶地下水具高分类新颖度与CPR依赖型结构。

3.4 MAGs的功能潜力

KEGG模块显示中心碳代谢常不完整，BI-06代谢功能较多样，BI-28功能谱受限且偏硫酸盐还原。基因组解析结果反映氧化与还原生境中微生物代谢广度不同，硫相关途径是研究焦点。

3.5 硫循环功能类群

3.5.1 异化硫代谢的基因组潜力

以KEGG M00596模块（sat/aprAB/dsrAB）评估，BI-04有4个MAG完整编码，BI-06有10个且含Desulfobulbia等，BI-28仅2个；因DsrAB可还原也可氧化，完整M00596仅示硫转化潜力而非确定硫酸盐还原。

3.5.2 通过SOX系统的硫代硫酸盐氧化

M00595模块（soxA/B/C/D/X/Z）仅存于BI-04与BI-06；BI-28无完整SOX，与其低ORP、高H₂S及较低硫酸盐一致，说明氧化界面功能在还原样中缺失或未被恢复。

3.5.3 DSR与SOX模块共编码

六个MAGs（BI-04两个、BI-06四个）同时完整编码M00596与M00595；Thiothrix属MAG更宜解读为多硫氧化系统共存。BI-28无共编码类群，反映其还原生境中硫氧化基因组潜力下降。

3.6 CPR与硫循环MAGs的功能分化

对比Patescibacteriota（89个）与Pseudomonadota/Nitrospirota/Desulfobacterota（28个），CPR基因组呼吸、脂类、维生素与氨基酸合成模块极贫，硫循环完整模块为0%；代谢多样MAGs中呼吸75%、硫代谢11.8%。结果说明CPR类群代谢流线化、生态上依赖其他微生物，而硫循环由少数代谢完整类群承担。

讨论与结论翻译

讨论指出，湿季样品通常更氧化、硫酸盐可检、H₂S低于检测；干季BI-28与BI-262氧化还原电位更低、H₂S可测、硫酸盐略低，SO₄^2-/H₂S比跨约四个数量级。补给动态中，湿季大气水入渗带来富氧水，干季补给减少与滞留时间延长促离子累积和氧耗。空间异质性使同季样品群落仍分异，说明岩溶水道、裂隙与栖位共同塑造微生物结构。

研究结论：在Irecê盆地硫酸富集岩溶含水层中，水化学与氧化还原异质性显著；16S与MAGs共同表明微生物群落具高分类新颖度与空间分层；M00596广泛、SOX限于氧化样、部分基因组共载还原与氧化硫模块；CPR类群代谢极度精简，硫循环由Pseudomonadota、Nitrospirota与Desulfobacterota等代谢完整类群支撑。该工作从基因组解析视角证实硫转化是此类非热液岩溶地下水生物地球化学中枢，对区域元素循环、碳酸盐岩溶蚀及地下水管理具有理论与应用意义。