，， ，，，

(1.西南林業(yè)大學(xué),云南省高校林木遺傳改良與繁育重點(diǎn)實(shí)驗(yàn)室，昆明 650224；2.北京林業(yè)大學(xué),林木育種國家工程實(shí)驗(yàn)室，北京 100083)

云南松(PinusyunnanensisFranch.)是我國云貴高原的主要鄉(xiāng)土樹種[1-3]，具有適應(yīng)性強(qiáng)、耐干旱瘠薄、木材用途廣等特點(diǎn)，在分布區(qū)域林業(yè)生產(chǎn)和生態(tài)經(jīng)濟(jì)建設(shè)中占有重要的地位[3-5]，天然分布于東經(jīng)96°～108°、北緯23°～30°之間，分布范圍廣，生境復(fù)雜多樣[6]，其中云南省是云南松的集中分布區(qū)[3,7-8]。目前云南松現(xiàn)存林分低矮、彎曲、扭曲等不良個體的比例逐漸增加，表現(xiàn)出林分衰退問題日益突出現(xiàn)象，亟需對云南松種質(zhì)資源開展保護(hù)工作，而遺傳多樣性研究是開展種質(zhì)資源保護(hù)的基礎(chǔ)[9]。

表型多樣性是遺傳多樣性與環(huán)境多樣性的綜合體現(xiàn)，主要研究群體在其分布區(qū)內(nèi)各種環(huán)境下的表型變異，是遺傳多樣性和生物系統(tǒng)學(xué)研究的重要內(nèi)容[10-11]，其表型性狀形態(tài)學(xué)結(jié)構(gòu)特征也是植物遺傳多樣性研究形態(tài)學(xué)標(biāo)記的重要內(nèi)容[12-16]，同時通過表型來檢測群體的變異和分化又是一種簡便、快捷、有效的方法[11,17]，特別是在短期內(nèi)對種質(zhì)資源進(jìn)行評價時，常作為一種重要的評價手段。球果作為裸子植物繁殖系統(tǒng)的重要特征[18]，在選擇壓力下會發(fā)生適應(yīng)性演化[19-20]，且果實(shí)形態(tài)往往是較穩(wěn)定的遺傳特征[20]，其表型性狀不僅決定物種的擴(kuò)散能力，也影響到物種的萌發(fā)和幼苗定居，進(jìn)而影響到種群的分布格局[21]。云南松種實(shí)性狀在不同海拔梯度群體或不同莖干類型群體的相關(guān)研究中初步表明存在一定的變異[22-24]，但針對不同分布區(qū)域群體間的球果表型變異未見研究報道。鑒于此，本研究在前期針葉表型變異研究基礎(chǔ)上[25]，以云南松主分布區(qū)不同分布區(qū)域天然群體的球果為材料，分析云南松群體間和群體內(nèi)球果性狀的表型變異程度與變異規(guī)律，以期為云南松資源保護(hù)與合理利用提供依據(jù)。

1 材料與方法

1.1 群體的分布及球果的采集

在云南松自然分布區(qū)內(nèi)選擇天然群體9個(表1)，樣株的年齡均在25年生以上[14]。為了最大化地代表各群體的遺傳變異狀況，采集的樣品均來自于不同植株，且各樣株不相鄰，間隔5倍以上樹高[23]，從樹冠南向中上部采集當(dāng)年成熟的球果，用于表型性狀測定。

表1 云南松天然群體信息[23]

表2 云南松球果性狀的方差分析

注：**表示在0.01%水平上差異顯著。

1.2 球果性狀的測定

包括球果長、球果粗、球果重量、球果長/球果粗等4個性狀。用電子游標(biāo)卡尺測量球果長(兩端的直線長度)和球果粗(最寬處直徑)[15]，精確到0.01 mm，用電子天平測量球果重，精確到0.01 g，計算球果長/球果粗，即球果的形狀指數(shù)，精確到0.01，每單株測定5～10個球果。

1.3 土壤樣品的采集與測定

土壤樣品的采集與測定同前期研究所述[25]，分別測定含水量、pH值、有機(jī)質(zhì)、水解氮、速效磷、速效鉀、全氮、全磷、全鉀及質(zhì)地共10個指標(biāo)，其中質(zhì)地采用定性描述，未用于后續(xù)的分析中。

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

對球果性狀采用巢式方差分析[10,26-27]計算群體間和群體內(nèi)的方差分量，線性模型為：

Yijk=μ+τi+δj(i)+εk(ij)，

2 結(jié)果與分析

2.1 群體間和群體內(nèi)的球果性狀變異特征

將4個球果表型性狀進(jìn)行巢式方差分析，各性狀的均方及F值見表2。云南松4個球果性狀在群體間和群體內(nèi)均呈極顯著差異(p<0.01)，表明云南松的球果表型性狀在群體間和群體內(nèi)存在廣泛的變異[23]。

2.2 球果性狀群體間的形態(tài)分化

對各個性狀成對群體間進(jìn)一步分析形態(tài)分化PST(表3和圖1)。球果重、球果長、球果粗和球果長/球果粗在各成對群體間的形態(tài)分化波動于0.000 0～0.462 9、0.000 0～0.213 1、0.000 0～0.605 7和0.000 0～0.145 5，平均分別為0.124 1、0.051 5、0.162 1和0.026 4，以球果長/球果粗有形態(tài)分化最低，即該性狀最穩(wěn)定?？傮w來看，球果性狀在各群體間的分化均較低，其中僅有一對群體間在球果粗方面表現(xiàn)出分化系數(shù)大于0.5。

2.3 球果性狀的變異特征

對同一群體的不同性狀或同一性狀的不同群體進(jìn)行比較[23]，9個群體球果表型性狀的相對極差見表4。云南松球果重、球果長、球果粗和球果長/球果粗的相對極差波動于0.42～0.96、0.43～0.90、0.46～1.00和0.31～0.96，平均分別為0.66、0.71、0.61和0.62。綜合各群體的相對極差來看，BH和MD群體變異程度比較高，而TC群體的變異比較低。

表3 云南松天然群體球果性狀的形態(tài)分化(PST)

表4 云南松天然群體球果性狀的相對極差

注：群體名稱的縮寫見表1。

注：A:PST=0；B:0

2.4 球果性狀與土壤因子間的相關(guān)性

對云南松球果性狀與土壤因子間相關(guān)分析表明(表5)，除球果重與有機(jī)質(zhì)含量、全氮含量，及球果長/球果粗與全鉀含量之間的相關(guān)性為顯著正相關(guān)關(guān)系外(p<0.05)，其它因子間的相關(guān)性均未達(dá)到顯著水平，相關(guān)性弱，變異可能是隨機(jī)的分布[23]。

表5 云南松天然群體球果性狀與土壤因子間的相關(guān)分析 (Spearman’s)

注：*表示在0.05%水平上相關(guān)性顯著。

2.5 基于球果性狀云南松天然群體的聚類

基于云南松4個球果性狀，對9個群體進(jìn)行主成分分析，結(jié)果見表6和圖2。第一主成分主要解釋球果長、球果粗和球果重，即球果的大小及其重量；第二主成分主要解釋球果長/球果粗，即球果的形狀。其中第一主成分將BH群體明顯區(qū)分于其它群體，第三主成分使MD與其它群體區(qū)分開來[23]。

將4個球果性狀在各群體的均值標(biāo)準(zhǔn)化后，計算各群體間的歐式距離，結(jié)果見表7。各距離波動于0.761～5.255之間，以BH群體和SJ群體間的距離最大，而以HQ群體和LX群體間的距離最小。

表6 云南松天然群體球果性狀的主成分分析

表7 基于球果性狀云南松天然群體間的歐式距離

注：群體名稱的縮寫見表1。

圖2 基于球果性狀云南松天然群體的主成分分析三維圖

圖3 基于球果性狀云南松天然群體的UPGMA聚類圖

基于球果性狀各成對群體間的歐氏距離，繪制云南松9個群體的UPGMA聚類圖(圖3)，結(jié)果表明，當(dāng)以歐氏距離2.8為閾值時，9個群體被分為3大類：其中BH單獨(dú)形成一類(Ⅰ)，SJ和EY聚成另一類(Ⅱ)，其余6個群體形成第三類(Ⅲ)。當(dāng)以歐氏距離1.5為閾值時，可將第三大類再分為3個亞類：LX和HQ 2個群體形成的第1亞類(Ⅲ-1)，KM、TC和MD 3個群體構(gòu)成的第2亞類(Ⅲ-2)，JS獨(dú)自形成一個亞類(Ⅲ-3)[23]。

3 討論與結(jié)論

遺傳變異是群體生物逐漸適應(yīng)環(huán)境和不斷進(jìn)化的產(chǎn)物，是遺傳信息的重要表征[35-36]。云南松球果表型性狀在群體間和群體內(nèi)均存在極顯著差異，在白皮松[10]、川西云杉[37]、云杉[38-40]、油松[41]、紅松[42]、華北落葉松[43]等研究中均揭示表型性狀變異的廣泛性，廣泛的形態(tài)變異能適應(yīng)不同環(huán)境[44]。不同海拔梯度云南松群體針葉、種實(shí)等性狀多樣性指數(shù)分析表明，云南松群體表現(xiàn)出豐富的表型多樣性[22,45]。此外，在不同地理分布區(qū)域群體間、不同莖干類型表型多樣性的研究均揭示云南松表型變異豐富，各表型性狀表現(xiàn)出不一樣的變異式樣[23-25,46-48]。豐富的表型變異是云南松適應(yīng)性廣、抗逆性強(qiáng)的基礎(chǔ)[25]。云南松分布范圍廣，且為異交的多年生植物，這些特性都有利于產(chǎn)生豐富的變異，從而適應(yīng)不同環(huán)境條件[23,25]。

各球果性狀在群體間存在極顯著差異，存在一定的形態(tài)分化，球果重、球果長、球果粗和球果長/球果重的形態(tài)分化分別為0.124 1、0.051 5、0.162 1和0.026 4，以球果粗和球果重的較高，總體來看，球果性狀形態(tài)分化都較低，群體內(nèi)變異是云南松球果變異的主要來源[23,25,49]，表型變異分析與遺傳多樣性分析一致[49]。在白皮松[10]、云杉[40]、川西云杉[37]、油松[41]、青海云杉[38-39]等研究中均揭示變異主要存在于群體內(nèi)。這可能與繁育系統(tǒng)有關(guān)，云南松為風(fēng)媒異花授粉植物，種子具翅，花粉具氣囊，可遠(yuǎn)距離擴(kuò)散[23]，從而導(dǎo)致群體間的基因流較大[49]，削弱群體間的遺傳分化，從而使群體間的變異降低。

不同性狀群體內(nèi)的變異程度有所差異，球果重的變異大于球果形狀大小的變異，其中又以球果長的變異高于球果粗的變異，云南松球果大小在縱向上的變異較橫向上的變異大，以球果粗最為穩(wěn)定，而球果重的遺傳變異最為豐富。對紅松的研究也報道相似的結(jié)果，張恒慶等對不同緯度和海拔下3個生境的天然紅松林研究表明，3個群體的變異系數(shù)均表現(xiàn)為球果重>球果長>球果粗，平均變異系數(shù)分別為30.77%、14.94%和10.60%，群體內(nèi)以球果粗最為穩(wěn)定[42]。研究還表明，云南松各球果性狀比針葉性狀變異幅度小[25]，即球果性狀比針葉性狀穩(wěn)定，這可能是因?yàn)槿~子性狀更容易受環(huán)境條件的影響[14]，相反，繁殖性狀更保守不易受環(huán)境條件的影響[50]，果實(shí)大小的性狀更為穩(wěn)定[20]。此外，不同群體的各個球果性狀的變異也不一樣，在9個群體中，BH和MD群體變異程度比較高，而TC群體的變異比較低。不同群體的相同性狀變異程度存在差異，這可能是由環(huán)境異質(zhì)性導(dǎo)致的[20]。

球果性狀與各土壤因子相關(guān)不顯著或弱相關(guān)，變異可能是隨機(jī)的分布。類似的如油松[41]、白皮松[10]、川西云杉[37]、華北落葉松[43]等研究中發(fā)現(xiàn)，球果性狀地理變異趨勢不明顯。聚類分析表明，有些地理上相近的群體并沒有優(yōu)先聚合成一類，在針葉性狀的分析中也表明，針葉性狀隨地理分布的變異趨勢不明顯，而與其生態(tài)因子分異間的關(guān)系更明顯[25]，這可能與群體間立地條件類型存在差異有關(guān)[20]，即基于球果性狀的云南松群體并不完全按地理分布進(jìn)行聚類，不宜利用地理分布區(qū)域?qū)υ颇纤汕蚬硇瓦M(jìn)行劃分。

綜上所述，云南松9個群體球果性狀在群體間和群體內(nèi)均存在極顯著差異，表明變異存在的廣泛性。各個性狀的變異主要存在于群體內(nèi)，開展種質(zhì)資源收集與保護(hù)時，可適當(dāng)增加群體內(nèi)個體數(shù)、減少群體取樣數(shù)。由各個球果性狀形態(tài)分化較低，兩兩群體間的形態(tài)分化說明群體間的分化較弱。相關(guān)分析表明，各球果性狀與土壤因子間的相關(guān)性較弱，變異可能是隨機(jī)的?；谇蚬誀畹闹鞒煞址治龊途垲惙治霰砻?，不同群體間的聚類并不完全按照地理分布聚類。

[1]中國科學(xué)院昆明植物研究所編著.云南植物志(第四卷):種子植物[M].北京:科學(xué)出版社,1986:54-57.

[2]中國科學(xué)院中國植物志編輯委員會.中國植物志(第7卷)[M].北京:科學(xué)出版社,1978:122-282.

[3]金振洲,彭鑒.云南松[M].昆明:云南科技出版社,2004:1-66.

[4]鄧喜慶,皇寶林,溫慶忠,等.云南松林資源動態(tài)研究[J].自然資源學(xué)報,2014,29(8):1 411-1 419.

[5]Zhang L,Xu WH,Ouyang ZY,et al.Determination of priority nature conservation areas and human disturbances in the Yangtze River Basin,China[J].Journal for Nature Conservation,2014,22(4):326-336.

[6]李本德,劉中天.云南松(Pinusyunnanensis)林的分布及其分區(qū)[J].云南大學(xué)學(xué)報,1984(1):33-46.

[7]陳飛,王健敏,陳曉鳴,等.基于Kira指標(biāo)的云南松氣候適宜性分析[J].林業(yè)科學(xué)研究,2012,25(5):576-581.

[8]陳飛,王健敏,孫寶剛,等.云南松的地理分布與氣候關(guān)系[J].林業(yè)科學(xué)研究,2012,25(2):163-168.

[9]SevΙk H,Ayan S,Turna I,et al.Genetic diversity among populations in Scotch pine (PinussilvestrisL.)seed stands of Western Black Sea Region in Turkey[J].African Journal of Biotechnology,2010,9(43):7 266-7 272.

[10]李斌,顧萬春,盧寶明.白皮松天然群體種實(shí)性狀表型多樣性研究[J].生物多樣性,2002,10(2):181-188.

[11]賈子瑞,張守攻,王軍輝.林芝云杉天然群體針葉與種實(shí)的變異及其地理趨勢[J].林業(yè)科學(xué)研究,2011,24(4):428-436.

[12]廖柏勇,陳麗君,王芳,等.苦楝種源果核和種子性狀地理變異的趨勢面分析[J].林業(yè)科學(xué)研究,2016,29(5):784-792.

[13]Boulli A,Baaziz M,M’Hirit O.Polymorphismof natural populations ofPinushalepensisMill.in Morocco as revealed by morphological characters[J].Euphytica,2001,119(3):309-316.

[14]Dangasuk OG,Panetsos KP.Altitudinal and longitudinal variations inPinusbrutia(Ten.)of Crete Island,Greece:some needle,cone and seed traits under natural habitats[J].New Forest,2004,27(3):269-284.

[15]Donahue JK,Upton JL.Geographical variation in leaves,cone and seeds of Pinus gregii in native forests[J].Forest Ecology and Management,1996,82(1):145-157.

[16]Wahid N,González-Martínez SC,El Hadrami I,et al.Variation of morphological traits in natural populations of maritime pine (PinuspinasterAit.) in Morocco[J].Annals of forest Science,2006,63(1):83-92.

[17]尚帥斌,郭俊杰,王春勝,等.海南島青梅天然居群表型變異[J].林業(yè)科學(xué),2015,51(2):154-162.

[18]劉貴峰,臧潤國,劉華,等.天山云杉種子形態(tài)性狀的地理變異[J].應(yīng)用生態(tài)學(xué)報,2012,23(6):1 455-1 461.

[19]楊玲,沈海龍,梁立東,等.不同產(chǎn)區(qū)野生花楸果實(shí)和種子的表型多樣性[J].東北林業(yè)大學(xué)學(xué)報,2009,37(2):8-10.

[20]董昕,王磊,魯儀增,等.山東稀有植物小果白刺天然群體表型變異研究[J].林業(yè)科學(xué)研究,2017,30(2):293-299.

[21]李帥鋒,蘇建榮,劉萬德,等.思茅松天然群體種實(shí)表型變異[J].植物生態(tài)學(xué)報,2013,37(11):998-1 009.

[22]徐楊,鄧麗麗,周麗,等.云南松不同海拔天然群體種實(shí)性狀表型多樣性研究[J].種子,2015,34(11):70-74,79.

[23]許玉蘭.云南松天然群體遺傳變異研究[D].北京林業(yè)大學(xué),2015.

[24]鄧麗麗,朱霞,和潤喜,等.云南松不同莖干類型種實(shí)性狀表型多樣性比較[J].種子,2017,36(3):4-9.

[25]許玉蘭,蔡年輝,陳詩,等.基于針葉性狀云南松天然群體表型分化研究[J].西南林業(yè)大學(xué)學(xué)報,2016,36(5):1-9.

[26]葛頌,王明麻,陳岳武.用同工酶研究馬尾松群體的遺傳結(jié)構(gòu)[J].林業(yè)科學(xué),1998,24(11):399-409.

[27]劉永紅,高桂琴,崔嵬,等.油松天然群體種實(shí)性狀表型多樣性分析[J].種子,2010,29(9):44-49.

[28]López R,López de Heredia U,Collada C,et al.Vulnerability to cavitation,hydraulic efficiency,growth and survival in an insular pine (Pinuscanariensis)[J].Annals of Botany,2013,111(6):1 167-1 179.

[29]Raeymaekers JAM,Van Houdt JKJ,Larmuseau MHD,et al.Divergent selection as revealed by PST and QTL-based FST in three-spined stickleback(Gasterosteus aculeatus)populations along a coastal-inland gradient[J].Molecular Ecology,2007,16(4):891-905.

[30]Leinonen T,O'hara RB,Cano JM,et al.Comparative studies of quantitative trait and neutral marker divergence:a meta-analysis[J].Journal of Evolutionary Biology,2008,21(1):1-17.

[31]Spitze K.Population-structure in Daphnia obtusa:quantitative genetic and allozymic variation[J].Genetics,1993,135(2):367-374.

[32]Boratyńska K,Jasińska AK,Ciepuch E.Effect of tree age on needle morphology and anatomy of Pinus uliginosa and Pinus silvestris-species-specific character separation during ontogenesis[J].Flora,2008,203(8):617-626.

[33]Steel RG,Torrie JH,Dickey DA.Principles and procedures of statistics.New York-Toronto-London:McGrow-Hill Book Company Inc.,1960:481.

[34]Rohlf F J.NTSYS-PC:Numerical Taxonomy and Multivariate Analysis System Version 1.80.Setauket New York:Distribution by Exeter SoftWare.1994.

[35]張振,張含國,周宇,等.紅松多無性系群體的種實(shí)性狀變異研究[J].北京林業(yè)大學(xué)學(xué)報,2015,37(2):67-78.

[36]柴春山,蘆娟,蔡國軍,等.文冠果人工種群的果實(shí)表型多樣性及其變異[J].林業(yè)科學(xué)研究,2013,26(2):181-191.

[37]辜云杰,羅建勛,吳遠(yuǎn)傳,等.川西云杉天然種群表型多樣性[J].植物生態(tài)學(xué)報,2009,33(2):291-301.

[38]王婭麗,李毅,陳曉陽.祁連山青海云杉天然群體的種實(shí)性狀遺傳多樣性研究[J].林業(yè)科學(xué),2008,44(2):70-77.

[39]王婭麗,李毅.祁連山青海云杉天然群體的種實(shí)性狀表型多樣性[J].植物生態(tài)學(xué)報,2008,32(2):355-362.

[40]羅建勛,顧萬春.云杉天然群體表型多樣性研究[J].林業(yè)科學(xué),2005,41(2):66-83.

[41]劉永紅,高桂琴,崔嵬,等.油松天然群體種實(shí)性狀表型多樣性分析[J].種子,2010,29(9):44-49.

[42]張恒慶,安利佳,祖元剛.天然紅松種群形態(tài)特征地理變異的研究[J].生態(tài)學(xué)報,1999,19(6):932-938.

[43]李文榮,齊力旺,韓有志.山西華北落葉松天然林的地理分布和種群變異規(guī)律的研究[J].林業(yè)科學(xué),1992,28(6):493-501.

[44]Ji M,Zhang Q,Deng J,et al.Intra-versus inter-population variation of cone and seed morphological traits ofPinustabulaeformisCarr.in northern China:impact of climate-related conditions[J].Polish Journal of Ecology,2011,59(4):717-727.

[45]徐楊,周麗,蔡年輝,等.云南松不同海拔群體的針葉性狀表型多樣性研究[J].云南農(nóng)業(yè)大學(xué)學(xué)報(自然科學(xué)),2016,31(1):109-114.

[46]Xu YL,Woeste K,Cai NH,et al.Variation in needle and cone traits in natural populations ofPinusyunnanensis[J].Journal of Forestry Research,2016,27(1):41-49.

[47]鄧麗麗,孫琪,許玉蘭,等.云南松不同莖干類型群體針葉性狀表型多樣性比較[J].西南林業(yè)大學(xué)學(xué)報,2016,36(3):30-37.

[48]鄧麗麗,周麗,蔡年輝,等.基于針葉性狀的云南松不同莖干類型遺傳變異分析[J].西南農(nóng)業(yè)學(xué)報,2017,30(3):530-534.

[49]Xu YL,Cai NH,Woeste K,et al.Genetic diversity and population structure ofPinusyunnanensisby simple sequence repeat markers[J].Forest Science,2016,62(1):38-47.

[50]Vasudeva R,Hanumantha M,Gunaga RP.Genetic variation for floral traits among teak (TectonagrandisLinn.f.) clones:Implications to seed orchard fertility[J].Current Science,2004,87(3):358-362.