Yunlong Dang1,2†
Qiao Dong1,2†Bowei Wu1,2Shuhua Yang1,2Jiaming Sun1,2Gengyuan Cui1,2Weixiang Xu1,2Meiling Zhao1,2Yunxuan Zhang1,2
Peng Li1,2*
Lin Li1,2*- 1College of Animal Science and Veterinary Medicine, Shenyang Agricultural University, Shenyang, China
- 2Key Laboratory of Ruminant Infectious Disease Prevention and Control (East), Ministry of Agriculture and Rural Affairs, Beijing, China
Liaoyu white cattle (LYWC) is a local breed in Liaoning Province, China. It has the advantages of grow quickly, high slaughter ratew, high meat quality and strong anti-stress ability. N6 methyladenosine (m6A) is a methylation modification of N6 position of RNA adenine, which is an important modification mechanism affecting physiological phenomena. In this study, we used the longissimus dorsi muscle of LYWC and SIMC for m6A-seq and RNA-seq high-throughput sequencing, and identified the key genes involved in muscle growth and m6A modification development by bioinformatics analysis. There were 31532 m6A peaks in the whole genome of LYWC and 47217 m6A peaks in the whole genome of SIMC. Compared with Simmental cattle group, LYWC group had 17,351 differentially expressed genes: 10,697 genes were up-regulated, 6,654 genes were down regulated, 620 differentially expressed genes were significant, while 16,731 differentially expressed genes were not significant. Among the 620 significantly differentially expressed genes, 295 genes were up-regulated and 325 genes were down regulated. In order to explore the relationship between m6A and mRNA expression in the muscles of LYWC and SIMC, the combined analysis of MeRIP-seq and RNA-seq revealed that 316 genes were m6A modified with mRNA expression. To identify differentially methylated genes related to muscle growth, four related genes were selected for quantitative verification in LYWC and SIMC. GO enrichment and KEGG analysis showed that the differentially expressed genes modified by m6A are mainly involved in skeletal muscle contraction, steroid biosynthesis process, redox process, PPAR pathway and fatty acid metabolism, and galactose metabolism. These results provide a theoretical basis for further research on the role of m6A in muscle growth and development.
Introduction
To date, more than 150 types of posttranscriptional modifications have been identified in the RNA of all living organisms (Boccaletto et al., 2018). The N6-methyladenosine (m6A) modification was discovered in the 1970s and was originally considered to be an abundant nucleotide modification of mRNA in eukaryotic cells (Jia et al., 2013; Yue et al., 2015). Biological functions of m6A modification are mediated by special binding proteins, including methyltransferases, demethylases, and effectors (Zhang et al., 2019). It is involved in various biological processes, such as lipid production and energy metabolism (Zhao et al., 2014; Wang et al., 2015a). In addition, m6A methylation regulates the splicing, expression, decay and translation of mRNA (Wang et al., 2014; Wang et al., 2015b; Xiao et al., 2016). Until recently, little was known about the specific function and mechanism of m6A. Similar to DNA and histone methylation, m6A methylation is also dynamic and reversible in mammals (Wang et al., 2019). It is modulated by several genes, including methyltransferases (METTL3, METTL4 and WTAP) (Liu et al., 2014; Ping et al., 2014), demethylases (ALKBH5 and FTO), (Jia et al., 2011; Zheng et al., 2013) and reading proteins (YTHDF, eIF3 and HNRNPC) (Cao et al., 2016). m6A modification is co installed by a variety of protein complexes (Roundtree et al., 2017). For example, YTHDF2 binds to m6A in mRNA to degrade target genes, while YTHDF1, YTHDF3 and eIF3 promote the translation of m6A containing transcripts (Wang et al., 2014; Wang et al., 2015a; Meyer et al., 2015). As the transferase of writers, METTL3 is composed of catalytic subunit and many other auxiliary subunits. This protein is very important for embryonic growth and development. Embryos lacking METTL3 show pluripotent degradation and damage (Geula et al., 2015). The distribution of mettl3 varies with the type of cell line. In some cases, the change of cell state will lead to the change of its distribution (Knuckles et al., 2017; Xiang et al., 2017). In the cytoplasm, METTL3 itself recognizes the 3′UTR m6A site on mRNA and promotes the formation of translation loop through the interaction with eif3h, so as to promote the protein translation of transcripts (Shi et al., 2019). METTL3 can be functionally regulated by PTM or its protein interaction. It is reported that METTL14 in human cells is phosphorylated at the residue serine399 site on the protein interface with METTL3, indicating that mettl3 has a regulatory function (Wang et al., 2016; Shi et al., 2019). m6A readers protein contains two kinds: one is a stable and direct protein containing YT521-B homology (YTH) domain, and the other is the common RNA binding domain (RBD) (Shi et al., 2019). Both the YTH domain family 1–3 (YTHDF1-3) and the YTH domain containing 1–2 (YTHDC1-2) in humans are stable and directly exercise the reading function of m6A.YTHDF1 and YTHDF3 translation initiation factors promote the translation of target transcripts in cells, and YTHDC2 mediates mRNA stability and translation and regulates cell development (Hsu et al., 2017). The other uses the common RNA binding domain (RBD), such as K homology (KH) domain, RNA recognition motif (RRM) domain and arginine/glycine rich (RGG) domain to preferentially bind the m6A containing region in RNA and exercise the function of m6A reader by regulating the surrounding RNA protein interaction (Shi et al., 2019). Most studies on m6A modification have focused on humans and mice (Dominissini et al., 2012). The m6A methylation is related to obesity (Wang et al., 2020). FTO was the first m6A mRNA demethylase that was discovered. It mediates DNA and RNA demethylation (Jia et al., 2011). The m6A demethylase FTO plays a key role in regulating postpartum growth and energy consumption. A study reported that AMPK regulates lipid accumulation in skeletal muscle cells by regulating FTO expression and FTO-dependent demethylation of m6A (Wu et al., 2017). Research reports on mouse animal models have shown that FTO plays an important role in the regulation of fat mass, adipogenesis, and body weight (Church et al., 2009; Fischer et al., 2009; Church et al., 2010; Gao et al., 2010; McMurray et al., 2013; Merkestein et al., 2015; Ronkainen et al., 2016).
LYWC are excellent beef cattle based on Charolais, breeding the fourth-generation hybrid herd with Liaoning local cattle as the female parent in Liaoning Province, China. A stable population with 93.75% Charolais pedigree and 6.25% local cattle pedigree was formed. LYWC grow quickly, and the slaughter rate was also higher than that of other beef cattle breeds. Due to the large market demand for beef in China, most farms choose LYWC for its excellent production performance. LYWC has wide adaptability, strong stress resistance, and outstanding cold resistance ability and can withstand a low-temperature environment of −30°C. Although LYWC has a better growth rate and slaughter rate, its rough myofiber always influences beef quality directly compared to Simmental and other beef cattle (Jing et al., 2011; Shuangyong et al., 2011). A large number studies have shown that m6A modification plays an important role in regulating lipid production and energy metabolism, inflammatory mechanisms and tumor formation. Based on the necessary functions of m6A modification in regulating gene expression and involving various biological processes, we speculate that m6A modification is involved in beef cattle muscle growth and development. In this study, we aimed to explore the global landscape of differentially expressed m6A methylation genes in muscle tissue between LYWC and SIMC and provide a theoretical basis for further research on the specific regulatory mechanism of unique meat quality and the optimization and selection of LYWC breeds. However, the effect, mechanism, and function of m6A modification on muscle growth and development still needs further research in the future.
Material and Method
Sample Collection and Ethics Statement
Three healthy male Liaoyu white cattles and three Simmental cattles were selected for this study and provided the same feed and drinking water during the breeding period. The breeding environment conditions were identical. The average birth weight of LYWC is 40.0 ± 2.0 kg, the average weight at 6 months is 218 ± 5.0 kg, the average weight at 12 months is 366.8 ± 5.0 kg, and the average weight at 24 months is 624.5 ± 5.0 kg. SIMC have an average birth weight of 41 ± 2.0 kg, an average weight of 200 ± 5.0 kg at 6 months, an average weight of 324 ± 5.0 kg at 12 months, and an average weight of 600 ± 5.0 kg at 24 months. They were raised from birth and slaughtered after 24 months.
The longissimus dorsi muscle samples of two breeds of beef cattle were collected after slaughter. A 1 cm3 muscle sample was taken from the inner side of the spine near the shoulder area. After that, muscle samples were immediately stored in liquid nitrogen. All experimental procedures were approved and performed according to the guidelines of the Laboratory Animal Management Committee of Shenyang Agricultural University.
Experimental Procedure
Total RNA was extracted using TRIzol reagent (Invitrogen, CA, United States). The quality and quantity of total RNA were analyzed by Bioanalyzer 2100 and RNA 6000 Nano Labchip kits (Angelon, CA, United States) with value of RIN >7.0. Oligo-dT magnetic beads were used to enrich total RNA with poly(A) mRNA. Approximately 200 µg of total RNA was subjected to isolation of poly(A) mRNA with poly-T oligo-attached magnetic beads (Invitrogen). The lysed RNA fragments were then incubated with m6A-specific antibodies (No. 202003, Synaptic Systems, Germany) in IP buffer (50 mM Tris-HCl, 750 mM NaCl and 0.5% Igepal CA-630) at 4°C for 2 h with BSA (0.5 μg/μl) (1 ml). The mixture was then incubated with protein A beads and eluted with elution buffer (1×IP buffer and 6.7 mM m6A). The eluted RNA was precipitated with 75% ethanol. According to the chain-specific library prepared by the dUTP method, the eluted m6A fragment (IP) and the unprocessed input control fragment were converted into the final cDNA library. The average insert size of the paired-end library was 100 ± 50 bp. We performed paired-end 2 × 150 bp sequencing on the Illumina NovaSeq™ 6000 platform of LC-BIO Biotech Ltd. (Hangzhou, China) according to the protocol recommended by the supplier.
Bioinformatics Analysis Process
First, Cutadapt and local Perl scripts were used to process the data obtained from sequencing to remove low-quality sequences, contaminated sequences, and linker sequences generated by the sequencer to obtain CleanData (Martin, 2011). Fastp (v0.12.6, data quality control doi: 10.1093/bioinformatics/btp616) was used to verify sequence quality, and HISAT2(v2.0.4, alignment reference sequence: doi: 10.1038/nmeth. 3317)) was used to map the read data to the Bos taurus genome of cattle with default parameters (Bos taurus_NCBI genome version NA) (USDA ARS) (Kim et al., 2015). The threshold settings of differential peak and differential expression are generally | log2Fc | ≥ 1 and P-val < 0.05. At the same time, qval/fdr is corrected for P-val. Exome Peak (v2.13.2, call peak and diff peak: doi: 10.1093/bioinformatics/btt171) read the IP and input data obtained in the experiment (Meng et al., 2014). This program uses bed or bam format files to identify the m6A peak and visualize it in the UCSC genome browser or with IGV software (http://www.igv.org). MEME and HOMER were used to discover known motifs and locate the peak of the obtained motifs using a Perl script (diff_peak |log2FC|≥1, pval<0.05; call peak log2FC ≥ 1, pval<0.05) (Bailey et al., 2009; Heinz et al., 2010). ChIPseeker analyzes the scanned peak calling and annotates the peak genes (Yu et al., 2015). Then, StringTie (v1.3.4d, assembly quantity: doi: 10.1038/NBT. 3122) was used to perform expression operations on all mRNA in the input database to calculate FPKM (FPKM = [total exon fragments/mapped exon readings (million) × exon length (kB)]) (Pertea et al., 2015). Using the R language package edgeR (v3.20.9, difference analysis: doi: 10.1093/bioinformatics/bt p616), differentially expressed mRNAs with log2 (fold change) > 1 or log2 (fold change)<−1 and p-value<0.05 were selected, the reference genome was ARS-UCD1.2 (https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/002/263/795/GCF_002263795.1_ARS -UCD1.2/) (Robinson et al., 2010). The edger input file was raw counts, and we used edger to analyze the PVAL and qval of the results, calculate fpkm values to measure the expression levels of genes, and compare the fold difference obtained by means of fpkm expression compared to the fold change. Functional enrichment we mapped differential gene functional annotations into different GO term/KEGG pathways by writing our own script, embodying the difference test as a hypergeometric test. The integration of MeRIP-seq and RNA-seq data is related through the annotation of peak, and the qualitative correlation is determined through the up/down of the two parts of regulation. Because exomepeak cannot output the quantification of peak level, it cannot calculate the correlation with the expression.
Real-Time Fluorescence Quantitative PCR
We tested four different genes with m6A methylation modification for qRT-PCR analysis, they are related to muscle growth and development (Kee and Hardeman, 2008; Otten et al., 2012; Flix et al., 2013; Siddique et al., 2016). We validated the methylation-modified differential genes and used a qRT-PCR kit (Takara, Dalian, China) to reverse-transcribe the total RNA extracted from the muscle into cDNA. SYBR Green (Vazyme-Q711, China) was used to perform real-time fluorescent quantitative PCR according to the instructions. The ACTB gene was used as an internal reference gene to standardize the expression level of genes. Three trials were performed on three LYWC and three SIM muscle samples. Primer 5 was used to design four pairs of primers, the primer list is shown in Table 1. All primers span the end of the gene. The relative expression of differentially expressed genes was calculated by the 2−△△Ct method. The data are expressed as the mean ± standard error (sample number n = 3). The t-test in SPSS statistical software (version 22.0, Chicago, IL, United States) was used to perform the statistical analyses in the two groups, and the difference was significant when p < 0.05.
TABLE 1. The primer list.
Results and Analysis
Sequence Statistics and Quality Control
First, the raw data generated by sequencing needed to be preprocessed. Cutadapt filtered out unqualified sequences and removed reads with an adaptor, low quality, and unsure base information. The original sequencing volume, effective sequencing volume, Q20, Q30, and GC content were counted, and appropriate evaluation was conducted. Effective data (Clean Data) was prepared for analysis. In the MeRIP-seq library, we obtained two sets of muscle sample data reads. Three biological replicates were performed in each group, and the effective reading data were as follows: LYWC group: 75018648, 63985000 and 73353244; SIM group: 34991588, 42841030 and 50385720. The percentages of valid data (Clean Data) in the two groups of data were 71.85, 71.45, and 73.96% in the LYWC group and 90.24, 91.73, and 96.68% in the SIM group, respectively. In the RNA-seq library, we obtained two sets of muscle sample data reads, each of which was subjected to three biological replicates: the effective read data were as follows: LYWC group: 53535846, 73285174, and 42315990; SIM group: 70936986, 96166358 and 66577584. The percentages of valid data in the two groups of data were 96.52, 96.42, and 96.24% in the LYWC group and 89.09, 98.52, and 90.23% in the SIM group, respectively (Table 2). In Table 2, the proportion of Q20% bases with a quality value ≥20 (sequencing error rate is less than 0.01) and the proportion of Q30% bases with a quality value ≥30 (the sequencing error rate is less than 0.001) are shown.
TABLE 2. Summary of reads quality control.
Map Data to Genome
We used HISAT2 for reference genome comparison of the preprocessed valid data and mapped reads to the Bos taurus cattle (Bos taurus_NCBI, version NA) genome with default parameters. By comparing the obtained reads with the reference genome sequence, we can perform detailed statistics on the data obtained by sequencing and its distribution in the genome. In the m6A-seq library, the IP samples of the longissimus dorsi muscle of beef cattle are LYWC_IP and SIM_IP, and we performed three replicates for each set of samples. The LYWC_IP effective data mapping read rates were 94.13, 94.87, and 95.20%; the SIM_IP effective data mapping read rates were 92.92, 90.97, and 92.30%. In the RNA-seq library, the longissimus dorsi samples of beef cattle are LYWC_input and SIM_input, and we performed three replicates for each set of samples. The effective data mapping read rates of LYWC_input are 96.95, 96.89, and 96.84%, respectively; the effective data mapping read rates of SIM_input are 94.95, 97.79, and 93.68%, respectively. Unique mapped reads are shown in Table 3. According to the regional distribution information of the reference genome, it can be defined as alignment to three parts of exon (exon), intron (intron) and intergenic (intergenic region).In general, the percentage of the sequence alignment to the exon region should be the highest. The results of this experiment showed that the IP samples of Liaoyu white cattle accounted for 97.54, 97.24 and 96.96% in the exon region; the ratios of the input samples were 98.08, 98.10 and 97.79%, respectively. The IP samples of Simmental cattle accounted for 96.81, 97.21 and 96.68% in the exon region; the ratios of the input samples were 96.86, 97.71 and 96.94%, and the results are shown in Figure 1.
TABLE 3. Summary of reads mapped to the cattle reference genome.
FIGURE 1. Comparison of the regional distribution with reference to the genome.
Identification of m6A Modification Sites and Analysis of Differentially Methylated Genes
Use peak-calling software, the R language toolkit exomePeak, was used to scan the m6A peak in the entire genetic dataset. Based on the identification of the IP and input libraries, biological information such as the position and length of the peak in the gene can be obtained. The call peak portion we choose P-val < 0.05, and the diff peak and diff expression portions generally choose|log2 fc| ≥ 1 and P-val < 0.05. We counted and combined all the samples and the degree of enrichment of the reads near the gene transcription start site (TSS). The peaks that could be combined near the TSS are represented in the form of a heat map, as shown in Figure 2A. ChIP seeker software was used to annotate the different peaks and perform GO and KEGG enrichment analyses. In general, the default p < 0.05 was the filter condition of the peak. Compared with Simmental cattle group, we screened 5631 difference peaks in Liaoyu white cattle group, of which 4,059 m6A peaks were significantly up-regulated and 66 m6A peaks were significantly down regulated, as shown in Figure 2B. The distribution of m6A peak in transcripts was analyzed. Analysis of the distribution of m6Apeak in the transcript was performed. We divided the transcript into four parts: 5′-UTR, 3′-UTR, first exon and other exons. It was used to analyze the distribution of different peaks in the original gene function, as shown in Figure 2C. As shown in Table 4, the m6A modification was mainly enriched in the 3′-UTR, and we report the top 20 differences m6Apeak. Difference factor <1 indicates hypomethylation, and difference factor >1 indicates hypermethylation.
FIGURE 2. (A) Heat map of the enrichment of reads near the TSS and the peak distribution near the TSS at the start site of gene transcription. (B) The distribution of peaks of original differences in gene function. (C) Distribution of differential peak on gene functional elements.
TABLE 4. The top 20 differentially expressed m6A peaks.
Motif Analysis
As a dynamic modification phenomenon, RNA methylation is mainly accomplished by the combined action of multiple methylases and methylation binding site motifs. A motif is a nucleotide sequence pattern with biological significance, and the sequence has a high degree of conservation. The methylases involved in the process of RNA methylation recognize the motifs in the gene to generate methylation and regulate gene expression. The motif software MEME was used to search for more credible motifs in the peak area and obtain information about the width, E-value, and location of each motif. We performed motif prediction on each set of samples, and the results are shown in Figure 3. A motif structure that is reported commonly in RNA modifications are RRACH (where R = A or G, H = A, C or U).
FIGURE 3. Sequence showing the motifs with significant differences in muscle samples at the m6A peak.
Whole Gene Analysis and Differential Gene Analysis
The expression level of genes is mainly measured by RPKM (reads per kilobase of exon model per million mapped reads) or FPKM (fragments per kilobase of exon model per million mapped reads) to measure the abundance value of gene expression. In our research, we chose FPKM to report the expression abundance values of different samples of known genes. Compared with Simmental cattle, Liao yu white cattle detected 17,351 differentially expressed genes, 620 genes were significantly different and 16,731 genes were not significantly different (|log2fc| ≥ 1 and p < 0.05). Among the differentially expressed genes, 10,697 genes were upregulated and 6,654 genes were downregulated. Among the 620 significantly differentially expressed genes, 295 were up-regulated and 325 were down regulated. Table 5 shows the top 20 differentially expressed genes we screened. Among the top 20 differentially expressed genes, there are 10 up-regulated genes and 10 down-regulated genes. We used Figures 4A,B to show the gene expression and expression density. We plotted the overall distribution statistics of differentially expressed genes, as shown in Figure 4C. Figure 4D shows the gene heat map between LYWC and SIM samples.
TABLE 5. The top 20 differentially expressed genes.
FIGURE 4. (A) Gene expression cassette diagram. (B) Gene expression density diagram. (C) Differential gene expression volcano diagram. In the figure, log2 of the fold change is the horizontal coordinate, and −log10 (p-value) is used as the vertical coordinate. The horizontal coordinate represents the gene expression in different samples; the vertical coordinate represents the significant difference in gene expression. Among them, red represents upregulated significantly differentially expressed genes, blue represents downregulated significantly differentially expressed genes, and gray represents nonsignificantly differentially expressed genes. (D) LYWC and SIMC gene heat map. Using zscore standardization, the expression levels of genes in different samples can be compared horizontally. From blue to red, the expression amount of genes ranges from low to high.
Joint Analysis of Differentially Expressed Genes and Differentially Expressed Genes
In the entire transcriptome sequencing, we found that there were upregulated and downregulated genes. In the MeRIP-seq sequencing results, according to the changes in abundance, we found that the methylation of the gene itself was upregulated and downregulated. Therefore, we combined the correlation analysis of the two sequencing results to compare and analyze the transcription level and methylation level. In the samples of the LYWC group, 13,624 genes have been modified by m6A, and in the samples of the SIMC group, 24,522 genes have been modified by m6A. We found that among the differentially expressed genes, 620 genes were significantly expressed. Based on this, we screened 316 genes with significant differential expression and m6A methylation modification. The result is shown in Figure 5. Since this experiment mainly explored the regulation of muscle growth and development, we screened four candidate genes related to muscle cell growth and development, as shown in Table 6. The m6A regulation of these genes was upregulated, while gene regulation was downregulated. At the same time, the difference between the two sets of samples were significant.
FIGURE 5. The result obtained by taking the intersection of the gene where the significant difference m6A peak is located and the significant difference expression gene is used, and a stricter screening threshold is used.
TABLE 6. M6A -modified candidate genes related to muscle cell growth and development.
GO Analysis and KEGG Pathway Analysis of Differentially Expressed m6A Methylation Genes
To deeply study the significance of m6A modification in physiological and biochemical processes, we conducted GO (http://www.geneontology.org/) and KEGG (http://www.kegg.jp/) analyses of the different peaks of m6A. The peaks selected were enriched with 898 GO items and 162 pathways. Figure 6A shows the top 25 items in biological processes, the top 15 items in cell components, and the top 10 items in molecular functions. GO analysis (Figure 6B) showed that differentially methylated genes significantly enriched fibers in fat granule tissue, skeletal muscle contraction, and muscle contraction. KEGG pathway analysis (Figure 6C) showed that differentially methylated genes were related to the p53 signaling pathway and PPAR pathway. At the same time, they are also involved in biological processes such as galactose metabolism, fatty acid metabolism, adipocytokine pathway, nitrogen metabolism, arginine synthesis, etc.
FIGURE 6. M6A differential peak gene ontology enrichment analysis and KEGG pathway analysis. (A) Main enrichment 3 and meaningful GO entries of m6Apeak. (B) The first 20 items have significantly enrichment GO terms. (C) The first 20 enriched pathways of the m6A peak.
Verification of Differentially Expressed Genes
To study the function of gene m6A modification and determine the key genes that regulate muscle growth and development in beef cattle muscle cells, We used qRT-PCR for experimental verification. RNA-seq results showed that among the screened differential genes, the expression of Liaoyu white cattle group was lower than that of Simmental cattle group. The qRT-PCR results also confirmed that the m6A-modified gene is indeed present in the Liaoyu white cattle muscle. The trends of these genes are consistent with the RNA-seq results (Figure 7).
FIGURE 7. QPCR results of four different m6A-modified genes in LYWC and SIMC.
Discussion
The modification of m6A is involved in many physiological processes, such as: Mediates mRNA output and synthesis, affects cell maturation, lipogenesis, maintains embryonic development stability, affects cell circadian rhythm, regulates stem cell differentiation, maintains Tregs stability, participates in inflammatory response, apoptosis, muscle production, cell Physiological and biochemical processes such as division. At the same time, the modification of RNA methylation is a dynamic change (Jia et al., 2011), (Zheng et al., 2013)- (Feng et al., 2018).m6A RNA modification is dynamically regulated by methyltransferases (writers) and demethylases (erasers). Since m6A bases cannot be directly detected by sequencing, the dissection of the m6A landscape is hindered; they do not change the base pairing properties and cannot be distinguished from conventional bases by reverse transcription (Brocard et al., 2017). Recently, new methods based on m6A immunoprecipitation or modified selective RNA chemistry methods to isolate modified RNA fragments coupled with high-throughput sequencing, namely, m6A-seq and MeRIP-seq, have identified thousands of hundred-nucleotide fragments containing modifications in the transcriptomes of mammalian cells (Dominissini et al., 2012; Meyer et al., 2012). Modification of m6A has been successively discovered in many animals, plants, bacteria and other microorganisms.
Based on numerous research reports, it is found that the majority of cows mammary gland and lymphocytes undergo m6A modification phenomenon (Burtseva et al., 1979; Horowitz et al., 1984; Wu et al., 2021). However, there are still few reports on beef cattle. We speculate that m6A is involved in the muscle growth and development process of beef cattle. Our data show that there are a large number of methylation modifications during muscle growth and development. It may have an important effect on the types of muscle fibers, the maturation of muscle cells, and the changes in muscle structure.
Through laser-induced microdamage of zebrafish muscles combined with cell repair, it was found that the XIRP1 gene is abundant in skeletal muscle and involved in cell repair, cells and new myofibrils, and the repair of damage does not involve cell proliferation (Otten et al., 2012). Troponin T (TNNT1) exists as a group of homologous proteins in the striated muscle of vertebrates and invertebrates. Mutations in the TNNT1 gene can cause rod-shaped myopathy. From animal model experiments, it was found that lack of TNNT1 reduced the content of slow fibers, accompanied by type II fiber hypertrophic growth and increased muscle fatigue (Feng et al., 2009; Wei et al., 2014). Myosin 2 (MYOM2) is the main component of the myofibril M-band of the sarcomere and the central gene in the interaction of sarcomere genes (Auxerre-Plantie et al., 2020). Research by Auxerre-Plantié found that loss of function and moderate knockdown of this gene can lead to myocardial expansion, and severe knockdown can lead to increased sarcomeric myosin (Auxerre-Plantie et al., 2020). Research by Andrei found that in hypothyroid rats, MYOM2 expression increased 3.4 times. Through small-molecule interference of RNA with MYOM2, it was found that the contraction speed of cardiomyocytes were severely reduced (Rozanski et al., 2013). Myosin heavy chain 6 (MYH6) is widely found mainly in the heart and smooth muscle. This gene is mainly expressed in type I fibers. The presence or absence of MYH6 and its family gene MYH7 determines the slow or fast-twitch phenotype of skeletal muscle (Stuart et al., 2016).
Based on GO enrichment and KEGG pathway analysis, we speculate that m6A modification in genes has potentially important functions and may play a vital role in certain pathways are involved in cell growth and development. Peroxisome proliferator-activated receptors (PPARs) are nuclear hormone receptors activated by fatty acids and their derivatives. They are ligand-activated receptors in the nuclear hormone receptor family. Three subtypes have been found in different species, which control many intracellular metabolic processes. The subtypes include PPARα (also known as NR1C1). PPARα participates in the liver and skeletal muscle through regulation and expression of lipid metabolism genes. PPARβ/δ participates in lipid oxidation and cell proliferation. PPARγ promotes the differentiation of adipocytes to enhance blood glucose uptake. PPAR transcriptional activity can be regulated by nongene crosstalk with phosphatases and kinases, including ERK1/2, p38-MAPK, PKA, PKC, AMPK and GSK3. At the same time, nuclear receptor coactivator (coactivator) and PPAR-RXR act synergistically and complement and stabilize the active transcription complex, which can regulate lipid metabolism and fat formation, maintain metabolic homeostasis and the expression of inflammation genes, and induce anticancer effects in a variety of human tumors.
The m6A regulation level of the differentially expressed genes screened in this study was negatively correlated with the transcription level. The RNA-seq results showed that the differential gene expression in LYWC was lower than that in SIM. The results of qRT-PCR showed that the differentially expressed genes for m6A methylation were all present in the muscle tissue of beef cattle. Therefore, this indicates that m6A modification not only participates in the process of muscle growth and development but may also regulate mRNA degradation.
Skeletal muscle development is a complex biological process. The regulatory role of myogenic regulatory factors and the study of apparent modifications, including DNA methylation and histone modification, in the regulation of skeletal muscle development have given us a preliminary understanding of the regulatory network of skeletal muscle development. Based on the involvement of m6A in the regulation of mouse brain development, fat formation, and other tissue development processes, we speculate that m6A is also involved in the regulation of skeletal muscle development (Zhao et al., 2014; Ma et al., 2018; Wang et al., 2018). Studies have shown that the regulation of METTL3 gene expression and regulation of m6A levels in myoblasts affect the differentiation process of myocytes and the expression of key regulatory genes (Kudou et al., 2017; Chen et al., 2019). This shows that m6A is involved in the regulation of muscle cell differentiation. During the development of animal embryos and the growth and development of brain tissue after birth, neural stem cells are required for differentiation and self-renewal. Studies have shown that knocking out METTL14 in mouse embryos will interrupt the cycle of radial glial cells in the nerves, which will eventually lead to a decrease in the thickness of the cerebral cortex and even death after birth. Overexpression and specific knockout of the FTO and METTL3 genes in pig adipocytes revealed that FTO expression levels are negatively correlated with m6A levels and positively correlated with adipogenesis, while METTL3 expression levels are positively correlated with m6A levels and negatively correlated with adipogenesis (Wang et al., 2015b). Mice lacking FTO function experience increased energy expenditure, growth retardation, lean body size after birth and deformity (Boissel et al., 2009). Through transcriptome sequencing of the muscle tissues of three different breeds of wild boar, Landrace pig and Rongchang pig, a complete transcriptome map of m6A was drawn. It was found that m6A is widely distributed in muscle tissue, and m6A is mainly enriched in related gene stop codons, 3′UTRs, and protein coding-regions. In addition, data show that there is a clear m6A peak around the stop codon of the cAMP response element-binding protein CREB and zinc finger protein ZNF-related genes, indicating that m6A is enriched here. CREB was first discovered as a transcriptional regulator of cell metabolism that regulates the cAMP response. It is an important gene regulating somatostatin. ZNF has also been considered one of the most important eukaryotic transcription factors and plays an important role in gene regulation.
In conclusion, this study analyzed m6A methylation modification in the muscle tissue of Liaoyu white cattle and Simmental cattle. Based on our experimental results, we speculate that m6A modification plays an important role in muscle growth and development. This study shows that TNNT1, XIRP1, MYOM2, and MYH6 are likely to play a key role in muscle growth and muscle differentiation. In addition, the data obtained through high-throughput sequencing provide a theoretical basis for further exploring the function of m6A modification on muscle growth and development. At the same time, the regulatory mechanism of m6A modification in muscle still needs to be studied in depth in the future.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI [accession: PRJNA778440].
Ethics Statement
The animal study was reviewed and approved by Laboratory Animal Management Committee of Shenyang Agricultural University.
Author Contributions
PL and LL conceived and designed the study. YD and QD performed the main experiments. BW, JS, GC, and WX Sample collection. MZ, YZ, and SY performed data analysis. YD, QD, PL, and LL drafted and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This research was funded by the National Natural Science Foundation of China (grants No. 31872538) and Scientific Research Funding Project of Liaoning Province (No. 2021JH1/10400033), Education Department of Liaoning Province (No. LSNFW201901, LSNJC202012). Liaoning Provincial joint fund for innovation capability improvement (2021-NLTS-11-05).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Auxerre-Plantie, E., Nielsen, T., Grunert, M., Olejniczak, O., Perrot, A., Özcelik, C., et al. (2020). Identification of MYOM2 as a Candidate Gene in Hypertrophic Cardiomyopathy and Tetralogy of Fallot and its Functional Evaluation in the Drosophila Heart. Dis. Model. Mech. 13, dmm045377. doi:10.1242/dmm.045377
Bailey, T. L., Boden, M., Buske, F. A., Frith, M., Grant, C. E., Clementi, L., et al. (2009). MEME SUITE: Tools for Motif Discovery and Searching. Nucleic Acids Res. 37, W202. doi:10.1093/nar/gkp335
Boccaletto, P., Machnicka, M. A., Purta, E., Piątkowski, P., Bagiński, B., Wirecki, T. K., et al. (2018). MODOMICS: a Database of RNA Modification Pathways. 2017 Update. Nucleic Acids Res. 46 (D1), D303–D307. doi:10.1093/nar/gkx1030
Boissel, S., Reish, O., Proulx, K., Kawagoe-Takaki, H., Sedgwick, B., Yeo, G. S. H., et al. (2009). Loss-of-function Mutation in the Dioxygenase-Encoding FTO Gene Causes Severe Growth Retardation and Multiple Malformations. Am. J. Hum. Genet. 85 (1), 106–111. doi:10.1016/j.ajhg.2009.06.002
Brocard, M., Ruggieri, A., and Locker, N. (2017). m6A RNA Methylation, a New Hallmark in Virus-Host Interactions. J. Gen. Virol. 98 (9), 2207–2214. doi:10.1099/jgv.0.000910
Burtseva, N. N., Azizov, Iu. M., and Vaniushin, B. F. (1979). Tissue Specificity of the Decrease of Cattle Lymphocyte DNA Methylation during Chronic Lymphoid Leukemia. Biokhimiia 44 (7), 1296–1302.
Cao, G., Li, H.-B., Yin, Z., and Flavell, R. A. (2016). Recent Advances in Dynamic M 6 A RNA Modification. Open Biol. 6 (4), 160003. doi:10.1098/rsob.160003
Chen, J. N., Chen, Y., Wei, Y. Y., Raza, M. A., Zou, Q., Xi, X. Y., et al. (2019). Regulation of M^(6)A RNA Methylation and its Effect on Myogenic Differentiation in Murine Myoblasts. Mol. Biol. (Mosk) 53 (3), 436–445. doi:10.1134/S0026898419030042
Church, C., Lee, S., Bagg, E. A. L., McTaggart, J. S., Deacon, R., Gerken, T., et al. (2009). A Mouse Model for the Metabolic Effects of the Human Fat Mass and Obesity Associated FTO Gene. Plos Genet. 5 (8), e1000599. doi:10.1371/journal.pgen.1000599
Church, C., Moir, L., McMurray, F., Girard, C., Banks, G. T., Teboul, L., et al. (2010). Overexpression of Fto Leads to Increased Food Intake and Results in Obesity. Nat. Genet. 42 (12), 1086–1092. doi:10.1038/ng.713
Dominissini, D., Moshitch-Moshkovitz, S., Schwartz, S., Salmon-Divon, M., Ungar, L., Osenberg, S., et al. (2012). Topology of the Human and Mouse m6A RNA Methylomes Revealed by m6A-Seq. Nature 485 (7397), 201–206. doi:10.1038/nature11112
Feng, H.-Z., Wei, B., and Jin, J.-P. (2009). Deletion of a Genomic Segment Containing the Cardiac Troponin I Gene Knocks Down Expression of the Slow Troponin T Gene and Impairs Fatigue Tolerance of Diaphragm Muscle. J. Biol. Chem. 284 (46), 31798–31806. doi:10.1074/jbc.m109.020826
Feng, Z., Li, Q., Meng, R., Yi, B., and Xu, Q. (2018). METTL 3 Regulates Alternative Splicing of MyD88 upon the Lipopolysaccharide‐induced Inflammatory Response in Human Dental Pulp Cells. J. Cel Mol Med 22 (5), 2558–2568. doi:10.1111/jcmm.13491
Fischer, J., Koch, L., Emmerling, C., Vierkotten, J., Peters, T., Brüning, J. C., et al. (2009). Inactivation of the Fto Gene Protects from Obesity. Nature 458 (7240), 894–898. doi:10.1038/nature07848
Flix, B., de la Torre, C., Castillo, J., Casal, C., Illa, I., and Gallardo, E. (2013). Dysferlin Interacts with Calsequestrin-1, Myomesin-2 and Dynein in Human Skeletal Muscle. Int. J. Biochem. Cel Biol. 45 (8), 1927–1938. doi:10.1016/j.biocel.2013.06.007
Gao, X., Shin, Y.-H., Li, M., Wang, F., Tong, Q., and Zhang, P. (2010). The Fat Mass and Obesity Associated Gene FTO Functions in the Brain to Regulate Postnatal Growth in Mice. PLoS One 5 (11), e14005. doi:10.1371/journal.pone.0014005
Geula, S., Moshitch-Moshkovitz, S., Dominissini, D., Mansour, A. A., Kol, N., Salmon-Divon, M., et al. (2015). m 6 A mRNA Methylation Facilitates Resolution of Naïve Pluripotency toward Differentiation. Science 347 (6225), 1002–1006. doi:10.1126/science.1261417
Heinz, S, Benner, C., Spann, N., Bertolino, E., Lin, Y. C., Laslo, P., et al. (2010). Simple Combinations of Lineage-Determining Transcription Factors Prime Cis-Regulatory Elements Required for Macrophage and B Cell Identities - ScienceDirect. Mol. Cel 38 (4), 576–589. doi:10.1016/j.molcel.2010.05.004
Horowitz, S., Horowitz, A., Nilsen, T. W., Munns, T. W., and Rottman, F. M. (1984). Mapping of N6-Methyladenosine Residues in Bovine Prolactin mRNA. Proc. Natl. Acad. Sci. 81 (18), 5667–5671. doi:10.1073/pnas.81.18.5667
Hsu, P. J., Zhu, Y., Ma, H., Guo, Y., Shi, X., Liu, Y., et al. (2017). Ythdc2 Is an N6-Methyladenosine Binding Protein that Regulates Mammalian Spermatogenesis. Cell Res 27 (9), 1115–1127. doi:10.1038/cr.2017.99
Jia, G., Fu, Y., and He, C. (2013). Reversible RNA Adenosine Methylation in Biological Regulation. Trends Genet. 29 (2), 108–115. doi:10.1016/j.tig.2012.11.003
Jia, G., Fu, Y., Zhao, X., Dai, Q., Zheng, G., Yang, Y., et al. (2011). N6-methyladenosine in Nuclear RNA Is a Major Substrate of the Obesity-Associated FTO. Nat. Chem. Biol. 7 (12), 885–887. doi:10.1038/nchembio.687
Jing, L., Huaiye, L., and Maobin, Q. (2011). A Contrast Experimental Report about slaughter Performance of Liaoyu White Cattle and 5 Varieties (Hybrid) Cattle. Mod. Anim. Husbandry Vet. (2), 30–33. doi:10.3969/j.issn.1672-9692.2011.02.011
Kee, A. J., and Hardeman, E. C. (2008). Tropomyosins in Skeletal Muscle Diseases. Adv. Exp. Med. Biol. 644, 143–157. doi:10.1007/978-0-387-85766-4_12
Kim, D., Langmead, B., and Salzberg, S. L. (2015). HISAT: A Fast Spliced Aligner with Low Memory Requirements. Nat. Methods 12 (4), 357–360. doi:10.1038/nmeth.3317
Knuckles, P., Carl, S. H., Musheev, M., Niehrs, C., Wenger, A., and Bühler, M. (2017). RNA Fate Determination through Cotranscriptional Adenosine Methylation and Microprocessor Binding. Nat. Struct. Mol. Biol. 24 (7), 561–569. doi:10.1038/nsmb.3419
Kudou, K., Komatsu, T., Nogami, J., Maehara, K., Harada, A., Saeki, H., et al. (2017). The Requirement of Mettl3-Promoted MyoD mRNA Maintenance in Proliferative Myoblasts for Skeletal Muscle Differentiation. Open Biol. 7 (9), 170119. doi:10.1098/rsob.170119
Liu, J., Yue, Y., Han, D., Wang, X., Fu, Y., Zhang, L., et al. (2014). A METTL3-METTL14 Complex Mediates Mammalian Nuclear RNA N6-Adenosine Methylation. Nat. Chem. Biol. 10 (2), 93–95. doi:10.1038/nchembio.1432
Ma, C., Chang, M., Lv, H., Zhang, Z.-W., Zhang, W., He, X., et al. (2018). RNA m6A Methylation Participates in Regulation of Postnatal Development of the Mouse Cerebellum. Genome Biol. 19 (1), 68. doi:10.1186/s13059-018-1435-z
Martin, M. (2011). Cutadapt Removes Adapter Sequences from High-Throughput Sequencing Reads. Embnet J. 17 (1). doi:10.14806/ej.17.1.200
McMurray, F., Church, C. D., Larder, R., Nicholson, G., Wells, S., Teboul, L., et al. (2013). Adult Onset Global Loss of the Fto Gene Alters Body Composition and Metabolism in the Mouse. Plos Genet. 9 (1), e1003166. doi:10.1371/journal.pgen.1003166
Meng, J., Lu, Z., Liu, H., Zhang, L., Zhang, S., Chen, Y., et al. (2014). A Protocol for RNA Methylation Differential Analysis with MeRIP-Seq Data and exomePeak R/Bioconductor Package. Methods A Companion Methods Enzymol. 69, 274. doi:10.1016/j.ymeth.2014.06.008
Merkestein, M., Laber, S., McMurray, F., Andrew, D., Sachse, G., Sanderson, J., et al. (2015). FTO Influences Adipogenesis by Regulating Mitotic Clonal Expansion. Nat. Commun. 6, 6792. doi:10.1038/ncomms7792
Meyer, K. D., Patil, D. P., Zhou, J., Zinoviev, A., Skabkin, M. A., Elemento, O., et al. (2015). 5′ UTR m6A Promotes Cap-independent Translation. Cell 163 (4), 999–1010. doi:10.1016/j.cell.2015.10.012
Meyer, K. D., Saletore, Y., Zumbo, P., Elemento, O., Mason, C. E., and Jaffrey, S. R. (2012). Comprehensive Analysis of mRNA Methylation Reveals Enrichment in 3′ UTRs and Near Stop Codons. Cell 149 (7), 1635–1646. doi:10.1016/j.cell.2012.05.003
Otten, C., van der Ven, P. F., Lewrenz, I., Paul, S., Steinhagen, A., Busch-Nentwich, E., et al. (2012). Xirp Proteins Mark Injured Skeletal Muscle in Zebrafish. PLoS One 7 (2), e31041. doi:10.1371/journal.pone.0031041
Pertea, M., Pertea, G. M., Antonescu, C. M., Chang, T.-C., Mendell, J. T., and Salzberg, S. L. (2015). StringTie Enables Improved Reconstruction of a Transcriptome from RNA-Seq Reads. Nat. Biotechnol. 33 (3), 290–295. doi:10.1038/nbt.3122
Ping, X.-L., Sun, B.-F., Wang, L., Xiao, W., Yang, X., Wang, W.-J., et al. (2014). Mammalian WTAP Is a Regulatory Subunit of the RNA N6-Methyladenosine Methyltransferase. Cel Res 24 (2), 177–189. doi:10.1038/cr.2014.3
Robinson, M. D., Mccarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data. Bioinformatics 26, 139. doi:10.1093/bioinformatics/btp616
Ronkainen, J., Mondini, E., Cinti, F., Cinti, S., Sebért, S., Savolainen, M. J., et al. (2016). Fto-Deficiency Affects the Gene and MicroRNA Expression Involved in Brown Adipogenesis and Browning of White Adipose Tissue in Mice. Int. J. Mol. Sci. 17 (11). doi:10.3390/ijms17111851
Roundtree, I. A., Evans, M. E., Pan, T., and He, C. (2017). Dynamic RNA Modifications in Gene Expression Regulation. Cell 169 (7), 1187–1200. doi:10.1016/j.cell.2017.05.045
Rozanski, A., Takano, A. P. C., Kato, P. N., Soares, A. G., Lellis-Santos, C., Campos, J. C., et al. (2013). M-protein Is Down-Regulated in Cardiac Hypertrophy Driven by Thyroid Hormone in Rats. Mol. Endocrinol. 27 (12), 2055–2065. doi:10.1210/me.2013-1018
Shi, H., Wei, J., and He, C. (2019). Where, when, and How: Context-dependent Functions of RNA Methylation Writers, Readers, and Erasers. Mol. Cel. 74 (4), 640–650. doi:10.1016/j.molcel.2019.04.025
Shuangyong, J., Huaiye, L., and Maobin, Q. (2011). A Contrast Experimental Report about Beef Quality of Liaoyu White Cattle and 5 Varieties (Hybrid) Cattle. Mod. J. Anim. Husbandry Vet. Med 2011 (3), 22–25. doi:10.3969/j.issn.1672-9692.2011.03.009
Siddique, B. S., Kinoshita, S., Wongkarangkana, C., Asakawa, S., and Watabe, S. (2016). Evolution and Distribution of Teleost myomiRNAs: Functionally Diversified myomiRs in Teleosts. Mar. Biotechnol. 18 (3), 436–447. doi:10.1007/s10126-016-9705-9
Stuart, C. A., Stone, W. L., Howell, M. E. A., Brannon, M. F., Hall, H. K., Gibson, A. L., et al. (2016). Myosin Content of Individual Human Muscle Fibers Isolated by Laser Capture Microdissection. Am. J. Physiology-Cell Physiol. 310 (5), C381–C389. doi:10.1152/ajpcell.00317.2015
Wang, C.-X., Cui, G.-S., Liu, X., Xu, K., Wang, M., Zhang, X.-X., et al. (2018). METTL3-mediated m6A Modification Is Required for Cerebellar Development. Plos Biol. 16 (6), e2004880. doi:10.1371/journal.pbio.2004880
Wang, J.-y., Chen, L.-j., and Qiang, P. (2020). The Potential Role of N6-Methyladenosine (m6A) Demethylase Fat Mass and Obesity-Associated Gene (FTO) in Human Cancers. Ott 13, 12845–12856. doi:10.2147/ott.s283417
Wang, X., Feng, J., Xue, Y., Guan, Z., Zhang, D., Liu, Z., et al. (2016). Structural Basis of N6-Adenosine Methylation by the METTL3-METTL14 Complex. Nature 534 (7608), 575–578. doi:10.1038/nature18298
Wang, X., Lu, Z., Gomez, A., Hon, G. C., Yue, Y., Han, D., et al. (2014). N6-methyladenosine-dependent Regulation of Messenger RNA Stability. Nature 505 (7481), 117–120. doi:10.1038/nature12730
Wang, X., Zhao, B. S., Roundtree, I. A., Lu, Z., Han, D., Ma, H., et al. (2015). N6-methyladenosine Modulates Messenger RNA Translation Efficiency. Cell 161 (6), 1388–1399. doi:10.1016/j.cell.2015.05.014
Wang, X., Zhu, L., Chen, J., and Wang, Y. (2015). mRNA m6A Methylation Downregulates Adipogenesis in Porcine Adipocytes. Biochem. Biophysical Res. Commun. 459 (2), 201–207. doi:10.1016/j.bbrc.2015.02.048
Wang, Y., Zheng, Y., Guo, D., Zhang, X., Guo, S., Hui, T., et al. (2019). m6A Methylation Analysis of Differentially Expressed Genes in Skin Tissues of Coarse and Fine Type Liaoning Cashmere Goats. Front. Genet. 10, 1318. doi:10.3389/fgene.2019.01318
Wei, B., Lu, Y., and Jin, J.-P. (2014). Deficiency of Slow Skeletal Muscle Troponin T Causes Atrophy of Type I Slow Fibres and Decreases Tolerance to Fatigue. J. Physiol. 592 (6), 1367–1380. doi:10.1113/jphysiol.2013.268177
Wu, K., Jia, S., Zhang, J., Zhang, C., Wang, S., Rajput, S. A., et al. (2021). Transcriptomics and Flow Cytometry Reveals the Cytotoxicity of Aflatoxin B1 and Aflatoxin M1 in Bovine Mammary Epithelial Cells. Ecotoxicology Environ. Saf. 209, 111823. doi:10.1016/j.ecoenv.2020.111823
Wu, W., Feng, J., Jiang, D., Zhou, X., Jiang, Q., Cai, M., et al. (2017). AMPK Regulates Lipid Accumulation in Skeletal Muscle Cells through FTO-dependent Demethylation of N6-Methyladenosine. Sci. Rep. 7, 41606. doi:10.1038/srep41606
Xiang, Y., Laurent, B., Hsu, C.-H., Nachtergaele, S., Lu, Z., Sheng, W., et al. (2017). RNA m6A Methylation Regulates the Ultraviolet-Induced DNA Damage Response. Nature 543 (7646), 573–576. doi:10.1038/nature21671
Xiao, W., Adhikari, S., Dahal, U., Chen, Y.-S., Hao, Y.-J., Sun, B.-F., et al. (2016). Nuclear M 6 A Reader YTHDC1 Regulates mRNA Splicing. Mol. Cel 61 (4), 507–519. doi:10.1016/j.molcel.2016.01.012
Yu, G., Wang, L. G., and He, Q. Y. (2015). ChIPseeker: an R/Bioconductor Package for ChIP Peak Annotation, Comparison and Visualization. Bioinformatics 31 (14), 2382–2383. doi:10.1093/bioinformatics/btv145
Yue, Y., Liu, J., and He, C. (2015). RNA N6-Methyladenosine Methylation in post-transcriptional Gene Expression Regulation. Genes Dev. 29 (13), 1343–1355. doi:10.1101/gad.262766.115
Zhang, C., Fu, J., and Zhou, Y. (2019). A Review in Research Progress Concerning m6A Methylation and Immunoregulation. Front. Immunol. 10, 922. doi:10.3389/fimmu.2019.00922
Zhao, X., Yang, Y., Sun, B.-F., Shi, Y., Yang, X., Xiao, W., et al. (2014). FTO-dependent Demethylation of N6-Methyladenosine Regulates mRNA Splicing and Is Required for Adipogenesis. Cel Res 24 (12), 1403–1419. doi:10.1038/cr.2014.151
Keywords: m6A methylation, RNA-seq, muscle growth and development, genetic modification, species
Citation: Dang Y, Dong Q, Wu B, Yang S, Sun J, Cui G, Xu W, Zhao M, Zhang Y, Li P and Li L (2022) Global Landscape of m6A Methylation of Differently Expressed Genes in Muscle Tissue of Liaoyu White Cattle and Simmental Cattle. Front. Cell Dev. Biol. 10:840513. doi: 10.3389/fcell.2022.840513
Received: 21 December 2021; Accepted: 22 February 2022;
Published: 10 March 2022.
Edited by:
Krzysztof Flisikowski, Technical University of Munich, GermanyReviewed by:
Haojie Zhang, Guangxi University, ChinaJiangbo Wei, University of Chicago, United States
Copyright © 2022 Dang, Dong, Wu, Yang, Sun, Cui, Xu, Zhao, Zhang, Li and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Peng Li, Lipeng2018@syau.edu.cn; Lin Li, Lilin619619@syau.edu.cn
†These authors have contributed equally to this work