王玉濤, 潘建飛,黃翠翠,胡 平,韓建林,,*

1 喀什大學(xué),生命與地理科學(xué)學(xué)院,喀什 844000 2 喀什大學(xué)新疆維吾爾自治區(qū)葉爾羌綠洲生態(tài)與生物資源研究高校重點(diǎn)實(shí)驗(yàn)室,喀什 844000 3 中國(guó)農(nóng)業(yè)科學(xué)院-國(guó)際家畜研究所畜禽牧草遺傳資源聯(lián)合實(shí)驗(yàn)室,北京 100073

雪雞屬鳥綱(Aves),雞形目(Galliformes),雉科(Phasianidae),雪雞屬(Tetraogallus),包括高加索雪雞(T.caucasicus)、里海雪雞(T.caspius)、阿爾泰雪雞(T.altaicus)、喜馬拉雅雪雞(T.himalayensis)和藏雪雞(T.tibetanus)5個(gè)種[1-2],后3種在中國(guó)境內(nèi)均有分布,其中喜馬拉雅雪雞又稱暗腹雪雞,藏雪雞又稱淡腹雪雞。喜馬拉雅雪雞是雪雞屬中體形最大的種[3],在國(guó)內(nèi)主要分布在昆侖山西部、祁連山脈和喜馬拉雅山西部等地,按照分布格局不同又分為有4個(gè)亞種：即指明亞種(T.h.himalayensis)分布于我國(guó)新疆西部；疆南亞種(T.h.grombszewskii)主要分布于昆侖山西部向東到新疆于田附近；青海亞種(T.h.koslowi)分布于自阿爾金山西部向東經(jīng)祁連山至青海湖南部地區(qū)；天山亞種(T.h.sewerzowi)分布于新疆天山山脈[3]。在國(guó)外,阿富汗、土耳其至尼泊爾、巴基斯坦和哈薩克斯坦均有分布[4-6]。昆侖山-帕米爾高原是喜馬拉雅雪雞的分布中心,然而由于全球氣候變化、人類活動(dòng)干擾、天敵侵害和非法狩獵,喜馬拉雅雪雞賴以生存環(huán)境不斷惡化,群體數(shù)量急劇下降。喜馬拉雅雪雞現(xiàn)為國(guó)家Ⅱ級(jí)重點(diǎn)保護(hù)動(dòng)物,已被IUCN紅色目錄列為低危(LC)物種[7],中國(guó)脊椎動(dòng)物紅色名錄列為近危(NV)[8]。

雪雞既具有較高的藥用和觀賞價(jià)值,又具有一定的生態(tài)價(jià)值。關(guān)于喜馬拉雅雪雞的研究主要集中在野外生境、巢址選擇、食性和繁殖習(xí)性,圈養(yǎng)馴化和繁育以及起源進(jìn)化方面[9-29]。在起源進(jìn)化方面,劉強(qiáng)等利用雪雞mtDNACytb基因全長(zhǎng)序列(1143bp)與雉科中其他19個(gè)屬的Cytb基因進(jìn)行同源性比較,結(jié)果表明雪雞屬在雉科中與鵪鶉屬親緣關(guān)系最近,其次是石雞屬,雪雞屬與鵪鶉屬的分歧時(shí)間為1.5百萬(wàn)—2.0百萬(wàn)年,與石雞屬分歧時(shí)間為1.7百萬(wàn)—2.4百萬(wàn)年[30]。Ruan等利用mtDNACytb基因部分序列(535bp)分析了6個(gè)喜馬拉雅雪雞群體,表明新疆天山地區(qū)雪雞群體遺傳多樣性豐富,且群體間存在基因交流[31-32]。王小立通過(guò)對(duì)喜馬拉雅雪雞的線粒體DNA控制區(qū)(mtDNAD-loop)序列和微衛(wèi)星進(jìn)行標(biāo)記檢測(cè),研究了我國(guó)喜馬拉雅雪雞11個(gè)地理種群的分子系統(tǒng)地理學(xué),結(jié)果表明11個(gè)地理種群分為4個(gè)進(jìn)化分支,且近期經(jīng)歷過(guò)種群分化、基因交流受到限制并且在冰期發(fā)生過(guò)擴(kuò)張[33-34]。但以上研究中,涉及東帕米爾高原境內(nèi)的喜馬拉雅雪雞樣本數(shù)量較少或不足,不足以代表整個(gè)地理種群。因此,關(guān)于雪雞地理學(xué)及遺傳多樣性研究,急需擴(kuò)大樣本數(shù)量以及地理覆蓋度。

帕米爾高原是地球上兩條巨大山帶 (阿爾卑斯-喜馬拉雅山帶和帕米爾-楚科奇山帶)的山結(jié),也是亞洲大陸南部和中部地區(qū)主要山脈的匯集處,這里群山起伏、連綿逶迤、雪峰群立、聳入云天,因此人口稀少,人類活動(dòng)較少,是野生動(dòng)物理想的棲息地。但隨著社會(huì)經(jīng)濟(jì)的飛速發(fā)展,人類活動(dòng)加劇,帕米爾高原植被退化嚴(yán)重,荒漠化進(jìn)程加速,野生動(dòng)物棲息地生境破碎化程度嚴(yán)重[35],邊緣效應(yīng)加劇,物種內(nèi)基因交流受阻,群體遺傳多樣性降低,導(dǎo)致該區(qū)域面臨生物多樣性喪失、生態(tài)系統(tǒng)服務(wù)低下和貧困的惡性循環(huán)[36],因此開展區(qū)域內(nèi)生物多樣性評(píng)估具有重要的意義。關(guān)于東帕米爾高原喜馬拉雅雪雞群體遺傳多樣性評(píng)估及相關(guān)研究尚未見報(bào)道,因此為闡明東帕米爾高原喜馬拉雅雪雞遺傳多樣性水平和系統(tǒng)進(jìn)化地位,本文選擇mtDNAD-loop序列信息,利用生物信息學(xué)軟件分析雪雞mtDNAD-loop序列特征和遺傳多樣性水平；選擇GenBank中已提交的雪雞序列,利用最大似然法構(gòu)建系統(tǒng)進(jìn)化樹以及中介網(wǎng)絡(luò)構(gòu)建系統(tǒng)發(fā)育關(guān)系,這可為喜馬拉雅雪雞遺傳資源保護(hù)和利用提供理論依據(jù)。

1 材料與方法

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

喜馬拉雅雪雞樣本是從2008年開始在新疆帕米爾高原的塔什庫(kù)爾干野生動(dòng)物自然保護(hù)區(qū)(紅其拉甫、達(dá)布達(dá)爾、恰拉奇古等區(qū)域)、阿克陶縣和帕米爾雪雞繁育基地收集雪雞尸體的肌肉組織,合計(jì)25個(gè)。收集的肌肉組織放入75%酒精中-80℃保存,肌肉組織樣本保存在喀什大學(xué)葉爾羌綠洲生態(tài)與生物資源研究高校重點(diǎn)實(shí)驗(yàn)室。

1.2 基因組DNA提取

肌肉組織基因組DNA提取參照《分子克隆》(第三版)的提取方法,采用常規(guī)的酚·氯仿抽提法,并結(jié)合實(shí)驗(yàn)室條件進(jìn)行了優(yōu)化和改進(jìn)。

1.3 引物設(shè)計(jì)和PCR擴(kuò)增

利用Primer 5軟件,根據(jù)GenBank上公布的雪雞mtDNA序列(登錄號(hào)：NC.027279)設(shè)計(jì)特異性擴(kuò)增引物。正向引物L(fēng)15725x-15705：5′-CCTAGCAGCCTCAGTACTCATCCTAC-3′,反向引物H1530x-1498：5′-CCGCGGTGGCTGGCACAAGATTTACC-3′,測(cè)序引物L(fēng)15725-seq1：5′-CCTCTATTTTTCCCACAAAACTAGG-3′。 引物由北京三博遠(yuǎn)志生物技術(shù)有限責(zé)任公司合成。

PCR反應(yīng)體系：總體積25μL,其中TaqDNA酶(500U)0.5μL,dNTPs(10mmol/L)0.5μL,Buffer(2.5mmol/L)12.5μL,上下游引物(10pmol/μL)1μL,DNA模板(100ng/μL)1μL,ddH2O 8.5μL。

擴(kuò)增反應(yīng)條件：94℃預(yù)變性3min,35個(gè)循環(huán)(94℃變性30s, 63℃退火30s, 72℃72延伸1.5min),72℃延伸7min,4℃保存。PCR產(chǎn)物采用0.8%的瓊脂糖凝膠電泳檢測(cè),擴(kuò)增產(chǎn)物送北京三博遠(yuǎn)志生物技術(shù)有限責(zé)任公司測(cè)序。

1.4 序列分析

Chromas Version 1.45對(duì)喜馬拉雅雪雞mtDNAD-loop區(qū)序列進(jìn)行人工校對(duì)。用MEGA 5.0比對(duì)序列以及統(tǒng)計(jì)序列堿基組成,并構(gòu)建系統(tǒng)發(fā)育樹。DnaSPV5用于統(tǒng)計(jì)單倍型,計(jì)算單倍型和核苷酸多樣性指標(biāo),利用Adobe Illustrator軟件對(duì)圖片進(jìn)行處理。利用Network4.6.1.1軟件對(duì)本研究所得序列和GenBank下載的喜馬拉雅雪雞和藏雪雞的同源mtDNAD-loop序列進(jìn)行網(wǎng)絡(luò)關(guān)系及其單倍型分布頻率分析。喜馬拉雅雪雞和藏雪雞序列信息見表1。

表1 雪雞mtDNA D-loop序列信息表Table 1 mtDNA D-loop sequence information of Tetraogallus

2 結(jié)果與分析

2.1 PCR擴(kuò)增

圖1 mtDNA PCR擴(kuò)增產(chǎn)物Fig.1 Amplifications of mtDNA gene

對(duì)25份喜馬拉雅雪雞mtDNA目的片段進(jìn)行PCR擴(kuò)增,擴(kuò)增產(chǎn)物用0.8%的瓊脂糖凝膠電泳檢測(cè)(圖1)。檢測(cè)產(chǎn)物大小與目的片段長(zhǎng)度一致(2496bp),且條帶清晰,無(wú)雜帶或拖帶,無(wú)非特異擴(kuò)增產(chǎn)物,可以用于PCR產(chǎn)物直接測(cè)序。

2.2 mtDNA D-loop 序列分析

2.2.1 mtDNA D-loop核苷酸組成

對(duì)25個(gè)喜馬拉雅雪雞mtDNAD-loop區(qū)進(jìn)行序列測(cè)定,測(cè)序結(jié)果經(jīng)比較、矯正后獲得了1157bp的全長(zhǎng)序列。與參照序列NC027279進(jìn)行對(duì)比編輯進(jìn)行分析。結(jié)果表明：mtDNAD-loop中T、C、A、G 4種堿基平均比例分別為33.0%、26.2%、26.8%、14.0%,A+T的平均含量為59.8%,顯著高于G+C平均含量40.2%。因此雪雞mtDNAD-loop區(qū)富含A和T。進(jìn)而說(shuō)明堿基組成具有偏倚性,符合禽類線粒體DNA序列特征。

2.2.2 mtDNA D-loop核苷酸變異與單倍型多樣性分析

利用DnaSP5.10.01軟件對(duì)25個(gè)喜馬拉雅雪雞的mtDNAD-loop高變區(qū)進(jìn)行單核苷酸多態(tài)性分析,以Hap1為標(biāo)準(zhǔn),共檢測(cè)到64個(gè)多態(tài)位點(diǎn),占所測(cè)核苷酸的5.5%。其中有29個(gè)單一多態(tài)位點(diǎn),位于第32、165、199、290、295、299、304、349、392、393、397、400、404、407、409、452、478、504、524、544、566、571、581、585、598、646、697、723和725位；33個(gè)簡(jiǎn)約信息位點(diǎn),位于第146、166、168、173、189、205、208、217、218、220、222、223、231、236、237、242、243、245、256、289、383、554、573、576、578、606、609、714、740、780、826、838和1003位；240和953位存在兩個(gè)插入/缺失位點(diǎn)。本實(shí)驗(yàn)檢測(cè)到轉(zhuǎn)換、顛換、插入或缺失4種核苷酸變異類型,其中轉(zhuǎn)換42次,占變異位點(diǎn)數(shù)的65.63%,A/G之間轉(zhuǎn)換8次,T/C之間轉(zhuǎn)換34次；顛換20次,占變異為點(diǎn)數(shù)的31.25%；存在兩處插入或缺失,在240處插入或缺一個(gè)C,在953處插入或缺失一個(gè)T。

通過(guò)DnsSP5.10.01軟件對(duì)喜馬拉雅雪雞的25個(gè)個(gè)體mtDNAD-loop區(qū)序列進(jìn)行多態(tài)性分析,結(jié)果共檢測(cè)到64個(gè)多態(tài)位點(diǎn),存在23種單倍型,命名為Hap1—Hap23(表2),25個(gè)雪雞中有2個(gè)個(gè)體共享單倍型Hap11,2個(gè)個(gè)體共享單倍型Hap19,其余23個(gè)個(gè)體均為獨(dú)立單倍型。平均單倍型多樣度(Hd)為0.92±0.0001,平均核苷酸多樣度(π)為0.00958,平均核苷酸差異度(K)為11.067,從而顯示東帕米爾高原喜馬拉雅雪雞核苷酸多樣性較低,單倍型多樣性高,具有較為豐富的遺傳多樣性。

表2 25個(gè)喜馬拉雅雪雞mtDNA D-loop 的23個(gè)單倍型變異序列Table 2 Sequence variation of 23 haplotypes in 25 Tetraogallus himalayensis group mtDNA D-loop

2.2.3 mtDNA D-loop區(qū)序列(1157bp)堿基變異和核苷酸多樣性分布

利用DnaSP5.10.01軟件分析高變區(qū)序列變異的位置(圖2)。結(jié)果表明：喜馬拉雅雪雞mtDNAD-loop序列堿基變異的分布有很大的差異,在1—80bp之間較保守,在80—300bp之間變異較大,300—500bp相對(duì)保守,在500—700bp之間變異增大,700—1157bp相對(duì)比較保守,其中在230bp處的堿基變異最豐富。這一結(jié)果符合mtDNAD-loop的結(jié)構(gòu)組成特點(diǎn),5′末端為高變區(qū)I,突變率較高的中間區(qū)域?yàn)楸Ｊ貐^(qū)域II,3′末端為高變區(qū)III。從而驗(yàn)證了擴(kuò)增序列是來(lái)自mtDNAD-loop區(qū)域。對(duì)mtDNAD-loop序列核苷酸多樣性分布分析(圖3),結(jié)果表明：在堿基1—80bp之間保守,80—270bp之間核苷酸多樣度較高,在270—500bp之間變異較少,500—700bp之間略微有所增加,在750bp之后逐漸趨于保守之勢(shì),在230bp處核苷酸多樣度的值最大。核苷酸變異分布與核苷酸多樣性分布基本一致。對(duì)mtDNAD-loop區(qū)序列(1157bp)錯(cuò)配分布分析(圖4),表明東帕米爾高原喜馬拉雅雪雞經(jīng)歷過(guò)群體擴(kuò)張。

2.3 系統(tǒng)發(fā)育與網(wǎng)絡(luò)關(guān)系分析

2.3.1 系統(tǒng)發(fā)育分析

結(jié)合本研究獲得數(shù)據(jù),下載GenBank已提交的雪雞屬mtDNAD-loop同源序列(886bp),比對(duì)分析共獲得84個(gè)單倍型雪雞(表1),利用MEGA6.1軟件,采用Kimura-2模型,以鄰接法構(gòu)建其NJ系統(tǒng)發(fā)育樹(圖5),并對(duì)拓?fù)鋱D進(jìn)行自展檢驗(yàn),重復(fù)抽樣次數(shù)為1000,然后以2條喜馬拉雅雪雞(GQ343551、GQ343550)的mtDNAD-loop同源序列作為根序列。

圖2 mtDNA D-loop區(qū)序列堿基變異分布Fig.2 mtDNA D-loop base variation distribution

圖3 mtDNA D-loop序列核苷酸多樣度分布Fig.3 mtDNA D-loop nucleotide diversity distribution

圖4 mtDNA D-loop堿基錯(cuò)配分布 Fig.4 Mismatch distributes for all Tetraogallus himalayensis samples

圖5 84個(gè)雪雞mtDNA D-loop單倍型構(gòu)建的NJ系統(tǒng)發(fā)育進(jìn)化樹Fig.5 Evolution tree Tetraogallus constructs 84 mtDNA D-loop haplotypes in NJ tree

從圖5系統(tǒng)發(fā)育樹可以發(fā)現(xiàn),84個(gè)單倍型共分為兩大支。一支為喜馬拉雅雪雞,一支為藏雪雞。而喜馬拉雅雪雞支又分為3支,本實(shí)驗(yàn)確定的23個(gè)單倍型喜馬拉雅雪雞聚在2個(gè)進(jìn)化支上,其中19個(gè)單倍型與Wang等[35],Luzhang等[33]采集自塔什庫(kù)爾干縣樣本(H32,H33,H34,AFXI-83)聚一個(gè)進(jìn)化支上,可見來(lái)自東帕米爾高原喜馬拉雅雪雞獨(dú)享進(jìn)化支,但單倍型Hap2、Hap17、Hap5、Hap15聚在另一個(gè)進(jìn)化支上,進(jìn)而說(shuō)明該區(qū)域喜馬拉雅雪雞出現(xiàn)種群分化。

2.3.2 網(wǎng)絡(luò)關(guān)系分析

Network網(wǎng)絡(luò)分析法是系統(tǒng)進(jìn)化分析最常用的方法之一,主要用于闡述種內(nèi)水平的系統(tǒng)發(fā)育關(guān)系,能清晰地表明單倍型之間的關(guān)系、各單倍型在各群體中的分布頻率及其基因流的方向。利用Network4.6.1.1軟件對(duì)雪雞mtDNAD-loop84個(gè)單倍型的網(wǎng)絡(luò)關(guān)系進(jìn)行分析(圖6)。結(jié)果顯示網(wǎng)絡(luò)關(guān)系分析與其系統(tǒng)發(fā)育NJ樹的結(jié)果完全一致,所有的單倍型聚為2個(gè)簇群,分別為喜馬拉雅雪雞簇群和藏雪雞簇群。按照動(dòng)物地理分布分析,本研究的19個(gè)單倍型與前人研究的4個(gè)單倍型聚在以東帕米爾高原為中心的喜馬拉雅雪雞簇中,且形成單獨(dú)的一支,不與其他區(qū)域雪雞共享進(jìn)化支；其余4個(gè)單倍型與其他單倍型聚在以昆侖山、天山、阿爾金山和祁連山為分布區(qū)域的喜馬拉雅雪雞簇中。由此可見喜馬拉雅雪雞遺傳特征具有明顯的地理分布特征,并由此形成不同的亞種,說(shuō)明雪雞進(jìn)化具有明顯的地理起源特征。所有藏雪雞均聚在以青藏高原為中心的簇中。

圖6 84個(gè)單倍型網(wǎng)絡(luò)關(guān)系Fig.6 Network relationship of the 84 haplotypes

3 討論

3.1 遺傳多樣性

本研究喜馬拉雅雪雞25個(gè)樣本獲得1157bp長(zhǎng)的mtDNAD-loop區(qū)序列,與其他鳥類控制區(qū)長(zhǎng)度差異不大。喜馬拉雅雪雞堿基A+T的含量為59.8%,高于G+C的含量40.2%,證明雪雞mtDNAD-loop區(qū)序列存在堿基偏倚性,符合鳥類線粒體DNA堿基組成偏倚性的特點(diǎn)。這與脊椎動(dòng)物mtDNA的特點(diǎn)一樣,存在一條重鏈和一條輕鏈,就是由于這兩條鏈的A+T的含量差異造成的。共發(fā)現(xiàn)64個(gè)變異位點(diǎn),占所檢測(cè)核苷酸的5.53%,其中有29個(gè)單一多態(tài)位點(diǎn),33個(gè)簡(jiǎn)約信息位點(diǎn)和兩個(gè)插入/缺失位點(diǎn),變異較高。

單倍型多樣度(Hd)和核苷酸多樣度(π值)是衡量一個(gè)群體mtDNA 變異程度的兩個(gè)重要指標(biāo)。π值越大遺傳多樣性越高。本實(shí)驗(yàn)喜馬拉雅雪雞π值為0.00958,平均單倍型多樣度為0.92±0.0001,說(shuō)明東帕米爾喜馬拉雅雪雞核苷酸多樣性低,而單倍型多樣性較高,這與王小立對(duì)喜馬拉雅雪雞的11個(gè)地理種群67個(gè)樣品的mtDNAD-loop的平均核苷酸多樣性為0.0103±0.00523,單倍型多樣性為0.977±0.00690結(jié)果一致[34]。安蓓對(duì)西藏雪雞部分mtDNAD-loop控制區(qū)序列進(jìn)行研究,結(jié)果表明藏雪雞的核苷酸多樣性(π值)為0.0036[37-38],由此可見,喜馬拉雅雪雞的核苷酸多樣性高于藏雪雞。研究認(rèn)為降低物種遺傳多樣性的主要原因有：(1)種群間隔離:長(zhǎng)時(shí)間缺乏基因交流；(2)棲息地破壞:自然種群棲息地片段化和嚴(yán)重近親繁殖；(3)小種群遺傳效應(yīng):遺傳漂變、瓶頸效應(yīng)和奠基者效應(yīng)等[39]。由此推測(cè)出東帕米爾高原喜馬拉雅雪雞可能由于缺乏基因交流和遺傳漂變導(dǎo)致核苷酸多樣性降低。鑒于此,喜馬拉雅雪雞作為國(guó)家Ⅱ級(jí)重點(diǎn)保護(hù)動(dòng)物,國(guó)家林業(yè)部門應(yīng)加強(qiáng)該物種的保護(hù),建議：首先擴(kuò)大位于東帕米爾高原的塔什庫(kù)爾干野生動(dòng)物自然保護(hù)區(qū)范圍并建立國(guó)家自然保護(hù)區(qū),恢復(fù)生態(tài)環(huán)境,提高雪雞棲息地的適宜性；其次在東帕米爾高原原地建立雪雞人工馴化繁育基地,借助分子遺傳標(biāo)記建立人工繁殖群體,以加強(qiáng)物種保護(hù)和利用；第三建立雪雞捕殺懲戒制度,對(duì)于人為捕殺者偷食者進(jìn)行嚴(yán)重懲戒。

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

對(duì)喜馬拉雅雪雞和GenBank中雪雞mtDNAD-loop共84個(gè)單倍型進(jìn)行系統(tǒng)發(fā)育和網(wǎng)絡(luò)關(guān)系分析,兩者聚類分析結(jié)果一致,84個(gè)單倍型聚為喜馬拉雅雪雞簇和藏雪雞簇。關(guān)于喜馬拉雅雪雞和藏雪雞起源研究,包新康從種內(nèi)亞種間分子系統(tǒng)發(fā)生關(guān)系、喜馬拉雅雪雞與藏雪雞分歧時(shí)間的估算以及高原及其周邊地質(zhì)環(huán)境演變方面認(rèn)定,喜馬拉雅雪雞和藏雪雞是在高原隆升前(也就是說(shuō)在3.6百萬(wàn)年前)就已經(jīng)分化形成了,在青藏高原隆升初期的多瑙河冰期(3.5百萬(wàn)—2.6百萬(wàn)年),淡腹組雪雞(藏雪雞又叫淡腹雪雞)的一支由天山東部地區(qū)向東南擴(kuò)散至青藏高原東部地區(qū)的一些高山,伴隨著青藏高原的隆升而形成藏雪雞；暗腹組雪雞(喜馬拉雅雪雞又稱暗腹雪雞)的一支形成喜馬拉雅雪雞后在高原隆升中期(大約0.88百萬(wàn)年前)開始向高原上擴(kuò)散并形成各亞種[40]。王小立研究了我國(guó)11個(gè)地理種群喜馬拉雅雪雞系統(tǒng)發(fā)育,結(jié)果表明喜馬拉雅雪雞存在4個(gè)進(jìn)化支,即昆侖山進(jìn)化支、青藏高原進(jìn)化支、天山進(jìn)化支和喀喇昆侖山進(jìn)化支,而喀喇昆侖山進(jìn)化支為東帕米爾高原(塔什庫(kù)爾干縣)所獨(dú)有,為獨(dú)立進(jìn)化支,不與其他地理群體雪雞共享進(jìn)化支[34],而本研究23個(gè)單倍型分布中19個(gè)單倍型(占84%)聚在以帕米爾高原為中心的進(jìn)化支上,這與王小立研究結(jié)果一致。但本研究有4個(gè)單倍型不與喀喇昆侖山進(jìn)化枝共享進(jìn)化支,暗示了東帕米爾高原喜馬拉雅雪雞經(jīng)歷過(guò)群體擴(kuò)張,具有兩個(gè)地理群體。帕米爾高原的群山起伏,連綿逶迤,冰川匯聚,雪峰群立,聳入云天,形成了獨(dú)特的地理區(qū)域,作為生長(zhǎng)在高海拔低溫地區(qū)的鳥類,這里形成的地理障礙使得雪雞基因交流受到限制,即形成獨(dú)立的地理類群或亞種,這一點(diǎn)在本研究得到進(jìn)一步證實(shí)。但東帕米爾高原喜馬拉雅雪雞是否經(jīng)歷過(guò)群體擴(kuò)張,由于受樣本數(shù)量的限制,尚需要進(jìn)一步驗(yàn)證。

[1] 李湘濤. 野生雞類的利用與保護(hù). 北京: 金盾出版社, 1997: 99-106.

[2] 鄭作新. 中國(guó)動(dòng)物志:鳥綱.第四卷,雞形目. 科學(xué)出版社, 1978.

[3] 馬鳴, 周永恒, 馬力. 新疆雪雞的分布及生態(tài)觀察. 野生動(dòng)物, 1991, (4): 15-16.

[4] 劉乃發(fā), 王香亭. 甘肅雞類分類區(qū)系研究. 甘肅科學(xué)學(xué)報(bào), 1992, 4(4): 29-38.

[5] 約翰·馬敬能, 卡倫·菲力普斯. 中國(guó)鳥類野外手冊(cè). 盧何芬, 譯. 長(zhǎng)沙: 湖南教育出版社, 2000.

[6] 沈孝宙, 王家駿. 中國(guó)雪雞的分類、地理分布和生態(tài). 動(dòng)物學(xué)雜志, 1963, 5(2): 67-68.

[7] IUCN. The IUCN Red List of Threatened SpeciesTM. Gland, Switzerland: IUCN. [2015-01-26]. http://www.iucnredlist.org.

[8] 蔣志剛, 江建平, 王躍招, 等. 中國(guó)脊椎動(dòng)物紅色名錄. 生物多樣性, 2016, 24(5): 500-551.

[9] 黃人鑫, 杜春華, 邵紅光, 張衛(wèi)紅. 中國(guó)高山雪雞的地理分布. 動(dòng)物學(xué)研究, 1992, 13(1): 66-66, 72-72.

[10] 黃人鑫, 馬力, 邵紅光, 姜濤. 新疆高山雪雞的生態(tài)和生物學(xué)的初步研究. 新疆大學(xué)學(xué)報(bào): 自然科學(xué)版, 1990, 7(3): 71-76.

[11] 盧汰春. 中國(guó)珍稀瀕危雉類. 大自然, 1989, (1): 37-39.

[12] 劉迺發(fā), 王香亭. 高山雪雞繁殖生態(tài)研究. 動(dòng)物學(xué)研究, 1990, 11(4): 299-302

[13] 常城, 劉迺發(fā), 王香亭. 暗腹雪雞的繁殖及食性. 動(dòng)物學(xué)報(bào), 1993, 39(1): 107-108.

[14] 阿德里, 馬鳴, 海肉拉, 王正已. 天山博格達(dá)南部雪雞的生態(tài)習(xí)性. 干旱區(qū)研究, 1997, 14(1): 84-87.

[15] 魏建輝, 陳玉平. 暗腹雪雞的生態(tài)習(xí)性初探. 甘肅林業(yè)科技, 2004, 29(4): 1-4, 8-8.

[16] 佚名. 新疆的珍貴野生動(dòng)物的飼養(yǎng)方法. 新疆林業(yè), 1975, (3): 34-34.

[17] 阿德力·麥地, 王德忠, 馬鳴, 海肉拉·麥得力汗. 天山博格達(dá)高山雪雞的生態(tài)習(xí)性與人工養(yǎng)殖技術(shù). 新疆畜牧業(yè), 1999, (4): 27-29.

[18] 阿德力·麥地, 王德忠, 馬鳴. 博格達(dá)高山雪雞的生物學(xué)特性及人工馴養(yǎng). 干旱區(qū)研究, 2000, 17(1): 75-77.

[19] 米來(lái)·夏日甫, 烏熱力哈孜, 賽務(wù)加甫, 冷青文, 海肉拉. 雪雞的馴化及人工養(yǎng)殖技術(shù). 中國(guó)家禽, 1999, 24(6): 31-32.

[20] 李世霞, 魏建輝. 暗腹雪雞的人工孵化技術(shù). 四川動(dòng)物, 2002, 21(1): 48-50.

[21] 王加福. 高山雪雞人工捕抓馴養(yǎng)繁殖適應(yīng)性研究. 甘肅畜牧獸醫(yī), 2001, 31(3): 12-13.

[22] 阿德力·麥地, 胡玉昆, 李凱輝. 暗腹雪雞人工飼養(yǎng)繁育初報(bào). 四川動(dòng)物, 2004, 23(2): 149-152.

[23] 常城, 劉迺發(fā), 王香亭. 暗腹雪雞青海亞種活動(dòng)規(guī)律及雛鳥羽毛生長(zhǎng)和成體秋季換羽. 甘肅科學(xué)學(xué)報(bào), 1994, 6(1): 77-81.

[24] 吾熱力哈孜, 米來(lái), 周自動(dòng), 海肉拉. 高山雪雞的內(nèi)臟器官解剖. 四川畜牧獸醫(yī), 2001, 28(7): 26-26, 28-28.

[25] 王輪, 王正己, 陳寧, 哈孜, 桂萍. 雪雞血液和生理參數(shù)的測(cè)定與分析(簡(jiǎn)報(bào)). 石河子農(nóng)學(xué)院學(xué)報(bào), 1995, (1): 63-64.

[26] 王正已, 王燕, 哈孜, 米來(lái), 許振營(yíng). 雪雞肉幾種氨基酸和微量元素的測(cè)定及對(duì)比分析. 石河子農(nóng)學(xué)院學(xué)報(bào), 1995, (1): 49-52.

[27] 文麗榮, 李國(guó)強(qiáng), 丁樹峰, 唐旭利. 新疆雪雞微量元素及氨基酸含量的測(cè)定. 新疆師范大學(xué)學(xué)報(bào): 自然科學(xué)版, 1996, 15(3): 44-46.

[28] 阿德力·麥地, 海茹拉. 新疆雪雞. 烏魯木齊: 新疆科技衛(wèi)生出版社, 1999: 110-125.

[29] 李紅燕, 王金富, 劉凱. 暗腹雪雞染色體的核型分析. 石河子大學(xué)學(xué)報(bào): 自然科學(xué)版, 2006, 24(3): 294-298.

[30] 劉強(qiáng), 吳敏, 張琳麟, 李悅. 暗腹雪雞細(xì)胞色素b基因的克隆及其在雉科中親緣關(guān)系的分析. 浙江大學(xué)學(xué)報(bào): 理學(xué)版, 2006, 33(1): 89-94.

[31] Ruan L Z, Zhang L X, Sun Q W, Liu N F. Genetic structure of Himalayan Snowcock (Tetraogallushimalayensis) populations in China. Biochemical Genetics, 2006, 44(9/10): 457-463.

[32] Ruan L Z, Zhang L X, Wen L Y, Sun Q W, Liu N F. Phylogeny and molecular evolution ofTetraogallusin China. Biochemical Genetics, 2005, 43(9/10): 507-518.

[33] 王小立. 喜馬拉雅雪雞的分子系統(tǒng)地理學(xué)研究[D]. 蘭州: 蘭州大學(xué), 2012.

[34] Wang X L, Qu J Y, Liu N F, Bao X K, Song S. Limited gene flow and partial isolation phylogeography of Himalayan SnowcockTetraogallushimalayensisbased on part mitochondrial D-loop sequences. Current Zoology, 2011, 57(6): 758-767.

[35] 王玉濤, 戴志剛, 楊世杰, 羅玉柱. 東帕米爾高原盤羊分布與棲息地植被覆蓋時(shí)空變化. 生態(tài)學(xué)報(bào), 2016, 36(1): 209-217.

[36] 傅伯杰, 于丹丹, 呂楠. 中國(guó)生物多樣性與生態(tài)系統(tǒng)服務(wù)評(píng)估指標(biāo)體系. 生態(tài)學(xué)報(bào), 2017, 37(2): 341-348.

[37] 安蓓. 西藏雪雞種群遺傳結(jié)構(gòu)與分子系統(tǒng)地理學(xué)研究[D]. 蘭州: 蘭州大學(xué), 2009.

[38] An B, Zhang L X, Browne S, Liu N F, Ruan L Z, Song S. Phylogeography of Tibetan Snowcock (Tetraogallustibetanus) in Qinghai-Tibetan Plateau. Molecular Phylogenetics and Evolution, 2009, 50(3): 526-533.

[39] Hunter, Malcolm L. Maintaining biodiversity in forest ecosystems. Cambridge University Press, 1999.

[40] 包新康. 我國(guó)高山鶉類分子系統(tǒng)發(fā)生研究[D]. 蘭州: 蘭州大學(xué), 2010.