Marginal effect of radial growth of Cryptomeria japonica to climate in Lushan Mountain of southwestern China
-
摘要:目的 深入认知树轮−气候要素之间的关系，揭示树木径向生长的主要限制性因素的相对贡献率及边际效应，以提升树轮在树轮生态学、树轮气候学研究中的应用价值。方法 基于庐山地区日本柳杉树轮宽度资料和庐山气象站气象资料，采用树轮气候学方法研制树轮宽度年表，基于相关分析方法初步识别树木径向生长的主要限制性因素，进而利用增强回归树分析方法揭示了庐山日本柳杉径向生长的主要影响因素相对贡献及边际效应。结果 正向影响庐山日本柳杉径向生长的因素按重要程度从大到小依次为当年1—3月平均最低气温（20.66%）、上年2—11月相对湿度（15.4%）、当年2—3月平均水汽压（9.47%），负向影响庐山日本柳杉径向生长的因素按重要程度从大到小依次为上年11月日照时数（20.81%）、上年5月最大日降水量（20.54%）、当年7月平均气温（13.11%）；树轮−气候之间的关系在阈值范围之内具有较好的线性关系，阈值范围外则不具有线性关系。结论 庐山日本柳杉径向生长受多种气候要素的综合影响，任一要素的影响均不是简单的线性关系，均存在明显的阈值效应。在分析树木径向生长对气候要素响应及进行树轮气候重建时，对边际效应问题应予以重视，以增强树轮−气候间关系的可信度及气候重建的可靠性。Abstract:Objective This paper aims to further understand the relationship between tree radial growth and climatic factors, reveal the relative importance as well as marginal effect of main driving climatic factors of tree radial growth, so as to enhance the application value of tree ring in the research of dendroecology and dendroclimatology.Method Based on the tree ring width data of Cryptomeria japonica in Lushan Mountain area and the meteorological data of Lushan Mountain meteorological station, the tree ring width chronology was developed by the dendroclimatology method. The main driving climatic factors of tree radial growth were initially identified based on the correlation analysis method, and then the relative importance and marginal effect of the main influencing factors on Cryptomeria japonica radial growth were revealed by utilizing boosted regression tree (BRT) method.Result The research results showed that the climatic factors that positively affected the radial growth of Cryptomeria japonica in Lushan Mountain, in descending order of importance, were the average minimum temperature from January to March of the current year (20.66%), the relative humidity from February to November of the previous year (15.4%), and the average vapor pressure from February to March of the current year (9.47%); the climatic factors that negatively affected the radial growth of Cryptomeria japonica in Lushan Mountain, in descending order of importance, were the sunshine hours in November of the previous year (20.81%), the maximum daily precipitation in the May of the previous year (20.54%) and the average temperature in July of the current year (13.11%). The relationship between tree ring and climate had a good linear relationship within the threshold range, and there was no linear relationship outside the threshold range.Conclusion The radial growth of Cryptomeria japonica in Lushan Mountain is affected by many climatic factors, and the influence of each factor has obvious marginal effect. It is important to pay close attention to the problem of marginal effects when performing tree ring based climate reconstruction, which will enhance the reliability of tree ring climate relationship and climate reconstruction.
-
随着全球气候变暖，极端气候频现，干旱、盐害等非生物胁迫已成为威胁全球农林业生产和发展的一个重大挑战[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′) AtACTIN2 GGTAACATTGTGCTCAGTGGTGG AACGACCTTAATCTTCATGCTGC AtSOD AGGAAACATCACTGTTGGAGAT GAGTTTGGTCCAGTAAGAGGAA AtCAT AGGATCAAACTTTGAGGGGTAG CTTGTGGTTCCTGGAATCTACT AtPOD CGTGCCCTTCATATTGTTGG GACGCCATCAACAACGAGTC 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 leaves2.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 proteins2.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 rewatering2.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 rewatering3. 讨 论
渗透胁迫处理诱导了胡杨叶片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基因对植物抗旱性的影响，研究结果有助于阐明胡杨的抗旱机制。
-
![]()
图 3 差值年表与气候要素相关关系
TMINC1-C3为 当年1—3月平均最低气温；TEMC7为当年7月平均气温；SUNP11为上年11月日照时数；WVPC2-C3为当年2—3月平均水汽压；RHP2-P11为上年2—11月相对湿度；PREP5为上年5月最大日降水量。下同。TMINC1-C3, average minimum temperature from January to March of the current year; TEMC7, average temperature in July of the current year; SUNP11, sunshine hours in November of the previous year; WVPC2-C3, average water vapor pressure from February to March of the current year; RHP2-P11, relative humidity from February to November of the previous year; PREP5, maximum daily precipitation in May of the previous year. The same below.
Figure 3. Relationship between difference chronology and climate factors
表 1 年表统计特征及公共区间分析结果
Table 1 Descriptive statistics of tree ring width chronology and the results of common interval analysis
年表统计量
Chronological statistics (1932—2018)公共区间统计量
Common interval statistics (1950—2016)平均敏感度
Mean
sensitivity标准差
Standard
deviation一阶自相关系数
First-order
autocorrelation
coefficient平均相关系数
Mean
correlation
coefficient树内相关系数
Within-tree
correlation
coefficient树间相关系数
Between-trees
correlation
coefficient信噪比
Signal-to-noise
ratio样本总体代表性
Total
representativeness
of sample第1主成分
解释方差量
Variance explained
by the first
principal component0.133 0.175 0.474 0.363 0.557 0.359 24.46 0.961 39.3% 
下载: 导出CSV
-
[1] Zhang T, Liu Y, Zhang R, et al. Tree-ring width based streamflow reconstruction for the Kaidu River originating from the central Tianshan Mountains since A. D. 1700[J]. Dendrochronologia, 2020, 61: 1−13.
[2] Chen Y, Gagen M H, Chen F, et al. Precipitation variations recorded in tree rings from the upper Salween and Brahmaputra River valleys, China[J]. Ecological Indicators, 2020, 113: 1−10.
[3] 蔡秋芳, 刘禹. 湖北麻城马尾松树轮宽度对气候的响应及1879年以来6—9月平均最高气温重建[J]. 科学通报, 2013, 58(增刊 I):169−177. Cai Q F, Liu Y. The June–September maximum mean temperature reconstruction from Masson pine (Pinus massoniana Lamb.) tree rings in Macheng, southeast China since 1879 AD[J]. Chinese Science Bulletin (Chinese Version), 2013, 58(Suppl. I): 169−177.
[4] Tei S, Sugimoto A, Yonenobu H, et al. Growth and physiological responses of larch trees to climate changes deduced from tree-ring widths and δ13C at two forest sites in eastern Siberia[J]. Polar Science, 2014, 8(2): 183−195.
[5] Bhuyan U, Zang C, Menzel A. Different responses of multispecies tree ring growth to various drought indices across Europe[J]. Dendrochronologia, 2017, 44: 1−8.
[6] 张辰华, 李书恒, 白红英, 等. 太白山地区7月NDVI多尺度周期变化及其对气候因子的响应[J]. 资源科学, 2019, 41(11):2131−2143. doi: 10.18402/resci.2019.11.15 Zhang C H, Li S H, Bai H Y, et al. Multi-scale periodic variation of NDVI in July and its response to climatic factors in the Taibai Mountain area[J]. Resources Science, 2019, 41(11): 2131−2143. doi: 10.18402/resci.2019.11.15
[7] Briffa K R, Jones P D, Schweingruber F H, et al. Unusual 20th-century summer warmth in a 1 000-year temperature record from Siberia[J]. Nature, 1995, 376: 156−159.
[8] Savva Y, Oleksyn J, Reich P B, et al. Interannual growth response of Norway spruce to climate along an altitudinal gradient in the Tatra Mountains, Poland[J]. Trees-Structure and Function, 2006, 20(6): 735−746.
[9] 白天军, 刘苑秋, 温林生, 等. 庐山日本柳杉早材与晚材年轮宽度对气候变化的响应[J]. 北京林业大学学报, 2020, 42(9):61−69. doi: 10.12171/j.1000-1522.20190439 Bai T J, Liu Y Q, Wen L S, et al. Response of earlywood and latewood ring width of Cryptomeria japonica to climate change in Lushan Mountain, eastern China[J]. Journal of Beijing Forestry University, 2020, 42(9): 61−69. doi: 10.12171/j.1000-1522.20190439
[10] Ridgeway G. Generalized boosted models: a guide to the GBM package[J]. Compute, 2005, 1: 1−12.
[11] Lewin W C, Mehner T, Ritterbusch D, et al. The influence of anthropogenic shoreline changes on the littoral abundance of fish species in German lowland lakes varying in depth as determined by boosted regression trees[J]. Hydrobiologia, 2014, 724(1): 293−306.
[12] 焦琳琳, 常禹, 申丹, 等. 利用增强回归树分析中国野火空间分布格局的影响因素[J]. 生态学杂志, 2015, 34(8):2288−2296. Jiao L L, Chang Y, Shen D, et al. Using boosted regression trees to analyze the factors affecting the spatial distribution pattern of wildfire in China[J]. Chinese Journal of Ecology, 2015, 34(8): 2288−2296.
[13] 葛跃, 王明新, 孙向武, 等. 基于增强回归树的城市PM2.5日均值变化分析: 以常州为例[J]. 环境科学, 2017, 38(2):485−494. Ge Y, Wang M X, Sun X W, et al. Variation analysis of daily PM2.5 concentrations based on boosted regression tree: a case study in Changzhou[J]. Environmental Science, 2017, 38(2): 485−494.
[14] Gu H L, Wang J, Ma L J, et al. Insights into the BRT (boosted regression trees) method in the study of the climate-growth relationship of Masson pine in subtropical China[J]. Forests, 2019, 10(3): 1−20.
[15] Cook E R, Kairiukstis L A K. Methods of dendrochronology[M]. Dordrecht: Kluwer Academic Publishers, 1990.
[16] 曹仁杰, 尹定财, 田昆, 等. 丽江老君山海拔上限长苞冷杉(Abies georgei)和云南铁杉(Tsuga dumosa)径向生长对气候变化的响应[J]. 生态学报, 2020, 40(17):6067−6076. Cao R J, Yin D C, Tian K, et al. Response of radial growth of Abies georgei and Tsuga dumosa to climate change at upper distributional limits on Laojun Mountain, Lijiang, Yunnan, China[J]. Acta Ecologica Sinica, 2020, 40(17): 6067−6076.
[17] 郑壮鹏, 赵思媛, 周非飞, 等. 亚洲树轮平均敏感度分布特征及其影响因素[J]. 应用生态学报, 2019, 30(3):805−813. Zheng Z P, Zhao S Y, Zhou F F, et al. Variations of mean sensitivity of tree rings in Asia and their influencing factors[J]. Chinese Journal of Applied Ecology, 2019, 30(3): 805−813.
[18] 乔晶晶, 王童, 潘磊, 等. 不同海拔和坡向马尾松树轮宽度对气候变化的响应[J]. 应用生态学报, 2019, 30(7):2231−2240. Qiao J J, Wang T, Pan L, et al. Responses of radial growth to climate change in Pinus massoniana at different altitudes and slopes[J]. Chinese Journal of Applied Ecology, 2019, 30(7): 2231−2240.
[19] 阮超越, 董志鹏, 李颖俊, 等. 闽北地区柳杉树轮宽度变化对气候因子的响应[J]. 亚热带资源与环境学报, 2020, 15(1):11−18. doi: 10.3969/j.issn.1673-7105.2020.01.003 Ruan C Y, Dong Z P, Li Y J, et al. Responses of tree-ring width of Cryptomeria to climate factors in northern Fujian[J]. Journal of Subtropical Resources and Environment, 2020, 15(1): 11−18. doi: 10.3969/j.issn.1673-7105.2020.01.003
[20] 王童, 孙玉军, 乔晶晶. 将乐林场马尾松树轮宽度对气候变化的响应[J]. 北京林业大学学报, 2019, 41(9):30−39. Wang T, Sun Y J, Qiao J J. Response of Pinus massoniana tree-ring width in the Jiangle Area of Fujian Province to climate change[J]. Journal of Beijing Forestry University, 2019, 41(9): 30−39.
[21] 陈秋艳, 勾晓华, 张军周, 等. 树轮宽度指示的神农架地区过去172年的冬春季温度变化[J]. 第四纪研究, 2015, 35(5):1145−1154. doi: 10.11928/j.issn.1001-7410.2015.05.10 Chen Q Y, Gou X H, Zhang J Z, et al. Tree-ring indication of winter-spring temperature in Shennongjia Mountain area of central china over the past 172 years[J]. Quaternary Sciences, 2015, 35(5): 1145−1154. doi: 10.11928/j.issn.1001-7410.2015.05.10
[22] 李宗善, 刘国华, 傅伯杰, 等. 川西卧龙国家级自然保护区树木生长对气候响应的时间稳定性评估[J]. 植物生态学报, 2010, 34(9):1045−1057. doi: 10.3773/j.issn.1005-264x.2010.09.005 Li Z S, Liu G H, Fu B J, et al. Evaluation of temporal stability in tree growth-climate response in Wolong National Natural Reserve, western Sichuan, China[J]. Chinese Journal of Plant Ecology, 2010, 34(9): 1045−1057. doi: 10.3773/j.issn.1005-264x.2010.09.005
[23] 钱恒君, 蔡秋芳, 刘禹. 湖南慈利马尾松树轮宽度对气候的响应及1854年以来冬季极端低温重建[J]. 地球环境学报, 2019, 10(6):543−555. Qian H J, Cai Q F, Liu Y. Response of tree-ring width of Pinus massoniana to climate change and winter extreme minimum temperature reconstruction since 1854 in Cili, Hunan Province of subtropical China[J]. Journal of Earth Environment, 2019, 10(6): 543−555.
[24] 尚华明, 尹仔锋, 魏文寿, 等. 基于树木年轮宽度重建塔里木盆地西北缘水汽压变化[J]. 中国沙漠, 2015, 35(5):1283−1290. doi: 10.7522/j.issn.1000-694X.2014.00116 Shang H M, Yin Z F, Wei W S, et al. Tree-ring based water vapor pressure reconstruction in northwestern edge of the Tarim Basin[J]. Journal of Desert Research, 2015, 35(5): 1283−1290. doi: 10.7522/j.issn.1000-694X.2014.00116
[25] 赵伯阳, 刘禹, 宋慧明, 等. 基于油松树轮宽度重建河北青龙过去123年平均相对湿度[J]. 地球环境学报, 2016, 7(5):509−520. doi: 10.7515/JEE201605008 Zhao B Y, Liu Y, Song H M, et al. Mean relative humidity reconstruction based on the tree-ring width from Chinese pines since 1890 from the Qinglong region, China[J]. Journal of Earth Environment, 2016, 7(5): 509−520. doi: 10.7515/JEE201605008
[26] 黄小梅, 肖丁木, 秦宁生. 利用树木年轮重建公元1639−2013年青南高原5—9月相对湿度变化[J]. 干旱区地理, 2018, 41(5):1002−1008. Huang X M, Xiao D M, Qin N S. Tree-ring based reconstruction of relative humidity from May to September in southern Qinghai Plateau during AD 1639−2013[J]. Arid Land Geography, 2018, 41(5): 1002−1008.
[27] 王亚军, 张永, 邵雪梅. 河北太行山南段树木年轮指示的167年来相对湿度变化研究[J]. 生态学报, 2019, 39(12):4570−4578. Wang Y J, Zhang Y, Shao X M. The study of variation in relative humidity from May to July during the past 167 years in the southern Taihang Mountains in Hebei Province[J]. Acta Ecologica Sinica, 2019, 39(12): 4570−4578.
[28] Sun C, Liu Y, Song H, et al. Sunshine duration reconstruction in the southeastern Tibetan Plateau based on tree-ring width and its relationship to volcanic eruptions[J]. Science of the Total Environment, 2018, 628−629: 707−714.
[29] 郑广宇, 王文杰, 王晓春, 等. 帽儿山地区兴安落叶松人工林树木年轮气候学研究[J]. 植物研究, 2012, 32(2):191−197. doi: 10.7525/j.issn.1673-5102.2012.02.011 Zheng G Y, Wang W J, Wang X C, et al. Tree-ring climatology of Larix gmelinii in Maoershan Region, Northeastern China[J]. Bulletin of Botanical Research, 2012, 32(2): 191−197. doi: 10.7525/j.issn.1673-5102.2012.02.011
[30] 蒋延玲, 周广胜, 王玉辉, 等. 内蒙古地带性针茅植物对CO2和气候变化的适应性研究进展[J]. 生态学报, 2015, 35(14):4559−4569. Jiang Y L, Zhou G S, Wang Y H, et al. Advances in the adaptability of zonal Stipa plants to CO2 and climate change in Inner Mongolia[J]. Acta Ecologica Sinica, 2015, 35(14): 4559−4569.
[31] 曹春福, 周非飞, 董志鹏, 等. 福建戴云山台湾松树轮−气候响应的线性和非线性模式研究[J]. 亚热带资源与环境学报, 2016, 11(1):44−51. doi: 10.3969/j.issn.1673-7105.2016.01.007 Cao C F, Zhou F F, Dong Z P, et al. Study on nonlinear climate-growth patterns of Pinus taiwanens in Daiyun Mountain, Fujian Province[J]. Journal of Subtropical Resources and Environment, 2016, 11(1): 44−51. doi: 10.3969/j.issn.1673-7105.2016.01.007
-
期刊类型引用(3)
1. 张欢欢，范挺秀，高双成，范丙友，史国安. 2个牡丹切花品种核型、花药形态与花粉特性的比较. 浙江农业学报. 2025(02): 329-337 . 
百度学术
2. 孙榕泽，宋焕芝，蒋珈琦，于晓南. 三倍体芍药品种染色体制片优化和核型分析. 河北农业大学学报. 2024(03): 46-55 . 百度学术
3. 段英姿，客绍英，王晓英，张胜珍，马艳芝. 15个油用紫斑牡丹品种核型及亲缘关系分析. 种子. 2023(07): 96-104 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 1315
- HTML全文浏览量: 365
- PDF下载量: 95
- 被引次数: 3
