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.

MicroRNA-181a Targets GNAI2 and Affects the Proliferation and Induction Ability of Dermal Papilla Cells: The Potential Involvement of the Wnt/ß-Catenin Signaling Pathway.

Wool is generated by hair follicles (HFs), which are crucial in defining the length, diameter, and morphology of wool fibers. However, the regulatory mechanism of HF growth and development remains largely unknown. Dermal papilla cells (DPCs) are a specialized cell type within HFs that play a crucial role in governing the growth and development of HFs. This study aims to investigate the proliferation and induction ability of ovine DPCs to enhance our understanding of the potential regulatory mechanisms underlying ovine HF growth and development. Previous research has demonstrated that microRNA-181a (miR-181a) was differentially expressed in skin tissues with different wool phenotypes, which indicated that miR-181a might play a crucial role in wool morphogenesis. In this study, we revealed that miR-181a inhibited the proliferation and induction ability of ovine DPCs by quantitative Real-time PCR (qRT-PCR), cell counting Kit-8 (CCK-8), 5-ethynyl-2'-deoxyuridine (EdU), flow cytometry, and alkaline phosphatase staining. Then, we also confirmed G protein subunit alpha i2 (GNAI2) is a target gene of miR-181a by dual luciferase reporter assay, qRT-PCR, and Western blot, and that it could promote the proliferation and induction ability of ovine DPCs. In addition, GNAI2 could also activate the Wnt/ß-Catenin signaling pathway in ovine DPCs. This study showed that miR-181a can inhibit the proliferation and induction ability of ovine DPCs by targeting GNAI2 through the Wnt/ß-Catenin signaling pathway.

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.

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.

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.

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.

Expression characteristics of pineal miRNAs at ovine different reproductive stages and the identification of miRNAs targeting the AANAT gene.

BACKGROUND: Many recent studies have shown that miRNAs play important roles in the regulation of animal reproduction, including seasonal reproduction. The pineal gland is a crucial hub in the regulation of seasonal reproduction. However, little is known about the expression characteristics of pineal miRNAs in different reproductive seasons (anestrus and breeding season). Therefore, the expression profiles and regulatory roles of ovine pineal miRNAs were investigated during different reproductive stages using Solexa sequencing technology and dual luciferase reporter assays. RESULTS: A total of 427 miRNAs were identified in the sheep pineal gland. Significant differences in miRNA expression were demonstrated between anestrus and the breeding season in terms of the frequency distributions of miRNA lengths, number of expressed miRNAs, and specifically and highly expressed miRNAs in each reproductive stage. KEGG analysis of the differentially expressed (DE) miRNAs between anestrus and the breeding season indicated that they are significantly enriched in pathways related to protein synthesis, secretion and uptake. Furthermore, transcriptome analysis revealed that many target genes of DE miRNAs in the ribosome pathway showed relatively low expression in the breeding season. On the other hand, analyses combining miRNA-gene expression data with target relationship validation in vitro implied that miR-89 may participate in the negative regulation of aralkylamine N-acetyltransferase (AANAT) mRNA expression by targeting its 3'UTR at a unique binding site. CONCLUSIONS: Our results provide new insights into the expression characteristics of sheep pineal miRNAs at different reproductive stages and into the negative regulatory effects of pineal miRNAs on AANAT mRNA expression.

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.

Expression Profiling and Functional Characterization of miR-26a and miR-130a in Regulating Zhongwei Goat Hair Development via the TGF-ß/SMAD Pathway.

The Zhongwei goat is an important and unique goat breed indigenous to China. It has a natural hair curling phenotype at birth, but the degree of curling gradually decreases with growth. The molecular mechanism underlying the dynamic changes in the wool curvature in Zhongwei goats is poorly understood. MicroRNAs (miRNAs) play important roles in many biological processes, including hair growth and development. In this study, we selected skins from Zhongwei goats at different ages (45 and 108 days) that exhibited different levels of hair curvature and performed miRNA sequencing to explore the molecular mechanism of hair bending. In total, 28 significantly differentially expressed miRNAs (DE miRNAs) were identified in the three groups of samples between the two developmental stages. An analysis of the target genes of the above-mentioned DE miRNAs by the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses indicated that the DE miRNAs were involved in signal pathways which were previously associated with hair bending and hair follicle development, such as the TGF-ß/SMAD, PI3K-Akt, JAK-STAT, and MAPK pathways. A comprehensive analysis of the correlations between the miRNA-seq results and issued transcriptional findings indicated that SMAD1 was a target gene of miR-26a and SMAD5 was a target gene of miR-130a. Furthermore, goat dermal papilla cells were successfully isolated and purified to determine the role of miRNAs in follicle development in vitro. The study results demonstrated that miR-130a and miR-26a had significant effects on the proliferation of dermal papilla cells. In addition, the detection results of mRNA and protein levels indicate that the overexpression of miR-26a can promote the expression of related genes in the TGF-ß/SMAD pathway, while miR-130a has the opposite substitution effect. The dual luciferase report test showed that miR-26a targeted the SMAD1 gene and reduced the expression of the SMAD1 protein in hair papillary cells. Our results identified DE microRNAs which perhaps change at the time of hair straightening in Zhongwei goats and explore the role of miR-26a and miR-130a in dermal papilla cells proliferation. The present study provided a theoretical basis to explore the mechanisms underlying the Zhongwei hair growth and curly phenotype.

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.

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.

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.