Synchronous profiling and analysis of mRNAs and ncRNAs in the dermal papilla cells from cashmere goats

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

CRABP1 Enhances the Proliferation of the Dermal Papilla Cells of Hu Sheep through the Wnt/β-catenin Pathway

BACKGROUND

The homologous proteins identified as cellular retinoic acid-binding proteins I and II (CRABP-I and CRABP-II) belong to a subset of intracellular proteins characterized by their robust affinity for retinoic acid, which plays an indispensable role in the development of hair follicle, including differentiation, proliferation, and apoptosis in keratinocytes. Previous research on Hu sheep hair follicles revealed the specific expression CRABP1 in dermal papilla cells (DPCs), suggesting that CRABP1 has a potential role in regulating the DPC population. Therefore, the main purpose of this study is to expose the performance of the CRABP1 genes in the development and proliferation of DPCs.

METHODS

Initially, overexpression and inhibition of CRABP1 in the DPCs were conducted through overexpression vector and siRNA. CCK-8, EDU, and RT-PCR cell cycle assays and immunostaining were performed to evaluate the proliferation and cell cycle of dermal papilla cells (DPCs). Although, the influence of CRABP1 upon β-catenin in dermal papilla cells (DPCs) was found using immunofluorescence labeling. Finally, RT-PCR was conducted to assess the impact of CRABP1 on the expression levels of CTNNB1, TCF4, and LEF1 in DPCs involved in the Wnt/β-catenin signaling pathway.

RESULTS

The results showed that CRABP1 overexpression promotes the growth rates of DPCs and significantly enhances the proportion of S-phase cells compared with the control group (p < 0.05). The results were the opposite when CRABP1 was a knockdown. In contrast, there was a significant decline in the mRNA expression levels of CTNNβ1, LEF1 (p < 0.05), and TCF4 (p < 0.01) by CRABP1 knockdown.

CONCLUSIONS

This study found that CRABP1 influences the expression of important genes within the Wnt/β-catenin signaling pathway and promotes DPC proliferation. This investigation provides a theoretical framework to explain the mechanisms that control hair follicle morphogenesis and development.

Genetic insights into fiber quality, coat color and adaptation in Changthangi and Muzzafarnagri sheep: A comparative skin transcriptome analysis

Changthangi sheep, which inhabit the high-altitude regions of Ladakh, are known for their fine fiber production and are characterized by grey skin and either black or white coats. In contrast, Muzzafarnagri sheep from the plains of Uttar Pradesh produce coarse wool and have white skin and coats. We conducted comparative global gene expression profiling on four biological replicates of skin from each breed. Notably, our analysis identified 149 up-regulated genes and 2,139 down-regulated genes in Changthangi sheep compared to Muzzafarnagri sheep, with a p-adjusted value (padj) of ≤0.05 and a Log2 fold change of ≥1.5. Gene Ontology analysis of the up-regulated genes revealed an enrichment of terms related to melanin biosynthesis and developmental pigmentation. Additionally, enriched KEGG pathways included tyrosine metabolism and metabolic pathways. Among the melanogenesis-related genes that exhibited higher expression in Changthangi sheep were TYR, TYRP1, DCT, SLC45A2, PMEL, MLANA, and OCA2. These findings confirm melanin's role in both the animals' black coat color and UV protection at high-altitude. Furthermore, we observed more pronounced expression of genes related to fiber quality, namely KRTAP6, KRTAP7, KRTAP13, and KRTAP2, in the fine wool-producing sheep from Ladakh. The results of the RNA sequencing were validated using real-time PCR on 10 genes governing fiber quality and coat color, with ACTB and PPIB serving as reference genes. In conclusion, our comparative skin transcriptome analysis of Changthangi and Muzzafarnagri sheep sheds light on the genetic differences associated with distinct phenotypic traits and environmental adaptability, offering valuable insights into the underlying mechanisms.

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.

Nanotechnology-based techniques for hair follicle regeneration

The hair follicle (HF) is a multicellular complex structure of the skin that contains a reservoir of multipotent stem cells. Traditional hair repair methods such as drug therapies, hair transplantation, and stem cell therapy have limitations. Advances in nanotechnology offer new approaches for HF regeneration, including controlled drug release and HF-specific targeting. Until recently, embryogenesis was thought to be the only mechanism for forming hair follicles. However, in recent years, the phenomenon of wound-induced hair neogenesis (WIHN) or de novo HF regeneration has gained attention as it can occur under certain conditions in wound beds. This review covers HF-specific targeting strategies, with particular emphasis on currently used nanotechnology-based strategies for both hair loss-related diseases and HF regeneration. HF regeneration is discussed in several modalities: modulation of the hair cycle, stimulation of progenitor cells and signaling pathways, tissue engineering, WIHN, and gene therapy. The HF has been identified as an ideal target for nanotechnology-based strategies for hair regeneration. However, some regulatory challenges may delay the development of HF regeneration nanotechnology based-strategies, which will be lastly discussed.

Transcriptomic Analysis Reveals Candidate Ligand-Receptor Pairs and Signaling Networks Mediating Intercellular Communication between Hair Matrix Cells and Dermal Papilla Cells from Cashmere Goats

Hair fiber growth is determined by the spatiotemporally controlled proliferation, differentiation, and apoptosis of hair matrix cells (HMCs) inside the hair follicle (HF); however, dermal papilla cells (DPCs), the cell population surrounded by HMCs, manipulate the above processes via intercellular crosstalk with HMCs. Therefore, exploring how the mutual commutations between the cells are molecularly achieved is vital to understanding the mechanisms underlying hair growth. Here, based on our previous successes in cultivating HMCs and DPCs from cashmere goats, we combined a series of techniques, including in vitro cell coculture, transcriptome sequencing, and bioinformatic analysis, to uncover ligand-receptor pairs and signaling networks mediating intercellular crosstalk. Firstly, we found that direct cellular interaction significantly alters cell cycle distribution patterns and changes the gene expression profiles of both cells at the global level. Next, we constructed the networks of ligand-receptor pairs mediating intercellular autocrine or paracrine crosstalk between the cells. A few pairs, such as LEP-LEPR, IL6-EGFR, RSPO1-LRP6, and ADM-CALCRL, are found to have known or potential roles in hair growth by acting as bridges linking cells. Further, we inferred the signaling axis connecting the cells from transcriptomic data with the advantage of CCCExplorer. Certain pathways, including INHBA-ACVR2A/ACVR2B-ACVR1/ACVR1B-SMAD3, were predicted as the axis mediating the promotive effect of INHBA on hair growth via paracrine crosstalk between DPCs and HMCs. Finally, we verified that LEP-LEPR and IL1A-IL1R1 are pivotal ligand-receptor pairs involved in autocrine and paracrine communication of DPCs and HMCs to DPCs, respectively. Our study provides a comprehensive landscape of intercellular crosstalk between key cell types inside HF at the molecular level, which is helpful for an in-depth understanding of the mechanisms related to hair growth.

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.

miR-29a-5p Inhibits Prenatal Hair Placode Formation Through Targeting EDAR by ceRNA Regulatory Network

Hair placode formation is an important stage of hair follicle morphogenesis and it is a complex process facilitated by non-coding RNAs. In this study, we conducted whole transcriptome sequencing analysis of skin, heart, liver, lung, and kidney tissues of day 41 (E41) normal and hairless pig embryos, and respectively detected 15, 8, and 515 skin-specific differentially expressed (DE) lncRNAs, miRNAs, and mRNAs. Furthermore, 18 competing endogenous RNA (ceRNA) networks were constructed. Following weighted gene co-expression network analysis (WGCNA) of stages E39, E41, E45, E52, and E60, between normal and hairless pig embryos, only two ceRNAs (lncRNA2162.1/miR-29a-5p/BMPR1b and lncRNA627.1/miR-29a-5p/EDAR) that showed period-specific differential expression in E41 skin were retained. Dual-luciferase reporter assays further indicated that EDAR was a direct, functioning target of miR-29a-5p and that no binding site was found in BMPR1b. Moreover, miR-29a-5p overexpression inhibited the mRNA and protein expression of EDAR while no significant differential expression of BMPR1b was detected. In addition, over-expressed lncRNA627.1 reduces the expression of miR-29a-5p and increase EDAR expression while inhibits lncRNA627.1 resulted in a opposite expression trend. Cell proliferation result demonstrated that lower expression of EDAR and lncRNA627.1 inhibited hair placode precursor cells (HPPCs) proliferation in a manner similar to that shown by over-expressed miR-29a-5p. This study identified that miR-29a-5p inhibited HPPCs proliferation via the suppression of EDAR expression in the EDA/EDAR signaling pathway, while lncRNA627.1 rescues EDAR expression. Our study provides a basis for a better understanding of the mechanisms underlying the ceRNA complex, miR29a-5p/EDAR/lncRNA627.1, that could regulate hair placode formation, which may help decipher diseases affecting human hair.

Integrated analysis of lncRNAs and mRNAs by RNA-Seq in secondary hair follicle development and cycling (anagen, catagen and telogen) of Jiangnan cashmere goat (Capra hircus)

BACKGROUND

Among the world's finest natural fiber composites is derived from the secondary hair follicles (SHFs) of cashmere goats yield one of the world's best natural fibres. Their development and cycling are characterized by photoperiodism with diverse, well-orchestrated stimulatory and inhibitory signals. Long non-coding RNA (lncRNAs) and mRNAs play important roles in hair follicle (HF) development. However, not many studies have explored their specific functions in cashmere development and cycling. This study detected mRNAs and lncRNAs with their candidate genes and related pathways in SHF development and cycling of cashmere goat. We utilized RNA sequencing (RNA-Seq) and bioinformatics analysis on lncRNA and mRNA expressions in goat hair follicles to discover candidate genes and metabolic pathways that could affect development and cycling (anagen, catagen, and telogen).

RESULTS

We identified 228 differentially expressed (DE) mRNAs and 256 DE lncRNA. For mRNAs, catagen and anagen had 16 upregulated and 35 downregulated DEGs, catagen and telogen had 18 upregulated and 9 downregulated DEGs and telogen and anagen had 52 upregulated and 98 downregulated DEGs. LncRNA witnessed 22 upregulated and 39 downregulated DEGs for catagen and anagen, 36 upregulated and 29 downregulated DEGs for catagen and telogen as well as 66 upregulated and 97 downregulated DEGs for telogen and anagen. Several key genes, including MSTRG.5451.2, MSTRG.45465.3, MSTRG.11609.2, CHST1, SH3BP4, CDKN1A, GAREM1, GSK-3β, DEFB103A KRTAP9-2, YAP1, S100A7A, FA2H, LOC102190037, LOC102179090, LOC102173866, KRT2, KRT39, FAM167A, FAT4 and EGFL6 were shown to be potentially important in hair follicle development and cycling. They were related to, WNT/β-catenin, mTORC1, ERK/MAPK, Hedgehog, TGFβ, NFkB/p38MAPK, caspase-1, and interleukin (IL)-1a signaling pathways.

CONCLUSION

This work adds to existing understanding of the regulation of HF development and cycling in cashmere goats via lncRNAs and mRNAs. It also serves as theoretical foundation for future SHF research in cashmere goats.

Single-Cell Sequencing Reveals Differential Cell Types in Skin Tissues of Liaoning Cashmere Goats and Key Genes Related Potentially to the Fineness of Cashmere Fiber

Cashmere fineness is one of the important factors determining cashmere quality; however, our understanding of the regulation of cashmere fineness at the cellular level is limited. Here, we used single-cell RNA sequencing and computational models to identify 13 skin cell types in Liaoning cashmere goats. We also analyzed the molecular changes in the development process by cell trajectory analysis and revealed the maturation process in the gene expression profile in Liaoning cashmere goats. Weighted gene co-expression network analysis explored hub genes in cell clusters related to cashmere formation. Secondary hair follicle dermal papilla cells (SDPCs) play an important role in the growth and density of cashmere. ACTA2, a marker gene of SDPCs, was selected for immunofluorescence (IF) and Western blot (WB) verification. Our results indicate that ACTA2 is mainly expressed in SDPCs, and WB results show different expression levels. COL1A1 is a highly expressed gene in SDPCs, which was verified by IF and WB. We then selected CXCL8 of SDPCs to verify and prove the differential expression in the coarse and fine types of Liaoning cashmere goats. Therefore, the CXCL8 gene may regulate cashmere fineness. These genes may be involved in regulating the fineness of cashmere in goat SDPCs; our research provides new insights into the mechanism of cashmere growth and fineness regulation by cells.

The Roles of Non-coding RNA in the Development and Regeneration of Hair Follicles: Current Status and Further Perspectives

Alopecia is a common problem that affects almost every age group and is considered to be an issue for cosmetic or psychiatric reasons. The loss of hair follicles (HFs) and hair caused by alopecia impairs self-esteem, thermoregulation, tactile sensation and protection from ultraviolet light. One strategy to solve this problem is HF regeneration. Many signalling pathways and molecules participate in the morphology and regeneration of HF, such as Wnt/β-catenin, Sonic hedgehog, bone morphogenetic protein and Notch. Non-coding RNAs (ncRNAs), especially microRNAs and long ncRNAs, have significant modulatory roles in HF development and regeneration via regulation of these signalling pathways. This review provides a comprehensive overview of the status and future prospects of ncRNAs in HF regeneration and could prompt novel ncRNA-based therapeutic strategies.

Cx43 promotes SHF-DPCs proliferation in the hair follicle of Albas cashmere goats from anagen to telogen

Connexin 43 (Cx43), known to form gap junction transmembrane channels between the cytoplasm of two adjacent cells, plays a key role in physiological functions, such as regulating cell growth, differentiation, and maintaining tissue homeostasis. Cashmere goat is an important farm animal that provides cashmere, which was produced by secondary hair follicles (SHF), for human consumption; however, there is no report about the role of Cx43 on the growth and development of SHF in cashmere goat. In this study, we investigated the effect of Cx43 on proliferation secondary hair follicle dermal papilla cells (SHF-DPCs) in Albas cashmere goat. In SHF-DPCs, Cx43 overexpression promoted cell proliferation and upregulated the expression of IGF-1, whereas Cx43 knockdown was associated with the opposite effects. These results suggested that Cx43 may promote cell proliferation by inducing IGF-1. Overall, our research not only contributes to a better understanding of the mechanism of the growth and development of SHF in cashmere goat, but also shed light on cashmere quality control in the future.

Cultivation of Hair Matrix Cells from Cashmere Goat Skins and Exemplified Applications

Simple Summary

A large scale of sequencing data pertaining to cashmere growth on cashmere goats have not been cost-effectively used due to the lack of in vitro cellular models, especially for hair matrix cells (HMCs)—the precursors of hair-forming keratinocytes, causing an enormous waste of data resources. Herein, we successfully isolated and cultivated previously unreported HMCs from cashmere goat skins and identified them morphologically and molecularly via their distinct appearance and signature genes’ expression from spatially adjacent dermal papilla cells. Through monitoring the effects of calcium and all-trans retinoic acid on HMCs using various biological techniques, we displayed that the cells are useful models to explore unsolved issues in hair fiber growth on goats. Therefore, our present success paves the road for further utilizing currently deposited data to unveil the secrets of cashmere growth and, ultimately, improve the quantity and quality of animal fibers.

Abstract

A functional interpretation of filtered candidates and predicted regulatory pathways related to cashmere growth from sequencing trials needs available cell models, especially for hair matrix cells (HMCs), whose continual proliferation and differentiation result in rapid hair growth. To fulfill such goals, we herein obtained primary goat HMCs via a microdissection-based method; optimized the selection of the culture medium and coating substances for better cell maintenance; and exemplified their usefulness through examining the effects of calcium and all-trans retinoic acid (ATRA) on cells using immunoblotting, flow cytometry, and other techniques. As a result, we successfully acquired primary and passaged goat HMCs with typical keratinocyte morphology. Calcium-free RPMI (Roswell Park Memorial Institute) 1640 and MEM (minimum Eagle’s medium) outperformed normal DMEM/F12 (Dulbecco’s modified Eagle’s medium/Nutrient Mixture F-12) on long-term cell maintenance, whereas serum-free media K-SFM and EpiLife failed to support cell growth. HMCs differed molecularly and morphologically from their neighbor dermal papilla cells on expressions of feature genes, such as HOXC13, and on characteristic keratinocyte-like appearances versus fibroblast shapes, respectively. Higher calcium concentrations significantly stimulated the expression of the genes (e.g., KRT1 and IVL) involved in keratinocyte differentiation and, promoted cell proliferation. Moreover, 10⁻⁵ M ATRA obviously boosted goat HMC expansions and changed their cell cycle distributions compared to the controls. Our study shines a light on researches exploring the mechanisms underlying the growth of cashmere.

Shedding light on cashmere goat hair follicle biology: from morphology analyses to transcriptomic landascape

BACKGROUND

Cashmere goat is known for its precious undercoat. Being photoperiod-dictated, cashmere growth has been studied focusing mainly on hair follicle cycle phases (anagen, catagen and telogen). An accurate molecular knowledge of the goat hair follicle cycle, disentangling gene expression changes during phases and recognizing timing boundaries, could be useful to improve cashmere goat management and ultimately cashmere production.

RESULTS

To better describe goat's hair follicle transcriptome we applied RNA-sequencing to isolated hair follicles from five Italian cashmere goats, during the anagen and catagen phase, identifying total of 214 differentially expressed genes (DEGs): 97 were up-regulated while 117 were down-regulated in catagen with respect to anagen. Gene Ontology and pathway analysis were performed. We detected 144 significant pathways spanning from estrogen, pluripotency of stem cells, thermogenesis and fatty acid metabolism that were strongly expressed during the hair follicle phases analysed. Finally, we validated promising DEGs by RT-qPCR in the same set of samples as well as in hair follicles and entire skin biopsies of another cashmere goats cohort accounting for early anagen, anagen, early catagen, and catagen phases.

CONCLUSIONS

As in the isolated hair follicles, some target genes were homogenously modulated during the four hair follicle phases. Ceruloplasmin (CP) and Keratin 4 (K4), confirmed their clear cut expression between growing and resting phase. In fact, K4 was almost absent in catagen phases while CP was barely expressed in anagen phases. In particular, the strong expression of K4 in early anagen makes it an eligible marker to track the beginning of a new hair cycle, and therefore defining the optimum time for cashmere harvesting.

Thymosin β4 Identified by Transcriptomic Analysis from HF Anagen to Telogen Promotes Proliferation of SHF-DPCs in Albas Cashmere Goat

Increasing cashmere yield is one of the important goals of cashmere goat breeding. To achieve this goal, we screened the key genes that can improve cashmere performance. In this study, we used the RNA raw datasets of the skin and dermal papilla cells of secondary hair follicle (SHF-DPCs) samples of hair follicle (HF) anagen and telogen of Albas cashmere goats and identified a set of significant differentially expressed genes (DEGs). To explore potential associations between gene sets and SHF growth features and to identify candidate genes, we detected functional enrichment and constructed protein-protein interaction (PPI) networks. Through comprehensive analysis, we selected Thymosin β4 (Tβ4), Rho GTPase activating protein 6 (ARHGAP6), ADAM metallopeptidase with thrombospondin type 1 motif 15, (ADAMTS15), Chordin (CHRD), and SPARC (Osteonectin), cwcv and kazal-like domains proteoglycan 1 (SPOCK1) as candidate genes. Gene set enrichment analysis (GSEA) for these genes revealed Tβ4 and ARHGAP6 have a close association with the growth and development of SHF-DPCs. However, the expression of Tβ4 in the anagen was higher than that in the telogen, so we finally chose Tβ4 as the ultimate research object. Overexpressing Tβ4 promoted and silencing Tβ4 inhibited the proliferation of SHF-DPCs. These findings suggest that Tβ4 can promote the growth and development of SHF-DPCs and indicate that this molecule may be a valuable target for increasing cashmere production.

Comprehensive analysis of circRNAs from cashmere goat skin by next generation RNA sequencing (RNA-seq)

Circular RNA (circRNA) is endogenous non-coding RNA (ncRNA) with a covalently closed circular structure. It is mainly generated through RNA alternative splicing or back-splicing. CircRNA is known in the majority of eukaryotes and very stable. However, knowledge of the circRNA involved in regulating cashmere fineness is limited. Skin samples were collected from Liaoning cashmere goats (LCG) and Inner Mongolia cashmere goats (MCG) during the anagen period. For differentially expressed circRNAs, RNA sequencing was performed, and the analysis led to an identification of 17 up-regulated circRNAs and 15 down-regulated circRNAs in LCG compared with MCG skin samples. In order to find the differentially expressed circRNAs in LCG, we carried out qPCRs on 10 candidate circRNAs in coarse type skin of LCG (CT-LCG) and fine type skin of LCG (FT-LCG). Four circRNAs: ciRNA128, circRNA6854, circRNA4154 and circRNA3620 were confirmed to be significantly differential expression in LCG. Also, a regulatory network of circRNAs-miRNAs was bioinformatically deduced and may help to understand molecular mechanisms of potential circRNA involvement in regulating cashmere fineness.