基于核基因和微衛(wèi)星的葉城沙蜥遺傳多樣性分析

謝慧　蔣羅雪　趙偉　祁玥　段嘉寶　李孟純　劉筱燁　李鈾

摘要 為了研究我國(guó)特有物種葉城沙蜥（Phrynocephalus axillaris）的遺傳多樣性，采用微衛(wèi)星和核基因分子標(biāo)記方法，對(duì)新疆7個(gè)地區(qū)所采集的96只葉城沙蜥進(jìn)行了遺傳多樣性分析?；趯?duì)6個(gè)微衛(wèi)星位點(diǎn)和核基因BDNF的PCR擴(kuò)增及瓊脂糖凝膠電泳檢測(cè)，利用軟件進(jìn)行分析，結(jié)果表明群體單倍型多樣性為0.707和平均核酸多樣性為0.003 14，觀察到的雜合度平均值和預(yù)期雜合度平均值分別為0.52和0.65，總Fst值為0.175，發(fā)現(xiàn)試驗(yàn)研究所采集的葉城沙蜥群體存在明顯的遺傳差異及種群間存在較大的遺傳分化。綜合分析說明采集的葉城沙蜥群群體間遺傳交流水平低，存在較高程度的遺傳分化，但由于樣本較少仍需對(duì)葉城沙蜥種群遺傳特征進(jìn)一步研究。補(bǔ)充和豐富了新疆部分地區(qū)葉城沙蜥的基礎(chǔ)遺傳數(shù)據(jù)，為后續(xù)合理利用及保護(hù)該地區(qū)內(nèi)其他物種的種質(zhì)資源奠定了基礎(chǔ)。

關(guān)鍵詞 葉城沙蜥；微衛(wèi)星; 腦源性神經(jīng)營(yíng)養(yǎng)因子（BDNF）;遺傳;種群

中圖分類號(hào) Q953 文獻(xiàn)標(biāo)識(shí)碼 A

文章編號(hào) 0517-6611（2023）06-0081-07

doi：10.3969/j.issn.0517-6611.2023.06.021

Genetic Diversity Analysis of Phrynocephalus axillaris Based on Nuclear Gene and Microsatellite

XIE Hui1，JIANG Luo－xue1，ZHAO Wei2 et al

（1.Life Science and Engineering College，Northwest University for Nationalities，Lanzhou，Gansu 730100;2.College of Life Sciences，Lanzhou University，Gansu Provincial Key Laboratory of Environmental Biomonitoring and Restoration，Lanzhou，Gansu 730030）

Abstract In order to study the genetic diversity of the endemic species Phrynocephalus axillaris in China，the genetic diversity of 96 P.axillaris collected from seven regions of Xinjiang was analyzed by microsatellite and nuclear gene molecular markers.Based on the PCR amplification and agarose gel electrophoresis detection of six microsatellite loci and nuclear gene BDNF，the software was used for analysis.The results showed that the haplotype diversity of the population was 0.707 and the average nucleic acid diversity was 0.003 14.The average observed heterozygosity and expected heterozygosity were 0.52 and 0.65，respectively，and the total Fst value was 0.175.It was found that there were obvious genetic differences and large genetic differentiation among the populations collected in the experiment.Comprehensive analysis showed that the level of genetic communication among the collected populations was low，and there was a high degree of genetic differentiation.However，due to the lack of samples，it is still necessary to further study the genetic characteristics of the population.This study supplemented and enriched the basic genetic data of P.axillaris in some regions of Xinjiang，and laid the foundation for the subsequent rational utilization and protection of the germplasm resources of other species in the region.

Key words Phrynocephalus axillaris;SSR;BDNF;Heredity;Populations

葉城沙蜥（Phrynocephalus axillaris）是鬣蜥科（Agamidae）沙蜥屬（Phrynocephalus）的爬行動(dòng)物［1］，屬于我國(guó)特有的一種小型爬行動(dòng)物，已被列入《世界自然保護(hù)聯(lián)盟》（IUCN）2010年瀕危物種紅色名錄。廣泛分布于我國(guó)甘肅、西藏、新疆天山山脈南部地區(qū)，包括塔里木盆地及周邊的吐魯番—哈密盆地和敦煌盆地［2］。棲息在戈壁荒漠或沙漠邊緣地帶以及固定沙丘的丘間平地，有少數(shù)植被覆蓋，垂直分布的海拔高度為65～1 500 m，在阿爾金山南麓甚至可達(dá)3 000 m左右［3］。利用線粒體基因?qū)τ谌~城沙蜥種群遺傳結(jié)構(gòu)的研究未得到很好的解析，有可能與其分化時(shí)間太短或不完全譜系分選有關(guān)［4-5］。而基于微衛(wèi)星位點(diǎn)變異速率快、多態(tài)性高等特點(diǎn)，有望更全面地解析葉城沙蜥的種群遺傳結(jié)構(gòu)［6］。

微衛(wèi)星標(biāo)記，又稱簡(jiǎn)短串聯(lián)重復(fù)（short tandem repeats，STRs）或簡(jiǎn)單序列重復(fù)（simple sequence repeats，SSRs），是一種基于DNA長(zhǎng)度多態(tài)性的分子標(biāo)記技術(shù)。微衛(wèi)星DNA由于具有特異性的PCR擴(kuò)增、多態(tài)信息容量（PIC）高、引物通用性好、突變率高、共顯性、符合孟德爾遺傳、易于檢測(cè)等優(yōu)勢(shì)，已被廣泛應(yīng)用于種群遺傳、譜系地理和親子鑒定等領(lǐng)域［7-11］。微衛(wèi)星標(biāo)記較線粒體基因更具豐富的多態(tài)性，也是野生動(dòng)物遺傳學(xué)分析常用的分子標(biāo)記之一，研究者可以利用分子標(biāo)記研究種群的擴(kuò)散模式，從而更有效地保護(hù)和管理野生動(dòng)物。仲光啟等［12］篩選了在黑鸛上可擴(kuò)增的18個(gè)微衛(wèi)星位點(diǎn)并進(jìn)行了遺傳特點(diǎn)分析，王萌等［13］探討了微衛(wèi)星標(biāo)記方法在大型貓科動(dòng)物保護(hù)群遺傳學(xué)研究中的應(yīng)用與挑戰(zhàn)，劉旭等［14］概述了微衛(wèi)星技術(shù)在珍稀瀕危野生動(dòng)物中的相關(guān)應(yīng)用以及分子生物學(xué)手段在珍稀瀕危野生動(dòng)物擴(kuò)散生態(tài)學(xué)的研究情況。對(duì)于沙蜥屬的研究主要與環(huán)境氣候變化有關(guān)，有紅尾沙蜥超氧化物歧化酶適應(yīng)高原的分子機(jī)制［15］、青海沙蜥遺傳和表型分化的環(huán)境相關(guān)性［16］、增溫對(duì)青海沙蜥繁殖生活史及后代適合度的影響［17-18］等。目前關(guān)于葉城沙蜥中的研究相對(duì)較少，主要在兩性異形和形態(tài)特征差異、線粒體和基于高通量測(cè)序微衛(wèi)星特征分析［19］，而關(guān)于微衛(wèi)星分析研究較少。該研究利用高質(zhì)量的SSR標(biāo)記評(píng)估96只葉城沙蜥的遺傳多樣性，解析不同群體遺傳結(jié)構(gòu)，為葉城沙蜥及新疆地區(qū)的物種保護(hù)提供理論依據(jù)。該研究利用6個(gè)微衛(wèi)星位點(diǎn)和1對(duì)核基因引物，分析了7個(gè)地理種群96個(gè)樣本的遺傳多樣性、遺傳分化和系統(tǒng)地理結(jié)構(gòu)，對(duì)葉城沙蜥進(jìn)行系統(tǒng)發(fā)育分析，補(bǔ)充和豐富新疆地區(qū)葉城沙蜥的基礎(chǔ)遺傳數(shù)據(jù)，為保護(hù)該地區(qū)內(nèi)其他物種提供參考。

1 材料與方法

1.1 試驗(yàn)材料 試驗(yàn)研究對(duì)象為新疆7個(gè)地區(qū)（且末、民豐、策勒、和田、皮山、葉城、喀什）采集的96只葉城沙蜥，將組織樣本均置于無(wú)水乙醇與生理鹽水（1∶1配比）組織保存液中，-20 ℃保存?zhèn)溆茫ū?）。

1.2 基因組DNA的提取

使用通用型基因組DNA提取試劑盒進(jìn)行基因組DNA的提取，并用紅色熒光核酸染料染色的1.4%瓊脂糖凝膠電泳進(jìn)行DNA的完整性和濃度檢測(cè)，并于放置于-20 ℃冰柜保存?zhèn)溆谩?/p>

1.3 引物篩選及PCR擴(kuò)增

試驗(yàn)隨機(jī)選取4只葉城沙蜥組織作為篩選階段的材料，根據(jù)已發(fā)表的沙蜥文獻(xiàn)中引物信息及重新設(shè)計(jì)的葉城引物進(jìn)行篩選［6］，選出3對(duì)南疆沙蜥（Pvms33、Pvms39、Phr51）［20-21］和新設(shè)計(jì)的3對(duì)葉城沙蜥（PA11、PA36、PA37）微衛(wèi)星引物（表2）。PCR反應(yīng)體系為（12 μL）：5×buffer（Mg2+ free） 2.4 μL、5 U/μL Immolase DNA聚合酶0.06 μL、10 μmol/L熒光標(biāo)記上游引物tag F（Fam、Ned、Vic、Pet）0.09 μL、10 μmol/L熒光標(biāo)記下游引物tag R 0.09 μL、0.4 μmol/L Primers 2.4 μL、DNA模板2 μL、其余由dd H2O補(bǔ)足12 μL。PCR擴(kuò)增反應(yīng)條件：95 ℃預(yù)變性10 min;92 ℃變性60 s，50 ℃退火90 s，72 ℃延伸60 s，5個(gè)循環(huán)；92 ℃變性30 s，63 ℃退火90 s，72 ℃延伸60 s，20個(gè)循環(huán);92 ℃變性15 s，54 ℃退火60 s，72 ℃延伸60 s，40個(gè)循環(huán);最后72 ℃延伸10 min。擴(kuò)增產(chǎn)物進(jìn)行瓊脂糖凝膠電泳，并在凝膠成像儀上觀察結(jié)果。將成功擴(kuò)增的合格PCR產(chǎn)物送至蘇州金唯智生物科技有限公司進(jìn)行毛細(xì)管電泳分析。

根據(jù)GenBank公布的葉城沙蜥核基因（腦源性神經(jīng)營(yíng)養(yǎng)因子，brain-derived neurotrophic factor，BDNF）部分序列（序列號(hào)為KJ363410），使用Geneious Prime 2022.0.1（www.geneious.com）進(jìn)行引物設(shè)計(jì)，上游引物：BDNF-11F（5′-TGCATGAAAGCTGCCCCAAT-3′），下游引物BDNF-600R（5′-CTGGGTAGTTCGGCACTGAT-3′）。引物由艾科瑞生物（湖南）公司合成。使用設(shè)計(jì)的引物對(duì)葉城沙蜥樣本進(jìn)行擴(kuò)增。通過試驗(yàn)，建立最佳PCR擴(kuò)增體系（25 μL）：2×Pro Taq Master Mix（dye plus）12.5 μL、10 μmol/L上下游引物各0.5 μL、模板DNA 2.0 μL、dd H2O 9.5 μL。PCR反應(yīng)條件：94? ℃預(yù)變性3 min；94? ℃變性45 s，58 ℃退火45 s，72? ℃延伸1 min，共循環(huán) 35 次；72 ℃延伸10 min，4 ℃保存。擴(kuò)增產(chǎn)物經(jīng)1.4%瓊脂糖凝膠電泳進(jìn)行純度檢測(cè)。檢測(cè)合格的擴(kuò)增產(chǎn)物由蘭州天啟基因生物科技有限公司進(jìn)行測(cè)序。

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

試驗(yàn)中所得到微衛(wèi)星數(shù)據(jù)用GeneMarker［22］對(duì)測(cè)序儀得到的原始數(shù)據(jù)進(jìn)行基因分型分析。用GenAlex V 6.5［23］和Popgene32［24］軟件分析等位基因數(shù)（number of allele，Na）、有效等位基因數(shù)（effective number of allele，Ne）、觀測(cè)雜合度（observed heterozygosity，Ho）、期望雜合度（expected heterozygosity，He）、群體間Nei’s遺傳距離和遺傳相似度、群體內(nèi)近交系數(shù)（population inbreeding coefficient，F(xiàn)is）、群體間的遺傳分化指數(shù)（genetic differentiation index，F(xiàn)st）、基因流（gene flow，Nm）、Shannon多樣性指數(shù)（I）及方差分析（analysis of molecular variance，AMOVA）。用Cervus3.0軟件［25］計(jì)算各位點(diǎn)的多態(tài)信息含量（Poly-morphism information content，PIC），并進(jìn)行哈迪-溫伯格平衡（Hardy-Weinberg equilibrium，HWE）檢驗(yàn)。用MEGA X［26］軟件繪制基于群體間的Nei’s遺傳距離的UPGMA（unweighted pair－group method with arithmetic mean，UPGMA）聚類樹。

試驗(yàn)中所選取的核基因BDNF片段長(zhǎng)度在550～650 bp，PCR產(chǎn)物測(cè)序方式為正反向。采用Geneious 2020.2.5（http：//www.geneious.com）處理測(cè)序結(jié)果原始ABI文件，經(jīng)過人工校對(duì)后獲得96條葉城沙蜥BDNF基因樣本，序列長(zhǎng)度為524 bp?；诤嘶駼DNF分析葉城沙蜥7個(gè)種群的遺傳多樣性，主要包括單倍型多樣性（haplotype diversity，h）和核苷酸多樣性（nucleotide diversity，π）。下載Genbank上已發(fā)表的紅原沙蜥（P.hongyuanensis）的BDNF全基因（序列號(hào)為AF497714.1）作為外群的物種序列。使用MEGA-X軟件［26］構(gòu)建NJ（neighbor-joining）和ML（maximum likelihood）系統(tǒng)發(fā)育樹，進(jìn)化樹Bootstrap置信值為重復(fù)抽樣1 000次。

2 結(jié)果與分析

2.1 PCR擴(kuò)增結(jié)果

試驗(yàn)篩選出6對(duì)多態(tài)性高且穩(wěn)定的微衛(wèi)星及1對(duì)核基因引物，進(jìn)行后續(xù)遺傳多樣性分析。對(duì)選取微衛(wèi)星引物和核基因引物（其目的條帶在500～750 bp）在96只樣本上進(jìn)行PCR擴(kuò)增（核基因有1樣本未成功擴(kuò)增），采用凝膠電泳檢測(cè)，目的條帶符合預(yù)期大小且清晰，核基因序列未提交至Genbank，將成功擴(kuò)增的產(chǎn)物送公司進(jìn)行測(cè)序。以部分樣品在引物PA37和BDNF的電泳結(jié)果為例（圖1），電泳圖和測(cè)序結(jié)果一致，并且能準(zhǔn)確讀取具體的片段大小。

2.2 遺傳多樣性分析

6對(duì)SSR引物對(duì)7個(gè)葉城沙蜥群體進(jìn)行擴(kuò)增，各位點(diǎn)的平均等位基因數(shù)（Na）和有效等位基因數(shù)（Ne）分別為5.452和3.585個(gè)，平均觀測(cè)雜合度（Ho）、期望雜合度（He）、無(wú)偏預(yù)期雜合度（uHe）、群內(nèi)近郊系數(shù)（Fis）、群體間分化系數(shù)（Fst）和基因流（Nm）分別為0.515、0.650、0.734、0.082、0.175和3.473（表3）。7個(gè)葉城沙蜥群體的等位基因數(shù)存在一定差異，其中Na最多的為QM群體（7.000個(gè)），CL和KS群體較少（表4）。分析葉城沙蜥種群在微衛(wèi)星位點(diǎn)的遺傳多樣性參數(shù)，就位點(diǎn)而言，可知PA11位點(diǎn)上的平均Na、Ne、I、He較高，最低的是在Phr51（表3）；就種群而言，QM最高，KS最低，在雜合度方面，QM表現(xiàn)出較高的雜合性且有最為豐富的遺傳多樣性（表4）。

基于核基因BDNF對(duì)葉城沙蜥7個(gè)種群的遺傳多樣性分析，主要包括單倍型多樣性和核苷酸多樣性（表1）。首先分析單倍型多樣性，95只蜥蜴中共定義14個(gè)單倍型，群體平均單倍型多樣性為0.523，群體單倍型多樣性范圍為0.345～0.806，其中策勒（CL）單倍型多樣性最高，為0.806；葉城（YC）最低，為0.345。在核酸多樣性方面，葉城沙蜥平均核苷酸多樣性為0.001 72，核苷酸多樣性范圍為0.001 04～0.002 66，種群內(nèi)平均單倍型以策勒（CL）最高，為0.002 66，葉城（YC）最低，為0.001 04，表明葉城沙蜥整體遺傳多樣性為高單倍型多樣性和低核苷酸多樣性模式。

2.2.1 基于微衛(wèi)星標(biāo)記的遺傳分化。

葉城沙蜥7個(gè)群體間的Nei’s遺傳相似度和遺傳距離見表5。7個(gè)群體間的遺傳距離介于0.283～0.728，Nei’s遺傳相似度介于0.483～0.753。YC和MF之間遺傳距離最短，KS和CL之間的遺傳距離最長(zhǎng)。

基于群體間Nei’s遺傳距離的聚類樹結(jié)果（圖2），7個(gè)葉城沙蜥種群分為兩大類，第一大類包括PS和KS 2個(gè)種群，其他種群為第二大類。第二大類再分為3個(gè)亞類，第Ⅰ類包括QM，第Ⅱ類包括CL和HT，第Ⅲ類包括MF和YC。AMOVA分析結(jié)果顯示，7個(gè)群體內(nèi)的遺傳變異占總體變異的62%。群體間的遺傳變異占38%，分子方差A(yù)MOVA數(shù)據(jù)表明的變異大部分存在于群體內(nèi)，說明這7個(gè)群體之間遺傳分化較小。

2.2.2 BDNF基因葉城沙蜥的系統(tǒng)發(fā)育分析。

基于核基因NJ系統(tǒng)發(fā)育分析，由圖3a可知，95個(gè)葉城沙蜥群體樣本均在同一大分支，同時(shí)KS、YC和PS聚在同一枝上。紅原沙蜥作為外群，獨(dú)立于另一分支上?；诤嘶蚪ǔ傻腗L進(jìn)化樹，通過計(jì)算找到此次最大似然法構(gòu)建進(jìn)化樹的最佳模型為Kimura 2-paprameter model+Gamma Distributed，由圖3b可知，95個(gè)葉城沙蜥群體樣本均在同一大分支，同時(shí)KS、YC和PS聚在同一枝上，而MF、CL、HT和QM分別聚在兩支上。紅原沙蜥作為外群，獨(dú)立于另一分支上，結(jié)果基本一致。

3 討論

3.1 7個(gè)葉城沙蜥群體遺傳多樣性

遺傳多樣性研究是生物多樣性研究的重要組成部分，是生物進(jìn)化和物種分化的基礎(chǔ)，研究物種遺傳多樣性對(duì)于評(píng)估物種資源狀況及在制定野生物種資源保護(hù)策略等方面具有重要意義［27］。等位基因數(shù)、觀測(cè)雜合度、期望雜合度和多態(tài)信息含量等參數(shù)均反映群體遺傳基礎(chǔ)多樣性水平，其值越大，說明基因豐富度越高、生物對(duì)環(huán)境的適應(yīng)能力越強(qiáng)。趙丹等［28］利用微衛(wèi)星標(biāo)記對(duì)我國(guó)胎生蜥蜴的3個(gè)種群的18只母本和67只子代進(jìn)行遺傳多樣性分析，檢測(cè)到各種群母本和子代的觀測(cè)等位基因數(shù)在3.90～8.13。駱?biāo)?9］基于微衛(wèi)星對(duì)西藏高原的80個(gè)紅尾沙蜥和34個(gè)西藏沙蜥樣本進(jìn)行多樣性分析，分別檢測(cè)出44和42個(gè)等位基因，平均觀測(cè)雜合度分別為0.360 7和0.399 6。該研究利用6對(duì)微衛(wèi)星分子標(biāo)記對(duì)7個(gè)葉城沙蜥群體進(jìn)行遺傳多樣性分析，發(fā)現(xiàn)平均觀測(cè)雜合度、平均期望雜合度分別為0.515和0.650。哈迪-溫伯格平衡檢驗(yàn)可以作為衡量群體遺傳結(jié)構(gòu)的標(biāo)準(zhǔn)，該研究中，各個(gè)群體偏離HWE平衡，雜合體缺失和無(wú)效等位基因的存在是各個(gè)群體顯著偏離HWE平衡的原因。

該研究的7個(gè)葉城沙蜥種群在6個(gè)所用的微衛(wèi)星位點(diǎn)上檢測(cè)到的平均Na/Ne值大于1，檢測(cè)到的平均等位基因數(shù)高于平均有效等位基因數(shù)，說明等位基因分布不勻，比值較大可能是由近交的原因及樣本數(shù)量不夠?qū)е碌?。衡量群體遺傳變異的一個(gè)重要參數(shù)就是基因雜合度，一般期望雜合度（He）越高說明群體的遺傳變異越大；觀測(cè)雜合度（Ho）越接近期望雜合度（He），說明群體受外來(lái)因素影響越小，處于遺傳平衡狀態(tài)。而基于雜合度計(jì)算的近交系數(shù)（Fis）可以反映出群體內(nèi)部的近交程度，它的數(shù)值大小反映出了群體中雜合子的缺失或者過剩程度，該研究中表明葉城沙蜥群體間存在一定程度的基因流，這也可以解釋基于BDNF構(gòu)建的NJ和ML進(jìn)化樹中，葉城沙蜥的各種群沒有完全獨(dú)立于一支的原因。

核基因在系統(tǒng)發(fā)育、生物地理等方面研究的巨大潛力逐漸被發(fā)現(xiàn)，由于核基因進(jìn)化速率相對(duì)較慢，能夠反映更豐富的生物學(xué)信息，因此選擇BDNF基因作為研究葉城沙蜥種群進(jìn)化的分子標(biāo)記，能夠得到更豐富的結(jié)果?；贐DNF基因采用鄰接法（NJ）和最大似然法（ML）構(gòu)建的2種進(jìn)化樹結(jié)果幾乎一致，通過分析可知，在葉城沙蜥種群中，7個(gè)地區(qū)共計(jì)95個(gè)樣本均在1大分支上，喀什地區(qū)與同市的葉城和皮山的葉城沙蜥親緣關(guān)系最近，且末、民豐、策勒和部分和田地區(qū)的葉城沙蜥群體在同一大支，與基于遺傳距離構(gòu)建的UPGMA圖大致一致，但其之間是否存在基因交流的個(gè)例無(wú)法確定，有待進(jìn)一步研究。

3.2 葉城沙蜥群體遺傳分化和遺傳結(jié)構(gòu)

遺傳分化參數(shù)Fst最大值為1表明兩個(gè)群體完全分化，最小值為0表明群體間無(wú)分化，在實(shí)際研究中，F(xiàn)st在0～0.05時(shí)，表明物種群體間的遺傳分化水平很低，可以不做考慮；在0.05～0.15時(shí)，表明群體間的遺傳分化水平處于中等程度；在0.15～0.25時(shí)群體間存在較大的遺傳分化；0.25以上時(shí)群體間就存在很大的遺傳分化［30］。該研究揭示的7個(gè)葉城沙蜥群體間遺傳分化指數(shù)Fst為0.175，分子方差A(yù)MOVA數(shù)據(jù)表明38%的變異存在于群體間，62%的變異存在于群體內(nèi)，表明群體間存在中等偏高程度的遺傳分化。

4 結(jié)論

研究使用多重?zé)晒釶CR及瓊脂糖凝膠電泳的方法，使用6對(duì)微衛(wèi)星引物和BDNF基因?qū)π陆貐^(qū)7個(gè)葉城沙蜥群體進(jìn)行遺傳多樣性分析。研究發(fā)現(xiàn)7個(gè)葉城沙蜥群體在6個(gè)微衛(wèi)星位點(diǎn)上的近交系數(shù)均值為0.082，Shannon多樣性指數(shù)為1.336；遺傳分化指數(shù)均值為0.175，基因流均值為3.473，表明群體間遺傳交流水平低，存在較高程度的遺傳分化；群體單倍型多樣性范圍為0.345～0.806及核苷酸多樣性范圍為0.001 04～0.002 66。目前我國(guó)葉城沙蜥群體的遺傳多樣性處于中等水平，群體間的遺傳分化受自身遷移能力和地理隔離等因素影響，亟需對(duì)我國(guó)葉城沙蜥進(jìn)行種質(zhì)資源保護(hù)，防止其遺傳多樣性繼續(xù)下降。

哈迪-溫伯格平衡是評(píng)價(jià)種群遺傳平衡的參數(shù)，處于遺傳平衡的種群不容易受到近交及外來(lái)種群的干擾［31］，種群相對(duì)穩(wěn)定。瀕危物種常常表現(xiàn)HWE平衡偏離，造成偏離的原因很多，如近親交配、種群遷移和遺傳漂變等。該研究中有6個(gè)位點(diǎn)偏離HWE平衡，均為雜合不足，且表現(xiàn)近親繁殖。對(duì)新疆7個(gè)地區(qū)葉城沙蜥的微衛(wèi)星和核基因進(jìn)行擴(kuò)增測(cè)序進(jìn)行分析，研究其遺傳多樣性，但目前的數(shù)據(jù)無(wú)法證明葉城沙蜥種群表現(xiàn)的遺傳特征及原因，有待進(jìn)一步研究，該研究補(bǔ)充和豐富了新疆部分地區(qū)葉城沙蜥的基礎(chǔ)遺傳數(shù)據(jù)，可為后續(xù)合理利用及保護(hù)該地區(qū)內(nèi)其他物種的種質(zhì)資源奠定基礎(chǔ)。

參考文獻(xiàn)

［1］ LI J，GUO X G，WANG Y Z.Population genetic differentiation of Phrynocephalus axillaris in east of Xinjiang Uygur Autonomous Region based on sequence variation of mitochondrial ND4－tRNALeu gene［J］.Zool Res，2010，31（4）：361-369.

［2］ 趙爾宓，趙肯堂，周開亞，等.中國(guó)動(dòng)物志：爬行綱 第二卷 有鱗目 蜥蜴亞目［M］.北京：科學(xué)出版社，1999.

［3］ 張鵬，袁國(guó)映.新疆兩棲爬行動(dòng)物［M］.烏魯木齊：新疆科學(xué)技術(shù)出版社，2005.

［4］ ZHANG Q，XIA L，HE J B，et al.Comparison of phylogeographic structure and population history of two Phrynocephalus species in the Tarim Basin and adjacent areas［J］.Mol Phylogenet Evol，2010，57（3）：1091-1104.

［5］ 丁新華.基于線粒體和核基因分子標(biāo)記的黑鰭鳈親緣地理研究［D］.上海：上海海洋大學(xué)，2020.

［6］ 宋琪，劉金龍，郭憲光.基于Roche 454 GS FLX高通量測(cè)序的葉城沙蜥基因組微衛(wèi)星特征分析［J］.四川動(dòng)物，2019，38（1）：62-67.

［7］ 張敏瑩，方弟安，周彥鋒，等.基于微衛(wèi)星標(biāo)記的錢塘江中下游三角魴親本群體和自然捕撈群體遺傳多樣性及遺傳結(jié)構(gòu)［J］.大連海洋大學(xué)學(xué)報(bào)，2022，37（5）：775-783.

［8］ 趙蒙蒙，滕暉，吳?；?，等.用于番鴨親子鑒定的微衛(wèi)星分子標(biāo)記初步篩選［J］.黑龍江畜牧獸醫(yī)，2021（5）：1-7，156.

［9］ 陳政昆，韋聰聰，朱雯，等.安徽白山羊保種群微衛(wèi)星標(biāo)記遺傳多樣性分析［J］.黑龍江畜牧獸醫(yī)，2020（9）：54-56，63.

［10］ 張代臻，唐伯平，孫紅英，等.DNA分子標(biāo)記在蝦蟹類系統(tǒng)地理學(xué)研究中的應(yīng)用［J］.水產(chǎn)科學(xué)，2008，27（5）：262-265.

［11］ VASHISTHA G，DEEPIKA S，DHAKATE P M，et al.The effectiveness of microsatellite DNA as a genetic tool in crocodilian conservation［J］.Conserv Genet Resour，2020，12（4）：733-744.

［12］ 仲光啟，由玉巖，徐艷春，等.黑鸛微衛(wèi)星位點(diǎn)的篩選及遺傳特點(diǎn)分析［J］.野生動(dòng)物學(xué)報(bào)，2019，40（3）：729-36.

［13］ 王萌，朱思雨，薛茂盛，等.微衛(wèi)星標(biāo)記方法在大型貓科動(dòng)物保護(hù)群遺傳學(xué)研究中的應(yīng)用與挑戰(zhàn)［J］.野生動(dòng)物學(xué)報(bào)，2019，40（3）：780-786.

［14］ 劉旭，丁由中.分子生物學(xué)技術(shù)應(yīng)用于野生動(dòng)物保護(hù)及擴(kuò)散研究概況［J］.野生動(dòng)物學(xué)報(bào)，2019，40（2）：497-501.

［15］ 魯松松，黃宇哲，張夢(mèng)洋.紅尾沙蜥超氧化物歧化酶適應(yīng)高原的分子機(jī)制［J］.蘭州大學(xué)學(xué)報(bào)（自然科學(xué)版），2022，58（1）：69-78，88.

［16］ 謝凡.青海沙蜥遺傳和表型分化的環(huán)境相關(guān)性［D］.南京：南京師范大學(xué)，2021.

［17］ 王坤.增溫對(duì)青海沙蜥繁殖生活史及后代適合度的影響［D］.保定： 河北大學(xué)，2020.

［18］ 祁玥.南疆沙蜥分子系統(tǒng)地理學(xué)［D］.蘭州：蘭州大學(xué)，2019.

［19］ 阿爾孜古·衣米提，徐峰，李俊.塔里木盆地西南部葉城沙蜥的兩性異形研究［J］.干旱區(qū)資源與環(huán)境，2021，35（7）：190-195.

［20］ URQUHART J，BI K，GOZDZIK A，et al.Isolation and characterization of microsatellite DNA loci in the toad－h(huán)eaded lizards，Phrynocephalus przewalskii complex［J］.Mol Ecol Notes，2005，5（4）：928-930.

［21］ ZHAN A B，F(xiàn)U J Z.Microsatellite DNA markers for three toad－h(huán)eaded lizard species （Phrynocephalus vlangalii，P.przewalskii and P.guttatus）［J］.Mol Ecol Resour，2009，9（2）：535-538.

［22］ HOLLAND M M，PARSON W.GeneMarker? HID：A reliable software tool for the analysis of forensic STR data［J］.J Forensic Sci，2011，56（1）：29-35.

［23］ PEAKALL R，SMOUSE P E.GenAIEx 6.5：Genetic analysis in Excel.Population genetic software for teaching and research：An update［J］.Bioinformatics，2012，28（19）：2537-2539.

［24］ ROUSSET F.GENEPOP’ 007：A complete re－implementation of the genepop software for Windows and Linux［J］.Mol Ecol Resour，2008，8（1）：103-106.

［25］ KALINOWSKI S T，TAPER M L，MARSHALL T C.Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment［J］.Mol Ecol，2007，16（5）：1099-1106.

［26］ KUMAR S，STECHER G，LI M，et al.MEGA X：Molecular evolutionary genetics analysis across computing platforms［J］.Mol Biol Evol，2018，35（6）：1547-1549.

［27］ 魏輔文，馬天笑，胡義波.中國(guó)瀕危獸類保護(hù)遺傳學(xué)研究進(jìn)展與展望［J］.獸類學(xué)報(bào)，2021，41（5）：571-580.

［28］ 趙丹，劉玉芬，趙文閣，等.基于微衛(wèi)星標(biāo)記的胎生蜥蜴母本與子代的遺傳多樣性研究［J］.哈爾濱師范大學(xué)自然科學(xué)學(xué)報(bào)，2017，33（4）：50-53.

［29］ 駱?biāo)?西藏高原兩種沙蜥遺傳結(jié)構(gòu)及三種沙晰親子鑒定［D］.蘭州：蘭州大學(xué)，2009.

［30］ 馮一寧，李因剛，祁銘，等.基于SSR標(biāo)記的福建省閩楠代表性群體遺傳多樣性分析［J］.南京林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），2022，46（4）：102-108.

［31］ 楊建寶，張艷萍，蘇軍虎，等.引進(jìn)美洲紅點(diǎn)鮭群體遺傳多樣性微衛(wèi)星的分析［J］.四川動(dòng)物，2012，31（3）：348-352.

基金項(xiàng)目 國(guó)家自然科學(xué)基金地區(qū)科學(xué)基金項(xiàng)目（32060311）；甘肅省自然科學(xué)基金項(xiàng)目（20JR10RA124）；西北民族大學(xué)中央高校基本科研業(yè)務(wù)費(fèi)項(xiàng)目（31920210143，31920200004，31920190083）；西北民族大學(xué)研究生科研創(chuàng)新項(xiàng)目（Yxm2021076）。

作者簡(jiǎn)介 謝慧（1997—），女，河南商丘人，碩士研究生，研究方向：動(dòng)物遺傳育種與繁殖。*通信作者，副教授，博士，碩士生導(dǎo)師，從事分子生態(tài)與進(jìn)化生物學(xué)研究。

收稿日期 2022-07-27