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.

Stromal cells (SCs) establish the compartmentalization of lymphoid tissues critical to the immune response. However, the full diversity of lymph node (LN) SCs remains undefined. Using droplet-based single-cell RNA sequencing, we identified nine peripheral LN non-endothelial SC clusters. Included are the established subsets, Ccl19 T-zone reticular cells (TRCs), marginal reticular cells, follicular dendritic cells (FDCs), and perivascular cells. We also identified Ccl19 TRCs, likely including cholesterol-25-hydroxylase cells located at the T-zone perimeter, Cxcl9 TRCs in the T-zone and interfollicular region, CD34 SCs in the capsule and medullary vessel adventitia, indolethylamine N-methyltransferase SCs in the medullary cords, and Nr4a1 SCs in several niches. These data help define how transcriptionally distinct LN SCs support niche-restricted immune functions and provide evidence that many SCs are in an activated state.

Summary Responses to anti-PD-1 immunotherapy occur but are infrequent in bladder cancer. The specific T cells that mediate tumor rejection are unknown. T cells from human bladder tumors and non-malignant tissue were assessed with single-cell RNA and paired T cell receptor (TCR) sequencing of 30,604 T cells from 7 patients. We find that the states and repertoires of CD8 + T cells are not distinct in tumors compared with non-malignant tissues. In contrast, single-cell analysis of CD4 + T cells demonstrates several tumor-specific states, including multiple distinct states of regulatory T cells. Surprisingly, we also find multiple cytotoxic CD4 + T cell states that are clonally expanded. These CD4 + T cells can kill autologous tumors in an MHC class II-dependent fashion and are suppressed by regulatory T cells. Further, a gene signature of cytotoxic CD4 + T cells in tumors predicts a clinical response in 244 metastatic bladder cancer patients treated with anti-PD-L1.