Nuclear stabilization of p53 requires a functional nucleolar surveillance pathway.

The nucleolar surveillance pathway monitors nucleolar integrity and responds to nucleolar stress by mediating binding of ribosomal proteins to MDM2, resulting in p53 accumulation. Inappropriate pathway activation is implicated in the pathogenesis of ribosomopathies, while drugs selectively activating the pathway are in trials for cancer. Despite this, the molecular mechanism(s) regulating this process are poorly understood. Using genome-wide loss-of-function screens, we demonstrate the ribosome biogenesis axis as the most potent class of genes whose disruption stabilizes p53. Mechanistically, we identify genes critical for regulation of this pathway, including HEATR3. By selectively disabling the nucleolar surveillance pathway, we demonstrate that it is essential for the ability of all nuclear-acting stresses, including DNA damage, to induce p53 accumulation. Our data support a paradigm whereby the nucleolar surveillance pathway is the central integrator of stresses that regulate nuclear p53 abundance, ensuring that ribosome biogenesis is hardwired to cellular proliferative capacity.

DNA methylation is required to maintain both DNA replication timing precision and 3D genome organization integrity.

DNA replication timing and three-dimensional (3D) genome organization are associated with distinct epigenome patterns across large domains. However, whether alterations in the epigenome, in particular cancer-related DNA hypomethylation, affects higher-order levels of genome architecture is still unclear. Here, using Repli-Seq, single-cell Repli-Seq, and Hi-C, we show that genome-wide methylation loss is associated with both concordant loss of replication timing precision and deregulation of 3D genome organization. Notably, we find distinct disruption in 3D genome compartmentalization, striking gains in cell-to-cell replication timing heterogeneity and loss of allelic replication timing in cancer hypomethylation models, potentially through the gene deregulation of DNA replication and genome organization pathways. Finally, we identify ectopic H3K4me3-H3K9me3 domains from across large hypomethylated domains, where late replication is maintained, which we purport serves to protect against catastrophic genome reorganization and aberrant gene transcription. Our results highlight a potential role for the methylome in the maintenance of 3D genome regulation.

BRG1 knockdown inhibits proliferation through multiple cellular pathways in prostate cancer.

BACKGROUND: BRG1 (encoded by SMARCA4) is a catalytic component of the SWI/SNF chromatin remodelling complex, with key roles in modulating DNA accessibility. Dysregulation of BRG1 is observed, but functionally uncharacterised, in a wide range of malignancies. We have probed the functions of BRG1 on a background of prostate cancer to investigate how BRG1 controls gene expression programmes and cancer cell behaviour. RESULTS: Our investigation of SMARCA4 revealed that BRG1 is over-expressed in the majority of the 486 tumours from The Cancer Genome Atlas prostate cohort, as well as in a complementary panel of 21 prostate cell lines. Next, we utilised a temporal model of BRG1 depletion to investigate the molecular effects on global transcription programmes. Depleting BRG1 had no impact on alternative splicing and conferred only modest effect on global expression. However, of the transcriptional changes that occurred, most manifested as down-regulated expression. Deeper examination found the common thread linking down-regulated genes was involvement in proliferation, including several known to increase prostate cancer proliferation (KLK2, PCAT1 and VAV3). Interestingly, the promoters of genes driving proliferation were bound by BRG1 as well as the transcription factors, AR and FOXA1. We also noted that BRG1 depletion repressed genes involved in cell cycle progression and DNA replication, but intriguingly, these pathways operated independently of AR and FOXA1. In agreement with transcriptional changes, depleting BRG1 conferred G1 arrest. CONCLUSIONS: Our data have revealed that BRG1 promotes cell cycle progression and DNA replication, consistent with the increased cell proliferation associated with oncogenesis.

Genome-wide RNAi screen for genes regulating glycolytic response to vemurafenib in BRAFV600 melanoma cells.

Identification of mechanisms underlying sensitivity and response to targeted therapies, such as the BRAF inhibitor vemurafenib, is critical in order to improve efficacy of these therapies in the clinic and delay onset of resistance. Glycolysis has emerged as a key feature of the BRAF inhibitor response in melanoma cells, and importantly, the metabolic response to vemurafenib in melanoma patients can predict patient outcome. Here, we present a multiparameter genome-wide siRNA screening dataset of genes that when depleted improve the viability and glycolytic response to vemurafenib in BRAFV600 mutated melanoma cells. These datasets are suitable for analysis of genes involved in cell viability and glycolysis in steady state conditions and following treatment with vemurafenib, as well as computational approaches to identify gene regulatory networks that mediate response to BRAF inhibition in melanoma.

Genes regulating membrane-associated E-cadherin and proliferation in adenomatous polyposis coli mutant colon cancer cells: High content siRNA screen.

Truncating mutations in the tumour suppressor gene APC occur frequently in colorectal cancers and result in the deregulation of Wnt signalling as well as changes in cell-cell adhesion. Using quantitative imaging based on the detection of membrane-associated E-cadherin, we undertook a protein coding genome-wide siRNA screen to identify genes that regulate cell surface E-cadherin in the APC-defective colorectal cancer cell line SW480. We identified a diverse set of regulators of E-cadherin that offer new insights into the regulation of cell-cell adhesion, junction formation and genes that regulate proliferation or survival of SW480 cells. Among the genes whose depletion promotes membrane-associated E-cadherin, we identified ZEB1, the microRNA200 family, and proteins such as a ubiquitin ligase UBE2E3, CDK8, sorting nexin 27 (SNX27) and the matrix metalloproteinases, MMP14 and MMP19. The screen also identified 167 proteins required for maintaining E-cadherin at cell-cell adherens junctions, including known junctional proteins, CTNND1 and CTNNA1, as well as signalling enzymes, DUSP4 and MARK2, and transcription factors, TEAD3, RUNX2 and TRAM2. A better understanding of the post-translational regulation of E-cadherin provides new opportunities for restoring cell-cell adhesion in APC-defective cells.

Epigenetic reprogramming at estrogen-receptor binding sites alters 3D chromatin landscape in endocrine-resistant breast cancer.

Endocrine therapy resistance frequently develops in estrogen receptor positive (ER+) breast cancer, but the underlying molecular mechanisms are largely unknown. Here, we show that 3-dimensional (3D) chromatin interactions both within and between topologically associating domains (TADs) frequently change in ER+ endocrine-resistant breast cancer cells and that the differential interactions are enriched for resistance-associated genetic variants at CTCF-bound anchors. Ectopic chromatin interactions are preferentially enriched at active enhancers and promoters and ER binding sites, and are associated with altered expression of ER-regulated genes, consistent with dynamic remodelling of ER pathways accompanying the development of endocrine resistance. We observe that loss of 3D chromatin interactions often occurs coincidently with hypermethylation and loss of ER binding. Alterations in active A and inactive B chromosomal compartments are also associated with decreased ER binding and atypical interactions and gene expression. Together, our results suggest that 3D epigenome remodelling is a key mechanism underlying endocrine resistance in ER+ breast cancer.

DNA Hypermethylation Encroachment at CpG Island Borders in Cancer Is Predisposed by H3K4 Monomethylation Patterns.

Promoter CpG islands are typically unmethylated in normal cells, but in cancer a proportion are subject to hypermethylation. Using methylome sequencing we identified CpG islands that display partial methylation encroachment across the 5' or 3' CpG island borders. CpG island methylation encroachment is widespread in prostate and breast cancer and commonly associates with gene suppression. We show that the pattern of H3K4me1 at CpG island borders in normal cells predicts the different modes of cancer CpG island hypermethylation. Notably, genetic manipulation of Kmt2d results in concordant alterations in H3K4me1 levels and CpG island border DNA methylation encroachment. Our findings suggest a role for H3K4me1 in the demarcation of CpG island methylation borders in normal cells, which become eroded in cancer.

Integrated epigenomic analysis stratifies chromatin remodellers into distinct functional groups.

BACKGROUND: ATP-dependent chromatin remodelling complexes are responsible for establishing and maintaining the positions of nucleosomes. Chromatin remodellers are targeted to chromatin by transcription factors and non-coding RNA to remodel the chromatin into functional states. However, the influence of chromatin remodelling on shaping the functional epigenome is not well understood. Moreover, chromatin remodellers have not been extensively explored as a collective group across two-dimensional and three-dimensional epigenomic layers. RESULTS: Here, we have integrated the genome-wide binding profiles of eight chromatin remodellers together with DNA methylation, nucleosome positioning, histone modification and Hi-C chromosomal contacts to reveal that chromatin remodellers can be stratified into two functional groups. Group 1 (BRG1, SNF2H, CHD3 and CHD4) has a clear preference for binding at 'actively marked' chromatin and Group 2 (BRM, INO80, SNF2L and CHD1) for 'repressively marked' chromatin. We find that histone modifications and chromatin architectural features, but not DNA methylation, stratify the remodellers into these functional groups. CONCLUSIONS: Our findings suggest that chromatin remodelling events are synchronous and that chromatin remodellers themselves should be considered simultaneously and not as individual entities in isolation or necessarily by structural similarity, as they are traditionally classified. Their coordinated function should be considered by preference for chromatin features in order to gain a more accurate and comprehensive picture of chromatin regulation.

Replication timing and epigenome remodelling are associated with the nature of chromosomal rearrangements in cancer.

DNA replication timing is known to facilitate the establishment of the epigenome, however, the intimate connection between replication timing and changes to the genome and epigenome in cancer remain largely uncharacterised. Here, we perform Repli-Seq and integrated epigenome analyses and demonstrate that genomic regions that undergo long-range epigenetic deregulation in prostate cancer also show concordant differences in replication timing. A subset of altered replication timing domains are conserved across cancers from different tissue origins. Notably, late-replicating regions in cancer cells display a loss of DNA methylation, and a switch in heterochromatin features from H3K9me3-marked constitutive to H3K27me3-marked facultative heterochromatin. Finally, analysis of 214 prostate and 35 breast cancer genomes reveal that late-replicating regions are prone to cis and early-replication to trans chromosomal rearrangements. Together, our data suggests that the nature of chromosomal rearrangement in cancer is related to the spatial and temporal positioning and altered epigenetic states of early-replicating compared to late-replicating loci.

Discovering cancer vulnerabilities using high-throughput micro-RNA screening.

Micro-RNAs (miRNAs) are potent regulators of gene expression and cellular phenotype. Each miRNA has the potential to target hundreds of transcripts within the cell thus controlling fundamental cellular processes such as survival and proliferation. Here, we exploit this important feature of miRNA networks to discover vulnerabilities in cancer phenotype, and map miRNA-target relationships across different cancer types. More specifically, we report the results of a functional genomics screen of 1280 miRNA mimics and inhibitors in eight cancer cell lines, and its presentation in a sophisticated interactive data portal. This resource represents the most comprehensive survey of miRNA function in oncology, incorporating breast cancer, prostate cancer and neuroblastoma. A user-friendly web portal couples this experimental data with multiple tools for miRNA target prediction, pathway enrichment analysis and visualization. In addition, the database integrates publicly available gene expression and perturbation data enabling tailored and context-specific analysis of miRNA function in a particular disease. As a proof-of-principle, we use the database and its innovative features to uncover novel determinants of the neuroblastoma malignant phenotype.

Acetylated histone variant H2A.Z is involved in the activation of neo-enhancers in prostate cancer.

Acetylation of the histone variant H2A.Z (H2A.Zac) occurs at active promoters and is associated with oncogene activation in prostate cancer, but its role in enhancer function is still poorly understood. Here we show that H2A.Zac containing nucleosomes are commonly redistributed to neo-enhancers in cancer resulting in a concomitant gain of chromatin accessibility and ectopic gene expression. Notably incorporation of acetylated H2A.Z nucleosomes is a pre-requisite for activation of Androgen receptor (AR) associated enhancers. H2A.Zac nucleosome occupancy is rapidly remodeled to flank the AR sites to initiate the formation of nucleosome-free regions and the production of AR-enhancer RNAs upon androgen treatment. Remarkably higher levels of global H2A.Zac correlate with poorer prognosis. Altogether these data demonstrate the novel contribution of H2A.Zac in activation of newly formed enhancers in prostate cancer.

Genome-wide functional analysis reveals central signaling regulators of lymphatic endothelial cell migration and remodeling.

Lymphatic vessels constitute a specialized vasculature that is involved in development, cancer, obesity, and immune regulation. The migration of lymphatic endothelial cells (LECs) is critical for vessel growth (lymphangiogenesis) and vessel remodeling, processes that modify the lymphatic network in response to developmental or pathological demands. Using the publicly accessible results of our genome-wide siRNA screen, we characterized the migratome of primary human LECs and identified individual genes and signaling pathways that regulate LEC migration. We compared our data set with mRNA differential expression data from endothelial and stromal cells derived from two in vivo models of lymphatic vessel remodeling, viral infection and contact hypersensitivity-induced inflammation, which identified genes selectively involved in regulating LEC migration and remodeling. We also characterized the top candidates in the LEC migratome in primary blood vascular endothelial cells to identify genes with functions common to lymphatic and blood vascular endothelium. On the basis of these analyses, we showed that LGALS1, which encodes the glycan-binding protein Galectin-1, promoted lymphatic vascular growth in vitro and in vivo and contributed to maintenance of the lymphatic endothelial phenotype. Our results provide insight into the signaling networks that control lymphangiogenesis and lymphatic remodeling and potentially identify therapeutic targets and biomarkers in disease specific to lymphatic or blood vessels.

Comprehensive evaluation of genome-wide 5-hydroxymethylcytosine profiling approaches in human DNA.

BACKGROUND: The discovery that 5-methylcytosine (5mC) can be oxidized to 5-hydroxymethylcytosine (5hmC) by the ten-eleven translocation (TET) proteins has prompted wide interest in the potential role of 5hmC in reshaping the mammalian DNA methylation landscape. The gold-standard bisulphite conversion technologies to study DNA methylation do not distinguish between 5mC and 5hmC. However, new approaches to mapping 5hmC genome-wide have advanced rapidly, although it is unclear how the different methods compare in accurately calling 5hmC. In this study, we provide a comparative analysis on brain DNA using three 5hmC genome-wide approaches, namely whole-genome bisulphite/oxidative bisulphite sequencing (WG Bis/OxBis-seq), Infinium HumanMethylation450 BeadChip arrays coupled with oxidative bisulphite (HM450K Bis/OxBis) and antibody-based immunoprecipitation and sequencing of hydroxymethylated DNA (hMeDIP-seq). We also perform loci-specific TET-assisted bisulphite sequencing (TAB-seq) for validation of candidate regions. RESULTS: We show that whole-genome single-base resolution approaches are advantaged in providing precise 5hmC values but require high sequencing depth to accurately measure 5hmC, as this modification is commonly in low abundance in mammalian cells. HM450K arrays coupled with oxidative bisulphite provide a cost-effective representation of 5hmC distribution, at CpG sites with 5hmC levels >~10%. However, 5hmC analysis is restricted to the genomic location of the probes, which is an important consideration as 5hmC modification is commonly enriched at enhancer elements. Finally, we show that the widely used hMeDIP-seq method provides an efficient genome-wide profile of 5hmC and shows high correlation with WG Bis/OxBis-seq 5hmC distribution in brain DNA. However, in cell line DNA with low levels of 5hmC, hMeDIP-seq-enriched regions are not detected by WG Bis/OxBis or HM450K, either suggesting misinterpretation of 5hmC calls by hMeDIP or lack of sensitivity of the latter methods. CONCLUSIONS: We highlight both the advantages and caveats of three commonly used genome-wide 5hmC profiling technologies and show that interpretation of 5hmC data can be significantly influenced by the sensitivity of methods used, especially as the levels of 5hmC are low and vary in different cell types and different genomic locations.

Systematic high-content genome-wide RNAi screens of endothelial cell migration and morphology.

Many cell types undergo migration during embryogenesis and disease. Endothelial cells line blood vessels and lymphatics, which migrate during development as part of angiogenesis, lymphangiogenesis and other types of vessel remodelling. These processes are also important in wound healing, cancer metastasis and cardiovascular conditions. However, the molecular control of endothelial cell migration is poorly understood. Here, we present a dataset containing siRNA screens that identify known and novel components of signalling pathways regulating migration of lymphatic endothelial cells. These components are compared to signalling in blood vascular endothelial cells. Further, using high-content microscopy, we captured a dataset of images of migrating cells following transfection with a genome-wide siRNA library. These datasets are suitable for the identification and analysis of genes involved in endothelial cell migration and morphology, and for computational approaches to identify signalling networks controlling the migratory response and integration of cell morphology, gene function and cell signaling. This may facilitate identification of protein targets for therapeutically modulating angiogenesis and lymphangiogenesis in the context of human disease.

Three-dimensional disorganization of the cancer genome occurs coincident with long-range genetic and epigenetic alterations.

A three-dimensional chromatin state underpins the structural and functional basis of the genome by bringing regulatory elements and genes into close spatial proximity to ensure proper, cell-type-specific gene expression profiles. Here, we performed Hi-C chromosome conformation capture sequencing to investigate how three-dimensional chromatin organization is disrupted in the context of copy-number variation, long-range epigenetic remodeling, and atypical gene expression programs in prostate cancer. We find that cancer cells retain the ability to segment their genomes into megabase-sized topologically associated domains (TADs); however, these domains are generally smaller due to establishment of additional domain boundaries. Interestingly, a large proportion of the new cancer-specific domain boundaries occur at regions that display copy-number variation. Notably, a common deletion on 17p13.1 in prostate cancer spanning the TP53 tumor suppressor locus results in bifurcation of a single TAD into two distinct smaller TADs. Change in domain structure is also accompanied by novel cancer-specific chromatin interactions within the TADs that are enriched at regulatory elements such as enhancers, promoters, and insulators, and associated with alterations in gene expression. We also show that differential chromatin interactions across regulatory regions occur within long-range epigenetically activated or silenced regions of concordant gene activation or repression in prostate cancer. Finally, we present a novel visualization tool that enables integrated exploration of Hi-C interaction data, the transcriptome, and epigenome. This study provides new insights into the relationship between long-range epigenetic and genomic dysregulation and changes in higher-order chromatin interactions in cancer.

Systematic Screening Identifies Dual PI3K and mTOR Inhibition as a Conserved Therapeutic Vulnerability in Osteosarcoma.

PURPOSE: Osteosarcoma is the most common cancer of bone occurring mostly in teenagers. Despite rapid advances in our knowledge of the genetics and cell biology of osteosarcoma, significant improvements in patient survival have not been observed. The identification of effective therapeutics has been largely empirically based. The identification of new therapies and therapeutic targets are urgently needed to enable improved outcomes for osteosarcoma patients. EXPERIMENTAL DESIGN: We have used genetically engineered murine models of human osteosarcoma in a systematic, genome-wide screen to identify new candidate therapeutic targets. We performed a genome-wide siRNA screen, with or without doxorubicin. In parallel, a screen of therapeutically relevant small molecules was conducted on primary murine- and primary human osteosarcoma-derived cell cultures. All results were validated across independent cell cultures and across human and mouse osteosarcoma. RESULTS: The results from the genetic and chemical screens significantly overlapped, with a profound enrichment of pathways regulated by PI3K and mTOR pathways. Drugs that concurrently target both PI3K and mTOR were effective at inducing apoptosis in primary osteosarcoma cell cultures in vitro in both human and mouse osteosarcoma, whereas specific PI3K or mTOR inhibitors were not effective. The results were confirmed with siRNA and small molecule approaches. Rationale combinations of specific PI3K and mTOR inhibitors could recapitulate the effect on osteosarcoma cell cultures. CONCLUSIONS: The approaches described here have identified dual inhibition of the PI3K-mTOR pathway as a sensitive, druggable target in osteosarcoma, and provide rationale for translational studies with these agents.

Synthetic Lethal Screens Identify Vulnerabilities in GPCR Signaling and Cytoskeletal Organization in E-Cadherin-Deficient Cells.

The CDH1 gene, which encodes the cell-to-cell adhesion protein E-cadherin, is frequently mutated in lobular breast cancer (LBC) and diffuse gastric cancer (DGC). However, because E-cadherin is a tumor suppressor protein and lost from the cancer cell, it is not a conventional drug target. To overcome this, we have taken a synthetic lethal approach to determine whether the loss of E-cadherin creates druggable vulnerabilities. We first conducted a genome-wide siRNA screen of isogenic MCF10A cells with and without CDH1 expression. Gene ontology analysis demonstrated that G-protein-coupled receptor (GPCR) signaling proteins were highly enriched among the synthetic lethal candidates. Diverse families of cytoskeletal proteins were also frequently represented. These broad classes of E-cadherin synthetic lethal hits were validated using both lentiviral-mediated shRNA knockdown and specific antagonists, including the JAK inhibitor LY2784544, Pertussis toxin, and the aurora kinase inhibitors alisertib and danusertib. Next, we conducted a 4,057 known drug screen and time course studies on the CDH1 isogenic MCF10A cell lines and identified additional drug classes with linkages to GPCR signaling and cytoskeletal function that showed evidence of E-cadherin synthetic lethality. These included multiple histone deacetylase inhibitors, including vorinostat and entinostat, PI3K inhibitors, and the tyrosine kinase inhibitors crizotinib and saracatinib. Together, these results demonstrate that E-cadherin loss creates druggable vulnerabilities that have the potential to improve the management of both sporadic and familial LBC and DGC.

Whole-genome multiparametric screening to identify modulators of epithelial-to-mesenchymal transition.

Metastasis accounts for the poor prognosis of the majority of solid tumors. The phenotypic transition of nonmotile epithelial tumor cells to migratory and invasive "mesenchymal" cells (epithelial-to-mesenchymal transition [EMT]) enables the transit of cancer cells from the primary tumor to distant sites. There is no single marker of EMT; rather, multiple measures are required to define cell state. Thus, the multiparametric capability of high-content screening is ideally suited for the comprehensive analysis of EMT regulators. The aim of this study was to generate a platform to systematically identify functional modulators of tumor cell plasticity using the bladder cancer cell line TSU-Pr1-B1 as a model system. A platform enabling the quantification of key EMT characteristics, cell morphology and mesenchymal intermediate filament vimentin, was developed using the fluorescent whole-cell-tracking reagent CMFDA and a fluorescent promoter reporter construct, respectively. The functional effect of genome-wide modulation of protein-coding genes and miRNAs coupled with those of a collection of small-molecule kinase inhibitors on EMT was assessed using the Target Activation Bioapplication integrated in the Cellomics ArrayScan platform. Data from each of the three screens were integrated to identify a cohort of targets that were subsequently examined in a validation assay using siRNA duplexes. Identification of established regulators of EMT supports the utility of this screening approach and indicated capacity to identify novel regulators of this plasticity program. Pathway analysis coupled with interrogation of cancer-related expression profile databases and other EMT-related screens provided key evidence to prioritize further experimental investigation into the molecular regulators of EMT in cancer cells.

Functional genomics down under: RNAi screening in the Victorian Centre for Functional Genomics.

The Victorian Centre for Functional Genomics (VCFG) is an RNAi screening facility housed at the Peter MacCallum Cancer Centre in Melbourne, Australia. The Peter Mac is Australia's largest dedicated Cancer Research Institute, home to a team of over 520 scientists that focus on understanding the genetic risk of cancer, the molecular events regulating cancer growth and dissemination and improving detection through new diagnostic tools (www.petermac.org). Peter Mac is a well recognised technology leader and established the VCFG with a view to enabling researchers Australia and New Zealand-wide access to cutting edge functional genomics technology, infrastructure and expertise. This review documents the technology platforms operated within the VCFG and provides insight into the workflows and analysis pipelines currently in operation.

Genome-wide functional genomic and transcriptomic analyses for genes regulating sensitivity to vorinostat.

Identification of mechanisms of resistance to histone deacetylase inhibitors, such as vorinostat, is important in order to utilise these anticancer compounds more efficiently in the clinic. Here, we present a dataset containing multiple tiers of stringent siRNA screening for genes that when knocked down conferred sensitivity to vorinostat-induced cell death. We also present data from a miRNA overexpression screen for miRNAs contributing to vorinostat sensitivity. Furthermore, we provide transcriptomic analysis using massively parallel sequencing upon knockdown of 14 validated vorinostat-resistance genes. These datasets are suitable for analysis of genes and miRNAs involved in cell death in the presence and absence of vorinostat as well as computational biology approaches to identify gene regulatory networks.