白桦<i>BpSPL</i>6基因启动子的克隆及表达分析

李豆; 苏功博; 胡晓晴; 宋婷婷; 孙庆斌; 徐昭; 王宏伟; 刘雪梅

doi:10.12171/j.1000-1522.20200174

白桦BpSPL6基因启动子的克隆及表达分析

李豆^1,,
苏功博¹,
胡晓晴¹,
宋婷婷¹,
孙庆斌¹,
徐昭¹,
王宏伟^1, ,,
刘雪梅^{1, 2, ,}

1.
东北林业大学生命科学学院，黑龙江哈尔滨 150040
2.
林木遗传育种学科创新引进基地（“111基地”），黑龙江哈尔滨 150040

基金项目: 黑龙江头雁创新团队计划（林木遗传育种创新研究团队），中央高校基本科研业务费专项基金E类项目（2572015EA05）

详细信息

作者简介:
李豆。主要研究方向：林木遗传育种。Email：1171057158@qq.com 　地址：150040 黑龙江省哈尔滨市香坊区和兴路26号东北林业大学生命科学学院

责任作者:
王宏伟，博士，副教授。主要研究方向：生物活性物质制备与应用。Email：84970486@qq.com　地址：同上

刘雪梅，博士，教授。主要研究方向：林木遗传育种。Email：695898040@qq.com　地址：同上

中图分类号: S718.46
计量
- 文章访问数: 2350
- HTML全文浏览量: 827
- PDF下载量: 141
出版历程
- 收稿日期: 2020-06-04
- 修回日期: 2020-07-15
- 网络出版日期: 2021-05-06
- 发布日期: 2022-02-24

Cloning and expression analysis of BpSPL6 promoter from Betula platyphylla

1.
College of Life Science, Northeast Forestry University, Harbin 150040, Heilongjiang, China
2.
Overseas Expertise Introduction Center for Discipline Innovation (“111Center”), Harbin 150040, Heilongjiang, China

摘要

摘要:
目的 SPL（SQUAMOSA promoter binding protein-like）是植物特有的转录因子，参与植物幼年期向成年期的转变、营养生长向生殖生长的转变、花发育、孢子发生、叶片和根发育、逆境响应等多个过程，在植物的生长发育过程中起着非常重要的作用。探究白桦中BpSPL6基因启动子区的顺式作用元件，以及该启动子在正常和胁迫条件下的表达模式，可为进一步研究BpSPL6基因的功能提供参考，也可为了解白桦的抗逆机制提供依据。
方法以本实验室组培白桦的总DNA为模板，经PCR克隆了BpSPL6基因上游1 703 bp的启动子序列，用PLACE和Plant CARE在线软件分析启动子区的顺式作用元件。构建BpSPL6基因启动子驱动GUS报告基因的植物表达载体并转化拟南芥，探究其组织表达特性和胁迫条件下的表达模式。
结果 PCR成功克隆了BpSPL6基因上游1 703 bp的启动子序列，对启动子区的顺式作用元件预测发现除了含有核心启动元件TATA-box、CAAT-box外，还包括2种特异组织表达元件（根、花粉），10种激素响应元件（生长素、赤霉素、水杨酸、脱落酸），4种脱水响应元件等。对转基因拟南芥进行GUS染色结果表明，BpSPL6基因启动子驱动的GUS基因在转基因拟南芥中的表达具有时空特异性。在拟南芥的整个发育过程中，BpSPL6基因启动子驱动GUS基因在真叶叶片中表达，但是表达部位不同。随着叶片的生长，首先在叶片的顶端表达，随后扩展到叶片的叶脉并直至整个叶片，并且表达量逐渐升高。同时BpSPL6基因启动子驱动的 GUS 基因在拟南芥营养生长时期的根部都有表达。并且经氯化钠和甘露醇胁迫后其表达量降低。对比两种胁迫，受到氯化钠胁迫后GUS基因的表达量变化更大，说明对氯化钠胁迫的响应更加强烈。
结论 BpSPL6基因可能参与了植物的叶片、根发育以及对盐和干旱胁迫的响应。
- 白桦 /
- SPL /
- 启动子 /
- 顺式作用元件 /
- GUS染色 /
- 盐胁迫 /
- 干旱胁迫
Abstract:
Objective SPL (SQUAMOSA promoter binding protein-like) is a plant-specific transcription factor and plays crucial roles in the growth and development of plants. It participates in multiple processes such as plant phase transformation, flower development, sporogenesis, leaf and root development, and stress response. Exploring the cis-acting elements of the promoter region of the BpSPL6 gene, as well as the expression pattern under normal and stress conditions, this paper provides a reference for further study of the function of BpSPL6 gene and also provides a basis for understanding the stress resistance mechanism of Betula platyphylla.
Method Using the total DNA of tissue-cultures B. platyphylla in this laboratory as a template, the 1 703 bp promoter sequence of BpSPL6 gene was cloned by PCR. The cis-elements of the promoter region were analyzed using PLACE and Plant CARE web tools. We constructed the BpSPL6 gene promoter-driven GUS expression vector, which was transformed into Arabidopsis thaliana. Transgenic A. thaliana was subjected to stress experiments, and GUS staining was used to analyze the expression pattern of the BpSPL6 gene promoter.
Result A 1 703 bp promoter sequence of BpSPL6 gene was cloned from B. platyphylla genomic DNA. Sequence analysis showed that the BpSPL6 gene promoter included core promoter elements TATA-box and CAAT-box, as well as cis-elements for specific parts (root and pollen), hormone response elements (auxin, gibberellin, salicylic acid, abscisic acid) and dehydration response elements. The GUS staining results of transgenic A. thaliana showed that the expression of GUS gene driven by BpSPL6 gene promoter had spatiotemporal specificity. Throughout the development of A. thaliana, the GUS gene driven by BpSPL6 gene promoter was expressed in euphylla, but the expression sites were different. As the leaf grew, it was first expressed at the tip of the leaf, then expanded to the leaf vein and reached the entire leaf, and the expression level gradually increased. At the same time, the GUS gene driven by BpSPL6 gene promoter was also expressed in the roots during vegetative growth. When transgenic A. thaliana was subjected to sodium chloride and mannitol stress, the expression level of GUS gene driven by BpSPL6 gene promoter decreased compared with control. Comparing the two stresses, the expression level of GUS gene changed more after being subjected to sodium chloride stress, indicating that the response to sodium chloride stress was stronger.
Conclusion BpSPL6 gene may be involved in plant leaf and root development and response to salt and drought stress.
- Betula platyphylla /
- SPL /
- promoter /
- cis-element /
- GUS staining /
- salt stress /
- drought stress

HTML全文

森林是全球三大生态系统之一，属于陆地生态系统的重要组分，是极为重要的自然资源，同时在水土保持、改善气候、维护生物多样性及碳汇等方面起着重要作用，是应对全球气候变化的指示剂和调节剂，具有特殊的生态效益^[1-4]。IPCC第四次评估报告指出随着全球工业化进程的加快，气候变暖对森林生长的影响日渐加剧^[5]。探究气候变化背景下森林生长的变化，对评估森林生态系统服务功能、预测未来森林演替方向以及科学指导营林活动具有重要意义^[6-7]。

研究表明北半球中高纬度地区变暖速率远远高于低纬度地区^[8-9]，已有研究利用野外调查、模型分析等手段发现气候变暖将对中高纬度森林产生重要影响^[10]。以往研究多采用模型估测^[11]或野外实地采取树木年轮^[12]来研究气候变化对中高纬度地区树木生长的影响，但模型的结论精度较低，而利用树木年轮的研究尺度较小。随着技术进步，卫星遥感已成为大时空尺度反映植被动态变化及监测其对气候变化响应特征的主要技术手段之一^[2,13]。国内外广泛利用归一化植被指数（NDVI）在大尺度范围反映植被生长情况及其对全球变化的响应，已获得良好效果^[14-15]。我国东北地区属于全球中高纬度变暖较显著区域之一，拥有独特的植被分布格局，其森林分布面积广泛，且对气候变化拥有较强敏感性^[10,16]。此前已有研究探究了东北植被NDVI变化及其与气候因子的关系，并取得一定进展：刘向培等^[17]结合1982—2006年GIMMS NDVI数据，发现东北地区植被年均NDVI在东北部高于西南部；神祥金等^[16]发现东北植被年均NDVI受降水和温度的影响均较为显著，但温度影响更强；史娜娜等^[18]采用相关分析发现东北森林区植被生长对降水的响应存在时滞性；罗新兰等^[19]利用标准化降水蒸散指数（SPEI）探究了东北植被NDVI对干旱的响应，得出森林NDVI同SPEI大体呈负相关关系；王小霞等^[2]分析了人类活动对东北森林NDVI的影响，认为人类活动对森林NDVI趋势的影响主要集中于2 km的缓冲区内。以上研究或是未针对未来东北森林NDVI变化趋势进行预测或是仅利用平均气温研究植被NDVI对气候变化的响应，而未考虑最高气温和最低气温可能对NDVI变化存在不同的作用。已有研究在其他区域发现最高和最低气温对森林生长具有差异性影响，如赵杰等^[15]发现夜晚增温对新疆森林植被活动的影响强于白天增温，Wang等^[20]发现欧洲温带树种春季物候的温度敏感性下降是由于白天增温导致的。白天和夜晚增温速率具有不对称性已成为当前全球气候变化的主要特征之一^[21-22]，因而在此背景下，本文分析比较最高气温和最低气温对东北森林NDVI的作用，以期更加全面地反映气候变化对东北森林的影响。

本文基于2000—2019年SPOT Vegetation NDVI数据和气象数据，分析东北森林生长季（4—10月）^[19]NDVI时空变化特征，探究生长季NDVI与主要气象因子（降水量、平均气温、平均最高气温和平均最低气温）的响应关系及时滞性特征，着重探讨白天、夜晚温度对东北森林生长季NDVI影响是否存在差异性特征，并利用Hurst指数预测东北森林未来生长趋势，为东北森林建设、保护与生态管理提供科学理论与数据。

1. 研究区概况与研究方法

1.1 研究区概况

东北地区位于我国东北部（38°72′ ~ 53°55′N、115°52′ ~ 135°09′E）^[17]，由黑龙江、吉林和辽宁三省以及内蒙古东部三市一盟构成^[11]，总面积约为127万km²。该地区森林总面积约40万km²，主要分布在大兴安岭、小兴安岭以及长白山，其中大部分位于温带大陆性季风气候区，夏季高温多雨、阳光充沛，冬季寒冷漫长，也有部分森林分布在寒温带^[2,23-24]，针叶林、阔叶林和针阔叶混交林是东北森林植被的3种主要森林类型^[25]。东北森林是我国面积最大的天然林区，是重要碳汇区及北方生态屏障，属于三北防护林体系及天然林保护工程的重点领域，具有难以估量的生态效益^[26]。本研究根据中国科学院资源环境科学数据中心（http://www.resdc.cn）提供的中国1∶1 000 000植被类型空间分布数据将东北森林分为针叶林、针阔叶混交林（混交林）和阔叶林进行对比分析（图1）。3种林型空间分布差异明显，其中，针叶林主要分布在大兴安岭，零星分布于小兴安岭和长白山区域；阔叶林广泛分布于小兴安岭、长白山、大兴安岭南部区域和西南部林草交错区；混交林分布范围很小，主要零星分布在小兴安岭和长白山，对本文研究结果影响较少。

图 1 研究区概况

审图号：GS(2021)5448号。下同。Base drawing review No.: GS(2021)5448. The same below.

Figure 1. Survey of the study area

下载: 全尺寸图片幻灯片

1.2 数据来源与处理

本文NDVI数据来源于中国科学院资源环境科学数据中心（http://www.resdc.cn）提供的“中国月度1 km植被指数（NDVI）空间分布数据集”。数据集基于SPOT-VEGETATION NDVI卫星遥感数据，在每10 d数据基础上利用最大合成法合成月数据，空间分辨率为1 km × 1 km，且在此基础上提出2000—2019年共20年东北森林生长季（4—10月）NDVI平均值，并根据前人研究^[27]，将多年生长季NDVI值大于0.1的区域定义为研究区。

本文采用的气象数据来源于中国气象数据网（http://data.cma.cn）提供的2000—2019年东北地区逐月降水量、平均气温、平均最高气温和平均最低气温数据，经过数据质量完整性及一致性检验，剔除存在缺失值和异常值的站点，最终选取97个气象站点数据进行分析。应用ArcGIS软件的插值分析工具，采取克里金插值法^[28-29]，基于各站点经纬度对气象数据进行插值，得到与NDVI数据投影、分辨率一致的栅格气象数据^[18]。

1.3 研究方法

1.3.1 趋势分析

本文利用一元线性回归分析方法逐像元计算2000—2019年东北森林生长季NDVI和各气象要素变化趋势^[30]，趋势计算公式如下^[31]：

${\theta _{{\rm{slope }}}} = \frac{{\left( {n \times \displaystyle \sum\nolimits_{i = 1}^n i \times {\rm{NDV}}{{\rm{I}}_i}} \right) - \left( {\displaystyle \sum\nolimits_{i = 1}^n i\displaystyle \sum\nolimits_{i = 1}^n {{\rm{NDV}}{{\rm{I}}_i}} } \right)}}{{n \times \displaystyle \sum\nolimits_{i = 1}^n {{i^2}} - {{\left( {\displaystyle \sum\nolimits_{i = 1}^n i } \right)}^2}}}$

式中： ${\theta _{{\rm{slope}}}}$ 表示各像元NDVI的变化趋势（斜率）；n为监测年数即为20年；i为年序号；NDVI_i为第i年的生长季NDVI平均值。 ${\theta _{{\rm{slope}}}}$ 为正值，则NDVI的年际变化呈现增加趋势，反之说明其呈现减少趋势。

1.3.2 Hurst指数

Hurst指数是用来定量描述研究对象的自相关性和长期依赖性的一种有效方法^[32]。有多种计算方法用来获取Hurst指数，而目前被认为可靠性最强、应用最广泛的方法是重际极差（R/S）分析法^[33-34]。本文利用基于R/S分析法的Hurst指数研究东北森林生长季 NDVI变化的持续性特征。基本原理如下^[35]：

对时间序列 ${\rm{NDVI}}_t$ ， $t=1,2, \cdots, n$ ，定义均值序列：

${\overline {{\mathop{\rm NDVI}\nolimits} } _{\tau }} = \frac{1}{\tau }\sum\nolimits_{t = 1}^\tau {{\rm{NDV}}{{\rm{I}}_{\tau }}},\;\;\;\;\;\tau = 1,2, \cdots,n$

积累离差：

${X_{(t,\tau )}} = \sum\nolimits_{t = 1}^{\tau} { {{\rm{NDV}}{{\rm{I}}_{t}}}} - {\overline {{\rm{NDVI}}} _{\tau }},\;\;\;\;1 \leqslant t \leqslant \tau$

极差：

${R_{\tau }} = \mathop {\max }\limits_{1 \leqslant t \leqslant \tau } {{{X}}_{(t,\tau )}} - \mathop {\min }\limits_{1 \leqslant t \leqslant \tau } {{{X}}_{(t,\tau )}},\;\;\;\;\tau = 1,2, \cdots,n$

标准差：

${S_{\tau }} = {\left[ {\frac{1}{\tau }\sum\nolimits_{t = 1}^\tau {{{\left( {{\text{NDV}}{{\text{I}}_{t}} - {\text{NDV}}{{\text{I}}_{\tau }}} \right)}^2}} } \right]^{\frac{1}{2}}},\;\;\;\;\tau = 1,2, \cdots,n$

若存在 $R/S \propto {\tau ^H}$ ，则说明该时间序列存在Hurst现象，Hurst指数即H值可通过双对数坐标系用最小二乘法拟合计算求得。当0.5 < H < 1，表明该时间序列具有正向持续性，即未来变化状况与过去20年变化趋势一致，且H值越接近1，正向持续性越强；当H = 0.5时，表明该时间序列为一个随机时间序列，即未来变化状况与过去变化趋势无关；当0 < H < 0.5，表明该时间序列具有反向持续性，即未来变化状况与过去20年变化趋势相反，且H值越接近0，反向持续性越强。

1.3.3 相关分析

运用相关分析方法来研究东北森林生长季NDVI和各气象要素之间的关系^[36]，从而分析东北森林生长季NDVI对气候变化的响应。计算公式如下^[37]：

${R_{xy}} = \dfrac{{\displaystyle \sum\nolimits_{i = 1}^n {\left( {{x_i} - \bar x} \right)} \left( {{y_i} - \bar y} \right)}}{{\sqrt {\displaystyle \sum\nolimits_{i = 1}^n {{{\left( {{x_i} - \bar x} \right)}^2}} \displaystyle \sum\nolimits_{i = 1}^n {{{\left( {{y_i} - \bar y} \right)}^2}} } }}$

式中：R_xy为相关系数；n指时间序列累计年数即20年；x_i和y_i分别为某年份生长季NDVI和某气象要素的数值； $\overline x$ 和 $\overline y$ 分别表示20年来生长季NDVI和某气象要素的均值。

2. 结果与分析

2.1 东北森林生长季NDVI时空变化特征

2.1.1 NDVI空间分布特征

利用2000—2019年东北森林生长季NDVI数据，计算20年平均NDVI的空间分布图（图2），从图中可以看出，总体上东北森林生长季平均NDVI值集中分布0.40 ~ 0.78，且在长白山、小兴安岭及大兴安岭等地区NDVI值偏高，这些地区都属于森林大面积成片分布区，人为干扰少^[2]；NDVI值较低的区域主要分布在西南部区域，这些区域属于林草交错区，人为干扰较大且自然气候条件较差^[38]，所以NDVI值较低，但该区域森林面积占整体东北森林面积较小。不同林型之间相比，混交林分布区生长季平均NDVI值最高（0.724），针叶林次之（0.655），阔叶林最低（0.647）。

图 2 东北森林生长季NDVI空间分布特征

Figure 2. Spatial distribution characteristics of growing season NDVI in forest zone of Northeast China

下载: 全尺寸图片幻灯片

2.1.2 NDVI变化趋势

东北森林生长季平均NDVI在2000—2019年期间总体呈上升趋势，趋势为每10年增加0.035，表明东北森林生长状况在20年来有着明显提升（图3），这和前人研究基本一致^[23,39]。针叶林、混交林和阔叶林生长季NDVI 均呈现上升趋势，其中混交林NDVI增长趋势最大（每10年增加0.037），阔叶林NDVI增长趋势最小（每10年增加0.035）。各类型森林的最低值均出现在2001年，最高值出现在2014年。2001—2006年、2011—2014年增速较快，2014年后所有林型森林生长季NDVI增长明显放缓。总体来说，这3种林型NDVI时间序列变化规律差异不大，在过去20年呈现出相对一致的增长过程。

图 3 东北森林生长季NDVI变化

Figure 3. Changes of growing season NDVI of forest zone in Northeast China

下载: 全尺寸图片幻灯片

图4表示东北森林生长季NDVI空间变化，根据杨尚武等^[40]提出的划分标准，将东北森林生长季NDVI变化斜率θ_slope划分为5个变化区间，即：中度退化、轻度退化、基本不变、轻度改善、中度改善。由图4和表1可知，在过去20年东北森林整体上以轻度改善为主，其面积为257 103 km²，占总面积的55.22%，其次基本不变的面积为204 138 km²，占总面积的43.84%，上述两种变化区间占总面积的99.06%。基本不变的主要分布在东北森林边缘的过渡带，其他地区主要呈现轻度改善，另外3种变化区间则零星分布，占比极小，这与该区域类似研究^[2,16]的结果基本一致。各林型中，混交林改善面积占自身面积比例最大（65.21%），阔叶林次之（55.62%），针叶林最小（54.32%）；阔叶林退化面积占自身面积比例最大（0.69%），混交林最小（0.27%）。东北森林近20年以改善为主，国家近年来实行的天然林保护工程减少了对东北森林的人为扰动表明这一变化主要由气候因素导致^[23,41]，我们接着试图结合气候变化特征和生长对气候因子响应特征分析东北森林生长变化的原因。

图 4 东北森林生长季NDVI空间变化

Figure 4. NDVI spatial variation in the growing season of forests in Northeast China

下载: 全尺寸图片幻灯片

表 1 2000—2019年东北森林生长季NDVI变化趋势统计

Table 1. Statistics in NDVI changing trend in forest growing season in Northeast China from 2000 to 2019

θ_slope	变化程度 Degree of change	所有森林 All forests		针叶林 Coniferous forest		混交林 Mixed forest		阔叶林 Broadleaved forest
θ_slope	变化程度 Degree of change	面积 Area/km²	比例 Proportion/%	面积 Area/km²	比例 Proportion/%	面积 Area/km²	比例 Proportion/%	面积 Area/km²	比例 Proportion/%
−0.0258 ~ −0.0090	中度退化 Moderately degraded	520	0.11	143	0.11	9	0.05	368	0.12
−0.0090 ~ −0.0034	轻度退化 Lightly degraded	2 141	0.46	325	0.25	37	0.22	1 779	0.57
−0.0034 ~ 0.0034	基本不变 Basically unchanged	204 138	43.84	59 623	45.32	5 852	34.53	137 265	43.70
0.0034 ~ 0.0097	轻度改善 Lightly improved	257 103	55.22	71 020	53.98	11 046	65.17	173 460	55.22
0.0097 ~ 0.0208	中度改善 Moderately improved	1 723	0.37	448	0.34	6	0.04	1 243	0.40

下载: 导出CSV

| 显示表格

2.1.3 NDVI 未来趋势分析

对近20年东北森林生长季NDVI进行Hurst指数分析，结果见图5。东北森林生长季NDVI的Hurst指数介于0.10 ~ 0.98之间，平均值为0.4516，Hurst指数小于0.5的像元占比为71.73%，表明未来东北森林生长季NDVI变化的反持续性强于正持续性，即东北森林生长季NDVI变化趋势在未来将发生逆转。各森林类型Hurst指数由小到大的顺序为：混交林（0.433 9） < 阔叶林（0.447 8） < 针叶林（0.463 0），均在整体上呈现反持续性。

图 5 东北森林生长季NDVI的Hurst指数及未来趋势

Figure 5. Hurst index and future trend of growing season NDVI of forest zone in Northeast China

下载: 全尺寸图片幻灯片

结合东北森林生长季NDVI变化斜率θ_slope以及Hurst指数，东北森林生长季NDVI未来趋势主要有4种情况^[42]：θ_slope小于0同时Hurst指数小于0.5，未来趋势为由退化到改善；θ_slope大于0同时Hurst指数小于0.5，未来趋势为由改善到退化；θ_slope小于0同时Hurst指数大于0.5，未来趋势为持续退化；θ_slope大于0同时Hurst指数大于0.5，未来趋势为持续改善。结合图5和表2可知，未来将出现大面积的森林由改善变为退化，各种森林类型均面临退化趋势，占总面积的71.94%，这同全球范围内出现的北方林森林衰退现象相吻合^[43]。其中混交林区域面临的退化风险最大，占其自身面积的76.23%；针叶林最小，占其自身面积的69.60%。同时，有占自身面积30.4%的针叶林、27.33%的阔叶林和23.77%的混交林未来将改善。气候变暖已经导致世界各地普遍出现森林衰退^[44]，并已有研究开始关注中国北方乃至东北部分区域出现的林木衰退^[45-46]。我国东北地区对气候变化较为敏感^[47]，此前有研究利用模型预测未来气候变化将影响东北地区森林演替^[11,48]，但目前许多机制仍未研究清楚，北方森林生态系统未来变化仍具有不可确定性^[49]，需要加强这一方面的研究为制定针对性策略提供支撑。

表 2 东北森林生长季NDVI的未来趋势统计

Table 2. Statistics of future trend of growing season NDVI of forest zone in Northeast China

未来趋势 Future trend	所有森林 All forests		针叶林 Coniferous forest		混交林 Mixed forest		阔叶林 Broadleaved forest
未来趋势 Future trend	面积 Area/km²	比例 Proportion/%	面积 Area/m²	比例 Proportion/%	面积 Area/m²	比例 Proportion/%	面积 Area/m²	比例 Proportion/%
退化到改善 Degradation to improvement	7 772	1.67	909	0.69	138	0.82	6710	2.14
改善到退化 Improvement to degradation	325 658	70.06	89 843	68.39	12 794	75.59	220 848	70.43
持续退化 Continuous degradation	8 737	1.88	1 583	1.21	108	0.64	7 024	2.24
持续改善 Continuous improvement	122 668	26.39	39 025	29.71	3 885	22.95	78 971	25.19

下载: 导出CSV

| 显示表格

2.2 东北森林生长季NDVI对气象因子的响应

2.2.1 东北森林生长季气象因子的变化特征

2000—2019年间，东北森林生长季主要气象因子变化趋势见表3。生长季降水量呈现显著上升趋势，增长速率为每10年10.37 mm。平均气温和平均最高气温上升趋势不明显，平均最低气温显著上升，这和全球大部分地区昼夜升温不同步的现象相一致^[21]。

表 3 2000—2019年东北森林生长季气象因子的变化趋势

Table 3. Changing trend of meteorological factors in growing season of forests in Northeast China from 2000 to 2019

森林类型 Forest type	降水量/（mm·10⁻¹ a⁻¹） Precipitation/ (mm·10⁻¹ year⁻¹)	平均气温/（℃·10⁻¹ a⁻¹） Mean temperature/ (℃·10⁻¹ year⁻¹)	平均最高气温/（℃·10⁻¹ a⁻¹） Mean max. temperature/ (℃·10⁻¹ year⁻¹)	平均最低气温/（℃·10⁻¹ a⁻¹） Mean min. temperature/ (℃·10⁻¹ year⁻¹)
所有森林 All forests	10.37^**	0.07	0.01	0.23^*
针叶林 Coniferous forest	8.73^*	0.12	0.01	0.35^**
混交林 Mixed forest	10.95^*	0.08	0.07	0.15
阔叶林 Broadleaved forest	11.03^**	0.05	0.01	0.18⁺
注：⁺代表P < 0.10，^代表P < 0.05，^代表P < 0.01。下同。Notes: ⁺ represents P < 0.10; ^represents P < 0.05; ^**represent P < 0.01. The same below.

下载: 导出CSV

| 显示表格

从图6可以看出，最高气温和最低气温变化的存在明显的空间差异特征，所以仅使用平均气温作为温度指标不能全面反映温度变化对植物生长的影响。研究区范围内，生长季降水量呈现增加趋势，且中部和北部区域增加速率强于南部区域，仅辽宁东部部分区域降水量为减少的趋势（图6）。平均最低气温增长趋势明显，但大兴安岭东部、小兴安岭东部和长白山南部存在下降趋势；平均最高气温增加区域和降低区域均有大面积分布，呈增加趋势的区域主要分布在大兴安岭北部和长白山区域，大兴安岭南部和小兴安岭则呈下降趋势；平均气温受最高温和最低温共同控制，在总体上呈现增长趋势，但增加速率低于最低气温。总的来说，2000—2019年东北中北部森林区域生长季暖湿化趋势明显，而南部区域趋势稍弱。

图 6 2000—2019年东北森林生长季气象因子变化趋势空间分布

Figure 6. Spatial distribution of meteorological factors in the growing season of forests in Northeast China from 2000 to 2019

下载: 全尺寸图片幻灯片

2.2.2 生长季NDVI与气候变化的关系

东北森林生长季平均NDVI与气象因子的相关性如表4所示，生长季平均NDVI与降水和平均最低气温总体呈现较强的正相关，表明降水和平均最低气温是影响整体东北森林生长的主要气象因子。影响针叶林生长的主要气象因子是平均最低气温和平均气温，其中最低气温通过了0.10水平的显著性检验，说明降水对针叶林生长的影响较小，东北针叶林生长主要受温度控制，这和郭金停^[39]的研究结果基本一致。针叶林主要分布在中高纬度地区，天气寒冷，升温可能增强其光合作用，促进生长^[2,50]。混交林生长与降水的相关性最高；阔叶林生长同降水和最低气温都具有显著的相关性（P < 0.05），说明其受降水和最低气温共同控制。不同林型的上述差异，可能是因为生理生态特性存在差异及水热条件空间差异所导致的。

表 4 东北森林生长季平均NDVI与气象因子的相关性

Table 4. Correlations between mean NDVI and meteorological factors in the growing season of forests in Northeast China

平均NDVI Mean NDVI	降水量 Precipitation	平均气温 Mean temperature	平均最高气温 Mean max. temperature	平均最低气温 Mean min. temperature
所有森林 All forests	0.44⁺	0.26	0.17	0.47^*
针叶林 Coniferous forest	0.24	0.37	0.30	0.42⁺
混交林 Mixed forest	0.35	0.22	0.11	0.27
阔叶林 Broadleaved forest	0.50^*	0.20	0.10	0.48^*

下载: 导出CSV

| 显示表格

图7为东北森林生长季平均NDVI与生长季降水、平均气温、平均最高气温和平均最低气温的相关系数空间分布，且根据显著性水平为0.05（相关系数为±0.42）阀值划分，可以看出各气象因子对东北森林NDVI的影响存在明显的空间异质性。表5显示出逐像元分类统计结果。

图 7 东北森林生长季NDVI与气象因子相关系数空间分布

Figure 7. Spatial distribution of correlation coefficients between NDVI in growing season and meteorological factors of forest zone in Northeast China

下载: 全尺寸图片幻灯片

表 5 东北森林生长季NDVI与气象因子相关性所占比例

Table 5. Proportion of correlation between NDVI and meteorological factors in the growing season of forests in Northeast China %

气象因子 Meteorological factor	相关性 Correlation	所有森林 All forests	针叶林 Coniferous forest	混交林 Mixed forest	阔叶林 Broadleaved forest
降水 Precipitation	显著负相关 Significantly negative correlation	0.51	0.13	0.63	0.67
	不显著负相关 Not significantly negative correlation	16.39	17.56	12.22	16.18
	不显著正相关 Not significantly positive correlation	69.46	77.90	80.74	65.19
	显著正相关 Significantly positive correlation	13.64	4.41	6.40	17.96
平均气温 Mean temperature	显著负相关 Significantly negative correlation	1.98	0.97	0.24	2.51
	不显著负相关 Not significantly negative correlation	21.24	10.03	16.93	26.24
	不显著正相关 Not significantly positive correlation	62.84	64.80	78.97	61.13
	显著正相关 Significantly positive correlation	13.94	24.20	3.87	10.12
平均最高气温 Mean max. temperature	显著负相关 Significantly negative correlation	2.02	0.95	0.51	2.56
	不显著负相关 Not significantly negative correlation	29.06	13.62	38.76	35.10
	不显著正相关 Not significantly positive correlation	59.61	72.43	52.11	54.48
	显著正相关 Significantly positive correlation	9.31	13.00	8.63	7.86
平均最低气温 Mean min. temperature	显著负相关 Significantly negative correlation	3.36	1.66	0.28	4.25
	不显著负相关 Not significantly negative correlation	21.60	18.03	29.19	22.75
	不显著正相关 Not significantly positive correlation	46.89	55.05	55.30	43.05
	显著正相关 Significantly positive correlation	28.15	25.26	15.22	29.95

下载: 导出CSV

| 显示表格

生长季NDVI与降水呈正相关的区域占整个东北森林面积的83.13%，除长白山南部区域外其他区域均大面积分布，尤其是西南部较为干旱的林草交错区对降水最为敏感。阔叶林分布区的北部区域（大兴安岭南部、小兴安岭和长白山东北部）对降水敏感性较高，说明这些区域的水分条件对制约了当地树木生长；生长季NDVI与降水和最低气温均呈现不显著负相关的区域主要集中在长白山南部，这里是东北森林水热条件较好的地区^[51]，生长季水热供应较为充足。

东北森林生长季NDVI与生长季平均气温、最高气温和最低气温均呈现正相关关系的区域面积占整个森林面积的比例分别为76.78%、68.92%和75.04%。研究区内，生长季平均最高气温和最低气温同NDVI相关系数的空间分布具有明显的空间异质性：最高气温同NDVI呈显著相关的区域主要分布在大兴安岭北部、长白山东北部；最低气温同NDVI呈显著相关的区域主要分布在大兴安岭北部和南部、小兴安岭西部。除占整个森林比例较低的混交林（约占整个森林面积的3.6%）外，阔叶林生长季NDVI同最低气温的相关性通过了显著性水平为0.05的检验，针叶林通过了显著性水平为0.10的检验，这与Peng等^[52]发现北半球干旱半干旱区域广泛存在着生长季NDVI与最低气温具有显著正相关性的结论基本一致。上述现象表明，最高气温和最低气温变化对东北森林生长具有非对称影响，夜晚升温比白天升温对东北森林生长的促进作用更强，这可能是由于夜晚增温促进植物白天光合作用的补偿效应比促进植物夜晚呼吸强度的影响更大^[53]。补偿作用是植物应对逆境的一种适应方式，广泛存在于各类植物中^[54-55]，当植物夜晚呼吸作用增强而消耗过多有机物后，第二天白天通过增强光合作用进而补充更多有机物，使生长得到恢复，并已有研究发现较高的立地条件也可以促进植物的补偿效应^[56-57]。东北土壤肥沃，生长季水热条件较好，可能有利于补偿效应发生。有研究者提出，夜晚升温有利于减少霜冻发生的频率，并因此提升植物的生产力^[15]。同时，昼夜不对称增温对植物地下根系、土壤微生物活性、土壤理化性质等均会产生影响，进而影响物质循环和能量流动乃至植物地上部分的生长^[58-59]。

总体而言，近20年东北森林生长季降水量和最低气温升高与东北森林生长季NDVI的增加相吻合，尤其是中部和北部地区，其相较于南部区域来说对气候响应更敏感，而气候变化也较南部更剧烈；但根据Hurst指数及未来趋势（图5）也可发现，东北森林未来面临衰退的区域主要分布在大兴安岭西部、小兴安岭及长白山东北部等气候变化较为剧烈的区域，结合林型分布推测，这可能是因为气候持续变化下未来水热耦合条件将导致生境不利于当地树木生长需求，导致森林未来可能出现衰退。

2.2.3 逐月NDVI与气象因子的相关性与时滞性分析

树木生长不仅受同期气候因子的影响，同时对其响应还存在一定程度的滞后性^[60-61]。本文将东北森林生长季逐月NDVI与同期及前1—3月的气象因子的响应进行分析，结果见表6。4月NDVI与同期以及前1月的平均气温、最高气温和最低气温均呈现显著正相关性，说明该时期NDVI不仅与同期温度显著相关，对前一月的温度也具有显著的滞后响应；5月NDVI与同期3个温度因子呈现出显著相关性，但已无响应滞后现象；而6月NDVI仅与同期最低气温关系显著，进入7月后则NDVI同温度关系不再显著。说明生长季初期，NDVI对温度的响应较强且存在一定时滞性，而随着温度升高，这种响应滞后现象逐步消失，这可能是因为前期升温有利于积累热量加速冻土、积雪融化，增加水分供应，促进生长季初期植物萌发及进行光合作用^[10,62]，随着温度条件持续向好这种响应滞后逐步消失，6月NDVI只与同期最低气温具有显著相关性，说明此时白天温度已经不再成为影响植物生长的主要温度因子，随着温度继续上升，生长较为旺盛的7—8月^[12]热量供应充足，植物生长对温度不再敏感。根据表7可知，各月温度变化均未达到显著，这说明单独月份的温度长期变化趋势对植被的生长改变较小。而6—8月的NDVI分别同前1—3月（即5月份）的降水呈显著正相关，说明在生长较为旺盛的6—8月，植物生长同生长旺盛期前的水分积累关系密切，同时东北森林地区5月份的降水上升趋势显著，这有助于森林植被生长。8月和9月降水呈显著上升的趋势，同期温度指标变化均不显著，生长季末期9—10月NDVI均和前1个月的最低气温显著正相关，同时和前1月降水的正相关性分别通过了显著性水平为0.05和0.10的检验，降水变化有利于提升NDVI值，这可能是在生长期末期前的低温偏高会导致植物光合作用增强，对水分需求提升^[46,63]，同时推迟植物枯黄的时间，即延长生长期^[64]。由此可见，东北森林NDVI在生长季前期对温度存在滞后响应，生长季中期对降水响应存在滞后现象，而生长季末期对最低气温和降水同时存在时滞性响应。在生长季前期，东北森林植被生长和最高气温及最低气温关系都较为显著，而生长季末期则只与最低气温关系密切。

表 6 东北森林逐月NDVI与气象因子的相关性

Table 6. Correlations between monthly NDVI and meteorological factors of forest zone in Northeast China

气象因子 Meteorological factor	项目 Item	4月 Apr.	5月 May	6月 Jun.	7月 Jul.	8月 Aug.	9月 Sep.	10月 Oct.
降水 Precipitation	同期 Same period	−0.351	0.136	0.098	−0.014	0.241	−0.321	−0.190
	前1月 Previous one month	−0.236	−0.059	0.683^**	0.261	−0.187	0.489^*	0.420⁺
	前2月 Previous two months	−0.231	−0.109	−0.275	0.608^**	0.200	−0.128	0.223
	前3月 Previous three months	−0.221	0.180	−0.141	−0.235	0.601^**	−0.078	0.042
平均气温 Mean temperature	同期 Same period	0.835^**	0.518^*	0.180	0.296	0.201	−0.131	0.179
	前1月 Previous one month	0.576^**	0.126	−0.201	−0.063	0.267	0.069	0.300
	前2月 Previous two months	0.239	0.072	0.014	−0.110	0.009	0.111	−0.074
	前3月 Previous three months	0.286	−0.110	0.079	0.058	−0.369	0.094	−0.005
平均最高气温 Mean max. temperature	同期 Same period	0.857^**	0.456^*	0.074	0.339	−0.011	0.044	0.215
	前1月 Previous one month	0.627^**	0.132	−0.285	−0.183	0.309	0.001	−0.052
	前2月 Previous two months	0.334	0.125	0.064	−0.177	−0.093	0.270	0.139
	前3月 Previous three months	0.307	−0.076	0.125	−0.110	−0.382	0.076	0.021
平均最低气温 Mean min. temperature	同期 Same period	0.633^**	0.567^**	0.453^*	0.180	0.123	−0.230	0.095
	前1月 Previous one month	0.512^*	0.033	0.113	0.294	0.177	0.473^*	0.478^*
	前2月 Previous two months	0.184	−0.007	−0.178	0.120	0.249	−0.075	0.093
	前3月 Previous three months	0.256	−0.143	0.007	−0.133	−0.217	0.180	−0.030

下载: 导出CSV

| 显示表格

表 7 东北森林逐月气象因子变化趋势

Table 7. Monthly changing trend of meteorological factors in the forest of Northeast China

气象因子 Meteorological factor	1月 Jan.	2月 Feb.	3月 Mar.	4月 Apr.	5月 May	6月 Jun.	7月 Jul.	8月 Aug.	9月 Sep.	10月 Oct.
降水量/（mm·10⁻¹ a⁻¹） Precipitation/(mm·10⁻¹ year⁻¹)	−3.85^**	1.51	−1.91	−5.90	16.68^*	18.49	1.30	29.62^*	15.62^*	−1.59
平均气温/（℃·10⁻¹ a⁻¹） Mean temperature/(℃·10⁻¹ year⁻¹)	0.61	0.12	0.71	0.31	0.20	−0.45	0.35	0.48	0.01	0.09
平均最高气温/（℃·10⁻¹ a⁻¹） Mean max. temperature/(℃·10⁻¹ year⁻¹)	0.61	0.08	0.99	0.69	0.19	−0.83	0.44	−0.33	−0.33	0.25
平均最低气温/（℃·10⁻¹ a⁻¹） Mean min. temperature/(℃·10⁻¹ year⁻¹)	0.58	0.04	0.40	−0.23	0.28	0.20	0.38	0.62	0.42	−0.05
注：负值表示降低。Note: the negative values indicate decrease.

下载: 导出CSV

| 显示表格

2.3 不足与展望

本文所做研究目前可能存在一些不足之处。首先，植被生长发育不仅仅和温度及降水有关，还受到诸如辐射、风速、云量等其他气候因子的影响，且人为活动也是影响植被生长的因素之一，而本文未能将上述影响考虑在内；其次，当前遥感卫星监测到的NDVI影像分辨率还比较低，会影响实验精度，对结果有一定影响；此外，极端天气对植物生长发育也存在一定影响，而本文受制于气象数据为月值数据集而未能深入分析。因此，未来研究可综合其他气候因子探讨其对森林植被生长的影响，利用更高精度的遥感影像进行分析，并考虑极端天气对植被生长动态的干扰作用，以期为科学营林和生态管护提供更多方面的参考。

3. 结　　论

本文利用SPOT Vegetation NDVI数据多角度分析了2000—2019年东北森林变化动态，并结合气象数据探究了气候变化条件下的响应特征，得出以下结论：

（1）近20年东北森林生长季NDVI值集中于0.40 ~ 0.78，并且整体呈现上升趋势，变化趋势为每10年增加0.035，其中混交林上升趋势最强，每10年达0.037，针叶林次之，阔叶林上升趋势最小，东北森林出现轻度改善和中度改善的面积占总面积的55.59%，基本不变的面积占43.84%，表明东北森林在过去20年主要以改善为主。

（2）在未来趋势上，东北森林生长季NDVI变化反持续性强于正持续性，70.06%的森林将面临由改善变为退化的趋势，26.39%将持续改善，预示未来东北森林呈退化态势，尤其是最低气温和降水变化较明显的中部和北部地区气候变化有可能由促进变为抑制。

（3）东北森林生长季NDVI与降水和最低气温均呈正相关关系，与最高气温关系不密切。生长季前期NDVI主要和最高气温和最低气温关系显著且具有一定程度的时滞性，生长季中期NDVI主要受前期水分条件影响；生长季后期NDVI主要受降水和最低气温影响且具有时滞性。最低气温对NDVI的影响在生长季前后期均有出现，而最高气温作用仅在生长季前期较为明显。不同林型对气象因子的响应存在差别，针叶林生长主要受温度因子影响，尤其是最低气温；阔叶林生长则主要受降水和最低气温的共同控制；混交林生长对降水最为敏感。

（4）东北森林地区在过去20年白天夜晚温度呈现不同步变化，夜晚温度上升比白天温度上升快，平均气温上升主要是由夜晚温度上升导致的。同时东北森林NDVI与最低气温的关系要显著强于最高气温。温度变化的昼夜不对称和响应的昼夜不对称是此区域NDVI上升的重要原因。此外，同期降水的上升也对NDVI的增加起到促进作用。这提示我们在未来的研究中，不应仅关注平均气温对植物生长的影响，也要考虑最高气温和最低气温影响的差异。

图 1 BpSPL6基因启动子的克隆

M. DL-5000 Marker；1. BpSPL6基因启动子扩增产物。M, DL-5000 Marker; 1, amplified production of BpSPL6 promoter.

Figure 1. Cloning of BpSPL6 promoter by PCR

下载: 全尺寸图片幻灯片

图 2 BpSPL6基因启动子核苷酸序列多重比对

1-BpSPL6-promoter. 测序结果；2-BpSPL6-promoter. 白桦基因组数据库中启动子序列；3-BpSPL6-promoter. 欧洲白桦基因组数据库中启动子序列。1-BpSPL6-promoter, sequencing results; 2-BpSPL6-promoter, BpSPL6 promoter sequences in Betula platyphylla genome database; 3-BpSPL6-promoter, BpSPL6 promoter sequences in Betula pendula genome database.

Figure 2. Multi-alignment of nucleotide sequences of BpSPL6 promoter

下载: 全尺寸图片幻灯片

图 3 白桦BpSPL6基因启动子序列及顺式作用元件分布

方框及彩色表示预测的顺势作用元件序列；方框下文字表示对应元件名称；图片下部注释了彩色序列对应的元件名称；下划线表示启动子上游未克隆序列；双下划线并加粗表示起始密码子位置。The box and color represent the predicted cis-elements; text under the box represents the corresponding cis-element name; the cis-element names corresponding to the color sequence are annotated at the bottom of the picture; the underline indicates the uncloned sequence upstream of the promoter; the double underline and bold indicate the position of start codon.

Figure 3. Promoter sequence and cis-element distribution of BpSPL6 gene promoter for B. platyphylla

下载: 全尺寸图片幻灯片

图 4 大肠杆菌重组质粒的菌液PCR检测

M. DL-5000 Marker；1. 阴性对照（水对照）；2 ~ 3. 大肠杆菌转化子。M, DL-5000 Marker; 1, negative control (water control); 2−3, Escherichia coli transformant.

Figure 4. PCR detection results of recombinant plasmid of Escherichia coli

下载: 全尺寸图片幻灯片

图 5 转基因拟南芥的检测

M. DL-5000 Marker；1. 阳性对照；2 ~ 6. 转基因株系；7. 阴性对照（野生型拟南芥）。M, DL-5000 Marker; 1, positive control; 2−6, transgenic lines; 7, negative control (wild type A. thaliana).

Figure 5. PCR detection of transgenic A. thaliana

下载: 全尺寸图片幻灯片

图 6 营养生长期的转基因拟南芥GUS染色

A ~ C. 生长7 d的拟南芥（A. 整株幼苗；B. 真叶；C. 根）。D ~ G. 生长10 d的拟南芥（D. 整株幼苗；E. 真叶；F. 下胚轴；G. 根）。A−C, A. thaliana seedling grown for 7 days (A, whole seedling; B, euphylla; C, root). D−G, A. thaliana seedling grown for 10 days (D, whole seedling; E, euphylla; F, hypocotyl; G, root).

Figure 6. Staining of GUS in transgenic A. thaliana in vegetative phase

下载: 全尺寸图片幻灯片

图 7 生殖生长期的转基因拟南芥GUS染色

A ~ G. 生长25 d的拟南芥；A. 整株幼苗；B. 成熟叶片；C. 幼嫩叶片；D. 幼嫩叶片和花；E. 花；F. 幼果；G. 根。A−G, A. thaliana seedling grown for 25 days; A, whole seedling; B, mature leaf; C, young leaf; D, young leaf and flower; E, flower; F, young fruit; G, root.

Figure 7. Staining of GUS in transgenic A. thaliana in reproductive phase

下载: 全尺寸图片幻灯片

图 8 不同非生物胁迫处理下转基因拟南芥的GUS染色

A. 整株幼苗；B. 成熟叶片；C. 幼嫩叶片；D. 根；E. 下胚轴。A, whole seedling; B, mature leaf; C, young leaf; D, root; E, hypocotyl.

Figure 8. Staining of GUS in transgenic A. thaliana under different abiotic stress treatments

下载: 全尺寸图片幻灯片

表 1 启动子顺式作用元件分析

Table 1 Analysis of promoter cis-elements

元件名称 Element name	基序序列 Motif sequence	个数 Number	生物学功能 Biological function
TATA-box	TATA	6	核心启动元件 Core actuating element
CAAT-box	CAAT	22	启动子增强区保守元件 Conservative element of promoter enhance region
GT1-motif	GGTTAAT	5	光响应元件 Element involved in light responsiveness
Box 4	ATTAAT	3	光响应元件 Element involved in light responsiveness
GA-motif	ATAGATAA	2	光响应元件 Element involved in light responsiveness
AT1-motif	AATTATTTTTTATT	1	光响应元件 Element involved in light responsiveness
AuxRR-core	GGTCCAT	2	生长素响应元件 Element involved in auxin responsiveness
TGA-element	AACGAC	1	生长素响应元件 Element involved in auxin responsiveness
CATATGGMSAUR	CATATG	4	生长素响应元件 Auxin response element
NTBBF1ARROLB	ACTTTA	4	生长素响应元件 Auxin response element
GARE2OSREP1	TAACGTA	1	赤霉素响应元件 Gibberellin response element
WRKY71OS	TGAC	12	赤霉素响应元件 Gibberellin response element
P-box	CCTTTTG	1	赤霉素响应元件 Gibberellin response element
WBOXATNPR1	TTGAC	4	水杨酸响应元件 Element involved in salicylic acid responsiveness
DPBFCOREDCDC3	ACACNNG	3	脱落酸响应元件 Element involved in abscisic acid responsiveness
ABRE	ACGTG	1	脱落酸响应元件 Element involved in abscisic acid responsiveness
MYB1AT	WAACCA	5	脱水响应元件 Dehydration response element
MYCCONSENSUSAT	CANNTG	12	脱水及寒冷响应元件 Dehydration and cold response element
ACGTATERD1	ACGT	4	脱水响应元件 Dehydration response element
MYBCORE	CNGTTR	1	脱水响应元件 Dehydration response element
WBOXNTERF3	TGACY	7	伤害响应元件 Damage response element
DOFCOREZM	AAAG	33	伤害响应元件 Damage response element
CCAATBOX1	CCAAT	4	热激信号响应元件 Heat shock signal response element
POLLEN1LELAT52	AGAAA	14	花粉特异表达的顺式作用元件 cis-acting element for pollen specific expression
ROOTMOTIFTAPOX1	ATATT	10	与根相关的顺式作用元件 cis-acting element in root
注：基序序列：N = A/T/G/C，R = A/G，W = A/T，Y = C/T。Notes: motif sequence: N = A/T/G/C, R = A/G, W = A/T, Y = C/T.

下载: 导出CSV

参考文献(39)

[1]	Pino M T, Skinner J S, Park E J, et al. Use of a stress inducible promoter to drive ectopic AtCBF expression improves potato freezing tolerance while minimizing negative effects on tuber yield[J]. Plant Biotechnology Journal, 2007, 5(5): 591−604. doi: 10.1111/j.1467-7652.2007.00269.x
[2]	Butler J E, Kadonaga J T. The RNA polymerase II core promoter: a key component in the regulation of gene expression[J]. Genes & Development, 2002, 16(20): 2583−2592.
[3]	Klein J, Saedler H, Huijser P. A new family of DNA binding proteins includes putative transcriptional regulators of the Antirrhinum majus floral meristem identity gene SQUAMOSA [J]. Molecular & General Genetics, 1996, 250(1): 7−16.
[4]	Yamasaki K, Kigawa T, Inoue M, et al. A novel zinc-binding motif revealed by solution structures of DNA-binding domains of Arabidopsis SBP-family transcription factors[J]. Journal of Molecular Biology, 2004, 337(1): 49−63. doi: 10.1016/j.jmb.2004.01.015
[5]	Usami T, Horiguchi G, Yano S, et al. The more and smaller cells mutants of Arabidopsis thaliana identify novel roles for SQUAMOSA PROMOTER BINDING PROTEIN-LIKE genes in the control of heteroblasty[J]. Development, 2009, 136(6): 955−964. doi: 10.1242/dev.028613
[6]	Cardon G H, Hohmann S, Nettesheim K, et al. Functional analysis of the Arabidopsis thaliana SBP-box gene SPL3: a novel gene involved in the floral transition[J]. Plant Journal, 1997, 12(2): 367−377. doi: 10.1046/j.1365-313X.1997.12020367.x
[7]	Wu G, Poethig R S. Temporal regulation of shoot development in Arabidopsis thaliana by miR156 and its target SPL3[J]. Development, 2006, 133(18): 3539−3547. doi: 10.1242/dev.02521
[8]	Jung J H, Lee H J, Ryu J Y, et al. SPL3/4/5 integrate developmental aging and photoperiodic signals into the FT-FD module in Arabidopsis flowering[J]. Molecular Plant, 2016, 9(12): 1647−1659. doi: 10.1016/j.molp.2016.10.014
[9]	Wang J W, Czech B, Weigel D. miR156-regulated SPL transcription factors define an endogenous flowering pathway in Arabidopsis thaliana[J]. Cell, 2009, 138(4): 738−749. doi: 10.1016/j.cell.2009.06.014
[10]	Unte U S, Sorensen A M, Pesaresi P, et al. SPL8, an SBP-box gene that affects pollen sac development in Arabidopsis[J]. Plant Cell, 2003, 15(4): 1009−1019. doi: 10.1105/tpc.010678
[11]	Preston J C, Jorgensen S A, Orozco R, et al. Paralogous SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes differentially regulate leaf initiation and reproductive phase change in petunia[J]. Planta, 2016, 243(2): 429−440. doi: 10.1007/s00425-015-2413-2
[12]	Yu N, Niu Q W, Ng K H, et al. The role of miR156/SPLs modules in Arabidopsis lateral root development[J]. Plant Journal, 2015, 83(4): 673−685. doi: 10.1111/tpj.12919
[13]	Zhang X H, Dou L H, Pang C Y, et al. Genomic organization, differential expression, and functional analysis of the SPL gene family in Gossypium hirsutum[J]. Molecular Genetics and Genomics, 2015, 290(1): 115−126. doi: 10.1007/s00438-014-0901-x
[14]	Filichkin S A, Ansariola M, Fraser V N, et al. Identification of transcription factors from NF-Y, NAC, and SPL families responding to osmotic stress in multiple tomato varieties[J]. Plant Science, 2018, 274: 441−450. doi: 10.1016/j.plantsci.2018.06.021
[15]	Stief A, Altmann S, Hoffmann K, et al. Arabidopsis miR156 regulates tolerance to recurring environmental stress through SPL transcription factors[J]. Plant Cell, 2014, 26(4): 1792−1807. doi: 10.1105/tpc.114.123851
[16]	Wang J W, Ye Y J, Xu M, et al. Roles of the SPL gene family and miR156 in the salt stress responses of tamarisk (Tamarix chinensis)[J]. BMC Plant Biology, 2019, 19(1): 370. doi: 10.1186/s12870-019-1977-6
[17]	Hou H M, Jia H, Yan Q, et al. Overexpression of a SBP-Box gene (VpSBP16) from Chinese wild Vitis species in Arabidopsis improves salinity and drought stress tolerance[J]. International Journal of Molecular Sciences, 2018, 19(4): 940.
[18]	Gou J Q, Debnath S, Sun L, et al. From model to crop: functional characterization of SPL8 in M. truncatula led to genetic improvement of biomass yield and abiotic stress tolerance in alfalfa[J]. Plant Biotechnology Journal, 2018, 16(4): 951−962. doi: 10.1111/pbi.12841
[19]	Ning K, Chen S, Huang H J, et al. Molecular characterization and expression analysis of the SPL gene family with BpSPL9 transgenic lines found to confer tolerance to abiotic stress in Betula platyphylla Suk[J]. Plant Cell, Tissue and Organ Culture (PCTOC), 2017, 130(3): 469−481. doi: 10.1007/s11240-017-1226-3
[20]	Wei H, Meilan R, Brunner A M, et al. Transgenic sterility in Populus: expression properties of the poplar PTLF, Agrobacterium NOS and two minimal 35S promoters in vegetative tissues[J]. Tree Physiology, 2006, 26(4): 401−410. doi: 10.1093/treephys/26.4.401
[21]	Cook M, Thilmony R. The OsGEX2 gene promoter confers sperm cell expression in transgenic rice[J]. Plant Molecular Biology Reporter, 2012, 30(5): 1138−1148. doi: 10.1007/s11105-012-0429-3
[22]	Wang J Y, Wang J P, Yang H F. Identification and functional characterization of the NAC gene promoter from Populus euphratica[J]. Planta, 2016, 244(2): 417−427. doi: 10.1007/s00425-016-2511-9
[23]	Chen Y, Chen S, Chen F, et al. Functional characterization of a Chrysanthemum dichrum stress-related promoter[J]. Molecular Biotechnology, 2012, 52(2): 161−169. doi: 10.1007/s12033-011-9483-6
[24]	丁雪峰, 刘鸿艳, 罗利军. 水稻OsGRAS1启动子的克隆及多样性分析[J]. 上海农业学报, 2010, 26(4): 8−14. doi: 10.3969/j.issn.1000-3924.2010.04.003 Ding X F, Liu H Y, Luo L J. Cloning and diversity analysis of the OsGRAS1 promoter in rice[J]. Acta Agriculturae Shanghai, 2010, 26(4): 8−14. doi: 10.3969/j.issn.1000-3924.2010.04.003
[25]	张勇, 胡晓晴, 李豆, 等. 白桦BpSPL8启动子的克隆及异源过表达BpSPL8对拟南芥耐旱性的影响[J]. 北京林业大学学报, 2019, 41(8): 67−75. Zhang Y, Hu X Q, Li D, et al. Cloning the promoter of BpSPL8 from Betula platyphylla and overexpression of BpSPL8 gene affecting drought tolerance in Arabidopsis thaliana[J]. Journal of Beijing Forestry University, 2019, 41(8): 67−75.
[26]	Higo K, Ugawa Y, Iwamoto M, et al. Plant cis-acting regulatory DNA elements (PLACE) database: 1999[J]. Nucleic Acids Research, 1999, 27(1): 297−300. doi: 10.1093/nar/27.1.297
[27]	Lescot M, Dehais P, Thijs G, et al. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences[J]. Nucleic Acids Research, 2002, 30(1): 325−327. doi: 10.1093/nar/30.1.325
[28]	Clough S J, Bent A F. Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana[J]. Plant Journal, 1998, 16(6): 735−743. doi: 10.1046/j.1365-313x.1998.00343.x
[29]	Chen X, Zhang Z, Liu D, et al. SQUAMOSA promoter-binding protein-like transcription factors: star players for plant growth and development[J]. Journal of Integrative Plant Biology, 2010, 52(11): 946−951. doi: 10.1111/j.1744-7909.2010.00987.x
[30]	Shikata M, Koyama T, Mitsuda N, et al. Arabidopsis SBP-box genes SPL10, SPL11 and SPL2 control morphological change in association with shoot maturation in the reproductive phase[J]. Plant & Cell Physiology, 2009, 50(12): 2133−2145.
[31]	Arshad M, Feyissa B A, Amyot L, et al. MicroRNA156 improves drought stress tolerance in alfalfa (Medicago sativa) by silencing SPL13[J]. Plant Science, 2017, 258: 122−136.
[32]	Xu M, Hu T, Zhao J, et al. Developmental functions of miR156-regulated SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes in Arabidopsis thaliana[J/OL]. PLoS Genetics, 2016, 12(8): e1006263 [2020−01−23]. https://doi.org/10.1371/journal.pgen.1006263.
[33]	Barrera-Rojas C H, Rocha G H B, Polverari L, et al. miR156-targeted SPL10 controls Arabidopsis root meristem activity and root-derived de novo shoot regeneration via cytokinin responses[J]. Journal of Experimental Botany, 2020, 71(3): 934−950. doi: 10.1093/jxb/erz475
[34]	Ye B B, Shang G D, Pan Y, et al. AP2/ERF transcription factors integrate age and wound signals for root regeneration[J]. Plant Cell, 2020, 32(1): 226−241. doi: 10.1105/tpc.19.00378
[35]	Zhang Y, Schwarz S, Saedler H, et al. SPL8, a local regulator in a subset of gibberellin-mediated developmental processes in Arabidopsis[J]. Plant Molecular Biology, 2007, 63(3): 429−439. doi: 10.1007/s11103-006-9099-6
[36]	Galvao V C, Horrer D, Kuttner F, et al. Spatial control of flowering by DELLA proteins in Arabidopsis thaliana[J]. Development, 2012, 139(21): 4072−4082. doi: 10.1242/dev.080879
[37]	Li J, Hou H M, Li X Q, et al. Genome-wide identification and analysis of the SBP-box family genes in apple (Malus × domestica Borkh.)[J]. Plant Physiology and Biochemistry, 2013, 70: 100−114.
[38]	Mao H D, Yu L J, Li Z J, et al. Genome-wide analysis of the SPL family transcription factors and their responses to abiotic stresses in maize[J]. Plant Gene, 2016, 6: 1−12.
[39]	Cui L G, Shan J X, Shi M, et al. The miR156-SPL9-DFR pathway coordinates the relationship between development and abiotic stress tolerance in plants[J]. Plant Journal, 2014, 80(6): 1108−1117. doi: 10.1111/tpj.12712

施引文献(24)

期刊类型引用(16)

1.	康向阳. 关于我国林木育种向智能分子设计育种发展的思考. 北京林业大学学报. 2024(03): 1-7 . 本站查看
2.	Junfei Hao，Na Chen，Pingyu Yan，Kaiyuan Xu，Lei Zhang，Hanguo Zhang. Study on the variation in and selection of Fraxinus mandshurica provenances and families in northeast China. Journal of Forestry Research. 2023(02): 519-529 . 必应学术
3.	朱艳，王久亮，王芳，安文娜，糜加轩，万雪琴，杨汉波. 泡核桃无性系油脂脂肪酸组成和营养成分的比较分析. 中国粮油学报. 2023(06): 84-91 . 百度学术
4.	杨涛，邱勇斌，沈汉，郑成忠，张振，王文月，金国庆，周志春. 柏木无性系和家系含碳量的早期评价与优良品系选择. 林业科学. 2023(09): 85-94 . 百度学术
5.	杜康，王加焕，李俊恒，李超，张平冬，康向阳. 毛白杨耐寒种质资源遗传鉴定及评价. 北京林业大学学报. 2023(12): 59-67 . 本站查看
6.	张超，张建林，张淑琴. 林木育种方法在园林植物培育中的应用. 分子植物育种. 2022(02): 499-502 . 百度学术
7.	闫卓，秦金凤，瞿辉，俞菲. 砧木和接穗对嫁接月季生长的影响. 江苏农业科学. 2022(04): 126-130 . 百度学术
8.	史月冬，郑宏，叶代全，施季森，边黎明. 杉木生长性状的空间与竞争效应及其对遗传参数估计的影响. 林业科学. 2022(05): 75-84 . 百度学术
9.	刘春鹏，张晓华，徐振华，杜克久，李计达，李向军，李新利. 毛白杨雄花发育特征及少粉型种质选择与成因分析. 浙江林业科技. 2021(04): 10-14 . 百度学术
10.	高本旺，欧阳芳群，高晗，李薇，雷华，田开春，祁松，王军辉. 鄂西地区欧洲云杉幼龄无性系生长差异及早期评价与选择. 林业科学研究. 2021(05): 88-94 . 百度学术
11.	杨晓伟，朱宁华，韩志强，刘天宇，李恒，杨捷，朱新传. 铁心杉半同胞子代优株父本鉴定. 中南林业科技大学学报. 2021(12): 26-35 . 百度学术
12.	康向阳. 林木三倍体育种研究进展及展望. 中国科学:生命科学. 2020(02): 136-143 . 百度学术
13.	康向阳. 林木遗传育种研究进展. 南京林业大学学报(自然科学版). 2020(03): 1-10 . 百度学术
14.	何丽玲. 杉木优良无性系组培繁育技术的研究. 农业技术与装备. 2020(07): 114-115 . 百度学术
15.	康向阳. 关于林木育种策略的思考. 北京林业大学学报. 2019(12): 15-22 . 本站查看
16.	吕运舟，蒋泽平，梁珍海，董筱昀，杨勇，黄利斌，施士争. 榔榆无性系苗期测定与评价. 江苏林业科技. 2019(06): 13-16 . 百度学术

其他类型引用(8)

资源附件(0)

图(8) / 表(1)

计量

文章访问数: 2350
HTML全文浏览量: 827
PDF下载量: 141
被引次数: 24

1. 研究区概况与研究方法
1.1 研究区概况
1.2 数据来源与处理
1.3 研究方法
1.3.1 趋势分析
1.3.2 Hurst指数
1.3.3 相关分析
2. 结果与分析
2.1 东北森林生长季NDVI时空变化特征
2.1.1 NDVI空间分布特征
2.1.2 NDVI变化趋势
2.1.3 NDVI 未来趋势分析
2.2 东北森林生长季NDVI对气象因子的响应
2.2.1 东北森林生长季气象因子的变化特征
2.2.2 生长季NDVI与气候变化的关系
2.2.3 逐月NDVI与气象因子的相关性与时滞性分析
2.3 不足与展望
3. 结　　论

1. 研究区概况与研究方法
1.1 研究区概况
1.2 数据来源与处理
1.3 研究方法
1.3.1 趋势分析
1.3.2 Hurst指数
1.3.3 相关分析
2. 结果与分析
2.1 东北森林生长季NDVI时空变化特征
2.1.1 NDVI空间分布特征
2.1.2 NDVI变化趋势
2.1.3 NDVI 未来趋势分析
2.2 东北森林生长季NDVI对气象因子的响应
2.2.1 东北森林生长季气象因子的变化特征
2.2.2 生长季NDVI与气候变化的关系
2.2.3 逐月NDVI与气象因子的相关性与时滞性分析
2.3 不足与展望
3. 结　　论

参考文献(39)

施引文献(24)

资源附件(0)

白桦BpSPL6基因启动子的克隆及表达分析

作者简介: 李豆。主要研究方向：林木遗传育种。Email：1171057158@qq.com 地址：150040 黑龙江省哈尔滨市香坊区和兴路26号东北林业大学生命科学学院

责任作者: 王宏伟，博士，副教授。主要研究方向：生物活性物质制备与应用。Email：84970486@qq.com 地址：同上 刘雪梅，博士，教授。主要研究方向：林木遗传育种。Email：695898040@qq.com 地址：同上

计量

出版历程