崇明東灘濕地蘆葦種群內功能性狀及生態對策變異分析

計陽晴, 李德志,1b,2,3*, 王虓, 李羚菱, 王瑩, 陳婧, Rasool Samreen Ghulam, Tehreem Rafique

崇明東灘濕地蘆葦種群內功能性狀及生態對策變異分析

計陽晴1a, 李德志1a,1b,2,3*, 王虓1a, 李羚菱1a, 王瑩1a, 陳婧1a, Rasool Samreen Ghulam1a, Tehreem Rafique1a

(1. 華東師范大學, a. 生態與環境科學學院; b. 上海市城市化生態過程與生態恢復重點實驗室, 上海 200241; 2. 崇明生態研究院, 上海 202162; 3. 自然資源部大都市區國土空間生態修復工程技術創新中心, 上海 200062)

植物種群功能性狀及生態對策變異對其自身環境適應能力及群落結構具有重要影響。該文基于崇明東灘濕地蘆葦()種群內3個小種群不同發育時期分株集合的多項性狀指標，分析了其性狀及生態對策變異，并基于CSR、LHS生態對策體系，對各小種群的生態對策進行了排序。結果表明, 3個蘆葦小種群的21項性狀指標及C、R對策得分均存在顯著差異。在LHS體系中，以比葉面積(11.5 mm2/mg)、株高(175 cm)和花序干質量(4 g)為分界線，可對3個小種群進行區分。處于生殖生長期的蘆葦分株具有更大的植株大小、更高的葉片干物質含量及S對策得分，而處于營養生長期的分株則具有較小的植株大小、較大的比葉面積及R對策得分。總體而言，崇明東灘濕地3個蘆葦小種群間及各小種群內處于不同生長發育階段的分株集合間，都存在顯著的性狀差異及生態對策差異，且二者的差異模式基本吻合。蘆葦種內性狀及生態對策變異特征，可為該區蘆葦種群的生態保育和恢復實踐提供重要的理論參考和技術支撐。

種內性狀差異；生態對策；CSR體系；LHS體系；蘆葦；濕地

種群內功能性狀變異指的是同一種群內不同的小種群間或個體間的功能性狀差異，屬于種內變異。Daniel等[1]總結了種內性狀變異改變群落結構或動態的6種機制，認為忽略種內性狀變異，可能會影響群落研究結果的正確性。此后的一項研究量化了種內性狀變異對群落總體性狀變異的貢獻程度，證實了種內性狀變異，對群落植物功能類型研究結果的潛在影響[2]。

植物生態對策是植物在適應環境脅迫和干擾、維持種群存續的過程中所產生的一系列特征[3]。隨著全球尺度植物生態學研究的不斷推進，科學家們提出了多種生態對策體系，通過刻畫植物對多維度環境因子的適應性權衡來劃分植物功能類型，從而達到簡化群落研究的目的[4]。然而，無論生態對策分類體系如何發展，植物性狀始終都是表征植物生態對策的基礎指標，因此，同一物種的不同個體, 在性狀上若存在顯著差異，即意味著不同個體在生態對策上也可能存在較為明顯的區別。一些研究表明，同一物種的不同種群，在性狀上的顯著差異, 的確使得不同種群在生態對策體系中的分布有所區別[5–6]。但也有的研究認為，現有的生態對策體系，無法反映出種內性狀變異的真實模式及變異水平[7–8]。這些研究結果提醒我們，在利用生態對策體系對群落結構進行簡化時，若簡單地將某一物種當成一個整體來討論，可能會對該物種在群落中的功能產生錯誤評估，但一個具有高度種內性狀變異性的物種，是否也會具有種內生態對策多樣性還有待驗證。

CSR體系(competitors-stress tolerators-ruderal species)和LHS體系(leaf-height-seed)是2種在植物生態學研究中應用最為廣泛的生態對策體系。前者認為，所有的植物在生存環境中，都將面臨競爭、脅迫及干擾3種環境因子的挑戰。植物生態對策應當是植物對三者適應性的綜合體現，并以3種對策間的權衡關系進行表征[9]。該體系自Grime首次提出以來，經歷了3次分類方法的革新[10]，被逐步應用、擴展，并分別向上和向下延伸至全球水平[11]和個體水平[12]。目前已有研究案例，應用Pierce等[11]所開發的電子表格工具“StrateFy”，成功揭示出種內不同種群間，以及同一物種不同個體發育階段間的生態對策差異[5–13]。LHS體系與“StrateFy”的建立基礎不同，除葉片性狀之外，它還包括植株整體大小(株高)及種子(質量)的相關性狀。直接利用比葉面積、株高及種子質量3項指標的性狀值進行三維投射，即可比較目標對象在LHS體系中的分布差異[3]。

蘆葦()是一種在世界范圍內廣泛分布的禾本科(Poaceae)植物, 由多種單倍型組成(haplotype)[14]。我國境內的蘆葦種群具有多樣化的單倍型[15]，是常用于濕地治理和造紙的、具有生態和經濟雙重價值的本土物種[16]。自19世紀蘆葦被引入北美洲之后[17]，單倍型M便迅速入侵了當地生態系統。在加拿大某些地區，該單倍型甚至取代了當地蘆葦單倍型[18]。于北美洲而言，蘆葦種下水平的多樣性使其具有本土種和入侵種的雙重身份，深刻地影響著當地生態系統的結構與功能。蘆葦在基因-生態系統水平上的多樣性，使得其成為研究植物環境適應對策的模式生物[19]。作為崇明島濕地的優勢種之一，蘆葦不僅對整個濕地的結構與功能具有重要影響，還能為遷徙鳥類提供越冬及繁殖的生境和食物[20]。據觀察，分布在崇明東灘濕地的蘆葦種群，沿東西方向呈現出較為明顯的形態學差異。此外，在生長季末期時，蘆葦種群中同時存在處于繁殖期的個體(分株)和少部分仍處在營養生長期的個體(分株)，且二者的形態學特征也存在較為明顯的區別。崇明東灘濕地蘆葦種群結構所展現出的多樣性，使其成為研究種內性狀及生態對策變異的合適樣地。

目前國外利用生態對策理論體系探討植物種群生態對策的案例相對較多[5–8]，而國內在此方面的研究則較為稀缺[10]，只有少數針對種間水平的研究[21–22]。有研究曾對黃河三角洲的蘆葦種內CSR對策差異進行過分析[23]，但并沒有關注蘆葦種群內不同發育階段分株間的生態對策差異，也沒有比較CSR體系和LHS體系在蘆葦種群內功能性狀及生態對策變異研究方面的差異和實際效果。

為探究崇明東灘濕地蘆葦種群的種群內功能性狀及生態對策變異，本文依據形態學性狀和分布特征，將該蘆葦種群劃分為3個小種群(或稱亞種群),將處于繁殖生長期的分株與處于營養生長期的分株視為2個不同的分株集合，比較了3個小種群與2個分株集合的25項和20項功能性狀，及3個小種群在2種生態對策體系中的特征，并提出以下假設：(1) 崇明東灘蘆葦種群存在較高的種群內功能性狀變異，表現為不同小種群及不同發育階段的分株集合在所調查的功能性狀上存在顯著差異；(2) 蘆葦種群具有生態對策變異，且不同小種群在CSR和LHS體系中的表現存在顯著差異，不同發育階段分株集合的CSR得分也存在顯著差異。本研究旨在為深刻認識崇明東灘濕地沿環境梯度分布的蘆葦種群的適應機制提供嶄新視角；揭示的蘆葦種群內功能性狀及生態對策變異特征，有望為該區蘆葦種群的生態保育和恢復實踐，提供重要的理論參考和技術支撐。

1 材料和方法

1.1 研究區域及采樣設置

研究地點位于上海市崇明東灘濕地鳥類保護區的東部。沿東西方向，把廣泛分布于東灘濕地的蘆葦種群，劃分為3個小種群，從東至西分別為: 前沿近海小種群(31°30′37″ N, 121°59′15″ E)、中部小種群(31°30′10″ N, 121°58′5″ E)及后沿近陸小種群(31°30′24″ N, 121°56′55″ E)。其中，前沿近海小種群呈斑塊狀分布，中部小種群及靠近98大堤的后沿近陸小種群沿南北方向呈帶狀分布。

于2021年10—11月對3個小種群進行采樣。在中部小種群、后沿近陸小種群內，沿98大堤(南北)方向隨機選取3個50 cm×50 cm的小樣方，且小樣方間相隔20 m以上。在前沿近海小種群中，隨機選取3個50 cm×50 cm的小樣方，以代表整個斑塊小種群的生長狀況。計數每個小樣方內，蘆葦的植株密度及結實率(即結實株數占植株總數的比例)，然后，在每個小樣方內，隨機選取5株處于繁殖生長階段的植株和5株處于營養生長階段的植株,測量其株高和基莖。每株植株擇莖干中上部的一片完整葉片剪下，放入黑色塑料封口袋中，供后續指標的測量；其余葉片及花序剪下后，放至信封中, 并將植株齊地收獲帶回實驗室。

1.2 植株性狀及測量方法

針對崇明東灘鹽沼濕地的生境特征，本研究選擇了25項功能性狀指標對蘆葦種群潛在的種群內功能性狀差異進行分析。將采集的葉片放入裝滿水的、黑色不透光的密封袋內，于4 ℃冷藏1晚后用紙巾擦凈表面水分，進行飽和鮮重(LSW)的測定; 隨后用佳能掃描儀LIDE 110對葉片進行掃描，再利用Image J (2.1.0.)測量葉長(Lleaf)、葉寬(Wleaf)與葉面積(LA)；最后將上述葉片放入烘箱，80 ℃烘干后稱量葉片干重(LDW)，并計算葉片干物質含量(LDMC)及比葉面積(SLA)[12]。干重測量完畢后將葉片粉碎用于后續生理生態學指標的測量。將信封中的葉片一并放入80 ℃烘箱中烘干，葉片總干重(ALDW)等于信封中的葉片干重與密封袋中單葉的干重之和。依據葉痕判斷葉片總數(LN)及葉片死亡數(DL); 未離體的葉片則以50%的葉面積枯萎為標準進行死亡判定，葉片死亡率(DLR)為葉片死亡數占葉片總數的百分比。用卷尺測量花序長度(Linf)，再將花序及采集的莖干部分放入80 ℃烘箱中烘干進行花序干重(DWinf)與莖干重(SDW)的測量。地上生物量為各器官干重的總和。對于處于繁殖時期的分株，利用各器官干重占地上生物量的百分比計算葉片(rleaf)、莖(rstem)及花序生物量(rinf)分配率。

采用蒽酮比色法測定非結構性碳水化合物含量(NSC)。取0.05 g干葉樣品，加入7.5 mL蒸餾水，在沸水浴中提取20 min。冷卻過濾后收集濾液，用蒸餾水定容至50 mL。取1 mL提取產物與4 mL 0.1% ()的蒽酮溶液混合均勻(蒽酮溶液由體積分數為80%的H2SO4溶液制備)。在沸水浴中浸泡10 min后，于620 nm波長下測定吸光值。同時制備0、10、20、40、60、80、100g/L的葡萄糖溶液以繪制標準曲線。NSC=(×)/(×)×100%, 式中，為在標準曲線上查出的糖含量(g)，為提取液總體積(50 mL)，為樣品重量(50g)，為測定時取用的體積(1 mL)。

采用福林酚試劑法[24]測定總酚含量, 取0.2 g粉碎的干葉片樣品，放入40 mL的蒸餾水中，于80 ℃的水浴鍋內，提取30 min，迅速冷卻過濾，濾液定容至50 mL后備用。吸取1 mL待測樣品，加入5 mL蒸餾水、1 mL福林酚試劑及3 mL質量濃度為7.5%的Na2CO3溶液，混合均勻后，于暗處顯色2 h，在765 nm波長下，測定反應體系的吸光值。配置濃度為0、10、20、30、40、50 mg/L的沒食子酸溶液，以繪制標準曲線。沒食子酸標準溶液吸光值的測定步驟與樣品相同?？偡雍繀⒄誎apali等[25]的公式計算。

參照Tong等[26]的方法測量葉綠素含量，將0.05 g干葉片樣品放入體積分數為95%的乙醇溶液中，浸提24 h，計算葉綠素a (Chl a)、b (Chl b)及總葉綠素含量(Chl a+b)。

1.3 蘆葦分株在CSR及LHS生態對策體系中的分布

利用Pierce等[11]提出的Excel工具“StrateFy”對蘆葦分株的CSR對策得分進行計算并分類；通過輸入葉面積、葉片飽和鮮重及葉片干重3項原始數據，即可返回CSR對策得分結果及分株對策類型。

本研究選擇花序干重代替種子質量，進行LHS體系的構建。此外，由于種群內水平功能性狀變異的范圍要遠小于種間水平，因此，本研究將采用比葉面積、株高和花序干重數據的絕對值來比較3個小種群在體系中的分布差異。

1.4 數據處理與分析

生長末期前沿近海、中部及后沿近陸小種群的結實率分別為87.94%、73.56%及77.59%，這意味著3個小種群內的大部分植株都處于繁殖階段，因此，應選擇繁殖分株的性狀數據集，來代表3個小種群的性狀變異情況，以進行后續分析。

對3個小種群的25項功能性狀進行單因素方差分析。若數據既符合正態分布，也符合方差齊次性，則用R 4.2.0中的“aov”函數，進行顯著性檢驗，同時，使用“agricolae”包進行事后多重比較。當數據不符合正態分布或方差齊次性時，使用“PMCM- Rplus”包，進行方差分析及多重比較。對于葉片總數、葉片枯萎數等不連續型性狀，可以將其視為有序分類變量，分析方法與不符合正態分布或方差齊次性的數據相同。對“StrateFy”返回的C、S、R對策得分，以同樣的方法，進行單因素方差分析，同時，利用R中的“ggtern”函數，繪制各樣本在CSR體系中的分布圖。

用R中自帶的“t.test”函數，對處于繁殖生長期的分株集合與處于營養生長期的分株集合的性狀差異，進行獨立樣本檢驗。對于不滿足正態分布或方差齊次性的數據，以及離散型的性狀指標，如葉片數、葉片枯萎數，則需使用“wilcox.test”函數,進行非秩序性檢驗。不同個體發育時期分株集合, 在CSR對策得分上的差異，以同樣的方式進行分析; CSR分布圖采用“ggtern”函數繪制。

為了更好地展示蘆葦種群處于繁殖生長期的個體在LHS體系中的分布情況，本研究選擇對比葉面積、株高和花序干重性狀值進行二維投影，利用“cor.test”函數，對3個性狀間的相關關系進行檢驗，并用“ggplot2”包繪制二維體系分布圖。

2 結果和分析

2.1 功能性狀的差異分析

生長末期3個小種群(以處于繁殖生長時期的分株代表整個小種群)在25項功能性狀指標中，僅有4項未見顯著差異，分別是Wleaf、LDMC、DLR和NSC (表1)。后沿近陸小種群繁殖分株在葉片、莖、花序及植株整體大小上(形態學大小包括長度和重量兩個方面)都占據優勢地位，而中部小種群繁殖分株雖然在葉片、莖、植株整體大小以及花序長度上都顯著小于后沿近陸小種群的繁殖分株，但其花序重量與后沿近陸小種群間，不存在顯著差異。這與3個小種群在地上生物量分配方面的特征相符。分析結果表明，前沿近海小種群在葉片生物量分配上，顯著高于其他2個小種群，而后沿近陸小種群則傾向于將大部分生物量分配給莖部，以支撐植株高度。與前沿近海小種群不同的是，盡管中部小種群的地上總生物量較低，但其仍傾向于將較多生物量分配于花序生長。生理生態學指標數據表明，前沿近海小種群具有較高的總酚含量，而中部小種群的葉綠素含量顯著高于其他2個小種群。在與CSR對策體系直接相關的3項葉片性狀指標中，LDMC在各小種群間不存在顯著差異，而前沿近海小種群和后沿近陸小種群分別在SLA和LA值上占據優勢。

表1 蘆葦小種群的功能性狀特征

同行數據后不同小寫字母表示差異顯著(<0.05)。

Data followed different letters in the same line indicate significant differences at 0.05 level.

前沿近海、中部和后沿近陸3個小種群分別有18、12和15項功能性狀指標在不同發育階段的分株集合間存在顯著差異(表2)。3個小種群內不同發育階段分株在Wleaf和NSC2項指標上均無顯著差異。形態學指標分析結果顯示，各小種群的生殖生長分株，在莖干及植株大小性狀值、LDMC、DL和DLR上均顯著高于營養生長分株。生理生態學指標的分析結果則表明，營養生長分株的葉片總酚含量及葉綠素含量顯著高于生殖生長分株(僅后沿近陸小種群的營養生長分株與生殖生長分株在葉片總酚含量上不存在顯著差異)。前沿近海與后沿近陸小種群生殖生長分株的LA值顯著大于營養生長個體；各小種群生殖生長分株的SLA值均顯著小于營養生長分株，而各小種群生殖生長分株的LDMC值則顯著高于營養生長分株。

2.2 CSR體系的排序

生長末期各小種群在CSR的3項對策得分(即各小種群生殖生長分株的CSR得分)如下：前沿近海小種群(C=34.93±3.42, S=58.44±2.39, R=6.63±2.38),中部小種群(C=35.34±5.40, S=60.38±5.96, R=4.29± 2.30)，后沿近陸小種群(C=39.11±4.62, S=58.18±4.86, R=2.72±2.33)。其中后沿近陸小種群的C對策得分顯著高于其他2個小種群(<0.05)，而前沿近海小種群的R對策得分則在3個小種群中占據極顯著優勢(<0.001)，各小種群的S對策得分不存在顯著差異。生長末期種群內的分株主要屬于CS、S/CS和S/CSR等3種類型。其中前沿近海和中部小種群的分株主要為前兩類對策，所占比例分別為66.67%、33.33%和46.67%、53.33%。后沿近陸小種群3種對策類型的分株比例分別為60%、33.33%和6.67%。綜合來看，該濕地蘆葦種群整體偏向CS對策類型，其中中部小種群比其他2個小種群更偏向S對策(圖1: D)。

不同發育階段分株集合在CSR體系中的排序如圖1: A～C所示。前沿近海、中部及后沿近陸小種群營養生長分株的CSR得分分別為(C=32.19±2.52, S= 47.80±3.49, R=20.02±3.01)、(C=37.05±7.59, S=47.13± 11.75, R=15.82±11.30)和(C=35.82±4.26, S=53.04± 5.39, R=11.15±4.07)。前沿近海小種群10—11月營養生長分株共包含3種對策類型，分別為S/CSR(46.67%)、CS/CSR (46.67%)及CSR (6.67%)，中部小種群營養分株共包含S/CSR、CS/CSR、CSR、CS、CR/CSR和S/CS共5種對策類型(分別占26.67%、13.33%、6.67%、40%、6.67%和6.67%)，而后沿近陸小種群的營養分株則有S/CSR (40%)、CS/CSR(26.67%)和CS (33.33%)共3種對策類型。綜合來看，各小種群營養生長分株的R對策得分均顯著高于生殖生長分株(<0.001)，而其S對策得分則顯著低于后者(前沿近海和中部小種群的顯著水平為=0.001,后沿近陸小種群的顯著水平為=0.05)。在C對策得分上，僅有前沿近海小種群的生殖生長分株與營養生長分株間存在顯著差異(<0.05)。

表2 蘆葦不同發育階段分株集合的功能性狀

*:<0.05; **:<0.01; ***:<0.001

圖1 蘆葦各小種群和不同發育階段分株集合的CSR排序。A: 前沿近海; B: 中部; C: 后沿近陸; D: 生長末期。

2.3 各小種群在LHS體系中的排序

生長末期3個小種群在LHS體系中的二維排序如圖2所示。在崇明東灘濕地蘆葦種群中，LHS體系的3個性狀指標間存在顯著的相關性，具體表現為SLA與其他2個性狀呈負相關，而株高和花序干重則呈顯著正相關關系。各小種群在3項指標的搭配上存在差異，使得LHS體系可以把不同小種群區分開來(大致以SLA=11.5 mm2/mg、株高=175 cm和DWinf=4 g為分界線)。前沿近海小種群個體傾向于高SLA值、較低的株高和較小的DWinf；中部小種群個體傾向于低SLA值、較低的株高和較大的DWinf；后沿近陸小種群個體則傾向于較低的SLA值、較高的株高和較大的DWinf。較大的SLA值與較小的DWinf是前沿近海小種群有別于其他2個小種群的特征，而中部小種群與后沿近陸小種群則可以依靠株高進行區分。

3 結論和討論

3.1 崇明東灘濕地蘆葦種群內功能性狀及生態對策變異

通過性狀對比分析、CSR體系排序及LHS體系排序研究，崇明東灘濕地3個蘆葦小種群間及各小種群內不同發育階段的分株集合間都存在較為顯著的功能性狀差異及生態對策差異。

沿東灘濕地東西方向，蘆葦種群呈現出C對策得分逐漸上升，而R對策得分逐漸降低的趨勢。各小種群間CSR對策得分的差異模式基本符合蘆葦種群水平的功能性狀變異特征，且CSR對策的差異, 主要由植株大小性狀的差異導致。3個小種群在以SLA、株高及DWinf構建的LHS體系中，具有明顯的排序差異。

在崇明東灘濕地，生境脅迫因子為塑造蘆葦功能性狀和生態對策的主導因素，而干擾因子則是蘆葦種群內功能性狀和生態對策變異產生的原因。水鹽因子是鹽沼濕地的主要限制性因子[27]，二者會對植物生長產生脅迫作用。有研究表明，NSC是植物抗逆性的重要指標之一，環境脅迫會促使植物體內的可溶性糖積累[28–30]。而本研究中，蘆葦各小種群的NSC含量未見顯著差異，這說明盡管在該生境內存在鹽分梯度[31]，但鹽脅迫對各小種群的影響程度大致相同。此外，很多植物還能夠通過增加葉片周轉率來抵抗脅迫因子[32]。崇明東灘蘆葦種群盡管在葉片死亡數和葉片總數上都存在顯著的種內差異, 但是各小種群的葉片死亡率卻保持基本一致，這一結果同樣能夠驗證上述觀點。CSR對策分析結果表明，蘆葦各小種群的S對策得分無顯著差異，且蘆葦種群主要呈現CS對策，蘆葦種群內功能性狀及生態對策的差異模式說明鹽沼濕地以脅迫為主導環境因素，但不同脅迫程度所產生的影響無顯著差異。同時，在該生境內，離海更近的區域受潮汐以及波浪干擾作用明顯[31]，因此，前沿近海小種群呈現出較高的R對策及總酚含量(總酚含量可被用來指征干擾因子的影響程度[33–34])。而后沿近陸小種群由于受到的干擾較少，生物量積累時間較多，因此株高一般較大、C對策得分較高。

有研究表明，植物在個體發育過程中，將經歷由R對策向S或C對策的轉變，且植物幼苗的SLA及LA值往往大于成熟的植株個體，而LDMC則小于后者[13]。蘆葦顯然也符合這一規律，蘆葦幼苗(處于營養生長的分株)在植株大小上顯著小于生殖生長期的成熟分株，同時，較高的SLA值說明幼苗更傾向于將資源和能量投入至與資源獲取相關的性狀中[12]，這樣的特征使蘆葦幼苗比處于生殖生長期的分株，具有更高的R對策得分及較低的S對策得分。此外，在本研究中，不同發育階段分株集合處在同一時間節點、同一小種群內，因此，環境因子可能不是二者產生差異的主要原因。

3.2 基于CSR和LHS體系的蘆葦種群內功能性狀及生態對策比較

有研究認為，種內性狀的高度變異性不會導致CSR對策差異[8]，原因可能是因為其僅僅對CSR對策類型進行了探討，而沒有對3項對策的具體得分進行比較分析。本研究結果表明，“StrateFy”表格工具確實能夠揭示同一生境內蘆葦的種群內功能性狀變異特征。

盡管應用CSR體系確實檢驗了東灘濕地蘆葦種群內功能性狀變異的存在，但該分類方法在種內水平應用時，仍然表現出一些不足之處，具體體現在CSR體系無法反映蘆葦花序性狀的變異特征。一些研究發現，繁殖性狀(如種子大小)的變異特征可能獨立于葉片經濟譜系及植株大小的變異特征，并成為植物應對干擾因子的主要對策[35–36]。在本研究中，蘆葦中部小種群具有較小的植株及較高的花序生物量分配比，這些都是R對策植物所擁有的典型特征[9]，但是CSR對策分析結果卻顯示中部小種群R對策得分與后沿近陸小種群無顯著差異，這表明基于葉片性狀指標建立的CSR分類方法，可能較難揭示蘆葦種群內水平上繁殖性狀的變異特征。

與CSR體系相比，LHS體系既包含葉片性狀，也包含莖干及繁殖性狀，這使其能夠更好地刻畫種群內功能性狀變異及生態對策差異。在初始的LHS生態對策體系中，Westoby[3]選擇種子質量這一繁殖性狀來表征物種成功定植于其他生境或斑塊的機會。然而，有研究表明，比起種子質量，植株高度才是決定種子傳播距離的關鍵特征[37]，這質疑了種子質量對物種傳播對策的解釋。相比之下，花序干重作為種子和其他傳播結構的質量總和，比種子質量更能代表物種的傳播能力。本研究通過花序干重、比葉面積及株高的性狀組合，成功區分出了3個小種群，說明以花序干重代替種子質量所構建的LHS體系符合蘆葦種群內功能性狀變異特征。此外,在種群內水平下，株高的變異程度要高于葉片化學計量特征性狀及葉片形態學性狀的變異程度[2], 這說明較之以葉片形態學指標為基礎的CSR分類方法，LHS體系可能更能展示出種群內生態對策變異的真實水平。

綜上，崇明東灘濕地3個蘆葦小種群間及各小種群內不同發育階段的分株集合間，都存在較為顯著的功能性狀差異及生態對策差異，且二者的差異模式基本吻合。在該濕地生境內，脅迫因子為塑造蘆葦功能性狀和生態對策的主導因素，而干擾因子和個體發育過程則是蘆葦種群內功能性狀變異和生態對策變異產生的原因。3個蘆葦小種群在21項功能性狀指標及C、R對策得分上均存在顯著差異。在LHS體系中，以比葉面積、株高和花序干重的特定取值為分界線，可對3個小種群進行區分。由SLA、株高及DWinf所構建的LHS體系能夠兼顧葉片、植株整體及繁殖性狀的變異特征，因此，該體系在蘆葦種群內功能性狀及生態對策差異研究方面，可能要優于以葉片性狀為基礎的CSR分類方法。該研究結果深化了我們此前對崇明東灘濕地沿環境梯度分布的蘆葦種群適應機制的認識；基于蘆葦種群內功能性狀及生態對策變異特征，可為該區蘆葦種群的生態保育和恢復實踐，提供重要的理論參考和技術支撐。

[1] Bolnick D I, Amarasekare P, Araújo M S, et al. Why intraspecific trait variation matters in community ecology [J]. Trends Ecol Evol, 2011, 26(4): 183–192. doi: 10.1016/j.tree.2011.01.009.

[2] Siefert A, Violle C, Chalmandrier L, et al. A global meta- analysis of the relative extent of intraspecific trait variation in plant communities [J]. Ecol Lett, 2015, 18(12): 1406–1419. doi: 10.1111/ele. 12508.

[3] Westoby M. A leaf-height-seed (LHS) plant ecology strategy scheme [J]. Plant Soil, 1998, 199(2): 213–227. doi: 10.1023/a:100432 7224729.

[4] Westoby M, Falster D S, Moles A T, et al. Plant ecological strategies: Some leading dimensions of variation between species [J]. Annu Rev Ecol Syst, 2002, 33(1): 125–159. doi: 10.1146/annurev.ecol sys.33.010802.150452.

[5] Astuti G, Ciccarelli D, Roma-Marzio F, et al. Narrow endemic speciesshows significant intraspecific variation in tertiary CSR strategy [J]. Plant Biosyst, 2019, 153(1): 12– 18. doi: 10.1080/11263504.2018.1435576.

[6] De Frenne P, Graae B J, Kolb A, et al. An intraspecific appli- cation of the leaf-height-seed ecology strategy scheme to forest herbs along a latitudinal gradient [J]. Ecography, 2011, 34(1): 132–140. doi: 10.1111/j.1600-0587.2010.06399.x.

[7] Kurze S, Bareither N, Metz J. Phenology, roots and repro- ductive allocation, but not the LHS scheme, shape ecotypes along an aridity gradient [J]. Perspect Plant Ecol Evol Syst, 2017, 29: 20–29. doi: 10.1016/j.ppees.2017.09.004.

[8] Leva?i? D, Jelaska S D. Consistent grime’s CSR strategy of daisy fleabane(L.) Pers. despite its high morphological variability-a case study from Zagreb and Medvednica Mt., Croatia [J]. Diversity, 2022, 14(1): 45. doi: 10.3390/d14010045.

[9] Grime J P. Vegetation classification by reference to strategies [J]. Nature, 1974, 250(5461): 26–31. doi: 10.1038/250026a0.

[10] YAO Z, WU Y M, YOU H L, et al. Research progress of plant competitor-stress tolerator-ruderal (CSR) theory and its application [J]. Acta Ecol Sin, 2022, 42(1): 24–36. [姚忠, 吳永明, 游海林, 等. 競爭-耐脅迫-雜草型植物對策理論及其應用研究進展 [J]. 生態學報, 2022, 42(1): 24–36. doi: 10.5846/stxb202101060055.]

[11] Pierce S, Negreiros D, Cerabolini B E L, et al. A global method for calculating plant CSR ecological strategies applied across biomes world-wide [J]. Funct Ecol, 2017, 31(2): 444–457. doi: 10. 1111/1365-2435.12722.

[12] Pierce S, Brusa G, Vagge I, et al. Allocating CSR plant func- tional types: The use of leaf economics and size traits to classify woody and herbaceous vascular plants [J]. Funct Ecol, 2013, 27(4): 1002– 1010. doi: 10.1111/1365-2435.12095.

[13] Dayrell R L C, Arruda A J, Pierce S, et al. Ontogenetic shifts in plant ecological strategies [J]. Funct Ecol, 2018, 32(12): 2730–2741. doi: 10.1111/1365-2435.13221.

[14] Saltonstall K. The naming ofhaplotypes [J]. Biol Invasions, 2016, 18(9): 2433–2441. doi: 10.1007/s10530-016-1192-4.

[15] An J X, Wang Q, Yang J, et al. Phylogeographic analyses ofin China: Native distribution and habitat prefe- rence of the haplotype that invaded North America [J]. J Syst Evol, 2012, 50(4): 334–340. doi: 10.1111/j.1759-6831.2012.00192.x.

[16] Ji Y H, Zhou G S, Lv G H, et al. Expansion ofin the Liaohe Delta, north-east China [J]. Weed Res, 2009, 49(6): 613– 620. doi: 10.1111/j.1365-3180.2009.00727.x.

[17] Saltonstall K. Cryptic invasion by a non-native genotype of the common reed,, into North America [J]. Proc Natl Acad Sci USA, 2002, 99(4): 2445–2449. doi: 10.1073/pnas.032477999.

[18] Lelong B, Lavoie C, Jodoin Y, et al. Expansion pathways of the exotic common reed (): A historical and genetic analysis [J]. Divers Distrib, 2007, 13(4): 430–437. doi: 10.1111/j.1472- 4642.2007.00351.x.

[19] Saltonstall K, Meyerson L A.: From genes to ecosystems [J]. Biol Invasions, 2016, 18(9): 2415–2420. doi: 10.1007/s10530-016-1240-0.

[20] Gan X J, Cai Y T, Choi C, et al. Potential impacts of invasiveon spring bird communities at Chongming Dongtan, a Chinese wetland of international importance [J]. Estuar Coast Shelf Sci, 2009, 83(2): 211–218. doi: 10.1016/j.ecss.2009.03.026.

[21] YANG Y T, ZENG H W, WEN Z M, et al. Classification of plant functional types based on LHS strategy in Yanhe River Basin [J]. J NW For Univ, 2019, 34(2): 84–91. [楊玉婷, 曾鴻文, 溫仲明, 等. 基于LHS策略的延河流域植物功能型劃分 [J]. 西北林學院學報, 2019, 34(2): 84–91. doi: 10.3969/j.issn.1001-7461.2019.02.13.]

[22] XU P C, YOU Z T, JI Y H, et al. Study on ecological strategies of 22 common woody plants inNature Reserve [J]. J Subtrop Resour Environ, 2019, 14(4): 23–29. [徐鵬程, 游章湉, 紀宇皝, 等. 格氏栲自然保護區22種常見木本植物的生態策略 [J]. 亞熱帶資源與環境學報, 2019, 14(4): 23–29. doi: 10.19687/j.cnki.1673- 7105.2019.04.004.]

[23] Zhou D y, Ni Y H, Yu X N, et al. Trait-based adaptability ofto the effects of soil water and salinity in the Yellow River Delta [J]. Ecol Evol, 2021, 11(16): 11352–11361. doi: 10.1002/ece3.7925.

[24] LI J, NIE J Y, LI H F, et al. On determination conditions for total polyphenols in fruits and its derived products by Folin-phenol methods [J]. J Fruit Sci, 2008, 25(1): 126–131. [李靜, 聶繼云, 李海飛, 等. Folin-酚法測定水果及其制品中總多酚含量的條件 [J]. 果樹學報, 2008, 25(1): 126–131. doi: 10.3969/j.issn.1009-9980.2008.01.026.]

[25] Kapali J, Sharma K R. Estimation of phytochemicals, antioxidant, antidiabetic and brine shrimp lethality activities of some medicinal plants growing in Nepal [J]. J Med Plants, 2021, 20(80): 102–116. doi: 10.52547/jmp.20.80.102.

[26] Tong Y, Gao L J, Xiao G M, et al. Microwave pretreatment- assisted ethanol extraction of chlorophylls from[J]. J Food Process Eng, 2012, 35(5): 792–799. doi: 10.1111/j.1745-4530. 2010.00629.x.

[27] CUI L L, LI G S, OUYANG N L, et al. Response of variation of net primary productivity to edaphic factors in the coastal salt marshes [J]. Acta Ecol Sin, 2020, 40(19): 7018–7029. [崔林林, 李國勝, 歐陽寧雷,等. 濱海鹽沼濕地植被凈初級生產力變化對土壤因子的響應 [J]. 生態學報, 2020, 40(19): 7018–7029. doi: 10.5846/stxb201905301134.]

[28] Ben Dkhil B, Denden M. Salt stress induced changes in germi- nation, sugars, starch and enzyme of carbohydrate metabolism inL. (Moench.) seeds [J]. Afr J Agric Res, 2010, 5(12): 1412–1418.

[29] Malone J G, Mittova V, Ratcliffe R G, et al. The response of carbohydrate metabolism in potato tubers to low temperature [J]. Plant Cell Physiol, 2006, 47(9): 1309–1322. doi: 10.1093/pcp/pcj101.

[30] Schmitz J, Heinrichs L, Scossa F, et al. The essential role of sugar metabolism in the acclimation response ofto high light intensities [J]. J Exp Bot, 2014, 65(6): 1619–1636. doi: 10. 1093/jxb/eru027.

[31] DING W H, JIANG J Y, LI X Z, eta al. Spatial distribution of species and influencing factors across salt marsh in southern Chongming Dongtan [J]. Acta Phytoecol Sin, 2015, 39(7): 704–716. [丁文慧, 姜俊彥, 李秀珍, 等. 崇明東灘南部鹽沼植被空間分布及影響因素分析 [J]. 植物生態學報, 2015, 39(7): 704–716. doi: 10.17521/cjpe.2015. 0067.]

[32] Shaar-Moshe L, Hayouka R, Roessner U, et al. Phenotypic and metabolic plasticity shapes life-history strategies under combinations of abiotic stresses [J]. Plant Direct, 2019, 3(1): e00113. doi: 10.1002/ pld3.113.

[33] Feeny P. Seasonal changes in oak leaf tannins and nutrients as a cause of spring feeding by winter moth caterpillars [J]. Ecology, 1970, 51(4): 565–581. doi: 10.2307/1934037.

[34] Forkner R E, Marquis R J, Lill J T. Feeny revisited: condensed tannins as anti-herbivore defences in leaf-chewing herbivore commu- nities of[J]. Ecol Entomol, 2004, 29(2): 174–187. doi: 10. 1111/j.1365-2311.2004.0590.x.

[35] Lavorel S, Mcintyre S, Landsberg J, et al. Plant functional classifications: From general groups to specific groups based on response to disturbance [J]. Trends Ecol Evol, 1997, 12(12): 474–478. doi: 10.1016/s0169-5347(97)01219-6.

[36] Pierce S, Bottinelli A, Bassani I, et al. How well do seed production traits correlate with leaf traits, whole-plant traits and plant ecological strategies? [J]. Plant Ecol, 2014, 215(11): 1351–1359. doi: 10.1007/s11258-014-0392-1.

[37] Thomson F J, Moles A T, Auld T D, et al. Seed dispersal distance is more strongly correlated with plant height than with seed mass [J]. J Ecol, 2011, 99(6): 1299–1307. doi: 10.1111/j.1365-2745. 2011.01867.x.

Analysis on Intra-population Variations of Functional Traits and Ecological Strategies ofin Chongming Dongtan Wetland Based on CSR and LHS System

JI Yangqing1, LI Dezhi1a,1b,2,3*, WANG Xiao1a, LI Lingling1a, WANG Ying1a, CHEN Jing1a, Rasool Samreen Ghulam1a, Tehreem Rafique1a

(1a. School of Ecological and Environmental Sciences; 1b. Key Laboratory of Urbanization and Ecological Restoration of Shanghai, East China Normal University,Shanghai 200241, China; 2. Institute of Eco-Chongming (IEC), Shanghai 202162, China; 3. Technology Innovation Center for Land Spatial Eco-restoration in Metropolitan Area, Ministry of Natural Resources,Shanghai 200062, China)

The variations of functional traits and ecological strategies of population may greatly affect its adaptive capacity and the structure of a community. Based on the multiple functional traits of the ramet groups at the different developmental stages of 3 subpopulations ofin Chongming Dongtan wetland, this paper analyzed its intra-population variations of functional traits and ecological strategies, and made an ordination of the ecological strategies of each subpopulation based on the CSR and LHS schemes. The results showed that 21 functional traits, C and R strategy scores among the 3 subpopulations ofwere significantly different. In the LHS scheme, 3 subpopulations were separated along the boundary line of the specific leaf area (11.5 mm2/mg), plant height (175 cm) and inflorescence dry weight (4 g). Ramets at the reproductive growth stage were characterized by larger plant sizes, higher leaf dry matter contents and S scores, while those at the vegetative growth stage, the smaller plant sizes, larger specific leaf areas and R scores. In general, there were significant differences in functional traits and ecological strategies among 3 subpopulations and between the ramet groups at the different developmental stages, and the divergent pattern of functional traits and ecological strategies were basically consistent. These might provide important theoretical reference and technical support for the ecological conservation and restoration ofpopulation in this area.

Intraspecific trait variation; Ecological strategy; CSR scheme; LHS scheme;; Wetland

10.11926/jtsb.4783

2023-03-01

2023-03-16

上海市科委國際科技合作一帶一路戰略青年科學家交流項目(19230742600)；自然資源部大都市區國土空間生態修復工程技術創新中心開放研究項目(CXZX2021B01, CXZX2022B01)；國家重點研發計劃(2017 YFC 0506002, 2016 YFC 0503102)資助

This work was supported by the Project for Young Scientist Exchange of One Belt and One Road Strategy in the International Science and Technology Cooperation of Shanghai Science and Technology Commission (Grant No. 19230742600), the Open Research Project for the Technology Innovation Center for Land Spatial Eco-restoration in Metropolitan Area, Ministry of Natural Resources (Grant No. CXZX2021B01, CXZX2022B01), and the National Key Research & Development Program of China (Grant No. 2017 YFC 0506002, 2016 YFC 0503102).

計陽晴(1998年生)，女，在讀碩士研究生，主要從事植物生態學研究。E-mail: 1213266277@qq.com

通訊作者 Corresponding author.E-mail: dzli@des.ecnu.edu.cn