Cathryn M. Gould

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. Clin Cancer Res; 21(14); 3216–29. ©2015 AACR.

A functional siRNA screen identifies genes modulating angiotensin II-mediated EGFR transactivation

Summary The angiotensin type 1 receptor (AT1R) transactivates the epidermal growth factor receptor (EGFR) to mediate cellular growth, however, the molecular mechanisms involved have not yet been resolved. To address this, we performed a functional siRNA screen of the human kinome in human mammary epithelial cells that demonstrate a robust AT1R–EGFR transactivation. We identified a suite of genes encoding proteins that both positively and negatively regulate AT1R–EGFR transactivation. Many candidates are components of EGFR signalling networks, whereas others, including TRIO, BMX and CHKA, have not been previously linked to EGFR transactivation. Individual knockdown of TRIO, BMX or CHKA attenuated tyrosine phosphorylation of the EGFR by angiotensin II stimulation, but this did not occur following direct stimulation of the EGFR with EGF, indicating that these proteins function between the activated AT1R and the EGFR. Further investigation of TRIO and CHKA revealed that their activity is likely to be required for AT1R–EGFR transactivation. CHKA also mediated EGFR transactivation in response to another G protein-coupled receptor (GPCR) ligand, thrombin, indicating a pervasive role for CHKA in GPCR–EGFR crosstalk. Our study reveals the power of unbiased, functional genomic screens to identify new signalling mediators important for tissue remodelling in cardiovascular disease and cancer.

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

BackgroundThe 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.ResultsWe 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.ConclusionsWe 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.

Genome-wide siRNA Screening at Biosafety Level 4 Reveals a Crucial Role for Fibrillarin in Henipavirus Infection.

Hendra and Nipah viruses (genus Henipavirus, family Paramyxoviridae) are highly pathogenic bat-borne viruses. The need for high biocontainment when studying henipaviruses has hindered the development of therapeutics and knowledge of the viral infection cycle. We have performed a genome-wide siRNA screen at biosafety level 4 that identified 585 human proteins required for henipavirus infection. The host protein with the largest impact was fibrillarin, a nucleolar methyltransferase that was also required by measles, mumps and respiratory syncytial viruses for infection. While not required for cell entry, henipavirus RNA and protein syntheses were greatly impaired in cells lacking fibrillarin, indicating a crucial role in the RNA replication phase of infection. During infection, the Hendra virus matrix protein co-localized with fibrillarin in cell nucleoli, and co-associated as a complex in pulldown studies, while its nuclear import was unaffected in fibrillarin-depleted cells. Mutagenesis studies showed that the methyltransferase activity of fibrillarin was required for henipavirus infection, suggesting that this enzyme could be targeted therapeutically to combat henipavirus infections.

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. Mol Cancer Ther; 14(5); 1213–23. ©2015 AACR.

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.

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.Acetylation of the histone variant H2A.Z at gene promoters is associated with oncogene activation; however, it is unclear if such modification has a role in regulating the function of enhancers. Here the authors show that acetylated H2A.Z is redistributed at cancer neo-enhancers and regulates the activity of specific enhancers of cancer-related genes.

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.

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

The signaling proteins and pathways involved in lymphatic endothelial cell migration are identified through a functional screen. Moving lymphatic endothelial cells about Lymphatic vessels return fluid and immune cells from peripheral tissues back to the circulation. The growth of new lymphatic vessels and their remodeling are critical for clearing infection and for metastasis of many cancer subtypes. Williams et al. compared the results of their functional siRNA screens with previously published mRNA data sets to identify genes that regulated lymphatic endothelial cell migration, a process critical for lymphatic vessel growth and remodeling, and genes that functioned in both lymphatic and blood endothelial cell migration. One of the top candidates to emerge from these analyses, the glycan-binding protein Galectin-1 not only promoted lymphatic vessel growth but also was important for maintaining lymphatic endothelial cell identity. Further analyses of the authors’ results may reveal lymphatic vessel–associated proteins that could be targeted to prevent edema, improve infection outcomes, or limit metastasis. 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.