Droplet single-cell RNA-sequencing (dscRNA-seq) has enabled rapid, massively parallel profiling of transcriptomes. However, assessing differential expression across multiple individuals has been hampered by inefficient sample processing and technical batch effects. Here we describe a computational tool, demuxlet, that harnesses natural genetic variation to determine the sample identity of each cell and detect droplets containing two cells. These capabilities enable multiplexed dscRNA-seq experiments in which cells from unrelated individuals are pooled and captured at higher throughput than in standard workflows. Using simulated data, we show that 50 SNPs per cell are sufficient to assign 97% of singlets and identify 92% of doublets in pools of up to 64 individuals. Given genotyping data for each of 8 pooled samples, demuxlet correctly recovers the sample identity of >99% of singlets and identifies doublets at rates consistent with previous estimates. We apply demuxlet to assess cell type-specific changes in gene expression in 8 pooled lupus patient samples treated with IFN-β and perform eQTL analysis on 23 pooled samples.

T-cell genome engineering holds great promise for cell-based therapies for cancer, HIV, primary immune deficiencies, and autoimmune diseases, but genetic manipulation of human T cells has been challenging. Improved tools are needed to efficiently “knock out” genes and “knock in” targeted genome modifications to modulate T-cell function and correct disease-associated mutations. CRISPR/Cas9 technology is facilitating genome engineering in many cell types, but in human T cells its efficiency has been limited and it has not yet proven useful for targeted nucleotide replacements. Here we report efficient genome engineering in human CD4+ T cells using Cas9:single-guide RNA ribonucleoproteins (Cas9 RNPs). Cas9 RNPs allowed ablation of CXCR4, a coreceptor for HIV entry. Cas9 RNP electroporation caused up to ∼40% of cells to lose high-level cell-surface expression of CXCR4, and edited cells could be enriched by sorting based on low CXCR4 expression. Importantly, Cas9 RNPs paired with homology-directed repair template oligonucleotides generated a high frequency of targeted genome modifications in primary T cells. Targeted nucleotide replacement was achieved in CXCR4 and PD-1 (PDCD1), a regulator of T-cell exhaustion that is a validated target for tumor immunotherapy. Deep sequencing of a target site confirmed that Cas9 RNPs generated knock-in genome modifications with up to ∼20% efficiency, which accounted for up to approximately one-third of total editing events. These results establish Cas9 RNP technology for diverse experimental and therapeutic genome engineering applications in primary human T cells.

Immunotherapies with chimeric antigen receptor (CAR) T cells and checkpoint inhibitors (including antibodies that antagonize programmed cell death protein 1 [PD-1]) have both opened new avenues for cancer treatment, but the clinical potential of combined disruption of inhibitory checkpoints and CAR T cell therapy remains incompletely explored. Here we show that programmed death ligand 1 (PD-L1) expression on tumor cells can render human CAR T cells (anti-CD19 4-1BBζ) hypo-functional, resulting in impaired tumor clearance in a sub-cutaneous xenograft model. To overcome this suppressed anti-tumor response, we developed a protocol for combined Cas9 ribonucleoprotein (Cas9 RNP)-mediated gene editing and lentiviral transduction to generate PD-1 deficient anti-CD19 CAR T cells. Pdcd1 (PD-1) disruption augmented CAR T cell mediated killing of tumor cells in vitro and enhanced clearance of PD-L1+ tumor xenografts in vivo. This study demonstrates improved therapeutic efficacy of Cas9-edited CAR T cells and highlights the potential of precision genome engineering to enhance next-generation cell therapies.

The majority of genetic variants associated with common human diseases map to enhancers, non-coding elements that shape cell-type-specific transcriptional programs and responses to extracellular cues1–3. Systematic mapping of functional enhancers and their biological contexts is required to understand the mechanisms by which variation in non-coding genetic sequences contributes to disease. Functional enhancers can be mapped by genomic sequence disruption4–6, but this approach is limited to the subset of enhancers that are necessary in the particular cellular context being studied. We hypothesized that recruitment of a strong transcriptional activator to an enhancer would be sufficient to drive target gene expression, even if that enhancer was not currently active in the assayed cells. Here we describe a discovery platform that can identify stimulus-responsive enhancers for a target gene independent of stimulus exposure. We used tiled CRISPR activation (CRISPRa)7 to synthetically recruit a transcriptional activator to sites across large genomic regions (more than 100 kilobases) surrounding two key autoimmunity risk loci, CD69 and IL2RA. We identified several CRISPRa-responsive elements with chromatin features of stimulus-responsive enhancers, including an IL2RA enhancer that harbours an autoimmunity risk variant. Using engineered mouse models, we found that sequence perturbation of the disease-associated Il2ra enhancer did not entirely block Il2ra expression, but rather delayed the timing of gene activation in response to specific extracellular signals. Enhancer deletion skewed polarization of naive T cells towards a pro-inflammatory T helper (TH17) cell state and away from a regulatory T cell state. This integrated approach identifies functional enhancers and reveals how non-coding variation associated with human immune dysfunction alters context-specific gene programs.

Group 2 innate lymphoid cells (ILC2s) and CD4+ T helper type 2 (Th2) cells are defined by their similar effector cytokines, which together mediate the features of allergic immunity. Here, we show that tissue ILC2s and Th2 cells differentiate independently but share overlapping effector function programs mediated by exposure to the tissue-derived cytokines interleukin (IL)-25, IL-33, and thymic stromal lymphopoietin (TSLP). Loss of these three tissue signals does not affect lymph node priming but abrogates terminal differentiation of effector Th2 cells and adaptive lung inflammation in a T cell-intrinsic manner. Our findings suggest how diverse perturbations activate type 2 immunity, uncovering a shared local tissue-elicited checkpoint that can be exploited to control both innate and adaptive allergic inflammation.

Over 90% of genetic variants associated with complex human traits map to non-coding regions, but little is understood about how they modulate gene regulation in health and disease. One possible mechanism is that genetic variants affect the activity of one or more cis-regulatory elements leading to gene expression variation in specific cell types. To identify such cases, we analyzed ATAC-seq and RNA-seq profiles from stimulated primary CD4 T cells in up to 105 healthy donors. We found that regions of accessible chromatin (ATAC-peaks) are co-accessible at kilobase and megabase resolution, consistent with the three-dimensional chromatin organization measured by in situ Hi-C in T cells. Fifteen percent of genetic variants located within ATAC-peaks affected the accessibility of the corresponding peak (local-ATAC-QTLs). Local-ATAC-QTLs have the largest effects on co-accessible peaks, are associated with gene expression and are enriched for autoimmune disease variants. Our results provide insights into how natural genetic variants modulate cis-regulatory elements, in isolation or in concert, to influence gene expression.

Kv7 (KCNQ) voltage-gated potassium channels control excitability in the brain, heart, and ear. Calmodulin (CaM) is crucial for Kv7 function, but how this calcium sensor affects activity has remained unclear. Here, we present X-ray crystallographic analysis of CaM:Kv7.4 and CaM:Kv7.5 AB domain complexes that reveal an Apo/CaM clamp conformation and calcium binding preferences. These structures, combined with small-angle X-ray scattering, biochemical, and functional studies, establish a regulatory mechanism for Kv7 CaM modulation based on a common architecture in which a CaM C-lobe calcium-dependent switch releases a shared Apo/CaM clamp conformation. This C-lobe switch inhibits voltage-dependent activation of Kv7.4 and Kv7.5 but facilitates Kv7.1, demonstrating that mechanism is shared by Kv7 isoforms despite the different directions of CaM modulation. Our findings provide a unified framework for understanding how CaM controls different Kv7 isoforms and highlight the role of membrane proximal domains for controlling voltage-gated channel function. VIDEO ABSTRACT.

While genetic variants are known to be associated with overall gene abundance in stimulated immune cells, less is known about their effects on alternative isoform usage. By analyzing RNA-seq profiles of monocyte-derived dendritic cells from 243 individuals, we uncovered thousands of unannotated isoforms synthesized in response to influenza infection and type 1 interferon stimulation. We identified more than a thousand quantitative trait loci (QTLs) associated with alternate isoform usage (isoQTLs), many of which are independent of expression QTLs (eQTLs) for the same gene. Compared with eQTLs, isoQTLs are enriched for splice sites and untranslated regions, but depleted of sequences upstream of annotated transcription start sites. Both eQTLs and isoQTLs explain a significant proportion of the disease heritability attributed to common genetic variants. At the ERAP2 locus, we shed light on the function of the gene and how two frequent, highly differentiated haplotypes with intermediate frequencies could be maintained by balancing selection. At baseline and following type 1 interferon stimulation, the major haplotype is associated with low ERAP2 expression caused by nonsense-mediated decay, while the minor haplotype, known to increase Crohn's disease risk, is associated with high ERAP2 expression. In response to influenza infection, we found two uncharacterized isoforms expressed from the major haplotype, likely the result of multiple perfectly linked variants affecting the transcription and splicing at the locus. Thus, genetic variants at a single locus could modulate independent gene regulatory processes in innate immune responses and, in the case of ERAP2, may confer a historical fitness advantage in response to virus.