Dengue and Zika RNA-RNA interactomes reveal pro- and anti-viral RNA in human cells.

BACKGROUND: Identifying host factors is key to understanding RNA virus pathogenicity. Besides proteins, RNAs can interact with virus genomes to impact replication. RESULTS: Here, we use proximity ligation sequencing to identify virus-host RNA interactions for four strains of Zika virus (ZIKV) and one strain of dengue virus (DENV-1) in human cells. We find hundreds of coding and non-coding RNAs that bind to DENV and ZIKV viruses. Host RNAs tend to bind to single-stranded regions along the virus genomes according to hybridization energetics. Compared to SARS-CoV-2 interactors, ZIKV-interacting host RNAs tend to be downregulated upon virus infection. Knockdown of several short non-coding RNAs, including miR19a-3p, and 7SK RNA results in a decrease in viral replication, suggesting that they act as virus-permissive factors. In addition, the 3'UTR of DYNLT1 mRNA acts as a virus-restrictive factor by binding to the conserved dumbbell region on DENV and ZIKV 3'UTR to decrease virus replication. We also identify a conserved set of host RNAs that interacts with DENV, ZIKV, and SARS-CoV-2, suggesting that these RNAs are broadly important for RNA virus infection. CONCLUSIONS: This study demonstrates that host RNAs can impact virus replication in permissive and restrictive ways, expanding our understanding of host factors and RNA-based gene regulation during viral pathogenesis.

m6A governs length-dependent enrichment of mRNAs in stress granules.

Stress granules are biomolecular condensates composed of protein and mRNA. One feature of stress granule-enriched mRNAs is that they are often longer than average. Another feature of stress granule-enriched mRNAs is that they often contain multiple N6-methyladenosine (m6A) residues. m6A is bound by the YTHDF proteins, creating mRNA-protein complexes that partition into stress granules in mammalian cells. Here we show that length-dependent enrichment of mRNAs in stress granules is mediated by m6A. Long mRNAs often contain one or more long exons, which are preferential sites of m6A formation. In mammalian cells lacking m6A, long mRNAs no longer show preferential stress granule enrichment. Furthermore, we show that m6A abundance more strongly predicts which short or long mRNAs are enriched in stress granules, rather than length alone. Thus, mRNA length correlates with mRNA enrichment in stress granules owing to the high prevalence of m6A in long mRNAs.

BTG1 mutation yields supercompetitive B cells primed for malignant transformation.

Multicellular life requires altruistic cooperation between cells. The adaptive immune system is a notable exception, wherein germinal center B cells compete vigorously for limiting positive selection signals. Studying primary human lymphomas and developing new mouse models, we found that mutations affecting BTG1 disrupt a critical immune gatekeeper mechanism that strictly limits B cell fitness during antibody affinity maturation. This mechanism converted germinal center B cells into supercompetitors that rapidly outstrip their normal counterparts. This effect was conferred by a small shift in MYC protein induction kinetics but resulted in aggressive invasive lymphomas, which in humans are linked to dire clinical outcomes. Our findings reveal a delicate evolutionary trade-off between natural selection of B cells to provide immunity and potentially dangerous features that recall the more competitive nature of unicellular organisms.

Alternative splicing of METTL3 explains apparently METTL3-independent m6A modifications in mRNA.

N6-methyladenosine (m6A) is a highly prevalent mRNA modification that promotes degradation of transcripts encoding proteins that have roles in cell development, differentiation, and other pathways. METTL3 is the major methyltransferase that catalyzes the formation of m6A in mRNA. As 30% to 80% of m6A can remain in mRNA after METTL3 depletion by CRISPR/Cas9-based methods, other enzymes are thought to catalyze a sizable fraction of m6A. Here, we reexamined the source of m6A in the mRNA transcriptome. We characterized mouse embryonic stem cell lines that continue to have m6A in their mRNA after Mettl3 knockout. We show that these cells express alternatively spliced Mettl3 transcript isoforms that bypass the CRISPR/Cas9 mutations and produce functionally active methyltransferases. We similarly show that other reported METTL3 knockout cell lines express altered METTL3 proteins. We find that gene dependency datasets show that most cell lines fail to proliferate after METTL3 deletion, suggesting that reported METTL3 knockout cell lines express altered METTL3 proteins rather than have full knockout. Finally, we reassessed METTL3's role in synthesizing m6A using an exon 4 deletion of Mettl3 and found that METTL3 is responsible for >95% of m6A in mRNA. Overall, these studies suggest that METTL3 is responsible for the vast majority of m6A in the transcriptome, and that remaining m6A in putative METTL3 knockout cell lines is due to the expression of altered but functional METTL3 isoforms.

The RNA m6A Reader YTHDF2 Maintains Oncogene Expression and Is a Targetable Dependency in Glioblastoma Stem Cells.

Glioblastoma is a universally lethal cancer driven by glioblastoma stem cells (GSC). Here, we interrogated N 6-methyladenosine (m6A) mRNA modifications in GSCs by methyl RNA immunoprecipitation followed by sequencing and transcriptome analysis, finding transcripts marked by m6A often upregulated compared with normal neural stem cells (NSC). Interrogating m6A regulators, GSCs displayed preferential expression, as well as in vitro and in vivo dependency, of the m6A reader YTHDF2, in contrast to NSCs. Although YTHDF2 has been reported to destabilize mRNAs, YTHDF2 stabilized MYC and VEGFA transcripts in GSCs in an m6A-dependent manner. We identified IGFBP3 as a downstream effector of the YTHDF2-MYC axis in GSCs. The IGF1/IGF1R inhibitor linsitinib preferentially targeted YTHDF2-expressing cells, inhibiting GSC viability without affecting NSCs and impairing in vivo glioblastoma growth. Thus, YTHDF2 links RNA epitranscriptomic modifications and GSC growth, laying the foundation for the YTHDF2-MYC-IGFBP3 axis as a specific and novel therapeutic target in glioblastoma. SIGNIFICANCE: Epitranscriptomics promotes cellular heterogeneity in cancer. RNA m6A landscapes of cancer and NSCs identified cell type-specific dependencies and therapeutic vulnerabilities. The m6A reader YTHDF2 stabilized MYC mRNA specifically in cancer stem cells. Given the challenge of targeting MYC, YTHDF2 presents a therapeutic target to perturb MYC signaling in glioblastoma.This article is highlighted in the In This Issue feature, p. 211.

Structure mapping of dengue and Zika viruses reveals functional long-range interactions.

Dengue (DENV) and Zika (ZIKV) viruses are clinically important members of the Flaviviridae family with an 11 kb positive strand RNA genome that folds to enable virus function. Here, we perform structure and interaction mapping on four DENV and ZIKV strains inside virions and in infected cells. Comparative analysis of SHAPE reactivities across serotypes nominates potentially functional regions that are highly structured, conserved, and contain low synonymous mutation rates. Interaction mapping by SPLASH identifies many pair-wise interactions, 40% of which form alternative structures, suggesting extensive structural heterogeneity. Analysis of shared interactions between serotypes reveals a conserved macro-organization whereby interactions can be preserved at physical locations beyond sequence identities. We further observe that longer-range interactions are preferentially disrupted inside cells, and show the importance of new interactions in virus fitness. These findings deepen our understanding of Flavivirus genome organization and serve as a resource for designing therapeutics in targeting RNA viruses.