Joram M. Mwacharo

Multiple genomic signatures of selection in goats and sheep indigenous to a hot arid environment.

Goats and sheep are versatile domesticates that have been integrated into diverse environments and production systems. Natural and artificial selection have shaped the variation in the two species, but natural selection has played the major role among indigenous flocks. To investigate signals of natural selection, we analyzed genotype data generated using the caprine and ovine 50K SNP BeadChips from Barki goats and sheep that are indigenous to a hot arid environment in Egypt’s Coastal Zone of the Western Desert. We identify several candidate regions under selection that spanned 119 genes. A majority of the genes were involved in multiple signaling and signal transduction pathways in a wide variety of cellular and biochemical processes. In particular, selection signatures spanning several genes that directly or indirectly influenced traits for adaptation to hot arid environments, such as thermo-tolerance (melanogenesis) (FGF2, GNAI3, PLCB1), body size and development (BMP2, BMP4, GJA3, GJB2), energy and digestive metabolism (MYH, TRHDE, ALDH1A3), and nervous and autoimmune response (GRIA1, IL2, IL7, IL21, IL1R1) were identified. We also identified eight common candidate genes under selection in the two species and a shared selection signature that spanned a conserved syntenic segment to bovine chromosome 12 on caprine and ovine chromosomes 12 and 10, respectively, providing, most likely, the evidence for selection in a common environment in two different but closely related species. Our study highlights the importance of indigenous livestock as model organisms for investigating selection sweeps and genome-wide association mapping.

Signatures of positive selection in African Butana and Kenana dairy zebu cattle

Butana and Kenana are two types of zebu cattle found in Sudan. They are unique amongst African indigenous zebu cattle because of their high milk production. Aiming to understand their genome structure, we genotyped 25 individuals from each breed using the Illumina BovineHD Genotyping BeadChip. Genetic structure analysis shows that both breeds have an admixed genome composed of an even proportion of indicine (0.75 ± 0.03 in Butana, 0.76 ± 0.006 in Kenana) and taurine (0.23 ± 0.009 in Butana, 0.24 ± 0.006 in Kenana) ancestries. We also observe a proportion of 0.02 to 0.12 of European taurine ancestry in ten individuals of Butana that were sampled from cattle herds in Tamboul area suggesting local crossbreeding with exotic breeds. Signatures of selection analyses (iHS and Rsb) reveal 87 and 61 candidate positive selection regions in Butana and Kenana, respectively. These regions span genes and quantitative trait loci (QTL) associated with biological pathways that are important for adaptation to marginal environments (e.g., immunity, reproduction and heat tolerance). Trypanotolerance QTL are intersecting candidate regions in Kenana cattle indicating selection pressure acting on them, which might be associated with an unexplored level of trypanotolerance in this cattle breed. Several dairy traits QTL are overlapping the identified candidate regions in these two zebu cattle breeds. Our findings underline the potential to improve dairy production in the semi-arid pastoral areas of Africa through breeding improvement strategy of indigenous local breeds.

Genomic footprints of dryland stress adaptation in Egyptian fat-tail sheep and their divergence from East African and western Asia cohorts

African indigenous sheep are classified as fat-tail, thin-tail and fat-rump hair sheep. The fat-tail are well adapted to dryland environments, but little is known on their genome profiles. We analyzed patterns of genomic variation by genotyping, with the Ovine SNP50K microarray, 394 individuals from five populations of fat-tail sheep from a desert environment in Egypt. Comparative inferences with other East African and western Asia fat-tail and European sheep, reveal at least two phylogeographically distinct genepools of fat-tail sheep in Africa that differ from the European genepool, suggesting separate evolutionary and breeding history. We identified 24 candidate selection sweep regions, spanning 172 potentially novel and known genes, which are enriched with genes underpinning dryland adaptation physiology. In particular, we found selection sweeps spanning genes and/or pathways associated with metabolism; response to stress, ultraviolet radiation, oxidative stress and DNA damage repair; activation of immune response; regulation of reproduction, organ function and development, body size and morphology, skin and hair pigmentation, and keratinization. Our findings provide insights on the complexity of genome architecture regarding dryland stress adaptation in the fat-tail sheep and showcase the indigenous stocks as appropriate genotypes for adaptation planning to sustain livestock production and human livelihoods, under future climates.

Polymorphisms in MC1R and ASIP Genes are Associated with Coat Color Variation in the Arabian Camel

Abstract Pigmentation in mammals is primarily determined by the distribution of eumelanin and pheomelanin, the ratio of which is mostly controlled by the activity of melanocortin 1 receptor (MC1R) and agouti signaling protein (ASIP) genes. Using 91 animals from 10 Arabian camel populations, that included the 4 predominant coat color phenotypes observed in the dromedary (light brown, dark brown, black, and white), we investigated the effects of the MC1R and ASIP sequence variants and identified candidate polymorphisms associated with coat color variation. In particular, we identified a single nucleotide polymorphism (SNP), found in the coding region of MC1R (901C/T), linked to the white coat color, whereas a 1-bp deletion (23delT/T) and a SNP (25G/A) in exon 2 of ASIP are associated with both black and dark-brown coat colors. Our results also indicate support that the light-brown coat color is likely the ancestral coat color for the dromedary. These sequence variations at the MC1R and ASIP genes represent the first documented evidence of candidate polymorphisms associated with Mendelian traits in the dromedary.

Mitochondrial DNA variation reveals maternal origins and demographic dynamics of Ethiopian indigenous goats

Abstract The Horn of Africa forms one of the two main historical entry points of domestics into the continent and Ethiopia is particularly important in this regard. Through the analysis of mitochondrial DNA (mtDNA) d‐loop region in 309 individuals from 13 populations, we reveal the maternal genetic variation and demographic dynamics of Ethiopian indigenous goats. A total of 174 variable sites that generated 231 haplotypes were observed. They defined two haplogroups that were present in all the 13 study populations. Reference haplotypes from the six globally defined goat mtDNA haplogroups show the two haplogroups present in Ethiopia to be A and G, the former being the most predominant. Although both haplogroups are characterized by an increase in effective population sizes (N e) predating domestication, they also have experienced a decline in N e at different time periods, suggesting different demographic histories. We observed seven haplotypes, six were directly linked to the central haplotypes of the two haplogroups and one was central to haplogroup G. The seven haplotypes were common between Ethiopia, Kenya, Egypt, and Saudi Arabia populations, suggesting common maternal history and the introduction of goats into East Africa via Egypt and the Arabian Peninsula, respectively. While providing new mtDNA data from a historically important region, our results suggest extensive intermixing of goats mediated by human socio‐cultural and economic interactions. These have led to the coexistence of the two haplogroups in different geographic regions in Ethiopia resulting in a large caprine genetic diversity that can be exploited for genetic improvement.

Population structure, genetic diversity and selection signatures within seven indigenous Pakistani goat populations

Goat farming in Pakistan depends on indigenous breeds that have adapted to specific agro-ecological conditions. Pakistan has a rich resource of goat breeds, and the genetic diversity of these goat breeds is largely unknown. In this study, genetic diversity and population structure were characterized from seven indigenous goat breeds using the goat 50K SNP chip. The genetic diversity analysis showed that Bugi toori goats have the highest inbreeding level, consistent with the highest linkage disequilibrium, lowest diversity and long run of heterozygosity segments. This indicates that this breed should be prioritized in future conservation activities. The population structure analysis revealed four fairly distinct clusters (including Bugi toori, Bari, Black Tapri and some Kamori) and three other breeds that are seemingly the results of admixture between these or related groups (some Kamori, Pateri, Tapri and White Tapri). The selection signatures were evaluated in each breed. A total of 2508 putative selection signals were reported. The 26 significant windows were identified in more than four breeds, and selection signatures spanned several genes that directly or indirectly influence traits included coat colour variation (KIT), reproduction (BMPR1B, GNRHR, INSL6, JAK2 and EGR4), body size (SOCS2), ear size (MSRB3) and milk composition (ABCG2, SPP1, CSN1S2, CSN2, CSN3 and PROLACTIN).