生态学报  2018, Vol. 38 Issue (11): 3816-3824

文章信息

高静, MUHANMMAD Said, 岳琳艳, 何永涛, 斯确多吉, 张宪州, 孔维栋.
GAO Jing, MUHANMMAD Said, YUE Linyan, HE Yongtao, TSECHOE Dorji, ZHANG Xianzhou, KONG Weidong.
藏北高原草甸土壤固碳微生物群落特征随海拔和季节的变化
Changes in CO2-fixing microbial community characteristics with elevation and season in alpine meadow soils on the northern Tibetan Plateau
生态学报. 2018, 38(11): 3816-3824
Acta Ecologica Sinica. 2018, 38(11): 3816-3824
http://dx.doi.org/10.5846/stxb201705160904

文章历史

收稿日期: 2017-05-16
网络出版日期: 2018-03-02
藏北高原草甸土壤固碳微生物群落特征随海拔和季节的变化
高静1,2 , MUHANMMAD Said1 , 岳琳艳1,2 , 何永涛3 , 斯确多吉1 , 张宪州3 , 孔维栋1,2     
1. 中国科学院青藏高原研究所, 高寒生态学与生物多样性重点实验室, 北京 100101;
2. 中国科学院大学资源与环境学院, 北京 100049;
3. 中国科学院地理与资源科学研究所, 生态系统网络观测与模拟重点实验室, 北京 100101
摘要: 研究土壤固碳微生物丰度、群落结构、多样性差异及其影响因子对了解青藏高原土壤碳循环和固碳潜力具有重要意义。采用定量PCR(qPCR)、末端限制性片段分析(T-RFLP)、克隆文库和测序方法,研究了青藏高原草甸土壤固碳微生物丰度与群落结构随海拔和季节的变化,主要结果如下:1)随海拔升高高寒草甸土壤固碳微生物丰度显著升高,但季节变化不明显,不同类别微生物固碳基因cbbL丰度依次为:Form IC > Form IAB > Form ID,其中Form IC类固碳微生物可达108拷贝数/g土壤,cbbL基因丰度与海拔、土壤含水量和铵态氮含量(NH4+-N)呈正相关关系,与土壤温度和pH值负相关;2)固碳微生物多样性和丰富度随海拔升高而升高,在4800m达到最大,且二者受季节影响较小,其群落结构随海拔升高而逐渐变化,主要受土壤pH值、海拔和土壤水分影响;3)Form IC类固碳微生物主要包括放线菌门和和变形菌门,其中α变形菌门是高寒草甸土壤优势固碳微生物类群。本研究有助于理解土壤微生物群落功能及其在土壤碳循环过程中的作用,为更准确评估高寒草甸土壤碳循环过程提供了科学依据。
关键词: 高寒草甸     固碳微生物     RubisCO酶     海拔梯度     青藏高原    
Changes in CO2-fixing microbial community characteristics with elevation and season in alpine meadow soils on the northern Tibetan Plateau
GAO Jing 1,2, MUHANMMAD Said 1, YUE Linyan 1,2, HE Yongtao 3, TSECHOE Dorji 1, ZHANG Xianzhou 3, KONG Weidong 1,2     
1. Key Laboratory of Alpine Ecology and Biodiversity, Institute of Tibetan Plateau Research, Chinese Academy of Sciences, Beijing 100101, China;
2. College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100049, China;
3. Key Laboratory of Ecosystem Network Observation and Modeling, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
Abstract: Soil microbial autotrophs play an important role in CO2 fixation in terrestrial ecosystems, particularly in vegetation-constrained ecosystems with environmental stresses, such as the Tibetan Plateau, which is characterized by low temperature, drought, and high UV. However, soil microbial autotrophic communities and their driving factors remain less appreciated in these terrestrial ecosystems. To understand carbon sequestration by soil microbial autotrophs and the carbon cycle in alpine meadow soils on the Tibetan Plateau, autotrophic microbial abundance, community structure, diversity, and their driving environmental factors were explored along an elevation gradient from 4400m to 5200m on the Tibetan Plateau. Additionally, the seasonal change in soil microbial autotrophs was explored at each elevation. The autotrophic microbial communities were characterized by quantitative PCR, T-RFLP, and cloning/sequencing methods targeting four types of cbbL gene, encoding the large subunit for the CO2-fixing protein ribulose 1, 5-bisphosphate carboxylase/oxygenase (RubisCO), Form IA/B, IC, and ID. High cbbL gene abundance was observed, which significantly increased with increasing elevations, with minimum seasonal change. Among the four types of microbial autotrophs, Form IC gene abundance was the highest at the level of 108 copies/g soil, followed by Form IA/B and Form ID. Form IC cbbL gene abundance was positively correlated with elevation, soil water content (SWC), and NH4+ content, and negatively correlated with soil temperature and pH (P < 0.01). Form IC Shannon diversity and richness significantly increased with increasing elevation, peaking at 4800m, but showed little effect of season. Form IC community structure exhibited a gradual shift along the elevation gradient, and was mainly driven by SWC, soil pH, and elevation. Form IC cbbL genes were identified to be Actinobacteria and Proteobacteria (including α-proteobacteria and β-proteobacteria) using clone library and sequencing, and the dominant phyla were α-proteobacteria in alpine meadow soils, including Bradyrhizobium, Rhizobium, and Mesorhizobium. Our results demonstrated high abundance of microbial autotrophs and high CO2 fixation potential in alpine grassland soils, and provide a novel model to identify dominant drivers of soil microbial communities and their ecological functions. These findings have important implications for understanding the functional role of microorganisms in the soil carbon cycle and could be helpful to more accurately assess carbon storage in alpine grassland soils on the Tibetan Plateau.
Key words: alpine meadow     cbbL     RubisCO     elevation     the Tibetan Plateau    

卡尔文循环是自然界中最主要的生物固碳途径[1]。核酮糖-1, 5-二磷酸羧化酶/加氧酶(RubisCO)是卡尔文循环的限速酶和关键酶[2]。目前自然界存在4种类型RubisCO (Ⅰ, Ⅱ, Ⅲ, Ⅳ), 其中RubisCOⅠ是广泛存在于光合和化能自养生物体中[3]。RubisCOⅠ又可细分为四类, 分别是Form IA, IB, IC, ID。RubisCO酶大亚基由cbbL基因编码, 95%以上的固碳生物含有该基因[4]cbbL基因可作为固碳微生物标记物, 已在海洋[5]、湖泊[6]和土壤[7]固碳微生物群落结构和生物多样性研究中广泛应用。海洋和湖泊中固碳微生物群落及其固碳能力已经得到很多研究, 但在土壤固碳微生物群落特征研究很少[7]

青藏高原因海拔高、年均温低、干旱和强紫外等等恶劣环境, 导致植物生长受到极大限制, 是典型生态脆弱区[8]和气候变化敏感区[9]。高寒草甸是青藏高原典型地带性植被生态系统[10], 该生态系统土壤中碳储量巨大[11]。近年来高寒草甸土壤微生物研究主要集中在土壤微生物量、土壤微生物活性[12]、土壤微生物群落结构变化[13]和土壤微生物多样性[14]等方面, 但土壤固碳微生物特征及其对高寒草甸土壤碳储量的贡献尚未得到深入研究。与低海拔地区相比, 青藏高原植被生长受到极大限制, 因此, 开展土壤微生物固碳研究对我们理解青藏高原土壤固碳微生物功能及土壤碳循环过程具有重要意义。

海拔梯度变化综合了多种环境因子的梯度效应[15]。随着海拔梯度的上升, 土壤和气候等环境因子均发生改变, 从而引起土壤微生物群落多样性的变化[16]。因此, 海拔梯度变化是研究土壤固碳微生物群落特征及其驱动因子的理想自然实验室。本文采用分子生物学方法(分别为定量PCR、T-RFLP、克隆文库和测序等方法), 研究了青藏高原不同海拔草甸土壤固碳微生物群落结构和多样性变化规律, 旨在探讨不同海拔梯度土壤固碳微生物群落特征及其环境影响因子。该研究可为深入理解高寒草甸土壤固碳微生物群落特征及土壤碳循环过程提供理论依据。

1 研究区域与方法 1.1 研究地区概况

本研究区域位于拉萨市北部当雄县辖区, 距离拉萨市区约170km, 坐标为91°05′E, 30°51′N, 平均海拔约4300—4700m, 背靠念青唐古拉山脉[17]。该地区属于高原性季风气候, 水热同期, 夏季短暂多雨;冬季寒冷干燥, 昼夜温差大。年均温1.5℃(1960—2000)[18], 年降水量476.8mm, 其中约80%的降水都集中在3月份至8月份, 降水季节性十分明显[19]。土壤类型为高寒草甸土, 生态系统为典型高山草甸, 植被主要包括建群种矮嵩草, 夹杂委陵菜、垫状点地梅等, 随海拔升高植被类型变化不大, 仅地上部生物量增加[20]

1.2 样品采集与处理

在念青唐古拉山南坡采样, 在2014年4个季节采样, 时间分别为2014年3月26日、6月21日、8月26日和10月23日, 样点分布在4400、4650、4800m和5200m的4个海拔, 每个海拔在取样点周围9m2内取3个重复, 每个重复土壤由直径2.5cm的土钻取3钻混合而成, 取样深度为0—10cm。土样采集后装塑封袋, 过2.0mm孔径筛子后一部分用无菌袋封装并迅速冷冻, 用于微生物群落实验, 于-80℃冷冻保存;另一部分风干后用于土壤理化因子测定。

1.3 土壤理化因子测定

测定土壤的碳氮相关理化因子, 主要为土壤总有机碳(TOC)和铵态氮(NH4+-N)和硝态氮(NO3--N)。TOC采用总有机碳分析仪测定(TOC-L, SHIMADZU, Japan), 铵态氮和硝态氮采用流动分析仪测定(AQ2, SEAL Analytical Inc., UK)。土壤含水量采用烘干法(24h)测定, 土壤pH值(土水比1: 2.5)采用酸度计(Sartorius PB-10, 德国)测定, 每个海拔土壤温度由HOBO自动气象站测定(Onset Inc., Bourne, MA, USA)。

1.4 土壤DNA提取

所有DNA样品均采用Powersoil® DNA Isolation Kit(MOBIO, USA)试剂盒提取, 土壤称取冷冻土样0.25g。操作流程按照试剂盒说明书进行。DNA提取完毕后取2μL DNA溶液检测DNA浓度和纯度(Nanodrop® ND-1000, USA), 并用琼脂糖凝胶电泳检测DNA完整性。

1.5 基因丰度测定

采用定量PCR(简称qPCR)方法检测样品中和固碳微生物cbbL基因丰度, PCR扩增引物均为特异性引物, 其序列见表 1。基因丰度检测使用罗氏定量PCR仪(Roche LightCycler480 Ⅱ, Roche, USA)进行, 标线采用含目的片段的已知数量质粒10倍系列稀释液[6]。qPCR扩增酶使用SYBR green kit(TaKaRa, 大连), 反应体系为10μL, 包括3μL缓冲液, 5μL的酶, 引物各0.5μL, 1μL DNA模板。qPCR反应条件为预变性95/94℃ 2 min, 94℃变性30 s, 退火温度分别是:Form IA/B基因61℃, Form IC基因63℃, Form ID基因53℃ 35个循环。在72℃收集荧光信号。

表 1 本研究采用的引物 Table 1 Primer sets used in the study
引物名称Primer 引物序列Sequence(5′-3′)
Form IC-F GAA CAT CAA YTC KCA GCC CTT[21]
Form IC-R TGG TGC ATC TGV CCG GCR TG[21]
Form ID-F GAT GAT GAR AAY ATT AAC TC[22]
Form ID-R ATT TGD CCA CAG TGD ATA CCA[22]
Form IA/B-F TCI GCI TGR AAC TAY GGT CG[23]
Form IA/B-R GGC ATR TGC CAI ACR TGR AT[23]
1.6 固碳微生物群落结构分析

草甸土壤中Form IC类固碳微生物丰度最高, 代表主要土壤固碳微生物类群, 远高于Form IAB和ID类固碳微生物。因此, 本研究采用末端限制性片段长度多态性(Terminal Restriction Fragment Length Polymorphism, 简称T-RFLP)技术深入分析了Form IC类固碳微生物群落结构和多样性。利用5′端标记有荧光的引物对功能基因进行扩增, 而后用Axygen公司生产的胶回收试剂盒对PCR产后进行切胶纯化, 将纯化产物进行限制性酶切, 本研究中采用MspI酶切Form IC目的基因片段, 酶切产物片段大小和相对量采用ABI 3730xl DNA测序仪(Applied Biosystems, CA, USA)测定, 不同末端限制片段(T-RFs)相对量用来计算固碳微生物群落结构和多样性指数。

1.7 固碳微生物克隆文库、测序与DNA序列系统发育分析

我们对4800m海拔处3月和6月土样各建立了Form IC类固碳微生物克隆文库, 并测序和进行系统发育分析。IC固碳基因PCR扩增采用与qPCR相同的引物, PCR反应体系与T-RFLP相同。将每个克隆文库的PCR样品按照的胶回收试剂盒说明书进行目的片段的回收纯化, 纯化产物用Promage公司生产的PEGM-T(Promega, USA)连接载体试剂盒进行连接反应并转化到感受态细胞中, 感受态细胞采用康为世纪公司生产的DH-5α菌株;通过LB平板蓝白斑筛选随机挑取白斑(40—45个);对随机挑选的白斑进行阳性克隆鉴定和测序, 鉴定结果为阳性、片段长度符合目标片段要求的克隆子进行Sanger测序(ABI 3730xl DNA测序仪, Applied Biosystems, CA, USA);采用MEGA 6.0对获得的目的基因DNA序列进行系统发育分析。

1.8 数据分析

克隆文库测序所得结果使用Mothur v.1.33.3选取相似性大于97%的DNA序列进行OTU划分, 使用邻接法(neighbor-joining)在MEGA 6.0软件中构建系统发育树。群落结构分析主要利用Canoco 5软件进行, 群落多样性利用R 3.1.3中的“vegan”软件包进行计算。相关性和差异显著性分析采用SPSS. 23软件, 作图采用Sigmaplot 10.0和Excel软件。

2 结果与分析 2.1 cbbL基因丰度随海拔和季节变化

土壤中各月份Form IAB基因丰度在106 拷贝数/g到108拷贝数/g之间(图 1)。随海拔升高, Form IAB的基因丰度逐渐增加, 且由4400m至4800m Form增加幅度较大, 4800m以上基本持平或稍下降。3月份Form IAB的基因丰度随海拔变化幅度最小。对海拔和季节进行双因素分析可知, 季节及其与海拔的交互作用对Form IAB基因丰度无显著影响(P>0.631)。

图 1 Form IAB基因丰度随海拔的变化图 Fig. 1 Changes in the abundance of Form IA/B with altitude

Form IC基因丰度各月份在108拷贝数/g到109拷贝数/g之间(图 2)。随海拔升高, 各月份Form IC基因丰度显著增加, 自4400m至4800m增加幅度较快, 4800m以上基本持平或有所下降。与其他月份相比, 3月份Form IC基因丰度随海拔升高变化最小。对海拔和季节进行双因素分析可知, 海拔升高和季节均显著影响Form IC基因丰度(P < 0.05), 且二者交互作用也显著影响Form IC基因丰度(P < 0.01)。

图 2 Form IC基因丰度随海拔的变化图 Fig. 2 Changes in the abundance of Form IC with altitude

Form ID基因丰度各月份在105拷贝数/g到106拷贝数/g之间(图 3)。随海拔升高, 各月份Form ID基因丰度无明显变化, 且10月份Form ID基因丰度随海拔升高变化最小。

图 3 Form ID基因丰度随海拔的变化图 Fig. 3 Changes in the abundance of Form ID with altitude

综上所述, 随海拔升高不同类型固碳微生物丰度变化趋势不同, Form IC和Form IAB基因丰度随海拔升高而逐渐增加, Form ID随海拔变化趋势不明显。其中Form IC类固碳微生物丰度在各个海拔均远高于其他两类固碳微生物。对比不同季节间基因丰度的变化情况, 在生长季(6月和8月)cbbL基因丰度变化较大, 而冬季(3月)变化幅度较小。

Pearson相关分析表明, 土壤Form IC基因丰度与海拔呈正相关关系, 即海拔越高, Form IC基因丰度越大(表 2), 说明高海拔土壤条件有利于IC类固碳微生物生长。Form IAB和Form IC类固碳微生物丰度与土壤温度和pH值呈负相关关系, 与土壤含水量和NH4+-N呈正相关关系。Form IC基因丰度与土壤TOC含量呈正相关关系。Form ID类固碳微生物丰度与海拔呈负相关关系, 即随海拔升高而降低, 与土壤温度呈正相关关系。

表 2 cbbL基因丰度与环境因子的相关性 Table 2 Pearson correlations between cbbL gene abundance and environmental factors
环境因子
Environmental factors
Form IC Form ID Form IAB
土壤含水量(SWC)
Soil water content
0.484** -0.064 0.341**
NH4+-N 0.316* -0.053 0.336**
NO3--N -0.027 0.050 -0.129
pH -0.529** 0.157 -0.370**
海拔
Altitude
0.327* -0.275* 0.133
土壤温度
Soil temperature
-0.371** 0.308* -0.289*
总有机碳
TOCTotal organic carbon
0.312* -0.175 0.120
*和**分别表示显著性达到0.05和0.01水平, 单尾检验
2.2 cbbL基因多样性指数随海拔和季节变化

随海拔升高各月份Form IC cbbL基因香农多样性指数均显著升高, 4800m出达到最高, 之后稍下降(图 4)。与香农多样性指数变化趋势相似, 随海拔显著升高各月份Form IC基因群落丰富度显著升高, 在4800m达到最大, 之后有所下降(图 5)。与其他月份相比, 3月份香农多样性指数和丰富度变化均最小, 可能是由于3月份土壤还在冻融期, 土壤温度与土壤水含量相对较低, 故对海拔变化不敏感。

图 4 Form IC基因香农多样性指数随海拔的变化 Fig. 4 Changes in the Shannon diversity of Form IC with altitude

图 5 Form IC基因群落丰富度随海拔的变化 Fig. 5 Changes in the richness of Form IC with altitude
2.3 草甸土壤固碳微生物群落结构随海拔和季节变化

Form IC类土壤固碳微生物群落结构随海拔升高而发生明显变化, 但随季节变化不明显(图 6)。根据海拔变化, 固碳微生物群落结构主要分为两个类群, 海拔4400m和4650m聚合为一类, 4800m和5200m聚合一类。RDA1和RDA2两个分量可解释30.36%群落结构变异, RDA1可单独解释变异的23.65%, 说明海拔升高导致的土壤水分、有机质含量、pH和土壤温度等土壤理化因子是驱动固碳微生物群落结构变化的主要因子。RDA2仅解释固碳微生物群落结构变异的6.71%, 其中季节是驱动RDA2的主要因子, 说明群落结构在一定程度存在明显季节差异。

图 6 Form IC基因群落结构RDA分析图 Fig. 6 Redundancy analysis of T-RFLP of Form IC genes E:海拔, altitude;SWC:土壤含水量, soil water content;TOC:总有机碳, total organic carbon;T_5:5cm土壤温度, 5cm soil temperature;Season:季节, season; NH4+:铵态氮, ammonium nitrogen;NO3-:硝态氮, nitrate nitrogen
2.4 固碳微生物群落组成及其系统发育分析

本研究对4800m海拔处3月和6月土壤进行克隆文库构建和cbbL基因片段序列分析, 共得到52个OTUs。本研究土壤中Form IC类固碳微生物属于3个门(图 7), 分别是放线菌(Actinobacteria)、α变形菌(α-proteobacteria)和β变形菌(β-proteobacteria)。放线菌主要PseudomocardiaNocardia等微生物, 与之前发现的一些青藏高原草地土壤克隆序列比较相似。α变形菌主要包括一些典型根瘤菌, 如Bradyrhizobium, MesorhizobiumRhizobium等, β变形菌主要包括一些化能自养型微生物, 如RubrivivaxVariovorvax等固碳微生物。

图 7 Form IC基因系统发育树 Fig. 7 Neighbor-joining phylogenetic tree of representative form IC cbbL gene sequences

α变形菌占总克隆文库的比例最高(48.1%), 放线菌占38.4%, β变形菌比例最低为13.5%。对比不同月份3种类群所占比例, 3月份α变形菌与放线菌所占比例相当分别占45%;6月份α变形菌比例从45%升高至52%, 而放线菌比例从45%降低至30%, 说明主要包括根瘤菌的α变形菌在生长季比例稍高于非生长季, 固碳微生物群落组成随季节变化。

3 讨论与结论

海拔梯度变化综合了多种环境因子的梯度效应[15]。在高寒草甸生态系统中, 随海拔梯度变化, 植物指数对气候的敏感性[24]、微生物的群落组成及多样性都会发生改变[25]。本研究发现, 随海拔升高, 草甸土壤温度与pH值逐渐降低, 土壤含水率与铵态氮(NH4+-N)含量逐渐升高, 这些结果与我们之前的发现一致[26]。草甸土壤中类型Form I固碳微生物中, Form IC基因丰度最高, 可达108拷贝数/L以上, 且随海拔升高而增加。cbbL基因丰度与土壤温度和pH值呈负相关关系, 与土壤含水量、NH4+-N和TOC呈正相关关系, 表明高海拔土壤中丰富的养分(包括TOC和NH4+-N)刺激了固碳微生物生长, 同时这些固碳微生物可忍受或适应较低温度。因此, 适应低温可能是青藏高原草甸土壤固碳微生物关键生存策略之一。丛静等[27]发现神农架自然保护区土壤微生物生物量随海拔升高而显著增加, 随海拔变化的土壤温度、土壤湿度、pH值和有机碳是重要的影响因子, 与本文研究结果相似。对比不同季节的cbbL基因丰度可发现, 生长季(6月和8月)cbbL基因丰度随海拔变化较大, 而非生长季(3月和10月)随海拔变化幅度较小。3月份各类固碳微生物丰度在各海拔梯度丰度基本不变, 说明3月份土壤处仍处于冷冻或冻融交替阶段, 土壤微生物还没有生长, 故各海拔梯度固碳微生物丰度基本相似。这些研究结果与王琴等[28]的发现一致, 即土壤微生物类群的海拔梯度效应存在季节变化特征。

海拔及其导致的土壤理化因子变化和季节均影响固碳微生物结构多样性。土壤pH值、海拔和水分是驱动群落结构变化的关键环境因子, 季节在一定程度上也会影响固碳微生物群落结构, 但其影响程度稍弱于土壤理化因子。不同海拔和季节均使土壤理化和养分发生改变, 从而影响固碳微生物生长和代谢活动, 最终导致固碳微生物种群结构的变化。肖可青等研究发现在稻田土壤中, 土壤有机碳和pH值对固碳微生物群落结构具有显著影响[29]。本文研究表明, 随海拔升高草甸土壤固碳微生物多样性和丰富度逐渐升高, 在4800m达到最高, 说明海拔梯度会影响固碳微生物丰度和多样性变化, 这与我们之前的研究结果一致[26]。宋贤冲等人[16]认为土壤含水量可能是造成不同海拔土壤微生物群落功能多样性差异的主要原因。土壤克隆文库和测序结果表明, 在生长季(3月)和非生长季(6月)草甸土壤Form IC类固碳微生物主要包括3个菌门, 分别是放线菌、α和β变形菌。其中α变形菌为优势菌门, 主要包括根瘤菌, 这与我们之前的研究发现相似[26]。刘琼等人[30]也发现稻田土壤中固碳微生物主要为变形菌和放线菌。有研究表明, β变形菌可以耐受各种环境压力, 尤其对于酸性土壤中低pH值具有良好的耐受性[31]

综上所述, 在本研究海拔梯度内, 草甸土壤固碳微生物丰度和群落结构主要受海拔变化影响, 而受季节变化影响较小, 其中, 随海拔梯度变化的土壤温度、土壤含水量、pH值和有机碳是重要的影响因子。固碳微生物多样性主要受土壤pH值和土壤含水量影响。该研究揭示了高寒草甸土壤固碳微生物丰度、群落结构特征及其环境驱动因子, 为深入理解高寒草地土壤微生物功能和土壤碳循环过程及其动态变化提供理论依据。

参考文献
[1] Tabita F R. Molecular and cellular regulation of autotrophic carbon dioxide fixation in microorganisms. Microbiological Reviews, 1988, 52(2): 155–189.
[2] 陈为钧, 赵贵文, 顾月华. RubisCO的研究进展. 生物化学与生物物理进展, 1999, 26(5): 433–436.
[3] Tabita F R, Satagopan S, Hanson T E, Kreel N E, Scott S S. Distinct form Ⅰ, Ⅱ, Ⅲ, and Ⅳ RubisCO proteins from the three kingdoms of life provide clues about RubisCO evolution and structure/function relationships. Journal of Experimental Botany, 2008, 59(7): 1515–1524.
[4] Uchino Y, Yokota A. "Green-like" and "red-like" RubisCO cbbL genes in Rhodobacter azotoformans. Molecular Biology and Evolution, 2003, 20(5): 821–830. DOI:10.1093/molbev/msg100
[5] John D E, Wang Z A, Liu X W, Byrne R H, Corredor J E, López J M, Cabrera A, Bronk D A, Tabita F R, Paul J H. Phytoplankton carbon fixation gene (RuBisCO) transcripts and air-sea CO2 flux in the Mississippi River plume. The ISME Journal, 2007, 1(6): 517–531.
[6] Kong W D, Ream D C, Priscu J C, Morgan-Kiss R M. Diversity and expression of RubisCO genes in a perennially ice-covered antarctic lake during the polar night transition. Applied and Environmental Microbiology, 2012, 78(12): 4358–4366. DOI:10.1128/AEM.00029-12
[7] 陈晓娟, 吴小红, 简燕, 袁红朝, 周萍, 葛体达, 童成立, 邹冬生, 吴金水. 农田土壤自养微生物碳同化潜力及其功能基因数量、关键酶活性分析. 环境科学, 2014, 35(3): 1144–1150.
[8] 孙鸿烈, 郑度, 姚檀栋, 张镱锂. 青藏高原国家生态安全屏障保护与建设. 地理学报, 2012, 67(1): 3–12. DOI:10.11821/xb201201001
[9] 姚檀栋, 朱立平. 青藏高原环境变化对全球变化的响应及其适应对策. 地球科学进展, 2006, 21(5): 459–464.
[10] Shen Z X, Fu G, Yu C Q, Sun W, Zhang X Z. Relationship between the growing season maximum enhanced vegetation index and climatic factors on the Tibetan Plateau. Remote Sensing, 2014, 6(8): 6765–6789. DOI:10.3390/rs6086765
[11] 王根绪, 程国栋, 沈永平. 青藏高原草地土壤有机碳库及其全球意义. 冰川冻土, 2002, 24(6): 693–700.
[12] 任佐华, 张于光, 李迪强, 肖启明, 蔡重阳. 三江源地区高寒草原土壤微生物活性和微生物量. 生态学报, 2011, 31(11): 3232–3238.
[13] 斯贵才, 王光鹏, 雷天柱, 张更新, 夏燕青, 陈芳荣. 青藏高原东北缘土壤微生物群落结构变化. 干旱区研究, 2015, 32(5): 849–855.
[14] 周宁一. 青藏高原微生物多样性研究. 微生物学通报, 2014, 41(11): 2378–2378.
[15] Djukic I, Zehetner F, Mentler A, Gerzabek M H. Microbial community composition and activity in different alpine vegetation zones. Soil Biology and Biochemistry, 2010, 42(2): 155–161. DOI:10.1016/j.soilbio.2009.10.006
[16] 宋 贤冲, 郭 丽梅, 田 红灯, 邓 小军, 赵 连生, 曹 继钊. 猫儿山不同海拔植被带土壤微生物群落功能多样性. 生态学报, 2017, 37(16): 1–8.
[17] 付刚, 沈振西. 藏北高原不同海拔高度高寒草甸蒸散与环境温湿度的关系. 中国草地学报, 2015, 37(3): 67–73.
[18] Ohtsuka T, Hirota M, Zhang X Z, Shimono A, Senga Y, Du M G, Yonemura S, Kawashima S, Tang Y H. Soil organic carbon pools in alpine to nival zones along an altitudinal gradient (4400-5300m) on the Tibetan Plateau. Polar Science, 2008, 2(4): 277–285. DOI:10.1016/j.polar.2008.08.003
[19] Fu G, Shen Z X, Zhang X Z, Shi P L, He Y T, Zhang Y J, Sun W, Wu J S, Zhou Y T, Pan X. Calibration of MODIS-based gross primary production over an alpine meadow on the Tibetan Plateau. Canadian Journal of Remote Sensing, 2012, 38(2): 157–168. DOI:10.5589/m12-023
[20] 周启龙. 西藏藏北地区不同草原类型区土壤理化性质特征研究. 安徽农业科学, 2016, 44(20): 124-125, 142-142. DOI:10.3969/j.issn.0517-6611.2016.20.040
[21] Alfreider A, Vogt C, Geiger-Kaiser M, Psenner R. Distribution and diversity of autotrophic bacteria in groundwater systems based on the analysis of RubisCO genotypes. Systematic and Applied Microbiology, 2009, 32(2): 140–150. DOI:10.1016/j.syapm.2008.11.005
[22] Paul J H, Alfreider A, Wawrik B. Micro-and macrodiversity in rbcL sequences in ambient phytoplankton populations from the southeastern Gulf of Mexico. Marine Ecology Progress Series, 2000, 198: 9–18. DOI:10.3354/meps198009
[23] Corredor J E, Wawrik B, Paul J H, Tran H, Kerkhof L, López J M, Dieppa A, Cárdenas O. Geochemical rate-RNA integration study:ribulose-1, 5-bisphosphate carboxylase/oxygenase gene transcription and photosynthetic capacity of planktonic photoautotrophs. Applied and Environmental Microbiology, 2004, 70(9): 5459–5468. DOI:10.1128/AEM.70.9.5459-5468.2004
[24] 沈振西, 孙维, 李少伟, 何永涛, 付刚, 张宪洲, 王江伟. 藏北高原不同海拔高度高寒草甸植被指数与环境温湿度的关系. 生态环境学报, 2015, 24(10): 1591–1598.
[25] 彭岳林, 蔡晓布. 丛枝菌根真菌群落沿高寒草原海拔梯度的变化特征. 生态学报, 2015, 35(22): 7475–7484.
[26] Guo G X, Kong W D, Liu J B, Zhao J X, Du H D, Zhang X Z, Xia P H. Diversity and distribution of autotrophic microbial community along environmental gradients in grassland soils on the Tibetan Plateau. Applied Microbiology and Biotechnology, 2015, 99(20): 8765–8776. DOI:10.1007/s00253-015-6723-x
[27] 丛静, 刘晓, 卢慧, 薛亚东, 刘学端, 尹华群, 李迪强, 张于光. 神农架自然保护区土壤微生物生物量碳、氮沿海拔梯度的变化及其影响因素. 生态学杂志, 2014, 33(12): 3381–3387.
[28] 王琴, 李菊, 孙辉. 海拔梯度上西南亚高山-高山土壤微生物生物量碳季节动态. 四川农业大学学报, 2013, 31(4): 386–392.
[29] Xiao K Q, Bao P, Bao Q L, Jia Y, Huang F Y, Su J Q, Zhu Y G. Quantitative analyses of ribulose-1, 5-bisphosphate carboxylase/oxygenase (RubisCO) large-subunit genes (cbbL) in typical paddy soils. Fems Microbiology Ecology, 2014, 87(1): 89–101. DOI:10.1111/1574-6941.12193
[30] 刘琼, 魏晓梦, 吴小红, 袁红朝, 王久荣, 李裕元, 葛体达, 吴金水. 稻田土壤固碳功能微生物群落结构和数量特征. 环境科学, 2017, 38(2): 760–768.
[31] Stopnisek N, Bodenhausen N, Frey B, Fierer N, Eberl L, Weisskopf L. Genus-wide acid tolerance accounts for the biogeographical distribution of soil Burkholderia populations. Environmental Microbiology, 2014, 16(6): 1503–1512. DOI:10.1111/emi.2014.16.issue-6