Comparison of gene expression profiles in core biopsies and corresponding surgical breast cancer samples
Breast Cancer Research volume 8, Article number: R51 (2006)
Gene expression profiling has been successfully used to classify breast cancer into clinically distinct subtypes, and to predict the risk of recurrence and treatment response. The aim of this study was to investigate whether the gene expression profile (GEP) detected in a core biopsy (CB) is representative for the entire tumor, since CB is an important tool in breast cancer diagnosis. Moreover, we investigated whether performing CBs prior to the surgical excision could influence the GEP of the respective tumor.
We quantified the RNA expression of 60 relevant genes by quantitative real-time PCR in paired CBs and surgical specimens from 22 untreated primary breast cancer patients. Subsequently, expression data were compared with independent GEPs obtained from tumors of 317 patients without preceding CB.
In 82% of the cases the GEP detected in the CB correlated very well with the corresponding profile in the surgical sample (rs ≥ 0.95, p < 0.001). Gene-by-gene analysis revealed four genes significantly elevated in the surgical sample compared to the CB; these comprised genes mainly involved in inflammation and the wound repair process as well as in tumor invasion and metastasis.
A GEP detected in a CB are representative for the entire tumor and is, therefore, of clinical relevance. The observed alterations of individual genes after performance of CB deserve attention since they might impact the clinical interpretation with respect to prognosis and therapy prediction of the GEP as detected in the surgical specimen following CB performance.
Gene expression profiling by parallel detection of thousands of genes permits the molecular signature (phenotype) of a tissue sample to be read and can, therefore, individually characterize a patient's tumor at the molecular level. Based on the gene expression profile (GEP) of a tumor, a molecular classification for breast cancer was proposed  and several molecular signatures were reported to predict the risk of recurrence and treatment response [2–5]. Such molecular analyses require only small amounts of material, such as tissue samples obtained by minimal invasive methods, for example, core biopsy (CB), which are used to assess the nature of palpable and non-palpable breast lesions to confirm or exclude the diagnosis of breast cancer [6–9].
Although CB investigations have become more and more important in the early workup of breast lesions, there are only a few investigations regarding the reliability of GEPs as detected in CBs . However, wound healing subsequent to CB can potentially induce gene expression alterations in the injured tissue. Many of these normally occurring reparative processes share molecular characteristics with an aggressive cancer phenotype, such as cell proliferation and survival, angiogenesis, and extracellular matrix remodeling; these biological hallmarks were shown to predict the clinical course in cancer [3, 11–13]. Therefore, potential biological changes induced by CB require further study since they may have important consequences for clinical outcome prediction and treatment decisions as well as the interpretation of GEP changes investigated in neo-adjuvant studies.
We explored the molecular expression levels of 60 genes using quantitative real-time PCR (qrt-PCR), a highly sensitive and reproducible method, in paired CB and surgical samples. These genes were selected according to their known links to malignant cell behavior in breast cancer and their importance in major cancer hallmarks, such as proliferation, survival, invasiveness and angiogenic potential, and in the wound healing process . Our first objective was to investigate whether the molecular profile of a CB is representative for the whole tumor. The second objective was to study if the repair process following CB alters the GEP and if this is influenced by the timeframe between the CB and the surgical excision.
Materials and methods
Tumor specimen acquisition
Between June 2004 and June 2005, 22 consecutive breast cancer patients entered this study, for which both CBs and surgical tumor tissue specimens (STs) were available. This study was undertaken at the Women's University Hospital Basel, Switzerland, and approved by the local institutional review board (EKBB permission Nr. 81/04). Written informed consent was obtained from all patients.
All CBs (14-gauge needle, Magnum® Core high speed, Bard Medica, Karlsruhe, Germany) were obtained under sonographic guidance (HDI 5'000 Sono CT®, Philips, Zurich, Switzerland) under local anesthesia through a skin incision in a sterile field. Five biopsy specimens were taken routinely for each patient. Two biopsy specimens were divided longitudinally; two halves to be used for molecular examination were, within one minute, stored in RNAlater®-solution (Qiagen, Basel, Switzerland), while the other two halves as well as the other three biopsies were immediately put in formaldehyde and sent for histological examination to the Institute of Pathology, University Hospital Basel.
After establishing the diagnosis of breast cancer, all patients underwent breast surgery with sentinel node lymphonodectomy or axillary lymphonodectomy if indicated. All surgical samples were examined by the same pathologist. If the tumor tissue was larger than 0.5 cm in diameter at the intra-operative frozen section, a representative piece containing more than 60% tumor cells was cryopreserved within five minutes and made available to the Stiftung Tumorbank Basel for molecular examination. The rest of the tumor tissue was embedded in paraffin for routine histological examination.
Reference study population
RNA expression levels of all 60 genes were detected using the same qrt-PCR method in 317 surgically excised breast cancer specimens  from patients undergoing primary surgery in 1992 to 1996 without previous examination by CB. All tissue samples were prepared by the pathologists as described above, all samples contained more than 60% tumor cells and they were cryoperserved within five minutes. The Stiftung Tumorbank Basel was subsequently responsible for collection, storage at -70°C and analysis.
RNA extraction and qrt-PCR
Detailed procedures have been published elsewhere . In brief, RNA was extracted using the RNeasy Mini Kit (Qiagen, Inc., Valencia, CA, USA), quantified and quality-checked on a Bioanalyzer 2100 (RNA 6000 Nano LabChip-Kit, Agilent Technologies, New Castle, DE, USA). High quality RNA samples were reverse-transcribed (10 mM DDT, 1 μg of hexamer primers, 2 U of MMLV Reverse Transcriptase (Invitrogen, Basel, Switzerland), 40 U of RNasin (Promega, Wallisellen, Switzerland), 0.5 mM of each dNTP (Promega), 1× reaction buffer). PCR primers were designed to be cDNA specific and ordered at GeneScan Europe (Freiburg, Germany). PCR was performed in 40 cycles on an ABI Prism 7000 using 2× SYBR Green I Master Mix (Applied Biosystems, Forster City, CA, USA) in a final volume of 25 μl. Relative quantities (ΔCt) were obtained by normalization against ribosomal 18S RNA, and standardization was achieved with Human Universal Standard RNA (Stratagene Europe, Amsterdam, The Netherlands). The 60 genes quantitatively assessed are listed in Table 1.
The same amount of RNA was used for the GEP analysis of each sample. For statistical analysis, ΔCt expression values of each gene were obtained by normalizing the raw gene values to 18S rRNA as a reference gene.
Cluster and TreeView programs were used to perform unsupervised hierarchical clustering of samples and genes (Spearman correlation, average linkage) . Spearman correlation coefficients were calculated to compare the GEPs of all paired samples. Differentially expressed genes were identified with the paired two-sample t-test. The Mann-Whitney U-test was used to compare median expression values of genes among different subgroups of patients. All statistical analyses were carried out at 5% level of significance and performed with S-Plus software (Version 6.1, Insightful Corporation, Seattle, WA, USA).
The mean age of the 22 patients was 63 years (range: 34 to 81 years). The period between CB and final surgery ranged from 1 to 23 days. The mean tumor diameter was 2.9 cm (range: 0.9 to 9 cm). Clinicopathological characteristics of each patient are listed in Table 1. Of note, one paired tissue sample was not further evaluated due to poor RNA quality, and RNA extracted from the two halves of the CB was not pooled in three cases but analyzed separately.
Comparison of the gene expression profile in paired samples
As shown in Figure 1, unsupervised hierarchical clustering revealed that paired CB and ST generally clustered together; in only four cases (patients 4, 10, 11, and 12) did the GEP of the CB not agglomerate with the profile of the respective ST. Interestingly, the two separate CBs taken from patient 11 were very similar to each other, although they differed from their ST GEP. The gene dendrogram of the cluster analysis also revealed that samples agglomerated in two main groups according to their respective estrogen receptor status, reconfirming the representative value of this study population.
Subsequent analysis of paired CBs and STs confirmed the high correlation between all samples (rs from 0.86 to 0.98, all p < 0.001) for each patient. A scatter plot of two representative examples of a paired CB/ST is displayed in Figure 2. The differences in paired GEPs does not seem to be related to the timeframe between CB and surgery or to any other clinicopathological parameters (Table 1).
Gene-by-gene analysis in core biopsies and paired surgical specimens
The comparison of the expression levels of individual genes by means of paired t-test showed no significant difference between CB and ST with the exception of four genes. Plasminogen activator inhibitor 1 (PAI-1; also known as SERPINE1) was significantly higher expressed (lower ΔCt values) in STs compared to CBs (p < 0.001, Table 1). Similar differences, although less pronounced, were observed for cyclooxygenase 2 (COX-2; also known as PTGS2; p < 0.001), urokinase plasminogen activator receptor (uPAR; also known as PLAUR; p = 0.003) and matrix metalloproteinase 1 (MMP1; p = 0.03). The increase in the expression of these genes was not related to the timeframe between CB and surgery. All other genes were very similarly expressed in paired CB/ST as shown in Figure 3. Table 1 lists differences in RNA expression values of PAI-1, COX-2 and ERBB2 for paired samples of each patient.
Histological re-examination of cryocuts of the surgical specimens revealed a certain amount of inflammation and fibrolysis. Whether these observations are due to cancerogenesis or to a de novo induced wound repair process can not be determined.
Comparison of the expression levels of selected genes with a reference study population of surgical tumor tissue specimens
To verify whether the higher expression levels of PAI-1 and COX-2 observed in STs could have been induced by the preceding CB procedure, we compared the expression levels of the same genes in an independent population of 317 primary breast cancer patients. These samples were investigated with the same qrt-PCR technique but were from patients from which no CB had been taken prior to surgery. As shown in Figure 4, the expression levels of PAI-1 and COX-2 measured in the independent STs without CB were found to be very similar to the levels detected in CBs and significantly different from those detected in STs after CB. Moreover, no variation at all was observed for the remaining genes as illustrated by the expression levels of ERBB2 as an example.
Ultrasound-guided CB is a well established method to diagnose breast cancer in women, since it is a reliable, and time- and cost-saving method. The clinical utility of the information gained by CB depends on whether the CB is representative for the whole tumor. Our data demonstrate that the quantitative expression levels of 60 genes detected in CBs were highly comparable to their paired STs in 17 out of the 21 cases investigated. Even in the cases where GEPs of a CB and ST did not agglomerate, the expression levels of the ER and progesterone receptor as well as ERBB2 measured in the CB were also representative for the whole tumor. This is important since today's therapy decisions are based on these markers obtained either by CB or ST (St Gallen consensus recommendations ). In addition, our results reconfirm previous observations reported by immunohistochemistry [19, 20] or semi-qrt-PCR .
Tissue sampling by CB causes a local injury, inducing wound healing that is characterized by recruitment of inflammatory cells, stimulation of stromal and epithelial cell proliferation, cell migration and increased angiogenesis. Analysis on a gene-by-gene basis demonstrated higher expression levels of PAI-1, COX-2, uPAR and MMP1 in STs compared to their paired CBs, whereas no changes were observed for all other genes. These results are not surprising since proteinases (such as PAI-1 and uPAR) are known to be essentially involved in the wound healing process [22, 23] and COX-2 plays roles in inflammation and angiogenesis [24–27].
However, many of these reparation processes show parallels with cancerogenesis [28–30]; while proteinases, their inhibitors, cyto-/chemokines and growth factors are essential for wound healing and tissue repair, they also play central roles in cancer progression. For example, uPA, uPAR and its inhibitor PAI-1 are responsible for the degradation and remodeling of the extracellular matrix, and are further involved in angiogenesis, cell adhesion and migration necessary for tumor cell invasion and metastasis [31, 32]. COX-2 can be induced by cytokines and growth factors during the inflammatory repair process as well as in cancer [24–26, 30] resulting in COX-2 overexpression observed in human malignancies [25–27].
Therefore, increased levels of these markers in the tumor specimen could suggest a more aggressive cancer phenotype. Indeed, elevated levels of uPA and PAI-1 are associated with poor clinical outcome in breast cancer and also have predictive value [33–35]. Moreover, a previously identified 'wound-response signature' turned out to be prognostic in several carcinomas, including breast cancer [12, 13]. Although COX-2 has been associated with increased Vascular Endothelial Growth Factor A (VEGF), estrogen synthesis, proliferation, apoptosis and invasion [25, 27, 36], in our study, higher levels of COX-2 were not accompanied by changes in the expression of genes involved in these processes, indicating that the observed molecular alterations influence data interpretation but not tumor aggressiveness.
Our study demonstrates that expression levels of ER, progesterone receptor, ERBB2 and other genes relevant for the management of breast cancer as detected in CBs are representative for the whole tumor. However, increased expression levels of proteinases (e.g. PAI-1, uPAR, MMP1) and COX-2 in STs compared to their paired CBs suggest induction of theses genes during the repair process following tissue injury caused by CBs. This observation is important since such molecular alterations may have an impact on the clinical interpretation of GEPs detected in STs with respect to the prediction of risk assessment and treatment response.
gene expression profile
quantitative real-time PCR
plasminogen activator inhibitor
Plasminogen Activator Inhibitor-1
surgical tumor tissue specimen
urokinase plasminogen activator receptor
Vascular Endothelial Growth Factor A.
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 USA. 2001, 98: 10869-10874. 10.1073/pnas.191367098.
Chang JC, Wooten EC, Tsimelzon A, Hilsenbeck SG, Gutierrez MC, Elledge R, Mohsin S, Osborne CK, Chamness GC, Allred DC, et al: Gene expression profiling for the prediction of therapeutic response to docetaxel in patients with breast cancer. Lancet. 2003, 362: 362-369. 10.1016/S0140-6736(03)14023-8.
Chang JC, Weiss H, Mohsin S: Apoptosis induction by the core biopsy procedure. J Clin Oncol. 2005, 23: 7240-10.1200/JCO.2005.02.5197.
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. 10.1056/NEJMoa021967.
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. 10.1038/415530a.
Fornage BD, Sneige N, Edeiken BS: Interventional breast sonography. Eur J Radiol. 2002, 42: 17-31. 10.1016/S0720-048X(01)00481-8.
Liberman L: Percutaneous image-guided core breast biopsy. Radiol Clin North Am. 2002, 40: 483-500. 10.1016/S0033-8389(01)00011-2. vi
Parker SH, Burbank F, Jackman RJ, Aucreman CJ, Cardenosa G, Cink TM, Coscia JL, Eklund GW, Evans WP, Garver PR, et al: Percutaneous large-core breast biopsy: a multi-institutional study. Radiology. 1994, 193: 359-364.
Philpotts LE: Controversies in core-needle breast biopsy. Semin Roentgenol. 2001, 36: 270-283. 10.1053/sroe.2001.25121.
Symmans WF, Ayers M, Clark EA, Stec J, Hess KR, Sneige N, Buchholz TA, Krishnamurthy S, Ibrahim NK, Buzdar AU, et al: Total RNA yield and microarray gene expression profiles from fine-needle aspiration biopsy and core-needle biopsy samples of breast carcinoma. Cancer. 2003, 97: 2960-2971. 10.1002/cncr.11435.
Menard S, Pupa SM, Campiglio M, Taglibue E, Balsari A, Fagnoni F, Costa A: Apoptosis induction by trastuzumab: possible role of the core biopsy intervention. J Clin Oncol. 2005, 23: 7238-7240. 10.1200/JCO.2005.02.4679.
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-10.1371/journal.pbio.0020007.
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 USA. 2005, 102: 3738-3743. 10.1073/pnas.0409462102.
Hanahan D, Weinberg RA: The hallmarks of cancer. Cell. 2000, 100: 57-70. 10.1016/S0092-8674(00)81683-9.
Urban P, Vuaroqueaux V, Labuhn M, Delorenzi M, Wirapati P, Wight E, Senn H-J, Benz C, Eppenberger U, Eppenberger-Castori S: Increased expression of urokinase-type plasminogen activator (uPA) mRNA determines adverse prognosis in ErbB2-positive primary breast cancer. J Clin Oncol.
Labuhn M, Vuaroqueaux V, Fina F, Schaller A, Nanni-Metellus I, Küng W, Eppenberger-Castori S, Martin PM, Eppenberger U: Simultaneous quantitative detection of relevant biomarkers in breast cancer by quantitative real-time PCR. Int J Biol Markers. 2006, 21: 30-39.
Eisen M, Spellman P, Brown P, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.
Goldhirsch A, Glick JH, Gelber RD, Coates AS, Thurlimann B, Senn HJ: Meeting highlights: international expert consensus on the primary therapy of early breast cancer 2005. Ann Oncol. 2005, 16: 1569-1583. 10.1093/annonc/mdi326.
Cavaliere A, Sidoni A, Scheibel M, Bellezza G, Brachelente G, Vitali R, Bucciarelli E: Biopathologic profile of breast cancer core biopsy: is it always a valid method?. Cancer Lett. 2005, 218: 117-121. 10.1016/j.canlet.2004.07.041.
Badoual C, Maruani A, Ghorra C, Lebas P, Avigdor S, Michenet P: Pathological prognostic factors of invasive breast carcinoma in ultrasound-guided large core biopsies-correlation with subsequent surgical excisions. Breast. 2005, 14: 22-27. 10.1016/j.breast.2004.07.005.
Pohlodek K, Galbavy S, Bartosova M, Mucha V, Holoman K, Pastorek J, Pastorekova S: Semi-quantitative RT-PCR assessment of molecular markers in breast large-core needle biopsies. Neoplasma. 2004, 51: 415-421.
Baker EA, Leaper DJ: Proteinases, their inhibitors, and cytokine profiles in acute wound fluid. Wound Repair Regen. 2000, 8: 392-398.
Stadelmann WK, Digenis AG, Tobin GR: Physiology and healing dynamics of chronic cutaneous wounds. Am J Surg. 1998, 176: 26S-38S. 10.1016/S0002-9610(98)00183-4.
Williams CS, Mann M, DuBois RN: The role of cyclooxygenases in inflammation, cancer, and development. Oncogene. 1999, 18: 7908-7916. 10.1038/sj.onc.1203286.
Howe LR, Dannenberg AJ: COX-2 inhibitors for the prevention of breast cancer. J Mammary Gland Biol Neoplasia. 2003, 8: 31-43. 10.1023/A:1025731204719.
Corcoran CA, He Q, Huang Y, Sheikh MS: Cyclooxygenase-2 interacts with p53 and interferes with p53-dependent transcription and apoptosis. Oncogene. 2005, 24: 1634-1640. 10.1038/sj.onc.1208353.
Ma X, Yang Q, Wilson KT, Kundu N, Meltzer SJ, Fulton AM: Promoter methylation regulates cyclooxygenase expression in breast cancer. Breast Cancer Res. 2004, 6: R316-321. 10.1186/bcr793.
Coussens LM, Werb Z: Inflammation and cancer. Nature. 2002, 420: 860-867. 10.1038/nature01322.
Dvorak HF: Tumors: wounds that do not heal. Similarities between tumor stroma generation and wound healing. N Engl J Med. 1986, 315: 1650-1659.
Balkwill F, Mantovani A: Inflammation and cancer: back to Virchow?. Lancet. 2001, 357: 539-545. 10.1016/S0140-6736(00)04046-0.
Harbeck N, Kates RE, Gauger K, Willems A, Kiechle M, Magdolen V, Schmitt M: Urokinase-type plasminogen activator (uPA) and its inhibitor PAI-I: novel tumor-derived factors with a high prognostic and predictive impact in breast cancer. Thromb Haemost. 2004, 91: 450-456.
Duffy MJ, Duggan C: The urokinase plasminogen activator system: a rich source of tumour markers for the individualised management of patients with cancer. Clin Biochem. 2004, 37: 541-548. 10.1016/j.clinbiochem.2004.05.013.
Decock J, Paridaens R, Cufer T: Proteases and metastasis: clinical relevance nowadays?. Curr Opin Oncol. 2005, 17: 545-550. 10.1097/01.cco.0000180435.39614.63.
Cianfrocca M, Goldstein LJ: Prognostic and predictive factors in early-stage breast cancer. Oncologist. 2004, 9: 606-616. 10.1634/theoncologist.9-6-606.
Look M, van Putten W, Duffy M, Harbeck N, Christensen IJ, Thomssen C, Kates R, Spyratos F, Ferno M, Eppenberger-Castori S, et al: Pooled analysis of prognostic impact of uPA and PAI-1 in breast cancer patients. Thromb Haemost. 2003, 90: 538-548.
Singh B, Berry JA, Shoher A, Ramakrishnan V, Lucci A: COX-2 overexpression increases motility and invasion of breast cancer cells. Int J Oncol. 2005, 26: 1393-1399.
This work was in part supported by the generous donation of Mrs Bärbel Girgner. We are grateful to Monika Heinemann, Sabine Ehret and Ulrike Weissenstein for technical assistance and data management and to Linda Herberich for proof reading.
The authors declare that they have no competing interests.
RS and SE designed the study, presented it to the ethical committee and drafted the manuscript. RS performed the sonographic guided core biopsies and coordinated the clinical part of the study. VV and ML contributed to the selection of the genes, selected primers and supervised RNA extraction and qrt-PCR. VV and PU contributed to statistical analysis, data interpretation and manuscript revision. GS performed the pathological analysis of all samples, prepared the surgical samples for molecular analysis and participated in data interpretation. EW, UE and WH participated in designing the study and writing the manuscript. SE coordinated the molecular investigation and performed statistical analysis. All authors read and approved the final manuscript.
About this article
Cite this article
Zanetti-Dällenbach, R., Vuaroqueaux, V., Wight, E. et al. Comparison of gene expression profiles in core biopsies and corresponding surgical breast cancer samples. Breast Cancer Res 8, R51 (2006). https://doi.org/10.1186/bcr1542