Gene expression profiling of peripheral blood cells for early detection of breast cancer
Breast Cancer Research volume 12, Article number: R7 (2010)
Early detection of breast cancer is key to successful treatment and patient survival. We have previously reported the potential use of gene expression profiling of peripheral blood cells for early detection of breast cancer. The aim of the present study was to refine these findings using a larger sample size and a commercially available microarray platform.
Blood samples were collected from 121 females referred for diagnostic mammography following an initial suspicious screening mammogram. Diagnostic work-up revealed that 67 of these women had breast cancer while 54 had no malignant disease. Additionally, nine samples from six healthy female controls were included. Gene expression analyses were conducted using high density oligonucleotide microarrays. Partial Least Squares Regression (PLSR) was used for model building while a leave-one-out (LOO) double cross validation approach was used to identify predictors and estimate their prediction efficiency.
A set of 738 probes that discriminated breast cancer and non-breast cancer samples was identified. By cross validation we achieved an estimated prediction accuracy of 79.5% with a sensitivity of 80.6% and a specificity of 78.3%. The genes deregulated in blood of breast cancer patients are related to functional processes such as defense response, translation, and various metabolic processes, such as lipid- and steroid metabolism.
We have identified a gene signature in whole blood that classifies breast cancer patients and healthy women with good accuracy supporting our previous findings.
Cancer of the breast is the most common cancer among women worldwide with an estimated 1,300,000 new cases and 465,000 deaths annually . In Norway, the age-adjusted incidence rate for breast cancer has more than doubled from 36.7 per 100,000 in the period 1953 to 1957 to 75.6 per 100,000 in the period 2003 to 2007 . To reduce breast cancer mortality, early detection and appropriate treatment play a key role . The five-year survival rate for stage I breast cancer in Norway in the period 1998 to 2002 was 95%, and 16.8% for stage IV metastatic breast cancer . This emphasizes the importance of early detection so that treatment can be initiated as early as possible during tumor development. Mammographic screening, physical examination and self examination are the main modalities for breast cancer detection today, but only mammography screening has been shown to reduce mortality. When a tumor is detectable in the breast, either by palpation or mammography, the tumor might have been present for several years and have had the ability to spread to distant organs. The growth rate of breast tumors varies considerably between subjects . Some tumors grow so rapidly that they escape a biannual screening program and hence show clinical symptoms before detection by mammography. In addition, mammographic sensitivity is significantly reduced in women with dense breast tissue, often seen in pre-menopausal women or those receiving menopausal hormone therapy . Due to the low sensitivity of mammography in women with dense breast tissue, other imaging modalities have been introduced in breast cancer screening including ultrasonography and magnetic resonance imaging (MRI). However, ultrasound is very operator-dependent, time-consuming, and is associated with many false positive results. MRI is expensive, and the high false positive rate, limited resources and lack of universally accepted imaging guidelines restrict the use of MRI in a screening setting. The need for improved methods to accurately detect breast cancer at an early stage is highly desirable.
Previous studies have found that use of peripheral blood cells for transcriptome analysis is valuable to assess disease-associated [6–10] and drug-response related gene signatures . We have previously demonstrated the potential use of gene expression profiling of peripheral blood cells for early detection of breast cancer . Blood samples are easily available, minimally invasive and can be collected at low cost making them an attractive alternative modality for diagnostic purposes. The rationale for using blood as a clinical sample is that breast cancer triggers a response in circulating blood cells, leading to a traceable change in the whole blood gene expression signature. In this study we aimed to refine our previous findings  with a different sample set, using a larger sample size and a commercially available microarray platform.
Materials and methods
Subject information and blood sampling for microarray experiments
Two hundred blood samples were collected between 2002 and 2004 at two Norwegian hospitals (Ullevål University Hospital and Haukeland University Hospital) after written informed consent under approval from the Regional Ethical Committee of Norway (Ref. no. 416-01151). The subjects included were randomly selected among women called in for a second look after a first suspect screening mammogram. The samples were collected prior to a clinical examination that includes diagnostic mammography and biopsy or fine needle aspiration in the case of a positive mammographic finding. Cytology revealed whether the findings were of malignant or benign origin. For the subjects with no abnormal mammographic findings, the standard of truth was mammography alone. From each woman, 2.5 ml blood was collected in PAXgene™ tubes (PreAnalytiX, Hombrechtikon, Switzerland) and left overnight at room temperature before storing at -80°C until use. As a result of method development and testing of various gene expression platforms, only 121 of the 200 samples initially collected were included in this study. The diagnostic mammograms and histopathology reports revealed that out of these 121 women, 57 had invasive breast cancer, 10 had ductal carcinoma in situ (DCIS) and 54 had no sign of malignant disease. Of these latter 54, 12 had benign findings including fibroadenomas, cysts and some unspecified findings (Table 1). Regarding the breast cancer subjects, tumor stage, grade and other relevant clinical data were recorded (Tables 1 and 2). The individuals in the case and control groups are balanced in relation to age, menopausal status and previous menopausal hormone therapy (Table 3). In addition to the 121 samples, five blood samples were collected from two healthy women at multiple time points (biological replicates), three blood samples from pregnant women, and one sample from a breast feeding healthy woman were collected, leaving 130 samples from 127 individuals for gene expression analysis (Table 1).
To control for technical variability such as different microarray production batches, lot variations of reagents and kits, day to day variations and effects related to different laboratory operators, a strict experimental design was followed. Samples were randomly divided into batches of 10, containing equal numbers of samples from women with breast cancer and those with no sign of the disease. All samples within each batch were handled together through each experimental step by one operator alone and the operators were blinded to cancer status. Two control samples were included in each batch following the same experimental procedures as the other 10. These control samples were composed of total RNA isolated from one healthy female. The order of the samples within each batch was randomized. In order to correct for any batch variations, we used the batch adjustment method described by Tibshirani . A total of 13 batches including 130 samples and 26 technical controls were thus analyzed.
PAXgene™ tubes were thawed over night in batches of 12 tubes and total RNA was extracted according to the manufacturer's protocol. Total RNA was stored at -80°C prior to analyses. RNA quality and quantity measures were conducted using the 2100 Bioanalyzer (Agilent Technologies, Santa Clara, California, USA) and the NanoDrop ND-1000 spectrophotometer (Thermo Scientific, Wilmington, Delaware, USA) respectively.
Microarray gene expression studies were conducted using single channel Applied Biosystems Human Genome Survey microarrays v2.0 containing 32,878 probes representing 29,098 genes. From each sample, 500 ng total RNA was amplified and labeled according to the NanoAmp RT-IVT Labeling Kit Protocol and hybridized onto the array for 16 hours at 55°C. Following hybridization, slides were manually washed and prepared according to the manufacturer's recommendation before image capturing using the AB1700 reader. Identification and quantification of gene expression signals, signal-to-noise ratios and flagging of failed spots were conducted using the Applied Biosystems Expression System software. Raw data files were exported for further analysis.
Data analysis was performed using R  and tools from the Bioconductor project , adapted to our needs. Data was preprocessed in the following way: data were log2 transformed while individual measurements with signal-to-noise <3 or flag values >8,191 were set as missing. Probes with more than 5% missing values over all 156 arrays were excluded. Preprocessing left 156 samples and 11,217 probes for further analyses. Data were standardized (that is, centered and scaled) and missing values were imputed with k-nearest neighbors imputation  using k = 10. Principal components analysis and ANOVA tests for each gene revealed that there were large batch-effects present in the data. Similar batch effects have previously been reported for the same type of data (Dumeaux V, et al., under revision). Each probe was individually treated for batch effects using a one way ANOVA procedure as described by Tibshirani . The 26 technical control samples were then excluded. For the biological replicates (multiple samples from one subject), signal intensities were averaged for each probe. Thus, 127 arrays, one from each individual remained for analysis. Finally, within-array normalization was conducted by global mean subtraction. The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus  and are accessible through GEO Series accession number [GEO:GSE16443].
Feature selection and classifier construction
The gene expression data served as predictors for predicting a dummy-coded response vector. The response vector was given the value -1 or 1 for each sample depending on it being a healthy control or a breast cancer case, respectively. A new gene expression sample was classified as diseased if the predicted value was larger than zero and as healthy otherwise.
Partial Least Squares Regression (PLSR) [18, 19] with double cross-validation was used to construct and test our classifier. PLSR with leave-one-out cross-validation (LOO-CV) was used in combination with Jackknife testing [20, 21] to select significant probes. In more detail, LOO-CV gives the optimal number of components and a set of regression coefficients associated to each probe and jackknife feature selection is used to select probes with regression coefficients different from 0 (P-value ≤ 0.05). A PLSR model is rebuilt on these significant probes and LOO-CV is again used to select the optimal number of components. Finally, the analysis described above is incorporated in an independent loop of LOO-CV in order to test classifier accuracy .
Functional enrichment analysis and biological interpretation
Reducing significant genes to core subsets is a useful step towards understanding biological mechanisms underlying the gene-set association with the phenotype of interest: a smaller number of genes are easier to understand and facilitate biological insight into disease processes. Global test  was used to identify the core probes most strongly explaining the difference between cases and controls. A Global test gene plot illustrates the influence of each individual probe on the significance result. The number of standard deviation of influence on the global test P-value above the reference line under the null hypothesis is termed the z-score. We identify probes with high z-scores (>2) as the core probes. Global test is not testing any specific null hypothesis. It is simply a useful analytical tool to reduce genes that have previously been found differentially expressed, to a core set, by gradually exploring the association of remaining genes as a set with a phenotype.
To explore functional enrichment and possible biological interactions among the genes identified we used the Database for Annotation, Visualization and Integrated Discovery (DAVID) , Human Experimental/Functional Mapper (HEFalMp)  and Graphle . DAVID is a functional annotation tool able to extract biological information out of a large list of genes, while Graphle is an interactive tool displaying relationships between genes predicted by HEFalMp. HEFalMp predicts interactions between genes based on data integration of a vast number of experimental results publicly available and reduce all findings to a single measurement of relatedness . Genes predicted to relate to each other often have a tendency to be co-regulated or are believed to carry out similar cellular tasks.
Construction and characterization of the 738 classifier
Partial Least Squares Regression (PLSR) was used for model building while a leave-one-out cross validation (LOO-CV) was used to evaluate the use of PLS with LOO-CV and Jackknife testing for feature selection. We observed a high number of latent components necessary in the PLS model (N = 21) to achieve a cross validated minimum in the error rate.
By using PLSR with Jackknife testing on all samples (N = 127), a set of 738 probes were identified as significant for disease classification (P < 0.05) between breast cancer patients and women not having the disease. The 738 probe list predicted cases and controls with an estimated accuracy of 79.5% based on LOO-CV with a sensitivity of 80.6% and specificity of 78.3%. Of the 67 breast cancer samples, 54 were predicted correctly, while 47 of the 60 healthy samples were assigned to the correct class (Figure 1a). When plotting the sensitivity versus 1-specificity in a receiver operating characteristics (ROC) curve (Figure 1b), we observe a good separation of the two groups with an area under curve (AUC) of 0.88. Of note, a permutation test (k = 2,000) of the response variable gave a maximum accuracy of 60.6% and AUC of 0.68 (see Figures S1 and S2 in Additional file 1). To assess whether the results could be further improved by a larger sample size, random balanced subsets of samples were analyzed repeating the classifier building process. Results from the analyses indicate that accuracy could be increased with more samples. These results could in theory have been used to estimate a higher achievable AUC if more samples had been available. Attempts to model the results have however proven unsuccessful, consequently no upper AUC has been estimated (see Additional file 2). Using a Fisher exact test, we analyzed whether any of the clinical characteristics were significantly overrepresented among the subjects incorrectly predicted (Figure 1a). Of the false negatives, 11 out of 13 were samples from women having small lesions (<2 cm) including two DCIS, being significantly overrepresented with P = 0.04. Four out of 10 subjects with DCIS are incorrectly predicted as healthy, although not significantly (P = 0.09) overrepresented among the false negatives. Parameters such as tumor grade, estrogen receptor (ER) status or menopausal status do not seem to affect the prediction of the cases in this study. Samples from pregnant women or from women with benign lesions were not overrepresented among the false positives.
Global test gene plot (see Additional file 3) illustrates the influence of each individual probe in the 738 list on the significance result (P = 0.001). Approximately equal numbers of probes are up-regulated (n = 395) and down-regulated (n = 343) in blood of breast cancer patients, with the median z-score equal to 0.55 (sd 1.70) and 0.84 (sd 2.72) respectively. Z-score filtering (Z >2) left 89 core up-regulated probes and 119 core down-regulated probes. We used the core probes for gene interaction prediction.
Functional enrichment analysis
Using DAVID  functional enrichment of the up- and down-regulated genes were investigated separately. Out of the 738 probes, a significant number were not annotated (n = 143) or had limited biological information and these were removed from the list along with duplicate gene symbols. Four hundred and ninety-three gene symbols from the total list were recognized by DAVID and included in the analysis. As background for the functional enrichment analyses, the 11,217 probes left after preprocessing were used. When analyzing the up-regulated genes alone we identified biological processes such as translation, defense response to bacterium, cellular biosynthetic process and response to external stimulus as enriched with false discovery rate (FDR) below 20% (Table 4), while processes involving various metabolic processes were enriched among the genes that were lower expressed in breast cancer patient compared to healthy controls (Table 5).
Graphle/HEFalMp  was used to predict interactions between the genes within each group. When including hundreds of genes in such analyses, giant hairballs of predicted interactions are generated making the results hard to interpret. To reduce the complexity of the interaction maps we selected only the core genes (z-score >2) from the global test analysis (see Additional file 3). After removing probes without annotation, Graphle recognized 47 of the up-regulated core genes and 95 of the down-regulated core genes and predicted their interactions (Figure 2 and Additional file 4). Further, we submitted only the core genes to DAVID to look at functional enrichment within the core genes of each group in particular (see Additional file 5).
The interaction map for the 47 core up-regulated genes identifies two main networks and many of the genes within each network seem to be connected to each other with high interaction confidence (Figure 2). One cluster includes mainly genes coding for ribosomal proteins, playing different roles in the translation machinery. The other cluster contains among others, genes involved in defense response to bacterium. Ten genes are not connected to either of the clusters using edge filter cutoff 0.648 (interaction confidence). The 95 core down-regulated genes do not appear to be as strongly related to each other (see Additional file 4). We observe one main cluster with genes predicted to relate to each other with edge filter cutoff set to 0.643. Many genes cluster in small, more vague interaction networks. No biological processes were enriched among the 95 genes. Edge weights for the genes with highest relatedness are listed in Additional file 6.
Finally, we compared the 738 gene list to the 37 (29 unique) genes published in our previous study . We applied the global test to our data to see whether the 37 gene set published in the initial study were differentially expressed between cases and controls. Twenty of the 29 unique genes were found in the filtered data of the present study, and this set of genes was not significantly differentially expressed between the cases and controls (see Additional file 7). Only two genes were overlapping between the two gene lists (RPS2 and RPL14), both coding for ribosomal proteins.
The biological signal from breast tumors recapitulated in whole blood does not appear to be very strong, reflected by the high number of latent components necessary in the PLS model. Other methods such as prediction analysis for microarray data (PAM) and support vector machines (SVM) were applied but did not improve classification accuracy (data not shown). Nonetheless, our results indicate that gene expression in whole blood serves as a possible diagnostic tool for early detection of breast cancer. We have identified a gene signature that separates breast cancer patients from healthy women with good accuracy. These results are in agreement with the findings in the pilot study, reporting a prediction accuracy of 82%  although for a different predictor. We use a rather liberal cut-off (P- value < 0.05) in the classifier construction and consider the probe list in biological terms, that is, several genes with moderate changes acting in concert within a pathway. The genes identified seem to reflect a biological response related to breast tumor growth. We also reduced the number of selected probes to a set of core genes more likely to be true positives and observe that similar biological processes are enrichment among the core genes up-regulated in blood of breast cancer patients.
False negatives and false positives
The size of the mammary lesion is the only clinical feature that is significantly overrepresented among the falsely predicted samples. Lesions (including DCIS) with size below 2 cm were found significantly overrepresented among the false negatives. It is reasonable that a lower tumor burden will give a weaker response in blood affecting the prediction efficacy.
In our previous study all three pregnant subjects included were predicted as having breast cancer. In this study only one of the samples from the three pregnant women are predicted as having breast cancer.
Since mammography is the standard of truth, we can not exclude the possibility that some of the false positives have very early stage breast cancer or other occult tumors not detectable by existing technology. Follow-up data of these women are unavailable so we can not verify or falsify such a hypothesis.
It is known that growing tumors communicate with the tissue in which they thrive, and also with the cells of the immune system of the host. The high rate of spontaneous occurring tumors in immunocompromised animals  and humans  reflects the inhibitory role of the immune system on tumor growth. The blood-tumor dialogue involves a broad spectrum of signaling molecules and such active cellular crosstalk seems to be reflected in the molecular blood signature of breast cancer patients discussed below.
A cancer-related gene expression signature in whole blood might reflect this communication. An increase or decrease of certain blood cell populations and their activities as a response to the tumor growth may also contribute to the observed difference.
Four biological processes are enriched with FDR below 20% when analyzing the genes up-regulated in blood of breast cancer patients (n = 243), including translation (GO:0006412), defense response to bacterium (GO:0042742), cellular biosynthetic process (GO:0044249) and response to external stimuli (GO:0009605). Among the genes down-regulated we identify processes involving lipid-, steroid-, catecholamine- and phenol metabolism (GO:0044255, 0008202, 0006629, 0006584, 0018958) as enriched.
Translation is a ribosome-mediated process where messenger RNAs (mRNAs) are translated into proteins. Translation is a process taking place in all cells, and it is difficult to draw any firm conclusions from this finding. However, in the pilot study we observed reduced expression of transcripts involved in protein synthesis among the breast cancer patients .
A defense related response observed in breast cancer patients is in agreement with our previous findings . The five genes involved in defense response to bacterium are DEFA3, LTF, CAMP, PPBP and S100A12, genes that all are either highly expressed in neutrophil granulocytes or activators of such. Neutrophil granulocytes are the most abundant type of leukocytes (approximately 60%), whose role is to recognize and kill microorganisms, but also tumor cells . Increased number of neutrophils (neutrophilia) is a sign of acute bacterial infection, but has also been reported in cancer patients, along with reduced lymphocyte counts (lymphocytopenia), referred to as an elevated neutrophil-lymphocyte ratio [30, 31]. Whether such a shift in blood cell populations is due to defense related mechanisms or as a response to tumor derived signals is still not well understood. It has been proposed that tumor cells can attract neutrophils by secreting interleukin 8 (IL8) and that the neutrophils, in a similar manner as in wounds, enhance angiogenesis, tumor growth and progression, and finally cell migration through the ECM . In contrast, one of the genes secreted by neutrophils; lactotransferrin (LTF) has been shown to have an inhibitory effect on tumor growth and metastasis via regulation of natural killer (NK) cell activity, modulation of expression of G1 proteins, inhibition of angiogenesis and enhancement of apoptosis [33, 34]. Interestingly, the gene cystatin A (CSTA), a cystein proteinase inhibitor, which is among the 49 core up-regulated genes has been proposed as a prognostic marker for breast cancer [35, 36]. Elevated lipocalin 2 (LCN2) levels has also been reported in tissue- and urine samples from patients with invasive breast cancer  and is proposed as a noninvasive biomarker for advanced breast cancer. It is believed that LCN2 promotes breast cancer progression by inducing epithelial to mesenchymal transition (EMT) and by increasing cell motility and invasiveness through down-regulation of E-cadherin.
Enrichment of genes involved in various metabolic processes among down-regulated genes suggests a change in the metabolism of breast cancer patients. Tumor growth often leads to dramatic metabolic changes in the host . Several studies have shown altered systemic lipid metabolism in cancer patients , often leading to cachexia. Although cancer cachexia is most common in patients with terminal malignancies, it has also been observed in patients with a relatively small tumor burden . The deregulation of lipid metabolism between cases and controls might reflect an early shift in the metabolism of the tumor bearer.
The gene interaction prediction analyses conducted using Graphle indicates that many of the core up-regulated genes seem to be linked to each other (Figure 2). When looking at the functional enrichment of the core up-regulated genes separately (see Additional file 5), we identify defense response to bacterium as the most significant process. This indicates that the core up-regulated genes carry much of the biological information that seems relevant in a blood-tumor dialogue context discussed above. We also identify taxis (GO:0042330, 0006935) as enriched among the core up-regulated genes alone. Taxis refers to movement of cells in response to external stimulus, possibly reflecting the movement of immune cells towards the growing tumor.
The signature identified in this study is being further refined to improve the diagnostic accuracy. A TaqMan based clinical test, BCtect®  has been developed in part based on the results from this study. This tool could constitute a fast and painless supplement to existing diagnostic technology, and offer a breast cancer test in areas where mammography screening is insufficient.
area under curve
cathelicidin antimicrobial peptide
Database for Annotation, Visualization and Integrated Discovery
ductal carcinoma in situ
Defensin, alpha 3, neutrophil-specific
epithelial to mesenchymal transition
false discovery rate
Human Experimental/Functional Mapper
invasive ductal carcinoma
invasive lobular carcinoma
magnetic resonance imaging
messenger ribonucleic acids
- NK cells:
natural killer cells
Partial Least Squares Regression
prediction analysis for microarray data
pro-platelet basic protein (chemokine (C-X-C motif) ligand 7)
receiver operating characteristics
S100 calcium binding protein A12
support vector machines
Kamangar F, Dores GM, Anderson WF: Patterns of cancer incidence, mortality, and prevalence across five continents: defining priorities to reduce cancer disparities in different geographic regions of the world. J Clin Oncol. 2006, 24: 2137-2150. 10.1200/JCO.2005.05.2308.
Cancer Registry of Norway. Cancer in Norway 2007: Cancer incidence, mortality, survival and prevalence in Norway. Bray, F. 2008, Ref Type 2007: Report
Chu KC, Tarone RE, Kessler LG, Ries LA, Hankey BF, Miller BA, Edwards BK: Recent trends in U.S. breast cancer incidence, survival, and mortality rates. J Natl Cancer Inst. 1996, 88: 1571-1579. 10.1093/jnci/88.21.1571.
Weedon-Fekjaer H, Lindqvist BH, Vatten LJ, Aalen OO, Tretli S: Breast cancer tumor growth estimated through mammography screening data. Breast Cancer Res. 2008, 10: R41-10.1186/bcr2092.
Kolb TM, Lichy J, Newhouse JH: Comparison of the performance of screening mammography, physical examination, and breast US and evaluation of factors that influence them: an analysis of 27,825 patient evaluations. Radiology. 2002, 225: 165-175. 10.1148/radiol.2251011667.
Achiron A, Gurevich M, Friedman N, Kaminski N, Mandel M: Blood transcriptional signatures of multiple sclerosis: unique gene expression of disease activity. Ann Neurol. 2004, 55: 410-417. 10.1002/ana.20008.
Gladkevich A, Kauffman HF, Korf J: Lymphocytes as a neural probe: potential for studying psychiatric disorders. Prog Neuropsychopharmacol Biol Psychiatry. 2004, 28: 559-576. 10.1016/j.pnpbp.2004.01.009.
Tang Y, Nee AC, Lu A, Ran R, Sharp FR: Blood genomic expression profile for neuronal injury. J Cereb Blood Flow Metab. 2003, 23: 310-319. 10.1097/00004647-200303000-00006.
Twine NC, Stover JA, Marshall B, Dukart G, Hidalgo M, Stadler W, Logan T, Dutcher J, Hudes G, Dorner AJ, Slonim DK, Trepicchio WL, Burczynski ME: Disease-associated expression profiles in peripheral blood mononuclear cells from patients with advanced renal cell carcinoma. Cancer Res. 2003, 63: 6069-6075.
Han M, Liew CT, Zhang HW, Chao S, Zheng R, Yip KT, Song ZY, Li HM, Geng XP, Zhu LX, Lin JJ, Marshall KW, Liew CC: Novel blood-based, five-gene biomarker set for the detection of colorectal cancer. Clin Cancer Res. 2008, 14: 455-460. 10.1158/1078-0432.CCR-07-1801.
Burczynski ME, Twine NC, Dukart G, Marshall B, Hidalgo M, Stadler WM, Logan T, Dutcher J, Hudes G, Trepicchio WL, Strahs A, Immermann F, Slonim DK, Dorner AJ: Transcriptional profiles in peripheral blood mononuclear cells prognostic of clinical outcomes in patients with advanced renal cell carcinoma. Clin Cancer Res. 2005, 11: 1181-1189.
Sharma P, Sahni NS, Tibshirani R, Skaane P, Urdal P, Berghagen H, Jensen M, Kristiansen L, Moen C, Sharma P, Zaka A, Arnes J, Sauer T, Akslen LA, Schlichting E, Borresen-Dale AL, Lonneborg A: Early detection of breast cancer based on gene-expression patterns in peripheral blood cells. Breast Cancer Res. 2005, 7: R634-R644. 10.1186/bcr1203.
Tibshirani R, Hastie T, Narasimhan B, Chu G: Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc Natl Acad Sci USA. 2002, 99: 6567-6572. 10.1073/pnas.082099299.
R Development Core Team R: A Language and Environment for Statistical Computing. 2009, Ref Type: Computer Program
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.
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: 520-525. 10.1093/bioinformatics/17.6.520.
Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.
Nguyen DV, Rocke DM: Partial least squares proportional hazard regression for application to DNA microarray survival data. Bioinformatics. 2002, 18: 1625-1632. 10.1093/bioinformatics/18.12.1625.
Wold H: Estimation of principal components and related models by iterative least squares. Multivariate Analysis. Edited by: Krishnaiah PR. 1966, New York: Academic Press, 391-420.
Gidskehaug L, Anderssen E, Flatberg A, Alsberg BK: A framework for significance analysis of gene expression data using dimension reduction methods. BMC Bioinformatics. 2007, 8: 346-10.1186/1471-2105-8-346.
Wu CFJ: Jackknife, bootstrap and other resampling plans in regression analysis. The Annals of Statistics. 1986, 14: 1261-1350. 10.1214/aos/1176350142.
Varma S, Simon R: Bias in error estimation when using cross-validation for model selection. BMC Bioinformatics. 2006, 7: 91-10.1186/1471-2105-7-91.
Goeman JJ, Geer van de SA, de Kort F, van Houwelingen HC: A global test for groups of genes: testing association with a clinical outcome. Bioinformatics. 2004, 20: 93-99. 10.1093/bioinformatics/btg382.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for annotation, visualization, and integrated discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.
Huttenhower C, Haley EM, Hibbs MA, Dumeaux V, Barrett DR, Coller HA, Troyanskaya OG: Exploring the human genome with functional maps. Genome Res. 2009, 19: 1093-1106. 10.1101/gr.082214.108.
Graphle. 2009, Ref Type: Computer Program, [http://sonorus.princeton.edu/graphle/]
Dunn GP, Bruce AT, Ikeda H, Old LJ, Schreiber RD: Cancer immunoediting: from immunosurveillance to tumor escape. Nat Immunol. 2002, 3: 991-998. 10.1038/ni1102-991.
Jain A, Patil VP, Fung J: Incidence of de novo cancer and lymphoproliferative disorders after liver transplantation in relation to age and duration of follow-up. Liver Transpl. 2008, 14: 1406-1411. 10.1002/lt.21609.
Koga Y, Matsuzaki A, Suminoe A, Hattori H, Hara T: Neutrophil-derived TNF-related apoptosis-inducing ligand (TRAIL): a novel mechanism of antitumor effect by neutrophils. Cancer Res. 2004, 64: 1037-1043. 10.1158/0008-5472.CAN-03-1808.
Cho H, Hur HW, Kim SW, Kim SH, Kim JH, Kim YT, Lee K: Pre-treatment neutrophil to lymphocyte ratio is elevated in epithelial ovarian cancer and predicts survival after treatment. Cancer Immunol Immunother. 2008, 58: 15-23. 10.1007/s00262-008-0516-3.
Walsh SR, Cook EJ, Goulder F, Justin TA, Keeling NJ: Neutrophil-lymphocyte ratio as a prognostic factor in colorectal cancer. J Surg Oncol. 2005, 91: 181-184. 10.1002/jso.20329.
De Larco JE, Wuertz BR, Furcht LT: The potential role of neutrophils in promoting the metastatic phenotype of tumors releasing interleukin-8. Clin Cancer Res. 2004, 10: 4895-4900. 10.1158/1078-0432.CCR-03-0760.
Bezault J, Bhimani R, Wiprovnick J, Furmanski P: Human lactoferrin inhibits growth of solid tumors and development of experimental metastases in mice. Cancer Res. 1994, 54: 2310-2312.
Rodrigues L, Teixeira J, Schmitt F, Paulsson M, Mansson HL: Lactoferrin and cancer disease prevention. Crit Rev Food Sci Nutr. 2009, 49: 203-217. 10.1080/10408390701856157.
Kuopio T, Kankaanranta A, Jalava P, Kronqvist P, Kotkansalo T, Weber E, Collan Y: Cysteine proteinase inhibitor cystatin A in breast cancer. Cancer Res. 1998, 58: 432-436.
Levicar N, Kos J, Blejec A, Golouh R, Vrhovec I, Frkovic-Grazio S, Lah TT: Comparison of potential biological markers cathepsin B, cathepsin L, stefin A and stefin B with urokinase and plasminogen activator inhibitor-1 and clinicopathological data of breast carcinoma patients. Cancer Detect Prev. 2002, 26: 42-49. 10.1016/S0361-090X(02)00015-6.
Yang J, Bielenberg DR, Rodig SJ, Doiron R, Clifton MC, Kung AL, Strong RK, Zurakowski D, Moses MA: Lipocalin 2 promotes breast cancer progression. Proc Natl Acad Sci USA. 2009, 106: 3913-3918. 10.1073/pnas.0810617106.
Tisdale MJ: Cachexia in cancer patients. Nat Rev Cancer. 2002, 2: 862-871. 10.1038/nrc927.
Gercel-Taylor C, Doering DL, Kraemer FB, Taylor DD: Aberrations in normal systemic lipid metabolism in ovarian cancer patients. Gynecol Oncol. 1996, 60: 35-41. 10.1006/gyno.1996.0008.
McAndrew PF: Fat metabolism and cancer. Surg Clin North Am. 1986, 66: 1003-1012.
Mackay J, Burford A, Sajun S, Pithers A: Clinical implications of gene expression profiling in cancer. Asia Pacific J of Oncology and Hematology. 2009, 1: 13-17.
This study was supported by the Functional Genomics (FUGE) program from the Norwegian Research Council (NFR-FUGE 159188/S10).
We thank Ole-Christian Lingjære, Einar Rødland and Robert Tibshirani for critically reviewing the statistics section and Simen Myhre for extracting RNA from all blood samples used in this study.
Contributors acknowledged were funded by the Norwegian Research Council (OCL and ER), National Science Foundation and National Institutes of Health (RT) and Oslo University Hospital (SM).
Torbjørn Lindahl, Derek Tobin, Nina Hagen, Anders Lönneborg and Praveen Sharma are employed by DiaGenic ASA and receive their salaries from the company. DiaGenic ASA develops diagnostic products for early detection of various diseases (including breast cancer). Anders Lönneborg and Praveen Sharma are co-founders of DiaGenic ASA and have substantial stocks in the company. Torbjørn Lindahl, Derek Tobin, Nina Hagen also have stocks/options in the company. Anders Lönneborg and Praveen Sharma are inventors of a gene-expression based method to detect disease, conditions or stages thereof (including breast cancer) using samples obtained from an area distant to the site of the disease (including peripheral blood). They protected the method by filing a patent in 1997. The company now holds several patents, and several applications are in process which combined covers the commercial use of the method and the products. The results disclosed in the present work are covered by DiaGenic's patent portfolio. The other authors declare that they have no competing interests.
JA carried out the laboratory work (together with NH), participated in the discussion of the analyses, carried out the functional analyses, prepared the majority of the figures and wrote the manuscript. TL carried out the statistical analyses (together with SS) and wrote the manuscript sections concerning the statistical analyses. VD carried out the global test analysis, participated in the data pre-processing, prepared figures, participated in the discussion of the analyses and critical revision of the manuscript. SS carried out the statistical analyses (together with TL) and participated in the writing of the sections concerning the statistical analyses. DT participated in the discussion and critical revision of the manuscript. NH carried out the laboratory work (together with JA). PSk was responsible for collection of the blood samples, provided clinical information and critical revision of the manuscript. AL conceived and coordinated the study (together with PS and ALBD) and critical reading of the manuscript. PSh conceived and coordinated the study (together with AL and ALBD), participated in the discussion of the statistical analyses and critical reading of the manuscript. ALBD conceived and coordinated the study (together with AL and PS), and participated in the discussion of all analyses and critical revision of the manuscript.
Electronic supplementary material
Additional file 1: Contains two figures, presenting the results from permutation tests (k = 2,000) of the response variables. Figure S1 shows a histogram of permuted accuracy values. The red line indicates the result presented in this study (79.5%) and is evidently better than that achieved by chance. Figure S2 shows a histogram of permuted AUC values. The red line indicates the result presented in this study (0.88) and is evidently better than that achieved by chance. (PDF 304 KB)
Additional file 2: A figure presenting a learning curve - AUC improvement with increasing sample size. The figure shows the prediction accuracy with random balanced sample subsets, using an increasing number of samples and repeating the classifier building and testing process. The blue line indicates the mean AUC, while the light blue lines indicate one standard deviation from the mean. The red dot indicates result reported in this study. Extrapolation of the results does not indicate that the upper limit has been reached. The variance of the AUC decreases with higher percentages, this is an expected result from using more samples to validate the classifier. (PDF 315 KB)
Additional file 3: A figure showing a ranked view of the 738 probes and their influence on the global test P-value. Probes with green bars show higher expression in blood of controls, while probes with red bars show higher expression in blood from women having breast cancer. The blue line indicates the influence of each probe on the global test P-value under the null hypothesis of no association. Black horizontal lines indicate one standard deviation of influence on the global test p-value above the reference line under the null hypothesis. The number of standard deviations is termed the z-score. Probes with high z-scores are the ones that most strongly explain the differences between cases and controls. The 208 core probes (z >2) are highlighted to the left. (PDF 177 KB)
Additional file 4: A figure showing the biological network prediction of the 95 core down-regulated genes in blood of breast cancer patients compared to controls, using edge weight cutoff 0.643 (interaction confidence). (PDF 2 MB)
Additional file 5: A table listing functional enrichment of core genes up-regulated in blood of breast cancer patients compared to healthy subjects. No biological processes were enriched among the 95 core genes down-regulated blood of breast cancer patients compared to healthy subjects. (XLS 14 KB)
Additional file 6: A table listing the interaction confidence predicted by HEFalMp/Graphle between the core genes in each group (z-score >2). (XLS 38 KB)
Additional file 7: A figure showing the influence of twenty of the annotated genes (some represented by multiple probes) from the 37 gene list published in the pilot study on the global test P-value in the present dataset. As illustrated by this plot, the enrichment of this set of 20 genes was not significant in relation to disease status in the present study. Only two of these genes are common with the 738 candidate gene identified; RPL14 and RPS2 (purple). (PDF 110 KB)
About this article
Cite this article
Aarøe, J., Lindahl, T., Dumeaux, V. et al. Gene expression profiling of peripheral blood cells for early detection of breast cancer. Breast Cancer Res 12, R7 (2010). https://doi.org/10.1186/bcr2472
- Breast Cancer
- Breast Cancer Patient
- Area Under Curve
- Partial Little Square Regression
- Core Gene