Conserved A-to-I RNA editing with non-conserved recoding expands the candidates of functional editing sites.

Adenosine-to-inosine (A-to-I) RNA editing recodes the genome and confers flexibility for the organisms to adapt to the environment. It is believed that RNA recoding sites are well suited for facilitating adaptive evolution by increasing the proteomic diversity in a temporal-spatial manner. The function and essentiality of a few conserved recoding sites are recognized. However, the experimentally discovered functional sites only make up a small corner of the total sites, and there is still the need to expand the repertoire of such functional sites with bioinformatic approaches. In this study, we define a new category of RNA editing sites termed 'conserved editing with non-conserved recoding' and systematically identify such sites in Drosophila editomes, figuring out their selection pressure and signals of adaptation at inter-species and intra-species levels. Surprisingly, conserved editing sites with non-conserved recoding are not suppressed and are even slightly overrepresented in Drosophila. DNA mutations leading to such cases are also favoured during evolution, suggesting that the function of those recoding events in different species might be diverged, specialized, and maintained. Finally, structural prediction suggests that such recoding in potassium channel Shab might increase ion permeability and compensate the effect of low temperature. In conclusion, conserved editing with non-conserved recoding might be functional as well. Our study provides novel aspects in considering the adaptive evolution of RNA editing sites and meanwhile expands the candidates of functional recoding sites for future validation.

BiLSTM- and CNN-Based m6A Modification Prediction Model for circRNAs.

m6A methylation, a ubiquitous modification on circRNAs, exerts a profound influence on RNA function, intracellular behavior, and diverse biological processes, including disease development. While prediction algorithms exist for mRNA m6A modifications, a critical gap remains in the prediction of circRNA m6A modifications. Therefore, accurate identification and prediction of m6A sites are imperative for understanding RNA function and regulation. This study presents a novel hybrid model combining a convolutional neural network (CNN) and a bidirectional long short-term memory network (BiLSTM) for precise m6A methylation site prediction in circular RNAs (circRNAs) based on data from HEK293 cells. This model exploits the synergy between CNN's ability to extract intricate sequence features and BiLSTM's strength in capturing long-range dependencies. Furthermore, the integrated attention mechanism empowers the model to pinpoint critical biological information for studying circRNA m6A methylation. Our model, exhibiting over 78% prediction accuracy on independent datasets, offers not only a valuable tool for scientific research but also a strong foundation for future biomedical applications. This work not only furthers our understanding of gene expression regulation but also opens new avenues for the exploration of circRNA methylation in biological research.

Non­m6A RNA modifications in haematological malignancies.

Dysregulated RNA modifications, stemming from the aberrant expression and/or malfunction of RNA modification regulators operating through various pathways, play pivotal roles in driving the progression of haematological malignancies. Among RNA modifications, N6-methyladenosine (m6A) RNA modification, the most abundant internal mRNA modification, stands out as the most extensively studied modification. This prominence underscores the crucial role of the layer of epitranscriptomic regulation in controlling haematopoietic cell fate and therefore the development of haematological malignancies. Additionally, other RNA modifications (non-m6A RNA modifications) have gained increasing attention for their essential roles in haematological malignancies. Although the roles of the m6A modification machinery in haematopoietic malignancies have been well reviewed thus far, such reviews are lacking for non-m6A RNA modifications. In this review, we mainly focus on the roles and implications of non-m6A RNA modifications, including N4-acetylcytidine, pseudouridylation, 5-methylcytosine, adenosine to inosine editing, 2'-O-methylation, N1-methyladenosine and N7-methylguanosine in haematopoietic malignancies. We summarise the regulatory enzymes and cellular functions of non-m6A RNA modifications, followed by the discussions of the recent studies on the biological roles and underlying mechanisms of non-m6A RNA modifications in haematological malignancies. We also highlight the potential of therapeutically targeting dysregulated non-m6A modifiers in blood cancer.

MeRIP-Seq for Identifying Stress-Responsive Transcriptome-Wide m6A Profiles in Plants.

Recent advancements in detection and mapping methods have enabled researchers to uncover the biological importance of RNA chemical modifications, which play a vital role in post-transcriptional gene regulation. Although numerous types of RNA modifications have been identified in higher eukaryotes, only a few have been extensively studied for their biological functions. Of these, N6-methyladenosine (m6A) is the most prevalent and important mRNA modification that influences various aspects of RNA metabolism, including mRNA stability, degradation, splicing, alternative polyadenylation, export, and localization, as well as translation. Thus, they have implications for a variety of biological processes, including growth, development, and stress responses. The m6A deposition or removal on transcripts is dynamic and is altered in response to internal and external cues. Because this mark can alter gene expression under stress conditions, it is essential to identify the transcripts that can acquire or lose this epitranscriptomic mark upon exposure to stress conditions. Here we describe a step-by-step protocol for identifying stress-responsive transcriptome-wide m6A changes using RNA immunoprecipitation followed by high-throughput sequencing (MeRIP-seq).

Reovirus infection induces transcriptome-wide unique A-to-I editing changes in the murine fibroblasts.

The conversion of Adenosine (A) to Inosine (I), by Adenosine Deaminases Acting on RNA or ADARs, is an essential post-transcriptional modification that contributes to proteome diversity and regulation in metazoans including humans. In addition to its transcriptome-regulating role, ADARs also play a major part in immune response to viral infection, where an interferon response activates interferon-stimulated genes, such as ADARp150, in turn dynamically regulating host-virus interactions. A previous report has shown that infection from reoviruses, despite strong activation of ADARp150, does not influence the editing of some of the major known editing targets, while likely editing others, suggesting a potentially nuanced editing pattern that may depend on different factors. However, the results were based on a handful of selected editing sites and did not cover the entire transcriptome. Thus, to determine whether and how reovirus infection specifically affects host ADAR editing patterns, we analyzed a publicly available deep-sequenced RNA-seq dataset, from murine fibroblasts infected with wild-type and mutant reovirus strains that allowed us to examine changes in editing patterns on a transcriptome-wide scale. To the best of our knowledge, this is the first transcriptome-wide report on host editing changes after reovirus infection. Our results demonstrate that reovirus infection induces unique nuanced editing changes in the host, including introducing sites uniquely edited in infected samples. Genes with edited sites are overrepresented in pathways related to immune regulation, cellular signaling, metabolism, and growth. Moreover, a shift in editing targets has also been observed, where the same genes are edited in infection and control conditions but at different sites, or where the editing rate is increased for some and decreased for other differential targets, supporting the hypothesis of dynamic and condition-specific editing by ADARs.

Functions and mechanisms of A-to-I RNA editing in filamentous ascomycetes.

Although lack of ADAR (adenosine deaminase acting on RNA) orthologs, genome-wide A-to-I editing occurs specifically during sexual reproduction in a number of filamentous ascomycetes, including Fusarium graminearum and Neurospora crassa. Unlike ADAR-mediated editing in animals, fungal A-to-I editing has a strong preference for hairpin loops and U at -1 position, which leads to frequent editing of UAG and UAA stop codons. Majority of RNA editing events in fungi are in the coding region and cause amino acid changes. Some of these editing events have been experimentally characterized for providing heterozygote and adaptive advantages in F. graminearum. Recent studies showed that FgTad2 and FgTad3, 2 ADAT (adenosine deaminase acting on tRNA) enzymes that normally catalyze the editing of A34 in the anticodon of tRNA during vegetative growth mediate A-to-I mRNA editing during sexual reproduction. Stage specificity of RNA editing is conferred by stage-specific expression of short transcript isoforms of FgTAD2 and FgTAD3 as well as cofactors such as AME1 and FIP5 that facilitate the editing of mRNA in perithecia. Taken together, fungal A-to-I RNA editing during sexual reproduction is catalyzed by ADATs and it has the same sequence and structural preferences with editing of A34 in tRNA.

METTL3 promotes the progression of osteosarcoma through the N6-methyladenosine modification of MCAM via IGF2BP1.

BACKGROUND: The molecular mechanisms of osteosarcoma (OS) are complex. In this study, we focused on the functions of melanoma cell adhesion molecule (MCAM), methyltransferase 3 (METTL3) and insulin like growth factor 2 mRNA binding protein 1 (IGF2BP1) in OS development. METHODS: qRT-PCR assay and western blot assay were performed to determine mRNA and protein expression of MCAM, METTL3, IGF2BP1 and YY1. MTT assay and colony formation assay were conducted to assess cell proliferation. Cell apoptosis, invasion and migration were evaluated by flow cytometry analysis, transwell assay and wound-healing assay, respectively. Methylated RNA Immunoprecipitation (MeRIP), dual-luciferase reporter, Co-IP, RIP and ChIP assays were performed to analyze the relationships of MCAM, METTL3, IGF2BP1 and YY1. The functions of METTL3 and MCAM in tumor growth were explored through in vivo experiments. RESULTS: MCAM was upregulated in OS, and MCAM overexpression promoted OS cell growth, invasion and migration and inhibited apoptosis. METTL3 and IGF2BP1 were demonstrated to mediate the m6A methylation of MCAM. Functionally, METTL3 or IGF2BP1 silencing inhibited OS cell progression, while MCAM overexpression ameliorated the effects. Transcription factor YY1 promoted the transcription level of METTL3 and regulated METTL3 expression in OS cells. Additionally, METTL3 deficiency suppressed tumor growth in vivo, while MCAM overexpression abated the effect. CONCLUSION: YY1/METTL3/IGF2BP1/MCAM axis aggravated OS development, which might provide novel therapy targets for OS.

An association study of m6A methylation with major depressive disorder.

OBJECTIVE: To find the relationship between N6-methyladenosine (m6A) genes and Major Depressive Disorder (MDD). METHODS: Differential expression of m6A associated genes between normal and MDD samples was initially identified. Subsequent analysis was conducted on the functions of these genes and the pathways they may affect. A diagnostic model was constructed using the expression matrix of these differential genes, and visualized using a nomogram. Simultaneously, an unsupervised classification method was employed to classify all patients based on the expression of these m6A associated genes. Following this, common differential genes among different clusters were computed. By analyzing the functions of the common differential expressed genes among clusters, the role of m6A-related genes in the pathogenesis of MDD patients was elucidated. RESULTS: Differential expression was observed in ELAVL1 and YTHDC2 between the MDD group and the control group. ELAVL1 was associated with comorbid anxiety in MDD patients. A linear regression model based on these two genes could accurately predict whether patients in the GSE98793 dataset had MDD and could provide a net benefit for clinical decision-making. Based on the expression matrix of ELAVL1 and YTHDC2, MDD patients were classified into three clusters. Among these clusters, there were 937 common differential genes. Enrichment analysis was also performed on these genes. The ssGSEA method was applied to predict the content of 23 immune cells in the GSE98793 dataset samples. The relationship between these immune cells and ELAVL1, YTHDC2, and different clusters was analyzed. CONCLUSION: Among all the m6A genes, ELAVL1 and YTHDC2 are closely associated with MDD, ELAVL1 is related to comorbid anxiety in MDD. ELAVL1 and YTHDC2 have opposite associations with immune cells in MDD.

Direct Analysis of HIV mRNA m6A Methylation by Nanopore Sequencing.

The post-transcriptional processing and chemical modification of HIV RNA are understudied aspects of HIV virology, primarily due to the limited ability to accurately map and quantify RNA modifications. Modification-specific antibodies or modification-sensitive endonucleases coupled with short-read RNA sequencing technologies have allowed for low-resolution or limited mapping of important regulatory modifications of HIV RNA such as N6-methyladenosine (m6A). However, a high-resolution map of where these sites occur on HIV transcripts is needed for detailed mechanistic understanding. This has recently become possible with new sequencing technologies. Here, we describe the direct RNA sequencing of HIV transcripts using an Oxford Nanopore Technologies sequencer and the use of this technique to map m6A at near single nucleotide resolution. This technology also provides the ability to identify splice variants with long RNA reads and thus, can provide high-resolution RNA modification maps that distinguish between overlapping splice variants. The protocols outlined here for m6A also provide a powerful paradigm for studying any other RNA modifications that can be detected on the nanopore platform.

Role of N6-methyladenosine-related lncRnas in pseudoexfoliation glaucoma.

To explore the role of lncRNA m6A methylation modification in aqueous humour (AH) of patients with pseudoexfoliation glaucoma (PXG). Patients with open-angle PXG under surgery from June 2021 to December 2021 were selected. Age- and gender-matched patients with age-related cataract (ARC) were chosen as control. Patients underwent detailed ophthalmic examinations. 0.05-0.1 ml AH were extracted during surgery for MeRIP-Seq and RNA-Seq. Joint analysis was used to screen lncRNAs with differential m6A methylation modification and expression. Online software tools were used to draw lncRNA-miRNA-mRNA network (ceRNA). Expression of lncRNAs and mRNAs was confirmed using quantitative real-time PCR. A total of 4151 lncRNAs and 4386 associated m6A methylation modified peaks were identified in the PXG group. Similarly, 2490 lncRNAs and 2595 associated m6A methylation modified peaks were detected in the control. Compared to the ARC group, the PXG group had 234 hypermethylated and 402 hypomethylated m6A peaks, with statistically significant differences (| Fold Change (FC) |≥2, p < 0.05). Bioinformatic analysis revealed that these differentially methylated lncRNA enriched in extracellular matrix formation, tight adhesion, TGF- ß signalling pathway, AMPK signalling pathway, and MAPK signalling pathway. Joint analysis identified 10 lncRNAs with differential m6A methylation and expression simultaneously. Among them, the expression of ENST000000485383 and ROCK1 were confirmed downregulated in the PXG group by RT-qPCR. m6A methylation modification may affect the expression of lncRNA and participate in the pathogenesis of PXG through the ceRNA network. ENST000000485383-hsa miR592-ROCK1 May be a potential target pathway for further investigation in PXG m6A methylation.

The metabolic baton: conducting the dance of N6-methyladenosine writing and erasing.

The modification N6-methyladenosine (m6A) plays an important role in determining the functional output of gene expression programs. Throughout the transcriptome, the levels of m6A are tightly regulated by the opposing activities of methyltransferases and demethylases, as well as the interaction of modified transcripts with m6A-dependent RNA-binding proteins that modulate transcript stability, often referred to as writers, erasers, and readers. The enzymatic activities of both writers and erasers are tightly linked to the cellular metabolic environment, as these enzymatic reactions rely on metabolism intermediaries as cofactors. In this review, we highlight the examples of intersection between metabolism and m6A-dependent gene regulation and discuss the different contexts where this interaction plays important roles.

Characterization of the m6A gene family in switchgrass and functional analysis of PvALKBH10 during flowering.

N6-methyladenosine (m6A), a nucleotide modification that is frequently seen in RNA, plays a crucial role in plant growth, development and stress resistance. However, the m6A regulatory machinery in switchgrass (Panicum virgatum L.), a model plant for cellulose-to-ethanol conversion, remains largely unknown. In this study, we identified 57 candidate genes involved in m6A-regulation in the switchgrass genome, and analyzed their chromosomal distribution, evolutionary relationships, and functions. Notably, we observed distinct gene expression patterns under salt and drought stress, with salt stress inducing writer and eraser genes, alongside drought stress predominantly affecting reader genes. Additionally, we knocked out PvALKBH10, an m6A demethylase gene, via CRISPR/Cas9 and found its potential function in controlling flowering time. This study provides insight into the genomic organization and evolutionary features of m6A-associated putative genes in switchgrass, and therefore serves as the basis for further functional studies.

Quantitative profiling N1-methyladenosine (m1A) RNA methylation from Oxford nanopore direct RNA sequencing data.

With the recent advanced direct RNA sequencing technique that proposed by the Oxford Nanopore Technologies, RNA modifications can be detected and profiled in a simple and straightforward manner. Majority nanopore-based modification studies were devoted to those popular types such as m6A and pseudouridine. To address current limitations on studying the crucial regulator, m1A modification, we conceived this study. We have developed an integrated computational workflow designed for the detection of m1A modifications from direct RNA sequencing data. This workflow comprises a feature extractor responsible for capturing signal characteristics (such as mean, standard deviations, and length of electric signals), a single molecule-level m1A predictor trained with features extracted from the IVT dataset using classical machine learning algorithms, a confident m1A site selector employing the binomial test to identify statistically significant m1A sites, and an m1A modification rate estimator. Our model achieved accurate molecule-level prediction (Average AUC = 0.9689) and reliable m1A site detection and quantification. To show the feasibility of our workflow, we conducted a study on in vivo transcribed human HEK293 cell line, and the results were carefully annotated and compared with other techniques (i.e., Illumina sequencing-based techniques). We believed that this tool will enabling a comprehensive understanding of the m1A modification and its functional mechanisms within cells and organisms.

Crosstalk between histone/DNA modifications and RNA N6-methyladenosine modification.

N6-methyladenosine (m6A) is the most prevalent internal RNA modification in eukaryotic messenger RNAs (mRNAs), regulating gene expression at the transcription and post-transcription levels. Complex interplay between m6A and other well-studied epigenetic modifications, including histone modifications and DNA modification, has been extensively reported in recent years. The crosstalk between RNA m6A modification and histone/DNA modifications plays a critical role in establishing the chromatin state for the precise and specific fine-tuning of gene expression and undoubtedly has profound impacts on both physiological and pathological processes. In this review, we discuss the crosstalk between RNA m6A modification and histone/DNA modifications, emphasizing their sophisticated communications and the mechanisms underlying to gain a comprehensive view of the biological relevance of m6A-based epigenetic network.

METTL3-mediated m6A methylation of SLC38A1 stimulates cervical cancer growth.

The objective of this study was to better characterize the role of the glutamine transporter SLC38A1 in cervical cancer and explore the underlying mechanisms. Data from public databases and clinical cervical cancer tissue samples were used to assess the expression of SLC38A1 and its prognostic significance. Immunohistochemical staining, qRT-PCR, and Western blotting were used to evaluate the expression of relevant genes and proteins. Cell viability, cell cycle, apoptosis, and intracellular glutamine content were measured using CCK-8, flow cytometry, and biochemical assays. Additionally, the RNA immunoprecipitation (RIP) assay was used to examine the impact of METTL3/IGF2BP3 on the m6A modification of the SLC38A1 3'UTR. Both cervical cancer specimens and cells showed significantly increased expression of SLC38A1 and its expression correlated with an unfavorable prognosis. Knockdown of SLC38A1 inhibited cell viability and cell cycle progression, induced apoptosis, and suppressed tumor growth in vivo. Glutaminase-1 inhibitor CB-839 reversed the effects of SLC38A1 overexpression. METTL3 promoted m6A modification of SLC38A1 and enhanced its mRNA stability through IGF2BP3 recruitment. Moreover, METTL3 silencing inhibited cell viability, cell cycle progression, intracellular glutamine content, and induced apoptosis, but these effects were reversed by SLC38A1 overexpression. In conclusion, METTL3-mediated m6A methylation of SLC38A1 stimulates cervical cancer progression. SLC38A1 inhibition is a potential therapeutic strategy for cervical cancer.

The N6-methyladenosine methylation landscape stratifies breast cancer into two subtypes with distinct immunological characteristics.

N6-methyladenosine (m6A) methylation modification affects the tumorigenesis and metastasis of breast cancer (BC). This study investigated the association between m6A regulator-mediated methylation modification patterns and characterization of the tumour microenvironment in BC, as well as their prognostic importance. Public gene expression data and clinical annotations were collected from The Cancer Genome Atlas (TCGA) database, the Gene Expression Omnibus website and the METABRIC program. We analysed the genetic expression, gene-gene interactions, gene mutations and copy number variations using R software. The data were screened for risk genes using the Cox risk regression model, and we developed an algorithm for risk score and its predictive value. Compared to adjacent normal tissue, we identified 16 differentially expressed m6A regulators in BC, including six writers and 10 readers. Under unsupervised clustering, two distinguished modification patterns were identified, cluster C1 and C2. Compared to m6A cluster C2, cluster C1 was found to be more involved in immune-related pathways, with a relatively higher immune score and stromal score (P < 0.05). Patients were divided into two groups based on their risk scores for survival analysis. The patients in the high-risk score group had significantly worse overall survival than patients in the low-risk score group, (P < 0.0001). The TCGA database validation revealed the same prognostic tendency. In summary, our study showed distinct m6A regulator modification patterns contribute to the immunological heterogeneity and diversity of BC. The development of m6A gene signatures and the m6A score aid in the prognostic prediction of patients with BC.

DDX21 mediates co-transcriptional RNA m6A modification to promote transcription termination and genome stability.

N6-methyladenosine (m6A) is a crucial RNA modification that regulates diverse biological processes in human cells, but its co-transcriptional deposition and functions remain poorly understood. Here, we identified the RNA helicase DDX21 with a previously unrecognized role in directing m6A modification on nascent RNA for co-transcriptional regulation. DDX21 interacts with METTL3 for co-recruitment to chromatin through its recognition of R-loops, which can be formed co-transcriptionally as nascent transcripts hybridize onto the template DNA strand. Moreover, DDX21's helicase activity is needed for METTL3-mediated m6A deposition onto nascent RNA following recruitment. At transcription termination regions, this nexus of actions promotes XRN2-mediated termination of RNAPII transcription. Disruption of any of these steps, including the loss of DDX21, METTL3, or their enzymatic activities, leads to defective termination that can induce DNA damage. Therefore, we propose that the R-loop-DDX21-METTL3 nexus forges the missing link for co-transcriptional modification of m6A, coordinating transcription termination and genome stability.

N6-methyladenosine-modified circTEAD1 stabilizes Yap1 mRNA to promote chordoma tumorigenesis.

BACKGROUND: Chordoma, a rare bone tumour with aggressive local invasion and high recurrence rate with limited understanding of its molecular mechanisms. Circular RNAs (circRNAs) have been extensively implicated in tumorigenesis, yet their involvement in chordoma remains largely unexplored. N6-methyladenosine (m6A) modification holds a crucial function in regulating protein translation, RNA degradation and transcription. METHODS: Initially, screening and validation of circTEAD1 in chordoma were conducted by high-throughput sequencing. Subsequently, sh-circTEAD1 and an overexpression plasmid were constructed. Colony formation assays, cell counting kit-8, Transwell and wound healing assays were utilized to validate the function of circTEAD1 in vitro. RNA pull-down assays identified the binding proteins of circTEAD1, which underwent verification through RNA immunoprecipitation (RIP). Methylated RIP assays were conducted to detect the m6A binding sites. Following this, luciferase assay, RT-qPCR, RIP and Western blotting analyses were conducted, revealing that Yap1 was the direct target of circTEAD1. Afterwards, the same methods were utilized for the validation of the function of Yap1 in chordoma in vitro. Finally, the regulatory relationship between circTEAD1 and Yap1 in chordoma was verified by an in vivo tumour formation assay. RESULTS: CircTEAD1 was identified as an upregulated circRNA in chordoma specimens, with heightened circTEAD1 expression emerging as a prognostic indicator. In vitro experiments convincingly demonstrated that circTEAD1 significantly promoted chordoma cell invasion, migration and aggressiveness. Furthermore, the analysis revealed that methyltransferase-like 3-mediated m6A modification facilitated the cytoplasmic export of circTEAD1. The circTEAD1/IGF2BP3/Yap1 mRNA RNA-protein ternary complex not only bolstered the stability of Yap1 mRNA but also exerted a pivotal role in driving chordoma tumorigenesis. CONCLUSIONS: In this study, the role of m6A-modified circTEAD1 in chordoma was identified. The findings offer novel insights into the potential molecular targets for chordoma therapy, shedding light on the intricate interplay between circRNAs, m6A modification and Yap1 mRNA in chordoma pathogenesis.

Circ_0053943 complexed with IGF2BP3 drives uveal melanoma progression via regulating N6-methyladenosine modification of Epidermal growth factor receptor.

Numerous studies have characterized the critical role of circular RNAs (circRNAs) as regulatory factors in the progression of multiple cancers. However, the biological functions of circRNAs and their underlying molecular mechanisms in the progression of uveal melanoma (UM) remain enigmatic. In this study, we identified a novel circRNA, circ_0053943, through re-analysis of UM microarray data and quantitative RT-PCR. Circ_0053943 was found to be upregulated in UM and to promote the proliferation and metastatic ability of UM cells in both in vitro and in vivo settings. Mechanistically, circ_0053943 was observed to bind to the KH1 and KH2 domains of insulin-like growth factor 2 mRNA-binding protein 3 (IGF2BP3), thereby enhancing the function of IGF2BP3 by stabilizing its target mRNA. RNA sequencing assays identified epidermal growth factor receptor (EGFR) as a target gene of circ_0053943 and IGF2BP3 at the transcriptional level. Rescue assays demonstrated that circ_0053943 exerts its biological function by stabilizing EGFR mRNA and regulating the downstream mitogen-activated protein kinase/extracellular signal-regulated kinase (MAPK/ERK) signaling pathway. Collectively, circ_0053943 may promote UM progression by stabilizing EGFR mRNA and activating the MAPK/ERK signaling pathway through the formation of a circ_0053943/IGF2BP3/EGFR RNA-protein ternary complex, thus providing a potential biomarker and therapeutic target for UM.

Single-molecule epitranscriptomic analysis of full-length HIV-1 RNAs reveals functional roles of site-specific m6As.

Although the significance of chemical modifications on RNA is acknowledged, the evolutionary benefits and specific roles in human immunodeficiency virus (HIV-1) replication remain elusive. Most studies have provided only population-averaged values of modifications for fragmented RNAs at low resolution and have relied on indirect analyses of phenotypic effects by perturbing host effectors. Here we analysed chemical modifications on HIV-1 RNAs at the full-length, single RNA level and nucleotide resolution using direct RNA sequencing methods. Our data reveal an unexpectedly simple HIV-1 modification landscape, highlighting three predominant N6-methyladenosine (m6A) modifications near the 3' end. More densely installed in spliced viral messenger RNAs than in genomic RNAs, these m6As play a crucial role in maintaining normal levels of HIV-1 RNA splicing and translation. HIV-1 generates diverse RNA subspecies with distinct m6A ensembles, and maintaining multiple of these m6As on its RNAs provides additional stability and resilience to HIV-1 replication, suggesting an unexplored viral RNA-level evolutionary strategy.