-
水青树(Tetracentron sinense Oliv.)隶属水青树科(Tetracentraceae)水青树属(Tetracentron),是水青树科仅有的一种高大落叶乔木[1];主要分布在我国西南及华中地区的900~3 500 m的深山、峡谷、陡坡悬崖等处;其树形伸展大气,叶形美观,颜色翠绿,可作为供观赏的行道树;其茎干中含有HIV和有关白细胞活性的化学成分[2-3],在医学研究方面有着重要的价值;其化石发现于中生代中期的始新世地层,是现存的较古老被子植物,其对于研究古代植物区系以及被子植物的系统演化和起源等方面具有重要的科学价值[4]。由于受自然条件的变化和人为因素干扰等的影响,水青树的种群和数量日益缩减,目前野外分布零散且稀少,已是国家Ⅱ级重点保护濒危植物[5]。国内外科学家已开始密切关注水青树的种质资源怎样才能得到有效保护和利用。
对于濒危植物的保护关键之处在于如何较好的保护其遗传多样性, 而人类保护珍稀濒危植物重要的措施之一就是对其进行迁地保护[6]。虽然水青树分布范围比较广泛,且种质资源数量丰富,但是如何对其进行有效的评价、研究和利用是水青树种质资源保护急需解决的难题[7]。同时,在对水青树进行迁地保护时,选择怎样的取样策略才能保证用尽可能少的种质保存尽可能多的现有遗传资源,这是在具体保护措施制定时面临的又一难题。“核心种质”的出现成为解决这一问题的新思路[8-10]。Frankel提出, 运用一定方法建立植物的核心种质,把最小的资源数量以及最大程度的遗传重复从原种质中筛选出来从而代表整个遗传资源的多样性[8]。因此,核心种质资源库在植物遗传资源的利用以及管理方面起了很大的作用[11];而在构建其核心种质资源库时,不同的取样策略对于不同种类植物的效果是不一样的,选择有效的取样策略就得尤为重要[12]。目前,刘娟等运用位点优先取样策略构建了野杏的种质资源库[13];赵宁等对山茱萸的种质资源库构建时提出采用位点优先取样策略是构建该物种的最适方法[14];胡建斌等运用比例法构建甜瓜的核心库[15]。我国陆续对一些其它木本植物进行了核心种质资源库的构建,如新疆野核桃[16]、新疆野苹果[17]等, 为它们的种质资源保护奠定了基础。然而,有关濒危木本植物种质库构建方面的研究还很少[18]。
经全面考察水青树的分布现状发现,水青树多呈零星散生,通常在400 m2的范围内仅有2~3棵成年植株。结合水青树实际分布情况,笔者于2015—2016年在四川、重庆、贵州、云南、陕西、甘肃、湖南、湖北这8个省市挑选出26个水青树居群共161个株系(图 1),构建了水青树的种质资源库。但是构建的种质资源库的多样性如何,且应该运用什么方法、什么材料来进一步构建水青树的核心种质,迄今为止,还没有相关研究。由于濒危植物水青树种子千粒质量很轻,它的平均值低于很多木本植物[4],野外收集有一定难度。相较于植物其它器官而言,叶对环境变化的响应比较明显,即环境一旦发生变化,叶的性状会有较大几率发生变化。换句话说,叶形态的变异也能在一定程度上反映其与环境之间的相互作用[19-20]。目前,通过叶表型变异来探讨植物种群的遗传变异已有较多研究[21-25]。因此,本文对161个水青树种质的叶表型变异情况进行评价,并利用多种方法比较分析了其核心种质资源库构建的策略,为水青树种质资源的进一步保护提供科学依据。
HTML
-
利用软件SPSS计算构建的161个种质的水青树表型性状各指标值见表 1。各表型性状在一定程度上存在差异。由于各表型性状存在量纲差别,而变异系数正好消除了各性状之间的量纲区别,因此,变异系数具有可比性和代表性。根据表 1结果可知:叶先端尖削指数的变异系数最大,为50%,其它性状的变异系数在10%~50%间(上表皮细胞弯曲度指数除外,它的变异系数最小)。整体而言,本实验所选的水青树叶表型性状都存在较大的变异,其遗传多样性较高,证明这161个种质资源能较好的反应水青树的遗传变异。
变异系数
variable coefficient标准差
standard deviation平均值
mean方差
variance极差
range上表皮细胞数目Number of epidermal cells/个 0.43 52.53 122.92 2 758.94 355.68 气孔器数目Number of stmata/个 0.39 4.95 12.82 24.54 23.91 气孔器长The stomata length/μm 0.16 7.04 42.78 49.57 29.89 气孔器宽The stomata width/μm 0.18 6.54 35.64 42.81 25.25 上表皮细胞面积Upper epidermal cell area/μm2 0.40 734.99 1 835.44 540 208.74 3 181.18 上表皮细胞周长Perimeter of epidermal cells/μm 0.27 51.77 192.99 2 680.64 248.67 上表皮细胞弯曲度指数Upper epidermal cell curvature index 0.08 0.11 1.29 0.01 0.53 上表皮细胞密度Density of epidermal cells/(个·mm-2) 0.42 374.07 889.40 139 930.13 2 370.14 气孔器密度Stomatal density/(个·mm-2) 0.39 35.90 92.92 1 288.96 173.26 叶面积Leaf area/cm2 0.25 12.72 51.03 161.75 73.50 叶周长Leaf perimeter/cm 0.12 3.59 28.88 12.91 21.37 叶柄长Petiole length/cm 0.18 0.47 2.53 0.22 2.44 叶长Length of leaf/cm 0.14 1.39 9.80 1.93 12.37 叶最长The longest leaf/cm 0.14 1.51 11.10 2.30 8.64 叶宽Leaf width/cm 0.18 1.30 7.12 1.69 11.79 叶基Phyllo podium/cm 0.14 0.51 3.61 0.26 3.69 叶基尖削指数Leaf base sharpening index 0.20 0.27 1.34 0.07 2.84 叶先端尖削指数Leaf apex sharpening index 0.50 0.32 0.65 0.10 4.27 叶全长Leaf length/cm 0.12 1.54 12.35 2.36 11.46 叶形指数Leaf index 0.34 0.49 1.46 0.24 6.06 Table 1. The phenotypic traits values of T. sinense
-
分析水青树20个叶表型性状的主成分,获得了5个能够代表其性状的特征值(分析时,从每个原始变量里提取的信息,SPSS默认保留其>1的值)大于1的主成分,且其累积贡献率达到了87.692%(表 2,自第6主成分开始特征值小于1)。对以上所提取出的主成分分别计算各主成分值,并进一步利用主成分值对各个体进行聚类分析,以此来研究各个体间的亲缘关系。
特征值
Eigenvalue贡献率/%
Contribution rate累积贡献率/%
Accumulating contribution rate1 7.598 33.036 33.036 2 5.959 25.910 58.946 3 3.614 15.715 74.661 4 1.747 7.595 82.257 5 1.250 5.436 87.692 6 0.652 2.836 90.528 Table 2. Principal Component Analysis (PCA)
-
采用随机取样法抽取10个比例的种质资源。随着抽样比例的增大,除CR外,叶表型性状的各参数均呈先逐渐增大后逐渐减小的趋势;其中,VR和VD在取样比例为35%时达到最大,此时CR和MD也相对较高(表 3)。因此,R6最能代表原种质,并得到56份核心种质。
种质资源库
Germplasm resources bank取样比例
Sampling ratio/%均值差异百分率
MD/%方差差异百分率
VD/%变异系数变化率
CR/%极差符合率
VR/%R1 10 1 21 50 85 R2 15 2 17 54 86 R3 20 2 16 58 87 R4 25 2 16 61 87 R5 30 0 19 64 86 R6 35 1 34 89 113 R7 40 0 26 89 110 R8 45 0 24 91 110 R9 50 0 19 91 108 R10 55 0 13 91 106 Table 3. Percentage difference between germplasm resources and original germplasm characters by random sampling strategy
-
用位点优先取样法来抽取10个比例的种质资源时发现,各取样比例的MD均较小,绝大部分为0;VD和CR均随取样比例的增大而呈现出相似的变化趋势,在取样比例为45%时达到最大(表 4),而此时VR的值相对较高。因此,S8最能反映原种质各情况,并得到72份核心种质。
种质资源库
Germplasm resources bank取样比例
Sampling ratio/%均值差异百分率
MD/%方差差异百分率
VD/%变异系数变化率
CR/%极差符合率
VR/%S1 10 1 73 82 47 S2 15 21 136 87 58 S3 20 1 80 88 58 S4 25 0 74 83 60 S5 30 0 71 82 60 S6 35 0 71 82 62 S7 40 0 71 82 66 S8 45 0 111 103 89 S9 50 0 109 103 92 S10 55 1 106 102 92 Table 4. Percentage difference between the germplasm resources and the original germplasm characters sampled by allele preferred sampling strategy
-
用偏离度取样策略构建水青树种质资源库时发现:MD的值大多接近0%;CR和VR值较大,但是VD值较小。从结果中可以得到,该种取样策略下,筛选出的最适种质资源应是45%取样比例下的种质(表 5),即,D8为最适种质资源,共72份。
种质资源库
Germplasm resources bank取样比例
Sampling ratio/%均值差异百分率
MD/%方差差异百分率
VD/%变异系数变化率
CR/%极差符合率
VR/%D1 10 2 20 63 107 D2 15 2 7 66 101 D3 20 1 7 67 94 D4 25 0 10 68 92 D5 30 0 17 69 88 D6 35 0 16 70 88 D7 40 1 12 71 91 D8 45 0 22 91 109 D9 50 0 13 91 106 D10 55 0 13 93 106 Table 5. Percentage difference between germplasm resources and original germplasm characters sampled by deviation sampling
-
比较不同策略构建的核心种质(表 6),位点优先取样策略构建得来的种质,其各表型性状的极差、均值和方差等值均大于随机和偏离度这2种策略;不仅如此,叶表型性状的变异系数在位点优先取样时最小,随机取样时却是最大的;MD在位点优先和偏离度取样时均为0;位点优先取样的VD远大于随机以及偏离度取样策略的VD;虽然位点优先取样法下的VR略小于其他2种方法,但综合考虑,位点优先取样策略更适合用来构建濒危植物水青树的核心种质资源库。因此,在相同的聚类方法下,结合位点优先取样策略一共构建了72份水青树核心种质资源。
核心种质资源库
Core collection极差
Range均值
Mean标准差
Standard deviation方差
Variance变异系数
Coefficient of variation均值差异百分率MD/% 方差差异百分率VD/% 变异系数变化率CR/% 极差符合率VR/% D8 251.59 147.51 53.93 27 856.49 0.27 0 22 91 109 R6 218.79 148.69 57.34 32 664.20 0.28 1 34 89 113 S8 260.66 149.40 52.74 28 184.87 0.25 0 111 103 89 Table 6. The values of each parameter under the three sampling strategies
-
由表 7可知:对初步构建的核心种质进行t检验,结果均不显著,这说明本次实验所构建的核心种质与原种质在遗传多样性上相差不大。检测种质各表型性状平均值、变异系数与原种质符合度发现,均值符合率达到94.4%以上;变异系数符合率也高达84.8%以上。因此,本次构建的核心种质资源库具有原种质中的各性状变异,且这72份种质是有效的。
编号
Code性状
Traits均值Mean 变异系数Coefficient of variation t值
t-value原种质
Primary germplasm核心种质
Core collection符合率
Coincidence rate原种质
Primary germplasm核心种质
Core collection符合率
Coincidence rate1 上表皮细胞数目Number of epidermal cells/个 122.92 116.05 0.94 0.43 0.36 0.85 0.98 2 气孔器数目Number of stomata/个 12.82 12.15 0.95 0.39 0.34 0.89 1.01 3 气孔器长The stomata length/μm 42.78 43.43 1.02 0.16 0.16 0.95 -0.65 4 气孔器宽The stomata width/μm 35.64 36.12 1.01 0.18 0.18 0.96 -0.52 5 上表皮细胞面积Upper epidermal cell area/μm2 1 835.44 1 953.04 1.06 0.40 0.38 0.95 -1.13 6 上表皮细胞周长Perimeter of epidermal cells/μm 192.99 200.12 1.04 0.27 0.25 0.93 -0.98 7 上表皮细胞弯曲度指数
Upper epidermal cell curvature index1.29 1.29 1.00 0.08 0.07 0.88 -0.42 8 上表皮细胞密度Density of epidermal cells/(个·mm-2) 889.40 840.96 0.95 0.42 0.36 0.86 0.97 9 气孔器密度Stomatal density/(个·mm-2) 92.92 88.03 0.95 0.39 0.34 0.89 1.01 10 叶面积Leaf area/cm2 51.03 51.50 1.01 0.25 0.24 0.95 -0.27 11 叶周长Leaf perimeter/cm 28.88 28.97 1.00 0.12 0.12 0.96 -0.17 12 叶柄长Petiole length/cm 2.53 2.53 1.00 0.18 0.18 0.96 0.08 13 叶长Length of leaf/cm 9.80 9.73 0.99 0.14 0.17 1.18 0.33 14 叶最长The longest leaf/cm 11.10 11.17 1.01 0.14 0.13 0.98 -0.35 15 叶宽Leaf width/cm 7.12 7.27 1.02 0.18 0.21 1.13 -0.78 16 叶基Phyllo podium/cm 3.61 3.65 1.01 0.14 0.14 1.00 -0.54 17 叶基尖削指数Leaf base sharpening index 1.34 1.35 1.01 0.20 0.25 1.23 -0.38 18 叶先端尖削指数Leaf apex sharpening index 0.65 0.68 1.05 0.50 0.70 1.40 -0.61 19 叶全长Leaf length/cm 12.35 12.29 1.00 0.12 0.14 1.09 0.26 20 叶形指数Leaf index 1.46 1.49 1.02 0.34 0.47 1.40 -0.39 注:t0.05=1.960,t0.01=2.576。 Table 7. To compare the characters of the core collection and the primary germplasm
-
利用主成分分析法确认最终得到的72份核心种质。表 8显示了各主成分解释原种质中原始变量总方差的情况。核心种质在第4个主成分之后特征值开始小于1,其4个主成分几乎包含了原种质各个原始变量80%以上的遗传信息,因此,保留4个主成分最为合适。构建的核心种质中,4个主成分的累积贡献率达到83.366%高于原种质的82.257%,可见,核心种质的贡献率高于原种质,这样有效的避免了种质的冗余。
主成分
Principal component原种质The original germplasm 核心种质The core collection 特征值
Eigenvalue贡献率/%
Contribution rate累积贡献率/%
Accumulating contribution rate特征值
Eigenvalue贡献率/%
Contribution rate累积贡献率/%
Accumulating contribution rate1 7.598 33.036 33.036 7.214 31.365 31.365 2 5.959 25.910 58.946 5.752 25.008 56.373 3 3.614 15.700 74.661 4.355 18.935 75.308 4 1.747 7.596 82.257 1.853 8.058 83.366 5 1.250 5.435 87.692 0.766 3.832 87.435 6 0.652 2.836 90.528 Table 8. 8 Principal component analysis of the germplasm and the core collection resources in the germplasm bank of T. sinense