Nanoencapsulation of MDM2 Inhibitor RG7388 and Class-I HDAC Inhibitor Entinostat Enhances their Therapeutic Potential Through Synergistic Antitumor Effects and Reduction of Systemic Toxicity.

Inhibitors of the p53-MDM2 interaction such as RG7388 have been developed to exploit latent tumor suppressive properties in p53 in 50% of tumors in which p53 is wild-type. However, these agents for the most part activate cell cycle arrest rather than death, and high doses in patients elicit on-target dose-limiting neutropenia. Recent work from our group indicates that combination of p53-MDM2 inhibitors with the class-I HDAC inhibitor Entinostat (which itself has dose-limiting toxicity issues) has the potential to significantly augment cell death in p53 wild-type colorectal cancer cells. We investigated whether coencapsulation of RG7388 and Entinostat within polymeric nanoparticles (NPs) could overcome efficacy and toxicity limitations of this drug combination. Combinations of RG7388 and Entinostat across a range of different molar ratios resulted in synergistic increases in cell death when delivered in both free drug and nanoencapsulated formats in all colorectal cell lines tested. Importantly, we also explored the in vivo impact of the drug combination on murine blood leukocytes, showing that the leukopenia induced by the free drugs could be significantly mitigated by nanoencapsulation. Taken together, this study demonstrates that formulating these agents within a single nanoparticle delivery platform may provide clinical utility beyond use as nonencapsulated agents.

SynLeGG: analysis and visualization of multiomics data for discovery of cancer 'Achilles Heels' and gene function relationships.

Achilles' heel relationships arise when the status of one gene exposes a cell's vulnerability to perturbation of a second gene, such as chemical inhibition, providing therapeutic opportunities for precision oncology. SynLeGG (www.overton-lab.uk/synlegg) identifies and visualizes mutually exclusive loss signatures in 'omics data to enable discovery of genetic dependency relationships (GDRs) across 783 cancer cell lines and 30 tissues. While there is significant focus on genetic approaches, transcriptome data has advantages for investigation of GDRs and remains relatively underexplored. SynLeGG depends upon the MultiSEp algorithm for unsupervised assignment of cell lines into gene expression clusters, which provide the basis for analysis of CRISPR scores and mutational status in order to propose candidate GDRs. Benchmarking against SynLethDB demonstrates favourable performance for MultiSEp against competing approaches, finding significantly higher area under the Receiver Operator Characteristic curve and between 2.8-fold to 8.5-fold greater coverage. In addition to pan-cancer analysis, SynLeGG offers investigation of tissue-specific GDRs and recovers established relationships, including synthetic lethality for SMARCA2 with SMARCA4. Proteomics, Gene Ontology, protein-protein interactions and paralogue information are provided to assist interpretation and candidate drug target prioritization. SynLeGG predictions are significantly enriched in dependencies validated by a recently published CRISPR screen.

The pseudo-caspase FLIP(L) regulates cell fate following p53 activation.

p53 is the most frequently mutated, well-studied tumor-suppressor gene, yet the molecular basis of the switch from p53-induced cell-cycle arrest to apoptosis remains poorly understood. Using a combination of transcriptomics and functional genomics, we unexpectedly identified a nodal role for the caspase-8 paralog and only human pseudo-caspase, FLIP(L), in regulating this switch. Moreover, we identify FLIP(L) as a direct p53 transcriptional target gene that is rapidly up-regulated in response to Nutlin-3A, an MDM2 inhibitor that potently activates p53. Genetically or pharmacologically inhibiting expression of FLIP(L) using siRNA or entinostat (a clinically relevant class-I HDAC inhibitor) efficiently promoted apoptosis in colorectal cancer cells in response to Nutlin-3A, which otherwise predominantly induced cell-cycle arrest. Enhanced apoptosis was also observed when entinostat was combined with clinically relevant, p53-activating chemotherapy in vitro, and this translated into enhanced in vivo efficacy. Mechanistically, FLIP(L) inhibited p53-induced apoptosis by blocking activation of caspase-8 by the TRAIL-R2/DR5 death receptor; notably, this activation was not dependent on receptor engagement by its ligand, TRAIL. In the absence of caspase-8, another of its paralogs, caspase-10 (also transcriptionally up-regulated by p53), induced apoptosis in Nutlin-3A-treated, FLIP(L)-depleted cells, albeit to a lesser extent than in caspase-8-proficient cells. FLIP(L) depletion also modulated transcription of canonical p53 target genes, suppressing p53-induced expression of the cell-cycle regulator p21 and enhancing p53-induced up-regulation of proapoptotic PUMA. Thus, even in the absence of caspase-8/10, FLIP(L) silencing promoted p53-induced apoptosis by enhancing PUMA expression. Thus, we report unexpected, therapeutically relevant roles for FLIP(L) in determining cell fate following p53 activation.

A chromatin-independent role of Polycomb-like 1 to stabilize p53 and promote cellular quiescence.

Polycomb-like proteins 1-3 (PCL1-3) are substoichiometric components of the Polycomb-repressive complex 2 (PRC2) that are essential for association of the complex with chromatin. However, it remains unclear why three proteins with such apparent functional redundancy exist in mammals. Here we characterize their divergent roles in both positively and negatively regulating cellular proliferation. We show that while PCL2 and PCL3 are E2F-regulated genes expressed in proliferating cells, PCL1 is a p53 target gene predominantly expressed in quiescent cells. Ectopic expression of any PCL protein recruits PRC2 to repress the INK4A gene; however, only PCL2 and PCL3 confer an INK4A-dependent proliferative advantage. Remarkably, PCL1 has evolved a PRC2- and chromatin-independent function to negatively regulate proliferation. We show that PCL1 binds to and stabilizes p53 to induce cellular quiescence. Moreover, depletion of PCL1 phenocopies the defects in maintaining cellular quiescence associated with p53 loss. This newly evolved function is achieved by the binding of the PCL1 N-terminal PHD domain to the C-terminal domain of p53 through two unique serine residues, which were acquired during recent vertebrate evolution. This study illustrates the functional bifurcation of PCL proteins, which act in both a chromatin-dependent and a chromatin-independent manner to regulate the INK4A and p53 pathways.

Activation of innate-adaptive immune machinery by poly(I:C) exposes a therapeutic vulnerability to prevent relapse in stroma-rich colon cancer.

OBJECTIVE: Stroma-rich tumours represent a poor prognostic subtype in stage II/III colon cancer (CC), with high relapse rates and limited response to standard adjuvant chemotherapy. DESIGN: To address the lack of efficacious therapeutic options for patients with stroma-rich CC, we stratified our human tumour cohorts according to stromal content, enabling identification of the biology underpinning relapse and potential therapeutic vulnerabilities specifically within stroma-rich tumours that could be exploited clinically. Following human tumour-based discovery and independent clinical validation, we use a series of in vitro and stroma-rich in vivo models to test and validate the therapeutic potential of elevating the biology associated with reduced relapse in human tumours. RESULTS: By performing our analyses specifically within the stroma-rich/high-fibroblast (HiFi) subtype of CC, we identify and validate the clinical value of a HiFi-specific prognostic signature (HPS), which stratifies tumours based on STAT1-related signalling (High-HPS v Low-HPS=HR 0.093, CI 0.019 to 0.466). Using in silico, in vitro and in vivo models, we demonstrate that the HPS is associated with antigen processing and presentation within discrete immune lineages in stroma-rich CC, downstream of double-stranded RNA and viral response signalling. Treatment with the TLR3 agonist poly(I:C) elevated the HPS signalling and antigen processing phenotype across in vitro and in vivo models. In an in vivo model of stroma-rich CC, poly(I:C) treatment significantly increased systemic cytotoxic T cell activity (p<0.05) and reduced liver metastases (p<0.0002). CONCLUSION: This study reveals new biological insight that offers a novel therapeutic option to reduce relapse rates in patients with the worst prognosis CC.

classifieR a flexible interactive cloud-application for functional annotation of cancer transcriptomes.

BACKGROUND: Transcriptionally informed predictions are increasingly important for sub-typing cancer patients, understanding underlying biology and to inform novel treatment strategies. For instance, colorectal cancers (CRCs) can be classified into four CRC consensus molecular subgroups (CMS) or five intrinsic (CRIS) sub-types that have prognostic and predictive value. Breast cancer (BRCA) has five PAM50 molecular subgroups with similar value, and the OncotypeDX test provides transcriptomic based clinically actionable treatment-risk stratification. However, assigning samples to these subtypes and other transcriptionally inferred predictions is time consuming and requires significant bioinformatics experience. There is no "universal" method of using data from diverse assay/sequencing platforms to provide subgroup classification using the established classifier sets of genes (CMS, CRIS, PAM50, OncotypeDX), nor one which in provides additional useful functional annotations such as cellular composition, single-sample Gene Set Enrichment Analysis, or prediction of transcription factor activity. RESULTS: To address this bottleneck, we developed classifieR, an easy-to-use R-Shiny based web application that supports flexible rapid single sample annotation of transcriptional profiles derived from cancer patient samples form diverse platforms. We demonstrate the utility of the " classifieR" framework to applications focused on the analysis of transcriptional profiles from colorectal (classifieRc) and breast (classifieRb). Samples are annotated with disease relevant transcriptional subgroups (CMS/CRIS sub-types in classifieRc and PAM50/inferred OncotypeDX in classifieRb), estimation of cellular composition using MCP-counter and xCell, single-sample Gene Set Enrichment Analysis (ssGSEA) and transcription factor activity predictions with Discriminant Regulon Expression Analysis (DoRothEA). CONCLUSIONS: classifieR provides a framework which enables labs without access to a dedicated bioinformation can get information on the molecular makeup of their samples, providing an insight into patient prognosis, druggability and also as a tool for analysis and discovery. Applications are hosted online at https://generatr.qub.ac.uk/app/classifieRc and https://generatr.qub.ac.uk/app/classifieRb after signing up for an account on https://generatr.qub.ac.uk .

A revised model of TRAIL-R2 DISC assembly explains how FLIP(L) can inhibit or promote apoptosis.

The long FLIP splice form FLIP(L) can act as both an inhibitor and promoter of caspase-8 at death-inducing signalling complexes (DISCs) formed by death receptors such as TRAIL-R2 and related intracellular complexes such as the ripoptosome. Herein, we describe a revised DISC assembly model that explains how FLIP(L) can have these opposite effects by defining the stoichiometry (with respect to caspase-8) at which it converts from being anti- to pro-apoptotic at the DISC. We also show that in the complete absence of FLIP(L), procaspase-8 activation at the TRAIL-R2 DISC has significantly slower kinetics, although ultimately the extent of apoptosis is significantly greater. This revised model of DISC assembly also explains why FLIP's recruitment to the TRAIL-R2 DISC is impaired in the absence of caspase-8 despite showing that it can interact with the DISC adaptor protein FADD and why the short FLIP splice form FLIP(S) is the more potent inhibitor of DISC-mediated apoptosis.

Identification of a BRCA1-mRNA splicing complex required for efficient DNA repair and maintenance of genomic stability.

Mutations within BRCA1 predispose carriers to a high risk of breast and ovarian cancers. BRCA1 functions to maintain genomic stability through the assembly of multiple protein complexes involved in DNA repair, cell-cycle arrest, and transcriptional regulation. Here, we report the identification of a DNA damage-induced BRCA1 protein complex containing BCLAF1 and other key components of the mRNA-splicing machinery. In response to DNA damage, this complex regulates pre-mRNA splicing of a number of genes involved in DNA damage signaling and repair, thereby promoting the stability of these transcripts/proteins. Further, we show that abrogation of this complex results in sensitivity to DNA damage, defective DNA repair, and genomic instability. Interestingly, mutations in a number of proteins found within this complex have been identified in numerous cancer types. These data suggest that regulation of splicing by the BRCA1-mRNA splicing complex plays an important role in the cellular response to DNA damage.

Tumor suppressor p53: from engaging DNA to target gene regulation.

The p53 transcription factor confers its potent tumor suppressor functions primarily through the regulation of a large network of target genes. The recent explosion of next generation sequencing protocols has enabled the study of the p53 gene regulatory network (GRN) and underlying mechanisms at an unprecedented depth and scale, helping us to understand precisely how p53 controls gene regulation. Here, we discuss our current understanding of where and how p53 binds to DNA and chromatin, its pioneer-like role, and how this affects gene regulation. We provide an overview of the p53 GRN and the direct and indirect mechanisms through which p53 affects gene regulation. In particular, we focus on delineating the ubiquitous and cell type-specific network of regulatory elements that p53 engages; reviewing our understanding of how, where, and when p53 binds to DNA and the mechanisms through which these events regulate transcription. Finally, we discuss the evolution of the p53 GRN and how recent work has revealed remarkable differences between vertebrates, which are of particular importance to cancer researchers using mouse models.

Chronic loss of STAG2 leads to altered chromatin structure contributing to de-regulated transcription in AML.

BACKGROUND: The cohesin complex plays a major role in folding the human genome into 3D structural domains. Mutations in members of the cohesin complex are known early drivers of myelodysplastic syndromes (MDS) and acute myeloid leukaemia (AML), with STAG2 the most frequently mutated complex member. METHODS: Here we use functional genomics (RNA-seq, ChIP-seq and HiChIP) to investigate the impact of chronic STAG2 loss on three-dimensional genome structure and transcriptional programming in a clinically relevant model of chronic STAG2 loss. RESULTS: The chronic loss of STAG2 led to loss of smaller loop domains and the maintenance/formation of large domains that, in turn, led to altered genome compartmentalisation. These changes in genome structure resulted in altered gene expression, including deregulation of the HOXA locus and the MAPK signalling pathway, resulting in increased sensitivity to MEK inhibition. CONCLUSIONS: The altered genomic architecture driven by the chronic loss of STAG2 results in altered gene expression that may contribute to leukaemogenesis and may be therapeutically targeted.

Recommendations for determining HPV status in patients with oropharyngeal cancers under TNM8 guidelines: a two-tier approach.

BACKGROUND: TNM8 staging for oropharyngeal squamous cell carcinomas (OPSCC) surrogates p16 immunohistochemistry for HPV testing. Patients with p16+ OPSCC may lack HPV aetiology. Here, we evaluate the suitability of TNM8 staging for guiding prognosis in such patients. METHODS: HPV status was ascertained using p16 immunohistochemistry and high-risk HPV RNA and DNA in situ hybridisation. Survival by stage in a cohort of OPSCC patients was evaluated using TNM7/TNM8 staging. Survival of p16+/HPV- patients was compared to p16 status. RESULTS: TNM8 staging was found to improve on TNM7 (log rank p = 0·0190 for TNM8 compared with p = 0·0530 for TNM7) in p16+ patients. Patients who tested p16+ but were HPV- (n = 20) had significantly reduced five-year survival (33%) compared to p16+ patients (77%) but not p16- patients (35%). Cancer stage was reduced in 95% of p16+/HPV- patients despite having a mortality rate twice (HR 2.66 [95% CI: 1.37-5.15]) that of p16+/HPV+ patients under new TNM8 staging criteria. CONCLUSION: Given the significantly poorer survival of p16+/HPV- OPSCCs, these data provide compelling evidence for use of an HPV-specific test for staging classification. This has particular relevance in light of potential treatment de-escalation that could expose these patients to inappropriately reduced treatment intensity as treatment algorithms evolve.

HPV16 Down-Regulates the Insulin-Like Growth Factor Binding Protein 2 to Promote Epithelial Invasion in Organotypic Cultures.

Cervical cancer is a multi-stage disease caused by human papillomaviruses (HPV) infection of cervical epithelial cells, but the mechanisms regulating disease progression are not clearly defined. Using 3-dimensional organotypic cultures, we demonstrate that HPV16 E6 and E7 proteins alter the secretome of primary human keratinocytes resulting in local epithelial invasion. Mechanistically, absence of the IGF-binding protein 2 (IGFBP2) caused increases in IGFI/II signalling and through crosstalk with KGF/FGFR2b/AKT, cell invasion. Repression of IGFBP2 is mediated by histone deacetylation at the IGFBP2 promoter and was reversed by treatment with histone deacetylase (HDAC) inhibitors. Our in vitro findings were confirmed in 50 invasive cancers and 79 cervical intra-epithelial neoplastic lesions caused by HPV16 infection, where IGFBP2 levels were reduced with increasing disease severity. In summary, the loss of IGFBP2 is associated with progression of premalignant disease, and sensitises cells to pro-invasive IGF signalling, and together with stromal derived factors promotes epithelial invasion.

Genome-wide characterization reveals complex interplay between TP53 and TP63 in response to genotoxic stress.

In response to genotoxic stress the TP53 tumour suppressor activates target gene expression to induce cell cycle arrest or apoptosis depending on the extent of DNA damage. These canonical activities can be repressed by TP63 in normal stratifying epithelia to maintain proliferative capacity or drive proliferation of squamous cell carcinomas, where TP63 is frequently overexpressed/amplified. Here we use ChIP-sequencing, integrated with microarray analysis, to define the genome-wide interplay between TP53 and TP63 in response to genotoxic stress in normal cells. We reveal that TP53 and TP63 bind to overlapping, but distinct cistromes of sites through utilization of distinctive consensus motifs and that TP53 is constitutively bound to a number of sites. We demonstrate that cisplatin and adriamycin elicit distinct effects on TP53 and TP63 binding events, through which TP53 can induce or repress transcription of an extensive network of genes by direct binding and/or modulation of TP63 activity. Collectively, this results in a global TP53-dependent repression of cell cycle progression, mitosis and DNA damage repair concomitant with activation of anti-proliferative and pro-apoptotic canonical target genes. Further analyses reveal that in the absence of genotoxic stress TP63 plays an important role in maintaining expression of DNA repair genes, loss of which results in defective repair.

Functional and genetic analysis of the colon cancer network.

Cancer is a complex disease that has proven to be difficult to understand on the single-gene level. For this reason a functional elucidation needs to take interactions among genes on a systems-level into account. In this study, we infer a colon cancer network from a large-scale gene expression data set by using the method BC3Net. We provide a structural and a functional analysis of this network and also connect its molecular interaction structure with the chromosomal locations of the genes enabling the definition of cis- and trans-interactions. Furthermore, we investigate the interaction of genes that can be found in close neighborhoods on the chromosomes to gain insight into regulatory mechanisms. To our knowledge this is the first study analyzing the genome-scale colon cancer network.

Genome-wide analysis of p63 binding sites identifies AP-2 factors as co-regulators of epidermal differentiation.

The p63 transcription factor (TP63) is critical in development, growth and differentiation of stratifying epithelia. This is highlighted by the severity of congenital abnormalities caused by TP63 mutations in humans, the dramatic phenotypes in knockout mice and de-regulation of TP63 expression in neoplasia altering the tumour suppressive roles of the TP53 family. In order to define the normal role played by TP63 and provide the basis for better understanding how this network is perturbed in disease, we used chromatin immunoprecipitation combined with massively parallel sequencing (ChIP-seq) to identify >7500 high-confidence TP63-binding regions across the entire genome, in primary human neonatal foreskin keratinocytes (HFKs). Using integrative strategies, we demonstrate that only a subset of these sites are bound by TP53 in response to DNA damage. We identify a role for TP63 in transcriptional regulation of multiple genes genetically linked to cleft palate and identify AP-2alpha (TFAP2A) as a co-regulator of a subset of these genes. We further demonstrate that AP-2gamma (TFAP2C) can bind a subset of these regions and that acute depletion of either TFAP2A or TFAP2C alone is sufficient to reduce terminal differentiation of organotypic epidermal skin equivalents, indicating overlapping physiological functions with TP63.

Dual-Hit Strategy for Therapeutic Targeting of Pancreatic Cancer in Patient-Derived Xenograft Tumors.

PURPOSE: Paracrine activation of pro-fibrotic hedgehog (HH) signaling in pancreatic ductal adenocarcinoma (PDAC) results in stromal amplification that compromises tumor drug delivery, efficacy, and patient survival. Interdiction of HH-mediated tumor-stroma crosstalk with smoothened (SMO) inhibitors (SHHi) "primes" PDAC patient-derived xenograft (PDX) tumors for increased drug delivery by transiently increasing vascular patency/permeability, and thereby macromolecule delivery. However, patient tumor isolates vary in their responsiveness, and responders show co-induction of epithelial-mesenchymal transition (EMT). We aimed to identify the signal derangements responsible for EMT induction and reverse them and devise approaches to stratify SHHi-responsive tumors noninvasively based on clinically-quantifiable parameters. EXPERIMENTAL DESIGN: Animals underwent diffusion-weighted magnetic resonance (DW-MR) imaging for measurement of intratumor diffusivity. In parallel, tissue-level deposition of nanoparticle probes was quantified as a marker of vascular permeability/perfusion. Transcriptomic and bioinformatic analysis was employed to investigate SHHi-induced gene reprogramming and identify key "nodes" responsible for EMT induction. RESULTS: Multiple patient tumor isolates responded to short-term SHH inhibitor exposure with increased vascular patency and permeability, with proportionate increases in tumor diffusivity. Nonresponding PDXs did not. SHHi-treated tumors showed elevated FGF drive and distinctly higher nuclear localization of fibroblast growth factor receptor (FGFR1) in EMT-polarized tumor cells. Pan-FGFR inhibitor NVP-BGJ398 (Infigratinib) reversed the SHHi-induced EMT marker expression and nuclear FGFR1 accumulation without compromising the enhanced permeability effect. CONCLUSIONS: This dual-hit strategy of SMO and FGFR inhibition provides a clinically-translatable approach to compromise the profound impermeability of PDAC tumors. Furthermore, clinical deployment of DW-MR imaging could fulfill the essential clinical-translational requirement for patient stratification.

Spatial Effects of Infiltrating T cells on Neighbouring Cancer Cells and Prognosis in Stage III CRC patients.

Colorectal cancer (CRC) is one of the most frequently occurring cancers, but prognostic biomarkers identifying patients at risk of recurrence are still lacking. In this study, we aimed to investigate in more detail the spatial relationship between intratumoural T cells, cancer cells, and cancer cell hallmarks, as prognostic biomarkers in stage III colorectal cancer patients. We conducted multiplexed imaging of 56 protein markers at single cell resolution on resected fixed tissue from stage III CRC patients who received adjuvant 5-fluorouracil-based chemotherapy. Images underwent segmentation for tumour, stroma and immune cells, and cancer cell 'state' protein marker expression was quantified at a cellular level. We developed a Python package for estimation of spatial proximity, nearest neighbour analysis focusing on cancer cell - T cell interactions at single-cell level. In our discovery cohort (MSK), we processed 462 core samples (total number of cells: 1,669,228) from 221 adjuvant 5FU-treated stage III patients. The validation cohort (HV) consisted of 272 samples (total number of cells: 853,398) from 98 stage III CRC patients. While there were trends for an association between percentage of cytotoxic T cells (across the whole cancer core), it did not reach significance (Discovery cohort: p = 0.07, Validation cohort: p = 0.19). We next utilized our region-based nearest neighbourhood approach to determine the spatial relationships between cytotoxic T cells, helper T cells and cancer cell clusters. In the both cohorts, we found that lower distance between cytotoxic T cells, T helper cells and cancer cells was significantly associated with increased disease-free survival. An unsupervised trained model that clustered patients based on the median distance between immune cells and cancer cells, as well as protein expression profiles, successfully classified patients into low-risk and high-risk groups (Discovery cohort: p = 0.01, Validation cohort: p = 0.003).

p63 maintains keratinocyte proliferative capacity through regulation of Skp2-p130 levels.

p63 is a master regulator of proliferation and differentiation in stratifying epithelia, and its expression is frequently altered in carcinogenesis. However, its role in maintaining proliferative capacity remains unclear. Here, we demonstrate that hypoproliferation and loss of differentiation in organotypic raft cultures of primary neonatal human foreskin keratinocytes (HFKs) depleted of the α and ß isoforms of p63 result from p53-p21-mediated accumulation of retinoblastoma (Rb) family member p130. Hypoproliferation in p63-depleted HFKs can be rescued by depletion of p53, p21(CIP1) or p130. Furthermore, we identified the gene encoding S-phase kinase-associated protein 2 (Skp2), the recognition component of the SCF(Skp2) E3 ubiquitin ligase, as a novel target of p63, potentially influencing p130 levels. Expression of Skp2 is maintained by p63 binding to a site in intron 2 and mRNA levels are downregulated in p63-depleted cells. Hypoproliferation in p63-depleted cells can be restored by re-expression of Skp2. Taken together, these results indicate that p63 plays a multifaceted role in maintaining proliferation in the mature regenerating epidermis, in addition to being required for differentiation.

A whole-genome massively parallel sequencing analysis of BRCA1 mutant oestrogen receptor-negative and -positive breast cancers.

BRCA1 encodes a tumour suppressor protein that plays pivotal roles in homologous recombination (HR) DNA repair, cell-cycle checkpoints, and transcriptional regulation. BRCA1 germline mutations confer a high risk of early-onset breast and ovarian cancer. In more than 80% of cases, tumours arising in BRCA1 germline mutation carriers are oestrogen receptor (ER)-negative; however, up to 15% are ER-positive. It has been suggested that BRCA1 ER-positive breast cancers constitute sporadic cancers arising in the context of a BRCA1 germline mutation rather than being causally related to BRCA1 loss-of-function. Whole-genome massively parallel sequencing of ER-positive and ER-negative BRCA1 breast cancers, and their respective germline DNAs, was used to characterize the genetic landscape of BRCA1 cancers at base-pair resolution. Only BRCA1 germline mutations, somatic loss of the wild-type allele, and TP53 somatic mutations were recurrently found in the index cases. BRCA1 breast cancers displayed a mutational signature consistent with that caused by lack of HR DNA repair in both ER-positive and ER-negative cases. Sequencing analysis of independent cohorts of hereditary BRCA1 and sporadic non-BRCA1 breast cancers for the presence of recurrent pathogenic mutations and/or homozygous deletions found in the index cases revealed that DAPK3, TMEM135, KIAA1797, PDE4D, and GATA4 are potential additional drivers of breast cancers. This study demonstrates that BRCA1 pathogenic germline mutations coupled with somatic loss of the wild-type allele are not sufficient for hereditary breast cancers to display an ER-negative phenotype, and has led to the identification of three potential novel breast cancer genes (ie DAPK3, TMEM135, and GATA4).

Profiling of the BRCA1 transcriptome through microarray and ChIP-chip analysis.

A role for BRCA1 in the direct and indirect regulation of transcription is well established. However, a comprehensive view of the degree to which BRCA1 impacts transcriptional regulation on a genome-wide level has not been defined. We performed genome-wide expression profiling and ChIP-chip analysis, comparison of which revealed that although BRCA1 depletion results in transcriptional changes in 1294 genes, only 44 of these are promoter bound by BRCA1. However, 27% of these transcripts were linked to transcriptional regulation possibly explaining the large number of indirect transcriptional changes observed by microarray analysis. We show that no specific consensus sequence exists for BRCA1 DNA binding but rather demonstrate the presence of a number of known and novel transcription factor (TF)- binding sites commonly found on BRCA1 bound promoters. Co-immunoprecipitations confirmed that BRCA1 interacts with a number of these TFs including AP2-α, PAX2 and ZF5. Finally, we show that BRCA1 is bound to a subset of promoters of genes that are not altered by BRCA1 loss, but are transcriptionally regulated in a BRCA1-dependent manner upon DNA damage. These data suggest a model, whereby BRCA1 is present on defined promoters as part of an inactive complex poised to respond to various genotoxic stimuli.