基于球果和種子杉木品質(zhì)遺傳多樣性分析及綜合評(píng)價(jià)

程琳 戴俊 賀錦鋒 謝樂(lè) 梁機(jī) 董利軍 藍(lán)肖 黃開(kāi)勇 白天道 黃鴻飛

摘要：【目的】明確不同種質(zhì)杉木球果和種子品質(zhì)性狀的遺傳多樣性，并篩選優(yōu)質(zhì)杉木資源，為杉木的良種選育及種質(zhì)資源的高效利用提供科學(xué)依據(jù)。【方法】以廣西2代杉木種子園40個(gè)無(wú)性系的球果和種子為研究對(duì)象，采用方差分析、聚類分析、主成分分析等分析方法，分別從球果表型性狀、種子發(fā)芽特性和種子優(yōu)良度性狀對(duì)15個(gè)品質(zhì)指標(biāo)進(jìn)行遺傳多樣性分析和綜合評(píng)價(jià)?！窘Y(jié)果】40個(gè)杉木無(wú)性系球果和種子15個(gè)品質(zhì)性狀的遺傳多樣性指數(shù)介于1.4792～2.0803，均值為1.9090，其中遺傳多樣性指數(shù)最大的是空粒率，最小的是霉變率;變異系數(shù)介于4.55%～120.72%，均值為29.60%，其中變異系數(shù)最大的是霉變率，最小的是果形指數(shù);遺傳力介于0.366～0.969，均值為0.843，其中遺傳力最大的是球果橫徑，最小的是空粒率。除空粒率在無(wú)性系間差異顯著（P<0.05）外，其他14個(gè)品質(zhì)性狀在無(wú)性系間差異均達(dá)極顯著（P<0.01）。40個(gè)杉木無(wú)性系可劃分為四大類群，其中，類群Ⅱ的綜合性狀最優(yōu)，類群Ⅰ的種子優(yōu)良性狀最好。對(duì)15個(gè)杉木球果和種子品質(zhì)性狀進(jìn)行主成分分析，綜合性狀排名前8位的無(wú)性系編號(hào)分別是14、34、28、31、38、9、29、27?！窘Y(jié)論】廣西2代杉木種子園40個(gè)無(wú)性系球果和種子品質(zhì)性狀的遺傳多樣性豐富，篩選出的8個(gè)無(wú)性系可作為優(yōu)良材料用于杉木種質(zhì)資源開(kāi)發(fā)、保存和利用。

關(guān)鍵詞： 杉木;球果;種子;品質(zhì)性狀;遺傳多樣性;綜合評(píng)價(jià)

中圖分類號(hào)： S791.27? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼： A 文章編號(hào)：2095-1191（2021）102805-09

Abstract：【Objective】The genetic diversity of cone and seed quality traits of different Cunninghamia lanceolata germplasms was clarified， and the high quality resources were screened to provide scientific basis for the breeding of C. lanceolata varieties and the efficient utilization of germplasm resources. 【Method】The cones and seeds of 40 clones in the second-generation seed orchard of C. lanceolata in Guangxi were selected as the research objects in this study. Using variance analysis， cluster analysis， principal component analysis， 15 quality traits from three aspects（cone phenotypic traits， seed germination characteristics and seed excellence traits） were analyzed and comprehensively evaluated.【Result】The genetic diversity index of 15 quality traits of 40 C. lanceolata clones ranged from 1.4792 to 2.0803， with an average of 1.9090. The largest genetic diversity index was empty seed rate， and the smallest was moldy seed rate. The coefficient of variation ranged from 4.55% to 120.72%， with an average of 29.60%. The largest coefficient of variation was moldy seed rate， and the smallest was fruit shape index. The heritability was between 0.366 and 0.969， with an average value of 0.843. The largest heritability was cone horizontal diameter， and the smallest was empty seed rate. Empty seed rate had significant differences among clones（P<0.05）， the other 14 quality traits had extremely significant differences among clones（P<0.01）. The 40 clones of C. lanceolata could be divided into four groups. Among them， group II had the best comprehensive traits， and group I had the best seed excellence traits. Principal component analysis was conducted on the quality traits of 15 C. lanceolata cones and seeds. The top eight clones with comprehensive traits were numbered as 14， 34， 28， 31， 38， 9， 29 and 27， respectively.【Conclusion】The genetic diversity of cone and seed quality traits of 40 clones in second-generation seed orchard of C. lanceolata in Guangxi is rich， and 8 clones screened can be used as excellent materials for the development， preservation and utilization of germplasm resources.

Key words：Cunninghamia lanceolata; cone;seed;quality traits; genetic diversity; comprehensive evaluation

Foundation item：Fundamental Research Funds of Central Public Welfare Research Innstitues of China Acadmy of Forestry（CAFYBB2017A001-1-4） Funding Project of Guangxi Key Laboratory of Excellent Timber Forest Resources Cultivation（2020-A-01-03）;Guangxi Key Science and Technology Project（Guike AA17204087-2）;Fundamental Research Funds of Guangxi Forestry Research Institute（Linke 202111）

0 引言

【研究意義】杉木[Cunninghamia lanceolata（Lamb.） Hook.]是我國(guó)特有的樹(shù)種，也是南方重要的商品用材樹(shù)種之一，具有較高的經(jīng)濟(jì)和生態(tài)價(jià)值，為全國(guó)林業(yè)發(fā)展做出了重要貢獻(xiàn)（陳代喜等，2015;黃開(kāi)勇等，2015）。據(jù)第九次全國(guó)森林資源清查結(jié)果顯示，廣西人工林面積733.53萬(wàn)ha，森林覆蓋率60.17%，分別位居全國(guó)第1位和第3位。廣西作為國(guó)內(nèi)杉木中心產(chǎn)區(qū)，現(xiàn)有杉木林面積150.38萬(wàn)ha，位列全區(qū)主要人工林樹(shù)種前列，十三五規(guī)劃期間，杉木人工林面積年均增加近5.33萬(wàn)ha，年需良種苗木逾1.4億株（陳代喜等，2018;覃歡蘭等，2020）。廣西林業(yè)科學(xué)研究院在區(qū)內(nèi)營(yíng)建杉木種子園超800 ha，為杉木產(chǎn)業(yè)化發(fā)展提供種質(zhì)保障（戴俊等，2019;程琳等，2020）。隨著人民生活水平的提高，社會(huì)對(duì)木材需求量的不斷增加，近年來(lái)市場(chǎng)上的杉木優(yōu)良種苗依然供不應(yīng)求。種子園球果和種子品質(zhì)是良種壯苗的保障，也是林木達(dá)到優(yōu)質(zhì)優(yōu)產(chǎn)的基礎(chǔ)，對(duì)提高種質(zhì)利用率及推動(dòng)良種化進(jìn)程具有重要意義?！厩叭搜芯窟M(jìn)展】遺傳多樣性是種質(zhì)評(píng)價(jià)的重要因子，對(duì)于優(yōu)良種質(zhì)選育、種質(zhì)創(chuàng)新和資源的開(kāi)發(fā)利用具有重要意義（齊明等，2018;陳越等，2019;丁明亮等，2020）。目前，國(guó)內(nèi)外很多學(xué)者已通過(guò)遺傳多樣性研究方法對(duì)桉（葉紹明等，2010）、松（馮源恒等，2019）、杉（李霞等，2020）等樹(shù)種進(jìn)行了林木對(duì)比試驗(yàn)，并闡述了遺傳變異規(guī)律，利于下一步種質(zhì)資源的保存與利用。陳曉陽(yáng)等（1996，1997）對(duì)杉木種子園球果和種子進(jìn)行解剖性狀和品質(zhì)性狀的研究，探討球果和種子表型多樣性規(guī)律，為良種選育提供了科學(xué)依據(jù)。張卓文（2004）基于表型性狀對(duì)1.5代杉木種子園球果和種子進(jìn)行了多樣性分析，綜合球果產(chǎn)量和單個(gè)球果出籽率2個(gè)指標(biāo)篩選出2個(gè)高產(chǎn)無(wú)性系7415和7431，并揭示了無(wú)性系組成對(duì)種子園球果和種子品質(zhì)的重要性。黃開(kāi)勇等（2015）通過(guò)對(duì)母樹(shù)截桿的方式提高了種子園球果和種子品質(zhì)，球果重增益60%以上，為種子園高產(chǎn)開(kāi)發(fā)提供了科學(xué)方法。鄭新華（2015）分析3代杉木種子園種子質(zhì)量，選擇出總體表現(xiàn)較好的5個(gè)家系作為種子園入園材料，豐富了杉木種質(zhì)資源。瞿虹等（2020）研究得出，杉木2代種子園球果采摘時(shí)間與種子品質(zhì)緊密相關(guān)，采收球果的最佳時(shí)間為每年10月24日之后?！颈狙芯壳腥朦c(diǎn)】球果和種子品質(zhì)性狀多樣性是進(jìn)行遺傳多樣性研究最重要的組成部分，而表型性狀、發(fā)芽特性及優(yōu)良程度是研究球果和種子品質(zhì)多樣性重要的評(píng)定指標(biāo)，但關(guān)于2代杉木球果和種子表型性狀、發(fā)芽特性、優(yōu)良程度的遺傳多樣性分析及種質(zhì)篩選的綜合性研究至今尚未有報(bào)道。【擬解決的關(guān)鍵問(wèn)題】通過(guò)對(duì)廣西2代杉木種子園40個(gè)無(wú)性系的球果和種子品質(zhì)性狀進(jìn)行測(cè)定和分析，明確種子園球果和種子品質(zhì)性狀的遺傳多樣性和遺傳變異程度，同時(shí)通過(guò)綜合評(píng)價(jià)篩選出具有綜合優(yōu)勢(shì)的無(wú)性系，為進(jìn)一步開(kāi)展杉木種質(zhì)資源保護(hù)利用和遺傳改良工作提供優(yōu)良種質(zhì)材料，為杉木產(chǎn)業(yè)化進(jìn)程提供科學(xué)理論依據(jù)。

1 材料與方法

1. 1 試驗(yàn)地概況

試驗(yàn)地位于廣西融水自治縣國(guó)家杉木良種基地（東經(jīng)109°16′，北緯25°05′），地處杉木適生區(qū)，海拔300～1000 m，地貌以中、低山地為主。屬溫潤(rùn)亞熱帶季風(fēng)氣候，年平均氣溫18.8 ℃，年均降水量1824.8 mm，雨量充沛，氣溫適中。

1. 2 球果和種子品質(zhì)檢測(cè)

供試40個(gè)杉木無(wú)性系均由廣西融水自治縣國(guó)家杉木良種基地2代杉木種子園提供，無(wú)性系編號(hào)、來(lái)源地和種質(zhì)特征見(jiàn)表1。于2019年12月對(duì)40個(gè)杉木無(wú)性系分號(hào)進(jìn)行球果和種子采摘和測(cè)量。試驗(yàn)材料設(shè)40個(gè)無(wú)性系號(hào)，每個(gè)無(wú)性系取40個(gè)球果，用游標(biāo)卡尺測(cè)量球果縱徑和橫徑。同時(shí)，每個(gè)無(wú)性系采用隨機(jī)選種的方法，共選300粒種子，設(shè)3個(gè)重復(fù)，每個(gè)重復(fù)100粒種子;將種子均勻擺放在鋪有濾紙的培養(yǎng)皿中，放入恒溫培養(yǎng)箱（26 ℃）中培養(yǎng)，每日加入去離子水以保證培養(yǎng)皿濕潤(rùn)，當(dāng)連續(xù)5 d均無(wú)種子發(fā)芽時(shí)結(jié)束該無(wú)性系發(fā)芽試驗(yàn);每日統(tǒng)計(jì)發(fā)芽情況。

1. 3 指標(biāo)測(cè)定及計(jì)算

參考胡瑞陽(yáng)等（2018）的方法，對(duì)杉木球果和種子品質(zhì)進(jìn)行遺傳多樣性研究。試驗(yàn)測(cè)定和計(jì)算的指標(biāo)具體如下：

（1）球果表型性狀（Cone phenotypic traits，CPT）：

球果縱徑（Cone vertical diameter，CVD，mm）;

球果橫徑（Cone horizontal diameter，CHD，mm）;

果形指數(shù)（Fruit shape index，F(xiàn)SI）=球果縱徑/球果橫徑;

出籽率（Seed rate，SR，%）=種子重/球果重×100。

（2）種子發(fā)芽特性（Seed germination characte-ristics，SGC）測(cè)定：

發(fā)芽率（Germination rate，GR，%）=發(fā)芽種子數(shù)/試驗(yàn)種子數(shù)×100;

發(fā)芽勢(shì)（Germination energy，GE，%）=發(fā)芽初期7 d內(nèi)正常發(fā)芽種子數(shù)/供試種子數(shù)×100;

發(fā)芽指數(shù)（Germination index，GI）=ΣGi/i（Gi表示第i天發(fā)芽數(shù)，i表示天數(shù)）;

發(fā)芽持續(xù)時(shí)間（Germination lasting time，GLT，d）=種子開(kāi)始發(fā)芽到發(fā)芽完成所需時(shí)間;

發(fā)芽起始時(shí)間（Germination initial time，GIT，d）=發(fā)芽試驗(yàn)開(kāi)始到種子開(kāi)始發(fā)芽所需時(shí)間;

發(fā)芽時(shí)間（Germination time，GT，d）=起始發(fā)芽時(shí)間+發(fā)芽持續(xù)時(shí)間。

（3）種子優(yōu)良性狀（Seed excellence traits，SET）測(cè)定：

優(yōu)良度（Excellent degree，ED，%）=優(yōu)良粒/測(cè)定種子粒數(shù)×100;

千粒重（Thousand-grain weight，TGW，g）：在氣干狀態(tài)下的1000粒純凈種子的重量;

空粒率（Empty seed rate，ESR，%）=空粒數(shù)/測(cè)定種子粒數(shù)×100;

澀粒率（Death seed rate，DSR，%）=澀粒數(shù)/測(cè)定種子粒數(shù)×100;

霉變率（Moldy seed rate，MSR，%）=發(fā)霉種子數(shù)/測(cè)定種子粒數(shù)×100。

（4）遺傳力（Heritability，h2）測(cè)定：

h2=1-1/F

式中，h2為廣義遺傳力，F(xiàn)為方差分析中的F值。

（5）多樣性指數(shù)（H'）測(cè)定：

品質(zhì)性狀參照史京京等（2019）的方法，劃分為10級(jí)，即1級(jí)<（x-2s）至10級(jí)>（x+2s），每0.5s為一級(jí)，x為平均值，s為標(biāo)準(zhǔn)差，每一級(jí)的分布頻率用于計(jì)算遺傳多樣性指數(shù)。

1. 4 統(tǒng)計(jì)分析

采用Excel 2007進(jìn)行數(shù)據(jù)匯總處理和制圖，以SPSS 19.0進(jìn)行方差分析、聚類分析和主成分分析。杉木球果和種子品質(zhì)性狀多樣性分析：采用變異系數(shù)進(jìn)行離散程度分析;采用Shannon-wiener多樣性指數(shù)（H'）研究球果和種子品質(zhì)性狀多樣性程度（謝向譽(yù)等，2017;董勝君等，2020）。

2 結(jié)果與分析

2. 1 杉木球果和種子品質(zhì)性狀多樣性分析

2. 1. 1 品質(zhì)性狀變異分析 對(duì)40個(gè)杉木無(wú)性系品質(zhì)性狀進(jìn)行多樣性分析，由表2可知，杉木無(wú)性系15個(gè)品質(zhì)性狀的變異系數(shù)介于4.55%～120.72%，總均值為29.60%。其中，球果表型性狀、種子發(fā)芽特性和種子優(yōu)良度性狀均值分別為12.03%、28.72%和44.71%。球果縱徑等14個(gè)品質(zhì)性狀在無(wú)性系間均存在極顯著差異（P<0.01，下同），僅空粒率在無(wú)性系間存在顯著差異（P<0.05，下同）。

2. 1. 2 品質(zhì)性狀遺傳變異分析 對(duì)15個(gè)杉木球果和種子無(wú)性系品質(zhì)性狀進(jìn)行遺傳力計(jì)算和分析，結(jié)果（表2）表明，杉木15個(gè)參試性狀的遺傳力在0.366～0.969。

2. 1. 3 品質(zhì)性狀分布頻率分析 通過(guò)對(duì)杉木球果和種子15個(gè)品質(zhì)性狀進(jìn)行變異分析，統(tǒng)計(jì)40個(gè)無(wú)性系品質(zhì)性狀分布頻率，并列出40個(gè)無(wú)性系15個(gè)品質(zhì)性狀的H'（表2）。15個(gè)品質(zhì)性狀的H'范圍為1.4792～2.0803，均值為1.9090。其中，空粒率（2.0803）的H'值最大，霉變率（1.4792）的H'值最小。H'>2.0000的性狀有空粒率、千粒重和球果縱徑，占參試性狀的20%;H'介于1.9000～2.0000的性狀占參試性狀的53%。15個(gè)參試性狀的H'均大于1.4000，表明杉木球果和種子品質(zhì)具有豐富的遺傳多樣性。

2. 2 杉木球果和種子品質(zhì)性狀聚類分析

采用離差平方和法，以歐氏距離為遺傳距離，對(duì)杉木無(wú)性系15個(gè)品質(zhì)性狀進(jìn)行聚類，結(jié)果表明，在歐式距離12處，40份杉木種質(zhì)可分為四大類群（圖1）。其中，類群Ⅰ包括22個(gè)無(wú)性系，占參試無(wú)性系的55%，總體特征為種子優(yōu)良性狀最佳、種子發(fā)芽特性好;類群Ⅱ包括6個(gè)無(wú)性系，占參試無(wú)性系的15%，總體特征為種子發(fā)芽率、發(fā)芽勢(shì)及發(fā)芽指數(shù)等發(fā)芽特性最佳，優(yōu)良性狀好，球果果形指數(shù)最大;類群Ⅲ包括2個(gè)無(wú)性系，占參試無(wú)性系的5%，屬于球果大而均勻、種子優(yōu)良度欠佳的類型;類群Ⅳ包括10個(gè)無(wú)性系，占參試無(wú)性系的25%，屬于種子優(yōu)良性狀較好、發(fā)芽特性欠佳、球果果形指數(shù)最小的類型（表3）。

2. 3 杉木球果和種子品質(zhì)性狀的綜合評(píng)價(jià)和指標(biāo)篩選

2. 3. 1 杉木球果和種子品質(zhì)性狀的相關(guān)分析 由表4可知，球果縱徑與球果橫徑呈極顯著正相關(guān)，與果形指數(shù)呈顯著正相關(guān);球果橫徑與千粒重呈極顯著正相關(guān);出籽率與發(fā)芽起始時(shí)間呈顯著負(fù)相關(guān);發(fā)芽率與發(fā)芽勢(shì)、發(fā)芽指數(shù)、發(fā)芽持續(xù)時(shí)間、發(fā)芽時(shí)間呈極顯著正相關(guān)，與霉變率呈極顯著負(fù)相關(guān);發(fā)芽勢(shì)與發(fā)芽指數(shù)呈極顯著正相關(guān)，與發(fā)芽起始時(shí)間呈極顯著負(fù)相關(guān);發(fā)芽指數(shù)與發(fā)芽持續(xù)時(shí)間呈極顯著正相關(guān)，與發(fā)芽時(shí)間和千粒重呈顯著正相關(guān)，與發(fā)芽起始時(shí)間呈顯著負(fù)相關(guān);發(fā)芽持續(xù)時(shí)間與發(fā)芽時(shí)間呈極顯著正相關(guān)，與霉變率呈顯著負(fù)相關(guān);發(fā)芽時(shí)間與霉變率呈極顯著負(fù)相關(guān);優(yōu)良度與澀粒率呈極顯著負(fù)相關(guān)?？梢?jiàn)，各品質(zhì)指標(biāo)間有一定的相關(guān)性，表明篩選出球果表型、種子發(fā)芽特性及優(yōu)良性狀兼優(yōu)的無(wú)性系存在可能性。

2. 3. 2 杉木球果和種子品質(zhì)性狀的綜合評(píng)價(jià) 40個(gè)杉木無(wú)性系15個(gè)品質(zhì)性狀的主成分分析結(jié)果見(jiàn)表5。前6個(gè)主成分的累積貢獻(xiàn)率達(dá)82.767%，包含不同無(wú)性系品質(zhì)性狀的絕大部分信息。第1主成分（PC1）特征值為3.745，貢獻(xiàn)率達(dá)24.968%，其中貢獻(xiàn)率最大的是發(fā)芽率，其次是發(fā)芽指數(shù);第2主成分（PC2）特征值為2.632，貢獻(xiàn)率達(dá)17.549%，其中貢獻(xiàn)率最大的是球果橫徑，其次是球果縱徑;第3主成分（PC3）特征值為2.029，貢獻(xiàn)率達(dá)13.529%，其中貢獻(xiàn)率最大的是發(fā)芽起始時(shí)間，其次是發(fā)芽時(shí)間;第4主成分（PC4）特征值為1.504，貢獻(xiàn)率達(dá)10.026%，其中貢獻(xiàn)率最大的是優(yōu)良度;第5主成分（PC5）特征值為1.400，貢獻(xiàn)率達(dá)9.335%，其中貢獻(xiàn)率最大的是果形指數(shù);第6主成分（PC6）特征值為1.104，貢獻(xiàn)率達(dá)7.361%，其中貢獻(xiàn)率最大的是出籽率。

將40個(gè)杉木無(wú)性系15個(gè)品質(zhì)性狀代入6個(gè)主成分，可得到每個(gè)無(wú)性系6個(gè)主成分的得分。以各主成分所占總貢獻(xiàn)率的比例作為權(quán)重系數(shù)（0.302、0.212、0.163、0.121、0.113、0.089），計(jì)算得到每個(gè)無(wú)性系球果和種子品質(zhì)性狀的綜合得分（U值）及排名（表6）。40個(gè)杉木無(wú)性系的綜合得分介于19.245～47.365，均值為31.786;按20%的入選率選出綜合排序前8的無(wú)性系編號(hào)，依次為14、34、28、31、38、9、29、27，篩選出的無(wú)性系綜合得分均值為41.043。

2. 3. 3 杉木球果和種子品質(zhì)性狀評(píng)價(jià)指標(biāo)的篩選 對(duì)綜合得分（U值）與15個(gè)品質(zhì)性狀進(jìn)行相關(guān)性分析（表7），結(jié)果表明，15個(gè)品質(zhì)性狀中，發(fā)芽特性與U值的相關(guān)性最高，其次是種子優(yōu)良程度，球果表型性狀與U值無(wú)顯著相關(guān)性。其中，U值與發(fā)芽率、發(fā)芽指數(shù)、發(fā)芽持續(xù)時(shí)間、發(fā)芽時(shí)間和發(fā)芽勢(shì)均呈極顯著正相關(guān)，與千粒重呈顯著正相關(guān)，與霉變率呈極顯著負(fù)相關(guān)。因此，綜合得分可作為球果和種子品質(zhì)綜合評(píng)價(jià)的重要指標(biāo)。

3 討論

3. 1 杉木無(wú)性系品質(zhì)性狀的遺傳差異

球果和種子品質(zhì)綜合評(píng)價(jià)的最終目標(biāo)是推廣利用，表型性狀是遺傳多樣性分析最直接的指標(biāo)，發(fā)芽特性是種子萌芽能力的重要評(píng)價(jià)指標(biāo)，優(yōu)良性狀是質(zhì)量評(píng)價(jià)的重要因子（周增亮等，2020）。本研究中15個(gè)杉木品質(zhì)性狀的遺傳力、變異系數(shù)和多樣性指數(shù)代表了40個(gè)無(wú)性系品質(zhì)的多樣性差異，15個(gè)品質(zhì)性狀遺傳力介于0.366～0.969，最大的是球果橫徑;變異系數(shù)范圍為4.55%～120.72%，最大的為霉變率;多樣性指數(shù)在1.4792～2.0803，空粒率、千粒重和球果長(zhǎng)的多樣性指數(shù)高（H'>2），有80%性狀的多樣性指數(shù)高于1.80。杉木種子發(fā)芽率（h2=0.902）和澀粒率（h2=0.918）遺傳力較高，受遺傳控制強(qiáng)，具有較好的選育效果，與陳仕昌（2019）研究一致。種子霉變率的變異系數(shù)最大，表明霉變程度受環(huán)境影響大;種子發(fā)芽率變異程度高于優(yōu)良度，與陳琴等（2019）、陳仕昌等（2019）的研究結(jié)果相同。本研究表明，40個(gè)杉木無(wú)性系遺傳多樣性豐富，可用于杉木種質(zhì)資源開(kāi)發(fā)和利用。

3. 2 杉木無(wú)性系品質(zhì)性狀的聚類分析

利用聚類分析將40個(gè)杉木無(wú)性系分為四大類群。其中，類群Ⅰ的種子優(yōu)良性狀最優(yōu);類群Ⅱ的種子發(fā)芽特性、優(yōu)良度和球果表型性狀俱佳;類群Ⅲ的球果大而均勻;類群Ⅳ的種子優(yōu)良度較好。來(lái)自廣西的無(wú)性系（編號(hào)14、28、31、34、38）和來(lái)自貴州的9號(hào)聚在類群Ⅱ，說(shuō)明種質(zhì)來(lái)源與品質(zhì)具有一定相關(guān)性，與NGUYEN（2016）的研究結(jié)果相同，可為種質(zhì)篩選和資源利用提供理論參考。

3. 3 杉木無(wú)性系品質(zhì)性狀的綜合評(píng)價(jià)

本研究采用主成分分析法把15個(gè)品質(zhì)性狀轉(zhuǎn)換為6個(gè)主成分，代表大部分指標(biāo)信息，其中，發(fā)芽率、球果橫徑、發(fā)芽起始時(shí)間、優(yōu)良度、果形指數(shù)和出籽率是主要構(gòu)成因子，可為綜合評(píng)價(jià)提供評(píng)定參考。結(jié)果顯示，綜合排序前8的無(wú)性系編號(hào)依次為14、34、28、31、38、9、29、27，可為進(jìn)一步種質(zhì)選育提供科學(xué)依據(jù)。

研究杉木球果和種子遺傳多樣性，探索遺傳變異豐富度，掌握其變異規(guī)律和特點(diǎn)，對(duì)于杉木優(yōu)良種質(zhì)資源的收集保存、挖掘利用及合理開(kāi)發(fā)具有極其重要的意義（李明等，2016）。種子園球果和種子品質(zhì)是有效提高林木產(chǎn)量和質(zhì)量的前提，各相關(guān)指標(biāo)間的聯(lián)系則為良種選育提供可能（Lambeth et al.，2001）。分子生物技術(shù)也逐漸應(yīng)用于杉木種質(zhì)多樣性分析，目前已有多種分子標(biāo)記應(yīng)用于杉木（王新民等，2014;文亞峰等，2015）。本研究選取廣西2代杉木種子園40個(gè)杉木無(wú)性系為試材，從球果表型性狀、種子發(fā)芽特性、種子優(yōu)良性狀進(jìn)行分析，下一步需要從分子水平對(duì)杉木種質(zhì)進(jìn)行深入分析，同時(shí)針對(duì)不同需求對(duì)杉木種質(zhì)進(jìn)行擴(kuò)繁和評(píng)價(jià)，選育出杉木新品種。

4 結(jié)論

40個(gè)杉木球果和種子的4個(gè)球果表型性狀、6個(gè)種子發(fā)芽特性和5個(gè)種子優(yōu)良性狀變異程度較高，遺傳多樣性豐富;其中，千粒重、球果縱徑、發(fā)芽率具有較強(qiáng)的遺傳變異性，且各性狀間具有一定相關(guān)性，可為杉木球果和種子品質(zhì)評(píng)定提供重要的參考依據(jù)。同時(shí)，篩選出的8個(gè)杉木無(wú)性系可為種子園營(yíng)建及良種繁育提供優(yōu)質(zhì)種質(zhì)材料。

參考文獻(xiàn)：

陳代喜，陳琴，蒙躍環(huán)，黃開(kāi)勇，黃海仲，劉凡勝，韋家國(guó)，莫亮榮. 2015. 杉木大徑材高效培育技術(shù)探討[J]. 南方農(nóng)業(yè)學(xué)報(bào)，46（2）： 293-298. [Chen D X，Chen Q，Meng Y H，Huang K Y，Huang H Z，Liu F S，Wei J G，Mo L R. 2015. High-efficient management technology of big-dia-meter timber of Cunninghamia lanceolata[J]. Journal of Southern Agriculture，46（2）：293-298.] doi：10.3969/jissn. 2095-1191.2015.2.293.

陳代喜，程琳，黃開(kāi)勇，戴俊，羅啟亮，莫宗恒，陳曉明，藍(lán)肖，唐紅亮. 2018. 30年生杉木半同胞子代遺傳特性測(cè)定及優(yōu)良家系和單株選擇[J]. 南方農(nóng)業(yè)學(xué)報(bào)，49（2）：307-312. [Chen D X，Cheng L，Huang K Y，Dai J，Luo Q L，Mo Z H，Chen X M，Lan X，Tang H L. 2018. Genetic characterization test and selection of superior family and individual of 30-year-old Chinese fir half-sib progeny[J]. Journal of Southern Agriculture，49（2）：307-312.] doi：10. 3969/j.issn.2095-1191.2018.02.16.

陳琴，董利軍，黃開(kāi)勇，戴俊，譚文婧，程琳，李魁鵬，黃鴻飛. 2019. 紅心杉種子園不同無(wú)性系種子品質(zhì)差異性分析[J]. 廣西林業(yè)科學(xué)，48（2）：152-156. [Chen Q，Dong L J，Huang K Y，Dai J，Tan W J，Cheng L，Li K P，Huang H F. 2019. Analysis on seed quality of different clones in red-heart Chinese fir seed orchard[J]. Guangxi Forest Science，48（2）：152-156.] doi：10.19692/j.cnki.gfs.2019.02.003.

陳仕昌. 2019. 1.5代杉木種子園種實(shí)差異及遺傳多樣性研究[D]. 南寧： 廣西大學(xué). [ Chen S C. 2019. Study on cone and seed difference and genetic diversity in 1.5 generation Cunninghamia lanceolata seed orchard[D]. Nanning： Guangxi University. ]

陳仕昌，李魁鵬，黃開(kāi)勇，戴俊，董利軍，譚文婧，梁機(jī). 2019. 1.5代杉木種子園不同無(wú)性系結(jié)實(shí)及種子差異[J]. 廣西林業(yè)科學(xué)，48（1）：27-33. [Chen S C，Li K P，Huang K Y，Dai J，Dong L J，Tan W J，Liang J. 2019. Different clone fruiting and seed differences in 1.5-generation seed orchard of Cunninghamia lanceolata[J]. Guangxi Forest Science，48（1）：27-33.] doi：10.3969/j.issn.1006-1126. 2019.01.005.

陳曉陽(yáng)，李文剛，潘奇敏，張定鉞. 1996. 杉木種子園球果解剖性狀的研究[J]. 北京林業(yè)大學(xué)學(xué)報(bào)，18（3）：42-47. [Chen X Y，Li W G，Pan Q M，Zhang D Y. 1996. The studies on anatomical traits of cones in Chinese-fir seed orchard[J]. Journal of Beijing Forestry University，18（3）：42-47.] doi：10.13332/j.1000-1522.1996.03.007.

陳曉陽(yáng)，潘奇敏，劉朝祿，楊明勝. 1997. 杉木種子園種子品質(zhì)的觀測(cè)和分析[J]. 北京林業(yè)大學(xué)學(xué)報(bào)，19（1）：32-38. [Chen X Y，Pan Q M，Liu C L，Yang M S. 1997. Observation and analysis on the traits of seeds from Chinese-fir seed orchards[J]. Journal of Beijing Forestry University，19（1）：32-38.]

陳越，丁明亮，張敦宇，付堅(jiān)，鐘巧芳，肖素勤，柯學(xué)，程在全. 2019. 云南水稻種質(zhì)資源農(nóng)藝性狀表型多樣性分析及綜合評(píng)價(jià)[J]. 南方農(nóng)業(yè)學(xué)報(bào)，50（9）：1922-1930. [Chen Y，Ding M L，Zhang D Y，F(xiàn)u J，Zhong Q F，Xiao S Q，Ke X，Cheng Z Q. 2019. Phenotypic diversity analysis and comprehensive evaluation of agronomic characters of rice germplasm resources in Yunnan Province[J]. Journal of Southern Agriculture，50（9）：1922-1930.] doi：10.3969/j.issn.2095-1191.2019.09.05.

程琳，陳代喜，黃開(kāi)勇，戴俊，陳曉明，李魁鵬，譚文婧，董利軍，唐紅亮，覃孟哲，韋連尤，王家妍. 2020. 福建第3代杉木種子園子代在廣西造林的早期表現(xiàn)[J]. 廣西林業(yè)科學(xué)，49（1）：113-116. [Cheng L，Chen D X，Huang K Y，Dai J，Chen X M，Li K P，Tan W J，Dong L J，Tang H L，Qin M Z，Wei L Y，Wang J Y. 2020. Early afforestation performance of Fujian third generation seed orchard pro-geny of Cunninghamia lanceolata in Guangxi[J]. Guangxi Forest Science，49（1）：113-116.] doi：10.19692/j.cnki.gfs. 2020.01.021.

戴俊，譚文婧，黃開(kāi)勇，陳琴，李魁鵬，董利軍，黃鴻飛，梁丹鳳. 2019. 紅心杉不同無(wú)性系結(jié)實(shí)特性及球果類型分析[J]. 廣西林業(yè)科學(xué)，48（2）：141-146. [Dai J，Tan W J，Huang K Y，Chen Q，Li K P，Dong L J，Huang H F，Liang D F. 2019. Analysis on fruiting characteristics and cone types of different clones of red-heart Chinese fir[J]. Guang-xi Forest Science，48（2）：141-146.] doi：10.19692/j.cnki.gfs.2019.02.001.

丁明亮，林麗萍，李明菊，胡欣，何遲，廖合發(fā)，趙紅，李紹祥，劉琨. 2020. 云南育成小麥品種（系）品質(zhì)性狀遺傳多樣性分析及綜合評(píng)價(jià)[J]. 南方農(nóng)業(yè)學(xué)報(bào)，51（2）：255-266. [Ding M L，Lin L P，Li M J，Hu X，He C，Liao H F，Zhao H，Li S X，Liu K. 2020. Genetic diversity analysis and comprehensive evaluation of quality traits of wheat varieties（lines）bred in Yunnan[J]. Journal of Southern Agriculture，51（2）：255-266.] doi：10.3969/j.issn.2095-1191.2020.02.002.

董勝君，孫永強(qiáng)，陳建華，盧彩云，劉權(quán)鋼，劉立新. 2020. 野杏無(wú)性系表型性狀多樣性分析及綜合評(píng)價(jià)[J]. 植物遺傳資源學(xué)報(bào)，21（5）：1156-1166. [Dong S J，Sun Y Q，Chen J H，Lu C Y，Liu Q G，Liu L X. 2020. Phenotypic traits diversity analysis and comprehensive evaluation of Armeniaca vulgaris var. ansu clones[J]. Journal of Plant Genetic Resources，21（5）：1156-1166.] doi：10.13430/j.cnki.jpgr.20200305001.

馮源恒，李火根，楊章旗，黃永利，陳虎. 2019. 廣西馬尾松三個(gè)優(yōu)良種源的遺傳多樣性及生長(zhǎng)性狀變異分析[J].南京林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），43（6）：67-72. [Feng Y H，Li H G，Yang Z Q，Huang Y L，Chen H. 2019. Analysis of genetic diversity and variation of growth traits of three superior provenances of Pinus massoniana（masson pine） in Guangxi[J]. Journal of Nanjing Forestry University（Natural Sciences Edition），43（6）：67-72.] doi：10. 3969/j. issn.1000-2006.201806039.

胡瑞陽(yáng)，孫宇涵，吳博，段紅靜，林華忠，方祿明，余小龍，李云. 2018. NaN3處理對(duì)杉木種子發(fā)芽及幼苗生長(zhǎng)的影響[J]. 東北林業(yè)大學(xué)學(xué)報(bào)，46（5）：6-11. [Hu R Y，Sun Y H，Wu B，Duan H J，Lin H Z，F(xiàn)ang L M，Yu X L，Li Y. 2018. Effect of NaN3 on Cunninghamia lanceolata seed germination and seedlings growth[J]. Journal of Northeast Forestry University，46（5）：6-11.] doi：10.13759/j.cnki. dlxb.2018.05.002.

黃開(kāi)勇，陳琴，唐文，戴俊，蔣華，黃鵬艷. 2015. 截桿處理對(duì)大齡杉木種子園種子產(chǎn)量與品質(zhì)的影響研究[J]. 西部林業(yè)科學(xué)，44（2）：29-35. [Huang K Y，Chen Q，Tang W，Dai J，Jiang H，Huang P Y. 2015. The influences of trunk-cutting treatment on seed yield and quality from elder Cunninghamia lanceolata seed orchard[J]. Journal of West China Forestry Science，44（2）：29-35.]

李明，黃木生，蘇爍爍，陳細(xì)貞，馬祥慶. 2016. 福建杉木王種質(zhì)資源的遺傳多樣性[J]. 森林與環(huán)境學(xué)報(bào)，36（3）：312-318. [Li M，Huang M S，Su S S，Chen X Z，Ma X Q. 2016. Genetic diversity of germplasm resources of the king of Chinese fir in Fujian provenances[J]. Journal of Forest and Environment，36（3）：312-318.] doi：10.13324/j.cnki.jfcf.2016.03.010.

李霞，王利寶，文亞峰，林軍，武星彤，袁美靈，張?jiān)趺羟?，李鑫? 2020. 杉木不同世代育種群體的遺傳多樣性[J]. 林業(yè)科學(xué)，56（11）：53-61. [Li X，Wang L B，Wen Y F，Lin J，Wu X T，Yuan M L， Zhang Y，Wang M Q，Li X Y. 2020. Genetic diversity of Chinese fir（Cunninghamia lanceolata） breeding populations among different generations[J]. Scientia Silvae Sinicae，56（11）：53-61.] doi：10.11707/j.1001-7488.20201106.

NGUYEN BA TRIEU（阮伯趙）. 2016. 越南杉木優(yōu)良種源早期選擇與中國(guó)杉木種源的比較研究[D]. 福州：福建農(nóng)林大學(xué). [NGUYEN BA TRIEU. 2016. Research on selection of good Chinese-fir provenances in Vietnam and compare with the provenances in China[D]. Fuzhou： Fujian Agriculture and Forestry University.]

齊明，王海蓉，彭九生，周建革，翁春媚，何貴平. 2018. 杉木育種園的遺傳多樣性和雜種優(yōu)勢(shì)群的劃分研究[J]. 南方林業(yè)科學(xué)，46（6）：17-21. [Qi M，Wang H R，Peng J S，Zhou J G，Weng C M，He G P. 2018. Genetic diversity and heterotic grouping of Cunninghamia lanceolata bree-ding garden[J]. South China Forestry Science，46（6）：17-21.] doi：10.16259/j.cnki.36-1342/s.2018.06.004.

覃歡蘭，程琳，池勇權(quán)，陳曉璋，甘炳幸. 2020. 沙床育苗技術(shù)對(duì)杉木苗木生長(zhǎng)影響的研究[J]. 林業(yè)科技，45（3）：1-4. [Qin H L，Cheng L，Chi Y Q，Chen X Z，Gan B X. 2020. Study on the effect of sand-bed seedling technology on seedling growth of Cunninghamia lanceolata seedlings[J]. Forestry Science & Technology，45（3）：1-4.] doi：10.19750/j.cnki.1001-9499.2020.03.001.

瞿虹，肖紀(jì)軍，季如春，曾甜甜，金霖芳. 2020. 杉木種子園球果采收時(shí)間與種子品質(zhì)的關(guān)系[J]. 園藝與種苗，40（4）：22-23. [Qu H，Xiao J J，Ji R C，Zeng T T，Jin L F. 2020. Relationship between cone harvest time and seed quality in Chinese fir seed orchard[J]. Horticulture & Seed，40（4）：22-23.] doi：10.16530/j.cnki.cn21-1574/s.2020.04.009.

史京京，薛盈文，郭偉，于崧，陸旺，于立河，胡慧影，郭燦，侯楠. 2019. 黑龍江西部地區(qū)引進(jìn)的飼用燕麥種質(zhì)資源遺傳多樣性分析及綜合評(píng)價(jià)[J]. 南方農(nóng)業(yè)學(xué)報(bào)，50（3）：515-523. [Shi J J，Xue Y W，Guo W，Yu S，Lu W，Yu L H，Hu H Y，Guo C，Hou N. 2019. Genetic diversity analysis and comprehensive evaluation of introduced Avena sativa L. germplasm resources in the western area of Heilongjiang[J]. Journal of Southern Agriculture，50（3）：515-523.] doi：10.3969/j.issn.2095-1191.2019.03.11.

王新民，鄭仁華，陳金慧，趙亞琦，徐陽(yáng)，王穎，施季森. 2014. 杉木SRAP-PCR體系優(yōu)化[J]. 南京林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），38（1）：15-20. [Wang X M，Zheng R H，Chen J H，Zhao Y Q，Xu Y， Wang Y， Shi J S. 2014. Optimization of Chinese fir SRAP-PCR system[J]. Journal of Nanjing Forestry University（Natural Sciences Edition），38（1）：15-20.] doi：10.3969/j.issn.1000-2006.2014.01.003.

文亞峰，韓文軍，周宏，徐剛標(biāo). 2015. 杉木轉(zhuǎn)錄組SSR挖掘及EST-SSR標(biāo)記規(guī)?；_(kāi)發(fā)[J]. 林業(yè)科學(xué)，51（11）：40-49. [Wen Y F，Han W J，Zhou H，Xu G B. 2015. SSR mining and development of EST-SSR markers for Cunninghamia lanceolata based on transcriptome sequences[J]. Scientia Silvae Sinicae，51（11）：40-49.] doi：10.11707/j.1001-7488. 20151106.

謝向譽(yù)，陸柳英，曾文丹，賴大欣，嚴(yán)華兵. 2017. 31份木薯種質(zhì)資源的鑒定評(píng)價(jià)及遺傳多樣性分析[J]. 南方農(nóng)業(yè)學(xué)報(bào)，48（3）：393-400. [Xie X Y，Lu L Y，Zeng W D，Lai D X，Yan H B. 2017. Identification， evaluation and genetic diversity analysis of 31 cassava germplasm resour-ces[J]. Journal of Southern Agriculture，48（3）：393-400.] doi：10.3969/j：issn.2095-1191.2017.03.003.

葉紹明，溫遠(yuǎn)光，楊梅，梁宏溫. 2010. 連栽桉樹(shù)人工林植物多樣性與土壤理化性質(zhì)的關(guān)聯(lián)分析[J]. 水土保持學(xué)報(bào)，24（4）：246-250. [Ye S M，Wen Y G，Yang M，Liang H W. 2010. Correlation analysis on biodiversity and soil physical & chemical properties of Eucalyptus spp. plantations under successive rotation[J]. Journal of Soil and Water Conservation，24（4）：246-250.] doi：10.13870/j.cnki.stbcxb.2010.04.051.

張卓文. 2004. 杉木種子園無(wú)性系球果及其種子分布特征[J]. 福建林學(xué)院學(xué)報(bào)，24（1）：45-49. [Zhang Z W. 2004. The female cone and its seed distribution among clones of Chinese fir in Chongyang seed orchard[J]. Journal of Fujian college of forestry，24（1）：45-49.] doi：10.13324/j.cnki.jfcf.2004.01.012.

鄭新華. 2015. 杉木第3代種子園不同家系種子質(zhì)量評(píng)價(jià)[J]. 亞熱帶農(nóng)業(yè)研究，11（2）：95-100. [Zheng X H. 2015. Eva-luation of different families’ seed quality in third-generation seed orchard of Cunninghamia lanceolata[J]. Subtropical Agriculture Research，11（2）：95-100.] doi：10. 13321/j.cnki.subtrop.agric.res.2015.02.005.

周增亮，胡冬南，晏雨鴻，何萬(wàn)存，陳夢(mèng)義，張旖旎，程離，劉娟. 2020. 馬尾松第二代種子園不同無(wú)性系生長(zhǎng)性狀和種實(shí)性狀差異研究[J]. 河南農(nóng)業(yè)大學(xué)學(xué)報(bào)， 54（4）： 566-574. [Zhou Z L， Hu D N，Yan Y H，He W C， Chen M Y，Zhang Y N，Cheng L，Liu J. 2020. Study on the differences of growth and seed characters among differentclones in the second-generation seed orchard of Pinus massoniana[J]. Journal of Henan Agricultural University， 54（4）： 566-574.] doi：10.16445/j.cnki.1000-2340.2020.04. 024.

Lambeth C，Lee B C，O’Malley，Wheeler N. 2001. Polymix breeding with parental analysis of progeny：An alternative to full-sib breeding and testing[J]. Theoretical Applied Genetics，103（6-7）：903-943. doi：10.1007/s001220100627.

（責(zé)任編輯 鄧慧靈）