Phenotypic variation and diversity of natural Pinus koraiensis populations in Jilin Province of northern China
-
摘要:目的 为揭示天然红松居群的表型分化程度及变异模式，以吉林省的6个天然红松居群为研究对象，方法 采用方差分析、主成分分析、聚类分析及多性状综合评价等方法对红松居群的13个表型性状（针叶性状和种实性状）进行系统分析和综合评价。结果 （1）除针叶束粗/针叶厚在居群间差异未达显著水平外，其余性状在居群内和居群间均达极显著差异水平（ P < 0.01）。（2）在居群间红松的平均表型分化系数为59.33%，其变异大于群体内（40.67%）。（3）6个居群的平均表型变异系数为11.30%，针叶性状和种实性状的表型变异系数依次为针叶性状（14.56%） > 球果性状（10.48%） > 种子性状（5.87%）。（4）主成分分析结果表明，红松居群表型多样性基本来源为针叶性状 > 种子性状 > 球果性状。（5）利用欧氏距离对红松居群进行聚类分析，将红松6个居群划分为3大类群，分别为P1与P6、P3与P5、P2与P4。（6）利用多性状综合评价法，分别以种实性状（种长、种宽、百粒质量和球果质量）和针叶性状（针叶长、针叶宽、针叶厚、针叶束粗）对天然红松居群进行评价，分别筛选出1个居群。结论 红松居群具有中等的表型多样性，居群间和居群内均存在丰富的表型变异，研究结果为今后红松种质资源的保护和利用提供基础，为红松优良居群的构建提供材料。Abstract:Objective In order to reveal the phenotypic differentiation and variation pattern of natural Pinus koraiensis populations, 6 natural populations in Jilin Province of northeastern China were taken as the research object.Method The 13 phenotypic traits (needles traits and fruit traits) of the Pinus koraiensis population were systematically analyzed and comprehensively assessed by means of variance analysis, principal component analysis, cluster analysis and comprehensive evaluation of multiple traits.Result (1) Except that there was no significant difference in the ratio of fascicle width to needle thick among populations, the other traits reached extremely significant differences within and among populations (P < 0.01). (2) The average phenotypic differentiation coefficient of Pinus koraiensis among the populations was 59.33%, and its variation was much smaller than within the population (40.67%). (3) The average phenotypic variation coefficient of the 6 populations was 11.30%, the phenotypic variation coefficient decreased from needle traits (14.56%) to cone traits (10.48%) and seed traits (5.87%). (4) The results of principal component analysis indicated that, the basic source of phenotypic diversity of Pinus koraiensis population was needle traits > seed traits > cone traits. (5) The results of cluster analysis showed that the 6 populations of Pinus koraiensis were divided into 3 groups. (6) Using comprehensive evaluation of multiple traits method, the natural
populations were evaluated by fruit traits (seed length, seed diameter, hundred seed mass, cone mass) and needle traits (needle length, needle diameter, needle thick, fascicle width), and 1 population was screened out, respectively. Conclusion Pinus koraiensis populations have moderate phenotypic diversity, and there are abundant phenotypic variations between and within groups. The research results can provide a basis for the protection and utilization of Pinus koraiensis germplasm resources, and it can provide materials for the construction of superior populations of Pinus koraiensis.-
Keywords:
- Pinus koraiensis /
- population /
- phenotypic trait /
- phenotypic differentiation
-
云南金钱槭(Dipteronia dyerana)是隶属于槭树科(Aceraceae)金钱槭属的一种落叶乔木，仅分布于云南东南部及贵州西南部海拔2 000~2 500 m的疏林中，为国家Ⅱ级珍惜濒危保护植物[1]。长期以来，云南金钱槭的种级系统位置界定，主要依据其独特的形态学特征(翅果成圆环形等)，并与金钱槭(Dipteronia sinensis)共同组成我国特有属——金钱槭属。
金钱槭属作为槭属(Acer)的姐妹属，常作为槭属等相关研究的外类群，但它常常不能形成单系群或属内物种嵌于槭属中。Tian等[2]采用ITS和trnL-F片段对槭属进行系统学研究中，云南金钱槭内嵌于槭属，导致金钱槭属没有形成单系群；Pfosser等[3]在关于亚洲槭属物种的亲缘关系重建中，金钱槭内嵌于槭属之内，导致金钱槭属在分析中不能作为外群；Grimm等[4]使用ITS片段，利用最大似然法构建槭属系统发育树时，云南金钱槭出现在槭属内部；Li等[5]研究发现，在最大简约法(MP)和最大似然法(ML)构建系统树时，金钱槭属的物种不能形成单系群。此外，相对早期的研究也存在类似情况[6-7]。Yang等[8]虽对此进行了讨论，但目标基因片段选择相对较少，且对于槭属属内分类组的辅助研究材料选择较为局限。因此，对云南金钱槭在其近缘物种中进行系统位置的确定研究，具有重要的理论与实践意义。
近年来，在物种系统位置相关分析中，基因片段的选择多集中于核ITS片段以及叶绿体trnL-trnF、rbcL、rpl16等片段[9-15]。此外，在金钱槭属及槭属的研究中，也常使用叶绿体psbM-trnD、trnD-trnT片段[5]以及psbA-trnH片段等(未发表数据)。因此，本研究选取8个基因片段，包括6个叶绿体基因片段(psbM-trnD、rbcL、trnD-trnT、rpl16、trnL-trnF和psbA-trnH)与2个核基因片段(ITS与CHS)，研究金钱槭属及其近缘属(槭属)的系统发育关系。通过多基因组合系统学分析，尝试理清云南金钱槭的系统位置及其与近缘类群的系统关系，并为其姐妹属(例如：槭属、七叶树属、掌叶木属等)的研究提供帮助。
1. 材料与方法
1.1 实验材料
本研究以云南金钱槭为材料，并在槭树属内14个分类组中[16]，每组各选取一个物种作为辅助材料(其中psbM-trnD、trnD-trnT与CHS选取不足14个)(表 1)。此外，七叶树属(Aesculus)作外群。叶绿体基因psbM-trnD、rbcL、trnD-trnT片段与核基因ITS片段序列全部选自GenBank，叶绿体基因rpl16、trnL-trnF与psbA-trnH片段大部分为本实验所得，核基因CHS片段全部为本实验所得(表 1)。每个实验物种采集其幼嫩叶片，分别放入置有变色硅胶的取样袋中干燥、备用。同时，压制凭证标本经由北京林业大学博士生高健鉴定并存放于北京林业大学林学院。
表 1 试验材料及其序列号Table 1. Species used in this study and their GenBank accession numbers属
Genera组
Section来源
Source编号
Code采集者
Collector种
Species序列登录号GenBank accession No. ITS CHS psbM- trnD rbcL trnD- trnT rpl16 trnH-PsbA trnL-F 七叶树属
Aesculus青秀山，南宁Qingxiu Mountain, Nanning k028 陆志成
LU Zhi-chengA. wangii AF406968 KX000286* NCBI A. glabra DQ659840 AY968623 NCBI A. flava DQ978441 AY968613 DQ978509 DQ978642 DQ978573 NCBI A. parviflora DQ978448 DQ978510 DQ978643 DQ978574 上海植物园，上海Shanghai Botanical Garden, Shanghai i008 杜宝明
DU Bao-mingA. chinensis KX000285* 金钱槭属
Diteronia堵河源保护区，十堰
Duheyuan National Nature Reserve, Shiyany010 郑德国
ZHENG De-guoD. dyerana DQ238333 KX000287* DQ659838 DQ978443 DQ659778 DQ978512 DQ978645 DQ978576 青秀山，南宁Qingxiu
Mountain, Nanningk029b 陆志成
LU Zhi-chengD. sinensis AY605290 KX000288* DQ659839 DQ978444 DQ659779 KU500460* KU500519* KU522513* 槭属
Acer尖齿枫组
sect. Arguta安图县，延边Antu
County, Yanbianq004 刘琪璟
LIU Qi-jingA. barbinerve AJ634569 KX000274* DQ978395 KU500434* KU500513* KU522490* NCBI A. acuminatum DQ659841 DQ659781 全缘叶枫组
sect. Oblonga上海植物园，上海
Shanghai Botanical Garden, Shanghaii002 杜宝明
DU Bao-mingA. buergerianum AY605466 KX000275* DQ659843 DQ659783 KU500435* KU500514* KU522491* NCBI A. buergerianum var. formosanum DQ978396 枫组sect. Acer NCBI A. caesium subsp. giraldii AY605295 DQ659893 DQ978397 DQ659833 DQ978453 DQ978584 DQ978519 穗状枫组
sect. SpicataNCBI A. caudatum subsp. multiserratum AY605432 DQ978401 华南植物园，广州South
China Botanical Garden, Guangzhoud005 简曙光
lJIAN Shu-guangA. caudatum KX000276* DQ659887 DQ659827 KU500444* KU500509* KU522508* 大花枫组
sect. MacranthaNCBI A. davidii subsp. davidii AY605392 DQ978406 浙江农林大学，临安
Zhejiang Agriculture and Forestry University, Lin’ano001 陈秋夏
CHEN Qiu-xiaA. davidii KX000277* DQ659848 DQ659788 浙江农林大学，临安
Zhejiang Agriculture and Forestry University, Lin'ano006 陈秋夏
CHEN Qiu-xiaA. davidii subsp. grosseri KU500459* KU500550* KU522512* 鸡爪枫组
sect. PalmataNCBI A. palmatum subsp. palmatum AY605425 中科院植物所，北京
Institute of Botany, Chinese Academy of Sciences, Beijingb015 廖培钧
LIAO Pei-junA. palmatum KX000280* DQ659867 DQ978421 DQ659807 KU500472* KU500503* KU522530* 三小叶枫组
sect. Trifoliata堵河源保护区，十堰
Duheyuan National Nature Reserve, Shiyany005 郑德国
ZHENG De-guoA. griseum AY605467 KX000278* DQ659857 DQ978411 DQ659797 KU500463* KU500530* KU522514* 扁果枫组
sect. HyptiocarpaNCBI A. laurinum AM113541 DQ978413 DQ659794 DQ978473 DQ978604 DQ978536 复叶枫组
sect. Negundo华南植物园，广州South
China Botanical Garden, Guangzhoud001a 简曙光
JIAN Shu-guangA. negundo AY605406 KX000279* DQ659864 DQ978417 DQ659804 KU500471* KU500535* KU522527* 五小叶枫组
sect. Pentaphtlla茂县，阿坝藏族羌族自治州Maoxian County, Aba Tibetan and Qiang Autonomous Prefecture u003 包维楷
BAO Wei-kaiA. pentaphyllum DQ238478 KX000281* DQ659870 DQ978422 DQ659810 KU500474* KU500538* KU522534* 疏毛枫组sect.
Pubescentia中科院植物所，北京Institute of Botany, Chinese Academy of Sciences, Beijing b013 廖培钧
LIAO Pei-junA. pilosum DQ238345 DQ978423 KU500491* KU500546* KU522535* NCBI A. pilosum var. stenolobum KX000282* 桐状枫组sect.
Platanoidea黑龙江森林植物园，哈尔滨Heilongjiang Forest Botanical Garden, Harbin h004 毛子军
MAO Zi-junA. pictum subsp. mono LK022665 KU500431* KU500534* KU522525* NCBI A. pictum subsp. pictum AB872557 茶条枫组
sect. Ginnala黑龙江森林植物园，哈尔滨Heilongjiang Forest Botanical Garden, Harbin h002 毛子军
MAO Zi-junA. tataricum subsp. ginnala AY605363 KX000284* DQ659855 DQ978436 DQ659795 KU500479* KU500521* KU522544* NCBI A. sterculiaceum subsp. franchetii DQ366145 DQ978435 坚果枫组
sect. Lithocarpa紫金山，南京Zijin Mountain, Nanjing m008 葛之葳
GE Zhi-weiA. sinopurpurascens KX000283* DQ659878 DQ659818 KU500490* KU500545* KU522538* 注:*代表本研究中获得的序列。Notes: * represents sequences obtained in present study; NCBI, National Center for Biotechnology Information. 1.2 DNA提取、聚合酶链式反应(PCR)及测序
每份样品取约0.2 g干燥叶片，采用植物基因组DNA小型提取试剂盒(Plant Genomic DNA kit)基因组DNA。叶绿体基因组片段rpl16、trnL-trnF、psbA-trnH和核基因组片段CHS的PCR扩增均在Labnet MultiGeneTM 96-well Gradient Thermal Cycler上完成。PCR反应体系为20 μL，主要包含DNA模板10 ~ 30 ng、50 mmol/L Tris-HCl、1.5 mmol/L MgCl2、1 mmol/L dNTPs、各0.3 μL的正反引物、0.2 U Taq DNA聚合酶(博尔纳德股份有限公司，台湾)。PCR反应程序为：94 ℃预变性3 min；94 ℃变性30 s，52~58 ℃退火30 s，72 ℃延伸60 s，共30个循环；最后72 ℃延伸10 min。PCR扩增产物经1%琼脂糖凝胶电泳检测后，送往台湾基龙米克斯生物科技股份有限公司并在ABI3730XL遗传分析仪上进行测序。
1.3 数据分析
利用Bioedit 7.1.11软件[17]对获得的基因组片段序列进行整理、拼接以及校正，将经手工矫正后的各基因组片段的序列提交至GenBank，其登录号见表 1。采用MP法与贝叶斯法(BI)分别构建8个单独基因以及6个叶绿体联合基因的系统发育关系。PAUP 4.0b10软件进行MP法分析时，空位(Gap)作为缺失状态，采用启发式搜索(Heuristic search)进行1 000次随机加入，以TBR(tree-bisection-reconnection)进行枝长交换，得到的系统树分枝的可靠性评价使用自展分析(bootstrap, BS)。BI法分析使用MrBayes ver.3.2，利用Modeltest ver.2.2进行模型和参数估计，选择最适碱基替代模型。贝叶斯推断以随机树(random tree)开始分析，2条马尔可夫链(Markov Chain Monte Carlo, MCMC)同时运行，温度参数值设为0.2，共运行1 000 000代，每100代保留1棵树。前25%代作为动态老化样本(burn-in samples)被舍弃，仅对剩余的静态抽样计算一致性树(consensus tree)并计算相关参数，一致性树上的各分支可信度评价用后验概率(posterior probability, PP)表示。
2. 结果与分析
2.1 6个叶绿体基因单独及联合数据分析
在6个叶绿体基因单独分析中，云南金钱槭与其属内姐妹种金钱槭形成单系群的有rbcL(BS=87%，PP=100%)与rpl16(BS=90%，PP=98%)片段；psbA-trnH、trnD-trnT与trnL-trnF片段能把金钱槭属与槭属植物很好分开，但不支持属内两个物种形成单系群；而psbM-trnD片段的结果显示金钱槭属与槭属混在一起，形成并系群。
在6个叶绿体基因联合分析中，覆盖云南金钱槭、金钱槭、槭树属内11个分类组以及外群七叶树属。七叶树属以及部分槭树属内分类组物种的基因片段来源于同一属及同一组的不同物种(表 1)。其中Aesculus flara缺psbM-trnD片段由Aesculus glabra的psbM-trnD片段代替，来象征七叶树属；尖齿枫组中Acer barbinerve缺psbM-trnD与trnD-trnT片段由Acer acuminatum的psbM-trnD与trnD-trnT片段代替，来象征尖齿枫组(Acer sect. Arguta)；坚果枫组中Acer sinopurpurascens缺rbcL片段由Acer sterculiaceum subsp. franchetii的rbcL片段代替，来象征坚果枫组(Acer sect. Lithocarpa)，此外，扁果枫组(Acer sect. Hyptiocarpa)、疏毛枫组(Acer sect. Pubescentia)、桐状枫组(Acer sect. Platanoidea)因序列不全，不列入6个叶绿体基因联合分析中。6个叶绿体基因联合矩阵的序列总长度为4 860 bp，变异位点为685个，其中信息位点为97个。联合数据的最优进化模型为GTR+G。最大简约法分析得出1棵最简约树，树长为913，一致性指数(consistency index)为0.863 1，保持性指数(retension index)为0.362 2。在叶绿体基因联合最简约树中(图 1)，槭属物种(BS = 99%，PP = 100%)被强烈支持形成单系群，而金钱槭属物种虽也被认为是单系群，但其支持度相对较低(BS = 52%，PP = 89%)。此外，金钱槭属与槭属被强烈支持为姐妹属(BS = 100%，PP = 100%)。最大简约法与贝叶斯分析结果显示，两种分析方法所构建的系统树拓扑结构较为一致，在自展支持率上，贝叶斯分析略高于最大简约法。
![]()
图 1 基于叶绿体基因合并序列分析构建的最大简约树菱形◆表示云南金钱槭的系统位置。分支上部数值表示MP分析的抽样自展值(≥50%), 分支下部数值是贝叶斯分析的后验概率值(≥50%)。以下各图同此。Figure 1. Single most parsimonious tree generated from combined chloroplast data◆ represents the phylogenetic status of D. dyeriana. Values at nodes indicate bootstrap support under maximum likelihood and parsimony. Bootstrap values (≥50%) are shown above the branches, and Bayesian posterior values (≥50%) are indicated below the branches. The same as below.2.2 CHS序列分析
CHS序列分析中涵盖了云南金钱槭、金钱槭、槭树属内11个物种(分别代表 11个组)及2个七叶树属物种(表 1)。CHS矩阵的序列总长度为526 bp, 变异位点为78个，其中信息位点为64个。数据的最优进化模型为HKY+G。最大简约法分析得出1棵最简约树，树长为114，一致性指数(CI)为0.736 8，保持性指数(RI)为0.687 5。在CHS序列数据分析中(图 2)，云南金钱槭与其姐妹种金钱槭形成高支持度的单系群(BS = 99%，PP = 99%)，但与6个叶绿体基因分析结果不同，金钱槭属与槭树属间关系较为模糊。
2.3 ITS序列分析
ITS序列分析中涵盖了云南金钱槭、金钱槭、槭树属内14个物种(分别代表 14个组)及1个七叶树属物种(表 1)。ITS矩阵的序列总长度为559 bp，变异位点为178个，其中信息位点为81个。联合数据的最优进化模型为GTR+G。最大简约法分析得出1棵最简约树，树长为415，一致性指数(CI)为0.566 3，保持性指数(RI)为0.174 3。在ITS序列数据分析中(图 3)，金钱槭属物种与槭属物种混在一起，为并系群。其中，云南金钱槭内嵌于槭属内部，与梣叶槭(Acer negundo)聚合在系统树的末端。
3. 讨论与结论
方文培[1]、Xu等[16]根据云南金钱槭冬芽裸露，叶系为羽状复叶(7 ~ 15小叶)以及果实具有圆形翅等特征，与具有类似形态的金钱槭一并归入金钱槭属。而以上的形态差异也正是金钱槭属与槭属主要的形态分类依据。在以往关于槭属的系统学研究中，出现过多次云南金钱槭的系统位置内嵌于槭属内部的现象[2, 4-5]。更有甚者，因金钱槭属物种不能与槭属物种区分，从而把2个属作为一个整体进行讨论[3]。在本研究中，通过叶绿体基因组合、CHS片段和ITS片段的分析显示，的确存在云南金钱槭系统位置不一致的现象。在6个叶绿体基因组合构建的系统发育树中，云南金钱槭与金钱槭形成单系群，且金钱槭属与槭属互为姐妹属(BS = 100%, PP = 100%)(图 1)，此结果支持Renner等[18]基于叶绿体基因构建的系统发育树。值得注意的是，以往基于叶绿体基因的研究中，曾出现过金钱槭属内物种内嵌于槭属内部[2-3]以及金钱槭属与槭属之间形成并系群的状况[5]，但这类状况多出现在基于少量叶绿体基因的研究结果(单个或两个)。而在本研究中增加叶绿体基因数量之后，云南金钱槭所在的金钱槭属与槭属之间被高度支持为姐妹属(图 1)。所以，从基于叶绿体基因遗传分析的角度来讲，支持前人形态分类学[1, 16]的结果。
相较于母系遗传的叶绿体基因，双亲遗传的核基因可能蕴含着更加丰富的系统发育信息[19]。在CHS基因构建的系统发育树中，云南金钱槭所在的金钱槭属为单系群并被强烈支持(BS=99%, PP=99%)，但与槭属间的亲缘关系并不明确(图 2)。而且，核CHS基因对于槭属属内物种间亲缘关系的解析力较差(图 2)。这可能是由于核CHS基因缺乏足够的变异信息位点所导致。而在核ITS基因构建的系统发育树中，云南金钱槭的系统位置位于槭属内部(图 3)，此前也有类似的研究结果[4, 8]。由于核基因相较于叶绿体基因具有较快的进化速率[20]，以及金钱槭属与槭属的分化时间约为62百万年前[18]。所以，在核ITS基因中，云南金钱槭与槭属物种可能因长时间的适应性进化导致序列出现趋同，而单亲遗传的叶绿体基因的变异相对保守，反应出较为原始的系统关系。需要指出的是，核ITS基因存在多拷贝的情况，况且槭树科的物种间杂交较为常见[3]，协同进化的时间可能不足以使多个拷贝变为相同的基因型。所以，克隆挑选出直系同源拷贝进行系统发育重建可能是下一步的研究方向。
综上所述，6个叶绿体基因组合与核CHS基因的序列分析结果，均支持云南金钱槭所在的金钱槭属为单系群，并且6个叶绿体基因组合高度支持金钱槭属与槭属互为姐妹属。核ITS基因虽然支持云南金钱槭归入槭属之内，但多拷贝的问题值得注意。结合形态学证据，尤其是云南金钱槭与槭属植物翅果形状的明显差异[1, 16]，建议云南金钱槭的系统位置继续放在金钱槭属之内。此外，完善槭属物种的采样种类以及增加实验所选择的基因数量，将有助于进一步准确界定云南金钱槭的系统位置。
致谢 感谢台湾师范大学廖培钧副教授、黄秉宏博士和内蒙古大学柴诗瑶老师对论文初稿的修改。
-
![]()
图 1 天然红松居群各表型性状的聚类热图
A. 种长；B. 种宽；C. 百粒质量；D. 球果长；E. 球果宽；F. 球果质量；G. 每层鳞片数量；H. 针叶长；I. 针叶宽；J. 针叶厚；K. 针叶束粗；L. 针叶长/针叶厚；M. 针叶束粗/针叶厚。A, seed length; B, seed diameter; C, hundred seed mass; D, cone length; E, cone diameter; F, cone mass; G, number of scale per layer; H, needle length; I, needle diameter; J, needle thick; K, fascicle width; L, ratio of needle length to needle thick; M, ratio of fascicle width to needle thick.
Figure 1. Clustering heat map of phenotypic traits of natural P. koraiensis populations
表 1 天然红松6个居群的地理位置
Table 1 Geographical location of 6 populations of natural Pinus koraiensis
居群
Population采集地
Collection place经度
Longitude纬度
Latitude海拔
Altitude/m无霜期
Frost-free season/dP1 靖宇县龙湾林场 Longwan Forest Farm, Jingyu County 126°36′09″E 42°05′05″N 836 110 P2 舒兰市平安村 Ping’an Village, Shulan City 127°25′12″E 44°00′36″N 207 125 P3 长白朝鲜自治县十八道沟村
Eighteen Daogou Village, Changbai Korean Autonomous County128°13′05″E 41°44′44″N 837 140 P4 抚松县岗西村 Gangxi Village, Fusong County 127°28′48″E 42°25′42″N 726 110 P5 汪清县丽城村 Licheng Village, Wangqing County 129°45′68″E 43°19′59″N 216 170 P6 白山市江源区振兴村 Zhenxing Village, Jiangyuan District, Baishan City 127°13′12″E 42°00′36″N 624 135 
下载: 导出CSV
表 2 天然红松居群各表型性状的平均值
Table 2 Average values of phenotypic traits in natural P. koraiensis populations
性状 Trait P1 P2 P3 P4 P5 P6 平均 Average 种长
Seed length (SL)/mm15.70 ± 0.50 14.85 ± 0.21 15.52 ± 0.87 15.06 ± 0.22 15.35 ± 0.91 14.83 ± 0.56 15.22 ± 0.69 种宽
Seed diameter (SD)/mm9.76 ± 0.21 10.10 ± 0.13 10.90 ± 0.26 9.82 ± 0.16 10.50 ± 0.51 10.16 ± 0.48 10.21 ± 0.51 百粒质量
Hundred seed mass (HSM)/g61.67 ± 6.48 65.00 ± 5.72 65.33 ± 9.37 59.83 ± 7.01 60.67 ± 8.98 59.67 ± 6.69 62.03 ± 7.74 球果长
Cone length (CL)/cm12.75 ± 0.9 13.26 ± 1.12 12.62 ± 1.10 12.49 ± 0.89 12.73 ± 1.08 13.00 ± 0.63 12.81 ± 0.99 球果宽
Cone diameter (CD)/cm7.29 ± 0.77 7.57 ± 1.07 6.84 ± 0.49 6.36 ± 0.41 6.89 ± 0.57 7.48 ± 0.49 7.07 ± 0.79 球果质量
Cone mass (CM)/g165.87 ± 24.02 194.91 ± 52.10 216.89 ± 44.89 164.66 ± 25.05 229.59 ± 53.89 201.76 ± 31.48 195.61 ± 46.66 每层鳞片数量
Number of scale per layer (NSPL)5.62 ± 0.35 5.34 ± 0.48 5.27 ± 0.39 5.61 ± 0.26 4.76 ± 0.30 5.71 ± 0.19 5.39 ± 0.47 针叶长
Needle length (NL)/mm107.10 ± 11.18 92.46 ± 15.59 90.36 ± 11.15 88.22 ± 14.12 101.30 ± 12.30 105.90 ± 14.33 97.56 ± 15.06 针叶宽
Needle diameter (ND)/mm1.12 ± 0.23 1.25 ± 0.17 1.17 ± 0.13 1.20 ± 0.16 1.31 ± 0.21 1.13 ± 0.13 1.20 ± 0.18 针叶厚
Needle thick (NT)/mm0.81 ± 0.19 0.90 ± 0.13 0.84 ± 0.10 0.85 ± 0.09 0.95 ± 0.16 0.82 ± 0.09 0.86 ± 0.14 针叶束粗
Fascicle width (FW)/mm1.66 ± 0.32 1.82 ± 0.24 1.60 ± 0.19 1.73 ± 0.24 1.88 ± 0.26 1.63 ± 0.18 1.72 ± 0.26 针叶长/针叶厚
NL/NT140.19 ± 43.07 103.65 ± 13.51 108.33 ± 18.77 105.58 ± 21.14 105.91 ± 27.91 131.04 ± 24.60 115.78 ± 29.84 针叶束粗/针叶厚
FW/NT2.10 ± 0.38 2.05 ± 0.21 1.90 ± 0.11 2.05 ± 0.18 1.98 ± 0.20 1.99 ± 0.15 2.01 ± 0.23
下载: 导出CSV
表 3 天然红松居群各表型性状的方差分析
Table 3 Variance analysis of phenotypic traits of natural P. koraiensis populations
性状
Trait均方Mean square (MS) 随机误差
Random
errorF 居群间
Among
populations居群内
Within
population居群间
Among
populations居群内
Within
population种长 SL 116.338 8 11.261 6 0.535 0 10.33** 21.05** 种宽SD 114.939 4 5.231 2 0.554 8 21.97** 9.43** 百粒质量 HSM 494.666 7 150.632 2 4.166 7 3.28** 36.15** 球果长 CL 8.339 1 2.648 0 0.368 3 3.15** 7.19** 球果宽 CD 17.325 7 1.382 4 0.249 8 12.53** 5.53** 球果质量 CM 62 592.360 0 4 516.614 0 663.124 1 13.86** 6.81** 每层鳞片数量 NSPL 11.238 8 0.386 4 0.127 9 29.08** 3.02** 针叶长 NL 6 109.012 0 523.820 2 58.217 2 11.66** 9.00** 针叶宽 ND 0.448 4 0.092 2 0.013 7 4.86** 6.74** 针叶厚 NT 0.255 7 0.051 4 0.007 4 4.97** 6.94** 针叶束粗 FW 0.998 0 0.176 9 0.020 0 5.64** 8.83** 针叶长/针叶厚 NL/NT 21 116.530 0 1 966.087 0 189.983 1 10.74** 10.35** 针叶束粗/针叶厚 FW/NT 0.408 1 0.152 1 0.050 2 2.68 3.03** 注：**表示差异极显著(P < 0.01)。下同。Notes: ** means difference is significant at the 0.01 level. Same as below.
下载: 导出CSV
表 4 天然红松居群各表型性状的方差分量分析
Table 4 Variance component analysis of phenotypic traits of natural P. koraiensis populations
性状
Trait方差分量
Variance component方差分量百分比
Percentage of variance component/%表型分化系数
Differentiation
coefficient of phenotypic
trait (Vst)/%居群间
Among populations居群内
Within population随机误差
Random error居群间
Among populations居群内
Within population种长 SL 3.502 6 1.787 8 0.535 0 60.13 30.69 66.21 种宽 SD 3.656 9 0.779 4 0.554 8 73.27 15.62 82.43 百粒质量 HSM 11.467 8 24.410 9 4.166 7 28.64 60.96 31.96 球果长 CL 0.189 7 0.380 0 0.368 3 20.23 40.51 33.30 球果宽 CD 0.531 4 0.188 8 0.249 8 54.79 19.46 73.79 球果质量 CM 1 935.858 0 642.248 3 663.124 1 59.73 19.82 75.09 每层鳞片数量 NSPL 0.361 7 0.430 1 0.127 9 67.90 8.09 89.36 针叶长 NL 186.173 1 77.600 5 58.217 2 57.82 24.10 70.58 针叶宽 ND 0.011 9 0.013 1 0.013 7 30.71 33.85 47.58 针叶厚 NT 0.006 8 0.007 3 0.007 4 31.61 34.04 48.15 针叶束粗 FW 0.027 4 0.026 2 0.020 0 37.23 35.57 51.14 针叶长/针叶厚 NL/NT 638.348 1 296.017 3 189.983 1 56.77 26.33 68.32 针叶束粗/针叶厚 FW/NT 0.008 5 0.017 0 0.050 2 11.27 22.43 33.44 平均 Average 45.39 28.57 59.33
下载: 导出CSV
表 5 天然红松居群各表型性状的表型变异系数
Table 5 Phenotypic variation coefficients of phenotypic traits in natural P. koraiensis populations
% 性状 Trait P1 P2 P3 P4 P5 P6 平均值 Average 种长 SL 3.18 1.44 1.46 1.46 5.96 3.78 2.88 种宽 SD 2.19 1.29 2.36 1.67 4.82 4.75 2.85 百粒质量 HSM 10.50 8.81 14.34 11.71 14.80 11.21 11.89 球果长 CL 7.31 8.42 8.72 7.15 8.52 4.85 7.50 球果宽 CD 10.58 14.11 7.13 6.49 8.25 6.61 8.86 球果质量 CM 14.48 26.73 20.70 15.21 23.47 15.60 19.37 每层鳞片数量 NSPL 6.29 9.03 7.45 4.60 6.28 3.38 6.17 针叶长 NL 10.44 16.86 12.34 16.01 12.14 13.53 13.55 针叶宽 ND 20.79 13.23 10.96 13.13 15.71 11.20 14.17 针叶厚 NT 23.28 14.10 11.32 10.69 16.69 11.41 14.58 针叶束粗 FW 18.99 13.34 11.98 14.06 13.62 11.21 13.87 针叶长/针叶厚 NL/NT 30.72 13.04 17.33 20.02 26.35 18.78 21.04 针叶束粗/针叶厚 FW/NT 18.23 10.49 5.80 8.90 10.15 7.29 10.14 平均 Average 13.61 11.61 10.15 10.09 12.83 9.51 11.30
下载: 导出CSV
表 6 天然红松居群各表型性状的相关性
Table 6 Correlations of phenotypic traits of natural P. koraiensis populations
性状
Trait种长
SL种宽
SD百粒质量
HSM球果长
CL球果宽
CD球果质量
CM每层鳞片
数量
NSPL针叶长
NL针叶宽
ND针叶厚
NT针叶束粗
FW针叶长/
针叶厚
NL/NT种宽 SD 0.424** 百粒质量 HSM 0.537** 0.404** 球果长 CL −0.011 0.039 0.039 球果宽 CD −0.036 −0.048 0.040 0.541** 球果质量 CM −0.077 0.393** 0.065 0.561** 0.484** 每层鳞片数 NSPL −0.074 −0.337** −0.019 0.313** 0.485** 0.053 针叶长 NL −0.070 −0.164* −0.271** 0.010 0.119 −0.059 −0.011 针叶宽 ND −0.128 −0.012 −0.077 −0.041 −0.089 0.018 −0.193** 0.098 针叶厚 NT −0.138 0.003 −0.070 −0.082 −0.067 0.005 −0.206** 0.128 0.861** 针叶束粗 FW −0.145 −0.110 −0.057 −0.014 −0.020 −0.007 −0.158* 0.128 0.847** 0.809** 针叶长/针叶厚
NL/NT0.096 −0.126 −0.087 0.091 0.172* −0.020 0.176* 0.528** −0.600** −0.684** −0.513** 针叶束粗/针叶厚
FW/NT0.062 −0.177* 0.041 0.075 0.093 −0.043 0.102 −0.008 −0.194** −0.459** 0.095 0.477** 注:*表示相关性达显著水平（P< 0. 05）。下同。
Note: * means correlation is significant at P <0. 05 level. Same as below.
下载: 导出CSV
表 7 天然红松居群各表型性状与环境因素的相关性
Table 7 Correlations between phenotypic traits and environmental factors in natural P. koraiensis populations
环境因素
Environment factor种长
SL种宽
SD百粒质量
HSM球果长
CL球果宽
CD球果质量
CM每层鳞片
数量
NSPL针叶长
NL针叶宽
ND针叶厚
NT针叶束粗
FW针叶长/
针叶厚
NL/NT针叶束粗/
针叶厚
FW/NT经度 Longitude 0.113 0.671 0.006 −0.225 −0.361 0.811 −0.932** −0.154 0.797 0.815* 0.568 −0.610 −0.592 纬度 Latitude −0.422 −0.107 0.208 0.580 0.264 0.195 −0.499 −0.179 0.803 0.807 0.903* −0.545 0.280 海拔 Altitude 0.474 −0.098 −0.076 −0.594 −0.312 −0.480 0.628 0.011 −0.844* −0.882* −0.889* 0.507 −0.041 无霜期
Frost-free season0.086 0.733 −0.040 −0.033 −0.051 0.934** −0.827* 0.176 0.581 0.652 0.373 −0.309 −0.672
下载: 导出CSV
表 8 天然红松居群各表型性状的主成分分析
Table 8 Principal component analysis of phenotypic traits of natural P. koraiensis populations
项目
Item主成分Ⅰ
Component Ⅰ主成分Ⅱ
Component Ⅱ主成分Ⅲ
Component Ⅲ主成分Ⅳ
Component Ⅳ主成分Ⅴ
Component Ⅴ特征值 Eigenvalue 3.45 2.22 2.19 1.31 1.18 贡献率 Contribution rate/% 26.50 17.10 6.88 10.11 9.10 累计贡献率 Cumulative contribution rate/% 26.50 43.60 60.48 70.59 79.69 种长 SL −0.19 0.61 −0.29 0.38 0.21 种宽 SD 0.02 0.78 −0.21 0.20 −0.30 百粒质量 HSM −0.10 0.69 −0.22 0.13 0.41 球果长 CL −0.22 0.30 0.74 −0.03 −0.03 球果宽 CD −0.28 0.18 0.79 0.00 0.05 球果质量 CM −0.10 0.56 0.59 0.05 −0.29 每层鳞片数量 NSPL −0.36 −0.14 0.52 −0.33 0.33 针叶长 NL −0.04 −0.43 0.23 0.70 −0.38 针叶宽 ND 0.89 0.00 0.23 0.20 0.13 针叶厚 NT 0.94 0.03 0.20 0.09 −0.05 针叶束粗 FW 0.81 −0.09 0.28 0.31 0.34 针叶长/针叶厚 NL/NT −0.80 0.30 0.03 0.47 −0.11 针叶束粗/针叶厚 FW/NT −0.41 −0.19 0.06 0.38 0.63
下载: 导出CSV
表 9 天然红松居群的多性状综合评价
Table 9 Comprehensive evaluation of multiple traits of natural P. koraiensis populations
种实性状 Fruit trait 针叶性状 Needle trait 居群 Population Qi 居群 Population Qi P3 1.983 P5 1.986 P5 1.973 P2 1.932 P2 1.928 P1 1.895 P6 1.916 P6 1.893 P1 1.895 P4 1.885 P4 1.872 P3 1.863
下载: 导出CSV
-
[1] 周雪燕, 高海燕, 李召珉, 等. 基于生长与结实评价红松种子园亲本[J]. 植物研究, 2020, 40(3):376−385. doi: 10.7525/j.issn.1673-5102.2020.03.008 Zhou X Y, Gao H Y, Li Z M, et al. Evaluating parents of Pinus koraienis seeds orchard with growth and fruiting[J]. Bulletin of Botanical Research, 2020, 40(3): 376−385. doi: 10.7525/j.issn.1673-5102.2020.03.008
[2] 丁宝永, 张树森, 张世英. 红松人工林季节周期生长规律的研究[J]. 东北林学院学报, 1981(4):19−32. Ding B Y, Zhang S S, Zhang S Y. Research on the growth law of seasonal period of Pinus koraiensis plantations[J]. Journal of North-eastern Forestry Institute, 1981(4): 19−32.
[3] 马建路, 庄丽文, 陈动, 等. 红松的地理分布[J]. 东北林业大学学报, 1992(5):40−48. Ma J L, Zhuang L W, Chen D, et al. Geographic distribution of Pinus koraienis in the world[J]. Journal of Northeast Forestry University, 1992(5): 40−48.
[4] Wang F, Liang D Y, Pei X N, et al. Study on the physiological indices of Pinus sibirica and Pinus koraiensis seedlings under cold stress[J]. Journal of Forestry Research, 2018, 30(4): 1255−1265.
[5] 于大炮, 周旺明, 包也, 等. 天保工程实施以来东北阔叶红松林的可持续经营[J]. 生态学报, 2015, 35(1):10−17. Yu D P, Zhou W M, Bao Y, et al. Forest management of Korean pine and broadleaf mixed forest in Northeast China since the implementation of Natural Forest Protection Project[J]. Acta Ecologica Sinica, 2015, 35(1): 10−17.
[6] 程春香, 毛子军, 靳世波, 等. 小兴安岭北部原始阔叶红松林红松结实气候敏感性及其种子年机制探讨[J]. 植物研究, 2017, 37(1):118−127. Cheng C X, Mao Z J, Jin S B, et al. Sensitivity of fruiting for Pinus koraiensis to climate change and mechanisms of masting in the original broad-leaved Korean pine forest in north Xiaoxing’an Mountain, China[J]. Bulletin of Botanical Research, 2017, 37(1): 118−127.
[7] 梁德洋, 金允哲, 赵光浩, 等. 50个红松无性系生长与木材性状变异研[J]. 北京林业大学学报, 2016, 38(6):51−59. Liang D Y, Jin Y Z, Zhao G H, et al. Variance analyses of growth and wood characteristics of 50 Pinus koraiensis clones[J]. Journal of Beijing Forestry University, 2016, 38(6): 51−59.
[8] 王卫, 史绍林, 刘洋. 人工红松生殖生长研究现状与展望[J]. 防护林科技, 2015(3):96−98. Wang W, Shi S L, Liu Y. Current status and prospects of research on reproductive growth of artificial Pinus koraiensis[J]. Protection Forest Science and Technology, 2015(3): 96−98.
[9] 刘晓婷, 李嘉琪, 李峪曦, 等. 红松半同胞家系变异分析及优良家系选择[J]. 分子植物育种, 2020, 18(13):4473−4482. Liu X T, Li J J, Li Y X, et al. Variation analysis and elite family selection of Pinus koraiensis half-sib families[J]. Molecular Plant Breeding, 2020, 18(13): 4473−4482.
[10] 朱立南. 辽宁地区天然红松种群遗传多样性研究[D]. 大连: 辽宁师范大学, 2011. Zhu L N. The study on genetic diversity of natural Korean pine populations from Liaoning Area[D]. Dalian: Liaoning Normal University, 2011.
[11] 张秦徽, 王洪武, 姜国云, 等. 红松半同胞家系变异分析及选择研究[J]. 植物研究, 2019, 39(4):557−567. doi: 10.7525/j.issn.1673-5102.2019.04.010 Zhang Q H, Wang H W, Jiang G Y, et al. Variation analysis and selection of Pinus koraiensis half-sib families[J]. Bulletin of Botanical Research, 2019, 39(4): 557−567. doi: 10.7525/j.issn.1673-5102.2019.04.010
[12] 贾俊玲. 吉林天然红松林遗传多样性的ISSR分析[D]. 大连: 辽宁师范大学, 2011. Jia J L. Analysis of genetic diversity of Jilin natural Pinus koraiensis population by ISSR marker[D]. Dalian: Liaoning Normal University, 2011.
[13] 王振宇, 陈小强. 红松籽油调节血脂作用研究[J]. 特产研究, 2004(1):1−10. doi: 10.3969/j.issn.1001-4721.2004.01.001 Wang Z Y, Chen X Q. Study on effects of Pinus koraiensis seed oil on blood lipid[J]. Special Wild Economic Animal and Plant Research, 2004(1): 1−10. doi: 10.3969/j.issn.1001-4721.2004.01.001
[14] 姜媛秀. 黑龙江地区天然红松种群遗传多样性的ISSR分析[D]. 大连: 辽宁师范大学, 2010. Jiang Y X. Genetic diversity of Piuns koraiensis populations from Heilongjiang Area revealed by ISSR analysis[D]. Dalian: Liaoning Normal University, 2010.
[15] 张巍, 王清君, 郭兴. 红松不同种源的遗传多样性分析[J]. 森林工程, 2017, 33(2):17−21. doi: 10.3969/j.issn.1006-8023.2017.02.005 Zhang W, Wang Q J, Guo X. Study on the genetic diversity of Pinus koraiensis in different provenances[J]. Forest Engineering, 2017, 33(2): 17−21. doi: 10.3969/j.issn.1006-8023.2017.02.005
[16] 欧汉彪, 林建勇, 李娟, 等. 不同种源楠木叶片表型性状变异分析[J]. 西南农业学报, 2020, 33(3):637−644. Ou H B, Lin J Y, Li J, et al. Analysis of variation of leaf phenotypic traits in different provenances of Phoebe bournei[J]. Southwest China Journal of Agricultural Sciences, 2020, 33(3): 637−644.
[17] Zhao X Y, Hou W, Zheng H Q, et al. Analyses of genotypic variation in white poplar clones at four sites in China[J]. Silvae Genetica, 2013, 62(4−5): 187−195.
[18] Zhao X Y, Li Y, Zheng M, et al. Comparative analysis of growth and photosynthetic characteristics of (Populus simonii × P. nigra) × (P. nigra × P. simonii) hybrid clones of different ploidides[J]. PLoS ONE, 2015, 10(4): e0119259. doi: 10.1371/journal.pone.0119259
[19] 邓丽丽, 蔡年辉, 孙琪, 等. 云南松不同茎干类型群体针叶性状表型多样性比较[J]. 西南林业大学学报, 2016, 36(3):30−37. Deng L L, Cai N H, Sun Q, et al. Compare on the needle phenotypic variations between the different type of trunk populations of Pinus yunnanensis[J]. Journal of Southwest forestry University, 2016, 36(3): 30−37.
[20] Yin S P, Xiao Z H, Zhao G H, et al. Variation analyses of growth and wood properties of Larix olgensis clones in China[J]. Journal of Forestry Research, 2017, 28(4): 687−697. doi: 10.1007/s11676-016-0359-2
[21] 包文泉, 乌云塔娜, 杜红岩, 等. 西藏光核桃表型性状遗传多样性分析[J]. 分子植物育种, 2018, 16(16):5463−5473. Bao W Q, Wuyuntana, Du H Y, et al. Genetic diversity analysis of the phenotypic traits of Amygdalus mira from Tibet Plateau[J]. Molecular Plant Breeding, 2018, 16(16): 5463−5473.
[22] Xu Y L, Woeste K, Cai N H, et al. Variation in needle and cone traits in natural populations of Pinus yunnanensis[J]. Journal of Forestry Research, 2016, 27(1): 41−49. doi: 10.1007/s11676-015-0153-6
[23] 李斌, 顾万春, 卢宝明. 白皮松天然群体种实性状表型多样性研究[J]. 生物多样性, 2002, 10(2):181−188. doi: 10.3321/j.issn:1005-0094.2002.02.008 Li B, Gu W C, Lu B M. A study on phenoytpic diversity of seeds and cones characteristies in Pinus bungeana[J]. Biodiversity Science, 2002, 10(2): 181−188. doi: 10.3321/j.issn:1005-0094.2002.02.008
[24] 张含国, 高士新, 张敏莉, 等. 长白落叶松天然群体遗传结构的研究[J]. 东北林业大学学报, 1995(6):21−31. Zhang H G, Gao S X, Zhang M L, et al. Study on natural population genetic structure of Larix olgensis Henry[J]. Journal of Northeast Forestry University, 1995(6): 21−31.
[25] 辜云杰, 罗建勋, 吴远伟, 等. 川西云杉天然种群表型多样性[J]. 植物生态学报, 2009, 33(2):291−301. Gu Y J, Luo J X, Wu Y W, et al. Phenotypic diversity in natural populations of Picea balfouriana in Sichuan, China[J]. Chinese Journal of Plant Ecology, 2009, 33(2): 291−301.
[26] 葛颂, 王明庥, 陈岳武. 用同工酶研究马尾松群体的遗传结构[J]. 林业科学, 1988(4):399−409. Ge S, Wang M X, Chen Y W. An analysis of population genetic structure of Masson pine by sozyme technique[J]. Scientia Silvae Sinicae, 1988(4): 399−409.
[27] Metougui M L, Mokhtari M, Maughan P J, et al. Morphological variability, heritability and correlation studies within an argan tree population (Argania spinosa (L.) Skeels) preserved in situ[J]. International Journal of Agriculture and Forestry, 2017, 7(2): 42−51.
[28] 张翠琴, 姬志峰, 林丽丽, 等. 五角枫种群表型多样性[J]. 生态学报, 2015, 35(16):5343−5352. Zhang C Q, Ji Z F, Lin L L, et al. Phenotypic diversity of Acer mono Maxim. population[J]. Acta Ecologica Sinica, 2015, 35(16): 5343−5352.
[29] 邱黛玉, 蔺海明, 陈垣, 等. 经纬度和海拔对当归成药期植株长势和早期抽薹的影响[J]. 草地学报, 2010, 18(6):838−843. Qiu D Y, Lin H M, Chen Y, et al. Effects of latitude, longitude and altitude on angelica growth and early bolting in medicine formation period[J]. Acta Agrestia Sinica, 2010, 18(6): 838−843.
[30] 刘永良. 油松地理种群针叶性状指标的变异与趋势[D]. 北京: 北京林业大学, 2011. Liu Y L. The variation and trend of needle trait indexes of Pinus tabulaeformis Carr. geographic populations[D]. Beijing: Beijing Forestry University, 2011.
[31] 林丽丽. 茶条槭种群表型差异研究及潜在分布区预测[D]. 太谷: 山西师范大学, 2015. Lin L L. Phenotypic variations and potential distribution area prediction of Acer ginnala populations[D]. Taigu: Shanxi Normal University, 2015.
[32] 侯维海, 王建林, 旦巴, 等. 不同生态因子条件下西藏青稞种子表型性状的相关分析[J]. 核农学报, 2017, 31(10):2063−2071. doi: 10.11869/j.issn.100-8551.2017.10.2063 Hou W H, Wang J L, Danba, et al. Phenotypic correlation analysis of Hullessbarley kernel traits form Tibet Plateau region under the condition of different ecological environment[J]. Journal of Nuclear Agricultural Sciences, 2017, 31(10): 2063−2071. doi: 10.11869/j.issn.100-8551.2017.10.2063
[33] 童跃伟, 唐杨, 陈红, 等. 红松种子园种群表型多样性研究[J]. 生态学报, 2019, 39(17):6341−6348. Tong Y W, Tang Y, Chen H, et al. Phenotypic diversity of Pinus koraiensis populations in a seed orchard[J]. Acta Ecologica Sinica, 2019, 39(17): 6341−6348.
[34] Feng F J, Han S J, Wang H M. Genetic diversity and genetic differentiation of natural Pinus koraiensis population[J]. Journal of Forestry Research, 2006, 17(1): 21−24. doi: 10.1007/s11676-006-0005-5
[35] 冯富娟, 隋心, 张冬东. 不同种源红松遗传多样性的研究[J]. 林业科技, 2008, 33(1):1−4. Feng F J, Sui X, Zhang D D. Studies on the genetic diversity of Pinus koraiensis in different provenance[J]. Forestry Science & Technology, 2008, 33(1): 1−4.
-
期刊类型引用(0)
其他类型引用(4)
计量
- 文章访问数: 1684
- HTML全文浏览量: 483
- PDF下载量: 58
- 被引次数: 4
