Scirpy: a Scanpy extension for analyzing single-cell T-cell receptor-sequencing data

Clonal landscape of autoantibody-secreting plasmablasts in COVID-19 patients

Whereas severe COVID-19 is often associated with elevated autoantibody titers, the underlying mechanism behind their generation has remained unclear. Here we report clonal composition and diversity of autoantibodies in humoral response to SARS-CoV-2. Immunoglobulin repertoire analysis and characterization of plasmablast-derived monoclonal antibodies uncovered clonal expansion of plasmablasts producing cardiolipin (CL)-reactive autoantibodies. Half of the expanded CL-reactive clones exhibited strong binding to SARS-CoV-2 antigens. One such clone, CoV1804, was reactive to both CL and viral nucleocapsid (N), and further showed anti-nucleolar activity in human cells. Notably, antibodies sharing genetic features with CoV1804 were identified in COVID-19 patient-derived immunoglobulins, thereby constituting a novel public antibody. These public autoantibodies had numerous mutations that unambiguously enhanced anti-N reactivity, when causing fluctuations in anti-CL reactivity along with the acquisition of additional self-reactivities, such as anti-nucleolar activity, in the progeny. Thus, potentially CL-reactive precursors may have developed multiple self-reactivities through clonal selection, expansion, and somatic hypermutation driven by viral antigens. Our results revealed the nature of autoantibody production during COVID-19 and provided novel insights into the origin of virus-induced autoantibodies.

Neoself-antigens are the primary target for autoreactive T cells in human lupus

Major histocompatibility complex class II (MHC-II) is the most significant genetic risk factor for systemic lupus erythematosus (SLE), but the nature of the self-antigens that trigger autoimmunity remains unclear. Unusual self-antigens, termed neoself-antigens, are presented on MHC-II in the absence of the invariant chain essential for peptide presentation. Here, we demonstrate that neoself-antigens are the primary target for autoreactive T cells clonally expanded in SLE. When neoself-antigen presentation was induced by deleting the invariant chain in adult mice, neoself-reactive T cells were clonally expanded, leading to the development of lupus-like disease. Furthermore, we found that neoself-reactive CD4+ T cells were significantly expanded in SLE patients. A high frequency of Epstein-Barr virus reactivation is a risk factor for SLE. Neoself-reactive lupus T cells were activated by Epstein-Barr-virus-reactivated cells through downregulation of the invariant chain. Together, our findings imply that neoself-antigen presentation by MHC-II plays a crucial role in the pathogenesis of SLE.

ThymoSpheres culture: A model to study human polyclonal unconventional T cells

In vitro cultures remain crucial for studying the fundamental mechanisms of human T-cell development. Here, we introduce a novel in vitro cultivation system based on ThymoSpheres (TS): dense spheroids consisting of DLL4-expressing stromal cells and human hematopoietic precursor cells, in the absence of thymic epithelial cells. These spheroids are subsequently cultured at the air-liquid interphase. TS generate large numbers of mature T cells, are easy to manipulate, scalable, and can be repeatably sampled to monitor T-cell differentiation. The mature T cells generated from primary human hematopoietic precursor cells were extensively characterized using single-cell RNA and combined T-cell receptor (TCR) sequencing. These predominantly CD8α T cells exhibit transcriptional and TCR CDR3 characteristics similar to the recently described human polyclonal αβ unconventional T cell (UTC) lineage. This includes the expression of hallmark genes associated with agonist selection, such as IKZF2 (Helios), and the expression of various natural killer receptors. The TCR repertoire of these UTCs is polyclonal and enriched for CDR3-associated autoreactive features and early rearrangements of the TCR-α chain. In conclusion, TS cultures offer an intriguing platform to study the development of this human polyclonal UTC lineage and its inducing selection mechanisms.

Identifying ADGRG1 as a specific marker for tumor-reactive T cells in acute myeloid leukemia

Besides chemotherapy and hematopoietic stem cell transplantation (HSCT), autologous T cells can also serve as a new treatment approach for AML patients. However, the features of tumor-reactive T cells and their distinctive markers still lack full description. To evaluate the characteristics of tumor-reactive T cells, we collected bone marrow (BM) T cells from newly diagnosed AML patients with RUNX1::RUNX1T1 as examples for paired single-cell RNA sequencing and single-cell V(D)J sequencing. Based on the STARTRAC-like algorithm, we defined bystander T cells and tumor-reactive T cells. Compared with bystander T cells, tumor-reactive T cells presented as senescent-like cytotoxic terminally differentiated T cells (Temra) with upregulated NK-related markers. Additionally, we found ADGRG1 could serve as the specific marker of CD8+ T tumor-reactive T cell and validated it through the Runx1Runx1t1/+; Mx1-Cre mouse model. In chimeric antigen receptor (CAR)-T and target cell system, ADGRG1 was selectively upregulated upon antigen-TCR encounter. Moreover, ADGRG1+CD8+ T cells released a higher level of IFN-γ and showed higher cell-killing ability when exposed to matched AML blasts. Together, our findings depict the single-cell profile of tumor-reactive T cells in AML BM and propose that ADGRG1 can act as an indicator of T cell tumor reactivity in AML, which may be further harnessed for adoptive cell therapy and tumor-reactive TCR enrichment.

Neoadjuvant PARPi or chemotherapy in ovarian cancer informs targeting effector Treg cells for homologous-recombination-deficient tumors

Homologous recombination deficiency (HRD) is prevalent in cancer, sensitizing tumor cells to poly (ADP-ribose) polymerase (PARP) inhibition. However, the impact of HRD and related therapies on the tumor microenvironment (TME) remains elusive. Our study generates single-cell gene expression and T cell receptor profiles, along with validatory multimodal datasets from >100 high-grade serous ovarian cancer (HGSOC) samples, primarily from a phase II clinical trial (NCT04507841). Neoadjuvant monotherapy with the PARP inhibitor (PARPi) niraparib achieves impressive 62.5% and 73.6% response rates per RECIST v.1.1 and GCIG CA125, respectively. We identify effector regulatory T cells (eTregs) as key responders to HRD and neoadjuvant therapies, co-occurring with other tumor-reactive T cells, particularly terminally exhausted CD8+ T cells (Tex). TME-wide interferon signaling correlates with cancer cells upregulating MHC class II and co-inhibitory ligands, potentially driving Treg and Tex fates. Depleting eTregs in HRD mouse models, with or without PARP inhibition, significantly suppresses tumor growth without observable toxicities, underscoring the potential of eTreg-focused therapeutics for HGSOC and other HRD-related tumors.

Neoantigen-specific cytotoxic Tr1 CD4 T cells suppress cancer immunotherapy

CD4+ T cells can either enhance or inhibit tumour immunity. Although regulatory T cells have long been known to impede antitumour responses1-5, other CD4+ T cells have recently been implicated in inhibiting this response6,7. Yet, the nature and function of the latter remain unclear. Here, using vaccines containing MHC class I (MHC-I) neoantigens (neoAgs) and different doses of tumour-derived MHC-II neoAgs, we discovered that whereas the inclusion of vaccines with low doses of MHC-II-restricted peptides (LDVax) promoted tumour rejection, vaccines containing high doses of the same MHC-II neoAgs (HDVax) inhibited rejection. Characterization of the inhibitory cells induced by HDVax identified them as type 1 regulatory T (Tr1) cells expressing IL-10, granzyme B, perforin, CCL5 and LILRB4. Tumour-specific Tr1 cells suppressed tumour rejection induced by anti-PD1, LDVax or adoptively transferred tumour-specific effector T cells. Mechanistically, HDVax-induced Tr1 cells selectively killed MHC-II tumour antigen-presenting type 1 conventional dendritic cells (cDC1s), leading to low numbers of cDC1s in tumours. We then documented modalities to overcome this inhibition, specifically via anti-LILRB4 blockade, using a CD8-directed IL-2 mutein, or targeted loss of cDC2/monocytes. Collectively, these data show that cytotoxic Tr1 cells, which maintain peripheral tolerance, also inhibit antitumour responses and thereby function to impede immune control of cancer.

Comparison of uridine and N1-methylpseudouridine mRNA platforms in development of an Andes virus vaccine

The rodent-borne Andes virus (ANDV) causes a severe disease in humans. We developed an ANDV mRNA vaccine based on the M segment of the viral genome, either with regular uridine (U-mRNA) or N1-methylpseudouridine (m1Ψ-mRNA). Female mice immunized by m1Ψ-mRNA developed slightly greater germinal center (GC) responses than U-mRNA-immunized mice. Single cell RNA and BCR sequencing of the GC B cells revealed similar levels of activation, except an additional cluster of cells exhibiting interferon response in animals vaccinated with U-mRNA but not m1Ψ-mRNA. Similar immunoglobulin class-switching and somatic hypermutations were observed in response to the vaccines. Female Syrian hamsters were immunized via a prime-boost regimen with two doses of each vaccine. The titers of glycoprotein-binding antibodies were greater for U-mRNA construct than for m1Ψ-mRNA construct; however, the titers of ANDV-neutralizing antibodies were similar. Vaccinated animals were challenged with a lethal dose of ANDV, along with a naïve control group. All control animals and two animals vaccinated with a lower dose of m1Ψ-mRNA succumbed to infection whereas other vaccinated animals survived without evidence of virus replication. The data demonstrate the development of a protective vaccine against ANDV and the lack of a substantial effect of m1Ψ modification on immunogenicity and protection in rodents.

Cytotoxic Lymphocyte-Monocyte Complex Reflects the Dynamics of Coronavirus Disease 2019 Systemic Immune Response

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection causes a variety of clinical manifestations, many of which originate from altered immune responses, either locally or systemically. Immune cell cross-talk occurs mainly in lymphoid organs. However, systemic cell interaction specific to coronavirus disease 2019 has not been well characterized. Here, by employing single-cell RNA sequencing and imaging flow cytometry analysis, we unraveled, in peripheral blood, a heterogeneous group of cell complexes formed by the adherence of CD14+ monocytes to different cytotoxic lymphocytes, including SARS-CoV-2-specific CD8+ T cells, γδ T cells, and natural killer T cells. These lymphocytes attached to CD14+ monocytes that showed enhanced inflammasome activation and pyroptosis-induced cell death in progression stage; in contrast, in the convalescent phase, CD14+ monocytes with elevated antigen presentation potential were targeted by cytotoxic lymphocytes, thereby restricting the excessive immune activation. Collectively, our study reports previously unrecognized cell-cell interplay in the SARS-CoV-2-specific immune response, providing new insight into the intricacy of dynamic immune cell interaction representing antiviral defense.

Computational estimation of clonal diversity in autoimmunity

Diversity is the cornerstone of the adaptive immune system, crucial for its effectiveness against constantly evolving pathogens that pose threats to higher vertebrates. Accurately measuring and interpreting this diversity presents challenges for immunologists, as changes in diversity and clonotype composition can tip the balance between protective immunity and autoimmunity. In this review, we present the current methods commonly used to measure diversity from single-cell T-cell receptor and B-cell receptor sequencing. We also discuss two case studies where single-cell sequencing and diversity estimations have led to breakthroughs in autoimmune disease discovery and therapeutic innovation, and reflect upon the necessity and importance of accurately defining and measuring lymphocyte diversity in these contexts.

Panpipes: a pipeline for multiomic single-cell and spatial transcriptomic data analysis

Single-cell multiomic analysis of the epigenome, transcriptome, and proteome allows for comprehensive characterization of the molecular circuitry that underpins cell identity and state. However, the holistic interpretation of such datasets presents a challenge given a paucity of approaches for systematic, joint evaluation of different modalities. Here, we present Panpipes, a set of computational workflows designed to automate multimodal single-cell and spatial transcriptomic analyses by incorporating widely-used Python-based tools to perform quality control, preprocessing, integration, clustering, and reference mapping at scale. Panpipes allows reliable and customizable analysis and evaluation of individual and integrated modalities, thereby empowering decision-making before downstream investigations.

Multi-modal generative modeling for joint analysis of single-cell T cell receptor and gene expression data

Recent advances in single-cell immune profiling have enabled the simultaneous measurement of transcriptome and T cell receptor (TCR) sequences, offering great potential for studying immune responses at the cellular level. However, integrating these diverse modalities across datasets is challenging due to their unique data characteristics and technical variations. Here, to address this, we develop the multimodal generative model mvTCR to fuse modality-specific information across transcriptome and TCR into a shared representation. Our analysis demonstrates the added value of multimodal over unimodal approaches to capture antigen specificity. Notably, we use mvTCR to distinguish T cell subpopulations binding to SARS-CoV-2 antigens from bystander cells. Furthermore, when combined with reference mapping approaches, mvTCR can map newly generated datasets to extensive T cell references, facilitating knowledge transfer. In summary, we envision mvTCR to enable a scalable analysis of multimodal immune profiling data and advance our understanding of immune responses.

nf-core/airrflow: An adaptive immune receptor repertoire analysis workflow employing the Immcantation framework

Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) is a valuable experimental tool to study the immune state in health and following immune challenges such as infectious diseases, (auto)immune diseases, and cancer. Several tools have been developed to reconstruct B cell and T cell receptor sequences from AIRR-seq data and infer B and T cell clonal relationships. However, currently available tools offer limited parallelization across samples, scalability or portability to high-performance computing infrastructures. To address this need, we developed nf-core/airrflow, an end-to-end bulk and single-cell AIRR-seq processing workflow which integrates the Immcantation Framework following BCR and TCR sequencing data analysis best practices. The Immcantation Framework is a comprehensive toolset, which allows the processing of bulk and single-cell AIRR-seq data from raw read processing to clonal inference. nf-core/airrflow is written in Nextflow and is part of the nf-core project, which collects community contributed and curated Nextflow workflows for a wide variety of analysis tasks. We assessed the performance of nf-core/airrflow on simulated sequencing data with sequencing errors and show example results with real datasets. To demonstrate the applicability of nf-core/airrflow to the high-throughput processing of large AIRR-seq datasets, we validated and extended previously reported findings of convergent antibody responses to SARS-CoV-2 by analyzing 97 COVID-19 infected individuals and 99 healthy controls, including a mixture of bulk and single-cell sequencing datasets. Using this dataset, we extended the convergence findings to 20 additional subjects, highlighting the applicability of nf-core/airrflow to validate findings in small in-house cohorts with reanalysis of large publicly available AIRR datasets.

Human SARS-CoV-2 challenge uncovers local and systemic response dynamics

The COVID-19 pandemic is an ongoing global health threat, yet our understanding of the dynamics of early cellular responses to this disease remains limited1. Here in our SARS-CoV-2 human challenge study, we used single-cell multi-omics profiling of nasopharyngeal swabs and blood to temporally resolve abortive, transient and sustained infections in seronegative individuals challenged with pre-Alpha SARS-CoV-2. Our analyses revealed rapid changes in cell-type proportions and dozens of highly dynamic cellular response states in epithelial and immune cells associated with specific time points and infection status. We observed that the interferon response in blood preceded the nasopharyngeal response. Moreover, nasopharyngeal immune infiltration occurred early in samples from individuals with only transient infection and later in samples from individuals with sustained infection. High expression of HLA-DQA2 before inoculation was associated with preventing sustained infection. Ciliated cells showed multiple immune responses and were most permissive for viral replication, whereas nasopharyngeal T cells and macrophages were infected non-productively. We resolved 54 T cell states, including acutely activated T cells that clonally expanded while carrying convergent SARS-CoV-2 motifs. Our new computational pipeline Cell2TCR identifies activated antigen-responding T cells based on a gene expression signature and clusters these into clonotype groups and motifs. Overall, our detailed time series data can serve as a Rosetta stone for epithelial and immune cell responses and reveals early dynamic responses associated with protection against infection.

Targeting T cell plasticity in kidney and gut inflammation by pooled single-cell CRISPR screening

Pro-inflammatory CD4+ T cells are major drivers of autoimmune diseases, yet therapies modulating T cell phenotypes to promote an anti-inflammatory state are lacking. Here, we identify T helper 17 (TH17) cell plasticity in the kidneys of patients with antineutrophil cytoplasmic antibody-associated glomerulonephritis on the basis of single-cell (sc) T cell receptor analysis and scRNA velocity. To uncover molecules driving T cell polarization and plasticity, we established an in vivo pooled scCRISPR droplet sequencing (iCROP-seq) screen and applied it to mouse models of glomerulonephritis and colitis. CRISPR-based gene targeting in TH17 cells could be ranked according to the resulting transcriptional perturbations, and polarization biases into T helper 1 (TH1) and regulatory T cells could be quantified. Furthermore, we show that iCROP-seq can facilitate the identification of therapeutic targets by efficient functional stratification of genes and pathways in a disease- and tissue-specific manner. These findings uncover TH17 to TH1 cell plasticity in the human kidney in the context of renal autoimmunity.

CXCR6+CD69+ CD8+ T cells in ascites are associated with disease severity in patients with cirrhosis

Background & Aims

Patients with advanced cirrhosis often develop hepatic decompensation, which is accompanied by systemic inflammation and may eventually lead to acute-on-chronic liver failure. One important cause of systemic hyperinflammation is a dysregulated overshooting immune response in ascites in the abdominal cavity. In this study, we analyzed the role of CD8+ T cells in the ascites immune compartment.

Methods

Peripheral blood and ascites fluid were collected from 50 patients with decompensated cirrhosis. Phenotype and functional responses of CD8+ T cells were analyzed, and obtained data were compared with each other as well as with healthy controls and patients with compensated cirrhosis.

Results

High-dimensional flow cytometry revealed that CD8+ T cells are abundant in the ascites of patients with cirrhosis and exhibit a chronically activated bystander phenotype with innate-like functions. Indeed, we identified distinct CXCR6+CD69+ clusters of late effector memory CD8+ T cells that were rarely found in blood and correlated with clinical parameters of disease severity. Moreover, this CD8+ T-cell population was hyperresponsive to innate cytokines and exhibited cytokine-mediated bystander activation. Interestingly, the Janus kinase (JAK) inhibitor tofacitinib was able to effectively block bystander-activated CXCR6+CD69+ CD8+ T cells and significantly suppress effector molecule production.

Conclusions

The results indicate that CXCR6+CD69+ CD8+ T cells in ascites are associated with disease severity and may contribute to inflammation in patients with decompensated cirrhosis, suggesting that targeted inhibition of this immune cell subset may be a viable therapeutic option.

Impact and Implications

Patients with advanced cirrhosis often develop hepatic decompensation, which is accompanied by systemic inflammation and eventually leads to acute-on-chronic liver failure. One important cause of systemic hyperinflammation is a dysregulated overshooting immune response in ascites in the abdominal cavity. In this study, we demonstrate that CXCR6+CD69+ CD8+ T cells are abundant in the ascites of patients with cirrhosis, exhibit a chronically activated bystander phenotype, and correlate with clinical parameters of disease severity. Moreover, we show that the Janus kinase (JAK) inhibitor tofacitinib can effectively block these bystander-activated CXCR6+CD69+ CD8+ T cells, suggesting that targeted inhibition of this immune cell subset may be a potential therapeutic strategy.

Clinical trial number

Prospective registry: INFEKTA (DRKS00010664).

Single-cell RNA sequencing reveals the expansion of circulating tissue-homing B cell subsets in secondary acute dengue viral infection

The roles of antibodies secreted by subsets of B cells in dengue virus (DENV) infection have been extensively studied, yet, the contribution of tissue-homing B cells to antiviral immunity remains unclear. In this study, we performed a comprehensive analysis of B cell subpopulations in peripheral blood samples from DENV-infected patients using single-cell RNA-sequencing (scRNA-seq) datasets and flow cytometry. We showed that plasma cells (PCs) and plasmablasts (PBs) were the predominant B cell populations during the acute phase of secondary natural DENV infection, but not in convalescent phase nor in healthy controls. Interestingly, these cells expressed proliferation, adhesion, and tissue-homing genes, including SELPLG, a homing marker of the skin, the initial infected site of DENV. Flow cytometry analysis confirmed a significant upregulation of cell surface expression of a cutaneous lymphocyte-associated antigen (CLA) encoded by SELPLG in PCs and PBs, compared to naive and memory B cells from the same patients. The analysis of an independent single-cell B-cell receptor sequencing (scBCR-seq) dataset of DENV-infected patients revealed that the peripheral blood PCs and PBs exhibited the highest clonal expansion in secondary DENV infection compared to other B cell subsets. These clonally expanded cells also expressed the highest levels of tissue-homing genes, including SELPLG. In addition, by utilizing a public scRNA-seq dataset of SARS-CoV2 infection, we demonstrated the upregulation of several tissue-homing genes in PCs and PBs. Our study provides evidence for the potential roles of tissue-homing B cell subsets in the context of immune responses against viral infections in humans.

Invariant γδTCR natural killer-like effector T cells in the naked mole-rat

The naked mole-rat (Heterocephalus glaber) is a long-lived rodent species showing resistance to the development of cancer. Although naked mole-rats have been reported to lack natural killer (NK) cells, γδ T cell-based immunity has been suggested in this species, which could represent an important arm of the immune system for antitumor responses. Here, we investigate the biology of these unconventional T cells in peripheral tissues (blood, spleen) and thymus of the naked mole-rat at different ages by TCR repertoire profiling and single-cell gene expression analysis. Using our own TCR annotation in the naked mole-rat genome, we report that the γδ TCR repertoire is dominated by a public invariant Vγ4-2/Vδ1-4 TCR, containing the complementary-determining-region-3 (CDR3)γ CTYWDSNYAKKLF / CDR3δ CALWELRTGGITAQLVF that are likely generated by short-homology-repeat-driven DNA rearrangements. This invariant TCR is specifically found in γδ T cells expressing genes associated with NK cytotoxicity and is generated in both the thoracic and cervical thymus of the naked mole-rat until adult life. Our results indicate that invariant Vγ4-2/Vδ1-4 NK-like effector T cells in the naked mole-rat can contribute to tumor immunosurveillance by γδ TCR-mediated recognition of a common molecular signal.

Single-cell immune repertoire analysis

Single-cell T cell and B cell antigen receptor-sequencing data analysis can potentially perform in-depth assessments of adaptive immune cells that inform on understanding immune cell development to tracking clonal expansion in disease and therapy. However, it has been extremely challenging to analyze and interpret T cells and B cells and their adaptive immune receptor repertoires at the single-cell level due to not only the complexity of the data but also the underlying biology. In this Review, we delve into the computational breakthroughs that have transformed the analysis of single-cell T cell and B cell antigen receptor-sequencing data.

Deep learning predictions of TCR-epitope interactions reveal epitope-specific chains in dual alpha T cells

T cells have the ability to eliminate infected and cancer cells and play an essential role in cancer immunotherapy. T cell activation is elicited by the binding of the T cell receptor (TCR) to epitopes displayed on MHC molecules, and the TCR specificity is determined by the sequence of its α and β chains. Here, we collect and curate a dataset of 17,715 αβTCRs interacting with dozens of class I and class II epitopes. We use this curated data to develop MixTCRpred, an epitope-specific TCR-epitope interaction predictor. MixTCRpred accurately predicts TCRs recognizing several viral and cancer epitopes. MixTCRpred further provides a useful quality control tool for multiplexed single-cell TCR sequencing assays of epitope-specific T cells and pinpoints a substantial fraction of putative contaminants in public databases. Analysis of epitope-specific dual α T cells demonstrates that MixTCRpred can identify α chains mediating epitope recognition. Applying MixTCRpred to TCR repertoires from COVID-19 patients reveals enrichment of clonotypes predicted to bind an immunodominant SARS-CoV-2 epitope. Overall, MixTCRpred provides a robust tool to predict TCRs interacting with specific epitopes and interpret TCR-sequencing data from both bulk and epitope-specific T cells.

Evolution of T cells in the cancer-resistant naked mole-rat

Naked mole-rats (NMRs) are best known for their extreme longevity and cancer resistance, suggesting that their immune system might have evolved to facilitate these phenotypes. Natural killer (NK) and T cells have evolved to detect and destroy cells infected with pathogens and to provide an early response to malignancies. While it is known that NMRs lack NK cells, likely lost during evolution, little is known about their T-cell subsets in terms of the evolution of the genes that regulate their function, their clonotypic diversity, and the thymus where they mature. Here we find, using single-cell transcriptomics, that NMRs have a large circulating population of γδT cells, which in mice and humans mostly reside in peripheral tissues and induce anti-cancer cytotoxicity. Using single-cell-T-cell-receptor sequencing, we find that a cytotoxic γδT-cell subset of NMRs harbors a dominant clonotype, and that their conventional CD8 αβT cells exhibit modest clonotypic diversity. Consistently, perinatal NMR thymuses are considerably smaller than those of mice yet follow similar involution progression. Our findings suggest that NMRs have evolved under a relaxed intracellular pathogenic selective pressure that may have allowed cancer resistance and longevity to become stronger targets of selection to which the immune system has responded by utilizing γδT cells.

Immune microniches shape intestinal Treg function

The intestinal immune system is highly adapted to maintaining tolerance to the commensal microbiota and self-antigens while defending against invading pathogens1,2. Recognizing how the diverse network of local cells establish homeostasis and maintains it in the complex immune environment of the gut is critical to understanding how tolerance can be re-established following dysfunction, such as in inflammatory disorders. Although cell and molecular interactions that control T regulatory (Treg) cell development and function have been identified3,4, less is known about the cellular neighbourhoods and spatial compartmentalization that shapes microorganism-reactive Treg cell function. Here we used in vivo live imaging, photo-activation-guided single-cell RNA sequencing5-7 and spatial transcriptomics to follow the natural history of T cells that are reactive towards Helicobacter hepaticus through space and time in the settings of tolerance and inflammation. Although antigen stimulation can occur anywhere in the tissue, the lamina propria-but not embedded lymphoid aggregates-is the key microniche that supports effector Treg (eTreg) cell function. eTreg cells are stable once their niche is established; however, unleashing inflammation breaks down compartmentalization, leading to dominance of CD103+SIRPα+ dendritic cells in the lamina propria. We identify and validate the putative tolerogenic interaction between CD206+ macrophages and eTreg cells in the lamina propria and identify receptor-ligand pairs that are likely to govern the interaction. Our results reveal a spatial mechanism of tolerance in the lamina propria and demonstrate how knowledge of local interactions may contribute to the next generation of tolerance-inducing therapies.

Expansion of memory Vδ2 T cells following SARS-CoV-2 vaccination revealed by temporal single-cell transcriptomics

γδ T cells provide rapid cellular immunity against pathogens. Here, we conducted matched single-cell RNA-sequencing and γδ-TCR-sequencing to delineate the molecular changes in γδ T cells during a longitudinal study following mRNA SARS-CoV-2 vaccination. While the first dose of vaccine primes Vδ2 T cells, it is the second administration that significantly boosts their immune response. Specifically, the second vaccination uncovers memory features of Vδ2 T cells, shaped by the induction of AP-1 family transcription factors and characterized by a convergent central memory signature, clonal expansion, and an enhanced effector potential. This temporally distinct effector response of Vδ2 T cells was also confirmed in vitro upon stimulation with SARS-CoV-2 spike-peptides. Indeed, the second challenge triggers a significantly higher production of IFNγ by Vδ2 T cells. Collectively, our findings suggest that mRNA SARS-CoV-2 vaccination might benefit from the establishment of long-lasting central memory Vδ2 T cells to confer protection against SARS-CoV-2 infection.

Universal recording of immune cell interactions in vivo

Immune cells rely on transient physical interactions with other immune and non-immune populations to regulate their function1. To study these 'kiss-and-run' interactions directly in vivo, we previously developed LIPSTIC (labelling immune partnerships by SorTagging intercellular contacts)2, an approach that uses enzymatic transfer of a labelled substrate between the molecular partners CD40L and CD40 to label interacting cells. Reliance on this pathway limited the use of LIPSTIC to measuring interactions between CD4+ T helper cells and antigen-presenting cells, however. Here we report the development of a universal version of LIPSTIC (uLIPSTIC), which can record physical interactions both among immune cells and between immune and non-immune populations irrespective of the receptors and ligands involved. We show that uLIPSTIC can be used, among other things, to monitor the priming of CD8+ T cells by dendritic cells, reveal the steady-state cellular partners of regulatory T cells and identify germinal centre-resident T follicular helper cells on the basis of their ability to interact cognately with germinal centre B cells. By coupling uLIPSTIC with single-cell transcriptomics, we build a catalogue of the immune populations that physically interact with intestinal epithelial cells at the steady state and profile the evolution of the interactome of lymphocytic choriomeningitis virus-specific CD8+ T cells in multiple organs following systemic infection. Thus, uLIPSTIC provides a broadly useful technology for measuring and understanding cell-cell interactions across multiple biological systems.

Single-cell transcriptome landscape of circulating CD4+ T cell populations in autoimmune diseases

CD4+ T cells are key mediators of various autoimmune diseases; however, their role in disease progression remains unclear due to cellular heterogeneity. Here, we evaluated CD4+ T cell subpopulations using decomposition-based transcriptome characterization and canonical clustering strategies. This approach identified 12 independent gene programs governing whole CD4+ T cell heterogeneity, which can explain the ambiguity of canonical clustering. In addition, we performed a meta-analysis using public single-cell datasets of over 1.8 million peripheral CD4+ T cells from 953 individuals by projecting cells onto the reference and cataloging cell frequency and qualitative alterations of the populations in 20 diseases. The analyses revealed that the 12 transcriptional programs were useful in characterizing each autoimmune disease and predicting its clinical status. Moreover, genetic variants associated with autoimmune diseases showed disease-specific enrichment within the 12 gene programs. The results collectively provide a landscape of single-cell transcriptomes of CD4+ T cell subpopulations involved in autoimmune disease.

An atlas of cells in the human tonsil

Palatine tonsils are secondary lymphoid organs (SLOs) representing the first line of immunological defense against inhaled or ingested pathogens. We generated an atlas of the human tonsil composed of >556,000 cells profiled across five different data modalities, including single-cell transcriptome, epigenome, proteome, and immune repertoire sequencing, as well as spatial transcriptomics. This census identified 121 cell types and states, defined developmental trajectories, and enabled an understanding of the functional units of the tonsil. Exemplarily, we stratified myeloid slan-like subtypes, established a BCL6 enhancer as locally active in follicle-associated T and B cells, and identified SIX5 as putative transcriptional regulator of plasma cell maturation. Analyses of a validation cohort confirmed the presence, annotation, and markers of tonsillar cell types and provided evidence of age-related compositional shifts. We demonstrate the value of this resource by annotating cells from B cell-derived mantle cell lymphomas, linking transcriptional heterogeneity to normal B cell differentiation states of the human tonsil.

Adipose-tissue regulatory T cells are a consortium of subtypes that evolves with age and diet

Foxp3+CD4+ regulatory T (Treg) cells found within tissues regulate local immunity, inflammation, and homeostasis. Tregs in epididymal visceral adipose tissue (eVAT) are critical regulators of local and systemic inflammation and metabolism. During aging and under obesogenic conditions, eVAT Tregs undergo transcriptional and phenotypic changes and are important for containing inflammation and normalizing metabolic indices. We have employed single-cell RNA sequencing, single-cell Tra and Trb sequencing, adoptive transfers, photoconvertible mice, cellular interaction analyses, and in vitro cultures to dissect the evolving heterogeneity of eVAT Tregs with aging and obesity. Distinct Treg subtypes with distinguishable gene expression profiles and functional roles were enriched at differing ages and with differing diets. Like those in lean mice, eVAT Tregs in obese mice were not primarily recruited from the circulation but instead underwent local expansion and had a distinct and diversified T cell receptor repertoire. The different eVAT-Treg subtypes were specialized in different functions; for example, the subtypes enriched in lean, but not obese, mice suppressed adipogenesis. The existence of functionally divergent eVAT-Treg subtypes in response to obesogenic conditions presents possibilities for precision therapeutics in the context of obesity.

Characterization of CD8+ T cells in immune-privileged organs of ZIKV-infected Ifnar1-/- mice

Zika virus (ZIKV) infection caused neurological complications and male infertility, leading to the accumulation of antigen-specific immune cells in immune-privileged organs (IPOs). Thus, it is important to understand the immunological responses to ZIKV in IPOs. We extensively investigated the ZIKV-specific T cell immunity in IPOs in Ifnar1-/- mice, based on an immunodominant epitope E294-302 tetramer. The distinct kinetics and functions of virus-specific CD8+ T cells infiltrated into different IPOs were characterized, with late elevation in the brain and spinal cord. Single epitope E294-302-specific T cells can account for 20-60% of the total CD8+ T cells in the brain, spinal cord, and testicle and persist for at least 90 days in the brain and spinal cord. The E294-302-specific TCRαβs within the IPOs are featured with the majority of clonotypes utilizing TRAV9N-3 paired with diverse TRBV chains, but with distinct αβ paired clonotypes in 7 and 30 days post-infection. Specific chemokine receptors, Ccr2 and Ccr5, were selectively expressed in the E294-302-specific CD8+ T cells within the brain and testicle, indicating an IPO-oriented migration of virus-specific CD8+ T cells after infection. Overall, this study adds to the understanding of virus-specific CD8+ T cell responses for controlling and clearing ZIKV infection in IPOs.IMPORTANCEThe immune-privileged organs (IPOs), such as the central nervous system and testicles, presented pathogenicity and inflammation after Zika virus (ZIKV) infection with infiltrated CD8+ T cells. Our data show that CD8+ T cells keep up with virus increases and decreases in immune-privileged organs. Furthermore, our study provides the first ex vivo comparative analyses of the composition and diversity related to TCRα/β clonotypes across anatomical sites and ZIKV infection phases. We show that the vast majority of TCRα/β clonotypes in tissues utilize TRAV9N-3 with conservation. Specific chemokine expression, including Ccr2 and Ccr5, was found to be selectively expressed in the E294-302-specific CD8+ T cells within the brain and testicle, indicating an IPO-oriented migration of the virus-specific CD8+ T cells after the infection. Our study adds insights into the anti-viral immunological characterization and chemotaxis mechanism of virus-specific CD8+ T cells after ZIKV infection in different IPOs.

Macro CD5L+ deteriorates CD8+T cells exhaustion and impairs combination of Gemcitabine-Oxaliplatin-Lenvatinib-anti-PD1 therapy in intrahepatic cholangiocarcinoma

Intratumoral immune status influences tumor therapeutic response, but it remains largely unclear how the status determines therapies for patients with intrahepatic cholangiocarcinoma. Here, we examine the single-cell transcriptional and TCR profiles of 18 tumor tissues pre- and post- therapy of gemcitabine plus oxaliplatin, in combination with lenvatinib and anti-PD1 antibody for intrahepatic cholangiocarcinoma. We find that high CD8 GZMB+ and CD8 proliferating proportions and a low Macro CD5L+ proportion predict good response to the therapy. In patients with a poor response, the CD8 GZMB+ and CD8 proliferating proportions are increased, but the CD8 GZMK+ proportion is decreased after the therapy. Transition of CD8 proliferating and CD8 GZMB+ to CD8 GZMK+ facilitates good response to the therapy, while Macro CD5L+-CD8 GZMB+ crosstalk impairs the response by increasing CTLA4 in CD8 GZMB+. Anti-CTLA4 antibody reverses resistance of the therapy in intrahepatic cholangiocarcinoma. Our data provide a resource for predicting response of the combination therapy and highlight the importance of CD8+T-cell status conversion and exhaustion induced by Macro CD5L+ in influencing the response, suggesting future avenues for cancer treatment optimization.

Dandelion uses the single-cell adaptive immune receptor repertoire to explore lymphocyte developmental origins

Assessment of single-cell gene expression (single-cell RNA sequencing) and adaptive immune receptor (AIR) sequencing (scVDJ-seq) has been invaluable in studying lymphocyte biology. Here we introduce Dandelion, a computational pipeline for scVDJ-seq analysis. It enables the application of standard V(D)J analysis workflows to single-cell datasets, delivering improved V(D)J contig annotation and the identification of nonproductive and partially spliced contigs. We devised a strategy to create an AIR feature space that can be used for both differential V(D)J usage analysis and pseudotime trajectory inference. The application of Dandelion improved the alignment of human thymic development trajectories of double-positive T cells to mature single-positive CD4/CD8 T cells, generating predictions of factors regulating lineage commitment. Dandelion analysis of other cell compartments provided insights into the origins of human B1 cells and ILC/NK cell development, illustrating the power of our approach. Dandelion is available at https://www.github.com/zktuong/dandelion .

Single-cell multi-omics analysis of COVID-19 patients with pre-existing autoimmune diseases shows aberrant immune responses to infection

In COVID-19, hyperinflammatory and dysregulated immune responses contribute to severity. Patients with pre-existing autoimmune conditions can therefore be at increased risk of severe COVID-19 and/or associated sequelae, yet SARS-CoV-2 infection in this group has been little studied. Here, we performed single-cell analysis of peripheral blood mononuclear cells from patients with three major autoimmune diseases (rheumatoid arthritis, psoriasis, or multiple sclerosis) during SARS-CoV-2 infection. We observed compositional differences between the autoimmune disease groups coupled with altered patterns of gene expression, transcription factor activity, and cell-cell communication that substantially shape the immune response under SARS-CoV-2 infection. While enrichment of HLA-DRlow CD14+ monocytes was observed in all three autoimmune disease groups, type-I interferon signaling as well as inflammatory T cell and monocyte responses varied widely between the three groups of patients. Our results reveal disturbed immune responses to SARS-CoV-2 in patients with pre-existing autoimmunity, highlighting important considerations for disease treatment and follow-up.

Systems immunology spanning tumors, lymph nodes, and periphery

The immune system defines a complex network of tissues and cell types that orchestrate responses across the body in a dynamic manner. The local and systemic interactions between immune and cancer cells contribute to disease progression. Lymphocytes are activated in lymph nodes, traffic through the periphery, and impact cancer progression through their interactions with tumor cells. As a result, therapeutic response and resistance are mediated across tissues, and a comprehensive understanding of lymphocyte dynamics requires a systems-level approach. In this review, we highlight experimental and computational methods that can leverage the study of leukocyte trafficking through an immunomics lens and reveal how adaptive immunity shapes cancer.

Sustained Liver HBsAg Loss and Clonal T- and B-Cell Expansion upon Therapeutic DNA Vaccination Require Low HBsAg Levels

BACKGROUND

Suppression of HBV DNA, inhibition of HBV surface (HBsAg) production and therapeutic vaccination to reverse HBV-specific T-cell exhaustion in chronic HBV patients are likely required to achieve a functional cure. In the AAV-HBV mouse model, therapeutic vaccination can be effective in clearing HBV when HBsAg levels are low. Using a single-cell approach, we investigated the liver immune environment with different levels of HBsAg and sustained HBsAg loss through treatment with a GalNAc-HBV-siRNA followed by therapeutic vaccination.

METHODS

AAV-HBV-transduced C57BL/6 mice were treated with GalNAc-HBV-siRNA to lower HBsAg levels and then vaccinated using a DNA vaccine. We used single-cell RNA and V(D)J sequencing to understand liver immune microenvironment changes.

RESULTS

GalNAc-HBV-siRNA, followed by therapeutic vaccination, achieved sustained HBsAg loss in all mice. This was accompanied by CD4 follicular helper T-cell induction, polyclonal activation of CD8 T cells and clonal expansion of plasma cells that were responsible for antibody production.

CONCLUSIONS

This study provides novel insights into liver immune changes at the single-cell level, highlighting the correlation between induced reduction of HBsAg levels and clonal expansion of CD4, CD8 T cells and plasma cells in the liver upon HBV siRNA and subsequent therapeutic vaccination.

Scaling up single-cell RNA-seq data analysis with CellBridge workflow

SUMMARY

Single-cell RNA sequencing (scRNA-seq) has revolutionized the study of gene expression at the individual cell level, unraveling unprecedented insights into cellular heterogeneity. However, the analysis of scRNA-seq data remains a challenging and time-consuming task, often demanding advanced computational expertise, rendering it impractical for high-volume environments and applications. We present CellBridge, an automated workflow designed to simplify the standard procedures entailed in scRNA-seq data analysis, eliminating the need for specialized computational expertise. CellBridge utilizes state-of-the-art computational methods, integrating a range of advanced functionalities, covering the entire process from raw unaligned sequencing reads to cell type annotation. Hence, CellBridge accelerates the pace of discovery by seamlessly enabling insights into vast volumes of scRNA-seq data, without compromising workflow control and reproducibility.

AVAILABILITY AND IMPLEMENTATION

The source code, detailed documentation, and materials required to reproduce the results are available on GitHub and archived in Zenodo. For the CellBridge pre-processing step (v1.0.0), access the GitHub repository at https://github.com/Sanofi-Public/PMCB-ToBridge and the Zenodo archive at https://zenodo.org/records/10246161. For the CellBridge processing step (v1.0.0), visit the GitHub repository at https://github.com/Sanofi-Public/PMCB-CellBridge and the Zenodo archive at https://zenodo.org/records/10246046.

Single cell dynamics of tumor specificity vs bystander activity in CD8+ T cells define the diverse immune landscapes in colorectal cancer

CD8+ T cell activation via immune checkpoint blockade (ICB) is successful in microsatellite instable (MSI) colorectal cancer (CRC) patients. By comparison, the success of immunotherapy against microsatellite stable (MSS) CRC is limited. Little is known about the most critical features of CRC CD8+ T cells that together determine the diverse immune landscapes and contrasting ICB responses. Hence, we pursued a deep single cell mapping of CRC CD8+ T cells on transcriptomic and T cell receptor (TCR) repertoire levels in a diverse patient cohort, with additional surface proteome validation. This revealed that CRC CD8+ T cell dynamics are underscored by complex interactions between interferon-γ signaling, tumor reactivity, TCR repertoire, (predicted) TCR antigen-specificities, and environmental cues like gut microbiome or colon tissue-specific 'self-like' features. MSI CRC CD8+ T cells showed tumor-specific activation reminiscent of canonical 'T cell hot' tumors, whereas the MSS CRC CD8+ T cells exhibited tumor unspecific or bystander-like features. This was accompanied by inflammation reminiscent of 'pseudo-T cell hot' tumors. Consequently, MSI and MSS CRC CD8+ T cells showed overlapping phenotypic features that differed dramatically in their TCR antigen-specificities. Given their high discriminating potential for CD8+ T cell features/specificities, we used the single cell tumor-reactive signaling modules in CD8+ T cells to build a bulk tumor transcriptome classification for CRC patients. This "Immune Subtype Classification" (ISC) successfully distinguished various tumoral immune landscapes that showed prognostic value and predicted immunotherapy responses in CRC patients. Thus, we deliver a unique map of CRC CD8+ T cells that drives a novel tumor immune landscape classification, with relevance for immunotherapy decision-making.

Progenitor-like exhausted SPRY1+CD8+ T cells potentiate responsiveness to neoadjuvant PD-1 blockade in esophageal squamous cell carcinoma

Neoadjuvant immune checkpoint blockade (ICB) demonstrates promise in operable esophageal squamous cell carcinoma (ESCC), but lacks available efficacy biomarkers. Here, we perform single-cell RNA-sequencing of tumors from patients with ESCC undergoing neoadjuvant ICB, revealing a subset of exhausted CD8+ T cells expressing SPRY1 (CD8+ Tex-SPRY1) that displays a progenitor exhausted T cell (Tpex) phenotype and correlates with complete response to ICB. We validate CD8+ Tex-SPRY1 cells as an ICB-specific predictor of improved response and survival using independent ICB-/non-ICB cohorts and demonstrate that expression of SPRY1 in CD8+ T cells enforces Tpex phenotype and enhances ICB efficacy. Additionally, CD8+ Tex-SPRY1 cells contribute to proinflammatory phenotype of macrophages and functional state of B cells, which thereby promotes antitumor immunity by enhancing CD8+ T cell effector functions. Overall, our findings unravel progenitor-like CD8+ Tex-SPRY1 cells' role in effective responses to ICB for ESCC and inform mechanistic biomarkers for future individualized immunotherapy.

Blood transcriptomics analysis offers insights into variant-specific immune response to SARS-CoV-2

Bulk RNA sequencing (RNA-seq) of blood is typically used for gene expression analysis in biomedical research but is still rarely used in clinical practice. In this study, we argue that RNA-seq should be considered a routine diagnostic tool, as it offers not only insights into aberrant gene expression and splicing but also delivers additional readouts on immune cell type composition as well as B-cell and T-cell receptor (BCR/TCR) repertoires. We demonstrate that RNA-seq offers vital insights into a patient's immune status via integrative analysis of RNA-seq data from patients infected with various SARS-CoV-2 variants (in total 240 samples with up to 200 million reads sequencing depth). We compare the results of computational cell-type deconvolution methods (e.g., MCP-counter, xCell, EPIC, quanTIseq) to complete blood count data, the current gold standard in clinical practice. We observe varying levels of lymphocyte depletion and significant differences in neutrophil levels between SARS-CoV-2 variants. Additionally, we identify B and T cell receptor (BCR/TCR) sequences using the tools MiXCR and TRUST4 to show that - combined with sequence alignments and pBLAST - they could be used to classify a patient's disease. Finally, we investigated the sequencing depth required for such analyses and concluded that 10 million reads per sample is sufficient. In conclusion, our study reveals that computational cell-type deconvolution and BCR/TCR methods using bulk RNA-seq analyses can supplement missing CBC data and offer insights into immune responses, disease severity, and pathogen-specific immunity, all achievable with a sequencing depth of 10 million reads per sample.

Immune Surveillance of Acute Myeloid Leukemia Is Mediated by HLA-Presented Antigens on Leukemia Progenitor Cells

Therapy-resistant leukemia stem and progenitor cells (LSC) are a main cause of acute myeloid leukemia (AML) relapse. LSC-targeting therapies may thus improve outcome of patients with AML. Here we demonstrate that LSCs present HLA-restricted antigens that induce T-cell responses allowing for immune surveillance of AML. Using a mass spectrometry-based immunopeptidomics approach, we characterized the antigenic landscape of patient LSCs and identified AML- and AML/LSC-associated HLA-presented antigens absent from normal tissues comprising nonmutated peptides, cryptic neoepitopes, and neoepitopes of common AML driver mutations of NPM1 and IDH2. Functional relevance of shared AML/LSC antigens is illustrated by presence of their cognizant memory T cells in patients. Antigen-specific T-cell recognition and HLA class II immunopeptidome diversity correlated with clinical outcome. Together, these antigens shared among AML and LSCs represent prime targets for T cell-based therapies with potential of eliminating residual LSCs in patients with AML.

SIGNIFICANCE

The elimination of therapy-resistant leukemia stem and progenitor cells (LSC) remains a major challenge in the treatment of AML. This study identifies and functionally validates LSC-associated HLA class I and HLA class II-presented antigens, paving the way to the development of LSC-directed T cell-based immunotherapeutic approaches for patients with AML. See related commentary by Ritz, p. 430 . This article is featured in Selected Articles from This Issue, p. 419.

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.

Barcoding intracellular reverse transcription enables high-throughput phenotype-coupled T cell receptor analyses

Assays linking cellular phenotypes with T cell or B cell antigen receptor sequences are crucial for characterizing adaptive immune responses. Existing methodologies are limited by low sample throughput and high cost. Here, we present INtraCEllular Reverse Transcription with Sorting and sequencing (INCERTS), an approach that combines molecular indexing of receptor repertoires within intact cells and fluorescence-activated cell sorting (FACS). We demonstrate that INCERTS enables efficient processing of millions of cells from pooled human peripheral blood mononuclear cell (PBMC) samples while retaining robust association between T cell receptor (TCR) sequences and cellular phenotypes. We used INCERTS to discover antigen-specific TCRs from patients with cancer immunized with a novel mutant KRAS peptide vaccine. After ex vivo stimulation, 28 uniquely barcoded samples were pooled prior to FACS into peptide-reactive and non-reactive CD4+ and CD8+ populations. Combining complementary patient-matched single-cell RNA sequencing (scRNA-seq) data enabled retrieval of full-length, paired TCR alpha and beta chain sequences for future validation of therapeutic utility.

Combinatorial Immunotherapy with Agonistic CD40 Activates Dendritic Cells to Express IL12 and Overcomes PD-1 Resistance

Checkpoint inhibitors have revolutionized cancer treatment, but resistance remains a significant clinical challenge. Myeloid cells within the tumor microenvironment can modulate checkpoint resistance by either supporting or suppressing adaptive immune responses. Using an anti-PD-1-resistant mouse melanoma model, we show that targeting the myeloid compartment via CD40 activation and CSF1R blockade in combination with anti-PD-1 results in complete tumor regression in a majority of mice. This triple therapy combination was primarily CD40 agonist-driven in the first 24 hours after therapy and showed a similar systemic cytokine profile in human patients as was seen in mice. Functional single-cell cytokine secretion profiling of dendritic cells (DC) using a novel microwell assay identified a CCL22+CCL5+ IL12-secreting DC subset as important early-stage effectors of triple therapy. CD4+ and CD8+ T cells are both critical effectors of treatment, and systems analysis of single-cell RNA sequencing data supported a role for DC-secreted IL12 in priming T-cell activation and recruitment. Finally, we showed that treatment with a novel IL12 mRNA therapeutic alone was sufficient to overcome PD-1 resistance and cause tumor regression. Overall, we conclude that combining myeloid-based innate immune activation and enhancement of adaptive immunity is a viable strategy to overcome anti-PD-1 resistance.

Selection of cross-reactive T cells by commensal and food-derived yeasts drives cytotoxic TH1 cell responses in Crohn's disease

Aberrant CD4+ T cell reactivity against intestinal microorganisms is considered to drive mucosal inflammation in inflammatory bowel diseases. The disease-relevant microbial species and the corresponding microorganism-specific, pathogenic T cell phenotypes remain largely unknown. In the present study, we identified common gut commensal and food-derived yeasts, as direct activators of altered CD4+ T cell reactions in patients with Crohn's disease (CD). Yeast-responsive CD4+ T cells in CD display a cytotoxic T helper cell (TH1 cell) phenotype and show selective expansion of T cell clones that are highly cross-reactive to several commensal, as well as food-derived, fungal species. This indicates cross-reactive T cell selection by repeated encounter with conserved fungal antigens in the context of chronic intestinal disease. Our results highlighted a role of yeasts as drivers of aberrant CD4+ T cell reactivity in patients with CD and suggest that both gut-resident fungal commensals and daily dietary intake of yeasts might contribute to chronic activation of inflammatory CD4+ T cell responses in patients with CD.

Immune profiling after allogeneic hematopoietic cell transplantation in pediatric acute myeloid leukemia

Although allogeneic hematopoietic cell transplant (allo-HCT) is curative for high-risk pediatric acute myeloid leukemia (AML), disease relapse remains the primary cause of posttransplant mortality. To identify pressures imposed by allo-HCT on AML cells that escape the graft-versus-leukemia effect, we evaluated immune signatures at diagnosis and posttransplant relapse in bone marrow samples from 4 pediatric patients using a multimodal single-cell proteogenomic approach. Downregulation of major histocompatibility complex class II expression was most profound in progenitor-like blasts and accompanied by correlative changes in transcriptional regulation. Dysfunction of activated natural killer cells and CD8+ T-cell subsets at relapse was evidenced by the loss of response to interferon gamma, tumor necrosis factor α signaling via NF-κB, and interleukin-2/STAT5 signaling. Clonotype analysis of posttransplant relapse samples revealed an expansion of dysfunctional T cells and enrichment of T-regulatory and T-helper cells. Using novel computational methods, our results illustrate a diverse immune-related transcriptional signature in posttransplant relapses not previously reported in pediatric AML.

Intradermal skin test with mRNA vaccines as a surrogate marker of T cell immunity in immunocompromised patients

OBJECTIVES

Intradermal skin test (IDT) with mRNA vaccines may represent a simple, reliable, and affordable tool to measure T cell response in immunocompromised patients who failed to mount serological responses following vaccination with mRNA covid-19 vaccines.

METHODS

We compared anti-SARS-CoV-2 antibodies and cellular responses in vaccinated immunocompromised patients (n = 58), healthy seronegative naive controls (NC, n = 8), and healthy seropositive vaccinated controls (VC, n = 32) by Luminex, spike-induced IFN-γ Elispot and an IDT. A skin biopsy 24 h after IDT and single-cell RNAseq was performed in three vaccinated volunteers.

RESULTS

Twenty-five percent of seronegative NC had a positive Elispot (2/8) and IDT (1/4), compared to 95% (20/21) and 93% (28/30) in seropositive VC, respectively. Single-cell RNAseq data in the skin of VC showed a predominant mixed population of effector helper and cytotoxic T cells. The TCR repertoire revealed 18/1064 clonotypes with known specificities against SARS-CoV-2, among which six were spike-specific. Seronegative immunocompromised patients with positive Elispot and IDT were in 83% (5/6) treated with B cell-depleting reagents, while those with negative IDT were all transplant recipients.

CONCLUSIONS

Our results indicate that delayed local reaction to IDT reflects vaccine-induced T-cell immunity opening new perspectives to monitor seronegative patients and elderly populations with waning immunity.

Induction of bronchus-associated lymphoid tissue is an early life adaptation for promoting human B cell immunity

Infants and young children are more susceptible to common respiratory pathogens than adults but can fare better against novel pathogens like severe acute respiratory syndrome coronavirus 2. The mechanisms by which infants and young children mount effective immune responses to respiratory pathogens are unknown. Through investigation of lungs and lung-associated lymph nodes from infant and pediatric organ donors aged 0-13 years, we show that bronchus-associated lymphoid tissue (BALT), containing B cell follicles, CD4+ T cells and functionally active germinal centers, develop during infancy. BALT structures are prevalent around lung airways during the first 3 years of life, and their numbers decline through childhood coincident with the accumulation of memory T cells. Single-cell profiling and repertoire analysis reveals that early life lung B cells undergo differentiation, somatic hypermutation and immunoglobulin class switching and exhibit a more activated profile than lymph node B cells. Moreover, B cells in the lung and lung-associated lymph nodes generate biased antibody responses to multiple respiratory pathogens compared to circulating antibodies, which are mostly specific for vaccine antigens in the early years of life. Together, our findings provide evidence for BALT as an early life adaptation for mobilizing localized immune protection to the diverse respiratory challenges during this formative life stage.

Best practices for single-cell analysis across modalities

Recent advances in single-cell technologies have enabled high-throughput molecular profiling of cells across modalities and locations. Single-cell transcriptomics data can now be complemented by chromatin accessibility, surface protein expression, adaptive immune receptor repertoire profiling and spatial information. The increasing availability of single-cell data across modalities has motivated the development of novel computational methods to help analysts derive biological insights. As the field grows, it becomes increasingly difficult to navigate the vast landscape of tools and analysis steps. Here, we summarize independent benchmarking studies of unimodal and multimodal single-cell analysis across modalities to suggest comprehensive best-practice workflows for the most common analysis steps. Where independent benchmarks are not available, we review and contrast popular methods. Our article serves as an entry point for novices in the field of single-cell (multi-)omic analysis and guides advanced users to the most recent best practices.

Multimodal immune cell phenotyping in GI biopsies reveals microbiome-related T cell modulations in human GvHD

Acute graft-versus-host disease (aGvHD) is a significant complication after allogeneic hematopoietic stem cell transplantation (aHSCT), but major factors determining disease severity are not well defined yet. By combining multiplexed tissue imaging and single-cell RNA sequencing on gastrointestinal biopsies from aHSCT-treated individuals with fecal microbiome analysis, we link high microbiome diversity and the abundance of short-chain fatty acid-producing bacteria to the sustenance of suppressive regulatory T cells (Tregs). Furthermore, aGvHD severity strongly associates with the clonal expansion of mainly CD8 T cells, which we find distributed over anatomically distant regions of the gut, persistent over time, and inversely correlated with the presence of suppressive Tregs. Overall, our study highlights the pathophysiological importance of expanded CD8 T cell clones in the progression of aGvHD toward more severe clinical manifestations and strongly supports the further development of microbiome interventions as GvHD treatment via repopulation of the gut Treg niche to suppress inflammation.

Clonally expanded PD-1-expressing T cells are enriched in synovial fluid of juvenile idiopathic arthritis patients

Juvenile idiopathic arthritis (JIA) is the most common chronic rheumatic condition in childhood. The disease etiology remains largely unknown; however, a key role in JIA pathogenesis is surely mediated by T cells. T-lymphocytes activity is controlled via signals, known as immune checkpoints. Delivering an inhibitory signal or blocking a stimulatory signal to achieve immune suppression is critical in autoimmune diseases. However, the role of immune checkpoints in chronic inflammation and autoimmunity must still be deciphered. In this study, we investigated at the single-cell level the feature of T cells in JIA chronic inflammation, both at the transcriptome level via single-cell RNA sequencing and at the protein level by flow cytometry. We found that despite the heterogeneity in the composition of synovial CD4+ and CD8+ T cells, those characterized by PD-1 expression were clonally expanded tissue-resident memory (Trm)-like cells and displayed the highest proinflammatory capacity, suggesting their active contribution in sustaining chronic inflammation in situ. Our data support the concept that novel therapeutic strategies targeting PD-1 may be effective in the treatment of JIA. With this approach, it may become possible to target overactive T cells regardless of their cytokine production profile.

A current and future perspective on T cell receptor repertoire profiling

T cell receptors (TCR) play a vital role in the immune system's ability to recognize and respond to foreign antigens, relying on the highly polymorphic rearrangement of TCR genes. The recognition of autologous peptides by adaptive immunity may lead to the development and progression of autoimmune diseases. Understanding the specific TCR involved in this process can provide insights into the autoimmune process. RNA-seq (RNA sequencing) is a valuable tool for studying TCR repertoires by providing a comprehensive and quantitative analysis of the RNA transcripts. With the development of RNA technology, transcriptomic data must provide valuable information to model and predict TCR and antigen interaction and, more importantly, identify or predict neoantigens. This review provides an overview of the application and development of bulk RNA-seq and single-cell (SC) RNA-seq to examine the TCR repertoires. Furthermore, discussed here are bioinformatic tools that can be applied to study the structural biology of peptide/TCR/MHC (major histocompatibility complex) and predict antigenic epitopes using advanced artificial intelligence tools.

Single-cell landscape of primary central nervous system diffuse large B-cell lymphoma

Understanding tumor heterogeneity and immune infiltrates within the tumor-immune microenvironment (TIME) is essential for the innovation of immunotherapies. Here, combining single-cell transcriptomics and chromatin accessibility sequencing, we profile the intratumor heterogeneity of malignant cells and immune properties of the TIME in primary central nervous system diffuse large B-cell lymphoma (PCNS DLBCL) patients. We demonstrate diverse malignant programs related to tumor-promoting pathways, cell cycle and B-cell immune response. By integrating data from independent systemic DLBCL and follicular lymphoma cohorts, we reveal a prosurvival program with aberrantly elevated RNA splicing activity that is uniquely associated with PCNS DLBCL. Moreover, a plasmablast-like program that recurs across PCNS/activated B-cell DLBCL predicts a worse prognosis. In addition, clonally expanded CD8 T cells in PCNS DLBCL undergo a transition from a pre-exhaustion-like state to exhaustion, and exhibit higher exhaustion signature scores than systemic DLBCL. Thus, our study sheds light on potential reasons for the poor prognosis of PCNS DLBCL patients, which will facilitate the development of targeted therapy.

Personalized RNA neoantigen vaccines stimulate T cells in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is lethal in 88% of patients1, yet harbours mutation-derived T cell neoantigens that are suitable for vaccines 2,3. Here in a phase I trial of adjuvant autogene cevumeran, an individualized neoantigen vaccine based on uridine mRNA-lipoplex nanoparticles, we synthesized mRNA neoantigen vaccines in real time from surgically resected PDAC tumours. After surgery, we sequentially administered atezolizumab (an anti-PD-L1 immunotherapy), autogene cevumeran (a maximum of 20 neoantigens per patient) and a modified version of a four-drug chemotherapy regimen (mFOLFIRINOX, comprising folinic acid, fluorouracil, irinotecan and oxaliplatin). The end points included vaccine-induced neoantigen-specific T cells by high-threshold assays, 18-month recurrence-free survival and oncologic feasibility. We treated 16 patients with atezolizumab and autogene cevumeran, then 15 patients with mFOLFIRINOX. Autogene cevumeran was administered within 3 days of benchmarked times, was tolerable and induced de novo high-magnitude neoantigen-specific T cells in 8 out of 16 patients, with half targeting more than one vaccine neoantigen. Using a new mathematical strategy to track T cell clones (CloneTrack) and functional assays, we found that vaccine-expanded T cells comprised up to 10% of all blood T cells, re-expanded with a vaccine booster and included long-lived polyfunctional neoantigen-specific effector CD8+ T cells. At 18-month median follow-up, patients with vaccine-expanded T cells (responders) had a longer median recurrence-free survival (not reached) compared with patients without vaccine-expanded T cells (non-responders; 13.4 months, P = 0.003). Differences in the immune fitness of the patients did not confound this correlation, as responders and non-responders mounted equivalent immunity to a concurrent unrelated mRNA vaccine against SARS-CoV-2. Thus, adjuvant atezolizumab, autogene cevumeran and mFOLFIRINOX induces substantial T cell activity that may correlate with delayed PDAC recurrence.