荒漠沙蜥遺傳多樣性的地理變異

劉昌景, 趙 偉, 周 蓉, 劉廼發(fā)

（蘭州大學(xué) 生命科學(xué)學(xué)院動(dòng)物學(xué)研究所, 甘肅 蘭州 730000）

荒漠沙蜥遺傳多樣性的地理變異

劉昌景, 趙 偉, 周 蓉, 劉廼發(fā)*

（蘭州大學(xué) 生命科學(xué)學(xué)院動(dòng)物學(xué)研究所,甘肅 蘭州730000）

荒漠沙蜥(Phrynocephalus przewalskii)是一種進(jìn)化地位至今存在爭(zhēng)議的沙漠蜥蜴, 該文基于12個(gè)不同地理種群119號(hào)荒漠沙蜥樣本線粒體基因組ND2序列, 在分析其系統(tǒng)進(jìn)化歷史的同時(shí), 研究了不同環(huán)境因素對(duì)其遺傳多樣性的影響。結(jié)果發(fā)現(xiàn), 12個(gè)地理種群聚類成3個(gè)明顯的分枝, 變色沙蜥(P. versicolor)和A枝處于與B、C兩枝并列的位置。所調(diào)查種群的遺傳多樣性存在顯著的地理變異, 分析顯示緯度、降雨和降雨變異對(duì)種群的遺傳多樣性都存在影響：隨緯度升高、降雨增加, 或降雨變異增加, 種群遺傳多樣性均有下降趨勢(shì)。推測(cè)緯度可能通過(guò)影響降雨和降雨變異對(duì)遺傳多樣性產(chǎn)生作用。

荒漠沙蜥; 系統(tǒng)進(jìn)化分析; 遺傳多樣性; 環(huán)境因素

生物在進(jìn)化過(guò)程中, 其表型特征受到遺傳基因和環(huán)境因素的共同作用。在自然選擇作用下, 個(gè)體遺傳型是表型進(jìn)化的基礎(chǔ), 而環(huán)境因素是表型分化的原因。研究表明，盡管表型進(jìn)化與分子進(jìn)化的速率不同, 但自然選擇對(duì)某特定表型的偏好經(jīng)常伴隨著該表型對(duì)應(yīng)的基因型發(fā)生頻率上的改變(Bruna et al, 1996)。許多生態(tài)遺傳學(xué)上的經(jīng)典研究, 如白樺尺蛾(Biston betularia)的“工業(yè)黑化”與其表型的選擇逆轉(zhuǎn)；深紅虎蛾(Panaxia dominula)的平衡多態(tài)性和生態(tài)因素對(duì)蝸牛(Fruticicolidaespp.)顏色變化的影響(Merrell, 1981), 它們已經(jīng)為環(huán)境因素可以通過(guò)影響個(gè)體表型進(jìn)而影響基因型頻率的觀點(diǎn)提供了經(jīng)驗(yàn)上的證據(jù)。在過(guò)去的20年里, 分子遺傳標(biāo)記的快速發(fā)展, 使得直接對(duì)基因型的地理變異進(jìn)行研究成為可能。同一物種生活在不同環(huán)境下的種群, 遺傳多樣性存在明顯地理變異, 而且這種變異與地理緯度、氣溫及降水等環(huán)境因子存在緊密聯(lián)系(Bernatchez & Wilson, 1998; Huang et al, 2005, 2007; Jin & Liu, 2008; Randi & Alkon, 1994)。對(duì)同一物種不同地理種群進(jìn)行遺傳多樣性變異的研究, 能夠?yàn)樯锒鄻有赃M(jìn)化理論提供重要的驗(yàn)證。

對(duì)特定物種進(jìn)行種群遺傳學(xué)研究過(guò)程中, 往往受到近緣種、伴生種或隱種的干擾。由于源自同一祖先的不同物種, 其表型差異更可能起源于祖先之間的表型差異, 而不是環(huán)境因素作用的結(jié)果(Harvey et al, 1995)。因此, 明確所研究物種的系統(tǒng)分類地位是種群遺傳學(xué)研究需要首先解決的難題。目前我們所熟知的物種分類地位的確定大多是基于傳統(tǒng)分類方法, 這種分類方法完全以物種的表型特征作為依據(jù)。由于表型特征在進(jìn)化過(guò)程中受到了自然選擇的洗禮, 不能完全反應(yīng)物種間的親緣度。因此, 以中性分子標(biāo)記為分類依據(jù)的分子系統(tǒng)學(xué)越來(lái)越受到分類學(xué)家的青睞(Hey, 2001)。

荒漠沙蜥(Phrynocephalus przewalskii)隸屬于蜥蜴目鬣蜥科沙蜥屬。它是棲息于荒漠或半荒漠地區(qū)的一種蜥蜴, 廣泛分布于我國(guó)北部及蒙古的荒漠地區(qū)。由于其表型特征變化多樣, 明顯表現(xiàn)出適應(yīng)各種荒漠生活的特征。因此, 它為研究同一物種不同種群遺傳多樣性的地理變異提供了很好的材料。然而, 當(dāng)前荒漠沙蜥的分類地位還存在廣泛爭(zhēng)議,主要表現(xiàn)為與變色沙蜥(Phrynocephalus versicolor)和草原沙蜥(Phrynocephalus frontalis)的親緣關(guān)系不甚明確。這3種蜥蜴在分布區(qū)上有較大重疊。最早, Zhao(1979)和Zhao et al(1999)根據(jù)形態(tài)學(xué)特征將這3種蜥蜴劃分為3個(gè)獨(dú)立種。Wang & Wang (1993)基于形態(tài)特征對(duì)分布在騰格里沙漠的3種沙蜥分析認(rèn)為：同一物種不同地理種群可能不是連續(xù)的, 盡管它們形態(tài)上還是相似的, 這種分布格局起因于青藏高原隆升引發(fā)的干旱化, 而與黃河形成無(wú)關(guān)。然而, 近年來(lái)基于線粒體分子標(biāo)記的系統(tǒng)學(xué)研究推翻了這種分類結(jié)果。例如, Pang et al (2003)基于線粒體(12S rRNA, 16S rRNA, 細(xì)胞色素b和ND4-tRNALEU)基因的分析結(jié)果表明，草原沙蜥、變色沙蜥和荒漠沙蜥聚類成同一個(gè)枝系, 其分類地位不甚明確。Wang & Fu(2004)分析認(rèn)為：3種沙蜥并未形成相互獨(dú)立的種, 而是依據(jù)地理分布形成東部、西北和西南3枝。分布在騰格里沙漠、鄂爾多斯高原及更東的種群可以命名為荒漠沙蜥(P. przewalskii), 而分布在西北地區(qū)的種群是變色沙蜥(P. versicolor), 西南部的可能是一個(gè)新種, 稱為“新荒漠沙蜥”(P. frontalis)。Gozdzik & Fu (2009)對(duì)“新荒漠沙蜥”的進(jìn)一步研究發(fā)現(xiàn), 線粒體數(shù)據(jù)支持將該種分成兩個(gè)獨(dú)立枝, 而形態(tài)和核基因的數(shù)據(jù)支持將其列為一個(gè)獨(dú)立枝。這些研究表明, 荒漠沙蜥和變色沙蜥及草原沙蜥的分類關(guān)系還沒(méi)有達(dá)成一致的結(jié)論。在研究荒漠沙蜥不同地理種群遺傳變異時(shí)需要排除另外兩種蜥蜴的干擾。

在本研究中, 我們首先應(yīng)用分子系統(tǒng)學(xué)方法對(duì)3種沙漠蜥蜴——荒漠沙蜥、草原沙蜥和變色沙蜥的系統(tǒng)進(jìn)化關(guān)系進(jìn)行了分析, 并與已發(fā)表的資料進(jìn)行了比對(duì)以檢驗(yàn)一致性。接下來(lái)采集荒漠沙蜥采樣地的環(huán)境因子, 包括：經(jīng)緯度、海拔、年均溫、年均降雨量、年均溫變異系數(shù)和年降雨量, 探討環(huán)境因子對(duì)不同荒漠沙蜥種群遺傳多樣性的影響。

1 材料和方法

1.1 研究地點(diǎn)和野外采樣

為了避免采樣過(guò)程中受到變色沙蜥和草原沙蜥的干擾, 我們的采樣地點(diǎn)集中在荒漠沙蜥分布區(qū)域的東部, 包括巴丹吉林沙漠、騰格里沙漠, 以及河西走廊東段。在三個(gè)研究區(qū)域選取了12個(gè)不同的地理種群(圖1)。各地區(qū)詳細(xì)的環(huán)境數(shù)據(jù)見(jiàn)表1,一部分?jǐn)?shù)據(jù)為采樣時(shí)測(cè)量, 包括經(jīng)緯度、海拔; 其余從國(guó)家氣象數(shù)據(jù)共享服務(wù)網(wǎng)獲取, 包括過(guò)去30年間年均溫和年均降雨量。在這12個(gè)不同的地理種群共采集到119號(hào)樣本。樣本采集后測(cè)量身體參數(shù)后處死, 然后浸泡在酒精中, 常溫保存。

圖1 荒漠沙蜥取樣圖與系統(tǒng)發(fā)生分支圖Fig. 1 The sites of sampled populations of Phrynocephalus przewalskii

1.2 分子生物學(xué)實(shí)驗(yàn)

選取動(dòng)物肝臟作為實(shí)驗(yàn)材料, 每個(gè)樣本取100 mg?？侱NA按照傳統(tǒng)酚-氯仿抽提法提取(Sambrook & Russell, 2001), DNA用200 μL TE緩沖液溶解, 并用紫外分光光度計(jì)檢測(cè)終濃度。兩條線粒體基因ND2的PCR引物H5617b/L5002(Macey & Schulte, 1998)用來(lái)進(jìn)行擴(kuò)增實(shí)驗(yàn)。PCR反應(yīng)體系：50 μL反應(yīng)液中含10 mmol/L Tris-HCl, pH 8.3;50 mmol/L KCl, 1.5 mmol/L MgCl2; Taq酶1U; 四種dNTP各150 mmol/L; 兩引物各10 pmol/L; DNA模板約100 ng。PCR反應(yīng)程序：95 ℃預(yù)變性4 min, 35個(gè)循環(huán)(94 ℃變性50 s, 48 ℃退火45 s, 72 ℃延伸50 s), 最后72 ℃延伸10 min（PE9700型PCR儀)。1.0%瓊脂糖凝膠電泳檢測(cè)PCR產(chǎn)物。線粒體ND2基因序列由上海生工生物工程公司進(jìn)行雙向測(cè)序。每個(gè)樣本雙向測(cè)序結(jié)果用CONTIG軟件進(jìn)行拼接,序列的排列和剪切由Clustal X軟件完成(Thompson et al, 1997)。

表1 12個(gè)荒漠沙蜥地理種群所在地點(diǎn)氣候與環(huán)境因子Tab. 1 Environmental factors at twelve sample locations

1.3 遺傳多樣性分析

種群的遺傳多樣性, 包括單倍型多樣性(h)和核苷酸多樣性(π)用軟件DNASP(5.0版, Librado & Rozas, 2009)來(lái)完成。MrBayesian系統(tǒng)發(fā)生分析用MrBayes 3.1.2軟件(Huelsenbeck & Ronquist, 2001)進(jìn)行。用青海沙蜥(Phrynocephalus valangalii)、貴德沙蜥(Phrynocephalus valangalii)和變色沙蜥作為外群。用MrModeltest(Nylander, 2004)進(jìn)行模型檢測(cè),用MDIV計(jì)算分歧時(shí)間, 采用的進(jìn)化速率為0.0156 substitutions/site/Ma。用ARELEQUIN軟件(3.1版)進(jìn)行AMOVA分析。

1.4 統(tǒng)計(jì)分析

我們將年均溫和年降雨量的變異系數(shù)作為統(tǒng)計(jì)分析中的環(huán)境變量(表1)。 由于氣候因子的絕對(duì)值及其穩(wěn)定性可能同時(shí)與遺傳多樣性相關(guān)(Huang et al, 2005, 2007), 故采用偏相關(guān)分析法控制其中一項(xiàng)(如絕對(duì)值)分析遺傳多樣性與另一項(xiàng)(如變異性)的關(guān)系。統(tǒng)計(jì)分析過(guò)程用軟件SPSS(16.0版)來(lái)完成,所有提供的數(shù)據(jù)用Mean±SE來(lái)表示, 顯著性水平設(shè)置為P=0.05。在進(jìn)行統(tǒng)計(jì)分析前, 連續(xù)型變量首先進(jìn)行正態(tài)性和方差齊性檢驗(yàn), 對(duì)于不滿足正態(tài)性和方差齊性的數(shù)據(jù)用非參數(shù)檢驗(yàn)的統(tǒng)計(jì)方法。

2 結(jié) 果

2.1 系統(tǒng)發(fā)育

文中測(cè)定的119條序列GenBank登錄號(hào)為：JQ688119-JQ688237。在119號(hào)荒漠沙蜥線粒體ND2基因序列中總共鑒定出37個(gè)單倍型。每個(gè)地理種群平均有9.8個(gè)(±3.9, 2～15,n=12)。安西種群?jiǎn)伪缎蛿?shù)目最少(2個(gè)), 山丹種群?jiǎn)伪缎蛿?shù)據(jù)最多(15個(gè))?；?7個(gè)線粒體單倍型所做的Bayesian系統(tǒng)樹(shù)顯示, 12個(gè)地理種群分成三個(gè)支持率非常高的枝系(A、B和C, 圖2)。其中A枝處于巴丹吉林沙漠的西邊, 主要包括額濟(jì)納旗、安西和山丹3個(gè)種群; B枝處于烏蘭布和沙漠周邊, 包括吉蘭泰、平羅和烏拉特后旗種群; C枝處于采樣的南部, 包括阿拉善左旗和右旗、民勤、大靖、甘塘和欒井。AMOVA的結(jié)果(表2)顯示, 大部分單倍型變異都存在于與這3個(gè)地理分支相對(duì)應(yīng)的3個(gè)地理群組之間, 而且地理群組間的遺傳分化極顯著。

Midv計(jì)算分析表明，A枝與B、C兩枝的分歧時(shí)間約為0.797 Ma前, B枝與C枝的分歧時(shí)間為0.618 Ma前。一個(gè)很明顯的結(jié)果是變色沙蜥和荒漠沙蜥的A枝聚類到一起, 它們處于與B、C兩枝并列的位置上, 說(shuō)明變色沙蜥和荒漠沙蜥在系統(tǒng)進(jìn)化上并沒(méi)有完全達(dá)到獨(dú)立種的水平。另外, 我們采集到的荒漠沙蜥樣本和兩個(gè)外群都處于并列位置,表明我們采集到的樣本可以視為同一個(gè)種, 滿足了對(duì)同一物種不同種群的遺傳多樣性進(jìn)行研究的前提。

圖2 根據(jù)37個(gè)單倍型616個(gè)堿基的mtDNA序列構(gòu)建的荒漠沙蜥貝葉斯進(jìn)化樹(shù)Fig. 2 Bayesian phylogenetic tree for 37 haplotypes using a 616 base pair fragment of the mitochondrial DNA of Phrynocephalus przewalskii

2.2 遺傳多樣性的地理變異

荒漠沙蜥各地理種群線粒體DNA單倍型數(shù)目、單倍型多樣性和核苷酸多樣性見(jiàn)表3。每個(gè)采樣點(diǎn)樣本量和單倍型數(shù)目具有顯著的正相關(guān)性(R2=0.725,n=12,P=0.008), 表明各地理種群?jiǎn)伪缎蛿?shù)目明顯受到樣本量的影響; 而單倍型多樣性(R2=0.0.26,n=12,P=0.437)和核苷酸多樣性(R2=0.08,n=12,P=0.824)不受樣本量的影響, 說(shuō)明可以用單倍型多樣性和核苷酸多樣性來(lái)表示該地理種群的遺傳多樣性高低。

表2 ND2部分序列分子變異分析(AMOVA)結(jié)果Tab. 2 Analysis of molecular variance (AMOVA) of Phrynocephalus przewalskii based on ND2 partial sequences

表3 荒漠沙蜥不同種群的單倍型多樣性和核苷酸多樣性Tab. 3 Haplotype and nucleotide diversities of Phrynocephalus przewalskii from different populations

Pearson相關(guān)分析結(jié)果顯示, 在經(jīng)緯度和海拔與種群的核苷酸多樣性均無(wú)顯著相關(guān)性; 僅緯度與種群的單倍型多樣性有顯著負(fù)相關(guān)(R2=0.630,n=12,P=0.038)。偏相關(guān)分析表明年均溫、溫度變異或年均降雨、降雨變異均對(duì)單倍型多樣性或核苷酸多樣性無(wú)顯著影響；但遺傳多樣性有隨降雨或降雨變異性增加而下降的趨勢(shì)(表4)。

表4 荒漠沙蜥遺傳多樣性和環(huán)境因素的關(guān)系Tab. 4 The relationship coefficients between genetic diversity and environmental factors

當(dāng)僅考慮“新荒漠沙蜥”種內(nèi), 即B和C支時(shí),偏相關(guān)表明單倍型多樣性與降雨變異系數(shù)有顯著負(fù)相關(guān)關(guān)系(R2=0.743,P=0.035); 而無(wú)論溫度或溫度變異仍與遺傳多樣性無(wú)關(guān)。緯度依然顯示與單倍型多樣性負(fù)相關(guān), 并且?guī)缀踹_(dá)到顯著水平(R2=0.658,P=0.054)。

3 討 論

一般認(rèn)為種群遺傳多樣性是在不斷的適應(yīng)環(huán)境和進(jìn)化中形成的, 是研究生物與環(huán)境相適應(yīng)的重要內(nèi)容(Merrell, 1981)。對(duì)適應(yīng)于局部環(huán)境的個(gè)體的自然選擇就會(huì)導(dǎo)致地理種群間的遺傳分化。氣候因子是導(dǎo)致生物進(jìn)化的重要因子, 其中溫度和降水是影響陸生脊椎動(dòng)物進(jìn)化的重要生態(tài)因子。Nevo(1981)發(fā)現(xiàn), 斑鬣蜥(Agama stellio)的遺傳多樣性與降水量相關(guān), 其次是溫度; Jin & Liu(2008)則發(fā)現(xiàn), 青海沙蜥的遺傳多樣性與降水量負(fù)相關(guān), 與溫度正相關(guān)。蜥蜴作為外溫動(dòng)物, 溫度和降雨對(duì)其生存有重要影響是顯而易見(jiàn)的, 但是也有研究表明溫度會(huì)對(duì)內(nèi)溫動(dòng)物的遺傳多樣性有重要影響(如石雞, Huang et al, 2005)?；哪瞅岬倪z傳多樣性存在顯著的地理變異, 但多樣性與溫度和溫度變異均無(wú)顯著相關(guān), 且相關(guān)系數(shù)較低。與分布在高海拔的青海沙蜥不同, 荒漠沙蜥分布在低海拔荒漠中, 溫度往往高于其生理極限, 故無(wú)論是年均溫度還是溫度變異性對(duì)荒漠沙蜥的影響應(yīng)該不會(huì)特別明顯。

在干旱、半干旱區(qū), 降雨往往成為驅(qū)動(dòng)和限制物種生命活動(dòng)的重要因子(Noy-Meir, 1973)。Jin et al (2008) 對(duì)青海沙蜥的研究表明, 青海沙蜥遺傳多樣性跟降水顯著負(fù)相關(guān), 并認(rèn)為只有少數(shù)單倍型被高降水所選擇。本文對(duì)荒漠沙蜥的研究也發(fā)現(xiàn), 荒漠沙蜥遺傳多樣性與降雨量之間的負(fù)相關(guān)性, 但并未達(dá)到顯著水平。這表明適于干旱環(huán)境的沙蜥會(huì)被高降雨量所限制, 導(dǎo)致其遺傳多樣性下降。

除環(huán)境因子的絕對(duì)值外, 其穩(wěn)定性亦會(huì)對(duì)多樣性產(chǎn)生重要影響(Huang et al, 2005, 2007)。用于解釋群落多樣性的氣候穩(wěn)定性學(xué)說(shuō)同樣可以用于解釋種群的遺傳多樣性。在穩(wěn)定環(huán)境中選擇相對(duì)松弛,會(huì)有更多的遺傳型被保存下來(lái); 在波動(dòng)環(huán)境中僅有那些適應(yīng)環(huán)境的遺傳型才會(huì)被保。由此可以預(yù)計(jì)物種的遺傳多樣性會(huì)隨環(huán)境的穩(wěn)定性增加而上升。Jin & Liu (2008) 發(fā)現(xiàn)青海沙蜥的遺傳多樣性隨溫度變異系數(shù)和降雨變異系數(shù)的增加有下降趨勢(shì)。本文亦發(fā)現(xiàn)雖然荒漠沙蜥的遺傳多樣性與降雨變異系數(shù)之間的相關(guān)未達(dá)到顯著水平, 但仍表明隨降雨變異系數(shù)的增加遺傳多樣性有下降趨勢(shì)。

地理因素可能會(huì)通過(guò)影響環(huán)境條件, 如溫度和降雨對(duì)生物遺傳多樣性產(chǎn)生影響。青藏高原地區(qū)海拔對(duì)溫度的影響要高于緯度影響, 而經(jīng)度變異往往伴隨著降雨量的改變, 故Jin & Liu (2008)認(rèn)為對(duì)青海沙蜥遺傳多樣性有影響的海拔與經(jīng)度可能是通過(guò)影響氣溫或降雨來(lái)起作用的。本文發(fā)現(xiàn), 荒漠沙蜥遺傳多樣性隨緯度增加有下降的趨勢(shì)。Pearson相關(guān)分析表明，緯度對(duì)降雨量和降雨變異系數(shù)有顯著影響而對(duì)溫度無(wú)顯著影響, 具體表現(xiàn)為隨緯度升高降雨量減少而降雨量變異系數(shù)增加。由此可見(jiàn),緯度對(duì)荒漠沙蜥遺傳多樣性的影響可能是通過(guò)影響降雨而產(chǎn)生作用的。

Pang et al (2003)對(duì)沙蜥屬物種線粒體數(shù)據(jù)的分析表明, 變色沙蜥、草原沙蜥和荒漠沙蜥聚成支持率較高的一支。Wang & Fu (2004)針對(duì)該三種進(jìn)行的研究發(fā)現(xiàn), 它們并未形成相互獨(dú)立的單系, 而是依據(jù)地理分布形成東部、西北和西南三支; 并提出東部支(騰格里沙漠、毛烏素沙漠及其東部)應(yīng)保留種名P. przewalskii(包含原P. frontalis), 西北支為P. versicolor, 西南部可能是一個(gè)新種。Gozdzik & Fu(2009)也通過(guò)分子和形態(tài)研究證明分布在內(nèi)蒙東部和蒙古的變色沙蜥和荒漠沙蜥應(yīng)為一個(gè)種。本文所構(gòu)建的系統(tǒng)樹(shù)表明, 變色沙蜥(P. versicolor)與荒漠沙蜥(P. przewalskii)并沒(méi)有形成相互獨(dú)立的單系,說(shuō)明它們很可能是一個(gè)種。系統(tǒng)樹(shù)中所有單倍型聚類成了三個(gè)支持率較高的地理支, 并且AMOVA的結(jié)果(表2)顯示, 大部分遺傳變異都存在于與這三個(gè)地理分支相對(duì)應(yīng)的三個(gè)地理群組之間, 而且地理群組間的遺傳分化極顯著, 結(jié)果支持Wang & Fu(2004)的分類結(jié)論。與Gozdzik & Fu(2009)不同,我們的研究表明，3個(gè)分支的分歧是在最近1 Ma內(nèi)發(fā)生的, 而非2～5 Ma或更久, 這可能與采取的計(jì)算方式不同有關(guān)。

Bernatchez L, Wilson CC. 1998. Comparative phylogeography of Nearctic and Palearctic fishes[J]. Mol Ecol,7(4): 431-452.

Bruna EM, Fisher RN, Case TJ. 1996. Morphological and genetic evolution appear decoupled in Pacific skinks (Squamata: Scincidae: Emoia)[J]. Proc Biol Sci,263(1371): 681-688.

Gozdzik A, Fu JZ. 2009. Are toad-headed lizards Phrynocephalus przewalskii and P. frontalis (family agamidae) the same species? Defining species boundaries with morphological and molecular data[J]. Russian J Herpetol,16(2): 107-118.

Harvey PH, Read AF, Nee S. 1995. Why ecologists need to be phylogenetically challenged[J]. J Ecol,83(3): 535-536.

Hey J. 2001. The mind of the species problem[J]. Trends Ecol Evol,16(7): 326-329.

Huang ZH, Liu NF, Luo SX, Long J, Xiao YA. 2007. Ecological genetics of rusty-necklaced partridge (Alectoris magna): Environmental factors and population genetic variability correlations[J]. Korean J Genet,29(2): 115-120.

Huang ZH, Liu NF, Zhou T, Ju B. 2005. Effects of environmental factors on the population genetic structure in chukar partridge (Alectoris chukar)[J]. J Arid Environ,62(3): 427-434.

Huelsenbeck JP, Ronquist F. 2001. MRBAYES: Bayesian inference of phylogenetic trees[J]. Bioinformatics,17(8): 754-755.

Jin YT, Brown RP, Liu NF. 2008. Cladogenesis and phylogeography of the lizard Phrynocephalus vlangalii (Agamidae) on the Tibetan plateau[J]. Mol Ecol,17(8): 1971-1982.

Jin YT, Liu NF. 2008. Ecological genetics of Phrynocephalus vlangalii on the North Tibetan (Qinghai) Plateau: Correlation between environmental factors and population genetic variability[J]. Biochem Genet,46(9-10): 598-604.

Librado P, Rozas J. 2009. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data[J]. Bioinformatics,25(11): 1451-1452.

Macey JR, Schulte JA. 1998. Phylogenetic relationships of toads in the Bufo bufo species group from the eastern escarpment of the Tibetan Plateau: A case of vicariance and dispersal[J]. Mol Phylogenet Evol,9(1):80-87.

Merrell DJ. 1981. Ecological Genetics[M]. Minneapolis: University of Minnesota Press.

Nevo E. 1981. Genetic variation and climatic selection in the Lizard Agama stellio in Israel and Sinai[J]. Theor Appl Genet,60(6): 369-380.

Noy-Meir I. 1973. Desert ecosystems: environment and producers[J]. Annu Rev Ecol Syst,4(1): 25-51.

Nylander JAA. 2004. MrModeltest Version 2. 0. Program Distributed by the Author[M]. Uppsala: Evolutionary Biology Centre, Uppsala University

Pang JF, Wang YZ, Zhong Y, Hoelzel AR, Papenfuss TJ, Zeng XM, Ananjeva NB, Zhang YP. 2003. A phylogeny of Chinese species in the genus Phrynocephalus (Agamidae) inferred from mitochondrial DNA sequences[J]. Mol Phylogenet Evol,27(3): 398-409.

Randi E, Alkon PU. 1994. Genetic structure of Chukar (Alectoris chukar) populations in Israel[J]. Auk,111(2): 416-426.

Sambrook J, Russell DW. 2001. Molecular Cloning: A Laboratory Manual[M]. New York: Cold Spring Harbor Laboratory Press.

Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. 1997. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools[J]. Nucleic Acids Res,25(24): 4876-4882.

Wang YZ, Fu JZ. 2004. Cladogenesis and vicariance patterns in the toad-headed lizard Phrynocephalus versicolor species complex[J]. Copeia,2004(2): 199-206.

Wang YZ, Wang HZ. 1993. Geographic variation and diversity in three species of Phrynocephalus in the Tengger Desert, western China[J]. Asia Herpetol Res,5: 65-73.

Zhao EM, Zhao KT, Zhou KY. 1999. Fauna Sinica: Reptilia[M]. Beijing: Science Press, 237-258. [趙爾宓,趙肯堂,周開(kāi)亞. 1999.中國(guó)動(dòng)物志:爬行綱[M].北京:科學(xué)出版社, 237-258.]

Zhao KT. 1979. A survey of the classilication and distribution of the toad-headed Agamids (Pnrynocephalus) in China[J]. Acta Sci Nat Univ Neimongol, (2): 111-121.[趙肯堂. 1979.中國(guó)沙蜥屬的分類和分布研究[J].內(nèi)蒙古大學(xué)學(xué)報(bào):自然科學(xué)版, (2): 111-121.]

Phylogeny and genetic diversity ofPhrynocephalus przewalskii

LIU Chang-Jing, ZHAO Wei, ZHOU Rong, LIU Nai-Fa*

(Institute of Zoology, School of Life Science, Lanzhou University, Lanzhou Gansu 730000)

Phrynocephalus przewalskiiis one kind of lizards inhabiting inclusively in the desert that has controversial viewpoints on its phylogeny. Based on mitochondrial ND2 gene of 119 samples from 12 geographic populations, we analyzed the effects of environmental factors on the variation of genetic diversity, as well as its relationship toP. versicolor. The results showed that these populations clustered into three major lineages, withP. versicolorembedded within one lineage. The twelve populations had great genetic diversity variation, which was tightly linked with local altitude, annual precipitation, and variation of annual precipitation. High latitudes, increased annual precipitation and great variation in annual precipitations may all have resulted in the decrease of genetic diversity. It thus assumed that altitude can change the genetic diversity of different geographic populations ofP. przewalskiiresulting from the effects of different local annual precipitation.

Phrynocephalus przewalskii; Phylogenetic analysis; Genetic diversity; Environmental factor

Q959.62; Q349.1; Q951.3

A

0254-5853-(2012)02-0127-06

10.3724/SP.J.1141.2012.02127

2011-12-15；接受日期：2012-02-15

國(guó)家自然科學(xué)基金資助項(xiàng)目(41071031)

?通信作者(Corresponding author)，E-mail: naifaliu@sohu.com