Whole genome structural analysis of Caribbean hair sheep reveals quantitative link to West African ancestry

Mitochondrial DNA control-region and coding-region data highlight geographically structured diversity and post-domestication population dynamics in worldwide donkeys

Donkeys (Equus asinus) have been used extensively in agriculture and transportations since their domestication, ca. 5000-7000 years ago, but the increased mechanization of the last century has largely spoiled their role as burden animals, particularly in developed countries. Consequently, donkey breeds and population sizes have been declining for decades, and the diversity contributed by autochthonous gene pools has been eroded. Here, we examined coding-region data extracted from 164 complete mitogenomes and 1392 donkey mitochondrial DNA (mtDNA) control-region sequences to (i) assess worldwide diversity, (ii) evaluate geographical patterns of variation, and (iii) provide a new nomenclature of mtDNA haplogroups. The topology of the Maximum Parsimony tree confirmed the two previously identified major clades, i.e. Clades 1 and 2, but also highlighted the occurrence of a deep-diverging lineage within Clade 2 that left a marginal trace in modern donkeys. Thanks to the identification of stable and highly diagnostic coding-region mutational motifs, the two lineages were renamed as haplogroup A and haplogroup B, respectively, to harmonize clade nomenclature with the standard currently adopted for other livestock species. Control-region diversity and population expansion metrics varied considerably between geographical areas but confirmed North-eastern Africa as the likely domestication center. The patterns of geographical distribution of variation analyzed through phylogenetic networks and AMOVA confirmed the co-occurrence of both haplogroups in all sampled populations, while differences at the regional level point to the joint effects of demography, past human migrations and trade following the spread of donkeys out of the domestication center. Despite the strong decline that donkey populations have undergone for decades in many areas of the world, the sizeable mtDNA variability we scored, and the possible identification of a new early radiating lineage further stress the need for an extensive and large-scale characterization of donkey nuclear genome diversity to identify hotspots of variation and aid the conservation of local breeds worldwide.

Refining the genomic profiles of North African sheep breeds through meta-analysis of worldwide genomic SNP data

Introduction

The development of reproducible tools for the rapid genotyping of thousands of genetic markers (SNPs) has promoted cross border collaboration in the study of sheep genetic diversity on a global scale.

Methods

In this study, we collected a comprehensive dataset of 239 African and Eurasian sheep breeds genotyped at 37,638 filtered SNP markers, with the aim of understanding the genetic structure of 22 North African (NA) sheep breeds within a global context.

Results and discussion

We revealed asubstantial enrichment of the gene pool between the north and south shores of the Mediterranean Sea, which corroborates the importance of the maritime route in the history of livestock. The genetic structure of North African breeds mirrors the differential composition of genetic backgrounds following the breed history. Indeed, Maghrebin sheep stocks constitute a geographically and historically coherent unit with any breed-level genetic distinctness among them due to considerable gene flow. We detected a broad east-west pattern describing the most important trend in NA fat-tailed populations, exhibited by the genetic closeness of Egyptian and Libyan fat-tailed sheep to Middle Eastern breeds rather than Maghrebin ones. A Bayesian FST scan analysis revealed a set of genes with potentially key adaptive roles in lipid metabolism (BMP2, PDGFD VEGFA, TBX15, and WARS2), coat pigmentation (SOX10, PICK1, PDGFRA, MC1R, and MTIF) and horn morphology RXFP2) in Tunisian sheep. The local ancestry method detected a Merino signature in Tunisian Noire de Thibar sheep near the SULF1gene introgressed by Merino's European breeds. This study will contribute to the general picture of worldwide sheep genetic diversity.

Genome-wide association with footrot in hair and wool sheep

Ovine footrot is an infectious disease with important contributions from Dichelobacter nodosus and Fusobacterium necrophorum. Footrot is characterized by separation of the hoof from underlying tissue, and this causes severe lameness that negatively impacts animal wellbeing, growth, and profitability. Large economic losses result from lost production as well as treatment costs, and improved genetic tools to address footrot are a valuable long-term goal. Prior genetic studies had examined European wool sheep, but hair sheep breeds such as Katahdin and Blackbelly have been reported to have increased resistance to footrot, as well as to intestinal parasites. Thus, footrot condition scores were collected from 251 U.S. sheep including Katahdin, Blackbelly, and European-influenced crossbred sheep with direct and imputed genotypes at OvineHD array (>500,000 single nucleotide polymorphism) density. Genome-wide association was performed with a mixed model accounting for farm and principal components derived from animal genotypes, as well as a random term for the genomic relationship matrix. We identified three genome-wide significant associations, including SNPs in or near GBP6 and TCHH. We also identified 33 additional associated SNPs with genome-wide suggestive evidence, including a cluster of 6 SNPs in a peak near the genome-wide significance threshold located near the glutamine transporter gene SLC38A1. These findings suggest genetic susceptibility to footrot may be influenced by genes involved in divergent biological processes such as immune responses, nutrient availability, and hoof growth and integrity. This is the first genome-wide study to investigate susceptibility to footrot by including hair sheep and also the first study of any kind to identify multiple genome-wide significant associations with ovine footrot. These results provide a foundation for developing genetic tests for marker-assisted selection to improve resistance to ovine footrot once additional steps like fine mapping and validation are complete.

Recent advances in the genomic resources for sheep

Sheep (Ovis aries) provide a vital source of protein and fibre to human populations. In coming decades, as the pressures associated with rapidly changing climates increase, breeding sheep sustainably as well as producing enough protein to feed a growing human population will pose a considerable challenge for sheep production across the globe. High quality reference genomes and other genomic resources can help to meet these challenges by: (1) informing breeding programmes by adding a priori information about the genome, (2) providing tools such as pangenomes for characterising and conserving global genetic diversity, and (3) improving our understanding of fundamental biology using the power of genomic information to link cell, tissue and whole animal scale knowledge. In this review we describe recent advances in the genomic resources available for sheep, discuss how these might help to meet future challenges for sheep production, and provide some insight into what the future might hold.

Differentiation of Hair Sheep Breeds Based on the Physiological and Blood Biochemical Changes in Response to Different Stressors Using Multivariate Analysis Techniques

Physiological and blood measurement changes due to high heat load, restricted feed intake, and limited drinking water availability in 135 animals of three hair sheep breeds (Dorper, Katahdin, and St. Croix) were subjected to multivariate analysis techniques. The objective of this analysis was to evaluate the ability of these variables to separate individual hair sheep into groups based on adaptation characteristics in response to three physiological stressors and identify variables with greater discriminatory power. There were 16, 8, and 13 physiological and blood variables obtained from high heat load, restricted feed consumption, and water intake studies, respectively, for multivariate analysis. Physiological variables such as respiration rate, rectal and skin temperature, and panting score were measured only in the heat stress study. The results of the cluster and canonical discriminant analyses showed the presence of wide divergence (p < 0.05) between St. Croix and other breeds in their responses to high heat loads and restricted-feed- and -water-intake conditions. Dorper and Katahdin were grouped (p > 0.05) together based on the changes in physiological variables, which were separated (p < 0.05) from those of St. Croix as a resilient group. The stepwise discriminant analysis indicated that skin temperature, panting score, rectal temperature, respiration rate, and blood urea nitrogen and oxygen concentrations were the significant (p < 0.05) discriminating variables in clustering individual sheep into groups based on their responses to the high-heat-stress condition. Under the limited feed intake condition, the significant (p < 0.05) traits responsible for the separation of St. Croix from Dorper and Katahdin were blood triglyceride and cholesterol concentrations, whereas blood hemoglobin, osmolality, protein, and albumin were most important discriminating variables under the limited water intake condition. In conclusion, the results of the present study suggest that the stress responses of Dorper and Katahdin are similar and different from that of St. Croix. This finding can be useful information for future decisions in developing climate-resilient sheep through selective breeding.

High-density SNP markers elucidate the genetic divergence and population structure of Noticiana sheep breed in the Mediterranean context

Among livestock species, sheep have played an early major role in the Mediterranean area. Italy has a long history of sheep breeding and, despite a dramatic contraction in numbers, still raise several local populations that may represent a unique source of genetic diversity. The Noticiana is a breed of the south-eastern part of Sicily appreciated both for its dairy products and for its resistance to harsh environment. In this study, the high-density Illumina Ovine SNP600K BeadChip array was used for the first genome-wide characterization of 48 individuals of Noticiana sheep to investigate its diversity, the genome structure and the relationship within the context of worldwide and Italian breeds. Moreover, the runs of homozygosity (ROH) pattern and the pairwise F_ST-outliers were examined. Noticiana reported moderate levels of genetic diversity. The high percentage of short and medium length ROH segments (93% under 4 Mb) is indicative of a within breed relatedness dating back to ancient times, despite the absence of management for the mating plans and the reduced population size. In the worldwide context, the Southern Italian, Spanish and Albanian breeds overlapped in a macro cluster which also included the Noticiana sheep. The results highlighted ancestral genetic components of Noticiana shared with Comisana breed, and showed the clear separation from the other Italian sheep. This is likely the consequence of the combined effects of genetic drift, small population size and reproductive isolation. ROH islands and F_ST-outliers approaches in Noticiana identified genes and QTLs involved in milk and meat production, as well as related to the local adaptation, and therefore are consistent with the phenotypic traits of the studied breed. Although a wider sampling could be useful to deepen the genomic survey on Noticiana, these results represent a crucial starting point for the characterization of an important local genetic resource, with a view of supporting the local economy and preserving the biodiversity of the sheep species.

The ancestral origin of the critically endangered Quadricorna sheep as revealed by genome-wide analysis

Livestock European diffusion followed different human migration waves from the Fertile Crescent. In sheep, at least two diffusion waves have shaped the current breeds' biodiversity generating a complex genetic pattern composed by either primitive or fine-wool selected breeds. Among primitive breeds, aside from sharing common ancestral genomic components, they also show several traits such as the policeraty, large horns in the ram, short tail, and a moulting fleece, considered as ancestral. Although most of the primitive breeds characterized by these traits are confined on the very edge of Northern Europe, several residual populations are also scattered in the Mediterranean region. In fact, although in Italy a large number of local breeds are already extinct, others are listed as critically endangered, and among these there is the Quadricorna breed which is a four-horned sheep characterized by several ancestral traits. In this context we genotyped 47 individuals belonging to the Quadricorna sheep breed, a relict and endangered breed, from Central and Southern Italy. In doing so we used the Illumina OvineSNP50K array in order to explore its genetic diversity and to compare it with other 33 primitive traits-related, Mediterranean and Middle-East breeds, with the specific aim to reconstruct its origin. After retaining 35,680 SNPs following data filtering, the overall genomic architecture has been explored by using genetic diversity indices, Principal Component Analysis (PCA) and admixture analysis, while the genetic relationships and migration events have been inferred using a neighbor-joining tree based on Reynolds' distances and by the maximum likelihood tree as implemented in treemix. Multiple convergent evidence from all our population genetics analyses, indicated that the two Quadricorna populations differ from all the other Italian breeds, while they resulted to be very close to the Middle Eastern and primitive European breeds. In addition, the genetic diversity indices highlighted values comparable with those of most of the other analyzed breeds, despite the two populations exhibit slightly different genetic indices suggesting different levels of genomic inbreeding and drift (FIS and FROH). The admixture analysis does not suggest any signal of recent gene exchange with other Italian local breeds, highlighting a rather ancestral purity of the two populations, while on the other hand the treemix analysis seems to suggest an ancient admixture with other primitive European breeds. Finally, all these evidences seem to trace back the residual Quadricorna sheep to an early Neolithic spread, probably following a Mediterranean route and that urgent conservation actions are needed in order to keep the breed and all related cultural products alive.

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.

New Insight into the Genome-Wide Diversity and Admixture of Six Colombian Sheep Populations

Creole sheep represent a strategic genetic resource for populations living in marginal areas under financial restrictions on the American continent. Six Colombian sheep breeds (two wool (BCL-Boyacá and NCL-Nariño, 12 and 14 samples) and four hair (OPCE-Ethiopian, 54 samples; OPCS-Sudan, 74 samples; OPCP-Pelibeuy, 59 samples; OPCW-Wayúu, 24 samples) were genotyped using the Illumina Ovine SNP50 BeadChip. Data was also included from international 44 breeds from International Sheep Genomics Consortium (ISGC) and from data published in previous a previous work on the Caribbean and African breeds. Although geographically separated, wool (NCL, BCL) and hair types (OPCE, OPCS, OPCW) presented little genetic differentiation (FST 0.05) at a global level but several groups of animals separated suggesting local clustering due to geographical isolation. The OPCP underwent a recent crossing with Mexican Pelibuey, explaining its differentiation. Findings in this work such as the proximity to West African Djallonké (WAD) and Barbados Black Belly (BBB), suggest different introductions of African type animals from the Caribbean region on a pre-existing genetic basis formed by animals deriving from the first importations coming from Europe in colonial times. As expected, Colombian wool breeds showed, in particular in Admixture software results, a greater genomic component in common with European breeds and in particular with Iberian ones (Churra). This study provides a basis for future research into the genetic diversity within and between the Colombian sheep breeds analysed, and scientific data for policy decisions on Farm Animal Genetic Resources (FAnGR).

Analysis of the Genetic Diversity and Population Structure of Four Senegalese Sheep Breeds Using Medium-Density Single-Nucleotide Polymorphisms

Simple Summary

This paper reported genetic parameters of four Senegalese sheep breeds, in relation to inbreeding, diversity and genetic proximity. The results provide informations on genetic conservation and adaptability of the breeds in the Senegalese context.

Abstract

In Senegal, sheep breeds have adapted to their environment and play a key socio-economic role. This study aimed to explore the genetic diversity and structure of four Senegalese sheep breeds (Peul-peul, Djallonke, Touabire, and Ladoum) and their relationships with global sheep breeds. To that end, forty-seven sheep were genotyped using the OvineSNP50 BeadChip, and these genotypic data were analysed with those of 73 sheep breeds representative of worldwide ovine diversity (2729 animals). The average observed heterozygosity (Ho) ranged from 0.293 in Djallonke sheep to 0.339 in Touabire sheep. The estimated Fis values were low, ranging from 0.019 for Ladoum to 0.034 for Peul-peul sheep. The estimated Fst values were low (0.003–0.044) among the trypanosusceptible breeds (Peul-peul, Touabire, and Ladoum) but high between the previous breeds and the trypanotolerant Djallonke breed (0.075–0.116), indicating better genetic conservation of the Djallonke sheep. A principal component analysis revealed clustering of the Senegalese sheep breeds according to their geographic distribution. However, owing to genetic improvement practices, the introgression of Touabire sheep blood seems to have reshaped the genetic landscape of the trypanosusceptible sheep breeds in Senegal. The Senegalese sheep breeds showed lower genetic diversity than their presumed ancestral sheep breeds of the Middle East. They also presented some relatedness with Caribbean sheep breeds, which reveals their contribution to the global genetic diversity and to the development of Caribbean sheep breeds.

Sheep Post-Domestication Expansion in the Context of Mitochondrial and Y Chromosome Haplogroups and Haplotypes

Mitochondrial DNA and nonrecombinant parts of Y-chromosome DNA are a great tool for looking at a species’ past. They are inherited for generations almost unaffected because they do not participate in recombination; thus, the time of occurrence of each mutation can be estimated based on the average mutation rate. Thanks to this, male and female haplogroups guide confirming events in the distant past (potential centers of domestication, settlement of areas, trade connections) as well as in modern breeding (crossbreeding, confirmation of paternity). This research focuses mainly on the development of domestic sheep and its post-domestication expansion, which has occurred through human trade from one continent to another. So far, five mitochondrial and five Y-chromosome haplogroups and dozens of their haplotypes have been detected in domestic sheep through studies worldwide. Mitochondrial DNA variability is more or less correlated with distance from the domestication center, but variability on the recombinant region of the Y chromosome is not. According to available data, central China shows the highest variability of male haplogroups and haplotypes.

Variability in Global DNA Methylation Rate Across Tissues and Over Time in Sheep

Recent studies highlighted the influence of epigenetic marks in the variability of many complex traits, both in plants and animals. These studied focused only on specific sites of the genome having differentially methylated profiles among individuals and/or tissues. In contrast, we recently used the methylation rate of the entire genome as a unique measure considered as a novel quantitative phenotype in sheep. This phenotype named global DNA methylation rate (GDMR), measured by luminometric assay, integrates the methylation level of each CpG dinucleotide within the 6 million of CCGG sites along the ovine genome. GDMR measured in blood previously showed moderate heritability of 0.20 and provided evidence for a genetic determinism. The main objective of the present study was to better characterize the GDMR phenotype in various tissues and investigate its variability in several breeds of sheep reared in the same environment. GDMR was measured on blood samples collected monthly from 59 growing male and female lambs (24 Romane, 23 Blackbelly and 12 Charollais), between birth and 4 months of age. Blood GDMR was on average around 80% and was influenced by the sampling date (p &lt; 0.001), the breed (p = 0.002) and the sex (p = 0.002). In addition, GDMR was determined in 12 somatic (frontal lobe, pituitary gland, heart, lung, sub cutaneous and perirenal adipose tissue, skeletal muscle, liver, spleen, adrenal gland, medulla and cortical kidney) and 6 reproductive tissues (ovary, oviduct, uterus, testis, epididymis and seminal vesicle). GDMR was on average 70% in somatic tissues but marked variation was observed depending on the tissue. The GDMR measured in blood was higher than that measured in other somatic tissues, and is not a good proxy of less accessible tissues. Female reproductive tissues had a 10% higher GDMR than male reproductive tissues. We demonstrated a significant influence of the breed on blood GDMR, certainly reflecting the influence of different genetic backgrounds. The effect of the breed on GDMR may be related to their specific abilities to adapt to and live in different conditions.

Whole-genome resequencing of worldwide wild and domestic sheep elucidates genetic diversity, introgression and agronomically important loci

Domestic sheep and their wild relatives harbor substantial genetic variants that can form the backbone of molecular breeding, but their genome landscapes remain understudied. Here, we present a comprehensive genome resource for wild ovine species, landraces and improved breeds of domestic sheep, comprising high-coverage (∼16.10×) whole genomes of 810 samples from 7 wild species and 158 diverse domestic populations. We detected, in total, ∼121.2 million single nucleotide polymorphisms, ∼61 million of which are novel. Some display significant (P < 0.001) differences in frequency between wild and domestic species, or are private to continent-wide or individual sheep populations. Retained or introgressed wild gene variants in domestic populations have contributed to local adaptation, such as the variation in the HBB associated with plateau adaptation. We identified novel and previously reported targets of selection on morphological and agronomic traits such as stature, horn, tail configuration, and wool fineness. We explored the genetic basis of wool fineness and unveiled a novel mutation (chr25: T7,068,586C) in the 3′-UTR of IRF2BP2 as plausible causal variant for fleece fiber diameter. We reconstructed prehistorical migrations from the Near Eastern domestication center to South-and-Southeast Asia and found two main waves of migrations across the Eurasian Steppe and the Iranian Plateau in the Early and Late Bronze Ages. Our findings refine our understanding of genome variation as shaped by continental migrations, introgression, adaptation, and selection of sheep.

Origin and population structure of Brazilian hair sheep breeds

Brazilian hair sheep constitute a genetic diversity hotspot. These animals are found in the harsh environments of the Brazilian Northwest (semi-arid) region. Genotypes (50K SNP chip) from seven Brazilian sheep breeds (five hair and two coarse wool types) and 87 worldwide breeds were used to test for population structure, admixture and genetic diversity. Moreover, phylogenetic trees evaluating migration events between genetic groups were built. Brazilian Somali, a fat-tailed breed, had a close relationship with East African breeds and clustered distinctly from other Brazilian breeds. Brazilian Blackbelly and Barbados Blackbelly had a close relationship. The Morada Nova breed did not show close relationships with European or African breeds, revealing a single migration event from an Algerian hair breed. Brazilian Fat-tail and Morada Nova share a common ancestor, but the former showed introgressions from Brazilian Somali and Afrikaner breeds, explaining the fat-tail phenotype. The Santa Inês breed received a substantial contribution from Brazilian Bergamasca and showed an admixed origin with recent introgressions from other breeds, mainly from Suffolk. Furthermore, Brazilian Somali and Brazilian Fat-tail are the most endangered sheep genetic resources in Brazil and should be the focus for ex situ conservation programs. In conclusion, Brazilian hair sheep show an African origin and are characterized by diverse genetic composition, reinforcing the need for conservation of these genetic resources, and at the same time, this highly diverse group has variability that can be used in breeding programs.

Mitochondrial DNA Variations in Colombian Creole Sheep Confirm an Iberian Origin and Shed Light on the Dynamics of Introduction Events of African Genotypes

The genetic origins and diversity of Creole sheep from five regions of Colombia were investigated based on mitochondrial DNA (mtDNA) variations across 89 sequences from five breeds: one wool Creole sheep (CL) and four hair Creole sheep, including Ethiopian (OPCE), Sudan (OPCS), Pelibuey (OPCP) and Wayúu (OPCW). A global comparison was done using 62 haplotypes from Iberian, African, Indian, Caribbean, Mexican, Caucasian and European sheep based on sequences retrieved from GenBank. This study aimed to identify the maternal origin of Colombian Creole sheep and their genetic relationships at a global level. The results showed 31 different haplotypes from Colombian Creole sheep, which can be assigned to maternal lineage B, the most common lineage found in European sheep breeds and the only one found in several Iberian breed (e.g., Churra, Spanish Merino) that most likely participated in the Creole formation. Additional analyses showed that wool and hair sheep retained a broad genetic identity despite being geographically separated. The global-level phylogenetic analysis revealed that Colombian Creole sheep belong to a distinct and defined genetic lineage that is likely the result of a founder effect with ecotypes of Iberian descent and the subsequent introduction of foreign breeds. This is consistent with historical reports on the presence of sheep in South America and, particularly, Colombia.

Genetic diversity evolution of a sheep breed reintroduced after extinction: Tracing back Christopher Columbus' first imported sheep

New World's hair sheep breeds may genetically stem from West African introgression into established ecotypes of Spanish descent presumably extinct in the XIX Century. However, present Canary non-wooled breeds have presumably regressively resulted from the absorption of primitive individuals through Venezuelan descendants. We studied genetic diversity, structure, and evolution of the Canary hair sheep since its reintroduction in the 1950s. Demographic and genetic variability were evaluated using ENDOG (v4.8). Effective population size based on individual inbreeding rate was around one third higher than when based on individual coancestry rate. Nei's distances and equivalent subpopulations number indicated a highly-structured population. Although genetic diversity loss since the founder generations could be considered small, narrower pedigree bottlenecks could result from intraflock breeding policies and excessive contribution of few ancestors. Long generation intervals could be considered when reducing inbreeding. Wright's fixation statistics indicated slight interflock inbreeding. Pedigree completeness suggested genetic parameters were reliable, hence controlling inbreeding negative effects, could indeed, be crucial preserving these animal resources, consolidating the population in the archipelago after reintroduction.

On the origin of European sheep as revealed by the diversity of the Balkan breeds and by optimizing population-genetic analysis tools

Background

In the Neolithic, domestic sheep migrated into Europe and subsequently spread in westerly and northwesterly directions. Reconstruction of these migrations and subsequent genetic events requires a more detailed characterization of the current phylogeographic differentiation.

Results

We collected 50 K single nucleotide polymorphism (SNP) profiles of Balkan sheep that are currently found near the major Neolithic point of entry into Europe, and combined these data with published genotypes from southwest-Asian, Mediterranean, central-European and north-European sheep and from Asian and European mouflons. We detected clines, ancestral components and admixture by using variants of common analysis tools: geography-informative supervised principal component analysis (PCA), breed-specific admixture analysis, across-breed \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{4}$$\end{document}f4 profiles and phylogenetic analysis of regional pools of breeds. The regional Balkan sheep populations exhibit considerable genetic overlap, but are clearly distinct from the breeds in surrounding regions. The Asian mouflon did not influence the differentiation of the European domestic sheep and is only distantly related to present-day sheep, including those from Iran where the mouflons were sampled. We demonstrate the occurrence, from southeast to northwest Europe, of a continuously increasing ancestral component of up to 20% contributed by the European mouflon, which is assumed to descend from the original Neolithic domesticates. The overall patterns indicate that the Balkan region and Italy served as post-domestication migration hubs, from which wool sheep reached Spain and north Italy with subsequent migrations northwards. The documented dispersal of Tarentine wool sheep during the Roman period may have been part of this process. Our results also reproduce the documented 18th century admixture of Spanish Merino sheep into several central-European breeds.

Conclusions

Our results contribute to a better understanding of the events that have created the present diversity pattern, which is relevant for the management of the genetic resources represented by the European sheep population.

The global diversity of Haemonchus contortus is shaped by human intervention and climate

Haemonchus contortus is a haematophagous parasitic nematode of veterinary interest. We have performed a survey of its genome-wide diversity using single-worm whole genome sequencing of 223 individuals sampled from 19 isolates spanning five continents. We find an African origin for the species, together with evidence for parasites spreading during the transatlantic slave trade and colonisation of Australia. Strong selective sweeps surrounding the β-tubulin locus, a target of benzimidazole anthelmintic drug, are identified in independent populations. These sweeps are further supported by signals of diversifying selection enriched in genes involved in response to drugs and other anthelmintic-associated biological functions. We also identify some candidate genes that may play a role in ivermectin resistance. Finally, genetic signatures of climate-driven adaptation are described, revealing a gene acting as an epigenetic regulator and components of the dauer pathway. These results begin to define genetic adaptation to climate in a parasitic nematode.

Based on single worm whole genome sequencing, the authors here characterise the global evolution of the gastrointestinal parasite Haemonchus contortus and identify genes that play a role in drug resistance as well as climate-driven adaptations involving an epigenetic regulator.

Inferring the population structure of the Maghreb sheep breeds using a medium-density SNP chip

North Africa has a great diversity of indigenous sheep breeds whose origin is linked to its environmental characteristics and to certain historical events that took place in the region. To date, few genome-wide studies have been conducted to investigate the population structure of North African indigenous sheep. The objective of the present study was to provide a detailed assessment of the genetic structure and admixture patterns of six Maghreb sheep populations using the Illumina 50K Ovine BeadChip and comparisons with 22 global populations of sheep and mouflon. Regardless of the method of analysis used, patterns of multiple hybridization events were observed within all North African populations, leading to a heterogeneous genetic architecture that varies according to the breed. The Barbarine population showed the lowest genetic heterogeneity and major southwest Asian ancestry, providing additional support to the Asian origin of the North African fat-tailed sheep. All other breeds presented substantial Merino introgression ranging from 15% for D'man to 31% for Black Thibar. We highlighted several signals of ancestral introgression between North African and southern European sheep. In addition, we identified two opposite gradients of ancestry, southwest Asian and central European, occurring between North Africa and central Europe. Our results provide further evidence of the weak global population structure of sheep resulting from high levels of gene flow among breeds occurring worldwide. At the regional level, signs of recent admixture among North African populations, resulting in a change of the original genomic architecture of minority breeds, were also detected.

Genome-wide SNP profiling of worldwide goat populations reveals strong partitioning of diversity and highlights post-domestication migration routes

BACKGROUND

Goat populations that are characterized within the AdaptMap project cover a large part of the worldwide distribution of this species and provide the opportunity to assess their diversity at a global scale. We analysed genome-wide 50 K single nucleotide polymorphism (SNP) data from 144 populations to describe the global patterns of molecular variation, compare them to those observed in other livestock species, and identify the drivers that led to the current distribution of goats.

RESULTS

A high degree of genetic variability exists among the goat populations studied. Our results highlight a strong partitioning of molecular diversity between and within continents. Three major gene pools correspond to goats from Europe, Africa and West Asia. Dissection of sub-structures disclosed regional gene pools, which reflect the main post-domestication migration routes. We also identified several exchanges, mainly in African populations, and which often involve admixed and cosmopolitan breeds. Extensive gene flow has taken place within specific areas (e.g., south Europe, Morocco and Mali-Burkina Faso-Nigeria), whereas elsewhere isolation due to geographical barriers (e.g., seas or mountains) or human management has decreased local gene flows.

CONCLUSIONS

After domestication in the Fertile Crescent in the early Neolithic era (ca. 12,000 YBP), domestic goats that already carried differentiated gene pools spread to Europe, Africa and Asia. The spread of these populations determined the major genomic background of the continental populations, which currently have a more marked subdivision than that observed in other ruminant livestock species. Subsequently, further diversification occurred at the regional level due to geographical and reproductive isolation, which was accompanied by additional migrations and/or importations, the traces of which are still detectable today. The effects of breed formation were clearly detected, particularly in Central and North Europe. Overall, our results highlight a remarkable diversity that occurs at the global scale and is locally partitioned and often affected by introgression from cosmopolitan breeds. These findings support the importance of long-term preservation of goat diversity, and provide a useful framework for investigating adaptive introgression, directing genetic improvement and choosing breeding targets.

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.