
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Human Cancer Biology |
Authors' Affiliations: 1 Institut Curie, Department of Tumor Biology, 2 Institut National de la Sante et de la Recherche Medicale Unit 830, Institut Curie, 3 Institut Curie, Translational Research Department, 4 Institut Curie, UMR144 Centre National de la Recherche Scientifique, 5 Institut Curie, Department of Surgery, 6 Head of the Breast Cancer Study Group of the Institut Curie, and 7 Institut Curie, Department of Radiation Oncology, Paris, France
Requests for reprints: Olivier Delattre, Institut National de la Sante et de la Recherche Medicale Unit 830, Institut Curie, 26 rue d'Ulm, 75248 Paris Cedex 05, France. Phone: 33-1-42-34-66-81; Fax: 33-1-42-34-66-30; E-mail: olivier.delattre{at}curie.fr.
| Abstract |
|---|
|
|
|---|
Experimental Design: We did a combined phenotypic and genomic analysis of a series of 57 DCIS integrated with gene expression profile analysis for 26 of the 57 cases.
Results: Thirty-two DCIS exhibited a luminal phenotype; 21 were ERBB2 positive, and 4 were ERBB2/estrogen receptor (ER) negative with 1 harboring a bona fide basal-like phenotype. Based on a CGH analysis, genomic types were identified in this series of DCIS with the 1q gain/16q loss combination observed in 3 luminal DCIS, the mixed amplifier pattern including all ERBB2, 12 luminal and 2 ERBB2-/ER- DCIS, and the complex copy number alteration profile encompassing 14 luminal and 1 ERBB2-/ER- DCIS. Eight cases (8 of 57; 14%) presented a TP53 mutation, all being amplifiers. Unsupervised analysis of gene expression profiles of 26 of the 57 DCIS showed that luminal and ERBB2-amplified, ER-negative cases clustered separately. We further investigated the effect of high and low copy number changes on gene expression. Strikingly, amplicons but also low copy number changes especially on 1q, 8q, and 16q in DCIS regulated the expression of a subset of genes in a very similar way to that recently described in invasive ductal carcinomas.
Conclusions: These combined approaches show that the molecular heterogeneity of breast ductal carcinomas exists already in in situ lesions and further indicate that DCIS and invasive ductal carcinomas share genomic alterations with a similar effect on gene expression profile.
DCIS correspond to a proliferation of malignant epithelial cells that spare the basal membrane with no stromal invasion. They can be classified into three categories (low-, intermediate-, and high-grade) according to nuclear grade, differentiation, and presence of central necrosis (2). However, this classification of DCIS is poorly reproducible among pathologists (3). Molecular studies, still based on a limited number of cases, have revealed some characteristics of high-grade DCIS, including the high frequency of ERBB2 overexpression, TP53 mutations, chromosomal 8p loss, and 1q gain (4–9). In contrast, positive estrogen receptors (ER), BCL2 expression, and 16q loss have been reported to characterize low-grade DCIS (10–12).
Treatment, regardless of DCIS nuclear grade, is based on mastectomy or lumpectomy in combination with radiotherapy (1, 13). Factors known to be associated with a higher risk of recurrence are young age, positive margins, necrosis, and high nuclear grade (14, 15). Some emerging data suggest that ERBB2-amplified DCIS could present a higher risk of recurrence (16, 17).
Recent gene expression profiling studies have reported that consistent differences in expression of a subset of genes can be identified between low-grade and high-grade DCIS (18–20). Interestingly, in situ lesions preferentially cluster with invasive lesions of the same grade, suggesting that DCIS are precursors of their similar grade invasive counterpart (20). In support of this hypothesis, DCIS and invasive components of the same tumors exhibit similar patterns of genetic alterations (4, 7).
Gene expression profile analysis classified invasive breast carcinomas into five groups: luminal A or B, basal like, ERBB2, and normal like carcinomas (21, 22). These analyses have led to the identification of specific immunophenotypic markers, which allow this classification to be used in clinical practice (23, 24). It is well-recognized that ERBB2 is up-regulated more frequently in DCIS than in invasive carcinomas (6). The identification of a basal-like entity is emerging with four recent studies describing DCIS lesions with this phenotype (18, 25–27). It nevertheless remains to be determined whether basal-like DCIS show specific genomic/transcriptomic alterations.
Genomic profiles also provide criteria for a clinically relevant classification, as combined genomic and gene expression studies in invasive cancers have highlighted the prognostic effect of high level amplifications associated with gene overexpression at chromosome regions 8p12, 11q13, 17q12, and 20q13 independently of the previously described major gene expression classes. These cases are defined as amplifiers. Among nonamplifiers, one group with a good prognosis harbors few genetic alterations but recurrent 16q loss associated with 1q gain (1q/16q). Another group, which includes most of basal-like but also a subset of luminal cases, is characterized by numerous low-level copy number alterations (CNA; refs. 28, 29). In invasive breast cancer, this integrated genomic approach, based on genomic and transcriptomic data, is very promising to achieve better clinical management of patients.
In the present study, we analyzed the immunophenotype, the array CGH profiles, the TP53 gene sequence, and the gene expression profile of DCIS frozen samples. We identified DCIS subtypes according to the molecular classification described for invasive ductal carcinomas (IDC), allowing comparison of the genomic, transcriptomic, and phenotypic profiles of ERBB2, luminal, and basal-like DCIS. Finally, correlations of genomic and transcriptomic data were used to identify candidate genes, the altered expression of which may play a critical role in early breast carcinogenesis via copy number changes.
| Materials and Methods |
|---|
|
|
|---|
|
Immunohistochemistry. Immunostaining was done according to previously published protocols (31). The expression of ER (clone 6F11; 1/200; Novocastra), progesterone receptor (PR; clone 1A6; 1/200; Novocastra), ERBB2 (clone CB11; 1/1,000; Novocastra), epidermal growth factor receptor (HER1; clone 31G7; 1/40; Zymed; Clinisciences), cytokeratin 5/6 (clone D5/16B4; 1/50; Dako), and cytokeratin 8/18 (clone DC10; 1/100; Zymed; Clinisciences) were evaluated. For each antibody, internal and external controls were included in the experiments.
TP53 sequencing. Exons 4 to 10 of TP53 were PCR amplified and bidirectionally sequenced using Big Dye Terminator chemistry (Applied Biosystems) with an ABI PRISM 3700 DNA Analyzer (primer sequences previously published in ref. 32).
Array CGH. The 3.5K BAC array together with experimental procedures used for hybridization and washing were as previously published (32). All BAC and PAC were verified by end sequencing before spotting.
Analysis of array CGH data. Data analysis was based on the normalized ratios of Cy5/Cy3 signals observed for each BAC clone that previously passed the flag assessment procedure (32, 33). For autosomal chromosomes, the status of each BAC locus was defined according to the VAMP (Vizualization and Analysis of Molecular Profiles) analysis procedure (34) using the GLAD (Gain and Loss Analysis of DNA) algorithm that, based on the Adaptive Weights Smoothing procedure, can be used to assign a status (gained, lost, or normal) to each chromosome region (33).
Regions of gain or loss and amplicons were defined as previously published (32). BACs known to harbor copy number variation in public databases (Database of Genomic variants)8 or identified in the laboratory by array CGH experiments using normal/normal DNAs hybridizations were excluded from the analysis.
Gene expression profile analysis. The DNA microarrays used in this study were the Human Genome U133 Set (HG-U133; Affymetrix), consisting of one GeneChip array (A) and containing almost 20,000 probe sets. Microarray data were simultaneously normalized using the GCRMA package version 1.1.4 in the R environment (R Development Core Team). Unsupervised hierarchical clustering of the tumor sample was done using the dChip software (35) with standard Pearson correlation as similarity measure and centroid as linkage criteria. Clustering was done on a subset of probe sets with an interquartile range of >2. Each probe set was preliminary standardized with respect to mean and SD. The comparison of mean gene expression between ERBB2-amplified and luminal DCIS was done using SAM (36) and Welch t tests. FDR for Welch t test was adjusted using the Benjamini-Hochberg corrections. Finally, for each probe set, the correlation with Cy5 to Cy3 ratios of the corresponding BAC clone was calculated using Pearson and Spearman correlation tests (for details see Supplementary Data 1).
| Results |
|---|
|
|
|---|
|
CNA (amplicons, gains, and losses) were identified using a one-megabase resolution BAC array. All samples but four exhibited altered profiles, indicating that the amount of stromal or other normal cells in the samples was not sufficient to mask copy number changes. Fifty-nine different amplicons were observed in the 57 DCIS (Supplementary Data 2). Fourteen amplicons were recurrent (observed in at least two tumors), including ERBB2 on chromosome 17q12, Cyclin D1 at 11q13, FGFR1 at 8p12, BMP7/BCAS1 at 20q13, but also five 17q amplicons distinct from ERBB2. Interestingly, the chromosomal localization of amplicons (distinct from ERBB2) differed between the ERBB2-amplified and the luminal DCIS groups, as half of the amplicons (13 of 26) of the ERBB2 group associated with ERBB2 were also localized on chromosome 17q, but only 5 of 29 amplicons of the luminal group were localized on this chromosome arm (P = 2.10–2).
The ERBB2 amplicon size differed from one tumor to another, ranging from 32.9 106 bp to 35.6 106 bp with a minimal common region of amplification located between 34.9 106 bp and 35.2 106 bp centered on the ERBB2 gene sequence but also containing other genes including STARD3, GRB7, and PERLD1.
Regions of gains and losses observed with a frequency of >20% in the luminal and ERBB2-amplified groups were identified. Common regions of gain were observed in both groups (Fig. 1C and D). They were located on chromosomes 1p, 1q, 5p, 7p, 7q, 8q, 14q, 16p, 19p, 19q, 17q, and 20q. Only one common and recurrent region of loss on 11q was identified. The ERBB2-amplified group presented specific regions of 3p, 4p, 4q, and 8p losses and 17q gains compared with the luminal group that presented specific regions of gains on 1q, 8p, and 17q and losses on 16q (Supplementary Data 3).
According to the criteria recently proposed by Fridlyand and Chin in 2006 (28, 29), DCIS cases could be classified into three categories: tumors characterized by the presence of amplicons (mixed amplifiers), tumors characterized by 1q gain and 16q loss (1q/16q group), and tumors with complex CNA. ERBB2-amplified DCIS never showed 16q loss. The luminal group exhibited genetic complexity with 12 cases classified as amplifiers, 3 cases harboring a 1q/16q profile, 14 cases exhibiting a complex pattern of CNA, and 3 cases characterized by a flat profile. Within the ERBB2-/ER-negative DCIS group, two belonged to the amplifier category with amplicons located on chromosomes 5p, 8p, and 22q, one harbored CNA and one case was charade rised by a flat profile (Fig. 1A and B).
Gene expression analysis. We also investigated transcriptomic differences between ERBB2-amplified and luminal DCIS. A subset of 26 cases for which good quality RNAs were available, including 9 ERBB2-amplified and 17 luminal DCIS, were analyzed on Affymetrix U133A. A total of 502 genes exhibited an interquartile range of >2 (Supplementary Data 4). These genes discriminated two major clusters of tumors mainly organized according to ERBB2 and hormonal receptor status rather than nuclear grade. One branch included most (8 of 9) ERBB2-amplified, ER-negative cases, and the other branch comprised all (17 of 17) luminal cases and one ERBB2-amplified, ER-negative/PR-positive case. As shown in Fig. 2 , genes coamplified with ERBB2 (GRB7, PERLD1, and STARD3) and genes of the ER pathway (ESR1, GATA3, STC2, LIV1, IGF1R...) were major contributors of this clustering. The most discriminant genes between ERBB2-amplified and luminal DCIS (Welch and Sam tests results, P < 0.01) are indicated in Supplementary Data 5. Unsupervised and supervised analyses showed that genes involved in estrogens receptor pathways, extracellular matrix remodeling (syndecan, ezrin, ...), immune response (TRBV19, ...), antiapoptotic pathways, cell adhesion (CXCR4, syndecan, and ITGB6), or fatty acid metabolism (LDRL) were differentially expressed between the ERBB2-amplified and luminal groups (Fig. 2). Finally, genes encoding important signal transduction and transcription factors molecules were differentially expressed. FGFR2, RET, KIF5, and IGF1R were overexpressed in luminal cases, whereas ERBB2, GRB7, PTPRC were overexpressed in ERBB2-amplified cases.
|
Genomic and transcriptomic correlations. Given the strong relationship between the genomic classification and the unsupervised clustering of gene expression profiles and to determine the effect of genomic changes on gene expression levels, we combined the analyses of the Affymetrix signals and array CGH ratios for the 26 DCIS analyzed with both approaches. Using the Pearson and Spearman correlation tests, we observed that the expression level of 416 genes was significantly correlated with the ratio of their corresponding BAC locus (372 BAC clones; Fig. 3 ). Correlations corresponding to amplicon were mostly located on 17q and corresponding to either the variable extent of the ERBB2 amplicon from one tumor to another or to other 17q amplicons, distinct from but co-occurring with ERBB2 (Supplementary Data 2). Other correlations due to amplicons mainly involved 8p (from 37 to 40 Mb), 11p12-p13 (from 34 to 36 Mb), 11q13 (from 69 to 78 Mb), and 20q13.32 (from 46 to 55 Mb; Supplementary Data 6). The main correlations due to low copy gains or losses were observed on chromosomes 1q, 8q, and 16q (Table 2 ; Supplementary Data 6). The analysis of genome expression correlations indicates that the amount of stromal or other normal cells in the samples was not sufficient to hamper the detection of a strong copy number effect upon gene expression.
|
|
| Discussion |
|---|
|
|
|---|
DCIS can be classified into the luminal, ERBB2, and basal-like groups previously described in IDC (21, 22). In our series of 57 DCIS, most tumors (56%) belonged to the luminal group based on ER expression and the absence of ERBB2 overexpression. The second group of tumors comprised ERBB2 overexpressing DCIS (37%) that were mostly high grade (81%) and ER negative/PR negative (90%) in agreement with previous reports showing that ERBB2 overexpression is rare in ER-positive DCIS (37). A third small group was composed of three triple-negative DCIS (5%) with only one case demonstrating a clear basal-like immunophenotype (2%) with positive staining for keratin 5/6 and epidermal growth factor receptor, confirming that the recently recognized basal-like DCIS group is a rare entity among pure DCIS (6-10% of cases; refs. 25, 26, 38, 39). This was also shown by a recent gene expression profile analysis of a series of DCIS (18).
When compared with the luminal group, the ERBB2-amplified cases presented a particular array CGH profile with more frequent 17q gains, 3p, 4p, 4q, and 8p losses. In contrast, the luminal DCIS group harbored significantly more frequent 16q losses. These observations refine previously published classic CGH data showing a significantly higher number of 17q gains in high nuclear grade DCIS (which correspond to most cases of our ERBB2-amplified, ER-negative group) and 16q losses in low-grade DCIS (which are represented in the present series by the luminal cases; refs. 4, 7, 11).
We also show that DCIS can be classified according to the genomic criteria recently proposed for IDC (28), as a first group of 34 tumors presented high-level amplification of genomic regions (amplifiers). All TP53 mutations were observed in this group. Remarkably, the recurrent amplicons of DCIS were located on 8p12 (FGFR1), 11q12 (Cyclin D1), 17q12 (ERBB2), and 20q13 (BMP7/BCAS1), corresponding also to the most frequent amplicons of IDC (28). Compared with other amplifier cases and as recently observed in IDC (28), half amplicons associated with ERBB2 were also located on chromosome 17q. In contrast, in the luminal group, only 17% of amplicons involved 17q. A second group comprised tumors with a 16q deletion associated with 1q gain but no additional CNA. All these cases presented a luminal phenotype. The last group comprised cases with more complex CNA. In this series, the three triple-negative DCIS did not present any recurrent genomic profile, but two were amplifiers including the bona fide basal-like DCIS, which also harbored a TP53 gene mutation.
Unsupervised analysis of gene expression profiles supports the observation that ERBB2-amplified and luminal DCIS constitute two distinct entities. Gene ontology analysis of the ERBB2 DCIS signature showed that the genes up-regulated in this group belonged to ERBB2 signal transduction pathways but also to fatty acid metabolism and other protein kinase and transmembrane signaling pathways. Genes down-regulated mainly belonged to apoptotic and ER regulation pathways. The list of genes differentially expressed between ERBB2 and luminal DCIS in this series was compared with the recently published list of genes that discriminate DCIS according to their differentiation status (18). Although the two comparisons were based on slightly different criteria, they both highlighted major genes involved in breast carcinogenesis such as ESR1, GRB7, and ERBB2. A link between the ERBB2 pathway and fatty acid metabolism has been previously reported in cell lines (40) and invasive breast carcinomas (41). The present study also showed that these pathways are activated right from the earliest stages of ERBB2-driven breast carcinogenesis. Interestingly, as reported in IDC (42, 43), overexpression of 17q12 genes STARD3, GRB7, and PERLD1 is also associated with ERBB2 overexpression in DCIS, which raises the question of their putative synergistic role during early breast carcinogenesis. Finally, Integrin β6, a receptor for fibronectin and cytotactin, is overexpressed in the ERBB2 DCIS group. Recent data have highlighted molecular interactions between ERBB2 receptors and integrin networks (41, 44). Moreover, experiments conducted with ITGB6 knockout mice indicate that ITGB6 may regulate transforming growth factor β and metalloprotease, in particular matrix metalloproteinase 12 expression (45). Further investigations are needed to understand the role of ITGB6 in ERBB2-amplified DCIS.
It is noteworthy that genes recently involved in angiogenesis or invasion such as SDC1 (syndecan; refs. 46–48) or CXCR4 are up-regulated in ERBB2-positive in situ lesions (49). Reciprocally, FGFR2, a gene recently linked to in susceptibility locus in breast cancer (50, 51), is up-regulated in luminal DCIS.
The correlation between array CGH and expression data shows that amplifications but also more subtle copy number changes lead to modifications of the gene expression profile. Strikingly, a large part of the genes, for which a strong correlation between expression level and genomic status was observed, also exhibited a similar correlation in IDC (28), indicating that genome alterations have similar consequences on gene expression in both IDC and DCIS. In line with the data reported by Chin et al. (28), these results identify genes located in frequently altered regions in breast cancer, strongly influenced by copy number variations, and whose gene dosage effects may play a role in early phases of breast carcinogenesis. Two genes on chromosome 16 seem to be of particular interest: ASCIZ (ATM/ATR-substrate ChK2-interacting Zinc finger; ref. 52), which encodes a protein involved in lesion-specific Rad51 focus formation, a process that also involves BRCA1 and BRCA2, two master genes in breast carcinogenesis; and CNOT1, which encodes a member of the CCR4-NOT complex, a global regulator of transcription that may act as a repressor of endogenous estrogens target genes (53). Interestingly, CNOT7, another member of the CCR4-NOT complex, is regulated by gene copy number variation on chromosome 8p, a chromosome frequently altered in ERBB2 DCIS tumors. Based on these observation, the hypothesis of haploinsufficiency of the CCR4-NOT complex playing a role in early breast tumorigenesis may now be tested.
Altogether, these integrated genomic and transcriptomic data confirm previous results showing strong similarities between IDC and DCIS (20, 54). Although we cannot exclude the possibility that subtle genetic alterations that may have escaped detection in this study could also play a role in the evolution from DCIS to IDC, this study raises the hypothesis that the difference between DCIS and IDC is not based on specific genetic alterations or on their consequences on gene expression profile. Because a recent report suggested that a 35-gene expression classifier might distinguish IDC and DCIS (18), the role of epigenetic modifications in the progression from DCIS to IDC will need to be evaluated.
In conclusion, our data show that DCIS already displays the molecular diversity observed in IDC and, therefore, can be classified according to molecular criteria distinguishing ERBB2-amplified DCIS, usually high-grade, ER negative with frequent TP53 mutations, from luminal DCIS corresponding to low/intermediate-grade, ER positive with a very low rate of TP53 mutations. In our series, only three cases were classified as triple negative and only one was classified as basal-like DCIS, which confirms that this last entity is rare among DCIS. Further studies are needed to define whether a classification of DCIS based on molecular markers may help to more accurately define cases associated with a higher risk of recurrence. And finally, genomic/transcriptomic correlations represent a promising tool to identify new genes and pathways important in early breast carcinogenesis.
| Acknowledgments |
|---|
Members of the Breast Cancer Study group in addition to coauthors: Bernard Asselain, Alain Aurias, Emmanuel Barillot, Marc Bollet, François Campana, Paul Cottu, Patricia de Cremoux, Véronique Diéras, Laurent Mignot, Jean-Yves Pierga, Marie-France Poupon, Dominique Stoppa-Lyonnet, Anne Tardivon, Fabienne Thibault, and Pascale This.
| Footnotes |
|---|
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.
Note: Supplementary data for this article are available at Clinical Cancer Research Online (http://clincancerres.aacrjournals.org/).
8 http://projects.tcag.ca/variation/ ![]()
Received 6/14/07; revised 10/28/07; accepted 11/29/07.
| References |
|---|
|
|
|---|
(v)β6-mediated TGF-β activation causes Mmp12-dependent emphysema. Nature 2003;422:169–73.[CrossRef][Medline]This article has been cited by other articles:
![]() |
A. Vincent-Salomon and O. Delattre Role of Genetic, Genomic, and Transcriptomic Factors Am. Assoc. Cancer Res. Educ. Book, April 18, 2009; 2009(1): 97 - 103. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |