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.     

The Balkans and the colonization of Europe: the post‐glacial range expansion of the wild boar,Sus scrofa

Aim We focus on the biogeographical role of the Balkan Peninsula as a glacial refugium and source of northward post‐glacial dispersal for many European taxa. Specifically, we analysed the genetic structure and variation of wild boar (Sus scrofa) samples primarily from Greece, a region that has repeatedly served as a glacial refugium within the Balkan Peninsula. Location Continental Greece, the Aegean island of Samos and Bulgaria. Methods We analysed wild boar samples from 18 localities. Samples from common domestic breeds were also examined to take into account interactions between wild and domesticated animals. Phylogenetic analyses were carried out on a 637‐bp fragment of the mitochondrial DNA control region in 200 wild boar and 27 domestic pigs. The sequences were also compared with 791 Eurasian wild boar and domestic pig D‐loop sequences obtained from GenBank. Results Ninety‐four haplotypes were identified in the European wild boar data set, of which 68 were found in the Balkan samples and assigned to two previously described clades: the E1 European and Near Eastern clades. All of the continental samples clustered in the E1 clade and the samples from Samos fell into the Near Eastern clade, consistent with the island’s proximity to Asia Minor. Intriguingly, 62 novel haplotypes were identified and are found exclusively in the Balkans. Only six haplotypes were shared between wild boar and domestic pigs. Main conclusions Our data reveal numerous novel and geographically restricted haplotypes in wild boar populations, suggesting the presence of separate refugia in the Balkans. Our analyses support the hypothesis of a post‐glacial wild boar expansion consistent with the leading edge model, north and west from modern day Greece, and suggest little maternal introgression of Near Eastern and domestic haplotypes into wild Balkan populations.     

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.     

Y‐chromosome and mtDNA variation confirms independent domestications and directional hybridization in South American camelids

Investigations of genetic diversity and domestication in South American camelids (SAC) have relied on autosomal microsatellite and maternally‐inherited mitochondrial data. We present the first integrated analysis of domestic and wild SAC combining male and female sex‐specific markers (male specific Y‐chromosome and female‐specific mtDNA sequence variation) to assess: (i) hypotheses about the origin of domestic camelids, (ii) directionality of introgression among domestic and/or wild taxa as evidence of hybridization and (iii) currently recognized subspecies patterns. Three male‐specific Y‐chromosome markers and control region sequences of mitochondrial DNA are studied here. Although no sequence variation was found in SRY and ZFY, there were seven variable sites in DBY generating five haplotypes on the Y‐chromosome. The haplotype network showed clear separation between haplogroups of guanaco–llama and vicuña–alpaca, indicating two genetically distinct patrilineages with near absence of shared haplotypes between guanacos and vicuñas. Although we document some examples of directional hybridization, the patterns strongly support the hypothesis that llama (Lama glama) is derived from guanaco (Lama guanicoe) and the alpaca (Vicugna pacos) from vicuña (Vicugna vicugna). Within male guanacos we identified a haplogroup formed by three haplotypes with different geographical distributions, the northernmost of which (Peru and northern Chile) was also observed in llamas, supporting the commonly held hypothesis that llamas were domesticated from the northernmost populations of guanacos (L. g. cacilensis). Southern guanacos shared the other two haplotypes. A second haplogroup, consisting of two haplotypes, was mostly present in vicuñas and alpacas. However, Y‐chromosome variation did not distinguish the two subspecies of vicuñas.     

Maternal and paternal genetic diversity of ancient sheep in Estonia from the Late Bronze Age to the post‐medieval period and comparison with other regions in Eurasia

Sheep were among the first domesticated animals to appear in Estonia in the late Neolithic and became one of the most widespread livestock species in the region from the Late Bronze Age onwards. However, the origin and historical expansion of local sheep populations in Estonia remain poorly understood. Here, we analysed fragments of the hypervariable D‐loop of mitochondrial DNA (mtDNA; 213 bp) and the Y‐chromosome SRY gene (130 bp) extracted from 31 archaeological sheep bones dated from approximately 800 BC to 1700 AD. The ancient DNA data of sheep from Estonia were compared with ancient sheep from Finland as well as a set of contemporary sheep breeds from across Eurasia in order to place them in a wider phylogeographical context. The analysis shows that: (i) 24 successfully amplified and analysed mtDNA sequences of ancient sheep cluster into two haplogroups, A and B, of which B is predominant; (ii) four of the ancient mtDNA haplotypes are novel; (iii) higher mtDNA haplotype diversity occurred during the Middle Ages as compared to other periods, a fact concordant with the historical context of expanding international trade during the Middle Ages; (iv) the proportion of rarer haplotypes declined during the expansion of sheep from the Near Eastern domestication centre to the northern European region; (v) three male samples showed the presence of the characteristic northern European haplotype, SNP G‐oY1 of the Y‐chromosome, and represent the earliest occurrence of this haplotype. Our results provide the first insight into the genetic diversity and phylogeographical background of ancient sheep in Estonia and provide basis for further studies on the temporal fluctuations of ancient sheep populations.     

Mitochondrial DNA D‐loop sequence variation in maternal lineages of Iranian native horses

To understand the origin and genetic diversity of Iranian native horses, mitochondrial DNA (mtDNA) D‐loop sequences were generated for 95 horses from five breeds sampled in eight geographical locations in Iran. Sequence analysis of a 247‐bp segment revealed a total of 27 haplotypes with 38 polymorphic sites. Twelve of 19 mtDNA haplogroups were identified in the samples. The most common haplotypes were found within haplogroup X2. Within‐population haplotype and nucleotide diversities of the five breeds ranged from 0.838 ± 0.056 to 0.974 ± 0.022 and 0.011 ± 0.002 to 0.021 ± 0.001 respectively, indicating a relatively high genetic diversity in Iranian horses. The identification of several ancient sequences common between the breeds suggests that the lineage of the majority of Iranian horse breeds is old and obviously originated from a vast number of mares. We found in all native Iranian horse breeds lineages of the haplogroups D and K, which is concordant with the previous findings of Asian origins of these haplogroups. The presence of haplotypes E and K in our study also is consistent with a geographical west–east direction of increasing frequency of these haplotypes and a genetic fusion in Iranian horse breeds.     

Origin of mitochondrial DNA diversity of domestic yaks

Background  

The domestication of plants and animals was extremely important anthropologically. Previous studies have revealed a general tendency for populations of livestock species to include deeply divergent maternal lineages, indicating that they were domesticated in multiple, independent events from genetically discrete wild populations. However, in water buffalo, there are suggestions that a similar deep maternal bifurcation may have originated from a single population. These hypotheses have rarely been rigorously tested because of a lack of sufficient wild samples. To investigate the origin of the domestic yak (Poephagus grunnies), we analyzed 637 bp of maternal inherited mtDNA from 13 wild yaks (including eight wild yaks from a small population in west Qinghai) and 250 domesticated yaks from major herding regions.     

Five ovine mitochondrial lineages identified from sheep breeds of the near East

Archaeozoological evidence indicates that sheep were first domesticated in the Fertile Crescent. To search for DNA sequence diversity arising from previously undetected domestication events, this survey examined nine breeds of sheep from modern-day Turkey and Israel. A total of 2027 bp of mitochondrial DNA (mtDNA) sequence from 197 sheep revealed a total of 85 haplotypes and a high level of genetic diversity. Six individuals carried three haplotypes, which clustered separately from the known ovine mtDNA lineages A, B, and C. Analysis of genetic distance, mismatch distribution, and comparisons with wild sheep confirmed that these represent two additional mtDNA lineages denoted D and E. The two haplogroup E sequences were found to link the previously identified groups A and C. The single haplogroup D sequence branched with the eastern mouflon (Ovis orientalis), urial (O. vignei), and argali (O. ammon) sheep. High sequence diversity (K = 1.86%, haplogroup D and O. orientalis) indicates that the wild progenitor of this domestic lineage remains unresolved. The identification in this study of evidence for additional domestication events adds to the emerging view that sheep were recruited from wild populations multiple times in the same way as for other livestock species such as goat, cattle, and pig.     

Strong and stable geographic differentiation of swamp buffalo maternal and paternal lineages indicates domestication in the China/Indochina border region

The swamp type of the Asian water buffalo is assumed to have been domesticated by about 4000 years BP, following the introduction of rice cultivation. Previous localizations of the domestication site were based on mitochondrial DNA (mtDNA) variation within China, accounting only for the maternal lineage. We carried out a comprehensive sampling of China, Taiwan, Vietnam, Laos, Thailand, Nepal and Bangladesh and sequenced the mtDNA Cytochrome b gene and control region and the Y‐chromosomal ZFY, SRY and DBY sequences. Swamp buffalo has a higher diversity of both maternal and paternal lineages than river buffalo, with also a remarkable contrast between a weak phylogeographic structure of river buffalo and a strong geographic differentiation of swamp buffalo. The highest diversity of the swamp buffalo maternal lineages was found in south China and north Indochina on both banks of the Mekong River, while the highest diversity in paternal lineages was in the China/Indochina border region. We propose that domestication in this region was later followed by introgressive capture of wild cows west of the Mekong. Migration to the north followed the Yangtze valley as well as a more eastern route, but also involved translocations of both cows and bulls over large distances with a minor influence of river buffaloes in recent decades. Bayesian analyses of various migration models also supported domestication in the China/Indochina border region. Coalescence analysis yielded consistent estimates for the expansion of the major swamp buffalo haplogroups with a credibility interval of 900 to 3900 years BP. The spatial differentiation of mtDNA and Y‐chromosomal haplotype distributions indicates a lack of gene flow between established populations that is unprecedented in livestock.     

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.     

Baboon phylogeny as inferred from complete mitochondrial genomes

Baboons (genus Papio) are an interesting phylogeographical primate model for the evolution of savanna species during the Pleistocene. Earlier studies, based on partial mitochondrial sequence information, revealed seven major haplogroups indicating multiple para‐ and polyphylies among the six baboon species. The most basal splits among baboon lineages remained unresolved and the credibility intervals for divergence time estimates were rather large. Assuming that genetic variation within the two studied mitochondrial loci so far was insufficient to infer the apparently rapid early radiation of baboons we used complete mitochondrial sequence information of ten specimens, representing all major baboon lineages, to reconstruct a baboon phylogeny and to re‐estimate divergence times. Our data confirmed the earlier tree topology including the para‐ and polyphyletic relationships of most baboon species; divergence time estimates are slightly younger and credibility intervals narrowed substantially, thus making the estimates more precise. However, the most basal relationships could not be resolved and it remains open whether (1) the most southern population of baboons diverged first or (2) a major split occurred between southern and northern clades. Our study shows that complete mitochondrial genome sequences are more effective to reconstruct robust phylogenies and to narrow down estimated divergence time intervals than only short portions of the mitochondrial genome, although there are also limitations in resolving phylogenetic relationships. Am J Phys Anthropol, 2013. © 2012 Wiley Periodicals, Inc.     

Origin and post-glacial dispersal of mitochondrial DNA haplogroups C and D in northern Asia

More than a half of the northern Asian pool of human mitochondrial DNA (mtDNA) is fragmented into a number of subclades of haplogroups C and D, two of the most frequent haplogroups throughout northern, eastern, central Asia and America. While there has been considerable recent progress in studying mitochondrial variation in eastern Asia and America at the complete genome resolution, little comparable data is available for regions such as southern Siberia--the area where most of northern Asian haplogroups, including C and D, likely diversified. This gap in our knowledge causes a serious barrier for progress in understanding the demographic pre-history of northern Eurasia in general. Here we describe the phylogeography of haplogroups C and D in the populations of northern and eastern Asia. We have analyzed 770 samples from haplogroups C and D (174 and 596, respectively) at high resolution, including 182 novel complete mtDNA sequences representing haplogroups C and D (83 and 99, respectively). The present-day variation of haplogroups C and D suggests that these mtDNA clades expanded before the Last Glacial Maximum (LGM), with their oldest lineages being present in the eastern Asia. Unlike in eastern Asia, most of the northern Asian variants of haplogroups C and D began the expansion after the LGM, thus pointing to post-glacial re-colonization of northern Asia. Our results show that both haplogroups were involved in migrations, from eastern Asia and southern Siberia to eastern and northeastern Europe, likely during the middle Holocene.     

Multiple maternal origins of native modern and ancient horse populations in China

To obtain more knowledge of the origin and genetic diversity of domestic horses in China, this study provides a comprehensive analysis of mitochondrial DNA (mtDNA) D-loop sequence diversity from nine horse breeds in China in conjunction with ancient DNA data and evidence from archaeological and historical records. A 247-bp mitochondrial D-loop sequence from 182 modern samples revealed a total of 70 haplotypes with a high level of genetic diversity. Seven major mtDNA haplogroups (A–G) and 16 clusters were identified for the 182 Chinese modern horses. In the present study, nine 247-bp mitochondrial D-loop sequences of ancient remains of Bronze Age horse from the Chifeng region of Inner Mongolia in China ( c. 4000–2000a bp ) were used to explore the origin and diversity of Chinese modern horses and the phylogenetic relationship between ancient and modern horses. The nine ancient horses carried seven haplotypes with rich genetic diversity, which were clustered together with modern individuals among haplogroups A, E and F. Modern domestic horse and ancient horse data support the multiple origins of domestic horses in China. This study supports the argument that multiple successful events of horse domestication, including separate introductions of wild mares into the domestic herds, may have occurred in antiquity, and that China cannot be excluded from these events. Indeed, the association of Far Eastern mtDNA types to haplogroup F was highly significant using Fisher's exact test of independence ( P  = 0.00002), lending support for Chinese domestication of this haplogroup. High diversity and all seven mtDNA haplogroups (A–G) with 16 clusters also suggest that further work is necessary to shed more light on horse domestication in China.     

Diversity and evolution of rice progenitors in Australia

In the thousands of years of rice domestication in Asia, many useful genes have been lost from the gene pool. Wild rice is a key source of diversity for domesticated rice. Genome sequencing has suggested that the wild rice populations in northern Australia may include novel taxa, within the AA genome group of close (interfertile) wild relatives of domesticated rice that have evolved independently due to geographic separation and been isolated from the loss of diversity associated with gene flow from the large populations of domesticated rice in Asia. Australian wild rice was collected from 27 sites from Townsville to the northern tip of Cape York. Whole chloroplast genome sequences and 4,555 nuclear gene sequences (more than 8 Mbp) were used to explore genetic relationships between these populations and other wild and domesticated rices. Analysis of the chloroplast and nuclear data showed very clear evidence of distinctness from other AA genome Oryza species with significant divergence between Australian populations. Phylogenetic analysis suggested the Australian populations represent the earliest‐branching AA genome lineages and may be critical resources for global rice food security. Nuclear genome analysis demonstrated that the diverse O. meridionalis populations were sister to all other AA genome taxa while the Australian O. rufipogon‐like populations were associated with the clade that included domesticated rice. Populations of apparent hybrids between the taxa were also identified suggesting ongoing dynamic evolution of wild rice in Australia. These introgressions model events similar to those likely to have been involved in the domestication of rice.     

Genetic relationship amongst the major non-coding regions of mitochondrial DNAs in wild boars and several breeds of domesticated pigs

We completed phylogenetic analysis of the major non-coding region of the mitochondrial DNA (mtDNA) from 159 animals of eight Euro-American and six East Asian domesticated pig breeds and 164 Japanese and five European wild boars. A total of 62 mtDNA haplotypes were detected. Alignment of these regions revealed nucleotide variations (including gaps) at 73 positions, including 58 sites with transition nucleotide substitutions, and two transversion substitutions. Phylogenetic analysis of the sequences could not organize domestic pig breeds into discrete clusters. In addition, many of the haplotypes found in members of diverged clustering groups were found primarily in Euro-American pig breeds, indicating extensive introgression of Asian domestic pigs into European breeds. Furthermore, phylogenetic analysis allocated the DNA sequences of non-coding regions into two different groups, and the deepest branchpoint of this porcine phylogeny corresponded to 86 000-136 000 years before present. This time of divergence would predate the historical period when the pig is thought to have been domesticated from the wild boar.     

The Enigmatic Origin of Bovine mtDNA Haplogroup R: Sporadic Interbreeding or an Independent Event of Bos primigenius Domestication in Italy?

Background

When domestic taurine cattle diffused from the Fertile Crescent, local wild aurochsen (Bos primigenius) were still numerous. Moreover, aurochsen and introduced cattle often coexisted for millennia, thus providing potential conditions not only for spontaneous interbreeding, but also for pastoralists to create secondary domestication centers involving local aurochs populations. Recent mitochondrial genomes analyses revealed that not all modern taurine mtDNAs belong to the shallow macro-haplogroup T of Near Eastern origin, as demonstrated by the detection of three branches (P, Q and R) radiating prior to the T node in the bovine phylogeny. These uncommon haplogroups represent excellent tools to evaluate if sporadic interbreeding or even additional events of cattle domestication occurred.

Methodology

The survey of the mitochondrial DNA (mtDNA) control-region variation of 1,747 bovine samples (1,128 new and 619 from previous studies) belonging to 37 European breeds allowed the identification of 16 novel non-T mtDNAs, which after complete genome sequencing were confirmed as members of haplogroups Q and R. These mtDNAs were then integrated in a phylogenetic tree encompassing all available P, Q and R complete mtDNA sequences.

Conclusions

Phylogenetic analyses of 28 mitochondrial genomes belonging to haplogroups P (N = 2), Q (N = 16) and R (N = 10) together with an extensive survey of all previously published mtDNA datasets revealed major similarities between haplogroups Q and T. Therefore, Q most likely represents an additional minor lineage domesticated in the Near East together with the founders of the T subhaplogroups. Whereas, haplogroup R is found, at least for the moment, only in Italy and nowhere else, either in modern or ancient samples, thus supporting an origin from European aurochsen. Haplogroup R could have been acquired through sporadic interbreeding of wild and domestic animals, but our data do not rule out the possibility of a local and secondary event of B. primigenius domestication in Italy.     

Comprehensive analysis of the mitochondrial DNA diversity in Chinese cattle

Complete mitochondrial DNA D‐loop sequences of 1105 individuals were used to assess the diversity of maternal lineages of cattle populations in China. In total, 250 taurine and 88 zebu haplotypes were identified. Five main haplogroups—T1a, T2, T3, T4 and T5—were identified in Bos taurus, whereas Bos indicus harbored two haplogroups—I1 and I2. Our results suggest that the distribution of T1a in Asia was concentrated mainly in the northeast region (northeast China, Korea and Japan); haplogroups T2, T3 and T4 were predominant in Chinese cattle; and T5 was sporadically detected in Mongolian and Pingwu cattle. In contrast to the widespread presence of I1, I2 was distributed only in southwestern China (Yunnan‐Guizhou Plateau and the Tibet Autonomous Region) and Xinjiang Uygur Autonomous Region. This is the first time that all five taurine haplogroups and two zebu haplogroups have been found in Mongolian cattle. In addition, eight individuals in Tibetan cattle carried the Bos grunniens mtDNA type. The high mtDNA diversity (H = 0.904 ± 0.008) and the weak genetic structure among the 57 Chinese cattle breeds/populations are consistent with their complex historical background, migration route and ecological environment.     

Mitochondrial DNA variation of Nigerian domestic helmeted guinea fowl

We analyzed genetic diversity of 215 mitochondrial DNA (mtDNA) D‐loop sequences from seven populations of domesticated helmeted guinea fowl (Numida meleagris) in Nigeria and compared that with results of samples collected in Kenya (n = 4) and China (n = 22). In total, 241 sequences were assigned to 22 distinct haplotypes. Haplotype diversity in Nigeria was 0.693 ± 0.022. The network grouped most matrilines into two main haplogroups: A and B. There was an absence of a geographic signal, and two haplotypes dominated across all locations with the exception of the Kebbi population in the northwest of the country; AMOVA also confirmed this observation (F_ST = 0.035). The low genetic diversity may be a result of recent domestication, whereas the lack of maternal genetic structure likely suggests the extensive genetic intermixing within the country. Additionally, the differentiation of the Kebbi population may be due to a certain demographic history and/or artificial selection that shaped its haplotype profile. The current data do not permit us to make further conclusions; therefore, more research evidence from genetics and archaeology is still required.     

Mitochondrial genome primers for Lake Malawi cichlids

Resolving the evolutionary history of rapidly diversifying lineages like the Lake Malawi Cichlid Flock demands powerful phylogenetic tools. Although this clade of over 500 species of fish likely diversified in less than two million years, the availability of extensive sequence data sets, such as complete mitochondrial genomes, could help resolve evolutionary patterns in this group. Using a large number of newly developed primers, we generated whole mitochondrial genome sequences for 14 Lake Malawi cichlids. We compared sequence divergence across protein‐coding regions of the mitochondrial genome and also compared divergence in the mitochondrial loci to divergence at two nuclear protein‐coding loci, Mitfb and Dlx2. Despite the widespread sharing of haplotypes of identical sequences at individual loci, the combined use of all protein‐coding mitochondrial loci provided a bifurcating phylogenetic hypothesis for the exemplars of major lineages within the Lake Malawi cichlid radiation. The primers presented here could have substantial utility for evolutionary analyses of mitochondrial evolution and hybridization within this diverse clade.     

Independent mitochondrial origin and historical genetic differentiation in North Eastern Asian cattle

In order to clarify the origin and genetic diversity of cattle in North Eastern Asia, this study examined mitochondrial displacement loop sequence variation and frequencies of Bos taurus and Bos indicus Y chromosome haplotypes in Japanese, Mongolian, and Korean native cattle. In mitochondrial analyses, 20% of Mongolian cattle carried B. indicus mitochondrial haplotypes, but Japanese and Korean cattle carried only B. taurus haplotypes. In contrast, all samples revealed B. taurus Y chromosome haplotypes. This may be due to the import of zebu and other cattle during the Mongol Empire era with subsequent crossing with native taurine cattle. B. taurus mtDNA sequences fall into several geographically distributed haplogroups and one of these, termed here T4, is described in each of the test samples, but has not been observed in Near Eastern, European or African cattle. This may have been locally domesticated from an East Eurasian strain of Bos primigenius.