Screening the key genes of hair follicle growth cycle in Inner Mongolian Cashmere goat based on RNA sequencing

Genetic diversity analysis of Inner Mongolia cashmere goats (Erlangshan subtype) based on whole genome re-sequencing

BACKGROUND

Inner Mongolia cashmere goat (IMCG), renowned for its superior cashmere quality, is a Chinese indigenous goat breed that has been developed through natural and artificial selection over a long period. However, recently, the genetic resources of IMCGs have been significantly threatened by the introduction of cosmopolitan goat breeds and the absence of adequate breed protection systems.

RESULTS

In order to assess the conservation effectiveness of IMCGs and efficiently preserve and utilize the purebred germplasm resources, this study analyzed the genetic diversity, kinship, family structure, and inbreeding of IMCGs utilizing resequencing data from 225 randomly selected individuals analyzed using the Plink (v.1.90), GCTA (v.1.94.1), and R (v.4.2.1) software. A total of 12,700,178 high-quality SNPs were selected through quality control from 34,248,064 SNP sites obtained from 225 individuals. The average minor allele frequency (MAF), polymorphic information content (PIC), and Shannon information index (SHI) were 0.253, 0.284, and 0.530, respectively. The average observed heterozygosity (Ho) and the average expected heterozygosity (He) were 0.355 and 0.351, respectively. The analysis of the identity by state distance matrix and genomic relationship matrix has shown that most individuals' genetic distance and genetic relationship are far away, and the inbreeding coefficient is low. The family structure analysis identified 10 families among the 23 rams. A total of 14,109 runs of homozygosity (ROH) were identified in the 225 individuals, with an average ROH length of 1014.547 kb. The average inbreeding coefficient, calculated from ROH, was 0.026 for the overall population and 0.027 specifically among the 23 rams, indicating a low level of inbreeding within the conserved population.

CONCLUSIONS

The IMCGs exhibited moderate polymorphism and a low level of kinship with inbreeding occurring among a limited number of individuals. Simultaneously, it is necessary to prevent the loss of bloodline to guarantee the perpetuation of the IMCGs' germplasm resources.

Genetic diversity and population structure in five Inner Mongolia cashmere goat populations using whole-genome genotyping

OBJECTIVE

As a charismatic species, cashmere goats have rich genetic resources. In the Inner Mongolia Autonomous Region, there are three cashmere goat varieties named and approved by the state. These goats are renowned for their high cashmere production and superior cashmere quality. Therefore, it is vitally important to protect their genetic resources as they will serve as breeding material for developing new varieties in the future.

METHODS

Three breeds including Inner Mongolia cashmere goats (IMCG), Hanshan White cashmere goats (HS), and Ujimqin white cashmere goats (WZMQ) were studied. IMCG were of three types: Aerbas (AEBS), Erlangshan (ELS), and Alashan (ALS). Nine DNA samples were collected for each population, and they were genomically re-sequenced to obtain high-depth data. The genetic diversity parameters of each population were estimated to determine selection intensity. Principal component analysis, phylogenetic tree construction and genetic differentiation parameter estimation were performed to determine genetic relationships among populations.

RESULTS

Samples from the 45 individuals from the five goat populations were sequenced, and 30,601,671 raw single nucleotide polymorphisms (SNPs) obtained. Then, variant calling was conducted using the reference genome, and 17,214,526 SNPs were retained after quality control. Individual sequencing depth of individuals ranged from 21.13× to 46.18×, with an average of 28.5×. In the AEBS, locus polymorphism (79.28) and expected heterozygosity (0.2554) proportions were the lowest, and the homologous consistency ratio (0.1021) and average inbreeding coefficient (0.1348) were the highest, indicating that this population had strong selection intensity. Conversely, ALS and WZMQ selection intensity was relatively low. Genetic distance between HS and the other four populations was relatively high, and genetic exchange existed among the other four populations.

CONCLUSION

The Inner Mongolia cashmere goat (AEBS type) population has a relatively high selection intensity and a low genetic diversity. The IMCG (ALS type) and WZMQ populations had relatively low selection intensity and high genetic diversity. The genetic distance between HS and the other four populations was relatively high, with a moderate degree of differentiation. Overall, these genetic variations provide a solid foundation for resource identification of Inner Mongolia Autonomous Region cashmere goats in the future.

Genomic Inbreeding and Runs of Homozygosity Analysis of Cashmere Goat

Cashmere goats are valuable genetic resources which are famous worldwide for their high-quality fiber. Runs of homozygosity (ROHs) have been identified as an efficient tool to assess inbreeding level and identify related genes under selection. However, there is limited research on ROHs in cashmere goats. Therefore, we investigated the ROH pattern, assessed genomic inbreeding levels and examined the candidate genes associated with the cashmere trait using whole-genome resequencing data from 123 goats. Herein, the Inner Mongolia cashmere goat presented the lowest inbreeding coefficient of 0.0263. In total, we identified 57,224 ROHs. Seventy-four ROH islands containing 50 genes were detected. Certain identified genes were related to meat, fiber and milk production (FGF1, PTPRM, RERE, GRID2, RARA); fertility (BIRC6, ECE2, CDH23, PAK1); disease or cold resistance and adaptability (PDCD1LG2, SVIL, PRDM16, RFX4, SH3BP2); and body size and growth (TMEM63C, SYN3, SDC1, STRBP, SMG6). 135 consensus ROHs were identified, and we found candidate genes (FGF5, DVL3, NRAS, KIT) were associated with fiber length or color. These findings enhance our comprehension of inbreeding levels in cashmere goats and the genetic foundations of traits influenced by selective breeding. This research contributes significantly to the future breeding, reservation and use of cashmere goats and other goat breeds.

Screening of hair follicle telogen-associated circRNAs in sheep and construction of their ceRNA network

Sheep breeds with hair-shedding traits have many advantages over non-shedding sheep breeds, not only because of reduced shearing labor and feeding management costs but also because it reduces in vitro parasites and improves adaptability to summer heat stress. The wool of Dorper sheep naturally sheds in spring due to the periodic growth of hair follicles. CircRNAs primarily regulate the morphogenesis of hair follicles through the ceRNA mechanism. In this study, five 2-year-old Dorper ewes with extreme hair-shedding phenotype (S) and three Dorper ewes with non-shedding (N) phenotype were selected for subsequent analyses. For RNA extraction, skin tissues were collected on 27th September 2019 (S1, N1), 3rd January 2020 (S2, N2), and 17th March 2020 (S3, N3), which were then subjected to RNA-seq. RNA-seq technology revealed 20,185 novel circRNAs in the hair follicles of Dorper sheep. Among them, 1450 circRNAs were differentially expressed (DE). Clustering heatmap and expression pattern analyses were performed on DE circRNAs, which indicated 78 circRNAs with T pattern (Telogen, highly expressed in telogen), and the source genes for candidate circRNAs were further screened by functional enrichment analysis, which identified 13 crucial genes enriched in pathways associated with hair follicle development. Additionally, a ceRNA regulatory network comprising 4 circRNAs, 11 miRNAs, and 13 target genes was constructed. Overall, this study screened circRNAs that may be associated with the telogen phase of hair follicles in sheep, providing a relevant theoretical basis for wool shedding in sheep and for breeding Dorper sheep with automatic wool shedding.

Single-cell profiling reveals a potent role of quercetin in promoting hair regeneration

Hair loss affects millions of people at some time in their life, and safe and efficient treatments for hair loss are a significant unmet medical need. We report that topical delivery of quercetin (Que) stimulates resting hair follicles to grow with rapid follicular keratinocyte proliferation and replenishes perifollicular microvasculature in mice. We construct dynamic single-cell transcriptome landscape over the course of hair regrowth and find that Que treatment stimulates the differentiation trajectory in the hair follicles and induces an angiogenic signature in dermal endothelial cells by activating HIF-1α in endothelial cells. Skin administration of a HIF-1α agonist partially recapitulates the pro-angiogenesis and hair-growing effects of Que. Together, these findings provide a molecular understanding for the efficacy of Que in hair regrowth, which underscores the translational potential of targeting the hair follicle niche as a strategy for regenerative medicine, and suggest a route of pharmacological intervention that may promote hair regrowth.

Identification and Characterization of Circular RNAs (circRNAs) Using RNA-Seq in Two Breeds of Cashmere Goats

Circular RNA (circRNA) is a type of non-coding RNA generated from back-splicing the reactions of linear RNA. It plays an important role in various cellular and biological processes. However, there are few studies about the regulatory effect of circRNAs on cashmere fiber traits in cashmere goats. In this study, the expression profiles of circRNAs in skin tissue were compared between Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats, with a significant difference in cashmere fiber yield, cashmere fiber diameter, and cashmere fiber color, using RNA-seq. A total of 11,613 circRNAs were expressed in the caprine skin tissue, and their type, chromosomal distribution, and length distribution were characterized. A total of 115 up-regulated circRNAs and 146 down-regulated circRNAs in LC goats were screened compared to ZB goats. The authenticity of 10 differentially expressed circRNAs was validated by detecting their expression levels and the head-to-tail splice junction using RT-PCR and DNA sequencing, respectively. The parent genes of differentially expressed circRNA were mainly enriched in some Gene Ontology (GO) terms and pathways related to cashmere fiber traits, such as the canonical Wnt signaling pathway, which is involved in the regulation of cell promotion, stem cell proliferation, Wnt signaling pathway regulation, epithelial morphogenesis, MAPK signaling pathway, and cell adhesion molecules pathway. Eight differentially expressed circRNAs were further selected to construct a circRNA-miRNA network, and some miRNAs that were previously reported as related to fiber traits were found in the network. This study provides a deep understanding of the roles of circRNAs in the regulation of cashmere fiber traits in cashmere goats and the involvement of differential splicing in phenotypic expression according to breed and region.

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.

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.

Identification and Molecular Analysis of m6A-circRNAs from Cashmere Goat Reveal Their Integrated Regulatory Network and Putative Functions in Secondary Hair Follicle during Anagen Stage

Simple Summary

Cashmere is a natural, high-end textile material. It is derived from the secondary hair follicle (SHFs) tissue in the skin of cashmere goats. Previous studies have indicated that m⁶A modifications in circRNA molecules play important roles in a variety of cells through multiple mechanisms. However, little information is available on the expression profile and functional regulatory characteristics of m⁶A-modified circRNA (m⁶A-circRNA) in SHFs of cashmere goats. In this study, a total of 15 m⁶A-circRNAs were identified. Six of these m⁶A-circRNAs were revealed to have significantly higher expression in skin at anagen than at telogen. To gain insight into the potential regulatory mechanisms of the anagen up-regulated m⁶A-circRNAs, we constructed the regulatory networks along with related pathways in SHFs of cashmere goats. In addition, we found that the expression trends of four m⁶A-circRNAs in the SHFs during SHF cycles were highly similar to their host genes. However, the expression patterns of two m⁶A-circRNAs were inconsistent with the linear RNAs from their host genes in the SHFs of cashmere goats. These results will provide new insights to elucidate the biological functions and regulatory features of m⁶A-circRNA in SHF development and cashmere growth in goats.

Abstract

N⁶-methyladenosine (m⁶A) is the most abundant modification in linear RNA molecules. Over the last few years, interestingly, many circRNA molecules are also found to have extensive m⁶A modification sites with temporal and spatial specific expression patterns. To date, however, little information is available concerning the expression profiling and functional regulatory characteristics of m⁶A modified circRNAs (m⁶A-circRNAs) in secondary hair follicles (SHFs) of cashmere goats. In this study, a total of fifteen m⁶A-circRNAs were identified and characterized in the skin tissue of cashmere goats. Of these, six m⁶A-circRNAs were revealed to have significantly higher expression in skin at anagen compared with those at telogen. The constructed ceRNA network indicated a complicated regulatory relationship of the six anagen up-regulated m⁶A-circRNAs through miRNA mediated pathways. Several signaling pathways implicated in the physiological processes of hair follicles were enriched based on the potential regulatory genes of the six anagen up-regulated m⁶A-circRNAs, such as TGF-beta, axon guidance, ribosome, and stem cell pluripotency regulatory pathways, suggesting the analyzed m⁶A-circRNAs might be essentially involved in SHF development and cashmere growth in cashmere goats. Further, we showed that four m⁶A-circRNAs had highly similar expression trends to their host genes in SHFs of cashmere goats including m⁶A-circRNA-ZNF638, -TULP4, -DNAJB6, and -CAT. However, the expression patterns of two m⁶A-circRNAs (m⁶A-circRNA-STAM2 and -CAAP1) were inconsistent with the linear RNAs from their host genes in the SHFs of cashmere goats. These results provide novel information for eluci-dating the biological function and regulatory characteristics of the m⁶A-circRNAs in SHF development and cashmere growth in goats.

Five SNPs Within the FGF5 Gene Significantly Affect Both Wool Traits and Growth Performance in Fine-Wool Sheep (Ovis aries)

Fibroblast growth factor 5 (FGF5) gene, a member of fibroblast growth factor superfamily, plays significant roles in the regulation of the hair growth cycle during the development of mammalian hair follicles as well as the skeletal muscle development. In this study, DNA sequencing was used to scan the putative SNPs within the full-length of FGF5 gene, and SNPscan high-throughput technique was applied in the individual genotyping of 604 crossbred sheep. 10 SNPs were identified within FGF5 gene while five of them located in intron 1 could be genotyped, namely SNP1 (g. 105914953 G > A), SNP2 (g. 105922232 T > C), SNP3 (g. 105922244 A > G), SNP4 (g. 105922334 A > T) and SNP5 (g. 105922340 G > T). All these SNPs were in accord with the Hardy-Weinberg equilibrium (P > 0.05), and displayed the moderate polymorphism with PIC values ranging from 0.302 to 0.374. Thereafter, the correlation analysis between each SNP locus and economic traits including wool length, greasy wool weight and growth performance of sheep was systematically implemented. In our results, SNP1, SNP3, SNP4 and SNP5 were significantly associated with wool length, greasy wool weight and growth traits of SG sheep (P < 0.05); SNP1, SNP2, SNP3, and SNP4 were significantly correlated with wool length and growth traits of SSG sheep (P < 0.05). Meanwhile, our study revealed a strong linkage disequilibrium (LD) relationship among these SNPs (r² > 0.33), except for SNP3 and SNP4 sites (r² = 0.30). Combination genotype analysis showed that combination genotypes were significantly associated with mean fiber diameter of SG (P < 0.05), and body weight trait of SSG (P < 0.01). The above findings suggested that these SNP loci might affect economic traits synergistically and could be regarded as potential molecular markers for improving both wool production and growth performance of fine-wool sheep, which lay a molecular foundation for the breeding of fine dual-purpose sheep thereby accelerating the pace of sheep breeding.

The Dynamic Change of Gene-Regulated Networks in Cashmere Goat Skin with Seasonal Variation

The Cashmere goat (Capra hircus) is renowned for its high-quality fiber production trait. The hair cycle in Cashmere goat has an annual rhythm. To deepen the understanding of the molecular foundation of annual rhythm in the skin of Cashmere goat, we did a comparative analysis of the Cashmere goat skin transcriptome all year round. 4002 Differentially expressed genes (DEGs) were identified with seasonal variations. 12 months transcriptome were divided into four developmental stages: Jan-Mar, Apr-Jul, Aug-Oct, and Nov-Dec based on gene expression patterns. 13 modules of highly correlated genes in skin were identified using WGCNA. Ten of these modules were consistent with the development stages. The gene function of those genes in each module was analyzed by functional enrichment. The results indicated that Wnt and Hedgehog signaling pathways were inhibited from January to March and activated from April to July. The cutaneous immune system of Cashmere goats has high activity from August to October. Fatty acid metabolism dominates goat skin from November to December. This study provides new information related to the annual skin development cycle, which could provide molecular biological significance for understanding the seasonal development and response to the annual rhythm of skin.

MicroRNA-mRNA Regulatory Networking Fine-Tunes Polyunsaturated Fatty Acid Synthesis and Metabolism in the Inner Mongolia Cashmere Goat

Fatty acid composition is an important aspect of meat quality in ruminants. Improving the beneficial fatty acid level in cashmere goat meat is important to its economic value. To investigate microRNAs (miRNAs) and mRNAs that regulate or coregulate polyunsaturated fatty acid (PUFA) synthesis and metabolism in the Inner Mongolia cashmere goat, we used longissimus dorsi muscle (WLM) and biceps femoris muscle (WBM) for transcript-level sequencing. RT-qPCR was used to evaluate the expression of mRNAs and miRNAs associated with PUFA synthesis and metabolism. The total PUFA content in the WBM was significantly higher than that in the WLM (P < 0.05). Our study is the first to systematically report miRNAs in cashmere goat meat. At the mRNA level, 20,375 genes were identified. ACSL1, CD36 and TECRL were at the center of a gene regulatory network and contributed significantly to the accumulation and metabolic regulation of fatty acids. At the miRNA level, 426 known miRNAs and 30 novel miRNAs were identified. KEGG analysis revealed that the miRNA target genes were involved mainly in the PPAR signaling pathway. The mRNA-miRNA coregulation analysis showed that ACSL1 was negatively targeted by nine miRNAs: chi-miR-10a-5p, chi-miR-10b-5p, chi-miR-130b-5p, chi-miR-15a-5p_R-1, chi-miR-15b-5p, chi-miR-16a-5p, chi-miR-16b-5p, chi-miR-181c-5p_R+1, and chi-miR-26b-5p. Finally, we speculated that the simultaneous silencing of ACSL1 by one or more of these nine miRNAs through PPAR signaling led to low ACSL1 expression in the WLM and, ultimately to high PUFA content in the WBM. Our study helps elucidate the metabolic regulation of fatty acids in Inner Mongolia cashmere goats.