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

Deoxyshikonin from Arnebiae Radix promotes hair growth by targeting the Wnt/β-catenin signaling pathway

BACKGROUND

Alopecia is a common skin condition with limited effective treatments. Arnebiae Radix (AR) is a traditional Asian herb used for various skin disorders. However, its specific components and the mechanisms underlying its hair growth-promoting effects remain elusive. Therefore, this study aimed to investigate the primary active components in AR that are responsible for hair growth as well as determine the molecular mechanisms responsible for treating alopecia.

METHODS

Alopecia areata mice were employed to assess the influence of AR extracts on hair growth. The active AR components were identified via High-Performance Liquid Chromatography (HPLC). Furthermore, network pharmacology and molecular docking were carried out to predict the key targets of the main active AR compound, which were validated by Surface Plasma Resonance (SPR) analysis. Moreover, the mechanism of action of the identified active AR compound on human dermal papilla cells (HDPCs) and alopecia areata mice was investigated to determine its effects on relevant signaling pathways.

RESULTS

It was found that Deoxyshikonin (Ds) was the active component in AR with hair growth-promoting potential. Furthermore, it was predicted that Ds targeted 112 alopecia-related targets, including biological processes (such as positive modulation of cell migration and protein phosphorylation) and pathways (including cell cycle- and Wnt signaling pathway-related genes pathways). Network analysis revealed that CCND1 and GSK3β were the hub targets of Ds when treating alopecia. Molecular docking showed a strong binding affinity between Ds and GSK3β, which was validated by SPR results. Moreover, Ds improved HDPC's proliferation ability and promoted hair regeneration in alopecia mice. Similarly, Ds increased VEGF and IGF-1, reduced TGF-β1 content and GSK3β expression, and enhanced the p-GSK3β and β-catenin expression in the Wnt/β-catenin pathway.

CONCLUSION

This study showed that Ds was the main active AR component with promising potential as a hair growth stimulant. Mechanistically, Ds primarily targets GSK3β to promote the Wnt/β-catenin signaling pathway. This suggests that Ds could be an innovative therapeutic candidate for promoting hair regeneration.

Melatonin's Role in Hair Follicle Growth and Development: A Cashmere Goat Perspective

Hair follicles, unique skin appendages, undergo cyclic phases (anagen, catagen, telogen) governed by melatonin and associated molecular pathways. Melatonin, synthesized in the pineal gland, skin, and gut, orchestrates these cycles through antioxidant activity and signaling cascades (e.g., Wnt, BMP). This review examines melatonin's biosynthesis across tissues, its regulation of cashmere growth patterns, and its interplay with non-coding RNAs and the gut-skin axis. Recent advances highlight melatonin's dual role in enhancing antioxidant capacity (via Keap1-Nrf2) and modulating gene expression (e.g., Wnt10b, CTNNB1) to promote hair follicle proliferation. By integrating multi-omics insights, we construct a molecular network of melatonin's regulatory mechanisms, offering strategies to improve cashmere yield and quality while advancing therapies for human alopecia.

Deciphering the molecular drivers for cashmere/pashmina fiber production in goats: a comprehensive review

Cashmere, also known as pashmina, is derived from the secondary hair follicles of Cashmere/Changthangi goats. Renowned as the world's most luxurious natural fiber, it holds significant economic value in the textile industry. This comprehensive review enhances our understanding of the complex biological processes governing cashmere/pashmina fiber development and quality, enabling advancements in selective breeding and fiber enhancement strategies. The review specifically examines the molecular determinants influencing fiber development, with an emphasis on keratins (KRTs) and keratin-associated proteins (KRTAPs). It also explores the roles of key molecular pathways, including Wnt, Notch, BMP, NF-kappa B, VEGF, cAMP, PI3K-Akt, ECM, cell adhesion, Hedgehog, MAPK, Ras, JAK-STAT, TGF-β, mTOR, melanogenesis, FoxO, Hippo, and Rap1 signaling. Understanding these intricate molecular cascades provides valuable insights into the mechanisms orchestrating hair follicle growth, further advancing the biology of this coveted natural fiber. Expanding multi-omics approaches will enhance breeding precision and deepen our understanding of molecular pathways influencing cashmere production. Future research should address critical gaps, such as the impact of environmental factors, epigenetic modifications, and functional studies of genetic variants. Collaboration among breeders, researchers, and policymakers is essential for translating genomic advancements into practical applications. Such efforts can promote sustainable practices, conserve biodiversity, and ensure the long-term viability of high-quality cashmere production. Aligning genetic insights with conservation strategies will support the sustainable growth of the cashmere industry while preserving its economic and ecological value.

Unlocking the genetic secrets of Dorper sheep: insights into wool shedding and hair follicle development

Dorper sheep is popular among farming enterprises with strong adaptability, disease resistance, and roughage tolerance, and an unique characteristic of natural shedding of wool. In a large number of observations on experimental sheep farms, it was found that the wool of some sheep still had not shed after May, thus manual shearing was required. Therefore, understanding the molecular mechanisms of normal hair follicles (HFs) development is crucial to revealing the improvement of sheep wool-related traits and mammalian skin-related traits. In this study, transcriptome analysis was performed on skin tissues of adult Dorper ewes in the shedding (S) and non-shedding (N) groups in September 2019, January 2020, and March 2020, respectively. The results identified 3,278 differentially expressed transcripts (DETs) in the three comparison groups within the S group, 720 DETs in the three comparison groups within the N group, and 1,342 DETs in the three comparison groups between the S-vs-N groups. Time-series expression analysis revealed 2 unique expression patterns in HF development, namely, elevated expression in the anagen phase (A pattern) and the telogen phase (T pattern). DETs with stage-specific expression had a significant presence in processes related to the hair cycle and skin development, and several classic signaling pathways involved in sheep HF development, such as Rap1, estrogen, PI3K-Akt, and MAPK, were detected. Combined analysis of DETs, time-series expression data, and weighted gene co-expression network analysis identified core genes and their transcripts influencing HF development, such as DBI, FZD3, KRT17, ZDHHC21, TMEM79, and HOXC13. Additionally, alternative splicing analysis predicted that the isoforms XM_004004383.4 and XM_012125926.3 of ZDHHC21 might play a crucial role in sheep HF development. This study is a valuable resource for explaining the morphology of normal growth and development of sheep HFs and the genetic foundation of mammalian skin-related traits. It also offers potential insights into factors influencing human hair advancement.

Transcriptome profiles of the skin associated with the color of the black and white coat of Angora goats

The coat color of mammals, determined by the distribution of melanin, particularly eumelanin and pheomelanin, reflects intricate genetic and molecular processes. However, our understanding of the relationship between coat color, gene expression, and polymorphisms in goats remains nascent. This study investigates transcriptomic differences between black and white Angora goats across three distinct hair growth phases. Skin tissue samples from both colored and white Angora goats were analyzed using mRNA expression profiling. Three skin samples were taken from each goat as biological duplicates at every stage of hair follicle growth (September, January, and March). In total, 36 samples were analyzed in this study, including samples from two Angora goat varieties, three developmental stages (three biological replicates), and two technical replicates for RNA sequencing. Significant differences in gene expression were observed between black and white goats at each growth phase, particularly in genes associated with the melanogenesis pathway. Specifically, several pigmentation genes were identified solely in black goats, indicating phase-specific and breed-specific regulation. Noteworthy genes, such as SLC2A1, STAR, and SLC7A5, exhibited differential expression patterns across growth phases in black goats, further highlighting the complexity of melanogenesis regulation. This is the first study to use mRNA expression profiling of skin tissues to analyze coat color differences between black and white coated Angora goats at the anagen, catagen, and telogen stages. The identification of phase-specific and black goat-specific pigmentation genes provides valuable insights into the complex mechanisms governing coat color formation.

The molecular anatomy of cashmere goat hair follicle during cytodifferentiation stage

BACKGROUND

Cashmere, named as "soft gold", derives from the secondary hair follicles (SHFs) of cashmere goat which is vital to Northwest China's economy. The cytodifferentiation stage (E120), mirroring the complete hair follicle (HF) structure of adult goats and marking a critical phase in SHF development. Therefore, this study aims to enhance the understanding of SHF development and its impact on fiber quality, informing breeding strategies.

RESULTS

From the scRNA-seq data analysis, the intricate processes and transcriptional dynamics of inner layer cell differentiation of HFs were unveiled in this study. we identified nine cell populations during cytodifferentiation and key structures such as the hair shaft and inner root sheath. And we discovered three main inner layer lineages and seven subpopulations, clarifying their roles in specialization and signaling. Pseudotime mapping analysis showed cell evolution from early stage to mature stages marked by unique gene expressions, and the intermediate stage on the differentiation of each lineage was revealed. The identification and spatial localization of specific transcription factors, such as GATA3, LEF1 and PRDM1, as well as keratin genes highlight regulatory pathways involved in HF development, which was further validated by immunofluorescence. These findings suggested the potential strategies to improve fiber quality, and the discovery of diverse cell types and their developmental molecular mechanisms, particularly in this species-specific context, offered a nuanced view of the regulatory mechanisms driving HF development in cashmere goats.

CONCLUSION

Overall, these findings provide a systematic molecular atlas of skin, defining three major branches and cell states of inner layer cells of HF, and determining how the branch-specific transcription factors, keratins, and signals coordinate HF morphogenesis during cytodifferentiation stage. This research not only advances skin tissue research in goats but also holds broader implications for the understanding of HF regeneration and development across various species.

Chromosome-level genome assembly of the cashmere goat

The goat, an early domesticated ruminant, is a reliable source of cashmere, meat and milk in global agricultural production. Despite this, the genome of cashmere-rich goats has yet to be characterized. Here, we assembled the nearly complete genome of a cashmere goat from a highly economically valuable Inner Mongolian Cashmere buck, utilizing a combination of PacBio HiFi, ONT ultra-long reads, and Hi-C technologies. The size of this genome is 2.76 Gb, with a contig N50 of 95.22 Mb. All assembled sequences were anchored onto 29 autosomes and both sex chromosomes, with only two gaps present on the X chromosome. We identified 1,333.29 Mb (48.26%) of repetitive sequences and predicted 22,480 protein-coding genes. Assembly quality assessment of the genome demonstrated that our assembled cashmere goat genome surpasses the continuity, completeness, and accuracy of other published goat genomes. Taken together, we provided the first cashmere goat assembly, bridging the gap in the genome of important economic breeds of domestic goats, and providing a valuable reference resource for goat genetics and genome research.

Genetic characterization of cashmere goat (Capra hircus) populations in Mongolia

Objective

Characterization studies of the phenotypic and genetic diversity of Mongolian goats are limited, despite several goat breeds being registered in the country. This study aimed to evaluate the phenotypic and genetic diversity of 14 cashmere goat populations in Mongolia, consisting largely of identified goat breeds.

Methods

Body weight, cashmere quality, and coat color were the phenotypic traits considered in this study. A linear model was used to fit body weight and cashmere traits, and least squares means (LSMs) were estimated for the region and location classes. Genetic diversity and structure were assessed using a goat 50K SNP array.

Results

The studied populations exhibited greater phenotypic diversity at the regional level. A very small overall differentiation index (Fst: 0.017) was revealed by Wright's Fst and a very small overall inbreeding index (F ROH1 :0.019) was revealed based on runs of homozygosity. Genetic clustering of populations by principal components showed large variances for the two goat populations of the Russian admixture (Gobi Gurvan Saikhan and Uuliin Bor), and smaller but differentiated clusters for the remaining populations. Similar results were observed in the admixture analysis, which identified populations with the highest (Govi Gurvan Saikhan and Uuliin Bor) and lowest (Tsagaan Ovoo Khar) exotic admixtures. A genomewide association study (GWAS) of body weight and cashmere traits identified a few significant variants on chromosomes 2, 4, 5, 9, and 15, with the strongest variant for cashmere yield on chromosome 4. The GWAS on coat color yielded nine significant variants, with the strongest variants located on chromosomes 6, 13, and 18 and potential associations with KIT, ASIP, and MC1R genes. These signals were also found in other studies on coat color and patterns in goats.

Conclusion

Mongolian cashmere goats showed relatively low genetic differentiation and low inbreeding levels, possibly caused by the traditional pastoral livestock management system and the practice of trading breeding bucks across provinces, along with a recent increase in the goat population. Further investigation of cashmere traits using larger samples and alternative methods may help identify the genes or genomic regions underlying cashmere quality in goats.

Stress adaptation in Tibetan cashmere goats is governed by inherent metabolic differences and manifested through variable cashmere phenotypes

Tibetan cashmere goats are not only served as a valuable model for studying adaptation to hypoxia and high-altitude conditions but also playing a pivotal role in bolstering local economies through the provision of premium quality cashmere yarn. In this study, we performed an integration and network analysis of metabolomic, transcriptomic and proteomic to elucidate the role of differentially expressed genes, important metabolites, and relevant cellular and metabolic pathways between the fine (average 12.04 ± 0.03 μm of mean fiber diameter) and coarse cashmere (average 14.88 ± 0.05 μm of mean fber diameter) producing by Tibetan cashmere goats. We identified a distinction of 56 and 71 differential metabolites (DMs) between the F and C cashmere groups under positive and negative ion modes, respectively. The KEGG pathway enrichment analysis of these DMs highlighted numerous pathways predominantly involved in amino acid and protein metabolism, as indicated by the finding that the most impactful pathway was the mammalian target of rapamycin (mTOR) signalling pathway. In the F group, we identified a distinctive metabolic profile where amino acid metabolites including serine, histidine, asparagine, glutamic acid, arginine, valine, aspartic acid, tyrosine, and methionine were upregulated, while lysine, isoleucine, glutamine, tryptophan, and threonine were downregulated. The regulatory network and gene co-expression network revealed crucial genes, metabolites, and metabolic pathways. The integrative omics analysis revealed a high enrichment of several pathways, notably encompassing protein digestion and absorption, sphingolipid signalling, and the synaptic vesicle cycle. Within the sphere of our integrative analysis, DNMT3B was identified as a paramount gene, intricately associated with significant proteins such as HMCN1, CPB2, GNG12, and LRP1. Our present study delineated the molecular underpinnings governing the variations in cashmere characteristics by conducting comprehensive analyses across metabolomic, transcriptomic, and proteomic dimensions. This research provided newly insights into the mechanisms regulating cashmere traits and facilitated the advancement of selective breeding programs aimed at cultivating high-quality superfine Tibetan cashmere goats.

The demographic history and adaptation of Canarian goat breeds to environmental conditions through the use of genome-wide SNP data

BACKGROUND

The presence of goats in the Canary Islands dates back to the late 1st millennium BC, which coincides with the colonization by the Amazigh settlers. However, the exact geographic origin of Canarian goats is uncertain since the Amazigh peoples were distributed over a wide spatial range. Nowadays, three Canarian breeds (Palmera, Majorera and Tinerfeña) are officially recognized, along with two distinct South and North Tinerfeña ecotypes, with the South Tinerfeña and Majorera goats thriving in arid and dry semi-desertic environments and the Palmera and North Tinerfeña goats are adapted to humid and temperate areas that are influenced by trade winds. Genotypes for 224 Canarian goats were generated using the Illumina Goat single nucleotide polymorphism (SNP)50 BeadChip. By merging these data with the genotypes from 1007 individuals of African and Southern European ancestry, our aim was to ascertain the geographic origin of the Canarian goats and identify genes associated with adaptation to diverse environmental conditions.

RESULTS

The diversity indices of the Canarian breeds align with most of those of the analyzed local breeds from Africa and Europe, except for the Palmera goats that showed lower levels of genetic variation. The Canarian breeds demonstrate a significant genetic differentiation compared to other populations, which indicates a history of prolonged geographic isolation. Moreover, the phylogenetic reconstruction indicated that the ancestry of the Canarian goats is fundamentally North African rather than West African. The ADMIXTURE and the TreeMix analyses showed no evidence of gene flow between Canarian goats and other continental breeds. The analysis of runs of homozygosity (ROH) identified 13 ROH islands while the window-based FST method detected 25 genomic regions under selection. Major signals of selection were found on Capra hircus (CHI) chromosomes 6, 7, and 10 using various comparisons and methods.

CONCLUSIONS

This genome-wide analysis sheds new light on the evolutionary history of the four breeds that inhabit the Canary Islands. Our findings suggest a North African origin of the Canarian goats. In addition, within the genomic regions highlighted by the ROH and FST approaches, several genes related to body size and heat tolerance were identified.

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.

Effect of Fibroblast Growth Factor 10 and an Interacting Non-Coding RNA on Secondary Hair Follicle Dermal Papilla Cells in Cashmere Goats' Follicle Development Assessed by Whole-Transcriptome Sequencing Technology

Cashmere, a keratinised product of secondary hair follicles (SHFs) in cashmere goats, holds an important place in international high-end textiles. However, research on the complex molecular and signal regulation during the development and growth of hair follicles (HFs), which is essential for the development of the cashmere industry, is limited. Moreover, increasing evidence indicates that non-coding RNAs (ncRNAs) participate in HF development. Herein, we systematically investigated a competing endogenous RNA (ceRNA) regulatory network mediated by circular RNAs (circRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs) in skin samples of cashmere goat embryos, using whole-transcriptome sequencing technology. We obtained 6468, 394, and 239 significantly differentially expressed mRNAs, circRNAs, and miRNAs, respectively. These identified RNAs were further used to construct a ceRNA regulatory network, mediated by circRNAs, for cashmere goats at a late stage of HF development. Among the molecular species identified, miR-184 and fibroblast growth factor (FGF) 10 exhibited competitive targeted interactions. In secondary HF dermal papilla cells (SHF-DPCs), miR-184 promotes proliferation, inhibits apoptosis, and alters the cell cycle via the competitive release of FGF10. This study reports that FGF10 and its interaction with ncRNAs significantly affect SHF-DPCs, providing a reference for research on the biology of HFs in cashmere goats and other mammals.

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.

Transcriptome analysis of pika heart tissue reveals mechanisms underlying the adaptation of a keystone species on the roof of the world

High-altitude environments impose intense stresses on living organisms and drive striking phenotypic and genetic adaptations, such as hypoxia resistance, cold tolerance, and increases in metabolic capacity and body mass. As one of the most successful and dominant mammals on the Qinghai-Tibetan Plateau (QHTP), the plateau pika (Ochotona curzoniae) has adapted to the extreme environments of the highest altitudes of this region and exhibits tolerance to cold and hypoxia, in contrast to closely related species that inhabit the peripheral alpine bush or forests. To explore the potential genetic mechanisms underlying the adaptation of O. curzoniae to a high-altitude environment, we sequenced the heart tissue transcriptomes of adult plateau pikas (comparing specimens from sites at two different altitudes) and Gansu pikas (O. cansus). Differential expression analysis and weighted gene co-expression network analysis (WGCNA) were used to identify differentially expressed genes (DEGs) and their primary functions. Key genes and pathways related to high-altitude adaptation were identified. In addition to the biological processes of signal transduction, energy metabolism and material transport, the identified plateau pika genes were mainly enriched in biological pathways such as the negative regulation of smooth muscle cell proliferation, the apoptosis signalling pathway, the cellular response to DNA damage stimulus, and ossification involved in bone maturation and heart development. Our results showed that the plateau pika has adapted to the extreme environments of the QHTP via protection against cardiomyopathy, tissue structure alterations and improvements in the blood circulation system and energy metabolism. These adaptations shed light on how pikas thrive on the roof of the world.

Lef1 and Dlx3 May Facilitate the Maturation of Secondary Hair Follicles in the Skin of Gansu Alpine Merino

Lymphatic enhancer factor 1 (Lef1) and distal-less homeobox 3 (Dlx3) are the transcription factors involved in regulating hair follicle development in mice, goats, and other animals. Their deletion can lead to hair follicle deficiency. In this study, hematoxylin−eosin staining (HE), real-time quantitative PCR (RT-qPCR), immunohistochemistry, and immunofluorescence were used to analyze the expression, location, and biological functions of Lef1 and Dlx3 in the lateral skin of Gansu Alpine Merino aged 1, 30, 60, and 90 days. The results revealed that the number of hair follicles decreased with age and was significantly higher at 1 day than in the other three age groups (p < 0.05). The mRNA levels of Lef1 and Dlx3 in the skin of 30-day old Gansu Alpine Merino were significantly higher than those in the other three age groups (p < 0.05). Protein expression of Lef1 and Dlx3 was lowest at 1 day (p < 0.05) and peaked at 60 days. Lef1 and Dlx3 exhibited a high density and strong positive expression in the dermal papillae; additionally, Dlx3 exhibited a high density and strong positive expression in the inner and outer root sheaths. Collectively, Lef1 and Dlx3 may facilitate the maturation of secondary hair follicles, which is mainly achieved through the dermal papillae and inner and outer root sheaths.