Deciphering regulation of FOXP3 expression in human conventional T cells.

FOXP3 is a lineage-defining transcription factor that controls differentiation and maintenance of suppressive function of regulatory T cells (Tregs). Foxp3 is exclusively expressed in Tregs in mice. However, in humans, FOXP3 is not only constitutively expressed in Tregs; it is also transiently expressed in stimulated CD4+CD25- conventional T cells (Tconvs)1-3. Mechanisms governing the expression of FOXP3 in human Tconvs are not understood. Here, we performed CRISPR interference (CRISPRi) screens using a 15K-member gRNA library tiling 39 kb downstream of the FOXP3 transcriptional start site (TSS) to 85 kb upstream of the TSS in Treg and Tconvs. The FOXP3 promoter and conserved non-coding sequences (CNS0, CNS1, CNS2 and CNS3), characterized as enhancer elements in murine Tregs, were required for maintenance of FOXP3 in human Tregs. In contrast, FOXP3 in human Tconvs depended on regulation at CNS0 and a novel Tconv-specific noncoding sequence (TcNS+) located upstream of CNS0. Arrayed validations of these sites identified an additional repressive cis-element overlapping with the PPP1R3F promoter (TcNS-). Pooled CRISPR knockouts revealed multiple transcription factors required for proper expression of FOXP3 in Tconvs, including GATA3, STAT5, IRF4, ETS1 and DNA methylation-associated regulators DNMT1 and MBD2. Analysis of ChIP-seq and ATAC-seq paired with knock-out (KO) of GATA3, STAT5, IRF4, and ETS1 revealed regulation of CNS0 and TcNS+ accessibility. Collectively, this work identified Treg-shared and Tconv-specific cis-elements and the trans-factors that interact with them, building a network of regulators controlling FOXP3 expression in human Tconvs.

Phenotypic profiling of human induced regulatory T cells at early differentiation: insights into distinct immunosuppressive potential.

Regulatory T cells (Tregs) play a key role in suppressing systemic effector immune responses, thereby preventing autoimmune diseases but also potentially contributing to tumor progression. Thus, there is great interest in clinically manipulating Tregs, but the precise mechanisms governing in vitro-induced Treg (iTreg) differentiation are not yet fully understood. Here, we used multiparametric mass cytometry to phenotypically profile human iTregs during the early stages of in vitro differentiation at single-cell level. A panel of 25 metal-conjugated antibodies specific to markers associated with human Tregs was used to characterize these immunomodulatory cells. We found that iTregs highly express the transcription factor FOXP3, as well as characteristic Treg-associated surface markers (e.g. CD25, PD1, CD137, CCR4, CCR7, CXCR3, and CD103). Expression of co-inhibitory factors (e.g. TIM3, LAG3, and TIGIT) increased slightly at late stages of iTreg differentiation. Further, CD103 was upregulated on a subpopulation of iTregs with greater suppressive capacity than their CD103- counterparts. Using mass-spectrometry-based proteomics, we showed that sorted CD103+ iTregs express factors associated with immunosuppression. Overall, our study highlights that during early stages of differentiation, iTregs resemble memory-like Treg features with immunosuppressive activity, and provides opportunities for further investigation into the molecular mechanisms underlying Treg function.

NKp30 and NKG2D contribute to natural killer recognition of HIV-infected cells.

Natural killer (NK) cells respond rapidly in early HIV-1 infection. HIV-1 prevention and control strategies harnessing NK cells could be enabled by mechanistic understanding of how NK cells recognize HIV-infected T cells. Here, we profiled the phenotype of human primary NK cells responsive to autologous HIV-1-infected CD4 + T cells in vitro. We characterized the patterns of NK cell ligand expression on CD4 + T cells at baseline and after infection with a panel of transmitted/founder HIV-1 strains to identify key receptor-ligand pairings. CRISPR editing of CD4 + T cells to knockout the NKp30 ligand B7-H6, or the NKG2D ligands MICB or ULBP2 reduced NK cell responses to HIV-infected cells in some donors. In contrast, overexpression of NKp30 or NKG2D in NK cells enhanced their targeting of HIV-infected cells. Collectively, we identified receptor-ligand pairs including NKp30:B7-H6 and NKG2D:MICB/ULBP2 that contribute to NK cell recognition of HIV-infected cells.

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.

Bidirectional epigenetic editing reveals hierarchies in gene regulation.

CRISPR perturbation methods are limited in their ability to study non-coding elements and genetic interactions. In this study, we developed a system for bidirectional epigenetic editing, called CRISPRai, in which we apply activating (CRISPRa) and repressive (CRISPRi) perturbations to two loci simultaneously in the same cell. We developed CRISPRai Perturb-seq by coupling dual perturbation gRNA detection with single-cell RNA sequencing, enabling study of pooled perturbations in a mixed single-cell population. We applied this platform to study the genetic interaction between two hematopoietic lineage transcription factors, SPI1 and GATA1, and discovered novel characteristics of their co-regulation on downstream target genes, including differences in SPI1 and GATA1 occupancy at genes that are regulated through different modes. We also studied the regulatory landscape of IL2 (interleukin-2) in Jurkat T cells, primary T cells and chimeric antigen receptor (CAR) T cells and elucidated mechanisms of enhancer-mediated IL2 gene regulation. CRISPRai facilitates investigation of context-specific genetic interactions, provides new insights into gene regulation and will enable exploration of non-coding disease-associated variants.

A proximal enhancer regulates RORA expression during early human Th17 cell differentiation.

Gene regulatory elements, such as enhancers, greatly influence cell identity by tuning the transcriptional activity of specific cell types. Dynamics of enhancer landscape during early human Th17 cell differentiation remains incompletely understood. Leveraging ATAC-seq-based profiling of chromatin accessibility and comprehensive analysis of key histone marks, we identified a repertoire of enhancers that potentially exert control over the fate specification of Th17 cells. We found 23 SNPs associated with autoimmune diseases within Th17-enhancers that precisely overlapped with the binding sites of transcription factors actively engaged in T-cell functions. Among the Th17-specific enhancers, we identified an enhancer in the intron of RORA and demonstrated that this enhancer positively regulates RORA transcription. Moreover, CRISPR-Cas9-mediated deletion of a transcription factor binding site-rich region within the identified RORA enhancer confirmed its role in regulating RORA transcription. These findings provide insights into the potential mechanism by which the RORA enhancer orchestrates Th17 differentiation.

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.

Long noncoding RNA LIRIL2R modulates FOXP3 levels and suppressive function of human CD4+ regulatory T cells by regulating IL2RA.

Regulatory T cells (Tregs) are central in controlling immune responses, and dysregulation of their function can lead to autoimmune disorders or cancer. Despite extensive studies on Tregs, the basis of epigenetic regulation of human Treg development and function is incompletely understood. Long intergenic noncoding RNAs (lincRNA)s are important for shaping and maintaining the epigenetic landscape in different cell types. In this study, we identified a gene on the chromosome 6p25.3 locus, encoding a lincRNA, that was up-regulated during early differentiation of human Tregs. The lincRNA regulated the expression of interleukin-2 receptor alpha (IL2RA), and we named it the lincRNA regulator of IL2RA (LIRIL2R). Through transcriptomics, epigenomics, and proteomics analysis of LIRIL2R-deficient Tregs, coupled with global profiling of LIRIL2R binding sites using chromatin isolation by RNA purification, followed by sequencing, we identified IL2RA as a target of LIRIL2R. This nuclear lincRNA binds upstream of the IL2RA locus and regulates its epigenetic landscape and transcription. CRISPR-mediated deletion of the LIRIL2R-bound region at the IL2RA locus resulted in reduced IL2RA expression. Notably, LIRIL2R deficiency led to reduced expression of Treg-signature genes (e.g., FOXP3, CTLA4, and PDCD1), upregulation of genes associated with effector T cells (e.g., SATB1 and GATA3), and loss of Treg-mediated suppression.

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.

Ultra-high efficiency T cell reprogramming at multiple loci with SEED-Selection.

Multiplexed reprogramming of T cell specificity and function can generate powerful next-generation cellular therapies. However, current manufacturing methods produce heterogenous mixtures of partially engineered cells. Here, we develop a one-step process to enrich for unlabeled cells with knock-ins at multiple target loci using a family of repair templates named Synthetic Exon/Expression Disruptors (SEEDs). SEED engineering associates transgene integration with the disruption of a paired endogenous surface protein, allowing non-modified and partially edited cells to be immunomagnetically depleted (SEED-Selection). We design SEEDs to fully reprogram three critical loci encoding T cell specificity, co-receptor expression, and MHC expression, with up to 98% purity after selection for individual modifications and up to 90% purity for six simultaneous edits (three knock-ins and three knockouts). These methods are simple, compatible with existing clinical manufacturing workflows, and can be readily adapted to other loci to facilitate production of complex gene-edited cell therapies.

CRISPR-Cas9 screen of E3 ubiquitin ligases identifies TRAF2 and UHRF1 as regulators of HIV latency in primary human T cells.

During HIV infection of CD4+ T cells, ubiquitin pathways are essential to viral replication and host innate immune response; however, the role of specific E3 ubiquitin ligases is not well understood. Proteomics analyses identified 116 single-subunit E3 ubiquitin ligases expressed in activated primary human CD4+ T cells. Using a CRISPR-based arrayed spreading infectivity assay, we systematically knocked out 116 E3s from activated primary CD4+ T cells and infected them with NL4-3 GFP reporter HIV-1. We found 10 E3s significantly positively or negatively affected HIV infection in activated primary CD4+ T cells, including UHRF1 (pro-viral) and TRAF2 (anti-viral). Furthermore, deletion of either TRAF2 or UHRF1 in three JLat models of latency spontaneously increased HIV transcription. To verify this effect, we developed a CRISPR-compatible resting primary human CD4+ T cell model of latency. Using this system, we found that deletion of TRAF2 or UHRF1 initiated latency reactivation and increased virus production from primary human resting CD4+ T cells, suggesting these two E3s represent promising targets for future HIV latency reversal strategies. IMPORTANCE: HIV, the virus that causes AIDS, heavily relies on the machinery of human cells to infect and replicate. Our study focuses on the host cell's ubiquitination system which is crucial for numerous cellular processes. Many pathogens, including HIV, exploit this system to enhance their own replication and survival. E3 proteins are part of the ubiquitination pathway that are useful drug targets for host-directed therapies. We interrogated the 116 E3s found in human immune cells known as CD4+ T cells, since these are the target cells infected by HIV. Using CRISPR, a gene-editing tool, we individually removed each of these enzymes and observed the impact on HIV infection in human CD4+ T cells isolated from healthy donors. We discovered that 10 of the E3 enzymes had a significant effect on HIV infection. Two of them, TRAF2 and UHRF1, modulated HIV activity within the cells and triggered an increased release of HIV from previously dormant or "latent" cells in a new primary T cell assay. This finding could guide strategies to perturb hidden HIV reservoirs, a major hurdle to curing HIV. Our study offers insights into HIV-host interactions, identifies new factors that influence HIV infection in immune cells, and introduces a novel methodology for studying HIV infection and latency in human immune cells.

Multiplexed CRISPR gene editing in primary human islet cells with Cas9 ribonucleoprotein.

Successful genome editing in primary human islets could reveal features of the genetic regulatory landscape underlying ß cell function and diabetes risk. Here, we describe a CRISPR-based strategy to interrogate functions of predicted regulatory DNA elements using electroporation of a complex of Cas9 ribonucleoprotein (Cas9 RNP) and guide RNAs into primary human islet cells. We successfully targeted coding regions including the PDX1 exon 1, and non-coding DNA linked to diabetes susceptibility. CRISPR-Cas9 RNP approaches revealed genetic targets of regulation by DNA elements containing candidate diabetes risk SNPs, including an in vivo enhancer of the MPHOSPH9 gene. CRISPR-Cas9 RNP multiplexed targeting of two cis-regulatory elements linked to diabetes risk in PCSK1, which encodes an endoprotease crucial for Insulin processing, also demonstrated efficient simultaneous editing of PCSK1 regulatory elements, resulting in impaired ß cell PCSK1 regulation and Insulin secretion. Multiplex CRISPR-Cas9 RNP provides powerful approaches to investigate and elucidate human islet cell gene regulation in health and diabetes.

Base-editing mutagenesis maps alleles to tune human T cell functions.

CRISPR-enabled screening is a powerful tool for the discovery of genes that control T cell function and has nominated candidate targets for immunotherapies1-6. However, new approaches are required to probe specific nucleotide sequences within key genes. Systematic mutagenesis in primary human T cells could reveal alleles that tune specific phenotypes. DNA base editors are powerful tools for introducing targeted mutations with high efficiency7,8. Here we develop a large-scale base-editing mutagenesis platform with the goal of pinpointing nucleotides that encode amino acid residues that tune primary human T cell activation responses. We generated a library of around 117,000 single guide RNA molecules targeting base editors to protein-coding sites across 385 genes implicated in T cell function and systematically identified protein domains and specific amino acid residues that regulate T cell activation and cytokine production. We found a broad spectrum of alleles with variants encoding critical residues in proteins including PIK3CD, VAV1, LCP2, PLCG1 and DGKZ, including both gain-of-function and loss-of-function mutations. We validated the functional effects of many alleles and further demonstrated that base-editing hits could positively and negatively tune T cell cytotoxic function. Finally, higher-resolution screening using a base editor with relaxed protospacer-adjacent motif requirements9 (NG versus NGG) revealed specific structural domains and protein-protein interaction sites that can be targeted to tune T cell functions. Base-editing screens in primary immune cells thus provide biochemical insights with the potential to accelerate immunotherapy design.

PIM kinases regulate early human Th17 cell differentiation.

The serine/threonine-specific Moloney murine leukemia virus (PIM) kinase family (i.e., PIM1, PIM2, and PIM3) has been extensively studied in tumorigenesis. PIM kinases are downstream of several cytokine signaling pathways that drive immune-mediated diseases. Uncontrolled T helper 17 (Th17) cell activation has been associated with the pathogenesis of autoimmunity. However, the detailed molecular function of PIMs in human Th17 cell regulation has yet to be studied. In the present study, we comprehensively investigated how the three PIMs simultaneously alter transcriptional gene regulation during early human Th17 cell differentiation. By combining PIM triple knockdown with bulk and scRNA-seq approaches, we found that PIM deficiency promotes the early expression of key Th17-related genes while suppressing Th1-lineage genes. Further, PIMs modulate Th cell signaling, potentially via STAT1 and STAT3. Overall, our study highlights the inhibitory role of PIMs in human Th17 cell differentiation, thereby suggesting their association with autoimmune phenotypes.

The lncRNA Malat1 inhibits miR-15/16 to enhance cytotoxic T cell activation and memory cell formation.

Proper activation of cytotoxic T cells via the T cell receptor and the costimulatory receptor CD28 is essential for adaptive immunity against viruses, intracellular bacteria, and cancers. Through biochemical analysis of RNA:protein interactions, we uncovered a non-coding RNA circuit regulating activation and differentiation of cytotoxic T cells composed of the long non-coding RNA Malat1 (Metastasis Associated Lung Adenocarcinoma Transcript 1) and the microRNA family miR-15/16. miR-15/16 is a widely and highly expressed tumor suppressor miRNA family important for cell proliferation and survival. miR-15/16 play important roles in T cell responses to viral infection, including the regulation of antigen-specific T cell expansion and memory. Comparative Argonaute-2 high-throughput sequencing of crosslinking immunoprecipitation (AHC) combined with gene expression profiling in normal and miR-15/16-deficient mouse T cells revealed a large network of hundreds of direct miR-15/16 target mRNAs, many with functional relevance for T cell activation, survival and memory formation. Among these targets, Malat1 contained the largest absolute magnitude miR-15/16-dependent AHC peak. This binding site was among the strongest lncRNA:miRNA interactions detected in the T cell transcriptome. We used CRISPR targeting with homology directed repair to generate mice with a 5-nucleotide mutation in the miR-15/16-binding site in Malat1. This mutation interrupted Malat1:miR-15/16 interaction, and enhanced the repression of other miR-15/16 target genes, including CD28. Interrupting Malat1 interaction with miR-15/16 decreased cytotoxic T cell activation, including the expression of interleukin 2 (IL-2) and a broader CD28-responsive gene program. Accordingly, Malat1 mutation diminished memory cell persistence in mice following LCMV Armstrong and Listeria monocytogenes infection. This study marks a significant advance in the study of long non-coding RNAs in the immune system by ascribing cell-intrinsic, sequence-specific in vivo function to Malat1. These findings have implications for T cell-mediated autoimmune diseases, antiviral and anti-tumor immunity, as well as lung adenocarcinoma and other malignancies where Malat1 is overexpressed.

HIC1 interacts with FOXP3 multi protein complex: Novel pleiotropic mechanisms to regulate human regulatory T cell differentiation and function.

Transcriptional repressor, hypermethylated in cancer 1 (HIC1) participates in a range of important biological processes, such as tumor repression, immune suppression, embryonic development and epigenetic gene regulation. Further to these, we previously demonstrated that HIC1 provides a significant contribution to the function and development of regulatory T (Treg) cells. However, the mechanism by which it regulates these processes was not apparent. To address this question, we used affinity-purification mass spectrometry to characterize the HIC1 interactome in human Treg cells. Altogether 61 high-confidence interactors were identified, including IKZF3, which is a key transcription factor in the development of Treg cells. The biological processes associated with these interacting proteins include protein transport, mRNA processing, non-coding (ncRNA) transcription and RNA metabolism. The results revealed that HIC1 is part of a FOXP3-RUNX1-CBFB protein complex that regulates Treg signature genes thus improving our understanding of HIC1 function during early Treg cell differentiation.

Modular pooled discovery of synthetic knockin sequences to program durable cell therapies.

Chronic stimulation can cause T cell dysfunction and limit the efficacy of cellular immunotherapies. Improved methods are required to compare large numbers of synthetic knockin (KI) sequences to reprogram cell functions. Here, we developed modular pooled KI screening (ModPoKI), an adaptable platform for modular construction of DNA KI libraries using barcoded multicistronic adaptors. We built two ModPoKI libraries of 100 transcription factors (TFs) and 129 natural and synthetic surface receptors (SRs). Over 30 ModPoKI screens across human TCR- and CAR-T cells in diverse conditions identified a transcription factor AP4 (TFAP4) construct that enhanced fitness of chronically stimulated CAR-T cells and anti-cancer function in vitro and in vivo. ModPoKI's modularity allowed us to generate an ∼10,000-member library of TF combinations. Non-viral KI of a combined BATF-TFAP4 polycistronic construct enhanced fitness. Overexpressed BATF and TFAP4 co-occupy and regulate key gene targets to reprogram T cell function. ModPoKI facilitates the discovery of complex gene constructs to program cellular functions.

Multiplexed CRISPR gene editing in primary human islet cells with Cas9 ribonucleoprotein.

Successful genome editing in primary human islets could reveal features of the genetic regulatory landscape underlying ß cell function and diabetes risk. Here, we describe a CRISPR-based strategy to interrogate functions of predicted regulatory DNA elements using electroporation of a complex of Cas9 ribonucleoprotein (Cas9 RNP) and guide RNAs into primary human islet cells. We successfully targeted coding regions including the PDX1 exon 1, and non-coding DNA linked to diabetes susceptibility. CRISPR/Cas9 RNP approaches revealed genetic targets of regulation by DNA elements containing candidate diabetes risk SNPs, including an in vivo enhancer of the MPHOSPH9 gene. CRISPR/Cas9 RNP multiplexed targeting of two cis-regulatory elements linked to diabetes risk in PCSK1, which encodes an endoprotease crucial for insulin processing, also demonstrated efficient simultaneous editing of PCSK1 regulatory elements, resulting in impaired ß cell PCSK1 regulation and insulin secretion. Multiplex CRISPR/Cas9 RNP provides powerful approaches to investigate and elucidate human islet cell gene regulation in health and diabetes.

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.