High-density resolution of the Kaposi's sarcoma associated herpesvirus transcriptome identifies novel transcript isoforms generated by long-range transcription and alternative splicing.

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.

Viral reprogramming of host transcription initiation.

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.

SpliceTools, a suite of downstream RNA splicing analysis tools to investigate mechanisms and impact of alternative splicing.

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.

Alterations in gene expression and microbiome composition upon calcium-sensing receptor deletion in the mouse esophagus.

The calcium-sensing receptor (CaSR), a G protein-coupled receptor, regulates Ca2+ concentration in plasma by regulating parathyroid hormone secretion. In other tissues, it is reported to play roles in cellular differentiation and migration and in secretion and absorption. We reported previously that CaSR can be conditionally deleted in the mouse esophagus. This conditional knockout (KO) (EsoCaSR-/-) model showed a significant reduction in the levels of adherens and tight junction proteins and had a marked buildup of bacteria on the luminal esophageal surface. To further examine the role of CaSR, we used RNA sequencing to determine gene expression profiles in esophageal epithelia of control and EsoCaSR-/-mice RNA Seq data indicated upregulation of gene sets involved in DNA replication and cell cycle in EsoCaSR-/-. This is accompanied by the downregulation of gene sets involved in the innate immune response and protein homeostasis including peptide elongation and protein trafficking. Ingenuity pathway analysis (IPA) demonstrated that these genes are mapped to important biological networks including calcium and Ras homologus A (RhoA) signaling pathways. To further explore the bacterial buildup in EsoCaSR-/- esophageal tissue, 16S sequencing of the mucosal-associated bacterial microbiome was performed. Three bacterial species, g_Rodentibacter, s_Rodentibacter_unclassified, and s_Lactobacillus_hilgardi were significantly increased in EsoCaSR-/-. Furthermore, metagenomic analysis of 16S sequences indicated that pathways related to oxidative phosphorylation and metabolism were downregulated in EsoCaSR-/- tissues. These data demonstrate that CaSR impacts major pathways of cell proliferation, differentiation, cell cycle, and innate immune response in esophageal epithelium. The disruption of these pathways causes inflammation and significant modifications of the microbiome.NEW & NOTEWORTHY Calcium-sensing receptor (CaSR) plays a significant role in maintaining the barrier function of esophageal epithelium. Using RNA sequencing, we show that conditional deletion of CaSR from mouse esophagus causes upregulation of genes involved in DNA replication and cell cycle and downregulation of genes involved in the innate immune response, protein translation, and cellular protein synthesis. Pathway analysis shows disruption of signaling pathways of calcium and actin cytoskeleton. These changes caused inflammation and esophageal dysbiosis.

SON and Its Alternatively Spliced Isoforms Control MLL Complex-Mediated H3K4me3 and Transcription of Leukemia-Associated Genes.

Dysregulation of MLL complex-mediated histone methylation plays a pivotal role in gene expression associated with diseases, but little is known about cellular factors modulating MLL complex activity. Here, we report that SON, previously known as an RNA splicing factor, controls MLL complex-mediated transcriptional initiation. SON binds to DNA near transcription start sites, interacts with menin, and inhibits MLL complex assembly, resulting in decreased H3K4me3 and transcriptional repression. Importantly, alternatively spliced short isoforms of SON are markedly upregulated in acute myeloid leukemia. The short isoforms compete with full-length SON for chromatin occupancy but lack the menin-binding ability, thereby antagonizing full-length SON function in transcriptional repression while not impairing full-length SON-mediated RNA splicing. Furthermore, overexpression of a short isoform of SON enhances replating potential of hematopoietic progenitors. Our findings define SON as a fine-tuner of the MLL-menin interaction and reveal short SON overexpression as a marker indicating aberrant transcriptional initiation in leukemia.

EBV miRNAs are potent effectors of tumor cell transcriptome remodeling in promoting immune escape.

The Epstein Barr virus (EBV) contributes to the tumor phenotype through a limited set of primarily non-coding viral RNAs, including 31 mature miRNAs. Here we investigated the impact of EBV miRNAs on remodeling the tumor cell transcriptome. Strikingly, EBV miRNAs displayed exceptionally abundant expression in primary EBV-associated Burkitt's Lymphomas (BLs) and Gastric Carcinomas (GCs). To investigate viral miRNA targeting, we used the high-resolution approach, CLASH in GC and BL cell models. Affinity constant calculations of targeting efficacies for CLASH hits showed that viral miRNAs bind their targets more effectively than their host counterparts, as did Kaposi's sarcoma-associated herpesvirus (KSHV) and murine gammaherpesvirus 68 (MHV68) miRNAs. Using public BL and GC RNA-seq datasets, we found that high EBV miRNA targeting efficacies translates to enhanced reduction of target expression. Pathway analysis of high efficacy EBV miRNA targets showed enrichment for innate and adaptive immune responses. Inhibition of the immune response by EBV miRNAs was functionally validated in vivo through the finding of inverse correlations between EBV miRNAs and immune cell infiltration and T-cell diversity in BL and GC datasets. Together, this study demonstrates that EBV miRNAs are potent effectors of the tumor transcriptome that play a role in suppressing host immune response.

Profiling the activity of the para-caspase MALT1 in B-cell acute lymphoblastic leukemia for potential targeted therapeutic application.

B cell acute lymphoblastic leukemia (B-ALL) remains a hard-to-treat disease with a poor prognosis in adults. Mucosa-associated lymphoid tissue lymphoma translocation protein 1 (MALT1) is a para-caspase required for B-cell receptor (BCR)-mediated NF-κB activation. Inhibition of MALT1 in preclinical models has proven efficacious in many B-cell malignancies including chronic lymphocytic leukemia, mantle cell lymphoma and diffuse large B-cell lymphoma. We sought to examine the role of MALT1 in B-ALL and determine the biological consequences of its inhibition. Targeting MALT1 with both Z-VRPR-fmk and MI-2 efficiently kills B-ALL cells independent of the cell-of-origin (pro, pre, mature) or the presence of the Philadelphia chromosome, and spares normal B-cells. The mechanism of cell death was through apoptotic induction, mostly in cycling cells. The proteolytic activity of MALT1 can be studied by measuring its ability to cleave its substrates. Surprisingly, with the exception of mature B-ALL, we did not detect cleavage of MALT1 substrates at baseline, nor after proteasomal inhibition or following activation of pre-BCR. To explore the possibility of a distinct role for MALT1 in B-ALL, independent of signaling through BCR, we studied the changes in gene expression profiling following a 24-hour treatment with MI-2 in 12 B-ALL cell lines. Our transcriptome analysis revealed a strong inhibitory effect on MYC-regulated gene signatures, further confirmed by Myc protein downregulation, concomitant with an increase in the Myc degrader FBXW7. In conclusion, our evidence suggests a novel role for MALT1 in B-ALL through Myc regulation and provides support for clinical testing of MALT1 inhibitors in B-ALL.

Transcriptional signatures of Zika virus infection in astrocytes.

Astrocytes are an early and important target of Zika virus (ZIKV) infection in the developing brain, but the impacts of infection on astrocyte function remain controversial. Given that nonhuman primate (NHP) models of ZIKV infection replicate aspects of neurologic disease seen in human infections, we cultured primary astrocytes from the brain tissue of infant rhesus macaques and then infected the cells with Asian or African lineage ZIKV to identify transcriptional patterns associated with infection in these cells. The African lineage virus appeared to have greater infectivity and promote stronger antiviral signaling, but infection by either strain ultimately produced typical virus response patterns. Both viruses induced hypoxic stress, but the Asian lineage strain additionally had an effect on metabolic and lipid biosynthesis pathways. Together, these findings describe an NHP astrocyte model that may be used to assess transcriptional signatures following ZIKV infection.

Screen technical noise in single cell RNA sequencing data.

We proposed a data cleaning pipeline for single cell (SC) RNA-seq data, where we first screen genes (gene-wise screening) followed by screening cell libraries (library-wise screening). Gene-wise screening is based on the expectation that for a gene with a low technical noise, a gene's count in a library will tend to increase with the increase of library size, which was tested using negative binomial regression of gene count (as dependent variable) against library size (as independent variable). Library-wise screening is based on the expectation that across-library correlations for housekeeping (HK) genes is expected to be higher than the correlations for non-housekeeping (NHK) genes in those libraries with low technical noise. We removed those libraries, whose mean pairwise correlation for HK genes is NOT significantly higher than that for NHK genes. We successfully applied the pipeline to two large SC RNA-seq datasets. The pipeline was also developed into an R package.

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.

Assessment of viral RNA in idiopathic pulmonary fibrosis using RNA-seq.

BACKGROUND: Numerous publications suggest an association between herpes virus infection and idiopathic pulmonary fibrosis (IPF). These reports have employed immunohistochemistry, in situ hybridization and/or PCR, which are susceptible to specificity artifacts. METHODS: We investigated the possible association between IPF and viral RNA expression using next-generation sequencing, which has the potential to provide a high degree of both sensitivity and specificity. We quantified viral RNA expression for 740 viruses in 28 IPF patient lung biopsy samples and 20 controls. Key RNA-seq results were confirmed using Real-time RT-PCR for select viruses (EBV, HCV, herpesvirus saimiri and HERV-K). RESULTS: We identified sporadic low-level evidence of viral infections in our lung tissue specimens, but did not find a statistical difference for expression of any virus, including EBV, herpesvirus saimiri and HERV-K, between IPF and control lungs. CONCLUSIONS: To the best of our knowledge, this is the first publication that employs RNA-seq to assess whether viral infections are linked to the pathogenesis of IPF. Our results do not address the role of viral infection in acute exacerbations of IPF, however, this analysis patently did not support an association between herpes virus detection and IPF.

Maintenance of AP-2-Dependent Functional Activities of Nef Restricts Pathways of Immune Escape from CD8 T Lymphocyte Responses.

Nef-specific CD8+ T lymphocytes (CD8TL) are linked to extraordinary control of primate lentiviral replication, but the mechanisms underlying their efficacy remain largely unknown. The immunodominant, Mamu-B*017:01+-restricted Nef195-203MW9 epitope in SIVmac239 partially overlaps a sorting motif important for interactions with host AP-2 proteins and, hence, downmodulation of several host proteins, including Tetherin (CD317/BST-2), CD28, CD4, SERINC3, and SERINC5. We reasoned that CD8TL-driven evolution in this epitope might compromise Nef's ability to modulate these important molecules. Here, we used deep sequencing of SIV from nine B*017:01+ macaques throughout infection with SIVmac239 to characterize the patterns of viral escape in this epitope and then assayed the impacts of these variants on Nef-mediated modulation of multiple host molecules. Acute variation in multiple Nef195-203MW9 residues significantly compromised Nef's ability to downregulate surface Tetherin, CD4, and CD28 and reduced its ability to prevent SERINC5-mediated reduction in viral infectivity but did not impact downregulation of CD3 or major histocompatibility complex class I, suggesting the selective disruption of immunomodulatory pathways involving Nef AP-2 interactions. Together, our data illuminate a pattern of viral escape dictated by a selective balance to maintain AP-2-mediated downregulation while evading epitope-specific CD8TL responses. These data could shed light on mechanisms of both CD8TL-driven viral control generally and on Mamu-B*017:01-mediated viral control specifically.IMPORTANCE A rare subset of humans infected with HIV-1 and macaques infected with SIV can control the virus without aid of antiviral medications. A common feature of these individuals is the ability to mount unusually effective CD8 T lymphocyte responses against the virus. One of the most formidable aspects of HIV is its ability to evolve to evade immune responses, particularly CD8 T lymphocytes. We show that macaques that target a specific peptide in the SIV Nef protein are capable of better control of the virus and that, as the virus evolves to escape this response, it does so at a cost to specific functions performed by the Nef protein. Our results help show how the virus can be controlled by an immune response, which could help in designing effective vaccines.

Driver gene mutations based clustering of tumors: methods and applications.

Motivation: Somatic mutations in proto-oncogenes and tumor suppressor genes constitute a major category of causal genetic abnormalities in tumor cells. The mutation spectra of thousands of tumors have been generated by The Cancer Genome Atlas (TCGA) and other whole genome (exome) sequencing projects. A promising approach to utilizing these resources for precision medicine is to identify genetic similarity-based sub-types within a cancer type and relate the pinpointed sub-types to the clinical outcomes and pathologic characteristics of patients. Results: We propose two novel methods, ccpwModel and xGeneModel, for mutation-based clustering of tumors. In the former, binary variables indicating the status of cancer driver genes in tumors and the genes' involvement in the core cancer pathways are treated as the features in the clustering process. In the latter, the functional similarities of putative cancer driver genes and their confidence scores as the 'true' driver genes are integrated with the mutation spectra to calculate the genetic distances between tumors. We apply both methods to the TCGA data of 16 cancer types. Promising results are obtained when these methods are compared to state-of-the-art approaches as to the associations between the determined tumor clusters and patient race (or survival time). We further extend the analysis to detect mutation-characterized transcriptomic prognostic signatures, which are directly relevant to the etiology of carcinogenesis. Availability and implementation: R codes and example data for ccpwModel and xGeneModel can be obtained from http://webusers.xula.edu/kzhang/ISMB2018/ccpw_xGene_software.zip. Supplementary information: Supplementary data are available at Bioinformatics online.

A comprehensive approach to expression of L1 loci.

L1 elements represent the only currently active, autonomous retrotransposon in the human genome, and they make major contributions to human genetic instability. The vast majority of the 500 000 L1 elements in the genome are defective, and only a relatively few can contribute to the retrotransposition process. However, there is currently no comprehensive approach to identify the specific loci that are actively transcribed separate from the excess of L1-related sequences that are co-transcribed within genes. We have developed RNA-Seq procedures, as well as a 1200 bp 5΄ RACE product coupled with PACBio sequencing that can identify the specific L1 loci that contribute most of the L1-related RNA reads. At least 99% of L1-related sequences found in RNA do not arise from the L1 promoter, instead representing pieces of L1 incorporated in other cellular RNAs. In any given cell type a relatively few active L1 loci contribute to the 'authentic' L1 transcripts that arise from the L1 promoter, with significantly different loci seen expressed in different tissues.

Transferring knowledge of bacterial protein interaction networks to predict pathogen targeted human genes and immune signaling pathways: a case study on M. tuberculosis.

BACKGROUND: Bacterial invasive infection and host immune response is fundamental to the understanding of pathogen pathogenesis and the discovery of effective therapeutic drugs. However, there are very few experimental studies on the signaling cross-talks between bacteria and human host to date. METHODS: In this work, taking M. tuberculosis H37Rv (MTB) that is co-evolving with its human host as an example, we propose a general computational framework that exploits the known bacterial pathogen protein interaction networks in STRING database to predict pathogen-host protein interactions and their signaling cross-talks. In this framework, significant interlogs are derived from the known pathogen protein interaction networks to train a predictive l2-regularized logistic regression model. RESULTS: The computational results show that the proposed method achieves excellent performance of cross validation as well as low predicted positive rates on the less significant interlogs and non-interlogs, indicating a low risk of false discovery. We further conduct gene ontology (GO) and pathway enrichment analyses of the predicted pathogen-host protein interaction networks, which potentially provides insights into the machinery that M. tuberculosis H37Rv targets human genes and signaling pathways. In addition, we analyse the pathogen-host protein interactions related to drug resistance, inhibition of which potentially provides an alternative solution to M. tuberculosis H37Rv drug resistance. CONCLUSIONS: The proposed machine learning framework has been verified effective for predicting bacteria-host protein interactions via known bacterial protein interaction networks. For a vast majority of bacterial pathogens that lacks experimental studies of bacteria-host protein interactions, this framework is supposed to achieve a general-purpose applicability. The predicted protein interaction networks between M. tuberculosis H37Rv and Homo sapiens, provided in the Additional files, promise to gain applications in the two fields: (1) providing an alternative solution to drug resistance; (2) revealing the patterns that M. tuberculosis H37Rv genes target human immune signaling pathways.

Global transcript structure resolution of high gene density genomes through multi-platform data integration.

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.

Secreted Oral Epithelial Cell Membrane Vesicles Induce Epstein-Barr Virus Reactivation in Latently Infected B Cells.

UNLABELLED: In the oral epithelium, peripheral stores of Epstein-Barr virus (EBV) are transmitted from infiltrating B cells to epithelial cells. Once the virus is transmitted to epithelial cells, the highly permissive nature of this cell type for lytic replication allows virus amplification and exchange to other hosts. Since the initial transfer of EBV from B cells to epithelial cells requires transitioning of the B-cell to a state that induces virus reactivation, we hypothesized that there might be epithelium-specific signals that allow the infiltrating B cells to sense the appropriate environment to initiate reactivation and begin this exchange process. We previously found that the epithelium-specific miR-200 family of microRNAs promotes EBV lytic replication. Here we show that there are high levels of miR-200 family members in oral and tonsillar epithelia and in saliva. Analysis of cultured oral epithelial cells (OKF6) showed that they actively secrete membrane vesicles (exosomes) that are enriched with miR-200 family members. Coculturing of EBV-positive B cells with OKF6 cells induced viral reactivation. Further, treatment of EBV-positive B cells with OKF6 cell-derived membrane vesicles promoted reactivation. Using a cell system that does not naturally express miR-200 family members, we found that enforced expression of a miR-200 family member produced membrane vesicles that were able to induce the lytic cascade in EBV-positive B cells. We propose that membrane vesicles secreted by oral and tonsillar epithelial cells may serve as a tissue-specific environmental cue that initiates reactivation in B cells, promoting the transfer of virus from peripheral B-cell stores to the oral epithelium to facilitate virus amplification and exchange to other hosts. IMPORTANCE: Epstein-Barr virus (EBV) is an important human pathogen that is causally associated with several lymphomas and carcinomas. The switch from latency to the lytic cycle is critical for successful host infection and for EBV pathogenesis. Although the EBV lytic cycle can be triggered by certain agents in vitro, the mechanisms that signal reactivation in vivo are poorly understood. We previously reported that endogenously expressed miR-200 family members likely play a role in facilitating the lytic tendencies of EBV in epithelial cells. Here we show that membrane vesicles secreted from oral epithelial cells contain miR-200 family members and that they can be transmitted to proximal EBV-positive B cells, where they trigger reactivation. We propose that this intercellular communication pathway may serve as a sensor mechanism for infiltrating B cells to recognize an appropriate environment to initiate reactivation, thereby allowing the exchange of virus to the oral epithelium.

Arsenic trioxide inhibits EBV reactivation and promotes cell death in EBV-positive lymphoma cells.

BACKGROUND: Epstein-Barr Virus (EBV) is associated with hematopoietic malignancies, such as Burkitt's lymphoma, post-transplantation lymphoproliferative disorder, and diffuse large B-cell lymphoma. The current approach for EBV-associated lymphoma involves chemotherapy to eradicate cancer cells, however, normal cells may be injured and organ dysfunction may occur with currently employed regimens. This research is focused on employing arsenic trioxide (ATO) as EBV-specific cancer therapy takes advantage of the fact the EBV resides within the malignant cells. METHODS AND RESULTS: Our research reveals that low ATO inhibits EBV gene expression and genome replication. EBV spontaneous reactivation starts as early as 6 h after re-suspending EBV-positive Mutu cells in RPMI media in the absence of ATO, however this does not occur in Mutu cells cultured with ATO. ATO's inhibition of EBV spontaneous reactivation is dose dependent. The expression of the EBV immediate early gene Zta and early gene BMRF1 is blocked with low concentrations of ATO (0.5 nM - 2 nM) in EBV latency type I cells and EBV-infected PBMC cells. The combination of ATO and ganciclovir further diminishes EBV gene expression. ATO-mediated reduction of EBV gene expression can be rescued by co-treatment with the proteasome inhibitor MG132, indicating that ATO promotes ubiquitin conjugation and proteasomal degradation of EBV genes. Co-immunoprecipitation assays with antibodies against Zta pulls down more ubiquitin in ATO treated cell lysates. Furthermore, MG132 reverses the inhibitory effect of ATO on anti-IgM-, PMA- and TGF-ß-mediated EBV reactivation. Thus, mechanistically ATO's inhibition of EBV gene expression occurs via the ubiquitin pathway. Moreover, ATO treatment results in increased cell death in EBV-positive cells compared to EBV-negative cells, as demonstrated by both MTT and trypan blue assays. ATO-induced cell death in EBV-positive cells is dose dependent. ATO and ganciclovir in combination further enhances cell death specifically in EBV-positive cells. CONCLUSION: ATO-mediated inhibition of EBV lytic gene expression results in cell death selectively in EBV-positive lymphocytes, suggesting that ATO may potentially serve as a drug to treat EBV-related lymphomas in the clinical setting.

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.