鄒萌，蔣瑞平，李佳倫，段淇洋，廖海，周嘉裕

西南交通大學(xué)生命科學(xué)與工程學(xué)院，四川 成都 610031

太白貝母Fritillaria taipaiensisP.Y.Li 為百合科貝母屬植物，是川貝母藥材的基原植物之一，以鱗莖入藥[1]。太白貝母含有生物堿成分，具有潤(rùn)肺止咳、化痰平喘等功效，臨床廣泛應(yīng)用于肺熱燥咳、肺陰 虛等。由于人為大量采挖和生態(tài)環(huán)境的破壞，太白貝母的野生資源受到嚴(yán)重破壞，已瀕臨枯竭[2]。鑒于太白貝母野生資源所面臨的嚴(yán)峻形勢(shì)，《陜西省中藥材保護(hù)和發(fā)展實(shí)施方案（2016－2020）》（陜政辦發(fā)〔2016〕9 號(hào)）已將太白貝母列入中藥材重點(diǎn)保護(hù)品種目錄。為緩解太白貝母需求與野生資源匱乏的矛盾，對(duì)太白貝母開(kāi)展引種與人工栽培是一種可行的替代方法。

有關(guān)太白貝母的栽培報(bào)道，始載于清光緒《大寧縣志》“貝母，銀廠坪所產(chǎn)為佳”?！按髮幙h”指今重慶市巫溪縣地區(qū)，說(shuō)明從清代起，當(dāng)?shù)匾殉蔀樘棕惸傅脑耘鄥^(qū)。其后，陜西、甘肅、重慶與四川等地也陸續(xù)有太白貝母的野生資源及栽培區(qū)的報(bào)道[3]，但盲目引種及擴(kuò)大栽培區(qū)會(huì)嚴(yán)重影響太白貝母生產(chǎn)的合理布局，甚至破壞太白貝母的道地性，導(dǎo)致藥材品質(zhì)下降。因此，開(kāi)展太白貝母生態(tài)環(huán)境適應(yīng)性研究具有重大現(xiàn)實(shí)意義。

本課題組通過(guò)2019 年6 月－2020 年9 月野外實(shí)地勘查、查詢(xún)相關(guān)文獻(xiàn)與網(wǎng)絡(luò)數(shù)據(jù)庫(kù)信息，整理出太白貝母在中國(guó)的地理分布區(qū)資料，利用最大熵（MaxEnt）模型和地理信息系統(tǒng)ArcGIS 軟件分析影響太白貝母地理分布的主要環(huán)境因子及其影響程度，并以此為依據(jù)分類(lèi)預(yù)測(cè)太白貝母不同程度適生區(qū)劃圖，分析其潛在地理分布格局，以期為太白貝母資源的可持續(xù)利用提供參考。

1 數(shù)據(jù)來(lái)源與方法

1.1 太白貝母分布數(shù)據(jù)

太白貝母的分布數(shù)據(jù)主要來(lái)源于3 個(gè)部分：①實(shí)驗(yàn)室人員野外實(shí)地考察，以及中國(guó)數(shù)字植物標(biāo)本館（http://www.cvh.org.cv/）標(biāo)本植物的分布地經(jīng)緯度信息；②通過(guò)查閱文獻(xiàn)中的分布區(qū)數(shù)據(jù)整理得到經(jīng)緯度信息，其中對(duì)僅有物種名及物種詳細(xì)分布地名稱(chēng)的，通過(guò)地球在線（https://www.earthol.com/）查找具體分布點(diǎn)并記錄；③通過(guò)全球生物多樣性信息平臺(tái)（https://www.gbif.org/）整理太白貝母在中國(guó)范圍內(nèi)分布點(diǎn)，并找到相應(yīng)的經(jīng)緯度信息。最終整理得到中國(guó)境內(nèi)115 份包含20 個(gè)地點(diǎn)的分布點(diǎn)數(shù)據(jù)，將各分布點(diǎn)數(shù)據(jù)統(tǒng)一轉(zhuǎn)換為十進(jìn)制并以.csv 格式保存。詳見(jiàn)表1。

表1 太白貝母分布數(shù)據(jù)

1.2 環(huán)境變量

從世界氣象數(shù)據(jù)（https://www.worldclim.org/）下載32 個(gè)環(huán)境變量因子，涵蓋19 個(gè)基礎(chǔ)生物氣候變量（bio1～19），以及1、5、6、7 月的最低溫度（tmin1、tmin5、tmin6、tmin7）、最高溫度（tmax1、tmax5、tmax6、tmax7）、平均溫度（tavg1、tavg5、tavg6、tavg7）和海拔（alt）。采用ChinaMap 中4p 格式的1∶400 萬(wàn)中華人民共和國(guó)地圖和中華人民共和國(guó)行政區(qū)劃圖（https://www.csdn.net/）作為底圖得到太白貝母在中國(guó)境內(nèi)的分布點(diǎn)。詳見(jiàn)圖1。

圖1 太白貝母分布點(diǎn)

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

1.3.1 環(huán)境因子選取

由于環(huán)境因子間有一定的相關(guān)性，所以本研究對(duì)32 個(gè)環(huán)境因子進(jìn)行相關(guān)性檢驗(yàn)后用于MaxEnt 模型。利用ArcGIS12.0 軟件提取各采樣點(diǎn)的點(diǎn)插值，以SPSS20.0 軟件Spearman 相關(guān)系數(shù)查驗(yàn)相關(guān)性，選取相關(guān)性＜0.8 的環(huán)境因子，若2 個(gè)因子相關(guān)性≥0.8，只選擇任意1 個(gè)因子運(yùn)行模型。

1.3.2 最大熵模型預(yù)測(cè)

運(yùn)行MaxEnt3.4.1 軟件，在樣本欄（samples）添加上述.csv 文件，在環(huán)境圖層欄（environmental layers）添加所選取的低相關(guān)性環(huán)境因子，采用刀切法（jackknife）預(yù)測(cè)各環(huán)境因子的貢獻(xiàn)率，制作響應(yīng)曲線（response curve）并生成初步適生區(qū)預(yù)測(cè)圖像。預(yù)測(cè)模型的準(zhǔn)確性是利用受試者工作特征（ROC）曲線下方面積（AUC）值判定，其值越接近1，說(shuō)明與隨機(jī)分布相距越遠(yuǎn)，且環(huán)境因子與預(yù)測(cè)的物種地理分布之間的相關(guān)性越大，即MaxEnt 模型預(yù)測(cè)的結(jié)果越準(zhǔn)確。

本研究采用Roberto、Natural breaks 等研究方法，選取25%分布點(diǎn)作為測(cè)試集與75%分布點(diǎn)作為訓(xùn)練集，模型運(yùn)算次數(shù)設(shè)定為500 次，重復(fù)運(yùn)算設(shè)定為10 次，其余參數(shù)均為默認(rèn)模型的自動(dòng)設(shè)置。

1.4 適應(yīng)性分布區(qū)劃分

將MaxEnt 模型運(yùn)行得到的結(jié)果導(dǎo)入ArcGIS12.0軟件，利用重分類(lèi)進(jìn)行疊加分析和地圖制作，采用人工分級(jí)方法劃分出太白貝母適應(yīng)性分布等級(jí)，得到主要生態(tài)因子影響下的適應(yīng)性區(qū)劃圖。根據(jù)MaxEnt 模型運(yùn)行結(jié)果中各環(huán)境因子的響應(yīng)曲線得出主要因子適宜值范圍及太白貝母適生區(qū)的生態(tài)特征。

2 結(jié)果與分析

2.1 預(yù)測(cè)的有效性

近年來(lái)，ROC 曲線分析法在預(yù)測(cè)物種潛在適生區(qū)模型評(píng)價(jià)中應(yīng)用廣泛[4]。該方法是對(duì)曲線下方面積進(jìn)行計(jì)算，即得到AUC 值，由此判斷模型模擬的精準(zhǔn)度。一般認(rèn)為，AUC＜0.6 時(shí)判斷價(jià)值較低，AUC＞0.8 時(shí)判斷價(jià)值較高，AUC＝1 時(shí)為理想分布[5]。本研究在對(duì)各環(huán)境因子重復(fù)運(yùn)行10 次后得到訓(xùn)練數(shù)據(jù)集AUC＝0.958，說(shuō)明MaxEnt 模型對(duì)太白貝母潛在分布區(qū)的預(yù)測(cè)結(jié)果較好，可信度非常高。見(jiàn)圖2。

圖2 太白貝母分布MaxEnt 模型ROC 曲線及AUC 值

2.2 環(huán)境因子貢獻(xiàn)率

jackknife 是在預(yù)測(cè)物種潛在分布的適生區(qū)和分析環(huán)境因子對(duì)所預(yù)測(cè)結(jié)果的影響程度大小，以及通過(guò)運(yùn)算確定對(duì)物種分布貢獻(xiàn)率最大的環(huán)境因子過(guò)程中最常用的方法。本研究根據(jù)jackknife 分析得到的預(yù)測(cè)結(jié)果顯示，各環(huán)境因子在對(duì)太白貝母適應(yīng)度的影響所占的權(quán)重不同。見(jiàn)圖3。

圖3 環(huán)境因子在MaxEnt 模型中預(yù)測(cè)太白貝母分布刀切圖

通過(guò)模型篩選得到對(duì)太白貝母潛在分布貢獻(xiàn)率＞0 的環(huán)境因子有6 個(gè)，即晝夜溫差均值（41.21%）、海拔（26.42%）、年降水量（18.32%）、溫度季節(jié)性變化的標(biāo)準(zhǔn)差（13.43%）、最濕月降水量（0.63%）、6月最低氣溫（0.12%）。其中晝夜溫差均值和海拔的貢獻(xiàn)率最高，溫度、海拔與降水是影響太白貝母分布區(qū)的重要環(huán)境因素，見(jiàn)表2。

表2 對(duì)太白貝母潛在分布貢獻(xiàn)率＞0 的6 個(gè)環(huán)境因子

2.3 主要環(huán)境因子

通過(guò)運(yùn)行MaxEnt 模型可以得到所選取各環(huán)境因子的單變量響應(yīng)曲線，從而分析其最適區(qū)間值。為明確主要環(huán)境因子的參數(shù)值與太白貝母存在概率之間的關(guān)系，本研究在MaxEnt 模型運(yùn)行結(jié)果基礎(chǔ)上繪制僅包含單一環(huán)境因子的響應(yīng)曲線，其中包含bio2、alt、bio12 與bio4。結(jié)果表明，太白貝母的潛在分布概率隨環(huán)境的改變波動(dòng)值較大。選擇存在概率＞0.2 的區(qū)間作為其適生區(qū)范圍。

晝夜溫差均值和海拔影響最大，其中當(dāng)晝夜溫差均值約為8.1 ℃時(shí)太白貝母適生概率最大，達(dá)到66.28%，晝夜溫差2.7～11.9 ℃為其適生區(qū)范圍。當(dāng)海拔為241 m 時(shí)達(dá)到太白貝母適生范圍；241 m＜海拔＜1 318 m時(shí)呈增長(zhǎng)趨勢(shì)；1 318 m時(shí)到達(dá)最適概率，約為65.7%；1 318 m＜海拔＜3 828 m 時(shí)適生概率逐漸降低。因此，在晝夜溫差約為8.1 ℃及海拔約為1 318 m 時(shí)，最適宜太白貝母生長(zhǎng)。詳見(jiàn)圖4、圖5。

圖4 太白貝母晝夜溫差均值響應(yīng)曲線

圖5 太白貝母海拔響應(yīng)曲線

2.4 太白貝母的適應(yīng)性劃分

將MaxEnt 模型的預(yù)測(cè)結(jié)果導(dǎo)入ArcGIS 軟件中運(yùn)行，對(duì)總體適生區(qū)進(jìn)行重分類(lèi)，確定把適生等級(jí)分為4 類(lèi)：0＜適應(yīng)指數(shù)≤0.06 為非適生區(qū)；0.06＜適應(yīng)指數(shù)≤0.23 為低適生區(qū)，主要分布于四川、甘肅南部、云南南部、浙江南部，在青海東南部、新疆北部和西藏西部有零星分布；0.23＜適應(yīng)指數(shù)≤0.51 為適生區(qū)，主要分布于甘肅南部、貴州西北部、四川中部及東北部、重慶南部、陜西南部及湖北西部；0.51＜適生指數(shù)＜0.86 為高適生區(qū)，主要分布于甘肅南部、貴州西北部、四川中部及東北部、重慶南部、陜西南部及湖北西部等地區(qū)。見(jiàn)圖6。潛在分布區(qū)在各?。▍^(qū)）的分布面積統(tǒng)計(jì)情況見(jiàn)表3。

圖6 太白貝母潛在分布核心區(qū)域

表3 太白貝母潛在分布區(qū)域面積

3 討論

利用生態(tài)位模型的運(yùn)行對(duì)物種的潛在分布區(qū)進(jìn)行預(yù)測(cè)，進(jìn)而對(duì)物種建立保護(hù)區(qū)及規(guī)劃人工栽培地域是目前廣泛應(yīng)用的方法。MaxEnt 模型運(yùn)行簡(jiǎn)單，預(yù)測(cè)結(jié)果的準(zhǔn)確性高且可信度較強(qiáng)，已被廣泛用于物種的潛在分布區(qū)預(yù)測(cè)及適應(yīng)性評(píng)價(jià)[6]。

本研究結(jié)合MaxEnt 模型和ArcGIS，采用的分布點(diǎn)數(shù)據(jù)涵蓋了太白貝母的現(xiàn)有分布區(qū)，結(jié)果發(fā)現(xiàn)太白貝母主要集中在陜西、甘肅、重慶及四川等地，與陳士林等[7]關(guān)于太白貝母最大生境區(qū)域的研究結(jié)果一致。

基于收集分布點(diǎn)的32 個(gè)環(huán)境因子數(shù)據(jù)，本研究對(duì)太白貝母的潛在分布區(qū)進(jìn)行預(yù)測(cè)及適應(yīng)性等級(jí)劃分。結(jié)果MaxEnt 模型ROC 曲線訓(xùn)練集AUC＝0.958，表示該模型預(yù)測(cè)準(zhǔn)確度達(dá)到較高水平，可用于太白貝母潛在分布區(qū)的預(yù)測(cè)。較之現(xiàn)有分布區(qū)[8]，潛在分布區(qū)增加了貴州、湖南及湖北的部分地區(qū)，其中貴州尤其值得關(guān)注，不僅與重慶接壤，且其總面積的76%可作為太白貝母的潛在分布區(qū)，表明貴州極有可能成為未來(lái)太白貝母產(chǎn)業(yè)的重點(diǎn)發(fā)展區(qū)域，需提前開(kāi)展合理布局與規(guī)劃。潛在分布區(qū)中，高適生區(qū)是太白貝母生長(zhǎng)與繁殖的最適區(qū)域，主要分布于甘肅南部、貴州西北部、四川中部和東北部、重慶南部、陜西南部及湖北西部等地區(qū)。這些區(qū)域集中了太白貝母的一些傳統(tǒng)產(chǎn)區(qū)，如陜西太白縣晝夜溫差大，氣候較寒冷，一直是太白貝母種植最適宜的地區(qū)[9]，印證了本研究預(yù)測(cè)結(jié)果的準(zhǔn)確性。太白貝母的潛在核心分布區(qū)與川貝母Fritillaria cirrhosaD.Don、暗紫貝母Fritillaria unibracteataP.K.Hsiao &K.C.Hsia 和甘肅貝母Fritillaria przewalskiiMaxim.ex Batal 等品種的分布區(qū)有一定范圍的重疊，其中在甘肅南部與甘肅貝母Fritillaria przewalskiiMaxim.ex Batal 有重疊，在四川省與暗紫貝母Fritillaria unibracteataP.K.Hsiao &K.C.Hsia 和川貝母Fritillaria cirrhosaD.Don 重疊，在湖北西部和重慶市與湖北貝母Fritillaria hupehensisHsiao et K.C.Hsia 有重疊。這些重疊區(qū)的存在，一方面體現(xiàn)出不同品種貝母所需的生態(tài)環(huán)境因子高度相似，另一方面也表明這些區(qū)域擁有較好的生態(tài)優(yōu)越性，盡管出現(xiàn)了不同貝母的空間生態(tài)位重疊，但彼此間卻不產(chǎn)生生存競(jìng)爭(zhēng)。

在影響太白貝母分布的主導(dǎo)環(huán)境因子中，氣候因子對(duì)太白貝母的影響大于土壤與地形因子。晝夜溫差均值是影響太白貝母分布的最重要環(huán)境因子，當(dāng)晝夜溫差均值為8.1 ℃時(shí)，最有利于太白貝母生長(zhǎng)過(guò)程中營(yíng)養(yǎng)物質(zhì)的累積，并促進(jìn)太白貝母鱗莖的生長(zhǎng)。有報(bào)道，等溫性是影響川貝母Fritillaria cirrhosaD.Don分布的最重要環(huán)境因子[10]。太白貝母與川貝母Fritillaria cirrhosaD.Don 具有不同的最重要環(huán)境因子，推測(cè)其原因可能在于太白貝母屬于川貝母類(lèi)群中較為獨(dú)立的品種。太白貝母的獨(dú)立性也體現(xiàn)在DNA上，其ITS2 等條形碼序列與其他品種川貝母有較大差別[11]，在進(jìn)化樹(shù)上處于不同分支，表明太白貝母較早從川貝母類(lèi)群中分化出來(lái)，導(dǎo)致其DNA 與所處環(huán)境與其他川貝母品種產(chǎn)生較大差別。由此，我們推測(cè)植物的生長(zhǎng)環(huán)境適應(yīng)性可能作為植物系統(tǒng)進(jìn)化研究的一個(gè)參數(shù)，但有待于進(jìn)一步驗(yàn)證。

本研究顯示，海拔是影響太白貝母分布的第2 個(gè)重要因素，最適海拔為1 318 m。研究顯示，不同海拔高度影響藥材品質(zhì)，海拔梯度造成的水分、溫度及土壤肥力等環(huán)境差異性會(huì)對(duì)植物生長(zhǎng)發(fā)育產(chǎn)生極大的影響，進(jìn)而導(dǎo)致植物生理和生態(tài)適應(yīng)性的變化[12]，如生長(zhǎng)在海拔1 500 m 的重慶產(chǎn)桔梗藥材質(zhì)量最優(yōu)[13]。由于暗紫貝母等其余川貝母品種的最低海拔在3 000 m 左右[14]，因此同一地區(qū)不同海拔可種植不同類(lèi)型的川貝母品種，如低處種植太白貝母，高處種植川貝母，如此則能有效利用生境資源，實(shí)現(xiàn)空間利用最大化。

目前川貝母的野生資源不容樂(lè)觀，而栽培技術(shù)已較為成熟，所以對(duì)太白貝母的適生區(qū)進(jìn)行預(yù)測(cè)并以此指導(dǎo)進(jìn)行合理種植與栽培十分必要。相較其他品種，太白貝母有望成為低海拔地區(qū)川貝母栽培的優(yōu)良候選品種。在合理規(guī)劃前提下建立太白貝母中藥材生產(chǎn)質(zhì)量管理規(guī)范（GAP）生產(chǎn)基地，對(duì)太白貝母資源的可持續(xù)利用有重要意義。