A model for accurate quantification of CRISPR effects in pooled FACS screens.

CRISPR screens are powerful tools to identify key genes that underlie biological processes. One important type of screen uses fluorescence activated cell sorting (FACS) to sort perturbed cells into bins based on the expression level of marker genes, followed by guide RNA (gRNA) sequencing. Analysis of these data presents several statistical challenges due to multiple factors including the discrete nature of the bins and typically small numbers of replicate experiments. To address these challenges, we developed a robust and powerful Bayesian random effects model and software package called Waterbear. Furthermore, we used Waterbear to explore how various experimental design parameters affect statistical power to establish principled guidelines for future screens. Finally, we experimentally validated our experimental design model findings that, when using Waterbear for analysis, high power is maintained even at low cell coverage and a high multiplicity of infection. We anticipate that Waterbear will be of broad utility for analyzing FACS-based CRISPR screens.

Systematic decoding of cis gene regulation defines context-dependent control of the multi-gene costimulatory receptor locus in human T cells.

Cis-regulatory elements (CREs) interact with trans regulators to orchestrate gene expression, but how transcriptional regulation is coordinated in multi-gene loci has not been experimentally defined. We sought to characterize the CREs controlling dynamic expression of the adjacent costimulatory genes CD28, CTLA4 and ICOS, encoding regulators of T cell-mediated immunity. Tiling CRISPR interference (CRISPRi) screens in primary human T cells, both conventional and regulatory subsets, uncovered gene-, cell subset- and stimulation-specific CREs. Integration with CRISPR knockout screens and assay for transposase-accessible chromatin with sequencing (ATAC-seq) profiling identified trans regulators influencing chromatin states at specific CRISPRi-responsive elements to control costimulatory gene expression. We then discovered a critical CCCTC-binding factor (CTCF) boundary that reinforces CRE interaction with CTLA4 while also preventing promiscuous activation of CD28. By systematically mapping CREs and associated trans regulators directly in primary human T cell subsets, this work overcomes longstanding experimental limitations to decode context-dependent gene regulatory programs in a complex, multi-gene locus critical to immune homeostasis.

Gene regulatory network inference from CRISPR perturbations in primary CD4+ T cells elucidates the genomic basis of immune disease.

The effects of genetic variation on complex traits act mainly through changes in gene regulation. Although many genetic variants have been linked to target genes in cis, the trans-regulatory cascade mediating their effects remains largely uncharacterized. Mapping trans-regulators based on natural genetic variation, including eQTL mapping, has been challenging due to small effects. Experimental perturbation approaches offer a complementary and powerful approach to mapping trans-regulators. We used CRISPR knockouts of 84 genes in primary CD4+ T cells to perturb an immune cell gene network, targeting both inborn error of immunity (IEI) disease transcription factors (TFs) and background TFs matched in constraint and expression level, but without a known immune disease association. We developed a novel Bayesian structure learning method called Linear Latent Causal Bayes (LLCB) to estimate the gene regulatory network from perturbation data and observed 211 directed edges among the genes which could not be detected in existing CD4+ trans-eQTL data. We used LLCB to characterize the differences between the IEI and background TFs, finding that the gene groups were highly interconnected, but that IEI TFs were much more likely to regulate immune cell specific pathways and immune GWAS genes. We further characterized nine coherent gene programs based on downstream effects of the TFs and linked these modules to regulation of GWAS genes, finding that canonical JAK-STAT family members are regulated by KMT2A, a global epigenetic regulator. These analyses reveal the trans-regulatory cascade from upstream epigenetic regulator to intermediate TFs to downstream effector cytokines and elucidate the logic linking immune GWAS genes to key signaling pathways.

CRISPR screens decode cancer cell pathways that trigger γδ T cell detection.

γδ T cells are potent anticancer effectors with the potential to target tumours broadly, independent of patient-specific neoantigens or human leukocyte antigen background1-5. γδ T cells can sense conserved cell stress signals prevalent in transformed cells2,3, although the mechanisms behind the targeting of stressed target cells remain poorly characterized. Vγ9Vδ2 T cells-the most abundant subset of human γδ T cells4-recognize a protein complex containing butyrophilin 2A1 (BTN2A1) and BTN3A1 (refs. 6-8), a widely expressed cell surface protein that is activated by phosphoantigens abundantly produced by tumour cells. Here we combined genome-wide CRISPR screens in target cancer cells to identify pathways that regulate γδ T cell killing and BTN3A cell surface expression. The screens showed previously unappreciated multilayered regulation of BTN3A abundance on the cell surface and triggering of γδ T cells through transcription, post-translational modifications and membrane trafficking. In addition, diverse genetic perturbations and inhibitors disrupting metabolic pathways in the cancer cells, particularly ATP-producing processes, were found to alter BTN3A levels. This induction of both BTN3A and BTN2A1 during metabolic crises is dependent on AMP-activated protein kinase (AMPK). Finally, small-molecule activation of AMPK in a cell line model and in patient-derived tumour organoids led to increased expression of the BTN2A1-BTN3A complex and increased Vγ9Vδ2 T cell receptor-mediated killing. This AMPK-dependent mechanism of metabolic stress-induced ligand upregulation deepens our understanding of γδ T cell stress surveillance and suggests new avenues available to enhance γδ T cell anticancer activity.

Systematic discovery and perturbation of regulatory genes in human T cells reveals the architecture of immune networks.

Gene regulatory networks ensure that important genes are expressed at precise levels. When gene expression is sufficiently perturbed, it can lead to disease. To understand how gene expression disruptions percolate through a network, we must first map connections between regulatory genes and their downstream targets. However, we lack comprehensive knowledge of the upstream regulators of most genes. Here, we developed an approach for systematic discovery of upstream regulators of critical immune factors-IL2RA, IL-2 and CTLA4-in primary human T cells. Then, we mapped the network of the target genes of these regulators and putative cis-regulatory elements using CRISPR perturbations, RNA-seq and ATAC-seq. These regulators form densely interconnected networks with extensive feedback loops. Furthermore, this network is enriched for immune-associated disease variants and genes. These results provide insight into how immune-associated disease genes are regulated in T cells and broader principles about the structure of human gene regulatory networks.

G1/S restriction point coordinates phasic gene expression and cell differentiation.

Pluripotent embryonic stem cells have a unique cell cycle structure with a suppressed G1/S restriction point and little differential expression across the cell cycle phases. Here, we evaluate the link between G1/S restriction point activation, phasic gene expression, and cellular differentiation. Expression analysis reveals a gain in phasic gene expression across lineages between embryonic days E7.5 and E9.5. Genetic manipulation of the G1/S restriction point regulators miR-302 and P27 respectively accelerates or delays the onset of phasic gene expression in mouse embryos. Loss of miR-302-mediated p21 or p27 suppression expedites embryonic stem cell differentiation, while a constitutive Cyclin E mutant blocks it. Together, these findings uncover a causal relationship between emergence of the G1/S restriction point with a gain in phasic gene expression and cellular differentiation.

Cell-cell adhesions in embryonic stem cells regulate the stability and transcriptional activity of ß-catenin.

E-cadherin (CDH1) is involved in maintaining cell-cell adhesions in embryonic stem cells (ESCs). However, its function in the context of cell fate decisions is largely unknown. Using mouse ESCs (mESCs), we demonstrate that E-cadherin and ß-catenin interact at the membrane and continue to do so upon internalization within the cell. Cdh1-/- mESCs failed to form tight colonies, with altered differentiation, marker expression and retention of pluripotency factors during differentiation. Interestingly, Cdh1-/- mESCs showed dramatically reduced ß-catenin levels. Transcriptional profiling of Cdh1-/- mESCs displayed a significant alteration in the expression of a subset of ß-catenin targets in a cell state- and GSK3ß-dependent manner. Our findings hint at hitherto unknown roles played by E-cadherin in regulating the activity of ß-catenin in ESCs.

CRISPR activation and interference screens decode stimulation responses in primary human T cells.

Regulation of cytokine production in stimulated T cells can be disrupted in autoimmunity, immunodeficiencies, and cancer. Systematic discovery of stimulation-dependent cytokine regulators requires both loss-of-function and gain-of-function studies, which have been challenging in primary human cells. We now report genome-wide CRISPR activation (CRISPRa) and interference (CRISPRi) screens in primary human T cells to identify gene networks controlling interleukin-2 (IL-2) and interferon-γ (IFN-γ) production. Arrayed CRISPRa confirmed key hits and enabled multiplexed secretome characterization, revealing reshaped cytokine responses. Coupling CRISPRa screening with single-cell RNA sequencing enabled deep molecular characterization of screen hits, revealing how perturbations tuned T cell activation and promoted cell states characterized by distinct cytokine expression profiles. These screens reveal genes that reprogram critical immune cell functions, which could inform the design of immunotherapies.

High-throughput, Efficient, and Unbiased Capture of Small RNAs from Low-input Samples for Sequencing.

MicroRNAs hold great promise as biomarkers of disease. However, there are few efficient and robust methods for measuring microRNAs from low input samples. Here, we develop a high-throughput sequencing protocol that efficiently captures small RNAs while minimizing inherent biases associated with library production. The protocol is based on early barcoding such that all downstream manipulations can be performed on a pool of many samples thereby reducing reagent usage and workload. We show that the optimization of adapter concentrations along with the addition of nucleotide modifications and random nucleotides increases the efficiency of small RNA capture. We further show, using unique molecular identifiers, that stochastic capture of low input RNA rather than PCR amplification influences the biased quantitation of intermediately and lowly expressed microRNAs. Our improved method allows the processing of tens to hundreds of samples simultaneously while retaining high efficiency quantitation of microRNAs in low input samples from tissues or bodily fluids.

GRHL2-Dependent Enhancer Switching Maintains a Pluripotent Stem Cell Transcriptional Subnetwork after Exit from Naive Pluripotency.

The enhancer landscape of pluripotent stem cells undergoes extensive reorganization during early mammalian development. The functions and mechanisms behind such reorganization, however, are unclear. Here, we show that the transcription factor GRHL2 is necessary and sufficient to activate an epithelial subset of enhancers as naive embryonic stem cells (ESCs) transition into formative epiblast-like cells (EpiLCs). Surprisingly, many GRHL2 target genes do not change in expression during the ESC-EpiLC transition. Instead, enhancers regulating these genes in ESCs diminish in activity in EpiLCs while GRHL2-dependent alternative enhancers become activated to maintain transcription. GRHL2 therefore assumes control over a subset of the naive network via enhancer switching to maintain expression of epithelial genes upon exit from naive pluripotency. These data evoke a model where the naive pluripotency network becomes partitioned into smaller, independent networks regulated by EpiLC-specific transcription factors, thereby priming cells for lineage specification.

Decoupling the impact of microRNAs on translational repression versus RNA degradation in embryonic stem cells.

Translation and mRNA degradation are intimately connected, yet the mechanisms that link them are not fully understood. Here, we studied these mechanisms in embryonic stem cells (ESCs). Transcripts showed a wide range of stabilities, which correlated with their relative translation levels and that did not change during early ESC differentiation. The protein DHH1 links translation to mRNA stability in yeast; however, loss of the mammalian homolog, DDX6, in ESCs did not disrupt the correlation across transcripts. Instead, the loss of DDX6 led to upregulated translation of microRNA targets, without concurrent changes in mRNA stability. The Ddx6 knockout cells were phenotypically and molecularly similar to cells lacking all microRNAs (Dgcr8 knockout ESCs). These data show that the loss of DDX6 can separate the two canonical functions of microRNAs: translational repression and transcript destabilization. Furthermore, these data uncover a central role for translational repression independent of transcript destabilization in defining the downstream consequences of microRNA loss.

Expression of Alternative Ago2 Isoform Associated with Loss of microRNA-Driven Translational Repression in Mouse Oocytes.

Mouse oocyte maturation, fertilization, and reprogramming occur in the absence of transcription, and thus, changes in mRNA levels and translation rate are regulated through post-transcriptional mechanisms [1]. Surprisingly, microRNA function, which is a major form of post-transcriptional regulation, is absent during this critical period of mammalian development [2, 3]. Here, we investigated the mechanisms underlying the global suppression of microRNA activity. In both mouse and frogs, microRNA function was active in growing oocytes but then absent during oocyte maturation. RNA sequencing (RNA-seq) of mouse oocytes uncovered that the microRNA effector protein AGO2 is predominantly expressed as an alternative isoform that encodes a truncated protein lacking all of the known essential domains. Full-length Ago2 as well as the related Argonautes (Ago1, Ago3, and Ago4) were lowly expressed in maturing mouse oocytes. Reintroduction of full-length AGO2 together with an exogenous microRNA in either mouse or frog oocytes restored translational repression of a target reporter. However, levels of endogenous transcripts remained unchanged. Consistent with a lack of microRNA activity, analysis of transcripts with alternative polyadenylation sites showed increased stability of transcripts with a longer 3' UTR during oocyte maturation. Redundant mechanisms protecting endogenous transcripts and the conserved loss of microRNA activity suggest a strong selection for suppressing microRNA function in vertebrate oocytes.