摘 要：【目的】通過分析38個滇山茶Camellia reticulata品種的12個花表型性狀在數(shù)量、大小和形狀三個維度上的分布和差異，探討滇山茶不同品種花表型變異的（橫向和縱向）規(guī)律，為滇山茶新種質(zhì)的選育和高效利用提供參考?！痉椒ā繉Φ嵘讲璧?2個花表型性狀進(jìn)行方差分析、相關(guān)性分析、主成分分析、綜合分析和聚類分析，箱線圖和應(yīng)用頻率分布函數(shù)等方法進(jìn)行分析?！窘Y(jié)果】滇山茶不同品種花的12個表型性狀存在豐富的變異，變異系數(shù)為13.00%～67.00%，平均變異系數(shù)為27.24%；變異系數(shù)法的權(quán)重計算結(jié)果顯示，雌蕊數(shù)（20.485%）的變異性影響最大，花瓣長寬比（3.99%）的變異性影響最?。幌嚓P(guān)性分析結(jié)果表明，12個性狀間都存在極顯著相關(guān)關(guān)系，花徑與花瓣寬度相關(guān)性最大（r=0.822）；主成分分析結(jié)果表明，前4個主成分提取到的特征值都均大于1.00，累計貢獻(xiàn)率達(dá)到82.129%；通過聚類分析將其分為3大類。綜合評價結(jié)果表明‘大理茶’的綜合得分最高（F=3.230），在38個滇山茶品種中綜合特性最好。對聚類分析中的3類滇山茶進(jìn)行箱形圖和頻率分布函數(shù)分析，其中9個性狀表現(xiàn)出顯著差異，其余3個性狀無顯著差異，第Ⅰ類各性狀分布較向上，第Ⅱ類分布較為居中，第Ⅲ類分布較向下。而在頻率分布函數(shù)分析中，各類的性狀均表現(xiàn)出明顯的偏置，變化程度依次為數(shù)量相關(guān)特征>尺寸相關(guān)特征>形狀相關(guān)特征?！窘Y(jié)論】滇山茶的12個花表型性狀存在豐富的變異現(xiàn)象，主要以雌雄蕊群逐漸瓣化為演化方向，研究結(jié)果確定了滇山茶花表型演化的趨勢，為后期滇山茶定向育種提供一定的理論依據(jù)。

關(guān)鍵詞：滇山茶；花表型性狀；不同品種；多樣性分析

中圖分類號：S794.4 文獻(xiàn)標(biāo)志碼：A 文章編號：1673-923X（2024）10-0181-15

基金項目：國家自然科學(xué)基金項目（32260416）；云南省教育廳科學(xué)研究基金項目（2022J0520）；云南省農(nóng)業(yè)基礎(chǔ)研究聯(lián)合專項項目（202101BD070001-095）；云南省“興滇英才支持計劃”青年人才專項項目（YNQR-QNRC-2019-104）。

Diversity analysis of flowers phenotypic traits of Camellia reticulata cultivars

XIAO Hongxinga，b，c， MIN Jiachena，b，c， WU Xueqina，b，c， HUANG Dia，b，c， CHEN Shengtonga，b，c， CHEN Longqinga，b，c， GENG Fanga，b，c

（a.College of Landscape Architecture and Horticulture； b. Southest Landscape Architecture Engineering Technology Research Center， National Forestry and Grassland Administration； c. Yunnan Functional Flower Resources and Industrialization Technology Engineering Research Center， Southwest Forestry University， Kunming 650224， Yunnan， China）

Abstract：【Objective】Phenotypic variation （including directions and degrees） of 12 floral traits from 38 cultivars of Camellia reticulata were explored by comparing their distributions and deviations in three different dimensions： floral organ number， size， and the shape.【Method】Analysis of variance， correlation analysis， principal component analysis， comprehensive analysis， cluster analysis， as well as box-plot and application of frequency distribution functions were conducted on the 12 floral phenotypic traits of C. reticulata.【Result】There were rich variations in 12 floral traits of different cultivars， with a coefficient of variation （CV） ranging from 13.00% to 67.00%， an average CV of 27.24%；The weight calculation using the CV showed that the variability of pistil number （20.485%） had the greatest impact， while the variability of petal length /petal width （3.99%） had the smallest impact； Correlation analysis showed that there was a highly significant correlation among all 12 traits， with the greatest correlation between flower diameter and petal width （r=0.822）. Principal component analysis showed that the feature values extracted from the first four principal components were all greater than 1.00， with a cumulative contribution rate of 82.129%； Cluster analysis divided 38 cultivars into three groups. The comprehensive evaluation indicated that ‘Dali Cha’ received the highest comprehensive score （F=3.230）， making it the best among the 38 cultivars of C. reticulata in terms of comprehensive characteristics. The three groups of 38 cultivars were also analyzed by using box plots and frequency distribution functions. There were only nine traits showed significant differences. Group I of floral traits all showed upward distributed box bodies in box plot analysis and Group III was downward distributed. However， all the 12 traits in three groups exhibited significant differences in frequency distribution functions. The variation degrees were quantitatively characterized by numbering traits>sizing traits>shaping traits. 【Conclusion】The 12 floral phenotypic traits of C. reticulata exhibit rich variability， primarily evolving towards gradual petalization of the stamen and pistil groups. The research findings establish the trend of floral phenotype evolution in C. reticulata， providing a theoretical basis for targeted breeding of C. reticulata in the future.

Keywords： Camellia reticulata； floral phenotypic traits； different cultivars； diversity analysis

滇山茶Camellia reticulata，屬于山茶科Theaceae山茶屬Camellia，常綠喬木或灌木，其觀賞性極高，品種豐富多彩、形態(tài)變異多樣、花型豐富多樣、花色艷而不妖[1-2]，位居云南八大名花之首，同時也是云南省的省花和昆明市、楚雄市和騰沖市的市花[3-4]。滇山茶是西部地區(qū)尤其是云南特有的木本觀花植物，但其應(yīng)用沒有體現(xiàn)資源優(yōu)勢。首先，據(jù)世界山茶屬品種注冊中心記載（截至2023年11月），滇山茶的品種僅有841個品種，遠(yuǎn)遠(yuǎn)低于位居首位的華東山茶C. japonica（18 271個品種）和茶梅C. sasanqua（1 334個品種）。其次，花色是滇山茶重要的觀賞性狀[5]，但滇山茶花色種類相比于其他木本觀花植物，如月季[6-7]、杜鵑花[8]等的花顏色比較單一，主要以紅色系為主，稀有淺色和白色花，更無黃色、藍(lán)色及其他顏色的品種；花型的變異包括單瓣型、半重瓣型、重瓣型，以重瓣型和半重瓣型為主，稀有單瓣型[9-10]。如今，滇山茶雖經(jīng)過長達(dá)1500多年的自然雜交和人工栽培選育，已育出形態(tài)變異多樣，色彩豐富的品種[11-12]，但這些品種的形態(tài)變異規(guī)律及演化趨勢目前鮮有報道。本研究將對滇山茶不同品種的花表型性狀特征進(jìn)行分析，以探索滇山茶花表型演化的趨勢，為后期滇山茶定向育種提供一定的理論依據(jù)。

表型性狀分析是基于形態(tài)學(xué)對植物性狀進(jìn)行分析，為植物的育種研究及分子生物學(xué)等方面提供一定的研究基礎(chǔ)[13]。其花表型性狀是評價植物種質(zhì)資源的重要指標(biāo)，被認(rèn)為是不同分類群復(fù)雜表型鑒定以及評價進(jìn)化與發(fā)育偏差相互作用的傳統(tǒng)性狀[14]，由于花形態(tài)特征所獲得的信息是由大量的定性和定量性狀數(shù)據(jù)集組成的，因此表型多樣性分析被認(rèn)為是對其進(jìn)行評價的最合適的分析工具[15-16]。目前國內(nèi)外山茶屬植物表型研究的報道多集中在杜鵑葉紅山茶C. azlea[17-18]、金花茶C. nitidissima[19]、香花油茶C. osmantha[20]、西南紅山茶C. pitardii[21]等。而關(guān)于滇山茶的研究主要集中在對花型進(jìn)行遺傳多樣性、遺傳關(guān)系的AFLP分析[22]和騰沖紅花油茶C. reticulata f. simplex花芽分化及雌配子體發(fā)育的形態(tài)學(xué)觀察[23]，對花表型性狀研究較少。

本研究通過對選定38個健壯生長、性狀穩(wěn)定的滇山茶品種的花器官表型性狀進(jìn)行分析，探討38個滇山茶品種12個花器官的表型變異：1）確定3個類群在花器官數(shù)量、大小和形狀3個不同維度上花器官變化的方向和程度；2）比較不同分析方法對生成變異規(guī)律的影響。從而分析滇山茶品種的形態(tài)學(xué)多樣性，以獲得滇山茶品種花表型性狀的現(xiàn)狀及變異的規(guī)律，為后期滇山茶定向育種提供理論依據(jù)及分子生物學(xué)等方面提供一定的研究基礎(chǔ)。

1 材料與方法

1.1 試驗材料

本試驗以昆明市植物園的山茶園（102°75′E，25°14′N，海拔1 992 m）和昆明市宜良縣正揚茶花園林綠化有限公司（103°23′E，24°88′N，海拔1 950 m）的健壯滇山茶品種為研究對象，選取生長狀況良好、健康無病害且處于盛花期（中期）的38個滇山茶品種的花（表1）作為試驗材料。

1.2 研究方法

本研究綜合參照中國茶花品種分類、測試指南及已知品種數(shù)據(jù)庫[24]，以及茶樹種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn)[25]。本研究對每個滇山茶品種隨機(jī)選擇4株植株，從每株茶花的中部以上分支采集4朵典型、標(biāo)準(zhǔn)一致的盛花期花，共計16朵，取其平均值進(jìn)行統(tǒng)計分析。然后，立即觀察并記錄所采每朵茶花的12個花部形態(tài)學(xué)指標(biāo)，包括3個數(shù)量特征（X）、7個尺寸特征（C）和2個形狀特征（S）（圖1、表2）。

研究方法包括三個組成部分：一是MS（測量調(diào)查），即對研究對象滇山茶花的各個器官的個體測量的收集，主要是數(shù)量性數(shù)據(jù)的收集；二是VG（視覺觀察調(diào)查），即對研究對象滇山茶花的各個器官的整體視覺觀察，主要是捕捉定性或視覺特征；三是結(jié)合了VS（具體觀察調(diào)查），即對研究對象滇山茶花的各個器官的個體觀察。

1.3 數(shù)據(jù)分析

箱線圖的繪制采用以下方式：箱體覆蓋了觀察區(qū)域的50%，上下橫線之間的區(qū)域覆蓋了觀察區(qū)域的90%，箱外的數(shù)值被視為異常值，本研究的重點在于箱體和均值。平均值用于鄧肯多重比較檢驗確定3大類之間的差異顯著性（P<0.05和P<0.01分別被視為顯著和極顯著），異常值則用于反映超越信息[31]，這對于滇山茶種質(zhì)創(chuàng)新非常重要。頻率分布函數(shù)分析使用Excel 2019、Origin 2022和SPSS 22.0軟件，擬合12個花器官表型性狀的頻率分布函數(shù)，探討變異規(guī)律（包括方向和程度），12個花部表型性狀均服從正態(tài)函數(shù)分布。

2 結(jié)果與分析

2.1 滇山茶不同品種花表型性狀間的多樣性分析

采用F檢驗分析滇山茶不同品種間的差異（表3），12個花表型性狀在38個滇山茶品種間都存在極顯著差異。12個花表型性狀的多樣性指數(shù)范圍為0.40～2.20。其中，雄蕊數(shù)的多樣性指數(shù)最高為2.20，其次為花瓣數(shù)和花徑，多樣性指數(shù)分別為2.18和2.02。而雌蕊長度的多樣性指數(shù)最低為0.40；其次是萼片長度和花瓣長寬比，多樣性指數(shù)分別為0.51和0.56。12個花表型性狀的變異系數(shù)范圍為13.00%～67.00%，平均變異系數(shù)為27.24%，說明38個滇山茶品種在花表型性狀間存在較大程度的變異。涉及12個花表型性狀變異系數(shù)高低依次為雌蕊數(shù)（67.00%）>雄蕊數(shù)（42.60%）>雌蕊長度（42.00%）>萼片長度（27.10%）>萼片長寬比（23.70%）>花瓣寬度（20.90%）>花瓣數(shù)（20.10%）>萼片寬度（19.90%）>雌蕊長度（19.80%）>花瓣長度（15.70%）>花徑（15.10%）>花瓣長寬比（13.00%）。說明各性狀中花瓣長寬比（13.00%）相對比較穩(wěn)定，而雌蕊數(shù)（67.00%）變異程度最大。

變異系數(shù)法的權(quán)重計算結(jié)果顯示（圖2），涉及12個花表型性狀變異系數(shù)權(quán)重高低依次為雌蕊數(shù)（20.485%）>雄蕊數(shù)（13.022%）>雌蕊長度（12.851%）>萼片長度（8.296%）>萼片長寬比（7.254%）>花瓣寬度（6.390%）>花瓣數(shù)（6.148%）>萼片寬度（6.100%）>雄蕊長度（6.056%）>花瓣長度（4.802%）>花徑（4.606%）>花瓣長寬比（3.990%），其中指標(biāo)權(quán)重最大值為雌蕊數(shù)（20.485%），最小值為花瓣長寬比（3.990%）?？傮w而言，雌蕊數(shù)、雄蕊數(shù)和雌蕊長度的變異性在涉及12個花表型性狀影響較大，而其余性狀的影響最小，但仍然對滇山茶的花表型和生殖特征有一定的影響。

2.2 滇山茶不同品種花表型性狀間相關(guān)性分析

38個滇山茶品種的12個花表型性狀相關(guān)性分析結(jié)果（圖2）表明，12個花表型性狀間都存在極顯著相關(guān)關(guān)系（P<0.05），其中雌蕊數(shù)與雄蕊數(shù)、雄蕊長度、雌蕊長度存在極顯著負(fù)相關(guān)，與花瓣長寬比存在顯著正相關(guān)；花瓣數(shù)除了與雌蕊數(shù)和雄蕊數(shù)無顯著性相關(guān)性外，與其他性狀都存在顯著性正相關(guān)或負(fù)相關(guān)；雄蕊數(shù)與雌蕊長度和雄蕊長度存在極顯著正相關(guān)；花徑和花瓣長度除了與雌蕊數(shù)、雄蕊數(shù)、花瓣長寬比和萼片長寬比無顯著性相關(guān)性外，與其他性狀都存在極顯著正相關(guān)；花瓣寬度花徑除了與雌蕊數(shù)、雄蕊數(shù)和萼片長寬比無顯著性相關(guān)性外與其他性狀都存在極顯著正相關(guān)；雄蕊長度與雌蕊長度、萼片長度、萼片寬度存在極顯著正相關(guān)；雌蕊長度與萼片長度、萼片寬度、萼片長寬比存在極顯著正相關(guān)；萼片長度與萼片寬度、萼片長寬比存在極顯著正相關(guān)；萼片長度與花瓣長寬比和萼片長寬比存在顯著負(fù)相關(guān)。這些結(jié)果對于了解滇山茶的花部特征及其遺傳基礎(chǔ)具有一定的重要性。

2.3 滇山茶不同品種花表型性狀主成分分析

對38個滇山茶品種的12個花表型性狀進(jìn)行主成分分析，結(jié)果如表4所示，前4個主成分提取到的特征值均大于1.00，累計貢獻(xiàn)率達(dá)到82.129%，表明這4個主成分能夠反映12個表型性狀的大部分信息。其中，第1主成分的特征值最大，為4.992，對應(yīng)的特征向量中花瓣長度的絕對值最大，為0.854，貢獻(xiàn)率達(dá)到41.602%。第2主成分的特征值為1.871，對應(yīng)的特征向量中雄蕊數(shù)的絕對值最大，為0.657，累計貢獻(xiàn)率達(dá)15.595%。第3主成分的特征值為1.724，對應(yīng)的特征向量中萼片長寬比的絕對值最大，為0.595，貢獻(xiàn)率達(dá)14.365%。第4主成分的特征值最小，為1.268，對應(yīng)的特征向量中花瓣長寬比的絕對值最大，為0.839，貢獻(xiàn)率達(dá)10.567%。這些結(jié)果揭示了前4個主成分在總體上捕捉了大部分表型性狀的變異情況。

2.4 滇山茶不同品種花表型性狀聚類分析

對38個滇山茶品種的12個花表型性狀進(jìn)行聚類分析，在歐氏距離（Euclidean distance）為15.0處將38個滇山茶品種資源分為3大類（圖3）。

本研究中3大類的表型性狀多樣性指數(shù)（H′）變化范圍為0.90～1.10，平均變異系數(shù)分別為27.24%、36.68%和21.52%，說明這3大類的花表型具有豐富的遺傳變異。

第Ⅰ類包含‘楚雄茶’、‘靖安茶’、‘大理茶’（圖3）、‘西林’、‘雪嬌’、‘節(jié)節(jié)高’等21個滇山茶品種，與第Ⅱ類和第Ⅲ類相比（表5），該類的雌蕊數(shù)、花徑、花瓣長度、花瓣寬度、雄蕊長度、萼片長度、萼片寬度、花瓣長寬比和萼片長寬比的平均值最大，僅雄蕊數(shù)、雌蕊長度和花瓣數(shù)處于中等水平。此類群平均變異系數(shù)處于中等水平，為27.24%，其中，花瓣數(shù)的變異系數(shù)最大，為72.50%，花瓣長寬比的變異系數(shù)最小為11.90%。說明該類群花較大，且瓣化程度相對較低，多為半重瓣，在觀賞上具有一定優(yōu)勢。

第Ⅱ類包含‘恨天高’、‘童子面’、‘朱砂紫袍’（圖3）、‘菊瓣’、‘賽菊瓣’、‘紫袍’等11個滇山茶品種，與第Ⅰ類和第Ⅲ類相比，該類僅花瓣數(shù)均值最大，雌蕊數(shù)、花徑、萼片長度和萼片寬度均值處于中等水平，而雄蕊數(shù)、花瓣長度、花瓣寬度、雌蕊長度、雄蕊長度和萼片長寬比的均值相對較小。此類平均變異系數(shù)最大，為36.68%，其中雌蕊數(shù)的變異系數(shù)最大，為 72.50%，其次是雌蕊長度為89.00%，花瓣長寬比的變異系數(shù)最小，為11.90%。說明該類瓣化程度相對較高，多為重瓣花，在觀賞上極具優(yōu)勢。

第Ⅲ類包含‘粉蝴蝶’、‘張家茶’、‘銀粉牡丹’、‘滇池粉妍’、‘東林’和‘早桃紅’（圖3）共6個滇山茶品種，與第Ⅰ類和第Ⅱ類相比，該類僅雄蕊數(shù)和雌蕊長度均值最大，花瓣長度、花瓣寬度和雄蕊長度的均值處于中等水平，花瓣數(shù)、雌蕊數(shù)、花徑、萼片長度、萼片寬度和花瓣長寬比的均值最小。此類平均變異系數(shù)最小，為21.52%，說明該類相比于其他2個類群，穩(wěn)定性較好，其中雌蕊數(shù)的變異系數(shù)最大，為70.00%，雄蕊數(shù)的變異系數(shù)最小，為11.90%。說明該類瓣化程度相對較低，雄蕊保留較多，在結(jié)實和雜交育種上更具優(yōu)勢。

2.5 滇山茶不同品種花表型性狀的綜合評價

滇山茶品種的12個花表型性狀的得分系數(shù)如表6所示，構(gòu)建其4個主成分為參數(shù)的線性組合模型，并分別計算38個滇山茶品種在4個主成分中的綜合系數(shù)（表7）?；诟髦鞒煞值呢暙I(xiàn)率進(jìn)一步建立滇山茶不同品種的綜合評價模型，用于計算38個滇山茶品種綜合得分，模型為F=（41.602F1+15.595F2+14.365F3+10.567F4）/ 82.129，根據(jù)此模型對38個滇山茶品種進(jìn)行綜合評價。在38個滇山茶品種的綜合評價中，排在前6的品種分別為‘大理茶’（3.230）、‘靖安茶’（2.980）、‘咪依嚕’（1.607）、‘鳳山茶’（1.471）、‘柳葉銀紅’（1.350）和‘寶珠茶’（1.292），都為第Ⅰ類的滇山茶品種。整體評價結(jié)果顯示，類群Ⅰ的品種綜合評價較好，排名比較靠前，而第Ⅱ類和第Ⅲ類的排名波動較大，結(jié)合聚類分析結(jié)果，發(fā)現(xiàn)花徑大、雄蕊數(shù)多等特征的滇山茶排名總體靠前，而花徑小、花瓣少的滇山茶品種基本都排在末位，這與滇山茶以花大色艷而著名相吻合。

2.6 滇山茶不同品種花表型性狀箱線圖分析

38個滇山茶品種的12個花表型特征包括器官數(shù)量、大小和形狀，通過箱線圖呈現(xiàn)（圖4）。對這些品種進(jìn)行了聚類分析分成3大類，差異分析結(jié)果表明，有9個花表型性狀具有顯著差異，其余3個花表型性狀則沒有明顯差異。這些發(fā)現(xiàn)可為滇山茶品種間花表型差異提供有力的統(tǒng)計支持，特別是在花瓣數(shù)、雄蕊數(shù)等方面的顯著變異，有助于深入挖掘品種間的形態(tài)差異和遺傳多樣性。

對38個滇山茶品種聚類的3大類群進(jìn)行分布比較，相比于第Ⅰ類，除花瓣數(shù)、雌蕊數(shù)、雄蕊數(shù)、花瓣寬度、萼片寬度外，第Ⅱ類和第Ⅲ類的其他花表型性狀均呈向下分布。各表型性狀的具體分布如下：在花器官數(shù)方面（圖5），第Ⅱ類的花瓣數(shù)及第Ⅰ類的雌蕊數(shù)、雄蕊數(shù)分布范圍較廣，第Ⅰ類的花瓣數(shù)和雌蕊數(shù)存在一定的越界個體（比例分別為2.30%和1.20%），3大類的雄蕊數(shù)分布比較集中，但存在極顯著差異，說明該性狀變異性較大；雌蕊數(shù)分布范圍較廣，但無顯著差異，說明該性狀變異性較小。在花器官大小方面，第Ⅰ類的花瓣寬度、花瓣長度、雄蕊長度、雌蕊長度、萼片長度、萼片寬度和類群Ⅲ的花直徑分布范圍較大。其中，第Ⅰ類的花瓣寬度、雄蕊長度、萼片寬度和第Ⅱ類的雌蕊長度、萼片長度、萼片寬度存在一定的越界個體（比例分別為5.90%、1.20%、 1.20%、9.00%、6.80%和2.30%），3個類群的雌蕊長度分布較集中，且存在顯著差異，說明該性狀變異性較大；花瓣寬度分布范圍較廣，但無顯著差異，說明該性狀變異性較小。從花器官形狀來看，第Ⅰ類的花瓣長寬比和第Ⅱ類的萼片長寬比的分布范圍都較大，且存在一定的越界個體（比例分別為3.60%和4.50%）。

2.7 滇山茶不同品種花表型變異的頻率分布函數(shù)分析

為了進(jìn)一步分析滇山茶3大類花表型性狀的變異情況，對12個花表型性狀進(jìn)行了頻率分布函數(shù)擬合（圖6），所有花器官表型性狀均遵循正態(tài)分布（R2 =0.916 4～0.998 3）（表8）。對38個滇山茶品種的3大類進(jìn)行頻率分布函數(shù)比較，除雌蕊數(shù)、花瓣長度、花瓣寬度的冪分布函數(shù)比較集中之外，其余9個性狀的冪函數(shù)分布較離散，進(jìn)一步說明滇山茶不同品種花表型性狀變異性較大。其中，在數(shù)量性狀方面，3大類花瓣數(shù)和雌蕊數(shù)的正態(tài)分布函數(shù)呈右偏（增加趨勢），而雄蕊數(shù)的正態(tài)分布函數(shù)呈左偏（減少趨勢），這表明增加的花瓣數(shù)可能由雄蕊轉(zhuǎn)化產(chǎn)生；在尺寸特征方面，花徑、雄蕊長度、雌蕊長度和萼片長度顯示出縱向和橫向偏移的趨勢，說明其性狀變異性較大，而花瓣長度和花瓣寬度的冪分布函數(shù)沒有出現(xiàn)明顯的偏移，說明其性狀變異性較?。辉谛螤钐卣鞣矫?，第Ⅰ類的花瓣長寬比相比于第Ⅱ類和第Ⅲ類的正態(tài)分布函數(shù)呈右偏，說明該類群相比于其他2個類變異性較大，而萼片長寬比在橫向方向沒有出現(xiàn)明顯的偏移，說明該性狀變異性較小。

3 討 論

表型性狀對植物種類的鑒別、授粉生態(tài)學(xué)的研究以及植物進(jìn)化和生態(tài)方面的理解具有重要意義[32]。通過對大規(guī)模表型數(shù)據(jù)的分析，可以更全面地了解基因型和表型之間的關(guān)系以及各種中間環(huán)節(jié)對性狀表達(dá)的影響[33-34]，有助于揭示物種內(nèi)遺傳多樣性和遺傳規(guī)律，為進(jìn)化、遺傳學(xué)和生物多樣性研究提供重要的數(shù)據(jù)支持[35]。利用相關(guān)性、聚類以及主成分分析方法對種質(zhì)資源性狀遺傳分析研究[36]，廣泛應(yīng)用于育種、多樣性研究及種質(zhì)分類[37]。花是植物重要的組織器官，在生命周期中扮演著多重角色，既為觀賞帶來美感，又在植物繁殖中發(fā)揮著重要的生物學(xué)功能[17]。通過對滇山茶花表型特征的研究和分析，有助于篩選出更具優(yōu)良性狀的滇山茶品種，對花表型性狀的改良和選育過程中的優(yōu)化具有重要作用[12]。

本研究中，38個滇山茶品種的12個花表型性狀的多樣性指數(shù)在0.40～2.20，其中有9個花表型性狀的多樣性指數(shù)大于1.00，說明滇山茶品種具有較為豐富的花表型多樣性，這與吳貴進(jìn)等[12]的研究結(jié)果一致。變異系數(shù)是一種用于度量數(shù)據(jù)離散程度的統(tǒng)計指標(biāo)，反映了性狀變量的離散程度[38]。本研究中，12個花表型性狀的變異系數(shù)范圍為13.00%～67.00%，均大于10.00%，證實了滇山茶不同品種間花表型性狀的變異程度較大，這與前人的研究結(jié)果相吻合[39]，這也說明滇山茶品種資源多樣性豐富，與辛桐[40]的研究保持一致，其中雌蕊數(shù)的變異系數(shù)（67.00%）最高，這可能與選取的材料中有些品種雌蕊數(shù)已退化有關(guān)。通過比較分析多樣性指數(shù)和變異系數(shù)發(fā)現(xiàn)，變異系數(shù)與多樣性指數(shù)的變化趨勢相反，滇山茶變異系數(shù)大的性狀其多樣性指數(shù)相對較低，這與張海平等[41]對睡蓮、丁釋豐等[42]對黃心樹及蘇群等[43]對睡蓮的研究結(jié)論相同。

相關(guān)性分析對滇山茶品種改良和生物技術(shù)研究都具有重要的指導(dǎo)意義[44]。本研究中，12個花表型性狀間都存在極顯著的相關(guān)關(guān)系，說明這12個花表性狀可能在滇山茶品種的花朵的演化過程中發(fā)揮著密切的協(xié)同作用，如雌蕊和雄蕊存在瓣化，從而導(dǎo)致花瓣數(shù)的增加。聚類分析結(jié)果表明，第I類可作為觀賞或育種優(yōu)質(zhì)親本材料；第II類可用于培育新的品種，以改進(jìn)觀賞性和適應(yīng)性等特征；第Ⅲ類滇山茶品種在某些特定用途或研究中更可靠和可預(yù)測。植物種質(zhì)資源綜合評價在育種研究中具有重要意義[45]，綜合評價結(jié)果表明，‘大理茶’的綜合得分最高（F=3.230），說明該品種在38個滇山茶品種綜合特性最好，結(jié)合12個花部性狀，該品種相較于其他品種觀賞價值高，但瓣化現(xiàn)象嚴(yán)重，不利于繁育后代。

在箱線分析中，在各類群中的花器官數(shù)量、大小和形狀方面存在一定比例的越界個體。比如，在花器官數(shù)量和大小方面，某些滇山茶品種可能表現(xiàn)出與其所屬類群不一致的特征，這可能是由于突變、特定生態(tài)條件、遺傳重組或其他遺傳機(jī)制引起的，可為選育新的滇山茶品種提供一定的參考。本研究在三大類進(jìn)行應(yīng)用頻率分布函數(shù)分析中，清楚地揭示滇山茶的表型變異，其所有花器官表型性狀均遵循正態(tài)分布，擬合度較高（R2為0.916 4～0.998 3），這也與箱線分析的結(jié)果保持一致，這些結(jié)果為滇山茶品種特性的理解和進(jìn)一步的育種提供了重要信息，為后續(xù)的滇山茶相關(guān)形態(tài)以及機(jī)理研究提供了方向和思路。

本研究通過對38個滇山茶品種的12個花表型性狀的全面分析，展現(xiàn)了滇山茶在花的形態(tài)特征上的豐富多樣性和顯著的變異性，進(jìn)而為滇山茶品種的選育和生物多樣性保護(hù)提供了有力的科學(xué)依據(jù)。研究揭示了38個滇山茶12個花表型性狀之間的復(fù)雜關(guān)系和演化趨勢，對于理解其遺傳多樣性和生物學(xué)特性具有重要意義。盡管如此，本研究存在局限性，主要體現(xiàn)在樣本數(shù)量有限，未能完全捕捉到滇山茶的整體遺傳和表型多樣性；研究方法側(cè)重于傳統(tǒng)的統(tǒng)計分析，對于表型性狀與基因型之間深層次聯(lián)系的解析不足；同時，對于環(huán)境因素如生長地點、氣候條件等外部影響因素考慮不夠全面，這可能影響到對表型性狀表達(dá)影響的全面理解。

鑒于這些局限性，未來的研究方向應(yīng)當(dāng)包括：1）通過擴(kuò)大滇山茶品種的樣本量和樣本覆蓋范圍，增加不同生態(tài)環(huán)境下的樣本，以期獲得更加全面和具有代表性的數(shù)據(jù)集；2）結(jié)合現(xiàn)代生物技術(shù)如分子信息標(biāo)記、基因組學(xué)和轉(zhuǎn)錄組學(xué)等方法，探討表型性狀與遺傳信息之間的聯(lián)系，以揭示性狀變異的分子機(jī)理，特別是滇山茶花瓣化現(xiàn)象的研究機(jī)理，如在油茶[46-47]的應(yīng)用；3）考慮到環(huán)境因素對植物表型的影響，在今后的研究中應(yīng)系統(tǒng)評估環(huán)境變量對滇山茶表型性狀表達(dá)的影響，以及可能的遺傳與環(huán)境的互作效應(yīng)；4）采用先進(jìn)數(shù)據(jù)分析技術(shù)提升分析效率和精確度，從大規(guī)模數(shù)據(jù)中發(fā)掘更深層次的生物學(xué)信息；5）開展功能性研究，如研究特定表型性狀（瓣化現(xiàn)象）對滇山茶觀賞價值和生態(tài)適應(yīng)性的影響，為滇山茶的應(yīng)用和保護(hù)提供科學(xué)依據(jù)。通過這些綜合性和多角度的研究，進(jìn)一步深化對滇山茶遺傳和表型多樣性的認(rèn)識，為其種質(zhì)資源的保護(hù)、利用和持續(xù)育種提供更加堅實的科學(xué)基礎(chǔ)。

4 結(jié) 論

本研究對38個滇山茶品種的12個花表型性狀進(jìn)行了多樣性分析，結(jié)果表明38個滇山茶品種在12個花的表型性狀上展現(xiàn)出豐富的多樣性和較大的變異程度。證實了滇山茶作為一個種質(zhì)資源的多樣性和豐富性。綜合評價結(jié)果表明，滇山茶品種‘大理茶’在綜合特性上表現(xiàn)最佳，但其瓣化現(xiàn)象的存在也提示了對繁育后代可能的不利影響。箱線分析和應(yīng)用頻率分布函數(shù)分析結(jié)果進(jìn)一步揭示了滇山茶的表型變異規(guī)律。特別值得注意的是，花表型性狀間的顯著相關(guān)性揭示了它們在滇山茶花朵演化中可能存在的密切協(xié)同作用，主要以雌雄蕊群逐漸瓣化為演化方向。研究結(jié)果確定了滇山茶花表型演化的趨勢，為理解植物演化過程中的功能性狀協(xié)同變化提供了新的視角，為深入了解滇山茶品種特性及其育種提供了寶貴信息。

參考文獻(xiàn)：

[1] FU M， YANG X， ZHENG J， et al. Unraveling the regulatory mechanism of color diversity in Camellia japonica petals by integrative transcriptome and metabolome analysis[J]. Frontiers in Plant Science，2021，12：685136.

[2] GENG F， NIE R， YANG N， et al. Integrated transcriptome and metabolome profiling of Camellia reticulata reveal mechanisms of flower color differentiation[J]. Frontiers in Genetics，2022，13：1059717.

[3] 夏麗芳.山茶花[M].北京：中國建筑工業(yè)出版社，2002. XIA L F. Camellia flower[M]. Beijing： China Architecture and Building Press，2002.

[4] 陳世軍，張明澤，姚玉仙.基于SSR標(biāo)記的黔南茶樹種質(zhì)資源DNA指紋圖譜構(gòu)建[J].植物遺傳資源學(xué)報，2017，18（1）： 106-111. CHEN S J， ZHANG M Z， YAO Y X. Establishment of DNA fingerprinting for tea germplasm from Qiannan Prefecture by SSR markers[J]. Journal of Plant Genetic Resources，2017， 18（1）：106-111.

[5] VELA P， SALINERO C， SAINZ M J. Phenological growth stages of Camellia japonica[J]. Annals of Applied Biology，2013， 162（2）：182-190.

[6] 溫佳辛，王超林，馮慧，等.月季花色研究進(jìn)展[J].園藝學(xué)報， 2021，48（10）：2044-2056. WEN J X， WANG C L， FENG H， et al. Research progress on flower color of rose[J]. Acta Horticulturae Sinica，2021，48（10）：2044-2056.

[7] 尹世華，王康，黃曉霞，等.47份月季品種表型多樣性分析及綜合評價[J].江西農(nóng)業(yè)大學(xué)學(xué)報，2021，43（1）：94-105. YIN S H， WANG K， HUANG X X， et al. Phenotypic diversity analysis and comprehensive evaluation of 47 rose resources[J]. Acta Agriculturae Universitatis Jiangxiensis，2021，43（1）：94-105.

[8] 王禹，張廣輝，赫京生，等.杜鵑花色研究進(jìn)展[J].世界林業(yè)研究，2020，33（5）：19-24. WANG Y， ZHANG G H， HE J S， et al. Research progress of Rhododendron flower color[J]. World Forestry Research，2020， 33（5）：19-24.

[9] 沈德周，縱丹，周安佩，等.AFLP引物組合數(shù)對準(zhǔn)確分析馬纓杜鵑遺傳多樣性的影響[J].西南林業(yè)大學(xué)學(xué)報，2016，36（6）： 22-28. SHEN D Z， ZONG D， ZHOU A P， et al. Effect of the quantity of AFLP primer combinations on accurately analyzing Rhododendron delavayi genetic diversity[J]. Journal of Southwest Forestry University，2016，36（6）：22-28.

[10] 司馬永康，郭玉紅，徐德兵，等.云南紅花油茶基于花瓣數(shù)量的分類[J].西部林業(yè)科學(xué)，2018，47（5）：80-89. SIMA Y K， GUO Y H， XU D B， et al. Classification of Camellia reticulata Lindley based on its petal numbers[J]. Journal of West China Forestry Science，2018，47（5）：80-89.

[11] 沈立新，梁洛輝，王慶華，等.騰沖紅花油茶自然類型及其品種類群劃分[J].林業(yè)資源管理，2009（6）：75-80. SHEN L X， LIANG L H， WANG Q H， et al. The natural types and variety groups of Camellia reticulata. Lindley[J]. Forest Resources Management，2009（6）：75-80.

[12] 吳貴進(jìn)，陳龍清，耿芳，等.滇山茶登錄品種統(tǒng)計及表型性狀分析[J].園藝學(xué)報，2023，50（10）：2157-2170. WU G J， CHEN L Q， GENG F， et al. Statistics and phenotypic traits analysis of Camellia reticulata registered cultivars[J]. Acta Horticulturae Sinica，2023，50（10）：2157-2170.

[13] 林存學(xué)，楊曉華，劉海榮.東北寒地96份李種質(zhì)資源表型性狀遺傳多樣性分析[J].園藝學(xué)報，2020，47（10）：1917-1929. LIN C X， YANG X H， LIU H R. Genetic diversity analysis of 96 plum germplasm resources by phenotypic traits in northeast cold area[J]. Acta Horticulturae Sinica，2020，47（10）：1917-1929.

[14] WESSINGE C A， HILEMAN L C. Accessibility， constraint， and repetition in adaptive floral evolution[J]. Developmental Biology，2016，419（1）：175-183.

[15] KHADIVI-KHUBAND A， ANJAM K. Morphological characterization of Prunus scoparia using multivariate analysis[J]. Plant Systematics and Evolution，2014，300：1361-1372.

[16] OLIVEIRA E J D， DIAS N L P， DANTAS J L L. Selection of morpho-agronomic descriptors for characterization of papaya cultivars[J]. Euphytica，2012，185（2）：253-265.

[17] 唐毓瑋，彭綠春，段福軍，等.15個高山杜鵑品種的花表型性狀多樣性分析[J/OL].分子植物育種，https：//link.cnki.net/urlid/ 46.1068.S.20231025.1128.004. TANG Y W， PENG L C， DUAN F J， et al. Analysis of floral phenotypic trait diversity of 15 Rhododendrons cultivars[J/OL]. Molecular Plant Breeding，https：//link.cnki.net/urlid/46.1068. S.20231025.1128.004.

[18] 徐斌，彭莉霞，張方秋，等.杜鵑紅山茶葉片主要性狀的遺傳多樣性分析[J].植物研究，2015，35（5）：730-734. XU B， PENG L X， ZHANG F Q， et al. Genetic diversity analysis for leaf main traits of Camellia azalea[J]. Bulletin of Botanical Research，2015，35（5）：730-734.

[19] 伍思宇，梁曉靜，廖健明，等.中國金花茶與越南金花茶葉表型差異性分析[J].廣西林業(yè)科學(xué)，2019，48（3）：336-341. WU S Y， LIANG X J， LIAO J M， et al. Analysis on leaf phenotypic differences of Chinese and Vietnamese Camellia section Chrysantha Chang[J]. Guangxi Forestry Science，2019， 48（3）：336-341.

[20] 戴艷花，杜佩蓮，郝丙青，等.香花油茶不同家系群體葉片和果實表型性狀的比較[J].經(jīng)濟(jì)林研究，2023，41（4）：109-117. DAI Y H， DU P L， HAO B Q， et al. Comparison of phenotypic characters of leave and fruit in different family groups of Camellia osmantha[J]. Non-wood Forest Research，2023，41（4）：109-117.

[21] 李芳，霍達(dá)，羅揚，等.西南紅山茶天然居群種實表型性狀變異分析[J].種子，2021，40（4）：109-114. LI F， HUO D， LUO Y， et al. Variation analysis of phenotypic traits of seed of natural populations of Camellia pitardii in southwest China[J]. Seed，2021，40（4）：109-114.

[22] 甘沛華，沈德周，何承忠，等.不同花型滇山茶遺傳多樣性與遺傳關(guān)系的AFLP分析[J].植物生理學(xué)報，2019，55（8）：1222-1230. GAN P H， SHEN D Z， HE C Z， et al. AFLP analysis on genetic diversity and relationship of Camellia reticulata of different floral patterns[J]. Plant Physiology Journal，2019，55（8）：1222-1230.

[23] 韋開萍，楊茜，和潤喜，等.騰沖紅花油茶花芽分化及雌配子體發(fā)育的形態(tài)學(xué)觀察[J].經(jīng)濟(jì)林研究，2023，41（3）：244-251. WEI K P， YANG X， HE R X， et al. Morphological study on floral bud differentiation and the development of female gametophytes in Camellia reticulate Lindl f. simplex Sealy[J]. Non-wood Forest Research，2023，41（3）：244-251.

[ 24] 張曉慶.中國茶花品種分類、測試指南及已知品種數(shù)據(jù)庫構(gòu)建[D].北京：中國林業(yè)科學(xué)研究院，2008. ZHANG X Q. The classification of Chinese Camellia cultivars， test guideline and construction of database for the known cultivars[D]. Beijing： Chinese Academy of Forestry，2008.

[25] 中華人民共和國農(nóng)業(yè)部.農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn)NY/T 2943—2016 茶樹種質(zhì)資源描述規(guī)范[S].北京：中國農(nóng)業(yè)出版社，2016. Ministry of Agriculture of the PRC. Chinese Agricultural Standard NY/T 2943-2016 Specification for the description of tea plant germplasm resources[S]. Beijing： China Agriculture Press，2016.

[26] SUMARTINI， T. Evaluation of mungbean lines for resistance to Cercospora leaves spot and powdery mildew （Erysiphe polygoni） diseases in the green house[J]. Earth and Environmental Science，2021，743（1）：1-7.

[27] 張悅，索玉靜，李芳東，等.柿種質(zhì)資源果實形態(tài)多樣性分析[J].園藝學(xué)報，2022，49（7）：1473-1490. ZHANG Y， SUO Y J， LI F D， et al. Analysis on fruit morphological diversity of persimmon germplasm resources[J]. Acta Horticulturae Sinica，2022，49（7）：1473-1490.

[28] 喻華平，趙志常，高愛平，等.基于主成分分析和聚類分析的23份黃皮種質(zhì)資源的品質(zhì)評價[J].熱帶作物學(xué)報，2022，43（7）： 1357-1364. YU H P， ZHAO Z C， GAO A P， et al. Quality evaluation of 23 species of Clausena lansium （Lour.） Skeels germplasm resources based on principal component analysis and cluster analysis[J]. Chinese Journal of Tropical Crops，2022，43（7）：1357-1364.

[29] 郭松，李在留，陳晶晶，等.不同種源掌葉木果實和種子表型性狀多樣性分析及綜合評價[J].植物資源與環(huán)境學(xué)報，2018， 27（4）：11-20. GUO S， LI Z L， CHEN J J， et al. Diversity analysis and comprehensive evaluation on phenotypic traits of fruit and seed of Handeliodendron bodinieri from different provenances[J]. Journal of Plant Resources and Environment，2018，27（4）：11-20.

[30] SHAIJU P N， OMANAKUMARI N. Floral morphology and systematics of the genus Thottea Rottb. （Aristolochiaceae） from the western Ghats， India[J]. Plant Systematics and Evolution，2012，288：213-225.

[31] ZHOU T， FAN J J， ZHAO M M， et al. Phenotypic variation of floral organs in Malus using frequency distribution functions[J]. BMC Plant Biology，2019，19（1）：1-11.

[32] 向貴生，王其剛，邱顯欽，等.云南川滇薔薇天然居群表型多樣性分析[J].云南大學(xué)報，2018，40（4）：786-794. XIANG G S， WANG Q G， QIU X Q， et al. Phenotypic diversity of natural population of Rosa soulieana in Yunnan[J]. Journal of Yunnan University，2018，40（4）：786-794.

[33] 張瑩，曹玉芬，張小雙，等.基于花表型性狀的梨種質(zhì)資源多樣性研究[J].園藝學(xué)報，2016，43（7）：1245-1256. ZHANG Y， CAO Y F， ZHANG X S， et al. Research on diversity of pear germplasm resources based on flowers phenotype traits[J]. Acta Horticulturae Sinica，2016，43（7）：1245-1256.

[34] 頡剛剛，歐陽麗婷，謝軍，等.新疆地區(qū)歐洲李葉片表型性狀多樣性及親緣關(guān)系分析[J].植物資源與環(huán)境學(xué)報，2018，27（3）： 72-78. XIE G G， OUYANG L T， XIE J， et al. Analyses on diversity of leaf phenotypic traits and genetic relationships of Prunus domestica in Xinjiang region[J]. Journal of Plant Resources and Environment，2018，27（3）：72-78.

[35] 吳根松，孫麗丹，郝瑞杰，等.梅花種質(zhì)資源表型多樣性研究[J].安徽農(nóng)業(yè)科學(xué)，2011，39（20）：12008-12009. WU G S， SUN L D， HAO R J， et al. Study on the phenotypic diversity of P. mume Sieb. et Zucc. germplasm resources[J]. Journal of Anhui Agricultural Sciences，2011，39（20）：12008-12009.

[36] 陳文星，楊琳懿，唐軍榮，等. 10個新樟屬植物葉片表型性狀多樣性分析[J].中南林業(yè)科技大學(xué)學(xué)報，2023，43（12）：83-93. CHEN W X， YANG L Y， TANG J R， et al. Analysis of the diversity on leaf phenotypic traits of 10 Neocinnamomum genus[J]. Journal of Central South University of Forestry Technology，2023，43（12）：83-93.

[37] 孫澤碩，李永華，柳新紅，等.四十一份櫻花品種資源表型性狀遺傳多樣性分析[J].北方園藝，2022（22）：57-66. SUN Z S， LI Y H， LIU X H， et al. Genetic diversity analysis of phenotypic traits in 41 cultivars of oriental cherry[J]. Northern Horticulture，2022（22）：57-66.

[38] 柳江群，尹明宇，左絲雨，等.長柄扁桃天然種群表型變異[J].植物生態(tài)學(xué)報，2017，41（10）：1091-1102. LIU J Q， YIN M Y， ZUO S Y， et al. Phenotypic variations in natural populations of Amygdalus pedunculata[J]. Chinese Journal of Plant Ecology，2017，41（10）：1091-1102.

[39] 王輝，謝利娟，羅丹，等.深圳外來月季品種的數(shù)量分類研究[J].西南林業(yè)大學(xué)學(xué)報，2013，33（4）：81-87. WANG H， XIE L J， LUO D， et al. Mathematic classification of the introduced rose cultivars in Shenzhen[J]. Journal of Southwest Forestry University，2013，33（4）：81-87.

[40] 辛桐.云南山茶花Camellia reticulata多樣性的人文因素及遺傳背景[D].北京：中央民族大學(xué)，2015. XIN T. The cultural factors and genetic background of the diversity of Camellia reticulata[D]. Beijing： Central University for Nationalities，2015.

[41] 張海平，房偉民，陳發(fā)棣，等.部分睡蓮屬植物形態(tài)性狀的多樣性分析[J].南京農(nóng)業(yè)大學(xué)學(xué)報，2009，32（4）：47-52. ZHANG H P， FANG W M， CHEN F D， et al. Investigation on the morphological diversity of taxa in genus Nymphaea[J]. Journal of Nanjing Agricultural University，2009，32（4）：47-52.

[42] 丁釋豐，高筱鈺，吳林源，等.廣東黃心樹種質(zhì)表型與遺傳多樣性分析[J].南方農(nóng)業(yè)學(xué)報，2021，52（5）：1291-1299. DING S F， GAO X Y， WU L Y， et al. Phenotype and genetic diversity of Machilus gamblei King ex Hook. f. germplasm in Guangdong[J]. Journal of Southern Agriculture，2021，52（5）：1291-1299.

[43] 蘇群，楊亞涵，田敏，等.49份睡蓮資源表型多樣性分析及綜合評價[J].西南農(nóng)業(yè)學(xué)報，2019，32（11）：2670-2681. SU Q， YANG Y H， TIAN M， et al. Phenotypic diversity analysis and comprehensive evaluation of 49 waterlily resources[J]. Southwest China Journal of Agricultural Sciences，2019，32（11）：2670-2681.

[44] 陶萍，宋燚，張晴暉，等.云南普洱茶種古茶樹資源表型性狀數(shù)據(jù)分析研究[J].西部林業(yè)科學(xué)，2020，49（6）：144-151. TAO P， SONG Y， ZHANG Q H， et al. Data analysis of phenotypic characteristics of ancient tea resources of Yunnan Pu’er tea species[J]. Journal of West China Forestry Science，2020， 49（6）：144-151.

[45] KIM E J， SA K J， PARK K C. Study of genetic diversity and relationships among accessions of foxtail millet Setaria italica （L.） P. Beauv. in Korea， China， and Pakistan using SSR marker[J]. Genes Genomics，2012，34（5）：529-539.

[46] 童海浪，張華健，李健，等.四倍體油茶全基因組SSR位點開發(fā)與特征分析[J].中南林業(yè)科技大學(xué)學(xué)報，2024，44（3）：126-135. TONG H L， ZHANG H J， LI J， et al. Development and characterization of whole genome SSR in tetraploid Camellia oleifera[J]. Journal of Central South University of Forestry Technology，2024，44（3）：126-135.

[47] 付超，劉凱，張幸，等.廣西鳳山中籽茶群體種質(zhì)SNP標(biāo)記開發(fā)及遺傳進(jìn)化分析[J].經(jīng)濟(jì)林研究，2023，41（3）：169-175，186. FU C， LIU K， ZHANG X， et al. SNP markers development and genetic structure analysis of Camellia oleifera genetic populations germplasm in Fengshan of Guangxi[J]. Non-wood Forest Research，2023，41（3）：169-175，186.

[本文編校：吳 彬]