Genome sequencing analysis identifies Epstein-Barr virus subtypes associated with high risk of nasopharyngeal carcinoma - PubMed (original) (raw)

. 2019 Jul;51(7):1131-1136.

doi: 10.1038/s41588-019-0436-5. Epub 2019 Jun 17.

Youyuan Yao 1 3, Hui Chen 4, Shanshan Zhang 1, Su-Mei Cao 1, Zhe Zhang 5, Bing Luo 6, Zhiwei Liu 7, Zilin Li 2, Tong Xiang 1, Guiping He 1, Qi-Sheng Feng 1, Li-Zhen Chen 1, Xiang Guo 1 8, Wei-Hua Jia 1, Ming-Yuan Chen 1, Xiao Zhang 1, Shang-Hang Xie 1, Roujun Peng 1, Ellen T Chang 9 10, Vincent Pedergnana 4, Lin Feng 1, Jin-Xin Bei 1, Rui-Hua Xu 1, Mu-Sheng Zeng 1, Weimin Ye 7, Hans-Olov Adami 7 11, Xihong Lin 2, Weiwei Zhai 12 13 14, Yi-Xin Zeng 15, Jianjun Liu 16 17

Affiliations

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

Miao Xu et al. Nat Genet. 2019 Jul.

Abstract

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 identify two non-synonymous EBV variants within BALF2 that are strongly associated with the risk of NPC (odds ratio (OR) = 8.69, P = 9.69 × 10-25 for SNP 162476_C; OR = 6.14, P = 2.40 × 10-32 for SNP 163364_T). The cumulative effects of these variants contribute to 83% of the overall risk of NPC in southern China. Phylogenetic analysis of the risk variants reveals a unique origin in Asia, followed by clonal expansion in NPC-endemic regions. Our results provide novel insights into the NPC endemic in southern China and also enable the identification of high-risk individuals for NPC prevention.

PubMed Disclaimer

Conflict of interest statement

Competing interests

The authors have no conflicts of interest to disclose.

Figures

Figure 1

Figure 1. Principal component and phylogenetic analyses of EBV genomes.

(a) Principal component analysis of 270 EBV isolates sequenced in the current study and 97 published isolates. The first two principal-component scores (PC1 and PC2) are plotted. Explaining 26.9% of the total genomic variance, PC1 discriminates between East Asian and Western/African strains; PC2 explains 15.3% of the total variance. (b) Phylogeny of 230 EBV single strains sequenced in the current study and 97 published strains. Macacine herpesvirus 4 genome sequence (NC_006146) was used as the outgroup to root the tree. Type 1 and Type 2 EBV lineages are indicated. The red dot on the phylogeny indicates the lineage of the NPC-dominant EBV strains, where 22 of 37 strains from healthy controls in NPC-endemic regions in southern China were located. Dashed lines in (a) and (b) indicate the separation between East Asian and Western/African strains. (c) Geographical origins and phenotypes of samples from which EBV strains were sequenced are shown with colors as indicated. (d) The normalized values of PC1 and PC2 scores are shown from blue to red. (e) The genotypes of SNPs 162215C>A, 162476T>C, and 163364C>T in each isolate.

Figure 2

Figure 2. Genome-wide association analysis of EBV variants in 156 NPC cases and 47 controls.

(a) Manhattan plot of genome-wide P values from the association analysis using a generalized-linear mixed model. The −log10-transformed P values (y axis) of 1545 variants in 156 NPC cases and 47 controls are presented according to their positions in the EBV genome. The minimum P value (SNP 162507C>T) is 9.17×10−5. The red line is the suggestive genome-wide significance P value threshold of 4.07×10−4. The three SNPs 162507C>T, 162852G>T and 162215C>A reaching genome-wide significance are labeled as green. (b) The regional plot of the posterior probabilities of association. The EBV genome was partitioned into overlapping 20-variant bins with 10-variant overlap between adjacent bins. The sum of the posterior probabilities for variants was assigned to each region. The one region from position 160971 to 163629 with strong evidence (> 0.85) for association with NPC risk is shown in green. (c) Schematic of EBV genes. Repetitive regions in EBV genomes are masked by light blue.

Comment in

References

    1. Epstein MA, Achong BG & Barr YM Virus Particles in Cultured Lymphoblasts from Burkitt’s Lymphoma. Lancet 1, 702–3 (1964). -PubMed
    1. Epstein A Why and How Epstein-Barr Virus Was Discovered 50 Years Ago. Curr Top Microbiol Immunol 390, 3–15 (2015). -PubMed
    1. Kieff ED & Rickinson AB Epstein-Barr Virus and Its Replication. in Fields’ virology Vol. 68A (eds. Knipe DM & Howley PM) 2603–2654 (Lippincott Williams & Wilkins, Wolters Kluwer, Philadelphia, 2007).
    1. Zhang LF et al. Incidence trend of nasopharyngeal carcinoma from 1987 to 2011 in Sihui County, Guangdong Province, South China: an age-period-cohort analysis. Chin J Cancer 34, 350–7 (2015). -PMC -PubMed
    1. Bei JX et al. A genome-wide association study of nasopharyngeal carcinoma identifies three new susceptibility loci. Nat Genet 42, 599–603 (2010). -PubMed

Methods-only references

    1. Li H & Durbin R Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–60 (2009). -PMC -PubMed
    1. Li H et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–9 (2009). -PMC -PubMed
    1. DePristo MA et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet 43, 491–8 (2011). -PMC -PubMed
    1. Cingolani P et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin) 6, 80–92 (2012). -PMC -PubMed
    1. Raab-Traub N & Flynn K The structure of the termini of the Epstein-Barr virus as a marker of clonal cellular proliferation. Cell 47, 883–9 (1986). -PubMed

Publication types

MeSH terms

Substances

Grants and funding

LinkOut - more resources