Search

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.

Anthropological and biophysical processes have shaped livestock genomes over Millenia and can explain their current geographic distribution and genetic divergence. We analyzed 57 Ethiopian indigenous domestic goat genomes alongside 67 equivalents of east, west, and north-west African, European, South Asian, Middle East, and wild Bezoar goats. Cluster, ADMIXTURE (K = 4) and phylogenetic analysis revealed four genetic groups comprising African, European, South Asian, and wild Bezoar goats. The Middle Eastern goats had an admixed genome of these four genetic groups. At K = 5, the West African Dwarf and Moroccan goats were separated from East African goats demonstrating a likely historical legacy of goat arrival and dispersal into Africa via the coastal Mediterranean Sea and the Horn of Africa. F ST , XP-EHH, and Hp analysis revealed signatures of selection in Ethiopian goats overlaying genes for thermo-sensitivity, oxidative stress response, high-altitude hypoxic adaptation, reproductive fitness, pathogen defence, immunity, pigmentation, DNA repair, modulation of renal function and integrated fluid and electrolyte homeostasis. Notable examples include TRPV1 (a nociception gene); PTPMT1 (a critical hypoxia survival gene); RETREG (a regulator of reticulophagy during starvation), and WNK4 (a molecular switch for osmoregulation). These results suggest that human-mediated translocations and adaptation to contrasting environments are shaping indigenous African goat genomes., (© 2024. The Author(s).)

Previous studies have shown that the BMP7 gene is differentially expressed in Hu sheep lamb skin of different pattern types, and its expression level is significantly correlated with hair follicle indices of different pattern types, but the molecular mechanism of the differential expression of the BMP7 gene remains unclear. This study investigated the effect of DNA methylation on the transcriptional expression of BMP7 . Firstly, we found that the mRNA expression of the BMP7 gene and the activity of the core promoter of the BMP7 gene were upregulated after 5-Aza-Deoxycytidine-induced demethylation treatment using qRT-PCR and double luciferase reporter assay. Then, we found that the proliferation of Hu sheep DPCs in vitro was promoted after 5-Aza-Deoxycytidine-induced demethylation treatment through qRT-PCR, CCK-8, and EdU assay, and that the overexpression of DNMT1 in DPCs induced the opposite effect. In addition, the results of the cell cycle assay reveal that the percentage of cells in the S phase was increased after 5-Aza-Deoxycytidine-induced demethylation treatment, and that the percentage of cells in the S phase was decreased after overexpression of DNMT1 in DPCs. This study indicated that the differential expression of the BMP7 gene in different patterns of Hu sheep lamb skin may be regulated by DNA methylation modification. In addition, DNA methylation can regulate the proliferation and cell cycle of DPCs in Hu sheep.

Ethiopia has about 52 million indigenous goats with marked phenotypic variability, which is the outcome of natural and artificial selection. Here, we obtained whole-genome sequence data of three Ethiopian indigenous goat populations (Arab, Fellata, and Oromo) from northwestern Ethiopia and analyzed their genome-wide genetic diversity, population structure, and signatures of selection. We included genotype data from four other Ethiopian goat populations (Abergelle, Keffa, Gumuz, and Woyto-Guji) and goats from Asia; Europe; and eastern, southern, western, and northern Africa to investigate the genetic predisposition of the three Ethiopian populations and performed comparative genomic analysis. Genetic diversity analysis showed that Fellata goats exhibited the lowest heterozygosity values (Ho = 0.288 ± 0.005 and He = 0.334 ± 0.0001). The highest values were observed in Arab goats (Ho = 0.310 ± 0.010 and He = 0.347 ± 4.35e-05). A higher inbreeding coefficient (F ROH = 0.137 ± 0.016) was recorded for Fellata goats than the 0.105 ± 0.030 recorded for Arab and the 0.112 ± 0.034 recorded for Oromo goats. This indicates that the Fellata goat population should be prioritized in future conservation activities. The three goat populations showed the majority (∼63%) of runs of homozygosity in the shorter (100-150 Kb) length category, illustrating ancient inbreeding and/or small founder effects. Population relationship and structure analysis separated the Ethiopian indigenous goats into two distinct genetic clusters lacking phylogeographic structure. Arab, Fellata, Oromo, Abergelle, and Keffa represented one genetic cluster. Gumuz and Woyto-Guji formed a separate cluster and shared a common genetic background with the Kenyan Boran goat. Genome-wide selection signature analysis identified nine strongest regions spanning 163 genes influencing adaptation to arid and semi-arid environments ( HOXC12, HOXC13, HOXC4, HOXC6 , and HOXC9 , MAPK8IP2 ), immune response ( IL18, TYK2, ICAM3, ADGRG1, and ADGRG3 ), and production and reproduction ( RARG and DNMT1 ). Our results provide insights into a thorough understanding of genetic architecture underlying selection signatures in Ethiopian indigenous goats in a semi-arid tropical environment and deliver valuable information for goat genetic improvement, conservation strategy, genome-wide association study, and marker-assisted breeding., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Sheriff, Ahbara, Haile, Alemayehu, Han and Mwacharo.)

Diarrhea is the most common issue in sheep farms, typically due to pathogenic Escherichia coli ( E. coli ) infections, such as E. coli F17. microRNA, a primary type of non-coding RNA, has been shown to be involved in diarrhea caused by pathogenic E. coli . To elucidate the profound mechanisms of miRNA in E. coli F17 infections, methods such as E. coli F17 adhesion assay, colony counting assay, relative quantification of bacterial E. coli fimbriae gene expression, indirect immune fluorescence (IF), Cell Counting Kit-8 (CCK-8), 5-ethynyl-2'-deoxyuridine (EdU), Western blotting (WB), and scratch assay were conducted to investigate the effect of miR-329b-5p overexpression/knock-down on E. coli F17 susceptibility of sheep intestinal epithelial cells (IECs). The findings indicated that miR-329b-5p enhances the E. coli F17 resistance of sheep IECs to E.coli F17 by promoting adhesion between E. coli F17 and IEC, as well as IEC proliferation and migration. In summary, miR-329b-5p plays a crucial role in the defense of sheep IECs against E. coli F17 infection, providing valuable insights into its mechanism of action.

Domestic goats are distributed worldwide, with approximately 35% of the one billion world goat population occurring in Africa. Ethiopia has 52.5 million goats, ~99.9% of which are considered indigenous landraces deriving from animals introduced to the Horn of Africa in the distant past by nomadic herders. They have continued to be managed by smallholder farmers and semi-mobile pastoralists throughout the region. We report here 57 goat genomes from 12 Ethiopian goat populations sampled from different agro-climates. The data were generated through sequencing DNA samples on the Illumina NovaSeq 6000 platform at a mean depth of 9.71x and 150 bp pair-end reads. In total, ~2 terabytes of raw data were generated, and 99.8% of the clean reads mapped successfully against the goat reference genome assembly at a coverage of 99.6%. About 24.76 million SNPs were generated. These SNPs can be used to study the population structure and genome dynamics of goats at the country, regional, and global levels to shed light on the species' evolutionary trajectory., (© 2024. The Author(s).)

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., (© 2023. The Author(s).)

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.

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.

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.

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.

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., (© 2023 The Authors. Animal Genetics published by John Wiley & Sons Ltd on behalf of Stichting International Foundation for Animal Genetics.)

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.

Previously, NCAPG was identified as a candidate gene associated with sheep growth traits. This study aimed to investigate the direct role of NCAPG in regulating myogenesis in embryonic myoblast cells and to investigate the association between single-nucleotide polymorphisms (SNPs) in its promoter region and sheep growth traits. The function of NCAPG in myoblast proliferation and differentiation was detected after small interfering RNAs (siRNAs) knocked down the expression of NCAPG . Cell proliferation was detected using CCK-8 assay, EdU proliferation assay, and flow cytometry cell cycle analysis. Cell differentiation was detected via cell immunofluorescence and the quantification of myogenic regulatory factors (MRFs). SNPs in the promoter region were detected using Sanger sequencing and genotyped using the improved multiplex ligation detection reaction (iMLDR ® ) technique. As a result, a notable decrease ( p < 0.01) in the percentage of EdU-positive cells in the siRNA-694-treated group was observed. A significant decrease ( p < 0.01) in cell viability after treatment with siRNA-694 for 48 h and 72 h was detected using the CCK-8 method. The quantity of S-phase cells in the siRNA-694 treatment group was significantly decreased ( p < 0.01). After interfering with NCAPG in myoblasts during induced differentiation, the relative expression levels of MRFs were markedly ( p < 0.05 or p < 0.01) reduced compared with the control group on days 5-7. The myoblast differentiation in the siRNA-694 treatment group was obviously suppressed compared with the control group. SNP1, SNP2, SNP3, and SNP4 were significantly ( p < 0.05) associated with all traits except body weight measured at birth and one month of age. SNP5 was significantly ( p < 0.05) associated with body weight, body height, and body length in six-month-old sheep. In conclusion, interfering with NCAPG can inhibit the proliferation and differentiation of ovine embryonic myoblasts. SNPs in its promoter region can serve as potential useful markers for selecting sheep growth traits.

In our previous study of Hu sheep hair follicles, we found that CRABP2 was highly expressed in DPCs, which suggested that CRABP2 may influence the number of DPCs. In the present study, we aimed to understand the effect of CRABP2 in Hu sheep dermal papilla cells (DPCs). First, we explored the influence of CRABP2 on the ability of Hu sheep DPCs' proliferation. Based on the results obtained from some experiments, such as CCK-8, EDU, qPCR, and Western blot experiment, we found that the overexpression of CRABP2 facilitated the proliferation of DPCs compared to the negative control group. Then, we also detected the effect of CRABP2 on the Wnt/β-catenin pathway based on the important function of the Wnt/β-catenin pathway in hair follicles. The results showed that CRABP2 could activate the Wnt/β-catenin pathway in DPCs, and it rescues the proliferation of DPCs when the Wnt/β-catenin pathway was inhibited. In summary, our findings indicate that CRABP2 is a vital functional gene in the proliferation of Hu sheep DPCs. Our study will be of great use for revealing the roles of CRABP2 in the hair follicles of Hu sheep.

Breeding programs involving either centralized nucleus schemes and/or importation of exotic germplasm for crossbreeding were not successful and sustainable in most Africa countries. Community-based breeding programs (CBBPs) are now suggested as alternatives that aim to improve local breeds and concurrently conserve them. Community-based breeding program is unique in that it involves the different actors from the initial phase of design up until implementation of the programs, gives farmers the knowledge, skills and support they need to continue making improvements long into the future and is suitable for low input systems. In Ethiopia, we piloted CBBPs in sheep and goats, and the results show that they are technically feasible to implement, generate genetic gains in breeding goal traits and result in socio-economic impact. In Malawi, CBBPs were piloted in local goats, and results showed substantial gain in production traits of growth and carcass yields. CBBPs are currently being integrated into goat pass-on programs in few NGOs and is out-scaled to local pig production. Impressive results have also been generated from pilot CBBPs in Tanzania. From experiential monitoring and learning, their success depends on the following: 1) identification of the right beneficiaries; 2) clear framework for dissemination of improved genetics and an up/out scaling strategy; 3) institutional arrangements including establishment of breeders' cooperatives to support functionality and sustainability; 4) capacity development of the different actors on animal husbandry, breeding practices, breeding value estimation and sound financial management; 5) easy to use mobile applications for data collection and management; 6) long-term technical support mainly in data management, analysis and feedback of estimated breeding values from committed and accessible technical staff; 7) complementary services including disease prevention and control, proper feeding, and market linkages for improved genotypes and non-selected counterparts; 8) a system for certification of breeding rams/bucks to ensure quality control; 9) periodic program evaluation and impact assessment; and 10) flexibility in the implementation of the programs. Lessons relating to technical, institutional, community dynamics and the innovative approaches followed are discussed., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Haile, Getachew, Rekik, Abebe, Abate, Jimma, Mwacharo, Mueller, Belay, Solomon, Hyera, Nguluma, Gondwe and Rischkowsky.)

Escherichia coli ( E. coli ) F17 is one of the most common pathogens causing diarrhea in farm livestock. In the previous study, we accessed the transcriptomic and microbiomic profile of E. coli F17-antagonism (AN) and -sensitive (SE) lambs; however, the biological mechanism underlying E. coli F17 infection has not been fully elucidated. Therefore, the present study first analyzed the metabolite data obtained with UHPLC-MS/MS. A total of 1957 metabolites were profiled in the present study, and 11 differential metabolites were identified between E. coli F17 AN and SE lambs (i.e., FAHFAs and propionylcarnitine). Functional enrichment analyses showed that most of the identified metabolites were related to the lipid metabolism. Then, we presented a machine-learning approach (Random Forest) to integrate the microbiome, metabolome and transcriptome data, which identified subsets of potential biomarkers for E. coli F17 infection (i.e., GlcADG 18:0-18:2, ethylmalonic acid and FBLIM1 ); furthermore, the PCCs were calculated and the interaction network was constructed to gain insight into the crosstalk between the genes, metabolites and bacteria in E. coli F17 AN/SE lambs. By combing classic statistical approaches and a machine-learning approach, our results revealed subsets of metabolites, genes and bacteria that could be potentially developed as candidate biomarkers for E. coli F17 infection in lambs.

CUT-like homeobox 1 protein (CUX1), also called CUX, CUTL1, and CDP, is a member of the DNA-binding protein homology family. Studies have shown that CUX1 is a transcription factor that plays an important role in the growth and development of hair follicles. The aim of this study was to investigate the effect of CUX1 on the proliferation of Hu sheep dermal papilla cells (DPCs) to reveal the role of CUX1 in hair follicle growth and development. First, the coding sequence (CDS) of CUX1 was amplified by PCR, and then CUX1 was overexpressed and knocked down in DPCs. A Cell Counting Kit-8 (CCK8), 5-ethynyl-2-deoxyuridine (EdU), and cell cycle assays were used to detect the changes in the proliferation and cell cycle of DPCs. Finally, the effects of overexpression and knockdown of CUX1 in DPCs on the expression of WNT10 , MMP7 , C-JUN , and other key genes in the Wnt/β-catenin signaling pathway were detected by RT-qPCR. The results showed that the 2034-bp CDS of CUX1 was successfully amplified. Overexpression of CUX1 enhanced the proliferative state of DPCs, significantly increased the number of S-phase cells, and decreased the number of G0/G1-phase cells ( p < 0.05). CUX1 knockdown had the opposite effects. It was found that the expression of MMP7 , CCND1 (both p < 0.05), PPARD , and FOSL1 (both p < 0.01) increased significantly after overexpression of CUX1 in DPCs, while the expression of CTNNB1 ( p < 0.05), C-JUN , PPARD , CCND1 , and FOSL1 (all p < 0.01) decreased significantly. In conclusion, CUX1 promotes proliferation of DPCs and affects the expression of key genes of the Wnt/β-catenin signaling pathway. The present study provides a theoretical basis to elucidate the mechanism underlying hair follicle development and lambskin curl pattern formation in Hu sheep.

Sox18 is a developmental gene that encodes transcription factors. It has been indicated as be a key gene affecting the growth and development of hair follicles, in which dermal papilla cells (DPCs) have been demonstrated to play an important role through their ability to induce the formation of hair follicles. Pre-laboratory studies have found that Sox18 is differentially expressed in the dermal papilla cells of different pattern types of Hu sheep. We speculated that Sox18 plays an important role in the dermal papilla cells of Hu sheep. In our study, we analyzed the effect of Sox18 on the induction ability of DPCs in order to elucidate the function and molecular mechanism of Sox18 in the DPCs of Hu sheep. We first identified the expression of Sox18 in the DPCs of Hu sheep by immunofluorescence staining. We then used alkaline phosphatase staining, cell morphology observations and RT-PCR to detect the effect of Sox18 on the induction of DPCs after overexpression of or interference with Sox18 . We also used RT-PCR, WB and immunofluorescence staining to detect the effect of Sox18 on the Wnt/β-catenin signal pathway in DPCs. We found that Sox18 was specifically expressed in the DPCs of Hu sheep, and that Sox18 could enhance the alkaline phosphatase activity in the DPCs of Hu sheep and accelerate cell agglutination. The results of RT-PCR revealed that Sox18 promoted the mRNA expression of Versican , HHIP and FGFRI , and inhibited the mRNA expression of BMP4 and WIF1 . Further studies showed that Sox18 promoted the expression of β-catenin and activated the Wnt/β-catenin signal pathway in DPCs. When the Wnt/β-catenin signal pathway of DPCs was activated, the induction ability of DPCs was enhanced. Overall, we believe that Sox18 could enhance the induction ability of DPCs in Hu sheep and regulate the induction ability of DPCs through the Wnt/β-catenin signal pathway.

The present study aimed to investigate the activity dynamics of Ixodes ricinus group ticks in a forest located in north-western Tunisia (Aïn Draham, Jendouba District) and assess the variation of abiotic factors (temperature, Normalized Difference Vegetation Index and relative humidity) during one year survey from September 2016 to August 2017 using the dragging sampling method. A total of 116 questing ticks was collected from the vegetation consisting of 47 adults (19 females and 28 males, sex ratio M:F = 1.47), 45 nymphs and 24 larvae representing 40.5, 38.8 and 20.7% of the total collected specimens, respectively. Adult I. ricinus were collected during October-May, nymphs during May-August and larvae during July-September. There were statistically significant correlations between adult tick numbers and mean daily relative humidity (Pearson r = 0.77; p = 0.003) and mean daily temperature (r = -0.74; p = 0.006). The comparison of 16S rDNA sequences from 20 adult ticks of approximately 444 bp length showed variability among 11 sequences. There was a low genetic variability (<1%) among the I. ricinus isolates collected from the forest. The amplicons showed >99% identity with I. ricinus and Ixodes inopinatus sequences from different countries and published in GenBank. These results should be complemented by further surveys in other Tunisian regions to better understand the influence of environmental factors on the biology of I. ricinus and the occurrence of sympatric I. inopinatus ticks. Different molecular markers should be used for better understanding of their taxonomic status., Competing Interests: Declaration of Competing Interest The authors declare that they have no competing interests., (Copyright © 2022 Elsevier B.V. All rights reserved.)

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., Competing Interests: Declaration of Competing Interest The authors declare no conflicts of interest., (Crown Copyright © 2022. Published by Elsevier Inc. All rights reserved.)

Previous genome-wide association studies (GWAS) have found that LAP3 may have the potential function to impact sheep muscle development. In order to further explore whether LAP3 expression has an important role in the development of sheep embryonic myoblasts, we conducted the spatiotemporal expression profile analysis of LAP3 at the tissue and cellular level. Then we used small interfering RNA and eukaryotic recombinant vectors to perform gain/loss-of-function analysis of LAP3 . CCK-8 detection, EdU staining, and flow cytometry were used to investigate the impact of LAP3 knockdown or overexpression on the proliferation of embryonic myoblasts. In addition, cell phenotype observation, MyHC indirect immunofluorescence, and quantitative detection of the expression changes of myogenic regulatory factors (MRFs) were used to explore the effect of LAP3 on myogenic differentiation. The results showed that the LAP3 expression level in muscle tissue of fetuses was significantly higher than that in newborn lambs and adult sheep, and its expression level on day 3 of differentiation was also significantly higher than that in the proliferation phase and other differentiation time points. LAP3 silencing could significantly increase cell viability and EdU-positive cells, as well as prolonging the length of S phase of myoblasts to promote proliferation, while the results were reversed when LAP3 was overexpressed. Moreover, LAP3 silencing significantly hindered myotube formation and down-regulated the expression levels of MRFs from day 5 to day 7 of terminal differentiation, while the results were reversed when LAP3 was highly expressed. Overall, our results suggested that the expression of LAP3 impacts on the development of sheep embryonic myoblasts which provides an important theoretical basis for molecular breeding of meat production in sheep.

Early growth response factor 1 (EGR1) is a zinc-finger transcription factor that plays a vital role in the development of hair follicles. According to our previous studies, EGR1 is a transcriptional promoter of the bone morphogenetic protein 7 (BMP7), a candidate gene involved in the proliferation of dermal papilla cells. Since hair follicles are the basis of lambskin pattern formation and dermal papilla cells (DPCs) act on hair follicle growth, in order to elucidate the role of EGR1 and hair follicles, this study aimed to investigate the biological role of EGR1 in DPCs. In our study, the EGR1 coding sequence (CDS) region was firstly cloned by polymerase chain reaction, and bioinformatics analysis was performed. Then, the function of EGR1 was detected by 5-ethynyl-2’-deoxyuridine (EDU) and Cell Counting Kit-8 (CCK8), and Western blot (WB) was conducted to analyze the cellular effect of EGR1 on DPCs. The proliferative effect of EGR1 on DPCs was also further confirmed by detecting its expression by qPCR and WB on marker genes of proliferation, including PCNA and CDK2. The sequence of the EGR1 CDS region of a lamb was successfully cloned, and its nucleic acid sequence was analyzed and found to be highly homologous to Rattus norvegicus, Mus musculus, Bos taurus and Homo sapiens. Predictive analysis of the protein encoded by EGR1 revealed that it is an extra-membrane protein, and not a secretory protein, with subcellular localization in the nucleus and cytoplasm. The proliferative effect of DPCs was significantly stronger (p < 0.01) in EGR1 up-regulated DPCs compared to the controls, while the opposite result was observed in EGR1 down-regulated DPCs. Markers of proliferation including PCNA and CDK2 also appeared to be differentially upregulated in EGR1 gene overexpression compared to the controls, with the opposite result in EGR1 gene downregulation. In summary, our study revealed that EGR1 promotes the proliferation of DPCs, and we speculate that EGR1 may be closely associated with hair follicle growth and development.

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 , F ST , 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., (© 2022 Stichting International Foundation for Animal Genetics.)

It has long been recognized that enterotoxigenic Escherichia coli (ETEC) is the major pathogen responsible for vomiting and diarrhea. E. coli F17, a main subtype of ETEC, is characterized by high morbidity and mortality in young livestock. However, the transcriptomic basis underlying E. coli F17 infection has not been fully understood. In the present study, RNA sequencing was conducted to explore the expression profiles of mRNAs and long non-coding RNAs (lncRNAs) in the jejunum of lambs who were identified as resistant or sensitive to E. coli F17 that was obtained in a challenge experiment. A total of 772 differentially expressed (DE) mRNAs and 190 DE lncRNAs were detected between the E. coli F17-resistance and E. coli F17-sensitive lambs (i.e., TFF2, LOC105606142, OLFM4, LYPD8, REG4, APOA4 , TCONS&#95;00223467, and TCONS&#95;00241897). Then, a two-step machine learning approach (RX) combination Random Forest and Extreme Gradient Boosting were performed, which identified 16 mRNAs and 17 lncRNAs as potential biomarkers, within which PPP2R3A and TCONS&#95;00182693 were prioritized as key biomarkers involved in E. coli F17 infection. Furthermore, functional enrichment analysis showed that peroxisome proliferator-activated receptor (PPAR) pathway was significantly enriched in response to E. coli F17 infection. Our finding will help to improve the knowledge of the mechanisms underlying E. coli F17 infection and may provide novel targets for future treatment of E. coli F17 infection., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Chen, Lv, Zhang, Hu, Cao, Ren, Getachew, Mwacharo, Haile and Sun.)

MicroRNA (miRNA) is of great importance to muscle growth and development, including the regulation of the proliferation and differentiation of skeletal muscle satellite cells (SMSCs). In our research group's previous study, we found that miR-181a is differentially expressed in the longissimus dorsi muscle of Hu sheep at different stages. We speculated that miR-181a may participate in the growth and development process of Hu sheep. To understand the mechanism of miR-181a regulating the growth and development of Hu sheep skeletal muscle, we extracted skeletal muscle satellite cells from the longissimus dorsi muscle of 3-month-old Hu sheep fetuses and performed a series of experiments. Our results showed that miR-181a suppressed SMSCs' proliferation using QRT-PCR, Western blot, CCK-8, EDU, and Flow cytometry cycle tests. In addition, QRT-PCR, Western blot, and immunofluorescence indicated that miR-181a facilitated the differentiation of SMSCs. Then, we used dual-luciferase reporter gene detection, QRT-PCR, and Western blot to find that the Yes1-related transcription regulator (YAP1) is the target gene of miR-181a. Our study supplies a research basis for understanding the regulation mechanism of miR-181a on the growth of Hu sheep skeletal muscle.

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 F ST , 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., (© 2022. The Author(s).)

It has long been recognized that enterotoxigenic Escherichia coli (ETEC) is the major pathogen responsible for vomiting and diarrhea. E. coli F17, a main subtype of ETEC, is characterized by high morbidity and mortality in young livestock. However, the transcriptomic basis underlying E. coli F17 infection has not been fully understood. In this study, RNA sequencing was performed to explore the expression profiles of circRNAs and miRNAs in the jejunum of E. coli F17-antagonism (AN) and -sensitive (SE) lambs. A total of 16,534 circRNAs and 271 miRNAs (125 novel miRNAs and 146 annotated miRNAs) were screened, and 214 differentially expressed (DE) circRNAs and 53 DE miRNAs were detected between the AN and SE lambs (i.e., novel&#95;circ&#95;0025840, novel&#95;circ&#95;0022779, novel&#95;miR&#95;107, miR-10b). Functional enrichment analyses showed that source genes of DE circRNAs were mainly involved in metabolic-related pathways, while target genes of DE miRNAs were mainly enriched in the immune response pathways. Then, a two-step machine learning approach combining Random Forest (RF) and XGBoost (candidates were first selected by RF and further assessed by XGBoost) was performed, which identified 44 circRNAs and 39 miRNAs as potential biomarkers (i.e., novel&#95;circ&#95;0000180, novel&#95;circ&#95;0000365, novel&#95;miR&#95;192, oar-miR-496-3p) for E. coli infection. Furthermore, circRNA-related and lncRNA-related ceRNA networks were constructed, containing 46 circRNA-miRNA-mRNA competing triplets and 630 lncRNA-miRNA-mRNA competing triplets, respectively. By conducting a serious of bioinformatic analyses, our results revealed important circRNAs and miRNAs that could be potentially developed as candidate biomarkers for intestinal inflammatory response against E. coli F17 infection; our study can provide novel insights into the underlying mechanisms of intestinal immunity.

Hair follicle development is closely associated with wool curvature. Current studies reveal the crucial role of microRNAs (miRNAs) in hair follicle growth and development. However, few studies are known regarding their role in wool curvature. To reveal the potential roles of miRNAs in Hu sheep lambskin with different patterns, a total of 37 differentially expressed (DE) miRNAs were identified in hair follicles between small waves (SM) and straight wool (ST) groups using RNA-seq. Through functional enrichment and miRNA-mRNA co-expression analysis, some key miRNAs (oar-miR-143, oar-miR-200b, oar-miR-10a, oar-miR-181a, oar-miR-10b, oar-miR-125b, etc.) and miRNA-mRNA pairs (miR-125b target CD34 , miR-181a target FGF12 , LMO3 , miR-200b target ZNF536 , etc.) were identified. Though direct or indirect ways affecting hair follicle development, these miRNAs and mRNAs may have possible effects on wool curvature, and this study thus provides valuable insight on potential pattern formation.

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Bone morphogenetic proteins (BMPs) are the structurally similar and highly conserved type of functional proteins that play an important role in hair follicle growth and development. BMP7 was a differentially expressed gene in different patterns of Hu sheep lambskin identified using Agilent microarray. Since hair follicle is the basis of pattern formation of lambskin, and its growth and development is governed by dermal papilla cells (DPCs), to clarify the role of BMP7 and hair follicle, our study was designed to investigate the regulation between BMP7 and DPCs. Firstly, the CDS region of BMP7 was cloned by 3'Race and PCR in Hu sheep and performed serious of bioinformatic analysis. Then, the effects of BMP7 on DPCs were analyzed after overexpression and interference of BMP7 in dermal papilla cells by CCK8, EdU, and PI assay. Additionally, qPCR was also conducted to clarify the relationship between BMP7 and the TGF-β/Smad signaling pathway. A total of 1296 bp of the BMP7 CDS region sequence was sucessfully cloned in Hu sheep, encoding a signal peptide of 431 amino acids, molecular weight was 49,316.9 Da and the isoelectric point (Pi) was 7.75. Nucleotide sequencing analysis of BMP7 revealed that Hu sheep had high homology with Bos taurus, Homo sapiens, and Canis lupus familiaris . Structure domain prediction showed that TGF-β superfamily domain exist between 330th-431th amino acid, BMP7 protein is a secreted protein. In BMP7 up-regulated DPCs, DPCs proliferation rate and cell cycle were significantly higher than that of NC group ( p < 0.05). Meanwhile, the expression level of Smad3 , Smad4 , Samd6 , and TGF-β1 in TGF-β/Smad signaling pathway were significantly lower than that in NC group ( p < 0.05). In BMP7 down-regulated DPCs, it presented the opposite result. In conclusion, our study showed that BMP7 had a positive effect on DPCs by accelerating the proliferation and cell cycle of DPCs, and hypothesized that regulate hair follicles growth and development via TGF-β/Smad signaling pathway. These findings may provide a synergistic target for the subsequent research of hair follicle growth and development.

The growth and development of skeletal muscle require a series of regulatory factors. MiRNA is a non-coding RNA with a length of about 22 nt, which can inhibit the expression of mRNA and plays an important role in the growth and development of muscle cells. The role of miR-22-3p in C2C12 cells and porcine skeletal muscle has been reported, but it has not been verified in Hu sheep skeletal muscle. Through qPCR, CCK-8, EdU and cell cycle studies, we found that overexpression of miR-22-3p inhibited proliferation of skeletal muscle cells ( p < 0.01). The results of qPCR and immunofluorescence showed that overexpression of miR-22-3p promoted differentiation of skeletal muscle cells ( p < 0.01), while the results of inhibiting the expression of miR-22-3p were the opposite. These results suggested that miR-22-3p functions in growth and development of sheep skeletal muscle cells. Bioinformatic analysis with mirDIP, miRTargets, and RNAhybrid software suggested IGFBP3 was the target of miR-22-3p, which was confirmed by dual-luciferase reporter system assay. IGFBP3 is highly expressed in sheep skeletal muscle cells. Overexpression of IGFBP3 was found to promote proliferation of skeletal muscle cells indicated by qPCR, CCK-8, EdU, and cell cycle studies ( p < 0.01). The results of qPCR and immunofluorescence experiments proved that overexpression of IGFBP3 inhibited differentiation of skeletal muscle cells ( p < 0.01), while the results of interfering IGFBP3 with siRNA were the opposite. These results indicate that miR-22-3p is involved in proliferation and differentiation of skeletal muscle cells by targeting IGFBP3 .

Escherichia coli ( E. coli ) F17 is a member of enterotoxigenic Escherichia coli , which can cause massive diarrhea and high mortality in newborn lambs. β-defensin is mainly produced by the epithelial tissue of the gastrointestinal tract in response to microbial infection. However, the molecular mechanism of sheep β-defensin 2 ( SBD -2) against E. coli F17 remains unclear. This study aims to reveal the antibacterial ability of SBD -2 against E. coli F17 infection in sheep. Firstly, we established the culture system of ovine intestinal epithelial cells (OIECs) in vitro, treated with different concentrations of E. coli F17 for an indicated time. Secondly, we performed RNA interference and overexpression to investigate the effect of SBD -2 expression on E. coli F17 adhesion to OIECs. Finally, inhibitors of NF-κB and MAPK pathways were pre-treated to explore the possible relationship involving in E. coli F17 infection regulating SBD -2 expression. The results showed that E. coli F17 markedly ( p < 0.01) upregulated the expression levels of SBD-2 mRNA and protein in a concentration- and time-dependent manner. Overexpression of SBD -2 contributed to enhancing E. coli F17 resistance in OIECs, while silencing SBD -2 dramatically improved the adhesion of E. coli F17 to OIECs ( p < 0.05 or p < 0.01). Furthermore, E. coli F17 stimulated SBD -2 expression was obviously decreased by pre-treatment with NF-κB inhibitor PDTC, p38 MAPK inhibitor SB202190 and ERK1/2 MAPK inhibitor PD98095 ( p < 0.05 or p < 0.01). Interestingly, adhesion of E. coli F17 to OIECs were highly enhanced by pre-treated with PDTC, SB202190 and PD98095. Our data suggested that SBD -2 could inhibit E. coli F17 infection in OIECs, possibly through NF-κB and MAPK signaling pathways. Our results provide useful theoretical basis on developing anti-infective drug and breeding for E. coli diarrhea disease-resistant sheep.

Wool curvature is the determining factor for lambskin quality of Hu lambs. However, the molecular mechanism of wool curvature formation is not yet known. miRNA has been proved to play an important role in hair follicle development, and we have discovered a differentially expressed miRNA, miR-143, in hair follicles of different curl levels. In this study, we first examined the effects of miR-143 on the proliferation and cell cycle of dermal papilla cells using CCK8, EdU and flow cytometry and showed that miR-143 inhibited the proliferation of dermal papilla cells and slowed down the cell cycle. Bioinformatics analysis was performed to predict the target genes KRT71 and CUX1 of miR-143, and both two genes were expressed at significantly higher levels in small waves than in straight lambskin wool ( p < 0.05) as detected by qPCR and Western blot (WB). Then, the target relationships between miR-143 and KRT71 and CUX1 were verified through the dual-luciferase assay in 293T cells. Finally, after overexpression and suppression of miR-143 in dermal papilla cells, the expression trend of CUX1 was contrary to that of miR-143. Meanwhile, KRT71 was not detected because KRT71 was not expressed in dermal papilla cells. Therefore, we speculated that miR-143 can target CUX1 to inhibit the proliferation of dermal papilla cells, while miR-143 can target KRT71 to regulate the growth and development of hair follicles, so as to affect the development of hair follicles and ultimately affect the formation of wool curvature.

Theileria annulata , which causes tropical theileriosis, is a major impediment to improving cattle production in Sudan. Tropical theileriosis disease is prevalent in the north and central regions of Sudan. Outbreaks of the disease have been observed outside the known endemic areas, in east and west regions of the country, due to changes in tick vector distribution and animal movement. A live schizont attenuated vaccination based on tissue culture technology has been developed to control the disease. The parasite in the field as well as the vaccine strain need to be genotyped before the vaccinations are practiced, in order to be able to monitor any breakthrough or breakdown, if any, after the deployment of the vaccine in the field. Nine microsatellite markers were used to genotype 246 field samples positive for T. annulata DNA and the vaccine strain. North and central populations have a higher multiplicity of infection than east and west populations. The examination of principal components showed two sub-structures with a mix of all four populations in both clusters and the vaccine strain used being aligned with left-lower cluster. Only the north population was in linkage equilibrium, while the other populations were in linkage disequilibrium, and linkage equilibrium was found when all samples were regarded as single population. The genetic identity of the vaccine and field samples was 0.62 with the north population and 0.39 with west population. Overall, genetic investigations of four T. annulata populations in Sudan revealed substantial intermixing, with only two groups exhibiting regional origin independence. In the four geographically distant regions analyzed, there was a high level of genetic variation within each population. The findings show that the live schizont attenuated vaccine, Atbara strain may be acceptable for use in all Sudanese regions where tropical theileriosis occurs., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Salih, Ali, Njahira, Taha, Mohammed, Mwacharo, Mbole-Kariuki, El Hussein, Bishop and Skilton.)

Autochthonous Sudanese cattle breeds, namely Baggara for beef and Butana and Kenana for dairy, are characterized by their adaptive characteristics and high performance in hot and dry agro-ecosystems. They are thus used largely by nomadic and semi-nomadic pastoralists. We analyzed the diversity and genetic structure of the BoLA-DRB3 gene, a genetic locus linked to the immune response, for the indigenous cattle of Sudan and in the context of the global cattle repository. Blood samples (n = 225) were taken from three indigenous breeds (Baggara; n = 113, Butana; n = 60 and Kenana; n = 52) distributed across six regions of Sudan. Nucleotide sequences were genotyped using the sequence-based typing method. We describe 53 alleles, including seven novel alleles. Principal component analysis (PCA) of the protein pockets implicated in the antigen-binding function of the MHC complex revealed that pockets 4 and 9 (respectively) differentiate Kenana-Baggara and Kenana-Butana breeds from other breeds. Venn analysis of Sudanese, Southeast Asian, European and American cattle breeds with 115 alleles showed 14 were unique to Sudanese breeds. Gene frequency distributions of Baggara cattle showed an even distribution suggesting balancing selection, while the selection index (ω) revealed the presence of diversifying selection in several amino acid sites along the BoLA-DRB3 exon 2 of these native breeds. The results of several PCA were in agreement with clustering patterns observed on the neighbor joining (NJ) trees. These results provide insight into their high survival rate for different tropical diseases and their reproductive capacity in Sudan's harsh environment., (© 2021. The Author(s).)

With climate change bound to affect food and feed production, emphasis will shift to resilient and adapted indigenous livestock to sustain animal production. However, indigenous livestock comprise several varieties, strains and ecotypes whose genomes are poorly characterized. Here, we investigated genomic variation in an African thin-tailed Desert Sheep sampled in Sudan, using 600K genotype data generated from 92 individuals representing five ecotypes. We included data from 18 fat-tailed and 45 thin-tailed sheep from China, to investigate shared ancestry and perform comparative genomic analysis. We observed a clear genomic differentiation between the African thin-tailed Desert Sheep and the Chinese thin-tailed and fat-tailed sheep, suggesting a broad genetic structure between the fat-tailed and thin-tailed sheep in general, and that at least two autosomal gene pools comprise the genome profile of the thin-tailed sheep. Further analysis detected two distinct genetic clusters in both the African thin-tailed Desert Sheep and the Chinese thin-tailed sheep, suggesting a fine-scale and complex genome architecture in thin-tailed sheep. Selection signature analysis suggested differences in adaptation, production, reproduction and morphology likely underly the fine-scale genetic structure in the African thin-tailed Desert Sheep. This may need to be considered in designing breeding programs and genome-wide association studies., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Abied, Ahbara, Berihulay, Xu, Islam, El-Hag, Rekik, Haile, Han, Ma, Zhao and Mwacharo.)

Knowledge on how adaptive evolution and human socio-cultural and economic interests shaped livestock genomes particularly in sub-Saharan Africa remains limited. Ethiopia is in a geographic region that has been critical in the history of African agriculture with ancient and diverse human ethnicity and bio-climatic conditions. Using 52K genome-wide data analysed in 646 individuals from 13 Ethiopian indigenous goat populations, we observed high levels of genetic variation. Although runs of homozygosity (ROH) were ubiquitous genome-wide, there were clear differences in patterns of ROH length and abundance and in effective population sizes illustrating differences in genome homozygosity, evolutionary history, and management. Phylogenetic analysis incorporating patterns of genetic differentiation and gene flow with ancestry modelling highlighted past and recent intermixing and possible two deep ancient genetic ancestries that could have been brought by humans with the first introduction of goats in Africa. We observed four strong selection signatures that were specific to Arsi-Bale and Nubian goats. These signatures overlapped genomic regions with genes associated with morphological, adaptation, reproduction and production traits due possibly to selection under environmental constraints and/or human preferences. The regions also overlapped uncharacterized genes, calling for a comprehensive annotation of the goat genome. Our results provide insights into mechanisms leading to genome variation and differentiation in sub-Saharan Africa indigenous goats., Competing Interests: The authors declare no conflict of interest regarding this manuscript., (© 2020 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

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., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The identification of genome-wide selection signatures can provide insights on the mechanisms of natural and/or artificial selection and uncover genes related to biological functions and/or phenotypes. Tibetan sheep are an important livestock in Tibet, providing meat and wool for Tibetans who are renown for breeding livestock that adapt well to high altitudes. Using whole-genome sequences with an effective sequencing depth of 5×, we investigated the genomic diversity and structure and, identified selection signatures of White Tibetan, Oula and Poll Dorset sheep. We obtained 30,163,679 Single Nucleotide Polymorphisms (SNPs) and 5,388,372 indels benchmarked against the ovine Oar&#95;v4.0 genome assembly. Next, using F ST , ZHp and XP-EHH approaches, we identified selection signatures spanning a set of candidate genes, including HIF1A, CAPN3, PRKAA1, RXFP2, TRHR and HOXA10 that are associated with pathways and GO categories putatively related to hypoxia responses, meat traits and disease resistance. Candidate genes and GO terms associated with coat color were also identified. Finally, quantification of blood physiological parameters, revealed higher levels of mean corpuscular hemoglobin measurement and mean corpuscular hemoglobin concentration in Tibetan sheep compared with Poll Dorset, suggesting a greater oxygen-carrying capacity in the Tibetan sheep and thus better adaptation to high-altitude hypoxia. In conclusion, this study provides a greater understanding of genome diversity and variations associated with adaptive and production traits in sheep.

Homozygosity of long sequence genotypes are a result of parents transmitting identical haplotypes, which can be used to estimate their auto-zygosity. Therefore, we used high-density SNP Chip data to characterize the auto-zygosity of each breed according to the occurrence and distribution of runs of homozygosity (ROH). Subsequently, we identified the genomic regions with high runs of homozygosity frequencies within individuals of each breed. We selected 96 sheep samples from five local Chinese sheep breeds belonging to different geographical locations. We identified 3046 ROHs within the study breed individuals, among which the longer segments (>1-5 Mb) were dominant. On average, ROH segments covered about 12% of the genomes; the coverage rate of OAR20 was the lowest and that of OAR2 was the highest. The distribution analysis of runs of homozygosity showed that the detected ROH mainly distributed between >26 and 28 Mb. The Hetian and Hu sheep showed the lowest ROH distribution. The estimation of homozygosity level reflects the history of modern and ancient inbreeding, which may affect the genomes of Chinese indigenous sheep breeds and indicate that some animals have experienced recent self-pollination events (Yabuyi, Karakul and Wadi). In these sheep breeds, the genomic regions were assumed to be under selection signatures frequently in line with long ROH. These regions included candidate genes associated with disease resistance traits ( 5S&#95;rRNA ), the innate and adaptive immune response ( HERC2 and CYFIP1 ), digestion and metabolism ( CENPJ ), growth (SPP1), body size and developments ( GJB2 and GJA3 ). This study highlighted new insights into the ROH patterns and provides a basis for future breeding and conservation strategies of Chinese sheep breeds.

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.

In Tunisia, livestock plays an important role in the economy; unfortunately, Tunisian sheep population faces several health challenges. The aim of this trial was to study phenotypic variability of four local sheep breeds and strains to abomasum nematodes infection and to Babesia/Theileria parasites. Faeces, blood and abomasum contents were collected from 310 sheep slaughtered in eight commercial slaughterhouses across North Tunisia. Haematological and biochemical parameters were assessed. DNA was extracted and catch-all primers were used to detect both Theileria spp. and Babesia spp. DNA. Faecal egg counts (FEC) was quantitatively assessed using simple flotation technique followed by McMaster technique. Male and female worms were collected from all abomasum contents and counted under a stereomicroscope. The percentage of faeces samples positive for GIN's eggs was 30.82%. After worms' recovery, the infection prevalence was estimated to 75.90%. The overall infection prevalence by Babesia spp. and Theileria spp. was 4.21%. The dispersion of observations plots obtained by principal component analysis (PCA) showed two clusters of individuals. The first cluster contains animals having positive Babesia/Theileria PCR, presence of nematodes in the abomasum contents and relatively low total worm count (TWC < 500) expect one animal which was found bearing high TWC (>500). In this same group, with a suspected form of resistance, animals showed normal values of albumin and normal haematological parameters (red blood cell count [RBC], haemoglobin [Hb] and packed cell volume [PCV]). The second cluster represents all the other observations in which subgroups of animals were distinguished on the basis of their potential resistance to abomasum nematodes. Multiple correlations showed significant positive correlations between RBC/Hb, RBC/PCV, PCV/Hb and FEC/TWC. Significant negative correlations were observed between TWC/RBC and TWC/Hb. It is concluded that the phenotypic variability among local sheep breeds is essential for more advanced genetic and genomic studies., (© 2020 The Authors. Veterinary Medicine and Science published by John Wiley & Sons Ltd.)

The domestication and subsequent global dispersal of livestock are crucial events in human history, but the migratory episodes during the history of livestock remain poorly documented [1-3]. Here, we first developed a set of 493 novel ovine SNPs of the male-specific region of Y chromosome (MSY) by genome mapping. We then conducted a comprehensive genomic analysis of Y chromosome, mitochondrial DNA, and whole-genome sequence variations in a large number of 595 rams representing 118 domestic populations across the world. We detected four different paternal lineages of domestic sheep and resolved, at the global level, their paternal origins and differentiation. In Northern European breeds, several of which have retained primitive traits (e.g., a small body size and short or thin tails), and fat-tailed sheep, we found an overrepresentation of MSY lineages y-HC and y-HB, respectively. Using an approximate Bayesian computation approach, we reconstruct the demographic expansions associated with the segregation of primitive and fat-tailed phenotypes. These results together with archaeological evidence and historical data suggested the first expansion of early domestic hair sheep and the later expansion of fat-tailed sheep occurred ∼11,800-9,000 years BP and ∼5,300-1,700 years BP, respectively. These findings provide important insights into the history of migration and pastoralism of sheep across the Old World, which was associated with different breeding goals during the Neolithic agricultural revolution., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.)

Understanding the genetic changes underlying phenotypic variation in sheep (Ovis aries) may facilitate our efforts towards further improvement. Here, we report the deep resequencing of 248 sheep including the wild ancestor (O. orientalis), landraces, and improved breeds. We explored the sheep variome and selection signatures. We detected genomic regions harboring genes associated with distinct morphological and agronomic traits, which may be past and potential future targets of domestication, breeding, and selection. Furthermore, we found non-synonymous mutations in a set of plausible candidate genes and significant differences in their allele frequency distributions across breeds. We identified PDGFD as a likely causal gene for fat deposition in the tails of sheep through transcriptome, RT-PCR, qPCR, and Western blot analyses. Our results provide insights into the demographic history of sheep and a valuable genomic resource for future genetic studies and improved genome-assisted breeding of sheep and other domestic animals.

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.

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.