Impact of biospecimens handling on biomarker research in breast cancer
© De Cecco et al; licensee BioMed Central Ltd. 2009
Received: 6 July 2009
Accepted: 24 November 2009
Published: 24 November 2009
Gene expression profiling is moving from the research setting to the practical clinical use.
Gene signatures able to correctly identify high risk breast cancer patients as well as to predict response to treatment are currently under intense investigation. While technical issues dealing with RNA preparation, choice of array platforms, statistical analytical tools are taken into account, the tissue collection process is seldom considered.
The time elapsed between surgical tissue removal and freezing of samples for biological characterizations is rarely well defined and/or recorded even for recently stored samples, despite the publications of standard operating procedures for biological sample collection for tissue banks.
Breast cancer samples from 11 patients were collected immediately after surgical removal and subdivided into aliquots. One was immediately frozen and the others were maintained at room temperature for respectively 2, 6 and 24 hrs. RNA was extracted and gene expression profile was determined using cDNA arrays. Phosphoprotein profiles were studied in parallel.
Delayed freezing affected the RNA quality only in 3 samples, which were not subjected to gene profiling. In the 8 breast cancer cases with apparently intact RNA also in sample aliquots frozen at delayed times, 461 genes were modulated simply as a function of freezing timing. Some of these genes were included in gene signatures biologically and clinically relevant for breast cancer. Delayed freezing also affected detection of phosphoproteins, whose pattern may be crucial for clinical decision on target-directed drugs.
Time elapsed between surgery and freezing of samples appears to have a strong impact and should be considered as a mandatory variable to control for clinical implications of inadequate tissue handling.
Our understanding of the underlying molecular mechanisms in various human tumors has increased exponentially over the last decades due to the rapid development and application of technologies such as DNA microarrays and mass spectrometry-based proteomics. DNA microarray technology has markedly contributed to the comprehension of the complexity of pathways governing aggressiveness and treatment response of human neoplasias . Further developments are expected as techniques are improving and allow the use of tiny amounts of tissue both frozen or even fixed and paraffin-embedded for comprehensive molecular analyses .
When comparing results from published microarray studies, differences in patient cohorts, treatment regimens, type of gene expression platform employed are usually taken into account, while procedures and timing related to the processes encompassed between surgical excision and freezing and/or fixation of the biological samples are poorly controlled. Such procedures, applied during sample handling may however significantly affect microarray data. In particular ischemia combined with room temperature storage due to the prolongation of the time elapsed between surgical removal and snap-freezing in liquid nitrogen is likely to alter gene expression patterns  as well as protein expression . If this is the case the gene expression data may be modified by an external source of variability, and consequently represent the result of a complicated interplay between disease-associated gene and conditions of sample handling rather than a specific disease condition.
Despite the definition of strict operating procedures for collection of samples in tissue banks , pre-analytical procedures have been scarcely ever controlled during the daily routine. Such pre-analytical variation is probably not likely to impact results from comprehensive genome-wide profiling studies designed to select or discover genes linked to a particular pathological condition. In fact when employing whole genome arrays the pre-analytical noise may be compensated by the large number of investigated transcripts. However, in the case of validation of signatures or even more in the case of their use for clinical decision, according to FDA-approved commercially available tests as the OncoDx™ (Genomic Health, Redwood City, CA) and the MammaPrint® (Agendia, Netherlands), it is very important to try to build gene signatures containing robust genes not affected by handling procedures and therefore to define which are the genes particularly prone to be modified by inadequate pre-analytical processing. Indeed the effect of inappropriate tissue handling is a critical issue not only for frozen samples, but also for fixed samples where the elapsed time between surgical removal and fixation adds technical variability to the possible alterations induced by fixation procedure.
Some studies have already addressed the issue in a number of human, rat and mouse tissues. Using real-time RT-PCR quantification in mouse liver specimens, Almeida et al  assessed the expression of six genes and showed their modulation under ischemic conditions both at two different temperatures mimicking surgical ischemic conditions and at room temperature waiting time prior to pathological examination. Similarly using cDNA microarrays three separate groups, Huang , Blackhall  and Dash , analyzed respectively specimens from a human colon normal mucosa sample, a couple of lung tumors and four prostate samples. All these studies disclosed differential gene expression patterns related to delays in tissue processing. Miyatake et al  drew the same conclusion investigating the effects of ischemia in different tissues of rat (lung, liver, kidney and spleen). In addition, the authors demonstrated a tissue-dependent transcriptional response against warm ischemia stressing this way the importance of homogeneous surgical procedures for result comparison.
Breast tumors have been extensively profiled by gene expression analysis and encouraging data are available both in the area of prognosis and treatment response prediction . In order to gain awareness of the gene expression artifacts linked to improper tissue handling also for this particular type of neoplasia, we have investigated gene expression profiles in a set of primary breast tumors, grossly subdivided into 4 aliquots and kept at room temperature for 0, 2, 6 and 24 hours after resection before snap-freezing in liquid nitrogen. In addition, we concomitantly explored the time-related preservation of tyrosine phosphorylated proteins, considered to be crucial in cancer development and directly implicated on target-directed therapeutics. Our results show that elapsed time between surgery and freezing has a strong effect on functional genomic analysis and may affect the expression of genes relevant for predictive signatures as well as the phosphorylation status of molecular targets. Therefore pre-analytical conditions should be taken into account as a separate variable both when designing and interpreting comprehensive molecular analyses on archival specimens.
Eleven histologically confirmed primary breast tumors were obtained from the Tissue Bank of the Fondazione IRCCS Istituto Nazionale Tumori between October 2006 and March 2007. For each sample a written informed consent signed by the patient was available authorizing the use of leftover material for research purposes. The study was approved by Independent Ethics Committee and Institutional Review Board. Samples used in the study were selected based on the size, carefully avoiding necrotic areas, fat and normal tissue. An adjacent section was stained and used for defining the percentage of tumor cells. Only specimens with more than 70% of tumor cells were included in the study. Each tissue sample was divided into 4 aliquots. One was immediately frozen while the remaining three were frozen after 2, 6 and 24 hours at room temperature.
RNA isolation and expression profiling
Tissue was pulverized using a Mikrodismembrator (Braun Biotech International, Germany). Total RNA was extracted with the Trizol reagent (Invitrogen, Carlsbad, CA) according to manufacture instructions and an additional DNase digestion was performed using the RNeasy kit (Qiagen, Valencia, CA). After each extraction a small fraction of RNA was used for quality and yield assessment. RNA total concentration and purity were determined by UV spectrometry. Total RNA electrophoretic profile was analysed by the Agilent RNA 6000 NanoLabChip kit on the Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA) using the software provided by the manufacturer for determination of RIN (RNA integrity number).
Probe labelling and hybridization were performed as previously described . The samples and a reference RNA (Universal Human Reference RNA, Stratagene, La Jolla, CA) were labelled directly with Cy3-dCTP (reference RNA) or Cy5-dCTP (sample RNA) (Amersham Biosciences, Little Chalfont, UK) and indirectly with 3DNA Submicro Expression Array Detection kit (Genisphere, Montvala, NJ). Hybridization was carried out in a hybridization station (Genomic Solutions, Ann Arbour, MI), slides were scanned using the GenePix 4000B microarray scanner and quantified using GenePix Pro 18.104.22.168 (Axon Instruments/Molecular Devices). The RNAs were hybridized on two different cDNA microarrays containing a total of 17172 unique clones selected from the Human sequence-verified I.M.A.G.E. clone collection (Research Genetics/Invitrogen, Carlsbad, CA) and spotted in triplicate. Raw gene expression data have been deposited at the European Bioinformatics Institute (EBI) ArrayExpress and are accessible through accession no. E-MEXP-2035.
Quantitative real-time PCR
Total RNA isolated for the microarray analysis was used to verify the quantity of specific messengers by real-time PCR for 4 biologically relevant differentially expressed genes (i.e., FGFR4, ESR1, ERBB2, FBLN2). RNA was reverse transcribed using the High-Capacity cDNA Archive Kit (Applied Biosystems, Foster City, CA). Samples were amplified in multiplex PCR reactions using one of the assays of interest labelled with FAM. PPIA was used as housekeeping gene. Reactions were performed in a final reaction volume of 20 μl containing cDNA template, 10 μl 2× TaqMan Gene Expression Assay (Applied Biosystems, Foster City, CA) and thermal cycling was performed on an ABI PRISM 7700 Sequencer Detector (Applied Biosystems, Foster City, CA).
Protein Profile Analysis
Following tissue pulverization, tumor samples were processed as described in . Briefly, samples were lysed in an ice-cold buffer containing 50 mM HEPES (pH 7.6), 150 mM NaCl, 10% glycerol, 1% Triton X-100, 1.5 mM MgCl2, 1 mM EGTA, 10 mM Na4P2O7, 100 mM NaF, protease inhibitors and in presence or absence of 1 mM sodium orthovanadate. After 30 min incubation with gentle rocking at 4°C, lysates were cleared by centrifugation (10 min at 13,000 rpm). Supernatants were collected and, after protein concentration determination, tested for the presence of specific proteins by Western blotting. For Western blotting, cell lysates were resolved by 4-12% SDS-PAGE (precast gel NuPAGE, Invitrogen). Proteins were transferred to nitrocellulose membranes (Hybond-C Super, Amersham Bioscience, Little Chalfont, UK), checked for equal sample loading by Red Ponceau S staining, and probed with the appropriate antibodies. Immunoreactive bands were visualized using horseradish peroxidase-linked secondary antiserum (Amersham Biosciences, Little Chalfont, UK) and detected using an enhanced ECL system (Amersham Biosciences, Little Chalfont, UK). For Western blotting, the following antibodies were used: anti-p-Tyr 4G10 (Upstate Biotechnology, Millipore Corporation, Billerica, MA); p-Neu (Tyr 1248)-R (Santa Cruz, Biothecnology, Santa Cruz, CA) for ERBB2 and anti-FAK (Transduction Laboratories, BD, Franklin Lakes, NJ).
Data analysis and Statistics
Raw microarray images and quantifications were stored and processed in Bio Array Software Environment (BASE, Lund Sweden) . Poor signal quality of background-corrected Cy3 and Cy5 intensities were flagged and a lowess normalization  was applied to each slide. Replicated spots were averaged and their log (base2) expression ratios (tumor/reference) were downloaded from BASE and imported into BRB-ArrayTools version 3.4.0. (Bethesda, MD) developed by Richard Simon and Amy Peng Lam . ANOVA analysis was performed following their instructions. All I.M.A.G.E clones annotations were updated with the latest release of NCBI Unigene (Build No. 194) using SOURCE. In order to find the overlap between our list of genes and genes included in public available signatures, we employed the Merge function of the web-application MatchMiner . Hierarchical clustering was performed using centered correlation and average linkage method.
For gene function analyses, a score obtained calculating the -Log(p-value) of the Fisher Exact Test indicates the probability that a function is obtained by chance; scores equal or greater than 3 give a 99.9% confidence level of not being generated by chance.
The genes identified by microarray analyses were challenged for gene function. After having imported their accession numbers into the Ingenuity Pathway Analysis (IPA) software (Ingenuity Systems, Mountain View, CA), they were categorized according to the gene functions present in the database.
Data analysis for real-time PCR was done using the Sequence Detector version 1.9 software. Relative log-expression of the genes (-ΔCt) was obtained subtracting the number of cycle threshold observed for the 18S gene from that observed for the gene of interest. After being exported further data analysis was performed in a BRB_ArrayToolv3.4.0.
The association between gene expression evaluated by cDNA microarry and real-time PCR was evaluated by linear regression analysis.
Results and Discussion
Gene expression analysis
Clinico-pathological characteristics of samples
Notwithstanding the relatively limited RNA degradation, a time-dependent modulation of genes was observed in the different serial aliquots, likely expected as a cellular reaction to stress conditions following the surgical resection of blood vessels and the suboptimal pre-analytical conditions. The gene expression profile of these serial specimens was performed using cDNA arrays containing 17,172 unique clones and analyzed by the ANOVA model. Using such approaches, genes whose modulation varied over the time independently of the intrinsic patient-specific differences were disclosed. A separate ANOVA model at significant threshold of p < 0.01 was fitted for each gene to relate the log-ratios expression to the factors under investigation (i.e. the time points and the different patients, the latter considered as blocking factor). Our results pointed out a significant variation as a function of time. Considering the number of genes modulated at each single time point compared to the sample immediately frozen as baseline, we noticed an increase in the number of modulated genes especially after 24 hours: 121 genes (0.76% of genes present in our chip) were modulated after 2 hours (66/121 were up-regulated); 157 genes (0.98%) were modulated after 6 hours (110/157 were up-regulated), 657 genes (4.1%) were modulated after 24 hours (287/657 were up-regulated) following tissue resection. Moreover, taking into account the overall time course of gene expression, we found 461 genes regulated at the significant threshold of p < 0.01, representing 2.88% of genes present in our chip. Around 1.03% of 461 genes were expected to be false positives according to the multiple testing correction method of Benjamini-Hochberg .
This result is in keeping with published data [3, 4, 6–9] in which a certain degree of gene expression variation as a function of collection procedures is reported. However, the modulation we noticed in breast cancer specimens was not as remarkable as the one reported by in Spruessel et al , where 20% of all detectable genes are modulated after 30 minutes of colon resection. In our study the majority of the genes (about two-thirds) were down-regulated, and up-regulated genes represented less than 1% of all the genes in our chip. Similarly Dash et al reported less than 0.6% of over-expressed genes in their study on radical prostatectomy specimens .
To better assess the impact of sample handling on the expression of critical genes, some of the modulated genes with an established biological role in breast cancer were subjected to a technical validation by qPCR. The statistically significant correlations between expression evaluated by cDNA microarray and Real Time PCR for FGFR4 (r = 0.57, p = 0.0008), FBLN2 (r = 0.595, p = 0.0004), ESR1 (r = 0.87, p = 0.0000008), ERBB2 (r = 0.754, p = 0.000002) clearly supports the fact that such a modulation is not imputable to technical artifacts, but represents a real biological phenomenon.
Effect of sample collection procedure on gene pathways
Comparison with published breast cancer microarray studies
The use of gene signatures for predicting prognosis and response to treatment is quickly moving from a research area to the clinical practice. The reliability of such tools as well as they performance are strictly related to accurate gene expression determination which could be affected by sample handling and storage. To better understand to which extent commonly used signatures in breast cancer could be affected by artifacts linked to sample handling we verified how many of our 461 genes were included in such gene lists [2, 28–33].
Genes common with main breast cancer genes signatures
28 genes common with the intrinsec genes list (Sorlie et al.)
Fibroblast growth factor receptor 4
Selenium binding protein 1
Periostin, osteoblast specific factor
Serine (or cysteine) proteinase inhibitor, clade H (heat shock protein 47), member 1, (collagen binding protein 1)
DEAD (Asp-Glu-Ala-Asp) box polypeptide 17
Histone 2, H2be
Tumor suppressor candidate 3
Zinc finger protein 516
Polycomb group ring finger 4
Putative acyl-CoA dehydrogenase
Inositol polyphosphate-4-phosphatase, type II, 105 kDa
Estrogen receptor 1
CAP-GLY domain containing linker protein family, member 4
Hypothetical protein FLJ20481
Ral GEF with PH domain and SH3 binding motif 1
Proteasome (prosome, macropain) subunit, beta type, 10
LPS-responsive vesicle trafficking, beach and anchor containing
Ubiquitin-activating enzyme E1 (A1S9T and BN75 temperature sensitivity complementing)
Collagen, type VI, alpha 3
Transmembrane 9 superfamily member 2
Transporter 1, ATP-binding cassette, sub-family B (MDR/TAP)
Trophoblast-derived noncoding RNA
24 genes common with wound healing signature (Chang et al)
Eukaryotic translation initiation factor 4 gamma, 1
primase, polypeptide 2A, 58 kDa
P450 (cytochrome) oxidoreductase
Fanconi anemia, complementation group G
Exosome component 7
Farnesyl diphosphate synthase (farnesyl pyrophosphate synthetase, dimethylallyltranstransferase, geranyltranstransferase)
Chromosome 5 open reading frame 4
Cell division cycle associated 4
H2A histone family, member V
Fatty acid binding protein 3, muscle and heart (mammary-derived growth inhibitor)
Hypothetical protein MGC4308
Fatty acid desaturase 1
Surfactant, pulmonary-associated protein B
Phospholipase A2 receptor 1, 180 kDa
Selenium binding protein 1
Phenylalanine-tRNA synthetase-like, beta subunit
MADS box transcription enhancer factor 2, polypeptide D (myocyte enhancer factor 2D)
Nuclear transport factor 2
Inhibitor of DNA binding 3, dominant negative helix-loop-helix protein
Inhibitor of DNA binding 2, dominant negative helix-loop-helix protein
Proteoglycan 2, bone marrow
H2A histone family, member X
WD repeat and SOCS box-containing 1
17 genes common with new intrinsic gene list (Zhiyan et al.)
Tripartite motif-containing 33
RAB5C, member RAS oncogene family
primase, polypeptide 2A, 58 kDa
RAB21, member RAS oncogene family
Serine/threonine protein kinase SSTK
ST6 beta-galactosamide alpha-2,6-sialyltranferase 1
Interferon regulatory factor 4
LPS-responsive vesicle trafficking, beach and anchor containing
Interleukin 13 receptor, alpha 1
Estrogen receptor 1
Microfibrillar-associated protein 2
Signal transducer and activator of transcription 1, 91 kDa
Protein tyrosine phosphatase, receptor type, f polypeptide (PTPRF), interacting protein (liprin), alpha 1
Protein tyrosine phosphatase, receptor type, f polypeptide (PTPRF), interacting protein (liprin), alpha 4
RAB31, member RAS oncogene family
Inositol polyphosphate-4-phosphatase, type II, 105 kDa
2 genes common with Recurrence Score genes (Paik et al)
Estrogen receptor 1
Matrix metalloproteinase 11 (stromelysin 3)
We then assumed that the extent of the impact of tissue handling on the prognostic/predictive reliability of a given gene signature, could be indirectly evaluated looking at its capability to maintain the clustering of individual samples collected at different times.
Results are reported in Fig. 4 for various gene signatures. After clustering the sample aliquots with the genes derived from the signature of Chang et al (343 genes in common with our data set) aliquots from 5 of 8 samples clustered together and similar data were obtained using the signature of Sorlie et al (271 genes in common), where in 6 out of 8 samples aliquots collected at different times clustered together. As expected, clustering of samples using signatures with a low number of genes as that of Paik et al , or with a small number of genes after searching for common genes within our data set, as that of Ivshina et al , were instead affected by inappropriate tissue handling as tissue aliquots obtained from the same patient did not cluster together. This underlines the fact that tissue mishandling has a limited effect on signatures with high gene numbers as already suggested in the Background section.
Interestingly, using the signatures with an intermediate number of genes (respectively 27 and 19 common genes) as in the case of the signature by Wang et al  and by van't Veer , we still observed that in most samples (5 out of 8 and 4 out of 8, respectively) aliquots from the same patient obtained at different times, did in fact cluster together. In these dendrograms the relatively low number of genes was probably partially compensated by the lack of overlapping (only 2 overlapping genes with the genes by Wang et al and none with genes by van't Veer) between such genes and the 461 genes that we identified as genes modulated by temperature. This may be an indirect proof of the fact that our 461 genes may be considered as reliable indicators of tissue mishandling
Impact of sample collection procedures on protein profile
Detection of phosphoepitopes is often very important for understanding whether a specific pathway is active, thus inadequate pre-analytical procedure may dramatically modify biological interpretation of data obtained on clinical samples. Therefore sample Pb1 and Pb2, for which ERBB2 immuhistochemistry positivity was detected (Table 1) and gene expression analysis showed no significant change as a function of the time delay in freezing, were tested in Western blotting for the level of Tyr 1248 phosphorylation of ERBB2 using a phospho-specific antibody (Fig 6C). The result was the same for both tumor specimens: signal decreased until disappearing within 24 hours post resection. When an anti-FAK antibody was used to representatively testing protein integrity during time course, it was found that the signal intensity disappeared at 24 hours of incubation (Fig 6C). Still considering Western blotting a non-quantitative assay it is clear that at 24 hours time point protein degradation was maximal.
The great opportunities offered by the development of the microarray technology have fueled many studies in cancer and more are underway exploiting proteomic tools. Unfortunately, when taking into consideration the reliability of results many factors are usually considered, eg. sample size, treatment homogeneity, type of follow up, while pre-analytical conditions are often neglected. In this study we demonstrate that the likelihood to detect transcriptional differences as well as protein patterns reflecting the biology of samples may be affected by a correct tissue handling.
RNA integrity is not a sufficient proof of a good quality sample and even in samples with high RIN gene expression may be affected by delays in freezing and hypoxia. The stress caused by resection of blood vessels leads in fact to rapid cellular changes typically activating hypoxia or apoptotic pathways. Genes included in signatures commonly used in breast cancer also may undergo transcriptional modulations linked to sample handling procedures. If this is the case misleading results may be obtained for single samples and probably evaluation of the performance of a single gene signature may also be influenced. Interestingly a down-regulation was consistently reported for ESR1, a gene playing a major role in clinical management of breast cancer. The latter finding can probably explain why when looking at messenger levels, the so-called ER gene cassettes are more predictive compared to ER itself . Tissue handling induced modulations are not only confined to messengers. In fact in the case of ESR1 itself also protein expression may be affected by prolonged hypoxia, although both down-regulations [36–38] and up-regulations have been observed .
Few studies exist on the stability of proteins and phosphoproteins in tissue after resection and in dependence of preanalytical handling. In this study it was taken into account the dynamic nature of cell protein composition and potentially of cell signaling in post excision setting. Generated data are limited to few specimens and further future in-depth investigation of the fluctuations of protein signaling and post translational modifications is needed. Notwithstanding some phosphoprotein profile fluctuations may potentially be due to tissue individuality, our non-quantitative data are sufficient to support the conclusion that orthovadate addition in lysis buffer is insufficient to block fluctuations in phosphoprotein profiles and also the recommendation to freeze tissue within few minutes of resection to preserve the state of phosphoproteins. In fact, the case of ERBB2 showed the a time-related phosphorylated status in the absence of significant changes at messenger level.
In summary all these data demonstrate that taking advantage of the new profiling techniques for designing tailored treatment strategies is a goal which cannot be achieved without an accurate standardization of tissue collection procedures. Tissue collection procedures should therefore be rigorously defined and controlled to avoid false results which may negatively affect the clinical outcome when attempting personalized treatment strategies. The simple obtainment of good quality RNA samples is not a guarantee for a reliable gene expression profile.
Ingenuity Pathway Analysis
RNA integrity number
Grant support: Associazione Italiana per la Ricerca sul Cancro (I Bongarzone), Cassa di Risparmio delle Province Lombarde Foundation and Regione Lombardia (MA Pierotti).
This work was also supported by funds (MG Daidone) obtained through a law by the Italian government which allowed Italian citizen to allocate the 5 × 1000 share of their tax payment to support a research or charitable institution of their choice. We wish to thank all those citizens who decided to donate their 5 × 1000 to IRCCS Fondazione Istituto Nazionale dei Tumori, Milano.
- Brennan DJ, Kelly C, Rexhepaj E, Dervan PA, Duffy MJ, Gallagher WM: Contribution of DNA and tissue microarray technology to the identification and validation of biomarkers and personalized medicine in breast cancer. Cancer Genomics Proteomics. 2007, 4: 121-134.PubMed
- Paik S, Shak S, Tang G, Kim C, Baker J, Cronin M, Baehner FL, Walker MG, Watson D, Park T, Hiller W, Fisher ER, Wickerham DL, Bryant J, Wolmark N: A multigene assay to predict recurrence of tamoxifen-treated, node-negative breast cancer patients. New Engl J Med. 2004, 351: 2817-2826. 10.1056/NEJMoa041588.View ArticlePubMed
- Almeida A, Paul Thiery J, Magdelénat H, Radvanyi F: Gene expression analysis by real-time reverse transcription polymerase chain reaction: influence of tissue handling. Anala Biochem. 2004, 328: 101-108. 10.1016/j.ab.2004.02.004.View Article
- Spruessel A, Steimann G, Jung M, Lee SA, Carr T, Fentz AK, Spangenberg J, Zornig C, Juhl HH, David KA: Tissue ischemia time affects gene and protein expression patterns within minutes following surgical tumor excision. Biotechniques. 2004, 36: 1030-1037.PubMed
- Mager SR, Oomen MH, Morente MM, Ratcliffe C, Knox K, Kerr DJ, Pezzella F, Riegman PH: Standard operating procedure for the collection of fresh frozen tissue samples. Eur J Cancer. 2007, 43: 828-834. 10.1016/j.ejca.2007.01.002.View ArticlePubMed
- Huang J, Qi R, Quackenbush J, Dauway E, Lazaridis E, Yeatman T: Effects of ischemia on gene expression. J Surg Res. 2001, 99: 222-227. 10.1006/jsre.2001.6195.View ArticlePubMed
- Blackhall FH, Pintilie M, Wigle DA, Jurisica I, Liu N, Radulovich N, Johnston MR, Keshavjee S, Tsao MS: Stability and heterogeneity of expression profiles in lung cancer specimens harvested following surgical resection. Neoplasia. 2004, 6: 761-767. 10.1593/neo.04301.PubMed CentralView ArticlePubMed
- Dash A, Maine IP, Varambally S, Shen R, Chinnaiyan AM, Rubin MA: Changes in differential gene expression because of warm ischemia time of radical prostatectomy specimens. Am J Pathol. 2002, 161: 1743-1748.PubMed CentralView ArticlePubMed
- Miyatake Y, Ikeda H, Michimata R, Koizumi S, Ishizu A, Nishimura N, Yoshiki T: Differential modulation of gene expression among rat tissues with warm ischemia. Exp Mol Pathol. 2004, 77: 222-230. 10.1016/j.yexmp.2004.07.004.View ArticlePubMed
- Haibe-Kains B, Desmedt C, Piette F, Buyse M, Cardoso F, Van't Veer L, Piccart M, Bontempi G, Sotiriou C: Comparison of Prognostic Gene Expression signatures for Breast Cancer. BMC Genomics. 2008, 9: 394-10.1186/1471-2164-9-394.PubMed CentralView ArticlePubMed
- Cappelletti V, Gariboldi M, De Cecco L, Toffanin S, Reid JF, Lusa L, Bajetta E, Celio L, Greco M, Fabbri A, Pierotti MA, Daidone MG: Patterns and changes in gene expression following neo-adjuvant anti-estrogen treatment in estrogen receptor-positive breast cancer. Endocr Relat Cancer. 2008, 15: 439-449. 10.1677/ERC-07-0274.View ArticlePubMed
- Gorla L, Cantu M, Miccichè F, Patelli C, Mondellini P, Pierotti MA, Bongarzone I: RET oncoproteins induce tyrosine phosphorylation changes of proteins involved in RNA metabolism. Cell Signal. 2006, 18: 2272-2282. 10.1016/j.cellsig.2006.05.016.View ArticlePubMed
- Saal LH, Troein C, Vallon-Christersson J, Gruvberger S, Borg A, Peterson C: BioArray Software Environment (BASE): A Platform for Comprehensive Management and Analysis of Microarray Data. Genome Biology. 2002, 3: software0003.1-0003.6. 10.1186/gb-2002-3-8-software0003.View Article
- 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, 3: e15-10.1093/nar/30.4.e15.View Article
- Simon R, Lam A, Li MC, Ngan M, Menenzes S, Zhao Y: Analysis of Gene Expression Data Using BRB-Array Tools. Cancer Informatics. 2007, 3: 11-17.PubMed CentralPubMed
- Diehn M, Sherlock G, Binkley G, Jin H, Metese JC, Hernandez-Boussard T, Rees CA, Cherry JM, Botstein D, Brown PO, Alizadeh AA: SOURCE: a unified genomic resource of functional annotations, ontologies, and gene expression data. Nucleic Acids Res. 2003, 31: 219-223. 10.1093/nar/gkg014.PubMed CentralView ArticlePubMed
- Bussey KJ, Kane D, Sunshine M, Narasimhan S, Nishizuka S, Reinhold WC, Zeeberg B, Ajay W, Weinstein JN: MatchMiner: a tool for batch navigation among gene and gene product identifiers. Genome Biology. 2003, 4: R27-10.1186/gb-2003-4-4-r27.PubMed CentralView ArticlePubMed
- Beekman JM, Reischl J, Henderson D, Bauer D, Ternes R, Peña C, Lathia C, Heubach JF: Recovery of microarray-quality RNA from frozen EDTA blood samples. J Pharmacol Toxicol Methods. 2009, 59: 44-49. 10.1016/j.vascn.2008.10.003.View ArticlePubMed
- Fleige S, Walf V, Huch S, Prgomet C, Sehm J, Pfaffl MW: Comparison of relative mRNA quantification models and the impact of RNA integrity in quantitative real-time RT-PCR. Biotechnol Lett. 2006, 28: 1601-1613. 10.1007/s10529-006-9127-2.View ArticlePubMed
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: A practical and powerful approach to multiple testing. J Roy Stat Soc, Ser B. 1995, 57: 289-300.
- Rommel C, Camps M, Ji H: PI3K delta and PI3K gamma: partners in crime in inflammation in rheumatoid arthritis and beyond?. Nat Rev Immunol. 2007, 3: 191-201. 10.1038/nri2036.View Article
- Lee KS, Park SJ, Kim SR, Min KH, Jin SM, Puri KD, Lee YC: Phosphoinositide 3-kinase-delta inhibitor reduces vascular permeability in a murine model of asthma. J Allergy Clin Immunol. 2006, 118: 403-9. 10.1016/j.jaci.2006.04.041.View ArticlePubMed
- Ivanov SV, Salnikow K, Ivanova AV, Bai L, Lerman MI: Hypoxic repression of STAT1 and its downstream genes by a pVHL/HIF-1 target DEC1/STRA13. Oncogene. 2007, 26: 802-12. 10.1038/sj.onc.1209842.View ArticlePubMed
- Huang Q, Gumireddy K, Schrier M, le Sage C, Nagel R, Nair S, Egan DA, Li A, Huang G, Klein-Szanto AJ, Gimotty PA, Katsaros D, Coukos G, Zhang L, Puré E, Agami R: The microRNAs miR-373 and miR-520c promote tumour invasion and metastasis. Nat Cell Biol. 2008, 10: 202-210. 10.1038/ncb1681.View ArticlePubMed
- Diaz LK, Zhou X, Wright ET, Cristofanilli M, Smith T, Yang Y, Sneige N, Sahin A, Gilcrease MZ: CD44 expression is associated with increased survival in node-negative invasive breast carcinoma. Clin Cancer Res. 2005, 11: 3309-3314. 10.1158/1078-0432.CCR-04-2184.View ArticlePubMed
- Lopez JI, Camenisch TD, Stevens MV, Sands BJ, McDonald J, Schroeder JA: CD44 attenuates metastatic invasion during breast cancer progression. Cancer Res. 2005, 65: 6755-6763. 10.1158/0008-5472.CAN-05-0863.View ArticlePubMed
- Meijer D, Sieuwerts AM, Look MP, van Agthoven T, Foekens JA, Dorssers LC: Fibroblast growth factor receptor 4 predicts failure on tamoxifen therapy in patients with recurrent breast cancer. Endocr Relat Cancer. 2008, 15: 101-111. 10.1677/ERC-07-0080.View ArticlePubMed
- Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, Rijn van de M, Jeffrey SS, Thorsen T, Quist H, Matese JC, Brown PO, Botstein D, Eystein Lønning P, Børresen-Dale AL: Gene expression patterns of breast carcinomas distinguish tumour subclasses with clinical implications. Proc Natl Acad Sci. 2001, 98: 10869-10874. 10.1073/pnas.191367098.PubMed CentralView ArticlePubMed
- Loi S, Haibe-Kains B, Desmedt C, Lallemand F, Tutt AM, Gillet C, Ellis P, Harris A, Bergh J, Foekens JA, Klijn JG, Larsimont D, Buyse M, Bontempi G, Delorenzi M, Piccart MJ, Sotiriou C: Definition of clinically distinct molecular subtypes in estrogen receptor-positive breast carcinomas through genomic grade. J Clin Oncol. 2007, 25: 1239-1246. 10.1200/JCO.2006.07.1522.View ArticlePubMed
- Chang HY, Sneddon JB, Alizadeh AA, Sood R, West RB, Montgomery K, Chi JT, Rijn van de 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.PubMed CentralView ArticlePubMed
- Ivshina AV, George J, Senko O, Mow B, Putti TC, Smeds J, Lindahl T, Pawitan Y, Hall P, Nordgren H, Wong JE, Liu ET, Bergh J, Kuznetsov VA, Miller LD: Genetic reclassification of histologic grade delineates new clinical subtypes of breast cancer. Cancer Res. 2006, 66: 10292-10301. 10.1158/0008-5472.CAN-05-4414.View ArticlePubMed
- van't Veer LJ, Dai H, Vijver van de MJ, He YD, Hart AA, Mao M, Peterse HL, Kooy van der 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 Article
- Hu Z, Fan C, Oh DS, Marron JS, He X, Qaqish BF, Livasy C, Carey LA, Reynolds E, Dressler L, Nobel A, Parker J, Ewend MG, Sawyer LR, Wu J, Liu Y, Nanda R, Tretiakova M, Ruiz Orrico A, Dreher D, Palazzo JP, Perreard L, Nelson E, Mone M, Hansen H, Mullins M, Quackenbush JF, Ellis MJ, Olopade OI, Bernard PS, Perou CM: The molecular portraits of breast tumors are conserved across microarray platforms. BMC Genomics. 2006, 7: 96-10.1186/1471-2164-7-96.PubMed CentralView ArticlePubMed
- Wang Y, Klijn JG, Zhang Y, Sieuwerts AM, Look MP, Yang F, Talantov D, Timmermans M, Meijer-van Gelder ME, Yu J, Jatkoe T, Berns EM, Atkins D, Foekens JA: Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer. Lancet. 2005, 365: 671-679.View ArticlePubMed
- Paik S, Tang G, Shak S, Kim C, Baker J, Kim W: Gene expression and benefit of chemotherapy in women with node-negative, estrogen receptor-positive breast cancer. J Clin Oncol. 2006, 24: 3726-3734. 10.1200/JCO.2005.04.7985.View ArticlePubMed
- Stoner M, Saville B, Wormke M, Dean D, Burghardt R, Safe S: Hypoxia induces proteasome-dependent degradation of estrogen receptor alpha in ZR-75 breast cancer cells. Mol Endocrinol. 2002, 16: 2231-2242. 10.1210/me.2001-0347.View ArticlePubMed
- Cooper C, Liu GY, Niu YL, Santos S, Murphy LC, Watson PH: Intermittent hypoxia induces proteasome-dependent down-regulation of estrogen receptor alpha in human breast carcinoma. Clin Cancer Res. 2004, 10: 8720-8727. 10.1158/1078-0432.CCR-04-1235.View ArticlePubMed
- Helczynska K, Kronblad A, Jogi A, Nilsson E, Beckman S, Landberg G, Palman S: Hypoxia promotes a dedifferentiated phenotype in ductal breast carcinoma in situ. Cancer Res. 2003, 63: 1441-1444.PubMed
- Bos R, Zhong H, Hanrahan CF, Mommers EC, Semenza GL, Pinedo HM, Abeloff MD, Simons JW, van Diest PJ, Wall van der E: Levels of hypoxia-inducible factor-1 alpha during breast carcinogenesis. J Natl Cancer Inst. 2001, 93: 309-314. 10.1093/jnci/93.4.309.View ArticlePubMed
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/9/409/prepub
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.