| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Molecular Biology, Pathobiology and Genetics |
1 Department of Oncology, University Hospital; 2 Lund Strategic Research Center for Stem Cell Biology and Cell Therapy and 3 Department of Theoretical Physics, Lund University; and 4 Department of Pathology, University Hospital, Lund, Sweden; 5 Abramson Family Cancer Research Institute, University of Pennsylvania, Philadelphia, Pennsylvania; and 6 Section of Cancer Genetics, Institute of Cancer Research, Sutton, Surrey, United Kingdom
Requests for reprints: Åke Borg, Department of Oncology, Lund University, SE-221 85 Lund, Sweden, SE 22185. Phone: 46-46-177502; Fax: 46-46-147327; E-mail: ake.borg{at}onk.lu.se.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Earlier studies have suggested different somatic genetic pathways in progression of inherited and sporadic tumors, and several histopathologic and clinical features differ among BRCA1, BRCA2, and sporadic breast cancers. For instance, BRCA1 tumors are often estrogen and progesterone receptor (ER/PR) negative, whereas BRCA2 tumors predominantly are ER and PR positive (47). BRCA1 tumors are primarily of high histologic grade and manifest high lymphocyte infiltration and continuous "pushing" tumor margins. BRCA2 tumors are more heterogeneous but have significantly less tubule formation compared with sporadic tumors, which account for an overall higher histologic grade (7). Using global gene expression profiling, we have previously identified a set of genes to distinguish among BRCA1, BRCA2, and sporadic tumors (8). In a parallel fashion, we used conventional metaphase comparative genomic hybridization (CGH) to characterize chromosomal aberrations in hereditary breast tumors (911). We found a higher frequency of copy number alterations among BRCA1 and BRCA2 tumors when compared with sporadic cases. BRCA1 tumors harbor frequent loss at 2q, 4p, 4q, 5q, and 12q, whereas BRCA2 tumors are characterized by a higher frequency of 6q and 13q losses as well as gains on 17q22-24 and 20q13. Others have also used metaphase CGH analysis to classify BRCA1 tumors based on patterns of genomic alterations and suggested gain of 3q and loss of 3p and 5q to distinguish BRCA1 from sporadic tumors (12).
Taken together, these findings imply that tumor profiling can be useful as diagnostics tools and encouraged us to extend our previous studies with an array-based CGH technique, providing higher consistency and resolution, to build classifiers based on genomic aberration patterns in tumor tissue. In the familial cancer clinic, candidates for BRCA1 and BRCA2 mutation screening are selected based on family history, age at onset, bilateral disease, and cases of ovarian cancers in the family. Because screening of these genes is laborious and costly, more precise criteria for selection and diagnosis is needed (e.g., a rationale corresponding to the tests for microsatellite instability and mismatch repair protein expression in diagnosis of hereditary nonpolyposis colon cancer; ref. 13). In the present study, we established DNA copy number profiles for BRCA1, BRCA2, and sporadic breast tumors using a microarray containing
5,000 individual BAC clones providing an average resolution of 1 Mbp. Furthermore, by comparing DNA copy number and gene expression in 11 BRCA1/2 tumors, we found novel genes affected by genomic aberrations in hereditary breast cancer. Our results form the basis for a new classification system of inherited breast cancer and reveal new molecular targets in tumor progression of inherited tumors. The approach of using array CGH in diagnostics may be more advantageous than gene expression profiling because paraffin-embedded material can be used and thus more feasible in the clinical setting.
| Materials and Methods |
|---|
|
|
|---|
|
3,600 fluorescence in situ hybridization and sequence-tagged site (STS)mapped BAC clones obtained from the RPCI-11 and CalTech A and B libraries, and
1,400 from the end-sequenced collection at The Institute for Genomic Research. Genome coverage was assessed by the UCSC human genome assembly Build 34 (http://www.genome.ucsc.edu), revealing an average resolution of
1 Mbp, with no gaps >2 Mbp. BAC DNA was amplified using degenerate oligonucleotide PCR (DOP) primers. Arrays were constructed with at least two replicates per clone on each slide using a Molecular Dynamics Gen3 spotter and a spotting solution of 50% DMSO (14). Subsequently, microarrays with tiling coverage and high resolution of chromosome 5 (>2,600 clones) were produced from BAC clones included in the 32K set of CHORI BACPAC Resources (http://bacpac.chori.org/genomicRearrays.php). Here, DOP-PCR products were obtained from 6 ng BAC DNA template and purified using filter based 96-well plates (PALL), dried, and resuspended in 50% DMSO to a concentration of 500 to 1,000 ng/µL. Arrays were printed using a MicroGrid II spotter (Biorobotics, Cambridge, MA) as described in detail elsewhere.7 DNA extraction and array comparative genomic hybridization hybridization. Genomic DNA was extracted from frozen tissue sections using a proteinase K treatment followed by phenol chloroform purification. DNA extraction from paraffin-embedded tumors was done using a xylene treatment before proteinase K digestion and phenol chloroform purification. For all samples, 1 µg of genomic DNA was labeled according to published protocols (14) using a random labeling kit (Invitrogen Life Technologies, Carlsbad, CA). Tumor DNA and male reference DNA was differentially labeled, pooled, mixed with human COT-1 DNA, dried down, and resuspended in a formamide-based buffer (14). The hybridization reactions were applied to arrays and the arrays were incubated under coverslips for 48 to 72 hours at 37°C. Slides were washed according to published protocols (15) and scanned using an Axon 4000A scanner (Axon Instruments, Wheatherford, TX). All experiments were done with "dye-swap" to account for variations in dye labeling efficiency and fluorescence emission and all BACs are spotted on the array in duplicate; thus, each data point is represented by four replicates.
RNA extraction and gene expression analysis. RNA was extracted from freshly frozen tumor tissue using Trizol reagent (Invitrogen Life Technologies) followed by RNeasy Midi kit (Qiagen, Chatsworth, CA). RNA quality was assessed using the RNA 6000 Nano LabChip Kit for Agilent 2100 bioanalyzer (Agilent Technologies, Palo Alto, CA) and concentration was determined using a NanoDrop Spectrophotometer (NanoDrop Technologies, Wilmington, DE). Arrays were produced by the SWEGENE DNA Microarray Resource Centre, Department of Oncology at Lund University. Human Genome Oligo Set Version 2.1 (containing 21,329 70-mer probes) and Human Genome Oligo Set Version 2.1 Upgrade (containing 5,462 70-mer probes) were obtained from Operon Biotechnologies, Inc. (Huntsville, AL). cDNA synthesis, labeling, and cleanup were overall done according to manufacturers' instructions using Pronto! Plus System 6 (Corning, Inc., Corning, NY) with 5 µg of total RNA as starting material. As reference RNA a commercial pool of cell lines (Stratagene Universal Reference, La Jolla, CA) was used. Arrays were scanned at two wavelengths using an Agilent G2505A DNA microarray scanner (Agilent Technologies).
Data analysis. Identification of individual spots on scanned arrays was done with Gene Pix Pro 4.0 (Axon Instruments, Union City, CA), and the quantified data matrix was loaded into Bio Array Software Environment BASE (16). Background-correction of Cy3 and Cy5 intensities was calculated using the median feature and median local background intensities provided in the quantified data matrix. Within arrays, intensity ratios for individual probes were calculated as background corrected intensity of tumor sample divided by background corrected intensity of reference sample. For the gene expression data spots for which background-corrected Cy3 or Cy5 intensities were <1 or >65,000 were removed from further analysis. In addition, spots that had been flagged during image analysis, either by the GenePix software or after manual inspection, were removed. Data within individual arrays were then adjusted by an implementation of the intensity-dependent normalization method based on a lowess fit, as described by Yang et al. (17). A filter was applied on remaining spots excluding all spots having signal-to-noise ratio of <2 and intensity less than background intensity. Next, data were mean centered using an iterative process so that all row-wise and column-wise mean log2 ratio values are close to 0.
For the CGH data, a signal-to-noise filter of
2 for the tumor channel and
1.5 for the reference channel was applied to the data and spots that failed to pass these criterions were excluded from further analysis and regarded as missing values. Average intensity ratios were calculated for spots present in both dye-swap hybridizations after filters had been applied and used in subsequent analysis. Furthermore, BAC clones with >10 missing values across the 52 tumors were excluded from further analysis (81% presence required). The filtered data was, for each array separately, centralized to a median ratio of unity. All filtering and normalization was done in BASE (16). The CGH-Plotter software, as an R (http://www.r-project.org) implementation in BASE, was used to identify regions of gains and losses (18). In CGH-Plotter, a moving median sliding window of three clones and a constant variable value of 15 were applied. Cutoff ratios for gains and losses were set to 1.15 and 0.87, respectively, corresponding to log2 (ratio) of ±0.2. All clones were designated gained, lost, or not changed giving us a ternary scale. Using the values given from CGH-Plotter, we calculated the percentage of altered clones, the number of altered regions in each tumor, and the mean value within each group. We used a two-sided Mann-Whitney test to calculate significant differences in copy number frequency between groups. To identify clones discriminating the three tumor groups, we did a Kruskal-Wallis test (19) as implemented in the Statistics Toolbox in Matlab (Mathworks, Natick, MA), on the ternary data. To determine the correlation coefficients, we used a Pearson correlation estimate using the results obtained from the CGH-Plotter software. For the Kruskal-Wallis test, only primary tumors were included in the analysis. Hierarchical clustering was done in BASE using a bottom-up approach with center of mass linkage and Pearson correlation coefficient distance measure.
Correlation of gene expression and copy number data was determined as follows. After normalization and quality filters had been applied, we selected oligonucleotide probes (henceforth referred to as reporter) for which expression data was present in 80% of 11 tumor samples. The data set was then filtered for variance across all hybridizations for each reporter individually, removing reporters having a log2 ratio SD of <0.5. To map reporters to genomic location of corresponding transcript target, each reporter sequence was blasted against Genbank sequences in Unigene build 176. We used 100% sequence similarity as cutoff criteria for associating a reporter with accession sequences. The ACID database (http://bioinfo.thep.lu.se/acid.html; ref. 20) was used to extract full gene information for accession numbers associated to each reporter through the blast analysis (Unigene build 176). Each gene was then mapped to a chromosomal position based on its gene symbol trough the UCSC Genome Browser (Hg16). For each BAC clone present in the BAC data set, we extended its stop position to cover the full base pair range between its adjacent BAC neighbors. The new BAC mapping was then connected to the previously mapped oligonucleotide reporters, creating reporter-BAC pairs. Standard Pearson correlation (henceforth referred to as correlation) for each reporter-BAC pair, in total 7,410, was then calculated to find pairs where gene expression and genomic copy number behaved concordantly. Next, sample assignments for the copy number profiles were randomly permuted and correlations recalculated keeping the expression data and reporter-BAC pair mapping intact. This procedure was repeated 10,000 times. The permutations allowed us to calculate P values for the different correlation bins (each bin represent a correlation range of 0.05). A P value cutoff of 0.01 corresponded to a correlation cut off of 0.75; that is, by chance, one would expect 1% of all pairs to have a correlation of
0.75. To determine whether the gene is "overexpressed" or "underexpressed," we used an SD method based on using a sliding window through an M-A plot. The method uses a defined window size for calculation of mean and SD for spots within the window. Each spot is assigned a window and is then compared with mean ratio and SD to see how many SD the spot is from the mean (21).
For classification of BRCA and sporadic breast tumors, we used supervised support vector machines in a leave-one-out iterative validation approach included in the TIGR MultiExperimentViewer (TMEV) software (22, 23). In the classification, whole genome copy number profiles obtained from the array CGH experiment were used. In the leave-one-out procedure, one tumor is left out and the remaining samples are used to construct a classifier that is validated on the left-out sample. The procedure is iterated for each sample, such that every sample gets a classification. A linear kernel with default values for the variables, except the diagonal factor set to 4, was used within the TMEV software (22). The significance of classification results was calculated using Fisher's exact test on 2 x 2 contingency tables using R (http://www.r-project.org).
| Results |
|---|
|
|
|---|
5,000 individual BAC clones evenly distributed over the genome, corresponding to an average resolution of
1 Mbp. Gains, losses, and high level amplifications were readily detected as shown by visualization in overview plots (Fig. 1A-B). The sizes of regions manifesting gain or loss were determined using the distance between the two closest, not altered clones and ranged from 1 Mbp to whole chromosome gains/losses. Overall, we found a higher frequency of copy number changes in BRCA1 (47.8 ± 12.7% altered clones) than in BRCA2 (34.9 ± 10.5% altered clones, P = 0.029) and sporadic tumors (28.1 ± 14.7% altered clones, P = 0.00078). The difference between BRCA2 and sporadic tumors was not significant (P = 0.18). In addition, we found a higher number of distinct altered regions in BRCA1 (91.1 ± 20.8) than in BRCA2 (79.4 ± 45.7, P = 0.064) and sporadic tumors (72.6 ± 32.1, P = 0.091). There was no apparent association between type or location of the BRCA1 or BRCA2 mutations and copy number changes, although it should be pointed out that all mutations but one were of protein truncating type and clustered within the central regions of the genes. There was, however, a concordance in genomic profiles of tumors occurring in the same patient (see below). The two male BRCA2 tumors had an overall genomic profile similar to BRCA2 tumors from female patients (see below).
|
3.5 Mbp and including >50 genes. Another region with frequent (83%) gain in BRCA1 tumors was found at 7q36.1-q36.3, including the cyclin-dependent kinase 5 (CDK5) gene. Narrow amplifications were detected on 6q, including the MYB oncogene and HBS1L gene. In BRCA2 tumors, overlapping regions with frequent (>75%) gains included 1q32.1-q41, 8q21.13, 8q22.1-q24.3, 17q23.3-q24.2, 17q25.1-qter, and 20q12-q13.12. The 17q23.3-q24.2 and 20q13.13 amplicons were shown to be most specific for BRCA2 tumors. A narrow amplification peak on 12q14.2-q21.1 including the MDM2 oncogene was revealed in one BRCA2 tumor. In addition, corroborating earlier findings, BRCA2 tumors exhibited recurrent gains on 20q, including high frequency of gain on 20q13.13. Among other frequently altered genes in breast cancer, the ERBB2 gene on 17q12 was amplified in seven (27%) of sporadic tumors but in none of the BRCA1 or BRCA2 tumors. The CCND1 gene on 11q13.3 was amplified in only one of the sporadic tumors, whereas another sporadic tumor (TB088) manifested discrete amplification of the EMSY gene, located 6 Mbp distal to CCND1. Regions with deletions. Regions showing recurrent loss in BRCA1 tumors (>75% of samples) included 4p15.32-p14, 4q31.3, 4q32.1-q34, 5q11.2-q23.3, 8pter-p12, 13q13.3-q21.32, 15q12-13.1, 15q15.3-q21.1, and 17p13.2-p12. Some of these regions (4q, 8p, and 13q) were commonly deleted (35-61%) also in sporadic tumors, although not to the same extent as in the inherited breast cancers. In addition to frequent heterozygous losses on chromosome 5q, a putative homozygous deletion on 5q12.1 was detected in a primary tumor and a local recurrence obtained from a BRCA1 mutation carrier. This deletion comprised two adjacent BAC clones present on the 1 Mb array resulting in log2 (ratios) below 1.5. The size of the deleted region was estimated to 2.8 Mbp and included nine known genes. Analysis of these samples using a BAC array with complete coverage and high-resolution (in average 50 kbp) of chromosome 5, narrowed down this homozygous deletion to a 700-kbp region including only three genes (PART, FKSG52, and DEPDC1B) and the 3' flanking regions of PDE4D (Fig. 2A). Primers for STSs mapped to the putative homozygous deletion were used for PCR analysis, which confirmed the homozygote loss in both cases (Fig. 2B). In 72% of the BRCA1 tumors investigated, a heterozygous loss of 5q12.1 was present, suggesting that this region harbors genes of importance for BRCA1 tumor development. Regions showing recurrent loss (>75%) in BRCA2 tumors included 8p23.3-p21.2, 11q14.3-q21, 11q24.2-q25, 13q13.3-q14.13, 16q22.2-23.2, and 17p13.1-p12. Of these regions, losses on 11q were the most specific alterations for BRCA2 tumors. The region comprising BRCA1 and adjacent clones on 17q21.31 was deleted in six BRCA1, three BRCA2, and two sporadic tumors, whereas the BRCA2 gene on 13q13.1 was deleted in nine BRCA1, eight BRCA2, and 12 sporadic tumors.
|
Correlation of gene expression and DNA copy number. To determine the effect of DNA alterations on gene expression, we did genome wide expression analysis on four BRCA2- and seven BRCA1-associated breast tumors from which array CGH results had been obtained. By matching BAC and oligonucleotide probes (in total, 7,410 BAC-oligonucleotide probe pairs), a significant correlation between genomic alterations and gene expression was found for 746 genes (P < 0.01), including several relevant candidate genes such as RAD50, HBS1L, BCAS4, CDK5, and TM4SF1. Moreover, 30 of the 746 genes were mapped to the chromosomal regions (see below) discriminating BRCA1, BRCA2, and sporadic tumors (Fig. 3).
|
5,000 BAC clones on our arrays, 3,591 passed our missing value requirements in these experiments. First, we used CGH-Plotter to assign a discrete state (gained, lost, or unchanged) to each copy number estimate for these 3,591 clones. Second, unsupervised hierarchical clustering using 3,591 BAC clones was conducted where a trend of subgrouping was apparent (Fig. 4A). Next, we used a Kruskal-Wallis test and supervised model to identify BAC clones that discriminate among the three tumor groups, including only primary tumors in the analysis. We found 169 significant BAC clones (P < 0.001), whereas only four clones would be expected by chance alone. Several regions confirmed findings from other studies using techniques with lower resolution, including large regions on chromosomes 4p, 4q, 5q, and 12q we previously showed as being frequently deleted in BRCA1 tumors. In addition, deletions on chromosome 15q were distinct features of BRCA1 tumors as well as gains on chromosome 3q and 7q. Moreover, gains on 7p and 17q24 and losses on 8p were found to be more prevalent in BRCA2 cases. Using hierarchical clustering, we could discriminate BRCA1, BRCA2, and sporadic tumors based on the 169 BAC clones (Fig. 4B). In the cluster analysis, 14 of 14 BRCA1 tumors were tightly clustered and separated from sporadic cases. BRCA2 tumors showed a somewhat higher similarity with the sporadic breast tumors but still displayed a distinct genomic profile of their own. Eight of 12 BRCA2 tumors were grouped in a subcluster whereas one outlier case clustered together with the sporadic cases and three showed high similarity to the BRCA1 tumors. In addition, two sporadic tumors (TB088 and TB026) showed high similarity to BRCA1 tumors.
|
|
| Discussion |
|---|
|
|
|---|
Here, we report genomic profiles for BRCA1 and BRCA2 tumors using DNA arrays comprised of
5,000 individual BAC clones providing an average resolution of 0.9 Mb. Genomic gains, heterozygous losses, high-level amplifications, and homozygous deletions were mapped. Our results confirm that BRCA1 and BRCA2 tumors exhibit a higher frequency of copy number alterations, although mainly shown here for the BRCA1 tumors compared with sporadic breast tumors with regard to percentage of altered clones. Although the number of discrete chromosomal regions altered was higher in hereditary than in sporadic tumors, this difference was not significant, indicating the presence of large alterations which signifies the gross chromosomal instability and aggressive phenotype generally displayed by BRCA1 and BRCA2 tumors. Both gene products are, among a variety of functions, important key players in homologous recombination mediated DNA repair; thus, loss of function increase genomic instability (3). Moreover, these tumors may evolve along distinct genomic programs as suggested by the occurrence of specific recurrent chromosomal alterations. For instance, the majority of all BRCA1 tumors showed losses on chromosomes 4p, 4q, and 5q, often accompanied by copy number gains or amplifications on 3q23-qter. The PIK3CA oncogene is one of several candidate genes in the relatively large 3q region. Interestingly, a potential role of CDK5 in the 7q36.1 amplicon is corroborated by gene expression analysis, where one of the BRCA1 tumors with 7q36-qter amplification was shown to have an elevated expression of CDK5 in comparison with a large number of other tumors (data not shown).
To further elucidate the effect of copy number changes to gene expression, we did gene expression analysis on 11 BRCA1 and BRCA2 tumors and found 746 genes for which DNA copy number and RNA expression were significantly correlated (P < 0.01). Interestingly, TM4SF1 located on 3q25.1, a region known to be commonly amplified in BRCA1 tumors, showed a very high correlation with gene expression data. TM4SF1 is a member of the transmembrane 4 superfamily, also known as the tetraspanin family. These are putative regulators of cancer cell motility and metastasis (26) which indicates that TM4SF1 is a candidate oncogene for BRCA1 tumors. Furthermore, several regions on 5q are frequently deleted in BRCA1 tumors, including 5q23.3 and RAD50, another gene for which gene and transcript copy numbers were significant correlated. RAD50 interacts with BRCA1 and is part of the "BRCA1-associated genome surveillance complex" that serves as a sensor of abnormal DNA structures (27). Moreover, BRCA2 tumors have a high frequency of DNA copy number gains on 20q13, which was confirmed in the present study suggesting a role of BCAS4 as one of the target genes. Recently, BCAS4 was found to be amplified and overexpressed in 9 of 13 breast cancer cell lines, suggesting an oncogenic function in sporadic breast cancer as well (28).
Intriguingly, we found a striking concordance regarding genomic profiles for two different tumors obtained from one patient carrying a BRCA2 mutation. These tumors were regarded as independent primary cancers as determined by histopathologic variables and arose in opposite breasts with a considerable time interval. Both tumors were invasive and included extensive DCIS, further arguing for an independent origin and against a metastatic spread. Recently, a study using metaphase CGH analysis of bilateral sporadic breast cancer revealed widely different alterations in the contralateral tumors, indicating that most bilateral cancers indeed have different origin (29). Although we cannot exclude the possibility that a contralateral breast cancer arise via hematogenic metastasis, our results indicate that a predisposing BRCA gene mutation and other constitutional factors, possibly in combination with a specific progenitor cell type of origin, may influence tumor development and genomic alterations along distinct genetic pathways or programs.
Heterozygous 5q deletions have previously been found to be common in BRCA1 tumors indicating the presence of one or several tumor suppressor genes (9). Moreover, CGH-targeted linkage analysis revealed a possible BRCA1 modifier locus on chromosome on 5q33-q34 (30). In our set of 12 independent primary tumors from BRCA1 mutation carriers, we found nine heterozygous and one homozygous deletion on 5q12.1. This novel homozygous deletion was found in both the primary tumor and its local recurrence. To determine its exact boundaries, we used a BAC array with complete coverage of chromosome 5, which allowed narrowing the region to approximately 700 kbp. The region fully comprises two known genes, PART1 that regulated by androgens in prostate cancer cells (31) and DEPDC1B that encodes a Rho GTPase activation protein and a putative tumor suppressor protein. The region is situated in the proximity of a proposed familial prostate cancer susceptibility locus (32) and also includes a single exon gene (FKSG52 or AF336879) of unknown function and the 3' flanking region of PDE4D, which encodes a phosphodiesterase associated with ischemic stroke (33). The role of PART1, FKSG52, PDE4D, or DEPDC1B in BRCA1 tumor development remains to be determined.
Using hierarchical clustering methods to group tumors based on ternary data for the 169 discriminating BAC clones, all BRCA1 tumors were tightly clustered and separated from sporadic cases. BRCA2 tumors showed a somewhat higher similarity with the sporadic tumors but still displayed a distinct genomic profile of their own. This is a remarkable finding in light of the inclusion of DNA extracts from both fresh-frozen and formalin-fixed tumors. Eight of 12 BRCA2 tumors were tightly grouped in a subcluster, whereas one BRCA2 outlier clustered together with the sporadic cases and three outlier cases clustered with BRCA1 tumors. This was not due to inclusion of male BRCA2 breast cancers, which had an overall similar genomic profile as female tumors, as shown earlier (10). It might, however, reflect that BRCA2 tumors constitute a more heterogeneous group, as also shown by a number of different histopathologic and clinical variables (7). It could also be a result of the relatively small number of samples in the present study. In addition, two sporadic cases seemed highly similar to the BRCA1 tumors with regard to their genomic profiles. Sporadic tumors with methylation of the BRCA1 gene promoter have been shown to display an expression profile similar to BRCA1-mutated tumors (8), and further analysis will reveal if BRCA1 methylation can influence also the genomic profiles of tumors. On the other hand, BRCA-like sporadic tumors could also represent missed BRCA1 germ line mutation carriers. Interestingly, one of the sporadic tumors (TB088) in the BRCA1 cluster exhibited amplification of the EMSY oncogene, which encodes a novel BRCA2-binding protein that suppress the chromatin remodeling function of BRCA2 (34). It is tempting to speculate in the role of EMSY activation to promote a tumor development similar to inherited tumors, but this will clearly need further studies to prove.
To further elucidate the biological mechanism behind the distinct copy number profiles for BRCA1/2-associated breast cancers, we extracted 30 of 746 significant genes correlating in gene expression and copy number data mapped to the discriminating regions pinpointing genes that might be important in the distinct tumor progression patterns in BRCA1 and BRCA2 breast cancers (Fig. 3). In addition to hierarchical clustering methods, a supervised leave-one-out iterative classification of the samples was conducted based on whole genome copy number data. Using this model, only one BRCA1 breast tumor was misclassified. Interestingly, this tumor was obtained from the single missense mutation (Cys64Gly) carrier included in the study, whereas all other BRCA tumors were from nonsense or frameshift mutation carriers. Although this mutation is known to be deleterious and disease-associated, transcript with Cys64Gly may escape nonsense-mediated decay, become translated, and retain some of the many functions of wild-type BRCA1 protein or acquiring new ones. In addition, the three sporadic breast tumors that were classified in the BRCA1 group were all ER and ERBB2 negative and of histologic grade 3 thus resembling tumors of basal-like phenotype that are similar to BRCA1 breast tumors with regard to gene expression profiles (35). In the BRCA2 classification model, 9 of 12 BRCA2 tumors were correctly classified indicating that these tumors constitute a more heterogeneous group than BRCA1 tumors. In addition, four sporadic tumors were classified as BRCA2-associated tumors. Nonetheless, even if the results of the present study must await confirmation in a larger and independent tumor material, our overall findings strongly suggest that BRCA1 and BRCA2 tumors can be classified based on copy number alterations. In addition, using CGH arrays with higher resolution, such as tiling DNA microarrays with complete coverage of the human genome, the regions discriminating the different tumor groups will be refined (36). Such a BRCA tumor classifier, applicable to formalin-fixed, paraffin-embedded tumor specimens, would be a valuable tool for test laboratories and oncogenetic clinics. It may also become of value for selection of "BRCAness" tumors for targeted treatment (37, 38).
| Acknowledgments |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
We thank Kazutoyo Osoegawa and Pieter de Jong (BACPAC Resources, Children's Hospital Oakland Research Institute, Oakland, CA; http://bacpac.chori.org/genomicRearrays.php) for obtaining the BAC clone DNA aliquots from the 32K set.
| Footnotes |
|---|
Received 2/18/05. Revised 6/11/05. Accepted 6/24/05.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
L. Melchor and J. Benitez An integrative hypothesis about the origin and development of sporadic and familial breast cancer subtypes Carcinogenesis, August 1, 2008; 29(8): 1475 - 1482. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. P. M. Suijkerbuijk, M. J. Fackler, S. Sukumar, C. H. van Gils, T. van Laar, E. van der Wall, M. Vooijs, and P. J. van Diest Methylation is less abundant in BRCA1-associated compared with sporadic breast cancer Ann. Onc., July 22, 2008; (2008) mdn409v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. V. Iakovlev, N. C.R. Arneson, V. Wong, C. Wang, S. Leung, G. Iakovleva, K. Warren, M. Pintilie, and S. J. Done Genomic Differences Between Pure Ductal Carcinoma In Situ of the Breast and that Associated with Invasive Disease: a Calibrated aCGH Study Clin. Cancer Res., July 15, 2008; 14(14): 4446 - 4454. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. A. Rakha, J. S. Reis-Filho, and I. O. Ellis Basal-Like Breast Cancer: A Critical Review J. Clin. Oncol., May 20, 2008; 26(15): 2568 - 2581. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Melchor, E. Honrado, J. Huang, S. Alvarez, T. L. Naylor, M. J. Garcia, A. Osorio, D. Blesa, M. R. Stratton, B. L. Weber, et al. Estrogen Receptor Status Could Modulate the Genomic Pattern in Familial and Sporadic Breast Cancer Clin. Cancer Res., December 15, 2007; 13(24): 7305 - 7313. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. K. Bodvarsdottir, L. Vidarsdottir, and J. E. Eyfjord AURKA and Breast Cancer in BRCA1/2 Mutation Carriers Cancer Epidemiol. Biomarkers Prev., November 1, 2007; 16(11): 2517 - 2517. [Full Text] [PDF] |
||||
![]() |
J. Huang, A. Gusnanto, K. O'Sullivan, J. Staaf, A. Borg, and Y. Pawitan Robust smooth segmentation approach for array CGH data analysis Bioinformatics, September 15, 2007; 23(18): 2463 - 2469. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Vincent-Salomon, C. Ganem-Elbaz, E. Manie, V. Raynal, X. Sastre-Garau, D. Stoppa-Lyonnet, M.-H. Stern, and E. Heard X Inactive-Specific Transcript RNA Coating and Genetic Instability of the X Chromosome in BRCA1 Breast Tumors Cancer Res., June 1, 2007; 67(11): 5134 - 5140. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Y. Hehir-Kwa, M. Egmont-Petersen, I. M. Janssen, D. Smeets, A. G. van Kessel, and J. A. Veltman Genome-wide Copy Number Profiling on High-density Bacterial Artificial Chromosomes, Single-nucleotide Polymorphisms, and Oligonucleotide Microarrays: A Platform Comparison based on Statistical Power Analysis DNA Res, March 15, 2007; (2007) dsm002v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Harris, M. Lowden, I. Clejan, M. Tzoneva, J. H. Thomas, J. Hodgkin, and S. Ahmed Mutator Phenotype of Caenorhabditis elegans DNA Damage Checkpoint Mutants Genetics, October 1, 2006; 174(2): 601 - 616. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||