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

Advances in single-cell transcriptomics in animal research

Understanding biological mechanisms is fundamental for improving animal production and health to meet the growing demand for high-quality protein. As an emerging biotechnology, single-cell transcriptomics has been gradually applied in diverse aspects of animal research, offering an effective method to study the gene expression of high-throughput single cells of different tissues/organs in animals. In an unprecedented manner, researchers have identified cell types/subtypes and their marker genes, inferred cellular fate trajectories, and revealed cell‒cell interactions in animals using single-cell transcriptomics. In this paper, we introduce the development of single-cell technology and review the processes, advancements, and applications of single-cell transcriptomics in animal research. We summarize recent efforts using single-cell transcriptomics to obtain a more profound understanding of animal nutrition and health, reproductive performance, genetics, and disease models in different livestock species. Moreover, the practical experience accumulated based on a large number of cases is highlighted to provide a reference for determining key factors (e.g., sample size, cell clustering, and cell type annotation) in single-cell transcriptomics analysis. We also discuss the limitations and outlook of single-cell transcriptomics in the current stage. This paper describes the comprehensive progress of single-cell transcriptomics in animal research, offering novel insights and sustainable advancements in agricultural productivity and animal health.

Single-Cell Sequencing Technology in Ruminant Livestock: Challenges and Opportunities

Advancements in single-cell sequencing have transformed the genomics field by allowing researchers to delve into the intricate cellular heterogeneity within tissues at greater resolution. While single-cell omics are more widely applied in model organisms and humans, their use in livestock species is just beginning. Studies in cattle, sheep, and goats have already leveraged single-cell and single-nuclei RNA-seq as well as single-cell and single-nuclei ATAC-seq to delineate cellular diversity in tissues, track changes in cell populations and gene expression over developmental stages, and characterize immune cell populations important for disease resistance and resilience. Although challenges exist for the use of this technology in ruminant livestock, such as the precise annotation of unique cell populations and spatial resolution of cells within a tissue, there is vast potential to enhance our understanding of the cellular and molecular mechanisms underpinning traits essential for healthy and productive livestock. This review intends to highlight the insights gained from published single-cell omics studies in cattle, sheep, and goats, particularly those with publicly accessible data. Further, this manuscript will discuss the challenges and opportunities of this technology in ruminant livestock and how it may contribute to enhanced profitability and sustainability of animal agriculture in the future.

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.

A single-cell transcriptome atlas of Lueyang black-bone chicken skin

As the largest organ of the body, the skin participates in various physiological activities, such as barrier function, sensory function, and temperature regulation, thereby maintaining the balance between the body and the natural environment. To date, compositional and transcriptional profiles in chicken skin cells have not been reported. Here, we report detailed transcriptome analyses of cell populations present in the skin of a black-feather chicken and a white-feather chicken using single-cell RNA sequencing (scRNA-seq). By analyzing cluster-specific gene expression profiles, we identified 12 cell clusters, and their corresponding cell types were also characterized. Subsequently, we characterized the subpopulations of keratinocytes, myocytes, mesenchymal cells, fibroblasts, and melanocytes. It is worth noting that we have identified a subpopulation of keratinocytes involved in pigment granule capture and a subpopulation of melanocytes involved in pigment granule deposition, both of which have a higher cell abundance in black-feather chicken compared to white-feather chicken. Meanwhile, we also compared the cellular heterogeneity features of Lueyang black-bone chicken skin with different feather colors. In addition, we also screened out 12 genes those could be potential markers of melanocytes. Finally, we validated the specific expression of SGK1, WNT5A, CTSC, TYR, and LAPTM5 in black-feather chicken, which may be the key candidate genes determining the feather color differentiation of Lueyang black-bone chicken. In summary, this study first revealed the transcriptome characteristics of chicken skin cells via scRNA-seq technology. These datasets provide valuable information for the study of avian skin characteristics and have important implications for future poultry breeding.

Single-cell RNA sequencing and UPHLC-MS/MS targeted metabolomics offer new insights into the etiological basis for male cattle-yak sterility

The variety of species can be efficiently increased by interspecific hybridization. However, because the males in the hybrid progeny are usually sterile, this heterosis cannot be employed when other cattle and yaks are hybridized. While some system-level studies have sought to explore the etiological basis for male cattle-yak sterility, no systematic cellular analyses of this phenomenon have yet been performed. Here, single-cell RNA sequencing and UPHLC-MS/MS targeted metabolomics methods were used to study the differences in testicular tissue between 4-year-old male yak and 4-year-old male cattle-yak, providing new and comprehensive insights into the causes of male cattle-yak sterility. Cattle-yak testes samples detected 6 somatic cell types and one mixed germ cell type. Comparisons of these cell types revealed the more significant differences in Sertoli cells (SCs) and [Leydig cells and myoid cells (LCs_MCs)] between yak and cattle-yak samples compared to other somatic cell clusters. Even though the LCs and MCs from yaks and cattle-yaks were derived from the differentiation of the same progenitor cells, a high degree of overlap between LCs and MCs was observed in yak samples. Still, only a small overlap between LCs and MCs was observed in cattle-yak samples. Functional enrichment analyses revealed that genes down-regulated in cattle-yak SCs were primarily enriched in biological activity, whereas up-regulated genes in these cells were enriched for apoptotic activity. Furthermore, the genes of up-regulated in LCs_MCs of cattle-yak were significantly enriched in enzyme inhibitor and molecular function inhibitor activity. On the other hand, the genes of down-regulated in these cells were enriched for signal receptor binding, molecular function regulation, positive regulation of biological processes, and regulation of cell communication activity. The most significant annotated differences between yak and cattle-yak LCs_MCs were associated with cell-to-cell communication. While yak LCs_MCs regulated spermatogenic cells at spermatogonia, spermatocyte, and spermatid levels, no such relationships were found between cattle-yak LCs_MCs and germ cells. This may suggest that the somatic niche in male cattle-yak testes is a microenvironment that is ultimately not favorable for spermatogenesis.

Interpretation of the Yak Skin Single-Cell Transcriptome Landscape

The morphogenesis of hair follicle structure is accompanied by the differentiation of skin tissue. Mammalian coats are produced by hair follicles. The formation of hair follicles requires signal transmission between the epidermis and dermis. However, knowledge of the transcriptional regulatory mechanism is still lacking. We used single-cell RNA sequencing to obtain 26,573 single cells from the scapular skin of yaks at hair follicle telogen and anagen stages. With the help of known reference marker genes, 11 main cell types were identified. In addition, we further analyzed the DP cell and dermal fibroblast lineages, drew a single-cell map of the DP cell and dermal fibroblast lineages, and elaborated the key genes, signals, and functions involved in cell fate decision making. The results of this study provide a very valuable resource for the analysis of the heterogeneity of DP cells and dermal fibroblasts in the skin and provide a powerful theoretical reference for further exploring the diversity of hair follicle cell types and hair follicle morphogenesis.

SOX18 Promotes the Proliferation of Dermal Papilla Cells via the Wnt/β-Catenin Signaling Pathway

SRY-box transcription factor 18 (SOX18) is known to play a crucial role in the growth and development of hair follicles (HF) in both humans and mice. However, the specific effect of SOX18 on sheep hair follicles remains largely unknown. In our previous study, we observed that SOX18 was specifically expressed within dermal papilla cells (DPCs) in ovine hair follicles, leading us to investigate its potential role in the growth of hair follicles in sheep. In the present study, we aimed to examine the effect of SOX18 in DPCs and preliminarily study its regulatory mechanism through RNA-seq. We initially found that the overexpression of SOX18 promoted the proliferation of DPCs compared to the negative control group, while the interference of SOX18 had the opposite effect. To gain further insight into the regulatory mechanism of SOX18, we conducted RNA-seq analysis after knocking down SOX18 in Hu sheep DPCs. The result showed that the Wnt/β-Catenin signaling pathway was involved in the growth process of DPC after SOX18 knockdown. Subsequently, we investigated the effect of SOX18 on the Wnt/β-Catenin signaling pathway in DPCs using TOP/FOP-flash, qRT-PCR, and Western blot (WB) analysis. Our data demonstrated that SOX18 could activate the Wnt/β-Catenin signaling pathway in DPCs. Additionally, we observed that SOX18 could rescue the proliferation of DPCs after inhibiting the Wnt/β-Catenin signaling pathway. These findings underscore the essential role of SOX18 as a functional molecule governing the proliferation of DPCs. Additionally, these findings also greatly enhance our understanding of the role of SOX18 in the proliferation of DPCs and the growth of wool in Hu sheep.

Developmental Mapping of Hair Follicles in the Embryonic Stages of Cashmere Goats Using Proteomic and Metabolomic Construction

The hair follicle (HF) is the fundamental unit for fleece and cashmere production in cashmere goats and is crucial in determining cashmere yield and quality. The mechanisms regulating HF development in cashmere goats during the embryonic period remain unclear. Growing evidence suggests that HF development involves complex developmental stages and critical events, and identifying the underlying factors can improve our understanding of HF development. In this study, samples were collected from embryonic day 75 (E75) to E125, the major HF developmental stages. The embryonic HFs of cashmere goats were subjected to proteomic and metabolomic analyses, which revealed dynamic changes in the key factors and signalling pathways controlling HF development at the protein and metabolic levels. Gene ontology and the Kyoto Encyclopaedia of Genes and Genomes were used to functionally annotate 1784 significantly differentially expressed proteins and 454 significantly differentially expressed metabolites enriched in different HF developmental stages. A joint analysis revealed that the oxytocin signalling pathway plays a sustained role in embryonic HF development by activating the MAPK and Ca2+ signalling pathways, and a related regulatory network map was constructed. This study provides a global perspective on the mechanism of HF development in cashmere goats and enriches our understanding of embryonic HF development.

Transcriptomic Profiling of Meat Quality Traits of Skeletal Muscles of the Chinese Indigenous Huai Pig and Duroc Pig

The Huai pig is a well-known indigenous pig breed in China. The main advantages of Huai pigs over Western commercial pig breeds include a high intramuscular fat (IMF) content and good meat quality. There are significant differences in the meat quality traits of the same muscle part or different muscle parts of the same variety. To investigate the potential genetic mechanism underlying the meat quality differences in different pig breeds or muscle groups, longissimus dorsi (LD), psoas major (PM), and biceps femoris (BF) muscle tissues were collected from two pig breeds (Huai and Duroc). There were significant differences in meat quality traits and amino acid content. We assessed the muscle transcriptomic profiles using high-throughput RNA sequencing. The IMF content in the LD, PM, and BF muscles of Huai pigs was significantly higher than that in Duroc pigs (p < 0.05). Similarly, the content of flavor amino acids in the three muscle groups was significantly higher in Huai pigs than that in Duroc pigs (p < 0.05). We identified 175, 110, and 86 differentially expressed genes (DEGs) between the LD, PM, and BF muscles of the Huai and Duroc pigs, respectively. The DEGs of the different pig breeds and muscle regions were significantly enriched in the biological processes and signaling pathways related to muscle fiber type, IMF deposition, lipid metabolism, PPAR signaling, cAMP signaling, amino acid metabolism, and ECM-receptor interaction. Our findings might help improve pork yield by using the obtained DEGs for marker-assisted selection and providing a theoretical reference for evaluating and improving pork quality.

Transcriptomic landscape reveals germline potential of porcine skin-derived multipotent dermal fibroblast progenitors

According to estimations, approximately about 15% of couples worldwide suffer from infertility, in which individuals with azoospermia or oocyte abnormalities cannot be treated with assisted reproductive technology. The skin-derived stem cells (SDSCs) differentiation into primordial germ cell-like cells (PGCLCs) is one of the major breakthroughs in the field of stem cells intervention for infertility treatment in recent years. However, the cellular origin of SDSCs and their dynamic changes in transcription profile during differentiation into PGCLCs in vitro remain largely undissected. Here, the results of single-cell RNA sequencing indicated that porcine SDSCs are mainly derived from multipotent dermal fibroblast progenitors (MDFPs), which are regulated by growth factors (EGF/bFGF). Importantly, porcine SDSCs exhibit pluripotency for differentiating into three germ layers and can effectively differentiate into PGCLCs through complex transcriptional regulation involving histone modification. Moreover, this study also highlights that porcine SDSC-derived PGCLCs specification exhibit conservation with the human primordial germ cells lineage and that its proliferation is mediated by the MAPK signaling pathway. Our findings provide substantial novel insights into the field of regenerative medicine in which stem cells differentiate into germ cells in vitro, as well as potential therapeutic effects in individuals with azoospermia and/or defective oocytes.

Combination of Transcriptomics and Proteomics Reveals Differentially Expressed Genes and Proteins in the Skin of EDAR Gene-Targeted and Wildtype Cashmere Goats

Cashmere goats play a pivotal role in the animal hair industry and are economically valuable. Cashmere is produced through the periodic growth of secondary hair follicles. To improve their yield of cashmere, the regulatory mechanisms of cashmere follicle growth and development need to be analysed. Therefore, in this study, EDAR gene-targeted cashmere goats were used as an animal model to observe the phenotypic characteristics of abnormal hair growth and development at the top of the head. Transcriptomic and proteomic techniques were used to screen for differentially expressed genes and proteins. In total, 732 differentially expressed genes were identified, including 395 upregulated and 337 downregulated genes. In addition, 140 differentially expressed proteins were identified, including 69 upregulated and 71 downregulated proteins. These results provide a research target for elucidating the mechanism through which EDAR regulates hair follicle growth in cashmere goats. It also enriches the available data on the regulatory network involved in hair follicle growth.

Hair Follicle Development and Cashmere Traits in Albas Goat Kids

The objectives of this trial were to study the growth and development of hair follicles and cashmere traits in cashmere goats and to provide a theoretical basis for the regulation of secondary hair follicle development and the scientific breeding selection of cashmere goats. Twelve single-fetal female kids were selected as research objects. A long-term tracking plan was created to regularly determine their growth performance, cashmere performance, and hair follicle traits. The results showed no significant difference in live weight after the first and second combing. The cashmere yield and unit yield of the first combing were significantly higher than those of the second combing (p < 0.05). Sections of hair follicles showed that the primary hair follicles are almost fully developed by 1 month, and the secondary hair follicles are fully developed by 5-6 months after birth. The primary hair follicle density (PFD) and secondary hair follicle density (SFD) were highest at birth and decreased within 1 month; and SFD was stable at 5-6 months of age. The change of MSFD took a maximum time of 2 to 3 months. The S:P increase reached its peak at 6 months. BMP4 expression increased with time. FGF2, FGF21 and BMP7 were higher at 3 months old than at the other two-time points. In conclusion, this study determined the total development time of primary and secondary hair follicles from morphology and speculated that FGF2, FGF21, and BMP7 may play a regulatory role in developing secondary hair follicles. Therefore, the period from birth to 6 months of age was the best time to regulate secondary hair follicle development in cashmere goats kids. The traits of the hair follicle and cashmere at 6 months of age could be breeding selection indicators for cashmere goats.

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.

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.

Selection signatures for fiber production in commercial species: A review

Natural fibers derived from diverse animal species have gained increased attention in recent years due to their favorable environmental effects, long-term sustainability benefits, and remarkable physical and mechanical properties that make them valuable raw materials used for textile and non-textile production. Domestication and selective breeding for the economically significant fiber traits play an imperative role in shaping the genomes and, thus, positively impact the overall productivity of the various fiber-producing species. These selection pressures leave unique footprints on the genome due to alteration in the allelic frequencies at specific loci, characterizing selective sweeps. Recent advances in genomics have enabled the discovery of selection signatures across the genome using a variety of methods. The increased demand for 'green products' manufactured from natural fibers necessitates a detailed investigation of the genomes of the various fiber-producing plant and animal species to identify the candidate genes associated with important fiber attributes such as fiber diameter/fineness, color, length, and strength, among others. The objective of this review is to present a comprehensive overview of the concept of selection signature and selective sweeps, discuss the main methods used for its detection, and address the selection signature studies conducted so far in the diverse fiber-producing animal species.

Defining ovine dermal papilla cell markers and identifying key signaling pathways regulating its intrinsic properties

Dermal papilla cell (DPC), one of the key cell types during hair follicle development and regeneration, specifies hair size, shape and cycling. It is also an important in vitro screening model for hair growth. Although some characteristics of DPCs, such as agglutinative growth and marker genes, have been studied in mice and humans, the intrinsic properties of ovine DPCs and the regulatory mechanism of the intrinsic properties during continued culture in vitro remained unknown. In this study, based on our previous single-cell transcriptome sequencing on sheep lambskin, we verified SOX18 and PDGFRA as the novel marker genes of ovine DPCs through immunofluorescence staining on skin sections and cultured DPCs. Using continued cell culture and alkaline phosphatase staining, we found that different from mice and humans, ovine DPCs exhibit particularly robust and stable aggregation with unbated alkaline phosphatase activity till 30 passages during continued culture in vitro. Also, we found that the expression of some marker genes and the activity of Wnt/β-catenin signaling differ between early passaged DPCs and multiple passaged DPCs. Further, using Wnt/β-catenin agonist and antagonist, we demonstrated that Wnt/β-catenin signaling could regulate cell aggregation and alkaline phosphatase activity of ovine DPCs through regulating FGF and IGF signaling. This study provides the basis for isolating ovine DPCs and defines their intrinsic properties, which contribute to improving wool performance and medicine of hair regeneration.

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.

miR-381-3p Inhibits Intramuscular Fat Deposition through Targeting FABP3 by ceRNA Regulatory Network

Intramuscular fat (IMF) deposition is an important determinant of pork quality and a complex process facilitated by non-coding ceRNAs. In this study, 52 Berkshire × Anqing Sixwhite crossbred pigs were slaughtered to measure eight carcass and pork quality traits. Whole-transcriptome sequencing analysis was performed using longissimus dorsi samples of six low- and high-IMF samples; 34 ceRNA networks, based on 881, 394, 158 differentially expressed (DE) lncRNAs, miRNAs, and mRNAs, were constructed. Following weighted gene co-expression network analysis between the low and high IMF, only one ceRNA, lncRNA4789/miR-381-3p/FABP3, that showed similar DE trend in longissimus dorsi tissue was retained. Dual-luciferase reporter assays further indicated that FABP3 was a direct, functional target of miR-381-3p, where miR-381-3p overexpression inhibited the mRNA and protein expression of FABP3. In addition, overexpressed lncRNA4789 attenuated the effect of miR-381-3p on FABP3 by sponging miR-381-3p. Cell function verification experiment demonstrated that miR-381-3p suppressed IMF deposition by inhibiting preadipocyte cell differentiation and lipid droplet deposition via the suppression of FABP3 expression in the peroxisome proliferator-activated receptor signalling pathway, whereas lncRNA4789 rescued FABP3 expression by sponging miR-381-3p. Our study may aid in identifying novel molecular markers for its optimization in IMF which is of importance in breeding for improving pork quality.

Understanding Mammalian Hair Follicle Ecosystems by Single-Cell RNA Sequencing

Simple Summary

Single-cell sequencing technology can reflect cell population heterogeneity at the single-cell level, leading to a better understanding of the role of individual cells in the microenvironment. Over the past few years, single-cell sequencing technology has not only made more new discoveries in the study of cellular heterogeneity of other rare cells such as stem cells, but has also become the most powerful research method for embryonic development, organ differentiation, cancer occurrence, and cell mapping. In this review, we outline the use of scRNA-seq in hair follicles. In particular, by focusing on landmark studies and the recent discovery of novel subpopulations of hair follicles, we summarize the phenotypic diversity of hair follicle cells and their links to hair follicle morphogenesis. Enhancing our understanding of the progress of hair follicle research will help to elucidate the regulatory mechanisms that determine the fate of different types of cells in the hair follicle, thereby guiding hair loss treatment and hair-producing economic animal breeding research.

Abstract

Single-cell sequencing technology can fully reflect the heterogeneity of cell populations at the single cell level, making it possible for us to re-recognize various tissues and organs. At present, the sequencing study of hair follicles is transiting from the traditional ordinary transcriptome level to the single cell level, which will provide diverse insights into the function of hair follicle cells. This review focuses on research advances in the hair follicle microenvironment obtained from scRNA-seq studies of major cell types in hair follicle development, with a special emphasis on the discovery of new subpopulations of hair follicles by single-cell techniques. We also discuss the problems and current solutions in scRNA-seq observation and look forward to its prospects.

Comprehensive Transcriptome Analysis of Hair Follicle Morphogenesis Reveals That lncRNA-H19 Promotes Dermal Papilla Cell Proliferation through the Chi-miR-214-3p/β-Catenin Axis in Cashmere Goats

Cashmere is initiated and develops in the fetal stages and the number and density of secondary hair follicles (SHFs) determine cashmere production and quality. Growing evidence indicates that both microRNA (miRNA) and long non-coding RNA (lncRNA) play an indispensable role in hair follicle (HF) growth and development. However, little is known about miRNAs, lncRNAs, and their functions as well as their interactions during cashmere initiation and development. Here, based on lncRNA and miRNA high-throughput sequencing and bioinformatics analysis, we identified 10,485 lncRNAs, 40,639 mRNAs, and 605 miRNAs in cashmere goat skin during HF induction, organogenesis, and cytodifferentiation stages. Among them, 521 lncRNAs, 5976 genes, and 204 miRNAs were differentially expressed (DE). KEGG analysis of DE genes indicated that ECM–receptor interaction and biosynthesis of amino acids were crucial for HF development. Notch, TGF-beta, and Wnt signaling pathways were also identified, which are conventional pathways associated with HF growth and development. Then, the ceRNA regulatory network was constructed, and the impact of lncRNA H19 was investigated in dermal papilla (DP) cells. The MTT, CCK-8, and EdU assays showed that the viability and proliferation of DP cells were promoted by H19, and mechanistic studies suggested that H19 performed its function through the chi-miR-214-3p/β-catenin axis. The present study created a resource for lncRNA, miRNA, and mRNA studies in cashmere morphogenesis. It could contribute to a better understanding of the molecular mechanism of ncRNAs involved in the regulation of HF growth and development.

Neural, adipocyte and hepatic differentiation potential of primary and secondary hair follicle stem cells isolated from Arbas Cashmere goats

Background

Arbas Cashmere goats are excellent domestic breeds with high yields of wool and cashmere. Their wool and cashmere can bring huge benefits to the livestock industry. Our studies intend to more fully understand the biological characteristics of hair follicle stem cells (HFSCs) in order to further explore the mechanisms of wool and cashmere regular regeneration. And they have been increasingly considered as promising multipotent cells in regenerative medicine because of their capacity to self-renew and differentiate. However, many aspects of the specific growth characteristics and differentiation ability of HFSCs remain unknown. This study aimed to further explore the growth characteristics and pluripotency of primary hair follicle stem cells (PHFSCs) and secondary hair follicle stem cells (SHFCs).

Results

We obtained PHFSCs and SHFSCs from Arbas Cashmere goats using combined isolation and purification methods. The proliferation and vitality of the two types of HFSCs, as well as the growth patterns, were examined. HFSC-specific markers and genes related to pluripotency, were subsequently identified. The PHFSCs and SHFSCs of Arbas Cashmere goat have a typical cobblestone morphology. Moreover, the PHFSCs and SHFSCs express HFSC surface markers, including CD34, K14, K15, K19 and LGR5. We also identified pluripotency-associated gene expression, including SOX2, OCT4 and SOX9, in PHFSCs and SHFSCs. Finally, PHFSCs and SHFSCs displayed multipotent abilities. PHFSCs and SHFSCs can be directed to differentiate into adipocyte-like, neural-like, and hepatocyte-like cells.

Conclusions

In conclusion, this study confirmed that the biological characteristics and differentiation potential of PHFSCs and SHFSCs from Arbas Cashmere goats. These findings broaden and refine our knowledge of types and characteristics of adult stem cells.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-022-03420-3.

Single-cell transcriptomics reveals male germ cells and Sertoli cells developmental patterns in dairy goats

Spermatogenesis holds considerable promise for human-assisted reproduction and livestock breeding based on stem cells. It occurs in seminiferous tubules within the testis, which mainly comprise male germ cells and Sertoli cells. While the developmental progression of male germ cells and Sertoli cells has been widely reported in mice, much less is known in other large animal species, including dairy goats. In this study, we present the data of single cell RNA sequencing (scRNA-seq) for 25,373 cells from 45 (pre-puberty), 90 (puberty), and 180-day-old (post-puberty) dairy goat testes. We aimed to identify genes that are associated with key developmental events in male germ cells and Sertoli cells. We examined the development of spermatogenic cells and seminiferous tubules from 15, 30, 45, 60, 75, 90, 180, and 240-day-old buck goat testes. scRNA-seq clustering analysis of testicular cells from pre-puberty, puberty, and post-puberty goat testes revealed several cell types, including cell populations with characteristics of spermatogonia, early spermatocytes, spermatocytes, spermatids, Sertoli cells, Leydig cells, macrophages, and endothelial cells. We mapped the timeline for male germ cells development from spermatogonia to spermatids and identified gene signatures that define spermatogenic cell populations, such as AMH, SOHLH1, INHA, and ACTA2. Importantly, using immunofluorescence staining for different marker proteins (UCHL1, C-KIT, VASA, SOX9, AMH, and PCNA), we explored the proliferative activity and development of male germ cells and Sertoli cells. Moreover, we identified the expression patterns of potential key genes associated with the niche-related key pathways in male germ cells of dairy goats, including testosterone, retinoic acid, PDGF, FGF, and WNT pathways. In summary, our study systematically investigated the elaborate male germ cells and Sertoli cells developmental patterns in dairy goats that have so far remained largely unknown. This information represents a valuable resource for the establishment of goat male reproductive stem cells lines, induction of germ cell differentiation in vitro, and the exploration of sequential cell fate transition for spermatogenesis and testicular development at single-cell resolution.

Identification of Genes Related to Hair Follicle Cycle Development in Inner Mongolia Cashmere Goat by WGCNA

Cashmere goat from Inner Mongolia is an excellent local breed in China, and the related cashmere product is a kind of precious textile raw material with high price. Cashmere is generated from secondary hair follicles, which has obvious annual periodicity and includes three different stages: anagen, catagen, and telogen. Therefore, we investigated skin transcriptome data for 12 months using weighted gene co-expression network analysis (WGCNA) to explore essential modules, pathways, and genes responsible for the periodic growth and development of secondary hair follicles. A total of 17 co-expression modules were discovered by WGCNA, and there is a strong correlation between steelblue module and month (0.65, p = 3E−09), anagen (0.52, p = 1E−05), telogen (−0.6, p = 8E−08). Gene expression was generally high during late anagen to catagen (June to December), while expression was downregulated from telogen to early anagen (January–May), which is similar to the growth rule of hair follicle cycle. KEGG pathway enrichment analyses of the genes of steelblue module indicated that genes are mainly enriched in Cell cycle, Wnt signaling pathway, p53 signaling pathway and other important signal pathways. These genes were also significantly enriched in GO functional annotation of the cell cycle, microtubule movement, microtubule binding, tubulin binding, and so on. Ten genes (WIF1, WNT11, BAMBI, FZD10, NKD1, LEF1, CCND3, E2F3, CDC6, and CDC25A) were selected from these modules, and further identified as candidate biomarkers to regulate periodic development of hair follicles using qRT-PCR. The Wnt signaling pathway and Cell cycle play an important role in the periodic development of hair follicles. Ten genes were identified as essential functional molecules related to periodic development of hair follicle. These findings laid a foundation for understanding molecular mechanisms in biological functions such as hair follicle development and hair growth in cashmere goats.

Identification of the Key Genes Associated with Different Hair Types in the Inner Mongolia Cashmere Goat

The Inner Mongolia cashmere goat is an excellent local breed in China. According to the characteristics of wool quilts, the Inner Mongolia cashmere goat can be divided into three types: a long-hair type (hair length of >22 cm), a short-hair type (hair length of ≤13 cm), and an intermediate type (hair length of >13 cm and ≤22 cm). It is found that hair length has a certain reference value for the indirect selection of other important economic traits of cashmere. In order to explore the molecular mechanisms and related regulatory genes of the different hair types, a weighted gene coexpression network analysis (WGCNA) was carried out on the gene expression data and phenotypic data of 12-month-old Inner Mongolia cashmere goats with a long-hair type (LHG) and a short-hair type (SHG) to explore the coexpression modules related to different coat types and nine candidate genes, and detect the relative expression of key candidate genes. The results showed that the WGCNA divided these genes into 19 coexpression modules and found that there was a strong correlation between one module and different hair types. The expression trends of this module’s genes were different in the two hair types, with high expression in the LHG and low expression in the SHG. GO functions are mainly concentrated in cellular components, including intermediate filaments (GO:0005882), intermediate filament cytoskeletons (GO:0045111), and cytoskeletal parts (GO:0044430). The KEGG pathway is mainly enriched in arginine as well as proline metabolism (chx00330) and the MAPK signaling pathway (chx04010). The candidate genes of the different hair types, including the KRT39, KRT74, LOC100861184, LOC102177231, LOC102178767, LOC102179881, LOC106503203, LOC108638293, and LOC108638298 genes, were screened. Through qRT-PCR, it was found that there were significant differences in these candidate genes between the two hair types, and most of them had a significant positive correlation with hair length. It was preliminarily inferred that these candidate genes could regulate the different hair types of cashmere goats and provide molecular markers for hair growth.

Single-cell sequencing reveals the new existence form of dermal papilla cells in the hair follicle regeneration of cashmere goats

The problem of human hair loss has caused widespread concern, however, such research is difficult because the periodicity is not obvious and the deeper levels knowledge of dermal papilla (DP) stem cells' differentiation are limited. Here, cashmere goats which have obvious periodicity of hair follicles were used, based on unbiased scRNA sequencing, we constructed DP cell lineage differentiation trajectory and revealed the key genes, signals and functions involved in cell fate decisions. And then we revealed the molecular landscape of hair follicle on regeneration. Revealed that DP cells differentiate into four intermediate cell states at different periodicity: Intermediate-cell-10 showed important functions in the growth and maintenance of cashmere; intermediate-cell-1 acting on apoptosis and cashmere shedding; intermediate-cell-0 initiated new follicular cycles, the migration of hair follicles and the occurrence of cashmere; and intermediate-cell-15 are suggested to be DP progenitor cells. In general, we provide new insights for hair regrowth. At the same time, it provides a new research ideas, directions and molecular landscape for the mechanism of dermal papilla cells.

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.