Genetic diversity and differentiation of five Cuban cattle breeds using 30 microsatellite loci

Conservation and improvement strategies in farm animals should be based on a combination of genetic and phenotypic characteristics. Genotype data from 30 microsatellites were used to assess the genetic diversity and relationships among five Cuban cattle breeds (Siboney de Cuba, Criollo Cubano, Cebú Cubano, Mambí de Cuba and Taíno de Cuba). All microsatellite markers were highly polymorphic in all the breeds. The expected heterozygosity ranged from 0.67 ± 0.02 in the Taíno de Cuba breed to 0.75 ± 0.02 in the Mambí de Cuba breed, and the observed heterozygosity ranged from 0.66 ± 0.03 in the Cebú Cubano breed to 0.73 ± 0.02 in the Siboney de Cuba breed. The genetic differentiation between the breeds was significant (p < 0.01) based on the infinitesimal model (F_ST). The exact test for Hardy‐Weinberg equilibrium within breeds showed a significant deviation in each breed (p < 0.0003) for one or more loci. The genetic distance and structure analysis showed that a significant amount of genetic variation is maintained in the local cattle population and that all breeds studied could be considered genetically distinct. The Siboney de Cuba and Mambí de Cuba breeds seem to be the most genetically related among the studied five breeds.     

柴达木绒山羊八个微卫星位点遗传多样性分析

以柴达木绒山羊为研究对象,对其8个微卫星位点进行遗传多样性分析,计算各位点的等位基因频率(P)、杂合度(H)、多态信息含量(PIC)和有效等位基因数(N)。结果表明:8个微卫星位点在柴达木绒山羊均呈高度的多态性(PIC>0.5);位点LSCV 24基因杂合度、有效等位基因和多态信息含量都最高(H=0.796,N=6.916,PIC=0.766);位点BM 3413的基因杂合度、多态信息含量都最小(H=0.676,PIC=0.615);位点MAF-70的有效等位基因最小(N=0·7672)。因此,柴达木绒山羊群体中遗传背景复杂,遗传多样性较丰富;这8个微卫星位点也可用于柴达木绒山羊羊绒性状的进一步研究。     

10个云南地方猪种微卫星遗传多样性分析

为阐明保山猪、迪庆藏猪、高贡黎山猪、丽江猪、明光小耳猪、滇南小耳猪、撒坝猪、大河猪、昭通猪、大河乌猪10个云南省优良地方猪种的群体遗传多样性和遗传结构,本试验参考ArkDB数据库中家猪14条染色体上的15对微卫星引物,对549个云南地方猪个体的耳组织样品进行了检测分析,计算各微卫星座位等位基因数（Na）、有效等位基因数（Ne）、Shannon's信息指数（I）、表观杂合度（Ho）、期望杂合度（He）、F-统计量（Fis、Fit、Fst）、基因流（Nm）、群体遗传距离等相关参数及多态信息含量（PIC）,采用NTsys 2.10软件构建聚类树,并采用STRUCTRUE 2.3.3软件评估群体遗传结构。结果显示：15个微卫星座位共检测到293个等位基因,各座位等位基因数在5~38个之间,平均等位基因数为19个,平均有效等位基因数8.3682个;10个云南地方猪群体的Shannon多样性、表观杂合度、期望杂合度和多态信息含量分别在1.4374~2.0317、0.6389~0.7756、0.7021~0.8281和0.6647~0.7993之间。各座位单个群体内近交系数（Fis）、总群体近交系数（Fit）、群体间遗传分化系数（Fst）的变化范围分别为-0.0678~0.4594、0.0618~0.5567和0.0713~0.1801;群体间Fst平均值为0.1101,处于中度遗传分化状态,有11.01%的遗传变异来自于群体间,大部分遗传变异来自于群体内;每个座位基因流程度较高,变化范围在1.6392~3.2551之间,平均值为2.0214。基于Nei氏遗传距离构建UPGMA系统发生树显示,高黎贡山猪与迪庆藏猪聚为一支,并与丽江猪、保山猪和明光小耳猪聚为一大支,其结果得到了STRUCTURE分析的验证。综上所述,10个云南地方猪种遗传多样性丰富,群体内有近交现象,群体间处于中度遗传分化,存在着一定的基因交流。     

闽南黄牛的微卫星多态性研究

：用8个微卫星标记对闽南黄牛的等位基因频率、多态信息含量、杂合度和有效等位基因数等指标进行统计分析,并在此基础上对其进行聚类分析和分类研究。结果：8个微卫星标记在所检测的闽南黄牛群体中都表现出较丰富的多态性,在4个群体8个微卫星座位共有58个等位基因,每个微卫星标记平均检测到7．2个等位基因（3～11）;8个微卫星标记的平均多态信息含量（PIC）为0．6471,各群体的平均多态信息含量（PIC）差并不显著;全部群体平均杂合度为0．2930;各群体平均有效等位基因数为3．37。此表明闽南黄牛在微卫星位点上具有丰富的遗传多样性;根据奈氏标准遗传距离进行UPGMA聚类分析,结果表明4个群体间的遗传变异不大。     

利用基因组重测序技术鉴定黑安格斯牛的纯度

【目的】基于10头黑安格斯牛和60头对照组牛的全基因组重测序数据,分析黑安格斯牛的纯度、遗传多样性及群体结构。【方法】全基因组重测序技术和生物信息学方法。【结果】通过对10头黑安格斯牛全基因组数据进行分析,共检测到15,064,459个SNP位点,其核苷酸多样性（pi）为0.0015,观测杂合度（Ho）为0.2381,期望杂合度（He）为0.2430,表明黑安格斯牛的遗传多样性较低;主成分分析和群体遗传结构分析发现,黑安格斯牛与欧洲普通牛聚为一类,其中有4头黑安格斯牛存在偏离情况,表明这4头牛主要与中国瘤牛与东亚普通牛韩牛之间存在杂交。【结论】10头黑安格斯牛中,6头为纯种黑安格斯牛,4头为杂种牛。     

4个肉牛品种微卫星多态性分析

本研究探讨了引进牛种（西门塔尔牛、利木赞牛）和山东地方黄牛（鲁西黄牛、利鲁牛）在20个微卫星位点上的多态性。试验从不同地区的种公牛站采集鲁西黄牛、利鲁牛、利木赞牛和西门塔尔牛的血液或精液样本,分别为56、27、31和30个,共计144个样本,血液样本采用Lab-Aid基因组DNA分离试剂盒提取基因组DNA,精液样本采用高盐法提取基因组DNA。利用20个微卫星标记,采用荧光标记毛细管电泳方法,对4个肉牛品种的DNA多态性进行分析,主要研究了每个群体的遗传变异指标（等位基因数、多态信息含量和杂合度）及群体间的遗传关系（F-统计量和基因流）。结果表明,肉牛基因组DNA条带整齐,无拖尾现象,质量较好,可用于本试验。20个微卫星座位中共检测到211个等位基因,平均等位基因数为10.6个。群体的平均观测杂合度在0.6147~0.7176之间,平均期望杂合度在0.6735~0.7862之间。在所有群体中各位点的多态信息含量在0.4853~0.8714之间,平均为0.7284,但在各个群体内的差异较大。近交系数及基因流分析结果表明,4个肉牛品种近交系数较低,群体间平均近交系数为0.085,每个微卫星位点的基因流均>1。总体来看,所选用的20个微卫星座位多态信息含量高,可作为有效遗传标记用于肉牛品种间遗传多样性分析。     

新疆13个绵羊群体遗传多样性及遗传分化的研究

利用10个微卫星标记对新疆13个绵羊群体的遗传多样性和遗传分化进行了分析。通过计算多态信息含量（PIC）、平均杂合度（H）、群体内近交系数（Fis）、遗传分化系数（Fst）、总群体近交系数（Fit）和基因流（Nm）等参数,评估各品种遗传多样性和品种间遗传分化。结果表明：13个绵羊群体10个微卫星标记的平均多态信息含量为0.6803;平均杂合度为0.7192,其中中国美利奴羊平均杂合度最高为0.7530,山区和田羊最低为0.6754;13个绵羊群体总近交系数为0.1066,群体内近交系数为0.0565,群体间基因分化系数为0.0531,说明5.65%的遗传变异来自群体间,94.35%的遗传变异来自群体内个体间的差异;基因流Nm平均值为4.4621（3.0808~7.4589）,说明各种群间存在或者在过去某个时期发生基因交流;各种群间平均遗传距离较低,为0.1702,平均遗传一致度数值较高,为0.8448,也再次证明了新疆13个绵羊群体间分化程度不高。基于Nei′s标准遗传距离运用UPGMA聚类法获得的新疆13个绵羊群体的聚类图与其品种育成史和地理分布基本相一致。     

Genetic structure of wild and domesticated grasscutters (Thryonomys swinderianus) from south-western Nigeria

Grasscutter (Thryonomys swinderianus) is a large rodent distributed across sub-Saharan Africa that is highly valued as a source of protein. There has been little effort to evaluate the genetic structure of grasscutters despite long-term harvesting pressure and over 40 years of grasscutter domestication in West Africa. Our objectives were to quantify the genetic structure of wild grasscutters, and to compare genetic variation from wild samples to those from various farmed samples within south-western Nigeria. We genotyped 145 wild and 88 domesticated individuals at 11 microsatellite loci and present results quantifying regional genetic structuring and the relative patterns of diversity among wild and domesticated grasscutter populations. Our data reflect high differentiation between wild and domesticated grasscutters, and significantly greater allelic richness and gene diversity in the former. Despite this, domesticated populations appear to have similar levels of observed heterozygosity and comparable levels of differentiation among domesticated samples relative to wild samples. This may be the result of high turnover within captive colonies, or frequent infusion of new animals. More detailed molecular and quantitative genetic studies are recommended on this species to be able to understand their natural variation, degree of connectivity and to improve strategies for domestication.     

中国6个双峰驼群体微卫星遗传多样性分析

为了解中国双峰驼群体的遗传多样性及不同种群间的遗传进化关系,本研究采用微卫星标记技术,对中国阿拉善驼、青海驼、南疆驼、北疆驼、肃北驼、苏尼特驼6个双峰驼群体进行了遗传多样性分析。通过计算杂合度（H）、多态信息含量（PIC）、有效等位基因（Ne）、Shannon信息指数等分析群体内遗传变异,通过计算F-统计量、基因流、遗传分化系数、遗传距离等分析群体间遗传进化关系。结果显示,10个微卫星位点共检测到了89个等位基因,平均每个位点检测到8.9个等位基因;所有位点均属中高度多态位点（YWLL08除外）,平均PIC值在0.488~0.752之间;6个群体观测杂合度值（0.355~0.448）都低于期望杂合度值（0.643~0.703）;几乎所有位点的Shannon指数都>1,且处于哈代-温伯格不平衡状态（P< 0.05）。群体间遗传分化系数Fst值为0.059,处于较低程度的中等分化状态; 6个双峰驼群体的平均Fis值均为正值,说明6个双峰驼群体都存在不同程度的近交。基于标准遗传距离D_S和遗传距离D_A进行聚类分析,南疆驼和北疆驼聚为一支,阿拉善驼、青海驼、肃北驼、苏尼特驼4个群体聚为一支。研究表明,中国双峰驼遗传多样性丰富,群体内遗传变异较大,存在一定近交现象;群体间存在着一定的基因流动,群体间的分化主要由群体内的遗传变异造成;6个群体分为2个类群。     

Genetic analysis of Mexican hairless pig populations.

The main Creole pig population in Mexico, the hairless Mexican pig, remains as an unimproved and endangered genetic resource. In order to learn more about the genetic characteristics of this pig population, we assessed the allele frequency of 10 microsatellite loci in 177 unrelated hairless pigs from seven regions at Mexico and in 111 pigs of four commercial breeds (Landrace, Large White, Hampshire, and Duroc). Genetic diversity in each population was estimated by the unbiased average heterozygosity and the allele number. Nei's standard genetic distances and a neighbor-joining dendrogram were used to reveal the genetic relationships among these populations. In this report, we present data showing that the level of the genetic diversity in Mexican hairless pigs is high compared with previous reports, and that they belong to a genetic lineage divergent from commercial breeds. Furthermore, Mexican hairless pigs seem to have developed several genetically distinct lines associated with their geographic location. We conclude that the Mexican Creole pig populations may be a reservoir of genetic diversity that is important to preserve and evaluate as a source of new alleles for the future improvement of commercial pig lines.     

Analysis on the Genetic Diversity of Six Bactrian Camel Populations in China Using Microsatellite Marker

In order to understand the genetic diversity of Bactrian camel and the genetic evolutionary relationships between different populations,the genetic diversity of Alashan camel, Qinghai camel, Nanjiang camel, Beijiang camel, Subei camel and Sunite camel populations in China were analyzed using 10 microsatellite markers. By calculating heterozygosity (H), polymorphism information content (PIC), effective allele (Ne), Shannon information index, the genetic variation within populations were analyzed. By calculating the F-statistics, gene flow, genetic differentiation coefficient and genetic distance to analyze the genetic relationship between populations. The results showed that 89 alleles were tested at the 10 microsatellite loci,8.9 alleles were tested at every locus in average. All loci were medium-highly polymorphic loci (except YWLL08),the average PIC of the Bactrian camel population was between 0.488 to 0.752. The observed heterozygosity (0.355 to 0.448) of 6 Bactrian camel populations was lower than the expected heterozygosity (0.643 to 0.703). Almost all loci Shannon index was greater than 1,and most of loci were in Hardy-Weinberg disequilibrium. The genetic differentiation coefficient between groups (Fst) value was 0.059 and in the low degree of moderately differentiated state. The average Fis values of 6 Bactrian camel populations were positive which suggested 6 Bactrian camel populations had different levels of inbreeding. The standard genetic distance (D_S) and genetic distance T (D_A) clustering analysis showed that the Nanjiang camel and Beijiang camel jointed as one group, the Alxa camel, Qinghai camel, Subei camel, Sunite camel jointed as the other group. Research showed that Chinese Bactrian camel was abundant in genetic diversity, genetic variation within population was larger, and there was a phenomenon of inbreeding. There was a certain gene flow between populations, the differentiation of populations were mainly caused by genetic variation within population. 6 Chinese Bactrian camel populations were divided into two groups.     

The Genetic Polymorphism Analysis in Tree Shrew(T.belangeri chinensis) based on Fluorescent Microsatellite

A total of 10 fluorescent labeling microsatellite markers detected by capillary electrophoresis were selected to analyze the genetic diversity of tree shrew(T.belangeri chinensis).Microsatellite Toolkit and PopGene softwares were used to calculate the genetic correlation indexes to assess the genetic diversity of the population.64 alleles with 6.4 alleles per locus were identified among 10 microsatellite loci.The average effective number of alleles, observed heterozygosity, expected heterozygosity and polymorphic information content were 3.7892, 0.6156, 0.6963, and 0.6367, respectively.Seven of all loci were high polymorphic, and three were moderate polymorphic.Five loci extremely significantly deviated from Hardy-Weinberg equilibrium(P<0.01).These novel polymorphic microsatellite markers indicated high genetic diversity, which will be used to evaluate genetic diversity in tree shrew.     

荧光标记微卫星分析广西树鼩群体的遗传多态性

选取10个微卫星位点,结合荧光标记(FAM) PCR扩增技术,运用毛细管电泳检测扩增产物,利用Microsatellite Toolkit、PopGene软件计算遗传相关指标,探讨广西树鼩群体的遗传多态性。结果发现,在10个微卫星基因座中,共检测到64个等位基因,平均等位基因数为6.4个,平均有效等位基因数为3.7892,平均观察杂合度和平均期望杂合度为0.6156和0.6963,平均多态信息含量为0.6367,其中7个是高度多态位点,3个属中度多态位点。此外,在10个微卫星基因座中,有5个位点偏离Hardy-Weinberg平衡(P<0.01)。结果表明广西树鼩遗传多态性较丰富,所选微卫星位点基本可用于评估广西树鼩群体遗传多样性。     

Evaluation of the genetic structure of sika deer (Cervus nippon) in Japan's Kanto and Tanzawa mountain areas,based on microsatellite markers

The browsing habits of sika deer (Cervus nippon) in Japan have caused serious ecological problems. Appropriate management of sika deer populations requires understanding the different genetic structures of local populations. In the present study, we used 10 microsatellite polymorphisms to explore the genetic structures of sika deer populations (162 individuals) living in the Kanto region. The expected heterozygosity of the Tanzawa mountain range population (Group I) was lower than that of the populations in the Kanto mountain areas (Group II). Our results suggest that moderate gene flow has occurred between the sika deer populations in the Kanto mountain areas (Group II), but not to or from the Tanzawa mountain range population (Group I). Also, genetic structure analysis showed that the Tanzawa population was separated from the other populations. This is probably attributable to a genetic bottleneck that developed in the Tanzawa sika deer population in the 1950s. However, we found that the Tanzawa population has since recovered from the bottleneck situation and now exhibits good genetic diversity. Our results show that it is essential to periodically evaluate the genetic structures of deer populations to develop conservation strategies appropriate to the specific structures of individual populations at any given time.     

Detection of occurrence of a recent genetic bottleneck event in Indian hill cattle breed Bargur using microsatellite markers

The effective number of breedable individuals is a crucial determinant for maintaining genetic variability within a population. The population of Bargur, the hill cattle of South India, has gone down drastically by more than 93 % in the past three decades, and only a few thousand animals are available at present. The present study was undertaken to evaluate Bargur cattle for mutation drift equilibrium and to detect the occurrence of recent genetic bottleneck event, if any, in this population. About 50 unrelated animals, true to the type, were sampled and genotyped at 25 microsatellite loci. The mean observed heterozygosity (0.808?±?0.023) was higher than the mean expected heterozygosity (0.762?±?0.008) with 15 out of 25 microsatellite loci exhibiting heterozygosity excess when assumed under Hardy–Weinberg equilibrium. To evaluate Bargur cattle for mutation drift equilibrium, three tests were performed under three different mutation models, viz., infinite allele model (IAM), stepwise mutation model (SMM) and two-phase model (TPM). The observed gene diversity (H _e) and expected equilibrium gene diversity (H _eq) were estimated under different models of microsatellite evolution. All the 25 loci were found to exhibit gene diversity excess under IAM and TPM, while 22 loci were having gene diversity excess under SMM. All the three statistical tests, viz., sign test, standardized differences test, and Wilcoxon sign rank test, revealed significant (P?<?0.01) deviation of Bargur cattle population from mutation-drift equilibrium under all the three models of mutation. Furthermore, the qualitative test of allele frequency distribution in Bargur cattle population revealed a strong mode shift from the normal L-shaped form suggesting that the population had experienced genetic bottleneck in the recent past. The occurrence of genetic bottleneck might have led to the loss of several rare alleles in the population, which point towards the need for efforts to conserve this important cattle germplasm. The present study is the first report in demonstrating the genetic basis of demographic bottleneck in an Indian cattle population.     

湘西黄牛遗传多样性的初步研究

利用6个微卫星标记对湘西黄牛的遗传多样性进行了初步分析.6对微卫星引物在湘西黄牛中共检测到65个等位基因,平均每对引物上检测到9.3个等位基因;平均有效等位基因7.8759个,等位基因在湘西黄牛群体中的分布很均匀;平均观测杂合度为0.8669,平均期望杂合度为0.8829.6个微卫星座位在湘西黄牛中平均多态信息含量(PIC)为0.8297,均大于0.5,都是高度多态性的.湘西黄牛杂合子的比例较大,基因纯合率都偏低,群体的遗传变异度高,表明湘西黄牛遗传多样性丰富,具有较大的选择潜力,可承受的选择压比较大,这是开展新品种培育工作的有利条件.     

塔里木马鹿微卫星遗传多样性与产茸量的相关性研究

利用7个微卫星遗传标记,采用PCR扩增,12%非变性聚丙烯酰胺凝胶电泳、Sanguinetti银染法显色,对新疆塔里木马鹿亚种遗传多样性进行检测,统计各群体的等位基因组成、平均有效等位基因数(E)和平均基因纯合率(Rh),根据等位基因频率计算出各群体的平均遗传杂合度(h)、多态信息含量(PIC)和群体间的遗传距离,利用分子进化遗传分析软件(MEGA),采用邻结法(NJ)重建系统发生树。根据等位基因频率,利用PHYLIP(3.6)分析软件,采用最大似然法(ML)构建系统发生树,应用自举检验(bootstrap test)估计系统树中节点的自引导值(bootstrap value),并进行了系统发生分析。同时对部分马鹿群体个体基因型与产茸量之间进行了相关分析。结果表明:7个微卫星位点在3个塔里木马鹿中的多态信息含量除BMS2508和Celjp0023没有多态性,以及BM5004为中度多态外,其他4个微卫星均为高度多态,可作为有效的遗传标记用于3个塔里木马鹿群体遗传多样性和系统发生关系的分析。所有马鹿群体的平均PIC为(0.5196)、h(0.5552)和E(2.45),其基因多态性和遗传多样性相对丰富。塔里木马鹿35团群体与沙雅群体的血缘关系及遗传距离近于阿拉尔群体,3个塔里木马鹿群体的系统发育关系基本符合其地理分布和育成史。微卫星位点BM4208的166bp/185bp基因型和微卫星BM888的208bp/208bp基因型可以对塔里木马鹿35团群体高产茸量进行分子遗传标记。     

利用微卫星分析康西草原马匹种群的遗传多样性和群体结构

旨在揭示康西草原马匹的品种来源以及群体遗传结构。本研究采集了当地60匹健康状况良好、体重适中、3～6岁之间不区分公母马的血液样品,试验分为6个种群并以纯血马等4个国外马品种、1个我国地方马品种作为对照， 通过设计引物并提取基因组DNA结合微卫星扩增法对染色体上12个微卫星位点的等位基因进行检测，并对各种群遗传多样性参数、群体分化指数、瓶颈效应以及群体结构进行分析。试验按照群体个数为依据分为6组，每组1个重复，每个重复以该组样本量为标准。分析6个种群马的遗传参数结果表明，其中共检测到221个等位基因,平均Shannon指数为1.691；有效等位基因数(Ne)为3.649～5.397；期望杂合度(He)为0.681～0.798；观测杂合度(Ho)为0.632～0.780；平均多态信息含量(PIC)为0.643～0.772。这些结果表明，6个群体都具有很高的遗传多样性。通过微卫星DNA位点的连锁不平衡及哈迪-温伯格平衡分析发现， 大多数位点的P值都大于校正值(P=0.006 5)，表示较多位点都处于哈迪-温伯格及连锁平衡状态。各群体间群体分化系数(Fst)表明，群体间分化程度较低(Fst<0.15)。基于Nei’s遗传距离与地理距离之间Mantel相关性检验显示二者之间呈现正相关性(R=0.502)，但未达到显著水平(P=0.310)；各马匹种群瓶颈效应分析表明，P值极显著(P_(IAM)<0.01)说明该地区的马群体受到过瓶颈效应的影响且该马群体历史上数量有大规模的减少。通过 FCA因子分析、基于Nei’s标准遗传距离UPGMA树状图以及Structure贝叶斯聚类分析发现，该地区马种群与纯血马和温血马聚为一类，说明该种群最有可能主要来源于温血马和纯血马。这些结果揭示，康西草原地区马匹具有较高的遗传多样性，该地马群体的血统与纯血马和温血马相近，这为当地马遗传资源评估以及开发利用打下了坚实的基础。     

The legacy of Columbus in American horse populations assessed by microsatellite markers

Criollo horse populations descend from horses brought from the Iberian Peninsula over the period of colonization (15th to 17th century). They are spread throughout the Americas and have potentially undergone genetic hybridization with other breeds in the recent past. In this study, 25 autosomal microsatellites were genotyped in 50 horse breeds representing Criollo populations from 12 American countries (27 breeds), breeds from the Iberian Peninsula (19), one breed each from France and Morocco and two cosmopolitan horse breeds (Thoroughbred and Arabian). The genetic relationships among breeds identified five clusters: Celtic; Iberian; North American with Thoroughbred influence; most Colombian breeds; and nearly all other Criollo breeds. The group of “all other Criollo breeds” had the closest genetic relationship with breeds originating from the Iberian Peninsula, specifically with the Celtic group. For the whole set of Criollo breeds analysed, the estimated genetic contribution from other breeds was approximately 50%, 30% and 20% for the Celtic, Iberian and Arab-Thoroughbred groups, respectively. The spatial distribution of genetic diversity indicates that hotspots of genetic diversity are observed in populations from Colombia, Ecuador, Brazil, Paraguay and western United States, possibly indicating points of arrival and dispersion of Criollo horses in the American continent. These results indicate that Criollo breeds share a common ancestry, but that each breed has its own identity.     

利用微卫星遗传标记评估5个西藏牦牛群体遗传多样性及群体结构

【目的】 评估当前西藏主要牦牛群体遗传多样性，梳理不同地区群体间遗传结构，明确西藏5个牦牛群体（阿里牦牛、斯布牦牛、娘亚牦牛、类乌齐牦牛和帕里牦牛）的保种情况和种群间系统发育关系。【方法】 利用13个微卫星标记（SSR）对5个牦牛群体共计195个个体进行基因分型，并对各群体的等位基因数量、基因杂合度、多态信息含量及群体间遗传距离等遗传参数进行评估。【结果】 阿里牦牛群体等位基因数最多（6.43），类乌齐牦牛等位基因数最少（5.00）；观测杂合度范围为0.5311（娘亚牦牛）~0.5995（类乌齐牦牛）。各群体内偏离Hardy-Weinberg平衡的位点数为5（类乌齐牦牛）~9（阿里牦牛）个；群体内近交系数最高为0.172（阿里牦牛），且4个群体（阿里牦牛、娘亚牦牛、斯布牦牛和帕里牦牛）存在显著近交风险（P<0.05）。从遗传结构来看，所有群体间均为显著遗传分歧（P<0.05），STRUCTURE分析结果显示，5个牦牛群体划分为3个簇，其中阿里牦牛较其他牦牛群体具有更为丰富的遗传背景。系统发育树结果表明，不同群体间系统发育关系相对独立，且与种群栖息地分布不一致。主坐标分析（PCoA）结果显示，不同群体间部分个体存在较近亲缘关系，表明不同群体间存在遗传物质交流。【结论】 5个西藏牦牛群体具有丰富的遗传多样性水平，且种群系统发育关系相对独立，但多数群体存在群体事件风险。本研究不仅有助于明确西藏牦牛地方种群遗传多样性水平，同时可为今后的保种策略提升提供理论参考。