Global Landscape of m6A Methylation of Differently Expressed Genes in Muscle Tissue of Liaoyu White Cattle and Simmental Cattle

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. N(6) methyladenosine (m6A) is a methylation modification of N(6) 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.