Aging impacts transcriptomes but not genomes of hormone-dependent breast cancers
© Yau et al., licensee BioMed Central Ltd. 2007
Received: 20 July 2007
Accepted: 12 September 2007
Published: 12 September 2007
Age is one of the most important risk factors for human malignancies, including breast cancer; in addition, age at diagnosis has been shown to be an independent indicator of breast cancer prognosis. Except for inherited forms of breast cancer, however, there is little genetic or epigenetic understanding of the biological basis linking aging with sporadic breast cancer incidence and its clinical behavior.
DNA and RNA samples from matched estrogen receptor (ER)-positive sporadic breast cancers diagnosed in either younger (age ≤ 45 years) or older (age ≥ 70 years) Caucasian women were analyzed by array comparative genomic hybridization and by expression microarrays. Array comparative genomic hybridization data were analyzed using hierarchical clustering and supervised age cohort comparisons. Expression microarray data were analyzed using hierarchical clustering and gene set enrichment analysis; differential gene expression was also determined by conditional permutation, and an age signature was derived using prediction analysis of microarrays.
Hierarchical clustering of genome-wide copy-number changes in 71 ER-positive DNA samples (27 younger women, 44 older women) demonstrated two age-independent genotypes; one with few genomic changes other than 1q gain/16q loss, and another with amplifications and low-level gains/losses. Age cohort comparisons showed no significant differences in total or site-specific genomic breaks and amplicon frequencies. Hierarchical clustering of 5.1 K genes variably expressed in 101 ER-positive RNA samples (53 younger women, 48 older women) identified six transcriptome subtypes with an apparent age bias (P < 0.05). Samples with higher expression of a poor outcome-associated proliferation signature were predominantly (65%) younger cases. Supervised analysis identified cancer-associated genes differentially expressed between the cohorts; with younger cases expressing more cell cycle genes and more than threefold higher levels of the growth factor amphiregulin (AREG), and with older cases expressing higher levels of four different homeobox (HOX) genes in addition to ER (ESR1). An age signature validated against two other independent breast cancer datasets proved to have >80% accuracy in discerning younger from older ER-positive breast cancer cases with characteristic differences in AREG and ESR1 expression.
These findings suggest that epigenetic transcriptome changes, more than genotypic variation, account for age-associated differences in sporadic breast cancer incidence and prognosis.
Age is the strongest demographic risk factor for most human malignancies, including breast cancer . About 80% of all breast cancers occur in women older than age 50; the 10-year probability of developing invasive breast cancer increases from <1.5% at age 40, to about 3% at age 50 and to >4% by age 70, resulting in a cumulative lifetime risk of 13.2% (one in eight) and a near ninefold higher incidence rate in women older than age 50 as compared with their younger counterparts [2, 3]. Despite awareness that breast cancer and other cancers are primarily age-related diseases, molecular and cellular hypotheses explaining the cancer–aging relationship have only recently emerged and remain clinically unproven .
At the subcellular level, normal human aging has been linked to increased genomic instability [5, 6], to global and promoter-specific epigenetic changes [7, 8], and to altered expression of genes involved in cell division and extracellular matrix remodeling [5, 6]. These associations have led to the hypothesis that the cancer-prone phenotype of an older individual results from the combined effects of cumulative mutational load, increased epigenetic gene silencing, telomere dysfunction and altered stromal milieu . Given the worrisome social, economic and medical consequences of an aging worldwide population, proposed biological mechanisms linking cancer with aging must be established in order to develop effective interventions.
As with normal organs and tissues, tumor biology can also change with aging [10, 11]. For sporadic breast cancer in particular, correlations between patient age at diagnosis, tumor biology and clinical prognosis have long been appreciated if not fully understood [12–16]. Younger age at diagnosis (≤ 45 years old) is associated with more aggressive breast cancer biomarkers, including overexpression of ERBB2/HER2 and ERBB1/HER1 growth factor receptors , abnormal p53 expression [13, 15], estrogen receptor (ER) negativity [12–16], higher nuclear grade and higher Ki-67 proliferation index [12–14, 16]. These breast cancer biomarkers are also interdependent, however; in particular, ER expression is inversely correlated with abnormal p53 , overexpression of ERBB2 , high Ki-67 and nuclear grade, and poor patient prognosis . It therefore remains unclear whether the age-specific biomarker features of breast cancer reflect the pleotropic background effects of aging on the normal mammary gland or age-specific differences in breast tumorigenesis; also, since most age-specific biomarkers strongly associate with the ER status, the effects of aging must be studied in histologically similar breast cancer phenotypes controlled for ER status.
The molecular and cellular effects of aging on both normal and malignant breast tissue are superimposed on a continuum of developmental changes that normally occur between puberty and menopause, heavily influenced by menstrual history and parity. In general, the normal mammary gland ER content (fmol receptor/g tissue) as well as the proportion of ER-expressing (ER-positive) ductal epithelial cells increase with each decade of age, and reach a plateau with menopause at about age 50 [18, 19]. In contrast, breast cancer ER expression continues to rise beyond menopause, reaching a near 25-fold differential between normal and malignant mammary gland ER expression in patients by age 70 .
Curiously, expression of some ER-inducible gene markers, such as progesterone receptor (PR), pS2, Bcl2 and cathepsin D, does not show any significant relationship with the age at diagnosis [13, 18], while other markers show increased expression in breast cancers arising earlier in life  – suggesting that the effects of aging may in part be attributed to age-related differences in estrogen-inducible ER pathways. Important in this regard is the age-related change in PR coexpression within ER-positive breast cancers, since PR has long been used as a clinical indicator for a functioning ER pathway in tumors likely to respond to endocrine therapy . Among all ethnic patient groups, ER-positive/PR-negative breast cancers show the greatest age-related increase in incidence after age 40 . Potentially relevant to this ER-positive/PR-negative phenotype is the fact that growth-factor-activated pathways downregulate PR expression [22–25], and that the inverse correlation between overexpression of the ERBB2 growth factor receptor and PR positivity is only seen in breast cancers arising after age 40 . Surprisingly, the natural perimenopausal decline in ovarian-produced estrogen serum levels do not fully account for age-related changes in ER-regulated mammary epithelial pathways, since the marked age-related increase in stromal and epithelial aromatase expression produces postmenopausal mammary gland estrogen levels comparable with those measured in premenopausal women .
To better understand the molecular and cellular influences of aging on breast cancer biology and clinical behavior, we performed a detailed study of phenotypically similar breast cancers arising in two disparate patient age groups. The DNA and the RNA were prospectively extracted from cryobanked samples of stage-matched and histology-matched ER-positive breast cancers diagnosed in either younger (age ≤ 45 years) or older (age ≥ 70 years) Caucasian women. These samples were analyzed by array comparative genomic hybridization (CGH) and by high-throughput expression microarrays to look for genetic and epigenetic differences between the age cohorts. Unsupervised hierarchical clustering of the combined data from both cohorts was used to search for age biases in clustered subsets, and this was followed by supervised comparisons between the two cohorts to delineate potential age-related genomic and transcriptome differences. Finally, a predictive analysis of microarrays (PAM) performed on the two age cohorts produced an age-specific expression signature that proved to have >80% predictive accuracy when validated against two other independent breast cancer datasets.
Materials and methods
Breast cancer samples extracted for DNA and RNA
Cryobanked breast cancer specimens excised from newly diagnosed Caucasian females were obtained from the University of California San Francisco (UCSF) Comprehensive Cancer Center Breast Oncology Program Tissue Core (n = 66) and from the National Cancer Institute – Bari (NCI-Bari) (n = 71), following multi-institutional review board approvals. Tumor specimen selection criteria included sporadic incidence (no first-degree relatives with breast cancer), >50% invasive cancer cellularity, a frozen wet weight of at least 100 mg, ER-positivity (>10% nuclear immunohistochemical stain), and patient age at diagnosis either ≤ 45 years (younger cases) or ≥ 70 years (older cases). The 66 UCSF cases were all node-negative and were of predominantly ductal histology (60/66); 54 cases were associated with outcome annotation. The 71 NCI-Bari cases were all of ductal histology with mixed nodal status and without any outcome annotation. The two age cohorts within both tumor sets showed no significant imbalances in the stage, the grade, and the PR status or ERBB2 status. PR-positivity was defined as ≥ 10% nuclear immunohistochemical staining, and ERBB2-positivity was defined as gene amplification.
Cryobanked specimens were pulverized under liquid nitrogen prior to nucleic acid extraction. DNA was purified from frozen tumor powders using the High Pure PCR Preparation Kit (Roche Diagnostics, Indianapolis, IN, USA) and its quality was verified by gel electrophoresis. Total RNA was purified using Trizol reagent as per the manufacturer's protocol (Invitrogen, Carlsbad, CA, USA). Some RNA samples initially stored in formamide were further purified through RNeasy columns according to the manufacturer's protocol (QIAGEN, Valencia, CA, USA). All RNA samples were quality verified on a bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). DNA from all 71 NCI-Bari specimens were used for array-based CGH analysis, but only 35 of these specimens also yielded sufficient RNA for expression microarray analysis. The 35 NCI-Bari RNA samples (RNA sample set 1) were therefore combined with 66 RNA samples prepared from the UCSF specimens (RNA sample set 2) to yield 101 total RNA samples (53 younger women, 48 older women) for expression microarray analysis.
Array comparative genomic hybridization and data processing
Array CGH and data processing were carried out as described in reference . Test and reference genomic DNAs (500–1,000 ng) were labeled with Cy3 and Cy5, respectively, in a random priming reaction using 25–50 μl reaction volumes. The two labeled DNAs together with Cot-1 DNA (100 μg) were hybridized for 48–72 hours at 37°C onto arrays of 2,464 BAC clones, each printed in triplicate (HumArray1.14 and Hum-Array2.0; UCSF Comprehensive Cancer Center Microarray Core). Data from both array versions were combined only for the BAC clones present on both arrays; duplicate clones were averaged and the final dataset contained 2,240 unique BACs. Images were acquired and the data were processed as previously described . For each tumor, data were plotted in genome order as the mean log2 ratio of the replicate spots for each clone normalized to the genome median log2 ratio. Array CGH data were deposited in the public Gene Expression Omnibus database (GSE8801).
Array CGH data were analyzed using circular binary segmentation  with default parameters to translate intensity measurements into regions of equal copy number, as implemented in the DNAcopy R/Bioconductor package . Missing values were imputed using the maximum value of two flanking segments, producing smoothed values. The gain and loss status for each probe was assigned using the mergeLevel procedure . Tumor profiles were clustered using smoothed imputed data with outliers present; and agglomerative hierarchical clustering was performed using the Euclidean distance as a similarity measure and using the Ward method to minimize the sum of variances to produce compact spherical clusters. A clone-wise comparison of the phenotypic groups was made by t-testing smoothed values and controlling for the false discovery rate, with a false-discovery-rate-adjusted P value ≤ 0.05 for significance. The Kruskal–Wallis rank sum test was used to analyze phenotypic associations with the following autosomal genomic parameters: the number of break points and chromosomes with break points, the number of amplifications and chromosomes with amplifications, the number of whole chromosome changes and the fraction of genomes altered.
The fractions of genome gained and lost were computed for each tumor, and the frequency of alterations at each clone was computed as the proportion of samples showing an alteration at that locus. The amplification status for a clone was determined by considering the width of the segment to which that clone belonged (0, if an outlier) and a minimum difference between the smoothed value of the clone and the segment means of the neighboring segments. A clone was declared amplified if it belonged to the segment spanning less than 20 Mb and the minimum difference was greater than exp(-x3), where x is the final smoothed value for the clone. This procedure allowed clones with small log2-ratio values to be declared amplified if they were high relative to surrounding clones. To calculate the number of chromosomes with amplifications, a chromosome was said to be amplified if at least one of its clones was amplified. Whole chromosome changes were assigned to chromosomes without identified breakpoints if a chromosomal segment mapped to a gain or a loss level. The number of chromosomal break points was calculated as the number of copy number levels within each chromosome across the genome minus the number of chromosomes. To calculate the number of chromosomes with break points, at least one break point per chromosome was necessary.
Microarray expression profiles and data processing
The total RNA (3–5 μg per sample) was labeled and analyzed using Affymetrix (Santa Clara, CA, USA) HT-HG_U133A Early Access Arrays with 22.9 K probes representing ~13 K unique UniGenes. Analyses were performed by standard Affymetrix procedures within the Lawrence Berkeley National Laboratory and Life Science Divison's Molecular Profiling Laboratory . Probe set measurements were generated from quantified Affymetrix image files (.CEL files) using the RMA algorithm in Bioconductor R. Array data were deposited in the public Gene Expression Omnibus database (GSE7378 and GSE8193).
Gene expression values were mean centered, with a low variation filter applied to exclude probe sets that did not have at least 10 observations exhibiting a twofold change from the mean. Filtered probes were annotated (GeneTraffic annotation file, March 2006) and those with unknown UniGene symbols were omitted, yielding a final significant probe set of 6,632 annotated probes representing 5,109 unique genes. Unsupervised hierarchical clustering of the mean centered significant probe set was performed using Cluster , and was visualized with Java TreeView . Phenotypes (for example, age cohort, PR status and ERBB2 status) of the resulting clusters were compared by chi-square test.
Gene set enrichment analysis (GSEA) was performed using GSEA software (version 2.01) [35, 36] to assess tumor phenotypes (for example, age cohort, PR status or ERBB2 status) with respect to specific gene signatures including a mitogen-activated protein kinase (MAPK)-regulated gene set , a luminal tumor subtype gene set  and a proliferation-associated gene set . The enrichment of all curated gene sets (c2) within the Molecular Signature Database , satisfying the gene set size-filtering criteria (maximum = 500, minimum = 10), were also evaluated. Ranking of UniGenes within each phenotype was based on a signal-to-noise metric; an enrichment score for each signature was derived as a function of the likelihood of that gene set being among the most highly ranked genes within the phenotype. For genes represented by more than one probe, the median expression level of all corresponding gene probes was used; the significance of the enrichment score was estimated by 1,000 random permutations of the sample labels in the tumor dataset. For gene sets showing statistically significant enrichment with respect to the age cohort (P < 0.05), unsupervised clustering was performed on the entire set of RNA samples using those genes; the outcome association of high gene set expressors versus low gene set expressors was also tested by Kaplan–Meier analysis of the 54 cases with known recurrence events.
A conditional permutation test was performed to identify differentially expressed genes with respect to age and conditional on the RNA sample sources (NCI-Bari and UCSF). For each gene, a linear-fit model of expression level versus age cohort was adjusted for the categorical sample source, and the P value was based on the t statistic on the slope of age. The null distribution used, however, was not the standard normal but was based on permuting the age cohort distribution within each sample source, which provides an exact test of the conditional independence of age given the sample source. This was expected to control for possible different age effects within each sample source, and to adjust for potential spurious associations based on the different sample sources . Reported associations between the differentially expressed genes and either cancer or aging were determined using the MEDGENE database [41, 42], entering the following search terms: Neoplasm, Breast Neoplasms, Carcinomas, Breast Carcinomas, and Aging, Premature. Functional annotation of the differentially expressed genes was performed by gene ontology (GO) analysis ; significant enrichment for specific biological functional categories (≥ 5 probes within a process class, Expression Analysis Systematic Explorer  score <0.05) was identified from the DAVID database [43, 45].
An age cohort gene signature was obtained by training PAM software  on RNA sample set 2 (UCSF), minimizing the cross-validation error for the individual age cohorts. The resulting PAM classifier was used to predict the age cohort of RNA sample set 1 (NCI-Bari). For external validation purposes, expression microarray data from Miller and colleagues  and from Sotiriou and colleagues  were downloaded from the NCBI public Gene Expression Omnibus database (GSE3494 and GSE2990). ER-positive cases with age characteristics matching our age cohorts were selected, resulting in 102 cases from Miller and colleagues' study and 47 cases from Sotiriou and colleagues' study. As the training set (RNA sample set 2) comprised only early-stage tumors, PAM validation was also performed on the external datasets restricted to only ER-positive, node-negative cases (35 cases from Sotiriou and colleagues , and 64 cases from Miller and colleagues ). Since these microarray studies used Affymetrix U133A platforms, probe set measurements were generated using the RMA algorithm in Bioconductor R and the resulting data were mapped to our significant gene set by Affymetrix probe identifiers, to which the PAM classifier was applied. The significance of the prediction accuracy was determined using Fisher's exact test.
Age and ER-positive breast cancer genomic profiles
Unsupervised analysis of ER-positive breast cancer expression profiles
When the cluster compositions with respect to age, PR status and ERBB2 status were statistically compared, only the age cohort distribution was found to be significantly different (P < 0.05, chi-square test). Patient outcome data were available on 54 (30 younger, 24 older) of the 101 cases, scattered evenly among the six transcriptome clusters. Kaplan–Meier probability curves for recurrence-free survival indicated that younger age and PR-negative status were associated with earlier relapse, but these outcome differences did not quite reach statistical significance (age, P = 0.09; PR status, P = 0.08; log-rank analyses). Kaplan–Meier curves for cases representing each of the six transcriptome clusters, however, did achieve significant separation (P = 0.025, log-rank analysis), with the predominantly younger group 2B cases showing the shortest survival (median recurrence-free survival = 2.5 years) and the predominantly older group 2A cases showing significantly more prolonged survival (median recurrence-free survival = 6.2 years).
Gene set enrichment analysis
Gene set enrichment analysis results for the gene list tested
Old over young cohort
Progesterone receptor-negative over progesterone receptor-positive tumors
ERBB2-negative over ERBB2-positive tumors
Nominal P value
FWER P value
Nominal P value
FWER P value
Nominal P value
FWER P value
Differential gene expression between age cohorts
Differentially expressed genes between the young and old cohorts
Genes with higher expression in the young cohort
Amphiregulin (schwannoma-derived growth factor)
Protease, serine, 2 (trypsin 2)
Parathyroid hormone-like hormone
Hydroxyprostaglandin dehydrogenase 15-(NAD)
Sperm protein associated with the nucleus, X-linked, family member A1///SPANX family, member B1///SPANX family, member A2///SPANX family, member C///SPANX family, member B2
Laminin, alpha 3
ATPase, H+ transporting, lysosomal 56/58 kDa, V1 subunit B, isoform 1 (renal tubular acidosis with deafness)
S100 calcium binding protein A2
Deiodinase, iodothyronine, type II
Protease, serine, 1 (trypsin 1)///protease, serine, 2 (trypsin 2)///protease, serine, 3 (mesotrypsin)///trypsinogen C
Fibroblast growth factor receptor 1
Tumor protein p73-like
Protease, serine, 1 (trypsin 1)
Chromosome 20 open reading frame 59
Discs, large homolog 7 (Drosophila)
Elongation of very long chain fatty acids (FEN1/Elo2, SUR4/Elo3, yeast)-like 2
Kinesin family member 2C
Serine/threonine kinase 6
CDC14 cell division cycle 14 homolog A (S. cerevisiae)
Elongation factor RNA polymerase II-like 3
RAD54 homolog B (S. cerevisiae)
Integrin, alpha 2 (CD49B, α2 subunit of VLA-2 receptor)
Genes with higher expression in the old cohort
Transmembrane channel-like 5
ST6 (α-N-acetyl-neuraminyl-2,3-β-galactosyl-1,3)-N-acetylgalactosaminide α-2,6-sialyltransferase 5///ST6 (α-N-acetyl-neuraminyl-2,3-β-galactosyl-1,3)-N-acetylgalactosaminide α-2,6-sialyltransferase 5
Phosphorylase, glycogen; liver (Hers disease, glycogen storage disease type VI)
Tumor necrosis factor (ligand) superfamily, member 10///tumor necrosis factor (ligand) superfamily, member 10
Golgi phosphoprotein 2
Dermatan sulfate proteoglycan 3
Hypothetical protein FLJ20152
Glycine amidinotransferase (l-arginine:glycine amidinotransferase)
Ectonucleoside triphosphate diphosphohydrolase 5
SAM and SH3 domain containing 1
Inositol 1,4,5-triphosphate receptor, type 1
Angiogenin, ribonuclease, RNase A family, 5///ribonuclease, RNase A family, 4
IQ motif containing GTPase activating protein 2
MANSC domain containing 1
Fumarylacetoacetate hydrolase (fumarylacetoacetase)
Rho GDP dissociation inhibitor (GDI) beta
TAP binding protein-like
Calmin (calponin-like, transmembrane)
Estrogen receptor 1
p8 protein (candidate of metastasis 1)
Sterol-C5-desaturase (ERG3 δ-5-desaturase homolog, fungal)-like
C-type lectin domain family 5, member A
Ras homolog gene family, member B
Phosphodiesterase 4A, cAMP-specific (phosphodiesterase E2 dunce homolog, Drosophila)
Chromosome 21 open reading frame 25
Chemokine (C–C motif) ligand 3///chemokine (C–C motif) ligand 3-like 1///chemokine (C–C motif) ligand 3-like 3
Coiled-coil domain containing 28A
Calmodulin 3 (phosphorylase kinase, delta)
PTPRF interacting protein, binding protein 2 (liprin β2)
Diazepam binding inhibitor (GABA receptor modulator, acyl-Coenzyme A binding protein)
Solute carrier family 25 (mitochondrial carrier, Aralar), member 12
Membrane-associated ring finger (C3HC4) 8
Solute carrier family 12 (potassium/chloride transporters), member 8
Fucosidase, alpha-L-1, tissue
RAN binding protein 2
trans-Golgi network protein 2
von Willebrand factor A domain containing 1
A MEDGENE database search (disease terms 'neoplasms', 'breast neoplasms', 'carcinomas', or 'carcinoma, ductal, breast') indicated that 29 of the 75 differentially expressed genes had a published first-degree association with cancer; in contrast, none of the 75 genes had any published association with aging (disease terms 'aging, premature'). GO analysis was used to determine whether the 75 genes differentially expressed between the age cohorts were enriched for specific biological processes. Several functional categories were significantly over-represented (Expression Analysis Systematic Explorer score < 0.05) in this gene set , including development, cell cycle, M-phase, morphogenesis and reproduction. As noted in Table 2, 17 of the genes were associated with development and nine genes were associated with the cell cycle, including five genes specifically associated with the M-phase.
Defining and validating an age cohort signature in ER-positive breast cancers
ESR1 and AREG were among the 44 genes in common between the age signature gene set and the differentially expressed gene set; and, as shown in Figure 4c, age-signature-defined subsets from all four sample datasets showed similar differences in mean expression levels (log2 intensity) for both ESR1 and AREG. Of note, a PR signature was similarly derived from RNA sample set 2 (103 probes; Additional file 4), but showed only 67% accuracy (P = 0.1) in predicting the PR status of RNA sample set 1.
Although there have been numerous studies of clinical factors addressing the relationship between age at diagnosis and breast cancer prognosis [12, 14, 16, 53–55], few studies have comprehensively investigated the age dependency of the many well-established prognostic breast cancer biomarkers, and no studies have used a prospective study design [13, 18]. Concerned about the established inverse relationship between the ER status and poor-risk biomarker surrogates of breast cancer proliferation and genomic instability [13, 18], the present study aimed to identify genomic and transcriptome changes associated with aging using DNA and RNA prospectively collected from stage-matched and histology-matched ER-positive breast cancers from younger women (age ≤ 45 years) and older women (age ≥ 70 years), analyzed by array CGH and high-throughput expression microarrays.
Similar bioinformatics-based approaches have been used to characterize aging effects in human fibroblasts [5, 6], lymphocytes  and myoblasts ; however, comparable efforts to investigate aging influences on human cancer biology have not been reported. Moreover, while ER-positive breast cancers have been well studied as a subgroup within unselected breast cancer phenotypes using array CGH [28, 57] or expression profiling [38, 39, 49–51], the present study represents the largest study reported to date using these powerful techniques to subset ER-positive breast cancers, while employing a statistical design powered to detect age-specific differences.
Array CGH analysis of 71 DNA samples confirmed that our ER-positive breast cancers were composed of two basic genotypes : a simple subtype characterized by few genomic copy number changes other than gain of 1q and loss of 16q, and a mixed amplifier subtype characterized by recurrent amplifications but otherwise low levels of genomic gains and losses. A third genomic subtype of breast cancer, referred to as complex, known to be almost exclusively composed of ER-negative breast cancers , was not observed in either of the two age cohorts studied. Neither the simple nor the mixed amplifier genomic subtypes of ER-positive breast cancer showed any particular age bias. Direct comparison of the two age cohorts for multiple array CGH parameters also revealed no significant differences in the fraction of genome altered, in whole chromosome changes or in total or site-specific amplicon frequencies. Although nonsignificant trends suggested slightly fewer oncogene amplifications within the older cohort, overall amplification frequencies for the most common oncogenes were as expected for ER-positive breast cancers [51, 58]: MYC (27%), CCND1 (23%), ZNF217 (17%), AIB1 (16%), MDM2 (8%), ESR1 (7%), ERBB2 (7%), and TOPO2A (7%). At the level of genomic resolution (~1 MB) achievable by BAC-based array CGH, there appeared to be few if any genetic differences between ER-positive breast cancers arising in women whose ages differ by more than 25 years. Future studies employing higher density genomic arrays are warranted to confirm this conclusion.
Microarray profiling of 101 RNA samples showed an average 65-fold range in ESR1 transcript levels across the entire collection of ER-positive breast cancers, with the older cohort showing significantly higher ESR1 levels as compared with the younger cohort, consistent with earlier biomarker studies . There was the expected close correlation between the ESR1 transcript levels and commonly observed ESR1 coexpressed genes (for example, GATA3) as well as other genes (for example, KRT8, KRT18) that characteristically define luminal-type breast cancer, although this tumor collection also contained several ERBB2-positive cases (10/101) that are not characteristically found in microarray-defined clusters of luminal-type breast cancer [38, 39, 49–51]. Hierarchical clustering of the ~5.1 K variably expressed genes also identified six transcriptome subtypes of ER-positive breast cancer with significant age biases (P < 0.05) but not associated with differing PR status. Based on relapse-free survival analyses of the 54 cases with known clinical outcome (30 younger women, 24 older women), there was a trend supporting a less favorable prognosis for the younger age cases (P = 0.09) and PR-negative cases (P = 0.08). The six age-biased transcriptome clusters, however, showed significantly different relapse-free survival outcomes (P = 0.025, log-rank analysis), suggesting that these transcriptome subtypes represent clinically relevant phenotypes of ER-positive breast cancer. Previous expression array studies analyzing fewer ER-positive cases have identified no more than two or three subsets of luminal-type breast cancer [38, 39, 49–51].
Reported gene signatures representing luminal, proliferation and MAPK markers were tested for their enrichment in one or the other of the age-stratified cohorts, and only the proliferation gene signature showed any significant age bias when multiple testing was accounted for, being more highly expressed in the younger cohort. This finding is consistent with earlier studies showing higher tumor grade and proliferation markers (for example, mitotic index and Ki-67 positivity) in younger age breast cancer patients . While none of the >1,000 curated gene sets in the Molecular Signature Database that were similarly evaluated demonstrated any significant age biases when multiple testing was account for, a trend was observed for enrichment of cell cycle genes in the younger cohort cases. Nine genes common to both the GO biological process cell cycle set and the proliferation signature set (BUB1, CCNB1, CCNE2, CDC25A, CDC7, MAD2L1, MCM4, ORC6L, PTTG1) were also present in our significant probe set. Among these, four genes (BUB1, CCNE2, MAD2L1, ORC6L) have been previously associated with poor-prognosis ER-positive breast cancers in a well-established 70-gene prognostic signature ; these genes are therefore probably important contributors to the more aggressive tumor characteristics of ER-positive breast cancers arising in younger patients.
Using only the proliferation gene signature to perform unsupervised hierarchical clustering of the 101 cases generated two comparably sized ER-positive subsets, one with higher expression and another with lower expression of the proliferation genes; the higher expressing subset contained most of the younger age cases (34/52) and all but one of the ERBB2-positive cases. When this proliferation signature was also used to dichotomize the 54 cases with known clinical outcome, the higher expressing cases showed significantly worse disease-free survival as compared with the lower expressing cases, consistent with reports on the association of a similar proliferation signature with poor outcome in patients with ER-positive breast cancer . Interestingly, despite a presumed mechanistic link between activation of growth factor receptors, MAPK signaling and cell proliferation, there was minimal overlap between genes in the reported MAPK and proliferation signatures, and no significant association was observed between the MAPK signature, age and ERBB2 positivity.
Despite the observed positive association between the ESR1 expression level and older age, no age association was seen for the luminal gene signature that included ESR1, ESR1-associated genes and estrogen-inducible genes. This finding is consistent with our previous report showing increased breast cancer ER protein with aging without comparably increased levels of such estrogen-inducible markers as PR, pS2, Bcl2 and cathepsin D , and suggesting reduced estrogen signaling in breast tumors of older patients. In keeping with these protein biomarker observations, differential gene expression analysis in the present study did not identify any known estrogen-inducible genes such as TFF1, PGR, IRS1, IGFBP4, PCNA, MYC, CCNA2 or DLEU2 as being more highly expressed in the older cohort despite higher expression of ESR1 in this cohort. In contrast, two estrogen-inducible growth-regulating genes, GREB1 and AREG, showed significantly higher expression levels in the younger cohort, in keeping with a recent study demonstrating a negative correlation between these estrogen-inducible genes and age . As GREB1 and AREG are known to induce cell proliferation upon estrogen activation [60, 61], their increased expression in the younger cohort offers some mechanistic basis for increased proliferative activity and gene expression in the younger cohort.
Of the 75 unique genes differentially expressed between younger and older cohorts, 24 genes showed increased expression in younger cases relative to older cases (including GREB1 and AREG) while 51 genes showed increased expression in older cases relative to younger cases (including ESR1). Comparison with a well-studied estrogen-inducible gene signature set  revealed that ~25% (19/75) of these differentially expressed genes overlapped with known early or late estrogen-responsive genes, and thus potentially reflected hormonal changes associated with menopause rather than aging effects. While two-thirds (13/19) of these potential estrogen-responsive genes showed appropriate directional changes according to cohort menopausal status, supporting this possibility, at least 75% of the differentially expressed genes would appear to be independent of menopausal differences in circulating estrogen levels and, therefore, potentially informative of age-related differences in ER-positive breast cancer biology. A comprehensive database search confirmed that at least 40% of these differentially expressed genes have reported direct links with malignancy; and while none have reported links with premature aging, one of the differentially expressed genes (KIF2C) has been previously implicated in aging studies of lymphocytes and fibroblasts , while six other genes (COBLL1, HPGD, HOXB2, PDE4A, SLC25A12, TP73L) were recently reported as differentially expressed with age in human skeletal muscle .
A search for annotated enrichment of the differentially expressed genes for specific biological processes (GO Biological Processes, Expression Analysis Systematic Explorer score < 0.05) indicated that 'development' and 'cell cycle/M-phase' were the most overrepresented functional gene categories. In keeping with the GSEA observation indicating a trend for enrichment of cell-cycle-associated genes in the younger cohort cases, differentially expressed cell cycle/M-phase genes (including positive regulators such as STK6, FGFR1 and DLG7) represented 20% (5/25) of all genes overexpressed in the younger cohort but only 8% (4/51) of those overexpressed in the older cohort. In contrast, the older cohort cases showed differentially increased expression of negative cell cycle regulators (such as SASHI and RHOB) and four developmentally essential homeobox genes (HOXB2, HOXB5, HOXB6, HOXB7), the latter finding also in keeping with the GSEA observed trend showing enrichment in the older cohort of HOX-regulated (NUP90-HOXA9 repressed) genes. Two of the overexpressed HOXB genes (HOXB6, HOXB7) have been specifically linked to mammary gland development and are known to be expressed in ER-positive breast cancer cells . HOXB7, in particular, known to be dependent on stromal (extracellular matrix) signaling, is transcriptionally upregulated in breast cancers metastatic to bone (relative to primary tumors), and is thought to play a role in promoting angiogenesis, growth factor-independent proliferation and DNA double-strand break repair, conferring breast cancer resistance to the genome destabilizing effects of DNA damage .
PAM was used to derive an age signature that consisted of 128 unique genes, including 44 of the 75 differentially expressed genes determined by our conditional permutation approach. The age signature was independently validated against two other age-matched ER-positive breast cancer microarray datasets and proved to have >80% accuracy in distinguishing younger from older ER-positive breast cancer cases. ESR1 and AREG were among the genes in common between the age signature and the differentially expressed gene sets; it is therefore not surprising that the age-signature-defined subsets from the two independent databases showed similar differences in the mean expression levels of these two genes as found in our age-defined cohorts. Only 28% of the age signature genes overlap with known early or late estrogen-responsive genes, suggesting that this age signature largely reflects age-related differences in the phenotype of ER-positive breast cancer rather than differences in circulating estrogen levels associated with menopausal status.
The fact that a PAM-derived PR signature did not perform well upon validation implies substantial heterogeneity between ER-positive breast cancers with the same PR status, and possibly indicates that confounding age-related gene expression changes are of greater biological importance than PR-related gene expression differences. Misclassification errors using the age signature were more prevalent among the older cohort cases, also suggesting greater variation in expression of the age signature genes with aging. Of further interest, the 128 age signature gene set was unable to accurately subset ER-negative cases identified from the two independent breast cancer datasets [47, 48], consistent with expression-array-based conclusions that the biology of ER-positive and ER-negative breast cancers are fundamentally distinct, and supporting the likelihood that the PAM-derived age signature incorporates biological profiles specific to ER-positive breast cancers but not ER-negative breast cancers.
This prospectively designed study addresses a pressing need to evaluate molecular and cellular hypotheses proposed to explain age-related differences in breast cancer incidence and clinical behavior. It is hard to reconcile the evidence gathered in this study of ER-positive breast cancers with the more general cancer-aging postulate that the breast-cancer-prone phenotype of an older woman results from genomic instability and age-accumulated mutational loads secondary to telomeric dysfunction and/or progressive DNA damage . More consistent with the present evidence is the likelihood that ER-positive breast cancers arising in older women relative to younger women do so by a fundamentally different tumorigenic process, manifested more by epigenetic transcriptome differences such as those regulated by HOX genes, and less by genomic differences that were not detected using state-of-the-art BAC-based CGH analyses. More pronounced expression of cell cycle and proliferation-associated genes emerged as a strong defining feature of ER-positive breast cancers arising in younger women, perhaps even driving their earlier clinical appearance; this observation is certainly consistent with the more aggressive clinical nature of early-age-onset breast cancer.
Age cohort study designs of this type are needed to not only confirm the specific transcriptome differences noted here, but also to look for common age-associated differences in gene classes and functional pathways that may enable us to generalize about the age-related biological differences driving ER-negative breast tumorigenesis as well as the many other age-associated epithelial malignancies other than breast cancer.
comparative genomic hybridization
gene set enrichment analysis
mitogen-activated protein kinase
National Cancer Institute – Bari
prediction analysis of microarrays
University of California San Francisco.
The authors appreciate the considerable assistance of Heidi Feiler, PhD, Director of the Affymetrix HTArray Genomics Core in the Molecular Profiling Laboratory of the Lawrence Berkeley National Laboratory. CY participated in these studies in partial fulfillment of dissertation research within the Joint University of California Berkeley – University of California San Francisco Bioengineering Graduate Program. These studies were supported in part by National Institutes of Health grants R01-AG020521, R01-CA71468, P30-CA82103, and P50-CA58207, and Hazel P. Munroe memorial funding to the Buck Institute.
- Edwards B, Howe HL, Ries L, Thun M, Rosenberg H, Wingo P, Jemal A, Feigal E: Annual Report to the Nation on the Status of Cancer, 1973–1999 featuring implications of age and aging on US cancer burden. Cancer. 2000, 94: 2766-2792. 10.1002/cncr.10593.View ArticleGoogle Scholar
- Smigal C, Jemal A, Ward E, Cokkinides V, Smith R, Howe HL, Thun M: Trends in breast cancer by race and ethnicity: update 2006. CA Cancer J Clin. 2006, 56: 168-183.View ArticlePubMedGoogle Scholar
- DevCan – probability of developing or dying of cancer software. [http://www.srab.cancer.gov/devcan]
- Benz CC, Campisi J, Cohen HJ, Ershler WB, Haubein L, Irminger-Finger I: Meeting report: translational research at the aging and cancer interface. Cancer Res. 2007, 67: 4560-4563. 10.1158/0008-5472.CAN-07-1236.View ArticlePubMedGoogle Scholar
- Geigl JB, Langer S, Barwisch S, Pfleghaar K, Lederer G, Speicher MR: Analysis of gene expression patterns and chromosomal changes associated with aging. Cancer Res. 2004, 64: 8550-8557. 10.1158/0008-5472.CAN-04-2151.View ArticlePubMedGoogle Scholar
- Ly DH, Lockhart DJ, Lerner RA, Schultz PG: Mitotic misregulation and human aging. Science. 2000, 287: 2486-2492. 10.1126/science.287.5462.2486.View ArticlePubMedGoogle Scholar
- Issa J-P: Aging, DNA methylation and cancer. Crit Rev Oncol Hematol. 1999, 32: 31-43. 10.1016/S1040-8428(99)00019-0.View ArticlePubMedGoogle Scholar
- Richardson B: Impact of aging on DNA methylation. Ageing Res Rev. 2003, 2: 245-261. 10.1016/S1568-1637(03)00010-2.View ArticlePubMedGoogle Scholar
- DePinho RA: The age of cancer. Nature. 2000, 408: 248-254. 10.1038/35041694.View ArticlePubMedGoogle Scholar
- Ershler WB, Longo DL: Aging and cancer: issues of basic and clinical science. J Natl Cancer Inst. 1997, 89: 1489-1497. 10.1093/jnci/89.20.1489.View ArticlePubMedGoogle Scholar
- Balducci L, Ershler WB: Cancer and ageing: a nexus at several levels. Nature Rev Cancer. 2005, 5: 655-662. 10.1038/nrc1675.View ArticleGoogle Scholar
- El Saghir N, Seoud M, Khalil M, Charafeddine M, Salem Z, Geara F, Shamseddine A: Effects of young age at presentation on survival in breast cancer. BMC Cancer. 2006, 6: 194-201. 10.1186/1471-2407-6-194.View ArticlePubMedPubMed CentralGoogle Scholar
- Eppenberger-Castori S, Moore IIDH, Thor AD, Edgerton SM, Kueng W, Eppenberger U, Benz CC: Age-associated biomarker profiles of human breast cancer. Int J Biochem Cell Biol. 2002, 34: 1318-1330. 10.1016/S1357-2725(02)00052-3.View ArticlePubMedGoogle Scholar
- Figueiredo JC, Ennis M, Knight JA, McLaughlin JR, Hood N, O'Malley F, Andrulis IL, Goodwin PJ: Influence of young age at diagnosis and family history of breast or ovarian cancer on breast cancer outcomes in a population-based cohort study. Breast Cancer Res Treat. 2007, 105: 69-80. 10.1007/s10549-006-9433-3.View ArticlePubMedGoogle Scholar
- Rudolph P, Olsson H, Bonatz G, Ratjen V, Bolte H, Baldetorp B, Fernö M, Parwaresch R, Alm P: Correlation between p53, c-erb B-2, and topoisomerase IIa expression, DNA ploidy, hormonal receptor status and proliferation in 356 node-negative breast carcinomas: prognostic implications. J Pathol. 1999, 187: 207-216. 10.1002/(SICI)1096-9896(199901)187:2<207::AID-PATH223>3.0.CO;2-U.View ArticlePubMedGoogle Scholar
- Rapiti E, Fioretta G, Verkooijen HM, Vlastos G, Schafer P, Sappino A-P, Kurtz J, Neyroud-Caspar I, Bouchardy C: Survival of young and older breast cancer patients in Geneva from 1990 to 2001. Eur J Cancer. 2005, 41: 1446-1452. 10.1016/j.ejca.2005.02.029.View ArticlePubMedGoogle Scholar
- Harvey JM, Clark GM, Osborne CK, Allred DC: Estrogen receptor status by immunohistochemistry is superior to the ligand-binding assay for predicting response to adjuvant endocrine therapy in breast cancer. J Clin Oncol. 1999, 17: 1474-1481.PubMedGoogle Scholar
- Quong J, Eppenberger-Castori S, Moore D, Scott GK, Birrer MJ, Kueng W, Eppenberger U, Benz CC: Age-dependent changes in breast cancer hormone receptors and oxidant stress markers. Breast Cancer Res Treat. 2002, 76: 221-236. 10.1023/A:1020886801674.View ArticlePubMedGoogle Scholar
- Shoker BS, Jarvis C, Sibson DR, Walker C, Sloane JP: Oestrogen receptor expression in the normal and pre-cancerous breast. J Pathol. 1999, 188: 237-244. 10.1002/(SICI)1096-9896(199907)188:3<237::AID-PATH343>3.0.CO;2-8.View ArticlePubMedGoogle Scholar
- Creighton C, Cordero K, Larios J, Miller R, Johnson M, Chinnaiyan A, Lippman M, Rae J: Genes regulated by estrogen in breast tumor cells in vitro are similarly regulated in vivo in tumor xenografts and human breast tumors. Genome Biol. 2006, 7: R28-10.1186/gb-2006-7-4-r28.View ArticlePubMedPubMed CentralGoogle Scholar
- Horwitz KB, McGuire WL: Estrogen control of progesterone receptor in human breast cancer. J Biol Chem. 1978, 253: 2223-2228.PubMedGoogle Scholar
- Arpino G, Weiss H, Lee AV, Schiff R, De Placido S, Osborne CK, Elledge RM: Estrogen receptor-positive, progesterone receptor-negative breast cancer: association with growth factor receptor expression and tamoxifen resistance. J Natl Cancer Inst. 2005, 97: 1254-1261.View ArticlePubMedGoogle Scholar
- Cui X, Zhang P, Deng W, Oesterreich S, Lu Y, Mills GB, Lee AV: Insulin-like growth factor-I inhibits progesterone receptor expression in breast cancer cells via the phosphatidylinositol 3-kinase/Akt/mammalian target of rapamycin pathway: progesterone receptor as a potential indicator of growth factor activity in breast cancer. Mol Endocrinol. 2003, 17: 575-588. 10.1210/me.2002-0318.View ArticlePubMedGoogle Scholar
- Cui X, Schiff R, Arpino G, Osborne CK, Lee AV: Biology of progesterone receptor loss in breast cancer and its implications for endocrine therapy. J Clin Oncol. 2005, 23: 7721-7735. 10.1200/JCO.2005.09.004.View ArticlePubMedGoogle Scholar
- Kim H-J, Cui X, Hilsenbeck SG, Lee AV: Progesterone receptor loss correlates with human epidermal growth factor receptor 2 overexpression in estrogen receptor-positive breast cancer. Clin Cancer Res. 2006, 12: 1013s-1018s. 10.1158/1078-0432.CCR-05-2128.View ArticlePubMedGoogle Scholar
- Huei-Jean H, Patrick N, Maria D, Robert P, Hans W, Erik Van L, Patrick B, Frederic A, Marie Rose C, Ignace V: Association between HER-2/neu and the progesterone receptor in oestrogen-dependent breast cancer is age-related. Breast Cancer Res Treat. 2005, 91: 81-87. 10.1007/s10549-004-8235-8.View ArticleGoogle Scholar
- van Landeghem AA, Poortman J, Nabuurs M, Thijssen JH: Endogenous concentration and subcellular distribution of estrogens in normal and malignant human breast tissue. Cancer Res. 1985, 45: 2900-2906.PubMedGoogle Scholar
- Fridlyand J, Snijders A, Ylstra B, Li H, Olshen A, Segraves R, Dairkee S, Tokuyasu T, Ljung B, Jain A, et al: Breast tumor copy number aberration phenotypes and genomic instability. BMC Cancer. 2006, 6: 96-108. 10.1186/1471-2407-6-96.View ArticlePubMedPubMed CentralGoogle Scholar
- Olshen AB, Venkatraman ES, Lucito R, Wigler M: Circular binary segmentation for the analysis of array-based DNA copy number data. Biostatistics. 2004, 5: 557-572. 10.1093/biostatistics/kxh008.View ArticlePubMedGoogle Scholar
- Gentleman R, Carey V, Bates D, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.View ArticlePubMedPubMed CentralGoogle Scholar
- Willenbrock H, Fridlyand J: A comparison study: applying segmentation to array CGH data for downstream analyses. Bioinformatics. 2005, 21: 4084-4091. 10.1093/bioinformatics/bti677.View ArticlePubMedGoogle Scholar
- Lawrence Berkeley National Laboratory and Life Science Divison's Molecular Profiling Laboratory. [http://hta.lbl.gov/HTA/index.html]
- Eisen MB, Spellman PT, Brown PO, 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Saldanha AJ: Java Treeview – extensible visualization of microarray data. Bioinformatics. 2004, 20: 3246-3248. 10.1093/bioinformatics/bth349.View ArticlePubMedGoogle Scholar
- GSEA software. [http://www.broad.mit.edu/gsea/]
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al: From the cover: gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102: 15545-15550. 10.1073/pnas.0506580102.View ArticlePubMedPubMed CentralGoogle Scholar
- Creighton CJ, Hilger AM, Murthy S, Rae JM, Chinnaiyan AM, El-Ashry D: Activation of mitogen-activated protein kinase in estrogen receptor α-positive breast cancer cells in vitro induces an in vivo molecular phenotype of estrogen receptor α-negative human breast tumors. Cancer Res. 2006, 66: 3903-3911. 10.1158/0008-5472.CAN-05-4363.View ArticlePubMedGoogle Scholar
- Hu Z, Fan C, Oh D, Marron JS, He X, Qaqish B, Livasy C, Carey L, Reynolds E, Dressler L, et al: The molecular portraits of breast tumors are conserved across microarray platforms. BMC Genom. 2006, 7: 96-107. 10.1186/1471-2164-7-96.View ArticleGoogle Scholar
- Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, et al: Molecular portraits of human breast tumours. Nature. 2000, 406: 747-752. 10.1038/35021093.View ArticlePubMedGoogle Scholar
- Lehman EL: Testing Statistical Hypotheses. 1986, New York: Wiley, 2View ArticleGoogle Scholar
- MEDGENE database. [http://medgene.med.harvard.edu/MEDGENE/]
- Hu Y, Hines LM, Weng H, Zuo D, Rivera M, Richardson A, LaBaer J: Analysis of genomic and proteomic data using advanced literature mining. J Proteome Res. 2003, 2: 405-412. 10.1021/pr0340227.View ArticlePubMedGoogle Scholar
- Dennis G, Sherman B, Hosack D, Yang J, Gao W, Lane H, Lempicki R: DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003, 4: R60-10.1186/gb-2003-4-9-r60.View ArticlePubMed CentralGoogle Scholar
- Hosack D, Dennis G, Sherman B, Lane H, Lempicki R: Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4: P4-10.1186/gb-2003-4-6-p4.View ArticleGoogle Scholar
- DAVID database. [http://david.abcc.ncifcrf.gov/home.jsp]
- Predictive Analysis of Microarrays software. [http://www-stat.stanford.edu/~tibs/PAM/]
- Miller LD, Smeds J, George J, Vega VB, Vergara L, Ploner A, Pawitan Y, Hall P, Klaar S, Liu ET, et al: From the cover: an expression signature for p53 status in human breast cancer predicts mutation status, transcriptional effects, and patient survival. Proc Natl Acad Sci USA. 2005, 102: 13550-13555. 10.1073/pnas.0506230102.View ArticlePubMedPubMed CentralGoogle Scholar
- Sotiriou C, Wirapati P, Loi S, Harris A, Fox S, Smeds J, Nordgren H, Farmer P, Praz V, Haibe-Kains B, et al: Gene expression profiling in breast cancer: understanding the molecular basis of histologic grade to improve prognosis. J Natl Cancer Inst. 2006, 98: 262-272.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Sorlie T, Tibshirani R, Parker J, Hastie T, Marron JS, Nobel A, Deng S, Johnsen H, Pesich R, Geisler S, et al: Repeated observation of breast tumor subtypes in independent gene expression data sets. Proc Natl Acad Sci USA. 2003, 100: 8418-8423. 10.1073/pnas.0932692100.View ArticlePubMedPubMed CentralGoogle Scholar
- West M, Blanchette C, Dressman H, Huang E, Ishida S, Spang R, Zuzan H, Olson JA, Marks JR, Nevins JR: Predicting the clinical status of human breast cancer by using gene expression profiles. Proc Natl Acad Sci USA. 2001, 98: 11462-11467. 10.1073/pnas.201162998.View ArticlePubMedPubMed CentralGoogle Scholar
- Takeda A, Goolsby C, Yaseen NR: NUP98-HOXA9 induces long-term proliferation and blocks differentiation of primary human CD34+ hematopoietic cells. Cancer Res. 2006, 66: 6628-6637. 10.1158/0008-5472.CAN-06-0458.View ArticlePubMedGoogle Scholar
- Han W, Kim S, Ae Park I, Kang D, Kim S-W, Youn Y-K, Oh S, Choe K, Noh D-Y: Young age: an independent risk factor for disease-free survival in women with operable breast cancer. BMC Cancer. 2004, 4: 82-89. 10.1186/1471-2407-4-82.View ArticlePubMedPubMed CentralGoogle Scholar
- Jayasinghe UW, Taylor R, Boyages J: Is age at diagnosis an independent prognostic factor for survival following breast cancer?. ANZ J Surg. 2005, 75: 762-767. 10.1111/j.1445-2197.2005.03515.x.View ArticlePubMedGoogle Scholar
- Oh JL, Bonnen M, Outlaw ED, Schechter NR, Perkins GH, Strom EA, Babiera G, Oswald MJ, Allen PK, Thames HD, et al: The impact of young age on locoregional recurrence after doxorubicin-based breast conservation therapy in patients 40 years old or younger: how young is 'young'?. Int J Radiat Oncol Biol Phys. 2006, 65: 1345-1352. 10.1016/j.ijrobp.2006.03.028.View ArticlePubMedGoogle Scholar
- Bortoli S, Renault V, Eveno E, Auffray C, Butler-Browne G, Pietu G: Gene expression profiling of human satellite cells during muscular aging using cDNA arrays. Gene. 2003, 321: 145-154. 10.1016/j.gene.2003.08.025.View ArticlePubMedGoogle Scholar
- Bergamaschi A, Kim YH, Wang P, Sorlie T, Hernandez-Boussard T, Lonning PE, Tibshirani R, Borresen-Dale AL, Pollack JR: Distinct patterns of DNA copy number alteration are associated with different clinicopathological features and gene-expression subtypes of breast cancer. Genes Chromosomes Cancer. 2006, 45: 1033-1040. 10.1002/gcc.20366.View ArticlePubMedGoogle Scholar
- van 't Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AAM, 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.View ArticlePubMedGoogle Scholar
- Dai H, van't Veer L, Lamb J, He YD, Mao M, Fine BM, Bernards R, van de Vijver M, Deutsch P, Sachs A, et al: A cell proliferation signature is a marker of extremely poor outcome in a subpopulation of breast cancer patients. Cancer Res. 2005, 65: 4059-4066. 10.1158/0008-5472.CAN-04-3953.View ArticlePubMedGoogle Scholar
- Rae JM, Johnson MD, Scheys JO, Cordero KE, Larios JM, Lippman ME: GREB1 is a critical regulator of hormone dependent breast cancer growth. Breast Cancer Res Treat. 2005, 92: 141-149. 10.1007/s10549-005-1483-4.View ArticlePubMedGoogle Scholar
- Ciarloni L, Mallepell S, Brisken C: Amphiregulin is an essential mediator of estrogen receptor α function in mammary gland development. Proc Natl Acad Sci USA. 2007, 104: 5455-5460. 10.1073/pnas.0611647104.View ArticlePubMedPubMed CentralGoogle Scholar
- Melov S, Tarnopolsky M, Beckman K, Felkey K, Hubbard A: Resistance exercise reverses aging in human skeletal muscle. PLoS ONE. 2007, 2: e465-10.1371/journal.pone.0000465.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen H, Sukumar S: Role of homeobox genes in normal mammary gland development and breast tumorigenesis. J Mammary Gland Biol Neoplasia. 2003, 8: 159-175. 10.1023/A:1025996707117.View ArticlePubMedGoogle Scholar
- Rubin E, Wu X, Zhu T, Cheung JCY, Chen H, Lorincz A, Pandita RK, Sharma GG, Ha HC, Gasson J, et al: A role for the HOXB7 homeodomain protein in DNA repair. Cancer Res. 2007, 67: 1527-1535. 10.1158/0008-5472.CAN-06-4283.View ArticlePubMedGoogle Scholar
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.