Inhibition of EZH2 transactivation function sensitizes solid tumors to genotoxic stress.

Drugs that block the activity of the methyltransferase EZH2 are in clinical development for the treatment of non-Hodgkin lymphomas harboring EZH2 gain-of-function mutations that enhance its polycomb repressive function. We have previously reported that EZH2 can act as a transcriptional activator in castration-resistant prostate cancer (CRPC). Now we show that EZH2 inhibitors can also block the transactivation activity of EZH2 and inhibit the growth of CRPC cells. Gene expression and epigenomics profiling of cells treated with EZH2 inhibitors demonstrated that in addition to derepressing gene expression, these compounds also robustly down-regulate a set of DNA damage repair (DDR) genes, especially those involved in the base excision repair (BER) pathway. Methylation of the pioneer factor FOXA1 by EZH2 contributes to the activation of these genes, and interaction with the transcriptional coactivator P300 via the transactivation domain on EZH2 directly turns on the transcription. In addition, CRISPR-Cas9-mediated knockout screens in the presence of EZH2 inhibitors identified these BER genes as the determinants that underlie the growth-inhibitory effect of EZH2 inhibitors. Interrogation of public data from diverse types of solid tumors expressing wild-type EZH2 demonstrated that expression of DDR genes is significantly correlated with EZH2 dependency and cellular sensitivity to EZH2 inhibitors. Consistent with these findings, treatment of CRPC cells with EZH2 inhibitors dramatically enhances their sensitivity to genotoxic stress. These studies reveal a previously unappreciated mechanism of action of EZH2 inhibitors and provide a mechanistic basis for potential combination cancer therapies.

Cistrome Explorer: an interactive visual analysis tool for large-scale epigenomic data.

SUMMARY: The regulation of genes by cis-regulatory elements (CREs) is complex and differs between cell types. Visual analysis of large collections of chromatin profiles across diverse cell types, integrated with computational methods, can reveal meaningful biological insights. We developed Cistrome Explorer, a web-based interactive visual analytics tool for exploring thousands of chromatin profiles in diverse cell types. Integrated with the Cistrome Data Browser database which contains thousands of ChIP-seq, DNase-seq and ATAC-seq samples, Cistrome Explorer enables the discovery of patterns of CREs across cell types and the identification of transcription factor binding underlying these patterns. AVAILABILITY AND IMPLEMENTATION: Cistrome Explorer and its source code are available at http://cisvis.gehlenborglab.org/ and released under the MIT License. Documentation can be accessed via http://cisvis.gehlenborglab.org/docs/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Chromosome segregation in Archaea: SegA- and SegB-DNA complex structures provide insights into segrosome assembly.

Genome segregation is a vital process in all organisms. Chromosome partitioning remains obscure in Archaea, the third domain of life. Here, we investigated the SegAB system from Sulfolobus solfataricus. SegA is a ParA Walker-type ATPase and SegB is a site-specific DNA-binding protein. We determined the structures of both proteins and those of SegA-DNA and SegB-DNA complexes. The SegA structure revealed an atypical, novel non-sandwich dimer that binds DNA either in the presence or in the absence of ATP. The SegB structure disclosed a ribbon-helix-helix motif through which the protein binds DNA site specifically. The association of multiple interacting SegB dimers with the DNA results in a higher order chromatin-like structure. The unstructured SegB N-terminus plays an essential catalytic role in stimulating SegA ATPase activity and an architectural regulatory role in segrosome (SegA-SegB-DNA) formation. Electron microscopy results also provide a compact ring-like segrosome structure related to chromosome organization. These findings contribute a novel mechanistic perspective on archaeal chromosome segregation.

Deciphering essential cistromes using genome-wide CRISPR screens.

Although millions of transcription factor binding sites, or cistromes, have been identified across the human genome, defining which of these sites is functional in a given condition remains challenging. Using CRISPR/Cas9 knockout screens and gene essentiality or fitness as the readout, we systematically investigated the essentiality of over 10,000 FOXA1 and CTCF binding sites in breast and prostate cancer cells. We found that essential FOXA1 binding sites act as enhancers to orchestrate the expression of nearby essential genes through the binding of lineage-specific transcription factors. In contrast, CRISPR screens of the CTCF cistrome revealed 2 classes of essential binding sites. The first class of essential CTCF binding sites act like FOXA1 sites as enhancers to regulate the expression of nearby essential genes, while a second class of essential CTCF binding sites was identified at topologically associated domain (TAD) boundaries and display distinct characteristics. Using regression methods trained on our screening data and public epigenetic profiles, we developed a model to predict essential cis-elements with high accuracy. The model for FOXA1 essentiality correctly predicts noncoding variants associated with cancer risk and progression. Taken together, CRISPR screens of cis-regulatory elements can define the essential cistrome of a given factor and can inform the development of predictive models of cistrome function.

Cistrome Data Browser: expanded datasets and new tools for gene regulatory analysis.

The Cistrome Data Browser (DB) is a resource of human and mouse cis-regulatory information derived from ChIP-seq, DNase-seq and ATAC-seq chromatin profiling assays, which map the genome-wide locations of transcription factor binding sites, histone post-translational modifications and regions of chromatin accessible to endonuclease activity. Currently, the Cistrome DB contains approximately 47,000 human and mouse samples with about 24,000 newly collected datasets compared to the previous release two years ago. Furthermore, the Cistrome DB has a new Toolkit module with several features that allow users to better utilize the large-scale ChIP-seq, DNase-seq, and ATAC-seq data. First, users can query the factors which are likely to regulate a specific gene of interest. Second, the Cistrome DB Toolkit facilitates searches for factor binding, histone modifications, and chromatin accessibility in any given genomic interval shorter than 2Mb. Third, the Toolkit can determine the most similar ChIP-seq, DNase-seq, and ATAC-seq samples in terms of genomic interval overlaps with user-provided genomic interval sets. The Cistrome DB is a user-friendly, up-to-date, and well maintained resource, and the new tools will greatly benefit the biomedical research community. The database is freely available at http://cistrome.org/db, and the Toolkit is at http://dbtoolkit.cistrome.org.

Genome-wide CRISPR-based gene knockout screens reveal cellular factors and pathways essential for nasopharyngeal carcinoma.

Early diagnosis of nasopharyngeal carcinoma (NPC) is difficult because of a lack of specific symptoms. Many patients have advanced disease at diagnosis, and these patients respond poorly to treatment. New treatments are therefore needed to improve the outcome of NPC. To better understand the molecular pathogenesis of NPC, here we used an NPC cell line in a genome-wide CRISPR-based knockout screen to identify the cellular factors and pathways essential for NPC (i.e. dependence factors). This screen identified the Moz, Ybf2/Sas3, Sas2, Tip60 histone acetyl transferase complex, NF-κB signaling, purine synthesis, and linear ubiquitination pathways; and MDM2 proto-oncogene as NPC dependence factors/pathways. Using gene knock out, complementary DNA rescue, and inhibitor assays, we found that perturbation of these pathways greatly reduces the growth of NPC cell lines but does not affect growth of SV40-immortalized normal nasopharyngeal epithelial cells. These results suggest that targeting these pathways/proteins may hold promise for achieving better treatment of patients with NPC.

Improved design and analysis of CRISPR knockout screens.

Motivation: Genome-wide clustered, regularly interspaced, short palindromic repeat (CRISPR)-Cas9 screen has been widely used to interrogate gene functions. However, the rules to design better libraries beg further refinement. Results: We found single guide RNA (sgRNA) outliers are characterized by higher G-nucleotide counts, especially in regions distal from the PAM motif and are associated with stronger off-target activities. Furthermore, using non-targeting sgRNAs as negative controls lead to strong bias, which can be mitigated by using sgRNAs targeting multiple 'safe harbor' regions. Custom-designed screens confirmed our findings and further revealed that 19 nt sgRNAs consistently gave the best signal-to-noise ratio. Collectively, our analysis motivated the design of a new genome-wide CRISPR/Cas9 screen library and uncovered some intriguing properties of the CRISPR-Cas9 system. Availability and implementation: The MAGeCK workflow is available open source at https://bitbucket.org/liulab/mageck_nest under the MIT license. Supplementary information: Supplementary data are available at Bioinformatics online.

Sequence determinants of improved CRISPR sgRNA design.

The CRISPR/Cas9 system has revolutionized mammalian somatic cell genetics. Genome-wide functional screens using CRISPR/Cas9-mediated knockout or dCas9 fusion-mediated inhibition/activation (CRISPRi/a) are powerful techniques for discovering phenotype-associated gene function. We systematically assessed the DNA sequence features that contribute to single guide RNA (sgRNA) efficiency in CRISPR-based screens. Leveraging the information from multiple designs, we derived a new sequence model for predicting sgRNA efficiency in CRISPR/Cas9 knockout experiments. Our model confirmed known features and suggested new features including a preference for cytosine at the cleavage site. The model was experimentally validated for sgRNA-mediated mutation rate and protein knockout efficiency. Tested on independent data sets, the model achieved significant results in both positive and negative selection conditions and outperformed existing models. We also found that the sequence preference for CRISPRi/a is substantially different from that for CRISPR/Cas9 knockout and propose a new model for predicting sgRNA efficiency in CRISPRi/a experiments. These results facilitate the genome-wide design of improved sgRNA for both knockout and CRISPRi/a studies.

Epigenetic response to environmental stress: Assembly of BRG1-G9a/GLP-DNMT3 repressive chromatin complex on Myh6 promoter in pathologically stressed hearts.

Chromatin structure is determined by nucleosome positioning, histone modifications, and DNA methylation. How chromatin modifications are coordinately altered under pathological conditions remains elusive. Here we describe a stress-activated mechanism of concerted chromatin modification in the heart. In mice, pathological stress activates cardiomyocytes to express Brg1 (nucleosome-remodeling factor), G9a/Glp (histone methyltransferase), and Dnmt3 (DNA methyltransferase). Once activated, Brg1 recruits G9a and then Dnmt3 to sequentially assemble repressive chromatin-marked by H3K9 and CpG methylation-on a key molecular motor gene (Myh6), thereby silencing Myh6 and impairing cardiac contraction. Disruption of Brg1, G9a or Dnmt3 erases repressive chromatin marks and de-represses Myh6, reducing stress-induced cardiac dysfunction. In human hypertrophic hearts, BRG1-G9a/GLP-DNMT3 complex is also activated; its level correlates with H3K9/CpG methylation, Myh6 repression, and cardiomyopathy. Our studies demonstrate a new mechanism of chromatin assembly in stressed hearts and novel therapeutic targets for restoring Myh6 and ventricular function. The stress-induced Brg1-G9a-Dnmt3 interactions and sequence of repressive chromatin assembly on Myh6 illustrates a molecular mechanism by which the heart epigenetically responds to environmental signals. This article is part of a Special Issue entitled: Cardiomyocyte Biology: Integration of Developmental and Environmental Cues in the Heart edited by Marcus Schaub and Hughes Abriel.

Partitioning the heart: mechanisms of cardiac septation and valve development.

Heart malformations are common congenital defects in humans. Many congenital heart defects involve anomalies in cardiac septation or valve development, and understanding the developmental mechanisms that underlie the formation of cardiac septal and valvular tissues thus has important implications for the diagnosis, prevention and treatment of congenital heart disease. The development of heart septa and valves involves multiple types of progenitor cells that arise either within or outside the heart. Here, we review the morphogenetic events and genetic networks that regulate spatiotemporal interactions between the cells that give rise to septal and valvular tissues and hence partition the heart.

The secondary heart field is a new site of calcineurin/Nfatc1 signaling for semilunar valve development.

Semilunar valve malformations are common human congenital heart defects. Bicuspid aortic valves occur in 2-3% of the population, and pulmonic valve stenosis constitutes 10% of all congenital heart disease in adults (Brickner et al., 2000) [1]. Semilunar valve defects cause valve regurgitation, stenosis, or calcification, leading to endocarditis or congestive heart failure. These complications often require prolonged medical treatment or surgical intervention. Despite the medical importance of valve disease, the regulatory pathways governing semilunar valve development are not entirely clear. In this report we investigated the spatiotemporal role of calcineurin/Nfatc1 signaling in semilunar valve development. We generated conditional knockout mice with calcineurin gene disrupted in various tissues during semilunar valve development. Our studies showed that calcineurin/Nfatc1 pathway signals in the secondary heart field (SHF) but not in the outflow tract myocardium or neural crest cells to regulate semilunar valve morphogenesis. Without SHF calcineurin/Nfatc1 signaling, the conal endocardial cushions-the site of prospective semilunar valve formation--first develop and then regress due to apoptosis, resulting in a striking phenotype with complete absence of the aortic and pulmonic valves, severe valve regurgitation, and perinatal lethality. This role of calcineurin/Nfatc1 signaling in the SHF is different from the requirement of calcineurin/Nfatc1 in the endocardium for semilunar valve formation (Chang et al., 2004) [2], indicating that calcineurin/Nfatc1 signals in multiple tissues to organize semilunar valve development. Also, our studies suggest distinct mechanisms of calcineurin/Nfat signaling for semilunar and atrioventricular valve morphogenesis. Therefore, we demonstrate a novel developmental mechanism in which calcineurin signals through Nfatc1 in the secondary heart field to promote semilunar valve morphogenesis, revealing a new supportive role of the secondary heart field for semilunar valve formation.

Determinants of transcription factor regulatory range.

Characterization of the genomic distances over which transcription factor (TF) binding influences gene expression is important for inferring target genes from TF chromatin immunoprecipitation followed by sequencing (ChIP-seq) data. Here we systematically examine the relationship between thousands of TF and histone modification ChIP-seq data sets with thousands of gene expression profiles. We develop a model for integrating these data, which reveals two classes of TFs with distinct ranges of regulatory influence, chromatin-binding preferences, and auto-regulatory properties. We find that the regulatory range of the same TF bound within different topologically associating domains (TADs) depend on intrinsic TAD properties such as local gene density and G/C content, but also on the TAD chromatin states. Our results suggest that considering TF type, binding distance to gene locus, as well as chromatin context is important in identifying implicated TFs from GWAS SNPs.

CRISPR Screens Identify Essential Cell Growth Mediators in BRAF Inhibitor-resistant Melanoma.

BRAF is a serine/threonine kinase that harbors activating mutations in ∼7% of human malignancies and ∼60% of melanomas. Despite initial clinical responses to BRAF inhibitors, patients frequently develop drug resistance. To identify candidate therapeutic targets for BRAF inhibitor resistant melanoma, we conduct CRISPR screens in melanoma cells harboring an activating BRAF mutation that had also acquired resistance to BRAF inhibitors. To investigate the mechanisms and pathways enabling resistance to BRAF inhibitors in melanomas, we integrate expression, ATAC-seq, and CRISPR screen data. We identify the JUN family transcription factors and the ETS family transcription factor ETV5 as key regulators of CDK6, which together enable resistance to BRAF inhibitors in melanoma cells. Our findings reveal genes contributing to resistance to a selective BRAF inhibitor PLX4720, providing new insights into gene regulation in BRAF inhibitor resistant melanoma cells.

Integrative analysis of pooled CRISPR genetic screens using MAGeCKFlute.

Genome-wide screening using CRISPR coupled with nuclease Cas9 (CRISPR-Cas9) is a powerful technology for the systematic evaluation of gene function. Statistically principled analysis is needed for the accurate identification of gene hits and associated pathways. Here, we describe how to perform computational analysis of CRISPR screens using the MAGeCKFlute pipeline. MAGeCKFlute combines the MAGeCK and MAGeCK-VISPR algorithms and incorporates additional downstream analysis functionalities. MAGeCKFlute is distinguished from other currently available tools by its comprehensive pipeline, which contains a series of functions for analyzing CRISPR screen data. This protocol explains how to use MAGeCKFlute to perform quality control (QC), normalization, batch effect removal, copy-number bias correction, gene hit identification and downstream functional enrichment analysis for CRISPR screens. We also describe gene identification and data analysis in CRISPR screens involving drug treatment. Completing the entire MAGeCKFlute pipeline requires ~3 h on a desktop computer running Linux or Mac OS with R support.

CRISPR-FOCUS: A web server for designing focused CRISPR screening experiments.

The recently developed CRISPR screen technology, based on the CRISPR/Cas9 genome editing system, enables genome-wide interrogation of gene functions in an efficient and cost-effective manner. Although many computational algorithms and web servers have been developed to design single-guide RNAs (sgRNAs) with high specificity and efficiency, algorithms specifically designed for conducting CRISPR screens are still lacking. Here we present CRISPR-FOCUS, a web-based platform to search and prioritize sgRNAs for CRISPR screen experiments. With official gene symbols or RefSeq IDs as the only mandatory input, CRISPR-FOCUS filters and prioritizes sgRNAs based on multiple criteria, including efficiency, specificity, sequence conservation, isoform structure, as well as genomic variations including Single Nucleotide Polymorphisms and cancer somatic mutations. CRISPR-FOCUS also provides pre-defined positive and negative control sgRNAs, as well as other necessary sequences in the construct (e.g., U6 promoters to drive sgRNA transcription and RNA scaffolds of the CRISPR/Cas9). These features allow users to synthesize oligonucleotides directly based on the output of CRISPR-FOCUS. Overall, CRISPR-FOCUS provides a rational and high-throughput approach for sgRNA library design that enables users to efficiently conduct a focused screen experiment targeting up to thousands of genes. (CRISPR-FOCUS is freely available at http://cistrome.org/crispr-focus/).

Integrative Analysis Reveals the Transcriptional Collaboration between EZH2 and E2F1 in the Regulation of Cancer-Related Gene Expression.

UNLABELLED: Overexpression of EZH2 is frequently linked to the advanced and metastatic stage of cancers. The mechanisms of its oncogenic function can be context specific, and may vary depending on the protein complexes that EZH2 interacts with. To identify novel transcriptional collaborators of EZH2 in cancers, a computational approach was developed that integrates protein-DNA binding data, cell perturbation gene expression data, and compendiums of tumor expression profiles. This holistic approach identified E2F1, a known mediator of the Rb tumor suppressor, as a transcriptional collaborator of EZH2 in castration-resistant prostate cancer. Subsequent analysis and experimental validation found EZH2 and E2F1 cobind to a subset of chromatin sites lacking H3K27 trimethylation, and activate genes that are critical for prostate cancer progression. The collaboration of EZH2 and E2F1 in transcriptional regulation is also observed in diffuse large B-cell lymphoma cell lines, where activation of the transcriptional network is concordant with the cellular response to the EZH2 inhibitor. IMPLICATIONS: The direct collaboration between EZH2 and Rb/E2F1 pathway provides an innovative mechanism underlying the cascade of tumor progression, and lays the foundation for the development of new anticancer targets/strategies.

Genome-scale deletion screening of human long non-coding RNAs using a paired-guide RNA CRISPR-Cas9 library.

CRISPR-Cas9 screens have been widely adopted to analyze coding-gene functions, but high-throughput screening of non-coding elements using this method is more challenging because indels caused by a single cut in non-coding regions are unlikely to produce a functional knockout. A high-throughput method to produce deletions of non-coding DNA is needed. We report a high-throughput genomic deletion strategy to screen for functional long non-coding RNAs (lncRNAs) that is based on a lentiviral paired-guide RNA (pgRNA) library. Applying our screening method, we identified 51 lncRNAs that can positively or negatively regulate human cancer cell growth. We validated 9 of 51 lncRNA hits using CRISPR-Cas9-mediated genomic deletion, functional rescue, CRISPR activation or inhibition and gene-expression profiling. Our high-throughput pgRNA genome deletion method will enable rapid identification of functional mammalian non-coding elements.

Quality control, modeling, and visualization of CRISPR screens with MAGeCK-VISPR.

High-throughput CRISPR screens have shown great promise in functional genomics. We present MAGeCK-VISPR, a comprehensive quality control (QC), analysis, and visualization workflow for CRISPR screens. MAGeCK-VISPR defines a set of QC measures to assess the quality of an experiment, and includes a maximum-likelihood algorithm to call essential genes simultaneously under multiple conditions. The algorithm uses a generalized linear model to deconvolute different effects, and employs expectation-maximization to iteratively estimate sgRNA knockout efficiency and gene essentiality. MAGeCK-VISPR also includes VISPR, a framework for the interactive visualization and exploration of QC and analysis results. MAGeCK-VISPR is freely available at http://bitbucket.org/liulab/mageck-vispr .

Tumor size matters differently in pulmonary adenocarcinoma and squamous cell carcinoma.

Little about primary tumor size and nodal/distant metastases among different cell types in non-small cell lung cancer (NSCLC) was discussed. This study aimed to investigate distinct associations between tumor size and nodal/distant metastases in pulmonary adenocarcinoma and squamous cell carcinoma. The study also aimed to clarify the cutoff size relating to a higher likelihood of metastases. We retrospectively evaluated 932 NSCLC patients over a 3-year period and focused on cases with primary tumors less than 4.0 cm in size. Our data showed that 2.5 cm was the critical cutoff size regarding increased nodal/distant metastases in adenocarcinoma (p<0.001), but not in squamous cell carcinoma (p>0.05). In addition, the incidence of nodal/distant metastases reached a plateau of more than 80% in adenocarcinoma when the tumor size exceeded 2.5 cm. In contrast, there was no such correlation observed in squamous cell carcinoma. This study showed that tumor size mattered differently in pulmonary adenocarcinoma and squamous cell carcinoma.