Transcriptome-wide identification and analysis reveals m6A regulation of metabolic reprogramming in shrimp (Marsupenaeus japonicus) under virus infection

It has been reported that the most common post-transcriptional modification of eukaryotic RNA is N6-methyladenosine (m6A). Previous studies show m6A is a key regulator for viral infection and immune response. However, whether there is a pathogen stimulus-dependent m6A regulation in invertebrate shrimp has not been studied. In this study, we performed a transcriptome-wide profiling of mRNA m6A methylation in shrimp (Marsupenaeus japonicus) after white spot syndrome virus (WSSV) infection by methylated RNA immunoprecipitation sequencing (MeRIP-seq). A total of 15,436 m6A peaks were identified in the shrimp, distributed in 8,108 genes, mainly enriched in the CDS, 3′ UTR region and near the stop codon. After WSSV infection, we identified 2,260 m6A peaks with significantly changes, of which 1,973 peaks were significantly up-regulated and 287 peaks were significantly down-regulated. 1,795 genes were identified as differentially methylated genes. GO and KEGG analysis showed that hyper-methylated genes or hypo-methylated genes were highly associated with innate immune process and related to metabolic pathways including HIF-1 signaling pathway, lysine degradation and Wnt signaling pathway. Combined analysis showed a positive correlation between m6A methylation levels and mRNA expression levels. In addition, computational predictions of protein-protein interaction indicated that genes with altered levels of m6A methylation and mRNA expression clustered in metabolism, DNA replication, and protein ubiquitination. ZC3H12A and HIF-1 were two hub genes in protein-protein interaction (PPI) network that involved in immune and metabolism processes, respectively. Our study explored the m6A methylation pattern of mRNA in shrimp after WSSV infection, exhibited the first m6A map of shrimp at the stage of WSSV induced metabolic reprogramming. These findings may reveal the possible mechanisms of m6A-mediated innate immune response in invertebrates.