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.

Children diagnosed with endemic Burkitt lymphoma (eBL) are deficient in interferon-gamma (IFN-gamma) responses to Epstein-Barr Nuclear Antigen1 (EBNA1), the viral protein that defines the latency I pattern in this B cell tumor. However, the contributions of immune-regulatory cytokines and phenotypes of the EBNA1-specific T cells have not been characterized for eBL. Using a bespoke flow cytometry assay we measured intracellular IFN-gamma, IL-10, IL-17A expression and phenotyped CD4(+) and CD8(+) T cell effector memory subsets specific to EBNA1 for eBL patients compared to two groups of healthy children with divergent malaria exposures. In response to EBNA1 and a malaria antigen (PfSEA-1A), the three study groups exhibited strikingly different cytokine expression and T cell memory profiles. EBNA1-specific IFN-gamma-producing CD4(+) T cell response rates were lowest in eBL (40%) compared to children with high malaria (84%) and low malaria (66%) exposures (p < 0.0001 and p = 0.0004, respectively). However, eBL patients did not differ in CD8(+) T cell response rates or the magnitude of IFN-gamma expression. In contrast, eBL children were more likely to have EBNA1-specific CD4(+) T cells expressing IL-10, and less likely to have polyfunctional IFN-gamma(+)IL-10(+) CD4(+) T cells (p = 0.02). They were also more likely to have IFN-gamma(+)IL-17A(+), IFN-gamma(+) and IL-17A(+) CD8(+) T cell subsets compared to healthy children. Cytokine-producing T cell subsets were predominantly CD45RA(+)CCR7(+) TNAIVE-LIKE cells, yet PD-1, a marker of persistent activation/exhaustion, was more highly expressed by the central memory (TCM) and effector memory (TEM) T cell subsets. In summary, our study suggests that IL-10 mediated immune regulation and depletion of IFN-gamma(+) EBNA1-specific CD4(+) T cells are complementary mechanisms that contribute to impaired T cell cytotoxicity in eBL pathogenesis.