VDJtools: Unifying Post-analysis of T Cell Receptor Repertoires

Immunopipe: a comprehensive and flexible scRNA-seq and scTCR-seq data analysis pipeline

Single-cell sequencing technologies provide us with information at the level of individual cells. Combining single-cell RNA-seq and single-cell TCR-seq profiling enables the exploration of cell heterogeneity and T-cell receptor repertoires simultaneously. Integrating both types of data can play a crucial role in enhancing our understanding of T-cell-mediated immunity and, in turn, facilitate the advancement of immunotherapy. Here, we present immunopipe, a comprehensive and flexible pipeline to perform integrated analysis of scRNA-seq and scTCR-seq data. In addition to the command line tool, we provide a user-friendly web interface for pipeline configuration and execution monitoring, benefiting researchers without extensive programming experience. With its comprehensive functionality and ease of use, immunopipe empowers researchers to uncover valuable insights from scRNA-seq and scTCR-seq data, ultimately advancing the understanding of immune responses and immunotherapy development.

Cholinergic regulation of thymocyte negative selection

The immune and nervous systems use a common chemical language for communication, namely, the cholinergic signaling involving acetylcholine (ACh) and its receptors (AChRs). Whether and how this language also regulates the development of the immune system is poorly understood. Here, we show that mouse CD4+CD8+ double-positive thymocytes express high levels of α9 nicotinic AChR (nAChR) and that this receptor controls thymic negative selection. α9 nAChR-deficient mice show an altered T cell receptor (TCR) repertoire and reduced CD4+ and CD8+ T cells in a mixed bone marrow chimera setting. α9 nAChR-mediated signaling regulates TCR strength and thymocyte survival. Thymic tuft cells, B cells and some T cells express choline acetyltransferase and are potential ACh sources, with ACh derived from T cells having the most important role. Furthermore, α9 nAChR deficiency during thymocyte development contributes to the altered development of autoimmune diseases in mice. Our results thus reveal a mechanism controlling immune cell development that involves cholinergic signaling.

Immunosequencing identifies signatures of T cell responses for early detection of nasopharyngeal carcinoma

To identify nasopharyngeal carcinoma (NPC)-relevant T cell receptors (TCRs), we profile the repertoires of peripheral blood TCRβ chains from 228 NPC patients, 241 at-risk controls positive for serum Epstein-Barr virus (EBV) VCA-IgA antibody, and 251 seronegative controls. We develop a TCR-based signature (T-score) based on 208 NPC-enriched CDR3β sequences, which accurately diagnoses NPC in both the original and independent validation cohorts. Notably, a higher T-score, associated with a shorter time interval to NPC diagnosis, effectively identifies early-stage NPC among EBV-seropositive at-risk individuals prior to clinical diagnosis. These NPC-enriched TCRs react against not only EBV-specific antigens but also non-EBV antigens expressed by NPC cells, indicating a broad range of specificities. Moreover, the abundance of NPC-enriched CD8+ T cells in blood correlates with the infiltration of non-exhausted T cell counterparts in tumors and predicts prolonged survival, suggesting that these NPC-enriched T cells have significant potential for disease monitoring and therapeutic applications.

Multi-omics sequencing of gastroesophageal junction adenocarcinoma reveals prognosis-relevant key factors and a novel immunogenomic classification

BACKGROUND

Gastroesophageal junction adenocarcinoma (GEJAC) exhibits distinct molecular characteristics due to its unique anatomical location. We sought to investigate effective and reliable molecular classification of GEJAC to guide personalized treatment.

METHODS

We analyzed the whole genomic, transcriptomic, T-cell receptor repertoires, and immunohistochemical data in 92 GEJAC patients and delineated the landscape of genetic and immune alterations. In addition to COSMIC nomenclature, the de novo nomenclature was also utilized to define signatures and investigate their correlation with survival. A novel molecular subtype was developed and validated in other cohorts.

RESULTS

We found 30 mutated driver genes, 7 novel genomic signatures, 3 copy-number variations, and 2 V-J gene usages related to prognosis that were not identified in previous study. A high frequency of COSMIC-SBS-384-1 and De novo-SV-32-A was associated with more neoantigen generation and a better survival. Using 19 molecular features, we identified three immune-related subtypes (immune inflamed, intermediate, and deserted) with discrete profiles of genomic signatures, immune status, and clinical outcome. The immune deserted subtype (27.2%) was characterized by an earlier KRAS mutation, worse immune reaction, and prognosis than the other two subtypes. The immune inflamed subtypes exhibited the highest levels of neoantigens, TCR/pMHC-binding strength, CD8 + T-cell infiltration, IFN-α/γ response pathways, and survival rate.

CONCLUSIONS

These results emphasize the immune reaction and prognostic value of novel molecular classifications based on multi-omics data and provide a solid basis for better management of GEJAC.

Two distinct durable human class-switched memory B cell populations are induced by vaccination and infection

Memory lymphocytes are durable cells that persist in the absence of antigen, but few human B cell subsets have been characterized in terms of durability. The relative durability of eight non-overlapping human B cell sub-populations covering 100% of all human class-switched B cells was interrogated. Only two long-lived B cell populations persisted in the relative absence of antigen. In addition to canonical germinal center-derived switched-memory B cells with an IgD-CD27+CXCR5+ phenotype, a second, non-canonical, but distinct memory population of IgD-CD27-CXCR5+ DN1 B cells was also durable, exhibited a unique TP63-linked transcriptional and anti-apoptotic signature, had low levels of somatic hypermutation, but was more clonally expanded than canonical switched-memory B cells. DN1 B cells likely evolved to preserve immunological breadth and may represent the human counterparts of rodent extrafollicular memory B cells that, unlike canonical memory B cells, can enter germinal centers and facilitate B cell and antibody evolution.

Efficient boosting of Omicron-reactive memory B cells after breakthrough infection protects from repeated exposure

Exploring the impact of persistent mutations in SARS-CoV-2 variants and reduced immunity on breakthrough infections (BTIs) is crucial, particularly in understanding how antigen-specific memory B cells (MBCs) respond to new variants. We followed 107 participants who received the ancestral inactivated vaccine and experienced one or two Omicron BTIs over six months. Using flow cytometry, SARS-CoV-2 antigen probes, single-cell RNA sequencing, and B cell receptor (BCR) profiling, we assessed MBCs and immune diversity. Our findings revealed that although neutralizing antibody levels decreased over time, the number of specific MBCs remained stable and matured progressively. Notably, pre-existing Omicron-specific MBCs played a key role in preventing secondary Omicron infections. Differential gene analysis showed enrichment in antigen processing and immune regulation pathways, while clonal lineage analysis revealed more B cell expansion and V(D)J gene-specific rearrangements in high neutralization samples. These results emphasize MBCs' critical role in long-term immunity and inform future vaccination strategies.

Repertoire-based mapping and time-tracking of T helper cell subsets in scRNA-Seq

Introduction

The functional programs of CD4+ T helper (Th) cell clones play a central role in shaping immune responses to different challenges. While advances in single-cell RNA sequencing (scRNA-Seq) have significantly improved our understanding of the diversity of Th cells, the relationship between scRNA-Seq clusters and the traditionally characterized Th subsets remains ambiguous.

Methods

In this study, we introduce TCR-Track, a method leveraging immune repertoire data to map phenotypically sorted Th subsets onto scRNA-Seq profiles.

Results and discussion

This approach accurately positions the Th1, Th1-17, Th17, Th22, Th2a, Th2, T follicular helper (Tfh), and regulatory T-cell (Treg) subsets, outperforming mapping based on CITE-Seq. Remarkably, the mapping is tightly focused on specific scRNA-Seq clusters, despite 4-year interval between subset sorting and the effector CD4+ scRNA-Seq experiment. These findings highlight the intrinsic program stability of Th clones circulating in peripheral blood. Repertoire overlap analysis at the scRNA-Seq level confirms that the circulating Th1, Th2, Th2a, Th17, Th22, and Treg subsets are clonally independent. However, a significant clonal overlap between the Th1 and cytotoxic CD4+ T-cell clusters suggests that cytotoxic CD4+ T cells differentiate from Th1 clones. In addition, this study resolves a longstanding ambiguity: we demonstrate that, while CCR10+ Th cells align with a specific Th22 scRNA-Seq cluster, CCR10-CCR6+CXCR3-CCR4+ cells, typically classified as Th17, represent a mixture of bona fide Th17 cells and clonally unrelated CCR10low Th22 cells. The clear distinction between the Th17 and Th22 subsets should influence the development of vaccine- and T-cell-based therapies. Furthermore, we show that severe acute SARS-CoV-2 infection induces systemic type 1 interferon (IFN) activation of naive Th cells. An increased proportion of effector IFN-induced Th cells is associated with a moderate course of the disease but remains low in critical COVID-19 cases. Using integrated scRNA-Seq, TCR-Track, and CITE-Seq data from 122 donors, we provide a comprehensive Th scRNA-Seq reference that should facilitate further investigation of Th subsets in fundamental and clinical studies.

A comprehensive review and evaluation of species richness estimation

MOTIVATION

The statistical problem of estimating the total number of distinct species in a population (or distinct elements in a multiset), given only a small sample, occurs in various areas, ranging from the unseen species problem in ecology to estimating the diversity of immune repertoires. Accurately estimating the true richness from very small samples is challenging, in particular for highly diverse populations with many rare species. Depending on the application, different estimation strategies have been proposed that incorporate explicit or implicit assumptions about either the species distribution or about the sampling process. These methods are scattered across the literature, and an extensive overview of their assumptions, methodology, and performance is currently lacking.

RESULTS

We comprehensively review and evaluate a variety of existing methods on real and simulated data with different compositions of rare and abundant species. Our evaluation shows that, depending on species composition, different methods provide the most accurate richness estimates. Simple methods based on the observed number of singletons yield accurate asymptotic lower bounds for several of the tested simulated species compositions, but tend to underestimate the true richness for heterogeneous populations and small samples containing 1% to 5% of the population. When the population size is known, upsampling (extrapolating) estimators such as PreSeq and RichnEst yield accurate estimates of the total species richness in a sample that is up to 10 times larger than the observed sample.

AVAILABILITY

Source code for data simulation and richness estimation is available at https://gitlab.com/rahmannlab/speciesrichness.

Single-cell analysis of aplastic anemia reveals a convergence of NK and NK-like CD8+ T cells with a disease-associated TCR signature

Immune aplastic anemia (AA) is a life-threatening bone marrow failure disorder driven by an autoimmune T cell attack against hematopoietic stem and progenitor cells (HSPCs). However, the exact autoantigen targets and role of other immune cells in the pathogenesis of AA are unknown. Here, we analyzed a cohort of 218 patients with AA using single-cell RNA and T cell receptor (TCR) αβ sequencing, TCRβ sequencing, flow cytometry, and plasma cytokine profiling. We identified natural killer (NK) cells and CD8+ terminally differentiated effector T (TEMRA) cells expressing NK receptors with AA-associated TCRβ motifs as the most dysregulated immune cell populations in AA bone marrow. Functional coculture experiments using primary HSPCs and immune cells showed that NK cells cannot kill HSPCs alone but may sensitize HSPCs to CD8+ T cell-mediated killing through production of interferons. Furthermore, HSPCs induced activation of T cell clones with CD8+ TEMRA NK-like phenotype in coculture. Our results reveal a convergent phenotype of innate and adaptive immune cells that may drive AA.

Comprehensive Analysis of TCR and BCR Repertoires: Insights into Methodologies, Challenges, and Applications

The diversity of T-cell receptors (TCRs) and B-cell receptors (BCRs) underpins the adaptive immune system's ability to recognize and respond to a wide array of antigens. Recent advancements in RNA sequencing have expanded its application beyond transcriptomics to include the analysis of immune repertoires, enabling the exploration of TCR and BCR sequences across various physiological and pathological contexts. This review highlights key methodologies and considerations for TCR and BCR repertoire analysis, focusing on the technical aspects of receptor sequence extraction, data processing, and clonotype identification. We compare the use of bulk and single-cell sequencing, discuss computational tools and pipelines, and evaluate the implications of examining specific receptor regions such as CDR3. By integrating immunology, bioinformatics, and clinical research, immune repertoire analysis provides valuable insights into immune function, therapeutic responses, and precision medicine approaches, advancing our understanding of health and disease.

IL-4 induces CD22 expression to restrain the effector program of virtual memory T cells

Parasitic helminths induce the production of interleukin-4 (IL-4), which causes the expansion of virtual memory CD8+ T cells (TVM cells), a cell subset that contributes to the control of coinfection with intracellular pathogens. However, the mechanisms regulating IL-4-dependent TVM cell activation and expansion remain ill defined. Here, we used single-cell RNA sequencing of CD8+ T cells to identify pathways that control IL-4-dependent TVM cell responses. Gene signature analysis of CD8+ T cells identified a cell cluster marked by CD22, a canonical regulator of B cell activation, as a selective surface marker of IL-4-induced TVM cells. CD22+ TVM cells were enriched for interferon-γ and granzyme A and retained a diverse TCR repertoire while enriched in self-reactive CDR3 sequences. CD22 intrinsically regulated the IL-4-induced CD8+ T cell effector program, resulting in reduced responsiveness of CD22+ TVM cells and regulatory functions to infection and inflammation. Thus, helminth-induced IL-4 drives the expansion and activation of TVM cells that is counterinhibited by CD22.

Alzheimer's disease may develop from changes in the immune system, cell cycle, and protein processing following alterations in ribosome function

The prevalence of Alzheimer's disease (AD) is increasing as society ages. The details of AD pathogenesis have not been fully elucidated, and a comprehensive gene expression analysis of the process leading up to the onset of AD would be helpful for understanding the mechanism. We performed an RNA sequencing analysis on a cohort of 1227 Japanese blood samples, representing 424 AD patients, 543 individuals with mild cognitive impairment (MCI), and 260 cognitively normal (CN) individuals. A total of 883 and 1169 statistically significant differentially expressed genes (DEGs) were identified between CN and MCI (CN-MCI) and between MCI and AD (MCI-AD), respectively. Pathway analyses using these DEGs, followed by protein-protein interaction network analysis, revealed key roles of ribosomal function in MCI progression, whereas immune responses, cell cycle, and protein processing in endoplasmic reticulum were involved in AD progression. Our findings indicate that the onset of AD might be associated with gene expression changes in the immune system, cell cycle, and protein processing following alterations in the expression of ribosomal protein genes during the MCI stage, although validation using brain tissue samples will be necessary in the future. Given the known effectiveness of delaying MCI progression in preventing AD, the genes related to ribosomal function might emerge as biomarkers for early diagnosis.

The compositional behavior of the human T cell receptor repertoire in ovarian cancer compared to healthy donors

The distinctive characteristics of an individual's T cell receptor repertoire are crucial in recognizing and responding to a diverse array of antigens, contributing to immune specificity and adaptability. The repertoire, famously vast due to a series of cellular mechanisms, can be quantified using repertoire sequencing. In this study, we sampled the repertoire of 85 women: ovarian cancer patients (OC) and healthy donors (HD), generating a dataset of T cell clones and their abundance. For the alpha chain we obtained 6.4·106 reads, with an average of 75936 clones per sample, and an average of 30607 clonotypes per sample. For the beta chain we obtained 13.6·106 reads, with an average of 160400 clones per sample, and an average of 70071 clonotypes per sample. The changes in dynamics of the repertoire can be observed in response to disease, with specific clones undergoing clonal expansion and contraction. The data provided here offers a unique view of immune system behavior in health and disease and can be used to stratify OC and HD.

Identification and validation of a T cell receptor targeting KRAS G12V in HLA-A*11:01 pancreatic cancer patients

T cells targeting a KRAS mutation can induce durable tumor regression in some patients with metastatic epithelial cancer. It is unknown whether T cells targeting mutant KRAS that are capable of killing tumor cells can be identified from peripheral blood of patients with pancreatic cancer. We developed an in vitro stimulation approach and identified HLA-A*11:01-restricted KRAS G12V-reactive CD8+ T cells and HLA-DRB1*15:01-restricted KRAS G12V-reactive CD4+ T cells from peripheral blood of 2 out of 6 HLA-A*11:01-positive patients with pancreatic cancer whose tumors expressed KRAS G12V. The HLA-A*11:01-restricted KRAS G12V-reactive T cell receptor (TCR) was isolated and validated to specifically recognize the KRAS G12V8-16 neoepitope. While T cells engineered to express this TCR specifically recognized all 5 tested human HLA-A*11:01+ and KRAS G12V+ pancreatic cancer organoids, the recognition was often modest, and tumor cell killing was observed in only 2 out of 5 organoids. IFN-γ priming of the organoids enhanced the recognition and killing by the TCR-engineered T cells. The TCR-engineered T cells could significantly slow the growth of an established organoid-derived xenograft in immunodeficient mice. Our data suggest that this TCR has potential for use in TCR-gene therapy, but additional strategies that enhance tumor recognition by the TCR-engineered T cells likely will be required to increase clinical activity.

Systematic evaluation of intratumoral and peripheral BCR repertoires in three cancers

The current understanding of humoral immune response in cancer patients suggests that tumors may be infiltrated with diffuse B cells of extra-tumoral origin or may develop organized lymphoid structures, where somatic hypermutation and antigen-driven selection occur locally. These processes are believed to be significantly influenced by the tumor microenvironment through secretory factors and biased cell-cell interactions. To explore the manifestation of this influence, we used deep unbiased immunoglobulin profiling and systematically characterized the relationships between B cells in circulation, draining lymph nodes (draining LNs), and tumors in 14 patients with three human cancers. We demonstrated that draining LNs are differentially involved in the interaction with the tumor site, and that significant heterogeneity exists even between different parts of a single lymph node (LN). Next, we confirmed and elaborated upon previous observations regarding intratumoral immunoglobulin heterogeneity. We identified B cell receptor (BCR) clonotypes that were expanded in tumors relative to draining LNs and blood and observed that these tumor-expanded clonotypes were less hypermutated than non-expanded (ubiquitous) clonotypes. Furthermore, we observed a shift in the properties of complementarity-determining region 3 of the BCR heavy chain (CDR-H3) towards less mature and less specific BCR repertoire in tumor-infiltrating B-cells compared to circulating B-cells, which may indicate less stringent control for antibody-producing B cell development in tumor microenvironment (TME). In addition, we found repertoire-level evidence that B-cells may be selected according to their CDR-H3 physicochemical properties before they activate somatic hypermutation (SHM). Altogether, our work outlines a broad picture of the differences in the tumor BCR repertoire relative to non-tumor tissues and points to the unexpected features of the SHM process.

A subunit vaccine Ag85A-LpqH focusing on humoral immunity provides substantial protection against tuberculosis in mice

The importance of humoral immunity in combating TB has gained extensive recognition. In this study, a subunit vaccine named Ag85A-LpqH (AL) was prepared by fusing the antigen Ag85A proved to induce robust T cell immune responses, and LpqH was shown to produce protective antibodies. The prevention and BCG prime-boost mouse models were established to test the vaccine efficacy. The results indicate that Ag85A-LpqH can induce substantial protection by reducing bacterial loads and pathological lesions. This vaccine can induce robust antibody responses, as well as T cell immune responses especially strong CD8+ T cell responses. Moreover, the serum from AL-immunized mice can reduce the bacterial load and lung pathology in mice. B cell receptor (BCR) sequencing revealed a notable rise in BCR diversity among mice immunized with AL. These results indicate that Ag85A-LpqH can be a promising vaccine candidate for tuberculosis prevention and control.

MHC-related protein 1-restricted recognition of cancer via a semi-invariant TCR-α chain

The T cell antigen presentation platform MR1 consists of 6 allomorphs in humans that differ by no more than 5 amino acids. The principal function of this highly conserved molecule involves presenting microbial metabolites to the abundant mucosal-associated invariant T (MAIT) cell subset. Recent developments suggest that the role of MR1 extends to presenting antigens from cancer cells, a function dependent on the K43 residue in the MR1 antigen binding cleft. Here, we successfully cultured cancer-activated, MR1-restricted T cells from multiple donors and confirmed that they recognized a wide range of cancer types expressing the most common MR1*01 and/or MR1*02 allomorphs (over 95% of the population), while remaining inert to healthy cells including healthy B cells and monocytes. Curiously, in all but one donor these T cells were found to incorporate a conserved TCR-α chain motif, CAXYGGSQGNLIF (where X represents 3-5 amino acids), because of pairing between 10 different TRAV genes and the TRAJ42 gene segment. This semi-invariance in the TCR-α chain is reminiscent of MAIT cells and suggests recognition of a conserved antigen bound to K43.

ATLAS-seq: a microfluidic single-cell TCR screen for antigen-reactive TCRs

Discovering antigen-reactive T cell receptors (TCRs) is central to developing effective engineered T cell immunotherapies. However, the conventional technologies for isolating antigen-reactive TCRs (i.e., major histocompatibility complex (MHC) multimer staining) focus on high-affinity interactions between the TCR and MHC-antigen complex, and may fail to identify TCRs with high efficacy for activating T cells. Here, we develop a microfluidic single-cell screening method for antigen-reactive T cells named ATLAS-seq (Aptamer-based T Lymphocyte Activity Screening and SEQuencing). This technology isolates and characterizes activated T cells via an aptamer-based fluorescent molecular sensor, which monitors the cytotoxic cytokine IFNγ secretion from single T cells upon antigen stimulation, followed by single-cell RNA and single-cell TCR sequencing. We use ATLAS-seq to screen TCRs reactive to cytomegalovirus (CMV) or prostate specific antigen (PSA) from peripheral blood mononuclear cells (PBMCs). ATLAS-seq identifies distinct TCR clonotype populations with higher T cell activation levels compared to TCRs recovered by MHC multimer staining. Select TCR clonotypes from ATLAS-seq are more efficient in target cell killing than those from MHC multimer staining. Collectively, ATLAS-seq provides an efficient and broadly applicable technology to screen antigen-reactive TCRs for engineered T cell immunotherapy.

Local treatment of HVJ-E with T cell costimulatory molecule stimulation elicits systemic anti-tumor effects

The tumor-infiltrating lymphocyte (TIL) is a crucial factor in controlling tumor growth. A therapeutic method activating TIL is desired for treating patients with metastatic tumors. Here, we show that treating a local tumor with a combination therapy of UV-irradiated hemagglutinating virus of Japan envelope (HVJ-E) plus agonist antibodies, including OX40, against T cell costimulatory molecules induces systemic anti-tumor effects in a T cell-dependent manner in multiple cancer cell lines. Transcriptome and T cell receptor repertoire analyses revealed that HVJ-E + anti-OX40 antibody treatment activates CD4 and CD8 T cells and promotes T cell trafficking between tumors. These systemic anti-tumor effects required an association between Nkg2d and Nkg2d ligands. Our findings provide insights into how systemic anti-tumor effects are induced and may help the development of therapeutic strategies for eliciting such effects.

Cracking the code of adaptive immunity: The role of computational tools

In recent years, the advances in high-throughput and deep sequencing have generated a diverse amount of adaptive immune repertoire data. This surge in data has seen a proportional increase in computational methods aimed to characterize T cell receptor (TCR) repertoires. In this perspective, we will provide a brief commentary on the various domains of TCR repertoire analysis, their respective computational methods, and the ongoing challenges. Given the breadth of methods and applications of TCR analysis, we will focus our perspective on sequence-based computational methods.

SLAM/SAP signaling regulates discrete γδ T cell developmental checkpoints and shapes the innate-like γδ TCR repertoire

During thymic development, most γδ T cells acquire innate-like characteristics that are critical for their function in tumor surveillance, infectious disease, and tissue repair. The mechanisms, however, that regulate γδ T cell developmental programming remain unclear. Recently, we demonstrated that the SLAM/SAP signaling pathway regulates the development and function of multiple innate-like γδ T cell subsets. Here, we used a single-cell proteogenomics approach to identify SAP-dependent developmental checkpoints and to define the SAP-dependent γδ TCR repertoire in mice. SAP deficiency resulted in both a significant loss of an immature Gzma+Blk+Etv5+Tox2+ γδT17 precursor population and a significant increase in Cd4+Cd8+Rorc+Ptcra+Rag1+ thymic γδ T cells. SAP-dependent diversion of embryonic day 17 thymic γδ T cell clonotypes into the αβ T cell developmental pathway was associated with a decreased frequency of mature clonotypes in neonatal thymus, and an altered γδ TCR repertoire in the periphery. Finally, we identify TRGV4/TRAV13-4(DV7)-expressing T cells as a novel, SAP-dependent Vγ4 γδT1 subset. Together, the data support a model in which SAP-dependent γδ/αβ T cell lineage commitment regulates γδ T cell developmental programming and shapes the γδ TCR repertoire.

Evolution of genome and immunogenome in esophageal squamous cell carcinomas driven by neoadjuvant chemoradiotherapy

Neoadjuvant chemoradiotherapy (NCRT) followed by surgery is a standard treatment for locally advanced esophageal squamous cell carcinomas (ESCCs). However, the evolution of genome and immunogenome in ESCCs driven by NCRT remains incompletely elucidated. We performed whole-exome sequencing of 51 ESCC tumors collected before and after NCRT, 36 of which were subjected to transcriptome sequencing. Clonal analysis identified clonal extinction in 13 ESCC patients wherein all pre-NCRT clones disappeared after NCRT, and clonal persistence in 9 patients wherein clones endured following NCRT. The clone-persistent patients showed higher pre-NCRT genomic intratumoral heterogeneity and worse prognosis than the clone-extinct ones. In contrast to the clone-extinct patients, the clone-persistent patients demonstrated a high proportion of subclonal neoantigens within pre-treatment specimens. Transcriptome analysis revealed increased immune infiltrations and up-regulated immune-related pathways after NCRT, especially in the clone-extinct patients. The number of T cell receptor-neoantigen interactions was higher in the clone-extinct patients than in the clone-persistent ones. The decrease in T cell repertoire evenness positively correlated to the decreased number of clonal neoantigens after NCRT, especially in the clone-extinct patients. In conclusion, we identified two prognosis-related clonal dynamic modes driven by NCRT in ESCCs. This study extended our knowledge of the ESCC genome and immunogenome evolutions driven by NCRT.

The T-cell receptor β chain CDR3 insights of bovine liver immune repertoire under heat stress

OBJECTIVE

The liver plays a dual role in regulating temperature and immune responses. Examining the influence of heat stress (HS) on liver T cells contributes significantly to understanding the intricate interplay between the immune system and hepatic tissues under thermal stress. This study focused on investigating the characteristics of the T-cell receptor (TCR) β chain CDR3 repertoire in bovine liver samples under both HS and pairfed (PF) environmental conditions.

METHODS

Sequencing data from six samples sourced from the GEO database underwent annotation. Utilizing immunarch and VDJtool software, the study conducted comprehensive analyses encompassing basic evaluation, clonality assessment, immune repertoire comparison, diversity estimation, gene usage profiling, VJ gene segment pairing scrutiny, clonal tracking, and Kmers analysis.

RESULTS

All four TCR chains, namely α, β, γ, and δ, were detected, with the α chains exhibiting the highest detection frequency, followed closely by the β chains. The prevalence of αβ TCRs in bovine liver samples underscored their crucial role in governing hepatic tissue's physiological functions. The TCR β CDR3 repertoire showcased substantial inter-individual variability, featuring diverse clonotypes exhibiting distinct amino acid lengths. Intriguingly, HS cattle displayed heightened diversity and clonality, suggesting potential peripheral T cell migration into the liver under environmental conditions. Notably, differential VJ gene pairings were observed in HS cattle compared to the PF, despite individual variations in V and J gene utilization. Additionally, while most high-frequency amino acid 5-mers remained consistent between the HS and PF, GELHF, and YDYHF were notably prevalent in the HS group. Across all samples, a prevalent trend of high-frequency 5mers skewed towards polar and hydrophobic amino acids was evident.

CONCLUSION

This study elucidates the characteristics of liver TCR β chain CDR3 repertoire under HS conditions, enhancing our understanding of HS implications.

Bulk T cell repertoire sequencing (TCR-Seq) is a powerful technology for understanding inflammation-mediated diseases

Multiple alterations in the T cell repertoire were identified in many chronic inflammatory diseases such as inflammatory bowel disease, multiple sclerosis, and rheumatoid arthritis, suggesting that T cells might, directly or indirectly, be implicated in these pathologies. This has sparked a deep interest in characterizing disease-associated T cell clonotypes as well as in identifying and quantifying their contribution to the pathophysiology of different autoimmune and inflammation-mediated diseases. Bulk T cell repertoire sequencing (TCR-Seq) has emerged as a powerful method to profile the T cell repertoire of a sample in a high throughput fashion. Given the increasing utilization of TCR-Seq, we aimed here to provide a comprehensive, up-to-date review of the method, its extensions, and its ability to investigate chronic and autoimmune diseases. Specifically, we started by introducing the immunological basis of TCR repertoire generation and features, followed by discussing different experimental approach to perform TCR-Seq, then we describe different methods and frameworks for analyzing the generated datasets. Subsequently, different experimental techniques for investigating the antigenicity of T cell clonotypes are described. Lastly, we discuss recent studies that utilized TCR-Seq to understand different inflammation-mediated diseases, discuss fallbacks of the technology and potential future directions to overcome current limitations.

Multiple myeloma long-term survivors exhibit sustained immune alterations decades after first-line therapy

The long-term consequences of cancer and its therapy on the patients' immune system years after cancer-free survival remain poorly understood. Here, we present an in-depth characterization of the bone marrow immune ecosystem of multiple myeloma long-term survivors, from initial diagnosis up to 17 years following a single therapy line and cancer-free survival. Using comparative single-cell analyses combined with molecular, genomic, and functional approaches, we demonstrate that multiple myeloma long-term survivors exhibit pronounced alterations in their bone marrow microenvironment associated with impaired immunity. These immunological alterations were frequently linked to an inflammatory immune circuit fueled by the long-term persistence or resurgence of residual myeloma cells. Notably, even in the complete absence of any detectable residual disease for decades, sustained changes in the immune system were observed, suggesting an irreversible 'immunological scarring' caused by the initial exposure to the cancer and therapy. Collectively, our study provides key insights into the molecular and cellular bone marrow ecosystem of long-term survivors of multiple myeloma, revealing both reversible and irreversible alterations in the immune compartment.

Early-life thymectomy leads to an increase of granzyme-producing γδ T cells in children with congenital heart disease

Congenital heart disease (CHD) is the most common birth defect in newborns, often requiring cardiac surgery with concomitant thymectomy that is known to increase disease susceptibility later in life. Studies of γδ T cells, which are one of the dominant T cells in the early fetal human thymus, are rare. Here, we provide a comprehensive analysis of the γδ T cell compartment via flow cytometry and next-generation sequencing in children and infants with CHD, who underwent cardiac surgery shortly after birth. A perturbation of the γδ T cell repertoire is evident, and Vδ1 T cell numbers are reduced. However, those cells that are present, do retain cytotoxicity. In contrast, GZMA+CD28+CD161hi innate effector Vγ9Vδ2 T cells are found in higher proportions. TCR-seq identifies an increase in TRDJ3+ γδ T cell clones in children with CHD, but not in a confirmatory group of neonates prior to CHD surgery, which overall points to a persistence of fetal-derived effector γδ T cells in children with CHD.

TCRosetta: An Integrated Analysis and Annotation Platform for T-cell Receptor Sequences

T cells and T-cell receptors (TCRs) are essential components of the adaptive immune system. Characterization of the TCR repertoire offers a promising and highly informative source for understanding the functions of T cells in the immune response and immunotherapy. Although TCR repertoire studies have attracted much attention, there are few online servers available for TCR repertoire analysis, especially for TCR sequence annotation or advanced analyses. Therefore, we developed TCRosetta, a comprehensive online server that integrates analytical methods for TCR repertoire analysis and visualization. TCRosetta combines general feature analysis, large-scale sequence clustering, network construction, peptide-TCR binding prediction, generation probability calculation, and k-mer motif analysis for TCR sequences, making TCR data analysis as simple as possible. The TCRosetta server accepts multiple input data formats and can analyze ∼ 20,000 TCR sequences in less than 3 min. TCRosetta is the most comprehensive web server available for TCR repertoire analysis and is freely available at https://guolab.wchscu.cn/TCRosetta/.

CXCR4 orchestrates the TOX-programmed exhausted phenotype of CD8+ T cells via JAK2/STAT3 pathway

Evidence from clinical trials suggests that CXCR4 antagonists enhance immunotherapy effectiveness in several cancers. However, the specific mechanisms through which CXCR4 contributes to immune cell phenotypes are not fully understood. Here, we employed single-cell transcriptomic analysis and identified CXCR4 as a marker gene in T cells, with CD8+PD-1high exhausted T (Tex) cells exhibiting high CXCR4 expression. By blocking CXCR4, the Tex phenotype was attenuated in vivo. Mechanistically, CXCR4-blocking T cells mitigated the Tex phenotype by regulating the JAK2-STAT3 pathway. Single-cell RNA/TCR/ATAC-seq confirmed that Cxcr4-deficient CD8+ T cells epigenetically mitigated the transition from functional to exhausted phenotypes. Notably, clinical sample analysis revealed that CXCR4+CD8+ T cells showed higher expression in patients with a non-complete pathological response. Collectively, these findings demonstrate the mechanism by which CXCR4 orchestrates CD8+ Tex cells and provide a rationale for combining CXCR4 antagonists with immunotherapy in clinical trials.

Sepsis shapes the human γδ TCR repertoire in an age- and pathogen-dependent manner

Sepsis affects 25 million children per year globally, leading to 2.9 million deaths and substantial disability in survivors. Extensive characterization of interactions between the host and bacteria in children is required to design novel preventive and therapeutic strategies tailored to this age group. Vγ9Vδ2 T cells are the first T cells generated in humans. These cells are defined by the expression of Vγ9Vδ2 T-cell receptors (TCRs, using the TRGV9 and TRDV2 gene segments), which react strongly against the prototypical bacterial phosphoantigen HMBPP. We investigated this reactivity by analyzing the TCR δ (TRD) repertoire in the blood of 76 children (0-16 years) with blood culture-proven bacterial sepsis caused by HMBPP-positive Escherichia coli or by HMBPP-negative Staphylococcus aureus or by HMBPP-negative Streptococcus pneumoniae. Strikingly, we found that S. aureus, and to a lesser extent E. coli but not S. pneumoniae, shaped the TRDV2 repertoire in young children (<2 years) but not in older children or adults. This dichotomy was due to the selective expansion of a fetal TRDV2 repertoire. Thus, young children possess fetal-derived Vγ9Vδ2 T cells that are highly responsive toward specific bacterial pathogens.

Age-associated clonal B cells drive B cell lymphoma in mice

Although cancer is an age-related disease, how the processes of aging contribute to cancer progression is not well understood. In this study, we uncovered how mouse B cell lymphoma develops as a consequence of a naturally aged system. We show here that this malignancy is associated with an age-associated clonal B cell (ACBC) population that likely originates from age-associated B cells. Driven by c-Myc activation, promoter hypermethylation and somatic mutations, IgM+ ACBCs clonally expand independently of germinal centers and show increased biological age. ACBCs become self-sufficient and support malignancy when transferred into young recipients. Inhibition of mTOR or c-Myc in old mice attenuates pre-malignant changes in B cells during aging. Although the etiology of mouse and human B cell lymphomas is considered distinct, epigenetic changes in transformed mouse B cells are enriched for changes observed in human B cell lymphomas. Together, our findings characterize the spontaneous progression of cancer during aging through both cell-intrinsic and microenvironmental changes and suggest interventions for its prevention.

Bioinformatics tools and resources for cancer and application

ABSTRACT

Tumor bioinformatics plays an important role in cancer research and precision medicine. The primary focus of traditional cancer research has been molecular and clinical studies of a number of fundamental pathways and genes. In recent years, driven by breakthroughs in high-throughput technologies, large-scale cancer omics data have accumulated rapidly. How to effectively utilize and share these data is particularly important. To address this crucial task, many computational tools and databases have been developed over the past few years. To help researchers quickly learn and understand the functions of these tools, in this review, we summarize publicly available bioinformatics tools and resources for pan-cancer multi-omics analysis, regulatory analysis of tumorigenesis, tumor treatment and prognosis, immune infiltration analysis, immune repertoire analysis, cancer driver gene and driver mutation analysis, and cancer single-cell analysis, which may further help researchers find more suitable tools for their research.

Development and characterization of human T-cell receptor (TCR) alpha and beta clones' library as biological standards and resources for TCR sequencing and engineering

Characterization of T-cell receptors (TCRs) repertoire was revolutionized by next-generation sequencing technologies; however, standardization using biological controls to facilitate precision of current alignment and assembly tools remains a challenge. Additionally, availability of TCR libraries for off-the-shelf cloning and engineering TCR-specific T cells is a valuable resource for TCR-based immunotherapies. We established nine human TCR α and β clones that were evaluated using the 5'-rapid amplification of cDNA ends-like RNA-based TCR sequencing on the Illumina platform. TCR sequences were extracted and aligned using MiXCR, TRUST4, and CATT to validate their sensitivity and specificity and to validate library preparation methods. The correlation between actual and expected TCR ratios within libraries confirmed accuracy of the approach. Our findings established the development of biological standards and library of TCR clones to be leveraged in TCR sequencing and engineering. The remaining human TCR clones' libraries for a more diverse biological control will be generated.

The T-cell repertoire of Spanish patients with COVID-19 as a strategy to link T-cell characteristics to the severity of the disease

BACKGROUND

The architecture and dynamics of T cell populations are critical in orchestrating the immune response to SARS-CoV-2. In our study, we used T Cell Receptor sequencing (TCRseq) to investigate TCR repertoires in 173 post-infection COVID-19 patients.

METHODS

The cohort included 98 mild and 75 severe cases with a median age of 53. We amplified and sequenced the TCR β chain Complementary Determining Region 3 (CDR3b) and performed bioinformatic analyses to assess repertoire diversity, clonality, and V/J allelic usage between age, sex and severity groups. CDR3b amino acid sequence inference was performed by clustering structural motifs and filtering validated reactive CDR3b to COVID-19.

RESULTS

Our results revealed a pronounced decrease in diversity and an increase in clonal expansion in the TCR repertoires of severe COVID-19 patients younger than 55 years old. These results reflect the observed trends in patients older than 55 years old (both mild and severe). In addition, we identified a significant reduction in the usage of key V alleles (TRBV14, TRBV19, TRBV15 and TRBV6-4) associated with disease severity. Notably, severe patients under 55 years old had allelic patterns that resemble those over 55 years old, accompanied by a skewed frequency of COVID-19-related motifs.

CONCLUSIONS

Present results suggest that severe patients younger than 55 may have a compromised TCR repertoire contributing to a worse disease outcome.

T-Cell Receptors Cross-Reactive to Coronaviral Epitopes Homologous to the SPR Peptide

The COVID-19 pandemic caused by the rapid spread of the novel coronavirus SARS-CoV-2, has promoted an interest in studying the T-cell immune response. It was found that the polyclonal and cross-reactive T-cell response against seasonal coronaviruses and other SARS-CoV-2 strains reduced disease severity. We investigated the immunodominant T-cell epitope SPRWYFYYYL from the nucleocapsid protein of SARS-CoV-2. The immune response to this epitope is characterized by the formation of highly homologous (convergent) receptors that have been found in the T-cell receptor (TCR) repertoires of different individuals. This epitope belongs to a group of highly conserved peptides that are rarely mutated in novel SARS-CoV-2 strains and are homologous to the epitopes of seasonal coronaviruses. It has been suggested that the cross-reactive response to homologous peptides contributes to the reduction of COVID-19 severity. However, some investigators have questioned this hypothesis, suggesting that the low affinity of the cross-reactive receptors reduces the strength of the immune response. The aim of this study was to evaluate the effect of amino acid substitutions in the SPR epitope on its binding affinity to specific TCRs. For this, we performed antigen-dependent cellular expansions were performed using samples from four COVID-19-transfected donors and sequenced their TCR repertoires. The resulting SPR-specific repertoire of β-chains in TCRs had a greater sequence diversity than the repertoire of α-chains. However, the TCR repertoires of all four donors contained public receptors, three of which were cloned and used to generate the Jurkat E6-1 TPR cell line. Only one of these receptors was activated by the SPR peptide and recognized with the same affinity by its mutant homologue LPRWYFYYY from seasonal coronaviruses. This indicates that the presence of the mutation did not affect the strength of the immune response, which may explain why the cross-reactive response to the SPR epitope is so frequent and contributes positively to COVID-19 infection.

Influenza vaccination stimulates maturation of the human T follicular helper cell response

The differentiation and specificity of human CD4+ T follicular helper cells (TFH cells) after influenza vaccination have been poorly defined. Here we profiled blood and draining lymph node (LN) samples from human volunteers for over 2 years after two influenza vaccines were administered 1 year apart to define the evolution of the CD4+ TFH cell response. The first vaccination induced an increase in the frequency of circulating TFH (cTFH) and LN TFH cells at week 1 postvaccination. This increase was transient for cTFH cells, whereas the LN TFH cells further expanded during week 2 and remained elevated in frequency for at least 3 months. We observed several distinct subsets of TFH cells in the LN, including pre-TFH cells, memory TFH cells, germinal center (GC) TFH cells and interleukin-10+ TFH cell subsets beginning at baseline and at all time points postvaccination. The shift toward a GC TFH cell phenotype occurred with faster kinetics after the second vaccine compared to the first vaccine. We identified several influenza-specific TFH cell clonal lineages, including multiple responses targeting internal influenza virus proteins, and found that each TFH cell state was attainable within a clonal lineage. Thus, human TFH cells form a durable and dynamic multitissue network.

Integrative alternative splicing analysis reveals new prognosis signature in B-cell acute lymphoblastic leukemia

The dysregulation of alternative splicing (AS) is increasingly recognized as a pivotal player in the pathogenesis, progression, and treatment resistance of B-cell acute lymphoblastic leukemia (B-ALL). Despite its significance, the clinical implications of AS events in B-ALL remain largely unexplored. This study developed a prognostic model based on 18 AS events (18-AS), derived from a meticulous integration of bioinformatics methodologies and advanced machine learning algorithms. The 18-AS signature observed in B-ALL distinctly categorized patients into different groups with significant differences in immune infiltration, V(D)J rearrangement, drug sensitivity, and immunotherapy outcomes. Patients classified within the high 18-AS group exhibited lower immune infiltration scores, poorer chemo- and immune-therapy responses, and worse overall survival, underscoring the model's potential in refining therapeutic strategies. To validate the clinical applicability of the 18-AS, we established an SF-AS regulatory network and identified candidate drugs. More importantly, we conducted in vitro cell proliferation assays to confirm our analysis, demonstrating that the High-18AS cell line (SUP-B15) exhibited significantly enhanced sensitivity to Dasatinib, Dovitinib, and Midostaurin compared to the Low-18AS cell line (REH). These findings reveal AS events as novel prognostic biomarkers and therapeutic targets, advancing personalized treatment strategies in B-ALL management.

Sculpting the tumour microenvironment by combining radiotherapy and ATR inhibition for curative-intent adjuvant immunotherapy

The combination of radiotherapy/chemoradiotherapy and immune checkpoint blockade can result in poor outcomes in patients with locally advanced head and neck squamous cell carcinoma (HNSCC). Here, we show that combining ATR inhibition (ATRi) with radiotherapy (RT) increases the frequency of activated NKG2A+PD-1+ T cells in animal models of HNSCC. Compared with the ATRi/RT treatment regimen alone, the addition of simultaneous NKG2A and PD-L1 blockade to ATRi/RT, in the adjuvant, post-radiotherapy setting induces a robust antitumour response driven by higher infiltration and activation of cytotoxic T cells in the tumour microenvironment. The efficacy of this combination relies on CD40/CD40L costimulation and infiltration of activated, proliferating memory CD8+ and CD4+ T cells with persistent or new T cell receptor (TCR) signalling, respectively. We also observe increased richness in the TCR repertoire and emergence of numerous and large TCR clonotypes that cluster based on antigen specificity in response to NKG2A/PD-L1/ATRi/RT. Collectively, our data point towards potential combination approaches for the treatment of HNSCC.

Genomic analysis and identification of a novel superantigen, SargEY, in Staphylococcus argenteus isolated from atopic dermatitis lesions

During surveillance of Staphylococcus aureus in lesions from patients with atopic dermatitis (AD), we isolated Staphylococcus argenteus, a species registered in 2011 as a new member of the genus Staphylococcus and previously considered a lineage of S. aureus. Genome sequence comparisons between S. argenteus isolates and representative S. aureus clinical isolates from various origins revealed that the S. argenteus genome from AD patients closely resembles that of S. aureus causing skin infections. We previously reported that 17%-22% of S. aureus isolated from skin infections produce staphylococcal enterotoxin Y (SEY), which predominantly induces T-cell proliferation via the T-cell receptor (TCR) Vα pathway. Complete genome sequencing of S. argenteus isolates revealed a gene encoding a protein similar to superantigen SEY, designated as SargEY, on its chromosome. Population structure analysis of S. argenteus revealed that these isolates are ST2250 lineage, which was the only lineage positive for the SEY-like gene among S. argenteus. Recombinant SargEY demonstrated immunological cross-reactivity with anti-SEY serum. SargEY could induce proliferation of human CD4+ and CD8+ T cells, as well as production of TNF-α and IFN-γ. SargEY showed emetic activity in a marmoset monkey model. SargEY and SET (a phylogenetically close but uncharacterized SE) revealed their dependency on TCR Vα in inducing human T-cell proliferation. Additionally, TCR sequencing revealed other previously undescribed Vα repertoires induced by SEH. SargEY and SEY may play roles in exacerbating the respective toxin-producing strains in AD.

IMPORTANCE

Staphylococcus aureus is frequently isolated from active lesions of atopic dermatitis (AD) patients. We reported that 17%-22% of S. aureus isolated from AD patients produced a novel superantigen staphylococcal enterotoxin Y (SEY). Unlike many S. aureus superantigens that activate T cells via T-cell receptor (TCR) Vß, SEY activates T cells via TCR Vα and stimulates cytokine secretion. Staphylococcus argenteus was isolated from AD patients during the surveillance for S. aureus. Phylogenetic comparison of the genome indicated that the isolate was very similar to S. aureus causing skin infections. The isolate encoded a SEY-like protein, designated SargEY, which, like SEY, activated T cells via the TCR Vα. ST2250 is the only lineage positive for SargEY gene. ST2250 S. argenteus harboring a superantigen SargEY gene may be a novel staphylococcal clone that infects human skin and is involved in the exacerbation of AD.

Microbiota dictate T cell clonal selection to augment graft-versus-host disease after stem cell transplantation

Allogeneic T cell expansion is the primary determinant of graft-versus-host disease (GVHD), and current dogma dictates that this is driven by histocompatibility antigen disparities between donor and recipient. This paradigm represents a closed genetic system within which donor T cells interact with peptide-major histocompatibility complexes (MHCs), though clonal interrogation remains challenging due to the sparseness of the T cell repertoire. We developed a Bayesian model using donor and recipient T cell receptor (TCR) frequencies in murine stem cell transplant systems to define limited common expansion of T cell clones across genetically identical donor-recipient pairs. A subset of donor CD4+ T cell clonotypes differentially expanded in identical recipients and were microbiota dependent. Microbiota-specific T cells augmented GVHD lethality and could target microbial antigens presented by gastrointestinal epithelium during an alloreactive response. The microbiota serves as a source of cognate antigens that contribute to clonotypic T cell expansion and the induction of GVHD independent of donor-recipient genetics.

Selection and characterization of human scFvs targeting the SARS-CoV-2 nucleocapsid protein isolated from antibody libraries of COVID-19 patients

In 2019, the novel SARS-CoV-2 coronavirus emerged in China, causing the pneumonia named COVID-19. At the beginning, all research efforts were focused on the spike (S) glycoprotein. However, it became evident that the nucleocapsid (N) protein is pivotal in viral replication, genome packaging and evasion of the immune system, is highly immunogenic, which makes it another compelling target for antibody development alongside the spike protein. This study focused on the construction of single chain fragments variable (scFvs) libraries from SARS-CoV-2-infected patients to establish a valuable, immortalized and extensive antibodies source. We used the Intracellular Antibody Capture Technology to select a panel of scFvs against the SARS-CoV-2 N protein. The whole panel of scFv was expressed and characterized both as intrabodies and recombinant proteins. ScFvs were then divided into 2 subgroups: those that exhibited high binding activity to N protein when expressed in yeast or in mammalian cells as intrabodies, and those purified as recombinant proteins, displaying affinity for recombinant N protein in the nanomolar range. This panel of scFvs against the N protein represents a novel platform for research and potential diagnostic applications.

Restoration of T and B Cell Differentiation after RAG1 Gene Transfer in Human RAG1 Defective Hematopoietic Stem Cells

Recombinase-activating gene (RAG)-deficient SCID patients lack B and T lymphocytes due to the inability to rearrange immunoglobulin and T cell receptor genes. The two RAG genes act as a required dimer to initiate gene recombination. Gene therapy is a valid treatment alternative for RAG-SCID patients who lack a suitable bone marrow donor, but developing such therapy for RAG1/2 has proven challenging. Using a clinically approved lentiviral vector with a codon-optimized RAG1 gene, we report here preclinical studies using CD34+ cells from four RAG1-SCID patients. We used in vitro T cell developmental assays and in vivo assays in xenografted NSG mice. The RAG1-SCID patient CD34+ cells transduced with the RAG1 vector and transplanted into NSG mice led to restored human B and T cell development. Together with favorable safety data on integration sites, these results substantiate an ongoing phase I/II clinical trial for RAG1-SCID.

SLAM/SAP signaling regulates discrete γδ T cell developmental checkpoints and shapes the innate-like γδ TCR repertoire

During thymic development, most γδ T cells acquire innate-like characteristics that are critical for their function in tumor surveillance, infectious disease, and tissue repair. The mechanisms, however, that regulate γδ T cell developmental programming remain unclear. Recently, we demonstrated that the SLAM-SAP signaling pathway regulates the development and function of multiple innate-like γδ T cell subsets. Here, we used a single-cell proteogenomics approach to identify SAP-dependent developmental checkpoints and to define the SAP-dependent γδ TCR repertoire. SAP deficiency resulted in both a significant loss of an immature Gzma + Blk + Etv5 + Tox2 + γδT17 precursor population, and a significant increase in Cd4 + Cd8+ Rorc + Ptcra + Rag1 + thymic γδ T cells. SAP-dependent diversion of embryonic day 17 thymic γδ T cell clonotypes into the αβ T cell developmental pathway was associated with a decreased frequency of mature clonotypes in neonatal thymus, and an altered γδ TCR repertoire in the periphery. Finally, we identify TRGV4/TRAV13-4(DV7)-expressing T cells as a novel, SAP-dependent Vγ4 γδT1 subset. Together, the data suggest that SAP-dependent γδ/αβ T cell lineage commitment regulates γδ T cell developmental programming and shapes the γδ TCR repertoire.

Harnessing HetHydrogel: A Universal Platform to Dropletize Single-Cell Multiomics

A universal platform is developed for dropletizing single cell plate-based multiomic assays, consisting of three main pillars: a miniaturized open Heterogeneous Hydrogel reactor (abbreviated HetHydrogel) for multi-step biochemistry, its tunable permeability that allows Tn5 tagmentation, and single cell droplet barcoding. Through optimizing the HetHydrogel manufacturing procedure, the chemical composition, and cell permeation conditions, simultaneous high-throughput mitochondrial DNA genotyping and chromatin profiling at the single-cell level are demonstrated using a mixed-species experiment. This platform offers a powerful way to investigate the genotype-phenotype relationships of various mtDNA mutations in biological processes. The HetHydrogel platform is believed to have the potential to democratize droplet technologies, upgrading a whole range of plate-based single cell assays to high throughput format.

The Immunogenomic Landscape of Peripheral High-Dose IL-2 Pharmacodynamics in Patients with Metastatic Renal Cell Carcinoma: A Benchmark for Next-Generation IL-2-Based Immunotherapies

High-dose (HD) IL-2 was the first immuno-oncology agent approved for treating advanced renal cell carcinoma and metastatic melanoma, but its use was limited because of substantial toxicities. Multiple next-generation IL-2 agents are being developed to improve tolerability. However, a knowledge gap still exists for the genomic markers that define the target pharmacology for HD IL-2 itself. In this retrospective observational study, we collected PBMC samples from 23 patients with metastatic renal cell carcinoma who were treated with HD IL-2 between 2009 and 2015. We previously reported the results of flow cytometry analyses. In this study, we report the results of our RNA-sequencing immunogenomic survey, which was performed on bulk PBMC samples from immediately before (day 1), during (day 3), and after treatment (day 5) in cycle 1 and/or cycle 2 of the first course of HD IL-2. As part of a detailed analysis of immunogenomic response to HD IL-2 treatment, we analyzed the changes in individual genes and immune gene signatures. By day 3, most lymphoid cell types had transiently decreased, whereas myeloid transcripts increased. Although most genes and/or signatures generally returned to pretreatment expression levels by day 5, certain ones representative of B cell, NK cell, and T cell proliferation and effector functions continued to increase, along with B cell (but not T cell) oligoclonal expansion. Regulatory T cells progressively expanded during and after treatment. They showed strong negative correlation with myeloid effector cells. This detailed RNA-sequencing immunogenomic survey of IL-2 pharmacology complements results of prior flow cytometry analyses. These data provide valuable pharmacological context for assessing PBMC gene expression data from patients dosed with IL-2-related compounds that are currently in development.

Long-lived central memory γδ T cells confer protection against murine cytomegalovirus reinfection

The involvement of γδ TCR-bearing lymphocytes in immunological memory has gained increasing interest due to their functional duality between adaptive and innate immunity. γδ T effector memory (TEM) and central memory (TCM) subsets have been identified, but their respective roles in memory responses are poorly understood. In the present study, we used subsequent mouse cytomegalovirus (MCMV) infections of αβ T cell deficient mice in order to analyze the memory potential of γδ T cells. As for CMV-specific αβ T cells, MCMV induced the accumulation of cytolytic, KLRG1+CX3CR1+ γδ TEM that principally localized in infected organ vasculature. Typifying T cell memory, γδ T cell expansion in organs and blood was higher after secondary viral challenge than after primary infection. Viral control upon MCMV reinfection was prevented when masking γδ T-cell receptor, and was associated with a preferential amplification of private and unfocused TCR δ chain repertoire composed of a combination of clonotypes expanded post-primary infection and, more unexpectedly, of novel expanded clonotypes. Finally, long-term-primed γδ TCM cells, but not γδ TEM cells, protected T cell-deficient hosts against MCMV-induced death upon adoptive transfer, probably through their ability to survive and to generate TEM in the recipient host. This better survival potential of TCM cells was confirmed by a detailed scRNASeq analysis of the two γδ T cell memory subsets which also revealed their similarity to classically adaptive αβ CD8 T cells. Overall, our study uncovered memory properties of long-lived TCM γδ T cells that confer protection in a chronic infection, highlighting the interest of this T cell subset in vaccination approaches.

γδ T cell profiling in a cohort of preterm infants reveals elevated frequencies of CD83+ γδ T cells in sepsis

Preterm infants are at high risk of developing neonatal sepsis. γδ T cells are thought to be an important set of effector cells in neonates. Here, γδ T cells were investigated in a longitudinal cohort of preterm neonates using next-generation sequencing, flow cytometry, and functional assays. During the first year of life, the Vγ9Vδ2 T cell subset showed dynamic phenotypic changes and elevated levels of fetal-derived Vγ9Vδ2 T cells were evident in infants with sepsis. Single-cell transcriptomics identified HLA-DRhiCD83+ γδ T cells in neonatal sepsis, which expressed genes related to antigen presentation. In vitro assays showed that CD83 was expressed on activated Vγ9Vδ2 T cells in preterm and term neonates, but not in adults. In contrast, activation of adult Vγ9Vδ2 T cells enhanced CD86 expression, which was presumably the key receptor to induce CD4 T cell proliferation. Together, we provide a map of the maturation of γδ T cells after preterm birth and highlight their phenotypic diversity in infections.

Deciphering Abdominal Aortic Diseases Through T-Cell Clonal Repertoire of Perivascular Adipose Tissue

BACKGROUND

Recent studies suggest that immune-mediated inflammation of perivascular adipose tissue of abdominal aortic aneurysms (AAAs) contributes to disease development and progression. Whether the perivascular adipose tissue of AAA is characterized by a specific adaptive immune signature remains unknown.

METHODS AND RESULTS

To investigate this hypothesis, we sequenced the T-cell receptor β-chain in the perivascular adipose tissue of patients with AAA and compared it with patients with aortic occlusive disease, who share the former anatomical site of the lesion and risk factors but differ in pathogenic mechanisms. Our results demonstrate that patients with AAA have a lower repertoire diversity than those with aortic occlusive disease and significant differences in variable/joining gene segment usage. Furthermore, we identified a set of 7 public T-cell receptor β-chain clonotypes that distinguished AAA and aortic occlusive disease with very high accuracy. We also found that the T-cell receptor β-chain repertoire differentially characterizes small and large AAAs (aortic diameter<55 mm and ≥55 mm, respectively).

CONCLUSIONS

This work supports the hypothesis that T cell-mediated immunity is fundamental in AAA pathogenesis and opens up new clinical perspectives.

Immunotoxin-mediated depletion of Gag-specific CD8+ T cells undermines natural control of SIV

Antibody-mediated depletion studies have demonstrated that CD8+ T cells are required for effective immune control of SIV. However, this approach is potentially confounded by several factors, including reactive CD4+ T cell proliferation, and provides no information on epitope specificity, a likely determinant of CD8+ T cell efficacy. We circumvented these limitations by selectively depleting CD8+ T cells specific for the Gag epitope CTPYDINQM (CM9) via the administration of immunotoxin-conjugated tetrameric complexes of CM9/Mamu-A*01. Immunotoxin administration effectively depleted circulating but not tissue-localized CM9-specific CD8+ T cells, akin to the bulk depletion pattern observed with antibodies directed against CD8. However, we found no evidence to indicate that circulating CM9-specific CD8+ T cells suppressed viral replication in Mamu-A*01+ rhesus macaques during acute or chronic progressive infection with a pathogenic strain of SIV. This observation extended to macaques with established infection during and after continuous antiretroviral therapy. In contrast, natural controller macaques experienced dramatic increases in plasma viremia after immunotoxin administration, highlighting the importance of CD8+ T cell-mediated immunity against CM9. Collectively, these data showed that CM9-specific CD8+ T cells were necessary but not sufficient for robust immune control of SIV in a nonhuman primate model and, more generally, validated an approach that could inform the design of next-generation vaccines against HIV-1.

Somatic mutations associate with clonal expansion of CD8+ T cells

Somatic mutations in T cells can cause cancer but also have implications for immunological diseases and cell therapies. The mutation spectrum in nonmalignant T cells is unclear. Here, we examined somatic mutations in CD4+ and CD8+ T cells from 90 patients with hematological and immunological disorders and used T cell receptor (TCR) and single-cell sequencing to link mutations with T cell expansions and phenotypes. CD8+ cells had a higher mutation burden than CD4+ cells. Notably, the biggest variant allele frequency (VAF) of non-synonymous variants was higher than synonymous variants in CD8+ T cells, indicating non-random occurrence. The non-synonymous VAF in CD8+ T cells strongly correlated with the TCR frequency, but not age. We identified mutations in pathways essential for T cell function and often affected lymphoid neoplasia. Single-cell sequencing revealed cytotoxic TEMRA phenotypes of mutated T cells. Our findings suggest that somatic mutations contribute to CD8+ T cell expansions without malignant transformation.

KSHV infection of B cells primes protective T cell responses in humanized mice

Kaposi sarcoma associated herpesvirus (KSHV) is associated with around 1% of all human tumors, including the B cell malignancy primary effusion lymphoma (PEL), in which co-infection with the Epstein Barr virus (EBV) can almost always be found in malignant cells. Here, we demonstrate that KSHV/EBV co-infection of mice with reconstituted human immune systems (humanized mice) leads to IgM responses against both latent and lytic KSHV antigens, and expansion of central and effector memory CD4+ and CD8+ T cells. Among these, KSHV/EBV dual-infection allows for the priming of CD8+ T cells that are specific for the lytic KSHV antigen K6 and able to kill KSHV/EBV infected B cells. This suggests that K6 may represent a vaccine antigen for the control of KSHV and its associated pathologies in high seroprevalence regions, such as Sub-Saharan Africa.