陕南水牛微卫星DNA遗传多样性研究

陕南水牛微卫星DNA遗传多样性研究王斌1,2，昝林森1,3[收稿日期]2023-01-17[基金项目]国家“十一五”科技支撑计划项目（2023BAD14B07）；西北农林科技大学拔尖人才支持计划项目[作者简介]王斌（1979—）男，陕西城固人，讲师，硕士，重要从事生物技术与家畜育种研究。E－mail：[[收稿日期]2023-01-17[基金项目]国家“十一五”科技支撑计划项目（2023BAD14B07）；西北农林科技大学拔尖人才支持计划项目[作者简介]王斌（1979—）男，陕西城固人，讲师，硕士，重要从事生物技术与家畜育种研究。E－mail：[通信作者]昝林森（1963—），陕西扶风人，专家，博士生导师，重要从事牛遗传育种与繁殖研究。E－mail：（1西北农林科技大学动物科技学院，陕西杨凌712100；2汉中职业技术学院农林系，陕西汉中723000；3国家肉牛改良中心，中国杨凌7121004；4河南师范大学生命科学学院，河南新乡453007）[摘要]【目的】分析4个陕南水牛群体的遗传多样性，为陕南水牛品种的资源保护和开发运用提供支持。【方法】运用10对微卫星引物，采用PCR扩增和非变性聚丙稀酰胺凝胶电泳技术，对192头陕南水牛（分别采自城固、南郑、宁强和汉滨4个群体，每群体48头）进行了遗传多样性检测，记录了各群体的等位基因频率、等位基因数、有效等位基因数（Ne）、遗传杂合度（H）、各群体间的奈氏标准遗传距离及多态信息含量（PIC），根据遗传距离进行了聚类分析。【结果】4个陕南水牛群体在10个微卫星位点共发现104个等位基因，其中有8个特有等位基因，等位基因频率为0.0052～0.3490，总群体各位点平均有效等位基因数为4.0174～13.1469；10个微卫星位点平均多态信息含量为0.7007～0.8932，均为高度多态；平均杂合度为0.7550～0.9288。聚类分析结果显示，城固群体与南郑群体一方面聚在一起，然后依次同宁强、汉滨水牛群体聚在一起。【结论】10对微卫星标记可作为有效的遗传标记用于陕南水牛遗传多样性的分析；陕南水牛群体变异限度多样性丰富，选育限度较低；聚类分析结果符合陕南水牛的地理分布格局和产区的自然环境，4个群体之间遗传差异较小,遗传一致性限度较大。[关键词]陕南水牛；微卫星；遗传多样性[中图分类号][文献标记码]A[文章编号]StudyongeneticdiversityofmicrosatelliteDNAinShaannanbuffaloWANGBin1,2,ZANLin-sen1,3,YANGYan-jie4,ZHONGXin1,HUANGLei1,WANGHong-bao1,3(1CollegeofAnimalScienceandTechnology，NorthwestA&FUniversity，Yangling，Shaanxi712100,China；2A&FDepartment，HanzhongVocationandTechnologyCollege，Hanzhong,Shaanxi723000,China；3NationalBeefCattleImprovementCenterinChina,Yangling,Shaanxi712100,China;4CollegeofLifeSciences,HenanNormalUniversity，Xinxiang，Henan453007,Abstract:【Objective】ToanalysegeneticdiversityoffourlandracesbuffaloinsouthernShaanxiandtoprovidesupportfortheprotectionandutilizationofresourcesofbuffalospeciesinsouthernShaanxi.【Method】192ShaannanBuffalo（wereselectedfromChenggu,Nanzheng,NingqiangandHanbin）wereusedtodetectthegeneticdiversity,andcountedallelefrequencyofeachgroup,numberofalleles,effectivenumberofalleles(Ne),heterozygosity(H),betweengroupsNei'sstandardgeneticdistanceandthepolymorphicinformationcontent(PIC),analysizegeneticdistancethroughclusteranalysis.【Result】Theresultsindicatethatthereare104allelesand8alleles-specificat10microsatellitelocusin4localbuffalopopulations.Theallelefrequencies;averageeffectivenumberofallelesandtheaveragepolymorphicinformationcontentrangingfrom0.0052－0.3490;4.0174－13.1469and0.7007－0.8932,meanwhile,theaverageheterozygosityrangefrom0.7550－0.9288.theywerehighlypolymorphic.ClusteranalysisshowedthatChenggugroupsfirstofalltogetherwithNanzheng,followedwithNingqiang,Hanbin.【Conclution】10pairsofmicrosatellitemarkerscanbeusedasaneffectivegeneticmarkerfortheanalysisofthegeneticdiversityofbuffaloinsouthernShaanxi.TherearerichdiversityofvariationinShaannanbuffalo.AnanlysisofclusteringshowedthattheseconsistentwiththegeographicaldistributionpatternsandnaturalenvironmentofShaananBuffalo.TherearebasicallysameatgeneticbasisbetweenFourlocalgroups(ChengguNanzhengNingqiangandHanbin)Keywords：Shaannanbuffalo；Microsatellite；geneticdiversity微卫星DNA(MicrosatelliteDNA)，又称短串联反复(Shorttandemrepeats，STR)或简朴反复序列(Simplesequencerepeat，SSR)，是以1～6个核苷酸为反复单位串联组成的长达几十个核苷酸的反复序列[1]，如（TG）n、（AAT）n、（GATA）n等。微卫星广泛存在于各种真核生物的基因组中，并且分布比较均匀，平均每10kb的DNA序列中就会出现1个微卫星序列[2]。微卫星具有数量大、分布广而均匀、多态性丰富、共显性遗传以及分析方法简便、快捷、结果反复性好等优点，目前已在动植物的基因组作图、群体遗传变异分析、群体起源分析、亲子鉴定、QTL检测与定位等研究中得到广泛应用[3-5]。陕南水牛是主产于汉中和安康盆地的地方水牛品种，重要分布于汉中、安康两市的十三县（区）[6],具有体质健壮、力大、役用性能好、耐粗饲、耐苦、耐劳、抗病性强等特点，有一定的乳、肉生产潜力，是适宜于陕南水田密集区饲养的一个重要畜种资源。目前关于陕南水牛分子水平遗传多样性的研究还比较少。张伟[7]对6个水牛类群54个个体（其中6头为陕南水牛）的mtDNAD-loop区全序列进行测序分析，从母系角度研究了中国水牛的起源。齐国强等[8]对9头陕南水牛的mtDNAD-loop区序列(930bp左右)、类群之间的系统进化关系以及遗传多样性进行了分析。张毅[9]运用线粒体D-loop区序列和Y染色体SNP两类标记,分别从母系和父系两个方面，对涉及陕南水牛在内的水牛群体进行了遗传多样性和起源驯化分析。截止目前，尚未见有关运用微卫星遗传标记从分子水平对陕南水牛进行遗传多样性和发生关系的研究报道。为此，本研究运用10个微卫星位点，对陕南水牛进行遗传学检测，初步探讨其遗传多样性和遗传分化关系，以期为陕南水牛在中国水牛中系统地位的拟定及进一步开展陕南水牛品种资源保护和开发运用工作提供参考。1材料与方法 1.1材料从陕南水牛中心产区城固（CG）、南郑（NZ）、宁强（NQ）、汉滨（HB）4地分别选择48头品种特性明显、没有与国外水牛杂交历史，并且个体间在3代或2代内没有血缘关系的个体，颈静脉采血（10mL/头），ACD抗凝（V(ACD)︰V(血液)＝1︰6），轻轻晃动摇匀，置于冰盒带回实验室，于－80℃1.2重要试剂蛋白酶K、TaqDNA聚合酶、10×Buffer、dNTPs、DNAMarker，均购自北京鼎国生物工程公司。分子量标记物pBR322/MspⅠ，由北京华美生物公司提供。1.3陕南水牛基因组DNA的制备参照《分子克隆实验指南》[10]，采用传统的酚、氯仿抽提法并略作改善，提取陕南水牛血样基因组DNA。将DNA用TE缓冲液溶解，然后用15g/L琼脂糖进行凝胶电泳检测，并用紫外分光光度计测定其浓度，于－80℃保存备用。1.4PCR扩增从2023年国际动物遗传学会（ISAG）和联合国粮农组织（FAO），推荐的用于家养水牛遗传多样性研究的微卫星引物中，筛选出10对引物（表1）用于实验，引物由上海生物工程技术服务有限公司合成。PCR反映体系为15μL：10×Buffer(含15mmol/LMg2+)1.5μL，2mmol/LdNTPs1.5μL，2U/μLTaqDNA聚合酶0.25μL，10pmol/μL上、下游引物各0.6μL，模板DNA(50ng/μL)1.0μL，ddH2O10.15μL。PCR扩增程序：95℃预变性5min；94℃30s，55～65.8℃（退火温度因引物不同而异，表1）退火35s，72℃30s，32个循环；72℃表110对微卫星引物的PCR条件Table1PCRrectionconditionsof10pairsprimerofmicrosatellite位点Locus染色体位置Chromosomalposition引物序列（5’→3’）Primersequence退火温度/℃Annealingtemeperature产物大小/bpProductsCSSM03317CACTGTGAATGCATGTGTGTGAGC65.8155～209CCCATGATAAGAGTGCAGATGACTILSTS01211TCTACCACCGATACAGATGG61.1109～125GAAGTAGGTAGTGCTGGAGGILSTS03124AATTCTAGGTGAACAGCAGC62.3276～304AAGACATACTCTCAGACTCCILSTS06820TGCAAAGAGTTGGATGGGGC61.5146～212AAGGTTGCTGCCTCATATCCHEL1311TAAGGACTTGAGATAAGGAG55.0137～193CCATCTACCTCCATCTTAACILSTS06115AAATTATAGGGGCCATACGG55.0133～171TGGCCTACCCTACCATTTCCILSTS09521GAAAGATGTTGCTAGTGGGG58.6207～247ATTCTCCTGTGAACCTCTCCCSSMO191TTGTCAGCAACTTCTTGTATCTTT59.9137～171TGTTTTAAGCCACCCAATTATTTGCSSM03627GGATAACTCAACCACACGTCTCTG59.9166～180AAGAAGTACTGGTTGCCAATCGTGCSSM04122AATTTCAAAGAACCGTTACACAGC56.0136～154AAGGGACTTGCAGGGACTAAAACA1.5PCR产物的检测用20g/L琼脂糖凝胶电泳检测PCR产物，将检测阳性的PCR产物用120g/L非变性聚丙烯酰胺凝胶于120V条件下电泳7～8h，银染显带并用紫外凝胶成像系统照相分析。以pBR322/MspⅠ为分子量标记物，用Gel-ProAnalyzer3.1分析软件鉴定个体等位基因片段大小及其基因型。1.6数据记录分析等位基因频率(Alelefrequency)、等位基因数（NumberofAlleles）、有效等位基因数(EffectiveNumberofAlleles，Ne)、遗传杂合度(Heterozygosity，H)、奈氏标准遗传距离(Standardgeneticdistance，Ds)的计算和聚类分析，均运用Popgen32软件进行；多态信息含量(Polymorphisminformationcontent，PIC)用PIC分析软件进行计算。2结果与分析2.1陕南水牛群体基因组DNA的检测将DNA用TE缓冲液溶解后，进行15g/L琼脂糖凝胶电泳检测，结果（图1）显示，DNA条带较亮，比较整齐，没有拖尾现象，可以进行下一步实验。图图1陕南水牛血样基因组DNA的检测注：图中1－12分别为12头水牛个体血样基因组DNAFig.1GelelectrophoresisofShaannanbuffalogenomicDNANote:1－12indicategenomicDNAof12buffaloinFig.RespectivelyFig1.2.2陕南水牛群体等位基因的组成及等位基因频率的计算城固、南郑、宁强及汉滨4个水牛群体在10个微卫星位点共发现104个等位基因，等位基因频率为0.0052～0.3490，其中8个为特有等位基因。8个特有等位基因中，城固水牛群体ILSTS012位点上有1个、ILSTS068位点上有3个，ILSTS095位点上有1个，该群体共有5个特有等位基因；南郑水牛群体在ILSTS068和ILSTS061位点上各有1个特有等位基因；宁强水牛群体在CSSM033位点上有1个特有等位基因。其中陕南汉滨水牛群体19头水牛个体微卫星ILSTS095的PCR扩增结果如图2所示。12123456789101112131415M16171819404bp404bp309bp242bp217bp207bp180bp图2汉滨（HB）水牛群体19头水牛个体微卫星ILSTS095PCR产物的电泳图谱M.Marker;1~19.19个汉滨水牛个体Fig.2PolyacrylamidegelelectrophoresisofmicrosatelliteILSTS095PCRproductsofBuffaloSouthernHanBin(HB)groupsof19buffalobuffaloM.Marker;1~19.19buffalosofHanbingroup2.3陕南水牛群体的多态信息含量（PIC）PIC是等位基因频率和等位基因数的变化函数，通常用来表达微卫星座位的杂合限度。PIC不仅是衡量等位基因片段多态性的抱负指标，即可以用来表达微卫星DNA变异限度的高低，同时也关系到一个座位的可用性和效率。PIC越大，该座位杂合子的比例越大，提供的遗传信息就越高。由表2可知，4个陕南水牛群体10个位点的多态信息含量为0.7007～0.8932，平均为0.8104。表24个陕南水牛群体在10个微卫星座位上的多态信息含量Table2Polymorphisminfomraitoncontentat10microsatellitelocusoffourlocalgroups位点LocusCGNZNQHB均值MeanCSSM0330.77540.75640.69850.77430.7512ILSTS0120.72130.64060.74430.69660.7007ILSTS0310.82780.81900.81590.78050.8108ILSTS0680.91060.91450.84700.90080.8932HEL130.85340.83380.85560.84700.8475ILSTS0610.84230.84640.85560.84290.8468ILSTS0950.84570.83690.83390.81440.8327CSSM0190.86240.87270.87330.84300.8629CSSM0360.77220.72590.74120.75160.7477CSSM0410.78820.73970.77030.77390.7680均值Mean0.82350.80510.80730.80570.81042.4陕南水牛群体的等位基因数与有效等位基因数（Ne）Ne是反映群体遗传变异大小的一个指标，一般用纯合度的倒数来表达，反映了等位基因间的互相影响。等位基因在群体中分布越均匀，有效等位基因数越接近所检测到的等位基因的绝对数。陕南水牛各群体在10个微卫星座位上的等位基因数见表3。由表3可知，ILSTS068座位的等位基因数最多，达成23个；CSSM036座位的等位基因数最少,只有6个;其余座位的等位基因数为8～11个。城固、南郑、宁强和汉滨4个水牛群体10个位点的平均等位基因数依次分别为9.3，8.9，9.1和8.5个。陕南水牛总群体在10个微卫星座位上的有效等位基因数见表4。由表3和表4可知，陕南水牛各群体的有效等位基因数与检测到的等位基因数的绝对数接近限度一致，表白等位基因在各群体中分布均匀。表34个陕南水牛群体在10个微卫星座位上的等位基因数Table3Numberofallelesat10microsatellitelocusof4localbuffalopopulationsinShaannan位点LocusCGNZNQHB总群体TotalpopulationCSSM03377879ILSTS01276657ILSTS03188878ILSTS0681616161523HEL131110111111ILSTS0611011101011ILSTS0951089811CSSMO191010101010CSSM03666556CSSM04187878均值Mean9.38.99.18.510.4表44个陕南水牛群体在10个微卫星座位上的有效等位基因数Table4Effectivenumberofallelesat10microsatellitelocusof4localbuffalopopulationsinShaannan位点LocusCGNZNQHB总群体TotalpopulationCSSM0335.05264.66403.66885.05265.0693ILSTS0124.14393.29144.51763.82724.0174ILSTS0316.50856.22706.12775.12006.3801ILSTS06812.000012.521710.867910.867913.1469HEL137.43236.69777.15537.20007.4714ILSTS0617.06757.20007.68007.06757.6008ILSTS0957.20006.85716.73686.06327.3143CSSMO198.00008.59708.66177.06758.7108CSSM0365.05264.21984.48254.66404.9034CSSM0415.30884.38024.82014.94425.1214均值Mean6.77666.46566.47186.18746.97362.5陕南水牛群体的遗传杂合度遗传杂合度（Heterozygosity，H)是某一座位上杂合子的盼望频率，可用于衡量某座位上遗传变异量的大小。某座位上的等位基因数目越多，各等位基因频率互相越接近，则遗传杂合度值越大，说明该座位上的遗传变量越大、越丰富。平均遗传杂合度（Meanheterozygosity)是某一群体各座位遗传杂合度的算术平均值。由表5可知，4个陕南水牛群体在各位点均表现出丰富的遗传多样性。表54个陕南水牛群体在10个微卫星座位上的遗传杂合度和平均遗传杂合度Table5heterozygosityandaverageterozygosityat10microsatellitelocusofgroupsandtotalpopulationofShaannanbuffalo位点Locus遗传杂合度Heterozygosity平均遗传杂合度AverageheterozygosityCGNZNQHBCSSM0330.81910.80230.74290.81910.8069ILSTS0120.77480.71100.79520.75440.7550ILSTS0310.86440.85730.85460.82180.8477ILSTS0680.93620.93970.92730.92730.9288HEL130.88390.86880.87850.87940.8707ILSTS0610.87680.87940.88830.87680.8730ILSTS0950.87940.87230.86970.85280.8678CSSM0190.89360.90250.90340.87680.8898CSSM0360.81910.77930.79340.80230.8002CSSM0410.82890.78810.80940.81470.8090均值Mean0.85760.84010.84630.84260.84492.6陕南水牛群体间的遗传距离应用GENEPOP软件计算城固、南郑、宁强及汉滨4个陕南水牛群体的奈氏标准遗传距离，所得结果见表6。表64个陕南水牛群体间的奈氏标准遗传距离Table6Nei’sstandardgeneticdistanceof4localbuffalopopulationsinShaannan群体GroupsCGNZNQHBCG****0.93870.92070.9154NZ0.0632****0.90280.8990NQ0.08260.1022****0.8898HB0.08840.10640.1168****由表6可知，城固、南郑、宁强和汉滨4个陕南水牛群体的奈氏标准遗传距离从大到小依次为宁强与汉滨(0.1168)、南郑与汉滨(0.1064)、南郑与宁强(0.1022)、城固与汉滨(0.0884)、城固与宁强(0.0826)及城固与南郑(0.0632)。2.7陕南水牛群体的聚类分析采用POPGENE软件，根据遗传资料计算城固、南郑、宁强和汉滨4个陕南水牛群体在10个微卫星座位上的遗传距离,绘制UPGMA聚类图，其结果见图3。由图3可知,城固水牛群体与南郑水牛群体一方面聚在一起，然后依次同宁强、汉滨水牛群体聚在一起。+------------------------------城固CG+------------1+---2+------------------------------南郑NZ!!--3+-------------------------------------------宁强NQ!+------------------------------------------------汉滨HB图3图34个陕南水牛群体间的UPGMA亲缘关系树状图Fig.3UPGMAdendrogramamong4localbuffalopopulationsinShaannan3讨论3.1样本的含量样本量是实验结果可靠性的主线保证，太多则会导致不必要的浪费，太少则无法得到对的的结果。Mathur等[11]指出，应用DNA指纹标记时，每群体或品系有13～20个个体就可以获得遗传变异的合理估计值。本实验因水牛无品种之分，各地方良种水牛有共同的遗传性和生物学特性，只有在大量标记、遗传杂合度较低即品种间遗传距离较大的情况下，小样本才可接受。张继全等[12]指出，遗传距离估计精度随取样个体数的增长而提高。对此，每个群体采集48头样本，有效保证了实验结果的可靠性，特别是对发现特有等位基因具有重要意义。3.2等位基因和特有等位基因本实验10个位点192个样本发现等位基因104个，特有等位基因8个，与叶鼎承[13]的报道结果（8个位点130个样本，发现等位基因78个，特有等位基因28个）相比，特有等位基因较少，表白陕南水牛各群体选育方向一致。从特有等位基因分布来看，城固群体有5个、南郑群体有2个、宁强群体有1个。各群体中只有南郑群体的2个个体在ILSTS068位点具有186bp的特有等位基因，其余各特有等位基因均只有1个个体具有，即均为稀有等位基因。该结果表白，陕南水牛城固群体受其他地方良种水牛的影响更大，但并不影响遗传距离的测定。张继全等[12]认为，稀有等位基因对品种即遗传关系确立的奉献极小。3.3群体遗传多样性分析Botstein等[14]提出了衡量基因变异限度高低的多态信息含量指标，当某位点的多态信息含量>0.5时，表白该位点为高度多态位点；当其值为0.25～0.5时，该位点为中度多态位点；当其值<0.25时为低度多态位点。本研究10个位点的多态信息含量为0.7007～0.8932，平均为0.8104，说明本研究所选的10个微卫星位点均为高度多态位点，可以用于遗传多样性的分析；同时也表白，陕南水牛多样性丰富，有很大的运用、选择潜力。此外，各位点的多态信息含量差异不大，相对较为集中，标准差为0.0614。由此可见，陕南水牛各位点相关性状均未受到强大选择压力的影响，且选择方向一致。许明[15]对中国中东部地区5个地方（峡江、海子、西林、恩施、江汉）水牛群体群体遗传多样性的研究表白，其多态信息含量为0.6425～0.6675，平均为0.6514±0.0101。刘建[16]对贵州6个地方水牛群体遗传多样性的研究显示，其多态信息含量为0.6476～0.7767，平均为0.7196±0.0270。本研究中4个陕南水牛群体的多态信息含量为0.8051～0.8235，平均0.8104±0.0088，说明与中国中东部5个地方水牛群体及贵州6个地方水牛群体遗传多样性相比，4个陕南水牛群体的遗传多样性更为丰富。此外，陕南水牛各群体多态信息含量的标准差仅为0.0088，互相之间高度接近，表白陕南水牛4个群体均未受到人工选择的强大影响，且所受到的选择压力强度基本一致。有效等位基因数越接近于所检测到的等位基因的绝对数，表白选育限度越高，人工选择压力越大。杨国锋等[17]报道，波尔山羊为高度选育品种，其有效等位基因数和检测到的等位基因数仅相差0.68，而川东白山羊为2.66。叶鼎承[13]、许明[15]、刘健[16]分别对福安水牛、中国中东部5个地方（峡江、海子、西林、恩施、江汉）水牛群体和贵州6个地方水牛群体进行了研究，各群体检测等位基因数与有效等位基因数之差的平均值分别为2.6275，4.2838和1.5135，而本研究为各水牛群体检测到的等位基因数与有效等位基因数之差的平均值为2.5996，说明陕南水牛有效等位基因数高于贵州水牛群体而低于中国中东部5个地方水牛群体，与福安水牛的有效等位基因数基本一致。该结果表白，陕南水牛的选育限度及人工选择压力同福安水牛基本一致，但高于贵州水牛而低于中国中东部5个地方水牛群体。本研究表白，陕南水牛各位点的遗传杂合度为0.7550～0.9288，平均遗传杂合度为0.8449±0.0515，各位点的遗传杂合度均大于0.5，说明陕南水牛具有较高的变异度，群体近交限度弱，具有丰富的多样性。据王冬蕾[18]报道，德昌、兴隆、富钟、温州、东流、福安水牛群体的平均遗传杂合度为0.6763～0.7194，平均为0.6890±0.0163；许明[15]研究发现，中国中东部地区5个地方水牛群体的遗传杂合度为0.7029～0.7580，平均为0.7304±0.0259；刘建[16]发现贵州6个水牛群体的杂合度为0.7999～0.7450，平均为0.7730±0.0152。本研究各群体遗传杂合度为0.8401～0.8576，平均为0.8467±0.0031。由此可见，陕南水牛的变异限度高于以上各地方水牛群体。这也许是由于陕南虽位于中国水牛分布的最北部，但整体沿汉江为中轴线呈东西分布，因而受到多个地方水牛群体的影响。从地理位置上来看，其同涪陵水牛、恩施水牛、江汉水牛中心产区较为接近，但是否受到它们的影响，尚有待进一步研究。3.4群体间遗传距离与聚类分析遗传距离是研究物种遗传多样性的基础，反映了所研究群体的系统发生关系，可用来描述群体的遗传结构和品种间的差异，其中Nei氏提出的标准遗传距离被广泛应用于群体遗传变异的分析。有研究者比较了几种不同的聚类方法,认为UPGMA是最为稳健的一种[19]。本研究表白，4个陕南水牛群体的遗传距离为0.0632～0.1168，这与王冬蕾[18]对我国6个水牛群体及张毅[9]对我国18个水牛群体的研究结果相似。本研究表白，城固水牛群体与南郑水牛群体一方面聚在一起，然后依次同宁强、汉滨水牛群体聚在一起，上述聚类关系基本符合陕南水牛的地理分布格局和产区的自然环境。本研究得到的多态信息含量、有效等位基因数及杂合度均一致表白，陕南水牛多样性丰富，选育限度较低。聚类分析表白，城固、南郑、宁强和汉滨4个陕南水牛群体之间的遗传距离基本处在同一水平,表白这4个水牛群体之间的遗传差异较小,遗传一致性限度较大。[参考文献][1]ParkerPG,SmowAA,SchugMD,etal.Whatmoleculescantellusaboutpoupulations:Choosingandusingamolecularmarker[J].Ecology,1998,79(2):361.[2]TautzD,HypervariabilityofsimplesequencesasageneralsourceforpolymorphicDNAmarkers[J].NuclAcidsRes,1989,17:6463.[3]ArranzeJJ,BayonY,PrimitivoFS,etal.Comparisionofproteinmarkersandmicrosatellitesindifferentiationofcattlepopulation[J].AnimGenet,1996,27:415-419.[4]BarkerJSF，MooreSS,HetzelDJS,etal.Geneticdiversityofasianwaterbuffalo(Bubalusbubalis):Microsatellitevarietyonandacomparisonwithproteincodingloci[J].AnimGenet,1997,28:103-115.[5]GoudarizKM,LsloeD,FuretJP,etal.Analysisofgeneticrelationshipbetween10cattlebreedswith17microsatellite[J].AnimalGenetics,1997,8:338-345.[6]邓银才,李锡章,禚振宇，等.秦巴山区家畜家禽及经济动物品种志[M].北京:中国农业科技出版社,1990.DengYC,LiXZ,XiangZY.etal.PoultyandlvestockeconomicanimalspeciesofQin-baMountainareas[M].Beijing：ChinaAgricultureScienceandTechnologyPress,1990.(inChinese)[7]张伟.中国6个地方水牛类群mtDNA遗传多样性与分子系统进化研究[D].陕西杨凌：西北农林科技大学，2023.ZhangW.China’ssixlocalbuffalomtDNAgeneticdiversityandmolecularphylogeny[D].Yangling，Shaanxi:NorthwestA&FUniversity，2023.(inChinese)[8]齐国强，昝林森，张桂芳，等.中国部分地方水牛品种mtDNAD-loop区遗传多样性与起源研究[J].畜牧兽医学报，2023,39(1):7-11.QiGQ,ZanLSZhangGF,etal.StudyongeneticdiversityandoriginofChineseBuffalomtDNAD-looparea[J].VeterinaryJournal，2023,39(1):7-11.(inChinese)[9]张毅.中国水牛分子遗传多样性及其起源分化的研究[D].北京:中国农业大学，2023.ZhangY.StudyongeneticdiversityandoriginofChinese[D].Beijing:ChinaAgricultureUniversity,2023.([10]萨姆布鲁克等.分子克隆实验指南[M].黄培堂,译.3版.北京：科学出版社,2023.Sombrook,J.etal.ExperimentGuideofMolicularColon[M].Huangpei-tang,Translated.3edition.Beijing:SciencePress,2023.[11]MathurPK,PonsuksillS,GroenAF,etal.Estimationofgeneticvariabilitywithinand