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...

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 doublestranded 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. IMPORTANCERecent 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. The Epstein-Barr virus (EBV) is an oncogenic human herpesvirus that infects 95% of the world's population (1). Although infection is typically asymptomatic, the virus will persist in the host for the duration of its life. In most immunocompetent individuals, the virus poses little serious health risk and its presence is unremarkable. Nevertheless, in...

Androgen receptor splice variants are alternatively spliced variants of androgen receptor that are C-terminally truncated and lack the canonical ligand-binding domain. Accumulating evidence has indicated a significant role of androgen receptor splice variants in mediating resistance of castration-resistant prostate cancer to current therapies and in predicting therapeutic responses. As such, there is an urgent need to target androgen receptor splicing variants for more effective treatment of castration-resistant prostate cancer. Identification of precise and critical targeting points to deactivate androgen receptor splicing variants relies on a deep understanding of how they are generated and the mechanisms of their action. In this review, we will focus on the emerging data on their generation, clinical significance, and mechanisms of action as well as the therapeutic impact of these findings.

Epigenetic regulation of gene expression is critical to phenotypic maintenance and transition of human breast cancer cells. HOX antisense intergenic RNA (HOTAIR) is a long intergenic non-coding RNA that epigenetically represses gene expression via recruitment of enhancer of zeste homolog 2 (EZH2), a histone methyltransferase. Elevated expression of HOTAIR promotes progression of breast cancer. In the current study we examined the expression and function of HOTAIR in MCF-7-TNR cells, a derivative of the luminal-like breast cancer cell line MCF-7 that acquired resistance to TNF-α-induced cell death. The expression of HOTAIR, markers of the luminal-like and basal-like subtypes, and growth were compared between MCF-7 and MCF-7-TNR cells. These variables were further assessed upon inhibition of HOTAIR, EZH2, p38 MAPK, and SRC kinase in MCF-7-TNR cells. When compared with MCF-7 cells, MCF-7-TNR cells exhibited an increase in the expression of HOTAIR, which correlated with characteristics of a luminal-like to basal-like transition as evidenced by dysregulated gene expression and accelerated growth. MCF-7-TNR cells exhibited reduced suppressive histone H3 lysine27 trimethylation on the HOTAIR promoter. Inhibition of HOTAIR and EZH2 attenuated the luminal-like to basal-like transition in terms of gene expression and growth in MCF-7-TNR cells. Inhibition of p38 and SRC diminished HOTAIR expression and the basal-like phenotype in MCF-7-TNR cells. HOTAIR was robustly expressed in the native basal-like breast cancer cells and inhibition of HOTAIR reduced the basal-like gene expression and growth. Our findings suggest HOTAIR-mediated regulation of gene expression and growth associated with the basal-like phenotype of breast cancer cells.