Cytochrome b sequences of ancient cattle and wild ox support phylogenetic complexity in the ancient and modern bovine populations

Mitochondrial DNA has been the traditional marker for the study of animal domestication, as its high mutation rate allows for the accumulation of molecular diversity within the time frame of domestic history. Additionally, it is exclusively maternally inherited and haplotypes become part of the domestic gene pool via actual capture of a female animal rather than by interbreeding with wild populations. Initial studies of British aurochs identified a haplogroup, designated P, which was found to be highly divergent from all known domestic haplotypes over the most variable portion of the D-loop. Additional analysis of a large and geographically representative sample of aurochs from northern and central Europe found an additional, separate aurochs haplotype, E. Until recently, the European aurochs appeared to have no matrilinear descendants among the publicly available modern cattle control regions sequenced; if aurochs mtDNA was incorporated into the domestic population, aurochs either formed a very small proportion of modern diversity or had been subsequently lost. However, a haplogroup P sequence has recently been found in a modern sample, along with a new divergent haplogroup called Q. Here we confirm the outlying status of the novel Q and E haplogroups and the modern P haplogroup sequence as a descendent of European aurochs, by retrieval and analysis of cytochrome b sequence data from twenty ancient wild and domesticated cattle archaeological samples.     

Y-SNPs do not indicate hybridisation between European aurochs and domestic cattle

Background

Previous genetic studies of modern and ancient mitochondrial DNA have confirmed the Near Eastern origin of early European domestic cattle. However, these studies were not able to test whether hybridisation with male aurochs occurred post-domestication. To address this issue, Götherström and colleagues (2005) investigated the frequencies of two Y-chromosomal haplotypes in extant bulls. They found a significant influence of wild aurochs males on domestic populations thus challenging the common view on early domestication and Neolithic stock-rearing. To test their hypothesis, we applied these Y-markers on Neolithic bone specimens from various European archaeological sites.

Methods and Findings

Here, we have analysed the ancient DNA of 59 Neolithic skeletal samples. After initial molecular sexing, two segregating Y-SNPs were identified in 13 bulls. Strikingly, our results do not support the hypothesis that these markers distinguish European aurochs from domesticated cattle.

Conclusions

The model of a rapid introduction of domestic cattle into Central Europe without significant crossbreeding with local wild cattle remains unchallenged.     

Complete mitochondrial genomes of Bos taurus and Bos indicus provide new insights into intra-species variation, taxonomy and domestication

The taurine and zebuine cattle breeds comprise the majority of the world cattle population but their taxonomic status is still controversial. The two forms of cattle are currently classified as Bos taurus and Bos indicus species and are differentiated primarily by the presence or absence of a hump. However, these two species hybridize readily, producing fully fertile offspring. We have determined and analyzed complete B. taurus and B. indicus mitochondrial genome sequences to investigate the extent of sequence divergences and to study their taxonomic status by molecular dating. The sequences encompassed 16,338 and 16,339 nucleotides, respectively, and differed at 237 positions. Estimated divergence times indicated that the two cattle lineages separated 1.7-2.0 million years ago. Combined phylogenetic analyses of 18 new and 130 previously reported extant B. taurus and B. indicus control region sequences with data from 32 archaeological specimens of the extinct wild aurochs (Bos primigenius) identified four major maternal lineages. B. primigenius haplotypes were present in all but the B. indicus lineage, and one B. taurus sequence clustered with B. primigenius P haplotypes that were not previously linked with domestic cattle. The B. indicus cluster and a recently reported new B. primigenius haplotype that represents a new lineage were approximately equidistant from the B. taurus cluster. These data suggest domestications from several differentiated populations of B. primigenius and a subspecies status for taurine (B. primigenius taurus) and zebuine (B. primigenius indicus) cattle.     

Cattle domestication in the Near East was followed by hybridization with aurochs bulls in Europe

Domesticated cattle were one of the cornerstones of European Neolithisation and are thought to have been introduced to Europe from areas of aurochs domestication in the Near East. This is consistent with mitochondrial DNA (mtDNA) data, where a clear separation exists between modern European cattle and ancient specimens of British aurochsen. However, we show that Y chromosome haplotypes of north European cattle breeds are more similar to haplotypes from ancient specimens of European aurochsen, than to contemporary cattle breeds from southern Europe and the Near East. There is a sharp north-south gradient across Europe among modern cattle breeds in the frequencies of two distinct Y chromosome haplotypes; the northern haplotype is found in 20 out of 21 European aurochsen or early domestic cattle dated 9500-1000 BC. This indicates that local hybridization with male aurochsen has left a paternal imprint on the genetic composition of modern central and north European breeds. Surreptitious mating between aurochs bulls and domestic cows may have been hard to avoid, or may have occurred intentionally to improve the breeding stock. Rather than originating from a few geographical areas only, as indicated by mtDNA, our data suggest that the origin of domestic cattle may be far more complex than previously thought.     

Early history of European domestic cattle as revealed by ancient DNA

We present an extensive ancient DNA analysis of mainly Neolithic cattle bones sampled from archaeological sites along the route of Neolithic expansion, from Turkey to North-Central Europe and Britain. We place this first reasonable population sample of Neolithic cattle mitochondrial DNA sequence diversity in context to illustrate the continuity of haplotype variation patterns from the first European domestic cattle to the present. Interestingly, the dominant Central European pattern, a starburst phylogeny around the modal sequence, T3, has a Neolithic origin, and the reduced diversity within this cluster in the ancient samples accords with their shorter history of post-domestic accumulation of mutation.     

Mitogenomes from Two Uncommon Haplogroups Mark Late Glacial/Postglacial Expansions from the Near East and Neolithic Dispersals within Europe

The current human mitochondrial (mtDNA) phylogeny does not equally represent all human populations but is biased in favour of representatives originally from north and central Europe. This especially affects the phylogeny of some uncommon West Eurasian haplogroups, including I and W, whose southern European and Near Eastern components are very poorly represented, suggesting that extensive hidden phylogenetic substructure remains to be uncovered. This study expanded and re-analysed the available datasets of I and W complete mtDNA genomes, reaching a comprehensive 419 mitogenomes, and searched for precise correlations between the ages and geographical distributions of their numerous newly identified subclades with events of human dispersal which contributed to the genetic formation of modern Europeans. Our results showed that haplogroups I (within N1a1b) and W originated in the Near East during the Last Glacial Maximum or pre-warming period (the period of gradual warming between the end of the LGM, ∼19 ky ago, and the beginning of the first main warming phase, ∼15 ky ago) and, like the much more common haplogroups J and T, may have been involved in Late Glacial expansions starting from the Near East. Thus our data contribute to a better definition of the Late and postglacial re-peopling of Europe, providing further evidence for the scenario that major population expansions started after the Last Glacial Maximum but before Neolithic times, but also evidencing traces of diffusion events in several I and W subclades dating to the European Neolithic and restricted to Europe.     

Mitochondrial haplotypes of European wild boars with 2n = 36 are closely related to those of European domestic pigs with 2n = 38

Wild boars from Western Europe have a 2n = 36 karyotype, in contrast to a karyotype of 2n = 38 in wild boars from Central Europe and Asia and in all domestic pigs. The phylogenetic status of this wild boar population is unclear, and it is not known if it has contributed to pig domestication. We have now sequenced the mtDNA control region from 30 European wild boars (22 with a confirmed 2n = 36 karyotype) and six Asian wild boars (two Hainan and four Dongbei wild boars) to address this question. The results revealed a close genetic relationship between mtDNA haplotypes from wild boars with 2n = 36 to those from domestic pigs with 2n = 38. Thus, we cannot exclude the possibility that wild boars with 2n = 36 may have contributed to pig domestication despite the karyotype difference. One of the European wild boars carried an Asian mtDNA haplotype, and this most likely reflects gene flow from domestic pigs to European wild boars. However, this gene flow does not appear to be extensive because the frequency of Asian haplotypes detected among European wild boars (c. 3%) were 10-fold lower than among European domestic pigs (c. 30%). Previous studies of mtDNA haplotypes have indicated that pig populations in Europe and Asia have experienced a population expansion, but it is not clear if the expansion occurred before or after domestication. The results of the present study are consistent with an expansion that primarily occurred prior to domestication because the mtDNA haplotypes found in European and Asian wild boars did not form their own clusters but were intermingled with haplotypes found in domestic pigs, indicating that they originated from the same population expansion.     

Mitochondrial DNA signals of late glacial recolonization of Europe from near eastern refugia

Human populations, along with those of many other species, are thought to have contracted into a number of refuge areas at the height of the last Ice Age. European populations are believed to be, to a large extent, the descendants of the inhabitants of these refugia, and some extant mtDNA lineages can be traced to refugia in Franco-Cantabria (haplogroups H1, H3, V, and U5b1), the Italian Peninsula (U5b3), and the East European Plain (U4 and U5a). Parts of the Near East, such as the Levant, were also continuously inhabited throughout the Last Glacial Maximum, but unlike western and eastern Europe, no archaeological or genetic evidence for Late Glacial expansions into Europe from the Near East has hitherto been discovered. Here we report, on the basis of an enlarged whole-genome mitochondrial database, that a substantial, perhaps predominant, signal from mitochondrial haplogroups J and T, previously thought to have spread primarily from the Near East into Europe with the Neolithic population, may in fact reflect dispersals during the Late Glacial period, ～19-12 thousand years (ka) ago.     

Mitochondrial DNA variation and evolution of Japanese black cattle (Bos taurus).

This article describes complete mitochondrial DNA displacement loop sequences from 32 Japanese Black cattle and the analysis of these data in conjunction with previously published sequences from African, European, and Indian subjects. The origins of North East Asian domesticated cattle are unclear. The earliest domestic cattle in the region were Bos taurus and may have been domesticated from local wild cattle (aurochsen; B. primigenius), or perhaps had an origin in migrants from the early domestic center of the Near East. In phylogenetic analyses, taurine sequences form a dense tree with a center consisting of intermingled European and Japanese sequences with one group of Japanese and another of all African sequences, each forming distinct clusters at extremes of the phylogeny. This topology and calibrated levels of sequence divergence suggest that the clusters may represent three different strains of ancestral aurochs, adopted at geographically and temporally separate stages of the domestication process. Unlike Africa, half of Japanese cattle sequences are topologically intermingled with the European variants. This suggests an interchange of variants that may be ancient, perhaps a legacy of the first introduction of domesticates to East Asia.     

Ancient DNA from South-East Europe Reveals Different Events during Early and Middle Neolithic Influencing the European Genetic Heritage

The importance of the process of Neolithization for the genetic make-up of European populations has been hotly debated, with shifting hypotheses from a demic diffusion (DD) to a cultural diffusion (CD) model. In this regard, ancient DNA data from the Balkan Peninsula, which is an important source of information to assess the process of Neolithization in Europe, is however missing. In the present study we show genetic information on ancient populations of the South-East of Europe. We assessed mtDNA from ten sites from the current territory of Romania, spanning a time-period from the Early Neolithic to the Late Bronze Age. mtDNA data from Early Neolithic farmers of the Starčevo Criş culture in Romania (Cârcea, Gura Baciului and Negrileşti sites), confirm their genetic relationship with those of the LBK culture (Linienbandkeramik Kultur) in Central Europe, and they show little genetic continuity with modern European populations. On the other hand, populations of the Middle-Late Neolithic (Boian, Zau and Gumelniţa cultures), supposedly a second wave of Neolithic migration from Anatolia, had a much stronger effect on the genetic heritage of the European populations. In contrast, we find a smaller contribution of Late Bronze Age migrations to the genetic composition of Europeans. Based on these findings, we propose that permeation of mtDNA lineages from a second wave of Middle-Late Neolithic migration from North-West Anatolia into the Balkan Peninsula and Central Europe represent an important contribution to the genetic shift between Early and Late Neolithic populations in Europe, and consequently to the genetic make-up of modern European populations.     

TYR基因外显子1的序列变异

为了分析家猪与野猪的遗传多样性及起源,测定了来自１２个中国地方家猪品种、３个欧洲引进猪品种以及8个中国野猪和2个越南野猪共36个个体的酪氨酸酶基因（TYR）外显子1的序列,共检出６个单核苷酸多态性位点（SNPs）,且这6个位点的变异均为同义突变,根据这些变异可将酪氨酸酶基因DNA序列归结为４种单倍型。结合已发表的数据,构建了简约中介网络图。 在网络图中,单倍型TYR*2主要为欧洲家猪与欧洲野猪和三条亚洲家猪染色体。大部分亚洲家猪和野猪共享单倍型TYR*1,表明这是一个亚洲类型的单倍型;同时也有部分欧洲家猪与野猪携带这一单倍型。 而单倍型TYR*3和TYR*4为本研究检测到的稀有单倍型,这两种单倍型主要由中国家猪与亚洲野猪组成。这种网络图结构支持家猪的欧洲和亚洲独立起源学说,同时也表明相当部分的欧洲家猪品种受到亚洲猪的基因渗透,而少量中国家猪和日本野猪也受到了欧洲猪的基因渗透。     

Adapt or die—Response of large herbivores to environmental changes in Europe during the Holocene

Climate warming and human landscape transformation during the Holocene resulted in environmental changes for wild animals. The last remnants of the European Pleistocene megafauna that survived into the Holocene were particularly vulnerable to changes in habitat. To track the response of habitat use and foraging of large herbivores to natural and anthropogenic changes in environmental conditions during the Holocene, we investigated carbon (δ¹³C) and nitrogen (δ¹⁵N) stable isotope composition in bone collagen of moose (Alces alces), European bison (Bison bonasus) and aurochs (Bos primigenius) in Central and Eastern Europe. We found strong variations in isotope compositions in the studied species throughout the Holocene and diverse responses to changing environmental conditions. All three species showed significant changes in their δ¹³C values reflecting a shift of foraging habitats from more open in the Early and pre‐Neolithic Holocene to more forest during the Neolithic and Late Holocene. This shift was strongest in European bison, suggesting higher plasticity, more limited in moose, and the least in aurochs. Significant increases of δ¹⁵N values in European bison and moose are evidence of a diet change towards more grazing, but may also reflect increased nitrogen in soils following deglaciation and global temperature increases. Among the factors explaining the observed isotope variations were time (age of samples), longitude and elevation in European bison, and time, longitude and forest cover in aurochs. None of the analysed factors explained isotope variations in moose. Our results demonstrate the strong influence of natural (forest expansion) and anthropogenic (deforestation and human pressure) changes on the foraging ecology of large herbivores, with forests playing a major role as a refugial habitat since the Neolithic, particularly for European bison and aurochs. We propose that high flexibility in foraging strategy was the key for survival of large herbivores in the changing environmental conditions of the Holocene.     

Novel Y‐chromosome short tandem repeats in Sus scrofa and their variation in European wild boar and domestic pig populations

Y‐chromosome markers are important tools for studying male‐specific gene flow within and between populations, hybridization patterns and kinship. However, their use in non‐human mammals is often hampered by the lack of Y‐specific polymorphic markers. We identified new male‐specific short tandem repeats (STRs) in Sus scrofa using the available genome sequence. We selected four polymorphic loci (5–10 alleles per locus), falling in one duplicated and two single‐copy regions. A total of 32 haplotypes were found by screening 211 individuals from eight wild boar populations across Europe and five domestic pig populations. European wild boar were characterized by significantly higher levels of haplotype diversity compared to European domestic pigs (H_D = 0.904 ± 0.011 and H_D = 0.491 ± 0.077 respectively). Relationships among STR haplotypes were investigated by combining them with single nucleotide polymorphisms at two linked genes (AMELY and UTY) in a network analysis. A differentiation between wild and domestic populations was observed (F_ST = 0.229), with commercial breeds sharing no Y haplotype with the sampled wild boar. Similarly, a certain degree of geographic differentiation was observed across Europe, with a number of local private haplotypes and high diversity in northern populations. The described Y‐chromosome markers can be useful to track male inheritance and gene flow in wild and domestic populations, promising to provide insights into evolutionary and population genetics in Sus scrofa.     

Modern taurine cattle descended from small number of near-eastern founders

Archaeozoological and genetic data indicate that taurine cattle were first domesticated from local wild ox (aurochs) in the Near East some 10,500 years ago. However, while modern mitochondrial DNA (mtDNA) variation indicates early Holocene founding event(s), a lack of ancient DNA data from the region of origin, variation in mutation rate estimates, and limited application of appropriate inference methodologies have resulted in uncertainty on the number of animals first domesticated. A large number would be expected if cattle domestication was a technologically straightforward and unexacting region-wide phenomenon, while a smaller number would be consistent with a more complex and challenging process. We report mtDNA sequences from 15 Neolithic to Iron Age Iranian domestic cattle and, in conjunction with modern data, use serial coalescent simulation and approximate Bayesian computation to estimate that around 80 female aurochs were initially domesticated. Such a low number is consistent with archaeological data indicating that initial domestication took place in a restricted area and suggests the process was constrained by the difficulty of sustained managing and breeding of the wild progenitors of domestic cattle.     

Between the Baltic and Danubian Worlds: The Genetic Affinities of a Middle Neolithic Population from Central Poland

For a long time, anthropological and genetic research on the Neolithic revolution in Europe was mainly concentrated on the mechanism of agricultural dispersal over different parts of the continent. Recently, attention has shifted towards population processes that occurred after the arrival of the first farmers, transforming the genetically very distinctive early Neolithic Linear Pottery Culture (LBK) and Mesolithic forager populations into present-day Central Europeans. The latest studies indicate that significant changes in this respect took place within the post-Linear Pottery cultures of the Early and Middle Neolithic which were a bridge between the allochthonous LBK and the first indigenous Neolithic culture of north-central Europe—the Funnel Beaker culture (TRB). The paper presents data on mtDNA haplotypes of a Middle Neolithic population dated to 4700/4600–4100/4000 BC belonging to the Brześć Kujawski Group of the Lengyel culture (BKG) from the Kuyavia region in north-central Poland. BKG communities constituted the border of the “Danubian World” in this part of Europe for approx. seven centuries, neighboring foragers of the North European Plain and the southern Baltic basin. MtDNA haplogroups were determined in 11 individuals, and four mtDNA macrohaplogroups were found (H, U5, T, and HV0). The overall haplogroup pattern did not deviate from other post-Linear Pottery populations from central Europe, although a complete lack of N1a and the presence of U5a are noteworthy. Of greatest importance is the observed link between the BKG and the TRB horizon, confirmed by an independent analysis of the craniometric variation of Mesolithic and Neolithic populations inhabiting central Europe. Estimated phylogenetic pattern suggests significant contribution of the post-Linear BKG communities to the origin of the subsequent Middle Neolithic cultures, such as the TRB.     

The Genetic Integrity of the Ex Situ Population of the European Wildcat (Felis silvestris silvestris) Is Seriously Threatened by Introgression from Domestic Cats (Felis silvestris catus)

Studies on the genetic diversity and relatedness of zoo populations are crucial for implementing successful breeding programmes. The European wildcat, Felis s. silvestris, is subject to intensive conservation measures, including captive breeding and reintroduction. We here present the first systematic genetic analysis of the captive population of Felis s. silvestris in comparison with a natural wild population. We used microsatellites and mtDNA sequencing to assess genetic diversity, structure and integrity of the ex situ population. Our results show that the ex situ population of the European wildcat is highly structured and that it has a higher genetic diversity than the studied wild population. Some genetic clusters matched the breeding lines of certain zoos or groups of zoos that often exchanged individuals. Two mitochondrial haplotype groups were detected in the in situ populations, one of which was closely related to the most common haplotype found in domestic cats, suggesting past introgression in the wild. Although native haplotypes were also found in the captive population, the majority (68%) of captive individuals shared a common mtDNA haplotype with the domestic cat (Felis s. catus). Only six captive individuals (7.7%) were assigned as wildcats in the STRUCTURE analysis (at K = 2), two of which had domestic cat mtDNA haplotypes and only two captive individuals were assigned as purebred wildcats by NewHybrids. These results suggest that the high genetic diversity of the captive population has been caused by admixture with domestic cats. Therefore, the captive population cannot be recommended for further breeding and reintroduction.     

Untangling the hybrid nature of modern pig genomes: a mosaic derived from biogeographically distinct and highly divergent Sus scrofa populations

The merging of populations after an extended period of isolation and divergence is a common phenomenon, in natural settings as well as due to human interference. Individuals with such hybrid origins contain genomes that essentially form a mosaic of different histories and demographies. Pigs are an excellent model species to study hybridization because European and Asian wild boars diverged ~1.2 Mya, and pigs were domesticated independently in Europe and Asia. During the Industrial Revolution in England, pigs were imported from China to improve the local pigs. This study utilizes the latest genomics tools to identify the origin of haplotypes in European domesticated pigs that are descendant from Asian and European populations. Our results reveal fine‐scale haplotype structure representing different ancient demographic events, as well as a mosaic composition of those distinct histories due to recently introgressed haplotypes in the pig genome. As a consequence, nucleotide diversity in the genome of European domesticated pigs is higher when at least one haplotype of Asian origin is present, and haplotype length correlates negatively with recombination frequency and nucleotide diversity. Another consequence is that the inference of past effective population size is influenced by the background of the haplotypes in an individual, but we demonstrate that by careful sorting based on the origin of haplotypes, both distinct demographic histories can be reconstructed. Future detailed mapping of the genomic distribution of variation will enable a targeted approach to increase genetic diversity of captive and wild populations, thus facilitating conservation efforts in the near future.     

Ancestral matrilineages and mitochondrial DNA diversity of the Lidia cattle breed

To clarify the genetic ancestry and the mitochondrial DNA (mtDNA) diversity of the Lidia cattle breed, a 521-bp D-loop fragment was sequenced in 527 animals belonging to 70 herds distributed across 29 lineages. The mtDNA diversity recorded was similar to that seen for Middle Eastern breeds and greater than that recorded for the majority of European breeds. Haplotype T3 was the most common (81%), followed by the African T1 haplotype (17%); very low frequencies were recorded for haplotypes T and T2. The results agree with there being two major ancestral lines for the Lidia breed, European and African, similar to that seen for other Mediterranean breeds. A wide range of variation in haplotype frequencies was seen between the examined lineages. Haplotype T3 was present in all those analysed; in five it was the only one present, and in only one lineage (Miura) was its frequency lower than that of T1. T1*, a haplotype reported in Criollo breeds and to date in only a single European breed (the Retinta breed from Spain), was found in a single animal belonging to the Concha y Sierra lineage. Network analysis of the Lidia breed revealed the presence of two major haplotypes: T3 and T1. The Lidia breed appears to be more closely related to prehistoric Iberian and Italian than to British aurochs.     

A genetic method to distinguish crossbred Inobuta from Japanese wild boars

To distinguish pig-wild boar crossbred Inobuta from Japanese wild boar populations, a genetic method by using mitochondrial DNA (mtDNA) haplotypes and the nuclear glucosephosphate isomerase-processed pseudogene (GPIP) was developed. Sixteen mtDNA haplotypes from 152 wild boars from Kyushu, Shikoku and Honshu islands of Japan were distinct from those from Asian and European domestic pigs. Five alleles of GPIP were classified into two groups: 1). alleles GPIP*1, GPIP*3 and GPIP*3a from Japanese wild boars, Asian wild boars and domestic pigs; 2). alleles GPIP*4 and GPIP*4a from European wild boars and domestic pigs. An extensive genetic survey was done to distinguish the crossbred Inobuta from 60 wild boars hunted on Tsushima Island, Goto Islands, and Nagasaki and Ooita Prefectures. The mtDNA haplotypes from the 60 samples showed Japanese wild boars, but four wild boar samples from Nagasaki Prefecture had the European GPIP allele, GPIP*4. These results showed that nuclear DNA polymorphism analysis is useful, in addition to mtDNA haplotype assay, to detect "Inobuta" having the European genotype from Japanese wild boar populations.     

The Origins of Domesticated Cattle

The origin of domestic cattle has perplexed archaeologists for more than a century. Researchers have proposed various theories, which offer alternative spatial and chronological models for the origin and spread of domesticated cattle. One point of discussion is whether domestic cattle had a single or multiple origins; however, most authorities considered that the first steps towards cattle domestication were taken in southwest Asia and that domesticated cattle entered Europe with pastoralists migrating from this region. Domesticated taurine cattle were thought to have entered Africa in successive waves from southwest Asia, while zebu cattle migrated into Africa at a later date from Arabia and the Indian subcontinent. Analysis of mitochondrial DNA (mtDNA) shows that taurine and zebu cattle divergence before the Holocene and were probably domesticated independently. Recent mtDNA sequence data shows that African and European taurine cattle were probably domesticated independently, but that there was a process of genetic introgression between taurine and zebu cattle in Africa. Ancient DNA studies over the last 10 years suggest that Northern European aurochsens apparently contributed little or nothing to domestic cattle while Southern European aurochsens apparently made a significant input. However, Middle Eastern aurochsen, unfortunately not typed yet, are expected to be to be very similar to European breeds as well, both because archeological data suggest that the major center of domestication for European Bos taurus breeds was the Fertile Crescent (9), and also because a mtDNA sequence from a Syrian specimen dated at 8,000–9,000 years ago shows a typical European haplotype found both in modern breeds and the Italian aurochsen. Evidence seems to suggest that small to moderate levels of local gene flow from wild Bos primigenius females in selected breeds were either accepted or may be reinforced by Neolithic breeders.