Single-cell Transcriptome Profiling reveals Dermal and Epithelial cell fate decisions during Embryonic Hair Follicle Development

Integrated-omics profiling unveils the disparities of host defense to ECM scaffolds during wound healing in aged individuals

Extracellular matrix (ECM) scaffold membranes have exhibited promising potential to better the outcomes of wound healing by creating a regenerative microenvironment around. However, when compared to the application in younger individuals, the performance of the same scaffold membrane in promoting re-epithelialization and collagen deposition was observed dissatisfying in aged mice. To comprehensively explore the mechanisms underlying this age-related disparity, we conducted the integrated analysis, combing single-cell RNA sequencing (scRNA-Seq) with spatial transcriptomics, and elucidated six functionally and spatially distinctive macrophage groups and lymphocytes surrounding the ECM scaffolds. Through intergroup comparative analysis and cell-cell communication, we characterized the dysfunction of Spp1+ macrophages in aged mice impeded the activation of the type Ⅱ immune response, thus inhibiting the repair ability of epidermal cells and fibroblasts around the ECM scaffolds. These findings contribute to a deeper understanding of biomaterial applications in varied physiological contexts, thereby paving the way for the development of precision-based biomaterials tailored specifically for aged individuals in future therapeutic strategies.

Application of multi-omics techniques to androgenetic alopecia: Current status and perspectives

The rapid advancement of sequencing technologies has enabled the generation of vast datasets, allowing for the in-depth analysis of sequencing data. This analysis has facilitated the validation of novel pathogenesis hypotheses for understanding and treating diseases through ex vivo and in vivo experiments. Androgenetic alopecia (AGA), a common hair loss disorder, has been a key focus of investigators attempting to uncover its underlying mechanisms. Abnormal changes in mRNA, proteins, and metabolites have been identified in individuals with AGA, and future developments in sequencing technologies may reveal new biomarkers for AGA. By integrating multiple omics analysis datasets such as genomics, transcriptomics, proteomics, and metabolomics-along with clinical phenotype data-we can achieve a comprehensive understanding of the molecular underpinnings of AGA. This review summarizes the data-mining studies conducted on various omics analysis datasets as related to AGA that have been adopted to interpret the biological data obtained from different omics layers. We herein discuss the challenges of integrative omics analyses, and suggest that collaborative multi-omics studies can enhance the understanding of the complete pathomechanism(s) of AGA by focusing on the interaction networks comprising DNA, RNA, proteins, and metabolites.

Unveiling microglia heterogeneity in intracerebral hemorrhage

Microglia are important innate immune cells in the brain, and a rich diversity of subtypes has recently been discovered that expand beyond the traditional classification of traditional M1 (pro-inflammatory) and M2 (anti-inflammatory) classifications. Intracerebral hemorrhage (ICH) is a devastating form of stroke, and the understanding of its later-stage pathological mechanisms remains incomplete. In this study, through the analysis of single-cell transcripts from mice brains 14 days post-ICH, three disease-associated expression patterns of microglia were identified. These include a lipid metabolism and phagocytosis phenotype reminiscent of Disease-Associated Microglia (DAM) initially discovered in Alzheimer's disease models, a phenotype associated with angiogenesis, and a relatively independent phenotype similar to the pro-inflammatory M1 state. These findings were further validated through immunofluorescence in both mouse and human specimens. In addition, analysis of single-cell transcripts from mice brains 3 days post-ICH suggested that microglia involved in lipid metabolism and phagocytosis likely emerge from early proliferating populations. Given the distinct origins and phenotypic characteristics of pro-inflammatory and reparative microglia, interventions targeting these cells hold the potential to modulate the delicate balance between injury and repair during the pathophysiological process of ICH, highlighting a pivotal direction for future therapeutic strategies.

GPC3-targeted CAR-T cells expressing GLUT1 or AGK exhibit enhanced antitumor activity against hepatocellular carcinoma

Chimeric antigen receptor-expressing T (CAR-T) cells induce robust antitumor responses in patients with hematologic malignancies. However, CAR-T cells exhibit only limited efficacy against solid tumors such as hepatocellular carcinoma (HCC), partially due to their limited expansion and persistence. CD8+ T cells, as key components of the adaptive immune response, play a central role in antitumor immunity. Aerobic glycolysis is the main metabolic feature of activated CD8+ T cells. In the tumor microenvironment, however, the uptake of large amounts of glucose by tumor cells and other immunosuppressive cells can impair the activation of T cells. Only when tumor-infiltrating lymphocytes (TILs) in the tumor microenvironment have a glycolytic advantage might the effector function of T cells be activated. Glucose transporter type 1 (GLUT1) and acylglycerol kinase (AGK) can boost glycolytic metabolism and activate the effector function of CD8+ T cells, respectively. In this study, we generated GPC3-targeted CAR-T cells overexpressing GLUT1 or AGK for the treatment of HCC. GPC3-targeted CAR-T cells overexpressing GLUT1 or AGK specifically and effectively lysed GPC3-positive tumor cells in vitro in an antigen-dependent manner. Furthermore, GLUT1 or AGK overexpression protected CAR-T cells from apoptosis during repeated exposures to tumor cells. Compared with second-generation CAR-T cells, GPC3-targeted CAR-T cells overexpressing GLUT1 or AGK exhibited greater CD8+ T-cell persistence in vivo and better antitumor effects in HCC allograft mouse models. Finally, we revealed that GLUT1 or AGK maintained anti-apoptosis ability in CD8+ T cells via activation of the PI3K/Akt pathway. This finding might identify a therapeutic strategy for advanced HCC.

TLR9 activation in large wound induces tissue repair and hair follicle regeneration via γδT cells

The mechanisms underlying tissue repair in response to damage have been one of main subjects of investigation. Here we leverage the wound-induced hair neogenesis (WIHN) models in adult mice to explore the correlation between degree of damage and the healing process and outcome. The multimodal analysis, in combination with single-cell RNA sequencing help to explore the difference in wounds of gentle and heavy damage degrees, identifying the potential role of toll-like receptor 9 (TLR9) in sensing the injury and regulating the immune reaction by promoting the migration of γδT cells. The TLR9 deficient mice or wounds injected with TLR9 antagonist have greatly impaired healing and lower WIHN levels. Inhibiting the migration of γδT cells or knockout of γδT cells also suppress the wound healing and regeneration, which can't be rescued by TLR9agonist. Finally, the amphiregulin (AREG) is shown as one of most important effectors secreted by γδT cells and keratinocytes both in silicon or in the laboratory, whose expression influences WIHN levels and the expression of stem cell markers. In total, our findings reveal a previously unrecognized role for TLR9 in sensing skin injury and influencing the tissue repair and regeneration by modulation of the migration of γδT cells, and identify the TLR9-γδT cells-areg axis as new potential targets for enhancing tissue regeneration.

Single-Cell and Spatial Transcriptome Analysis of Dermal Fibroblast Development in Perinatal Mouse Skin: Dynamic Lineage Differentiation and Key Driver Genes

Several single-cell RNA studies of developing mouse skin have elucidated the molecular and cellular processes involved in skin development. However, they have primarily focused on either the fetal or early postnatal period, leaving a gap in our understanding of skin development. In this study, we conducted a comprehensive time-series analysis by combining single-cell RNA-sequencing datasets collected at different stages of development (embryonic days 13.5, 14.5, and 16.5 and postnatal days 0, 2, and 4) and validated our findings through multipanel in situ spatial transcriptomics. Our analysis indicated that embryonic fibroblasts exhibit heterogeneity from a very early stage and that the rapid determination of each lineage occurs within days after birth. The expression of putative key driver genes, including Hey1, Ebf1, Runx3, and Sox11 for the dermal papilla trajectory; Lrrc15 for the dermal sheath trajectory; Zfp536 and Nrn1 for the papillary fibroblast trajectory; and Lrrn4cl and Mfap5 for the fascia fibroblast trajectory, was detected in the corresponding, spatially identified cell types. Finally, cell-to-cell interaction analysis indicated that the dermal papilla lineage is the primary source of the noncanonical Wnt pathway during skin development. Together, our study provides a transcriptomic reference for future research in the field of skin development and regeneration.

Single-cell transcriptome profiling reveals cell type-specific variation and development in HLA expression of human skin

Skin, the largest organ of body, is a highly immunogenic tissue with a diverse collection of immune cells. Highly polymorphic human leukocyte antigen (HLA) molecules have a central role in coordinating immune responses as recognition molecules. Nevertheless, HLA gene expression patterns among diverse cell types within a specific organ, like the skin, have yet to be thoroughly investigated, with stromal cells attracting much less attention than immune cells. To illustrate HLA expression profiles across different cell types in the skin, we performed single-cell RNA sequencing (scRNA-seq) analyses on skin datasets, covering adult and fetal skin, and hair follicles as the skin appendages. We revealed the variation in HLA expression between different skin populations by examining normal adult skin datasets. Moreover, we evaluated the potential immunogenicity of multiple skin populations based on the expression of classical HLA class I genes, which were well represented in all cell types. Furthermore, we generated scRNA-seq data of developing skin from fetuses of 15 post conception weeks (PCW), 17 PCW, and 22 PCW, delineating the dynamic expression of HLA genes with cell type-dependent variation among various cell types during development. Notably, the pseudotime trajectory analysis unraveled the significant variance in HLA genes during the evolution of vascular endothelial cells. Moreover, we uncovered the immune-privileged properties of hair follicles at single-cell resolution. Our study presents a comprehensive single-cell transcriptomic landscape of HLA genes in the skin, which provides new insights into variation in HLA molecules and offers a clue for allogeneic skin transplantation.

Spatially resolved transcriptomic profiling of placental development in dairy cow

The placenta plays a crucial role in successful mammalian reproduction. Ruminant animals possess a semi-invasive placenta characterized by a highly vascularized structure formed by maternal endometrial caruncles and fetal placental cotyledons, essential for full-term fetal development. The cow placenta harbors at least two trophoblast cell populations: uninucleate (UNC) and binucleate (BNC) cells. However, the limited capacity to elucidate the transcriptomic dynamics of the placental natural environment has resulted in a poor understanding of both the molecular and cellular interactions between trophoblast cells and niches, and the molecular mechanisms governing trophoblast differentiation and functionalization. To fill this knowledge gap, we employed Stereo-seq to map spatial gene expression patterns at near single-cell resolution in the cow placenta at 90 and 130 days of gestation, attaining high-resolution, spatially resolved gene expression profiles. Based on clustering and cell marker gene expression analyses, key transcription factors, including YBX1 and NPAS2, were shown to regulate the heterogeneity of trophoblast cell subpopulations. Cell communication and trajectory analysis provided a framework for understanding cell-cell interactions and the differentiation of trophoblasts into BNCs in the placental microenvironment. Differential analysis of cell trajectories identified a set of genes involved in regulation of trophoblast differentiation. Additionally, spatial modules and co-variant genes that help shape specific tissue structures were identified. Together, these findings provide foundational insights into important biological pathways critical to the placental development and function in cows.

Integrating Single-Cell and Spatial Transcriptomics Reveals Heterogeneity of Early Pig Skin Development and a Subpopulation with Hair Placode Formation

The dermis and epidermis, crucial structural layers of the skin, encompass appendages, hair follicles (HFs), and intricate cellular heterogeneity. However, an integrated spatiotemporal transcriptomic atlas of embryonic skin has not yet been described and would be invaluable for studying skin-related diseases in humans. Here, single-cell and spatial transcriptomic analyses are performed on skin samples of normal and hairless fetal pigs across four developmental periods. The cross-species comparison of skin cells illustrated that the pig epidermis is more representative of the human epidermis than mice epidermis. Moreover, Phenome-wide association study analysis revealed that the conserved genes between pigs and humans are strongly associated with human skin-related diseases. In the epidermis, two lineage differentiation trajectories describe hair follicle (HF) morphogenesis and epidermal development. By comparing normal and hairless fetal pigs, it is found that the hair placode (Pc), the most characteristic initial structure in HFs, arises from progenitor-like OGN+/UCHL1+ cells. These progenitors appear earlier in development than the previously described early Pc cells and exhibit abnormal proliferation and migration during differentiation in hairless pigs. The study provides a valuable resource for in-depth insights into HF development, which may serve as a key reference atlas for studying human skin disease etiology using porcine models.

Gynostemma pentaphyllum Hydrodistillate and Its Major Component Damulin B Promote Hair Growth-Inducing Properties In Vivo and In Vitro via the Wnt/β-Catenin Pathway in Dermal Papilla Cells

Alopecia, a prevalent yet challenging condition with limited FDA-approved treatments which is accompanied by notable side effects, necessitates the exploration of natural alternatives. This study elucidated the hair growth properties of Gynostemma pentaphyllum leaf hydrodistillate (GPHD) both in vitro and in vivo. Furthermore, damulin B, a major component of GPHD, demonstrated hair growth-promoting properties in vitro. Beyond its established anti-diabetic, anti-obesity, and anti-inflammatory attributes, GPHD exhibited hair growth induction in mice parallel to minoxidil. Moreover, it upregulated the expression of autocrine factors associated with hair growth, including VEGF, IGF-1, KGF, and HGF. Biochemical assays revealed that minoxidil, GPHD, and damulin B induced hair growth via the Wnt/β-catenin pathway through AKT signaling, aligning with in vivo experiments demonstrating improved expression of growth factors. These findings suggest that GPHD and damulin B contribute to the hair growth-inducing properties of dermal papilla cells through the AKT/β-catenin signaling pathway.

MDIC3: Matrix decomposition to infer cell-cell communication

Crosstalk among cells is vital for maintaining the biological function and intactness of systems. Most existing methods for investigating cell-cell communications are based on ligand-receptor (L-R) expression, and they focus on the study between two cells. Thus, the final communication inference results are particularly sensitive to the completeness and accuracy of the prior biological knowledge. Because existing L-R research focuses mainly on humans, most existing methods can only examine cell-cell communication for humans. As far as we know, there is currently no effective method to overcome this species limitation. Here, we propose MDIC3 (matrix decomposition to infer cell-cell communication), an unsupervised tool to investigate cell-cell communication in any species, and the results are not limited by specific L-R pairs or signaling pathways. By comparing it with existing methods for the inference of cell-cell communication, MDIC3 obtained better performance in both humans and mice.

Single-Cell Transcriptome Sequence Profiling on the Morphogenesis of Secondary Hair Follicles in Ordos Fine-Wool Sheep

The Ordos fine-wool sheep is a high-quality breed in China that produces superior natural textiles and raw materials such as wool and lamb meat. However, compared to the Australian Merino sheep, there is still a gap in terms of the wool fiber fineness and wool yield. The hair follicle is the main organ that controls the type of wool fiber, and the morphological changes in the secondary hair follicle are crucial in determining wool quality. However, the process and molecular mechanisms of hair follicle morphogenesis in Ordos fine-wool sheep are not yet clear. Therefore, analyzing the molecular mechanisms underlying the process of follicle formation is of great significance for improving the fiber diameter and wool production of Ordos fine-wool sheep. The differential expressed genes, APOD, POSTN, KRT5, and KRT15, which related to primary hair follicles and secondary hair follicles, were extracted from the dermal papillae. Based on pseudo-time analysis, the differentiation trajectories of dermal lineage cells and epidermal lineage cells in the Ordos fine-wool sheep were successfully constructed, providing a theoretical basis for breeding research in Ordos fine-wool sheep.

Integrating scRNA-seq and bulk RNA-seq to characterize infiltrating cells in the colorectal cancer tumor microenvironment and construct molecular risk models

Colorectal cancer (CRC) is a malignancy that is both highly lethal and heterogeneous. Although the correlation between intra-tumoral genetic and functional heterogeneity and cancer clinical prognosis is well-established, the underlying mechanism in CRC remains inadequately understood. Utilizing scRNA-seq data from GEO database, we re-isolated distinct subsets of cells, constructed a CRC tumor-related cell differentiation trajectory, and conducted cell-cell communication analysis to investigate potential interactions across cell clusters. A prognostic model was built by integrating scRNA-seq results with TCGA bulk RNA-seq data through univariate, LASSO, and multivariate Cox regression analyses. Eleven distinct cell types were identified, with Epithelial cells, Fibroblasts, and Mast cells exhibiting significant differences between CRC and healthy controls. T cells were observed to engage in extensive interactions with other cell types. Utilizing the 741 signature genes, prognostic risk score model was constructed. Patients with high-risk scores exhibited a significant correlation with unfavorable survival outcomes, high-stage tumors, metastasis, and low responsiveness to chemotherapy. The model demonstrated a strong predictive performance across five validation cohorts. Our investigation involved an analysis of the cellular composition and interactions of infiltrates within the microenvironment, and we developed a prognostic model. This model provides valuable insights into the prognosis and therapeutic evaluation of CRC.

Association analysis for SNPs of BAAT and COL1A1 genes with cashmere production performance and other production traits in Liaoning cashmere goats

This study aimed to investigate the relationship between the polymorphism of bile acid-CoA: amino acid N-acyltransferase (BAAT) and collagen type I alpha 1 chain (COL1A1) genes and the production performance of Liaoning Cashmere goat (LCG). The potential single nucleotide polymorphisms (SNPs) of LCG were detected by sequence comparison of BAAT and COL1A1 genes and PCR-Seq polymorphism, and the effect of SNPs on production performance was analyzed by SPSS software. The results showed that three SNPs loci were detected in BAAT gene: G7900A, T7967C, C7998T, and one SNP locus T6716C was detected in COL1AL gene. At G7900A locus, the dominant genotype for cashmere performance was GG, and the dominant genotype for body measurement traits and milk production traits was AG. At T7967C locus, the dominant genotype for cashmere performance was TT, and the dominant genotype for body measurement traits and milk production traits was CC. At C7998T locus, TT was the dominant genotype for cashmere performance, body measurement traits, and milk production traits. At the T6716C locus, TT was the dominant genotype for cashmere performance, body measurement traits, and milk production traits. H1H1: AACC is the dominant haplotype combination. Therefore, this study will provide a reliable reference for future research on cashmere production performance, body measurement traits, and milk production traits of LCG.

Induced differentiation of primordial germ cell like cells from SOX9+ porcine skin derived stem cells

Understanding the mechanisms behind porcine primordial germ cell like cells (pPGCLCs) development, differentiation, and gametogenesis is crucial in the treatment of infertility. In this study, SOX9+ skin derived stem cells (SOX9+ SDSCs) were isolated from fetal porcine skin and a high-purity SOX9+ SDSCs population was obtained. The SOX9+ SDSCs were induced to transdifferentiate into PGCLCs during 8 days of cultured. The results of RNA-seq, western blot and immunofluorescence staining verified SDSCs have the potential to transdifferentiate into PGCLCs from aspects of transcription factor activation, germ layer differentiation, energy metabolism, and epigenetic changes. Both adherent and suspended cells were collected. The adherent cells were found to be very similar to early porcine primordial germ cells (pPGCs). The suspended cells resembled late stage pPGCs and had a potential to enter meiotic process. This SDSCs culture-induced in vitro model is expected to provide suitable donor cells for stem cell transplantation in the future.

Molecular and spatial landmarks of early mouse skin development

A wealth of specialized cell populations within the skin facilitates its hair-producing, protective, sensory, and thermoregulatory functions. How the vast cell-type diversity and tissue architecture develops is largely unexplored. Here, with single-cell transcriptomics, spatial cell-type assignment, and cell-lineage tracing, we deconstruct early embryonic mouse skin during the key transitions from seemingly uniform developmental precursor states to a multilayered, multilineage epithelium, and complex dermal identity. We identify the spatiotemporal emergence of hair-follicle-inducing, muscle-supportive, and fascia-forming fibroblasts. We also demonstrate the formation of the panniculus carnosus muscle (PCM), sprouting blood vessels without pericyte coverage, and the earliest residence of mast and dendritic immune cells in skin. Finally, we identify an unexpected epithelial heterogeneity within the early single-layered epidermis and a signaling-rich periderm layer. Overall, this cellular and molecular blueprint of early skin development-which can be explored at https://kasperlab.org/tools-establishes histological landmarks and highlights unprecedented dynamic interactions among skin cells.

Glutamate-aspartate transporter dysfunction enhances aminoglycoside-induced cochlear hair cell death via NMDA receptor activation

Glutamate is a crucial neurotransmitter for hearing transduction in the cochlea, but excess glutamate is detrimental to the survival of cochlear sensory cells. Glutamate-aspartate transporter (GLAST) is the major transporter for glutamate removal; however, its role in aminoglycoside-induced hair cell loss is not well studied. In the present study, we first investigated the localization and expression of GLAST over the course of development of the mouse cochlea, and we found that inhibition of GLAST increased hair cell death. However, when the glutamate receptor NMDAR was inhibited by D-AP5, hair cell death was no longer increased by the GLAST inhibitor. Our results indicate that GLAST inhibition aggravates damage to cochlear hair cells, which may occur via NMDAR, and this suggests new clinical strategies for ameliorating the ototoxicity associated with the dysfunction of glutamate metabolism.

Temporal transcriptome features identify early skeletal commitment during human epiphysis development at single-cell resolution

Human epiphyseal development has been mainly investigated through radiological and histological approaches, uncovering few details of cellular temporal genetic alternations. Using single-cell RNA sequencing, we investigated the dynamic transcriptome changes during post-conception weeks (PCWs) 15-25 of human distal femoral epiphysis cells. We find epiphyseal cells contain multiple subtypes distinguished by specific markers, gene signatures, Gene Ontology (GO) enrichment analysis, and gene set variation analysis (GSVA). We identify the populations committed to cartilage or ossification at this time, although the secondary ossification centers (SOCs) have not formed. We describe the temporal alternation in transcriptional expression utilizing trajectories, transcriptional regulatory networks, and intercellular communication analyses. Moreover, we find the emergence of the ossification-committed population is correlated with the COL2A1-(ITGA2/11+ITGB1) signaling. NOTCH signaling may contribute to the formation of cartilage canals and ossification via NOTCH signaling. Our findings will advance the understanding of single-cell genetic changes underlying fetal epiphysis development.

Skin fibroblast functional heterogeneity in health and disease

Fibroblasts are the major cell population of connective tissue, including the skin dermis, and are best known for their function in depositing and remodelling the extracellular matrix. Besides their role in extracellular matrix homeostasis, fibroblasts have emerged as key players in many biological processes ranging from tissue immunity and wound healing to hair follicle development. Recent advances in single-cell RNA-sequencing technologies have revealed an astonishing transcriptional fibroblast heterogeneity in the skin and other organs. A key challenge in the field is to understand the functional relevance and significance of the identified new cell clusters in health and disease. Here, we discuss the functionally distinct fibroblast subtypes identified in skin homeostasis and repair and how they evolve in fibrotic disease conditions, in particular keloid scars and cancer. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

A single-cell transcriptome atlas of pig skin characterizes anatomical positional heterogeneity

Different anatomical locations of the body skin show differences in their gene expression patterns depending on different origins, and the inherent heterogeneous information can be maintained in adults. However, highly resolvable cellular specialization is less well characterized in different anatomical regions of the skin. Pig is regarded as an excellent model animal for human skin research in view of its similar physiology to human. In this study, single-cell RNA sequencing was performed on pig skin tissues from six different anatomical regions of Chenghua (CH) pigs, with a superior skin thickness trait, and the back site of large white (LW) pigs. We obtained 233,715 cells, representing seven cell types, among which we primarily characterized the heterogeneity of the top three cell types, including smooth muscle cells (SMCs), endothelial cells (ECs), and fibroblasts (FBs). Then, we further identified several subtypes of SMCs, ECs, and FBs, and discovered the expression patterns of site-specific genes involved in some important pathways such as the immune response and extracellular matrix (ECM) synthesis in different anatomical regions. By comparing differentially expressed genes of skin FBs among different anatomical regions, we considered TNN, COL11A1, and INHBA as candidate genes for facilitating ECM accumulation. These findings of heterogeneity in the main skin cell types from different anatomical sites will contribute to a better understanding of inherent skin information and place the potential focus on skin generation, transmission, and transplantation, paving the foundation for human skin priming.

An optimized force-triggered density gradient sedimentation method for isolation of pelage follicle dermal papilla cells from neonatal mouse skin

BACKGROUND

The dermal papilla cells are a specialized population of mesenchymal cells located at the base of the hair follicle (HF), which possess the capacity to regulate HF morphogenesis and regeneration. However, lack of cell-type specific surface markers restricts the isolation of DP cells and application for tissue engineering purposes.

METHODS

We describe a novel force-triggered density gradient sedimentation (FDGS) method to efficiently obtain purified follicular DP-spheres cells from neonatal mouse back skin, utilizing only centrifugation and optimized density gradients.

RESULTS

Expression of characteristic DP cell markers, alkaline phosphatase, β-catenin, versican, and neural cell adhesion molecules, were confirmed by immunofluorescence. Further, the patch assays demonstrated that DP cells maintained their hair regenerative capacity in vivo. Compared with current methods, including microdissection and fluorescence-activated cell sorting, the FDGS technique is simpler and more efficient for isolating DP cells from neonatal mouse skin.

CONCLUSIONS

The FDGS method will improve the research potential of neonatal mouse pelage-derived DP cells for tissue engineering purposes.

Murine skin-derived multipotent papillary dermal fibroblast progenitors show germline potential in vitro

BACKGROUND

Many laboratories have described the in vitro isolation of multipotent cells with stem cell properties from the skin of various species termed skin-derived stem cells (SDSCs). However, the cellular origin of these cells and their capability to give rise, among various cell types, to male germ cells, remain largely unexplored.

METHODS

SDSCs were isolated from newborn mice skin, and then differentiated into primordial germ cell-like cells (PGCLCs) in vitro. Single-cell RNA sequencing (scRNA-seq) was then applied to dissect the cellular origin of SDSCs using cells isolated from newborn mouse skin and SDSC colonies. Based on an optimized culture strategy, we successfully generated spermatogonial stem cell-like cells (SSCLCs) in vitro.

RESULTS

Here, using scRNA-seq and analyzing the profile of 7543 single-cell transcriptomes from newborn mouse skin and SDSCs, we discovered that they mainly consist of multipotent papillary dermal fibroblast progenitors (pDFPs) residing in the dermal layer. Moreover, we found that epidermal growth factor (EGF) signaling is pivotal for the capability of these progenitors to proliferate and form large colonies in vitro. Finally, we optimized the protocol to efficiently generate PGCLCs from SDSCs. Furthermore, PGCLCs were induced into SSCLCs and these SSCLCs showed meiotic potential when cultured with testicular organoids.

CONCLUSIONS

Our findings here identify pDFPs as SDSCs derived from newborn skin and show for the first time that such precursors can be induced to generate cells of the male germline.

Single-nucleus RNA and ATAC sequencing uncovers the molecular and cellular characteristics in the musk gland of Chinese forest musk deer (Moschus berezovskii)

The Chinese forest musk deer (FMD; Moschus berezovskii) is an endangered artiodactyl mammal. Musk secreted by the musk gland of male has extremely high economic and medicinal value. However, the molecular and cellular characteristics of the musk gland have not been studied. Here, we investigated the diversity and transcriptional composition of musk gland cell types and the effect of cell type-specific chromatin accessibility on gene expression using single-nucleus RNA sequencing (snRNA-seq) and single-nucleus ATAC sequencing (snATAC-seq) association analysis. Based on uniform manifold approximation and projection (UMAP) analysis, we identified 13 cell types from the musk gland, which included two different acinar cells (cluster 0 and cluster 10). Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that many pathways related to musk secretion were enriched in acinar cells. Our analysis also revealed acinar cell core transcription factors and core target genes, and further constructed acinar cell-specific regulatory networks. In cluster 0, 11 core target genes (Nedd4l, Adcy9, Akr1c1, Vapb, Me1, Acsl1, Acss3, Srd5a1, Scnn1a, Acadm, and Nceh1) possibly related to musk secretion were regulated by 24 core transcription factors (SP3, NFIC, NR6A1, EHF, RUNX1, TFAP2A, RREB1, GRHL2, NFIB, ELF1, MAX, KLF5, REL, HES1, POU2F3, TFDP1, NR2C1, ATF7, MEIS1, NR4A2, NFIA, PBX1, ZNF652, and NFKB1). In cluster 10, four core target genes (Akr1c1, Pcca, Atp1b1, and Sgk1) possibly related to musk secretion were regulated by 10 core transcription factors (BARX2, EHF, PBX1, RUNX1, NFIB, FOXP1, KLF3, KLF6, ETV6, and NR3C2). Moreover, the credibility of snRNA-seq and snATAC-seq data was verified by fluorescence in situ hybridization and immunohistochemistry. Finally, cell communication analysis demonstrated that the two types of acinar cells mainly have communications in musk secretion-related processes. In conclusion, we provided important insights and invaluable resources for the molecular and cellular characteristics of the musk gland, which will lay a foundation for the study of musk secretion mechanism in the future.

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.

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.

Single-cell chromatin landscapes of mouse skin development

The coat of mammals is produced by hair follicles, and hair follicle is an important and complex accessory organ of skin. As a complex physiological regulation process, hair follicle morphogenesis is regulated by a series of signal pathway factors, involves the interaction of multiple cell types and begins in the early embryonic stage. However, its transcriptional regulatory mechanism is unclear. We have therefore utilized single-cell ATAC sequencing to obtain the chromatin accessibility landscapes of 6,928, 6,961 and 7,374 high-quality cells from the dorsal skins of E13.5, E16.5 and P0 mice (Mus musculus), respectively. Based on marker gene activity clustering, we defined 6, 8 and 5 distinct cell types in E13.5, E16.5 and P0 stages, respectively. Furtherly, we integrated the fibroblasts and keratinocytes clusters, performed further analysis and re-clustered. The single cell map of the chromatin open area was drawn from each cell type and the mechanism of cell transcription regulation was explored. Collectively, our data provide a reference for deeply exploring the epigenetic regulation mechanism of mouse hair follicles development.

Investigation of hub genes and immune infiltration in androgenetic alopecia using bioinformatics analysis

Background

Androgenetic alopecia (AGA) is a type of non-scarring hair loss. Current drugs for AGA are accompanied by adverse reactions and a high recurrence rate. Thus, the discovery of diagnostic biomarkers and therapeutic targets for AGA remains imperatively warranted.

Methods

The GSE90594 dataset, which contained scalp skin biopsies from 14 male AGA cases and healthy volunteers, was used to identify the differentially expressed genes (DEGs). Functional enrichment analysis was subsequently performed. Next, the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database combined with the cytoHubba plugin of Cytoscape were used to obtain the key genes of AGA. Thereafter, the Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) algorithm was performed to evaluate the relative abundance of immune cells between male AGA patients and healthy controls. The correlation between key genes and infiltrating immune cells was analyzed to obtain the significant immune-cell related genes (IRGs), then intersected with the DEGs between immortalized balding and non-balding human dermal papilla cells (DPCs) of the GSE93766 dataset as well as the DEGs obtained by the GSE90594 dataset, thus obtaining the hub genes of AGA. Finally, the hub genes were validated using GSE36169, which contained expression profiling of tissues biopsied from haired and bald scalps of five individuals with AGA.

Results

A total of 234 DEGs were obtained from the GSE90594 dataset, which were mainly enriched in the extracellular matrix (ECM)-related pathways and immune-related activities. The STRING database and ten algorithms in the cytoHubba plugin of Cytoscape disclosed 21 key DEGs. The results of the CIBERSORT algorithm revealed the relative abundances of 20 kinds of immune cells between diseased and healthy individuals, and yielded 15 IRGs involved in the pathogenesis of AGA. Next, the intersection analysis identified four hub genes of AGA, comprising COL1A2, PCOLCE, ITGAX, and LOX. The GSE36169 dataset validated the expression pattern of hub genes in the haired scalp of AGA patients.

Conclusions

We discovered that the hub genes identified are closely linked with the causative factors of AGA, which could be used as the viable diagnostic and therapeutic target in the clinical applications.

Wound healing, fibroblast heterogeneity, and fibrosis

Fibroblasts are highly dynamic cells that play a central role in tissue repair and fibrosis. However, the mechanisms by which they contribute to both physiologic and pathologic states of extracellular matrix deposition and remodeling are just starting to be understood. In this review article, we discuss the current state of knowledge in fibroblast biology and heterogeneity, with a primary focus on the role of fibroblasts in skin wound repair. We also consider emerging techniques in the field, which enable an increasingly nuanced and contextualized understanding of these complex systems, and evaluate limitations of existing methodologies and knowledge. Collectively, this review spotlights a diverse body of research examining an often-overlooked cell type-the fibroblast-and its critical functions in wound repair and beyond.

Hedgehog signaling reprograms hair follicle niche fibroblasts to a hyper-activated state

Hair follicle stem cells are regulated by dermal papilla fibroblasts, their principal signaling niche. Overactivation of Hedgehog signaling in the niche dramatically accelerates hair growth and induces follicle multiplication in mice. On single-cell RNA sequencing, dermal papilla fibroblasts increase heterogeneity to include new Wnt5a^high states. Transcriptionally, mutant fibroblasts activate regulatory networks for Gli1, Alx3, Ebf1, Hoxc8, Sox18, and Zfp239. These networks jointly upregulate secreted factors for multiple hair morphogenesis and hair-growth-related pathways. Among these is non-conventional TGF-β ligand Scube3. We show that in normal mouse skin, Scube3 is expressed only in dermal papillae of growing, but not in resting follicles. SCUBE3 protein microinjection is sufficient to induce new hair growth, and pharmacological TGF-β inhibition rescues mutant hair hyper-activation phenotype. Moreover, dermal-papilla-enriched expression of SCUBE3 and its growth-activating effect are partially conserved in human scalp hair follicles. Thus, Hedgehog regulates mesenchymal niche function in the hair follicle via SCUBE3/TGF-β mechanism.

3D printing of skin equivalents with hair follicle structures and epidermal-papillary-dermal layers using gelatin/hyaluronic acid hydrogels

Recent advances in three-dimensional (3D) bioprinting technologies enabled the fabrication of sophisticated live 3D tissue analogs. Although various hydrogel-based bioink has been reported, the development of advanced bioink materials that can reproduce the composition of native extracellular matrix (ECM) accurately and mimic the intrinsic property of laden cells is still challenging. In this work, 3D printed skin equivalents incorporating hair follicle structures and epidermal-papillary-dermal layers are fabricated with gelatin methacryloyl (GelMA)/hyaluronic acid (HA) MA (HAMA) hydrogel (GelMA/HAMA) bioink. The composition of collagen and glycosaminoglycan (GAG) of native skin was recapitulated by adjusting the combination of GelMA and HAMA. The GelMA/HAMA bioink was proven to have excellent viscoelastic and physicochemical properties, 3D printability, cytocompatibility, and functionality to maintain the hair inductive potency and facilitated spontaneous hair pore development. Overall, we suggest that the GelMA/HAMA hydrogels can be promising candidates as bioinks for the 3D printing of skin equivalents with epidermal-papillary-dermal multi-layers and hair follicle structures, and they might serve as a useful model in skin tissue engineering and regeneration.

Cellular Heterogeneity Facilitates the Functional Differences Between Hair Follicle Dermal Sheath Cells and Dermal Papilla Cells: A New Classification System for Mesenchymal Cells within the Hair Follicle Niche

Mesenchymal stem cells (MSCs) are known for their self-renewal and multi-lineage differentiation potential, with these cells often being evaluated in the regulation and maintenance of specific cellular niches including those of the hair follicle. Most mesenchymal stem cells in the hair follicles are housed in the dermal papilla (DP) and dermal sheath (DS), with both niches characterized by a broad variety of cellular subsets. However, while most previous studies describing the hair follicle mesenchymal niche treated all DP and DS cells as Hair Follicle Mesenchymal Stem Cells (HF-MSCs), the high number of cellular subsets would suggest that these cells are actually too heterogenous for such a broad definition. Given this we designed this study to evaluate the differentiation processes in these cells and used this data to create a new set of classifications for DP and DS cells, dividing them into "hair follicle mesenchymal stem cells (HF-MSCs)", "hair follicle mesenchymal progenitor cells (HF-MPCs)", and "hair follicle mesenchymal functional cells (HF-MFCs)". In addition, those cells that possess self-renewal and differentiation were re-named hair follicle derived mesenchymal multipotent cells (HF-MMCs). This new classification may help to further our understanding of the heterogeneity of hair follicle dermal cells and provide new insights into their evaluation.

Gene network analysis reveals candidate genes related with the hair follicle development in sheep

Background

Merino sheep are the most famous fine wool sheep in the world. They have high wool production and excellent wool quality and have attracted worldwide attention. The fleece of the Merino sheep is composed predominantly of wool fibers grown from secondary wool follicles. Therefore, it is necessary to study the development of hair follicles to understand the mechanism of wool production. The hair follicle is a complex biological system involved in a dynamic process governed by gene regulation. The hair follicle development process is very complex and poorly understood. The purpose of our research is to identify candidate genes related to hair follicle development, provide a theoretical molecular breeding basis for the cultivation of fine wool sheep, and provide a reference for the problems of hair loss and alopecia areata that affect human beings.

Results

We analyzed mRNAs data in skin tissues of 18 Merino sheep at four embryonic days (E65, E85, E105 and E135) and two postnatal days (P7 and P30). G1 to G6 represent hair follicles developmental at six stages (i.e. E65 to P30). We identified 7879 differentially expressed genes (DEGs) and 12623 novel DEGs, revealed different expression patterns of these DEGs at six stages of hair follicle development, and demonstrated their complex interactions. DEGs with stage-specific expression were significantly enriched in epidermal differentiation and development, hair follicle development and hair follicle morphogenesis and were enriched in many pathways related to hair follicle development. The key genes (LAMA5, WNT10A, KRT25, SOSTDC1, ZDHHC21, FZD1, BMP7, LRP4, TGFβ2, TMEM79, SOX10, ITGB4, KRT14, ITGA6, and GLI2) affecting hair follicle morphogenesis were identified by network analysis.

Conclusion

This study provides a new reference for the molecular basis of hair follicle development and lays a foundation for further improving sheep hair follicle breeding. Candidate genes related to hair follicular development were found, which provided a theoretical basis for molecular breeding for the culture of fine wool sheep. These results are a valuable resource for biological investigations of fleece evolution in animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08552-2.

Fibroblast Heterogeneity in Healthy and Wounded Skin

Fibroblasts are the main cell type in the dermis. They are responsible for the synthesis and deposition of structural proteins such as collagen and elastin, which are integrated into the extracellular matrix (ECM). Mouse and human studies using flow cytometry, cell culture, skin reconstitution, and lineage tracing experiments have shown the existence of different subpopulations of fibroblasts, including papillary fibroblasts, reticular fibroblasts, and fibroblasts comprising the dermal papilla at the base of the hair follicle. In recent years, the technological advances in single-cell sequencing have allowed researchers to study the repertoire of cells present in full-thickness skin including the dermis. Multiple groups have confirmed that distinct fibroblast populations can be identified in mouse and human dermis on the basis of differences in the transcriptional profile. Here, we discuss the current state of knowledge regarding dermal fibroblast heterogeneity in healthy mouse and human skin, highlighting the similarities and differences between mouse and human fibroblast subpopulations. We also discuss how fibroblast heterogeneity may provide insights into physiological wound healing and its dysfunction in pathological states such as hypertrophic and keloid scars.

Morphogenesis, Growth Cycle and Molecular Regulation of Hair Follicles

As one of the main appendages of skin, hair follicles play an important role in the process of skin regeneration. Hair follicle is a tiny organ formed by the interaction between epidermis and dermis, which has complex and fine structure and periodic growth characteristics. The hair growth cycle is divided into three continuous stages, growth (anagen), apoptosis-driven regression (catagen) and relative quiescence (telogen). And The Morphogenesis and cycle of hair follicles are regulated by a variety of signal pathways. When the signal molecules in the pathways are abnormal, it will affect the development and cycle of hair follicles, which will lead to hair follicle-related diseases.This article will review the structure, development, cycle and molecular regulation of hair follicles, in order to provide new ideas for solving diseases and forming functional hair follicle.

Epithelial-Mesenchymal Interaction in Hair Regeneration and Skin Wound Healing

Interactions between epithelial and mesenchymal cells influence hair follicles (HFs) during embryonic development and skin regeneration following injury. Exchanging soluble molecules, altering key pathways, and extracellular matrix signal transduction are all part of the interplay between epithelial and mesenchymal cells. In brief, the mesenchyme contains dermal papilla cells, while the hair matrix cells and outer root sheath represent the epithelial cells. This study summarizes typical epithelial–mesenchymal signaling molecules and extracellular components under the control of follicular stem cells, aiming to broaden our current understanding of epithelial–mesenchymal interaction mechanisms in HF regeneration and skin wound healing.

Integrative analysis of circRNAs from Yangtze River Delta white goat neck skin tissue by high-throughput sequencing (circRNA-seq)

The Yangtze River Delta white goat is a unique goat species that can produce superior-quality brush hair. The formation of this brush hair is controlled by a series of critical genes and related signaling pathways. Circular RNAs (circRNAs), are ubiquitous endogenous non-coding RNAs that regulate many biological and physiological processes in mammals. However, little is known about the potential regulatory role of circRNAs on superior-quality brush hair formation in Yangtze River Delta white goat. In this study, high-throughput sequencing technology was used to only detect circRNAs in the neck skin tissue of normal-quality goats (NHQs) and superior-quality goats (HQs). A total of 61 803 circRNAs were identified and 32 of them were differentially expressed in the NHQ group vs. the HQ group. Functional enrichment analysis showed that the source gene of differentially expressed circRNAs (DE-circRNAs) was enriched mostly in platelet activation and the focal adhesion signal pathway. Action mechanism analysis revealed that DE-circRNAs could sponge to many identified miRNAs, including miR-31, miR-125b, miR-let-7a and miR-149-5p, which have important roles in goat hair follicle stem cell growth, hair follicle development and morphogenesis. Altogether, our findings provide a valuable basis for studying circRNAs involved in superior-quality brush hair traits and meanwhile advance our understanding of circRNA complex regulation mechanisms in Yangtze River Delta white goat skin hair follicle development.

Single-Cell Transcriptomics Reveals the Molecular Anatomy of Sheep Hair Follicle Heterogeneity and Wool Curvature

Wool is the critical textile raw material which is produced by the hair follicle of sheep. Therefore, it has important implications to investigate the molecular mechanism governing hair follicle development. Due to high cellular heterogeneity as well as the insufficient cellular, molecular, and spatial characterization of hair follicles on sheep, the molecular mechanisms involved in hair follicle development and wool curvature of sheep remains largely unknown. Single-cell RNA sequencing (scRNA-seq) technologies have made it possible to comprehensively dissect the cellular composition of complex skin tissues and unveil the differentiation and spatial signatures of epidermal and hair follicle development. However, such studies are lacking so far in sheep. Here, single-cell suspensions from the curly wool and straight wool lambskins were prepared for unbiased scRNA-seq. Based on UAMP dimension reduction analysis, we identified 19 distinct cell populations from 15,830 single-cell transcriptomes and characterized their cellular identity according to specific gene expression profiles. Furthermore, novel marker gene was applied in identifying dermal papilla cells isolated in vitro. By using pseudotime ordering analysis, we constructed the matrix cell lineage differentiation trajectory and revealed the dynamic gene expression profiles of matrix progenitors' commitment to the hair shaft and inner root sheath (IRS) cells. Meanwhile, intercellular communication between mesenchymal and epithelial cells was inferred based on CellChat and the prior knowledge of ligand–receptor pairs. As a result, strong intercellular communication and associated signaling pathways were revealed. Besides, to clarify the molecular mechanism of wool curvature, differentially expressed genes in specific cells between straight wool and curly wool were identified and analyzed. Our findings here provided an unbiased and systematic view of the molecular anatomy of sheep hair follicle comprising 19 clusters; revealed the differentiation, spatial signatures, and intercellular communication underlying sheep hair follicle development; and at the same time revealed the potential molecular mechanism of wool curvature, which will give important new insights into the biology of the sheep hair follicle and has implications for sheep breeding.

A tissue-bioengineering strategy for modeling rare human kidney diseases in vivo

The lack of animal models for some human diseases precludes our understanding of disease mechanisms and our ability to test prospective therapies in vivo. Generation of kidney organoids from Tuberous Sclerosis Complex (TSC) patient-derived-hiPSCs allows us to recapitulate a rare kidney tumor called angiomyolipoma (AML). Organoids derived from TSC2-/- hiPSCs but not from isogenic TSC2+/- or TSC2+/+ hiPSCs share a common transcriptional signature and a myomelanocytic cell phenotype with kidney AMLs, and develop epithelial cysts, replicating two major TSC-associated kidney lesions driven by genetic mechanisms that cannot be consistently recapitulated with transgenic mice. Transplantation of multiple TSC2-/- renal organoids into the kidneys of immunodeficient rats allows us to model AML in vivo for the study of tumor mechanisms, and to test the efficacy of rapamycin-loaded nanoparticles as an approach to rapidly ablate AMLs. Collectively, our experimental approaches represent an innovative and scalable tissue-bioengineering strategy for modeling rare kidney disease in vivo.

A Single-cell Transcriptome Atlas of Cashmere Goat Hair Follicle Morphogenesis

Cashmere, also known as soft gold, is produced from the secondary hair follicles (SHFs) of cashmere goats. The number of SHFs determines the yield and quality of cashmere; therefore, it is of interest to investigate the transcriptional profiles present during cashmere goat hair follicle development. However, mechanisms underlying this development process remain largely unexplored, and studies regarding hair follicle development mostly use a murine research model. In this study, to provide a comprehensive understanding of cellular heterogeneity and cell fate decisions, single-cell RNA sequencing was performed on 19,705 single cells of the dorsal skin from cashmere goat fetuses at induction (embryonic day 60; E60), organogenesis (E90), and cytodifferentiation (E120) stages. For the first time, unsupervised clustering analysis identified 16 cell clusters, and their corresponding cell types were also characterized. Based on lineage inference, a detailed molecular landscape was revealed along the dermal and epidermal cell lineage developmental pathways. Notably, our current data also confirmed the heterogeneity of dermal papillae from different hair follicle types, which was further validated by immunofluorescence analysis. The current study identifies different biomarkers during cashmere goat hair follicle development and has implications for cashmere goat breeding in the future.

A hypothetical role for autophagy during the day/night rhythm-regulated melatonin synthesis in the rat pineal gland

Melatonin is a highly conserved molecule that regulates day/night rhythms; it is associated with sleep improvement, reactive oxygen species (ROS) scavenging, anti-aging effects, and seasonal and circadian rhythms and has been a hot topic of research for decades. Using single-cell RNA sequencing, a recent study describes a single-cell transcriptome atlas for the rat pineal gland. Based on a more comprehensive analysis of the retrieved data (Mays et al., PLoS One, 2018, 13, e0205883), results from the current study unveiled the underappreciated gene regulatory network behind different cell populations in the pineal gland. More importantly, our study here characterized, for the first time, the day/night activation of autophagy flux in the rat pineal gland, indicating a potential role of autophagy in regulating melatonin synthesis in the rat pineal gland. These findings emphasized a hypothetical role of day/night autophagy in linking the biological clock with melatonin synthesis. Furthermore, ultrastructure analysis of pinealocytes provided fascinating insights into differences in their intracellular structure between daytime and nighttime. In addition, we also provide a preliminary description of cell-cell communication in the rat pineal gland. In summary, the current study unveils the day/night regulation of autophagy in the rat pineal gland, raising a potential role of autophagy in day/night-regulated melatonin synthesis.

High-resolution single-cell transcriptomics reveals heterogeneity of self-renewing hair follicle stem cells

Multipotent bulge stem cells (SCs) fuel the hair follicle (HF) cyclic growth during adult skin homeostasis, but their intrinsic molecular heterogeneity is not well understood. These hair follicle stem cells (HFSCs) engage in bouts of self-renewal, migration and differentiation during the hair cycle. Here, we perform high-resolution single-cell RNA sequencing (scRNA-seq) of HFSCs sorted as CD34+ /K14-H2BGFP+ from mouse skin at mid-anagen, the self-renewal stage. We dissect the transcriptomic profiles and unravel population-specific transcriptional heterogeneity. Unsupervised clustering reveals five major HFSC populations, which distinguished by known markers associated with both the bulge and the outer root sheath (ORS) underneath. These populations include quiescent bulge, ORS cellular states and proliferative cells. Lineage trajectory analysis predicted the prospective differentiation path of these cellular states and their corresponding self-renewing subpopulations. The bulge population itself can be further sub-divided into distinct subpopulations that can be mapped to the upper, mid and lower bulge regions, and present a decreasing quiescence score. Gene set enrichment analysis (GSEA) revealed new markers and suggested potentially distinct functions of the ORS and bulge subpopulations. This included communications between the upper bulge subpopulation and sensory nerves and between the upper ORS and skin vasculature, as well as enrichment of a bulge subset in cell migratory functions. The lower ORS enriched genes may potentially enable nutrients passing from the surrounding fat and vasculature cells towards the proliferating hair matrix cells. Thus, we provide a comprehensive account of HFSC molecular heterogeneity during their self-renewing stage, which enables future HF functional studies.

Infertility network and hub genes for nonobstructive azoospermia utilizing integrative analysis

Non-obstructive azoospermia (NOA) is the most severe form of male infertility owing to the absence of sperm during ejaculation as a result of failed spermatogenesis. The molecular mechanisms of NOA have not been well studied. Here, we revealed the dysregulated differentially expressed genes in NOA and related signaling pathways or biological processes. Cluster features of biological processes include spermatogenesis, fertilization, cilium movement, penetration of zona pellucida, sperm chromatin condensation, and being significantly enriched metabolic pathways in proximal tubule bicarbonate reclamation, aldosterone synthesis and secretion, glycolysis and glycogenesis pathways in NOA using Gene Ontology analysis and pathway enrichment analysis. The NOA gene co-expression network was constructed by weighted gene co-expression network analysis to identify the hub genes (CHD5 and SPTBN2). In addition, we used another Gene Expression Omnibus dataset (GSE45887) to validate these hub genes. Furthermore, we used the Seurat package to classify testicular tissue cells from NOA patients and to characterize the differential expression of hub genes in different cell types from different adult males based on the scRNA-seq dataset (GSE106487). These results provide new insights into the pathogenesis of NOA. Of particular note, CHD5 and SPTBN2 may be potential biomarkers for the diagnosis and treatment of NOA.

Inference and analysis of cell-cell communication using CellChat

Understanding global communications among cells requires accurate representation of cell-cell signaling links and effective systems-level analyses of those links. We construct a database of interactions among ligands, receptors and their cofactors that accurately represent known heteromeric molecular complexes. We then develop CellChat, a tool that is able to quantitatively infer and analyze intercellular communication networks from single-cell RNA-sequencing (scRNA-seq) data. CellChat predicts major signaling inputs and outputs for cells and how those cells and signals coordinate for functions using network analysis and pattern recognition approaches. Through manifold learning and quantitative contrasts, CellChat classifies signaling pathways and delineates conserved and context-specific pathways across different datasets. Applying CellChat to mouse and human skin datasets shows its ability to extract complex signaling patterns. Our versatile and easy-to-use toolkit CellChat and a web-based Explorer ( http://www.cellchat.org/ ) will help discover novel intercellular communications and build cell-cell communication atlases in diverse tissues.

Identification of Key Pathways and Genes Related to the Development of Hair Follicle Cycle in Cashmere Goats

The development of hair follicle in cashmere goats shows significant periodic change, as with mice and humans. However, for cashmere goat with double-coat, the periodic change may be due to other regulatory molecules and signal pathways. To understand the mechanism of periodic development of hair follicle, we performed a weighted gene coexpression network analysis (WGCNA) to mine key genes and establish an interaction network by utilizing the NCBI public dataset. Ten coexpression modules, including 7689 protein-coding genes, were constructed by WGCNA, six of which are considered to be significantly related to the development of the hair follicle cycle. A functional enrichment analysis for each model showed that they are closely related to ECM- receptor interaction, focal adhesion, PI3K-Akt signaling pathway, estrogen signaling pathway, and so on. Combined with the analysis of differential expressed genes, 12 hub genes from coexpression modules were selected as candidate markers, i.e., COL1A1, C1QTNF6, COL1A2, AQP3, KRTAP3-1, KRTAP11-1, FA2H, NDUFS5, DERL2, MRPL14, ANTKMT and XAB2, which might be applied to improve cashmere production.