-
东北红豆杉(Taxus cuspidata)隶属于红豆杉科(Taxaceae),是第三纪孑遗乔木物种,也是国家Ⅰ级重点保护濒危植物,被列为《全国极小种群野生植物拯救保护工程规划(2011—2015年)》的典型极小种群野生植物[1]。在我国,东北红豆杉主要分布在东北三省,近10年来,由于其能够提炼出的紫杉醇具有抗癌功效,而被人类活动不同程度破坏和干扰,导致种群稳定性不断下降,如今又面临气候变化的威胁,其保护形势十分严峻[2]。万基中等[3]用MaxEnt模型预测东北红豆杉潜在分布区可能会发生变化,所以需要采取更有效的保护措施来应对气候变化对东北红豆杉造成的干扰。迁地保护是珍稀频危植物重返自然生境的回归引种及野生居群恢复重建的基础和原始材料的重要保障措施,同时也降低了物种灭绝的风险。在迁地保护的基础上,通过人工繁殖将植物引入到其原有的自然或半自然的栖息地,植物本身可以通过建立充足的遗传资源来适应并形成一个可以自然维持和更新的新种群[5]。然而当前制约珍稀濒危植物迁地保护的很重要的一个原因是引种的植物种类虽然多,但是成活率低,每年都要大量引进植物[6]。造成这种现象的原因是引种植物难以适应新生气候条件以及人为管理不善等[7]。因此,在开展东北红豆杉引种回归前进行相关的产地稳定性试验是十分必要的。
国内外专家对东北红豆杉其自身的生物学特性及濒危机制进行了很多的探索与研究,主要集中在东北红豆杉繁殖生长[8],种群特征及其与生境关系[9],紫杉醇的提取[10-11],以及光合生理特征[12]。然而对不同生境下不同产地东北红豆杉幼苗其生长特性以及稳定性的研究尚未见报道。近年来,植物与环境的关系一直是植物生理生态学研究的热点问题[13-14]。基因型和环境相互作用促使大部分树种受环境不同的影响而表现各不相同,有效的G×E分析对正确评估树种适应性和稳定性有重要作用,为树种的保护提供直接的依据[15]。目前,广泛用于评估植物稳定性和适应性的方法有很多,如秩次分析法、Franics和高稳系数法、Shukla模型、Kannenberg模型及AMMI模型(Additive main effects and multiplicative interaction)等[16]。其中AMMI模型是目前国际上流行的分析品种区域试验数据非常有效的模型,该模型与方差分析模型、线性回归模型相比,应用范围更广更有效,它的特点是把主成分分析与方差分析相结合,不仅能分析互作效应的显著性,还能借助双标图和互作效应值直观而定量地描绘出互作效应的特点,目前已广泛用于研究区域试验中农作物的基因型与环境互作(G×E)效应和产量性状分析等[17]。如AMMI模型已被广泛应用到水稻(Oryza sativa)[18]、旱地冬小麦(Triticum aestivum)[19]、大豆(Glycine max)[20]等作物的大范围、区域试验分析研究领域。AMMI模型也可以应用于树木的区域试验稳定性分析[21-22]。
相关研究表明,相同种源不同无性系东北红豆杉幼苗个体表现出不同的变异程度,且不同生长指标变异的规律性不明显[23]。本研究采用完全随机区组试验,进一步揭示5个不同纬度的东北红豆杉产地幼苗在4个生境下的变异规律,应用AMMI模型开展产地与地点间的互作效应分析,产地间稳定性分析,进而筛选出东北红豆杉适合不同试验点迁地保护的生产力高、稳定性强的优良产地幼苗。以推动濒危东北红豆杉野生资源的保护,为将来东北红豆杉的回归引种以及不同生境下的响应方式和生存策略研究提供理论基础。
-
东北红豆杉分布于中国吉林老爷岭、张广才岭及长白山区,其自然地理分布范围为32°30′~53°00′N、123°00′~155°00′E,垂直分布于海拔250~1 200 m之间[24]。山东、江苏等省有栽培。遵循东北红豆杉分布区的特点,从南到北,沿纬度由低到高的顺序,选择山东烟台、北京海淀、辽宁本溪、吉林汪清4个地点作为野外移栽试验地点。试验地设在山东烟台大季家镇极小种群示范基地、北京林业大学八家苗圃科研基地、辽宁本溪森林经营研究所和吉林汪清百草沟镇新田村。样地的地理气候条件概况如表 1。
表 1 各试验点的地理气候条件
Table 1. Geographical and climatic conditions of the test sites
序号
No.试验地点
Test site纬度
Latitude经度
Longitude海拔
Altitude/m年均温
Annual average
temperature/℃年降水量
Annual
precipitation/mm年日照时数
Annual sunshine
hour/hⅠ 山东Shandong Province 37°39′N 121°03′E 51 13.4 627.6 2 656 Ⅱ 北京Beijing 40°00′N 116°20′E 49 12.5 612.0 2 720 Ⅲ 辽宁Liaoning Province 40°52′N 123°54′E 239 8.3 784.4 2 678 Ⅳ 吉林Jilin Province 43°12′N 129°37′E 168 4.5 549.4 2 357 -
2017年3月,在东北红豆杉分布区选取具有代表性并且能确认产地的3年生实生苗。试验材料本着生长健壮、长势一致的植株进行选取。5个产地概况如表 2。
表 2 东北红豆杉幼苗产地概况
Table 2. Overview of the origin of Taxus cuspidata seedlings
产地
Source纬度
Latitude经度
Longitude苗木来源
Seedling source辽宁宽甸
Kuandian County of Liaoning
Province40°03′N 124°20′E 辽宁省丹东市宽甸满族自治县鸭绿江生态苗木繁育基地
Ecological seedling breeding base of Yalu River in Kuandian Manchu Autonomous
County, Dandong City, Liaoning Province辽宁本溪
Benxi City of Liaoning Province41°19′N 123°47′E 辽宁省本溪市桓仁满族自治县桓仁宏创苗木基地
Huanren Hongchuang seedling breeding base of Huanren Manchu Autonomous
County, Benxi City, Liaoning Province吉林和龙
Helong City of Jilin Province42°31′N 129°00′E 吉林省和龙甑峰苗木有限公司
Helongzhenfeng Seedling Breeding Limited Company in Jilin Province吉林汪清
Wangqing County of Jilin Province43°18′N 129°47′E 吉林省延边朝鲜族自治州汪清县东北红豆杉生物科技有限公司
Taxus cuspidata Biotechnology Limited Company of Wangqing County, Yanbian
Korean Autonomous Prefecture, Jilin Province黑龙江穆棱
Muling City of Heilongjiang
Province44°36′N 130°12′E 黑龙江省穆棱市振才伟业红豆杉苗木基地
Zhencai Weiye Taxus cuspidata seedling breeding base in Muling City,
Heilongjiang Province -
2017年4月初将5个产地东北红豆杉3年生幼苗移栽到4个试验地点。4个地点的试验设置均采用完全随机区组设计,每个试验地随机划分为4个区组,每组种植各产地幼苗40株(共计200株),株间距35 cm×35 cm,并在试验区四周分别种植2行保护行作为补充。生长期间只进行补水和除草管理。为保证数据可靠性,各试验地于2017年5月底和10月底分别进行2次调查测定,在各产地的每个组中随机抽取20株幼苗,分别测量幼苗的苗高、冠幅(植株两个相互垂直方向的冠辐直径的平均值)、当年抽梢长(每产地16株每株选5支当年抽梢进行测量)、侧枝数、最长侧枝长等5个生长指标,精确到0.1 cm,用游标卡尺测量地径,精确到0.01 mm。地栽当年5月测量其生长性状指标并统计苗木缓苗率,10月底再次测量并统计苗木保存率,各生长性状生长量为其当年生长季净增加量即10月底测量值减去5月底测量值。
-
采用Excel 2010进行数据整理,运用SPSS22.0进行方差分析和多重比较处理,对多重比较结果用Originpro2017制图;并利用DPS7.05软件中的AMMI模型进行基因与环境的互作效应分析[25],然后以AMMI模型分析结果为基础,参照吴为人[26]的方法计算稳定性参数Di指标,定量分析各产地东北红豆杉幼苗苗高生长量的稳定性。
(1) AMMI模型分析方法
$$ {{Y}_{ge}}=\mu +{{\alpha }_{g}}+{{\beta }_{e}}+\sum\limits_{r=1}^{n}{{{\lambda }_{r}}{{\varphi }_{gr}}{{\sigma }_{er}}+{{\rho }_{ge}}+{{\varepsilon }_{ger}}} $$ 式中:Yge是产地g在地点e中的的苗高平均观测值;μ代表苗高总平均值;αg为产地的平均偏差(各个产地的苗高平均值减去总的平均值);βe为地点的平均偏差(各个地点的苗高平均值减去总的平均值);λr为第r个主成分分析的特征值;φgr为第r个主成分的产地主成分得分;σer为第r个主成分的地点主成分得分;r为保留在模型中的主成分因子轴的总个数;ρge为提取n个IPCA轴后留下的残差;εger为误差。
(2) 稳定性参数计算
$$ {{D}_{i}}=\sqrt{\sum\limits_{i=1}^{n}{(\text{IPCA})_{_{ge}}^{2}}} $$ 式中:n为显著的IPCA个数;Di为第g个产地或第e个地点在n个IPCA上的主成分得分,Di大小为第g个产地或第e个地点的稳定性定量指标。
-
由表 3可知:山东试验点各产地幼苗地径、侧枝数、最长侧枝长和当年抽梢长生长量表现最好,分别高达2.60 mm、1.6个、2.38 cm和6.07 cm,高于北京试验点地径、侧枝数、最长侧枝长和当年抽梢长均值的3.77、1.88、10.35和1.30倍。辽宁试验点的苗高生长量表现最好,高于北京试验点苗高的2.97倍。4个试验点中,当年抽梢长变异范围都最大,从南到北变异范围略有减小,依次为3.3~10.2 cm、1.9~8.3 cm、3.3~9.2 cm、2.2~6.7 cm,分别相差近3.1、4.4、2.8、3.0倍。从变异系数来看,山东试验点的苗高、冠幅、最长侧枝长、当年抽梢长的变异系数最小,分别为29.15%、39.95%、36.32%和30.63%。北京试验点的地径、冠幅、侧枝数、当年抽梢长变异系数明显高于其他试验点,表明东北红豆杉对北京当年环境比较敏感,在北京选择优良单株方面相对具有较大的挖掘潜力。
表 3 不同试验点东北红豆杉幼苗生长性状的变异参数
Table 3. Variation parameters for growth traits of Taxus cuspidata seedlings at different sites
试验点
Test site生长性状
Growth trait均值
Mean标准差
Standard deviation变异范围
Range of variation变异系数
CV/%山东Shandong Province 苗高Seedling height/cm 4.14 1.21 2.50~6.80 29.15 地径Basal diameter/mm 2.60 0.85 1.56~5.56 32.56 冠幅Crown width/cm 1.93 0.77 0.70~3.50 39.95 侧枝数Number of lateral branch 1.60 1.23 0~4.00 76.94 最长侧枝长Maximum length of lateral branch/cm 2.38 0.86 1.10~4.30 36.32 当年抽梢长Length of shoot in current year/cm 6.07 1.86 3.30~10.20 30.63 北京Beijing 苗高Seedling height/cm 1.44 0.50 0.50~2.20 35.12 地径Basal diameter/mm 0.69 0.26 0.11~1.12 37.23 冠幅Crown width/cm 0.81 0.52 0.10~1.80 64.80 侧枝数Number of lateral branch 0.85 0.75 0~2.00 87.65 最长侧枝长Maximum length of lateral branch/cm 0.23 0.17 0~0.60 74.78 当年抽梢长Length of shoot in current year/cm 4.67 1.64 1.90~8.30 35.19 辽宁Liaoning Province 苗高Seedling height/cm 4.28 1.52 2.10~8.00 35.62 地径Basal diameter/mm 1.51 0.27 0.96~2.14 17.76 冠幅Crown width/cm 2.19 1.09 0.10~4.10 49.96 侧枝数Number of lateral branch 0.45 0.31 0~2.00 67.78 最长侧枝长Maximum length of lateral branch/cm 1.54 1.23 0.20~4.50 80.24 当年抽梢长Length of shoot in current year/cm 5.26 1.63 3.30~9.20 31.06 吉林Jilin Province 苗高Seedling height/cm 1.59 0.71 0.50~3.40 44.68 地径Basal diameter/mm 1.48 0.47 0.71~2.38 32.02 冠幅Crown width/cm 1.40 0.61 0.60~2.80 43.97 侧枝数Number of lateral branch 0.53 0.34 0~2.00 63.40 最长侧枝长Maximum length of lateral branch/cm 1.24 0.49 0.50~2.40 39.43 当年抽梢长Length of shoot in current year/cm 4.21 1.33 2.20~6.70 31.58 -
由表 4的联合方差分析结果可知:地径净生长量在产地间和地点间表现出极显著(P<0.01)的差异,而在产地与地点的互作效应上不显著,苗高、冠幅、侧枝数、最长侧枝长、当年抽梢长的净生长量在产地间和地点间以及产地与地点的互作效应上均表现出极显著(P<0.01)的差异,说明各性状受环境影响较强烈。根据方差分析结果, 计算各差异来源的方差分量组成比,用以说明地点、产地等各差异来源的变异在总变异中的比例。从方差分量来看,苗高、地径、冠幅、当年抽梢长生长量在地点间的方差分量分别为62.89%、64.37%、32.65%、52.34%,远大于其他差异来源,说明地点效应对东北红豆杉生长影响最大。侧枝数生长量在产地间的方差分量为58.17%,远大于其他差异来源,表明受遗传控制较大。最长侧枝长生长量在地点效应、产地与地点互作效应上方差分量较高为22.76%和26.69%,表明最长侧枝长由产地遗传基础与生长环境共同控制。
表 4 不同试验点东北红豆杉幼苗生长性状联合方差分析
Table 4. Joint analysis of variance for growth traits of Taxus cuspidata seedlings at different sites
生长性状
Growth trait地点Site 区组(地点)Block(site) 产地Source 产地×地点Source×site 随机误差
Random errordf F 方差分量
Variance
component/%df F 方差分量
Variance
component/%df F 方差分量
Variance
component/%df F 方差分量
Variance
component/%df 方差分量
Variance
component/%苗高Seedling height/cm 3 98.448** 62.89 3 2.457 1.57 4 15.329** 13.06 12 4.046** 10.34 57 12.14 地径Basal diameter/mm 3 67.622** 64.37 3 2.219 2.11 4 7.607** 9.65 12 1.518** 5.78 57 18.09 冠幅Crown width/cm 3 23.217** 32.65 3 0.415 0.58 4 6.735** 12.63 12 4.874** 27.42 57 26.72 侧枝数Number of lateral branch 3 13.304** 16.05 3 1.403 1.05 4 6.854** 58.17 12 3.900** 10.54 57 14.19 最长侧枝长Maximum length of
lateral branch/cm3 67.760** 22.76 3 1.508 2.40 4 10.848** 15.64 12 6.682** 26.69 57 32.51 当年抽梢长Length of shoot in
current year/cm3 21.493** 52.34 3 1.409 1.16 4 58.424** 11.17 12 3.527** 20.65 57 14.68 注:*差异显著,P<0.05;**差异极显著,P<0.01。Notes: * means significant difference,P<0.05;** means extremely significant difference,P<0.01. -
图 1结果显示:4个试验点间,东北红豆杉幼苗产地间在当年抽梢长生长量Ducan多重比较差异最大,其次是冠幅生长量;差异最小的是苗高生长量。对于不同产地的幼苗来说,汪清的幼苗在各试验点的生长表现均很突出,例如苗高生长量在山东、北京和辽宁试验点均为最大地径生长量在北京和辽宁2个试验点均最大,当年抽梢生长量在山东、北京和辽宁试验点均最大,本溪幼苗在各试验点生长表现次之,例如苗高生长量在吉林试验点最大,地径生长量在山东和吉林试验点均最大。另外,4个试验点间,苗高、地径、冠幅和当年抽梢长生长量最好的是辽宁和山东试验点,最低的是北京试验点。
-
由表 5可知,4个地点的平均缓苗率由高到低分别为98.55%(吉林)>98%(山东)>96.82%(辽宁)>88.79%(北京);平均保存率由高到低分别为97.49%(吉林)>92.16%(辽宁)>89.92%(山东)>85.64%(北京),从4个试验点整体来看,保存率呈随着纬度增加而增加的趋势。从缓苗率和保存率对比来看山东试验点的苗子在试验期间死亡率最高,尤其是宽甸产地的幼苗成活率下降19.75%,在调查过程中发现宽甸产地苗木干枯死亡比较多, 对北京地区进一步分析表明,宽甸产地幼苗成活率也最低。从整体来看,各试验点内保存率最低的都是宽甸产地幼苗,北京试验点最高的是和龙产地幼苗,而其他试验点最高的都是汪清产地幼苗。
表 5 不同试验点东北红豆杉幼苗缓苗率和保存率
Table 5. Prolonged seedling conservation and preservation rate of Taxus cuspidata seedlings at different sites
% 试验点
Test site产地
Source缓苗率Seedling conservation rate 保存率Preservation rate 均值±标准差
Mean± SD变异范围
Range of variation均值±标准差
Mean± SD变异范围
Range of variation山东Shandong Province KD 91.60±5.15 85.11~95.83 71.85±15.16 51.02~86.96 BX 100.00±0.00 100.00~100.00 91.40±3.07 88.64~95.65 HL 99.49±1.02 97.96~100.00 95.98±4.97 89.80~100.00 WQ 98.91±2.18 95.65~100.00 97.27±2.11 95.56~100.00 ML 100.00±0.00 100.00~100.00 93.09±5.77 84.78~97.87 均值Mean 98.00±4.01 85.11~100.00 89.92±11.75 51.02~100.00 北京Beijing KD 60.44±15.01 48.89~82.22 56.41±12.43 46.67~73.33 BX 86.51±4.26 82.50~92.50 82.94±4.13 76.74~85.00 HL 98.88±1.31 97.50~100.00 98.88±1.31 97.50~100.00 WQ 100.00±0.00 100.00~100.00 94.38±9.66 80.00~100.00 ML 98.13±2.39 95.00~100.00 95.63±5.15 90.00~100.00 均值Mean 88.79±16.62 48.89~100.00 85.64±17.37 46.67~100.00 辽宁Liaoning Province KD 90.49±5.08 85.42~95.83 78.40±6.32 72.92~87.50 BX 96.88±2.69 93.75~100.00 91.15±3.56 87.50~95.83 HL 99.48±1.04 97.92~100.00 98.43±2.00 95.83~100.00 WQ 100.00±0.00 100.00~100.00 100.00±0.00 100.00~100.00 ML 97.28±3.21 93.75~100.00 92.83±5.97 86.05~97.87 均值Mean 96.82±4.37 85.42~100.00 92.16±8.72 72.92~100.00 吉林Jilin Province KD 97.47±0.98 96.00~98.00 95.93±6.88 85.71~100.00 BX 95.78±7.24 85.00~100.00 96.52±4.72 90.00~100.00 HL 100.00±0.00 100.00~100.00 97.00±4.76 90.00~100.00 WQ 100.00±0.00 100.00~100.00 99.50±1.00 98.00~100.00 ML 99.51±0.98 98.04~100.00 98.49±1.98 95.83~100.00 均值Mean 98.55±3.39 85.00~100.00 97.49±4.14 85.71~100.00 -
苗高是评价苗木质量好坏的主要指标[27]。在苗高净生长量在基因和环境间以及两者互作效应上均达到极显著(P<0.01)水平的基础上,采用回归模型和AMMI模型对其稳定性进行分析(表 6)。结果如下:苗高生长量在产地、地点以及两者互作效应的方差分量为13.06%、62.89%和10.34%,说明对苗高生长量总变异起作用的大小顺序依次为环境、基因型和基因与环境互作。AMMI模型分析表明,前两个主成分因子IPCA的F值均达到极显著水平,将剩余的不显著主成分值合并为残差。IPCA1、IPCA2的平方和分别占互作平方和的90.85%、8.08%,即前2项互作总成分解释了98.93%,说明用AMMI模型进行稳定性分析具有代表性,IPCA1和IPCA2代表的互作部分能对东北红豆杉幼苗苗高的稳定性做出判断。
表 6 东北红豆杉苗高性状的方差分析及AMMI模型分析
Table 6. ANOVA and AMMI analysis for the seeding height growth of Taxus cuspidata
分析方法
Analytical method变异来源
Source of variationdf 平方和
Sum of
square方差分量
Variance
component/%均方
Mean
squareF P 方差分析Analysis
of variance总的Total 79 232.034 2.937 地点内区组Location within the block 3 3.642 1.57 1.214 2.457 0.072 0 总处理Total treatment 19 200.226 10.538 19.879 0.000 1 产地Source 4 30.297 13.06 7.574 15.329 0.000 1 地点Site 3 145.937 62.89 48.646 98.448 0.000 1 产地×地点Source× site 12 23.992 10.34 1.999 4.046 0.000 3 误差Error 57 28.165 0.494 AMMI模型AMMI
analysisPCA1 6 21.798 90.85 3.633 28.506 0.000 1 PCA2 4 1.940 8.08 0.485 3.805 0.008 0 残差Residual 2 0.255 0.127 误差Error 60 31.808 0.530 为进一步对东北红豆杉产地和地点的稳定性进行分析,以苗高平均生长量为x坐标轴,以IPCA1值为y坐标轴做图 2如下。在AMMI标图中,产地与试验点同时位于x坐标轴一侧的互作效应为积极促进,与之不同侧的互作则为反向减弱[17]。结果显示在x轴方向试验点间要比产地间的分散,说明同一产地东北红豆杉幼苗在不同试验点的苗高生长表现差异较大。在y轴方向上,山东和辽宁试验点对宽甸、和龙、汪清产地幼苗具有正向的促进作用,相应的北京和吉林试验点对本溪、穆棱产地幼苗也具有正向的促进作用。
图 2 不同产地和地点间苗高生长量的双标图
Figure 2. AMMI plots for seeding height growth of Taxus cuspidata sources at different sites
为定量辨别产地稳定性与试验点分辨力,依据式(2)计算方法求得稳定性参数Di值,其中Di值越小,说明产地的稳定性越好;Di值越大,表明试验点对产地来源的分辨力即越好[17]。由表 7值可知,在苗高平均生长量上,试验的5个产地稳定性排名按序为和龙>穆棱>宽甸>本溪>汪清,4个试验点的稳定性排名按序为:北京>山东>吉林>辽宁。
表 7 东北红豆杉幼苗在显著互作效应轴上的得分及稳定性参数
Table 7. Scores of significant IPCA axis and stability parameter for Taxus cuspidata seedlings
项目
Item变量
Variable苗高生长量平均值
Mean growth of
seedling height/cm离差
DeviationIPCA1 IPCA2 Di 产地Source KD 2.531 -0.328 0.325 0.691 0.764 BX 2.456 -0.403 -0.962 -0.130 0.970 HL 3.069 0.210 0.136 0.013 0.136 WQ 3.969 1.110 1.012 -0.423 1.097 ML 2.269 -0.590 -0.511 -0.150 0.533 试验点Test place 山东Shandong Province 4.140 1.281 0.473 -0.665 0.816 北京Beijing 1.435 -1.424 -0.637 0.018 0.637 辽宁Liaoning Province 4.275 1.416 1.002 0.472 1.107 吉林Jilin Province 1.585 -1.274 -0.838 0.175 0.856 注:表中IPCA栏的数字表示各个主成分值。Note: data under IPCA in the table represent the principal component value. 综合对图 2和表 7的数据进行分析,表明东北红豆杉的5个产地中选出生长良好且稳定的产地为和龙;生长良好但不够稳定的产地是汪清;生长量不高但适应能力强的产地是穆棱;生长较差且不够稳定的产地是本溪。
-
相关研究表明,南方红豆杉(Taxus wallichiana var. mairei)生长与温度呈显著正相关,在一定范围内温度越高对生长越有利[28]。山东试验点内东北红豆杉各指标的生长量方面均优于其他3个试验点,山东烟台试验点地处温带海洋性季风气候,本身生长期长且水热条件良好,土壤适宜,促进了东北红豆杉生长。而北京试验点的东北红豆杉生长量普遍低于另外3个试验点,这与该试验点所在地土壤板结,2017年降雨量小等有关,在不利的环境条件下东北红豆杉生长量必然受到影响,但是这为后续抗旱产地的选择提供可能。本试验把采自不同环境条件下的东北红豆杉种植在4个相同环境的试验点下,发现试验点内的表型性状变异系数平均高达30%以上,说明东北红豆杉的表型变异不仅仅来源于表型可塑性,其更具有很高的遗传分化,南方红豆杉相关研究也验证了苗期性状为高度遗传[29-30]。各性状在地点间与产地间均存在极显著差异(P<0.01),说明同一产地东北红豆杉幼苗在不同的环境条件下的性状生长表现都明显不同,不同产地幼苗在同一试点内性状生长也各不一样,产地与试验间存在明显的互作效应。这一结果与王艺等[31]对南方红豆杉不同产地苗期性状产地变异的研究结果相一致,说明通过产地间在不同生境表现筛选出适应当地速生优良产地的潜力很大。
-
交互移植试验是揭示环境与遗传因素对植物性状作用机制的重要手段;也是预测植物在气候变化背景下分布范围变迁的重要研究方法之一[32]。经分析,辽宁和汪清试验点属于交互移植试验,且产地与试验点的互作效应对东北红豆杉的各项生长指标都存在显著影响,在辽宁本溪试验点中,当地本溪产地幼苗没有外地产地在生长中表现出显著的优势,同样在吉林汪清试验点中,当地汪清产地幼苗也没有外地产地在生长中表现出显著的优势,那么可以认为东北红豆杉没有存在局域适应现象,说明东北红豆杉可以尝试往南迁地保护。与南方红豆杉迁地保护有相似的现象,在地理上超出了自然分布区,说明东北红豆杉原生境正在或已经变成了“致危生境”[33]。
-
不同产地的幼苗保存率也是评价其稳定性大小的一项重要指标。经分析,各试验点内产地间保存率最低的都是宽甸产地,而在山东、辽宁、吉林试验点,最高的都是汪清产地,北京试验点最高的则是和龙产地,说明宽甸产地是成活率表现最差的幼苗,汪清产地相比是表现最好的幼苗。经分析,幼苗的保存率高低受试验点的地理气候条件影响较大,宽甸产地的幼苗本身十分脆弱,且缓苗后紧接着迎来夏季气温较高,光照较强,未及时实施浇水和遮阳等管理措施来减轻幼苗受到的胁迫程度。因此,适当遮阳处理有助于幼苗成活。另外,从4个试验点内各产地幼苗的保存率和生长量分析来看,保存率低的幼苗产地并没有因生长空间更大而更快,保存率高的产地幼苗也没有因空间竞争而生长缓慢。因此,保存率是幼苗适应气候的具体表现。在4个试验点中,山东试验点总体表现最好,该试验点在生长期内宽甸产地苗子死亡率最高,说明宽甸产地不适合引种保护,引种要尽量选取苗粗壮成活率高的汪清产地幼苗。
-
区域试验是基因型与环境互作研究中的热点,是对树种进行客观评价的理论基础及确定品种适应范围的重要依据[35]。过去有关林木稳定性分析多采用回归分析模型来评价,然而基因型与环境的互作效应是比较复杂且不是简单的线性叠加,所以分析效果不理想。AMMI模型把方差分析和主成分分析两者优点结合在同一模型,并且借助稳定性参数能定量直观对稳定性进行描述,从而提高产地幼苗稳定性分析的精确性。从图 2可知,在横轴方向上,汪清产地的苗高生长量均值最大,穆棱产地的苗高生长量均值最小;在纵轴方向上汪清、本溪产地的IPCA1值最大,说明这2个产地与环境间的互作效应最强,对环境条件较为敏感,稳定性弱,而和龙产地的IPCA1值最小,说明其受环境条件的影响小,稳定性较高。
-
(1) 东北红豆杉幼苗在山东试验点稳定较好且各生长性状生长量较大;在北京和吉林试验点的生长稳定性一般,但各生长性状生长量最小;在辽宁试验点的各生长性状生长量大但是稳定性最差。山东试验点适合种植东北红豆杉,建议可以往南试迁地保护。(2)由于保护目标不同,我们对于东北红豆杉的保护方向也有所侧重,应根据东北红豆杉产地特性制定相应的保护措施。对于生长快且稳定的产地,如和龙产地幼苗适合在4个试验点进行大面积迁地保护;对于生长快而相对不稳定性的保护主要侧重于其在生长适宜地区,如汪清产地幼苗可在山东烟台和辽宁本溪地区进行重点迁地保护;而对于生长量不高的产地保护则根据其生长适应能力强,对不良环境的抵御能力高等特点,如穆棱产地幼苗建议在环境条件相对恶劣的北京地区进行试迁地保护。而本溪产地幼苗生长较差且稳定性较差,适合在其正向互作效应的吉林生长,建议在吉林汪清迁地保护;宽甸产地幼苗在4个试验点都表现最差且保存率低不适合迁地保护。在引种时要尽量利用有利的互作,避免负互作效应带来的不利影响。
由于本试验只对4年生的东北红豆杉进行研究,存在环境变化造成的偶然性,同时试验点的选择不够密,覆盖范围不够广,为保证数据可靠,未来进行多年多点试验以减少误差显得尤为必要。
Growth stability analysis of ex situ conservation of Taxus cuspidata seedlings from different sources
-
摘要:
目的通过控制试验研究东北红豆杉宽甸、本溪、和龙、汪清、穆棱5个产地幼苗的生长特征,旨在探讨东北红豆杉在潜在分布区的引种适应性,为其迁地保护提供理论基础。 方法试验分别在山东、北京、辽宁、吉林4个试验点观测3年生苗,对其苗高、地径、冠幅等7个性状进行多地点联合分析和多重比较。另外,用AMMI模型对东北红豆杉苗高进行基因型(G)、环境(E)和基因型与环境(G×E)互作分析,比较不同产地东北红豆杉幼苗在不同生境下生长的稳定性。 结果(1) 山东试验点各产地幼苗地径、侧枝数、最长侧枝长和当年抽梢长生长量表现最好,分别高于北京试验点的3.77、1.88、10.35和1.30倍。各试验点内的表型性状变异系数平均高达30%以上,说明幼苗更具有很高的遗传分化。(2)从方差分量来看,苗高、地径、冠幅、当年抽梢长生长量在地点间的方差分量分别为62.89%、64.37%、32.65%、52.34%,远大于其他差异来源,说明地点效应对东北红豆杉生长影响最大。平均保存率由高到低分别为97.49%(吉林)>92.16%(辽宁)>89.92%(山东)>85.64%(北京)。 结论(1) 东北红豆杉生长性状除地径生长量在产地与试验地互作效应上差异不显著外,其他都具有极显著的产地、地点以及产地×地点的互作效应。山东、北京、辽宁试验点生长较好的是汪清产地幼苗;吉林试验点生长较好的是本溪产地幼苗。(2)和龙产地幼苗生长良好且稳定,适合在4个试验点进行大面积迁地保护;宽甸产地幼苗在4个试验点都表现最差且保存率低,不适合迁地保护。山东生长比北京好,表明可以对东北红豆杉进行往南迁地保护,前提必须做好产地幼苗筛选试验。 Abstract:ObjectiveThe seedling growth characteristics about five sources of Kuandian, Benxi, Helong, Wangqing and Muling of Taxus cuspidata were studied through the control experiment, aiming at exploring the adaptability of Taxus cuspidata in potential distribution areas, which will provide theoretical basis for ex situ conservation. MethodExperiments were conducted at four pilot sites in Shandong Province of eastern China, Beijing City, Liaoning Province and Jilin Province of northeastern China. Three-year-old seedlings were planted. Multi-site joint analysis and multiple comparisons of seven traits (including seedling height, basal diameter, and crown width, etc.) were conducted. In addition, the interactions of genotype (G), environment (E), and genotype with environment (G×E) were used to analyze the height of Taxus cuspidata in the AMMI model to compare the growth stability of Taxus cuspidata from different sources in different habitats. Result(1) The diameter of the seedlings, the number of lateral branch, the maximum length of lateral branch and the length of shoot in current year were the best in Shandong Province, which were 3.77 times, 1.88 times, 10.35 times and 1.30 times higher than the Beijing test site. The coefficient of variation of phenotypic traits averaged over 30% in each test site, indicating that seedlings had higher genetic differentiation. (2) From the aspect of variance, the variance components of height, basal diameter, crown width and length of shoot in current year were 62.89%, 64.37%, 32.65%, and 52.34%, which were much larger than the other sources. It is indicated that the location effect had the greatest impact on the growth of Taxus cuspidata. The average preservation rate was 97.49% (Jilin)>92.16% (Liaoning)>89.92% (Shandong)>85.64% (Beijing). Conclusion(1) The growth of basal diameter was not significantly different in the interaction between sources and experimental sites. The growth of other growth traits showed significant differences (P < 0.01) between sources and locations, as well as the interactions between sources and locations. Good performance of growth in Shandong Province, Beijing City, and Liaoning Province is the source of Wangqing, and the better growing in Jilin pilot is Benxi source. (2) The Helong source which grow well and stably is suitable for large-scale ex situ conservation in 4 pilot sites. The Kuandian sources have the worst performance in four pilot sites and the preservation rate is lower than other sources, which is not suitable for ex-situ conservation. The distribution area in Shandong Province is better than Beijing, indicating that the Taxus cuspidata can be moved to the south, and prerequisites must be done in the seed source screening test. -
图 1 不同试验点东北红豆杉幼苗生长性状差异的多重比较
Ⅰ.山东;Ⅱ.北京;Ⅲ.辽宁;Ⅳ.吉林。KD.宽甸产地;BX.本溪产地;HL.和龙产地;WQ.汪清产地;ML.穆棱产地。下同。不同小写字母表示同一试验点不同产地之间差异显著(P<0.05)。
Figure 1. Multiple comparisons for growth trait of Taxus cuspidata seedlings at different sites
Ⅰ. Shandong; Ⅱ. Beijing; Ⅲ. Liaoning; Ⅳ. Jilin. KD. Kuandian Source; BX. Benxi Source; HL. Helong Source; WQ. Wangqing Source; ML. Muling Source. The same below. Different letters indicated significant difference among different sources in the same test site (P < 0.05).
表 1 各试验点的地理气候条件
Table 1. Geographical and climatic conditions of the test sites
序号
No.试验地点
Test site纬度
Latitude经度
Longitude海拔
Altitude/m年均温
Annual average
temperature/℃年降水量
Annual
precipitation/mm年日照时数
Annual sunshine
hour/hⅠ 山东Shandong Province 37°39′N 121°03′E 51 13.4 627.6 2 656 Ⅱ 北京Beijing 40°00′N 116°20′E 49 12.5 612.0 2 720 Ⅲ 辽宁Liaoning Province 40°52′N 123°54′E 239 8.3 784.4 2 678 Ⅳ 吉林Jilin Province 43°12′N 129°37′E 168 4.5 549.4 2 357 表 2 东北红豆杉幼苗产地概况
Table 2. Overview of the origin of Taxus cuspidata seedlings
产地
Source纬度
Latitude经度
Longitude苗木来源
Seedling source辽宁宽甸
Kuandian County of Liaoning
Province40°03′N 124°20′E 辽宁省丹东市宽甸满族自治县鸭绿江生态苗木繁育基地
Ecological seedling breeding base of Yalu River in Kuandian Manchu Autonomous
County, Dandong City, Liaoning Province辽宁本溪
Benxi City of Liaoning Province41°19′N 123°47′E 辽宁省本溪市桓仁满族自治县桓仁宏创苗木基地
Huanren Hongchuang seedling breeding base of Huanren Manchu Autonomous
County, Benxi City, Liaoning Province吉林和龙
Helong City of Jilin Province42°31′N 129°00′E 吉林省和龙甑峰苗木有限公司
Helongzhenfeng Seedling Breeding Limited Company in Jilin Province吉林汪清
Wangqing County of Jilin Province43°18′N 129°47′E 吉林省延边朝鲜族自治州汪清县东北红豆杉生物科技有限公司
Taxus cuspidata Biotechnology Limited Company of Wangqing County, Yanbian
Korean Autonomous Prefecture, Jilin Province黑龙江穆棱
Muling City of Heilongjiang
Province44°36′N 130°12′E 黑龙江省穆棱市振才伟业红豆杉苗木基地
Zhencai Weiye Taxus cuspidata seedling breeding base in Muling City,
Heilongjiang Province表 3 不同试验点东北红豆杉幼苗生长性状的变异参数
Table 3. Variation parameters for growth traits of Taxus cuspidata seedlings at different sites
试验点
Test site生长性状
Growth trait均值
Mean标准差
Standard deviation变异范围
Range of variation变异系数
CV/%山东Shandong Province 苗高Seedling height/cm 4.14 1.21 2.50~6.80 29.15 地径Basal diameter/mm 2.60 0.85 1.56~5.56 32.56 冠幅Crown width/cm 1.93 0.77 0.70~3.50 39.95 侧枝数Number of lateral branch 1.60 1.23 0~4.00 76.94 最长侧枝长Maximum length of lateral branch/cm 2.38 0.86 1.10~4.30 36.32 当年抽梢长Length of shoot in current year/cm 6.07 1.86 3.30~10.20 30.63 北京Beijing 苗高Seedling height/cm 1.44 0.50 0.50~2.20 35.12 地径Basal diameter/mm 0.69 0.26 0.11~1.12 37.23 冠幅Crown width/cm 0.81 0.52 0.10~1.80 64.80 侧枝数Number of lateral branch 0.85 0.75 0~2.00 87.65 最长侧枝长Maximum length of lateral branch/cm 0.23 0.17 0~0.60 74.78 当年抽梢长Length of shoot in current year/cm 4.67 1.64 1.90~8.30 35.19 辽宁Liaoning Province 苗高Seedling height/cm 4.28 1.52 2.10~8.00 35.62 地径Basal diameter/mm 1.51 0.27 0.96~2.14 17.76 冠幅Crown width/cm 2.19 1.09 0.10~4.10 49.96 侧枝数Number of lateral branch 0.45 0.31 0~2.00 67.78 最长侧枝长Maximum length of lateral branch/cm 1.54 1.23 0.20~4.50 80.24 当年抽梢长Length of shoot in current year/cm 5.26 1.63 3.30~9.20 31.06 吉林Jilin Province 苗高Seedling height/cm 1.59 0.71 0.50~3.40 44.68 地径Basal diameter/mm 1.48 0.47 0.71~2.38 32.02 冠幅Crown width/cm 1.40 0.61 0.60~2.80 43.97 侧枝数Number of lateral branch 0.53 0.34 0~2.00 63.40 最长侧枝长Maximum length of lateral branch/cm 1.24 0.49 0.50~2.40 39.43 当年抽梢长Length of shoot in current year/cm 4.21 1.33 2.20~6.70 31.58 表 4 不同试验点东北红豆杉幼苗生长性状联合方差分析
Table 4. Joint analysis of variance for growth traits of Taxus cuspidata seedlings at different sites
生长性状
Growth trait地点Site 区组(地点)Block(site) 产地Source 产地×地点Source×site 随机误差
Random errordf F 方差分量
Variance
component/%df F 方差分量
Variance
component/%df F 方差分量
Variance
component/%df F 方差分量
Variance
component/%df 方差分量
Variance
component/%苗高Seedling height/cm 3 98.448** 62.89 3 2.457 1.57 4 15.329** 13.06 12 4.046** 10.34 57 12.14 地径Basal diameter/mm 3 67.622** 64.37 3 2.219 2.11 4 7.607** 9.65 12 1.518** 5.78 57 18.09 冠幅Crown width/cm 3 23.217** 32.65 3 0.415 0.58 4 6.735** 12.63 12 4.874** 27.42 57 26.72 侧枝数Number of lateral branch 3 13.304** 16.05 3 1.403 1.05 4 6.854** 58.17 12 3.900** 10.54 57 14.19 最长侧枝长Maximum length of
lateral branch/cm3 67.760** 22.76 3 1.508 2.40 4 10.848** 15.64 12 6.682** 26.69 57 32.51 当年抽梢长Length of shoot in
current year/cm3 21.493** 52.34 3 1.409 1.16 4 58.424** 11.17 12 3.527** 20.65 57 14.68 注:*差异显著,P<0.05;**差异极显著,P<0.01。Notes: * means significant difference,P<0.05;** means extremely significant difference,P<0.01. 表 5 不同试验点东北红豆杉幼苗缓苗率和保存率
Table 5. Prolonged seedling conservation and preservation rate of Taxus cuspidata seedlings at different sites
% 试验点
Test site产地
Source缓苗率Seedling conservation rate 保存率Preservation rate 均值±标准差
Mean± SD变异范围
Range of variation均值±标准差
Mean± SD变异范围
Range of variation山东Shandong Province KD 91.60±5.15 85.11~95.83 71.85±15.16 51.02~86.96 BX 100.00±0.00 100.00~100.00 91.40±3.07 88.64~95.65 HL 99.49±1.02 97.96~100.00 95.98±4.97 89.80~100.00 WQ 98.91±2.18 95.65~100.00 97.27±2.11 95.56~100.00 ML 100.00±0.00 100.00~100.00 93.09±5.77 84.78~97.87 均值Mean 98.00±4.01 85.11~100.00 89.92±11.75 51.02~100.00 北京Beijing KD 60.44±15.01 48.89~82.22 56.41±12.43 46.67~73.33 BX 86.51±4.26 82.50~92.50 82.94±4.13 76.74~85.00 HL 98.88±1.31 97.50~100.00 98.88±1.31 97.50~100.00 WQ 100.00±0.00 100.00~100.00 94.38±9.66 80.00~100.00 ML 98.13±2.39 95.00~100.00 95.63±5.15 90.00~100.00 均值Mean 88.79±16.62 48.89~100.00 85.64±17.37 46.67~100.00 辽宁Liaoning Province KD 90.49±5.08 85.42~95.83 78.40±6.32 72.92~87.50 BX 96.88±2.69 93.75~100.00 91.15±3.56 87.50~95.83 HL 99.48±1.04 97.92~100.00 98.43±2.00 95.83~100.00 WQ 100.00±0.00 100.00~100.00 100.00±0.00 100.00~100.00 ML 97.28±3.21 93.75~100.00 92.83±5.97 86.05~97.87 均值Mean 96.82±4.37 85.42~100.00 92.16±8.72 72.92~100.00 吉林Jilin Province KD 97.47±0.98 96.00~98.00 95.93±6.88 85.71~100.00 BX 95.78±7.24 85.00~100.00 96.52±4.72 90.00~100.00 HL 100.00±0.00 100.00~100.00 97.00±4.76 90.00~100.00 WQ 100.00±0.00 100.00~100.00 99.50±1.00 98.00~100.00 ML 99.51±0.98 98.04~100.00 98.49±1.98 95.83~100.00 均值Mean 98.55±3.39 85.00~100.00 97.49±4.14 85.71~100.00 表 6 东北红豆杉苗高性状的方差分析及AMMI模型分析
Table 6. ANOVA and AMMI analysis for the seeding height growth of Taxus cuspidata
分析方法
Analytical method变异来源
Source of variationdf 平方和
Sum of
square方差分量
Variance
component/%均方
Mean
squareF P 方差分析Analysis
of variance总的Total 79 232.034 2.937 地点内区组Location within the block 3 3.642 1.57 1.214 2.457 0.072 0 总处理Total treatment 19 200.226 10.538 19.879 0.000 1 产地Source 4 30.297 13.06 7.574 15.329 0.000 1 地点Site 3 145.937 62.89 48.646 98.448 0.000 1 产地×地点Source× site 12 23.992 10.34 1.999 4.046 0.000 3 误差Error 57 28.165 0.494 AMMI模型AMMI
analysisPCA1 6 21.798 90.85 3.633 28.506 0.000 1 PCA2 4 1.940 8.08 0.485 3.805 0.008 0 残差Residual 2 0.255 0.127 误差Error 60 31.808 0.530 表 7 东北红豆杉幼苗在显著互作效应轴上的得分及稳定性参数
Table 7. Scores of significant IPCA axis and stability parameter for Taxus cuspidata seedlings
项目
Item变量
Variable苗高生长量平均值
Mean growth of
seedling height/cm离差
DeviationIPCA1 IPCA2 Di 产地Source KD 2.531 -0.328 0.325 0.691 0.764 BX 2.456 -0.403 -0.962 -0.130 0.970 HL 3.069 0.210 0.136 0.013 0.136 WQ 3.969 1.110 1.012 -0.423 1.097 ML 2.269 -0.590 -0.511 -0.150 0.533 试验点Test place 山东Shandong Province 4.140 1.281 0.473 -0.665 0.816 北京Beijing 1.435 -1.424 -0.637 0.018 0.637 辽宁Liaoning Province 4.275 1.416 1.002 0.472 1.107 吉林Jilin Province 1.585 -1.274 -0.838 0.175 0.856 注:表中IPCA栏的数字表示各个主成分值。Note: data under IPCA in the table represent the principal component value. -
[1] 臧润国, 董鸣, 李俊清, 等.典型极小种群野生植物保护与恢复技术研究[J].生态学报, 2016, 36(22):7130-7135. http://d.old.wanfangdata.com.cn/Periodical/stxb201622022 Zang R G, Dong M, Li J Q, et al. Conservation and restoration for typical critically endangered wild plants with extremely small population[J]. Acta Ecologica Sinica, 2016, 36(22):7130-7135. http://d.old.wanfangdata.com.cn/Periodical/stxb201622022 [2] 周志强, 刘彤, 胡林林, 等.穆棱东北红豆杉年轮—气候关系及其濒危机制[J].生态学报, 2010, 30(9):2304-2310. http://d.old.wanfangdata.com.cn/Periodical/stxb201009009 Zhou Z Q, Liu T, Hu L L, et al. Tree ring-climate response and endangered mechanism deliberation of Japanese yew (Taxus cuspidata) in Muling Nature Reserve[J]. Acta Ecologica Sinica, 2010, 30(9):2304-2310. http://d.old.wanfangdata.com.cn/Periodical/stxb201009009 [3] 万基中, 王春晶, 韩士杰, 等.气候变化压力下建立东北红豆杉优先保护区的模拟规划[J].沈阳农业大学学报, 2014, 45(1):28-32. http://d.old.wanfangdata.com.cn/Periodical/synydxxb201401007 Wan J Z, Wang C J, Han S J, et al. The planning of priority protection area for Taxus cuspidata under climate change[J]. Journal of Shenyang Agricultural University, 2014, 45(1):28-32. http://d.old.wanfangdata.com.cn/Periodical/synydxxb201401007 [4] 吴雪琴, 徐刚标, 梁艳, 等.南岭地区观光木自然和人工迁地保护种群的遗传多样性[J].生物多样性, 2013, 21(1):71-79. http://d.old.wanfangdata.com.cn/Periodical/swdyx201301010 Wu X Q, Xu G B, Liang Y, et al. Genetic diversity of natural and planted populations of Tsoongiodendron odorum from the Nanling Mountains[J]. Biodiversity Science, 2013, 21(1):71-79. http://d.old.wanfangdata.com.cn/Periodical/swdyx201301010 [5] 任海.植物园与植物回归[J].生物多样性, 2017, 25(9):945-950. doi: 10.17520/biods.2017181 Ren H. The role of botanical gardens in reintroduction of plants[J]. Biodiversity Science, 2017, 25(9):945-950. doi: 10.17520/biods.2017181 [6] 黄宏文, 张征.中国植物引种栽培及迁地保护的现状与展望[J].生物多样性, 2012, 20(5): 559-571. http://d.old.wanfangdata.com.cn/Periodical/swdyx201205003 Huang H W, Zhang Z. Current status and prospects of ex situ cultivation and conservation of plants in China[J]. Biodiversity Science, 2012, 20(5): 559-571. http://d.old.wanfangdata.com.cn/Periodical/swdyx201205003 [7] 李思锋, 祁云枝, 张莹, 等.西安植物园国家重点保护野生植物资源引种研究[J].中国农学通报, 2009, 25(17):227-232. http://d.old.wanfangdata.com.cn/Periodical/zgnxtb200917048 Li S F, Qi Y Z, Zhang Y, et al. A brief study on the introduction of national key preserved wild plant resources in Xi'an Botanical Garden[J]. Chinese Agricultural Science Bulletin, 2009, 25(17):227-232. http://d.old.wanfangdata.com.cn/Periodical/zgnxtb200917048 [8] 姬少玲.东北红豆杉繁育技术研究进展[J].辽宁林业科技, 2015(4):61-63. doi: 10.3969/j.issn.1001-1714.2015.04.023 Ji S L. Research progress on the breeding techniques of Taxus cuspidate[J]. Liaoning Forestry Science and Technology, 2015(4):61-63. doi: 10.3969/j.issn.1001-1714.2015.04.023 [9] 刘彤, 李云灵, 周志强, 等.天然东北红豆杉(Taxus cuspidata)种内和种间竞争[J].生态学报, 2007, 27(3):924-929. doi: 10.3321/j.issn:1000-0933.2007.03.012 Liu T, Li Y L, Zhou Z Q, et al. Intraspecific and interspecific competition of Japanese yew (Taxus cuspidata)[J]. Acta Ecologica Sinica, 2007, 27(3):924-929. doi: 10.3321/j.issn:1000-0933.2007.03.012 [10] 常醉, 郭娜, 刘彤, 等.天然东北红豆杉中紫杉醇和三尖杉宁碱含量的分布及变化规律[J].中国中药杂志, 2011, 36(3):294-298. http://d.old.wanfangdata.com.cn/Periodical/zgzyzz201103015 Chang Z, Guo N, Liu T, et al. Distribution and variation of paclitaxel and cephalomannine contents in wild Taxus cuspidata[J]. China Journal of Chinese Material Medica, 2011, 36(3):294-298. http://d.old.wanfangdata.com.cn/Periodical/zgzyzz201103015 [11] Rodney C, Raymond E B, Ketchum R M L, et al. Taxol biosynthesis and molecular genetics[J]. Phytochemistry Reviews, 2006(5):75-97. doi: 10.1007/s11101-005-3748-2 [12] 刘彤, 崔海娇, 吴淑杰, 等.东北红豆杉幼苗光合和荧光特性对不同光照条件的响应[J].北京林业大学学报, 2013, 35(3):65-70. http://j.bjfu.edu.cn/article/id/9925 Liu T, Cui H J, Wu S J, et al. Response of photosynthetic and fluorescence characteristics of Japanese yew seedlings to different light conditions[J]. Journal of Beijing Forestry University, 2013, 35(3):65-70. http://j.bjfu.edu.cn/article/id/9925 [13] Damesin C. Respiration and photosynthesis characteristics of current-year stems of Fagus sylvatica: from the seasonal pattern to an annual balance[J]. New Phytologist, 2003, 158(3): 465-475. doi: 10.1046/j.1469-8137.2003.00756.x [14] 王小菲, 高文强, 刘建锋, 等.不同生境对栓皮栎幼苗光合生理特性的影响[J].生态学报, 2016, 36(24):8062-8070. http://d.old.wanfangdata.com.cn/Periodical/stxb201624018 Wang X F, Gao W Q, Liu J F, et al. Effects of habitat change on the photosynthetic characteristics of Quercus variabilis seedlings[J]. Acta Ecologica Sinica, 2016, 36(24):8062-8070. http://d.old.wanfangdata.com.cn/Periodical/stxb201624018 [15] 宋淑芳, 周冀衡, 李强, 等.基因与环境互作对云南保山烤烟主要潜香型物质的影响[J].应用生态报, 2014, 25(11):3223-3228. http://d.old.wanfangdata.com.cn/Periodical/yystxb201411021 Song S F, Zhou J H, Li Q, et al. Effects of interaction of genotypes with environments on major latently fragrant substances of flue-cured tobacco in Baoshan of Yunnan[J]. Southwest China. Chinese Journal of Applied Ecology, 2014, 25(11):3223-3228. http://d.old.wanfangdata.com.cn/Periodical/yystxb201411021 [16] 祝遵凌, 施曼, 周琦, 等.欧洲鹅耳枥区域化试验及适应性初步评价[J].浙江农林大学学报, 2016, 33(5):834-840. http://d.old.wanfangdata.com.cn/Periodical/zjlxyxb201605015 Zhu Z L, Shi M, Zhou Q, et al. Regional adaptability trials for Carpinus betulus in China[J]. Journal of Zhejiang A & F University, 2016, 33(5):834-840. http://d.old.wanfangdata.com.cn/Periodical/zjlxyxb201605015 [17] 刘宇, 徐焕文, 尚福强, 等.3个地点白桦种源试验生长稳定性分析[J].北京林业大学学报, 2016, 38(5):50-57. doi: 10.13332/j.1000-1522.20150463 Liu Y, Xu H W, Shang F Q, et al. Growth stability of Betula platyphylla sources from three sites[J]. Journal of Beijing Forestry University, 2016, 38(5):50-57. doi: 10.13332/j.1000-1522.20150463 [18] 郑桂萍, 蔡永盛, 赵洋, 等.利用AMMI模型进行寒地水稻品质分析[J].核农学报, 2015, 29(2):296-303. http://d.old.wanfangdata.com.cn/Periodical/hnxb201502012 Zheng G P, Cai Y S, Zhao Y, et al. Application of AMMI model in the analysis of rice quality at cold area[J]. Journal of Nuclear Agricultural Sciences, 2015, 29(2):296-303. http://d.old.wanfangdata.com.cn/Periodical/hnxb201502012 [19] 施万喜.利用AMMI模型分析陇东旱地冬小麦新品种(系)丰产稳产性[J].干旱地区农业研究, 2009, 27(3):37-43. http://d.old.wanfangdata.com.cn/Periodical/ghdqnyyj200903008 Shi W X. Analysis on high yield and stability of winter wheat cultivars in dryland of east Gansu by AMMI model[J]. Agricultural Research in the Arid Areas, 2009, 27(3):37-43. http://d.old.wanfangdata.com.cn/Periodical/ghdqnyyj200903008 [20] 刘博, 卫玲, 樊云茜, 等.基于AMMI模型的黄淮海夏大豆国家区试产量分析[J].中国农学通报, 2015, 31(27):69-74. doi: 10.11924/j.issn.1000-6850.casb15030053 Liu B, Wei L, Pan Y Q, et al. Yield analysis of national regional trial of soybean in Huang-Huai-Hai Area based on AMMI model[J]. Chinese Agricultural Science Bulletin, 2015, 31(27):69-74. doi: 10.11924/j.issn.1000-6850.casb15030053 [21] 乔转运, 王军辉, 顾万春, 等.桤木种源区域试验的AMMI模型分析[J].河南农业大学学报, 2005, 39(2):171-177. doi: 10.3969/j.issn.1000-2340.2005.02.012 Qiao Z Y, Wang J H, Gu W C, et al. AMMI model analysis on regional alder provenance trials of Alnas cremastogyn[J]. Journal of Henan Agricultural University, 2005, 39(2):171-177. doi: 10.3969/j.issn.1000-2340.2005.02.012 [22] 徐焕文, 刘宇, 李志新, 等.5年生白桦杂种子代多点稳定性分析及优良家系选择[J].北京林业大学学报, 2015, 37(12):24-31. doi: 10.13332/j.1000-1522.20140466 Xu H W, Li Y, Li Z X, et al.Analysis of the stability and superiority of five-year-old birch crossbreed families based on a multi-site test[J]. Journal of Beijing Forestry University, 2015, 37(12):24-31. doi: 10.13332/j.1000-1522.20140466 [23] 张华, 王金国.相同种源不同无性系东北红豆杉生长规律的比较研究[J].林业科技, 2011, 36(6):10-12. http://d.old.wanfangdata.com.cn/Periodical/lykj201106006 Zhang H, Wang J G. Comparison of growth rhythms to different Taxus cuspidata clones from same provenance[J]. Forestry Science & Technology, 2011, 36(6):10-12. http://d.old.wanfangdata.com.cn/Periodical/lykj201106006 [24] 吴榜华, 戚继忠.东北红豆杉植物地理学研究[J].应用与环境生物学报, 1995, 1(3):219-225. doi: 10.3321/j.issn:1006-687X.1995.03.002 Wu B H, Qi J Z. Study on phytogeography of Taxus cuspidata[J]. Chinese Journal of Applied and Environmental Biology, 1995, 1(3):219-225. doi: 10.3321/j.issn:1006-687X.1995.03.002 [25] 唐启义. DPS数据处理系统[M].北京:科学出版社, 2013. Tang Q Y.DPS data processing system[M]. Beijing: Science Press, 2013. [26] 吴为人.对基于AMMI模型的品种稳定性分析方法的一点改进[J].遗传, 2000, 22 (1): 31-32. http://d.old.wanfangdata.com.cn/Periodical/yc200001010 Wu W R. An improvement on the method of variety stability analysis based on the AMMI model[J]. Hereditas, 2000, 22(1): 31-32. http://d.old.wanfangdata.com.cn/Periodical/yc200001010 [27] 李志新, 赵曦阳, 杨成君, 等.转TaLEA基因小黑杨株系变异及生长稳定性分析[J].北京林业大学学报, 2013, 35(2):57-62. http://j.bjfu.edu.cn/article/id/9878 Li Z X, Zhao X Y, Yang C J, et al. Variation and growth adaptability analysis of transgenic Populus simonii×P. nigra lines carrying TaLEA gene[J]. Journal of Beijing Forestry University, 2013, 35(2):57-62. http://j.bjfu.edu.cn/article/id/9878 [28] 肖遥, 张蕊, 楚秀丽, 等.24个产地南方红豆杉在两试验点的生长差异及其选择[J].林业科学研究, 2017, 30(2):342-348. http://d.old.wanfangdata.com.cn/Periodical/lykxyj201702022 Xiao Y, Zhang R, Chu X L, et al. Growth variation and selection of Taxus wallichiana var. mairei from 24 locations in two trial plots[J]. Forest Research, 2017, 30(2):342-348. http://d.old.wanfangdata.com.cn/Periodical/lykxyj201702022 [29] 肖遥, 楚秀丽, 徐肇友, 等.南方红豆杉2年生容器苗多点试验的生长节律家系变异[J].林业科学研究, 2016, 29(2):238-244. doi: 10.3969/j.issn.1001-1498.2016.02.014 Xiao Y, Chu X L, Xu Z Y, et al. Family variation on growth rhythm of two-year-old container seedlings of Taxus wallichiana var. mairei in multi-location trials[J]. Forest Research, 2016, 29(2):238-244. doi: 10.3969/j.issn.1001-1498.2016.02.014 [30] 曾志光, 肖复明, 包国华, 等.木荷地理种源苗期性状遗传变异研究[J].林业科学研究, 2005, 18(1):27-30. doi: 10.3321/j.issn:1001-1498.2005.01.005 Zeng Z G, Xiao F M, Bao G H, et al. A study on variations of seedling traits in Schima superba geographic provenances[J]. Forest Research, 2005, 18(1):27-30. doi: 10.3321/j.issn:1001-1498.2005.01.005 [31] 王艺, 张蕊, 冯建国, 等.不同种源南方红豆杉生长差异分析及早期速生优良产地筛选[J].植物资源与环境学报, 2012, 21(4):41-47. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zwzyyhjxb201204006 Wang Y, Zhang R, Feng J G, et al. Analysis on growth difference of Taxus wallichiana var. mairei from different sources and selection of fast-growing sources at early age[J]. Journal of Plant Resources and Environment, 2012, 21(4):41-47. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zwzyyhjxb201204006 [32] 周蒙, 刘文耀, 马文章, 等.不同地理种源紫茎泽兰的生态适应性比较[J].应用生态学报, 2009, 20(7):1643-1649. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yystxb200907019 Zhou M, Liu W Y, Ma W Z, et al. Ecological adaptability of different provenance Eupatorium adenophorum Spreng[J]. Chinese Journal of Applied Ecology, 2009, 20(7):1643-1649. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yystxb200907019 [33] 贺善安, 李新华, 彭峰, 等.南方红豆杉迁地保护小种群动态的研究[J].植物资源与环境学报, 2007, 16(1):35-39. doi: 10.3969/j.issn.1674-7895.2007.01.008 He S A, Li X H, Peng F, et al. Study on dynamic of small ex-situ conservation population of Taxus chinensis var. mairei[J]. Journal of Plant Resources and Environment, 2007, 16(1):35-39. doi: 10.3969/j.issn.1674-7895.2007.01.008 [34] Nagamitsu T, Nagasak A K, Yoshimaru H, et al. Provenance tests for survival and growth of 50-year-old Japanese larch (Larix kaempferi) trees related to climatic conditions in central Japan[J]. Tree Genetics & Genomes, 2014, 10(1):87-99. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=a4f051a8dc88e25359e2b4c365bde397 -