Endangered wild yak (Bos grunniens) in the Tibetan plateau and adjacent regions: Population size,distribution, conservation perspectives and its relation to the domestic subspecies

The yak (Bos grunniens) is a long-haired bovid, endemic to the Tibetan Plateau and the adjacent high-altitude regions. The domesticated subspecies of yak (B. grunniens grunniens) are abundant and closely associated with the livelihoods of herders, while the wild subspecies of yak (B. grunniens mutus) are endangered due primarily to anthropogenic effects. The endangered status of wild yaks calls for consideration, if we are to secure its long term survival, hence this study. Here we hope to provide baseline information necessary for further research and protection of the wild yak resources. We use published data to discuss their evolution, their characteristics as well as their distribution in the Tibetan Plateau and the adjacent high-altitude regions. We were able to come up with a world wild yak distribution map, which may be useful for establishing protected areas, as well as updating the species IUCN Red List Status. From the data available, we were also able to provide an estimate of the wild yak population in China (∼22,000 wild yaks living in China), corresponding to 90% of the total world population. We further discuss the major threats to yaks, and we give some suggestions for future and sustainable conservation.     

Structural Variants Selected during Yak Domestication Inferred from Long-Read Whole-Genome Sequencing

Structural variants (SVs) represent an important genetic resource for both natural and artificial selection. Here we present a chromosome-scale reference genome for domestic yak (Bos grunniens) that has longer contigs and scaffolds (N50 44.72 and 114.39 Mb, respectively) than reported for any other ruminant genome. We further obtained long-read resequencing data for 6 wild and 23 domestic yaks and constructed a genetic SV map of 372,220 SVs that covers the geographic range of the yaks. The majority of the SVs contains repetitive sequences and several are in or near genes. By comparing SVs in domestic and wild yaks, we identified genes that are predominantly related to the nervous system, behavior, immunity, and reproduction and may have been targeted by artificial selection during yak domestication. These findings provide new insights in the domestication of animals living at high altitude and highlight the importance of SVs in animal domestication.     

野牦牛和家牦牛粪便菌群与短链脂肪酸关系的研究

短链脂肪酸(SCFA)是反刍动物吸收饲草、饲料中营养物质的重要形式。肠道菌群能够降解食物生成SCFA并影响其比例。本文通过16S r DNA测序和气相色谱质谱联用仪,分别测定了野牦牛(Bos mutus)和家牦牛(Bos grunniens)粪便菌群组成及SCFA含量,通过比较分析两种牦牛肠道菌群与SCFA的关系,筛选出野牦牛肠道中与SCFA高浓度有正相关关系的菌群。结果显示,野牦牛粪便菌群主要有厚壁菌门(Firmicutes)(66. 47%)、拟杆菌门(Bacteroidetes)(26. 00%)和变形菌门(Proteobacteria)(3. 48%),主要的科有瘤胃球菌科(Ruminococcaceae)(55. 18%)、拟杆菌科(Bacteroidaceae)(8. 75%)和毛螺菌科(Lachnospiraceae)(7. 57%),家牦牛的菌群结构和组成与野牦牛相似。野牦牛粪便中SCFA以乙酸和丙酸为主,乙酸、丙酸、异丁酸、正丁酸和正戊酸的含量均显著高于家牦牛(P <0. 01)。Spearman相关分析显示,野牦牛粪便菌群中紫单胞菌科(Porphyromonadaceae)、拟杆菌科(Bacteroidaceae)、普雷沃氏菌科(Paraprevotellaceae)、理研菌科(Rikenellaceae)和韦荣球菌科(Erysipelotrichaceae)与SCFA具有较强相关性(r> 0. 4),而家牦牛仅有弱相关性(r <0. 3)。说明牦牛后肠道具有丰富的能够促进SCFA生成的益生菌群,进而提高食物的转化效率。     

Associations of single nucleotide polymorphisms in candidate genes with the polled trait in Datong domestic yaks

The domestic yak (Bos grunniens) is an iconic symbol of animal husbandry at high altitudes. Yaks exhibit unique external characteristics including long hair and large horns. However, hornless yaks can be found in different breeds and different populations. The hornless trait is also known as polled, and the POLL locus has been fine‐mapped to chromosome 1 in cattle (Bos taurus), although the underlying genetic basis of the polled trait is still unclear in the yak. Thus, we performed an association study to identify the genetic polymorphisms responsible for the polled trait in the yak. Fifty polled Datong domestic yaks and 51 horned individuals were selected randomly from a huge herd and were used as the case and control groups respectively for the association analysis. Twelve genes located in the candidate region of the POLL locus in cattle were used as references to detect DNA polymorphisms related to yak polledness, which were analyzed by sequencing and a high‐resolution melting test. We applied Fisher's exact test and haplotype analysis to show that a 147‐kb segment that included three protein‐coding genes C1H21orf62, GCFC1 and SYNJ1 was the most likely location of the POLL mutation in domestic yaks.     

No silver bullet? Snow leopard prey selection in Mt. Kangchenjunga,Nepal

In this study, we investigated the impact of domestic and wild prey availability on snow leopard prey preference in the Kangchenjunga Conservation Area of eastern Nepal—a region where small domestic livestock are absent and small wild ungulate prey are present. We took a comprehensive approach that combined fecal genetic sampling, macro‐ and microscopic analyses of snow leopard diets, and direct observation of blue sheep and livestock in the KCA. Out of the collected 88 putative snow leopard scat samples from 140 transects (290 km) in 27 (4 × 4 km²) sampling grid cells, 73 (83%) were confirmed to be from snow leopard. The genetic analysis accounted for 19 individual snow leopards (10 males and 9 females), with a mean population size estimate of 24 (95% CI: 19–29) and an average density of 3.9 snow leopards/100 km² within 609 km². The total available prey biomass of blue sheep and yak was estimated at 355,236 kg (505 kg yak/km² and 78 kg blue sheep/km²). From the available prey biomass, we estimated snow leopards consumed 7% annually, which comprised wild prey (49%), domestic livestock (45%), and 6% unidentified items. The estimated 47,736 kg blue sheep biomass gives a snow leopard‐to‐blue sheep ratio of 1:59 on a weight basis. The high preference of snow leopard to domestic livestock appears to be influenced by a much smaller available biomass of wild prey than in other regions of Nepal (e.g., 78 kg/km² in the KCA compared with a range of 200–300 kg/km² in other regions of Nepal). Along with livestock insurance scheme improvement, there needs to be a focus on improved livestock guarding, predator‐proof corrals as well as engaging and educating local people to be citizen scientists on the importance of snow leopard conservation, involving them in long‐term monitoring programs and promotion of ecotourism.     

野牦牛线粒体基因组序列测定及其系统进化

野牦牛属高寒地区的特有物种,是我国最珍贵的野生动物遗传资源之一,已被列为国家一级重点保护动物。对野牦牛mtDNA进行全序列测定和结构分析,并基于线粒体基因组序列对其系统发生进行了探讨。结果表明:(1)野牦牛线粒体基因组全序列的大小为16 322 bp,整个基因组由37个编码基因和D-loop区组成;22个tRNA基因序列长度为1 524 bp、2个RNA基因序列长度为2 528 bp、13个编码蛋白基因序列长度为11420 bp、D-loop区长度为892 bp。基因组中无间隔序列,基因间排列紧密,基因内无内含子。(2)野牦牛具有较丰富的遗传多样性。(3)分子系统发生关系显示牦牛为牛亚科中的一个独立属,即牦牛属(Poephagus),牦牛属包括家牦牛(Poephagus grunniens)和野牦牛(Poephagus mutus)2个种。野牦牛线粒体基因组全序列的获得和结构解析对研究牦牛的起源、演化和分类,以及野牦牛遗传资源的保护、开发和利用均具有重要的理论和实际意义。     

Accumulation of deleterious mutations in the domestic yak genome

Deleterious mutations play an important functional role, affecting trait phenotypes in ways that decrease the fitness of organisms. Estimating the frequency of occurrence and abundance has been a topic of much interest, especially in crops and livestock. The processes of domestication and breeding allow deleterious mutations to persist at high frequency, and identifying such deleterious mutations is particularly important for breed improvement. Here, we assessed genome‐wide patterns of deleterious variation in 59 domestic and 13 wild yaks using genome resequencing data. Based on the intersection of results given by three methods (provean , polyphen 2 and sift 4g ), we identified 3187 putative deleterious mutation sites affecting 2586 genes in domestic yaks and 2067 affecting 1701 genes in wild yaks. Multiple lines of evidence indicate a significant increase in the load of deleterious mutations in domesticated yaks compared to wild yaks. Private deleterious genes were found to be associated with the perception of smell and detection of chemical stimulus. We also identified 36 genes related to Mendelian genetic diseases involved in sensory perception, skeletal development and the nervous and immune systems. This study not only adds to the understanding of the genetic basis of yak domestication but also provides a rich catalog of variants that will facilitate future breeding‐related research on the yak genome and on other bovid species.     

Main roads and land cover shaped the genetic structure of a Mediterranean island wild boar population

Patterns of genetic differentiation within and among animal populations might vary due to the simple effect of distance or landscape features hindering gene flow. An assessment of how landscape connectivity affects gene flow can help guide management, especially in fragmented landscapes. Our objective was to analyze population genetic structure and landscape genetics of the native wild boar (Sus scrofa meridionalis) population inhabiting the island of Sardinia (Italy), and test for the existence of Isolation‐by‐Distance (IBD), Isolation‐by‐Barrier (IBB), and Isolation‐by‐Resistance (IBR). A total of 393 Sardinian wild boar samples were analyzed using a set of 16 microsatellite loci. Signals of genetic introgression from introduced non‐native wild boars or from domestic pigs were revealed by a Bayesian cluster analysis including 250 reference individuals belonging to European wild populations and domestic breeds. After removal of introgressed individuals, genetic structure in the population was investigated by different statistical approaches, supporting a partition into five discrete subpopulations, corresponding to five geographic areas on the island: north‐west (NW), central west (CW), south‐west (SW), north‐central east (NCE), and south‐east (SE). To test the IBD, IBB, and IBR hypotheses, we optimized resistance surfaces using genetic algorithms and linear mixed‐effects models with a maximum likelihood population effects parameterization. Landscape genetics analyses revealed that genetic discontinuities between subpopulations can be explained by landscape elements, suggesting that main roads, urban settings, and intensively cultivated areas are hampering gene flow (and thus individual movements) within the Sardinian wild boar population. Our results reveal how human‐transformed landscapes can affect genetic connectivity even in a large‐sized and highly mobile mammal such as the wild boar, and provide crucial information to manage the spread of pathogens, including the African Swine Fever virus, endemic in Sardinia.     

西藏牦牛mtDNA D-loop区的遗传多样性及其遗传分化

通过测定和分析西藏11个牦牛类群114个个体的mtDNA D-loop区全序列,对西藏牦牛的遗传多样性、类群间的亲缘关系及其遗传分化进行了研究。结果表明:①西藏牦牛mtDNA D-loop区全序列长度为890—896 bp,4种核苷酸T、C、A、G的平均比例分别为28.5%、25.3%、32.4%、13.8%,西藏牦牛mtDNA D-loop区富含碱基A+T,表现出一定的碱基偏好性。②共检测到130个变异位点,占分析总位点数的14.33%;其中单一多态位点85个,占多态位点总数的65.38%,简约信息位点45个,占多态位点总数的34.62%。序列变异中碱基缺失、插入和碱基替换等均有,其中碱基替换变异类型中转换114次,颠换12次,在转换变异类型中以A/G、T/C为主,占95.61%,在颠换变异类型中以A/T为主,占75%。③在114个个体中鉴定出90种单倍型,单倍型多样性为0.981±0.008,核苷酸多样性为0.01056±0.00701,均说明西藏牦牛具有丰富的单倍型类型。④90种单倍型分为2个聚类簇(Ⅰ、Ⅱ),聚类簇Ⅰ包含80种单倍型,占全部单倍型的88.89%,涵盖本研究中所有的西藏牦牛类群;聚类簇Ⅱ中有10种单倍型,占单倍型总数的11.11%,涉及的类群有工布江达、帕里、丁青、巴青、江达、类乌齐、桑桑、桑日、斯布,说明西藏牦牛可能有2个母系起源。⑤西藏牦牛类群间核苷酸分歧度(Dxy)在0.503%—1.416%之间,聚类分析和AMOVA分析显示西藏牦牛可分为两大类,康布牦牛、嘉黎牦牛为一类,其余的牦牛类群为另一类。     

The genetics of brown coat color and white spotting in domestic yaks (Bos grunniens)

Domestic yaks (Bos grunniens) exhibit two major coat color variations: a brown vs. wild‐type black pigmentation and a white spotting vs. wild‐type solid color pattern. The genetic basis for these variations in color and distribution remains largely unknown and may be complicated by a breeding history involving hybridization between yaks and cattle. Here, we investigated 92 domestic yaks from China using a candidate gene approach. Sequence variations in MC1R, PMEL and TYRP1 were surveyed in brown yaks; TYRP1 was unassociated with the coloration and excluded. Recessive mutations from MC1R, or p.Gln34*, p.Met73Leu and possibly p.Arg142Pro, are reported in bovids for the first time and accounted for approximately 40% of the brown yaks in this study. The remaining 60% of brown individuals correlated with a cattle‐derived deletion mutation from PMEL (p.Leu18del) in a dominant manner. Degrees of white spotting found in yaks vary from color sidedness and white face, to completely white. After examining the candidate gene KIT, we suggest that color‐sided and all‐white yaks are caused by the serial translations of KIT (Cs₆ or Cs₂₉) as reported for cattle. The white‐faced phenotype in yaks is associated with the KIT haplotype S^wf. All KIT mutations underlying the serial phenotypes of white spotting in yaks are identical to those in cattle, indicating that cattle are the likely source of white spotting in yaks. Our results reveal the complex genetic origins of domestic yak coat color as either native in yaks through evolution and domestication or as introduced from cattle through interspecific hybridization.     

贺兰山牦牛冬春季的生境选择

在2009年12月－2010年1月和2010年4－5月,采用样线法结合直接观察法对贺兰山牦牛的冬春季生境选择进行了研究。结果表明,牦牛冬季对11种生境因子有选择性,偏爱山地针叶林带,海拔小于2 000 m,优势乔木为灰榆,坡度小于10,下坡位,距水源距离大于1 200 m,人为干扰距离2 000－4 000 m,隐蔽级大于70 %;春季牦牛对13种利用生境生态因子有选择性,偏爱于亚高山灌丛和草甸带,海拔大于3 000 m,乔木密度小于1株,乔木高度小于3 m,乔木距离大于3 m,灌木密度大于4 0株,灌木距离小于1 m,植被盖度大于7 0 %,上坡位,距水源的距离小于8 00 m,人为干扰距离大于4 000 m,隐蔽级大于7 0 %。冬春季牦牛在海拔、植被类型、地形特征、优势乔木、灌木种类、坡位、坡向、人为干扰距离、距水源距离上存在显著差异。主成分分析表明冬季第一主成分的贡献率21.100 %,其中绝对值较大的相关系数是乔木距离、优势乔木、乔木高度和乔木密度;春季第一主成分的贡献率是31.247 %,其中绝对值较大的相关系数是植被类型、海拔高度、地形特征和灌木密度。与其他分布地区的牦牛相比,贺兰山地区的牦牛能适应当地的地理特征和气候环境。     

Domestication relaxed selective constraints on the yak mitochondrial genome

Wild yaks roaming in the high-altitude Qinghai-Tibetan plateau have to maintain high metabolic efficiency. However, domestic yaks do not require such high efficiency because of their limited activity. Hence, domestication may have caused the relaxation of selective constraints on the yak mitochondrial genome because mitochondrial mutations are extremely sensitive to energy-related selective pressures. We have tested this hypothesis by analyzing the mitochondrial genomes of 51 domestic yaks and 21 wild yaks. The results show that the ratio of nonsynonymous/synonymous substitutions in mitochondrial protein-coding genes is significantly higher in domestic yak lineages than those of wild yaks. This genetic difference suggests that the relaxation of selective constraints following the domestication in addition to bottlenecks has allowed faster accumulation of nonsilent substitutions in the yak mitochondrial genome, despite its short domestication history.     

Comprehensive investigation of nucleotide diverdity in yaks

To understand the maternal genetic diversity of Tianzhu white yak better, we analyzed mtDNA D‐loop sequences of 209 Tianzhu white yaks, which are from the central region of Tianzhu white yak habitat. Accordingly, a total of 45 haplotypes were identified in Tianzhu white yaks in this study, and 18 of them were unique. The nucleotide diversity and haplotype diversity of population studied were 0.024 ± 0.003 and 0.946 ± 0.007 respectively, revealing that Tianzhu white yak possess a relatively high genetic diversity. The phylogenetic analysis, combining D‐loop sequences in this study with 533 previous published D‐loop sequences of 13 yak breeds, indicated that Tianzhu white yaks fell mainly into haplogroup A and that a small portion belonged to haplogroups B, C, D and E. Moreover, six haplotypes of 20 individuals identified in Tianzhu white yak were in the taurine haplogroup, indicating hybridization between Bos taurus and Tianzhu white yaks. In summary, this study supplies a comprehensive maternal genetic pattern for Tianzhu white yak and provides a basic reference for future breeding programs to conserve the purebred Tianzhu white yak.     

Phylogeographical analyses of domestic and wild yaks based on mitochondrial DNA: new data and reappraisal

Aim We aimed to examine the phylogeographical structure and demographic history of domestic and wild yaks (Bos grunniens) based on a wide range of samples and complete mitochondrial genomic sequences. Location The Qinghai‐Tibetan Plateau (QTP) of western China. Methods All available D‐loop sequences for 405 domesticated yaks and 47 wild yaks were examined, including new sequences from 96 domestic and 34 wild yaks. We further sequenced the complete mitochondrial genomes of 48 domesticated and 21 wild yaks. Phylogeographical analyses were performed using the mitochondrial D‐loop and the total genome datasets. Results We recovered a total of 123 haplotypes based on the D‐loop sequences in wild and domestic yaks. Phylogenetic analyses of this dataset and the mitochondrial genome data suggested three well‐supported and divergent lineages. Two lineages with six D‐loop haplogroups were recovered for all morphological breeds of domestic yaks across their distributions in the QTP, while one more lineage and more endemic haplogroups or haplotypes were found for wild yaks. Based on the mitochondrial genome data, the divergences of the three lineages were estimated to have occurred around 420,000 and 580,000 years ago, consistent with the geological records of two large glaciation events experienced in the QTP. Main conclusions There are distinct phylogeographical differences between wild and domestic yaks. However, there is no apparent geographical correlation between identified haplogroups and distributions of domestic yaks. Three differentiated lineages of yaks probably evolved allopatrically in different regions during the Pleistocene glaciation events, then reunited into a single gene pool during post‐glacial population expansion and migrations before the start of the domestication of yaks in the Holocene.     

Assessment of cattle genetic introgression into domestic yak populations using mitochondrial and microsatellite DNA markers

Hybridization between yak Poephagus grunniens and taurine Bos taurus or indicine B. indicus cattle has been widely practiced throughout the yak geographical range, and gene flow is expected to have occurred between these species. To assess the impact of cattle admixture on domestic yak, we examined 1076 domestic yak from 29 populations collected in China, Bhutan, Nepal, India, Pakistan, Kyrgyzstan, Mongolia and Russia using mitochondrial DNA and 17 autosomal microsatellite loci. A cattle diagnostic marker‐based analysis reveals cattle‐specific mtDNA and/or autosomal microsatellite allele introgression in 127 yak individuals from 22 populations. The mean level of cattle admixture across the populations, calculated using allelic information at 17 autosomal microsatellite loci, remains relatively low (mY_cattle = 2.66 ± 0.53% and Q_cattle = 0.69 ± 2.58%), although it varies a lot across populations as well as among individuals within population. Although the level of cattle admixture shows a clear geographical structure, with higher levels of admixture in the Qinghai‐Tibetan Plateau and Mongolian and Russian regions, and lower levels in the Himalayan and Pamir Plateau region, our results indicate that the level of cattle admixture is not significantly correlated with the altitude across geographical regions as well as within geographical region. Although yak‐cattle hybridization is primarily driven to produce F₁ hybrids, our results show that the subsequent gene flow between yak and cattle took place and has affected contemporary genetic make‐up of domestic yak. To protect yak genetic integrity, hybridization between yak and cattle should be tightly controlled.     

Genome‐wide variation within and between wild and domestic yak

The yak is one of the few animals that can thrive in the harsh environment of the Qinghai‐Tibetan Plateau and adjacent Alpine regions. Yak provides essential resources allowing Tibetans to live at high altitudes. However, genetic variation within and between wild and domestic yak remain unknown. Here, we present a genome‐wide study of the genetic variation within and between wild and domestic yak. Using next‐generation sequencing technology, we resequenced three wild and three domestic yak with a mean of fivefold coverage using our published domestic yak genome as a reference. We identified a total of 8.38 million SNPs (7.14 million novel), 383 241 InDels and 126 352 structural variants between the six yak. We observed higher linkage disequilibrium in domestic yak than in wild yak and a modest but distinct genetic divergence between these two groups. We further identified more than a thousand of potential selected regions (PSRs) for the three domestic yak by scanning the whole genome. These genomic resources can be further used to study genetic diversity and select superior breeds of yak and other bovid species.     

Prevalence,Associated Risk Factors,and Phylogenetic Analysis of Toxocara vitulorum Infection in Yaks on the Qinghai Tibetan Plateau,China

Toxocara vitulorum has been rarely reported in yaks at high altitudes and remote areas of Sichuan Province of Tibetan Plateau of China. The current study was designed to investigate the prevalence, associated risk factors, and phylogenetic characteristics of T. vitulorum in yak calves on the Qinghai Tibetan plateau. Fecal samples were collected from 891 yak calves and were examined for the presence of T. vitulorum eggs by the McMaster technique. A multivariable logistic regression model was employed to explore variables potentially associated with exposure to T. vitulorum infection. T. vitulorum specimens were collected from the feces of yaks in Hongyuan of Sichuan Province, China. DNA was extracted from ascaris. After PCR amplification, the sequencing of ND1 gene was carried out and phylogenetic analyses was performed by MEGA 6.0 software. The results showed that 64 (20.1%; 95% CI 15.8–24.9%), 75 (17.2; 13.8–21.1), 29 (40.9; 29.3–53.2), and 5 (7.6; 2.5–16.8) yak calves were detected out to excrete T. vitulorum eggs in yak calve feces in Qinghai, Tibet, Sichuan, and Gansu, respectively. The present study revealed that high infection and mortality by T. vitulorum is wildly spread on the Qinghai Tibetan plateau, China by fecal examination. Geographical origin, ages, and fecal consistencies are the risk factors associated with T. vitulorum prevalence by logistic regression analysis. Molecular detection and phylogenetic analysis of ND1 gene of T. vitulorum indicated that T. vitulorum in the yak calves on the Qinghai Tibetan plateau are homologous to preveiously studies reported.     

The molecular phylogenetic signature of Bali cattle revealed by maternal and paternal markers

Bali cattle is a domestic cattle breed that can be found in Malaysia. It is a domestic cattle that was purely derived from a domestication event in Banteng (Bos javanicus) around 3,500 BC in Indonesia. This research was conducted to portray the phylogenetic relationships of the Bali cattle with other cattle species in Malaysia based on maternal and paternal lineage. We analyzed the cytochrome c oxidase I (COI) mitochondrial gene and SRY of Y chromosome obtained from five species of the Bos genus (B. javanicus, Bos gaurus, Bos indicus, Bos taurus, and Bos grunniens). The water buffalo (Bubalus bubalis) was used as an outgroup. The phylogenetic relationships were observed by employing several algorithms: Neighbor-Joining (PAUP version 4.0), Maximum parsimony (PAUP version 4.0) and Bayesian inference (MrBayes 3.1). Results from the maternal data showed that the Bali cattle formed a monophyletic clade, and together with the B. gaurus clade formed a wild cattle clade. Results were supported by high bootstrap and posterior probability values together with genetic distance data. For the paternal lineage, the sequence variation is low (with parsimony informative characters: 2/660) resulting an unresolved Neighbor-Joining tree. However, Bali cattle and other domestic cattle appear in two monophyletic clades distinct from yak, gaur and selembu. This study expresses the potential of the COI gene in portraying the phylogenetic relationships between several Bos species which is important for conservation efforts especially in decision making since cattle is highly bred and hybrid breeds are often formed. Genetic conservation for this high quality beef cattle breed is important by maintaining its genetic characters to prevent extinction or even decreased the genetic quality.     

牦牛的分类学地位及起源研究:mtDNA D-loop序列的分析

牦牛的起源与属级分类学地位至今仍然存在一定的争议。我们测定了家养牦牛和野生牦牛线粒体控制区(D-loop)序列,并以此构建牦牛和牛属、野牛属、水牛属以及非洲水牛属相关种的系统发育树。研究结果表明线粒体D-loop区与Cytb基因序列在构建牛族的系统发育具有同样重要的价值。系统发育关系显示野牛属的灭绝种草原野牛与现存种美洲野牛先聚合为一单系群,然后再和牦牛形成一单系分支,表明牦牛与野牛属的草原野牛、美洲野牛亲缘关系最近,具有最近的共同祖先,而与牛属的其它亚洲物种亲缘关系较远。因此,本研究不支持将牦牛独立为牦牛属———Poephagus,牛属与野牛属在分类上也应合并为一个属。基于上述研究结果和化石证据,我们进一步对牦牛起源的历史背景进行了讨论,认为牦牛与野牛属的分化是由于第四纪气候变化在欧亚大陆发生的,野牛通过白令陆桥进入北美;冰期结束后,由于欧亚大陆其它地区温度升高,牦牛只能局限分布在较为寒冷的青藏高原;而野牛属在北美先后分化为草原野牛和美洲野牛,前者可能是后者的直接祖先。     

Genetic Diversity and Divergence in Chinese Yak (Bos grunniens) Populations Inferred from Blood Protein Electrophoresis

In 6 Chinese yak (Bos. Grunniens) populations including 177 yaks, 34 blood protein loci were studied by horizontal starch gel electrophoresis, four of these loci (AKP, ALB, LDH-1, TF) were found to be polymorphic. The percentage of polymorphic loci(P) is 0.118, the mean individual heterozygosity(H) is 0.015, which means a low level of genetic diversity in the whole Chinese yak population. The coefficient of gene differentiation (G _ST ) is 0.0625, which indicated an almost-indistinguishable divergence among different populations at the level of blood protein electrophoresis.