• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于譜系信號的寶天曼省級保護植物物種維持機制研究

    2020-10-14 07:15:25閆滿玉白兵勇
    林業(yè)科技 2020年2期

    閆滿玉 姚 松 白兵勇

    摘要: ?親緣關(guān)系相近的物種在形態(tài)特征等方面具有更多的相似性,譜系信號分析對生態(tài)位保守、物種共存機制方面的研究具有重要意義。本文選用可計算雙變量性狀的D值法和EM-mantel方法,計算寶天曼保護區(qū)內(nèi)省級重點保護被子植物不同生活型的譜系信號。結(jié)果顯示:喬木型具有強烈的譜系信號,灌木和草本型兩種方法的計算結(jié)果具有一定的差異;分析表明:喬木型的主要致危機制是生境過濾,而灌木和草本型的主要致危機制為人為干擾;最后,對保護區(qū)植物保護工作提出建議。

    關(guān)鍵詞: ?保護植物; ?譜系信號; ?物種共存; ?寶天曼自然保護區(qū)

    中圖分類號: ? S 759. 9 ? ? ? ? ? ? ? 文獻標(biāo)識碼: ? A ? ? ? ? ? ? ? ?文章編號:1001 - 9499(2020)02 - 0026 - 06

    物種在生物學(xué)和其他統(tǒng)計單元(statistic units)上都不是獨立的,而是具有相似的進化歷史[ 1 ]。親緣關(guān)系近的物種,在形態(tài)特征、生理性狀、行為特征等方面具有更多的相似性[ 2 ]。譜系信息主要反映生物進化過程中的分化時間差異,同樣的分化時間,相近的物種具有較近的親緣關(guān)系,分化時間相異的物種親緣關(guān)系較遠[ 3 ]。研究表明,物種滅絕不是隨機的,而是在分類學(xué)和譜系學(xué)上呈現(xiàn)出一定的聚集效應(yīng),親緣關(guān)系近的物種滅絕風(fēng)險相似[ 4 - 5 ]。

    研究物種間譜系和性狀間的關(guān)系(也就是譜系信號 Phylogenetic Signal)對認(rèn)識譜系生態(tài)位保守具有重要意義[ 6 - 8 ],有助于提高群落物種共存機制方面的認(rèn)識[ 9 - 11 ]。近年來有關(guān)譜系信號檢驗的研究層出不窮[ 12 - 15 ],截至到目前,主要有以下5種方法:1)Pagel,s λ法,利用布朗進化模型檢驗譜系信號的早期方法[ 16 ];2)K值法,是以布朗運動為進化模型的系統(tǒng)發(fā)育信號的檢驗方法,目前該方法最為常用,但僅能分析連續(xù)變量性狀[ 17 ];3) PVR法,該方法利用特征向量進行譜系信號檢驗,能較好的解釋布郎運動模型在進化速率方面的偏差[ 18 ];4) D值法,該方法最先應(yīng)用在動物雙變量性狀譜系信號檢測,但一次僅能處理一個性狀[ 19 ];5)EM-mantel檢驗法,結(jié)合進化模型的mantel檢驗法,可以處理多種性狀類型[ 1 ]。前3種方法在以往的研究中使用較多,但是僅僅適用于連續(xù)變量性狀,并且一次僅能計算一個性狀。鑒于上述因素限制,F(xiàn)ritz和Purvis于2010年提出了一種可以計算雙變量和連續(xù)變量性狀的譜系信號方法,該方法廣泛應(yīng)用在生態(tài)學(xué)和動物學(xué)領(lǐng)域[ 20 - 21 ],但是D值法同樣存在一次僅能計算一個性狀的劣勢。最近Vanderlei和 Leandro提出了一種結(jié)合進化模型的譜系信號矩陣相關(guān)檢驗方法(EM-mantel test),該方法不僅能夠處理連續(xù)型變量性狀,還能夠計算雙變量型性狀,同時還具備一次分析多個性狀的優(yōu)勢[ 1 ]。

    寶天曼自然保護區(qū)地處暖溫帶和北亞熱帶交匯區(qū),野生動植物資源豐富,具有重要的經(jīng)濟、社會和科研價值[ 22 - 23 ]。以往該區(qū)的珍稀瀕危植物的研究主要關(guān)注物種科屬分布、區(qū)系組分、生態(tài)分布、生長特點、生活型[ 23 - 25 ];與相鄰保護區(qū)的對比研究分析[ 26 ]、保育恢復(fù)措施[ 27 ]等,尚未有從進化的角度進行深層次的研究探討。因此,本研究結(jié)合2005年河南省人民政府頒布的《河南省重點保護植物名錄》[ 28 ]和本區(qū)已有研究結(jié)果[ 21 , 29 - 30 ],比對篩選得到寶天曼保護區(qū)省級重點保護被子植物名錄,利用Phylomatic軟件[ 31 ]建立帶有支長的譜系樹;以植物生活型為檢驗性狀(雙變量性狀),選用能夠計算雙變量性狀譜系信號的D值法和EM-mantel法,檢驗不同生活型的省級重點保護植物是否存在譜系信號及強度;探討珍稀瀕危植物可能存在的維持機制,為保護區(qū)今后開展重點保護植物保育工作提供科學(xué)依據(jù)。

    1 方 法

    1. 1 研究區(qū)概況

    寶天曼國家級自然保護區(qū)(33.34°~33.60°N,111.78°~112.06° E)位于河南省西南部內(nèi)鄉(xiāng)縣境內(nèi),山勢呈東南——西北走向,總面積5 412 hm2,最高峰海拔1 830 m。保護區(qū)地處南暖溫帶向北亞熱帶過渡區(qū)域,是我國同緯度中森林生態(tài)系統(tǒng)保存較為完整、物種最豐富的地區(qū)之一。該區(qū)屬大陸性季風(fēng)氣候,四季分明,年均溫度15.2 ℃,年降雨量885.6 mm,年蒸發(fā)量991.6 mm,相對濕度68%,極端最高溫41.2 ℃,極端最低溫-14.8℃。植被區(qū)劃屬于暖溫帶落葉闊葉林地帶,垂直分布比較明顯[ 22 , 32 ]。寶天曼古老的地質(zhì)和復(fù)雜的地形,為植物的生存繁衍提供了得天獨厚的環(huán)境條件,共有種子植物150科783屬2 147種(含種以下等級)[ 23 ]。

    1. 2 譜系樹構(gòu)建

    首先根據(jù)2005年河南省人民政府頒布的《河南省重點保護植物名錄》[ 28 ]、省級重點保護植物名校正[ 23 ]和本保護區(qū)相關(guān)研究 [ 29 - 30 ],比對篩選得到寶天曼保護區(qū)省級重點保護被子植物64種;其次利用被子植物分類系統(tǒng)III(APG III)[ 33 ]校正物種的科屬信息;最后,將物種學(xué)名以及相關(guān)科屬信息輸入到網(wǎng)頁版Phylomatic V3軟件[ 31 ],利用其數(shù)據(jù)庫zanne2014構(gòu)建帶有分支長度的譜系樹(圖1)。

    1. 3 譜系信號檢測

    本研究生活型性狀數(shù)據(jù)屬性為雙變量(用1和0表示,1表示具有某種性狀,0表示無此種性狀),綜合考量上述5種譜系信號檢測方法,選用能夠計算雙變量性狀譜系信號的D值法和EM-mantel法。

    1. 3. 1 D 值法

    D值法是通過比較譜系樹所有姊妹分支差異之和的大小來反映是否存在譜系信號及其強度,性狀聚集的姊妹分支差異之和較小,性狀隨機的姊妹分支的差異之和較大。公式為:

    其中dobs為姊妹分支長度的真實差異,db是模擬布朗運動[ 19 , 34 ]情況下的姊妹分支長度的差異,dr表示模擬隨機情況下的姊妹分支長度的差異。如果D=1,說明性狀隨機分布,沒有譜系信號;如果D=0,說明性狀聚集,具有譜系信號;如果D<0,說明性狀高度聚集,譜系信號強烈。

    1. 3. 2 EM-mantel 方法

    傳統(tǒng)上多性狀譜系信號檢測廣泛采用矩陣相關(guān)性檢驗(mantel test)[ 7 , 35 - 36 ],但矩陣相關(guān)性檢驗的零假設(shè)是兩矩陣間差異為零,而物種卻有相似的進化歷程[ 1 ],因此該方法存在解釋率較低和Ⅰ類錯誤(TypeⅠerror)率較高的弊端[ 37 ]。結(jié)合進化模型的矩陣相關(guān)性檢驗(EM-mantel test)之所以能夠解決傳統(tǒng)方法存在的諸多問題,是由于EM-mantel方法的零假設(shè)是在某一特定的進化模型下,譜系距離矩陣和性狀相似性矩陣兩者間不存在差異。布朗運動模型在解釋性狀進化、內(nèi)稟性狀值等方面不存在其他假設(shè)機制[ 1 ],同時D方法選用布朗進化模型,為便于兩種方法間的比較分析,EM- ?mantel法選用布朗進化模型。

    EM-mamtel法首先根據(jù)布朗進化模型生成模擬矩陣,其次利用傳統(tǒng)矩陣相關(guān)性檢驗方法分析布朗進化模型矩陣和真實矩陣間是否存在差異,最后根據(jù)p值確定是否有譜系信號及其強度。計算步驟簡述如下:1)以譜系樹數(shù)據(jù)為依據(jù)生成布朗進化模型下的性狀值(用1和0表示,1表示具有某種性狀,0表示無此種性狀);2)根據(jù)譜系樹支長數(shù)據(jù)模擬性狀分布;3)利用傳統(tǒng)矩陣相關(guān)性檢驗方法計算真實矩陣和模擬矩陣兩者的差異;4)重復(fù)上述計算多次,最后得到雙尾型p值。該方法的詳細介紹見Vanderlei和 Leandro的文章[ 1 ]。

    本研究所有計算均由R 3.4.0軟件(https://cran.r-pro-ject.org/)實現(xiàn),其中物種APG III分類系統(tǒng)科屬校正使用plantlist軟件包,D值法使用caper軟件包,EM-mantel法使用EM.mantel函數(shù) [ 1 ],兩種方法均模擬999次。

    2 結(jié) 果

    2. 1 寶天曼保護區(qū)省級重點保護被子植物名錄

    通過比較篩選該區(qū)省級重點保護被子植物共64種(表1)。其中,喬木34種,占53.13%,如杈葉槭(Acer ?robustum )、米心水青岡(Fagus ?engleriana)、紫莖(Stewartia sinensis)等;灌木9種,占14.06%,如蝟實(Kolkwitzia ?amabilis)、河南杜鵑(Rhododendron ?henanense)、玉鈴花(Styrax ?obassis)等;草本20種,占33.33%,如多花蘭(Cymbidium ?floribundum)、絞股藍(Gynostemma ?pentaphyllum)、延齡草(Trillium ?tschonoskii)等;藤本僅黑蕊獼猴桃(Actinidia ?melanandra)1種??茖匐`屬關(guān)系,蘭科10種,占15.63%;樟科8種占12.50%;毛茛科和槭樹科各4種,各占6.25%;其他22科共占65.62%。區(qū)系組成分析,北溫帶分布19種,占29.69%;東亞——北美洲際間斷分布和中國特有分布均為7種,共占21.08%;熱帶亞洲類型分布共18種,共占27.54%;其他4種類型共13種,共占21.69%。

    2. 2 省級重點保護被子植物不同生活型譜系信號分析

    區(qū)內(nèi)省級重點保護藤本植物僅黑蕊獼猴桃1種,無法計算此生活型的譜系信號,故本研究僅選取喬木、灌木和草本3種生活型(表2)。由此可知:兩種方法的計算結(jié)果存在一定差異: D值法3種生活型(喬木、灌木和草本)均譜系信號顯著,其D值分別為-1.07、-0.22和-1.34;而EM-mantel法僅喬木類型譜系信號顯著(p=0.001),其他兩種生活型譜系信號不顯著。

    3 分析與討論

    3. 1 兩種方法計算結(jié)果差異分析

    本研究選取能夠計算雙變量性狀譜系信號的兩種方法,計算結(jié)果顯示灌木和草本兩種生活型存在差異:D值法計算結(jié)果顯示兩種生活型譜系信號均顯著,而EM-mantel法結(jié)果顯示譜系信號不顯著。Hardy 和Pavoine研究表明,Mantel檢驗結(jié)果與物種的功能性狀矩陣和譜系距離矩陣相關(guān)[ 38 ],本研究中草本共20種隸屬9科17屬,其中50%為蘭科植物,而喬木生活型共34種隸屬15科26屬,樟科占比最高(僅占23.53%),因此兩者性狀矩陣的差異導(dǎo)致EM-mantle法計算結(jié)果不顯著。另外,本研究為方便比較分析,EM-mantel法選用布朗運進化模型,然而該方法尚有其他復(fù)雜模型,如Ornstein Uhlenbeck 模型[ 39 ]、Burst模型[ 40 ]等,模型選擇差異可能對計算結(jié)果產(chǎn)生一定的影響。此外,物種庫太大或太小均可影響生態(tài)過程的解釋能力,譜系分析中物種占區(qū)域物種庫30%~60%均在可接受范圍[ 43 ],本研究中灌木僅有9種(占物種庫14.06%),也可能對譜系分析結(jié)果造成一定的影響。

    3. 2 省級重點保護植物致危機制探討

    物種固然有自然消亡的現(xiàn)象,自然選擇的結(jié)果使地球上眾多物種淘汰,但是這種消亡過程是極其緩慢的。本研究區(qū)珍稀瀕危植物數(shù)量稀少,分布區(qū)域狹窄,一些種類正處于瀕危狀態(tài),究其原因有歷史、人為和植物自身等多方面原因[ 23 ]。古植物學(xué)研究表明,我國有過多次植物群的南北進退,特別是在第四紀(jì)冰川襲擊,而本區(qū)正處于遷移的中間地帶[ 22 ],一些珍稀植物在局部優(yōu)越的生境下得以生存和繁衍。本區(qū)重點保護植物主要分布于封閉山谷,避風(fēng)坡林地、溝壑雜木林及其林下、灌叢和林緣山坡等,除水曲柳、領(lǐng)春木、青檀和金錢槭等種類有散生小片群落分布外,絕大部分為零星植株分布,如銀鵲樹16株、水青樹和天竺桂等不超過10株[ 22 ]。因此生境過濾是喬木生活型省級重點保護植物的主要維持機制。

    本研究中,喬木生活型在兩種方法下均表現(xiàn)出強烈的譜系信號,這也證明了如果呈現(xiàn)出強烈的譜系聚集格局,那么性狀在決定物種滅絕風(fēng)險之中就起著直接或間接的作用,也就是說生境過濾是其主要維持機制這一研究結(jié)論的正確性。絕大部分灌木和草本為傳統(tǒng)藥用或觀賞植物,比如紫斑牡丹、延齡草、大葉三七、天麻等,另外草本類型中僅蘭科占15.63%,因此人為干擾因素占主要作用。同時EM-mantle法的計算結(jié)果顯示,這兩種生活型譜系信號不顯著,因此人為因素干擾是灌木和草本生活型重點保護植物的主要致危因子。研究結(jié)果也在一定層面上證實了:如果譜系呈現(xiàn)隨機的分布格局,那么性狀不是導(dǎo)致物種滅絕的限制因子,也就是說生態(tài)位理論為其主要維持機制的正確性。

    綜上,本區(qū)喬木生活型重點保護植物的主要致危因子為生境過濾,而灌木和草本生活型的主要致危因子是人為干擾。

    3. 3 寶天曼保護區(qū)重點保護植物保護措施

    自然生態(tài)環(huán)境是植物生存繁衍的基礎(chǔ),特定的生態(tài)環(huán)境和植物之間構(gòu)成了相互依存的群落關(guān)系,一旦立地生境遭到破壞,重點保護植物將面臨滅頂之災(zāi)。為了更好地保護和發(fā)展本區(qū)重點保護植物資源,根據(jù)該研究分析結(jié)果和以往研究[ 26 ],提出針對不同生活型采取分類保護的措施,即對森林生態(tài)系統(tǒng)的優(yōu)勢種——喬木型重點保護植物,應(yīng)該重點保護其立地生境。對灌木和草本型植物來說,人為因素干擾是其重要致危因子,其保護措施應(yīng)重點著眼于減少人為干擾。

    參考文獻

    [1] Debastiani V J, da Silva Duarte L. Evolutionary models and phylogenetic signal assessment via Mantel test[J]. Evolutionary Biology, 2017, 44(1): 135 - 143.

    [2] Harvey P H, Pagel M D. The comparative method in evolutionary biology[M]. Oxford: Oxford university press, 1991.

    [3] Webb C O, Losos J B, Agrawal A A. Integrating phylogenies into community ecology[J]. Ecology, 2006, 87(sp7).

    [4] Purvis A. Phylogenetic approaches to the study of extinction [J]. Annual Review of Ecology, Evolution, and Systematics, 2008, 39.

    [5] Roy K, Hunt G, Jablonski D. Phylogenetic conservatism of extinctions in marine bivalves[J]. Science, 2009, 325(5941): 733 - 737.

    [6] Wiens J J, Graham C H. Niche conservatism: integrating evolution, ecology, and conservation biology[J]. Annu. Rev. Ecol. Evol. Syst., 2005, 36: 519 - 539.

    [7] Losos J B. Phylogenetic niche conservatism, phylogenetic signal and the relationship between phylogenetic relatedness and ecological similarity among species[J]. Ecology letters, 2008, 11(10): 995 - 1 003.

    [8] Cooper N, Jetz W, Freckleton R P. Phylogenetic comparativeapproaches for studying nicheconservatism[J]. Journal of evolutionary biology, 2010, 23(12): 2 529 - 2 539.

    [9] Webb C O, Ackerly D D, McPeek M A, et al. Phylogenies and community ecology[J]. Annual review of ecology and systematics, 2002, 33(1): 475 - 505.

    [10] Mouquet N, Devictor V, Meynard C N, et al. Ecophylogenetics:advances and perspectives[J]. Biological reviews, 2012, 87(4): 769 - 785.

    [11] Gerhold P, Cahill J F, Winter M, et al. Phylogenetic patterns are not proxies of community assembly mechanisms (they are far better)[J]. Functional Ecology, 2015, 29(5): 600 - 614.

    [12] Gilbert G S, Webb C O. Phylogenetic signal in plant pathogen- host range[J]. Proceedings of the National Academy of Sciences, 2007, 104(12): 4 979 - 4 983.

    [13] Ackerly D. Conservatism and diversification of plant functional traits: evolutionary rates versus phylogenetic signal[J]. Proceedings of the National Academy of Sciences, 2009, 106(Supplement 2): 19 699 - 19 706.

    [14] Rafferty N E, Nabity P D. A global test for phylogenetic signal in shifts in flowering time under climate change[J]. Journal of Ecology, 2017, 105(3): 627 - 633.

    [15] Frishkoff L O, Valpine P, M'Gonigle L K. Phylogenetic occupancy models integrate imperfect detection and phylogenetic signal to analyze community structure[J]. Ecology, 2017, 98(1): 198 - 210.

    [16] Pagel M. Inferring the historical patterns of biological evolution[J]. Nature, 1999, 401(6756): 877.

    [17] Blomberg S P, Garland Jr T, Ives A R. Testing for phylogenetic signal in comparative data: behavioral traits are more labile[J]. Evolution, 2003, 57(4): 717 - 745.

    [18] DinizFilho F, Alexandre J, Rangel T F, et al. Exploring patterns of interspecific variation in quantitative traits using sequential phylogenetic eigenvector regressions[J]. Evolution, 2012, 66(4): 1 079 - 1 090.

    [19] Fritz S A, Purvis A. Selectivity in mammalian extinction risk and threat types: a new measure of phylogenetic signal strength in binary traits[J]. Conservation Biology, 2010, 24(4): 1 042 - 1 051.

    [20] Díaz S, Purvis A, Cornelissen J H C, et al. Functional traits, the phylogeny of function, and ecosystem service vulnerability[J]. Ecology and evolution, 2013, 3(9): 2 958 - 2 975.

    [21] Jetz W, Thomas G H, Joy J B, et al. Global distribution and conservation of evolutionary distinctness in birds[J]. Current Biology, 2014, 24(9): 919 - 930.

    [22] 宋朝樞. ?寶天曼自然保護區(qū)科學(xué)考察集[J]. ?北京: 中國林業(yè)出版社, 1994, 1(996): 2.

    [23] 朱從波, ?王萬里, ?劉曉靜, ?等. ?寶天曼自然保護區(qū)珍稀瀕危植物研究[J]. ?安徽農(nóng)業(yè)科學(xué), ?2011a, 39(16): ?9 467 - 9 470.

    [24] 盧炯林, ?王磐基, ?王正用. ?河南省寶天曼自然保護區(qū)珍稀瀕危植物資源的初步研究[J]. ?河南大學(xué)學(xué)報 (自然科學(xué)版),1988, ?3: ?35 - 41.

    [25] 王磐基, ?楊永芳, ?胡彧. ?河南伏牛山區(qū)植物多樣性及保護研究[J]. ?面向 21 世紀(jì)的中國生物多樣性保護——第三屆全國生物多樣性保護與持續(xù)利用研討會論文集, ?1998, 198 - 207.

    [26] 朱從波, ?李丹丹, ?劉曉靜, ?等. ?河南省 5 個自然保護區(qū)珍稀瀕危植物的比較研究[J]. ?福建林業(yè)科技, ?2011b, ?38(3): ?127 - 130.

    [27] 董云嵐, ?李愛琴, ?丁向陽, ?等. ?寶天曼自然保護區(qū)生物資源及保護利用[J]. ?河南科學(xué), ?2002, ?20(1): ?61 - 65.

    [28] 河南省林業(yè)廳. ?河南省重點保護植物名錄[R]. ?鄭州: ?河南省人民政府, ?2005.

    [29] 王文靜, ?宋曉潔. ?寶天曼自然保護區(qū)珍稀瀕危植物的概況及區(qū)系分析[J]. ?河南師范大學(xué)學(xué)報 (自然科學(xué)版), ?2001, ?29(2): ?115 - 117.

    [30] 孫偉娟, ?梅光朋, ?李家美. ?河南省重點保護植物名錄修訂[J]. 河南農(nóng)業(yè)大學(xué)學(xué)報, ?2011, ?45(3): ?358 - 361.

    [31] Webb C O, Donoghue M J. Phylomatic: tree assembly forapplied phylogenetics[J]. Molecular Ecology Resources, 2005, 5(1): 181 - 183.

    [32] 閆滿玉, ?杜曉軍, ?趙愛花, ?等. ?河南寶天曼落葉闊葉林木本植物單物種-面積關(guān)系[J]. ?生物多樣性, ?2015, ?23(5): 630 - 640.

    [33] Angiosperm Phylogeny Group. An update of the Angiosperm Phylogeny Group classification for the orders and families of flowering plants: APG III[J]. Botanical Journal of the Linnean Society, 2009, 161(2): 105 - 121.

    [34] Ackerly, David. "Conservatism and diversification of plantfunctional traits: evolutionary rates versus phylogenetic signal. "Proceedings of the National Academy of Sciences 106.Supplement 2 (2009): 19 699 - 19 706.

    [35] Mantel N. The detection of disease clustering and a generalized regression approach[J]. Cancer research, 1967, 27(2 Part 1): 209 - 220.

    [36] Pillar V D, Duarte L S. A framework for metacommunity ?analysis of phylogenetic structure[J]. Ecology letters, 2010, 13(5): 587 - 596.

    [37] Legendre P, FORTIN M J. Comparison of the Mantel test andalternative approaches for detecting complex multivariate relationships in the spatial analysis of genetic data[J]. Molecular ecology resources, 2010, 10(5): 831 - 844.

    [38] Hardy O J, Pavoine S. Assessing phylogenetic signal withmeasurement error: a comparison of mantel tests, blomberg et al.'sk, and phylogenetic distograms[J]. Evolution, 2012, 66(8): 2 614 - 2 621.

    [39] Hansen T F. Stabilizing selection and the comparative analysis of adaptation[J]. Evolution, 1997, 51(5): 1 341 - 1 351.

    [40] Harmon L J, Glor R E. Poor statistical performance of the Mantel test in phylogenetic comparative analyses[J]. Evolution, 2010, 64(7): 2 173 - 2 178.

    [41] Kraft N J B, Cornwell W K, Webb C O, et al. Trait evolution, community assembly, and the phylogenetic structure of ecological communities[J]. The American Naturalist, 2007, 170(2): 271 - 283.

    第1作者簡介: ?閆滿玉(1988-), ?男, ?碩士, ?研究方向為植被生態(tài)學(xué)。

    通訊作者: ?齊 光(1983-), ?男, ?博士, ?副教授, ?研究方向為植被生態(tài)學(xué)。

    收稿日期: 2020 - 01 - 20

    (責(zé)任編輯: ? 李 丹)

    Abstract Species closely related to be more similar to each other than to more distantly related ones, the calculation of phylogenetic signal play a important role in phylogenetic conservatism and species assembly patterns. Based on D and EM-mantel that could estimate phylogenetic signal for binary traits, we estimated the phylogenetic signal of the key preserved angiosperm's different life forms in Baotianman Nature Reserve, Henan Province. Our results showed that the arbor has significant phylogenetic signal(p=0.001), while the values of shrub and herbaceous were very different in terms of D and EM-mantle. Our results supported the extinction risk of arbor was dominated by habitat filtering, and those of shrub and herbaceous were human disturbance.

    Key words Endangered plants; ?Phylogenetic signal; ?Species coexistence; ?Baotianman Nature Reserve (BNR)

    热re99久久国产66热| 一本久久精品| www.av在线官网国产| 伦理电影大哥的女人| 尾随美女入室| 一级毛片aaaaaa免费看小| 国产欧美日韩一区二区三区在线 | 桃花免费在线播放| 七月丁香在线播放| 一边亲一边摸免费视频| 一本—道久久a久久精品蜜桃钙片| 美女主播在线视频| 亚洲欧美成人综合另类久久久| 亚洲av中文av极速乱| 黄色日韩在线| 久久女婷五月综合色啪小说| 人人妻人人爽人人添夜夜欢视频 | 自拍偷自拍亚洲精品老妇| 男女啪啪激烈高潮av片| 高清黄色对白视频在线免费看 | 熟女人妻精品中文字幕| 乱系列少妇在线播放| 国内精品宾馆在线| 精品一区二区三区视频在线| 国产美女午夜福利| 99热6这里只有精品| 国产视频内射| 欧美精品高潮呻吟av久久| 久久女婷五月综合色啪小说| 国产成人免费观看mmmm| 卡戴珊不雅视频在线播放| 超碰97精品在线观看| 久久婷婷青草| 亚洲精品自拍成人| 一级毛片电影观看| 国产视频内射| 欧美激情极品国产一区二区三区 | 91久久精品电影网| 久久人人爽人人爽人人片va| 2021少妇久久久久久久久久久| 极品少妇高潮喷水抽搐| 国产熟女欧美一区二区| 亚洲伊人久久精品综合| 国模一区二区三区四区视频| 国产日韩欧美亚洲二区| 国产免费福利视频在线观看| 国产精品免费大片| 亚洲欧美日韩另类电影网站| 日韩在线高清观看一区二区三区| 如日韩欧美国产精品一区二区三区 | 日韩欧美一区视频在线观看 | 少妇 在线观看| 成人综合一区亚洲| 国产成人精品无人区| 各种免费的搞黄视频| 欧美日韩国产mv在线观看视频| 国产爽快片一区二区三区| 国产伦在线观看视频一区| 蜜臀久久99精品久久宅男| 久久久精品94久久精品| 热re99久久精品国产66热6| 亚洲精华国产精华液的使用体验| 深夜a级毛片| 一级二级三级毛片免费看| 国产日韩欧美视频二区| 亚洲在久久综合| av在线观看视频网站免费| 中国国产av一级| 久久人人爽人人片av| 蜜桃在线观看..| 亚洲一区二区三区欧美精品| 国产亚洲最大av| 狂野欧美白嫩少妇大欣赏| 秋霞伦理黄片| 91aial.com中文字幕在线观看| 人人澡人人妻人| 天堂俺去俺来也www色官网| 91精品一卡2卡3卡4卡| a级毛片在线看网站| 久久久精品94久久精品| 香蕉精品网在线| av有码第一页| 成人国产av品久久久| 久久久久精品性色| 青春草视频在线免费观看| 国产免费一级a男人的天堂| 久久久久久久久大av| 亚洲国产最新在线播放| 日韩欧美一区视频在线观看 | 中文字幕亚洲精品专区| 亚洲经典国产精华液单| 日本vs欧美在线观看视频 | 汤姆久久久久久久影院中文字幕| 啦啦啦中文免费视频观看日本| 天堂俺去俺来也www色官网| av在线老鸭窝| 久久久久精品性色| 丰满迷人的少妇在线观看| 免费久久久久久久精品成人欧美视频 | 亚洲综合精品二区| 99热全是精品| 日韩视频在线欧美| 永久网站在线| 日日撸夜夜添| 六月丁香七月| 国产亚洲一区二区精品| 精品视频人人做人人爽| 中文字幕人妻丝袜制服| 欧美日韩精品成人综合77777| 亚洲精品国产色婷婷电影| 国模一区二区三区四区视频| 国产91av在线免费观看| 九草在线视频观看| 欧美日韩视频精品一区| 成人特级av手机在线观看| 一本大道久久a久久精品| 高清在线视频一区二区三区| 亚洲国产精品999| 我的老师免费观看完整版| 精华霜和精华液先用哪个| 草草在线视频免费看| 亚洲精品国产色婷婷电影| 国产精品嫩草影院av在线观看| 蜜桃久久精品国产亚洲av| 色5月婷婷丁香| 色网站视频免费| 狠狠精品人妻久久久久久综合| 久久这里有精品视频免费| 黑丝袜美女国产一区| 中文天堂在线官网| 在线观看三级黄色| 亚洲欧美清纯卡通| 亚洲成人一二三区av| 九九爱精品视频在线观看| 日韩在线高清观看一区二区三区| 伊人亚洲综合成人网| 久久久久精品久久久久真实原创| 欧美变态另类bdsm刘玥| 男女无遮挡免费网站观看| 亚洲av二区三区四区| 高清黄色对白视频在线免费看 | 曰老女人黄片| 天堂8中文在线网| 亚洲va在线va天堂va国产| 国产视频内射| 久久久精品免费免费高清| 国产av精品麻豆| a级毛片在线看网站| 国产精品成人在线| 亚洲精品乱久久久久久| 在线观看免费高清a一片| 黄色视频在线播放观看不卡| 中文字幕亚洲精品专区| 青春草视频在线免费观看| 久久久欧美国产精品| 久久女婷五月综合色啪小说| 男女免费视频国产| 国产精品国产三级专区第一集| 亚洲国产av新网站| 国产精品免费大片| videossex国产| 一级黄片播放器| 国产欧美亚洲国产| 人妻夜夜爽99麻豆av| 亚洲精品,欧美精品| 99久久中文字幕三级久久日本| 一区二区三区四区激情视频| 亚洲欧美一区二区三区国产| 欧美区成人在线视频| 我的老师免费观看完整版| 日韩精品有码人妻一区| 中国国产av一级| 国产 精品1| 最新中文字幕久久久久| 亚洲精品乱久久久久久| 丰满迷人的少妇在线观看| 老司机亚洲免费影院| 曰老女人黄片| 日日啪夜夜爽| 丰满人妻一区二区三区视频av| 国产日韩一区二区三区精品不卡 | av播播在线观看一区| 婷婷色综合大香蕉| 亚洲国产精品国产精品| 在线观看免费视频网站a站| 国产白丝娇喘喷水9色精品| 久久久久精品性色| 一级av片app| 晚上一个人看的免费电影| 国产精品久久久久久精品古装| 午夜视频国产福利| 麻豆成人av视频| 亚洲av欧美aⅴ国产| 26uuu在线亚洲综合色| 男女边吃奶边做爰视频| 亚洲欧美一区二区三区黑人 | 校园人妻丝袜中文字幕| av女优亚洲男人天堂| 国产在线一区二区三区精| 国产在线免费精品| 麻豆成人午夜福利视频| 欧美日韩亚洲高清精品| 在线观看免费视频网站a站| 精品国产一区二区三区久久久樱花| 又大又黄又爽视频免费| 一区二区av电影网| 三级国产精品欧美在线观看| 菩萨蛮人人尽说江南好唐韦庄| 欧美国产精品一级二级三级 | 亚洲中文av在线| 只有这里有精品99| 国产极品天堂在线| 桃花免费在线播放| 久久精品久久精品一区二区三区| 夜夜骑夜夜射夜夜干| 久久韩国三级中文字幕| 啦啦啦啦在线视频资源| 日韩精品免费视频一区二区三区 | 精品一区二区免费观看| 亚洲精品乱久久久久久| 一区二区三区免费毛片| 国产乱来视频区| 午夜视频国产福利| 亚洲av日韩在线播放| 午夜影院在线不卡| 亚洲无线观看免费| 高清黄色对白视频在线免费看 | 国产视频内射| 又黄又爽又刺激的免费视频.| 最近中文字幕高清免费大全6| 午夜老司机福利剧场| 99久久人妻综合| 欧美人与善性xxx| 如日韩欧美国产精品一区二区三区 | av福利片在线观看| 国产黄频视频在线观看| 成人特级av手机在线观看| 欧美日韩国产mv在线观看视频| 久久影院123| 婷婷色综合www| 内地一区二区视频在线| 日本黄色片子视频| 麻豆成人av视频| 少妇高潮的动态图| 亚洲精品亚洲一区二区| 一级a做视频免费观看| 久久韩国三级中文字幕| 亚洲不卡免费看| 老司机影院毛片| 亚洲av男天堂| 午夜福利影视在线免费观看| 日本欧美视频一区| 男女啪啪激烈高潮av片| 亚洲国产精品一区三区| 男人狂女人下面高潮的视频| 国产精品人妻久久久久久| 国产探花极品一区二区| 国产黄频视频在线观看| 亚洲情色 制服丝袜| 麻豆精品久久久久久蜜桃| 蜜桃在线观看..| 国产午夜精品久久久久久一区二区三区| 成人特级av手机在线观看| 99热这里只有是精品在线观看| 人妻一区二区av| 免费人妻精品一区二区三区视频| 激情五月婷婷亚洲| 亚洲精品成人av观看孕妇| 久久久久久久久久成人| 国产熟女午夜一区二区三区 | av国产精品久久久久影院| 丝瓜视频免费看黄片| 人人妻人人澡人人爽人人夜夜| 久久国产乱子免费精品| 夫妻午夜视频| 黄色视频在线播放观看不卡| 在线观看免费视频网站a站| 肉色欧美久久久久久久蜜桃| 丰满人妻一区二区三区视频av| av国产久精品久网站免费入址| 中文在线观看免费www的网站| 天堂中文最新版在线下载| 高清毛片免费看| 亚洲欧美成人精品一区二区| 亚洲av日韩在线播放| 高清av免费在线| 少妇 在线观看| kizo精华| 亚洲美女黄色视频免费看| 成人特级av手机在线观看| 国产视频内射| 女性生殖器流出的白浆| 色婷婷av一区二区三区视频| 高清欧美精品videossex| 赤兔流量卡办理| 中国三级夫妇交换| 97在线人人人人妻| 人人妻人人澡人人看| 日本黄色片子视频| 一区二区三区精品91| 国产老妇伦熟女老妇高清| 三级经典国产精品| 国产亚洲最大av| 免费黄网站久久成人精品| 女性生殖器流出的白浆| 亚洲av成人精品一二三区| 成人无遮挡网站| 国产精品人妻久久久久久| 啦啦啦啦在线视频资源| av视频免费观看在线观看| 久久99一区二区三区| 成人特级av手机在线观看| 国产精品久久久久久精品古装| 黄色一级大片看看| 少妇人妻久久综合中文| 成人无遮挡网站| 中文字幕久久专区| 一区二区三区精品91| 午夜91福利影院| 国产伦在线观看视频一区| 亚洲精品乱久久久久久| 五月天丁香电影| 国产免费又黄又爽又色| 老司机亚洲免费影院| 在线观看免费视频网站a站| 亚洲真实伦在线观看| 亚洲欧美日韩卡通动漫| av黄色大香蕉| 国产在线免费精品| 女人精品久久久久毛片| 熟女人妻精品中文字幕| 欧美人与善性xxx| 国产一区二区在线观看av| 日本欧美国产在线视频| 亚洲国产精品999| 最黄视频免费看| 日韩欧美精品免费久久| 男女无遮挡免费网站观看| 久久99蜜桃精品久久| 99久久综合免费| 精品一品国产午夜福利视频| 亚洲熟女精品中文字幕| 男人爽女人下面视频在线观看| av.在线天堂| 精品国产国语对白av| 日韩,欧美,国产一区二区三区| 在线观看国产h片| h视频一区二区三区| 久久久久视频综合| 国产一区有黄有色的免费视频| 99久久精品热视频| 国产爽快片一区二区三区| 日本午夜av视频| 免费看av在线观看网站| 如日韩欧美国产精品一区二区三区 | 黄色日韩在线| av免费在线看不卡| 亚洲欧美日韩另类电影网站| 日日摸夜夜添夜夜添av毛片| 欧美日韩视频精品一区| 精品少妇久久久久久888优播| 亚洲国产精品国产精品| 午夜福利网站1000一区二区三区| 有码 亚洲区| 香蕉精品网在线| 成人亚洲精品一区在线观看| 18禁裸乳无遮挡动漫免费视频| 男人爽女人下面视频在线观看| 少妇精品久久久久久久| 午夜福利在线观看免费完整高清在| xxx大片免费视频| 亚洲精品第二区| 久久久久久久久久人人人人人人| 激情五月婷婷亚洲| 丰满饥渴人妻一区二区三| 人人妻人人澡人人爽人人夜夜| 亚洲内射少妇av| 特大巨黑吊av在线直播| 永久免费av网站大全| 国产综合精华液| 亚洲av成人精品一区久久| 深夜a级毛片| 欧美日韩av久久| 人妻制服诱惑在线中文字幕| 亚洲内射少妇av| 婷婷色av中文字幕| 国产成人精品无人区| 免费观看性生交大片5| 国产色爽女视频免费观看| 中文欧美无线码| 国产精品久久久久成人av| 亚洲,欧美,日韩| 在线观看三级黄色| 妹子高潮喷水视频| 看十八女毛片水多多多| 最近中文字幕2019免费版| 亚洲精品国产色婷婷电影| 少妇人妻精品综合一区二区| 精品少妇黑人巨大在线播放| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲精品第二区| 久久久久国产网址| 一区二区av电影网| 欧美xxxx性猛交bbbb| 久久精品国产亚洲av涩爱| 亚洲av福利一区| 久久精品夜色国产| 少妇高潮的动态图| 大码成人一级视频| 特大巨黑吊av在线直播| 纵有疾风起免费观看全集完整版| av国产精品久久久久影院| 两个人免费观看高清视频 | 下体分泌物呈黄色| .国产精品久久| 亚洲精品国产成人久久av| 蜜桃久久精品国产亚洲av| 美女大奶头黄色视频| 亚洲成人手机| 久久av网站| 九色成人免费人妻av| 久久久久久久久久久丰满| 在线观看免费日韩欧美大片 | 亚洲精品国产av成人精品| 黄色日韩在线| 亚洲av免费高清在线观看| 国产免费福利视频在线观看| 午夜激情福利司机影院| 国产熟女午夜一区二区三区 | .国产精品久久| 婷婷色综合www| 我的老师免费观看完整版| 最近的中文字幕免费完整| 日本黄大片高清| 最近中文字幕高清免费大全6| 80岁老熟妇乱子伦牲交| 五月伊人婷婷丁香| 水蜜桃什么品种好| 国产av精品麻豆| 人人澡人人妻人| 少妇精品久久久久久久| 日韩不卡一区二区三区视频在线| 水蜜桃什么品种好| 插逼视频在线观看| av在线播放精品| 久久国产精品男人的天堂亚洲 | 久久热精品热| 男人和女人高潮做爰伦理| 精品人妻熟女毛片av久久网站| 国产精品人妻久久久影院| 自线自在国产av| 日本午夜av视频| 欧美日韩亚洲高清精品| 中文天堂在线官网| 热re99久久精品国产66热6| 成人免费观看视频高清| 精品少妇久久久久久888优播| 男女无遮挡免费网站观看| 色吧在线观看| 纯流量卡能插随身wifi吗| 日日啪夜夜撸| 国产 精品1| 在线观看免费日韩欧美大片 | 十八禁网站网址无遮挡 | 亚洲av电影在线观看一区二区三区| 欧美xxxx性猛交bbbb| 深夜a级毛片| 久久精品国产鲁丝片午夜精品| 亚洲国产欧美在线一区| 国产精品不卡视频一区二区| 男男h啪啪无遮挡| 日本免费在线观看一区| 日韩 亚洲 欧美在线| 久久99蜜桃精品久久| 熟妇人妻不卡中文字幕| 久久久久视频综合| 免费看不卡的av| 国产av国产精品国产| 人体艺术视频欧美日本| 51国产日韩欧美| 亚洲精品日韩av片在线观看| 老女人水多毛片| 亚洲精品国产成人久久av| 成年人免费黄色播放视频 | 伊人久久精品亚洲午夜| 最近手机中文字幕大全| 亚洲精品乱码久久久久久按摩| 国产一级毛片在线| 18禁裸乳无遮挡动漫免费视频| 80岁老熟妇乱子伦牲交| 亚洲欧美日韩卡通动漫| 国产成人精品久久久久久| 青春草国产在线视频| 26uuu在线亚洲综合色| 欧美在线黄色| 久久ye,这里只有精品| 国产在线免费精品| www日本在线高清视频| 日韩大码丰满熟妇| 成年av动漫网址| 亚洲第一欧美日韩一区二区三区 | 亚洲中文av在线| 日韩免费高清中文字幕av| 免费看十八禁软件| 黑人巨大精品欧美一区二区蜜桃| 老司机福利观看| 亚洲成av片中文字幕在线观看| 咕卡用的链子| 亚洲精品日韩在线中文字幕| 91字幕亚洲| 婷婷丁香在线五月| 十八禁网站网址无遮挡| 亚洲va日本ⅴa欧美va伊人久久 | 欧美日韩av久久| 国产欧美日韩一区二区精品| 亚洲欧美色中文字幕在线| 国产高清国产精品国产三级| 国产成人精品在线电影| 亚洲专区国产一区二区| 国产精品秋霞免费鲁丝片| 久久亚洲精品不卡| 国产成人免费无遮挡视频| 国产成人精品久久二区二区91| 啦啦啦啦在线视频资源| 少妇的丰满在线观看| 亚洲av美国av| 高清欧美精品videossex| 在线观看免费日韩欧美大片| 丁香六月天网| 亚洲国产精品一区三区| 久久中文字幕一级| 激情视频va一区二区三区| 女人精品久久久久毛片| 精品乱码久久久久久99久播| 亚洲欧美一区二区三区久久| 十八禁人妻一区二区| 久久人人97超碰香蕉20202| 午夜久久久在线观看| 亚洲中文av在线| 黄色视频在线播放观看不卡| bbb黄色大片| 日韩精品免费视频一区二区三区| 久久 成人 亚洲| 天天操日日干夜夜撸| 午夜精品国产一区二区电影| 亚洲avbb在线观看| 国产精品.久久久| 777久久人妻少妇嫩草av网站| bbb黄色大片| 80岁老熟妇乱子伦牲交| 99久久99久久久精品蜜桃| 久久久久国内视频| 精品亚洲成a人片在线观看| 91大片在线观看| 欧美激情高清一区二区三区| 亚洲精品国产区一区二| 老司机福利观看| 国产一区二区激情短视频 | 国产av又大| 日韩视频在线欧美| 国产成人免费无遮挡视频| av超薄肉色丝袜交足视频| 男人舔女人的私密视频| 午夜久久久在线观看| 91成年电影在线观看| 一二三四在线观看免费中文在| 精品国产一区二区久久| 国产91精品成人一区二区三区 | 欧美av亚洲av综合av国产av| xxxhd国产人妻xxx| 交换朋友夫妻互换小说| 亚洲全国av大片| 伦理电影免费视频| 香蕉丝袜av| 天堂俺去俺来也www色官网| 国产区一区二久久| 国产伦理片在线播放av一区| 亚洲va日本ⅴa欧美va伊人久久 | 一二三四在线观看免费中文在| 老司机午夜十八禁免费视频| kizo精华| 99国产精品一区二区三区| 欧美人与性动交α欧美精品济南到| 亚洲国产看品久久| 丁香六月欧美| 丝袜人妻中文字幕| 91成人精品电影| 脱女人内裤的视频| 亚洲精华国产精华精| 乱人伦中国视频| 9色porny在线观看| 亚洲一码二码三码区别大吗| 成人国产一区最新在线观看| 青春草亚洲视频在线观看| 日本a在线网址| 午夜视频精品福利| 国产精品欧美亚洲77777| 国产免费现黄频在线看| 欧美成人午夜精品| 亚洲成国产人片在线观看| 黄色视频不卡| 日韩 亚洲 欧美在线| 十分钟在线观看高清视频www| 高清视频免费观看一区二区| 欧美激情久久久久久爽电影 | 真人做人爱边吃奶动态| 亚洲精品自拍成人| 欧美精品高潮呻吟av久久| 国产成人免费无遮挡视频| 欧美亚洲 丝袜 人妻 在线| 日韩视频在线欧美| 王馨瑶露胸无遮挡在线观看| 久久九九热精品免费| 久久精品aⅴ一区二区三区四区| 久久精品国产综合久久久|