Abstract
Using EpiLymph case-control data, we found that chronic lymphocytic leukemia patients were more likely to have abnormal reactive serological patterns to Epstein Barr virus than controls. Here, we aimed to assess whether this association is modified by genetic variants. We examined 1,305 Single Nucleotide Polymorphisms from 300 selected genes related to various pathways in 240 cases and 513 controls from five European centers. In a recessive model, patients positive to aberrant antibody pattern and homozygous for rare genotypes in rs8113877T>G or rs17576A>G of the MMP9 gene were at highest risk of chronic lymphocytic leukemia. In a dominant model, TP73 showed the highest risk in patients positive to aberrant antibody pattern and homozygous for the wild-type genotype in rs1885859G>C or rs3765701A>T. All interactions were additive and no main effect was observed. The strong interactions observed may be indicative of a specific pathway in cancer genesis. Confirmation of these results is warranted.Introduction
Chronic lymphocytic leukemia (CLL), one of the commonest forms of lymphoid malignancies, is characterized by an accumulation of abnormal B-cells in the peripheral blood. Chronic lymphocytic leukemia patients were three times more likely to have aberrant anti-Epstein Barr virus antibody (ab_EBV) patterns than controls.1 Recent prospective data supported this association.2 Monoclonal B-CLL tumors are generally EBV negative and frequently arise from the subgroup of CD5+ B cells, which are resistant to EBV infection. Chronic inflammation or antigen triggering may lead to loss of homeostatic B-cell control and outgrowth of a malignant B-cell clone, particularly in the elderly. Epstein Barr virus, which is normally well controlled by host immune responses but can be sub-clinically reactivated by a variety of stress signals, might trigger B-cell activation.3 Such Epstein Barr virus reactivation is reflected by aberrant Epstein Barr virus seroreactivity. This process might become functional in the background of “minor” genetic defects that increase the chance of erroneous events over time. The increased chronic lymphocytic leukemia risk in individuals with antibody signature of Epstein Barr virus reactivation, as reflected in an increased antibody diversity1,2 might reflect such an indirect role of Epstein Barr virus in chronic lymphocytic leukemia pathogenesis.
To further analyze the role of genetic susceptibility markers, we selected 1,305 single nucleotide polymorphisms (SNPs) from 300 genes associated with inflammation, apoptosis, immunoregulation, cell-cycle, metabolic and DNA-repair pathways, and examined their interaction with abnormal reactive serological patterns to Epstein Barr virus among 240 chronic lymphocytic leukemia patients and 513 controls from EpiLymph.
Design and Methods
Study population
The design of the EpiLymph study has been described in detail elsewhere.1 Briefly, the Czech Republic, France, Ireland, Italy, Germany and Spain took part from 1998 to 2003, using a comparable questionnaire and study protocol. Overall, the study includes 2,362 incident lymphoma cases and 2,465 controls. Chronic lymphocytic leukemia and small lymphocytic lymphoma were grouped and classified according to the World Health Organization (WHO) classification.4 Patients with organ transplantation or HIV infection were excluded from analysis. The large majority of chronic lymphocytic leukemia cases were untreated at the time of blood collection. The study was approved by the Research ethics committee and all patients gave their consent to participate in the study.
Laboratory methods
Detailed information on Epstein Barr virus serology has been reported elsewhere.1 Genomic DNA was quantified using PicoGreen and diluted to a final concentration of 50 ng/AL. Genotyping was performed using a customized Illumina Golden Gate genotyping assay (Illumina Inc., San Diego, USA) and all SNP genotype data were manually clustered using BeadStudio software (Illumina).
Candidate gene choice and SNP selection
In selecting candidate genes, we aimed at SNPs located in genes coding for innate immunity, specifically in NFKB and TLR signaling, other immunoregulatory and apoptosis pathways, and cell cycle regulation. Also genes predisposing for diseases associated with an elevated risk of lymphoma and susceptibility to infectious agents were considered. A list of all investigated polymorphisms along with their MAP and Hardy-Weinberg Equilibrium (HWE) P values can be found in the Online Supplementary Appendix 1.
SNP genotyping and exclusions
Overall, 241 chronic lymphocytic leukemia cases and 518 controls from EpiLymph were genotyped for this study. Following the exclusion of 3 samples that had less than 95% individual call rate and 3 subjects due to discrepancy between sex call from X chromosome genotype data and self-reported data, 240 chronic lymphocytic leukemia cases (including 6 small lymphocytic lymphoma) and 513 controls were included in the analysis. There was no difference between the proportion of cases and controls excluded due to low call rate (0.04% vs. 1.00%, P value=0.43, for cases and controls, respectively).
Of the 1,419 SNPs that were selected for genotyping, 114 were excluded (one monomorphic SNP, 27 SNPs from X chromosome, 82 SNPs with call rate less than 95% and 4 SNPs with HWE less than 1×10 among controls). For all control patients, minor allele frequency was over 1%.
Statistical analysis
Using STATA10 (StataCorp, USA), genotypic distributions of each SNP were assessed for Hardy-Weinberg Equilibrium among controls using exact χ goodness of fit test. Logistic regression models adjusted for age (tertiles based on the control age distribution), sex and country were used to estimate the odds ratio for chronic lymphocytic leukemia by each SNP, haplotypes, ab_Epstein Barr virus status and for the main effects of each SNP with detected interaction. All statistical tests were two-sided.
Dominant and recessive inheritance models were examined. Statistical interactions using multiplicative and additive scales were considered. Two tests of significance, HA0 and HM0, were used to assess departure towards multiplicative or additive interactions, respectively.5 Data on interactions were presented in a two-by-four table format. Based on these two derived P values, interactions were defined as additive, multiplicative, both or neither. Interactions were considered as genuine if HA0 and/or HM0 were rejected at the 5% level in Spain and Germany separately (the 2 countries with the most data). Matched analysis on age (−/+ 2 years), sex and center was also performed as a sensitivity analysis. The haplotype block structure was determined with Haploview 4.1 software. Haplotypes and haplotype pairs were reconstructed using the haplo.stats package for R statistical software (http://www.r-project.org/).
Results and Discussion
The characteristics of the 240 chronic lymphocytic leukemia patients and 513 controls are shown in Table 1. A higher proportion of cases than controls were men (P= 0.10) and the mean age was 63.8 years (standard deviation (SD) 9.5) in cases and 63.1 years (SD 11.7 years) in controls. The proportion of individuals with aberrant reaction to antibodies against Epstein Barr virus was 20% in controls and 39% in cases.
Four SNPs were found to interact with ab_EBV in chronic lymphocytic leukemia risk (Table 2). All interactions were additive. In a recessive model, two MMP9 SNPs (rs8113877 and rs17576) were found to modify the effect of ab_EBV on chronic lymphocytic leukemia (P values=0.002 and 0.004 for rs8113877 and rs17576, respectively). Overall, ab_EBV positive individuals were more likely to develop chronic lymphocytic leukemia compared to those who were negative, independently of their genotype. Compared to ab_EBV positive carriers of TT/TG and AA/AG for the MMP9 SNPs rs8113877 and rs17576, respectively, individuals carrying the GG genotype and being positive to ab_EBV were more likely to have been diagnosed with chronic lymphocytic leukemia. By contrast, among ab_EBV negative patients, the presence of both rare alleles (GG) decreased significantly the odds ratio for chronic lymphocytic leukemia. In a dominant model, two TP73 SNPs (rs3765701 and rs1885859) modified the association between ab_EBV and chronic lymphocytic leukemia (P values=0.005 and 0.003 for rs3765701 and rs1885859, respectively). Odds ratio of chronic lymphocytic leukemia were significantly higher among all groups compared to ab_EBV negative carriers of wild-type homozygous genotypes. Antagonist effects were observed among ab_EBV positive patients with carriers of the wild-type homozygous genotypes having higher significant odds ratio of chronic lymphocytic leukemia than those with other genotypes. Restricting the analyses to chronic lymphocytic leukemia patients with lower Rai stage did not materially change the results suggesting that clinical stage does not confound the observed association (Online Supplementary Appendix 2). Matched analysis decreased the power of the study but did not substantially alter the results (Online Supplementary Appendix 3). A list of SNPs with HA0 and/or HM0 rejected at the 10% level in Spain and Germany separately are shown in the Online Supplementary Appendix 4.
The same data as Table 2 were displayed based on multiple contingency tables (Figure 1). Here, odds ratios of chronic lymphocytic leukemia in ab_EBV positive compared to ab_EBV negative patients are shown according to the genotype of the four SNPs. Overall, the odds ratio of chronic lymphocytic leukemia in all ab_EBV positive patients compared to all ab_EBV negative patients was 2.76 (95% CI=1.91 to 3.98). The odds ratios of chronic lymphocytic leukemia associated with ab_EBV positivity were significantly higher among carriers of the GG homozygous genotypes for the 2 MMP9 SNPs than in those without this genetic characteristic. In relation to the two TP73 polymorphisms and based on a dominant model, odds ratio of chronic lymphocytic leukemia associated with ab_EBV was higher among carriers of the homozygous genotypes AA and GG compared to those with the AT/TT or GC/CC for rs3765701 and rs1885859.
Haplotype analyses and LD plots of the genotyped SNPs of the MMP9 and TP73 genes are shown in the Online Supplementary Appendices 5, 6 and 7, respectively. All five MMP9 gene variants were in strong linkage disequilibrium. Overall, the haplotype pair analysis for MMP9 supported the results observed in single SNP analysis. In particular, a significant greater odds ratio of chronic lymphocytic leukemia in relation to ab_EBV was found with the haplotype pair Hap2-Hap3 where each haplotype includes two copies of the rare allele G (Wald test P=0.00005). Using Hapmap (Phase II, NCBI B36, dbSNP b126) we found that the SNPs rs1885859 and rs3765701 of the TP73 were in strong linkage disequilibrium (D′=0.93; r²=0.453) and decided to include the borderline results for SNP rs3765701 with the main results.
A role for Epstein Barr virus activation in the etiology and progression of chronic lymphocytic leukemia is increasingly discussed.6,7 Epstein Barr virus control requires a constant lifelong state of alert of the immune system, involving approximately 1% of all T cells, which may decrease with increasing age.8 A raised EBV-reactivation with age may be reflected in an increased incidence of EBV-driven lymphoproliferative diseases and B-cell lymphomas in the elderly.9 Various stress factors can trigger Epstein Barr virus reactivation,8 evoking immune responses that are reflected in a broadened spectrum of anti-EBV antibodies. Transcriptional activity of the Epstein Barr virus latent membrane protein 1 (LMP1) with known transforming potential in vitro was recently detected in a significant proportion of chronic lymphocytic leukemia patients.7 Atypical Epstein Barr virus serological patterns1,2 appear to be associated with chronic lymphocytic leukemia and its clinical grade independently of age (Online Supplementary Appendix 8). Monoclonal B-CLL are generally EBV negative and the role of Epstein Barr virus, if any, would probably be indirect via B-cell triggering and facilitation of malignant growth.1 By considering that a major effect in the carcinogenic process was determined by presence of ab_EBV, genetic factors were evaluated as risk modifiers.
MMP9 belongs to a family of zinc-dependent endopeptidases which degrade type IV collagen, a major component of the basement membrane which is involved in tumor growth, apoptosis, angiogenesis, invasion and metastasis.10 The SNP rs17576 has been previously associated with chronic lymphocytic leukemia in a large population-based case-control study11 using recessive models. MMP9 is activated and over-expressed in B-cell chronic lymphocytic leukemia.12 Importantly, MMP9 is induced by EBV-encoded oncogene Latent Membrane Protein 1 (LMP1)13 and is further stimulated by Epstein Barr virus lytic switch protein ZEBRA that also stimulates LMP1 re-expression.14 LMP1 of Epstein Barr virus has been found to induce MMP9 expression in lymphoblastoid cells14 and in EBV-immortalized B lymphocytes.15 In particular, nuclear antigen 3C (EBNA3C) proteins up-regulate MMP9 in B cells from Burkitt’s lymphoma16 and multiple myeloma cells.17 MMP9 proteins might interact with Epstein Barr virus proteins to maintain a microenvironment supporting the proliferation of B cells18 and, in particular, angiogenesis through migration, proliferation and invasion of B-CLL cells might play an important indirect role in malignant growth of chronic lymphocytic leukemia disease.19
TP73 belongs to the transcription factor p53 family (together with p63) and maps to the p36 region of chromosome 1, a region frequently deleted in neuroblastoma and other tumors20 and affected by chromosomal aberrations, including insertions, deletions, and Copy Number Variation.21 To our knowledge, none of the TP73 polymorphisms that we reported here have been examined in previous epidemiological studies. In relation to Epstein Barr virus, molecular studies have shown that silencing of TP73 by hypermethylation is associated with EBV-related gastric cancer while no consistent association is observed in EBV-negative gastric cancer.22 EBV-associated lymphoid malignancies such as Burkitt’s lymphoma23 have also been associated with methylation of TP73.
It has to be stressed that our results might simply reflect the play of chance given the large number of performed tests. Further potential limitations of this study include hospital-based controls which could introduce some selection biases. The strengths of the EpiLymph study lie in its multi-center dimension that allows validation across countries with large differences in ab_EBV positivity. Unfortunately, none of the SNPs that have been associated with chronic lymphocytic leukemia in recent Genome Wide Association studies was genotyped for this study.24 None of the four SNPs was found to be associated with chronic lymphocytic leukemia as a main effect (Online Supplementary Appendix 9), suggesting that the presence of both risk factors (genetic and environmental) might be necessary to act on the disease. Analyses of interactions are crucial to better understand the complexity of the etiology of chronic lymphocytic leukemia. Our preliminary findings, therefore, suggest novel pathways in lymphoma-genesis and deserve to be validated.
Acknowledgments
The authors would like to thank all those who took part in this study, providing questionnaire data and blood samples. We thank the International Agency for Research on Cancer for the DNA preparation.
Footnotes
- Funding: this EpiLymph work was supported by a European Commission Grant n. QLK4-CT-2000-00422 (EpiLymph studies); Association pour la Recherche contre le Cancer (ARC no.5111); the Fondation de France (1999 008471; EpiLymphFrance); Compagnia di San Paolo di Torino, Programma Oncologia 2001 (EpiLymph-Italy); Italian Ministry of Education, University and Research (PRIN 2007, 2007WEJLZB); Spanish Ministry of Health grant PI040091 and RCESP 03/09 (EpiLymph-Spain); Marato TV3 (051210) (EpiLymph-Spain); the German Federal Office for Radiation Protection (StSch4261 and StSch4420; EpiLymph-Germany). Anna González-Neira was supported by the Foundation Genoma España.
- The online version of this article has a Supplementary Appendix.
- Authorship and Disclosures The information provided by the authors about contributions from persons listed as authors and in acknowledgments is available with the full text of this paper at www.haematologica.org. Financial and other disclosures provided by the authors using the ICMJE (www.icmje.org) Uniform Format for Disclosure of Competing Interests are also available at www.haematologica.org.
- Received July 29, 2010.
- Revision received October 28, 2010.
- Accepted November 3, 2010.
References
- De Sanjosé SS, Benavente Y, Nieters A, Foretova L, Maynadié M, Cocco PL. Epstein-Barr virus infection and risk of lymphoma: immunoblot analysis of antibody responses against EBV-related proteins in a large series of lymphoma subjects and matched controls. Int J Cancer. 2007; 121(8):1806-12. PubMedhttps://doi.org/10.1002/ijc.22857Google Scholar
- Bertrand KA, Birmann BM, Chang ET, Spiegelman D, Aster JC, Zhang SM, Laden F. A prospective study of Epstein-Barr virus antibodies and risk of non-Hodgkin lymphoma. Blood. 2010. Google Scholar
- Glaser R, Padgett DA, Litsky ML, Baiocchi RA, Yang EV, Chen M. Stress-associated changes in the steady-state expression of latent Epstein-Barr virus: implications for chronic fatigue syndrome and cancer. Brain Behav Immun. 2005; 19(2):91-103. PubMedhttps://doi.org/10.1016/j.bbi.2004.09.001Google Scholar
- Jaffe ES. WHO Classification of Tumours-Tumours of Haematopoietic and Lymphoid Tissues. IARC Press: Lyon; 2001. Google Scholar
- Berrington de González A, Cox DR. Additive and multiplicative models for the joint effect of two risk factors. Biostatistics. 2005; 6(1):1-9. PubMedhttps://doi.org/10.1093/biostatistics/kxh024Google Scholar
- Kostareli E, Hadzidimitriou A, Stavroyianni N, Darzentas N, Athanasiadou A, Gounari M. Molecular evidence for EBV and CMV persistence in a subset of patients with chronic lymphocytic leukemia expressing stereotyped IGHV4-34 B-cell receptors. Leukemia. 2009; 23(5):919-24. PubMedhttps://doi.org/10.1038/leu.2008.379Google Scholar
- Tarrand JJ, Keating MJ, Tsimberidou AM, O’Brien S, LaSala RP, Han XY. Epstein-Barr virus latent membrane protein 1 mRNA is expressed in a significant proportion of patients with chronic lymphocytic leukemia. Cancer. 2010; 116(4):880-7. PubMedhttps://doi.org/10.1002/cncr.24839Google Scholar
- Stowe RP, Kozlova EV, Yetman DL, Walling DM, Goodwin JS, Glaser R. Chronic her-pesvirus reactivation occurs in aging. Exp Gerontol. 2007; 42(6):563-70. PubMedhttps://doi.org/10.1016/j.exger.2007.01.005Google Scholar
- Shimoyama Y, Asano N, Kojima M, Morishima S, Yamamoto K, Oyama T. Age-related EBV-associated B-cell lympho-proliferative disorders: diagnostic approach to a newly recognized clinicopathological entity. Pathol Int. 2009; 59(12):835-43. PubMedhttps://doi.org/10.1111/j.1440-1827.2009.02466.xGoogle Scholar
- Egeblad M, Werb Z. New functions for the matrix metalloproteinases in cancer progression. Nat Rev Cancer. 2002; 2(3):161-74. PubMedhttps://doi.org/10.1038/nrc745Google Scholar
- Skibola CF, Bracci PM, Halperin E, Nieters A, Hubbard A, Paynter RA. Polymorphisms in the estrogen receptor 1 and vitamin C and matrix metalloproteinase gene families are associated with susceptibility to lymphoma. PLoS One. 2008; 3(7):e2816. PubMedhttps://doi.org/10.1371/journal.pone.0002816Google Scholar
- Bauvois B, Dumont J, Mathiot C, Kolb JP. Production of matrix metalloproteinase-9 in early stage B-CLL: suppression by interferons. Leukemia. 2002; 16(5):791-8. PubMedhttps://doi.org/10.1038/sj.leu.2402472Google Scholar
- Yoshizaki T, Sato H, Furukawa M, Pagano JS. The expression of matrix metalloproteinase 9 is enhanced by Epstein-Barr virus latent membrane protein 1. Proc Natl Acad Sci USA. 1998; 95(7):3621-6. PubMedhttps://doi.org/10.1073/pnas.95.7.3621Google Scholar
- Yoshizaki T, Sato H, Murono S, Pagano JS, Furukawa M. Matrix metalloproteinase 9 is induced by the Epstein-Barr virus BZLF1 transactivator. Clin Exp Metastasis. 1999; 17(5):431-6. PubMedhttps://doi.org/10.1023/A:1006699003525Google Scholar
- Trocmé C, Gaudin P, Berthier S, Barro C, Zaoui P, Morel F. Human B lymphocytes synthesize the 92-kDa gelatinase, matrix metalloproteinase-9. J Biol Chem. 1998; 273 (32):20677-84. PubMedhttps://doi.org/10.1074/jbc.273.32.20677Google Scholar
- Kuppers DA, Lan K, Knight JS, Robertson ES. Regulation of matrix metalloproteinase 9 expression by Epstein-Barr virus nuclear antigen 3C and the suppressor of metastasis Nm23-H1. J Virol. 2005; 79(15):9714-24. PubMedhttps://doi.org/10.1128/JVI.79.15.9714-9724.2005Google Scholar
- Barillé S, Akhoundi C, Collette M, Mellerin MP, Rapp MJ, Harousseau JL. Metalloproteinases in multiple myeloma: production of matrix metalloproteinase-9 (MMP-9), activation of proMMP-2, and induction of MMP-1 by myeloma cells. Blood. 1997; 90(4):1649-55. PubMedGoogle Scholar
- Noël A, Jost M, Maquoi E. Matrix metalloproteinases at cancer tumor-host interface. Semin Cell Dev Biol. 2008; 19(1):52-60. PubMedhttps://doi.org/10.1016/j.semcdb.2007.05.011Google Scholar
- Letilovic T, Vrhovac R, Verstovsek S, Jaksic B, Ferrajoli A. Role of angiogenesis in chronic lymphocytic leukemia. Cancer. 2006; 107(5):925-34. PubMedhttps://doi.org/10.1002/cncr.22086Google Scholar
- Kaghad M, Bonnet H, Yang A, Creancier L, Biscan JC, Valent A. Monoallelically expressed gene related to p53 at 1p36, a region frequently deleted in neuroblastoma and other human cancers. Cell. 1997; 90(4):809-19. PubMedhttps://doi.org/10.1016/S0092-8674(00)80540-1Google Scholar
- Iafrate AJ, Feuk L, Rivera MN, Listewnik ML, Donahoe PK, Qi Y. Detection of large-scale variation in the human genome. Nat Genet. 2004; 36(9):949-51. PubMedhttps://doi.org/10.1038/ng1416Google Scholar
- Ushiku T, Chong JM, Uozaki H, Hino R, Chang MS, Sudo M. p73 gene promoter methylation in Epstein-Barr virus-associated gastric carcinoma. Int J Cancer. 2007; 120(1):60-6. PubMedhttps://doi.org/10.1002/ijc.22275Google Scholar
- Corn PG, Kuerbitz SJ, van Noesel MM, Esteller M, Compitello N, Baylin SB. Transcriptional silencing of the p73 gene in acute lymphoblastic leukemia and Burkitt’s lymphoma is associated with 5′ CpG island methylation. Cancer Res. 1999; 59(14):3352-6. PubMedGoogle Scholar
- Di Bernardo MC, Crowther-Swanepoel D, Broderick P, Webb E, Sellick G, Wild R. A genome-wide association study identifies six susceptibility loci for chronic lymphocytic leukemia. Nat Genet. 2008; 40(10):1204-10. PubMedhttps://doi.org/10.1038/ng.219Google Scholar