Skip to main content

Breast tumor copy number aberration phenotypes and genomic instability

Abstract

Background

Genomic DNA copy number aberrations are frequent in solid tumors, although the underlying causes of chromosomal instability in tumors remain obscure. Genes likely to have genomic instability phenotypes when mutated (e.g. those involved in mitosis, replication, repair, and telomeres) are rarely mutated in chromosomally unstable sporadic tumors, even though such mutations are associated with some heritable cancer prone syndromes.

Methods

We applied array comparative genomic hybridization (CGH) to the analysis of breast tumors. The variation in the levels of genomic instability amongst tumors prompted us to investigate whether alterations in processes/genes involved in maintenance and/or manipulation of the genome were associated with particular types of genomic instability.

Results

We discriminated three breast tumor subtypes based on genomic DNA copy number alterations. The subtypes varied with respect to level of genomic instability. We find that shorter telomeres and altered telomere related gene expression are associated with amplification, implicating telomere attrition as a promoter of this type of aberration in breast cancer. On the other hand, the numbers of chromosomal alterations, particularly low level changes, are associated with altered expression of genes in other functional classes (mitosis, cell cycle, DNA replication and repair). Further, although loss of function instability phenotypes have been demonstrated for many of the genes in model systems, we observed enhanced expression of most genes in tumors, indicating that over expression, rather than deficiency underlies instability.

Conclusion

Many of the genes associated with higher frequency of copy number aberrations are direct targets of E2F, supporting the hypothesis that deregulation of the Rb pathway is a major contributor to chromosomal instability in breast tumors. These observations are consistent with failure to find mutations in sporadic tumors in genes that have roles in maintenance or manipulation of the genome.

Peer Review reports

Background

Genomic DNA copy number aberrations are frequent in solid tumors [1]. The wide range in the number and types of chromosome level alterations are likely to reflect the different solutions taken by individual tumors to escape normal protective mechanisms. Thus, the spectrum of alterations is likely to reflect a composite of selection and particular failures in genome surveillance mechanism(s). The interplay between selection and genetic instability in shaping tumor genomes is currently most clearly established in tumors with defects in mismatch repair. These tumors have a high frequency of nucleotide sequence level aberrations, fewer DNA copy number alterations and characteristic histological phenotype [1]. On the other hand, less is known about specific gene defects that give rise to chromosome level aberrations in tumors. Mutations in genes encoding proteins involved in mitosis and DNA damage sensing and repair mechanisms, which are associated with chromosomal level instability have been identified in cancer-prone syndromes, including ATM, TP53, BRCA1, BRCA2, NBS1 and BUB1B, however they are rarely mutated in sporadic tumors [2, 3]. Similarly, searches for mutations in genes that participate in maintenance or manipulation of the genome (e.g. genes involved in DNA repair, replication, spindle checkpoints etc.) have found only a small number of mutations in tumors [3]. Nevertheless, deregulation of functions that maintain genome stability appears to occur early in tumors, as activation of the DNA damage checkpoint, possibly in response to DNA replication stress, is evident in pre-malignant lesions [4, 5]. Similarly, telomere shortening is observed in pre-malignant lesions, supporting a role for telomere dysfunction early in tumor development [6]. Other proposed routes to instability include deregulation of CCNE1 and AURKA expression through loss of function of FBXW7 (hCdc4) [7] and more global alteration in gene expression due to deregulation of the Rb pathway [8]. The foregoing discussion suggests that failures in a number of different processes that maintain genome integrity could contribute to the wide variety of genomic alterations in solid tumors. Often these aberrations include net gain or loss of whole chromosomes (aneuploidy) or parts of chromosomes. Gene amplification, defined as a copy number increase of a restricted region of a chromosome arm may also occur. Here we investigated the numbers and types of copy number alterations in tumors and whether they were associated with differential expression of genes likely to play a role in manipulation or maintenance of the genome. These studies found three subtypes of breast tumors distinguished by copy number aberrations. Telomere dysfunction was implicated in the propensity to amplify, since shorter telomeres and differential expression of genes involved in telomere maintenance were associated with the numbers of amplicons and the presence of at least one amplicon, respectively. On the other hand, the number of lower magnitude gains and losses of chromosomal segments was associated with differential expression of genes involved in processes maintaining or manipulating the genome. These genes are significantly enriched for the known targets of E2F. Furthermore, we observed enhanced expression of most E2F target genes, indicating that over expression rather than deficiency was associated with genetic instability. These observations support the hypothesis that deregulation of the Rb/E2F pathway is a major contributor to chromosomal instability in breast tumors.

Methods

Specimens

Frozen tumor tissue was obtained from the University of California San Francisco Comprehensive Cancer Center Breast Oncology Program Tissue Bank. All specimens were collected under approved protocols from UCSF with patient consent. Patient characteristics are provided in Supplementary Table 1 (Additional file 1). Expression and copy number data from a second set of ductal invasive breast tumors were used and patient characteristics are given in Chin et al. (submitted). The patient groups in both sets were similar in terms of their genomic and pathological characterization.

Extraction of nucleic acids

Nucleic acids were extracted from tumor blocks as described previously [9, 10]. Briefly, blocks were trimmed with a razor blade to remove normal tissue and cryosections were obtained from either side of the block to ascertain that tumor cells comprised greater than 70% of the specimen. DNA was extracted using the QUIamp tissue kits (29304, Qiagen).

TP53sequencing

Exons 5–8 of TP53 were amplified from genomic DNA and cycle sequencing was carried out as described previously [11].

Array CGH and data processing

Array CGH, imaging and data analysis were carried out as described previously using arrays of 2464 genomic clones (BAC or P1) each printed in triplicate (HumArray1.14 and HumArray2.0) [11, 12]. Data processing is described in detail in the Supplementary Methods (Additional file 2) and the array data are available in Supplementary Table 2 (Additional file 3).

Telomere length assessment

The mean TRF length was measured using the TeloTAGGG telomere length assay kit (Roche Applied Science). Briefly, 1 μg genomic DNA was digested with Hinf I and Rsa I restriction enzymes and electrophoretically resolved on 0.8% agarose/1X TAE. The gels were blotted to a nylon membrane (Positive charged, Roche) and fixed by UV-crosslinking. After hybridization with digoxigenin labeled telomere specific probe, the signals were visualized with an alkaline phosphatase – CDP-Star chemiluminescent system. The filters were exposed to X-ray film and the mean TRF length was calculated using Quantity One software.

Statistical methods

A detailed description of the methods used for all aspects of the data analysis is provided in the Supplementary Methods (Additional file 2).

Results

Genomic analysis of breast tumors

Application of array CGH to the analysis of copy number aberrations in 62 sporadic ductal invasive breast tumors and five BRCA1 mutant tumors revealed a number of frequent low level gains and losses (Supplementary Tables 1 and 2, Additional files 1 and 3) and 12 regions of recurrent amplification (Table 1). We observed significant differences in the spectrum of aberrations with respect to estrogen receptor (ER) status (Figure 1) consistent with other published reports [13]. We also found that tumors with mutations in exons 5–8 of TP53 showed a higher frequency of alterations, as well as significant differences in the frequency of alteration of certain regions of the genome compared to tumors without mutations as indicated by the t-statistic for each clone (Figure 2). Moreover, we found a highly significant concordance between the test statistic for association of the particular chromosomal regions with TP53 mutation status in our data and an independent set of more than 100 primary breast tumors (Fedele et al., personal communication) (Pearson correlation of 0.53 corresponding to the p-value << 0.001).

Table 1 Recurrent amplicons in breast tumors with examples of some candidate oncogenes
Figure 1
figure 1

Frequency plot of copy number alterations in ER positive and negative tumors. The top two panels show the frequency of gains, indicated by the green bars ranging from 0 to 1, and losses, indicated by the red bars ranging from 0 to -1, in 62 sporadic breast tumors for each clone. The bottom panel displays the magnitude of the t-statistic for each clone computed based on the smoothed data as described in the Methods. The horizontal dotted lines indicate the statistic cut-off corresponding to the FDR-adjusted p-value of 0.05 (blue) and 0.1 (green).

Figure 2
figure 2

Analysis of TP53 mutation in breast tumors. Frequency plot of copy number changes in TP53 mutant and wild type tumors. The top two panels show the frequency of gains, indicated by the green bars ranging from 0 to 1, and losses, indicated by the red bars ranging from 0 to -1, in 62 sporadic breast tumors for each clone. The bottom panel displays the magnitude of the t-statistic for each clone computed based on the smoothed data as described in the Methods. The horizontal dotted lines indicate the statistic cut-off corresponding to the FDR-adjusted p-value of 0.01 (red), 0.05 (blue) and 0.1 (green). Thus, copy number alterations occurring more frequently in TP53 mutant tumors included losses of regions on 3p, 4q, 5q, 15q, 17q and gain of a small region on distal 8q.

Hierarchical clustering of tumors according to their genome-wide DNA copy number profiles revealed three main branches. Tumors within each of the branches also differed in the number of copy number changes that were present, as well as the frequency of particular aberrations (Figures 3 and 4). The same three clusters were observed in an independent set of breast tumors (Chin et al., submitted), thus confirming our initial observation. The groups in both sets agreed in terms of their genomic and pathological characterization.

Figure 3
figure 3

Genomic analysis of breast tumors reveals three subtypes. A. Hierarchical clustering of 62 ductal invasive breast tumors and five BRCA1 mutant tumors based on their genome-wide DNA copy number profiles. Individual clones are represented as rows, ordered by chromosome and genome position according to the July 2003 freeze of the human genome. Clones on the p-arm and q-arm of chromosomes are indicated in shades of dark gray (odd numbered chromosomes) or light gray (even numbered chromosomes). Acrocentric chromosomes are shown in dark or light gray. Columns represent individual tumor samples. The estrogen receptor status of the tumors is shown in shades of blue (dark blue = ER negative, light blue = ER positive), BRCA1 mutant tumors are indicated in orange, and TP53 mutation status is indicated with a maroon box for TP53 mutant tumors, a gray box for tumors with no detected mutation and a black box if the TP53 status is unknown. Copy number losses are indicated in red, gains in green and amplifications as yellow dots. Three main clusters are evident. B-F. Genome-wide copy number aberrations profiles of sporadic and hereditary (BRCA1) breast tumors are plotted as the normalized log2ratio for each clone sorted by chromosome and ordered according to genome position from the p-arm to the q-arm. Normalized copy number ratios of genomic DNA are shown for a tumor from the 1q/16q cluster with few copy number changes including gain of 1q and loss of 16q (B), a tumor from the ER negative, complex cluster showing many low level chromosome changes and few amplifications (C), a tumor from the amplifier cluster with low level gains and losses and amplifications (D) and tumors from patients with mutations in BRCA1 (E and F). G. Numbers and types of copy number aberrations in breast tumor subtypes. The mean numbers of whole chromosome copy number changes, copy number transitions and amplifications were determined for the tumors within each subtype. H. Numbers and types of copy number aberrations in breast tumor subtypes. The mean numbers of chromosomes showing no copy number change, whole chromosome copy number changes, copy number transitions, copy number transitions at centromeres and amplifications were determined for the tumors within each subtype. I. Association with disease-specific survival. Significance of the log-rank test was used to assess the association between a genomic subclass and survival phenotypes. The significance was declared at p < 0.05. Patients with complex tumors experienced significantly worse outcome compared to the other groups.

Figure 4
figure 4

Copy number changes more frequently associated with one subtype. Frequency for each clone of gains and losses, which were uniquely present in more than 50% of samples of one subtype and in less than 30% of samples in other subtypes. Gains are indicated by green bars, ranging from 0 to 1, and losses, by the red bars ranging from 0 to -1 for each clone.

At the low end of chromosomal level instability are ER positive tumors (n = 7 tumors, Figure 1A, left branch), designated 1q/16q, as their genomes showed very few copy number changes other than gain of 1q and loss of 16q (Figure 3A, B, G and 3H, Figure 4). Tumors in this group were exclusively of moderately or well differentiated grade, stage II, and did not recur. These tumors had very high within group similarity with average pairwise Pearson correlation of 0.76.

At the other extreme of genome instability are sporadic tumors (n = 16, Figure 3A, middle branch) in which the mean fraction of the genome at altered copy number is greater than 0.6 due to the presence of many low level copy number aberrations (Figure 3A, C, G and 3H). Copy number losses involving chromosomes 3p, 4, 5q, 11p, 14q, 15q, 17q and 18q were more prevalent in this subtype than in others (Figure 4). All but one tumor in this group are ER negative, all were of high grade and patients experienced significantly worse outcome as compared to other groups (Figure 3I). Four had mutations in exons 5–8 of TP53, accompanied by a copy number loss encompassing the locus (Table 2). In addition this cluster contains all of the familial BRCA1 mutant tumors (Figure 3E and 3F) in our dataset. Similar to BRCA1 mutant tumors, they show a relatively high degree of within group similarity in regions of aberration in spite of the presence of many alterations (see Methods for discussion of statistical analysis, Additional file 2). We refer to this group as "complex" in recognition of their many low level copy number alterations.

Table 2 TP53 mutations in breast tumors

The third group (n = 39, Figure 3A, right branch) comprised of both ER positive and negative tumors is characterized by the presence of low level gains and losses and recurrent amplifications (Figure 3A, D, G and 3H). Gains involving chromosome 8q and 12p were more frequent in this group than the others (Figure 4). We refer to this group as "mixed amplifiers." The more frequently occurring amplifications in this group, which occurred predominantly in the ER positive tumors involved 8p, including FGFR1 (6 ER positive tumors/6 tumors with 8p amplification), 11q13, including CCND1 (8 ER positive tumors/9 tumors with 11q13 amplification) and regions of 20q including ZNF217 (6 ER positive tumors/9 tumors with 20q amplification). Within this subtype, amplification of 17q (ERBB2) was present in both ER positive (n = 5) and ER negative tumors (n = 2) (Figure 5A and Table 1).

Figure 5
figure 5

Overview of most frequent DNA amplifications on chromosomes 8, 11, 17 and 20 in 62 breast tumors determined by genome-wide array CGH. A. Heat map depiction of aberrations on chromosomes 8, 11, 17 and 20 and typical chromosome copy number profiles showing amplifications of 8p (including FGFR1), 11q13 (including CCND1), 17q (ERBB2) and regions on 20q (including ZNF217). Note that the chromosome 11 copy number profiles vary depending on whether amplification of chromosome 8 is also present. In both cases copy number losses distal to the amplicon are observed, however in the absence of chromosome 8 amplification (right), the region of loss extends distally from the amplified region, whereas, when chromosome 8 was amplified (left), the copy number loss includes regions proximal as well as distal to CCND1. B. For each tumor, the numbers of copy number transitions is compared to the number of chromosome arms with at least one amplification.

Association of copy number aberration types with alterations in processes/genes involved in maintenance and manipulation of the genome

The discrimination of breast tumor subtypes based on copy number aberrations led us to investigate possible associations of copy number aberration types with alterations in processes/genes involved in maintenance of genome stability. Over expression or depletion of such genes in vitro results in a variety of genome instability phenotypes, including disruptions of chromosome integrity, aberrant mitoses, aberrant cell division, etc. As telomere dysfunction has been widely proposed as a source of genetic instability in tumors, we first investigated the possible association of telomere attrition with copy number aberrations. We determined average telomere length in 28 breast tumors using Southern blotting. We found an inverse correlation between telomere length and number of chromosome arms with amplification (Figure 6, Spearman correlation = -0.42, p = 0.02). Moreover we observed an inverse association between telomere length and the presence of at least one amplicon (median length in amplified samples of 6.3 compared to unamplified, 7.4), but the comparison was underpowered and statistically not significant (Wilcoxon rank sum test, p = 0.25). These observations suggest a role for telomere attrition in promoting amplification in breast tumors.

Figure 6
figure 6

Telomere length measurements in 28 breast tumor samples. Plotted is the telomere length determined by Southern blotting relative to the number of chromosome arms with amplification. A significant inverse correlation is evident.

Next, we investigated whether expression levels of genes that play a role in maintenance or manipulation of the genome varied among tumors with greater or lesser numbers of copy number aberrations. To carry out this analysis we used a second independent set of 101 ductal invasive breast tumors for which copy number profiles and Affymetrix High Throughput Array (HTA) GeneChip® expression data were available (Chin et al., submitted). We determined the number and type of copy number changes in each tumor by counting three types of copy number alterations; copy number changes involving whole chromosomes, low level gains and losses affecting extended portions of chromosomes, and amplifications defined as focal regions of increased copy number [14]. Specifically, a clone was declared amplified if it belonged to a copy number segment <20 Mb and the increase in ratio exceeded the criterion described in the Statitical Methods. The distinction between gains and amplifications can be seen in the copy number profiles in Figure 5A. A copy number gain spanning 8q can be seen in the top left profile, while the wide variety in amplicon profiles is evident by comparison of all the profiles. We enumerate low level changes by counting "copy number transitions," the number of changes in the CGH profile from one copy number level to another that occur within chromosomes (see Supplementary Methods for further discussion of aberration finding, Additional file 2). Since the spacing between clones is ~1.5 Mb, focal aberrations that fall between clones on the array will be missed. On the other hand, all copy number transitions will be recorded, but the precision with which they will be located on the genome will depend on clone spacing. We note that these copy number analyses found that the number of copy number transitions associated with amplifications varied over a wide range in tumors of all subtypes in both datasets, however the greatest number of amplifications did not occur in the samples with either the smallest or largest number of copy number transitions (Figure 5B).

We tested for associations between gene expression and copy number aberrations by developing a list of 426 genes assigned to functional categories, "DNA replication," "DNA damage/repair," "cell cycle," "mitosis," "centrosome" (centrosome and centrosome cycle) and "telomere" using Gene Ontology Annotation (GOA) terms and reference to the literature (Supplementary Methods and Supplementary Table 3, Additional files 2 and 4). Many of the genes were assigned to more than one process. Expression of 350 of the 426 stability genes could be analyzed in the breast tumor data set. Controlling for estrogen receptor status as a possible confounder, we observed that the telomere functional class was significantly associated with the presence of at least one amplicon, with half of the genes showing positive and half showing negative association (Table 3). On the other hand, we found that mitosis, cell cycle, DNA replication and DNA damage/repair functional classes were highly significantly enriched for association with copy number transitions (Table 3). Moreover these associations held when only tumors of the complex subtype were considered, indicating that the results are not due to confounding between the large number of low level transitions and complex subtype. Expression of 146 stability genes was significantly associated with the number of copy number transitions (false discovery rate, FDR < 0.05). Most associations were positive (120/146 genes), indicating that enhanced expression of these genes was associated with greater numbers of copy number transitions (Figure 7). The number of amplifications was associated with mitosis, cell cycle and to a lesser extent, DNA replication categories. Again, this relation held when considering all samples, as well as only samples within the mixed amplifier subtype. Twenty-five individual genes were associated with number of amplifications (FDR < 0.05) and 21 were in common with the group of genes associated with copy number transitions (Figure 7). Here too, we observed that most associations were positive (21/25 genes), indicating that increased expression of the genes was observed in tumors with more amplifications. Finally, we investigated how the GOA categories represented by our list of 426 stability genes ranked among all known GOA categories with respect to associations with frequency of copy number aberrations by considering all probes measured by the expression analysis. In this subsequent unsupervised analysis, we found enrichment for genes associated with copy number transitions (FDR < 0.05) in the same functional classes, e.g. mitosis, cell cycle, cell division and DNA replication (Holm adjusted p-value < 0.005), providing further support for these associations.

Table 3 Association of expression of functional classes with copy number aberration types
Figure 7
figure 7

Signaling pathway and E2F responsive genes associated with numbers of copy number transitions and amplifications. The receptor signaling cascade impinging on Rb includes a number of up-stream genes that are frequently altered in cancers and result in deregulation of the E2F family of transcription factors through inhibition or loss of Rb function. Here, we show genes enriched for association with numbers of copy number transitions (FDR < 0.05). Genes, which have been identified as E2F targets (asterisk) and/or are correlated with E2F1 expression levels (|Pearson correlation| > 0.3) are shown in the top box. Genes which showed less correlation with E2F1 are shown in the bottom box. Genes also associated with numbers of amplifications (FDR < 0.05) are highlighted in yellow. Four additional genes associated with amplifications (FDR < 0.05) are ASCIZ, MSH5, RAE1 and TDG. Red, increased expression. Blue, decreased expression.

We noted that the 146 stability genes associated with numbers of copy number transitions included E2F1, and they are significantly enriched for genes known to be targets of E2F1 (p < 2 × 10-6, Fisher exact test, Figure 7). Moreover the expression levels of known E2F1 target genes were highly correlated with E2F1 expression (p < 2 × 10-10, Supplementary Table 4, Additional file 5). These observations provide in vivo validation of the in vitro determinations of E2F1 target genes. They are also consistent with deregulation of E2F being a major contributor to genomic instability affecting numbers of copy number transitions and amplifications. Taken together these observations suggest that telomere attrition and deregulated expression of genes in the other functional classes, particularly those that are targets of E2F, contribute to the numbers of chromosomal alterations.

Discussion

Our analysis of large numbers of breast tumors by array CGH revealed variety in the numbers and types of copy number alterations in the tumor genomes. In the ductal invasive breast tumors reported here, three subtypes were distinguished by copy number alterations. The subtypes differed with respect to the numbers and types of aberrations, as well as patient survival. The1q/16q subtype with very few copy number alterations in addition to gain of 1q and loss of 16q was associated with the best patient outcome, consistent with other studies. Searches for tumor suppressor gene(s) on 16q have failed to find mutations in candidate genes in the region in ductal invasive breast cancer, although mutations in E cadherin and loss of 16q are characteristic of lobular breast tumors. Two genes involved in telomere maintenance, TERF2 and TERF2IP were among those ruled out as tumor suppressors on 16q, as was E2F4 [15–17]. The stability of the genome of these tumors also suggests that copy number alterations of these and other stability genes mapping within the aberrant regions, +1q and -16q are less likely to contribute to chromosomal level instability in breast cancer.

Complex tumors with extensive chromosomal level instability were associated with poor patient survival. They are similar to BRCA1 hereditary tumors in their copy number alterations [18, 19] (Figure 3). BRCA1 participates in a number of cell functions that maintain genome integrity either directly through double strand break repair or indirectly through maintenance of checkpoints at G1, S and mitosis [20–22]. Thus, it is possible that BRCA1 [23, 24] or the genes/pathways that interact with BRCA1 are defective in this subtype either through mutation, silencing or copy number mediated dosage effects. We note that the copy number loss on 17q associated with this subtype includes the BRCA1 locus (9/16 tumors, Figure 4).

The discrimination of breast tumor subtypes based on copy number aberrations led us to investigate possible associations of copy number aberration types with alterations in processes/genes involved in maintenance of genome stability. We observed shorter telomeres in tumors with greater numbers of amplifications, consistent with telomere attrition promoting this type of copy number aberration in breast tumors. Telomere dysfunction, often referred to as "telomere crisis" has been implicated in amplification, particularly by breakage-fusion-bridge processes. On the other hand, our analyses of stability gene expression in relation to copy number aberration types found that expression of genes in the functional classes; "mitosis," "cell cycle," "replication," and "DNA damage/repair" were associated with greater numbers of copy number transitions. Furthermore, a subsequent analysis found significant enrichment for these same classes among all GOA groups when analyzed with GOStats [25]. The number of amplicons was associated with similar functional groups, "mitosis" and "cell cycle." Many of these genes are E2F targets [26–36] and therefore potentially coordinately deregulated due to Rb pathway defects [37]. Abrogation of Rb pathway function is frequent in breast tumors by loss of expression of Rb or altered expression of inhibitors of Rb activity (e.g. loss/silencing of CDKN2A (p16) and amplification and/or over expression of CCND1, CDK4, CDK6) (Figure 7). It is interesting to note that whereas E2F1 is up-regulated in breast tumors, its expression is low in prostate tumors [38], which typically have genomes with fewer copy number changes than most ductal invasive breast cancers [39]. For example, in an array CGH dataset of 64 primary prostate tumor samples [39], the median number of copy number transitions was 13 per tumor genome compared to 30 in our primary breast tumor samples (p < 5 × 10-9, Wilcoxon rank sum test). Mechanistic support for a central role of E2F1 in genomic instability comes from a recent report that elevated numbers of DNA double strand breaks are present in cell lines with deregulated E2F1 and Rb deficiency [40].

Chromosomal instability has been observed in vitro when many of these E2F target genes (Figure 7) associated with replication, DNA repair, cell cycle control and the mitotic checkpoint are mutated, knocked out or knocked down using siRNA [8, 41, 42]. Contrary to expectation, we observed that greater chromosomal instability in breast tumors is associated with increased expression levels of many of these genes, even though they have loss of function instability phenotypes. These assays further demonstrate that loss of a single copy of some of the genes results instability or cancer prone phenotypes. Genes that have been shown to be haploinsufficient in this way and that are among those we identified as showing significant association with the number of copy number aberrations in our tumors (FDR < 0.05) include RAD17, ATM and RB1, which are expressed at lower levels in tumors with more copy number changes. These genes are also negatively correlated with E2F1 expression. Other genes showing haploinsufficiency in vitro, MAD2L1, PLK4, BUB1B and CHEK1 show enhanced expression in association with number of chromosomal changes and are positively correlated with E2F1 expression (Supplementary Table 4, Additional file 5). As all seven of the above mentioned genes with haploinsufficiency phenotypes map to regions of frequent loss in breast tumors and genetic instability phenotypes are associated with deficiency in these genes, we asked whether loss of function might play a role in the subset of tumors in which there is a copy number loss of the locus. Specifically, we asked if their expression levels were down regulated when there is a copy number loss. Although 118 of the genome stability genes showed highly significant reduction in expression in tumors in which the locus was lost (FDR < 0.05, one-sided Wilcoxon rank sum test), we found little difference in expression level with copy number loss for MAD2L1, PLK4, ATM and RB1, whereas BUB1B was increased in expression in tumors with loss of the locus (Supplementary Table 4, Additional file 5). Only expression of RAD17 was significantly reduced when lost (unadjusted p = 8 × 10-4, Wilcoxon rank sum test), suggesting that RAD17 might be haploinsufficient in tumors with copy number loss of the locus at 5q13.

Our observations in tumors support the hypothesis that global alteration of expression of genes involved in processes such as chromosome segregation and maintenance of genome integrity, driven by deregulation of E2F, underlies much of the chromosomal instability in breast tumors. Furthermore gene expression appears to be relatively up-regulated. On the one hand, this observation seems contradictory in light of the phenotypes resulting from mutational analyses of genes involved in maintenance of genome stability. Such in vitro studies have generally assessed the consequences of functional deficiency one gene at a time and have found that individually many genes have loss of function instability phenotypes. On the other hand, as many of these genes participate in multi-protein complexes that depend on proper stoichiometry for function, alterations resulting in overproduction or deficiency are likely to have similar or related phenotypes (reviewed in [43]). Indeed, in mammalian cells, instability phenotypes have been reported in association with both up and down regulation of genes such as MAD2L1 [8, 41], ATR [44, 45], PLK4 [46] and AURKA [47]. Further studies will be required not only to assess instability phenotypes when expression levels are increased, but also how phenotypes might vary when multiple genes are up-regulated.

In tumors, changes in gene dosage due to low level copy number alterations may also lead to small alterations in expression of multiple genes, which together could contribute to dysfunction of processes manipulating the genome, resulting in more error prone cell division cycles. Thus, during tumor progression, genome instability may be enhanced not only by deregulation of E2F, but also by the acquisition of greater numbers of copy number changes encompassing more genes involved in genome maintenance. Since genetic instability is an on-going feature of tumors, allowing them to evolve resistance to therapy, the ability to recognize the active mechanisms of instability in tumors may help to guide therapeutic decisions.

Conclusion

Application of array CGH to the study of breast tumors found three subtypes. Investigation of the numbers and types of copy number alterations in tumors and their association with differential expression of genes likely to play a role in manipulation or maintenance of the genome implicated telomere dysfunction in the propensity to amplify. On the other hand, the number of lower magnitude gains and losses of chromosomal segments was associated with differential expression of genes which were significantly enriched for the known targets of E2F, supporting the hypothesis that deregulation of E2F underlies much of the chromosomal instability in breast tumors. Furthermore, we observed enhanced expression of most E2F target genes, indicating that over expression rather than deficiency was associated with genetic instability. These observations provide a possible explanation for the failure to find mutations in sporadic tumors in genes that have roles in maintenance or manipulation of the genome.

References

  1. Albertson DG, Collins C, McCormick F, Gray JW: Chromosome aberrations in solid tumors. Nat Genet. 2003, 34: 369-376. 10.1038/ng1215.

    Article  CAS  PubMed  Google Scholar 

  2. Cahill DP, Lengauer C, Yu J, Riggins GJ, Willson JK, Markowitz SD, Kinzler KW, Vogelstein B: Mutations of mitotic checkpoint genes in human cancers. Nature. 1998, 392: 300-303. 10.1038/32688.

    Article  CAS  PubMed  Google Scholar 

  3. Wang Z, Cummins JM, Shen D, Cahill DP, Jallepalli PV, Wang TL, Parsons DW, Traverso G, Awad M, Silliman N, Ptak J, Szabo S, Willson JK, Markowitz SD, Goldberg ML, Karess R, Kinzler KW, Vogelstein B, Velculescu VE, Lengauer C: Three classes of genes mutated in colorectal cancers with chromosomal instability. Cancer Res. 2004, 64: 2998-3001. 10.1158/0008-5472.CAN-04-0587.

    Article  CAS  PubMed  Google Scholar 

  4. Bartkova J, Horejsi Z, Koed K, Kramer A, Tort F, Zieger K, Guldberg P, Sehested M, Nesland JM, Lukas C, Orntoft T, Lukas J, Bartek J: DNA damage response as a candidate anti-cancer barrier in early human tumorigenesis. Nature. 2005, 434: 864-870. 10.1038/nature03482.

    Article  CAS  PubMed  Google Scholar 

  5. Gorgoulis VG, Vassiliou LV, Karakaidos P, Zacharatos P, Kotsinas A, Liloglou T, Venere M, Ditullio RA, Kastrinakis NG, Levy B, Kletsas D, Yoneta A, Herlyn M, Kittas C, Halazonetis TD: Activation of the DNA damage checkpoint and genomic instability in human precancerous lesions. Nature. 2005, 434: 907-913. 10.1038/nature03485.

    Article  CAS  PubMed  Google Scholar 

  6. Feldser DM, Hackett JA, Greider CW: Telomere dysfunction and the initiation of genome instability. Nat Rev Cancer. 2003, 3: 623-627. 10.1038/nrc1142.

    Article  CAS  PubMed  Google Scholar 

  7. Mao JH, Perez-Losada J, Wu D, Delrosario R, Tsunematsu R, Nakayama KI, Brown K, Bryson S, Balmain A: Fbxw7/Cdc4 is a p53-dependent, haploinsufficient tumour suppressor gene. Nature. 2004, 432: 775-779. 10.1038/nature03155.

    Article  CAS  PubMed  Google Scholar 

  8. Hernando E, Nahle Z, Juan G, Diaz-Rodriguez E, Alaminos M, Hemann M, Michel L, Mittal V, Gerald W, Benezra R, Lowe SW, Cordon-Cardo C: Rb inactivation promotes genomic instability by uncoupling cell cycle progression from mitotic control. Nature. 2004, 430: 797-802. 10.1038/nature02820.

    Article  CAS  PubMed  Google Scholar 

  9. Albertson DG, Ylstra B, Segraves R, Collins C, Dairkee SH, Kowbel D, Kuo WL, Gray JW, Pinkel D: Quantitative mapping of amplicon structure by array CGH identifies CYP24 as a candidate oncogene. Nat Genet. 2000, 25: 144-146. 10.1038/75985.

    Article  CAS  PubMed  Google Scholar 

  10. Pinkel D, Segraves R, Sudar D, Clark S, Poole I, Kowbel D, Collins C, Kuo WL, Chen C, Zhai Y, Dairkee SH, Ljung BM, Gray JW, Albertson DG: High resolution analysis of DNA copy number variation using comparative genomic hybridization to microarrays. Nat Genet. 1998, 20: 207-211. 10.1038/2524.

    Article  CAS  PubMed  Google Scholar 

  11. Snijders AM, Schmidt BL, Fridlyand J, Dekker N, Pinkel D, Jordan RC, Albertson DG: Rare amplicons implicate frequent deregulation of cell fate specification pathways in oral squamous cell carcinoma. Oncogene. 2005, 24: 4232-4242. 10.1038/sj.onc.1208601.

    Article  CAS  PubMed  Google Scholar 

  12. Snijders AM, Nowak N, Segraves R, Blackwood S, Brown N, Conroy J, Hamilton G, Hindle AK, Huey B, Kimura K, Law S, Myambo K, Palmer J, Ylstra B, Yue JP, Gray JW, Jain AN, Pinkel D, Albertson DG: Assembly of microarrays for genome-wide measurement of DNA copy number. Nat Genet. 2001, 29: 263-264. 10.1038/ng754.

    Article  CAS  PubMed  Google Scholar 

  13. Loo LW, Grove DI, Williams EM, Neal CL, Cousens LA, Schubert EL, Holcomb IN, Massa HF, Glogovac J, Li CI, Malone KE, Daling JR, Delrow JJ, Trask BJ, Hsu L, Porter PL: Array comparative genomic hybridization analysis of genomic alterations in breast cancer subtypes. Cancer Res. 2004, 64: 8541-8549. 10.1158/0008-5472.CAN-04-1992.

    Article  CAS  PubMed  Google Scholar 

  14. Snijders AN, Fridlyand J, Mans D, Segraves R, Jain AN, Pinkel D, Albertson DG: Shaping tumor and drug resistant genomes by instability and selection. Oncogene. 2003, 22: 4370-4379. 10.1038/sj.onc.1206482.

    Article  CAS  PubMed  Google Scholar 

  15. Rakha EA, Armour JA, Pinder SE, Paish CE, Ellis IO: High-resolution analysis of 16q22.1 in breast carcinoma using DNA amplifiable probes (multiplex amplifiable probe hybridization technique) and immunohistochemistry. Int J Cancer. 2005, 114: 720-729. 10.1002/ijc.20738.

    Article  CAS  PubMed  Google Scholar 

  16. Rakha EA, Pinder SE, Paish EC, Robertson JF, Ellis IO: Expression of E2F-4 in invasive breast carcinomas is associated with poor prognosis. J Pathol. 2004, 203: 754-761. 10.1002/path.1573.

    Article  CAS  PubMed  Google Scholar 

  17. van Wezel T, Lombaerts M, van Roon EH, Philippo K, Baelde HJ, Szuhai K, Cornelisse CJ, Cleton-Jansen AM: Expression analysis of candidate breast tumour suppressor genes on chromosome 16q. Breast Cancer Res. 2005, 7: R998-R1004. 10.1186/bcr1337.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Jonsson G, Naylor TL, Vallon-Christersson J, Staaf J, Huang J, Ward MR, Greshock JD, Luts L, Olsson H, Rahman N, Stratton M, Ringner M, Borg A, Weber BL: Distinct genomic profiles in hereditary breast tumors identified by array-based comparative genomic hybridization. Cancer Res. 2005, 65: 7612-7621.

    PubMed  Google Scholar 

  19. Wessels LF, van Welsem T, Hart AA, van't Veer LJ, Reinders MJ, Nederlof PM: Molecular classification of breast carcinomas by comparative genomic hybridization: a specific somatic genetic profile for BRCA1 tumors. Cancer Res. 2002, 62: 7110-7117.

    CAS  PubMed  Google Scholar 

  20. Scully R, Xie A, Nagaraju G: Molecular functions of BRCA1 in the DNA damage response. Cancer Biol Ther. 2004, 3: 521-527.

    Article  CAS  PubMed  Google Scholar 

  21. Wang RH, Yu H, Deng CX: A requirement for breast-cancer-associated gene 1 (BRCA1) in the spindle checkpoint. Proc Natl Acad Sci U S A. 2004, 101: 17108-17113. 10.1073/pnas.0407585101.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Yoshida K, Miki Y: Role of BRCA1 and BRCA2 as regulators of DNA repair, transcription, and cell cycle in response to DNA damage. Cancer Sci. 2004, 95: 866-871. 10.1111/j.1349-7006.2004.tb02195.x.

    Article  CAS  PubMed  Google Scholar 

  23. Staff S, Isola J, Tanner M: Haplo-insufficiency of BRCA1 in sporadic breast cancer. Cancer Res. 2003, 63: 4978-4983.

    CAS  PubMed  Google Scholar 

  24. Wei M, Grushko TA, Dignam J, Hagos F, Nanda R, Sveen L, Xu J, Fackenthal J, Tretiakova M, Das S, Olopade OI: BRCA1 promoter methylation in sporadic breast cancer is associated with reduced BRCA1 copy number and chromosome 17 aneusomy. Cancer Res. 2005, 65: 10692-10699. 10.1158/0008-5472.CAN-05-1277.

    Article  CAS  PubMed  Google Scholar 

  25. Beissbarth T, Speed TP: GOstat: find statistically overrepresented Gene Ontologies within a group of genes. Bioinformatics. 2004, 20: 1464-1465. 10.1093/bioinformatics/bth088.

    Article  CAS  PubMed  Google Scholar 

  26. Berkovich E, Ginsberg D: ATM is a target for positive regulation by E2F-1. Oncogene. 2003, 22: 161-167. 10.1038/sj.onc.1206144.

    Article  CAS  PubMed  Google Scholar 

  27. Carrassa L, Broggini M, Vikhanskaya F, Damia G: Characterization of the 5'flanking region of the human Chk1 gene: identification of E2F1 functional sites. Cell Cycle. 2003, 2: 604-609.

    Article  CAS  PubMed  Google Scholar 

  28. Ishida S, Huang E, Zuzan H, Spang R, Leone G, West M, Nevins JR: Role for E2F in control of both DNA replication and mitotic functions as revealed from DNA microarray analysis. Mol Cell Biol. 2001, 21: 4684-4699. 10.1128/MCB.21.14.4684-4699.2001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Iwanaga R, Komori H, Ishida S, Okamura N, Nakayama K, Nakayama KI, Ohtani K: Identification of novel E2F1 target genes regulated in cell cycle-dependent and independent manners. Oncogene. 2005

    Google Scholar 

  30. Jiang Y, Saavedra HI, Holloway MP, Leone G, Altura RA: Aberrant regulation of survivin by the RB/E2F family of proteins. J Biol Chem. 2004, 279: 40511-40520. 10.1074/jbc.M404496200.

    Article  CAS  PubMed  Google Scholar 

  31. Kel AE, Kel-Margoulis OV, Farnham PJ, Bartley SM, Wingender E, Zhang MQ: Computer-assisted identification of cell cycle-related genes: new targets for E2F transcription factors. J Mol Biol. 2001, 309: 99-120. 10.1006/jmbi.2001.4650.

    Article  CAS  PubMed  Google Scholar 

  32. Muller H, Bracken AP, Vernell R, Moroni MC, Christians F, Grassilli E, Prosperini E, Vigo E, Oliner JD, Helin K: E2Fs regulate the expression of genes involved in differentiation, development, proliferation, and apoptosis. Genes Dev. 2001, 15: 267-285. 10.1101/gad.864201.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Polager S, Kalma Y, Berkovich E, Ginsberg D: E2Fs up-regulate expression of genes involved in DNA replication, DNA repair and mitosis. Oncogene. 2002, 21: 437-446. 10.1038/sj.onc.1205102.

    Article  CAS  PubMed  Google Scholar 

  34. Ren B, Cam H, Takahashi Y, Volkert T, Terragni J, Young RA, Dynlacht BD: E2F integrates cell cycle progression with DNA repair, replication, and G(2)/M checkpoints. Genes Dev. 2002, 16: 245-256. 10.1101/gad.949802.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Vernell R, Helin K, Muller H: Identification of target genes of the p16INK4A-pRB-E2F pathway. J Biol Chem. 2003, 278: 46124-46137. 10.1074/jbc.M304930200.

    Article  CAS  PubMed  Google Scholar 

  36. Markey M, Siddiqui H, Knudsen ES: Geminin is targeted for repression by the retinoblastoma tumor suppressor pathway through intragenic E2F sites. J Biol Chem. 2004, 279: 29255-29262. 10.1074/jbc.M313482200.

    Article  CAS  PubMed  Google Scholar 

  37. Black EP, Huang E, Dressman H, Rempel R, Laakso N, Asa SL, Ishida S, West M, Nevins JR: Distinct gene expression phenotypes of cells lacking Rb and Rb family members. Cancer Res. 2003, 63: 3716-3723.

    CAS  PubMed  Google Scholar 

  38. Zacharatos P, Kotsinas A, Evangelou K, Karakaidos P, Vassiliou LV, Rezaei N, Kyroudi A, Kittas C, Patsouris E, Papavassiliou AG, Gorgoulis VG: Distinct expression patterns of the transcription factor E2F-1 in relation to tumour growth parameters in common human carcinomas. J Pathol. 2004, 203: 744-753. 10.1002/path.1582.

    Article  CAS  PubMed  Google Scholar 

  39. Paris PL, Andaya A, Fridlyand J, Jain AN, Weinberg V, Kowbel D, Brebner JH, Simko J, Watson JE, Volik S, Albertson DG, Pinkel D, Alers JC, van der Kwast TH, Vissers KJ, Schroder FH, Wildhagen MF, Febbo PG, Chinnaiyan AM, Pienta KJ, Carroll PR, Rubin MA, Collins C, van Dekken H: Whole genome scanning identifies genotypes associated with recurrence and metastasis in prostate tumors. Hum Mol Genet. 2004, 13: 1303-1313. 10.1093/hmg/ddh155.

    Article  CAS  PubMed  Google Scholar 

  40. Pickering MT, Kowalik TF: Rb inactivation leads to E2F1-mediated DNA double-strand break accumulation. Oncogene. 2005

    Google Scholar 

  41. Michel LS, Liberal V, Chatterjee A, Kirchwegger R, Pasche B, Gerald W, Dobles M, Sorger PK, Murty VV, Benezra R: MAD2 haplo-insufficiency causes premature anaphase and chromosome instability in mammalian cells. Nature. 2001, 409: 355-359. 10.1038/35053094.

    Article  CAS  PubMed  Google Scholar 

  42. Wang X, Zou L, Zheng H, Wei Q, Elledge SJ, Li L: Genomic instability and endoreduplication triggered by RAD17 deletion. Genes Dev. 2003, 17: 965-970. 10.1101/gad.1065103.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Ouspenski II, Elledge SJ, Brinkley BR: New yeast genes important for chromosome integrity and segregation identified by dosage effects on genome stability. Nucleic Acids Res. 1999, 27: 3001-3008. 10.1093/nar/27.15.3001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Cortez D, Guntuku S, Qin J, Elledge SJ: ATR and ATRIP: partners in checkpoint signaling. Science. 2001, 294: 1713-1716. 10.1126/science.1065521.

    Article  CAS  PubMed  Google Scholar 

  45. Smith L, Liu SJ, Goodrich L, Jacobson D, Degnin C, Bentley N, Carr A, Flaggs G, Keegan K, Hoekstra M, Thayer MJ: Duplication of ATR inhibits MyoD, induces aneuploidy and eliminates radiation-induced G1 arrest. Nat Genet. 1998, 19: 39-46.

    Article  CAS  PubMed  Google Scholar 

  46. Habedanck R, Stierhof YD, Wilkinson CJ, Nigg EA: The Polo kinase Plk4 functions in centriole duplication. Nat Cell Biol. 2005, 7: 1140-1146. 10.1038/ncb1335.

    Article  CAS  PubMed  Google Scholar 

  47. Bischoff JR, Plowman GD: The Aurora/Ipl1p kinase family: regulators of chromosome segregation and cytokinesis. Trends Cell Biol. 1999, 9: 454-459. 10.1016/S0962-8924(99)01658-X.

    Article  CAS  PubMed  Google Scholar 

  48. Olivier M, Eeles R, Hollstein M, Khan MA, Harris CC, Hainaut P: The IARC TP53 database: new online mutation analysis and recommendations to users. Hum Mutat. 2002, 19: 607-614. 10.1002/humu.10081.

    Article  CAS  PubMed  Google Scholar 

Pre-publication history

Download references

Acknowledgements

We thank members of the UCSF Comprehensive Cancer Center Genome Analysis Shared Resource Facility, Sonia Mirza, Julie Weng, Maimie Yu, and Facility Manager, David Ginzinger for carrying out the TP53 sequencing. This work was supported by NIH grants CA90421 and CA101359 (JF, AMS, BY, HL, DGA), CA78913 (RS, TT, ANJ, DP), and CA58207, the Office of Health and Environmental Research of the U.S. Department of Energy (Contract DE-AC03-76SF00098) and the Avon foundation (SD, BML, JM, JZ, KC, SD, HF, JWG, FW).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Donna G Albertson.

Additional information

Competing interests

The author(s) declare that they have no competing interests.

Authors' contributions

JF performed the statistical analysis and contributed to writing the manuscript, AMS and BY performed the array CGH and contributed to writing the manuscript, HL performed the telomere length measurements, AO performed statistical analysis, RS performed array CGH, SD directed the dissection of tumor tissue and extraction of nucleic acids, TT contributed to the development of analytical procedures for image analysis and performed the image analysis, BML was responsible for pathologic diagnosis, AJN contributed to the development of analytical procedures, JM and JZ manage and direct the UCSF Comprehensive Cancer Center Risk Program and provided access to BRCA1 mutant tumors, KC and SD performed the array CGH for breast study II, HF was responsible for the expression array measurements, JG and FW directed breast study II and shared their unpublished data, DP participated in technology development that made this study possible, DGA directed the study, participated in technology development and was responsible for preparation of the manuscript.

Jane Fridlyand, Antoine M Snijders, Bauke Ylstra contributed equally to this work.

Electronic supplementary material

12885_2006_448_MOESM1_ESM.xls

Additional File 1: Additional file 1 – Supplementary Table 1. Patient CharacteristicsProvides information on patient samples. (XLS 38 KB)

12885_2006_448_MOESM2_ESM.doc

Additional File 2: Additional file 2 – Supplementary Statistical Methods Provides a detailed description of statistical methods. (DOC 54 KB)

12885_2006_448_MOESM3_ESM.xls

Additional File 3: Additional file 3 – Supplementary Table 2. Array CGH data Provides array CGH log2ratios for clones and genome order for all samples. (XLS 3 MB)

12885_2006_448_MOESM4_ESM.xls

Additional File 4: Additional file 4 – Supplementary Table 3. Genome stability genes and associations with copy number aberrations Provides information on 426 genome stability genes and associations with copy number aberration types. (XLS 662 KB)

12885_2006_448_MOESM5_ESM.xls

Additional File 5: Additional file 5 – Supplementary Table 4. Genome stability genes and their relation to E2F. Table provides information on whether the genes are known targets of E2F1, correlation of their expression with E2F1, correlation of their expression with copy number in samples in which copy number of the locus is reduced. (XLS 394 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Fridlyand, J., Snijders, A.M., Ylstra, B. et al. Breast tumor copy number aberration phenotypes and genomic instability. BMC Cancer 6, 96 (2006). https://doi.org/10.1186/1471-2407-6-96

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2407-6-96

Keywords