Single-cell transcriptomic analyses reveal distinct immune cell contributions to epithelial barrier dysfunction in checkpoint inhibitor colitis.

Immune checkpoint inhibitor (ICI) therapy has revolutionized oncology, but treatments are limited by immune-related adverse events, including checkpoint inhibitor colitis (irColitis). Little is understood about the pathogenic mechanisms driving irColitis, which does not readily occur in model organisms, such as mice. To define molecular drivers of irColitis, we used single-cell multi-omics to profile approximately 300,000 cells from the colon mucosa and blood of 13 patients with cancer who developed irColitis (nine on anti-PD-1 or anti-CTLA-4 monotherapy and four on dual ICI therapy; most patients had skin or lung cancer), eight controls on ICI therapy and eight healthy controls. Patients with irColitis showed expanded mucosal Tregs, ITGAEHi CD8 tissue-resident memory T cells expressing CXCL13 and Th17 gene programs and recirculating ITGB2Hi CD8 T cells. Cytotoxic GNLYHi CD4 T cells, recirculating ITGB2Hi CD8 T cells and endothelial cells expressing hypoxia gene programs were further expanded in colitis associated with anti-PD-1/CTLA-4 therapy compared to anti-PD-1 therapy. Luminal epithelial cells in patients with irColitis expressed PCSK9, PD-L1 and interferon-induced signatures associated with apoptosis, increased cell turnover and malabsorption. Together, these data suggest roles for circulating T cells and epithelial-immune crosstalk critical to PD-1/CTLA-4-dependent tolerance and barrier function and identify potential therapeutic targets for irColitis.

Blood immunophenotyping identifies distinct kidney histopathology and outcomes in patients with lupus nephritis.

Lupus nephritis (LN) is a frequent manifestation of systemic lupus erythematosus, and fewer than half of patients achieve complete renal response with standard immunosuppressants. Identifying non-invasive, blood-based pathologic immune alterations associated with renal injury could aid therapeutic decisions. Here, we used mass cytometry immunophenotyping of peripheral blood mononuclear cells in 145 patients with biopsy-proven LN and 40 healthy controls to evaluate the heterogeneity of immune activation in patients with LN and to identify correlates of renal parameters and treatment response. Unbiased analysis identified 3 immunologically distinct groups of patients with LN that were associated with different patterns of histopathology, renal cell infiltrates, urine proteomic profiles, and treatment response at one year. Patients with enriched circulating granzyme B+ T cells at baseline showed more severe disease and increased numbers of activated CD8 T cells in the kidney, yet they had the highest likelihood of treatment response. A second group characterized primarily by a high type I interferon signature had a lower likelihood of response to therapy, while a third group appeared immunologically inactive by immunophenotyping at enrollment but with chronic renal injuries. Main immune profiles could be distilled down to 5 simple cytometric parameters that recapitulate several of the associations, highlighting the potential for blood immune profiling to translate to clinically useful non-invasive metrics to assess immune-mediated disease in LN.

Cytotoxic CD4+ T cells eliminate senescent cells by targeting cytomegalovirus antigen.

Senescent cell accumulation has been implicated in the pathogenesis of aging-associated diseases, including cancer. The mechanism that prevents the accumulation of senescent cells in aging human organs is unclear. Here, we demonstrate that a virus-immune axis controls the senescent fibroblast accumulation in the human skin. Senescent fibroblasts increased in old skin compared with young skin. However, they did not increase with advancing age in the elderly. Increased CXCL9 and cytotoxic CD4+ T cells (CD4 CTLs) recruitment were significantly associated with reduced senescent fibroblasts in the old skin. Senescent fibroblasts expressed human leukocyte antigen class II (HLA-II) and human cytomegalovirus glycoprotein B (HCMV-gB), becoming direct CD4 CTL targets. Skin-resident CD4 CTLs eliminated HCMV-gB+ senescent fibroblasts in an HLA-II-dependent manner, and HCMV-gB activated CD4 CTLs from the human skin. Collectively, our findings demonstrate HCMV reactivation in senescent cells, which CD4 CTLs can directly eliminate through the recognition of the HCMV-gB antigen.

Differentiation of exhausted CD8+ T cells after termination of chronic antigen stimulation stops short of achieving functional T cell memory.

T cell exhaustion is associated with failure to clear chronic infections and malignant cells. Defining the molecular mechanisms of T cell exhaustion and reinvigoration is essential to improving immunotherapeutic modalities. Here we confirmed pervasive phenotypic, functional and transcriptional differences between memory and exhausted antigen-specific CD8+ T cells in human hepatitis C virus (HCV) infection before and after treatment. After viral cure, phenotypic changes in clonally stable exhausted T cell populations suggested differentiation toward a memory-like profile. However, functionally, the cells showed little improvement, and critical transcriptional regulators remained in the exhaustion state. Notably, T cells from chronic HCV infection that were exposed to antigen for less time because of viral escape mutations were functionally and transcriptionally more similar to memory T cells from spontaneously resolved HCV infection. Thus, the duration of T cell stimulation impacts exhaustion recovery, with antigen removal after long-term exhaustion being insufficient for the development of functional T cell memory.

Genome-wide enhancer maps link risk variants to disease genes.

Genome-wide association studies (GWAS) have identified thousands of noncoding loci that are associated with human diseases and complex traits, each of which could reveal insights into the mechanisms of disease1. Many of the underlying causal variants may affect enhancers2,3, but we lack accurate maps of enhancers and their target genes to interpret such variants. We recently developed the activity-by-contact (ABC) model to predict which enhancers regulate which genes and validated the model using CRISPR perturbations in several cell types4. Here we apply this ABC model to create enhancer-gene maps in 131 human cell types and tissues, and use these maps to interpret the functions of GWAS variants. Across 72 diseases and complex traits, ABC links 5,036 GWAS signals to 2,249 unique genes, including a class of 577 genes that appear to influence multiple phenotypes through variants in enhancers that act in different cell types. In inflammatory bowel disease (IBD), causal variants are enriched in predicted enhancers by more than 20-fold in particular cell types such as dendritic cells, and ABC achieves higher precision than other regulatory methods at connecting noncoding variants to target genes. These variant-to-function maps reveal an enhancer that contains an IBD risk variant and that regulates the expression of PPIF to alter the membrane potential of mitochondria in macrophages. Our study reveals principles of genome regulation, identifies genes that affect IBD and provides a resource and generalizable strategy to connect risk variants of common diseases to their molecular and cellular functions.

Disabled Homolog 2 Controls Prometastatic Activity of Tumor-Associated Macrophages.

Tumor-associated macrophages (TAM) are regulators of extracellular matrix (ECM) remodeling and metastatic progression, the main cause of cancer-associated death. We found that disabled homolog 2 mitogen-responsive phosphoprotein (DAB2) is highly expressed in tumor-infiltrating TAMs and that its genetic ablation significantly impairs lung metastasis formation. DAB2-expressing TAMs, mainly localized along the tumor-invasive front, participate in integrin recycling, ECM remodeling, and directional migration in a tridimensional matrix. DAB2+ macrophages escort the invasive dissemination of cancer cells by a mechanosensing pathway requiring the transcription factor YAP. In human lobular breast and gastric carcinomas, DAB2+ TAMs correlated with a poor clinical outcome, identifying DAB2 as potential prognostic biomarker for stratification of patients with cancer. DAB2 is therefore central for the prometastatic activity of TAMs. SIGNIFICANCE: DAB2 expression in macrophages is essential for metastasis formation but not primary tumor growth. Mechanosensing cues, activating the complex YAP-TAZ, regulate DAB2 in macrophages, which in turn controls integrin recycling and ECM remodeling in 3-D tissue matrix. The presence of DAB2+ TAMs in patients with cancer correlates with worse prognosis.This article is highlighted in the In This Issue feature, p. 1611.

SLAMF6​ deficiency augments tumor killing and skews toward an effector phenotype revealing it as a novel T cell checkpoint.

SLAMF6 is a homotypic receptor of the Ig-superfamily whose exact role in immune modulation has remained elusive. Its constitutive expression on resting and activated T cells precludes it from being a bona fide exhaustion marker. By breeding Pmel-1 mice with SLAMF6 -/- mice, we generated donors for T cells lacking SLAMF6 and expressing a transgenic TCR for gp100-melanoma antigen. Activated Pmel-1xSLAMF6 -/- CD8+ T cells displayed improved polyfunctionality and strong tumor cytolysis. T-bet was the dominant transcription factor in Pmel-1 x SLAMF6 -/- cells, and upon activation, they acquired an effector-memory phenotype. Adoptive transfer of Pmel-1 x SLAMF6 -/- T cells to melanoma-bearing mice resulted in lasting tumor regression in contrast to temporary responses achieved with Pmel-1 T cells. LAG-3 expression was elevated in the SLAMF6 -/- cells, and the addition of the LAG-3-blocking antibody to the adoptive transfer protocol improved the SLAMF6 -/- T cells and expedited the antitumor response even further. The results from this study support the notion that SLAMF6 is an inhibitory immune receptor whose absence enables powerful CD8+ T cells to eradicate tumors.

The immune system helps to protect our bodies from illnesses and infections. Immunotherapies are medicines designed to treat diseases, such as cancer, by boosting the immune system against the condition. This is a powerful approach but so far immunotherapies have only had partial success and there is a need for further improvements. One protein called SLAMF6 is found on cells from the immune system that attack and kill cancer cells. Immunotherapies that suppress SLAMF6 on immune cells called killer T cells could increase immune system activity helping to treat cancers, particularly melanoma skin cancers. So far the potential for SLAMF6 as a target for immunotherapy has not been fully explored. Hajaj et al. created mice with killer T cells that recognized skin cancer cells and lacked SLAMF6. These modified cells were better at fighting cancer, producing more anti-cancer chemicals called cytokines and killing more cancer cells. The modified cells had a lasting effect on tumors and helped the mice to live longer. The effects could be further boosted by treating the mice in combination with other immunotherapies. SLAMF6 is a possible new target for skin cancer immunotherapy that could help more people to live longer following cancer diagnosis. The next step is to create a drug to target SLAMF6 in humans and to test it in clinical trials.

An immune-cell signature of bacterial sepsis.

Dysregulation of the immune response to bacterial infection can lead to sepsis, a condition with high mortality. Multiple whole-blood gene-expression studies have defined sepsis-associated molecular signatures, but have not resolved changes in transcriptional states of specific cell types. Here, we used single-cell RNA-sequencing to profile the blood of people with sepsis (n = 29) across three clinical cohorts with corresponding controls (n = 36). We profiled total peripheral blood mononuclear cells (PBMCs, 106,545 cells) and dendritic cells (19,806 cells) across all subjects and, on the basis of clustering of their gene-expression profiles, defined 16 immune-cell states. We identified a unique CD14+ monocyte state that is expanded in people with sepsis and validated its power in distinguishing these individuals from controls using public transcriptomic data from subjects with different disease etiologies and from multiple geographic locations (18 cohorts, n = 1,467 subjects). We identified a panel of surface markers for isolation and quantification of the monocyte state and characterized its epigenomic and functional phenotypes, and propose a model for its induction from human bone marrow. This study demonstrates the utility of single-cell genomics in discovering disease-associated cytologic signatures and provides insight into the cellular basis of immune dysregulation in bacterial sepsis.

A large peptidome dataset improves HLA class I epitope prediction across most of the human population.

Prediction of HLA epitopes is important for the development of cancer immunotherapies and vaccines. However, current prediction algorithms have limited predictive power, in part because they were not trained on high-quality epitope datasets covering a broad range of HLA alleles. To enable prediction of endogenous HLA class I-associated peptides across a large fraction of the human population, we used mass spectrometry to profile >185,000 peptides eluted from 95 HLA-A, -B, -C and -G mono-allelic cell lines. We identified canonical peptide motifs per HLA allele, unique and shared binding submotifs across alleles and distinct motifs associated with different peptide lengths. By integrating these data with transcript abundance and peptide processing, we developed HLAthena, providing allele-and-length-specific and pan-allele-pan-length prediction models for endogenous peptide presentation. These models predicted endogenous HLA class I-associated ligands with 1.5-fold improvement in positive predictive value compared with existing tools and correctly identified >75% of HLA-bound peptides that were observed experimentally in 11 patient-derived tumor cell lines.

The immune cell landscape in kidneys of patients with lupus nephritis.

Lupus nephritis is a potentially fatal autoimmune disease for which the current treatment is ineffective and often toxic. To develop mechanistic hypotheses of disease, we analyzed kidney samples from patients with lupus nephritis and from healthy control subjects using single-cell RNA sequencing. Our analysis revealed 21 subsets of leukocytes active in disease, including multiple populations of myeloid cells, T cells, natural killer cells and B cells that demonstrated both pro-inflammatory responses and inflammation-resolving responses. We found evidence of local activation of B cells correlated with an age-associated B-cell signature and evidence of progressive stages of monocyte differentiation within the kidney. A clear interferon response was observed in most cells. Two chemokine receptors, CXCR4 and CX3CR1, were broadly expressed, implying a potentially central role in cell trafficking. Gene expression of immune cells in urine and kidney was highly correlated, which would suggest that urine might serve as a surrogate for kidney biopsies.

Multiplexed enrichment and genomic profiling of peripheral blood cells reveal subset-specific immune signatures.

Specialized immune cell subsets are involved in autoimmune disease, cancer immunity, and infectious disease through a diverse range of functions mediated by overlapping pathways and signals. However, subset-specific responses may not be detectable in analyses of whole blood samples, and no efficient approach for profiling cell subsets at high throughput from small samples is available. We present a low-input microfluidic system for sorting immune cells into subsets and profiling their gene expression. We validate the system's technical performance against standard subset isolation and library construction protocols and demonstrate the importance of subset-specific profiling through in vitro stimulation experiments. We show the ability of this integrated platform to identify subset-specific disease signatures by profiling four immune cell subsets in blood from patients with systemic lupus erythematosus (SLE) and matched control subjects. The platform has the potential to make multiplexed subset-specific analysis routine in many research laboratories and clinical settings.

A secreted PD-L1 splice variant that covalently dimerizes and mediates immunosuppression.

Targeting immune checkpoint pathways, such as programmed death ligand-1 (PD-L1, also known as CD274 or B7-H1) or its receptor programmed cell death-1 (PD-1) has shown improved survival for patients with numerous types of cancers, not limited to lung cancer, melanoma, renal cell carcinoma, and Hodgkin lymphoma. PD-L1 is a co-inhibitory molecule whose expression on the surface of tumor cells is associated with worse prognosis in many tumors. Here we describe a splice variant (secPD-L1) that does not splice into the transmembrane domain, but instead produces a secreted form of PD-L1 that has a unique 18 amino acid tail containing a cysteine that allows it to homodimerize and more effectively inhibit lymphocyte function than monomeric soluble PD-L1. We show that recombinant secPD-L1 can dimerize and inhibit T-cell proliferation and IFN-gamma production in vitro. The secPD-L1 variant is expressed by malignant cells in vitro that also express high levels of full-length PD-L1. Transcriptomic analysis of gene expression across The Cancer Genome Atlas found the strongest association of secPD-L1 with full-length PD-L1, but also with subsets of immunologic genes, such as in myeloid-derived suppressor cells. Moreover, the splice variant is also expressed in normal tissues and within normal peripheral blood cells it is preferentially expressed in activated myeloid cells. This is the first report of a form of secreted PD-L1 that homodimerizes and is functionally active. SecPD-L1 may function as a paracrine negative immune regulator within the tumor, since secPD-L1 does not require a cell-to-cell interaction to mediate its inhibitory effect.

Methods for high-dimensional analysis of cells dissociated from cryopreserved synovial tissue

BACKGROUND: Detailed molecular analyses of cells from rheumatoid arthritis (RA) synovium hold promise in identifying cellular phenotypes that drive tissue pathology and joint damage. The Accelerating Medicines Partnership RA/SLE Network aims to deconstruct autoimmune pathology by examining cells within target tissues through multiple high-dimensional assays. Robust standardized protocols need to be developed before cellular phenotypes at a single cell level can be effectively compared across patient samples. METHODS: Multiple clinical sites collected cryopreserved synovial tissue fragments from arthroplasty and synovial biopsy in a 10% DMSO solution. Mechanical and enzymatic dissociation parameters were optimized for viable cell extraction and surface protein preservation for cell sorting and mass cytometry, as well as for reproducibility in RNA sequencing (RNA-seq). Cryopreserved synovial samples were collectively analyzed at a central processing site by a custom-designed and validated 35-marker mass cytometry panel. In parallel, each sample was flow sorted into fibroblast, T-cell, B-cell, and macrophage suspensions for bulk population RNA-seq and plate-based single-cell CEL-Seq2 RNA-seq. RESULTS: Upon dissociation, cryopreserved synovial tissue fragments yielded a high frequency of viable cells, comparable to samples undergoing immediate processing. Optimization of synovial tissue dissociation across six clinical collection sites with ~ 30 arthroplasty and ~ 20 biopsy samples yielded a consensus digestion protocol using 100 µg/ml of Liberase™ TL enzyme preparation. This protocol yielded immune and stromal cell lineages with preserved surface markers and minimized variability across replicate RNA-seq transcriptomes. Mass cytometry analysis of cells from cryopreserved synovium distinguished diverse fibroblast phenotypes, distinct populations of memory B cells and antibody-secreting cells, and multiple CD4+ and CD8+ T-cell activation states. Bulk RNA-seq of sorted cell populations demonstrated robust separation of synovial lymphocytes, fibroblasts, and macrophages. Single-cell RNA-seq produced transcriptomes of over 1000 genes/cell, including transcripts encoding characteristic lineage markers identified. CONCLUSIONS: We have established a robust protocol to acquire viable cells from cryopreserved synovial tissue with intact transcriptomes and cell surface phenotypes. A centralized pipeline to generate multiple high-dimensional analyses of synovial tissue samples collected across a collaborative network was developed. Integrated analysis of such datasets from large patient cohorts may help define molecular heterogeneity within RA pathology and identify new therapeutic targets and biomarkers.

Corrigendum: An immunogenic personal neoantigen vaccine for patients with melanoma.

This corrects the article DOI: 10.1038/nature22991.

An immunogenic personal neoantigen vaccine for patients with melanoma.

Effective anti-tumour immunity in humans has been associated with the presence of T cells directed at cancer neoantigens, a class of HLA-bound peptides that arise from tumour-specific mutations. They are highly immunogenic because they are not present in normal tissues and hence bypass central thymic tolerance. Although neoantigens were long-envisioned as optimal targets for an anti-tumour immune response, their systematic discovery and evaluation only became feasible with the recent availability of massively parallel sequencing for detection of all coding mutations within tumours, and of machine learning approaches to reliably predict those mutated peptides with high-affinity binding of autologous human leukocyte antigen (HLA) molecules. We hypothesized that vaccination with neoantigens can both expand pre-existing neoantigen-specific T-cell populations and induce a broader repertoire of new T-cell specificities in cancer patients, tipping the intra-tumoural balance in favour of enhanced tumour control. Here we demonstrate the feasibility, safety, and immunogenicity of a vaccine that targets up to 20 predicted personal tumour neoantigens. Vaccine-induced polyfunctional CD4+ and CD8+ T cells targeted 58 (60%) and 15 (16%) of the 97 unique neoantigens used across patients, respectively. These T cells discriminated mutated from wild-type antigens, and in some cases directly recognized autologous tumour. Of six vaccinated patients, four had no recurrence at 25 months after vaccination, while two with recurrent disease were subsequently treated with anti-PD-1 (anti-programmed cell death-1) therapy and experienced complete tumour regression, with expansion of the repertoire of neoantigen-specific T cells. These data provide a strong rationale for further development of this approach, alone and in combination with checkpoint blockade or other immunotherapies.

Mass Spectrometry Profiling of HLA-Associated Peptidomes in Mono-allelic Cells Enables More Accurate Epitope Prediction.

Identification of human leukocyte antigen (HLA)-bound peptides by liquid chromatography-tandem mass spectrometry (LC-MS/MS) is poised to provide a deep understanding of rules underlying antigen presentation. However, a key obstacle is the ambiguity that arises from the co-expression of multiple HLA alleles. Here, we have implemented a scalable mono-allelic strategy for profiling the HLA peptidome. By using cell lines expressing a single HLA allele, optimizing immunopurifications, and developing an application-specific spectral search algorithm, we identified thousands of peptides bound to 16 different HLA class I alleles. These data enabled the discovery of subdominant binding motifs and an integrative analysis quantifying the contribution of factors critical to epitope presentation, such as protein cleavage and gene expression. We trained neural-network prediction algorithms with our large dataset (>24,000 peptides) and outperformed algorithms trained on datasets of peptides with measured affinities. We thus demonstrate a strategy for systematically learning the rules of endogenous antigen presentation.

A Genome-wide CRISPR Screen in Primary Immune Cells to Dissect Regulatory Networks.

Finding the components of cellular circuits and determining their functions systematically remains a major challenge in mammalian cells. Here, we introduced genome-wide pooled CRISPR-Cas9 libraries into dendritic cells (DCs) to identify genes that control the induction of tumor necrosis factor (Tnf) by bacterial lipopolysaccharide (LPS), a key process in the host response to pathogens, mediated by the Tlr4 pathway. We found many of the known regulators of Tlr4 signaling, as well as dozens of previously unknown candidates that we validated. By measuring protein markers and mRNA profiles in DCs that are deficient in known or candidate genes, we classified the genes into three functional modules with distinct effects on the canonical responses to LPS and highlighted functions for the PAF complex and oligosaccharyltransferase (OST) complex. Our findings uncover new facets of innate immune circuits in primary cells and provide a genetic approach for dissection of mammalian cell circuits.

Immunogenetics. Dynamic profiling of the protein life cycle in response to pathogens.

Protein expression is regulated by the production and degradation of messenger RNAs (mRNAs) and proteins, but their specific relationships remain unknown. We combine measurements of protein production and degradation and mRNA dynamics so as to build a quantitative genomic model of the differential regulation of gene expression in lipopolysaccharide-stimulated mouse dendritic cells. Changes in mRNA abundance play a dominant role in determining most dynamic fold changes in protein levels. Conversely, the preexisting proteome of proteins performing basic cellular functions is remodeled primarily through changes in protein production or degradation, accounting for more than half of the absolute change in protein molecules in the cell. Thus, the proteome is regulated by transcriptional induction for newly activated cellular functions and by protein life-cycle changes for remodeling of preexisting functions.

CRISPR-Cas9 knockin mice for genome editing and cancer modeling.

CRISPR-Cas9 is a versatile genome editing technology for studying the functions of genetic elements. To broadly enable the application of Cas9 in vivo, we established a Cre-dependent Cas9 knockin mouse. We demonstrated in vivo as well as ex vivo genome editing using adeno-associated virus (AAV)-, lentivirus-, or particle-mediated delivery of guide RNA in neurons, immune cells, and endothelial cells. Using these mice, we simultaneously modeled the dynamics of KRAS, p53, and LKB1, the top three significantly mutated genes in lung adenocarcinoma. Delivery of a single AAV vector in the lung generated loss-of-function mutations in p53 and Lkb1, as well as homology-directed repair-mediated Kras(G12D) mutations, leading to macroscopic tumors of adenocarcinoma pathology. Together, these results suggest that Cas9 mice empower a wide range of biological and disease modeling applications.

A lentiviral RNAi library for human and mouse genes applied to an arrayed viral high-content screen.

To enable arrayed or pooled loss-of-function screens in a wide range of mammalian cell types, including primary and nondividing cells, we are developing lentiviral short hairpin RNA (shRNA) libraries targeting the human and murine genomes. The libraries currently contain 104,000 vectors, targeting each of 22,000 human and mouse genes with multiple sequence-verified constructs. To test the utility of the library for arrayed screens, we developed a screen based on high-content imaging to identify genes required for mitotic progression in human cancer cells and applied it to an arrayed set of 5,000 unique shRNA-expressing lentiviruses that target 1,028 human genes. The screen identified several known and approximately 100 candidate regulators of mitotic progression and proliferation; the availability of multiple shRNAs targeting the same gene facilitated functional validation of putative hits. This work provides a widely applicable resource for loss-of-function screens, as well as a roadmap for its application to biological discovery.