Lambskin of the Hu sheep exhibits high economic value due to its water-wave pattern. Wool curvature is the key factor of the pattern types and quality of lambskin, and it is formed by the interaction between dermal papilla cells and hair matrix cells in the hair follicle, which is regulated by various genes and signaling pathways. Herein, three full-sibling pairs of two-day-old healthy lambs (n = 6) were divided into a straight wool group (ST) and small waves group (SM) with three repetitions. RNA-seq was applied to determine the expression profile of mRNAs and lncRNAs in Hu sheep hair follicles. 25 differentially expressed mRNAs and 75 differentially expressed lncRNAs were found between SM and ST. FGF12, ATP1B4, and TCONS_00279168 were probably associated with hair follicle development. Then, Gene Ontology (GO) and KEGG enrichment analysis were implemented for the functional annotation of target genes of differentially expressed lncRNAs. The results showed that many genes, such as FGF12 and ATP1B4, were found enriched in PI3K-Akt signaling, MAPK signaling, and Ras signaling pathway associated with hair follicle growth and development. In addition, the interaction network of differentially expressed lncRNAs and mRNAs showed that a total of 6 differentially expressed lncRNAs were associated with 12 differentially expressed mRNAs, which may be as candidate mRNAs and lncRNAs. TCONS_00279168 may target ATP1B4 and FGF12 to regulate MAPK, PI3K-Akt, Ras signaling pathways involved in the sheep hair follicle development process. These results will provide the basis for exploring hair follicle development.

Sheep show characteristics of phenotypic diversity and adaptation to diverse climatic regions. Previous studies indicated associations between copy number variations (CNVs) and climate-driven adaptive evolution in humans and other domestic animals. Here, we constructed a genomic landscape of CNVs (n = 39,145) in 47 old autochthonous populations genotyped at a set of high-density (600 K) SNPs to detect environment-driven signatures of CNVs using a multivariate regression model. We found 136 deletions and 52 duplications that were significantly (P_adj. < 0.05) associated with climatic variables. These climate-mediated selective CNVs are involved in functional candidate genes for heat stress and cold climate adaptation (e.g., B3GNTL1, UBE2L3, and TRAF2), coat and wool-related traits (e.g., TMEM9, STRA6, RASGRP2, and PLA2G3), repairing damaged DNA (e.g., HTT), GTPase activity (e.g., COPG), fast metabolism (e.g., LMF2 and LPIN3), fertility and reproduction (e.g., SLC19A1 and CCDC155), growth-related traits (e.g., ADRM1 and IGFALS), and immune response (e.g., BEGAIN and RNF121) in sheep. In particular, we identified significant (P_adj. < 0.05) associations between probes in deleted/duplicated CNVs and solar radiation. Enrichment analysis of the gene sets among all the CNVs revealed significant (P_adj. < 0.05) enriched gene ontology terms and pathways related to functions such as nucleotide, protein complex, and GTPase activity. Additionally, we observed overlapping between the CNVs and 140 known sheep QTLs. Our findings imply that CNVs can serve as genomic markers for the selection of sheep adapted to specific climatic conditions.

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.

Hu sheep, a unique Chinese breed with high reproductive performance, are also well known for their rare white lambskin in China. The quality of lambskin is affected by hair follicles, and dermal papilla cells are an important component of hair follicles that plays a key role in hair follicle growth and development. This study helps elucidate the effect of miR-148a and miR-10a on hair follicle growth and development. Based on the results of gene chip and high-throughput sequencing, bone morphogenetic protein 7 (BMP7) was used as a research object. Bioinformatics analysis and the dual-luciferase reporter system indicated that, along with Western blot and quantitative real-time polymerase chain reaction (qRT-PCR) that miR-148a and miR-10a target relationships with BMP7. BMP7 was the target gene both for miR-148a and miR-10a by the dual-luciferase reporter system and Western blot. Hu sheep dermal papilla cells were successfully isolated and purified, and after transfecting miR-148a/miR-10a mimics and inhibitors into dermal papilla cells, a Cell Counting Kit-8 (CCK-8) was used to determine that miR-148a/miR-10a inhibited the proliferation of Hu sheep dermal papilla cells. In addition, after the overexpression of miR-148a, the expression levels of Smad3 (P < 0.05), Smad6 (P < 0.05), Smad4 (P < 0.01), and Smad5 (P < 0.01) were significantly higher than those of the control groups. After the inhibition of miR-148a, the expression levels of Smad3 (P < 0.05), Smad4 (P < 0.05), and TGF-β (P < 0.01) were significantly lower than those of the control groups. After the overexpression of miR-10a, the expression levels of Smad1 (P < 0.01), Smad2 (P < 0.05), Smad4 (P < 0.01), Smad5 (P < 0.01), and TGF-β (P < 0.05) were significantly lower than those of the control groups. After the inhibition of miR-10a, the expression levels of Smad1 (P < 0.01) and Smad2 (P < 0.05) were significantly lower than those of the control groups. These results revealed the target relationship between miR-148a, miR-10a and BMP7, and the effect of miR-148a and miR-10a on the proliferation of dermal papilla cells. They will provide the basis for a follow-up study on how miR-148a, and miR-10a mediate BMP7 regulation 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.

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.