A gene expression profile for detection of sufficient tumour cells in breast tumour tissue: microarray diagnosis eligibility
© Roepman et al; licensee BioMed Central Ltd. 2009
Received: 23 January 2009
Accepted: 12 August 2009
Published: 12 August 2009
Microarray diagnostics of tumour samples is based on measurement of prognostic and/or predictive gene expression profiles. Typically, diagnostic profiles have been developed using bulk tumour samples with a sufficient amount of tumour cells (usually >50%). Consequentially, a diagnostic results depends on the minimal percentage of tumour cells within a sample. Currently, tumour cell percentage is assessed by conventional histopathological review. However, even for experienced pathologists, such scoring remains subjective and time consuming and can lead to ambiguous results.
In this study we investigated whether we could use transcriptional activity of a specific set of genes instead of histopathological review to identify samples with sufficient tumour cell content. Genome-wide gene expression measurements were used to develop a transcriptional gene profile that could accurately assess a sample's tumour cell percentage.
Supervised analysis across 165 breast tumour samples resulted in the identification of a set of 13 genes which expression correlated with presence of tumour cells. The developed gene profile showed a high performance (AUC 0.92) for identification of samples that are suitable for microarray diagnostics. Validation on 238 additional breast tumour samples indicated a robust performance for correct classification with an overall accuracy of 91 percent and a kappa score of 0.63 (95%CI 0.47–0.73).
The developed 13-gene profile provides an objective tool for assessment whether a breast cancer sample contains sufficient tumour cells for microarray diagnostics. It will improve the efficiency and throughput for diagnostic gene expression profiling as it no longer requires histopathological analysis for initial tumour percentage scoring. Such profile will also be very use useful for assessment of tumour cell percentage in biopsies where conventional histopathology is difficult, such as fine needle aspirates.
Microarray diagnostics of tumour specimens is based on gene expression measurement of a specific set of predictive or prognostic genes. Bulk tumour samples that consist of tumour cells admixed with surrounding stromal tissue are commonly used for microarray analysis. Although tumour stroma likely plays an important role in tumour development and metastasis [1–3], gene expression profiles are typically generated using tissue that contains sufficient amount of tumour cells, not stroma. Most prognostic gene profiles were originally identified using samples containing at least 50% tumour cells and are, therefore, likely based on gene expression of the neoplastic tissue in question. However, some profiles, e.g. MammaPrint [4, 5] have been shown to be accurate with 30% tumour in the diagnostic specimen (AM Glas, unpublished data).
Currently, tumour cell percentage assessment is assessed using heamatoxilin-eosine (H/E) stained specimen for histopathological review. However, even for experienced pathologists, histopathological tumour scoring remains subjective and time consuming and can lead to inconclusive results [6–8] and variable tumour cell percentage scoring. Moreover, tumour cell scoring can be more difficult for core biopsies and especially fine-needle aspirates that are too small or unsuitable for H/E analysis. A tumour percentage scoring method based on tumour cell mRNA transcription levels would provide an additional method to more reliably determine tumour cell content in a reduced time frame.
Herein we report the development of a molecular profile that can accurately identify breast cancer samples with sufficient tumour cell content for diagnostic purpose. This assessment is based on transcription levels and is able to reduce the time that is needed for a microarray experiment as no initial pathological tumour cell percentage (TCP) scoring will be required before sample processing. Moreover, the identified tumour percentage profile would facilitate microarray diagnostics of small specimens for which tumour sections can not be generated for pathological scoring.
Four hundred and three frozen tumour samples or tumour samples preserved in RNALater from breast cancer patients were used. At the time of initial diagnosis, all patients had provided consent in the sense that their tumour samples could be used for investigational purposes. The study was carried out in accordance with the ethical standards of the Helsinki Declaration and was approved by the Medical ethical Board of the participating medical centers and hospitals. All samples were de-identified, analyzed anonymously and were part of research implementation studies for MammaPrint.
Genes used for building a tumour cell percentage associated profile.
Homo sapiens cDNA clone IMAGE:1473008 3'
Homo sapiens cDNA FLJ21777 fis, clone HEP00173
Homo sapiens cDNA FLJ37541 fis, clone BRCAN2026340
Homo sapiens anaphase promoting complex subunit 7
Homo sapiens cDNA clone IMAGE:4837645
Homo sapiens mRNA for hypothetical protein FLJ20694 variant, clone: COL06209
Homo sapiens HEAT repeat containing 3
Homo sapiens hypothetical protein LOC147804
Homo sapiens O-linked N-acetylglucosamine (GlcNAc) transferase
Homo sapiens proline rich 13 (PRR13)
Homo sapiens synaptosomal-associated protein, 29 kDa
Homo sapiens HSPC084 mRNA, partial cds
Homo sapiens tropomyosin 3 (TPM3)
Biological function analysis was performed using Ingenuity Pathway Analysis (IPA version 6.3, Ingenuity Systems Inc, Redwood City, CA). Full genome gene expression measurement of all samples analysed in this study is publicly available at the Gene Expression Omnibus (GEO) with accession number GSE16201.
Pathological tumour cell percentage related gene expression in breast cancer
To determine the variation in histopathological tumour cell percentage scoring, H/E stained tumour sections were repeatably analysed in a period of six months by five different pathologists (Figure 1). Although the variation between pathologists (10–15%) was larger than the variation observed in repeated scoring over time by the same pathologist (5–10%), the average variation in tumour cell percentage scoring was 10 percent (range 0 to 60%). This observed pathological variation (10%) is taken into account for classifier development as described in detail in the materials and methods.
Using 95 breast tumour samples that ranged in tumour cell content from zero to 85 percent, based on analysis of HE slides, a tumour percentage gene expression profile was identified. A 3-fold cross validation (CV) method was used to identify a gene expression profile (nearest-mean classifier) that showed a strong association with pathological TCP. This procedure was repeated five-hundred times (multiple sampling approach) and genes were ranked according to their CV performance. The multiple sampling CV performance indicated a strong association of the selected nearest-mean classifier model with pathological TCP (R2 = 0.42, P < 0.001, Wilcoxon rank-sum test). To identify the optimal/minimal number of genes to be used in the classifier model, the top-ranked gene set was gradually expanded and the classifier performance (area under ROC curve, AUC) was determined for each gene set size.
Identification of a 13-gene profile for assessment of TCP
Although optimal performance was reached with the complete set of 35 genes, a substantial number of genes from the 35-gene set showed a significant single-gene performance (AUC >0.80, 12 genes; AUC >0.80, 15 genes). A minimal combination of just two random genes from within the 35-gene set resulted in an AUC of greater than 0.80. This performance increased to >0.85 for random set of least 5 genes (data not shown).
Tumour cell percentage profile (TCP) performance.
Median TCP profile index
Training cohorts (n = 165)
low pathological TCP (<30%)
medium/high pathological TCP (≥30%)
P < 0.0001
Kappa score: 0.68 (95% 0.50 – 0.72)
Validation cohort (n = 238)
low pathological TCP (<30%)
medium/high pathological TCP (≥30%)
P < 0.0001
Kappa score: 0.63 (95%CI 0.47–0.73)
Biological function analysis indicated that the 35 TCP associated gene set was enriched for genes associated with cancer related processes such as cellular morphology (MGAT3, AQP1, TPM3, MYLK, SFRP1, NTRK2: p = 1.7e-3), cell-to-cell signalling (JAM2, LOC338328, OGT, MYLK, NTRK2, SNAP29: P = 1.7e-3), cellular movement and invasion (MGAT3, TPM3, MYLK, SFRP1, NTRK2, CCL15: P = 1.7e-3) and cellular survival and apoptosis (TPM3, PRR13, OGT, MYLK, SFRP1: P = 3.6e-5).
Validation of the 13-gene profile on a cohort of 238 breast tumour samples
Analysis of continuous TCP profile index and continuous pathological tumour cell percentage scoring resulted in a significant association between both methods (R2 = 0.48, Wilcoxon P < 0.001, AUC 0.90) and indicated that the profile might also be useful for indication of accurate tumour cell percentage instead of a low-, medium- or high-TCP classification. Statistical analysis of the profile indexes indicated a significant difference in outcome between samples with a low and medium TCP (P = 0.01, Student's T-test) and between samples with a low and high TCP (P < 0.0001) (Figure 3B). Although the difference in index between medium and high TCP samples was borderline significant (P = 0.02), the majority of medium TCP samples (78%) were classified as high TCP by the molecular profile. This result indicates that both medium and high TCP harbour a different tumour cell related gene profile compared to low TCP samples. The validation results confirm that the molecular profile can accurately determine samples with a sufficient tumour percentage for microarray diagnostics based on transcriptional analysis of 13 identified genes.
This study reports the development of a tumour cell percentage (TCP) assessment method for breast cancer samples based on transcriptional analysis of 13 genes. The 13-gene molecular profile has been validated on an independent cohort of 238 samples and can accurately identify breast tumour samples with a sufficient number of tumour cells for microarray diagnostics. Tumour percentage scoring based on the molecular profile is identical to the pathologist scoring for more than 90 percent of all analysed tumour samples. Although the variation in pathologist TCP scoring is in agreement with the number of discrepancies between pathological and molecular profile classification, this inconsistency is likely caused by the difference between the number of tumour cells present in a tumour sample and the actual tumour cell specific mRNA levels. In addition to the known variation and inconsistency in pathological scoring of tumour slides (e.g. Figure 1 and [6, 7]), Hsu et al. indicated that formalin fixation of slides can result in tumour cell shrinkage  and can lead to an underestimation of tumour cell content. Since the TCP profile is based on transcriptional levels of high TCP related gene expression within fresh or frozen tumour tissue, we believe that the developed gene profile likely gives a better indication whether a sample is suitable for microarray diagnostics compared to a pathological tumour cell percentage scoring on formalin-fixed H/E stained slides.
The utility of the gene profile lies in its capability to identify tumours with a high percentage of tumour cells compared to tumours with insufficient tumour content for subsequent microarray diagnostics. Conventional histopathological review results in tumour cell percentage scorings up to 10% increments but is laborious and requires an experienced (in-house) pathologist. The profile, on the other hand, is able to distinguish samples with low-, medium- or high-tumour cell content. Although the 13-gene profile provides a more qualitative measurement compared to the quantitative pathological assessment, it is a value tool for an objective TCP scoring based on transcriptional levels that can quickly identify samples suitable for diagnostics.
Since the gene profile was developed in such a way to mimic a pathological scoring which has been described as inconsistent and subjective, one might argue that the TCP profile also suffers from these factors. The main goal of this study, however, was not to develop a more accurate tool for TCP assessment but an objective method that can be performed independent of pathological expertise and which is based on transcriptional gene levels instead of the number of tumour cells. The use of a robust cross validation procedure that included a mimicked pathological variation was therefore included in the selection model that should, in principle, select a set of genes that are robust to this variation. Nevertheless, the 10 percent misclassification between the gene profile and pathological scoring might partly be attributed to this phenomenon.
Future application of the developed TCP profile could significantly improve the throughput of microarray diagnostics. Currently, after sample arrival, RNA processing and expression analysis cannot proceed until a histopathological TCP analysis confirms that the tumour sample contains sufficient tumour cells for analysis. As indicated above, pathological analysis will remain necessary for detection of ductal carcinoma in situ, necrosis and a detailed assessment about the percentage of tumour cells to define the suitability of the specimen. However, replacement of initial pathologist TCP scoring by a faster transcription based analysis that is able to identify whether a sample likely contains sufficient tumour tissue will shorten the processing time for microarray diagnostics. More importantly, the 13-gene molecular profile for tumour cell percentage (TCP) enables gene expression diagnostics for small samples or those on which no H/E evaluation is possible.
While gene expression profiles tend to be more robust with inclusion of a larger set of genes , we decided to limit the developed profile to a relatively small number of genes with optimal performance. The rationale behind this strategy was that transcriptional TCP assessment is preferably done before diagnostic microarray analysis as this assessment will indicate whether a sample is qualified for gene expression profiling. Future development of the 13-gene molecular profile into a RT-qPCR based assay, will allow a less subjective qualification of breast tumour samples as suitable for microarray diagnostics. This way, only samples with a sufficient TCP will be used for microarray diagnostics, saving time, money, and eliminating the need for a pathologist to score TCP on qualified specimens.
Whether one uses a qPCR assay for specimen selection or uses parallel microarray readout for final approval, the developed 13-gene TCP profile will provide an additional tool for objective assessment of sufficient tumour cell content in breast cancer tissue. It will improve the efficiency and throughput for diagnostic gene expression profiling as it eliminates the need for histopathological analysis in initial tumour percentage scoring. Moreover, it will also allows likely qualify small clinical samples, such as fine-needle aspirates, for transcriptional diagnostics.
The authors like to thank Diederik Wehkamp and Tako Bruinsma for bioinformatics support and Johan Westerga for pathological review of analyzed tumour samples. Also do we want to thank Femke de Snoo, Richard Bender, Frits Michiels, René Bernards and Laura van't Veer for fruitful discussions and critical review of this manuscript.
- Mueller MM, Fusenig NE: Friends or foes – bipolar effects of the tumour stroma in cancer. Nat Rev Cancer. 2004, 4: 839-849. 10.1038/nrc1477.View ArticlePubMedGoogle Scholar
- Christofori G: New signals from the invasive front. Nature. 2006, 441: 444-450. 10.1038/nature04872.View ArticlePubMedGoogle Scholar
- Kalluri R, Zeisberg M: Fibroblasts in cancer. Nat Rev Cancer. 2006, 6: 392-401. 10.1038/nrc1877.View ArticlePubMedGoogle Scholar
- 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, Schreiber GJ, Kerkhoven RM, Roberts C, Linsley PS, Bernards R, Friend SH: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415: 530-536. 10.1038/415530a.View ArticleGoogle Scholar
- Glas AM, Floore A, Delahaye LJ, Witteveen AT, Pover RC, Bakx N, Lahti-Domenici JS, Bruinsma TJ, Warmoes MO, Bernards R, Wessels LF, Van't Veer LJ: Converting a breast cancer microarray signature into a high-throughput diagnostic test. BMC Genomics. 2006, 7: 278-10.1186/1471-2164-7-278.View ArticlePubMedPubMed CentralGoogle Scholar
- Furness PN, Taub N, Assmann KJ, Banfi G, Cosyns JP, Dorman AM, Hill CM, Kapper SK, Waldherr R, Laurinavicius A, Marcussen N, Martins AP, Nogueira M, Regele H, Seron D, Carrera M, Sund S, Taskinen EI, Paavonen T, Tihomirova T, Rosenthal R: International variation in histologic grading is large, and persistent feedback does not improve reproducibility. Am J Surg Pathol. 2003, 27: 805-810. 10.1097/00000478-200306000-00012.View ArticlePubMedGoogle Scholar
- Ross JS, Symmans WF, Pusztai L, Hortobagyi GN: Standardizing slide-based assays in breast cancer: hormone receptors, HER2, and sentinel lymph nodes. Clin Cancer Res. 2007, 13: 2831-2835. 10.1158/1078-0432.CCR-06-2522.View ArticlePubMedGoogle Scholar
- Kirkegaard T, Edwards J, Tovey S, McGlynn LM, Krishna SN, Mukherjee R, Tam L, Munro AF, Dunne B, Bartlett JM: Observer variation in immunohistochemical analysis of protein expression, time for a change?. Histopathology. 2006, 48: 787-794. 10.1111/j.1365-2559.2006.02412.x.View ArticlePubMedGoogle Scholar
- Roepman P, Kemmeren P, Wessels LF, Slootweg PJ, Holstege FC: Multiple robust signatures for detecting lymph node metastasis in head and neck cancer. Cancer Res. 2006, 66: 2361-2366. 10.1158/0008-5472.CAN-05-3960.View ArticlePubMedGoogle Scholar
- Hsu PK, Huang HC, Hsieh CC, Hsu HS, Wu YC, Huang MH, Hsu WH: Effect of formalin fixation on tumor size determination in stage I non-small cell lung cancer. Ann Thorac Surg. 2007, 84: 1825-1829. 10.1016/j.athoracsur.2007.07.016.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/2/52/prepub