Multiple genomic solutions for local adaptation in two closely related species (sheep and goats) facing the same climatic constraints.

The question of how local adaptation takes place remains a fundamental question in evolutionary biology. The variation of allele frequencies in genes under selection over environmental gradients remains mainly theoretical and its empirical assessment would help understanding how adaptation happens over environmental clines. To bring new insights to this issue we set up a broad framework which aimed to compare the adaptive trajectories over environmental clines in two domesticated mammal species co-distributed in diversified landscapes. We sequenced the genomes of 160 sheep and 161 goats extensively managed along environmental gradients, including temperature, rainfall, seasonality and altitude, to identify genes and biological processes shaping local adaptation. Allele frequencies at putatively adaptive loci were rarely found to vary gradually along environmental gradients, but rather displayed a discontinuous shift at the extremities of environmental clines. Of the 430 candidate adaptive genes identified, only 6 were orthologous between sheep and goats and those responded differently to environmental pressures, suggesting different putative mechanisms involved in local adaptation in these two closely related species. Interestingly, the genomes of the 2 species were impacted differently by the environment, genes related to signatures of selection were most related to altitude, slope and rainfall seasonality for sheep, and summer temperature and spring rainfall for goats. The diversity of candidate adaptive pathways may result from a high number of biological functions involved in the adaptations to multiple eco-climatic gradients, and a differential role of climatic drivers on the two species, despite their co-distribution along the same environmental gradients. This study describes empirical examples of clinal variation in putatively adaptive alleles with different patterns in allele frequency distributions over continuous environmental gradients, thus showing the diversity of genetic responses in adaptive landscapes and opening new horizons for understanding genomics of adaptation in mammalian species and beyond.

Broad maternal geographic origin of domestic sheep in Anatolia and the Zagros.

We investigated the controversial origin of domestic sheep (Ovis aries) using large samples of contemporary and ancient domestic individuals and their closest wild relatives: the Asiatic mouflon (Ovis gmelini), the urial (Ovis vignei) and the argali (Ovis ammon). A phylogeny based on mitochondrial DNA, including 213 new cytochrome-b sequences of wild Ovism confirmed that O. gmelini is the maternal ancestor of sheep and precluded mtDNA contributions from O. vignei (and O. gmelini × O. vignei hybrids) to domestic lineages. We also produced 54 new control region sequences showing shared haplogroups (A, B, C and E) between domestic sheep and wild O. gmelini which localized the domestication center in eastern Anatolia and central Zagros, excluding regions further east where exclusively wild haplogroups were found. This overlaps with the geographic distribution of O. gmelini gmelini, further suggesting that the maternal origin of domestic sheep derives from this subspecies. Additionally, we produced 57 new CR sequences of Neolithic sheep remains from a large area covering Anatolia to Europe, showing the early presence of at least three mitochondrial haplogroups (A, B and D) in Western colonization routes. This confirmed that sheep domestication was a large-scale process that captured diverse maternal lineages (haplogroups).

Genomic Uniqueness of Local Sheep Breeds From Morocco.

Sheep farming is a major source of meat in Morocco and plays a key role in the country's agriculture. This study aims at characterizing the whole-genome diversity and demographic history of the main Moroccan sheep breeds, as well as to identify selection signatures within and between breeds. Whole genome data from 87 individuals representing the five predominant local breeds were used to estimate their level of neutral genetic diversity and to infer the variation of their effective population size over time. In addition, we used two methods to detect selection signatures: either for detecting selective sweeps within each breed separately or by detecting differentially selected regions by contrasting different breeds. We identified hundreds of genomic regions putatively under selection, which related to several biological terms involved in local adaptation or the expression of zootechnical performances such as Growth, UV protection, Cell maturation or Feeding behavior. The results of this study revealed selection signatures in genes that have an important role in traits of interest and increased our understanding of how genetic diversity is distributed in these local breeds. Thus, Moroccan local sheep breeds exhibit both a high genetic diversity and a large set of adaptive variations, and therefore, represent a valuable genetic resource for the conservation of sheep in the context of climate change.

VarGoats project: a dataset of 1159 whole-genome sequences to dissect Capra hircus global diversity.

BACKGROUND: Since their domestication 10,500 years ago, goat populations with distinctive genetic backgrounds have adapted to a broad variety of environments and breeding conditions. The VarGoats project is an international 1000-genome resequencing program designed to understand the consequences of domestication and breeding on the genetic diversity of domestic goats and to elucidate how speciation and hybridization have modeled the genomes of a set of species representative of the genus Capra. FINDINGS: A dataset comprising 652 sequenced goats and 507 public goat sequences, including 35 animals representing eight wild species, has been collected worldwide. We identified 74,274,427 single nucleotide polymorphisms (SNPs) and 13,607,850 insertion-deletions (InDels) by aligning these sequences to the latest version of the goat reference genome (ARS1). A Neighbor-joining tree based on Reynolds genetic distances showed that goats from Africa, Asia and Europe tend to group into independent clusters. Because goat breeds from Oceania and Caribbean (Creole) all derive from imported animals, they are distributed along the tree according to their ancestral geographic origin. CONCLUSIONS: We report on an unprecedented international effort to characterize the genome-wide diversity of domestic goats. This large range of sequenced individuals represents a unique opportunity to ascertain how the demographic and selection processes associated with post-domestication history have shaped the diversity of this species. Data generated for the project will also be extremely useful to identify deleterious mutations and polymorphisms with causal effects on complex traits, and thus will contribute to new knowledge that could be used in genomic prediction and genome-wide association studies.

Local adaptations of Mediterranean sheep and goats through an integrative approach.

Small ruminants are suited to a wide variety of habitats and thus represent promising study models for identifying genes underlying adaptations. Here, we considered local Mediterranean breeds of goats (n = 17) and sheep (n = 25) from Italy, France and Spain. Based on historical archives, we selected the breeds potentially most linked to a territory and defined their original cradle (i.e., the geographical area in which the breed has emerged), including transhumant pastoral areas. We then used the programs PCAdapt and LFMM to identify signatures of artificial and environmental selection. Considering cradles instead of current GPS coordinates resulted in a greater number of signatures identified by the LFMM analysis. The results, combined with a systematic literature review, revealed a set of genes with potentially key adaptive roles in relation to the gradient of aridity and altitude. Some of these genes have been previously implicated in lipid metabolism (SUCLG2, BMP2), hypoxia stress/lung function (BMPR2), seasonal patterns (SOX2, DPH6) or neuronal function (TRPC4, TRPC6). Selection signatures involving the PCDH9 and KLH1 genes, as well as NBEA/NBEAL1, were identified in both species and thus could play an important adaptive role.

Genetic Variations and Differential DNA Methylation to Face Contrasted Climates in Small Ruminants: An Analysis on Traditionally-Managed Sheep and Goats.

The way in which living organisms mobilize a combination of long-term adaptive mechanisms and short-term phenotypic plasticity to face environmental variations is still largely unknown. In the context of climate change, understanding the genetic and epigenetic bases for adaptation and plasticity is a major stake for preserving genomic resources and the resilience capacity of livestock populations. We characterized both epigenetic and genetic variations by contrasting 22 sheep and 21 goats from both sides of a climate gradient, focusing on free-ranging populations from Morocco. We produced for each individual Whole-Genome Sequence at 12X coverage and MeDIP-Seq data, to identify regions under selection and those differentially methylated. For both species, the analysis of genetic differences (FST) along the genome between animals from localities with high vs. low temperature annual variations detected candidate genes under selection in relation to environmental perception (5 genes), immunity (4 genes), reproduction (8 genes) and production (11 genes). Moreover, we found for each species one differentially methylated gene, namely AGPTA4 in goat and SLIT3 in sheep, which were both related, among other functions, to milk production and muscle development. In both sheep and goats, the comparison between genomic regions impacted by genetic and epigenetic variations suggests that climatic variations impacted similar biological pathways but different genes.

Search for Selection Signatures Related to Trypanosomosis Tolerance in African Goats.

Livestock is heavily affected by trypanosomosis in Africa. Through strong selective pressure, several African indigenous breeds of cattle and small ruminants have acquired varying degrees of tolerance against this disease. In this study, we combined LFMM and PCAdapt for analyzing two datasets of goats from West-Central Africa and East Africa, respectively, both comprising breeds with different assumed levels of trypanotolerance. The objectives were (i) to identify molecular signatures of selection related to trypanotolerance; and (ii) to guide an optimal sampling for subsequent studies. From 33 identified signatures, 18 had been detected previously in the literature as being mainly associated with climatic adaptations. The most plausible signatures of trypanotolerance indicate the genes DIS3L2, COPS7B, PD5A, UBE2K, and UBR1. The last gene is of particular interest since previous literature has already identified E3-ubiquitin ligases as playing a decisive role in the immune response. For following-up on these findings, the West-Central African area appears particularly relevant because of (i) a clear parasitic load gradient related to a humidity gradient, and (ii) still restricted admixture levels between goat breeds. This study illustrates the importance of protecting local breeds, which have retained unique allelic combinations conferring their remarkable adaptations.

Genome-Wide Detection of Structural Variations Reveals New Regions Associated with Domestication in Small Ruminants.

During domestication processes, changes in selective pressures induce multiple phenotypical, physiological, and behavioral changes in target species. The rise of next-generation sequencing has provided a chance to study the genetics bases of these changes, most of the time based on single nucleotide polymorphisms (SNPs). However, several studies have highlighted the impact of structural variations (SVs) on individual fitness, particularly in domestic species. We aimed at unraveling the role of SVs during the domestication and later improvement of small ruminants by analyzing whole-genome sequences of 40 domestic sheep and 11 of their close wild relatives (Ovis orientalis), and 40 goats and 18 of their close wild relatives (Capra aegagrus). Using a combination of detection tools, we called 45,796 SVs in Ovis and 15,047 SVs in Capra genomes, including insertions, deletions, inversions, copy number variations, and chromosomal translocations. Most of these SVs were previously unreported in small ruminants. 69 and 45 SVs in sheep and goats, respectively, were in genomic regions with neighboring SNPs highly differentiated between wilds and domestics (i.e., putatively related to domestication). Among them, 25 and 20 SVs were close to or overlapping with genes related to physiological and morpho-anatomical traits linked with productivity (e.g., size, meat or milk quality, wool color), reproduction, or immunity. Finally, several of the SVs differentiated between wilds and domestics would not have been detected by screening only the differentiation of SNPs surrounding them, highlighting the complementarity of SVs and SNPs based approaches to detect signatures of selection.

Altitudinal Zonation of Green Algae Biodiversity in the French Alps.

Mountain environments are marked by an altitudinal zonation of habitat types. They are home to a multitude of terrestrial green algae, who have to cope with abiotic conditions specific to high elevation, e.g., high UV irradiance, alternating desiccation, rain and snow precipitations, extreme diurnal variations in temperature and chronic scarceness of nutrients. Even though photosynthetic green algae are primary producers colonizing open areas and potential markers of climate change, their overall biodiversity in the Alps has been poorly studied so far, in particular in soil, where algae have been shown to be key components of microbial communities. Here, we investigated whether the spatial distribution of green algae followed the altitudinal zonation of the Alps, based on the assumption that algae settle in their preferred habitats under the pressure of parameters correlated with elevation. We did so by focusing on selected representative elevational gradients at distant locations in the French Alps, where soil samples were collected at different depths. Soil was considered as either a potential natural habitat or temporary reservoir of algae. We showed that algal DNA represented a relatively low proportion of the overall eukaryotic diversity as measured by a universal Eukaryote marker. We designed two novel green algae metabarcoding markers to amplify the Chlorophyta phylum and its Chlorophyceae class, respectively. Using our newly developed markers, we showed that elevation was a strong correlate of species and genus level distribution. Altitudinal zonation was thus determined for about fifty species, with proposed accessions in reference databases. In particular, Planophila laetevirens and Bracteococcus ruber related species as well as the snow alga Sanguina genus were only found in soil starting at 2,000 m above sea level. Analysis of environmental and bioclimatic factors highlighted the importance of pH and nitrogen/carbon ratios in the vertical distribution in soil. Capacity to grow heterotrophically may determine the Trebouxiophyceae over Chlorophyceae ratio. The intensity of freezing events (freezing degree days), proved also determinant in Chlorophyceae distribution. Guidelines are discussed for future, more robust and precise analyses of environmental algal DNA in mountain ecosystems and address green algae species distribution and dynamics in response to environmental changes.

An evaluation of sequencing coverage and genotyping strategies to assess neutral and adaptive diversity.

Whole genome sequences (WGS) greatly increase our ability to precisely infer population genetic parameters, demographic processes, and selection signatures. However, WGS may still be not affordable for a representative number of individuals/populations. In this context, our goal was to assess the efficiency of several SNP genotyping strategies by testing their ability to accurately estimate parameters describing neutral diversity and to detect signatures of selection. We analysed 110 WGS at 12× coverage for four different species, i.e., sheep, goats and their wild counterparts. From these data we generated 946 data sets corresponding to random panels of 1K to 5M variants, commercial SNP chips and exome capture, for sample sizes of five to 48 individuals. We also extracted low-coverage genome resequencing of 1×, 2× and 5× by randomly subsampling reads from the 12× resequencing data. Globally, 5K to 10K random variants were enough for an accurate estimation of genome diversity. Conversely, commercial panels and exome capture displayed strong ascertainment biases. Besides the characterization of neutral diversity, the detection of the signature of selection and the accurate estimation of linkage disequilibrium (LD) required high-density panels of at least 1M variants. Finally, genotype likelihoods increased the quality of variant calling from low coverage resequencing but proportions of incorrect genotypes remained substantial, especially for heterozygote sites. Whole genome resequencing coverage of at least 5× appeared to be necessary for accurate assessment of genomic variations. These results have implications for studies seeking to deploy low-density SNP collections or genome scans across genetically diverse populations/species showing similar genetic characteristics and patterns of LD decay for a wide variety of purposes.

Genetic homogenization of indigenous sheep breeds in Northwest Africa.

Northwest-African sheep represent an ideal case-study for assessing the potential impact of genetic homogenization as a threat to the future of traditional breeds that are adapted to local conditions. We studied ten Algerian and Moroccan breeds of sheep, including three transboundary breeds, distributed over a large part of the Maghreb region, which represents a geographically and historically coherent unit. Our analysis of the dataset that involved carrying out Genome-wide SNP genotyping, revealed a high level of homogenization (ADMIXTURE, NetView, fineSTRUCTURE and IBD segments analyses), in such a way that some breeds from different origins appeared genetically undistinguished: by grouping the eight most admixed populations, we obtained a mean global FST value of 0.0024. The sPCA analysis revealed that the major part of Morocco and the Northern part of Algeria were affected by the phenomenon, including most of the breeds considered. Unsupervised cross-breeding with the popular Ouled-Djellal breed was identified as a proximate cause of this homogenization. The issue of transboundary breeds was investigated, and the Hamra breed in particular was examined via ROH fragments analysis. Genetic diversity was considered in the light of historical archives and anthropological works. All of these elements taken together suggest that homogenization as a factor affecting the Maghrebin sheep stock, has been particularly significant over the last few decades, although this process probably started much earlier. In particular, we have identified the policies set by the French administration during the colonial period of the region's history as a causal factor that probably contributed significantly to this process. The genetic homogenization that we have observed calls into question the integrity of the farm animal genomic resources represented by these local breeds, whose conservation is of critical importance to the future of the livestock sector.

Old origin of a protective endogenous retrovirus (enJSRV) in the Ovis genus.

Sheep, the Jaagsiekte sheep retrovirus (JSRV) and its endogenous forms (enJSRVs) are a good model to study long-time relationships between retroviruses and their hosts. Taking advantage of 76 whole genome resequencing data of wild and domestic Ovis, we investigated the evolution of this relationship. An innovative analysis of re-sequencing data allowed characterizing 462 enJSRVs insertion sites (including 435 newly described insertions) in the Ovis genus. We focused our study on endogenous copies inserted in the q13 locus of chromosome 6 (6q13). Those copies are known to confer resistance against exogenous JSRV thanks to alleles bearing a mutation in the gag gene. We characterized (i) the distribution of protective and non-protective alleles across Ovis species and (ii) the copy number variation of the 6q13 locus. Our results challenged the previous hypothesis of fixation and amplification of the protective copies in relation with domestication, and allowed building a new model for the evolution of the 6q13 locus. JSRV would have integrated the 6q13 locus after the Ovis-Capra divergence (5-11 MYA) and before the Ovis diversification (2.4-5 MYA). The protective mutation in the enJSRV 6q13 copy appeared shortly after its insertion and was followed by genomic amplifications, after the divergence between Pachyform lineage on one side and the Argaliform and moufloniform lineages on the other (2.4-5 MYA). Considering the potential selective advantage of the protective mutation, its fixation in both sheep and its closest wild relative Ovis orientalis may be due to natural selection before domestication from O. orientalis populations.

Genetic homogeneity of North-African goats.

North Africa represents a rich and early reservoir of goat genetic diversity, from which the main African breeds have been derived. In this study, the genetic diversity of four indigenous Algerian goat breeds (i.e., Arabia, Makatia, M'Zabite and Kabyle, with n = 12 for each breed) has been investigated for the first time by genome-wide SNP genotyping; moreover in a broader context, genetic structuration of Algerian and Moroccan goats was explored (via FST, MDS, STRUCTURE, FineSTRUCTURE, BAPS, sPCA and DAPC analyses). At national level, the study revealed high level of genetic diversity and a significant phenomenon of admixture affecting all the Algerian breeds. At broader scale, clear global genetic homogeneity appeared considering both Algerian and Moroccan stocks. Indeed, genetic structuration was almost nonexistent among Arabia (from Algeria), Draa, Black and Nord (from Morocco), while the ancestral Kabyle and M'Zabite breeds, reared by Berber peoples, showed genetic distinctness. The study highlighted the threat to the Maghrebin stock, probably induced by unsupervised cross-breeding practices which have intensified in recent centuries. Moreover, it underlined the necessity to deepen our understanding of the genetic resources represented by the resilient North-African goat stock.

Ancient goat genomes reveal mosaic domestication in the Fertile Crescent.

Current genetic data are equivocal as to whether goat domestication occurred multiple times or was a singular process. We generated genomic data from 83 ancient goats (51 with genome-wide coverage) from Paleolithic to Medieval contexts throughout the Near East. Our findings demonstrate that multiple divergent ancient wild goat sources were domesticated in a dispersed process that resulted in genetically and geographically distinct Neolithic goat populations, echoing contemporaneous human divergence across the region. These early goat populations contributed differently to modern goats in Asia, Africa, and Europe. We also detect early selection for pigmentation, stature, reproduction, milking, and response to dietary change, providing 8000-year-old evidence for human agency in molding genome variation within a partner species.

Convergent genomic signatures of domestication in sheep and goats.

The evolutionary basis of domestication has been a longstanding question and its genetic architecture is becoming more tractable as more domestic species become genome-enabled. Before becoming established worldwide, sheep and goats were domesticated in the fertile crescent 10,500 years before present (YBP) where their wild relatives remain. Here we sequence the genomes of wild Asiatic mouflon and Bezoar ibex in the sheep and goat domestication center and compare their genomes with that of domestics from local, traditional, and improved breeds. Among the genomic regions carrying selective sweeps differentiating domestic breeds from wild populations, which are associated among others to genes involved in nervous system, immunity and productivity traits, 20 are common to Capra and Ovis. The patterns of selection vary between species, suggesting that while common targets of selection related to domestication and improvement exist, different solutions have arisen to achieve similar phenotypic end-points within these closely related livestock species.

Sheep genome functional annotation reveals proximal regulatory elements contributed to the evolution of modern breeds.

Domestication fundamentally reshaped animal morphology, physiology and behaviour, offering the opportunity to investigate the molecular processes driving evolutionary change. Here we assess sheep domestication and artificial selection by comparing genome sequence from 43 modern breeds (Ovis aries) and their Asian mouflon ancestor (O. orientalis) to identify selection sweeps. Next, we provide a comparative functional annotation of the sheep genome, validated using experimental ChIP-Seq of sheep tissue. Using these annotations, we evaluate the impact of selection and domestication on regulatory sequences and find that sweeps are significantly enriched for protein coding genes, proximal regulatory elements of genes and genome features associated with active transcription. Finally, we find individual sites displaying strong allele frequency divergence are enriched for the same regulatory features. Our data demonstrate that remodelling of gene expression is likely to have been one of the evolutionary forces that drove phenotypic diversification of this common livestock species.

Whole mitochondrial genomes unveil the impact of domestication on goat matrilineal variability.

BACKGROUND: The current extensive use of the domestic goat (Capra hircus) is the result of its medium size and high adaptability as multiple breeds. The extent to which its genetic variability was influenced by early domestication practices is largely unknown. A common standard by which to analyze maternally-inherited variability of livestock species is through complete sequencing of the entire mitogenome (mitochondrial DNA, mtDNA). RESULTS: We present the first extensive survey of goat mitogenomic variability based on 84 complete sequences selected from an initial collection of 758 samples that represent 60 different breeds of C. hircus, as well as its wild sister species, bezoar (Capra aegagrus) from Iran. Our phylogenetic analyses dated the most recent common ancestor of C. hircus to ~460,000 years (ka) ago and identified five distinctive domestic haplogroups (A, B1, C1a, D1 and G). More than 90 % of goats examined were in haplogroup A. These domestic lineages are predominantly nested within C. aegagrus branches, diverged concomitantly at the interface between the Epipaleolithic and early Neolithic periods, and underwent a dramatic expansion starting from ~12-10 ka ago. CONCLUSIONS: Domestic goat mitogenomes descended from a small number of founding haplotypes that underwent domestication after surviving the last glacial maximum in the Near Eastern refuges. All modern haplotypes A probably descended from a single (or at most a few closely related) female C. aegagrus. Zooarchaelogical data indicate that domestication first occurred in Southeastern Anatolia. Goats accompanying the first Neolithic migration waves into the Mediterranean were already characterized by two ancestral A and C variants. The ancient separation of the C branch (~130 ka ago) suggests a genetically distinct population that could have been involved in a second event of domestication. The novel diagnostic mutational motifs defined here, which distinguish wild and domestic haplogroups, could be used to understand phylogenetic relationships among modern breeds and ancient remains and to evaluate whether selection differentially affected mitochondrial genome variants during the development of economically important breeds.