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

Evolutionary analysis of LMP-1 genetic diversity in EBV-associated nasopharyngeal carcinoma: Bioinformatic insights into oncogenic potential

EBV latent membrane protein 1 (LMP-1) is an important oncogene involved in the induction and maintenance of EBV infection and the activation of several cell survival and proliferative pathways. The genetic diversity of LMP-1 has an important role in immunogenicity and tumorigenicity allowing escape from host cell immunity and more metastatic potential of LMP-1 variants. This study explored the evolutionary of LMP-1 in EBV-infected patients at an advanced stage of nasopharyngeal carcinoma (NPC). Detection of genetic variability in LMP-1 genes was carried out using Sanger sequencing. Bioinformatic analysis was conducted for translation and nucleotide alignment. Phylogenetic analysis was used to construct a Bayesian tree for a deeper understanding of the genetic relationships, evolutionary connections, and variations between sequences. Genetic characterization of LMP-1 in NPC patients revealed the detection of polymorphism in LMP-1 Sequences. Motifs were identified within three critical LMP-1 domains, such as PQQAT within CTAR1 and YYD within CTAR2. The presence of the JACK3 region at specific sites within CTAR3, as well as repeat regions at positions (122-132) and (133-143) within CTAR3, was also annotated. Additionally, several mutations were detected including 30 and 69 bp deletions, 33 bp repeats, and 15 bp insertion. Although LMP-1 strains appear to be genetically diverse, they are closely related to 3 reference strains: prototype B95.8, Med- 30 bp deletion, and Med + 30 bp deletion. In our study, one of the strains harboring the 30 bp deletion had both bone and bone marrow metastasis which could be attributed to the fact that LMP-1 is involved in tumor metastasis, evasion and migration of NPC cells. This study provided valuable insights into genetic variability in LMP-1 sequences of EBV in NPC patients. Further functional studies would provide a more comprehensive understanding of the molecular characteristics, epidemiology, and clinical implications of LMP-1 polymorphisms in EBV-related malignancies.

The impact of HLA polymorphism on herpesvirus infection and disease

Human Leukocyte Antigens (HLA) are cell surface molecules, central in coordinating innate and adaptive immune responses, that are targets of strong diversifying natural selection by pathogens. Of these pathogens, human herpesviruses have a uniquely ancient relationship with our species, where coevolution likely has reciprocating impact on HLA and viral genomic diversity. Consistent with this notion, genetic variation at multiple HLA loci is strongly associated with modulating immunity to herpesvirus infection. Here, we synthesize published genetic associations of HLA with herpesvirus infection and disease, both from case/control and genome-wide association studies. We analyze genetic associations across the eight human herpesviruses and identify HLA alleles that are associated with diverse herpesvirus-related phenotypes. We find that whereas most HLA genetic associations are virus- or disease-specific, HLA-A*01 and HLA-A*02 allotypes may be more generally associated with immune susceptibility and control, respectively, across multiple herpesviruses. Connecting genetic association data with functional corroboration, we discuss mechanisms by which diverse HLA and cognate receptor allotypes direct variable immune responses during herpesvirus infection and pathogenesis. Together, this review examines the complexity of HLA-herpesvirus interactions driven by differential T cell and Natural Killer cell immune responses.

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.

Genomic landscape of Epstein–Barr virus in familial nasopharyngeal carcinoma

To better understand the genomic characteristics of Epstein–Barr virus (EBV) in familial nasopharyngeal carcinoma (NPC), we sequenced the EBV genomes by whole-genome capture in 38 unrelated patients with NPC family history in first-degree relatives and 47 healthy controls, including 13 with family history and 34 without. Compared with type 1 reference genome, mutation hotspots were observed in the latent gene regions of EBV in familial NPC cases. Population structure analysis showed that one cluster has a higher frequency in familial cases than in controls (OR=5.33, 95 % CI 2.50–11.33, P=1.42×10−5), and similar population structure composition was observed among familial and sporadic NPC cases in high-endemic areas. By genome-wide association analysis, four variants were found to be significantly associated with familial NPC. Consistent results were observed in the meta-analysis integrating two published case-control EBV sequencing studies in NPC high-endemic areas. High-risk haplotypes of EBV composed of 34 variants were associated with familial NPC risk (OR=13.85, 95 % CI 4.13–46.44, P=2.06×10−5), and higher frequency was observed in healthy blood-relative controls with NPC family history (9/13, 69.23 %) than those without family history (16/34, 47.06%). This study suggested the potential contribution of EBV high-risk subtypes to familial aggregation of NPC.

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 (EBV) deletions as biomarkers of response to treatment of chronic active EBV

Chronic active Epstein-Barr virus (CAEBV) disease is a rare condition characterised by persistent EBV infection in previously healthy individuals. Defective EBV genomes were found in East Asian patients with CAEBV. In the present study, we sequenced 14 blood EBV samples from three UK patients with CAEBV, comparing the results with saliva CAEBV samples and other conditions. We observed EBV deletions in blood, some of which may disrupt viral replication, but not saliva in CAEBV. Deletions were lost overtime after successful treatment. These findings are compatible with CAEBV being associated with the evolution and persistence of EBV+ haematological clones that are lost on successful treatment.

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.

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.

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.

Analysis of latent T-cell epitopes in Epstein-Barr virus isolated from extranodal nasal-type natural killer/T-cell lymphoma in Taiwanese population

Extranodal nasal-type natural killer (NK)/T-cell lymphoma (NKTCL) is an aggressive lymphoma that is prevalent among East Asian and South American populations. Although Epstein-Barr virus (EBV) is commonly detected in NKTCL, there are limited studies that have analyzed the EBV genomic variations in NKTCL. In this study, 8 EBV latent genes were analyzed using targeted gene sequencing in 23 formalin-fixed paraffin-embedded tissues derived from 18 patients with NKTCL. Five cases with paired samples were comparatively analyzed. The consistency of EBV sequencing data between tissue samples was high (96.3%-98.7%), whereas that of variant calling among the tissue samples and plasma samples (74.3%-79.2%) was low. The highest densities of non-synonymous variants were detected in the EBNA3B gene. Among the 74 known T-cell epitopes, 363 non-synonymous variants were identified in 32 (43.2%) epitopes. Additionally, the AVFDRKSDAK (A1S/P and V2F/M/L) and YHLIVDTDSL (I4L and L10R/V/G/H) epitopes were associated with 5 patterns of amino acid changes in EBNA3B and EBNA-2, respectively. The frequency of variation in the human leukocyte antigen (HLA)-restricted epitopes with corresponding HLA types common among Taiwanese population was significantly low (P = 0.011), whereas that in anchor residues was significantly high (P = 0.012). In conclusion, this study demonstrated the genomic diversity of EBV in NKTCL and its correlation with the HLA-restricted epitope variations in Taiwanese population. The findings of this study provide useful insights for the development of novel therapeutic strategies for NKTCL.

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.

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.

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.