Updated Epstein-Barr virus (EBV) DNA sequence and analysis of a promoter for the BART (CST, BARF0) RNAs of EBV

Exploring the anti-EBV potential of suberoylanilide hydroxamic acid: Induction of apoptosis in infected cells through suppressing BART gene expression and inducing lytic infection

Epstein-Barr virus (EBV) is linked to lymphoma and epithelioma but lacks drugs specifically targeting EBV-positive tumors. BamHI A Rightward Transcript (BART) miRNAs are expressed in all EBV-positive tumors, suppressing both lytic infection and host cell apoptosis. We identified suberoylanilide hydroxamic acid (SAHA), an inhibitor of histone deacetylase enzymes, as an agent that suppresses BART promoter activity and transcription of BART miRNAs. SAHA treatment demonstrated a more pronounced inhibition of cell proliferation in EBV-positive cells compared to EBV-negative cells, affecting both p53 wild-type and mutant gastric epithelial cells. SAHA treatment enhanced lytic infection in wild-type EBV-infected cells, while also enhancing cell death in BZLF1-deficient EBV-infected cells. It reduced BART gene expression by 85% and increased the expression of proapoptotic factors targeted by BART miRNAs. These findings suggest that SAHA not only induces lytic infection but also leads to cell death by suppressing BART miRNA transcription and promoting the apoptotic program.

Epstein-Barr Viruses: Their Immune Evasion Strategies and Implications for Autoimmune Diseases

Epstein-Barr virus (EBV), a member of the γ-herpesvirus family, is one of the most prevalent and persistent human viruses, infecting up to 90% of the adult population globally. EBV's life cycle includes primary infection, latency, and lytic reactivation, with the virus primarily infecting B cells and epithelial cells. This virus has evolved sophisticated strategies to evade both innate and adaptive immune responses, thereby maintaining a lifelong presence within the host. This persistence is facilitated by the expression of latent genes such as EBV nuclear antigens (EBNAs) and latent membrane proteins (LMPs), which play crucial roles in viral latency and oncogenesis. In addition to their well-known roles in several types of cancer, including nasopharyngeal carcinoma and B-cell lymphomas, recent studies have identified the pathogenic roles of EBV in autoimmune diseases such as multiple sclerosis, rheumatoid arthritis, and systemic lupus erythematosus. This review highlights the intricate interactions between EBV and the host immune system, underscoring the need for further research to develop effective therapeutic and preventive strategies against EBV-associated diseases.

Epstein-Barr virus evades restrictive host chromatin closure by subverting B cell activation and germinal center regulatory loci

Chromatin accessibility fundamentally governs gene expression and biological response programs that can be manipulated by pathogens. Here we capture dynamic chromatin landscapes of individual B cells during Epstein-Barr virus (EBV) infection. EBV+ cells that exhibit arrest via antiviral sensing and proliferation-linked DNA damage experience global accessibility reduction. Proliferative EBV+ cells develop expression-linked architectures and motif accessibility profiles resembling in vivo germinal center (GC) phenotypes. Remarkably, EBV elicits dark zone (DZ), light zone (LZ), and post-GC B cell chromatin features despite BCL6 downregulation. Integration of single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq), single-cell RNA sequencing (scRNA-seq), and chromatin immunoprecipitation sequencing (ChIP-seq) data enables genome-wide cis-regulatory predictions implicating EBV nuclear antigens (EBNAs) in phenotype-specific control of GC B cell activation, survival, and immune evasion. Knockouts validate bioinformatically identified regulators (MEF2C and NFE2L2) of EBV-induced GC phenotypes and EBNA-associated loci that regulate gene expression (CD274/PD-L1). These data and methods can inform high-resolution investigations of EBV-host interactions, B cell fates, and virus-mediated lymphomagenesis.

Epstein-Barr Virus and the Pathogenesis of Diffuse Large B-Cell Lymphoma

Epstein-Barr virus (EBV), defined as a group I carcinogen by the World Health Organization (WHO), is present in the tumour cells of patients with different forms of B-cell lymphoma, including Burkitt lymphoma, Hodgkin lymphoma, post-transplant lymphoproliferative disorders, and, most recently, diffuse large B-cell lymphoma (DLBCL). Understanding how EBV contributes to the development of these different types of B-cell lymphoma has not only provided fundamental insights into the underlying mechanisms of viral oncogenesis, but has also highlighted potential new therapeutic opportunities. In this review, we describe the effects of EBV infection in normal B-cells and we address the germinal centre model of infection and how this can lead to lymphoma in some instances. We then explore the recent reclassification of EBV+ DLBCL as an established entity in the WHO fifth edition and ICC 2022 classifications, emphasising the unique nature of this entity. To that end, we also explore the unique genetic background of this entity and briefly discuss the potential role of the tumour microenvironment in lymphomagenesis and disease progression. Despite the recent progress in elucidating the mechanisms of this malignancy, much work remains to be done to improve patient stratification, treatment strategies, and outcomes.

Wildebeest-Derived Malignant Catarrhal Fever: A Bovine Peripheral T Cell Lymphoma Caused by Cross-Species Transmission of Alcelaphine Gammaherpesvirus 1

Gammaherpesviruses (γHVs) include viruses that can induce lymphoproliferative diseases and tumors. These viruses can persist in the long term in the absence of any pathological manifestation in their natural host. Alcelaphine gammaherpesvirus 1 (AlHV-1) belongs to the genus Macavirus and asymptomatically infects its natural host, the wildebeest (Connochaetes spp.). However, when transmitted to several susceptible species belonging to the order Artiodactyla, AlHV-1 is responsible for the induction of a lethal lymphoproliferative disease, named wildebeest-derived malignant catarrhal fever (WD-MCF). Understanding the pathogenic mechanisms responsible for the induction of WD-MCF is important to better control the risks of transmission and disease development in susceptible species. The aim of this review is to synthesize the current knowledge on WD-MCF with a particular focus on the mechanisms by which AlHV-1 induces the disease. We discuss the potential mechanisms of pathogenesis from viral entry into the host to the maintenance of viral genomes in infected CD8⁺ T lymphocytes, and we present current hypotheses to explain how AlHV-1 infection induces a peripheral T cell lymphoma-like disease.

Sequence analysis of Epstein-Barr virus RPMS1 gene in malignant hematopathy of Northern China

The RPMS1 gene is the only member of the BamHI-A rightward transcripts (BARTs) family for which a full-length complementary DNA has been identified, and RPMS1 transcript has been confirmed in many Epstein-Barr virus (EBV)-positive malignancies. However, the effects of sequence variations of RPMS1 in hematological malignancies and their biological significance are unclear. To explore the association between RPMS1 gene variations and hematological malignancy, the RPMS1 gene of 391 EBV-positive samples from patients with EBV-positive leukemia, myelodysplastic syndromes and lymphoma in northern China were sequenced. On the basis of phylogenetic tree and mutation characteristics of RPMS1, all the sequences were divided into five major types: RPMS1-A, RPMS1-B, RPMS1-C, RPMS1-E, and RPMS1-F. RPMS1-A type, similar to the prototype B95-8, was identified in 71.87% (281/391) of samples and was the major type in all subpopulations. The frequency of RPMS1-F type was significantly higher in all malignant hematopathy groups than in healthy donors. The Hodgkin lymphoma group contained more RPMS1-F than other malignant hematopathy groups, and acute myeloid leukemia contained more RPMS1-C type than other malignant hematopathy groups. Therefore, RPMS1-A is the main type of RPMS1 gene in northern China, and RPMS1-F may be associated with hematologic malignancies.

Accurate reconstruction of viral genomes in human cells from short reads using iterative refinement

BACKGROUND

After an infection, human cells may contain viral genomes in the form of episomes or integrated DNA. Comparing the genomic sequences of different strains of a virus in human cells can often provide useful insights into its behaviour, activity and pathology, and may help develop methods for disease prevention and treatment. To support such comparative analyses, the viral genomes need to be accurately reconstructed from a large number of samples. Previous efforts either rely on customized experimental protocols or require high similarity between the sequenced genomes and a reference, both of which limit the general applicability of these approaches. In this study, we propose a pipeline, named ASPIRE, for reconstructing viral genomes accurately from short reads data of human samples, which are increasingly available from genome projects and personal genomics. ASPIRE contains a basic part that involves de novo assembly, tiling and gap filling, and additional components for iterative refinement, sequence corrections and wrapping.

RESULTS

Evaluated by the alignment quality of sequencing reads to the reconstructed genomes, these additional components improve the assembly quality in general, and in some particular samples quite substantially, especially when the sequenced genome is significantly different from the reference. We use ASPIRE to reconstruct the genomes of Epstein Barr Virus (EBV) from the whole-genome sequencing data of 61 nasopharyngeal carcinoma (NPC) samples and provide these sequences as a resource for EBV research.

CONCLUSIONS

ASPIRE improves the quality of the reconstructed EBV genomes in published studies and outperforms TRACESPipe in some samples considered.

Epstein-Barr Virus in Burkitt Lymphoma in Africa Reveals a Limited Set of Whole Genome and LMP-1 Sequence Patterns: Analysis of Archival Datasets and Field Samples From Uganda, Tanzania, and Kenya

Epstein-Barr virus (EBV) is associated with endemic Burkitt lymphoma (eBL), but the contribution of EBV variants is ill-defined. Studies of EBV whole genome sequences (WGS) have identified phylogroups that appear to be distinct for Asian versus non-Asian EBV, but samples from BL or Africa, where EBV was first discovered, are under-represented. We conducted a phylogenetic analysis of EBV WGS and LMP-1 sequences obtained primarily from BL patients in Africa and representative non-African EBV from other conditions or regions using data from GenBank, Sequence Read Archive, or Genomic Data Commons for the Burkitt Lymphoma Genome Sequencing Project (BLGSP) to generate data to support the use of a simpler biomarker of geographic or phenotypic associations. We also investigated LMP-1 patterns in 414 eBL cases and 414 geographically matched controls in the Epidemiology of Burkitt Lymphoma in East African children and minors (EMBLEM) study using LMP-1 PCR and Sanger sequencing. Phylogenetic analysis revealed distinct genetic patterns of African versus Asian EBV sequences. We identified 281 single nucleotide variations (SNVs) in LMP-1 promoter and coding region, which formed 12 unique patterns (A to L). Nine patterns (A, AB, C, D, F, I, J, K and L) predominated in African EBV, of which four were found in 92% of BL samples (A, AB, D, and H). Predominant patterns were B and G in Asia and H in Europe. EBV positivity in peripheral blood was detected in 95.6% of EMBLEM eBL cases versus 79.2% of the healthy controls (odds ratio [OR] =3.83; 95% confidence interval 2.06-7.14). LMP-1 was successfully sequenced in 66.7% of the EBV DNA positive cases but in 29.6% of the controls (ORs ranging 5-11 for different patterns). Four LMP-1 patterns (A, AB, D, and K) were detected in 63.1% of the cases versus 27.1% controls (ORs ranges: 5.58-11.4). Dual strain EBV infections were identified in WGS and PCR-Sanger data. In conclusion, EBV from Africa is phylogenetically separate from EBV in Asia. Genetic diversity in LMP-1 formed 12 patterns, which showed promising geographic and phenotypic associations. Presence of multiple strain infection should be considered in efforts to refine or improve EBV markers of ancestry or phenotype.

Lay Summary

Epstein-Barr virus (EBV) infection, a ubiquitous infection, contributes to the etiology of both Burkitt Lymphoma (BL) and nasopharyngeal carcinoma, yet their global distributions vary geographically with no overlap. Genomic variation in EBV is suspected to play a role in the geographical patterns of these EBV-associated cancers, but relatively few EBV samples from BL have been comprehensively studied. We sought to compare phylogenetic patterns of EBV genomes obtained from BL samples in Africa and from tumor and non-tumor samples from elsewhere. We concluded that EBV obtained from BL in Africa is genetically separate from EBV in Asia. Through comprehensive analysis of nucleotide variations in EBV's LMP-1 gene, we describe 12 LMP-1 patterns, two of which (B and G) were found mostly in Asia. Four LMP-1 patterns (A, AB, D, and F) accounted for 92% of EBVs sequenced from BL in Africa. Our results identified extensive diversity of EBV, but BL in Africa was associated with a limited number of variants identified, which were different from those identified in Asia. Further research is needed to optimize the use of PCR and sequencing to study LMP-1 diversity for classification of EBV variants and for use in epidemiologic studies to characterize geographic and/or phenotypic associations of EBV variants with EBV-associated malignancies, including eBL.

Genome-Wide Profiling of Epstein-Barr Virus (EBV) Isolated from EBV-Related Malignancies

Epstein–Barr virus (EBV) is the cause of certain cancers, such as Burkitt lymphoma, Hodgkin lymphoma, NK/T cell lymphoma, nasopharyngeal carcinoma, and a subset of gastric carcinomas. The genome-wide characteristics of EBV are essential to understand the diversity of strains isolated from EBV-related malignancies, provide the first opportunity to test the general validity of the EBV genetic map and explore recombination, geographic variation, and the major features of variation in this virus. Moreover, understanding more about EBV sequence variations isolated from EBV-related malignancies might give important implications for the development of effective prophylactic and therapeutic vaccine approaches targeting the personalized or geographic-specific EBV antigens in these aggressive diseases. In this chapter, we will mainly focus on the EBV genome-wide profiling in three common EBV-related cancers in Asia, including nasopharyngeal carcinoma, EBV-associated gastric carcinoma, and NK/T-cell lymphoma.

Genetic Patterns Found in the Nuclear Localization Signals (NLSs) Associated with EBV-1 and EBV-2 Provide New Insights into Their Contribution to Different Cell-Type Specificities

Simple Summary

The Epstein–Barr virus (EBV) has been implicated in several human neoplastic diseases. The EBV-1 can transform B cells into LCL more efficiently than EBV-2, and EBV-2 preferentially infects T-cell lymphocytes. The EBNA3A oncoprotein has an essential role in B-cell transformation. The six peptide motifs called nuclear localization signals (NLSs) from EBNA3A ensure nucleocytoplasmic protein trafficking. Multiple NLSs have been suggested to enhance EBNA3 function or different specificities to different cell types; however, a comprehensive assessment of their genetic variability has not been addressed. Our objective was to study the NLSs’ variability and their relationship with EBV types. Based on a comprehensive analysis of over a thousand EBNA3A sequences from different clinical manifestations and geographic locations, we found that EBNA3A from EBV-2 has two of the six NLSs altered, and genetic patterns in the NLSs are associated with EBV-1 and EBV-2.

Abstract

The Epstein–Barr virus (EBV) is a globally dispersed pathogen involved in several human cancers of B-cell and non-B-cell origin. EBV has been classified into EBV-1 and EBV-2, which have differences in their transformative ability. EBV-1 can transform B-cells into LCL more efficiently than EBV-2, and EBV-2 preferentially infects T-cell lymphocytes. The EBNA3A oncoprotein is a transcriptional regulator of virus and host cell genes, and is required in order to transform B-cells. EBNA3A has six peptide motifs called nuclear localization signals (NLSs) that ensure nucleocytoplasmic protein trafficking. The presence of multiple NLSs has been suggested to enhance EBNA3 function or different specificities in different cell types. However, studies about the NLS variability associated with EBV types are scarce. Based on a systematic sequence analysis considering more than a thousand EBNA3A sequences of EBV from different human clinical manifestations and geographic locations, we found differences in NLSs’ nucleotide structures among EBV types. Compared with the EBNA3A EBV-1, EBNA3A EBV-2 has two of the six NLSs altered, and these mutations were possibly acquired by recombination. These genetic patterns in the NLSs associated with EBV-1 and EBV-2 provide new information about the traits of EBNA3A in EBV biology.

A global phylogenetic analysis of Japanese tonsil-derived Epstein-Barr virus strains using viral whole-genome cloning and long-read sequencing

Epstein-Barr virus (EBV) establishes lifelong latent infection in the majority of healthy individuals, while it is a causative agent for various diseases, including some malignancies. Recent high-throughput sequencing results indicate that there are substantial levels of viral genome heterogeneity among different EBV strains. However, the extent of EBV strain variation among asymptomatically infected individuals remains elusive. Here, we present a streamlined experimental strategy to clone and sequence EBV genomes derived from human tonsillar tissues, which are the reservoirs of asymptomatic EBV infection. Complete EBV genome sequences, including those of repetitive regions, were determined for seven tonsil-derived EBV strains. Phylogenetic analyses based on the whole viral genome sequences of worldwide non-tumour-derived EBV strains revealed that Asian EBV strains could be divided into several distinct subgroups. EBV strains derived from nasopharyngeal carcinoma-endemic areas constitute different subgroups from a subgroup of EBV strains from non-endemic areas, including Japan. The results could be consistent with biased regional distribution of EBV-associated diseases depending on the different EBV strains colonizing different regions in Asian countries.

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.

Expanding the Geographic Characterisation of Epstein-Barr Virus Variation through Gene-Based Approaches

The Epstein-Barr Virus (EBV) infects the vast majority of human individuals worldwide (~90%) and is associated with several diseases, including different types of cancer and multiple sclerosis, which show wide variation in incidence among global geographical regions. Genetic variants in EBV genomic sequences have been used to determine the geographical structure of EBV isolates, but our understanding of EBV diversity remains highly incomplete. We generated sequences for 13 pivotal EBV genes derived from 103 healthy individuals, expanding current EBV diversity datasets with respect to both geographic coverage and number of isolates per region. These newly generated sequences were integrated with the more than 250 published EBV genomes, generating the most geographically comprehensive data set of EBV strains to date. We report remarkable variation in single-gene phylogenies that, when analysed together, show robust signals of population structure. Our results not only confirm known major global patterns of geographic variation, such as the clear separation of Asian isolates from the rest, and the intermixed relationships among African, European and Australian isolates, but yield novel phylogenetic relationships with previously unreported populations. We provide a better understanding of EBV's population structure in South America, Africa and, by the inclusion of Turkey and Georgia, we also gain insight into EBV diversity in Western Asia, a crossroads connecting Europe, Africa and Asia. In summary, our results provide a detailed world-wide characterisation of EBV genetic clusters, their enrichment in specific geographic regions, novel inter-population relationships, and a catalogue of geographically informative EBV genetic variants.

Formation of a DNA triple helical structure at BOLF1 gene of human herpesvirus 4 (HH4) genome

Eukaryotic genomes contain a large number of pyrimidine-purine rich regions and such regions can assume varied DNA conformations, including triple-stranded structures. These structures have fascinated scientists because of their considerable therapeutic applications. These structures have also profound implications in the field of nanotechnology as they can be used to develop DNA-based nanostructures and materials. Therefore, for any application, it is important to understand the formation of triplex structures, both in quantitative and qualitative terms. A combination of gel electrophoresis, UV-thermal denaturation and circular dichroism (CD) spectroscopy was used to investigate the formation of inter- as well as intramolecular triplex, in pyrimidine motif at BOLF1 gene of human herpesvirus 4 (HH4) genome. This gene codes for inner tegument protein, which plays crucial roles in viral replication. The said oligopurine•oligopyrimidine duplex was targeted via a designed triple helix forming oligopyrimidine nucleotide (TFO) in intermolecular as well as intramolecular fashion. Our studies revealed that intramolecular triplex formation takes place at acidic as well as at neutral pH; whereas low pH is required for its intermolecular version. This comparative study between inter- and intramolecular triplex allowed us to demonstrate that intramolecular structure is more stable to its intermolecular counterpart. Numerous models for mono-, bi- and trimolecular structures adopted by these DNA sequences have been suggested. This report adds to our existing knowledge about DNA triple helical structures.

Role of Viral and Host microRNAs in Immune Regulation of Epstein-Barr Virus-Associated Diseases

Epstein-Barr virus (EBV) is an oncogenic human herpes virus that was discovered in 1964. Viral non-coding RNAs, such as BamHI-A rightward fragment-derived microRNAs (BART miRNAs) or BamHI-H rightward fragment 1-derived miRNAs (BHRF1 miRNA) in EBV-infected cells have been recently reported. Host miRNAs are also upregulated upon EBV infection. Viral and host miRNAs are important in maintaining viral infection and evasion of host immunity. Although miRNAs in EBV-infected cells often promote cell proliferation by targeting apoptosis or cell cycle, this review focuses on the regulation of the recognition of the host immune system. This review firstly describes the location and organization of two clusters of viral miRNAs, then describes evasion from host immune surveillance systems by modulating viral gene expression or inhibiting innate and acquired immunity by viral miRNAs as well as host miRNAs. Another topic is the enigmatic depletion of viral miRNAs in several types of EBV-infected tumor cells. Finally, this review introduces the strong correlation of nasopharyngeal cancer cases with a newly identified single nucleotide polymorphism that enhances BART miRNA promoter activity.

Translational profiling of B cells infected with the Epstein-Barr virus reveals 5' leader ribosome recruitment through upstream open reading frames

The Epstein-Barr virus (EBV) genome encodes several hundred transcripts. We have used ribosome profiling to characterize viral translation in infected cells and map new translation initiation sites. We show here that EBV transcripts are translated with highly variable efficiency, owing to variable transcription and translation rates, variable ribosome recruitment to the leader region and coverage by monosomes versus polysomes. Some transcripts were hardly translated, others mainly carried monosomes, showed ribosome accumulation in leader regions and most likely represent non-coding RNAs. A similar process was visible for a subset of lytic genes including the key transactivators BZLF1 and BRLF1 in cells infected with weakly replicating EBV strains. This suggests that ribosome trapping, particularly in the leader region, represents a new checkpoint for the repression of lytic replication. We could identify 25 upstream open reading frames (uORFs) located upstream of coding transcripts that displayed 5′ leader ribosome trapping, six of which were located in the leader region shared by many latent transcripts. These uORFs repressed viral translation and are likely to play an important role in the regulation of EBV translation.

Co-evolution of sites under immune selection shapes Epstein-Barr Virus population structure

Epstein–Barr virus (EBV) is one of the most common viral infections in humans and persists within its host for life. EBV therefore represents an extremely successful virus that has evolved complex strategies to evade the host’s innate and adaptive immune response during both initial and persistent stages of infection. Here, we conducted a comparative genomics analysis on 223 whole genome sequences of worldwide EBV strains. We recover extensive genome-wide linkage disequilibrium (LD) despite pervasive genetic recombination. This pattern is explained by the global EBV population being subdivided into three main subpopulations, one primarily found in East Asia, one in Southeast Asia and Oceania, and the third including most of the other globally distributed genomes we analyzed. Additionally, sites in LD were overrepresented in immunogenic genes. Taken together, our results suggest that host immune selection and local adaptation to different human host populations has shaped the genome-wide patterns of genetic diversity in EBV.

Epstein-Barr virus strain variation and cancer

Epstein-Barr virus (EBV) is a human tumor virus and is etiologically linked to various malignancies. Certain EBV-associated diseases, such as Burkitt lymphomas and nasopharyngeal carcinomas, are endemic and exhibit biased geographic distribution worldwide. Recent advances in deep sequencing technology enabled high-throughput sequencing of the EBV genome from clinical samples. Rapid cloning and sequencing of cancer-derived EBV genomes, followed by reconstitution of infectious virus, have also become possible. These developments have revealed that various EBV strains are differentially distributed throughout the world, and that the behavior of cancer-derived EBV strains is different from that of the prototype EBV strain of non-cancerous origin. In this review, we summarize recent progress and future perspectives regarding the association between EBV strain variation and cancer.

Pathogenesis and biomarkers of natural killer T cell lymphoma (NKTL)

Natural killer T cell lymphoma (NKTL) is an aggressive disease with very poor treatment outcomes in the advanced stages. With chemotherapy, initial response rates to treatment are high but responses are short lived. A better understanding of the complex molecular pathogenesis of this disease is essential in order to design and develop better therapeutics with improved efficacy. This review aims to summarise the key pathogenic mechanisms in NKTL which may have significant prognostic and therapeutic implications.

The Phylogeographic Diversity of EBV and Admixed Ancestry in the Americas⁻Another Model of Disrupted Human-Pathogen Co-Evolution

Epstein-Barr virus (EBV) is an etiological agent for gastric cancer with significant worldwide variations. Molecular characterizations of EBV have shown phylogeographical variations among healthy populations and in EBV-associated diseases, particularly the cosegregated BamHI-I fragment and XhoI restriction site of exon 1 of the LMP-1 gene. In the Americas, both cosegregated variants are present in EBV carriers, which aligns with the history of Asian and European human migration to this continent. Furthermore, novel recombinant variants have been found, reflecting the genetic makeup of this continent. However, in the case of EBV-associated gastric cancer (EBV-associated GC), the cosegregated European BamHI-“i” fragment and XhoI restriction site strain prevails. Thus, we propose that a disrupted coevolution between viral phylogeographical strains and mixed human ancestry in the Americas might explain the high prevalence of this particular gastric cancer subtype. This cosegregated region contains two relevant transcripts for EBV-associated GC, the BARF-1 and miR-BARTs. Thus, genome-wide association studies (GWAS) or targeted sequencing of both transcripts may be required to clarify their role as a potential source of this disrupted coevolution.

The roles of EBV-encoded microRNAs in EBV-associated tumors

Epstein-Barr virus (EBV) is believed to be a pathogen causing a number of human cancers, but the pathogenic mechanisms remain unclear. An increasing number of studies have indicated that EBV-encoded microRNAs (EBV miRNAs) are expressed in a latency type- and tumor type-dependent manner, playing important roles in the development and progression of EBV-associated tumors. By targeting one or more genes of the virus and the host, EBV miRNAs are responsible for the deregulation of a variety of viral and host cell biological processes, including viral replication, latency maintenance, immune evasion, cell apoptosis and metabolism, and tumor proliferation and metastasis. In addition, some EBV miRNAs can be used as excellent diagnostic, prognostic and treatment efficacy predictive biomarkers for EBV-associated tumors. More importantly, EBV miRNA-targeting therapeutics have emerged and have been developing rapidly, which may open a new era in the treatment of EBV-associated tumors in the near future.

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.

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.

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.

Heterogeneity of the Epstein-Barr Virus (EBV) Major Internal Repeat Reveals Evolutionary Mechanisms of EBV and a Functional Defect in the Prototype EBV Strain B95-8

Epstein-Barr virus (EBV) is a ubiquitous pathogen of humans that can cause several types of lymphoma and carcinoma. Like other herpesviruses, EBV has diversified through both coevolution with its host and genetic exchange between virus strains. Sequence analysis of the EBV genome is unusually challenging because of the large number and lengths of repeat regions within the virus. Here we describe the sequence assembly and analysis of the large internal repeat 1 of EBV (IR1; also known as the BamW repeats) for more than 70 strains. The diversity of the latency protein EBV nuclear antigen leader protein (EBNA-LP) resides predominantly within the exons downstream of IR1. The integrity of the putative BWRF1 open reading frame (ORF) is retained in over 80% of strains, and deletions truncating IR1 always spare BWRF1. Conserved regions include the IR1 latency promoter (Wp) and one zone upstream of and two within BWRF1. IR1 is heterogeneous in 70% of strains, and this heterogeneity arises from sequence exchange between strains as well as from spontaneous mutation, with interstrain recombination being more common in tumor-derived viruses. This genetic exchange often incorporates regions of <1 kb, and allelic gene conversion changes the frequency of small regions within the repeat but not close to the flanks. These observations suggest that IR1-and, by extension, EBV-diversifies through both recombination and breakpoint repair, while concerted evolution of IR1 is driven by gene conversion of small regions. Finally, the prototype EBV strain B95-8 contains four nonconsensus variants within a single IR1 repeat unit, including a stop codon in the EBNA-LP gene. Repairing IR1 improves EBNA-LP levels and the quality of transformation by the B95-8 bacterial artificial chromosome (BAC).IMPORTANCE Epstein-Barr virus (EBV) infects the majority of the world population but causes illness in only a small minority of people. Nevertheless, over 1% of cancers worldwide are attributable to EBV. Recent sequencing projects investigating virus diversity to see if different strains have different disease impacts have excluded regions of repeating sequence, as they are more technically challenging. Here we analyze the sequence of the largest repeat in EBV (IR1). We first characterized the variations in protein sequences encoded across IR1. In studying variations within the repeat of each strain, we identified a mutation in the main laboratory strain of EBV that impairs virus function, and we suggest that tumor-associated viruses may be more likely to contain DNA mixed from two strains. The patterns of this mixing suggest that sequences can spread between strains (and also within the repeat) by copying sequence from another strain (or repeat unit) to repair DNA damage.

The conundrum of the Epstein-Barr virus-associated gastric carcinoma in the Americas

Epstein-Barr virus-associated gastric carcinoma shows a higher prevalence in the Americas than Asia. We summarize all studies of Epstein Barr virus-associated gastric carcinoma in the Americas, focusing on host characteristics, environmental associations and phylogeographic diversity of Epstein-Barr virus strains. In the Americas, the prevalence of Epstein Barr virus-associated gastric carcinoma is 11.4%, more frequent in males and portray predominantly diffuse-type histology. EBERs, EBNAs, BARTs and LMP are the highest expressed genes; their variations in healthy individuals may explain the phylogeographic diversity of Epstein-Barr virus across the region. Gastric cancer cases harbor exclusively the western genotype (subtype D and kept Xho I site), suggesting a disrupted co-evolution between the pathogen and its host. Epstein-Barr virus-associated gastric carcinoma molecular subtype cases from The Cancer Genome Atlas display PIK3CA gene mutations, amplification of JAK2, PD-L1 and PD-L2 and CpG island methylator phenotype, leading to more extensive methylation of host and viral genomes than any other subtypes from the study. Environmental conditions include negative- and positive- associations with being firstborn child and smoking, respectively. A marginal association with H. pylori has also been reported. Lymphoepithelioma-like carcinoma is associated with Epstein Barr virus in 80%-86% of cases, most of which have been included as part of Epstein Barr virus-associated gastric carcinoma series (prevalence 1.1%-7.6%). Whether these cases represent a variant of Epstein-Barr virus-associated gastric carcinoma is discussed. We propose novel research strategies to solve the conundrum of the high prevalence of Epstein-Barr virus-associated gastric carcinoma in the Americas.

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.

Epstein-Barr virus strains and variations: Geographic or disease-specific variants?

The Epstein-Barr Virus (EBV) is associated with the development of several diseases, including infectious mononucleosis (IM), Burkitt's Lymphoma (BL), Nasopharyngeal Carcinoma, and other neoplasias. The publication of EBV genome 1984 led to several studies regarding the identification of different viral strains. Currently, EBV is divided into EBV type 1 (B95-8 strain) and EBV type 2 (AG876 strain), also known as type A and type B, which have been distinguished based upon genetic differences in the Epstein-Barr nuclear antigens (EBNAs) sequence. Several other EBV strains have been described in the past 10 years considering variations on EBV genome, and many have attempted to clarify if these variations are ethnic or geographically correlated, or if they are disease related. Indeed, there is an increasing interest to describe possible specific disease associations, with emphasis on different malignancies. These studies aim to clarify if these variations are ethnic or geographically correlated, or if they are disease related, thus being important to characterize the epidemiologic genetic distribution of EBV strains on our population. Here, we review the current knowledge on the different EBV strains and variants and its association with different diseases. J. Med. Virol. 89:373-387, 2017. © 2016 Wiley Periodicals, Inc.

Epstein-Barr virus: more than 50 years old and still providing surprises

It is more than 50 years since the Epstein-Barr virus (EBV), the first human tumour virus, was discovered. EBV has subsequently been found to be associated with a diverse range of tumours of both lymphoid and epithelial origin. Progress in the molecular analysis of EBV has revealed fundamental mechanisms of more general relevance to the oncogenic process. This Timeline article highlights key milestones in the 50-year history of EBV and discusses how this virus provides a paradigm for exploiting insights at the molecular level in the diagnosis, treatment and prevention of cancer.

Genome-wide analysis of Epstein-Barr virus (EBV) isolated from EBV-associated gastric carcinoma (EBVaGC)

Epstein-Barr virus (EBV) is linked to the development of a variety of malignancies, including EBV-associated gastric carcinoma (EBVaGC). In this study, EBVaGC was detected in 15 (7.3%) of 206 GC cases. To identify the EBV genomic variation, EBV genomic sequences isolated from 9 EBVaGC biopsy specimens were successfully retrieved, designated EBVaGC1 to EBVaGC9. By comparative analysis of these strains with another 6 completely sequenced EBV strains, EBV-wild type, B95–8, AG876, GD1, GD2, and HKNPC1, it was demonstrated that EBVaGC1 to 9 were most closely related to the GD1 strain. Phylogenetic analysis of the GC biopsy specimen-derived EBV (GC-EBV) genomes was subsequently performed to assess their genomic diversity and it exhibited the greatest divergence from the type 2 strain, AG876. Compared with the reference EBV strain GD1, they harbored 961 variations in total, including 919 substitutions, 23 insertions, and 19 deletions. Single nucleotide polymorphism (SNP) density varied substantially across all known open reading frames and was highest in latency-associated genes. Moreover, we identified 2 interstrain recombinants at the EBNA1 locus, which provided a further mechanism for the generation of diversity. Some T-cell epitope sequences in EBNA1 and LMP2A genes showed extensive variation across strains, which implied their importance in the development of vaccines and T-cell therapy. In conclusion, we reported the first genome-wide view of sequence variation of EBV isolated from primary EBVaGC biopsy specimens, which might serve as an effective method for further understanding the genomic variations contribute to EBVaGC carcinogenesis and treatment.

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.

Highly Efficient CRISPR/Cas9-Mediated Cloning and Functional Characterization of Gastric Cancer-Derived Epstein-Barr Virus Strains

The Epstein-Barr virus (EBV) is etiologically linked to approximately 10% of gastric cancers, in which viral genomes are maintained as multicopy episomes. EBV-positive gastric cancer cells are incompetent for progeny virus production, making viral DNA cloning extremely difficult. Here we describe a highly efficient strategy for obtaining bacterial artificial chromosome (BAC) clones of EBV episomes by utilizing a CRISPR/Cas9-mediated strand break of the viral genome and subsequent homology-directed repair. EBV strains maintained in two gastric cancer cell lines (SNU719 and YCCEL1) were cloned, and their complete viral genome sequences were determined. Infectious viruses of gastric cancer cell-derived EBVs were reconstituted, and the viruses established stable latent infections in immortalized keratinocytes. While Ras oncoprotein overexpression caused massive vacuolar degeneration and cell death in control keratinocytes, EBV-infected keratinocytes survived in the presence of Ras expression. These results implicate EBV infection in predisposing epithelial cells to malignant transformation by inducing resistance to oncogene-induced cell death.

IMPORTANCE

Recent progress in DNA-sequencing technology has accelerated EBV whole-genome sequencing, and the repertoire of sequenced EBV genomes is increasing progressively. Accordingly, the presence of EBV variant strains that may be relevant to EBV-associated diseases has begun to attract interest. Clearly, the determination of additional disease-associated viral genome sequences will facilitate the identification of any disease-specific EBV variants. We found that CRISPR/Cas9-mediated cleavage of EBV episomal DNA enabled the cloning of disease-associated viral strains with unprecedented efficiency. As a proof of concept, two gastric cancer cell-derived EBV strains were cloned, and the infection of epithelial cells with reconstituted viruses provided important clues about the mechanism of EBV-mediated epithelial carcinogenesis. This experimental system should contribute to establishing the relationship between viral genome variation and EBV-associated diseases.

From Conventional to Next Generation Sequencing of Epstein-Barr Virus Genomes

Genomic sequences of Epstein–Barr virus (EBV) have been of interest because the virus is associated with cancers, such as nasopharyngeal carcinoma, and conditions such as infectious mononucleosis. The progress of whole-genome EBV sequencing has been limited by the inefficiency and cost of the first-generation sequencing technology. With the advancement of next-generation sequencing (NGS) and target enrichment strategies, increasing number of EBV genomes has been published. These genomes were sequenced using different approaches, either with or without EBV DNA enrichment. This review provides an overview of the EBV genomes published to date, and a description of the sequencing technology and bioinformatic analyses employed in generating these sequences. We further explored ways through which the quality of sequencing data can be improved, such as using DNA oligos for capture hybridization, and longer insert size and read length in the sequencing runs. These advances will enable large-scale genomic sequencing of EBV which will facilitate a better understanding of the genetic variations of EBV in different geographic regions and discovery of potentially pathogenic variants in specific diseases.

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.

Characterization of Epstein-Barr virus strains and LMP1-deletion variants in Portugal

Variations in the genome sequence of Epstein-Barr Virus (EBV) are thought to lead to differential interaction with host cells, immune evasion, and transformation. The discussion regarding EBV strains as having a geographic or disease-association has been increasing and the majority of studies are performed in Asiatic populations. We developed a case-control study with 139 individuals, including 96 subjects with different malignancies and 43 healthy individuals, from the North region of Portugal. We have used PCR protocols for the characterization of EBV strains (type A or B) based on EBNA3C genome variation and for the LMP1 30bp-deletion variants (wt-LMP1 or del-LMP1). Our study showed that type A is the most prevalent in our population (100% of healthy controls, 96.9% of aHSCT patients, 90.8% of HNSCC patients, and 94.9% of NPC patients) and that type B was significantly associated with NPC (P = 0.019; RR = 8.90). Regarding the LMP1 30bp-deletion, we found a similar distribution of both wt- and del-LMP1 variants in controls and dispare results in cases: del-LMP1 was more frequent in aHSCT and HNSCC patients (64.7% and 63.2%, respectively) and wt-LMP1 in NPC patients (100%). In fact, the study reveals that wt-LPM1 was significantly associated with NPC (P < 0.001; RR = 18.4). Hence, our study showed that EBV type B and wt-LMP1 variant seem to be associated with NPC in our population, with a clear disease-association for wt-LMP1. These results contribute for the knowledge of EBV genetic diversity among Caucasian populations.

The Epstein-Barr virus BamHI C promoter is not essential for B cell immortalization in vitro, but it greatly enhances B cell growth transformation

Epstein-Barr virus (EBV) infection of B cells leads to the sequential activation of two viral promoters, Wp and Cp, resulting in the expression of six EBV nuclear antigens (EBNAs) and the viral Bcl2 homologue BHRF1. The viral transactivator EBNA2 is required for this switch from Wp to Cp usage during the initial stages of infection. EBNA2-dependent Cp transcription is mediated by the EBNA2 response element (E2RE), a region that contains at least two binding sites for cellular factors; one of these sites, CBF1, interacts with RBP-JK, which then recruits EBNA2 to the transcription initiation complex. Here we demonstrate that the B cell-specific transcription factor BSAP/Pax5 binds to a second site, CBF2, in the E2RE. Deletion of the E2RE in the context of a recombinant virus greatly diminished levels of Cp-initiated transcripts during the initial stages of infection but did not affect the levels of Wp-initiated transcripts or EBNA mRNAs. Consistent with this finding, viruses deleted for the E2RE were not markedly impaired in their ability to induce B cell transformation in vitro. In contrast, a larger deletion of the entire Cp region did reduce EBNA mRNA levels early after infection and subsequently almost completely ablated lymphoblastoid cell line (LCL) outgrowth. Notably, however, rare LCLs could be established following infection with Cp-deleted viruses, and these were indistinguishable from wild-type-derived LCLs in terms of steady-state EBV gene transcription. These data indicate that, unlike Wp, Cp is dispensable for the virus' growth-transforming activity.

IMPORTANCE

Epstein-Barr virus (EBV), a B lymphotropic herpesvirus etiologically linked to several B cell malignancies, efficiently induces B cell proliferation leading to the outgrowth of lymphoblastoid cell lines (LCLs). The initial stages of this growth-transforming infection are characterized by the sequential activation of two viral promoters, Wp and Cp, both of which appear to be preferentially active in target B cells. In this work, we have investigated the importance of Cp activity in initiating B cell proliferation and maintaining LCL growth. Using recombinant viruses, we demonstrate that while Cp is not essential for LCL outgrowth in vitro, it enhances transformation efficiency by >100-fold. We also show that Cp, like Wp, interacts with the B cell-specific activator protein BSAP/Pax5. We suggest that EBV has evolved this two-promoter system to ensure efficient colonization of the host B cell system in vivo.

Epstein-Barr Virus: From the Detection of Sequence Polymorphisms to the Recognition of Viral Types

The Epstein-Barr virus is etiologically linked with the development of benign and malignant diseases, characterized by their diversity and a heterogeneous geographic distribution across the world. The virus possesses a 170-kb-large genome that encodes for multiple proteins and non-coding RNAs. Early on there have been numerous attempts to link particular diseases with particular EBV strains, or at least with viral genetic polymorphisms. This has given rise to a wealth of information whose value has been difficult to evaluate for at least four reasons. First, most studies have looked only at one particular gene and missed the global picture. Second, they usually have not studied sufficient numbers of diseased and control cases to reach robust statistical significance. Third, the functional significance of most polymorphisms has remained unclear, although there are exceptions such as the 30-bp deletion in LMP1. Fourth, different biological properties of the virus do not necessarily equate with a different pathogenicity. This was best illustrated by the type 1 and type 2 viruses that markedly differ in terms of their transformation abilities, yet do not seem to cause different diseases. Reciprocally, environmental and genetic factors in the host are likely to influence the outcome of infections with the same virus type. However, with recent developments in recombinant virus technology and in the availability of high throughput sequencing, the tide is now turning. The availability of 23 complete or nearly complete genomes has led to the recognition of viral subtypes, some of which possess nearly identical genotypes. Furthermore, there is growing evidence that some genetic polymorphisms among EBV strains markedly influence the biological and clinical behavior of the virus. Some virus strains are endowed with biological properties that explain crucial clinical features of patients with EBV-associated diseases. Although we now have a better overview of the genetic diversity within EBV genomes, it has also become clear that defining phenotypic traits evinced by cells infected by different viruses usually result from the combination of multiple polymorphisms that will be difficult to identify in their entirety. However, the steadily increasing number of sequenced EBV genomes and cloned EBV BACS from diseased and healthy patients will facilitate the identification of the key polymorphisms that condition the biological and clinical behavior of the viruses. This will allow the development of preventative and therapeutic approaches against highly pathogenic viral strains.

Epstein-Barr Virus Strain Variation

What is wild-type Epstein-Barr virus and are there genetic differences in EBV strains that contribute to some of the EBV-associated diseases? Recent progress in DNA sequencing has resulted in many new Epstein-Barr virus (EBV) genome sequences becoming available. EBV isolates worldwide can be grouped into type 1 and type 2, a classification based on the EBNA2 gene sequence. Type 1 transforms human B cells into lymphoblastoid cell lines much more efficiently than type 2 EBV and molecular mechanisms that may account for this difference in cell transformation are now becoming understood. Study of geographic variation of EBV strains independent of the type 1/type 2 classification and systematic investigation of the relationship between viral strains, infection and disease are now becoming possible. So we should consider more directly whether viral sequence variation might play a role in the incidence of some EBV-associated diseases.

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.

Clustered microRNAs of the Epstein-Barr virus cooperatively downregulate an epithelial cell-specific metastasis suppressor

The Epstein-Barr virus (EBV) encodes its own microRNAs (miRNAs); however, their biological roles remain elusive. The commonly used EBV B95-8 strain lacks a 12-kb genomic region, known as BamHI A rightward transcripts (BART) locus, where a number of BART miRNAs are encoded. Here, bacterial artificial chromosome (BAC) technology was used to generate an EBV B95-8 strain in which the 12-kb region was fully restored at its native locus [BART(+) virus]. Epithelial cells were stably infected with either the parental B95-8 virus or the BART(+) virus, and BART miRNA expression was successfully reconstituted in the BART(+) virus-infected cells. Microarray analyses of cellular gene expression identified N-myc downstream regulated gene 1 (NDRG1) as a putative target of BART miRNAs. The NDRG1 protein was barely expressed in B cells, highly expressed in epithelial cells, including primary epithelial cells, and strongly downregulated in the BART(+) virus-infected epithelial cells of various origins. Although in vitro reporter assays identified BART22 as being responsible for the NDRG1 downregulation, EBV genetic analyses revealed that BART22 was not solely responsible; rather, the entire BART miRNA cluster 2 was responsible for the downregulation. Immunohistochemical analyses revealed that the expression level of the NDRG1 protein was downregulated significantly in EBV-positive nasopharyngeal carcinoma specimens. Considering that NDRG1 encodes an epithelial differentiation marker and a suppressor of metastasis, these data implicate a causative relationship between BART miRNA expression and epithelial carcinogenesis in vivo.

IMPORTANCE

EBV-related epithelial cancers, such as nasopharyngeal carcinomas and EBV-positive gastric cancers, encompass more than 80% of EBV-related malignancies. Although it is known that they express high levels of virally encoded BART miRNAs, how these miRNAs contribute to EBV-mediated epithelial carcinogenesis remains unknown. Although a number of screenings have been performed to identify targets of viral miRNAs, many targets likely have not been identified, especially in case of epithelial cell infection. This is the first study to use EBV genetics to perform unbiased screens of cellular genes that are differentially expressed in viral miRNA-positive and -negative epithelial cells. The result indicates that multiple EBV-encoded miRNAs cooperatively downregulate NDRG1, an epithelial differentiation marker and suppressor of metastasis. The experimental system described in this study should be useful for further clarifying the mechanism of EBV-mediated epithelial carcinogenesis.

Genomic diversity of Epstein-Barr virus genomes isolated from primary nasopharyngeal carcinoma biopsy samples

Undifferentiated nasopharyngeal carcinoma (NPC) has a 100% association with Epstein-Barr virus (EBV). However, only three EBV genomes isolated from NPC patients have been sequenced to date, and the role of EBV genomic variations in the pathogenesis of NPC is unclear. We sought to obtain the sequences of EBV genomes in multiple NPC biopsy specimens in the same geographic location in order to reveal their sequence diversity. Three published EBV (B95-8, C666-1, and HKNPC1) genomes were first resequenced using the sequencing workflow of target enrichment of EBV DNA by hybridization, followed by next-generation sequencing, de novo assembly, and joining of contigs by Sanger sequencing. The sequences of eight NPC biopsy specimen-derived EBV (NPC-EBV) genomes, designated HKNPC2 to HKNPC9, were then determined. They harbored 1,736 variations in total, including 1,601 substitutions, 64 insertions, and 71 deletions, compared to the reference EBV. Furthermore, genes encoding latent, early lytic, and tegument proteins and glycoproteins were found to contain nonsynonymous mutations of potential biological significance. Phylogenetic analysis showed that the HKNPC6 and -7 genomes, which were isolated from tumor biopsy specimens of advanced metastatic NPC cases, were distinct from the other six NPC-EBV genomes, suggesting the presence of at least two parental lineages of EBV among the NPC-EBV genomes. In conclusion, much greater sequence diversity among EBV isolates derived from NPC biopsy specimens is demonstrated on a whole-genome level through a complete sequencing workflow. Large-scale sequencing and comparison of EBV genomes isolated from NPC and normal subjects should be performed to assess whether EBV genomic variations contribute to NPC pathogenesis.

IMPORTANCE

This study established a sequencing workflow from EBV DNA capture and sequencing to de novo assembly and contig joining. We reported eight newly sequenced EBV genomes isolated from primary NPC biopsy specimens and revealed the sequence diversity on a whole-genome level among these EBV isolates. At least two lineages of EBV strains are observed, and recombination among these lineages is inferred. Our study has demonstrated the value of, and provided a platform for, genome sequencing of EBV.

Fast selection of miRNA candidates based on large-scale pre-computed MFE sets of randomized sequences

Background

Small RNAs are important regulators of genome function, yet their prediction in genomes is still a major computational challenge. Statistical analyses of pre-miRNA sequences indicated that their 2D structure tends to have a minimal free energy (MFE) significantly lower than MFE values of equivalently randomized sequences with the same nucleotide composition, in contrast to other classes of non-coding RNA. The computation of many MFEs is, however, too intensive to allow for genome-wide screenings.

Results

Using a local grid infrastructure, MFE distributions of random sequences were pre-calculated on a large scale. These distributions follow a normal distribution and can be used to determine the MFE distribution for any given sequence composition by interpolation. It allows on-the-fly calculation of the normal distribution for any candidate sequence composition.

Conclusion

The speedup achieved makes genome-wide screening with this characteristic of a pre-miRNA sequence practical. Although this particular property alone will not be able to distinguish miRNAs from other sequences sufficiently discriminative, the MFE-based P-value should be added to the parameters of choice to be included in the selection of potential miRNA candidates for experimental verification.

Identification and characterization of EBV genomes in spontaneously immortalized human peripheral blood B lymphocytes by NGS technology

Background

We conducted genomic sequencing to identify Epstein Barr Virus (EBV) genomes in 2 human peripheral blood B lymphocytes that underwent spontaneous immortalization promoted by mycoplasma infections in culture, using the high-throughput sequencing (HTS) Illumina MiSeq platform. The purpose of this study was to examine if rapid detection and characterization of a viral agent could be effectively achieved by HTS using a platform that has become readily available in general biology laboratories.

Results

Raw read sequences, averaging 175 bps in length, were mapped with DNA databases of human, bacteria, fungi and virus genomes using the CLC Genomics Workbench bioinformatics tool. Overall 37,757 out of 49,520,834 total reads in one lymphocyte line (# K4413-Mi) and 28,178 out of 45,335,960 reads in the other lymphocyte line (# K4123-Mi) were identified as EBV sequences. The two EBV genomes with estimated 35.22-fold and 31.06-fold sequence coverage respectively, designated K4413-Mi EBV and K4123-Mi EBV (GenBank accession number KC440852 and KC440851 respectively), are characteristic of type-1 EBV.

Conclusions

Sequence comparison and phylogenetic analysis among K4413-Mi EBV, K4123-Mi EBV and the EBV genomes previously reported to GenBank as well as the NA12878 EBV genome assembled from database of the 1000 Genome Project showed that these 2 EBVs are most closely related to B95-8, an EBV previously isolated from a patient with infectious mononucleosis and WT-EBV. They are less similar to EBVs associated with nasopharyngeal carcinoma (NPC) from Hong Kong and China as well as the Akata strain of a case of Burkitt’s lymphoma from Japan. They are most different from type 2 EBV found in Western African Burkitt’s lymphoma.

CD40L expression permits CD8+ T cells to execute immunologic helper functions

CD8(+) T cells play an essential role in immunity against intracellular pathogens, with cytotoxicity being considered their major effector mechanism. However, we here demonstrate that a major part of central and effector memory CD8(+) T cells expresses CD40L, one key molecule for CD4(+) T-cell-mediated help. CD40L(+) CD8(+) T cells are detectable among human antigen-specific immune responses, including pathogens such as influenza and yellow fever virus. CD40L(+) CD8(+) T cells display potent helper functions in vitro and in vivo, such as activation of antigen-presenting cells, and exhibit a cytokine expression signature similar to CD4(+) T cells and unrelated to cytotoxic CD8(+) T cells. The broad occurrence of CD40L(+) CD8(+) T cells in cellular immunity implicates that helper functions are not only executed by major histocompatibility complex (MHC) class II-restricted CD4(+) helper T cells but are also a common feature of MHC class I-restricted CD8(+) T cell responses. Due to their versatile functional capacities, human CD40L(+) CD8(+) T cells are promising candidate cells for immune therapies, particularly when CD4(+) T-cell help or pathogen-associated molecular pattern signals are limited.

Epstein-barr virus sequence variation-biology and disease

Some key questions in Epstein-Barr virus (EBV) biology center on whether naturally occurring sequence differences in the virus affect infection or EBV associated diseases. Understanding the pattern of EBV sequence variation is also important for possible development of EBV vaccines. At present EBV isolates worldwide can be grouped into Type 1 and Type 2, a classification based on the EBNA2 gene sequence. Type 1 EBV is the most prevalent worldwide but Type 2 is common in parts of Africa. Type 1 transforms human B cells into lymphoblastoid cell lines much more efficiently than Type 2 EBV. Molecular mechanisms that may account for this difference in cell transformation are now becoming clearer. Advances in sequencing technology will greatly increase the amount of whole EBV genome data for EBV isolated from different parts of the world. Study of regional variation of EBV strains independent of the Type 1/Type 2 classification and systematic investigation of the relationship between viral strains, infection and disease will become possible. The recent discovery that specific mutation of the EBV EBNA3B gene may be linked to development of diffuse large B cell lymphoma illustrates the importance that mutations in the virus genome may have in infection and human disease.

The pre-NH(2)-terminal domain of the herpes simplex virus 1 DNA polymerase catalytic subunit is required for efficient viral replication

The catalytic subunit of herpes simplex virus 1 DNA polymerase (HSV-1 Pol) has been extensively studied; however, its full complement of functional domains has yet to be characterized. A crystal structure has revealed a previously uncharacterized pre-NH(2)-terminal domain (residues 1 to 140) within HSV-1 Pol. Due to the conservation of the pre-NH(2)-terminal domain within the herpesvirus Pol family and its location in the crystal structure, we hypothesized that this domain provides an important function during viral replication in the infected cell distinct from 5'-3' polymerase activity. We identified three pre-NH(2)-terminal Pol mutants that exhibited 5'-3' polymerase activity indistinguishable from that of wild-type Pol in vitro: deletion mutants PolΔN43 and PolΔN52 that lack the extreme N-terminal 42 and 51 residues, respectively, and mutant PolA(6), in which a conserved motif at residues 44 to 49 was replaced with alanines. We constructed the corresponding pol mutant viruses and found that the polΔN43 mutant displayed replication kinetics similar to those of wild-type virus, while polΔN52 and polA(6) mutant virus infection resulted in an 8-fold defect in viral yield compared to that achieved with wild type and their respective rescued derivative viruses. Additionally, both polΔN52 and polA(6) viruses exhibited defects in viral DNA synthesis that correlated with the observed reduction in viral yield. These results strongly indicate that the conserved motif within the pre-NH(2)-terminal domain is important for viral DNA synthesis and production of infectious virus and indicate a functional role for this domain.

DNA of a circular minichromosome linearized by restriction enzymes or other reagents is resistant to further cleavage: an influence of chromatin topology on the accessibility of DNA

The accessibility of DNA in chromatin is an essential factor in regulating its activities. We studied the accessibility of the DNA in a ∼170 kb circular minichromosome to DNA-cleaving reagents using pulsed-field gel electrophoresis and fibre-fluorescence in situ hybridization on combed DNA molecules. Only one of several potential sites in the minichromosome DNA was accessible to restriction enzymes in permeabilized cells, and in growing cells only a single site at an essentially random position was cut by poisoned topoisomerase II, neocarzinostatin and γ-radiation, which have multiple potential cleavage sites; further sites were then inaccessible in the linearized minichromosomes. Sequential exposure to combinations of these reagents also resulted in cleavage at only a single site. Minichromosome DNA containing single-strand breaks created by a nicking endonuclease to relax any unconstrained superhelicity was also cut at only a single position by a restriction enzyme. Further sites became accessible after ≥95% of histones H2A, H2B and H1, and most non-histone proteins were extracted. These observations suggest that a global rearrangement of the three-dimensional packing and interactions of nucleosomes occurs when a circular minichromosome is linearized and results in its DNA becoming inaccessible to probes.

Diversity of Epstein-Barr virus BamHI-A rightward transcripts and their expression patterns in lytic and latent infections

Epstein-Barr virus (EBV) BamHI-A rightward transcripts (BARTs; also designated complementary strand transcripts or CSTs) have been demonstrated to contain several splicing forms in EBV-infected cells. To date, however, little is known about the actual full-length splicing form and its functions. In the present study, we proved that six forms of BARTs were present in EBV-positive cell lines and various tissue specimens with different EBV infection patterns. Of the BART-encoded genes, mRNA of four major splicing forms, including BARF0, RPMS1, RPMS1A and A73, were expressed in all EBV-infected cells. On the other hand, mRNA of two minor splicing forms, RK-BARF0 and RB3, was rarely detected, or if at all, at very low expression levels. Both RPMS1A and RPMS1 mRNA was transcribed at higher levels in EBV-infected cells. In particular, RPMS1 mRNA was expressed abundantly in epithelial carcinoma cells, including gastric carcinoma and nasopharyngeal carcinoma, in association with a lytic infection signal, BZLF1 mRNA. The four major splicing forms were expressed much less in B-cell lines with an integrated EBV genome than in those with episomal EBV genomes. These data indicate that at least six splicing forms can be expressed by EBV-infected cells or tissues, although the expression patterns or levels differ for different infection states such as lytic and latent infections.