本报告介绍一项关于咸海沿岸土壤微生物群落对地球化学和矿物学变化响应的原创性研究。该研究由以中国科研人员为主体的国际团队完成，通讯作者为中国地质大学（北京）的姜红春博士，合作单位包括中国科学院新疆生态与地理研究所、中国地质大学（武汉）、德国莱布尼茨农业景观研究中心、乌兹别克斯坦国立大学、亚美尼亚埃里温州立大学、挪威卑尔根大学、美国太平洋西北国家实验室以及中山大学等多国科研机构。研究成果以《Onshore soil microbes and endophytes respond differently to geochemical and mineralogical changes in the Aral Sea》为题，于2021年发表在环境科学领域知名期刊《Science of the Total Environment》（第765卷）上。
本研究属于环境微生物学、地球化学与生态学的交叉领域。核心科学问题围绕湖泊持续干涸这一全球性环境变化背景下，新暴露的湖床土壤中微生物群落是如何形成和演替的。
研究背景与动机： 咸海自20世纪60年代以来，由于气候变化和水资源无序利用，经历了剧烈的萎缩和盐化，导致了严重的生态和环境灾难，如沙漠化、盐尘暴和生物多样性丧失。随着水位下降，大面积湖床暴露，形成了不同时期（从1970年代到2018年）暴露的土壤带。这些土壤继承了不同时期湖水的盐度和地球化学特征，形成了一个从远岸（暴露时间早、盐度低）到近岸（暴露时间晚、盐度高）的连续环境梯度。土壤微生物作为生态系统功能的关键驱动者，其在这种人为-自然双重压力下的形成机制和响应模式尚不清楚。同时，植物内生菌（endophytes）在帮助宿主应对环境胁迫（如盐度）方面扮演重要角色，但其在植物地上部分对土壤地球化学变化的响应也知之甚少。因此，本研究旨在填补这些知识空白，为理解湖泊干涸对沿岸土壤生态系统的影响提供科学依据。
研究目标： 具体研究目标包括：1）揭示咸海西岸不同暴露时期湖床土壤的地球化学和矿物学变化规律；2）阐明沿岸土壤中细菌和古菌群落的多样性、组成及预测功能沿环境梯度的变化；3）明确地球化学、矿物学和植物多样性等环境因子对土壤微生物群落构建的相对贡献；4）对比分析优势植物地上部分内生菌群落与土壤微生物群落在响应环境变化方面的差异；5）验证“土壤微生物比植物地上部分内生菌更容易受到盐度及其共变因素的影响”这一假设。
本研究是一项系统的野外采样与室内分析相结合的实证研究，主要包含以下五个关键步骤：
第一步：研究地点选择与样本采集。 研究区域位于咸海西岸一片长约1.5公里、宽约50米的暴露湖床上。根据卡拉卡尔帕克斯坦科学院自然研究所的资料，该区域自1970年代起逐步暴露，形成了六个不同暴露年代的土壤带（E48: 暴露约48年，至2018年；E38: 约38年；E28: 约28年；E18: 约18年；E9: 约9年；E1: 约1年）。此外，还采集了当前湖岸附近的沉积物（E0）作为对照。在每个土壤带，研究人员采集了靠近优势植物根部的土壤（深度3-10厘米）以及健康优势植物的地上部分组织样本。所有样本均按照分子生物学和地球化学分析的要求进行无菌处理、保存并运输至实验室。
第二步：地球化学与矿物学分析。 对采集的土壤和沉积物样本进行了一系列分析。首先，测量了总可溶性盐（Total Soluble Salts， TSS）和总有机碳（Total Organic Carbon， TOC）含量。TSS通过将土壤与去离子水混合后，利用离子色谱法测量上清液中的主要阴阳离子（Na⁺， K⁺， Ca²⁺， Mg²⁺， Cl⁻， SO₄²⁻）浓度并加和得到。TOC则使用总有机碳分析仪测定。其次，使用X射线衍射（XRD）技术对样本进行矿物学分析，鉴定并半定量了石英、方解石、白云石、石膏、石盐、粘土矿物等组分的相对含量。此外，研究还采用了一项先进的分析技术——飞行时间二次离子质谱（Time-of-Flight Secondary Ion Mass Spectrometry， ToF-SIMS），用于表征样本表面的分子特征。这是本研究方法学上的一个亮点。样本被冻干研磨后压制成薄片，使用TOF.SIMS5仪器进行正负离子谱采集。获得的质谱数据通过主成分分析（PCA）进行处理，以可视化不同样本间表面化学组成的差异，并识别导致差异的关键离子信号（如Ca⁺、Na⁺、SO₄²⁻等），这为理解矿物表面微环境和有机质分布提供了分子级别的信息。
第三步：DNA提取与高通量测序。 对于土壤样本，直接使用土壤DNA提取试剂盒提取总DNA。对于植物样本，在进行DNA提取前，实施了严格的组织表面消毒程序（包括用Tween 20、乙醇和次氯酸钠依次处理），并通过将最后一次漂洗液涂布平板培养来验证表面灭菌效果，确保所获DNA来自植物组织内部的内生菌而非附生菌。提取的DNA使用针对古菌和细菌16S rRNA基因V4区的通用引物进行PCR扩增。古菌引物为Arch519F/Arch915R，细菌引物为515F/806R。PCR产物经纯化、等量混合建库后，在Illumina HiSeq2500平台上进行双端250 bp测序。
第四步：测序数据分析与功能预测。 对获得的原始测序数据进行质量过滤、拼接，并在97%的相似度水平下聚类为操作分类单元（OTUs）。利用SILVA数据库对OTUs进行物种分类学注释。随后，将OTU表进行标准化（抽平）处理，用于计算α多样性指数（如观测OTU数、香农指数）。为了比较群落结构差异，计算了样本间的Bray-Curtis相异性矩阵，并进行了非度量多维尺度分析（NMDS）和相似性分析（ANOSIM）。为了探究微生物的潜在功能，研究采用了FAPROTAX（原核生物分类功能注释）管道，将16S rRNA基因序列映射到已知的功能类群数据库，预测了与碳、氮、硫循环等相关的代谢功能。
第五步：统计学分析与假设检验。 研究运用了多种统计方法验证环境因子与微生物群落的关系。通过Mantel检验评估了微生物群落结构与TSS、TOC、植物物种数以及各矿物成分之间的相关性。使用变异分区分析（Variation Partitioning Analysis， VPA）量化了地球化学因子（Env， 包括TSS和TOC）、植物多样性和矿物学这三组变量对微生物群落变异的解释比例。通过线性回归分析了特定微生物类群（纲水平）相对丰度与TSS的相关性，以及微生物群落相异性与石膏、方解石含量差异的相关性。所有分析均在R语言环境中完成。
结果一：土壤地球化学与矿物学的系统性梯度变化。 分析结果显示，从远岸（E48）到近岸（E1/E0），土壤的盐度（TSS）从0.5 g/L（淡水水平）急剧增加至71.3 g/L（超盐水平）。TOC含量则呈现相反趋势，在远岸土壤中较高。矿物组成发生了规律性演替：远岸老土壤（E48， E38）以碳酸盐矿物（方解石、白云石）为主；在中等盐度的E28和E18土壤中，蒸发岩矿物石膏开始出现并成为主要组分；而在最高盐度的E9和E1土壤中，石盐开始大量出现。ToF-SIMS的PCA结果进一步证实了这一化学梯度：远岸样品富含Ca⁺物种和含氮有机物信号，而近岸样品则富含NaCl簇和SO₄²⁻信号。这些结果清晰地描绘了咸海湖床土壤在持续干涸过程中，盐分积累和矿物沉淀的序列，为后续微生物研究提供了明确的环境背景。
结果二：土壤微生物多样性、组成及功能沿梯度的显著变化。 测序分析共获得902个古菌OTU和8753个细菌OTU。无论是古菌还是细菌，其α多样性（观测OTU数和香农指数）均随盐度增加而显著降低，符合高盐压力降低微生物多样性的生态学原理。群落组成也发生显著改变。例如，古菌中的盐杆菌纲（Halobacteria）相对丰度随盐度增加而上升，在近岸高盐土壤中占主导；而奇古菌门下的Nitrososphaeria等类群则在低盐土壤中更丰富。细菌中，γ-变形菌纲（Gammaproteobacteria）、拟杆菌纲（Bacteroidia）和红嗜热菌纲（Rhodothermia）的丰度向近岸增加，而放线菌纲（Actinobacteria）在远岸更丰富。FAPROTAX功能预测表明，潜在代谢功能也发生分异：远岸低盐土壤中，古菌的硝化作用和细菌的烃类降解等功能相对活跃；而近岸高盐土壤中，与硫酸盐还原、硫化合物呼吸相关的细菌功能更为突出。
结果三：矿物学是塑造土壤微生物群落的最主要因素。 Mantel检验表明，土壤微生物群落结构与TSS、大多数矿物成分以及植物物种数均存在显著相关性。VPA分析进一步揭示，在所有测量的环境变量中，矿物学对古菌和细菌群落变异的解释贡献最大（分别占可解释变异的很大部分），超过了盐度（TSS）和植物多样性的单独贡献。线性回归分析具体显示，古菌群落相异性与石膏含量差异显著相关，而细菌群落相异性则与石膏和方解石含量差异均显著相关。这一结果强调了在咸海这种特殊的地球化学背景下，矿物表面提供的微生境和养分条件对微生物群落构建起到了至关重要的选择作用。
结果四：植物内生菌群落对土壤环境变化不敏感。 对七个优势植物地上部分的内生菌分析显示，其群落多样性和组成与土壤的TSS、矿物学或植物物种数均无显著相关性。更重要的是，同一种植物（如Chenopodium album， Haloxylon ammodendron）即使生长在不同盐度的土壤带中，其地上部分的内生菌群落组成也高度相似。相反，在同一土壤带内，不同植物物种之间的内生菌群落则存在差异。这表明，植物地上部分的内生菌群落主要受宿主植物种类（植物基因型或生理特性）控制，而对土壤地球化学和矿物学梯度表现出很强的“缓冲”或“隔离”能力。这支持了研究最初的假设，即植物可能通过内部的脱盐和水分保持机制，维持了相对稳定的内生微环境，使得内生菌免受外部剧烈环境变化的直接影响。
本研究系统揭示了咸海持续干涸过程中，新暴露湖床土壤的地球化学与矿物学演化规律，及其对土壤微生物和植物内生菌群落截然不同的塑造作用。
核心结论： 1）沿岸土壤形成了一个盐度和矿物组成（从碳酸盐到蒸发岩）的连续梯度。2）土壤微生物（古菌和细菌）的多样性、组成和潜在功能沿此梯度发生显著变化，且群落差异与石膏、方解石等矿物含量密切相关，其中矿物学是驱动群落变异的最主要环境因子。3）与之形成鲜明对比的是，优势植物地上部分的内生菌群落不受土壤地球化学和矿物学变化的影响，其组成主要取决于宿主植物种类。
科学意义与应用价值： - 理论价值： 深化了对极端环境（干旱盐碱化）下微生物群落构建机制的理解，首次在咸海地区明确了矿物学在塑造沿岸土壤微生物群落中的主导作用。同时，通过对比土壤微生物与植物内生菌的响应差异，凸显了植物在缓冲环境胁迫、维持内部微生物稳态方面的重要生态功能，为研究植物-微生物互作提供了新视角。 - 应用价值： 研究成果有助于评估和预测湖泊萎缩、土地盐碱化等全球性环境问题对陆地生态系统，特别是土壤微生物生态功能的长期影响。对耐盐内生菌资源的认识，也为未来开发生物菌剂，用于盐碱地植被恢复和农业可持续发展提供了潜在的微生物种质资源信息和理论指导。
研究在讨论部分也指出，尽管测量了多种环境变量，但仍有超过60%的微生物群落变异未能得到解释。这暗示了除确定性环境筛选（如盐度、矿物）外，随机性过程（如扩散限制、随机出生死亡）以及未测量的生物相互作用（如物种间竞争、互利共生）可能在群落构建中同样扮演重要角色。这为未来更深入的研究指明了方向，例如需要结合网络分析、中性模型等工具来进一步解析群落组装机制。此外，研究也承认基于16S rRNA基因的功能预测（FAPROTAX）存在局限性，未来的宏基因组或宏转录组研究将能更准确地揭示微生物的真实功能活性。