miR-143 Targeting CUX1 to Regulate Proliferation of Dermal Papilla Cells in Hu Sheep

Chromosome-level genome sequencing and multi-omics of the Hungarian White Goose (Anser anser domesticus) reveals novel miRNA-mRNA regulation mechanism of waterfowl feather follicle development

The Hungarian White Goose (Anser anser domesticus) is an excellent European goose breed, with high feather and meat production. Despite its importance in the poultry industry, no available genome assembly information has been published. This study aimed to present Chromosome-level and functional genome sequencing of the Hungarian White Goose. The results showed that the genome assembly has a total length of 1115.82 Mb, 39 pairs of chromosomes, 92.98% of the BUSCO index, and contig N50 and scaffold N50 were up to 2.32 Mb and 60.69 Mb, respectively. Annotation of the genome assembly revealed 19550 genes, 286 miRNAs, etc. We identified 235 expanded and 1,167 contracted gene families in this breed compared with the other 16 species. We performed a positive selection analysis between this breed and four species of Anatidae to uncover the genetic information underlying feather follicle development. Further, we detected the function of miR-199-x, miR-143-y, and miR-23-z on goose embryonic skin fibroblast. In summary, we have successfully generated a highly complete genome sequence of the Hungarian white goose, which will provide a great resource to improve our understanding of gene functions and enhance the studies on feather follicle development at the genomic level.

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.

Association between DNA Methylation in the Core Promoter Region of the CUT-like Homeobox 1 (CUX1) Gene and Lambskin Pattern in Hu Sheep

CUT-like homeobox 1 (CUX1) has been proven to be a key regulator in sheep hair follicle development. In our previous study, CUX1 was identified as a differential expressed gene between Hu sheep lambskin with small wave patterns (SM) and straight wool patterns (ST); however, the exact molecular mechanism of CUX1 expression has been obscure. As DNA methylation can regulate the gene expression, the potential association between CUX1 core promotor region methylation and lambskin pattern in Hu sheep was explored in the present study. The results show that the core promoter region of CUX1 was present at (-1601-(-1) bp) upstream of the transcription start site. A repressive region (-1151-(-751) bp) was also detected, which had a strong inhibitory effect on CUX1 promoter activity. Bisulfite amplicon sequencing revealed that no significant difference was detected between the methylation levels of CUX1 core promoter region in SM tissues and ST tissues. Although the data demonstrated the differential expression of CUX1 between SM and ST probably has no association with DNA methylation, the identification of the core region and a potential repressive region of CUX1 promoter can enrich the role of CUX1 in Hu sheep hair follicle development.

Effect of CUX1 on the Proliferation of Hu Sheep Dermal Papilla Cells and on the Wnt/β-Catenin Signaling Pathway

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.

RNA-Seq Reveals the Roles of Long Non-Coding RNAs (lncRNAs) in Cashmere Fiber Production Performance of Cashmere Goats in China

Long non-coding RNAs (lncRNAs) are a kind of non-coding RNA being >200 nucleotides in length, and they are found to participate in hair follicle growth and development and wool fiber traits regulation. However, there are limited studies reporting the role of lncRNAs in cashmere fiber production in cashmere goats. In this study, Liaoning cashmere (LC) goats (n = 6) and Ziwuling black (ZB) goats (n = 6) with remarkable divergences in cashmere yield, cashmere fiber diameter, and cashmere color were selected for the construction of expression profiles of lncRNAs in skin tissue using RNA sequencing (RNA-seq). According to our previous report about the expression profiles of mRNAs originated from the same skin tissue as those used in the study, the cis and trans target genes of differentially expressed lncRNAs between the two caprine breeds were screened, resulting in a lncRNA-mRNA network. A total of 129 lncRNAs were differentially expressed in caprine skin tissue samples between LC goats and ZB goats. The presence of 2 cis target genes and 48 trans target genes for the differentially expressed lncRNAs resulted in 2 lncRNA-cis target gene pairs and 93 lncRNA-trans target gene pairs. The target genes concentrated on signaling pathways that were related to fiber follicle development, cashmere fiber diameter, and cashmere fiber color, including PPAR signaling pathway, metabolic pathways, fatty acid metabolism, fatty acid biosynthesis, tyrosine metabolism, and melanogenesis. A lncRNA-mRNA network revealed 22 lncRNA-trans target gene pairs for seven differentially expressed lncRNAs selected, of which 13 trans target genes contributed to regulation of cashmere fiber diameter, while nine trans target genes were responsible for cashmere fiber color. This study brings a clear explanation about the influences of lncRNAs over cashmere fiber traits in cashmere goats.

Discovery and Functional Analysis of Secondary Hair Follicle miRNAs during Annual Cashmere Growth

Secondary hair follicles (SHFs) produce the thermoregulatory cashmere of goats. MicroRNAs (miRNAs) play indispensable roles in hair follicle formation and growth. However, most studies examining miRNAs related to cashmere have been performed on goat skin. It remains unclear which miRNAs are highly expressed in SHFs or how miRNAs affect cashmere growth. In the present study, we isolated the SHFs under a dissecting microscope and analyzed the miRNA signatures during annual cashmere growth. Small-RNA sequencing followed by genome-wide expression analysis revealed that early anagen is a crucial phase for miRNA regulation of the cashmere growth, as revealed by two predominant groups of miRNAs. Although they exhibited opposite expression patterns, both groups demonstrated sharp changes of expression when in transit from early anagen to mid-anagen. In addition, we identified 96 miRNA signatures that were differentially expressed between different phases among 376 miRNAs. Functional analysis of the predicted target genes of highly expressed or differentially expressed miRNAs indicated that these miRNAs were involved in signal pathways associated with SHF development, regeneration, and regression. Furthermore, miR-143-3p was preferentially expressed in SHFs and Itga6 was identified as one of targets. The dual-luciferase and in situ hybridization assay demonstrated that miR-143-3p directly repressed the expression of Itga6, suggesting a possible novel role for miR-143-3p in cashmere growth.

Regulation of secondary hair follicle cycle in cashmere goats by miR-877-3p targeting IGFBP5 gene

Cashmere, a highly valuable animal product derived from cashmere goats, holds significant economic importance. MiRNAs serve as crucial regulators in the developmental processes of mammalian hair follicles. Understanding the regulation of miRNAs during the hair follicle cycle is essential for enhancing cashmere quality. In this investigation, we employed high-throughput sequencing technology to analyze the expression profiles of miRNAs in the secondary hair follicles of Jiangnan cashmere goats at different stages. Through bioinformatics analysis, we identified differentially expressed miRNAs (DE miRNAs). The regulatory relationships between miRNAs and their target genes were verified using multiple techniques, including RT-qPCR, western blot, Dual-Luciferase Reporter, and CKK-8 assays. Our findings revealed the presence of 193 DE miRNAs during various stages of the hair follicle cycle in Jiangnan cashmere goats. Based on the previously obtained mRNA data, the target genes of DE miRNA were predicted, and 1,472 negative regulatory relationships between DE miRNAs and target genes were obtained. Notably, the expression of chi-miR-877-3p was down-regulated during the telogen (Tn) phase compared to the anagen (An) and catagen (Cn) phases, while the IGFBP5 gene exhibited up-regulation. Further validation experiments confirmed that overexpression of chi-miR-877-3p in dermal papilla cells suppressed IGFBP5 gene expression and facilitated cell proliferation. The results of this study provide novel insights for analyzing the hair follicle cycle.

Effect of Sox18 on the Induction Ability of Dermal Papilla Cells in Hu Sheep

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

Comparative Analysis of mRNA and miRNA Expression between Dermal Papilla Cells and Hair Matrix Cells of Hair Follicles in Yak

The interaction between the dermal papilla cells (DPCs) and epidermal hair matrix cells (HMCs) of hair follicles (HFs) is crucial for the growth and development of HFs, but the molecular mechanism is complex and remains unclear. MicroRNAs (miRNAs) are the key signaling molecules for cellular communication. In this study, the DPCs and HMCs of yak were isolated and cultured, and the differentially expressed mRNA and miRNA were characterized to analyze the molecular basis of the interaction between DPCs and HMCs during hair follicle (HF) development in yak. The mRNA differential expression and functional enrichment analysis revealed that there were significant differences between DPCs and HMCs, and they showed the molecular functional characteristics of dermal cells and epidermal cells, respectively. Multiple KEGG pathways related to HF development were enriched in the highly expressed genes in DPCs, while the pathways associated with microbiota and immunity were significantly enriched in the highly expressed genes in HMCs. By combining analysis with our previous 10× genomics single-cell transcriptome data, 39 marker genes of DPCs of yak were identified. A total of 123 relatively specifically expressed miRNAs were screened; among these, the miRNAs associated with HF development such as miR-143, miR-214, miR-125b, miR-31, and miR-200 were presented. In conclusion, the large changes in yak DPCs and HMCs for both mRNA and miRNA expression were revealed, and numerous specifically expressed mRNAs and miRNAs in DPCs or HMCs were identified, which may contribute to the interaction and cellular communication between DPCs and HMCs during HF development in yak.

An integrative analysis of the lncRNA-miRNA-mRNA competitive endogenous RNA network reveals potential mechanisms in the murine hair follicle cycle

Alopecia is a common progressive disorder associated with abnormalities of the hair follicle cycle. Hair follicles undergo cyclic phases of hair growth (anagen), regression (catagen), and rest (telogen), which are precisely regulated by various mechanisms. However, the specific mechanism associated with hair follicle cycling, which includes noncoding RNAs and regulation of competitive endogenous RNA (ceRNA) network, is still unclear. We obtained data from publicly available databases and performed real-time quantitative polymerase chain reaction validations. These analyses revealed an increase in the expression of miRNAs and a decrease in the expression of target mRNAs and lncRNAs from the anagen to telogen phase of the murine hair follicle cycle. Subsequently, we constructed the ceRNA networks and investigated their functions using enrichment analysis. Furthermore, the androgenetic alopecia (AGA) microarray data analysis revealed that several novel alopecia-related genes were identified in the ceRNA networks. Lastly, GSPT1 expression was detected using immunohistochemistry. Our analysis revealed 11 miRNAs (miR-148a-3p, miR-146a-5p, miR-200a-3p, miR-30e-5p, miR-30a-5p, miR-27a-3p, miR-143-3p, miR-27b-3p, miR-126a-3p, miR-378a-3p, and miR-22-3p), 9 target mRNAs (Atp6v1a, Cdkn1a, Gadd45a, Gspt1, Mafb, Mitf, Notch1, Plk2, and Slc7a5), and 2 target lncRNAs (Neat1 and Tug1) were differentially expressed in hair follicle cycling. The ceRNA networks were made of 12 interactive miRNA-mRNA pairs and 13 miRNA-lncRNA pairs. The functional enrichment analysis revealed the enrichment of hair growth–related signaling pathways. Additionally, GSPT1 was downregulated in androgenetic alopecia patients, possibly associated with alopecia progression. The ceRNA network identified by our analysis could be involved in regulating the hair follicle cycle.

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.