High-resolution genomic and expression analyses of copy number alterations in HER2-amplified breast cancer

Introduction

HER2 gene amplification and protein overexpression (HER2+) define a clinically challenging subgroup of breast cancer with variable prognosis and response to therapy. Although gene expression profiling has identified an ERBB2 molecular subtype of breast cancer, it is clear that HER2+ tumors reside in all molecular subtypes and represent a genomically and biologically heterogeneous group, needed to be further characterized in large sample sets.

Methods

Genome-wide DNA copy number profiling, using bacterial artificial chromosome (BAC) array comparative genomic hybridization (aCGH), and global gene expression profiling were performed on 200 and 87 HER2+ tumors, respectively. Genomic Identification of Significant Targets in Cancer (GISTIC) was used to identify significant copy number alterations (CNAs) in HER2+ tumors, which were related to a set of 554 non-HER2 amplified (HER2-) breast tumors. High-resolution oligonucleotide aCGH was used to delineate the 17q12-q21 region in high detail.

Results

The HER2-amplicon was narrowed to an 85.92 kbp region including the TCAP, PNMT, PERLD1, HER2, C17orf37 and GRB7 genes, and higher HER2 copy numbers indicated worse prognosis. In 31% of HER2+ tumors the amplicon extended to TOP2A, defining a subgroup of HER2+ breast cancer associated with estrogen receptor-positive status and with a trend of better survival than HER2+ breast cancers with deleted (18%) or neutral TOP2A (51%). HER2+ tumors were clearly distinguished from HER2- tumors by the presence of recurrent high-level amplifications and firestorm patterns on chromosome 17q. While there was no significant difference between HER2+ and HER2- tumors regarding the incidence of other recurrent high-level amplifications, differences in the co-amplification pattern were observed, as shown by the almost mutually exclusive occurrence of 8p12, 11q13 and 20q13 amplification in HER2+ tumors. GISTIC analysis identified 117 significant CNAs across all autosomes. Supervised analyses revealed: (1) significant CNAs separating HER2+ tumors stratified by clinical variables, and (2) CNAs separating HER2+ from HER2- tumors.

Conclusions

We have performed a comprehensive survey of CNAs in HER2+ breast tumors, pinpointing significant genomic alterations including both known and potentially novel therapeutic targets. Our analysis sheds further light on the genomically complex and heterogeneous nature of HER2+ tumors in relation to other subgroups of breast cancer.     

Comprehensive copy number profiles of breast cancer cell model genomes

Introduction

Breast cancer is the most commonly diagnosed cancer in women worldwide and consequently has been extensively investigated in terms of histopathology, immunochemistry and familial history. Advances in genome-wide approaches have contributed to molecular classification with respect to genomic changes and their subsequent effects on gene expression. Cell lines have provided a renewable resource that is readily used as model systems for breast cancer cell biology. A thorough characterization of their genomes to identify regions of segmental DNA loss (potential tumor-suppressor-containing loci) and gain (potential oncogenic loci) would greatly facilitate the interpretation of biological data derived from such cells. In this study we characterized the genomes of seven of the most commonly used breast cancer model cell lines at unprecedented resolution using a newly developed whole-genome tiling path genomic DNA array.

Methods

Breast cancer model cell lines MCF-7, BT-474, MDA-MB-231, T47D, SK-BR-3, UACC-893 and ZR-75-30 were investigated for genomic alterations with the submegabase-resolution tiling array (SMRT) array comparative genomic hybridization (CGH) platform. SMRT array CGH provides tiling coverage of the human genome permitting break-point detection at about 80 kilobases resolution. Two novel discrete alterations identified by array CGH were verified by fluorescence in situ hybridization.

Results

Whole-genome tiling path array CGH analysis identified novel high-level alterations and fine-mapped previously reported regions yielding candidate genes. In brief, 75 high-level gains and 48 losses were observed and their respective boundaries were documented. Complex alterations involving multiple levels of change were observed on chromosome arms 1p, 8q, 9p, 11q, 15q, 17q and 20q. Furthermore, alignment of whole-genome profiles enabled simultaneous assessment of copy number status of multiple components of the same biological pathway. Investigation of about 60 loci containing genes associated with the epidermal growth factor family (epidermal growth factor receptor, HER2, HER3 and HER4) revealed that all seven cell lines harbor copy number changes to multiple genes in these pathways.

Conclusion

The intrinsic genetic differences between these cell lines will influence their biologic and pharmacologic response as an experimental model. Knowledge of segmental changes in these genomes deduced from our study will facilitate the interpretation of biological data derived from such cells.     

Identification of specific gene copy number changes in asbestos-related lung cancer

Asbestos is a well-known lung cancer-causing mineral fiber. In vitro and in vivo experiments have shown that asbestos can cause chromosomal damage and aberrations. Lung tumors, in general, have several recurrently amplified and deleted chromosomal regions. To investigate whether a distinct chromosomal aberration profile could be detected in the lung tumors of heavily asbestos-exposed patients, we analyzed the copy number profiles of 14 lung tumors from highly asbestos-exposed patients and 14 matched tumors from nonexposed patients using classic comparative genomic hybridization (CGH). A specific profile could lead to identification of the underlying genes that may act as mediators of tumor formation and progression. In addition, array CGH analyses on cDNA microarrays (13,000 clones) were carried out on 20 of the same patients. Classic CGH showed, on average, more aberrations in asbestos-exposed than in nonexposed patients, and an altered region in chromosome 2 seemed to occur more frequently in the asbestos-exposed patients. Array CGH revealed aberrations in 18 regions that were significantly associated with either of the two groups. The most significant regions were 2p21-p16.3, 5q35.3, 9q33.3-q34.11, 9q34.13-q34.3, 11p15.5, 14q11.2, and 19p13.1-p13.3 (P < 0.005). Furthermore, 11 fragile sites coincided with the 18 asbestos-associated regions (P = 0.08), which may imply preferentially caused DNA damage at these sites. Our findings are the first evidence, indicating that asbestos exposure may produce a specific DNA damage profile.     

A genomic copy number signature predicts radiation exposure in post‐Chernobyl breast cancer

Breast cancer is the second leading cause of cancer death among women worldwide and besides life style, age and genetic risk factors, exposure to ionizing radiation is known to increase the risk for breast cancer. Further, DNA copy number alterations (CNAs), which can result from radiation‐induced double‐strand breaks, are frequently occurring in breast cancer cells. We set out to identify a signature of CNAs discriminating breast cancers from radiation‐exposed and non‐exposed female patients. We analyzed resected breast cancer tissues from 68 exposed female Chernobyl clean‐up workers and evacuees and 68 matched non‐exposed control patients for CNAs by array comparative genomic hybridization analysis (aCGH). Using a stepwise forward–backward selection approach a non‐complex CNA signature, that is, less than ten features, was identified in the training data set, which could be subsequently validated in the validation data set (p value < 0.05). The signature consisted of nine copy number regions located on chromosomal bands 7q11.22‐11.23, 7q21.3, 16q24.3, 17q21.31, 20p11.23‐11.21, 1p21.1, 2q35, 2q35, 6p22.2. The signature was independent of any clinical characteristics of the patients. In all, we identified a CNA signature that has the potential to allow identification of radiation‐associated breast cancer at the individual level.     

Estimating the number of rate limiting genomic changes for human breast cancer

We used multistage models that incorporate the age dependent dynamics of normal breast tissue, clonal expansion of intermediate cells and mutational events to fit data for the age-specific incidence of breast cancers in the surveillance, epidemiology, and end results (SEER) registry. Our results suggest that two or three rate limiting events occurring at rates characteristic of point mutation rates for normal mammalian cells set in motion a sequence of other genomic changes that lead with high probability to breast carcinoma.     

High-resolution analyses of copy number changes in disseminated tumor cells of patients with breast cancer

The presence of disseminated tumor cells (DTCs) in bone marrow (BM) identifies breast cancer patients with less favorable outcome. Furthermore, molecular characterization is required to investigate the malignant potential of these cells. This study presents a single-cell array comparative genomic hybridization (SCaCGH) method providing molecular analysis of immunomorphologically detected DTCs. The resolution limit of the method was estimated using the cancer cell line SK-BR-3 on 44 and 244k arrays. The technique was further tested on 28 circulating tumor cells and four hematopoietic cells (HCs) from peripheral blood (n = 8 patients). The SCaCGH method was finally applied to 24 DTCs, three immunopositive cells morphologically classified as probable HCs from breast cancer patients and five HC controls from BM (n = 7 patients plus n = 1 healthy donor). The frequency of copy number changes of the DTCs revealed similarities with primary breast tumor samples. Three of the patients had available profiles for DTCs and the corresponding tumor tissue from primary surgery. More than two-third of the analyzed DTCs disclosed equivalent changes, both to each other and to the corresponding primary disease, whereas the rest of the cells showed balanced profiles. The probable HCs revealed either balanced profiles (n = 2) or changes comparable to the tumor tissue and DTCs (n = 1), indicating morphological overlap between HCs and DTCs. Similar aberration patterns were visible in DTCs collected at diagnosis and at 3 years relapse-free follow-up. SCaCGH may be a powerful tool for the molecular characterization of DTCs.     

Germline copy number variants are not associated with globally acquired copy number changes in familial breast tumours

A characteristic of sporadic and familial breast tumours is genomic instability, resulting from either inherited mutations in genes that control genome integrity or mutations that are acquired in somatic cells during development. It is well established that abnormal chromosome number and structural changes to chromosomes play an important role in the cause and progression of breast cancer. Familial BRCA1 breast tumours are characterised by basal-like phenotype and high-histological grade which are typically associated with increased genomic instability. Consistent with previous studies, the genomes with the greatest number of base pairs covered by copy number change were typically found in basal-like and/or high-histological grade breast tumours within our cohort. Moreover, we show that luminal A tumours that are high grade had significantly less copy number variant (CNV) coverage than the more clinically aggressive high-grade luminal B tumours, suggesting that chromosomal instability rather than cellular differentiation contributes to the aggressive nature of luminal B tumours. It has previously been proposed that germline CNVs may contribute to somatically acquired chromosome changes in the tumour, but this is the first study to address this idea in breast cancer. By comparing germline CNVs and tumour-specific CNVs in matched breast tumour and normal tissue using data from the Illumina Human CNV370 duo beadarray, we provide evidence that germline CNVs do not tend to act as a foundation on which larger chromosome copy number aberrations develop in tumour cells. Further studies are required with increased sequence resolution that will detect smaller CNVs and define CNV breakpoints to comprehensively assess the relationship between inherited genomic variation and genome evolution in breast cancer.     

Protein expression and gene copy number changes of receptor tyrosine kinase in thymomas and thymic carcinomas

BackgroundInsulin-like growth factor-1 receptor (IGF-1R), epidermal growth factor receptor (EGFR), human epidermal growth factor receptor-type 2 (HER2), and c-Met are members of the receptor tyrosine kinases (RTKs). The associations between the RTK status [protein expression and gene copy number (GCN)] and patient characteristics and between the RTK status and prognosis remain undetermined.Materials and methodsThe study included 140 patients who underwent surgery for thymic tumors. Protein expression was evaluated by immunohistochemistry (IHC) and GCN was evaluated by bright-field in situ hybridization (BISH). The correlations between the RTK status and clinicopathological findings were examined.ResultsIGF-1R protein was frequently detected in thymic carcinoma (83.8%) and EGFR in thymic tumors (91.4%). Thirty-six and 39 tumors were BISH high for IGF-1R and EGFR, respectively: 28 and 25 exhibited high polysomy; 8 and 14 exhibited gene amplification. No tumor was positive for HER2 or c-Met by IHC and BISH. Multivariate analysis revealed that IGF-1R gene amplification (P = 0.027), thymic carcinoma histology, and higher tumor stage were significantly correlated with an adverse prognosis.ConclusionsThymic epithelial tumors frequently express IGF-1R and/or EGFR proteins. IGF-1R gene amplification is suggested to define an unfavorable subset for thymic epithelial tumors.     

High-resolution analysis of genomic copy number alterations in bladder cancer by microarray-based comparative genomic hybridization

We have screened 22 bladder tumour-derived cell lines and one normal urothelium-derived cell line for genome-wide copy number changes using array comparative genomic hybridization (CGH). Comparison of array CGH with existing multiplex-fluorescence in situ hybridization (M-FISH) results revealed excellent concordance. Regions of gain and loss were defined more accurately by array CGH, and several small regions of deletion were detected that were not identified by M-FISH. Numerous genetic changes were identified, many of which were compatible with previous results from conventional CGH and loss of heterozygosity analyses on bladder tumours. The most frequent changes involved complete or partial loss of 4q (83%) and gain of 20q (78%). Other frequent losses were of 18q (65%), 8p (65%), 2q (61%), 6q (61%), 3p (56%), 13q (56%), 4p (52%), 6p (52%), 10p (52%), 10q (52%) and 5p (43%). We have refined the localization of a region of deletion at 8p21.2-p21.3 to an interval of approximately 1 Mb. Five homozygous deletions of tumour suppressor genes were confirmed, and several potentially novel homozygous deletions were identified. In all, 15 high-level amplifications were detected, with a previously reported amplification at 6p22.3 being the most frequent. Real-time PCR analysis revealed a novel candidate gene with consistent overexpression in all cell lines with the 6p22.3 amplicon.     

Clinical significance of high focal adhesion kinase gene copy number and overexpression in invasive breast cancer

Increased expression and activation of focal adhesion kinase (FAK) was reported to be an unfavorable factor in various human cancers, including breast cancer. In order to study FAK gene status as a prognostic factor, we evaluated FAK gene copy number and FAK protein expression in invasive breast cancer. Tumors from 435 patients with invasive breast cancer were evaluated for FAK gene status using fluorescence in situ hybridization (FISH) and immunohistochemistry (IHC) based on tissue array method. Survival analyses were performed using the Kaplan–Meier method in 267 patients. 42 out of 362 evaluable cases (11.6%) showed high polysomy and 22 cases (6.1%) had gene amplification by FISH. 108 out of 393 evaluable cases (27.5%) showed FAK overexpression by IHC. FAK FISH positivity was significantly associated with higher histologic grade, higher T stage, negative estrogen receptor expression, negative progesterone receptor expression and triple-negative phenotype; FAK overexpression with higher histologic grade and triple-negative phenotype. FAK overexpression was noted in 57.8% (37 of 64) of FAK FISH+ cases. The concordance of FISH and IHC results for FAK gene was observed in 74.9% (271 of 362). Survival analyses revealed the patients with FAK FISH+ breast cancer had significantly shorter overall survival and relapse-free survival compared to those with FISH− breast cancer. In node-positive breast cancer patients who received postoperative systemic treatment, the patients with FAK FISH+ showed significantly shorter 5-year survival rates. Despite of high significant concordance between the results of FISH and IHC for FAK gene in invasive breast cancer, only FAK FISH positivity was an unfavorable prognostic factor.     

Germline DNA copy number variation in familial and early-onset breast cancer

Introduction

Genetic factors predisposing individuals to cancer remain elusive in the majority of patients with a familial or clinical history suggestive of hereditary breast cancer. Germline DNA copy number variation (CNV) has recently been implicated in predisposition to cancers such as neuroblastomas as well as prostate and colorectal cancer. We evaluated the role of germline CNVs in breast cancer susceptibility, in particular those with low population frequencies (rare CNVs), which are more likely to cause disease."

Methods

Using whole-genome comparative genomic hybridization on microarrays, we screened a cohort of women fulfilling criteria for hereditary breast cancer who did not carry BRCA1/BRCA2 mutations.

Results

The median numbers of total and rare CNVs per genome were not different between controls and patients. A total of 26 rare germline CNVs were identified in 68 cancer patients, however, a proportion that was significantly different (P = 0.0311) from the control group (23 rare CNVs in 100 individuals). Several of the genes affected by CNV in patients and controls had already been implicated in cancer.

Conclusions

This study is the first to explore the contribution of germline CNVs to BRCA1/2-negative familial and early-onset breast cancer. The data suggest that rare CNVs may contribute to cancer predisposition in this small cohort of patients, and this trend needs to be confirmed in larger population samples.     

Increased copy number at 3p14 in breast cancer

Introduction  

The present study was conducted to investigate if chromosome band 3p14 is of any pathogenic significance in the malignant process of breast cancer. Genetic studies have implicated a tumour suppressor gene on chromosome arm 3p and we have proposed LRIG1 at 3p14 as a candidate tumour suppressor. The LRIG1 gene encodes an integral membrane protein that counteracts signalling by receptor tyrosine kinases belonging to the ERBB family. LRIG1 mRNA and protein are expressed in many tissues, including breast tissue.     

Insulin-like growth factor-1 receptor gene expression is associated with survival in breast cancer: a comprehensive analysis of gene copy number, mRNA and protein expression

Insulin-like growth factor-1 receptor (IGF1R) plays a key role in the initiation and progression of breast cancer. However, its prognostic relevance to breast cancer patients has long been a matter of debate. In a series of 325 primary invasive breast cancer patients, we performed a comprehensive analysis of IGF1R at the levels of gene copy number, mRNA expression and protein expression. The relationship between the IGF1R status and the clinicopathological characteristics and prognosis was evaluated. IGF1R mRNA levels not only correlated with protein expression, but also were significantly associated with several clinicopathological parameters and prognosis. Patients with low nuclear grade, negative axillary lymph nodes, positive hormone receptor, negative Her2, negative Ki67, and luminal subtype tumors showed higher expression levels of IGF1R mRNA, which was shown to be a significant univariate parameter for both relapse-free survival and breast cancer-specific survival (BCSS) as well as a significant multivariate parameter for BCSS. IGF1R protein expression showed an association with a prolonged BCSS in univariate analysis. In contrast, IGF1R gene copy number was not correlated with mRNA and protein expression, and harbored no prognostic value. When studied in the luminal tumor subtype groups, IGF1R mRNA level was still significantly associated with a better BCSS. Overall, our data indicated a correlation between IGF1R mRNA expression and protein expression in primary breast cancer. In particular, IGF1R mRNA expression appeared to be a good prognostic marker both in the entire cohort and in the luminal subtype group. These data may serve as background information for IGF1R-targeted therapy.     

Tree models for dependent copy number changes in bladder cancer

We analyzed comparative genomic hybridization data on a collection of 237 bladder cancer tumors with the aim of identifying sets of copy number aberrations (CNAs) that tend to occur together. A test based on Fisher's exact test for pairs, but taking into account multiple testing, showed strong dependencies amongst several pairs of aberrations including (+1q, -11p), (+17q, +20q), (+10p, -17p), (-8p, -17p), (+5p, +10p). To determine whether co-occurrence of CNAs may characterize tumor subtypes, we used two recently proposed methods to construct tree models of tumor progression. We constructed tree models for all the tumors, the tumors of stage pT1, and the tumors of stages pT2-4. The tree models confirmed that most of the non-random events and the associations are the same for different stages. We conclude that the combination of large data sets and tree models provide a useful approach to systematically identifying tumor subgroups characterized by more than a single chromosomal aberration.