Viral gene sequences from an enlarged set of about 200 Epstein-Barr virus (EBV) strains, including many primary isolates, have been used to investigate variation in key viral genetic regions, particularly LMP1, Zp, gp350, EBNA1, and the BART microRNA (miRNA) cluster 2. Determination of type 1 and type 2 EBV in saliva samples from people from a wide range of geographic and ethnic backgrounds demonstrates a small percentage of healthy white Caucasian British people carrying predominantly type 2 EBV. Linkage of Zp and gp350 variants to type 2 EBV is likely to be due to their genes being adjacent to the EBNA3 locus, which is one of the major determinants of the type 1/type 2 distinction. A novel classification of EBNA1 DNA binding domains, named QCIGP, results from phylogeny analysis of their protein sequences but is not linked to the type 1/type 2 classification. The BART cluster 2 miRNA region is classified into three major variants through single-nucleotide polymorphisms (SNPs) in the primary miRNA outside the mature miRNA sequences. These SNPs can result in altered levels of expression of some miRNAs from the BART variant frequently present in Chinese and Indonesian nasopharyngeal carcinoma (NPC) samples. The EBV genetic variants identified here provide a basis for future, more directed analysis of association of specific EBV variations with EBV biology and EBV-associated diseases.IMPORTANCE Incidence of diseases associated with EBV varies greatly in different parts of the world. Thus, relationships between EBV genome sequence variation and health, disease, geography, and ethnicity of the host may be important for understanding the role of EBV in diseases and for development of an effective EBV vaccine. This paper provides the most comprehensive analysis so far of variation in specific EBV genes relevant to these diseases and proposed EBV vaccines. By focusing on variation in LMP1, Zp, gp350, EBNA1, and the BART miRNA cluster 2, new relationships with the known type 1/type 2 strains are demonstrated, and a novel classification of EBNA1 and the BART miRNAs is proposed.

Undifferentiated nasopharyngeal carcinoma (NPC) is 100% associated with Epstein‐Barr virus (EBV) as oncogenic driver. NPC is often diagnosed late due to initial vague complaints and obscured location. Prior studies suggest that measurement of EBV DNA load and RNA transcripts in nasopharyngeal (NP) brushings is useful for minimally invasive NPC diagnosis. However, whether these EBV markers relate to local virus replication or reflect tumor origin remains to be demonstrated. To resolve this, we analysed EBV‐DNA characteristics and quantified latent and lytic viral RNA transcripts in NP brushings and matching frozen NP‐biopsy specimens from patients suspected of having NPC.We observed non‐fragmented and Cp‐promotor methylated EBV‐DNA in both NP brushings and biopsies suggestive of tumor origin. Using quantitative RT‐PCR we determined expression levels of 7 critical latent (EBER1, Qp‐EBNA1, EBNA2, BART, LMP1, LMP2, BARF1) and 5 lytic (Zta, Rta, TK, PK and VCA‐p18) RNA transcripts. Although latent and early lytic RNA transcripts were frequently detected in conjunction with high EBV viral load, in both brushings and biopsies the latent transcripts prevailed and reflected a typical NPC‐associated latency‐II transcription profile without EBNA2. Late lytic RNA transcripts were rare and detected at low levels mainly in NP brushings, suggestive of abortive viral reactivation rather than complete virus replication. EBV‐IgA serology (EBNA1, VCA, Zta) did not correlate to the level of viral reactivation in situ.Overall, viral RNA profiling, DNA fragmentation and methylation analysis in NP brushings and parallel biopsies indicate that NP brush sampling provides a true and robust indicator of NPC tumor presence.

Cell‐free microRNA (miRNA) in biofluids released by tumors in either protein or vesicle‐bound form, represent promising minimally‐invasive cancer biomarkers. However, a highly abundant non‐tumor background in human plasma and serum complicates the discovery and detection of tumor‐selective circulating miRNAs. We performed small RNA sequencing on serum and plasma RNA from Nasopharyngeal Carcinoma (NPC) patients. Collectively, Epstein Barr virus‐encoded miRNAs, more so than endogenous miRNAs, signify presence of NPC. However, RNAseq‐based EBV miRNA profiles differ between NPC patients, suggesting inter‐tumor heterogeneity or divergent secretory characteristics. We determined with sensitive qRT‐PCR assays that EBV miRNAs BART7‐3p, BART9‐3p and BART13‐3p are actively secreted by C666.1 NPC cells bound to extracellular vesicles (EVs) and soluble ribonucleoprotein complexes. Importantly, these miRNAs are expressed in all primary NPC tumor biopsies and readily detected in nasopharyngeal brushings from both early and late‐stage NPC patients. Increased levels of BART7‐3p, BART9‐3p and particularly BART13‐3p, distinguish NPC patient sera from healthy controls. Receiver operating characteristic curve analysis using sera from endemic NPC patients, other head and neck cancers and individuals with asymptomatic EBV‐infections reveals a superior diagnostic performance of EBV miRNAs over anti‐EBNA1 IgA serology and EBV‐DNA load (AUC 0.87–0.96 vs 0.86 and 0.66 respectively). The high specificity of circulating EBV‐BART13‐3p (97%) for NPC detection is in agreement with active secretion from NPC tumor cells. We conclude EV‐bound BART13‐3p in circulation is a promising, NPC‐selective, biomarker that should be considered as part of a screening strategy to identify NPC in endemic regions.

BackgroundEpstein-Barr virus (EBV) commonly infects the general population and has been associated with nasopharyngeal carcinoma (NPC), which has a high incidence in certain regions. This study aimed to address how EBV variations contribute to the risk of NPC.MethodsUsing logistic regression analysis and based on the sequence variations at EBV-encoded RPMS1, a multi-stage association study was conducted to identify EBV variations associated with NPC risk. A protein degradation assay was performed to characterize the functional relevance of the RPMS1 variations.ResultsBased on EBV-encoded RPMS1 variations, a single nucleotide polymorphism (SNP) in the EBV genome (locus 155391: G>A, named G155391A) was associated with NPC in 157 cases and 319 healthy controls from an NPC endemic region in South China [P < 0.001, odds ratio (OR) = 4.47, 95% confidence interval (CI) 2.71–7.37]. The results were further validated in three independent cohorts from the NPC endemic region (P < 0.001, OR = 5.20, 95% CI 3.18–8.50 in 168 cases vs. 241 controls, and P < 0.001, OR = 5.27, 95% CI 4.06–6.85 in 726 cases vs. 880 controls) and a non-endemic region (P < 0.001, OR = 7.52, 95% CI 3.69–15.32 in 58 cases vs. 612 controls). The combined analysis in 1109 cases and 2052 controls revealed that the SNP G155391A was strongly associated with NPC (Pcombined < 0.001, OR = 5.27, 95% CI 4.31–6.44). Moreover, the frequency of the SNP G155391A was associated with NPC incidence but was not associated with the incidences of other EBV-related malignancies. Furthermore, the protein degradation assay showed that this SNP decreased the degradation of the oncogenic RPMS1 protein.ConclusionsOur study identified an EBV variation specifically and significantly associated with a high risk of NPC. These findings provide insights into the pathogenesis of NPC and strategies for prevention.