Natural adaptation and human selection of northeast African sheep genomes.

African sheep manifest diverse but distinct physio-anatomical traits, which are the outcomes of natural- and human-driven selection. Here, we generated 34.8 million variants from 150 indigenous northeast African sheep genomes sequenced at an average depth of ∼54× for 130 samples (Ethiopia, Libya) and ∼20× for 20 samples (Sudan). These represented sheep from diverse environments, tail morphology and post-Neolithic introductions to Africa. Phylogenetic and model-based admixture analysis provided evidence of four genetic groups corresponding to altitudinal geographic origins, tail morphotypes and possible historical introduction and dispersal of the species into and across the continent. Running admixture at higher levels of K (6 ≤ K ≤ 25), revealed cryptic levels of genome intermixing as well as distinct genetic backgrounds in some populations. Comparative genomic analysis identified targets of selection that spanned conserved haplotype structures overlapping clusters of genes and gene families. These were related to hypoxia responses, ear morphology, caudal vertebrae and tail skeleton length, and tail fat-depot structures. Our findings provide novel insights underpinning morphological variation and response to human-driven selection and environmental adaptation in African indigenous sheep.

Integrated Environmental and Genomic Analysis Reveals the Drivers of Local Adaptation in African Indigenous Chickens.

Breeding for climate resilience is currently an important goal for sustainable livestock production. Local adaptations exhibited by indigenous livestock allow investigating the genetic control of this resilience. Ecological niche modeling (ENM) provides a powerful avenue to identify the main environmental drivers of selection. Here, we applied an integrative approach combining ENM with genome-wide selection signature analyses (XPEHH and Fst) and genotype-environment association (redundancy analysis), with the aim of identifying the genomic signatures of adaptation in African village chickens. By dissecting 34 agro-climatic variables from the ecosystems of 25 Ethiopian village chicken populations, ENM identified six key drivers of environmental challenges: One temperature variable-strongly correlated with elevation, three precipitation variables as proxies for water availability, and two soil/land cover variables as proxies of food availability for foraging chickens. Genome analyses based on whole-genome sequencing (n = 245), identified a few strongly supported genomic regions under selection for environmental challenges related to altitude, temperature, water scarcity, and food availability. These regions harbor several gene clusters including regulatory genes, suggesting a predominantly oligogenic control of environmental adaptation. Few candidate genes detected in relation to heat-stress, indicates likely epigenetic regulation of thermo-tolerance for a domestic species originating from a tropical Asian wild ancestor. These results provide possible explanations for the rapid past adaptation of chickens to diverse African agro-ecologies, while also representing new landmarks for sustainable breeding improvement for climate resilience. We show that the pre-identification of key environmental drivers, followed by genomic investigation, provides a powerful new approach for elucidating adaptation in domestic animals.

The maternal origin of indigenous domestic chicken from the Middle East, the north and the horn of Africa.

BACKGROUND: Indigenous domestic chicken represents a major source of protein for agricultural communities around the world. In the Middle East and Africa, they are adapted to hot dry and semi-dry areas, in contrast to their wild ancestor, the Red junglefowl, which lives in humid and sub-humid tropical areas. Indigenous populations are declining following increased demand for poultry meat and eggs, favouring the more productive exotic commercial breeds. In this paper, using the D-loop of mitochondrial DNA as a maternally inherited genetic marker, we address the question of the origin and dispersal routes of domestic chicken of the Middle East (Iraq and Saudi Arabia), the northern part of the African continent (Algeria and Libya) and the Horn of Africa (Ethiopia). RESULTS: The analysis of the mtDNA D-loop of 706 chicken samples from Iraq (n = 107), Saudi Arabia (n = 185), Algeria (n = 88), Libya (n = 23), Ethiopia (n = 211) and Pakistan (n = 92) show the presence of five haplogroups (A, B, C, D and E), suggesting more than one maternal origin for the studied populations. Haplogroup E, which occurred in 625 samples, was the most frequent in all countries. This haplogroup most likely originates from the Indian subcontinent and probably migrated following a terrestrial route to these different countries. Haplotypes belonging to haplogroup D were present in all countries except Algeria and Libya, it is likely a legacy of the Indian Ocean maritime trading network. Haplogroup A was present in all countries and may be of commercial origin. Haplogroup B was found only in Ethiopia. Haplogroup C was only detected in the South-Western region of Saudi Arabia and in Ethiopia. CONCLUSION: The results support a major influence of the Indian subcontinent on the maternal diversity of the today's chicken populations examined here. Most of the diversity occurs within rather than between populations. This lack of phylogeographic signal agrees with both ancient and more recent trading networks having shaped the modern-day diversity of indigenous chicken across populations and countries.

Diversity of endogenous avian leukosis virus subgroup E (ALVE) insertions in indigenous chickens.

BACKGROUND: Avian leukosis virus subgroup E (ALVE) insertions are endogenous retroviruses (ERV) that are restricted to the domestic chicken and its wild progenitor. In commercial chickens, ALVE are known to have a detrimental effect on productivity and provide a source for recombination with exogenous retroviruses. The wider diversity of ALVE in non-commercial chickens and the role of these elements in ERV-derived immunity (EDI) are yet to be investigated. RESULTS: In total, 974 different ALVE were identified from 407 chickens sampled from village populations in Ethiopia, Iraq, and Nigeria, using the recently developed obsERVer bioinformatics identification pipeline. Eighty-eight percent of all identified ALVE were novel, bringing the known number of ALVE integrations to more than 1300 across all analysed chickens. ALVE content was highly lineage-specific and populations generally exhibited a large diversity of ALVE at low frequencies, which is typical for ERV involved in EDI. A significantly larger number of ALVE was found within or near coding regions than expected by chance, although a relative depletion of ALVE was observed within coding regions, which likely reflects selection against deleterious integrations. These effects were less pronounced than in previous analyses of chickens from commercial lines. CONCLUSIONS: Identification of more than 850 novel ALVE has trebled the known diversity of these retroviral elements. This work provides the basis for future studies to fully quantify the role of ALVE in immunity against exogenous ALV, and development of programmes to improve the productivity and welfare of chickens in developing economies.

Comparative analysis of the chicken IFITM locus by targeted genome sequencing reveals evolution of the locus and positive selection in IFITM1 and IFITM3.

BACKGROUND: The interferon-induced transmembrane (IFITM) protein family comprises a class of restriction factors widely characterised in humans for their potent antiviral activity. Their biological activity is well documented in several animal species, but their genetic variation and biological mechanism is less well understood, particularly in avian species. RESULTS: Here we report the complete sequence of the domestic chicken Gallus gallus IFITM locus from a wide variety of chicken breeds to examine the detailed pattern of genetic variation of the locus on chromosome 5, including the flanking genes ATHL1 and B4GALNT4. We have generated chIFITM sequences from commercial breeds (supermarket-derived chicken breasts), indigenous chickens from Nigeria (Nsukka) and Ethiopia, European breeds and inbred chicken lines from the Pirbright Institute, totalling of 206 chickens. Through mapping of genetic variants to the latest chIFITM consensus sequence our data reveal that the chIFITM locus does not show structural variation in the locus across the populations analysed, despite spanning diverse breeds from different geographic locations. However, single nucleotide variants (SNVs) in functionally important regions of the proteins within certain groups of chickens were detected, in particular the European breeds and indigenous birds from Ethiopia and Nigeria. In addition, we also found that two out of four SNVs located in the chIFITM1 (Ser36 and Arg77) and chIFITM3 (Val103) proteins were simultaneously under positive selection. CONCLUSIONS: Together these data suggest that IFITM genetic variation may contribute to the capacities of different chicken populations to resist virus infection.

Altitude-dependent agro-ecologies impact the microbiome diversity of scavenging indigenous chicken in Ethiopia.

BACKGROUND: Scavenging indigenous village chickens play a vital role in sub-Saharan Africa, sustaining the livelihood of millions of farmers. These chickens are exposed to vastly different environments and feeds compared to commercial chickens. In this study, we analysed the caecal microbiota of 243 Ethiopian village chickens living in different altitude-dependent agro-ecologies. RESULTS: Differences in bacterial diversity were significantly correlated with differences in specific climate factors, topsoil characteristics, and supplemental diets provided by farmers. Microbiota clustered into three enterotypes, with one particularly enriched at high altitudes. We assembled 9977 taxonomically and functionally diverse metagenome-assembled genomes. The vast majority of these were not found in a dataset of previously published chicken microbes or in the Genome Taxonomy Database. CONCLUSIONS: The wide functional and taxonomic diversity of these microbes highlights their importance in the local adaptation of indigenous poultry, and the significant impacts of environmental factors on the microbiota argue for further discoveries in other agro-ecologies. Video Abstract.

Whole genome sequences of 234 indigenous African chickens from Ethiopia.

Indigenous chickens predominate poultry production in Africa. Although preferred for backyard farming because of their adaptability to harsh tropical environments, these populations suffer from relatively low productivity compared to commercial lines. Genome analyses can unravel the genetic potential of improvement of these birds for both production and resilience traits for the benefit of African poultry farming systems. Here we report whole-genome sequences of 234 indigenous chickens from 24 Ethiopian populations distributed under diverse agro-climatic conditions. The data represents over eight terabytes of paired-end sequences from the Ilumina HiSeqX platform with an average coverage of about 57X. Almost 99% of the sequence reads could be mapped against the chicken reference genome (GRCg6a), confirming the high quality of the data. Variant calling detected around 15 million SNPs, of which about 86% are known variants (i.e., present in public databases), providing further confidence on the data quality. The dataset provides an excellent resource for investigating genetic diversity and local environmental adaptations with important implications for breed improvement and conservation purposes.

Phenotypic and genetic parameter estimates for early growth, growth rate and growth efficiency-related traits of Fogera cattle in Ethiopia.

BACKGROUND: Understanding the phenotypic and genetic parameter estimates of growth traits is important for an effective livestock genetic improvement programme. OBJECTIVES: In this study, we evaluated the phenotypic performances and estimated genetic parameters for birthweight (BWT), weaning weight (WWT), pre-weaning average daily gain (PADG), pre-weaning Kleiber ratio (PKR), pre-weaning growth efficiency (PGE) and pre-weaning relative growth rate (PRGR) of Fogera cattle in Ethiopia. METHODS: Growth data collected from 2000 to 2018 in Andassa Livestock Research Center were used for the study. General linear model of SAS 9.1 was used to estimate the least squares mean (LSM) ± standard error (SE) for phenotypic performances, and AI-REML of Wombat software combined with a series of five single-trait animal models to estimate phenotypic variance and its direct, maternal and residual components. Calf sex, calf birth season and calf birth year were the fixed effects considered. RESULTS: The overall LSM ± SE BWT, WWT, PADG, PKR, PGE and PRGR were 21.28 ± 0.05 kg, 97.99 ± 0.67 kg, 320.29 ± 2.79 g, 10.10 ± 0.04, 3.51 ± 0.35 and 1.95 ± 0.00, respectively. All the fixed effects considered significantly (p < 0.001) affected all the traits. The direct heritability estimates for BWT, WWT, PADG, PKR, PGE and PRGR were 0.21 ± 0.07, 0.26 ± 0.01, 0.55 ± 0.19, 0.53 ± 0.18, 0.33 ± 0.00 and 0.50 ± 0.00, respectively. The genetic correlations among the traits ranged from negative (-0.20 ± 0.04; BWT-PKR) to positive (0.99 ± 0.00; BW-PGE, BW-GR, WWT-PGE, WWT-PGR, ADG-PGR, PKR-PGR, PKR-PGE and PGE-PGR). Similarly, the phenotypic correlations ranged from -0.03 ± 0.20 to 0.99 ± 0.01; BWT-PGE, BWT-PRGR, WWT-PGE, WWT-PRGR, PKR-PGE, PKR-PRGR and PGE-PRGR). CONCLUSION: The positive and larger phenotypic and genetic correlations between most of the traits implied that selection based on one trait could improve the other traits. However, the negative phenotypic and genetic correlation between BWT-PKRA implies that selection of Fogera calves based on either of the traits has an adverse effect on the other. Therefore, caution should be taken when designing the selection criteria for growth improvement.

Whole-Genome Sequence Data Suggest Environmental Adaptation of Ethiopian Sheep Populations.

Great progress has been made over recent years in the identification of selection signatures in the genomes of livestock species. This work has primarily been carried out in commercial breeds for which the dominant selection pressures are associated with artificial selection. As agriculture and food security are likely to be strongly affected by climate change, a better understanding of environment-imposed selection on agricultural species is warranted. Ethiopia is an ideal setting to investigate environmental adaptation in livestock due to its wide variation in geo-climatic characteristics and the extensive genetic and phenotypic variation of its livestock. Here, we identified over three million single nucleotide variants across 12 Ethiopian sheep populations and applied landscape genomics approaches to investigate the association between these variants and environmental variables. Our results suggest that environmental adaptation for precipitation-related variables is stronger than that related to altitude or temperature, consistent with large-scale meta-analyses of selection pressure across species. The set of genes showing association with environmental variables was enriched for genes highly expressed in human blood and nerve tissues. There was also evidence of enrichment for genes associated with high-altitude adaptation although no strong association was identified with hypoxia-inducible-factor (HIF) genes. One of the strongest altitude-related signals was for a collagen gene, consistent with previous studies of high-altitude adaptation. Several altitude-associated genes also showed evidence of adaptation with temperature, suggesting a relationship between responses to these environmental factors. These results provide a foundation to investigate further the effects of climatic variables on small ruminant populations.

Functional evolution of the colony-stimulating factor 1 receptor (CSF1R) and its ligands in birds.

Macrophage colony-stimulating factor (CSF1 or M-CSF) and interleukin 34 (IL34) are secreted cytokines that control macrophage survival and differentiation. Both act through the CSF1 receptor (CSF1R), a type III transmembrane receptor tyrosine kinase. The functions of CSF1R and both ligands are conserved in birds. We have analyzed protein-coding sequence divergence among avian species. The intracellular tyrosine kinase domain of CSF1R was highly conserved in bird species as in mammals but the extracellular domain of avian CSF1R was more divergent in birds with multiple positively selected amino acids. Based upon crystal structures of the mammalian CSF1/IL34 receptor-ligand interfaces and structure-based alignments, we identified amino acids involved in avian receptor-ligand interactions. The contact amino acids in both CSF1 and CSF1R diverged among avian species. Ligand-binding domain swaps between chicken and zebra finch CSF1 confirmed the function of variants that confer species specificity on the interaction of CSF1 with CSF1R. Based upon genomic sequence analysis, we identified prevalent amino acid changes in the extracellular domain of CSF1R even within the chicken species that distinguished commercial broilers and layers and tropically adapted breeds. The rapid evolution in the extracellular domain of avian CSF1R suggests that at least in birds this ligand-receptor interaction is subjected to pathogen selection. We discuss this finding in the context of expression of CSF1R in antigen-sampling and antigen-presenting cells.

Genome-wide insights of Ethiopian indigenous sheep populations reveal the population structure related to tail morphology and phylogeography.

BACKGROUND: Ethiopian sheep living in different climatic zones and having contrasting morphologies are a most promising subject of molecular-genetic research. Elucidating their genetic diversity and genetic structure is critical for designing appropriate breeding and conservation strategies. OBJECTIVE: The study was aimed to investigate genome-wide genetic diversity and population structure of eight Ethiopian sheep populations. METHODS: A total of 115 blood samples were collected from four Ethiopian sheep populations that include Washera, Farta and Wollo (short fat-tailed) and Horro (long fat-tailed). DNA was extracted using Quick-DNA™ Miniprep plus kit. All DNA samples were genotyped using Ovine 50 K SNP BeadChip. To infer genetic relationships of Ethiopian sheep at national, continental and global levels, genotype data on four Ethiopian sheep (Adilo, Arsi-Bale, Menz and Black Head Somali) and sheep from east, north, and south Africa, Middle East and Asia were included in the study as reference. RESULTS: Mean genetic diversity of Ethiopian sheep populations ranged from 0.352 ± 0.14 for Horro to 0.379 ± 0.14 for Arsi-Bale sheep. Population structure and principal component analyses of the eight Ethiopian indigenous sheep revealed four distinct genetic cluster groups according to their tail phenotype and geographical distribution. The short fat-tailed sheep did not represent one genetic cluster group. Ethiopian fat-rump sheep share a common genetic background with the Kenyan fat-tailed sheep. CONCLUSION: The results of the present study revealed the principal component and population structure follows a clear pattern of tail morphology and phylogeography. There is clear signature of admixture among the study Ethiopian sheep populations.

Genome-Wide Variation, Candidate Regions and Genes Associated With Fat Deposition and Tail Morphology in Ethiopian Indigenous Sheep.

Variations in body weight and in the distribution of body fat are associated with feed availability, thermoregulation, and energy reserve. Ethiopia is characterized by distinct agro-ecological and human ethnic farmer diversity of ancient origin, which have impacted on the variation of its indigenous livestock. Here, we investigate autosomal genome-wide profiles of 11 Ethiopian indigenous sheep populations using the Illumina Ovine 50 K SNP BeadChip assay. Sheep from the Caribbean, Europe, Middle East, China, and western, northern and southern Africa were included to address globally, the genetic variation and history of Ethiopian populations. Population relationship and structure analysis separated Ethiopian indigenous fat-tail sheep from their North African and Middle Eastern counterparts. It indicates two main genetic backgrounds and supports two distinct genetic histories for African fat-tail sheep. Within Ethiopian sheep, our results show that the short fat-tail sheep do not represent a monophyletic group. Four genetic backgrounds are present in Ethiopian indigenous sheep but at different proportions among the fat-rump and the long fat-tail sheep from western and southern Ethiopia. The Ethiopian fat-rump sheep share a genetic background with Sudanese thin-tail sheep. Genome-wide selection signature analysis identified eight putative candidate regions spanning genes influencing growth traits and fat deposition (NPR2, HINT2, SPAG8, INSR), development of limbs and skeleton, and tail formation (ALX4, HOXB13, BMP4), embryonic development of tendons, bones and cartilages (EYA2, SULF2), regulation of body temperature (TRPM8), body weight and height variation (DIS3L2), control of lipogenesis and intracellular transport of long-chain fatty acids (FABP3), the occurrence and morphology of horns (RXFP2), and response to heat stress (DNAJC18). Our findings suggest that Ethiopian fat-tail sheep represent a uniquely admixed but distinct genepool that presents an important resource for understanding the genetic control of skeletal growth, fat metabolism and associated physiological processes.