紅錐葉片功能性狀對地理環(huán)境因子的響應(yīng)

中圖分類號：Q948 文獻(xiàn)標(biāo)識碼：A 文章編號：1000-3142（2025）04-0773-14

Response of functional traits of Castanopsis hystrix leaves to geographical environment factors

LI Qiguo'， JIANG Weixin'， TAN Zhangqiang2，BAI Tiandao1*， SHEN Wenhui（1.KeyLboofalesdssdstoouiaofstmbinentluiKeyLaboratoryofForestEcologyndConservationCoegeofFrestryGuangxiUniersityNanningo，hin;.GuagiAcademyofForestrySciences/GuangxiCenterforImprovedCastanopsishystrixBreeding，Nanning53o2，China）

收稿日期：2024-07-09 接受日期：2024-09-18

基金項目：中央引導(dǎo)地方科技發(fā)展資金項目（桂科ZY21195049）。

第一作者：厲啟國（1998—），碩士研究生，主要從事林木遺傳育種研究，（E-mail）2397096665 @ qq.com。通信作者：白天道，博士，副教授，主要從事林木遺傳資源保育研究，（E-mail） btd@ gxu.edu.cn。

Abstract：Leafisanimportantnutrientorganforphysiologicalandbiochemicalactivitiesof plants，anditsfunctional traitscanreflecttheenvironmentaladaptabilityof plants toacertainextent.Inordertorevealtherelationshipbetween thefunctional traitsofCastanopsis hystrix leavesand geographicalenvironmentfactors，andtounderstand itsadaptive response to environmental changes，in this study，six populations of c .hystrix in Guangxi were used as the research objects.Various statistical methods including nested analysis of variance，principal component analysis，and cluster analysiswere employedto analyze the variationof leaf functionaltraitsandtheircorelationswith geographical environment factors，toexplorethe patemof geographicaland environmentalvariationintheleaffunctional traitsof C. hystrix. The results were as follows：（1）The 13 leaf functional traits of C .hystrix showed highly significant differences among and within populations （ Plt;0 . 0 0 1 ），with high diversity.（2）The leaf functional traitswere dominated bywithin population variation[average trait differentiation coefficient ，andthe differentiationbetween populations of morphological traits （ 1 3 . 1 1 % - 4 6 . 9 3 % ）washigher than that of physiological and biochemical traits （10.95%- 22.59 % ），indicating that the morphological traitsof leaves responded more significantly to the environment.（3）Theleaffunctionaltraits were positivelycorelated withaltitude，latitudeandlongitude，butnegatively correlated with annual precipitation，indicating that the leaf functional traits of C ：hystrixwere affectedbyaltitude， geographic distance and precipitation.（4）Cluster analysis divided the participating C .hystrix populations into two groups，and Mantel testshowedthatleaffunctional traitswere significqantlypositivelycorelated withboth geographic distance （ P=0 . 0 2 8 ）and altitude （ P=0 . 0 1 3 ），and the effects of altitude on leaf functional traits were relatively more significant.In summary，the functional traits of c .hystrix leavesare rich in variation，which mainly originates from within populations and showa patern of geographic variation dominatedbychanges in altitude，latitude，longitude and precipitation. The results of this study provides a theoretical basis for the conservation of c . hystrix germplasm resources and the selection and utilization of seed sources.

Keywords： Castanopsishystrix，functional traits，geographical environment factors， variationpaterns， association analysis

植物功能性狀是指植物在長期適應(yīng)外界環(huán)境的過程中形成的形態(tài)、生理和物候特征（Cornelissenetal.，2003；宋帥帥等，2023），與外界環(huán)境具有相互適應(yīng)性和協(xié)同演化性（王超等，2021）。葉片是植物進行光合、呼吸、蒸騰等生理活動的重要營養(yǎng)器官（Maetal.，2023），是高靈敏度的環(huán)境傳感器，最能體現(xiàn)植物的功能以及其對環(huán)境的適應(yīng)性（Garnieretal.，2004；Sandoval-Granilloamp;Meave，2023）。因此，葉片功能性狀是植物功能性狀的重要組成部分，也是植物響應(yīng)環(huán)境變化的可靠指標(biāo)（熊玲等，2024）。一般認(rèn)為，葉片的形態(tài)性狀和生理生化性狀是葉片功能性狀的重要組成部分，共同影響著植物對環(huán)境的適應(yīng)能力（孫梅等，2017），如葉長、葉寬、葉面積、葉厚等葉片形態(tài)性狀可以反映不同種群間植株的生長差異及葉片變化規(guī)律（孟益德等，2022）；可溶性糖、可溶性蛋白、丙二醛和脯氨酸等生理生化性狀與植物的抗逆性有關(guān)（黃相玲等，2018）。近年來，隨著對全球氣候變化關(guān)注度的提高，植物葉片功能性狀與地理環(huán)境因子（海拔、年均溫、降水等）的響應(yīng)研究越來越受關(guān)注。研究表明，斑皮檸檬桉（Corymbiacitriodorassp.variegata）（唐慶蘭等，2022）苦楝（Meliaazedarach）（陳麗君等，2016）、毛白楊（Populustomentosa）（Du etal.，2014）等樹木的葉長、葉寬、葉面積等功能性狀隨著海拔的升高、緯度的增加，呈現(xiàn)逐漸變大的趨勢，而蒙古櫟（Quercusmongolica）（尹杰，2023）、青楊（Populuscathayana）（曹德美等，2021）等植物則與之相反，體現(xiàn)了不同種類植物葉片對環(huán)境的適應(yīng)性策略不同（羅欣語等，2024）。因此，探明特定植物葉功能性狀與環(huán)境因子的關(guān)系，對于深人了解環(huán)境因子在驅(qū)動植物葉功能性狀變異方面的貢獻(xiàn)，以及制定合理的資源保育和經(jīng)營策略具有重要意義（劉樂樂等，2017）。

紅錐（Castanopsishystrix），為殼斗科（Fagaceae）錐屬（Castanopsis）常綠喬木，是我國華南地區(qū)的優(yōu)良鄉(xiāng)土用材林樹種，在廣西、廣東、云南等省（區(qū)）均有分布，主要集中分布在廣西的東部和南部（申文輝等，2014）。紅錐生長迅速，材質(zhì)良好，是優(yōu)質(zhì)的家具、船只、農(nóng)業(yè)用具材料，也是良好的造紙原料，果實富含淀粉，可做飼料（龔循勝等，2022），具有較高的經(jīng)濟價值和生態(tài)價值。盡管前人針對紅錐生態(tài)和生理適應(yīng)性開展了一些研究（梁艷紅和滕維超，2021；李娜等，2023；歐陽子龍等，2024），但不同地理種源紅錐葉片功能性狀在多大程度上受地理環(huán)境因子選擇作用，目前尚不十分清楚。鑒于此，本研究以廣西區(qū)內(nèi)6個代表性紅錐種群的115個自由授粉子代家系所營建的同質(zhì)園為研究對象，采用方差分析、主成分分析、聚類分析等方法，擬探討以下問題：（1）同質(zhì)園不同紅錐種群葉片功能性狀的變異、分化情況；（2）主要地理環(huán)境因子對紅錐葉片功能性狀塑造作用以及其地理變異規(guī)律。

1材料與方法

1.1試驗地概況

試驗地位于廣西南寧市西鄉(xiāng)塘區(qū)廣西壯族自治區(qū)林業(yè)科學(xué)研究院（以下簡稱廣西林科院）的老虎嶺試驗林場（ E ），平均海拔1 8 0 m ，屬南亞熱帶季風(fēng)氣候，水熱充沛，年均溫 ，年平均降水量 ，平均相對濕度

8 0 % ，年無霜期 3 6 0 d 。地形為低丘陵，地勢平緩，土壤以紅壤（磚紅壤）為主， 值5＼～6，土壤肥力均勻，適宜紅錐等造林樹種生長發(fā)育（黃志玲等，2016）。

1.2試驗材料

本試驗材料來自2012年營建于廣西林科院老虎嶺試驗林場的紅錐自由授粉家系子代試驗林，建林材料來自廣西區(qū)內(nèi)6個代表性種群的115個單株上的自由授粉種子（家系）（表1），其中老虎嶺（Laohuling，LHL）種群的13個家系來自1982年營建于該試驗林場紅錐試驗林優(yōu)良單株子代，該試驗林為廣西博白縣、天等縣等地點的天然種群收集的優(yōu)樹混系種子育苗營建；其余5個種群[博白縣（BobaiCounty，BB）、岑溪市（CenxiCity，CX）、浦北縣（PubeiCounty，PB）、容縣（RongCounty，RX）、天等縣（TiandengCounty，TD）］的試驗材料均采自各種群地天然林單株。

樣本采集：葉片樣本采集工作于2023年7月進行，在試驗林中以家系為單位進行葉片采集，每個家系隨機選擇生長良好、無明顯病蟲害的2個植株，每株選取樹冠中部正常生長的15片健康成熟葉片（每家系共計采集30片），放入密封袋內(nèi)，寫好標(biāo)簽，置于有冰袋的保溫箱內(nèi)帶回實驗室測定相關(guān)指標(biāo)。

1.3葉片功能性狀測定

每家系隨機選取采集的15片葉片進行形態(tài)性狀（表2）測定：吸干葉表水分，使用 電子天平測量鮮重，電子游標(biāo)卡尺測量葉厚；使用葉面積掃描儀（美國CID公司，型號：LI-3000C）掃描葉片形態(tài)，利用ImageJ軟件測量葉長、葉寬、葉面積等指標(biāo)；分別將每家系的15片葉片單獨用信封包裹起來置于烘箱內(nèi) 烘干至恒重，使用1/10000電子天平測量干質(zhì)重。基于上述數(shù)據(jù)計算葉形指數(shù)（葉長/葉寬）和比葉面積（葉面積/葉干重）。

生理生化性狀測定：參照《植物生理學(xué)實驗技術(shù)》（張治安和陳展宇，2008）及時測量每家系采取的另外15片葉片的總?cè)~綠素、可溶性糖、可溶性蛋白等重要生理生化指標(biāo)。具體測定方法：采用乙醇浸泡比色法測定總?cè)~綠素含量，硫代巴比妥酸法測定葉片丙二醛含量，酸性芘三酮法測定脯氨酸，蒽酮比色法測定可溶性糖，考馬斯亮藍(lán)G-250（CoomassiebrilliantblueG-25O）法測定可溶性蛋白。葉片功能性狀的單位及縮寫見表2。

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

性狀變異分析：利用統(tǒng)計軟件R4.2.2（RCoreTeam，2023）的aov函數(shù)進行巢式方差分析，比較種群間和種群內(nèi)葉片功能性狀指標(biāo)的差異，線性模型： 。式中， 為第 i 個種群第 j 個家系第 k 個觀測值， 為總平均值， 為種源效應(yīng)值， 為種群內(nèi)家系隨機效應(yīng)值， 為隨機誤差（Logan，2010）。采用性狀分化系數(shù)（traitdifferentiation coefficient， ）來評估紅錐種源間性狀分化水平，計算公式： 0 式中， 為種群間方差分量， 為種群內(nèi)（家系）方差分量（葛頌等，1988）。利用agricolae1.37數(shù)據(jù)包（deMendiburu，2023）進行Duncan's多重比較。

性狀和地理環(huán)境因子的關(guān)聯(lián)分析：基于采樣區(qū)域及其周邊氣象站觀測數(shù)據(jù)，結(jié)合相關(guān)文獻(xiàn)記載（張鳳良等，2011；李娜等，2023），篩選出海拔、經(jīng)度、緯度、年均溫和年降水量5個地理環(huán)境因子，對參試性狀進行主成分分析，獲得綜合反映葉片功能性狀的主成分，并篩選能解釋紅錐葉片變異的主要性狀?；赑earson相關(guān)分析探究13個葉片功能性狀與環(huán)境因子間的相關(guān)性?；趨⒃囆誀顚ΨN群進行層次聚類分析（Ward.D2法），并利用Mantel檢驗探究種群葉功能性狀地理變異模式。所有分析均采用R4.2.2及其相關(guān)擴展包執(zhí)行（Kassambaraamp; Mundt，2020；Weiamp; Simko，2021；Harrell，2024）。

2 結(jié)果與分析

2.1紅錐葉片功能性狀變異分析

參試紅錐種群葉片功能性狀在種群間和種群內(nèi)（家系間）2個層次的方差分析結(jié)果（表3）表明，參試的13個葉片功能性狀在種群間和種群內(nèi)存在極顯著性差異（ Plt;0 . 0 0 1 ）。葉片功能性狀的變異系數(shù)結(jié)果（表4）顯示，LA性狀的平均變異系數(shù)最大（ 4 6 . 7 1 % ），MDA性狀次之（ 3 5 . 4 6 % ，SP性狀最小 （ 7 . 1 9 % ）。此外，參試種群間的形態(tài)性狀（ 1 5 . 8 8 % ～ 4 6 . 7 1 % 的變異系數(shù)高于生理生化性狀 （ 7 . 1 9 % ～ 3 5 . 4 6 % 。綜合分析葉片功能性狀變異情況發(fā)現(xiàn)，紅錐葉片功能性狀的變異具有顯著的種群差異，形態(tài)性狀的種群分化高于生理生化性狀。

表3參試紅錐種群間和種群內(nèi)的葉片功能性狀巢氏方差分析

紅錐葉片功能性狀的多重比較分析（表5）顯示，13個葉片功能性狀在種群間存在顯著性差異。在形態(tài)性狀方面，CX種群的LL、LW、FW、DW、LA、SLA和LI等性狀均較大，其次為RX種群，而LT性狀的最大值出現(xiàn)在BB種群中。在生理生化性狀方面，RX種群MDA和SP2個指標(biāo)高于其他種群，而TC、Pro和SS3個生理性狀最大值分別在LHL、TD和PB3個種群內(nèi)。綜上表明，紅錐葉形態(tài)與生理生化性狀在群體間存在分化差異；參試種群中，CX種群與其他種群在形態(tài)性狀上的差異更加明顯。

表5紅錐種群葉片功能性狀的多重比較分析

2.2紅錐種群間的葉片功能性狀分化

13個紅錐葉片功能性狀的方差分量和性狀分化系數(shù) 結(jié)果（表6）顯示，紅錐葉片的性狀分化系數(shù)介于 1 0 . 9 5 % （SS）至 4 6 . 9 3 % （LI）之間，表明不同性狀在種群間和種群內(nèi)的方差分量占比存在顯著性差異。在葉形態(tài)性狀中，種群間的分化尤為顯著，其中LW、LT、DW、LA和LI等性狀的分化系數(shù)顯著高于均值（ 2 2 . 2 8 % ）。然而，在生理生化性狀中，僅有SP（ 2 2 . 5 9 % ）性狀的分化系數(shù)略高于均值。整體而言，紅錐葉生理生化性狀的分化系數(shù)普遍低于形態(tài)性狀。

2.3紅錐葉功能性狀和地理環(huán)境因子的主成分分析

葉片功能性狀與地理環(huán)境因子的主成分分析結(jié)果（圖1）顯示，主成分1的貢獻(xiàn)率為 4 0 . 6 3 % ，主成分2的為 2 6 . 8 4 % 。葉形態(tài)性狀與主成分1呈現(xiàn)較高的正相關(guān)關(guān)系，其中LL、FW、LW和LA等形態(tài)性狀分布集中，并且彼此間呈顯著性正相關(guān)；結(jié)合Pearson相關(guān)分析（圖2）可知，葉片形態(tài)性狀（LT除外）主要受海拔影響，隨海拔增加而增加，經(jīng)緯度的影響次之；LT則主要受年均溫的調(diào)控。與形態(tài)性狀相比，生理生化性狀分布較為分散，從向量夾角來看，Pro、TC與MDA、SP相關(guān)性不顯著，與SS呈負(fù)相關(guān);結(jié)合Pearson相關(guān)分析（圖2）可知，MDA、SP和SS 等生理性狀隨年降水量及經(jīng)緯度增加而增加，但與海拔呈負(fù)相關(guān)；Pro和TC的變化趨勢與其他生理性狀相反。綜上可知，紅錐葉片功能性狀主要受海拔、經(jīng)緯度和年降水量的協(xié)同調(diào)控。

表6紅錐葉片功能性狀的方差分量和分化系數(shù)

Table 6Variance components and differentiation coefficients of leaf functional traits of Castanopsis hystrix

2.4聚類分析與Mantel檢驗

基于13個葉片功能性狀對6個紅錐種群進行聚類分析（圖3），發(fā)現(xiàn)當(dāng)候選聚類數(shù)為2時，平均輪廓系數(shù)最大。以最佳聚類數(shù)2為基準(zhǔn)，當(dāng)表型距離大于3.45時可將6個種群分為2類，其中第1類包括博白縣（BB）、浦北縣（PB）、容縣（RX）、天等縣（TD）和老虎嶺（LHL）5個種群，而第2類僅有岑溪市（CX）1個種群。

為進一步分析紅錐種群葉片功能性狀變異與地理距離、海拔之間的關(guān)系，分別將6個種群葉片功能性狀的歐式距離與地理距離、海拔差進行Mantel檢驗，結(jié)果（圖4）顯示，功能性狀與地理距離（ 0 . 1 3 5 ， P = 0 . 0 2 8 ）、海拔差（ 均具有顯著相關(guān)性，表明紅錐種群的葉片功能性狀變異受地理隔離和海拔等因素的綜合影響。

3 討論與結(jié)論

3.1紅錐葉片功能性狀變異分析

不同種群間葉片性狀分化系數(shù)（ ）的差異可以反映葉片性狀的變異來源以及植物對種源地原生境的適應(yīng)能力（李洪果等，2019；Gomoryetal.，2023）。本研究結(jié)果表明，紅錐葉片形態(tài)性狀在種群間和種群內(nèi)均存在極顯著性差異（ Plt;0 . 0 0 1 ），種群間變異占比為 2 6 . 7 8 % （ ），略高于李娜等（2023）對17個種群紅錐的229個無性系葉片形態(tài)性狀研究結(jié)果（ ），但低于板栗（Castnanea mollisima， ）（劉亞斌和郭素娟，2020）、蒙古櫟（Quercusmongolica， 6 3 . 1 6 % ）（尹杰，2023）等其他殼斗科植物的種群分化。這可能與樹種特性、種群分布范圍等有關(guān)，種群地理來源越廣泛（生境差異越大），葉片性狀變異程度越大（ X u etal.，2023）。本試驗材料主要來源于廣西東南部，地理分布較上述2種殼斗科植物更為狹窄，種群間的基因交流可能更加頻繁，在一定程度上降低了種群間的遺傳分化；同時，在植物生活周期中幼苗階段對外界環(huán)境的變化最為敏感（Fieldsetal.，2017；劉從等，2018），上述蒙古櫟的試驗材料均來自林木的幼苗期（2年生幼苗），而本研究的試驗材料來自11年生的紅錐人工林，葉片的形態(tài)性狀可能更趨于穩(wěn)定。

Alt.海拔；Lon.經(jīng)度；Lat.緯度；MAT.年均溫；MAP.年均降水量；LL.葉長；LW.葉寬；LT.葉厚；FW.鮮重；DW.干重； LA.葉面積；SLA.比葉面積；LI.葉形指數(shù)；TC.總?cè)~綠素；MDA.丙二醛；Pro.脯氨酸；SS.可溶性糖；SP.可溶性蛋白。 下同。 Alt.Altitude;LonLongitude;Lat.Latitude;MAT.Meananualtemperature;MAP.Meanaualprecipitation;LL.Leaflength;LW.Leaf width;LT.Leafthickness；FW.Fresh weight；DW.Dryweight；LA.Leafarea；SLA.Specficleafarea；LILeafindex；TC.Total chlorophyll；MDA.Malondialdehyde;Pro.Proline；SS.Soluble sugar；SP.Soluble protein.The same below.

Fig.1Principal componentanalysisof leaffunctional traitsofCastanopsis hystrixand itsgeographicalenvironmentfact

表示 P lt; 0 . 0 5 ，***表示 Plt;0 . 0 1 。（204號 indicates P lt; 0 . 0 5 ，and indicates Plt;0 . 0 1

Fig.2Pearson correlations between leaf functional traits of Castanopsis hystrix geographical environment factors

葉片的生理生化性狀可以反映植物對逆境的適應(yīng)性（Tuveretal.，2022）。本研究結(jié)果顯示，紅錐不同種群間葉片生理生化性狀均存在極顯著性差異（ P lt; 0 . 0 0 1 ），不同種群間的變異占比為1 5 . 0 6 % ，與殼斗科栗屬（Castanea）雜交F1代葉片生理生化性狀的變異分析結(jié)果一致（章平生等，2021）。值得注意的是，盡管紅錐葉片形態(tài)性狀和生理生化性狀在種群間存在顯著性差異，但生理生化性狀種群間的變異占比低于形態(tài)性狀，與北京百花山23種落葉木本植物（高永龍等，2024）、中國北部小葉楊（Populussimonii）（李炳虹等，2020）檫木（Sassafras tzumu）（蔣艾平等，2016）等研究結(jié)果一致。推測其原因可能與紅錐的葉片形態(tài)性狀和生理生化性狀的遺傳穩(wěn)定性和功能差異有關(guān)。葉形性狀決定了植物在特定生境下的資源獲取效率，可能更多受到長期生境選擇作用，遺傳穩(wěn)定較高，而生理生化性狀（如內(nèi)二醛、脯氨酸、可溶性糖等與脅迫相關(guān)的生理生化指標(biāo)）可能對短期環(huán)境波動更為敏感（謝福春等，2008；Xieetal.，2014）。本研究中，老虎嶺林場的地理環(huán)境條件相對一致，其生境適宜紅錐林木的生長發(fā)育（黃志玲等，2016），使得參試紅錐種源受到相對較小的環(huán)境脅迫，生理生化響應(yīng)趨于一致。葉形態(tài)性狀因受原生境的長期選擇所用，具有較高的遺傳穩(wěn)定性，使得同質(zhì)園條件下種源間仍維持著較高分化（Dareusetal.，2021）。

研究表明，葉片的形態(tài)性狀（葉面積、比葉面積等）和生理生化性狀（葉綠素含量、可溶性糖等）等功能性狀與植物對外界資源的獲取、利用和保存息息相關(guān)（Reich，2014；Yangetal.，2019），通過分析葉片性狀之間的差異和葉壽命的長短，可以推測植物葉片性狀的“投資-收益”策略（龐世龍等，2021）。紅錐為常綠闊葉林樹種，受其生命周期的影響，葉片具有較長的壽命（Wangetal.，2023）；本研究中，紅錐具有較小的葉面積（ ）、比葉面積（ 和中等的葉綠素（ ）、可溶性糖（57.752 ）含量，說明參試紅錐種群在同質(zhì)園條件下具有較低的光合速率?；谏鲜鼋Y(jié)果推測，紅錐葉片在“投資-收益\"權(quán)衡中更趨向于選擇壽命長、比葉面積小、光合速率小的“緩慢投資-收益”型（保守型）策略（范宏坤等，2024），與同屬植物米（Castanopsiscarlesii）（張耀藝等，2021）的研究結(jié)論一致。

3.2葉片功能性狀的地理變異規(guī)律

在植物生長發(fā)育過程中，受海拔、經(jīng)緯度、年降水量等地理環(huán)境因子的驅(qū)動，其葉片不同功能性狀間協(xié)同作用，逐漸形成了特定的生態(tài)適應(yīng)策略（程雯等，2019；Lietal.，2021）。本研究中，除LT以外，形態(tài)性狀之間具有顯著正相關(guān)性，這意味著參試紅錐形態(tài)性狀在生長和適應(yīng)過程中具有較高的協(xié)同性，與蒙古櫟（尹杰，2023）、米櫧（仲小茹等，2022）、麻櫟（Quercusacutissima）（張元燕等，2014）等殼斗科植物的研究結(jié)果類似。結(jié)合Pearson相關(guān)分析結(jié)果，LL、LW、FW、LA等性狀與海拔和經(jīng)緯度呈正相關(guān)，其原因可能是LL、LW和LA等形態(tài)性狀與植物光照輻射吸收能力和儲存生物量有關(guān)（趙夏緯等，2019；蔡繼醇等，2024）。隨著海拔、緯度的升高，紅錐的LL、LW和LA等形態(tài)性狀不斷增長，在一定程度上可以增強紅錐林木的光合效率，以適應(yīng)溫度的降低和維持自身生長發(fā)育的需求，提高在林分中的競爭力。在生理生化性狀中，MDA、SS和SP等生理生化性狀與海拔呈負(fù)相關(guān)，與年降水量呈正相關(guān)，這可能與MDA、SP和SS的生理特性和紅錐的生長習(xí)性有關(guān)。研究表明，滲透調(diào)節(jié)是植物響應(yīng)外界環(huán)境變化的重要生理機制（Kouetal.，2023），MDA、SP和SS是植物生長發(fā)育的重要調(diào)節(jié)因子，MDA含量可以反映植物細(xì)胞膜脂的過氧化程度和植物對逆境條件反應(yīng)的強弱（范蘇魯?shù)龋?011），SS和SP在植物體內(nèi)相互作用，共同維持植物的滲透壓平衡（Liuetal.，2023）。紅錐主要分布在我國南部地區(qū)，喜溫暖濕潤的生態(tài)環(huán)境（洪偉等，2001），其適宜生長溫度約為 （劉春苑，2022）。因此，溫度對紅錐林木的生長發(fā)育至關(guān)重要（丘小軍等，2006），而本研究樣品于7月采集，月均溫大于28C （賈艷紅等，2023），基于此推測，低海拔地區(qū)的紅錐葉片生理生化性狀（MDA、SP和SS等）在一定程度上受到夏季溫度脅迫的影響，葉片生理生化性狀含量有所增加，以維持細(xì)胞內(nèi)外滲透壓的平衡，保障紅錐林木的生長發(fā)育。

基于13個功能性狀將6個紅錐種群分為2類，其中地理距離較遠(yuǎn)的TD種群與BB和PB等種群聚為1類，而平均海拔最高的CX（ 5 6 2 m 種群單獨為1類，說明海拔造成的生境異質(zhì)性可能對紅錐種群生態(tài)適應(yīng)性有重要影響。盡管Mantel檢驗結(jié)果顯示，參試紅錐種群葉片性狀的表型距離與地理距離、海拔差均具有顯著關(guān)系，符合距離隔離（isolationbydistance，IBD）模式，但值得注意的是，紅錐葉片性狀的表型距離與海拔差（ P =

0.013）的相關(guān)性高于地理距離（ ，說明其受海拔的影響更大。因此，開展紅錐種質(zhì)資源收集時，除了考慮種群間的地理隔離外，還需特別關(guān)注不同海拔梯度下的自然種群。綜上所述，本研究中紅錐葉片功能性狀呈現(xiàn)“海拔-經(jīng)緯度-降水”三者相互作用的地理變異模式。

3.3結(jié)論

紅錐種群間葉功能性狀變異豐富，變異主要來源于種群內(nèi)；基于主成分分析發(fā)現(xiàn)，紅錐葉片功能性狀間的相關(guān)性關(guān)系比較復(fù)雜，形態(tài)性狀之間呈正相關(guān)，與生理生化性狀間呈負(fù)相關(guān)；綜合地理環(huán)境因子和葉片功能性狀可知，紅錐葉片功能性狀的變異受海拔、經(jīng)緯度、降水等地理環(huán)境因子的綜合驅(qū)動。因此，本研究通過研究紅錐葉片功能性狀特征及其與地理環(huán)境因子的關(guān)系，有助于了解不同種群間紅錐葉片功能性狀的地理格局及其種內(nèi)變異的環(huán)境驅(qū)動力，為紅錐種質(zhì)資源的引種馴化、選育工作及開展其遺傳改良等提供了理論指導(dǎo)。

參考文獻(xiàn)：

CAI JC，LI ZJ，ZHOU GY，et al.，2024.Effects and actionmechanisms of increased light availability to understoryCinnamomum camphora seedlings[J].Journal of CentralSouthUniversityofForestryamp;Technology，44（7）：1-11.［蔡繼醇，李兆佳，周光益，等，2024.增加光照對林下樟樹苗木的影響及作用機理［J].中南林業(yè)科技大學(xué)學(xué)報，44（7）：1-11.]

CAODM，ZHANG YH，CHENG XQ，et al.，2021.Geneticvariation ofleaf phenotypic traitsindifferentpopulationsofPopuluscathayana[J].ScientiaSilvaeSinicae，57（8）：56-67.［曹德美，張亞紅，成星奇，等，2021.青楊不同種群葉片表型性狀的遺傳變異［J].林業(yè)科學(xué)，57（8）：56-67.]

CHENLJ，LIUMQ，LIAOBY，etal.，2O16.Variationofleafmorphological charactersofMeliaazedarach from differentprovenances[J]. Journal of Southwest Forestry University，36（5）：16-20.[陳麗君，劉明騫，廖柏勇，等，2016.苦楝不同種源葉片性狀變異研究［J］.西南林業(yè)大學(xué)學(xué)報，36（5） ：16-20.]

CHENGW，YUYH，XIONGKN，etal.，2O19.Leaf functionaltraitsof dominant species in karst plateau-canyon areas[J].Guihaia，39（8）：1039-1049.[程雯，喻陽華，熊康寧，等，2019.喀斯特高原峽谷優(yōu)勢種葉片功能性狀分析[J].廣西植物，39（8）：1039-1049.]

CORNELISSEN JHC，LAVOREL S，GARNIER E，et al.，2003.A handbook of protocols for standardised and easymeasurement ofplant functional traitsworldwide [J].Australian Journal of Botany，51（4） ：335-380.

DAREUS R，ACHARYA JP，PAUDEL DR，et al.，2021.Phenotypic diversity for phenological and agronomic traits inthe UC-Riverside cowpea （Vigna unguiculata L. Walp） mini-core collection [J]. Crop Science，61（5）：3551-3563.

DE MENDIBURU F， 2023.Agricolae： statistical procedures foragricultural research.R package version 1.3- 7[EB/OL].[2024 - 04 - 30]. https：//CRAN.R - project. org/package τ= agricolae.

DUQ，XU B，GONG C，et al.，2014.Variation in growth，leaf，and wood property traits of Chinese white poplar（Populus tomentosa），a major industrial tree species inNorthern China [J]. Canadian Journal of Forest Research，44（4）：326-339.

FAN HK，ZENG T，JINGZ，et al.，2O24.Leaf trait variationand trade-offs among growth types of broadleaf plants in XiaoHinggan Mountains [J].Chinese Journal of Plant Ecolgy，48（3）：364-376.［范宏坤，曾濤，金光澤，等，2024.小興安嶺不同生長型闊葉植物葉性狀變異及權(quán)衡［J]．植物生態(tài)學(xué)報，48（3）：364-376.]

FAN SL，YUAN ZH，F(xiàn)ENG LJ，et al.，2011 Effects of droughtstress on physiological and biochemical parameters of Dahliapinnata[J]. Journal of Applied Ecology，22（3）：651-657.［范蘇魯，苑兆和，馮立娟，等，2011.干旱脅迫對大麗花生理生化指標(biāo)的影響［J].應(yīng)用生態(tài)學(xué)報，22（3）：651-657.]

FIELDS WR，GRANT EHC，LOWE WH，2017. Detectingspatial ontogenetic niche shiftsin complex dendriticecological networks[J].Ecosphere，8（2）：e01662.

GAO YL，SUN YL，XU MZ，et al.，2024.Variationcharacteristics in leaf functional traits of woody plants indeciduous broadleaved forest community in Baihua Mountainof Beijing[J]. Journal of Beijing Forestry University，46（4）：40-51.［高永龍，孫艷麗，徐銘澤，等，2024.北京百花山落葉闊葉林群落內(nèi)木本植物的葉片功能性狀變異特征［J].北京林業(yè)大學(xué)學(xué)報，46（4）：40-51.]

GARNIER E，CORTEZ J，BILLES G，et al.，2004. Plantfunctional markers capture ecosystem properties duringsecondary succession[J].Ecology，85（9）：2630-2637.

GE S，WANG MX，CHEN YW，1988.An analysisofpopulation genetic structure of masson pine by isozymetechnique[J]. Scientia Silvae Sinicae，24（4）：399 -409.［葛頌，王明麻，陳岳武，1988.用同工酶研究馬尾松群體的遺傳結(jié)構(gòu)［J].林業(yè)科學(xué)，24（4）：399-409.]

GOMORYD，DITMAROVA L，HRIVNAK M，et al.，2023.Inconsistent phenotypic differentiation at physiological traitsin Norway spruce （Picea abies Karst.） provenances undercontrasting water regimes[J].Central European ForestryJournal，69（4）：214-223.

GONG XS，WANG YN，LIU LT，et al.，2022. Analysis on thedevelopment status of Castanopsis hystrix based on Cite Space[J].South China Forestry Science，50（3）： 41-46.[龔循勝，汪雁楠，劉麗婷，等，2022.基于Cite Space 的紅錐發(fā)展現(xiàn)狀分析［J].南方林業(yè)科學(xué)，50（3）：41-46.]

HARRELL JR FE，2024.Hmisc：Harrell misecellaneous.Rpackage version 5.1-3 [EB/OL].[2024-04-30]. https：//CRAN.R-project.org/package τ= τ Hmisc.

HONG W，LIU J，WU CZ，2001. A study on structure andspatlal distribution pattern of Castanopsis hystrix population[J].Scientia Silvae Sinicae，37（Sp.1）：6-10.[洪偉，柳江，吳承禎，2001.紅錐種群結(jié)構(gòu)和空間分布格局的研究[J].林業(yè)科學(xué)，37（Sp.1）：6-10.]

HUANG XL， ZHANG MY， ZHU LQ，et al.，2018. Comparisonofmorphological traits and physiological and biochemicalcharacteristics of leaves of four species in Podocarpus[J].Acta Agriculturae Jiangxi，30（8）：12-15.[黃相玲，張明月，朱栗瓊，等，2018.4種羅漢松葉片形態(tài)性狀及生理生化特性比較［J].江西農(nóng)業(yè)學(xué)報，30（8）：12-15.]

HUANG ZL， SHEN WH， ZHU JY，et al.，2016. Early growthperformance of ten rare tree species in southern Guangxi[J]. Journal of West China Forestry Science，45（5）： 84-88.［黃志玲，申文輝，朱積余，等，2016.桂南地區(qū)10種珍貴樹種人工林早期生長表現(xiàn)［J]．西部林業(yè)科學(xué)，45（5） ： 84-88.]

JIA YH，SU XQ，LI HR，2023. Characteristic analysis ofclimate change in Nanning from 1991 to 2020[J]. Journal ofGuangxi Normal University（Natural Science Edition），41（2）：190-200.［賈艷紅，蘇筱茜，李鶴冉，2023.1991—2020 年南寧市氣候變化特征分析［J].廣西師范大學(xué)學(xué)報（自然科學(xué)版），41（2）：190-200.]

JIANG AP， JIANG JM，LIU J， 2016. Responses of leaf traits ofSassafras tsumu（Hemsl.）Hemsl．a(chǎn)long an altitudinalgradient[J]. Chinese Journal of Ecology，35（6）：1467-1474.［蔣艾平，姜景民，劉軍，2016.檫木葉片性狀沿海拔梯度的響應(yīng)特征［J].生態(tài)學(xué)雜志，35（6）：1467-1474.]

KASSAMBARAA，MUNDTF，2O2O.Factoextra：Extract andvisualize the results of multivariate data analyses.R packageversion 1.0-7[EB/OL].[2024-04－30]. https：//CRAN.R-project.org/package=factoextra.

KOU J，YAN DH，QIN BT，et al.，2023.Physiologicalresponse mechanism of European birch（Betula pendulaRoth）to PEG-induced drought stress and hydration [J].Frontiers in Plant Science，14：1226456.

LI BH，SONG YP，ZHANG DQ，et al.，2O2O.Variationanalysis of photosynthetic physiology and leaf phenotype onPopulus simonii in northern China [J]. Journal of FujianAgriculture andForestry University（Natural ScienceEdition），49（1）：45-50.[李炳虹，宋躍朋，張德強，等，2020.中國北部小葉楊光合生理與葉片表型變異分析[J].福建農(nóng)林大學(xué)學(xué)報（自然科學(xué)版），49（1）：45-50.]

LI HG，CHEN DZ，XU JS，et al.，2019. Phenotypic diversityand variation in natural populations of Erythrophleum fordi，an endangered plant species [J]. Scientia Silvae Sinicae，55（4）：69-83.［李洪果，陳達(dá)鎮(zhèn)，許靖詩，等，2019.瀕危植物格木天然種群的表型多樣性及變異［J]．林業(yè)科學(xué)，55（4） ： 69-83.]

LI N，YANG YM，LI ZY，et al.，2O23.Variation analysis ofleaf phenotypic in different provenance of Castanopsis hystrix[J].Subtropical Plant Science，52（4）：318-326.[李娜，楊袁木，李子玥，等，2023.不同種源紅錐葉片表型變異分析［J].亞熱帶植物科學(xué)，52（4）：318-326.]

LI Y，LI S，LU X，et al.，2021.Leaf phenotypic variation ofendangered plant Tetracentron sinense Oliv.and influence ofgeographical andclimatic factors[J]. Journalof ForestryResearch，32（2） ： 623-636.

LIANG YH， TENG WC，2021. Effect of manganese stress ongrowth and physiological indexes of Castanopsis hystrix youngseedling［J].Guangxi Forestry Science，50（6）：685 -689.［梁艷紅，滕維超，2021.錳脅迫對紅錐幼苗生長和生理指標(biāo)的影響[J].廣西林業(yè)科學(xué)，50（6）：685-689.]

LIU C，TIAN T，LI S，et al.，2O18. Growth response of Chinesewoodyplant seedlings to different light intensities［J].ActaEcologica Sinica，38（2）：518-527.[劉從，田甜，李珊，等，2018.中國木本植物幼苗生長對光照強度的響應(yīng)[J].生態(tài)學(xué)報，38（2）：518-527.]

LIUCY，2O22.Biological characteristicsandartificialcultivation technology of Castanopsis hystrix [J].AnhuiAgricultural Science Buletin，28（8）：97-98，112.[劉春苑，2022.紅錐生物學(xué)特性及其人工培育技術(shù)［J].安徽農(nóng)學(xué)通報，28（8）：97-98，112.]

LIU FH，ZHAO Y，WANG XR，et al.，2023.Physiologicalresponse and drought resistance evaluation ofGleditsiasinensis seedlings under drought-rehydration state[J].Scientific Reports，13（1） ：19963.

LIU LL，DU N，PEI CP，et al.，2017. Advances in populationepigenetics of plant [J]. Chinese Journal of Ecology，36（9）：2615-2622.［劉樂樂，杜寧，裴翠萍，等，2017.植物群體表觀遺傳學(xué)研究進展［J].生態(tài)學(xué)雜志，36（9）：2615-2622.]

LIU YB，GUO SJ，2020.Leaf functional traits of 39 chestnutgermplasms in Qianxiregion [J]. Journal of NortheastForestry University，48（4）：8-14.[劉亞斌，郭素娟，2020.遷西地區(qū)39份板栗種質(zhì)葉片的功能性狀多樣性[J]．東北林業(yè)大學(xué)學(xué)報，48（4）：8-14.]

LOGAN M， 2O1o. Biostatistical design and analysis using R：Apractical guide[M].West Sussex，UK：Wiley-Blackwell：283-312.

LUO XY，WANG Z， CHEN Y，et al.， 2024. Research progresson phenotypic variation of leaves and fruits among naturalpopulations of woody plants [J]. Contemporary Horticulture，47（5）：66-68.[羅欣語，王芷，陳飚，等，2024.木本植物天然種群間葉、果表型性狀變異研究進展［J]．現(xiàn)代園藝，47（5）： 66-68.]

MA M，ZHAO CZ，KANG MP，et al.，2023. Periploca sepiumleaf dry weight， area，and thickness scale diferently underdifrent lightenvironmentsin Taohe riparian forests[J].Flora，301：152255.

MENG YD，DU HY，WANG L，et al.，2O22.Diversity analysisof leaf phenotypic traits of Eucommia ulmoides germplasmresources[J].Forest Research，35（5）：103-112.[孟益德，杜紅巖，王璐，等，2022.杜仲種質(zhì)資源葉片表型性狀多樣性分析［J].林業(yè)科學(xué)研究，35（5）：103-112.]

OUYANG ZL，LEI K，JIA XL，et al.，2024. Effects of siliconand manganese on growth and physiology of Castanopsishystrix[J]. Chinese Journal of Tropical Agriculture，44（5）：41-49.[歐陽子龍，雷康，賈湘璐，等，2024.硅錳互作對紅錐幼苗生長生理的影響［J]．熱帶農(nóng)業(yè)科學(xué)，44（5）：41-49.]

PANG SL，OU ZY，LING FC，et al.，2021．Leaf traitvariations and leaf economic spectrum of 18 plant speciessuitable in a karst area of Southwest Guangxi [J]. ChineseJournal of Ecology，40（10）：3041-3049.［龐世龍，歐芷陽，凌福誠，等，2021.桂西南巖溶區(qū)18種適生植物葉性狀變異及經(jīng)濟譜［J].生態(tài)學(xué)雜志，40（10）：3041-3049.]

QIU XJ，ZHU JY，JIANG Y，et al.，2OO6. Study on naturaldistribution and adaptability of Castanopsis hystrix [J].Journal of Guangxi Agriculture and Biological Science，25（2）：175-179.［丘小軍，朱積余，蔣燚，等，2006.紅錐的天然分布與適生條件研究［J].廣西農(nóng)業(yè)生物科學(xué)，25（2）： 175-179.]

R CORE TEAM，2023. R：A language and environment forstatistical computing[EB/OL].[2024-04-30]. https：//www.mendeley. com/catalogue/c36f6d0e- 0456 - 3314 -

REICH PB，2014. The world-wide“fast-slow”plant economicsspectrum：A traits manifesto [J]. Journal of Ecology，102（2） ： 275-301.

SANDOVAL-GRANILLO V，MEAVE JA， 2O23.Leaf functionaldiversity and environmental filtering in a tropical dry forest：Comparison between two geological substrates[J]. Ecologyand Evolution，13（9）：e10491.

SHEN WH，ZHU JY，LIU X，et al.，2014. Regional testingand superior provenance selection of Castanopsis hystrixprovenances [J]. Journal of Central South University ofForestryamp; Technology，34（3）：11-17.［申文輝，朱積余，劉秀，等，2014.紅錐種源區(qū)域試驗與優(yōu)良種源選擇[J].中南林業(yè)科技大學(xué)學(xué)報，34（3）：11-17.]

SONG SS，WANG J，CAI CH，et al.，2023. Response andadaptation of leaf functional traits of an endangered lianaMonimopetalum chinense to environmental factors [J]. ActaEcologica Sinica，43（17）：7252-7262.[宋帥帥，王建，蔡朝暉，等，2023.瀕危植物永瓣藤葉片功能性狀對環(huán)境因子的響應(yīng)［J].生態(tài)學(xué)報，43（17）：7252-7262.]

SUNM，TIANK，ZHANGY，etal.，2O17.Research onleaffunctional traits and their environmental adaptation[J].Plant ScienceJournal，35（6）：940-949.［孫梅，田昆，張贊，等，2017.植物葉片功能性狀及其環(huán)境適應(yīng)研究[J].植物科學(xué)學(xué)報，35（6）：940-949.]

TANG QL，LIANG LH，DENG ZY，et al.，2O22.Variationanalysisonleafmorphologiccharactersofdifferentprovenance of Corymbia citriodora ssp. variegata [J].MolecularPlantBreeding，20（16）：5454-5462.[唐慶蘭，梁麗華，鄧紫宇，等，2022.不同種源斑皮檸檬桉葉片形態(tài)變異分析［J].分子植物育種，20（16）：5454-5462.]

TUVER GY， EKINCI M， YILDIRIM E， 2022. Morphological，physiologicalandbiochemical responsestocombinedcadmium and drought stress in radish （Raphanus sativus L.）[J].Rendiconti Lincei. Scienze Fisiche e Naturali， 33（2）：419-429.

WANG C，LU J， ZHOU CN，et al.，2021. Altitude distributionofleaf functional traitsof Quercusaquifolioidesinsoutheastern Tibet[J].Journal of Forest and Environment，41（4）：366-372.[王超，盧杰，周晨霓，等，2021.藏東南川滇高山櫟葉功能性狀海拔分布特征［J].森林與環(huán)境學(xué)報，41（4）：366-372.]

WANG H，PRENTICE IC，WRIGHT IJ，et al.，2023.Leafeconomics fundamentals explained by optimality principles[J]. Science Advances，9（3）： eadd5667.

WEI TY，SIMKOV，2O21.R package“corrplot”：visualizationof a Correlation Matrix （Version O.92）[EB/OL].[2024-04-30]. https：//github.com/taiyun/corrplot.

XIE FC， ZHANG WT，LIU FQ，et al.， 2008.Effect of soil NaClstresson physiological and biochemical characteristicsofClerodendrum trichotonum Thumb [J].Acta AgriculturaeUniversitatis Jiangxiensis，30（5）：839-844.［謝福春，張文婷，劉富強，等，2008.土壤鹽脅迫對海州常山生理生化特性的影響［J].江西農(nóng)業(yè)大學(xué)學(xué)報，30（5）：839-844.]

XIE LN，MA CC，GUO HY，et al.，2014. Distribution patternof Caragana species under the influence of climate gradientin the Inner Mongolia region，China[J]. Journal of AridLand，6（3） ： 311-323.

XIONG L，LONG CL，LIANG S，et al.，2O24. Response of leaffunctional traits of woody plants to soil characteristics inkarst forests [J].Journal of Tropical and SubtropicalBotany，32（3）：310-318.［熊玲，龍翠玲，梁盛，等，2024.喀斯特森林木本植物葉片功能性狀對土壤特性的響應(yīng)［J].熱帶亞熱帶植物學(xué)報，32（3）：310-318.]

XU TY，NIU X，WANGB，et al.，2O23.Variations in leaffunctionaltraitsandphotosyntheticparametersofCunninghamia lanceolata provenances[J]. Forests，14（9） ： 1708.

YANGY，DOUY，CHENGH，et al.，2019.Plant functionaldiversity drives carbon storage following vegetation restorationinLoess Plateau，China [J].Journal of EnvironmentalManagement，246： 668-678.

YIN J，2O23.Analysis of leaf phenotypic character variation ofQuercus mongolica seedlings from different provenances[J].Forest Science and Technology（8）：56-60.[尹杰，2023.不同種源蒙古櫟幼苗葉表型性狀變異分析［J].林業(yè)科技通訊（8）：56-60.]

ZHANG FL，ZHANG FQ，PAN W，et al.，2O11.Variationanalysis on leaf traits of 17 Castanopsis hystix provenances[J].Guangdong Forestry Science and Technology，27（3）：20-26.[張鳳良，張方秋，潘文，等，2011.17個紅錐種源葉片性狀變異分析[J].廣東林業(yè)科技，27（3）：20-26.]

ZHANG PS，JIANG XB，GONG BC，et al.，2O21.Leafphenotype variation and heterosis in F1 progeny of crossbetween Castanea mollissima and C.henryi[J].BulletinofBotanical Research，41（2）：281-293.［章平生，江錫兵，龔榜初，等，2021.板栗與錐栗雜交F1代葉片表型變異及雜種優(yōu)勢研究［J].植物研究，41（2）：281-293.]

ZHANG YY，NI XY，YANG J，et al.，2O21. Nitrogen andphosphorus resorption and stoichiometric characteristics ofdiffrent tree species in a mid-subtropical common-garden，China[J]. Chinese Journal of Applied Ecology，32（4）：1154-1162.［張耀藝，倪祥銀，楊靜，等，2021.中亞熱帶同質(zhì)園不同樹種氮磷重吸收及化學(xué)計量特征［J].應(yīng)用生態(tài)學(xué)報，32（4）：1154-1162.]

ZHANG YY，YU MK，F(xiàn)ANG YM，2014.Analysisonphenotypic trait variation of different provenances of Quercusacutissima [J]. Journal of Plant Resources and Environment，23（3）：36-44.[張元燕，虞木奎，方炎明，2014.麻櫟不同種源的表型性狀變異分析［J].植物資源與環(huán)境學(xué)報，23（3）：36-44.]

ZHANG ZA，CHEN ZY，2008. Experimental techniques ofplantphysiology[M]//CHEN ZY.Changchun：JilinUniversityPress：178-193.［張治安，陳展宇，2008.植物生理學(xué)實驗技術(shù)［M]//陳展宇.長春：吉林大學(xué)出版社：178-193.]

ZHAO XW，WANG YF，MA WM，2019.The relationship oftranspirationratewithleafcharactersofThermopsis lanceolatein different slopes of analpine meadow[J].Acta EcologicaSinica，39（7）：2494-2500.［趙夏緯，王一峰，馬文梅，2019.高寒草地不同坡向披針葉黃華蒸騰速率與葉性狀的關(guān)系［J].生態(tài)學(xué)報，39（7）：2494-2500.]

ZHONG XR， ZHANG L， PAN XH， et al.， 2022. Response andadaptation of leaf functional traits to different altitudesinevergreen broad-leaved forest of Castanopsis carlesii[J].Acta Agriculturae Universitatis Jiangxiensis，44（6）：1438-1447.［仲小茹，張露，潘昕昊，等，2022.常綠闊葉林米櫧葉片功能性狀對不同海拔梯度的響應(yīng)與適應(yīng)［J].江西農(nóng)業(yè)大學(xué)學(xué)報，44（6）：1438-1447.]

（責(zé)任編輯 李莉王登惠）