IBD sharing patterns as intra-breed admixture indicators in small ruminants.

In this study, we investigated how IBD patterns shared between individuals of the same breed could be informative of its admixture level, with the underlying assumption that the most admixed breeds, i.e. the least genetically isolated, should have a much more fragmented genome. We considered 111 goat breeds (i.e. 2501 individuals) and 156 sheep breeds (i.e. 3304 individuals) from Europe, Africa and Asia, for which beadchip SNP genotypes had been performed. We inferred the breed's level of admixture from: (i) the proportion of the genome shared by breed's members (i.e. "genetic integrity level" assessed from ADMIXTURE software analyses), and (ii) the "AV index" (calculated from Reynolds' genetic distances), used as a proxy for the "genetic distinctiveness". In both goat and sheep datasets, the statistical analyses (comparison of means, Spearman correlations, LM and GAM models) revealed that the most genetically isolated breeds, also showed IBD profiles made up of more shared IBD segments, which were also longer. These results pave the way for further research that could lead to the development of admixture indicators, based on the characterization of intra-breed shared IBD segments, particularly effective as they would be independent of the knowledge of the whole genetic landscape in which the breeds evolve. Finally, by highlighting the fragmentation experienced by the genomes subjected to crossbreeding carried out over the last few generations, the study reminds us of the need to preserve local breeds and the integrity of their adaptive architectures that have been shaped over the centuries.

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.

Candidate signatures of positive selection for environmental adaptation in indigenous African cattle: A review.

Environmental adaptation traits of indigenous African cattle are increasingly being investigated to respond to the need for sustainable livestock production in the context of unpredictable climatic changes. Several studies have highlighted genomic regions under positive selection probably associated with adaptation to environmental challenges (e.g. heat stress, trypanosomiasis, tick and tick-borne diseases). However, little attention has focused on pinpointing the candidate causative variant(s) controlling the traits. This review compiled information from 22 studies on signatures of positive selection in indigenous African cattle breeds to identify regions under positive selection. We highlight some key candidate genome regions and genes of relevance to the challenges of living in extreme environments (high temperature, high altitude, high infectious disease prevalence). They include candidate genes involved in biological pathways relating to innate and adaptive immunity (e.g. BoLAs, SPAG11, IL1RL2 and GFI1B), heat stress (e.g. HSPs, SOD1 and PRLH) and hypoxia responses (e.g. BDNF and INPP4A). Notably, the highest numbers of candidate regions are found on BTA3, BTA5 and BTA7. They overlap with genes playing roles in several biological functions and pathways. These include but are not limited to growth and feed intake, cell stability, protein stability and sweat gland development. This review may further guide targeted genome studies aiming to assess the importance of candidate causative mutations, within regulatory and protein-coding genome regions, to further understand the biological mechanisms underlying African cattle's unique adaption.

Expression analysis of TLR signaling pathway genes under lipopolysaccharide-induced and E. coli F17-infected sheep intestinal epithelial cells.

Escherichia coli (E. coli) F17 is one of the main pathogens causing diarrhea in young livestock. The specific F17 fimbriae and lipopolysaccharide (LPS) in the surface components of E. coli F17 induces immune activation via interacting with the intestinal epithelial cells (IECs)-expressed innate immune toll-like receptors (TLRs) signaling pathway. In this study, the expression patterns of eight canonical genes from the TLR signaling pathway (IL-6, IL-8, IL-1ß, TLR4, MyD88, CD14, TNF-α and TRAF6) were analyzed in LPS-induced IECs, E. coli F17-infected IECs and ileum tissue of E. coli F17-infected lambs. The results showed that increased expression levels of all the studied genes were observed following post-LPS-induced and E. coli F17-infected treatment, with TLR4 having the highest up-regulated expression multiple (compared to NC, fold change = 17.94 and 20.11, respectively), and CD14 having the lowest up-regulated expression multiple (fold change = 2.68 and 1.59, respectively), and higher expression levels of all the studied TLR signaling pathway genes were observed in ileum tissue of E. coli F17 antagonistic (AN) lambs than in E. coli F17 sensitive (SE) lambs. Furthermore, when compared to LPS-induced IECs, E. coli F17-infected IECs showed a more pronounced increase in the expression of IL6, TLR4 and TNF-α, indicating the different roles of these genes in the IECs resistance to E. coli F17 infection. Our results demonstrate that the TLR signaling pathway likely promotes immune activation and provide the first evidence that TLRs have a significant potential to protect against E. coli F17 infections.

Whole-genome resequencing of Dorper and Hu sheep to reveal selection signatures associated with important traits.

Dorper and Hu sheep exhibit different characteristics in terms of reproduction, growth, and meat quality. Comparison of the genomes of two breeds help to reveal important genomic information. In this study, whole genome resequencing of 30 individuals (Dorper, DB and Hu sheep, HY) identified 15,108,125 single nucleotide polymorphisms (SNPs). Population differentiation (Fst) and cross population extended haplotype homozygosity (XP-EHH) were performed for selective signal analysis. In total, 106 and 515 overlapped genes were present in both the Fst results and XP-EHH results in HY vs DB and in DB vs HY, respectively. In HY vs DB, 106 genes were enriched in 12 GO terms and 83 KEGG pathways, such as ATP binding (GO:0005524) and PI3K-Akt signaling pathway (oas04151). In DB vs HY, 515 genes were enriched in 109 GO terms and 215 KEGG pathways, such as skeletal muscle cell differentiation (GO:0035914) and MAPK signaling pathway (oas04010). According to the annotation results, we identified a series of candidate genes associated with reproduction (UNC5C, BMPR1B, and GLIS1), meat quality (MECOM, MEF2C, and MYF6), and immunity (GMDS, GALK1, and ITGB4). Our investigation has uncovered genomic information for important traits in sheep and provided a basis for subsequent studies of related traits.

The identification and validation of target genes of IGFBP3 protein in sheep skeletal muscle cells.

This study aimed to identify the target genes of IGFBP3（insulin growth factor binding protein）protein and to investigate its target genes effects on the proliferation and differentiation of Hu sheep skeletal muscle cells. IGFBP3 was an RNA-binding protein that regulates mRNA stability. Previous studies have reported that IGFBP3 promotes the proliferation of Hu sheep skeletal muscle cells and inhibits differentiation, but the downstream genes that bind to it have not been reported yet. We predicted the target genes of IGFBP3 through RNAct and sequencing data, and verified by qPCR and RIP（RNA Immunoprecipitation）experiments, and demonstrated GNAI2（G protein subunit alpha i2）as one of the target gene of IGFBP3. After interference with siRNA, we carried out qPCR, CCK8, EdU, and immunofluorescence experiments, and found that GNAI2 can promote the proliferation and inhibit differentiation of Hu sheep skeletal muscle cells. This study revealed the effects of GNAI2 and provided one of the regulatory mechanisms of IGFBP3 protein underlying sheep muscle development.

Effects of YAP1 on proliferation and differentiation of Hu sheep skeletal muscle satellite cells in vitro.

This study aimed to understand the expression level of YAP1 in the skeletal muscle of Hu sheep and to reveal the regulatory mechanism of YAP1 on Hu sheep skeletal muscle satellite cells (SMSCs). Previous research by our group has found that YAP1 may affect the growth and development of Hu sheep skeletal muscle. In the present study, we found the expression of YAP1 in the skeletal muscle is higher than in other tissues of Hu sheep. Then, we detected the effect of YAP1 on proliferation and differentiation in Hu sheep SMSCs. According to the results of qPCR, CCK-8, EDU, and Western blot, compared to the group of negative control, overexpression of YAP1 promoted the proliferation and inhibited the differentiation of SMSCs according to the results of qPCR, CCK-8, EDU, Western blot, while the interference of YAP1 was on the contrary. Overall, our study suggests that YAP1 is an important functional molecule in the growth and development of skeletal muscle by regulating the proliferation and differentiation of SMSCs. These findings are of great use for understanding the roles of YAP1 in the skeletal muscle of Hu sheep.

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.

Genome variation in tick infestation and cryptic divergence in Tunisian indigenous sheep.

BACKGROUND: Ticks are obligate haematophagous ectoparasites considered second to mosquitos as vectors and reservoirs of multiple pathogens of global concern. Individual variation in tick infestation has been reported in indigenous sheep, but its genetic control remains unknown. RESULTS: Here, we report 397 genome-wide signatures of selection overlapping 991 genes from the analysis, using ROH, LR-GWAS, XP-EHH, and FST, of 600 K SNP genotype data from 165 Tunisian sheep showing high and low levels of tick infestations and piroplasm infections. We consider 45 signatures that are detected by consensus results of at least two methods as high-confidence selection regions. These spanned 104 genes which included immune system function genes, solute carriers and chemokine receptor. One region spanned STX5, that has been associated with tick resistance in cattle, implicating it as a prime candidate in sheep. We also observed RAB6B and TF in a high confidence candidate region that has been associated with growth traits suggesting natural selection is enhancing growth and developmental stability under tick challenge. The analysis also revealed fine-scale genome structure indicative of cryptic divergence in Tunisian sheep. CONCLUSIONS: Our findings provide a genomic reference that can enhance the understanding of the genetic architecture of tick resistance and cryptic divergence in indigenous African sheep.

Historical Introgression from Wild Relatives Enhanced Climatic Adaptation and Resistance to Pneumonia in Sheep.

How animals, particularly livestock, adapt to various climates and environments over short evolutionary time is of fundamental biological interest. Further, understanding the genetic mechanisms of adaptation in indigenous livestock populations is important for designing appropriate breeding programs to cope with the impacts of changing climate. Here, we conducted a comprehensive genomic analysis of diversity, interspecies introgression, and climate-mediated selective signatures in a global sample of sheep and their wild relatives. By examining 600K and 50K genome-wide single nucleotide polymorphism data from 3,447 samples representing 111 domestic sheep populations and 403 samples from all their seven wild relatives (argali, Asiatic mouflon, European mouflon, urial, snow sheep, bighorn, and thinhorn sheep), coupled with 88 whole-genome sequences, we detected clear signals of common introgression from wild relatives into sympatric domestic populations, thereby increasing their genomic diversities. The introgressions provided beneficial genetic variants in native populations, which were significantly associated with local climatic adaptation. We observed common introgression signals of alleles in olfactory-related genes (e.g., ADCY3 and TRPV1) and the PADI gene family including in particular PADI2, which is associated with antibacterial innate immunity. Further analyses of whole-genome sequences showed that the introgressed alleles in a specific region of PADI2 (chr2: 248,302,667-248,306,614) correlate with resistance to pneumonia. We conclude that wild introgression enhanced climatic adaptation and resistance to pneumonia in sheep. This has enabled them to adapt to varying climatic and environmental conditions after domestication.

Genomic analysis of 10 years of artificial selection in community-based breeding programs in two Ethiopian indigenous sheep breeds.

In recent times, community-based breeding programs (CBBPs) have been advocated as the best strategy for genetic improvement of local breeds in smallholder farms in developing countries. Since 2009, CBBPs have been implemented for Ethiopian Bonga and Menz sheep to improve growth rates resulting in significant genetic gains in 6-month weights. With the hypothesis that selection could be impacting their genomes, we systematically screened for possible genome changes in the two breeds by analyzing 600K BeadChip genotype data of 151 individuals (with the highest breeding values for 6-month weights) from CBBP flocks against 98 individuals from non-CBBP flocks. We observed no differences in genetic diversity and demographic dynamics between CBBP and non-CBBP flocks. Selection signature analysis employing ROH, logistic regression genome-wide association study , FST , XP-EHH and iHS revealed 5 (Bonga) and 11 (Menz) overlapping regions under selection, that co-localized with QTLs for production (body size/weight, growth, milk yield), meat/milk quality, and health/parasite resistance, suggesting that the decade-long selection has likely started to impact their genomes. However, genome-wide genetic differentiation between the CBBP and non-CBBP flocks is not yet clearly evident.

Preliminary study of melatonin in the proliferation and apoptosis of Hu sheep dermal papilla cells in vitro.

Previous studies have shown that melatonin has a certain regulatory effect on the growth of sheep wool. However, the mechanism of melatonin action remains unknown. In the present study, we aimed to understand the role of exogenous melatonin in the dermal papilla cells of Hu sheep. To confirm the optimal melatonin treatment regimen for Hu sheep dermal papilla cells, we detected the cell viability by exposing them to different concentrations of melatonin and different treatment times. The results showed that cell viability was best when dermal papilla cells were treated with 1000 pg/ml of melatonin for 48 h. According to the results of qPCR, CCK-8, EDU, Western blot, and Flow cytometry analysis, we found that 1000 pg/ml melatonin promoted the proliferation and inhibited the apoptosis of dermal papilla cells compared with the exogenous melatonin blank group (control group). Furthermore, we also found that 1000 pg/ml of melatonin promoted the cell cycle progress of dermal papilla cells according to the results of qPCR and Flow cytometry analysis. Overall, our findings showed that melatonin plays an important role in the dermal papilla cells of Hu sheep.

The wild species genome ancestry of domestic chickens.

BACKGROUND: Hybridisation and introgression play key roles in the evolutionary history of animal species. They are commonly observed within several orders in wild birds. The domestic chicken Gallus gallus domesticus is the most common livestock species. More than 65 billion chickens are raised annually to produce meat and 80 million metric tons of egg for global human consumption by the commercial sector. Unravelling the origin of its genetic diversity has major application for sustainable breeding improvement programmes. RESULTS: In this study, we report genome-wide analyses for signatures of introgression between indigenous domestic village chicken and the four wild Gallus species. We first assess the genome-wide phylogeny and divergence time across the genus Gallus. Genome-wide sequence divergence analysis supports a sister relationship between the Grey junglefowl G. sonneratii and Ceylon junglefowl G. lafayettii. Both species form a clade that is sister to the Red junglefowl G. gallus, with the Green junglefowl G. varius the most ancient lineage within the genus. We reveal extensive bidirectional introgression between the Grey junglefowl and the domestic chicken and to a much lesser extent with the Ceylon junglefowl. We identify a single case of Green junglefowl introgression. These introgressed regions include genes with biological functions related to development and immune system. CONCLUSIONS: Our study shows that while the Red junglefowl is the main ancestral species, introgressive hybridisation episodes have impacted the genome and contributed to the diversity of the domestic chicken, although likely at different levels across its geographic range.

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.

Estimates of genetic parameters and trends for reproduction traits in Bonga sheep, Ethiopia.

Investigation was carried out to assess the effects of environmental factors and to estimate genetic parameters and trends for reproductive traits in Bonga sheep, Ethiopia. Animals used in this study were managed by two communities involved in a community-based breeding programs (CBBPs) from 2009 to 2018. The database consisted of 15,595 individual phenotypic information from about 1500 ewes with variable number of records for each trait. The traits analyzed were age at first lambing (AFL), lambing interval (LI), and litter size (LS). Fixed effect analysis was done using the general linear model procedures of SAS. The Average Information Restricted Maximum Likelihood method of WOMBAT, fitting univariate animalmodel, was used to estimate heritabilities, repeatabilities and breeding values. Results showed that Bonga sheep had overall mean AFL, LI, and LS of 453 ± 109 days, 254 ± 51 days, and 1.43 ± 0.008, respectively. All traits were influenced (p < 0.01) by lambing year, lambing season, breeder communities, and parity of ewes. Estimates of heritability for AFL, LI, and LS were 0.015 ± 0.143, 0.009 ± 0.070, and 0.085 ± 0.110, respectively. The low heritabilities for the traits are expected and indicate low possibility of achieving rapid genetic progress through phenotypic selection. The repeatability estimates for LI and LS were low (0.109 and 0.196, respectively) indicating that non-genetic factors had significant influence to the variation in these traits among parities; therefore, selection decision on ewes should consider repeated records. The genetic trend for AFL, LI, and LS over the years was significant (p < 0.01). Positive values for LS and negative for AFL and LI were recorded, implying that the well-structured CBBPs have resulted in measurable genetic gains for the reproductive traits.

Correction to: Genetic diversity and matrilineal genetic origin of fat-rumped sheep in Ethiopia.

The article Genetic diversity and matrilineal genetic origin of fat-rumped sheep in Ethiopia, written by Nigussie H., Mwacharo J.M., Osama S., Agaba M., Mekasha Y., Kebede K., Abegaz S., Pal S.K., was originally published Online First without Open Access.

Genome-wide scans identify known and novel regions associated with prolificacy and reproduction traits in a sub-Saharan African indigenous sheep (Ovis aries).

Maximizing the number of offspring born per female is a key functionality trait in commercial- and/or subsistence-oriented livestock enterprises. Although the number of offspring born is closely associated with female fertility and reproductive success, the genetic control of these traits remains poorly understood in sub-Saharan Africa livestock. Using selection signature analysis performed on Ovine HD BeadChip data from the prolific Bonga sheep in Ethiopia, 41 candidate regions under selection were identified. The analysis revealed one strong selection signature on a candidate region on chromosome X spanning BMP15, suggesting this to be the primary candidate prolificacy gene in the breed. The analysis also identified several candidate regions spanning genes not reported before in prolific sheep but underlying fertility and reproduction in other species. The genes associated with female reproduction traits included SPOCK1 (age at first oestrus), GPR173 (mediator of ovarian cyclicity), HB-EGF (signalling early pregnancy success) and SMARCAL1 and HMGN3a (regulate gene expression during embryogenesis). The genes involved in male reproduction were FOXJ1 (sperm function and successful fertilization) and NME5 (spermatogenesis). We also observed genes such as PKD2L2, MAGED1 and KDM3B, which have been associated with diverse fertility traits in both sexes of other species. The results confirm the complexity of the genetic mechanisms underlying reproduction while suggesting that prolificacy in the Bonga sheep, and possibly African indigenous sheep is partly under the control of BMP15 while other genes that enhance male and female fertility are essential for reproductive fitness.

Mitochondrial and nuclear multilocus phylogeny of Rhipicephalus ticks from Kenya.

The morphological diversity of African ticks of the genus Rhipicephalus and subgenus Boophilus have been studied in detail. However, their taxonomy remains poorly resolved with limited molecular studies performed to improve inter-species discrimination. Herein, ribosomal cytochrome c oxidase I (COI), 12S ribosomal DNA (12S rDNA) and nuclear ribosomal DNA internal transcriber spacer 2 (ITS2) were analyzed in Rhipicephalus tick populations in Kenya. While the morphological and molecular criteria separated R. e. evertsi, R. pulchellus and R. appendiculatus from other members of the genus, except the morphologically similar sibling species R. zambeziensis, this was not the case for other tick populations. COI sequences of Rhipicephalus ticks from Ruma National Park (RNP) in Southwestern Kenya, that were morphologically similar to R. praetextatus/R. simus, a formed distinct clade and barcode gap group. 12S rDNA haplotypes of this population were 99% identical to a GenBank accession of R. muhsamae which is thought to be endemic in West and Central Africa. However, the ITS2 locus indicated that the RNP samples were genetically closest to ticks identified morphologically as R. praetextatus. The COI and 12S rDNA haplotype sequences of R. praetextatus clustered closely with R. simus reference sequences though the two species occurred in distinct barcode gap groups. Our results suggest that the R. simus/R. praetextatus/R. muhsamae comprise a closely related tick species complex found across sub-Saharan Africa and includes the yet to be described RNP population. More studies on the biology, ecology and genomics of all life stages of tick species in the complex may clarify their taxonomic status. A continent-wide study that combines morphology, DNA marker sequencing and emerging methods, such as mass spectrometry and whole-genome resequencing may reveal the diversity and distribution of taxa within the genus Rhipicephalus in sub-Saharan Africa.

Ancient and modern DNA reveal dynamics of domestication and cross-continental dispersal of the dromedary.

Dromedaries have been fundamental to the development of human societies in arid landscapes and for long-distance trade across hostile hot terrains for 3,000 y. Today they continue to be an important livestock resource in marginal agro-ecological zones. However, the history of dromedary domestication and the influence of ancient trading networks on their genetic structure have remained elusive. We combined ancient DNA sequences of wild and early-domesticated dromedary samples from arid regions with nuclear microsatellite and mitochondrial genotype information from 1,083 extant animals collected across the species' range. We observe little phylogeographic signal in the modern population, indicative of extensive gene flow and virtually affecting all regions except East Africa, where dromedary populations have remained relatively isolated. In agreement with archaeological findings, we identify wild dromedaries from the southeast Arabian Peninsula among the founders of the domestic dromedary gene pool. Approximate Bayesian computations further support the "restocking from the wild" hypothesis, with an initial domestication followed by introgression from individuals from wild, now-extinct populations. Compared with other livestock, which show a long history of gene flow with their wild ancestors, we find a high initial diversity relative to the native distribution of the wild ancestor on the Arabian Peninsula and to the brief coexistence of early-domesticated and wild individuals. This study also demonstrates the potential to retrieve ancient DNA sequences from osseous remains excavated in hot and dry desert environments.

Genetic diversity and matrilineal genetic origin of fat-rumped sheep in Ethiopia.

Ethiopia is home to a diverse gene pool of indigenous sheep populations. Therefore, a better understanding of genetic variation holds the key to future utilization through conservation. Three of these breeds, Afar, Blackhead Somali, and Hararghe Highland, are found in eastern Ethiopia where they contribute significantly to the livelihood of most pastoralist, agro-pastoralist, and smallholder farmers. These indigenous sheep are recognized on the basis of morphotype and their genetic distinction remains unknown. Here, to assess genetic variation, and matrilineal genetic origin and relationship of fat-rumed sheep found in eastern Ethiopia, 300 individuals from the three breeds were genotyped for 22 microsatellite markers and sequenced for the mitochondrial DNA displacement loop (mtDNA d-loop) region. The overall HO and HE were 0.57 and 0.75, respectively. Differentiation statistics revealed that a high proportion (97%) of the total genetic variation was explained by differences between individuals within populations. Genotype assignment independent of the population of origin showed K = 2 to be the optimum number of genetic backgrounds present in the dataset. This result was further confirmed by mtDNA D-loop sequences comparison in which the matrilineal genetic origin of eastern Ethiopia sheep is from two haplotype groups (types A and B) among the five haplotypes globally observed. Taken together, our findings suggest that the sheep populations from three breeds originated from two ancestral genetic backgrounds that may have diverged prior to their introduction to Ethiopia. However, to obtain a complete picture of the evolutionary dynamics of Ethiopian indigenous sheep, more samples and populations from within and outside of the country will need to be analyzed.