2. Cellular m6A regulating proteins and their functions
The addition of a methyl group to the
N6 position of adenosine is catalyzed by a multi-subunit methyltransferase complex, also referred to as the m
6A writer complex. The catalytic core is formed by a heterodimer of methyltransferase-like 3 (METTL3) and methyltransferase-like 14 (METTL14), with Wilms tumor 1-associating protein (WTAP) serving as a scaffold protein that stabilizes the methyltransferase complex [
7,
8]. Additional components of the methyltransferase complex that may direct the complex to sites of m
6A modification include RNA binding motifs protein 15 (RBM15) and KIAA1429, also known as vir-like m
6A methyltransferase associated (VIRMA) [
9,
10]. However, the function of each protein associated with the methyltransferase complex is not fully understood.
The writer complex exhibits preference for modification of adenosine in the consensus sequence 5’-DRACH-3' (D=A/G/U, R=A/G, H= U/A/C) in mammalian mRNAs [
11]. Transcriptome-wide m
6A mapping techniques have revealed that m
6A is most common in 3’ UTRs and long exons of mRNA [
11]. Although DRACH motifs occur frequently, only ~10% of them are m
6A methylated in cells, suggesting that there are context-dependent signals for the methylation of specific DRACH motifs [
11].
m
6A modification is reversible through the action of demethylases, or erasers. These enzymes include alkB homolog 5 (ALKBH5) and fat mass and obesity-associated protein (FTO), and both enzymes can remove m
6A modification [
12,
13]. However, since the discovery of these enzymes, it has been shown that FTO exhibits a strong preference for demethylating
N6, 2’-O-methyladenosine (m
6Am) residues next to the 5’ cap, rather than m
6A at internal sites of mRNA [
14]. Therefore, while expression of the individual eraser enzymes can be manipulated to affect overall changes in m
6A levels, the role of each endogenous eraser in the context of virus infection in cells remain unclear.
The functional consequence of m
6A modification is determined by the binding of m
6A-specific RNA binding proteins, or readers. Readers that are most often studied in the context of m
6A function are the YTH domain-containing family of proteins, YTHDC1 and YTHDF1-3. YTHDC1 regulates the nuclear export of mature mRNA [
3]. The cytoplasmic YTHDF1-3 proteins were initially determined to be functionally distinct. YTHDF1 promotes translation, while YTHDF2 mediates degradation of m
6A-modified RNAs [
1]. YTHDF3 acts in synergy with YTHDF1 and YTHDF2 to enhance translation or mediate RNA decay [
15]. However, more recent studies suggest that YTHDF1-3 are functionally redundant and mediate RNA degradation [
16]. The discrepancies regarding YTHDF protein functions have been reviewed elsewhere [
17,
18,
19]. While YTHD family readers are often the focus of studies designed to elucidate the role of m
6A in viral replication, many other potential reader proteins have also been described (reviewed in [
17]). These less studied m
6A readers also deserve attention in investigating the mechanisms of m
6A in regulating cellular gene expression and virus replication.
A major challenge in determining the function of m
6A in regulating virus replication is the ability to specifically remove m
6A on viral RNA, especially transcripts, while leaving the modification of cellular RNA unperturbed. Almost all studies to date have utilized overexpression or knockdown/knockout of writers, readers, and erasers to infer the function of viral RNA m
6A during infection. However, these approaches will also lead to changes in the m
6A modification of host cell RNA, which may have indirect effects on virus replication. Such indirect effects must always be carefully considered during the interpretation of these studies. Additional approaches, including mutations of m
6A site of viral genomes, are often used in confirming specific effects [
20,
21,
22,
23,
24].
3. Mapping m6A modification sites on HIV-1 RNA
Several studies have reported the location of m
6A in HIV-1 RNA from infected cells or purified virions (
Figure 1A,B). The techniques, virus strains, and cell types used among studies are summarized in
Table 1. In the first report, methylated RNA immunoprecipitation sequencing (MeRIP-seq) was used to determine m
6A modified regions of HIV-1 RNA from infected MT4 CD4+ T cells [
25]. meRIP-seq utilizes fragmented RNA for m
6A-specific immunoprecipitation followed by high-throughput sequencing and has a resolution of ~100-200 nt [
11]. This analysis identified 14 regions predicted to contain m
6A modifications [
25] (
Figure 1A). In contrast to subsequent reports, this data analysis did not include the 5’ and 3’ UTR sequences, which led to the absence of input or m
6A signal mapping to these regions. Peaks mapped primarily to coding sequences and exhibited overlap with several splicing regulatory sequences. Unique to this study was the identification of a peak located in
env and
rev response element (RRE), which was the focus of subsequent functional analysis. Methylation of predicted sites in the RRE was confirmed by primer extension assay in infected cells [
25]. However, none of the subsequent studies reported by other groups found m
6A to be present in the RRE [
26,
27].
The next report identified sites of m
6A modification in HIV-1 RNA from CEM-SS CD4+ T cells infected with single-cycle HIV-1 [
26]. This study employed photo-crosslinking-assisted m
6A sequencing (PA-m
6A-seq), which modifies the meRIP protocol through the addition of 4-thiouridine (4SU) labeling of viral RNA followed by UV crosslinking of an m
6A-specific antibody to the labeled RNA prior to sequencing. This approach improves the resolution of sites of m
6A modification to ~20 nt due to the introduction of U to C mutations at the sites of crosslink between antibody and 4SU [
32]. In this analysis, all predicted sites of methylation mapped to the 3’ UTR ~15% of the genome (
Figure 1A). The PA-m
6A-seq results were compared to photoactivatable ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) performed to identify the binding sites for YTHDF reader proteins on the HIV-1 RNA. Comparison of reader binding sites in HIV-1 RNA from HEK293T cells transfected with an HIV-1 proviral plasmid and m
6A peaks from infected CEM-SS cells reveal four regions in common, located in
env/rev,
nef, and the 3’ UTR. The authors also sought to assess the topology of reader binding sites in different HIV-1 strains. PAR-CLIP data to assess reader binding sites in RNA of HIV-1
BaL and HIV-1
JR-CSF show that these primary isolates have the four binding sites in common with HIV-1
NL4-3, while also possessing additional reader binding peaks. This illustrates the importance of considering cell-type and HIV-1 strain-specific differences in m
6A modification of viral RNA. This group also reported the location of m
6A modifications in HIV-1 genomic RNA (gRNA) purified from virions compared to cell-associated RNA from CEM-SS cells. PA-m
6A-seq revealed major peaks in
env/rev,
nef, and 3’ UTR of gRNA, consistent with the results of m
6A mapping to cell-associated viral RNA (
Figure 1B) [
28]. Interestingly, additional peaks were detected in several coding sequences of virion RNA including
pol,
env, and
vpr.
Another study using meRIP-seq confirmed the presence of major m
6A peaks in
env/rev,
nef, and the 3’ UTR in HIV-1 RNA from infected Jurkat and primary CD4+ T cells (
Figure 1A) [
27]. This analysis also identified additional smaller but prominent peaks in the 5’ UTR and
gag. Cell-associated HIV-1 RNA collected from HEK293T cells transfected with an HIV-1 proviral plasmid showed similar sites of m
6A modification as viral RNA from infected CD4+ T cells, with an additional major peak in
tat. These results suggest that in addition to cell type differences, it is possible that patterns of HIV-1 m
6A modification differ in virus producer cells and target cells. Validation of m
6A modification at the predicted sites was performed by crosslinking-immunoprecipitation sequencing (CLIP-seq) to identify the binding sites of YTHDF1-3 proteins on HIV-1 RNA isolated from infected HeLa/CD4 cells. CLIP-seq analysis revealed reader binding sites in
env/
rev,
nef, and the 5’ UTR, providing further support for the presence of m
6A at these locations. Of note, reader binding was detected in several coding regions of the HIV-1 genome, sites that were not detected by meRIP-seq. This could be due to different cell types used for meRIP-seq and CLIP-seq or due to more sensitive and robust detection using a crosslinking approach.
Mapping sites of m
6A modification at a single timepoint does not provide information about how modification states may change over the course of infection. One group investigated whether temporal changes in m
6A modification of viral RNA occur at different times during a single-cycle HIV-1 infection [
29]. meRIP-seq was performed using cell-associated RNA from SupT1 cells infected with single-cycle HIV-1 for 12, 24, or 36 hrs post-infection (hpi). Each of these time points is earlier than the previously reported analyses, all performed at 96 hpi. m
6A peak mapping at these earlier time points again showed enrichment of m
6A in the 3’ end of the RNA, with peaks also detected in
pol,
pol/vif, and
vpu (
Figure 1A). One notable change in m
6A mapping over time revealed a large peak overlapping the RNA packaging signal ψ at 36 hpi but not earlier. RNA isolated from virions released from the same cells at 36 hpi retains m
6A modification in the
env/nef and 3’ UTR regions, as well at the peak overlapping ψ, but is largely lacking prominent peaks in the remaining coding regions (
Figure 1B). These results suggest that m
6A modification patterns differ between packaged gRNA and translated viral transcripts. The peak width defined for these analyses is quite broad and ranges from ~200-500 nt, while meRIP-seq involves fragmentation of RNA to 100-200 nt in length. This low resolution prevents meaningful prediction of specific DRACH motifs that may be modified.
Another report comparing m
6A modification of infected-cell HIV-1 RNA compared to RNA purified from viral particles also showed that the cell-associated RNA has a different pattern of methylation than that of packaged gRNA [
30]. RNA was prepared from HEK293T cells transfected with HIV-1 proviral plasmid or from supernatant virus. Consistent with all other reports, both viral RNA preparations exhibited peaks overlapping
env/rev,
nef, and the 3’ UTR (
Figure 1A, B). In contrast to the results reported by Cristinelli et al. [
29], a peak was observed in the 5’ UTR of intracellular RNA that was absent from virion RNA. The authors conclude that m
6A in the 5’ UTR regulates genome packaging and focus on this phenotype for subsequent functional studies.
Techniques such as meRIP-seq and PA-m
6A-seq reveal the population-level modification state of fragmented RNA, and therefore cannot provide information on the presence of m
6A in individual viral genomes or transcripts. Therefore, these data do not provide insight into whether there is differential m
6A deposition among specific splice isoforms. The most recent analysis of m
6A sites in HIV-1 RNA was performed using a modification of direct RNA sequencing (DRS) that resulted in a remarkable improvement of full-length genome reads from 0.01% to 34.9% [
31]. This approach allowed for full-length, single-molecule analysis at single-nucleotide resolution for both viral gRNA and splice isoforms. Additional stringency was provided by using three independent bioinformatics tools for the identification of m
6A signals. Seven common sites were identified across all three analyses. Interestingly, two of these sites, both in
pol, were not located in canonical DRACH motif sequences (
Figure 1A). Four sites in the 1.2 kb of the 3’ HIV-1 genome coincide with the common peaks identified in
env/
rev and the 3’ UTR in each of the previous studies [
26,
27,
28,
29]. These m
6A are located at 8079, 8110, 8975, and 8989 of HIV-1
NL4-3 [
31].
Almost all available data describing the location of m
6A in HIV-1 RNA is derived from cell lines infected
in vitro. One study demonstrated that the pattern of viral RNA modification is similar in a CD4+ T cell line and primary CD4+ T cells infected
ex vivo, suggesting that these results might be relevant to HIV-1 infected individuals [
27]. The only information regarding
in vivo HIV-1 infection comes from a transgenic rat model of HIV-1 infection [
33]. These animals harbor a
gag/
pol deleted provirus that expresses viral genes in microglia and astrocytes [
34]. meRIP-seq of RNA isolated from the hippocampus of these animals revealed m
6A in the 5’ UTR and
nef and 3’UTR, consistent with cells infected
in vitro and
ex vivo (
Figure 1A) [
27,
33].
Differences in m
6A site prediction among studies can arise for several reasons, including biological, technical, and computational differences in the experimental approach. Such variables include HIV-1 strain, cell type, virus purification method, sequencing approach, data processing, and bioinformatic tools. In addition to the issue of low resolution provided by meRIP-seq, concerns have been raised about reproducibility of the method, even when using the same experimental conditions [
35]. Future studies would benefit from the use of sequencing approaches that allow for the identification of m
6A at single-nucleotide resolution.
Despite some inconsistencies among studies, m
6A located in
env (2 sites) and the 3’ UTR (2 sites) have been reproducibly identified across several data sets. In addition, adenosine residues in viral DRACH motifs that are predicted to contain m
6A are highly conserved in the HIV-1 genome that is otherwise highly susceptible to mutation [
25,
26,
31,
36]. In addition, only a very small fraction of the nearly 250 DRACH motifs in the HIV-1 full-length RNA [
31] appear to contain m
6A, suggesting that modification of these sites is very highly specific. These data suggest that these m
6A modifications are especially important in facilitating efficient virus replication.
4. HIV-1 infection modulates the cellular RNA m6A profile
Epitranscriptomic modification of mRNA represents an important layer of gene expression regulation in mammalian cells. While much focus has been placed on determining the m6A profile of HIV-1 RNA, another area of interest is how viral infection influences changes in m6A modification of host cell transcripts. This can provide mechanistic insights into how HIV-1 induces changes in host cell gene expression, either directly through the action of viral proteins or because of the intrinsic antiviral response to infection.
Infection of various cell types including differentiated monocytic, microglial and CD4+ T cell lines as well as primary CD4+ T cells with HIV-1 induces an upregulation of m
6A levels in total cellular RNA, through mechanisms that remain unclear [
25,
29,
37,
38,
39]. Upregulation of m
6A in CD4+ T cells and microglial cells is not due to changes in the expression of m
6A writers or erasers [
37]. It is possible that methyltransferase complex activity is being regulated at the level of complex formation, stability, or enzymatic activity.
In CD4+ T cells, m
6A upregulation is also observed in cells treated with the reverse transcriptase inhibitor nevirapine, heat inactivated virus, and gp120, the surface unit of the HIV-1 envelope protein [
37]. These data demonstrate that m
6A upregulation is independent of virus replication. m
6A upregulation occurs after treatment of cells with recombinant gp120 proteins from CXCR4- and CCR5-tropic HIV-1 and is blocked by occlusion of the primary receptor CD4 with anti-CD4 antibodies [
37]. Whether HIV-1 co-receptor binding to gp120 is also required for m
6A upregulation of cellular RNA remains unknown. Analysis of m
6A levels in polyadenylated and non-coding RNA reveal upregulation of m
6A in both, including ribosomal RNA (rRNA) [
37]. The 18S and 28S rRNA subunits both contain a single m
6A that is modified by the methyltransferases METTL5 and ZCCHC4, respectively [
40,
41]. Ribosomes will form normally in the absence of METTL5 or ZCCHC4, but the overall translational efficiency is reduced [
40,
42]. This raises the interesting possibility that HIV-1 may enhance its protein translation through m
6A modification of the ribosome; however, this possible mechanism has not yet been explored.
High throughput sequencing approaches such as meRIP-seq allow for the identification of changes in the relative abundance of m
6A-modified host cell mRNAs. The first transcriptome-wide meRIP-seq defining the topology of m
6A in human cells showed an enrichment of m
6A near stop codons and in coding sequences, particularly long exons [
11]. Mapping of m
6A in cellular transcripts from HIV-1-infected CD4+ T cell lines compared to uninfected controls reveals no significant changes in the distribution of m
6A deposition to 5’ UTR, exons, introns, or 3’ UTR, with the majority of m
6A sites found in exons and the 3’ UTR as expected [
25,
27,
29].
Changes in relative abundance of m
6A-modified transcripts alone does not indicate differential methylation of the RNA, because these changes could simply be a result of differential expression of m
6A-modified transcripts. Therefore, true differential m
6A methylation of HIV-1 infected cells can only be determined by measuring modification level and transcript abundance and comparing those data to the same transcript in uninfected control cells. Applying this strategy, over 3,615 transcripts were found to be hypermethylated upon HIV-1 infection compared to hypermethylation of only 777 transcripts [
29]. These data support an overall upregulation of m
6A modification of cellular transcripts in SupT1 cells upon HIV-1 infection. Comparing m
6A methylation patterns at 12, 24, and 36 hrs post-infection, only 2% of differentially methylated transcripts were common among all three time points [
29]. This shows that m
6A modification patterns of cellular RNA are dynamic and temporally regulated during HIV-1 infection.
A new analysis combining both RNA-seq and meRIP-seq data showed that changes in transcript methylation occurred disproportionately in coding sequences compared to the 3’ UTR in a microglial cell line during a single-cycle HIV-1 infection [
38]. Over 2,000 transcripts were reported to be hypo- or hyper-methylated in infected cells, suggesting that HIV-1 infection induces changes in gene expression at the level of m
6A modification. Pathway analysis of these transcripts show significant enrichment of genes involved in signal transduction, including the Ras and MAPK signaling pathways [
38]. Further studies will be required to determine the functional effect of differential methylation on these signaling pathways during HIV-1 infection of microglia.
5. HIV-1 m6A suppresses the induction of type I interferon (IFN-I) in macrophages
Internal 2’-O-methylation of HIV-1 RNA serves as an immune evasion strategy during HIV-1 infection of monocyte-derived macrophages and dendritic cells [
43]. To test whether m
6A also suppresses IFN-I in response to HIV-1 infection, viruses were produced in the presence of eraser (ALKBH5 or FTO) overexpression or knockout [
44]. Differentiated monocytic U937 cells (macrophage-like cells) were then infected with HIV-1 containing low, normal, or high levels of m
6A in the incoming RNA genome. The important advantage of manipulating m
6A levels in producer cells rather than target cells is that it avoids the perturbation of cellular m
6A levels so that observed phenotypes could be attributed specifically to m
6A in the incoming HIV-1 gRNA. Using complementary approaches for modulating m
6A, the authors demonstrate that the level of m
6A in the HIV-1 genome is inversely correlated with the expression of IFN-α and IFN-β after transfection or infection of differentiated U937 cells [
44]. The phenotype was observed in primary monocyte-derived macrophages transfected with HIV-1 RNA [
44]. Taken together, the results suggest that m
6A suppresses the expression of IFN-I by HIV-1 RNA. Delivery of RNA to cells by transfection does not recapitulate the pathway by which gRNA is delivered to cells during infection. Therefore, future mechanistic studies should be done using HIV-1 infected cells, particularly macrophages.
Viral RNA is detected in the host cell by protein sensors that recognize pathogen-associated molecular patterns in viral molecules to discriminate self from non-self RNA [
45]. Binding of these sensors to their RNA ligands activates a signaling cascade to induce the expression of IFN-I. To determine which host cell cytosolic RNA sensor may be inhibited by m
6A-modified HIV-1 RNA, U937 knockout (KO) or knockdown (KD) cell lines were generated lacking expression of retinoic acid-induced gene I (RIG-I) or melanoma differentiation-associated gene 5 (MDA5), respectively. These cells were transfected with a 42-mer HIV-1 derived RNA oligonucleotide with or without a single m
6A. The sequence of the oligonucleotide corresponds to a portion of the HIV-1 5’ UTR of HIV-1, where potential m
6A sites were predicted by meRIP-seq studies [
25,
27]. While the expression of IFN-α and IFN-β was lower in both RIG-I KO and MDA5 KD cells compared to control, RIG-I KO cells were no longer able to discriminate between oligonucleotides with or without m
6A [
44]. These data suggest that m
6A may function in evading RIG-I sensing, and this would be consistent with a growing body of research in the field of innate immunity to RNA viruses [
46,
47,
48,
49]. Also, it was shown that the infection of differentiated U937 cells with m
6A deficient HIV-1 induces more phosphorylation of IRF3 and IRF7 as compared to control viral particles. Based on this, it appears that the presence of m
6A on HIV-1 RNA hinders IFN-1 production during the early stages of HIV-1 infection. However, a major limitation of this study is the use of RNA oligonucleotide transfection, especially given uncertainties regarding the exposure of HIV-1 genomes to the cytoplasm during transit to the nuclear pore [
50,
51,
52,
53]. It will be critical to determine whether m
6A in the HIV-1 genome avoids detection by RIG-I in the context of infected cells, and if so, what is the underlying mechanisms. It is possible that other cellular RNA sensors might be involved in sensing m
6A lower HIV-1 RNA.
7. m6A regulates HIV-1 RNA splicing and nuclear export
HIV-1 gene expression is driven by the viral 5’ long terminal repeat (LTR) of integrated provirus in a productively infected cell. The deposition of m
6A on nascent transcripts occurs co-transcriptionally in the nucleus [
57]. It is therefore possible for m
6A to affect RNA splicing and nuclear export. The first study reporting the presence of m
6A in HIV-1 RNA identified m
6A enrichment in the RRE, with two potential methylation sites in hairpin loop IIB [
25]. The RRE is a highly conserved and structured RNA element in HIV-1 RNA that is essential for the nuclear export of partially spliced or unspliced transcripts [
58]. The authors hypothesized that m
6A within the RRE may be required for efficient viral RNA export. Mutation of the predicted methylation site in the hairpin IIB bulge dramatically reduced the levels of
env mRNA in HEK293T cells transfected with HIV-1 proviral DNA plasmid. Subcellular fractionation prior to measurement of
env mRNA indicated a defect in nuclear export. However, Rev binding to these mutated RRE was not measured. A variety of
in vitro assays to assess the structure and Rev binding properties of synthesized RRE IIB with or without m
6A argued against a significant impact of m
6A modification on IIB structure [
59]. Further, a separate study made the same RRE mutations in a different HIV-1 laboratory strain and showed that when
env and
rev are provided
in trans, there is no difference in viral gene expression in either producer or target cells [
36]. The RRE was not identified as a site of m
6A modification or reader binding in any subsequent studies [
26,
27,
28]. Therefore, the function of m
6A in the RRE remains to be confirmed.
The nuclear m
6A reader YTHDC1 is known to be involved in RNA stability, splicing, and nuclear export of m
6A-modified viral and cellular RNA [
3,
60,
61]. Two independent studies assessed HIV-1 RNA splicing and nuclear export under conditions of overexpression or depletion of YTHDC1 during single-round infections [
26,
62]. Both studies found that YTHDC1 is required for not only maintaining overall viral RNA levels, but also for appropriate selection of splice sites. Splicing of HIV-1 transcripts is complex and involves the utilization of multiple splice donor and acceptor sites [
63]. Specifically, YTHDC1 increased the use of more 3’ splice acceptor sites with a concomitant decrease in the use of splice acceptor sites further 5’ in the HIV-1 genome [
56]. Manipulating the expression levels of YTHDC1 had no effect on the cellular NONO mRNA that is not m
6A-modified [
56]. This suggests that the observed changes in HIV-1 RNA upon YTHDC1 silencing or overexpression are due to the presence of m
6A (
Figure 2). Overexpressed and endogenous YTHDC1 bind to spliced and unspliced HIV-1 RNA in infected cells in a METTL3-dependent manner [
62] (
Figure 2). YTHDC1 binding sites were identified adjacent to HIV-1 splice acceptor (SA) sites 3 and 7 [
56]. However, when the predicted m
6A site near SA3 was mutated, there was no effect on the utilization of SA3 [
56]. Therefore, whether YTHDC1 is acting directly on viral m
6A is unclear. Both groups also reported that overexpression or depletion of YTHDC1 had no effect on the nuclear/cytoplasmic accumulation of unspliced, partially spliced, or completely spliced viral RNA [
56,
62]. These results argue against a role for YTHDC1 in HIV-1 viral RNA nuclear export, which is unexpected given the role of YTHDC1 in nuclear export of cellular mRNA [
3].
Several independent studies reported m
6A modification in the 3’ region of HIV-1 RNA overlap between
env and the second exon of
rev [
26,
27,
28,
29,
30,
33]. DRS identified the modified adenosine in a DRACH motif in the
env/rev overlap at position 8079 of HIV-1
NL4-3, downstream of SA7 [
31]. Two other high confidence m
6A are located in the 3’ UTR. These m
6A represent three of the four most prominent sites of m
6A modification, in agreement with other studies, and were mutated for functional studies. Interestingly, mutation of these sites individually had no significant effect on viral protein expression, the proportion of unspliced RNA, p24 release, or virion infectivity [
31]. Only when all three adenosines were mutated was there a significant reduction in virus infectivity. The reduction in HIV-1 replication is likely due to a decrease in the abundance of unspliced RNA, which is necessary to produce infectious virions. The triple mutation also results in modest but significant increases in the proportion of partially and completely spliced viral RNA relative to total RNA [
31]. This suggests functional redundancy of the individual m
6A and that the modifications are involved in the regulation of splicing. This idea is further supported by the observation that completely spliced RNAs are more m
6A-modified at these three predominant sites than unspliced RNAs [
31]. It will be interesting to uncover the mechanisms by which these three m
6A site modulate HIV-1 RNA splicing, whether other m
6A are also functionally redundant, and how a single nucleotide modification can substitute for the function of another.
8. m6A enhances post-integration HIV-1 RNA abundance, stability, and translation
One strategy to indirectly test whether m
6A inhibits or enhances HIV-1 replication is to manipulate the expression of writer complex components or m
6A erasers in target cells. Silencing or overexpression of writers or erasers would presumably only effect viral RNA that is made post-integration; however, this has not been assessed in detail. Several studies that have used this approach agree that in general, m
6A is required for efficient HIV-1 replication as measured by viral RNA, viral protein, and progeny particle release [
25,
27,
30,
56]. One exception to this general phenotype is a study using single-cycle HIV-1 transduction of HeLa cells [
62]. The exact reason for the discrepancy is unknown but could be due to the kinetics of viral gene expression in this experimental system compared to replication competent HIV-1 entry by plasma membrane fusion.
Another area of controversy in the field of HIV-1 and m
6A is the effect of readers on viral replication. One study measured HIV-1 RNA and protein levels in HEK293T producer cells overexpressing each YTHDF1-3 reader individually and found that in all cases, both HIV-1 RNA and protein levels were increased [
26] (
Figure 2). In addition, silencing of YTHDF2 decreases the half-life of HIV-1 RNA and cellular RNAs containing m
6A, yet had no significant effect on cellular RNAs that are not though to be modified [
56]. These results all suggest that m
6A readers enhance HIV-1 replication by stabilizing viral RNA. The overexpressed readers are assumed to be acting on viral RNA, however, enhanced replication can also be due to changes in cellular gene expression that indirectly affect viral gene expression. Consistent results were reported after YTHDF2 overexpression or knockout in CEM T-cells infected with single-cycle HIV-1 [
26,
56]. In these cells there was a reduction in Gag protein levels and p24 levels of released virions, again suggesting that YTHDF2 is required for efficient viral gene expression (
Figure 2). These results are in direct opposition to other reports that show a decrease in viral RNA and protein in infected HeLa, Jurkat, and primary CD4+ T cells [
27,
54]. As discussed above, this was concluded to be a result of inhibition of RT through binding of readers to incoming viral RNA leading ultimately to its degradation [
27,
54,
55]. One possibility for these seemingly opposing results is that m
6A has a negative effect on HIV-1 gRNA but a potential benefit at post-integration phases of infection. The use of different cell types, viruses, and length of infection make direct comparison of data from different groups imprecise.
Two potential functions for viral RNA m
6A are destabilization of RNA or enhanced translation efficiency mediated by m
6A readers [
1,
2]. To assess the functional significance of m
6A modifications in the 3’-UTR of HIV-1 RNA, Kennedy et al. placed the wild-type or m
6A-deficient (A to G mutated) forms of the HIV-1 3’-UTR downstream of a reporter gene [
26]. Their results revealed that substituting A with G significantly reduced both RNA and protein levels of the reporter. These data argue against RNA destabilization or enhanced translation and indicate that m
6A in the 3’ UTR of HIV-1 RNA enhance the abundance of an unrelated RNA (
Figure 2). Recruiting any of the YTHDF readers to the 3’UTR of a reporter, independent of m
6A, also resulted in enhanced expression of the reporter [
26]. Taken together, these results suggest that m
6A modification in the 3’UTR of HIV-1 mRNA enhances viral gene expression at the mRNA and protein levels.
9. m6A inhibits HIV-1 RNA packaging and reduces virion infectivity
Several groups reported the presence of m
6A in the 5’ UTR of HIV-1 RNA [
27,
29,
30]. The 5’ UTR is highly structured and contains several regulatory elements that are indispensable for efficient HIV-1 replication [
64]. One important aspect of HIV-1 replication that is mediated by RNA sequences in the 5’ UTR is the selection of full-length RNA genomes for packaging into progeny virions [
65]. There are two DRACH motifs of particular interest in the 5’UTR. One is in the primer binding site (PBS) required for the initiation of RT and another is near the dimer initiation sequence (DIS) that is involved in a monomer-dimer switch that regulates whether a given full-length HIV-1 RNA is translated (monomer) or selected for packaging (dimer) [
54].
Mutation of either of these two adenosine residues led to a slight increase in p24 levels in HEK293T producer cells, yet reduced infectivity of the release virions when used to infect target cells [
54]. RNA secondary structure prediction
in silico showed an altered structure of the PBS but not the DIS because of the mutations [
54]. However, due to the possibility of the disruption of primer binding for RT or Gag binding during particle assembly as a result of these mutations, the interpretation of these data with respect to m
6A modification can be multifaceted.
Several lines of evidence suggest that m
6A modification impairs packaging of gRNA. meRIP-seq data derived from intracellular viral RNA and virion RNA revealed a prominent peak representing m
6A in the 5’ UTR of intracellular RNA, but not RNA from purified virus particles [
30]. In addition, the m
6A/A level of intracellular viral RNA was significantly higher than that in packaged RNA [
30]. This raises the possibility that m
6A is playing a role in the preferential selection of gRNA molecules for packaging. Overexpression of METTL3/14, which presumably increases m
6A modification of viral RNA, resulted in higher levels of Gag protein and lower levels of RNA packaging, whereas silencing of the methyltransferase complex resulted in the opposite phenotype [
30]. These results are consistent with the observed decrease in Gag protein in other studies where METTL3 is silenced [
25,
27]. Interestingly, the interaction between Gag and full-length viral RNA in the cytoplasm was decreased under conditions of METTL3/14 overexpression, further supporting the idea that increased methylation of viral RNA inhibits is ability to be packaged [
30]. Unexpectedly, deletion of the two adenosine residues that are potential sites of methylation resulted in a significant increase in overall m
6A content of viral RNA and a significant decrease in packaging efficiency of gRNA [
30]. Due to deletion rather than mutation of these adenosines, it is unknown whether the proper secondary structure of the 5’ UTR was affected. Another study also demonstrated that virion RNA contains fewer m
6A modifications that viral transcripts, however this study only reported modification of three adenosines in the 3’ ~1.4 kb of the viral genome. In addition, YTHDF proteins bind to m
6A modified gRNA and form a complex with Gag in an RNA-dependent manner [
27,
54]. It is therefore possible that YTHDF binding prevents efficient binding of Gag and therefore inhibits packaging. Further experiments are needed to investigate this hypothesis.
HIV-1 Gag colocalizes in the nucleus with unspliced viral RNA [
66]. Interestingly, Gag also colocalizes with FTO exclusively in the nucleus and there is an increase in viral RNA m
6A levels in the absence of gag expression that can be rescued when Gag is provided
in trans [
30]. These data suggest that FTO mediates demethylation of viral RNA in a Gag-dependent manner (
Figure 2). However, while ALKBH5 and FTO can both demethylate HIV-1 RNA when each demethylase is overexpressed, it is not known which of these proteins acts on viral RNA under conditions of endogenous expression levels. Emerging evidence supports a predominant role for FTO in the removal of methyl groups from m
6Am and not m
6A in human cells [
67,
68]. Further studies are needed to determine the molecular mechanisms of the removal of m
6A from HIV-1 RNA and how this process is regulated to control translation or packaging of full-length gRNA.