RESUMEN
Kaposi's sarcoma-associated herpesvirus is the etiologic agent of Kaposi's sarcoma and two B-cell malignancies. Recent advancements in sequencing technologies have led to high resolution transcriptomes for several human herpesviruses that densely encode genes on both strands. However, for KSHV progress remained limited due to the overall low percentage of KSHV transcripts, even during lytic replication. To address this challenge, we have developed a target enrichment method to increase the KSHV-specific reads for both short- and long-read sequencing platforms. Furthermore, we combined this approach with the Transcriptome Resolution through Integration of Multi-platform Data (TRIMD) pipeline developed previously to annotate transcript structures. TRIMD first builds a scaffold based on long-read sequencing and validates each transcript feature with supporting evidence from Illumina RNA-Seq and deepCAGE sequencing data. Our stringent innovative approach identified 994 unique KSHV transcripts, thus providing the first high-density KSHV lytic transcriptome. We describe a plethora of novel coding and non-coding KSHV transcript isoforms with alternative untranslated regions, splice junctions and open-reading frames, thus providing deeper insights on gene expression regulation of KSHV. Interestingly, as described for Epstein-Barr virus, we identified transcription start sites that augment long-range transcription and may increase the number of latency-associated genes potentially expressed in KS tumors.
Asunto(s)
Empalme Alternativo , Herpesvirus Humano 8 , Transcriptoma , Herpesvirus Humano 8/genética , Humanos , Transcriptoma/genética , Transcripción Genética , Regulación Viral de la Expresión Génica , Sistemas de Lectura Abierta/genética , Secuenciación de Nucleótidos de Alto Rendimiento , Sarcoma de Kaposi/virología , Sarcoma de Kaposi/genética , ARN Viral/genética , ARN Viral/metabolismoRESUMEN
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.
Asunto(s)
Infecciones por Virus de Epstein-Barr , Herpesvirus Humano 4 , Iniciación de la Transcripción Genética , Replicación Viral , Humanos , Herpesvirus Humano 4/fisiología , Interacciones Huésped-Patógeno , Regiones Promotoras Genéticas , TATA Box , Factores de Transcripción/metabolismo , Sitio de Iniciación de la Transcripción , Transcripción Genética , Proteínas Virales/metabolismo , Proteínas Virales/genética , Infecciones por Virus de Epstein-Barr/metabolismo , Infecciones por Virus de Epstein-Barr/virologíaRESUMEN
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.
RESUMEN
As a fundamental aspect of normal cell signaling and disease states, there is great interest in determining alternative splicing (AS) changes in physiologic, pathologic, and pharmacologic settings. High throughput RNA sequencing and specialized software to detect AS has greatly enhanced our ability to determine transcriptome-wide splicing changes. Despite the richness of this data, deriving meaning from sometimes thousands of AS events is a substantial bottleneck for most investigators. We present SpliceTools, a suite of data processing modules that arms investigators with the ability to quickly produce summary statistics, mechanistic insights, and functional significance of AS changes through command line or through an online user interface. Utilizing RNA-seq datasets for 186 RNA binding protein knockdowns, nonsense mediated RNA decay inhibition, and pharmacologic splicing inhibition, we illustrate the utility of SpliceTools to distinguish splicing disruption from regulated transcript isoform changes, we show the broad transcriptome footprint of the pharmacologic splicing inhibitor, indisulam, we illustrate the utility in uncovering mechanistic underpinnings of splicing inhibition, we identify predicted neo-epitopes in pharmacologic splicing inhibition, and we show the impact of splicing alterations induced by indisulam on cell cycle progression. Together, SpliceTools puts rapid and easy downstream analysis at the fingertips of any investigator studying AS.
Asunto(s)
Empalme Alternativo , Empalme del ARN , Empalme Alternativo/genética , Sulfonamidas , Transcriptoma/genética , Análisis de Secuencia de ARN/métodosRESUMEN
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.
Asunto(s)
Empalme Alternativo , Cisplatino , Resistencia a Antineoplásicos , Proteínas de Unión al ARN , Factores de Transcripción , Empalme Alternativo/genética , Cisplatino/farmacología , Cisplatino/metabolismo , Daño del ADN , Proteínas Nucleares/metabolismo , Precursores del ARN/genética , Precursores del ARN/metabolismo , Proteínas de Unión al ARN/genética , Proteínas de Unión al ARN/metabolismo , Factores de Transcripción/metabolismo , Línea Celular Tumoral , Humanos , AnimalesRESUMEN
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.
Asunto(s)
Vesículas Extracelulares , MicroARNs , Neoplasias , Animales , Modelos Animales de Enfermedad , Células Endoteliales , Vesículas Extracelulares/genética , Ratones , MicroARNs/genética , Microambiente Tumoral , Macrófagos Asociados a TumoresRESUMEN
BACKGROUND: Extracellular vesicles (EVs) are released by nearly every cell type and have attracted much attention for their ability to transfer protein and diverse RNA species from donor to recipient cells. Much attention has been given so far to the features of EV short RNAs such as miRNAs. However, while the presence of mRNA and long noncoding RNA (lncRNA) transcripts in EVs has also been reported by multiple different groups, the properties and function of these longer transcripts have been less thoroughly explored than EV miRNA. Additionally, the impact of EV export on the transcriptome of exporting cells has remained almost completely unexamined. Here, we globally investigate mRNA and lncRNA transcripts in endothelial EVs in multiple different conditions. RESULTS: In basal conditions, long RNA transcripts enriched in EVs have longer than average half-lives and distinctive stability-related sequence and structure characteristics including shorter transcript length, higher exon density, and fewer 3' UTR A/U-rich elements. EV-enriched long RNA transcripts are also enriched in HNRNPA2B1 binding motifs and are impacted by HNRNPA2B1 depletion, implicating this RNA-binding protein in the sorting of long RNA to EVs. After signaling-dependent modification of the cellular transcriptome, we observed that, unexpectedly, the rate of EV enrichment relative to cells was altered for many mRNA and lncRNA transcripts. This change in EV enrichment was negatively correlated with intracellular abundance, with transcripts whose export to EVs increased showing decreased abundance in cells and vice versa. Correspondingly, after treatment with inhibitors of EV secretion, levels of mRNA and lncRNA transcripts that are normally highly exported to EVs increased in cells, indicating a measurable impact of EV export on the long RNA transcriptome of the exporting cells. Compounds with different mechanisms of inhibition of EV secretion affected the cellular transcriptome differently, suggesting the existence of multiple EV subtypes with different long RNA profiles. CONCLUSIONS: We present evidence for an impact of EV physiology on the characteristics of EV-producing cell transcriptomes. Our work suggests a new paradigm in which the sorting and packaging of transcripts into EVs participate, together with transcription and RNA decay, in controlling RNA homeostasis and shape the cellular long RNA abundance profile.
Asunto(s)
Vesículas Extracelulares , MicroARNs , ARN Largo no Codificante , Movimiento Celular , Vesículas Extracelulares/genética , Vesículas Extracelulares/metabolismo , MicroARNs/metabolismo , ARN Largo no Codificante/genética , ARN Largo no Codificante/metabolismo , ARN Mensajero/genética , ARN Mensajero/metabolismoRESUMEN
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.
Asunto(s)
Empalme Alternativo , Infecciones por Virus de Epstein-Barr , Humanos , Herpesvirus Humano 4/genética , ARN Mensajero/genética , Diferenciación Celular/genéticaRESUMEN
The U6 snRNA, the core catalytic component of the spliceosome, is extensively modified post-transcriptionally, with 2'-O-methylation being most common. However, how U6 2'-O-methylation is regulated remains largely unknown. Here we report that TFIP11, the human homolog of the yeast spliceosome disassembly factor Ntr1, localizes to nucleoli and Cajal Bodies and is essential for the 2'-O-methylation of U6. Mechanistically, we demonstrate that TFIP11 knockdown reduces the association of U6 snRNA with fibrillarin and associated snoRNAs, therefore altering U6 2'-O-methylation. We show U6 snRNA hypomethylation is associated with changes in assembly of the U4/U6.U5 tri-snRNP leading to defects in spliceosome assembly and alterations in splicing fidelity. Strikingly, this function of TFIP11 is independent of the RNA helicase DHX15, its known partner in yeast. In sum, our study demonstrates an unrecognized function for TFIP11 in U6 snRNP modification and U4/U6.U5 tri-snRNP assembly, identifying TFIP11 as a critical spliceosome assembly regulator.
Asunto(s)
Factores de Empalme de ARN/metabolismo , Empalme del ARN/fisiología , ARN Nuclear Pequeño/metabolismo , Ribonucleoproteína Nuclear Pequeña U4-U6/metabolismo , Ribonucleoproteína Nuclear Pequeña U5/metabolismo , Nucléolo Celular/metabolismo , Supervivencia Celular , Cuerpos Enrollados/metabolismo , Células HeLa , Humanos , Metilación , Mitosis , Proteínas Nucleares/metabolismo , Motas Nucleares/metabolismo , Unión Proteica , Estabilidad Proteica , Precursores del ARN/metabolismo , Factores de Empalme de ARN/genética , ARN Nucleolar Pequeño/metabolismo , Empalmosomas/metabolismoRESUMEN
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.
Asunto(s)
Factores de Transcripción con Cremalleras de Leucina de Carácter Básico/metabolismo , Productos del Gen tax/metabolismo , Infecciones por HTLV-I/metabolismo , Leucemia-Linfoma de Células T del Adulto/virología , Proteínas de los Retroviridae/metabolismo , Células HEK293 , Infecciones por HTLV-I/etiología , Virus Linfotrópico T Tipo 1 Humano , Humanos , Células Jurkat , Empalme del ARN , ARN Mensajero , Factor de Empalme U2AF/metabolismoRESUMEN
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.
Asunto(s)
Empalme Alternativo , Proteínas de Fusión Oncogénica/metabolismo , Proteína Proto-Oncogénica c-fli-1/metabolismo , Factores de Empalme de ARN/metabolismo , Proteína EWS de Unión a ARN/metabolismo , Proteínas Represoras/metabolismo , Proteínas de Unión a Calmodulina/genética , Proteínas de Unión a Calmodulina/metabolismo , Línea Celular , Línea Celular Tumoral , Células HeLa , Células Endoteliales de la Vena Umbilical Humana/metabolismo , Humanos , Dominios Proteicos , Sarcoma de Ewing/genética , Sarcoma de Ewing/metabolismo , Regulador Transcripcional ERG/química , Regulador Transcripcional ERG/metabolismoRESUMEN
Muscle formation is controlled by a number of key myogenic transcriptional regulators that govern stage-specific gene expression programs and act as terminal effectors of intracellular signaling pathways. To date, the role of phosphatases in the signaling cascades instructing muscle development remains poorly understood. Here, we show that a specific PP2A-B55δ holoenzyme is necessary for skeletal myogenesis. The primary role of PP2A-B55δ is to dephosphorylate histone deacetylase 4 (HDAC4) following myocyte differentiation and ensure repression of Myocyte enhancer factor 2D (MEF2D)-dependent gene expression programs during myogenic fusion. As a crucial HDAC4/MEF2D target gene that governs myocyte fusion, we identify ArgBP2, an upstream inhibitor of Abl, which itself is a repressor of CrkII signaling. Consequently, cells lacking PP2A-B55δ show upregulation of ArgBP2 and hyperactivation of CrkII downstream effectors, including Rac1 and FAK, precluding cytoskeletal and membrane rearrangements associated with myoblast fusion. Both in vitro and in zebrafish, loss-of-function of PP2A-B55δ severely impairs fusion of myocytes and formation of multinucleated muscle fibers, without affecting myoblast differentiation. Taken together, our results establish PP2A-B55δ as the first protein phosphatase to be involved in myoblast fusion and suggest that reversible phosphorylation of HDAC4 may coordinate differentiation and fusion events during myogenesis.
Asunto(s)
Histona Desacetilasas/metabolismo , Factores de Transcripción MEF2/metabolismo , Mioblastos/citología , Mioblastos/metabolismo , Animales , Fusión Celular , Línea Celular , Citoesqueleto/metabolismo , Embrión no Mamífero/metabolismo , Holoenzimas/metabolismo , Ratones , Morfogénesis , Desarrollo de Músculos , Fibras Musculares Esqueléticas/citología , Fenotipo , Fosforilación , ARN Mensajero/genética , ARN Mensajero/metabolismo , Transducción de Señal , Transcripción Genética , Pez Cebra/embriología , Proteína de Unión al GTP rac1/metabolismoRESUMEN
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.
Asunto(s)
Gammaherpesvirinae/metabolismo , Infecciones por Herpesviridae/metabolismo , Sistemas de Lectura Abierta , ARN Viral/biosíntesis , Transcriptoma , Animales , Estudio de Asociación del Genoma Completo , Ratones , Células 3T3 NIHRESUMEN
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.
RESUMEN
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.
Asunto(s)
Perfilación de la Expresión Génica , Herpesvirus Humano 4/genética , Secuenciación de Nucleótidos de Alto Rendimiento , Transcriptoma , Biología Computacional/métodos , Regulación Viral de la Expresión Génica , Genoma Viral , Genómica/métodos , Humanos , Análisis de Secuencia de ARN , Programas Informáticos , Navegador WebRESUMEN
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.
Asunto(s)
Genoma Viral , Herpesvirus Humano 4/genética , ARN Mensajero/genética , Estadística como Asunto , Empalme Alternativo/genética , Línea Celular , ADN Intergénico/genética , Exones/genética , Humanos , Anotación de Secuencia Molecular , Poliadenilación/genética , Regiones Promotoras Genéticas , ARN Mensajero/metabolismo , Secuencias Repetitivas de Ácidos Nucleicos/genética , Reproducibilidad de los Resultados , Análisis de Secuencia de ARN , Sitio de Iniciación de la Transcripción , Transcripción Genética , Transcriptoma/genéticaRESUMEN
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.
Asunto(s)
Regulación Viral de la Expresión Génica , Herpesvirus Humano 4/genética , Herpesvirus Humano 4/fisiología , Proteínas Asociadas a Matriz Nuclear/metabolismo , Factores de Transcripción de Octámeros/metabolismo , ARN no Traducido/metabolismo , Proteínas de Unión al ARN/metabolismo , Transcripción Genética , Replicación Viral , Línea Celular , Proteínas de Unión al ADN , Interacciones Huésped-Patógeno , Humanos , Inmunidad Innata , Unión Proteica , ARN no Traducido/genética , Origen de RéplicaRESUMEN
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.
Asunto(s)
Herpesviridae/aislamiento & purificación , Secuenciación de Nucleótidos de Alto Rendimiento , Linfoma/virología , Retroviridae/aislamiento & purificación , Línea Celular Tumoral , Perfilación de la Expresión Génica , Herpesviridae/genética , Herpesviridae/fisiología , Humanos , Retroviridae/genética , Retroviridae/fisiología , Integración ViralRESUMEN
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.