日本落葉松遺傳變異分析及優(yōu)良家系無性系選擇

王芳 馬茂 王佳興 張含國

摘 要：為探索日本落葉松遺傳變異規(guī)律及選擇優(yōu)良家系無性系，以內(nèi)蒙古自治區(qū)赤峰市旺業(yè)甸實(shí)驗(yàn)林場的59個(gè)日本落葉松無性系測定林和81個(gè)日本落葉松家系子代試驗(yàn)林為研究對象，測定其樹高，對其進(jìn)行方差分析和變異分析。結(jié)果表明，日本落葉松無性系和家系樹高性狀存在豐富的變異，變異系數(shù)范圍分別為17.84%～20.67%和24.68%～31.19%。方差分析表明，每一年，無性系間和家系間均存在顯著差異，重復(fù)力和遺傳力范圍分別為0.482～0.494和0.793～0.819，屬于中等重復(fù)力和高遺傳力。以20%的入選率為標(biāo)準(zhǔn)，最終選擇出12個(gè)優(yōu)良無性系和16個(gè)優(yōu)良家系，無性系現(xiàn)實(shí)增益和家系現(xiàn)實(shí)增益范圍分別為5.73%～16.48%和7.59%～15.35%。

關(guān)鍵詞：日本落葉松；樹高；家系；無性系；變異分析

中圖分類號：S791.223??? 文獻(xiàn)標(biāo)識碼：A?? 文章編號：1006-8023（2023）04-0048-10

Analysis of Genetic Variation and Selection of Superior

Families and Clones of Larix kaempferi

WANG Fang1， MA Mao1， WANG Jiaxing2， ZHANG Hanguo2*

（1.Wangyedian Experimental Forestry Farm of Karqin Qi， Chifeng 024423， China;

2.State Key Laboratory of Tree Genetics and Breeding， Northeast Forestry University， Harbin 150040， China）

Abstract：In order to explore the genetic variation pattern of Larix kaempferi and to select superior family and clones， 59 clones and 81 families of Larix kaempferi in Wangyedian Experimental Forestry Farm， Chifeng City， Inner Mongolia Autonomous Region， were taken as research objects， and their tree heights were measured， and their ANOVA and variance analysis were conducted. The results showed that there was abundant variation of tree height traits in Larix kaempferi clones and families， with coefficients of variation ranging from 17.84%-20.67% and 24.68%-31.19%， respectively. ANOVA showed the difference was highly significant between clones and families every year， with a range in 0.482-0.494 and 0.793-0.819 of repeatability and heritability， respectively， which were moderate repeatability and high heritability. Using a 20% entry rate， 12 superior clones and 16 superior families were finally selected， with realistic gains of clones and realistic gains of families ranging were 5.73% to 16.48% and 7.59% to 15.35%， respectively.

Keywords：Larix kaempferi; tree height; family; clone; analysis of variance

收稿日期：2022-09-17

基金項(xiàng)目：國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目研究任務(wù)（2022YFD2201001）

第一作者簡介：王芳，高級工程師。研究方向?yàn)榱帜具z傳育種。E-mail： 651387751@qq.com

通信作者：張含國，博士，教授。研究方向?yàn)榱帜具z傳育種。E-mail： hanguozhang1@ sina.com

引文格式：王芳，馬茂，王佳興，等.日本落葉松遺傳變異分析及優(yōu)良家系無性系選擇[J].森林工程，2023，39（4）：48-57.

WANG F， MA M， WANG J X， et al.Analysis of genetic variation and selection of superior families and clones of Larix kaempferi[J]. Forest Engineering， 2023， 39（4）：48-57.

0 引言

日本落葉松（Larix kaempferi）分布在日本本州中部地區(qū)和關(guān)東地區(qū)[1]。日本落葉松樹干端直，姿態(tài)優(yōu)美，葉色翠綠，適應(yīng)范圍廣，生長初期較快，抗病性較強(qiáng)，木材力學(xué)性能較高，有較好的耐腐性，可做建筑材料和工業(yè)用材的原料，并可從其木材中提取松節(jié)油、酒精和纖維素等化學(xué)物品，用途很廣[2-3]。日本落葉松早在20世紀(jì)中期就引入我國，東北、華北、西北乃至華中、西南地區(qū)都有栽植，其中山東、遼寧引種最早，已形成大面積人工林[4]。

變異系數(shù)大小可反映群體的變異程度，是衡量性狀遺傳變異潛力的有效指標(biāo)[5-6]。羅芊芊等[7]在對馬尾松無性系進(jìn)行研究時(shí)發(fā)現(xiàn)，馬尾松無性系生長和形質(zhì)性狀存在豐富的遺傳差異，各性狀變異系數(shù)變幅為12.44%～36.36%；肖德卿等[8]在對紅豆樹優(yōu)樹家系進(jìn)行研究時(shí)發(fā)現(xiàn)，3年生紅豆樹生長和形質(zhì)性狀在家系間的遺傳差異均達(dá)到極顯著水平（P＜0.01），兩地點(diǎn)生長和形質(zhì)性狀的變異系數(shù)變幅為23.3%～65.1%。

增益是衡量選擇效果的最重要的參數(shù)之一，反映了下一代比親本可能增加的收獲量[9]。聶林芽等[10]對杉木半同胞家系進(jìn)行研究得出結(jié)論，優(yōu)良家系地徑遺傳增益為10.17%～16.46%，家系增產(chǎn)效果相對顯著；黃云鵬等[11]對42個(gè)4年生的福建山櫻花嫁接無性系進(jìn)行研究，最終選出速生福建山櫻花無性系4個(gè)，樹高、胸徑、冠幅遺傳增益均值分別為11.12%、13.35%、9.52%。

在林木良種選育過程中，優(yōu)良家系選擇至關(guān)重要，在選擇過程中既要考慮廣泛的遺傳基礎(chǔ)，又應(yīng)考慮選擇的遺傳穩(wěn)定性[12]。葛清鋒等[13]以17年生日本落葉松半同胞家系為研究對象，選擇出優(yōu)良家系并對其早期選擇結(jié)果進(jìn)行驗(yàn)證；潘艷艷等[14]以67個(gè)日本落葉松全同胞家系為材料進(jìn)行分析選擇出優(yōu)良家系；張正剛等[15]以36個(gè)日本落葉松自由授粉家系為研究對象進(jìn)行研究，最終選擇出8個(gè)優(yōu)良家系。

無性系與傳統(tǒng)的實(shí)生苗相比，具有目的性強(qiáng)、集約經(jīng)營和生產(chǎn)周期短的特點(diǎn)，優(yōu)良遺傳材料的取得和大規(guī)模的無性增殖利用是無性系林業(yè)的發(fā)展基礎(chǔ)之一[16]。杜彥昌等[17]以小隴山林區(qū)18個(gè)日本落葉松無性系為研究對象，選擇出4個(gè)優(yōu)良適生無性系；杜超群等[18]以11年生日本落葉松混系超級苗無性系為研究對象，開展優(yōu)良無性系初步選擇，最終選擇出5個(gè)優(yōu)良無性系。

為了獲得高產(chǎn)穩(wěn)產(chǎn)、遺傳增益高的良種，對無性系和家系進(jìn)行評價(jià)尤為重要[19]。本研究以內(nèi)蒙古自治區(qū)赤峰市旺業(yè)甸林場的59個(gè)日本落葉松無性系測定林和81個(gè)日本落葉松家系子代測定林為研究對象，測量其樹高，對其樹高性狀進(jìn)行變異分析和方差分析，選擇出優(yōu)良家系和無性系，為日后日本落葉松優(yōu)良家系無性系的大面積推廣以及營造速生林提供優(yōu)良材料。

1 材料與方法

1.1 試驗(yàn)地條件情況

旺業(yè)甸實(shí)驗(yàn)林場位于內(nèi)蒙古自治區(qū)赤峰市喀喇沁旗西南部，地處燕山山脈北麓七老圖山支脈，地理位置東經(jīng)118°09′～118°30′，北緯41°21′～41°39′，海拔800～1 890 m。經(jīng)營區(qū)屬暖溫帶半干旱地區(qū)，為明顯的大陸性季風(fēng)氣候。年降水量400～600 mm，年平均氣溫4.2 ℃，無霜期117 d，土壤以典型的棕色森林土為主。

1.2 試驗(yàn)材料

日本落葉松無性系測定林位于旺業(yè)甸實(shí)驗(yàn)林場美林營林區(qū)前敖包溝61林班45小班，試驗(yàn)材料來自日本落葉松種子園，在種子園分無性系采集穗條并嫁接育苗。2013年4月嫁接育苗，2014年4月造林，參試無性系共59個(gè)（表1），單行2株小區(qū)5次重復(fù)，株行距1.5 m×2 m，隨機(jī)排列，整地方式為穴狀。日本落葉松家系子代測定林位于旺業(yè)甸實(shí)驗(yàn)林場古山營林區(qū)南岔林班21小班，自由授粉家系子代測定種子來自日本落葉松種子園。2015年4月分家系進(jìn)行田間試驗(yàn)育苗，2017年4月造林，參試家系共81個(gè)（表1），單行10株小區(qū)6次重復(fù)，株行距1.5 m×2 m。

1.3 試驗(yàn)方法

2019—2021年秋季當(dāng)日本落葉松停止生長后，用塔尺對家系和無性系測定林當(dāng)年樹高進(jìn)行測量。數(shù)據(jù)分析采用Excel和Spss 22.0統(tǒng)計(jì)分析軟件進(jìn)行分析。

以家系小區(qū)平均值為單元進(jìn)行方差分析，方差分析采用線性模型

Yijk=μ+Bi+Fj+BFij+eijk

式中：Yijk為第i個(gè)區(qū)組第j個(gè)家系第k個(gè)觀測值；μ為群體平均效應(yīng)；Bi為第i個(gè)區(qū)組效應(yīng)；Fj為第j家系的效應(yīng)；BFij為第i家系和第j區(qū)組的互作效應(yīng)；eijk為機(jī)誤。

變異系數(shù)（CV）計(jì)算公式

CV=S/X×100%。

式中：S表示表型標(biāo)準(zhǔn)差；X表示樹高性狀平均值。

家系遺傳力計(jì)算公式

h2=（MSf-MSe）/MSf=1-1/F。

式中：MSf表示家系均方；MSe表示誤差均方；F表示方差分析中的F值。

家系遺傳增益、現(xiàn)實(shí)增益計(jì)算公式[20]

ΔG=（Xi-X）Δh2/X×100%。

G=（Xi-X）/X×100%。

式中：ΔG為遺傳增益；G為現(xiàn)實(shí)增益；Xi為入選家系平均值；X為家系平均值；h2為家系遺傳力。

無性系重復(fù)力計(jì)算公式

R=1-1/F。

無性系遺傳增益、現(xiàn)實(shí)增益計(jì)算公式

ΔG=RW/X×100%。

G=W/X×100%。

式中：W為選擇差；R為無性系重復(fù)力。

2 結(jié)果與分析

2.1 日本落葉松無性系樹高變異分析及初步選擇

2.1.1 日本落葉松無性系樹高變異分析

對7—9年生日本落葉松無性系進(jìn)行變異分析，結(jié)果見表2。日本落葉松樹高性狀存在豐富的變異，并且隨著樹齡的增長，變異系數(shù)逐漸降低。變異系數(shù)范圍在17.84%～20.67%，平均變異系數(shù)為19.05%。

對不同年度的日本落葉松各無性系進(jìn)行變異分析，并列出變異系數(shù)較大前5位和較小后5位的無性系。

7年生日本落葉松無性系樹高平均變異系數(shù)為20.67%，其中變異系數(shù)較高的5個(gè)無性系為81、D2、87、73081和38，平均變異系數(shù)為33.82%，比總體平均變異系數(shù)高13.15%；變異系數(shù)較低的5個(gè)無性系為73118、93、RT1、2019和S2，平均變異系數(shù)為9.29%，比總體平均變異系數(shù)低11.38%。

8年生日本落葉松無性系樹高平均變異系數(shù)為18.67%，其中變異系數(shù)較高的5個(gè)無性系為95、81、D2、38和73096，平均變異系數(shù)為31.27%，比總體平均變異系數(shù)高12.60%；變異系數(shù)較低的5個(gè)無性系為41、88、94、93和RT1，平均變異系數(shù)為8.55%，比總體平均變異系數(shù)低10.12%。

9年生日本落葉松無性系樹高平均變異系數(shù)為17.84%，其中變異系數(shù)較高的5個(gè)無性系為81、73113、73108、86和38，平均變異系數(shù)為27.83%，比總體平均變異系數(shù)高9.99%；變異系數(shù)較低的5個(gè)無性系為76、RB、S9、93和41，平均變異系數(shù)為9.24%，比總體平均變異系數(shù)低8.60%。

可見，日本落葉松無性系存在著豐富的遺傳變異，除9年生變異較小的5個(gè)無性系平均變異系數(shù)高于8年生外，其他變異系數(shù)均隨著樹齡的變化而降低。3個(gè)年度共有的變異大的家系有81和38，變異小的家系有93。

2.1.2 日本無性系樹高方差分析及無性系初步選擇

對日本落葉松無性系各年度樹高進(jìn)行方差分析，結(jié)果見表3。3個(gè)年度在無性系間均表現(xiàn)出差異極顯著，3個(gè)年度重復(fù)力分別為0.494、0.482和0.489，屬于中等重復(fù)力水平，可以進(jìn)行無性系選擇。

入選率為10%時(shí)，遺傳增益分別為7.02%、6.4%和6.72%；入選率為20%時(shí)，遺傳增益分別為5.92%、5.04%和5%。通過估算遺傳增益可以看出，由于總體變異系數(shù)逐年減小，遺傳增益隨著林齡增長而略有降低，但林齡越大，選擇結(jié)果會越準(zhǔn)確。

通過樹高均值和無性系變異系數(shù)，以20%入選率為標(biāo)準(zhǔn)，對9年生日本落葉松無性系進(jìn)行選擇，最終選擇73109、73121、D81、2019、D61、76、72、73099、82、S4、79和44號12個(gè)優(yōu)良無性系，結(jié)果見表4?，F(xiàn)實(shí)增益范圍為5.73%～16.48%。

2.2 日本落葉松家系變異分析及初步選擇

2.2.1 日本落葉松家系樹高變異分析

對5—7年生的日本落葉松家系進(jìn)行變異分析，結(jié)果見表5。與無性系相似，家系樹高性狀也存在豐富的變異。并且隨著樹齡的增長，變異系數(shù)也逐漸降低。變異系數(shù)范圍為24.68%～31.19%，平均變異系數(shù)為28.16%。

對不同年度日本落葉松各家系進(jìn)行變異分析，并列出變異系數(shù)較高前5和較低后5位的家系。

5年生日本落葉松家系樹高平均變異系數(shù)為31.19%，其中變異系數(shù)較高的5個(gè)家系為99、44、87、88和S2，平均變異系數(shù)為37.63%，比總體平均變異系數(shù)高6.44%；變異系數(shù)較低的5個(gè)家系分別為73098、73108、76、D81和41，平均變異系數(shù)為21.44%，比總體平均變異系數(shù)低9.75%。

6年生日本落葉松家系樹高平均變異系數(shù)為28.61%，其中變異系數(shù)較高的5個(gè)家系為44、87、99、88和73094，平均變異系數(shù)為34.97%，比總體平均變異系數(shù)高6.36%；變異系數(shù)較低的5個(gè)家系分別為73098、76、73108、D81和41，平均變異系數(shù)為18.08%，比總體平均變異系數(shù)低10.53%。

7年生日本落葉松家系樹高平均變異系數(shù)為24.68%，其中變異系數(shù)較高的5個(gè)家系為44、87、99、73101和D41，平均變異系數(shù)為30.82%，比總體平均變異系數(shù)高6.14%；變異系數(shù)較低的5個(gè)家系分別為RT2、D81、82、S4和41，平均變異系數(shù)為15.12%，比總體平均變異系數(shù)低9.56%。

可見，日本落葉松家系存在著豐富的遺傳變異，并且隨著樹齡的增長，變異系數(shù)逐年降低。3個(gè)年度共有的變異大的家系有44、87和99；共有的變異小的家系有41和D81。

2.2.2 日本落葉松家系樹高方差分析及家系初步選擇

對日本落葉松家系樹高進(jìn)行方差分析，結(jié)果見表6。3個(gè)年份在家系以及區(qū)組家系互作間均差異極顯著。3個(gè)年度年樹高性狀遺傳力分別為0.819、0.793和0.804，受到較強(qiáng)遺傳控制，適合在家系水平進(jìn)行遺傳改良。

入選率為10%時(shí)，遺傳增益分別為13.75%、11.09%和9.66%；入選率為20%時(shí)，遺傳增益分別為11.74%、9.43%和8.51%。與無性系相似，總體遺傳增益系數(shù)逐年減小。林齡越大，選擇結(jié)果越準(zhǔn)確。

通過樹高均值和家系變異系數(shù)，以20%入選率為標(biāo)準(zhǔn)，對7年生日本落葉松家系進(jìn)行選擇，最終選擇82、77、D82、73098、D81、D61、76、73105、RT2、74、73109、S9、81、D12、94和73108號16個(gè)優(yōu)良家系，結(jié)果見表7。現(xiàn)實(shí)增益范圍為7.59%～15.35%。

2.3 無性系與家系對比分析

2.3.1 變異對比分析

對57個(gè)具有相同系號的日本落葉松家系和無性系變異進(jìn)行對比分析，結(jié)果見表8。由表8發(fā)現(xiàn)同一年齡下家系平均變異系數(shù)要大于無性系平均變異系數(shù)，符合有性繁殖變異較大的理論，7年生家系變異系數(shù)比無性系變異系數(shù)高2.72%。

對7年生各家系和無性系進(jìn)行變異分析，結(jié)果發(fā)現(xiàn)排名前5的無性系變異系數(shù)大于排名前5的家系變異系數(shù)，而排名后5的無性系變異系數(shù)要小于排名后5的家系變異系數(shù)，說明家系、無性系變異較大或較小系號是不同的。排名前5的無性系平均變異系數(shù)比家系高3.22%，排名后5的無性系平均變異系數(shù)比家系低5.68%。

2.3.2入選系號及增益對比分析

對7年生日本落葉松相同家系和無性系進(jìn)行方差分析，結(jié)果見表9。樹高性狀在家系和無性系間均差異極顯著。同一年度家系遺傳力大于無性系重復(fù)力。

通過樹高均值和變異系數(shù)，以20%入選率為標(biāo)準(zhǔn)，對日本落葉松無性系和家系進(jìn)行選擇，最終選擇出73121、2019、73109、D81、D61、S2、D24、72、82、S4和76號11個(gè)優(yōu)良無性系和82、77、73098、D81、D61、76、74、73109、S9、81和94號11個(gè)優(yōu)良家系，結(jié)果見表10和表11。無性系現(xiàn)實(shí)增益為7.62%～16.09%；家系現(xiàn)實(shí)增益范圍為6.30%～13.60%。

從結(jié)果上看，選出的優(yōu)良家系和優(yōu)良無性系相同系號有82、73109、D81、D61號4個(gè)系號，說明這4個(gè)系號有性繁殖和無性繁殖均能獲得較大增益。

3 結(jié)論和討論

林木的遺傳變異是遺傳改良的前提，是遺傳信息的重要表征。豐富的遺傳變異決定了物種在育種過程中的改良潛力，掌握林木的遺傳變異規(guī)律，是育種策略制定的基礎(chǔ)[21-23]。本研究結(jié)果顯示，日本落葉松無性系和家系樹高性狀存在豐富的變異，7—9年生無性系變異系數(shù)范圍為17.84%～20.67%，平均變異系數(shù)為19.05%；5—7年生家系變異系數(shù)范圍為24.68%～31.19%，平均變異系數(shù)為28.16%。研究發(fā)現(xiàn)，無論是無性系還是家系，樹齡越大，變異系數(shù)越小，這與馬茂等[24]對日本落葉松無性系和繆小飛等[25]對雜種落葉松家系的研究結(jié)果一致。同一年齡相同系號的日本落葉松家系和無性系變異對比分析，結(jié)果顯示家系平均變異系數(shù)要大于無性系平均變異系數(shù)，而且家系、無性系變異較大或較小系號是不同的。

在林木育種研究中，方差分析是評估變異幅度的重要方法[26]。方差分析結(jié)果顯示，3個(gè)年份樹高在無性系間和家系間差異極顯著。遺傳力和重復(fù)力是衡量某一性狀在不同時(shí)間、不同地點(diǎn)其表型特征是否可以持續(xù)穩(wěn)定遺傳的判定指標(biāo)，遺傳力（重復(fù)力）越高，說明受外界環(huán)境影響越小[27，28]。3個(gè)年度樹高性狀無性系重復(fù)力范圍為0.482～0.494，屬于中等重復(fù)力，表現(xiàn)較為穩(wěn)定。3個(gè)年度家系遺傳力范圍為0.793～0.819，屬于高遺傳力，受到較強(qiáng)遺傳控制。相同材料同一年度家系遺傳力大于無性系重復(fù)力。

日本落葉松無性系入選率為20%時(shí)，3個(gè)年度樹高性狀遺傳增益分別為5.92%、5.04%和5%。家系入選率為20%時(shí)，3個(gè)年度遺傳增益分別為11.74%、9.43%和8.51%。由于試驗(yàn)材料的遺傳差異隨著年齡增長會發(fā)生變化，因此選擇樹齡較大的材料，選擇的結(jié)果會相對準(zhǔn)確一些。

通過樹高平均值和變異系數(shù)，以20%入選率為標(biāo)準(zhǔn)，最終選擇出12個(gè)優(yōu)良無性系和16個(gè)優(yōu)良家系，其中D81、D61、76、73109、76、82號為無性系和家系共有系號，表明這6個(gè)系號有性繁殖和無性繁殖均能獲得較大增益。通過對比共有家系和無性系，最終選擇出包括S2、D24、71號在內(nèi)的15個(gè)無性系和16個(gè)家系。

【參 考 文 獻(xiàn)】

[1]中國科學(xué)院中國植物志編輯委員會.中國植物志[M].北京：科學(xué)出版社，1978.

Editorial Committee of China Plant Journal， Chinese Academy of Sciences. China plant journal[M]. Beijing： Science Press， 1978.

[2]LAI M， SUN X M， CHEN D S， et al. Age-related trends in genetic parameters for Larix kaempferi and their implications for early selection[J]. BMC Genetics， 2014， 15（1）：1-8.

[3]FUKATSU E， TSUBOMURA M， FUJISAWA Y， et al. Genetic improvement of wood density and radial growth in Larix kaempferi： result from a diallel mating test[J]. Annals of Forest Science， 2013， 70：451-459.

[4]滕秀蘭，王福森，李樹森，等.日本落葉松品種引進(jìn)及育苗造林技術(shù)探析[J].防護(hù)林科技，2018（1）：91-92.

TENG X L， WANG F S， LI S S， et.al. Introduction of Larix kaempferi species and seedling planting technology[J]. Protection Forest Science and Technology， 2018（1）：91-92.

[5]LEPOITTEVIN C， ROUSSEAU J P， GUILLEMIN A， et al. Genetic parameters of growth， straightness sand wood chemistry traits in Pinus pinaster[J]. Annals of Forest Science， 2011， 68（4）：873-884.

[6]高玉紅，魏萌，李凱迪，等.不同紅河橙居群果實(shí)表型的變異[J].經(jīng)濟(jì)林研究，2021，39（4）：88-96.

GAO Y H， WEI M， LI K D， et al. Phenotypic variation of Citrus hongheensis fruit from different provenances[J]. Non-wood Forest Research， 2021， 39（4）： 88-96.

[7]羅芊芊，肖德卿，徐洪輝，等.馬尾松無性系生長和形質(zhì)性狀變異分析和優(yōu)選[J].東北林業(yè)大學(xué)學(xué)報(bào)，2022，50（9）：40-44.

LUO Q Q， XIAO D Q， XU H H， et al.Variation analysis and optimization of growth and form-quality traits of Pinus massoniana clones[J]. Journal of Northeast Forestry University， 2022， 50（9）：40-44.

[8]肖德卿，鄧章文，羅芊芊，等.幼齡紅豆樹生長和形質(zhì)性狀家系變異分析[J].林業(yè)科學(xué)研究，2021，34（3）：152-157.

XIAO D Q， DENG Z W， LUO Q Q， et al.Growth and form-quality traits diversity analysis of young Ormosia hosiei family[J]. Forest Research， 2021， 34（3）：152-157.

[9]RUOTSALAINEN S， LINDGREN D. Predicting genetic gain of backward and forward selection in forest tree breeding[J]. Silvae Genetica，1998， 47（1）：42-50.

[10]聶林芽，曾素平，厲月橋，等.杉木半同胞家系生長性狀遺傳分析及優(yōu)良家系選擇[J].林業(yè)科技通訊，2022（10）：3-9.

NIE L Y， ZENG S P， LI Y Q， et al. Genetic analysis of growth traits and selection of excellent families in half sib families of Cunninghamia lanceolata[J]. Forest Science and Technology， 2022（10）：3-9.

[11]黃云鵬，丁珌，潘惠忠，等.福建山櫻花無性系生長性狀遺傳變異分析及選擇[J].福建林業(yè)科技，2022，49（2）：7-10，25.

HUANG Y P， DING M， PAN H Z， et al. Genetic variation of growth trait and superior clone selection of Cerasus campanulata[J]. Joural of Fujian Forestry Science and Technology， 2022， 49（2）：7-10， 25.

[12]潘艷艷.日本落葉松種子園親本及其子代變異研究[D].哈爾濱：東北林業(yè)大學(xué)，2019.

PAN Y Y. Study on genetic variation of parents and progenies in Larix kaempferi seed orchard[D]. Harbin： Northeast Forestry University， 2019.

[13]葛清鋒，曾美，侯義梅.日本落葉松半同胞家系生長變異研究[J].湖北林業(yè)科技，2022，51（3）：6-9，21.

GE Q F， ZENG M， HOU Y M. Study on growth variation of Larix kaempferi half-sib families[J]. Hubei Forestry Science and Technology， 2022， 51（3）：6-9， 21.

[14]潘艷艷，許貴友，董利虎，等.日本落葉松全同胞家系苗期生長性狀遺傳變異[J].南京林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），2019，43（2）：14-22.

PAN Y Y， XU G Y， DONG L H， et.al. Genetic variations of seedling growth traits among full-sib families of Larix kaempferi[J]. Journal of Nanjing Forestry University （Natural Sciences）， 2019， 43（2）：14-22.

[15]張正剛，馬建偉，靳新春，等.日本落葉松 [J].西北林學(xué)院學(xué)報(bào)，2013，28（4）：74-79.

ZHANG Z G， MA J W， JIN X C， et al. Analysis and selection of progeny test forest of the open pollinated family of Larix kaempferi[J]. Journal of Northwest Forestry University， 2013， 28（4）：74-79.

[16]DONALDSON L， BURDON R. Clonal variation and repeatability of micro fibrilangle in Pinus radiata[J]. New Zealand Journal of Forestry Science， 1995， 25（2）：164-174.

[17]杜彥昌，馬建偉，安三平，等.小隴山林區(qū)日本落葉松無性系遺傳變異及初選[J].林業(yè)科技通訊，2018（7）：11-15.

DU Y C， MA J W， AN S P， et al. Genetic variation and primary selection of clones of Larix kaempferi in Xiaolong mountain forest area[J]. Forest Science and Technology， 2018（7）：11-15.

[18]杜超群，許業(yè)洲，孫曉梅，等.鄂西亞高山區(qū)日本落葉松無性系生長性狀變異分析與早期選擇[J].華中農(nóng)業(yè)大學(xué)學(xué)報(bào)，2015，34（3）：19-23.

DU C Q， XU Y Z， SUN X M， et al. Variation of growth traits and early selection of Larix kaempferi clones in sub-alpine area of western Hubei Province[J]. Journal of Huazhong Agricultural University， 2015， 34（3）：19-23.

[19]李巖，朱嘉瑤，王喜和，等.紅松優(yōu)樹無性系及其子代的生長評價(jià)與選擇研究[J].北京林業(yè)大學(xué)學(xué)報(bào)，2021，43（10）：38-46.

LI Y， ZHU J Y， WANG X H， et al. Growth evaluation and selection study of elite clones and its offspring families in Pinus koraiensis[J]. Journal of Beijing Forestry University， 2021， 43（10）：38-46.

[20]歐陽天林，朱柯帆，邱建勛，等.木荷種子園自由授粉家系生長遺傳變異及初選[J/OL].中南林業(yè)科技大學(xué)學(xué)報(bào)，2022（9）：17-23.

OUYANG T L， ZHU K F， QIU J X， et al. Genetic variation and selection of the free-pollinated families in Schima superba seed orchard[J]. Journal of Central South University of Forestry & Technology， 2022（9）：17-23.

[21]王云鵬，張蕊，周志春，等.木荷優(yōu)樹自由授粉家系早期生長性狀遺傳變異動態(tài)規(guī)律[J].林業(yè)科學(xué)，2020，56（9）：77-86.

WANG Y P， ZHANG R， ZHOU Z C， et al. Dynamic patterns of genetic variation in early growth traits of the open-pollinated families of Schima superba plus tree[J]. Scientia Silvae Sinicae， 2020， 56（9）：77-86.

[22]王瑋槐，陳雨桐，黨乾順，等.紅松無性系和優(yōu)樹種子性狀及營養(yǎng)成分分析[J].森林工程，2019，35（2）：11-20，26.

WANG W H， CHEN Y T， DANG Q S， et al.Analysis of characters and nutritional components of Pinus koraiensis seeds on clones and superior trees[J]. Forest Engineering， 2019， 35（2）：11-20， 26.

[23]周燕，黃盛怡，楊孟晴，等.赤皮青岡不同家系苗期生長和葉片性狀遺傳變異[J].西部林業(yè)科學(xué)，2022，51（5）：75-80.

ZHOU Y， HUANG S Y， YANG M Q， et al. Genetic variation of growth and leaf traits of different Cyclobalanopsis gilva families at seedling stage[J]. Journal of West China Forestry Science， 2022， 51（5）： 75-80.

[24]馬茂，徐開源，張含國.日本落葉松無性系樹高變異與初步選擇[J].林業(yè)科技，2019，44（5）：5-7.

MA M， XU K Y， ZAHNG H G. Tree height variation and preliminary selection of clones of Larix kaempferi[J]. Forestry Science & Technology， 2019， 44（5）：5-7.

[25]繆小飛，張含國，侯丹，等.雜種落葉松家系遺傳變異及多點(diǎn)穩(wěn)定性[J].東北林業(yè)大學(xué)學(xué)報(bào)，2018，46（12）：1-8.

MIU X F， ZHANG H G， HOU D， et al. Genetic variation and multi-point stability of families of Hybrid larch[J]. Journal of Northeast Forestry University， 2018， 46（12）：1-8.

[26]SAFAVI S A， POURDAD S A， MOHMMAD T， et al. Assessment of genetic variation among safflower （Carthamus tinctorius L.） accessions using agro-morphological traits and molecular markers[J]. Journal of Food Agriculture ＆ Environment， 2010， 8（3）：616-625.

[27]劉明國，張欣，董勝君，等.西伯利亞杏優(yōu)選無性系間數(shù)量性狀的差異與重復(fù)力研究[J].沈陽農(nóng)業(yè)大學(xué)學(xué)報(bào)，2015，46（5）：548-554.

LIU M G， ZHANG X， DONG S J， et al. Differences and repeatabilities of quantitative characters of Prunus sibirica superior clones[J]. Journal of Shenyang Agricultural University， 2015， 46（5）：548-554.

[28]史加亮，李鳳瑞，趙文超，等.陸地棉特異種質(zhì)主要性狀配合力和遺傳力分析[J].山東農(nóng)業(yè)科學(xué)，2021，53（1）：14-19.

SHI J L， LI F R， ZHAO W C， et al. Analysis on combining ability and heritability of main traits for special heteroplasm in upland cotton[J]. Shandong Agricultural Sciences， 2021， 53（1）：14-19.