NRG1 fusions in breast cancer.

BACKGROUND: NRG1 gene fusions may be clinically actionable, since cancers carrying the fusion transcripts can be sensitive to tyrosine kinase inhibitors. The NRG1 gene encodes ligands for the HER2(ERBB2)-ERBB3 heterodimeric receptor tyrosine kinase, and the gene fusions are thought to lead to autocrine stimulation of the receptor. The NRG1 fusion expressed in the breast cancer cell line MDA-MB-175 serves as a model example of such fusions, showing the proposed autocrine loop and exceptional drug sensitivity. However, its structure has not been properly characterised, its oncogenic activity has not been fully explained, and there is limited data on such fusions in breast cancer. METHODS: We analysed genomic rearrangements and transcripts of NRG1 in MDA-MB-175 and a panel of 571 breast cancers. RESULTS: We found that the MDA-MB-175 fusion-originally reported as a DOC4(TENM4)-NRG1 fusion, lacking the cytoplasmic tail of NRG1-is in reality a double fusion, PPP6R3-TENM4-NRG1, producing multiple transcripts, some of which include the cytoplasmic tail. We hypothesise that many NRG1 fusions may be oncogenic not for lacking the cytoplasmic domain but because they do not encode NRG1's nuclear-localised form. The fusion in MDA-MB-175 is the result of a very complex genomic rearrangement, which we partially characterised, that creates additional expressed gene fusions, RSF1-TENM4, TPCN2-RSF1, and MRPL48-GAB2. We searched for NRG1 rearrangements in 571 breast cancers subjected to genome sequencing and transcriptome sequencing and found four cases (0.7%) with fusions, WRN-NRG1, FAM91A1-NRG1, ARHGEF39-NRG1, and ZNF704-NRG1, all splicing into NRG1 at the same exon as in MDA-MB-175. However, the WRN-NRG1 and ARHGEF39-NRG1 fusions were out of frame. We identified rearrangements of NRG1 in many more (8% of) cases that seemed more likely to inactivate than to create activating fusions, or whose outcome could not be predicted because they were complex, or both. This is not surprising because NRG1 can be pro-apoptotic and is inactivated in some breast cancers. CONCLUSIONS: Our results highlight the complexity of rearrangements of NRG1 in breast cancers and confirm that some do not activate but inactivate. Careful interpretation of NRG1 rearrangements will therefore be necessary for appropriate patient management.

Structural analysis of the genome of breast cancer cell line ZR-75-30 identifies twelve expressed fusion genes.

BACKGROUND: It has recently emerged that common epithelial cancers such as breast cancers have fusion genes like those in leukaemias. In a representative breast cancer cell line, ZR-75-30, we searched for fusion genes, by analysing genome rearrangements. RESULTS: We first analysed rearrangements of the ZR-75-30 genome, to around 10kb resolution, by molecular cytogenetic approaches, combining array painting and array CGH. We then compared this map with genomic junctions determined by paired-end sequencing. Most of the breakpoints found by array painting and array CGH were identified in the paired end sequencing-55% of the unamplified breakpoints and 97% of the amplified breakpoints (as these are represented by more sequence reads). From this analysis we identified 9 expressed fusion genes: APPBP2-PHF20L1, BCAS3-HOXB9, COL14A1-SKAP1, TAOK1-PCGF2, TIAM1-NRIP1, TIMM23-ARHGAP32, TRPS1-LASP1, USP32-CCDC49 and ZMYM4-OPRD1. We also determined the genomic junctions of a further three expressed fusion genes that had been described by others, BCAS3-ERBB2, DDX5-DEPDC6/DEPTOR and PLEC1-ENPP2. Of this total of 12 expressed fusion genes, 9 were in the coamplification. Due to the sensitivity of the technologies used, we estimate these 12 fusion genes to be around two-thirds of the true total. Many of the fusions seem likely to be driver mutations. For example, PHF20L1, BCAS3, TAOK1, PCGF2, and TRPS1 are fused in other breast cancers. HOXB9 and PHF20L1 are members of gene families that are fused in other neoplasms. Several of the other genes are relevant to cancer-in addition to ERBB2, SKAP1 is an adaptor for Src, DEPTOR regulates the mTOR pathway and NRIP1 is an estrogen-receptor coregulator. CONCLUSIONS: This is the first structural analysis of a breast cancer genome that combines classical molecular cytogenetic approaches with sequencing. Paired-end sequencing was able to detect almost all breakpoints, where there was adequate read depth. It supports the view that gene breakage and gene fusion are important classes of mutation in breast cancer, with a typical breast cancer expressing many fusion genes.

The role of tandem duplicator phenotype in tumour evolution in high-grade serous ovarian cancer.

High-grade serous ovarian carcinoma (HGSOC) is characterized by genomic instability, ubiquitous TP53 loss, and frequent development of platinum resistance. Loss of homologous recombination (HR) is a mutator phenotype present in 50% of HGSOCs and confers hypersensitivity to platinum treatment. We asked which other mutator phenotypes are present in HGSOC and how they drive the emergence of platinum resistance. We performed whole-genome paired-end sequencing on a model of two HGSOC cases, each consisting of a pair of cell lines established before and after clinical resistance emerged, to describe their structural variants (SVs) and to infer their ancestral genomes as the SVs present within each pair. The first case (PEO1/PEO4), with HR deficiency, acquired translocations and small deletions through its early evolution, but a revertant BRCA2 mutation restoring HR function in the resistant lineage re-stabilized its genome and reduced platinum sensitivity. The second case (PEO14/PEO23) had 216 tandem duplications and did not show evidence of HR or mismatch repair deficiency. By comparing the cell lines to the tissues from which they originated, we showed that the tandem duplicator mutator phenotype arose early in progression in vivo and persisted throughout evolution in vivo and in vitro, which may have enabled continual evolution. From the analysis of SNP array data from 454 HGSOC cases in The Cancer Genome Atlas series, we estimate that 12.8% of cases show patterns of aberrations similar to the tandem duplicator, and this phenotype is mutually exclusive with BRCA1/2 carrier mutations.

Single-molecule analysis of genome rearrangements in cancer.

Rearrangements of the genome can be detected by microarray methods and massively parallel sequencing, which identify copy-number alterations and breakpoint junctions, but these techniques are poorly suited to reconstructing the long-range organization of rearranged chromosomes, for example, to distinguish between translocations and insertions. The single-DNA-molecule technique HAPPY mapping is a method for mapping normal genomes that should be able to analyse genome rearrangements, i.e. deviations from a known genome map, to assemble rearrangements into a long-range map. We applied HAPPY mapping to cancer cell lines to show that it could identify rearrangement of genomic segments, even in the presence of normal copies of the genome. We could distinguish a simple interstitial deletion from a copy-number loss at an inversion junction, and detect a known translocation. We could determine whether junctions detected by sequencing were on the same chromosome, by measuring their linkage to each other, and hence map the rearrangement. Finally, we mapped an uncharacterized reciprocal translocation in the T-47D breast cancer cell line to about 2 kb and hence cloned the translocation junctions. We conclude that HAPPY mapping is a versatile tool for determining the structure of rearrangements in the human genome.

Large duplications at reciprocal translocation breakpoints that might be the counterpart of large deletions and could arise from stalled replication bubbles.

Reciprocal chromosome translocations are often not exactly reciprocal. Most familiar are deletions at the breakpoints, up to megabases in extent. We describe here the opposite phenomenon-duplication of tens or hundreds of kilobases at the breakpoint junction, so that the same sequence is present on both products of a translocation. When the products of the translocation are mapped on the genome, they overlap. We report several of these "overlapping-breakpoint" duplications in breast cancer cell lines HCC1187, HCC1806, and DU4475. These lines also had deletions and essentially balanced translocations. In HCC1187 and HCC1806, we identified five cases of duplication ranging between 46 kb and 200 kb, with the partner chromosome showing deletions between 29 bp and 31 Mb. DU4475 had a duplication of at least 200 kb. Breakpoints were mapped using array painting, i.e., hybridization of chromosomes isolated by flow cytometry to custom oligonucleotide microarrays. Duplications were verified by fluorescent in situ hybridization (FISH), PCR on isolated chromosomes, and cloning of breakpoints. We propose that these duplications are the counterpart of deletions and that they are produced at a replication bubble, comprising two replication forks with the duplicated sequence in between. Both copies of the duplicated sequence would go to one daughter cell, on different products of the translocation, while the other daughter cell would show deletion. These duplications may have been overlooked because they may be missed by FISH and array-CGH and may be interpreted as insertions by paired-end sequencing. Such duplications may therefore be quite frequent.

High-resolution analysis of genomic alteration on chromosome arm 8p in urothelial carcinoma.

Loss of chromosome arm 8p, sometimes in combination with amplification of proximal 8p, is found in urothelial carcinoma (UC) and other epithelial cancers and is associated with more advanced tumor stage. We carried out array comparative genomic hybridization on 174 UC and 33 UC cell lines to examine breakpoints and copy number. This was followed by a detailed analysis of the cell lines using fluorescence in situ hybridization (FISH) and, in some cases, M-FISH, to refine breakpoints and determine translocation partners, heterozygosity analysis, and analysis of expression of selected genes. We showed an overall pattern of 8p loss with reduced heterozygosity and reduced gene expression. Amplification was seen in some samples and shown in the cell line JMSU1 to correlate with overexpression of ZNF703, ERLIN2, PROSC, GPR124, and BRF2. Apart from the centromere, no single breakpoint was overrepresented, and we postulate that frequent complex changes without consistent breakpoints reflect the need for alterations of combinations of genes. The region around 2 Mb, which was homozygously deleted in one cell line and includes the gene ARHGEF10 and the micro-RNA hsa-mir-596, is one candidate tumor suppressor gene region.

Progressive 3q amplification consistently targets SOX2 in preinvasive squamous lung cancer.

RATIONALE: Amplification of distal 3q is the most common genomic aberration in squamous lung cancer (SQC). SQC develops in a multistage progression from normal bronchial epithelium through dysplasia to invasive disease. Identifying the key driver events in the early pathogenesis of SQC will facilitate the search for predictive molecular biomarkers and the identification of novel molecular targets for chemoprevention and therapeutic strategies. For technical reasons, previous attempts to analyze 3q amplification in preinvasive lesions have focused on small numbers of predetermined candidate loci rather than an unbiased survey of copy-number variation. OBJECTIVES: To perform a detailed analysis of the 3q amplicon in bronchial dysplasia of different histological grades. METHODS: We use molecular copy-number counting (MCC) to analyze the structure of chromosome 3 in 19 preinvasive bronchial biopsy specimens from 15 patients and sequential biopsy specimens from 3 individuals. MEASUREMENTS AND MAIN RESULTS: We demonstrate that no low-grade lesions, but all high-grade lesions, have 3q amplification. None of seven low-grade lesions progressed clinically, whereas 8 of 10 patients with high-grade disease progressed to cancer. We identify a minimum commonly amplified region on chromosome 3 consisting of 17 genes, including 2 known oncogenes, SOX2 and PIK3CA. We confirm that both genes are amplified in all high-grade dysplastic lesions tested. We further demonstrate, in three individuals, that the clinical progression of high-grade preinvasive disease is associated with incremental amplification of SOX2, suggesting this promotes malignant progression. CONCLUSIONS: These findings demonstrate progressive 3q amplification in the evolution of preinvasive SQC and implicate SOX2 as a key target of this dynamic process.

High-resolution array CGH clarifies events occurring on 8p in carcinogenesis.

BACKGROUND: Rearrangement of the short arm of chromosome 8 (8p) is very common in epithelial cancers such as breast cancer. Usually there is an unbalanced translocation breakpoint in 8p12 (29.7 Mb - 38.5 Mb) with loss of distal 8p, sometimes with proximal amplification of 8p11-12. Rearrangements in 8p11-12 have been investigated using high-resolution array CGH, but the first 30 Mb of 8p are less well characterised, although this region contains several proposed tumour suppressor genes. METHODS: We analysed the whole of 8p by array CGH at tiling-path BAC resolution in 32 breast and six pancreatic cancer cell lines. Regions of recurrent rearrangement distal to 8p12 were further characterised, using regional fosmid arrays. FISH, and quantitative RT-PCR on over 60 breast tumours validated the existence of similar events in primary material. RESULTS: We confirmed that 8p is usually lost up to at least 30 Mb, but a few lines showed focal loss or copy number steps within this region. Three regions showed rearrangements common to at least two cases: two regions of recurrent loss and one region of amplification. Loss within 8p23.3 (0 Mb - 2.2 Mb) was found in six cell lines. Of the genes always affected, ARHGEF10 showed a point mutation of the remaining normal copies in the DU4475 cell line. Deletions within 12.7 Mb - 19.1 Mb in 8p22, in two cases, affected TUSC3. A novel amplicon was found within 8p21.3 (19.1 Mb - 23.4 Mb) in two lines and one of 98 tumours. CONCLUSION: The pattern of rearrangements seen on 8p may be a consequence of the high density of potential targets on this chromosome arm, and ARHGEF10 may be a new candidate tumour suppressor gene.

PMC42, a breast progenitor cancer cell line, has normal-like mRNA and microRNA transcriptomes.

INTRODUCTION: The use of cultured cell lines as model systems for normal tissue is limited by the molecular alterations accompanying the immortalisation process, including changes in the mRNA and microRNA (miRNA) repertoire. Therefore, identification of cell lines with normal-like expression profiles is of paramount importance in studies of normal gene regulation. METHODS: The mRNA and miRNA expression profiles of several breast cell lines of cancerous or normal origin were measured using printed slide arrays, Luminex bead arrays, and real-time reverse transcription-polymerase chain reaction. RESULTS: We demonstrate that the mRNA expression profiles of two breast cell lines are similar to that of normal breast tissue: HB4a, immortalised normal breast epithelium, and PMC42, a breast cancer cell line that retains progenitor pluripotency allowing in-culture differentiation to both secretory and myoepithelial fates. In contrast, only PMC42 exhibits a normal-like miRNA expression profile. We identified a group of miRNAs that are highly expressed in normal breast tissue and PMC42 but are lost in all other cancerous and normal-origin breast cell lines and observed a similar loss in immortalised lymphoblastoid cell lines compared with healthy uncultured B cells. Moreover, like tumour suppressor genes, these miRNAs are lost in a variety of tumours. We show that the mechanism leading to the loss of these miRNAs in breast cancer cell lines has genomic, transcriptional, and post-transcriptional components. CONCLUSION: We propose that, despite its neoplastic origin, PMC42 is an excellent molecular model for normal breast epithelium, providing a unique tool to study breast differentiation and the function of key miRNAs that are typically lost in cancer.

Architectures of somatic genomic rearrangement in human cancer amplicons at sequence-level resolution.

For decades, cytogenetic studies have demonstrated that somatically acquired structural rearrangements of the genome are a common feature of most classes of human cancer. However, the characteristics of these rearrangements at sequence-level resolution have thus far been subject to very limited description. One process that is dependent upon somatic genome rearrangement is gene amplification, a mechanism often exploited by cancer cells to increase copy number and hence expression of dominantly acting cancer genes. The mechanisms underlying gene amplification are complex but must involve chromosome breakage and rejoining. We sequenced 133 different genomic rearrangements identified within four cancer amplicons involving the frequently amplified cancer genes MYC, MYCN, and ERBB2. The observed architectures of rearrangement were diverse and highly distinctive, with evidence for sister chromatid breakage-fusion-bridge cycles, formation and reinsertion of double minutes, and the presence of bizarre clusters of small genomic fragments. There were characteristic features of sequences at the breakage-fusion junctions, indicating roles for nonhomologous end joining and homologous recombination-mediated repair mechanisms together with nontemplated DNA synthesis. Evidence was also found for sequence-dependent variation in susceptibility of the genome to somatic rearrangement. The results therefore provide insights into the DNA breakage and repair processes operative in somatic genome rearrangement and illustrate how the evolutionary histories of individual cancers can be reconstructed from large-scale cancer genome sequencing.

Co-amplification of 8p12 and 11q13 in breast cancers is not the result of a single genomic event.

Epithelial cancers frequently have multiple amplifications, and particular amplicons tend to occur together. These co-amplifications have been suggested to result from amplification of pre-existing junctions between two chromosomes, that is, translocation junctions. We investigated this hypothesis for two amplifications frequent in breast cancer, at 8p12 and 11q13, which had been reported to be associated in Southern blot studies. We confirmed that both genomic amplification and expression of genes was correlated between the frequently-amplified regions of 8p and 11q, in array CGH and microarray expression data, supporting the importance of co-amplification. We examined by FISH the physical structure of co-amplifications that we had identified by array CGH, in five breast cancer cell lines (HCC1500, MDA-MB-134, MDA-MB-175, SUM44, and ZR-75-1), four breast tumors, and a pancreatic cancer cell line (SUIT2). We found a variety of arrangements: amplification of translocation junctions; entirely independent amplification of the two regions on separate chromosomes; and separate amplification of 8p and 11q sequences in distinct sites on the same rearranged chromosome. In this last arrangement, interphase nuclei often showed intermingling of FISH signals from 8p12 and 11q13, giving a false impression that the sequences were interdigitated. We conclude that co-amplification of the main 8p and 11q amplicons in breast tumors is not usually the result of a preceding translocation event but most likely reflects selection of clones that have amplified both loci. This article contains supplementary material available at http://www.interscience.wiley.com/jpages/1045-2257/suppmat.

Genomic analysis of the 8p11-12 amplicon in familial breast cancer.

Amplification of 8p11-12 has been recurrently reported in sporadic breast cancer. These studies define a complex molecular structure with a set of minimal amplified regions, and different putative oncogenes that show a strong correlation between amplification and over-expression such as ZNF703/FLJ14299, SPFH2/C8orf2, BRF2 and RAB11FIP. However, none of these studies were carried out on familial breast malignancies. We have studied the incidence, molecular features and clinical value of this amplification in familial breast tumors associated with BRCA1, BRCA2 and non-BRCA1/2 gene mutations. We detected 9 out of 80 familial tumors with this amplicon by chromosomal comparative genomic hybridization. Next, we used a high-resolution comparative genomic hybridization array covering the 8p11-12 region to characterize this chromosomal region. This approach allowed us to define 2 cores of common amplification that largely overlap with those reported in sporadic tumors. Our findings confirm the molecular complexity of this chromosomal region and indicate that this genomic event is a common alteration in breast cancer, present not only in sporadic but also in familial tumors. Finally, we found correlation between the 8p11-12 amplification and proliferation (Ki-67) and cyclin E expression, which further proves in familial tumors the poor prognosis association previously reported in sporadic breast cancer.

Genomic imbalances and patterns of karyotypic variability in mantle-cell lymphoma cell lines.

Mantle-cell lymphoma (MCL) is genetically characterized by 11q13 chromosomal translocations involving the CCND1 gene. We have characterized five MCL cell lines, JVM-2, GRANTA-519, REC-1, JEKO-1, and NCEB-1, combining metaphase and array comparative genomic hybridization, multicolor-FISH, and molecular analysis. Our results revealed common gained regions at 2p14, 9q31.2-qter, 11q13.1-q21, 13q14-q21.2, 13q34-qter and 18q21.1-q22.1, and losses at 1p21.2-p31.1, 2p11.2, 8p21.2-pter, 9p21.3-pter, 11q23.3-qter, 17p11.2-pter, and 17q21.2-q22.2. All cell lines except JVM-2, displayed moderate or high numerical chromosome instability. In addition, an ongoing level of chromosome rearrangements was observed in REC-1. Surprisingly, NCEB-1 carried several stable mouse chromosomes and showed expression of both human and murine bcl-2 protein. Our findings indicate that these cell lines represent three patterns of chromosome evolution in MCL and may be useful to understand the pathogenesis of this neoplasm.

A 1 Mb minimal amplicon at 8p11-12 in breast cancer identifies new candidate oncogenes.

Amplification of 8p11-12 is a well-known alteration in human breast cancers but the driving oncogene has not been identified. We have developed a high-resolution comparative genomic hybridization array covering 8p11-12 and analysed 33 primary breast tumors, 20 primary ovarian tumors and 27 breast cancer cell lines. Expression analysis of the genes in the region was carried out by using real-time quantitative PCR and/or oligo-microarray profiling. In all, 24% (8/33) of the breast tumors, 5% (1/20) of the ovary tumors and 15% (4/27) of the cell lines showed 8p11-12 amplification. We identified a 1 Mb segment of common amplification that excludes previously proposed candidate genes. Some of the amplified genes did not show overexpression, whereas for others, overexpression was not specifically attributable to amplification. The genes FLJ14299, C8orf2, BRF2 and RAB11FIP, map within the 8p11-12 minimal amplicon, two have a putative function consistent with an oncogenic role, these four genes showed a strong correlation between amplification and overexpression and are therefore the best candidate driver oncogenes at 8p12.

Gene expression changes induced by estrogen and selective estrogen receptor modulators in primary-cultured human endometrial cells: signals that distinguish the human carcinogen tamoxifen.

Tamoxifen has long been the endocrine treatment of choice for women with breast cancer and is now employed for prophylactic use in women at high risk from breast cancer. Other selective estrogen receptor modulators (SERMs), such as raloxifene, mimic some of tamoxifen's beneficial effects and, like tamoxifen, exhibit a complex mixture of organ-specific estrogen agonist and antagonistic properties. However, accompanying the positive effects of tamoxifen has been the emergence of evidence for an increased risk of endometrial cancer associated with its use. A more complete understanding of the mechanism(s) of SERM carcinogenicity and endometrial effects is therefore required. We have sought to compare and characterise the transcript profile of tamoxifen, raloxifene and the agonist estradiol in human endometrial cells. Using primary cultures of human endometria, to best emulate the in vivo responses in a manageable in vitro system, we have shown 230 significant changes in gene expression for epithelial cultures and 83 in stromal cultures, either specific to 17beta-estradiol, tamoxifen or raloxifene, or changed across more than one of the treatments. Considering the transcriptome as a whole, the endometrial responses to raloxifene or tamoxifen were more similar than either drug was to 17beta-estradiol. Treatment of endometrial cultures with tamoxifen resulted in the largest number of gene changes relative to control cultures and a high proportion of genes associated with regulation of gene transcription, cell-cycle control and signal transduction. Tamoxifen-specific changes that might point towards mechanisms for its proliferative response in the endometrium included changes in retinoblastoma and c-myc binding proteins, the APCL, dihydrofolate reductase (DHFR) and E2F1 genes and other transcription factors. Tamoxifen was also found to give rise to the highest number of gene expression changes common to those that characterise malignant endometria. It is anticipated that this study will provide leads for further and more focused investigation into SERM carcinogenicity.

A recurrent chromosome breakpoint in breast cancer at the NRG1/neuregulin 1/heregulin gene.

Most studies of genomic rearrangements in common cancers have focused on regional gains and losses, but some rearrangements may break within specific genes. We previously reported that five breast cancer cell lines have chromosome translocations that break in the NRG1 gene and that could cause abnormal NRG1 expression. NRG1 encodes the Neuregulins 1 (formerly the Heregulins), ligands for members of the ErbB/epidermal growth factor-receptor family, which includes ErbB2/HER2. We have now screened for breaks at NRG1 in paraffin sections of breast tumors. Tissue microarrays were screened by fluorescence in situ hybridization, with hybridization probes proximal and distal to the expected breakpoints. This screen detects breaks but does not distinguish between translocation or deletion breakpoints. The screen was validated with array-comparative genomic hybridization on a custom 8p12 high-density genomic array to detect a lower copy number of the sequences that were lost distal to the breaks. We also precisely mapped the breaks in five tumors with different hybridization probes. Breaks in NRG1 were detected in 6% (19 of 323) of breast cancers and in some lung and ovarian cancers. In an unselected series of 213 cases with follow-up, breast cancers where the break was detected tended to be high-grade (65% grade III compared with 28% of negative cases). They were, like breast tumors in general, mainly ErbB2 low (11 of 13 were low) and estrogen receptor positive (11 of 13 positive).