更吉卓瑪，李　彥，賈留坤，夏銘澤，高慶波

(1 中國科學(xué)院西北高原生物研究所 中國科學(xué)院高山植物適應(yīng)與進(jìn)化重點(diǎn)實(shí)驗(yàn)室，西寧 810001；2 中國科學(xué)院大學(xué)，北京 100039；3 中國科學(xué)院西北高原生物研究所 青海省作物分子育種重點(diǎn)實(shí)驗(yàn)室，西寧 810001)

由于復(fù)雜的地理拓?fù)浣Y(jié)構(gòu)、微環(huán)境以及更廣泛的生態(tài)位，高山地區(qū)往往具有較高的生物多樣性，如高加索山、安第斯山、阿爾卑斯山等[1-3]。青藏高原及其周邊的高山地區(qū)包含了3個(gè)全球生物多樣性熱點(diǎn)地區(qū)，即東部的橫斷山區(qū)、南部的喜馬拉雅山系及西部的阿爾泰和天山山脈[4]。一般認(rèn)為，與造山運(yùn)動(dòng)相關(guān)的地理隔離和異域分化是導(dǎo)致青藏高原地區(qū)植物多樣性較高的原因之一[5-6]?？焖佥椛浞只谇嗖馗咴貐^(qū)的很多植物類群中被檢測到，如馬先蒿屬(PedicularisL.)[7]、紅景天屬(RhodiolaL.)[8]、風(fēng)毛菊屬(SaussureaDC.)[9]、虎耳草屬(SaxifragaL.)[10-11]、槖吾屬-垂頭菊屬-蟹甲草屬復(fù)合群(LigulariaCass.-CremanthodiumBenth.-ParasenecioW. W. Sm. & J. Small)[12]等。然而，與較高分類階元的分化模式不同，青藏高原地區(qū)植物種內(nèi)居群間的分化可能更多的受到第四紀(jì)冰期氣候波動(dòng)的影響。更新世晚期以來，青藏高原經(jīng)歷了至少4次規(guī)模較大的冰期，其中聶聶雄拉冰期(第二冰期，約720～500 ka)強(qiáng)度最大，冰川擴(kuò)張最廣，此后冰期強(qiáng)度逐漸減弱[13-14]。而全球末次盛冰期時(shí)青藏高原冰川的擴(kuò)張程度非常有限，這與歐洲和北美大陸末次冰期形成統(tǒng)一的冰蓋不同[15]。第四紀(jì)冰期與間冰期交替所引起的氣候波動(dòng)對(duì)于現(xiàn)存生物類群的地理分布格局和遺傳結(jié)構(gòu)產(chǎn)生了極大的影響[16-17]。由于復(fù)雜的地理地貌特征、未被統(tǒng)一的大冰蓋覆蓋以及東亞季風(fēng)氣候和印度季風(fēng)氣候的影響，青藏高原地區(qū)的植物可能具有更復(fù)雜的第四紀(jì)冰期進(jìn)化歷史[18-19]。一方面，某些植物，如祁連圓柏(Juniperusprzewalskii)[20]、青海云杉(Piceacrassifolia)[21]等，在第四紀(jì)冰期時(shí)由高原臺(tái)面退縮到海拔較低的東南部邊緣避難所內(nèi)，間冰期或冰期后再從邊緣避難所回遷到高原臺(tái)面，從而形成現(xiàn)有的地理分布格局。另一方面，由于第四紀(jì)冰期時(shí)青藏高原未被統(tǒng)一的大冰蓋覆蓋，僅有不同程度的山地冰川發(fā)育，高原臺(tái)面可能存在適宜的局部微環(huán)境而成為某些耐寒植物的避難所。對(duì)銀露梅(Potentillaglabra)[22]、露蕊烏頭(Aconitumgymnandrum)[23]和西川紅景天(Rhodiolaalsia)[24]等的研究揭示，在第四紀(jì)冰期時(shí)高原臺(tái)面的居群并沒有退縮到高原邊緣避難所，而是在高原臺(tái)面上的一個(gè)或幾個(gè)較大的避難所中保留下來，間冰期或冰期后再由避難所向附近擴(kuò)散。此外，對(duì)西藏圓柏復(fù)合群(Juniperustibeticacomplex)[25]、西藏沙棘(Hippophaetibetana)[26]、菊葉紅景天(Rhodiolachrysanthemifolia)[27]等的研究發(fā)現(xiàn)，這些物種第四紀(jì)冰期時(shí)在高原臺(tái)面存在多個(gè)微型避難所，這些微型避難所在物種的分布范圍內(nèi)呈均勻分布。然而，具有不同分布范圍、生境和生活史特征的植物可能具有不同的第四紀(jì)冰期進(jìn)化歷史。對(duì)青藏高原地區(qū)更多的植物進(jìn)行譜系地理學(xué)研究，不但可為高山植物演化歷史與第四紀(jì)冰期的關(guān)系提供依據(jù)，同時(shí)對(duì)青藏高原生物多樣性保護(hù)具有重要啟示。

虎耳草屬是虎耳草科(Saxifragaceae)最大的屬，約有450～500種，主要分布在歐洲和亞洲的高山地區(qū)，某些物種環(huán)北極分布，另有一些物種擴(kuò)散到南美洲的落基山脈、安第斯山脈和火地島[10-11]。中國產(chǎn)虎耳草屬植物約220種，主要分布在青藏高原-喜馬拉雅地區(qū)[28]?；⒍輰僦参锉粡V泛地用于譜系地理學(xué)和系統(tǒng)發(fā)育學(xué)研究，來揭示北極和高山地區(qū)植物的分化式樣[10-11, 29-32]。然而，青藏高原作為虎耳草屬物種多樣性中心之一，該地區(qū)虎耳草屬植物種內(nèi)分化與第四紀(jì)冰期的關(guān)系尚不明確。唐古特虎耳草(S.tanguticaEngl.)為多年生草本植物，生于海拔2 900～4 600 m的高山草甸、灌叢、石隙，分布于青海、甘肅、四川、西藏等地，是青藏高原地區(qū)高山草甸的重要組成成分，對(duì)維持高寒草甸生態(tài)系統(tǒng)的穩(wěn)定具有重要意義。本研究選取葉綠體DNA(cpDNA)非編碼區(qū)trnL-trnF和rpl16片段對(duì)青藏高原地區(qū)的唐古特虎耳草進(jìn)行譜系地理學(xué)研究，目的在于揭示該物種居群內(nèi)與居群間的遺傳分布式樣，討論唐古特虎耳草對(duì)于第四紀(jì)氣候波動(dòng)的響應(yīng)，探討其現(xiàn)有地理分布格局的歷史成因。

1　材料和方法

1.1　實(shí)驗(yàn)材料

研究材料于2007年至2014年采自青海、西藏、四川、甘肅，共采集唐古特虎耳草18個(gè)居群209個(gè)個(gè)體(表1，圖1)。每個(gè)居群隨機(jī)采集2～24個(gè)個(gè)體，個(gè)體之間相隔至少10 m以上。采集生長良好植株的新鮮嫩葉后立即用硅膠干燥，帶回實(shí)驗(yàn)室存放于-20 ℃冰箱備用。憑證標(biāo)本由R. J. Gornall教授和高慶波副研究員共同鑒定，館藏于中國科學(xué)院西北高原植物研究所青藏高原生物標(biāo)本館(HNWP)。

表1　唐古特虎耳草18個(gè)居群的采集信息

餅狀圖表示每個(gè)居群中單倍型的頻率圖1　唐古特虎耳草18個(gè)居群的采樣圖以及葉綠體單倍型分布圖Pie charts show the proportions of haplotypes within each populationFig.1　Map of the 18 sampled populations of S. tangutica and the distribution of cpDNA haplotypes in the species

1.2　DNA提取與PCR擴(kuò)增

本研究采用改良的CTAB法從硅膠干燥的葉片中提取總DNA[33]。采用通用引物對(duì)所有個(gè)體的cpDNAtrnL-trnF和rpl16片段進(jìn)行PCR擴(kuò)增[34-35]。PCR擴(kuò)增反應(yīng)體系為50 μL：包含5 μL的10 × PCR Buffer(含1.5 mmol/L MgCl2)，2 μL的10 mmol/L dNTPs，正反引物各1.0 μL(5 pmol/L)，TaqDNA聚合酶(TaKaRa，大連)0.4 μL，1 μL的DNA模板，用雙蒸水補(bǔ)足到50 μL。PCR擴(kuò)增反應(yīng)程序：94 ℃預(yù)變性10 min；94 ℃變性1 min，58 ℃退火50 s，72 ℃延伸1 min，32循環(huán)；72 ℃延伸10 min。PCR產(chǎn)物用1%的瓊脂糖凝膠電泳進(jìn)行檢測，然后送測序公司進(jìn)行雙向測序。

1.3　數(shù)據(jù)分析

利用Chromas 2.6.2軟件(http://www.technelysium.com.au)對(duì)測序峰圖進(jìn)行肉眼核對(duì)。利用MEGA 7.0.26軟件對(duì)序列進(jìn)行對(duì)位排列并統(tǒng)計(jì)序列的堿基組成[36]，發(fā)現(xiàn)trnL-trnF序列上存在一段約20 bp的poly A結(jié)構(gòu)，為了確保堿基變異的可靠性，poly A及兩側(cè)約40 bp的堿基被刪除。同樣，對(duì)rpl16序列中出現(xiàn)的大于11 bp的poly A和poly T結(jié)構(gòu)也進(jìn)行了類似的刪除處理。利用DnaSP 5.10軟件統(tǒng)計(jì)cpDNA聯(lián)合序列的變異位點(diǎn)并確定單倍型[37]?？傔z傳多樣性(HT)、居群內(nèi)平均遺傳多樣性(Hs)、居群間遺傳分化系數(shù)GST和NST值均利用PERMUT軟件進(jìn)行計(jì)算[38]。其中GST僅考慮單倍型的頻率，而NST兼顧了單倍型頻率和單倍型之間的相似性。使用U-統(tǒng)計(jì)方法對(duì)GST和NST進(jìn)行比較(1 000次重復(fù)的置換檢驗(yàn))以檢測單倍型變異的地理分布模式。NST大于GST且顯著，表明居群在分布范圍內(nèi)具有顯著的譜系地理結(jié)構(gòu)。由于PERMUT軟件要求每個(gè)居群至少要有3個(gè)個(gè)體，所以，居群P13(2個(gè)個(gè)體)在進(jìn)行PERMUT分析時(shí)未包含在內(nèi)。

運(yùn)用Arlequin 3.5.2計(jì)算每個(gè)居群的單倍型多樣性(He)和核苷酸多樣性(π)[39]。利用Arlequin軟件包中的分子變異分析(AMOVA)檢測居群內(nèi)和居群間的遺傳變異水平，并對(duì)單倍型分布的FST進(jìn)行評(píng)價(jià)(1 000次置換檢驗(yàn))[40]。Tajima’sD、Fu & Li’sD*兩種無限突變位點(diǎn)模型的中性檢驗(yàn)及歧點(diǎn)分布分析都在DnaSP 5.10程序中完成[41-42]。

以最大簡約法為原則，利用NETWORK 4.6.1.3軟件構(gòu)建單倍型中央連接網(wǎng)狀圖[43-44]。在運(yùn)行該軟件時(shí)，對(duì)于多堿基的插入、缺失或突變位點(diǎn)，被認(rèn)為是通過一步突變形成。在單倍型Network圖中，原始的單倍型多存在于網(wǎng)絡(luò)圖的內(nèi)部，而近代起源的單倍型發(fā)生在網(wǎng)絡(luò)圖的邊緣位置[45-46]。

2　結(jié)果與分析

2.1　唐古特虎耳草的單倍型分布

對(duì)18個(gè)居群209個(gè)個(gè)體的trnL-trnF和rpl16片段進(jìn)行測序，MEGA軟件對(duì)位排列后的矩陣長度為1 626 bp，通過DnaSP軟件鑒定出74個(gè)單倍型，有25個(gè)簡約信息位點(diǎn)(表2)。在檢測出的74個(gè)單倍型中，有53個(gè)單倍型是居群特有單倍型，單倍型H5出現(xiàn)的頻率最高，有46個(gè)個(gè)體擁有此單倍型。表3顯示在18個(gè)居群中，16個(gè)居群擁有特有單倍型。居群的遺傳多樣性(He)范圍為0.490 9～1.000 0，平均為0.853 6；居群的核苷酸多樣性(π)范圍為0.007 506～0.046 481，平均為0.022 145。

2.2　居群遺傳結(jié)構(gòu)特征

通過PERMUT程序計(jì)算得出唐古特虎耳草居群內(nèi)平均遺傳多樣性HS值(0.843)、總遺傳多樣性HT值(0.933)、居群間遺傳分化GST(0.097)和NST值(0.109)。使用U-統(tǒng)計(jì)方法對(duì)唐古特虎耳草單倍型變異的地理結(jié)構(gòu)進(jìn)行檢驗(yàn)后發(fā)現(xiàn)NST大于GST，但不顯著(P> 0.05)，且居群間的遺傳分化水平很低，表明唐古特虎耳草在現(xiàn)有分布區(qū)域內(nèi)沒有明顯的譜系地理結(jié)構(gòu)。分子變異分析(AMOVA)的結(jié)果表明，居群間的遺傳變異僅為8.15%(FST= 0.081 45；P< 0.01)，而居群內(nèi)的遺傳變異為91.85%(表4)，進(jìn)一步揭示了唐古特虎耳草遺傳變異主要存在于居群內(nèi)。

2.3　基于cpDNA聯(lián)合序列的居群歷史動(dòng)態(tài)分析

對(duì)唐古特虎耳草cpDNA聯(lián)合序列所有個(gè)體進(jìn)行中性檢驗(yàn)，結(jié)果表明，Tajima’sD(-2.045 07，P< 0.05)和Fu & Li’D*(-3.629 27，P< 0.05)均為顯著的負(fù)值，推測其分布范圍或者個(gè)體數(shù)近期經(jīng)歷過擴(kuò)張。歧點(diǎn)分布分析結(jié)果顯示，在唐古特虎耳草整個(gè)分布范圍內(nèi)觀測到的歧點(diǎn)分布呈單峰曲線(圖2)，其SSD值(0.007 2)和Raggednes index值(0.007 6)均不顯著(P> 0.05)，說明該物種在整個(gè)分布范圍內(nèi)經(jīng)歷過近期擴(kuò)張現(xiàn)象，這與中性檢驗(yàn)結(jié)果一致。

2.4　單倍型之間的系統(tǒng)發(fā)育關(guān)系

最大簡約法(maximum parsimony, MP)、最大似然法(maximum likelihood, ML)、貝葉斯法(Bayesian inference, BI)構(gòu)建的單倍型系統(tǒng)發(fā)育樹均呈現(xiàn)“刷子”結(jié)構(gòu)，不能很好地揭示單倍型之間的系統(tǒng)發(fā)育關(guān)系，也表明了唐古特虎耳草單倍型的快速輻射分化(結(jié)果未顯示)。NETWORK軟件可更好揭示遺傳分化較淺序列之間的系統(tǒng)發(fā)育關(guān)系?；谧畲笏迫环?gòu)建74個(gè)單倍型的中央鏈接網(wǎng)狀圖，結(jié)果顯示(圖3)，居群共享單倍型，如H1、H2、H3、H5、H16等位于網(wǎng)狀圖的中央位置，而居群特有單倍型一般位于網(wǎng)狀圖的邊緣。根據(jù)溯祖理論，共享單倍型為較古老的單倍型，而居群特有單倍型可能是通過近期輻射分化衍生而來的較年輕的單倍型。

表2　唐古特虎耳草74個(gè)單倍型trnL-trnF和rpl16片段簡約信息位點(diǎn)

續(xù)表2Continued Table 2

單倍型Hap簡約信息位點(diǎn)ParsimonyinformativesitetrnL-trnFrpl16230355414456611792813928966988100010901124119213341352135913601414145114571465147915511562H53-------------------------H54--------------?A----------H55---------?C--?A------------H56--------?T--------?GC-----?TH57-----------?T-?A------?A---?TH58-------------------------H59-------------------------H60----------------?A--------H61----------?A--------------H62-------------------------H63---------?C--?A------------H64--------?T---------------?TH65-C-----C-?C---------------H66----C------?T-----------?GTH67--------?T--?T-------?AA-?A-?TH68-------------------------H69-C----A-------?A----------H70--------------?A----------H71-------------------------H72A----T------------------?TH73----------------?A--------H74----------?A--------------

注：-代表與單倍型H1相同的堿基

Note: - stands for the same base ashaplotypeH1

表3　唐古特虎耳草18個(gè)居群的單倍型組成、遺傳多樣性和核苷酸多樣性

表4　唐古特虎耳草18個(gè)居群的分子遺傳變異(AMOVA)分析結(jié)果

圖2　唐古特虎耳草209個(gè)個(gè)體葉綠體DNA2個(gè)片段聯(lián)合數(shù)據(jù)的歧點(diǎn)分布分析Fig.2　Mismatch distribution for 209 individuals of S. tangutica based on two cpDNA fragments

3　討　論

3.1　唐古特虎耳草的居群遺傳結(jié)構(gòu)

遺傳變異分析表明唐古特虎耳草cpDNA片段總遺傳多樣性水平較高，HT= 0.933。類似的結(jié)果也在青藏高原很多高山植物中被檢測到，如偏花報(bào)春(Primulasecundiflora,HT= 0.966)[47]、西藏沙棘(Hippophaetibetana,HT= 0.956)[26]、綿參(Eriophytonwallichii,HT= 0.979)[48]、西川紅景天(Rhodiolaalsia,HT= 0.950)[24]、菊葉紅景天(Rhodiolachrysanthemifolia,HT= 0.937)[27]等。然而，這些物種大多數(shù)都具有低的居群內(nèi)平均遺傳多樣性。在本研究中，唐古特虎耳草不僅具有較高的總遺傳多樣性，居群內(nèi)的平均遺傳多樣性也很高(HS= 0.843)，表明唐古特虎耳草較高的遺傳變異主要由居群內(nèi)遺傳變異所致。這一結(jié)果也得到了分子變異分析(AMOVA)的證實(shí)。AMOVA分析表明，唐古特虎耳草居群內(nèi)遺傳變異占總變異的91.85%，而居群間遺傳變異僅為8.15%，進(jìn)一步揭示了唐古特虎耳草遺傳變異主要存在于居群內(nèi)。如此高的居群內(nèi)遺傳變異在青藏高原的植物中并不多見，推測造成該遺傳結(jié)構(gòu)的原因如下：一方面，由于檢測到的74個(gè)單倍型中只有少數(shù)共享單倍型，但其個(gè)體數(shù)所占比例很高(67.46%)，并且這些共享單倍型在居群中隨機(jī)分布，因此降低了居群間的遺傳差異；另一方面，居群特有單倍型比例高(53/74)，某些居群固定多個(gè)特有單倍型，從而增加了居群內(nèi)的遺傳變異。特有單倍型比例較高的物種往往不能揭示出明顯的譜系地理格局[22, 24, 27]。 此外，對(duì)虎耳草屬山羊臭組(Saxifragasect.CiliataeHaw.)系統(tǒng)發(fā)育和物種形成的研究表明，唐古特虎耳草所在的唐古拉亞組(subsect.HirculoideaeEngl. & Irmsch.)具有較快的分化速率，該亞組約110個(gè)物種是在最近的2.12百萬年(Ma)分化而來[10]。由此推測，唐古特虎耳草具有更短的居群分化歷史。事實(shí)上，基于NETWORK構(gòu)建的單倍型系統(tǒng)發(fā)育關(guān)系表明，大量的特有單倍型位于NETWORK圖的邊緣，揭示了唐古特虎耳草單倍型的近期快速分化。在第四紀(jì)冰期前，唐古特虎耳草可能廣泛分布于青藏高原及其周邊地區(qū)，并廣泛固定某些古老單倍型(如H5)。隨著第四紀(jì)冰期氣候的波動(dòng)，原來連續(xù)的居群片段化，形成了相互隔離的居群，并經(jīng)歷了異域分化，產(chǎn)生大量特有單倍型。在此過程中，瓶頸效應(yīng)和遺傳漂變等作用使得特有單倍型在居群中固定下來，抹去了居群原有的部分遺傳結(jié)構(gòu)。由于單獨(dú)進(jìn)化時(shí)間較短，特有單倍型之間的差異大多表現(xiàn)在單個(gè)堿基的變異。

圓的大小表示單倍型頻率，實(shí)心方塊代表缺失單倍型圖3　唐古特虎耳草cpDNA聯(lián)合片段74種單倍型的中央連接網(wǎng)狀圖Circle size is proportional to haplotype frequencies, parallelograms represent missing haplotypesFig.3　Median-joining network of 74 haplotypes of S. tangutica based on two cpDNA fragments

3.2　冰期避難所及居群進(jìn)化歷史的推測

第四紀(jì)循環(huán)的氣候波動(dòng)對(duì)生物區(qū)系的分布和遺傳結(jié)構(gòu)產(chǎn)生了巨大的影響[49]。植物對(duì)第四紀(jì)冰期的反應(yīng)可以根據(jù)現(xiàn)有的遺傳結(jié)構(gòu)進(jìn)行追溯。于海彬等[50]綜述了36種高山植物的譜系地理格局，總結(jié)出高山植物面對(duì)青藏高原復(fù)雜的地理結(jié)構(gòu)、微氣候和生境條件所呈現(xiàn)出的譜系歷史。一部分高山植物在第四紀(jì)冰期時(shí)退縮到青藏高原邊緣避難所，冰期后回遷到高原臺(tái)面。符合這種“退縮-回遷”模式的物種，其遺傳結(jié)構(gòu)表現(xiàn)為遺傳多樣性較高的居群和古老單倍型往往集中在高原邊緣避難所內(nèi)，居群遺傳多樣性由高原邊緣向高原臺(tái)面遞減。此外，一些耐寒植物，在冰期時(shí)沒有完全退縮到低海拔的邊緣地區(qū)，而是在高原臺(tái)面形成了一個(gè)或幾個(gè)相互隔離的避難所，冰期后發(fā)生局部擴(kuò)張。符合這種“高原臺(tái)面避難所-局部擴(kuò)張”模式的物種，其遺傳結(jié)構(gòu)表現(xiàn)為特有單倍型或者遺傳多樣性較高的居群分布在高原臺(tái)面上的一個(gè)或幾個(gè)相互隔離的區(qū)域。另外，某些物種第四紀(jì)冰期時(shí)在高原臺(tái)面形成多個(gè)微型避難所，間冰期或冰期后發(fā)生小范圍的擴(kuò)張，其遺傳多樣性較高的居群和特有單倍型在物種的分布范圍內(nèi)均勻分布。

本研究沒有揭示出唐古特虎耳草遺傳多樣性較高的某一地理區(qū)域，也沒有發(fā)現(xiàn)遺傳多樣性由高原邊緣向高原臺(tái)面遞減。相反，本研究所涉及的唐古特虎耳草居群普遍具有較高的遺傳多樣性，特有單倍型在分布范圍內(nèi)均勻分布。我們推測，由于青藏高原地形的高度異質(zhì)性以及復(fù)雜的局部微氣候和生境條件，在第四紀(jì)冰期時(shí)，唐古特虎耳草在高原臺(tái)面上未被冰川覆蓋的多個(gè)微型避難所中存活下來。唐古特虎耳草居群對(duì)第四紀(jì)冰期氣候波動(dòng)的反應(yīng)可能更多地表現(xiàn)為垂直海拔高度的遷移，而非大規(guī)模的水平遷移，這與銀露梅(Potentillaglabra)[22]和西川紅景天(Rhodiolachrysanthemifolia)[27]的進(jìn)化歷史相似。居群在垂直遷移的過程中發(fā)生了瓶頸效應(yīng)和奠基者效應(yīng)，居群隨機(jī)固定少量古老單倍型。隨后的異域分化和遺傳漂變使得較年輕的衍生特有單倍型在居群中固定下來，形成了唐古特虎耳草現(xiàn)有的居群遺傳結(jié)構(gòu)。居群歷史動(dòng)態(tài)檢驗(yàn)結(jié)果表明，Tajima’sD和Fu & Li’D*均為顯著的負(fù)值，歧點(diǎn)分布分析呈單峰曲線，說明唐古特虎耳草經(jīng)歷過近期分布范圍或者個(gè)體數(shù)的擴(kuò)張[51-53]。物種經(jīng)歷分布范圍擴(kuò)張事件，其原分布中心應(yīng)具有較高的遺傳多樣性和特有單倍型，并且原分布中心的居群包含擴(kuò)散區(qū)域居群所具有的單倍型。但在本研究中，大量的特有單倍型在唐古特虎耳草的整個(gè)分布范圍內(nèi)均勻分布，遺傳多樣性和核苷酸多樣性較高的居群也呈現(xiàn)均勻分布，并未檢測到分布范圍擴(kuò)張的信號(hào)。由此推測，該物種可能經(jīng)歷過個(gè)體數(shù)快速擴(kuò)張事件，從而導(dǎo)致了中性檢驗(yàn)的負(fù)值和錯(cuò)配分布的單峰曲線。

綜上所述，唐古特虎耳草具有較高的遺傳變異，且主要由居群內(nèi)遺傳變異所致，大量特有單倍型是近期輻射分化的結(jié)果。唐古特虎耳草的進(jìn)化歷史與西藏圓柏復(fù)合群(Juniperustibeticacomplex)[25]、西藏沙棘(Hippophaetibetana)[26]、菊葉紅景天(Rhodiolachrysanthemifolia)[27]相似，第四紀(jì)冰期時(shí)在高原臺(tái)面存在多個(gè)微型避難所。第四紀(jì)冰期氣候波動(dòng)使得原來連續(xù)的居群片段化，形成了相互隔離的居群，并發(fā)生異域快速分化，形成大量特有單倍型。唐古特虎耳草居群對(duì)第四紀(jì)冰期氣候波動(dòng)的反應(yīng)可能更多地表現(xiàn)為垂直海拔高度的遷移，而非大規(guī)模的水平遷移。此研究為青藏高原地區(qū)譜系地理學(xué)研究提供了一個(gè)新的案例，為了解第四紀(jì)冰期高山植物的避難所和現(xiàn)有分布格局的成因提供了參考依據(jù)。

參考文獻(xiàn):

[1]MYERS N, MITTERMEIER R A, MITTERMEIER C G,etal. Biodiversity hotspots for conservation priorities [J].Nature, 2000,403(6 772): 853-858.

[2]HOORN C, MOSBRUGGER V, MULCH A,etal. Biodiversity from mountain building [J].NatureGeoscience, 2013,6(3): 154.

[3]HUGHES C E. The tropical Andean plant diversity powerhouse [J].NewPhytologist, 2016,210(4): 1 152-1 154.

[5]WEN J, ZHANG J-Q, NIE Z-L,etal. Evolutionary diversifications of plants on the Qinghai-Tibetan Plateau [J].FrontiersinGenetics, 2014,5: 4.

[6]XING Y, REE R H. Uplift-driven diversification in the Hengduan Mountains, a temperate biodiversity hotspot [J].ProceedingsoftheNationalAcademyofSciencesoftheUnitedStatesofAmerica, 2017,114(17): E3444-E3451.

[7]EATON D A R, FENSTER C B, HEREFORD J,etal. Floral diversity and community structure inPedicularis(Orobanchaceae) [J].Ecology, 2012,93(8): S182-S194.

[8]ZHANG J Q, MENG S Y, WEN J,etal. Phylogenetic relationships and character evolution ofRhodiola(Crassulaceae) based on nuclear ribosomal ITS and plastidtrnL-FandpsbA-trnHsequences [J].SystematicBotany, 2014,39(2): 441-451.

[9]WANG Y J, SUSANNA A, VON RAAB-STRAUBE E,etal. Island-like radiation ofSaussurea(Asteraceae: Cardueae) trigged by uplifts of the Qinghai-Tibetan Plateau [J].BotanicalJournaloftheLinneanSociety, 2009,97(4): 893-903.

[10]GAO Q B, LI Y H, GORNALL R J,etal. Phylogeny and speciation inSaxifragasect.Ciliatae(Saxifragaceae): evidence frompsbA-trnH,trnL-Fand ITS sequences [J].Taxon, 2015,64(4): 703-713.

[11]EBERSBACH J, MUELLNER-RIEHL A N, MICHALAK I,etal. In and out of the Qinghai-Tibet Plateau: divergence time estimation and historical biogeography of the large arctic-alpine genusSaxifragaL [J].JournalofBiogeography, 2017,44(4): 900-910.

[12]LIU J Q, WANG Y J, WANG A L,etal. Radiation and diversification within theLigularia-Cremanthodium-Paraseneciocomplex (Asteraceae) trigged by uplift of the Qinghai-Tibetan Plateau [J].MolecularPhylogeneticsandEvolution, 2006,38(1): 31-49.

[13]ZHENG B, XU Q, SHEN Y. The relationship between climate change and Quaternary glacial cycles on the Qinghai-Tibetan Plateau: review and speculation [J].QuaternaryInternational, 2002,97-98(1): 93-101.

[14]OWEN L A, BENN D I. Equilibrium-line altitudes of the Last Glacial Maximum for the Himalaya and Tibet: an assessment and evaluation of results [J].QuaternaryInternational, 2005,138-139(3): 55-78.

[15]LEHMKUHL F, OWEN L A. Late Quaternary glaciation of Tibet and the bordering mountains: a review [J].Boreas, 2005,34(2): 87-100.

[16]HEWITT G M. Post-glacial re-colonization of European biota [J].BiologicalJournaloftheLinneanSociety, 1999,68(1/2): 87-112.

[17]HEWITT G M. Speciation, hybrid zones and phylogeography-or seeing genes in space and time [J].MolecularEcology, 2001,10(3): 537-549.

[18]ZHANG D, FENGQUAN L, JIANMIN B. Eco-environmental effects of the Qinghai-Tibet Plateau uplift during the Quaternary in China [J].EnvironmentalGeology, 2000,39(12): 1 352-1 358.

[19]HEWITT G M. Genetic consequences of climatic oscillations in the Quaternary [J].PhilosophicalTransactionsoftheRoyalSocietyofLondonB:BiologicalSciences, 2004,359(1 442): 183-195.

[20]ZHANG Q, CHIANG T Y, GEORGE M,etal. Phylogeography of the Qinghai-Tibetan Plateau endemicJuniperusprzewalskii(Cupressaceae) inferred from chloroplast DNA sequence variation [J].MolecularEcology, 2005,14(11): 3 513-3 524.

[21]MENG L, YANG R, ABBOTT R J,etal. Mitochondrial and chloroplast phylogeography ofPiceacrassifoliaKom. (Pinaceae) in the Qinghai-Tibetan Plateau and adjacent highlands [J].MolecularEcology, 2007,16(19): 4 128-4 137.

[22]WANG L Y, IKEDA H, LIU T L,etal. Repeated range expansion and glacial endurance ofPotentillaglabra(Rosaceae) in the Qinghai-Tibetan Plateau [J].JournalofIntegrativePlantBiology, 2009,51(7): 698-706.

[23]WANG L, ABBOTT R J, ZHENG W,etal. History and evolution of alpine plants endemic to the Qinghai-Tibetan Plateau:Aconitumgymnandrum(Ranunculaceae) [J].MolecularEcology, 2009,18(4): 709-721.

[24]GAO Q, ZHANG D, DUAN Y,etal. Intraspecific divergences ofRhodiolaalsia(Crassulaceae) based on plastid DNA and internal transcribed spacer fragments [J].BotanicalJournaloftheLinneanSociety, 2012,168(2): 204-215.

[25]OPGENOORTH L, VENDRAMIN G G, MAO K,etal. Tree endurance on the Tibetan Plateau marks the world’s highest known tree line of the Last Glacial Maximum [J].NewPhytologist, 2010,185(1): 332-342.

[26]WANG H, QIONG L, SUN K,etal. Phylogeographic structure ofHippophaetibetana(Elaeagnaceae) highlights the highest microrefugia and the rapid uplift of the Qinghai-Tibetan Plateau [J].MolecularEcology, 2010,19(14): 2 964-2 979.

[27]GAO Q B, ZHANG F Q, XING R,etal. Phylogeographic study revealed microrefugia for an endemic species on the Qinghai-Tibetan Plateau:Rhodiolachrysanthemifolia(Crassulaceae) [J].PlantSystematicsandEvolution, 2016,302(9):1 179-1 193.

[28]WU C Y, RAVEN P H. Flora of China [M]. Beijing: Science Press; St. Louis: Missouri Botanical Garden Press, 2001,8: 208-344.

[29]ABBOTT R J, COMES H P. Evolution in the Arctic: a phylogeographic analysis of the circumarctic plant,Saxifragaoppositifolia(Purple saxifrage) [J].NewPhytologist, 2003,161(1): 211-224.

[30]DECHAINE E G, ANDERSON S A, MCNEW J M,etal. On the evolutionary and biogeographic history ofSaxifragasect.Trachyphyllum(Gaud.) Koch (Saxifragaceae Juss.) [J].PLoSOne, 2013,8(7): e69814.

[31]EBERSBACH J, SCHNITZLER J, FAVRE A,etal. Evolutionary radiations in the species-rich mountain genusSaxifragaL [J].BMCEvolutionaryBiology, 2017,17(1): 119.

[32]GAO Q B, LI Y, GENGJI Z M,etal. Population genetic differentiation and taxonomy of three closely related species ofSaxifraga(Saxifragaceae) from southern Tibet and the Hengduan Mountains [J].FrontiersinPlantScience, 2017,8: 1 325.

[33]DOYLE J J, DOYLE J L. A rapid DNA isolation procedure for small quantities of fresh leaf material [J].PhytochemicalBulletin,BotanicalSocietyofAmerica, 1987,19(1): 11-15.

[34]TABERLET P, GIELLY L, PAUTOU G,etal. Universal primers for amplification of three non-coding regions of chloroplast [J].PlantMolecularBiology, 1991,17(5): 1 105-1 109.

[35]SCARCELLI N, BARNAUD A, EISERHARDT W,etal. A set of 100 chloroplast DNA primer pairs to study population genetics and phylogeny in Monocotyledons [J].PLoSOne, 2011,6(5): e19954.

[36]KUMAR S, STECHER G, TAMURA K. MEGA 7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets [J].MolecularBiologyandEvolution, 2016,33(7): 1 870-1 874.

[37]LIBRADO P, ROZAS J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data [J].Bioinformatics, 2009,25(11): 1 451-1 452.

[38]PONS O, PETIT R J. Measuring and testing genetic differentiation with ordered versus unordered alleles [J].Genetics, 1996,144(3): 1 237-1 245.

[39]EXCOFFIER L, LISCHER H E L. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows [J].MolecularEcologyResources, 2010,10(3):564-567.

[40]WEIR B S, COCKERHAM C C. Estimating F-Statistics for the analysis of population structure [J].Evolution, 1984,38(6): 1 358-1 370.

[41]TAJIMA F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism [J].Genetics, 1989,123(3): 585-595.

[42]FU Y X, LI W H. Statistical tests of neutrality of mutations [J].Genetics, 1993,133(3): 693-709.

[43]BANDELT H J, FORSTER P, R?HL A. Median-joining networks for inferring intraspecific phylogenies [J].MolecularBiologyandEvolution, 1999,16(1): 37-48.

[44]POLZIN T, DANESHMAND S V. On Steiner trees and minimum spanning trees in hypergrahs [J].OperationsResearchLetters, 2003,31(1): 12-20.

[45]GOLDING G B. The detection of deleterious selection using ancestors inferred from a phylogenetic history [J].GeneticsResearch, 1987,49(1): 71-82.

[46]CRANDALL K A, TEMPLETON A R. Empirical tests of some predictions from coalescent theory with applications to intraspecific phylogeny reconstruction [J].Genetics, 1993,134(3): 959-969.

[47]WANG F Y, GONG X, HU C M,etal. Phylogeography of an alpine speciesPrimulasecundiflorainferred from the chloroplast DNA sequence variation [J].JournalofSystematicsandEvolution, 2008,46(1): 13-22.

[48]王曉雄,樂霽培,孫航,等.青藏高原高山流石灘特有植物綿參的譜系地理學(xué)研究 [J]. 植物分類與資源學(xué)報(bào), 2011,33(6): 605-614.

WANG X X, YUE J P, SUN H,etal. Phylogeographical study onEriophytonwallichii(Labiatae) from Alpine scree of Qinghai-Tibetan Plateau [J].PlantDiversity&Resources, 2011,33(6): 605-614.

[49]DYNESIUS M, JANSSON R. Evolutionary consequences of changes in species' geographical distributions driven by Milankovitch climate oscillations [J].ProceedingsoftheNationalAcademyofSciencesoftheUnitedStatesofAmerica, 2000,97(16): 9 115-9 120.

[50]于海彬, 張鐿鋰. 青藏高原及其周邊地區(qū)高山植物譜系地理學(xué)研究進(jìn)展 [J]. 西北植物學(xué)報(bào), 2013,33(6): 1 268-1 278.

YU H B, ZHANG Y L. Advances in phylogeography of alpine plants in the Tibetan Plateau and adjacent regions [J].ActaBotanicaBoreali-OccidentaliaSinica, 2013,33(6): 1 268-1 278.

[51]RAY N, CURRAT M, EXCOFFIER L. Intra-deme molecular diversity in spatially expanding populations [J].MolecularBiologyandEvolution, 2003,20(1): 76-86.

[52]EXCOFFIER L. Patterns of DNA sequence diversity and genetic structure after a range expansion: lessons from the infinite-island model [J].MolecularEcology, 2004,13(4): 853-864.

[53]HARPENDING H C, BATZER M A, GURVEN M,etal. Genetic traces of ancient demography [J].ProceedingsoftheNationalAcademyofSciencesoftheUnitedStatesofAmerica, 1998,95(4): 1 961-1 967.