Novel viral splicing events and open reading frames revealed by long-read direct RNA sequencing of adenovirus transcripts.

Adenovirus is a common human pathogen that relies on host cell processes for transcription and processing of viral RNA and protein production. Although adenoviral promoters, splice junctions, and polyadenylation sites have been characterized using low-throughput biochemical techniques or short read cDNA-based sequencing, these technologies do not fully capture the complexity of the adenoviral transcriptome. By combining Illumina short-read and nanopore long-read direct RNA sequencing approaches, we mapped transcription start sites and RNA cleavage and polyadenylation sites across the adenovirus genome. In addition to confirming the known canonical viral early and late RNA cassettes, our analysis of splice junctions within long RNA reads revealed an additional 35 novel viral transcripts that meet stringent criteria for expression. These RNAs include fourteen new splice junctions which lead to expression of canonical open reading frames (ORFs), six novel ORF-containing transcripts, and 15 transcripts encoding for messages that could alter protein functions through truncation or fusion of canonical ORFs. In addition, we detect RNAs that bypass canonical cleavage sites and generate potential chimeric proteins by linking distinct gene transcription units. Among these chimeric proteins we detected an evolutionarily conserved protein containing the N-terminus of E4orf6 fused to the downstream DBP/E2A ORF. Loss of this novel protein, E4orf6/DBP, was associated with aberrant viral replication center morphology and poor viral spread. Our work highlights how long-read sequencing technologies combined with mass spectrometry can reveal further complexity within viral transcriptomes and resulting proteomes.

DRUMMER-rapid detection of RNA modifications through comparative nanopore sequencing.

MOTIVATION: The chemical modification of ribonucleotides regulates the structure, stability and interactions of RNAs. Profiling of these modifications using short-read (Illumina) sequencing techniques provides high sensitivity but low-to-medium resolution i.e. modifications cannot be assigned to specific transcript isoforms in regions of sequence overlap. An alternative strategy uses current fluctuations in nanopore-based long read direct RNA sequencing (DRS) to infer the location and identity of nucleotides that differ between two experimental conditions. While highly sensitive, these signal-level analyses require high-quality transcriptome annotations and thus are best suited to the study of model organisms. By contrast, the detection of RNA modifications in microbial organisms which typically have no or low-quality annotations requires an alternative strategy. Here, we demonstrate that signal fluctuations directly influence error rates during base-calling and thus provides an alternative approach for identifying modified nucleotides. RESULTS: DRUMMER (Detection of Ribonucleic acid Modifications Manifested in Error Rates) (i) utilizes a range of statistical tests and background noise correction to identify modified nucleotides with high confidence, (ii) operates with similar sensitivity to signal-level analysis approaches and (iii) correlates very well with orthogonal approaches. Using well-characterized DRS datasets supported by independent meRIP-Seq and miCLIP-Seq datasets we demonstrate that DRUMMER operates with high sensitivity and specificity. AVAILABILITY AND IMPLEMENTATION: DRUMMER is written in Python 3 and is available as open source in the GitHub repository: https://github.com/DepledgeLab/DRUMMER. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Adenovirus prevents dsRNA formation by promoting efficient splicing of viral RNA.

Eukaryotic cells recognize intracellular pathogens through pattern recognition receptors, including sensors of aberrant nucleic acid structures. Sensors of double-stranded RNA (dsRNA) are known to detect replication intermediates of RNA viruses. It has long been suggested that annealing of mRNA from symmetrical transcription of both top and bottom strands of DNA virus genomes can produce dsRNA during infection. Supporting this hypothesis, nearly all DNA viruses encode inhibitors of dsRNA-recognition pathways. However, direct evidence that DNA viruses produce dsRNA is lacking. Contrary to dogma, we show that the nuclear-replicating DNA virus adenovirus (AdV) does not produce detectable levels of dsRNA during infection. In contrast, abundant dsRNA is detected within the nucleus of cells infected with AdV mutants defective for viral RNA processing. In the presence of nuclear dsRNA, the cytoplasmic dsRNA sensor PKR is relocalized and activated within the nucleus. Accumulation of viral dsRNA occurs in the late phase of infection, when unspliced viral transcripts form intron/exon base pairs between top and bottom strand transcripts. We propose that DNA viruses actively limit dsRNA formation by promoting efficient splicing and mRNA processing, thus avoiding detection and restriction by host innate immune sensors of pathogenic nucleic acids.

Adenovirus Remodeling of the Host Proteome and Host Factors Associated with Viral Genomes.

Viral infections are associated with extensive remodeling of the cellular proteome. Viruses encode gene products that manipulate host proteins to redirect cellular processes or subvert antiviral immune responses. Adenovirus (AdV) encodes proteins from the early E4 region which are necessary for productive infection. Some cellular antiviral proteins are known to be targeted by AdV E4 gene products, resulting in their degradation or mislocalization. However, the full repertoire of host proteome changes induced by viral E4 proteins has not been defined. To identify cellular proteins and processes manipulated by viral products, we developed a global, unbiased proteomics approach to analyze changes to the host proteome during infection with adenovirus serotype 5 (Ad5) virus. We used whole-cell proteomics to measure total protein abundances in the proteome during Ad5 infection. Since host antiviral proteins can antagonize viral infection by associating with viral genomes and inhibiting essential viral processes, we used Isolation of Proteins on Nascent DNA (iPOND) proteomics to identify proteins associated with viral genomes during infection with wild-type Ad5 or an E4 mutant virus. By integrating these proteomics data sets, we identified cellular factors that are degraded in an E4-dependent manner or are associated with the viral genome in the absence of E4 proteins. We further show that some identified proteins exert inhibitory effects on Ad5 infection. Our systems-level analysis reveals cellular processes that are manipulated during Ad5 infection and points to host factors counteracted by early viral proteins as they remodel the host proteome to promote efficient infection. IMPORTANCE Viral infections induce myriad changes to the host cell proteome. As viruses harness cellular processes and counteract host defenses, they impact abundance, post-translational modifications, interactions, or localization of cellular proteins. Elucidating the dynamic changes to the cellular proteome during viral replication is integral to understanding how virus-host interactions influence the outcome of infection. Adenovirus encodes early gene products from the E4 genomic region that are known to alter host response pathways and promote replication, but the full extent of proteome modifications they mediate is not known. We used an integrated proteomics approach to quantitate protein abundance and protein associations with viral DNA during virus infection. Systems-level analysis identifies cellular proteins and processes impacted in an E4-dependent manner, suggesting ways that adenovirus counteracts potentially inhibitory host defenses. This study provides a global view of adenovirus-mediated proteome remodeling, which can serve as a model to investigate virus-host interactions of DNA viruses.

Comparative proteomics identifies Schlafen 5 (SLFN5) as a herpes simplex virus restriction factor that suppresses viral transcription.

Intrinsic antiviral host factors confer cellular defence by limiting virus replication and are often counteracted by viral countermeasures. We reasoned that host factors that inhibit viral gene expression could be identified by determining proteins bound to viral DNA (vDNA) in the absence of key viral antagonists. Herpes simplex virus 1 (HSV-1) expresses E3 ubiquitin-protein ligase ICP0 (ICP0), which functions as an E3 ubiquitin ligase required to promote infection. Cellular substrates of ICP0 have been discovered as host barriers to infection but the mechanisms for inhibition of viral gene expression are not fully understood. To identify restriction factors antagonized by ICP0, we compared proteomes associated with vDNA during HSV-1 infection with wild-type virus and a mutant lacking functional ICP0 (ΔICP0). We identified the cellular protein Schlafen family member 5 (SLFN5) as an ICP0 target that binds vDNA during HSV-1 ΔICP0 infection. We demonstrated that ICP0 mediates ubiquitination of SLFN5, which leads to its proteasomal degradation. In the absence of ICP0, SLFN5 binds vDNA to repress HSV-1 transcription by limiting accessibility of RNA polymerase II to viral promoters. These results highlight how comparative proteomics of proteins associated with viral genomes can identify host restriction factors and reveal that viral countermeasures can overcome SLFN antiviral activity.

Direct RNA sequencing reveals m6A modifications on adenovirus RNA are necessary for efficient splicing.

Adenovirus is a nuclear replicating DNA virus reliant on host RNA processing machinery. Processing and metabolism of cellular RNAs can be regulated by METTL3, which catalyzes the addition of N6-methyladenosine (m6A) to mRNAs. While m6A-modified adenoviral RNAs have been previously detected, the location and function of this mark within the infectious cycle is unknown. Since the complex adenovirus transcriptome includes overlapping spliced units that would impede accurate m6A mapping using short-read sequencing, here we profile m6A within the adenovirus transcriptome using a combination of meRIP-seq and direct RNA long-read sequencing to yield both nucleotide and transcript-resolved m6A detection. Although both early and late viral transcripts contain m6A, depletion of m6A writer METTL3 specifically impacts viral late transcripts by reducing their splicing efficiency. These data showcase a new technique for m6A discovery within individual transcripts at nucleotide resolution, and highlight the role of m6A in regulating splicing of a viral pathogen.

Adenovirus-mediated ubiquitination alters protein-RNA binding and aids viral RNA processing.

Viruses promote infection by hijacking the ubiquitin machinery of the host to counteract or redirect cellular processes. Adenovirus encodes two early proteins, E1B55K and E4orf6, that together co-opt a cellular ubiquitin ligase complex to overcome host defences and promote virus production. Adenovirus mutants lacking E1B55K or E4orf6 display defects in viral RNA processing and protein production, but previously identified substrates of the redirected ligase do not explain these phenotypes. Here, we used a quantitative proteomics approach to identify substrates of E1B55K/E4orf6-mediated ubiquitination that facilitate RNA processing. While all currently known cellular substrates of E1B55K and E4orf6 are degraded by the proteasome, we uncovered RNA-binding proteins as high-confidence substrates that are not decreased in overall abundance. We focused on two RNA-binding proteins, RALY and hnRNP-C, which we confirm are ubiquitinated without degradation. Knockdown of RALY and hnRNP-C increased levels of viral RNA splicing, protein abundance and progeny production during infection with E1B55K-deleted virus. Furthermore, infection with E1B55K-deleted virus resulted in an increased interaction of hnRNP-C with viral RNA and attenuation of viral RNA processing. These data suggest that viral-mediated ubiquitination of RALY and hnRNP-C relieves a restriction on viral RNA processing and reveal an unexpected role for non-degradative ubiquitination in the manipulation of cellular processes during virus infection.

Identification of Host Biomarkers of Epstein-Barr Virus Latency IIb and Latency III.

Deciphering the molecular pathogenesis of virally induced cancers is challenging due, in part, to the heterogeneity of both viral gene expression and host gene expression. Epstein-Barr virus (EBV) is a ubiquitous herpesvirus prevalent in B-cell lymphomas of immune-suppressed individuals. EBV infection of primary human B cells leads to their immortalization into lymphoblastoid cell lines (LCLs), serving as a model of these lymphomas. In previous studies, reports from our laboratory have described a temporal model for immortalization with an initial phase characterized by expression of Epstein-Barr nuclear antigens (EBNAs), high levels of c-Myc activity, and hyperproliferation in the absence of the latent membrane proteins (LMPs), called latency IIb. This is followed by the long-term outgrowth of LCLs expressing the EBNAs along with the LMPs, particularly NFκB-activating LMP1, defining latency III. However, LCLs express a broad distribution of LMP1 such that a subset of these cells express LMP1 at levels similar to those seen in latency IIb, making it difficult to distinguish these two latency states. In this study, we performed mRNA sequencing (mRNA-Seq) on early EBV-infected latency IIb cells and latency III LCLs sorted by NFκB activity. We found that latency IIb transcriptomes clustered independently from latency III independently of NFκB. We identified and validated mRNAs defining these latency states. Indeed, we were able to distinguish latency IIb cells from LCLs expressing low levels of LMP1 using multiplex RNA-fluorescence in situ hybridization (RNA-FISH) targeting EBV EBNA2 or LMP1 and human CCR7 or MGST1 This report defines latency IIb as a bona fide latency state independent from latency III and identifies biomarkers for understanding EBV-associated tumor heterogeneity.IMPORTANCE EBV is a ubiquitous pathogen, with >95% of adults harboring a life-long latent infection in memory B cells. In immunocompromised individuals, latent EBV infection can result in lymphoma. The established expression profile of these lymphomas is latency III, which includes expression of all latency genes. However, single-cell analysis of EBV latent gene expression in these lymphomas suggests heterogeneity where most cells express the transcription factor, EBNA2, and only a fraction of the cells express membrane protein LMP1. Our work describes an early phase after infection where the EBNAs are expressed without LMP1, called latency IIb. However, LMP1 levels within latency III vary widely, making these states hard to discriminate. This may have important implications for therapeutic responses. It is crucial to distinguish these states to understand the molecular pathogenesis of these lymphomas. Ultimately, better tools to understand the heterogeneity of these cancers will support more-efficacious therapies in the future.

c-Myc Represses Transcription of Epstein-Barr Virus Latent Membrane Protein 1 Early after Primary B Cell Infection.

Recent evidence has shown that the Epstein-Barr virus (EBV) oncogene LMP1 is not expressed at high levels early after EBV infection of primary B cells, despite its being essential for the long-term outgrowth of immortalized lymphoblastoid cell lines (LCLs). In this study, we found that expression of LMP1 increased 50-fold between 7 days postinfection and the LCL state. Metabolic labeling of nascent transcribed mRNA indicated that this was primarily a transcription-mediated event. EBNA2, the key viral transcription factor regulating LMP1, and CTCF, an important chromatin insulator, were recruited to the LMP1 locus similarly early and late after infection. However, the activating histone H3K9Ac mark was enriched at the LMP1 promoter in LCLs relative to that in infected B cells early after infection. We found that high c-Myc activity in EBV-infected lymphoma cells as well as overexpression of c-Myc in an LCL model system repressed LMP1 transcription. Finally, we found that chemical inhibition of c-Myc both in LCLs and early after primary B cell infection increased LMP1 expression. These data support a model in which high levels of endogenous c-Myc activity induced early after primary B cell infection directly repress LMP1 transcription.IMPORTANCE EBV is a highly successful pathogen that latently infects more than 90% of adults worldwide and is also causally associated with a number of B cell malignancies. During the latent life cycle, EBV expresses a set of viral oncoproteins and noncoding RNAs with the potential to promote cancer. Critical among these is the viral latent membrane protein LMP1. Prior work suggests that LMP1 is essential for EBV to immortalize B cells, but our recent work indicates that LMP1 is not produced at high levels during the first few weeks after infection. Here we show that transcription of the LMP1 gene can be negatively regulated by a host transcription factor, c-Myc. Ultimately, understanding the regulation of EBV oncogenes will allow us to better treat cancers that rely on these viral products for survival.

Take your PIKK: tumour viruses and DNA damage response pathways.

Viruses regulate cellular processes to facilitate viral replication. Manipulation of nuclear proteins and pathways by nuclear replicating viruses often causes cellular genome instability that contributes to transformation. The cellular DNA damage response (DDR) safeguards the host to maintain genome integrity, but DNA tumour viruses can manipulate the DDR to promote viral propagation. In this review, we describe the interactions of DNA tumour viruses with the phosphatidylinositol 3-kinase-like protein kinase (PIKK) pathways, which are central regulatory arms of the DDR. We review how signalling through the ataxia telangiectasia mutated (ATM), ataxia telangiectasia and Rad3 related (ATR), and DNA-dependent protein kinases (DNA-PK) influences viral life cycles, and how their manipulation by viral proteins may contribute to tumour formation.This article is part of the themed issue 'Human oncogenic viruses'.

The Epstein-Barr virus miR-BHRF1 microRNAs regulate viral gene expression in cis.

The Epstein-Barr virus (EBV) miR-BHRF1 microRNA (miRNA) cluster has been shown to facilitate B-cell transformation and promote the rapid growth of the resultant lymphoblastoid cell lines (LCLs). However, we find that expression of physiological levels of the miR-BHRF1 miRNAs in LCLs transformed with a miR-BHRF1 null mutant (∆123) fails to increase their growth rate. We demonstrate that the pri-miR-BHRF1-2 and 1-3 stem-loops are present in the 3'UTR of transcripts encoding EBNA-LP and that excision of pre-miR-BHRF1-2 and 1-3 by Drosha destabilizes these mRNAs and reduces expression of the encoded protein. Therefore, mutational inactivation of pri-miR-BHRF1-2 and 1-3 in the ∆123 mutant upregulates the expression of not only EBNA-LP but also EBNA-LP-regulated mRNAs and proteins, including LMP1. We hypothesize that this overexpression causes the reduced transformation capacity of the ∆123 EBV mutant. Thus, in addition to regulating cellular mRNAs in trans, miR-BHRF1-2 and 1-3 also regulate EBNA-LP mRNA expression in cis.

Epstein-Barr virus ensures B cell survival by uniquely modulating apoptosis at early and late times after infection.

Latent Epstein-Barr virus (EBV) infection is causally linked to several human cancers. EBV expresses viral oncogenes that promote cell growth and inhibit the apoptotic response to uncontrolled proliferation. The EBV oncoprotein LMP1 constitutively activates NFκB and is critical for survival of EBV-immortalized B cells. However, during early infection EBV induces rapid B cell proliferation with low levels of LMP1 and little apoptosis. Therefore, we sought to define the mechanism of survival in the absence of LMP1/NFκB early after infection. We used BH3 profiling to query mitochondrial regulation of apoptosis and defined a transition from uninfected B cells (BCL-2) to early-infected (MCL-1/BCL-2) and immortalized cells (BFL-1). This dynamic change in B cell survival mechanisms is unique to virus-infected cells and relies on regulation of MCL-1 mitochondrial localization and BFL-1 transcription by the viral EBNA3A protein. This study defines a new role for EBNA3A in the suppression of apoptosis with implications for EBV lymphomagenesis.

Mitogen-induced B-cell proliferation activates Chk2-dependent G1/S cell cycle arrest.

B-cell activation and proliferation can be induced by a variety of extracellular stimuli. The fate of an activated B cell following mitogen stimulation can be dictated by the strength or duration of the signal, the expression of downstream signaling components necessary to promote proliferation, and the cell intrinsic sensors and regulators of the proliferative program. Previously we have identified the DNA damage response (DDR) signaling pathway as a cell intrinsic sensor that is activated upon latent infection of primary human B cells by Epstein-Barr virus (EBV). Here we have assessed the role of the DDR as a limiting factor in the proliferative response to non-viral B-cell mitogens. We report that TLR9 activation through CpG-rich oligonucleotides induced B-cell hyper-proliferation and an ATM/Chk2 downstream signaling pathway. However, B-cell activation through the CD40 pathway coupled with interleukin-4 (IL-4) promoted proliferation less robustly and only a modest DDR. These two mitogens, but not EBV, modestly induced intrinsic apoptosis that was independent from the DDR. However, all three mitogens triggered a DDR-dependent G1/S phase cell cycle arrest preventing B-cell proliferation. The extent of G1/S arrest, as evidenced by release through Chk2 inhibition, correlated with B-cell proliferation rates. These findings have implications for the regulation of extra-follicular B-cell activation as it may pertain to the development of auto-immune diseases or lymphoma.

Dynamic Epstein-Barr virus gene expression on the path to B-cell transformation.

Epstein-Barr virus (EBV) is an oncogenic human herpesvirus in the γ-herpesvirinae subfamily that contains a 170-180kb double-stranded DNA genome. In vivo, EBV commonly infects B and epithelial cells and persists for the life of the host in a latent state in the memory B-cell compartment of the peripheral blood. EBV can be reactivated from its latent state, leading to increased expression of lytic genes that primarily encode for enzymes necessary to replicate the viral genome and structural components of the virion. Lytic cycle proteins also aid in immune evasion, inhibition of apoptosis, and the modulation of other host responses to infection. In vitro, EBV has the potential to infect primary human B cells and induce cellular proliferation to yield effectively immortalized lymphoblastoid cell lines, or LCLs. EBV immortalization of B cells in vitro serves as a model system for studying EBV-mediated lymphomagenesis. While much is known about the steady-state viral gene expression within EBV-immortalized LCLs and other EBV-positive cell lines, relatively little is known about the early events after primary B-cell infection. It was previously thought that upon latent infection, EBV only expressed the well-characterized latency-associated transcripts found in LCLs. However, recent work has characterized the early, but transient, expression of lytic genes necessary for efficient transformation and delayed responses in the known latency genes. This chapter summarizes these recent findings that show how dynamic and controlled expression of multiple EBV genes can control the activation of B cells, entry into the cell cycle, the inhibition of apoptosis, and innate and adaptive immune responses.

Analysis of Epstein-Barr virus-regulated host gene expression changes through primary B-cell outgrowth reveals delayed kinetics of latent membrane protein 1-mediated NF-κB activation.

Epstein-Barr virus (EBV) is an oncogenic human herpesvirus that dramatically reorganizes host gene expression to immortalize primary B cells. In this study, we analyzed EBV-regulated host gene expression changes following primary B-cell infection, both during initial proliferation and through transformation into lymphoblastoid cell lines (LCLs). While most EBV-regulated mRNAs were changed during the transition from resting, uninfected B cells through initial B-cell proliferation, a substantial number of mRNAs changed uniquely from early proliferation through LCL outgrowth. We identified constitutively and dynamically EBV-regulated biological processes, protein classes, and targets of specific transcription factors. Early after infection, genes associated with proliferation, stress responses, and the p53 pathway were highly enriched. However, the transition from early to long-term outgrowth was characterized by genes involved in the inhibition of apoptosis, the actin cytoskeleton, and NF-κB activity. It was previously thought that the major viral protein responsible for NF-κB activation, latent membrane protein 1 (LMP1), is expressed within 2 days after infection. Our data indicate that while this is true, LCL-level LMP1 expression and NF-κB activity are not evident until 3 weeks after primary B-cell infection. Furthermore, heterologous NF-κB activation during the first week after infection increased the transformation efficiency, while early NF-κB inhibition had no effect on transformation. Rather, inhibition of NF-κB was not toxic to EBV-infected cells until LMP1 levels and NF-κB activity were high. These data collectively highlight the dynamic nature of EBV-regulated host gene expression and support the notion that early EBV-infected proliferating B cells have a fundamentally distinct growth and survival phenotype from that of LCLs.