Transcriptome-Wide, Unbiased Profiling of Ribonuclease Targeting Chimeras.

Various approaches have been developed to target RNA and modulate its function with modes of action including binding and cleavage. Herein, we explored how small molecule binding is correlated with cleavage induced by heterobifunctional ribonuclease targeting chimeras (RiboTACs), where RNase L is recruited to cleave the bound RNA target, in a transcriptome-wide, unbiased fashion. Only a fraction of bound targets was cleaved by RNase L, induced by RiboTAC binding. Global analysis suggested that (i) cleaved targets generally form a region of stable structure that encompasses the small molecule binding site; (ii) cleaved targets have preferred RNase L cleavage sites nearby small molecule binding sites; (iii) RiboTACs facilitate a cellular interaction between cleaved targets and RNase L; and (iv) the expression level of the target influences the extent of cleavage observed. In one example, we converted a binder of LGALS1 (galectin-1) mRNA into a RiboTAC. In MDA-MB-231 cells, the binder had no effect on galectin-1 protein levels, while the RiboTAC cleaved LGALS1 mRNA, reduced galectin-1 protein abundance, and affected galectin-1-associated oncogenic cellular phenotypes. Using LGALS1, we further assessed additional factors including the length of the linker that tethers the two components of the RiboTAC, cellular uptake, and the RNase L-recruiting module on RiboTAC potency. Collectively, these studies may facilitate triangulation of factors to enable the design of RiboTACs.

The RNA secondary structure of androgen receptor-FL and V7 transcripts reveals novel regulatory regions.

The androgen receptor (AR) is a ligand-dependent nuclear transcription factor belonging to the steroid hormone nuclear receptor family. Due to its roles in regulating cell proliferation and differentiation, AR is tightly regulated to maintain proper levels of itself and the many genes it controls. AR dysregulation is a driver of many human diseases including prostate cancer. Though this dysregulation often occurs at the RNA level, there are many unknowns surrounding post-transcriptional regulation of AR mRNA, particularly the role that RNA secondary structure plays. Thus, a comprehensive analysis of AR transcript secondary structure is needed. We address this through the computational and experimental analyses of two key isoforms, full length (AR-FL) and truncated (AR-V7). Here, a combination of in-cell RNA secondary structure probing experiments (targeted DMS-MaPseq) and computational predictions were used to characterize the static structural landscape and conformational dynamics of both isoforms. Additionally, in-cell assays were used to identify functionally relevant structures in the 5' and 3' UTRs of AR-FL. A notable example is a conserved stem loop structure in the 5'UTR of AR-FL that can bind to Poly(RC) Binding Protein 2 (PCBP2). Taken together, our results reveal novel features that regulate AR expression.

ScanFold 2.0: a rapid approach for identifying potential structured RNA targets in genomes and transcriptomes.

A major limiting factor in target discovery for both basic research and therapeutic intervention is the identification of structural and/or functional RNA elements in genomes and transcriptomes. This was the impetus for the original ScanFold algorithm, which provides maps of local RNA structural stability, evidence of sequence-ordered (potentially evolved) structure, and unique model structures comprised of recurring base pairs with the greatest structural bias. A key step in quantifying this propensity for ordered structure is the prediction of secondary structural stability for randomized sequences which, in the original implementation of ScanFold, is explicitly evaluated. This slow process has limited the rapid identification of ordered structures in large genomes/transcriptomes, which we seek to overcome in this current work introducing ScanFold 2.0. In this revised version of ScanFold, we no longer explicitly evaluate randomized sequence folding energy, but rather estimate it using a machine learning approach. For high randomization numbers, this can increase prediction speeds over 100-fold compared to ScanFold 1.0, allowing for the analysis of large sequences, as well as the use of additional folding algorithms that may be computationally expensive. In the testing of ScanFold 2.0, we re-evaluate the Zika, HIV, and SARS-CoV-2 genomes and compare both the consistency of results and the time of each run to ScanFold 1.0. We also re-evaluate the SARS-CoV-2 genome to assess the quality of ScanFold 2.0 predictions vs several biochemical structure probing datasets and compare the results to those of the original ScanFold program.

Thermodynamic and structural characterization of an EBV infected B-cell lymphoma transcriptome.

Epstein-Barr virus (EBV) is a widely prevalent human herpes virus infecting over 95% of all adults and is associated with a variety of B-cell cancers and induction of multiple sclerosis. EBV accomplishes this in part by expression of coding and noncoding RNAs and alteration of the host cell transcriptome. To better understand the structures which are forming in the viral and host transcriptomes of infected cells, the RNA structure probing technique Structure-seq2 was applied to the BJAB-B1 cell line (an EBV infected B-cell lymphoma). This resulted in reactivity profiles and secondary structural analyses for over 10000 human mRNAs and lncRNAs, along with 19 lytic and latent EBV transcripts. We report in-depth structural analyses for the human MYC mRNA and the human lncRNA CYTOR. Additionally, we provide a new model for the EBV noncoding RNA EBER2 and provide the first reported model for the EBV tandem terminal repeat RNA. In-depth thermodynamic and structural analyses were carried out with the motif discovery tool ScanFold and RNAfold prediction tool; subsequent covariation analyses were performed on resulting models finding various levels of support. ScanFold results for all analyzed transcripts are made available for viewing and download on the user-friendly RNAStructuromeDB.

Expansion of the RNAStructuromeDB to include secondary structural data spanning the human protein-coding transcriptome.

RNA plays vital functional roles in almost every component of biology, and these functional roles are often influenced by its folding into secondary and tertiary structures. An important role of RNA secondary structure is in maintaining proper gene regulation; therefore, making accurate predictions of the structures involved in these processes is important. In this study, we have expanded on our previous work that led to the creation of the RNAStructuromeDB. Unlike this previous study that analyzed the human genome at low resolution, we have now scanned the protein-coding human transcriptome at high (single nt) resolution. This provides more robust structure predictions for over 100,000 isoforms of known protein-coding genes. Notably, we also utilize the motif identification tool, ScanFold, to model structures with high propensity for ordered/evolved stability. All data have been uploaded to the RNAStructuromeDB, allowing for easy searching of transcripts, visualization of data tracks (via the Integrative Genomics Viewer or IGV), and download of ScanFold data-including unique highly-ordered motifs. Herein, we provide an example analysis of MAT2A to demonstrate the utility of ScanFold at finding known and novel secondary structures, highlighting regions of potential functionality, and guiding generation of functional hypotheses through use of the data.

Transcriptome-Wide Mapping of Small-Molecule RNA-Binding Sites in Cells Informs an Isoform-Specific Degrader of QSOX1 mRNA.

The interactions between cellular RNAs in MDA-MB-231 triple negative breast cancer cells and a panel of small molecules appended with a diazirine cross-linking moiety and an alkyne tag were probed transcriptome-wide in live cells. The alkyne tag allows for facile pull-down of cellular RNAs bound by each small molecule, and the enrichment of each RNA target defines the compound's molecular footprint. Among the 34 chemically diverse small molecules studied, six bound and enriched cellular RNAs. The most highly enriched interaction occurs between the novel RNA-binding compound F1 and a structured region in the 5' untranslated region of quiescin sulfhydryl oxidase 1 isoform a (QSOX1-a), not present in isoform b. Additional studies show that F1 specifically bound RNA over DNA and protein; that is, we studied the entire DNA, RNA, and protein interactome. This interaction was used to design a ribonuclease targeting chimera (RIBOTAC) to locally recruit Ribonuclease L to degrade QSOX1 mRNA in an isoform-specific manner, as QSOX1-a, but not QSOX1-b, mRNA and protein levels were reduced. The RIBOTAC alleviated QSOX1-mediated phenotypes in cancer cells. This approach can be broadly applied to discover ligands that bind RNA in cells, which could be bioactive themselves or augmented with functionality such as targeted degradation.

Analyses of human cancer driver genes uncovers evolutionarily conserved RNA structural elements involved in posttranscriptional control.

Experimental breakthroughs have provided unprecedented insights into the genes involved in cancer. The identification of such cancer driver genes is a major step in gaining a fuller understanding of oncogenesis and provides novel lists of potential therapeutic targets. A key area that requires additional study is the posttranscriptional control mechanisms at work in cancer driver genes. This is important not only for basic insights into the biology of cancer, but also to advance new therapeutic modalities that target RNA-an emerging field with great promise toward the treatment of various cancers. In the current study we performed an in silico analysis on the transcripts associated with 800 cancer driver genes (10,390 unique transcripts) that identified 179,190 secondary structural motifs with evidence of evolutionarily ordered structures with unusual thermodynamic stability. Narrowing to one transcript per gene, 35,426 predicted structures were subjected to phylogenetic comparisons of sequence and structural conservation. This identified 7,001 RNA secondary structures embedded in transcripts with evidence of covariation between paired sites, supporting structure models and suggesting functional significance. A select set of seven structures were tested in vitro for their ability to regulate gene expression; all were found to have significant effects. These results indicate potentially widespread roles for RNA structure in posttranscriptional control of human cancer driver genes.

A survey of RNA secondary structural propensity encoded within human herpesvirus genomes: global comparisons and local motifs.

There are nine herpesviruses known to infect humans, of which Epstein-Barr virus (EBV) is the most widely distributed (>90% of adults infected). This ubiquitous virus is implicated in a variety of cancers and autoimmune diseases. Previous analyses of the EBV genome revealed numerous regions with evidence of generating unusually stable and conserved RNA secondary structures and led to the discovery of a novel class of EBV non-coding (nc)RNAs: the stable intronic sequence (sis)RNAs. To gain a better understanding of the roles of RNA structure in EBV biology and pathogenicity, we revisit EBV using recently developed tools for genome-wide motif discovery and RNA structural characterization. This corroborated previous results and revealed novel motifs with potential functionality; one of which has been experimentally validated. Additionally, since many herpesviruses increasingly rival the seroprevalence of EBV (VZV, HHV-6 and HHV-7 being the most notable), analyses were expanded to include all sequenced human Herpesvirus RefSeq genomes, allowing for genomic comparisons. In total 10 genomes were analyzed, for EBV (types 1 and 2), HCMV, HHV-6A, HHV-6B, HHV-7, HSV-1, HSV-2, KSHV, and VZV. All resulting data were archived in the RNAStructuromeDB (https://structurome.bb.iastate.edu/herpesvirus) to make them available to a wide array of researchers.

Plasmodium falciparum translational machinery condones polyadenosine repeats.

Plasmodium falciparum is a causative agent of human malaria. Sixty percent of mRNAs from its extremely AT-rich (81%) genome harbor long polyadenosine (polyA) runs within their ORFs, distinguishing the parasite from its hosts and other sequenced organisms. Recent studies indicate polyA runs cause ribosome stalling and frameshifting, triggering mRNA surveillance pathways and attenuating protein synthesis. Here, we show that P. falciparum is an exception to this rule. We demonstrate that both endogenous genes and reporter sequences containing long polyA runs are efficiently and accurately translated in P. falciparum cells. We show that polyA runs do not elicit any response from No Go Decay (NGD) or result in the production of frameshifted proteins. This is in stark contrast to what we observe in human cells or T. thermophila, an organism with similar AT-content. Finally, using stalling reporters we show that Plasmodium cells evolved not to have a fully functional NGD pathway.

Differential expression analysis comparing EBV uninfected to infected human cell lines identifies induced non-micro small non-coding RNAs.

Epstein-Barr virus (EBV) is a ubiquitous human herpes virus, which is implicated in cancer and various autoimmune diseases. This study profiles non-micro small non-coding RNA expression changes induced by latent EBV infection. Using small RNA-Seq, 346 non-micro small RNAs were identified as being significantly differentially expressed between EBV(+) BJAB-B1 and EBV(-) BJAB cell lines. Select small RNA expression changes were experimentally validated in the BJAB-B1 cell line as well as the EBV-infected Raji and Jijoye cell lines. This latter analysis recapitulated the previously identified induction of vault RNA1, while also finding novel evidence for the deregulation of several tRNAs and a snoRNA.

Translation of the intrinsically disordered protein α-synuclein is inhibited by a small molecule targeting its structured mRNA.

Many proteins are refractory to targeting because they lack small-molecule binding pockets. An alternative to drugging these proteins directly is to target the messenger (m)RNA that encodes them, thereby reducing protein levels. We describe such an approach for the difficult-to-target protein α-synuclein encoded by the SNCA gene. Multiplication of the SNCA gene locus causes dominantly inherited Parkinson's disease (PD), and α-synuclein protein aggregates in Lewy bodies and Lewy neurites in sporadic PD. Thus, reducing the expression of α-synuclein protein is expected to have therapeutic value. Fortuitously, the SNCA mRNA has a structured iron-responsive element (IRE) in its 5' untranslated region (5' UTR) that controls its translation. Using sequence-based design, we discovered small molecules that target the IRE structure and inhibit SNCA translation in cells, the most potent of which is named Synucleozid. Both in vitro and cellular profiling studies showed Synucleozid directly targets the α-synuclein mRNA 5' UTR at the designed site. Mechanistic studies revealed that Synucleozid reduces α-synuclein protein levels by decreasing the amount of SNCA mRNA loaded into polysomes, mechanistically providing a cytoprotective effect in cells. Proteome- and transcriptome-wide studies showed that the compound's selectivity makes Synucleozid suitable for further development. Importantly, transcriptome-wide analysis of mRNAs that encode intrinsically disordered proteins revealed that each has structured regions that could be targeted with small molecules. These findings demonstrate the potential for targeting undruggable proteins at the level of their coding mRNAs. This approach, as applied to SNCA, is a promising disease-modifying therapeutic strategy for PD and other α-synucleinopathies.

RNA structural analysis of the MYC mRNA reveals conserved motifs that affect gene expression.

The MYC gene encodes a human transcription factor and proto-oncogene that is dysregulated in over half of all known cancers. To better understand potential post-transcriptional regulatory features affecting MYC expression, we analyzed secondary structures in the MYC mRNA using a program that is optimized for finding small locally-folded motifs with a high propensity for function. This was accomplished by calculating folding metrics across the MYC sequence using a sliding analysis window and generating unique consensus base pairing models weighted by their lower-than-random predicted folding energy. A series of 30 motifs were identified, primarily in the 5' and 3' untranslated regions, which show evidence of structural conservation and compensating mutations across vertebrate MYC homologs. This analysis was able to recapitulate known elements found within an internal ribosomal entry site, as well as discover a novel element in the 3' UTR that is unusually stable and conserved. This novel motif was shown to affect MYC expression, potentially via the modulation of miRNA target accessibility or other trans-regulatory factors. In addition to providing basic insights into mechanisms that regulate MYC expression, this study provides numerous, potentially druggable RNA targets for the MYC gene, which is considered "undruggable" at the protein level.

Analysis of a structured intronic region of the LMP2 pre-mRNA from EBV reveals associations with human regulatory proteins and nuclear actin.

OBJECTIVE: The pre-mRNA of the Epstein-Barr virus LMP2 (latent membrane protein 2) has a region of unusual RNA structure that partially spans two consecutive exons and the entire intervening intron; suggesting RNA folding might affect splicing-particularly via interactions with human regulatory proteins. To better understand the roles of protein associations with this structured intronic region, we undertook a combined bioinformatics (motif searching) and experimental analysis (biotin pulldowns and RNA immunoprecipitations) of protein binding. RESULT: Characterization of the ribonucleoprotein composition of this region revealed several human proteins as interactors: regulatory proteins hnRNP A1 (heterogeneous nuclear ribonucleoprotein A1), hnRNP U, HuR (human antigen R), and PSF (polypyrimidine tract-binding protein-associated splicing factor), as well as, unexpectedly, the cytoskeletal protein actin. Treatment of EBV-positive cells with drugs that alter actin polymerization specifically showed marked effects on splicing in this region. This suggests a potentially novel role for nuclear actin in regulation of viral RNA splicing.

Comparative Analysis of Gammaherpesvirus Circular RNA Repertoires: Conserved and Unique Viral Circular RNAs.

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.

The Epstein Barr virus circRNAome.

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.

Host miRNA degradation by Herpesvirus saimiri small nuclear RNA requires an unstructured interacting region.

Herpesvirus saimiri, an oncogenic herpesvirus, during latency produces seven small nuclear RNAs, called the Herpesvirus saimiri U RNAs (HSUR1-7). HSUR1 mediates degradation of the host microRNA, miR-27, via a process that requires imperfect base-pairing. The decreased levels of miR-27 lead to prolonged T-cell activation and likely contribute to oncogenesis. To gain insight into HSUR1-mediated degradation of miR-27, we probed the in vivo secondary structure of HSUR1 and coupled this with bioinformatic structural analyses. The results suggest that HSUR1 adopts a conformation different than previously believed and that the region complementary to miR-27 lacks stable structure. To determine whether HSUR1 structural flexibility is important for its ability to mediate miR-27 degradation, we performed structurally informative mutagenic analyses of HSUR1. HSUR1 mutants in which the miR-27 binding site sequence is preserved, but sequestered in predicted helices, lose their ability to decrease miR-27 levels. These results indicate that the HSUR1 miR27-binding region must be available in a conformationally flexible segment for noncoding RNA function.

In silico discovery and modeling of non-coding RNA structure in viruses.

This review covers several computational methods for discovering structured non-coding RNAs in viruses and modeling their putative secondary structures. Here we will use examples from two target viruses to highlight these approaches: influenza A virus-a relatively small, segmented RNA virus; and Epstein-Barr virus-a relatively large DNA virus with a complex transcriptome. Each system has unique challenges to overcome and unique characteristics to exploit. From these particular cases, generically useful approaches can be derived for the study of additional viral targets.

EBV noncoding RNA binds nascent RNA to drive host PAX5 to viral DNA.

EBER2 is an abundant nuclear noncoding RNA expressed by the Epstein-Barr virus (EBV). Probing its possible chromatin localization by CHART revealed EBER2's presence at the terminal repeats (TRs) of the latent EBV genome, overlapping previously identified binding sites for the B cell transcription factor PAX5. EBER2 interacts with PAX5 and is required for the localization of PAX5 to the TRs. EBER2 knockdown phenocopies PAX5 depletion in upregulating the expression of LMP2A/B and LMP1, genes nearest the TRs. Knockdown of EBER2 also decreases EBV lytic replication, underscoring the essential role of the TRs in viral replication. Recruitment of the EBER2-PAX5 complex is mediated by base-pairing between EBER2 and nascent transcripts from the TR locus. The interaction is evolutionarily conserved in the related primate herpesvirus CeHV15 despite great sequence divergence. Using base-pairing with nascent RNA to guide an interacting transcription factor to its DNA target site is a previously undescribed function for a trans-acting noncoding RNA.

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.

RNA families in Epstein-Barr virus.

Epstein-Barr virus (EBV) is a tumorigenic human γ-herpesvirus, which produces several known structured RNAs with functional importance: two are implicated in latency maintenance and tumorigenic phenotypes, EBER1 and EBER2; a viral small nucleolar RNA (v-snoRNA1) that may generate a small regulatory RNA; and an internal ribosomal entry site in the EBNA1 mRNA. A recent bioinformatics and RNA-Seq study of EBV identified two novel EBV non-coding (nc)RNAs with evolutionary conservation in lymphocryptoviruses and likely functional importance. Both RNAs are transcribed from a repetitive region of the EBV genome (the W repeats) during a highly oncogenic type of viral latency. One novel ncRNA can form a massive (586 nt) hairpin, while the other RNA is generated from a short (81 nt) intron and is found in high abundance in EBV-infected cells.