• Scopus收录期刊
  • CSCD（核心库）来源期刊
  • 中文核心期刊
  • 中国科技核心期刊
  • F5000顶尖学术来源期刊
  • RCCSE中国核心学术期刊
高级检索

基于多层感知机的长白落叶松人工林林分生物量模型

徐奇刚, 雷相东, 国红, 李海奎, 李玉堂

徐奇刚, 雷相东, 国红, 李海奎, 李玉堂. 基于多层感知机的长白落叶松人工林林分生物量模型[J]. 北京林业大学学报, 2019, 41(5): 97-107. DOI: 10.13332/j.1000-1522.20190035
引用本文: 徐奇刚, 雷相东, 国红, 李海奎, 李玉堂. 基于多层感知机的长白落叶松人工林林分生物量模型[J]. 北京林业大学学报, 2019, 41(5): 97-107. DOI: 10.13332/j.1000-1522.20190035
shu
Xu Qigang, Lei Xiangdong, Guo Hong, Li Haikui, Li Yutang. Stand biomass model of Larix olgensis plantations based on multi-layer perceptron networks[J]. Journal of Beijing Forestry University, 2019, 41(5): 97-107. DOI: 10.13332/j.1000-1522.20190035
Citation: Xu Qigang, Lei Xiangdong, Guo Hong, Li Haikui, Li Yutang. Stand biomass model of Larix olgensis plantations based on multi-layer perceptron networks[J]. Journal of Beijing Forestry University, 2019, 41(5): 97-107. DOI: 10.13332/j.1000-1522.20190035
shu

基于多层感知机的长白落叶松人工林林分生物量模型

  • 1.

    中国林业科学研究院资源信息研究所，北京 100091

  • 2.

    国家林业和草原局森林经营与生长模拟实验室，北京 100091

  • 3.

    吉林省林业调查规划院，吉林 长春 130022

基金项目: 林业行业公益性科研项目“我国主要林区林地立地质量和生产力评价研究”（201504303）
详细信息
    作者简介:

    徐奇刚，博士生。主要研究方向：森林生长收获模型。Email：adslxqg@126.com　地址：100091 北京市海淀区香山路东小府1号中国林业科学研究院资源信息研究所

    责任作者:

    雷相东，研究员，博士生导师。主要研究方向：森林生长模型与模拟。Email：xdlei@caf.ac.cn　地址：同上

Stand biomass model of Larix olgensis plantations based on multi-layer perceptron networks

  • 1.

    Institute of Forest Resource Information Techniques, Chinese Academy of Forestry, Beijing 100091, China

  • 2.

    Key Labotary of Forest Management and Growth Modelling, National Forestry and Grassland Administration, Beijing 100091, China

  • 3.

    Jilin Forestry Inventory and Planning Institute, Changchun 130022, Jilin, China

HTML全文

    摘要
  • 摘要:
    目的神经网络模型能避免林分生物量模型建模时自变量共线性与异方差问题，研究多层感知机在林分生物量模型中的应用，为森林经营单位、区域生物量和碳储量的估算提供方法和依据。
    方法以长白落叶松人工林为研究对象，利用吉林省一类清查固定样地的917组数据，分别建立了基于传统的对数转化后线性模型和神经网络多层感知机的地上生物量和总生物量模型。使用AIC、决定系数（R2）、均方根误差（RMSE）、相对均方根误差（RMSEr）和平均绝对误差（MAE）来评价模型。
    结果估计精度最高的模型是输入单元为林分平均胸径（D）−平均高（H）−林分密度指数（S）−海拔（HB）−坡向（PX）−坡位（PW）、2个隐藏层、隐单元数为40−20的神经网络模型，与传统对数转换线性回归模型相比，地上生物量和总生物量模型的调整决定系数（Adj.R2）分别从0.902 1提高到了0.914 1，从0.897 9提高到了0.908 9；RMSEr分别从6.330 5%降低到了5.992 2%，从6.490 1%降低到了6.153 6%。包含立地因子的神经网络模型比未包含立地因子的神经网络模型估计精度略有提升，地上生物量与总生物量的Adj.R2分别提高了0.88%和0.99%，RMSEr分别降低了5.33%和5.46%。
    结论多层感知机生物量模型的估计精度比传统回归模型略有提高，但它可以避免模型选型和违背传统统计假设的处理等问题，且能够一次性计算地上生物量和总生物量模型，有一定优势。
    Abstract:
    ObjectiveNeural network model can avoid the collinearity and heteroscedasticity of variables in modeling forest stand biomass. This paper aims to apply multi-layer perceptron networks to forest biomass model to provide methods for the calculation and prediction of forest biomass and carbon stocks at forest management unit and regional levels.
    MethodBased on 917 observations from the sample plots of larch plantations from national forest inventory in Jilin Province of northeastern China, the aboveground and total biomass models by log-transformed linear regression, and multi-layer perceptron networks with and without site factors were established. AIC, adjusted R2, RMSE, RMSEr and MAE were used to evaluate the models.
    ResultThe model with the highest prediction accuracy was the neural network one with the input unit quadratic mean diameters (D), stand mean height (H), stand density index (S), altitude (HB), slope (PD), slope aspect (PX), slope position (PW), two hidden layers and hidden unit number of 40−20. Compared with traditional log-transformed linear regression model, the adjusted R2 of the aboveground and total biomass models was increased from 0.902 1 to 0.914 1, and from 0.897 9 to 0.908 9, RMSEr was decreased from 6.330 5% to 5.992 2%, and from 6.490 1% to 6.153 6%, respectively. The neural network model with site factors had slightly higher estimation accuracy than that without site factors. The adjusted R2 of the aboveground and total biomass models was increased by 0.88% and 0.99%, and RMSEr decreased by 5.33% and 5.46%, respectively.
    ConclusionBiomass model based on multi-layer perceptron networks had similar performance in terms of model accuracy, but it could avoid treating the traditional assumptions such as the collinearity and heteroscedasticity of variables, and had the ability to calculate aboveground and total biomass models at one time.
    • HTML全文

  • 随着全球气候变暖，极端气候频现，干旱、盐害等非生物胁迫已成为威胁全球农林业生产和发展的一个重大挑战[1]。植物为了适应逆境环境并得以生存，进化出相关基因来维持细胞内结构和功能机制的稳定[2]。迄今为止，已发现许多与干旱和盐胁迫防御相关的功能基因，其中包括膜联蛋白基因。膜联蛋白是一个超级基因家族的多功能脂结合蛋白，到目前已从超过65个物种中分离和克隆了编码基因，包括真菌、原生生物、植物、高等脊椎动物和原核生物 [3]。Creutz等人[4]1978年首次从牛肾上腺髓质嗜铬细胞中分离纯化出一种蛋白质Annexin A7（synexin），这种蛋白促使嗜铬颗粒通过依赖Ca2+来结合一些膜或者内膜成分。而国内学者对植物膜联蛋白的研究较晚，第1个植物膜联蛋白是从悬浮培养的番茄 (Solanum lycopersicum) 细胞中发现和分离的[5]，且与动物膜联蛋白在序列上有较高的同源性[6]。随着对植物膜联蛋白的进一步了解，发现膜联蛋白广泛分布于各种植物中，如棉花（Gossypium spp.）[7]、烟草（Nicotiana tabacum）[8]、玉米（Zea mays）[9]、拟南芥（Arabidopsis thaliana）[10]、水稻（Oryza sativa）[11]等。不同植物中的膜联蛋白以基因家族的形式存在，其家族成员各自具有不同的功能和时空表达模式，表达具有组织特异性。据报道膜联蛋白参与重要的生物过程，如膜运输、细胞骨架组织、细胞内稳态和离子运输，具有Ca2+结合、离子通道和过氧化物酶活性等不同的生物学功能，在植物生长发育、信号转导、适应环境和响应非生物逆境胁迫中发挥着重要作用[12]。在渗透胁迫、脱落酸（ABA）或干旱等非生物胁迫时紫花苜蓿（Medicago sativa）中的膜联蛋白被激活上调表达[13]。拟南芥中膜联蛋白的表达能响应各种非生物胁迫（如盐胁迫、渗透胁迫和低温胁迫）而发生变化[14]。Cantero等人[15]评估了拟南芥中的8种膜联蛋白基因表达模式，并建立了它们在盐、干旱和其他非生物胁迫下的差异调控机制。拟南芥t-dna突变体表明，AtAnn1和AtAnn4在萌发过程中能应答渗透胁迫[16]。拟南芥AtAnn1基因是植物中研究较多的一个膜联蛋白，过表达AtAnn1导致耐旱性增强，而AtAnn1缺失突变株抗旱能力降低，对干旱更敏感[17]。印度芥菜（Brassica juncea）中，过表达AnnBj1分别增强了转基因烟草和棉花的耐旱性和耐盐性[18]。小麦 经低温胁迫可诱导膜联蛋白p39和p22.5的表达，嵌插在质膜上感受或转导钙信号的同时还调节胞质中的钙离子浓度适应胁迫环境[19]。虽然膜联蛋白在植物胁迫反应中发挥重要作用，但不同植物膜联蛋白的分子性能、生理和生物学功能是否有所不同，还属未知。

    胡杨（Populus euphratica）适应非生物胁迫环境能力极强，其抗逆机制受到国内外广泛关注。胡杨耐旱性的生理及分子机制虽有较多的研究，但目前未有胡杨膜联蛋白与水分胁迫相关性的报道。本研究通过对比胡杨Annexin1转基因拟南芥、野生型、突变体（atann1）在渗透胁迫、干旱胁迫及复水后的表型及生理生化反应，明确PeAnn1在植物耐受水分胁迫中的作用，研究结果有助于揭示树木的抗旱性机制。

    1.   材料与方法

    1.1   实验材料

    将来自新疆的胡杨1年生实生苗栽植在北京林业大学苗圃的温室，培养3个月后，利用250 mmol/L甘露醇处理12 h，在渗透胁迫处理0、4、8、12 h后采集叶片，进行RT-qPCR分析。

    将野生型拟南芥（WT）、突变体（atann1）和2个过表达胡杨PeAnn1拟南芥（PeAnn1-OE1和PeAnn1-OE2）种子在超净台中用质量分数为1%的次氯酸钠消毒10 min后，无菌水冲洗5次，播种于1/2 MS固体培养基（含有10 g/L蔗糖和质量分数为0.3%的植物凝胶，pH 5.8）上，4 ℃低温暗培养春化2 d，取出置于22 ℃、16 h光照/8 h黑暗的人工气候培养箱中生长。生长7 d的幼苗移栽于草炭土和蛭石的混合比例为 2∶1(w∶w)的土壤中，转入22 ℃、16 h光照/8 h黑暗，55% ~ 65%相对湿度的温室中生长。

    胡杨苗木培养，PeAnn1基因克隆，RT-qPCR分析，以及转基因植物的培育及分子检测参见文献[20]。RNA提取试剂盒购自康为世纪生物科技有限公司；cDNA反转录试剂盒购自普洛麦格生物技术有限公司；其余试剂由北京拜尔迪生物科技有限公司提供。引物合成由北京睿博兴科生物技术有限公司完成。RNA提取、纯化、RT-qPCR参见文献[21]。

    1.2   PeAnn1序列与进化树分析

    通过NCBI（http://blast.ncbi.nlm.nih.gov/Blast.cgi）在线获取不同植物Annexin同源基因序列。利用Mega 6软件（http://www.megasoftware.net/index.php）对胡杨PeAnn1、毛白杨（Populus tomentosa)PtAnn1、拟南芥AtAnn1、野生大豆（Glycine soja）GsANN和水稻OsAnn8等基因序列进行多重序列比对分析，并对不同物种同源Annexins构建系统进化树。

    1.3   渗透胁迫下拟南芥耐旱性检测

    1.3.1   种子萌发率

    在1/2 MS固体培养基中添加不同浓度的甘露醇模拟渗透胁迫，甘露醇浓度梯度设为150、200、250、300 mmol/L 4个梯度，以不加甘露醇的培养基作为对照。拟南芥种子在超净台中用质量分数为1%的次氯酸钠消毒10 min后，无菌水冲洗5次，依次点播在各个培养基上，每个株系点播50粒种子，每种处理做3次重复。4 ℃低温暗培养2 d后置于22 ℃、16 h光照/8 h黑暗的恒温培养箱中培养，10 d后记录萌发率。

    1.3.2   根　长

    选取在1/2MS固体培养基中萌发5 d，长势相似的拟南芥幼苗，移栽到正常对照培养基和甘露醇浓度梯度150、200、250、300 mmol/L的固体培养基上，观察拍照并测量记录拟南芥竖直生长7 d后的根长。

    1.4   渗透胁迫下转基因拟南芥H2O2含量的测定

    采用荧光探针H2DCF-DA（molecular probe）染色测定H2O2含量[22]。将在1/2MS固体培养基上萌发、生长7 d的各株系拟南芥幼苗（WT、atann1、PeAnn1-OE1和PeAnn1-OE2）移到含有250 mmol/L甘露醇的LMS培养基中，12 h后取出幼苗，将幼苗浸入10 μmol/L H2DCF-DA的探针溶液，室温避光孵育15 min。取出幼苗，用蒸馏水冲洗3 ~ 4次，去除表面残留探针染料。利用Leica SP5激光共聚焦显微镜（Wetzlar, 德国）测定H2DCF-DA绿色荧光，激发光波长488 nm，发射光波长510 ~ 530 nm，相对荧光强度由Image-Pro Plus version 4.5软件计算。

    1.5   渗透胁迫下转基因拟南芥抗氧化酶活性的测定

    采用紫外吸收法[23]测定过氧化氢酶（CAT）活性，以1 min内ΔA240减少0.01的酶量为1个过氧化氢酶活性单位（U）。采用NBT（氮蓝四唑）光化还原法[24]测定超氧化物歧化酶（SOD）活性，以抑制NBT光还原50%的酶量为1个酶活力单位（U）。采用愈创木酚法[25]测定过氧化物酶（POD）活性，以每分钟OD值变化（升高）0.01为1个酶活性单位（U）。按照考马斯亮蓝G250法进行蛋白含量的测定。

    1.6   转基因拟南芥中抗逆相关基因表达水平的检测

    将野生型拟南芥和转基因拟南芥播种于1/2MS培养基上进行培养，萌发5 d后长势相似的拟南芥移至250 mmol/L甘露醇的1/2MS培养基上进行处理。处理10 d后取样，利用TRIzon试剂（康为世纪）提取总RNA并反转录合成第1条cDNA，通过荧光定量PCR检测抗逆相关基因SOD（超氧化物歧化酶基因）、CAT（过氧化氢酶基因）、POD（过氧化物酶基因）的转录水平表达情况，以拟南芥AtACTIN2（At3g18780）为内参，各抗氧化酶基因荧光定量引物见表1。

    表  1  本文实验中所用到的引物序列
    Table  1.  Gene-specific primer sequences used in this study
    引物名称 Primer name上游引物 Forward primer (5′−3′)下游引物 Reverse primer (5′−3′)
    AtACTIN2GGTAACATTGTGCTCAGTGGTGGAACGACCTTAATCTTCATGCTGC
    AtSODAGGAAACATCACTGTTGGAGATGAGTTTGGTCCAGTAAGAGGAA
    AtCATAGGATCAAACTTTGAGGGGTAGCTTGTGGTTCCTGGAATCTACT
    AtPODCGTGCCCTTCATATTGTTGGGACGCCATCAACAACGAGTC
    下载: 导出CSV 
    | 显示表格

    1.7   自然干旱条件下拟南芥叶绿素荧光参数的测定

    选取生长20 d，状况良好且一致的野生型、突变体和转基因株系通过控制浇水模拟土壤干旱，浇足水后当天作为干旱的第1天，设置5个重复。分别在正常浇水、干旱处理8 d和复水3 d后采样进行相关指标的测定。

    1.7.1   叶绿素荧光参数

    每次早上10：00开始测定，选取暗适应30 min后的拟南芥叶片。利用Junior PAM便携式荧光仪（WALZ，德国）测定叶片正常浇水，干旱处理8 d和复水3 d后的荧光参数，包括PSⅡ最大光量子效率（Fv/Fm）、实际光合量子产量（ΦPSⅡ）和相对电子传递速率（ETR）。各个参数直接从仪器导出，每个株系重复3 ~ 5次，取平均值。

    1.7.2   叶绿素相对含量

    利用日本Minolta公司生产的便携式叶绿素计（SPAD-502 Plus）原位测定叶绿素相对含量，以SPAD值表示[26]。测定时选取长势一致的叶片、叶缘和叶脉中间的部位。要注意避开叶脉集中的部位，同时适当遮挡直射阳光，以保证测量的准确性。每个株系重复3 ~ 5次，取平均值。

    1.8   数据分析

    实验数据采用Excel和SPSS 19.0软件进行数据统计、分析并绘图。用SPSS 19.0软件中的最小显著性差异法（LSD）进行单因素方差分析，用Duncan法进行多重比较分析，显著性水平设定为α = 0.05。

    2.   结果与分析

    2.1   渗透胁迫诱导胡杨PeAnn1基因表达

    胡杨苗木经250 mmol/L甘露醇处理后，PeAnn1基因表达发生明显变化，PeAnn1基因相对表达量在渗透胁迫处理4 h后显著提高，之后下降，在渗透胁迫处理12 h后降至对照水平（图1）。

    图  1  甘露醇处理对胡杨PeAnn1基因表达的影响
    胡杨1年生苗木经250 mmol/L甘露醇处理12 h，渗透胁迫处理0、4、8、12 h后采集叶片进行RT-qPCR分析，每个数值均为3次生物样本重复，*表示差异显著，P < 0.05。One-year-old seedlings of P. euphratica were subjected to 250 mmol/L mannitol and stressed for 12 hours. Then leaves were sampled after being treated for 0, 4, 8, 12 hours under osmotic stress and used for RT-qPCR analysis. Each column is the mean of three biological repeats. Asterisk represents significant differences at P < 0.05 level.
    Figure  1.  Effects of mannitol stress on expression of PeAnn1 in Populus euphratica leaves
    下载: 全尺寸图片 幻灯片

    2.2   胡杨PeAnn1序列与进化树分析

    胡杨PeAnn1与毛白杨PtAnn1、拟南芥AtAnn1、大豆GsANN和水稻OsAnn8同源基因的相似度分别为85%、72.56%、56%、29.51%，胡杨PeAnn1与毛白杨PtAnn1的相似性最高（图2A）。系统进化树分析结果表明，PeAnn1与毛白杨PtAnn1的亲缘关系最为接近（图2B）。

    图  2  胡杨PeAnn1基因序列比对与系统进化树分析
    A. 胡杨 PeAnn1基因与其他植物Annexin同源基因多重序列比对。B. 不同植物Annexin蛋白的系统进化树分析。胡杨PeAnn1（XM_011018027.1）、毛白杨PtAnn1（JX986594.1）、拟南芥AtAnn1（NM_103274.4）、野生大豆GsANN（GU474544.1）、水稻OsAnn8（LOC_Os09g20330.1）。A, multiple alignment of the deduced gene sequences of PeAnn1 with other Annexins from different plant species. B, phylogenetic tree analysis of Annexin proteins. P. euphratica Anneixn1（XM_011018027.1）, P. tomentosa PtAnn1（JX986594.1）, Arabidopsis thaliana AtAnn1（NM_103274.4）, Glycine soja GsANN（GU474544.1）, and Oryza sative OsAnn8（LOC_Os09g20330.1）.
    Figure  2.  Multiple alignment of PeAnn1 gene sequences with other Annexins from different plant species and phylogenetic tree analysis of Annexin proteins
    下载: 全尺寸图片 幻灯片

    2.3   渗透胁迫下拟南芥的萌发率

    将各株系拟南芥种子播种于含有不同浓度甘露醇的1/2 MS培养基中培养，观察萌发率和生长状态。从图3A中可以看出，未经甘露醇处理时，各个基因型的拟南芥均正常生长，长势一致，无显著差异。但在含有甘露醇的培养基上，各株系的拟南芥的生长均受到抑制，且随着甘露醇浓度的升高，生长抑制作用逐渐加强。当甘露醇浓度达到250 mmol/L时，野生型、突变体种子萌发率急剧降低，降至75.8%、78.7%，转基因拟南芥OE1和OE2的种子受到的抑制作用更严重，萌发率降至69.0%、71.6%（图3B）。当甘露醇浓度达到300 mmol/L时，转基因株系OE1和OE2的种子萌发率分别为46.9%、48.9%，仅是对照组的47.8%、49.1%，与野生型和突变体间差异显著（P < 0.05, 图3B）。结果表明，过表达PeAnn1拟南芥株系OE1和OE2在甘露醇培养基上的长势弱于野生型和突变体，对渗透胁迫更敏感。

    图  3  不同浓度甘露醇对野生型（WT）、Annexin1突变体（atann1）和过表达PeAnn1（PeAnn1-OE1和PeAnn1-OE2）拟南芥萌发率的影响
    A. 甘露醇对种子萌发生长的影响（拟南芥播种在含有不同浓度甘露醇的1/2MS培养基上，生长10 d后的照片）；B. 种子萌发率分析。不同字母表示在P < 0.05水平上差异显著，下同。A, effects of mannitol stress on seed germination rate (seeds of Arabidopsis thaliana are allowed to germinate on 1/2 MS medium supplemented with different concentrations of mannitol. Representative photographs are taken after 10 days of treatment); B, analysis on seed germination rate. Different letters denote significant differences at P < 0.05 level, the same below.
    Figure  3.  Effects of different concentrations mannitol on seed germination rate of WT, atann1 and PeAnn1-transgenic Arabidopsis thaliana (PeAnn1-OE1, PeAnn1-OE2)
    下载: 全尺寸图片 幻灯片

    2.4   渗透胁迫下拟南根长的变化

    根长结果显示，在无甘露醇的对照培养基上，转基因拟南芥OE2根长生长低于其他株系（图4A、4B）。在含有不同浓度梯度（150、200、250、300 mmol/L）的甘露醇培养基上，各拟南芥的根长生长均受到抑制，且随着处理浓度的增加抑制作用逐渐加强，其长度逐渐缩短（图4A）。250 mmol/L浓度的甘露醇处理后，野生型和突变体株系根长分别降低了42.5%和46.3%，转基因株系OE1和OE2根的长度下降更显著，达到58.8%和57.7%（图4B）。当甘露醇的浓度达到300 mmol/L时，转基因OE1和OE2的根长分别降低了72.4%和70.7%（图4B），说明高浓度的甘露醇对拟南芥根生长产生极大的抑制作用。而渗透胁迫下野生型与突变体根长显著高于转基因拟南芥OE1和OE2，说明过表达PeAnn1影响了转基因株系的根长生长。根据甘露醇梯度实验结果，确定250 mmol/L甘露醇为拟南芥渗透胁迫处理的最适浓度，并进行H2O2含量、抗氧化酶活性和基因表达等指标的测定。

    图  4  不同浓度甘露醇对野生型（WT）、Annexin1突变体（atann1）和过表达PeAnn1（PeAnn1-OE1和PeAnn1-OE2）拟南芥根长生长的影响
    A. 甘露醇对根长生长的影响（拟南芥播种在1/2MS培养基萌发5 d后，移植到含有不同浓度甘露醇的1/2MS培养基垂直生长7 d后的照片）, 标尺为1 cm；B. 根系生长分析。A, effects of mannitol stress on root growth (seeds of Arabidopsis thaliana are allowed to germinate on 1/2 MS medium for 5 days, then transplanted to 1/2 MS medium supplemented with different concentrations of mannitol. Representative photographs are taken after 7 days of treatment). Scale bar = 1 cm; B, root growth analysis.
    Figure  4.  Effects of different concentrations mannitol on root growth of WT, atann1 and PeAnn1-transgenicArabidopsis thaliana (PeAnn1-OE1, PeAnn1-OE2)
    下载: 全尺寸图片 幻灯片

    2.5   渗透胁迫下拟南芥H2O2含量的变化

    H2DCFDA是一种脂溶性荧光探针，具有细胞渗透性且自身无荧光的氧化反应指示剂。进入细胞后被细胞酯酶水解成非脂溶性的DCFH且不被质膜排出，之后在过氧化物酶的催化下可被H2O2氧化转变成为强荧光性的DCF，细胞内DCF的荧光强度与胞内H2O2水平呈显著正相关，因此被广泛用来检测细胞内活性氧（ROS）的含量和监测细胞氧化还原反应过程[27]。H2DCFDA的测定结果显示，在正常条件下，转基因株系OE1和OE2根细胞H2O2水平较高，与野生型、突变体拟南芥有显著差异。经250 mmol/L甘露醇渗透胁迫后，除转基因株系OE2外，各拟南芥株系的H2O2水平均有显著增加（图5），说明渗透胁迫下，植物体内积累了大量的活性氧，而且250 mmol/L甘露醇对转基因株系的作用更加明显：OE1和OE2 H2O2的水平分别是野生型拟南芥、突变体的1.19 ~ 1.36倍和1.37 ~ 1.55倍（图5），说明转基因拟南芥在渗透逆境下难以清除体内过多的H2O2，抗氧化防御系统能力弱于野生型拟南芥和突变体。

    图  5  甘露醇胁迫下野生型（WT）、Annexin1突变体（atann1）和过表达PeAnn1（PeAnn1-OE1和PeAnn1-OE2）拟南芥根细胞H2O2荧光强度的变化
    甘露醇（250 mmol/L）胁迫下拟南芥根细胞H2O2荧光强度的变化。H2O2 fluorescence intensity variation of Arabidopsis thaliana under mannitol stress (250 mmol/L).
    Figure  5.  Effects of mannitol stress on H2O2 fluorescence intensity in root cells of wild-type (WT), atann1 and transgenic lines(PeAnn1-OE1, PeAnn1-OE2)
    下载: 全尺寸图片 幻灯片

    2.6   渗透胁迫下拟南芥抗氧化物酶活性的变化

    当植物处于不利于自身生长的条件时，活性氧增多，SOD、CAT、POD等抗氧化酶清除超氧阴离子和H2O2，从而保持活性氧的动态平衡，减轻ROS对植物的危害[28]。通过对250 mmol/L甘露醇胁迫处理后各株系中SOD活性的检测，发现甘露醇处理后各拟南芥体内的SOD活性均显著增加，野生型和突变体的SOD值变化幅度较大，分别较对照组上升了33.1%和52.8%，说明拟南芥幼苗的SOD对其在干旱胁迫下自由基的清除具有重要作用。转基因株系OE1和OE2的SOD活性增幅较低，且在渗透胁迫后，SOD活性显著低于野生型和突变体（图6A）。POD能清除对植物体有毒害作用的氧自由基，防止细胞衰老，提高植物的抗逆能力。由图6B可知，渗透胁迫后，各拟南芥的POD活性均呈不同程度的下降，其中野生型、转基因株系OE1和OE2与对照组间差异显著，分别下降了15.2%、33.4%、35.7%；突变体拟南芥无显著差异。CAT作为植物体内分解、清除H2O2的重要保护酶类之一。在250 mmol/L甘露醇胁迫条件下，水分的缺失诱导植物体内CAT活性提高，从而加速H2O2的分解，避免羟基、自由基的产生，对植物细胞膜的过氧化有一定减缓作用。图6C中结果显示，正常条件下各个拟南芥株系的CAT活性差异显著。胁迫处理后，野生型和突变体的CAT活性显著增加，而转基因株系OE1的CAT活性升高幅度较小，OE2的CAT活性显著降低，较处理前下降了29.3%。

    图  6  甘露醇胁迫下野生型（WT）、Annexin1突变体（atann1）和过表达PeAnn1（PeAnn1-OE1和PeAnn1-OE2）拟南芥中SOD酶活性、POD酶活性和CAT酶活性的变化
    A. 甘露醇胁迫下SOD活性；B. 甘露醇胁迫下POD活性；C. 甘露醇胁迫下CAT活性。A, SOD activity under mannitol stress; B, POD activity under mannitol stress; C, CAT activity under mannitol stress.
    Figure  6.  Effects of mannitol stress on the activities of SOD, POD and CAT in wild-type (WT), atann1 and transgenic lines (PeAnn1-OE1, PeAnn1-OE2)
    下载: 全尺寸图片 幻灯片

    2.7   渗透胁迫下抗氧化酶基因表达量的变化

    为了探究PeAnn1转基因株系在渗透胁迫下抗氧化酶活性下降的原因，本文检测了渗透胁迫后各基因型抗氧化酶基因的表达。图7A结果表明，渗透胁迫处理后，AtSOD基因的表达水平在野生型和突变体拟南芥中显著升高，与对照组相比分别上升了50%、38.5%；而在转基因株系OE1和OE2中显著下降。图7B中显示，渗透胁迫处理的野生型和突变体拟南芥AtPOD基因表达水平几乎没有变化，而转基因株系OE1和OE2与对照组间差异显著，AtPOD表达量分别下降了33.4%、35.7%。同样，在渗透胁迫后，各个拟南芥的AtCAT表达水平也均显著降低，转基因株系OE1和OE2的AtCAT基因表达水平下降幅度最大，分别较对照组下降了58.2%、66.1%（图7C）。综上可知，250 mmol/L甘露醇处理条件下转基因株系OE1和OE2体内AtSOD、AtPOD和AtCAT基因的表达水平显著均低于野生型和突变体（图7A~7C），说明过表达胡杨PeAnn1降低了拟南芥ROS清除酶基因的表达，导致不能清除体内过量的H2O2，抗旱性较弱。

    图  7  甘露醇胁迫对野生型（WT）、Annexin1突变体（atann1）和过表达PeAnn1（PeAnn1-OE1和PeAnn1-OE2）拟南芥抗氧化酶基因AtSOD、AtPOD和AtCAT表达的影响
    A. 甘露醇胁迫下AtSOD基因表达；B. 甘露醇胁迫下AtPOD基因表达；C. 甘露醇胁迫下AtCAT基因表达；内参基因为AtACTIN2，n = 3。A, AtSOD expression under mannitol stress; B, AtPOD expression under mannitol stress; C, AtCAT expression under mannitol stress; reference gene: AtACTIN2, n = 3.
    Figure  7.  Effects of mannitol stress on relative expression of antioxidant-enzyme genes (AtSOD, AtPOD and AtCAT) in wild-type (WT), atann1 and transgenic lines (PeAnn1-OE1, PeAnn1-OE2)
    下载: 全尺寸图片 幻灯片
    图  8  干旱和复水后野生型（WT）、Annexin1突变体（atann1）和过表达PeAnn1（PeAnn1-OE1和PeAnn1-OE2）拟南芥的光合参数变化
    A. 干旱及复水后拟南芥的生长状况；B~D. 正常浇水、干旱8 d、复水3 d拟南芥的PSⅡ最大光量子效率（Fv/Fm）、相对电子传递效率（ETR）、实际光合量子产量（ΦPSⅡ）的变化。A, plant performance of Arabidopsis thaliana after drought stress and rewatering; B−D, maximum photon efficiency of PSⅡ（Fv/Fm）, relative electron transfer efficiency（ETR）, and actual photosynthetic quantum yield（ΦPSⅡ） of Arabidopsis thaliana under normal watering, 8 days of drought stress and 3 days of rewatering.
    Figure  8.  Changes of photosynthetic parameters in wild-type (WT), atann1, and transgenic lines (PeAnn1-OE1, PeAnn1-OE2) under drought stress and rewatering
    下载: 全尺寸图片 幻灯片

    2.8   干旱胁迫下拟南芥叶绿素荧光参数的变化

    对各株系拟南芥进行土壤干旱处理，并测定叶绿素含量及荧光参数。干旱处理前，拟南芥生长状况没有明显差异，干旱胁迫8 d后植株失水枯黄，野生型和突变体拟南芥长势较好。复水3 d后，转基因株系死亡率较高（图8A）。

    叶绿素荧光参数反映植物光反应中心活性和光合生理状况。胁迫条件下PSⅡ的最大光量子效率（Fv/Fm）的变化可反映PSⅡ的光抑制程度[29]。如图8B所示，在干旱胁迫前，野生型、突变体和转基因株系OE1和OE2的Fv/Fm为0.706 ~ 0.747，其中转基因株系OE2的Fv/Fm最低，与其他拟南芥株系有显著差异（P < 0.05）。干旱胁迫8 d后，各株系的Fv/Fm有不同程度的降低，转基因株系的Fv/Fm下降更为显著；复水3 d后，PSⅡ最大光化学效率均迅速恢复，但转基因株系的恢复幅度不及野生型和突变体（图8B）。说明PeAnn1转基因株系的光反应中心在干旱胁迫条件下受到逆境的伤害较大。

    干旱处理前，转基因株系和突变体的相对电子传递效率（ETR）、实际光合量子产量（ΦPSⅡ）均低于野生型（图8C、8D）。干旱8 d后，转基因株系的ETR和ΦPSⅡ值均大幅度降低，明显高于野生型、突变体株系的下降幅度（图8C、8D）。复水3 d后，各个拟南芥株系的相对电子传递效率（ETR）和实际光合量子产量（ΦPSⅡ）均快速恢复，但转基因株系的恢复程度远低于野生型、突变体拟南芥（图8C、8D）。这说明，与光反应中心相似，PeAnn1转基因株系的捕光蛋白复合体在干旱胁迫下被严重破坏，电子传递受阻，光能转化效率降低。

    2.9   干旱胁迫下拟南芥叶绿素含量的变化

    叶绿素相对含量采用叶绿素仪（SPAD 仪）原位测定[22]。由图9可以看出，在干旱前，突变体atann1的SPAD值显著高于野生型和转基因株系OE1和OE2（P < 0.05）。干旱8 d后，SPAD值均显著降低，其中转基因株系下降幅度最大。复水后第3天，各株系拟南芥的SPAD值均有所增加，但仍远低于干旱前的水平，特别是转基因株系。这表明干旱胁迫使转基因拟南芥叶绿素含量显著下降。与转基因株系相比，突变体atann1能在干旱和复水期间维持较高的叶绿素含量。

    图  9  干旱和复水后野生型（WT）、Annexin1突变体（atann1）和过表达PeAnn1（PeAnn1-OE1和PeAnn1-OE2）拟南芥的叶绿素相对含量
    拟南芥正常浇水、干旱8 d、复水3 d叶绿素SPAD值的变化。Chlorophyll SPAD values of Arabidopsis thaliana under normal watering, 8 days of drought stress and 3 days of rewatering.
    Figure  9.  Changes of relative chlorophyll content in wild-type (WT), atann1 and transgenic plant (PeAnn1-OE1, PeAnn1-OE2) under drought stress and rewatering
    下载: 全尺寸图片 幻灯片

    3.   讨　　论

    渗透胁迫处理诱导了胡杨叶片PeAnn1基因瞬时的上调表达（图1），表明PeAnn1参与了胡杨对水分胁迫的响应。为了明确胡杨PeAnn1在植物耐受水分胁迫中的作用，本文对野生型（WT）、Annexin1突变体（atann1）和PeAnn1转基因株系OE1和OE2进行渗透胁迫和干旱胁迫处理，研究各株系的表型差异。结果表明，甘露醇处理后，WT、atann1和转基因株系OE1和OE2的种子萌发率和根长均受到抑制，且抑制作用随甘露醇浓度提高而增强（图3、4）。在高渗透胁迫（250 ~ 300 mmol/L甘露醇）条件下，转基因植株的萌发率、根长显著低于野生型和突变体。干旱胁迫下，转基因株系的光合参数下降幅度高于野生型和突变体（图8、9），这些表型结果说明，过表达胡杨PeAnn1使转基因植株对水分胁迫的敏感性提高。

    PeAnn1转基因株系对水分胁迫耐受性的下降与光合参数的下降有关。叶绿素含量、PSⅡ最大光量子效率（Fv/Fm）、相对电子传递效率（ETR）和实际光合量子产量（ΦPSⅡ）的高低反映植物的光合能力和生长状况[30]。SPAD 值是一项表征叶片光合活性的指标，它与单位叶面积的叶绿素含量及叶绿素的密度显著正相关[31]，因此，常用SPAD 值的大小来衡量叶片中叶绿素含量。干旱胁迫下SPAD显著降低（图9），说明在干旱条件下叶绿素蛋白质的合成受到抑制，叶绿素逐渐分解，光合色素含量下降。干旱还会破坏捕光蛋白复合体，使电子传递受阻，光能转化效率降低，导致Fv/Fm、ETR、ΦPSⅡ显著下降（图8）。复水后，光合参数有所回升，说明光系统并未完全受到损伤。但与野生型和突变体相比，转基因株系修复损伤能力最弱，这与其干旱期间光系统受损的程度有关（图8、9）。

    在逆境胁迫下,植株体内会生成大量的超氧自由基和过氧化氢，SOD作为植物抗氧化系统的第一道防线，通过歧化超氧阴离子生成H2O2，生成的H2O2进一步被CAT、POD等分解消除，保护细胞免受氧化伤害[32]。渗透胁迫下各株系H2O2水平提高（图5），野生型和突变体维持较高抗氧化酶活性，有助于活性氧的清除（图6）。但转基因株系中SOD、CAT、POD酶活较低（图6），不能有效清除活性氧（图5），减轻氧化胁迫对细胞膜的伤害，难以保持光合系统的相对稳定。这与在印度芥菜[18]和小麦[19]中的研究结果不完全一致，表明不同物种膜联蛋白在抗氧化防御中的作用有所差异。众多研究表明，抗氧化酶活性与编码基因表达的表达量相关[33-34]。在渗透胁迫后PeAnn1转基因株系OE1和OE2抗氧化酶SOD、POD和CAT的活性和基因表达量均显著低于野生型和突变体（图6、7），表明过表达PeAnn1株系在干旱胁迫下不能及时清除体内大量的活性氧，破坏了细胞内蛋白质的结构和功能。

    关于膜联蛋白在植物耐受水分逆境中的作用有不同的研究报道。Lee等人[16]用甘露醇模拟渗透胁迫，发现过表达AtAnn1表现强耐旱性。这与Konopka-Postupolska等人[17]报道的结果类似，但不同于Huh等人[35]研究的结果。PeAnn1基因序列与拟南芥AtAnn1有较大差别（图2），因此，PeAnn1在调控植物抗旱性中作用可能有所不同。王希等人[36]通过超表达GsANN发现转基因拟南芥提高了对盐胁迫和干旱胁迫的敏感性，却志群等人[37]研究水稻OsAnn8基因可能通过负调控方式参与水稻对干旱胁迫的响应。这些研究与本文结论类似。另外，Liao等[38]发现拟南芥MYB30转录因子调控植物胁迫响应的信号转导机制。拟南芥MYB30通过结合ANN1和ANN4基因的启动子，抑制ANN1和ANN4基因表达，从而影响ANNs介导的特异胞内钙信号的产生，以此调控植物的氧化胁迫和热胁迫响应过程。胡杨PeAnn1是否也与相关转录因子互作，通过调控胞内钙信号影响抗氧化酶基因的表达，有待深入研究。

    综上所述，本文从生理和分子转录水平上揭示了胡杨PeAnn1基因对植物抗旱性的影响，研究结果有助于阐明胡杨的抗旱机制。

  • 图  1   各个林分变量与地上生物量和总生物量的散点图

    Figure  1.   Scatter plots of stand variables and aboveground and total biomass

    下载: 全尺寸图片 幻灯片

    图  2   对数转换线性回归模型残差图（n = 917）

    Figure  2.   Residual plot of predictions and observations based on the logarithmic transformation linear regression model (n = 917)

    下载: 全尺寸图片 幻灯片

    图  3   人工神经网络模型（模型25）残差图（n = 184）

    Figure  3.   Scatter plots of predictions and observations based on artificial neural network model (model 25) (n = 184)

    下载: 全尺寸图片 幻灯片

    图  4   人工神经网络模型（模型25）训练集与验证集预测结果对比图

    Figure  4.   Artificial neural network model (model 25) prediction comparison of training set and validation set

    下载: 全尺寸图片 幻灯片

    图  5   人工神经网络模型（模型34）残差图（n = 184）

    Figure  5.   Residual plot of predictions and observations based on artificial neural network model (model 34) (n = 184)

    下载: 全尺寸图片 幻灯片

    图  6   人工神经网络模型（模型34）训练集与验证集预测结果对比图

    Figure  6.   Artificial neural network model (model 34) prediction comparison of training set and validation set

    下载: 全尺寸图片 幻灯片

    图  7   人工神经网络模型（模型34）残差图（n = 917）

    Figure  7.   Residual plot of predictions and observations based on artificial neural network model (model 34) (n = 917)

    下载: 全尺寸图片 幻灯片

    表  1   样地基本因子与生物量统计量

    Table  1   Summary statistics of sample plot basic variables and biomass

    林分因子　　　
    Stand factor　　　    		最大值
    Max.    		最小值
    Min.    		平均值
    Mean    		标准差
    Standard deviation    		变异系数
    Coefficient of variation/%
    年龄/a Age/year57 5 27 11 40.10
    郁闭度 Crown density 1.00 0.20 0.67 0.2030.09
    株数密度/（株·hm− 2） Stand density/(tree·ha− 1)4 033 267 1 257 675 53.64
    林分断面积/（m2·hm− 2） Stand basal area/(m2·ha− 1) 34.63 3.2513.46 6.6749.55
    蓄积量/（m3·hm− 2） Stand volume/(m3·ha− 1)284.3112.9286.2051.1659.35
    平均胸径 Quadratic mean diameter (D)/cm27.75.812.2 3.629.38
    平均树高 Mean tree height (H)/m23.08.415.9 2.717.16
    海拔 Altitude (HB)/m1 190 10 525 209 39.88
    坡度 Slope (PD)/(°)36 0 11 7 68.93
    地上生物量/（t·hm− 2） Aboveground biomass/(t·ha− 1)174.62 9.9060.0033.4055.67
    总生物量/（t·hm− 2） Total biomass/(t·ha− 1)213.2712.0373.6340.9655.63
    下载: 导出CSV

    表  2   地上生物量对数转换线性回归模型参数估计值

    Table  2   Estimated parameters of aboveground biomass based on log-linear regression model

    编号 No.异速方程 Allometric equationb0b1b2b3
    1lnY = b0 + b1lnD 0.319 8**1.463 8***−−
    2lnY = b0 + b1ln（D2 × H） − 0.676 8***0.599 2***−−
    3lnY = b0 + b1lnD + b2lnH− 3.188 2***0.148 2ns 2.450 3***−
    4lnY =b0 + b1ln（D2 × H × S）− 2.989 5***0.502 7***−−
    5lnY = b0 + b1lnD + b2lnH + b3lnS− 3.865 1***0.404 9***0.750 0***0.778 1***
    6lnY = b0 + b1ln（D2 × S）− 2.554 1***0.589 0***−−
    7lnY = b0 + b1lnD + b2lnS− 2.993 9***0.764 8***0.828 7***−
    8lnY =b0 + b1ln（D2 × H） + b2lnS− 3.372 0***0.314 4***0.803 6***−
    9lnY =b0 +b1ln（D2 × S） + b3lnH− 2.307 6***0.628 8***− 0.248 5** −
    注：D为平均胸径；H为平均树高；S为林分密度指数; b0、b1、b2、b3为模型参数；ns代表参数在0.05水平上不显著；**代表在0.01水平上显著；***代表在0.001水平上显著。下同。Y为地上生物量；Notes: D represents quadratic mean diameter; H represents mean tree height; S represents stand density index; b0, b1, b2, b3 are parameters of models; ns represents no significant difference at 0.05 level; ** represents significant difference at 0.01 level; *** represents significant difference at 0.001 level. The same below. Y represents aboveground biomass.
    下载: 导出CSV

    表  3   地上生物量对数线性回归模型拟合优度与精度检验

    Table  3   Fitting statistics of aboveground biomass based on log-linear regression model

    编号
    No.    		AIC决定系数
    R2    		调整决定系数
    Adj.R2    		均方根误差/(t·hm− 2)
    RMSE/(t·ha− 1)    		相对均方根误差
    RMSEr/%    		平均绝对误差/(t·hm− 2)
    MAE/(t·ha− 1)
    18 303.537 90.448 60.448 024.790 415.049 9 19.134 8
    28 219.147 20.500 90.500 423.584 114.317 6 18.190 9
    38 091.547 90.572 60.571 621.825 613.250 0 16.851 1
    47 301.056 00.826 20.826 013.918 38.449 610.278 4
    56 888.135 90.902 40.902 110.427 76.330 56.850 3
    67 240.699 10.835 90.835 713.525 08.210 89.940 4
    76 955.132 70.894 20.894 010.859 06.592 37.358 0
    86 910.318 90.900 30.900 110.542 26.400 07.005 9
    97 235.057 20.836 70.836 413.489 48.189 29.888 6
    下载: 导出CSV

    表  4   总生物量对数转换线性回归模型参数估计值

    Table  4   Estimated parameters of total biomass based on log-linear regression model

    编号 No.异速方程 Allometric equationb0b1b2b3
    10lnY = b0 + b1lnD 0.534 8***1.459 6***−−
    11lnY = b0 + b1ln（D2 × H）− 0.460 4** 0.597 7** −−
    12lnY = b0 + b1lnD + b2lnH− 2.986 0***0.139 2ns 2.459 3***−
    13lnY = b0 + b1ln（D2 × H × S）− 2.770 5***0.501 7***−−
    14lnY = b0 + b1lnD + b2lnH + b4lnS− 3.661 8***0.395 4***0.762 0***0.776 7***
    15lnY = b0 + b1ln（D2 × S）− 2.335 5***0.587 7***−−
    16lnY = b0 + b1lnD + b2lnS− 2.776 5***0.761 1***0.828 2***−
    17lnY = b0 + b1ln（D2 × H） + b2lnS− 3.153 5***0.313 1***0.803 0***−
    18lnY = b0 + b1ln（D2 × S） + b3lnH− 2.095 2***0.626 5***− 0.242 3** −
    注：Y为总生物量。Note： Y represents total biomass.
    下载: 导出CSV

    表  5   总生物量对数转换线性回归拟合优度与精度检验

    Table  5   Fitting statistics of total biomass based on log-linear regression model

    编号 No.AICR2Adj.R2RMSE/(t·hm− 2)
    RMSE/(t·ha− 1)    		RMSEr/%MAE/(t·hm− 2)
    MAE/(t·ha− 1)
    108 684.928 00.444 50.443 930.510 515.161 023.547 6
    118 600.925 10.497 00.496 429.033 914.427 322.397 4
    128 472.676 50.569 70.568 726.853 213.343 720.732 3
    137 696.929 90.821 50.821 417.292 88.593 012.781 9
    147 292.576 70.898 20.897 913.060 86.490 18.778 7
    157 638.875 60.831 10.830 916.825 88.360 912.362 3
    167 359.598 60.889 70.889 513.594 16.755 09.369 3
    177 315.456 50.895 90.895 713.207 46.562 98.973 9
    187 633.789 10.831 90.831 516.785 28.340 812.309 4
    下载: 导出CSV

    表  6   未加入立地因子的人工神经网络模型结果与拟合优度检验（基于测试集，n = 184）

    Table  6   Results and goodness of fitting statistics based on artificial neural network models without site factors (based on test set, n = 184)

    编号
    No.    		输入单元
    Input unit    		隐藏层数
    Number of hidden layers    		隐单元数
    Number of hidden units    		训练轮数
    Epoch    		R2Adj.R2
    AGBTotal BAGBTotal B
    19D2H240−202040.523 60.505 80.521 00.503 1
    20D2HS240−20960.721 80.712 40.720 30.710 8
    21D2HS340−40−201120.793 60.782 30.792 50.781 1
    22D2HS440−40−40−201310.790 40.780 70.789 20.779 5
    23D-H240−20920.626 30.617 20.622 20.612 9
    24D-S240−203500.900 90.896 60.899 80.895 4
    25D-H-S240−202910.910 90.907 20.909 40.905 6
    26D-H-S340−40−202590.909 90.906 40.908 40.904 9
    27D-H-S440−40−40−201780.895 50.892 60.893 70.890 8
    28D2HS-D-H-S240−203750.909 30.904 90.907 30.902 8
    29D2HS-D-H-S340−40−202620.908 30.903 60.906 20.901 4
    30D2HS-D-H-S440−40−40−201380.899 80.895 40.897 50.893 1
    注：AGB代表地上生物量（t/hm2），Total B代表总生物量（t/hm2）。下同。Notes: AGB stands for aboveground biomass (t/ha) and Total B stands for total biomass (t/ha). The same below.
    下载: 导出CSV

    表  7   未加入立地因子的人工神经网络模型精度检验（基于测试集，n = 184）

    Table  7   Accuracy test statistics based on artificial neural network models without site factors (based on test set, n = 184)

    编号 No.RMSE/(t·hm− 2)
    RMSE/(t·ha− 1)    		RMSEr/%MAE/(t·hm− 2)
    MAE/(t·ha− 1)
    AGBTotal BAGBTotal BAGBTotal B
    1924.874 530.791 015.458 815.735 59.559 111.827 7
    2019.006 923.491 811.812 312.004 16.946 6 8.630 2
    2114.785 218.412 3 9.832 6 9.870 25.236 9 6.507 1
    2216.501 120.514 310.255 010.482 75.998 0 7.459 6
    2322.029 727.102 013.690 813.848 98.709 110.745 9
    2411.345 514.087 4 7.050 9 7.198 53.999 8 5.012 9
    2510.758 913.345 2 6.686 4 6.819 33.692 8 4.736 1
    2610.815 713.400 0 6.721 6 6.847 33.719 4 4.710 5
    2711.650 414.353 8 7.240 4 7.334 73.746 1 4.744 5
    2810.851 413.508 8 6.743 8 6.902 93.668 1 4.687 1
    2910.912 613.601 0 6.781 9 6.950 03.855 4 4.887 5
    3011.409 914.166 4 7.091 0 7.238 94.286 3 5.376 1
    下载: 导出CSV

    表  8   加入立地因子的人工神经网络模型结果与拟合优度检验（基于测试集，n = 184）

    Table  8   Results and goodness of fit statistics based on artificial neural network models with site factors (based on test set, n = 184)

    编号
    No.    		输入单元
    Input unit    		隐藏层数
    Number of hidden layer    		隐单元数
    Number of hidden units    		训练轮数
    Epoch    		R2Adj.R2
    AGBTotal BAGBTotal B
    31D-H-S-HB-PD-PX-PW240−201130.905 90.901 40.893 00.887 8
    32D-H-S-HB-PD-PX-PW340−40−201600.898 40.895 10.884 50.880 7
    33D-H-S-HB-PD-PX-PW440−40−40−201280.881 50.880 30.865 20.863 9
    34D-H-S-HB-PX-PW240−201260.910 90.906 60.909 40.905 1
    35D-H-S-HB-PX-PW340−40−201790.898 30.895 60.885 10.882 0
    36D-H-S-HB-PX-PW440−40−40−201360.907 10.903 20.894 90.890 6
    37D-H-S-PX-PW240−201080.908 00.904 40.896 70.892 6
    38D-H-S-PX-PW340−40−20650.901 70.897 30.889 60.884 6
    39D-H-S-PX-PW440−40−40−201260.900 70.897 90.888 40.885 3
    注：HB代表海拔（m），PD代表坡度（°），PX代表坡向，PW代表坡位。Notes: HB stands for altitude, PD stands for slope (°), PX stands for slope aspect, and PW stands for slope position.
    下载: 导出CSV

    表  9   加入立地因子的人工神经网络模型精度检验（基于测试集，n = 184）

    Table  9   Accuracy test statistics based on artificial neural network models with site factors (based on test set, n = 184)

    编号 No.RMSE/(t·hm− 2)
    RMSE/(t·ha− 1)    		RMSEr/%MAE/(t·hm− 2)
    MAE/(t·ha− 1)
    AGBTotal BAGBTotal BAGBTotal B
    3111.053 313.758 06.870 67.030 23.760 54.818 1
    3211.484 914.185 07.137 67.248 43.883 24.943 2
    3312.406 315.153 97.710 27.743 54.216 15.212 7
    3410.758 013.385 16.685 86.839 73.611 94.662 4
    3511.491 614.152 07.141 77.231 53.731 34.704 2
    3610.987 313.625 36.828 36.962 43.903 24.921 7
    3710.929 913.541 26.792 66.919 43.707 04.742 8
    3811.298 614.039 97.021 87.174 33.756 44.809 0
    3911.358 613.998 27.059 17.152 93.846 94.933 7
    注：AGB代表地上生物量（t/hm2），Total B代表总生物量（t/hm2）。Notes: AGB stands for aboveground biomass (t/ha) and Total B stands for total biomass (t/ha).
    下载: 导出CSV

    表  10   各组最优模型的比较（n = 917）

    Table  10   Comparison of the 4 best models for each group (n = 917)

    模型 Model地上生物量 Aboveground biomass总生物量 Total biomass
    R2Adj.R2RMSE/(t·hm− 2)
    RMSE/(t·ha− 1)    		RMSEr/%MAE/(t·hm− 2)
    MAE/(t·ha− 1)    		R2Adj.R2RMSE/(t·hm− 2)
    RMSE/(t·ha− 1)    		RMSEr/%MAE/(t·hm− 2)
    MAE/(t·ha− 1)
    50.902 40.902 110.427 76.330 56.850 3
    140.898 20.897 913.060 86.490 18.778 7
    250.906 40.906 110.426 06.329 56.985 10.900 40.900 013.099 06.509 08.990 2
    340.916 10.914 1 9.870 55.992 26.737 70.911 00.908 912.383 76.153 68.625 0
    下载: 导出CSV
    • 参考文献(22)
  • [1] 罗云建. 森林生物量的估算方法及其研究进展[J]. 林业科学, 2009, 45(8):129−134. doi: 10.3321/j.issn:1001-7488.2009.08.023

    Luo Y J. Forest biomass estimation methods and their prospects[J]. Scientia Silvae Sinicae, 2009, 45(8): 129−134. doi: 10.3321/j.issn:1001-7488.2009.08.023
    [2] 雷相东, 张会儒, 牟惠生. 东北过伐林区蒙古栎林分相容性生物量模型研究[J]. 第四纪研究, 2010, 30(3):559−565. doi: 10.3969/j.issn.1001-7410.2010.03.14

    Lei X D, Zhang H R, Mu H S. Compatible stand biomass models of Mongolia oak forests in over logged forest regions, Northeast China[J]. Quaternary Sciences, 2010, 30(3): 559−565. doi: 10.3969/j.issn.1001-7410.2010.03.14
    [3] 程堂仁, 冯菁, 马钦彦, 等. 基于森林资源清查资料的林分生物量相容性线性模型[J]. 北京林业大学学报, 2007, 29(5):110−113. doi: 10.3321/j.issn:1000-1522.2007.05.022

    Cheng T R, Feng J, Ma Q Y, et al. Linear models with compatibility of stand biomass based on the forest resource inventory data[J]. Journal of Beijing Forestry University, 2007, 29(5): 110−113. doi: 10.3321/j.issn:1000-1522.2007.05.022
    [4] 董利虎, 李凤日. 大兴安岭东部主要林分类型乔木层生物量估算模型[J]. 应用生态学报, 2018, 29(9):2825−2834.

    Dong L H, Li F R. Stand-level biomass estimation models for the tree layer of main forest types in East Daxing ’an Mountains, China.[J]. Chinese Journal of Applied Ecology, 2018, 29(9): 2825−2834.
    [5] 欧光龙, 胥辉, 王俊峰, 等. 思茅松天然林林分生物量混合效应模型构建[J]. 北京林业大学学报, 2015, 37(3):101−110.

    Ou G L, Xu H, Wang J F, et al. Building mixed effect models of stand biomass for Simao pine (Pinus kesiya var. langbianensis) natural forest[J]. Journal of Beijing Forestry University, 2015, 37(3): 101−110.
    [6] 赵嘉诚, 李海奎. 杉木单木和林分水平地下生物量模型的构建[J]. 林业科学, 2018, 54(2):81−89.

    Zhao J C, Li H K. Establishment of below-ground biomass equations for Chinese fir at tree and stand level[J]. Scientia Silvae Sinicae, 2018, 54(2): 81−89.
    [7]

    Vahedi A A. Artificial neural network application in comparison with modeling allometric equations for predicting above-ground biomass in the Hyrcanian mixed-beech forests of Iran[J]. Biomass and Bioenergy, 2016, 88: 66−76. doi: 10.1016/j.biombioe.2016.03.020
    [8]

    Sileshi G W. A critical review of forest biomass estimation models, common mistakes and corrective measures[J]. Forest Ecology and Management, 2014, 329: 237−254. doi: 10.1016/j.foreco.2014.06.026
    [9] 董利虎. 东北林区主要树种及林分类型生物量模型研究[D]. 哈尔滨: 东北林业大学, 2015.

    Dong L H. Developing individual and stand-level biomass equations in Northeast China forest area[D]. Harbin: Northeast Forest University, 2015.
    [10]

    Nandy S, Singh R, Ghosh S, et al. Neural network-based modelling for forest biomass assessment[J]. Carbon Management, 2017, 8(4): 305−317. doi: 10.1080/17583004.2017.1357402
    [11]

    Stas S M, Rutishauser E, Chave J, et al. Estimating the aboveground biomass in an old secondary forest on limestone in the Moluccas, Indonesia: comparing locally developed versus existing allometric models[J]. Forest ecology and management, 2017, 389: 27−34. doi: 10.1016/j.foreco.2016.12.010
    [12]

    Ercanlı İ, Günlü A, Şenyurt M, et al. Artificial neural network models predicting the leaf area index: a case study in pure even-aged Crimean pine forests from Turkey[J]. Forest Ecosystems, 2018, 5(1): 29. doi: 10.1186/s40663-018-0149-8
    [13]

    Özçelik R, Diamantopoulou M J, Eker M, et al. Artificial neural network models: an alternative approach for reliable aboveground Pine tree biomass prediction[J]. Forest Science, 2017, 63(3): 291−302.
    [14] 王轶夫, 孙玉军, 郭孝玉. 基于BP神经网络的马尾松立木生物量模型研究[J]. 北京林业大学学报, 2013, 35(2):17−21.

    Wang Y F, Sun Y J, Guo X Y. Single-tree biomass modeling of Pinus massoniana based on BP neural network[J]. Journal of Beijing Forestry University, 2013, 35(2): 17−21.
    [15] 解雅麟, 王海燕, 雷相东. 基于3-PG模型的长白落叶松人工林生长和生物量模拟[J]. 南京林业大学学报(自然科学版), 2018, 42(1):141−148.

    Xie Y L, Wang H Y, Lei X D. Growth and biomass simulation of Larix olgensis plantations based on 3-PG model[J]. Journal of Nanjing Forestry University(Natural Science Edition), 2018, 42(1): 141−148.
    [16] 国家林业与草原局. 立木生物量模型及碳计量参数落叶松[S]. 北京: 中国标准出版社, 2016.

    National Forestry and Grassland Administration of the People ’s Republic of China. Tree biomass models and related parameters to carbon accounting for Larix[S]. Beijing: Standards Press of China, 2016.
    [17]

    Goodfellow I, Bengio Y, Courville A, et al. Deep learning[M]. Cambridge: MIT Press, 2016.
    [18]

    Uzun H, Yıldız Z, Goldfarb J L, et al. Improved prediction of higher heating value of biomass using an artificial neural network model based on proximate analysis[J]. Bioresource Technology, 2017, 234: 122−130. doi: 10.1016/j.biortech.2017.03.015
    [19]

    Tieleman T, Hinton G. Lecture 6.5-rmsprop: divide the gradient by a running average of its recent magnitude[J]. COURSERA: Neural Networks for Machine Learning, 2012, 4(2): 26−31.
    [20]

    Xiao X, White E P, Hooten M B, et al. On the use of log-transformation vs. nonlinear regression for analyzing biological power laws[J]. Ecology, 2011, 92(10): 1887−1894. doi: 10.1890/11-0538.1
    [21] 符利勇, 雷渊才, 孙伟, 等. 不同林分起源的相容性生物量模型构建[J]. 生态学报, 2014, 34(6):1464−1470.

    Fu L Y, Lei Y C, Sun W, el al. Development of compatible biomass models for trees from different stand origin[J]. Acta Ecologica Sinica, 2014, 34(6): 1464−1470.
    [22]

    Özçelik R, Diamantopoulou M J, Brooks J R, et al. Estimating tree bole volume using artificial neural network models for four species in Turkey[J]. Journal of Environmental Management, 2010, 91(3): 742−753.
    • 相关文章
    • 施引文献(3)

  • 期刊类型引用(3)

    1. 张欢欢，范挺秀，高双成，范丙友，史国安. 2个牡丹切花品种核型、花药形态与花粉特性的比较. 浙江农业学报. 2025(02): 329-337 . 百度学术
    2. 孙榕泽，宋焕芝，蒋珈琦，于晓南. 三倍体芍药品种染色体制片优化和核型分析. 河北农业大学学报. 2024(03): 46-55 . 百度学术
    3. 段英姿，客绍英，王晓英，张胜珍，马艳芝. 15个油用紫斑牡丹品种核型及亲缘关系分析. 种子. 2023(07): 96-104 . 百度学术

    其他类型引用(0)

    • 资源附件(0)
图(7)  /  表(10)
计量
  • 文章访问数:  4894
  • HTML全文浏览量:  875
  • PDF下载量:  131
  • 被引次数: 3
出版历程
  • 收稿日期:  2019-01-14
  • 修回日期:  2019-03-12
  • 网络出版日期:  2019-04-29
  • 发布日期:  2019-04-30

目录

摘要
HTML全文

  • 1.   材料与方法

  • 1.1   实验材料

  • 1.2   PeAnn1序列与进化树分析

  • 1.3   渗透胁迫下拟南芥耐旱性检测

  • 1.3.1   种子萌发率
  • 1.3.2   根　长

  • 1.4   渗透胁迫下转基因拟南芥H2O2含量的测定

  • 1.5   渗透胁迫下转基因拟南芥抗氧化酶活性的测定

  • 1.6   转基因拟南芥中抗逆相关基因表达水平的检测

  • 1.7   自然干旱条件下拟南芥叶绿素荧光参数的测定

  • 1.7.1   叶绿素荧光参数
  • 1.7.2   叶绿素相对含量

  • 1.8   数据分析

  • 2.   结果与分析

  • 2.1   渗透胁迫诱导胡杨PeAnn1基因表达

  • 2.2   胡杨PeAnn1序列与进化树分析

  • 2.3   渗透胁迫下拟南芥的萌发率

  • 2.4   渗透胁迫下拟南根长的变化

  • 2.5   渗透胁迫下拟南芥H2O2含量的变化

  • 2.6   渗透胁迫下拟南芥抗氧化物酶活性的变化

  • 2.7   渗透胁迫下抗氧化酶基因表达量的变化

  • 2.8   干旱胁迫下拟南芥叶绿素荧光参数的变化

  • 2.9   干旱胁迫下拟南芥叶绿素含量的变化

  • 3.   讨　　论

参考文献(22)
相关文章
施引文献(3)
资源附件(0)

/

返回文章
返回

版权所有 © 2013 《北京林业大学学报》编辑部

地址：北京市海淀区清华东路35号邮政编码：100083

电话：010-62337673Email：bldxed@bjfu.edu.cn 京ICP备05066833号-1

本系统由北京仁和汇智信息技术有限公司开发  百度统计