Integrated Hair Follicle Profiles of microRNAs and mRNAs to Reveal the Pattern Formation of Hu Sheep Lambskin

Identification of miR-141 as a Regulator of Epidermal Homeostasis

MicroRNAs, small endogenous noncoding RNAs, are involved in the regulation of epidermal homeostasis. Among them, miR-203 was the most described and expressed in human epidermis, promoting keratinocyte (KC) differentiation by repressing genes involved in proliferation. To identify other microRNAs involved in this process, the miRNomes of normal human KCs cultured in monolayer (2-dimensional) or in 3-dimensional reconstructed skin were compared. Besides miR-203, miR-141 was one of the most expressed microRNAs in 3-dimensional culture and was overexpressed in 3-dimensional versus 2-dimensional condition, that is, during KC differentiation. Functional experiments revealed that, mostly expressed in the basal layer, miR-141 decreased KC proliferation and clonogenicity while promoting differentiation. Target prediction algorithm coupled with transcriptomic data of KCs overexpressing miR-141 as well as 3' untranslated region luciferase assays highlighted CCND2 mRNA as a direct target of miR-141, leading to its downregulation by miR-141 overexpression. Finally, CCND2 silencing decreased KC proliferation and induced differentiation, revealing that miR-141 action was mediated by CCND2. MiR-141 features were also compared with those of miR-203 in parallel experiments. Although miR-141 displayed functions similar to those of miR-203, it exhibited different localization and targets, suggesting a joint participation of miR-141 and miR-203 in engaging and maintaining KC toward differentiation, respectively.

Potential Involvement of miR-144 in the Regulation of Hair Follicle Development and Cycle Through Interaction with Lhx2

BACKGROUND

Cashmere, known as "soft gold", is a highly prized fiber from Cashmere goats, produced by secondary hair follicles. Dermal papilla cells, located at the base of these follicles, regulate the proliferation and differentiation of hair matrix cells, which are essential for hair growth and cashmere formation. Recent studies emphasize the role of microRNAs (miRNAs) in controlling gene expression within these processes.

METHODS

This study centered on exploring the targeted regulatory interaction between miR-144 and the Lhx2 gene. Utilizing methodologies like miRNA target prediction, luciferase reporter assays, and quantitative PCR, they assessed the interplay between miR-144 and Lhx2. Dermal papilla cells derived from Cashmere goats were cultured and transfected with either miR-144 mimics or inhibitors to observe the subsequent effects on Lhx2 expression.

RESULTS

The results demonstrated that miR-144 directly targets the Lhx2 gene by binding to its mRNA, leading to a decrease in Lhx2 expression. This modulation of Lhx2 levels influenced the behavior of dermal papilla cells, affecting their ability to regulate hair matrix cell proliferation and differentiation. Consequently, the manipulation of miR-144 levels had a significant impact on the growth cycle of cashmere wool.

CONCLUSIONS

The findings suggest miR-144 regulates hair follicle dynamics by targeting Lhx2, offering insights into hair growth mechanisms. This could lead to innovations in enhancing cashmere production, fleece quality, and addressing hair growth disorders. Future research may focus on adjusting miR-144 levels to optimize Lhx2 expression and promote hair follicle activity.

Melatonin promotes the proliferation of dermal papilla cells in cashmere goats via activation of chi-let-7d-5p/WNT2 axis

Exogenous melatonin promotes the differentiation of secondary hair follicles in Cashmere goats, thereby improving cashmere production. MicroRNAs (miRNAs) play a crucial role in regulating post-transcriptional gene expression and influence hair follicle growth. However, the mechanism through which melatonin regulates hair follicle development via miRNA mediation remains unclear. In this study, we used RNA-seq to identify differentially expressed (DE) miRNAs during melatonin-induced growth of secondary hair follicles in inner Mongolian Cashmere goats. In total, 170 DE miRNAs were identified. Enrichment analysis revealed that the target genes of these DE miRNAs were related to biological processes such as protein modification; cytoskeletal components; and the Notch, Wnt, and MAPK signaling pathways. The miRNA-mRNA regulatory network suggested that the DE miRNA chi-let-7d-5p negatively regulates WNT2 expression. Mechanistic studies revealed that melatonin promotes the proliferation of DP cells in Cashmere goats via the chi-let-7d-5p/WNT2 axis.

Population structure and selective signature of Kirghiz sheep by Illumina Ovine SNP50 BeadChip

Objective

By assessing the genetic diversity and associated selective traits of Kirghiz sheep (KIR), we aim to uncover the mechanisms that contribute to sheep's adaptability to the Pamir Plateau environment.

Methods

This study utilized Illumina Ovine SNP50 BeadChip data from KIR residing in the Pamir Plateau, Qira Black sheep (QBS) inhabiting the Taklamakan Desert, and commonly introduced breeds including Dorper sheep (DOR), Suffolk sheep (SUF), and Hu sheep (HU). The data was analyzed using principal component analysis, phylogenetic analysis, population admixture analysis, kinship matrix analysis, linkage disequilibrium analysis, and selective signature analysis. We employed four methods for selective signature analysis: fixation index (Fst), cross-population extended homozygosity (XP-EHH), integrated haplotype score (iHS), and nucleotide diversity (Pi). These methods aim to uncover the genetic mechanisms underlying the germplasm resources of Kirghiz sheep, enhance their production traits, and explore their adaptation to challenging environmental conditions.

Results

The test results unveiled potential selective signals associated with adaptive traits and growth characteristics in sheep under harsh environmental conditions, and annotated the corresponding genes accordingly. These genes encompass various functionalities such as adaptations associated with plateau, cold, and arid environment (ETAA1, UBE3D, TLE4, NXPH1, MAT2B, PPARGC1A, VEGFA, TBX15 and PLXNA4), wool traits (LMO3, TRPS1, EPHA5), body size traits (PLXNA2, EFNA5), reproductive traits (PPP3CA, PDHA2, NTRK2), and immunity (GATA3).

Conclusion

Our study identified candidate genes associated with the production traits and adaptation to the harsh environment of the Pamir Plateau in Kirghiz sheep. These findings provide valuable resources for local sheep breeding programs. The objective of this study is to offer valuable insights for the sustainable development of the Kirghiz sheep industry.

Analysis of the Long Non-Coding and Messenger RNA Expression Profiles in the Skin Tissue of Super Merino and Small-Tailed Han Sheep

Wool quality and yield are two important economic livestock traits. However, there are relatively few molecular studies on lncRNA for improving sheep wool, so these require further exploration. In this study, we examined skin tissue from the upper scapula of Super Merino (SM) and Small-Tailed Han (STH) sheep during the growing period. The apparent difference was verified via histological examination. High-throughput RNA sequencing identified differentially expressed (DE) long non-coding (lncRNAs) and messenger RNAs (mRNAs). The target gene of DE lncRNA and DE genes were enrichment analyzed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). A Reverse Transcription quantitative Polymerase Chain Reaction (RT-qPCR) was used to verify randomly selected DE lncRNAs and mRNAs. Finally, the DE, RAC2, WNT11, and FZD2 genes, which were enriched in the Wnt signaling pathway, were detected via immunohistochemistry. The results showed that a total of 20,888 lncRNAs and 31,579 mRNAs were identified in the skin tissues of the two sheep species. Among these, 56 lncRNAs and 616 mRNAs were differentially expressed. Through qRT-PCR, the trends in the randomly selected DE genes' expression were confirmed to be aligned with the RNA-seq results. GO and KEGG enrichment analysis showed that DE lncRNA target genes were enriched in GO terms as represented by epidermal and skin development and keratin filature and in KEGG terms as represented by PI3K-Akt, Ras, MAPK, and Wnt signaling pathways, which were related to hair follicle growth and development. Finally, immunohistochemistry staining results indicated that RAC2, WNT11, and FZD2 were expressed in dermal papilla (DP). The lncRNAs MSTRG.9225.1 and MSTRG.98769.1 may indirectly participate in the regulation of hair follicle growth, development, and fiber traits by regulating their respective target genes, LOC114113396(KRTAP15-1), FGF1, and IGF1. In addition, MSTRG.84658.1 may regulate the Wnt signaling pathway involved in the development of sheep hair follicles by targeting RAC2. This study provides a theoretical reference for improving sheep breeding in the future and lays a foundation for further research on the effects of MSTRG.84658.1 and the target gene RAC2 on dermal papilla cells (DPC).

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

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

Identification of Potential miRNA-mRNA Regulatory Network Associated with Growth and Development of Hair Follicles in Forest Musk Deer

In this study, sRNA libraries and mRNA libraries of HFs of FMD were constructed and sequenced using an Illumina HiSeq 2500, and the expression profiles of miRNAs and genes in the HFs of FMD were obtained at the anagen and catagen stages. In total, 565 differentially expressed unigenes (DEGs) were identified, 90 of which were upregulated and 475 of which were downregulated. In the BP category of GO enrichment, the DEGs were enriched in the processes related to HF development and differentiation, including the hair cycle regulation and processes, HF development, skin epidermis development, regulation of HF development, skin development, the Wnt signaling pathway, and the BMP signaling pathway. Through KEGG analysis it was found that DEGs were significantly enriched in pathways associated with HF development and growth. A total of 186 differentially expressed miRNAs (DEmiRNAs) were screened (p < 0.05) in the HFs of FMD at the anagen stage vs. the catagen stage, 33 of which were upregulated and 153 of which were downregulated. Through DEmiRNA-mRNA association analysis, we found DEmiRNAs and target genes that mainly play regulatory roles in HF development and growth. The enrichment analysis of DEmiRNA target genes revealed similarities with the enrichment results of DEGs associated with HF development. Notably, both sets of genes were enriched in key pathways such as the Notch signaling pathway, melanogenesis, the cAMP signaling pathway, and cGMP-PKG. To validate our findings, we selected 11 DEGs and 11 DEmiRNAs for experimental verification using RT-qPCR. The results of the experimental validation were consistent with the RNA-Seq results.

Imprinted Dlk1-Gtl2 cluster miRNAs are potential epigenetic regulators of lamb fur quality

BACKGROUND

Tan and Hu sheep are well-known local breeds in China, producing lamb fur with unique ornamental and practical values highly appreciated by consumers worldwide. Fur quality is optimal at one month of age and gradually declines with time. Despite active research on its genetic mechanism using transcriptomic and whole genome bisulfite sequencing analysis, the main effective gene locus has not been found, and its regulatory mechanism is still unclear, which limits the breeding and improvement of fur traits.

RESULTS

Scapular skin samples from newborn (1-month old) and adult (24-month old) Tan sheep were utilized for small ribonucleic acid (RNA) sequencing Principal Component Analysis (PCA) showed that the newborn and adult groups were completely separated. Differential expression analysis of micro-RNAs (miRNAs) identified 32 up-regulated miRNAs and 48 down-regulated miRNAs in the newborn groups. All up-regulated miRNAs were located in the imprinted. Dlk1-Gtl2 locus on chromosome 18, whereas all down-regulated miRNAs were distributed across the sheep chromosomes, without a clear pattern of positional consistency. Further, by systematically analyzing the target genes and signaling pathways of all 32 up-regulated miRNAs, we found that the PI3K-AKT signaling pathway has the potential to be targeted and regulated by most of the miRNAs in the Dlk1-Gtl2 region. In addition, we also re-analyzed miRNA sequencing data from public databases on Hu lambs (full sibling Hu lambs with high- and low-quality fur characteristics). Again, it was found that most of the up-regulated miRNAs in lambs with high-quality fur were also located in the Dlk1-Gtl2 region, whereas this patter was not present for down-regulated miRNAs.

CONCLUSION

Sequencing of miRNAs in conjunction with public databases was employed to identify miRNAs within the imprinted Dlk1-Gtl2 region on chromosome 18, suggesting their potential roles as epigenetic regulators of fur traits. Small RNAs located at the Dlk1-Gtl2 locus were identified as having the potential to systematically regulate the PI3K-AKT signaling pathway, thereby indicating the relevance of the Dlk1-Gtl2/PI3K-AKT axis in the context of fur traits. Selection of parental specific expressed imprinted genes in the process of conserving and exploiting lamb fur traits should be emphasized.

Integrated analysis of miRNAs and mRNA profiling reveals the potential roles of miRNAs in sheep hair follicle development

BACKGROUND

Merino sheep exhibit high wool production and excellent wool quality. The fleece of Merino sheep is predominantly composed of wool fibers grown from hair follicles (HFs). The HF is a complex biological system involved in a dynamic process governed by gene regulation, and gene expression is regulated by microRNAs (miRNAs). miRNA inhibits posttranscriptional gene expression by specifically binding to target messenger RNA (mRNA) and plays an important role in regulating gene expression, the cell cycle and biological development sequences. The purpose of this study was to examine mRNA and miRNA binding to identify key miRNAs and target genes related to HF development. This will provide new and important insights into fundamental mechanisms that regulate cellular activity and cell fate decisions within and outside of the skin.

RESULTS

We analyzed miRNA data in skin tissues collected from 18 Merino sheep on four embryonic days (E65, E85, E105 and E135) and two postnatal days (D7 and D30) and identified 87 differentially expressed miRNAs (DE-miRNAs). These six stages were further divided into two longer developmental stages based on heatmap cluster analysis, and the results showed that DE-mRNAs in Stage A were closely related to HF morphogenesis. A coanalysis of Stage A DE-mRNAs and DE-miRNAs revealed that 9 DE-miRNAs and 17 DE-mRNAs presented targeting relationships in Stage A. We found that miR-23b and miR-133 could target and regulate ACVR1B and WNT10A. In dermal fibroblasts, the overexpression of miR-133 significantly reduced the mRNA and protein expression levels of ACVR1B. The overexpression of miR-23b significantly reduced the mRNA and protein expression levels of WNT10A.

CONCLUSION

This study provides a new reference for understanding the molecular basis of HF development and lays a foundation for further improving sheep HF breeding. miRNAs and target genes related to hair follicular development were found, which provided a theoretical basis for molecular breeding for the culture of fine-wool sheep.

Genetics of the phenotypic evolution in sheep: a molecular look at diversity-driving genes

BACKGROUND

After domestication, the evolution of phenotypically-varied sheep breeds has generated rich biodiversity. This wide phenotypic variation arises as a result of hidden genomic changes that range from a single nucleotide to several thousands of nucleotides. Thus, it is of interest and significance to reveal and understand the genomic changes underlying the phenotypic variation of sheep breeds in order to drive selection towards economically important traits.

REVIEW

Various traits contribute to the emergence of variation in sheep phenotypic characteristics, including coat color, horns, tail, wool, ears, udder, vertebrae, among others. The genes that determine most of these phenotypic traits have been investigated, which has generated knowledge regarding the genetic determinism of several agriculturally-relevant traits in sheep. In this review, we discuss the genomic knowledge that has emerged in the past few decades regarding the phenotypic traits in sheep, and our ultimate aim is to encourage its practical application in sheep breeding. In addition, in order to expand the current understanding of the sheep genome, we shed light on research gaps that require further investigation.

CONCLUSIONS

Although significant research efforts have been conducted in the past few decades, several aspects of the sheep genome remain unexplored. For the full utilization of the current knowledge of the sheep genome, a wide practical application is still required in order to boost sheep productive performance and contribute to the generation of improved sheep breeds. The accumulated knowledge on the sheep genome will help advance and strengthen sheep breeding programs to face future challenges in the sector, such as climate change, global human population growth, and the increasing demand for products of animal origin.

Transcriptional Characteristics Showed That miR-144-y/FOXO3 Participates in Embryonic Skin and Feather Follicle Development in Zhedong White Goose

Simple Summary

Feather is one of the most valuable and economical products in goose farming and plays a crucial physiological role in birds. For avian biology and the poultry industry, it is essential to comprehend and regulate how skin and feather follicles develop during embryogenesis. This study showed that several key regulatory genes (FOXO3, CTGF, and PTCH1, among others) and miRNAs (miR-144-y) participated in the developmental process of the skin and feather follicles in Zhedong white goose. Our findings are particularly important because they will serve as a valuable resource for upcoming studies on down feathers in agricultural economic growth regarding complex molecular mechanisms and breeding techniques.

Abstract

Skin and feather follicle development are essential processes for goose embryonic growth. Transcriptome and next-generation sequencing (NGS) network analyses were performed to improve the genome of Zhedong White goose and discover the critical genes, miRNAs, and pathways involved in goose skin and feather follicle morphogenesis. Sequencing output generated 6,002,591,668 to 8,675,720,319 clean reads from fifteen libraries. There were 1234, 3024, 4416, and 5326 different genes showing differential expression in four stages, E10 vs. E13, E10 vs. E18, E10 vs. E23, and E10 vs. E28, respectively. The differentially expressed genes (DEGs) were found to be implicated in multiple biological processes and pathways associated with feather growth and development, such as the Wnt signaling pathway, cell adhesion molecules, ECM–receptor interaction signaling pathways, and cell cycle and DNA replication pathways, according to functional analysis. In total, 8276 DEGs were assembled into twenty gene profiles with diverse expression patterns. The reliability of transcriptome results was verified by real-time quantitative PCR by selecting seven DEGs and five miRNAs. The localization of forkhead box O3 (FOXO3), connective tissue growth factor (CTGF), protein parched homolog1 (PTCH1), and miR-144-y by in situ hybridization showed spatial-temporal expression patterns and that FOXO3 and miR-144-y have an antagonistic targeting relationship. The correlation coefficient of FOXO3 and miR-144-y was -0.948, showing a strong negative correlation. Dual-luciferase reporter assay results demonstrated that miR-144-y could bind to the expected location to suppress the expression of FOXO3, which supports that there is a targeting relationship between them. The detections in this report will provide critical insight into the complex molecular mechanisms and breeding practices underlying the developmental characteristics of skin and feather follicles in Zhedong white geese.

Identification of photoperiod-induced specific miRNAs in the adrenal glands of Sunite sheep (Ovis aries)

In seasonal estrus, it is well known that melatonin-regulated biorhythm plays a key role. Some studies indicate that the adrenal gland plays an important role in reproduction in mammals, but the molecular mechanism is not clear. This study used an artificially controlled light photoperiod model, combined with RNA-seq technology and bioinformatics analysis, to analyze the messenger RNA (mRNA) and microRNA (miRNA) of ewe (Sunite) adrenal glands under different photoperiod treatments. After identification, the key candidate genes GRHL2, CENPF, FGF16 and SLC25A30 that photoperiod affects reproduction were confirmed. The miRNAs (oar-miR-544-3p, oar-miR-411b-5p, oar-miR-376e-3p, oar-miR-376d, oar-miR-376b-3p, oar-miR-376a-3p) were specifically expressed in the adrenal gland. The candidate mRNA-miRNA pairs (e.g., SLC25A30 coagulated by novel miRNA554, novel miRNA555 and novel miRNA559) may affect seasonal estrus. In summary, we constructed relation network of the mRNAs and miRNAs of sheep adrenal glands using RNA sequencing and bioinformatics analysis, thereby, providing a valuable genetic variation resource for sheep genome research, which will contribute to the study of complex traits in sheep.