Ancient cattle genomics, origins, and rapid turnover in the Fertile Crescent

See allHide authors and affiliations

Science  12 Jul 2019:
Vol. 365, Issue 6449, pp. 173-176
DOI: 10.1126/science.aav1002

How cow genomes have moo-ved

Cattle were domesticated ∼10,000 years ago, but analysis of modern breeds has not elucidated their origins. Verdugo et al. performed genome-wide analysis of 67 ancient Near Eastern Bos taurus DNA samples. Several populations of ancient aurochs were progenitors of domestic cows. These genetic lineages mixed ∼4000 years ago in a region around the Indus Valley. Interestingly, mitochondrial analysis indicated that genetic material likely derived from arid-adapted Bos indicus (zebu) bulls was introduced by introgression.

Science, this issue p. 173


Genome-wide analysis of 67 ancient Near Eastern cattle, Bos taurus, remains reveals regional variation that has since been obscured by admixture in modern populations. Comparisons of genomes of early domestic cattle to their aurochs progenitors identify diverse origins with separate introgressions of wild stock. A later region-wide Bronze Age shift indicates rapid and widespread introgression of zebu, Bos indicus, from the Indus Valley. This process was likely stimulated at the onset of the current geological age, ~4.2 thousand years ago, by a widespread multicentury drought. In contrast to genome-wide admixture, mitochondrial DNA stasis supports that this introgression was male-driven, suggesting that selection of arid-adapted zebu bulls enhanced herd survival. This human-mediated migration of zebu-derived genetics has continued through millennia, altering tropical herding on each continent.

The extinct Eurasian aurochs (Bos primigenius) was domesticated circa 10,500 years before present (yr B.P.) within the restricted locality of the Upper Euphrates and Tigris drainages of the Fertile Crescent (1, 2). However, the true extent and nature of interactions between humans and aurochs resulting in modern day domestic cattle are obscure.

Mitochondrial DNA (mtDNA) diversity in modern Bos taurus cattle suggests a highly restricted initial domestic pool of ~80 females (35). However, a more complex relationship with wild populations is evidenced by introgression from local aurochs into British cattle and the genomic divergence of B. indicus (zebu) cattle from the Indus Valley region (6, 7). Zebu genomic influence is pervasive in modern Near Eastern herds (8). Two theories account for this: one suggests an origin from genomically intermediate Near Eastern aurochs, whereas a second hypothesizes that these Near Eastern herds resulted from an introgression of domestic zebu genomes into the region from the east, either in a discrete active process—perhaps responding to climate fluctuation—or a passive diffusion over many millennia (9).

To analyze now-obscured early cattle genome strata from the region of B. taurus domestication, we retrieved genome-wide data from 67 ancient bovines (including six aurochs). These date from Mesolithic to early Islamic periods, and despite poor preservation, which is typical of the region, we obtained an average genome coverage of 0.9× (table S1).

The pattern of genetic variation in extant cattle is well established. European B. taurus, West African B. taurus, and B. indicus of South Asian origin represent three distinct apices in plotted principal components (PCs) (Fig. 1A). Geographically intermediate populations, such as Near Eastern and East African animals, fall in genetically intermediate positions (7, 8, 10). Projecting ancient cattle genomes (provenance shown in Fig. 1B) against this genetic landscape (Fig. 1A), we observe that to the left of PC1, earlier (Neolithic and Bronze Age) genomes fall in three geographical clusters (a, Balkans; b, Anatolia/Iran; and c, southern Levant) along with modern European and African B. taurus, whereas B. indicus breeds are separated and represented on the far right of the PC plot (Fig. 1A). This suggests that cattle origins included two divergent aurochs populations that formed the basis of the B. indicusB. taurus divide.

Fig. 1 Procrustes projection principal components analysis of ancient cattle.

(A) Ancient animals are projected on modern 770K Bovine single-nucleotide polymorphism (SNP) genotypes, shown as background gray asterisks (figs. S1 and S2). Four clusters are highlighted: Neolithic Balkans, which plot with modern Europeans (a); a group of mainly Anatolian and Iranian cattle close to four aurochs from the Near East (b); and Levantine cattle that fall into two groups, a cluster of earlier samples (c) and a cluster of later samples (d) close to contemporary Near Eastern cattle with B. indicus admixture. (B) Approximate geographical distribution of ancient sample sites.

Six ancient aurochs genomes, including four from the greater Near East, provide additional context: two ~9000-year-old samples from the Levantine Aceramic village of Abu Ghosh (Abu1 and Abu2), a 7500-year-old sample from the early Anatolian settlement Çatalhöyük (Ch22), and a 7000-year-old Armenian aurochs (Gyu2) (11). These four genomes fall close to the Anatolia and Iran ancient domestic cattle cluster (Fig. 1A, cluster b) and reveal this as the oldest ancestral stratum of B. taurus. The genomic signature of this earliest population has been obscured in modern Near Eastern cattle by later admixture. From this group, we sequenced a well-preserved 8000-year-old Anatolian genome (Sub1) (11) to 13.5× coverage and use this in D statistics testing for zebu introgression in other ancient individuals (Fig. 2).

Fig. 2 Zebu introgression through time in the greater Near East.

Whole genome D statistic calculations with a gaur (B. gaurus) as an outgroup and the Neolithic Anatolian domestic genome Sub1 as a representative nonadmixed individual (see inset). Lower coverage samples with <200 informative sites were excluded. A step change in zebu introgression is apparent circa 4000 yr B.P. MtDNA counts of taurine (gray) and the single zebu (black) mtDNA in ancient domestic cattle are graphed (bottom, with a shared time axis) in 250-year intervals.

B. indicus cattle are adapted to, and predominate in, modern arid and tropical regions of the world (11). Zebu cattle originated circa 8000 yr B.P. (12). However, despite archaeological evidence for contact between civilizations of the Fertile Crescent region and the Indus Valley (9), the influence of the zebu genome is detectable in ancient Southwest Asian cattle only 4000 years later (Fig. 2). However, after ~4000 yr B.P., hybrid animals (median 35% indicine ancestry) are found across the Near East, from Central Asia and Iran to the Caucasus and Mediterranean shores of the southern Levant (table S2 and fig. S1). During this period, depictions and osteological evidence for B. indicus also appear in the region (9, 13). In contrast to autosomal data, but similar to earlier work (14), we find persistence of B. taurus mitochondria, suggesting introgression may have been mediated by bulls (Fig. 2).

This sharp influx may have been stimulated by the onset of a period of increased aridity known as the 4.2-thousand-year abrupt climate change event (9, 1517). This multicentury drought coincided with empire collapse in both Mesopotamia and Egypt as well as a decline in the Indus civilization and has been accepted as the boundary defining the onset of our current geological age, the Meghalayan (18).

Three features of this zebu influx after ~4000 yr B.P. attest that the influx was likely driven by adaptation and/or human agency rather than passive diffusion. First, the extent of indicine introgression does not follow a simple east-to-west gradient; for example, it is pronounced in Levantine genomes from the western edge of the Near East. Second, the introgression was widespread and took place in a relatively restricted time interval after four millennia of barely detectable B. indicus influence. Third, it was plausibly driven by bull choice, as we observe up to ~70% autosomal genome change but a retained substratum of B. taurus mtDNA haplotypes (Fig. 2 and table S3). Hybrid B. taurus–B. indicus herds may have enabled the survival of communities under stress and perhaps facilitated expansion of herding into more-peripheral regions. Restocking after herd decline may have also been a factor. Westward human migration has been documented around this time (19, 20) along with archaeological evidence for the appearance of other South Asian taxa such as water buffalo and Asian elephants in the Near East (21), suggesting the movement of large animals by people.

Before zebu admixture, ancient southern Levantine animals occupy a distinctive space within the PC plot (Fig. 1A, cluster c), toward modern African cattle and adjacent to a 9000-yr-B.P. Epipalaeolithic Moroccan aurochs (Th7). A 7000-yr-B.P. Mesolithic British aurochs genome (CPC98) (6) also plots away from the core Anatolia/Iran ancestral Near Eastern cluster and close to Neolithic Balkan (cluster a) and modern European cattle. These genetic affinities in ancient cattle suggest an early secondary recruitment from diverse wild populations.

Concordantly, D statistic tests of allele sharing by cattle population pairs with three divergent aurochs confirm that early cattle exhibit asymmetric relationships with different wild populations (Fig. 3). The most extreme deviations are found in comparisons featuring the B. taurus Levantine population (Fig. 1, cluster c); these share the least affinity with the British and Armenian aurochs (z-score > 5.67; P < 10−5) but more with the Moroccan Epipalaeolithic sample. We infer that a distinct strain of aurochs, probably from the Levant and similar to those ranging across North Africa, had considerable input into early cattle in the southern Levant. The Mesolithic British aurochs also shows asymmetric affinity with the Neolithic Balkans samples, implying that the hybridization of European aurochs (6) was initiated more than 7000 years ago, close to the onset of human herding of cattle in Europe. These findings are supported by a qpgraph analysis (figs. S2 and S3) and cannot be explained by cattle-into-aurochs admixture, as both the British and Moroccan aurochs have securely predomestic dates. Although each of these three aurochs have divergent mtDNA haplotypes falling outside normal B. taurus variation, ancient domesticates display typical modern domestic haplotypes (fig. S4). This points toward common matrilineal origins for domestic taurine cattle and away from an archaeologically less parsimonious interpretation that our observed ancient genetic structure may have arisen from separate domestications; it also suggests that introgression may have been via mating with wild males. Sexually mature bulls, because of size and aggression, were likely the most dangerous stock in Neolithic villages, and thus unsupervised field insemination by aurochs bulls may have played a role in early herd management (22).

Fig. 3 Clade integrity of ancient population pairs with respect to aurochs introgression.

(A) The test D(gaur, aurochs; ancient group1, ancient group2) reveals asymmetric affinities of aurochs genomes with pre-4000-yr-B.P. cattle populations. Levantine cattle show reduced allele sharing relative to other populations with the Armenian (Gyu2) and British (CPC98) aurochs but more with the Moroccan aurochs (Th7). Balkan cattle show asymmetric affinities with the British aurochs. Bars denote two standard errors. (B) Geographical location of aurochs tested. (C) Distribution of ancient domestic cattle groups tested; post-4000-yr-B.P. Near Eastern samples were excluded because of their zebu admixture.

Distinct genotypes and phenotypes in B. taurus cattle native to Africa, such as tolerance of tropical infections, have been attributed to either local domestication or introgression from African aurochs (10, 23). However, ancient Levantine genome affinity with North African aurochs hints that this distinctiveness may have origins in the southern Fertile Crescent. Supporting this, the B. taurus mtDNA haplogroup (T1), which is almost fixed in African cattle populations (24), is the most frequent in the southern Levant, including earliest samples, but was not found among other ancient domesticates (table S3).

B. taurus were initially derived from a restricted northern Fertile Crescent genetic background, but early domestic cattle outside this region gained heterogeneous inputs from diverse aurochs strains, including contributions specific to European and African cattle ancestors. After ~4200 yr B.P., gross genome turnover reflecting the spread of B. indicus, and likely associated with climate change, was effected by cattle herders throughout southwest and central Asia, representing the start of a global B. indicus genome diaspora (25) that continues today.

Supplementary Materials

Materials and Methods

Figs. S1 to S14

Tables S1 to S13

References (27200)

  • Independent researcher

References and Notes

  1. See the supplementary materials.
Acknowledgments: We thank L. Frantz, G. Larson, and M. Djamali for valuable suggestions. We thank L. Cassidy for helpful discussions and scripts, and R. Verdugo for contributions in figure design. We thank excavators, archaeozoologists, and museums who permitted sampling of bones from their excavations and collections without which this project would not have been possible, including E. Galili, T. Levy, C. Grigson, R. Gophna, A. Maeir, S. Gitin, A. Ben-Tor, D. Master, E. van den Brink, S. Cohen, E. Rosenstock, P. Biehl, I. Hodder, J.J. Roodenberg, S. Alpaslan Roodenberg, D. Helmer, H. Greenfield, J. Vuković, M. Radivojević, B. Roberts, M. Marić, G. Grabež, J. Kuzmanović-Cvetković, A. Bulatović, M. Spasić, E. Rosenstock, K. Kaniuth, and B. Eivers. We thank the Israel Antiquities Authority for permitting sampling of the Israeli sites (under permit). We thank the Çatalhöyük Research Project and Ministry of Culture and Tourism of the Republic of Turkey. We thank the Iranian Cultural Heritage Handicraft and Tourism organization and the National Museum of Iran (NMI) (J. Nokandeh, director, and F. Biglari, cultural deputy). We are grateful to H. Laleh and A. Aliyari, directors of the Bioarchaeology Laboratory, Central Laboratory, University of Tehran. The ATM Project of MNHN supported sampling of several sites as well as the LIA HAOMA CNRS project. The participation of ZIN RAS (state assignment AAAA-A17-117022810195-3) to this research is acknowledged. Funding: Supported by ERC Investigator grant 295729-CodeX. Additional support from Science Foundation Ireland Award 12/ERC/B2227, Trinseq and the SFI/HEA Irish Centre for High-End Computing (ICHEC). M.D.T. was partially supported by the Marie Skłodowska-Curie Individual Fellowship H2020-MSCA-IF-2016 747424 (SCRIBE). V.E.M. was partially supported by the NERC Grant NE/P012574/1. M.J.C. acknowledges support from DNRF128. P.M.D. was supported by the HERA Joint Research Programme “Uses of the Past” (CitiGen) and the European Union’s Horizon 2020 research and innovation program (649307). Author contributions: D.G.B. conceived of the project and designed research, with input from J.B. and M.J.C.; M.P.V., V.E.M., A.S., V.M., K.G.D., M.D.T., and A.J.H. performed ancient DNA laboratory work; M.D.T. and V.M. performed modern DNA laboratory work; M.P.V. and V.E.M. performed the computational analyses with input from D.G.B., M.D.T., P.M.D., and K.G.D.; other coauthors provided data and samples. D.G.B., M.P.V., V.E.M., A.S., and M.D.T. wrote the paper and supplementary information with input from all coauthors. Competing interests: The authors declare that they have no competing interests. Data and materials availability: ​Raw reads have been deposited at the European Nucleotide Archive (ENA) with project number PRJEB31621. The code used in this study for manipulating sequence files can be found at Zenodo (26). Mitochondrial phylogenies and modern SNP data are available at

Stay Connected to Science

Navigate This Article