Time-course RNA-seq analysis reveals stage-specific and melatonin-triggered gene expression patterns during the hair follicle growth cycle in Capra hircus

Transcriptome meta-analysis reveals the hair genetic rules in six animal breeds and genes associated with wool fineness

Wool plays an irreplaceable role in the lives of livestock and the textile industry. The variety of hair quality and shape leads to the diversity of its functions and applications, and the finer wool has a higher economic value. In this study, 10 coarse and 10 fine ordos fine wool sheep skin samples were collected for RNA-seq, and coarse and fine skin/hair follicle RNA-seq datasets of other five animal breeds were obtained from NCBI. Weighted gene co-expression network analysis showed that the common genes were clustered into eight modules. Similar gene expression patterns in sheep and rabbits with the same wool types, different gene expression patterns in animal species with different hair types, and brown modules were significantly correlated with species and breeds. GO and KEGG enrichment analyses showed that, most genes in the brown module associated with hair follicle development. Hence, gene expression patterns in skin tissues may determine hair morphology in animal. The analysis of differentially expressed genes revealed that 32 highly expressed candidate genes associated with the wool fineness of Ordos fine wool sheep. Among them, KAZALD1 (grey module), MYOC (brown module), C1QTNF6 (brown module), FOS (tan module), ITGAM, MX2, MX1, and IFI6 genes have been reported to be involved in the regulation of the hair follicle cycle or hair loss. Additionally, 12 genes, including KAZALD1, MYOC, C1QTNF6, and FOS, are differentially expressed across various animal breeds and species. The above results suggest that different sheep breeds share a similar molecular regulatory basis of wool fineness. Finally, the study provides a theoretical reference for molecular breeding of sheep breeds as well as for the investigation of the origin and evolution of animal hair.

Genomic Inbreeding and Runs of Homozygosity Analysis of Cashmere Goat

Cashmere goats are valuable genetic resources which are famous worldwide for their high-quality fiber. Runs of homozygosity (ROHs) have been identified as an efficient tool to assess inbreeding level and identify related genes under selection. However, there is limited research on ROHs in cashmere goats. Therefore, we investigated the ROH pattern, assessed genomic inbreeding levels and examined the candidate genes associated with the cashmere trait using whole-genome resequencing data from 123 goats. Herein, the Inner Mongolia cashmere goat presented the lowest inbreeding coefficient of 0.0263. In total, we identified 57,224 ROHs. Seventy-four ROH islands containing 50 genes were detected. Certain identified genes were related to meat, fiber and milk production (FGF1, PTPRM, RERE, GRID2, RARA); fertility (BIRC6, ECE2, CDH23, PAK1); disease or cold resistance and adaptability (PDCD1LG2, SVIL, PRDM16, RFX4, SH3BP2); and body size and growth (TMEM63C, SYN3, SDC1, STRBP, SMG6). 135 consensus ROHs were identified, and we found candidate genes (FGF5, DVL3, NRAS, KIT) were associated with fiber length or color. These findings enhance our comprehension of inbreeding levels in cashmere goats and the genetic foundations of traits influenced by selective breeding. This research contributes significantly to the future breeding, reservation and use of cashmere goats and other goat breeds.

Transcriptome Profiling Unveils Key Genes Regulating the Growth and Development of Yangzhou Goose Knob

Goose is one of the most economically valuable poultry species and has a distinct appearance due to its possession of a knob. A knob is a hallmark of sexual maturity in goose (Anser cygnoides) and plays crucial roles in artificial selection, health status, social signaling, and body temperature regulation. However, the genetic mechanisms influencing the growth and development of goose knobs remain completely unclear. In this study, histomorphological and transcriptomic analyses of goose knobs in D70, D120, and D300 Yangzhou geese revealed differential changes in tissue morphology during the growth and development of goose knobs and the key core genes that regulate goose knob traits. Observation of tissue sections revealed that as age increased, the thickness of the knob epidermis, cuticle, and spinous cells gradually decreased. Additionally, fat cells in the dermis and subcutaneous connective tissue transitioned from loose to dense. Transcriptome sequencing results, analyzed through differential expression, Weighted Gene Co-expression Network Analysis (WGCNA), and pattern expression analysis methods, showed D70-vs.-D120 (up-regulated: 192; down-regulated: 423), D70-vs.-D300 (up-regulated: 1394; down-regulated: 1893), and D120-vs.-D300 (up-regulated: 1017; down-regulated: 1324). A total of 6243 differentially expressed genes (DEGs) were identified, indicating varied expression levels across the three groups in the knob tissues of D70, D120, and D300 Yangzhou geese. These DEGs are significantly enriched in biological processes (BP) such as skin morphogenesis, the regulation of keratinocyte proliferation, and epidermal cell differentiation. Furthermore, they demonstrate enrichment in pathways related to goose knob development, including ECM-receptor interaction, NF-kappa B, and PPAR signaling. Through pattern expression analysis, three gene expression clusters related to goose knob traits were identified. The joint analysis of candidate genes associated with goose knob development and WGCNA led to the identification of key core genes influencing goose knob development. These core genes comprise WNT4, WNT10A, TCF7L2, GATA3, ADRA2A, CASP3, SFN, KDF1, ERRFI1, SPRY1, and EVPL. In summary, this study provides a reference for understanding the molecular mechanisms of goose knob growth and development and provides effective ideas and methods for the genetic improvement of goose knob traits.

Light: An Extrinsic Factor Influencing Animal-based Research

Light is an environmental factor that is extrinsic to animals themselves and that exerts a profound influence on the regulation of circadian, neurohormonal, metabolic, and neurobehavioral systems of all animals, including research animals. These widespread biologic effects of light are mediated by distinct photoreceptors-rods and cones that comprise the conventional visual system and melanopsin-containing intrinsically photosensitive retinal ganglion cells (ipRGCs) of the nonvisual system that interact with the rods and cones. The rods and cones of the visual system, along with the ipRGCs of the nonvisual system, are species distinct in terms of opsins and opsin concentrations and interact with one another to provide vision and regulate circadian rhythms of neurohormonal and neurobehavioral responses to light. Here, we review a brief history of lighting technologies, the nature of light and circadian rhythms, our present understanding of mammalian photoreception, and current industry practices and standards. We also consider the implications of light for vivarium measurement, production, and technological application and provide simple recommendations on artificial lighting for use by regulatory authorities, lighting manufacturers, designers, engineers, researchers, and research animal care staff that ensure best practices for optimizing animal health and well-being and, ultimately, improving scientific outcomes.

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.

Identification of N6-Methyladenosine-Related Factors and the Prediction of the Regulatory Mechanism of Hair Follicle Development in Rex and Hycole Rabbits

Hair follicle development directly affects the development of the rabbit fur industry. The growth and development of a hair follicle is modified and regulated by many genes and mechanisms. M6A is an important RNA modification. However, there are few studies on the effects of the regulation of m6A on hair follicle growth and development. In this study, hematoxylin-eosin (HE) staining was used to explore the difference in hair follicle development between Rex rabbits and Hycole rabbits, and we performed m6A sequencing to identify the key genes with m6A modification in hair follicle growth. The results showed that the hair length, coarse hair percentage, primary hair follicle ratio, and skin thickness of Hycole rabbits were significantly higher than those of Rex rabbits. However, the proportion of secondary hair follicles in Hycole rabbits was significantly lower than that in Rex rabbits. In addition, we found five differential methylases, 20 differential genes, and 24 differential signaling pathways related to hair growth and development. The results of the Sankey diagram showed that 12 genes were related to 13 signal pathways. Finally, we found that five methylases regulated the development of hair follicles through differential genes/signal pathways. These findings laid a molecular foundation for the function of m6A modification in hair development.

Hedgehog signaling is required for endometrial remodeling and myometrial homeostasis in the cycling mouse uterus

Decades of work demonstrate that the mammalian estrous cycle is controlled by cycling steroid hormones. However, the signaling mechanisms that act downstream, linking hormonal action to the physical remodeling of the cycling uterus, remain unclear. To address this issue, we analyzed gene expression at all stages of the mouse estrous cycle. Strikingly, we found that several genetic programs well-known to control tissue morphogenesis in developing embryos displayed cyclical patterns of expression. We find that most of the genetic architectures of Hedgehog signaling (ligands, receptors, effectors, and transcription factors) are transcribed cyclically in the uterus, and that conditional disruption of the Hedgehog receptor smoothened not only elicits a failure of normal cyclical thickening of the endometrial lining but also induces aberrant deformation of the uterine smooth muscle. Together, our data shed light on the mechanisms underlying normal uterine remodeling specifically and cyclical gene expression generally.

LncRNAs regulate the cyclic growth and development of hair follicles in Dorper sheep

Introduction

Hair follicles in Dorper sheep are characterized by seasonal cyclic growth and development, consequently resulting in hair shedding during spring. The cyclic growth and development of hair follicles are regulated by several influencing factors such as photoperiods, hormones, age of the animal, genes, long non-coding RNAs (lncRNAs), and signaling pathways.

Methods

In the present study, skin samples of five shedding sheep (S), used as experimental animals, and three non-shedding sheep (N), used as controls, were collected at three time points (September 27, 2019; January 3, 2020; and March 17, 2020) for RNA sequencing (RNA-seq) technology. Nine different groups (S1-vs-S2, S1-vs-S3, S2-vs-S3, N1- vs-N2, N1-vs-N3, N2-vs-N3, S1-vs-N1, S2-vs-N2, and S3-vs-N3) were compared using FDR < 0.05 and log 21 FC >as thresholds to assess the differences in the expression of lncRNAs.

Results and discussion

In total, 395 differentially expressed (DE) lncRNAs were screened. Cluster heatmap analysis identified two types of expression patterns, namely, high expression during the anagen phase (A pattern) and high expression during the telogen phase (T pattern). Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the target genes were largely enriched in the Estrogen signaling pathway, PI3K-Akt signaling pathway, Fc gamma R-mediated phagocytosis, and cell adhesion molecules (CAMs), which are associated with hair follicle cyclic growth and development-related pathways. In addition, 17 pairs of lncRNAs-target genes related to hair follicle cyclic growth and development were screened, and a regulatory network was constructed. Altogether, candidate lncRNAs and their regulated target genes were screened that contributed to sheep hair follicle cyclic growth and development. We believe these findings will provide useful insights into the underlying regulatory mechanisms.

Altered hair root gene expression profiles highlight calcium signaling and lipid metabolism pathways to be associated with curly hair initiation and maintenance in Mangalitza pigs

Hair types have been under strong targeted selection in domestic animals for their impact on skin protection, thermoregulation and exterior morphology, and subsequent economic importance. In pigs, a very special hair phenotype was observed in Mangalitza, who expresses a thick coat of curly bristles and downy hair. Two breed-specific missense variants in TRPM2 and CYP4F3 were suggested to be associated with the Mangalitza pig's hair shape due to their role in hair follicle morphogenesis reported for human and mice. However, the mechanism behind this expression of a curly hair type is still unclear and needs to be explored. In our study, hair shafts were measured and investigated for the curvature of the hair in Mangalitza and crossbreeds in comparison to straight-coated pigs. For molecular studies, hair roots underwent RNA sequencing for a differential gene expression analysis using DESeq2. The output matrix of normalized counts was then used to construct weighted gene co-expression networks. The resulting hair root gene expression profiles highlighted 454 genes to be significantly differentially expressed for initiation of curly hair phenotype in newborn Mangalitza piglets versus post-initiation in later development. Furthermore, 2,554 genes showed a significant differential gene expression in curly hair in comparison to straight hair. Neither TRPM2 nor CYP4F3 were identified as differentially expressed. Incidence of the genes in weighted co-expression networks associated with TRPM2 and CYP4F3, and prominent interactions of subsequent proteins with lipids and calcium-related pathways suggested calcium signaling and/or lipid metabolism as essential players in the induction of the curly hair as well as an ionic calcium-dependency to be a prominent factor for the maintenance of this phenotype. Subsequently, our study highlights the complex interrelations and dependencies of mutant genes TRPM2 and CYP4F3 and associated gene expression patterns, allowing the initiation of curly hair type during the development of a piglet as well as the maintenance in adult individuals.

Hair Follicle Transcriptome Analysis Reveals Differentially Expressed Genes That Regulate Wool Fiber Diameter in Angora Rabbits

Wool fiber diameter (WFD) is an important index of wool traits and the main determinant of wool quality and value. However, the genetic determinants of fiber diameter have not yet been fully elucidated. Here, coarse and fine wool of Wan strain Angora rabbits and their hair follicle traits were characterized. The results indicated significant differences in the diameters of wool fibers and their hair follicles. The RNA sequencing (RNA-Seq) technique was used to identify differences in gene expression in hair follicles between coarse and fine wool. In total, 2574 differentially expressed genes (DEGs) were found between the two hair follicle groups. Transcription factors, keratin-associated protein (KAP) and keratin (KRT) families, and ECM-related genes may control the structure of fine fibers in rabbits. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that skin development, epidermal cell and keratinocyte differentiation, epithelium development, and Notch and ribosome signaling pathways were significantly enriched, respectively. GSEA further filtered six important pathways and related core genes. PPI analysis also mined functional DEGs associated with hair structure, including LEF1, FZD3, SMAD3, ITGB6, and BMP4. Our findings provide valuable information for researching the molecular mechanisms regulating wool fiber and could facilitate enhanced selection of super-fine wool rabbits through gene-assisted selection in the future.

Melatonin Promotes the Development of Secondary Hair Follicles in Adult Cashmere Goats by Activating the Keap1-Nrf2 Signaling Pathway and Inhibiting the Inflammatory Transcription Factors NFκB and AP-1

Exogenous melatonin (MT) has been used to promote the growth of secondary hair follicles and improve cashmere fiber quality, but the specific cellular-level mechanisms involved are unclear. This study was carried out to investigate the effect of MT on the development of secondary hair follicles and on cashmere fiber quality in cashmere goats. The results showed that MT improved secondary follicle numbers and function as well as enhanced cashmere fiber quality and yield. The MT-treated goat groups had high secondary-to-primary ratios (S:P) for hair follicles, greater in the elderly group (p < 0.05). Antioxidant capacities of secondary hair follicles improved fiber quality and yield in comparison with control groups (p < 0.05/0.01). Levels of reactive oxygen and nitrogen species (ROS, RNS) and malondialdehyde (MDA) were lowered (p < 0.05/0.01) by MT. There was significant upregulation of antioxidant genes (for SOD-3; GPX-1; NFE2L2) and the protein of nuclear factor (Nrf2), and downregulation of the Keap1 protein. There were significant differences in the expression of genes for secretory senescence-associated phenotype (SASP) cytokines (IL-1β, IL-6, MMP-9, MMP-27, CCL-21, CXCL-12, CXCL-14, TIMP-1,2,3) plus their protein of key transcription factors, nuclear factor kappa B (NFκB) and activator protein-1 (AP-1), in comparison with the controls. We concluded that MT could enhance antioxidant capacity and reduce ROS and RNS levels of secondary hair follicles through the Keap1-Nrf2 pathway in adult cashmere goats. Furthermore, MT reduced the expression of the SASP cytokines genes by inhibiting the protein of NFκB and AP-1 in the secondary hair follicles in older cashmere goats, thus delaying skin aging, improving follicle survival, and increasing the number of secondary hair follicles. Collectively, these effects of exogenous MT enhanced the quality and yield of cashmere fibers, especially at 5-7 years old.

Construction of transcriptome atlas of white yak hair follicle during anagen and catagen using single-cell RNA sequencing

BACKGROUND

As the direct organ of villus, hair follicles have obvious seasonal cycles. The hair follicle cycle is orchestrated by multiple cell types that together direct cell renewal and differentiation. But the regulation property of hair follicle cells from anagen to catagen in yak is still unknown.

RESULTS

In this study, single-cell RNA sequencing was performed on 24,124 single cells of the scapular skin from white yak. Based on tSNE cluster analysis, the cell types of IFE-DC, epidermal cell lines, fibroblasts, keratinocytes, IRS, DS, INFU, and other cells in yak hair follicles during anagen and catagen were successfully identified, and the gene expression profiles were described. The GO enrichment analysis indicated the different cells characteristic genes to be mainly enriched in the epidermal development, epithelial cell differentiation and wound healing pathways. The pseudotime trajectory analysis described the differentiation trajectory of the epidermal lineage and dermal lineage of the hair follicle during anagen and catagen. Moreover, the dynamic changes of the genes like LHX2, KRT25, and KRT71 were found to be highly expressed in HS and IRS, but not in the IFE-DC, INFU, and keratinocyte during differentiation.

CONCLUSIONS

Our results analyzed the time-varying process of gene expression in the dermal cell lineage and epidermal cell lineage of hair follicles during anagen and catagen during fate differentiation was expounded at the single cell level, revealing the law of fate specialization of different types of cells. In addition, based on the enrichment analysis, the transcriptional regulatory factors involved in the different cell fates were also revealed. These results will help to enhance our understanding of yak hair follicle cycle and promote the development and utilization of yak villus.