Viral reprogramming of host transcription initiation.

Viruses are master remodelers of the host cell environment in support of infection and virus production. For example, viruses typically regulate cell gene expression through modulating canonical cell promoter activity. Here, we show that Epstein Barr virus (EBV) replication causes 'de novo' transcription initiation at 29674 new transcription start sites throughout the cell genome. De novo transcription initiation is facilitated in part by the unique properties of the viral pre-initiation complex (vPIC) that binds a TATT[T/A]AA, TATA box-like sequence and activates transcription with minimal support by additional transcription factors. Other de novo promoters are driven by the viral transcription factors, Zta and Rta and are influenced by directional proximity to existing canonical cell promoters, a configuration that fosters transcription through existing promoters and transcriptional interference. These studies reveal a new way that viruses interact with the host transcriptome to inhibit host gene expression and they shed light on primal features driving eukaryotic promoter function.

Transcriptomic analysis identifies B-lymphocyte kinase as a therapeutic target for desmoplastic small round cell tumor cancer stem cell-like cells.

Desmoplastic small round cell tumor (DSRCT) is an aggressive pediatric cancer caused by the EWSR1-WT1 fusion oncoprotein. The tumor is refractory to treatment with a 5-year survival rate of only 15-25%, necessitating the development of novel therapeutics, especially those able to target chemoresistant subpopulations. Novel in vitro cancer stem cell-like (CSC-like) culture conditions increase the expression of stemness markers (SOX2, NANOG) and reduce DSRCT cell line susceptibility to chemotherapy while maintaining the ability of DSRCT cells to form xenografts. To gain insights into this chemoresistant model, RNA-seq was performed to elucidate transcriptional alterations between DSRCT cells grown in CSC-like spheres and normal 2-dimensional adherent state. Commonly upregulated and downregulated genes were identified and utilized in pathway analysis revealing upregulation of pathways related to chromatin assembly and disassembly and downregulation of pathways including cell junction assembly and extracellular matrix organization. Alterations in chromatin assembly suggest a role for epigenetics in the DSRCT CSC-like state, which was further investigated with ATAC-seq, identifying over 10,000 differentially accessible peaks, including 4444 sphere accessible peaks and 6,120 adherent accessible peaks. Accessible regions were associated with higher gene expression, including increased accessibility of the CSC marker SOX2 in CSC-like culture conditions. These analyses were further utilized to identify potential CSC therapeutic targets, leading to the identification of B-lymphocyte kinase (BLK) as a CSC-enriched, EWSR1-WT1-regulated, druggable target. BLK inhibition and knockdown reduced CSC-like properties, including abrogation of tumorsphere formation and stemness marker expression. Importantly, BLK knockdown reduced DSRCT CSC-like cell chemoresistance, making its inhibition a promising target for future combination therapy.

The transcription factor c-Jun inhibits RBM39 to reprogram pre-mRNA splicing during genotoxic stress.

Genotoxic agents, that are used in cancer therapy, elicit the reprogramming of the transcriptome of cancer cells. These changes reflect the cellular response to stress and underlie some of the mechanisms leading to drug resistance. Here, we profiled genome-wide changes in pre-mRNA splicing induced by cisplatin in breast cancer cells. Among the set of cisplatin-induced alternative splicing events we focused on COASY, a gene encoding a mitochondrial enzyme involved in coenzyme A biosynthesis. Treatment with cisplatin induces the production of a short isoform of COASY lacking exons 4 and 5, whose depletion impedes mitochondrial function and decreases sensitivity to cisplatin. We identified RBM39 as a major effector of the cisplatin-induced effect on COASY splicing. RBM39 also controls a genome-wide set of alternative splicing events partially overlapping with the cisplatin-mediated ones. Unexpectedly, inactivation of RBM39 in response to cisplatin involves its interaction with the AP-1 family transcription factor c-Jun that prevents RBM39 binding to pre-mRNA. Our findings therefore uncover a novel cisplatin-induced interaction between a splicing regulator and a transcription factor that has a global impact on alternative splicing and contributes to drug resistance.

Endothelial extracellular vesicles promote tumour growth by tumour-associated macrophage reprogramming.

Tumour-derived extracellular vesicles (EVs) participate in tumour progression by deregulating various physiological processes including angiogenesis and inflammation. Here we report that EVs released by endothelial cells in a mammary tumour environment participate in the recruitment of macrophages within the tumour, leading to an immunomodulatory phenotype permissive for tumour growth. Using RNA-Seq approaches, we identified several microRNAs (miRNAs) found in endothelial EVs sharing common targets involved in the regulation of the immune system. To further study the impact of these miRNAs in a mouse tumour model, we focused on three miRNAs that are conserved between humans and mouse, that is, miR-142-5p, miR-183-5p and miR-222-3p. These miRNAs are released from endothelial cells in a tumour microenvironment and are transferred via EVs to macrophages. In mouse mammary tumour models, treatment with EVs enriched in these miRNAs leads to a polarization of macrophages toward an M2-like phenotype, which in turn promotes tumour growth.

Reversal of splicing infidelity is a pre-activation step in B cell differentiation.

Introduction: B cell activation and differentiation is central to the adaptive immune response. Changes in exon usage can have major impacts on cellular signaling and differentiation but have not been systematically explored in differentiating B cells. Methods: We analyzed exon usage and intron retention in RNA-Seq data from subsets of human B cells at various stages of differentiation, and in an in vitro laboratory model of B cell activation and differentiation (Epstein Barr virus infection). Results: Blood naïve B cells were found to have an unusual splicing profile, with unannotated splicing events in over 30% of expressed genes. Splicing changed substantially upon naïve B cell entry into secondary lymphoid tissue and before activation, involving significant increases in exon commitment and reductions in intron retention. These changes preferentially involved short introns with weak splice sites and were likely mediated by an overall increase in splicing efficiency induced by the lymphoid environment. The majority of transcripts affected by splicing changes showed restoration of encoded conserved protein domains and/or reduced targeting to the nonsense-mediated decay pathway. Affected genes were enriched in functionally important immune cell activation pathways such as antigen-mediated signaling, cell cycle control and mRNA processing and splicing. Discussion: Functional observations from donor B cell subsets in progressive states of differentiation and from timecourse experiments using the in vitro model suggest that these widespread changes in mRNA splicing play a role in preparing naïve B cells for the decisive step of antigen-mediated activation and differentiation.

The HTLV-1 viral oncoproteins Tax and HBZ reprogram the cellular mRNA splicing landscape.

Viral infections are known to hijack the transcription and translation of the host cell. However, the extent to which viral proteins coordinate these perturbations remains unclear. Here we used a model system, the human T-cell leukemia virus type 1 (HTLV-1), and systematically analyzed the transcriptome and interactome of key effectors oncoviral proteins Tax and HBZ. We showed that Tax and HBZ target distinct but also common transcription factors. Unexpectedly, we also uncovered a large set of interactions with RNA-binding proteins, including the U2 auxiliary factor large subunit (U2AF2), a key cellular regulator of pre-mRNA splicing. We discovered that Tax and HBZ perturb the splicing landscape by altering cassette exons in opposing manners, with Tax inducing exon inclusion while HBZ induces exon exclusion. Among Tax- and HBZ-dependent splicing changes, we identify events that are also altered in Adult T cell leukemia/lymphoma (ATLL) samples from two independent patient cohorts, and in well-known cancer census genes. Our interactome mapping approach, applicable to other viral oncogenes, has identified spliceosome perturbation as a novel mechanism coordinated by Tax and HBZ to reprogram the transcriptome.

ERG transcription factors have a splicing regulatory function involving RBFOX2 that is altered in the EWS-FLI1 oncogenic fusion.

ERG family proteins (ERG, FLI1 and FEV) are a subfamily of ETS transcription factors with key roles in physiology and development. In Ewing sarcoma, the oncogenic fusion protein EWS-FLI1 regulates both transcription and alternative splicing of pre-messenger RNAs. However, whether wild-type ERG family proteins might regulate splicing is unknown. Here, we show that wild-type ERG proteins associate with spliceosomal components, are found on nascent RNAs, and induce alternative splicing when recruited onto a reporter minigene. Transcriptomic analysis revealed that ERG and FLI1 regulate large numbers of alternative spliced exons (ASEs) enriched with RBFOX2 motifs and co-regulated by this splicing factor. ERG and FLI1 are associated with RBFOX2 via their conserved carboxy-terminal domain, which is present in EWS-FLI1. Accordingly, EWS-FLI1 is also associated with RBFOX2 and regulates ASEs enriched in RBFOX2 motifs. However, in contrast to wild-type ERG and FLI1, EWS-FLI1 often antagonizes RBFOX2 effects on exon inclusion. In particular, EWS-FLI1 reduces RBFOX2 binding to the ADD3 pre-mRNA, thus increasing its long isoform, which represses the mesenchymal phenotype of Ewing sarcoma cells. Our findings reveal a RBFOX2-mediated splicing regulatory function of wild-type ERG family proteins, that is altered in EWS-FLI1 and contributes to the Ewing sarcoma cell phenotype.

Genome-wide Transcript Structure Resolution Reveals Abundant Alternate Isoform Usage from Murine Gammaherpesvirus 68.

The gammaherpesviruses, including Epstein-Barr virus (EBV), Kaposi's sarcoma-associated herpesvirus (KSHV), and murine gammaherpesvirus 68 (MHV68, MuHV-4, γHV68), are etiologic agents of a wide range of lymphomas and non-hematological malignancies. These viruses possess large and highly dense dsDNA genomes that feature >80 bidirectionally positioned open reading frames (ORFs). The abundance of overlapping transcripts and extensive splicing throughout these genomes have until now prohibited high throughput-based resolution of transcript structures. Here, we integrate the capabilities of long-read sequencing with the accuracy of short-read platforms to globally resolve MHV68 transcript structures using the transcript resolution through integration of multi-platform data (TRIMD) pipeline. This approach reveals highly complex features, including: (1) pervasive overlapping transcript structures; (2) transcripts containing intra-gene or trans-gene splices that yield chimeric ORFs; (3) antisense and intergenic transcripts containing ORFs; and (4) noncoding transcripts. This work sheds light on the underappreciated complexity of gammaherpesvirus transcription and provides an extensively revised annotation of the MHV68 transcriptome.

Gammaherpesvirus Readthrough Transcription Generates a Long Non-Coding RNA That Is Regulated by Antisense miRNAs and Correlates with Enhanced Lytic Replication In Vivo.

Gammaherpesviruses, including the human pathogens Epstein⁻Barr virus (EBV) and Kaposi's sarcoma-associated herpesvirus (KSHV) are oncogenic viruses that establish lifelong infections in hosts and are associated with the development of lymphoproliferative diseases and lymphomas. Recent studies have shown that the majority of the mammalian genome is transcribed and gives rise to numerous long non-coding RNAs (lncRNAs). Likewise, the large double-stranded DNA virus genomes of herpesviruses undergo pervasive transcription, including the expression of many as yet uncharacterized lncRNAs. Murine gammaperherpesvirus 68 (MHV68, MuHV-4, HV68) is a natural pathogen of rodents, and is genetically and pathogenically related to EBV and KSHV, providing a highly tractable model for studies of gammaherpesvirus biology and pathogenesis. Through the integrated use of parallel data sets from multiple sequencing platforms, we previously resolved transcripts throughout the MHV68 genome, including at least 144 novel transcript isoforms. Here, we sought to molecularly validate novel transcripts identified within the M3/M2 locus, which harbors genes that code for the chemokine binding protein M3, the latency B cell signaling protein M2, and 10 microRNAs (miRNAs). Using strand-specific northern blots, we validated the presence of M3-04, a 3.91 kb polyadenylated transcript that initiates at the M3 transcription start site and reads through the M3 open reading frame (ORF), the M3 poly(a) signal sequence, and the M2 ORF. This unexpected transcript was solely localized to the nucleus, strongly suggesting that it is not translated and instead may function as a lncRNA. Use of an MHV68 mutant lacking two M3-04-antisense pre-miRNA stem loops resulted in highly increased expression of M3-04 and increased virus replication in the lungs of infected mice, demonstrating a key role for these RNAs in regulation of lytic infection. Together these findings suggest the possibility of a tripartite regulatory relationship between the lncRNA M3-04, antisense miRNAs, and the latency gene M2.

Analysis of EBV Transcription Using High-Throughput RNA Sequencing.

High-throughput sequencing of RNA is used to analyze the transcriptomes of viruses and cells, providing information about transcript structure and abundance. A wide array of programs and pipelines has been created to manage and interpret the abundance of data generated from high-throughput RNA sequencing experiments. This protocol details the use of free and open-source programs to align RNA-Seq reads to a reference genome, visualize read coverage and splice junctions, estimate transcript abundance, and evaluate differential expression of transcripts in different conditions. Particular concerns related to EBV and viral transcriptomics are addressed and access to EBV reference files is provided.

Global transcript structure resolution of high gene density genomes through multi-platform data integration.

Annotation of herpesvirus genomes has traditionally been undertaken through the detection of open reading frames and other genomic motifs, supplemented with sequencing of individual cDNAs. Second generation sequencing and high-density microarray studies have revealed vastly greater herpesvirus transcriptome complexity than is captured by existing annotation. The pervasive nature of overlapping transcription throughout herpesvirus genomes, however, poses substantial problems in resolving transcript structures using these methods alone. We present an approach that combines the unique attributes of Pacific Biosciences Iso-Seq long-read, Illumina short-read and deepCAGE (Cap Analysis of Gene Expression) sequencing to globally resolve polyadenylated isoform structures in replicating Epstein-Barr virus (EBV). Our method, Transcriptome Resolution through Integration of Multi-platform Data (TRIMD), identifies nearly 300 novel EBV transcripts, quadrupling the size of the annotated viral transcriptome. These findings illustrate an array of mechanisms through which EBV achieves functional diversity in its relatively small, compact genome including programmed alternative splicing (e.g. across the IR1 repeats), alternative promoter usage by LMP2 and other latency-associated transcripts, intergenic splicing at the BZLF2 locus, and antisense transcription and pervasive readthrough transcription throughout the genome.

New Noncoding Lytic Transcripts Derived from the Epstein-Barr Virus Latency Origin of Replication, oriP, Are Hyperedited, Bind the Paraspeckle Protein, NONO/p54nrb, and Support Viral Lytic Transcription.

UNLABELLED: We have previously shown that the Epstein-Barr virus (EBV) likely encodes hundreds of viral long noncoding RNAs (vlncRNAs) that are expressed during reactivation. Here we show that the EBV latency origin of replication (oriP) is transcribed bi-directionally during reactivation and that both leftward (oriPtLs) and rightward (oriPtRs) transcripts are largely localized in the nucleus. While the oriPtLs are most likely noncoding, at least some of the oriPtRs contain the BCRF1/vIL10 open reading frame. Nonetheless, oriPtR transcripts with long 5' untranslated regions may partially serve noncoding functions. Both oriPtL and oriPtR transcripts are expressed with late kinetics, and their expression is inhibited by phosphonoacetic acid. RNA sequencing (RNA-seq) analysis showed that oriPtLs and oriPtRs exhibited extensive "hyperediting" at their Family of Repeat (FR) regions. RNA secondary structure prediction revealed that the FR region of both oriPtLs and oriPtRs may form large evolutionarily conserved and thermodynamically stable hairpins. The double-stranded RNA-binding protein and RNA-editing enzyme ADAR was found to bind to oriPtLs, likely facilitating editing of the FR hairpin. Further, the multifunctional paraspeckle protein, NONO, was found to bind to oriPt transcripts, suggesting that oriPts interact with the paraspeckle-based innate antiviral immune pathway. Knockdown and ectopic expression of oriPtLs showed that it contributes to global viral lytic gene expression and viral DNA replication. Together, these results show that these new vlncRNAs interact with cellular innate immune pathways and that they help facilitate progression of the viral lytic cascade. IMPORTANCE: Recent studies have revealed that the complexity of lytic herpesviral transcriptomes is significantly greater than previously appreciated with hundreds of viral long noncoding RNAs (vlncRNAs) being recently discovered. Work on cellular lncRNAs over the past several years has just begun to give us an initial appreciation for the array of functions they play in complex formation and regulatory processes in the cell. The newly identified herpesvirus lncRNAs are similarly likely to play a variety of different functions, although these functions are likely tailored to specific needs of the viral infection cycles. Here we describe novel transcripts derived from the EBV latency origin of replication. We show that they are hyperedited, that they interact with a relatively newly appreciated antiviral pathway, and that they play a role in facilitating viral lytic gene expression. These investigations are a starting point to unraveling the complex arena of vlncRNA function in herpesvirus lytic replication.

High-throughput RNA sequencing-based virome analysis of 50 lymphoma cell lines from the Cancer Cell Line Encyclopedia project.

UNLABELLED: Using high-throughput RNA sequencing data from 50 common lymphoma cell culture models from the Cancer Cell Line Encyclopedia project, we performed an unbiased global interrogation for the presence of a panel of 740 viruses and strains known to infect human and other mammalian cells. This led to the findings of previously identified infections by Epstein-Barr virus (EBV), Kaposi's sarcoma herpesvirus (KSHV), and human T-lymphotropic virus type 1 (HTLV-1). In addition, we also found a previously unreported infection of one cell line (DEL) with a murine leukemia virus. High expression of murine leukemia virus (MuLV) transcripts was observed in DEL cells, and we identified four transcriptionally active integration sites, one being in the TNFRSF6B gene. We also found low levels of MuLV reads in a number of other cell lines and provided evidence suggesting cross-contamination during sequencing. Analysis of HTLV-1 integrations in two cell lines, HuT 102 and MJ, identified 14 and 66 transcriptionally active integration sites with potentially activating integrations in immune regulatory genes, including interleukin-15 (IL-15), IL-6ST, STAT5B, HIVEP1, and IL-9R. Although KSHV and EBV do not typically integrate into the genome, we investigated a previously identified integration of EBV into the BACH2 locus in Raji cells. This analysis identified a BACH2 disruption mechanism involving splice donor sequestration. Through viral gene expression analysis, we detected expression of stable intronic RNAs from the EBV BamHI W repeats that may be part of long transcripts spanning the repeat region. We also observed transcripts at the EBV vIL-10 locus exclusively in the Hodgkin's lymphoma cell line, Hs 611.T, the expression of which were uncoupled from other lytic genes. Assessment of the KSHV viral transcriptome in BCP-1 cells showed expression of the viral immune regulators, K2/vIL-6, K4/vIL-8-like vCCL1, and K5/E2-ubiquitin ligase 1 that was significantly higher than expression of the latency-associated nuclear antigen. Together, this investigation sheds light into the virus composition across these lymphoma model systems and provides insights into common viral mechanistic principles. IMPORTANCE: Viruses cause cancer in humans. In lymphomas the Epstein-Barr virus (EBV), Kaposi's sarcoma herpesvirus (KSHV) and human T-lymphotropic virus type 1 are major contributors to oncogenesis. We assessed virus-host interactions using a high throughput sequencing method that facilitates the discovery of new virus-host associations and the investigation into how the viruses alter their host environment. We found a previously unknown murine leukemia virus infection in one cell line. We identified cellular genes, including cytokine regulators, that are disrupted by virus integration, and we determined mechanisms through which virus integration causes deregulation of cellular gene expression. Investigation into the KSHV transcriptome in the BCP-1 cell line revealed high-level expression of immune signaling genes. EBV transcriptome analysis showed expression of vIL-10 transcripts in a Hodgkin's lymphoma that was uncoupled from lytic genes. These findings illustrate unique mechanisms of viral gene regulation and to the importance of virus-mediated host immune signaling in lymphomas.

Global bidirectional transcription of the Epstein-Barr virus genome during reactivation.

Epstein-Barr virus (EBV) reactivation involves the ordered induction of approximately 90 viral genes that participate in the generation of infectious virions. Using strand-specific RNA-seq to assess the EBV transcriptome during reactivation, we found extensive bidirectional transcription extending across nearly the entire genome. In contrast, only 4% of the EBV genome is currently bidirectionally annotated. Most of the newly identified transcribed regions show little evidence of coding potential, supporting noncoding roles for most of these RNAs. Based on previous cellular long noncoding RNA size calculations, we estimate that there are likely hundreds more EBV genes expressed during reactivation than was previously known. Limited 5' and 3' rapid amplification of cDNA ends (RACE) experiments and findings of novel splicing events by RNA-seq suggest that the complexity of the viral genome during reactivation may be even greater. Further analysis of antisense transcripts at some of the EBV latency gene loci showed that they are "late" genes, they are nuclear, and they tend to localize in areas of the nucleus where others find newly synthesized viral genomes. This raises the possibility that these transcripts perform functions such as new genome processing, stabilization, organization, etc. The finding of a significantly more complex EBV transcriptome during reactivation changes our view of the viral production process from one that is facilitated and regulated almost entirely by previously identified viral proteins to a process that also involves the contribution of a wide array of virus encoded noncoding RNAs. Epstein-Barr virus (EBV) is a herpesvirus that infects the majority of the world's population, in rare cases causing serious disease such as lymphoma and gastric carcinoma. Using strand-specific RNA-seq, we have studied viral gene expression during EBV reactivation and have discovered hundreds more viral transcripts than were previously known. The finding of alternative splicing and the prevalence of overlapping transcripts indicate additional complexity. Most newly identified transcribed regions do not encode proteins but instead likely function as noncoding RNA molecules which could participate in regulating gene expression, gene splicing or even activities such as viral genome processing. These findings broaden the scope of what we need to consider to understand the viral manufacturing process. As more detailed studies are undertaken they will likely change the way we view this process as a whole.

Epstein-Barr virus and human herpesvirus 6 detection in a non-Hodgkin's diffuse large B-cell lymphoma cohort by using RNA sequencing.

Comprehensive virome analysis of RNA sequence (RNA-seq) data sets from 118 non-Hodgkin's B-cell lymphomas revealed a small subset that is positive for Epstein-Barr virus (EBV) or human herpesvirus 6B (HHV-6B), with one coinfection. EBV transcriptome analysis revealed expression of the latency genes RPMS1, LMP1, and LMP2, with one sample additionally showing a high level of early lytic expression and another sample showing a high level of EBNA2 expression. HHV-6B transcriptome analysis revealed that the majority of genes were transcribed.