The effect of the stromal component of breast tumours on prediction of clinical outcome using gene expression microarray analysis.

The aim of this study was to examine the effect of the cellular composition of biopsies on the error rates of multigene predictors of response of breast tumours to neoadjuvant adriamycin and cyclophosphamide (AC) chemotherapy. Core biopsies were taken from primary breast tumours of 43 patients prior to AC, and subsequent clinical response was recorded. Post-chemotherapy (day 21) samples were available for 16 of these samples. Frozen sections of each core were used to estimate the proportion of invasive cancer and other tissue components at three levels. Transcriptional profiling was performed using a cDNA array containing 4,600 elements. Twenty-three (53%) patients demonstrated a 'good' and 20 (47%) a 'poor' clinical response. The percentage invasive tumour in core biopsies collected from these patients varied markedly. Despite this, agglomerative clustering of sample expression profiles showed that almost all biopsies from the same tumour aggregated as nearest neighbours. SAM (significance analysis of microarrays) regression analysis identified 144 genes which distinguished high- and low-percentage invasive tumour biopsies at a false discovery rate of not more than 5%. The misclassification error of prediction of clinical response using microarray data from pre-treatment biopsies (on leave-one-out cross-validation) was 28%. When prediction was performed on subsets of samples which were more homogeneous in their proportions of malignant and stromal cells, the misclassification error was considerably lower (8%-13%, p < 0.05 on permutation). The non-tumour content of breast cancer samples has a significant effect on gene expression profiles. Consideration of this factor improves accuracy of response prediction by expression array profiling. Future gene expression array prediction studies should be planned taking this into account.