Molecular characterization, expression and methylation status analysis of BMP4 gene in skin tissue of Liaoning cashmere goat during hair follicle cycle

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.

CeRNA regulates network and expression and SNP effect on NFKBIA of cashmere fineness

In this study, a total of 1140 Liaoning Cashmere Goats (LCG) were genotyped for single nucleotide polymorphism (SNP) of NFKBIA gene. There are 15 SNPs and 7 genotypes have been found, and G1547A (GG) genotype has been associated with cashmere fineness and cashmere yield. An integrated ceRNA regulatory network of NFKBIA gene was made. To prove NFKBIA and these non-coding RNAs (ncRNAs) may be related to cashmere fineness, we performed qPCR on these ncRNA in LCG coarse type skin (CT-LCG) and LCG fine type skin (FT-LCG). The result of qPCR showed lncRNA XLOC_011060 and ciRNA452 are at high expression level in CT-LCG, all miRNAs appear high expressed in FT-LCG, and mir-93 was the most significant difference between CT-LCG and FT-LCG. In addition, five miRNAs were selected for qPCR in different genotypes. The qPCR results showed that mir-93 might negatively regulate cashmere fineness and mir-17-5p may play a positive role in regulating cashmere fineness of individuals with G1355A (AG) genotype. These results demonstrated that NFKBIA gene is associated with cashmere fineness of LCG and G1547A (GG) genotype is the preferred marker genotype for cashmere fineness.

Association between the cashmere production performance, milk production performance, and body size traits and polymorphism of COL6A5 and LOC102181374 genes in Liaoning cashmere goats

The purpose of this study was to analyze the relationship between COL6A5 (collagen type VI alpha 5 chain) and LOC102181374 (alcohol dehydrogenase 1) genes and the production performance of Liaoning cashmere goats by single nucleotide polymorphism (SNP). We have searched for SNP loci of COL6A5 and LOC102181374 genes through sequence alignment and PCR experiments, and have used SPSS and SHEsis software to analyze production data. We obtained five SNP loci in total, including three SNP loci (G50985A, G51140T, G51175A) in COL6A5 gene and two SNP loci (A10067G, T10108C) in LOC102181374 gene. The genotypes G50985A (AG), G51140T (GT), G51175A (AA), A10067G (AA), and T10108C (CC) of these loci have certain advantages in improving the production performance of Liaoning cashmere goats. The haplotype combinations that can improve production performance in COL6A5 gene were H1H5:AGGGAG, H4H4:GGGGAA, and H4H4:GGGGAA. H3H3:GGCC and H2H4:AGTT were the dominant combinations in LOC102181374 gene. At G51175A and A10067G loci, we found that H1H2:AAAG and H1H3:AGAA have dominant effects. These results may provide some support for the molecular breeding of production traits in Liaoning cashmere goats.

Multi-omics analysis of functional substances and expression verification in cashmere fineness

BACKGROUND

Numerous factors influence the growth and development of cashmere. Existing research on cashmere has predominantly emphasized a single omics level. Integrating multi-omics analyses can offer a more comprehensive understanding by encompassing the entire spectrum. This study more accurately and comprehensively identified the key factors influencing cashmere fineness using multi-omics analysis.

METHODS

This study used skin tissues of coarse cashmere type (CT_LCG) and fine cashmere type Liaoning cashmere goats (FT_LCG) for the analysis. This study employed an integrated approach involving transcriptomics, translatomics, proteomics, and metabolomics to identify substances associated with cashmere fineness. The findings were validated using parallel reaction monitoring (PRM) and multiple reaction monitoring (MRM) techniques.

RESULTS

The GO functional enrichment analysis identified three common terms: multicellular organismal process, immune system process, and extracellular region. Furthermore, the KEGG enrichment analysis uncovered the involvement of the arachidonic acid metabolic pathway. Protein expression trends were verified using PRM technology. The expression trends of KRT79, as confirmed by PRM, were consistent with those observed in TMT proteomics and exhibited a positive regulatory effect on cashmere fineness. Metabolite expression trends were confirmed using MRM technology. The expression trends of 9 out of 15 validated metabolites were in agreement with those identified in the non-targeted metabolomics analysis.

CONCLUSIONS

This study employed multi-omics analysis to identify key regulators of cashmere fineness, including PLA2G12A, KRT79, and prostaglandin B2. The findings of this study offer valuable data and establish a theoretical foundation for conducting comprehensive investigations into the molecular regulatory mechanisms and functional aspects of cashmere fineness.

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 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.

Selection of Cashmere Fineness Functional Genes by Translatomics

Cashmere fineness is an important index to evaluate cashmere quality. Liaoning Cashmere Goat (LCG) has a large cashmere production and long cashmere fiber, but its fineness is not ideal. Therefore, it is important to find genes involved in cashmere fineness that can be used in future endeavors aiming to improve this phenotype. With the continuous advancement of research, the regulation of cashmere fineness has made new developments through high-throughput sequencing and genome-wide association analysis. It has been found that translatomics can identify genes associated with phenotypic traits. Through translatomic analysis, the skin tissue of LCG sample groups differing in cashmere fineness was sequenced by Ribo-seq. With these data, we identified 529 differentially expressed genes between the sample groups among the 27197 expressed genes. From these, 343 genes were upregulated in the fine LCG group in relation to the coarse LCG group, and 186 were downregulated in the same relationship. Through GO enrichment analysis and KEGG enrichment analysis of differential genes, the biological functions and pathways of differential genes can be found. In the GO enrichment analysis, 491 genes were significantly enriched, and the functional region was mainly in the extracellular region. In the KEGG enrichment analysis, the enrichment of the human papillomavirus infection pathway was seen the most. We found that the COL6A5 gene may affect cashmere fineness.

BMP7 Functions to Regulate Proliferation of Dermal Papilla Cells in Hu Sheep

Bone morphogenetic proteins (BMPs) are the structurally similar and highly conserved type of functional proteins that play an important role in hair follicle growth and development. BMP7 was a differentially expressed gene in different patterns of Hu sheep lambskin identified using Agilent microarray. Since hair follicle is the basis of pattern formation of lambskin, and its growth and development is governed by dermal papilla cells (DPCs), to clarify the role of BMP7 and hair follicle, our study was designed to investigate the regulation between BMP7 and DPCs. Firstly, the CDS region of BMP7 was cloned by 3'Race and PCR in Hu sheep and performed serious of bioinformatic analysis. Then, the effects of BMP7 on DPCs were analyzed after overexpression and interference of BMP7 in dermal papilla cells by CCK8, EdU, and PI assay. Additionally, qPCR was also conducted to clarify the relationship between BMP7 and the TGF-β/Smad signaling pathway. A total of 1,296 bp of the BMP7 CDS region sequence was sucessfully cloned in Hu sheep, encoding a signal peptide of 431 amino acids, molecular weight was 49,316.9Da and the isoelectric point (Pi) was 7.75. Nucleotide sequencing analysis of BMP7 revealed that Hu sheep had high homology with Bos taurus, Homo sapiens, and Canis lupus familiaris. Structure domain prediction showed that TGF-β superfamily domain exist between 330th-431th amino acid, BMP7 protein is a secreted protein. In BMP7 up-regulated DPCs, DPCs proliferation rate and cell cycle were significantly higher than that of NC group (P < 0.05). Meanwhile, the expression level of Smad3, Smad4, Samd6, and TGF-β1 in TGF-β/Smad signaling pathway were significantly lower than that in NC group (P < 0.05). In BMP7 down-regulated DPCs, it presented the opposite result. In conclusion, our study showed that BMP7 had a positive effect on DPCs by accelerating the proliferation and cell cycle of DPCs, and hypothesized that regulate hair follicles growth and development via TGF-β/Smad signaling pathway. These findings may provide a synergistic target for the subsequent research of hair follicle growth and development.

m6A Methylation Analysis of Differentially Expressed Genes in Skin Tissues of Coarse and Fine Type Liaoning Cashmere Goats

N6-methyladenosine (m6A) is the most common internal modification in mRNAs of all higher eukaryotes. Here we perform two high-throughput sequencing methods, m6A-modified RNA immunoprecipitation sequence (MeRIP-seq) and RNA sequence (RNA-seq) to identify key genes with m6A modification in cashmere fiber growth. A total of 9,085 m6A sites were differentially RNA m6A methylated as reported from by MeRIP-seq, including 7,170 upregulated and 1,915 downregulated. In addition, by comparing m6A-modified genes between the fine-type Liaoning cashmere goat (FT-LCG) and coarse-type Liaoning Cashmere Goat (CT-LCG) skin samples, we obtain 1,170 differentially expressed genes. In order to identify the differently methylated genes related to cashmere fiber growth, 19 genes were selected to validate by performing qRT-PCR in FT-LCG and CT-LCG. In addition, GO enrichment analysis shows that differently methylated genes are mainly involved in keratin filament and intermediate filament. These findings provide a theoretical basis for future research on the function of m6A modification during the growth of cashmere fiber.

The Overexpression of Tβ4 in the Hair Follicle Tissue of Alpas Cashmere Goats Increases Cashmere Yield and Promotes Hair Follicle Development

Increased cashmere yield and improved quality are some goals of cashmere goat breeding. Thymosin beta-4 (Tβ4) plays a key role in the growth and development of hair follicles. For the past ten years, we have evaluated the role of Tβ4 by establishing a flock of 15 cashmere goats that specifically overexpress the Tβ4 gene in the hair follicles. These Tβ4 overexpression (Tβ4-OE) cashmere goats had more secondary hair follicles than the WT goats and produced more cashmere. Meanwhile, combined analysis of the skin transcriptome and proteome in cashmere goats suggested that Tβ4 may affect hair growth by interacting with keratin type II cytoskeletal 4 epidermal (KRT4) to mediate the extracellular signal-regulated protein kinase (ERK) signaling pathway, thereby promoting the development of secondary hair follicles, and consequently, increasing cashmere yield. Thus, the specific overexpression of Tβ4 in the hair follicles of cashmere goats effectively increased the cashmere yield.

Preliminary study on microR-148a and microR-10a in dermal papilla cells of Hu sheep

Background

Hu sheep, a unique Chinese breed with high reproductive performance, are also well known for their rare white lambskin in China. The quality of lambskin is affected by hair follicles, and dermal papilla cells are an important component of hair follicles that plays a key role in hair follicle growth and development. This study helps elucidate the effect of miR-148a and miR-10a on hair follicle growth and development.

Results

Based on the results of gene chip and high-throughput sequencing, bone morphogenetic protein 7 (BMP7) was used as a research object. Bioinformatics analysis and the dual-luciferase reporter system indicated that, along with Western blot and quantitative real-time polymerase chain reaction (qRT-PCR) that miR-148a and miR-10a target relationships with BMP7. BMP7 was the target gene both for miR-148a and miR-10a by the dual-luciferase reporter system and Western blot. Hu sheep dermal papilla cells were successfully isolated and purified, and after transfecting miR-148a/miR-10a mimics and inhibitors into dermal papilla cells, a Cell Counting Kit-8 (CCK-8) was used to determine that miR-148a/miR-10a inhibited the proliferation of Hu sheep dermal papilla cells. In addition, after the overexpression of miR-148a, the expression levels of Smad3 (P < 0.05), Smad6 (P < 0.05), Smad4 (P < 0.01), and Smad5 (P < 0.01) were significantly higher than those of the control groups. After the inhibition of miR-148a, the expression levels of Smad3 (P < 0.05), Smad4 (P < 0.05), and TGF-β (P < 0.01) were significantly lower than those of the control groups. After the overexpression of miR-10a, the expression levels of Smad1 (P < 0.01), Smad2 (P < 0.05), Smad4 (P < 0.01), Smad5 (P < 0.01), and TGF-β (P < 0.05) were significantly lower than those of the control groups. After the inhibition of miR-10a, the expression levels of Smad1 (P < 0.01) and Smad2 (P < 0.05) were significantly lower than those of the control groups.

Conclusions

These results revealed the target relationship between miR-148a, miR-10a and BMP7, and the effect of miR-148a and miR-10a on the proliferation of dermal papilla cells. They will provide the basis for a follow-up study on how miR-148a, and miR-10a mediate BMP7 regulation of hair follicle growth and development.

An Integrated Analysis of Cashmere Fineness lncRNAs in Cashmere Goats

Animal growth and development are regulated by long non-coding RNAs (lncRNAs). However, the functions of lncRNAs in regulating cashmere fineness are poorly understood. To identify the key lncRNAs that are related to cashmere fineness in skin, we have collected skin samples of Liaoning cashmere goats (LCG) and Inner Mongolia cashmere goats (MCG) in the anagen phase, and have performed RNA sequencing (RNA-seq) approach on these samples. The high-throughput sequencing and bioinformatics analyses identified 437 novel lncRNAs, including 93 differentially expressed lncRNAs. We also identified 3084 differentially expressed messenger RNAs (mRNAs) out of 27,947 mRNAs. Gene ontology (GO) analyses of lncRNAs and target genes in cis show a predominant enrichment of targets that are related to intermediate filament and intermediate filament cytoskeleton. According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, sphingolipid metabolism is a significant pathway for lncRNA targets. In addition, this is the first report to reveal the possible lncRNA–mRNA regulatory network for cashmere fineness in cashmere goats. We also found that lncRNA XLOC_008679 and its target gene, KRT35, may be related to cashmere fineness in the anagen phase. The characterization and expression analyses of lncRNAs will facilitate future studies on the potential value of fiber development in LCG.

Transcriptomic analysis reveals critical genes for the hair follicle of Inner Mongolia cashmere goat from catagen to telogen

There are two main types of hair follicle in Inner Mongolia Cashmere goats, the primary hair follicle (PHF) producing hair fibers and the secondary hair follicle (SHF) producing cashmere fibers. Of both fibers from cashmere-bearing goats in Aerbasi, Inner Mongolia, the timing of cyclical phases for the cashmere have been well clarified but hair fibers have been less noticeable. Herein, we evaluated transcriptome of PHF and SHF from the same three goats in Aerbasi at the catagen- and telogen phase of cashmere growth. We totally found 1977 DEGs between PHFs at the telogen and catagen phases of SHF, 1199 DEGs between telogen- and catagen SHF, 2629 DEGs between PHF at the catagen phase of SHF and catagen SHF, and 755 DEGs between PHF at the telogen phase of SHF and telogen SHF. By analyzing gene functions based on GO and KEGG database, we found that the DEGs have functions in muscle contraction and muscle filament sliding between catagen- and telogen SHF, indicating that arrector pilli muscles might play a role on the transition from catagen to telogen. Moreover, considering that the enriched GO and KEGG categories of the DEGs between PHF and SHF, we suggested that part of PHF might rest in their own anagen phase when SHF are at catagen, but PHF might enter into the telogen phase at SHF’s telogen. Finally, we highly recommended the several potential genes acting as the regulators of the transition between growth phases including IL17RB and eight members of ZNF. These results provide insight into molecular mechanisms on the transition of SHF from catagen to telogen together with PHF’s growth situation at SHF’s catagen and telogen in Inner Mongolia Cashmere goats.

A lncRNA-H19 transcript from secondary hair follicle of Liaoning cashmere goat: Identification, regulatory network and expression regulated potentially by its promoter methylation

The H19 transcript (imprinted maternally expressed transcript) is well-known as long noncoding RNA (lncRNA), and it is thought to be associated with the inductive capacity of dermal papilla cells for hair-follicle reconstruction. In this study, we isolated and characterized a lncRNA-H19 transcript from the secondary hair follicle of Liaoning cashmere goat. Also, we investigated its transcriptional pattern and methylation status of H19 gene in secondary hair follicle of this breed during different stages of hair follicle cycle. Nucleotide composition analysis indicated that guanine (G) and cytosine (C) are the dominant nucleotides in the lncRNA-H19 transcript of Liaoning cashmere goat with the highest frequency distribution (11.25%) of GG nucleotide pair. The regulatory network showed that lncRNA-H19 transcript appears to have remarkably diverse regulatory relationships with its related miRNAs and the potential target genes. In secondary hair follicle, the relative expression of lncRNA-H19 transcript at the anagen phase is significantly higher than that at both telogen and catagen phases suggesting that lncRNA-H19 transcript might play essential roles in the formation and growth of cashmere fiber of goat. Methylation analysis indicated that the methylation of the promoter region of H19 gene most likely participates in its transcriptional suppression in secondary hair follicle of Liaoning cashmere goat.