Abstract
Clustering analysis of microarray data is often criticized for giving ambiguous results because of sensitivity to data perturbation or clustering techniques used. In this paper, we describe a new method based on principal component analysis and ensemble consensus clustering that avoids these problems. We illustrate the method on a public microarray dataset from 36 breast cancer patients of whom 31 were diagnosed with at least two of three pathological stages of disease (atypical ductal hyperplasia (ADH), ductal carcinoma in situ (DCIS) and invasive ductal carcinoma (IDC). Our method identifies an optimum set of genes and divides the samples into stable clusters which correlate with clinical classification into Luminal, Basal-like and Her2+ subtypes. Our analysis reveals a hierarchical portrait of breast cancer progression and identifies genes and pathways for each stage, grade and subtype. An intriguing observation is that the disease phenotype is distinguishable in ADH and progresses along distinct pathways for each subtype. The genetic signature for disease heterogeneity across subtypes is greater than the heterogeneity of progression from DCIS to IDC within a subtype, suggesting that the disease subtypes have distinct progression pathways. Our method identifies six disease subtype and one normal clusters. The first split separates the normal samples from the cancer samples. Next, the cancer cluster splits into low grade (pathological grades 1 and 2) and high grade (pathological grades 2 and 3) while the normal cluster is unchanged. Further, the low grade cluster splits into two subclusters and the high grade cluster into four. The final six disease clusters are mapped into one Luminal A, three Luminal B, one Basal-like and one Her2+. We confirm that the cancer phenotype can be identified in early stage because the genes altered in this stage progressively alter further as the disease progresses through DCIS into IDC. We identify six subtypes of disease which have distinct genetic signatures and remain separated in the clustering hierarchy. Our findings suggest that the heterogeneity of disease across subtypes is higher than the heterogeneity of the disease progression within a subtype, indicating that the subtypes are in fact distinct diseases.
MeSH Terms
Algorithms
Artificial Intelligence
Biomarkers, Tumor/analysis
Breast Neoplasms/diagnosis,metabolism
Carcinoma, Ductal/diagnosis,metabolism
Diagnosis, Computer-Assisted/methods
Disease Progression
Female
Gene Expression Profiling/methods
Humans
Neoplasm Proteins/analysis
Oligonucleotide Array Sequence Analysis/methods
Pattern Recognition, Automated/methods
Principal Component Analysis
Reproducibility of Results
Sensitivity and Specificity
Chemicals
Biomarkers, Tumor
Neoplasm Proteins
Authors & Affiliations
8 authors, click to expand affiliations / ORCID
Alexe Gabriela
Scanfeld Daniel
Tamayo Pablo
Mesirov Jill P
Ganesan Shridar
DeLisi Charles
Bhanot Gyan
References (19)
19 references, click to expand
-
Molecular profiling in breast cancer.
Rev Endocr Metab Disord. 2007 Sep;8(3):185-98
PMID: 17464566
-
Patterns and emerging mechanisms of the angiogenic switch during tumorigenesis.
Cell. 1996 Aug 9;86(3):353-64
PMID: 8756718
-
Microarray analysis and tumor classification.
N Engl J Med. 2006 Jun 8;354(23):2463-72
PMID: 16760446
-
Molecular portraits of human breast tumours.
Nature. 2000 Aug 17;406(6797):747-52
PMID: 10963602
-
Estrogen receptor status in breast cancer is associated with remarkably distinct gene expression patterns.
Cancer Res. 2001 Aug 15;61(16):5979-84
PMID: 11507038
-
The hallmarks of cancer.
Cell. 2000 Jan 7;100(1):57-70
PMID: 10647931
-
High expression of lymphocyte-associated genes in node-negative HER2+ breast cancers correlates with lower recurrence rates.
Cancer Res. 2007 Nov 15;67(22):10669-76
PMID: 18006808
-
Gene expression profiles of human breast cancer progression.
Proc Natl Acad Sci U S A. 2003 May 13;100(10):5974-9
PMID: 12714683
-
A gene network for navigating the literature.
Nat Genet. 2004 Jul;36(7):664
PMID: 15226743
-
Survivin, Survivin-2B, and Survivin-deItaEx3 expression in medulloblastoma: biologic markers of tumour morphology and clinical outcome.
Br J Cancer. 2005 Jan 31;92(2):359-65
PMID: 15655550
-
Survivin, a novel anti-apoptosis inhibitor, expression in uterine cervical cancer and relationship with prognostic factors.
Int J Gynecol Cancer. 2005 Jan-Feb;15(1):113-9
PMID: 15670305
-
DAVID: Database for Annotation, Visualization, and Integrated Discovery.
Genome Biol. 2003;4(5):P3
PMID: 12734009
-
Data perturbation independent diagnosis and validation of breast cancer subtypes using clustering and patterns.
Cancer Inform. 2007 Feb 19;2:243-74
PMID: 19458770
-
Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer.
Lancet. 2005 Feb 19-25;365(9460):671-9
PMID: 15721472
-
Gene expression profiles do not consistently predict the clinical treatment response in locally advanced breast cancer.
Mol Cancer Ther. 2006 Nov;5(11):2914-8
PMID: 17121939
-
MatchMiner: a tool for batch navigation among gene and gene product identifiers.
Genome Biol. 2003;4(4):R27
PMID: 12702208
-
Basal phenotype identifies a poor prognostic subgroup of breast cancer of clinical importance.
Eur J Cancer. 2006 Dec;42(18):3149-56
PMID: 17055256
-
Distinct molecular mechanisms underlying clinically relevant subtypes of breast cancer: gene expression analyses across three different platforms.
BMC Genomics. 2006;7:127
PMID: 16729877
-
Molecular classification of cancer: class discovery and class prediction by gene expression monitoring.
Science. 1999 Oct 15;286(5439):531-7
PMID: 10521349