榧树天然群体中种子表型特征与化学成分的变异分析

董雷鸣，曾燕如，邬玉芬，黄银芝，武栋，戴文圣

（1.浙江农林大学亚热带森林培育国家重点实验室培育基地，浙江临安 311300；2.浙江省宁海林特技术推广总站，浙江宁海 315600；3.浙江省黄岩区头陀镇人民政府农业办公室，浙江台州 318026）

榧树天然群体中种子表型特征与化学成分的变异分析

董雷鸣1，曾燕如1，邬玉芬2，黄银芝3，武栋1，戴文圣1

（1.浙江农林大学亚热带森林培育国家重点实验室培育基地，浙江临安 311300；2.浙江省宁海林特技术推广总站，浙江宁海 315600；3.浙江省黄岩区头陀镇人民政府农业办公室，浙江台州 318026）

对安徽省黄山市天然榧树Torreya grandis群体内种子表型性状和化学成分的变异进行分析研究，为优良榧树种质资源收集、品种选育及榧树遗传研究奠定基础。结果表明：①种蒲质量、种蒲横径、种蒲纵径、种核横径、种核纵径、种核质量、种形指数、核形指数和出核率的变幅分别为5.24～12.23 g，1.963～3.018 cm，2.255～3.685 cm，1.438～2.270 cm，2.009～3.208 cm，2.30～5.19 g，0.695～0.902，0.555～0.823，34.80%～66.49%；②营养成分粗脂肪、淀粉和蛋白质的质量分数变幅分别为232.96～425.54 g·kg-1，59.95～114.40 g·kg-1和153.22～231.85 g·kg-1；矿质元素磷、钾、钙、镁、铜、锌、锰和镉的变幅分别为0.81～2.29 g·kg-1，23.90～53.72 g·kg-1，80.18～224.38mg·kg-1，765.84～798.57mg·kg-1，8.10～31.60 mg·kg-1，40.43～128.00 mg·kg-1，1.45～23.32 mg·kg-1，0.01～0.11 mg·kg-1；多数单株样品砷和铅呈现未检出状态；③种蒲质量、种蒲横径、种蒲纵径、种核横径、种核纵径及种核质量两两之间极显著正相关（P＜0.01）；种蒲和种核的形状与质量的关系无规律；种蒲质量与种形指数、种核质量与核形指数之间无明显的相关关系，而种核质量与种形指数之间存在显著的负相关关系（P＜0.05）；出核率与种蒲横径、种蒲纵径、种核横径、种核纵径均存在极显著的负相关（P＜0.01），与种核指数极显著正相关（P＜0.01），与种形指数无显著的相关关系；3种营养成分质量分数之间不存在任何的相关关系，除粗脂肪与种核质量间存在显著的负相关关系外（P＜0.05），营养成分与种子表型性状及矿质元素间不存在相关关系。钾与镁、铜与磷、铜与锰的质量分数间均存在显著的正相关（P＜0.01）；锰的质量分数与种形指数和核形指数间呈极显著正相关（P＜0.01）。表3参15

经济林学；榧树；种子表型；化学成分；变异分析；天然群体

榧树Torreya grandis隶属于红豆杉科Taxaceae榧属Torreya，是中国特有的珍稀树种，自然分布于中国江苏南部、浙江、福建北部、江西北部、安徽南部，南至湖南西南部及贵州松桃等中亚热带和北亚热带地区[1-2]。但是榧树资源分布不均，浙江省最多，安徽省黄山市也有较多分布[3]。榧树种子是经济价值很高的坚果，可以榨油、食用或入药；种仁含有丰富的油脂、蛋白质、矿物元素及特殊的维生素[4]，其脂肪酸的主要组成以亚油酸和油酸为主，其中不饱和脂肪酸占到脂肪酸总量的79%以上[5]。榧树多雌雄异株，异花授粉，不同单株之间基因交流的几率很大，加上分布范围广，分布区地理环境差异大，使得天然群体内有较高的遗传多样性，种内性状变异十分复杂[6-8]。沈登锋等[9]的研究发现，天然林里的榧树种子不仅大小差异大，且存在性状特别或良好的种质，不乏一些综合性状优良、品质达到或超过香榧的单株种质[3，9]，这些种质是优良品种选育的基础材料。迄今为止，榧树种下仅香榧T.grandis‘Merrilli，1个栽培类型，还有很大的育种潜力。因此，榧树资源天然的遗传变异是开展榧树育种工作的基础，但尚无对天然群体内种子理化性质变异开展系统研究的报道。本研究以安徽省黄山市天然榧树群体为研究对象，评估天然群体中雌株种子表型及营养成分与部分矿质元素含量的变异，结果可为优良单株的选育、榧树遗传变异研究提供参考。

1 材料与方法

1.1 材料

1.1.1 榧树天然林概况榧树天然林位于安徽省黄山市徽州区呈坎镇小容村（29°57＇N，118°14＇E）。该地属亚热带季风性湿润气候，年最高气温为40.0℃，年最低气温为-10.0℃，年平均气温16.0℃，年平均降水量1 450.0 mm。夏季气温稍低，而冬季气温稍高，日照稍短。全年无霜期240 d。该地分布有大量百年以上的野生榧树，适合开展榧树资源评价和优株选育工作。

1.1.2 样品采集2010年10月榧树种子成熟时节随机选择50株彼此至少相距50 m生长健壮、正常结实的雌株，分别编号、定位（GPS，eTrex Vista，Garmin），并在树干喷涂油漆。采集种子不少于100粒·株-1，用于种实性状分析。由于榧树的经济价值及经济产量主要体现在种实品质上，因此，在样品采集过程中注意与当地农户进行交流，收集农户认为比较好的单株种质。

1.2 方法

1.2.1 种子表型性状的测定从每棵树采集的种子中随机取30粒，用电子数显卡尺（中国上海量具刃具厂）及电子天平（PB1502-L，瑞士）测定带假种皮种蒲的纵径、横径、质量及去除假种皮后种核的纵径、横径、质量，并计算种形指数（种蒲横径/种蒲纵径）、核形指数（种核横径/种核纵径）及出核率[（种核质量/种蒲质量）×100%]。

1.2.2 种子营养成分与矿质元素质量分数的测定表型性状测定后，将种子洗净晾干装入封口袋。取20粒·样-1，去种壳后放入烘箱105℃干燥2 h，然后用粉碎机（DFY-400，温岭市林大机械有限公司）粉碎，并混合均匀，按编号装入封口袋待用。粗脂肪质量分数按照GB/T 5512-2008索氏抽提法进行测定；采用蒽酮比色法和紫外分光光度计（UV-2550，日本）测定种子的淀粉质量分数；利用定氮仪（FOSSKjeltec 2300，瑞典）及凯氏定氮法[10]测定种子的蛋白质质量分数。采用火焰原子吸收法在原子吸收光谱仪（ICE3000 SERIES，美国）上测定钾（K），钙（Ca），镁（Mg），铜（Cu），锌（Zn），锰（Mn）6种矿质元素的质量分数；采用等离子体质谱法在等离子体质谱仪（XSERIES2 ICP-MS，美国）上测定砷（As），镉（Cd），铅（Pb）3种矿质元素；磷（P）元素在紫外分光光度计（UV-2550，日本）上采用钼锑抗比色法单独测定。所有处理重复测定3次。

1.2.3 数据处理种子各化学成分实验测得的数据计算参照武栋[11]的描述。用统计软件SPSS 13.0对种蒲和种核表型性状，以及种子各化学成分数据分析结果进行统计性描述，并两两进行相关分析。

2 结果与分析

2.1 榧树种子表型性状和化学成分的变异

2.1.1 种蒲和种核大小的变异榧树种蒲样品中单粒质量最大值为12.23 g，最小值为5.24 g，平均质量8.96 g。种蒲质量的变异系数最高，为20.1%（表1），说明天然群体中种蒲大小存在一定变异。种核质量的变幅为2.30～5.19 g，变异系数为16.8%（表1）。种蒲横径最大值为3.018 cm，最小1.963 cm，均值为2.368 cm；种蒲纵径的变幅为2.255～3.685 cm，平均2.931 cm。种核横径变幅为1.438～2.270 cm，平均1.762 cm；种核纵径变幅为2.009～3.208 cm，平均2.507 cm。以上4个指标的变异系数相近，分别为10.3%，10.5%，10.2%和10.5%（表1）。

表1 榧树种子表型性状的描述性统计值Table 1 Phenotypic values of seed traits in Torreya grandis

2.1.2 种蒲和种核形状的变异种形指数和核型指数数值范围为0～1，能够反映榧树种蒲和种核的形状变化。值越接近于0，表明种蒲或种核越细长，反之则越接近于球形。种蒲样品多呈椭球形或球形，去除假种皮的种核则一端较为圆润，另一端渐尖，测量数据也证实了这些现象：种形指数变幅范围为0.695～0.902，平均0.811；核形指数最大值为0.823，最小为0.555，平均0.707；变异系数则分别为7.6%和15.3%（表1）。

2.1.3 出核率的变异榧树种子一般食用（如果品质好的话）或用作培育嫁接香榧的砧木，因此出核率是衡量榧树产量的一个重要指标，也是重要的经济性状之一。样品中最高出核率为66.49%，最低34.80%，平均为45.15%，变异系数为15.3%（表1）。

2.1.4 营养成分的变异从表2中可知，本研究中样品粗脂肪质量分数变幅为232.96～425.54 g·kg-1，平均值为337.25 g·kg-1，远低于香榧的平均值578.20 g·kg-1；淀粉质量分数变幅为59.95～114.40 g·kg-1，平均值为85.60 g·kg-1，高于香榧的53.20 g·kg-1；蛋白质质量分数平均174.81 g·kg-1，高于香榧的140.2 g·kg-1，变幅为153.22～231.85 g·kg-1（表2）。蛋白质变异最小，变异系数为8.3%，粗脂肪和淀粉质量分数的变异系数分别为13.2%和16.5%。

表2 榧树种子化学成分的描述性统计值Table 2 Statistically descriptive valuesof seed chemical compositions in Torreya grandis

2.1.5 矿质元素的变异部分矿质元素（如锌）是对人身体有益的元素，而过量有害元素（如砷、镉、铅）的摄入则是有害的。本研究中，磷质量分数的变幅为0.81～2.29 g·kg-1，平均为1.37 g·kg-1；钾元素的变幅为23.90～53.72 g·kg-1，平均为40.98 g·kg-1；钙的变幅为80.18～224.38 mg·kg-1，平均为139.76 mg· kg-1；镁元素的平均质量分数为779.57 mg·kg-1，变幅为765.84～798.57mg·kg-1；铜的变幅为8.10～31.60 mg·kg-1，平均为20.89 mg·kg-1；锌的变幅为40.43～128.00 mg·kg-1，平均为69.50 mg·kg-1；锰的变幅为1.45～23.32mg·kg-1，平均为9.69mg·kg-1。50份样品中分别有21份和49份砷和铅的质量分数低于检出限，最大值分别为0.11mg·kg-1，0.05mg·kg-1；镉的变幅为0.01～0.11 mg·kg-1，平均为0.04 mg·kg-1。各元素变异程度不一，镁的变异系数仅为0.95%；锰的变异系数最大，为50.05%；其他几种的变异系数为13.54%～25.36%。

2.2 种子表型性状与化学成分的相关性

相关分析发现，种蒲质量、种蒲横径、种蒲纵径、种核横径、种核纵径及种核质量两两之间极显著正相关（P＜0.01），Pearson相关系数均＞0.5；由于种形指数、种核指数和出核率由前6个原初指标计算得来，分别与相应的分子呈极显著的正相关，与分母则呈极显著的负相关关系（P＜0.01）；种蒲和种核的形状与质量的关系无规律；种蒲质量与种形指数、种核质量与种核指数之间无明显的相关关系，而种核质量与种形指数之间存在显著的负相关关系（P＜0.05）；出核率与种蒲横径、种蒲纵径、种核横径、种核纵径均存在极显著的负相关，与核形指数极显著正相关（P＜0.01），与种形指数无显著的相关关系（表3）。

3种营养成分之间不存在任何的相关关系。除含油率与种核质量间存在显著的负相关关系外，营养成分与种子表型性状及矿质元素间不存在相关关系。钾与镁、铜与磷、铜与锰间均存在显著的正相关（P＜0.05）；锰元素与种形指数和核形指数极显著正相关（P＜0.01）（表3）。

3 讨论

3.1 榧树种子表型性状的变异

生物的变异和多样性是开展育种工作的前提。榧树的分布范围广，分布区内有极为丰富的变异。种子是榧树主要的经济产出，榧树分布区当地百姓根据种子的外形将榧树分为不同的变异类型，如花生榧、米榧、芝麻榧、象牙榧、圆榧等[9]。本研究较为系统地研究了安徽省黄山市1个天然榧树群体内种子表形性状和化学成分，结果表明：种子表型性状存在一定变异，除了种蒲和种核质量的变异较大（变异系数分别为20.1%和16.8%），其余指标如横纵经、种形指数和核形指数都较小，这些性状的变异均小于黎章矩等[12]在多个天然群体中采集的不同个体种子样品表型性状的变异，他们统计的榧树种蒲质量最大为18.49 g，最小为4.63 g，变异系数为35.15%，这可能与群体内不同个体间存在较大的遗传差异有关。

3.2 榧树种子化学成分的变异

营养成分粗脂肪、淀粉和蛋白质质量分数的变异以淀粉最高，为16.5%；粗脂肪次之，为13.2%；蛋白质最小，仅8.3%，与前人[12]的研究结果相一致。香榧的风味及口感要优于普通的榧籽，主要依赖于其较高的含油率和较低的淀粉质量分数，所以实生榧树种子品质的好坏也取决于脂肪、淀粉和蛋白质的高低[12]。从含油率的测定结果看，该区域内的榧树普遍较低，而淀粉和蛋白质却高于香榧，说明品质近于或高于香榧的个体极少。目前，榧树产地一些居民食用一些风味较好的榧籽，一般的榧籽则出售用于砧木的繁育，榧籽的开发利用相对比较落后，建议开展榧树多目标选育，充分开发利用榧籽中的天然成分。

矿质元素的变异系数差别较大，镁元素的仅为0.95%，而锰元素却高达50.05%，可能与小地形的土壤成分有关。砷、镉、铅3种对人体有害元素的含量绝大部分均远小于香榧[13]及以基于植物源的进口保健品如人参、绿藻、银杏、茶类、大豆磷脂及其他中草药类[14]，多数单株呈现未检出状态，说明该区榧树的生长环境优异，较少受到污染，而香榧林人工集约经营的程度高，人工施肥及使用化学药品可能是导致香榧中这些元素含量高的因素。

3.3 种子表型性状与化学成分的相关性

相关分析的结果表明：多个种子表型性状两两之间极显著正相关。种子的形状与质量的关系无规律，说明种子质量不随种子长短胖瘦的变化而变化。出核率是描述产量的一个重要指标，它与种蒲的纵横径均存在极显著的负相关（P＜0.01），与种核指数极显著正相关（P＜0.01），与种形指数无显著的相关关系，说明种蒲越大，出核率越低，而去除假种皮的种核越接近球形的种子出核率越高，结合种核的质量，可以直接预测每株树大致的种核产量。已有研究表明，种核形状细长的往往品质较好，即可食性较强[3]，但本研究结果表明，除含油率与种核质量间存在显著的负相关关系外，其他营养成分指标与种子表型形状间无相关关系，不能证明种核形状与品质之间存在关联。含油率与种核质量负相关表明种核单位质量越大，其含油率越低，相应的品质也就较差，这与程晓建等[3]的报道一致，该结果可以用于筛选优良种质。此外我们的研究发现，种子越大，所形成的苗木生长越快[13]，那么可以选择相应的母本通过嫁接方式营建种子园，生产的种子可以为香榧提供优良砧本。部分矿质元素如钾与镁、铜与磷、铜与锰的质量分数间均存在显著的正相关（P＜0.05），表明这些元素在榧树种子代谢过程中存在一定的依存关系。锰元素与种形指数和核形指数极显著正相关（P＜0.01），说明外形细长的种子中锰元素质量分数显著高于偏圆的种子，那么在育种中可以选择合适的亲本来提高锰元素的质量分数。

总体上看，该群体中种子的表型性状和营养成分的变异程度不是特别大。曾用SRAP标记对该群体雌性榧树进行遗传多样性分析，结果表明：Nei，s基因多样性指数H和Shannon信息指数I分别为0.166 6和0.262 3[15]，也表明变异程度不是特别高。母本对其半同胞种子的影响如何有待进一步研究。另外，在香榧产业大发展之前，许多地方的农民不了解榧树雌雄异株的特性，砍伐了许多不结果的雄株，以致于目前许多地方雄株花粉不足，需从异地采集花粉进行人工授粉。雄株的减少对群体内变异的变化也起着一定的作用。此外，如果仅仅以含油率作为衡量优株的唯一指标，那么满足这个要求的个体数很少。然而榧树的育种工作需要从多个方面展开，如产量、抗性，尤其是用于香榧嫁接的优良砧木的选育，包括优良雄株的选育，以充分利用现有的且越来越少的天然种质资源。

[1]郑万钧，傅立国.中国植物志：第7卷[M].北京：科学出版社，1978：459.

[2]程晓建，黎章矩，喻卫武，等.榧树的资源分布与生态习性[J].浙江林学院学报，2007，24（4）：383-388.

CHENG Xiaojian，LIZhangju，YUWeiwu，et al.Distribution and ecological characteristics ofTorreya grandisin China [J].JZhejiang For Coll，2007，24（4）：383-388.

[3]程晓建，黎章矩，戴文圣，等.榧树种质资源调查与评价[J].果树学报，2009，26（5）：654-658.

CHENG Xiaojian，LIZhangju，DAIWensheng，et al.Investigation and evaluation ofTorreya grandisgermplasm resource[J].JFruit Sci，2009，26（5）：654-658.

[4]黎章矩，骆成方，程晓建，等.香榧种子成分分析及营养评价[J].浙江林学院学报，2005，26（5）：540-544.

LIZhangju，LUO Chengfang，CHENG Xiaojian，et al.Component analysis and nutrition evaluation of seeds ofTorreya grandis‘Merrillii，[J].JZhejiang For Coll，2005，26（5）：540-544.

[5]陈振德，郑汉臣，傅秋华，等.国产榧属植物种子油含量及其脂肪酸测定[J].中国中药杂志，1998，23（8）：456-457.

CHEN Zhende，ZHENG Hanchen，FU Qiuhua，et al.Determination of oil contents and fatty acids in seeds ofTorreyaArn.in China[J].China JChin Materia Med，1998，23（8）：456-457.

[6]闵会.浙江省香榧天然群体的遗传多样性研究[D].临安：浙江农林大学，2009.

MIN Hui.Genetic Diversity and Germplasm Resources Conservation of Different Torreya grandisFort ex Lindl.[D]. Lin，an：Zhejiang A&F University，2009.

[7]李扬，姚小华，王开良，等.野生香榧种实性状变异研究[J].浙江林业科技，2009，29（3）：35-38.

LIYang，YAO Xiaohua，WANG Kailiang，et al.Variation and correlativity analysis of fruit characters of wild Torreya grandis‘Merrillii，[J].JZhejiang For Sci Technol，2009，29（3）：35-38.

[8]董雷鸣，沈登锋，喻卫武，等.榧树雄株若干性状变异初探[J].浙江农林大学学报，2012，29（5）：715-721.

DONG Leiming，SHEN Dengfeng，YUWeiwu，et al.Trait variations ofmaleTorreya grandistrees[J].JZhejiang A& FUniv，2012，29（5）：715-721.

[9]沈登锋，曾燕如，喻卫武，等.榧树种质资源的收集与种子理化性质的初步分析[J].浙江农林大学学报，2011，28（5）：747-752.

SHEN Dengfeng，ZENG Yanru，YU Weiwu，et al.Collection ofTorreya grandisgermplasm and analysis of seeds，physico-chemical characteristics[J].JZhejiang A&F Univ，2011，28（5）：747-752.

[10]黄晓荣，曹承富，杜世州，等.全自动定氮仪测定小麦籽粒蛋白质[J].安徽农业科学，2009，37（19）：8823-8824.

HUANG Xiaorong，CAO Chengfu，DU Shizhou，et al.Determination on protein in wheat grains by fully automatic azotometer[J].JAnhui Agri Sci，2009，37（19）：8823-8824.

[11]武栋.安徽徽州榧树种群内种子理化性质的分析[D].临安：浙江农林大学，2012.

WU Dong.Analysis of Physico-chemical Properties of the Seeds From a Natural Torreya population in Huizhou，Anhui[D].Lin，an：Zhejiang A&F University，2012.

[12]黎章矩，戴文圣.中国香榧[M].北京：科学出版社，2007.

[13]戴文圣，黎章矩，曹福亮，等.香榧林地土壤及其种子矿物元素含量分析[J].浙江林业科技，2006，26（2）：14-18.

DAIWensheng，LIZhangju，CAO Fuliang，et al.Mineral elements inTorreya grandis‘Merrillii，seeds and their forestland soils[J].JZhejiang For Sci Technol，2006，26（2）：14-18.

[14]王永芳，赵馨，韩宏伟.进口保健食品中铅含量监测与评价[J].中国食品卫生杂志，1999，11（5）：25-26.

WANG Yongfang，ZHAO Xin，HAN Hongwei.Monitoring and evaluation of lead in impported helth food[J].Chin JFood Hygi，1999，11（5）：25-26.

[15]董雷鸣.榧树群体遗传结构及与苗期生长性状关联的SRAP标记分析[D].临安：浙江农林大学，2013.

DONG Leiming.Population Genetic Structure and Association Analysis of Juvenile Growth Traits Using SRAPMarkers in Torreya grandis[D].Lin，an：Zhejiang A&F University，2013.

Variations in phenotypic traits and chemical compositions of seeds from a natural population in Torreya grandis

DONG Leiming1，ZENG Yanru1，WU Yufen2，HUANG Yinzhi3，WU Dong1，DAIWensheng1
（1.The Nurturing Station for the State Key Laboratory of Subtropical Silviculture，Zhejiang A&F University，Lin，an 311300，Zhejiang，China；2.Ninghai Forest Technological Extension Station，Ninghai 315600，Zhejiang，China；3. Agricultural Office，People，s Government of Toutuo Town，Huangyan District，Taizhou 318026，Zhejiang，China）

Variations in phenotypic traits and chemical compositions of Torreya grandis seeds sampled from a natural population located in Xiaorong Village of Chengkan Town，Huangshan City，Anhui Province were analyzed with a view to laying a foundation for elite germplasm collection，selection and breeding of new cultivars，and genetic studies.Results showed as follows：（1）For seedswith and withoutaril，their weight，transverse diameter，vertical diameter shape index and the weight percentage of seedswith aril to those without aril ranged from 5.24-12.23 g，2.30-5.19 g，1.963-3.018 cm，1.438-2.270 cm，2.255-3.685 cm，2.009-3.208 cm，0.695-0.902，0.555-0.823 and 34.80%-66.49%，respectively.（2）The fat，starch and protein contents changed from232.96 to 425.54 g·kg-1，59.95-114.40 g·kg-1，and 153.22-231.85 g·kg-1，respectively.The content of P，K，Ca，Mg，Cu，Zn，Mn and Cd varied from 0.81 to 2.29 g·kg-1，23.90-53.72 g·kg-1，80.18-224.38 mg·kg-1，765.84-798.57 mg·kg-1，8.10-31.60mg·kg-1，40.43-128.00 mg·kg-1，1.45-23.32 mg·kg-1，and 0.01-0.11mg· kg-1，respectively.As and Pb were not detectable in most of samples.And（3）The weight，transverse diameter，and vertical diameter of seedswith and without aril had amost significantly positive linear correlation between each other（P＜0.01），but the relationship between the seed shape and weight of both seedswith or without arilwas irregular.There was no significant correlation between weight and shape index of seedswith and without aril.Butweight of seeds without aril was negatively correlated with shape index of seeds with aril at a significant level（P＜0.05）.The weight percentage of seedswith aril to those withoutarilwasmost significantly in negative correlation with transverse and vertical diameter of seedswith and without aril and positively correlated with shape index of seeds without aril at a most significant level，but not correlated with shape index of seeds with aril.There was no significant correlation between 3 nutritive ingredients，which had no correlation with both phenotypic traits of seeds and chemical composition，either，except the significantly negative correlation between fat percentage and weight of seeds without aril.There was a significant positive correlation between K and Mg，Cu and P as well as Cu and Mn contents.The Mn content was positively correlated with shape index of seedswith and withoutaril atamost significant level.[Ch，3 tab.15 ref.]

cash forestry；Torreya grandis；seed phenotype；chemical composition；variation analysis；natural population

S722.1

A

2095-0756（2014）02-0224-07

2013-05-14；

2013-05-27

浙江省科技计划项目（2008C22003，2012C12904-12）；浙江农林大学研究生科研创新基金资助项目（31220113240135）；浙江省宁波市科技计划项目（2011Y16）

董雷鸣，从事经济林培育与利用研究。E-mail：dongleiming2008@126.com。通信作者：曾燕如，教授，博士，从事经济林培育与利用研究。E-mail：yrzeng@zafu.edu.com