From A to m6A: The Emerging Viral Epitranscriptome

Recent insights into N6-methyladenosine during viral infection

The RNA modification of N6-methyladenosine (m6A) controls many aspects of RNA function that impact biological processes, including viral infection. In this review, we highlight recent work that shapes our current understanding of the diverse mechanisms by which m6A can regulate viral infection by acting on viral or cellular mRNA molecules. We focus on emerging concepts and understanding, including how viral infection alters the localization and function of m6A machinery proteins, how m6A regulates antiviral innate immunity, and the multiple roles of m6A in regulating specific viral infections. We also summarize the recent studies on m6A during SARS-CoV-2 infection, focusing on points of convergence and divergence. Ultimately, this review provides a snapshot of the latest research on m6A during viral infection.

Elucidation of the Epitranscriptomic RNA Modification Landscape of Chikungunya Virus

The genomes of positive-sense (+) single-stranded RNA (ssRNA) viruses are believed to be subjected to a wide range of RNA modifications. In this study, we focused on the chikungunya virus (CHIKV) as a model (+) ssRNA virus to study the landscape of viral RNA modification in infected human cells. Among the 32 distinct RNA modifications analysed by mass spectrometry, inosine was found enriched in the genomic CHIKV RNA. However, orthogonal validation by Illumina RNA-seq analyses did not identify any inosine modification along the CHIKV RNA genome. Moreover, CHIKV infection did not alter the expression of ADAR1 isoforms, the enzymes that catalyse the adenosine to inosine conversion. Together, this study highlights the importance of a multidisciplinary approach to assess the presence of RNA modifications in viral RNA genomes.

Binding Evolution of the Dengue Virus Envelope Against DC-SIGN: A Combined Approach of Phylogenetics and Molecular Dynamics Analyses Over 30 Years of Dengue Virus in Brazil

The Red Queen Hypothesis (RQH), derived from Lewis Carroll's "Through the Looking-Glass", postulates that organisms must continually adapt in response to each other to maintain relative fitness. Within the context of host-pathogen interactions, the RQH implies an evolutionary arms race, wherein viruses evolve to exploit hosts and hosts evolve to resist viral invasion. This study delves into the dynamics of the RQH in the context of virus-cell interactions, specifically focusing on virus receptors and cell receptors. We observed multiple virus-host systems and noted patterns of co-evolution. As viruses evolved receptor-binding proteins to effectively engage with cell receptors, cells countered by altering their receptor genes. This ongoing mutual adaptation cycle has influenced the molecular intricacies of receptor-ligand interactions. Our data supports the RQH as a driving force behind the diversification and specialization of both viral and host cell receptors. Understanding this co-evolutionary dance offers insights into the unpredictability of emerging viral diseases and potential therapeutic interventions. Future research is crucial to dissect the nuanced molecular changes and the broader ecological consequences of this ever-evolving battle. Here, we combine phylogenetic inferences, structural modeling, and molecular dynamics analyses to describe the epidemiological characteristics of major Brazilian DENV strains that circulated from 1990 to 2022 from a combined perspective, thus providing us with a more detailed picture on the dynamics of such interactions over time.

Ramifications of m6A Modification on ncRNAs in Cancer

N6-methyladenosine (m6A) is an RNA modification wherein the N6-position of adenosine is methylated. It is one of the most prevalent internal modifications of RNA and regulates various aspects of RNA metabolism. M6A is deposited by m6A methyltransferases, removed by m6A demethylases, and recognized by reader proteins, which modulate splicing, export, translation, and stability of the modified mRNA. Recent evidence suggests that various classes of non- coding RNAs (ncRNAs), including microRNAs (miRNAs), circular RNAs (circRNAs), and long con-coding RNAs (lncRNAs), are also targeted by this modification. Depending on the ncRNA species, m6A may affect the processing, stability, or localization of these molecules. The m6A- modified ncRNAs are implicated in a number of diseases, including cancer. In this review, the author summarizes the role of m6A modification in the regulation and functions of ncRNAs in tumor development. Moreover, the potential applications in cancer prognosis and therapeutics are discussed.

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.

Challenges to mapping and defining m6A function in viral RNA

Viral RNA molecules contain multiple layers of regulatory information. This includes features beyond the primary sequence, such as RNA structures and RNA modifications, including N6-methyladenosine (m6A). Many recent studies have identified the presence and location of m6A in viral RNA and have found diverse regulatory roles for this modification during viral infection. However, to date, viral m6A mapping strategies have limitations that prevent a complete understanding of the function of m6A on individual viral RNA molecules. While m6A sites have been profiled on bulk RNA from many viruses, the resulting m6A maps of viral RNAs described to date present a composite picture of m6A across viral RNA molecules in the infected cell. Thus, for most viruses, it is unknown if unique viral m6A profiles exist throughout infection, nor if they regulate specific viral life cycle stages. Here, we describe several challenges to defining the function of m6A in viral RNA molecules and provide a framework for future studies to help in the understanding of how m6A regulates viral infection.

N6-methyladenosine modification is not a general trait of viral RNA genomes

Despite the nuclear localization of the m6A machinery, the genomes of multiple exclusively-cytoplasmic RNA viruses, such as chikungunya (CHIKV) and dengue (DENV), are reported to be extensively m6A-modified. However, these findings are mostly based on m6A-Seq, an antibody-dependent technique with a high rate of false positives. Here, we address the presence of m6A in CHIKV and DENV RNAs. For this, we combine m6A-Seq and the antibody-independent SELECT and nanopore direct RNA sequencing techniques with functional, molecular, and mutagenesis studies. Following this comprehensive analysis, we find no evidence of m6A modification in CHIKV or DENV transcripts. Furthermore, depletion of key components of the host m6A machinery does not affect CHIKV or DENV infection. Moreover, CHIKV or DENV infection has no effect on the m6A machinery's localization. Our results challenge the prevailing notion that m6A modification is a general feature of cytoplasmic RNA viruses and underscore the importance of validating RNA modifications with orthogonal approaches.

Utilization of nanopore direct RNA sequencing to analyze viral RNA modifications

Modifications on viral RNAs (vRNAs), either genomic RNAs or RNA transcripts, have complex effects on the viral life cycle and cellular responses to viral infection. The advent of Oxford Nanopore Technologies Direct RNA Sequencing provides a new strategy for studying RNA modifications. To this end, multiple computational tools have been developed, but a systemic evaluation of their performance in mapping vRNA modifications is lacking. Here, 10 computational tools were tested using the Sindbis virus (SINV) RNAs isolated from infected mammalian (BHK-21) or mosquito (C6/36) cells, with in vitro-transcribed RNAs serving as modification-free control. Three single-mode approaches were shown to be inapplicable in the viral context, and three out of seven comparative methods required cutoff adjustments to reduce false-positive predictions. Utilizing optimized cutoffs, an integrated analysis of comparative tools suggested that the intersected predictions of Tombo_com and xPore were significantly enriched compared with the background. Consequently, a pipeline integrating Tombo_com and xPore was proposed for vRNA modification detection; the performance of which was supported by N6-methyladenosine prediction in severe acute respiratory syndrome coronavirus 2 RNAs using publicly available data. When applied to SINV RNAs, this pipeline revealed more intensive modifications in subgenomic RNAs than in genomic RNAs. Modified uridines were frequently identified, exhibiting substantive overlapping between vRNAs generated in different cell lines. On the other hand, the interpretation of other modifications remained unclear, underlining the limitations of the current computational tools despite their notable potential.IMPORTANCEComputational approaches utilizing Oxford Nanopore Technologies Direct RNA Sequencing data were almost exclusively designed to map eukaryotic epitranscriptomes. Therefore, extra caution must be exercised when using these tools to detect vRNA modifications, as in most cases, vRNA modification profiles should be regarded as unknown epitranscriptomes without prior knowledge. Here, we comprehensively evaluated the performance of 10 computational tools in detecting vRNA modification sites. All tested single-mode methods failed to differentiate native and in vitro-transcribed samples. Using optimized cutoff values, seven tested comparative tools generated very different predictions. An integrated analysis showed significant enrichment of Tombo_com and xPore predictions against the background. A pipeline for vRNA modification detection was proposed accordingly and applied to Sindbis virus RNAs. In conclusion, our study underscores the need for the careful application of computational tools to analyze viral epitranscriptomics. It also offers insights into alphaviral RNA modifications, although further validation is required.

Modifying the antiviral innate immune response by selective writing, erasing, and reading of m6A on viral and cellular RNA

Viral infection and the antiviral innate immune response are regulated by the RNA modification m6A. m6A directs nearly all aspects of RNA metabolism by recruiting RNA-binding proteins that mediate the fate of m6A-containing RNA. m6A controls the antiviral innate immune response in diverse ways, including shielding viral RNA from detection by antiviral sensors and influencing the expression of cellular mRNAs encoding antiviral signaling proteins, cytokines, and effector proteins. While m6A and the m6A machinery are important for the antiviral response, the precise mechanisms that determine how the m6A machinery selects specific viral or cellular RNA molecules for modification during infection are not fully understood. In this review, we highlight recent findings that shed light on how viral infection redirects the m6A machinery during the antiviral response. A better understanding of m6A targeting during viral infection could lead to new immunomodulatory and therapeutic strategies against viral infection.

Mapping m6A Sites on HIV-1 RNA Using Oligonucleotide LC-MS/MS

The biological significance of chemical modifications to the ribonucleic acid (RNA) of human immunodeficiency virus type-1 (HIV-1) has been recognized. However, our understanding of the site-specific and context-dependent roles of these chemical modifications remains limited, primarily due to the absence of nucleotide-resolution mapping of modification sites. In this study, we present a method for achieving nucleotide-resolution mapping of chemical modification sites on HIV-1 RNA using liquid chromatography and tandem mass spectrometry (LC-MS/MS). LC-MS/MS, a powerful tool capable of directly analyzing native RNAs, has proven effective for mapping RNA modifications in small RNA molecules, including ribosomal RNA and transfer RNA. However, longer RNAs have posed challenges, such as the 9 Kb HIV-1 virion RNA, due to the complexity of and ambiguity in mass differences among RNase T1-cleaved RNA fragments in LC-MS/MS data. Here, we introduce a new target RNA enrichment method to isolate small local RNA fragments of HIV-1 RNA that potentially harbor site-specific N6-methyladenosine (m6A) modifications. In our initial trial, we used target-specific DNA probes only and encountered insufficient RNA fragmentation due to inefficient S1 digestion near the target site. Recognizing that inefficient S1 digestion by HIV-1 RNA is likely due to the formation of secondary structures in proximity to the target site, we designed multiple DNA probes annealing to various sites of HIV-1 RNA to better control the structures of RNA substrates for S1 digestion. The use of these non-target DNA probes significantly improved the isolation of more homogeneous target RNA fragments of approximately 50 bases in length. Oligonucleotide LC-MS/MS analysis of these isolated target RNA fragments successfully separated and detected both m6A-methylated and non-methylated oligomers at the two m6A-predicted sites. The principle of this new target enrichment strategy holds promise and should be broadly applicable to the analysis of any lengthy RNA that was previously deemed infeasible for investigation using oligonucleotide LC-MS/MS.

Plant YTHDF proteins are direct effectors of antiviral immunity against an N6-methyladenosine-containing RNA virus

In virus-host interactions, nucleic acid-directed first lines of defense that allow viral clearance without compromising growth are of paramount importance. Plants use the RNA interference pathway as a basal antiviral immune system, but additional RNA-based mechanisms of defense also exist. The infectivity of a plant positive-strand RNA virus, alfalfa mosaic virus (AMV), relies on the demethylation of viral RNA by the recruitment of the cellular N6-methyladenosine (m6 A) demethylase ALKBH9B, but how demethylation of viral RNA promotes AMV infection remains unknown. Here, we show that inactivation of the Arabidopsis cytoplasmic YT521-B homology domain (YTH)-containing m6 A-binding proteins ECT2, ECT3, and ECT5 is sufficient to restore AMV infectivity in partially resistant alkbh9b mutants. We further show that the antiviral function of ECT2 is distinct from its previously demonstrated function in the promotion of primordial cell proliferation: an ect2 mutant carrying a small deletion in its intrinsically disordered region is partially compromised for antiviral defense but not for developmental functions. These results indicate that the m6 A-YTHDF axis constitutes a novel branch of basal antiviral immunity in plants.

m6A Regulates the Stability of Cellular Transcripts Required for Efficient KSHV Lytic Replication

The epitranscriptomic modification N⁶-methyladenosine (m⁶A) is a ubiquitous feature of the mammalian transcriptome. It modulates mRNA fate and dynamics to exert regulatory control over numerous cellular processes and disease pathways, including viral infection. Kaposi's sarcoma-associated herpesvirus (KSHV) reactivation from the latent phase leads to the redistribution of m⁶A topology upon both viral and cellular mRNAs within infected cells. Here we investigate the role of m⁶A in cellular transcripts upregulated during KSHV lytic replication. Our results show that m⁶A is crucial for the stability of the GPRC5A mRNA, whose expression is induced by the KSHV latent-lytic switch master regulator, the replication and transcription activator (RTA) protein. Moreover, we demonstrate that GPRC5A is essential for efficient KSHV lytic replication by directly regulating NFκB signalling. Overall, this work highlights the central importance of m⁶A in modulating cellular gene expression to influence viral infection.

Cells at early and late stages of infection with Senecavirus A: Comparative analysis of N6-methyladenosine modification on mRNAs

Infection with Senecavirus A (SVA) causes differential phenotypes in cells. In this study, cells were inoculated with SVA for culture. At 12 and 72 h post infection, cells were independently harvested for high-throughput RNA sequencing, and further methylated RNA immunoprecipitation sequencing. The resultant data were comprehensively analyzed for mapping N⁶-methyladenosine (m⁶A)-modified profiles of SVA-infected cells. More importantly, m⁶A-modified regions were identified in the SVA genome. A dataset of m⁶A-modified mRNAs was generated for screening out differentially m⁶A-modified mRNAs, further subjected to a series of in-depth analyses. This study not only showed statistical differentiation of m⁶A-modified sites between two SVA-infected groups, but also demonstrated that SVA genome, as a positive-sense, single-stranded mRNA, itself could be modified through the m⁶A pattern. Out of the six samples of SVA mRNAs, only three were identified to be m⁶A-modified, implying that the epigenetic effect might not be a crucial driving force for SVA evolution.

m6A reader proteins: the executive factors in modulating viral replication and host immune response

N6-Methyladenosine (m6A) modification is the most abundant covalent modification of RNA. It is a reversible and dynamic process induced by various cellular stresses including viral infection. Many m6A methylations have been discovered, including on the genome of RNA viruses and on RNA transcripts of DNA viruses, and these methylations play a positive or negative role on the viral life cycle depending on the viral species. The m6A machinery, including the writer, eraser, and reader proteins, achieves its gene regulatory role by functioning in an orchestrated manner. Notably, data suggest that the biological effects of m6A on target mRNAs predominantly depend on the recognition and binding of different m6A readers. These readers include, but are not limited to, the YT521-B homology (YTH) domain family, heterogeneous nuclear ribonucleoproteins (HNRNPs), insulin-like growth factor 2 mRNA-binding proteins (IGF2BPs), and many others discovered recently. Indeed, m6A readers have been recognized not only as regulators of RNA metabolism but also as participants in a variety of biological processes, although some of these reported roles are still controversial. Here, we will summarize the recent advances in the discovery, classification, and functional characterization of m6A reader proteins, particularly focusing on their roles and mechanisms of action in RNA metabolism, gene expression, and viral replication. In addition, we also briefly discuss the m6A-associated host immune responses in viral infection.

Infection phase-dependent dynamics of the viral and host N6-methyladenosine epitranscriptome in the lifecycle of an oncogenic virus in vivo

Dynamic alteration of the epitranscriptome exerts regulatory effects on the lifecycle of oncogenic viruses in vitro. However, little is known about these effects in vivo because of the general lack of suitable animal infection models of these viruses. Using a model of rapid-onset Marek's disease lymphoma in chickens, we investigated changes in viral and host messenger RNA (mRNA) N6-methyladenosine (m⁶ A) modification during Marek's disease virus (MDV) infection in vivo. We found that the expression of major epitranscriptomic proteins varies among viral infection phases, reprogramming both the viral and the host epitranscriptomes. Specifically, the methyltransferase-like 3 (METTL3)/14 complex was suppressed during the lytic and reactivation phases of the MDV lifecycle, whereas its expression was increased during the latent phase and in MDV-induced tumors. METTL3/14 overexpression inhibits, whereas METTL3/14 knockdown enhances, MDV gene expression and replication. These findings reveal the dynamic features of the mRNA m⁶ A modification program during viral replication in vivo, especially in relation to key pathways involved in tumorigenesis.

Regulatory circular RNAs in viral diseases: applications in diagnosis and therapy

Circular RNA (circRNA) forms closed loops via back-splicing in precursor mRNA, resisting exonuclease degradation. In higher eukaryotes, protein-coding genes create circRNAs through exon back-splicing. Unlike mRNAs, circRNAs possess unique production and structural traits, bestowing distinct cellular functions and biomedical potential. In this review, we explore the pivotal roles of viral circRNAs and associated RNA in various biological processes. Analysing the interactions between viral circRNA and host cellular machinery yields fresh insights into antiviral immunity, catalysing the development of potential therapeutics. Furthermore, circRNAs serve as enduring biomarkers in viral diseases due to their stable translation within specific tissues. Additionally, a deeper understanding of translational circRNA could expedite the establishment of circRNA-based expression platforms, meeting the rising demand for broad-spectrum viral vaccines. We also highlight the applications of circular RNA in biomarker studies as well as circRNA-based therapeutics. Prospectively, we expect a technological revolution in combating viral infections using circRNA.

Analyzing viral epitranscriptomes using nanopore direct RNA sequencing

RNA modifications are a common occurrence across all domains of life. Several chemical modifications, including N6-methyladenosine, have also been found in viral transcripts and viral RNA genomes. Some of the modifications increase the viral replication efficiency while also helping the virus to evade the host immune system. Nonetheless, there are numerous examples in which the host's RNA modification enzymes function as antiviral factors. Although established methods like MeRIP-seq and miCLIP can provide a transcriptome- wide overview of how viral RNA is modified, it is difficult to distinguish between the complex overlapping viral transcript isoforms using the short read-based techniques. Nanopore direct RNA sequencing (DRS) provides both long reads and direct signal readings, which may carry information about the modifications. Here, we describe a refined protocol for analyzing the RNA modifications in viral transcriptomes using nanopore technology.

Role of Epitranscriptomic and Epigenetic Modifications during the Lytic and Latent Phases of Herpesvirus Infections

Herpesviruses are double-stranded DNA viruses occurring at a high prevalence in the human population and are responsible for a wide array of clinical manifestations and diseases, from mild to severe. These viruses are classified in three subfamilies (Alpha-, Beta- and Gammaherpesvirinae), with eight members currently known to infect humans. Importantly, all herpesviruses can establish lifelong latent infections with symptomatic or asymptomatic lytic reactivations. Accumulating evidence suggest that chemical modifications of viral RNA and DNA during the lytic and latent phases of the infections caused by these viruses, are likely to play relevant roles in key aspects of the life cycle of these viruses by modulating and regulating their replication, establishment of latency and evasion of the host antiviral response. Here, we review and discuss current evidence regarding epitranscriptomic and epigenetic modifications of herpesviruses and how these can influence their life cycles. While epitranscriptomic modifications such as m⁶A are the most studied to date and relate to positive effects over the replication of herpesviruses, epigenetic modifications of the viral genome are generally associated with defense mechanisms of the host cells to suppress viral gene transcription. However, herpesviruses can modulate these modifications to their own benefit to persist in the host, undergo latency and sporadically reactivate.

Impact of the Potential m6A Modification Sites at the 3′UTR of Alfalfa Mosaic Virus RNA3 in the Viral Infection

We have previously reported the presence of m6A in the AMV (Alfamovirus, Bromoviridae) genome. Interestingly, two of these putative m6A-sites are in hairpin (hp) structures in the 3’UTR of the viral RNA3. One site (2012AAACU2016) is in the loop of hpB, within the coat protein binding site 1 (CPB1), while the other (1900UGACC1904) is in the lower stem of hpE, a loop previously associated with AMV negative-strand RNA synthesis. In this work, we have performed in vivo experiments to assess the role of these two regions, containing the putative m6A-sites in the AMV cycle, by introducing compensatory point mutations to interfere with or abolish the m6A-tag of these sites. Our results suggest that the loop of hpB could be involved in viral replication/accumulation. Meanwhile, in the 1900UGACC1904 motif of the hpE, the maintenance of the adenosine residue and the lower stem hpE structure are necessary for in vivo plus-strand accumulation. These results extend our understanding of the requirements for hpE in the AMV infection cycle, indicating that both the residue identity and the base-pairing capacity in this structure are essential for viral accumulation.

Melatonin: Regulation of Viral Phase Separation and Epitranscriptomics in Post-Acute Sequelae of COVID-19

The relentless, protracted evolution of the SARS-CoV-2 virus imposes tremendous pressure on herd immunity and demands versatile adaptations by the human host genome to counter transcriptomic and epitranscriptomic alterations associated with a wide range of short- and long-term manifestations during acute infection and post-acute recovery, respectively. To promote viral replication during active infection and viral persistence, the SARS-CoV-2 envelope protein regulates host cell microenvironment including pH and ion concentrations to maintain a high oxidative environment that supports template switching, causing extensive mitochondrial damage and activation of pro-inflammatory cytokine signaling cascades. Oxidative stress and mitochondrial distress induce dynamic changes to both the host and viral RNA m6A methylome, and can trigger the derepression of long interspersed nuclear element 1 (LINE1), resulting in global hypomethylation, epigenetic changes, and genomic instability. The timely application of melatonin during early infection enhances host innate antiviral immune responses by preventing the formation of "viral factories" by nucleocapsid liquid-liquid phase separation that effectively blockades viral genome transcription and packaging, the disassembly of stress granules, and the sequestration of DEAD-box RNA helicases, including DDX3X, vital to immune signaling. Melatonin prevents membrane depolarization and protects cristae morphology to suppress glycolysis via antioxidant-dependent and -independent mechanisms. By restraining the derepression of LINE1 via multifaceted strategies, and maintaining the balance in m6A RNA modifications, melatonin could be the quintessential ancient molecule that significantly influences the outcome of the constant struggle between virus and host to gain transcriptomic and epitranscriptomic dominance over the host genome during acute infection and PASC.

The Epitranscriptome in miRNAs: Crosstalk, Detection, and Function in Cancer

The epitranscriptome encompasses all post-transcriptional modifications that occur on RNAs. These modifications can alter the function and regulation of their RNA targets, which, if dysregulated, result in various diseases and cancers. As with other RNAs, miRNAs are highly modified by epitranscriptomic modifications such as m6A methylation, 2'-O-methylation, m5C methylation, m7G methylation, polyuridine, and A-to-I editing. miRNAs are a class of small non-coding RNAs that regulates gene expression at the post-transcriptional level. miRNAs have gathered high clinical interest due to their role in disease, development, and cancer progression. Epitranscriptomic modifications alter the targeting, regulation, and biogenesis of miRNAs, increasing the complexity of miRNA regulation. In addition, emerging studies have revealed crosstalk between these modifications. In this review, we will summarize the epitranscriptomic modifications-focusing on those relevant to miRNAs-examine the recent crosstalk between these modifications, and give a perspective on how this crosstalk expands the complexity of miRNA biology.

Lessons Learned and Yet-to-Be Learned on the Importance of RNA Structure in SARS-CoV-2 Replication

SARS-CoV-2, the etiological agent responsible for the COVID-19 pandemic, is a member of the virus family Coronaviridae, known for relatively extensive (~30-kb) RNA genomes that not only encode for numerous proteins but are also capable of forming elaborate structures. As highlighted in this review, these structures perform critical functions in various steps of the viral life cycle, ultimately impacting pathogenesis and transmissibility. We examine these elements in the context of coronavirus evolutionary history and future directions for curbing the spread of SARS-CoV-2 and other potential human coronaviruses. While we focus on structures supported by a variety of biochemical, biophysical, and/or computational methods, we also touch here on recent evidence for novel structures in both protein-coding and noncoding regions of the genome, including an assessment of the potential role for RNA structure in the controversial finding of SARS-CoV-2 integration in “long COVID” patients. This review aims to serve as a consolidation of previous works on coronavirus and more recent investigation of SARS-CoV-2, emphasizing the need for improved understanding of the role of RNA structure in the evolution and adaptation of these human viruses.

Nanopore-Based Detection of Viral RNA Modifications

The chemical modification of ribonucleotides plays an integral role in the biology of diverse viruses and their eukaryotic host cells. Mapping the precise identity, location, and abundance of modified ribonucleotides remains a key goal of many studies aimed at characterizing the function and importance of a given modification. While mapping of specific RNA modifications through short-read sequencing approaches has powered a wealth of new discoveries in the past decade, this approach is limited by inherent biases and an absence of linkage information. Moreover, in viral contexts, the challenge is increased due to the compact nature of viral genomes giving rise to many overlapping transcript isoforms that cannot be adequately resolved using short-read sequencing approaches. The recent emergence of nanopore sequencing, specifically the ability to directly sequence native RNAs from virus-infected host cells, provides not just a new methodology for mapping modified ribonucleotides but also a new conceptual framework for what can be derived from the resulting sequencing data. In this minireview, we provide a detailed overview of how nanopore direct RNA sequencing works, the computational approaches applied to identify modified ribonucleotides, and the core concepts underlying both. We further highlight recent studies that have applied this approach to interrogating viral biology and finish by discussing key experimental considerations and how we predict that these methodologies will continue to evolve.

EBV Exploits RNA m6A Modification to Promote Cell Survival and Progeny Virus Production During Lytic Cycle

N6-methyladenosine (m⁶A) mediates various biological processes by affecting RNA stability, splicing, and translational efficiency. The roles of m⁶A modification in Epstein-Barr virus (EBV) infection in the lytic phase are unclear. Here, knockout of the m⁶A methyltransferase, N6-methyladenosine methyltransferase-like 3 (METTL3), or inhibition of methylation by DAA or UZH1a decreased the expression of viral lytic proteins and reduced progeny virion production. Interestingly, cell growth and viability were decreased by induction of the lytic cycle in METTL3-knockout or inhibitor-treated cells. Apoptosis was induced in those conditions possibly because of a decreased level of the anti-apoptotic viral protein, BHRF1. Therefore, m⁶A shows potential as a target of lytic induction therapy for EBV-associated cancers, including Burkitt lymphoma.

Pathogenic roles of m6A modification in viral infection and virus-driven carcinogenesis

N6-methyladenosine (m6A) is a prevalent modification of RNA in eukaryotes, bacteria, and viruses. It is highly conserved and can affect the structure, localization, and biology functions of RNA. In recent years, multiple m6A methylation sites have been identified in the viral RNA genome and transcripts of DNA viruses. This modification occurs commonly during the primary infection and is dynamically regulated by a methyltransferase (writers), demethylase (eraser) and m6A-binding proteins (readers) within the host cells. The abnormal m6A modification not only affects the replication of pathogenic viruses and host immune response but also contributes to the pathogenesis of virus-induced cancers. In this review, we highlight recent advances on the mechanism of m6A modification on viral replication, host immune response and carcinogenesis to provide a novel insight for epigenetic prevention of viral infection and virus-driven carcinogenesis.

Chemical biology and medicinal chemistry of RNA methyltransferases

RNA methyltransferases (MTases) are ubiquitous enzymes whose hitherto low profile in medicinal chemistry, contrasts with the surging interest in RNA methylation, the arguably most important aspect of the new field of epitranscriptomics. As MTases become validated as drug targets in all major fields of biomedicine, the development of small molecule compounds as tools and inhibitors is picking up considerable momentum, in academia as well as in biotech. Here we discuss the development of small molecules for two related aspects of chemical biology. Firstly, derivates of the ubiquitous cofactor S-adenosyl-l-methionine (SAM) are being developed as bioconjugation tools for targeted transfer of functional groups and labels to increasingly visible targets. Secondly, SAM-derived compounds are being investigated for their ability to act as inhibitors of RNA MTases. Drug development is moving from derivatives of cosubstrates towards higher generation compounds that may address allosteric sites in addition to the catalytic centre. Progress in assay development and screening techniques from medicinal chemistry have led to recent breakthroughs, e.g. in addressing human enzymes targeted for their role in cancer. Spurred by the current pandemic, new inhibitors against coronaviral MTases have emerged at a spectacular rate, including a repurposed drug which is now in clinical trial.

Identification and Characterization of BmNPV m6A Sites and Their Possible Roles During Viral Infection

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the most serious pathogens and causes serious economic losses in sericulture. At present, there is no epigenetic modification of BmNPV transcripts, especially of m6A, and this modification mediates diverse cellular and viral functions. This study showed that m6A modifications are widespread in BmNPV transcripts in virally infected cells and the identified m6A peaks with a conserved RRACH sequence. m6A sites predominantly appear in the coding sequences (CDS) and the 3'-end of CDS. About 37% of viral genes with m6A sites deleted from the viral genome did not produce any infectious virions in KOV-transfected cells. Among the viral genes related to replication and proliferation, ie-1 mRNA was identified with a higher m6A level than other viral genes. The m6A sites in the ie-1 mRNA may be negatively related to the protein expression. Viral replication was markedly inhibited in cells overexpressed with BmYTHDF3 in a dose-dependent manner, and a contrary effect was found in si-BmYTHDF3-transfected cells. Collectively, the identification of putative m6A modification in BmNPV transcripts provides a foundation for comprehensively understanding the viral infection, replication, and pathobiology in silkworms.

Comparative Virus-Host Protein Interactions of the Bluetongue Virus NS4 Virulence Factor

Bluetongue virus (BTV) is the etiologic agent of a non-contagious arthropod-borne disease transmitted to wild and domestic ruminants. BTV induces a large panel of clinical manifestations ranging from asymptomatic infection to lethal hemorrhagic fever. Despite the fact that BTV has been studied extensively, we still have little understanding of the molecular determinants of BTV virulence. In our report, we have performed a comparative yeast two-hybrid (Y2H) screening approach to search direct cellular targets of the NS4 virulence factor encoded by two different serotypes of BTV: BTV8 and BTV27. This led to identifying Wilms’ tumor 1-associated protein (WTAP) as a new interactor of the BTV-NS4. In contrast to BTV8, 1, 4 and 25, NS4 proteins from BTV27 and BTV30 are unable to interact with WTAP. This interaction with WTAP is carried by a peptide of 34 amino acids (NS422−55) within its putative coil-coiled structure. Most importantly, we showed that binding to WTAP is restored with a chimeric protein where BTV27-NS4 is substituted by BTV8-NS4 in the region encompassing residue 22 to 55. We also demonstrated that WTAP silencing reduces viral titers and the expression of viral proteins, suggesting that BTV-NS4 targets a cellular function of WTAP to increase its viral replication.

Shaping the Innate Immune Response Through Post-Transcriptional Regulation of Gene Expression Mediated by RNA-Binding Proteins

Innate immunity is the frontline of defense against infections and tissue damage. It is a fast and semi-specific response involving a myriad of processes essential for protecting the organism. These reactions promote the clearance of danger by activating, among others, an inflammatory response, the complement cascade and by recruiting the adaptive immunity. Any disequilibrium in this functional balance can lead to either inflammation-mediated tissue damage or defense inefficiency. A dynamic and coordinated gene expression program lies at the heart of the innate immune response. This expression program varies depending on the cell-type and the specific danger signal encountered by the cell and involves multiple layers of regulation. While these are achieved mainly via transcriptional control of gene expression, numerous post-transcriptional regulatory pathways involving RNA-binding proteins (RBPs) and other effectors play a critical role in its fine-tuning. Alternative splicing, translational control and mRNA stability have been shown to be tightly regulated during the innate immune response and participate in modulating gene expression in a global or gene specific manner. More recently, microRNAs assisting RBPs and post-transcriptional modification of RNA bases are also emerging as essential players of the innate immune process. In this review, we highlight the numerous roles played by specific RNA-binding effectors in mediating post-transcriptional control of gene expression to shape innate immunity.

Direct RNA Sequencing Reveals SARS-CoV-2 m6A Sites and Possible Differential DRACH Motif Methylation among Variants

The causative agent of COVID-19 pandemic, SARS-CoV-2, has a 29,903 bases positive-sense single-stranded RNA genome. RNAs exhibit about 150 modified bases that are essential for proper function. Among internal modified bases, the N6-methyladenosine, or m6A, is the most frequent, and is implicated in SARS-CoV-2 immune response evasion. Although the SARS-CoV-2 genome is RNA, almost all genomes sequenced thus far are, in fact, reverse transcribed complementary DNAs. This process reduces the true complexity of these viral genomes because the incorporation of dNTPs hides RNA base modifications. Here, we present an initial exploration of Nanopore direct RNA sequencing to assess the m6A residues in the SARS-CoV-2 sequences of ORF3a, E, M, ORF6, ORF7a, ORF7b, ORF8, N, ORF10 and the 3'-untranslated region. We identified fifteen m6A methylated positions, of which, six are in ORF N. Additionally, because m6A is associated with the DRACH motif, we compared its distribution in major SARS-CoV-2 variants. Although DRACH is highly conserved among variants, we show that variants Beta and Eta have a fourth position C > U change in DRACH at 28,884b that could affect methylation. This is the first report of direct RNA sequencing of a Brazilian SARS-CoV-2 sample coupled with the identification of modified bases.

The m6A landscape of polyadenylated nuclear (PAN) RNA and its related methylome in the context of KSHV replication

Polyadenylated nuclear (PAN) RNA is a long noncoding transcript involved in Kaposi's sarcoma-associated herpesvirus (KSHV) lytic reactivation and regulation of cellular and viral gene expression. We have previously shown that PAN RNA has dynamic secondary structure and protein binding profiles that can be influenced by epitranscriptomic modifications. N6-methyladenosine (m6A) is one of the most abundant chemical signatures found in viral RNA genomes and virus-encoded RNAs. Here, we combined antibody-independent next-generation mapping with direct RNA sequencing to address the epitranscriptomic status of PAN RNA in KSHV infected cells. We showed that PAN m6A status is dynamic, reaching the highest number of modifications at the late lytic stages of KSHV infection. Using a newly developed method, termed selenium-modified deoxythymidine triphosphate (SedTTP)-reverse transcription (RT) and ligation assisted PCR analysis of m6A (SLAP), we gained insight into the fraction of modification at identified sites. By applying comprehensive proteomic approaches, we identified writers and erasers that regulate the m6A status of PAN, and readers that can convey PAN m6A phenotypic effects. We verified the temporal and spatial subcellular availability of the methylome components for PAN modification by performing confocal microscopy analysis. Additionally, the RNA biochemical probing (SHAPE-MaP) outlined local and global structural alterations invoked by m6A in the context of full-length PAN RNA. This work represents the first comprehensive overview of the dynamic interplay that takes place between the cellular epitranscriptomic machinery and a specific viral RNA in the context of KSHV infected cells.