Loading [MathJax]/jax/output/SVG/jax.js
高级检索+

中国红豆杉TcNPR1基因的克隆与功能研究

谢莎, 许想平, 张蒙, 付春华, 余龙江

谢莎, 许想平, 张蒙, 付春华, 余龙江. 中国红豆杉TcNPR1基因的克隆与功能研究[J]. 植物科学学报, 2014, 32(4): 383-393. DOI: 10.3724/SP.J.1142.2014.40383
引用本文: 谢莎, 许想平, 张蒙, 付春华, 余龙江. 中国红豆杉TcNPR1基因的克隆与功能研究[J]. 植物科学学报, 2014, 32(4): 383-393. DOI: 10.3724/SP.J.1142.2014.40383
XIE Sha, XU Xiang-Ping, ZHANG Meng, FU Chun-Hua, YU Long-Jiang. Cloning and Function Analysis of the TcNPR1 Gene from Taxus chinesis[J]. Plant Science Journal, 2014, 32(4): 383-393. DOI: 10.3724/SP.J.1142.2014.40383
Citation: XIE Sha, XU Xiang-Ping, ZHANG Meng, FU Chun-Hua, YU Long-Jiang. Cloning and Function Analysis of the TcNPR1 Gene from Taxus chinesis[J]. Plant Science Journal, 2014, 32(4): 383-393. DOI: 10.3724/SP.J.1142.2014.40383
谢莎, 许想平, 张蒙, 付春华, 余龙江. 中国红豆杉TcNPR1基因的克隆与功能研究[J]. 植物科学学报, 2014, 32(4): 383-393. CSTR: 32231.14.SP.J.1142.2014.40383
引用本文: 谢莎, 许想平, 张蒙, 付春华, 余龙江. 中国红豆杉TcNPR1基因的克隆与功能研究[J]. 植物科学学报, 2014, 32(4): 383-393. CSTR: 32231.14.SP.J.1142.2014.40383
XIE Sha, XU Xiang-Ping, ZHANG Meng, FU Chun-Hua, YU Long-Jiang. Cloning and Function Analysis of the TcNPR1 Gene from Taxus chinesis[J]. Plant Science Journal, 2014, 32(4): 383-393. CSTR: 32231.14.SP.J.1142.2014.40383
Citation: XIE Sha, XU Xiang-Ping, ZHANG Meng, FU Chun-Hua, YU Long-Jiang. Cloning and Function Analysis of the TcNPR1 Gene from Taxus chinesis[J]. Plant Science Journal, 2014, 32(4): 383-393. CSTR: 32231.14.SP.J.1142.2014.40383

中国红豆杉TcNPR1基因的克隆与功能研究

基金项目: 

国家自然科学基金面上项目(31270342)。

详细信息
    作者简介:

    谢莎(1988-),女,硕士研究生,研究方向为植物次生代谢转录调控(E-mail:402765362@qq.com);许想平(1968-),女,中级技师,主要从事植物次生代谢研究(E-mail:506589862@qq.com)。

    通讯作者:

    付春华(1974-),女,副教授,主要从事植物次生代谢与调控研究(E-mail:fuch2003@126.com)。

  • 中图分类号: Q943.2

Cloning and Function Analysis of the TcNPR1 Gene from Taxus chinesis

  • 摘要: 水杨酸(SA)可诱导红豆杉细胞中紫杉醇的合成,病程相关基因非表达子基因(NPR1)是SA介导的植物系统获得性抗性(SAR)信号途径中的关键信号分子。本研究从中国红豆杉(Taxus chinesis)细胞中克隆了TcNPR1基因,该基因完整的开放阅读框大小为1857 bp,编码蛋白序列由619个氨基酸组成,含有典型的NPR1保守结构域、锌指结构域和锚蛋白结构域。进化树分析发现,TcNPR1蛋白与拟南芥AtNPR3和AtNPR4类蛋白的亲源关系较近。TcNPR1基因表达响应SA、干旱和NaCl胁迫处理,且随着处理时间增加,TcNPR1基因表达量逐渐增大。超表达TcNPR1的转基因烟草植株对干旱的耐受力增强,叶片内可溶性糖含量达43.8 mg/g(FW),比对照高出2倍左右,说明TcNPR1基因具有表达功能,此结果为探究中国红豆杉细胞响应SA诱导大量合成紫杉醇的信号途径提供了依据。
    Abstract: Salicylic acid (SA), a common plant hormone, can increase the biosynthesis of various compounds, such as taxol in Taxus cells. NPR1 is the key gene in SA mediated plant system acquired resistance (SAR) signal pathway. To clarify the regulating mechanism of SA induced taxol biosynthesis, the TcNPR1 gene was isolated from Taxus cells, and was found to contain a 1857 bp ORF (open reading frame), which may encode a 619 aa protein. Sequence alignment showed TcNPR1 contained typical conserve domains such as Broad-Complex, Tramtrack, and Bric-a-brac/Pox virus and Zinc Finger (BTB/POZ domain), NPR1/NIM1-like-C terminal and ankyrin repeats. Phylogenetic analysis showed that TcNPR1 had a close relationship with AtNPR3 and AtNPR4 of Arabidopsis. Expression analysis revealed that the expression of the TcNPR1 gene was significantly induced by SA, drought and NaCl. In comparison to the wild type, transgenic tobacco plants over-expressing TcNPR1 showed elevated drought tolerance with soluble sugar increasing to 43.8 mg/g(FW) after 4 d of drought treatment, about two times more than the control, revealing that TcNPR1 was a functional protein. The available evidence verified that TcNPR1 may respond to SA and may be involved in systemic acquired resistance (SAR) in Taxus species. Further studies on the TcNPR1 gene are essential for clarifying its function in taxol biosynthesis.
  • 人类活动的影响使大气、海洋和陆地显著变暖[1]。森林约占陆地面积的三分之一,是陆地生态系统的重要组成部分,其固定的碳约占陆地生物圈碳储存量的一半[2]。更重要的是,森林生态系统通过吸收和储存大气CO2,在碳汇方面发挥重要作用,对于缓解全球变暖具有不可忽视的意义。以往研究一般使用净生态系统生产力(NEP)表示生态系统的碳净储存量,反映森林生态系统与大气系统之间的CO2交换量[3]。NEP对气候变化敏感,会随着气候因素的变化发生时空异质性的改变[4]。然而,未来气候变化条件下,森林NEP如何响应尚不清楚。因此,通过实验观测或模型模拟来准确地评估森林生态系统NEP与气候的关系对于陆地生态系统碳循环研究具有重要意义。

    目前,国内外开展森林生态系统碳汇估算方法主要有:样地清查法、涡度相关法、碳同化反演法和模型模拟法[5]。近年来,基于过程的生态模型模拟方法,具有可定量区分不同驱动因子对陆地碳汇变化的贡献,可预测陆地碳汇未来变化,并可开展长时间序列指定区域生态系统碳汇评估等优点,已成为区域或全球森林碳汇评估的主要手段。目前已有多种估算陆地生态系统NEP时空变化的生态模型,不同模型在原理、结构和参数上存在差别,导致在评估中存在很大的不确定性[6]。国内外对NEP的研究主要借助CASA和TEC模型,且集中在对其时间动态、影响因素及不确定性的分析,认为气温、降水、太阳辐射以及人类活动等因素是NEP变化的主要驱动因子[7-12]。近20年中国陆地生态系统NEP的年平均值为42.79~61.91 g/m2,呈上升趋势,并表现出明显的空间分异特征,碳汇作用占68.83%,主要集中在东北、华东、华南和华中地区[9]。中国陆地生态系统的整体碳汇能力呈增加趋势,各区域时间变化不一[10]。如,榆林地区荒漠生态系统的碳汇能力在2000-2020年呈增加趋势,2001年之前为碳源,2001年后变为碳汇;2000-2020年碳储量呈增加趋势[13];渭河流域在2001-2020年间也呈显著增加趋势[8]。2000年以来黄河流域年均NEP为92.7 g/m2,总体上表现为碳汇;NEP总体呈增加趋势,其中62.4%的地区增加趋势达到显著水平,植被固碳能力提升明显[11]。从气象因子来看,温度、降水和太阳辐射是影响NEP的主要驱动因子[14]。据报道,荒漠生态系统21年来碳汇能力的显著增强主要是由降水引起的[12]。近20年来,我国陆地生态系统NEP碳汇能力持续上升,降水与太阳辐射是其主要影响因素[9]。曹云等[11]在对黄河流域净生态系统生产力驱动因素的分析中发现,NEP多与降水呈正相关,与日照时数呈负相关,且降水的影响范围超过日照时数。

    由于综合校准和验证数据的可用性不足,关于气候对森林碳循环影响的大规模研究往往依赖于生态系统过程模型,从而将经验模型限制在野外测量的地块级别上,在区域或全球尺度上具有不确定性[15, 16]。考虑到植被物候和相关生长过程,基于叶物候和生长过程的森林生态系统碳通量过程模型FORCCHN2(FORest Ecosystem Carbon Budget Model For CHiNa version 2.0)被广泛应用。在北美进行测试时,与其他6个模型相比,该模型降低了不确定性,且提高了GPP的预测性能[17]。该模型已广泛应用于区域森林生态系统的碳收支模拟,且模拟效果得到了观测数据验证[18-20]

    本研究以南水北调中线工程水源地(以下简称水源地)为对象,研究区包括退耕还林(草)计划、天然林保护以及水土保持等生态建设工程[21],生态服务功能突出。森林覆盖率高,主要以落叶阔叶林和年轻森林为主,对该地区水源涵养能力和生态安全具有重要作用。目前在水源地已经开展了部分研究,但主要集中在植被覆盖变化[22]、径流变化[23]、土地利用格局[24]、水土流失效应[25]、水质变化[26]等方面,关于森林固碳功能变化及其气候响应机制的报道很少。本文利用第二代基于个体树木的森林生态系统碳循环模型FORCCHN2.0估算水源地2002-2020年植被的NEP、GPP、土壤呼吸(Soil Respiration,SR)、植物呼吸(Vegetation Respiration,VR),并分析森林NEP的时空变化特征,同时,探讨气温、降水、辐射等气候因子对NEP、GPP、SR、VR的影响,以揭示森林生态系统生产力变化的影响因素。研究结果旨在为水源地森林生态系统的固碳能力评估以及应对气候变化的政策制定提供科学依据。

    南水北调中线工程水源地位于(31°70′~34°10′N,105°50′~112°10′E)秦岭山脉与大巴山山脉之间,涵盖湖北西北部、陕西东南部、河南西南部等3省44个县市,总面积9.52万km2,海拔为100~3 177 m[27]图1)。研究区属北亚热带季风气候区,年平均降水量923 mm,降水年际变化较大,空间分布不均匀,75%的降水发生在5−10月[23]。年平均气温为14 ℃~16 ℃,最冷季和最热季平均温度分别为2.4 ℃和28.4 ℃。研究区主要包括汉江和丹江两大流域,地貌由山地、丘陵、平原组成[28]。水源地西部和南部植被以森林为主,自然植被主要有常绿针叶林、落叶阔叶林和混交林等,植被覆盖度较高,是南水北调中线工程的水源涵养区[21]

    图  1  研究区域高程(A)、土地覆盖类型(B)与树龄(C)
    非森林区域被掩膜为白色。底图审图号为GS(2020)4619,下同。
    Figure  1.  Elevation (A), land cover types (B), and forest age (C) of the study area
    Non-forest cover areas are masked in white. Map review number is GS (2020) 4619, same below.

    土地覆盖使用MODIS的MCD12Q1数据,该数据集包含17个主要土地覆盖类型,包括11个自然植被类型,3个土地开发和镶嵌地类,以及3个非草木土地类型定义类(https://modis.gsfc.nasa.gov/data/dataprod/mod12.php)。数字高程模型数据(DEM)来源于谷歌地球引擎(GEE)NASA NASADEM Digital Elevation 30 m产品(https://code.earthengine.google.com/9de757ff2bb2ef83b79db2d1181ef544)。叶面积指数(LAI)数据提取自全球陆地表面卫星(GLASS)产品(http://www.glass.umd.edu/Download.html)。气象数据包括每日最高、最低和平均气温(℃)、降水量(mm)、风速(m/s)、太阳总辐射量(W/m2)、相对湿度(%)和大气压(hPa),数据来源于ERA5数据集(https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-land?tab=overview)。土壤质地类型数据来源于世界土壤数据库(HWSD,Harmonized World Soil Database),该数据提供了各个格网点的土壤类型、土壤相位、土壤理化性质等信息(http://www.fao.org/soils-portal/soil-survey/soil-maps-and-databases/harmonized-world-soil-database-v12/en/)。通量验证数据来自秦岭森林生态系统国家野外科学观测研究站(50°56′16″N,121°30′38″E,海拔2 150 m),该站代表了暖温带云杉林生态类型[29]。林龄来自于中国2020年30 m林龄数据集[30]。为方便后续分析,所有数据(除DEM)均重采样至0.05°地理分辨率。

    本文使用的FORCCHN2为基于个体生长过程的中国森林生态系统碳收支模型FORCCHN的第二代模型,耦合了树木的非结构性碳水化合物库(Non-structural carbohydrates,NSC),并将叶物候和相关生长过程加入到了框架内[10]。该模型由每日气象条件驱动逐个计算一定面积斑块上每株树木的碳收支,通过求和及耦合土壤碳循环模型所计算的土壤碳收支,得到生态系统在单位面积上的碳收支。模型建立在两个不同的时间尺度上,包括每日和年度过程。日常过程由光合作用、植物呼吸、光合产物分配、凋落物产生以及土壤呼吸和转移组成[31]。模型假设每日光合作用所产生的碳用于叶和细根的生长,呼吸作用所需消耗的碳以及多余碳的进行储存。光合作用采取与光照条件有关的光饱和曲线进行模拟,生长速率由静态参数决定,呼吸作用包括维持呼吸(RM)和生长呼吸(RG)(即VR=RG+RM)。年度过程则分为林分同化、树高和胸径形式的生长以及粗木屑的生产。上一年储存的碳和当年树木碳的年净增加量用于木质部分茎和枝的尺寸生长,以维持树木对有利资源的竞争,剩余部分的碳被存入NSC库。当产生和储存的碳不足以树木的生长和呼吸时,模型认为树木已经死亡。死亡的树木最终通过凋落作用加入土壤枯枝落叶库中。FORCCHN2模型计算碳收支的主要过程和计算流程如下:

    GPPi=min(fc×fdry×fT×GPPMi,an×aNS)

    式中,GPPi为第i棵树的个体总光合(kg C/d);GPPMi为树木的个体最大总光合;fcfdryfT、an×aNS分别为CO2浓度、土壤水分条件、气温及土壤有效氮对总光合的影响,其中aNS为土壤有效氮量(kg/m2),an为同化作用碳氮比参数,an=150。

    GPPMi=2×Amj×DLKljln1+1+Klj×Slj×PARiAmj1+Klj×Slj×PARi×exp(Klj×LAIi)Amj

    式中,DL为当天可照时数;PARi为该树木的个体正午冠顶光合有效辐射;Amj为第j类树木类型的单叶最大可能光合;Klj为第j类树木类型的消光系数,无量纲;Slj为第j类树木类型的光强-光合作用曲线初始斜率;LAIi为该树木的叶面积指数。

    单株树木的日维持呼吸由叶、枝、干、根和细根的维持呼吸构成。叶、枝、干呼吸构成地上部分呼吸量,根和细根呼吸构成地下部分呼吸量。它们的日呼吸量(kg/d)表示为:

    RMik=tresp×rkCik

    式中,i为第i棵树,k为各个器官,tresp表示气温对植物呼吸的影响系数,rk为在气温为15 ℃时,叶、枝、干、根和细根的相对呼吸速率;Cik为相应碳库的库量。

    生长呼吸指的是初级光合产物转化为结构干物质时的消耗,在树木碳平衡中起重要作用,在日光合作用大于维持呼吸时才有净积累。本模型采用与总光合-维持呼吸成正比的生长呼吸计算公式,即单株树木的生长呼吸为:

    RGi=tresp×rg×(GPPiRMi)

    式中,rg是生长呼吸系数,为一常数,一般取值0.2~0.3,本模型取0.25,即每合成1 g干物质需要1.38 g的葡萄糖。

    土壤呼吸(SR)过程逐日模拟,采用改进后的土壤碳收支模型CENTURY[32]进行模拟和预测。

    NEP的计算取决于GPP、生长呼吸(RG)、维持呼吸(RM)以及土壤呼吸(SR):

    dNEPdt=dGPPidtdRMidtdRGidtdSRdt

    式中,i是地块中第i棵树,t为时间(d)。

    其他详细模型计算方程参考Fang等[17]的研究。

    结合Sen斜率估计法(Theil-Sen)趋势分析和趋势检验法(Mann-Kendall检验)分析NEP增减的空间格局。采用Theil-Sen中值趋势分析方法和Mann-Kendall检验来解释NEP的年际变化趋势。Sen斜率估计是一种非参数统计的趋势计算方法,适用于长时间序列数据的趋势分析。Mann-Kendall检验是一种非参数的时间序列趋势性检验方法,不需要测量值服从正态分布,不受缺失值和异常值的影响,适用于长时间序列数据的趋势显著检验。两种方法相结合已成为确定长期时间序列数据趋势的重要手段[33, 34]

    采用相关分析方法研究植被与气候因子的关系。计算2002-2020年水源地NEP与对应年份日均总辐射、年降水量和年平均气温的相关系数,并进行显著性检验,反映不同气候因子与NEP的相关程度[35, 36]

    为评估模型对森林生态系统碳通量的模拟性能,本文分别采用模拟结果与观测站点观测值的相关系数(R)和均方根误差(RMSE)来进行分析。

    将本研究中FORCCHN2模型的预测结果与秦岭生态站2010-2013年的通量数据进行对比(图2),结果显示,GPP、ER(生态系统呼吸,SR+VR)、NEP的模拟值与实测值均显著相关(P<0.001),说明FORCCHN2模型对水源地NEP、GPP、ER具有良好的模拟能力。

    图  2  总初级生产力(A)、生态系统呼吸(B)、净生态系统生产力(C)模拟值与观测值比较
    GPP:总初级生产力;ER:生态系统呼吸;NEP:净生态系统生产力;OB:观测值;SIM:模拟值。
    Figure  2.  Comparison between simulated and observed values of gross primary productivity (A), ecosystem respiration (B), and net ecosystem productivity (C)
    GPP: Gross primary productivity; ER: Ecosystem respiration; NEP: Net ecosystem productivity; OB: Observed value; SIM: Simulated value.

    图3所示,水源地GPP年均值为1 320.75 g/m2,最低值和最高值分别为2003年的1 177.51和2013年的1 432.98 g/m2。VR年均值为596.94 g/m2,2003年最低,2013年最高。2002-2020年,水源地SR增加趋势显著(P<0.001),增速为4.21 g·m−2·a−1,年均值为235.13 g/m2,2003年最低,2019年最高,分别为196.25、266.29 g/m2。水源地NEP年际变化主要由GPP与ER决定,由于研究期间ER增加的速率大于GPP,导致NEP呈下降趋势(P<0.05),下降速度约为8.85 g·m−2·a−1,年平均值为488.68 g/m2,总的来说,水源地NEP大于0,仍为碳汇功能区域。

    图  3  2002-2020年水源地GPP(A)、植被呼吸(B)、土壤呼吸(C)、NEP(D)年际变化特征
    Figure  3.  Interannual variation characteristics of GPP (A), vegetation respiration (B), soil respiration (C), and NEP (D) in the water source area from 2002 to 2020

    总体来说,水源地西北部的秦岭山脉地区、西南部的大巴山脉地区为成熟林区域(林龄60年以上),森林NEP显著降低,水源地区域东部的丹江口水库周边地区为低龄林区域(林龄60年以下,图1:C),森林NEP无显著变化(图4)。水源地森林NEP在2002-2020年像元尺度上的平均值为308.88~761.22 g/m2,10年变化趋势为−148.76~14.09 g/m2

    图  4  2002-2020年水源地NEP平均值(A)、10年变化趋势(B)和变化显著性(C)空间分布特征
    S1:极显著减少;S2:显著减少;S3:微减少;S4:无显著变化。
    Figure  4.  Spatial distribution characteristics of multi-year average of NEP (A), decade trend of NEP (B), and significance of changes in NEP (C) in the water source area from 2002 to 2020
    S1: Extremely significant reduction; S2: Significant reduction; S3: Slight decrease; S4: Not significant.

    气象因子的年际变化特征如图5所示。总体上,2002-2020年平均气温在11.93 ℃~13.44 ℃,2003年最低,2013年最高,年际变化呈显著上升趋势(P<0.05),上升率为0.062 ℃/10 a。年降水量为939.08~1 481.42 mm,最低值出现在2016年,最高值为2003年,年际变化不显著。日均总辐射为122.85~145.78 W/m2,最低值在2003年,最高值在2013年,年际变化不显著。

    图  5  2002-2020年水源地气候要素年际变化趋势
    P:年总降水;T:年均温;R:日均总辐射。
    Figure  5.  Interannual variation trends in climate elements in the water source area from 2002 to 2020
    P: Annual total precipitation; T: Mean annual temperature; R: Daily total radiation. Same below.

    本研究对水源地NEP与气候要素进行了像元尺度相关分析,并进行了显著性检验,结果如图6所示。分析表明,日均总辐射和年均温分别与水源地98.41%(主要位于西部与西南部)和19.55%区域的NEP呈显著正相关,而年降水与61.57%区域的NEP呈显著负相关。GPP与日均总辐射和年均温呈显著正相关,通过显著性检验(95%)的区域分别占99.71%和98.38%,与日均总辐射相关系数较高的区域位于西南部与西北部地区,而GPP与年降水在97.20%的区域呈显著负相关。SR与57.53%区域的年均温呈显著正相关(主要位于东部),与26.16%区域的年降水呈显著负相关。日均总辐射和年均温分别与83.53%和92.37%区域的VR呈显著正相关,而年降水与83.18%区域的VR呈显著负相关。同时,对气候要素进行年空间变化Theil-Sen趋势分析(图7)。结果表明,年均温的空间变化格局为水源地东部升高速率比西部快;年总降水变化为中部地区、东北部地区呈下降趋势,而南部与西部地区为增加趋势;日均总辐射的变化格局为东部下降速率大于西北部与西南部。

    图  6  水源地不同指标相关系数的空间分布
    结果已通过95%显著性检验,NEP、GPP、VR与日均总辐射(R)相关分析图通过99.9%的显著性检验。
    Figure  6.  Spatial distribution of correlation coefficients of different indicators in the water source area
    Results have passed the 95% significance test. Correlation analysis figures of NEP, GPP, and VR with daily total radiation R have passed the 99.9% significance test.
    图  7  水源地年均温(A)、年总降水(B)和日均总辐射(C)的空间变化趋势
    Figure  7.  Trends of mean annual temperature (A), annual total precipitation (B), and daily total radiation (C) in the water source area

    本研究分析了2002-2020年SR及GPP年际变化趋势的空间格局(图8)。结果显示,水源地东部丹江口水库周边区域的SR增加速度高于GPP;西北与西南部秦岭、大巴山脉区域的GPP呈下降趋势。

    图  8  水源地GPP(A)和SR(B)的年际变化趋势
    Figure  8.  Annual trends of GPP (A) and SR (B) in the water source area from 2002 to 2020

    FORCCHN2模型在全球以及北美森林已被广泛使用[18-20, 35],本研究利用该模型模拟了水源地森林NEP的变化,将模拟数据与野外生态站观测数据对比发现,二者结果较为一致,表明该模型的模拟结果可靠。

    本研究发现,2002-2020年间,水源地森林的碳汇能力呈下降趋势,且具有空间异质性。水源地NEP整体呈显著下降趋势,可能是土壤呼吸增加所导致的。由于水源地处于我国退耕还林政策先行试验区,丹江口水库周边区域的植被在研究期间处于幼龄林到中龄林状态,土壤裸露面积较大,水分散失较快,土壤呼吸较高,从而导致生态系统的碳汇能力降低[37]。相关分析表明,土壤呼吸与温度相关显著区域位于丹江口水库周边区域,与该结果一致。因此,在全球气候变暖背景下,森林树木生长与气候变化使得土壤呼吸的增长比GPP快,但随着树木的生长,土壤呼吸增长减慢,GPP的增加快于生态系统的呼吸作用,水源地的碳汇能力也会逐渐增加。

    此外,水源地NEP的分布及其受气候因子的影响,具有空间差异性。研究结果发现,秦岭和大巴山脉森林的NEP下降比丹江口水库周边区域更为明显,这可能与海拔有一定关系[38]。不同海拔的水热环境、太阳能量通量、CO2浓度等均不同,这些都会对植被NEP产生影响。逐像元相关分析结果表明,日均总辐射和年均温分别与水源地98.41%和19.55%区域的NEP呈显著正相关,而年降水与61.57%区域的NEP呈显著负相关,这说明NEP受太阳辐射的影响最大,其次是降水和温度。这与Wang等[39]对中国森林生态系统的研究结果一致。研究表明,自1999年开始,随着中国退耕还林工程的实施,受降水影响的植被比例逐渐增加,而受气温影响的植被比例下降[40, 41]。研究发现,与丹江口水库周边区相比,秦岭、大巴山区域的太阳总辐射下降速率较快,且GPP在该区域与太阳辐射呈显著正相关,最终导致NEP的下降更明显。前人研究认为,秦岭和大巴山及其毗邻地区的太阳辐射量在年际间呈明显下降趋势,可能是云量和气溶胶的增加所导致[42, 43]。综上所述,水源地NEP受水、光、热条件的共同影响,且气候变化对西北与西南部区域固碳能力的负面影响大于东部丹江口水库的周边区域。

    尽管气候变化在水源地森林生态系统的NEP变化中起重要作用,但森林生态系统固碳受到多种因素的综合影响。如,土地利用变化可能导致森林覆盖的变化,进而影响GPP、ER[44],最终影响NEP的变化特征。未来可进一步提高模型模拟的精度,并加强NEP的多元影响因素分析,如CO2施肥效应、土地利用变化和人为活动影响等。

    2002-2020年间,水源地森林NEP的多年平均值为308.88~761.22 g/m2,固碳能力呈下降趋势,且秦岭和大巴山脉森林的NEP下降快于丹江口水库周边区域。此外,水源地NEP分别与日均总辐射和年均温呈显著正相关,而与年降水呈显著负相关;GPP、SR、VR均与温度及辐射呈显著正相关,与降水呈显著负相关。气候变化对水源地西北、西南地区森林固碳能力的负面影响大于东部。

  • [1]

    Zhang W, Curtin C, Franco C. Towards manipulation of post-biosynthetic events in secondary metabolism of plant cell cultures[J]. Enzyme Microb Tech, 2002, 30(6): 688-696.

    [2]

    Kim BJ, Gibson DM, Shuler ML. Effect of subculture and elicitation on instability of taxol production in Taxus sp. suspension cultures[J]. Biotechnol Prog, 2004, 20(6): 1666-1673.

    [3]

    Croteau R, Ketchum REB, Long RM, Kaspera R. Taxol biosynthesis and molecular genetics[J]. Phytochem Rev, 2006, 5(1): 75-97.

    [4]

    Wilson SA, Roberts SC. Recent advances towards development and commercialization of plant cell culture processes for the synthesis of biomolecules[J]. Plant Biotech J, 2012, 10(3): 249-268.

    [5]

    Wang YD, Wu JC, Yuan YJ. Salicylic acid-induced taxol production and isopentenyl pyrophosphate biosynthesis in suspension cultures of Taxus chinensis var. mairei[J]. Cell Biol Int, 2007, 31(10):1179-83.

    [6]

    Park SW, Kaimoyo E, Kumar D, Mosher S, Klessig DF. Methyl salicylate is a critical mobile signal for plant systemic acquired resistance[J]. Science, 2007, 318(5847): 113-116.

    [7]

    Johnson C, Boden E, Arias J. Salicylic acid and NPR1 induce the recruitment of trans-activating TGA factors to a defense gene promoter in Arabidopsis[J]. Plant Cell, 2003, 15(8):1846-1858.

    [8]

    Zhou J, Trifa Y, Silva H, Pontier D, Lam E, Shah J, Klessig, DF. NPR1 differentially interacts with members of the TGA/OBF family of transcription factors that bind an element of the PR-1 gene required for induction by salicylic acid[J]. Mol Plant Microbe In, 2000, 13(2):191-202.

    [9]

    Mou Z, Fan W, Dong X. Inducers of plant syste-mic acquired resistance regulate NPR1 function through redox changes[J]. Cell, 2003, 113(7):935-944.

    [10]

    Fu ZQ, Yan SP, Saleh A, Wang W, Ruble J, Oka N, Mohan R, Spoel S, Tada Y, Zheng N, Dong XN. NPR3 and NPR4 are receptors for the immune signal salicylic acid in plants[J]. Nature, 2012, 486(7402): 228-232.

    [11] 郭丽红, 陈善娜, 龚明. 根癌农杆菌法转化烟草的条件探索[J]. 云南大学学报:自然科学版, 2003, 25(2): 148-152.
    [12]

    Caldwell KS, Michelmore RW. Arabidopsis thaliana genes encoding defense signaling and recognition proteins exhibit contrasting evolutionary dynamics[J]. Genetics, 2009, 181(2):671-684.

    [13]

    Zhang Y, Cheng YT, Qu N, Zhao QG, Bi DL, LX. Negative regulation of defense responses in Arabidopsis by two NPR1 paralogs[J]. Plant J, 2006, 48(5): 647-656.

    [14]

    Ha CM, Jun JH, Nam HG, Fletcher JC. Blade-on-petiole1 encodes a BTB/POZ domain protein required for leaf morphogenesis in Arabidopsis tha-liana[J]. Plant Cell Physiol, 2004, 45(10):1361-1370.

    [15]

    Hepworth SR, Zhang Y, McKim S, Li X, Haughn GW. Blade-on-petiole-dependent signaling controls leaf and floral patterning in Arabidopsis[J]. Plant Cell, 2005, 17(5): 1434-1448.

    [16]

    Norberg M, Holmlund M, Nilsson O. The blade on petiole genes act redundantly to control the growth and development of lateral organs[J]. Development, 2005, 132(9): 2203-2213.

    [17]

    Bork P. Hundreds of ankyrin-like repeats in functionally diverse proteins: mobile modules that cross phyla horizontally[J]. Proteins, 1993, 17(4): 363-374.

    [18]

    Sedgwick SG, Smerdon SJ. The ankyrin repeat: a diversity of interactions on a common structural framework[J]. Trends Biochem Sci, 1999, 24(8): 311-316.

    [19]

    Niki T, Mitsuhara I, Seo S, Ohtsubo N, Ohashi Y. Antagonistic effect of salicylic acid and jasmonic acid on the expression of pathogenesis-related (PR) protein genes in wounded mature tobacco leaves[J]. Plant Cell Physiol, 1998, 39(5): 500-507.

    [20]

    Feys BJ, Parker JE. Interplay of signaling pathways in plant disease resistance[J]. Trends Ge-net, 2000, 16(10): 449-455.

计量
  • 文章访问数:  1404
  • HTML全文浏览量:  0
  • PDF下载量:  1835
  • 被引次数: 0
出版历程
  • 收稿日期:  2014-05-26
  • 修回日期:  2014-06-14
  • 网络出版日期:  2022-11-01
  • 发布日期:  2014-08-29

目录

/

返回文章
返回