Evolution and domestication of the Bovini species

Domestication of the Bovini species (taurine cattle, zebu, yak, river buffalo and swamp buffalo) since the early Holocene (ca. 10 000 BCE) has contributed significantly to the development of human civilization. In this study, we review recent literature on the origin and phylogeny, domestication and dispersal of the three major Bos species – taurine cattle, zebu and yak – and their genetic interactions. The global dispersion of taurine and zebu cattle was accompanied by population bottlenecks, which resulted in a marked phylogeographic differentiation of the mitochondrial and Y-chromosomal DNA. The high diversity of European breeds has been shaped through isolation-by-distance, different production objectives, breed formation and the expansion of popular breeds. The overlapping and broad ranges of taurine and zebu cattle led to hybridization with each other and with other bovine species. For instance, Chinese gayal carries zebu mitochondrial DNA; several Indonesian zebu descend from zebu bull × banteng cow crossings; Tibetan cattle and yak have exchanged gene variants; and about 5% of the American bison contain taurine mtDNA. Analysis at the genomic level indicates that introgression may have played a role in environmental adaptation.     

牛科动物HSL基因序列分析及其分子进化研究

在对牛科中4种动物即牦牛、瘤牛、普通牛和水牛HSL基因外显子Ⅰ部分核苷酸序列进行测定的基础上,与Gen-Bank中其他物种相应基因核苷酸序列、氨基酸序列进行了比对分析,并构建了牦牛与其他物种间分子系统进化树。结果表明:牦牛与普通牛、瘤牛、水牛、猪、人、小鼠、大鼠7个物种HSL基因外显子Ⅰ部分核苷酸序列间保守性较高,同源性大小依次为99.8%、99.6%、97.4%、90.6%、88.4%、83.5%、82.3%。相应氨基酸序列间保守性更高,同源性分别为100%、100%、98.2%、94.0%、92.2%、89.8%、89.8%。牦牛与各物种该基因部分核苷酸序列间碱基变异类型主要表现为碱基转换和颠换,无碱基插入和缺失发生,碱基转换的频率高于颠换的频率;在核苷酸水平上的多数碱基替换都是同义替换;序列间单碱基变异位点大多出现在同一位点,多发生在密码子第3位,其次是第1位,最少发生在第2位,符合分子进化的中性学说。HSL基因外显子Ⅰ部分核苷酸序列进行多序列对位排列构建的各物种间分子系统进化树结果表明,普通牛和瘤牛首先聚为一类,再分别与牦牛、水牛、猪、人聚类,最后与大鼠、小鼠聚为一类。该聚类结果与动物学上的分类结果一致,表明HSL基因外显子Ⅰ部分核苷酸序列适合于构建物种间分子系统进化树。研究表明,牦牛、普通牛和瘤牛3个物种间的遗传距离大小相近,牦牛和水牛间的遗传距离与普通牛、瘤牛和水牛间的遗传距离大小相当。牦牛、普通牛和瘤牛3个物种间的遗传距离远小于它们各自与水牛这一物种的遗传距离,它们三者之间的亲缘关系也相对于它们各自与水牛间的亲缘关系都较近,故将牦牛、普通牛和瘤牛划分在同一个属——牛属(Bos)更为合理。     

Complete sequence of the yak (Bos grunniens) mitochondrial genome and its evolutionary relationship with other ruminants

The yak (Bos grunniens) is the most important domesticated species in the Qinhai-Tibetan Plateau. In present study, the complete sequence of the yak mitochondrial genome was determined. Sequence analysis revealed that there are no differences with cattle in the yak mitochondrial genome organization. Interestingly, within the D-loop, the conserved sequence blocks are less conserved than surrounding regions. Neighbor-Joining (NJ) trees based on single genes, gene sets and concatenated genes of mitochondrial genome were constructed. The analysis identified the yak as a sister group of a cattle/zebu clade. Based on substitutions in 22 tRNA genes, 12S rRNA gene and 16S rRNA gene, the dating of divergence between yak and cattle/zebu, and yak and water buffalo, was proposed to have occurred 4.38-5.32 and 10.54-13.85 million years before present, respectively. This is consistent with the paleontologyical data. Yak and sheep/goat divergent dating predicts that their divergence occurred at 13.14-27.99 million years before the present day.     

Mitochondrial DNA variation in cattle of South China: origin and introgression

Ten restriction endonucleases were used to investigate the mitochondrial DNA restriction fragment length polymorphism (mtDNA RFLP) of 11 native cattle breeds and one cultivated cattle breed in South China. Twenty-three restriction morphs were detected, which can be sorted into five haplotypes. A phylogenetic tree of the haplotypes was constructed by using the 'upgMa' method. Our study showed that haplotype I and II are identical to the zebu (Bos indicus) and taurine (Bos taurus) haplotypes, respectively. Zebu and taurine were the two major origins of cattle populations in South China, and the zebu probably had more influence on the native cattle population than taurine did. Haplotype III is identical to haplotype I of yak (Bos grunniens), which was only detected in the Diqing cattle breed. Haplotype IV was detected for the first time. This haplotype, found only in Dehong cattle, might be from an independent domestication event, probably from another Bos indicus population. Divergence of haplotypes I and IV occurred about 268,000-535,000 years ago, much earlier than the 10,000-year history of cattle husbandry. Our results also suggest a secondary introgession of mtDNA from yak to Diqing cattle.     

中国部分地方黄牛品种线粒体D-loop区的遗传变异与进化分析

测定了13个黄牛品种125个个体的线粒体D-loop区段的全序列,包括12个中国地方黄牛品种的123个个体和德国黄牛2个个体,并进行了分析。结果显示,共检测到93个变异位点,57个单倍型,平均核苷酸差异(average number ofnucleotide differences,k)为22.708,核苷酸多样度(nucleotide diversity,π)为0.0251±0.00479,单倍型多样度(haplotypediversity,Hd)为0.888±0.026,表明我国黄牛品种遗传多样性非常丰富。构建的Neighbor-Joining进化树显示这13个品种主要分成两大类型:普通牛和瘤牛;新发现的特殊类型Ⅲ只有一个西藏阿沛甲咂牛的个体,它与牦牛D-loop序列最相近,证明西藏地区的黄牛与牦牛之间存在基因渗入现象。普通牛和瘤牛在日喀则驼峰牛中占的比例分别是64.3%和35.7%,在阿沛甲咂牛中占的比例分别是50.0%和50.0%,证明了西藏的黄牛也有瘤牛类型。云南牛品种的单倍型非常丰富证明了云南在中国黄牛起源上的重要地位;在27个中国黄牛品种中(本研究11个品种以及GenBank上的16个品种)找到了中国瘤牛的核心单倍型i1,并且对它进行了讨论。同时证明了西藏瘤牛独立于中国瘤牛核心类群的特殊性。     

Genetic Structure and Differentiation of Three Chinese Indigenous Cattle Populations

Levels of genetic differentiation, gene flow, and genetic structure of three indigenous cattle populations (Luxi, Bohai, and Minnan) and two reference cattle populations (Chinese Holstein and Qinhai yak) in China were estimated using the information from 12 microsatellites, and 141 microsatellite alleles were identified. The mean number of alleles per locus ranged from 2.9005 in yak to 4.9722 in Holstein. The observed heterozygosity ranged from 0.5325 (yak) to 0.7719 (Holstein); 29 private alleles were detected. The global heterozygote deficit across all populations amounted to 58.5% (p < 0.001). The overall significant (p < 0.001) deficit of heterozygotes because of inbreeding within breeds amounted to 43.2%. The five cattle populations were highly differentiated (F _st = 26.9%, p < 0.001) at all loci. The heterozygote deficit within the population was highest in Luxi cattle and lowest in yak. The average number of effective migrants exchanged per generation was highest (1.149) between Luxi and Holstein, and lowest (0.509) between Luxi and yak. With the application of prior population information, cluster analysis achieved posterior probabilities from 91% to 98% of correctly assigning individuals to populations. Combining the information of cluster analysis, gene flow, and Structure analysis, the five cattle populations belong to three genetic clusters, a taurine (Luxi and Chinese Holstein), a zebu (Bohai and Minnan), and a yak cluster. This indicates that Bohai black is closer to Bos indicus than Luxi cattle. The evolution and development of three indigenous cattle populations are discussed.     

A microsatellite survey of cattle from a centre of origin: the Near East

Eight humpless cattle breeds from the Near East, three from Europe, one from West Africa and two zebu breeds from India were screened with 20 microsatellite loci. Breeds from the Near East revealed considerable levels of introgression from zebu cattle, which was apparent most in populations from the East and which declined in populations further West. This nonrandom pattern is suggestive of the introduction of zebu cattle from the East. Notwithstanding the overlay of zebu alleles, it was possible to demonstrate that Near Eastern cattle exhibited significantly higher levels of allelic diversity than breeds from other regions, which is consistent with the view that this region represents a primary domestication centre for Bos taurus cattle. The hypothesis that B. taurus and B. indicus cattle have separate domestic origins is also supported by the survey, a large genetic divergence being apparent between the nonhybrid taurine and zebu groups.     

Diversity and phylogeny of mitochondrial DNA isolated from mithun Bos frontalis located in Bhutan

We sequenced the 16S rRNA gene in mitochondrial DNA to characterize mithun located in Bhutan and to increase our understanding of its origin. We compared mithun with yak, European cattle, Bhutanese zebu and Indian zebu. Sequencing revealed low nucleotide diversity within the mithun population and their phylogenetic proximity to gaur. A close relationship between Bhutanese mithun and gaur was confirmed by an additional comparison with wild gaur specimens from three locations in Bhutan. Direct domestication of mithun from gaur was supported, while maternal contribution from the cattle lineage during domestication was not supported.     

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.     

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

野牦牛属高寒地区的特有物种,是我国最珍贵的野生动物遗传资源之一,已被列为国家一级重点保护动物。对野牦牛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个种。野牦牛线粒体基因组全序列的获得和结构解析对研究牦牛的起源、演化和分类,以及野牦牛遗传资源的保护、开发和利用均具有重要的理论和实际意义。     

A new insight into cattle's maternal origin in six Asian countries

The domestication of cattle fuelled the development of agricultural society in the history of human being.The evolution and genetic relationship of cattle can be elucidated by investigating the variation of mitochondrial DNA (mtDNA) D-loop sequence.In this study,we built a cattle phylogeny with a pool of 856 individual D-loop sequences,of which 264 Chinese cattle D-loop sequences were obtained in this study (141 ones were first analyzed,and 123 were first submitted) and the rest sequences of cattle from six Asian countries (Japan,Korea,Mongolia,Nepal,India and China) were retrieved from GenBank.Our results indicated that cattle from six Asian countries fell into three clades,Bos taurus (taurine),Bos indicus (zebu) and yak.Four main haplogroups T1A,T2,T3 (including T3A and T3B) and T5 were found in taurine,and two haplogroups I1 and I2 in zebu.Furthermore,we found that I1 and I2 haplogroups were separated by four variable sites rather than five ones and four haplogroups or sub-haplogroups of T1A,T3A,T3B and T5 were found for the first time in these Asian cattle.These data brought us a new insight into cattle's genetic structure in these six Asian countries.The geographical distribution of haplogroups was also outlined to provide systematic information on cattle genetic resources.     

MspI Allelic Pattern of Bovine Growth Hormone Gene in Indian Zebu Cattle (Bos indicus) Breeds

The MspI allelic variation in intron III of the bovine growth hormone (bGH) gene was explored using PCR-RFLP in 750 animals belonging to 17 well-recognized breeds of Indian zebu cattle (Bos indicus) reared in different geographic locations of the country. Restriction digestion analysis of a 329-bp PCR fragment of the bGH intron III region with MspI restriction enzyme revealed two alleles (MspI− and MspI+) and two genotypes (−/− and +/−) across the 17 cattle breeds studied. The allelic frequency varied from 0.67 to 0.94 for MspI (−) and from 0.06 to 0.33 for MspI (+) across the 17 breeds, with a combined average frequency of 0.87 and 0.13, respectively. No animal with +/+ genotype was detected across the samples analyzed. The chi-square test showed that the difference in MspI allelic frequency was not significant (p > 0.05), regardless of the geographic origin, coat color, or utility of the cattle breed. The high MspI (−) allele frequencies obtained for Indian zebu cattle in this study are in sharp contrast to those reported for taurine breeds from northern Europe, Mediterranean countries, and America. Findings of this study further substantiate the hypothesis that the MspI (−) allele has an Indian origin.     

Correlating Bayesian date estimates with climatic events and domestication using a bovine case study

The tribe Bovini contains a number of commercially and culturally important species, such as cattle. Understanding their evolutionary time scale is important for distinguishing between post-glacial and domestication-associated population expansions, but estimates of bovine divergence times have been hindered by a lack of reliable calibration points. We present a Bayesian phylogenetic analysis of 481 mitochondrial D-loop sequences, including 228 radiocarbon-dated ancient DNA sequences, using a multi-demographic coalescent model. By employing the radiocarbon dates as internal calibrations, we co-estimate the bovine phylogeny and divergence times in a relaxed-clock framework. The analysis yields evidence for significant population expansions in both taurine and zebu cattle, European aurochs and yak clades. The divergence age estimates support domestication-associated expansion times (less than 12 kyr) for the major haplogroups of cattle. We compare the molecular and palaeontological estimates for the Bison-Bos divergence.     

Genetic diversity and admixture of indigenous cattle from North Ethiopia: implications of historical introgressions in the gateway region to Africa

Microsatellite variation was surveyed to determine the genetic diversity, population structure and admixture of seven North Ethiopian cattle breeds by combining multiple microsatellite data sets of Indian and West African zebu, and European, African and Near-Eastern taurine in genetic analyses. Based on allelic distribution, we identified four diagnostic alleles (HEL1-123 bp, CSSM66-201 bp, BM2113-150 bp and ILSTS6-285 bp) specific to the Near-Eastern taurine. Results of genetic relationship and population structure analyses confirmed the previously established marked genetic distinction between taurine and zebu, and indicated further divergence among the bio-geographical groupings of breeds such as North Ethiopian, Indian and West African zebu, and African, European and Near-Eastern taurine. Using the diagnostic alleles for bio-geographical groupings and a Bayesian method for population structure inference, we estimated the genetic influences of major historical introgressions in North Ethiopian cattle. The breeds have been heavily (>90%) influenced by zebu, followed by African, European and the Near-Eastern taurine. Overall, North Ethiopian cattle show a high level of within-population genetic variation (e.g. observed heterozygosity = 0.659-0.687), which is in the upper range of that reported for domestic cattle and indicates their potential for future breeding applications, even in a global context. Rather low but significant population differentiation (F(ST) = 1.1%, P < 0.05) was recorded as a result of multiple introgression events and strong genetic exchanges among the North Ethiopian breeds.     

Phylogeny of bovine species based on AFLP fingerprinting

The Bovini species comprise both domestic and wild cattle species. Published phylogenies of this tribe based on mitochondrial DNA contain anomalies, while nuclear sequences show only low variation. We have used amplified fragment length polymorphism (AFLP) fingerprinting in order to detect variation in loci distributed over the nuclear genome. Computer-assisted scoring of electrophoretic fingerprinting patterns yielded 361 markers, which provided sufficient redundancy to suppress stochastic effects of intraspecies polymorphisms and length homoplasies (comigration of non-homologous fragments). Tree reconstructions reveal three clusters: African buffalo with water buffalo, ox with zebu, and bison with wisent. Similarity values suggest a clustering of gaur and banteng, but bifurcating clustering algorithms did not assign consistent positions to these species and yak. We propose that because of shared polymorphisms and reticulations, tree topologies are only partially adequate to represent the phylogeny of the Bovini. Principal-coordinate analysis positions zebu between a gaur/banteng cluster and taurine cattle. This correlates with the region of origin of these species and suggests that genomic distances between the cattle species have been influenced by genetic exchange between neighbouring ancestral populations.     

Genetic diversity,introgression and relationships among West/Central African cattle breeds

Genetic diversity, introgression and relationships were studied in 521 individuals from 9 African Bos indicus and 3 Bos taurus cattle breeds in Cameroon and Nigeria using genotype information on 28 markers (16 microsatellite, 7 milk protein and 5 blood protein markers). The genotypes of 13 of the 16 microsatellite markers studied on three European (German Angus, German Simmental and German Yellow) and two Indian (Nelore and Ongole) breeds were used to assess the relationships between them and the African breeds. Diversity levels at microsatellite loci were higher in the zebu than in the taurine breeds and were generally similar for protein loci in the breeds in each group. Microsatellite allelic distribution displayed groups of alleles specific to the Indian zebu, African taurine and European taurine. The level of the Indian zebu genetic admixture proportions in the African zebus was higher than the African taurine and European taurine admixture proportions, and ranged from 58.1% to 74.0%. The African taurine breed, Muturu was free of Indian zebu genes while its counter Namchi was highly introgressed (30.2%). Phylogenic reconstruction and principal component analysis indicate close relationships among the zebu breeds in Cameroon and Nigeria and a large genetic divergence between the main cattle groups – African taurine, European taurine and Indian zebu, and a central position for the African zebus. The study presents the first comprehensive information on the hybrid composition of the individual cattle breeds of Cameroon and Nigeria and the genetic relationships existing among them and other breeds outside of Africa. Strong evidence supporting separate domestication events for the Bos species is also provided.     

雷琼牛GH基因第五外显子遗传多样性及其分化

对16头雷琼牛GH基因第5外显子序列进行分析,发现了1个变异位点,定义了2种单倍型。引用巴州牦牛2个个体GH基因同源区序列并结合GenBank中牛属普通牛、其它瘤牛和牦牛3个种群与水牛1个远缘种GH基因同源区序列,分别采用邻接（NJ）法和最大简约（MP）法构建分子系统发育树,得到基本一致的拓扑结构,结果显示GH基因的分化早于雷琼牛（瘤牛）、其它瘤牛、普通牛、牦牛和水牛的分化,瘤牛物种内存在多型,同时证实了GH基因第5外显子区有着较高的突变率。     

Assessing introgression of Sahelian zebu genes into native Bos taurus breeds in Burkina Faso

A total of 350 samples were analyzed to estimate zebu gene proportions into two different taurine cattle breeds of Burkina Faso (Lobi and N’Dama) using 38 microsatellites and various statistical methodologies. West African and East African zebu samples were sequentially used as reference parental populations. Furthermore, N’Dama cattle from Congo, the composite South African Bonsmara cattle breed and a pool of European cattle were used successively as second parental populations. Independently of the methodology applied: (a) the use of West African zebu samples gave higher admixture coefficients than the East African zebu; (b) the higher zebu proportions were estimated when the European cattle was used as parental population 2; and (c) the use of the N’Dama population from Congo as parental population 2 gave the more consistent zebu proportion estimates for both the Lobi and the N’Dama breeds. In any case, the zebu admixture proportions estimated were not negligible and were always higher in the N’Dama cattle than in the Lobi cattle of Burkina Faso. This suggested that the introgression of Sahelian zebu genes into the taurine cattle of Southern West Africa can follow a complex pattern that can depend on local agro-ecological features. The current research pointed out that the estimation of admixture coefficients is highly dependent on both the assumptions underlying the methodologies applied and the selection of parental populations. Our analyses suggest that either too high or nil genetic identity between the parental and the expectedly derived populations must be avoided.     

Genetic characterization of Argentine and Bolivian Creole cattle breeds assessed through microsatellites

In the present report, the polymorphisms from 9 microsatellites were used to assess genetic diversity and relationships in 4 Creole cattle breeds from Argentina and Bolivia, 4 European taurine breeds, and 2 American zebu populations. The Creole populations display a relatively high level of genetic variation as estimated by allelic diversity and heterozygosity, whereas the British breeds displayed reduced levels of genetic diversity. The analysis of molecular variance indicated that 7.8% of variance can be explained by differences among taurine and zebu breeds. Consistent with these results, the first principal component (PC), which comprised the 40% of the total variance, clearly distinguishes these 2 groups. In addition, all constructed phylogenetic trees cluster together Nelore and Brahman breeds with robust bootstrap values. Only 1% of variance was due to difference between American Creole and European taurine cattle. Although this secondary split was supported by the classical genetic distance and the second PC (15%), the topology of trees is not particularly robust. The presence of zebu-specific alleles in Creole cattle allowed estimating a moderate degree of zebu admixture. When these data were compared with mitochondrial and Y chromosomal studies, a clear pattern of male-mediated introgression was revealed. The results presented here contribute to the understanding of origin and history of the American Creole cattle.     

Microsatellite DNA Variation and the Evolution, Domestication and Phylogeography of Taurine and Zebu Cattle (Bos Taurus and Bos Indicus)

Genetic variation at 20 microsatellite loci was surveyed to determine the evolutionary relationships and molecular biogeography of 20 different cattle populations from Africa, Europe and Asia. Phylogenetic reconstruction and multivariate analysis highlighted a marked distinction between humpless (taurine) and humped (zebu) cattle, providing strong support for a separate origin for domesticated zebu cattle. A molecular clock calculation using bison (Bison sp.) as an outgroup gave an estimated divergence time between the two subspecies of 610,000-850,000 years. Substantial differences in the distribution of alleles at 10 of these loci were observed between zebu and taurine cattle. These markers subsequently proved very useful for investigations of gene flow and admixture in African populations. When these data were considered in conjunction with previous mitochondrial and Y chromosomal studies, a distinctive male-mediated pattern of zebu genetic introgression was revealed. The introgression of zebu-specific alleles in African cattle afforded a high resolution perspective on the hybrid nature of African cattle populations and also suggested that certain West African populations of valuable disease-tolerant taurine cattle are under threat of genetic absorption by migrating zebu herds.