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

MiR-23b and miR-133 Cotarget TGFβ2/NOTCH1 in Sheep Dermal Fibroblasts, Affecting Hair Follicle Development

Wool is produced and controlled by hair follicles (HFs). However, little is known about the mechanisms involved in HF development and regulation. Sheep dermal fibroblasts (SDFs) play a key role in the initial stage of HF development. Analyzing the molecular mechanism that regulates early HF development in superfine wool sheep is of great importance for better understanding the HF morphogenesis process and for the breeding of fine wool sheep. Here, we show that two microRNAs (miRNAs) affect the development of HFs by targeting two genes that are expressed by SDFs. Meanwhile, the overexpression and inhibition of oar-miR-23b and oar-miR-133 in SDFs cells and cell proliferation, apoptosis, and migration were further detected using a CCK-8 assay, an Annexin V-FITC assay, a Transwell assay, and flow cytometry. We found that oar-miR-23b, oar-miR-133, and their cotarget genes TGFβ2 and NOTCH1 were differentially expressed during the six stages of HF development in superfine wool sheep. Oar-miR-23b and oar-miR-133 inhibited the proliferation and migration of SDFs and promoted the apoptosis of SDFs through TGFβ2 and NOTCH1. oar-miR-23b and oar-miR-133 inhibited the proliferation and migration of SDFs by jointly targeting TGFβ2 and NOTCH1, thereby inhibiting the development of superfine wool HFs. Our research provides a molecular marker that can be used to guide the breeding of ultrafine wool sheep.

Understanding Circular RNAs in Health, Welfare, and Productive Traits of Cattle, Goats, and Sheep

Circular RNAs (circRNAs) are unique noncoding RNA molecules, notable for their covalent closed-loop structures, which play a crucial role in regulating gene expression across a variety of biological processes. This review comprehensively synthesizes the existing knowledge of circRNAs in three key livestock species: Bos taurus (cattle), Ovis aries (sheep), and Capra hircus (goats). It focuses on their functional importance and emerging potential as biomarkers for disease detection, stress response, and overall physiological health. Specifically, it delves into the expression and functionality of circRNAs in these species, paying special attention to traits critical to livestock productivity such as milk production, meat quality, muscle development, wool production, immune responses, etc. We also address the current challenges faced in circRNA research, including the need for standardized methodologies and broader studies. By providing insights into the molecular mechanisms regulated by circRNAs, this review underscores their scientific and economic relevance in the livestock industry. The potential of circRNAs to improve animal health management and the quality of animal-derived products aligns with growing consumer concerns for animal welfare and sustainability. Thus, this paper aims to guide future research directions while supporting the development of innovative strategies in livestock management and breeding.

CircERCC6 Positively Regulates the Induced Activation of SHF Stem Cells in Cashmere Goats via the miR-412-3p/BNC2 Axis in an m6A-Dependent Manner

The cashmere, a kind of nature protein fiber, is one of the main use of cashmere goats. The induced activation of secondary hair follicle (SHF) stem cells by the dermal papilla cell-derived signals is a key biological process for the morphogenesis and growth of cashmere fiber in cashmere goats. Previously, the circRNA-ERCC6 (circERCC6) was identified from cashmere goat SHFs; however, its biological significance is unclear in the SHF physiology process of cashmere goats. In this study, we found that circERCC6 exhibited significantly higher expression at anagen SHF bulge compared with the counterpart of telogen and harbored three m6A modified sites (named m6A-685, m6A-862, and m6A-995) through methylation immunoprecipitation using a real-time quantitative polymerase chain reaction (Me-RIP-qPCR) technique. The knockdown experiments of circERCC6 in SHF stem cells showed that circERCC6 positively regulates the induced activation of SHF stem cells in cashmere goats. Through a dual-luciferase reporter assay, we demonstrated that m6A-modified circERCC6 (m6A-circERCC6) sponged miR-412-3p to upregulate the expression of BNC2 mRNA in SHFstem cells. Through m6A-deficient mutant assay in circERCC6 knockdown SHF stem cells, we further showed that m6A modification within circERCC6 is required to mediate the miR-412-3p/BNC2 axis to finally promote the proper induced activation of SHF stem cells in cashmere goats.

Genome-wide DNA methylation and transcriptome analyses reveal the key gene for wool type variation in sheep

BACKGROUND

Wool fibers are valuable materials for textile industry. Typical wool fibers are divided into medullated and non-medullated types, with the former generated from primary wool follicles and the latter by either primary or secondary wool follicles. The medullated wool is a common wool type in the ancestors of fine wool sheep before breeding. The fine wool sheep have a non-medullated coat. However, the critical period determining the type of wool follicles is the embryonic stage, which limits the phenotypic observation and variant contrast, making both selection and studies of wool type variation fairly difficult.

RESULTS

During the breeding of a modern fine (MF) wool sheep population with multiple-ovulation and embryo transfer technique, we serendipitously discovered lambs with ancestral-like coarse (ALC) wool. Whole-genome resequencing confirmed ALC wool lambs as a variant type from the MF wool population. We mapped the significantly associated methylation locus on chromosome 4 by using whole genome bisulfite sequencing signals, and in turn identified the SOSTDC1 gene as exons hypermethylated in ALC wool lambs compare to their half/full sibling MF wool lambs. Transcriptome sequencing found that SOSTDC1 was expressed dozens of times more in ALC wool lamb skin than that of MF and was at the top of all differentially expressed genes. An analogy with the transcriptome of coarse/fine wool breeds revealed that differentially expressed genes and enriched pathways at postnatal lamb stage in ALC/MF were highly similar to those at the embryonic stage in the former. Further experiments validated that the SOSTDC1 gene was specifically highly expressed in the nucleus of the dermal papilla of primary wool follicles.

CONCLUSION

In this study, we conducted genome-wide differential methylation site association analysis on differential wool type trait, and located the only CpG locus that strongly associated with primary wool follicle development. Combined with transcriptome analysis, SOSTDC1 was identified as the only gene at this locus that was specifically overexpressed in the primary wool follicle stem cells of ALC wool lamb skin. The discovery of this key gene and its epigenetic regulation contributes to understanding the domestication and breeding of fine wool sheep.

N6-Methyladenosine modification (m6A) of circRNA-ZNF638 contributes to the induced activation of SHF stem cells through miR-361-5p/Wnt5a axis in cashmere goats

OBJECTIVE

The objective of this study was to investigate the effects of N6-Methyladenosine modification-circRNA-zinc finger protein 638 (m6A-circRNA-ZNF638) on the induced activation of secondary hair follicle (SHF) stem cells with its potential mechanisms in cashmere goats.

METHODS

The m6A modification of ZNF638 was analyzed using methylation immunoprecipitation with real-time quantitative polymerase chain reaction technique in SHF stem cells. The effects of circRNA-ZNF638 on the induced activation of SHF stem cells in m6A dependence were evaluated through the overexpression of circRNA-ZNF638/its m6Adeficient mutants in circRNA-ZNF638 knockdown SHF stem cells. The competitive binding of miR-361-5p to circRNA-ZNF638/Wnt5a 3'- untranslated region was analyzed through Dual-luciferase reporter assay.

RESULTS

The m6A-circRNA-ZNF638 had significantly higher transcription at anagen SHF bulge of cashmere goats compared with that at telogen, as well as it positively regulated the induced activation of SHF-stem cells in cashmere goats. Mechanismly, m6A-circRNA-ZNF638 sponged miR-361-5p to heighten the transcriptional expression of Wnt5a gene in SHFstem cells. We further demonstrated that the internal m6A modification within circRNAZNF638 is required for mediating the miR-361-5p/Wnt5a pathway to regulate the induced activation of SHF stem cells through an introducing of m6A-deficient mutant of circRNAZNF638.

CONCLUSION

The circRNA-ZNF638 contributes the proper induced activation of SHF-stem cells in cashmere goats in m6A-dependent manner through miR-361-5p/Wnt5a axis.

The Important Role of m6A-Modified circRNAs in the Differentiation of Intramuscular Adipocytes in Goats Based on MeRIP Sequencing Analysis

Intramuscular fat contributes to the improvement of goat meat quality. N⁶-Methyladenosine (m6A)-modified circular RNAs play important roles in adipocyte differentiation and metabolism. However, the mechanisms by which m6A modifies circRNA before and after differentiation of goat intramuscular adipocytes remain poorly understood. Here, we performed methylated RNA immunoprecipitation sequencing (MeRIP-seq) and circRNA sequencing (circRNA-seq) to determine the distinctions in m6A-methylated circRNAs during goat adipocyte differentiation. The profile of m6A-circRNA showed a total of 427 m6A peaks within 403 circRNAs in the intramuscular preadipocytes group, and 428 peaks within 401 circRNAs in the mature adipocytes group. Compared with the intramuscular preadipocytes group, 75 peaks within 75 circRNAs were significantly different in the mature adipocytes group. Furthermore, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of intramuscular preadipocytes and mature adipocytes showed that the differentially m6A-modified circRNAs were enriched in the PKG signaling pathway, endocrine and other factor-regulated calcium reabsorption, lysine degradation, etc. m6A-circRNA-miRNA-mRNA interaction networks predicted the potential m6A-circRNA regulation mechanism in different goat adipocytes. Our results indicate that there is a complicated regulatory relationship between the 12 upregulated and 7 downregulated m6A-circRNAs through 14 and 11 miRNA mediated pathways, respectively. In addition, co-analysis revealed a positive association between m6A abundance and levels of circRNA expression, such as expression levels of circRNA_0873 and circRNA_1161, which showed that m6A may play a vital role in modulating circRNA expression during goat adipocyte differentiation. These results would provide novel information for elucidating the biological functions and regulatory characteristics of m6A-circRNAs in intramuscular adipocyte differentiation and could be helpful for further molecular breeding to improve meat quality in goats.

Identification and functional analysis of m6A in the mammary gland tissues of dairy goats at the early and peak lactation stages

N6-methyladenosine (m⁶A) is the most common reversible epigenetic RNA modification in the mRNA of all higher eukaryotic organisms and plays an important role in the regulation of gene expression and cell function. In this study, m⁶A-modified methylated RNA immunoprecipitation sequencing (MeRIP-seq) and transcriptome sequencing (RNA-seq) were used to identify the key genes with m⁶A modification during mammary gland development and lactation in dairy goats. The results showed that m⁶A methylation occurred at 3,927 loci, which were significantly enriched in the 3′ untranslated region (3′UTR) and the termination codon region. In the early stage and peak stage of lactation, m⁶A methylation occurred extensively in mammary tissues, and a total of 725 differentially expressed m⁶A-modified genes were obtained, all negatively correlated with mRNA expression. In addition, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that different methylated genes were mainly involved in the growth and apoptosis of mammary epithelial cells through signaling pathways, such as the mitogen-activated protein kinase (MAPK) and phospholipase D pathways, and then affected the development and lactation of mammary gland. All in all, we identified and analyzed the methylation events related to the development and lactation regulation of mammary gland at the early and peak lactation stages, and provided a theoretical basis to reveal the physiological regulatory system of mammary gland development and lactation in dairy goats.