Integrated Profiling of Basal and Luminal Breast Cancers

Basal and luminal are two molecular subtypes of breast cancer with opposite histoclinical features. We report a combined, high-resolution analysis of genome copy number and gene expression in primary basal and luminal breast cancers. First, we identified and compared genomic alterations in 45 basal and 48 luminal tumors by using 244K oligonucleotide array comparative genomic hybridization (aCGH). We found various genome gains and losses and rare high-level gene amplifications that may provide therapeutic targets. We show that gain of 10p is a new alteration in basal breast cancer and that a subregion of the 8p12 amplification is specific of luminal tumors. Rare high-level amplifications contained BCL2L2, CCNE, EGFR, FGFR2, IGF1R, NOTCH2 , and PIK3CA . Potential gene breaks involved ETV6 and FLT3 . Second, we analyzed both aCGH and gene expression profiles for 42 basal and 32 luminal breast cancers. The results support the existence of specific oncogenic pathways in basal and luminal breast cancers, involving several potential oncogenes and tumor suppressor genes (TSG). In basal tumors, 73 candidate oncogenes were identified in chromosome regions 1q21-23, 10p14, and 12p13 and 28 candidate TSG in regions 4q32-34 and 5q11-23. In luminal breast cancers, 33 potential oncogenes were identified in 1q21-23, 8p12-q21, 11q13, and 16p1213 and 61 candidate TSG in 16q12-13, 16q22-24, and 17p13. HORMAD1 (P = 6.5 10 ) and ZNF703 (P = 7 10 ) were the most significant basal and luminal potential oncogenes, respectively. Finally, among 10p candidate oncogenes associated with basal subtype, we validated CDC123/C10orf7 protein as a basal marker. [Cancer Res 2007;67(24):11565–75]


Introduction
Breast cancer is a heterogeneous disease whose evolution is difficult to predict.Consequently, treatment is not as adapted as it should be.Breast cancer can be defined at the clinical, histologic, cellular, and molecular levels.Progress and efforts to integrate these definitions allow a better management of the disease (1).Five major breast cancer molecular subtypes (luminal A and B, basal, ERBB2 overexpressing, and normal-like) have been identified (2)(3)(4)(5).They are associated with different prognosis.Luminal A breast cancers express hormone receptors, have an overall good prognosis, and can be treated by hormone therapy.ERBB2-overexpressing breast cancers have a poor prognosis and are treated by targeted therapy using trastuzumab or lapatinib (see ref. 6 for review).No specific therapy is available for the other subtypes, although the prognosis of basal and luminal B tumors is poor.
Progress can be made in several directions.First, by identifying, among good prognosis tumors such as luminal A breast cancers, the cases that will relapse and metastasize.Second, by developing new drugs that will allow a better management of poor-prognosis breast cancers, such as basal tumors.The identification of new markers and therapeutic targets is necessary to meet these objectives.This relies on the study of tumor samples at the genomic (7)(8)(9)(10)(11)(12)(13)(14), gene expression (2,4,15,16), or proteomic (17)(18)(19) levels.
We have established here the genomic and gene expression profiles of a series of breast cancer samples by using highresolution array comparative genomic hybridization (aCGH) and DNA microarrays, respectively.We have focused our analysis on two breast cancer subtypes with overall opposite features: luminal and basal.

Materials and Methods
Breast cancer samples.Tumor tissues were collected from 202 patients with invasive adenocarcinoma who underwent initial surgery at the Institut Paoli-Calmettes (Marseilles, France) between 1992 and 2004 ( from a cohort of 1,185 patients with frozen tumor sample).Each patient gave written informed consent and the study was approved by our institutional review board.Samples were macrodissected and frozen in liquid nitrogen within 30 min of removal.Several sets of tumors were studied.Gene expression was analyzed on a set of 202 tumors.Gene expression and genome profilings were established on a subset of 74 tumors (42 basal and 32 luminal A; as defined in ref. 4) for which RNA and DNA samples were both available.Genome profiles were established on 93 tumors, including these 74 tumors and an additional set of 19 cases selected from immunohistochemical data on estrogen receptor (ER), progesterone receptor (PR), ERBB2, and other basal and luminal markers (20): 2 ER À /PR À /ERBB2 À (''triple negative'') and 17 ER + /ERBB2 À (21) tumors associated with basallike and luminal-like immunohistochemical phenotypes, respectively.Although a third of the luminal tumors were not profiled on expression microarrays, immunohistochemical-based phenotypes were established for all studied tumors: 98% and 88% of the luminal samples used in this study were ER and PR positive, respectively (see Supplementary Table S1).The set studied by aCGH thus comprised 44 basal and 49 luminal tumors (mainly luminal A tumors).All tumor sections were reviewed by two pathologists (J.Jacquemier and E. Charafe-Jauffret) before analysis.All specimens contained >60% of tumor cells (as assessed before RNA extraction using frozen sections adjacent to the profiled samples).The main histoclinical characteristics of samples are listed in Supplementary Table S1.Immunohistochemical data included status for the following: ER, PR, and P53 (positivity cutoff value of 1%); ERBB2 (0-3+ score, DAKO HercepTest kit scoring guidelines, with >1+ defined as positive); and Ki67 (positivity cutoff value of 20%).Among basal tumors were 14 medullary breast cancers (MBC) defined on Ridolfi's criteria (22).MBCs displayed characteristics similar to series in literature with 100% of samples Scarff-Bloom-Richardson grade 3, ER negative, ERBB2 negative, Ki67 positive, and 49% P53 positive.After surgery, patients were treated using a multimodal approach according to standard guidelines.Eleven normal breast tissue samples pooled in four RNA samples were also profiled on gene expression microarrays.
DNA and RNA extraction.DNA and RNA were extracted from frozen samples by using guanidium isothiocynanate and cesium chloride gradient as described previously (23).RNA integrity was controlled on Agilent Bioanalyzer (Agilent Technologies).
Gene expression profiling with DNA microarrays.Expression profiles were established for 202 breast cancers and 4 normal breast samples with Affymetrix U133 Plus 2.0 human oligonucleotide microarrays following a described protocol (24).Scanning was done with Affymetrix GeneArray scanner and quantification with Affymetrix GeneChip Operating Software.Data were analyzed by the robust multichip average method in R using Bioconductor and associated packages (25).Robust multichip average did background adjustment, quantile normalization, and summarization of 11 oligonucleotides per gene.Before analysis, a filtering process removed from the data set the genes with low and poorly measured expression as defined by an expression value inferior to 100 units in all samples.All data were then log 2 transformed for display and analysis.Before combined analysis, gene expression levels for each tumor sample were centered by the average expression levels of the four normal breast samples.The up-regulation or down-regulation status for each probe set was assigned using a threshold of |1| corresponding to twice the expression found in the normal breast pool.For multiple probe sets mapping to the same gene, the probe sets with an extension ''_at'', next ''s_at'', and followed by all other extensions were preferentially kept.When several probe sets with the best extension were available, the one with the highest median value was retained.Probe sets were mapped to the human genome according to the hg17 Build 35 from the National Center for Biotechnology Information (NCBI) Ensembl database and the University of California at Santa Cruz Genome Bioinformatics database Genome Browser.
Array-comparative genomic hybridization.Genomic imbalances of 93 breast tumors were analyzed by aCGH using 244K CGH Microarrays (Hu-244A, Agilent Technologies) following a described protocol (26).Scanning was done with Agilent Autofocus Dynamic Scanner (G2565BA, Agilent Technologies).Data analysis was done as described (27) and visualized with CGH Analytics 3.4 software (Agilent Technologies).Extraction of data (log 2 ratio) was done from CGH Analytics, whereas normalized and filtered log 2 ratio was obtained from ''Feature Extraction'' software (Agilent Technologies).Data generated by probes mapped to X and Y chromosomes were eliminated.The final data set contained 225,388 unique probes covering 22,509 genes and intergenic part following the hg17 human genome mapping.Copy number changes were characterized as reported previously (10,28).aCGH data were analyzed using circular binary segmentation (29) with default variables to translate intensity measurements in regions of equal copy number as implemented in the DNA copy R/Bioconductor package (25).Thus, each probe was assigned a segment value referred to as its ''smoothed'' value.The gain and loss status for each probe was assigned using threshold of |0.5| on at least five consecutive probes across the genome.Data were submitted to the Cluster program (30) using Pearson correlation as similarity metric and centroid linkage clustering.Results were displayed using TreeView program (30).The number of copy number transitions was computed based on the initial DNA copy segmentation by counting the number of copy number transitions in the genome (31).Identification of copy number variations was based on a previous study (32).The frequency of alterations was computed for each probe locus as the proportion of samples showing an aberration therein.Alteration frequencies were evaluated by Fisher's exact test and false discovery rate (FDR) was applied to correct the multiple testing hypothesis (33).
Combined genomic and expression data analysis.For each gene, data from two corresponding probe sets (Agilent and Affymetrix) were compared.Frequencies of combined alterations, overexpression associated with amplification and the opposite, were evaluated.Data were then log 2 transformed and submitted to the Cluster program (30) using data mediancentered on genes, Pearson correlation as similarity metric, and centroid linkage clustering.Results were displayed using the TreeView program (30). 7luorescence in situ hybridization analysis.Fluorescence in situ hybridization (FISH) analysis on tissue microarrays using bacterial artificial chromosome (BAC) probes has been described previously (34).
Immunohistochemistry on breast cancer tissues.Anti-epidermal growth factor receptor (EGFR) antibody (clone 31G7, 1:20 dilution; Zymed Laboratories) and anti-CDC123 antibody (clone 4B9, 1:300 dilution; Abnova Corp.) were used for immunohistochemistry.A pepsin treatment was done on tissue sections for 30 min before immunohistochemistry with EGFR antibody.Immunohistochemistry was done on 5-mm sections using DAKO LSAB2 kit in a DAKO Autostainer (DakoCytomation).Antibodies were incubated for 1 h in citrate buffer.After staining, slides were evaluated by two pathologists.EGFR and CDC123 (also called C10orf7) showed membrane and cytoplasmic signal, respectively.Intensity of EGFR and CDC123 staining was scored from 0 (no staining) to 3 (strong and diffuse staining).The results were estimated by the percentage of positive cells (P, from 0% to 100%) and the intensity of the staining (I, from 0 to 3) and expressed by the Quick score (Q = P Â I, from 0 to 300).
Statistical analyses.Differences in clinical and pathologic features between basal and luminal cases were assessed by Fisher's exact test for qualitative variables with discrete categories and Wilcoxon test for continuous variables.

Results
Identification of molecular subtypes of breast cancer by gene expression profiling.The five molecular subtypes of breast cancer have been identified using an intrinsic set of f500 genes (2)(3)(4).We looked for these subtypes in our series of 202 breast cancers profiled with Affymetrix DNA microarrays by using the 476 probe sets common to the intrinsic 500-gene set and our 29,623 filtered probe sets.We first defined the expression centroid of each subtype for these 476 probe sets in the 122 samples from Sorlie et al. (4) by computing the median expression for each gene in the corresponding core samples.We then measured the correlation of each of our 202 samples with each centroid: 80 samples were close to the luminal A centroid, 2 to the luminal B centroid, 60 to the basal centroid, 27 to the ERBB2 centroid, and 19 to the normallike centroid.Fourteen samples displayed a correlation inferior to 0.15 with any centroid and were not attributed to any subtype with this 476-gene set.Although basal and ERBB2 subtypes have been repeatedly recognized in independent data sets, identification of the other subtypes remains difficult (35); the intrinsic gene set is not optimal to classify all samples (2)(3)(4).In all sets, some samples remain unclassified because they have a low correlation with subtype centroids (4).Other studies have attempted to improve the stratification of breast cancer patients (see ref. 35 for details).In our study, we kept the standard definition given by the initial intrinsic gene of Sorlie et al. (4) to subtype tumors in the same conditions as in recent integrated studies.Thus, a total of 142 samples were either basal or luminal in our series (data not shown).
Genome analysis of basal and luminal breast tumors.To identify genomic alterations in basal and luminal breast cancers, we used high-resolution aCGH.DNA was available for 42 basal and 32 luminal cases as defined above.We selected 19 additional cases (2 basal and 17 luminal) from immunohistochemical data.We thus analyzed gene copy number aberrations (CNA) in 93 primary breast tumors corresponding to 45 basal and 48 luminal tumors.CNAs are reported on frequency plots in Supplementary Fig. S1A to D. We distinguished low-level CNAs, which can derive from aneuploidy, from high-level CNAs, which result from focal amplification of chromosome regions, by using two different threshold values (log 2 ratio >|0.5| and |1|) in basal (Supplementary Fig. S1A and C) and luminal (Supplementary Fig. S1B and D) tumors.
Thus, basal and luminal breast cancers show both common and specific alterations.
Gene alterations differentially associated with basal and luminal breast cancer.Unsupervised hierarchical clustering (Fig. 1) of the 93 samples and the 225,388 oligonucleotide probes (excluding X and Y probes) revealed significant differences in the spectrum of alterations in basal and luminal breast cancers.A total of 30,600 probes were differently represented in basal and luminal cases (P < 0.05, Fisher's exact test, FDR corrected).Specific regional aberrations (gains and losses defined with a threshold value of log 2 ratio >|0.5|) are shown for basal (red) and luminal (blue) tumors in Supplementary Fig. S2.Genes with significant differences of CNA frequencies are listed in Supplementary Tables S2A and B and summarized in Supplementary Table S2C.
We distinguished two types of gain alterations: copy number gains and high-level amplifications.
To characterize the presence of potential gene breakages, we identified genome alterations associated with copy number transition profiles within a chromosome as described previously (31).
Copy number transitions are likely to reflect DNA strand breakages that may lead to nonreciprocal translocations.We identified transition profiles with a threshold of aberration score (log 2 ratio) >|0.5| in 1,975 genes in at least one of the 93 tumors.A total of 501 genes showed copy number transitions in at least two different tumors (Supplementary Table S3A).Results were coherent with our previous analysis of NRG1 (38), FHIT (39), and PARK2 (40), the two latter being associated with fragile sites.Eight other large genes associated with fragile sites showed transitions: CSMD1 (5%), CSMD3 (3%), NBEA (2%), NRXN3 (2%), CDH13 (2%), ACCN1 (2%), DCC (2%), and DMD (2%; ref. 41).Copy number transitions in PTEN were detected in 4% of tumors.We also observed transition in ETV6 (Supplementary Fig. S9A), in agreement with previous FISH results (42), and in FLT3, suggesting that these genes frequently targeted in leukemia could also be affected in breast cancer (Supplementary Fig. S9B).The density of probes allowed to map the breakpoints (Supplementary Fig. S9C).Statistical analysis showed that frequent copy number transitions targeted 6 and 15 genes in basal and luminal tumors, respectively (Supplementary Table S3B).Copy number variations (32) were identified in 11 of them (Supplementary Table S3B).
Correlation between gene expression and genome CNA.To identify genes whose expression levels were significantly modified in relation to DNA copy number changes, we studied 74 samples (42 basal and 32 luminal) by both aCGH and RNA expression analyses.The genomic and gene expression profiles were then compared across subtype.To define subtype-specific candidates, aCGH data from 7,694 probes with CNA frequencies different between the two subtypes were compared with the Affymetrix data of their 1,834 corresponding independent genes (Fig. 4).The expression of 197 genes was deregulated in relation to CNA with differences between basal and luminal breast cancers (Supplementary Tables S4A and S5A).Combined genomic and gene expression status of these 197 candidate genes is shown for each individual tumor in Supplementary Tables S4B and S5B.The most significant candidates are shown in Fig. 5A and B.
A total of 106 genes showed gains correlated with up-regulated mRNA expression (Supplementary Table S4A): 73 and 33 genes (corresponding to 77 and 34 probes, respectively) were associated with basal and luminal subtypes, respectively.A total of 91 genes showed losses correlated with down-regulated mRNA expression (Supplementary Table S5A): 28 and 61 genes (corresponding to 30 and 61 probes, respectively) were associated with basal and luminal subtype, respectively.These results suggest the existence of potential oncogenes and tumor suppressor gene (TSG) differentially associated with each subtype.
Conversely, 237 candidates (152 potential oncogenes and 85 potential TSG; Supplementary Table S6) exhibited correlation between CNA and mRNA expression without specificity for basal or luminal subtype.The presence of common candidate genes in the two subtypes suggests the existence of common pathways of oncogenesis.
For each individual tumor exhibiting the most frequent copy number transitions, gene CNA and expression profiles for 12 of the 21 targeted genes are shown in Supplementary Table S7 (Affymetrix data were not available for the other 9 targeted genes).No basal tumor displayed multiple transition profiles, whereas luminal tumors showed 4 to 13 potential gene breakages.The effect of these alterations in the luminal tumors was not significant on the 5-year metastasis-free survival (MFS) and other clinical features (data not shown).Among the 12 genes exhibiting copy number transitions, none presented a down-regulated expression.Five genes, CHRD, TNXB, SPTBN5, RAB11FIP3, and BC012355, were up-regulated in 38%, 38%, 50%, 40%, and 72% of the corresponding tumors, respectively.
Data comparison with recent integrated genomic studies.The 197 basal-and luminal-specific candidate genes and the 237 common candidate genes were compared with those obtained in two recent integrated whole-genome analyses of breast cancers (Supplementary Table S8; refs.10,11).A total of 1,183 candidate genes (11) and 66 highly amplified and overexpressed genes (10) were compared with the here-defined candidates.
In conclusion, our study, using a higher-resolution approach, generated a list of 73 and 33 candidate oncogenes exclusively associated with basal or luminal tumors, respectively (see Fig. 5 and Supplementary Tables S4A and C).Among them, 65 and 25 genes were newly defined as candidate oncogenes associated with basal and luminal subtypes, respectively.HORMAD1 (P = 6.5 Â 10 À5 ) and ZNF703 (P = 7 Â 10 À4 ) were the most significant basal and luminal potential oncogenes, respectively.
Genome profiling of MBC.Gene expression profiling has defined MBCs as a subset of basal breast cancers (24).To determine whether genomic alterations may also distinguish MBCs (n = 14) from ductal grade 3 carcinomas of basal subtype (ductal carcinomas-BSGIII; n = 19), we compared CNA frequency plots obtained from genomic profiles in these two groups of breast cancers (Supplementary Fig. S10).Unsupervised hierarchical clustering separated the two groups and revealed differences in CNA frequencies (P < 0.05, Fisher's exact test, FDR corrected; data not shown).Gains of chromosomal regions 4q34 and 12p13 and loss of 4q32-33, 11q13, and 15q11 were associated with MBCs,  S4A).FISH analysis of basal (C ) and luminal (D ) cases, B8700 and L9721, respectively, detects this regional amplification found significantly associated with basal tumors.S9A and B. Among genes exhibiting copy number gains (log 2 ratio >0.5), 86 presented gain frequencies different between MBCs and ductal carcinomas-BSGIII (9 for MBCs and 79 for basal non-MBC; Supplementary Table S9A).Only four genes presented copy number losses correlated with MBCs (Supplementary Table S9B).
For the same samples (14 MBCs and 19 ductal carcinomas-BSGIII), gene expression was also analyzed and compared with genomic profiles to identify candidate genes.Only gains associated with overexpression of LOC440731 (1q42.2) were associated with MBCs.
Validation of candidate oncogenes.First, we focused our analysis on gain of 10p, a new alteration in basal breast cancer.Twenty 10p genes were deregulated in relation to copy number gain with a significant difference in relation with basal subtype (Supplementary Table S4A).To begin to validate 10p candidates in our samples (40 basal and 31 luminal), we measured by immunohistochemistry the protein expression of the CDC123/ C10orf7 growth promoter, for which a reliable antibody is available.CDC123 was expressed in 16 basal breast cancers and in none of the luminal samples (Fig. 6A).CDC123 protein expression was associated with basal breast cancer (P = 0.000203, m 2 test) and its level correlated with CDC123 mRNA expression (Spearman test rho = 0.296; P = 0.012), which is itself influenced by the genomic status of the gene (Supplementary Tables S4A and B).The effect of CDC123 protein expression in the 16 basal tumors was not significant on the 5-year MFS and other clinical features (data not shown).No correlation between the CDC123 protein expression and the MBC type was established (P = 0.2464, Fisher's exact test).
Second, we considered important to analyze protein expression from genes targeted by the high-level amplifications because, although rare, they may provide interesting therapeutic targets.As an example, we studied the EGFR amplification in greater details (Supplementary Fig. S8B).EGFR protein expression was analyzed in the same samples (42 basal and 32 luminal) by immunohistochemistry (Fig. 6B).Two basal tumors showed EGFR amplification and EGFR loss, respectively (Supplementary Fig. S8B).EGFR mRNA was expressed at high level in 3 basal cases (including the amplified case), at normal level in 14 basal (35.8%) and 1 luminal cases, and at low level (or inferior to twice the normal breast tissue expression level) in 25 basal (60%) and 31 luminal tumors (96.6%).Membrane positivity was observed in 31 basal cases (79.4%) with a quick score higher than 10, including 17 cases with >50% of positive cells.The basal case with EGFR gene loss did not show EGFR expression.The EGFR-amplified case had the strongest immunohistochemical signal, with the maximum quick score of 300, suggesting a high production of protein as a consequence of gene amplification (Fig. 6B).EGFR protein expression was associated with basal cases (P = 1.496eÀ9 , m 2 test) as reported previously (42).

Discussion
We have described a combined, high-resolution analysis of gene expression and genome copy number in primary human breast cancers of basal and luminal subtypes.It was performed to characterize genomic events and identify candidate cancer genes and potential therapeutic targets that may help identify and treat two major molecular breast cancer subtypes.We used the 244K aCGH matrix from Agilent Technologies, which allows a better resolution than BAC-based aCGH (43).
Genomic alterations and molecular subtypes.We observed the existence of multiple genomic alterations supporting the idea that breast cancers are highly heterogeneous.Several regions and genes targeted by these alterations were different in basal and luminal subtypes.We confirmed some strong associations between subtype and genome profiles reported in previous studies (9)(10)(11); good examples are gains at 12p13 and 11q13 regions associated with basal and luminal subtypes, respectively.We identified gain of 10p as a new alteration of basal breast cancer and found that gains of 8p11-12 targeting the previously defined A1 and A2 amplicons (37) are associated with luminal subtype.
We observed rare high-level amplifications in basal tumors, suggesting again that breast cancers are highly heterogeneous even within a given subtype.They affected small regions, including Left, unsupervised hierarchical clustering of genome copy number profiles measured for the 74 breast tumors by aCGH with 7,694 discriminant probes identified to significantly differentiate basal and luminal tumors and to cover genes commonly analyzed in genomic and gene expression studies.Vertical orange line, limits of the two major clusters.Data are displayed as described in Fig. 1.To the immediate right are plotted the corresponding significant CNA frequency differences observed for luminal (blue ) and basal (red) subtypes.Horizontal lines on frequency plots, chromosome boundaries; horizontal dashed lines, centromere locations.Gene expression of the 1,834 corresponding independent genes was analyzed on the same 74 tumors and allowed to draw an unsupervised hierarchical clustering of genes (shown to the right of the figure) with regard to their respective chromosomal location.The aCGH data from 7,694 probes that revealed significant differences of CNA frequencies between the two subtypes (genomic profiles) were compared with the expression of their 1,834 corresponding independent genes (expression profiles).To the immediate left of the expression profiles are plotted (to the right or to the left part , respectively) the candidate genes that showed significant differences in copy number gain correlated with up-regulated expression or in copy number loss correlated with down-regulated expression between basal (red) and luminal (blue ) subtypes.S4A and S5A, respectively.The most significant candidate oncogenes (A) and TSG (B ) are shown associated with their corresponding chromosome CNA frequency plot defined in basal and luminal tumors.For a better visualization of the targeted region, a threshold value of log 2 ratio >|0.3| was used to draw chromosome CNA frequency plots using CGH Analytics software.For each subtype, all the candidate genes were significantly defined with a FDR-adjusted t test P < 0.05, but some of them were associated with a significance of P < 0.001 (noted with asterisk).Underlined genes have been associated with oncogenesis (see Supplementary Tables S4A and S5A for references).

Cancer Research
Cancer Res 2007; 67: (24).December 15, 2007 PIK3CA (3q26), IGF1R (15q26), and CCNE1 (19q11-12), but also single genes, such as EGFR (7p11), FGFR2 (10q26), and BCL2L2 (14q11).EGFR, FGFR2, EGFR, and IGF1R are tyrosine kinase receptors and thus potential therapeutic targets.Knowledge of CNAs can have immediate clinical use in diagnosis and can, in some cases, provide useful prognostic information (see ref. 44 for review).The existence of these amplifications and such high degree of heterogeneity confirms that molecular profiling will be paramount to select the appropriate treatment.
Specific genomic losses were also detected in basal and luminal subtypes.Our results agree with loss of heterozygosity (LOH) studies: LOH at 4p and 5q defines a subclass of basal breast cancers, whereas LOH at 1p and 16q occurs preferentially in ER-positive breast cancers (45).Losses of 4p and 5q associated with basal breast cancers targeted several genes (Supplementary Tables S2B and C), including candidate or known TSG: SLIT2 (4p15.31),GPR125 (4p15.31),RASA1 (5q14.3),and APC (5q22.2).Genes coding for P53, its repressor PRDM1 (6q21), and its effector PERP (6q23.3)were among genes with losses associated with luminal subtype (Supplementary Tables S2B and C).However, expression of these three genes was not significantly deregulated in relation to copy number loss.
We also detected genomic alterations through analysis of copy number transitions, which could be associated with chromosome breaks.They targeted 501 genes in at least two different tumors.Results were consistent with previous studies (38)(39)(40), underlining the quality and the resolutive power of the aCGH analysis.Ten large genes associated with fragile sites were targeted, with the higher frequency for CSMD1 (41).Copy number transitions with frequencies different between subtypes targeted 6 and 15 genes in basal and luminal tumors, respectively.The difference in the number of involved genes may suggest a different mode of genomic instability in luminal comparatively to basal tumors.Among genes with frequent transitions in basal subtype (Supplementary Table S3B) was PUM2, a human homologue of Pumilio, which codes for a protein required to maintain germ-line stem cells in protostomians (46).Among genes with frequent copy number transitions in luminal samples (Supplementary Table S3B), few are characterized except GLI2 (GLI-Kruppel family member GLI2), CHRD (chordin), DOM3Z, and TNXB.
Genes associated with basal and luminal subtypes or common to the two subtypes.A total of 106 genes showed differences in copy number gain correlated with overexpression Figure 6.CDC123 and EGFR, two basal candidates.A, CDC123 expression was analyzed by immunohistochemistry and quantified.As illustrated with luminal L5342 sample, no luminal case expressed CDC123, whereas medullary B11614 was the most positive basal case.CDC123 was also found expressed in the basal B12478 sample.B, EGFR expression was analyzed by immunohistochemistry and quantified in 12 basal tumors.EGFR expression was negative in the case with EGFR gene loss (B11442), whereas the EGFR -amplified case (B8595) displayed the strongest immunohistochemical signal intensity, suggesting a high production of EGFR protein as a consequence of gene amplification (Supplementary Fig. S8B).Increased EGFR expression was also observed in cases without EGFR amplification (B9752 and B10033).EGFR and CDC123 showed membrane and cytoplasmic signal, respectively.CDC123 and EGFR scoring associated the percentage of positive cells (P , from 0% to 100%) and the intensity of the staining (I , from 0 to 3) as described in Materials and Methods.
between basal (73 genes) and luminal (33 genes) subtypes.The extent of concordance between high-level amplification and increased gene expression has been analyzed previously (47).In basal and luminal subtypes, several genes have been associated with oncogenesis or with the control of various cellular processes (see Supplementary Table S4A and Fig. 5 for details and references).Interestingly, among basal candidate genes, several code for proteins potentially involved in transcription or cell cycle regulation, whereas luminal candidates encode proteins involved in cellular differentiation.Among the eight genes with high expression in relation to high-level amplification in luminal tumors, six have been associated with cancer (Supplementary Table S4C).The genes targeted by low-level CNAs could be selected during oncogenesis to increase basal metabolism of the tumor cells, whereas genes deregulated by high-level amplifications are predominantly involved in signaling and cell cycle regulation (10).IKBKE (1q32.1)has been recently identified as an oncogene in breast cancers (48).In our work, IKBKE was not found among specific and common candidate oncogenes.This indicates that our list is not comprehensive and further supports the great heterogeneity of breast cancer.
Taken together, our and previous studies (10,11) support the existence of a combined activation of several potential oncogenes in 1q21-23, 8p12, 10p14, 11q13, and 12p13 regions involved in specific basal and luminal oncogenic pathways.Our highresolution study pointed to additional candidate oncogenes such as 10p genes and HORMAD1 and ZNF703.
A total of 89 genes with copy number losses correlated with down-regulated mRNA expression were associated with basal (28 genes) and luminal (61 genes) subtypes, respectively.In basal and luminal subtypes, only 7 and 11 genes have been thus far associated with cancer as potential TSG or as coding for proteins involved in growth inhibition pathways (see Supplementary Table S5A and Fig. 5 for details and references), respectively.A total of 237 common candidate cancer genes were associated with both molecular subtypes.These may participate in general oncogenic mechanisms.One possibility is that these genes are targeted in a stem cell background, whereas subtype-specific alterations occur in committed progenitors (49).
Validation of candidate oncogenes.We found that the 10p candidate oncogene CDC123/C10orf7, which is potentially involved in cell cycle regulation, is expressed in basal breast cancer at the protein level (see Supplementary Table S4A).Similarly, we confirmed the strong relation between EGFR protein level and basal subtype.Frequent EGFR positivity has been observed in metaplastic carcinomas, a subset of basal carcinomas (50).Pursuing these analyses with the list of potential candidate breast cancer genes may provide therapeutic targets.
In conclusion, we have studied genome copy number and gene expression in primary human breast cancers.We have identified genomic events that distinguish basal and luminal tumors and have identified various rare high-level amplifications that may provide therapeutic targets.We have characterized gain of 10p as a new alteration in basal breast cancer and have shown that a limited region of 8p11-12 amplicon is found in luminal rather than in basal tumors.Our study pointed to basal and luminal candidate oncogenes, such as HORMAD1 and ZNF703, defined as the best significant basal and luminal potential cancer genes, respectively.We have described 73 candidate oncogenes and 28 candidate TSG as well as 33 candidate oncogenes and 61 candidate TSG that may contribute to the pathophysiology of basal and luminal breast tumors, respectively.At last, among 10p candidate oncogenes associated with basal subtypes, we validated CDC123/C10orf7 protein as a potential basal marker (P = 2 Â 10 À4 ).
Integrated genomics help understand mammary oncogenesis and develop future treatments associating therapeutic targets with various molecular categories of tumors.Such treatments could target candidates associated with a specific subtype, candidates common to all subtypes, and the rare highly amplified candidates with high protein expression level.In this study, six potential oncogenes associated with basal or luminal subtype (Supplementary Table S5A), nine commonly found (Supplementary Table S8), as well as several rare highly amplified candidates are potentially druggable.

Figure 1 .
Figure 1.Recurrent abnormalities in 93 primary breast tumors.Left, unsupervised hierarchical clustering of genome copy number profiles measured for 93 primary breast tumors by aCGH on 225.388 probes (without X and Y).Red, increased genome copy number; green, decreased genome copy number.Bar to the left, chromosome locations with chromosome 1pter to the top and 22qter to the bottom.Vertical orange line, limits of the two major clusters.The location of the odd-numbered chromosomes is indicated.Below the dendrogram, the first seven rows indicate biological and clinical features of the tumors, whereas the last row indicates subtypes.Color codes and corresponding legends are indicated in the box located to the right at the top of the figure.MED, medullary; LOB, lobular; MIX, mixed lobular and ductal; DUC, ductal; NA, nonassigned.CNA frequencies discriminated basal (red ) and luminal (blue ) tumors.Right, frequencies of genome copy number gains and losses are plotted as a function of chromosome location with the same previous sequential order for, from the left to the right, both (global) basal and luminal subtypes.Horizontal lines, chromosome boundaries.Positive and negative values indicate frequencies of tumors showing copy number increase and decrease, respectively, with gains and losses as described in Materials and Methods.

Figure 2 .
Figure 2. 12p13.3 gain is a significant genomic signature of basal tumors.A, aCGH profiles of chromosome 12 show 12p13.3gain frequency in basal tumors.Profiles were established with CGH Analytics software.Vertical green bars, location of three FISH BAC probes mapped in the 12p13.31region centered on the NOL1 gene (B) and previously used in FISH analysis (33).B, from top to bottom and from telomere to centromere is shown the combination of BAC clones (RP5-940J5, RP11-433J6, and RP11-578M14; green boxes ) previously used as biotinylated FISH probe (revealed in green, FITC) on breast tumors organized on a tissue microarray; Mb scale of the corresponding 12p13.31region corresponds to cytogenetic bands; genes are taken from the Build 36.1 from NCBI (March 2006 version).Red-colored genes MRPL51, CNAP1, GAPDH, NOL1, MLF2, CDCA3 , and ENO2 exhibit gain copy number correlated with overexpression in basal tumors (see Supplementary TableS4A).FISH analysis of basal (C ) and luminal (D ) cases, B8700 and L9721, respectively, detects this regional amplification found significantly associated with basal tumors.

Figure 3 .
Figure 3. 10p gain is a new significant genomic signature of basal tumors.For each basal (left part ) and luminal (right part ) cases, the regional genomic profiles were established, with CGH Analytics software, within the genomic interval (13.1-0Mb) of the short arm of chromosome 10 (hg17 human genome mapping; Build 35 from NCBI, May 2004 version).Color profiles correspond to different cases.Several basal cases showed 10p13-15.3gain or amplification (A), whereas no luminal case displayed this alteration (B ).The various profiles associated with 10p high-level amplification of the three basal cases B8667, B11895, and B11442 suggested the presence of several amplicons.

Figure 4 .
Figure4.Correlation between gene expression and genome alterations.Genomic and gene expression profiles were established for 74 breast tumors (42 basal and 32 luminal).Left, unsupervised hierarchical clustering of genome copy number profiles measured for the 74 breast tumors by aCGH with 7,694 discriminant probes identified to significantly differentiate basal and luminal tumors and to cover genes commonly analyzed in genomic and gene expression studies.Vertical orange line, limits of the two major clusters.Data are displayed as described in Fig.1.To the immediate right are plotted the corresponding significant CNA frequency differences observed for luminal (blue ) and basal (red) subtypes.Horizontal lines on frequency plots, chromosome boundaries; horizontal dashed lines, centromere locations.Gene expression of the 1,834 corresponding independent genes was analyzed on the same 74 tumors and allowed to draw an unsupervised hierarchical clustering of genes (shown to the right of the figure) with regard to their respective chromosomal location.The aCGH data from 7,694 probes that revealed significant differences of CNA frequencies between the two subtypes (genomic profiles) were compared with the expression of their 1,834 corresponding independent genes (expression profiles).To the immediate left of the expression profiles are plotted (to the right or to the left part , respectively) the candidate genes that showed significant differences in copy number gain correlated with up-regulated expression or in copy number loss correlated with down-regulated expression between basal (red) and luminal (blue ) subtypes.

Figure 5 .
Figure 5.Most significant candidate cancer genes in basal and luminal tumors.One hundred six candidate oncogenes and 89 putative TSGs are listed in Supplementary TablesS4A and S5A, respectively.The most significant candidate oncogenes (A) and TSG (B ) are shown associated with their corresponding chromosome CNA frequency plot defined in basal and luminal tumors.For a better visualization of the targeted region, a threshold value of log 2 ratio >|0.3| was used to draw chromosome CNA frequency plots using CGH Analytics software.For each subtype, all the candidate genes were significantly defined with a FDR-adjusted t test P < 0.05, but some of them were associated with a significance of P < 0.001 (noted with asterisk).Underlined genes have been associated with oncogenesis (see Supplementary Tables S4A and S5A for references).