Our appreciation for the extent of Epstein Barr virus (EBV) transcriptome complexity continues to grow through findings of EBV encoded microRNAs, new long non-coding RNAs as well as the more recent discovery of over a hundred new polyadenylated lytic transcripts. Here we report an additional layer to the EBV transcriptome through the identification of a repertoire of latent and lytic viral circular RNAs. Utilizing RNase R-sequencing with cell models representing latency types I, II, and III, we identified EBV encoded circular RNAs expressed from the latency Cp promoter involving backsplicing from the W1 and W2 exons to the C1 exon, from the EBNA BamHI U fragment exon, and from the latency long non-coding RPMS1 locus. In addition, we identified circular RNAs expressed during reactivation including backsplicing from exon 8 to exon 2 of the LMP2 gene and a highly expressed circular RNA derived from intra-exonic backsplicing within the BHLF1 gene. While expression of most of these circular RNAs was found to depend on the EBV transcriptional program utilized and the transcription levels of the associated loci, expression of LMP2 exon 8 to exon 2 circular RNA was found to be cell model specific. Altogether we identified over 30 unique EBV circRNAs candidates and we validated and determined the structural features, expression profiles and nuclear/cytoplasmic distributions of several predominant and notable viral circRNAs. Further, we show that two of the EBV circular RNAs derived from the RPMS1 locus are detected in EBV positive clinical stomach cancer specimens. This study increases the known EBV latency and lytic transcriptome repertoires to include viral circular RNAs and it provides an essential foundation and resource for investigations into the functions and roles of this new class of EBV transcripts in EBV biology and diseases.

Increased expression of the full-length androgen receptor (AR-FL) and AR splice variants (AR-Vs) drives the progression of castration-resistant prostate cancer (CRPC). The levels of AR-FL and AR-V transcripts are often tightly correlated in individual CRPC samples, yet our understanding of how their expression is co-regulated is limited. Here, we report a role of c-Myc in accounting for coordinated AR-FL and AR-V expression. Analysis of gene expression data from 159 metastatic CRPC samples and 2142 primary prostate tumors showed that the level of c-Myc is positively correlated with that of individual AR isoforms. A striking positive correlation also exists between the activity of the c-Myc pathway and the level of individual AR isoforms, between the level of c-Myc and the activity of the AR pathway, and between the activities of the two pathways. Moreover, the c-Myc signature is highly enriched in tumors expressing high levels of AR, as is the AR signature in c-Myc-high-expressing tumors. Using shRNA knockdown, we confirmed c-Myc regulation of expression and activity of AR-FL and AR-Vs in cell models and a patient-derived xenograft model. Mechanistically, c-Myc promotes the transcription of the AR gene and enhances the stability of the AR-FL and AR-V proteins without altering AR RNA splicing. Importantly, inhibiting c-Myc sensitizes enzalutamide-resistant cells to growth inhibition by enzalutamide. Overall, this study highlights a critical role of c-Myc in regulating the coordinated expression of AR-FL and AR-Vs that is commonly observed in CRPC and suggests the utility of targeting c-Myc as an adjuvant to AR-directed therapy.

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. IMPORTANCE 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. E pstein-Barr virus (EBV) is a human gammaherpesvirus that is prevalent in all human populations. Although infection is frequently asymptomatic, it has been linked to a number of serious diseases such as Burkitt's lymphoma, nasopharyngeal carcinoma, and posttransplant lymphoproliferative disorder (1). Like other herpesviruses, EBV can exist in both lytic and latent phases, with primary lytic infections often progressing to lifelong latent infections with generally sporadic subsequent episodes o...

Using an enhanced RNA-Seq pipeline to analyze Epstein-Barr virus (EBV) transcriptomes, we investigated viral and cellular gene expression in the Akata cell line following B-cell-receptor-mediated reactivation. Robust induction of EBV gene expression was observed, with most viral genes induced >200-fold and with EBV transcripts accounting for 7% of all mapped reads within the cell. After induction, hundreds of candidate splicing events were detected using the junction mapper TopHat, including a novel nonproductive splicing event at the gp350/gp220 locus and several alternative splicing events at the LMP2 locus. A more detailed analysis of lytic LMP2 transcripts showed an overall lack of the prototypical type III latency splicing events. Analysis of nuclear versus cytoplasmic RNA-Seq data showed that the lytic forms of LMP2, EBNA-2, EBNA-LP, and EBNA-3A, -3B, and -3C have higher nuclear-to-cytoplasmic accumulation ratios than most lytic genes, including classic late genes. These data raise the possibility that at least some lytic transcripts derived from these latency gene loci may have unique, noncoding nuclear functions during reactivation. Our analysis also identified two previously unknown genes, BCLT1 and BCRT2, that map to the BamHI C-region of the EBV genome. Pathway analysis of cellular gene expression changes following B-cell receptor activation identified an inflammatory response as the top predicted function and ILK and TREM1 as the top predicted canonical pathways. E pstein-Barr virus (EBV) is a human pathogen that causes malignancies including Burkitt's lymphoma, Hodgkin's disease, and nasopharyngeal carcinoma (13). EBV has a complex infection cycle involving a number of different viral gene expression programs. These individual programs facilitate distinct tasks that are required for specific infection stages. Like all herpesviruses, EBV utilizes both latent gene expression programs, in which only limited numbers of viral genes are expressed, and a replicative gene expression program, in which the bulk of EBV genes are expressed to produce infectious virus.Efficient and synchronous virus reactivation can be modeled by activating the B-cell receptor (BCR) in the EBV-positive Burkitt's lymphoma cell line Akata (15). In this system, reactivation leads to an ordered induction of immediate-early (e.g., BZLF1 and BRLF1), early (e.g., BMRF1), and late genes, with immediate-early genes peaking at approximately 2 to 6 h and late genes peaking at approximately 6 to 24 h postinduction (6,15,20). Interestingly, Yuan et al. (20) found that EBV latency genes were induced following reactivation, suggesting a role for these genes in the lytic cycle (20).Second-generation RNA-Seq technology allows the simultaneous interrogation of gene expression and transcript structure at a high level of accuracy and at a single-nucleotide resolution. We have recently shown the application of RNA-Seq to the interrogation of EBV transcriptomes in two type I latency cell lines, Akata and Mutu I (11). Here we have improved our RNA-Se...

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 crosscontamination 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. IMPORTANCEViruses 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 sh...