秦之曠 劉娜 周霞 李志輝 曹福亮 李何

摘?要：赤皮青岡（Quercus gilva）是我國(guó)中亞熱帶地區(qū)極具經(jīng)濟(jì)價(jià)值的珍貴用材樹(shù)種。為揭示其表型多樣性與變異規(guī)律以及影響表型變異的關(guān)鍵地理氣候因子，該研究以14個(gè)赤皮青岡天然種群115個(gè)單株為材料，對(duì)其15個(gè)樹(shù)體和葉片性狀進(jìn)行測(cè)量，并采用巢式方差分析、表型分化系數(shù)分析、多樣性指數(shù)分析、相關(guān)性分析、主成分分析以及聚類分析等方法，探究赤皮青岡的表型變異規(guī)律及其與地理氣候因子的相關(guān)性。結(jié)果表明：（1）赤皮青岡表型多樣性水平較高，15個(gè)表型性狀的變異系數(shù)和Shannon-Wiener指數(shù)平均值分別為35.070%和1.998，14個(gè)種群的變異系數(shù)在14.94%（洞口）～35.56%（龍山）之間、Shannon-Wiener指數(shù)在1.127（松陽(yáng)）～1.980（常寧）之間。（2）15個(gè)表型性狀在種群間和種群內(nèi)均存在極顯著差異（P<0.01），種群平均表型分化系數(shù)為41.491%，表型變異主要來(lái)源于種群內(nèi)。（3）相關(guān)性分析發(fā)現(xiàn)，赤皮青岡部分性狀間存在顯著或極顯著相關(guān)性，表型性狀與地理因子間相關(guān)性不明顯，降水是影響赤皮青岡表型性狀的主要?dú)夂蛞蜃?。?）主成分分析發(fā)現(xiàn)，基于11個(gè)表型可提煉出4個(gè)主成分，累計(jì)貢獻(xiàn)率為82.961%。（5）根據(jù)聚類分析可將14個(gè)種群分為三類，分別為小葉片小葉柄類、大葉片大葉柄類和中葉片小葉柄類，種群的聚類結(jié)果與其地理位置間無(wú)明顯相關(guān)性。該研究結(jié)果為赤皮青岡種質(zhì)資源的保護(hù)和利用提供了科學(xué)依據(jù)，為其良種選育奠定了重要基礎(chǔ)。

關(guān)鍵詞： 赤皮青岡， 珍貴樹(shù)種， 種質(zhì)資源， 表型性狀， 變異規(guī)律

中圖分類號(hào)：Q948.12

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1000-3142（2023）09-1622-14

收稿日期：2022-09-15

基金項(xiàng)目：國(guó)家自然科學(xué)基金 （32201589）； 湖南省自然科學(xué)基金 （2021JJ41069）； 湖南省林業(yè)科技創(chuàng)新資金項(xiàng)目 （XLKY202218）； 福建省公益類科研院所專項(xiàng)（2022R1010007）。

第一作者： 秦之曠（1997-），碩士，研究方向?yàn)榱帜居N，（E-mail）1667830341@qq.com。

*通信作者：李何，博士，講師，碩士研究生導(dǎo)師，研究方向?yàn)榱帜居N，（E-mail）lihecsuft@126.com。

Phenotypic diversity of Quercus gilva natural

populations in middle subtropical China

QIN Zhikuang1， LIU Na1， ZHOU Xia1， LI Zhihui1， CAO Fuliang2， LI He1*

（ 1. College of Forestry， Central South University of Forestry and Technology， Changsha 410004， China;

2. College of Forestry， Nanjing Forestry University， Nanjing 210037， China ）

Abstract：Quercus gilva is a precious timber tree species in middle subtropical China. In order to investigate the phenotypic diversity， phenotypic variation pattern and the major geographic and climatic factors influencing phenotypic variation of Q. gilva populations， 15 growth and leaf traits of 115 individuals from 14 natural populations were measured， and then the phenotypic diversity， variation pattern and the correlation between phenotypic variation and geographic and climatic factors were studied using nested analysis of variance， phenotypic differentiation coefficient， diversity index， correlation analysis， principal component analysis and cluster analysis. The results were as follows： （1） The average variation coefficient and the average Shannon-Wiener index of 15 phenotypic traits were 35.070% and 1.998， respectively. Variation coefficient and Shannon-Wiener index of 14 populations ranged from 14.94% （Dongkou） to 35.56% （Longshan） and from 1.127 （Songyang） to 1.980 （Changning）， respectively. These results indicated that Q. gilva showed a relatively high level of phenotypic diversity. （2） Significant differences in 15 phenotypic traits were found among and within populations （P<0.01）， with an average phenotypic differentiation coefficient of 41.491%. This result elucidated that the phenotypic variation within populations was greater than that among populations. （3） Correlation analysis showed that there were significant or entremely significant correlations among parts of the traits， whereas no correlation was found between phenotypic traits and geographic factors. Precipitation was the major factor that influenced the phenotypic traits of Q. gilva. （4） Principle component analysis indicated that the first four principal components explained 82.961% of total phenotypic contribution rate of Q. gilva. （5） Cluster analysis revealed that 14 Q. gilva populations were divided into three groups according to their different leaf and petiole size， such as small-leaved micropetiole， large-leaved macropetiole and medium-leaved micropetiole， which was not related to the geographic distribution of populations. This study provides scientific basis for the protection and utilization of Q. gilva germplasm resources， and lays an important foundation for its improved variety breeding.

Key words： Quercus gilva， precious tree species， germplasm resources， phenotypic traits， variation pattern

赤皮青岡（Quercus gilva）是殼斗科青岡屬常綠喬木，其樹(shù)干通直，邊材呈黃褐色，心材呈暗紅褐色，紋理優(yōu)美具有光澤，強(qiáng)韌有彈性，有獨(dú)特香味，抗蟲(chóng)蛀，不易腐蝕，是我國(guó)亞熱帶地區(qū)極具經(jīng)濟(jì)價(jià)值和發(fā)展前景的珍貴用材樹(shù)種之一（端木，1995；Zaynab et al.， 2018）。其樹(shù)體高大，為深根性樹(shù)種，在水土保持和涵養(yǎng)水源方面有良好的表現(xiàn)，從而具有很高的生態(tài)價(jià)值（夏蟬，2021）。赤皮青岡在我國(guó)中亞熱帶地區(qū)分布較廣，主要分布于湖南、福建、浙江、江西、貴州、臺(tái)灣等?。╖eng et al.， 2019；歐陽(yáng)澤怡等，2021）。然而，從20世紀(jì)開(kāi)始，赤皮青岡天然群落受到明顯的人為干擾，絕大部分赤皮青岡天然林被過(guò)度采伐，采伐跡地用于城鄉(xiāng)發(fā)展、筑路修橋或人工林營(yíng)造，導(dǎo)致大多赤皮青岡處于散生和小群落生存狀態(tài)，其種質(zhì)資源大量流失（Deng et al.， 2018）。

種質(zhì)資源既是遺傳多樣性的載體，也是種質(zhì)創(chuàng)新和林木遺傳改良的物質(zhì)基礎(chǔ)，闡明種質(zhì)資源的多樣性是對(duì)其合理保護(hù)和利用的前提（Bai et al.， 2019）。植物表型是遺傳因子和環(huán)境在長(zhǎng)期過(guò)程中共同作用的結(jié)果，表型多樣性是遺傳多樣性的直接體現(xiàn)（李洪果等，2019；郭燕等，2022）。利用植物的表型性狀檢測(cè)遺傳多樣性具有簡(jiǎn)便、快捷的特點(diǎn)（尚帥斌等，2015），是長(zhǎng)期以來(lái)最基本和最直接的方法。近年來(lái)，種質(zhì)資源的表型性狀分析已在閩楠（Phoebe bournei）（Li et al.， 2018；歐漢彪等，2021）、黑楊（Populus nigra）（陳存等，2021）、板栗（Castanea mollissima）（郭燕等，2022）、茶條槭（Acer ginnala）（武艷虹等，2018）、唐古特白刺（Nitraria tanguticum）（吉德娟等，2021）等林木上開(kāi)展，通過(guò)分析各表型性狀，揭示了表型多樣性和變異，挖掘了影響表型變異的關(guān)鍵因子，為相關(guān)樹(shù)種種質(zhì)資源的保護(hù)和利用奠定了重要基礎(chǔ)。

目前，赤皮青岡天然種群仍處于持續(xù)減少的狀態(tài)，天然種質(zhì)資源不斷流失，調(diào)查與闡明天然種質(zhì)的多樣性以合理制定其保護(hù)和利用策略是當(dāng)前亟待開(kāi)展的重要研究。但是，關(guān)于赤皮青岡種質(zhì)資源的相關(guān)研究非常有限，僅謝?。?011）和林文俊等（2017）對(duì)赤皮青岡天然群落進(jìn)行了調(diào)查，并且調(diào)查群落僅分布在閩清和擎天巖村，很難具有代表性。在分子層面，雖然有基于ISSR分子標(biāo)記技術(shù)（朱品紅，2014）和葉綠體基因組測(cè)序技術(shù)（Sugiura et al.， 2015）的赤皮青岡遺傳多樣性和遺傳結(jié)構(gòu)的相關(guān)報(bào)道，但所調(diào)查的種群分布較集中，不足以揭示我國(guó)赤皮青岡種質(zhì)的遺傳多樣性；在表型層面，僅謝?。?021）對(duì)4個(gè)種源赤皮青岡幼林的生長(zhǎng)性狀進(jìn)行了分析。樹(shù)高、胸徑、材積和干型等是用材樹(shù)種的重要經(jīng)濟(jì)性狀，而葉片是植物進(jìn)行光合和蒸騰作用的重要媒介，其形態(tài)直接影響生長(zhǎng)發(fā)育和繁殖特性，并最終影響樹(shù)種的適應(yīng)和進(jìn)化（郭燕等，2022）。

本研究以具有代表性的、基本覆蓋赤皮青岡在大陸地區(qū)分布范圍的14個(gè)赤皮青岡天然種群115個(gè)單株為材料，對(duì)其15個(gè)樹(shù)體和葉片性狀進(jìn)行測(cè)量，采用巢式方差分析、表型分化系數(shù)分析、多樣性指數(shù)分析、相關(guān)性分析、主成分分析以及聚類分析等方法，通過(guò)分析表型多樣性和變異及其地理變化，擬探討：（1）赤皮青岡天然種群表型多樣性水平；（2）赤皮青岡表型多樣性呈現(xiàn)的變異規(guī)律；（3）赤皮青岡表型變異與地理氣候因子之間的相關(guān)性。以期為赤皮青岡種質(zhì)資源的收集、保護(hù)和利用提供科學(xué)依據(jù)，為其良種選育奠定重要基礎(chǔ)。

1?材料與方法

1.1 材料

基于前期對(duì)赤皮青岡天然分布區(qū)內(nèi)種質(zhì)資源的調(diào)查，2021年5月至2022年3月對(duì)14個(gè)赤皮青岡天然種群共115個(gè)單株進(jìn)行了調(diào)查和樣本采集，采樣種群基本覆蓋了我國(guó)中亞熱帶赤皮青岡天然分布區(qū)，各種群的地理氣候信息和樣本數(shù)量如表1所示。對(duì)樣本數(shù)大于10個(gè)的種群，隨機(jī)選取10株樣本，并確保樣本間較大的距離；由于部分種群樣本較少或采樣困難，因此對(duì)樣本數(shù)小于10株的種群，則采集種群內(nèi)的所有樣本。用GPS記錄各樣本的經(jīng)度和緯度，利用測(cè)高儀、胸徑尺等工具測(cè)量且記錄樹(shù)體性狀。用高枝剪等工具在當(dāng)?shù)?2：00—14：00時(shí)樹(shù)冠中下層的陽(yáng)面，獲取健康枝條，隨機(jī)采集成熟葉片（李珊等，2016），并選取10片無(wú)病蟲(chóng)害的完整葉片進(jìn)行集中拍照、記錄。

1.2 氣候數(shù)據(jù)的獲取

從全球氣候數(shù)據(jù)網(wǎng)WorldClim（https：//www.worldclim.org/）下載當(dāng)前時(shí)期氣候數(shù)據(jù)，共下載與溫度和降水相關(guān)的19個(gè)氣候因子，數(shù)據(jù)精度為2.5’，版本為2.1（Fick & Hijmans， 2017），使用R包Raster中extract函數(shù)提取種群所在位置的氣候數(shù)據(jù)。在19個(gè)氣候因子中，由于一些因子間可能存在高度相關(guān)性，因此用prcomp函數(shù)對(duì)種群所在地的氣候因子進(jìn)行主成分分析，提取前三個(gè)主成分中載荷大于0.8的氣候因子，得到5個(gè)氣候因子（表1），用于相關(guān)性分析。

1.3 表型性狀的測(cè)定

對(duì)各種質(zhì)的樹(shù)高、胸徑、冠徑等樹(shù)體性狀進(jìn)行測(cè)量，其中冠徑為東西冠徑和南北冠徑的平均值，每個(gè)性狀重復(fù)測(cè)量3次，取平均值作為實(shí)測(cè)值。

計(jì)算各種質(zhì)的單株材積（v），公式為v=0.4π×（d1.3/2）×h+3×10（式中，d1.3為胸徑，h為樹(shù)高）（林昌庚，1964）。

記錄樹(shù)干干型情況，對(duì)這一非數(shù)值性狀由3人分別進(jìn)行1～10賦值統(tǒng)計(jì)。根據(jù)劉雨等（2011）的研究且結(jié)合赤皮青岡實(shí)際情況進(jìn)行以下賦值：有明顯主干無(wú)分枝，樹(shù)干通直記為10；有明顯主干無(wú)分枝，但樹(shù)干輕微彎曲記為9；有明顯主干無(wú)分枝，但樹(shù)干彎曲明顯記為8；有明顯主干且出現(xiàn)1個(gè)分岔，但分岔和主干較為通直記為7；有明顯主干和1個(gè)分岔，樹(shù)干較為彎曲記為6；有明顯主干和多個(gè)分岔，主干分枝較直記為5；有主干但第1個(gè)分岔較早，枝下高小于1.3 m，主干分枝較直記為4；有主干但第1個(gè)分岔較早，枝下高小于1.3 m，主干分枝較彎曲記為3；無(wú)明顯主干，但主要枝條較直記為2；無(wú)明顯主干，枝條較彎曲記為1。

使用ImageJ 圖像分析軟件分析各種質(zhì)的葉片圖像，測(cè)量獲取葉片的葉長(zhǎng)、葉寬、葉面積、葉柄長(zhǎng)、葉柄寬、葉柄面積、鋸齒數(shù)等數(shù)據(jù)，重復(fù)3次，取平均值作為實(shí)測(cè)值。通過(guò)測(cè)量值計(jì)算葉長(zhǎng)寬比、葉柄長(zhǎng)寬比和葉脈間距。其中，葉脈間距計(jì)算公式為葉脈間距=葉長(zhǎng)/（單邊葉脈數(shù)+1）（郭燕等，2020）。

1.4 數(shù)據(jù)的處理和分析

1.4.1 巢式方差分析和多重比較?在SPSS 26軟件中，采用巢式設(shè)計(jì)模型對(duì)赤皮青岡的15個(gè)表型性狀進(jìn)行巢式方差分析，其線性模型為Yijk=m+ai+bj（i）+eijk。式中， Yijk為第 i個(gè)種群的第 j個(gè)單株的第 k個(gè)觀測(cè)值， m為總體平均值， ai為第 i個(gè)種群的效應(yīng)值，bj（i）為第 i個(gè)種群的第 j個(gè)單株的效應(yīng)值， eijk為第ijk個(gè)觀測(cè)值的試驗(yàn)誤差。

1.4.2 多樣性指數(shù)分析?在Excel軟件中，通過(guò)公式CV＝δ/x× 100%和VST=δ2t/（δ2t+δ2s）計(jì)算各性狀的變異系數(shù)（CV）及表型分化系數(shù)（VST）（葛頌等，1988）。式中，δ2t表示種群間方差，δ2s表示種群內(nèi)方差。在Past 4.09中計(jì)算各種群表型性狀的Shannon-Wiener多樣性指數(shù)，先根據(jù)平均值（mean，x）和標(biāo)準(zhǔn)差（SD，δ）將所有性狀數(shù)據(jù)分成 10級(jí)，轉(zhuǎn)換為質(zhì)量性狀，第一級(jí)為< x-2δ，中間每級(jí)增加0.5δ，直到第10級(jí)為≥ x+2δ，再計(jì)算Shannon-Wiener多樣性指數(shù)。

1.4.3 相關(guān)性分析、主成分分析與聚類分析?使用SPSS 26軟件，通過(guò)Pearson相關(guān)系數(shù)分析法，研究赤皮青岡各表型性狀間以及性狀與地理位置和氣候因子間的相關(guān)性。在R軟件中，先通過(guò)prcomp函數(shù)表型數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化，再進(jìn)行主成分分析?；谥鞒煞值梅窒禂?shù)矩陣，提取貢獻(xiàn)率最大的第一和第二主成分，采用平方歐式距離和組間連接法進(jìn)行種群間聚類分析。赤皮青岡天然種群受破壞較嚴(yán)重，大多種群內(nèi)樣本數(shù)量較少，并且種群間樹(shù)齡結(jié)構(gòu)差異較大，靖州和松陽(yáng)種群僅由3株古樹(shù)和少量次生幼樹(shù)組成，保靖、通道、正板、插柳、安化、龍山和周寧等種群內(nèi)僅存幾株古樹(shù)和中齡赤皮青岡，鳳凰、洞口、常寧、婺源和普陀種群為樹(shù)齡較大的同齡林。因此，在進(jìn)行主成分分析、聚類分析以及表型性狀與地理和氣候因子相關(guān)性分析時(shí)，剔除樹(shù)高、胸徑、冠徑和材積4個(gè)受樹(shù)齡影響較大的表型性狀，使用其余11個(gè)表型性狀進(jìn)行分析。

2?結(jié)果與分析

2.1 赤皮青岡表型性狀變異

赤皮青岡15個(gè)表型性狀的變異系數(shù)平均值為35.070%（表2）。其中，材積的變異系數(shù)最大（為152.606%），變異幅度為0.05～58.72 m3。葉片長(zhǎng)寬比變異幅度最?。?1.868%），變異幅度為2.540～4.534。15個(gè)表型性狀的變異系數(shù)由大到小依次為材積＞胸徑＞冠徑＞樹(shù)高＞葉柄面積＞葉面積＞葉柄長(zhǎng)寬比＞葉柄長(zhǎng)＞葉脈間距＞鋸齒數(shù)＞葉長(zhǎng)＞葉柄寬＞葉寬＞干型＞葉長(zhǎng)寬比。這說(shuō)明赤皮青岡樹(shù)體性狀變異更加豐富。

15個(gè)表型性狀的多樣性指數(shù)（H）在1.543～2.142之間，平均值為1.998（表2）。其中，材積的多樣性指數(shù)最低（為1.543），葉寬的多樣性指數(shù)最高（為2.142）。15個(gè)表型性狀的多樣性由高到低分別為葉寬＞葉脈間距＞葉長(zhǎng)＞葉長(zhǎng)寬比＞鋸齒數(shù)＞葉柄寬＞樹(shù)高＞葉柄長(zhǎng)＞葉面積＞葉柄長(zhǎng)寬比＞葉柄面積＞冠徑＞胸徑＞干型＞材積。這說(shuō)明赤皮青岡在個(gè)體水平上存在較為豐富的表型變異。

2.2 赤皮青岡表型變異來(lái)源及種群變異規(guī)律

赤皮青岡15個(gè)表型性狀在種群間和種群內(nèi)的方差分析結(jié)果見(jiàn)表3。由表3的F檢驗(yàn)結(jié)果可知，無(wú)論是在種群間還是種群內(nèi)15個(gè)表型性狀都存在極顯著差異（P<0.01），說(shuō)明赤皮青岡的樹(shù)體和葉片性狀在種群間和種群內(nèi)的變異都較為豐富。

由表4可知，保靖、插柳、安化、龍山、婺源、周寧、松陽(yáng)種群的樹(shù)高、胸徑、冠徑和干型等樹(shù)體指標(biāo)值都較大，其中松陽(yáng)種群的樹(shù)高、材積、冠徑、干型均為最大，表現(xiàn)為樹(shù)形高大通直的特征；靖州種群的胸徑最大；鳳凰、洞口、常寧、普陀種群的平均材積相對(duì)較小，與各種群實(shí)際樹(shù)齡結(jié)構(gòu)相符。靖州、常寧、普陀種群的葉長(zhǎng)寬比較小，葉片表現(xiàn)為短而寬的葉形；保靖、龍山、婺源種群的葉長(zhǎng)寬比較大，葉片呈細(xì)長(zhǎng)型；插柳種群的葉長(zhǎng)和葉面積最大；正板赤皮青岡的葉片鋸齒數(shù)最多，葉面積較大，表現(xiàn)為葉片大且鋸齒多的特點(diǎn)。

2.3 赤皮青岡種群內(nèi)表型變異特征

變異系數(shù)（CV）表示種群內(nèi)表型性狀的變異情況，變異系數(shù)越大表明表型的離散程度越大。由表5可知，14個(gè)種群的變異系數(shù)在14.94%（洞口）～35.56%（龍山）之間，平均值為23.07%，說(shuō)明在種群水平上赤皮青岡的表型性狀變異水平較高；14個(gè)種群中，僅洞口種群的15個(gè)表型性狀變異系數(shù)平均值未超過(guò)15%，說(shuō)明該種群赤皮青岡表型變異水平較低；15個(gè)表型性狀中干型、葉長(zhǎng)、葉寬、葉長(zhǎng)寬比、葉柄寬和鋸齒數(shù)的變異系數(shù)平均值都未超過(guò)15%，說(shuō)明這6個(gè)性狀較為穩(wěn)定。與表2結(jié)果一致，除干型外，樹(shù)高、胸徑、冠徑和材積等樹(shù)體性狀的變異系數(shù)平均值大于葉片性狀的變異系數(shù)平均值，說(shuō)明赤皮青岡天然種質(zhì)的樹(shù)體性狀受樹(shù)齡影響過(guò)高而表現(xiàn)出較高變異水平。

2.4 赤皮青岡種群內(nèi)表型多樣性

Shannon-Wiener多樣性指數(shù)反映不同表型性狀的豐富和均勻程度，14個(gè)種群的15個(gè)表型性狀的多樣性指數(shù)見(jiàn)表6。由表6可知，Shannon-Wiener指數(shù)在1.127（松陽(yáng)）～1.980（常寧）之間，14個(gè)種群的多樣性指數(shù)從高到低依次為常寧＞普陀＞洞口＞插柳＞保靖＞鳳凰＞周寧＞通道＞婺源＞龍山＞安化＞正板＞靖州＞松陽(yáng)。Shannon-Wiener多樣性指數(shù)與種群內(nèi)樣本數(shù)量大小有一定的關(guān)系，靖州、安化、松陽(yáng)等樣本較少的種群表現(xiàn)出較低的Shannon-Wiener指數(shù)，在種質(zhì)資源收集和保護(hù)工作中應(yīng)重視。

2.5 赤皮青岡表型性狀變異來(lái)源及種群間表型分化

表型分化系數(shù)（VST）是種群間表型變異占總表型的百分比，用于評(píng)估種群間的表型分化情況。由表7可知，15個(gè)赤皮青岡表型性狀的分化系數(shù)在23.420%（葉柄長(zhǎng)）～57.359%（胸徑）之間。赤皮青岡種群間表型方差分量占總表型變異的31.682%，種群內(nèi)占44.677%，表型分化系數(shù)為41.491%，說(shuō)明種群內(nèi)的變異在赤皮青岡表型變異中起主要作用。

2.6 赤皮青岡表型性狀間的相關(guān)性及其與地理和氣候因子間的相關(guān)性

15個(gè)赤皮青岡表型性狀間的相關(guān)性結(jié)果如表8所示。其中，部分性狀間存在顯著或極顯著相關(guān)性。樹(shù)高與胸徑、材積、干型、冠徑、葉長(zhǎng)寬比呈極顯著正相關(guān)性，與葉長(zhǎng)和葉脈間距呈顯著正相關(guān)，說(shuō)明較高的赤皮青岡往往具有較大的胸徑、較優(yōu)良的干型和冠形，葉片也較為狹長(zhǎng)。干型與胸徑、材積、冠徑均呈極顯著正相關(guān)，說(shuō)明優(yōu)良的干型在一定程度上可以提高赤皮青岡的生長(zhǎng)潛力。葉面積與葉長(zhǎng)寬比無(wú)明顯相關(guān)性，說(shuō)明不同種質(zhì)的葉長(zhǎng)寬比表現(xiàn)出一定程度的穩(wěn)定性。葉寬與葉柄長(zhǎng)、葉柄寬、葉柄面積均呈極顯著正相關(guān)，葉長(zhǎng)和葉面積與葉柄寬和葉柄面積呈極顯著正相關(guān)，說(shuō)明大的葉片通常伴隨葉柄粗長(zhǎng)的特征。葉脈間距僅與鋸齒數(shù)呈極顯著負(fù)相關(guān)，說(shuō)明鋸齒越多葉脈間距越小。

11個(gè)表型性狀與種群所在地的地理和氣候因子的相關(guān)性分析結(jié)果如表9所示。赤皮青岡表型性狀與地理氣候因子間的相關(guān)性不明顯。干型與最濕月份降水量、最濕季度降水量和最暖季度降水量呈顯著正相關(guān)，葉寬與最暖季度降水量呈顯著負(fù)相關(guān)，葉柄寬與最濕月份降水量和最暖季度降水量呈顯著負(fù)相關(guān)；經(jīng)度、緯度、最冷月份最低溫度和最濕季度平均溫度與赤皮青岡表型性狀間無(wú)相關(guān)性。這表明降水對(duì)赤皮青岡表型性狀的影響大于溫度和地理因子。

2.7 赤皮青岡表型性狀的主成分分析和聚類分析

對(duì)赤皮青岡11個(gè)表型性狀進(jìn)行主成分分析，得到4個(gè)特征根大于1的主成分，累計(jì)貢獻(xiàn)率為82.961% （表10）。由表10可知， PC1特征值為4.205，解釋總變異的38.231%， 主要包括葉長(zhǎng)、 葉寬、葉面積（載荷的絕對(duì)值＞0.8），主要表征葉片大??；PC2的特征值為2.447，解釋了22.248%的總變異，主要包括葉柄長(zhǎng)和葉柄長(zhǎng)寬比，與葉柄性狀相關(guān)；PC3特征值為1.344，主要表征干型，解釋了12.214%的總變異；PC4特征值為1.129，主要表征鋸齒數(shù)，解釋了10.268%的總變異。由圖1可知，赤皮青岡表型性狀未呈現(xiàn)出明顯的地理變異規(guī)律，其中湖南鳳凰和福建周寧種群以其較低的PC1得分與其他種群區(qū)分開(kāi)來(lái)。

以平方歐式距離為14個(gè)赤皮青岡種群間的距離，采用組間連接法進(jìn)行聚類（圖2）。當(dāng)平方歐式距離為15時(shí)，根據(jù)11個(gè)表型性狀可將14個(gè)種群分為三類（與主成分分析結(jié)果相符）：第一類為湖南鳳凰和福建周寧，主要特征為葉長(zhǎng)、葉寬、葉面積、葉柄長(zhǎng)、葉柄寬、葉柄面積和葉脈間距值均最?。ū?），表現(xiàn)為小葉片小葉柄類；第二類為湖南靖州、浙江松陽(yáng)和湖南插柳，其葉寬、葉面積、葉柄長(zhǎng)、葉柄寬、葉柄面積和葉脈間距的值最大，葉長(zhǎng)寬比值最低，說(shuō)明第二類呈現(xiàn)出葉片大、葉寬且葉柄粗長(zhǎng)的特征，表現(xiàn)為大葉片大葉柄類；第三類包括湖南保靖、湖南通道、湖南正板、湖南洞口、湖南安化、湖南龍山、湖南常寧、江西婺源和浙江普陀9個(gè)種群，屬于第一類和第二類的銜接部分，葉片形態(tài)接近第二類，葉長(zhǎng)、葉寬、葉面積的值普遍較大，但葉柄形態(tài)更接近第一類，普遍較小，總體表現(xiàn)為中葉片小葉柄類。通過(guò)比較種群的聚類結(jié)果和地理位置之間的關(guān)系，發(fā)現(xiàn)二者關(guān)聯(lián)并不顯著，說(shuō)明赤皮青岡表型性狀與地理因子間的相關(guān)性不明顯。

3?討論與結(jié)論

3.1 赤皮青岡表型性狀的多樣性

表型性狀的變異是環(huán)境和遺傳因子共同影響下的表達(dá)（曹德美等，2021），尤其是赤皮青岡這類珍稀植物，種質(zhì)資源流失嚴(yán)重且面臨復(fù)雜的生存環(huán)境和長(zhǎng)期的地理隔離，極易發(fā)生樹(shù)體和葉片等表型性狀的變異。

本研究結(jié)果表明，赤皮青岡15個(gè)表型性狀的Shannon-Wiener多樣性指數(shù)在1.543～2.142之間，平均值為1.998。高于瀕危植物云南黃連（1.634～1.839）（Coptis chinensis）（楊維澤等，2013），與瀕危植物格木（Erythrophleum fordii）（2.027 8）（李洪果等，2019）相近，低于四川珙桐（2.360）（Davidia involucrata）（張騰駒等，2019）。這表明赤皮青岡雖面臨種質(zhì)資源流失的問(wèn)題，仍表現(xiàn)出較高的多樣性水平，特別是葉長(zhǎng)、葉寬和葉脈間距等性狀。葉片的形態(tài)直接影響葉片生長(zhǎng)發(fā)育特性，間接影響植物的適應(yīng)能力，赤皮青岡葉片表型性狀的多樣性說(shuō)明其具有較強(qiáng)的環(huán)境適應(yīng)性。

本研究中，赤皮青岡15個(gè)表型性狀的極大值是極小值的1.78（葉長(zhǎng)寬比）～1 174.40（材積）倍，其中樹(shù)體性狀的極大值與極小值的比值普遍大于葉片性狀，表明樹(shù)體性狀較葉片性狀更易受到樹(shù)齡本身的影響。本研究中，赤皮青岡15個(gè)表型性狀的變異系數(shù)在11.868%～152.606%之間，天然種質(zhì)間表型性狀差異較大，尤其是樹(shù)高、胸徑、冠徑和材積等受樹(shù)齡影響大的樹(shù)體性狀；剔除以上4個(gè)性狀后，其余11個(gè)表型性狀的變異系數(shù)為11.868%～29.718%，與灰楸（Catalpa fargesii）（9.66%～30.39%） （于曉池等，2021） 和板栗（Castanea mollissima）（8.68%～32.73%）（郭燕等，2022）相近。這說(shuō)明赤皮青岡干型和葉片性狀變異豐富，在不同環(huán)境因子和遺傳因子的影響下做出相應(yīng)的改變，這與Shannon-Wiener多樣性指數(shù)的結(jié)論一致，表明赤皮青岡在復(fù)雜生長(zhǎng)環(huán)境中表現(xiàn)出較強(qiáng)的適應(yīng)能力。

3.2 赤皮青岡表型性狀的變異來(lái)源和分化

本研究中，赤皮青岡15個(gè)表型性狀在種群間和種群內(nèi)都表現(xiàn)出較豐富的變異，而種群內(nèi)的表型變異更為豐富（平均表型分化系數(shù)VST=41.491%）。這與大多數(shù)林木種質(zhì)資源的表型研究結(jié)果基本一致，如小果白刺（Nitraria sibirica）（40.71%）（董昕等，2017）、香椿（Toona sinensis）（40.209%）（吳軍等，2018）等樹(shù)種均存在較為豐富的表型變異。這表明赤皮青岡葉片等表型性狀在種群間的遺傳分化程度相對(duì)較低，可能是本研究所調(diào)查種群均位于我國(guó)亞熱帶地區(qū)，分布區(qū)內(nèi)土壤理化性質(zhì)和氣候環(huán)境等條件差別不大，并且赤皮青岡在自然分布范圍內(nèi)呈碎片化分布的時(shí)間較短，暫未因生殖隔離而形成較大的差異。

本研究中，赤皮青岡11個(gè)葉片及干型性狀無(wú)論是在種群間還是在種群內(nèi)都存在極顯著差異，而種群內(nèi)的表型變異為主要變異來(lái)源（表型分化系數(shù)在23.420%～56.757%之間）。同時(shí)，14個(gè)種群的11個(gè)葉片及干型性狀的變異系數(shù)平均值在7.35%（干型）～23.33%（葉柄面積）之間，并且同一表型性狀的Shannon-Wiener指數(shù)在種群間存在較大差異，提示該物種較易發(fā)生種群分化。本研究所調(diào)查的赤皮青岡天然種群大多零散分布于村落周圍，雖然個(gè)別天然種群保存較為完整，但大部分種群因受人為干擾較嚴(yán)重而造成天然生境的改變，加上生境不連續(xù)、樣本較少的零散分布，種群間的基因流受阻，從而引起種群分化和表型性狀變異。

3.3 赤皮青岡表型變異與地理氣候因子間的相關(guān)性

不同物種對(duì)環(huán)境因子的敏感度和適應(yīng)力有一定的差異，其表型性狀呈現(xiàn)出不同的地理變異規(guī)律，如板栗（郭燕等，2022）的葉片由南向北表現(xiàn)出長(zhǎng)寬比變小的趨勢(shì)，唐古特白刺（吉德娟等，2021）的葉片厚度與海拔呈正相關(guān)，青楊（Populus cathayana）（曹德美等，2021）的側(cè)脈數(shù)與經(jīng)緯度呈極顯著正相關(guān)。本研究中，在適宜赤皮青岡生長(zhǎng)的環(huán)境下，經(jīng)緯度和溫度對(duì)赤皮青岡表型性狀的影響不顯著，而降水與赤皮青岡表型性狀間存在一定的相關(guān)性。這與姜小龍（2020）的研究結(jié)果類似，降水是影響南嶺青岡（Quercus championii）和福建青岡（Q. chungii）分布最大的氣候因子。本研究發(fā)現(xiàn)，氣候因子中最濕月份降水量、最濕季度降水量以及最暖季度降水量?jī)H與赤皮青岡干型呈正相關(guān)，而最濕月份降水量和最暖季度降水量與葉寬和葉柄寬呈負(fù)相關(guān)。在適宜的降水范圍內(nèi)，降水量的增加會(huì)促進(jìn)維管形成層的活動(dòng)，有助于樹(shù)干的均勻增粗，從而表現(xiàn)出優(yōu)良的干型。隨著降水量的增加，櫟類樹(shù)種的葉片和葉柄傾向于長(zhǎng)條狀發(fā)展（李東勝等，2013），葉寬和葉柄寬表現(xiàn)出減小的趨勢(shì)，該現(xiàn)象與苦楝（Melia azedarach）（陳麗君等，2016）、板栗（郭燕等，2022）等的研究結(jié)果相似。因此，降水是影響赤皮青岡表型和分布最主要的因子，這為赤皮青岡的資源圃和人工林基地建設(shè)提供了科學(xué)依據(jù)。

3.4 赤皮青岡的保護(hù)與開(kāi)發(fā)利用策略

赤皮青岡天然種群的生境不連續(xù)、跨度較大、碎片化較為嚴(yán)重，目前大多天然種群處于被孤立狀態(tài)且包含的種質(zhì)數(shù)量較少。赤皮青岡種子為頑拗性種子，外形較大，自然擴(kuò)散距離十分有限，并且花粉的傳播受到我國(guó)南方復(fù)雜地形的阻礙，長(zhǎng)此以往，高片段化分布的種群間基因流受阻（姜小龍，2020），加上環(huán)境異質(zhì)性的影響，從而引起種群分化和表型性狀變異。因此，僅開(kāi)展就地保護(hù)既不能有效地為赤皮青岡種群間的基因交流提供廊道，也不利于種群的延續(xù)和遺傳多樣性的維持。為了更好地保護(hù)、管理和利用赤皮青岡種質(zhì)資源，應(yīng)加強(qiáng)對(duì)現(xiàn)有天然種群和單株的保護(hù)力度，特別是單株數(shù)量稀少的種群；加快開(kāi)展種質(zhì)資源收集工作，鑒于赤皮青岡種群內(nèi)表型變異水平較高，要注重對(duì)種群內(nèi)每個(gè)單株的收集；應(yīng)將赤皮青岡作為鄉(xiāng)土樹(shù)種在其天然分布區(qū)內(nèi)推廣造林，以適應(yīng)能力較強(qiáng)的大葉片大葉柄類的靖州、松陽(yáng)和插柳種群為種源進(jìn)行采種育苗和造林，既可提高赤皮青岡的適應(yīng)性，又可改善種群片段化分布的現(xiàn)狀，為種群間的交流提供廊道，幫助該物種的延續(xù)。

參考文獻(xiàn)：

BAI QS， CAI YL， HE BX， et al.， 2019. Core set construction and association analysis of Pinus massoniana from Guangdong Province in southern China using SLAF-seq［J］. Sci Rep， 9： 13157.

CAO DM， ZHANG YH， CHENG XQ， et al.， 2021. Genetic variation of leaf phenotypic traits in different populations of Populus cathayana ［J］. Sci Silv Sin， 57（8）： 56-67.［曹德美， 張亞紅， 成星奇， 等， 2021. 青楊不同種群葉片表型性狀的遺傳變異 ［J］. 林業(yè)科學(xué)， 57（8）： 56-67.］

CHEN C， DING CJ， HUANG QJ， et al.， 2021. Construction of phenotypic core collection of Populus deltoides ［J］. For Res， 34（2）： 1-11.［陳存， 丁昌俊， 黃秦軍， 等， 2021. 美洲黑楊表型核心種質(zhì)庫(kù)構(gòu)建 ［J］. 林業(yè)科學(xué)研究， 34（2）： 1-11.］

CHEN LJ， LIU MQ， LIAO BY， et al.， 2016. Variation of leaf morphological characters of Melia azedarach from different provenances ［J］. J SW For Univ， 36（5）： 16-20.［陳麗君， 劉明騫， 廖柏勇， 等， 2016. 苦楝不同種源葉片性狀變異研究 ［J］. 西南林業(yè)大學(xué)學(xué)報(bào)， 36（5）： 16-20.］

DENG M， JIANG XL， HIPP AL， et al.， 2018. Phylogeny and biogeography of East Asian evergreen oaks （Quercus section Quercus; Fagaceae）： Insights into the Cenozoic history of evergreen broad-leaved forests in subtropical Asia ［J］. Mol Phylogenet Evol， 119： 170-181.

DONG X， WANG L， LU YZ， et al.， 2017. Phenotypic variation of Nitraria sibirica natural populations in Shandong ［J］. For Res， 30（2）： 293-299.［董昕， 王磊， 魯儀增， 等， 2017. 山東稀有植物小果白刺天然群體表型變異研究 ［J］. 林業(yè)科學(xué)研究， 30（2）： 293-299.］

DUAN M， 1995. Comprehensive utilization of Quercus resources in China ［J］. J Beijing For Univ， 17（2）： 109-110.［端木， 1995. 我國(guó)青岡屬資源的綜合利用 ［J］. 北京林業(yè)大學(xué)學(xué)報(bào)， 17（2）： 109-110.］

FICK SE， HIJMANS RJ， 2017. WorldClim 2： new 1km spatial resolution climate surfaces for global land areas ［J］. Int J Clim， 37（12）： 4302-4315.

GE S， WANG MX， CHEN YW， 1988. An analysis of population genetic structure of Masson pine by isozyme tech-nique ［J］. Sci Silv Sin， 24（4）： 399-409． ［葛頌， 王明庥， 陳岳武， 1988. 用同工酶研究馬尾松群體的遺傳結(jié)構(gòu) ［J］. 林業(yè)科學(xué)， 24（4）： 399-409.］

GUO Y， ZHANG SH， LI Y， et al.， 2022. Diversity analysis of 36 leaf phenotypic traits of Chinese chestnut ［J］. Sci Agric Sin， 55（5）： 991-1009.［郭燕， 張樹(shù)航， 李穎， 等， 2022. 中國(guó)板栗36個(gè)葉片表型性狀的多樣性 ［J］. 中國(guó)農(nóng)業(yè)科學(xué)， 55（5）： 991-1009.］

GUO Y， ZHANG SH， LI Y， et al.， 2020. Studies on the leaf morphology， anatomical structure and drought resistance evaluation of 238 Chinese chestnut varieties （Strains） ［J］. Acta Hortic Sin， 47（6）： 1033-1046.［郭燕， 張樹(shù)航， 李穎， 等， 2020. 中國(guó)板栗238份品種（系）葉片形態(tài)、解剖結(jié)構(gòu)及其抗旱性評(píng)價(jià) ［J］. 園藝學(xué)報(bào)， 47（6）： 1033-1046.］

JIANG XL， 2020. Phylogenetic relationship and population genetic structure of Quercus chungii and Q. championii［D］. Changsha： Central South University of Forestry and Technology.［姜小龍， 2020. 福建青岡和嶺南青岡系統(tǒng)發(fā)育關(guān)系及居群遺傳結(jié)構(gòu) ［D］. 長(zhǎng)沙： 中南林業(yè)科技大學(xué).］

JI DJ， ZHANG DF， YU Q， 2021. Phenotypic diversity of Nitraria tangutorum Bobr. in Qaidam basin ［J］. J Centr S Univ For Technol， 41（10）： 57-66.［吉德娟， 張得芳， 于倩， 2021. 柴達(dá)木盆地唐古特白刺的表型多樣性 ［J］. 中南林業(yè)科技大學(xué)學(xué)報(bào)， 41（10）： 57-66.］

LIN CG， 1964. Stem control in forest volume measurement technology ［J］. Sci Silv Sin， （4）： 87-97.［林昌庚， 1964. 林木蓄積量測(cè)算技術(shù)中的干形控制問(wèn)題 ［J］. 林業(yè)科學(xué)， （4）： 87-97.］

LI YG， LIU XH， MA JW， et al.， 2018. Phenotypic variation in Phoebe bournei populations preserved in the primary distribution area ［J］. J For Res， 29（1）： 35-44.

LI DS， SHI ZM， FENG QH， et al.， 2013.Response of leaf morphometric traits of Quercus species to climate in the temperate zone of the North-South Transect of Eastern China ［J］. Acta Phytoecol Sin， 37（9）： 793-802.［李東勝， 史作民， 馮秋紅， 等， 2013. 中國(guó)東部南北樣帶暖溫帶區(qū)櫟屬樹(shù)種葉片形態(tài)性狀對(duì)氣候條件的響應(yīng) ［J］. 植物生態(tài)學(xué)報(bào)， 37（9）： 793-802.］

LI S， GAN XH， HAN HY， et al.， 2016. Leaf phenotypic traits of Tetracentron sinense， an endangered plant species ［J］. For Res， 29（5）： 687-697.［李珊， 甘小洪， 憨宏艷， 等， 2016. 瀕危植物水青樹(shù)葉的表型性狀變異 ［J］. 林業(yè)科學(xué)研究， 29（5）： 687－697．］

LIU Y， LI DW， QIN YS， et al.， 2011. Phenotypic diversity for natural populations of Juniperus rigida in the Helan Mountains of Ningxia ［J］. J Zhejiang A & F Univ， 28（4）： 619-627. ［劉雨， 李登武， 秦廷松， 等， 2011. 寧夏賀蘭山杜松天然群體的表型多樣性 ［J］. 浙江農(nóng)林大學(xué)學(xué)報(bào)， 28（4）： 619-627.］

LIN WJ， LIU QY， LUO X， et al.， 2017. Research on species diversity and niche of natural small population of Quercus glauca ［J］. Fujian For， （6）： 38-41. ［林文俊， 劉沁月， 羅蕭， 等， 赤皮青岡天然小居群物種多樣性與生態(tài)位研究 ［J］. 福建林業(yè)， （6）： 38-41.］

OU HB， DONG LJ， LI J， et al.， 2021. Relationship between phenotypic diversity of Phoebe bournei deed traits and environmental geographical factors ［J］. J NE For Univ， 49（8）： 45-50.［歐漢彪， 董利軍， 李娟， 等， 2021. 閩楠種實(shí)性狀表型多樣性與環(huán)境地理因子的關(guān)系 ［J］. 東北林業(yè)大學(xué)學(xué)報(bào)， 49（8）： 45-50.］

OUYANG ZY， OUYANG SL， WU JY， et al.， 2021. Research progress of precious commercial tree species Quercus gilva ［J］. Hunan For Sci Technol， 48（6）： 74-79.［歐陽(yáng)澤怡， 歐陽(yáng)碩龍， 吳際友， 等， 2021. 珍貴用材樹(shù)種赤皮青岡研究進(jìn)展 ［J］. 湖南林業(yè)科技， 48（6）： 74-79.］

SHANG SB， 2015. Genetic diversity in natural populations of Vatica mangachapoi on Hainan Island， China ［D］. Beijing： Chinese Academy of Forestry.［尚帥斌， 2015. 海南島青梅天然居群遺傳多樣性 ［D］. 北京： 中國(guó)林業(yè)科學(xué)研究院.］

SUGIURA N， TANG DQ， KUROKOCHI H， et al.， 2015. Genetic structure of Quercus gilva Blume in Japan as revealed by chloroplast DNA sequences ［J］. Botany， 93（12）： 873-880.

WU J， ZHONG ZZ， LOU J， et al.， 2018. A study on phenotypic diversity of cone and seed in natural populations of Toona sinensis ［J］. Acta Agric Univ Jiangxi， 40（2）： 248-256.［吳軍， 鐘志真， 樓軍， 等， 2018. 香椿天然群體種實(shí)性狀表型多樣性研究 ［J］. 江西農(nóng)業(yè)大學(xué)學(xué)報(bào)， 40（2）：248-256.］

WU YH， FAN ZL， LI J， et al.， 2018. Phenotypic diversity of seeds and fruits in natural populations of Acer ginnala in China ［J］. Guihaia， 38（6）： 795-803.［武艷虹， 樊澤璐， 李佳， 等， 2018. 茶條槭自然種群種子和果實(shí)表型多樣性研究 ［J］. 廣西植物， 38（6）： 795-803.］

XIA C， 2021. Effects of shading on the growth and physiological characteristics of Quercus gilva seedlings ［D］. Changsha： Central South University of Forestry and Technology.［夏嬋， 2021. 遮陰對(duì)赤皮青岡幼苗生長(zhǎng)及生理特性的影響 ［D］. 長(zhǎng)沙： 中南林業(yè)科技大學(xué).］

XIE J， 2021. Variations in growth traits and provenance selection of Quercus gilv ［J］. J Fujian Agric For Univ （Nat Sci Ed）， 50（5）： 636-640.［謝健， 2021. 不同種源赤皮青岡幼林生長(zhǎng)性狀變異及種源選擇 ［J］. 福建農(nóng)林大學(xué)學(xué)報(bào)（自然科學(xué)版）， 50（5）： 636-640.］

XIE J， 2011. Research on population ecology of Quercus gilvs ［D］. Fuzhou： Fujian Agriculture and Forestry University.［謝健， 2011. 赤皮青岡種群生態(tài)學(xué)研究 ［D］. 福州： 福建農(nóng)林大學(xué).］

YANG WZ， JIN H， LI WY， et al.， 2013. Phenotypic diversity of different populations of endangered plant Coptis yunnanensis ［J］. J Yunnan Univ （Nat Sci Ed）， 35（5）： 719-726.［楊維澤， 金航， 李晚誼， 等， 2013. 瀕危植物云南黃連不同居群表型多樣性研究 ［J］. 云南大學(xué)學(xué)報(bào)（自然科學(xué)版）， 35（5）： 719-726.］

YU XC， LI F， OU Y， et al.， 2021. Construction of core germplasm of Sorbus mandshurica based on phenotype ［J］. For Res， 34（6）： 38-45.［于曉池， 李鳳， 歐陽(yáng)， 等， 2021. 基于表型的灰楸核心種質(zhì)構(gòu)建 ［J］. 林業(yè)科學(xué)研究， 34（6）： 38-45.］

ZAYNAB M， PAN D， NOMAN A， et al.， 2018. Transcriptome approach to address low seed germination in Quercus gilva to save forest ecology ［J］. Biochem Syst Ecol， 81： 62-69.

ZENG QM， LIU B， LIU RQ， et al.， 2019. The complete chloroplast genome sequence of Quercus gilva （Fagaceae）［J］. Mitchondrial DNA B， 4（2）： 2493-2494.

ZHANG TJ， CHEN XH， KANG XK， et al.， 2019. Phenotypic diversity of leaf morphologic traits of Davidia involucrata natural populations in Sichuan Province ［J］. Sin J Acta Ecol， 38（1）： 35-43.［張騰駒， 陳小紅， 康喜坤， 等， 2019. 四川省珙桐天然種群葉表型多樣性 ［J］. 生態(tài)學(xué)雜志， 38（1）： 35-43.］

ZHU PH， 2014. Analysis on genetic diversity and genetic structure of Quercus gilva population ［D］. Changsha： Central South University of Forestry and Technology.［朱品紅， 2014. 赤皮青岡居群遺傳多樣性與遺傳結(jié)構(gòu)分析 ［D］. 長(zhǎng)沙： 中南林業(yè)科技大學(xué).］

（責(zé)任編輯?蔣巧媛?王登惠）