Natural Variation of Epstein-Barr Virus Genes, Proteins, and Primary MicroRNA

Lytic and Latent Genetic Diversity of the Epstein-Barr Virus Reveals Raji-Related Variants from Southeastern Brazil Associated with Recombination Markers

Epstein-Barr virus (EBV) is a ubiquitous gammaherpesvirus etiologically associated with benign and malignant diseases. Since the pathogenic mechanisms of EBV are not fully understood, understanding EBV genetic diversity is an ongoing goal. Therefore, the present work describes the genetic diversity of the lytic gene BZLF1 in a sampling of 70 EBV-positive cases from southeastern Brazil. Additionally, together with the genetic regions previously characterized, the aim of the present study was to determine the impact of viral genetic factors that may influence EBV genetic diversity. Accordingly, the phylogenetic analysis of the BZLF1 indicated two main clades with high support, BZ-A and BZ-B (PP > 0.85). Thus, the BZ-A clade was the most diverse clade associated with the main polymorphisms investigated, including the haplotype Type 1 + V3 (p < 0.001). Furthermore, the multigene phylogenetic analysis (MLA) between BZLF1 and the oncogene LMP1 showed specific clusters, revealing haplotypic segregation that previous single-gene phylogenies from both genes failed to demonstrate. Surprisingly, the LMP1 Raji-related variant clusters were shown to be more diverse, associated with BZ-A/B and the Type 2/1 + V3 haplotypes. Finally, due to the high haplotypic diversity of the Raji-related variants, the number of DNA recombination-inducing motifs (DRIMs) was evaluated within the different clusters defined by the MLA. Similarly, the haplotype BZ-A + Raji was shown to harbor a greater number of DRIMs (p < 0.001). These results call attention to the high haplotype diversity of EBV in southeast Brazil and strengthen the hypothesis of the recombinant potential of South American Raji-related variants via the LMP1 oncogene.

Meta-analysis of Epstein-Barr virus genomes in Southern Chinese identifies genetic variants and high risk viral lineage associated with nasopharyngeal carcinoma

Genetic variants in Epstein-Barr virus (EBV) have been strongly associated with nasopharyngeal carcinoma (NPC) in South China. However, different results regarding the most significant viral variants, with polymorphisms in EBER2 and BALF2 loci, have been reported in separate studies. In this study, we newly sequenced 100 EBV genomes derived from 61 NPC cases and 39 population controls. Comprehensive genomic analyses of EBV sequences from both NPC patients and healthy carriers in South China were conducted, totaling 279 cases and 227 controls. Meta-analysis of genome-wide association study revealed a 4-bp deletion downstream of EBER2 (coordinates, 7188-7191; EBER-del) as the most significant variant associated with NPC. Furthermore, multiple viral variants were found to be genetically linked to EBER-del forming a risk haplotype, suggesting that multiple viral variants might be associated with NPC pathogenesis. Population structure and phylogenetic analyses further characterized a high risk EBV lineage for NPC revealing a panel of 38 single nucleotide polymorphisms (SNPs), including those in the EBER2 and BALF2 loci. With linkage disequilibrium clumping and feature selection algorithm, the 38 SNPs could be narrowed down to 9 SNPs which can be used to accurately detect the high risk EBV lineage. In summary, our study provides novel insight into the role of EBV genetic variation in NPC pathogenesis by defining a risk haplotype of EBV for downstream functional studies and identifying a single high risk EBV lineage characterized by 9 SNPs for potential application in population screening of NPC.

Characterization of target gene regulation by the two Epstein-Barr virus oncogene LMP1 domains essential for B-cell transformation

IMPORTANCE

Epstein-Barr virus (EBV) causes multiple human cancers, including B-cell lymphomas. In cell culture, EBV converts healthy human B-cells into immortalized ones that grow continuously, which model post-transplant lymphomas. Constitutive signaling from two cytoplasmic tail domains of the EBV oncogene latent membrane protein 1 (LMP1) is required for this transformation, yet there has not been systematic analysis of their host gene targets. We identified that only signaling from the membrane proximal domain is required for survival of these EBV-immortalized cells and that its loss triggers apoptosis. We identified key LMP1 target genes, whose abundance changed significantly with loss of LMP1 signals, or that were instead upregulated in response to switching on signaling by one or both LMP1 domains in an EBV-uninfected human B-cell model. These included major anti-apoptotic factors necessary for EBV-infected B-cell survival. Bioinformatics analyses identified clusters of B-cell genes that respond differently to signaling by either or both domains.

Unique synapomorphies and high diversity in South American Raji-related Epstein-Barr virus genomes

BACKGROUND

Epstein-Barr virus (EBV) is a human gammaherpesvirus etiologically linked to several benign and malignant diseases. EBV-associated malignancies exhibit an unusual global distribution that might be partly attributed to virus and host genetic backgrounds.

OBJECTIVES

To assemble a new genome of EBV (CEMO3) from a paediatric Burkitt's lymphoma from Rio de Janeiro State (Southeast Brazil). In addition, to perform global phylogenetic analysis using complete EBV genomes, including CEMO3, and investigate the genetic relationship of some South American (SA) genomes through EBV subgenomic targets.

METHODS

CEMO3 was sequenced through next generation sequencing and its coverage and gaps were corrected through the Sanger method. CEMO3 and 67 EBV genomes representing diverse geographic regions were evaluated through maximum likelihood phylogenetic analysis. Further, the polymorphism of subgenomic regions of some SA EBV genomes were assessed.

FINDINGS

The whole bulk tumour sequencing yielded 23,217 reads related to EBV, which 172,713 base pairs of the newly EBV genome CEMO3 was assembled. The CEMO3 and most SA EBV genomes clustered within the SA subclade closely related to the African Raji strain, forming the South American/Raji clade. Notably, these Raji-related genomes exhibit significant genetic diversity, characterised by distinctive synapomorphies at some gene levels absent in the original Raji strain.

CONCLUSION

The CEMO3 represents a new South American EBV genome assembled. Albeit the majority of EBV genomes from SA are Raji-related, it harbours a high diversity different from the original Raji strain.

Exploring the Genetic Diversity of Epstein-Barr Virus among Patients with Gastric Cancer in Southern Chile

The Epstein-Barr virus (EBV) has been associated with gastric cancer (GC), one of the deadliest malignancies in Chile and the world. Little is known about Chilean EBV strains. This study aims to investigate the frequency and genetic diversity of EBV in GC in patients in southern Chile. To evaluate the prevalence of EBV in GC patients from the Chilean population, we studied 54 GC samples using the gold standard detection method of EBV-encoded small RNA (EBER). The EBV-positive samples were subjected to amplification and sequencing of the Epstein-Barr virus nuclear protein 3A (EBNA3A) gene to evaluate the genetic diversity of EBV strains circulating in southern Chile. In total, 22.2% of the GC samples were EBV-positive and significantly associated with diffuse-type histology (p = 0.003). Phylogenetic analyses identified EBV-1 and EBV-2 in the GC samples, showing genetic diversity among Chilean isolates. This work provides important information for an epidemiological follow-up of the different EBV subtypes that may cause GC in southern Chile.

Distribution of Epstein-Barr Virus LMP1 Variants in Patients with Infectious Mononucleosis and Association with Selected Biochemical and Hematological Parameters

The molecular diversity of Epstein-Barr virus (EBV) is exceptionally complex and based on the characterization of sequences coding for several viral genes. The aim of this study was to analyze the distribution of EBV types 1 and 2 and to characterize LMP1 variants in a cohort of 73 patients with infectious mononucleosis (IM), as well as to investigate a possible association between viral diversity and relevant clinical parameters. Population-based sequencing of EBNA-2 gene showed the presence of EBV type 1 in all IM patients. Analysis of LMP1 gene found a restricted repertoire of LMP1 variants with the predominance of wild-type B95-8, China1, Mediterranean and North Carolina variants with the presence of more than one LMP1 variant in 16.4% of patients. Co-infections with different LMP1 variants were associated with significantly higher levels of C-reactive protein and lower levels of maximal neutrophil counts and minimal platelet count. The results of this study have shown a narrow repertoire of LMP1 variants and an exclusive presence of EBV type 1 in a cohort of IM from Croatia, suggesting a characteristic local molecular pattern of this virus. The clinical importance of distinct immunobiological features of IM patients with LMP1 variant co-infections needs to be investigated further.

T cell-mediated immunity during Epstein-Barr virus infections in children

Epstein-Barr virus (EBV) infection is extremely common worldwide, with approximately 90% of adults testing positive for EBV antibodies. Human are susceptible to EBV infection, and primary EBV infection typically occurs early in life. EBV infection can cause infectious mononucleosis (IM) as well as some severe non-neoplastic diseases, such as chronic active EBV infection (CAEBV) and EBV-associated hemophagocytic lymphohistiocytosis (EBV-HLH), which can have a heavy disease burden. After primary EBV infection, individuals develop robust EBV-specific T cell immune responses, with EBV-specific CD8⁺ and part of CD4⁺ T cells functioning as cytotoxic T cells, defending against virus. Different proteins expressed during EBV's lytic replication and latent proliferation can cause varying degrees of cellular immune responses. Strong T cell immunity plays a key role in controlling infection by decreasing viral load and eliminating infected cells. However, the virus persists as latent infection in EBV healthy carriers even with robust T cell immune response. When reactivated, it undergoes lytic replication and then transmits virions to a new host. Currently, the relationship between the pathogenesis of lymphoproliferative diseases and the adaptive immune system is still not fully clarified and needs to be explored in the future. Investigating the T cell immune responses evoked by EBV and utilizing this knowledge to design promising prophylactic vaccines are urgent issues for future research due to the importance of T cell immunity.

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.

Molecular Characterisation of Epstein-Barr Virus in Classical Hodgkin Lymphoma

Hodgkin lymphomas (HLs) are a heterogeneous group of lymphoid neoplasia associated with Epstein-Barr virus (EBV) infection. EBV, considered to be an important etiological co-factor in approximately 1% of human malignancies, can be classified into two genotypes based on EBNA-2, EBNA-3A and EBNA-3C sequences, and into genetic variants based on the sequence variation of the gene coding for the LMP1 protein. Here, we present the results on the distribution of EBV genotypes 1 and 2 as well as LMP1 gene variants in 50 patients with EBV-positive classical HL selected from a cohort of 289 histologically verified cases collected over a 9-year period in a tertiary clinical center in the Southeast of Europe. The population-based sequencing of the EBNA-3C gene showed the exclusive presence of EBV genotype 1 in all cHL samples. The analysis of EBV LMP1 variant distribution showed a predominance of the wild-type strain B95-8 and the Mediterranean subtype with 30 bp deletion. These findings could contribute to the understanding of EBV immunobiology in cHL as well as to the development of a prophylactic and therapeutic vaccine.

Type 1 and Type 2 Epstein-Barr viruses induce proliferation, and inhibit differentiation, in infected telomerase-immortalized normal oral keratinocytes

Differentiated epithelial cells are an important source of infectious EBV virions in human saliva, and latent Epstein-Barr virus (EBV) infection is strongly associated with the epithelial cell tumor, nasopharyngeal carcinoma (NPC). However, it has been difficult to model how EBV contributes to NPC, since EBV has not been shown to enhance proliferation of epithelial cells in monolayer culture in vitro and is not stably maintained in epithelial cells without antibiotic selection. In addition, although there are two major types of EBV (type 1 (T1) and type 2 (T2)), it is currently unknown whether T1 and T2 EBV behave differently in epithelial cells. Here we inserted a G418 resistance gene into the T2 EBV strain, AG876, allowing us to compare the phenotypes of T1 Akata virus versus T2 AG876 virus in a telomerase-immortalized normal oral keratinocyte cell line (NOKs) using a variety of different methods, including RNA-seq analysis, proliferation assays, immunoblot analyses, and air-liquid interface culture. We show that both T1 Akata virus infection and T2 AG876 virus infection of NOKs induce cellular proliferation, and inhibit spontaneous differentiation, in comparison to the uninfected cells when cells are grown without supplemental growth factors in monolayer culture. T1 EBV and T2 EBV also have a similar ability to induce epithelial-to-mesenchymal (EMT) transition and activate canonical and non-canonical NF-κB signaling in infected NOKs. In contrast to our recent results in EBV-infected lymphoblastoid cells (in which T2 EBV infection is much more lytic than T1 EBV infection), we find that NOKs infected with T1 and T2 EBV respond similarly to lytic inducing agents such as TPA treatment or differentiation. These results suggest that T1 and T2 EBV have similar phenotypes in infected epithelial cells, with both EBV types enhancing cellular proliferation and inhibiting differentiation when growth factors are limiting.

Spatial Dispersal of Epstein-Barr Virus in South America Reveals an African American Variant in Brazilian Lymphomas

Epstein−Barr virus (EBV) is a saliva-borne ɣ-herpesvirus associated with benign and malignant lymphoproliferation. EBV-mediated tumorigenic mechanisms are not fully understood and may be related to viral genetic variations. In this work, we characterize the genetic diversity of EBV from Brazil, assessing 82 samples derived from saliva from asymptomatic carriers (n = 45), biopsies of benign reactive hyperplasia (n = 4), and lymphomas (n = 33). Phylogenetic and phylogeographic analysis of the entire coding region of the LMP-1 was performed. Additionally, type 1/type 2 distinction by the EBNA3C gene and Zp variants were evaluated. Our results revealed a high diversity of EBV in Brazil, with the co-circulation of four main clades, described here as: Mediterranean (40.2%, n = 33), Raji/Argentine (39%, n = 32), B95-8 (6.1%, n = 5), and Asian II (1.2%, n = 1). The Raji/Argentine and Mediterranean clades were the most prevalent in South America (45% and 28%, respectively). The Raji/Argentine clade was associated with polymorphisms I124V/I152L, del30 bp, and ins15 bp (p < 0.0001, to all clades) and with a high haplotype diversity related to EBV type and Zp variants. We found that a Raji/Argentine subclade spread primarily from Brazil and later to other South American countries. Although no LMP1 variant has been directly associated with disease, the Raji/Argentine clade was predominantly clustered with lymphomas (61%) and the Mediterranean clade with non-malignant cases (59%) (p = 0.1). These data highlight the high genetic diversity of EBV circulating in Brazil, calling attention to a Raji-related variant with great recombination potential in Brazilian lymphomas.

Epstein-Barr Virus Genome Deletions in Epstein-Barr Virus-Positive T/NK Cell Lymphoproliferative Diseases

The main target cells for Epstein-Barr virus (EBV) infection and persistence are B lymphocytes, although T and NK cells can also become infected. In this paper, we characterize the EBV present in 21 pediatric and adult patients who were treated in France for a range of diseases that involve infection of T or NK cells. Of these 21 cases, 5 pediatric patients (21%) and 11 adult patients (52%) were of Caucasian origin. In about 30% of the cases, some of the EBV genomes contain a large deletion. The deletions are different in every patient but tend to cluster near the BART region of the viral genome. Detailed investigation of a family in which several members have persistent T or NK cell infection by EBV indicates that the virus genome deletions arise or are selected independently in each individual patient. Genome sequence polymorphisms in the EBV in these T or NK cell diseases reflect the geographic origin of the patient and not a distinct type of EBV (the 21 cases studied included examples of both type 1 and type 2 EBV infection). Using virus produced from type 1 or type 2 EBV genomes cloned in bacterial artificial chromosome (BAC) vectors, we demonstrate infection of T cells in cord blood from healthy donors. Our results are consistent with transient infection of some T cells being part of normal asymptomatic infection by EBV in young children. IMPORTANCE EBV contributes to several types of human cancer. Some cancers and nonmalignant lymphoproliferative diseases involving T or NK cells contain EBV. These diseases are relatively frequent in Japan and China and have been shown sometimes to have deletions in the EBV genome in the disease cells. We identify further examples of deletions within the EBV genome associated with T or NK cell diseases, and we provide evidence that the virus genomes with these deletions are most likely selected in the individual cases, rather than being transmitted between people during infection. We demonstrate EBV infection of cord blood T cells by highly characterized, cloned EBV genomes and suggest that transient infection of T cells may be part of normal asymptomatic infection by EBV in young children.

Epstein-Barr Virus (EBV) Is Mostly Latent and Clonal in Angioimmunoblastic T Cell Lymphoma (AITL)

The Epstein-Barr virus (EBV) is associated with angioimmunoblastic T cell lymphoma (AITL), a peripheral T lymphoma of poor prognosis in at least 90% of cases. The role of EBV in this pathology is unknown. Using next-generation sequencing, we sequenced the entire EBV genome in biopsies from 18 patients with AITL, 16 patients with another EBV-associated lymphoma, and 2 controls. We chose an EBV target capture method, given the high specificity of this technique, followed by a second capture to increase sensitivity. We identified two main viral strains in AITL, one of them associated with the mutations BNRF1 S542N and BZLF1 A206S and with mutations in the EBNA-3 and LMP-2 genes. This strain was characterized in patients with short post-diagnosis survival. The main mutations found during AITL on the most mutated latency or tegument genes were identified and discussed. We showed that the virus was clonal in all the AITL samples, suggesting that it may be involved in this pathology. Additionally, EBV was latent in all the AITL samples; for one sample only, the virus was found to be latent and probably replicative, depending on the cells. These various elements support the role of EBV in AITL.

Human herpesvirus diversity is altered in HLA class I binding peptides

Herpesviruses are ubiquitous, genetically diverse DNA viruses, with long-term presence in humans associated with infrequent but significant pathology. Human leukocyte antigen (HLA) class I presents intracellularly derived peptide fragments from infected tissue cells to CD8+ T and natural killer cells, thereby directing antiviral immunity. Allotypes of highly polymorphic HLA class I are distinguished by their peptide binding repertoires. Because this HLA class I variation is a major determinant of herpesvirus disease, we examined if sequence diversity of virus proteins reflects evasion of HLA presentation. Using population genomic data from Epstein–Barr virus (EBV), human cytomegalovirus (HCMV), and Varicella–Zoster virus, we tested whether diversity differed between the regions of herpesvirus proteins that can be recognized, or not, by HLA class I. Herpesviruses exhibit lytic and latent infection stages, with the latter better enabling immune evasion. Whereas HLA binding peptides of lytic proteins are conserved, we found that EBV and HCMV proteins expressed during latency have increased peptide sequence diversity. Similarly, latent, but not lytic, herpesvirus proteins have greater population structure in HLA binding than nonbinding peptides. Finally, we found patterns consistent with EBV adaption to the local HLA environment, with less efficient recognition of EBV isolates by high-frequency HLA class I allotypes. Here, the frequency of CD8+ T cell epitopes inversely correlated with the frequency of HLA class I recognition. Previous analyses have shown that pathogen-mediated natural selection maintains exceptional polymorphism in HLA residues that determine peptide recognition. Here, we show that HLA class I peptide recognition impacts diversity of globally widespread pathogens.

Reduced IRF4 expression promotes lytic phenotype in Type 2 EBV-infected B cells

Humans are infected with two types of EBV (Type 1 (T1) and Type 2 (T2)) that differ substantially in their EBNA2 and EBNA 3A/B/C latency proteins and have different phenotypes in B cells. T1 EBV transforms B cells more efficiently than T2 EBV in vitro, and T2 EBV-infected B cells are more lytic. We previously showed that both increased NFATc1/c2 activity, and an NFAT-binding motif within the BZLF1 immediate-early promoter variant (Zp-V3) contained in all T2 strains, contribute to lytic infection in T2 EBV-infected B cells. Here we compare cellular and viral gene expression in early-passage lymphoblastoid cell lines (LCLs) infected with either T1 or T2 EBV strains. Using bulk RNA-seq, we show that T2 LCLs are readily distinguishable from T1 LCLs, with approximately 600 differentially expressed cellular genes. Gene Set Enrichment Analysis (GSEA) suggests that T2 LCLs have increased B-cell receptor (BCR) signaling, NFAT activation, and enhanced expression of epithelial-mesenchymal-transition-associated genes. T2 LCLs also have decreased RNA and protein expression of a cellular gene required for survival of T1 LCLs, IRF4. In addition to its essential role in plasma cell differentiation, IRF4 decreases BCR signaling. Knock-down of IRF4 in a T1 LCL (infected with the Zp-V3-containing Akata strain) induced lytic reactivation whereas over-expression of IRF4 in Burkitt lymphoma cells inhibited both NFATc1 and NFATc2 expression and lytic EBV reactivation. Single-cell RNA-seq confirmed that T2 LCLs have many more lytic cells compared to T1 LCLs and showed that lytically infected cells have both increased NFATc1, and decreased IRF4, compared to latently infected cells. These studies reveal numerous differences in cellular gene expression in B cells infected with T1 versus T2 EBV and suggest that decreased IRF4 contributes to both the latent and lytic phenotypes in cells with T2 EBV.

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.

Do Epstein–Barr Virus Mutations and Natural Genome Sequence Variations Contribute to Disease?

Most of the world’s population is infected by the Epstein–Barr virus (EBV), but the incidence of the diseases associated with EBV infection differs greatly in different parts of the world. Many factors may determine those differences, but variation in the virus genome is likely to be a contributing factor for some of the diseases. Here, we describe the main forms of EBV genome sequence variation, and the mechanisms by which variations in the virus genome are likely to contribute to disease. EBV genome deletions or polymorphisms can also provide useful markers for monitoring disease. If some EBV strains prove to be more pathogenic than others, this suggests the possible value of immunising people against infection by those pathogenic strains.

Epstein-Barr virus miR-BHRF1-3 targets the BZLF1 3'UTR and regulates the lytic cycle

Suppression of lytic viral gene expression is a key aspect of the Epstein-Barr virus (EBV) life cycle to facilitate the establishment of latent infection. Molecular mechanisms regulating transitions between EBV lytic replication and latency are not fully understood. Here, we investigated the impact of viral microRNAs on the EBV lytic cycle. Through functional assays, we found that miR-BHRF1-3 attenuates EBV lytic gene expression following reactivation. To understand the miRNA targets contributing to this activity, we performed Ago PAR-CLIP analysis on EBV-positive, reactivated Burkitt's lymphoma cells and identified multiple miR-BHRF1-3 interactions with viral transcripts. Using luciferase reporter assays, we confirmed a miRNA interaction site within the 3'UTR of BZLF1 which encodes the essential immediate early (IE) transactivator Zta. Comparison of >850 published EBV genomes identified sequence polymorphisms within the miR-BHRF1-3 locus that deleteriously affect miRNA expression and function. Molecular interactions between the homologous viral miRNA, miR-rL1-17, and IE transcripts encoded by rhesus lymphocryptovirus were further identified. Our data demonstrate that regulation of IE gene expression by a BHRF1 miRNA is conserved amongst lymphocryptoviruses, and further reveal virally-encoded genetic elements that orchestrate viral antigen expression during the lytic cycle. Importance Epstein-Barr virus infection is predominantly latent in healthy individuals, while periodic cycles of reactivation are thought to facilitate persistent lifelong infection. Lytic infection has been linked to development of certain EBV-associated diseases. Here, we demonstrate that EBV miR-BHRF1-3 can suppress lytic replication by directly inhibiting Zta expression. Moreover, we identify nucleotide variants that impact the function of miR-BHRF1-3, which may contribute to specific EBV pathologies.

Comprehensive Evolutionary Analysis of Complete Epstein-Barr Virus Genomes from Argentina and Other Geographies

The sequence variability of the Epstein–Barr virus has been extensively studied throughout previous years in isolates from various geographic regions and consequent variations at both genetic and genomic levels have been described. However, isolates from South America were underrepresented in these studies. Here, we sequenced 15 complete EBV genomes that we analyzed together with publicly available raw NGS data for 199 EBV isolates from other parts of the globe by means of a custom-built bioinformatic pipeline. The phylogenetic relations of the genomes, the geographic structure and variability of the data set, and the evolution rates for the whole genome and each gene were assessed. The present work contributes to overcoming the scarcity of complete EBV genomes from South America and is the most comprehensive geography-related variability study, which involved determining the actual contribution of each EBV gene to the geographic segregation of the entire genome. Moreover, to the best of our knowledge, we established for the first time the evolution rate for the entire EBV genome based on a host–virus codivergence-independent assumption and assessed their evolution rates on a gene-by-gene basis, which were related to the encoded protein function. Considering the evolution of dsDNA viruses with a codivergence-independent approach may lay the basis for future research on EBV evolution. The exhaustive bioinformatic analysis performed on this new dataset allowed us to draw a novel set of conclusions regarding the genome evolution of EBV.

Meet the Insidious Players: Review of Viral Infections in Head and Neck Cancer Etiology with an Update on Clinical Trials

Head and neck cancers (HNC) occur in the upper aerodigestive tract and are among the most common cancers. The etiology of HNC is complex, involving many factors, including excessive tobacco and alcohol consumption; over the last two decades, oncogenic viruses have also been recognized as an important cause of HNC. Major etiological agents of nasopharynx carcinoma and oropharyngeal carcinoma include Epstein-Barr virus (EBV) and human papillomaviruses (HPVs), both of which are able to interfere with cell cycle control. Additionally, the association of hepatitis C and hepatitis B infection was observed in oral cavity, oropharyngeal, laryngeal, and nasopharyngeal cancers. Overall prognoses depend on anatomic site, stage, and viral status. Current treatment options, including radiotherapy, chemotherapy, targeted therapies and immunotherapies, are distributed in order to improve overall patient prognosis and survival rates. However, the interplay between viral genome sequences and the health, disease, geography, and ethnicity of the host are crucial for understanding the role of viruses and for development of potential personalized treatment and prevention strategies. This review provides the most comprehensive analysis to date of a vast field, including HNC risk factors, as well as viral mechanisms of infection and their role in HNC development. Additionally, currently available treatment options investigated through clinical practice are emphasized in the paper.

T cell epitope screening of Epstein-Barr virus fusion protein gB

Glycoprotein B (gB) is an essential fusion protein for the Epstein-Barr virus (EBV) infection of both B cells and epithelial cells and is thus a promising target antigen for a prophylactic vaccine to prevent or reduce EBV-associated disease. T cell responses play key roles in the control of persistent EBV infection and in the efficacy of a vaccine. However, to date, T cell responses to gB have been characterized for only a limited number of human leukocyte antigen (HLA) alleles. Here, we screened gB T cell epitopes in 23 healthy EBV carriers and ten patients with nasopharyngeal cancer (NPC) using a peptide library spanning the entire gB sequence. We identified twelve novel epitopes in the context of seven new HLA restrictions that are common in Asian populations. Two epitopes, gB_214-223 and gB_840-849, restricted by HLA-B*58:01 and B*38:02, respectively, elicited specific CD8⁺ T cell responses to inhibit EBV-driven B cell transformation. Interestingly, gB-specific CD8⁺ T cells were more frequent in healthy viral carriers with EBV reactivation than in those without EBV reactivation, indicating that EBV reactivation in vivo stimulates both humoral (VCA-gp125-IgA) and cellular responses to gB. We further found that most gB epitopes are conserved among different EBV strains. Our study broadens the diversity and HLA restrictions of gB epitopes and suggests that gB is a common target of T cell responses in healthy viral carriers with EBV reactivation. In particular, the precisely mapped and conserved gB epitopes provide valuable information for prophylactic vaccine development.ImportanceT cells are crucial for the control of persistent EBV infection and the development of EBV-associated diseases. The EBV gB protein is essential for virus entry into B cells and epithelial cells and is thus a target antigen for vaccine development. Understanding T cell responses to gB is important for subunit vaccine design. Herein, we comprehensively characterized T cell responses to full-length gB. Our results expand the available gB epitopes and HLA restrictions, particularly those common in Asian populations. Furthermore, we showed that gB-specific CD8⁺ T cells inhibit B cell transformation ex vivo and that gB-specific CD8⁺ T cell responses in vivo may be associated with intermittent EBV reactivation in asymptomatic viral carriers. These gB epitopes are highly conserved among geographically separated EBV strains. Precisely mapped and conserved T cell epitopes may contribute to immune monitoring and to the development of a gB subunit vaccine.

The influence of human genetic variation on Epstein-Barr virus sequence diversity

Epstein-Barr virus (EBV) is one of the most common viruses latently infecting humans. Little is known about the impact of human genetic variation on the large inter-individual differences observed in response to EBV infection. To search for a potential imprint of host genomic variation on the EBV sequence, we jointly analyzed paired viral and human genomic data from 268 HIV-coinfected individuals with CD4 + T cell count < 200/mm3 and elevated EBV viremia. We hypothesized that the reactivated virus circulating in these patients could carry sequence variants acquired during primary EBV infection, thereby providing a snapshot of early adaptation to the pressure exerted on EBV by the individual immune response. We searched for associations between host and pathogen genetic variants, taking into account human and EBV population structure. Our analyses revealed significant associations between human and EBV sequence variation. Three polymorphic regions in the human genome were found to be associated with EBV variation: one at the amino acid level (BRLF1:p.Lys316Glu); and two at the gene level (burden testing of rare variants in BALF5 and BBRF1). Our findings confirm that jointly analyzing host and pathogen genomes can identify sites of genomic interactions, which could help dissect pathogenic mechanisms and suggest new therapeutic avenues.

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.

Cytomegalovirus, Epstein-Barr virus, and human herpesvirus-6 infections in patients with myalgic еncephalomyelitis/chronic fatigue syndrome

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a disabling multisystem chronic disease. The etiology and pathogenesis of ME/CFS are unknown. Infections of cytomegalovirus (CMV), Epstein-Barr virus (EBV), and human herpesvirus-6 (HHV-6) are suspected as etiological agents for ME/CFS. This study aims to estimate prevalence and type (active/latent) of EBV, CMV, and HHV-6 infections in Bulgarian ME/CFS patients. In the study were included 58 patients with ME/CFS and 50 healthy controls. Virus-specific antibodies were detected by enzyme-linked immunosorbent assay and viral genomic sequences in peripheral blood mononuclear cell (PBMCs) and plasma samples by nested polymerase chain reaction (PCR). We did not observe any significant differences in virus-specific immunoglobulin G and immunoglobulin M positivity rates between patients with ME/CFS and control group. In ME/CFS plasma samples, EBV DNA was found in 24.1%, CMV DNA in 3.4%, and HHV-6 DNA in 1.7% of samples. EBV DNA was detected in 4%, and CMV and HHV-6 DNA were not found in plasma samples of controls. The frequency of viral genome detection in PBMCs of patients and controls was 74% vs 78% for CMV, 81% vs 84% for EBV, and 82.8% vs 82% for HHV-6. The difference in frequency of EBV active infection in ME/CFS and control group was statistically significant (P = .0027). No ME/CFS and control individuals with active CMV and HHV-6 infection were observed. In conclusion, this study using both serological and PCR-based techniques for distinguishing between active and latent infection showed high rate of active EBV infection among patients with ME/CFS indicating that at least in a subset of cases, EBV is important factor for the development of disease.

Sequencing Analysis of Plasma Epstein-Barr Virus DNA Reveals Nasopharyngeal Carcinoma-Associated Single Nucleotide Variant Profiles

BACKGROUND

Nasopharyngeal carcinoma (NPC) is strongly associated with Epstein-Barr virus (EBV) infection. Plasma EBV DNA is a validated screening tool for NPC. In screening, there are some individuals who do not have NPC but carry EBV DNA in plasma. Currently it is not known from screening if there may be any genotypic differences in EBV isolates from NPC and non-NPC subjects. Also, low concentrations of EBV DNA in plasma could pose challenge to such EBV genotypic analysis through plasma DNA sequencing.

METHODS

In a training dataset comprised of plasma DNA sequencing data of NPC and non-NPC subjects, we studied the difference in the EBV single nucleotide variant (SNV) profiles between the two groups. The most differentiating SNVs across the EBV genome were identified. We proposed an NPC risk score to be derived from the genotypic patterns over these SNV sites. We subsequently analyzed the NPC risk scores in a testing set.

RESULTS

A total of 661 significant SNVs across the EBV genome were identified from the training set. In the testing set, NPC plasma samples were shown to have high NPC risk scores, which suggested the presence of NPC-associated EBV SNV profiles. Among the non-NPC samples, there was a wide range of NPC risk scores. These results support the presence of diverse SNV profiles of EBV isolates from non-NPC subjects.

CONCLUSION

EBV genotypic analysis is feasible through plasma DNA sequencing. The NPC risk score may be used to inform the cancer risk based on the EBV genome-wide SNV profile.

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.

Epstein-Barr Virus Genomes Reveal Population Structure and Type 1 Association with Endemic Burkitt Lymphoma

Endemic Burkitt lymphoma (eBL), the most prevalent pediatric cancer in sub-Saharan Africa, is distinguished by its inclusion of Epstein-Barr virus (EBV). In order to better understand the impact of EBV variation in eBL tumorigenesis, we improved viral DNA enrichment methods and generated a total of 98 new EBV genomes from both eBL cases (n = 58) and healthy controls (n = 40) residing in the same geographic region in Kenya. Using our unbiased methods, we found that EBV type 1 was significantly more prevalent in eBL patients (74.5%) than in healthy children (47.5%) (odds ratio = 3.24, 95% confidence interval = 1.36 to 7.71, P = 0.007), as opposed to similar proportions in both groups. Controlling for EBV type, we also performed a genome-wide association study identifying six nonsynonymous variants in the genes EBNA1, EBNA2, BcLF1, and BARF1 that were enriched in eBL patients. In addition, viruses isolated from plasma of eBL patients were identical to their tumor counterparts consistent with circulating viral DNA originating from the tumor. We also detected three intertypic recombinants carrying type 1 EBNA2 and type 2 EBNA3 regions, as well as one novel genome with a 20-kb deletion, resulting in the loss of multiple lytic and virion genes. Comparing EBV types, viral genes displayed differential variation rates as type 1 appeared to be more divergent, while type 2 demonstrated novel substructures. Overall, our findings highlight the complexities of the EBV population structure and provide new insight into viral variation, potentially deepening our understanding of eBL oncogenesis.IMPORTANCE Improved viral enrichment methods conclusively demonstrate EBV type 1 to be more prevalent in eBL patients than in geographically matched healthy controls, which previously underrepresented the prevalence of EBV type 2. Genome-wide association analysis between cases and controls identifies six eBL-associated nonsynonymous variants in EBNA1, EBNA2, BcLF1, and BARF1 genes. Analysis of population structure reveals that EBV type 2 exists as two genomic subgroups and was more commonly found in female than in male eBL patients.

CD21 (Complement Receptor 2) Is the Receptor for Epstein-Barr Virus Entry into T Cells

Epstein-Barr virus (EBV) is associated with a number of T-cell diseases, including some peripheral T-cell lymphomas, hemophagocytic lymphohistiocytosis, and chronic active EBV disease. The tropism of EBV for B cells and epithelial cell infection has been well characterized, but infection of T cells has been minimally explored. We have recently shown that the EBV type 2 (EBV-2) strain has the unique ability to infect mature T cells. Utilizing an ex vivo infection model, we sought to understand the viral glycoprotein and cellular receptor required for EBV-2 infection of T cells. Here, using a neutralizing-antibody assay, we found that viral gp350 and complement receptor 2 (CD21) are required for CD3⁺ T-cell infection. Using the HB5 anti-CD21 antibody clone but not the Bly-4 anti-CD21 antibody clone, we detected expression of CD21 on both CD4⁺ and CD8⁺ T cells, with the highest expression on naive CD4 and CD8⁺ T-cell subsets. Using CRISPR to knock out CD21, we demonstrated that CD21 is necessary for EBV entry into the Jurkat T-cell line. Together, these results indicate that EBV uses the same viral glycoprotein and cellular receptor for both T- and B-cell infection.IMPORTANCE Epstein-Barr virus (EBV) has a well-described tropism for B cells and epithelial cells. Recently, we described the ability of a second strain of EBV, EBV type 2, to infect mature peripheral T cells. Using a neutralizing antibody assay, we determined that EBV uses the viral glycoprotein gp350 and the cellular protein CD21 to gain entry into mature peripheral T cells. CRISPR-Cas9 deletion of CD21 on the Jurkat T-cell line confirmed that CD21 is required for EBV infection. This study has broad implications, as we have defined a function for CD21 on mature peripheral T cells, i.e., as a receptor for EBV. In addition, the requirement for gp350 for T-cell entry has implications for EBV vaccine studies currently targeting the gp350 glycoprotein to prevent EBV-associated diseases.

miRNAs: EBV Mechanism for Escaping Host's Immune Response and Supporting Tumorigenesis

Epstein-Barr virus (EBV) or human herpesvirus 4 (HHV-4) is a ubiquitous human oncogenic virus, and the first human virus found to express microRNAs (miRNAs). Its genome contains two regions encoding more than 40 miRNAs that regulate expression of both viral and human genes. There are numerous evidences that EBV miRNAs impact immune response, affect antigen presentation and recognition, change T- and B-cell communication, drive antibody production during infection, and have a role in cell apoptosis. Moreover, the ability of EBV to induce B-cell transformation and take part in mechanisms of oncogenesis in humans is well known. Although EBV infection is associated with development of various diseases, the role of its miRNAs is still not understood. There is abundant data describing EBV miRNAs in nasopharyngeal carcinoma and several studies that have tried to evaluate their role in gastric carcinoma and lymphoma. This review aims to summarize so far known data about the role of EBV miRNAs in altered regulation of gene expression in human cells in EBV-associated diseases.

Genomic variations in EBNA3C of EBV associate with posttransplant lymphoproliferative disorder

Epstein-Barr Virus (EBV) is a ubiquitous virus linked to a variety of lymphoid and epithelial malignancies. In solid organ and hematopoietic stem cell transplant recipients, EBV is causally associated with posttransplant lymphoproliferative disorder (PTLD), a group of heterogeneous lymphoid diseases. EBV+ B cell lymphomas that develop in the context of PTLD are generally attributed to the immunosuppression required to promote graft survival, but little is known regarding the role of EBV genome diversity in the development of malignancy. We deep-sequenced the EBV genome from the peripheral blood of 18 solid organ transplant recipients, including 6 PTLD patients. Sequences from 6 EBV+ spontaneous lymphoblastoid B cell lines (SLCL) were similarly analyzed. The EBV genome from PTLD patients had a significantly greater number of variations than EBV from transplant recipients without PTLD. Importantly, there were 15 nonsynonymous variations, including 8 in the latent cycle gene EBNA3C that were associated with the development of PTLD. One of the nonsynonymous variations in EBNA3C is located within a previously defined T cell epitope. These findings suggest that variations in the EBV genome can contribute to the pathogenesis of PTLD.

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

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

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

Structure-based design of small-molecule inhibitors of EBNA1 DNA binding blocks Epstein-Barr virus latent infection and tumor growth

Epstein-Barr virus (EBV) is a DNA tumor virus responsible for 1 to 2% of human cancers including subtypes of Burkitt's lymphoma, Hodgkin's lymphoma, gastric carcinoma, and nasopharyngeal carcinoma (NPC). Persistent latent infection drives EBV-associated tumorigenesis. Epstein-Barr nuclear antigen 1 (EBNA1) is the only viral protein consistently expressed in all EBV-associated tumors and is therefore an attractive target for therapeutic intervention. It is a multifunctional DNA binding protein critical for viral replication, genome maintenance, viral gene expression, and host cell survival. Using a fragment-based approach and x-ray crystallography, we identify a 2,3-disubstituted benzoic acid series that selectively inhibits the DNA binding activity of EBNA1. We characterize these inhibitors biochemically and in cell-based assays, including chromatin immunoprecipitation and DNA replication assays. In addition, we demonstrate the potency of EBNA1 inhibitors to suppress tumor growth in several EBV-dependent xenograft models, including patient-derived xenografts for NPC. These inhibitors selectively block EBV gene transcription and alter the cellular transforming growth factor-β (TGF-β) signaling pathway in NPC tumor xenografts. These EBNA1-specific inhibitors show favorable pharmacological properties and have the potential to be further developed for the treatment of EBV-associated malignancies.

B cells infected with Type 2 Epstein-Barr virus (EBV) have increased NFATc1/NFATc2 activity and enhanced lytic gene expression in comparison to Type 1 EBV infection

Humans are infected with two distinct strains (Type 1 (T1) and Type 2 (T2)) of Epstein-Barr virus (EBV) that differ substantially in their EBNA2 and EBNA 3A/B/C latency genes and the ability to transform B cells in vitro. While most T1 EBV strains contain the "prototype" form of the BZLF1 immediate-early promoter ("Zp-P"), all T2 strains contain the "Zp-V3" variant, which contains an NFAT binding motif and is activated much more strongly by B-cell receptor signalling. Whether B cells infected with T2 EBV are more lytic than cells infected with T1 EBV is unknown. Here we show that B cells infected with T2 EBV strains (AG876 and BL5) have much more lytic protein expression compared to B cells infected with T1 EBV strains (M81, Akata, and Mutu) in both a cord blood-humanized (CBH) mouse model and EBV-transformed lymphoblastoid cell lines (LCLs). Although T2 LCLs grow more slowly than T1 LCLs, both EBV types induce B-cell lymphomas in CBH mice. T1 EBV strains (M81 and Akata) containing Zp-V3 are less lytic than T2 EBV strains, suggesting that Zp-V3 is not sufficient to confer a lytic phenotype. Instead, we find that T2 LCLs express much higher levels of activated NFATc1 and NFATc2, and that cyclosporine (an NFAT inhibitor) and knockdown of NFATc2 attenuate constitutive lytic infection in T2 LCLs. Both NFATc1 and NFATc2 induce lytic EBV gene expression when combined with activated CAMKIV (which is activated by calcium signaling and activates MEF2D) in Burkitt Akata cells. Together, these results suggest that B cells infected with T2 EBV are more lytic due to increased activity of the cellular NFATc1/c2 transcription factors in addition to the universal presence of the Zp-V3 form of BZLF1 promoter.

Reactivation of Epstein-Barr virus by a dual-responsive fluorescent EBNA1-targeting agent with Zn2+-chelating function

EBNA1 is the only Epstein–Barr virus (EBV) latent protein responsible for viral genome maintenance and is expressed in all EBV-infected cells. Zn²⁺ is essential for oligomerization of the functional EBNA1. We constructed an EBNA1 binding peptide with a Zn²⁺ chelator to create an EBNA1-specific inhibitor (ZRL₅P₄). ZRL₅P₄ by itself is sufficient to reactivate EBV from its latent infection. ZRL₅P₄ is able to emit unique responsive fluorescent signals once it binds with EBNA1 and a Zn²⁺ ion. ZRL₅P₄ can selectively disrupt the EBNA1 oligomerization and cause nasopharyngeal carcinoma (NPC) tumor shrinkage, possibly due to EBV lytic induction. Dicer1 seems essential for this lytic reactivation. As can been seen, EBNA1 is likely to maintain NPC cell survival by suppressing viral reactivation.

Epstein–Barr nuclear antigen 1 (EBNA1) plays a vital role in the maintenance of the viral genome and is the only viral protein expressed in nearly all forms of Epstein–Barr virus (EBV) latency and EBV-associated diseases, including numerous cancer types. To our knowledge, no specific agent against EBV genes or proteins has been established to target EBV lytic reactivation. Here we report an EBNA1- and Zn²⁺-responsive probe (ZRL₅P₄) which alone could reactivate the EBV lytic cycle through specific disruption of EBNA1. We have utilized the Zn²⁺ chelator to further interfere with the higher order of EBNA1 self-association. The bioprobe ZRL₅P₄ can respond independently to its interactions with Zn²⁺ and EBNA1 with different fluorescence changes. It can selectively enter the nuclei of EBV-positive cells and disrupt the oligomerization and oriP-enhanced transactivation of EBNA1. ZRL₅P₄ can also specifically enhance Dicer1 and PML expression, molecular events which had been reported to occur after the depletion of EBNA1 expression. Importantly, we found that treatment with ZRL₅P₄ alone could reactivate EBV lytic induction by expressing the early and late EBV lytic genes/proteins. Lytic induction is likely mediated by disruption of EBNA1 oligomerization and the subsequent change of Dicer1 expression. Our probe ZRL₅P₄ is an EBV protein-specific agent that potently reactivates EBV from latency, leading to the shrinkage of EBV-positive tumors, and our study also suggests the association of EBNA1 oligomerization with the maintenance of EBV latency.

Antibodies to Human Herpesviruses in Myalgic Encephalomyelitis/Chronic Fatigue Syndrome Patients

Myalgic encephalomyelitis, also referred to as chronic fatigue syndrome (ME/CFS) is a debilitating disease characterized by myalgia and a sometimes severe limitation of physical activity and cognition. It is exacerbated by physical and mental activity. Its cause is unknown, but frequently starts with an infection. The eliciting infection (commonly infectious mononucleosis or an upper respiratory infection) can be more or less well diagnosed. Among the human herpesviruses (HHV-1-8), HHV-4 (Epstein-Barr virus; EBV), HHV-6 (including HHV-6A and HHV-6B), and HHV-7, have been implicated in the pathogenesis of ME/CFS. It was therefore logical to search for serological evidence of past herpesvirus infection/reactivation in several cohorts of ME/CFS patients (all diagnosed using the Canada criteria). Control samples were from Swedish blood donors. We used whole purified virus, recombinant proteins, and synthetic peptides as antigens in a suspension multiplex immunoassay (SMIA) for immunoglobulin G (IgG). The study on herpesviral peptides based on antigenicity with human sera yielded novel epitope information. Overall, IgG anti-herpes-viral reactivities of ME/CFS patients and controls did not show significant differences. However, the high precision and internally controlled format allowed us to observe minor relative differences between antibody reactivities of some herpesviral antigens in ME/CFS versus controls. ME/CFS samples reacted somewhat differently from controls with whole virus HHV-1 antigens and recombinant EBV EBNA6 and EA antigens. We conclude that ME/CFS samples had similar levels of IgG reactivity as blood donor samples with HHV-1-7 antigens. The subtle serological differences should not be over-interpreted, but they may indicate that the immune system of some ME/CFS patients interact with the ubiquitous herpesviruses in a way different from that of healthy controls.

F factor plasmid-mediated Epstein-Barr virus genome introduction establishes an EBV positive NPC cell model

Background

Most Epstein-Barr virus (EBV)-positive cells lose the EBV episomes upon prolonged propagation.

Purpose

The purposes of this study were to establish a simple cell model for nasopharyngeal carcinoma (NPC) research by introducing a plasmid with the EBV genome into NPC cells and then to investigate the resulting changes in malignant biological behaviour and NPC-associated signalling pathways.

Methods

HONE1 NPC cells were transfected with F-factor plasmids including the EBV genome (HONE1-EBV cells). Then cell proliferation, migration, cell cycle distribution and apoptosis were evaluated in vitro by using CCK8, transwell and flow cytometry assays respectively. EBV-encoded proteins and cell signal tranducting proteins were detected by western blot assays. EBV-encoded RNAs were detected by in situ hybridization. EBV particles were assayed by transmission electron microscope (TEM). The morphology of cells were detected by immunofluorescence assays for alpha-tubulin.

Results

Latent membrane protein 1 (LMP1), latent membrane protein 2A (LMP2A), Epstein-Barr nuclear antigen 1 (EBNA1) and EBV-encoded small RNAs (EBERs) were successfully expressed in HONE1-EBV cells. No EBV particles were founded by TEM. Introduction of the EBV genome significantly promoted proliferation, cell cycle progression and migration and inhibited apoptosis in HONE1 cells. Immunofluorescence assays showed that the morphology of HONE1-EBV cells changed into spindle. Furthermore, EBV genome introduction significantly inhibited the JAK/STAT signalling pathway, while it activated the PI3K-AKT and NF-κB signalling pathways in HONE1 cells.

Conclusion

These findings suggest that F-factor plasmid-mediated EBV genome introduction was successful in constructing an EBV positive cell model, which showed deteriorated biological behavior and activated NPC-associated signalling pathways. This model can serve as a good tool for studying EBV in NPC, but the subtle differences in cancer-associated pathways must be considered.

The Function and Therapeutic Potential of Epstein-Barr Virus-Encoded MicroRNAs in Cancer

Epstein-Barr virus (EBV) is a ubiquitous human γ-herpesvirus that infects over 90% of the global population. EBV is considered a contributory factor in a variety of malignancies including nasopharyngeal carcinoma, gastric carcinoma, Burkitt lymphoma, and Hodgkin’s lymphoma. Notably, EBV was the first virus found to encode microRNAs (miRNAs). Increasing evidence indicates that EBV-encoded miRNAs contribute to the carcinogenesis and development of EBV-associated malignancies. EBV miRNAs have been shown to inhibit the expression of genes involved in cell proliferation, apoptosis, invasion, and immune signaling pathways. Therefore, EBV miRNAs perform a significant function in the complex host-virus interaction and EBV-driven carcinogenesis. However, the integrated mechanisms underlying the roles of EBV miRNAs in carcinogenesis remain to be further explored. In this review, we describe recent advances regarding the involvement of EBV miRNAs in the pathogenesis of EBV-associated malignancies and discuss their potential utility as cancer biomarkers. An in-depth investigation into the pro-carcinogenic role of EBV miRNAs will expand our knowledge of the biological processes associated with virus-driven tumors and contribute to the development of novel therapeutic strategies for the treatment of EBV-associated malignancies.

A reliable Epstein-Barr Virus classification based on phylogenomic and population analyses

The Epstein-Barr virus (EBV) infects more than 90% of the human population, playing a key role in the origin and progression of malignant and non-malignant diseases. Many attempts have been made to classify EBV according to clinical or epidemiological information; however, these classifications show frequent incongruences. For instance, they use a small subset of genes for sorting strains but fail to consider the enormous genomic variability and abundant recombinant regions present in the EBV genome. These could lead to diversity overestimation, alter the tree topology and misinterpret viral types when classified, therefore, a reliable EBV phylogenetic classification is needed to minimize recombination signals. Recombination events occur 2.5-times more often than mutation events, suggesting that recombination has a much stronger impact than mutation in EBV genomic diversity, detected within common ancestral node positions. The Hierarchical Bayesian Analysis of Population Structure (hierBAPS) resulted in the differentiation of 12 EBV populations showed seven monophyletic and five paraphyletic. The populations identified were related to geographic location, of which three populations (EBV-p1/Asia/GC, EBV-p2/Asia II/Tumors and EBV-p4/China/NPC) were related to tumor development. Therefore, we proposed a new consistent and non-simplistic EBV classification, beneficial in minimizing the recombination signal in the phylogeny reconstruction, investigating geography relationship and even infer associations to human diseases. These EBV classifications could also be useful in developing diagnostic applications or defining which strains need epidemiological surveillance.

Clinical characteristics and prognostic value of pre-retreatment plasma epstein-barr virus DNA in locoregional recurrent nasopharyngeal carcinoma

PURPOSE

To define the clinical characteristics and prognostic value of pre-retreatment plasma Epstein-Barr virus (EBV) DNA, we investigated EBV status in locoregional recurrent nasopharyngeal carcinoma (lrNPC) patients.

METHODS

Between April 2008 and August 2016, the data of patients with nonmetastatic lrNPC were retrospectively reviewed. The survival indexes of patients between different pre-retreatment EBV status groups were compared.

RESULTS

A total of 401 patients with nonmetastatic lrNPC were enrolled, and 197 (49.1%) patients had detectable pre-retreatment plasma EBV DNA. Treatment included radiotherapy alone (n = 37 patients), surgery alone (n = 105), radiotherapy (n = 208), surgery combined with radiotherapy (n = 20), chemotherapy and targeted therapy (n = 31). Median follow-up was 32 months. The 3-year locoregional relapse-free survival (LRRFS), distant metastasis-free survival (DMFS), and overall survival (OS) rates for the entire cohort were 64.8%, 89.4%, and 58.8%, respectively. The estimated 3-year LRRFS, DMFS, and OS rates for the pre EBV-positive group vs the pre EBV-negative group were 54.2% vs 75.0% (P < 0.001), 86.6% vs 91.9% (P = 0.05), 51.6% vs 65.9% (P = 0.01), respectively. Among patients in the clinical stage rI/II, there were 17 patients in the radiotherapy alone group and 49 patients in the surgery alone group. And there was no significant difference in overall survival between radiotherapy and surgery, even among the different pre-EBV statuses (P > 0.05). In terms of long-term toxic and side effects, the incidence of radioactive temporal lobe injury in the radiotherapy group was higher than that in the surgery group (35.3% vs 8.2%, P < 0.001), and no statistically significant difference was found in other long-term toxic and side effects.

CONCLUSIONS

The positive rate of pre-retreatment plasma EBV DNA in lrNPC is lower than primary NPC. The prognosis of EBV DNA negative group is better than positive group. For locally early-stage lrNPC, regardless of EBV DNA status, radiotherapy and surgery are available options and both can achieve better long-term survival.

Genome sequencing analysis identifies Epstein-Barr virus subtypes associated with high risk of nasopharyngeal carcinoma

Epstein-Barr virus (EBV) infection is ubiquitous worldwide and is associated with multiple cancers, including nasopharyngeal carcinoma (NPC). The importance of EBV viral genomic variation in NPC development and its striking epidemic in southern China has been poorly explored. Through large-scale genome sequencing of 270 EBV isolates and two-stage association study of EBV isolates from China, we identified two non-synonymous EBV variants within BALF2 strongly associated with the risk of NPC (odds ratio (OR) = 8.69, P=9.69×10⁻²⁵ for SNP 162476_C; OR = 6.14, P=2.40×10⁻³² for SNP 163364_T). The cumulative effects of these variants contributed to 83% of the overall risk of NPC in southern China. Phylogenetic analysis of the risk variants revealed a unique origin in Asia, followed by clonal expansion in NPC-endemic regions. Our results provide novel insights into NPC endemic in southern China and also enable the identification of high-risk individuals for NPC prevention.

Whole-genome sequencing and association analysis of 270 Epstein-Barr virus (EBV) isolates from China identify two non-synonymous EBV variants within BALF2 strongly associated with the risk of nasopharyngeal carcinoma.

Prevalence of human herpesviruses infections in nonmalignant tonsils: The SPLIT study

OBJECTIVE

To assess the prevalence of all known human herpesviruses (HHV) in tonsils of an age-stratified large sample of immunocompetent children and adults.

METHODS

Patients undergoing tonsillectomy for benign indications were recruited in 19 French hospitals. After resection, the entire outer surfaces of right and left half tonsils were extensively brushed. A highly sensitive species-specific multiplex assay was used to detect herpes simplex virus 1 (HSV1), HSV2, Epstein-Barr virus (EBV; types 1 and 2), and human cytomegalovirus (CMV) DNA in 688, as well as varicella zoster virus (VZV), HHV6A, HHV6B, HHV7, and Kaposi's sarcoma-associated herpesvirus (KSHV) DNA in a subset of 440 tonsil brushings.

RESULTS

Overall 85% of tonsil brushing samples were infected with at least one HHV species. HHV7 and EBV were the most prevalent (≈70%), followed by HHV6B (≈50%), HSV1, CMV, VZV (≈2%), and KSHV and HSV2 (<1%), while HHV6A was not detected. EBV prevalence was significantly higher in adults than in children, whereas it was opposite for HHV6B and VZV. No difference in HHV prevalence was observed by sex. In multivariate analysis, EBV detection was associated with age greater than or equal to 15 years (prevalence ratio [PR] = 1.8; 95% confidence interval [CI]: 1.5-2.3) and marginally with tobacco smoking (PR = 1.2; 95% CI: 1.1-1.3).

CONCLUSION

Differing patterns of HHV infection in tonsils in a large age-stratified population were described. This study is by far the largest available and shows that EBV, HHV6B, and HHV7 are commonly detected in the tonsils in both men and women, in contrast to other HHVs.

Epstein-Barr virus latent gene EBNA-1 genetic diversity among transplant patients compared with patients with infectious mononucleosis

INTRODUCTION

As a step toward evaluating the association between Epstein-Barr virus genetic diversity and post-transplant lymphoproliferative disorder (PTLD), we conducted a preliminary study to compare the genetic diversity of the EBNA-1 gene among transplant patients and patients with infectious mononucleosis (IM).

METHODS

We sequenced the EBNA-1 gene in blood samples from study subjects using Sanger methodology. The sequences were aligned with a reference strain and compared with publicly available sequences.

RESULTS

We analyzed 33 study samples and 25 publicly available sequences along with the reference strain B95-8. The evaluable samples were from sixteen patients with IM (median age 14.0 years, range 2-24) and 17 transplant patients. There were six children without PTLD (median age 1.93 years, range 0.79-7.46) and 11 who developed PTLD (median age 5.67 years, range 0.96-17.45). A predominant EBNA-1 variant (P-thr) was identified across the study groups. Differences were observed between the samples from the IM patients compared with the transplant samples.

CONCLUSION

The predominant EBNA-1 strain is in contrast to reports of the predominant strain in North America. The results suggest differences between the EBNA-1 strains among the study groups. Further studies will examine the relationship between EBNA-1 strains and PTLD occurrence and outcomes.

IRF4 promotes Epstein-Barr virus activation in Burkitt's lymphoma cells

Epstein-Barr virus (EBV) establishes a life-long latency in memory B cells, whereas plasma cell differentiation is linked to EBV lytic reactivation from latently infected B cells. EBV lytic replication is mediated by the two immediate-early switch proteins Zta and RTA. Both plasma cell transcription factors XBP-1 and Blimp-1 have been shown to enable the triggering of EBV lytic reactivation by activating the transcription of Zta or RTA. Here we show that interferon regulatory factor 4 (IRF4), another plasma cell transcription factor that is either not expressed or expressed at a low level in EBV-positive Burkitt's lymphoma (BL) cells, can activate the promoters of EBV Zta and RTA, but is not sufficient to elicit EBV lytic reactivation in latently infected BL cells. However, ectopic IRF4 expression can augment EBV lytic gene expression induced by anti-immunoglobulin (anti-Ig) or sodium butyrate treatment in all tested lymphoma cells, whereas IRF4 knockout in Raji cells, the only BL cell line with detectable endogenous IRF4 expression, abolishes EBV lytic gene expression induced by anti-Ig, and this is accompanied by the reduction of Blimp-1 expression, whose overexpression, in turn, can rescue EBV lytic gene expression in IRF4 knockout Raji cells. Furthermore, IRF4 knockout impairs B cell receptor (BCR) signalling activation, which is required for BCR-mediated EBV reactivation. Altogether, these results demonstrate that IRF4 facilitates EBV lytic reactivation in BL cells, which involves the regulation of Blimp-1 expression and BCR signalling pathways.

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.

Epstein–Barr Virus and Cancer

Epstein–Barr virus (EBV) contributes to about 1.5% of all cases of human cancer worldwide, and viral genes are expressed in the malignant cells. EBV also very efficiently causes the proliferation of infected human B lymphocytes. The functions of the viral proteins and small RNAs that may contribute to EBV-associated cancers are becoming increasingly clear, and a broader understanding of the sequence variation of the virus genome has helped to interpret their roles. The improved understanding of the mechanisms of these cancers means that there are great opportunities for the early diagnosis of treatable stages of EBV-associated cancers and the use of immunotherapy to target EBV-infected cells or overcome immune evasion. There is also scope for preventing disease by immunization and for developing therapeutic agents that target the EBV gene products expressed in the cancers.

Natural Variations in BRLF1 Promoter Contribute to the Elevated Reactivation Level of Epstein-Barr Virus in Endemic Areas of Nasopharyngeal Carcinoma

Background

Epstein-Barr virus (EBV) infection is a crucial risk factor for nasopharyngeal carcinoma (NPC), but the mechanism for its elevated activation level in NPC endemic areas remains unclear. This study aims to identify the EBV natural variations contributed to the different reactivation potential between NPC endemic and non-endemic areas.

Methods

1030 subjects were recruited in China, including 303 healthy individuals from two NPC non-endemic areas, 483 healthy people from three endemic areas and 244 NPC patients. Among which, saliva DNA samples from 244 participants were sequenced for the EBV immediate early (IE) genes of BRLF1 and BZLF1, their promoters were included; the rest 786 subjects were used for the validation of significant variations among three different populations. Haplotype and population structure analysis were conducted. Dual-luciferase assay was used to detect the promoter activity.

Results

A total of 246 distinct variations were detected, 29 showed significant difference in the frequencies between healthy people from NPC endemic area and non-endemic area. Population structure analysis clustered EBV strains into 9 subgroups mostly in accordance with the geographical origin of samples. Interestingly, two EBV genotypes, Rp-V1 and Rp-V2, were identified according to the linkage relationship of the variations in BRLF1 promoter (Rp). Rp-V1 has higher frequency in NPC endemic areas than in non-endemic areas (52.38% vs 18.15%, P = 2.07 × 10⁻¹⁴), and was associated with higher oral EBV DNA levels (adjusted OR = 1.64, 95% CI = 1.21–2.24, P = .002), suggesting a more powerful activation ability of Rp-V1 than that of the prototype Rp-of the EBV strain; On the contrary, Rp-V2 has higher frequency in NPC non-endemic areas than in endemic areas (18.48% vs 0.38%, P = 1.17 × 10⁻⁷), might represent a reduced activation potential of EBV. Further dual-luciferase assay showed Rp-V1 has higher promoter activity while compared with Rp-V2 (P < .0001). Notably, Rp-V1 impaired the transcription repression effect of YY1 while Rp-V2 strengthened the transcription repression effect of EBF1 on Rp. In addition, significant differences of Rta 393–407 CTL epitope which may influence the recognition of Rta by CD8+ T cells were detected between healthy people from NPC endemic area and non-endemic area.

Conclusions

This study identified natural variations in cis-acting elements (YY1 and EBF1) of EBV Rp altering Rp transcription activities, which may contribute to the elevated EBV activation level in NPC endemic areas than non-endemic areas.

Sequence Variation of Epstein-Barr Virus: Viral Types, Geography, Codon Usage, and Diseases

Epstein-Barr virus causes most cases of infectious mononucleosis and posttransplant lymphoproliferative disease. It contributes to several types of cancer, including Hodgkin's lymphoma, Burkitt's lymphoma, diffuse large B cell lymphoma, nasopharyngeal carcinoma, and gastric carcinoma. EBV genome variation is important because some of the diseases associated with EBV have very different incidences in different populations and geographic regions, and differences in the EBV genome might contribute to these diseases. Some specific EBV genome alterations that appear to be significant in EBV-associated cancers are already known, and current efforts to make an EBV vaccine and antiviral drugs should also take account of sequence differences in the proteins used as targets.

One hundred thirty-eight new Epstein-Barr virus (EBV) genome sequences have been determined. One hundred twenty-five of these and 116 from previous reports were combined to produce a multiple-sequence alignment of 241 EBV genomes, which we have used to analyze variation within the viral genome. The type 1/type 2 classification of EBV remains the major form of variation and is defined mostly by EBNA2 and EBNA3, but the type 2 single-nucleotide polymorphisms (SNPs) at the EBNA3 locus extend into the adjacent gp350 and gp42 genes, whose products mediate infection of B cells by EBV. A small insertion within the BART microRNA region of the genome was present in 21 EBV strains. EBV from saliva of U.S. patients with chronic active EBV infection aligned with the wild-type EBV genome with no evidence of WZhet rearrangements. The V3 polymorphism in the Zp promoter for BZLF1 was found to be frequent in nasopharyngeal carcinoma cases from both Hong Kong and Indonesia. Codon usage was found to differ between latent and lytic cycle EBV genes, and the main forms of variation of the EBNA1 protein have been identified.

IMPORTANCE Epstein-Barr virus causes most cases of infectious mononucleosis and posttransplant lymphoproliferative disease. It contributes to several types of cancer, including Hodgkin's lymphoma, Burkitt's lymphoma, diffuse large B cell lymphoma, nasopharyngeal carcinoma, and gastric carcinoma. EBV genome variation is important because some of the diseases associated with EBV have very different incidences in different populations and geographic regions, and differences in the EBV genome might contribute to these diseases. Some specific EBV genome alterations that appear to be significant in EBV-associated cancers are already known, and current efforts to make an EBV vaccine and antiviral drugs should also take account of sequence differences in the proteins used as targets.

Recombination rates along the entire Epstein Barr virus genome display a highly heterogeneous landscape

Epstein Barr virus (EBV) has a large DNA genome assumed to be stable, but also subject to mutational processes such as nucleotide substitution and recombination, the latter explored to a lesser extent. Moreover, differences in the extent of recombination events across herpes sub-families were recently reported. Given the relevance of recombination in viral evolution and its possible impact in pathogenesis, we aimed to fully characterize and quantify its extension in all available EBV complete genome by assessing global and local recombination rate values (⍴/bp). Our results provide the first EBV recombination map based on recombination rates assessment, both at a global and gene by gene level, where the mean value for the entire genome was 0.035 (HPDI 0.020-0.062) ⍴/bp. We quantified how this evolutionary process changes along the EBV genome, and proved it to be non-homogeneous, since regulatory regions depicted the lowest recombination rate values while repetitive regions the highest signal. Moreover, GC content rich regions seem not to be linked to high recombination rates as previously reported. At an intragenic level, four genes (EBNA3C, EBNA3B, BRRF2 and BBLF2-BBLF3) presented a recombination rate above genome average. We specifically quantified the signal strength among different recombination-initiators previously described features and concluded that those which elicited the greatest amount of changes in ⍴/bp, TGGAG and CCCAG, were two well characterized recombination inducing motifs in eukaryotic cells. Strikingly, although TGGAG was not the most frequently detected DNA motif across the EBV genome (697 hits), it still induced a significantly greater proportion of initiation events (0.025 events/hits) than other more represented motifs, p-value = 0.04; one tailed proportion test. Present results support the idea that diversity and evolution of herpesviruses are impacted by mechanisms, such as recombination, which extends beyond the usual consideration of point mutations.