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.

Interwoven Trimeric Cage-Catenanes with Topological Chirality.

Catenanes have gained increasing attention for their unique features such as topological chirality. To date, the majority of works have focused on catenanes comprising monocyclic rings. Due to the lack of efficient synthetic strategy, catenanes of multiannulated monomers remain scarce. Here, we report the one-pot synthesis of an interwoven trimeric cage-catenane in high yield by dynamic imine condensation between diamine linkers of suitable length and trialdehyde panels in stoichiometry. The formation of cage-catenane is driven by the efficient 6-fold π-π stacking of panels. The monomeric cage and trimeric cage-catenane are interconvertible with reversible imine chemistry, with the latter thermodynamically being more favored. Using a topology-based statistical model, we first reveal that the formation probability of the interwoven catenane surpasses that of its chain-like isomer by 20%. When this pure mathematical model is refined by taking into account the strong template effect provided by the π-π stacking of aromatic panels, it shows that the interwoven structure emerges as the dominant species, almost ruling out the formation of the latter. Although composed of achiral cage monomers, the topological chirality of the interwoven trimeric catenane is unraveled by chiral-high-performance liquid chromatography (HPLC) and circular dichroism (CD) spectroscopy, and single-crystal X-ray diffraction (XRD) analysis of the interwoven cage-catenane also reveals a pair of two topological enantiomers. Our probability analysis-aided rationale would provide a design rationale for guiding the efficient synthesis of topologically sophisticated structures.

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.

Unraveling the Dynamic Molecular Motions of a Twin-Cavity Cage with Slow Configurational but Rapid Conformational Interconversions.

Featuring diverse structural motions/changes, dynamic molecular systems hold promise for executing complex tasks. However, their structural complexity presents formidable challenge in elucidating their kinetics, especially when multiple structural motions are intercorrelated. We herein introduce a twin-cavity cage that features interconvertible C3- and C1-configurations, with each configuration exhibiting interchangeable P- and M-conformations. This molecule is therefore composed of four interconnected chiral species (P)-C3, (M)-C3, (P)-C1, (M)-C1. We showcase an effective approach to decouple these sophisticated structural changes into two kinetically distinct pathways. Utilizing time-dependent 1H NMR spectroscopy at various temperatures, which disregards the transition between mirror-image conformations, we first determine the rate constant (kc) for the C3- to C1-configuration interconversion, while time-dependent circular dichroism spectroscopy at different temperatures quantifies the observed rate constant (kobs) of the ensemble of all the structural changes. As kobs ≫ ${{\rm { \gg }}}$ kc, it allows us to decouple the overall molecular motions into a slow configurational transformation and rapid conformational interconversions, with the latter further dissected into two independent conformational interchanges, namely (P)-C3 ← → ${ \mathbin{{\stackrel{\textstyle\rightarrow} { {\smash{\leftarrow}\vphantom{_{\vbox to.5ex{\vss}}}} } }} }$ (M)-C3 and (P)-C1 ← → ${ \mathbin{{\stackrel{\textstyle\rightarrow} { {\smash{\leftarrow}\vphantom{_{\vbox to.5ex{\vss}}}} } }} }$ (M)-C1. This work, therefore, sheds light on the comprehensive kinetic study of complex molecular dynamics, offering valuable insights for the rational design of smart dynamic materials for applications of sensing, separation, catalysis, molecular machinery, etc.

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.

CRISPR-DO for genome-wide CRISPR design and optimization.

MOTIVATION: Despite the growing popularity in using CRISPR/Cas9 technology for genome editing and gene knockout, its performance still relies on well-designed single guide RNAs (sgRNA). In this study, we propose a web application for the Design and Optimization (CRISPR-DO) of guide sequences that target both coding and non-coding regions in spCas9 CRISPR system across human, mouse, zebrafish, fly and worm genomes. CRISPR-DO uses a computational sequence model to predict sgRNA efficiency, and employs a specificity scoring function to evaluate the potential of off-target effect. It also provides information on functional conservation of target sequences, as well as the overlaps with exons, putative regulatory sequences and single-nucleotide polymorphisms (SNPs). The web application has a user-friendly genome-browser interface to facilitate the selection of the best target DNA sequences for experimental design. AVAILABILITY AND IMPLEMENTATION: CRISPR-DO is available at http://cistrome.org/crispr/ CONTACT: qiliu@tongji.edu.cn or hanxu@jimmy.harvard.edu or xsliu@jimmy.harvard.eduSupplementary information: Supplementary data are available at Bioinformatics online.

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.

Spheroid models to elaborate the broken symmetry and equivalent volume of molecules in crystalline phase.

Dense packing of particles has provided powerful models to elaborate the important structural features of matter in various systems such as liquid, glassy, and crystalline phases. The simplest sphere packing models can represent and capture salient properties of the building blocks for covalent, metallic, and ionic crystals; it, however, becomes insufficient to reflect the broken symmetry of the commonly anisotropic molecules in molecular crystals. Here, we develop spheroid models with a minimal degree of anisotropy, which serve as a simple geometrical representation for a rich spectrum of molecules-including both isotropic and anisotropic, convex and concave ones-in crystalline phases. Our models are determined via an inverse packing approach: Given a molecular crystal, an optimal spheroid model is constructed using a contact diagram, which depicts the packing relationship between neighboring molecules within the crystal. The spheroid models are capable of accurately capturing the broken symmetry and characterizing the equivalent volume of molecules in the crystalline phases. Moreover, our model retrieves such molecular information from low-quality x-ray diffraction data with poorly resolved structures, and by using soft spheroids, it can also describe the packing behavior in cocrystals.

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/).

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.