Hidden codes in mRNA: Control of gene expression by m6A.

Information in mRNA has largely been thought to be confined to its nucleotide sequence. However, the advent of mapping techniques to detect modified nucleotides has revealed that mRNA contains additional information in the form of chemical modifications. The most abundant modified nucleotide is N6-methyladenosine (m6A), a methyl modification of adenosine. Although early studies viewed m6A as a dynamic and tissue-specific modification, it is now clear that the mRNAs that contain m6A and the location of m6A in those transcripts are largely universal and are influenced by gene architecture, i.e., the size and location of exons and introns. m6A can affect nuclear processes such as splicing and epigenetic regulation, but the major effect of m6A on mRNAs is to promote degradation in the cytoplasm. m6A marks a functionally related cohort of mRNAs linked to certain biological processes, including cell differentiation and cell fate determination. m6A is also enriched in other cohorts of mRNAs and can therefore affect their respective cellular processes and pathways. Future work will focus on understanding how the m6A pathway is regulated to achieve control of m6A-containing mRNAs.

KIAA1429/VIRMA promotes breast cancer progression by m6 A-dependent cytosolic HAS2 stabilization.

N6 -methyladenosine (m6 A), the most abundant internal modification in eukaryotic mRNA, plays important roles in many physiological and pathological processes, including the development and progression of cancer. RNA modification by m6 A is regulated by methyltransferases, demethylases, and m6 A-binding proteins that function in large part by regulating mRNA expression and function. Here, we investigate the expression of m6 A regulatory proteins in breast cancer. We find that expression of KIAA1429/VIRMA, a component of the m6 A methyltransferase complex, is upregulated in breast cancer tissue and correlates positively with poor survival. KIAA1429/VIRMA is mislocalized to the cytosol of breast cancer tissues and cell lines, and shRNA-mediated knockdown inhibits breast cancer cell proliferation, migration, and invasion. Mechanistically, KIAA1429/VIRMA is shown to bind to the m6 A-dependent RNA-binding protein insulin-like growth factor 2 mRNA-binding protein 3 (IGF2BP3), leading to recruitment and stabilization of m6 A-modified hyaluronan synthase 2 (HAS2) mRNA. HAS2 mRNA and KIAA1429/VIRMA mRNA levels correlate positively in breast cancer tissues, suggesting that the KIAA1429/VIRMA-IGF2BP3-HAS2 axis promotes breast cancer growth and contributes to poor prognosis.

VIRMA promotes nasopharyngeal carcinoma, tumorigenesis, and metastasis by upregulation of E2F7 in an m6A-dependent manner.

The N6-methyladenosine (m6A) modification possesses new and essential roles in tumor initiation and progression by regulating mRNA biology. However, the role of aberrant m6A regulation in nasopharyngeal carcinoma (NPC) remains unclear. Here, through comprehensive analyses of NPC cohorts from the GEO database and our internal cohort, we identified that VIRMA, an m6A writer, is significantly upregulated in NPC and plays an essential role in tumorigenesis and metastasis of NPC, both in vitro and in vivo. High VIRMA expression served as a prognostic biomarker and was associated with poor outcomes in patients with NPC. Mechanistically, VIRMA mediated the m6A methylation of E2F7 3'-UTR, then IGF2BP2 bound, and maintained the stability of E2F7 mRNA. An integrative high-throughput sequencing approach revealed that E2F7 drives a unique transcriptome distinct from the classical E2F family in NPC, which functioned as an oncogenic transcriptional activator. E2F7 cooperated with CBFB-recruited RUNX1 in a non-canonical manner to transactivate ITGA2, ITGA5, and NTRK1, strengthening Akt signaling-induced tumor-promoting effect.

Overexpression of VIRMA confers vulnerability to breast cancers via the m6A-dependent regulation of unfolded protein response.

Virilizer-like m6A methyltransferase-associated protein (VIRMA) maintains the stability of the m6A writer complex. Although VIRMA is critical for RNA m6A deposition, the impact of aberrant VIRMA expression in human diseases remains unclear. We show that VIRMA is amplified and overexpressed in 15-20% of breast cancers. Of the two known VIRMA isoforms, the nuclear-enriched full-length but not the cytoplasmic-localised N-terminal VIRMA promotes m6A-dependent breast tumourigenesis in vitro and in vivo. Mechanistically, we reveal that VIRMA overexpression upregulates the m6A-modified long non-coding RNA, NEAT1, which contributes to breast cancer cell growth. We also show that VIRMA overexpression enriches m6A on transcripts that regulate the unfolded protein response (UPR) pathway but does not promote their translation to activate the UPR under optimal growth conditions. Under stressful conditions that are often present in tumour microenvironments, VIRMA-overexpressing cells display enhanced UPR and increased susceptibility to death. Our study identifies oncogenic VIRMA overexpression as a vulnerability that may be exploited for cancer therapy.

HNRNPC promotes collagen fiber alignment and immune evasion in breast cancer via activation of the VIRMA-mediated TFAP2A/DDR1 axis.

BACKGROUND: Cancers aggressively reorganize collagen in their microenvironment, leading to the evasion of tumor cells from immune surveillance. However, the biological significance and molecular mechanism of collagen alignment in breast cancer (BC) have not been well established. METHODS: In this study, BC-related RNA-Seq data were obtained from the TCGA database to analyze the correlation between DDR1 and immune cells. Mouse BC cells EO771 were selected for in vitro validation, and dual-luciferase experiments were conducted to examine the effect of TFAP2A on DDR1 promoter transcription activity. ChIP experiments were performed to assess TFAP2A enrichment on the DDR1 promoter, while Me-RIP experiments were conducted to detect TFAP2A mRNA m6A modification levels, and PAR-CLIP experiments were conducted to determine VIRMA's binding to TFAP2A mRNA and RIP experiments to investigate HNRNPC's recognition of m6A modification on TFAP2A mRNA. Additionally, an in vivo mouse BC transplant model and the micro-physiological system was constructed for validation, and Masson staining was used to assess collagen fiber arrangement. Immunohistochemistry was conducted to identify the number of CD8-positive cells in mouse BC tumors and Collagen IV content in ECM, while CD8 + T cell migration experiments were performed to measure CD8 + T cell migration. RESULTS: Bioinformatics analysis showed that DDR1 was highly expressed in BC and negatively correlated with the proportion of anti-tumor immune cell infiltration. In vitro cell experiments indicated that VIRMA, HNRNPC, TFAP2A, and DDR1 were highly expressed in BC cells. In addition, HNRNPC promoted TFAP2A expression and, therefore, DDR1 transcription by recognizing the m6A modification of TFAP2A mRNA by VIRMA. In vivo animal experiments further confirmed that VIRMA and HNRNPC enhanced the TFAP2A/DDR1 axis, promoting collagen fiber alignment, reducing anti-tumor immune cell infiltration, and promoting immune escape in BC. CONCLUSION: This study demonstrated that HNRNPC promoted DDR1 transcription by recognizing VIRMA-unveiled m6A modification of TFAP2A mRNA, which enhanced collagen fiber alignment and ultimately resulted in the reduction of anti-tumor immune cell infiltration and promotion of immune escape in BC.

Role of m6A methyltransferase component VIRMA in multiple human cancers (Review).

N6-Methyladenosine (m6A) modification is one of the most widely distributed RNA modifications in eukaryotes. It participates in various RNA functions and plays vital roles in tissue development, stem cell formation and differentiation, heat shock response control, and circadian clock controlling, particularly during tumor development. The reversible regulation of m6A modification is affected by the so-called 'reader', 'writer' and 'eraser'. As a required component and the largest methyltransferase, vir-like m6A methyltransferase associated (VIRMA) can promote the progression of cancer and is associated with poor survival in multiple types of cancer. The present review investigated the role of VIRMA in various types of cancer. In an m6A-dependent or -independent manner, VIRMA can play an oncogenic role by regulating cancer cell proliferation, migration and invasion, metastasis, apoptosis resistance and tumor growth in different pathways by targeting stem factors, CCAT1/2, ID2, GATA3, CDK1, c-Jun, etc. VIRMA can also predict better prognosis in kidney renal clear cell carcinoma (KIRC), kidney renal papillary cell carcinoma (KIRP) and papillary thyroid carcinoma by TCGA analysis. The obvious oncogenic roles of VIRMA observed in different types of cancer and the mechanisms of VIRMA promoting cancers provided the basis for potential therapeutic targeting for cancer treatments.

m6A RNA modification and its writer/reader VIRMA/YTHDF3 in testicular germ cell tumors: a role in seminoma phenotype maintenance.

BACKGROUND: Covalent RNA modifications, such as N-6-methyladenosine (m6A), have been associated with various biological processes, but their role in cancer remains largely unexplored. m6A dynamics depends on specific enzymes whose deregulation may also impact in tumorigenesis. Herein, we assessed the differential abundance of m6A, its writer VIRMA and its reader YTHDF3, in testicular germ cell tumors (TGCTs), looking for clinicopathological correlates. METHODS: In silico analysis of TCGA data disclosed altered expression of VIRMA (52%) and YTHDF3 (48%), prompting subsequent validation. Formalin-fixed paraffin-embedded tissues from 122 TGCTs (2005-2016) were selected. RNA extraction, cDNA synthesis and real-time qPCR (Taqman assays) for VIRMA and YTHDF3 were performed, as well as immunohistochemistry for VIRMA, YTHDF3 and m6A, for staining intensity assessment. Associations between categorical variables were assessed using Chi square and Fisher's exact test. Distribution of continuous variables between groups was compared using the nonparametric Mann-Whitney and Kruskal-Wallis tests. Biomarker performance was assessed through receiver operating characteristics (ROC) curve construction and a cut-off was established by Youden's index method. Statistical significance was set at p < 0.05. RESULTS: In our cohort, VIRMA and YTHDF3 mRNA expression levels differed among TGCT subtypes, with Seminomas (SEs) depicting higher levels than Non-Seminomatous tumors (NSTs) (p < 0.01 for both). A positive correlation was found between VIRMA and YTHDF3 expression levels. VIRMA discriminated SEs from NSTs with AUC = 0.85 (Sensitivity 77.3%, Specificity 81.1%, PPV 71.6%, NPV 85.3%, Accuracy 79.7%). Immunohistochemistry paralleled transcript findings, as patients with strong m6A immunostaining intensity depicted significantly higher VIRMA mRNA expression levels and stronger VIRMA immunoexpression intensity (p < 0.001 and p < 0.01, respectively). CONCLUSION: Abundance of m6A and expression of VIRMA/YTHDF3 were different among TGCT subtypes, with higher levels in SEs, suggesting a contribution to SE phenotype maintenance. VIRMA and YTHDF3 might cooperate in m6A establishment in TGCTs, and their transcript levels accurately discriminate between SEs and NSTs, constituting novel candidate biomarkers for patient management.

Glycolysis and tumor progression promoted by the m6A writer VIRMA via m6A-dependent upregulation of STRA6 in pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is one of the most aggressive and lethal malignancies, highlighting the urgent need to elucidate the underlying oncogenic mechanisms. VIRMA is a classic isoform of methyltransferases that participates in epigenetic transcriptomic modification in eukaryotic mRNAs. However, the exact roles of VIRMA in PDAC remain unclear. Here, we identified that VIRMA is highly expressed in PDAC, and histone modifications of the promoter may partly account for this dysregulation. Moreover, VIRMA is closely related to glycolysis and poor prognosis in PDAC. We further determined that STRA6 is a direct downstream target of VIRMA in PDAC by RNA sequencing (RNA-seq) and m6A sequencing (m6A-seq). VIRMA is involved in gene expression regulation via 3' UTR targeting of STRA6 mRNA. Furthermore, the m6A reader IGF2BP2 was shown to critically contribute to the stability of STRA6 mRNA. We describe the role of VIRMA in promoting signaling via the STRA6/STAT3 axis, which results in increased levels of HIF-1α, a key activator of glycolysis. In vivo and in vitro experiments reveal that the VIRMA-STRA6-STAT3-HIF-1α axis plays an instrumental role in glycolysis and tumor progression in PDAC. In conclusion, we demonstrate that VIRMA can increase glycolysis in PDAC by upregulating STRA6, a cell surface membrane protein that stimulates the STAT3 pathway, thereby activating HIF-1α and leading to pancreatic cancer malignancy. Overall, our data strongly suggest that the VIRMA-STRA6-STAT3-HIF-1α axis is a viable therapeutic target in PDAC.

VIRMA promotes neuron apoptosis via inducing m6A methylation of STK10 in spinal cord injury animal models.

BACKGROUND: Spinal cord injury (SCI) occurs as a devastating neuropathic disease. The role of serine-threonine kinase 10 (STK10) in the development of SCI remains unclear. OBJECTIVE: This study aimed to investigate the action of m6A methylation on STK10 in the apoptosis of spinal cord neurons in the pathogenesis of SCI and the possible underlying mechanisms. METHODS: Rat model of SCI was established and subsequently evaluated for motor function, pathological conditions, and apoptosis of spinal cord neurons. And the effects of overexpression of STK10 on neuronal cells in animal models of spinal cord injury and glyoxylate deprivation (OGD) cell models were evaluated. m6A2Target database and SRAMP database were used to predict the m6A methylation sites of STK10. The methylation kits were used to detect overall m6A methylation. Finally, the interaction between STK10 and vir like m6A methyltransferase associated (VIRMA) was explored in animal and cellular models. RESULTS: STK10 is markedly decreased in spinal cord injury models and overexpression of STK10 inhibits neuronal apoptosis. VIRMA can induce m6A methylation of STK10. VIRMA is over-expressed in spinal cord injury models and negatively regulates the expression of STK10. m6A methylation and apoptosis of neuronal cells are reduced by the knockdown of VIRMA and STK10 shRNA have shown the opposite effects. CONCLUSIONS: VIRMA promotes neuronal apoptosis in spinal cord injury by regulating STK10 m6A methylation.

VIRMA Facilitates Triple-Negative Breast Cancer Progression via Increasing m6A-Dependent KIF15 Expression.

BACKGROUND: Vir like N6-methyladenosine (m6A) methyltransferase associated protein (VIRMA) is associated with various tumors, but the specific role of VIRMA in triple-negative breast cancer (TNBC) and the mechanisms are still unclear. Thus, in this study, in addition to the effect of VIRMA on TNBC, the underlying mechanisms were also explored. METHODS: In vitro, VIRMA expression was detected by quantitative real-time polymerase chain reaction (qRT-PCR) and western blot; VIRMA lentiviral overexpression vector (LV-VIRMA) and lentiviral vector connected with the shRNA targeting VIRMA (LV-shVIRMA) were constructed to explore the functional role of VIRMA; RNA immunoprecipitation and qRT-PCR were performed to assess the relationship between VIRMA and kinesin family 15 (KIF15). In vivo, female Balb/C mice (n = 6) were subcutaneously injected with TNBC cells transfected with LV-shRNA + LV-NC (negative control), LV-shVIRMA + LV-NC, and LV-shVIRMA + LV-KIF15, tumor volume, weight and immunohistochemistry staining of Ki-67 were employed to assess breast tumor growth; immunohistochemistry of VIRMA and KIF15 were performed to examine VIRMA and KIF15 expression in breast tumor tissues. RESULTS: Compared to normal breast epithelial cells, VIRMA was increased in TNBC cells (p < 0.01 and p < 0.001). LV-VIRMA elevated proliferation, metastasis and invasion of TNBC cells in comparison with LV-NC (p < 0.001), while VIRMA knockdown resulted in the opposite effects in comparison with LV-shRNA NC (p < 0.01 and p < 0.001). Also, compared to LV-shRNA NC, LV-shVIRMA downregulated KIF15 expression by reducing KIF15 mRNA stability (p < 0.05 and p < 0.001), which was dependent on m6A. Furthermore, compared to LV-shVIRMA + LV-NC, LV-shVIRMA + LV-KIF15 not only reversed the reduced proliferation, metastasis and invasion of TNBC cells (p < 0.05, p < 0.01, and p < 0.001), but also reversed the decreased tumor weight and volume (p < 0.05, p < 0.01, and p < 0.001). CONCLUSIONS: The above results indicated that VIRMA promoted TNBC progression by upregulating m6A-dependent KIF15 expression, providing a better understanding of the pathogenesis of TNBC.

VIRMA facilitates intrahepatic cholangiocarcinoma progression through epigenetic augmentation of TMED2 and PARD3B mRNA stabilization.

BACKGROUND: N6-methyladenine modification of RNA, a critical component of the regulatory role at the post-transcriptional level, has a crucial effect on tumor development and progression. vir-Like m6A methyltransferase associated (VIRMA) has been recently discovered as an N6-methyladenine methyltransferase; however, its specific role in intrahepatic cholangiocarcinoma (ICC) remains to be investigated in-depth. METHODS: VIRMA expression and its association with clinicopathological characteristics were evaluated using The Cancer Genome Atlas (TCGA) dataset and tissue microarrays. In vivo and in vitro assays were performed to determine the role of VIRMA in ICC proliferation and metastasis. The underlying mechanism by which VIRMA influences ICC was clarified by RNA sequencing (RNA-seq), methylated RNA immunoprecipitation sequencing (MeRIP-seq), SLAM sequencing (SLAM-seq), RNA immunoprecipitation, a luciferase reporter assay, and chromatin immunoprecipitation assay. RESULTS: VIRMA showed high expression in ICC tissues, and this finding predicted a dismal prognostic outcome. The high expression of VIRMA in ICC was due to the demethylation of H3K27me3 modification in the promoter region. Functionally, VIRMA is required for the endothelial-mesenchymal transition (EMT) process in ICC cells, as shown by multiple ICC models in in vitro and in vivo experiments. Mechanistically, multi-omics analysis using ICC cells demonstrated that TMED2 and PARD3B were the direct downstream target of VIRMA. The methylated TMED2 and PARD3B transcripts were directly recognized by HuR, which exerted stabilizing effects on its bound RNA. VIRMA-induced expression of TMED2 and PARD3B activated the Akt/GSK/ß-catenin and MEK/ERK/Slug signaling pathways, thereby promoting ICC proliferation and metastasis. CONCLUSIONS: The present study showed that VIRMA plays a critical role in ICC development by stabilizing TMED2 and PARD3B expression through the m6A-HuR-mediated mechanism. Thus, demonstrating VIRMA and its pathway as candidate therapeutic targets for ICC treatment.

The yeast RNA methylation complex consists of conserved yet reconfigured components with m6A-dependent and independent roles.

N6-methyladenosine (m6A), the most abundant mRNA modification, is deposited in mammals/insects/plants by m6A methyltransferase complexes (MTC) comprising a catalytic subunit and at least five additional proteins. The yeast MTC is critical for meiosis and was known to comprise three proteins, of which two were conserved. We uncover three novel MTC components (Kar4/Ygl036w-Vir1/Dyn2). All MTC subunits, except for Dyn2, are essential for m6A deposition and have corresponding mammalian MTC orthologues. Unlike the mammalian bipartite MTC, the yeast MTC is unipartite, yet multifunctional. The mRNA interacting module, comprising Ime4, Mum2, Vir1, and Kar4, exerts the MTC's m6A-independent function, while Slz1 enables the MTC catalytic function in m6A deposition. Both functions are critical for meiotic progression. Kar4 also has a mechanistically separate role from the MTC during mating. The yeast MTC constituents play distinguishable m6A-dependent, MTC-dependent, and MTC-independent functions, highlighting their complexity and paving the path towards dissecting multi-layered MTC functions in mammals.

IGF2BP3 is an essential N6-methyladenosine biotarget for suppressing ferroptosis in lung adenocarcinoma cells.

A lack of promising targets leads to poor prognosis in patients with lung adenocarcinoma (LUAD). Therefore, it is urgent to identify novel therapeutic targets. The importance of the N6-methyladenosine (m6A) RNA modification has been demonstrated in various types of tumors; however, knowledge of m6A-related proteins in LUAD is still limited. Here, we found that insulin-like growth factor 2 mRNA binding protein 3 (IGF2BP3), an m6A reader protein, is highly expressed in LUAD and associated with poor prognosis. IGF2BP3 desensitizes ferroptosis (a new form of regulated cell death) in a manner dependent on its m6A reading domain and binding capacity to m6A-methylated mRNAs encoding anti-ferroptotic factors, including but not limited to glutathione peroxidase 4 (GPX4), solute carrier family 3 member 2 (SLC3A2), acyl-CoA synthetase long chain family member 3 (ACSL3), and ferritin heavy chain 1 (FTH1). After IGF2BP3 overexpression, expression levels and mRNA stabilities of these anti-ferroptotic factors were successfully sustained. Notably, significant correlations between SLC3A2, ACSL3, and IGF2BP3 were revealed in clinical LUAD specimens, further establishing the essential role of IGF2BP3 in desensitizing ferroptosis. Inducing ferroptosis has been gradually accepted as an alternative strategy to treat tumors. Thus, IGF2BP3 could be a potential target for the future development of new biomaterial-associated therapeutic anti-tumor drugs.

The component of the m6A writer complex VIRMA is implicated in aggressive tumor phenotype, DNA damage response and cisplatin resistance in germ cell tumors.

BACKGROUND: Germ cell tumors (GCTs) are developmental cancers, tightly linked to embryogenesis and germ cell development. The recent and expanding field of RNA modifications is being increasingly implicated in such molecular events, as well as in tumor progression and resistance to therapy, but still rarely explored in GCTs. In this work, and as a follow-up of our recent study on this topic in TGCT tissue samples, we aim to investigate the role of N6-methyladenosine (m6A), the most abundant of such modifications in mRNA, in in vitro and in vivo models representative of such tumors. METHODS: Four cell lines representative of GCTs (three testicular and one mediastinal), including an isogenic cisplatin resistant subline, were used. CRISPR/Cas9-mediated knockdown of VIRMA was established and the chorioallantoic membrane assay was used to study its phenotypic effect in vivo. RESULTS: We demonstrated the differential expression of the various m6A writers, readers and erasers in GCT cell lines representative of the major classes of these tumors, seminomas and non-seminomas, and we evidenced changes occurring upon differentiation with all-trans retinoic acid treatment. We showed differential expression also among cells sensitive and resistant to cisplatin treatment, implicating these players in acquisition of cisplatin resistant phenotype. Knockdown of VIRMA led to disruption of the remaining methyltransferase complex and decrease in m6A abundance, as well as overall reduced tumor aggressiveness (with decreased cell viability, tumor cell proliferation, migration, and invasion) and increased sensitivity to cisplatin treatment, both in vitro and confirmed in vivo. Enhanced response to cisplatin after VIRMA knockdown was related to significant increase in DNA damage (with higher γH2AX and GADD45B levels) and downregulation of XLF and MRE11. CONCLUSIONS: VIRMA has an oncogenic role in GCTs confirming our previous tissue-based study and is further involved in response to cisplatin by interfering with DNA repair. These data contribute to our better understanding of the emergence of cisplatin resistance in GCTs and support recent attempts to therapeutically target elements of the m6A writer complex.

VIRMA-Dependent N6-Methyladenosine Modifications Regulate the Expression of Long Non-Coding RNAs CCAT1 and CCAT2 in Prostate Cancer.

RNA methylation at position N6 in adenosine (m6A) and its associated methyltransferase complex (MTC) are involved in tumorigenesis. We aimed to explore m6A biological function for long non-coding RNAs (lncRNAs) in prostate cancer (PCa) and its clinical significance. m6A and MTC levels in PCa cells were characterized by ELISA and western blot. Putative m6A-regulated lncRNAs were identified and validated by lncRNA profiler qPCR array and bioinformatics analysis, followed by m6A/RNA co-immunoprecipitation. Impact of m6A depletion on RNA stability was assessed by Actinomycin D assay. The association of m6A-levels with PCa prognosis was examined in clinical samples. Higher m6A-levels and VIRMA overexpression were detected in metastatic castration-resistant PCa (mCRPC) cells (p < 0.05). VIRMA knockdown in PC-3 cells significantly decreased m6A-levels (p = 0.0317), attenuated malignant phenotype and suppressed the expression of oncogenic lncRNAs CCAT1 and CCAT2 (p < 0.00001). VIRMA depletion and m6A reduction decreased the stability and abundance of CCAT1/2 transcripts. Higher expression of VIRMA, CCAT1, and CCAT2 as a group variable was an independent predictor of poor prognosis (HR = 9.083, CI95% 1.911-43.183, p = 0.006). VIRMA is a critical factor sustaining m6A-levels in PCa cells. VIRMA downregulation attenuates the aggressive phenotype of PCa by overall reduction of m6A-levels decreasing stability and abundance of oncogenic lncRNAs.