Key Laboratory of Qinghai-Tibetan Plateau Animal Genetic Resource Reservation and Utilization, Ministry of Education, Southwest Minzu University, Chengdu 610041, China.
Key Laboratory of Qinghai-Tibetan Plateau Animal Genetic Resource Reservation and Exploitation of Sichuan Province, Southwest Minzu University, Chengdu 610041, China.
Int J Mol Sci. 2023 Mar 2;24(5):4817. doi: 10.3390/ijms24054817.
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.
肌内脂肪有助于改善羊肉品质。N6-甲基腺苷(m6A)修饰的环状 RNA 在脂肪细胞分化和代谢中发挥重要作用。然而,在山羊肌内脂肪细胞分化前后,m6A 修饰环状 RNA 的机制仍知之甚少。在这里,我们进行了甲基化 RNA 免疫沉淀测序(MeRIP-seq)和环状 RNA 测序(circRNA-seq),以确定在山羊脂肪细胞分化过程中 m6A 修饰的环状 RNA 之间的差异。m6A-circRNA 图谱在肌内前脂肪细胞组中总共显示了 403 个环状 RNA 中的 427 个 m6A 峰,在成熟脂肪细胞组中显示了 401 个环状 RNA 中的 428 个峰。与肌内前脂肪细胞组相比,成熟脂肪细胞组中 75 个环状 RNA 中的 75 个峰差异显著。此外,肌内前脂肪细胞和成熟脂肪细胞的基因本体论(GO)和京都基因与基因组百科全书(KEGG)分析表明,差异 m6A 修饰的环状 RNA 富集于蛋白激酶 G(PKG)信号通路、内分泌和其他因子调节的钙重吸收、赖氨酸降解等途径中。m6A-circRNA-miRNA-mRNA 相互作用网络预测了不同山羊脂肪细胞中潜在的 m6A-circRNA 调控机制。我们的结果表明,在 12 个上调和 7 个下调的 m6A-circRNA 之间,通过 14 个和 11 个 miRNA 介导的通路分别存在复杂的调控关系。此外,联合分析显示 m6A 丰度与环状 RNA 表达水平之间存在正相关关系,如环状 RNA_0873 和环状 RNA_1161 的表达水平,这表明 m6A 可能在调节山羊脂肪细胞分化过程中的环状 RNA 表达中发挥重要作用。这些结果将为阐明 m6A-circRNA 在山羊肌内脂肪细胞分化中的生物学功能和调控特征提供新的信息,并有助于进一步进行分子育种,以改善山羊的肉质。