Comprehensive analysis of the expression profiles of mRNA, lncRNA, circRNA, and miRNA in primary hair follicles of coarse sheep fetal skin.

BACKGROUND: The Qinghai Tibetan sheep, a local breed renowned for its long hair, has experienced significant deterioration in wool characteristics due to the absence of systematic breeding practices. Therefore, it is imperative to investigate the molecular mechanisms underlying follicle development in order to genetically enhance wool-related traits and safeguard the sustainable utilization of valuable germplasm resources. However, our understanding of the regulatory roles played by coding and non-coding RNAs in hair follicle development remains largely elusive. RESULTS: A total of 20,874 mRNAs, 25,831 circRNAs, 4087 lncRNAs, and 794 miRNAs were annotated. Among them, we identified 58 DE lncRNAs, 325 DE circRNAs, 924 DE mRNAs, and 228 DE miRNAs during the development of medullary primary hair follicle development. GO and KEGG functional enrichment analyses revealed that the JAK-STAT, TGF-ß, Hedgehog, PPAR, cGMP-PKG signaling pathway play crucial roles in regulating fibroblast and epithelial development during skin and hair follicle induction. Furthermore, the interactive network analysis additionally identified several crucial mRNA, circRNA, and lncRNA molecules associated with the process of primary hair follicle development. Ultimately, by investigating DEmir's role in the ceRNA regulatory network mechanism, we identified 113 circRNA-miRNA pairs and 14 miRNA-mRNA pairs, including IGF2BP1-miR-23-x-novel-circ-01998-MSTRG.7111.3, DPT-miR-370-y-novel-circ-005802-MSTRG.14857.1 and TSPEAR-oar-miR-370-3p-novel-circ-005802- MSTRG.10527.1. CONCLUSIONS: Our study offers novel insights into the distinct expression patterns of various transcription types during hair follicle morphogenesis, establishing a solid foundation for unraveling the molecular mechanisms that drive hair development and providing a scientific basis for selectively breeding desirable wool-related traits in this specific breed.

Genetic diversity and selection of Tibetan sheep breeds revealed by whole-genome resequencing.

OBJECTIVE: This study aimed to elucidate the underlying gene regions responsible for productive, phenotypic or adaptive traits in different ecological types of Tibetan sheep and the discovery of important genes encoding valuable traits. METHODS: We used whole-genome resequencing to explore the genetic relationships, phylogenetic tree, and population genetic structure analysis. In addition, we identified 28 representative Tibetan sheep single-nucleotide polymorphisms (SNPs) and genomic selective sweep regions with different traits in Tibetan sheep by fixation index (Fst) and the nucleotide diversity (θπ) ratio. RESULTS: The genetic relationships analysis showed that each breed partitioned into its own clades and had close genetic relationships. We also identified many potential breed-specific selective sweep regions, including genes associated with hypoxic adaptability (MTOR, TRHDE, PDK1, PTPN9, TMTC2, SOX9, EPAS1, PDGFD, SOCS3, TGFBR3), coat color (MITF, MC1R, ERCC2, TCF25, ITCH, TYR, RALY, KIT), wool traits (COL4A2, ERC2, NOTCH2, ROCK1, FGF5, SOX9), and horn phenotypes (RXFP2). In particular, a horn-related gene, RXFP2, showed the four most significantly associated SNP loci (g. 29481646 A>G, g. 29469024 T>C, g. 29462010 C>T, g. 29461968 C>T) and haplotypes. CONCLUSION: This finding demonstrates the potential for genetic markers in future molecular breeding programs to improve selection for horn phenotypes. The results will facilitate the understanding of the genetic basis of production and adaptive unique traits in Chinese indigenous Tibetan sheep taxa and offer a reference for the molecular breeding of Tibetan sheep.

Transcriptome Reveals Long Non-coding RNAs and mRNAs Involved in Primary Wool Follicle Induction in Carpet Sheep Fetal Skin.

Murine primary hair follicle induction is driven by the communication between the mesenchyme and epithelium and mostly governed by signaling pathways including wingless-related integration site (WNT), ectodysplasin A receptor (EDAR), bone morphogenetic protein (BMP), and fibroblast growth factor (FGF), as observed in genetically modified mouse models. Sheep skin may serve as a valuable system for hair research owing to the co-existence of sweat glands with wool follicles in trunk skin and asynchronized wool follicle growth pattern similar to that of human head hair follicles. However, the mechanisms underlying wool follicle development remain largely unknown. To understand how long non-coding RNAs (lncRNAs) and mRNAs function in primary wool follicle induction in carpet wool sheep, we conducted high-throughput RNA sequencing and revealed globally altered lncRNAs (36 upregulated and 26 downregulated), mRNAs (228 elevated and 225 decreased), and 80 differentially expressed novel transcripts. Several key signals in WNT (WNT2B and WNT16), BMP (BMP3, BMP4, and BMP7), EDAR (EDAR and EDARADD), and FGF (FGFR2 and FGF20) pathways, and a series of lncRNAs, including XLOC_539599, XLOC_556463, XLOC_015081, XLOC_1285606, XLOC_297809, and XLOC_764219, were shown to be potentially important for primary wool follicle induction. GO and KEGG analyses of differentially expressed mRNAs and potential targets of altered lncRNAs were both significantly enriched in morphogenesis biological processes and transforming growth factor-ß, Hedgehog, and PI3K-Akt signaling, as well as focal adhesion and extracellular matrix-receptor interactions. The prediction of mRNA-mRNA and lncRNA-mRNA interaction networks further revealed transcripts potentially involved in primary wool follicle induction. The expression patterns of mRNAs and lncRNAs of interest were validated by qRT-PCR. The localization of XLOC_297809 and XLOC_764219 both in placodes and dermal condensations was detected by in situ hybridization, indicating important roles of lncRNAs in primary wool follicle induction and skin development. This is the first report elucidating the gene network of lncRNAs and mRNAs associated with primary wool follicle early development in carpet wool sheep and will shed new light on selective wool sheep breeding.

Defining Key Genes Regulating Morphogenesis of Apocrine Sweat Gland in Sheepskin.

The apocrine sweat gland is a unique skin appendage in humans compared to mouse and chicken models. The absence of apocrine sweat glands in chicken and murine skin largely restrains further understanding of the complexity of human skin biology and skin diseases, like hircismus. Sheep may serve as an additional system for skin appendage investigation owing to the distributions and histological similarities between the apocrine sweat glands of sheep trunk skin and human armpit skin. To understand the molecular mechanisms underlying morphogenesis of apocrine sweat glands in sheepskin, transcriptome analyses were conducted to reveal 1631 differentially expressed genes that were mainly enriched in three functional groups (cellular component, molecular function and biological process), particularly in gland, epithelial, hair follicle and skin development. There were 7 Gene Ontology (GO) terms enriched in epithelial cell migration and morphogenesis of branching epithelium that were potentially correlated with the wool follicle peg elongation. An additional 5 GO terms were enriched in gland morphogenesis (20 genes), gland development (42 genes), salivary gland morphogenesis and development (8 genes), branching involved in salivary gland morphogenesis (6 genes) and mammary gland epithelial cell differentiation (4 genes). The enriched gland-related genes and two Kyoto Encyclopedia of Genes and Genomes pathway genes (WNT and TGF-ß) were potentially involved in the induction of apocrine sweat glands. Genes named BMPR1A, BMP7, SMAD4, TGFB3, WIF1, and WNT10B were selected to validate transcript expression by qRT-PCR. Immunohistochemistry was performed to localize markers for hair follicle (SOX2), skin fibroblast (PDGFRB), stem cells (SOX9) and BMP signaling (SMAD5) in sheepskin. SOX2 and PDGFRB were absent in apocrine sweat glands. SOX9 and SMAD5 were both observed in precursor cells of apocrine sweat glands and later in gland ducts. These results combined with the upregulation of BMP signaling genes indicate that apocrine sweat glands were originated from outer root sheath of primary wool follicle and positively regulated by BMP signaling. This report established the primary network regulating early development of apocrine sweat glands in sheepskin and will facilitate the further understanding of histology and pathology of apocrine sweat glands in human and companion animal skin.