Geographic Population Structure in Epstein-Barr Virus Revealed by Comparative Genomics

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.

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.

Whole genome sequence analysis of equid gammaherpesvirus -2 field isolates reveals high levels of genomic diversity and recombination

Background

Equid gammaherpesvirus 2 (EHV2) is a gammaherpesvirus with a widespread distribution in horse populations globally. Although its pathogenic significance can be unclear in most cases of infection, EHV2 infection can cause upper respiratory tract disease in foals. Co-infection of different strains of EHV2 in an individual horse is common. Small regions of the EHV2 genome have shown considerable genetic heterogeneity. This could suggest genomic recombination between different strains of EHV2, similar to the extensive recombination networks that have been demonstrated for some alphaherpesviruses. This study examined natural recombination and genome diversity of EHV2 field isolates.

Results

Whole genome sequencing analysis of 18 EHV2 isolates, along with analysis of two publicly available EHV2 genomes, revealed variation in genomes sizes (from 173.7 to 184.8 kbp), guanine plus cytosine content (from 56.7 to 57.8%) and the size of the terminal repeat regions (from 17,196 to 17,551 bp). The nucleotide sequence identity between the genomes ranged from 86.2 to 99.7%. The estimated average inter-strain nucleotide diversity between the 20 EHV2 genomes was 2.9%. Individual gene sequences showed varying levels of nucleotide diversity and ranged between 0 and 38.1%. The ratio of nonsynonymous substitutions, Ka, to synonymous substitutions, Ks, (Ka/Ks) suggests that over 50% of EHV2 genes are undergoing diversifying selection. Recombination analyses of the 20 EHV2 genome sequences using the recombination detection program (RDP4) and SplitsTree revealed evidence of viral recombination.

Conclusions

Analysis of the 18 new EHV2 genomes alongside the 2 previously sequenced genomes revealed a high degree of genetic diversity and extensive recombination networks. Herpesvirus genome diversification and virus evolution can be driven by recombination, and our findings are consistent with recombination being a key mechanism by which EHV2 genomes may vary and evolve.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08789-x.

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.

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.

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.

A comprehensive analysis of genetic diversity of EBV reveals potential high-risk subtypes associated with nasopharyngeal carcinoma in China

Epstein-Barr virus (EBV), a widespread oncovirus, is associated with multiple cancers including nasopharyngeal carcinoma (NPC), gastric cancer and diverse lymphoid malignancies. Recent studies reveal that specific EBV strains or subtypes are associated with NPC development in endemic regions. However, these NPC specific subtypes were only identified in a portion of infected individuals due possibly to the limited samples size studied or the complicated population structures of the virus. To identify additional high-risk EBV subtypes, we conducted a comprehensive genetic analysis of 22 critical viral proteins by using the largest dataset of 628 EBV genomes and 792 sequences of single target genes/proteins from GenBank. The phylogenetic, principal component and genetic structure analyses of these viral proteins were performed through worldwide populations. In addition to the general Asia-Western/Africa geographic segregation, population structure analysis showed a 'Chinese-unique' cluster (96.57% isolates from China) was highly enriched in the NPC patients, compared to the healthy individuals (89.6% vs. 44.5%, P < 0.001). The newly identified EBV subtypes, which contains four Chinese-specific NPC-associated amino acid substitutions (BALF2 V317M, BNRF1 G696R, V1222I and RPMS1 D51E), showed a robust positive association with the risk of NPC in China (Odds Ratio = 4.80, 20.00, 18.24 and 32.00 for 1, 2, 3 and 4 substitutions, respectively, P _trend <0.001). Interestingly, the coincidence of positively selected sites with NPC-associated substitutions suggests that adaptive nonsynonymous mutation on critical proteins, such as BNRF1, may interact with host immune system and contribute to the carcinogenesis of NPC. Our findings provide a comprehensive overview of EBV genetic structure for worldwide populations and offer novel clues to EBV carcinogenesis from the aspect of evolution.

When Viruses Cross Developmental Pathways

Aberrant regulation of developmental pathways plays a key role in tumorigenesis. Tumor cells differ from normal cells in their sustained proliferation, replicative immortality, resistance to cell death and growth inhibition, angiogenesis, and metastatic behavior. Often they acquire these features as a consequence of dysregulated Hedgehog, Notch, or WNT signaling pathways. Human tumor viruses affect the cancer cell hallmarks by encoding oncogenic proteins, and/or by modifying the microenvironment, as well as by conveying genomic instability to accelerate cancer development. In addition, viral immune evasion mechanisms may compromise developmental pathways to accelerate tumor growth. Viruses achieve this by influencing both coding and non-coding gene regulatory pathways. Elucidating how oncogenic viruses intersect with and modulate developmental pathways is crucial to understanding viral tumorigenesis. Many currently available antiviral therapies target viral lytic cycle replication but with low efficacy and severe side effects. A greater understanding of the cross-signaling between oncogenic viruses and developmental pathways will improve the efficacy of next-generation inhibitors and pave the way to more targeted antiviral therapies.

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.

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.

Molecular characterization of Epstein-Barr virus variants detected in the oral cavity of adolescents in Cali, Colombia

Introduction: The Epstein-Barr virus (EBV) is an ubiquitous and oncogenic virus associated with the development of diseases such as infectious mononucleosis, Burkitt’s lymphoma, nasopharyngeal carcinoma, and other neoplasms. Currently, two types are recognized: EBV-1 and EBV-2, which have genetic differences with their EBNA nuclear antigens. Likewise, due to the high degree of heterogeneity and variability found in the LMP1 protein of the virus, variants associated with pathogenesis or specific geographic regions have been described. Objective: To identify and characterize molecularly EBV variants detected in the oral cavity of 84 adolescents in Cali, Colombia. Materials and methods: Conventional PCR amplification, purification, and sequencing of the gen EBNA3C were carried out to typify the virus and the C-ter domain of the LMP1 protein to identify variants. We also conducted a phylogenetic and nucleotide variant analysis of the obtained sequences versus pathogenic or geographic variants reported in GenBank-NCBI. Results: The predominant viral subtype was EBV-1 (79%); 72.6% was grouped with the pathogenic variant Raji, derived from B lymphocytes of a patient with Burkitt›s lymphoma, 13.7% was related to a variant of Mediterranean origin, and 13.7% was not grouped with any of the reference variants. Conclusions: This is the first time that variants of LMP1-EBV have been identified in Cali, Colombia. Additional studies are necessary to characterize the unidentified variant and to determine if it is pathogenic or if it is just an isolate present in the city of Cali.

Evolution and Genetic Diversity of Primate Cytomegaloviruses

Cytomegaloviruses (CMVs) infect many mammals, including humans and non-human primates (NHPs). Human cytomegalovirus (HCMV) is an important opportunistic pathogen among immunocompromised patients and represents the most common infectious cause of birth defects. HCMV possesses a large genome and very high genetic diversity. NHP-infecting CMVs share with HCMV a similar genomic organization and coding content, as well as the course of viral infection. Recent technological advances have allowed the sequencing of several HCMV strains from clinical samples and provided insight into the diversity of NHP-infecting CMVs. The emerging picture indicates that, with the exclusion of core genes (genes that have orthologs in all herpesviruses), CMV genomes are relatively plastic and diverse in terms of gene content, both at the inter- and at the intra-species level. Such variability most likely underlies the strict species-specificity of these viruses, as well as their ability to persist lifelong and with relatively little damage to their hosts. However, core genes, despite their strong conservation, also represented a target of adaptive evolution and subtle changes in their coding sequence contributed to CMV adaptation to different hosts. Indubitably, important knowledge gaps remain, the most relevant of which concerns the role of viral genetics in HCMV-associated human disease.

Molecular Epidemiology of Epstein-Barr virus (EBV) in Patients with Hematologic Malignancies

Background:

Epstein-Barr virus (EBV) is associated with different malignant diseases, such as Hodgkin lymphoma (HL) and lymphoproliferative disorders. Patients with hematologic malignancies by variable severity could be suspected for the infection with different types of this virus. This preliminary study reported the genotyping and related viral load of Epstein-Barr virus in Iranian patients with hematologic malignancies for estimation of possible factors affecting malignancy.

Methods:

Peripheral blood mononuclear cells (PBMC) of HL (n=20), NHL (n=29), acute lymphocytic leukemia (ALL) (n=18) and chronic lymphocytic leukemia (CLL) (n=12) were obtained. After DNA extraction, a nested-PCR and a conventional-PCR targeting EBNA-2 and EBNA-3C genes were performed. A real-time PCR assay for viral load quantitation carried out. Standard curve analysis used for evaluation of amplification specificity.

Results:

Of 79 included patients, 34 (43%) were EBV positive. There were 23.5% (8/34), 38.2% (13/34), 23.5% (8/34), 14.8% (5/34) in HL, NHL, ALL and CLL groups, respectively. Also, the main genotype was genotype I (91.2%) which it follows by 8.8% (3/34) genotype II. The real-time PCR assay showed the mean viral load ± std. deviation was 2.75×10⁵ ± 1.202×10⁶ copies/μg DNA and the higher viral load was seen in NHL patients.

Conclusion:

This preliminary investigation in Iran shows that the main EBV genotype into our region probably is genotype I (91.2%) which it is similar to others. We could not find any statistically significant association between the virus infection and viral load with any specific disease and patients’ demographic data.

Human Herpesvirus Sequencing in the Genomic Era: The Growing Ranks of the Herpetic Legion

The nine human herpesviruses are some of the most ubiquitous pathogens worldwide, causing life-long latent infection in a variety of different tissues. Human herpesviruses range from mild childhood infections to known tumour viruses and 'trolls of transplantation'. Epstein-Barr virus was the first human herpesvirus to have its whole genome sequenced; GenBank now includes thousands of herpesvirus genomes. This review will cover some of the recent advances in our understanding of herpesvirus diversity and disease that have come about as a result of new sequencing technologies, such as target enrichment and long-read sequencing. It will also look at the problem of resolving mixed-genotype infections, whether with short or long-read sequencing methods; and conclude with some thoughts on the future of the field as herpesvirus population genomics becomes a reality.

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.

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.

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.

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.

A cancer-associated Epstein-Barr virus BZLF1 promoter variant enhances lytic infection

Latent Epstein-Barr virus (EBV) infection contributes to both B-cell and epithelial-cell malignancies. However, whether lytic EBV infection also contributes to tumors is unclear, although the association between malaria infection and Burkitt lymphomas (BLs) may involve excessive lytic EBV replication. A particular variant of the viral promoter (Zp) that controls lytic EBV reactivation is over-represented, relative to its frequency in non-malignant tissue, in EBV-positive nasopharyngeal carcinomas and AIDS-related lymphomas. To date, no functional differences between the prototype Zp (Zp-P) and the cancer-associated variant (Zp-V3) have been identified. Here we show that a single nucleotide difference between the Zp-V3 and Zp-P promoters creates a binding site for the cellular transcription factor, NFATc1, in the Zp-V3 (but not Zp-P) variant, and greatly enhances Zp activity and lytic viral reactivation in response to NFATc1-inducing stimuli such as B-cell receptor activation and ionomycin. Furthermore, we demonstrate that restoring this NFATc1-motif to the Zp-P variant in the context of the intact EBV B95.8 strain genome greatly enhances lytic viral reactivation in response to the NFATc1-activating agent, ionomycin, and this effect is blocked by the NFAT inhibitory agent, cyclosporine, as well as NFATc1 siRNA. We also show that the Zp-V3 variant is over-represented in EBV-positive BLs and gastric cancers, and in EBV-transformed B-cell lines derived from EBV-infected breast milk of Kenyan mothers that had malaria during pregnancy. These results demonstrate that the Zp-V3 enhances EBV lytic reactivation to physiologically-relevant stimuli, and suggest that increased lytic infection may contribute to the increased prevalence of this variant in EBV-associated malignancies.

Early Epstein-Barr Virus Genomic Diversity and Convergence toward the B95.8 Genome in Primary Infection

Over 90% of the world's population is persistently infected with Epstein-Barr virus. While EBV does not cause disease in most individuals, it is the common cause of acute infectious mononucleosis (AIM) and has been associated with several cancers and autoimmune diseases, highlighting a need for a preventive vaccine. At present, very few primary, circulating EBV genomes have been sequenced directly from infected individuals. While low levels of diversity and low viral evolution rates have been predicted for double-stranded DNA (dsDNA) viruses, recent studies have demonstrated appreciable diversity in common dsDNA pathogens (e.g., cytomegalovirus). Here, we report 40 full-length EBV genome sequences obtained from matched oral wash and B cell fractions from a cohort of 10 AIM patients. Both intra- and interpatient diversity were observed across the length of the entire viral genome. Diversity was most pronounced in viral genes required for establishing latent infection and persistence, with appreciable levels of diversity also detected in structural genes, including envelope glycoproteins. Interestingly, intrapatient diversity declined significantly over time (P < 0.01), and this was particularly evident on comparison of viral genomes sequenced from B cell fractions in early primary infection and convalescence (P < 0.001). B cell-associated viral genomes were observed to converge, becoming nearly identical to the B95.8 reference genome over time (Spearman rank-order correlation test; r = -0.5589, P = 0.0264). The reduction in diversity was most marked in the EBV latency genes. In summary, our data suggest independent convergence of diverse viral genome sequences toward a reference-like strain within a relatively short period following primary EBV infection.IMPORTANCE Identification of viral proteins with low variability and high immunogenicity is important for the development of a protective vaccine. Knowledge of genome diversity within circulating viral populations is a key step in this process, as is the expansion of intrahost genomic variation during infection. We report full-length EBV genomes sequenced from the blood and oral wash of 10 individuals early in primary infection and during convalescence. Our data demonstrate considerable diversity within the pool of circulating EBV strains, as well as within individual patients. Overall viral diversity decreased from early to persistent infection, particularly in latently infected B cells, which serve as the viral reservoir. Reduction in B cell-associated viral genome diversity coincided with a convergence toward a reference-like EBV genotype. Greater convergence positively correlated with time after infection, suggesting that the reference-like genome is the result of selection.

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.

A novel recombinant variant of latent membrane protein 1 from Epstein Barr virus in Argentina denotes phylogeographical association

Epstein Barr virus (EBV) infection in Argentina occurs at an early age and occasionally develops infectious mononucleosis (IM). EBV is also related with lymphomas. LMP1, the viral oncoprotein is polymorphic and is used to define viral variants.