Malignant Lymphomas |
1 INSERM U618, European Institute of Peptide Research (IFR23), Rouen;
2 Centre Henri Becquerel, Department of Hematology, Rouen;
3 INSERM U837, Université de Lille 2, IFR 114 and from the Institut de Recherches sur le Cancer, Lille;
4 INSERM U614, Institut Hospitalo-Universitaire de Recherche, Rouen, France
Correspondence: Fabrice Jardin, Department of Hematology, Centre Henri Becquerel, 76000 Rouen, France. E-mail:fabrice.jardin{at}rouen.fnclcc.fr
|
|
|---|
Design and Methods: Two polymerase chain reaction assays (multiplex polymerase chain reaction of short fluorescent fragments, QMPSF) were designed to detect gains or losses of c-REL, BCL6, SIM1, PTPRK, MYC, CDKN2A, MDM2, CDKN1B, TP53 and BCL2. Array comparative genomic hybridization was simultaneously performed to evaluate the sensitivity and predictive value of the QMPSF assay. The biological and clinical relevance of this assay were assessed.
Results: The predictive value of the QMPSF assay for detecting abnormal DNA copy numbers ranged between 88–97%, giving an overall concordance rate of 92% with comparative genomic hybridization results. In 77 cases of diffuse large B-cell lymphomas, gains of MYC, CDKN1B, c-REL and BCL2 were detected in 12%, 40%, 27% and 29%, respectively. TP53 and CDKN2A deletions were observed in 22% and 36% respectively. BCL2 and CDKN2A allelic status correlated with protein expression. TP53 mutations were associated with allelic deletions in 45% of cases. The prognostic value of a single QMPSF assay including TP53, MYC, CDKN2A, SIM1 and CDKN1B was predictive of the outcome independently of the germinal center B-cell like/non-germinal center B-cell like subtype or the International Prognostic Index.
Conclusions: QMPSF is a reliable and flexible method for detecting somatic quantitative genetic alterations in diffuse large B-cell lymphomas and could be integrated in future prognostic predictive models.
Key words: deletions, gains, diffuse large B-cell lymphoma, prognosis.
|
|
|---|
|
|
|---|
Array CGH
The CGH analysis was performed using a high resolution 60-mer oligonucleotide-based microarray that contains ~43,000 probes, with an average spatial resolution of 35 kB (Human genome CGH array 44B, Agilent Technologies, CA, USA). High molecular weight DNA was prepared using the standard method. Restriction was performed as recommended by the manufacturer of the arrays. Tumor DNA was labeled with cyanine-5 (Cy5) and reference DNA (pooled normal DNA, Promega, Madison, WI, USA) was labeled with cyanine-3 (Cy3). To increase the probability of detecting relevant genomic gains or losses involving candidate genes related to the outcome of DLBCL, CGH was performed in 17/77 cases selected on the basis of their particularly unfavorable outcome. The analyses of microarray images were performed with the Agilent CGH analytics 3.4.27 software. Classification as gain or loss was based on identification as such by the CGH plotter and visual inspection of the log2 ratios. Signal log2 ratios greater than 0.25 or less then –0.25 were considered to indicate gains and losses, respectively.
QMPSF assay
QMPSF is a sensitive method for detecting genomic deletions or duplications based on the simultaneous amplification of short genomic fragments using dye-labeled primers under quantitative conditions (patent FR 020924).13,14,16 Polymerase chain reaction (PCR) products were analyzed on a sequencing platform used in the fragment analysis mode in which both peaks heights and areas are proportional to the quantity of template present for each target sequence. We designed two distinct QMPSF assays which contain the following target genes: Assay 1 - MYC (8q24), TP53 (17p13), CDKN2A (9p21), SIM1 (6q16) and CDKN1B (12p13.1); Assay 2 - c-REL (2p13), BCL6 (3q27), PTPRK (6q22), BCL2 (18q21) and MDM2 (12q15). The CECR1 gene, located at 22q11 was chosen as a reference gene, considering the fact that it appears uncommonly affected by aneuploidy or focal gains or losses in our own cytogenetic database and in published DLBCL series.12,17,18 Primer pairs were designed for each of these 11 genes to generate PCR fragments ranging from 150 to 250 base pairs and chosen in a way that they do not encompass polymorphisms (Online Supplemental Table S1). PCR were run from 100 ng of genomic DNA in a final volume of 25 µL with 0.16 mmol/L of each deoxynucleoside triphosphate, 1.5 mmol/L MgCl2, 1 unit of thermoprime Plus DNA polymerase (AB gene, Epson, United Kingdom), 5% DMSO and 0.5 to 1.6 µmol/L of each primer, one primer of each pair carrying a 6-FAM label. After initial denaturation for 3 min at 94°C, 20 cycles were performed consisting of denaturation, 94°C for 15 sec, annealing 90°C for 15 sec (ramping 3°C/sec) and extension 70°C, 15 sec (ramping 3°C/sec, followed by a final extension step for 5 min at 70°C). Two control DNA were used (commercial DNA, Roche and a DNA extracted from a reactive lymph node) to calculate the mean normal/tumoral peak height ratio. Using this approach, we demonstrated that TP53 and MDM2 somatic defects could be reliably detected when the proportion of tumoral cells was as low as 20%.15 In addition, polymorphic gene copy number changes were excluded in some cases using matched non-tumoral DNA as control.
QMPSF validation
Considering array-CGH as the reference method, QMPSF and array CGH were both performed in 17 cases. A correlation between CGH log2 and QMPSF ratio was established and allowed the sensitivity, specificity, positive and negative predictive values of the QMPSF assay to be determined. To determine the reliable QMPSF ratio for detecting gene gains or losses, the equation of the regression curve obtained was used to deduce the QMPSF ratio cut-offs corresponding to a CGH log2 ratio of –0.25 (loss) and +0.25 (gain).
TP53 mutational status
To investigate the frequency of TP53 mutations, the highly conserved exons 5 to 8 (central core domain) were screened for the mutation, as described elsewhere.19
Immunohistochemistry
Immunohistochemical studies were performed on formalin-fixed, paraffin-embedded tissue sections of lymph node (n = 76) or spleen tissue (n = 1) using antibodies directed against BCL2, p53 (Dako), BCL6, p27 (Novacastra), p16 (Biocare Medical), and c-REL (Calbiochem). Cases were classed as expressing BCL2 and c-REL if the protein was detected in > 50% tumor cells, and p27 and p16 positive if the protein was detected in > 10% tumor cells. GCB and non-GCB phenotypes were defined using the decision tree established by Hans with the same cut-offs.20
Statistic analysis
The linear relationship between QMPSF ratio and CGH fluorescence ratio was established using Pearsons coefficient (R). Overall survival was measured from the time of diagnosis to the date of death or last follow-up alive. Progression-free survival was calculated from the initiation of the treatment to the date of relapse, progression or death from any cause. Progression-free survival and overall survival rates were estimated by the Kaplan-Meier method and statistical differences were assessed by the log-rank test. A Fishers exact test was used to evaluate the unequal distribution of the different genetic abnormalities and the GCB/non-GCB phenotype and to correlate protein expression and allelic status. A multivariate analysis using a Cox model was conducted to assess the independent prognostic influence of the International Prognostic Index and QMPSF score. Analyses were performed using StatView® and SEM software.21
|
|
|---|
QMPSF assays
QMPSF and CGH experiments were both performed in 17 cases (Figure 1). The linear correlation between CGH log2 ratio and QMPSF ratio is illustrated in Figure 2A. The R2 coefficient was 0.70, indicating a good concordance between the results of the two experiments. From the curve equation, a gene loss, defined by a mean CGH log2 ratio < –0.25, corresponds to a QMPSF ratio below 0.83. A gain detected by a mean CGH log2 ratio > +0.25 corresponds to a QMSPF ratio above 1.13. To maximize detection of true losses and gains, the ratio values finally used were 0.7 and 1.2, respectively. With these cut-offs, DNA copy number changes were confirmed in 156/170 amplicons, giving an overall concordance rate of 92%. The positive and negative predictive values of QMPSF for detecting gains were 88% and 97%, respectively. Similarly, the positive and negative predictive values of QMPSF for detecting gene losses were 90% and 97%, respectively.
![]() View larger version (47K): [in a new window] [Download PPT slide] |
Figure 1. Results of QMPSF and CGH experiments in a single DLBCL case (A) Array CGH shows (from the left to the right) a chromosome 8 in germline configuration, a gain of the long arm of chromosome 17, a short deletion of the 9p21 band, trisomy 12 and trisomy 6. Horizontal blue lines on the chromosome pictogram indicate the location of the five genes included in the QMPSF assay. (B) The corresponding QMPSF assay enables, in a single PCR assay, the detection of gains or losses of gene copy numbers of MYC (8q24), TP53 (17p13), CDKN2A (9p21), CDKN1B (12p13), and SIM1 (6q16). Tumor and normal DNA electropherograms are indicated in blue and orange, respectively. Amplicons are separated and identified by their respective expected sizes (indicated on the upper abscissa). Peak height ratios between tumor and normal DNA are indicated for each gene. CERC1, located on chromosome 22 is used as the reference gene to normalize peak ratios (1.00). MYC and TP53 QMPSF ratios indicated no gene copy number abnormality. A decrease of the QMPSF CDKN2A ratio below a cut-off of 0.7 (0.32) corresponds to the 9p21 loss detected by CGH. In contrast, trisomy 12 and 6 are detected by an increase of the QMPSF ratios (>1.2) of CDKN1B and SIM1, respectively.
|
![]() View larger version (64K): [in a new window] [Download PPT slide] |
Figure 2. Correlation between the QMPSF assay and CGH experiments in de novo DLBCL (A) The correlation curve obtained was used to deduce the most reliable minimal QMPSF ratio cut-offs corresponding to a CGH log2 ratio of –0.25 and +0.25 (gray areas) to detect gene losses and gains, respectively. To maximize detection of true gene copy number changes, QMPSF ratio cut-offs of 0.7 and 1.2 were finally used (horizontal red lines). (B) Correlation between the QMPSF assay and CGH experiments for the CDKN2A gene. Coding CDKN2A gene protein (p16ink4a) detected by immunohistochemistry is indicated for each individual cases (positive cases, gray dots; negative cases, black dots). (C) Distinct pattern of genomic losses involving the CDKN2A gene (p16ink4a) located on chromosome 9p21 in nine patients with DLBCL
|
Frequencies of gain and loss in the overall DLBCL population
Gain and loss frequencies of the ten genes analyzed by the two QMPSF assays are indicated in Table 1. Some target genes were mainly gained such as MYC, CDKN1B, MDM2, c-REL, or BCL2. By contrast, CDKN2A, TP53, SIM1 and PTPRK were almost exclusively deleted. CDKN2A loss was observed in 36% of DLBCL. In 13 cases, the QMPSF ratio was below 0.45, corresponding to a CGH log2 ratio < –0.5, indicating a homozygous deletion. SIM1 (6q16) and PTPRK (6q22) deletions were more frequently observed in MUM1-positive DLBCL (p=0.004 and 0.008, respectively) and in the non-GCB DLBCL subgroup. c-REL gains were observed in 13/31 (42%) GCB-DLBCL and in 8/46 (17%) non-GCB DLBCL (p=0.02). Gain of BCL2 copy was more frequently observed in the non-GCB subtype. Among the 13 cases with homozygous CDKN2A deletion, 11 belonged to the non-GCB subtype, and two to the GCB subtype (p=0.06). However, some simultaneous gene copy number abnormalities also occurred frequently in the same tumors, even if these gene are not located on the same chromosome. For instance, 11/16 cases (68%) with BCL6 (3q27) gains also had BCL2 (18q21) gains (p=0.0003). Multiple losses of tumor suppressor genes can be observed. For instance, in five cases (6%), concomitant loss of TP53 and CKND2A was observed. Furthermore, CDKN2A allelic loss was associated with SIM1 and PTPRK deletions (p=0.003). The details of allelic status of each gene are indicated in the supplemental data (Online Supplemental Table S2). Furthermore, using matched non-tumoral DNA as a control, we confirmed that gene copy number changes were not inherited polymorphisms (data not shown).
|
View this table: [in a new window] [Download PPT slide] |
Table 1. Allelic status assessed by QMPSF in DLBCL and correlation with protein expression assessed by immunohistochemistry and GCB/non GCB subtypes.
|
Prognostic significance of QMPSF assays
The prognostic value of the QMPSF assays was analyzed for each individual gene or using a scoring system established as the amount of gene number abnormalities detected by a single QMPSF assay. Only gain of CDKN1B (12p13.1) was related to significantly shorter progression-free and overall survival. The 3-year progression-free survival rate was 34% for patients with gain of CDKN1B, (16–46% CI-95%), and 67% (52–79%, CI-95%) for patients with a germline configuration (p=0.001). Given the fact that tumor behavior is considered to be the result of multiple gene alterations, we assessed the additive or synergic effect of each gene gain or loss determined by a single QMPSF assay. Each abnormality was scored as +1 and was included in a scoring system. The addition of four abnormalities, MYC and CDKN1B gains, TP53 and CDKN2A losses, was the most powerful scoring system to predict the outcome. The progression-free and overall survival were significantly shorter for patients with a QMPSF score >1 than for patients whose tumors were scored 0–1. This prognostic value held true for the GCB and the non-GCB subtypes. The score remained significantly predictive of a shorter overall survival in the high risk group and tended to be predictive in the low risk group (Figure 3). The prognostic value of the score also held true when only patients treated with CHOP/CHOP-like regimens (excluding patients treated with autologous stem cell transplantation or rituximab as frontline therapy) were considered. In this subgroup (n=58), the 3-year progression-free survival rate was 64% (48–77%, CI-95%) for patients with a score of 0–1 and 16% (5–38%, CI-95%) for patients with a high QMPSF score ( p <0.001). Similarly, the 3-year overall survival rate was 67% (51–79%, CI 95%) for the low score QMPSF group and 26% (12–49%, CI-95%) for patients with a score > 1 (p=0.001).
![]() View larger version (31K): [in a new window] [Download PPT slide] |
Figure 3. Kaplan-Meier survival curves according to the QMPSF score for overall and progression-free survival of patients with DLBCL (A) Kaplan-Meier analysis for the overall population (n=77). (B–C) Kaplan-Meier analysis for GCB (n=31) and non-GCB (n=46 ) DLBCL subtypes. (D–E) Kaplan-Meier analysis for patients defined by the International Prognostic Index as being at low risk (0–2) or high risk (3–5).
|
![]() View larger version (26K): [in a new window] [Download PPT slide] |
Figure 4. Impact of QMPSF/IPI combined score on survival of patients with DLBCL. International Prognostic Index (IPI) and QMPSF scores define three prognostic groups with different prognoses. The group with a favourable prognosis includes patients with both low QMPSF and IPI scores (n=34). The group with an unfavourable prognosis is defined by patients with both high QMPSF and IPI scores (n=17). The intermediate group comprises patients with low IPI/high QMPSF scores or high IPI/QMPSF scores (n=24).
|
|
|
|---|
Previous studies have reported comparative genome analyses of DLBCL, showing relevant differences in the genomic imbalance patterns of the activated B-cell-like and GCB subgroups.10,18,22 In the present study, we confirmed most of the genomic imbalances previously reported. For instance, it was recently shown that 6q21–22 losses, and 3q27 or 18q21-22 gains were more frequently observed in the activated B-cell-like subtypes, corresponding to PTPRK loss, and BCL6 or BCL2 gains, respectively, detected by QMPSF.10 A gain of c-REL copy number was observed in 29% of cases, a rate similar to that detected by CGH or by Southern blot.23–25 This gain is predominantly but not exclusively observed in the GCB subtype.2,22 Interestingly, c-REL copy number gain tends to be more frequently associated with protein expression only in the GCB subtype, indicating that c-REL protein deregulation pathways may be distinct in the two DLBCL subtypes. This observation was recently suggested by the correlation between chromosomal copy number changes and mRNA levels, revealing that genomic copy number gains in 2p14-16, 12q12-15, 3q27-qter and 18q21-q22 lead to subtype-specific up-regulation of genes located in these regions.10
Deletions of the CDKN2A gene were detected in 36% of the DLBCL cases. Tagawa and co-workers reported a more frequent loss of 9p21 in the activated B-cell-like group.11 Here we demonstrated at the gene resolution level that 9p21 loss detected by CGH involved the CDKN2A gene. The lack of p16ink4a protein expression most frequently observed in cases of gene deletion indicated that this mechanism contributes, possibly in combination with methylation, to the down-regulation of this tumor suppressor gene.
To determine the biological relevance of TP53 deletions detected by QMPSF, TP53 mutation status of the central core binding domain was simultaneously assessed. Interestingly, as reported in a large series of NHL of various histology, we observed that TP53 mutations were mainly present in tumors in which allelic loss had occurred.26 This observation indicates that p53 mutations play both recessive and dominant roles in lymphoma.
We demonstrated that a single PCR assay, based on the gene copy numbers of TP53, MYC, CDKN2A and CDKN1B, had a prognostic value independent of the International Prognostic Index and added to its predictive power. A validation of our QMPSF score in an independent set of patients treated uniformly with regimens including rituximab is now necessary. It is likely that genes that are predictive in DLBCL mainly treated first line without rituximab will have a different impact in patients treated with rituximab plus chemotherapy.27,28 Because our approach is very flexible, the QMPSF assay can be easily upgraded and adapted to incorporate additional genes more predictive in this setting. For instance, it was shown that patients whose tumors have 3p11-12 gains have a worse prognosis and new potential tumor suppressor genes have been recently identified.10,29 These data should be integrated to establish a second generation QMPSF assay to predict the outcome of lymphoma treated by immunochemotherapy.
The online version of this article contains a supplemental appendix.
All the authors made substantial contributions to the conception and design of the study, or acquisition of data, or analysis and interpretation of data, drafting the article or revising it critically for important intellectual content; all approved of the version to be published. The authors reported no potential conflicts of interest.
Received for publication September 17, 2007. Revision received October 29, 2007. Accepted for publication November 22, 2007.
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||