Identification of new viral genes and transcript isoforms during Epstein-Barr virus reactivation using RNA-Seq

Managing Viral Emerging Infectious Diseases via Current and Future Molecular Diagnostics

Emerging viral infectious diseases have been a constant threat to global public health in recent times. In managing these diseases, molecular diagnostics has played a critical role. Molecular diagnostics involves the use of various technologies to detect the genetic material of various pathogens, including viruses, in clinical samples. One of the most commonly used molecular diagnostics technologies for detecting viruses is polymerase chain reaction (PCR). PCR amplifies specific regions of the viral genetic material in a sample, making it easier to detect and identify viruses. PCR is particularly useful for detecting viruses that are present in low concentrations in clinical samples, such as blood or saliva. Another technology that is becoming increasingly popular for viral diagnostics is next-generation sequencing (NGS). NGS can sequence the entire genome of a virus present in a clinical sample, providing a wealth of information about the virus, including its genetic makeup, virulence factors, and potential to cause an outbreak. NGS can also help identify mutations and discover new pathogens that could affect the efficacy of antiviral drugs and vaccines. In addition to PCR and NGS, there are other molecular diagnostics technologies that are being developed to manage emerging viral infectious diseases. One of these is CRISPR-Cas, a genome editing technology that can be used to detect and cut specific regions of viral genetic material. CRISPR-Cas can be used to develop highly specific and sensitive viral diagnostic tests, as well as to develop new antiviral therapies. In conclusion, molecular diagnostics tools are critical for managing emerging viral infectious diseases. PCR and NGS are currently the most commonly used technologies for viral diagnostics, but new technologies such as CRISPR-Cas are emerging. These technologies can help identify viral outbreaks early, track the spread of viruses, and develop effective antiviral therapies and vaccines.

Integrative profiling of Epstein-Barr virus transcriptome using a multiplatform approach

BACKGROUND

Epstein-Barr virus (EBV) is an important human pathogenic gammaherpesvirus with carcinogenic potential. The EBV transcriptome has previously been analyzed using both Illumina-based short read-sequencing and Pacific Biosciences RS II-based long-read sequencing technologies. Since the various sequencing methods have distinct strengths and limitations, the use of multiplatform approaches have proven to be valuable. The aim of this study is to provide a more complete picture on the transcriptomic architecture of EBV.

METHODS

In this work, we apply the Oxford Nanopore Technologies MinION (long-read sequencing) platform for the generation of novel transcriptomic data, and integrate these with other's data generated by another LRS approach, Pacific BioSciences RSII sequencing and Illumina CAGE-Seq and Poly(A)-Seq approaches. Both amplified and non-amplified cDNA sequencings were applied for the generation of sequencing reads, including both oligo-d(T) and random oligonucleotide-primed reverse transcription. EBV transcripts are identified and annotated using the LoRTIA software suite developed in our laboratory.

RESULTS

This study detected novel genes embedded into longer host genes containing 5'-truncated in-frame open reading frames, which potentially encode N-terminally truncated proteins. We also detected a number of novel non-coding RNAs and transcript length isoforms encoded by the same genes but differing in their start and/or end sites. This study also reports the discovery of novel splice isoforms, many of which may represent altered coding potential, and of novel replication-origin-associated transcripts. Additionally, novel mono- and multigenic transcripts were identified. An intricate meshwork of transcriptional overlaps was revealed.

CONCLUSIONS

An integrative approach applying multi-technique sequencing technologies is suitable for reliable identification of complex transcriptomes because each techniques has different advantages and limitations, and the they can be used for the validation of the results obtained by a particular approach.

Pathogenic Role of Epstein-Barr Virus in Lung Cancers

Human oncogenic viruses account for at least 12% of total cancer cases worldwide. Epstein–Barr virus (EBV) is the first identified human oncogenic virus and it alone causes ~200,000 cancer cases and ~1.8% of total cancer-related death annually. Over the past 40 years, increasing lines of evidence have supported a causal link between EBV infection and a subgroup of lung cancers (LCs). In this article, we review the current understanding of the EBV-LC association and the etiological role of EBV in lung carcinogenesis. We also discuss the clinical impact of the knowledge gained from previous research, challenges, and future directions in this field. Given the high clinical relevance of EBV-LC association, there is an urgent need for further investigation on this topic.

Role of Virally Encoded Circular RNAs in the Pathogenicity of Human Oncogenic Viruses

Human oncogenic viruses are a group of important pathogens that etiologically contribute to at least 12% of total cancer cases in the world. As an emerging class of non-linear regulatory RNA molecules, circular RNAs (circRNAs) have gained increasing attention as a crucial player in the regulation of signaling pathways involved in viral infection and oncogenesis. With the assistance of current circRNA enrichment and detection technologies, numerous novel virally-encoded circRNAs (vcircRNAs) have been identified in the human oncogenic viruses, initiating an exciting new era of vcircRNA research. In this review, we discuss the current understanding of the roles of vcircRNAs in the respective viral infection cycles and in virus-associated pathogenesis.

Comprehensive Epstein-Barr Virus Transcriptome by RNA-Sequencing in Angioimmunoblastic T Cell Lymphoma (AITL) and Other Lymphomas

The Epstein-Barr virus (EBV) is associated with angioimmunoblastic T cell lymphoma (AITL) in more than 80% of cases. Few studies have focused on this association and it is not clear now what role the virus plays in this pathology. We used next-generation sequencing (NGS) to study EBV transcriptome in 14 AITLs compared to 21 other lymphoma samples and 11 cell lines including 4 lymphoblastoid cell lines (LCLs). Viral transcripts were recovered using capture probes and sequencing was performed on Illumina. Bam-HI A rightward transcripts (BARTs) were the most latency transcripts expressed in AITLs, suggesting they may play a role in this pathology. Thus, BARTs, already described as highly expressed in carcinoma cells, are also very present in AITLs and other lymphomas. They were poorly expressed in cell lines other than LCLs. AITLs showed a latency IIc, with BNLF2a gene expression. For most AITLs, BCRF1, which encodes a homologous protein of human interleukin 10, vIL-10, was in addition expressed. This co-expression can contribute to immune escape and survival of infected cells. Considering these results, it can be assumed that EBV plays a pathogenic role in AITLs.

Widespread Traces of Lytic Kaposi Sarcoma-Associated Herpesvirus in Primary Effusion Lymphoma at Single-Cell Resolution

Cancer cells of primary effusion lymphoma (PEL) often contain both Kaposi sarcoma-associated herpesvirus (KSHV) and Epstein-Barr virus (EBV). We measured the interplay of human, KSHV, and EBV transcription in a cell culture model of PEL using single-cell RNA sequencing. The data detect widespread trace expression of lytic KSHV genes.

Cancer cells of primary effusion lymphoma (PEL) often contain both Kaposi sarcoma-associated herpesvirus (KSHV) and Epstein-Barr virus (EBV). We measured the interplay of human, KSHV, and EBV transcription in a cell culture model of PEL using single-cell RNA sequencing. The data detect widespread trace expression of lytic KSHV genes.

Histone Loaders CAF1 and HIRA Restrict Epstein-Barr Virus B-Cell Lytic Reactivation

Epstein-Barr virus (EBV) infects 95% of adults worldwide and causes infectious mononucleosis. EBV is associated with endemic Burkitt lymphoma, Hodgkin lymphoma, posttransplant lymphomas, nasopharyngeal and gastric carcinomas. In these cancers and in most infected B-cells, EBV maintains a state of latency, where nearly 80 lytic cycle antigens are epigenetically suppressed. To gain insights into host epigenetic factors necessary for EBV latency, we recently performed a human genome-wide CRISPR screen that identified the chromatin assembly factor CAF1 as a putative Burkitt latency maintenance factor. CAF1 loads histones H3 and H4 onto newly synthesized host DNA, though its roles in EBV genome chromatin assembly are uncharacterized. Here, we found that CAF1 depletion triggered lytic reactivation and virion secretion from Burkitt cells, despite also strongly inducing interferon-stimulated genes. CAF1 perturbation diminished occupancy of histones 3.1 and 3.3 and of repressive histone 3 lysine 9 and 27 trimethyl (H3K9me3 and H3K27me3) marks at multiple viral genome lytic cycle regulatory elements. Suggestive of an early role in establishment of latency, EBV strongly upregulated CAF1 expression in newly infected primary human B-cells prior to the first mitosis, and histone 3.1 and 3.3 were loaded on the EBV genome by this time point. Knockout of CAF1 subunit CHAF1B impaired establishment of latency in newly EBV-infected Burkitt cells. A nonredundant latency maintenance role was also identified for the DNA synthesis-independent histone 3.3 loader histone regulatory homologue A (HIRA). Since EBV latency also requires histone chaperones alpha thalassemia/mental retardation syndrome X-linked chromatin remodeler (ATRX) and death domain-associated protein (DAXX), EBV coopts multiple host histone pathways to maintain latency, and these are potential targets for lytic induction therapeutic approaches.IMPORTANCE Epstein-Barr virus (EBV) was discovered as the first human tumor virus in endemic Burkitt lymphoma, the most common childhood cancer in sub-Saharan Africa. In Burkitt lymphoma and in 200,000 EBV-associated cancers per year, epigenetic mechanisms maintain viral latency, during which lytic cycle factors are silenced. This property complicated EBV's discovery and facilitates tumor immunoevasion. DNA methylation and chromatin-based mechanisms contribute to lytic gene silencing. Here, we identified histone chaperones CAF1 and HIRA, which have key roles in host DNA replication-dependent and replication-independent pathways, respectively, as important for EBV latency. EBV strongly upregulates CAF1 in newly infected B-cells, where viral genomes acquire histone 3.1 and 3.3 variants prior to the first mitosis. Since histone chaperones ATRX and DAXX also function in maintenance of EBV latency, our results suggest that EBV coopts multiple histone pathways to reprogram viral genomes and highlight targets for lytic induction therapeutic strategies.

The BHLF1 Locus of Epstein-Barr Virus Contributes to Viral Latency and B-Cell Immortalization

The Epstein-Barr virus (EBV) BHLF1 gene encodes an abundant linear and several circular RNAs believed to perform noncoding functions during virus replication, although an open reading frame (ORF) is retained among an unknown percentage of EBV isolates. Evidence suggests that BHLF1 is also transcribed during latent infection, which prompted us to investigate the contribution of this locus to latency. Analysis of transcripts transiting BHLF1 revealed that its transcription is widespread among B-cell lines supporting the latency I or III program of EBV protein expression and is more complex than originally presumed. EBV-negative Burkitt lymphoma cell lines infected with either wild-type or two different BHLF1 mutant EBVs were initially indistinguishable in supporting latency III. However, cells infected with BHLF1 ^- virus ultimately transitioned to the more restrictive latency I program, whereas cells infected with wild-type virus either sustained latency III or transitioned more slowly to latency I. Upon infection of primary B cells, which require latency III for growth in vitro, both BHLF1 ^- viruses exhibited variably reduced immortalization potential relative to the wild-type virus. Finally, in transfection experiments, efficient protein expression from an intact BHLF1 ORF required the EBV posttranscriptional regulator protein SM, whose expression is limited to the replicative cycle. Thus, one way in which BHLF1 may contribute to latency is through a mechanism, possibly mediated or regulated by a long noncoding RNA, that supports latency III critical for the establishment of EBV latency and lifelong persistence within its host, whereas any retained protein-dependent function of BHLF1 may be restricted to the replication cycle.IMPORTANCE Epstein-Barr virus (EBV) has significant oncogenic potential that is linked to its latent infection of B lymphocytes, during which virus replication is not supported. The establishment of latent infection, which is lifelong and can precede tumor development by years, requires the concerted actions of nearly a dozen EBV proteins and numerous small non-protein-coding RNAs. Elucidating how these EBV products contribute to latency is crucial for understanding EBV's role in specific malignancies and, ultimately, for clinical intervention. Historically, EBV genes that contribute to virus replication have been excluded from consideration of a role in latency, primarily because of the general incompatibility between virus production and cell survival. However, here, we provide evidence that the genetic locus containing one such gene, BHLF1, indeed contributes to key aspects of EBV latency, including its ability to promote the continuous growth of B lymphocytes, thus providing significant new insight into EBV biology and oncogenic potential.

Using Direct RNA Nanopore Sequencing to Deconvolute Viral Transcriptomes

The genomes of DNA viruses encode deceptively complex transcriptomes evolved to maximize coding potential within the confines of a relatively small genome. Defining the full range of viral RNAs produced during an infection is key to understanding the viral replication cycle and its interactions with the host cell. Traditional short-read (Illumina) sequencing approaches are problematic in this setting due to the difficulty of assigning short reads to individual RNAs in regions of transcript overlap and to the biases introduced by the required recoding and amplification steps. Additionally, different methodologies may be required to analyze the 5' and 3' ends of RNAs, which increases both cost and effort. The advent of long-read nanopore sequencing simplifies this approach by providing a single assay that captures and sequences full length RNAs, either in cDNA or native RNA form. The latter is particularly appealing as it reduces known recoding biases whilst allowing more advanced analyses such as estimation of poly(A) tail length and the detection of RNA modifications including N6 -methyladenosine. Using herpes simplex virus (HSV)-infected primary fibroblasts as a template, we provide a step-by-step guide to the production of direct RNA sequencing libraries suitable for sequencing using Oxford Nanopore Technologies platforms and provide a simple computational approach to deriving a high-quality annotation of the HSV transcriptome from the resulting sequencing data. © 2020 by John Wiley & Sons, Inc. Basic Protocol 1: Productive infection of primary fibroblasts with herpes simplex virus Support Protocol: Cell passage and plating of primary fibroblasts Basic Protocol 2: Preparation and sequencing of dRNA-seq libraries from virus-infected cells Basic Protocol 3: Processing, alignment, and analysis of dRNA-seq datasets.

High-Throughput Sequence Analysis of Peripheral T-Cell Lymphomas Indicates Subtype-Specific Viral Gene Expression Patterns and Immune Cell Microenvironments

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

Roles of Non-coding RNAs During Herpesvirus Infection

Non-coding RNAs (ncRNAs) play essential roles in multiple aspects of the life cycles of herpesviruses and contribute to lifelong persistence of herpesviruses within their respective hosts. In this chapter, we discuss the types of ncRNAs produced by the different herpesvirus families during infection, some of the cellular ncRNAs manipulated by these viruses, and the overall contributions of ncRNAs to the viral life cycle, influence on the host environment, and pathogenesis.

Detection of Epstein-Barr Virus Infection in Non-Small Cell Lung Cancer

Previous investigations proposed a link between the Epstein-Barr virus (EBV) and lung cancer (LC), but the results are highly controversial largely due to the insufficient sample size and the inherent limitation of the traditional viral screening methods such as PCR. Unlike PCR, current next-generation sequencing (NGS) utilizes an unbiased method for the global assessment of all exogenous agents within a cancer sample with high sensitivity and specificity. In our current study, we aim to resolve this long-standing controversy by utilizing our unbiased NGS-based informatics approaches in conjunction with traditional molecular methods to investigate the role of EBV in a total of 1127 LC. In situ hybridization analysis of 110 LC and 10 normal lung samples detected EBV transcripts in 3 LC samples. Comprehensive virome analyses of RNA sequencing (RNA-seq) data sets from 1017 LC and 110 paired adjacent normal lung specimens revealed EBV transcripts in three lung squamous cell carcinoma and one lung adenocarcinoma samples. In the sample with the highest EBV coverage, transcripts from the BamHI A region accounted for the majority of EBV reads. Expression of EBNA-1, LMP-1 and LMP-2 was observed. A number of viral circular RNA candidates were also detected. Thus, we for the first time revealed a type II latency-like viral transcriptome in the setting of LC in vivo. The high-level expression of viral BamHI A transcripts in LC suggests a functional role of these transcripts, likely as long non-coding RNA. Analyses of cellular gene expression and stained tissue sections indicated an increased immune cell infiltration in the sample expressing high levels of EBV transcripts compared to samples expressing low EBV transcripts. Increased level of immune checkpoint blockade factors was also detected in the sample with higher levels of EBV transcripts, indicating an induced immune tolerance. Lastly, inhibition of immune pathways and activation of oncogenic pathways were detected in the sample with high EBV transcripts compared to the EBV-low LC indicating the direct regulation of cancer pathways by EBV. Taken together, our data support the notion that EBV likely plays a pathological role in a subset of LC.

Interpretation of EBV infection in pan-cancer genome considering viral life cycle: LiEB (Life cycle of Epstein-Barr virus)

We report a novel transcriptomic analysis workflow called LiEB (Life cycle of Epstein-Barr virus) to characterize distributions of oncogenic virus, Epstein-Barr virus (EBV) infection in human tumors. We analyzed 851 The Cancer Genome Atlas whole-transcriptome sequencing (WTS) data to investigate EBV infection by life cycle information using three-step LiEB workflow: 1) characterize virus infection generally; 2) align transcriptome sequences against a hybrid human-EBV genome, and 3) quantify EBV gene expression. Our results agreed with EBV infection status of public cell line data. Analysis in stomach adenocarcinoma identified EBV-positive cases involving PIK3CA mutations and/or CDKN2A silencing with biologically more determination, compared to previous reports. In this study, we found that a small number of colorectal adenocarcinoma cases involved with EBV lytic gene expression. Expression of EBV lytic genes was also observed in 3% of external colon cancer cohort upon WTS analysis. Gene set enrichment analysis showed elevated expression of genes related to E2F targeting and interferon-gamma responses in EBV-associated tumors. Finally, we suggest that interpretation of EBV life cycle is essential when analyzing its infection in tumors, and LiEB provides high capability of detecting EBV-positive tumors. Observation of EBV lytic gene expression in a subset of colon cancers warrants further research.

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.

Electronic supplementary material

The online version of this article (10.1186/s13104-019-4070-1) contains supplementary material, which is available to authorized users.

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.

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.

Characterization of viral RNA splicing using whole-transcriptome datasets from host species

RNA alternative splicing (AS) is an important post-transcriptional mechanism enabling single genes to produce multiple proteins. It has been well demonstrated that viruses deploy host AS machinery for viral protein productions. However, knowledge on viral AS is limited to a few disease-causing viruses in model species. Here we report a novel approach to characterizing viral AS using whole transcriptome dataset from host species. Two insect transcriptomes (Acheta domesticus and Planococcus citri) generated in the 1,000 Insect Transcriptome Evolution (1KITE) project were used as a proof of concept using the new pipeline. Two closely related densoviruses (Acheta domesticus densovirus, AdDNV, and Planococcus citri densovirus, PcDNV, Ambidensovirus, Densovirinae, Parvoviridae) were detected and analyzed for AS patterns. The results suggested that although the two viruses shared major AS features, dramatic AS divergences were observed. Detailed analysis of the splicing junctions showed clusters of AS events occurred in two regions of the virus genome, demonstrating that transcriptome analysis could gain valuable insights into viral splicing. When applied to large-scale transcriptomics projects with diverse taxonomic sampling, our new method is expected to rapidly expand our knowledge on RNA splicing mechanisms for a wide range of viruses.

EBV and Apoptosis: The Viral Master Regulator of Cell Fate?

Epstein-Barr virus (EBV) was first discovered in cells from a patient with Burkitt lymphoma (BL), and is now known to be a contributory factor in 1-2% of all cancers, for which there are as yet, no EBV-targeted therapies available. Like other herpesviruses, EBV adopts a persistent latent infection in vivo and only rarely reactivates into replicative lytic cycle. Although latency is associated with restricted patterns of gene expression, genes are never expressed in isolation; always in groups. Here, we discuss (1) the ways in which the latent genes of EBV are known to modulate cell death, (2) how these mechanisms relate to growth transformation and lymphomagenesis, and (3) how EBV genes cooperate to coordinately regulate key cell death pathways in BL and lymphoblastoid cell lines (LCLs). Since manipulation of the cell death machinery is critical in EBV pathogenesis, understanding the mechanisms that underpin EBV regulation of apoptosis therefore provides opportunities for novel therapeutic interventions.

A Temporal Proteomic Map of Epstein-Barr Virus Lytic Replication in B Cells

Epstein-Barr virus (EBV) replication contributes to multiple human diseases, including infectious mononucleosis, nasopharyngeal carcinoma, B cell lymphomas, and oral hairy leukoplakia. We performed systematic quantitative analyses of temporal changes in host and EBV proteins during lytic replication to gain insights into virus-host interactions, using conditional Burkitt lymphoma models of type I and II EBV infection. We quantified profiles of >8,000 cellular and 69 EBV proteins, including >500 plasma membrane proteins, providing temporal views of the lytic B cell proteome and EBV virome. Our approach revealed EBV-induced remodeling of cell cycle, innate and adaptive immune pathways, including upregulation of the complement cascade and proteasomal degradation of the B cell receptor complex, conserved between EBV types I and II. Cross-comparison with proteomic analyses of human cytomegalovirus infection and of a Kaposi-sarcoma-associated herpesvirus immunoevasin identified host factors targeted by multiple herpesviruses. Our results provide an important resource for studies of EBV replication.

Current Trends and Alternative Scenarios in EBV Research

Epstein-Barr virus (EBV) infection is associated with several distinct hematological and epithelial malignancies, e.g., Burkitt lymphoma, Hodgkin lymphoma, nasopharyngeal carcinoma, gastric carcinoma, and others. The association with several malignant tumors of local and worldwide distribution makes EBV one of the most important tumor viruses. Furthermore, because EBV can cause posttransplant lymphoproliferative disease, transplant medicine has to deal with EBV as a major pathogenic virus second only to cytomegalovirus. In this review, we summarize briefly the natural history of EBV infection and outline some of the recent advances in the pathogenesis of the major EBV-associated neoplasms. We present alternative scenarios and discuss them in the light of most recent experimental data. Emerging research areas including EBV-induced patho-epigenetic alterations in host cells and the putative role of exosome-mediated information transfer in disease development are also within the scope of this review. This book contains an in-depth description of a series of modern methodologies used in EBV research. In this introductory chapter, we thoroughly refer to the applications of these methods and demonstrate how they contributed to the understanding of EBV-host cell interactions. The data gathered using recent technological advancements in molecular biology and immunology as well as the application of sophisticated in vitro and in vivo experimental models certainly provided deep and novel insights into the pathogenetic mechanisms of EBV infection and EBV-associated tumorigenesis. Furthermore, the development of adoptive T cell immunotherapy has provided a novel approach to the therapy of viral disease in transplant medicine and hematology.

Analysis of EBV Transcription Using High-Throughput RNA Sequencing

High-throughput sequencing of RNA is used to analyze the transcriptomes of viruses and cells, providing information about transcript structure and abundance. A wide array of programs and pipelines has been created to manage and interpret the abundance of data generated from high-throughput RNA sequencing experiments. This protocol details the use of free and open-source programs to align RNA-Seq reads to a reference genome, visualize read coverage and splice junctions, estimate transcript abundance, and evaluate differential expression of transcripts in different conditions. Particular concerns related to EBV and viral transcriptomics are addressed and access to EBV reference files is provided.

MHC class I-associated peptides derive from selective regions of the human genome

MHC class I-associated peptides (MAPs) define the immune self for CD8+ T lymphocytes and are key targets of cancer immunosurveillance. Here, the goals of our work were to determine whether the entire set of protein-coding genes could generate MAPs and whether specific features influence the ability of discrete genes to generate MAPs. Using proteogenomics, we have identified 25,270 MAPs isolated from the B lymphocytes of 18 individuals who collectively expressed 27 high-frequency HLA-A,B allotypes. The entire MAP repertoire presented by these 27 allotypes covered only 10% of the exomic sequences expressed in B lymphocytes. Indeed, 41% of expressed protein-coding genes generated no MAPs, while 59% of genes generated up to 64 MAPs, often derived from adjacent regions and presented by different allotypes. We next identified several features of transcripts and proteins associated with efficient MAP production. From these data, we built a logistic regression model that predicts with good accuracy whether a gene generates MAPs. Our results show preferential selection of MAPs from a limited repertoire of proteins with distinctive features. The notion that the MHC class I immunopeptidome presents only a small fraction of the protein-coding genome for monitoring by the immune system has profound implications in autoimmunity and cancer immunology.

Insights into the role of immunosenescence during varicella zoster virus infection (shingles) in the aging cell model

Varicella zoster virus (VZV) is the etiological agent of shingles, a painful skin rash that affects a significant proportion of the elderly population. In the present study, we used two aging cell models, Hutchinson-Gilford progeria syndrome (HGPS) fibroblasts and stress or replicative senescence-induced normal human dermal fibroblasts (NHDFs), to investigate age-associated susceptibility to VZV infection. VZV infectivity titers were significantly associated with donor age in HGPS fibroblasts and senescence induction in NHDFs. High throughput RNA-sequencing (RNA-seq) analysis was performed to assess global and dynamic changes in the host transcriptomes of VZV-infected aging cells. Analysis of differentially expressed genes (DEGs) indicated that VZV infection in aged HGPS fibroblasts resembled that in senescent NHDFs, particularly in terms of genes associated with pattern recognition receptors in virus sensing network, providing novel insights into the mechanisms of senescence-associated susceptibility to VZV infection. Additionally, we identified stimulator of interferon genes (STING) as a potential VZV sensing receptor. Knockdown of STING expression resulted in increased viral replication in primary fibroblasts, whereas STING overexpression led to suppression of VZV plaque formation. In conclusion, our findings highlight the important role of immunosenescence following VZV infection and provide significant insights into the mechanisms underlying cellular sensing of VZV infection and the induction of immune responses in aged skin cells.

Analysis of Epstein-Barr virus and cellular gene expression during the early phases of Epstein-Barr virus lytic induction

In order to develop novel host/pathogen real-time PCR assays for routine diagnostic use, early gene expression patterns from both Epstein-Barr virus (EBV) and Raji cells were examined after inducing the lytic life cycle using 12-O-tetradecanoyl-13-phorbol ester and sodium butyrate. Real-time PCR identified several highly induced (>90-fold) EBV lytic genes over a 48 h time course during the lytic induction phase. Latent genes were induced at low levels during this phase. The cellular response to lytic viral replication is poorly understood. Whole human genome microarray analysis identified 113 cellular genes regulated twofold or more by EBV, including 63 upregulated and 46 downregulated genes, over a 24 h time course post-induction. The most upregulated gene was CHI3L1, a chitinase-3-like 1 protein (18.1-fold; P<0.0084), and the most downregulated gene was TYMS, a thymidylate synthetase (-7.6-fold). Gene Ontology enrichment analysis using MetaCore software revealed cell cycle (core), cell cycle (role of anaphase-promoting complex) in cell cycle regulation) and lymphatic diseases as the most significantly represented biological network processes, canonical pathways and disease biomarkers, respectively. Chemotaxis, DNA damage and inflammation (IL-4 signalling) together with lymphoproliferative disorders and non-Hodgkin's lymphoma were significantly represented biological processes and disease biomarkers.

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.

Epstein-Barr virus latency type and spontaneous reactivation predict lytic induction levels

The human Epstein-Barr virus (EBV) evades the immune system by entering a transcriptionally latent phase in B cells. EBV in tumor cells expresses distinct patterns of genes referred to as latency types. Viruses in tumor cells also display varying levels of lytic transcription resulting from spontaneous reactivation out of latency. We measured this dynamic range of lytic transcription with RNA deep sequencing and observed no correlation with EBV latency types among genetically different viruses, but type I cell lines reveal more spontaneous reactivation than isogenic type III cultures. We further determined that latency type and spontaneous reactivation levels predict the relative amount of induced reactivation generated by cytotoxic chemotherapy drugs. Our work has potential implications for personalizing medicine against EBV-transformed malignancies. Identifying latency type or measuring spontaneous reactivation may provide predictive power in treatment contexts where viral production should be either avoided or coerced.

EBV Noncoding RNAs

EBV expresses a number of viral noncoding RNAs (ncRNAs) during latent infection, many of which have known regulatory functions and can post-transcriptionally regulate viral and/or cellular gene expression. With recent advances in RNA sequencing technologies, the list of identified EBV ncRNAs continues to grow. EBV-encoded RNAs (EBERs) , the BamHI-A rightward transcripts (BARTs) , a small nucleolar RNA (snoRNA) , and viral microRNAs (miRNAs) are all expressed during EBV infection in a variety of cell types and tumors. Recently, additional novel EBV ncRNAs have been identified. Viral miRNAs, in particular, have been under extensive investigation since their initial identification over ten years ago. High-throughput studies to capture miRNA targets have revealed a number of miRNA-regulated viral and cellular transcripts that tie into important biological networks. Functions for many EBV ncRNAs are still unknown; however, roles for many EBV miRNAs in latency and in tumorigenesis have begun to emerge. Ongoing mechanistic studies to elucidate the functions of EBV ncRNAs should unravel additional roles for ncRNAs in the viral life cycle. In this chapter, we will discuss our current knowledge of the types of ncRNAs expressed by EBV, their potential roles in viral latency, and their potential involvement in viral pathogenesis.

Identification of Epstein-Barr Virus Replication Proteins in Burkitt's Lymphoma Cells

The working model to describe the mechanisms used to replicate the cancer-associated virus Epstein-Barr virus (EBV) is partly derived from comparisons with other members of the Herpes virus family. Many genes within the EBV genome are homologous across the herpes virus family. Published transcriptome data for the EBV genome during its lytic replication cycle show extensive transcription, but the identification of the proteins is limited. We have taken a global proteomics approach to identify viral proteins that are expressed during the EBV lytic replication cycle. We combined an enrichment method to isolate cells undergoing EBV lytic replication with SILAC-labeling coupled to mass-spectrometry and identified viral and host proteins expressed during the EBV lytic replication cycle. Amongst the most frequently identified viral proteins are two components of the DNA replication machinery, the single strand DNA binding protein BALF2, DNA polymerase accessory protein BMRF1 and both subunits of the viral ribonucleoside-diphosphate reductase enzyme (BORF2 and BaRF1). An additional 42 EBV lytic cycle proteins were also detected. This provides proteomic identification for many EBV lytic replication cycle proteins and also identifies post-translational modifications.

EBV glycoproteins: where are we now?

Glycoproteins are critical to virus entry, to spread within and between hosts and can modify the behavior of cells. Many viruses carry only a few, most found in the virion envelope. EBV makes more than 12, providing flexibility in how it colonizes its human host. Some are dedicated to getting the virus through the cell membrane and on toward the nucleus of the cell, some help guide the virus back out and on to the next cell in the same or a new host. Yet others undermine host defenses helping the virus persist for a lifetime, maintaining a presence that is mostly tolerated and serves to perpetuate EBV as one of the most common infections of man.

Host Gene Expression Is Regulated by Two Types of Noncoding RNAs Transcribed from the Epstein-Barr Virus BamHI A Rightward Transcript Region

In Epstein-Barr virus-infected epithelial cancers, the alternatively spliced BamHI A rightward transcripts (BARTs) are the most abundant viral polyadenylated RNA. The BART introns form the template for the production of 44 microRNAs (miRNAs), and the spliced and polyadenylated exons form nuclear non-protein-coding RNAs. Analysis of host cell transcription by RNA-seq during latency in AGS cells identified a large number of reproducibly changed genes. Genes that were downregulated were enriched for BART miRNA targets. Bioinformatics analysis predicted activation of the myc pathway and downregulation of XBP1 as likely mediators of the host transcriptional changes. Effects on XBP1 activity were not detected in these cells; however, myc activation was confirmed through use of a myc-responsive luciferase reporter. To identify potential regulatory properties of the spliced, polyadenylated BART RNAs, a full-length cDNA clone of one of the BART isoforms was obtained and expressed in the Epstein-Barr virus (EBV)-negative AGS cells. The BART cDNA transcript remained primarily nuclear yet induced considerable and consistent changes in cellular transcription, as profiled by RNA-seq. These transcriptional changes significantly overlapped the transcriptional changes induced during latent EBV infection of these same cells, where the BARTs are exclusively nuclear and do not encode proteins. These data suggest that the nuclear BART RNAs are functional long noncoding RNAs (lncRNAs). The abundant expression of multiple forms of noncoding RNAs that contribute to growth regulation without expression of immunogenic proteins would be an important mechanism for viral oncogenesis in the presence of a functional immune system.

IMPORTANCE

Infection with Epstein-Barr virus (EBV) is nearly ubiquitous in the human population; however, it does contribute to the formation of multiple types of cancer. In immunocompromised patients, EBV causes multiple types of lymphomas by expressing viral oncogenes that promote growth and survival of infected B lymphocytes. EBV-positive gastric carcinoma does not require immune suppression, and the viral oncoproteins that are frequent targets for an immunological response are not expressed. This study demonstrates using transcriptional analysis that the expression of various classes of viral non-protein-coding RNAs likely contribute to the considerable changes in the host transcriptional profile in the AGS gastric cancer cell line. This is the first report to show that the highly expressed polyadenylated BamHI A rightward transcripts (BART) viral transcript in gastric carcinoma is in fact a functional viral long noncoding RNA. These studies provide new insight into how EBV can promote transformation in the absence of viral protein expression.

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.

Genome-wide analysis of wild-type Epstein-Barr virus genomes derived from healthy individuals of the 1,000 Genomes Project

Most people in the world (∼90%) are infected by the Epstein–Barr virus (EBV), which establishes itself permanently in B cells. Infection by EBV is related to a number of diseases including infectious mononucleosis, multiple sclerosis, and different types of cancer. So far, only seven complete EBV strains have been described, all of them coming from donors presenting EBV-related diseases. To perform a detailed comparative genomic analysis of EBV including, for the first time, EBV strains derived from healthy individuals, we reconstructed EBV sequences infecting lymphoblastoid cell lines (LCLs) from the 1000 Genomes Project. As strain B95-8 was used to transform B cells to obtain LCLs, it is always present, but a specific deletion in its genome sets it apart from natural EBV strains. After studying hundreds of individuals, we determined the presence of natural EBV in at least 10 of them and obtained a set of variants specific to wild-type EBV. By mapping the natural EBV reads into the EBV reference genome (NC007605), we constructed nearly complete wild-type viral genomes from three individuals. Adding them to the five disease-derived EBV genomic sequences available in the literature, we performed an in-depth comparative genomic analysis. We found that latency genes harbor more nucleotide diversity than lytic genes and that six out of nine latency-related genes, as well as other genes involved in viral attachment and entry into host cells, packaging, and the capsid, present the molecular signature of accelerated protein evolution rates, suggesting rapid host–parasite coevolution.

Latent Membrane Protein 2 (LMP2)

LMP2A is an EBV-encoded protein with three domains: (a) an N-terminal cytoplasmic domain, which has PY motifs that bind to WW domain-containing E3 ubiquitin ligases and an ITAM that binds to SH2 domain-containing proteins, (b) a transmembrane domain with 12 transmembrane segments that localizes LMP2A in cellular membranes, and (c) a 27-amino acid C-terminal domain which mediates homodimerization and heterodimerization of LMP2 protein isoforms. The most prominent two isoforms of the protein are LMP2A and LMP2B. The LMP2B isoform lacks the 19-amino acid N-terminal domain found in LMP2A, which modulates cellular signaling resulting in a baseline activation of B cells and degradation of cellular kinases leading to the downregulation of normal B cell signaling pathways. These two seemingly contradictory processes allow EBV to establish and maintain latency. LMP2 is expressed in many EBV-associated malignancies. While its antigenic properties may be useful in developing LMP2-specific immunity, the LMP2A N-terminal motifs also provide a basis to target LMP2A-modulated cellular kinases for the development of treatment strategies.

Unexpected patterns of Epstein-Barr virus transcription revealed by a high throughput PCR array for absolute quantification of viral mRNA

We have validated a flexible, high-throughput and relatively inexpensive RT-QPCR array platform for absolute quantification of Epstein-Barr virus transcripts in different latent and lytic infection states. Several novel observations are reported. First, during infection of normal B cells, Wp-initiated latent gene transcripts remain far more abundant following activation of the Cp promoter than was hitherto suspected. Second, EBNA1 transcript levels are remarkably low in all forms of latency, typically ranging from 1 to 10 transcripts per cell. EBNA3A, -3B and -3C transcripts are likewise very low in Latency III, typically at levels similar to or less than EBNA1 transcripts. Thirdly, a subset of lytic gene transcripts is detectable in Burkitt lymphoma lines at low levels, including: BILF1, which has oncogenic properties, and the poorly characterized LF1, LF2 and LF3 genes. Analysis of seven African BL biopsies confirmed this transcription profile but additionally revealed significant expression of LMP2 transcripts.

Oncogenes and RNA splicing of human tumor viruses

Approximately 10.8% of human cancers are associated with infection by an oncogenic virus. These viruses include human papillomavirus (HPV), Epstein–Barr virus (EBV), Merkel cell polyomavirus (MCV), human T-cell leukemia virus 1 (HTLV-1), Kaposi's sarcoma-associated herpesvirus (KSHV), hepatitis C virus (HCV) and hepatitis B virus (HBV). These oncogenic viruses, with the exception of HCV, require the host RNA splicing machinery in order to exercise their oncogenic activities, a strategy that allows the viruses to efficiently export and stabilize viral RNA and to produce spliced RNA isoforms from a bicistronic or polycistronic RNA transcript for efficient protein translation. Infection with a tumor virus affects the expression of host genes, including host RNA splicing factors, which play a key role in regulating viral RNA splicing of oncogene transcripts. A current prospective focus is to explore how alternative RNA splicing and the expression of viral oncogenes take place in a cell- or tissue-specific manner in virus-induced human carcinogenesis.

Pervasive transcription of a herpesvirus genome generates functionally important RNAs

Pervasive transcription is observed in a wide range of organisms, including humans, mice, and viruses, but the functional significance of the resulting transcripts remains uncertain. Current genetic approaches are often limited by their emphasis on protein-coding open reading frames (ORFs). We previously identified extensive pervasive transcription from the murine gammaherpesvirus 68 (MHV68) genome outside known ORFs and antisense to known genes (termed expressed genomic regions [EGRs]). Similar antisense transcripts have been identified in many other herpesviruses, including Kaposi’s sarcoma-associated herpesvirus and human and murine cytomegalovirus. Despite their prevalence, whether these RNAs have any functional importance in the viral life cycle is unknown, and one interpretation is that these are merely “noise” generated by functionally unimportant transcriptional events. To determine whether pervasive transcription of a herpesvirus genome generates RNA molecules that are functionally important, we used a strand-specific functional approach to target transcripts from thirteen EGRs in MHV68. We found that targeting transcripts from six EGRs reduced viral protein expression, proving that pervasive transcription can generate functionally important RNAs. We characterized transcripts emanating from EGRs 26 and 27 in detail using several methods, including RNA sequencing, and identified several novel polyadenylated transcripts that were enriched in the nuclei of infected cells. These data provide the first evidence of the functional importance of regions of pervasive transcription emanating from MHV68 EGRs. Therefore, studies utilizing mutation of a herpesvirus genome must account for possible effects on RNAs generated by pervasive transcription.

The fact that pervasive transcription produces functionally important RNAs has profound implications for design and interpretation of genetic studies in herpesviruses, since such studies often involve mutating both strands of the genome. This is a common potential problem; for example, a conservative estimate is that there are an additional 73,000 nucleotides transcribed antisense to annotated ORFs from the 119,450-bp MHV68 genome. Recognizing the importance of considering the function of each strand of the viral genome independently, we used strand-specific approaches to identify six regions of the genome encoding transcripts that promoted viral protein expression. For two of these regions, we mapped novel transcripts and determined that targeting transcripts from these regions reduced viral replication and the expression of other viral genes. This is the first description of a function for these RNAs and suggests that novel transcripts emanating from regions of pervasive transcription are critical for the viral life cycle.

RNA CoMPASS: a dual approach for pathogen and host transcriptome analysis of RNA-seq datasets

High-throughput RNA sequencing (RNA-seq) has become an instrumental assay for the analysis of multiple aspects of an organism's transcriptome. Further, the analysis of a biological specimen's associated microbiome can also be performed using RNA-seq data and this application is gaining interest in the scientific community. There are many existing bioinformatics tools designed for analysis and visualization of transcriptome data. Despite the availability of an array of next generation sequencing (NGS) analysis tools, the analysis of RNA-seq data sets poses a challenge for many biomedical researchers who are not familiar with command-line tools. Here we present RNA CoMPASS, a comprehensive RNA-seq analysis pipeline for the simultaneous analysis of transcriptomes and metatranscriptomes from diverse biological specimens. RNA CoMPASS leverages existing tools and parallel computing technology to facilitate the analysis of even very large datasets. RNA CoMPASS has a web-based graphical user interface with intrinsic queuing to control a distributed computational pipeline. RNA CoMPASS was evaluated by analyzing RNA-seq data sets from 45 B-cell samples. Twenty-two of these samples were derived from lymphoblastoid cell lines (LCLs) generated by the infection of naïve B-cells with the Epstein Barr virus (EBV), while another 23 samples were derived from Burkitt's lymphomas (BL), some of which arose in part through infection with EBV. Appropriately, RNA CoMPASS identified EBV in all LCLs and in a fraction of the BLs. Cluster analysis of the human transcriptome component of the RNA CoMPASS output clearly separated the BLs (which have a germinal center-like phenotype) from the LCLs (which have a blast-like phenotype) with evidence of activated MYC signaling and lower interferon and NF-kB signaling in the BLs. Together, this analysis illustrates the utility of RNA CoMPASS in the simultaneous analysis of transcriptome and metatranscriptome data. RNA CoMPASS is freely available at http://rnacompass.sourceforge.net/.

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.

Global bidirectional transcription of the Epstein-Barr virus genome during reactivation

Epstein-Barr virus (EBV) reactivation involves the ordered induction of approximately 90 viral genes that participate in the generation of infectious virions. Using strand-specific RNA-seq to assess the EBV transcriptome during reactivation, we found extensive bidirectional transcription extending across nearly the entire genome. In contrast, only 4% of the EBV genome is currently bidirectionally annotated. Most of the newly identified transcribed regions show little evidence of coding potential, supporting noncoding roles for most of these RNAs. Based on previous cellular long noncoding RNA size calculations, we estimate that there are likely hundreds more EBV genes expressed during reactivation than was previously known. Limited 5' and 3' rapid amplification of cDNA ends (RACE) experiments and findings of novel splicing events by RNA-seq suggest that the complexity of the viral genome during reactivation may be even greater. Further analysis of antisense transcripts at some of the EBV latency gene loci showed that they are "late" genes, they are nuclear, and they tend to localize in areas of the nucleus where others find newly synthesized viral genomes. This raises the possibility that these transcripts perform functions such as new genome processing, stabilization, organization, etc. The finding of a significantly more complex EBV transcriptome during reactivation changes our view of the viral production process from one that is facilitated and regulated almost entirely by previously identified viral proteins to a process that also involves the contribution of a wide array of virus encoded noncoding RNAs. Epstein-Barr virus (EBV) is a herpesvirus that infects the majority of the world's population, in rare cases causing serious disease such as lymphoma and gastric carcinoma. Using strand-specific RNA-seq, we have studied viral gene expression during EBV reactivation and have discovered hundreds more viral transcripts than were previously known. The finding of alternative splicing and the prevalence of overlapping transcripts indicate additional complexity. Most newly identified transcribed regions do not encode proteins but instead likely function as noncoding RNA molecules which could participate in regulating gene expression, gene splicing or even activities such as viral genome processing. These findings broaden the scope of what we need to consider to understand the viral manufacturing process. As more detailed studies are undertaken they will likely change the way we view this process as a whole.

A gammaherpesvirus uses alternative splicing to regulate its tropism and its sensitivity to neutralization

Human gammaherpesviruses are associated with the development of lymphomas and epithelial malignancies. The heterogeneity of these tumors reflects the ability of these viruses to route infection to different cell types at various stages of their lifecycle. While the Epstein Barr virus uses gp42 – human leukocyte antigen class II interaction as a switch of cell tropism, the molecular mechanism that orientates tropism of rhadinoviruses is still poorly defined. Here, we used bovine herpesvirus 4 (BoHV-4) to further elucidate how rhadinoviruses regulate their infectivity. In the absence of any gp42 homolog, BoHV-4 exploits the alternative splicing of its Bo10 gene to produce distinct viral populations that behave differently based on the originating cell. While epithelial cells produce virions with high levels of the accessory envelope protein gp180, encoded by a Bo10 spliced product, myeloid cells express reduced levels of gp180. As a consequence, virions grown in epithelial cells are hardly infectious for CD14+ circulating cells, but are relatively resistant to antibody neutralization due to the shielding property of gp180 for vulnerable entry epitopes. In contrast, myeloid virions readily infect CD14+ circulating cells but are easily neutralized. This molecular switch could therefore allow BoHV-4 to promote either, on the one hand, its dissemination into the organism, or, on the other hand, its transmission between hosts.

Gammaherpesviruses are highly prevalent human and animal pathogens. These viruses display sophisticated entry mechanisms, allowing them to infect different cell types inside a host but also to transmit between hosts in the presence of neutralizing antibodies. Here, we used bovine herpesvirus 4 (BoHV-4) to decipher how some gammaherpesviruses manage to do this. We found that, as function of the originating cell types, BoHV-4 is able to modify its tropism as well as its sensitivity to antibody neutralization just by controlling the alternative splicing of one of its genes. This virus therefore exploits post-transcriptional events to generate viral populations with distinct phenotypes.

Genome-wide analyses of Epstein-Barr virus reveal conserved RNA structures and a novel stable intronic sequence RNA

BACKGROUND

Epstein-Barr virus (EBV) is a human herpesvirus implicated in cancer and autoimmune disorders. Little is known concerning the roles of RNA structure in this important human pathogen. This study provides the first comprehensive genome-wide survey of RNA and RNA structure in EBV.

RESULTS

Novel EBV RNAs and RNA structures were identified by computational modeling and RNA-Seq analyses of EBV. Scans of the genomic sequences of four EBV strains (EBV-1, EBV-2, GD1, and GD2) and of the closely related Macacine herpesvirus 4 using the RNAz program discovered 265 regions with high probability of forming conserved RNA structures. Secondary structure models are proposed for these regions based on a combination of free energy minimization and comparative sequence analysis. The analysis of RNA-Seq data uncovered the first observation of a stable intronic sequence RNA (sisRNA) in EBV. The abundance of this sisRNA rivals that of the well-known and highly expressed EBV-encoded non-coding RNAs (EBERs).

CONCLUSION

This work identifies regions of the EBV genome likely to generate functional RNAs and RNA structures, provides structural models for these regions, and discusses potential functions suggested by the modeled structures. Enhanced understanding of the EBV transcriptome will guide future experimental analyses of the discovered RNAs and RNA structures.

Complexities of gammaherpesvirus transcription revealed by microarrays and RNAseq

Technological advances in genome-wide transcript analysis, referred to as the transcriptome, using microarrays and deep RNA sequencing methodologies are rapidly extending our understanding of the genetic content of the gammaherpesviruses (γHVs). These vast transcript analyses continue to uncover the complexity of coding transcripts due to alternative splicing, translation initiation and termination, as well as regulatory RNAs of the γHVs. A full assessment of the transcriptome requires that our analysis be extended to the virion and exosomes of infected cells since viral and host mRNAs, miRNAs, and other noncoding RNAs seem purposefully incorporated to exert function upon delivery to naïve cells. Understanding the regulation, biogenesis and function of the recently discovered transcripts will extend beyond pathogenesis and oncogenic events to offer key insights for basic RNA processes of the cell.

Differences in gastric carcinoma microenvironment stratify according to EBV infection intensity: implications for possible immune adjuvant therapy

Epstein-Barr virus (EBV) is associated with roughly 10% of gastric carcinomas worldwide (EBVaGC). Although previous investigations provide a strong link between EBV and gastric carcinomas, these studies were performed using selected EBV gene probes. Using a cohort of gastric carcinoma RNA-seq data sets from The Cancer Genome Atlas (TCGA), we performed a quantitative and global assessment of EBV gene expression in gastric carcinomas and assessed EBV associated cellular pathway alterations. EBV transcripts were detected in 17% of samples but these samples varied significantly in EBV coverage depth. In four samples with the highest EBV coverage (hiEBVaGC – high EBV associated gastric carcinoma), transcripts from the BamHI A region comprised the majority of EBV reads. Expression of LMP2, and to a lesser extent, LMP1 were also observed as was evidence of abortive lytic replication. Analysis of cellular gene expression indicated significant immune cell infiltration and a predominant IFNG response in samples expressing high levels of EBV transcripts relative to samples expressing low or no EBV transcripts. Despite the apparent immune cell infiltration, high levels of the cytotoxic T-cell (CTL) and natural killer (NK) cell inhibitor, IDO1, was observed in the hiEBVaGCs samples suggesting an active tolerance inducing pathway in this subgroup. These results were confirmed in a separate cohort of 21 Vietnamese gastric carcinoma samples using qRT-PCR and on tissue samples using in situ hybridization and immunohistochemistry. Lastly, a panel of tumor suppressors and candidate oncogenes were expressed at lower levels in hiEBVaGC versus EBV-low and EBV-negative gastric cancers suggesting the direct regulation of tumor pathways by EBV.

Epstein-Barr virus (EBV) is detected in roughly 10% of gastric carcinoma (GC) cases worldwide. Despite a strong link between EBV and gastric carcinoma, the contribution of EBV to the tumor environment in EBV associated gastric carcinoma is unclear. We performed a global assessment of EBV and host cell gene expression in gastric carcinoma tumors from 71 patients to link EBV genes (and expression intensities) to cell and microenvironmental changes. In addition to the finding that EBV is associated with down-regulated tumor regulatory genes, this study revealed that samples with high levels of EBV gene expression (hiEBVaGCs) displayed elevated immune cell infiltration with high interferon-gamma (IFNG) expression compared to samples with low or no EBV gene expression. Despite this evidence of increased immune posturing, hiEBVaGC samples also showed elevated expression of the potent immune cell inhibitor, IDO1. This finding may partly explain the persistence of these virus associated tumors in the face of local immune cell concentration. Importantly, the small molecule IDO inhibitor, 1MT (1-methyl Tryptophan), has been shown to reverse the tolerance inducing effects of IDO1 in other tumors. We propose that stratification of gastric carcinomas into EBV-negative, EBV-low and EBV-high may provide indicator value for the use of IDO1 inhibitors as adjuvant therapies against hiEBVaGCs.

Next-generation sequencing technologies in diagnostic virology

The data deluge produced by next-generation sequencing (NGS) technologies is an appealing feature for clinical virologists that are involved in the diagnosis of emerging viral infections, molecular epidemiology of viral pathogens, drug-resistance testing, and also like to do some basic and clinical research. Indeed, NGS platforms are being implemented in many clinical and research laboratories, as the costs of these platforms are progressively decreasing. We provide here some suggestions for virologists who are planning to implement a NGS platform in their clinical laboratory and an overview on the potential applications of these technologies in diagnostic virology.

Whole-genome sequencing of the Akata and Mutu Epstein-Barr virus strains

Using a simple viral genome enrichment approach, we report the de novo assembly of the Akata and Mutu Epstein-Barr virus (EBV) genomes from a single lane of next-generation sequencing (NGS) reads. The Akata and Mutu viral genomes are type I EBV strains of approximately 171 kb in length. Evidence for genome heterogeneity was found for the Akata but not for the Mutu strain. A comparative analysis of Akata with another four completely sequenced EBV strains, B95-8/Raji, AG876, Mutu, and GD1, demonstrated that the Akata strain is most closely related to the GD1 strain and exhibits the greatest divergence from the type II strain, AG876. A global comparison of latent and lytic gene sequences showed that the four latency genes, EBNA2, EBNA3A, EBNA3B, and EBNA3C, are uniquely defining of type I and type II strain differences. Within type I strains, LMP1, the latency gene, is among the most divergent of all EBV genes, with three insertion or deletion loci in its CTAR2 and CTAR3 signaling domains. Analysis of the BHLF1 and LF3 genes showed that the reading frames identified in the B95-8/Raji genome are not conserved in Akata (or Mutu, for BHLF1), suggesting a primarily non-protein-coding function in EBV's life cycle. The Akata and Mutu viral-genome sequences should be a useful resource for homology-based functional prediction and for molecular studies, such as PCR, RNA-seq, recombineering, and transcriptome studies. As an illustration, we identified novel RNA-editing events in ebv-miR-BART6 antisense transcripts using the Akata and Mutu reference genomes.

Alternative splicing in human tumour viruses: a therapeutic target?

Persistent infection with cancer risk-related viruses leads to molecular, cellular and immune response changes in host organisms that in some cases direct cellular transformation. Alternative splicing is a conserved cellular process that increases the coding complexity of genomes at the pre-mRNA processing stage. Human and other animal tumour viruses use alternative splicing as a process to maximize their transcriptomes and proteomes. Medical therapeutics to clear persistent viral infections are still limited. However, specific lessons learned in some viruses [e.g. HIV and HCV (hepatitis C virus)] suggest that drug-directed inhibition of alternative splicing could be useful for this purpose. The present review describes the basic mechanisms of constitutive and alternative splicing in a cellular context and known splicing patterns and the mechanisms by which these might be achieved for the major human infective tumour viruses. The roles of splicing-related proteins expressed by these viruses in cellular and viral gene regulation are explored. Moreover, we discuss some currently available drugs targeting SR (serine/arginine-rich) proteins that are the main regulators of constitutive and alternative splicing, and their potential use in treatment for so-called persistent viral infections.

Global impact of RNA splicing on transcriptome remodeling in the heart

In the eukaryotic transcriptome, both the numbers of genes and different RNA species produced by each gene contribute to the overall complexity. These RNA species are generated by the utilization of different transcriptional initiation or termination sites, or more commonly, from different messenger RNA (mRNA) splicing events. Among the 30,000+ genes in human genome, it is estimated that more than 95% of them can generate more than one gene product via alternative RNA splicing. The protein products generated from different RNA splicing variants can have different intracellular localization, activity, or tissue-distribution. Therefore, alternative RNA splicing is an important molecular process that contributes to the overall complexity of the genome and the functional specificity and diversity among different cell types. In this review, we will discuss current efforts to unravel the full complexity of the cardiac transcriptome using a deep-sequencing approach, and highlight the potential of this technology to uncover the global impact of RNA splicing on the transcriptome during development and diseases of the heart.

Human papillomavirus gene expression is controlled by host cell splicing factors

HPVs (human papillomaviruses) infect stratified epithelia and cause a variety of lesions ranging from benign warts to invasive tumours. The virus life cycle is tightly linked to differentiation of the keratinocyte it infects: papillomaviruses modulate host gene expression to ensure efficient virus replication. For example, the viral transcription factor E2 can directly up-regulate, in an epithelial differentiation-dependent manner, cellular SRSFs [SR (serine/arginine-rich) splicing factors] that control constitutive and alternative splicing. Changes in alternative splicing and the mechanisms controlling this for viral mRNAs have been the subject of intense exploration. However, to date experiments have only been carried out in model systems because the genetic systems suitable for studying alternative splicing of viral RNAs in the context of the virus life cycle are relatively recent and technically challenging. Now using these life cycle-supporting systems, our laboratory has identified SR proteins as important players in differentiation-dependent regulation of HPV gene expression. Better understanding of the role of cellular factors in regulating the virus life cycle is needed as it may help development of novel diagnostic approaches and antiviral therapies in the future.