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

History and genetic diversity of African sheep: Contrasting phenotypic and genomic diversity

Domesticated sheep have adapted to contrasting and extreme environments and continue to play important roles in local community-based economies throughout Africa. Here we review the Neolithic migrations of thin-tailed sheep and the later introductions of fat-tailed sheep into eastern Africa. According to contemporary pictorial evidence, the latter occurred in Egypt not before the Ptolemaic period (305-25 BCE). We further describe the more recent history of sheep in Egypt, the Maghreb, west and central Africa, central-east Africa, and southern Africa. We also present a comprehensive molecular survey based on the analysis of 50 K SNP genotypes for 59 African breeds contributed by several laboratories. We propose that gene flow and import of fat-tailed sheep have partially overwritten the diversity profile created by the initial migration. We found a genetic contrast between sheep north and south of the Sahara and a west-east contrast of thin- and fat-tailed sheep. There is no close relationship between African and central and east Asian fat-tailed breeds, whereas we observe within Africa only a modest effect of tail types on breed relationships.

The African introgression of Murciano Granadina goats has a Moroccan origin and displays remarkable levels of inter-individual variability

There is evidence that Murciano Granadina (MG), the most important caprine dairy breed in Spain, has been introgressed by African goats, but the precise geographic origin of such introgression has not been identified yet. Moreover, an accurate estimate of the magnitude of this African introgression is lacking, since current estimates are based on small numbers of sampled individuals. The aim of our work was to tackle these two issues by genotyping 500 MG goats with the Goat SNP50 BeadChip and comparing their genotypes with those of reference populations from Spain (Bermeya), France (Saanen), Morocco (Barcha, Draa, Ghazalia, Noire de Atlas, Nord, Moroccan), Egypt (Barki, Oasis, Saidi), Algeria (Arabia, Makatia, M'Zabite, Kabyle), Tunisia (Tunisian native breeds) and Sudan (Desert, Nilotic, Taggar). The population of 500 MG goats was subdivided into 10 datasets of 50 individuals to ensure that sample sizes of the target (MG) and reference populations are balanced. Performance of an unsupervised ADMIXTURE analysis demonstrated that MG goats have a North African ancestry, with an average proportion of 4.4 ± 2.3%. Next, we did a supervised ADMIXTURE analysis that revealed that the Moroccan genetic component reaches a proportion of 4.01 ± 3.9% in MG goats, while the Algerian (0.001 ± 0.001%), Egyptian (0.2 ± 0.1%), Sudanese (0.1 ± 0.1%) and Tunisian (0.3 ± 0.4%) components are present in extremely small proportions. The historical circumstances of this introgression event are currently unknown, but several plausible scenarios are outlined. Moreover, our results show considerable inter-individual heterogeneity regarding the magnitude of the Moroccan introgression of MG goats (0%- 12% depending on the MG data set under analysis). This result implies that reliable estimates about the introgression of autochthonous livestock by exotic breeds can only be obtained by extensively sampling target populations.

Investigating the footprint of post-domestication dispersal on the diversity of modern European, African and Asian goats

BACKGROUND

Goats were domesticated in the Fertile Crescent about 10,000 years before present (YBP) and subsequently spread across Eurasia and Africa. This dispersal is expected to generate a gradient of declining genetic diversity with increasing distance from the areas of early livestock management. Previous studies have reported the existence of such genetic cline in European goat populations, but they were based on a limited number of microsatellite markers. Here, we have analyzed data generated by the AdaptMap project and other studies. More specifically, we have used the geographic coordinates and estimates of the observed (Ho) and expected (He) heterozygosities of 1077 European, 1187 African and 617 Asian goats belonging to 38, 43 and 22 different breeds, respectively, to find out whether genetic diversity and distance to Ganj Dareh, a Neolithic settlement in western Iran for which evidence of an early management of domestic goats has been obtained, are significantly correlated.

RESULTS

Principal component and ADMIXTURE analyses revealed an incomplete regional differentiation of European breeds, but two genetic clusters representing Northern Europe and the British-Irish Isles were remarkably differentiated from the remaining European populations. In African breeds, we observed five main clusters: (1) North Africa, (2) West Africa, (3) East Africa, (4) South Africa, and (5) Madagascar. Regarding Asian breeds, three well differentiated West Asian, South Asian and East Asian groups were observed. For European and Asian goats, no strong evidence of significant correlations between Ho and He and distance to Ganj Dareh was found. In contrast, in African breeds we detected a significant gradient of diversity, which decreased with distance to Ganj Dareh.

CONCLUSIONS

The detection of a genetic cline associated with distance to the Ganj Dareh in African but not in European or Asian goat breeds might reflect differences in the post-domestication dispersal process and subsequent migratory movements associated with the management of caprine populations from these three continents.

The demographic history and adaptation of Canarian goat breeds to environmental conditions through the use of genome-wide SNP data

BACKGROUND

The presence of goats in the Canary Islands dates back to the late 1st millennium BC, which coincides with the colonization by the Amazigh settlers. However, the exact geographic origin of Canarian goats is uncertain since the Amazigh peoples were distributed over a wide spatial range. Nowadays, three Canarian breeds (Palmera, Majorera and Tinerfeña) are officially recognized, along with two distinct South and North Tinerfeña ecotypes, with the South Tinerfeña and Majorera goats thriving in arid and dry semi-desertic environments and the Palmera and North Tinerfeña goats are adapted to humid and temperate areas that are influenced by trade winds. Genotypes for 224 Canarian goats were generated using the Illumina Goat single nucleotide polymorphism (SNP)50 BeadChip. By merging these data with the genotypes from 1007 individuals of African and Southern European ancestry, our aim was to ascertain the geographic origin of the Canarian goats and identify genes associated with adaptation to diverse environmental conditions.

RESULTS

The diversity indices of the Canarian breeds align with most of those of the analyzed local breeds from Africa and Europe, except for the Palmera goats that showed lower levels of genetic variation. The Canarian breeds demonstrate a significant genetic differentiation compared to other populations, which indicates a history of prolonged geographic isolation. Moreover, the phylogenetic reconstruction indicated that the ancestry of the Canarian goats is fundamentally North African rather than West African. The ADMIXTURE and the TreeMix analyses showed no evidence of gene flow between Canarian goats and other continental breeds. The analysis of runs of homozygosity (ROH) identified 13 ROH islands while the window-based FST method detected 25 genomic regions under selection. Major signals of selection were found on Capra hircus (CHI) chromosomes 6, 7, and 10 using various comparisons and methods.

CONCLUSIONS

This genome-wide analysis sheds new light on the evolutionary history of the four breeds that inhabit the Canary Islands. Our findings suggest a North African origin of the Canarian goats. In addition, within the genomic regions highlighted by the ROH and FST approaches, several genes related to body size and heat tolerance were identified.

Historical Westward Migration Phases of Ovis aries Inferred from the Population Structure and the Phylogeography of Occidental Mediterranean Native Sheep Breeds

In this study, the genetic relationship and the population structure of western Mediterranean basin native sheep breeds are investigated, analyzing Maghrebian, Central Italian, and Venetian sheep with a highly informative microsatellite markers panel. The phylogeographical analysis, between breeds’ differentiation level (Wright’s fixation index), gene flow, ancestral relatedness measured by molecular coancestry, genetic distances, divergence times estimates and structure analyses, were revealed based on the assessment of 975 genotyped animals. The results unveiled the past introduction and migration history of sheep in the occidental Mediterranean basin since the early Neolithic. Our findings provided a scenario of three westward sheep migration phases fitting properly to the westward Neolithic expansion argued by zooarcheological, historical and human genetic studies.

Characterizing the Mitochondrial Diversity of Arbi Goats from Tunisia

Arbi is one of the main local goat breeds in Tunisia, representing an important economic resource in arid and hot areas where cattle and sheep cannot thrive successfully. In the current work, we have characterized the mitochondrial diversity of 26 Arbi goats by partially sequencing the mitochondrial D-loop region. These sequences plus 10 retrieved from GenBank were analyzed with the DnaSP v.5.10.1, evidencing the existence of 12 different haplotypes. Nucleotide and haplotype diversities were 0.02 and 0.96. Moreover, median-joining network analysis showed that all D-loop sequences from Arbi goats correspond to haplogroup A and that in general they do not cluster with sequences from other goat breeds. The high diversity that has been observed in North African goats is compatible with the maritime diffusion of the Neolithic package 10,000-7000 YBP. Moreover, there are evidences that local Tunisian breeds have been extensively crossed with highly productive transboundary breeds in order to improve meat and milk yields. These uncontrolled crossing practices may lead to the loss of alleles that play key roles in the adaptation of Tunisian local breeds to a harsh environment.

Genetic diversity of the North African population revealed by the typing of SNPs in the DRD2/ANKK1 genomic region

The dopamine - related genes, like dopamine D2 receptor (DRD2) gene and ankyrin repeat and kinase domain containing 1 (ANKK1) gene are implicated in neurological functions. Some polymorphisms of the DRD2/ANKK1 locus (TaqIA, TaqIB, TaqID) have been used to study genetic diversity and the evolution of human populations. The present investigation aims to assess the genetic diversity in seven North African populations in order to explore their genetic structure and to compare them to others worldwide populations studied for the same locus. Nine single nucleotide polymorphisms (SNPs) from the DRD2/ANKK1 locus (rs1800497 TaqIA, rs2242592, rs1124492, rs6277, rs6275, rs1079727, rs2002453, rs2234690 and rs1079597 TaqIB) were typed in 366 individuals from seven North African populations: six from Tunisia (Sousse, Smar, Kesra, Kairouan, Mehdia and Kerkennah) and one from Libya. The allelic frequencies of rs2002453 and rs2234690 were higher in the Smar population than in the other North African populations. More, the Smar population showed the lowest average heterozygosity (0.313). The principal component analysis (PCA) showed that the Smar population was clearly separated from others. Furthermore, linkage disequilibrium analysis shown a high linkage disequilibrium in the North African population and essentially in Smar population. Comparison with other world populations has shown that the heterozygosity of North African population was very close to that of the African and European populations. The PCA and the haplotypic analysis suggested the presence of an important Eurasian genetic component for the North African population. These results suggested that the Smar population was isolated from the others North Africans ones by its peculiar genetic structure because of isolation, endogamy and genetic drift. On the other hand, the North African population is characterized by a multi ancestral gene pool from Eurasia and sub-Saharan Africa due to human migration since prehistoric times.

Dissection of ancestral genetic contributions to Creole goat populations

Goats have played a key role as source of nourishment for humans in their expansion all over the world in long land and sea trips. This has guaranteed a place for this species in the important and rapid episode of livestock expansion triggered by Columbus' arrival in the Americas in the late 1400s. The aims of this study are to provide a comprehensive perspective on genetic diversity in American goat populations and to assess their origins and evolutionary trajectories. This was achieved by combining data from autosomal neutral genetic markers obtained in more than two thousand samples that encompass a wide range of Iberian, African and Creole goat breeds. In general, even though Creole populations differ clearly from each other, they lack a strong geographical pattern of differentiation, such that populations of different admixed ancestry share relatively close locations throughout the large geographical range included in this study. Important Iberian signatures were detected in most Creole populations studied, and many of them, particularly the Cuban Creole, also revealed an important contribution of African breeds. On the other hand, the Brazilian breeds showed a particular genetic structure and were clearly separated from the other Creole populations, with some influence from Cape Verde goats. These results provide a comprehensive characterisation of the present structure of goat genetic diversity, and a dissection of the Iberian and African influences that gave origin to different Creole caprine breeds, disentangling an important part of their evolutionary history. Creole breeds constitute an important reservoir of genetic diversity that justifies the development of appropriate management systems aimed at improving performance without loss of genomic diversity.

Differential distribution of Y-chromosome haplotypes in Swiss and Southern European goat breeds

The analysis of Y-chromosome variation has provided valuable clues about the paternal history of domestic animal populations. The main goal of the current work was to characterize Y-chromosome diversity in 31 goat populations from Central Eastern (Switzerland and Romania) and Southern Europe (Spain and Italy) as well as in reference populations from Africa and the Near East. Towards this end, we have genotyped seven single nucleotide polymorphisms (SNPs), mapping to the SRY, ZFY, AMELY and DDX3Y Y-linked loci, in 275 bucks from 31 populations. We have observed a low level of variability in the goat Y-chromosome, with just five haplotypes segregating in the whole set of populations. We have also found that Swiss bucks carry exclusively Y1 haplotypes (Y1A: 24%, Y1B1: 15%, Y1B2: 43% and Y1C: 18%), while in Italian and Spanish bucks Y2A is the most abundant haplotype (77%). Interestingly, in Carpathian goats from Romania the Y2A haplotype is also frequent (42%). The high Y-chromosome differentiation between Swiss and Italian/Spanish breeds might be due to the post-domestication spread of two different Near Eastern genetic stocks through the Danubian and Mediterranean corridors. Historical gene flow between Southern European and Northern African goats might have also contributed to generate such pattern of genetic differentiation.

Goat domestication and breeding: a jigsaw of historical, biological and molecular data with missing pieces

Domestic goats (Capra hircus) are spread across the five continents with a census of 1 billion individuals. The worldwide population of goats descends from a limited number of bezoars (Capra aegagrus) domesticated 10 000 YBP (years before the present) in the Fertile Crescent. The extraordinary adaptability and hardiness of goats favoured their rapid spread over the Old World, reaching the Iberian Peninsula and Southern Africa 7000 YBP and 2000 YBP respectively. Molecular studies have revealed one major mitochondrial haplogroup A and five less frequent haplogroups B, C, D, F and G. Moreover, the analysis of autosomal and Y-chromosome markers has evidenced an appreciable geographic differentiation. The implementation of new molecular technologies, such as whole-genome sequencing and genome-wide genotyping, allows for the exploration of caprine diversity at an unprecedented scale, thus providing new insights into the evolutionary history of goats. In spite of a number of pitfalls, the characterization of the functional elements of the goat genome is expected to play a key role in understanding the genetic determination of economically relevant traits. Genomic selection and genome editing also hold great potential, particularly for improving traits that cannot be modified easily by traditional selection.