湿地与森林、海洋并称为地球三大生态系统，其在调节生物多样性和保护自然环境等方面具有不可替代的功能与效益。西藏作为地球上独特的环境地域单元，孕育了非常独特的高寒生态环境，有着丰富的湿地资源。根据全国湿地资源第二次调查结果，西藏现有河流、湖泊、沼泽等各类湿地6.53 × 10⁴ km²，湿地率（即湿地面积占国土面积的比率）达5.31%，高居全国第二位^[1]。西藏独特的自然地理环境条件，也孕育了丰富的湿地生物多样性。近几十年来，西藏高原湿地的生物多样性受气候变化和人为活动等因素的威胁越来越严重，其生态状况也引起了世界范围的高度关注^[2]。

植物是湿地生态系统的重要组成部分，也是其主要的初级生产者，对维持生态系统结构和功能的稳定有重要作用^[3]。植物多样性是评价湿地生态服务价值和生态系统稳定性的重要指标，与土壤、海拔、气候等环境因子密切相关。一方面，环境因子可以影响植物的生长与群落分布，如有研究表明土壤养分含量直接制约植被的生长发育^[4]。海拔也被认为是决定植被群落分布的重要因素之一^[5]，二者的关系目前主要有两种理论。第一种是“中间高度膨胀理论”^{[6, 7]}，即植物多样性在中海拔位置达到峰值，而在高海拔和低海拔地区相对较低。第二种是“负相关理论”^[8]，即植物多样性随着海拔的增加而降低。另一方面，越来越多的研究表明植物的凋落物和根系分泌物可显著改变土壤质地和肥力，进而可能影响植物的多样性^{[9, 10]}。有研究发现，西藏地区水生植物物种丰富度和系统发育多样性随海拔的增加而显著降低^[11]。也有研究表明，西藏高寒沼泽湿地植物多样性的变异受到多种环境因子相互作用的影响^[12]。目前对西藏湿地植物的研究大多集中在植物物种调查^{[13, 14]}和保护^[15]上，而对湿地植物与高寒环境之间关系的探究相对较少^{[16, 17]}，因此，探讨西藏湿地植物多样性与环境因子之间关系具有重要意义。

土壤微生物作为湿地生态系统的分解者，在能量流动、物质循环、生态系统演替以及生物多样性维持等方面发挥重要的生态功能。细菌是最常见和最丰富的土壤微生物类群，驱动着绝大多数生物地球化学循环过程，如土壤有机物矿化、硝化、反硝化等^{[2, 18]}。土壤细菌多样性与生态系统稳定性及土壤肥力有关，且能在一定程度上指示土壤养分状况和土壤质量。研究表明，土壤性质^[19]、气候因素^{[20, 21]}、地理距离^[22]、植物多样性^[23]等因子能显著影响土壤细菌多样性。如在高山草地和森林土壤中，pH是土壤细菌多样性形成和维持的主要驱动因子^{[24, 25]}；而在石灰质沙漠土壤中，土壤营养盐含量显著改变细菌群落的组成和结构^[26]。李明家等^[27]对西藏横断山区怒江和澜沧江两个流域进行了研究，发现气候和环境因子是影响细菌多样性的重要因素。甄莉等^[28]探究了西藏热泉沉积物中硫氧化细菌的多样性，发现溶解有机碳含量显著影响硫氧化菌群的多样性。此外，植被特征差异也能不同程度地解释土壤微生物的多样性。如在旱地生态系统中，植物丰富度、盖度以及叶片面积等功能性状可以解释土壤细菌多样性的变异^[29]。近年来，土壤微生物多样性的相关研究已经引起了国内外的广泛关注，但目前对西藏高寒环境下湿地土壤细菌多样性的研究还较少。

基于此，本研究选取26个具有代表性的西藏高原湿地为对象，探讨湿地植物和土壤细菌多样性，明确影响植物多样性与土壤细菌多样性的关键因子，试图回答“高寒环境下，哪些环境因子能够解释湿地植物和土壤细菌多样性的空间变化？”这一科学问题。研究结果旨在为了解高寒湿地生态系统的植物、微生物和土壤的相互关系，保护湿地生态系统生物多样性提供理论支撑。

1.   材料与方法

1.1   研究地点

于2020年8月在西藏选取26个相互隔离的典型湿地（28.87° ～ 33.44°N，79.77° ～ 91.88°E）开展植被调查与土壤采样工作（图1）。所选湿地的海拔范围为3974 ～ 5206 m，年平均气温（Mean annual temperature, MAT）为−2.72℃ ～ 6.25℃，年平均降水量（Mean annual precipitation, MAP）为41 ～ 429 mm。所选湿地8月的平均降水量（Average precipitation）为67 mm，平均气温（Average temperature）为10.2℃。

图  1  西藏湿地采样位置图

Figure  1.  Locations of 26 Tibetan Plateau wetlands

下载: 全尺寸图片 幻灯片

依据典型性原则，在每个湿地内选择能够代表整个植被和土壤等特征的地段设置1个样地。记录样地位置（经度、纬度、海拔）及周围环境特点（放牧干扰程度、水源补给方式等）等信息。在样地内随机布设5 ～ 15个样方，样方面积为1 m²（1 m × 1 m），记录样方内的植物种类、个体数（株 / 丛数）和盖度。随机选取其中1个样方，样方内出现的所有植物各取10个叶片用信封保存，当天测定生物量鲜重和叶面积，并带回实验室测干重、碳氮磷含量等植物功能性状。在该样方内取10 g表层根际土放入冻存管中，带回实验室测定细菌多样性，另取150 g根际土放入自封袋中，带回实验室检测土壤理化性质。

1.2   植物多样性与功能性状测定

根据样方调查数据，计算植物的物种丰富度（Richness）、香农指数（Shannon index）等多样性参数。采样结束后，将植物地上部分与地下部分分开，吸干表面水分后测生物量鲜重。叶片在叶面积板上展开并拍照，使用Image J软件测叶面积（Leaf area, LA）。植物样品带回实验室放置于烘箱中，105℃杀青15 min后，立即降低烘箱温度至80℃，48 h后称量生物量干重，并计算干物质含量（Leaf dry matter content, LDMC）和比叶面积（Specific leaf area, SLA）。烘干后的植物样品用球磨仪粉碎，使用元素分析仪（Elementar Analysensysteme, Vario MACRO cube, 美国）测定叶片总碳（TC）和总氮（TN）含量，使用等离子体原子发射光谱仪（Perkin Elmer, Optima 8000 DV, 美国）测定叶片总磷（TP）含量。根据样方内各物种的相对多度（个体数），计算群落水平上各功能性状的群落加权平均数（Community weighted mean, CWM）。

1.3   土壤细菌多样性的测定

基于16S rRNA基因，采用扩增子高通量测序技术测定土壤的细菌多样性。DNA提取步骤如下：称取0.5 g土壤样品，分别加入PB（pH值 8.0）溶液和TNS（pH值 8.0）溶液，在样品破碎仪（Fastprep-24, 美国）下进行破碎，并离心得到上清液，然后分别用试剂PCI（苯酚∶氯仿∶异戊醇 = 25∶24∶1）、CI（氯仿∶异戊醇 = 24∶1）进行提取，离心得到的沉淀加入TE（pH值 8.0）溶液进行溶解，即DNA溶液。DNA溶液提纯后送至上海派森诺公司进行高通量测序。细菌群落DNA片段的双端（Paired-end）测序在Illumina MiSeq 300PE平台进行，对获得的序列使用DADA2方法进行去引物、质量过滤、去噪、拼接和去嵌合体等操作^[30]，获得的序列即为操作分类单元（Operational taxonomic units, OTUs）代表序列。使用QIIME 2软件中的“qiime feature-table rarefy”功能，抽平深度设为最低样本序列量的97%，获得抽平OTUs表。通过QIIME 2软件计算土壤细菌的Chao1和香农等多样性指数^{[31, 32]}。

1.4   土壤理化性质的测定

新鲜土壤样品使用烘干法测量含水率（Moisture）。土样经自然风干、碾碎、去除细根等预处理后，分别过10目和100目筛备用。过10目筛的土样采用5%的H₂O₂去除有机质，然后使用激光粒度粒形仪（Mastersizer 3000, 英国）分析土壤黏粒含量（Clay）。过100目筛的土样用哈纳水质分析仪（HANNA HI98194, 意大利）测定pH值和电导率（Electrical conductivity, EC）。土样（过100目筛）经稀HCl酸化后，用元素分析仪测定有机碳含量（TOC）。使用元素分析仪测定土样（过100目筛）TC和TN含量。土样（过100目筛）经微波消解仪消解后采用等离子体原子发射光谱仪测定TP含量。

1.5   气候数据获取

采用WorldClim v2中的Average temperature（℃）、Precipitation（mm）和Bioclimatic variables（分辨率均为30″）栅格数据集，将其导入到ArcGIS 10.6软件中，使用Spatial Analyst Tools工具箱中的Extract Multi Values To Points工具，提取26个采样位点的月平均温度（Average temperature）、月降水量（Average precipitation），并计算年平均温度（MAT）和年降水量（MAP）等两个气候参数。

1.6   统计分析

实验数据使用R 3.6.2软件进行分析。用Shapiro-Wilk法检测数据的正态性和方差齐性，对不符合正态分布的数据进行对数转换或者平方根转换。采用Mantel分析植物功能性状与细菌多样性的关系，采用相关性分析探讨环境因子与植物多样性、细菌多样性的关系。利用多元回归分析寻找植物多样性和细菌多样性的关键影响因子，并通过变差分解得到各因子的解释率。使用基于距离的冗余分析（Distance-based redundancy analysis, db-RDA）对细菌群落的组成情况进行可视化分析，并通过与排序轴的多元回归置换检验得出影响细菌群落组成的关键环境因子和植物因子。

2.   结果与分析

2.1   植物与土壤细菌的多样性

本研究对26个样地的调查结果表明，植物丰富度为1 ～ 6，平均值为3.04；植物香农指数在0 ~ 1.80，平均为1.06（表1）。土壤细菌的Chao1指数在3654 ～ 6509，平均值为4859；土壤细菌的香农指数最小值为8.72，最大值为11.30。

表  1  植物与土壤细菌多样性

Table  1.  Plant and soil bacterial diversity

指标 Index		最小值 Minimum	最大值 Maximum	平均值 Mean
	植物丰富度 Plant Richness	1	6	3.04
	植物香农指数 Plant Shannon	0	1.80	1.06
	细菌Chao1指数 Bacteria Chao1	3654	6509	4859
	细菌香农指数 Bacteria Shannon	8.72	11.30	10.41

下载: 导出CSV 

| 显示表格

2.2   环境因子与植物多样性的关系

相关性分析结果表明（图2），植物丰富度仅与土壤电导率显著相关，而与其他环境因子没有明显的相关性；植物香农指数与所有环境因子的相关性均不显著。回归分析结果显示（表2），海拔、年平均温度和土壤电导率仅能解释植物丰富度26.71%的变异，但能解释植物香农指数31.92%的变异，其中，土壤电导率是解释植物丰富度和香农指数变异的首要因子，解释量分别为18.41%和16.79%。

图  2  环境因子与植物多样性、土壤细菌多样性的相关性

Figure  2.  Correlations between environmental factors and plant and soil bacterial diversity

下载: 全尺寸图片 幻灯片

表  2  土壤细菌多样性、植物多样性与环境因子的回归分析解释率

Table  2.  Multiple linear regression analysis of plant and soil bacterial diversity and environmental factors

指标 Index	细菌Chao1指数 Bacterial Chao1 / %	细菌香农指数 Bacterial Shannon / %	植物丰富度指数 Plant Richness / %	植物香农指数 Plant Shannon / %
纬度 Latitude	8.50	6.00	−	−
海拔 Altitude	8.57	21.69	8.30	6.82
年降水 MAP	16.30	−	−	−
pH 值	22.57	16.52	−	−
年均温 MAT	−	−	−	8.31
电导率 EC	−	−	18.41	16.79
总解释率 R2	56.25	44.21	26.71	31.92

下载: 导出CSV 

| 显示表格

2.3   环境因子、植物与土壤细菌多样性的关系

相关性分析结果表明（图2），土壤细菌Chao1指数与土壤电导率呈负相关，土壤细菌Chao1、香农指数均与海拔及年降水量呈正相关，而与土壤pH呈负相关。同时，土壤细菌香农指数还与年均温、土壤总有机碳含量显著相关。进一步进行回归分析，结果显示，纬度、海拔、年降水量和土壤pH值能够解释土壤细菌Chao1指数56.25%的变异；同时，纬度、海拔和土壤pH能够解释土壤细菌香农指数44.21%的变异（表2）。

Mantel分析结果表明，土壤细菌Chao1多样性指数与群落水平的植物叶片总氮含量（CWM.TN）呈显著正相关关系（图3）。回归分析表明，植物叶片干物质含量（CWM.LDMC）与植物丰富度是解释土壤细菌Chao1指数变化的主要因子，分别贡献了18.21%和12.14%的解释度，且细菌Chao1指数随CWM.LDMC与植物丰富度的增加而增加。植物叶片总碳含量（CWM.TC）是影响土壤细菌香农指数变异的关键因子，贡献了14.01%的解释度（表3），且CWM.TC与土壤细菌香农指数呈正相关。

图  3  植物多样性、功能性状和土壤细菌多样性的Mantel相关分析

Figure  3.  Mantel correlation analysis of plant diversity, functional traits, and soil bacterial diversity

下载: 全尺寸图片 幻灯片

表  3  植物多样性、功能性状和土壤细菌多样性的回归分析解释率

Table  3.  Regression analysis of plant diversity, functional traits, and soil bacterial diversity

指标 Index	细菌Chao1指数 Bacterial Chao1 / %	细菌香农指数 Bacterial Shannon / %
群落加权平均叶干物质含量 CWM.LDMC	18.21	−
植物丰富度指数Richness	12.14	−
群落加权平均叶总碳含量CWM.TC	−	14.01
总解释率 R2	30.35	14.01

下载: 导出CSV 

| 显示表格

冗余分析结果显示（图4：A），群落水平的植物功能性状（叶片总氮含量CWM.TN和总碳含量CWM.TC）和土壤理化性质（总碳和pH值）对细菌群落组成影响最大，共解释了12.48%的变异，其中前两轴分别解释了5.58%和2.90%。土壤总碳和pH值与第1轴呈正相关，而叶片总氮和总碳含量与第1轴呈负相关。韦恩图分析结果表明（图4：B），叶片总氮和总碳、土壤总碳和pH值分别解释了细菌群落3.15%、1.73%、2.11%和7.12%的差异，而它们的共同解释率很低，仅为0.08%。此外，与植物因子相比，环境因子对细菌群落组成的影响更大（图4：B）。

图  4  植物、环境因子与土壤细菌群落组成的关系

A：冗余分析；B：韦恩分析。

Figure  4.  Relationships between plants, environmental factors, and soil bacterial community composition

A: Redundancy analysis; B: Venn analysis.

下载: 全尺寸图片 幻灯片

3.   讨论

植物多样性的变化受多种环境因素的影响，包括气候、海拔和土壤理化性质等。本研究结果表明，西藏高原湿地土壤电导率与植物丰富度呈显著负相关，且是植物多样性变异的首要驱动因子。土壤电导率直接反映土壤的盐度，过高的土壤电导率会改变植物根区的渗透压，影响植物的代谢路径，进而影响其生长速率^[33]。崔乔等^[34]对盐沼湿地的研究表明，植物群落组成的模式取决于土壤电导率。气候因子（年平均温度和年降水量）也是影响植物多样性的关键因素，一方面可直接制约植物的生长、分布和繁衍，另一方面决定了土壤的水热条件，从而间接影响植物的多样性^{[35, 36]}。本研究结果表明，年平均温度是影响植物多样性的重要因子之一，而年降水量与植物多样性并没有显著的相关性，这可能是因为青藏高原高寒湿地的植物多样性更易受温度的影响^[37]。Hong等^[38]在河岸带湿地的研究也发现温度而非降水是影响植物物种多样性的关键因子。海拔的变化通常伴随着气候条件的改变（如温度随海拔上升而下降），进而影响植物的多样性。在西藏高原湿地中，我们发现海拔也是影响植物多样性的一个重要因子，随海拔的升高植物多样性显著降低。

细菌多样性的变化往往与环境因子密切相关。本研究发现，环境因子能够解释土壤细菌Chao1指数56.25%的变异和香农指数44.21%的变异，是影响西藏高原湿地土壤细菌多样性变化的关键因素。此外，土壤pH是湿地细菌多样性变化的关键驱动因子，细菌多样性随土壤pH值的增加而显著降低，这与Kang等^[39]的研究结果一致。土壤pH会通过影响土壤基质组成、化学性质和养分利用效率来改变细菌群落的活性和多样性^[40]。同时，海拔和纬度也是细菌多样性变化的关键驱动因子，细菌的Chao1指数和香农指数均随海拔的增加而显著增加。Wang等^[41]对不同海拔梯度细菌群落的生物多样性模式进行了研究，结果发现细菌多样性与海拔呈显著正相关。虽然纬度也是影响细菌多样性变化的关键因子，但本研究发现二者之间不存在显著相关。这可能是因为纬度的变化主要是通过影响气候条件和土壤理化性质，进而改变细菌的多样性。细菌Chao1指数还与年均降水量呈显著正相关，与杨阳等^[42]的研究结果一致。

西藏高原湿地植物的功能性状和物种丰富度也是驱动土壤细菌多样性变化的关键因子。很多研究发现土壤细菌多样性随植物多样性的增加而增加^{[43, 44]}，但也有研究指出两者并无显著相关性^[45]。植物多样性很可能是通过增加凋落物和根系分泌物等的营养资源，从而促进了土壤细菌多样性的增加^{[46 - 48]}。此外，植物多样性的变化也会引起土壤理化特征的改变（如电导率、有机碳含量等），进而影响土壤细菌多样性。Cao等^[49]研究发现植物可以通过影响土壤有机碳、总氮等理化因子间接影响土壤细菌多样性。本研究结果表明，与环境因子相比，植物因子对土壤细菌多样性变化的贡献较小，这也说明西藏高原湿地土壤的细菌多样性主要由环境因子（气候、土壤理化性质）决定。