不同抗性短枝木麻黃種源苗木接種青枯病菌后酚類物質含量的變化

魏永成, 張勇, 仲崇祿, 孟景祥, 陳珍

?

不同抗性短枝木麻黃種源苗木接種青枯病菌后酚類物質含量的變化

魏永成, 張勇, 仲崇祿*, 孟景祥, 陳珍

(中國林業科學研究院熱帶林業研究所，廣州 510520)

為了解短枝木麻黃()抗青枯病的機理，對接種青枯病菌()后短枝木麻黃的單寧和黃酮含量變化進行了分析。結果表明，不同抗性短枝木麻黃種源小枝的總酚和單寧含量呈現不同的變化趨勢，高抗、中抗種源均呈現先升高后降低的變化趨勢，峰值均約為126 mg g–1，但中抗種源的峰值出現時間較晚，而易感種源則呈逐漸升高趨勢。抗、感種源木麻黃接種青枯菌后，小枝中縮合單寧含量均呈現逐漸升高的趨勢，但高抗種源的縮合單寧含量均顯著高于易感種源，增加70.33%。抗性種源黃酮含量呈S型上升趨勢，易感種源則持續緩慢升高。這表明接種青枯病菌后，抗、感短枝木麻黃種源表現出不同的防御特征，次生物質含量增幅越大，抑菌抗氧化能力越強，短枝木麻黃表現出的抗性越強。

短枝木麻黃；青枯病菌；接種；單寧；黃酮

木麻黃抗風、耐瘠薄，是我國華南地區營建沿海防護林的先鋒樹種，在一些立地條件較差的前緣沙質地帶，木麻黃甚至是唯一可造林的樹種[1]，具有顯著的生態、經濟價值[2–3]。青枯病()是一種危害嚴重、傳播廣泛的毀滅性土傳病害, 該病傳播速度快、發病死亡率高，可危害54科200余種植物，包括許多具有重要經濟和社會價值的植物, 如煙草()、番茄()、馬鈴薯()、香蕉()和桉樹(sp.)等[4]，是世界范圍內最難防治的細菌性重大病害之一[5]。青枯病是木麻黃的重要病害之一，嚴重影響我國華南沿海防護林的生態安全。青枯病菌作為一種土傳細菌，具有傳播速度快、變異大等特點，人們在化學防治和生物防治方面做了大量工作，但依然無法改變青枯病不斷蔓延的趨勢，如何有效防治青枯病已成為影響木麻黃發展的重要因素。植物在長期進化過程中形成的自身抗性是控制青枯病擴散的重要策略[6]。在影響植物抗病性的因素中，次生代謝物質是植物產生抗病性的重要生理生化基礎，能夠有效抑制病原菌在寄主體內的擴散與增殖，甚至直接殺滅病原菌[7]。

酚類物質是植物重要的次生代謝物質，參與許多生理生化過程，許多研究表明酚類物質在植物抵御病原物侵襲過程中發揮了重要作用[8]。受病原物侵襲后，可誘導植物產生大量可溶性酚類物質，抑制病原物的遷移和擴增[9]。細胞壁上的大量酚類物質也與植物的抗病性密切相關[10]。黃酮類化合物也可以抑制多酚氧化酶、細胞色素氧化酶等的活性，進而影響植物的抗病性[11–13]。綜上，植物酚類物質和黃酮類化合物與植物抗性有密切關系，研究植物受到病原菌侵襲后酚類物質、黃酮類化合物的含量變化規律，可以在一定程度上反映植株的抗病特征。

本研究以選育的抗青枯病短枝木麻黃()高抗(R)、中抗(M)和易感(S)為材料，探討感染青枯病后小枝中酚類物質和黃酮含量的變化規律，揭示短枝木麻黃對青枯病的抗性生理機理, 明確短枝木麻黃抗病與感病種源間的差異特征，為短枝木麻黃抗青枯病品種選育提供科學依據。

1 材料和方法

1.1 材料

于2017年5月在廣東省徐聞縣的種源試驗林(2014年造林)中收集短枝木麻黃()材料，包括3個抗性種源(R、M和S)，每種源選擇50株，每株取10支幼嫩小枝，水培生根后移至營養杯培養。供試菌株為中國林科院熱帶林業研究所前期獲得的高致病力青枯病菌() GL-3菌株[14]。

1.2 接種和取樣

參考劉勇等[15]的方法，將GL-3菌株在TTC (2,2,3-三苯基四唑氯)培養基上活化48 h，選擇單菌落接種到普通營養培養基中擴大培養24 h，然后接種到液體培養基中30℃搖床24 h。培養液靜置30 min收集菌體，用蒸餾水將菌液稀釋成3×109cfu mL–1用于接種。

2018年3月，大部分幼苗高度在40 cm以上，對幼苗進行接種處理。接種時將所有幼苗從土中取出,小心去掉土壤盡量避免幼苗根部損傷，用水洗凈根部。剪掉約1/3根系，一半植株浸入30℃細菌懸浮液中(I, inoculation) 30 min，另外一半植株浸入蒸餾水作為對照(W)。將接種后的植株移栽至裝有營養土(黃土∶泥炭土∶沙子=1∶1∶1)的塑料盆(長×寬×高= 40 cm×25 cm×8 cm)中，每盆20株。共6種組合(R-I、R-W、M-I、M-W、S-I、S-W)，采用隨機完全區組設計，每組合處理苗木4盆80株，重復3次，共12盆240株。接種后每天澆水，相對濕度控制在80%左右。

接種當天取中上部的功能小枝，以后每天固定時間取樣1次，直至有植株完全枯萎死亡，共取樣8次。每次取樣時每組合隨機選5株，-20℃冰箱保存備用。

1.3 指標測定

按王翠穎等[16–17]的方法提取和測定單寧和黃酮含量，總酚含量利用Folin-Ciocalteu試劑比色法測定[18]，縮合單寧用丁醇/鹽酸(95∶5)法測定[19]。3次重復，取平均值。

1.4 數據分析

利用SPSS 19.0和EXCEL 2007對數據進行統計分析，以均值±標準誤表示。采用Duncan多重比較法對數據進行方差分析。

2 結果和分析

2.1 總酚含量

接種青枯菌后，不同抗性短枝木麻黃種源苗木小枝中的總酚含量呈現不同的變化趨勢(表1)。高抗、中抗種源均呈現先升高后降低的趨勢，但中抗種源的峰值出現時間較晚；易感種源則呈逐漸升高趨勢，在感病末期雖略有下降，但整體水平較高。高抗短枝木麻黃接種青枯菌后，總酚含量迅速升高，在接種后2 d達最大值126.89 mg g–1，比對照高47.20%。隨后迅速降低，接種后7 d與對照的差異不顯著。中抗種源接種后，總酚含量也會迅速升高, 但直到接種后4 d才達最大值126.16 mg g–1, 比對照高39.22%，隨后迅速降低，接種后5~7 d與對照差異不顯著。易感種源受到青枯菌侵襲后，總酚含量持續上升，最大值(109.57 mg g–1)出現在第5天，接種后6 d時略微下降，在接種后7 d比對照高32.53%。

接種青枯菌后1 d，高抗、中抗、易感短枝木麻黃間的總酚含量差異不顯著。接種2 d后，總酚含量依次為高抗>中抗>易感，高抗與易感種源間的差異顯著。接種3~4 d后，中抗種源的總酚含量高于其余種源。接種后5~7 d，3種源間的差異不顯著，且易感種源的總酚含量維持在較高濃度。

2.2 單寧含量

不同抗性種源短枝木麻黃接種青枯菌后，小枝中單寧含量的變化趨勢與總酚含量的一致(表2)。高抗、中抗種源單寧含量呈現明顯的先升后降的變化趨勢，易感種源則持續升高，在接種后6~7 d略微下降。高抗種源接種青枯菌2 d后的單寧含量達最大(125.92 mg g–1)，比對照高47.05%。中抗種源在接種后4 d達最大(125.41 mg g–1)，比對照高39.58%。易感種源接種青枯菌5 d后達最大值，單寧含量比對照高19.56%。

表1 短枝木麻黃接種青枯菌后總酚含量(mg g–1)

R: 高抗種源; M: 中抗種源; S: 易感種源; I: 接種處理; W: 對照; 數據后不同字母表示差異顯著(<0.05)。下同。

R: High resistant provenance; M: Media resistant provenance; S: Susceptible provenance; I: Inoculation; W: Control; Data followed different letters indicate significant difference at 0.05 level. The same is following Tables.

接種青枯菌后1 d，高抗、中抗和易感短枝木麻黃間的單寧含量差異不顯著。接種2 d后的單寧含量依次為高抗>中抗>易感，且差異顯著。接種4 d后，中抗種源的單寧含量顯著高于高抗和易感。接種5~7 d后，3種源間的差異不顯著，易感種源的單寧含量維持在較高水平。這表明高抗短枝木麻黃受到青枯菌侵襲后，能迅速大量合成單寧并有效發揮抗菌作用，而易感種源的單寧合成反應緩慢且效用不明顯。

2.3 縮合單寧含量

不同抗性種源短枝木麻黃接種青枯菌后，小枝的縮合單寧含量均呈逐漸升高的變化趨勢(表3)。接種青枯菌后7 d，高抗種源的縮合單寧含量達最大值13.03 mg g–1，比對照高70.33%，中抗和易感種源分別為10.36和10.2 mg g–1，分別比對照高19.35%和42.06%。高抗種源顯著高于易感種源，這表明高抗、易感短枝木麻黃受到青枯菌侵襲后，均會誘導產生縮合單寧，但高抗種源的縮合單寧含量更高, 導致抵御青枯菌侵襲的效果不同。

2.4 黃酮含量的變化

不同抗性種源短枝木麻黃接種青枯菌后，小枝的黃酮含量呈逐漸升高的變化趨勢(表4)。其中，高抗、中抗種源的黃酮含量呈S型上升趨勢，易感種源則持續緩慢升高。高抗種源接種青枯菌后5 d的黃酮含量達最大值56.66 mg g–1，比對照高35.68%，中抗種源在接種后7 d達最大值54.08 mg g–1，比對照高27.49%，而易感種源接種青枯菌7 d后，黃酮含量比對照高14.54%。

表2 短枝木麻黃接種青枯菌后單寧含量(mg g–1)

表3 短枝木麻黃接種青枯菌后縮合單寧含量(mg g–1)

表4 短枝木麻黃接種青枯菌后黃酮含量(mg g–1)

接種青枯菌后1~2 d，高抗、中抗、易感種源間的黃酮含量差異不顯著。接種后3~7 d，高抗種源的黃酮含量為51~57 mg g–1，均顯著高于易感種源。這表明高抗短枝木麻黃受到青枯菌侵襲后，大量合成黃酮，起到抑制青枯菌的作用。

3 討論

木麻黃受到青枯病菌侵襲后，體內生理生化物質會發生系統性變化。已有研究主要集中在少數或者單一無性系，極大地限制了木麻黃抗病育種工作的進展。同時由于病原菌不斷發生變異，抗病品種選育和應用相對滯后，使得許多木麻黃林大面積染病死亡，急需新的抗病種質材料充實到沿海防護林帶。本研究從種源層次解析短枝木麻黃抗性機制, 擴大了遺傳資源篩選范圍，對短枝木麻黃抗病選育提供了良好的理論基礎。

林木抗病性是在長期演化過程中所形成的生物特性，根據抗性產生機制可以分為組成型抗性和誘導型抗性[20]。組成型抗性是植物不論受侵襲與否均產生與抗性有關的物質，誘導型抗性是植物所產生的一種更為主動的抗性反應方式，即在受到如病原菌等物的侵襲后，植物誘導產生一些物質來發揮防御作用[21]。本研究中，受到青枯菌侵襲后，高抗和易感短枝木麻黃的單寧和黃酮含量均會升高，高抗種源的單寧含量呈先升后降的變化趨勢，黃酮含量呈S型增長趨勢，而易感種源的單寧和黃酮含量均呈逐漸升高的變化趨勢，但高抗種源中單寧和黃酮含量增幅高于易感種源，這表明短枝木麻黃受青枯菌侵襲后表現出誘導型抗性。因此，誘導產生的單寧和黃酮含量可以作為評價短枝木麻黃抗青枯病的生理指標。同時，不同種源單寧和黃酮含量峰值出現時間并不一致，為快速有效地鑒定短枝木麻黃抗病與否，需要進一步明確各抗性種質材料合成單寧和黃酮的閾值。

單寧和黃酮類化合物是多酚中的代表性成分,是短枝木麻黃受到病原菌侵襲后產生的兩類重要次生代謝物質，來源于莽草酸途徑和苯丙氨酸代謝途徑，具有抑菌、抗氧化等多種功能，是植物體內主要的防御措施[22]。王翠穎等[16]報道抗病木麻黃無性系在接種青枯病菌前后的單寧含量均高于感病無性系。林秀琴等[23]在篩選木麻黃抗星天牛品系時，認為抗蟲等級越高，植物體內的單寧和黃酮等次生代謝物含量越高。這與本研究結果基本一致, 但本研究中不同種源間短枝木麻黃的單寧和黃酮含量在接種當天或未接種時的差異不顯著，這可能是由于不同種源短枝木麻黃遺傳特性造成的。而受侵襲后抗感種源的單寧含量變化趨勢不同則是因為某些基因差異表達的結果。

由于短枝木麻黃受到青枯病菌侵襲后會出現較長時間的潛伏期，因此本研究中采用高劑量菌液接種，以檢測短枝木麻黃感病后的生理生化變化。但在野外，由于降雨充沛、人工干擾等因素影響, 短枝木麻黃并不會表現出枯萎死亡的病癥，甚至出現部分枯萎植株逐漸恢復的情況。短枝木麻黃中的單寧和黃酮有組成型和誘導型兩類，其中誘導型的單寧和黃酮主要是受到病原菌侵襲后產生的，二者的變化規律能夠在一定程度上反映植物的抗性機理, 進一步研究其他防御物質的變化規律將有利于全面揭示抗性短枝木麻黃對青枯病的抗性機制。

[1] Zhong C L, Zhang Y, Chen Y, et al. Casuarina research and applications in China [J]. Symbiosis, 2010, 50(1/2): 107–114. doi: 10. 1007/s13199-009-0039-5.

[2] ZHONG C L, MANSOUR S, NAMBIAR-VEETIL M, et al.: A model tree for basic research in actinorhizal symbiosis [J]. J Bioscience, 2013, 38(4): 815–823. doi: 10.1007/s12038-013-9370-3.

[3] ZHANG Y, ZHONG C L, HAN Q, et al. Reproductive biology and breeding system in(Casuarinaceae): Impli- cation for genetic improvement [J]. Aust J Bot, 2016, 64(2): 120–128. doi: 10.1071/BT15184.

[4] QIAO J Q, CHEN Z Y, LIU Y Z, et al. Research progress on bacterial wilt of nightshade family [J]. Acta Phytopathol Sin, 2013, 43(1): 1–10. doi: 10.3969/j.issn.0412-0914.2013.01.001.喬俊卿, 陳志誼, 劉郵洲, 等. 茄科作物青枯病研究進展 [J]. 植物病理學報, 2013, 43(1): 1–10. doi: 10.3969/j.issn.0412-0914.2013.01.001.

[5] SALANOUBAT M, GENIN S, ARTIGUENAVE F, et al. Genome sequence of the plant pathogen[J]. Nature, 2002, 415(6871): 497–502. doi: 10.1038/415497a.

[6] HUANG G H, ZHONG C L, SU X H, et al. Genetic variation and structure of native and introduced(L. Johnson) provenances [J]. Silv Genet, 2009, 58(1/2/3/4/5/6): 79–85. doi: 10. 1515/sg-2009-0010.

[7] SUN T H, LI J, WANG T, et al. Difference in protein expression of apple leaves with different resistance response tof. sp.infection [J]. Acta Hort Sin, 2018, 45(3): 409–420. doi: 10.16420/j.issn.0513-353x.2017-0686.孫天驊, 李佳, 王濤, 等. 抗病與感病蘋果葉片應答輪紋病菌侵染的蛋白質表達差異分析 [J]. 園藝學報, 2018, 45(3): 409–420. doi: 10.16420/j.issn.0513-353x.2017-0686.

[8] TENG T, CAO F X, WANG M, et al. Effects of pine wood nematode (PWN) infection on the activity of PAL and the content of phenols in the host pine [J]. J CS Univ For Technol, 2007, 27(3): 124–127. doi: 10. 3969/j.issn.1673-923X.2007.03.013.滕濤, 曹福祥, 王猛, 等. 松材線蟲侵染對松樹苯丙氨酸解氨酶及酚類物質的影響 [J]. 中南林業科技大學學報, 2007, 27(3): 124–127. doi: 10.3969/j.issn.1673-923X.2007.03.013.

[9] GUO H L, CHENG G W, CHEN J, et al. The function of phenolic metabolism in resistant mechanism of gray leaf spot of corn [J]. Acta Phytopathol Sin, 2003, 33(4): 342–346. doi: 10.3321/j.issn:0412-0914. 2003.04.012.郭紅蓮, 程根武, 陳捷, 等. 玉米灰斑病抗性反應中酚類物質代謝作用的研究 [J]. 植物病理學報, 2003, 33(4): 342–346. doi: 10.3321/ j.issn:0412-0914.2003.04.012.

[10] CHEN J, GAO H M, ZHU Y G, et al. Effect of phenols and its meta- bolism on activity cell wall degrading enzymes produced by[J]. Acta Phytopathol Sin, 1996, 26(2): 171–176. doi: 10.13926/j.cnki.apps.1996.02.019.陳捷, 高洪敏, 朱有釭, 等. 酚類物質和代謝對瓜果腐霉菌產生的細胞壁降解酶活性的影響 [J]. 植物病理學報, 1996, 26(2): 171–176. doi: 10.13926/j.cnki.apps.1996.02.019.

[11] HU F, KONG C H, XU X H, et al. Inhibitory effect of flavones fromon the major pathogens in citrus orchard [J]. Chin J Appl Ecol, 2002, 13(9): 1166–1168.胡飛, 孔垂華, 徐效華, 等. 勝紅薊黃酮類物質對柑桔園主要病原菌的抑制作用 [J]. 應用生態學報, 2002, 13(9): 1166–1168.

[12] WU Y H, ZHONG L J, ZHAO X X. Relationship among total phenol, flavonoid and PPO activity with the necrosis of tobacco infected by PVYN [J]. Acta Phytopathol Sin, 2007, 37(4): 398–402. doi: 10.3321/ j.issn:0412-0914.2007.04.010.吳元華, 鐘麗娟, 趙秀香. 煙草感染PVY~N后葉脈壞死與總酚、類黃酮及PPO關系研究 [J]. 植物病理學報, 2007, 37(4): 398–402. doi: 10.3321/j.issn:0412-0914.2007.04.010.

[13] HUANG Q X, QU Y Y, YAO Z P, et al. Correlation between fusarium wilt resistance and expression levels of genes involved in flavonoid metabolism pathway inL. [J]. Acta Agron Sin, 2017, 43(12): 1791–1801. doi: 10.3724/SP.J.1006.2017.01791.黃啟秀, 曲延英, 姚正培, 等. 海島棉枯萎病抗性與類黃酮代謝途徑基因表達量的相關性 [J]. 作物學報, 2017, 43(12): 1791–1801. doi: 10.3724/SP.J.1006.2017.01791.

[14] WANG S K, WANG J, XU D P. Comparison of four different DNA extraction methods ofin eucalyptus and sensi- tivity of PCR detection [J]. For Pest Dis, 2007, 26(5): 4–7. doi: 10. 3969/j.issn.1671-0886.2007.05.002.王勝坤, 王軍, 徐大平. 四種桉樹青枯菌DNA提取方法及PCR檢測靈敏度比較 [J]. 中國森林病蟲, 2007, 26(5): 4–7. doi: 10.3969/j. issn.1671-0886.2007.05.002.

[15] LIU Y, QIN X Y, LI W Z, et al. The resistance evaluation to bacterial wilt of tobacco germplasm in Yunnan Province [J]. J Plant Genet Resour, 2010, 11(1): 10–16.劉勇, 秦西云, 李文正, 等. 抗青枯病煙草種質資源在云南省的評價 [J]. 植物遺傳資源學報, 2010, 11(1): 10–16.

[16] WANG C Y, SUN S, WANG J. Relationship between tannin content fromand their resistance to[J]. Shandong Chem Ind, 2015, 44(3): 77,82. doi: 10.3969/j. issn.1008-021X.2015.03.024.王翠穎, 孫思, 王軍. 木麻黃抗感青枯病品系單寧含量差異 [J]. 山東化工, 2015, 44(3): 77, 82. doi: 10.3969/j.issn.1008-021X.2015. 03.024.

[17] WANG C Y, SUN S, ZHOU E X, et al. Purification of total flavonoids fromby macroporous adsorption resin & inhibition of total flavonoids on[J]. Prot For Sci Technol, 2012(1): 11–14,32. doi: 10.3969/j.issn.1005-5215.2012. 01.005.王翠穎, 孫思, 周而勛, 等. 大孔吸附樹脂純化木麻黃總黃酮及對青枯菌的抑制 [J]. 防護林科技, 2012(1): 11–14,32. doi: 10.3969/j. issn.1005-5215.2012.01.005.

[18] MAKKAR H P S, BLUMMEL M, BOROWY N K, et al. Gravimetric determination of tannins and their correlations with chemical and protein precipitation methods [J]. J Sci Food Agric, 2010, 61(2): 161– 165. doi: 10.1002/jsfa.2740610205.

[19] MAKKAR H P S. Quantification of Tannins in Tree and Shrub Foliage: A Laboratory Manual [M]. Netherlands: Springer, 2003: 53–54.

[20] WANG X S, YANG C L, WANG S S, et al. Changes of phenols and lignin contents in alfalfa leaf damaged by[J]. Chin J Appl Ecol, 2014, 25(6): 1688–1692.王小珊, 楊成霖, 王森山, 等. 牛角花齒薊馬為害后苜蓿葉酚類物質和木質素含量的變化 [J]. 應用生態學報, 2014, 25(6): 1688– 1692.

[21] WU Y Q, GUO Q Y. Potential resistance of tannins: Flavoniods in upland cotton against(Hübner) [J]. Acta Ecol Sin, 2001, 21(2): 286–289. doi: 10.3321/j.issn:1000-0933.2001.02.018.武予清, 郭予元. 棉花單寧: 黃酮類化合物對棉鈴蟲的抗性潛力 [J]. 生態學報, 2001, 21(2): 286–289. doi: 10.3321/j.issn:1000-0933. 2001.02.018.

[22] XU X Y, ZHANG W Q, HUANG Y H, et al. Comparison of identification methods of bacterial wilt resistance inand screening of resistant germplasm resources [J]. J S China Agric Univ, 2017, 38(4): 87–94. doi: 10.7671/j.issn.1001-411X.2017.04.015.許秀玉, 張衛強, 黃鈺輝, 等. 木麻黃青枯病抗性鑒定方法比較及抗病種質篩選 [J]. 華南農業大學學報, 2017, 38(4): 87–94. doi: 10. 7671/j.issn.1001-411X.2017.04.015.

[23] LIN X Q, HUANG J S, CAI S P, et al. Correlation of inclusions contents of differentspp. strains and resistibility to[J]. J CS Univ For Technol, 2014, 34(5): 37–41. doi: 10.3969/j.issn.1673-923X.2014.05.008.林秀琴, 黃金水, 蔡守平, 等. 木麻黃內含物含量與對星天牛抗性關系的研究 [J]. 中南林業科技大學學報, 2014, 34(5): 37–41. doi: 10.3969/j.issn.1673-923X.2014.05.008.

Changes in Phenols Contents inProvenances with Different Resistibility Infected by

WEI Yong-cheng, ZHANG Yong, ZHONG Chong-lu*, MENG Jing-xiang, CHEN Zhen

(Research Institute of Tropical Forestry, Chinese Academy of Forestry, Guangzhou 510520, China)

In order to understand the resistance mechanism of casuarinas () against, the changes in contents of tannin and phenols in casuarinas infected bywere analyzed. The results showed that the change trends of tannin and phenols contents in provenances with different resistibility were different. The contents of tannin and phenols in high- and media-resistant provenances increased at first and then decreased with the peak of 126 mg g–1, and the peak of media-resistant provenances came later. The susceptible provenances showed rising trend. After infected with, the condensed tannin content increased gradually in three provenances, and that in the high-resistant provenances were 70.33% higher than that in susceptible provenances. The content of flavone in resistant provenances showed S type rising trend, while that in susceptible provenances increased slowly. Therefore, it was demonstrated that there were different defense mechanisms between resistant and susceptible provenances. The greater the content of secondary compounds increase, the stronger bacteriostasis and antioxidant activities,would show higher resistance to bacterial wilt.

;; Inoculate; Tannin; Flavone

10.11926/jtsb.3983

2018-08-14

2018-09-21

中央級公益性科研院所基本科研業務費專項資金(CAFYBB2017QA007)；廣東省林業科技創新項目(2014KJCX017)資助

This work was supported by the Special Project for Basic Scientific Research of Central Public Welfare Institution (Grant No. CAFYBB2017QA007), and the Project for Forestry Science and Technology Innovation of Guangdong Province of China (Grant No. 2014KJCX017).

魏永成(1988~ )，男，博士，助理研究員，主要從事林木遺傳改良研究。E-mail: yongcheng369@163.com

通信作者 Corresponding author.E-mail: zclritf@126.com