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.

Recent studies have identified circular RNAs (circRNAs) expressed from the Epstein-Barr virus (EBV) and Kaposi’s sarcoma herpesvirus (KSHV) human DNA tumor viruses. To gain initial insights into the potential relevance of EBV circRNAs in virus biology and disease, we assessed the circRNAome of the interspecies homologue rhesus macaque lymphocryptovirus (rLCV) in a naturally occurring lymphoma from a simian immunodeficiency virus (SIV)-infected rhesus macaque. This analysis revealed rLCV orthologues of the latency-associated EBV circular RNAs circRPMS1_E4_E3a and circEBNA_U. Also identified in two samples displaying unusually high lytic gene expression was a novel rLCV circRNA that contains both conserved and rLCV-specific RPMS1 exons and whose backsplice junctions flank an rLCV lytic origin of replication (OriLyt). Analysis of a lytic infection model for the murid herpesvirus 68 (MHV68) rhadinovirus identified a cluster of circRNAs near an MHV68 lytic origin of replication, with the most abundant of these, circM11_ORF69, spanning the OriLyt. Lastly, analysis of KSHV latency and reactivation models revealed the latency associated circRNA originating from the vIRF4 gene as the predominant viral circRNA. Together, the results of this study broaden our appreciation for circRNA repertoires in the Lymphocryptovirus and Rhadinovirus genera of gammaherpesviruses and provide evolutionary support for viral circRNA functions in latency and viral replication. IMPORTANCE Infection with oncogenic gammaherpesviruses leads to long-term viral persistence through a dynamic interplay between the virus and the host immune system. Critical for remodeling of the host cell environment after the immune responses are viral noncoding RNAs that modulate host signaling pathways without attracting adaptive immune recognition. Despite the importance of noncoding RNAs in persistent infection, the circRNA class of noncoding RNAs has only recently been identified in gammaherpesviruses. Accordingly, their roles in virus infection and associated oncogenesis are unknown. Here we report evolutionary conservation of EBV-encoded circRNAs determined by assessing the circRNAome in rLCV-infected lymphomas from an SIV-infected rhesus macaque, and we report latent and lytic circRNAs from KSHV and MHV68. These experiments demonstrate utilization of the circular RNA class of RNAs across 4 members of the gammaherpesvirus subfamily, and they identify orthologues and potential homoplastic circRNAs, implying conserved circRNA functions in virus biology and associated malignancies.

Deregulated expression of circular RNAs (circRNAs) is associated with various human diseases, including many types of cancer. Despite their growing links to cancer, there has been limited characterization of circRNAs in metastatic castration-resistant prostate cancer, the major cause of prostate cancer mortality. Here, through the analysis of an exome-capture RNA-seq dataset from 47 metastatic castration-resistant prostate cancer samples and ribodepletion and RNase R RNAsequencing of patient-derived xenografts (PDXs) and cell models, we identified 13 circRNAs generated from the key prostate cancer driver gene-androgen receptor (AR). We validated and characterized the top four most abundant, clinically relevant AR circRNAs. Expression of these AR circRNAs was upregulated during castration-resistant progression of PDXs. The upregulation was not due to global increase of circRNA formation in these tumors. Instead, the levels of AR circRNAs correlated strongly with that of the linear AR transcripts (both AR and AR-variants) in clinical samples and PDXs, indicating a transcriptional mechanism of regulation. In cultured cells, androgen suppressed the expression of these AR circRNAs and the linear AR transcripts, and the suppression was attenuated by an antiandrogen. Using nuclear/cytoplasmic fractionation and RNA Users may view, print, copy, and download text and data-mine the content in such documents, for the purposes of academic research, subject always to the full Conditions of use:

Certain peripheral T-cell lymphomas (PTCLs) have been associated with viral infection, particularly infection with Epstein-Barr virus (EBV). However, a comprehensive virome analysis across PTCLs has not previously been reported. Here we utilized published whole-transcriptome RNA sequencing (RNA-seq) data sets from seven different PTCL studies and new RNA-seq data from our laboratory to screen for virus association, to analyze viral gene expression, and to assess B- and T-cell receptor diversity paradigms across PTCL subtypes. In addition to identifying EBV in angioimmunoblastic T-cell lymphoma (AITL) and extranodal NK/T-cell lymphoma (ENKTL), two PTCL subtypes with well-established EBV associations, we also detected EBV in several cases of anaplastic large-cell lymphoma (ALCL), and we found evidence of infection by the oncogenic viruses Kaposi’s sarcoma-associated herpesvirus and human T-cell leukemia virus type 1 in isolated PTCL cases. In AITLs, EBV gene expression analysis showed expression of immediate early, early, and late lytic genes, suggesting either low-level lytic gene expression or productive infection in a subset of EBV-infected B-lymphocyte stromal cells. Deconvolution of immune cell subpopulations demonstrated a greater B-cell signal in AITLs than in other PTCL subtypes, consistent with a larger role for B-cell support in the pathogenesis of AITL. Reconstructed T-cell receptor (TCR) and B-cell receptor (BCR) repertoires demonstrated increased BCR diversity in AITLs, consistent with a possible EBV-driven polyclonal response. These findings indicate potential alternative roles for EBV in PTCLs, in addition to the canonical oncogenic mechanisms associated with EBV latent infection. Our findings also suggest the involvement of other viruses in PTCL pathogenesis and demonstrate immunological alterations associated with these cancers.IMPORTANCEIn this study, we utilized next-generation sequencing data from 7 different studies of peripheral T-cell lymphoma (PTCL) patient samples to globally assess viral associations, provide insights into the contributions of EBV gene expression to the tumor phenotype, and assess the unique roles of EBV in modulating the immune cell tumor microenvironment. These studies revealed potential roles for EBV replication genes in some PTCL subtypes, the possible role of additional human tumor viruses in rare cases of PTCLs, and a role for EBV in providing a unique immune microenvironmental niche in one subtype of PTCLs. Together, these studies provide new insights into the understudied role of tumor viruses in PTCLs.

Yeast genomes contain variable number tandem repeats (VNTRs) within coding regions of DNA. A significant number of these genes are involved in cell rescue, defence and virulence and are regulated by genetic elements associated with stress. Alleles that encode variable length, single amino acid tracts, are mainly associated with transcription and proteins localized within the nucleus. Alleles that encode proteins containing oligopeptide repeats or minisatellites are over-represented in cell wall and extracellular space locations. Functional analysis of the latter group reveals that these proteins are involved in biogenesis of cellular components and in interaction with the cellular environment, especially in relation to stress resistance, heat shock response, temperature perception and adhesion. A significantly high number of these proteins have regions rich in threonine and/or serine that contain repeated sequences, variable in length within yeast species. DNA sequences encoding serine-and/or threonine-rich regions give rise to polymorphic alleles and therefore may confer a selective advantage to cells. We propose that these regions are the focus of mutational and recombination events that, when coupled with directed selection, may contribute to genetic variation within stress-related genes.