Cell Simulation as Cell Segmentation.

Single-cell spatial transcriptomics promises a highly detailed view of a cell's transcriptional state and microenvironment, yet inaccurate cell segmentation can render this data murky by misattributing large numbers of transcripts to nearby cells or conjuring nonexistent cells. We adopt methods from ab initio cell simulation to rapidly infer morphologically plausible cell boundaries that preserve cell type heterogeneity. Benchmarking applied to datasets generated by three commercial platforms show superior performance and computational efficiency of this approach compared with existing methods. We show that improved accuracy in cell segmentation aids greatly in detection of difficult to accurately segment tumor infiltrating immune cells such as neutrophils and T cells. Lastly, through improvements in our ability to delineate subsets of tumor infiltrating T cells, we show that CXCL13-expressing CD8+ T cells tend to be more closely associated with tumor cells than their CXCL13-negative counterparts in data generated from renal cell carcinoma patient samples. Proseg is available under at open source license at https://github.com/dcjones/proseg.

MR1-restricted T cell clonotypes are associated with "resistance" to Mycobacterium tuberculosis infection.

T cells are required for protective immunity against Mycobacterium tuberculosis. We recently described a cohort of Ugandan household contacts of tuberculosis cases who appear to "resist" M. tuberculosis infection (resisters; RSTRs) and showed that these individuals harbor IFN-γ-independent T cell responses to M. tuberculosis-specific peptide antigens. However, T cells also recognize nonprotein antigens via antigen-presenting systems that are independent of genetic background, known as donor-unrestricted T cells (DURTs). We used tetramer staining and flow cytometry to characterize the association between DURTs and "resistance" to M. tuberculosis infection. Peripheral blood frequencies of most DURT subsets were comparable between RSTRs and latently infected controls (LTBIs). However, we observed a 1.65-fold increase in frequency of MR1-restricted T (MR1T) cells among RSTRs in comparison with LTBIs. Single-cell RNA sequencing of 18,251 MR1T cells sorted from 8 donors revealed 5,150 clonotypes that expressed a common transcriptional program, the majority of which were private. Sequencing of the T cell receptor α/T cell receptor δ (TCRα/δ) repertoire revealed several DURT clonotypes were expanded among RSTRs, including 2 MR1T clonotypes that recognized mycobacteria-infected cells in a TCR-dependent manner. Overall, our data reveal unexpected donor-specific diversity in the TCR repertoire of human MR1T cells as well as associations between mycobacteria-reactive MR1T clonotypes and resistance to M. tuberculosis infection.

Multi-omic profiling reveals the endogenous and neoplastic responses to immunotherapies in cutaneous T cell lymphoma.

Cutaneous T cell lymphomas (CTCLs) are skin cancers with poor survival rates and limited treatments. While immunotherapies have shown some efficacy, the immunological consequences of administering immune-activating agents to CTCL patients have not been systematically characterized. We apply a suite of high-dimensional technologies to investigate the local, cellular, and systemic responses in CTCL patients receiving either mono- or combination anti-PD-1 plus interferon-gamma (IFN-γ) therapy. Neoplastic T cells display no evidence of activation after immunotherapy. IFN-γ induces muted endogenous immunological responses, while anti-PD-1 elicits broader changes, including increased abundance of CLA+CD39+ T cells. We develop an unbiased multi-omic profiling approach enabling discovery of immune modules stratifying patients. We identify an enrichment of activated regulatory CLA+CD39+ T cells in non-responders and activated cytotoxic CLA+CD39+ T cells in leukemic patients. Our results provide insights into the effects of immunotherapy in CTCL patients and a generalizable framework for multi-omic analysis of clinical trials.

HA-1-targeted T cell receptor (TCR) T cell therapy for recurrent leukemia after hematopoietic stem cell transplantation.

Relapse is the leading cause of death after allogeneic hematopoietic stem cell transplantation (HCT) for leukemia. T cells engineered by gene transfer to express T cell receptors (TCR; TCR-T) specific for hematopoietic-restricted minor histocompatibility (H) antigens may provide a potent selective anti-leukemic effect post-HCT. We conducted a phase I clinical trial employing a novel TCR-T product targeting the minor H antigen HA-1 to treat or consolidate treatment of persistent or recurrent leukemia and myeloid neoplasms. The primary objective was to evaluate the feasibility and safety of administration of HA-1 TCR-T post-HCT. CD8+ and CD4+ T cells expressing the HA-1 TCR and a CD8-co-receptor were successfully manufactured from HA-1 disparate HCT donors. One or more infusions of HA-1 TCR-T following lymphodepleting chemotherapy were administered to nine HCT recipients who had developed disease recurrence post-HCT. TCR-T cells expanded and persisted in vivo after adoptive transfer. No dose-limiting toxicities occurred. Although the study was not designed to assess efficacy, four patients achieved or maintained complete remissions following lymphodepletion and HA-1 TCR-T, with one ongoing at >2 years. Single-cell RNA sequencing of relapsing/progressive leukemia after TCR-T therapy identified upregulated molecules associated with T cell dysfunction or cancer cell survival. HA-1 TCR-T therapy appears feasible and safe and shows preliminary signals of efficacy. This clinical trial is registered at clinicaltrials.gov as NCT03326921.

Modeling T cell temporal response to cancer immunotherapy rationalizes development of combinatorial treatment protocols.

Successful immunotherapy relies on triggering complex responses involving T cell dynamics in tumors and the periphery. Characterizing these responses remains challenging using static human single-cell atlases or mouse models. To address this, we developed a framework for in vivo tracking of tumor-specific CD8+ T cells over time and at single-cell resolution. Our tools facilitate the modeling of gene program dynamics in the tumor microenvironment (TME) and the tumor-draining lymph node (tdLN). Using this approach, we characterize two modes of anti-programmed cell death protein 1 (PD-1) activity, decoupling induced differentiation of tumor-specific activated precursor cells from conventional type 1 dendritic cell (cDC1)-dependent proliferation and recruitment to the TME. We demonstrate that combining anti-PD-1 therapy with anti-4-1BB agonist enhances the recruitment and proliferation of activated precursors, resulting in tumor control. These data suggest that effective response to anti-PD-1 therapy is dependent on sufficient influx of activated precursor CD8+ cells to the TME and highlight the importance of understanding system-level dynamics in optimizing immunotherapies.

Circulating cancer-specific CD8 T cell frequency is associated with response to PD-1 blockade in Merkel cell carcinoma.

Understanding cancer immunobiology has been hampered by difficulty identifying cancer-specific T cells. Merkel cell polyomavirus (MCPyV) causes most Merkel cell carcinomas (MCCs). All patients with virus-driven MCC express MCPyV oncoproteins, facilitating identification of virus (cancer)-specific T cells. We studied MCPyV-specific T cells from 27 patients with MCC using MCPyV peptide-HLA-I multimers, 26-color flow cytometry, single-cell transcriptomics, and T cell receptor (TCR) sequencing. In a prospective clinical trial, higher circulating MCPyV-specific CD8 T cell frequency before anti-PD-1 treatment was strongly associated with 2-year recurrence-free survival (75% if detectable, 0% if undetectable, p = 0.0018; ClinicalTrial.gov: NCT02488759). Intratumorally, such T cells were typically present, but their frequency did not significantly associate with response. Circulating MCPyV-specific CD8 T cells had increased stem/memory and decreased exhaustion signatures relative to their intratumoral counterparts. These results suggest that cancer-specific CD8 T cells in the blood may play a role in anti-PD-1 responses. Thus, strategies that augment their number or mobilize them into tumors could improve outcomes.

Crosstalk between CD64+MHCII+ macrophages and CD4+ T cells drives joint pathology during chikungunya.

Communications between immune cells are essential to ensure appropriate coordination of their activities. Here, we observed the infiltration of activated macrophages into the joint-footpads of chikungunya virus (CHIKV)-infected animals. Large numbers of CD64+MHCII+ and CD64+MHCII- macrophages were present in the joint-footpad, preceded by the recruitment of their CD11b+Ly6C+ inflammatory monocyte precursors. Recruitment and differentiation of these myeloid subsets were dependent on CD4+ T cells and GM-CSF. Transcriptomic and gene ontology analyses of CD64+MHCII+ and CD64+MHCII- macrophages revealed 89 differentially expressed genes, including genes involved in T cell proliferation and differentiation pathways. Depletion of phagocytes, including CD64+MHCII+ macrophages, from CHIKV-infected mice reduced disease pathology, demonstrating that these cells play a pro-inflammatory role in CHIKV infection. Together, these results highlight the synergistic dynamics of immune cell crosstalk in driving CHIKV immunopathogenesis. This study provides new insights in the disease mechanism and offers opportunities for development of novel anti-CHIKV therapeutics.

Merkel cell polyomavirus-specific and CD39+CLA+ CD8 T cells as blood-based predictive biomarkers for PD-1 blockade in Merkel cell carcinoma.

Merkel cell carcinoma is a skin cancer often driven by Merkel cell polyomavirus (MCPyV) with high rates of response to anti-PD-1 therapy despite low mutational burden. MCPyV-specific CD8 T cells are implicated in anti-PD-1-associated immune responses and provide a means to directly study tumor-specific T cell responses to treatment. Using mass cytometry and combinatorial tetramer staining, we find that baseline frequencies of blood MCPyV-specific cells correlated with response and survival. Frequencies of these cells decrease markedly during response to therapy. Phenotypes of MCPyV-specific CD8 T cells have distinct expression patterns of CD39, cutaneous lymphocyte-associated antigen (CLA), and CD103. Correspondingly, overall bulk CD39+CLA+ CD8 T cell frequencies in blood correlate with MCPyV-specific cell frequencies and similarly predicted favorable clinical outcomes. Conversely, frequencies of CD39+CD103+ CD8 T cells are associated with tumor burden and worse outcomes. These cell subsets can be useful as biomarkers and to isolate blood-derived tumor-specific T cells.

In-depth analysis of human virus-specific CD8+ T cells delineates unique phenotypic signatures for T cell specificity prediction.

Following viral infection, the human immune system generates CD8+ T cell responses to virus antigens that differ in specificity, abundance, and phenotype. A characterization of virus-specific T cell responses allows one to assess infection history and to understand its contribution to protective immunity. Here, we perform in-depth profiling of CD8+ T cells binding to CMV-, EBV-, influenza-, and SARS-CoV-2-derived antigens in peripheral blood samples from 114 healthy donors and 55 cancer patients using high-dimensional mass cytometry and single-cell RNA sequencing. We analyze over 500 antigen-specific T cell responses across six different HLA alleles and observed unique phenotypes of T cells specific for antigens from different virus categories. Using machine learning, we extract phenotypic signatures of antigen-specific T cells, predict virus specificity for bulk CD8+ T cells, and validate these predictions, suggesting that machine learning can be used to accurately predict antigen specificity from T cell phenotypes.

mRNA vaccination boosts S-specific T cell memory and promotes expansion of CD45RAint TEMRA-like CD8+ T cells in COVID-19 recovered individuals.

SARS-CoV-2 infection and mRNA vaccination both elicit spike (S)-specific T cell responses. To analyze how T cell memory from prior infection influences T cell responses to vaccination, we evaluated functional T cell responses in naive and previously infected vaccine recipients. Pre-vaccine S-specific responses are predictive of subsequent CD8+ T cell vaccine-response magnitudes. Comparing baseline with post-vaccination TCRß repertoires, we observed large clonotypic expansions correlated with the frequency of spike-specific T cells. Epitope mapping the largest CD8+ T cell responses confirms that an HLA-A∗03:01 epitope was highly immunodominant. Peptide-MHC tetramer staining together with mass cytometry and single-cell sequencing permit detailed phenotyping and clonotypic tracking of these S-specific CD8+ T cells. Our results demonstrate that infection-induced S-specific CD8+ T cell memory plays a significant role in shaping the magnitude and clonal composition of the circulating T cell repertoire after vaccination, with mRNA vaccination promoting CD8+ memory T cells to a TEMRA-like phenotype.

Single-cell protein expression profiling resolves circulating and resident memory T cell diversity across tissues and infection contexts.

Memory CD8+ T cells can be broadly divided into circulating (TCIRCM) and tissue-resident memory T (TRM) populations. Despite well-defined migratory and transcriptional differences, the phenotypic and functional delineation of TCIRCM and TRM cells, particularly across tissues, remains elusive. Here, we utilized an antibody screening platform and machine learning prediction pipeline (InfinityFlow) to profile >200 proteins in TCIRCM and TRM cells in solid organs and barrier locations. High-dimensional analyses revealed unappreciated heterogeneity within TCIRCM and TRM cell lineages across nine different organs after either local or systemic murine infection models. Additionally, we demonstrated the relative effectiveness of strategies allowing for the selective ablation of TCIRCM or TRM populations across organs and identified CD55, KLRG1, CXCR6, and CD38 as stable markers for characterizing memory T cell function during inflammation. Together, these data and analytical framework provide an in-depth resource for memory T cell classification in both steady-state and inflammatory conditions.

An information theoretic approach to detecting spatially varying genes.

A key step in spatial transcriptomics is identifying genes with spatially varying expression patterns. We adopt an information theoretic perspective to this problem by equating the degree of spatial coherence with the Jensen-Shannon divergence between pairs of nearby cells and pairs of distant cells. To avoid the notoriously difficult problem of estimating information theoretic divergences, we use modern approximation techniques to implement a computationally efficient algorithm designed to scale with in situ spatial transcriptomics technologies. In addition to being highly scalable, we show that our method, which we call maximization of spatial information (Maxspin), improves accuracy across several spatial transcriptomics platforms and a variety of simulations when compared with a variety of state-of-the-art methods. To further demonstrate the method, we generated in situ spatial transcriptomics data in a renal cell carcinoma sample using the CosMx Spatial Molecular Imager and used Maxspin to reveal novel spatial patterns of tumor cell gene expression.

Persistent serum protein signatures define an inflammatory subcategory of long COVID.

Long COVID or post-acute sequelae of SARS-CoV-2 (PASC) is a clinical syndrome featuring diverse symptoms that can persist for months following acute SARS-CoV-2 infection. The aetiologies may include persistent inflammation, unresolved tissue damage or delayed clearance of viral protein or RNA, but the biological differences they represent are not fully understood. Here we evaluate the serum proteome in samples, longitudinally collected from 55 PASC individuals with symptoms lasting ≥60 days after onset of acute infection, in comparison to samples from symptomatically recovered SARS-CoV-2 infected and uninfected individuals. Our analysis indicates heterogeneity in PASC and identified subsets with distinct signatures of persistent inflammation. Type II interferon signaling and canonical NF-κB signaling (particularly associated with TNF), appear to be the most differentially enriched signaling pathways, distinguishing a group of patients characterized also by a persistent neutrophil activation signature. These findings help to clarify biological diversity within PASC, identify participants with molecular evidence of persistent inflammation, and highlight dominant pathways that may have diagnostic or therapeutic relevance, including a protein panel that we propose as having diagnostic utility for differentiating inflammatory and non-inflammatory PASC.

Distinct immune responses associated with vaccination status and protection outcomes after malaria challenge.

Understanding immune mechanisms that mediate malaria protection is critical for improving vaccine development. Vaccination with radiation-attenuated Plasmodium falciparum sporozoites (PfRAS) induces high level of sterilizing immunity against malaria and serves as a valuable tool for the study of protective mechanisms. To identify vaccine-induced and protection-associated responses during malarial infection, we performed transcriptome profiling of whole blood and in-depth cellular profiling of PBMCs from volunteers who received either PfRAS or noninfectious mosquito bites, followed by controlled human malaria infection (CHMI) challenge. In-depth single-cell profiling of cell subsets that respond to CHMI in mock-vaccinated individuals showed a predominantly inflammatory transcriptome response. Whole blood transcriptome analysis revealed that gene sets associated with type I and II interferon and NK cell responses were increased in prior to CHMI while T and B cell signatures were decreased as early as one day following CHMI in protected vaccinees. In contrast, non-protected vaccinees and mock-vaccinated individuals exhibited shared transcriptome changes after CHMI characterized by decreased innate cell signatures and inflammatory responses. Additionally, immunophenotyping data showed different induction profiles of vδ2+ γδ T cells, CD56+ CD8+ T effector memory (Tem) cells, and non-classical monocytes between protected vaccinees and individuals developing blood-stage parasitemia, following treatment and resolution of infection. Our data provide key insights in understanding immune mechanistic pathways of PfRAS-induced protection and infective CHMI. We demonstrate that vaccine-induced immune response is heterogenous between protected and non-protected vaccinees and that inducted-malaria protection by PfRAS is associated with early and rapid changes in interferon, NK cell and adaptive immune responses. Trial Registration: ClinicalTrials.gov NCT01994525.

CD4+ T cells expressing CX3CR1, GPR56, with variable CD57 are associated with cardiometabolic diseases in persons with HIV.

Persons with HIV (PWH) on long-term antiretroviral therapy (ART) have a higher incidence and prevalence of cardiometabolic diseases attributed, in part, to persistent inflammation despite viral suppression. In addition to traditional risk factors, immune responses to co-infections such as cytomegalovirus (CMV) may play an unappreciated role in cardiometabolic comorbidities and offer new potential therapeutic targets in a subgroup of individuals. We assessed the relationship of CX3CR1+, GPR56+, and CD57+/- T cells (termed CGC+) with comorbid conditions in a cohort of 134 PWH co-infected with CMV on long-term ART. We found that PWH with cardiometabolic diseases (non-alcoholic fatty liver disease, calcified coronary arteries, or diabetes) had higher circulating CGC+CD4+ T cells compared to metabolically healthy PWH. The traditional risk factor most correlated with CGC+CD4+ T cell frequency was fasting blood glucose, as well as starch/sucrose metabolites. While unstimulated CGC+CD4+ T cells, like other memory T cells, depend on oxidative phosphorylation for energy, they exhibited higher expression of carnitine palmitoyl transferase 1A compared to other CD4+ T cell subsets, suggesting a potentially greater capacity for fatty acid ß-oxidation. Lastly, we show that CMV-specific T cells against multiple viral epitopes are predominantly CGC+. Together, this study suggests that among PWH, CGC+ CD4+ T cells are frequently CMV-specific and are associated with diabetes, coronary arterial calcium, and non-alcoholic fatty liver disease. Future studies should assess whether anti-CMV therapies could reduce cardiometabolic disease risk in some individuals.

Intraventricular B7-H3 CAR T Cells for Diffuse Intrinsic Pontine Glioma: Preliminary First-in-Human Bioactivity and Safety.

Diffuse intrinsic pontine glioma (DIPG) remains a fatal brainstem tumor demanding innovative therapies. As B7-H3 (CD276) is expressed on central nervous system (CNS) tumors, we designed B7-H3-specific chimeric antigen receptor (CAR) T cells, confirmed their preclinical efficacy, and opened BrainChild-03 (NCT04185038), a first-in-human phase I trial administering repeated locoregional B7-H3 CAR T cells to children with recurrent/refractory CNS tumors and DIPG. Here, we report the results of the first three evaluable patients with DIPG (including two who enrolled after progression), who received 40 infusions with no dose-limiting toxicities. One patient had sustained clinical and radiographic improvement through 12 months on study. Patients exhibited correlative evidence of local immune activation and persistent cerebrospinal fluid (CSF) B7-H3 CAR T cells. Targeted mass spectrometry of CSF biospecimens revealed modulation of B7-H3 and critical immune analytes (CD14, CD163, CSF-1, CXCL13, and VCAM-1). Our data suggest the feasibility of repeated intracranial B7-H3 CAR T-cell dosing and that intracranial delivery may induce local immune activation. SIGNIFICANCE: This is the first report of repeatedly dosed intracranial B7-H3 CAR T cells for patients with DIPG and includes preliminary tolerability, the detection of CAR T cells in the CSF, CSF cytokine elevations supporting locoregional immune activation, and the feasibility of serial mass spectrometry from both serum and CSF. This article is highlighted in the In This Issue feature, p. 1.

Bystander CD4+ T cells infiltrate human tumors and are phenotypically distinct.

Tumor-specific T cells likely underpin effective immune checkpoint-blockade therapies. Yet, most studies focus on Treg cells and CD8+ tumor-infiltrating lymphocytes (TILs). Here, we study CD4+ TILs in human lung and colorectal cancers and observe that non-Treg CD4+ TILs average more than 70% of total CD4+ TILs in both cancer types. Leveraging high dimensional analyses including mass cytometry, we reveal that CD4+ TILs are phenotypically heterogeneous, within each tumor and across patients. Consistently, we find different subsets of CD4+ TILs showing characteristics of effectors, tissue resident memory (Trm) or exhausted cells (expressing PD-1, CTLA-4 and CD39). In both cancer types, the frequencies of CD39- non-Treg CD4+ TILs strongly correlate with frequencies of CD39- CD8+ TILs, which we and others have previously shown to be enriched for cells specific for cancer-unrelated antigens (bystanders). Ex-vivo, we demonstrate that CD39- CD4+ TILs can be specific for cancer-unrelated antigens, such as HCMV epitopes. Overall, our findings highlight that CD4+ TILs can also recognize cancer-unrelated antigens and suggest measuring CD39 expression as a straightforward way to quantify or isolate bystander CD4+ T cells.

Single-cell analysis reveals clonally expanded tumor-associated CD57+ CD8 T cells are enriched in the periphery of patients with metastatic urothelial cancer responding to PD-L1 blockade.

BACKGROUND: A growing body of evidence suggests that T-cell responses against neoantigens are critical regulators of response to immune checkpoint blockade. We previously showed that circulating neoantigen-specific CD8 T cells in patients with lung cancer responding to anti-Programmed death-ligand 1 (PD-L1) (atezolizumab) exhibit a unique phenotype with high expression of CD57, CD244, and KLRG1. Here, we extended our analysis on neoantigen-specific CD8 T cells to patients with metastatic urothelial cancer (mUC) and further profiled total CD8 T cells to identify blood-based predictive biomarkers of response to atezolizumab. METHODS: We identified tumor neoantigens from 20 patients with mUC and profiled their peripheral CD8 T cells using highly multiplexed combinatorial tetramer staining. Another set of patients with mUC treated with atezolizumab (n=30) or chemotherapy (n=40) were selected to profile peripheral CD8 T cells by mass cytometry. Using single-cell transcriptional analysis (single-cell RNA sequencing (scRNA-seq)), together with CITE-seq (cellular indexing of transcriptomes and epitopes by sequencing) and paired T-cell receptor (TCR) sequencing, we further characterized peripheral CD8 T cells in a subset of patients (n=16). RESULTS: High frequency of CD57 was observed in neoantigen-specific CD8 T cells in patients with mUC responding to atezolizumab. Extending these findings to bulk CD8 T cells, we found higher frequency of CD57 expressing CD8 T cells before treatment in patients responding to atezolizumab (n=20, p<0.01) but not to chemotherapy. These findings were corroborated in a validation cohort (n=30, p<0.01) and notably were independent of known biomarkers of response. scRNA-seq analysis identified a clonally expanded cluster enriched within CD57+ CD8 T cells in responding patients characterized by higher expression of genes associated with activation, cytotoxicity, and tissue-resident memory markers. Furthermore, compared with CD57- CD8 T cells, TCRs of CD57+ CD8 T cells showed increased overlap with the TCR repertoire of tumor-infiltrating T cells. CONCLUSIONS: Collectively, we show high frequencies of CD57 among neoantigen-specific and bulk CD8 T cells in patients responding to atezolizumab. The TCR repertoire overlap between peripheral CD57+ CD8 T cells and tumor-infiltrating lymphocytes suggest that accumulation of peripheral CD57+ CD8 T cells is reflective of an ongoing antitumor T-cell response. Our findings provide evidence and rationale for using circulating CD8 T cells expressing CD57 as a readily accessible blood-based biomarker for selecting patients with mUC for atezolizumab therapy.

Single-cell immunology of SARS-CoV-2 infection.

Gaining a better understanding of the immune cell subsets and molecular factors associated with protective or pathological immunity against severe acute respiratory syndrome coronavirus (SARS-CoV)-2 could aid the development of vaccines and therapeutics for coronavirus disease 2019 (COVID-19). Single-cell technologies, such as flow cytometry, mass cytometry, single-cell transcriptomics and single-cell multi-omic profiling, offer considerable promise in dissecting the heterogeneity of immune responses among individual cells and uncovering the molecular mechanisms of COVID-19 pathogenesis. Single-cell immune-profiling studies reported to date have identified innate and adaptive immune cell subsets that correlate with COVID-19 disease severity, as well as immunological factors and pathways of potential relevance to the development of vaccines and treatments for COVID-19. For facilitation of integrative studies and meta-analyses into the immunology of SARS-CoV-2 infection, we provide standardized, download-ready versions of 21 published single-cell sequencing datasets (over 3.2 million cells in total) as well as an interactive visualization portal for data exploration.