Search

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_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_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.

To understand the ovarian basis for prolificacy of Bonga sheep, a total of 31 ewes were selected based on litter size (LS) records and divided into two groups: High Prolificacy (HP) (n = 20) with LS ≥ 2 and Low Prolificacy (LP) (n = 11) with LS = 1. At a synchronized estrus, follicular dynamics were determined using transrectal ultrasonography. Plasma estradiol concentrations were also monitored. In total 27 ewes were observed in estrus being 9/11 LP (82%) and 18/20 HP (90%). On the day of estrus (day 0), the mean number of large follicles was higher (p < .05) in HP (1.78 ± 0.19) than in LP (1.0 ± 0.28) ewes. Prior to estrus, more (p < .05) medium follicles were visible for HP compared to LP ewes. Plasma estradiol concentrations were higher in HP compared to LP ewes (18.91 ± 0.41 vs. 14.51 ± 0.65 pg/ml; p < .05) and similarly was ovulation number (2.3 ± 0.15 vs. 1.28 ± 0. 14; p < .05). Higher ovulation rates and litter size in Bonga sheep are evidenced by the previous presence of more large follicles and the existence of co-dominance effects as most likely medium follicles are selected to ovulate., (© 2020 Japanese Society of Animal Science.)

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.

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

Genomic regions subjected to selection frequently show signatures such as within-population reduced nucleotide diversity and outlier values of differentiation among differentially selected populations. In this study, we analyzed 50K SNP genotype data of 373 animals belonging to 23 sheep breeds of different geographic origins using the Rsb (extended haplotype homozygosity) and FST statistical approaches, to identify loci associated with the fat-tail phenotype. We also checked if these putative selection signatures overlapped with regions of high-homozygosity (ROH). The analyses identified novel signals and confirmed the presence of selection signature in genomic regions that harbor candidate genes known to affect fat deposition. Several genomic regions that frequently appeared in ROH were also identified within each breed, but only two ROH islands overlapped with the putative selection signatures. The results reported herein provide the most complete genome-wide study of selection signatures for fat-tail in African and Eurasian sheep breeds; they also contribute insights into the genetic basis for the fat tail phenotype in sheep, and confirm the great complexity of the mechanisms that underlie quantitative traits, such as the fat-tail., Competing Interests: The authors have declared that no competing interests exist.

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

Simple Summary: Dermal papilla cells (DPCs) play a key role in the growth and development of hair follicles, and their number can regulate the shape and size of hair. Previous studies found that the BMP7 gene is differentially expressed in Hu sheep lamb skin of different pattern types and plays a key role in the proliferation of Hu sheep DPCs. However, the molecular mechanism of the differential expression of the BMP7 gene remains unclear. Here, we aim to investigate the cause of the BMP7 gene differentially expressed in Hu sheep lamb skin of different pattern types. Our findings demonstrate that DNA methylation modification can affect the expression of the BMP7 gene and can regulate the proliferation and cell cycle of Hu sheep DPCs. 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. [ABSTRACT FROM AUTHOR]

Simple Summary: Sheep diarrhea is one of the most common illnesses in sheep farming and diarrhea in sheep caused by E. coli occupies a significant proportion of them. E. coli F17 is one of the main members of E. coli causing diarrhea in sheep. This study is devoted to improving the resistance of sheep itself to E. coli F17. This study provides theoretical support for the solution of E. coli F17-caused diarrhea in sheep. 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Simple Summary: Screening for polymorphisms in the promoter region of a functional gene is an effective way to identify useful markers for improving sheep growth and development. Non-SMC condensin I complex subunit G (NCAPG) is a candidate gene linked with sheep growth and development. Its explicit role in muscle development is still unclear, and markers in NCAPG's promoter region have not been explored yet. The goal of this study was to investigate the direct role of NCAPG in regulating myogenic development and the differentiation of myoblasts and explore potential markers in its promoter region in relation to sheep growth and development traits. To achieve this goal, cell proliferation and differentiation after RNA interference with NCAPG were investigated in embryonic myoblasts. In addition, the genetic markers in the promoter region of NCAPG were scanned, and association analysis between the markers and sheep growth and development traits was carried out. The results suggest that interfering with NCAPG inhibits the proliferation and differentiation of myoblasts. Five variants detected in the promoter region of NCAPG were significantly (p < 0.05) associated with sheep growth and development traits. These results provide direct evidence of NCAPG regulating myogenesis and provide useful genetic markers to increase the efficacy of the selection of sheep growth and development traits. 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Simple Summary: Dermal papilla cells are a vital cell type in the hair follicles. It has been reported that the number of dermal papilla cells could affect the growth and development of hair follicles. Our previous study found that the CRABP2 gene is highly expressed in dermal papilla cells, but its function in dermal papilla cells is still unclear. In this study, we detected the role of the CRABP2 gene in the proliferation of dermal papilla cells and its effect on the Wnt/β-Catenin pathway. Our findings indicate that CRABP2 could promote the proliferation of dermal papilla cells by activating the Wnt/β-Catenin pathway. 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. [ABSTRACT FROM AUTHOR]

Simple Summary: Escherichia coli (E. coli) F17 is one of the major pathogenic bacteria responsible for diarrhea in farm animals; however, little is known about the biological mechanism underlying E. coli F17 infection. The aim of our study was to reveal the interplay between intestinal genes, metabolites and bacteria in E. coli F17 infected sheep. Our results confirm that the intestinal differ significantly in sheep with different E. coli F17 susceptibility, and integrated omics analyses reveal subsets of potential biomarkers for E. coli F17 infection (i.e., GlcADG 18:0-18:2, ethylmalonic acid and FBLIM1). Our results can help in the development of new insight for the treatment of farm animals infected by E. coli F17. 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Simple Summary: The Sox18 gene is important for the growth and function of dermal papilla cells, but its effect on the induction ability of dermal papilla cells in Hu sheep is not clear. Through experiments, we analyzed the effect of Sox18 on the induction ability of dermal papilla cells in Hu sheep. The results showed that Sox18 could enhance the induction ability of dermal papilla cells, and regulate the induction ability of dermal papilla cells through the Wnt/β-catenin signal pathway. We believe that Sox18 is an important gene that affects the induction ability of Hu sheep's dermal papilla cells, and these findings may provide new avenues for studying the growth and development mechanism of hair. 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Simple Summary: Diarrhea and vomiting caused by Escherichia coli (E. coli) F17 are considered significant threats to animal farming. In the present study, RNA-Seq was performed to investigate the potential circRNA and miRNA biomarkers for E. coli F17-antagonism (AN) and -sensitive (SE) lambs. The results indicated that circRNA and miRNA expression is closely associated with the susceptibility of E. coli F17 in lambs. Numbers of circRNAs and miRNAs may serve as potential biomarkers for intestinal inflammatory response against E. coli F17 infection. Our study can provide a preliminary understanding of the underlying mechanisms of intestinal immunity. 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_circ_0025840, novel_circ_0022779, novel_miR_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_circ_0000180, novel_circ_0000365, novel_miR_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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

To understand the ovarian basis for prolificacy of Bonga sheep, a total of 31 ewes were selected based on litter size (LS) records and divided into two groups: High Prolificacy (HP) (n = 20) with LS = 2 and Low Prolificacy (LP) (n = 11) with LS = 1. At a synchronized estrus, follicular dynamics were determined using transrectal ultrasonography. Plasma estradiol concentrations were also monitored. In total 27 ewes were observed in estrus being 9/11 LP (82%) and 18/20 HP (90%). On the day of estrus (day 0), the mean number of large follicles was higher (p < .05) in HP (1.78 ± 0.19) than in LP (1.0 ± 0.28) ewes. Prior to estrus, more (p < .05) medium follicles were visible for HP compared to LP ewes. Plasma estradiol concentrations were higher in HP compared to LP ewes (18.91 ± 0.41 vs. 14.51 ± 0.65 pg/ml; p < .05) and similarly was ovulation number (2.3 ± 0.15 vs. 1.28 ± 0. 14; p < .05). Higher ovulation rates and litter size in Bonga sheep are evidenced by the previous presence of more large follicles and the existence of co-dominance effects as most likely medium follicles are selected to ovulate. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Simple Summary: The meat production of Hu sheep affects the mutton supply, and how to increase the meat production is the primary focus of genetic breeding. Therefore, in order to explore the key factors affecting the proliferation and differentiation of Hu sheep skeletal muscle cells, we performed functional verification of miRNA at the cellular level. Our findings are helpful to clarify the molecular regulation mechanism of proliferation and differentiation of skeletal muscle cells, which will benefit the molecular breeding of Hu sheep. 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. [ABSTRACT FROM AUTHOR]

Simple Summary: This study was conducted to explore the antibacterial ability of sheep β-defensin 2 (SBD-2) against E. coli F17 infection in ovine intestinal epithelial cells (OIECs). Our data revealed that E. coli F17 induces SBD-2 expression in OIECs in vitro, which appears to be mediated through the activation of the signaling pathways NF-κB and MAPK. Our results provide a novel insight for the functionality of SBD-2, which could be useful for developing anti-infective drugs and/or breeding for E. coli diarrhea disease-resistant sheep. 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Through long term natural and artificial selection, domestic sheep (Ovis aries) have become adapted to a diverse range of agro-ecological environments and display multiple phenotypic traits. Characterization of diversity and selection signature is essential for genetic improvement, understanding of environmental adaptation, as well as utilization and conservation of sheep genetic resources. Here, we aimed to assess genomic diversity, population structure, and genomic selection among five Chinese native sheep breeds using 600K high density SNP genotypes. A total of 96 animals of the five breeds were selected from different geographical locations with extremely dry or humid conditions. We found a high proportion of informative SNPs, ranging from 93.3% in Yabuyi to 95.5% in Wadi, Hu, and Hetian sheep. The average pairwise population differentiation (FST) between the breeds was 0.048%, ranging from 0.022% to 0.054%, indicating their low to moderate differentiation. PCA, ADMIXTURE, and phylogenetic tree analyses revealed a clustering pattern of the five Chinese sheep breeds according to their geographical distribution, tail type, coat color, body size, and breeding history. The genomic regions under putative selection identified by FST and XP-EHH approaches frequently overlapped across the breeds, and spanned genes associated with adaptation to extremely dry or humid environments, innate and adaptive immune responses, and growth, wool, milk, and reproduction traits. The present study offers novel insight into genomic adaptation to dry and humid climates in sheep among other domestic animals and provides a valuable resource for further investigation. Moreover, it contributes useful information to sustainable utilization and conservation of sheep genetic resources. [ABSTRACT FROM AUTHOR]

Simple Summary: Bone morphogenetic protein 7 (BMP7) is one of the largest secretory signal conductive molecules and is in the TGF-β superfamily. It plays an important role in the growth and development of hair follicles. We cloned the proximal promoter of the BMP7 gene for bioinformatics analysis. Dual-luciferase reporter system and overexpression were used to analyze the key regions and transcription factor binding sites. There was high activity between −758 bp and −545 bp in the core region of the gene and a possible binding site for transcription factors SP1 and EGR1. The BMP7 gene is involved in the growth and development of hair follicles but its regulation mechanism is unclear. We studied the regulation mechanism of the BMP7 promoter by cloning the proximal promoter of BMP7 for bioinformatics analysis. A series of missing vectors was then constructed for dual-fluorescein activity detection based on the bioinformatics analysis results. We tested transcription-factor binding-site mutations and transcription factor over-expression to analyze the transcriptional regulation principle of the BMP7 promoter region. The upstream transcriptional regulatory region of the BMP7 gene proximal promoter was predicted by bioinformatics. There were −1216 bp to −1166 bp and −632 bp to −582 bp transcription initiation sites in the upstream transcriptional regulatory region of the BMP7 gene proximal promoter. The CpG islands' distribution showed that there were many CpG islands at −549 bp to 1 bp. A dual-luciferase assay revealed high activity between −758 bp and −545 bp in the core region and a possible binding site for transcription factors SP1 and EGR1. The transcriptional activity of BMP7 was significantly decreased in the transcriptional regulatory region of the BMP7 after EGR1 and SP1 mutation. Transcription was significantly enhanced by over expression of the EGR1 transcription factor, which strongly suggests that EGR1 and SP1 play important roles in BMP7 regulation. [ABSTRACT FROM AUTHOR]