Evidence that the 5p12 Variant rs10941679 Confers Susceptibility to Estrogen-Receptor-Positive Breast Cancer through FGF10 and MRPS30 Regulation.

Genome-wide association studies (GWASs) have revealed increased breast cancer risk associated with multiple genetic variants at 5p12. Here, we report the fine mapping of this locus using data from 104,660 subjects from 50 case-control studies in the Breast Cancer Association Consortium (BCAC). With data for 3,365 genotyped and imputed SNPs across a 1 Mb region (positions 44,394,495-45,364,167; NCBI build 37), we found evidence for at least three independent signals: the strongest signal, consisting of a single SNP rs10941679, was associated with risk of estrogen-receptor-positive (ER+) breast cancer (per-g allele OR ER+ = 1.15; 95% CI 1.13-1.18; p = 8.35 × 10-30). After adjustment for rs10941679, we detected signal 2, consisting of 38 SNPs more strongly associated with ER-negative (ER-) breast cancer (lead SNP rs6864776: per-a allele OR ER- = 1.10; 95% CI 1.05-1.14; p conditional = 1.44 × 10-12), and a single signal 3 SNP (rs200229088: per-t allele OR ER+ = 1.12; 95% CI 1.09-1.15; p conditional = 1.12 × 10-05). Expression quantitative trait locus analysis in normal breast tissues and breast tumors showed that the g (risk) allele of rs10941679 was associated with increased expression of FGF10 and MRPS30. Functional assays demonstrated that SNP rs10941679 maps to an enhancer element that physically interacts with the FGF10 and MRPS30 promoter regions in breast cancer cell lines. FGF10 is an oncogene that binds to FGFR2 and is overexpressed in ∼10% of human breast cancers, whereas MRPS30 plays a key role in apoptosis. These data suggest that the strongest signal of association at 5p12 is mediated through coordinated activation of FGF10 and MRPS30, two candidate genes for breast cancer pathogenesis.

Polymorphisms in a Putative Enhancer at the 10q21.2 Breast Cancer Risk Locus Regulate NRBF2 Expression.

Genome-wide association studies have identified SNPs near ZNF365 at 10q21.2 that are associated with both breast cancer risk and mammographic density. To identify the most likely causal SNPs, we fine mapped the association signal by genotyping 428 SNPs across the region in 89,050 European and 12,893 Asian case and control subjects from the Breast Cancer Association Consortium. We identified four independent sets of correlated, highly trait-associated variants (iCHAVs), three of which were located within ZNF365. The most strongly risk-associated SNP, rs10995201 in iCHAV1, showed clear evidence of association with both estrogen receptor (ER)-positive (OR = 0.85 [0.82-0.88]) and ER-negative (OR = 0.87 [0.82-0.91]) disease, and was also the SNP most strongly associated with percent mammographic density. iCHAV2 (lead SNP, chr10: 64,258,684:D) and iCHAV3 (lead SNP, rs7922449) were also associated with ER-positive (OR = 0.93 [0.91-0.95] and OR = 1.06 [1.03-1.09]) and ER-negative (OR = 0.95 [0.91-0.98] and OR = 1.08 [1.04-1.13]) disease. There was weaker evidence for iCHAV4, located 5' of ADO, associated only with ER-positive breast cancer (OR = 0.93 [0.90-0.96]). We found 12, 17, 18, and 2 candidate causal SNPs for breast cancer in iCHAVs 1-4, respectively. Chromosome conformation capture analysis showed that iCHAV2 interacts with the ZNF365 and NRBF2 (more than 600 kb away) promoters in normal and cancerous breast epithelial cells. Luciferase assays did not identify SNPs that affect transactivation of ZNF365, but identified a protective haplotype in iCHAV2, associated with silencing of the NRBF2 promoter, implicating this gene in the etiology of breast cancer.

Intrinsic subtypes and benefit from postmastectomy radiotherapy in node-positive premenopausal breast cancer patients who received adjuvant chemotherapy - results from two independent randomized trials.

BACKGROUND: The study of the intrinsic molecular subtypes of breast cancer has revealed differences among them in terms of prognosis and response to chemotherapy and endocrine therapy. However, the ability of intrinsic subtypes to predict benefit from adjuvant radiotherapy has only been examined in few studies. METHODS: Gene expression-based intrinsic subtyping was performed in 228 breast tumors collected from two independent post-mastectomy clinical trials (British Columbia and the Danish Breast Cancer Cooperative Group 82b trials), where pre-menopausal patients with node-positive disease were randomized to adjuvant radiotherapy or not. All patients received adjuvant chemotherapy and a subgroup of patients underwent ovarian ablation. Tumors were classified into intrinsic subtypes: Luminal A, Luminal B, HER2-enriched, Basal-like and Normal-like using the research-based PAM50 classifier. RESULTS: In the British Columbia study, patients treated with radiation had an overall significant lower incidence of locoregional recurrence compared to the controls. For Luminal A tumors the risk of loco-regional recurrence was low and was further lowered by adjuvant radiation. These findings were validated in the DBCG 82b study. The individual data from the two cohorts were merged, the hazard ratio (HR) for loco-regional recurrence associated with giving radiation was 0.34 (0.19 to 0.61) overall and 0.12 (0.03 to 0.52) for Luminal A tumors. CONCLUSIONS: In both postmastectomy trials, patients with Luminal A tumors turned out to have a significant lower incidence of loco-regional recurrence when randomized to adjuvant radiotherapy, leaving no indication to omit postmastectomy adjuvant radiation in pre-menopausal high-risk patients with Luminal A tumors. It was not possible to evaluate the effect of radiotherapy among the other subtypes because of limited sample sizes.

Integrative clustering reveals a novel split in the luminal A subtype of breast cancer with impact on outcome.

BACKGROUND: Breast cancer is a heterogeneous disease at the clinical and molecular level. In this study we integrate classifications extracted from five different molecular levels in order to identify integrated subtypes. METHODS: Tumor tissue from 425 patients with primary breast cancer from the Oslo2 study was cut and blended, and divided into fractions for DNA, RNA and protein isolation and metabolomics, allowing the acquisition of representative and comparable molecular data. Patients were stratified into groups based on their tumor characteristics from five different molecular levels, using various clustering methods. Finally, all previously identified and newly determined subgroups were combined in a multilevel classification using a "cluster-of-clusters" approach with consensus clustering. RESULTS: Based on DNA copy number data, tumors were categorized into three groups according to the complex arm aberration index. mRNA expression profiles divided tumors into five molecular subgroups according to PAM50 subtyping, and clustering based on microRNA expression revealed four subgroups. Reverse-phase protein array data divided tumors into five subgroups. Hierarchical clustering of tumor metabolic profiles revealed three clusters. Combining DNA copy number and mRNA expression classified tumors into seven clusters based on pathway activity levels, and tumors were classified into ten subtypes using integrative clustering. The final consensus clustering that incorporated all aforementioned subtypes revealed six major groups. Five corresponded well with the mRNA subtypes, while a sixth group resulted from a split of the luminal A subtype; these tumors belonged to distinct microRNA clusters. Gain-of-function studies using MCF-7 cells showed that microRNAs differentially expressed between the luminal A clusters were important for cancer cell survival. These microRNAs were used to validate the split in luminal A tumors in four independent breast cancer cohorts. In two cohorts the microRNAs divided tumors into subgroups with significantly different outcomes, and in another a trend was observed. CONCLUSIONS: The six integrated subtypes identified confirm the heterogeneity of breast cancer and show that finer subdivisions of subtypes are evident. Increasing knowledge of the heterogeneity of the luminal A subtype may add pivotal information to guide therapeutic choices, evidently bringing us closer to improved treatment for this largest subgroup of breast cancer.

Assessment of variation in immunosuppressive pathway genes reveals TGFBR2 to be associated with prognosis of estrogen receptor-negative breast cancer after chemotherapy.

INTRODUCTION: Tumor lymphocyte infiltration is associated with clinical response to chemotherapy in estrogen receptor (ER) negative breast cancer. To identify variants in immunosuppressive pathway genes associated with prognosis after adjuvant chemotherapy for ER-negative patients, we studied stage I-III invasive breast cancer patients of European ancestry, including 9,334 ER-positive (3,151 treated with chemotherapy) and 2,334 ER-negative patients (1,499 treated with chemotherapy). METHODS: We pooled data from sixteen studies from the Breast Cancer Association Consortium (BCAC), and employed two independent studies for replications. Overall 3,610 single nucleotide polymorphisms (SNPs) in 133 genes were genotyped as part of the Collaborative Oncological Gene-environment Study, in which phenotype and clinical data were collected and harmonized. Multivariable Cox proportional hazard regression was used to assess genetic associations with overall survival (OS) and breast cancer-specific survival (BCSS). Heterogeneity according to chemotherapy or ER status was evaluated with the log-likelihood ratio test. RESULTS: Three independent SNPs in TGFBR2 and IL12B were associated with OS (P <10⁻³) solely in ER-negative patients after chemotherapy (267 events). Poorer OS associated with TGFBR2 rs1367610 (G > C) (per allele hazard ratio (HR) 1.54 (95% confidence interval (CI) 1.22 to 1.95), P = 3.08 × 10⁻4) was not found in ER-negative patients without chemotherapy or ER-positive patients with chemotherapy (P for interaction <10-3). Two SNPs in IL12B (r² = 0.20) showed different associations with ER-negative disease after chemotherapy: rs2546892 (G > A) with poorer OS (HR 1.50 (95% CI 1.21 to 1.86), P = 1.81 × 10⁻4), and rs2853694 (A > C) with improved OS (HR 0.73 (95% CI 0.61 to 0.87), P = 3.67 × 10⁻4). Similar associations were observed with BCSS. Association with TGFBR2 rs1367610 but not IL12B variants replicated using BCAC Asian samples and the independent Prospective Study of Outcomes in Sporadic versus Hereditary Breast Cancer Study and yielded a combined HR of 1.57 ((95% CI 1.28 to 1.94), P = 2.05 × 10⁻5) without study heterogeneity. CONCLUSIONS: TGFBR2 variants may have prognostic and predictive value in ER-negative breast cancer patients treated with adjuvant chemotherapy. Our findings provide further insights into the development of immunotherapeutic targets for ER-negative breast cancer.

A polymorphism in the base excision repair gene PARP2 is associated with differential prognosis by chemotherapy among postmenopausal breast cancer patients.

BACKGROUND: Personalized therapy considering clinical and genetic patient characteristics will further improve breast cancer survival. Two widely used treatments, chemotherapy and radiotherapy, can induce oxidative DNA damage and, if not repaired, cell death. Since base excision repair (BER) activity is specific for oxidative DNA damage, we hypothesized that germline genetic variation in this pathway will affect breast cancer-specific survival depending on treatment. METHODS: We assessed in 1,408 postmenopausal breast cancer patients from the German MARIE study whether cancer specific survival after adjuvant chemotherapy, anthracycline chemotherapy, and radiotherapy is modulated by 127 Single Nucleotide Polymorphisms (SNPs) in 21 BER genes. For SNPs with interaction terms showing p<0.1 (likelihood ratio test) using multivariable Cox proportional hazard analyses, replication in 6,392 patients from nine studies of the Breast Cancer Association Consortium (BCAC) was performed. RESULTS: rs878156 in PARP2 showed a differential effect by chemotherapy (p=0.093) and was replicated in BCAC studies (p=0.009; combined analysis p=0.002). Compared to non-carriers, carriers of the variant G allele (minor allele frequency=0.07) showed better survival after chemotherapy (combined allelic hazard ratio (HR)=0.75, 95% 0.53-1.07) and poorer survival when not treated with chemotherapy (HR=1.42, 95% 1.08-1.85). A similar effect modification by rs878156 was observed for anthracycline-based chemotherapy in both MARIE and BCAC, with improved survival in carriers (combined allelic HR=0.73, 95% CI 0.40-1.32). None of the SNPs showed significant differential effects by radiotherapy. CONCLUSIONS: Our data suggest for the first time that a SNP in PARP2, rs878156, may together with other genetic variants modulate cancer specific survival in breast cancer patients depending on chemotherapy. These germline SNPs could contribute towards the design of predictive tests for breast cancer patients.

Deregulation of cancer-related miRNAs is a common event in both benign and malignant human breast tumors.

MicroRNAs (miRNAs) are endogenous non-coding RNAs, which play an essential role in the regulation of gene expression during carcinogenesis. The role of miRNAs in breast cancer has been thoroughly investigated, and although many miRNAs are identified as cancer related, little is known about their involvement in benign tumors. In this study, we investigated miRNA expression profiles in the two most common types of human benign tumors (fibroadenoma/fibroadenomatosis) and in malignant breast tumors and explored their role as oncomirs and tumor suppressor miRNAs. Here, we identified 33 miRNAs with similar deregulated expression in both benign and malignant tumors compared with the expression levels of those in normal tissue, including breast cancer-related miRNAs such as let-7, miR-21 and miR-155. Additionally, messenger RNA (mRNA) expression profiles were obtained for some of the same samples. Using integrated mRNA/miRNA expression analysis, we observed that overexpression of certain miRNAs co-occurred with a significant downregulation of their candidate target mRNAs in both benign and malignant tumors. In support of these findings, in vitro functional screening of the downregulated miRNAs in non-malignant and breast cancer cell lines identified several possible tumor suppressor miRNAs, including miR-193b, miR-193a-3p, miR-126, miR-134, miR-132, miR-486-5p, miR-886-3p, miR-195 and miR-497, showing reduced growth when re-expressed in cancer cells. The finding of deregulated expression of oncomirs and tumor suppressor miRNAs in benign breast tumors is intriguing, indicating that they may play a role in proliferation. A role of cancer-related miRNAs in the early phases of carcinogenesis and malignant transformation can, therefore, not be ruled out.

Relationship between the prognostic and predictive value of the intrinsic subtypes and a validated gene profile predictive of loco-regional control and benefit from post-mastectomy radiotherapy in patients with high-risk breast cancer.

BACKGROUND: Breast cancer is characterized by great molecular heterogeneity demonstrated, e.g. by the intrinsic subtypes. Administration of post-mastectomy radiotherapy (PMRT) does, however, not reflect this heterogeneity. A gene profile (DBCG-RT profile) has recently been developed and validated, and has shown prognostic impact in terms of loco-regional failure and predictive impact for PMRT. Reports have also shown predictive value in terms of benefit of PMRT from intrinsic subtypes and derived approximations. The aim of this study was to examine: 1) the agreement between various methods for determining the intrinsic subtypes; and 2) the relationship between the prognostic and predictive impact of the DBCG-RT profile and the intrinsic subtypes. MATERIAL AND METHODS: Intrinsic subtypes and the DBCG-RT profile was determined from microarray analysis based on fresh frozen tissue from 191 patients included in the Danish Breast Cancer Cooperative Group (DBCG) 82bc trial. Corresponding formalin-fixed, paraffin-embedded tissue was available from 146 of these patients and from another 890 DBCG82bc patients. Estrogen receptor, progesterone receptor, HER2, CK5/6, Ki-67 and EGFR were combined into immunohistochemical approximations of the intrinsic subtypes. Endpoint considered was loco-regional recurrence (LRR). RESULTS: The DBCG-RT profile identified a group of patients with low risk of LRR and no additional benefit from PMRT among all subtypes. Combining six immunohistochemical markers identified a subgroup of triple negative patients with high risk of LRR and significant benefit from PMRT. Agreement in the different assignments of tumors to the subtypes was suboptimal, and the clinical outcome and predicted benefit from PMRT varied according to the method used for assignment. CONCLUSION: The prognostic and predictive information obtained from the DBCG-RT profile cannot be substituted by any approximation of the tumors intrinsic subtype. The predictive value of the intrinsic subtypes in terms of PMRT was influenced by the method used for assignment to the intrinsic subtypes.

Large-scale genotyping identifies a new locus at 22q13.2 associated with female breast size.

BACKGROUND: Individual differences in breast size are a conspicuous feature of variation in human females and have been associated with fecundity and advantage in selection of mates. To identify common variants that are associated with breast size, we conducted a large-scale genotyping association meta-analysis in 7169 women of European descent across three independent sample collections with digital or screen film mammograms. METHODS: The samples consisted of the Swedish KARMA, LIBRO-1 and SASBAC studies genotyped on iCOGS, a custom illumina iSelect genotyping array comprising of 211 155 single nucleotide polymorphisms (SNPs) designed for replication and fine mapping of common and rare variants with relevance to breast, ovary and prostate cancer. Breast size of each subject was ascertained by measuring total breast area (mm(2)) on a mammogram. RESULTS: We confirm genome-wide significant associations at 8p11.23 (rs10086016, p=1.3×10(-14)) and report a new locus at 22q13 (rs5995871, p=3.2×10(-8)). The latter region contains the MKL1 gene, which has been shown to impact endogenous oestrogen receptor α transcriptional activity and is recruited on oestradiol sensitive genes. We also replicated previous genome-wide association study findings for breast size at four other loci. CONCLUSIONS: A new locus at 22q13 may be associated with female breast size.

Overrepresentation of transcription factor families in the genesets underlying breast cancer subtypes.

BACKGROUND: The human genome contains a large amount of cis-regulatory DNA elements responsible for directing both spatial and temporal gene-expression patterns. Previous studies have shown that based on their mRNA expression breast tumors could be divided into five subgroups (Luminal A, Luminal B, Basal, ErbB2(+) and Normal-like), each with a distinct molecular portrait. Whole genome gene expression analysis of independent sets of breast tumors reveals repeatedly the robustness of this classification. Furthermore, breast tumors carrying a TP53 mutation show a distinct gene expression profile, which is in strong association to the distinct molecular portraits. The mRNA expression of 552 genes, which varied considerably among the different tumors, but little between two samples of the same tumor, has been shown to be sufficient to separate these tumor subgroups. RESULTS: We analyzed in silico the transcriptional regulation of genes defining the subgroups at 3 different levels: 1. We studied the pathways in which the genes distinguishing the subgroups of breast cancer may be jointly involved including upstream regulators (1st and 2nd level of regulation) as well as downstream targets of these genes. 2. Then we analyzed the promoter areas of these genes (-500 bp tp +100 bp relative to the transcription start site) for canonical transcription binding sites using Genomatix. 3. We looked for the actual expression levels of the identified TF and how they correlate with the overrepresentation of their TF binding sites in the separate groups. We report that promoter composition of the genes that most strongly predict the patient subgroups is distinct. The class-predictive genes showed a clearly different degree of overrepresentation of transcription factor families in their promoter sequences. CONCLUSION: The study suggests that transcription factors responsible for the observed expression pattern in breast cancers may lead us to important biological pathways.

The expression of the long NEAT1_2 isoform is associated with human epidermal growth factor receptor 2-positive breast cancers.

The long non-coding RNA NEAT1 locus is transcribed into two overlapping isoforms, NEAT1_1 and NEAT1_2, of which the latter is essential for the assembly of nuclear paraspeckles. NEAT1 is abnormally expressed in a wide variety of human cancers. Emerging evidence suggests that the two isoforms have distinct functions in gene expression regulation, and recently it was shown that NEAT1_2, but not NEAT1_1, expression predicts poor clinical outcome in cancer. Here, we report that NEAT1_2 expression correlates with HER2-positive breast cancers and high-grade disease. We provide evidence that NEAT1_1 and NEAT1_2 have distinct expression pattern among different intrinsic breast cancer subtypes. Finally, we show that NEAT1_2 expression and paraspeckle formation increase upon lactation in humans, confirming what has previously been demonstrated in mice.

AXL Is a Driver of Stemness in Normal Mammary Gland and Breast Cancer.

The receptor tyrosine kinase AXL is associated with epithelial plasticity in several solid tumors including breast cancer and AXL-targeting agents are currently in clinical trials. We hypothesized that AXL is a driver of stemness traits in cancer by co-option of a regulatory function normally reserved for stem cells. AXL-expressing cells in human mammary epithelial ducts co-expressed markers associated with multipotency, and AXL inhibition abolished colony formation and self-maintenance activities while promoting terminal differentiation in vitro. Axl-null mice did not exhibit a strong developmental phenotype, but enrichment of Axl + cells was required for mouse mammary gland reconstitution upon transplantation, and Axl-null mice had reduced incidence of Wnt1-driven mammary tumors. An AXL-dependent gene signature is a feature of transcriptomes in basal breast cancers and reduced patient survival irrespective of subtype. Our interpretation is that AXL regulates access to epithelial plasticity programs in MaSCs and, when co-opted, maintains acquired stemness in breast cancer cells.

Author Correction: Pan-cancer analysis of homozygous deletions in primary tumours uncovers rare tumour suppressors.

The original version of this Article omitted a declaration from the competing interests statement, which should have included the following: 'K.P.W. is President of Tempus Lab, Inc., Chicago, IL, USA'. This has now been corrected in both the PDF and HTML versions of the Article.

Breast cancer quantitative proteome and proteogenomic landscape.

In the preceding decades, molecular characterization has revolutionized breast cancer (BC) research and therapeutic approaches. Presented herein, an unbiased analysis of breast tumor proteomes, inclusive of 9995 proteins quantified across all tumors, for the first time recapitulates BC subtypes. Additionally, poor-prognosis basal-like and luminal B tumors are further subdivided by immune component infiltration, suggesting the current classification is incomplete. Proteome-based networks distinguish functional protein modules for breast tumor groups, with co-expression of EGFR and MET marking ductal carcinoma in situ regions of normal-like tumors and lending to a more accurate classification of this poorly defined subtype. Genes included within prognostic mRNA panels have significantly higher than average mRNA-protein correlations, and gene copy number alterations are dampened at the protein-level; underscoring the value of proteome quantification for prognostication and phenotypic classification. Furthermore, protein products mapping to non-coding genomic regions are identified; highlighting a potential new class of tumor-specific immunotherapeutic targets.

Basal-like breast cancer engages tumor-supportive macrophages via secreted factors induced by extracellular S100A4.

The tumor microenvironment (TME) may influence both cancer progression and therapeutic response. In breast cancer, particularly in the aggressive triple-negative/basal-like subgroup, patient outcome is strongly associated with the tumor's inflammatory profile. Tumor-associated macrophages (TAMs) are among the most abundant immune cells in the TME, shown to be linked to poor prognosis and therapeutic resistance. In this study, we investigated the effect of the metastasis- and inflammation-associated microenvironmental factor S100A4 on breast cancer cells (BCCs) of different subtypes and explored their further interactions with myeloid cells. We demonstrated that extracellular S100A4 activates BCCs, particularly the basal-like subtype, to elevate secretion of pro-inflammatory cytokines. The secreted factors promoted conversion of monocytes to TAM-like cells that exhibited protumorigenic activities: stimulated epithelial-mesenchymal transition, proliferation, chemoresistance, and motility in cancer cells. In conclusion, we have shown that extracellular S100A4 instigates a tumor-supportive microenvironment, involving a network of cytokines and TAM-like cells, which was particularly characteristic for basal-like BCCs and potentiated their aggressive properties. The S100A4-BCC-TAM interaction cascade could be an important contributor to the aggressive behavior of this subtype and should be further explored for therapeutic targeting.

Author Correction: Pan-cancer analysis of homozygous deletions in primary tumours uncovers rare tumour suppressors.

The original version of this Article contained an error in the author affiliations. The affiliation of Kevin P. White with Tempus Labs, Inc., Chicago, IL, USA was inadvertently omitted.This has now been corrected in both the PDF and HTML versions of the Article.

Time series analysis of neoadjuvant chemotherapy and bevacizumab-treated breast carcinomas reveals a systemic shift in genomic aberrations.

BACKGROUND: Chemotherapeutic agents such as anthracyclines and taxanes are commonly used in the neoadjuvant setting. Bevacizumab is an antibody which binds to vascular endothelial growth factor A (VEGFA) and inhibits its receptor interaction, thus obstructing the formation of new blood vessels. METHODS: A phase II randomized clinical trial of 123 patients with Her2-negative breast cancer was conducted, with patients treated with neoadjuvant chemotherapy (fluorouracil (5FU)/epirubicin/cyclophosphamide (FEC) and taxane), with or without bevacizumab. Serial biopsies were obtained at time of diagnosis, after 12 weeks of treatment with FEC ± bevacizumab, and after 25 weeks of treatment with taxane ± bevacizumab. A time course study was designed to investigate the genomic landscape at the three time points when tumor DNA alterations, tumor percentage, genomic instability, and tumor clonality were assessed. Substantial differences were observed with some tumors changing mainly between diagnosis and at 12 weeks, others between 12 and 25 weeks, and still others changing in both time periods. RESULTS: In both treatment arms, good responders (GR) and non-responders (NR) displayed significant difference in genomic instability index (GII) at time of diagnosis. In the combination arm, copy number alterations at 25 loci at the time of diagnosis were significantly different between the GR and NR. An inverse aberration pattern was also observed between the two extreme response groups at 6p22-p12 for patients in the combination arm. Signs of subclonal reduction were observed, with some aberrations disappearing and others being retained during treatment. Increase in subclonal amplification was observed at 6p21.1, a locus which contains the VEGFA gene for the protein which are targeted by the study drug bevacizumab. Of the 13 pre-treatment samples that had a gain at VEGFA, 12 were responders. Significant decrease of frequency of subclones carrying gains at 17q21.32-q22 was observed at 12 weeks, with the peak occurring at TMEM100, an ALK1 receptor signaling-dependent gene essential for vasculogenesis. This implies that cells bearing amplifications of VEGFA and TMEM100 are particularly sensitive to this treatment regime. CONCLUSIONS: Taken together, these results suggest that heterogeneity and subclonal architecture influence the response to targeted treatment in combination with chemotherapy, with possible implications for clinical decision-making and monitoring of treatment efficacy. TRIAL REGISTRATION: NCT00773695 . Registered 15 October 2008.

Publisher Correction: Evidence that breast cancer risk at the 2q35 locus is mediated through IGFBP5 regulation.

This corrects the article DOI: 10.1038/ncomms5999.

Pan-cancer analysis of homozygous deletions in primary tumours uncovers rare tumour suppressors.

Homozygous deletions are rare in cancers and often target tumour suppressor genes. Here, we build a compendium of 2218 primary tumours across 12 human cancer types and systematically screen for homozygous deletions, aiming to identify rare tumour suppressors. Our analysis defines 96 genomic regions recurrently targeted by homozygous deletions. These recurrent homozygous deletions occur either over tumour suppressors or over fragile sites, regions of increased genomic instability. We construct a statistical model that separates fragile sites from regions showing signatures of positive selection for homozygous deletions and identify candidate tumour suppressors within those regions. We find 16 established tumour suppressors and propose 27 candidate tumour suppressors. Several of these genes (including MGMT, RAD17, and USP44) show prior evidence of a tumour suppressive function. Other candidate tumour suppressors, such as MAFTRR, KIAA1551, and IGF2BP2, are novel. Our study demonstrates how rare tumour suppressors can be identified through copy number meta-analysis.

The Longitudinal Transcriptional Response to Neoadjuvant Chemotherapy with and without Bevacizumab in Breast Cancer.

Purpose: Chemotherapy-induced alterations to gene expression are due to transcriptional reprogramming of tumor cells or subclonal adaptations to treatment. The effect on whole-transcriptome mRNA expression was investigated in a randomized phase II clinical trial to assess the effect of neoadjuvant chemotherapy with the addition of bevacizumab.Experimental Design: Tumor biopsies and whole-transcriptome mRNA profiles were obtained at three fixed time points with 66 patients in each arm. Altogether, 358 specimens from 132 patients were available, representing the transcriptional state before treatment start, at 12 weeks and after treatment (25 weeks). Pathologic complete response (pCR) in breast and axillary nodes was the primary endpoint.Results: pCR was observed in 15 patients (23%) receiving bevacizumab and chemotherapy and 8 patients (12%) receiving only chemotherapy. In the estrogen receptor-positive patients, 11 of 54 (20%) treated with bevacizumab and chemotherapy achieved pCR, while only 3 of 57 (5%) treated with chemotherapy reached pCR. In patients with estrogen receptor-positive tumors treated with combination therapy, an elevated immune activity was associated with good response. Proliferation was reduced after treatment in both treatment arms and most pronounced in the combination therapy arm, where the reduction in proliferation accelerated during treatment. Transcriptional alterations during therapy were subtype specific, and the effect of adding bevacizumab was most evident for luminal-B tumors.Conclusions: Clinical response and gene expression response differed between patients receiving combination therapy and chemotherapy alone. The results may guide identification of patients likely to benefit from antiangiogenic therapy. Clin Cancer Res; 23(16); 4662-70. ©2017 AACR.