Prehistoric contacts over the Straits of Gibraltar indicated by genetic analysis of Iberian Bronze Age cattle

The Portuguese mertolenga cattle breed: a review

This review summarizes the state of the art regarding the phylogenetic origins, recent history and present-day main traits and uses of the Mertolenga cattle breed from Southern Portugal, particularly those related to production performances and product quality. Named after the historical city of Mértola, in southern Portugal, the Mertolenga is one of the fifteen autochthonous bovine breeds of Portugal. It is a cattle breed thoroughly adapted to the poor Mediterranean pastures of the southern regions of the Iberian Peninsula. It is used predominantly in an extensive to semi-intensive sylvopastoral production system called montado, where pastures are combined with helm and cork trees. Its productive traits allow for a good adaptation to the intense dry heat and pasture shortage during the summer, and compensatory growth in autumn and spring, when pastures regenerate. They are small to medium sized animals, with well-balanced bodies, roan, red, or, less often, red-spotted coats, and known for their nervous temperament. Although this breed experienced a severe decline in numbers in the 1970s and classified as endangered in the 1990s, the work of a few breeders led to the establishment of larger Mertolenga breed inventories, starting from a limited base. For this reason, the entire breed has today a strong influence from a very few herds and sires. Reproduction is still mostly achieved using natural mating, and the males are often kept with the breeding females all year long. It is a heterogeneous breed both phenotypically and genetically. Recent phylogenetic studies have revealed the Mertolenga as a one of the most genetically diverse breeds in the country and in the Iberian Peninsula and helped classify this breed, once believed to be a variety of the Alentejana breed. These studies also showed genetic relations with other breeds in the Iberian Peninsula. Mertolenga beef currently benefits from several certifications, the most important one being the PDO - Protected Denomination of Origin.

Late History of Cattle Breeds in Central Europe in Light of Genetic and Archaeogenetic Sources-Overview, Thoughts, and Perspectives

Although Europe was not a primary centre of cattle domestication, its expansion from the Middle East and subsequent development created a complex pattern of cattle breed diversity. Many isolated populations of local historical breeds still carry the message about the physical and genetic traits of ancient populations. Since the way of life of human communities starting from the eleventh millennium BP was strongly determined by livestock husbandry, the knowledge of cattle diversity through the ages is helpful in the interpretation of many archaeological findings. Historical cattle diversity is currently at the intersection of two leading directions of genetic research. Firstly, it is archaeogenetics attempting to recover and interpret the preserved genetic information directly from archaeological finds. The advanced archaeogenetic approaches meet with the population genomics of extant cattle populations. The immense amount of genetic information collected from living cattle, due to its key economic role, allows for reconstructing the genetic profiles of the ancient populations backwards. The present paper aims to place selected archaeogenetic, genetic, and genomic findings in the picture of cattle history in Central Europe, as suggested by archaeozoological and historical records. Perspectives of the methodical connection between the genetic approaches and the approaches of traditional archaeozoology, such as osteomorphology and osteometry, are discussed. The importance, actuality, and effectiveness of combining different approaches to each archaeological find, such as morphological characterization, interpretation of the historical context, and molecular data, are stressed.

Ancient DNA confirms diverse origins of early post-Columbian cattle in the Americas

Before the arrival of Europeans, domestic cattle (Bos taurus) did not exist in the Americas, and most of our knowledge about how domestic bovines first arrived in the Western Hemisphere is based on historical documents. Sixteenth-century colonial accounts suggest that the first cattle were brought in small numbers from the southern Iberian Peninsula via the Canary archipelago to the Caribbean islands where they were bred locally and imported to other circum-Caribbean regions. Modern American heritage cattle genetics and limited ancient mtDNA data from archaeological colonial cattle suggest a more complex story of mixed ancestries from Europe and Africa. So far little information exists to understand the nature and timing of the arrival of these mixed-ancestry populations. In this study we combine ancient mitochondrial and nuclear DNA from a robust sample of some of the earliest archaeological specimens from Caribbean and Mesoamerican sites to clarify the origins and the dynamics of bovine introduction into the Americas. Our analyses support first arrival of cattle from diverse locales and potentially confirm the early arrival of African-sourced cattle in the Americas, followed by waves of later introductions from various sources over several centuries.

Iron age genomic data from Althiburos - Tunisia renew the debate on the origins of African taurine cattle

The Maghreb is a key region for understanding the dynamics of cattle dispersal and admixture with local aurochs following their earliest domestication in the Fertile Crescent more than 10,000 years ago. Here, we present data on autosomal genomes and mitogenomes obtained for four archaeological specimens of Iron Age (∼2,800 cal BP-2,000 cal BP) domestic cattle from the Eastern Maghreb, i.e. Althiburos (El Kef, Tunisia). D-loop sequences were obtained for an additional eight cattle specimens from this site. Maternal lineages were assigned to the elusive R and ubiquitous African-T1 haplogroups found in two and ten Althiburos specimens, respectively. Our results can be explained by post-domestication hybridization of Althiburos cattle with local aurochs. However, we cannot rule out an independent domestication in North Africa considering the shared ancestry of Althiburos cattle with the pre-domestic Moroccan aurochs and present-day African taurine cattle.

Genetic Structure Analysis of 155 Transboundary and Local Populations of Cattle (Bos taurus, Bos indicus and Bos grunniens) Based on STR Markers

Every week, 1-2 breeds of farm animals, including local cattle, disappear in the world. As the keepers of rare allelic variants, native breeds potentially expand the range of genetic solutions to possible problems of the future, which means that the study of the genetic structure of these breeds is an urgent task. Providing nomadic herders with valuable resources necessary for life, domestic yaks have also become an important object of study. In order to determine the population genetic characteristics, and clarify the phylogenetic relationships of modern representatives of 155 cattle populations from different regions of the world, we collected a large set of STR data (10,250 individuals), including unique native cattle, 12 yak populations from Russia, Mongolia and Kyrgyzstan, as well as zebu breeds. Estimation of main population genetic parameters, phylogenetic analysis, principal component analysis and Bayesian cluster analysis allowed us to refine genetic structure and provided insights in relationships of native populations, transboundary breeds and populations of domestic yak. Our results can find practical application in conservation programs of endangered breeds, as well as become the basis for future fundamental research.

Recovery of mitogenomes from whole genome sequences to infer maternal diversity in 1883 modern taurine and indicine cattle

Maternal diversity based on a sub-region of mitochondrial genome or variants were commonly used to understand past demographic events in livestock. Additionally, there is growing evidence of direct association of mitochondrial genetic variants with a range of phenotypes. Therefore, this study used complete bovine mitogenomes from a large sequence database to explore the full spectrum of maternal diversity. Mitogenome diversity was evaluated among 1883 animals representing 156 globally important cattle breeds. Overall, the mitogenomes were diverse: presenting 11 major haplogroups, expanding to 1309 unique haplotypes, with nucleotide diversity 0.011 and haplotype diversity 0.999. A small proportion of African taurine (3.5%) and indicine (1.3%) haplogroups were found among the European taurine breeds and composites. The haplogrouping was largely consistent with the population structure derived from alternate clustering methods (e.g. PCA and hierarchical clustering). Further, we present evidence confirming a new indicine subgroup (I1a, 64 animals) mainly consisting of breeds originating from China and characterised by two private mutations within the I1 haplogroup. The total genetic variation was attributed mainly to within-breed variance (96.9%). The accuracy of the imputation of missing genotypes was high (99.8%) except for the relatively rare heteroplasmic genotypes, suggesting the potential for trait association studies within a breed.

How have sheep breeds differentiated from each other in Morocco? Genetic structure and geographical distribution patterns

Background

Based on the relatively homogeneous origin of the sheep breeds in Morocco that originate mainly from Iberia, it is highly relevant to address the question of how these very diverse sheep populations differentiated from each other. The Mountains of the High Atlas and Middle Atlas are expected to constitute North–South and West–East geographical barriers, respectively, which could have shaped the history of the differentiation of sheep breeds. The aim of this study was to test this hypothesis by considering the genetic structure and the spatial distribution of five major breeds (Sardi, Timahdite, Beni Guil, Boujaad and D’man) and one minor breed (Blanche de Montagne), by analysing the mtDNA control region, using 30 individuals per breed.

Results

Phylogenetic and network analyses did not indicate any clear separation among the studied breeds and discriminant component principal analysis showed some overlap between them, which indicates a common genetic background. The calculated pairwise F_ST values and Nei’s genetic distances revealed that most breeds showed a moderate genetic differentiation. The lowest and highest degrees of differentiation were retrieved in the Beni Guil and Boujaad breeds, respectively. Analysis of molecular variance (AMOVA) indicated that more than 95% of the genetic diversity occurs within individuals, while between- and within-population variabilities represent only 1.332% and 2.881%, respectively. Isolation-by-distance, spatial Principal Component Analysis (sPCA), and spatial AMOVA analyses evidenced clear examples of geographical structuration among the breeds, both between and within breeds. However, several enigmatic relationships remain, which suggest the occurrence of complex events leading to breed differentiation.

Conclusions

The approaches used here resulted in a convergent view on the hypothetic events that could have led to the progressive differentiation between the Moroccan breeds. The major split seems to be linked to the West–East barrier of the Middle Atlas, whereas the influence of the High Atlas is less obvious and incompletely resolved. The study of additional breeds that have settled near the High Atlas should clarify the relationships between the breeds of the West part of the country, in spite of their small population size.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00679-2.

Biomolecular insights into North African-related ancestry, mobility and diet in eleventh-century Al-Andalus

Historical records document medieval immigration from North Africa to Iberia to create Islamic al-Andalus. Here, we present a low-coverage genome of an eleventh century CE man buried in an Islamic necropolis in Segorbe, near Valencia, Spain. Uniparental lineages indicate North African ancestry, but at the autosomal level he displays a mosaic of North African and European-like ancestries, distinct from any present-day population. Altogether, the genome-wide evidence, stable isotope results and the age of the burial indicate that his ancestry was ultimately a result of admixture between recently arrived Amazigh people (Berbers) and the population inhabiting the Peninsula prior to the Islamic conquest. We detect differences between our sample and a previously published group of contemporary individuals from Valencia, exemplifying how detailed, small-scale aDNA studies can illuminate fine-grained regional and temporal differences. His genome demonstrates how ancient DNA studies can capture portraits of past genetic variation that have been erased by later demographic shifts-in this case, most likely the seventeenth century CE expulsion of formerly Islamic communities as tolerance dissipated following the Reconquista by the Catholic kingdoms of the north.

The genetic ancestry of American Creole cattle inferred from uniparental and autosomal genetic markers

Cattle imported from the Iberian Peninsula spread throughout America in the early years of discovery and colonization to originate Creole breeds, which adapted to a wide diversity of environments and later received influences from other origins, including zebu cattle in more recent years. We analyzed uniparental genetic markers and autosomal microsatellites in DNA samples from 114 cattle breeds distributed worldwide, including 40 Creole breeds representing the whole American continent, and samples from the Iberian Peninsula, British islands, Continental Europe, Africa and American zebu. We show that Creole breeds differ considerably from each other, and most have their own identity or group with others from neighboring regions. Results with mtDNA indicate that T1c-lineages are rare in Iberia but common in Africa and are well represented in Creoles from Brazil and Colombia, lending support to a direct African influence on Creoles. This is reinforced by the sharing of a unique Y-haplotype between cattle from Mozambique and Creoles from Argentina. Autosomal microsatellites indicate that Creoles occupy an intermediate position between African and European breeds, and some Creoles show a clear Iberian signature. Our results confirm the mixed ancestry of American Creole cattle and the role that African cattle have played in their development.

Genetic Diversity and Population Structure of Ethiopian Sheep Populations Revealed by High-Density SNP Markers

Sheep in Ethiopia are adapted to a wide range of environments, including extreme habitats. Elucidating their genetic diversity is critical for improving breeding strategies and mapping quantitative trait loci associated with productivity. To this end, the present study investigated the genetic diversity and population structure of five Ethiopian sheep populations exhibiting distinct phenotypes and sampled from distinct production environments, including arid lowlands and highlands. To investigate the genetic relationships in greater detail and infer population structure of Ethiopian sheep breeds at the continental and global levels, we analyzed genotypic data of selected sheep breeds from the Ovine SNP50K HapMap dataset. All Ethiopian sheep samples were genotyped with Ovine Infinium HD SNP BeadChip (600K). Mean genetic diversity ranged from 0.29 in Arsi-Bale to 0.32 in Menz sheep, while estimates of genetic differentiation among populations ranged from 0.02 to 0.07, indicating low to moderate differentiation. An analysis of molecular variance revealed that 94.62 and 5.38% of the genetic variation was attributable to differences within and among populations, respectively. Our population structure analysis revealed clustering of five Ethiopian sheep populations according to tail phenotype and geographic origin-i.e., short fat-tailed (very cool high-altitude), long fat-tailed (mid to high-altitude), and fat-rumped (arid low-altitude), with clear evidence of admixture between long fat-tailed populations. North African sheep breeds showed higher levels of within-breed diversity, but were less differentiated than breeds from Eastern and Southern Africa. When African breeds were grouped according to geographic origin (North, South, and East), statistically significant differences were detected among groups (regions). A comparison of population structure between Ethiopian and global sheep breeds showed that fat-tailed breeds from Eastern and Southern Africa clustered together, suggesting that these breeds were introduced to the African continent via the Horn and migrated further south.

Contacts in the last 90,000 years over the Strait of Gibraltar evidenced by genetic analysis of wild boar (Sus scrofa)

Contacts across the Strait of Gibraltar in the Pleistocene have been studied in different research papers, which have demonstrated that this apparent barrier has been permeable to human and fauna movements in both directions. Our study, based on the genetic analysis of wild boar (Sus scrofa), suggests that there has been contact between Africa and Europe through the Strait of Gibraltar in the Late Pleistocene (at least in the last 90,000 years), as shown by the partial analysis of mitochondrial DNA. Cytochrome b and the control region from North African wild boar indicate a close relationship with European wild boar, and even some specimens belong to a common haplotype in Europe. The analyses suggest the transformation of the wild boar phylogeography in North Africa by the emergence of a natural communication route in times when sea levels fell due to climatic changes, and possibly through human action, since contacts coincide with both the Last Glacial period and the increasing human dispersion via the strait.

Abundance and Characterization of Perfect Microsatellites on the Cattle Y Chromosome

Microsatellites or simple sequence repeats (SSRs) are found in most organisms and play an important role in genomic organization and function. To characterize the abundance of SSRs (1-6 base-pairs [bp]) on the cattle Y chromsome, the relative frequency and density of perfect or uninterrupted SSRs based on the published Y chromosome sequence were examined. A total of 17,273 perfect SSRs were found, with total length of 324.78 kb, indicating that approximately 0.75% of the cattle Y chromosome sequence (43.30 Mb) comprises perfect SSRs, with an average length of 18.80 bp. The relative frequency and density were 398.92 loci/Mb and 7500.62 bp/Mb, respectively. The proportions of the six classes of perfect SSRs were highly variable on the cattle Y chromosome. Mononucleotide repeats had a total number of 8073 (46.74%) and an average length of 15.45 bp, and were the most abundant SSRs class, while the percentages of di-, tetra-, tri-, penta-, and hexa-nucleotide repeats were 22.86%, 11.98%, 11.58%, 6.65%, and 0.19%, respectively. Different classes of SSRs varied in their repeat number, with the highest being 42 for dinucleotides. Results reveal that repeat categories A, AC, AT, AAC, AGC, GTTT, CTTT, ATTT, and AACTG predominate on the Y chromosome. This study provides insight into the organization of cattle Y chromosome repetitive DNA, as well as information useful for developing more polymorphic cattle Y-chromosome-specific SSRs.

Detecting the existence of gene flow between Spanish and North African goats through a coalescent approach

Human-driven migrations are one of the main processes shaping the genetic diversity and population structure of domestic species. However, their magnitude and direction have been rarely analysed in a statistical framework. We aimed to estimate the impact of migration on the population structure of Spanish and African goats. To achieve this goal, we analysed a dataset of 1,472 individuals typed with 23 microsatellites. Population structure of African and Spanish goats was moderate (mean FST = 0.07), with the exception of the Canarian and South African breeds that displayed a significant differentiation when compared to goats from North Africa and Nigeria. Measurement of gene flow with Migrate-n and IMa coalescent genealogy samplers supported the existence of a bidirectional gene flow between African and Spanish goats. Moreover, IMa estimates of the effective number of migrants were remarkably lower than those calculated with Migrate-n and classical approaches. Such discrepancies suggest that recent divergence, rather than extensive gene flow, is the main cause of the weak population structure observed in caprine breeds.

Genetic origin, admixture and population history of aurochs (Bos primigenius) and primitive European cattle

The domestication of taurine cattle initiated ~10 000 years ago in the Near East from a wild aurochs (Bos primigenius) population followed by their dispersal through migration of agriculturalists to Europe. Although gene flow from wild aurochs still present at the time of this early dispersion is still debated, some of the extant primitive cattle populations are believed to possess the aurochs-like primitive features. In this study, we use genome-wide single nucleotide polymorphisms to assess relationship, admixture patterns and demographic history of an ancient aurochs sample and European cattle populations, several of which have primitive features and are suitable for extensive management. The principal component analysis, the model-based clustering and a distance-based network analysis support previous works suggesting different histories for north-western and southern European cattle. Population admixture analysis indicates a zebu gene flow in the Balkan and Italian Podolic cattle populations. Our analysis supports the previous report of gene flow between British and Irish primitive cattle populations and local aurochs. In addition, we show evidence of aurochs gene flow in the Iberian cattle populations indicating wide geographical distribution of the aurochs. Runs of homozygosity (ROH) reveal that demographic processes like genetic isolation and breed formation have contributed to genomic variations of European cattle populations. The ROH also indicate recent inbreeding in southern European cattle populations. We conclude that in addition to factors such as ancient human migrations, isolation by distance and cross-breeding, gene flow between domestic and wild-cattle populations also has shaped genomic composition of European cattle populations.

A genome-wide perspective about the diversity and demographic history of seven Spanish goat breeds

Background

The main goal of the current work was to infer the demographic history of seven Spanish goat breeds (Malagueña, Murciano-Granadina, Florida, Palmera, Mallorquina, Bermeya and Blanca de Rasquera) based on genome-wide diversity data generated with the Illumina Goat SNP50 BeadChip (population size, N = 176). Five additional populations from Europe (Saanen and Carpathian) and Africa (Tunisian, Djallonké and Sahel) were also included in this analysis (N = 80) for comparative purposes.

Results

Our results show that the genetic background of Spanish goats traces back mainly to European breeds although signs of North African admixture were detected in two Andalusian breeds (Malagueña and Murciano-Granadina). In general, observed and expected heterozygosities were quite similar across the seven Spanish goat breeds under analysis irrespective of their population size and conservation status. For the Mallorquina and Blanca de Rasquera breeds, which have suffered strong population declines during the past decades, we observed increased frequencies of large-sized (ROH), a finding that is consistent with recent inbreeding. In contrast, a substantial part of the genome of the Palmera goat breed comprised short ROH, which suggests a strong and ancient founder effect.

Conclusions

Admixture with African goats, genetic drift and inbreeding have had different effects across the seven Spanish goat breeds analysed in the current work. This has generated distinct patterns of genome-wide diversity that provide new clues about the demographic history of these populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0229-6) contains supplementary material, which is available to authorized users.

Mitogenomes from Egyptian Cattle Breeds: New Clues on the Origin of Haplogroup Q and the Early Spread of Bos taurus from the Near East

BACKGROUND

Genetic studies support the scenario that Bos taurus domestication occurred in the Near East during the Neolithic transition about 10 thousand years (ky) ago, with the likely exception of a minor secondary event in Italy. However, despite the proven effectiveness of whole mitochondrial genome data in providing valuable information concerning the origin of taurine cattle, until now no population surveys have been carried out at the level of mitogenomes in local breeds from the Near East or surrounding areas. Egypt is in close geographic and cultural proximity to the Near East, in particular the Nile Delta region, and was one of the first neighboring areas to adopt the Neolithic package. Thus, a survey of mitogenome variation of autochthonous taurine breeds from the Nile Delta region might provide new insights on the early spread of cattle rearing outside the Near East.

METHODOLOGY

Using Illumina high-throughput sequencing we characterized the mitogenomes from two cattle breeds, Menofi (N = 17) and Domiaty (N = 14), from the Nile Delta region. Phylogenetic and Bayesian analyses were subsequently performed.

CONCLUSIONS

Phylogenetic analyses of the 31 mitogenomes confirmed the prevalence of haplogroup T1, similar to most African cattle breeds, but showed also high frequencies for haplogroups T2, T3 and Q1, and an extremely high haplotype diversity, while Bayesian skyline plots pointed to a main episode of population growth ~12.5 ky ago. Comparisons of Nile Delta mitogenomes with those from other geographic areas revealed that (i) most Egyptian mtDNAs are probably direct local derivatives from the founder domestic herds which first arrived from the Near East and the extent of gene flow from and towards the Nile Delta region was limited after the initial founding event(s); (ii) haplogroup Q1 was among these founders, thus proving that it underwent domestication in the Near East together with the founders of the T clades.

Genome-wide insights into population structure and genetic history of Tunisian local cattle using the illumina bovinesnp50 beadchip

Background

Tunisian local cattle populations are at risk of extinction as they were massively crossed with imported breeds. Preservation of indigenous livestock populations is important because each of them comprises a unique set of genes resulting from a local environment-driven selection that occurred over hundreds of years. The diversity and genetic structure of Tunisian local cattle populations are poorly understood. However, such information is crucial to the conservation and sustainable use of genetic resources.

In addition, comparing the genomic structure of population sets from different parts of the world could help yield insight into their origin and history.

In the present study, we provide a detailed assessment of the population structure of the three Tunisian local cattle populations using various methods, and we highlight their origin and history by investigating approximately ~38,000 SNPs in a broad panel of 878 individuals from 37 worldwide cattle breeds representative of African, European and indicine populations.

Results

Our study revealed a low level of divergence and high genetic diversity in Tunisian local cattle reflecting low levels of genetic drift. A Comparison with the worldwide cattle panel pinpointed the admixed origin of the genome of the three Tunisian populations with the two main European and African ancestries. Our results were in agreement with previous historical and archaeological reports about the past gene flow that existed between North African and South European breeds, in particular with Iberian cattle. We also detected a low-level indicine introgression in the three Tunisian populations and we inferred that indicine ancestry was inherited via African ancestors.

Conclusions

Our results represent the first study providing genetic evidence about the origin and history of Tunisian local cattle. The information provided by the fine-scale genetic characterization of our study will enhance the establishment of a national conservation strategy for these populations. These results may enable the identification of genetic variants involved in adaptation to harsh environmental conditions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1638-6) contains supplementary material, which is available to authorized users.

The genetic prehistory of domesticated cattle from their origin to the spread across Europe

Background

Cattle domestication started in the 9^th millennium BC in Southwest Asia. Domesticated cattle were then introduced into Europe during the Neolithic transition. However, the scarcity of palaeogenetic data from the first European domesticated cattle still inhibits the accurate reconstruction of their early demography. In this study, mitochondrial DNA from 193 ancient and 597 modern domesticated cattle (Bos taurus) from sites across Europe, Western Anatolia and Iran were analysed to provide insight into the Neolithic dispersal process and the role of the local European aurochs population during cattle domestication.

Results

Using descriptive summary statistics and serial coalescent simulations paired with approximate Bayesian computation we find: (i) decreasing genetic diversity in a southeast to northwest direction, (ii) strong correlation of genetic and geographical distances, iii) an estimated effective size of the Near Eastern female founder population of 81, iv) that the expansion of cattle from the Near East and Anatolia into Europe does not appear to constitute a significant bottleneck, and that v) there is evidence for gene-flow between the Near Eastern/Anatolian and European cattle populations in the early phases of the European Neolithic, but that it is restricted after 5,000 BCE.

Conclusions

The most plausible scenario to explain these results is a single and regionally restricted domestication process of cattle in the Near East with subsequent migration into Europe during the Neolithic transition without significant maternal interbreeding with the endogenous wild stock. Evidence for gene-flow between cattle populations from Southwestern Asia and Europe during the earlier phases of the European Neolithic points towards intercontinental trade connections between Neolithic farmers.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0203-2) contains supplementary material, which is available to authorized users.

On the History of Cattle Genetic Resources

Cattle are our most important livestock species because of their production and role in human culture. Many breeds that differ in appearance, performance and environmental adaptation are kept on all inhabited continents, but the historic origin of the diverse phenotypes is not always clear. We give an account of the history of cattle by integrating archaeological record and pictorial or written sources, scarce until 300 years ago, with the recent contributions of DNA analysis. We describe the domestication of their wild ancestor, migrations to eventually all inhabited continents, the developments during prehistory, the antiquity and the Middle Ages, the relatively recent breed formation, the industrial cattle husbandry in the Old and New World and the current efforts to preserve the cattle genetic resources. Surveying the available information, we propose three main and overlapping phases during the development of the present genetic diversity: (i) domestication and subsequent wild introgression; (ii) natural adaptation to a diverse agricultural habitat; and (iii) breed development.

Over the sands and far away: interpreting an Iberian mitochondrial lineage with ancient Western African origins

OBJECTIVES

There is an ongoing effort to characterize the genetic links between Africa and Europe, mostly using lineages and haplotypes that are specific to one continent but had an ancient origin in the other. Mitochondrial DNA has been proven to be a very useful tool for this purpose since a high number of putatively European-specific variants of the African L* lineages have been defined over the years. Due to their geographic locations, Spain and Portugal seem to be ideal places for searching for these lineages.

METHODS

Five members of a minor branch of haplogroup L3f were found in recent DNA samplings in the region of Asturias (Northern Spain), which is known for its historical isolation. The frequency of L3f in this population (≈1%) is unexpectedly high in comparison with other related lineages in Europe. Complete mitochondrial DNA sequencing of these L3f lineages, as well phylogenetic and phylogeographic comparative analyses have been performed.

RESULTS

The L3f variant found in Asturias seems to constitute an Iberian-specific haplogroup, distantly related to lineages in Northern Africa and with a deep ancestry in Western Africa. Coalescent algorithms estimate the minimum arrival time as 8,000 years ago, and a possible route through the Gibraltar Strait.

CONCLUSIONS

Results are concordant with a previously proposed Neolithic connection between Southern Europe and Western Africa, which might be key to the proper understanding of the ancient links between these two continents.

A Genomic Approach for Distinguishing between Recent and Ancient Admixture as Applied to Cattle

Genomic data facilitate opportunities to track complex population histories of divergence and gene flow. We developed a metric, scaled block size (SBS), which uses the nonrecombined block size of introgressed regions of chromosomes to differentiate between recent and ancient types of admixture, and applied it to the reconstruction of admixture in cattle. Cattle are descendants of 2 independently domesticated lineages, taurine and indicine, which diverged more than 200 000 years ago. Several breeds have hybrid ancestry between these divergent lineages. Using 47 506 single-nucleotide polymorphisms, we analyzed the genomic architecture of the ancestry of 1369 individuals. We focused on 4 groups with admixed ancestry, including 2 anciently admixed African breeds (n = 58; n = 43), New World cattle of Spanish origin (n = 51), and known recent hybrids (n = 46). We estimated the ancestry of chromosomal regions for each individual and used the SBS metric to differentiate the timing of admixture among groups and among individuals within groups. By comparing SBS values of test individuals with standards with known recent hybrid ancestry, we were able to differentiate individuals of recent hybrid origin from other admixed cattle. We also estimated ancestry at the chromosomal scale. The X chromosome exhibits reduced indicine ancestry in recent hybrid, New World, and western African cattle, with virtually no evidence of indicine ancestry in New World cattle.

Direct estimates of natural selection in Iberia indicate calcium absorption was not the only driver of lactase persistence in Europe

Lactase persistence (LP) is a genetically determined trait whereby the enzyme lactase is expressed throughout adult life. Lactase is necessary for the digestion of lactose--the main carbohydrate in milk--and its production is downregulated after the weaning period in most humans and all other mammals studied. Several sources of evidence indicate that LP has evolved independently, in different parts of the world over the last 10,000 years, and has been subject to strong natural selection in dairying populations. In Europeans, LP is strongly associated with, and probably caused by, a single C to T mutation 13,910 bp upstream of the lactase (LCT) gene (-13,910*T). Despite a considerable body of research, the reasons why LP should provide such a strong selective advantage remain poorly understood. In this study, we examine one of the most widely cited hypotheses for selection on LP--that fresh milk consumption supplemented the poor vitamin D and calcium status of northern Europe's early farmers (the calcium assimilation hypothesis). We do this by testing for natural selection on -13,910*T using ancient DNA data from the skeletal remains of eight late Neolithic Iberian individuals, whom we would not expect to have poor vitamin D and calcium status because of relatively high incident UVB light levels. None of the eight samples successfully typed in the study had the derived T-allele. In addition, we reanalyze published data from French Neolithic remains to both test for population continuity and further examine the evolution of LP in the region. Using simulations that accommodate genetic drift, natural selection, uncertainty in calibrated radiocarbon dates, and sampling error, we find that natural selection is still required to explain the observed increase in allele frequency. We conclude that the calcium assimilation hypothesis is insufficient to explain the spread of LP in Europe.

Ancient mtDNA analysis of early 16(th) century Caribbean cattle provides insight into founding populations of New World creole cattle breeds

The Columbian Exchange resulted in a widespread movement of humans, plants and animals between the Old and New Worlds. The late 15^th to early 16^th century transfer of cattle from the Iberian Peninsula and Canary Islands to the Caribbean laid the foundation for the development of American creole cattle (Bos taurus) breeds. Genetic analyses of modern cattle from the Americas reveal a mixed ancestry of European, African and Indian origins. Recent debate in the genetic literature centers on the ‘African’ haplogroup T1 and its subhaplogroups, alternatively tying their origins to the initial Spanish herds, and/or from subsequent movements of taurine cattle through the African slave trade. We examine this problem through ancient DNA analysis of early 16^th century cattle bone from Sevilla la Nueva, the first Spanish colony in Jamaica. In spite of poor DNA preservation, both T3 and T1 haplogroups were identified in the cattle remains, confirming the presence of T1 in the earliest Spanish herds. The absence, however, of “African-derived American” haplotypes (AA/T1c1a1) in the Sevilla la Nueva sample, leaves open the origins of this sub-haplogroup in contemporary Caribbean cattle.

New World cattle show ancestry from multiple independent domestication events

Previous archeological and genetic research has shown that modern cattle breeds are descended from multiple independent domestication events of the wild aurochs (Bos primigenius) ∼10,000 y ago. Two primary areas of domestication in the Middle East/Europe and the Indian subcontinent resulted in taurine and indicine lines of cattle, respectively. American descendants of cattle brought by European explorers to the New World beginning in 1493 generally have been considered to belong to the taurine lineage. Our analyses of 47,506 single nucleotide polymorphisms show that these New World cattle breeds, as well as many related breeds of cattle in southern Europe, actually exhibit ancestry from both the taurine and indicine lineages. In this study, we show that, although European cattle are largely descended from the taurine lineage, gene flow from African cattle (partially of indicine origin) contributed substantial genomic components to both southern European cattle breeds and their New World descendants. New World cattle breeds, such as Texas Longhorns, provide an opportunity to study global population structure and domestication in cattle. Following their introduction into the Americas in the late 1400s, semiferal herds of cattle underwent between 80 and 200 generations of predominantly natural selection, as opposed to the human-mediated artificial selection of Old World breeding programs. Our analyses of global cattle breed population history show that the hybrid ancestry of New World breeds contributed genetic variation that likely facilitated the adaptation of these breeds to a novel environment.

Ancient DNA studies: new perspectives on old samples

In spite of past controversies, the field of ancient DNA is now a reliable research area due to recent methodological improvements. A series of recent large-scale studies have revealed the true potential of ancient DNA samples to study the processes of evolution and to test models and assumptions commonly used to reconstruct patterns of evolution and to analyze population genetics and palaeoecological changes. Recent advances in DNA technologies, such as next-generation sequencing make it possible to recover DNA information from archaeological and paleontological remains allowing us to go back in time and study the genetic relationships between extinct organisms and their contemporary relatives. With the next-generation sequencing methodologies, DNA sequences can be retrieved even from samples (for example human remains) for which the technical pitfalls of classical methodologies required stringent criteria to guaranty the reliability of the results. In this paper, we review the methodologies applied to ancient DNA analysis and the perspectives that next-generation sequencing applications provide in this field.

Ancient DNA extracted from Danish aurochs (Bos primigenius): genetic diversity and preservation

We extracted DNA from 39 Danish aurochs specimens and successfully amplified and sequenced a 252 base pair long fragment of the multivariable region I of the mitochondrial control region from 11 specimens. The sequences from these specimens dated back to 9830-2865 14Cyr BP and represent the first study of genetic variation of Danish aurochs. In addition, for all specimens we address correlations between the ability to obtain DNA sequences and various parameters such as the age of the sample, the collagen content, the museum storage period, Danish geography and whether the specimens were found in an archeological or geological context. We find that aurochs from southern Scandinavia display a star-shaped population genetic structure, that is indicative of a local and relatively recent diversification from a few ancestral haplotypes that may have originated in the ancestral Western European population before migration northwards during the retreat of the glaciers. Scenarios suggesting several invasions of genetically distinct aurochs are not supported by these analyses. Rather, our results suggest that a single continuous migration northward occurred. Our findings also suggest, although with only limited support, that aurochs in Northwestern Europe underwent a population expansion beginning shortly after the retreat of the glacial ice from Denmark and had a stable population size until the population decline that must have occurred prior to extinction. The absence of haplotypes similar to modern domestic cattle in our aurochs suggests that introgression between these species must have been limited, if it occurred at all. We found that the successful recovery of genetic material for PCR amplification correlates with sample age and local geographic conditions. However, contrary to other studies, we found no significant correlation between length of time in museum storage or the type of the locality in which a specimen was discovered (archeological or geological) and amplification success. Finally, we found large variances in our estimates of collagen content preventing an evaluation of this as an indicator of preservation quality.

Conservation genetics of cattle, sheep, and goats

Cattle, sheep and goats were domesticated about 10,000 years ago, spread out of the domestication centers in Europe, Asia, and Africa during the next few thousands years, and gave many populations locally adapted. After a very long period of soft selection, the situation changed dramatically 200 years ago with the emergence of the breed concept. The selection pressure strongly increased, and the reproduction among breeds was seriously reduced, leading to the fragmentation of the initial gene pool. More recently, the selection pressure was increased again via the use of artificial insemination, leading to a few industrial breeds with very high performances, but with low effective population sizes. Beside this performance improvement of industrial breeds, genetic resources are being lost, because of the replacement of traditional breeds by high performance industrial breeds at the worldwide level, and because of the loss of genetic diversity in these industrial breeds. Many breeds are already extinct, and genetic resources in cattle, sheep, and goats are thus highly endangered, particularly in developed countries. The recent development of next generation sequencing technologies opens new avenues for properly characterizing the genetic resources, not only in the very diverse domestic breeds, but also in their wild relatives. Based on sound genetic characterization, urgent conservation measures must be taken to avoid an irremediable loss of farm animal genetic resources, integrating economical, sociological, and political parameters.

The complete mitochondrial genome of an 11,450-year-old aurochsen (Bos primigenius) from Central Italy

Background

Bos primigenius, the aurochs, is the wild ancestor of modern cattle breeds and was formerly widespread across Eurasia and northern Africa. After a progressive decline, the species became extinct in 1627. The origin of modern taurine breeds in Europe is debated. Archaeological and early genetic evidence point to a single Near Eastern origin and a subsequent spread during the diffusion of herding and farming. More recent genetic data are instead compatible with local domestication events or at least some level of local introgression from the aurochs. Here we present the analysis of the complete mitochondrial genome of a pre-Neolithic Italian aurochs.

Results

In this study, we applied a combined strategy employing both multiplex PCR amplifications and 454 pyrosequencing technology to sequence the complete mitochondrial genome of an 11,450-year-old aurochs specimen from Central Italy. Phylogenetic analysis of the aurochs mtDNA genome supports the conclusions from previous studies of short mtDNA fragments - namely that Italian aurochsen were genetically very similar to modern cattle breeds, but highly divergent from the North-Central European aurochsen.

Conclusions

Complete mitochondrial genome sequences are now available for several modern cattle and two pre-Neolithic mtDNA genomes from very different geographic areas. These data suggest that previously identified sub-groups within the widespread modern cattle mitochondrial T clade are polyphyletic, and they support the hypothesis that modern European breeds have multiple geographic origins.

Mitochondrial DNA patterns in the Iberian Northern plateau: population dynamics and substructure of the Zamora province

Several studies have shown the importance of recent events in the configuration of the genetic landscape of a specific territory. In this context, due to the phenomena of repopulation and demographic fluctuations that took place in recent centuries, the Iberian Northern plateau is a very interesting case study. The main aim of this work is to check if recent population movements together with existing boundaries (geographical and administrative) have influenced the current genetic composition of the area. To accomplish this general purpose, mitochondrial DNA variations of 214 individuals from a population located in the Western region of the Iberian Northern plateau (the province of Zamora) were analyzed. Results showed a typical Western European mitochondrial DNA haplogroup composition. However, unexpected high frequencies of U5, HV0, and L haplogroups were found in some regions. The analyses of microdifferentiation showed that there are differences between regions, but no geographic substructure organization can be noticed. It can be stated that the differences observed in the genetic pool of the sampled area at regional level results from the mixture of different populations carrying new lineages into this area at different points in history.

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.

High frequency of lactose intolerance in a prehistoric hunter-gatherer population in northern Europe

Background

Genes and culture are believed to interact, but it has been difficult to find direct evidence for the process. One candidate example that has been put forward is lactase persistence in adulthood, i.e. the ability to continue digesting the milk sugar lactose after childhood, facilitating the consumption of raw milk. This genetic trait is believed to have evolved within a short time period and to be related with the emergence of sedentary agriculture.

Results

Here we investigate the frequency of an allele (-13910*T) associated with lactase persistence in a Neolithic Scandinavian population. From the 14 individuals originally examined, 10 yielded reliable results. We find that the T allele frequency was very low (5%) in this Middle Neolithic hunter-gatherer population, and that the frequency is dramatically different from the extant Swedish population (74%).

Conclusions

We conclude that this difference in frequency could not have arisen by genetic drift and is either due to selection or, more likely, replacement of hunter-gatherer populations by sedentary agriculturalists.

Analysis of STR markers reveals high genetic structure in Portuguese native cattle

Genetic structure and diversity of 13 Portuguese native and 3 imported cattle breeds were assessed with 39 microsatellites. Allelic richness per locus was high, with an overall average of 8.3 +/- 2.5. The mean observed and expected heterozygosities were 0.673 +/- 0.043 and 0.691 +/- 0.034, respectively. The mean number of alleles per breed ranged between 5.36 +/- 1.27 and 7.87 +/- 2.66. Brava de Lide and Mirandesa breeds had the lowest genetic diversity, whereas Minhota, Arouquesa, and Mertolenga had the highest. Significant (P < 0.05) heterozygote deficit was detected in all breeds except Garvonesa, Marinhoa, Minhota, and Limousin. Hardy-Weinberg deviations are most probably due to inbreeding, particularly in Alentejana, Brava de Lide, Mertolenga, and Ramo Grande (F(is) > 0, P < 0.0001). Based on the principal component and the Neighbor-Net analyses, Mirandesa was the most genetically distinct breed. Even though admixture was detected across all breeds (6.7%, q < 0.800), the molecular structure was consistent with original breed designations, with the exception of Cachena that had a clear influence of Barrosã (K = 15). Mertolenga showed substructure with independent clustering of red speckled animals. The percentage animals correctly assigned was >or=90 in all breeds except Cachena, Garvonesa, and Preta (q >or= 0.800). The results obtained here confirmed that high levels of genetic diversity exist within Portuguese native cattle and that the breeds are highly structured. Conservation measures should be implemented for all native breeds to minimize inbreeding.

Origins and genetic diversity of New World Creole cattle: inferences from mitochondrial and Y chromosome polymorphisms

The ancestry of New World cattle was investigated through the analysis of mitochondrial and Y chromosome variation in Creoles from Argentina, Brazil, Mexico, Paraguay and the United States of America. Breeds that influenced the Creoles, such as Iberian native, British and Zebu, were also studied. Creoles showed high mtDNA diversity (H = 0.984 +/- 0.003) with a total of 78 haplotypes, and the European T3 matriline was the most common (72.1%). The African T1a haplogroup was detected (14.6%), as well as the ancestral African-derived AA matriline (11.9%), which was absent in the Iberian breeds. Genetic proximity among Creoles, Iberian and Atlantic Islands breeds was inferred through their sharing of mtDNA haplotypes. Y-haplotype diversity in Creoles was high (H = 0.779 +/- 0.019), with several Y1, Y2 and Y3 haplotypes represented. Iberian patrilines in Creoles were more difficult to infer and were reflected by the presence of H3Y1 and H6Y2. Y-haplotypes confirmed crossbreeding with British cattle, mainly of Hereford with Pampa Chaqueño and Texas Longhorn. Male-mediated Bos indicus introgression into Creoles was found in all populations, except Argentino1 (herd book registered) and Pampa Chaqueño. The detection of the distinct H22Y3 patriline with the INRA189-90 allele in Caracú suggests introduction of bulls directly from West Africa. Further studies of Spanish and African breeds are necessary to elucidate the origins of Creole cattle, and determine the exact source of their African lineages.

The multifaceted origin of taurine cattle reflected by the mitochondrial genome

A Neolithic domestication of taurine cattle in the Fertile Crescent from local aurochsen (Bos primigenius) is generally accepted, but a genetic contribution from European aurochsen has been proposed. Here we performed a survey of a large number of taurine cattle mitochondrial DNA (mtDNA) control regions from numerous European breeds confirming the overall clustering within haplogroups (T1, T2 and T3) of Near Eastern ancestry, but also identifying eight mtDNAs (1.3%) that did not fit in haplogroup T. Sequencing of the entire mitochondrial genome showed that four mtDNAs formed a novel branch (haplogroup R) which, after the deep bifurcation that gave rise to the taurine and zebuine lineages, constitutes the earliest known split in the mtDNA phylogeny of B. primigenius. The remaining four mtDNAs were members of the recently discovered haplogroup Q. Phylogeographic data indicate that R mtDNAs were derived from female European aurochsen, possibly in the Italian Peninsula, and sporadically included in domestic herds. In contrast, the available data suggest that Q mtDNAs and T subclades were involved in the same Neolithic event of domestication in the Near East. Thus, the existence of novel (and rare) taurine haplogroups highlights a multifaceted genetic legacy from distinct B. primigenius populations. Taking into account that the maternally transmitted mtDNA tends to underestimate the extent of gene flow from European aurochsen, the detection of the R mtDNAs in autochthonous breeds, some of which are endangered, identifies an unexpected reservoir of genetic variation that should be carefully preserved.

The genetic legacy of religious diversity and intolerance: paternal lineages of Christians, Jews, and Muslims in the Iberian Peninsula

Most studies of European genetic diversity have focused on large-scale variation and interpretations based on events in prehistory, but migrations and invasions in historical times could also have had profound effects on the genetic landscape. The Iberian Peninsula provides a suitable region for examination of the demographic impact of such recent events, because its complex recent history has involved the long-term residence of two very different populations with distinct geographical origins and their own particular cultural and religious characteristics—North African Muslims and Sephardic Jews. To address this issue, we analyzed Y chromosome haplotypes, which provide the necessary phylogeographic resolution, in 1140 males from the Iberian Peninsula and Balearic Islands. Admixture analysis based on binary and Y-STR haplotypes indicates a high mean proportion of ancestry from North African (10.6%) and Sephardic Jewish (19.8%) sources. Despite alternative possible sources for lineages ascribed a Sephardic Jewish origin, these proportions attest to a high level of religious conversion (whether voluntary or enforced), driven by historical episodes of social and religious intolerance, that ultimately led to the integration of descendants. In agreement with the historical record, analysis of haplotype sharing and diversity within specific haplogroups suggests that the Sephardic Jewish component is the more ancient. The geographical distribution of North African ancestry in the peninsula does not reflect the initial colonization and subsequent withdrawal and is likely to result from later enforced population movement—more marked in some regions than in others—plus the effects of genetic drift.

Y chromosome haplotype analysis in Portuguese cattle breeds using SNPs and STRs

DNA samples from 307 males of 13 Portuguese native cattle breeds, 57 males of the 3 major exotic breeds in Portugal (Charolais, Friesian, and Limousin), and 5 Brahman (Bos indicus) were tested for 5 single nucleotide polymorphisms, 1 "indel," and 7 microsatellites specific to the Y chromosome. The 13 Y-haplotypes defined included 3 previously described patrilines (Y1, Y2, and Y3) and 10 new haplotypes within Bos taurus. Native cattle contained most of the diversity with 7 haplotypes (H2Y1, H3Y1, H5Y1, H7Y2, H8Y2, H10Y2, and H12Y2) found only in these breeds. H6Y2 and H11Y2 occurred in high frequency across breeds including the exotics. Introgression of Friesian cattle into Ramo Grande was inferred through their sharing of haplotype H4Y1. Among the native breeds, Mertolenga had the highest haplotype diversity (0.68 +/- 0.07), Brava de Lide was the least differentiated. The analyses of molecular variance showed significant (P < 0.0001) differences between breeds with more than 64% of the total genetic variation found among breeds within groups and 33-35% within breeds. The detection of INRA189-104 allele in 8 native breeds suggested influence of African cattle in breeds of the Iberian Peninsula. The presence in Portuguese breeds of Y1 patrilines, also found in aurochs, could represent more ancient local haplotypes.

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.