RNA families in Epstein-Barr virus

Intron-capture RNA-seq reveals the landscape of intronic RNAs in Arabidopsis

Intronic RNAs have been overlooked for a long time: They are functional, but treated as "junk." In this work, we designed a new sequencing strategy to investigate intronic RNAs. By using intron-capture RNA-seq, we systematically analyzed the intronic RNAs in Arabidopsis by zooming into the intronic regions an order of magnitude deeper than in previous work. Our key findings include: (1) Intron-capture RNA-seq is a much more efficient approach to analyze intronic RNAs than total RNA-seq and mRNA-seq. (2) We identified three types of intronic RNAs, and found that the GC pattern differs significantly between the introns with and without intronic RNAs. (3) We detected many hidden elements in introns, including circular RNAs, splice junctions, and transcripts that have previously been overlooked. (4) The expression of these intronic RNAs varies during the time course of pathogen infection, which indicates that an unknown mechanism may exist for these RNAs. (5) We also demonstrated that most of intronic RNAs are detectable in both Arabidopsis and rice, suggesting that these non-coding molecules are conserved. Taken together, this work proposes an efficient strategy to analyze intronic RNAs, and provides an unprecedented view of this essential component in biological pathways.

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.

Epstein-Barr virus: Biology and clinical disease

Epstein-Barr virus (EBV) is a ubiquitous, oncogenic virus that is associated with a number of different human malignancies as well as autoimmune disorders. The expression of EBV viral proteins and non-coding RNAs contribute to EBV-mediated disease pathologies. The virus establishes life-long latency in the human host and is adept at evading host innate and adaptive immune responses. In this review, we discuss the life cycle of EBV, the various functions of EBV-encoded proteins and RNAs, the ability of the virus to activate and evade immune responses, as well as the neoplastic and autoimmune diseases that are associated with EBV infection in the human population.

Mono a Mano: ZBP1’s Love–Hate Relationship with the Kissing Virus

Z-DNA binding protein (ZBP1) very much represents the nuclear option. By initiating inflammatory cell death (ICD), ZBP1 activates host defenses to destroy infectious threats. ZBP1 is also able to induce noninflammatory regulated cell death via apoptosis (RCD). ZBP1 senses the presence of left-handed Z-DNA and Z-RNA (ZNA), including that formed by expression of endogenous retroelements. Viruses such as the Epstein–Barr “kissing virus” inhibit ICD, RCD and other cell death signaling pathways to produce persistent infection. EBV undergoes lytic replication in plasma cells, which maintain detectable levels of basal ZBP1 expression, leading us to suggest a new role for ZBP1 in maintaining EBV latency, one of benefit for both host and virus. We provide an overview of the pathways that are involved in establishing latent infection, including those regulated by MYC and NF-κB. We describe and provide a synthesis of the evidence supporting a role for ZNA in these pathways, highlighting the positive and negative selection of ZNA forming sequences in the EBV genome that underscores the coadaptation of host and virus. Instead of a fight to the death, a state of détente now exists where persistent infection by the virus is tolerated by the host, while disease outcomes such as death, autoimmunity and cancer are minimized. Based on these new insights, we propose actionable therapeutic approaches to unhost EBV.

Epstein-Barr Virus and Systemic Autoimmune Diseases

Epstein-Barr Virus (EBV) is an extremely successful human herpes virus, which infects essentially all human beings at some time during their life span. EBV infection and the associated immune response results in production of antibodies (seroconversion), which occurs mainly during the first years of life, but may also happen during adolescence or later in life. Infection of adolescents can result in infectious mononucleosis, an acute serious condition characterized by massive lymphocytosis. Transmission of EBV mainly occurs through saliva but can rarely be spread through semen or blood, e.g. through organ transplantations and blood transfusions. EBV transmission through oral secretions results in infection of epithelial cells of the oropharynx. From the epithelial cells EBV can infect B cells, which are the major reservoir for the virus, but other cell types may also become infected. As a result, EBV can shuttle between different cell types, mainly B cells and epithelial cells. Moreover, since the virus can switch between a latent and a lytic life cycle, EBV has the ability to cause chronic relapsing/reactivating infections. Chronic or recurrent EBV infection of epithelial cells has been linked to systemic lupus erythematosus and Sjögren’s syndrome, whereas chronic/recurrent infection of B cells has been associated with rheumatoid arthritis, multiple sclerosis and other diseases. Accordingly, since EBV can shuttle between epithelial cells and B cells, the systemic autoimmune diseases often occur as overlapping syndromes with symptoms and characteristic autoantibodies (e.g. antinuclear antibodies and rheumatoid factors) reflecting epithelial and/or B cell infection.

Epstein-Barr Virus and Multiple Sclerosis

Multiple sclerosis (MS) is a neurologic disease affecting myelinated nerves in the central nervous system (CNS). The disease often debuts as a clinically isolated syndrome, e.g., optic neuritis (ON), which later develops into relapsing-remitting (RR) MS, with temporal attacks or primary progressive (PP) MS. Characteristic features of MS are inflammatory foci in the CNS and intrathecal synthesis of immunoglobulins (Igs), measured as an IgG index, oligoclonal bands (OCBs), or specific antibody indexes. Major predisposing factors for MS are certain tissue types (e.g., HLA DRB1*15:01), vitamin D deficiency, smoking, obesity, and infection with Epstein-Barr virus (EBV). Many of the clinical signs of MS described above can be explained by chronic/recurrent EBV infection and current models of EBV involvement suggest that RRMS may be caused by repeated entry of EBV-transformed B cells to the CNS in connection with attacks, while PPMS may be caused by more chronic activity of EBV-transformed B cells in the CNS. In line with the model of EBV's role in MS, new treatments based on monoclonal antibodies (MAbs) targeting B cells have shown good efficacy in clinical trials both for RRMS and PPMS, while MAbs inhibiting B cell mobilization and entry to the CNS have shown efficacy in RRMS. Thus, these agents, which are now first line therapy in many patients, may be hypothesized to function by counteracting a chronic EBV infection.

The roles of miRNAs and lncRNAs in Epstein-Barr virus associated epithelial cell tumors

Epstein-Barr virus (EBV) infection is highly prevalent in the population and is known to be associated with a variety of human tumors, such as nasopharyngeal carcinoma, gastric cancer, and lymphoma; however, the mechanisms of EBV carcinogenesis remain unclear. Recent studies have revealed that many non-coding RNAs participate in the regulation of proliferation, migration, invasion, and other processes in EBV-associated tumor, and the interaction between ncRNAs and the potential target genes has gradually become a research hotspot. Therefore, here, we discuss the expression and roles of ncRNAs in EBV-associated epithelial tumors.

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.

Essential role of inverted repeat in Epstein-Barr virus IR-1 in B cell transformation; geographical variation of the viral genome

Many regions of the Epstein–Barr virus (EBV) genome, repeated and unique sequences, contribute to the geographical variation observed between strains. Here we use a large alignment of curated EBV genome sequences to identify major sites of variation in the genome of type 1 EBV strains; the CAO deletion in latent membrane protein 1 (LMP1) is the most frequent major indel present in the unique regions of EBV strains from various parts of the world. Principal component analysis was used to identify patterns of sequence variation and nucleotide positions in the sequences that can distinguish EBV from some different geographical regions. Viral genome sequence variation also affects interpretation of genetic content; known genes, origins of replication and gene expression control regions explain most of the viral genome but there are still a few sections of unknown function. One of these EBV genome regions contains a large inverted repeat sequence (invR) within the IR-1 major internal repeat array. We deleted this invR sequence and showed that this abolished the ability of the virus to transform human B cells into lymphoblastoid cell lines.

This article is part of the theme issue ‘Silent cancer agents: multi-disciplinary modelling of human DNA oncoviruses’.

Idiosyncrasies of Viral Noncoding RNAs Provide Insights into Host Cell Biology

Like their host cells, many viruses express noncoding RNAs (ncRNAs). Despite the technical challenge of ascribing function to ncRNAs, diverse biological roles for virally expressed ncRNAs have been described, including regulation of viral replication, modulation of host gene expression, host immune evasion, cellular survival, and cellular transformation. Insights into conserved interactions between viral ncRNAs and host cell machinery frequently lead to novel findings concerning host cell biology. In this review, we discuss the functions and biogenesis of ncRNAs produced by animal viruses. Specifically, we describe noncanonical pathways of microRNA (miRNA) biogenesis and novel mechanisms used by viruses to manipulate miRNA and messenger RNA stability. We also highlight recent advances in understanding the function of viral long ncRNAs and circular RNAs.

The RNA encoding the microtubule-associated protein tau has extensive structure that affects its biology

Tauopathies are neurodegenerative diseases that affect millions of people worldwide including those with Alzheimer’s disease. While many efforts have focused on understanding the role of tau protein in neurodegeneration, there has been little done to systematically analyze and study the structures within tau’s encoding RNA and their connection to disease pathology. Knowledge of RNA structure can provide insights into disease mechanisms and how to affect protein production for therapeutic benefit. Using computational methods based on thermodynamic stability and evolutionary conservation, we identified structures throughout the tau pre-mRNA, especially at exon-intron junctions and within the 5′ and 3′ untranslated regions (UTRs). In particular, structures were identified at twenty exon-intron junctions. The 5′ UTR contains one structured region, which lies within a known internal ribosome entry site. The 3′ UTR contains eight structured regions, including one that contains a polyadenylation signal. A series of functional experiments were carried out to assess the effects of mutations associated with mis-regulation of alternative splicing of exon 10 and to identify regions of the 3′ UTR that contain cis-regulatory elements. These studies defined novel structural regions within the mRNA that affect stability and pre-mRNA splicing and may lead to new therapeutic targets for treating tau-associated diseases.

Epstein-Barr virus infection as a predisposing factor for multiple sclerosis. An update from molecular biology, immunology and epidemiology

Introduction: Epstein-Barr virus is an infectious agent used to immortalize and induce polyclonal activation of B cells. It has been widely described that this virus produces changes in the cells it infects and in the immune response, and stimulates the development of autoimmune diseases.Objective: To characterize the association between Epstein-Barr virus and multiple sclerosis described in current scientific literature.Materials and methods: A 59-years range literature search was conducted in the PubMed, ScienceDirect, Redalyc and SciELO databases using the following MeSH terms: “Epstein-Barr virus, multiple sclerosis autoimmune diseases, autoimmune diseases of the nervous system”.Results: Many studies describe the association between Epstein-Barr virus and multiple sclerosis. It is believed that acute infection and viral reactivation promote the development of multiple sclerosis.Conclusions: It is necessary to conduct further research on the pathogenesis and morphophysiological and neuroimmunological changes –at the ecological, molecular, cellular, tissue, organic and systemic level– induced by the immune response and that favor the development of multiple sclerosis.

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.

A Genome-Wide Epstein-Barr Virus Polyadenylation Map and Its Antisense RNA to EBNA

Epstein-Barr virus (EBV) is a ubiquitous human pathogen associated with Burkitt's lymphoma and nasopharyngeal carcinoma. Although the EBV genome harbors more than a hundred genes, a full transcription map with EBV polyadenylation profiles remains unknown. To elucidate the 3' ends of all EBV transcripts genome-wide, we performed the first comprehensive analysis of viral polyadenylation sites (pA sites) using our previously reported polyadenylation sequencing (PA-seq) technology. We identified that EBV utilizes a total of 62 pA sites in JSC-1, 60 in Raji, and 53 in Akata cells for the expression of EBV genes from both plus and minus DNA strands; 42 of these pA sites are commonly used in all three cell lines. The majority of identified pA sites were mapped to the intergenic regions downstream of previously annotated EBV open reading frames (ORFs) and viral promoters. pA sites lacking an association with any known EBV genes were also identified, mostly for the minus DNA strand within the EBNA locus, a major locus responsible for maintenance of viral latency and cell transformation. The expression of these novel antisense transcripts to EBNA were verified by 3' rapid amplification of cDNA ends (RACE) and Northern blot analyses in several EBV-positive (EBV+) cell lines. In contrast to EBNA RNA expressed during latency, expression of EBNA-antisense transcripts, which is restricted in latent cells, can be significantly induced by viral lytic infection, suggesting potential regulation of viral gene expression by EBNA-antisense transcription during lytic EBV infection. Our data provide the first evidence that EBV has an unrecognized mechanism that regulates EBV reactivation from latency.IMPORTANCE Epstein-Barr virus represents an important human pathogen with an etiological role in the development of several cancers. By elucidation of a genome-wide polyadenylation landscape of EBV in JSC-1, Raji, and Akata cells, we have redefined the EBV transcriptome and mapped individual polymerase II (Pol II) transcripts of viral genes to each one of the mapped pA sites at single-nucleotide resolution as well as the depth of expression. By unveiling a new class of viral lytic RNA transcripts antisense to latent EBNAs, we provide a novel mechanism of how EBV might control the expression of viral latent genes and lytic infection. Thus, this report takes another step closer to understanding EBV gene structure and expression and paves a new path for antiviral approaches.

Stable Intronic Sequence RNAs (sisRNAs): An Expanding Universe

Intronic sequences are often regarded as 'nonsense' transcripts that are rapidly degraded. We highlight here recent studies on intronic sequences that play regulatory roles as long noncoding RNAs (lncRNAs) which are classified as sisRNAs. Interestingly, sisRNAs come in different forms and are produced via a variety of ways. They regulate genes at the DNA, RNA, and protein levels, and frequently engage in autoregulatory feedback loops to ensure cellular homeostasis under normal and stress conditions. Future directions, evolutionary insights, and potential implications of dysregulated sisRNAs are also discussed, especially in relation to human pathogenesis.

The Structure-To-Function Relationships of Gammaherpesvirus-Encoded Long Non-Coding RNAs and Their Contributions to Viral Pathogenesis

Advances in next-generation sequencing have facilitated the discovery of a multitude of long non-coding RNAs (lncRNAs) with pleiotropic functions in cellular processes, disease, and viral pathogenesis. It came as no surprise when viruses were also revealed to transcribe their own lncRNAs. Among them, gammaherpesviruses, one of the three subfamilies of the Herpesviridae, code their largest number. These structurally and functionally intricate non-coding (nc) transcripts modulate cellular and viral gene expression to maintain viral latency or prompt lytic reactivation. These lncRNAs allow for the virus to escape cytosolic surveillance, sequester, and re-localize essential cellular factors and modulate the cell cycle and proliferation. Some viral lncRNAs act as “messenger molecules”, transferring information about viral infection to neighboring cells. This broad range of lncRNA functions is achieved through lncRNA structure-mediated interactions with effector molecules of viral and host origin, including other RNAs, proteins and DNAs. In this review, we discuss examples of gammaherpesvirus-encoded lncRNAs, emphasize their unique structural attributes, and link them to viral life cycle, pathogenesis, and disease progression. We will address their potential as novel targets for drug discovery and propose future directions to explore lncRNA structure and function relationship.

Determination and genome-wide analysis of Epstein-Barr virus (EBV) sequences in EBV-associated gastric carcinoma from Guangdong, an endemic area of nasopharyngeal carcinoma

About 10 % of gastric carcinoma worldwide is associated with EBV, which is defined as EBV-associated gastric carcinoma (EBVaGC). To date, EBV sequence data from EBVaGC in Guangdong, China, an endemic area of nasopharyngeal carcinoma (NPC), are not available. In the present study, two EBV genomes from EBVaGC specimens from Guangdong (designated as GDGC1 and GDGC2) were determined by next-generation sequencing, de novo assembly and joining of contigs by Sanger sequencing. In addition, we sequenced EBV from two Korean EBVaGC cell lines, YCCEL1 and SNU-719. Genomic diversity, including single nucleotide polymorphisms (SNPs) and insertions and deletions (indels), phylogenetic analysis and rates of protein evolution, was performed using bioinformatics software. The four gastric carcinoma-derived EBV (GC-EBV) were all type I. Compared with the reference EBV genome, a total of 1815 SNPs (146 indels), 1519 SNPs (106 indels), 1812 SNPs (126 indels) and 1484 SNPs (106 indels) were found in GDGC1, GDGC2, YCCEL1 and SNU-719, respectively. These variations were distributed across the entire genome, especially in latent genes. In contrast, the sequences of promoters and non-coding RNAs were strictly conserved. Phylogenetic analyses suggested the presence of at least two parental lineages of EBV among the GC-EBV genomes. Rates of protein evolution analyses showed that lytic genes were under purifying selection; in contrast, latency genes were under positive selection. In conclusion, this study determined the EBV genomes in EBVaGC from Guangdong and performed a detailed genome-wide analysis of GC-EBV, which would be helpful for further understanding of the relationship between EBV genomic variation and EBVaGC carcinogenesis.

Adaptation by copy number variation in monopartite viruses

•

Viral gene amplification allows rapid adaptation, especially for large DNA viruses.

•

Amplifications often incur fitness costs for viral replication.

•

Adaptive mutations can arise in the expanded locus or elsewhere enabling collapse.

•

Genome amplifications provide genetic substrate for the evolution of modified or new functions.

Viruses evolve rapidly in response to host defenses and to exploit new niches. Gene amplification, a common adaptive mechanism in prokaryotes, archaea, and eukaryotes, has also contributed to viral evolution, especially of large DNA viruses. In experimental systems, gene amplification is one mechanism for rapidly overcoming selective pressures. Because the amplification generally incurs a fitness cost, emergence of adaptive point mutations within the amplified locus or elsewhere in the genome can enable collapse of the locus back to a single copy. Evidence of gene amplification followed by subfunctionalization or neofunctionalization of the copies is apparent by the presence of families of paralogous genes in many DNA viruses. These observations suggest that copy number variation has contributed broadly to virus evolution.

Generation of Drosophila sisRNAs by Independent Transcription from Cognate Introns

Although stable intronic sequence RNAs (sisRNAs) are conserved in plants and animals, their functional significance is still unclear. We identify a pool of polyadenylated maternally deposited sisRNAs in Drosophila melanogaster. These sisRNAs can be generated by independent transcription from the cognate introns. The ovary-specific poly(A) polymerase Wispy mediates the polyadenylation of maternal sisRNAs and confers their stability as maternal transcripts. A developmentally regulated sisRNA sisR-3 represses the expression of a long noncoding RNA CR44148 and is required during development. Our results expand the pool of sisRNAs and suggest that sisRNAs perform regulatory functions during development in Drosophila.

•

Identification of polyadenylated sisRNAs

•

sisRNAs can be produced from independent transcription

•

sisR-3 regulates a long noncoding RNA

•

sisR-3 is required during development

Rapid CRISPR/Cas9-Mediated Cloning of Full-Length Epstein-Barr Virus Genomes from Latently Infected Cells

Herpesviruses have relatively large DNA genomes of more than 150 kb that are difficult to clone and sequence. Bacterial artificial chromosome (BAC) cloning of herpesvirus genomes is a powerful technique that greatly facilitates whole viral genome sequencing as well as functional characterization of reconstituted viruses. We describe recently invented technologies for rapid BAC cloning of herpesvirus genomes using CRISPR/Cas9-mediated homology-directed repair. We focus on recent BAC cloning techniques of Epstein-Barr virus (EBV) genomes and discuss the possible advantages of a CRISPR/Cas9-mediated strategy comparatively with precedent EBV-BAC cloning strategies. We also describe the design decisions of this technology as well as possible pitfalls and points to be improved in the future. The obtained EBV-BAC clones are subjected to long-read sequencing analysis to determine complete EBV genome sequence including repetitive regions. Rapid cloning and sequence determination of various EBV strains will greatly contribute to the understanding of their global geographical distribution. This technology can also be used to clone disease-associated EBV strains and test the hypothesis that they have special features that distinguish them from strains that infect asymptomatically.

EBV-encoded miRNAs target ATM-mediated response in nasopharyngeal carcinoma

Nasopharyngeal carcinoma (NPC) is a highly invasive epithelial malignancy that is prevalent in southern China and Southeast Asia. It is consistently associated with latent Epstein-Barr virus (EBV) infection. In NPC, miR-BARTs, the EBV-encoded miRNAs derived from BamH1-A rightward transcripts, are abundantly expressed and contribute to cancer development by targeting various cellular and viral genes. In this study, we establish a comprehensive transcriptional profile of EBV-encoded miRNAs in a panel of NPC patient-derived xenografts and an EBV-positive NPC cell line by small RNA sequencing. Among the 40 miR-BARTs, predominant expression of 22 miRNAs was consistently detected in these tumors. Among the abundantly expressed EBV-miRNAs, BART5-5p, BART7-3p, BART9-3p, and BART14-3p could negatively regulate the expression of a key DNA double-strand break (DSB) repair gene, ataxia telangiectasia mutated (ATM), by binding to multiple sites on its 3'-UTR. Notably, the expression of these four miR-BARTs represented more than 10% of all EBV-encoded miRNAs in tumor cells, while downregulation of ATM expression was commonly detected in all of our tested sequenced samples. In addition, downregulation of ATM was also observed in primary NPC tissues in both qRT-PCR (16 NP and 45 NPC cases) and immunohistochemical staining (35 NP and 46 NPC cases) analysis. Modulation of ATM expression by BART5-5p, BART7-3p, BART9-3p, and BART14-3p was demonstrated in the transient transfection assays. These findings suggest that EBV uses miRNA machinery as a key mechanism to control the ATM signaling pathway in NPC cells. By suppressing these endogenous miR-BARTs in EBV-positive NPC cells, we further demonstrated the novel function of miR-BARTs in inhibiting Zta-induced lytic reactivation. These findings imply that the four viral miRNAs work co-operatively to modulate ATM activity in response to DNA damage and to maintain viral latency, contributing to the tumorigenesis of NPC. © 2017 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

RNA2DMut: a web tool for the design and analysis of RNA structure mutations

With the widespread application of high-throughput sequencing, novel RNA sequences are being discovered at an astonishing rate. The analysis of function, however, lags behind. In both the cis- and trans-regulatory functions of RNA, secondary structure (2D base-pairing) plays essential regulatory roles. In order to test RNA function, it is essential to be able to design and analyze mutations that can affect structure. This was the motivation for the creation of the RNA2DMut web tool. With RNA2DMut, users can enter in RNA sequences to analyze, constrain mutations to specific residues, or limit changes to purines/pyrimidines. The sequence is analyzed at each base to determine the effect of every possible point mutation on 2D structure. The metrics used in RNA2DMut rely on the calculation of the Boltzmann structure ensemble and do not require a robust 2D model of RNA structure for designing mutations. This tool can facilitate a wide array of uses involving RNA: for example, in designing and evaluating mutants for biological assays, interrogating RNA-protein interactions, identifying key regions to alter in SELEX experiments, and improving RNA folding and crystallization properties for structural biology. Additional tools are available to help users introduce other mutations (e.g., indels and substitutions) and evaluate their effects on RNA structure. Example calculations are shown for five RNAs that require 2D structure for their function: the MALAT1 mascRNA, an influenza virus splicing regulatory motif, the EBER2 viral noncoding RNA, the Xist lncRNA repA region, and human Y RNA 5. RNA2DMut can be accessed at https://rna2dmut.bb.iastate.edu/.

Genetics and Molecular Biology of Epstein-Barr Virus-Encoded BART MicroRNA: A Paradigm for Viral Modulation of Host Immune Response Genes and Genome Stability

Epstein-Barr virus, a ubiquitous human herpesvirus, is associated through epidemiologic evidence with common autoimmune syndromes and cancers. However, specific genetic mechanisms of pathogenesis have been difficult to identify. In this review, the author summarizes evidence that recently discovered noncoding RNAs termed microRNA encoded by Epstein-Barr virus BARF (BamHI A right frame) termed BART (BamHI A right transcripts) are modulators of human immune response genes and genome stability in infected and bystander cells. BART expression is apparently regulated by complex feedback loops with the host immune response regulatory NF-κB transcription factors. EBV-encoded BZLF-1 (ZEBRA) protein could also regulate BART since ZEBRA contains a terminal region similar to ankyrin proteins such as IκBα that regulate host NF-κB. BALF-2 (BamHI A left frame transcript), a viral homologue of the immunoglobulin and T cell receptor gene recombinase RAG-1 (recombination-activating gene-1), may also be coregulated with BART since BALF-2 regulatory sequences are located near the BART locus. Viral-encoded microRNA and viral mRNA transferred to bystander cells through vesicles, defective viral particles, or other mechanisms suggest a new paradigm in which bystander or hit-and-run mechanisms enable the virus to transiently or chronically alter human immune response genes as well as the stability of the human genome.

Comparative analysis of 22 Epstein-Barr virus genomes from diseased and healthy individuals

Epstein-Barr virus (EBV) infects most of the world's population and is causally associated with several human cancers, but little is known about how EBV genetic variations might influence EBV-associated diseases and their geographical patterns. In the present study, 22 EBV whole-genome sequences from diseased and healthy individuals were analysed to explore EBV sequence variations at the whole-genome level. We found that the 22 EBV genomes were generally highly similar to each other at the genome level. However, varying degrees of genetic diversity were detected across the entire genome, especially in the latent genes. In contrast, the sequences of promoters and non-coding RNAs were strictly conserved. These findings suggested that both latent genes and non-coding RNAs play important roles in the EBV life cycle. When we investigated changes in known T-cell epitopes in some latent and lytic proteins, we observed that some T-cell epitopes were changed, while others were conserved. These findings indicate that the effect of EBV variations in protein sequences that seem to have been selected by the host immune system should be considered when conducting EBV-targeted immunotherapy. Taken together, our results provide a global view of EBV genome sequence variation, which not only is important for designing vaccines and immunotherapy for EBV but also adds to the understanding of EBV biology and the relationships between viral sequence variation and EBV-associated diseases.

A proximity-dependent assay for specific RNA-protein interactions in intact cells

The proximity ligation assay (PLA) is an immune staining method that detects protein-protein interactions in fixed cells. We describe here RNA-PLA, a simple adaptation of this technology that allows the detection of specific RNA-protein interactions in fixed cells by using a DNA oligonucleotide that hybridizes to a target RNA in combination with an antibody that recognizes the protein bound to the target RNA. Stable and transient RNA-protein interactions can be readily detected by generation of a fluorescent signal in discrete compartments in intact fixed cells with high specificity. We demonstrate that this approach requires the colocalization of the binding protein and its RNA target in the same cellular compartment, use of an oligonucleotide complementary to the target RNA, and the presence of a binding site for the protein in the target RNA.

Identification of IgH gene rearrangement and immunophenotype in an animal model of Epstein-Barr virus-associated lymphomas

Epstein-Barr virus (EBV) is a human oncogenic herpesvirus associated with lymphoma and nasopharyngeal carcinoma. Because the susceptible hosts of EB virus are limited to human and cotton-top tamarins (Saguinus oedipus), there have been no appropriate animal models until the lymphoma model induced by EBV in human peripheral blood lymphocyte (hu-PBL)/SCID chimeric mice was reported. However, it is still controversial whether the EBV-associated lymphoma induced in hu-PBL/SCID mice is a monoclonal tumor. In this study, we transplanted normal human peripheral blood lymphocytes (hu-PBL) from six donors infected with EBV into SCID mice to construct hu-PBL/SCID chimeric mice. The induced tumors were found in the mediastinum or abdominal cavity of SCID mice. Microscopic observation exhibited tumor cells that were large and had a plasmablastic, centroblastic or immunoblastic-like appearance. Immunophenotyping assays showed the induced tumors were LCA-positive, CD20/CD79a-positive (markers of B cells), and CD3/CD45RO-negative (markers of T cells). A human-specific Alu sequence could be amplified by Alu-PCR. This confirmed that induced tumors were B-cell lymphomas originating from the transplanted human lymphocytes rather than mouse cells. EBER in situ hybridization detected positive signals in the nuclei of the tumor cells. Expression of EBV-encoded LMP1, EBNA-1, and EBNA-2 in the tumors was significantly positive. PCR-based capillary electrophoresis analysis of IgH gene rearrangement revealed a monoclonal peak and single amplification product in all six cases of induced tumors. This indicated that EBV can induce monoclonal proliferation of human B lymphocytes and promotes the development of lymphoma. J. Med. Virol. 88:1804-1813, 2016. © 2016 Wiley Periodicals, Inc.

Consideration of Epstein-Barr Virus-Encoded Noncoding RNAs EBER1 and EBER2 as a Functional Backup of Viral Oncoprotein Latent Membrane Protein 1

The Epstein-Barr virus (EBV)-encoded noncoding RNAs EBER1 and EBER2 are highly abundant through all four latency stages of EBV infection (III-II-I-0) and have been associated with an oncogenic phenotype when expressed in cell lines cultured in vitro. In vivo, EBV-infected B cells derived from freshly isolated lymphocytes show that EBER1/2 deletion does not impair viral latency. Based on published quantitative proteomics data from BJAB cells expressing EBER1 and EBER2, we propose that the EBERs, through their activation of AKT in a B-cell-specific manner, are a functionally redundant backup of latent membrane protein 1 (LMP1)-an essential oncoprotein in EBV-associated malignancies, with a main role in AKT activation. Our proposed model may explain the lack of effect on viral latency establishment in EBER-minus EBV infection.

Viral Noncoding RNAs in Cancer Biology

Over 12 % of all human cancers are caused by oncoviruses, primarily including Epstein-Barr virus (EBV), high-risk human papillomaviruses (HPVs), hepatitis B and C viruses (HBV and HCV, respectively), and Kaposi's sarcoma herpesvirus (KSHV). In addition to viral oncoproteins, a variety of noncoding RNAs (ncRNAs) produced by oncoviruses have been recognized as important cofactors that contribute to the oncogenic events. In this chapter, we will focus on the recent understanding of the long and short noncoding RNAs, as well as microRNAs of the viruses, and discuss their roles in the biology of multistep oncogenesis mediated by established human oncoviruses.

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.

New Noncoding Lytic Transcripts Derived from the Epstein-Barr Virus Latency Origin of Replication, oriP, Are Hyperedited, Bind the Paraspeckle Protein, NONO/p54nrb, and Support Viral Lytic Transcription

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

IMPORTANCE

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

Viral noncoding RNAs: more surprises

Eukaryotic cells produce several classes of long and small noncoding RNA (ncRNA). Many DNA and RNA viruses synthesize their own ncRNAs. Like their host counterparts, viral ncRNAs associate with proteins that are essential for their stability, function, or both. Diverse biological roles--including the regulation of viral replication, viral persistence, host immune evasion, and cellular transformation--have been ascribed to viral ncRNAs. In this review, we focus on the multitude of functions played by ncRNAs produced by animal viruses. We also discuss their biogenesis and mechanisms of action.

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

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.