MDM2 promotor polymorphism and disease characteristics in chronic lymphocytic leukemia: results of an individual patient data-based meta-analysis

Abstract

A number of single nucleotide polymorphisms have been associated with disease predisposition in chronic lymphocytic leukemia. A single nucleotide polymorphism in the MDM2 promotor region, MDM2SNP309, was shown to soothe the p53 pathway. In the current study, we aimed to clarify the effect of the MDM2SNP309 on chronic lymphocytic leukemia characteristics and outcome. We performed a meta-analysis of data from 2598 individual patients from 10 different cohorts. Patients’ data and genetic analysis for MDM2SNP309 genotype, immunoglobulin heavy chain variable region mutation status and fluorescence in situ hybridization results were collected. There were no differences in overall survival based on the polymorphism (log rank test, stratified by study cohort; P=0.76; GG genotype: cohort-adjusted median overall survival of 151 months; TG: 153 months; TT: 149 months). In a multivariable Cox proportional hazards regression analysis, advanced age, male sex and unmutated immunoglobulin heavy chain variable region genes were associated with inferior survival, but not the MDM2 genotype. The MDM2SNP309 is unlikely to influence disease characteristics and prognosis in chronic lymphocytic leukemia. Studies investigating the impact of individual single nucleotide polymorphisms on prognosis are often controversial. This may be due to selection bias and small sample size. A meta-analysis based on individual patient data provides a reasonable strategy for prognostic factor analyses in the case of small individual studies. Individual patient data-based meta-analysis can, therefore, be a powerful tool to assess genetic risk factors in the absence of large studies.

Introduction

A number of single nucleotide polymorphisms (SNP) of different genes have been associated with disease predisposition.41 While the potential effects of gene polymorphisms are widely acknowledged, most studies investigating the impact of individual SNPs on prognosis in chronic lymphocytic leukemia (CLL) have been controversial.115 None of the allelic variants described so far contribute to relevant current clinical risk models in CLL.1312

Sample sizes of studies to detect associations between SNP markers and clinical end points should be reasonably large, taking into account the very large number of candidate SNPs. In the absence of large studies, individual patient data (IPD)-based meta-analyses can provide a reliable assessment of possibly prognostic SNP markers. However, few IPD-based meta-analyses on SNP data have been performed so far.

A number of years ago, a single nucleotide polymorphism in the MDM2 promotor region (IVS1+309) was discovered and shown to soothe the p53 pathway by influencing MDM2 transcript and protein levels.1514 Patients with the Li-Fraumeni syndrome develop cancers a decade earlier, when carrying the G/G allele of MDM2. These important findings have led to the investigation of the role of the MDM2 polymorphism in a variety of cancers with conflicting results. A combined meta-analysis of breast, colorectal and lung cancers found no effect of the SNP on predisposition in colorectal and breast cancer, but showed a small but significant predisposition to lung cancer.16 Most positive findings on the MDM2 genotype were related to disease onset rather than outcome.16 The exceptionally well documented functional consequences of the MDM2 and the important prognostic role of p53,1817 as well as the particularly striking effects of the SNP in CLL in one study, has led to the study of the MDM2 in many centers, providing conflicting results.1911965

In the current study, we aimed to clarify the prognostic value of MDM2, as well as the correlation between this SNP and the characteristics of CLL, by using an IPD-based meta-analysis.20 To this end, we have collected data from 10 European centers contributing pre-determined clinical and genetic information on their patient cohorts. The combined data set is the largest cohort analyzed for the prognostic impact of genetic risk groups in CLL and thus constitutes a powerful data set to approach prognostic models in this disease.

Methods

Cohorts

This meta-analysis includes individual patient data from 10 cohorts (n=2598) (Figure 1). The cohorts have, in part, been published previously.1195 Two data sets were unavailable either because of regulatory issues21 or due to the decision of the investigator.2322 Tables 1 and 2 show the characteristics of the meta-analysis data set. Details are provided in the Online Supplementary Appendix.

Figure 1.Summary of individual studies as reported.

Table 1.Comparison of clinical and biological characteristics between the ten study cohorts.

Table 2.Comparison of time-to-event characteristics between the ten study cohorts.

Statistical analysis

Meta-analyses are applied to combine the results of related studies to produce an overall effect estimate. The large majority of meta-analyses use summary statistics supplied by the trial investigators or extracted from trial reports. Trial-level results are then combined to produce an overall effect estimate. A more powerful and less biased approach collects individual patient data (IPD) directly from the researchers responsible for each trial. Use of IPD has several advantages over the aggregate data approach, including standardization of statistical analyses, the assessment of potential causes of heterogeneity, and the investigation of interactions and non-linear effects. IPD-based meta-analysis is especially useful when time-to-event data are of interest. For survival analysis, the approach uses individual survival times and takes account of censoring. This provides a more powerful and flexible approach compared to the use of aggregate data since power is lost if time is ignored.

In our study, we apply an approach based on a single proportional hazards regression model stratified by trial. In order to incorporate random effects, a mixed-effect regression model is used. The primary end point of the IPD-based meta-analysis study presented here was overall survival (OS), and the second was time to first treatment (TFT). OS was defined as the time from diagnosis to death (event) or last follow up. TFT was defined as the time interval between the date of diagnosis and date of first CLL-specific treatment or death. TFT was censored at last treatment-free follow-up date. The log rank test, stratified by center, was used to compare survival curves with respect to MDM2 genotypes. For multivariable survival analysis, the proportional hazards regression model of Cox was applied including patients’ sex, age, IGVH gene mutation status and FISH results. For the MDM2, an additive genetic model was used with values 0, 1, and 2, indicating the number of copies of the variant G allele for each patient. A stratified approach was used, allowing each center to have a different base-line hazard. We applied mixed effects Cox proportional hazards models with a random center effect to consider between-center heterogeneity.24 Univariable survival distributions according to genotype are provided as direct adjusted survival curves, adjusted for cohort.25 Forest plots were used to present results for individual centers and the combined analysis. Cochran’s Q was applied to test heterogeneity, whereas the I index was used to quantify the degree of heterogeneity between cohorts. Bias and systematic heterogeneity was checked by contour-enhanced funnel plots, where the standard error is plotted against the log hazard ratio estimate.26 The contour enhanced funnel plots also indicate regions of statistical significance at levels of 0.01, 0.05 and 0.1, testing the null hypotheses that log hazard ratios are equal to zero.

Results

A total of 2598 patients were included in this study, gathered from 10 different single-center series from Europe. The comparison of the base-line characteristics of different cohorts is shown in Table 1, while time-to-event data (survival time, time to first treatment) are shown in Table 2. Median overall survival from diagnosis ranged from a median of 104 months to 301 months. Not surprisingly, there was a difference in outcome of the cohorts with respect to overall survival and time to first treatment reflecting the patient cohort differences and the respective referral patterns. As expected, the most important known genetic risk factors (17p-, 11q-, unmutated IGHV genes) showed hierarchical impact on OS (data not shown).

MDM2SNP309 genotype and disease characteristics

There was no significant difference in genotype distribution among cohorts (P=0.11). Furthermore, the MDM2 was not associated with particular disease characteristics, such as age of onset, 17p-, 11q-, +12q, or IGHV gene mutation status (Table 3). To assess the potential influence of the genotype on the patients’ prognosis, we analyzed TFT and OS using an additive genetic model. Results of the analyses of OS in individual cohorts and the summarizing IPD meta-analysis are illustrated by the forest plot shown in Figure 2. The funnel plot of individual studies was roughly symmetrical, indicating no evidence of publication bias in this study (Figure 3). Results for TFT were comparable. To account for heterogeneity, meta-analysis was carried out stratified for study cohorts. Stratified log rank tests showed no significant differences based on the polymorphism with respect to OS and TFT (P=0.76 and P=0.58, respectively). Specifically, the GG genotype group had a cohort-adjusted median OS of 151 months compared to 153 months for the TG genotype group and 149 months for the TT genotype group (Figure 4). Cohort-adjusted median TFT was 43 months for the GG genotype group, and 42 months for the TG and the TT genotype group.

Figure 2.Forest plot of the IPD-based meta-analysis for overall survival. The plot shows the individual results of univariable Cox regression models for each study cohort together with the meta-analysis result with respect to the MDM2^SNP309 genotype using an additive genetic model. The estimated hazard ratios, HR, together with their 95% confidence intervals are shown numerically and graphically. In addition the numbers of events, D, and observations, N, as well as the proportions of the three genotypes are provided. The results of individual studies are shown as squares centered on each cohort’s point estimate of the hazard ratio (HR). Size of a square represents the precision of the estimate (the inverse of the squared standard error). 95% confidence intervals are represented by horizontal lines. The overall HR estimate from the IPD meta-analysis and its confidence interval are shown at the bottom, represented by a diamond. Testing heterogeneity resulted in Q(df=9) = 14.9, P = 0.10, with I² =39%.

Figure 3.Contour-enhanced funnel plot of the OS analysis. The black dots indicate the cohort-specific estimates of the log hazard ratio and their standard errors; the vertical dashed line shows the summary estimate. Cohorts are indicated by the first two letters of their names. The shaded regions reveal the P value in each study relating to the test of the null hypothesis that the log (hazard ratio) is zero. The funnel plot of the OS analyses with respect to MDM2^SNP309 genotypes is roughly symmetrical, indicating no evidence of publication bias in this study.

Figure 4.Estimated distributions of overall survival. The solid lines show the direct adjusted estimates of survival curves corresponding to MDM2 genotype. Adjustment was made with respect to the different cohorts of the study. Dashed lines represent the unadjusted Kaplan-Meier estimates.

Table 3.Association of the MDM2^SNP309 genotype and patients’ characteristics.

Because of the known interaction of MDM2 with the p53 pathway, we also assessed the impact of the SNP in the group with and the group without 17p deletion. In the subset of patients without 17p deletion, the survival distribution in the sets defined by the three genotypes were comparable (stratified log rank test: P=0.58). Accordingly, no statistically significant impact of the genotype could be demonstrated in the subset with 17p deletion (stratified log rank test: P=0.56). Similar results were obtained for TFT.

Multivariable meta-analysis

We performed a cohort-stratified Cox proportional hazards regression analysis to assess the impact of the MDM2 genotype on overall and treatment-free survival considering known prognostic factors. The first model included MDM2 genotype, sex, age and IGHV gene mutation status as well as the stratification according to cohorts. As shown in Table 4, advanced age, male sex and unmutated IGHV genes were significantly associated with inferior survival but not the MDM2 genotype (P=0.31). Results for TFT were comparable.

Table 4.Cox regression models for overall survival, stratified by study center. Stratified Cox proportional hazards model including covariates age, sex, IGHV mutation status and MDM2^SNP309 (10 centers) and extended model, additionally including cytogenetic aberrations (17p-,11q- and +12q; 9 centers).

In a separate model, we included cytogenetic aberrations (17p-, 11q-, +12q). In this model, sex, age, IGHV gene mutation status, and cytogenetic aberrations were associated with increased risk of death. Again, the MDM2 genotype was not significantly associated with overall survival (P=0.28; Table 4). Considering treatment-free survival, age (P=0.10) and the MDM2 genotype (P=0.85) were the only factors that were not statistically significant.

Discussion

The current study used individual patient data to assess the role of a particular host genetic factor in CLL. In contrast to other large studies, we focused on the prognostic impact of the MDM2 polymorphism. The size of the cohort allowed us to perform subgroup analysis with sufficient statistical power to account for the biological heterogeneity of CLL.

In response to cellular stress, the p53 protein is stabilized and regulates the activity of key effectors of cellular processes, such as DNA repair, cell-cycle arrest, senescence, and apoptosis. These p53-mediated responses are crucial both in reducing cancer frequency and in mediating the response of commonly used cancer therapies. Not surprisingly, p53 loss or mutation leads to very poor outcome in hematologic malignancies, and the impact is particularly striking in CLL. There is a large body of evidence suggesting that the p53 pathway harbors functional inherited single-nucleotide polymorphisms (SNPs) that affect p53 signaling in cells, resulting in differences in cancer risk and clinical outcome in humans (reviewed by Grochola et al.27). Recent elegant work with mice carrying the MDM2 allele showed that MDM2 cells exhibit elevated mdm2 levels, reduced p53 levels, and decreased apoptosis with the MDM2 allele potentiating the tumor phenotype, and altered the tumor spectrum in mice inheriting a p53 hot-spot mutation.28

Although the number of studies on the prognostic impact of individual SNPs in cancer, including CLL, is rapidly increasing, the validation of the results obtained in these studies is often disappointing. Meta-analysis has become an essential part of the process of SNP evaluation. In CLL, IPD-based meta-analysis have so far been limited to clinical trial analysis focusing on clinical questions such as early versus late treatment,29 immunoglobulin prophylaxis30 or response to treatment.3231 This is, therefore, the first IPD based meta-analysis in which the clinical and prognostic significance of a biomarker such as MDM2 gene polymorphism is evaluated. To that purpose, we have compiled a unique data set of 2598 patients from 10 centers using IPD-based meta-analysis to analyze the association between MDM2 polymorphism, genetic and clinical features and clinical outcome in CLL. We observed a strong heterogeneity in the cohorts studied, likely to be due to referral patterns and health care differences across countries. The incidence of the MDM2 genotype in different studies argues for comparability of the different techniques used to detect it.

In contrast to single center studies, we did not find a correlation between MDM2, disease characteristics and outcome, including patients’ age at diagnosis, genetic lesions, IGHV mutations, TFT or overall survival. On the other hand, the current study provides a unique dataset to estimate the impact of genetic aberrations and IGHV mutational status for outcome of patients with CLL. There are additional smaller studies that we have not been able to include in the analysis for various reasons;332321 based on the limited number of cases involved the omission is unlikely to affect the results.

In our analysis, we have used an additive genetic model considering the number of copies of the variant G allele and stratification by study cohort. We have also analyzed the data using mixed effects Cox models with random cohort effect providing comparable results (data not shown). One shortcoming of our study is the lack of data on TP53 mutations, which are now accepted as markers of poor prognosis in CLL but which were not widely available when the series analyzed here were investigated. The analysis within the group of patients with 17p deletion did not show an impact of the MDM2 suggesting that, in CLL, interaction with p53 status is unlikely. While no sample exchange and validation experiments were performed, the incidence of the SNP allele distribution was similar across centers suggesting that technology reproducibility was not a confounding factor. While comprehensive data on the mRNA or Protein expression of mdm2 in CLL are missing, a larger study detected no impact of the MDM2 on mRNA levels.11 Based on our data, further meta-analyses are encouraged; this will help classify and develop the hierarchy of relevant prognostic factors in this disease. This is particularly important as, with the identification of novel mutations (e.g. SF3B1, NOTCH1, BIRC3, ATM, MYD88), the number of prognostic factors has increased, necessitating large cohorts to avoid over-fitted models.

Footnotes

• Authorship and Disclosures Information on authorship, contributions, and financial & other disclosures was provided by the authors and is available with the online version of this article at www.haematologica.org.

• Received November 20, 2013.
• Accepted May 19, 2014.