The molecular portraits of breast tumors are conserved across microarray platforms
- Zhiyuan Hu1, 2,
- Cheng Fan1,
- Daniel S Oh1, 2,
- JS Marron3,
- Xiaping He1, 2,
- Bahjat F Qaqish4,
- Chad Livasy5,
- Lisa A Carey6,
- Evangeline Reynolds6,
- Lynn Dressler6,
- Andrew Nobel3,
- Joel Parker7,
- Matthew G Ewend6,
- Lynda R Sawyer6,
- Junyuan Wu1,
- Yudong Liu1,
- Rita Nanda8,
- Maria Tretiakova8,
- Alejandra Ruiz Orrico9,
- Donna Dreher9,
- Juan P Palazzo9,
- Laurent Perreard10,
- Edward Nelson11,
- Mary Mone11,
- Heidi Hansen11,
- Michael Mullins12,
- John F Quackenbush12,
- Matthew J Ellis13,
- Olufunmilayo I Olopade8,
- Philip S Bernard12 and
- Charles M Perou1, 2, 5Email author
© Hu et al; licensee BioMed Central Ltd. 2006
Received: 15 February 2006
Accepted: 27 April 2006
Published: 27 April 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 . 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" . The "purest" test sets are comprised of samples not generated by the primary investigators to remove any possibility of bias . 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
Validation of the Intrinsic/UNC gene list
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 2D). The HER2+ intrinsic subtype (pink dendrogram branch in Figure 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 2F). Using the program EASE, 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 2C). GATA3 has recently been shown to be somatically mutated in some ER+ breast tumors, and some of the genes in Figure 2C are GATA3-regulated (FOXA1 and TFF3), 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 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 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
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, 1 = greater than 2 cm); node status was a binary variable (0 = no positive nodes, 1 = one or more positive nodes); age was a continuous variable formatted as decade-years. Hazard ratios for Intrinsic Subtypes were calculated relative to the Luminal A subtype. Variables found to be significant (p < 0.05) in the Cox proportional hazards model are shown in bold.
Hazard Ratio (95% CI)
Age, per decade
Tumor grade 2 vs. 1
Tumor grade 3 vs. 1
Hazard Ratio (95% CI)
Age, per decade
Tumor grade 2 vs. 1
Tumor grade 3 vs. 1
Basal-like vs. LumA
HER2+/ER- vs. LumA
LumB vs. LumA
IFN vs. LumA
Normal-like vs. LumA
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 1B 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 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
Association between tumor histologic grade and intrinsic subtype in the 315-sample combined test set.
Two-way contingency table
Statistics for two-way contingency table analysis
Cramer's V†† 0.42
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 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 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 . 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 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. . The tumor groups identified by the SSP showed significant differences in RFS (Figure 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 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, 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.  and predictive of outcomes on the ER+ tamoxifen-treated data set of Ma et al. . Many other gene expression based predictors for breast cancer patients have been developed, and in a complementary publication, 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. , and (4) the 2-gene (HOXB13:IL17BR) ratio predictor developed by Ma et al. . 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. , 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 . 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 . 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, and 2 μg of Cy5-labeled experimental sample. Microarrays were hybridized overnight, washed, dried, and scanned as described in . The image files were analyzed with GenePix Pro 4.1 and loaded into the UNC-CH Microarray Database where a Lowess normalization procedure was performed to adjust the Cy3 and Cy5 channels. All primary microarray data associated with this study are available at , 2006 #2192} and in the GEO 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 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 ).
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 and Sotiriou et al. 2003). 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. Next, missing values were k-NN imputed. Gene annotations from each dataset were converted into UniGene Cluster IDs (UCIDs, Build 161) using the SOURCE database, 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 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 . 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 Figure2 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 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.
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 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: Molecular portraits of human breast tumours. Nature. 2000, 406 (6797): 747-752. 10.1038/35021093.PubMedView Article
- Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, van de Rijn M, Jeffrey SS: Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci U S A. 2001, 98 (19): 10869-10874. 10.1073/pnas.191367098.PubMedPubMed CentralView Article
- Sorlie T, Tibshirani R, Parker J, Hastie T, Marron JS, Nobel A, Deng S, Johnsen H, Pesich R, Geisler S: Repeated observation of breast tumor subtypes in independent gene expression data sets. Proc Natl Acad Sci U S A. 2003, 100 (14): 8418-8423. 10.1073/pnas.0932692100.PubMedPubMed CentralView Article
- 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 (3): R3-10.1186/bcr433.PubMedPubMed CentralView Article
- 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: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415 (6871): 530-536. 10.1038/415530a.PubMedView Article
- Ma XJ, Wang Z, Ryan PD, Isakoff SJ, Barmettler A, Fuller A, Muir B, Mohapatra G, Salunga R, Tuggle JT: A two-gene expression ratio predicts clinical outcome in breast cancer patients treated with tamoxifen. Cancer Cell. 2004, 5 (6): 607-616. 10.1016/j.ccr.2004.05.015.PubMedView Article
- Huang E, Cheng SH, Dressman H, Pittman J, Tsou MH, Horng CF, Bild A, Iversen ES, Liao M, Chen CM: Gene expression predictors of breast cancer outcomes. Lancet. 2003, 361 (9369): 1590-1596. 10.1016/S0140-6736(03)13308-9.PubMedView Article
- Zhao H, Langerod A, Ji Y, Nowels KW, Nesland JM, Tibshirani R, Bukholm IK, Karesen R, Botstein D, Borresen-Dale AL: Different gene expression patterns in invasive lobular and ductal carcinomas of the breast. Mol Biol Cell. 2004, 15 (6): 2523-2536. 10.1091/mbc.E03-11-0786.PubMedPubMed CentralView Article
- Bertucci F, Finetti P, Rougemont J, Charafe-Jauffret E, Cervera N, Tarpin C, Nguyen C, Xerri L, Houlgatte R, Jacquemier J: Gene expression profiling identifies molecular subtypes of inflammatory breast cancer. Cancer Res. 2005, 65 (6): 2170-2178. 10.1158/0008-5472.CAN-04-4115.PubMedView Article
- Paik S, Shak S, Tang G, Kim C, Baker J, Cronin M, Baehner FL, Walker MG, Watson D, Park T: A multigene assay to predict recurrence of tamoxifen-treated, node-negative breast cancer. N Engl J Med. 2004, 351 (27): 2817-2826. 10.1056/NEJMoa041588.PubMedView Article
- van de Vijver MJ, He YD, van't Veer LJ, Dai H, Hart AA, Voskuil DW, Schreiber GJ, Peterse JL, Roberts C, Marton MJ: A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med. 2002, 347 (25): 1999-2009. 10.1056/NEJMoa021967.PubMedView Article
- Bhattacharjee A, Richards WG, Staunton J, Li C, Monti S, Vasa P, Ladd C, Beheshti J, Bueno R, Gillette M: Classification of human lung carcinomas by mRNA expression profiling reveals distinct adenocarcinoma subclasses. Proc Natl Acad Sci U S A. 2001, 98 (24): 13790-13795. 10.1073/pnas.191502998.PubMedPubMed CentralView Article
- Chung CH, Parker JS, Karaca G, Wu J, Funkhouser WK, Moore D, Butterfoss D, Xiang D, Zanation A, Yin X: Molecular classification of head and neck squamous cell carcinomas using patterns of gene expression. Cancer Cell. 2004, 5 (5): 489-500. 10.1016/S1535-6108(04)00112-6.PubMedView Article
- Garber ME, Troyanskaya OG, Schluens K, Petersen S, Thaesler Z, Pacyna-Gengelbach M, van de Rijn M, Rosen GD, Perou CM, Whyte RI: Diversity of gene expression in adenocarcinoma of the lung. Proc Natl Acad Sci U S A. 2001, 98 (24): 13784-13789. 10.1073/pnas.241500798.PubMedPubMed CentralView Article
- Michiels S, Koscielny S, Hill C: Prediction of cancer outcome with microarrays: a multiple random validation strategy. Lancet. 2005, 365 (9458): 488-492. 10.1016/S0140-6736(05)17866-0.PubMedView Article
- Jenssen TK, Hovig E: Gene-expression profiling in breast cancer. Lancet. 2005, 365 (9460): 634-635.PubMedView Article
- 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 (1): 14-18.PubMedView Article
- Ioannidis JP: Microarrays and molecular research: noise discovery?. Lancet. 2005, 365 (9458): 454-455.PubMedView Article
- 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 (18): 10393-10398. 10.1073/pnas.1732912100.PubMedPubMed CentralView Article
- Benito M, Parker J, Du Q, Wu J, Xiang D, Perou CM, Marron JS: Adjustment of systematic microarray data biases. Bioinformatics. 2004, 20 (1): 105-114. 10.1093/bioinformatics/btg385.PubMedView Article
- Hosack DA, Dennis G, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4 (10): R70-10.1186/gb-2003-4-10-r70.PubMedPubMed CentralView Article
- Usary J, Llaca V, Karaca G, Presswala S, Karaca M, He X, Langerod A, Karesen R, Oh DS, Dressler LG: Mutation of GATA3 in human breast tumors. Oncogene. 2004, 23 (46): 7669-7678. 10.1038/sj.onc.1207966.PubMedView Article
- Perou CM, Jeffrey SS, van de Rijn M, Rees CA, Eisen MB, Ross DT, Pergamenschikov A, Williams CF, Zhu SX, Lee JC: Distinctive gene expression patterns in human mammary epithelial cells and breast cancers. Proc Natl Acad Sci U S A. 1999, 96 (16): 9212-9217. 10.1073/pnas.96.16.9212.PubMedPubMed CentralView Article
- Chung CH, Bernard PS, Perou CM: Molecular portraits and the family tree of cancer. Nat Genet. 2002, 32 (Suppl): 533-540. 10.1038/ng1038.PubMedView Article
- Whitfield ML, Sherlock G, Saldanha AJ, Murray JI, Ball CA, Alexander KE, Matese JC, Perou CM, Hurt MM, Brown PO: Identification of genes periodically expressed in the human cell cycle and their expression in tumors. Mol Biol Cell. 2002, 13 (6): 1977-2000. 10.1091/mbc.02-02-0030..PubMedPubMed CentralView Article
- Bromberg JF, Horvath CM, Wen Z, Schreiber RD, Darnell JE: 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 (15): 7673-7678. 10.1073/pnas.93.15.7673.PubMedPubMed CentralView Article
- 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 (6): 1980-1991.PubMed
- Van Belle G, Fisher L: Biostatistics: a methodology for the health sciences. 2004, Hoboken, NJ: Wiley-Interscience John Wiley & Sons, 2View Article
- 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 (16): 1605-1616. 10.1056/NEJMoa031046.PubMedView Article
- Bair E, Tibshirani R: Semi-supervised methods to predict patient survival from gene expression data. PLoS Biol. 2004, 2 (4): E108-10.1371/journal.pbio.0020108.PubMedPubMed CentralView Article
- Chang HY, Nuyten DS, Sneddon JB, Hastie T, Tibshirani R, Sorlie T, Dai H, He YD, van't Veer LJ, Bartelink H: 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 (10): 3738-3743. 10.1073/pnas.0409462102.PubMedPubMed CentralView Article
- 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
- 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 (2): E7-10.1371/journal.pbio.0020007.PubMedPubMed CentralView Article
- Hu Z, Troester M, Perou CM: High reproducibility using sodium hydroxide-stripped long oligonucleotide DNA microarrays. Biotechniques. 2005, 38 (1): 121-124.PubMedView Article
- Novoradovskaya N, Whitfield ML, Basehore LS, Novoradovsky A, Pesich R, Usary J, Karaca M, Wong WK, Aprelikova O, Fero M: Universal Reference RNA as a standard for microarray experiments. BMC Genomics. 2004, 5 (1): 20-10.1186/1471-2164-5-20.PubMedPubMed CentralView Article
- UNC Microarray Database. [https://genome.unc.edu/]
- 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 (4): e15-10.1093/nar/30.4.e15.PubMedPubMed CentralView Article
- UNC Breast Tumor Data. [https://genome.unc.edu/pubsup/breastTumor/]
- Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/geo/]
- 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 (9): 5116-5121. 10.1073/pnas.091062498.PubMedPubMed CentralView Article
- 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 (6): 520-525. 10.1093/bioinformatics/17.6.520.PubMedView Article
- Diehn M, Sherlock G, Binkley G, Jin H, Matese JC, Hernandez-Boussard T, Rees CA, Cherry JM, Botstein D, Brown PO: SOURCE: a unified genomic resource of functional annotations, ontologies, and gene expression data. Nucleic Acids Res. 2003, 31 (1): 219-223. 10.1093/nar/gkg014.PubMedPubMed CentralView Article
- 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 (25): 14863-14868. 10.1073/pnas.95.25.14863.PubMedPubMed CentralView Article
- Hastie T, Tibshirani R, Friedman JH: The elements of statistical learning: data mining, inference, and prediction. 2001, New York: SpringerView Article
- 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 (10): 6567-6572. 10.1073/pnas.082099299.PubMedPubMed CentralView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.