基于毛细管电泳的日本落叶松MSAP反应体系的建立

刁 姝，孙 超，孙晓梅，张守攻

(中国林业科学研究院林业研究所，国家林业局林木培育重点实验室，北京 100091)

基于毛细管电泳的日本落叶松MSAP反应体系的建立

刁 姝，孙 超，孙晓梅，张守攻*

(中国林业科学研究院林业研究所，国家林业局林木培育重点实验室，北京 100091)

[目的]建立基于毛细管电泳的日本落叶松最佳MSAP反应体系，获得条带清晰，信号强度均一的DNA甲基化图谱，为今后开展日本落叶松表观遗传变异研究奠定基础。[方法]以日本落叶松未成熟木质部为试材，设置不同梯度的DNA模板用量以及选择扩增体系中引物浓度、Mg2+浓度、dNTP浓度和exTaq DNA聚合酶浓度，确定各组分不同用量对毛细管电泳效果的影响，从而确立最适的MSAP反应体系。[结果]结果表明模板DNA以及exTaq DNA聚合酶浓度对毛细管电泳结果影响不明显，而高的引物浓度会引起条带的减少。Mg2+及dNTP浓度是影响毛细管电泳结果最为关键的因素，浓度过低或过高会导致毛细管电泳图谱条带明显减少，杂峰显著增多。[结论] 最终确立的反应体系为：DNA 模板100 ng，引物0.75 pmol·μL-1，Mg2+2 mmol·L-1，dNTP 0.375 mmol·L-1和exTaq DNA聚合酶 0.05 U·μL-1。

DNA甲基化；甲基化敏感扩增多态性(MSAP)；毛细管电泳；体系建立；日本落叶松

DNA甲基化是指S-腺苷甲硫氨酸在DNA甲基化转移酶的作用下，将甲基基团共价结合到胞嘧啶5’碳位上的过程[1]。 DNA甲基化是目前研究最为深入的表观遗传修饰。DNA甲基化既发生于对称的CG，CHG序列上，也发生于非对称的CHH(H=A，C或T)序列上[2]。DNA甲基化参与调控植物体多种生长发育过程及表型可塑性，如拟南芥(Arabidopsisthaliana)花发育过程[3]、水稻(Oryzasativa)种子发育过程[4]及红树植物拉关木(Lagunculariaracemosa)对不同环境的适应性[5]等。

目前应用最为广泛的全基因组水平DNA甲基化检测方法主要分为两类，一类依托于高通量测序技术，包括全基因组亚硫酸盐测序法(Whole genome bisulfite sequencing，WGBs)、简化代表性亚硫酸盐测序法(Reduced reperesentation bisulfite sequencing，RRBS)和DNA甲基化免疫共沉淀测序法(Methylated DNA immunoprecipitation sequencing，MeDIP-seq)等[6]，可在全基因组水平获得DNA甲基化信息，但价格昂贵且需要物种的全基因组序列；另一类是基于限制性内切酶的甲基化敏感扩增多样性(methylation-sensitive amplification polymorphism，MSAP)检测技术。MSAP技术是在AFLP技术基础上发展起来的，将AFLP中的内切酶MseI替换成同裂酶HpaII及MspI[7]，通过比较两种同裂酶分别与EcoR I组合获得电泳条带的差异确定CCGG位点甲基化情况。该技术的优点是操作简单、成本低，无需已知参考基因组序列，缺点是只能获得全基因组部分DNA甲基化信息。

日本落叶松(Larixkaempferi)是我国重要的经济和生态树种，具有生长速度快、抗病性强和适应性广的特点[8]，开展日本落叶松DNA甲基化研究有助于在表观遗传水平解析其表型变异。由于参考基因组序列未知，李爱应用基于传统的聚丙烯酰胺凝胶电泳的MSAP技术在日本落叶松中检测了1447个5’-CCGG 位点，并发现日本落叶松杂种优势可能与DNA甲基化模式的大量调整密切相关[9]，但聚丙烯酰胺凝胶电泳耗时长，扩增条带难以确认，灵敏性差，不利于大样本的表观遗传变异研究。基于毛细管电泳(capillary electrophoresis，CE)的MSAP技术可以克服上述缺陷，但对反应体系要求比较高。本研究拟通过分析DNA模板用量及选择性扩增体系中主要成分浓度对毛细管电泳结果的影响，确立最佳反应体系，获得条带清晰、信号强度均一的DNA甲基化图谱，为今后开展日本落叶松表观遗传变异研究奠定基础。

1 材料与方法

1.1 DNA 提取、酶切、连接以及预扩增

试验材料为16年生日本落叶松1.3米胸径处未成熟木质部(图1)，采集于湖北省建始县长岭岗林场日本落叶松原生种源试验林，该试验林包括草津、浅间、富士、伊那、日光以及松本六个种源。随机选取草津种源中的一株作为本研究的试验材料。采用DNeasy Plant Mini 试剂盒(Qiagen，德国) 提取总DNA，使用NanoDrop8000分光光度计及琼脂糖凝胶电泳检测DNA浓度与质量。DNA质量应满足：1)1.8

图1 试验样品采集示意图Fig.1 The presentation of test sample collecting

用10 U同裂酶MspI和HpaII分别与10 UEcoR I(NEB，美国)组合酶切基因组DNA，37°C 酶切12 h。EcoR I/HpaII组合的酶切产物80°C失活20 min，EcoR I/MspI酶切产物65°C失活20 min。通过T4 DNA连接酶连上接头，连接体系为：2 μL 10×T4 DNA 连接酶Buffer, 10 μL酶切产物, 1μL T4 DNA 酶(400U·μL-1), 1 μLEcoR I 接头(5 pmol·μL-1), 1 μLHpaII/MspI接头(50 pmol·μL-1)，5 μL H2O。16°C 连接16 h，65°C失活20 min(表1)。

预扩增的反应体系为：exTaq10×buffer 2 μL，exTaq DNA聚合酶(U·μL-1) 0.2 μL，dNTP(2.5 mmol·L-1) 1.6 μL，EcoR I+0预扩增引物0.5 μL，HpaII/MspI+0预扩增引物0.5 μL(表1)，连接产物4 μL，灭菌水11.2 μL，总体积为20 μL。

1.2 选择性扩增体系的优化

稀释20倍的预扩增产物用于选择性扩增。选用E3-HM25为选择性扩增引物组合，优化DNA模板用量以及引物、Mg2+、dNTP和exTaq DNA聚合酶的浓度，设置梯度为：1)DNA模板100 ng，200 ng，300 ng以及400 ng； 2)引物0.25 pmol·μL-1，0.5 pmol·μL-1，0.75 pmol·μL-1，1 pmol·μL-1；3)Mg2+1 mmol·L-1，1.5 mmol·L-1，2 mmol·L-1，2.5 mmol·L-1；4)dNTP 0.1 mmol·L-1，0.2 mmol·L-1，0.3 mmol·L-1，0.4 mmol·L-1；5)exTaq DNA聚合酶0.05 U·μL-1，0.125 U·μL-1，0.375 U·μL-1，0.5 U·μL-1。利用E3-HM13、E3-HM33以及E8-HM12三对引物来验证优化选择性扩增体系的稳定性。

预扩增及选择性扩增程序均采用吕金辉的方法[10]。选择性扩增产物使用仪器ABI3730进行检测(由北京睿博兴科生物技术有限公司完成)，毛细管电泳图谱应用GenemarkerV2.2.0软件获得。

表1 接头以及引物序列

Table 1 The DNA sequence of adaptors and primers

接头/引物Adaptors/primers序列Sequences荧光修饰Fluorescentmodification接头AdaptorsEcoRI-adapterFCTCGTAGACTGCGTACC无EcoRI-adapterRAATTGGTACGCAGTCTAC无HpaII/MspI-adapterFGACGATGAGTCTCGAT无HpaII/MspI-adapterRCGATCGAGACTCAT无预扩Pre-amplificationEcoRI+0GACTGCGTACCAATTC无HpaII/MspI+0ATGAGTCTCGATCGG无选扩SelectiveamplificationE3-AACGACTGCGTACCAATTCAAC5’-FAME8-ACCGACTGCGTACCAATTCACC5’-FAMH/M13-CAAATGAGTCTCGATCGGCAA无H/M12-CATATGAGTCTCGATCGGCAT无H/M25-AGTATGAGTCTCGATCGGAGT无H/M33-CTTATGAGTCTCGATCGGCTT无

2 结果与分析

2.1 DNA模板用量对毛细管电泳结果的影响

模板DNA用量对毛细管电泳结果影响不明显。由图2可知，四种DNA模板用量(100 ng、200 ng、300 ng以及400 ng)得到毛细管电泳图谱在信号强度、信号均一性和杂峰干扰程度等方面没有显著差别。为了减少模板DNA的使用量，选择DNA模板用量为100 ng。

2.2 引物浓度对毛细管电泳结果的影响

引物浓度是影响毛细管电泳结果的关键因素(图3)，浓度过高时，片段的数量明显减少。引物浓度0.25 pmol·μL-1、0.5 pmol·μL-1和0.75 pmol·μL-1的毛细管电泳结果相似，但以0.75 pmol·μL-1时毛细管电泳的信号强度最高；而当引物浓度为1 pmol·μL-1时EcoR I/HpaII酶组合及EcoR I/MspI酶组合毛细管电泳结果虽然在200 bp和400 bp的主峰信号强度很强，但在400～600 bp片段的信号强度却非常弱。因此，最适引物浓度为0.75 pmol·μL-1。

H代表EcoR I/Hpa II酶切组合的毛细管电泳结果；M代表EcoR I/Msp I酶切组合的毛细管电泳结果(下同)，Ⅰ, Ⅱ, Ⅲ和Ⅳ分别代表DNA模板量为100 ng、200 ng、300 ng以及400 ng时的毛细管电泳结果。H represented the results of capillary electrophoresis by the endonuclease combination of EcoR I and Hpa II and M represented the results of capillary electrophoresis by the endonuclease combination of EcoR I and Msp I, which were the same meanings in the following figures; Ⅰ, Ⅱ, Ⅲ and Ⅳ stood for the results of capillary electrophoresis for 100 ng, 200 ng, 300 ng and 400 ng DNA template. 图2 不同DNA模板用量的毛细管电泳结果Fig. 2 The results of capillary electrophoresis for different DNA template usage

Ⅰ、Ⅱ、Ⅲ和Ⅳ分别代表引物浓度为0.25 pmol·μL-1， 0.5 pmol·μL-1，0.75 pmol·μL-1，1 pmol·μL-1时的毛细管电泳结果。Ⅰ, Ⅱ, Ⅲ and Ⅳ represented the results of capillary electrophoresis for 0.25 pmol·μL-1, 0.5 pmol·μL-1, 0.75 pmol·μL-1,1 pmol·μL-1 primer concentration respectively.图3 不同引物浓度的毛细管电泳结果Fig. 3 The results of capillary electrophoresis for different primer concentration

Ⅰ、Ⅱ、Ⅲ和Ⅳ分别代表Mg2+浓度为1 mmol·L-1，1.5 mmol·L-1，2 mmol·L-1，2.5 mmol·L-1时的毛细管电泳结果。Ⅰ, Ⅱ, Ⅲ and Ⅳ represented the results of capillary electrophoresis for 1 mmol·L-1, 1.5 mmol·L-1, 2 mmol·L-1, 2.5 mmol·L-1 Mg2+ concentration respectively.图4 不同Mg2+ 浓度的毛细管电泳结果Fig.4 The results of capillary electrophoresis for different Mg2+ concentration

2.3 Mg2+浓度对荧光检测结果的影响

选择性扩增体系中Mg2+浓度是影响毛细管电泳结果的关键因素(图4)。当Mg2+浓度为1 mmol·L-1时，EcoR I/HpaII酶组合及EcoR I/MspI酶组合的毛细管电泳结果主要集中在几个主带上；而当Mg2+浓度为2.5 mmol·L-1时，EcoR I/HpaII酶组合在400 bp的一个主峰信号很强，其余信号相对较弱，EcoR I/MspI酶组合杂峰相对较多；当Mg2+浓度为1.5 mmol·L-1或2 mmol·L-1时，均获得了非常理想的毛细管电泳结果，杂峰少，信号响度均一。可见，Mg2+浓度在1.5 mmol·L-1和2 mmol·L-1时毛细管电泳的效果最佳。

2.4 dNTP浓度对毛细管电泳结果的影响

选择性扩增体系中dNTP浓度是影响毛细管电泳结果的关键因素(图5)。当dNTP浓度为0.1 mmol·L-1时，EcoR I/HpaII酶组合及EcoR I/MspI酶组合的毛细管电泳结果均在两个主峰处表现出很强的信号值，其它信号峰很弱；当dNTP浓度为0.2 mmol·L-1时，EcoR I/MspI酶组合仍然表现在两个主峰处信号很强，其它信号峰很弱；当dNTP浓度为0.3 mmol·L-1及0.4 mmol·L-1时毛细管电泳的效果有了明显的提升，信号值均一，杂峰较少。可见，选择性扩增体系中的dNTP浓度在0.3 mmol·L-1和0.4 mmol·L-1时具有最佳的毛细管电泳效果。

Ⅰ、Ⅱ、Ⅲ和Ⅳ分别代表dNTP浓度为0.1 mmol·L-1, 0.2 mmol·L-1, 0.3 mmol·L-1 0.4 mmol·L-1时的毛细管电泳结果。Ⅰ, Ⅱ, Ⅲ and Ⅳ represented the results of capillary electrophoresis for 0.1 mmol·L-1, 0.2 mmol·L-1, 0.3 mmol·L-1, 0.4 mmol·L-1 dNTP concentration respectively.图5 不同dNTP浓度的毛细管电泳结果Fig. 5 The results of capillary electrophoresis for different dNTP concentration

2.5 exTaq DNA聚合酶浓度对毛细管电泳结果的影响

与普通Taq DNA聚合酶相比，虽然exTaq DNA聚合酶稳定性及准确性更高，但其成本也相对高很多。四种exTaq DNA聚合酶浓度的毛细管电泳结果(图6)，在信号强度的均一性、杂峰程度上没有显著差异。四种exTaq酶浓度获得的EcoR I/HpaII酶组合的毛细管电泳的信号强度是一致的，但EcoR I/MspI酶组合的毛细管电泳的信号强度表现出随浓度的增加而降低的趋势。出于成本的考量，最终选定exTaq酶浓度为0.05 U·μL-1。

2.6 选择性扩增体系的确立与验证

综上所述，本研究初步确立了基于毛细管电泳的MSAP反应体系：DNA模板用量100 ng，引物 0.75 pmol·μL-1，Mg2+1.5～2.0 mmol·L-1，dNTP 0.3～0.4 mmol·L-1，exTaq DNA聚合酶0.05 U·μL-1。为了进一步优化体系，又在1.5～2 mmol·L-1和dNTP 0.3～0.4 mmol·L-1之间设置了五个梯度，两两组合，检测了相应组合获得的毛细管电泳结果的质量。结果表明，Mg2+2.0 mmol·L-1、dNTP 0.375 mmol·L-1的组合具有最佳的条带扩增效果(图7)。选择三对选择性扩增引物(表1)，对优化后选择性扩增体系的稳定性进行了验证。由图8可知，三对引物均获得了亮度清晰、信号较强且相对均一、杂峰较少的毛细管电泳图谱，表明所建的反应体系稳定性较高。

Ⅰ、Ⅱ、Ⅲ和Ⅳ分别代表exTaq DNA聚合酶浓度为0.05 U·L-1, 0.125 U·L-1, 0.375 U·L-1, 0.5 U·L-1时的毛细管电泳结果。Ⅰ, Ⅱ, Ⅲ and Ⅳ represented the results of capillary electrophoresis for 0.05 U·L-1, 0.125 U·L-1, 0.375 U·L-1, 0.5 U·L-1 exTaq DNA polymerase concentration respectively.图6 不同exTaq DNA聚合酶浓度的毛细管电泳结果Fig. 6 The results of capillary electrophoresis for different exTaq DNA polymerase concentration

图7 优化后的毛细管电泳结果Fig. 7 The results of capillary electrophoresis for optimized systems

EX-HMX-H/M中EX-HMX代表选择性扩增使用的引物组合The ‘EX-HMX’ in the presentation of ‘EX-HMX-H/M’ presented the primer combinations in selective amplification.图8 优化体系稳定性验证Fig.8 The verification for the optimized systems

3 讨论

开展DNA甲基化研究对于解析表型变异的表观遗传机理具有重要意义[11-13]。MSAP检测技术是目前应用最为广泛的DNA甲基化检测技术[14]。相对于高通量的DNA甲基化检测方法，MSAP检测技术具有成本低，操作简单以及无需已经参考基因组的特点。基于毛细管电泳的MSAP检测技术较传统基于聚丙烯酰胺凝胶电泳的MSAP检测技术而言，省去了大量的制胶、跑板、染色、显影以及读带等繁琐的过程，极大地缩短了DNA甲基化检测的时间，显著地提高了效率，便于开展大样本的表观遗传变异研究，并且能够更加灵敏地识别浓度很低的扩增片段。本研究综合考虑了影响毛细管电泳扩增效果的多个因素，确立了基于毛细管电泳的日本落叶松MSAP技术最佳反应体系，为大批量分析日本落叶松DNA甲基化奠定了基础。

研究发现Mg2+及dNTP浓度是影响毛细管电泳结果最为关键的因素，浓度过低或过高均会导致毛细管电泳图谱条带明显减少，杂峰数量显著增多。该结果与宋春艳[15]、吕金辉[10]基于毛细管电泳的AFLP条件优化结果一致，同时也与李金龙[16]在北京油鸡中报道一致：Mg2+最佳浓度为1.25 mmol·L-1，浓度过低条带较少，过高时杂峰比较多；dNTP最佳浓度为 0.15 mmol·L-1，浓度过高使得条带减少。本研究最佳反应体系中Mg2+浓度为2 mmol·L-1，dNTP浓度为0.375 mmol·L-1，这可能是由于物种不同，酶切获得的片段多少不同造成的。

引物浓度也是影响毛细管电泳结果的重要因素，当引物浓度为1 pmol·μL-1时，毛细管图谱的信号会集中在几个主峰上。吕金辉[10]和李金龙[16]研究发现引物浓度是影响毛细管电泳结果的关键因素，过高会造成非特异性扩增，过低会造成片段的缺失。

DNA模板用量和exTaq DNA聚合酶浓度对毛细管电泳结果影响较小。为了降低成本，选定DNA模板用量100 ng和exTaq DNA聚合酶浓度0.05 U·μL-1为最佳配置，可获得较好的条带扩增效果。吕金辉[10]认为DNA模板用量会影响基于毛细管电泳的AFLP结果。我们推测差异的原因可能在于物种的差异或者是DNA模板用量梯度差异设置造成的。同时，吕金辉[10]也发现exTaq DNA聚合酶的使用量不是影响毛细管电泳结果的关键因素，这与我们的结果是一致的。

4 结论

本文以日本落叶松未成熟木质部为试材，探索不同DNA模板用量及选择扩增体系中引物浓度、Mg2+浓度、dNTP浓度和exTaq DNA聚合酶浓度等因素对毛细管电泳效果的影响，从而确立最适的MSAP反应体系。结果表明，DNA模板以及exTaq DNA聚合酶的用量对毛细管电泳结果的影响不明显，过高的引物浓度会减少扩增的条带数目，过低或过高的Mg2+和dNTP浓度会导致毛细管电泳图谱条带明显减少，杂峰显著增多。最终确立了最佳反应体系为：DNA 模板100 ng，引物0.75 pmol·μL-1，Mg2+2 mmol·L-1，dNTP 0.375 mmol·L-1和exTaq DNA聚合酶 0.05 U·μL-1。通过三对选择性扩增引物的验证性实验表明，该体系可获得亮度清晰、信号较强且强度相对均一以及杂峰较少的毛细管电泳图谱。本文首次建立了基于毛细管电泳的日本落叶松最佳MSAP反应体系，为今后开展日本落叶松表观遗传变异相关研究奠定基础。

[1] 谭建新, 孙玉洁. 表观基因组学研究方法进展与评价[J]. 遗传, 2009, 31(1):3-12.

[2] Cokus SJ, Feng S, Zhang X,etal. Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning[J]. Nature, 2008, 452(7184):215-219.

[3] Yang H, Chang F, You C,etal. Whole-genome DNA methylation patterns and complex associations with gene structure and expression during flower development in Arabidopsis[J]. The Plant Journal: for Cell and Molecular Biology, 2015, 81(2):268-281.

[4] Xing MQ, Zhang YJ, Zhou SR,etal. Global Analysis Reveals the Crucial Roles of DNA Methylation during Rice Seed Development[J]. Plant Physiology, 2015, 168(4):1417-1432.

[5] Lira-Medeiros CF, Parisod C, Fernandes RA,etal. Epigenetic variation in mangrove plants occurring in contrasting natural environment[J]. PloS One, 2010, 5(4):e10326.

[6] 沈 圣, 屈彦纯, 张 军. 下一代测序技术在表观遗传学研究中的重要应用及进展[J]. 遗传, 2014, 36(3):256-275.

[7] Schulz B, Eckstein RL, Durka W. Scoring and analysis of methylation-sensitive amplification polymorphisms for epigenetic population studies[J]. Molecular Ecology Resources, 2013, 13(4):642-653.

[8] Lai M, Sun X, Chen D,etal. Age-related trends in genetic parameters for Larix kaempferi and their implications for early selection[J]. BMC Genetics, 2014, 15 Suppl 1:S10.

[9] Li A, Song W-Q, Chen C-B,etal. DNA methylation status is associated with the formation of heterosis in Larix kaempferi intraspecific hybrids[J]. Molecular Breeding, 2013, 31(2):463-475.

[10] 吕金辉, 胡建军, 卢孟柱. 基于毛细管电泳的柳树 AFLP 分子标记研究[J]. 北京林业大学学报, 2012, 34(1):51-57.

[11] Ci D, Song Y, Du Q,etal. Variation in genomic methylation in natural populations of Populus simonii is associated with leaf shape and photosynthetic traits[J]. Journal of Experimental Botany, 2016, 67(3):723-737.

[12] Latzel V, Allan E, Bortolini Silveira A,etal. Epigenetic diversity increases the productivity and stability of plant populations[J]. Nature Communications, 2013, 4:2875.

[13] Medrano M, Herrera CM, Bazaga P. Epigenetic variation predicts regional and local intraspecific functional diversity in a perennial herb[J]. Molecular Ecology, 2014, 23(20):4926-4938.

[14] Alonso C, Pérez R, Bazaga P,etal. MSAP markers and global cytosine methylation in plants: a literature survey and comparative analysis for a wild growing species[J]. Molecular Ecology Resources, 2016, 16(1):80-90.

[15] 宋春艳, 楚 鹰, 陈明杰, 等. 基于 CEQ 遗传分析仪的 AFLP 条件的优化试验[J]. 食用菌, 2006, 28(3):14-15.

[16] 李金龙, 唐韶青, 邹智元, 等. 北京油鸡 MSAP 毛细管电泳荧光检测技术的建立[J]. 遗传, 2014, 36(5):495-502.

(责任编辑：张 研)

Establishment ofLarixkaempferiMSAP System Based on Capillary Electrophoresis

DIAO Shu, SUN Chao, SUN Xiao-mei, ZHANG Shou-gong

(Research Institute of Forestry, Chinese Academy of Forestry, Key Laboratory of Tree Breeding and Cultivation, State Forestry Administration, Beijing 100091, China)

[Objective]This study aims at constructing the excellent reaction system of MSAP technology and obtaining DNA methylation maps with clear bands and uniform signal intensity, so as to provide references for the study on genetic variation ofLarixkaempferiin the future.[Method] Taking immature xylem as experimental materials, the impact of five factors (dosage of DNA template, concentration of primers, Mg2+, dNTP and exTaq DNA polymerase) on the result of capillary electrophoresis was investigated. [Result] It was found that different concentrations of DNA template and the exTaq DNA polymerase did not affect the result of capillary electrophoresis. However, high concentration of primers would lead to sparse fragments in MSAP profile. Mg2+concentration and dNTP concentration were the most important factors which determined the result of capillary electrophoresis. When the concentration of Mg2+/dNTP was unsuitable, the amount of band would significantly decrease, and the amount of impurity peak would significantly increase. [Conclusion] The MASP system determined for capillary electrophoresis inL.kaempferiis: DNA template dosage 100 ng, primers 0.75 pmol·μL-1, Mg2+2 mmol·L-1, dNTP 0.375 mmol·L-1and exTaq DNA polymerase 0.05 U·μL-1. It can be widely used in detecting DNA methylation.

DNA methylation; methylation sensitive amplification polymorphism (MSAP); capillary electrophoresis; system establishment;Larixkaempferi

10.13275/j.cnki.lykxyj.2017.03.001

2016-09-21

林业公益性行业科研专项经费重点项目(20150414)。

刁姝(1987—)，女，博士生，研究方向为林木遗传育种。E-mail:diaoshu0802@163.com

* 责任作者：张守攻(1957—)，研究员，研究方向为针叶树遗传改良。Tel:010-62889685;E-mail:larch_rif@163.com

S791.223

A

1001-1498(2017)03-0361-07