Multiparameter DNA content analysis identifies distinct groups in primary breast cancer.

BACKGROUND: Multiparameter flow cytometry is a robust and reliable method for determining tumour DNA content applicable to formalin-fixed paraffin-embedded (FFPE) tissue. This study examined the clinical and pathological associations of DNA content in primary breast cancer using an improved multiparametric technique.
METHODS: The FFPE tissue from 201 primary breast cancers was examined and the cancers categorised according to their DNA content using multiparametric flow cytometry incorporating differential labelling of stromal and tumour cell populations. Mathematical modelling software (ModFit 3.2.1) was used to calculate the DNA index (DI) and percentage S-phase fraction (SPF%) for each tumour. Independent associations with clinical and pathological parameters were sought using backward stepwise Binary Logistic Regression (BLR) and Cox's Regression (CR) analysis.
RESULTS: Tumours were grouped into four categories based on the DI of the tumour cell population. Low DI tumours (DI=0.76-1.14) associated with progesterone receptor-positive status (P=0.012, BLR), intermediate DI (DI=1.18-1.79) associated with p53 mutant tumours (P=0.001, BLR), high DI (DI1.80) tumours with human epidermal growth factor receptor 2 (HER2)-positive status (P=0.004, BLR) and 'multiploid tumours' (two or more tumour DNA peaks) did not show any significant associations. Tumours with high SPF% (10%) independently associated with poor overall survival (P=0.027, CR).
CONCLUSION: Multiparametric flow analysis of FFPE tissue can accurately assess tumour DNA content. Tumour sub-populations associated with biomarkers of prognosis or likely response to therapy. The alterations in DNA content present the potential for greater understanding of the mechanisms underlying clinically significant biomarker changes in primary breast cancer.

Chromosomal instability (CIN) is a hallmark of cancer, and for the common epithelial tumours complex and clonal karyotypic evolution is the predominant driver of clinical disease. Breast carcinomas show gain and/or loss of either whole chromosomes (aneuploidy) or segments of chromosomes (segmental aneuploidy), which frequently results in abnormal DNA content (Fridlyand ). Flow cytometry provides a reliable estimate of the DNA content of tumour populations and provides a global assessment of genetic instability. Flow cytometric and cytogenetic evidence suggests that most common epithelial cancers have a chromosomal constitution in the near diploid or sub-tetraploid range (Cornelisse ; Olaharski ; Weaver ). It is hypothesised that tumours proceed from near diploid to high structural and numerical chromosomal aneuploidy by a process of endoreduplication. This is followed by further non-random gain and/or loss of specific chromosomes or chromosomal regions resulting in an intermediate DNA index (DI) with enhanced malignant potential (Fujiwara ; Ganem ).

In breast cancer, previous flow cytometry studies have correlated DNA content with tumour grade (Hitchcock ), p53 mutation status (Eissa ; Friedrich ), disease recurrence/relapse (Ewers ; Stal ), and survival (Shiao ; Utada ), suggesting a therapeutic and prognostic clinical relevance. Tumour DNA content, in conjunction with S-phase fraction (SPF), provides a reliable measure of breast cancer aggression (Ferno ; Bagwell ). Using Cox's proportional hazard modelling, low DI, node negative, breast cancers were associated with good prognosis (Bagwell ). However, the technical limitations of the single parametric flow cytometry technique, including the inability to differentiate between DNA diploid tumour and normal cell populations (Hedley ), does not allow accurate sub-grouping of breast tumours, thus limiting further statistical analysis. Additionally, previous attempts to improve DNA content analysis of formalin-fixed paraffin-embedded (FFPE) samples using the single parameter technique by introducing an external control as a reference were unsuccessful, in particular for the detection of tumours with a DI close to 1.0 (Schutte ; Darzynkiewicz, 2010). Further, high DI tumours in the DNA tetraploid range with <25% tumour cells could not be differentiated from low DI tumours with an increased G2/M phase (Bauer, 1993). These limitations of single parametric flow have made DNA content analysis controversial and prevented the technique from achieving clinical utility (Ross ).

Dual parameter DNA content flow cytometry of FFPE tissues allows detection of epithelial cells by anti-keratin antibodies (Nylander ; Leers ). This was further improved by Corver , who used anti-keratin and anti-vimentin antibodies to allow specific identification of the tumour epithelial cells along with diploid stromal cells used as an internal reference, thus enabling accurate grouping of tumours based on their DI. This also enabled identification of DNA aneuploid tumour populations with a DI close to 1.0. It has been previously suggested that carcinoma-associated fibroblasts can be genetically abnormal and can have an abnormal DNA content; however, more recent studies have shown that, whatever the gene expression features of carcinoma-associated fibroblasts, the chromosome complement of stromal cells is diploid (Allinen ; Hosein ; Corver ).

Additionally, the multiparametric technique facilitated the identification of high DI tumour (epithelial) cells overlapping the G2/M phase of the diploid stromal cells allowing characterisation of these tumours. These advantages of the multiparametric technique make it a robust platform for DNA content analysis of FFPE tissues and give it the potential to be a clinically useful tool. The aims of this study were to confirm the application of multiparameter flow cytometry on FFPE breast cancer specimens and to examine the clinical and biomarker associations of multiparametric flow analysis in primary breast cancer.

Materials and methods

Patients and tissues

The FFPE tissue blocks from 201 primary, previously untreated operable breast cancers were selected at random from patients diagnosed between 1997 and 2002. The mean age of the patients was 61 years (age range 28–88 years). All patients provided consent; the project received ethical permission from Tayside Tissue Bank under delegated authority from the Local Tayside Research Ethics Committee.

Blocks were reviewed to confirm the presence of >40% tumour cells. Tumours were graded according to the NHS guidelines for pathology reporting of breast disease (Bloom and Richardson, 1957) by a specialist breast pathologist (CAP or LBJ) (Table 1). Human epidermal growth factor receptor 2 (HER2) status was determined using immunohistochemistry (IHC) with confirmation of HER2-positive tumours by fluorescent in situ hybridisation (FISH) (Hsi and Tubbs, 2004; Purdie ). The Quickscore (Detre ) method was employed to determine oestrogen receptor (ER) and progesterone receptor (PgR) status. P53 mutation status was obtained using the p53 Amplichip test (Roche, Pleasanton, CA, USA) that detects single base pair substitutions and deletions (Baker ). Breast cancer-specific overall survival, disease-free survival and recurrence were recorded with a minimum of 5 years follow-up for patients included in the study.

Table 1

Summary of the clinico-histopathological variables for 201 breast cancers studied using flow cytometry for tumour DNA content

	N	%
All	201
DNA content
Low	79	39.3
Intermediate	42	20.8
High	23	11.4
Mixed	57	28.3
Tumour grade
1	23	11.4
2	79	39.3
3	94	46.8
Missing	5	2.5
ER status
Negative	48	23.9
Positive	152	75.6
Missing	1	0.5
PgR status
Negative	89	44.3
Positive	111	55.2
Missing	1	0.5
HER2 status
Negative	170	84.6
Positive	28	13.9
Missing	3	1.5
p53 status
Wild type	146	72.6
Mutant	55	27.4
Survival status
Alive	163	81.1
Dead	38	18.9
S-phase fraction
Low	105	52.2
Medium	41	20.4
High	55	27.4

Abbreviations: ER, oestrogen receptor; PgR, progesterone receptor; HER2, human epidermal growth factor receptor 2.

Tissue processing and staining for flow cytometry

Tissue sections were processed as described by Corver . Briefly, 3 × 60 μm thick sections from each block were de-waxed using xylene and rehydrated with decreasing concentrations of alcohol before antigen retrieval using 10 mℳ citric acid buffer (pH 6) for 1 h at 80 °C. Tissue was dissociated with a solution of 0.1% collagenase (Sigma-Aldrich, Dorset, UK) and dispase (Gibco; Invitrogen, Paisley, UK) for 1 h. One million cells were incubated overnight with anti-cytokeratin (2 μg ml−1 MNF116, DAKO (Glostrup, Denmark) and 5 μg ml−1 AE1AE3, Chemicon International (Temecula, CA, USA) MAB3412) and anti-vimentin (diluted 1:5 V9-2b, Antibodies for Research Applications BV, Gouda, The Netherlands; http://www.abforresearch.nl/) primary antibodies. Cells were washed with PBA (500 ml PBS/5% BSA/0.05% Tween-20) and treated with fluorescently conjugated secondary antibodies, FITC-IgG1 (Southern Biotech, Birmingham, AL, USA) for cytokeratin detection and PE-IgG2b (Southern Biotech) for vimentin detection for 30 min. Finally, the cells were resuspended in a solution of 10 μℳ propidium iodide (PI) (Sigma-Aldrich) and 0.1% RNase (Sigma-Aldrich) for 2–3 h or overnight before analysis.

Flow cytometry and analysis

Cell suspensions were analysed using a FACScan (BD Biosciences, San Jose, CA, USA) with CellQuest software (version 3; BD Biosciences). DNA content analysis was performed by acquiring 50 000 events from each sample and simultaneously measuring PI, FITC-labelled epithelial, and PE-labelled stromal cell populations (Corver ; Figure 1C). The keratin-positive tumour epithelial cell population detected using flow cytometry ranged from 6% to 90% of the total events analysed. Doublet discrimination was applied to remove cell clumps and enable the analysis of single cell events by using PI area vs PI width signal and noise/debris in front of the first G1 peak was also gated out (Figure 1A) during acquisition. The 530/30 band pass filter was used to detect FITC-labelled epithelial cells in the FL1 channel, PE-labelled stromal cells were detected using the 585/42 band pass filter in the FL2 channel and finally, PI was detected using the 650 long pass filter in the FL3 channel. According to the published guidelines for DNA content analysis from FFPE tissues, samples with a coefficient of variation (CV) of >8% should be excluded from data analysis to ensure robust data generation (Hedley, 1993). All breast cancers included in the present study yielded a CV ranging between 3.05% and 7.84% and the reproducibility of the technique confirmed with 15 randomly selected samples (data not shown).

Figure 1

The four DNA content categories based on multi-parametric analysis. (A) Dot plot showing doublet discrimination to ensure analysis of a single cell population. (B) PI-positive cells as a negative control. (C) Dot plot showing keratin-positive epithelial cell population (EP), and vimentin-positive stromal cell population (SP). (D) DNA content histograms showing an example of a low DI, intermediate DI, high DI, and a multiploid tumour. The epithelial tumour cell population (green) is shown overlapping the normal stromal cell population (red) is shown. Histograms in (E) show the mathematical modelling analysis using the ModFit algorithm. The color reproduction of this figure is available at the British Journal of Cancer journal online.

The DI and percentage SPF (SPF%) were calculated using ModFit 3.2.1 and WinList 6.0 (Verity Software House, Topsham, ME, USA). The two software packages were remotely linked allowing ModFit to use the median of the diploid stromal cells to calculate the DI of the tumour (epithelial) cell population (mean of G0/G1 tumour epithelial peak/mean of G0/G1 reference stromal peak) (Hedley, 1993; Dayal ; Corver and ter Haar, 2011). A data file contained all events except those at the lower and higher ends of the linear scale (Figure 1A), to allow the ModFit algorithm to correct for debris, clumps, and aggregates. ModFit analysis also included a gate to select all keratin-positive events in the FSC vs keratin dot plot, thus including doublets, noise and debris, allowing ModFit to handle data correction without manual interference.

Statistical analysis

Flow cytometry data were compared with clinical and histopathological parameters including histological tumour grade, p53 mutation status, HER2 status, ER status, and PgR status. An in-house data analysis tool (INSPIRE), which includes the two tail Fisher's exact test (FET) was used for preliminary statistical analysis (Quinlan ). Disease recurrence, disease-free survival, and overall survival were studied using the Kaplan–Meier (KM) log-rank test. Significant results were further analysed to confirm independent associations using the backward stepwise Binary Logistic Regression (BLR) and Cox's Regression (CR) statistical analysis. The strength of the associations was assessed by calculating the β coefficient. The P-values for each association were used to discriminate independent associations between the variables analysed. For all analyses, the null hypothesis was rejected at α level of 5% (P<0.05).

Results

Categorising primary breast tumours

In all, 201 primary breast cancers were grouped into four categories based on the DI of the tumour cell population (Figures 1D and E). Tumours with a DI between 0.76 and1.14 were observed to be in the ‘low DI' range (n=79), ‘intermediate DI' tumours (n=42) had a DI ranging from 1.18 to 1.79, whereas ‘high DI' tumours (n=23) had a DI of ⩾1.80, tumours with two or more peaks from one or more of the above three categories were classed as ‘multiploid' tumours (n=57). DNA index was consistent and no variations were observed when DNA content analysis was repeated for 15 randomly selected tumours ensuring the reliability of the data generated and DI cutoffs were carefully determined based on published recommended standards (Hedley, 1993; Shankey ). Figure 2 summarises the DI for all 201 samples used in this study

Figure 2

A bimodal DNA index distribution curve of 201 breast cancers. The number of samples is represented on the y axis and the DNA index is indicated on the x axis.

Based on the technical precision of multiparameter flow analysis (Corver ), tumours in the low DI range could further be categorised into three groups, (i) group 1 (n=26), with 0.76⩾DI⩽0.95, (ii) group 2 (n=25), with 0.96⩾DI⩽1.05, and (iii) group 3 (n=28), with 1.06⩾DI⩽1.14 (Figure 3).

Figure 3

Low DI sub-categories. Histograms showing a group 1 low DI tumour with overlapping stromal (SP) and epithelial (EP) cell population and DI=0.90, group 2 low DI=0.97, and group 3 low DI tumour with DI=1.11.

The SPF% was also obtained for each tumour using ModFit analysis. Tumours were separated into three categories based on SPF% as described previously (Stal ), low SPF% <5% (n=105), medium SPF%=5–10% (n=41), and high SPF% ⩾10% (n=55) tumours.

Clinical associations

Univariate statistical analysis

Low DI tumours showed significant associations with good prognostic features, including HER2-negative status (P=0.03, FET), ER-positive status (P=0.003, FET), PgR-positive status (P=0.013, FET), p53 wild type (P=0.03, FET), tumour grade 1–2 (P<0.0001, FET), disease-free survival (P=0.0027, KM), and overall survival (P=0.0023, KM). Further, analysis of the low DI cases (into subgroups 1, 2, and 3) indicated that group 2 low DI tumours associated with better prognostic markers when compared with the group 1 low DI tumours. The group 1 low DI tumour category (n=26) associated with poor prognostic features such as poor Nottingham prognostic index (NPI) (Galea ) (P=0.03, FET) and tumour grade 3 (P=0.02, FET). Whereas, group 2 low DI tumours (n=25) associated with lymph node-negative status (P=0.02, FET). Group 3 tumours showed no significant associations.

Intermediate DI tumours showed an association with markers of poor prognosis, including ER-negative status (P=0.02, FET), p53 mutant status (P=0.0016, FET), tumour grade 3 (P=0.0014, FET), and disease recurrence (P=0.04, KM). High DI tumours were only associated with HER2-positive status (P=0.006, FET) and multiploid tumours showed no significant associations to any of the clinical parameters examined.

The above tumour categories were further assessed using clinical markers associated with the luminal equivalent (ER, PgR), HER2-positive and triple receptor-negative breast cancers (TNBC). Low DI tumours were significantly associated with clinical markers (ER+ve and/or PgR+ve) that are characteristic of the luminal type cancers (P=0.004, FET), whereas intermediate DI tumours associated with TNBC (P=0.007, FET). The molecular subtyping of these tumours based on microarray expression analysis was outside the scope of this study.

Tumours with high SPF (⩾10%) associated with ER-negative status (P=0.001, FET), PgR-negative status (P=0.001, FET), p53 mutant tumours and intermediate-high DI tumours (P=0.001, FET).

Multivariate statistical analysis

Significant associations obtained from univariate analysis were tested using CR (Table 2) and BLR analyses (Table 3) to determine independent associations. An independent association between tumour DNA content and overall survival was not confirmed by multivariate analysis (Table 3). However, low DI tumours independently associated with PgR-positive status (P=0.012, BLR), intermediate DI tumours associated with p53 mutant status (P=0.001, BLR), and high DI tumours independently associated with HER2-positive status (P=0.004, BLR) (Table 3; Figure 4). High SPF (⩾10%) independently associated (P=0.027, CR) with low patient survival (Table 2; Figure 4).

Table 2

Cox's regression (CR) analysis

Response variable	Predictor variables	Hazards ratio	95% Low	95% High	P-value
Death	PgR positive	0.32	0.149	0.686	0.003
	p53 mutant	2.61	1.266	5.383	0.009
	Lymph node positive	2.782	1.239	6.249	0.013
	High SPF (⩾10)	2.218	1.094	4.499	0.027

Abbreviations: PgR=progesterone receptor; SPF=S-phase fraction.

All variables were analysed by backwards step-wise Cox's proportional hazards regression model to find independent associations with survival. Independent associations were confirmed with PgR-positive status, p53 mutation, positive lymph nodes, and high S-phase fraction (SPF3).

Table 3

Backward stepwise Binary Logistic Regression (BLR) statistical analysis

Response variable	Predictor variables	β	P-value	OR (95% CI)
PgR	p53 wild type	0.72	0.030	2.05 (1.07–3.93)
	Low DI tumours	−0.67	0.030	0.51 (0.28–0.94)
p53 wild type	ER	1.69	<0.001	5.40 (2.62–11.13)
	Intermediate DI tumours	−1.02	0.009	0.36 (0.17–0.78)
ER	p53 wild type	1.68	<0.001	5.37 (2.38–12.14)
	HER2	1.20	0.013	3.34 (1.29–8.63)
HER2	ER	1.38	0.002	3.96 (1.65–9.46)
	High DI tumours	−1.55	0.003	0.21 (0.08–0.59)
Low DI tumours	PgR	−0.76	0.012	0.47 (0.26–0.84)
Intermediate DI tumours	p53 wild type	−1.17	0.001	0.31 (0.15–0.63)
High DI tumours	HER2	−1.42	0.004	0.24 (0.09–0.64)

Abbreviations: CI=confidence interval; DI=DNA index; ER=oestrogen receptor; HER2=human epidermal growth factor receptor 2; OR=odds ratio; PgR=progesterone receptor.

Variables were analysed by backwards step-wise BLR. All significant associations between the predictor and response variables are shown and non-significant associations are omitted for brevity and clarity.

Figure 4

Independent associations of DNA content categories with prognostic markers. Independent associations between DNA content categories and molecular markers associated with breast cancer prognosis using BLR and CR analysis. The sense of association (negative or positive) and the corresponding level of confidence (P-value) are also indicated in the figure. The direction of the arrows indicates the direction of association found by BLR and CR. DI, DNA index; SPF%, S-phase fraction.

Discussion

This study has, for the first time, applied the multiparametric flow technique to FFPE tissue from 201 primary breast cancers. The technique allows assessment of tumour DNA content by differential staining of keratin expressing epithelial tumour cells (FITC labelled) from vimentin expressing normal stromal cells (PE labelled) (Corver ). The ability of this technique to utilise routinely processed FFPE clinical material of mixed tumour and stromal tissues is beneficial compared with other methods of DNA analysis, such as deep sequencing, that requires a minimum of 50–70% tumour cells from frozen tissue and relies on other labour intensive techniques. Reassuringly, the DI values obtained in this study from FFPE sections show a similar distribution (Figure 2) to previously published data on fresh/frozen tissue material (Kimmig ), underpinning the reliability of this multiparameter flow technique for DNA content analysis of FFPE tissues.

The association of low DI tumours with good prognostic indicators in breast cancer confirms previous findings (Cunningham ; Spiethoff ). Furthermore, reliable subgrouping of low DI FFPE tumours based on the DI, while not possible using the single parametric technique, demonstrated an association of group 2 low DI tumours with lymph node-negative status (P=0.02, FET), not previously reported using FFPE tissue. This suggests that the technical deficiencies of single parameter flow cytometry may have masked clinically relevant associations between lymph-node status and tumour DNA content (Beerman ; Bergers ; Spiethoff ).

Tumours belonging to the intermediate DI, high DI, and multiploid categories (Figure 1D and E) could not be appropriately grouped in previous studies due to the technical limitations of the single parameter technique to differentiate between epithelial and stromal cell populations. As a result, multiploid tumours may have been incorrectly allocated to either the intermediate or high DI tumour groups. A reliable separation of intermediate DI and high DI tumours could only be achieved as exemplified here with the multiparametric flow method. The association of DNA content categories with different markers of prognosis (Figure 4) in the present study indicates that these tumours clinically behave differently. This highlights the clinical potential of multiparametric flow analysis of breast cancer to determine sub-groups of cancers within tumour categories previously identified by DNA content using single parametric flow analysis.

Biomarker associations

The PgR expression and tumour p53 mutation have been postulated as independent markers of prognosis in breast cancer (Overgaard ; Liu ) as confirmed here (Table 3; Figure 4). An association of low DI tumours with ER-positive and PgR-positive status has been previously reported (Leers and Nap, 2001); however, for the first time we demonstrate multivariate analysis showing an independent association of low DI tumours with PgR, but not ER, positive IHC.

Tumours with aneuploid DNA have been associated with poor prognostic features (Friedrich ; Leers and Nap, 2001). On further interrogation these tumours revealed two distinct subgroups, (i) intermediate DI tumours, which associated with poor prognostic features including p53 mutation and triple negative status; and (ii) high DI tumours associated with HER2 positivity but not with any other markers of prognosis. Multivariate statistical analysis confirmed these associations suggesting the presence of distinct molecular categories within the aneuploid tumours of potential clinical significance.

High SPF% independently associated with poor survival using CR analysis. Previous studies have failed to identify SPF% as an independent prognostic marker using single parameter DNA content analysis in lymph node-positive breast cancer patients (Witzig , 1993) but have found an association with overall survival in lymph node-negative tumours (Cunningham ) and, recently, in both early onset and advanced stage breast cancer (Vielh ). In the present study, SPF% was not found to be significantly associated with lymph-node status; however, both lymph-node status and SPF% were independently associated with survival (Table 2; Figure 4). The SPF% data presented here are inherently different from previously reported studies as the tumour (epithelial) cell fraction SPF% has been specifically calculated, not the total SPF%. Interestingly, high SPF% also associated with intermediate and high DI tumours (P=0.001, FET), the tumour categories in this study associated with poor prognostic markers.

The identification of clinically distinct DNA content sub-groups in the present study suggests that a reliable DNA content analysis platform is required and achievable to investigate the complex biological mechanisms that drive cancer. The SPF% and DNA content analysis using multiparameter flow analysis associate significantly with prognosis and warrant further investigation. Moreover, an understanding of the basic biology responsible for tetraploidy and sub-tetraploidy in breast cancer is likely to have direct relevance to similar events in other common epithelial tumours.

Conclusion

Multiparameter DNA content analysis of FFPE primary breast cancers identifies categories of breast cancer, based on DI, associated with conventional clinical and pathological parameters of prognostic and therapeutic significance. Investigating the mechanisms involved in tumour cell DNA content may provide future therapeutic opportunities in breast cancer.

Ethical approval

Ethical approval was obtained from the Tayside Tissue Bank under delegated authority from the Tayside Local Research Ethics Committee for all patient samples included in this study.