吳東山，梁健敏，楊章旗，黃永利

（1. 廣西壯族自治區(qū)林業(yè)科學研究院，國家林業(yè)局馬尾松工程技術(shù)研究中心，廣西馬尾松工程技術(shù)研究中心，廣西 南寧，530002；2. 廣西貴港市覃塘林場，廣西 貴港，537121；3. 廣西南寧市林業(yè)科學研究所，廣西 武鳴 530100 ）

馬尾松Pinus massoniana是我國南方分布面積最大的針葉樹種，也是我國重要的速生用材及紙漿材樹種。馬尾松作為一種特殊的經(jīng)濟樹種，還可以生產(chǎn)另外一個重要的高經(jīng)濟價值產(chǎn)品——松脂。松脂是香精、香料、造紙等重要的工業(yè)原料，我國90%的松脂來自于馬尾松。從上世紀70年代就進行了馬尾松[1]、濕地松P. elliottii[2-6]，思茅松P. kesiya var. langbianensis[7]等主要采脂樹種產(chǎn)脂性狀遺傳研究，取得了一定的成就。但上世紀80年代營建的種子園多以用材為主要選育目標，對馬尾松產(chǎn)脂性狀選育的重視程度不夠。目前不論是生長性狀還是產(chǎn)脂性狀的研究，大多數(shù)仍以單一性狀選擇為主，這導致在進行半同胞子代家系測定與種子園的去劣疏伐時損失掉部分產(chǎn)脂性狀優(yōu)良的家系及無性系。因此，實現(xiàn)馬尾松材用和產(chǎn)脂性狀的同步改良，對兩者進行聯(lián)合選擇，達到高材、高產(chǎn)脂優(yōu)良性狀于一體的目的，滿足生產(chǎn)需求。

研究表明，根據(jù)林木不同的選育目標運用針對性的評價方法能夠選出性狀遺傳穩(wěn)定性優(yōu)良的家系[8-14]。目前，選擇產(chǎn)脂及生長兼優(yōu)的家系主要的評價方法有遺傳主成分、遺傳距離聚類、育種值綜合評分法及指數(shù)選擇法等。比較這4種方法在油松P. tabuliformis子代生長測定，結(jié)果表明指數(shù)選擇法被認為是評價多性狀優(yōu)良家系理想的方法[15]。綜合評價方法在馬尾松等松樹多性狀選擇上已有不少的應(yīng)用，但在馬尾松半同胞家系子代產(chǎn)脂性狀的聯(lián)合選擇研究中還未進行應(yīng)用。本研究對20年生馬尾松的53個半同胞家系子代采脂林生長及產(chǎn)脂性狀進行遺傳和相關(guān)分析，應(yīng)用性狀指數(shù)選擇法進行優(yōu)良家系選擇，以探討馬尾松半同胞家系多性狀綜合選擇的可行性，為更合理選擇利用優(yōu)良家系材料和種子園的去劣疏伐提供科學依據(jù)。

1 材料與方法

1.1 試驗材料

馬尾松子代測定林位于廣西壯族自治區(qū)南寧市林業(yè)科學研究所國家馬尾松良種基地，23°10′05′ N，108°00′01′ E。屬熱帶北緣季風氣候，干濕季節(jié)明顯，年平均氣溫21.5℃，1月平均氣溫12.5℃，年平均有霜日23 d，年降水量1 246 mm。土壤為第四紀紅土發(fā)育而成的中壤質(zhì)厚層赤紅壤，土層厚1 m以上，pH值4.5～5.0。試驗地海拔120 m，地勢平坦。1981年營建馬尾松無性系種子園，面積57 hm2，1992年進行半同胞子代家系造林，包含對照，均為同齡1年生苗自由授粉子代家系。2011年進行子代家系的采脂試驗。共53個參試號（進行采脂試驗時試驗林有69個家系，編號為1～69，去除達不到數(shù)量株數(shù)家系后共有53個參試號進行試驗），其中48個家系來自該基地馬尾松無性系種子園半同胞子代，另外5個為對照處理（包括本地馬尾松CK1，古蓬馬尾松CK2、桐棉馬尾松CK3，信宜馬尾松CK4，貴州馬尾松CK5）。試驗設(shè)計采用完全隨機區(qū)組，設(shè)6個重復，每小區(qū)10株，造林株行距2 m×2 m。自然整枝未進行間伐，采脂當年試驗林保留密度在1 125～1 350株·hm-2之間。

1.2 調(diào)查內(nèi)容與方法

2011年7月采脂前，測定試驗林樹高（H）、胸徑（D），單株材積（V）按照標準公式V = 0.714 265 437 ×104D1.867008H0.9014632計算材積[16]。2011年8－12月采用下降式采脂法[6]每天采割1次，每10 d稱1次松脂重量；分單株記錄產(chǎn)量并計算產(chǎn)脂力（M），產(chǎn)脂力為每天每10 cm割溝長松脂重量。采脂結(jié)束后測定割面寬度、割溝長度及割面處周長。割面寬與采割處周長均在第一刀“V”形底部測量。

1.3 統(tǒng)計分析方法

試驗中有關(guān)數(shù)據(jù)的方差統(tǒng)計分析均采用SAS 6.12[17]完成。表型變異系數(shù)（CV）和遺傳變異系數(shù)（GCV）、遺傳力、遺傳增益、遺傳相關(guān)分析參照相關(guān)文獻[18]進行計算，采用選擇指數(shù)方法[19-20]對家系進行綜合評價。

2 結(jié)果與分析

2.1 馬尾松家系半同胞子代生長與產(chǎn)脂性狀的遺傳變異

通過對馬尾松半同胞家系的生長及產(chǎn)脂性狀的表型變異分析可知（表 1），生長性狀與產(chǎn)脂性狀在半同胞子代家系間都存在豐富的變異。生長性狀在表型變異及遺傳變異上都稍高于產(chǎn)脂性狀，其中材積的遺傳變異系數(shù)最大，達到12.41%，產(chǎn)脂力的遺傳變異系數(shù)最小為4.74%?？梢婑R尾松半同胞子代家系的表型選擇有一定效果，通過選擇有望獲得遺傳品質(zhì)優(yōu)良的家系。

表1 馬尾松半同胞家系生長與產(chǎn)脂性狀表型變異Table 1 Phenotypic variation of growth and resin yield in half-sib families of P. massoniana

進一步方差分析結(jié)果表明（表2），20年生馬尾松自由授粉半同胞子代家系間的胸徑、材積、產(chǎn)脂力間均存在極顯著的差異，樹高性狀存在顯著差異。生長性狀遺傳方差分量高于產(chǎn)脂性狀，表明樹高、胸徑、材積生長受到遺傳控制的程度高于產(chǎn)脂性狀。參試家系各性狀家系遺傳力都較高，樹高、胸徑、材積的家系遺傳力都在0.52以上，產(chǎn)脂性狀的家系遺傳力也達到了0.468。這表明生長以及產(chǎn)脂性狀仍受中等強度以上的遺傳控制，較大的變異能夠提供豐富的選擇空間，為進一步選擇提供了可能。

表2 馬尾松半同胞家系生長與產(chǎn)脂性狀方差分析和遺傳估算Table 2 ANOVA and genetic estimation of growth and resin yield in half-sib families of P. massoniana

2.2 馬尾松家系半同胞子代生長與產(chǎn)脂性狀的相關(guān)分析

馬尾松子代半同胞家系各性狀的相關(guān)分析見表 3。樹高、胸徑、材積生長性狀之間及生長與產(chǎn)脂性狀之間均達顯著的表型正相關(guān)和遺傳正相關(guān)。生長性狀之間的相關(guān)性高于生長性狀與產(chǎn)脂性狀間的相關(guān)性，生長性狀中胸徑與產(chǎn)脂性狀間的遺傳正相關(guān)最緊密，達到0.972 8，表現(xiàn)出極顯著的正相關(guān)；樹高與產(chǎn)脂間的遺傳正相關(guān)最弱，為0.766 9；而材積與產(chǎn)脂力也表現(xiàn)出極顯著的遺傳正相關(guān)，因此，在改良生長性狀的同時，產(chǎn)脂性狀也得到一定的間接增益，可對兩者同時進行改良。

表3 馬尾松半同胞家系各性狀表型相關(guān)系數(shù)及遺傳相關(guān)系數(shù)Table 3 Phenotypic and genetic correlation coefficient of the traits in half-sib families of P. massoniana

2.3 生長量和產(chǎn)脂性狀評價及優(yōu)良家系選擇

參試家系的平均材積和產(chǎn)脂力均超過對照種源對應(yīng)的平均值（表4）。以材積為評價指標時，48個半同胞家系中材積超過對照平均值10%以上的家系占29.2%，超過對照平均值20%的占18.8%，超過對照平均值30%的占6.3%。以產(chǎn)脂力為評價指標時，產(chǎn)脂力超過對照平均值10%以上的家系占25%，超過對照平均值20%的占12.5%，超過對照平均值30%的占4.2%。表4同時表明，雖然以材積和產(chǎn)脂力作為評價指標是基于兩個不同的育種目的，2個單性狀分別進行選擇時家系前10%～30%的排名也存在一些差異，但是從兩者分別與對照的比較表明在選擇過程中材積與產(chǎn)脂力可以得到一定程度的兼顧。

根據(jù)生長以及產(chǎn)脂性狀的遺傳差異和相關(guān)分析結(jié)果，以20年生時馬尾松半同胞家系的H，D，V，M為變量，以家系均值標準差的倒數(shù)為權(quán)重配合建立指數(shù)選擇方程，得到的方程式為：

Ii= 0.522 8H + 0.300 7D + 9.586 1V + 0.437 7M。以高出材率、高產(chǎn)脂量為選擇目標，進行優(yōu)良家系選擇。

表4 不同選育目標下馬尾松半同胞家系材積和產(chǎn)脂力與對照均值比較Table 4 Comparison on volume and resin yield in half-sib families of P. massoniana with that of the control under the different selection target

根據(jù)群體均值以及各家系性狀均值代入指數(shù)方程計算得到群體及各家系聚合指數(shù)分別為 I0及II，其值分別為：I0=11.637 1，δI= 2.044 5。令各家系均值標準差為δI，便可得到區(qū)間，小于的家系綜合性狀差于，之間的家系可以作為綜合選擇對象，大于的家系為綜合表現(xiàn)優(yōu)良家系。

應(yīng)用選擇指數(shù)對半同胞家系進行綜合評價，按照前20%排名選擇出56號，54號，33號，6號等10個生長與產(chǎn)脂性狀優(yōu)良的家系。表5顯示，入選家系的各個性狀具有較高的遺傳改良潛力，樹高、胸徑、材積、產(chǎn)脂力的平均遺傳增益分別為16.99%，17.02%，18.39%，23.92%，分別比群體平均值大5.23%，14.66%，0.965%，15.57%。另外，按照產(chǎn)脂力遺傳增益的大小排序，發(fā)現(xiàn)對于家系的材積增益排序存在一定程度的差異，產(chǎn)脂力高的家系材積增益僅略高于群體均值。雖然綜合指數(shù)選擇存在一定的缺陷，但通過多性狀的綜合指數(shù)選擇，排名靠前的家系材積、產(chǎn)脂力得到的遺傳增益都獲得較高的遺傳增益，因此指數(shù)選擇在實際應(yīng)用中具有較強的可行性。

表5 優(yōu)良家系各性狀的遺傳增益及選擇指數(shù)Table 5 Genetic gain and selection index of traits in selected families

2.4 選擇方法及與對照的比較分析

在林木育種過程中，根據(jù)不同的選育目的采用針對性的方法會得到不同的結(jié)果。選擇方法有多種，或分析單一性狀選擇或多性狀綜合選擇，在分別按照單性狀方法選擇條件下取排名前30%家系及與對照的遺傳增益分析做作了比較，結(jié)果見表6。

表6 家系的單性狀選擇及與對照比較分析Table 6 Comparison and analysis of single sex selection in family and comparison with control

結(jié)合表5和表6分析，按照材積和產(chǎn)脂力2個單性狀選擇時排名前30%的家系在指數(shù)選擇方法結(jié)果中均出現(xiàn)，比較產(chǎn)脂力的遺傳增益及材積指標時，指數(shù)選擇方法與單一性狀選擇家系的相對順序排名差別不大，但是對同一個家系分析時，產(chǎn)脂力遺傳增益及材積的順序卻有較大差別，材積指標上，51號家系在兩種方法下都排名第一，但產(chǎn)脂力的遺傳增益均不高，排名靠后，56號，54號家系在用兩種方法比較時，產(chǎn)脂力遺傳增益及材積均排名前4名，但是在使用指數(shù)選擇方法下這2個家系的產(chǎn)脂力遺傳增益分別排在第1，第2名，單性狀選擇時排名則為第1，第3名。兩種方法選擇時材積的家系排名變化大于產(chǎn)脂力的遺傳增益家系排名。由此可見，利用綜合指數(shù)和單性狀選擇出的優(yōu)良家系的材積、產(chǎn)脂力增益是存在差別的，材積優(yōu)良家系排序、產(chǎn)脂力優(yōu)良家系排序也各不一致。綜合2種方法，指數(shù)選擇法下產(chǎn)脂力、材積的遺傳增益家系排名都排在前列，單性狀選擇法下前30%的入選家系不盡相同，且2種方法選擇出的相同家系排名也相對靠后。2種方法得出的家系排名與對照相比增益也不同，家系排名前30%的材積與產(chǎn)脂力遺傳增益均高于對照種源均值?？傮w上用單性狀選擇時材積的遺傳增益平均提高了6.22%，而產(chǎn)脂力則下降了3.46%。在家系間進行比較分析發(fā)現(xiàn)，用單一性狀選擇時有2個家系的材積遺傳增益下降，降幅在5.89%～11.72%，其余家系的材積遺傳增益均提高。產(chǎn)脂性狀卻表現(xiàn)不同，4個家系的產(chǎn)脂力遺傳增益下降，降幅在0.35%～14.28%。這表明采用2種不同方法對材積和產(chǎn)脂力選擇有明顯的影響，為獲得更高的遺傳增益，在強調(diào)材積性狀時宜采用單一性狀方法，在強調(diào)材積時兼顧產(chǎn)脂力性狀時宜采用指數(shù)選擇方法。另外，從表6中還可以看出，對照的平均產(chǎn)脂力低于家系產(chǎn)脂力，僅為家系產(chǎn)脂力的79.75%，其中 CK3產(chǎn)脂力最高，CK4最低。對照的平均遺傳增益低于家系平均遺傳增益，僅為家系遺傳增益的79.21%。通過指數(shù)選擇法對產(chǎn)脂性狀進行評價后發(fā)現(xiàn)，其遺傳增益均值均高于對照，這表明指數(shù)選擇結(jié)果更好些。

3 結(jié)論與討論

53參試個馬尾松半同胞子代家系的樹高、胸徑、材積、產(chǎn)脂力在家系間存在極顯著差異（P＜0.01）。生長性狀的家系遺傳力都在0.52以上，其中以胸徑家系遺傳力最高，產(chǎn)脂性狀家系遺傳力最低，但也達到0.46，表明生長性狀和產(chǎn)脂性狀都受到中等強度以上的遺傳控制，在進行家系選擇時具有較高的遺傳穩(wěn)定性。這樣的家系遺傳力與范林元等[21]、周志春[22]、秦國峰[23]、劉月蓉[1]的研究結(jié)果基本一致。在家系水平上，生長性狀與產(chǎn)脂性狀都表現(xiàn)出顯著的正遺傳相關(guān)，其中，胸徑與產(chǎn)脂具有極顯著的遺傳正相關(guān)。

若只考慮材積性狀或產(chǎn)脂力性狀時選擇出一批高材積與高產(chǎn)脂力優(yōu)樹，分別與對照的材積和產(chǎn)脂力進行比較，雖然單一性狀選擇方法與綜合指數(shù)選擇法各自強調(diào)的重點不同，從入選的結(jié)果看也有一定的差別，但是 2種方法下所選擇出的排名靠前的優(yōu)良家系在綜合指數(shù)選擇方法下相同的優(yōu)良家系排名更靠前，用綜合指數(shù)選擇法得出的是材積、產(chǎn)脂力遺傳增益都較高的優(yōu)良家系，它不一定是某個性狀特別突出的家系，但應(yīng)該是較高增益的各個性狀的平衡體。這將為后期高世代育種提供更豐富的遺傳材料。

用樹高、胸徑、材積和產(chǎn)脂力作多性狀指數(shù)選擇方程對家系進行生長性狀和產(chǎn)脂性狀的聯(lián)合評價，選擇出10個優(yōu)良家系，入選家系樹高、胸徑、材積、產(chǎn)脂力的平均遺傳增益分別為16.99%，17.02%，18.39%，23.92%，分別比群體平均值大增加5.23%，14.66%，0.965%，15.57%。通過多性狀指數(shù)選擇，彌補了單一性狀選擇的缺陷，說明是可行的選擇方法。

開展馬尾松高產(chǎn)脂家系選擇對于提高生產(chǎn)力十分重要，本試驗表明馬尾松半同胞家系的生長與產(chǎn)脂性狀選擇取得較好的效果。

[1] 劉月蓉. 高產(chǎn)脂馬尾松半同胞、全同胞子代產(chǎn)脂量測定[J]. 西南林學院學報，2005，25（3）：33－35.

[2] 呂鵬信，潘志剛. 濕地松高產(chǎn)脂家系的產(chǎn)脂量變異及增益[J]. 林業(yè)科學研究，1992，5（6）：700－705.

[3] 徐六一，虞沐奎. 濕地松高產(chǎn)脂品系早期性狀的研究[J]. 安徽農(nóng)業(yè)科學，2001，29（2）：228－230.

[4] 胡興宜，許業(yè)洲，楊代貴，等. 高產(chǎn)脂濕地松引種評價及優(yōu)樹選擇[J]. 林業(yè)科技開發(fā)，2006，20（1）：19－21.

[5] 莊偉瑛，張玉英，鄒元熹. 高產(chǎn)脂濕地松選擇和相關(guān)因子的分析[J]. 江西農(nóng)業(yè)大學學報：自然科學版，2008，29（1）：55－65.

[6] 冉澤文. 影響松脂采集增產(chǎn)的因素[J]. 林產(chǎn)化工通訊，2001，35（2）：35－36.

[7] 李思廣，付玉嬪，張快富，等. 高產(chǎn)脂思茅松半同胞子代測定[J]. 浙江林學院學報，2008，25（2）：158－162.

[8] 葉培忠. 杉木遺傳×環(huán)境互作和遺傳穩(wěn)定性研究[J]. 南京林學院學報，1980，3（3）：35－46.

[9] 王明庥. 黑楊新無性系生長的適應(yīng)性和遺傳穩(wěn)定性分析[J]. 南京林學院學報，1987（4）：13－25.

[10] 顧萬春. 毛白楊優(yōu)良無性系選育——生產(chǎn)力、遺傳穩(wěn)定性和適應(yīng)性評價[J]. 林業(yè)科學研究，1990，3（3）：222－227.

[11] 陳天華. 馬尾松種源試驗點效應(yīng)與穩(wěn)定性分析[J]. 南京林學院學報，1994，3（3）：19－25.

[12] 曾令海，王以珊，何波祥，等. 馬尾松高產(chǎn)脂力類型的選擇及子代測定的遺傳評價[J]. 廣東林業(yè)科技，1994（2）：35－38.

[13] MATZIRIS D. Variation in Growth, Flowering and cone production in a clonal seed orchard of Aleppo Pine Grown in Greece. [J]. Silvae Genet,1997，46，4：224－228.

[14] PAPAIOANNOC，I. The timits of Geographic Distribution of Pinus helepensis and Pinus brutia in NE Chalkidiki and Thein Associated Vegetation. [J]. Bull Nat Sc Athens，1995，1－13.

[15] 劉永紅，楊培華，樊軍鋒，等. 油松優(yōu)良家系多性狀選擇方法研究[J]. 西北農(nóng)林科技大學學報：自然科學版，2006，34（12）：115－120.

[16] 楊章旗. 馬尾松材性與產(chǎn)脂性狀遺傳改良研究[D]. 北京：北京林業(yè)大學，2012，55－67.

[17] 惠大豐，姜長鑒. 統(tǒng)計分析系統(tǒng)SAS軟件實用教程[M]. 北京：北京航空航天大學出版社，1996，46－54.

[18] 續(xù)九如. 林木數(shù)量遺傳學[M]. 北京：中國林業(yè)出版社，2006，61－78..

[19] 吳清金. 25年生馬尾松優(yōu)樹子代生產(chǎn)力評價及家系選擇[J]. 林業(yè)科技開發(fā)，2012，26（2）：73－76.

[20] 李光友，徐建民，陸釗華，等. 尾葉桉二代測定林家系的綜合評選[J]. 林業(yè)科學研究，2005，18（1）：57－61.

[21] 范林元，賴煥林，季孔庶，等. 馬尾松實生種子園建園家系遺傳值估算與優(yōu)良家系評選[J]. 東北林業(yè)大學報，2004，32（4）：3－11.

[22] 周志春. 馬尾松自由授粉家系生長和材質(zhì)的遺傳變異分析及聯(lián)合選擇[J]. 林業(yè)科學研究，1994，7（3）：263－268.

[23] 秦國峰. 馬尾松造紙材的變異、遺傳和選擇育種策略[J]. 林業(yè)科學研究，1995，8（5）：506－513.