Amin Izadpanah
Jay Rappaport
Prasun K. Datta- 1Division of Comparative Pathology, Tulane National Primate Center, Covington, LA, United States
- 2Department of Microbiology and Immunology, School of Medicine, Tulane University, New Orleans, LA, United States
Recent studies on the epitranscriptomic code of SARS-CoV-2 infection have discovered various RNA modifications, such as N6-methyladenosine (m6A), pseudouridine (Ψ), and 2′-O-methylation (Nm). The effects of RNA methylation on SARS-CoV-2 replication and the enzymes involved in this mechanism are emerging. In this review, we summarize the advances in this emerging field and discuss the role of various players such as readers, writers, and erasers in m6A RNA methylation, the role of pseudouridine synthase one and seven in epitranscriptomic modification Ψ, an isomer of uridine, and role of nsp16/nsp10 heterodimer in 2′-O-methylation of the ribose sugar of the first nucleotide of SARS-CoV-2 mRNA. We also discuss RNA expression levels of various enzymes involved in RNA modifications in blood cells of SARS-CoV-2 infected individuals and their impact on host mRNA modification. In conclusion, these observations will facilitate the development of novel strategies and therapeutics for targeting RNA modification of SARS-CoV-2 RNA to control SARS-CoV-2 infection.
Introduction
The etiological agent of coronavirus disease-2019 (COVID-19) is the novel beta coronavirus, severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) (Datta et al., 2020; Wu et al., 2020; Zhou et al., 2020). To date, this virus has infected 430 million people and killed more than five million people worldwide (https://coronavirus.jhu.edu/data). The life cycle of SARS-CoV-2 initiates with the viral envelope protein spike (S) binding predominantly the cellular receptor angiotensin-converting enzyme 2 (ACE2) and other potential alternate ACE2 independent receptors, such as CD147, Tyrosine-protein kinase receptor UFO (AXL), Kringle-containing protein marking the eye and the nose protein 1(KREMEN1) and Asialoglycoprotein receptor 1 (ASGR1) (Wang et al., 2020; Wang et al., 2021; Gu et al., 2022). Viral entry via ACE2 dependent mechanism also requires other cellular proteases such as serine proteases (TMPRSS2, TMPRSS11D, and TMPRSS13), furin, and cathepsin L (CTSL) that are involved in the activation of the spike protein through proteolytic cleavage (Bestle et al., 2020; Hashimoto et al., 2021; Hoffmann et al., 2020; Mellott et al., 2021). The release of the viral genome from the endosomal compartment into the cytoplasm is a prerequisite for initiating the viral replication cycle. During the life cycle of SARS-CoV-2, viral proteins and nucleic acid closely interact with many host proteins to regulate viral replication (Gordon et al., 2020; Daniloski et al., 2021; Schmidt et al., 2021). SARS-CoV-2 can alter host nucleic acids and proteins to promote successful viral replication and impair or shut cellular responses to infection (Banerjee et al., 2020).
Many studies have shown that most eukaryotic mRNAs, rRNAs, and tRNAs (Li and Mason, 2014; Roundtree et al., 2017), and viral RNAs (Dang et al., 2019; Tsai and Cullen, 2020; Imam et al., 2020) have multiple forms of RNA modifications that together are defined as the ‘Epitranscriptome’ (Saletore et al., 2012). To date, more than 170 RNA modifications have been identified (Wiener and Schwartz, 2021). N6-methyladenosine (m6A) is the most characterized RNA post-transcriptional modifications (PTM), the other modifications that have been identified are 5-methyl cytosine (m5C), 5-hydroxymethylcytosine (hm5C), 7-methylguanosine (m7G), 1-methylguanosine (m1G), pseudouridine (Ψ), N6, N6-dimethyladenosine (m6,2A), ribose-methylation (2′-O-Me), uridine (U), and inosine (I) (Shi et al., 2020). Elucidation of these modifications’ functional roles in RNAs has shown that they play a role in nuclear export of RNAs, pre-splicing of mRNA, stability of RNA, translation initiation, and viral infection (Shi et al., 2020). Concomitantly, several techniques have emerged that enable the identification and study of significant modifications, namely m6A m5C, and Ψ (Sarkar et al., 2021; Wiener and Schwartz, 2021).
SARS-CoV-2 Transcriptome
The novel SARS-CoV-2 virus contains a positive-sense, single-stranded RNA genome of ∼30 kb (Wu et al., 2020). The positive sense nature of the genome enables immediate translation, which produces two polypeptides named pp1a (440–500 kDa) or pp1ab (740–810 kDa). The polypeptide pp1a is derived from the open reading frame-1 (ORF1), while the pp1ab is derived from the ORF1ab. ORF1ab occurs due to a -1-ribosome frameshift signal upstream of the ORF1a stop codon, enabling continued translation and thus the larger pp1ab. The polypeptides pp1a and pp1ab are cleaved by viral proteases nsp3 and nsp5 to yield 11 or 15 non-structural proteins (nsps), respectively. Both structural proteins and Nsps depend on translation within the host, but Nsps are viral proteins that are not packaged inside the virus. The nsp12 encodes an RNA-dependent RNA polymerase (RdRp) essential to the coronavirus replication cycle. The RdRp uses the positive-sense genomic template to generate negative-strand RNA. These negative-strand RNA intermediates become the template for producing genomic RNA and positive-sense sub-genomic RNAs (sgRNAs). Genomic RNAs are copies of the original genome of the coronavirus and are thus packaged into the progeny virion (Figure 1A). The sgRNAs contain the same 3′ sequence as the genomic RNA; however, the 5′ ends are different from that of genomic RNA. In coronaviruses, a discontinuous transcription mode occurs during negative-sense RNA synthesis, referred to as leader-body fusion (Sola et al., 2015; Miller and Koev, 2000). Beginning from the 3′ end of the genomic RNA, RdRp transcribes the body and is halted at the transcription-regulatory sequence in the body (TRS-B). RdRp then resumes transcription at the TRS-L (transcription-regulatory sequence at the leader sequence), a more 5′ location on the genomic RNA. The TRS is located next to ORFs. Therefore, the positive sense mRNA (sgRNA) produced from the negative-sense RNA contains the leader-sequence fused to the distal ORF, for example, the S, E, M, and N proteins (Figure 1B). Importantly, sgRNAs encode structural proteins in SARS-CoV-2, including spike protein (S), an envelope protein (E), membrane protein (M), and nucleocapsid protein (N), and accessory proteins 3a, 6, 7a, 7b, 8, and 10.
FIGURE 1. (A) Schematic representation of the SARS-CoV-2 biogenesis. SARS-CoV-2 enters host cells interacting with the angiotensin-converting enzyme 2 (ACE2) receptor by the surface spike (S) protein. Upon entry of the virus into the host cell, viral genomic RNA is released into the cytoplasm, where it is translated into viral polymerase proteins. Here, sub-genomic (–) RNAs are synthesized and used as templates for sub-genomic (+) messenger RNAs (mRNAs). The nucleocapsid (N) structural protein and viral RNA are replicated, transcribed, and synthesized in the cytoplasm. In contrast, other viral structural proteins, including the S protein, membrane (M) protein, and envelope (E) protein, are transcribed and then translated in the endoplasmic reticulum (ER). The structural proteins traverse the ER-Golgi intermediate compartment for virion assembly, followed by the release of the nascent virion from the host cell via exocytosis. (B) Schematic representation of SARS-CoV-2 sgRNA synthesis. SARS CoV-2 sgRNAs are synthesized via discontinuous transcription. Beginning from the 3′ end of the genomic RNA, RdRp transcribes the body and is halted at the transcription-regulatory sequence in the body (TRS-B). RdRp then resumes transcription at the TRS-L (transcription-regulatory sequence at the leader sequence), a more 5′ location on the genomic RNA. Next, the newly generated negative strand is used as a template for positive-strand synthesis. The TRS is located next to ORFs. Therefore, the positive sense mRNA (sgRNA) produced from the negative-sense RNA contains the leader-sequence fused to the distal ORF. Created with BioRender.
Methods for Identification of RNA Modifications in SARS-CoV-2
Several methodologies have been developed to identify RNA modifications (mostly m6A) in eukaryotic cells and viruses. We only describe and compare the techniques used so far to detect SARS-CoV-2 modifications for brevity.
Antibody-Based Detection Methods
meRIP-seq/m6A-seq: (Dominissini et al., 2012; Meyer et al., 2012). In this method, m6A-specific antibodies are used to immunoprecipitate total RNA. Following isolation of m6A bound RNA, the RNA is reverse transcribed to cDNA and deep sequenced using next-generation sequencing (NGS) protocols to obtain high-resolution reads of m6A methylated RNA. The limitations of this technique are 1) issues with m6A antibody specificity, and 2) read lengths are ∼100–200 nucleotides wide, and the bioinformatic prediction of m6A residues is limited to one site per peak.
miCLIP (Linder et al., 2015). In this method, the RNA is first sheared and incubated with anti-m6A antibodies to cross-link RNA using UV light. m6A antibody bound-RNA complexes are recovered by protein A/G-affinity purification, followed by SDS-PAGE and transfer to nitrocellulose membrane. The RNA is then released from the membrane by proteinase K and reverse transcribed to generate a cDNA library using the iCLIP protocol (Hafner et al., 2010). The resulting cDNA is PCR-amplified before sequencing using NGS protocols. The limitations of this technique are 1) issues with m6A antibody specificity, 2) use of excess RNA, and 3) high probability for the introduction of mutation or truncation in relation to the position of the modified adenosine by m6A antibody.
Liquid chromatography-Tandem Mass Spectrometry (LC-MS/MS)
Oligonucleotide or nucleoside LC/MS (Jora et al., 2019); High-Resolution Ion Mobility Spectrometry-Mass Spectrometry (Kenderdine et al., 2020). In this technique, the total RNA or purified mRNA is digested into individual nucleotides, separated by LC, and quantified by MS. The MS peaks from the sample are compared with the MS peaks of standards to assess all nucleotide modifications in an RNA sample. The limitations of this technique are 1) it requires large amounts of input RNA, 2) no information can be obtained about the location of the modification in an RNA molecule. The advantage of this technique is that even if the levels of the modified nucleotides are low, the measurements are quantitative and reproducible between studies.
Direct RNA Nanopore Sequencing (DRS)
An alternative to sequencing-by-synthesis is the DRS platform developed by Oxford Nanopore Technologies (Loman et al., 2015). This technique involves ratcheting a single RNA molecule tethered to a motor protein through a protein nanopore sensor (E. coli CsgG-derived nanopore) embedded in a synthetic hydrophobic membrane. The sensor measures changes in an ionic current as RNA passes through the nanopore; information about changes in current and dwell time in the pore identifies the unmodified and modified nucleotide.
Compared to antibody-based methods, the strengths of this technique enable one to detect any RNA modifications in full-length native RNAs at single-nucleotide resolution without the need for reverse transcription or PCR amplification (Depledge and Wilson, 2020; Leger et al., 2021). The limitations are 1) high input RNA required and 2) low throughput.
RNA Modifications in SARS-CoV-2 RNA
The most predominant RNA modifications are N6-methyladenosine (m6A), N7-methylguanosine (m7G), 5-methylcytosine (m5C), 5-hydroxymethylcytosine (hm5C), pseudouridine (Ψ) and 2ʹ-O-methylations (Nm) (Figure 2). Three different RNA modifications were found in SARS-CoV-2 RNA, namely, m6A, Ψ, and 2′-O-methylation.
FIGURE 2. Schematic representation of the common RNA modifications in mRNA. Created with Biorender.
The RNA methyltransferase complex is responsible for m6A modification in RNA. This complex is mainly composed of m6A writers METTL3 (methyltransferase like 3), METTL14 (methyltransferase like 14), and WTAP (Wilm’s tumor 1-associated protein). METTL14 forms a heterodimer with METTL3; whereas WTAP is a part of the core complex and influences m6A deposition by METTL3-METTL14 (Little et al., 2000; Bujnicki et al., 2002; Liu et al., 2014; Ping et al., 2014). RBM15 (RNA-binding motif protein 15) plays a role as an adaptor protein essential for the initial recruitment of the writers onto pre-mRNAs (Patil et al., 2016). VIRMA (vir-like m6A methyltransferase associated), an adaptor protein, facilitates interactions with other proteins of the RNA methyltransferase complex (Schwartz et al., 2014). ZC3H13 (Zinc finger CCCH domain-containing protein 13) is a protein that interferes positively with the binding between the adaptor protein RB15 and WTAP (Knuckles et al., 2018).
The addition of Ψ to RNA or pseudouridylation is catalyzed by a class of proteins called pseudouridine synthases (PUS). In humans, there are 13 pseudouridine synthases (Borchardt et al., 2020). Studies have shown that PUS1, PUS7, and TRUB1 pseudouridylate mRNAs (Safra et al., 2017; Carlile et al., 2019).
Viral RNAs predominantly undergo 2′-O-methylation using methyltransferase encoded by the virus itself (Netzband and Pager, 2020).
The first study by Kim and co-workers used the nanopore-based direct RNA sequencing (DRS) approach to study the SARS-CoV-2 transcriptome (Kim et al., 2020). The DRS approach was combined with the complementary DNA nano ball sequencing method to enable a robust investigation of the SARS-CoV-2 transcriptome. Vero cells were infected with SARS-CoV-2 (BetaCoV/Korea/KCDC03/2020). RNA-seq revealed significant read abundance from the 3′ end of the genome, corresponding to the sgRNAs. The high sequencing depth of this method also enabled evaluation of the junctions between the 5′ and 3′ breakpoints, which confirmed the TRS mediated leader-body fusion mechanism of discontinuous transcription for sgRNA synthesis. The transcript frequencies in descending order are N, S, 7a, 3a, 8, M, E, 6, and 7b. Interestingly, Kim et al. identified at least 41 potential RNA modification sites with an AAGAA motif. It was also shown that long viral RNAs (such as gRNA, S, 3a, E, and M) were more frequently modified than the shorter ones (6, 7a, 7b, 8, and N), suggesting that the mechanism of RNA modification may be RNA specific (Kim et al., 2020). Polyadenylation analysis revealed that SARS-CoV-2 sgRNA transcripts have two distinct polyA tail lengths of 30 nucleotides (nt) and 45 nt. These varying lengths may reflect time post-infection. It was shown that polyA tails of bovine coronavirus mRNA increase from 45 nt immediately post-infection to 65 nt 6–9 h post-infection decrease to 30 nt at 120 h post-infection (Wu et al., 2013). Therefore, the authors theorized that the 30 nt polyA tail might reflect aged and decay-prone RNAs (Kim et al., 2020). Forty-one potential sites of base modification (detected by ionic current changes) were correlated with shorter polyA tails, although the type of RNA modification was unknown (Kim et al., 2020).
N6-Methyladenosine (m6A)
The most common post-transcriptional RNA modification is the N6-methyladenosine (m6A). The m6A modification is commonly found on the DRACH motif (D = A, G, U; R = A, G; H = A, C, U) (Linder et al., 2015). The studies by Liu and co-workers profiled the m6A methylome of SARS-CoV-2 RNA in Vero and Huh7 cells (Liu et al., 2021). To identify the m6A methylome, a refined RNA immunoprecipitation (RIP) strategy was used. Total RNA was subjected to m6A antibody treatment for immunoprecipitation and subsequently sequenced (RIP-seq). In Vero cells, RIP-seq revealed four m6A peaks at 24 h post-infection (hpi) and nine additional m6A peaks at 56hpi. These suggest that m6A modification increases at a later point of infection, supported by increased intensity of m6A at 56 hpi vs. 24 hpi. In Huh7 cells, RIP-seq detected six m6A peaks, corresponding with those in Vero cells at 56 hpi. To determine specific m6A sites, modified m6A individual-nucleotide-resolution crosslinking and miCLIP (Linder et al., 2015) were employed. The studies identified three m6A sites in ORF1ab, one m6A site in ORF 7a, three m6A sites in N, and one m6A site in ORF 10 (Liu et al., 2021). The propensity of the m6A sites to cluster at the 3′ end suggests the involvement of the sgRNAs. Further analysis revealed correlative mutations for each of the m6A sites, except for the second site in ORF1ab. Mutations localizing in the first and third m6A sites (both in ORF1ab) were mainly European strains, while mutations localizing in the fourth and sixth m6A sites (in ORF7a and N, respectively) were primarily in North American isolates. Importantly, m6A sites were also found in the negative-sense strand RNA; however, the ratio of negative-sense to positive sense m6A sites was <1%, and miCLIP could not yield accurate localization. m6A sites on negative-sense RNA suggest a dynamic role for methylation during viral infection and replication.
Using MeRIP-seq to detect m6A modification in the SARS-CoV-2 RNA isolated from infected Vero E6 cells (Zhang et al., 2021), five m6A peaks were identified in the 5′ end (nucleotide 36 to 753 and nucleotide 1,023–1,324) and the 3′ end (nucleotide 27,493 to 27,913, nucleotide 28,475 to 28,706, and nucleotide 28,944–29,751). The m6A residues were in the ORF1ab-, N-, and ORF10-coding regions of the SARS-CoV-2 genome (Zhang et al., 2021). Additional experiments using nanopore-based direct RNA sequencing (DRS) confirmed the specific m6A modification sites in RNAs extracted from Vero E6, A549-ACE2, and Huh7 cells infected with SARS-CoV-2. In brief, six m6A sites were found to be conserved in all the infected cell lines, and the m6A motif in the SARS-CoV-2 genome was mainly GGACA (Zhang et al., 2021). Analysis of SARS-CoV-2 RNA in A549-ACE2 cells and Vero E6 cells revealed nine common m6A sites. However, three m6A peaks in SARS-CoV-2 from Huh7 cells differed from those in Vero E6 and A549-ACE2 cells SARS-CoV-2 RNA. These results show that both conserved and different m6A sites exist in the SARS-CoV-2 genome in other cell lines.
Studies by Campos et al. (Campos et al., 2021) employed DRS to assess the location of m6A residues in SARS-CoV-2 RNA isolated from supernatants of SARS-CoV-2-infected Vero E6 cells that are enriched for genomic RNAs and not sgRNA. They identified fifteen m6A methylated positions, of which six are in ORF N, three in ORF3a, one each in E, M, and ORF7a, and three in ORF7b. In addition, the studies also showed that m6A is associated with the DRACH motif that is highly conserved among variants. However, since the variants Beta (B.1.351) and Eta (B.1.525) have a fourth position C > U change in DRACH at nucleotide position 28,884 of SARS-CoV-2 RNA, this may affect methylation (Campos et al., 2021).
Using MACS2 for peak-calling (Zhang et al., 2008) of meRIP-seq data of SARS-CoV-2 RNA, the studies by Burgess and co-workers (Burgess et al., 2021) identified 14 peak regions in genome length transcripts (gRNA) only, suggesting that both gRNA and sgRNAs harbor m6A residues.
The studies by Li and co-workers (Li et al., 2021) employed liquid chromatography-tandem mass spectrometry (LC-MS/MS and MS/MS/MS) to profile the RNA modifications in SARS-CoV-2 RNA isolated from Vero E6 infected cells. Numerous RNA modifications were detected in SARS-CoV-2 RNA that includes 2′-O-methylated derivatives of all four canonical nucleosides (Am, Cm, Um, and Gm), modified cytidine derivatives (ac4C, m3C, and m5C), two modified uridine derivatives (J and m5U), and two modified adenosines (m6A and m6,6A), and 2-thiocytidine. The authors estimated that m6A accounted for 0.096% of adenosine in the virus and estimated eight m6A-modified sites. The distribution of m6A residues was further validated using MeRIP-seq of full-length SARS-CoV-2 RNA purified from Vero E6 cells. The authors reported that m6A peaks were present in the ORF1ab and 3′ end of the SARS-CoV-2 genome, especially in the N region of SARS-CoV-2 RNA isolated from infected Vero and Caco-2 cells (Li et al., 2021).
Reanalysis of the 41 modified sites reported by Kim and co-workers (Kim et al., 2020) in subsequent studies (Li et al., 2021) found modified adenosine sites embedded in the DRACH motif (D = G/A/U, R = G/A, H = A/U/C) in more than half of the modified sites that were excluded from the analysis of AAGAA motif. This observation revealed the existence of m6A modification in the BetaCoV/Korea/KCDC03/2020 viral genome, which was not reported earlier. The results of the above investigations are summarized in Table 1.
TABLE 1. Summary of the studies on SARS-CoV-2 m6A modifications.
Pseudouridine (Ψ)
Pseudouridine (Ψ), residue formed by isomerization of uridine (U), is found at high levels (>1%) in eukaryotic rRNA and tRNA and at lower levels (<1%) in eukaryotic mRNA (Carlile et al., 2014; Borchardt et al., 2020). Ψ residues are also seen in viral RNAs (McIntyre et al., 2018; Furuse, 2021). Using DRS, five high confidence Ψ sites in TRS-S, five in TRS-3a, five in TRS-E, and five in TRS-M of SARS-CoV-2 sgRNAs were detected (Fleming et al., 2021). In TRS-3a, -E, and -M, five identified peaks mapped to the nucleotide (U27164, U28039, U28759, U28927, and U29418). Biochemical validation showed that recombinant PUS1 and PUS7 introduced Ψ in synthetic RNAs at positions U28927 and U29418, while recombinant PUS1 introduced Ψ in synthetic RNAs at position U29418 (Fleming et al., 2021).
2′-O-Methylation
Numerous studies on RNA capping in coronaviruses showed the involvement of several non-structural proteins (nsps): nsp13, a bifunctional helicase and RNA/NTP triphosphatase; nsp14, a bifunctional mRNA cap guanine-N7 methyltransferase and 3′→5′ mismatch exonuclease; nsp16, a cap ribose 2′-O methyltransferase; and a guanylyltransferase (Ivanov and Ziebuhr, 2004; Bouvet et al., 2010; Chen et al., 2011; Bouvet et al., 2012). 2′-O methylation occurs on the ribose sugar of the first nucleotide of SARS-CoV-2 mRNA. It is catalyzed by a complex of SARS-CoV-2 nsp16 and nsp10 in the presence of cognate RNA substrate analog and methyl donor, S-adenosyl methionine (Viswanathan et al., 2020). In addition, the high-resolution structure of the ternary complex of the SARS-CoV-2 RNA cap/nsp16/nsp10 complex showed that the ligand-binding site in nsp16/10 allows accommodation of small molecules outside of the catalytic pocket (Viswanathan et al., 2020). In a recent study using Nm-seq, a total of 130 2′-O methylation sites were identified in the SARS-CoV-2 genome. The 2′-O methylation was enriched in the 5′ and 3′ UTRs of SARS-CoV-2 (Yang et al., 2021).
m6A Modification in SARS-CoV-2: Role of Writers, Erasers, and Readers
The writer enzymes add the m6A modification, recognized by the reader enzymes and removed by the eraser enzymes. Readers play a critical role in identifying the m6A mark and regulating the fate of m6A-marked mRNA; for brevity, we only discuss the readers, writers, and erasers that play a role in SARS-CoV-2 RNA modifications. Characterization of the writer protein complex involved in m6A modification in mRNA is methyl-transferase-like three METTL3 and its homolog METTL14 (Bujnicki et al., 2002; Liu et al., 2014). Both proteins act synergistically by forming a stable heterodimer METTL3 and METTL14 (1:1 stoichiometry) to facilitate m6A addition to mRNA methylation in cells (Liu et al., 2014). More recently, a third critical component of this complex, Wilm’s tumor-associated protein (WTAP), which lacks methyltransferase activity, was essential for the methylation process (Little et al., 2000; Ping et al., 2014). WTAP interacts with the METTL3/METTL14 heterodimer to direct the localization of this m6A writer complex to nuclear speckles, where splicing occurs (Liu et al., 2014). Interestingly, RNA-binding proteins, RNA-binding motif protein 15 (RBM15), and its paralogue RBM15B were shown to direct the writer complex to the XIST long non-coding RNA and several cellular mRNAs to guide site-specific methylation of target RNA (Patil et al., 2016). Zinc finger CCCH domain-containing protein 13 (ZC3H13) was identified as a nuclear m6A writer gene (Wen et al., 2018).
Among the m6A erasers, the first protein identified was fat mass and obesity protein (FTO), an alpha-ketoglutarate-dependent dioxygenase member of the AlkB family that demethylates RNA (Jia et al., 2011). AlkB homolog 5, RNA demethylase (ALKBH5), another member of the AlkB family, also promotes mRNA demethylation rates similar to FTO (Zheng et al., 2013).
The cytosolic m6A readers are the YT521-B homology domain-containing proteins (YTHDF1-3) that regulate target mRNA fate. YTHDF1-3 have a highly conserved single-stranded RNA-binding domain at the carboxy terminus (the YTH domain) and a less conserved amino-terminal region (Zhang et al., 2010). YTHDF1 enhances ribosome assembly of m6A mRNA and interacts with the translation initiation factor 3 (eIF3) to promote translation of m6A containing mRNA (Wang et al., 2015). YTHDF2 promotes the degradation of non-translating m6A-modified mRNAs (Wang X et al., 2014). YTHDF3 promotes the translation and degradation of mRNA (Shi et al., 2017). YTH domain-containing 1 (YTHDC1) is a nuclear m6A reader protein that mediates mRNA metabolism and regulates mRNA splicing (Xiao et al., 2016). YTH domain-containing 2 (YTHDC2) is an m6A reader protein with a 3′-5′ RNA helicase activity, binds to m6A mRNA, and regulates mRNA translation and stability (Wojtas et al., 2017). Insulin-like growth factor-2 (IGF2) mRNA-binding proteins 1, 2, and 3 (IGF2BP1/2/3) are a new family of m6A readers that prevents decay of m6A-modified mRNAs (Nielsen et al., 1999; Bell et al., 2013; Huang et al., 2018).
To explain the role of writers, erasers, and readers in the SARS-CoV-2 biogenesis, the studies by Zhang and co-workers (Zhang et al., 2021) showed that METTL3 overexpression in Vero E6 cells induced an abundance of m6A containing viral RNA levels, while METLL3 or FTO knockdown significantly decreased or increased virus levels, respectively. This study also showed that SARS-CoV-2 infection induced the expression of METTL3 and altered distribution in both the nucleus and cytoplasm and co-localized with RdRp. In addition, METTL14, WTAP, ALKBH5, and FTO were shown to co-localize with the viral protein N. These observations taken together suggest that SARS-CoV-2 infection affects host cells m6A methyltransferase and demethylases.
Using A549 cells expressing human ACE2 receptor and a SARS-CoV-2 reporter virus expressing mNeonGreen (icSARS-CoV-2-mNG) as a model of SARS-CoV-2 infection, the effects of knockdown of METTL3, YTHDF1, YTHDF2, and YTHDF3 was used tested (Burgess et al., 2021). The studies showed that siRNAs against METTL3 reduced the percentage of infected cells by 78–81% compared to control cells. While, (knock-down of YTHDF2 or YTHDF3 reduced the percentage of infected cells by 42–66% or 75–76% of the control, respectively. In contrast, the knockdown of YTHDF1 reduced SARS-CoV-2 infection by 23–89%, dependent upon the targeting siRNAs used. Co-depleting all three YTHDF1, YTHDF2, and YTHDF3 reduced the number of infected cells by 72% (Burgess et al., 2021).
Liu and co-workers further investigated the role of METTL3, METTL14, ALKBH5, and YTHDF2 by generating individual knockdown cell lines (Liu et al., 2021). Knockdown of METTL3, METTL14, and YTHDF2 significantly increased viral infection and replication, while knockdown of ALKBH5 decreased viral infection and replication. YTHDF2 was shown to induce the decay of m6A transcripts. Therefore, the authors concluded that host m6A methylome regulators play a critical role in SARS-CoV-2 replication. In the studies by Li and co-workers (Li et al., 2021), knockdown of METTL3 using two different small hairpin RNAs in Caco-2 cells infected with SARS-CoV-2 resulted in a reduction of SARS-CoV-2 viral load and proviral gene expression. MeRIP-seq analysis of SARS-CoV-2 isolated from METTL3 knock-down (METTL3-KD) Caco-2 cells showed decreased m6A occupancy in the N region of the SARS-CoV-2 virus (Li et al., 2021). Li et al. also demonstrate that the presence of m6A on the SARS-CoV-2 RNA avoids sensing by RIG-1 (retinoic acid-inducible gene 1) (Li et al., 2021) (Figure 3). RIG-1 is a cytosolic pattern recognition receptor essential to innate immunity through its ability to activate the type 1 interferon response (Figure 3).
FIGURE 3. m6A modification in SARS-CoV-2 RNA and cellular mRNA by METLL3 and RBM15, respectively, and its effect on evasion of host innate immune response by SARS-CoV-2 and host cell. (Left) After entering the host cell, the viral genome enters the replication phase, and METTL3 introduces m6A residues in viral RNA. m6A methylated SARS-CoV-2 does not bind RIG-I and induces the expression of antiviral molecules such as IFNβ. Viral infection also enhances METTL3 and RBM15, leading to enhanced host mRNA m6A methylation resulting in activation of inflammatory gene expression and programmed cell death. (Right) When antiviral molecules such as IFNβ levels increase, it leads to activation JAK-STAT pathway and interferon-stimulated gene (ISG) expression. Created with BioRender.
Collectively, these studies demonstrate that a functioning m6A RNA modification pathway is beneficial to SARS-CoV-2 replication. Further, it provides proof of concept that targeting cellular components of this intensively studied RNA modification pathway could ultimately lead to new therapeutic opportunities to control these important viral pathogens.
Impact of SARS-CoV-2 on Host Cell m6A Methylome
Many studies have documented that RNA modification of viral RNAs through methylation plays a role in the evasion of the host’s innate immune response (Kariko et al., 2005; Durbin et al., 2016; Lu et al., 2020). Conversely, viral infection may also play a role in regulating host methylation. Indeed, in recent studies, SARS-CoV-2 infection influenced the host cell m6A methylome. RNA seq analysis of METTL3-KD SARS-CoV-2 infected Caco-2 cells showed differential expression of genes involved in regulating metabolic processes, immune response, RNA processing, ribosome function, protein modification, DNA repair, cell cycle, and cell differentiation (Li et al., 2021). However, knockdown of METTL3 in Calu-3 cells infected with SARS-CoV-2 significantly enhanced IFNβ and ISG expression (STAT1, STAT2, and IRF7) and cytokine/chemokine expression (IL-6, IL-8, TNF, CXCL10/IP10, CCL5, CXCL1, CXCL3, and CCL20), while METTL14 KD only increased IFNγ, IL-6, TNF, and CXCL3 levels (Li et al., 2021).
m6A epitranscriptomic microarray studies using peripheral blood samples of COVID-19 patients and healthy controls showed altered m6A modification levels in lymphocytes (Meng et al., 2021). Analysis of m6A regulators identified RBM15, WTAP, YTHD3 and IGF2BP1 overexpression and down regulation of METTL16, YTHDF2, YTHDC2 and IGF2BP2 in patients with COVID-19. Hypermethylation of caspase (CASP) 1, CASP5, and tribbles homolog one gene (TRIB1), thymic stromal lymphopoietin (TSLP), DEAD-box helicase 3 X-linked (DDX3X), and interleukin 17 receptor B (IL17RB) mRNAs and subsequent translation promotes programmed cell death and activation of an abnormal inflammatory response in severe patients (Meng et al., 2021).
Pseudouridine (Ψ) in SARS-CoV-2 Vaccines
Epitranscriptomic modifications may be leveraged for therapeutic mRNA vaccine design. Indeed, both Pfizer and Moderna mRNA vaccines utilize Ψ. The Pfizer-BioNTech vaccine (BNT162b2; trade name: Comirnaty; generic name: tozinameran) is an mRNA-based vaccine, whereby synthetic mRNA coding for the SARS-CoV-2 spike protein is encapsulated within a lipid nanoparticle. Upon intramuscular administration, the lipid nanoparticle facilitates delivery of the mRNA payload into local muscle cells or infiltrating immune cells. The mRNA is then translated to produce spike protein, which may be delivered to the plasma membrane as an anchored spike protein or processed and presented by the major histocompatibility complex (MHC), thus generating an immune response. Despite inducing a protective immune response, mRNA vaccine technology has a significant downside to forming undesirable RNA structures. This may cause uncontrolled immune activation, which is detrimental to the host and may also decrease protein translation by inhibiting the ribosome and activating mRNA-degrading ribonucleases, which impairs the efficacy of the mRNA vaccine. The beneficial effects of incorporating Ψ in mRNA vaccines maybe because of the reduction of undesirable RNA structures (Nance and Meier, 2021).
Ψ may reduce the synthesis of antisense RNA, thus improving the translational efficiency of mRNA vaccines1. Ψ may also reduce the formation of double-stranded RNA (dsRNA) (Nance and Meier, 2021). dsRNA is recognized by TLR3 and RIG-1 and initiates the expression of inflammatory genes through NF-κB and IRF activities (Kariko et al., 2004; Durbin et al., 2019; Mauger et al., 2019). Single-stranded poly-uridine has been shown to induce TLR7 activation, an effect that is reduced with Ψ because of the steric incompatibility (Heil et al., 2004; Mauger et al., 2019). N1-methyl-Ψ is the form used in the mRNA vaccines against SARS-CoV-2 (Morais et al., 2021). The methylated-N1 prevents irregular bonding with guanine, uracil, or cytosine and promotes faithful bonding of N1-methylated Ψ with adenine. Consistently, N1-methyl-Ψ inhibited activation of innate immune sensors such as TLR and showed improved protein translation in vivo (Andries et al., 2015).
Altogether, while not explicitly demonstrated robustly in current COVID-19 vaccines yet, these favorable properties of Ψ, which reduce unwanted immunostimulatory potential, justify its adoption in current mRNA vaccines and perhaps in future mRNA-based therapeutics. However, it is essential to note that the location of Ψ within the transcript can influence the net effect on translation. For example, Ψ within the 5’ UTR may block translation initiation, while Ψ within the coding sequence may increase the functional half-life of the mRNA and thus increase translation. Further, once the mRNA of the vaccines is degraded, it is unknown how the modified nucleosides may affect cellular physiology. Characterization of the breakdown and potential salvaging of modified nucleotides from vaccine mRNA is an area that requires further study and attention.
Exploiting base modifications in mRNA vaccines is critical to improving efficacy. This is exemplified by the Curevac COVID-19 mRNA vaccines (CVnCoV), which utilized unmodified mRNA, yielded only 48% efficacy against disease severity (Kremsner et al., 2021). And while these suboptimal results, compared to Pfizer and Moderna vaccines with ∼90% efficacy, maybe due to the lower dose used, it is still likely that using N1-methyl- Ψ mRNA would have increased effectiveness (Dolgin, 2021; Morais et al., 2021).
Conclusion and future Directions
In summary, RNA modifications and the enzymatic machinery are emerging as essential players in regulating SARS-CoV-2 mediated infection and subsequent COVID-19 disease progression. m6A modifications are differentially expressed in host mRNAs in COVID-19 patients, and some essential RNA modifications such as m6A, Ψ, and 2′-O-methylation in SARS-CoV-2 RNA have been identified that permits the virus to interfere with the host’s natural immune response and thereby create a haven for viral replication specifically in the lungs in infected individuals. In this regard, studies in diverse viral infections showed that m6A and 2′-O-methylation control antiviral type I IFN (Kong et al., 2019; Winkler et al., 2019) and interferon gene (ISG) responses (Zheng et al., 2017; Williams et al., 2020).
Future studies at the tissue level are needed to identify the actual writer, reader, and eraser enzyme (s) and PUS synthases in the course of SARS-CoV-2 infection in organs other than the lung. Studies are also needed to determine if pseudouridylation and 2′-O-methylation of SARS-CoV-2 RNA allow the virus to escape host immune response.
There are a plethora of vaccines available to prevent infection (Tatsi et al., 2021) and clinical treatment that includes potential antiviral drug candidates, monoclonal antibodies, and glucocorticoids (Gavriatopoulou et al., 2021) to mitigate the severity of COVID-19. The therapeutic modulation of RNA modification enzymes that are altered or dysregulated during SARS-CoV-2 infection could be a new strategy for attenuating SARS-CoV-2 replication. The strategy that can be employed is to use inhibitors or activators of METTL3, METTL14, YTHDF2, ALKBH5, and FTO to control SARS-CoV-2 replication by modulating the m6A methylome. In this context, selective inhibition of METTL3 catalytic activity using STM2457 inhibits HCoV-OC43 and SARS-CoV-2 infection (Burgess et al., 2021). Since YTHDFs are involved in orchestrating innate antiviral immunity, we can consider them as potential targets of viral countermeasures. Since the ligand-binding site in SARS-CoV-2 nsp16/10 can accommodate small molecules outside of the catalytic pocket, besides SAM- and RNA cap-binding pockets, it is prudent to screen for inhibitors as antiviral therapies to treat SARS-CoV-2 infections.
Author Contributions
Conceptualization, JR and PD; writing: original draft preparation and figures, AI and PD; review and editing, JR and PD. All authors have read and agreed to the content of the manuscript.
Funding
This work was supported by a TNPRC base grant (P51 OD01110459) and startup funds from Tulane University.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Andries, O., Mc Cafferty, S., De Smedt, S. C., Weiss, R., Sanders, N. N., and Kitada, T. (2015). N1-methylpseudouridine-incorporated mRNA Outperforms Pseudouridine-Incorporated mRNA by Providing Enhanced Protein Expression and Reduced Immunogenicity in Mammalian Cell Lines and Mice. J. Controlled Release 217, 337–344. doi:10.1016/j.jconrel.2015.08.051
Banerjee, A. K., Blanco, M. R., Bruce, E. A., Honson, D. D., Chen, L. M., Chow, A., et al. (2020). SARS-CoV-2 Disrupts Splicing, Translation, and Protein Trafficking to Suppress Host Defenses. Cell 183 (5), 1325–1339. doi:10.1016/j.cell.2020.10.004
Bell, J. L., Wächter, K., Mühleck, B., Pazaitis, N., Köhn, M., Lederer, M., et al. (2013). Insulin-like Growth Factor 2 mRNA-Binding Proteins (IGF2BPs): post-transcriptional Drivers of Cancer Progression? Cell. Mol. Life Sci. 70 (15), 2657–2675. doi:10.1007/s00018-012-1186-z
Bestle, D., Heindl, M. R., Limburg, H., Van Lam van, T., Pilgram, O., Moulton, H., et al. (2020). TMPRSS2 and Furin Are Both Essential for Proteolytic Activation of SARS-CoV-2 in Human Airway Cells. Life Sci. Alliance 3 (9), e202000786. doi:10.26508/lsa.202000786
Borchardt, E. K., Martinez, N. M., and Gilbert, W. V. (2020). Regulation and Function of RNA Pseudouridylation in Human Cells. Annu. Rev. Genet. 54, 309–336. doi:10.1146/annurev-genet-112618-043830
Bouvet, M., Debarnot, C., Imbert, I., Selisko, B., Snijder, E. J., Canard, B., et al. (2010). In Vitro Reconstitution Of SARS-Coronavirus Mrna Cap Methylation. PLoS Pathog. 6 (4), e1000863. doi:10.1371/journal.ppat.1000863
Bouvet, M., Imbert, I., Subissi, L., Gluais, L., Canard, B., and Decroly, E. (2012). RNA 327;-end Mismatch Excision By The Severe Acute Respiratory Syndrome Coronavirus Nonstructural Protein nsp10/nsp14 Exoribonuclease Complex. Proc. Natl. Acad. Sci. USA 109 (24), 9372–9377. doi:10.1073/pnas.1201130109
Bujnicki, J. M., Feder, M., Radlinska, M., and Blumenthal, R. M. (2002). Structure Prediction and Phylogenetic Analysis of a Functionally Diverse Family of Proteins Homologous to the MT-A70 Subunit of the Human mRNA:m6A Methyltransferase. J. Mol. Evol. 55 (4), 431–444. doi:10.1007/s00239-002-2339-8
Burgess, H. M., Depledge, D. P., Thompson, L., Srinivas, K. P., Grande, R. C., Vink, E. I., et al. (2021). Targeting the m6A RNA Modification Pathway Blocks SARS-CoV-2 and HCoV-OC43 Replication. Genes Dev. 35 (13-14), 1005–1019. doi:10.1101/gad.348320.121
Campos, J. H. C., Maricato, J. T., Braconi, C. T., Antoneli, F., Janini, L. M. R., and Briones, M. R. S. (2021). Direct RNA Sequencing Reveals SARS-CoV-2 m6A Sites and Possible Differential DRACH Motif Methylation Among Variants. Viruses 13 (11), 2108. doi:10.3390/v13112108
Carlile, T. M., Martinez, N. M., Schaening, C., Su, A., Bell, T. A., Zinshteyn, B., et al. (2019). mRNA Structure Determines Modification by Pseudouridine Synthase 1. Nat. Chem. Biol. 15 (10), 966–974. doi:10.1038/s41589-019-0353-z
Carlile, T. M., Rojas-Duran, M. F., Zinshteyn, B., Shin, H., Bartoli, K. M., and Gilbert, W. V. (2014). Pseudouridine Profiling Reveals Regulated mRNA Pseudouridylation in Yeast and Human Cells. Nature 515, 143–146. doi:10.1038/nature13802
Chen, Y., Su, C., Ke, M., Jin, X., Xu, L., Zhang, Z., et al. (2011). Biochemical and Structural Insights into the Mechanisms of SARS Coronavirus RNA Ribose 2′-O-Methylation by Nsp16/nsp10 Protein Complex. Plos Pathog. 7 (10), e1002294. doi:10.1371/journal.ppat.1002294
Daniloski, K., Jordan, T. X, Wessels, H. H., Hoagland, D. A., Kasela, S., Legut, M., et al. (2021). Identification of Required Host Factors for SARS-CoV-2 Infection in Human Cells. Cell 184 (1), 92–105.e16. doi:10.1016/j.cell.2020.10.030
Dang, W., Xie, Y., Cao, P., Xin, S., Wang, J., Li, S., et al. (2019). N6-Methyladenosine and Viral Infection. Front. Microbiol. 10, 417. doi:10.3389/fmicb.2019.00417
Datta, P. K., Liu, F., Fischer, T., Rappaport, J., and Qin, X. (2020). SARS-CoV-2 Pandemic and Research Gaps: Understanding SARS-CoV-2 Interaction with the ACE2 Receptor and Implications for Therapy. Theranostics 10 (16), 7448–7464. doi:10.7150/thno.48076
Depledge, D. P., and Wilson, A.-C. (2020). Using Direct RNA Nanopore Sequencing To Deconvolute Viral Transcriptomes. Curr. Protoc. Microbiol. 57 (1), e99. doi:10.1002/cpmc.99
Dolgin, E. (2021). Curevac COVID Vaccine Let-Down Spotlights Mrna Design Challenges. Nature 594 (7864), 483. doi:10.1038/d41586-021-01661-0
Dominissini, D., Moshitch-Moshkovitz, S., Schwartz, S., Salmon-Divon, M., Ungar, L., Osenberg, S., et al. (2012). Topology of the Human and Mouse m6A RNA Methylomes Revealed by m6A-Seq. Nature 485 (7397), 201–206. doi:10.1038/nature11112
Durbin, A. F., Wang, C., Marcotrigiano, J., and Gehrke, L. (2016). RNAs Containing Modified Nucleotides Fail to Trigger RIG-I Conformational Changes for Innate Immune Signaling. mBio 7 (5), e00833–16. doi:10.1128/mBio.00833-16
Durbin, J. K., Miller, D. K., Niekamp, J., and Khisamutdinov, E. F. (2019). Modulating Immune Response with Nucleic Acid Nanoparticles. Molecules 24 (20), 3740. doi:10.3390/molecules24203740
Fleming, A. M., Mathewson, N. J., Howpay Manage, S. A., and Burrows, C. J. (2021). Nanopore Dwell Time Analysis Permits Sequencing and Conformational Assignment of Pseudouridine in SARS-CoV-2. ACS Cent. Sci. 7 (10), 1707–1717. doi:10.1021/acscentsci.1c00788
Furuse, Y. (2021). RNA Modifications in Genomic RNA of Influenza A Virus and the Relationship between RNA Modifications and Viral Infection. Ijms 22 (17), 9127. doi:10.3390/ijms22179127
Gavriatopoulou, M., Ntanasis-Stathopoulos, I., Korompoki, E., Fotiou, D., Migkou, M., Tzanninis, I.-G., et al. (2021). Emerging Treatment Strategies for COVID-19 Infection. Clin. Exp. Med. 21 (2), 167–179. doi:10.1007/s10238-020-00671-y
Gordon, D. E., Jang, G. M., Bouhaddou, M., Xu, J., Obernier, K., White, K. M., et al. (2020). A SARS-CoV-2 Protein Interaction Map Reveals Targets for Drug Repurposing. Nature 583 (7816), 459–468. doi:10.1038/s41586-020-2286-9
Gu, Y., Cao, J., Zhang, X., Gao, H., Wang, Y., Wang, J., et al. (2022). Receptome Profiling Identifies KREMEN1 and ASGR1 as Alternative Functional Receptors of SARS-CoV-2. Cell Res 32 (1), 24–37. doi:10.1038/s41422-021-00595-6
Hafner, M., Landthaler, M., Burger, L., Khorshid, M., Hausser, J., Berninger, P., et al. (2010). Transcriptome-wide Identification of RNA-Binding Protein and microRNA Target Sites by PAR-CLIP. Cell 141 (1), 129–141. doi:10.1016/j.cell.2010.03.009
Hashimoto, R., Sakamoto, A., Deguchi, S., Yi, R., Sano, E., Hotta, A., et al. (2021). Dual Inhibition Of TMPRSS2 And Cathepsin Bprevents SARS-Cov-2 Infection In IPS Cells. Mol. Ther. Nucleic Acids. 26, 1107–1114. doi:10.1016/j.omtn.2021.10.016
Heil, F., Hemmi, H., Hochrein, H., Ampenberger, F., Kirschning, C., Akira, S., et al. (2004). Species-specific Recognition of Single-Stranded RNA via Toll-like Receptor 7 and 8. Science 303 (5663), 1526–1529. doi:10.1126/science.1093620
Hoffmann, M., Kleine-Weber, H., Schroeder, S., Krüger, N., Herrler, T., Erichsen, S., et al. (2020). SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell 181 (2), 271–280. e8. doi:10.1016/j.cell.2020.02.052
Huang, H., Weng, H., Sun, W., Qin, X., Shi, H., Wu, H., et al. (2018). Recognition of RNA N6-Methyladenosine by IGF2BP Proteins Enhances mRNA Stability and Translation. Nat. Cel Biol 20, 285–295. doi:10.1038/s41556-018-0045-z
Imam, H., Kim, G.-W., and Siddiqui, A. (2020). Epitranscriptomic(N6-methyladenosine) Modification of Viral RNA and Virus-Host Interactions. Front. Cel. Infect. Microbiol. 10, 584283. doi:10.3389/fcimb.2020.584283
Ivanov, K. A., and Ziebuhr, J. (2004). Human Coronavirus 229E Nonstructural Protein 13: Characterization of Duplex-Unwinding, Nucleoside Triphosphatase, and RNA 5′-Triphosphatase Activities. J. Virol. 78 (14), 7833–7838. doi:10.1128/JVI.78.14.7833-7838.2004
Jia, G., Fu, Y., Zhao, X., Dai, Q., Zheng, G., Yang, Y., et al. (2011). N6-methyladenosine in Nuclear RNA Is a Major Substrate of the Obesity-Associated FTO Erratum in. Nat. Chem. Biol. 7 (12), 885–887. doi:10.1038/nchembio.687
Jora, M., Lobue, P. A., Ross, R. L., Williams, B., and Addepalli, B. (2019). Detection of Ribonucleoside Modifications by Liquid Chromatography Coupled with Mass Spectrometry. Biochim. Biophys. Acta (Bba) - Gene Regul. Mech. 1862 (3), 280–290. doi:10.1016/j.bbagrm.2018.10.012
Kariko, K., Ni, H., Capodici, J., Lamphier, M., and Weissman, D. (2004). mRNA Is an Endogenous Ligand for Toll-like Receptor 3. J. Biol. Chem. 279 (13), 12542–12550. doi:10.1074/jbc.M310175200
Kariko, K., Buckstein, M., Ni, H., and Weissman, D. (2005). Suppression of RNA Recognition by Toll-like Receptors: the Impact of Nucleoside Modification and the Evolutionary Origin of RNA. Immunity 23 (2), 165–175. doi:10.1016/j.immuni.2005.06.008
Kenderdine, T., Nemati, R., Baker, A., Palmer, M., Ujma, J., FitzGibbon, M., et al. (2020). High-resolution Ion Mobility Spectrometry-Mass Spectrometry of Isomeric/isobaric Ribonucleotide Variants. J. Mass. Spectrom. 55 (2), e4465. doi:10.1002/jms.4465
Kim, D., Lee, J. Y., Yang, J. S., Kim, J. W., Kim, V. N., and Chang, H. (2020). The Architecture of SARS-CoV-2 Transcriptome. Cell 181 (4), 914–921.e10. doi:10.1016/j.cell.2020.04.011
Knuckles, P., Lence, T., Haussmann, I. U., Jacob, D., Kreim, N., Carl, S. H., et al. (2018). Zc3h13/Flacc Is Required for Adenosine Methylation by Bridging the mRNA-Binding Factor Rbm15/Spenito to the m6A Machinery Component Wtap/Fl(2)d. Genes Dev. 32 (5-6), 415–429. doi:10.1101/gad.309146.117
Kong, W., Rivera-Serrano, E. E., Neidleman, J. A., and Zhu, J. (2019). HIV-1 Replication Benefits from the RNA Epitranscriptomic Code. J. Mol. Biol. 431 (24), 5032–5038. doi:10.1016/j.jmb.2019.09.021
Kremsner, P. G., Ahuad Guerrero, R. A., Arana-Arri, E., Aroca Martinez, G. J., Bonten, M., Chandler, R., et al. (2021). Efficacy and Safety of the CVnCoV SARS-CoV-2 mRNA Vaccine Candidate in Ten Countries in Europe and Latin America (HERALD): a Randomised, Observer-Blinded, Placebo-Controlled, Phase 2b/3 Trial. Lancet Infect. Dis. S1473-3099 (21), 00677. doi:10.1016/S1473-3099(21)00677-0
Leger, A., Amaral, P. P., Pandolfini, L., Capitanchik, C., Capraro, F., Miano, V., et al. (2021). RNA Modifications Detection by Comparative Nanopore Direct RNA Sequencing. Nat. Commun. 12 (1), 7198. doi:10.1038/s41467-021-27393-3
Li, N., Hui, H., Bray, B., Gonzalez, G. M., Zeller, M., Anderson, K. G., et al. (2021). METTL3 Regulates Viral m6A RNA Modification and Host Cell Innate Immune Responses during SARS-CoV-2 Infection. Cel Rep. 35 (6), 109091. doi:10.1016/j.celrep.2021.109091
Li, S., and Mason, C. E. (2014). The Pivotal Regulatory Landscape of RNA Modifications. Annu. Rev. Genom. Hum. Genet. 15, 127–150. doi:10.1146/annurev-genom-090413-025405
Linder, B., Grozhik, A. V., Olarerin-George, A. O., Meydan, C., Mason, C. E., and Jaffrey, S. R. (2015). Single-nucleotide-resolution Mapping of m6A and m6Am throughout the Transcriptome. Nat. Methods 12 (8), 767–772. doi:10.1038/nmeth.3453
Little, N. A., Hastie, N. D., and Davies, R. C. (2000). Identification of WTAP, a Novel Wilms' Tumour 1-associating Protein. Hum. Mol. Genet. 9 (15), 2231–2239. doi:10.1093/oxfordjournals.hmg.a018914
Liu, J. e., Xu, Y.-P., Li, K., Ye, Q., Zhou, H.-Y., Sun, H., et al. (2021). The m6A Methylome of SARS-CoV-2 in Host Cells. Cel Res 31 (4), 404–414. doi:10.1038/s41422-020-00465-7
Liu, J., Yue, Y., Han, D., Wang, X., Fu, Y., Zhang, L., et al. (2014). A METTL3-METTL14 Complex Mediates Mammalian Nuclear RNA N6-Adenosine Methylation. Nat. Chem. Biol. 10 (2), 93–95. doi:10.1038/nchembio.1432
Loman, N. J., Quick, J., and Simpson, J. T. (2015). A Complete Bacterial Genome Assembled De Novo Using Only Nanopore Sequencing Data. Nat. Methods 12 (8), 733–735. doi:10.1038/nmeth.3444
Lu, M., Zhang, Z., Xue, M., Zhao, B. S., Harder, O., Li, A., et al. (2020). N6-methyladenosine Modification Enables Viral RNA to Escape Recognition by RNA Sensor RIG-I. Nat. Microbiol. 5 (4), 584–598. doi:10.1038/s41564-019-0653-9
Mauger, D. M., Cabral, B. J., Presnyak, V., Su, S. V., Reid, D. W., Goodman, B., et al. (2019). mRNA Structure Regulates Protein Expression through Changes in Functional Half-Life. Proc. Natl. Acad. Sci. USA 116 (48), 24075–24083. doi:10.1073/pnas.1908052116
McIntyre, W., Netzband, R., Bonenfant, G., Biegel, J. M., Miller, C., Fuchs, G., et al. (2018). Positive-sense RNA Viruses Reveal the Complexity and Dynamics of the Cellular and Viral Epitranscriptomes during Infection. Nucleic Acids Res. 46, 5776–5791. doi:10.1093/nar/gky029
Mellott, D. M., Tseng, C.-T., Drelich, A., Fajtová, P., Chenna, B. C., Kostomiris, D. H., et al. (2021). A Clinical-Stage Cysteine Protease Inhibitor Blocks SARS-CoV-2 Infection of Human and Monkey Cells. ACS Chem. Biol. 16 (4), 642–650. doi:10.1021/acschembio.0c00875
Meng, Y., Zhang, Q., Wang, K., Zhang, X., Yang, R., Bi, K., et al. (2021). RBM15-mediated N6-methyladenosine Modification affects COVID-19 severity By Regulating The Expression Of Multitarget Genes. Cell Death Dis. 12 (8), 732. doi:10.1038/s41419-021-04012-z
Meyer, K. D., Saletore, Y., Zumbo, P., Elemento, O., Mason, C. E., and Jaffrey, S. R. (2012). Comprehensive Analysis of mRNA Methylation Reveals Enrichment in 3′ UTRs and Near Stop Codons. Cell 149 (7), 1635–1646. doi:10.1016/j.cell.2012.05.003
Miller, W. A., and Koev, G. (2000). Synthesis of Subgenomic RNAs by Positive-Strand RNA Viruses. Virology 273 (1), 1–8. doi:10.1006/viro.2000.0421
Morais, P., Adachi, H., and Yu, Y.-T. (2021). The Critical Contribution of Pseudouridine to mRNA COVID-19 Vaccines. Front. Cel Dev. Biol. 9, 789427. doi:10.3389/fcell.2021.789427
Nance, K. D., and Meier, J. L. (2021). Modifications in an Emergency: The Role of N1-Methylpseudouridine in COVID-19 Vaccines. ACS Cent. Sci. 7 (5), 748–756. doi:10.1021/acscentsci.1c00197
Netzband, R., and Pager, C. T. (2020). Epitranscriptomic marks: Emerging Modulators of RNA Virus Gene Expression. WIREs RNA 11 (3), e1576. doi:10.1002/wrna.1576
Nielsen, J., Christiansen, J., Lykke-Andersen, J., Johnsen, A. H., Wewer, U. M., and Nielsen, F. C. A. (1999). A Family Of Insulin-Like Growth Factor II Mrna-Binding Proteins Represses Translation In Late Development. Mol. Cell Biol. 19 (2), 1262–1270. doi:10.1128/MCB.19.2.1262
Patil, D. P., Chen, C.-K., Pickering, B. F., Chow, A., Jackson, C., Guttman, M., et al. (2016). m6A RNA Methylation Promotes XIST-Mediated Transcriptional repressionA RNA Methylation Promotes XIST-Mediated Transcriptional Repression. Nature 537 (7620), 369–373. doi:10.1038/nature193427
Ping, X.-L., Sun, B.-F., Wang, L., Xiao, W., Yang, X., Wang, W.-J., et al. (2014). Mammalian WTAP Is a Regulatory Subunit of the RNA N6-Methyladenosine Methyltransferase. Cel Res 24 (2), 177–189. doi:10.1038/cr.2014.3
Roundtree, I. A., Evans, M. E., Pan, T., and He, C. (2017). Dynamic RNA Modifications in Gene Expression Regulation. Cell 169 (7), 1187–1200. doi:10.1016/j.cell.2017.05.045
Safra, M., Nir, R., Farouq, D., Vainberg Slutskin, I., and Schwartz, S. (2017). TRUB1 Is the Predominant Pseudouridine Synthase Acting on Mammalian mRNA via a Predictable and Conserved Code. Genome Res. 27 (3), 393–406. doi:10.1101/gr.207613.116
Saletore, Y., Meyer, K., Korlach, J., Vilfan, I. D., Jaffrey, S., and Mason, C. E. (2012). The Birth of the Epitranscriptome: Deciphering the Function of RNA Modifications. Genome Biol. 13, 175. doi:10.1186/gb-2012-13-10-175
Sarkar, A., Gasperi, W., Begley, U., Nevins, S., Huber, S. M., Dedon, P. C., et al. (2021). Detecting the Epitranscriptome. WIREs RNA 12 (6), e1663. doi:10.1002/wrna.1663
Schmidt, N., Lareau, C. A., Keshishian, H., Ganskih, S., Schneider, C., Hennig, T., et al. (2021). The SARS-CoV-2 RNA-Protein Interactome in Infected Human Cells. Nat. Microbiol. 6 (3), 339–353. doi:10.1038/s41564-020-00846-z
Schwartz, S., Mumbach, M. R., Jovanovic, M., Wang, T., Maciag, K., Bushkin, G. G., et al. (2014). Perturbation of m6A Writers Reveals Two Distinct Classes of mRNA Methylation at Internal and 5′ Sites. Cel Rep. 8 (1), 284–296. doi:10.1016/j.celrep.2014.05.048
Shi, H., Chai, P., Jia, R., and Fan, X. (2020). Novel Insight into the Regulatory Roles of Diverse RNA Modifications: Re-defining the Bridge between Transcription and Translation. Mol. Cancer 19 (1), 78. doi:10.1186/s12943-020-01194-6
Shi, H., Wang, X., Lu, Z., Zhao, B. S., Ma, H., Hsu, P. J., et al. (2017). YTHDF3 Facilitates Translation and Decay of N6-Methyladenosine-Modified RNA. Cel Res 27 (3), 315–328. doi:10.1038/cr.2017.15
Sola, I., Almazán, F., Zúñiga, S., Enjuanes, L., Wewer, U. M., and Nielsen, F. C. A. (2015). Continuous and Discontinuous RNA Synthesis in Coronaviruses. Annu. Rev. Virol. 2 (1), 265–288. doi:10.1146/annurev-virology-100114-055218
Tatsi, E.-B., Filippatos, F., and Michos, A. (2021). SARS-CoV-2 Variants and Effectiveness of Vaccines: a Review of Current Evidence. Epidemiol. Infect. 149, e237. doi:10.1017/S0950268821002430
Tsai, K., and Cullen, B. R. (2020). Epigenetic and Epitranscriptomic Regulation of Viral Replication. Nat. Rev. Microbiol. 18 (10), 559–570. doi:10.1038/s41579-020-0382-3
Viswanathan, T., Arya, S., Chan, S.-H., Qi, S., Dai, N., Misra, A., et al. (2020). Structural Basis of RNA Cap Modification by SARS-CoV-2. Nat. Commun. 11 (1), 3718. doi:10.1038/s41467-020-17496-8
Wang, K., Chen, W., Zhang, Z., Deng, Y., Lian, J.-Q., Du, P., et al. (2020). CD147-spike Protein Is a Novel Route for SARS-CoV-2 Infection to Host Cells. Sig Transduct Target. Ther. 5 (1), 283. doi:10.1038/s41392-020-00426-x
Wang, S., Qiu, Z., Hou, Y., Deng, X., Xu, W., Zheng, T., et al. (2021). AXL Is a Candidate Receptor for SARS-CoV-2 that Promotes Infection of Pulmonary and Bronchial Epithelial Cells. Cel Res 31 (2), 126–140. doi:10.1038/s41422-020-00460-y
Wang, X., Lu, Z., Gomez, A., Hon, G. C., Yue, Y., Han, D., et al. (2014). N6-methyladenosine-dependent Regulation of Messenger RNA Stability. Nature 505 (7481), 117–120. doi:10.1038/nature12730
Wang, X., Zhao, B. S., Roundtree, I. A., Lu, Z., Han, D., Ma, H., et al. (2015). N6-methyladenosine Modulates Messenger RNA Translation Efficiency. Cell 161 (6), 1388–1399. doi:10.1016/j.cell.2015.05.014
Wang, Y., Li, Y., Toth, J. I., Petroski, M. D., Zhang, Z., and Zhao, J. C. (2014). N6-methyladenosine Modification Destabilizes Developmental Regulators in Embryonic Stem Cells. Nat. Cel Biol 16 (2), 191–198. doi:10.1038/ncb2902
Wen, J., Lv, R., Ma, H., Shen, H., He, C., Wang, J., et al. (2018). Zc3h13 Regulates Nuclear RNA m6A Methylation and Mouse Embryonic Stem Cell Self-Renewal. Mol Cell. 69 (6), 1028–1038.e6. doi:10.1016/j.molcel.2018.02.015
Wiener, D., and Schwartz, S. (2021). The Epitranscriptome beyond m6A. Nat. Rev. Genet. 22 (2), 119–131. doi:10.1038/s41576-020-00295-8
Williams, G. D., Gokhale, N. S., Snider, D. L., and Horner, S. M. (2020). The mRNA Cap 2′- O -Methyltransferase CMTR1 Regulates the Expression of Certain Interferon-Stimulated Genes. mSphere 5 (3), e00202–20. doi:10.1128/mSphere.00202-20
Winkler, R., Gillis, E., Lasman, L., Safra, M., Geula, S., Soyris, C., et al. (2019). m6A Modification Controls the Innate Immune Response to Infection by Targeting Type I Interferons. Nat. Immunol. 20, 173–182. doi:10.1038/s41590-018-0275-z
Wojtas, M. N., Pandey, R. R., Mendel, M., Homolka, D., Sachidanandam, R., and Pillai, R. S. (2017). Regulation of m6A Transcripts by the 3ʹ→5ʹ RNA Helicase YTHDC2 Is Essential for a Successful Meiotic Program in the Mammalian Germline. Mol. Cel 68 (2), 374–387. e12. doi:10.1016/j.molcel.2017.09.021
Wu, F., Zhao, S., Yu, B., Chen, Y.-M., Wang, W., Song, Z.-G., et al. (2020). A New Coronavirus Associated with Human Respiratory Disease in China. Nature 579 (7798), 265–269. doi:10.1038/s41586-020-2008-3
Wu, H.-Y., Ke, T.-Y., Liao, W.-Y., and Chang, N.-Y. (2013). Regulation of Coronaviral Poly(A) Tail Length during Infection. PLoS One 8 (7), e70548. doi:10.1371/journal.pone.0070548
Xiao, W., Adhikari, S., Dahal, U., Chen, Y.-S., Hao, Y.-J., Sun, B.-F., et al. (2016). Nuclear M 6 A Reader YTHDC1 Regulates mRNA Splicing. Mol. Cel 61 (4), 507–519. doi:10.1016/j.molcel.2016.01.012
Yang, S. L., DeFalco, L., Anderson, D. E., Zhang, Y., Aw, S. Y., Lim, X. N., et al. (2021). Comprehensive Mapping of SARS-CoV-2 Interactions In Vivo Reveals Functional Virus-Host Interactions. Nat. Commun. 12 (1), 5113. doi:10.1038/s41467-021-25357-1
Zhang, X., Hao, H., Ma, L., Zhang, Y., Hu, X., Chen, Z., et al. (2021). Methyltransferase-like 3 Modulates Severe Acute Respiratory Syndrome Coronavirus-2 RNA N6-Methyladenosine Modification and Replication. mBio 12 (4), e0106721. doi:10.1128/mBio.01067-21
Zhang, Y., Liu, T., Meyer, C. A., Eeckhoute, J., Johnson, D. S., Bernstein, B. E., et al. (2008). Model-based Analysis of ChIP-Seq (MACS). Genome Biol. 9 (9), R137. doi:10.1186/gb-2008-9-9-r137
Zhang, Z., Theler, D., Kaminska, K. H., Hiller, M., de la Grange, P., Pudimat, R., et al. (2010). The YTH Domain Is a Novel RNA Binding Domain. J. Biol. Chem. 285 (19), 14701–14710. doi:10.1074/jbc.M110.104711
Zheng, G., Dahl, J. A., Niu, Y., Fedorcsak, P., Huang, C.-M., Li, C. J., et al. (2013). ALKBH5 Is a Mammalian RNA Demethylase that Impacts RNA Metabolism and Mouse Fertility. Mol. Cel 49 (1), 18–29. doi:10.1016/j.molcel.2012.10.015
Zheng, Q., Hou, J., Zhou, Y., Li, Z., and Cao, X. (2017). The RNA Helicase DDX46 Inhibits Innate Immunity by Entrapping m6A-Demethylated Antiviral Transcripts in the Nucleus. Nat. Immunol. 18, 1094–1103. doi:10.1038/ni.3830
Keywords: epitranscriptome, SARS-CoV-2, RNA modifications, lung, COVID-19
Citation: Izadpanah A, Rappaport J and Datta PK (2022) Epitranscriptomics of SARS-CoV-2 Infection. Front. Cell Dev. Biol. 10:849298. doi: 10.3389/fcell.2022.849298
Received: 05 January 2022; Accepted: 01 March 2022;
Published: 08 April 2022.
Edited by:
Ana Soares, University of Aveiro, PortugalReviewed by:
Ricardo Soto-Rifo, University of Chile, ChileJia Sheng, University at Albany, United States
Copyright © 2022 Izadpanah, Rappaport and Datta. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Prasun K. Datta, cGRhdHRhQHR1bGFuZS5lZHU=