摘 要：堡狀峰是類似城堡外貌而凸出于周圍環(huán)境之上的特殊山峰，位于中國(guó)閩南烏山風(fēng)景區(qū)堡狀峰地貌上的烏岡櫟（殼斗科）種群的表型發(fā)生了變異，其與原類型的顯著不同在于葉背具宿存性黃褐色星狀毛。為解釋這一現(xiàn)象，該研究對(duì)堡狀峰頂部的烏岡櫟種群分布進(jìn)行調(diào)查，包括其中的變異型、原始型、過(guò)渡型，以及堡狀峰山腳、另一異域的典型原始型種群被研究，基于葉表皮解剖，利用方差分析統(tǒng)計(jì)顯微性狀的差異，結(jié)合生態(tài)因子，探究變異的原因及演化的可能途徑。結(jié)果表明：（1）變異型在山頂西側(cè)形成以其為優(yōu)勢(shì)種的常綠闊葉灌叢群落，東側(cè)的小片灌叢以原始型為主，后者內(nèi)部還有少量過(guò)渡型。（2）各種群?jiǎn)渭?xì)胞毛狀體基（STB）與復(fù)合毛狀體基（CTB）在葉上表皮的密度很小甚至小于每平方毫米1個(gè)，宏觀的表現(xiàn)為無(wú)毛。原始型、過(guò)渡型、變異型的STB密度均值在葉下表皮遞增。變異型種群與原始種群相比，葉下表皮CTB或典型星狀毛（Tst）的密度由每平方毫米0.481～1.122個(gè)變異至每平方毫米57.2個(gè)，宏觀表現(xiàn)為形成密毛。氣孔長(zhǎng)寬比在原始型（1.187～1.205）至變異型（1.258）間遞增。（3）方差分析結(jié)果顯示，變異型葉下表皮STB與CTB密度及氣孔長(zhǎng)寬比較原始型顯著增加（大）（Plt;0.01），但氣孔密度與尺寸在種群間未表現(xiàn)出規(guī)律與差異性（P＞0.05）。烏山堡狀峰頂?shù)耐寥镭汃ぁ⑺疅釛l件差。綜合生境與統(tǒng)計(jì)的結(jié)果表明，變異型烏岡櫟葉背毛狀體及氣孔長(zhǎng)寬比的變異具抵抗干旱或極端氣溫的功能。由于花粉流未曾間斷，因此相比于地理隔離，特殊生境促使烏岡櫟發(fā)生適應(yīng)性演化應(yīng)是該變異產(chǎn)生的主要原因。該研究為亞熱帶特殊花崗巖地貌促進(jìn)種群分化和物種形成提供了直接的證據(jù)。

關(guān)鍵詞：陸地島嶼系統(tǒng)， 烏岡櫟， 葉顯微性狀， 自然選擇， 種群分化

中圖分類號(hào)：Q944" "文獻(xiàn)標(biāo)識(shí)碼：A" "文章編號(hào)：1000-3142（2024）03-0488-13

Phenotypic variation of Quercus phillyreoides populations on castle peak of South Fujian， China LIN Yizhe1， ZHU Yiyao2， FU Houhua1， MA Liang3， CHEN Shipin1*

（ 1. College of Forestry， Fujian Agriculture and Forestry University， Fuzhou 350002， China; 2. Yunxiao County Forestry Bureau，Zhangzhou 363399， Fujian， China; 3. Fujian Health College， Fuzhou 350101， China ）

Abstract： "Castle peak is a mountain that resembles the appearance of a castle and protrudes above its surroundings. The phenotypes of Quercus phillyreoides （Fagaceae） populations located on castle peak in the Wushan Scenic Spot of southern Fujian， China， are variegated. It differs significantly from the original form in the presence of persistent yellow-brown stellate hairs on the leaf abaxial surface. To explain this phenomenon， the distribution of Q. phillyreoides populations at the summit of castle peak was investigated， where variant， original， and intermediate forms， as well as typical original populations at the foot of castle peak and another allopatric area， were studied. Differences in micromorphological traits were counted using ANOVA based on leaf epidermal anatomy. Ecological factors were linked with a view to discovering the causes of variation and possible pathways of evolution. The results were as follows： （1） The variant form formed an evergreen broad-leaved scrub community with it as the dominant species on the west side of the summit， and the small patch of scrub on the east side is dominated by the original form， with a few intermediate form individuals within the latter. （2） The densities of single trichome base （STB） and compound trichome base （CTB） in the leaf adaxial epidermis of various populations was very low or even less than 1 No.·mm-2， and the macroscopic expression was glabrous. The mean values of STB density of original， intermediate， and variant forms were increasing in the leaf abaxial epidermis. The mean density of CTB or typical stellate trichome （Tst） on the leaf abaxial epidermis varied from 0.481-1.122 No.·mm-2 to 57.2 No.·mm-2 in the variant form population compared with the original form populations， which [JP3]macroscopically showed the formation of dense hairs. The stomatal length-width ratio increases from original form （1.187-[JP2]1.205） to variant form （1.258）. （3） ANOVA results showed that variant form leaf abaxial epidermal STB， CTB density and stomatal length-width ratio were significantly increased （larger） than the original form （Plt;0.01）， but stomatal density and size did not show regularity and variability among populations （Pgt;0.05）. Poor soil， water and heat conditions at the summit of the Wushan castle peak， and combining the background with the obtained results indicate that" the variation in the variant Q. phillyreoides leaf abaxial trichomes and stomatal aspect ratio were resistant to drought or temperature extremes. Since pollen flow has not been interrupted， the adaptive evolution of Q. phillyreoides due to specific habitats should be the main reason for this variation compared to geographical isolation. The microscopic traits showed no significant differences （P gt; 0.05） between the original form population on the top of castle peak and the rest original form populations confirms that the variation has been fixed by long-term selection and the short-term environmental influences play a limited role in shaping the variant traits. This study provides direct evidence that subtropical exceptional granitic landforms promote population differentiation and species formation.

Key words： continental island system， Quercus phillyreoides， leaf microscopic character， natural selection， population differentiation

花崗巖堡狀峰（castle peak）（下稱堡狀峰）是四周由陡直崖壁圍限呈城堡狀的一種花崗巖地貌類型（梁詩(shī)經(jīng)和文斐成，2010），可連綿成墻狀山體。中國(guó)東南沿海的堡狀峰以新近紀(jì)（Neogene）后地質(zhì)構(gòu)造發(fā)生中度抬升為基礎(chǔ)，經(jīng)化學(xué)風(fēng)化、流水下切侵蝕塑造而來(lái)（崔之久等，2007）。堡狀峰實(shí)際為介于石林（stone forest）與石蛋（pebbly）之間的一種花崗巖過(guò)渡地貌類型。其特點(diǎn)為四周陡峻，巖壁裸露，峰頂平坦，可沉積土壤。較周邊丘陵，堡狀峰突兀地佇于山間，生態(tài)特征與低地截然不同（Porembski， 2007）。

熱帶地區(qū)的另一種花崗巖地貌“島山”（inselberg）與堡狀峰相似，但區(qū)別在于前者構(gòu)造抬升輕微，經(jīng)化學(xué)風(fēng)化在平原上形成孤立的巖石露頭（崔之久等，2007）。島山上的植物群落常與周邊地區(qū)不同，被認(rèn)為是孤立山體特殊生境產(chǎn)生的效果；同時(shí)，其孤立程度取決于周圍植被類型，較于真實(shí)島嶼，生態(tài)隔離彌補(bǔ)了島山之間地理距離相對(duì)較短的不足（Porembski， 2007）。生態(tài)與距離隔離以及不同于山腳的特殊生境共同作用于島山的物種演化。就物種多樣性而言，有研究顯示其他因素（如面積、環(huán)境篩）而非隔離距離與物種豐富度間有更大關(guān)聯(lián)（Conceio et al.， 2007；Henneron et al.， 2019），表明在較小空間尺度內(nèi)，在異質(zhì)性生境下發(fā)生了適應(yīng)性輻射。大量研究表明，一些植物種群形態(tài)與遺傳結(jié)構(gòu)在或遠(yuǎn)或近的島山上發(fā)生變異，基因流發(fā)生不同程度的限制，種群內(nèi)遺傳多樣性降低，促進(jìn)了種群分化與物種形成（Barbará et al.， 2008；Byrne et al.， 2019；Feliciano et al.， 2022）。

島山作為陸地島嶼系統(tǒng)用于觀察生態(tài)現(xiàn)象與物種形成，在溫帶與熱帶受到長(zhǎng)期的關(guān)注（Porembski amp; Barthlott， 2012）。中國(guó)擁有不同的地質(zhì)史，喜馬拉雅山脈、阿爾泰-天山地區(qū)、中國(guó)西南山地這類稱作“天空之島”（空島，sky islands）的陸地島嶼近年的關(guān)注度較高，其中地形地貌及其改變與隨之產(chǎn)生的氣候變遷對(duì)譜系地理結(jié)構(gòu)的影響研究較多（He amp; Jiang， 2014）。西南地區(qū)的盆地、峽谷在更新世冰川期為物種提供相互間隔的小生境避難所，這些避難所便成為遺傳分化的中心（Xie et al.， 2012；Yang ZY et al.， 2012），青藏高原的抬升還讓異地分化的物種二次相遇發(fā)生偶然的雜交種化事件（Yang FS et al.， 2012；Zhang et al.， 2019），現(xiàn)有證據(jù)表明空島上的高寒植物共享了一連串地質(zhì)運(yùn)動(dòng)引發(fā)的譜系地理結(jié)構(gòu)與進(jìn)化歷史（Luo et al.， 2016）。然而，在中國(guó)亞熱帶的花崗巖堡狀峰、峰叢、峰墻、石林等地貌上開(kāi)展有關(guān)進(jìn)化與生態(tài)的研究甚為少見(jiàn)。

狹義櫟屬（Quercus）全世界約300種，中國(guó)約51種、14變種、1變型（中國(guó)科學(xué)院中國(guó)植物志編輯委員會(huì)，1998），由于復(fù)雜的網(wǎng)狀進(jìn)化歷史，具有極高物種多樣性（Hipp et al.， 2020），櫟屬植物受地理與古氣候的影響在各地發(fā)生適應(yīng)性演化（Yang et al.， 2021），葉表皮性狀常視為重要特征之一被用于種間界定（Deng et al.， 2017）。2019年11月，我們?cè)诟＝ㄊ∧喜空闹菔性葡龊驮t安兩縣交界的烏山風(fēng)景區(qū)開(kāi)展全國(guó)中藥資源普查（圖1：A， C），于烏山主峰紅旗巖（Hongqiyan）頂發(fā)現(xiàn)一葉柄、葉背顯著被星狀毛的烏岡櫟（Quercus phillyreoides）群體（其余性狀未表現(xiàn)出顯著不同）。[JP]經(jīng)過(guò)1年的持續(xù)觀察，發(fā)現(xiàn)其成熟葉片毛被不脫落。烏岡櫟是中國(guó)秦嶺淮河以南廣泛分布的殼斗科（Fagaceae）櫟屬冬青礫組（Section Ilex）常綠闊葉樹(shù)，日本也有分布，新葉有時(shí)疏被毛，老葉無(wú)毛或僅葉背中脈被疏柔毛（中國(guó)科學(xué)院中國(guó)植物志編輯委員會(huì)，1998）。[JP]對(duì)此，我們推測(cè)該種群葉表皮特征或多或少發(fā)生變異，故對(duì)烏山風(fēng)景區(qū)堡狀峰山頂、山腳及另一異域的烏岡櫟種群做顯微形態(tài)的比較研究，通過(guò)葉解剖獲取表皮性狀，所得性狀數(shù)據(jù)經(jīng)方差分析檢驗(yàn)與比較，擬探討以下問(wèn)題：（1）變異種群哪些葉表皮性狀發(fā)生了變異，變異的規(guī)律與趨勢(shì)如何；（2）變異產(chǎn)生與堡狀峰生境的關(guān)系，以及其是否具備生態(tài)合理性；（3）種群分化的歷史進(jìn)程，以及其是否具備演化合理性。通過(guò)對(duì)以上問(wèn)題的探討，以期為亞熱帶特殊花崗巖地貌促進(jìn)種群分化提供證據(jù)。

1 材料與方法

1.1 研究樣品

采集樣品之前，我們查閱了烏岡櫟模式標(biāo)本及中國(guó)數(shù)字標(biāo)本館（Chinese Virtual Herbarium， CVH， https：//www.cvh.ac.cn/）的所有普通標(biāo)本（有照片的標(biāo)本），2份模式標(biāo)本（Cavalerie， J. 3626）均為葉無(wú)毛被類型，209份普通標(biāo)本中有1份采自廣西的標(biāo)本與我們?cè)诩t旗巖頂所見(jiàn)相似（葉背顯著具毛，陳照宙50295），采集簽記錄的生境、生活型亦類同。然而，中國(guó)分布之北界未見(jiàn)該性狀出現(xiàn)，表明葉背具穩(wěn)定毛被的性狀不是環(huán)境直接作用的結(jié)果，亦非烏岡櫟表型可塑性可簡(jiǎn)單解釋。與此同時(shí)，我們?cè)诩t旗巖頂不但發(fā)現(xiàn)葉背具毛的種群，典型原始型亦有分布。

為了解釋這些現(xiàn)象，我們采集了紅旗巖頂原始型、變異型與假定的過(guò)渡型，以及紅旗巖腳原始型的葉片與憑證標(biāo)本。由表1可知，野外觀察到的過(guò)渡型是相對(duì)原始型葉背的幾乎完全無(wú)毛以及變異型的葉背幾乎完全為毛狀體覆蓋（葉下表皮不可見(jiàn)）而言，即過(guò)渡型雖然擁有宿存性的毛，但其葉下表皮在全株均明顯可見(jiàn)，介于原始型和變異型兩者之間。由于原始型的穩(wěn)定性，允許我們以任意異域環(huán)境的烏岡櫟種群代表全生境范圍的所有烏岡櫟群體，來(lái)說(shuō)明紅旗巖頂?shù)倪@種變異。5個(gè)（極小）種群被研究，每種群隨機(jī)選擇3個(gè)成年個(gè)體（個(gè)體間大于等于20 m或最大間隔距離），每個(gè)體東西南北重復(fù)取樣（共4片）；考慮到毛狀體對(duì)氣孔性狀的遮[JP]蓋將不利于后期觀察與統(tǒng)計(jì)，對(duì)變異型每個(gè)方向額外采集1片葉子，利用解剖刀去除毛被。葉片被切割成葉尖、葉中、葉基3段，保存于FAA固定液（38%福爾馬林∶冰醋酸∶70%乙醇=1∶1∶18，體積比）中。樣品采集在2 d內(nèi)完成，采集時(shí)間和天氣均相同。葉表皮切片的制作方法參考Deng等（2017）。

1.2 研究地概況

漳州烏山主體位于福建南部與廣東交界的詔安、云霄兩縣境內(nèi)，為東南海岸平原、東北漳江、西南東溪環(huán)繞；南北走向，南北約45 km，東西約20 km，面積約900 km2。烏山風(fēng)景區(qū)位于烏山中部，主峰海拔1 050.6 m，山體四周懸崖絕壁，稍有棱角，頂部平坦，與附近山谷高差達(dá)215 m（圖1：A）。烏山基巖為花崗巖， 新近紀(jì)后局部構(gòu)造中度抬升山谷，上坡，青云山峽谷瀑布風(fēng)景區(qū)，永泰，福建，中國(guó)。119°0′28.39″ E， 25°45′40.01″ N，海拔818 m，2021年12月9日Valley， uphill， Qingyun Mountain Canyon Waterfall Scenic Spot， Yongtai， Fujian， China. 119°0′28.39″ E， 25°45′40.01″ N， alt. 818 m， 9 Dec. 2021（圖1：B）。烏山主峰紅旗巖屬幼年（juvenile stage）至青年（adolescence stage）發(fā)育期花崗巖地貌（Chen et al.， 2009），即為典型的堡狀峰（圖1：C）。山頂保留中新世統(tǒng)一夷平面的殘留部分，山腳堆積深厚風(fēng)化殼，發(fā)育成紅壤（崔之久等，2007）。亞熱帶季風(fēng)氣候，年均降水1 519.6 mm，夏季均降水624 mm，冬季均降水121 mm；白天均溫27 ℃，夜間均溫20" ℃，夏季最高溫31 ℃，冬季最低溫12 ℃。地帶植被為常綠闊葉林，殼斗科、樟科（Lauraceae）、山茶科（Theaceae）等為該地帶森林的優(yōu)勢(shì)類群（吳中倫和蔣有緒，1997）。

1.3 統(tǒng)計(jì)和分析方法

利用BK5000型光學(xué)顯微鏡觀察葉上下表皮細(xì)胞垂周壁、毛狀體基、毛狀體及氣孔。葉表皮術(shù)語(yǔ)及特征描述參考Deng等（2017）。每切片選1個(gè)視野（200×）拍攝記錄。用ImageJ 1.53f51 （https：//imagej.nih.gov/ij/index.html）計(jì)測(cè)視野內(nèi)全部氣孔的數(shù)目，并測(cè)量尺寸（保衛(wèi)細(xì)胞外輪廓的長(zhǎng)徑及短徑）（Arena et al.， 2017）。毛狀體/毛狀體基數(shù)目由人工在每個(gè)視野中計(jì)數(shù)得到。隨后氣孔、毛狀體/毛狀體基數(shù)目轉(zhuǎn)化為密度，氣孔尺寸計(jì)算出長(zhǎng)寬比（長(zhǎng)徑與短徑的比值）。將每片葉子葉尖、葉中、葉基各性狀的數(shù)據(jù)求平均，以排除隨機(jī)誤差。由于4個(gè)方向的測(cè)量值并無(wú)明顯的規(guī)律，因此將其全部納入原始數(shù)據(jù)集，以減小抽樣個(gè)體數(shù)較少可能存在的系統(tǒng)誤差。用SPSS 25對(duì)原始數(shù)據(jù)集各性狀實(shí)施F檢驗(yàn)與多重比較（LSD），對(duì)結(jié)果具顯著性的性狀繪制箱型圖。Plt;0.01，認(rèn)為差異極顯著；Plt;0.05，則認(rèn)為差異顯著。

2 結(jié)果與分析

2.1 表型變異和空間分布

變異種群生于紅旗巖頂?shù)幕◢弾r風(fēng)化殼上（圖2：A-D），在山頂西側(cè)大片生長(zhǎng)，集群分布明顯，形成以其為優(yōu)勢(shì)種的常綠闊葉灌叢群落，伴生種有鳳凰潤(rùn)楠（Machilus phoenicis）、漳州越橘（Vaccinium zhangzhouense）、猴頭杜鵑（Rhododendron simiarum）[JP+1]等。東側(cè)的小片灌叢以原始型為主，數(shù)量較少；而在變異種群外圍亦見(jiàn)零星分布。北側(cè)狹長(zhǎng)、地勢(shì)較低地帶的變異型可能由山頂西擴(kuò)散而來(lái)。在原始種群內(nèi)部，發(fā)現(xiàn)葉背毛稍顯發(fā)育的個(gè)體，視其為過(guò)渡型。

唯獨(dú)葉背密生成熟時(shí)不脫落的土黃色或黃褐色星狀毛，一年生枝的星狀毛也較密集而與原類型不同（圖2：F-G）。[JP+1]紅旗巖山頂烏岡櫟變異型的花序、殼斗、堅(jiān)果（圖2：H）及其他葉性狀（形狀、大小、質(zhì)地、顏色及葉緣等）均與原類型無(wú)異。

2.2 葉表皮性狀和顯著性結(jié)果

2.2.1 葉上表皮（圖3：A-E） 葉上表皮細(xì)胞垂周壁厚度均勻，有4種形態(tài)，即平直、弧形、波狀、蜿蜒（圖3：A， D-E），在各類各地種群間無(wú)規(guī)律。除在紅旗巖腳原始型2（種群4）未在視野中觀察到單細(xì)胞毛狀體基（STB）以外，其余或多或少存在單細(xì)胞（圖3：A）及復(fù)合毛狀體基（CTB）（圖3：B，C， E）2種類型，但CTB和STB密度均值小于每平方毫米7個(gè)；在青云山原始型3（種群5）觀察到來(lái)源于STB的殘存的頭狀毛狀體（Ca），然而，均未發(fā)現(xiàn)與CTB相對(duì)應(yīng)的毛狀體類型。綜上表明，各類型種群STB與CTB在葉上表皮的密度很小甚至小于每平方毫米1個(gè)，然而即便存在毛狀體基，毛狀體也不一定發(fā)生或容易脫落，在宏觀上的表現(xiàn)即為無(wú)毛。特征統(tǒng)計(jì)見(jiàn)表2。

2.2.2 葉下表皮（圖3：F-L） 葉下表皮細(xì)胞垂周壁在各種群內(nèi)多少表現(xiàn)出平直、弧形、波狀之漸變， 厚度均勻（圖3：H， J）。全部種群均觀察到明

[JP]顯的頭狀毛狀體（Ca）及STB（圖3：I-L），后者的密度為每平方毫米0～10.125個(gè)。由表3可知，原始型、過(guò)渡型、變異型的STB密度均值遞增。雖然CTB各種群均存在（圖3：F， H-I， K），但是由其發(fā)生的典型星狀毛（Tst）僅見(jiàn)于紅旗巖頂變異型（種群1）、過(guò)渡型（種群2）及青云山原始型3（種群5）[HJ1.84mm]（圖3：F-G， I， L）。然而，青云山原始型3的葉下表皮星狀毛只局限于中脈的基部附近，對(duì)應(yīng)的CTB密度均值則小于每平方毫米1個(gè)；相比之，紅旗巖頂變異型的Tst存在于中脈、任意側(cè)脈或次級(jí)脈上，CTB密度均值達(dá)每平方毫米57.2個(gè)，中脈、側(cè)脈、次級(jí)脈上的CTB面積逐級(jí)減??；過(guò)渡型的CTB密度均值介于兩者之間，為每平方毫米26.076個(gè)。這表明變異型種群與其他種群相比，CTB或Tst的增加很明顯，宏觀表現(xiàn)為密被星狀毛，過(guò)渡型的葉背雖有星狀毛，但葉下表皮也明顯可見(jiàn)，原始型則幾乎完全無(wú)毛。另外，氣孔類型均為無(wú)規(guī)則型（圖3：H， J）。各種群氣孔密度均值在每平方毫米453～516個(gè)之間，氣孔尺寸之短徑為24.494～25.812 μm，長(zhǎng)徑介于28.906～31.047 μm，變異型與原始型在這3個(gè)性狀的變化無(wú)規(guī)律。氣孔長(zhǎng)寬比均值則在原始型（1.187～1.205）、過(guò)渡型（1.240）、變異型（1.258）間遞增。相關(guān)特征統(tǒng)計(jì)見(jiàn)表3。

2.2.3 顯著性檢驗(yàn) 方差分析結(jié)果（表2， 表3）顯示，葉上表皮STB（P=0.002 0）、葉下表皮STB（P=0.006 0）與CTB（P=0.000 5）密度及氣孔長(zhǎng)寬比（P=0.002 0）在不同種群之間存在極顯著性差異（Plt;0.01）。葉上表皮STB密度雖然不同種群差異極顯著，但經(jīng)事后檢驗(yàn)發(fā)現(xiàn)變異型與原始型2、3，原始型1與原始型2、3共同貢獻(xiàn)了該差異，表明該性狀在葉上表皮的發(fā)生較隨機(jī)。葉下表皮的事后檢驗(yàn)結(jié)果顯示葉下表皮STB與CTB密度及氣孔長(zhǎng)寬比的差異主要是由變異型引起（圖4）。變異型STB與CTB密度顯著大于原始型，無(wú)論原始型位于山頂、山腳或異地。變異型氣孔平均長(zhǎng)寬比為1.258，極大值1.333；原始型均值不大于1.205，極大值1.274。這一性狀在變異型中在一定程度上有增大，盡管在影像中看上去并不明顯。過(guò)渡型各性狀值均居中間水平。以上結(jié)果表明，變異型在一些葉下表皮性狀（STB與CTB密度、氣孔長(zhǎng)寬比）上與原始型相比發(fā)生了顯著性變異（Plt;0.01）。但是，對(duì)生境頗為敏感的氣孔密度與尺寸在山頂、山腳或異地的種群間未表現(xiàn)出差異性（P＞0.05）（表3）。

3 討論與結(jié)論

3.1 堡狀峰生境促進(jìn)了種群分化

同為花崗巖地貌的島山導(dǎo)致種群變異和物種形成的案例已經(jīng)屢見(jiàn)不鮮（Mota et al.， 2020；Ruas et al.， 2020；Tavares et al.， 2022），這通常歸因于距離隔離導(dǎo)致有限的基因流發(fā)生的遺傳漂變和近交（Barbará et al.， 2007；Mota et al.， 2020），然而有時(shí)異質(zhì)性生境被認(rèn)為起到了更大作用， 驅(qū)使同一物種在不同種群間發(fā)生分化（Milá et al.， 2010；Mallet et al.， 2014）。雖然島山生境通常不適合物種的生存，但是世界上一些地區(qū)的島山卻成為生物多樣性的中心，如巴西東南部、馬達(dá)加斯加、澳大利亞西南部，推測(cè)與隔離及物種長(zhǎng)期的堅(jiān)持促進(jìn)分化有關(guān)（Porembski， 2007），因此島山多特有種。

烏山紅旗巖是典型堡狀峰地貌，與島山類似，頂部生境與同緯南亞熱帶季風(fēng)氣候區(qū)迥異，與周邊地區(qū)相比，具有一定的特殊性：夏季光照強(qiáng)度大，日照時(shí)間長(zhǎng)；冬季夜間氣溫低；氣溫日較差大；降水滯留土壤時(shí)間短，雨水沖刷土壤造成流失；土壤本身貧瘠而??；盛行風(fēng)速大。一般情況下，同一物種在一定程度的干燥與適生條件下生長(zhǎng)相比氣孔密度增大，但當(dāng)干旱進(jìn)一步加強(qiáng)，密度將減小，氣孔尺寸則在干旱條件下持續(xù)減?。╔u amp; Zhou， 2008）。然而，烏山不同生境下的烏岡櫟并未表現(xiàn)出這一傾向，并且在更大空間尺度的櫟屬植物研究案例中也存在這一現(xiàn)象（Sánchez-Acevedo et al.， 2022）。由于氣孔密度與尺寸通常具有一定關(guān)聯(lián)性（Hetherington amp; Woodward， 2003），因此不同種群烏岡櫟的氣孔尺寸與密度一樣無(wú)顯著性差異。

為了應(yīng)對(duì)烏山山頂特殊生境的挑戰(zhàn)，在我們研究的中烏岡櫟似乎傾向于依賴毛狀體并調(diào)整氣孔的形狀，因?yàn)榻Y(jié)果顯示葉下表皮單細(xì)胞毛狀體基（STB）、復(fù)合毛狀體基（CTB）及氣孔長(zhǎng)寬比在不同種群間差異顯著。毛狀體有吸附空氣中水滴、減少蒸騰的功能（Moles et al.， 2020），在冬季具有一定保暖效果（Peng et al.， 2015）。就荒漠植物而言，在應(yīng)對(duì)寒、旱、鹽堿、高溫、強(qiáng)輻射等不同非生物脅迫時(shí)表皮毛的變化不一，在受機(jī)械損傷、干旱、長(zhǎng)日照、高光照、鹽脅迫時(shí)，通常會(huì)導(dǎo)致表皮毛長(zhǎng)度或密度的增大 （馬亞麗等， 2015）。烏山頂變異型烏岡櫟的葉背STB、CTB密度與各原始型相比均顯著，盡管毛狀體在極端環(huán)境中可能被催生，但山頂原始型1的存在，以及與原始型2、3之間無(wú)顯著差異，證實(shí)環(huán)境的短期作用對(duì)塑造烏岡櫟葉背毛狀體的能力有限。顯著增多葉背的星狀毛（Tst） 勢(shì)必有助其適應(yīng)山頂貧瘠土[JP+2]壤儲(chǔ)水少、高蒸散導(dǎo)致的干旱，以及冬季花崗巖山體比熱小、散熱快導(dǎo)致的冷害、凍害。頭狀毛狀體具備合成、貯藏和分泌諸多代謝物的作用（張繼偉等，2018），增多的頭狀毛或能增加烏岡櫟應(yīng)對(duì)山頂環(huán)境脅迫的抗性。關(guān)于環(huán)境對(duì)毛狀體的作用機(jī)制，就擬南芥（Arabidopsis thaliana）的研究顯示，擁有豐富表皮毛的突變體的較于表皮毛缺失突變體其莖、葉蠟質(zhì)層含有的C32+脂肪酸混合物含量顯著增加，并且表皮毛所覆蓋的蠟質(zhì)C32+脂肪酸含量亦被證實(shí)高于周圍細(xì)胞，表皮毛可聯(lián)合蠟質(zhì)層生物合成來(lái)利于減少增發(fā)和調(diào)控溫度（Hegebarth et al.， 2016）。對(duì)黃瓜（Cucumis sativus）表皮毛突變體的研究表明，GL1與TTG1為葉表皮發(fā)生的核心基因，GL3、EGL3為必要基因，這些基因的轉(zhuǎn)錄產(chǎn)物構(gòu)成正調(diào)控因子（Liu et al.， 2016）。但是，在TRY、CPC等基因轉(zhuǎn)錄的負(fù)調(diào)控因子介導(dǎo)下，表皮毛的發(fā)生受到抑制（Schellmann et al.， 2002）。毛狀體的發(fā)生發(fā)育機(jī)理研究表明其由多基因調(diào)控（張繼偉等，2018），變異型烏岡櫟葉背毛狀體作為山頂環(huán)境脅迫響應(yīng)的產(chǎn)物，可能經(jīng)歷了負(fù)調(diào)控因子突變、正調(diào)控因子轉(zhuǎn)座或增強(qiáng)等過(guò)程。堡狀峰頂環(huán)境因子與毛狀體發(fā)生的內(nèi)在聯(lián)系需進(jìn)一步研究方可明確。

氣孔長(zhǎng)寬比實(shí)際上反映了氣孔開(kāi)度（stomatal aperture），其變化依賴于受離子及糖控制膨壓的副衛(wèi)細(xì)胞、細(xì)胞骨架組織、跨膜運(yùn)輸及基因表達(dá)的聯(lián)合作用，大量證據(jù)顯示單一的刺激性信號(hào)通路不足以讓氣孔開(kāi)度發(fā)生變化（Hetherington amp; Woodward， 2003）。在本次研究的烏岡櫟種群間，與毛狀體同理，同一時(shí)刻原始型1的氣孔開(kāi)度似乎沒(méi)有因居處堡狀峰頂而與山腳或青云山的原始型2、3在這一性狀上表現(xiàn)出顯著差異。但是，在變異型中，氣孔長(zhǎng)寬比顯著增大，即氣孔形狀發(fā)生了變異，干旱、葉溫過(guò)高或過(guò)低經(jīng)信號(hào)傳遞最終在葉片上表現(xiàn)出氣孔開(kāi)度的減?。ń阎竦?，2015；段驊等，2019）。由此直觀地推測(cè)，在同等條件下變異型種群通過(guò)減少蒸騰以便應(yīng)對(duì)干旱或極端氣溫，并在大風(fēng)環(huán)境中關(guān)閉氣孔使機(jī)動(dòng)性得以增強(qiáng)。

綜上所述，堡狀峰地貌是變異產(chǎn)生的淵緣，該地貌生境，特別是貧瘠而薄的土壤，導(dǎo)致水熱條件差于周圍地區(qū)，STB、CTB及氣孔長(zhǎng)寬比的變異對(duì)適應(yīng)該生境有積極作用。這些性狀變異應(yīng)是穩(wěn)定的，它們并不受環(huán)境脅迫快速而直接地形成；經(jīng)過(guò)1年的觀察，在變異型或原始型內(nèi)，不同發(fā)育時(shí)期亦穩(wěn)定不變。

3.2 種群分化的可能途徑

烏山直接形成于新近紀(jì)后中國(guó)東南海岸地帶不同程度的隆起抬升與風(fēng)化侵蝕（Chen et al.，2009），更新世末次冰期導(dǎo)致高海拔植物下潛到附近“避難所”（陳冬梅等，2011），在間冰期或全新世（Holocene）很可能讓烏岡櫟重回山頂。嚙齒類是殼斗科種子傳播的主要貢獻(xiàn)者，但為抵達(dá)堡狀峰頂，最初的遷移可能依賴于鳥(niǎo)類（Oh， 2015）。原始型烏岡櫟抵達(dá)山頂后，若受生境脅迫發(fā)生適應(yīng)性演化，適應(yīng)性強(qiáng)的個(gè)體具有優(yōu)勢(shì)，進(jìn)而擴(kuò)大成種群。然而，山頂卻存在小范圍分布或零散于變異型周圍的原始型種群或個(gè)體。這表明原始型本身尚可忍耐干旱等極端生境。可以將其看作是基于自然選擇與競(jìng)爭(zhēng)排除原理種群殘存的一個(gè)剖面，同時(shí)，鳥(niǎo)類傳播這類偶然事件也可能對(duì)山頂群落演替產(chǎn)生重大影響，這均可導(dǎo)致當(dāng)前原始型與變異型共存，并且通過(guò)雜交讓過(guò)渡型產(chǎn)生。這表明過(guò)渡型之所以被觀察到，可能是物種形成過(guò)程殘留的隨機(jī)變異個(gè)體，抑或是雜交的產(chǎn)物。但是作為風(fēng)媒樹(shù)種，四周低地的原始型花粉流應(yīng)當(dāng)從未斷絕（Ashley， 2021）。紅旗巖堡狀峰的地理隔離作用可能相當(dāng)有限， 堡狀峰特殊生境顯然較地理隔離對(duì)變異的產(chǎn)生起到了更大的促進(jìn)作用。山頂?shù)倪^(guò)渡型及山腳未調(diào)查到變異型的現(xiàn)狀表明變異型可以與原始型融合而消失。但是就目前而言，堡狀峰的確產(chǎn)生了一定的生態(tài)學(xué)島嶼效應(yīng)。

3.3 祖先性狀和平行演化

較于氣孔長(zhǎng)寬比，毛狀體的變異在櫟屬內(nèi)有更多演化規(guī)律。首先需明確，烏山的變異型烏岡櫟與擁有類似葉背毛狀體的近緣種直線距離超過(guò)300 km，可排除雜交起源。

Deng等（2017）研究表明，STB是櫟屬的祖征，不同演化支或多或少存在，可在種間存在顯著的密度差異，烏山變異型與原始型烏岡櫟葉上下表皮STB即如此。然而，櫟屬分子系統(tǒng)結(jié)果顯示CTB在不同支系中多次起源或丟失（Yang et al.， 2017；Hipp et al.， 2020），冬青礫組不同種的情況不一，因此變異型發(fā)生CTB密度的變異是可能的。表皮毛發(fā)育受多基因直接與間接調(diào)控（普莉等，2003），選擇促進(jìn)基因互作與重排，古老祖先多態(tài)性分選也可使類似性狀反復(fù)產(chǎn)生。為解析它們的關(guān)系，需利用分子手段做各類型群體的研究。

另外，烏岡櫟的毛狀體變異在廣西的標(biāo)本（陳照宙50295），以及Global Biodiversity Information Facility （GBIF， https：//www.gbif.org/） 額外檢查的56份標(biāo)本中的2份來(lái)自琉球群島的標(biāo)本（Wright， 無(wú)號(hào); Type of Quercus wrightii = Q. phillyreoides）中發(fā)現(xiàn)，若這不是單株突變，那就表明它們之間可能發(fā)生了平行演化。

綜上得到以下結(jié)論：福建漳州烏山紅旗巖的烏岡櫟變異型的葉下表皮STB、CTB、氣孔長(zhǎng)寬比與原始型相比發(fā)生了自然變異，即STB、CTB顯著增多（Plt;0.01），氣孔長(zhǎng)寬比顯著增大（Plt;0.01）。烏山堡狀峰特殊生境促使變異產(chǎn)生，其具備一定生態(tài)與演化合理性。這為亞熱帶特殊花崗巖地貌促進(jìn)種群分化提供了有效證據(jù)。烏岡櫟廣泛分布于中國(guó)南方及日本，是毛狀體發(fā)生與演化、平行物種形成等研究的理想對(duì)象。中國(guó)是世界上花崗巖分布最廣的國(guó)家之一，又以廣東和福建兩省最多（占中國(guó)30%～40%）（崔之久等，2007）。構(gòu)造運(yùn)動(dòng)加速花崗巖出露地表，形成極其多樣的地貌景觀，中、高度隆起形成的堡狀峰、峰叢、石林、獨(dú)立巨峰等在大小空間尺度上廣泛存在。特殊花崗巖地貌上的植物區(qū)系、群落生態(tài)、極小種群、物種演化等有望被進(jìn)一步關(guān)注。

參考文獻(xiàn)：

ARENA ET， RUEDEN CT， HINER MC， et al.， 2017. Quantitating the cell： turning images into numbers with ImageJ ［J］. Wiley Interdiscip Rev Dev Biol， 6（2）：" e260.

ASHLEY MV， 2021. Answers blowing in the wind： A quarter century of genetic studies of pollination in oaks ［J］. Forests， 12（5）：" 575.

BARBARá"T， LEXER C， MARTINELLI G， et al.， 2008. Within-population spatial genetic structure in four naturally fragmented species of a neotropical inselberg radiation， Alcantarea imperialis， A-geniculata， A-glaziouana and A-regina （Bromeliaceae） ［J］. Heredity （Edinb）， 101（3）：" 285-296.

BARBARá T， MARTINELLI G， FAY MF， et al.， 2007. Population differentiation and species cohesion in two closely related plants adapted to neotropical high-altitude ‘inselbergs’， Alcantarea imperialis and Alcantarea geniculata （Bromeliaceae） ［J］. Mol Ecol， 16（10）：" 1981-1992.

BYRNE M， KRAUSS SL， MILLAR MA， et al.， 2019. Persistence and stochasticity are key determinants of genetic diversity in plants associated with banded iron formation inselbergs ［J］. Biol Rev Camb Philos Soc， 94（3）：" 753-772.

CHEN YX， CUI ZJ， YANG JQ， 2009. Influence of climate and tectonic movements on granite landforms in China ［J］. J Geogr Sci， 19（5）：" 587-599.

CHEN DM， KANG HZ， LIU CJ， 2011. An overview on the potential Quaternary glacial refugia of plants in China mainland ［J］. Bull Bot Res， 31（5）：" 623-632." ［陳冬梅， 康宏樟， 劉春江， 2011. 中國(guó)大陸第四紀(jì)冰期潛在植物避難所研究進(jìn)展 ［J］. 植物研究， 31（5）： 623-632.］

CONCEI??O AA， PIRANI JR， MEIRELLES ST， 2007. Floristics， structure and soil of insular vegetation in four quartzite-sandstone outcrops of “Chapada Diamantina”， Northeast Brazil ［J］. Rev Bras Bot， 30（4）：" 641-656.

CUI ZJ， YANG JQ， CHEN YX， 2007. The type and evolution of the granite landforms in China ［J］. Acta Geograph Sin， 62（7）：" 675-690." ［崔之久， 楊建強(qiáng)， 陳藝鑫， 2007. 中國(guó)花崗巖地貌的類型特征與演化 ［J］. 地理學(xué)報(bào)， 62（7）： 675-690.］

DENG M， JIANG XL， SONG YG， et al.， 2017. Leaf epidermal features of Quercus Group Ilex （Fagaceae） and their application to species identification ［J］. Rev Palaeobot Palynol， 237： 10-36.

DUAN H， TONG H， LIU YQ， et al.， 2019. Research advances in the effect of heat and drought on rice and its mechanism ［J］. Chin J Rice Sci， 33（3）：" 206-218." ［段驊， 佟卉， 劉燕清， 等， 2019. 高溫和干旱對(duì)水稻的影響及其機(jī)制的研究進(jìn)展 ［J］. 中國(guó)水稻科學(xué)， 33（3）： 206-218.］

Editorial Committee of Chinese Flora， Chinese Academy of Sciences， 1998. Flora of China：" Volume 22. Angiosperms， Dicotyledoneae， Fagaceae， Ulmaceae， Rhoipteleaceae ［M］. Beijing： Science Press： 255-257." ［中國(guó)科學(xué)院中國(guó)植物志編輯委員會(huì)， 1998. 中國(guó)植物志： 第二十二卷. 被子植物門 雙子葉植物綱 殼斗科 榆科 馬尾樹(shù)科 ［M］. 北京： 科學(xué)出版社： 255-257.］

FELICIANO DC， De GODOY SM， MARQUES DA SILVA JF， et al.， 2022. Landscape genetics reveal low diversity and adaptive divergence in Portulaca hatschbachii （Portulacaceae）：" an endangered species endemic to rocky outcrops of the Atlantic Forest ［J］. Bot J Linn Soc， 200（1）：" 116-141.

HE K， JIANG XL， 2014. Sky islands of southwest China. I： an overview of phylogeographic patterns ［J］. Chin Sci Bull， 59（7）：" 585-597.

HEGEBARTH D， BUSCHHAUS C， WU M， et al.， 2016. The composition of surface wax on trichomes of Arabidopsis thaliana differs from wax on other epidermal cells ［J］. Plant J， 88（5）： 762-774.

HENNERON L， SARTHOU C， DE MASSARY J， et al.， 2019. Habitat diversity associated to island size and environmental filtering control the species richness of rock-savanna plants in neotropical inselbergs ［J］. Ecography， 42（9）：" 1536-1547.

HETHERINGTON AM， WOODWARD FI， 2003. The role of stomata in sensing and driving environmental change ［J］. Nature， 424： 901-908.

HIPP AL， MANOS PS， HAHN M， et al.， 2020. Genomic landscape of the global oak phylogeny ［J］. New Phytol， 226（4）：" 1198-1212.

JIANG ZZ， ZHU HG， ZHANG Q， et al.， 2015. Progress of influence of low temperature on plant photosynthesis ［J］. Crops， （3）：" 23-28." ［姜籽竹， 朱恒光， 張倩， 等， 2015. 低溫脅迫下植物光合作用的研究進(jìn)展 ［J］. 作物雜志， （3）： 23-28.］

LIANG SJ， WEN FC， 2010. Morphologic characteristics of the miarolitic alkali feldspar granite in Taimushan Mountian Fujian Province ［J］. Geology Fujian， 29（4）：" 334-346." ［梁詩(shī)經(jīng)， 文斐成， 2010. 福建太姥山晶洞堿長(zhǎng)花崗巖地貌特征 ［J］. 福建地質(zhì)， 29（4）： 334-346.］

LIU XW， BARTHOLOMEW E， CAI YL， et al.， 2016. Trichome-related mutants provide a new perspective on multicellular trichome initiation and development in cucumber （Cucumis sativus L）" ［J］. Front Plant Sci， 7： 1187.

LUO D， YUE JP， SUN WG， et al.， 2016. Evolutionary history of the subnival flora of the Himalaya-Hengduan Mountains： first insights from comparative phylogeography of four perennial herbs ［J］. J Biogeogr， 43（1）： 31-43.

MA YL， WANG L， LIU YX， et al.， 2015. Uptates on stress tolerance of main accessory structures and their synergetic interaction in desert plants ［J］. Plant Physiol J， 51（11）： 1821-1836." ［馬亞麗， 王璐， 劉艷霞， 等， 2015. 荒漠植物幾種主要附屬結(jié)構(gòu)的抗逆功能及其協(xié)同調(diào)控的研究進(jìn)展 ［J］. 植物生理學(xué)報(bào)， 51（11）： 1821-1836.］

MALLET B， MARTOS F， BLAMBERT L， et al.， 2014. Evidence for isolation-by-habitat among populations of an epiphytic orchid species on a small oceanic island ［J］. PLoS ONE， 9（2）：" e87469.

MIL B， WARREN BH， HEEB P， et al.， 2010. The geographic scale of diversification on islands： genetic and morphological divergence at a very small spatial scale in the Mascarene grey white-eye （Aves： Zosterops borbonicus） ［J］. BMC Evol Biol， 10（1）：" 158.

MOLES AT， LAFFAN SW， KEIGHERY M， et al.， 2020. A hairy situation： Plant species in warm， sunny places are more likely to have pubescent leaves ［J］. J Biogeogr， 47（9）：" 1934-1944.

MOTA MR， PINHEIRO F， DOS SANTOS LEAL BS， et al.， 2020. From micro- to macro-evolution： insights from a Neotropical bromeliad with high population genetic structure adapted to rock outcrops ［J］. Heredity （Edinb）， 125（5）：" 353-370.

OH SH， 2015. Sea， wind， or bird： origin of Fagus multinervis （Fagaceae） inferred from chloroplast DNA sequences ［J］. Korean J Plant Taxon， 45（3）：" 213-220.

PENG DL， YANG N， SONG B， et al.， 2015. Woolly and overlapping leaves dampen temperature fluctuations in reproductive organ of an alpine Himalayan forb ［J］. J Plant Ecol， 8（2）：" 159-165.

POREMBSKI S， 2007. Tropical inselbergs： habitat types， adaptive strategies and diversity patterns ［J］. Rev Bras Bot， 30（4）：" 579-586.

POREMBSKI S， BARTHLOTT W， 2012. Inselbergs： biotic diversity of isolated rock outcrops in tropical and temperate regions ［M］. Berlin： Springer Science amp; Business Media： 103-176.

PU L， SUO JF， XUE YB， 2003. Molecular control of plant trichome development ［J］. Acta Genet Sin， 30（11）：" 1078-1084." ［普莉， 索金鳳， 薛勇彪， 2003. 植物表皮毛發(fā)育的分子遺傳控制 ［J］. 遺傳學(xué)報(bào)， 30（11）： 1078-1084.］

RUAS RDB， PAGGI GM， AGUIAR-MELO C， et al.， 2020. Strong genetic structure in Dyckia excelsa （Bromeliaceae）， an endangered species found on ironstone outcrops in Pantanal， Brazil ［J］. Bot J Linn Soc， 192（4）：" 691-705.

SáNCHEZ-ACEVEDOT V， TORRES-MIRANDA A， GUTIRREZ-ARELLANO C， et al.， 2022. Variability in leaf morphological traits of an endemic Mexican oak （Quercus mexicana Bonpl.） along an environmental gradientt ［J］. Bot Sci， 100（3）：" 579-599.

SCHELLMANN S， KIRIK V， WADA T， et al.， 2002.TRIPTYCHON and CAPRICE mediate lateral inhibition during trichome and root hair patterning in Arabidopsis ［J］. EMBO J， 21（19）：" 5036-5046.

TAVARES MM， FERRO M， LEAL BSS， et al.， 2022. Speciation with gene flow between two Neotropical sympatric species （Pitcairnia spp.： Bromeliaceae） ［J］. Ecol Evol， 12（5）：" e8834.

WU ZL， JIANG YX， 1997. Forest zoning in China ［M］//China Forest Editorial Committee. Forest in China. Beijing： China Forestry Publishing House： 437-512." ［吳中倫， 蔣有緒， 1997. 中國(guó)森林分區(qū) ［M］//中國(guó)森林編輯委員會(huì). 中國(guó)森林. 北京： 中國(guó)林業(yè)出版社： 437-512.］

XIE XF， YAN HF， WANG FY， et al.， 2012. Chloroplast DNA phylogeography of Primula ovalifolia in central and adjacent southwestern China： Past gradual expansion and geographical isolation ［J］. J Syst Evol， 50（4）：" 284-294.

XU ZZ， ZHOU GS， 2008. Responses of leaf stomatal density to water status and its relationship with photosynthesis in a grass ［J］. J Exp Bot， 59（12）：" 3317-3325.

YANG FS， QIN AL， LI YF， et al.， 2012. Great genetic differentiation among populations of Meconopsis integrifolia and its implication for plant speciation in the Qinghai-Tibetan Plateau ［J］. PLoS ONE， 7（5）：" e37196.

YANG J， VAZQUEZ L， CHEN XD， et al.， 2017. Development of chloroplast and nuclear DNA markers for Chinese oaks （Quercus Subgenus Quercus） and assessment of their utility as DNA barcodes ［J］. Front Plant Sci， 8： 816.

YANG YC， ZHOU T， QIAN ZQ，et al.， 2021. Phylogenetic relationships in Chinese oaks （Fagaceae， Quercus）：" Evidence from plastid genome using low-coverage whole genome sequencing ［J］. Genomics， 113（3）：" 1438-1447.

YANG ZY， YI TS， PAN YZ， et al.， 2012. Phylogeography of an alpine plant Ligularia vellerea （Asteraceae） in the Hengduan Mountains ［J］. J Syst Evol， 50（4）：" 316-324.

ZHANG D， HAO GQ， GUO XY， et al.， 2019. Genomic insight into “sky island” species diversification in a mountainous biodiversity hotspot ［J］. J Syst Evol， 57（6）：" 633-645.

ZHANG JW， ZHAO JJ， ZHOU Q， et al.， 2018. Progress in research of plant trichome ［J］. Chin Bull Bot， 53（5）：" 726-737." ［張繼偉， 趙杰才， 周琴， 等， 2018. 植物表皮毛研究進(jìn)展 ［J］. 植物學(xué)報(bào)， 53（5）： 726-737.］

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

基金項(xiàng)目： "國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目（2019YFD1000400）； 君子峰國(guó)家級(jí)自然保護(hù)區(qū)維管植物資源調(diào)查項(xiàng)目（閩財(cái)資環(huán)指 ［2019］15號(hào)）； 中醫(yī)藥公共衛(wèi)生服務(wù)補(bǔ)助專項(xiàng)（財(cái)社［2019］39號(hào)）。

第一作者： 林毅喆（1999—），碩士研究生，主要從事生物多樣性研究，（E-mail）linyizhe_0220@163.com。

*通信作者： "陳世品，博士，教授，主要從事植物分類學(xué)研究，（E-mail）fjcsp@126.com。