Abstract
We analysed the molecular genetic profiles of breast cancer samples before and after neoadjuvant chemotherapy with combination doxorubicin and cyclophosphamide (AC). DNA was obtained from microdissected frozen breast core biopsies from 44 patients before chemotherapy. Additional samples were obtained before the second course of chemotherapy (D21) and after the completion of the treatment (surgical specimens) in 17 and 21 patients, respectively. Microarray-based comparative genome hybridisation was performed using a platform containing ∼5800 bacterial artificial chromosome clones (genome-wide resolution: 0.9 Mb). Analysis of the 44 pretreatment biopsies revealed that losses of 4p, 4q, 5q, 12q13.11–12q13.12, 17p11.2 and 17q11.2; and gains of 1p, 2p, 7q, 9p, 11q, 19p and 19q were significantly associated with oestrogen receptor negativity. 16q21–q22.1 losses were associated with lobular and 8q24 gains with ductal types. Losses of 5q33.3–q4 and 18p11.31 and gains of 6p25.1–p25.2 and Xp11.4 were associated with HER2 amplification. No correlations between DNA copy number changes and clinical response to AC were found. Microarray-based comparative genome hybridisation analysis of matched pretreatment and D21 biopsies failed to identify statistically significant differences, whereas a comparison between matched pretreatment and surgical samples revealed a statistically significant acquired copy number gain on 11p15.2–11p15.5. The modest chemotherapy-driven genomic changes, despite profound loss of cell numbers, suggest that there is little therapeutic selection of resistant non-modal cell lineages.
Similar content being viewed by others
Main
Breast cancer is a heterogeneous disease comprising tumours with remarkably distinct clinical behaviour (Simpson et al, 2005). In the last 20 years, a multitude of prognostic and predictive markers have been tested, however only oestrogen receptor, progesterone receptor and HER2 are currently used to tailor the therapy of breast cancer patients (Goldhirsch et al, 2005). The use of expression arrays to derive molecular profiles that are predictive of clinical outcome has received great attention in the last 5 years (van’t Veer et al, 2005; Reis-Filho et al, 2006b), especially in the neoadjuvant setting (Davidson and Morrow, 2005). Predictive signatures for different neoadjuvant chemotherapy regimens have been reported by some investigators but with less success by others (Chang et al, 2003; Ayers et al, 2004; Cleator and Ashworth, 2004; Chang et al, 2005a; Hannemann et al, 2005; Iwao-Koizumi et al, 2005; Modlich et al, 2005; Rouzier et al, 2005; Dressman et al, 2006; Reis-Filho et al, 2006b). Although some of these results are promising, issues related to the instability of mRNA, experimental design and data analysis have led many to call into question the validity of current approaches (Ransohoff, 2004; Brenton et al, 2005; Ioannidis, 2005; Reis-Filho et al, 2006b).
Chromosomal comparative genomic hybridisation (cCGH) has been widely used to analyse the pattern of unbalanced genomic aberrations in breast cancer (Roylance et al, 1999; Buerger et al, 1999a, 1999b; Albertson, 2003; O’Connell, 2003; Reis-Filho et al, 2005a). More recent studies have employed higher resolution array-based CGH (aCGH) showing the enormous complexity of breast cancer genomes. Nonetheless, these studies have consistently reported the same regions of gain (1q, 8q, 11q, 17q, 20q) and loss (8p, 9p, 13q, 16q) (Rennstam et al, 2003; Loo et al, 2004; Naylor et al, 2005; Nessling et al, 2005). Unbalanced chromosomal aberrations and changes in gene copy numbers have been reported as a mechanism for acquired drug resistance to chemotherapy in cell line studies (Leyland-Jones et al, 1999; Yasui et al, 2004). However, no genome-wide analysis addressing response prediction by aCGH in breast cancer patients has been published.
The primary aim of this study was to determine whether aCGH could be used to identify unbalanced genomic changes predictive of response to preoperative doxorubicin and cyclophosphamide (AC) combination. We also sought to investigate the hypotheses that resistance to chemotherapy could be driven either by selection of chemotherapy resistant populations of neoplastic cells or by the induction of specific genomic aberrations.
Materials and methods
Patients and samples
Patients were recruited and treated at the Royal Marsden Hospital (Sutton and London, UK). Eligible patients were those undergoing neoadjuvant adriamycin and AC chemotherapy treatment at doses of 60 and 600 mg m−2, respectively, three weekly, for a clinically measurable breast cancer. Approval by the Royal Marsden Hospital Clinical Research and Ethics committees was granted and written consent was obtained in all cases. Patients were offered neoadjuvant treatment for one of several standard indications including locally advanced or inflammatory breast cancer, high tumour to breast size ratio and tumours located close to the nipple. The majority of the patients were from a previously published study on RNA expression profiling (Cleator et al, 2006).
Diagnosis was confirmed on core-cut biopsies subjected to routine histological analysis. Patients who demonstrated at least a partial clinical response received six cycles of treatment before local treatment. Patients with no or only marginal response after three or four cycles, proceeded directly to local treatment or were commenced on alternative systemic treatment (docetaxel).
Clinical size of tumour was recorded as the largest diameter and a diameter perpendicular to this. This was recorded before commencement and at completion of treatment. Clinical response was categorised as follows: no palpable abnormality after treatment, complete clinical response (cCR); >50% reduction in the product of the bidimensional measurements, partial response (PR); <50% reduction in the product of bidimensional measurements was recorded as no change (NC); residual ill-defined thickening after a good response, minimal residual disease (MRD) (Cleator et al, 2006). No patients in the study demonstrated progressive disease, an increase of more than 25% in the product of bidimensional measurements. Those cases in which there was no residual invasive carcinoma at surgery were classified as a complete pathological response (pCR). Good responders were defined as pCR, cCR, MRD; poor responders were defined as PR or NC. Eight patients undergoing a complete clinical and radiological (on ultrasound) response received radiation only as local treatment. Therefore some of the cCRs may represent undocumented pCRs.
Research 14-gauge core biopsies were collected before commencing treatment and again at 21 days after first treatment. Samples were snap frozen in liquid nitrogen. All samples were thereafter coded using a study number as an identifier. Frozen cores were embedded in optimum cutting temperature embedding compound and sectioned at −20°C in a cryostat. Sections (5 μm thickness) were taken for hematoxylin and eosin staining to assess histological character superficially from the core as soon as ‘full-face’ was reached. The percentage of cells comprising invasive malignant disease was recorded. For patients in whom multiple biopsies were available, that with the highest invasive content was used for microarray analysis. Samples with less than 20% invasive cancer content were excluded from the study. Eight micrometer thick, nuclear fast red-stained representative frozen sections from 44 patients with primary breast cancer obtained before a neoadjuvant combination of AC were microdissected. In 22 patients consent for a second research core needle biopsy was granted and the procedure was technically feasible; in 17 of these, samples obtained 3 weeks after first course of chemotherapy (before cycle two) contained proportions of neoplastic cells suitable for microdissection (i.e. >20% of neoplastic cells). Twenty-one surgical formalin-fixed paraffin-embedded resection samples taken after 4–6 cycles of AC and with sufficient amounts of residual neoplastic cells were retrieved from the pathology files of the Royal Marsden Hospital.
Fluorescent in situ hybridisation
Fluorescent in situ hybridisation (FISH) analysis was performed on representative 4-μm sections using PathVysion HER-2 DNA probe mixture containing a HER-2 DNA probe (190 kb Spectrum Orange directly labelled DNA probe) and the CEP 17 DNA probe (5.4 kb Spectrum Green directly labelled fluorescent DNA probe specific for the chromosome 17 α-satellite DNA sequence) as described previously (Dowsett et al, 2003). A total of 60 cells were scored for red and green signals for each section and results were expressed as a ratio of red to green signals, with a ratio ⩾2.0 being considered amplified (Dowsett et al, 2003).
Microdissection and DNA extraction
Frozen and formalin-fixed samples were subjected to aCGH analysis after microdissection of neoplastic cells with a sterile needle under a stereomicroscope (Olympus SZ61) from one to three consecutive 8 μm nuclear fast red-stained sections (median=1). Estimated purity of tumour cells ranged from 75 to 100% (median=90%) (Supplementary Figure 1). DNA was extracted according to standard methods (Reis-Filho et al, 2005b). The DNA yield and purity was assessed by spectrophotometry and the DNA quality (DNA fragment size range) by agarose gel electrophoresis as described previously (Reis-Filho et al, 2005b).
Microarray-based comparative genomic hybridisation
Microarray-based CGH was performed with in-house produced microarrays containing ∼5.8K bacterial artificial chromosome (BAC) clones, spaced out throughout the genome in ∼0.9 Mb intervals. The Breakthrough Breast Cancer Research Centre CGH microarray is composed of the same BAC clones as described by Reis-Filho et al (2005b).
Mapping of the BAC clones was retrieved from public sources and positioned according to the May 2004 build of the human genome sequence (hg17). When genomic positioning was dubious or conflicting, BAC end pair sequencing and FISH mapping was performed. Clones that either (i) showed poor quality end sequences or (ii) hybridised to multiple chromosomal locations or to a cytogenetic location inconsistent with their position in the sequence assembly were excluded from analysis. Bacterial artificial chromosome clones were spotted in triplicate onto Corning GAPSII-coated glass slides (Corning, NY, USA). Labelling of 250 ng of non-amplified DNA obtained from microdissected frozen sections or 1000 ng of DNA retrieved from microdissected formalin-fixed paraffin-embedded tissue sections, hybridisation and washes were carried out essentially as described previously (Reis-Filho et al, 2005b, 2006a; Natrajan et al, 2006).
Arrays were scanned with a GenePix 4000A scanner (Axon Instruments Inc., Union City, CA, USA); fluorescence data were processed with GenePix 4.1 image analysis software (Axon Instruments Inc.) as described previously (Reis-Filho et al, 2005b, 2006a; Natrajan et al, 2006).
Data analysis
The log2 ratios were normalised for spatial and intensity-dependent biases using a two-dimensional Loess local regression. Experiment replicates (‘dye-swaps’) were collated, BAC clone replicate spots averaged, and clones excluded with poor reproducibility between replicates (s.d. >0.2). Samples with >70% missing/poor values were excluded, as were clones with missing/poor values in >70% samples. Missing data were imputed by k-nearest neighbours (Troyanskaya et al, 2001) and clones with no mapping information (May 2004 build of the human genome sequence, hg17) were excluded. A total of 5001 different BAC clones were used in the final analysis.
Statistical analysis
All data transformation and statistical analysis were carried out in R 2.0.1 (http://www.r-project.org/), BioConductor 1.5 (http://www.bioconductor.org/), making extensive use of modified versions of the package aCGH in particular (Paris et al, 2004; Natrajan et al, 2006) and S-Plus (version 6.2.1). The log2 ratio of each BAC clone in each tumour sample was centred by the median of each case. Thresholds for determining copy number gains and losses were determined as exceeding ±3.0 × s.d. of the mean of these control experiments (log2 ratio of ±1). For identification of DNA copy number alterations, data were smoothed using a local polynomial adaptive weights smoothing procedure for regression problems with additive errors (Hupe et al, 2004; Natrajan et al, 2006).
Associations between genomic loci were assessed by calculating Pearson's correlations between thresholded values for each clone, assigned as 1, 0 or −1 for gain, NC, or loss in copy number. Thresholded data for each clone were also used for categorical analysis using a Fisher's exact test adjusted for multiple-testing with a correction for multiple testing using the step-down permutation procedure maxT, providing strong control of the family-wise type I error rate.
Results
Patient, tumour and biopsy characteristics
A total of 44 patients were included in the study. Patient and tumour characteristics are summarised in Table 1 according to clinical response. Median tumour diameter was 4 cm (range 2–10). Two patients had inflammatory breast cancer (T4d). Median follow-up was 24 months (range, 9–38 months).
Of the 44 patients, 24 (55%) demonstrated a ‘good’ and 20 (45%) a ‘poor’ clinical response. Four patients (9%) were documented as undergoing a pathological Complete Response (pCR), 17 (37%) underwent a cCR and three (7%) had ill-defined thickening (MRD) at end of treatment; eight (19%) had a PR and 12 (28%) showed NC.
Validation of aCGH method
DNA was extracted from 61 frozen samples. The median DNA yield obtained per sample was 660 ng (range 152–3188 ng). Owing to constraints in amount of tumour material available from frozen core biopsies after manual dissection, we assessed the reproducibility and quality of profiles obtained with 250, 500 and 1000 ng of tumour DNA. All DNA concentrations showed optimal results, with the Pearson's coefficient of correlation ranging from 0.93 to 0.96 (data not shown). Biological replicates of aCGH experiments also demonstrated high reproducibility: two independent research core biopsies were obtained before neoadjuvant chemotherapy for five patients. DNA extraction and CGH arrays were performed separately for these paired samples. For each of the matched pairs, the Pearson's correlation coefficient ranged from 0.64 to 0.86 (mean 0.78, data not shown). In six cases, HER2 amplification was defined by aCGH. Fluorescent in situ hybridisation analysis confirmed the results in all cases, providing further evidence to support the validity of the aCGH analysis methods employed in this study (data not shown).
For six cases, high-resolution cCGH was performed and the genetic profiles compared with those obtained with aCGH. Correlation was good-to-excellent, with correlations for low-level gains and deletions >20 Mb and any amplification ranging from 60 to 87.5% (median=79.5%, mean=77.8%, data not shown). However, aCGH was more sensitive in detecting small losses and low-level gains than HR-CGH.
Genomic alterations in 44 pre-chemotherapy breast cancer samples
To identify genomic regions harbouring recurrent unbalanced genomic changes, we plotted the frequency of tumours showing gain or loss for each BAC across the genome (Figure 1A). The most frequent (>30%) genomic changes (Table 2) comprised gains of 1q (66%), 5p (32%), 8q (70%), 16p (36%) and 20q (41%) with the smallest regions of genomic gain on 1q31.1–1q31.2 and 1q22–1q25.3, 5p15.31–5p15.33, 8q23.1–8q25.1, 16p11.2–16p12.2 and 20q13.13–20q13.33, respectively. Losses were observed on 4q (39%), 8p (50%), 9p (36%), 11q (32%), 13q (36%), 16q (52%) 17p (50%) and 18q (39%) with the smallest regions of deletion on 4q32.3–4q33, 8p21.1–p23.3, 9p22.2–9p24.3, 11q23.1–11q25, 13q14.11–13q14.3, 16q23.1–q24.3, 17p12–17p 13.3 and 18q22.1–q23. In addition to the large regional alterations, the resolution of the BAC array allowed us to map smaller regions of gain or loss. Bacterial artificial chromosome clones gained or deleted in >30% of the tumours are described in Supplementary Tables 1 and 2, respectively.
Comparison of genomic alterations in different phenotypes of breast cancer
On the basis of three distinguishing phenotypic characteristics (ER, HER2 and histological type), we characterised genetic alterations that might be associated with subtypes of breast cancer on the 44 pretreatment biopsies. Oestrogen receptor (ER)-negative tumours (N=15) were characterised by significantly more chromosomal changes than ER-positive tumours (Table 3 and Supplementary Figure 2). Subgroup analysis showed that the following chromosomal changes were more frequently associated with ER negativity: gains on 1p, 2p, 2q, 6q, 7q, 9p, 11q and 19q with large regions on 1p31–1p34.2, 7q32–7q36 and losses on 4p, 4q, 5p, 5q, 12q and 17q as well as large regions such as the whole long arm of chromosome 5 and whole chromosome 4.
Pearson's correlation analysis of the subset of ER-positive tumours (n=29) revealed chromosomal alterations that coexist in this group of tumours. The heat map in Figure 1B shows regions of positive (change in same direction: red) and negative (change in opposite direction: blue) correlations. Relatively large regions of correlation were seen between changes in copy number of BACs on chromosome 1p and 8p (P<0.01), 3p and 4p (P<0.001), 4q and 8q (P<0.01) and 8q and 12q (P<0.01) were observed, as were strong inverse correlations between 1p and 12q (P<0.0005), 7p and 22 (P<0.001) and 11q and 12q (P<0.001) (Figure 1B). Pearson's correlation matrix of ER-negative tumours revealed strong direct associations between 1p and 7q (P<0.001), 16p and 17q (P<0.001), and strong inverse correlations between 1q and 14q (P<0.005), 2q and 9q (P<0.001), 4q and 15q (P<0.001), 6q and 11q (P<0.001), 7p and 11q (P<0.001), 7p and 12q (P<0.01), 9p and 16p (P<0.001) and 13q and 15q (P<0.005) (Figure 1C).
Comparing invasive lobular carcinomas (n=7) with non-lobular cases (37 cases), 16q loss (from 16q12.1 to 16q24.3) was significantly associated with lobular phenotype (unadjusted Fisher's exact test <0.01) (Supplementary Table 3 and Supplementary Figure 3). This region contains a number of cadherin genes (CDH1, CDH3, CDH5, CDH 8, CDH 11). Gain of 8q21.11–q24.23 was significantly associated with the ductal phenotype (Supplementary Table 3 and Supplementary Figure 3).
Comparison between the six cases with HER2 amplification on 17q11.2 with the 38 HER2-negative cases revealed more frequent gains of 6p25.2–p24.2, 17q12 and Xp22.33–p11.36, and loss of 18p11.31 (Supplementary Table 3 and Supplementary Figure 4). On the other hand, loss of 16q23.3–q24.1 was more frequently observed in HER2 – tumours. This would be expected given that all HER2+cases were of histological grade III and loss of 16q is significantly less frequently found in this group of tumours (Reis-Filho et al, 2005a; Simpson et al, 2005).
Prediction of tumour response
Comparison of the molecular genetic profiles of objective clinical responders (complete and partial clinical responders, N=24) with those of non-responders revealed loss of a 13.3 Mb region on 13q31.1–13q32.2 to be the only significant difference (Figure 2). This region includes the genes SLITRK6, SLITRK5, GPC5, GPC6, DCT, TGDS, SOX21, ABCC4, CLDN10, DZIP1, DNAJC3, UGCGL2, HS6ST3, HSP90AB6P, OXGR1, MBNL2 and RAP2A. In addition, this region also encompasses the micro RNA cluster miR-17-92, which is reported to induce tumour cell growth and to be overexpressed and sometimes amplified in aggressive forms of lung cancer (Hayashita et al, 2005).
Comparison of pre-chemotherapy samples to D21 samples, before cycle two of AC
Unsupervised hierarchical clustering analysis based upon genetic alterations on all chromosomes showed that all pre- and D21 post-biopsies clustered together (N=17 cases) (Figure 3A). There were no significant differences in the profiles between poor (10 cases) and good responders (seven cases) (P=0.95, Kolmogorov–Smirnov test). The comparison of the pattern and frequency of unbalanced genomic changes detected by aCGH analysis in matched pretreatment and D21 biopsies revealed no significant differences (multi-Fisher's exact test, data not shown).
Comparison of pre-chemotherapy samples to surgical samples, after 4–6 cycles of AC
Twenty-one surgical samples out of 44 cases were available (48%). All but three segregated with the respective pre-chemotherapy sample by unsupervised hierarchical clustering (Figure 3B). For 12 cases, pre-, D21 and at surgery samples were available. Unsupervised hierarchical clustering analysis based upon genetic alterations on all chromosomes was performed. All 12 ‘triplicates’ clustered together.
Comparison between the genomic profiles obtained with DNA extracted from matched frozen and formalin-fixed pretreatment core biopsies were performed (n=2). The profiles showed a good concordance (r2>0.8). Subsequently, a comparison of the aCGH profiles of 21-matched pretreatment core biopsies and the respective excision specimens obtained after 4 months of neoadjuvant chemotherapy revealed the presence of gains of 4q13.1 (9/21, 43%), 11p15.2–p15.5 (14/21, 66%), 12q13.3 (14/21, 66%), 18p11.21 (12/21, 52%) and 19q13.2 (12/21 52%) in the excision specimens (Figure 4 and Supplementary Table 4). After FDR adjustment, only the gain on 11p15.2–p15.5 retained its statistical significance. Although we cannot completely rule out that this change could result from the comparison between DNA extracted from frozen and formalin-fixed samples, this 14.3 Mb telomeric region encompasses interesting oncogene candidates, such as H-RAS, a bona fide oncogene, and IGF2, which has been reported to show loss of imprinting in up to 30–60% of breast cancers (McCann et al, 1996; Wu et al, 1997). In the pre-chemotherapy samples, gain on 11p15.2–p15.5 was observed in five out of 10 good responders and in two out of 11 poor responders (not significant, Fisher's exact test).
Discussion
In the present study, the frequency and pattern of unbalanced genomic aberrations were similar to those described in previous studies where genome-wide changes were analysed by means of chromosomal CGH (Tirkkonen et al, 1998; Buerger et al, 1999a, 1999b; Roylance et al, 1999; Seute et al, 2001; Rennstam et al, 2003) and aCGH (Callagy et al, 2005; Naylor et al, 2005; Nessling et al, 2005). Recurrent gains on chromosome 1q, 8q, 11q, 17q and 20q and losses on 6q, 8p, 9p, 13q and 16q were the most prevalent changes. In addition, we confirmed the association between ER positivity and gain of 1q coupled with loss of 16q (Farabegoli et al, 2004; Reis-Filho et al, 2005a; Simpson et al, 2005) and the more prevalent deletions of 4p16 and 4p15, 5q and 17p11.2 in ER-negative tumours (Loo et al, 2004). In contrast to previous studies (Loo et al, 2004), gains of 8q24.1 (MYC) and 17q12 (HER2) were not significantly more frequent in ER-positive tumours. This is expected, given that only breast carcinomas of histological grades 2 and 3 were present in the population and gains of 8q are rather frequent in grade 2 and 3 breast cancers (Buerger et al, 1999a, 1999b; Roylance et al, 1999).
The comparison between the genomic profiles obtained for ductal and lobular carcinomas were also in agreement with previous studies (Buerger et al, 1999b; Shelley Hwang et al, 2004; Reis-Filho et al, 2005a; Simpson et al, 2005; Stange et al, 2006): gain of 1q and deletions of 16q were the most prevalent changes in lobular carcinomas, whereas gain of 8q was significantly more frequent in grade 2 and 3 ductal carcinomas. However, we could define the smallest region of overlap of the deletions of 16q, which mapped to 16q21–q22.1 and encompassed the region of the cadherin gene cluster, and the gain of 8q, which encompassed two regions 8q13.2–q21.13 and 8q21.3–qtel (Supplementary Table 3). The most significant gain comparing ductal carcinoma with lobular carcinomas mapped to 8q24.11 (117.8–118.0 Mb), which encompasses RAD21 and eucaryotic translation initiation factor 3, subunit 3 gamma (EIF3S3). RAD21 is believed to function in sister chromatid alignment as part of the cohesin complex and also in double-strand break repair and influences cellular proliferation (Atienza et al, 2005), whereas EIF3S3 is reported to be amplified and overexpressed in up to 20% of breast carcinomas (Nupponen et al, 1999). Our observations are consistent with previous studies on breast cancer, confirming the robustness of our aCGH protocol, the validity of our analysis method and the likelihood that we have a representative set of tumours.
Several attempts have been made to predict clinical or pathological response to neoadjuvant chemotherapy in breast cancer using gene expression arrays. Gene ‘signatures’ or predictors have been devised for several chemotherapy regimens, including paclitaxel followed by fluorouracil, AC, AC/doxorubicin-docetaxel and taxane only chemotherapy (Chang et al, 2003, 2005b; Ayers et al, 2004; Hannemann et al, 2005; Cleator et al, 2006). Although these results are promising, the exceedingly small sample size and limitations with the current technology and analysis methods have so far precluded definitive conclusions (Brenton et al, 2005; Reis-Filho et al, 2006b).
We have recently reported a similar expression profiling study (Cleator et al, 2006) that included the 44 patients studied here plus a small number of others for whom pretreatment tissue for DNA analysis was not available. Neither unsupervised nor supervised methods could separate the responders from non-responders. In the current study, using the same cohort of patients, which included 24 good clinical responders, aCGH analysis revealed a deletion of a large region from 13q31.1 to 13q34 as the only significant copy number change associated with response to chemotherapy. Although this statistical association lost its significance after correction for multiple comparisons, the correction method, which we adopted may be too conservative, given that unbalanced chromosomal aberrations usually encompass more than one BAC clone. One of the genes deleted in this region is a subtype of heat-shock protein 90, HSP90AB6. HSP 90 is a molecular chaperone whose association is required for stability and function of multiple signalling proteins that promote cancer cell growth and/or survival (Chen et al, 2005). Further studies with a larger sample size focusing on this particular region are warranted. We cannot exclude the possibility that analysis similar to ours of a larger series of samples might identify other DNA aberrations of importance for response to AC chemotherapy. It is however likely that if these were present in a large proportion of the population they would have been uncovered by the present analysis.
Sequential sampling of tumour during neoadjuvant treatment can be used to detect gene expression modification induced by therapy as we have demonstrated with endocrine therapy (Mackay et al, 2005). This was demonstrated for chemotherapy by Hannemann et al (2005), who observed that tumours that responded to neoadjuvant chemotherapy showed dramatic changes in their expression profiles when compared to the changes observed in non-responders (Hannemann et al, 2005). On the other hand, a comparison between the transcriptomic profiles of tumours subjected to taxane-based neoadjuvant chemotherapy before and 3 months after treatment revealed strikingly different patterns, independent of initial sensitivity or resistance (Chang et al, 2005b).
Cell line studies have demonstrated that changes in gene copy numbers may lead to acquired resistance to chemotherapy (Leyland-Jones et al, 1999; Shimizu et al, 2002; Yasui et al, 2004). Specific gains of genetic material mapping to multi-drug resistance (MDR) gene MDR1 locus have been reported in drug resistant cell lines (Shimizu et al, 2002; Kuwano et al, 2003). Recurring amplicon 7q11.2–q21 identified by CGH in doxorubicin-resistant hepatocellular carcinoma cell lines coincided with the localisation of MDR1 (Pang et al, 2005). Induction of DNA damage response genes such as p21 were obtained by doxorubicin in breast cancer cell lines (Troester et al, 2004). Both chromosome 7 alterations and several cytogenetic changes involving the 7q21 locus are associated with the development of MDR in sarcoma cells (Chen et al, 2002). Analysis of genomic amplifications and deletions revealed specific genetic alterations common to both intrinsic and acquired doxorubicin resistance including ABCB1, PGY3 (ABCB4) and BAK (Turton et al, 2001).
Most patients show some tumour shrinkage with neoadjuvant chemotherapy and for many there may be profound loss of malignant cells with some showing a pCR. Characterisation of the residual cells may be expected to a better understanding of the causes of resistance and allow the identification of the means to overcome the resistance. Given the molecular heterogeneity between breast carcinomas and the composition of most breast cancers being of multiple non-modal clones, we hypothesised that the cell loss might lead to selection of breast cancer cell lineages that were resistant to chemotherapy by virtue of specific DNA alterations. Few studies have previously addressed this issue. A very small study, involving CGH of just four tumours from breast cancer patients after neoadjuvant chemotherapy (Fazeny-Dorner et al, 2003) showed typical DNA imbalances for ductal breast cancer. Three patients showed involvement of several regions bearing genes of drug resistance (MDR1, BCRP, MRP1, RFC1); the fourth patient displayed an amplification in the region of MYC. In our study, the hierarchical clustering of sequential samples at baseline, 3 weeks and at surgery indicated that any changes induced in DNA profiles by therapy were modest compared with the differences that were present between the patients. This interpretation was supported by the observation that statistically significant chemotherapy-driven genomic changes were not detected within 3 weeks (i.e. after a single course of chemotherapy). Comparison of the pretreatment and the tumour specimens excised after 4 months of neoadjuvant chemotherapy revealed the appearance of genomic gains on 11p15.2–p15.5. Given the experimental design of our study, we cannot define whether these amplifications were chemotherapy induced or were present in a non-modal population of the primary tumour and selected by chemotherapy, as the majority of surgical samples by definition were obtained from poor responders (66%). Interestingly, pathological relaxation of the imprinting pattern in this region is reported to be found in 30–60% of breast cancers (McCann et al, 1996; Wu et al, 1997). However we did not find a correlation between amplification on 11p15.5 and clinical response (Han et al, 2006). Furthermore, these changes on surgical samples, need to be cautiously interpreted as surgical samples were subjected to different fixation conditions (Devries et al, 2005). Thus in contrast to observations with cancer cell lines, we did not find large or frequent chemotherapy-induced acquired genomic changes. Given that aCGH provides an average of the pattern of genomic gains and losses in the cell population studied, chemotherapy-driven losses or low-level gains occurring in non-modal clones are unlikely to be detected.
In conclusion, array CGH is a powerful method for the genome-wide detection of chromosomal imbalances and allowed us to detect molecular genetic aberrations associated with specific breast cancer subgroups (ductal vs lobular, ER-negative vs ER-positive tumours). A molecular genetic profile specific of good responders to neoadjuvant chemotherapy was not detectable in our series. Chemotherapy-driven genomic changes were not detected following 3 weeks of treatment and only a single change after completion of treatment. The hypothesis of resistance to neoadjuvant chemotherapy by the selection of non-modal cell lineages, which differ by gene amplifications or losses is not supported by our results.
Change history
16 November 2011
This paper was modified 12 months after initial publication to switch to Creative Commons licence terms, as noted at publication
References
Albertson DG (2003) Profiling breast cancer by array CGH. Breast Cancer Res Treat 78: 289–298
Atienza JM, Roth RB, Rosette C, Smylie KJ, Kammerer S, Rehbock J, Ekblom J, Denissenko MF (2005) Suppression of RAD21 gene expression decreases cell growth and enhances cytotoxicity of etoposide and bleomycin in human breast cancer cells. Mol Cancer Ther 4: 361–368
Ayers M, Symmans WF, Stec J, Damokosh AI, Clark E, Hess K, Lecocke M, Metivier J, Booser D, Ibrahim N, Valero V, Royce M, Arun B, Whitman G, Ross J, Sneige N, Hortobagyi GN, Pusztai L (2004) Gene expression profiles predict complete pathologic response to neoadjuvant paclitaxel and fluorouracil, doxorubicin, and cyclophosphamide chemotherapy in breast cancer. J Clin Oncol 22: 2284–2293
Brenton JD, Carey LA, Ahmed AA, Caldas C (2005) Molecular classification and molecular forecasting of breast cancer: ready for clinical application? J Clin Oncol 23: 7350–7360
Buerger H, Otterbach F, Simon R, Poremba C, Diallo R, Decker T, Riethdorf L, Brinkschmidt C, Dockhorn-Dworniczak B, Boecker W (1999a) Comparative genomic hybridization of ductal carcinoma in situ of the breast-evidence of multiple genetic pathways. J Pathol 187: 396–402
Buerger H, Otterbach F, Simon R, Schafer KL, Poremba C, Diallo R, Brinkschmidt C, Dockhorn-Dworniczak B, Boecker W (1999b) Different genetic pathways in the evolution of invasive breast cancer are associated with distinct morphological subtypes. J Pathol 189: 521–526
Callagy G, Pharoah P, Chin SF, Sangan T, Daigo Y, Jackson L, Caldas C (2005) Identification and validation of prognostic markers in breast cancer with the complementary use of array-CGH and tissue microarrays. J Pathol 205: 388–396
Chang JC, Hilsenbeck SG, Fuqua SAW (2005a) The promise of microarrays in the management and treatment of breast cancer. Breast Cancer Res 7: 100–104
Chang JC, Wooten EC, Tsimelzon A, Hilsenbeck SG, Gutierrez MC, Elledge R, Mohsin S, Osborne CK, Chamness GC, Allred DC, O’Connell P (2003) Gene expression profiling for the prediction of therapeutic response to docetaxel in patients with breast cancer. Lancet 362: 362–369
Chang JC, Wooten EC, Tsimelzon A, Hilsenbeck SG, Gutierrez MC, Tham YL, Kalidas M, Elledge R, Mohsin S, Osborne CK, Chamness GC, Allred DC, Lewis MT, Wong H, O’Connell P (2005b) Patterns of resistance and incomplete response to docetaxel by gene expression profiling in breast cancer patients. J Clin Oncol 23: 1169–1177
Chen B, Piel WH, Gui L, Bruford E, Monteiro A (2005) The HSP90 family of genes in the human genome: insights into their divergence and evolution. Genomics 86: 627–637
Chen GK, Lacayo NJ, Duran GE, Wang Y, Bangs CD, Rea S, Kovacs M, Cherry AM, Brown JM, Sikic BI (2002) Preferential expression of a mutant allele of the amplified MDR1 (ABCB1) gene in drug-resistant variants of a human sarcoma. Genes Chromosomes Cancer 34: 372–383
Cleator S, Ashworth A (2004) Molecular profiling of breast cancer: clinical implications. Br J Cancer 90: 1120–1124
Cleator S, Tsimelzon A, Ashworth A, Dowsett M, Dexter T, Powles T, Hilsenbeck S, Wong H, Osborne CK, O’Connell P, Chang JC (2006) Gene expression patterns for doxorubicin (Adriamycin) and cyclophosphamide (Cytoxan) (AC) response and resistance. Breast Cancer Res Treat 95: 229–233
Davidson NE, Morrow M (2005) Sometimes a great notion – an assessment of neoadjuvant systemic therapy for breast cancer. J Natl Cancer Inst 97: 159–161
Devries S, Nyante S, Korkola J, Segraves R, Nakao K, Moore D, Bae H, Wilhelm M, Hwang S, Waldman F (2005) Array-based comparative genomic hybridization from formalin-fixed, paraffin-embedded breast tumors. J Mol Diagn 7: 65–71
Dowsett M, Bartlett J, Ellis IO, Salter J, Hills M, Mallon E, Watters AD, Cooke T, Paish C, Wencyk PM, Pinder SE (2003) Correlation between immunohistochemistry (HercepTest) and fluorescence in situ hybridization (FISH) for HER-2 in 426 breast carcinomas from 37 centres. J Pathol 199: 418–423
Dressman H, Hans C, Bild A, Olson Jr JA (2006) Gene expression profiles of multiple breast cancer phenotypes and response to neoadjuvant chemotherapy. Clin Cancer Res 12: 819–826
Farabegoli F, Hermsen MA, Ceccarelli C, Santini D, Weiss MM, Meijer GA, van Diest PJ (2004) Simultaneous chromosome 1q gain and 16q loss is associated with steroid receptor presence and low proliferation in breast carcinoma. Mod Pathol 17: 449–455
Fazeny-Dorner B, Piribauer M, Wenzel C, Fakhrai N, Pirker C, Berger W, Sedivy R, Rudas M, Filipits M, Okamoto I, Marosi C (2003) Cytogenetic and comparative genomic hybridization findings in four cases of breast cancer after neoadjuvant chemotherapy. Cancer Genet Cytogenet 146: 161–166
Goldhirsch A, Glick JH, Gelber RD, Coates AS, Thurlimann B, Senn HJ (2005) Meeting highlights: international expert consensus on the primary therapy of early breast cancer 2005. Ann Oncol 16: 1569–1583
Han W, Han MR, Kang JJ, Bae JY, Lee JH, Bae YJ, Lee JE, Shin HJ, Hwang KT, Hwang SE, Kim SW, Noh DY (2006) Genomic alterations identified by array comparative genomic hybridization as prognostic markers in tamoxifen-treated estrogen receptor-positive breast cancer. BMC Cancer 6: 92
Hannemann J, Oosterkamp HM, Bosch CA, Velds A, Wessels LF, Loo C, Rutgers EJ, Rodenhuis S, van de Vijver MJ (2005) Changes in gene expression associated with response to neoadjuvant chemotherapy in breast cancer. J Clin Oncol 23: 3331–3342
Hayashita Y, Osada H, Tatematsu Y, Yamada H, Yanagisawa K, Tomida S, Yatabe Y, Kawahara K, Sekido Y, Takahashi T (2005) A polycistronic microRNA cluster, miR-17-92, is overexpressed in human lung cancers and enhances cell proliferation. Cancer Res 65: 9628–9632
Hupe P, Stransky N, Thiery JP, Radvanyi F, Barillot E (2004) Analysis of array CGH data: from signal ratio to gain and loss of DNA regions. Bioinformatics 20: 3413–3422
Ioannidis JP (2005) Microarrays and molecular research: noise discovery? Lancet 365: 454–455
Iwao-Koizumi K, Matoba R, Ueno N, Kim SJ, Ando A, Miyoshi Y, Maeda E, Noguchi S, Kato K (2005) Prediction of docetaxel response in human breast cancer by gene expression profiling. J Clin Oncol 23: 422–431
Kuwano M, Uchiumi T, Hayakawa H, Ono M, Wada M, Izumi H, Kohno K (2003) The basic and clinical implications of ABC transporters, Y-box-binding protein-1 (YB-1) and angiogenesis-related factors in human malignancies. Cancer Sci 94: 9–14
Leyland-Jones B, Kelland LR, Harrap KR, Hiorns LR (1999) Genomic imbalances associated with acquired resistance to platinum analogues. Am J Pathol 155: 77–84
Loo LW, Grove DI, Williams EM, Neal CL, Cousens LA, Schubert EL, Holcomb IN, Massa HF, Glogovac J, Li CI, Malone KE, Daling JR, Delrow JJ, Trask BJ, Hsu L, Porter PL (2004) Array comparative genomic hybridization analysis of genomic alterations in breast cancer subtypes. Cancer Res 64: 8541–8549
Mackay A, Dixon JM, Urruticoechea A, Dexter T, Iravani M, Fenwick K, Young O, White S, Miller WR, Evans DB, Ashworth A, Dowsett M (2005) Molecular determinants of aromatase inhibitor sensitivity in primary breast cancer. Breast Cancer Res Treat 94: S48 (abstract)
McCann AH, Miller N, O’Meara A, Pedersen I, Keogh K, Gorey T, Dervan PA (1996) Biallelic expression of the IGF2 gene in human breast disease. Hum Mol Genet 5: 1123–1127
Modlich O, Prisack HB, Munnes M, Audretsch W, Bojar H (2005) Predictors of primary breast cancers responsiveness to preoperative epirubicin/cyclophosphamide-based chemotherapy: translation of microarray data into clinically useful predictive signatures. J Transl Med 3: 32
Natrajan R, Williams RD, Hing SN, Mackay A, Reis-Filho JS, Fenwick K, Iravani M, Valgeirsson h, Grigoriadis A, Langford CF, Dovey O, Gregory SG, Weber BL, Ashworth A, Grundy PE, Pritchard-Jones K, Jones C (2006) Array CGH profiling of favourable histology wilms tumours reveals novel gains and losses associated with relapse. J Pathol 210: 49–58
Naylor TL, Greshock J, Wang Y, Colligon T, Yu QC, Clemmer V, Zaks TZ, Weber BL (2005) High resolution genomic analysis of sporadic breast cancer using array-based comparative genomic hybridization. Breast Cancer Res 7: R1186–R1198
Nessling M, Richter K, Schwaenen C, Roerig P, Wrobel G, Wessendorf S, Fritz B, Bentz M, Sinn HP, Radlwimmer B, Lichter P (2005) Candidate genes in breast cancer revealed by microarray-based comparative genomic hybridization of archived tissue. Cancer Res 65: 439–447
Nupponen NN, Porkka K, Kakkola L, Tanner M, Persson K, Borg A, Isola J, Visakorpi T (1999) Amplification and overexpression of p40 subunit of eukaryotic translation initiation factor 3 in breast and prostate cancer. Am J Pathol 154: 1777–1783
O’Connell P (2003) Genetic and cytogenetic analyses of breast cancer yield different perspectives of a complex disease. Breast Cancer Res Treat 78: 347–357
Pang E, Hu Y, Chan KY, Lai PB, Squire JA, Macgregor PF, Beheshti B, Albert M, Leung TW, Wong N (2005) Karyotypic imbalances and differential gene expressions in the acquired doxorubicin resistance of hepatocellular carcinoma cells. Lab Invest 85: 664–674
Paris PL, Andaya A, Fridlyand J, Jain AN, Weinberg V, Kowbel D, Brebner JH, Simko J, Watson JE, Volik S, Albertson DG, Pinkel D, Alers JC, van der Kwast TH, Vissers KJ, Schroder FH, Wildhagen MF, Febbo PG, Chinnaiyan AM, Pienta KJ, Carroll PR, Rubin MA, Collins C, van Dekken H (2004) Whole genome scanning identifies genotypes associated with recurrence and metastasis in prostate tumors. Hum Mol Genet 13: 1303–1313
Ransohoff DF (2004) Rules of evidence for cancer molecular-marker discovery and validation. Nat Rev Cancer 4: 309–314
Reis-Filho JS, Simpson PT, Gale T, Lakhani SR (2005a) The molecular genetics of breast cancer: the contribution of comparative genomic hybridization. Pathol Res Pract 201: 713–725
Reis-Filho JS, Simpson PT, Jones C, Steele D, Mackay A, Iravani M, Fenwick K, Valgeirsson H, Lambros M, Ashworth A, Palacios J, Schmitt F, Lakhani SR (2005b) Pleomorphic lobular carcinoma of the breast: role of comprehensive molecular pathology in characterization of an entity. J Pathol 207: 1–13
Reis-Filho JS, Simpson PT, Turner NC, Lambros MB, Jones C, Mackay A, Grigoriadis A, Sarrio D, Savage K, Dexter T, Iravani M, Fenwick K, Weber B, Hardisson D, Schmitt FC, Palacios J, Lakhani SR, Ashworth A (2006a) Comprehensive molecular genetic analysis reveals FGFR1 as a potential therapeutic target for lobular breast carcinomas. Clin Cancer Res 12: doi:10.1158/1078-0432.CCR-06-1164
Reis-Filho JS, Westbury C, Pierga JY (2006b) The impact of expression profiling on prognostic and predictive testing in breast cancer. J Clin Pathol 59: 225–231
Rennstam K, Ahlstedt-Soini M, Baldetorp B, Bendahl PO, Borg A, Karhu R, Tanner M, Tirkkonen M, Isola J (2003) Patterns of chromosomal imbalances defines subgroups of breast cancer with distinct clinical features and prognosis. A study of 305 tumors by comparative genomic hybridization. Cancer Res 63: 8861–8868
Rouzier R, Perou CM, Symmans WF, Ibrahim N, Cristofanilli M, Anderson K, Hess KR, Stec J, Ayers M, Wagner P, Morandi P, Fan C, Rabiul I, Ross JS, Hortobagyi GN, Pusztai L (2005) Breast cancer molecular subtypes respond differently to preoperative chemotherapy. Clin Cancer Res 11: 5678–5685
Roylance R, Gorman P, Harris W, Liebmann R, Barnes D, Hanby A, Sheer D (1999) Comparative genomic hybridization of breast tumors stratified by histological grade reveals new insights into the biological progression of breast cancer. Cancer Res 59: 1433–1436
Seute A, Sinn HP, Schlenk RF, Emig R, Wallwiener D, Grischke EM, Hohaus S, Dohner H, Haas R, Bentz M (2001) Clinical relevance of genomic aberrations in homogeneously treated high-risk stage II/III breast cancer patients. Int J Cancer 93: 80–84
Shelley Hwang E, Nyante SJ, Yi Chen Y, Moore D, DeVries S, Korkola JE, Esserman LJ, Waldman FM (2004) Clonality of lobular carcinoma in situ and synchronous invasive lobular carcinoma. Cancer 100: 2562–2572
Shimizu H, Fukuda T, Ghazizadeh M, Nagashima M, Kawanami O, Suzuki T (2002) Molecular cytogenetic characterization of drug-resistant leukemia cell lines by comparative genomic hybridization and fluorescence in situ hybridization. Jpn J Cancer Res 93: 902–910
Simpson PT, Reis-Filho JS, Gale T, Lakhani SR (2005) Molecular evolution of breast cancer. J Pathol 205: 248–254
Stange DE, Radlwimmer B, Schubert F, Traub F, Pich A, Toedt G, Mendrzyk F, Lehmann U, Eils R, Kreipe H, Lichter P (2006) High-resolution genomic profiling reveals association of chromosomal aberrations on 1q and 16p with histologic and genetic subgroups of invasive breast cancer. Clin Cancer Res 12: 345–352
Tirkkonen M, Tanner M, Karhu R, Kallioniemi A, Isola J, Kallioniemi OP (1998) Molecular cytogenetics of primary breast cancer by CGH. Genes Chromosomes Cancer 21: 177–184
Troester MA, Hoadley KA, Sorlie T, Herbert BS, Borresen-Dale AL, Lonning PE, Shay JW, Kaufmann WK, Perou CM (2004) Cell-type-specific responses to chemotherapeutics in breast cancer. Cancer Res 64: 4218–4226
Troyanskaya O, Cantor M, Sherlock G, Brown P, Hastie T, Tibshirani R, Botstein D, Altman RB (2001) Missing value estimation methods for DNA microarrays. Bioinformatics 17: 520–525
Turton NJ, Judah DJ, Riley J, Davies R, Lipson D, Styles JA, Smith AG, Gant TW (2001) Gene expression and amplification in breast carcinoma cells with intrinsic and acquired doxorubicin resistance. Oncogene 20: 1300–1306
van’t Veer LJ, Paik S, Hayes DF (2005) Gene expression profiling of breast cancer: a new tumor marker. J Clin Oncol 23: 1631–1635
Wu HK, Squire JA, Catzavelos CG, Weksberg R (1997) Relaxation of imprinting of human insulin-like growth factor II gene, IGF2, in sporadic breast carcinomas. Biochem Biophys Res Commun 235: 123–129
Yasui K, Mihara S, Zhao C, Okamoto H, Saito-Ohara F, Tomida A, Funato T, Yokomizo A, Naito S, Imoto I, Tsuruo T, Inazawa J (2004) Alteration in copy numbers of genes as a mechanism for acquired drug resistance. Cancer Res 64: 1403–1410
Author information
Authors and Affiliations
Corresponding author
Additional information
Supplementary Information accompanies the paper on British Journal of Cancer website (http://www.nature.com/bjc)
Rights and permissions
From twelve months after its original publication, this work is licensed under the Creative Commons Attribution-NonCommercial-Share Alike 3.0 Unported License. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc-sa/3.0/
About this article
Cite this article
Pierga, JY., Reis-Filho, J., Cleator, S. et al. Microarray-based comparative genomic hybridisation of breast cancer patients receiving neoadjuvant chemotherapy. Br J Cancer 96, 341–351 (2007). https://doi.org/10.1038/sj.bjc.6603483
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/sj.bjc.6603483
Keywords
This article is cited by
-
Integrative analysis of copy number and gene expression in breast cancer using formalin-fixed paraffin-embedded core biopsy tissue: a feasibility study
BMC Genomics (2017)
-
Annexin A6 is a scaffold for PKCα to promote EGFR inactivation
Oncogene (2013)
-
Genomic amplification and oncogenic properties of the GASC1 histone demethylase gene in breast cancer
Oncogene (2009)
-
A high-resolution integrated analysis of genetic and expression profiles of breast cancer cell lines
Breast Cancer Research and Treatment (2009)
-
Elevated chemokine receptor CXCR4 expression in primary tumors following neoadjuvant chemotherapy predicts poor outcomes for patients with locally advanced breast cancer (LABC)
Breast Cancer Research and Treatment (2009)