Genetic Screen for Cell Fitness in High or Low Oxygen Highlights Mitochondrial and Lipid Metabolism.

Human cells are able to sense and adapt to variations in oxygen levels. Historically, much research in this field has focused on hypoxia-inducible factor (HIF) signaling and reactive oxygen species (ROS). Here, we perform genome-wide CRISPR growth screens at 21%, 5%, and 1% oxygen to systematically identify gene knockouts with relative fitness defects in high oxygen (213 genes) or low oxygen (109 genes), most without known connection to HIF or ROS. Knockouts of many mitochondrial pathways thought to be essential, including complex I and enzymes in Fe-S biosynthesis, grow relatively well at low oxygen and thus are buffered by hypoxia. In contrast, in certain cell types, knockout of lipid biosynthetic and peroxisomal genes causes fitness defects only in low oxygen. Our resource nominates genetic diseases whose severity may be modulated by oxygen and links hundreds of genes to oxygen homeostasis.

High-Throughput Discovery and Characterization of Human Transcriptional Effectors.

Thousands of proteins localize to the nucleus; however, it remains unclear which contain transcriptional effectors. Here, we develop HT-recruit, a pooled assay where protein libraries are recruited to a reporter, and their transcriptional effects are measured by sequencing. Using this approach, we measure gene silencing and activation for thousands of domains. We find a relationship between repressor function and evolutionary age for the KRAB domains, discover that Homeodomain repressor strength is collinear with Hox genetic organization, and identify activities for several domains of unknown function. Deep mutational scanning of the CRISPRi KRAB maps the co-repressor binding surface and identifies substitutions that improve stability/silencing. By tiling 238 proteins, we find repressors as short as ten amino acids. Finally, we report new activator domains, including a divergent KRAB. These results provide a resource of 600 human proteins containing effectors and demonstrate a scalable strategy for assigning functions to protein domains.

Structural Basis of Teneurin-Latrophilin Interaction in Repulsive Guidance of Migrating Neurons.

Teneurins are ancient metazoan cell adhesion receptors that control brain development and neuronal wiring in higher animals. The extracellular C terminus binds the adhesion GPCR Latrophilin, forming a trans-cellular complex with synaptogenic functions. However, Teneurins, Latrophilins, and FLRT proteins are also expressed during murine cortical cell migration at earlier developmental stages. Here, we present crystal structures of Teneurin-Latrophilin complexes that reveal how the lectin and olfactomedin domains of Latrophilin bind across a spiraling beta-barrel domain of Teneurin, the YD shell. We couple structure-based protein engineering to biophysical analysis, cell migration assays, and in utero electroporation experiments to probe the importance of the interaction in cortical neuron migration. We show that binding of Latrophilins to Teneurins and FLRTs directs the migration of neurons using a contact repulsion-dependent mechanism. The effect is observed with cell bodies and small neurites rather than their processes. The results exemplify how a structure-encoded synaptogenic protein complex is also used for repulsive cell guidance.

Identification of RIOK2 as a master regulator of human blood cell development.

Anemia is a major comorbidity in aging, chronic kidney and inflammatory diseases, and hematologic malignancies. However, the transcriptomic networks governing hematopoietic differentiation in blood cell development remain incompletely defined. Here we report that the atypical kinase RIOK2 (right open reading frame kinase 2) is a master transcription factor (TF) that not only drives erythroid differentiation, but also simultaneously suppresses megakaryopoiesis and myelopoiesis in primary human stem and progenitor cells. Our study reveals the previously uncharacterized winged helix-turn-helix DNA-binding domain and two transactivation domains of RIOK2 that are critical to regulate key hematopoietic TFs GATA1, GATA2, SPI1, RUNX3 and KLF1. This establishes RIOK2 as an integral component of the transcriptional regulatory network governing human hematopoietic differentiation. Importantly, RIOK2 mRNA expression significantly correlates with these TFs and other hematopoietic genes in myelodysplastic syndromes, acute myeloid leukemia and chronic kidney disease. Further investigation of RIOK2-mediated transcriptional pathways should yield therapeutic approaches to correct defective hematopoiesis in hematologic disorders.

Disentangling Genetic and Environmental Effects on the Proteotypes of Individuals.

Proteotypes, like genotypes, have been found to vary between individuals in several studies, but consistent molecular functional traits across studies remain to be quantified. In a meta-analysis of 11 proteomics datasets from humans and mice, we use co-variation of proteins in known functional modules across datasets and individuals to obtain a consensus landscape of proteotype variation. We find that individuals differ considerably in both protein complex abundances and stoichiometry. We disentangle genetic and environmental factors impacting these metrics, with genetic sex and specific diets together explaining 13.5% and 11.6% of the observed variation of complex abundance and stoichiometry, respectively. Sex-specific differences, for example, include various proteins and complexes, where the respective genes are not located on sex-specific chromosomes. Diet-specific differences, added to the individual genetic backgrounds, might become a starting point for personalized proteotype modulation toward desired features.

Promoter-Intrinsic and Local Chromatin Features Determine Gene Repression in LADs.

It is largely unclear whether genes that are naturally embedded in lamina-associated domains (LADs) are inactive due to their chromatin environment or whether LADs are merely secondary to the lack of transcription. We show that hundreds of human promoters become active when moved from their native LAD position to a neutral context in the same cells, indicating that LADs form a repressive environment. Another set of promoters inside LADs is able to "escape" repression, although their transcription elongation is attenuated. By inserting reporters into thousands of genomic locations, we demonstrate that escaper promoters are intrinsically less sensitive to LAD repression. This is not simply explained by promoter strength but by the interplay between promoter sequence and local chromatin features that vary strongly across LADs. Enhancers also differ in their sensitivity to LAD chromatin. This work provides a general framework for the systematic understanding of gene regulation by repressive chromatin.

Pervasive Chromatin-RNA Binding Protein Interactions Enable RNA-Based Regulation of Transcription.

Increasing evidence suggests that transcriptional control and chromatin activities at large involve regulatory RNAs, which likely enlist specific RNA-binding proteins (RBPs). Although multiple RBPs have been implicated in transcription control, it has remained unclear how extensively RBPs directly act on chromatin. We embarked on a large-scale RBP ChIP-seq analysis, revealing widespread RBP presence in active chromatin regions in the human genome. Like transcription factors (TFs), RBPs also show strong preference for hotspots in the genome, particularly gene promoters, where their association is frequently linked to transcriptional output. Unsupervised clustering reveals extensive co-association between TFs and RBPs, as exemplified by YY1, a known RNA-dependent TF, and RBM25, an RBP involved in splicing regulation. Remarkably, RBM25 depletion attenuates all YY1-dependent activities, including chromatin binding, DNA looping, and transcription. We propose that various RBPs may enhance network interaction through harnessing regulatory RNAs to control transcription.

Hypoxia Rescues Frataxin Loss by Restoring Iron Sulfur Cluster Biogenesis.

Friedreich's ataxia (FRDA) is a devastating, multisystemic disorder caused by recessive mutations in the mitochondrial protein frataxin (FXN). FXN participates in the biosynthesis of Fe-S clusters and is considered to be essential for viability. Here we report that when grown in 1% ambient O2, FXN null yeast, human cells, and nematodes are fully viable. In human cells, hypoxia restores steady-state levels of Fe-S clusters and normalizes ATF4, NRF2, and IRP2 signaling events associated with FRDA. Cellular studies and in vitro reconstitution indicate that hypoxia acts through HIF-independent mechanisms that increase bioavailable iron as well as directly activate Fe-S synthesis. In a mouse model of FRDA, breathing 11% O2 attenuates the progression of ataxia, whereas breathing 55% O2 hastens it. Our work identifies oxygen as a key environmental variable in the pathogenesis associated with FXN depletion, with important mechanistic and therapeutic implications.

Visualizing Engrafted Human Cancer and Therapy Responses in Immunodeficient Zebrafish.

Xenograft cell transplantation into immunodeficient mice has become the gold standard for assessing pre-clinical efficacy of cancer drugs, yet direct visualization of single-cell phenotypes is difficult. Here, we report an optically-clear prkdc-/-, il2rga-/- zebrafish that lacks adaptive and natural killer immune cells, can engraft a wide array of human cancers at 37°C, and permits the dynamic visualization of single engrafted cells. For example, photoconversion cell-lineage tracing identified migratory and proliferative cell states in human rhabdomyosarcoma, a pediatric cancer of muscle. Additional experiments identified the preclinical efficacy of combination olaparib PARP inhibitor and temozolomide DNA-damaging agent as an effective therapy for rhabdomyosarcoma and visualized therapeutic responses using a four-color FUCCI cell-cycle fluorescent reporter. These experiments identified that combination treatment arrested rhabdomyosarcoma cells in the G2 cell cycle prior to induction of apoptosis. Finally, patient-derived xenografts could be engrafted into our model, opening new avenues for developing personalized therapeutic approaches in the future.

A Compendium of Genetic Modifiers of Mitochondrial Dysfunction Reveals Intra-organelle Buffering.

Mitochondrial dysfunction is associated with a spectrum of human conditions, ranging from rare, inborn errors of metabolism to the aging process. To identify pathways that modify mitochondrial dysfunction, we performed genome-wide CRISPR screens in the presence of small-molecule mitochondrial inhibitors. We report a compendium of chemical-genetic interactions involving 191 distinct genetic modifiers, including 38 that are synthetic sick/lethal and 63 that are suppressors. Genes involved in glycolysis (PFKP), pentose phosphate pathway (G6PD), and defense against lipid peroxidation (GPX4) scored high as synthetic sick/lethal. A surprisingly large fraction of suppressors are pathway intrinsic and encode mitochondrial proteins. A striking example of such "intra-organelle" buffering is the alleviation of a chemical defect in complex V by simultaneous inhibition of complex I, which benefits cells by rebalancing redox cofactors, increasing reductive carboxylation, and promoting glycolysis. Perhaps paradoxically, certain forms of mitochondrial dysfunction may best be buffered with "second site" inhibitors to the organelle.

Modulation of Protein-Interaction States through the Cell Cycle.

Global profiling of protein expression through the cell cycle has revealed subsets of periodically expressed proteins. However, expression levels alone only give a partial view of the biochemical processes determining cellular events. Using a proteome-wide implementation of the cellular thermal shift assay (CETSA) to study specific cell-cycle phases, we uncover changes of interaction states for more than 750 proteins during the cell cycle. Notably, many protein complexes are modulated in specific cell-cycle phases, reflecting their roles in processes such as DNA replication, chromatin remodeling, transcription, translation, and disintegration of the nuclear envelope. Surprisingly, only small differences in the interaction states were seen between the G1 and the G2 phase, suggesting similar hardwiring of biochemical processes in these two phases. The present work reveals novel molecular details of the cell cycle and establishes proteome-wide CETSA as a new strategy to study modulation of protein-interaction states in intact cells.

Universal Chimeric Antigen Receptors for Multiplexed and Logical Control of T Cell Responses.

T cells expressing chimeric antigen receptors (CARs) are promising cancer therapeutic agents, with the prospect of becoming the ultimate smart cancer therapeutics. To expand the capability of CAR T cells, here, we present a split, universal, and programmable (SUPRA) CAR system that simultaneously encompasses multiple critical "upgrades," such as the ability to switch targets without re-engineering the T cells, finely tune T cell activation strength, and sense and logically respond to multiple antigens. These features are useful to combat relapse, mitigate over-activation, and enhance specificity. We test our SUPRA system against two different tumor models to demonstrate its broad utility and humanize its components to minimize potential immunogenicity concerns. Furthermore, we extend the orthogonal SUPRA CAR system to regulate different T cell subsets independently, demonstrating a dually inducible CAR system. Together, these SUPRA CARs illustrate that multiple advanced logic and control features can be implemented into a single, integrated system.

Mapping the Genetic Landscape of Human Cells.

Seminal yeast studies have established the value of comprehensively mapping genetic interactions (GIs) for inferring gene function. Efforts in human cells using focused gene sets underscore the utility of this approach, but the feasibility of generating large-scale, diverse human GI maps remains unresolved. We developed a CRISPR interference platform for large-scale quantitative mapping of human GIs. We systematically perturbed 222,784 gene pairs in two cancer cell lines. The resultant maps cluster functionally related genes, assigning function to poorly characterized genes, including TMEM261, a new electron transport chain component. Individual GIs pinpoint unexpected relationships between pathways, exemplified by a specific cholesterol biosynthesis intermediate whose accumulation induces deoxynucleotide depletion, causing replicative DNA damage and a synthetic-lethal interaction with the ATR/9-1-1 DNA repair pathway. Our map provides a broad resource, establishes GI maps as a high-resolution tool for dissecting gene function, and serves as a blueprint for mapping the genetic landscape of human cells.

An Integrated Genome-wide CRISPRa Approach to Functionalize lncRNAs in Drug Resistance.

Resistance to chemotherapy plays a significant role in cancer mortality. To identify genetic units affecting sensitivity to cytarabine, the mainstay of treatment for acute myeloid leukemia (AML), we developed a comprehensive and integrated genome-wide platform based on a dual protein-coding and non-coding integrated CRISPRa screening (DICaS). Putative resistance genes were initially identified using pharmacogenetic data from 760 human pan-cancer cell lines. Subsequently, genome scale functional characterization of both coding and long non-coding RNA (lncRNA) genes by CRISPR activation was performed. For lncRNA functional assessment, we developed a CRISPR activation of lncRNA (CaLR) strategy, targeting 14,701 lncRNA genes. Computational and functional analysis identified novel cell-cycle, survival/apoptosis, and cancer signaling genes. Furthermore, transcriptional activation of the GAS6-AS2 lncRNA, identified in our analysis, leads to hyperactivation of the GAS6/TAM pathway, a resistance mechanism in multiple cancers including AML. Thus, DICaS represents a novel and powerful approach to identify integrated coding and non-coding pathways of therapeutic relevance.

Efficient Generation of Transcriptomic Profiles by Random Composite Measurements.

RNA profiles are an informative phenotype of cellular and tissue states but can be costly to generate at massive scale. Here, we describe how gene expression levels can be efficiently acquired with random composite measurements-in which abundances are combined in a random weighted sum. We show (1) that the similarity between pairs of expression profiles can be approximated with very few composite measurements; (2) that by leveraging sparse, modular representations of gene expression, we can use random composite measurements to recover high-dimensional gene expression levels (with 100 times fewer measurements than genes); and (3) that it is possible to blindly recover gene expression from composite measurements, even without access to training data. Our results suggest new compressive modalities as a foundation for massive scaling in high-throughput measurements and new insights into the interpretation of high-dimensional data.

In Situ Capture of Chromatin Interactions by Biotinylated dCas9.

Cis-regulatory elements (CREs) are commonly recognized by correlative chromatin features, yet the molecular composition of the vast majority of CREs in chromatin remains unknown. Here, we describe a CRISPR affinity purification in situ of regulatory elements (CAPTURE) approach to unbiasedly identify locus-specific chromatin-regulating protein complexes and long-range DNA interactions. Using an in vivo biotinylated nuclease-deficient Cas9 protein and sequence-specific guide RNAs, we show high-resolution and selective isolation of chromatin interactions at a single-copy genomic locus. Purification of human telomeres using CAPTURE identifies known and new telomeric factors. In situ capture of individual constituents of the enhancer cluster controlling human ß-globin genes establishes evidence for composition-based hierarchical organization. Furthermore, unbiased analysis of chromatin interactions at disease-associated cis-elements and developmentally regulated super-enhancers reveals spatial features that causally control gene transcription. Thus, comprehensive and unbiased analysis of locus-specific regulatory composition provides mechanistic insight into genome structure and function in development and disease.

DEL-1 promotes macrophage efferocytosis and clearance of inflammation.

Resolution of inflammation is essential for tissue homeostasis and represents a promising approach to inflammatory disorders. Here we found that developmental endothelial locus-1 (DEL-1), a secreted protein that inhibits leukocyte-endothelial adhesion and inflammation initiation, also functions as a non-redundant downstream effector in inflammation clearance. In human and mouse periodontitis, waning of inflammation was correlated with DEL-1 upregulation, whereas resolution of experimental periodontitis failed in DEL-1 deficiency. This concept was mechanistically substantiated in acute monosodium-urate-crystal-induced inflammation, where the pro-resolution function of DEL-1 was attributed to effective apoptotic neutrophil clearance (efferocytosis). DEL-1-mediated efferocytosis induced liver X receptor-dependent macrophage reprogramming to a pro-resolving phenotype and was required for optimal production of at least certain specific pro-resolving mediators. Experiments in transgenic mice with cell-specific overexpression of DEL-1 linked its anti-leukocyte-recruitment action to endothelial cell-derived DEL-1 and its efferocytic/pro-resolving action to macrophage-derived DEL-1. Thus, the compartmentalized expression of DEL-1 facilitates distinct homeostatic functions in an appropriate context that can be harnessed therapeutically.

Perturb-Seq: Dissecting Molecular Circuits with Scalable Single-Cell RNA Profiling of Pooled Genetic Screens.

Genetic screens help infer gene function in mammalian cells, but it has remained difficult to assay complex phenotypes-such as transcriptional profiles-at scale. Here, we develop Perturb-seq, combining single-cell RNA sequencing (RNA-seq) and clustered regularly interspaced short palindromic repeats (CRISPR)-based perturbations to perform many such assays in a pool. We demonstrate Perturb-seq by analyzing 200,000 cells in immune cells and cell lines, focusing on transcription factors regulating the response of dendritic cells to lipopolysaccharide (LPS). Perturb-seq accurately identifies individual gene targets, gene signatures, and cell states affected by individual perturbations and their genetic interactions. We posit new functions for regulators of differentiation, the anti-viral response, and mitochondrial function during immune activation. By decomposing many high content measurements into the effects of perturbations, their interactions, and diverse cell metadata, Perturb-seq dramatically increases the scope of pooled genomic assays.

CRISPR Inversion of CTCF Sites Alters Genome Topology and Enhancer/Promoter Function.

CTCF and the associated cohesin complex play a central role in insulator function and higher-order chromatin organization of mammalian genomes. Recent studies identified a correlation between the orientation of CTCF-binding sites (CBSs) and chromatin loops. To test the functional significance of this observation, we combined CRISPR/Cas9-based genomic-DNA-fragment editing with chromosome-conformation-capture experiments to show that the location and relative orientations of CBSs determine the specificity of long-range chromatin looping in mammalian genomes, using protocadherin (Pcdh) and ß-globin as model genes. Inversion of CBS elements within the Pcdh enhancer reconfigures the topology of chromatin loops between the distal enhancer and target promoters and alters gene-expression patterns. Thus, although enhancers can function in an orientation-independent manner in reporter assays, in the native chromosome context, the orientation of at least some enhancers carrying CBSs can determine both the architecture of topological chromatin domains and enhancer/promoter specificity. These findings reveal how 3D chromosome architecture can be encoded by linear genome sequences.

Improving prime editing with an endogenous small RNA-binding protein.

Prime editing enables the precise modification of genomes through reverse transcription of template sequences appended to the 3' ends of CRISPR-Cas guide RNAs1. To identify cellular determinants of prime editing, we developed scalable prime editing reporters and performed genome-scale CRISPR-interference screens. From these screens, a single factor emerged as the strongest mediator of prime editing: the small RNA-binding exonuclease protection factor La. Further investigation revealed that La promotes prime editing across approaches (PE2, PE3, PE4 and PE5), edit types (substitutions, insertions and deletions), endogenous loci and cell types but has no consistent effect on genome-editing approaches that rely on standard, unextended guide RNAs. Previous work has shown that La binds polyuridine tracts at the 3' ends of RNA polymerase III transcripts2. We found that La functionally interacts with the 3' ends of polyuridylated prime editing guide RNAs (pegRNAs). Guided by these results, we developed a prime editor protein (PE7) fused to the RNA-binding, N-terminal domain of La. This editor improved prime editing with expressed pegRNAs and engineered pegRNAs (epegRNAs), as well as with synthetic pegRNAs optimized for La binding. Together, our results provide key insights into how prime editing components interact with the cellular environment and suggest general strategies for stabilizing exogenous small RNAs therein.