馬 麗， 單立山， 解婷婷， 種培芳， 楊 潔， 師亞婷

(甘肅農(nóng)業(yè)大學(xué)林學(xué)院， 甘肅 蘭州 730070)

植物功能性狀是指植物在個(gè)體水平上所表現(xiàn)出的植物形態(tài)、生理等特征，能夠有效表達(dá)植物在長(zhǎng)期進(jìn)化過(guò)程中對(duì)外界環(huán)境變化的響應(yīng)[1-2]。葉片作為植物地上部分的主要器官，是植物進(jìn)行光合作用和物質(zhì)生產(chǎn)的主要場(chǎng)所[3]，其功能性狀如葉片含水量、比葉面積、葉片碳氮磷含量和解剖結(jié)構(gòu)等最能體現(xiàn)植物對(duì)環(huán)境因子的適應(yīng)性進(jìn)化[4-6]。目前，植物功能性狀研究已成為植物功能生態(tài)學(xué)研究中的前沿和熱點(diǎn)問(wèn)題之一[7-8]。雖有大量研究報(bào)道了植物葉功能性狀變異及其與環(huán)境因子的關(guān)系[9-11]，但并不能說(shuō)明植物葉功能性狀是受環(huán)境影響還是遺傳因素的影響，而同質(zhì)園實(shí)驗(yàn)?zāi)軌蛟谝欢ǔ潭壬峡刂骗h(huán)境因素所造成的影響，進(jìn)而說(shuō)明不同種源植物葉功能性狀的變異是由于原種源地環(huán)境差異造成的植物遺傳性基因變異。因此，基于同質(zhì)園實(shí)驗(yàn)探究植物葉功能性狀的變異及其與環(huán)境因子的關(guān)系具有重要意義，有助于我們更好地理解不同種源地植物對(duì)環(huán)境梯度變化的響應(yīng)和適應(yīng)機(jī)制。

植物葉片通過(guò)改變自身的形態(tài)和化學(xué)元素來(lái)適應(yīng)環(huán)境變化[12-13]。目前，關(guān)于植物葉功能性狀受環(huán)境因子影響的研究主要集中在全球等大尺度層面[14]，但這種研究會(huì)忽略一些小尺度層面的環(huán)境因子如微地形、生境等，對(duì)這些環(huán)境因子的研究也許能夠揭示大尺度下所無(wú)法解釋的環(huán)境過(guò)濾效應(yīng)[15]。植物為適應(yīng)環(huán)境會(huì)形成以該環(huán)境為主的特定地理種源，遺傳變異就是植物為適應(yīng)該環(huán)境變化的適應(yīng)結(jié)果，而由于短期環(huán)境因子如土壤質(zhì)量等引起的植物葉功能性狀的變異是植物適應(yīng)短期環(huán)境變化的結(jié)果[16]。此外，在以往試驗(yàn)中沿自然環(huán)境梯度采樣時(shí)，很難分辨植物表現(xiàn)型、基因型和環(huán)境之間的關(guān)系，如黃河上游的寧夏平原和黃河下游的黃河三角洲地區(qū)的蘆葦(Phragmitesaustralis)植物葉功能性狀的差異同時(shí)受到環(huán)境和遺傳兩個(gè)因素的共同影響[17]；全緣葉綠絨蒿(Meconopsisintegrifolia)隨著海拔升高，株型變矮，葉片組織結(jié)構(gòu)緊密度增大[18]。然而，同質(zhì)園實(shí)驗(yàn)將不同生境的物種栽培在環(huán)境因素相對(duì)一致的同一園地內(nèi)，進(jìn)行統(tǒng)一栽培管理，使其生境保持一致；該方法能夠減少立地條件和氣候變化造成的影響，在局部適應(yīng)性研究中發(fā)揮重要作用[19]。近年來(lái)，雖然有關(guān)同質(zhì)園的生態(tài)學(xué)研究已有諸多的報(bào)道，但主要集中在森林植物葉功能性狀變異和對(duì)環(huán)境的適應(yīng)性研究等方面[20-21]，而關(guān)于荒漠植物在其葉功能性狀變異及影響機(jī)制方面的研究則較少。

白刺(Nitrariatangutorum)和紅砂(Reaumuriasoongarica)是半干旱、干旱地區(qū)常見(jiàn)的鹽生、旱生荒漠植物，具有抗旱、抗寒、耐高溫和耐鹽堿等特性以及防風(fēng)固沙、水土保持等作用[22-23]。白刺和紅砂廣泛分布于我國(guó)西北荒漠地區(qū)，并在當(dāng)?shù)匦纬闪艘运鼈優(yōu)榻ㄈ悍N或優(yōu)勢(shì)種的荒漠群落[23]。以往對(duì)白刺和紅砂的研究多集中在自然環(huán)境梯度下葉片形態(tài)解剖、根系的形態(tài)和功能、對(duì)極端環(huán)境的適應(yīng)性及群落功能等方面[24-27]，但是采用同質(zhì)園實(shí)驗(yàn)方法對(duì)不同種源兩種典型荒漠植物葉功能性狀的研究相對(duì)較少。因此，本研究分析了不同地理種源兩種典型荒漠植物葉片的肉質(zhì)化程度(Degree of fleshiness，DOF)、含水量(Water content，WC)、組織密度(Tissue density，TD)、比葉面積(Specific leaf area，SLA)、碳(Carbon，C)含量、氮(Nitrogen，N)含量、磷(Phosphorus，P)含量、碳氮比(Carbon and nitrogen ratio，C∶N)、碳磷比(Carbon and phosphorus，C∶P)和氮磷比(Nitrogen and phosphorus ratio，N∶P)指標(biāo)，旨在探討不同地理種源兩種典型荒漠植物葉功能性狀在種源間是否存在差異以及其與種源環(huán)境因子的關(guān)系，以期為兩種典型荒漠植物選擇優(yōu)良種源、適宜的種植環(huán)境及兩種典型荒漠植物對(duì)環(huán)境適應(yīng)策略等提供理論依據(jù)。

1 材料與方法

1.1 研究地概況

同質(zhì)園位于甘肅省武威市涼州區(qū)羊下壩鎮(zhèn)甘肅武威市良種繁育中心(38°24′ N，103°9′ E)，地處河西走廊北部，屬于典型的溫帶大陸性荒漠氣候，干旱少雨，土壤性質(zhì)為典型砂壤土，pH值為8.2，有機(jī)碳含量為14.30 g·kg-1，全氮含量為0.68 g·kg-1，全磷含量為1.42 g·kg-1，速效鉀含量為152.6 mg·kg-1，0～20 cm土壤容重為1.57 g·cm-3，平均土壤含水量0.4%，海拔為1 378 m，年均降雨量為113.2 mm，年均蒸發(fā)量為2 604.3 mm，年均氣溫為6.9℃，7月平均溫為24.9℃，1月平均氣溫為—9.9℃，日最高氣溫為35.9℃，日最低氣溫為—25.1℃[28-29]。

1.2 試驗(yàn)材料的培養(yǎng)及選取

2011年在白刺和紅砂的成熟季節(jié)，沿自然環(huán)境降水梯度在我國(guó)甘肅酒泉(JQ)、張掖(ZY)、武威(WW)和蘭州(LZ)天然分布區(qū)分別選取每個(gè)種源中生長(zhǎng)良好、無(wú)病蟲(chóng)害的多個(gè)優(yōu)勢(shì)單株，采集種子，并進(jìn)行標(biāo)記裝袋；2012年1月在武威市涼州區(qū)羊下壩鎮(zhèn)甘肅武威市良種繁育中心苗圃溫室進(jìn)行穴盤(pán)育苗，基質(zhì)成分為沙土10%～20%，砂壤土60%～80%，腐熟羊糞10%～20%，播種后進(jìn)行日常的苗木管理。2013年4月，土壤解凍后，進(jìn)行移栽，建立白刺和紅砂種源試驗(yàn)地。移栽時(shí)采用隨機(jī)區(qū)組試驗(yàn)設(shè)計(jì)，株行距0.5 m×1.5 m。移栽后隨即灌水，保證苗木的成活率。待苗木成活后，適時(shí)除草，按日常管理方法對(duì)苗木進(jìn)行管理[29]。

2020年在苗木生長(zhǎng)旺季，以同質(zhì)園內(nèi)定植的4個(gè)種源地的白刺和紅砂作為研究對(duì)象(4個(gè)種源地的地理位置與環(huán)境因子見(jiàn)表1[30])，分別選取每個(gè)種源中生長(zhǎng)健壯、長(zhǎng)勢(shì)良好、無(wú)病蟲(chóng)害的10個(gè)成年植株，剪取每個(gè)植株直立莖中部的葉片放入打濕的信封袋中，并且裝入封口塑料袋內(nèi)封口，立即帶回實(shí)驗(yàn)室用于葉片功能指標(biāo)測(cè)定。

表1 白刺和紅砂4個(gè)種源地的地理位置與環(huán)境因子Table 1 Geographical location and environmental factors of the four provenances of N. tangutorum and R. soongarica

1.3 測(cè)定方法

1.3.1葉片形態(tài)性狀的測(cè)定 將采集葉片上的雜質(zhì)清理掉，稱(chēng)量葉鮮重，再將葉片用根系分析軟件(Win-RhIZO 2008a)測(cè)量葉片表面積、體積等形態(tài)特征。測(cè)量完成之后，將葉片置于75℃烘箱烘至恒重，稱(chēng)量干重并記錄。每種源各有3次重復(fù)。白刺和紅砂葉片的肉質(zhì)化程度、含水量、比葉面積和組織密度的計(jì)算公式如下：

(1)

(2)

(3)

(4)

1.3.2葉片碳(C)、氮(N)和磷(P)含量的測(cè)定 將以上測(cè)定完的植物樣品用碾槽碾碎，裝進(jìn)密封袋用于進(jìn)行葉片C，N，P含量的測(cè)定。C含量測(cè)定采用重鉻酸鉀氧化加熱法；N含量測(cè)定采用凱氏定氮法；P含量測(cè)定采用釩鉬黃比色法測(cè)定[31]。

1.4 氣象數(shù)據(jù)來(lái)源

本研究選擇了對(duì)植物葉功能性狀具有影響的環(huán)境因子。氣象數(shù)據(jù)以酒泉、張掖、武威和蘭州4個(gè)種源地的數(shù)據(jù)為主(http://www.weather.com.cn/)。通過(guò)對(duì)數(shù)據(jù)的整理，計(jì)算出每個(gè)種源附近各站點(diǎn)的經(jīng)度(Longitude，Lng)、緯度(Latitude，Lat)、海拔(Altitude，ALT)、年均溫度(Mean annual temperature，MAT)、年均降水量(Mean annual precipitation，MAP)、1月均溫(January annual temperature，JaT)、7月均溫(July annual temperature，JuT)等數(shù)據(jù)，作為本次研究的氣象數(shù)據(jù)。

1.5 數(shù)據(jù)處理

本試驗(yàn)所有數(shù)據(jù)統(tǒng)計(jì)分析和相關(guān)作圖均是使用Excel 2016，SPSS 24.0(America)，Origin 2018(America)和Canoco5.0(America)軟件完成。采用單因素方差分析法分析同一物種在不同種源間葉片各性狀差異顯著性，并用Duncan檢驗(yàn)法和最小顯著差異法(LSD)進(jìn)行多重比較；用Pearson相關(guān)系數(shù)檢驗(yàn)白刺和紅砂各葉功能性狀之間的相關(guān)性。采用冗余分析(RDA)法和逐步線性回歸法分析白刺和紅砂葉功能性狀對(duì)不同種源地環(huán)境因子的響應(yīng)關(guān)系；在RDA分析之前，先對(duì)葉功能性狀進(jìn)行去趨勢(shì)對(duì)應(yīng)分析(DCA)，確定白刺和紅砂是采用以線性模型排序的冗余分析(RDA)還是以單峰模型排序的典型對(duì)應(yīng)分析(CCA)。若最大排序軸長(zhǎng)小于3，則適合線性模型排序分析的RDA方法；若最大排序軸大于4，則適合單峰模型排序的CCA方法；若最大排序軸長(zhǎng)在3～4之間，則兩種方法都可以。本研究中最大排序軸長(zhǎng)為0.1，小于3，故采用RDA方法分析，且將種源地環(huán)境因子作為解釋變量(空心箭頭)，將葉功能性狀作為響應(yīng)變量(實(shí)心箭頭)。文中數(shù)據(jù)均為平均值±標(biāo)準(zhǔn)誤差。

2 結(jié)果與分析

2.1 白刺和紅砂葉片形態(tài)性狀的種源間差異

由圖1可知，白刺葉片肉質(zhì)化程度在JQ種源與其他地理種源差異顯著(P<0.05)，且在JQ種源最小(3.99)，與ZY種源、WW種源和LZ種源相比，JQ種源白刺葉片肉質(zhì)化程度減小8.62%，14.36%和10.95%；含水量在JQ種源與其他地理種源差異顯著(P<0.05)，且在JQ種源最小(74.88)，與ZY種源、WW種源和LZ種源相比，JQ種源白刺葉片含水量減小2.56%，4.47%和3.22%；組織密度在種源間差異顯著(P<0.05)，且在ZY種源最大(0.08 g·cm-3)，在WW種源最小(0.06 g·cm-3)，與WW種源相比，ZY種源白刺葉片組織密度顯著增加了30.79%。

圖1 不同地理種源白刺和紅砂葉片形態(tài)結(jié)構(gòu)Fig.1 Leaf morphological structure of N. tangutorum and R. soongarica from different geographic provenances注：JQ-酒泉、ZY-張掖、WW-武威、LZ-蘭州；不同小寫(xiě)字母表示同一物種在不同種源間差異顯著(P<0.05)。下同Note:JQ-Jiuquan,ZY-Zhangye,WW-Wuwei,LZ-Lanzhou;different lowercase letters indicate significant difference between different provenances under the same species at the 0.05 level. The same as below

紅砂葉片的比葉面積在JQ和WW種源間差異顯著(P<0.05)，與WW種源相比，JQ種源的紅砂葉片的比葉面積增加了22.24%；肉質(zhì)化程度、含水量和組織密度在各個(gè)種源均無(wú)顯著差異。

2.2 白刺和紅砂葉片C，N，P含量及其計(jì)量比的種源間差異

由圖2可知，白刺葉片N含量在種源間差異顯著(P<0.05)，在JQ種源最大，為3.79 g·kg-1，在LZ種源最小，為3.31 g·kg-1；葉片N∶P在JQ種源與其他種源間均差異顯著(P<0.05)，且在JQ種源比值最大，為10.37，與ZY種源、WW種源和LZ種源相比，JQ種源葉片N∶P增加了47.95%，24.72%和63.72%。

圖2 不同地理種源白刺和紅砂葉片C，N，P含量及其計(jì)量比Fig.2 Contents of C,N,and P in leaves of N. tangutorum and R. soongarica and their stoichiometric ratio from different geographic provenances

紅砂葉片P含量在JQ和WW種源間差異顯著(P<0.05)，在JQ種源最大，為0.37 g·kg-1，在WW種源最小，為0.28 g·kg-1，與WW種源相比，JQ種源紅砂葉片P含量增加了30.40%。

二者的C含量、C∶N和C∶P在各個(gè)種源均無(wú)顯著差異。

2.3 白刺和紅砂葉片功能性狀間的相關(guān)性

通過(guò)Pearson相關(guān)性分析得出(表2)，白刺葉片的組織密度與肉質(zhì)化程度、含水量和比葉面積呈極顯著負(fù)相關(guān)關(guān)系(P<0.01)；N含量與比葉面積呈顯著負(fù)相關(guān)關(guān)系(P<0.05)，C∶P與C含量呈極顯著正相關(guān)關(guān)系(P<0.01)；N∶P與N含量呈極顯著正相關(guān)關(guān)系(P<0.01)，與C∶N呈顯著負(fù)相關(guān)關(guān)系(P<0.05)。通過(guò)Pearson相關(guān)性分析可以看出(表3)，紅砂葉片的組織密度與比葉面積呈極顯著負(fù)相關(guān)關(guān)系(P<0.01)；C∶N與C含量呈極顯著正相關(guān)關(guān)系(P<0.01)，與N含量呈顯著負(fù)相關(guān)關(guān)系(P<0.05)；C∶P與C含量呈極顯著正相關(guān)關(guān)系(P<0.01)，與P含量呈顯著負(fù)相關(guān)關(guān)系(P<0.05)，與C∶N呈極顯著正相關(guān)關(guān)系(P<0.01)；N∶P與N含量呈顯著正相關(guān)關(guān)系(P<0.05)。

表2 不同地理種源白刺葉片功能性狀間的相關(guān)性Table 2 Correlation among leaf functional traits of N. tangutorum from different geographic provenances

表3 不同地理種源紅砂葉片功能性狀間的相關(guān)性Table 3 Correlation among leaf functional traits of R. soongarica from different geographic provenances

2.4 不同地理種源白刺和紅砂葉片功能性狀與種源地環(huán)境因子的關(guān)系

2.4.1葉片功能性狀與種源地環(huán)境因子的冗余分析 通過(guò)RDA分析(表4)，第1，2軸的白刺葉功能性狀與環(huán)境的相關(guān)系數(shù)分別為0.9163，0.9213，特征值分別為0.3265，0.1574，共解釋了48.40%的葉功能性狀和92.37%的葉功能性狀與種源地環(huán)境的關(guān)系，表明前兩個(gè)排序軸能夠較好地反映白刺葉功能性狀與種源地環(huán)境因子的關(guān)系。蒙特卡羅檢驗(yàn)第1排序軸和所有排序軸均達(dá)到極顯著水平(P<0.01)，表明結(jié)果具有統(tǒng)計(jì)學(xué)意義。由圖3可知，海拔與肉質(zhì)化程度、含水量呈正相關(guān)關(guān)系，與組織密度呈負(fù)相關(guān)關(guān)系；緯度與N∶P呈正相關(guān)關(guān)系，與P含量呈負(fù)相關(guān)關(guān)系；年均降水量與C含量、N含量、C∶P呈負(fù)相關(guān)關(guān)系。RDA分析(表5)表明，第1，2軸的特征值分別為0.804 6，0.006 9，共解釋了81.16%的紅砂葉功能性狀和99.69%的葉功能性狀與種源地環(huán)境的關(guān)系，表明前兩個(gè)排序軸能夠較好地反映紅砂葉功能性狀與種源地環(huán)境因子的關(guān)系。蒙特卡羅檢驗(yàn)第1排序軸和所有排序軸均達(dá)到極顯著水平(P<0.01)，表明結(jié)果具有統(tǒng)計(jì)學(xué)意義。由圖4可知，海拔、年降水量和年均溫度與第1排序軸呈正相關(guān)關(guān)系，且相關(guān)性依次為海拔>年降水量>年均溫度；海拔與紅砂比葉面積、P含量呈負(fù)相關(guān)關(guān)系。

圖4 紅砂葉功能性狀與種源地環(huán)境因子的RDA排序Fig.4 Redundancy analysis ordination of leaf functional traits of R. soongarica and geographic environmental factors

表4 白刺葉性狀與環(huán)境因子的RDA排序Table 4 Redundancy analysis ordination of leaf traits and environmental factors of N.tangutorum

表5 紅砂葉性狀與環(huán)境因子的RDA排序Table 5 Redundancy analysis ordination of leaf traits and environmental factors of R. soongarica

圖3 白刺葉功能性狀與種源地環(huán)境因子的RDA排序Fig.3 Redundancy analysis ordination of leaf functional traits of N. tangutorum and geographic environmental factors注：經(jīng)度-Lng,緯度-Lat,海拔-ALT,年均溫度-MAT,年均降水量-MAP,下同Note:Longitude-Lng,Latitude-Lat,Altitude-ALT,Annual average temperature-MAT,Annual average precipitation-MAP,the same as below

2.4.2葉功能性狀與種源地環(huán)境因子的逐步線性回歸分析 由表6可知，白刺葉片肉質(zhì)化程度、含水量與海拔呈極顯著正相關(guān)關(guān)系(P<0.01)；葉組織密度與海拔、緯度均呈極顯著負(fù)相關(guān)關(guān)系(P<0.01)，且海拔的標(biāo)準(zhǔn)化回歸系數(shù)的絕對(duì)值(|—0.697|)大于緯度的標(biāo)準(zhǔn)化回歸系數(shù)的絕對(duì)值(|—0.441|)，表明影響葉組織密度的主要環(huán)境因子是海拔，其次是緯度；葉片N含量與年均降水量呈顯著負(fù)相關(guān)關(guān)系(P<0.05)，與1月均溫呈顯著正相關(guān)關(guān)系(P<0.05)，且年均降水量的標(biāo)準(zhǔn)化回歸系數(shù)的絕對(duì)值(|—2.704|)大于1月均溫的標(biāo)準(zhǔn)化回歸系數(shù)(2.126)，表明影響葉片N含量的主要環(huán)境因子是年均降水量，其次是1月均溫；葉片P含量與緯度呈顯著負(fù)相關(guān)關(guān)系(P<0.05)，葉片N∶P與緯度呈顯著正相關(guān)關(guān)系(P<0.05)，與7月均溫呈顯著負(fù)相關(guān)關(guān)系(P<0.05)，且緯度的標(biāo)準(zhǔn)化回歸系數(shù)(0.708)大于7月均溫的標(biāo)準(zhǔn)化回歸系數(shù)的絕對(duì)值(|-0.441|)，這表明影響葉片N∶P的主要環(huán)境因子是緯度，其次是7月均溫。因此，白刺葉片肉質(zhì)化程度、含水量、組織密度主要受海拔的影響，葉片N含量主要受年均降水量的影響，葉片P含量、N∶P主要受緯度的影響。而紅砂葉片的比葉面積和P含量主要受海拔影響。

表6 白刺和紅砂葉功能性狀與環(huán)境因子的逐步回歸分析Table 6 Stepwise regression analysis of leaf functional traits and environmental factors of N. tangutorum and R. soongarica

3 討論

3.1 不同地理種源白刺和紅砂葉片形態(tài)性狀的變異特征

植物在長(zhǎng)期適應(yīng)種源環(huán)境下，必然會(huì)產(chǎn)生生理生態(tài)特性的分化，形成不同的地理種源，導(dǎo)致同一葉功能性狀在不同種源間差異顯著[32]。大量研究表明，植物長(zhǎng)期受種源地環(huán)境影響，導(dǎo)致其葉功能性狀產(chǎn)生遺傳分化[33-35]。本研究結(jié)果也顯示，不同地理種源白刺在同質(zhì)園環(huán)境中葉片肉質(zhì)化程度、含水量、組織密度表現(xiàn)出顯著性差異(圖1)，表明白刺受種源地環(huán)境的影響，其葉功能性狀差異可能與遺傳分化有關(guān)。有研究表明，在野外環(huán)境和同質(zhì)園環(huán)境中，植物由于受種源環(huán)境的影響，產(chǎn)生了遺傳穩(wěn)定性[17]。楊彪生等[36]研究發(fā)現(xiàn)，野外環(huán)境下白刺葉片肉化程度、含水量和組織密度均呈現(xiàn)顯著性差異，這一野外研究結(jié)果與本研究中同質(zhì)園實(shí)驗(yàn)的研究結(jié)果一致，說(shuō)明白刺葉功能性狀通過(guò)對(duì)種源環(huán)境的適應(yīng)可能形成了遺傳穩(wěn)定性。楊彪生等[36]對(duì)野外環(huán)境中紅砂葉片肉質(zhì)化程度、含水量和組織密度的研究發(fā)現(xiàn)這些指標(biāo)均呈顯著性差異。然而，本研究中同質(zhì)園內(nèi)不同地理種源紅砂葉片的肉質(zhì)化程度、含水量、組織密度均無(wú)顯著性差異(圖1)，說(shuō)明紅砂在同質(zhì)園環(huán)境中不受種源環(huán)境的影響。

比葉面積可以反映植物對(duì)資源的獲取能力，且也能夠反映出植物對(duì)不同環(huán)境的適應(yīng)性特征，是研究植物功能性狀的重要指標(biāo)之一[37-38]。有研究認(rèn)為，植物在資源豐富的環(huán)境中生長(zhǎng)比葉面積較大，在資源匱乏的環(huán)境中生長(zhǎng)往往比葉面積較小[39]。在本研究中，同質(zhì)園中紅砂的比葉面積(311.26 cm2·g-1)高于野外環(huán)境中紅砂的比葉面積(277.32 cm2·g-1)，這可能是同質(zhì)園環(huán)境有較充足的水分、光照等條件，與野外環(huán)境相比，同質(zhì)園中比葉面積利用光合的能力較強(qiáng)，能夠快速獲取所需資源。

3.2 不同地理種源白刺和紅砂葉片C，N，P含量化學(xué)計(jì)量及其計(jì)量比的變異特征

葉片C，N，P含量對(duì)異質(zhì)環(huán)境的適應(yīng)是表型分化和遺傳分化共同的結(jié)果，并在長(zhǎng)期的進(jìn)化過(guò)程中形成不同的生理生態(tài)和養(yǎng)分利用特征，若在同質(zhì)園環(huán)境中，這種生理生態(tài)的分化依然存在，則認(rèn)為植物為了適應(yīng)環(huán)境產(chǎn)生了遺傳分化[17]。有研究表明，同質(zhì)園試驗(yàn)植物葉片的C，N，P含量及其計(jì)量比在種源間表現(xiàn)出顯著性差異，說(shuō)明植物葉片中化學(xué)計(jì)量特征在對(duì)種源地環(huán)境條件的長(zhǎng)期進(jìn)化適應(yīng)過(guò)程中產(chǎn)生了遺傳分化[40-41]；本研究發(fā)現(xiàn)白刺葉片N含量，N∶P在不同地理種源間表現(xiàn)出的結(jié)果與上述結(jié)果相同(圖2)，說(shuō)明白刺受種源環(huán)境的影響可能使其葉片化學(xué)元素產(chǎn)生了遺傳分化。不同種源翅莢木(Zeniainsignis)葉片C含量、N含量及C∶N存在差異顯著，說(shuō)明其主要葉片養(yǎng)分性狀可能主要受長(zhǎng)期的種源地環(huán)境影響產(chǎn)生了遺傳分化[32]，這與本文研究中紅砂葉片C含量、N含量及C∶N無(wú)顯著差異的結(jié)論相悖(圖2)，這可能是因?yàn)槌崆v木研究地屬于中亞熱帶海洋性季風(fēng)氣候，雨水充沛，氣候溫暖，而本文研究地屬于干旱性大陸氣候，干旱少雨，兩個(gè)研究地的立地條件和氣候環(huán)境截然不同，所以對(duì)其葉功能性狀的影響程度也不同。

C，N，P含量化學(xué)計(jì)量比對(duì)植物營(yíng)養(yǎng)元素的平衡起限制性作用，并已經(jīng)廣泛運(yùn)用于植物養(yǎng)分循環(huán)和限制的研究[42-43]。N∶P常被看作植物在生長(zhǎng)過(guò)程中獲取環(huán)境養(yǎng)分的限制性因素[1，44]。有研究表明，植物生長(zhǎng)發(fā)育受N和P限制的N∶P比值為14和16，即當(dāng)N∶P<14時(shí)，植物生長(zhǎng)發(fā)育主要受N限制；當(dāng)N∶P>16時(shí)，植物生長(zhǎng)發(fā)育主要受P限制；當(dāng)14

3.3 白刺和紅砂葉功能性狀間的關(guān)系

植物葉片在長(zhǎng)期適應(yīng)地理環(huán)境的過(guò)程中，通過(guò)內(nèi)部不同功能間的協(xié)調(diào)，形成一系列適應(yīng)環(huán)境的葉功能性狀間的不同組合，并非葉片性狀孤立發(fā)揮作用[47-48]。通常認(rèn)為，植物葉片的比葉面積與組織密度呈負(fù)相關(guān)的關(guān)系，即比葉面積增大時(shí)，組織密度會(huì)變小。比葉面積和組織密度作為植物資源分配利用的重要指標(biāo)，均能夠反映植物對(duì)環(huán)境的適應(yīng)策略(水分利用策略)[49-51]。本研究發(fā)現(xiàn)，白刺和紅砂的比葉面積與組織密度呈極顯著負(fù)相關(guān)關(guān)系(表2、表3)，出現(xiàn)這種現(xiàn)象可能與植物地處干旱少雨的環(huán)境密切有關(guān)，即當(dāng)植物處于干旱環(huán)境時(shí)，植物為了適應(yīng)干旱環(huán)境，從周?chē)h(huán)境中不斷汲取養(yǎng)分，使得比葉面積增大，但此過(guò)程中植物對(duì)光照需求較大，而對(duì)葉片中儲(chǔ)存的養(yǎng)分和水分的需求較小，故養(yǎng)分資源的利用效率相對(duì)較低，因此植物葉片相對(duì)較薄，葉組織密度相對(duì)較小。因此，白刺和紅砂葉片通過(guò)比葉面積和組織密度的平衡形成對(duì)不同水分環(huán)境的適應(yīng)策略[47，52]。有研究表明，作為植物生長(zhǎng)中主要的營(yíng)養(yǎng)元素，N和P具有一定的相互促進(jìn)作用，N∶P可以反映N和P兩種營(yíng)養(yǎng)元素供給狀況的相對(duì)有效性[53]。白刺和紅砂葉片N含量與N∶P呈顯著正相關(guān)的關(guān)系(表2、表3)，這是由于植物對(duì)養(yǎng)分的利用受到了環(huán)境和自身需求的影響，植物通過(guò)調(diào)節(jié)自身葉片營(yíng)養(yǎng)元素及化學(xué)計(jì)量比來(lái)適應(yīng)環(huán)境中養(yǎng)分的供給情況[41]。葉片C∶N和C∶P表示植物吸收N，P元素所能同化碳的能力，是植物養(yǎng)分利用效率的重要指標(biāo)[54]。本研究中，紅砂葉片C∶N與C∶P呈極顯著正相關(guān)關(guān)系(表3)，說(shuō)明紅砂葉片具有較高的碳同化能力。

3.4 種源地環(huán)境對(duì)白刺和紅砂葉功能性狀的影響

廣泛分布的植物對(duì)其環(huán)境適應(yīng)的范圍較廣，可以通過(guò)表型可塑性、局域適應(yīng)等來(lái)面對(duì)環(huán)境變化，其選擇的策略將決定植物應(yīng)對(duì)氣候變化的能力[55]，一般來(lái)說(shuō)，葉性狀若與種源地環(huán)境因子顯著相關(guān)，則可認(rèn)為發(fā)生局域適應(yīng)[56]。結(jié)合RDA分析和逐步線性回歸法分析，白刺葉片N含量與年均降水量有相關(guān)性(圖3，表6)，這可能是因?yàn)楸狙芯?個(gè)種源地干旱少雨，除LZ種源外，其余種源年均降水量均不足200 mm，而且水分是植物生長(zhǎng)發(fā)育過(guò)程中所必需的因子，與其他因子相比，植物功能性狀對(duì)水分的響應(yīng)更為顯著[33]，因此說(shuō)明白刺葉片N含量主要受水分的影響。一般來(lái)說(shuō)，海拔對(duì)植物功能性狀的影響機(jī)理主要是隨著海拔的升高，植物生長(zhǎng)發(fā)育逐漸受到溫度和水分的限制，因而植物對(duì)環(huán)境的適應(yīng)性會(huì)發(fā)生顯著變化[57]。本研究結(jié)果顯示，白刺肉質(zhì)化程度和含水量與海拔呈極顯著正相關(guān)關(guān)系，組織密度與海拔呈極顯著負(fù)相關(guān)關(guān)系(圖3，表6)，并且海拔對(duì)植物比葉面積、組織密度、干物質(zhì)含量等影響顯著的結(jié)論已得到證實(shí)[57-58]，因此白刺肉質(zhì)化程度、含水量和組織密度隨海拔的變化受溫度和水分的限制。海拔對(duì)紅砂葉片比葉面積和P含量的影響呈顯著負(fù)相關(guān)關(guān)系(圖4，表6)，并且隨著海拔的上升，比葉面積和P含量趨向于降低(圖1，圖2)，這可能是由于本研究所采樣的4個(gè)種源地均屬于高海拔地區(qū)，使得植物生長(zhǎng)所需要的溫度和水分受到了影響，植物為了適應(yīng)高海拔的惡劣環(huán)境，其物質(zhì)代謝發(fā)生改變，導(dǎo)致植物獲取資源的能力和對(duì)環(huán)境的適應(yīng)能力發(fā)生了變化[59]。因不同種源地植物葉功能性狀的變異受到多個(gè)地理環(huán)境因子的共同影響，本文只選取了緯度、經(jīng)度、海拔、年均降水量、年均溫度等環(huán)境因子對(duì)葉片功能性狀進(jìn)行研究，沒(méi)有涉及到各種源地物種所在地的土壤養(yǎng)分、坡向、相對(duì)濕度、積溫等環(huán)境因子。因此，在今后的研究中，還需在這方面作進(jìn)一步研究。

4 結(jié)論

通過(guò)研究同質(zhì)園中各地理種源白刺和紅砂葉功能性狀的種間差異，發(fā)現(xiàn)白刺在適應(yīng)環(huán)境的過(guò)程中可能產(chǎn)生了遺傳分化，且受環(huán)境影響可能形成了遺傳穩(wěn)定性，而紅砂在同質(zhì)園中不受種源環(huán)境的影響；兩種植物葉片生長(zhǎng)發(fā)育主要受氮素限制；受環(huán)境的影響，兩種植物葉片內(nèi)部功能性狀的平衡變化形成了對(duì)不同環(huán)境的適應(yīng)策略；分析白刺和紅砂葉功能性狀與種源環(huán)境因子的關(guān)系，發(fā)現(xiàn)海拔是影響兩種植物的主導(dǎo)環(huán)境因子，故兩種植物生長(zhǎng)受溫度和水分的限制。