Determination and inference of eukaryotic transcription factor sequence specificity.

Transcription factor (TF) DNA sequence preferences direct their regulatory activity, but are currently known for only ∼1% of eukaryotic TFs. Broadly sampling DNA-binding domain (DBD) types from multiple eukaryotic clades, we determined DNA sequence preferences for >1,000 TFs encompassing 54 different DBD classes from 131 diverse eukaryotes. We find that closely related DBDs almost always have very similar DNA sequence preferences, enabling inference of motifs for ∼34% of the ∼170,000 known or predicted eukaryotic TFs. Sequences matching both measured and inferred motifs are enriched in chromatin immunoprecipitation sequencing (ChIP-seq) peaks and upstream of transcription start sites in diverse eukaryotic lineages. SNPs defining expression quantitative trait loci in Arabidopsis promoters are also enriched for predicted TF binding sites. Importantly, our motif "library" can be used to identify specific TFs whose binding may be altered by human disease risk alleles. These data present a powerful resource for mapping transcriptional networks across eukaryotes.

DNA Repair Profiling Reveals Nonrandom Outcomes at Cas9-Mediated Breaks.

The repair outcomes at site-specific DNA double-strand breaks (DSBs) generated by the RNA-guided DNA endonuclease Cas9 determine how gene function is altered. Despite the widespread adoption of CRISPR-Cas9 technology to induce DSBs for genome engineering, the resulting repair products have not been examined in depth. Here, the DNA repair profiles of 223 sites in the human genome demonstrate that the pattern of DNA repair following Cas9 cutting at each site is nonrandom and consistent across experimental replicates, cell lines, and reagent delivery methods. Furthermore, the repair outcomes are determined by the protospacer sequence rather than genomic context, indicating that DNA repair profiling in cell lines can be used to anticipate repair outcomes in primary cells. Chemical inhibition of DNA-PK enabled dissection of the DNA repair profiles into contributions from c-NHEJ and MMEJ. Finally, this work elucidates a strategy for using "error-prone" DNA-repair machinery to generate precise edits.

Author Correction: CPSF3-dependent pre-mRNA processing as a druggable node in AML and Ewing's sarcoma.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

CPSF3-dependent pre-mRNA processing as a druggable node in AML and Ewing's sarcoma.

The post-genomic era has seen many advances in our understanding of cancer pathways, yet resistance and tumor heterogeneity necessitate multiple approaches to target even monogenic tumors. Here, we combine phenotypic screening with chemical genetics to identify pre-messenger RNA endonuclease cleavage and polyadenylation specificity factor 3 (CPSF3) as the target of JTE-607, a small molecule with previously unknown target. We show that CPSF3 represents a synthetic lethal node in a subset of acute myeloid leukemia (AML) and Ewing's sarcoma cancer cell lines. Inhibition of CPSF3 by JTE-607 alters expression of known downstream effectors in AML and Ewing's sarcoma lines, upregulates apoptosis and causes tumor-selective stasis in mouse xenografts. Mechanistically, it prevents the release of newly synthesized pre-mRNAs, resulting in read-through transcription and the formation of DNA-RNA hybrid R-loop structures. This study implicates pre-mRNA processing, and specifically CPSF3, as a druggable target providing an avenue to therapeutic intervention in cancer.

Extensive rewiring and complex evolutionary dynamics in a C. elegans multiparameter transcription factor network.

Gene duplication results in two identical paralogs that diverge through mutation, leading to loss or gain of interactions with other biomolecules. Here, we comprehensively characterize such network rewiring for C. elegans transcription factors (TFs) within and across four newly delineated molecular networks. Remarkably, we find that even highly similar TFs often have different interaction degrees and partners. In addition, we find that most TF families have a member that is highly connected in multiple networks. Further, different TF families have opposing correlations between network connectivity and phylogenetic age, suggesting that they are subject to different evolutionary pressures. Finally, TFs that have similar partners in one network generally do not in another, indicating a lack of pressure to retain cross-network similarity. Our multiparameter analyses provide unique insights into the evolutionary dynamics that shaped TF networks.

Genome-wide CRISPR screen for PARKIN regulators reveals transcriptional repression as a determinant of mitophagy.

PARKIN, an E3 ligase mutated in familial Parkinson's disease, promotes mitophagy by ubiquitinating mitochondrial proteins for efficient engagement of the autophagy machinery. Specifically, PARKIN-synthesized ubiquitin chains represent targets for the PINK1 kinase generating phosphoS65-ubiquitin (pUb), which constitutes the mitophagy signal. Physiological regulation of PARKIN abundance, however, and the impact on pUb accumulation are poorly understood. Using cells designed to discover physiological regulators of PARKIN abundance, we performed a pooled genome-wide CRISPR/Cas9 knockout screen. Testing identified genes individually resulted in a list of 53 positive and negative regulators. A transcriptional repressor network including THAP11 was identified and negatively regulates endogenous PARKIN abundance. RNAseq analysis revealed the PARKIN-encoding locus as a prime THAP11 target, and THAP11 CRISPR knockout in multiple cell types enhanced pUb accumulation. Thus, our work demonstrates the critical role of PARKIN abundance, identifies regulating genes, and reveals a link between transcriptional repression and mitophagy, which is also apparent in human induced pluripotent stem cell-derived neurons, a disease-relevant cell type.

TMEM41B is a novel regulator of autophagy and lipid mobilization.

Autophagy maintains cellular homeostasis by targeting damaged organelles, pathogens, or misfolded protein aggregates for lysosomal degradation. The autophagic process is initiated by the formation of autophagosomes, which can selectively enclose cargo via autophagy cargo receptors. A machinery of well-characterized autophagy-related proteins orchestrates the biogenesis of autophagosomes; however, the origin of the required membranes is incompletely understood. Here, we have applied sensitized pooled CRISPR screens and identify the uncharacterized transmembrane protein TMEM41B as a novel regulator of autophagy. In the absence of TMEM41B, autophagosome biogenesis is stalled, LC3 accumulates at WIPI2- and DFCP1-positive isolation membranes, and lysosomal flux of autophagy cargo receptors and intracellular bacteria is impaired. In addition to defective autophagy, TMEM41B knockout cells display significantly enlarged lipid droplets and reduced mobilization and ß-oxidation of fatty acids. Immunostaining and interaction proteomics data suggest that TMEM41B localizes to the endoplasmic reticulum (ER). Taken together, we propose that TMEM41B is a novel ER-localized regulator of autophagosome biogenesis and lipid mobilization.

A gene-centered C. elegans protein-DNA interaction network provides a framework for functional predictions.

Transcription factors (TFs) play a central role in controlling spatiotemporal gene expression and the response to environmental cues. A comprehensive understanding of gene regulation requires integrating physical protein-DNA interactions (PDIs) with TF regulatory activity, expression patterns, and phenotypic data. Although great progress has been made in mapping PDIs using chromatin immunoprecipitation, these studies have only characterized ~10% of TFs in any metazoan species. The nematode C. elegans has been widely used to study gene regulation due to its compact genome with short regulatory sequences. Here, we delineated the largest gene-centered metazoan PDI network to date by examining interactions between 90% of C. elegans TFs and 15% of gene promoters. We used this network as a backbone to predict TF binding sites for 77 TFs, two-thirds of which are novel, as well as integrate gene expression, protein-protein interaction, and phenotypic data to predict regulatory and biological functions for multiple genes and TFs.

Transcription factor binding to Caenorhabditis elegans first introns reveals lack of redundancy with gene promoters.

Gene expression is controlled through the binding of transcription factors (TFs) to regulatory genomic regions. First introns are longer than other introns in multiple eukaryotic species and are under selective constraint. Here we explore the importance of first introns in TF binding in the nematode Caenorhabditis elegans by combining computational predictions and experimentally derived TF-DNA interaction data. We found that first introns of C. elegans genes, particularly those for families enriched in long first introns, are more conserved in length, have more conserved predicted TF interactions and are bound by more TFs than other introns. We detected a significant positive correlation between first intron size and the number of TF interactions obtained from chromatin immunoprecipitation assays or determined by yeast one-hybrid assays. TFs that bind first introns are largely different from those binding promoters, suggesting that the different interactions are complementary rather than redundant. By combining first intron and promoter interactions, we found that genes that share a large fraction of TF interactions are more likely to be co-expressed than when only TF interactions with promoters are considered. Altogether, our data suggest that C. elegans gene regulation may be additive through the combined effects of multiple regulatory regions.

Multiple transcription factors directly regulate Hox gene lin-39 expression in ventral hypodermal cells of the C. elegans embryo and larva, including the hypodermal fate regulators LIN-26 and ELT-6.

BACKGROUND: Hox genes encode master regulators of regional fate specification during early metazoan development. Much is known about the initiation and regulation of Hox gene expression in Drosophila and vertebrates, but less is known in the non-arthropod invertebrate model system, C. elegans. The C. elegans Hox gene lin-39 is required for correct fate specification in the midbody region, including the Vulval Precursor Cells (VPCs). To better understand lin-39 regulation and function, we aimed to identify transcription factors necessary for lin-39 expression in the VPCs, and in particular sought factors that initiate lin-39 expression in the embryo. RESULTS: We used the yeast one-hybrid (Y1H) method to screen for factors that bound to 13 fragments from the lin-39 region: twelve fragments contained sequences conserved between C. elegans and two other nematode species, while one fragment was known to drive reporter gene expression in the early embryo in cells that generate the VPCs. Sixteen transcription factors that bind to eight lin-39 genomic fragments were identified in yeast, and we characterized several factors by verifying their physical interactions in vitro, and showing that reduction of their function leads to alterations in lin-39 levels and lin-39::GFP reporter expression in vivo. Three factors, the orphan nuclear hormone receptor NHR-43, the hypodermal fate regulator LIN-26, and the GATA factor ELT-6 positively regulate lin-39 expression in the embryonic precursors to the VPCs. In particular, ELT-6 interacts with an enhancer that drives GFP expression in the early embryo, and the ELT-6 site we identified is necessary for proper embryonic expression. These three factors, along with the factors ZTF-17, BED-3 and TBX-9, also positively regulate lin-39 expression in the larval VPCs. CONCLUSIONS: These results significantly expand the number of factors known to directly bind and regulate lin-39 expression, identify the first factors required for lin-39 expression in the embryo, and hint at a positive feedback mechanism involving GATA factors that maintains lin-39 expression in the vulval lineage. This work indicates that, as in other organisms, the regulation of Hox gene expression in C. elegans is complicated, redundant and robust.

Enhanced yeast one-hybrid assays for high-throughput gene-centered regulatory network mapping.

A major challenge in systems biology is to understand the gene regulatory networks that drive development, physiology and pathology. Interactions between transcription factors and regulatory genomic regions provide the first level of gene control. Gateway-compatible yeast one-hybrid (Y1H) assays present a convenient method to identify and characterize the repertoire of transcription factors that can bind a DNA sequence of interest. To delineate genome-scale regulatory networks, however, large sets of DNA fragments need to be processed at high throughput and high coverage. Here we present enhanced Y1H (eY1H) assays that use a robotic mating platform with a set of improved Y1H reagents and automated readout quantification. We demonstrate that eY1H assays provide excellent coverage and identify interacting transcription factors for multiple DNA fragments in a short time. eY1H assays will be an important tool for mapping gene regulatory networks in Caenorhabditis elegans and other model organisms as well as in humans.

Yeast one-hybrid assays for gene-centered human gene regulatory network mapping.

Gateway-compatible yeast one-hybrid (Y1H) assays provide a convenient gene-centered (DNA to protein) approach to identify transcription factors that can bind a DNA sequence of interest. We present Y1H resources, including clones for 988 of 1,434 (69%) predicted human transcription factors, that can be used to detect both known and new interactions between human DNA regions and transcription factors.

Enhanced Y1H assays for Arabidopsis.

We present an Arabidopsis thaliana full-length transcription factor resource of 92% of root stele-expressed transcription factors and 74.5% of root-expressed transcription factors. We demonstrate its use with enhanced yeast one-hybrid (eY1H) screening for rapid, systematic mapping of plant transcription factor-promoter interactions. We identified 158 interactions with 13 stele-expressed promoters, many of which occur physically or are regulatory in planta.

Yeast one-hybrid assays: a historical and technical perspective.

Since its development about two decades ago, the yeast one-hybrid (Y1H) assay has become an important technique for detecting physical interactions between sequence-specific regulatory transcription factor proteins (TFs) and their DNA target sites. Multiple versions of the Y1H methodology have been developed, each with technical differences and unique advantages. We will discuss several of these technical variations in detail, and also provide some ideas for how Y1H assays can be further improved.

Discovery of Ligands for TRIM58, a Novel Tissue-Selective E3 Ligase.

Redirecting E3 ligases to neo-substrates, leading to their proteasomal disassembly, known as targeted protein degradation (TPD), has emerged as a promising alternative to traditional, occupancy-driven pharmacology. Although the field has expanded tremendously over the past years, the choice of E3 ligases remains limited, with an almost exclusive focus on CRBN and VHL. Here, we report the discovery of novel ligands to the PRY-SPRY domain of TRIM58, a RING ligase that is specifically expressed in erythroid precursor cells. A DSF screen, followed by validation using additional biophysical methods, led to the identification of TRIM58 ligand TRIM-473. A basic SAR around the chemotype was established by utilizing a competitive binding assay employing a short FP peptide probe derived from an endogenous TRIM58 substrate. The X-ray co-crystal structure of TRIM58 in complex with TRIM-473 gave insights into the binding mode and potential exit vectors for bifunctional degrader design.

Cell adhesion molecule KIRREL1 is a feedback regulator of Hippo signaling recruiting SAV1 to cell-cell contact sites.

The Hippo/YAP pathway controls cell proliferation through sensing physical and spatial organization of cells. How cell-cell contact is sensed by Hippo signaling is poorly understood. Here, we identified the cell adhesion molecule KIRREL1 as an upstream positive regulator of the mammalian Hippo pathway. KIRREL1 physically interacts with SAV1 and recruits SAV1 to cell-cell contact sites. Consistent with the hypothesis that KIRREL1-mediated cell adhesion suppresses YAP activity, knockout of KIRREL1 increases YAP activity in neighboring cells. Analyzing pan-cancer CRISPR proliferation screen data reveals KIRREL1 as the top plasma membrane protein showing strong correlation with known Hippo regulators, highlighting a critical role of KIRREL1 in regulating Hippo signaling and cell proliferation. During liver regeneration in mice, KIRREL1 is upregulated, and its genetic ablation enhances hepatic YAP activity, hepatocyte reprogramming and biliary epithelial cell proliferation. Our data suggest that KIRREL1 functions as a feedback regulator of the mammalian Hippo pathway through sensing cell-cell interaction and recruiting SAV1 to cell-cell contact sites.

DamID in C. elegans reveals longevity-associated targets of DAF-16/FoxO.

Insulin/IGF-1 signaling controls metabolism, stress resistance and aging in Caenorhabditis elegans by regulating the activity of the DAF-16/FoxO transcription factor (TF). However, the function of DAF-16 and the topology of the transcriptional network that it crowns remain unclear. Using chromatin profiling by DNA adenine methyltransferase identification (DamID), we identified 907 genes that are bound by DAF-16. These were enriched for genes showing DAF-16-dependent upregulation in long-lived daf-2 insulin/IGF-1 receptor mutants (P=1.4e(-11)). Cross-referencing DAF-16 targets with these upregulated genes (daf-2 versus daf-16; daf-2) identified 65 genes that were DAF-16 regulatory targets. These 65 were enriched for signaling genes, including known determinants of longevity, but not for genes specifying somatic maintenance functions (e.g. detoxification, repair). This suggests that DAF-16 acts within a relatively small transcriptional subnetwork activating (but not suppressing) other regulators of stress resistance and aging, rather than directly regulating terminal effectors of longevity. For most genes bound by DAF-16::DAM, transcriptional regulation by DAF-16 was not detected, perhaps reflecting transcriptionally non-functional TF 'parking sites'. This study demonstrates the efficacy of DamID for chromatin profiling in C. elegans.

The C. elegans Snail homolog CES-1 can activate gene expression in vivo and share targets with bHLH transcription factors.

Snail-type transcription factors (TFs) are found in numerous metazoan organisms and function in a plethora of cellular and developmental processes including mesoderm and neuronal development, apoptosis and cancer. So far, Snail-type TFs are exclusively known as transcriptional repressors. They repress gene expression by recruiting transcriptional co-repressors and/or by preventing DNA binding of activators from the basic helix-loop-helix (bHLH) family of TFs to CAGGTG E-box sequences. Here we report that the Caenorhabditis elegans Snail-type TF CES-1 can activate transcription in vivo. Moreover, we provide results that suggest that CES-1 can share its binding site with bHLH TFs, in different tissues, rather than only occluding bHLH DNA binding. Together, our data indicate that there are at least two types of CES-1 target genes and, therefore, that the molecular function of Snail-type TFs is more plastic than previously appreciated.

USP7 inhibits Wnt/ß-catenin signaling through promoting stabilization of Axin.

Axin is a key scaffolding protein responsible for the formation of the ß-catenin destruction complex. Stability of Axin protein is regulated by the ubiquitin-proteasome system, and modulation of cellular concentration of Axin protein has a profound effect on Wnt/ß-catenin signaling. Although E3s promoting Axin ubiquitination have been identified, the deubiquitinase responsible for Axin deubiquitination and stabilization remains unknown. Here, we identify USP7 as a potent negative regulator of Wnt/ß-catenin signaling through CRISPR screens. Genetic ablation or pharmacological inhibition of USP7 robustly increases Wnt/ß-catenin signaling in multiple cellular systems. USP7 directly interacts with Axin through its TRAF domain, and promotes deubiquitination and stabilization of Axin. Inhibition of USP7 regulates osteoblast differentiation and adipocyte differentiation through increasing Wnt/ß-catenin signaling. Our study reveals a critical mechanism that prevents excessive degradation of Axin and identifies USP7 as a target for sensitizing cells to Wnt/ß-catenin signaling.

Genome-wide CRISPR screening reveals genetic modifiers of mutant EGFR dependence in human NSCLC.

EGFR-mutant NSCLCs frequently respond to EGFR tyrosine kinase inhibitors (TKIs). However, the responses are not durable, and the magnitude of tumor regression is variable, suggesting the existence of genetic modifiers of EGFR dependency. Here, we applied a genome-wide CRISPR-Cas9 screening to identify genetic determinants of EGFR TKI sensitivity and uncovered putative candidates. We show that knockout of RIC8A, essential for G-alpha protein activation, enhanced EGFR TKI-induced cell death. Mechanistically, we demonstrate that RIC8A is a positive regulator of YAP signaling, activation of which rescued the EGFR TKI sensitizing phenotype resulting from RIC8A knockout. We also show that knockout of ARIH2, or other components in the Cullin-5 E3 complex, conferred resistance to EGFR inhibition, in part by promoting nascent protein synthesis through METAP2. Together, these data uncover a spectrum of previously unidentified regulators of EGFR TKI sensitivity in EGFR-mutant human NSCLC, providing insights into the heterogeneity of EGFR TKI treatment responses.