Search tips
Search criteria 


Logo of bmcgenoBioMed Centralsearchsubmit a manuscriptregisterthis articleBMC Genomics
BMC Genomics. 2006; 7: 96.
Published online Apr 27, 2006. doi:  10.1186/1471-2164-7-96
PMCID: PMC1468408
The molecular portraits of breast tumors are conserved across microarray platforms
Zhiyuan Hu,1,2 Cheng Fan,1 Daniel S Oh,1,2 JS Marron,3 Xiaping He,1,2 Bahjat F Qaqish,4 Chad Livasy,5 Lisa A Carey,6 Evangeline Reynolds,6 Lynn Dressler,6 Andrew Nobel,3 Joel Parker,7 Matthew G Ewend,6 Lynda R Sawyer,6 Junyuan Wu,1 Yudong Liu,1 Rita Nanda,8 Maria Tretiakova,8 Alejandra Ruiz Orrico,9 Donna Dreher,9 Juan P Palazzo,9 Laurent Perreard,10 Edward Nelson,11 Mary Mone,11 Heidi Hansen,11 Michael Mullins,12 John F Quackenbush,12 Matthew J Ellis,13 Olufunmilayo I Olopade,8 Philip S Bernard,12 and Charles M Peroucorresponding author1,2,5
1Lineberger Comprehensive Cancer Center, University of North Carolina, Chapel Hill, NC 27599, USA
2Department of Genetics, University of North Carolina, Chapel Hill, NC 27599, USA
3Department of Statistics and Operations Research, University of North Carolina, Chapel Hill, NC 27599, USA
4Department of Biostatistics, University of North Carolina, Chapel Hill, NC 27599, USA
5Department of Pathology and Laboratory Medicine, University of North Carolina, Chapel Hill, NC 27599, USA
6Department of Medicine, University of North Carolina, Chapel Hill, NC 27599, USA
7Constella Health Sciences, 2605 Meridian Parkway, Durham, NC 27713, USA
8Section of Hematology/Oncology, Department of Medicine, Committees on Genetics and Cancer Biology, University of Chicago, 5841 South Maryland Avenue, Chicago, IL 60637-1463, USA
9Department of Pathology, Thomas Jefferson University, 132 South 10th Street Philadelphia, PA 19107, USA
10The ARUP Institute for Clinical and Experimental Pathology, 500 Chipeta Way, Salt Lake City, Utah 84108, USA
11Department of Surgery, University of Utah School of Medicine, 30 N 1900 E, Salt Lake City, Utah 84132, USA
12Department of Pathology, University of Utah School of Medicine, 30 N 1900 E, Salt Lake City, Utah 84132, USA
13Department of Medicine, Division of Oncology, Washington University School of Medicine and Siteman Cancer Center, St Louis, Missouri, USA
corresponding authorCorresponding author.
Zhiyuan Hu: zhiyuan_hu/at/; Cheng Fan: cfan2004/at/; Daniel S Oh: daniel_oh/at/; JS Marron: marron/at/; Xiaping He: xiaping_he/at/; Bahjat F Qaqish: qaqish/at/; Chad Livasy: chad_livasy/at/; Lisa A Carey: Lisa_Carey/at/; Evangeline Reynolds: evangeline_reynolds/at/; Lynn Dressler: lynn_dressler/at/; Andrew Nobel: nobel/at/; Joel Parker: jparker/at/; Matthew G Ewend: ewend/at/; Lynda R Sawyer: lynda_sawyer/at/; Junyuan Wu: junyuan_wu/at/; Yudong Liu: ydliu/at/; Rita Nanda: rnanda/at/; Maria Tretiakova: mtretiakova/at/; Alejandra Ruiz Orrico: Alejandra.RuizOrrico/at/; Donna Dreher: donna.dreher/at/; Juan P Palazzo: juan.palazzo/at/; Laurent Perreard: laurent.perreard/at/; Edward Nelson: edward.nelson/at/; Mary Mone: mary.mone/at/; Heidi Hansen: heidi.hansen/at/; Michael Mullins: michael.mullins/at/; John F Quackenbush: jonquack/at/; Matthew J Ellis: mellis/at/; Olufunmilayo I Olopade: folopade/at/; Philip S Bernard: phil.bernard/at/; Charles M Perou: cperou/at/
Received February 15, 2006; Accepted April 27, 2006.
Validation of a novel gene expression signature in independent data sets is a critical step in the development of a clinically useful test for cancer patient risk-stratification. However, validation is often unconvincing because the size of the test set is typically small. To overcome this problem we used publicly available breast cancer gene expression data sets and a novel approach to data fusion, in order to validate a new breast tumor intrinsic list.
A 105-tumor training set containing 26 sample pairs was used to derive a new breast tumor intrinsic gene list. This intrinsic list contained 1300 genes and a proliferation signature that was not present in previous breast intrinsic gene sets. We tested this list as a survival predictor on a data set of 311 tumors compiled from three independent microarray studies that were fused into a single data set using Distance Weighted Discrimination. When the new intrinsic gene set was used to hierarchically cluster this combined test set, tumors were grouped into LumA, LumB, Basal-like, HER2+/ER-, and Normal Breast-like tumor subtypes that we demonstrated in previous datasets. These subtypes were associated with significant differences in Relapse-Free and Overall Survival. Multivariate Cox analysis of the combined test set showed that the intrinsic subtype classifications added significant prognostic information that was independent of standard clinical predictors. From the combined test set, we developed an objective and unchanging classifier based upon five intrinsic subtype mean expression profiles (i.e. centroids), which is designed for single sample predictions (SSP). The SSP approach was applied to two additional independent data sets and consistently predicted survival in both systemically treated and untreated patient groups.
This study validates the "breast tumor intrinsic" subtype classification as an objective means of tumor classification that should be translated into a clinical assay for further retrospective and prospective validation. In addition, our method of combining existing data sets can be used to robustly validate the potential clinical value of any new gene expression profile.
The classification of human tumors using microarray data has been an area of intense research, but it remains a daunting task to validate a new profile and generate a clinically useful test. Many different gene expression-based predictors have been developed for breast cancer [1-9], and two different gene expression predictors have reached the final step of prospective clinical trial testing [10,11]. Using cDNA microarrays, we previously identified five distinct subtypes of breast tumors arising from at least two distinct cell types (basal-like and luminal epithelial cells) [1-3]. This molecular taxonomy was based upon an "intrinsic" gene set, which was identified using a supervised analysis to select genes that showed little variance within repeated samplings of the same tumor, but which showed high variance across tumors [1]. We showed that an intrinsic gene set reflects the stable biological properties of tumors and typically identifies distinct tumor subtypes that have prognostic significance, even though no knowledge of outcome was used to derive this gene set [3,12-14].
A major challenge for microarray studies, especially those with clinical implications, is validation [15,16]. Due to the practical barriers of cost and access to large numbers of fresh frozen tumor samples with associated clinical information, very few microarray studies have analyzed enough samples to allow promising initial findings to be sufficiently validated to justify the major investment required for clinical testing. An efficient approach would be to use public gene expression data repositories as test sets; however, it has been difficult to compare and/or combine data sets from independent laboratories due to differences in sample preparation, experimental design, and microarray platforms. An accepted method for validation is to derive a prognostic/predictive gene set from a "training set" and then apply it to a completely independent "test set" [17]. The "purest" test sets are comprised of samples not generated by the primary investigators to remove any possibility of bias [18]. In this study, we illustrate the successful application of these principles by (1) deriving a new breast tumor intrinsic gene list that identifies the "intrinsic" biological features of breast tumors and (2) validating this predictor using a combined test set of 311 breast tumor samples compiled from the public domain. These analyses show that the breast tumor intrinsic subtypes are significant predictors of outcome when correcting for standard clinical parameters, and that common patterns of expression and outcome predictions can be identified in data sets generated by independent labs.
Identification of the Intrinsic/UNC gene set
Our goals were to (1) create a new breast tumor intrinsic list, (2) validate this list on an independent dataset to show the clinical significance of the "intrinsic" classifications, and (3) to derive an objective "intrinsic subtype" classifier that could be used clinically (see Figure Figure11 for overview of analyses performed). An intrinsic analysis is a "within class" versus "across classes" analysis that identifies genes that show low variability within a group (i.e. a tumor-metastasis pair), but which show high variation in expression across different tumors; in essence, one is selecting for genes that are consistently expressed when individual tumors are examined, but that vary in expression across different tumors. To develop a new breast tumor intrinsic gene set (Intrinsic/UNC), we assayed a training set of 105 breast tumor samples and 9 normal breast samples, which contained 26 sample pairs (See Additional file 2, 146 microarray experiments in total), using Agilent oligo microarrays. Using the intrinsic analysis method as described in Sorlie et al. 2003[3], we identified an intrinsic gene set of 1410 microarray elements representing 1300 genes. We felt it important to create a new intrinsic list because first, we wanted to take advantage of newer microarrays (Agilent arrays with 17,000 genes vs. 8,000 gene cDNA microarrays previously used[3]), and second, we wanted to use paired tumor samples that were not before-and-after chemotherapy pairs, but were instead pre-treatment tumor pairs. The Intrinsic/UNC gene set showed overlap with a previous breast tumor intrinsic gene set (108 genes in common with the Intrinsic/Stanford gene set of Sorlie et al. 2003[3]), but also showed a significant increase in gene number likely due to the greater number of genes present on current microarrays.
Figure 1
Figure 1
Overview of the analysis methods and datasets used in this paper.
Validation of the Intrinsic/UNC gene list
To evaluate the Intrinsic/UNC gene set on an independent test dataset, we applied it to a "combined test set" of 315 breast samples (311 tumors and 4 normal breast samples) using hierarchical clustering methods as have been done previously [1-3]. The "combined test set" of 315 breast samples was a single data set created by combining together the data from Sorlie et al. 2001 and 2003 (cDNA microarrays)[2,3], van't Veer et al. 2002 (custom Agilent oligo microarrays)[5] and Sotiriou et al.2003 (cDNA microarrays)[19]. We created a single data table of these three sets by first identifying the common genes present across all three microarray data sets (2800 genes). Next, we used Distance Weighted Discrimination (DWD) to combine these three data sets together [20]. DWD is a multivariate analysis tool that is able to identify systematic biases present in separate data sets and then make a global adjustment to compensate for these biases; in essence, each separate data set is a multi-dimensional cloud of data points, and DWD takes two points clouds and shifts one such that it more optimally overlaps the other. Finally, we determined that 306 of the 1300 unique Intrinsic/UNC genes were present in the combined test set and performed a hierarchical clustering analysis of these 306 genes and 315 samples (Figure (Figure2;2; see Additional file 1, for the complete cluster diagram). We analyzed the combined test set instead of analyzing each of the 3 datasets separately because we believed this would provide more statistical power to perform multivariate analysis, and would yield more meaningful results because any finding would need to be shared/present across all 3 datasets. Remarkably, despite the loss of genes in the Intrinsic/UNC list due to the requirement of having to be present on 4 different microarray platforms, the hierarchical clustering analysis in Figure Figure22 identified the five main subtypes/groups corresponding to the previously defined HER2+/ER-, Basal-like, LumA, LumB and Normal Breast-like tumor groups [2,3].
Figure 2
Figure 2
Hierarchical cluster analysis of the 315-sample combined test set using the Intrinsic/UNC gene set reduced to 306 genes. (A) Overview of complete cluster diagram. (B) Experimental sample-associated dendrogram. (C) Luminal/ER+ gene cluster with GATA3-regulated (more ...)
As shown in previous studies, a HER2+ expression cluster was observed in the cluster analysis of the "combined test set" and contained multiple genes from the 17q11 amplicon including HER2/ERBB2 and GRB7 (Figure (Figure2D).2D). The HER2+ intrinsic subtype (pink dendrogram branch in Figure Figure2B)2B) was predominantly ER-negative (i.e. HER2+/ER-) as previously shown. A Basal-like expression cluster was also present and contained genes (i.e. c-KIT, FOXC1 and P-Cadherin) previously identified to be characteristic of basal epithelial cells (Figure (Figure2F).2F). Using the program EASE[21], the Gene Ontology (GO) categories "extracellular space" and "extracellular region" were over-represented relative to chance in the Basal epithelial gene cluster. As shown in previous studies, a Luminal/ER+ expression cluster was present and contained ER, XBP1, FOXA1 and GATA3 (Figure (Figure2C).2C). GATA3 has recently been shown to be somatically mutated in some ER+ breast tumors, and some of the genes in Figure Figure2C2C are GATA3-regulated (FOXA1 and TFF3)[22], thus showing the functional clustering of a transcription factor and some of its direct targets. The Gene Ontology (GO) categories "transcription regulator activity" and "DNA binding" were over-represented relative to chance in the Luminal/ER+ gene cluster.
The most significant difference between the previous Intrinsic/Stanford gene lists and the new Intrinsic/UNC gene list was that the latter contained a large proliferation signature (Figure (Figure2G)2G) [23-25]. As expected, EASE analysis showed that the GO categories "mitotic cell cycle" and "M phase" were over-represented relative to chance in the proliferation signature. The inclusion of proliferation genes in the Intrinsic/UNC gene set, but not in the Intrinsic/Stanford gene set, is likely due to the fact that the Intrinsic/Stanford lists were based upon before-and-after chemotherapy paired samples of the same tumor, while the Intrinsic/UNC list was based upon paired samples taken at the same time point with respect to chemotherapy (22/26 were pre-treatment pairs). This finding suggests that tumor cell proliferation rates do vary before and after chemotherapy, but that proliferation is a reproducible and intrinsic feature of a tumor's expression profile.
A possible new tumor group (IFN) characterized by the high expression of Interferon (IFN)-regulated genes was observed in the combined test set analysis (Figure (Figure2E).2E). According to EASE, the GO categories "immune response" and "defense response" were over-represented relative to chance in the interferon-regulated gene cluster. This cluster contained STAT1, which is thought to be the transcription factor responsible for mediating IFN-regulation of gene expression [26,27]. Genes in the IFN cluster have been linked to lymph node metastasis and poor prognosis [7,13]. In summary, the Intrinsic/UNC list contained more genes than previous lists, encompasses most features of the Intrinsic/Stanford list (i.e. Basal, Luminal/ER+, and HER2-amplicon gene clusters) and adds the biologically and clinically relevant proliferation signature.
Tumor subtypes identified by the Intrinsic/UNC gene set are predictive of outcome
To determine how many biologically relevant tumor subtypes/groups might be present within the cluster in Figure Figure2,2, we used 3 criteria, which resulted in the identification of 6 potential subtypes/groups. The first criterion was the simple and obvious dendrogram branching pattern (Figure (Figure2B)2B) suggesting six groups. Second was the observation that each of the six groups uniquely expressed distinct sets of known biologically relevant genes including the basal, luminal/ER+, HER2-amplicon, IFN-regulated, and proliferation-associated signatures. Third was our knowledge of the previous classifications made by the Sorlie et al. 2003 Intrinsic/Stanford list of the Stanford/Norway samples (these samples are identified in Additional file 1): there was a high concordance (78%) between the classification of these samples made using either the Sorlie et al. 2003 Intrinsic/Stanford list or the Intrinsic/UNC list (excluding the IFN samples). Therefore, the 311 tumors/patients were stratified into six groups, and we proceeded to look for differences in outcomes and associations with other clinical parameters between these six groups. The Intrinsic/UNC gene set identified tumor groups/subtypes that were predictive of Relapse-Free Survival (RFS, Figure Figure3A)3A) and Overall Survival (OS, p = 0.000001, data not shown) in Kaplan-Meier survival analysis on the combined test set. As previously seen in Sorlie et al. (2001 and 2003), the LumA group had the best outcome while the HER2+/ER-, Basal-like, and LumB groups had significantly worse outcomes. The new IFN class had a Kaplan-Meier survival curve similar to that of LumB, and both showed elevated proliferation rates when compared to LumA (Figure (Figure2G2G).
Figure 3
Figure 3
Kaplan-Meier survival curves of breast tumors classified by intrinsic subtype. Survival curves are shown for (A) the 315-sample combined test set classified by hierarchical clustering using the Intrinsic/UNC gene set and (B) the 60-sample Ma et al., (more ...)
In the combined test set, the standard clinical parameters of ER status, node status, grade, and tumor size (note: data for clinical HER2 status was not available) were significant predictors of RFS using Kaplan-Meier analysis (Figure (Figure4),4), thus showing that the act of combining three different patient sets together did not destroy the prognostic abilities of these standard markers. In a multivariate Cox proportional hazards analysis of the combined test set using these standard clinical parameters, size, grade and ER status were significant predictors of RFS (Table (Table1A1A).
Figure 4
Figure 4
Kaplan-Meier survival curves using RFS as the endpoint, for the common clinical parameters present within the 315-sample combined test set. Survival curves are shown for (A) ER status, (B) node status, (C) histologic grade (1 = well-differentiated, 2 (more ...)
Table 1
Table 1
Multivariate Cox proportional hazards analysis of (A) standard clinical factors alone, or with (B) the Intrinsic Subtypes in relation to Relapse-Free Survival for the 315-sample combined test set. Size was a binary variable (0 = diameter of 2 cm or less, (more ...)
To further evaluate the prognostic/predictive value of the intrinsic subtype classification, we performed multivariate Cox proportional hazards analysis of the combined test set using the six intrinsic subtypes/groups defined above and the five standard clinical parameters with RFS, OS, or DSS as the endpoint (Table (Table1B1B shows analysis for RFS). The intrinsic subtypes, when added to the multivariate model containing the standard clinical variables, resulted in a model significantly more predictive of RFS, OS, and DSS (p = 0.01, 0.009, and 0.04 respectively, by the likelihood-ratio test). In multivariate analysis for RFS (Table (Table1B),1B), the Basal-like, LumB and HER2+/ER- subtypes had hazard ratios significantly greater than 1 (LumA served as the reference group), while the IFN and Normal Breast-like groups were not significant. Thus, the intrinsic subtypes classifications of LumA, LumB, Basal-like and HER2+/ER- add new and important prognostic information beyond what the standard clinical predictors provide.
Associations of the intrinsic subtypes with clinical and biological parameters
To further characterize and better understand the intrinsic subtypes, we determined whether an association existed between intrinsic subtype and grade, node status, ER status, age, and tumor size in the combined test set. Two-way contingency table analysis showed significant association between grade and subtype, with HER2+/ER- and Basal-like tumors more likely to be grade 3 (Table (Table2).2). The Cramer's V statistic[28], which measures the strength of association between two variables in a contingency table, indicated a substantial association (Cramer's V > 0.36) between grade and subtype. Two-way contingency table analysis did not show significant association between node status and subtype (p = 0.44), but did show significant association between ER status and subtype (p < 0.0001; Cramer's V = 0.72) and between tumor size and subtype (p = 0.01; Cramer's V = 0.17). As would be expected, ER+ tumors were more likely to be LumA or LumB. As indicated by the low Cramer's V (Cramer's V < 0.19 indicates a low relationship), tumor size and subtype were not strongly correlated.
Table 2
Table 2
Association between tumor histologic grade and intrinsic subtype in the 315-sample combined test set.
To determine association between age and subtype, we used an unpaired Student's t-test to compare the average ages of diagnosis of each tumor subtype. Interestingly, the average age of diagnosis for HER2+/ER- tumors was significantly less than that for all other tumor types. The average age of diagnosis for LumA tumors was significantly greater than that for LumB tumors.
Derivation and application of a Single Sample Predictor
A caveat to the above analyses is that our classifications were based upon hierarchical clustering, which is a powerful tool for intrinsic class discovery, but which is not suited for individual sample classification because to classify a new sample would require a reanalysis of all samples. Therefore, we wanted to create an unchanging and objective method to classify tumors according to intrinsic subtype that could be clinically applicable. To this end, we developed a Single Sample Predictor (SSP) using the combined test set hierarchically clustered using the 306 Intrinsic/UNC genes (Figure (Figure1).1). For the SSP, a mean expression profile (i.e. centroid) was created for each subtype that was significant in the multivariate analysis (LumA, LumB, Basal-like, HER2+/ER-) and for the Normal Breast-like group using the combined test set (Figure (Figure2).2). Next, any new sample is then compared to each Centroid and assigned by the SSP to the nearest subtype/centroid as determined by Spearman correlation (note: this SSP is based on methods developed by Tibshirani and colleagues[3,29,30]); thus, the SSP contains five different idealized profiles, and any new sample is compared to each of the five profiles and assigned a profile label (i.e. subtype name) based upon the single idealized profile it most resembled.
To validate the SSP, we tested it on two additional datasets not used previously. The first was the 60-patient Ma et al. dataset, which represents a group of early stage ER+ tamoxifen-treated patients [6]. The SSP classified these samples as follows: 2 Basal-like, 2 HER2+/ER-, 12 Normal Breast-like, 34 LumA, and 9 LumB. The 2 Basal-like and 2 HER2+/ER- assigned samples were excluded from a survival analysis because they were too few for a meaningful survival analysis and possibly were misclassified ER-negative tumors. Among the remaining samples the SSP classification was a significant predictor of RFS (p = 0.04, Figure Figure3B),3B), due to the poor outcome of the LumB group. Next, we applied the SSP to a 96-sample test set of local only (surgery)-treated patients from Chang et al. [31]. The tumor groups identified by the SSP showed significant differences in RFS (Figure (Figure3C,3C, p = 0.0006) and OS (p = 0.001, data not shown) in Kaplan-Meier analysis, with the poor outcome groups as expected: LumB, Basal-like, and HER2+/ER-. Thus, the SSP identified tumor groups that are truly prognostic and have significantly different outcomes as was seen before: namely, LumA always has the most favorable outcome, while LumB, Basal-like and HER2+/ER- do poorly[2,3,9,19].
We also applied the SSP onto the 105-sample dataset used to derive the Intrinsic/UNC gene list, which is technically not a test set for the SSP because it was used to derive the Intrinsic/UNC gene set. The tumor groups identified by the SSP showed significantly different RFS (Figure (Figure3D,3D, p = 0.02) and OS (p = 0.03, data not shown) in Kaplan-Meier analysis with the poor outcome groups again being LumB, Basal-like, and HER2+/ER-. A subset of the 105-sample dataset (48 in total) had been previously characterized using an immunohistochemical (IHC) analysis[32], which showed that (1) all 18 Basal-like tumors were ER-negative and HER2-negative (defined as not having a 3+ score on HER2 IHC analysis), (2) all 18 luminal subtype tumors were ER-positive and HER2-negative, and (3) all 12 HER2+/ER- subtype tumors were ER-negative and 11 out of these 12 showed HER2-overexpression (defined as having a 3+ score on HER2 IHC analysis). Thus, the SSP correlated with many standard clinical parameters, and was also able to identify clinically relevant groups (i.e. LumA vs. LumB) not identifiable using the standard clinical assays, thus indicating potential value as an objective classification method that should be developed further as a clinically applicable test.
The development and validation of gene sets for cancer patients requires significant resources because large training and test sets are required to achieve robust results. In fact, microarray studies are often criticized for a lack of rigorous validation due to small sample sizes [17,18]. Therefore, we utilized a previously described microarray data set combining method (Distance Weighted Discrimination) to create a large validation test set of over 300 tumors, and used it to validate a newly derived gene list for breast cancer prognostication and prediction. This approach allowed us to perform a multivariate analysis in which we show for the first time that the intrinsic subtype classification adds valuable information in the presence of five standard clinical parameters. We believe this combined test set is a valid test set for use in our analysis because after the multiple data sets were combined, the prognostic abilities of the standard clinical variables such as ER and grade remained intact.
The remarkable power of our DWD-based approach is indicated by the fact that although samples came from different platforms, hierarchical clustering analysis of the combined data set managed to group samples and genes based upon biology, and not some artifact caused by combining the data sets together. Evidence that this grouping reflected biology and not some artifact comes from (1) the finding that various Gene Ontology terms were significantly over-represented relative to chance in individual gene clusters seen in this analysis and (2) the groupings of the samples showed inter-dataset mixing and were significant predictors of outcome in univariate Kaplan-Meier and multivariate Cox analysis. It is also remarkable that this classification was successful in predicting outcome despite the fact that the Intrinsic/UNC gene set was reduced from 1300 genes to 306 genes in the combined test set; this indicates the robust nature of the intrinsic subtypes as defined by the new Intrinsic/UNC gene list.
One of the accomplishments of this manuscript was to develop an unchanging and objective intrinsic subtype predictor that could be used routinely in the clinical setting. This was accomplished by first identifying a new intrinsic gene set and then using this set to develop the Single Sample Predictor (SSP) that was shown here to be both prognostic on the local therapy-only patient subset from Chang et al. [31] and predictive of outcomes on the ER+ tamoxifen-treated data set of Ma et al. [6]. Many other gene expression based predictors for breast cancer patients have been developed, and in a complementary publication[33], we tested the intrinsic subtype SSP developed here, relative to those predictions made by four other previously published breast cancer prognostic/predictive gene sets using a single patient/tumor set of 295 cases; the four other expression-based predictors used were (1) the "70-gene" Good vs. Poor outcome predictor developed by van't Veer and colleagues[5,11], (2) the "Wound-Response" profile developed by Chang et al[31,34], (3) the "Recurrence Score (RS)" profile developed by Paik et al. [10], and (4) the 2-gene (HOXB13:IL17BR) ratio predictor developed by Ma et al. [6]. The results showed that of samples classified as Basal-like, HER2+/ER-, or LumB by the SSP, 93–100% were classified by the 70-gene, RS and Wound-Response predictors as being in each predictor's bad prognosis group. These data suggest that a high concordance exists across these multiple predictors, in particular the RS, 70-gene and Intrinsic Subtypes; thus, the new intrinsic gene list and classification method developed here, when compared to other predictors as accomplished in Fan et al. [33], showed that a high concordance across predictors exists, which provides additional validation for each predictor.
The results of this study advance our current knowledge of the intrinsic breast tumor subtypes and provides an objective method (SSP) for prospectively classifying tumors that could be used in the clinical setting. More broadly speaking, our findings show that while the individual brushstrokes (i.e. genes) may sometimes show discordance across data sets, the portraits created by the combined patterns of the individual brushstrokes is conserved and recognizable across datasets because of the similarities to the family portrait [24]. Moreover, these data show that the breast tumor intrinsic subtypes identified using the Intrinsic/UNC gene list can be generalized to many different patient sets, both treated and untreated.
Sample collection, RNA isolation and microarray hybridization
105 fresh frozen breast tumor samples and 9 normal breast samples were obtained using IRB-approved protocols at 4 institutions: the University of North Carolina at Chapel Hill (UNC-CH), The University of Utah, Thomas Jefferson University, and the University of Chicago. This sample set represents an ethnically and geographically diverse cohort. Additional file 2 contains clinical data for these samples. Patients were heterogeneously treated according to the standard of care dictated by disease stage, ER and HER2 status.
Total RNA was purified from each sample using the Qiagen RNAeasy Kit. RNA integrity was determined using the RNA 6000 Nano LabChip Kit and Agilent 2100 Bioanalyzer. Total RNA amplification and labeling were done as previously described in [35]. Microarray hybridizations were performed using Agilent Human oligonucleotide (1Av1, 1Av2 and custom designed 1Av1-based) microarrays using 2 μg of Cy3-labeled common reference sample that is a modified version of the Stratagene Human Universal Reference[36], and 2 μg of Cy5-labeled experimental sample. Microarrays were hybridized overnight, washed, dried, and scanned as described in [35]. The image files were analyzed with GenePix Pro 4.1 and loaded into the UNC-CH Microarray Database[37] where a Lowess normalization procedure was performed to adjust the Cy3 and Cy5 channels[38]. All primary microarray data associated with this study are available at [39], 2006 #2192} and in the GEO[40] under the accession number of GSE1992, series GSM34424-GSM34568.
Identification of the intrinsic gene set
We derived a new breast tumor intrinsic gene set, referred to as the "Intrinsic/UNC" list, using a training set composed of the 105 tumor samples described above, 9 normal breast samples, and 26 sample pairs (in total, represented by 146 microarrays). 15, 9, and 2 of the 26 sample pairs were different physical pieces of the same tumor (taken at the same time point), tumor-metastasis pairs and normal sample pairs, respectively. The background subtracted, Lowess normalized log2 ratio of Cy5 to Cy3 intensity values were first filtered to select genes that had a signal intensity of at least 30 units above background in both the Cy5 and Cy3 channels. Only genes that met these criteria in at least 70% of the 146 microarrays were included for subsequent analysis. Next, we performed an "intrinsic" analysis as described previously[3] using the 26 sample pairs and 86 additional microarrays. An intrinsic analysis identifies genes showing low variability in expression within paired samples but high variability in expression across different tumors; for each gene a ratio of "within-pair variance" to "between-subject variance" is computed. Genes with ratios below one standard deviation of the mean ratio were defined as "intrinsic". This analysis resulted in 1410 microarray elements representing 1300 genes being identified as "intrinsic". In order to obtain an estimate of the number of false-positive intrinsic genes, we permuted the sample labels to generate 26 random pairs and 86 non-paired samples. This permutation was performed 100 times and the intrinsic scores were calculated for each. These permuted scores were used to determine a threshold on the intrinsic score corresponding to a false discovery rate (FDR) less than 1%. The selected threshold resulted in 1410 microarray features being called significant with a median FDR = 0.3% and 90th percentile FDR = 0.5%. (See Tusher et al. for a complete description of this calculation [41]).
Creation and analyses of the combined test set
The independent test set was a 315-sample "combined test set" consisting of three DNA microarray datasets (Sorlie et al. 2001 and 2003[2,3], van't Veer et al. 2002[5] and Sotiriou et al. 2003[19]). To combine these datasets obtained from different microarray platforms, we performed the following pre-processing methods. First, the R/G ratios in each dataset were log2 transformed and Lowess normalized[38]. Next, missing values were k-NN imputed[42]. Gene annotations from each dataset were converted into UniGene Cluster IDs (UCIDs, Build 161) using the SOURCE database[43], and multiple occurrences of a UCID were collapsed by taking the median value for that ID within each experiment and platform, which resulted in ~2800 genes having expression data in all three datasets. Next, Distance Weighted Discrimination[20] was performed in a pair-wise fashion by first combining the Sorlie et al. and Sotiriou et al. datasets, and then combining this with the van't Veer et al. dataset to make a single dataset. In the final pre-processing step, each microarray experiment was normalized such that each column/experimental sample was standardized to N(0,1), and each row/gene was median centered. 306 of the 1300 Intrinsic/UNC genes had microarray data present in the combined test set and were used in a two-way average-linkage hierarchical cluster analysis [44]. Cluster results were visualized using the program "Treeview".
Derivation of the Single Sample Predictor
The Single Sample Predictor (SSP) is a Nearest Centroid-based method based upon the work of Hastie and Tibshirani [3,45,46]. Our SSP classifies an individual sample according to its nearest centroid as determined by Spearman correlation. To derive our SSP, we utilized the 315-sample combined test set from FigureFigure22 to create centroids for each of the five intrinsic subtypes (LumA, LumB, HER2+/ER-, Basal-like and Normal Breast-like). Please note that we did not create a centroid for the IFN group because it failed significance in multivariate testing, but did create a centroid for the Normal Breast-like group because we feel it is important to be able to identify true normal samples; an H&E examination of most tumor samples falling into the Normal Breast-like category shows that this is occurring mainly because of too much normal tissue contamination.
To create each intrinsic subtype centroid, we averaged the gene expression profiles for samples clearly assigned to each subtype (limiting the analysis to 249 of the 315 samples) using the hierarchical clustering dendrogram as a guide (Figure (Figure2).2). We then applied the SSP to two independent test datasets: (1) the Ma et al. 60-sample ER+ tamoxifen-treated tumor dataset and (2) the Chang et al. 96-sample local only-treated tumor dataset. By matching UCIDs, microarray data for as many as possible of the 306 Intrinsic/UNC genes was obtained from these 2 datasets. To remove microarray platform/source systematic biases, we applied DWD to the 2 test datasets relative to the combined test set. The SSP was then used to classify tumors by intrinsic subtype in these 2 test datasets. Using similar methods, the SSP was also applied to the 105-sample training set used to derive the intrinsic/UNC gene set.
Survival analyses
Kaplan-Meier survival plots were compared using the Cox-Mantel log-rank test in WinSTAT for Excel (R. Fitch Software). Two-way contingency table analysis and unpaired Student's t-test were done using WinSTAT. For the "combined test set", multivariate Cox proportional hazards analysis was performed using SAS (Cary, NC).
C.M.P. was the Principal Investigator and instigated the study, helped with design and wrote the paper, while J.P., O.I.O, and P.S.B. were the Principal Investigators at each of the three other participating institution and were involved in the study design, implementation and manuscript writing. C.F., J.S.M., B.F.Q., A.N., and J.P. were responsible for the statistical analyses and some writing. Z.H and X.H. performed all of the tumor RNA preparation and microarray experiments and were involved in the writing. J.W. and Y.L. were responsible for all data management and some data analysis. C.L. was responsible for the pathological assessment of most tumor samples and was involved in the writing. M.E. and D.O. were involved in data analysis, interpretation and writing. Tumor sample collection, clinical data acquisition and interpretation was accomplished by L.C., M.E., R.N., M.T., A.R.O., D.D., L.P., E.N., M.M., H.H., M.M., J.F.Q., L.R.S., E.R., and L.D., and it should be noted that this was separately accomplished at four institutions.
Supplementary Material
Additional File 2
Supplemental Table 1. Clinical and microarray information associated with each patient in the 105-sample training dataset.
Additional File 1
Supplemental Figure 1. Complete hierarchical cluster diagram of the 315-sample combined test set analyzed using the Intrinsic/UNC gene set, which was reduced to 306 genes based upon the gene overlap between datasets. Sorlie et al. sample names begin with the letters "BC", Sotiriou et al. sample names begin with "Exp", and van't Veer et al. sample names begin with "sample".
C.M.P. was supported by funds from the NCI Breast SPORE program to UNC-CH (P50-CA58223-09A1), by the National Institute of Environmental Health Sciences (U19-ES11391-03) and by NCI (RO1-CA-101227-01). P.S.B. was supported by NCI R33-CA97769-01, O.I.O. by the National Institute of Environmental Health Sciences (P50 ESO12382), A.N by NSF Grant DMS 0406361, J.S.M. by NSF Grant DMS-0308331, and L.A.C. by NIH M01RR00046. The Breast Cancer Research Foundation supported C.M.P., O.I.O and L.R.S.
  • Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, et al. Molecular portraits of human breast tumours. Nature. 2000;406:747–752. doi: 10.1038/35021093. [PubMed] [Cross Ref]
  • Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, van de Rijn M, Jeffrey SS, et al. Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci U S A. 2001;98:10869–10874. doi: 10.1073/pnas.191367098. [PubMed] [Cross Ref]
  • Sorlie T, Tibshirani R, Parker J, Hastie T, Marron JS, Nobel A, Deng S, Johnsen H, Pesich R, Geisler S, et al. Repeated observation of breast tumor subtypes in independent gene expression data sets. Proc Natl Acad Sci U S A. 2003;100:8418–8423. doi: 10.1073/pnas.0932692100. [PubMed] [Cross Ref]
  • Sotiriou C, Powles TJ, Dowsett M, Jazaeri AA, Feldman AL, Assersohn L, Gadisetti C, Libutti SK, Liu ET. Gene expression profiles derived from fine needle aspiration correlate with response to systemic chemotherapy in breast cancer. Breast Cancer Res. 2002;4:R3. doi: 10.1186/bcr433. [PMC free article] [PubMed] [Cross Ref]
  • van 't Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, Mao M, Peterse HL, van der Kooy K, Marton MJ, Witteveen AT, et al. Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002;415:530–536. doi: 10.1038/415530a. [PubMed] [Cross Ref]
  • Ma XJ, Wang Z, Ryan PD, Isakoff SJ, Barmettler A, Fuller A, Muir B, Mohapatra G, Salunga R, Tuggle JT, et al. A two-gene expression ratio predicts clinical outcome in breast cancer patients treated with tamoxifen. Cancer Cell. 2004;5:607–616. doi: 10.1016/j.ccr.2004.05.015. [PubMed] [Cross Ref]
  • Huang E, Cheng SH, Dressman H, Pittman J, Tsou MH, Horng CF, Bild A, Iversen ES, Liao M, Chen CM, et al. Gene expression predictors of breast cancer outcomes. Lancet. 2003;361:1590–1596. doi: 10.1016/S0140-6736(03)13308-9. [PubMed] [Cross Ref]
  • Zhao H, Langerod A, Ji Y, Nowels KW, Nesland JM, Tibshirani R, Bukholm IK, Karesen R, Botstein D, Borresen-Dale AL, et al. Different gene expression patterns in invasive lobular and ductal carcinomas of the breast. Mol Biol Cell. 2004;15:2523–2536. doi: 10.1091/mbc.E03-11-0786. [PMC free article] [PubMed] [Cross Ref]
  • Bertucci F, Finetti P, Rougemont J, Charafe-Jauffret E, Cervera N, Tarpin C, Nguyen C, Xerri L, Houlgatte R, Jacquemier J, et al. Gene expression profiling identifies molecular subtypes of inflammatory breast cancer. Cancer Res. 2005;65:2170–2178. doi: 10.1158/0008-5472.CAN-04-4115. [PubMed] [Cross Ref]
  • Paik S, Shak S, Tang G, Kim C, Baker J, Cronin M, Baehner FL, Walker MG, Watson D, Park T, et al. A multigene assay to predict recurrence of tamoxifen-treated, node-negative breast cancer. N Engl J Med. 2004;351:2817–2826. doi: 10.1056/NEJMoa041588. [PubMed] [Cross Ref]
  • van de Vijver MJ, He YD, van't Veer LJ, Dai H, Hart AA, Voskuil DW, Schreiber GJ, Peterse JL, Roberts C, Marton MJ, et al. A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med. 2002;347:1999–2009. doi: 10.1056/NEJMoa021967. [PubMed] [Cross Ref]
  • Bhattacharjee A, Richards WG, Staunton J, Li C, Monti S, Vasa P, Ladd C, Beheshti J, Bueno R, Gillette M, et al. Classification of human lung carcinomas by mRNA expression profiling reveals distinct adenocarcinoma subclasses. Proc Natl Acad Sci U S A. 2001;98:13790–13795. doi: 10.1073/pnas.191502998. [PubMed] [Cross Ref]
  • Chung CH, Parker JS, Karaca G, Wu J, Funkhouser WK, Moore D, Butterfoss D, Xiang D, Zanation A, Yin X, et al. Molecular classification of head and neck squamous cell carcinomas using patterns of gene expression. Cancer Cell. 2004;5:489–500. doi: 10.1016/S1535-6108(04)00112-6. [PubMed] [Cross Ref]
  • Garber ME, Troyanskaya OG, Schluens K, Petersen S, Thaesler Z, Pacyna-Gengelbach M, van de Rijn M, Rosen GD, Perou CM, Whyte RI, et al. Diversity of gene expression in adenocarcinoma of the lung. Proc Natl Acad Sci U S A. 2001;98:13784–13789. doi: 10.1073/pnas.241500798. [PubMed] [Cross Ref]
  • Michiels S, Koscielny S, Hill C. Prediction of cancer outcome with microarrays: a multiple random validation strategy. Lancet. 2005;365:488–492. doi: 10.1016/S0140-6736(05)17866-0. [PubMed] [Cross Ref]
  • Jenssen TK, Hovig E. Gene-expression profiling in breast cancer. Lancet. 2005;365:634–635. [PubMed]
  • Simon R, Radmacher MD, Dobbin K, McShane LM. Pitfalls in the use of DNA microarray data for diagnostic and prognostic classification. J Natl Cancer Inst. 2003;95:14–18. [PubMed]
  • Ioannidis JP. Microarrays and molecular research: noise discovery? Lancet. 2005;365:454–455. [PubMed]
  • Sotiriou C, Neo SY, McShane LM, Korn EL, Long PM, Jazaeri A, Martiat P, Fox SB, Harris AL, Liu ET. Breast cancer classification and prognosis based on gene expression profiles from a population-based study. Proc Natl Acad Sci U S A. 2003;100:10393–10398. doi: 10.1073/pnas.1732912100. [PubMed] [Cross Ref]
  • Benito M, Parker J, Du Q, Wu J, Xiang D, Perou CM, Marron JS. Adjustment of systematic microarray data biases. Bioinformatics. 2004;20:105–114. doi: 10.1093/bioinformatics/btg385. [PubMed] [Cross Ref]
  • Hosack DA, Dennis G, Jr, Sherman BT, Lane HC, Lempicki RA. Identifying biological themes within lists of genes with EASE. Genome Biol. 2003;4:R70. doi: 10.1186/gb-2003-4-10-r70. [PMC free article] [PubMed] [Cross Ref]
  • Usary J, Llaca V, Karaca G, Presswala S, Karaca M, He X, Langerod A, Karesen R, Oh DS, Dressler LG, et al. Mutation of GATA3 in human breast tumors. Oncogene. 2004;23:7669–7678. doi: 10.1038/sj.onc.1207966. [PubMed] [Cross Ref]
  • Perou CM, Jeffrey SS, van de Rijn M, Rees CA, Eisen MB, Ross DT, Pergamenschikov A, Williams CF, Zhu SX, Lee JC, et al. Distinctive gene expression patterns in human mammary epithelial cells and breast cancers. Proc Natl Acad Sci U S A. 1999;96:9212–9217. doi: 10.1073/pnas.96.16.9212. [PubMed] [Cross Ref]
  • Chung CH, Bernard PS, Perou CM. Molecular portraits and the family tree of cancer. Nat Genet. 2002;32:533–540. doi: 10.1038/ng1038. [PubMed] [Cross Ref]
  • Whitfield ML, Sherlock G, Saldanha AJ, Murray JI, Ball CA, Alexander KE, Matese JC, Perou CM, Hurt MM, Brown PO, et al. Identification of genes periodically expressed in the human cell cycle and their expression in tumors. Mol Biol Cell. 2002;13:1977–2000. doi: 10.1091/mbc.02-02-0030.. [PMC free article] [PubMed] [Cross Ref]
  • Bromberg JF, Horvath CM, Wen Z, Schreiber RD, Darnell JE., Jr Transcriptionally active Stat1 is required for the antiproliferative effects of both interferon alpha and interferon gamma. Proc Natl Acad Sci U S A. 1996;93:7673–7678. doi: 10.1073/pnas.93.15.7673. [PubMed] [Cross Ref]
  • Matikainen S, Sareneva T, Ronni T, Lehtonen A, Koskinen PJ, Julkunen I. Interferon-alpha activates multiple STAT proteins and upregulates proliferation-associated IL-2Ralpha, c-myc, and pim-1 genes in human T cells. Blood. 1999;93:1980–1991. [PubMed]
  • Van Belle G, Fisher L. Biostatistics: a methodology for the health sciences. 2. Hoboken, NJ: Wiley-Interscience John Wiley & Sons; 2004.
  • Bullinger L, Dohner K, Bair E, Frohling S, Schlenk RF, Tibshirani R, Dohner H, Pollack JR. Use of gene-expression profiling to identify prognostic subclasses in adult acute myeloid leukemia. N Engl J Med. 2004;350:1605–1616. doi: 10.1056/NEJMoa031046. [PubMed] [Cross Ref]
  • Bair E, Tibshirani R. Semi-supervised methods to predict patient survival from gene expression data. PLoS Biol. 2004;2:E108. doi: 10.1371/journal.pbio.0020108. [PMC free article] [PubMed] [Cross Ref]
  • Chang HY, Nuyten DS, Sneddon JB, Hastie T, Tibshirani R, Sorlie T, Dai H, He YD, van't Veer LJ, Bartelink H, et al. Robustness, scalability, and integration of a wound-response gene expression signature in predicting breast cancer survival. Proc Natl Acad Sci U S A. 2005;102:3738–3743. doi: 10.1073/pnas.0409462102. [PubMed] [Cross Ref]
  • Livasy CA, Karaca G, Nanda R, Tretiakova MS, Olopade OI, Moore DT, Perou CM. Phenotypic evaluation of the basal-like subtype of invasive breast carcinoma. Mod Pathol. 2005 [PubMed]
  • Fan C, Oh DS, Wessels L, Weigelt B, Nuyten DSA, Nobel AB, van't Veer LJ, Perou CM. Different gene expression-based predictors for breast cancer patients are concordant. N Engl J Med.
  • Chang HY, Sneddon JB, Alizadeh AA, Sood R, West RB, Montgomery K, Chi JT, van de Rijn M, Botstein D, Brown PO. Gene expression signature of fibroblast serum response predicts human cancer progression: similarities between tumors and wounds. PLoS Biol. 2004;2:E7. doi: 10.1371/journal.pbio.0020007. [PMC free article] [PubMed] [Cross Ref]
  • Hu Z, Troester M, Perou CM. High reproducibility using sodium hydroxide-stripped long oligonucleotide DNA microarrays. Biotechniques. 2005;38:121–124. [PubMed]
  • Novoradovskaya N, Whitfield ML, Basehore LS, Novoradovsky A, Pesich R, Usary J, Karaca M, Wong WK, Aprelikova O, Fero M, et al. Universal Reference RNA as a standard for microarray experiments. BMC Genomics. 2004;5:20. doi: 10.1186/1471-2164-5-20. [PMC free article] [PubMed] [Cross Ref]
  • UNC Microarray Database
  • Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP. Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002;30:e15. doi: 10.1093/nar/30.4.e15. [PMC free article] [PubMed] [Cross Ref]
  • UNC Breast Tumor Data
  • Gene Expression Omnibus
  • Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001;98:5116–5121. doi: 10.1073/pnas.091062498. [PubMed] [Cross Ref]
  • Troyanskaya O, Cantor M, Sherlock G, Brown P, Hastie T, Tibshirani R, Botstein D, Altman RB. Missing value estimation methods for DNA microarrays. Bioinformatics. 2001;17:520–525. doi: 10.1093/bioinformatics/17.6.520. [PubMed] [Cross Ref]
  • Diehn M, Sherlock G, Binkley G, Jin H, Matese JC, Hernandez-Boussard T, Rees CA, Cherry JM, Botstein D, Brown PO, et al. SOURCE: a unified genomic resource of functional annotations, ontologies, and gene expression data. Nucleic Acids Res. 2003;31:219–223. doi: 10.1093/nar/gkg014. [PMC free article] [PubMed] [Cross Ref]
  • Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998;95:14863–14868. doi: 10.1073/pnas.95.25.14863. [PubMed] [Cross Ref]
  • Hastie T, Tibshirani R, Friedman JH. The elements of statistical learning: data mining, inference, and prediction. New York: Springer; 2001.
  • Tibshirani R, Hastie T, Narasimhan B, Chu G. Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc Natl Acad Sci U S A. 2002;99:6567–6572. doi: 10.1073/pnas.082099299. [PubMed] [Cross Ref]
Articles from BMC Genomics are provided here courtesy of
BioMed Central