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

Studies in humans and mice have revealed that hair follicle morphogenesis relies on tightly coordinated ectodermal–mesodermal interactions, involving multiple signals and regulatory factors. DNA methylation and long non-coding RNA (lncRNA) play a critical role in early embryonic skin development by controlling gene expression. Acting as an indirect regulator, lncRNA could recruit DNA methyltransferases to specific genomic sites to methylate DNA. However, the molecular regulation mechanisms underlying hair follicle morphogenesis is unclear in cashmere goat. In this study, RNA-seq and whole-genome bisulfite sequencing (WGBS) in embryonic day 65 (E 65) and E 120 skin tissues of cashmere goat were used to reveal this complex regulatory process. The RNA-seq, qRT-PCR, and immunohistochemistry results showed that Wnt signaling played an important role in both hair follicle induction and differentiation stage; transcriptional factors (TFs), including HOXC13, SOX9, SOX21, JUNB, LHX2, VDR, and GATA3, participated in hair follicle differentiation via specific expression at E 120. Subsequently, the combination of WGBS and RNA-seq analysis showed that the expression of some hair follicle differentiation genes and TF genes were negatively correlated with the DNA methylation level generally. A portion of hair follicle differentiation genes were methylated and repressed in the hair follicle induction stage but were subsequently demethylated and expressed during the hair follicle differentiation stage, suggesting that DNA methylation plays an important role in hair morphogenesis by regulating associated gene expression. Furthermore, 45 upregulated and 147 downregulated lncRNAs in E 120 compared with E 65 were identified by lncRNA mapping, and then the potential differentially expressed lncRNAs associated with DNA methylation on the target gene were revealed. In conclusion, critical signals and genes were revealed during hair follicle morphogenesis in the cashmere goat. In this process, DNA methylation was lower in the hair follicle differentiation compared with the hair follicle induction stage and may play an important role in hair morphogenesis by regulating associated gene expression. Furthermore, potential lncRNAs associated with DNA methylation on target genes were delineated. This study enriches the regulatory network and molecular mechanisms on hair morphogenesis.

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.

The plectin (PLEC) gene is crucial in regulating muscle development and maintaining the cytoskeleton. An abnormal expression of PLEC can lead to muscle atrophy and muscular dystrophy. In a previous study, we found that Leizhou black goats exhibit abundant structural variations in the PLEC gene. However, the genetic effects of these variations on growth traits and meat quality in goats are not fully understood. In this study, three PLEC copy number variations (CNVs) were identified in a population of 417 Leizhou black goats, using quantitative polymerase chain reaction (qPCR) technology. Population distribution analysis revealed a high abundance of various types of these three CNVs. PLEC mRNA was found to be highly expressed in muscle tissue and remained consistently high from 1 month to 24 months after birth. Specifically, the gain type of CNV-1 (chr14: 81056401-81064800) showed a significant association with PLEC mRNA expression in muscle tissue (p < 0.01). The sequence of CNV-1 in PLEC shares similarities with three domain superfamilies associated with muscle development and skin disease. Furthermore, there were significant differences in chest circumference, body weight, carcass weight, the cross-sectional area of the longissimus dorsi lumbar muscle, and shear stress between different types of CNV-1 (p < 0.05). Notably, goats with the CNV-1 gain type demonstrated better phenotypic values compared to those with loss and normal types. These findings suggest that PLEC CNV-1 could play a crucial role in the growth and muscle development of Leizhou black goats, making it a potential marker for assisted selection in goat breeding.

In recent years, the poultry industry has embraced genome-wide assisted breeding techniques for broiler production, focusing primarily on single nucleotide polymorphisms (SNPs). However, the potential of genomic structural variations (SVs) in broiler traits remains largely unexplored due to cost and technology limitations. We developed a new chicken 21 K SV genotyping array and assessed its performance on 303 Wenchang chickens. Our evaluation showed a 99.20 % call rate with 57.13 % polymorphism across all microarray loci that identified SVs within the Wenchang chicken population. Genome-wide association studies (GWAS) of these SVs identified five significant variants associated with economically important traits, along with growth and development-related candidate genes including PLXAN4, ACBD3, and RNASET2. These findings demonstrate the value of SV-based approaches for improving genetic selection methods and enhancing productivity in poultry breeding programs. Our results show that high-quality 21 K SV genotyping arrays have significant potential for poultry breeding applications.