Transcriptome-wide N 6-Methyladenosine Methylome Profiling Reveals m6A Regulation of Skeletal Myoblast Differentiation in Cattle (Bos taurus)

Transcriptome-wide study of mRNAs modified by m6A RNA methylation in the testis development of dairy goats

N6-methyladenosine (m6A) is an important epigenetic modification in RNA, playing a crucial role in regulating the production and aging of animal testicular sperm. This study extracted mRNA from the testicular tissue of male goats before and after sexual maturity, generating a methylation map through preliminary experiments and methylation immunoprecipitation sequencing. The results showed that during the development of dairy goats, the expression levels of marker genes related to testicular development and methylation-related enzymes changed significantly. A total of 36,602 peaks and 11,223 genes were identified in the two groups, including 2989 differential peaks (427 upregulated and 2562 downregulated) and 1457 differentially expressed genes (833 upregulated and 624 downregulated). The abundance of m6A was positively correlated with gene expression levels. This study reports for the first time the mRNA profiles of m6A modifications across the entire transcriptome during testicular development in Guanzhong dairy goats, providing a new perspective for genetic improvement in goats.

Estimation of Genetic Parameters and GWAS on Water Efficiency Traits in the Senepol Cattle

The need for producing in environmentally resilient system drives new research to achieve sustainable beef production. Water footprint of the beef supply chain is a concern that must be addressed, aiming to improve water use within the production chain. One approach is genetic selection of beef cattle for water efficiency. However, it is essential to understand the genetic architecture and mechanisms involved in the expression of this phenotype to choose the best selection criteria. Thus, our study aimed to estimate genetic parameters for water efficiency traits, conduct a genome-wide association study (GWAS) and identify the genetic networks and biological processes involved. A population of 1762 purebred Senepol cattle was phenotyped for the following water efficiency traits: water intake (WI), gross water efficiency (GWE), water conversion ratio (WCR), residual water intake based on average daily gain (RWIADG) and residual water intake based on dry matter intake (RWIDMI). A subset of 1342 animals was genotyped using GGP Bovine 50 K SNP Chip with (734 animals) or 100 K (508 animals), and imputation from 50 K to 100 K was performed with Beagle software. The heritability estimates were 0.36 ± 0.06, 0.26 ± 0.05, 0.22 ± 0.05, 0.24 ± 0.05 and 0.20 ± 0.05 for WI, GWE, WCR, RWIADG and RWIDMI, respectively. Unlike the raw measures of WI, the phenotypic correlations between average daily gain (ADG) and the residuals (RWIDMI and RWIADG) were zero. All water efficiency traits were moderately to highly correlated with each other. GWAS were used to estimate the effect of 79,860 single nucleotide polymorphisms (SNPs), and significant SNPs were only observed for WCR. Enrichment analysis of genes in the significant regions revealed the involvement of different biological processes, such as saliva production, water transport, renal system and immune system. Genetic selection of Senepol cattle for water efficiency traits is feasible and can reduce water requirements for meat production. Water efficiency measures are polygenic traits, and different biological processes act simultaneously on the expression of related phenotypes.

The Spatial-Temporal Alternative Splicing Profile Reveals the Functional Diversity of FXR1 Isoforms in Myogenesis

Alternative splicing (AS) is a fundamental mechanism contributing to proteome diversity, yet its comprehensive landscape and regulatory dynamics during skeletal muscle development remain largely unexplored. Here, the temporal AS profiles are investigated during myogenesis in five vertebrates, conducting comprehensive profiling across 27 developmental stages in skeletal muscle and encompassing ten tissues in adult pigs. The analysis reveals a pervasive and evolutionarily conserved pattern of alternative exon usage throughout myogenic differentiation, with hundreds of skipped exons (SEs) showing developmental regulation, particularly within skeletal muscle. Notably, this study identifies a muscle-specific SE (exon 15) within the Fxr1 gene, whose AS generates two dynamically expressed isoforms with distinct functions: the isoform without exon 15 (Fxr1E15 -) regulates myoblasts proliferation, while the isoform incorporating exon 15 (Fxr1E15+) promotes myogenic differentiation and fusion. Transcriptome analysis suggests that specifically knocking-down Fxr1E15+ isoform in myoblasts modulates differentiation by influencing gene expression and splicing of specific targets. The increased inclusion of exon 15 during differentiation is mediated by the binding of Rbm24 to the intron. Furthermore, in vivo experiments indicate that the Fxr1E15+ isoform facilitates muscle regeneration. Collectively, these findings provide a comprehensive resource for AS studies in skeletal muscle development, underscoring the diverse functions and regulatory mechanisms governing distinct Fxr1 isoforms in myogenesis.

Comprehensive analysis of the expression patterns and function of the FTO-LINE1 axis in yak tissues and muscle satellite cells

Background

The long interspersed nuclear element 1 (LINE1) retrotransposon has been identified as a specific substrate for fat mass and obesity-related gene (FTO), which facilitates the removal of N6-methyladenosine modifications from its targeted RNAs.

Methods

This study examined the dynamic interaction between FTO and LINE1 in yak tissues and muscle satellite cells, utilizing RT-qPCR, RNA immunoprecipitation (RIP), immunofluorescence staining, and techniques involving overexpression and interference of FTO and LINE1 to elucidate the relationship between FTO and LINE1 in yak tissues and muscle satellite cells.

Results

Cloning and analysis of the FTO coding sequence in Jiulong yak revealed a conserved protein structure across various Bos breeds, with notable homology observed with domestic yak, domestic cattle, and Java bison. Comprehensive examination of FTO and LINE1 gene expression patterns in Jiulong yaks revealed consistent trends across tissues in both sexes. FTO mRNA levels were markedly elevated in the heart and kidney, while LINE1 RNA was predominantly expressed in the heart. Immunoprecipitation confirmed the direct interaction between the FTO protein and LINE1 RNA in yak tissues and muscle satellite cells. The FTO-LINE1 axis was confirmed by a significant decrease in LINE1 RNA enrichment following its expression interference in yak muscle satellite cells. Overexpression of FTO substantially reduced the expression of recombinant myogenic factor 5 (MYF5). However, FTO interference had no discernible effect on MYF5 and myoblast determination protein 1 (MYOD1) mRNA levels. Immunofluorescence analysis revealed no alterations in Ki-67 protein expression following FTO interference or overexpression. However, phalloidin staining demonstrated enhancement in the myotube fusion rate of yak muscle satellite cells upon LINE1 interference.

Conclusion

This comprehensive mapping of the FTO and LINE1 mRNA expression patterns establishes a direct interaction between the FTO protein and LINE1 RNA in yak. The findings suggest that FTO overexpression promotes muscle satellite cells differentiation, whereas LINE1 negatively regulates myotube fusion. The study provides fundamental insights into the role of the FTO-LINE1 axis in determining the fate of muscle satellite cells in yak, laying a solid theoretical foundation for future investigations.

Proteomic Analysis Reveals the Effects of Different Dietary Protein Levels on Growth and Development of Jersey-Yak

Jersey-yak is a hybrid offspring of Jersey cattle and yak (Bos grunniens). Changing the feeding system of Jersey-yak can significantly improve its growth performance. In this study, tandem mass tag (TMT) proteomics technology was used to determine the differentially expressed proteins (DEPs) of the longissimus lumborum (LL) muscle of Jersey-yak fed different protein levels of diet. The results showed that compared with the traditional grazing feeding, the growth performance of Jersey-yaks was significantly improved by crude protein supplementation after grazing. A total of 3368 proteins were detected in these muscle samples, of which 3365 were quantified. A total of 434 DEPs were identified. Through analyses, it was found that some pathways related to muscle growth and development were significantly enriched, such as Rap1 signaling pathway, mTOR signaling pathway, and TGF-beta signaling pathway. A number of DEPs enriched in these pathways are related to muscle cell development, differentiation, and muscle development, including integrin subunit alpha 7 (ITGA7), myosin heavy chain 8 (MYH8), and collagen type XII alpha 1 chain (COL12A1). In conclusion, the results of this study provide insights into the proteomics of different feeding patterns of Jersey-yak, providing a stronger basis for further understanding the biological mechanism of hybrid varieties.

Regulatory role of N6-methyladenosine in intramuscular fat deposition in chicken

Intramuscular fat (IMF) has a pivotal influence on meat quality, with its deposition being a multifaceted physiological interaction of several regulatory factors. N6-methyladenosine (m6A), the preeminent epigenetic alteration among eukaryotic RNA modifications, holds a crucial role in moderating post-transcriptional gene expression. However, there is a dearth of comprehensive understanding regarding the functional machinery of m6A modification in the context of IMF deposition in poultry. Our current study entails an analysis of the disparities in IMF within the breast and leg of 180-day-old Jingyuan chickens. We implemented methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) to delve into the distribution of m6A and its putative regulatory frameworks on IMF deposition in chickens. The findings demonstrated a markedly higher IMF content in leg relative to breast (P < 0.01). Furthermore, the expression of METTL14, WTAP, FTO, and ALKBH5 was significantly diminished in comparison to that of breast (P < 0.01). The m6A peaks in the breast and leg primarily populated 3'untranslated regions (3'UTR) and coding sequence (CDS) regions. The leg, when juxtaposed with the breast, manifested 176 differentially methylated genes (DMGs), including 151 hyper-methylated DMGs and 25 hypo-methylated DMGs. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed a pronounced enrichment of DMGs in the biosynthesis of amino acids, peroxisome, Fatty acid biosynthesis, fatty acid elongation, and cell adhesion molecules (CAMs) pathways. Key DMGs, namely ECH1, BCAT1, and CYP1B1 were implicated in the regulation of muscle lipid anabolism. Our study offers substantial insight and forms a robust foundation for further exploration of the functional mechanisms of m6A modification in modulating IMF deposition. This holds profound theoretical importance for improving and leveraging meat quality in indigenous chicken breeds.

Interactive regulation of DNA demethylase gene TET1 and m6A methyltransferase gene METTL3 in myoblast differentiation

DNA methylation (5mC) and mRNA N⁶-methyladenosine (m⁶A) play an essential role in gene transcriptional regulation. DNA methylation has been well established to be involved in skeletal muscle development. Interacting regulatory mechanisms between DNA methylation and mRNA m⁶A modification have been identified in a variety of biological processes. However, the effect of m⁶A on skeletal muscle differentiation and the underlying mechanisms are still unclear. It is also unknown whether there is an interaction between DNA methylation and mRNA m⁶A modification in skeletal myogenesis. In the present study, we used m⁶A-IP-qPCR, LC-MS/MS and dot blot assays to determine that the DNA demethylase gene, TET1, exhibited increased m⁶A levels and decreased mRNA expression during bovine skeletal myoblast differentiation. Dual-luciferase reporter assays and RIP experiments demonstrated that METTL3 suppressed TET1 expression by regulating TET1 mRNA stability in a m⁶A-YTHDF2-dependent manner. Furthermore, TET1 mediated DNA demethylation of itself, MYOD1 and MYOG, thereby stimulating their expression to promote myogenic differentiation. Ectopic expression of TET1 rescued the effect of METTL3 knockdown on reduced myotubes. In contrast, TET1 knockdown impaired the myogenic differentiation promoted by METTL3 overexpression. Moreover, ChIP experiments found that TET1 could bind and demethylate METTL3 DNA, which enhanced METTL3 expression. In addition, TET1 knockdown decreased m⁶A levels. ChIP assays also showed that TET1 knockdown contributed to the binding of H3K4me3 and H3K27me3 to METTL3 DNA. Our results revealed a negative feedback regulatory loop between TET1 and METTL3 in myoblast differentiation, which unveiled the interplay among DNA methylation, RNA methylation and histone methylation in skeletal myogenesis.

Transcriptome-Wide Study of mRNAs and lncRNAs Modified by m6A RNA Methylation in the Longissimus Dorsi Muscle Development of Cattle-Yak

Cattle-yak is a hybrid F₁ generation of cattle and yak, which has a history of more than 3000 years and has shown better production performance and higher economic benefits than those of yaks. However, up to now, there has been no study on the transcriptome-wide m⁶A methylation profile of bovine skeletal muscle and its potential biological function during muscle development. Here, we observed significant changes in the expression levels of muscle-related marker genes and methylation-related enzymes during the development of cattle-yak, and the overall m⁶A content in the Longissimus dorsi muscle of 18-month-old cattle-yak decreased significantly. A total of 36,602 peaks, 11,223 genes and 8388 lncRNAs were identified in the two groups, including 2989 differential peaks (427 up-regulated peaks and 2562 down-regulated peaks), 1457 differentially expressed genes (833 up-regulated genes and 624 down-regulated genes) and 857 differentially expressed lncRNAs (293 up-regulated lncRNAs and 564 down-regulated lncRNAs). GO and KEGG analysis revealed that they were significantly enriched in some muscle-related pathways (Wnt signaling pathway and MAPK signaling pathway) and high-altitude adaptation-related pathway (HIF-1 signaling pathway). Moreover, m⁶A abundance was positively correlated with gene expression levels, while it was negatively correlated with lncRNA expression levels. This indicates that m⁶A modification played an important role in the Longissimus dorsi muscle development of cattle-yak; however, the regulation mechanism of m⁶A-modified mRNA and lncRNA may be different. This study was the first report of transcriptome-wide m⁶A-modified mRNAs and lncRNAs atlas in the Longissimus dorsi muscle development of cattle-yak, one which will provide new perspectives for genetic improvement in bovines.

Roles of MEF2A and MyoG in the transcriptional regulation of bovine LATS2 gene

As an important downstream effector gene in the hippo signaling pathway, large tumor suppressor gene 2 (LATS2) is involved in cell proliferation and differentiation, organ size and tissue regeneration, and plays an important role in regulating the growth and development of animal muscles. The purpose of this study is to explore the temporal expression of bovine LATS2 gene, and determine the key transcription factors for regulating bovine LATS2 gene. The result showed that bovine LATS2 gene was highly expressed in liver and longissimus dorsi, and was up-regulated in infancy muscle. In addition, it was highly expressed on the 2th day during the differentiation stage of myoblast. The upstream 1.7 Kb sequence of the 5 'translation region of bovine LATS2 gene was cloned, and 7 different deletion fragments were amplified by the upstream primers. These fragments were constructed into double luciferase reporter vectors and transfected into myoblasts and myotubes cells, respectively to detect the core promoter regions. In addition, the key transcription factors of the core promoter sequence of the bovine LATS2 gene were analyzed and predicted by online software. Combining with site-directed mutations, siRNA interference and chromatin immunoprecipitation technology, it was identified that MEF2A and MyoG combined in core promoter region (-248/-56) to regulate the transcription activity of bovine LATS2 gene. The results have laid a theoretical foundation for exploring the molecular regulation mechanism of LATS2 gene in the process of muscle growth.

Research progress of m6A regulation during animal growth and development

Environmental factors, genetic factors, and epigenetics are involved in animal growth and development. Among them, methylation is one of the abundant modifications of epigenetics. N6-methyladenosine(m⁶A) is extensive in cellular RNA, of which mRNA is the most common internal modification. m⁶A modification regulates life activities dynamically and reversibly, including expressed genes, RNA metabolism, and protein translation. The m⁶A modifications are closely related to human diseases involving heart failure, tumors, and cancer. It is relatively in-depth in the medical field. However, there are few studies on its biochemical function in animals. We summarized the latest paper related to the chemical structure and role of the writers, the erasers, and the readers to study exerting dynamic regulation of m⁶A modification of animal growth and development. Furthermore, the key roles of m⁶A modification were reported in the process of RNA metabolism. Finally, the dynamic regulation of m⁶A modification in animal growth and development was reviewed, including brain development, fertility, fat deposition, and muscle production. It reveals the key roles of m⁶A modification and the regulation of gene expression, aiming to provide new ideas for m⁶A methylation in animal growth and development.

Transcriptome-wide study revealed m6A and miRNA regulation of embryonic breast muscle development in Wenchang chickens

N⁶-Methyladenosine (m6A) modification has been shown to play important role in skeletal muscle development. Wenchang chickens are commonly used as a high-quality animal model in researching meat quality. However, there have been no previous reports regarding the profile of m6A and its function in the embryonic breast muscle development of Wenchang chickens. In this paper, we identified different developmental stages of breast muscle in Wenchang chickens and performed m6A sequencing and miRNA sequencing in the breast muscle of embryos. Embryo breast muscles were weighed and stained with hematoxylin–eosin after hatching. We found that myofibers grew fast on the 10^th day after hatching (E10) and seldom proliferated beyond the 19^th day after hatching (E19). A total of 6,774 differentially methylated genes (DMGs) were identified between E10 and E19. For RNA-seq data, we found 5,586 differentially expressed genes (DEGs). After overlapping DEGs and DMGs, we recorded 651 shared genes (DEMGs). Subsequently, we performed miRNA-seq analysis and obtained 495 differentially expressed miRNAs (DEMs). Then, we overlapped DEMGs and the target genes of DEMs and obtained 72 overlapped genes (called miRNA-m6A-genes in this study). GO and KEGG results showed DEMGs enriched in many muscle development-related pathways. Furthermore, we chose WNT7B, a key regulator of skeletal muscle development, to perform IGV visualization analysis and found that the m6A levels on the WNT7B gene between E10 and E19 were significantly different. In conclusion, we found that miRNAs, in conjunction with m6A modification, played a key role in the embryonic breast muscle development of Wenchang chickens. The results of this paper offer a theoretical basis for the study of m6A function in muscle development and fat deposition of Wenchang chickens.

MEF2C Expression Is Regulated by the Post-transcriptional Activation of the METTL3-m6A-YTHDF1 Axis in Myoblast Differentiation

N⁶-methyladenosine (m⁶A) plays an essential role in regulating gene expression. However, the effect of m⁶A on skeletal myoblast differentiation and the underlying mechanisms are still unclear. Here, we ascertained mRNA m⁶A methylation exhibited declined changes during bovine skeletal myoblast differentiation, and both MEF2C mRNA expression and m⁶A levels were significantly increased during myoblast differentiation. We found that MEF2C with mutated m⁶A sites significantly inhibited myoblast differentiation compared with wild-type MEF2C. METTL3 promoted MEF2C protein expression through posttranscriptional modification in an m⁶A-YTHDF1-dependent manner. Moreover, MEF2C promoted the expression of METTL3 by binding to its promoter. These results revealed that there is a positive feedback loop between these molecules in myoblast differentiation. Our study provided new insights into skeletal muscle differentiation and fusion, which may provide an RNA methylation-based approach for molecular genetics and breeding in livestock as well as for the treatment of muscle-related diseases.

m6A Methylases Regulate Myoblast Proliferation, Apoptosis and Differentiation

Simple Summary

N⁶-methyladenosine (m⁶A) is the most prevalent methylation modification in eukaryotic mRNA, and it plays an important role in regulating gene expression. Previous studies found that m⁶A methylation plays a role in mammalian skeletal muscle development. Skeletal muscle is an important factor that regulates livestock muscle quality and maintains metabolic homeostasis, and skeletal myogenesis is regulated by a series of transcription factors. However, the role of m⁶A in bovine skeletal myogenesis is unclear. In this study, we examined the expression patterns of the m⁶A methylase genes METTL3, METTL14, WTAP, FTO and ALKBH5 in bovine skeletal muscle tissue and during myogenesis in myoblasts. Furthermore, we used bovine skeletal muscle myoblasts as the object of study to discover the regulatory role of these genes in the process of skeletal myogenesis in vitro. Our findings indicate that these five m⁶A methylases have pronounced and diverse functions in regulating bovine skeletal myoblast proliferation, apoptosis and myogenic differentiation, which can contribute to further understanding the roles of m⁶A in skeletal muscle development.

Abstract

N⁶-methyladenosine (m⁶A) plays an important role in regulating gene expression. Previous studies found that m⁶A methylation affects skeletal muscle development. However, the effect of m⁶A methylases on bovine skeletal myogenesis is still unclear. Here, we found that the expression of m⁶A demethylases (FTO and ALKBH5) was significantly higher in the longissimus dorsi muscle of adult cattle than in newborn cattle. In contrast, the expression of m⁶A methyltransferases (METTL3, METTL14 and WTAP) was reduced. The mRNA expression of all five genes was found to be increased during the myogenesis of myoblasts in vitro. Knockdown of FTO or METTL3 promoted myoblast proliferation, inhibited myoblast apoptosis and suppressed myogenic differentiation, whereas ALKBH5 knockdown had the opposite effect. METTL14 knockdown enhanced myoblast proliferation and impaired myogenic differentiation. WTAP knockdown attenuated proliferation and contributed to apoptosis but did not affect differentiation. Furthermore, the functional domains of these five m⁶A methylases are conserved across species. Our results suggest that m⁶A methylases are involved in regulating skeletal muscle development and that there may be a complex network of m⁶A methylation regulating skeletal myogenesis.