Dynamic N6-methyladenosine RNA methylation landscapes reveal epi-transcriptomic modulation induced by ammonia nitrogen exposure in the Pacific whiteleg shrimp Litopenaeus vannamei

Abstract

Despite the versatility of RNA m⁶A methylation in regulating various biological processes, its involvement in the physiological response to ammonia nitrogen toxicity in decapod crustaceans like shrimp remains enigmatic. Here, we provided the first characterization of dynamic RNA m⁶A methylation landscapes induced by toxic ammonia exposure in the Pacific whiteleg shrimp Litopenaeus vannamei. The global m⁶A methylation level showed significant decrease following ammonia exposure, and most of the m⁶A methyltransferases and m⁶A binding proteins were significantly repressed. Distinct from many well-studied model organisms, m⁶A methylated peaks in the transcriptome of L. vannamei were enriched not only near the termination codon and in the 3' untranslated region (UTR), but also around the start codon and in the 5' UTR. Upon ammonia exposure, 11,430 m⁶A peaks corresponding to 6113 genes were hypo-methylated, and 5660 m⁶A peaks from 3912 genes were hyper-methylated. The differentially methylated genes showing significant changes in expression were over-represented by genes associated with metabolism, cellular immune defense and apoptotic signaling pathways. Notably, the m⁶A-modified ammonia-responsive genes encompassed a subset of genes related to glutamine synthesis, purine conversion and urea production, implying that m⁶A methylation may modulate shrimp ammonia stress responses partly through these ammonia metabolic processes.

Keywords: Ammonia nitrogen; Decapod crustacean; M(6)A; RNA methylation; Toxicity.