Single-cell lineage mapping of a diverse virus-specific naive CD4 T cell repertoire

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.

Plasmodium infection induces phenotypic, clonal, and spatial diversity among differentiating CD4+ T cells

Naive CD4+ T cells must differentiate in order to orchestrate immunity to Plasmodium, yet understanding of their emerging phenotypes, clonality, spatial distributions, and cellular interactions remains incomplete. Here, we observe that splenic polyclonal CD4+ T cells differentiate toward T helper 1 (Th1) and T follicular helper (Tfh)-like states and exhibit rarer phenotypes not elicited among T cell receptor (TCR) transgenic counterparts. TCR clones present at higher frequencies exhibit Th1 skewing, suggesting that variation in major histocompatibility complex class II (MHC-II) interaction influences proliferation and Th1 differentiation. To characterize CD4+ T cell interactions, we map splenic microarchitecture, cellular locations, and molecular interactions using spatial transcriptomics at near single-cell resolution. Tfh-like cells co-locate with stromal cells in B cell follicles, while Th1 cells in red pulp co-locate with activated monocytes expressing multiple chemokines and MHC-II. Spatial mapping of individual transcriptomes suggests that proximity to chemokine-expressing monocytes correlates with stronger effector phenotypes in Th1 cells. Finally, CRISPR-Cas9 gene disruption reveals a role for CCR5 in promoting clonal expansion and Th1 differentiation. A database of cellular locations and interactions is presented: https://haquelab.mdhs.unimelb.edu.au/spatial_gui/.

Bhlhe40 Promotes CD4+ T Helper 1 Cell and Suppresses T Follicular Helper Cell Differentiation during Viral Infection

In response to acute infection, naive CD4+ T cells primarily differentiate into T helper 1 (Th1) or T follicular helper (Tfh) cells that play critical roles in orchestrating cellular or humoral arms of immunity, respectively. However, despite the well established role of T-bet and BCL-6 in driving Th1 and Tfh cell lineage commitment, respectively, whether additional transcriptional circuits also underlie the fate bifurcation of Th1 and Tfh cell subsets is not fully understood. In this article, we study how the transcriptional regulator Bhlhe40 dictates the Th1/Tfh differentiation axis in mice. CD4+ T cell-specific deletion of Bhlhe40 abrogates Th1 but augments Tfh differentiation. We also assessed an increase in germinal center B cells and Ab production, suggesting that deletion of Bhlhe40 in CD4+ T cells not only alters Tfh differentiation but also their capacity to provide help to B cells. To identify molecular mechanisms by which Bhlhe40 regulates Th1 versus Tfh lineage choice, we first performed epigenetic profiling in the virus specific Th1 and Tfh cells following LCMV infection, which revealed distinct promoter and enhancer activities between the two helper cell lineages. Furthermore, we identified that Bhlhe40 directly binds to cis-regulatory elements of Th1-related genes such as Tbx21 and Cxcr6 to activate their expression while simultaneously binding to regions of Tfh-related genes such as Bcl6 and Cxcr5 to repress their expression. Collectively, our data suggest that Bhlhe40 functions as a transcription activator to promote Th1 cell differentiation and a transcription repressor to suppress Tfh cell differentiation.

Single-cell RNA sequencing reveals the dynamics and heterogeneity of lymph node immune cells during acute and chronic viral infections

Introduction

The host immune response determines the differential outcome of acute or chronic viral infections. The comprehensive comparison of lymphoid tissue immune cells at the single-cell level between acute and chronic viral infections is largely insufficient.

Methods

To explore the landscape of immune responses to acute and chronic viral infections, single-cell RNA sequencing(scRNA-seq), scTCR-seq and scBCR-seq were utilized to evaluate the longitudinal dynamics and heterogeneity of lymph node CD45+ immune cells in mouse models of acute (LCMV Armstrong) and chronic (LCMV clone 13) viral infections.

Results

In contrast with acute viral infection, chronic viral infection distinctly induced more robust NK cells and plasma cells at the early stage (Day 4 post-infection) and acute stage (Day 8 post-infection), respectively. Moreover, chronic viral infection exerted decreased but aberrantly activated plasmacytoid dendritic cells (pDCs) at the acute phase. Simultaneously, there were significantly increased IgA+ plasma cells (MALT B cells) but differential usage of B-cell receptors in chronic infection. In terms of T-cell responses, Gzma-high effector-like CD8+ T cells were significantly induced at the early stage in chronic infection, which showed temporally reversed gene expression throughout viral infection and the differential usage of the most dominant TCR clonotype. Chronic infection also induced more robust CD4+ T cell responses, including follicular helper T cells (Tfh) and regulatory T cells (Treg). In addition, chronic infection compromised the TCR diversity in both CD8+ and CD4+ T cells.

Discussion

In conclusion, gene expression and TCR/BCR immune repertoire profiling at the single-cell level in this study provide new insights into the dynamic and differential immune responses to acute and chronic viral infections.

Machine learning analysis of the T cell receptor repertoire identifies sequence features of self-reactivity

The T cell receptor (TCR) determines specificity and affinity for both foreign and self-peptides presented by the major histocompatibility complex (MHC). Although the strength of TCR interactions with self-pMHC impacts T cell function, it has been challenging to identify TCR sequence features that predict T cell fate. To discern patterns distinguishing TCRs from naive CD4+ T cells with low versus high self-reactivity, we used data from 42 mice to train a machine learning (ML) algorithm that identifies population-level differences between TCRβ sequence sets. This approach revealed that weakly self-reactive T cell populations were enriched for longer CDR3β regions and acidic amino acids. We tested our ML predictions of self-reactivity using retrogenic mice with fixed TCRβ sequences. Extrapolating our analyses to independent datasets, we predicted high self-reactivity for regulatory T cells and slightly reduced self-reactivity for T cells responding to chronic infections. Our analyses suggest a potential trade-off between TCR repertoire diversity and self-reactivity. A record of this paper's transparent peer review process is included in the supplemental information.

A spatial sequencing atlas of age-induced changes in the lung during influenza infection

Influenza virus infection causes increased morbidity and mortality in the elderly. Aging impairs the immune response to influenza, both intrinsically and because of altered interactions with endothelial and pulmonary epithelial cells. To characterize these changes, we performed single-cell RNA sequencing (scRNA-seq), spatial transcriptomics, and bulk RNA sequencing (bulk RNA-seq) on lung tissue from young and aged female mice at days 0, 3, and 9 post-influenza infection. Our analyses identified dozens of key genes differentially expressed in kinetic, age-dependent, and cell type-specific manners. Aged immune cells exhibited altered inflammatory, memory, and chemotactic profiles. Aged endothelial cells demonstrated characteristics of reduced vascular wound healing and a prothrombotic state. Spatial transcriptomics identified novel profibrotic and antifibrotic markers expressed by epithelial and non-epithelial cells, highlighting the complex networks that promote fibrosis in aged lungs. Bulk RNA-seq generated a timeline of global transcriptional activity, showing increased expression of genes involved in inflammation and coagulation in aged lungs. Our work provides an atlas of high-throughput sequencing methodologies that can be used to investigate age-related changes in the response to influenza virus, identify novel cell-cell interactions for further study, and ultimately uncover potential therapeutic targets to improve health outcomes in the elderly following influenza infection.

Molecular, metabolic, and functional CD4 T cell paralysis in the lymph node impedes tumor control

CD4 T cells are central effectors of anti-cancer immunity and immunotherapy, yet the regulation of CD4 tumor-specific T (TTS) cells is unclear. We demonstrate that CD4 TTS cells are quickly primed and begin to divide following tumor initiation. However, unlike CD8 TTS cells or exhaustion programming, CD4 TTS cell proliferation is rapidly frozen in place by a functional interplay of regulatory T cells and CTLA4. Together these mechanisms paralyze CD4 TTS cell differentiation, redirecting metabolic circuits, and reducing their accumulation in the tumor. The paralyzed state is actively maintained throughout cancer progression and CD4 TTS cells rapidly resume proliferation and functional differentiation when the suppressive constraints are alleviated. Overcoming their paralysis established long-term tumor control, demonstrating the importance of rapidly crippling CD4 TTS cells for tumor progression and their potential restoration as therapeutic targets.

ePlatypus: an ecosystem for computational analysis of immunogenomics data

MOTIVATION

The maturation of systems immunology methodologies requires novel and transparent computational frameworks capable of integrating diverse data modalities in a reproducible manner.

RESULTS

Here, we present the ePlatypus computational immunology ecosystem for immunogenomics data analysis, with a focus on adaptive immune repertoires and single-cell sequencing. ePlatypus is an open-source web-based platform and provides programming tutorials and an integrative database that helps elucidate signatures of B and T cell clonal selection. Furthermore, the ecosystem links novel and established bioinformatics pipelines relevant for single-cell immune repertoires and other aspects of computational immunology such as predicting ligand-receptor interactions, structural modeling, simulations, machine learning, graph theory, pseudotime, spatial transcriptomics, and phylogenetics. The ePlatypus ecosystem helps extract deeper insight in computational immunology and immunogenomics and promote open science.

AVAILABILITY AND IMPLEMENTATION

Platypus code used in this manuscript can be found at github.com/alexyermanos/Platypus.

LRT: Integrative analysis of scRNA-seq and scTCR-seq data to investigate clonal differentiation heterogeneity

Single-cell RNA sequencing (scRNA-seq) data has been widely used for cell trajectory inference, with the assumption that cells with similar expression profiles share the same differentiation state. However, the inferred trajectory may not reveal clonal differentiation heterogeneity among T cell clones. Single-cell T cell receptor sequencing (scTCR-seq) data provides invaluable insights into the clonal relationship among cells, yet it lacks functional characteristics. Therefore, scRNA-seq and scTCR-seq data complement each other in improving trajectory inference, where a reliable computational tool is still missing. We developed LRT, a computational framework for the integrative analysis of scTCR-seq and scRNA-seq data to explore clonal differentiation trajectory heterogeneity. Specifically, LRT uses the transcriptomics information from scRNA-seq data to construct overall cell trajectories and then utilizes both the TCR sequence information and phenotype information to identify clonotype clusters with distinct differentiation biasedness. LRT provides a comprehensive analysis workflow, including preprocessing, cell trajectory inference, clonotype clustering, trajectory biasedness evaluation, and clonotype cluster characterization. We illustrated its utility using scRNA-seq and scTCR-seq data of CD8+ T cells and CD4+ T cells with acute lymphocytic choriomeningitis virus infection. These analyses identified several clonotype clusters with distinct skewed distribution along the differentiation path, which cannot be revealed solely based on scRNA-seq data. Clones from different clonotype clusters exhibited diverse expansion capability, V-J gene usage pattern and CDR3 motifs. The LRT framework was implemented as an R package 'LRT', and it is now publicly accessible at https://github.com/JuanXie19/LRT. In addition, it provides two Shiny apps 'shinyClone' and 'shinyClust' that allow users to interactively explore distributions of clonotypes, conduct repertoire analysis, implement clustering of clonotypes, trajectory biasedness evaluation, and clonotype cluster characterization.

COVID-19 relapse associated with SARS-CoV-2 evasion from CD4+ T-cell recognition in an agammaglobulinemia patient

A 25-year-old patient with a primary immunodeficiency lacking immunoglobulin production experienced a relapse after a 239-day period of persistent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Viral genetic sequencing demonstrated that SARS-CoV-2 had evolved during the infection period, with at least five mutations associated with host cellular immune recognition. Among them, the T32I mutation in ORF3a was found to evade recognition by CD4⁺ T cells. The virus found after relapse showed an increased proliferative capacity in vitro. SARS-CoV-2 may have evolved to evade recognition by CD4⁺ T cells and increased in its proliferative capacity during the persistent infection, likely leading to relapse. These mutations may further affect viral clearance in hosts with similar types of human leukocyte antigens. The early elimination of SARS-CoV-2 in immunocompromised patients is therefore important not only to improve the condition of patients but also to prevent the emergence of mutants that threaten public health.

Recruitment of epitope-specific T cell clones with a low-avidity threshold supports efficacy against mutational escape upon re-infection

Repetitive pathogen exposure leads to the dominant outgrowth of T cell clones with high T cell receptor (TCR) affinity to the relevant pathogen-associated antigens. However, low-affinity clones are also known to expand and form immunological memory. While these low-affinity clones contribute less immunity to the original pathogen, their role in protection against pathogens harboring immune escape mutations remains unclear. Based on identification of the TCR repertoire and functionality landscape of naive epitope-specific CD8⁺ T cells, we reconstructed defined repertoires that could be followed as polyclonal populations during immune responses in vivo. We found that selective clonal expansion is governed by clear TCR avidity thresholds. Simultaneously, initial recruitment of broad TCR repertoires provided a polyclonal niche from which flexible secondary responses to mutant epitopes could be recalled. Elucidating how T cell responses develop "from scratch" is informative for the development of enhanced immunotherapies and vaccines.

Themis suppresses the effector function of CD8+ T cells in acute viral infection

CD8+ T cells play a central role in antiviral immune responses. Upon infection, naive CD8+ T cells differentiate into effector cells to eliminate virus-infected cells, and some of these effector cells further differentiate into memory cells to provide long-term protection after infection is resolved. Although extensively investigated, the underlying mechanisms of CD8+ T-cell differentiation remain incompletely understood. Themis is a T-cell-specific protein that plays critical roles in T-cell development. Recent studies using Themis T-cell conditional knockout mice also demonstrated that Themis is required to promote mature CD8+ T-cell homeostasis, cytokine responsiveness, and antibacterial responses. In this study, we used LCMV Armstrong infection as a probe to explore the role of Themis in viral infection. We found that preexisting CD8+ T-cell homeostasis defects and cytokine hyporesponsiveness do not impair viral clearance in Themis T-cell conditional knockout mice. Further analyses showed that in the primary immune response, Themis deficiency promoted the differentiation of CD8+ effector cells and increased their TNF and IFNγ production. Moreover, Themis deficiency impaired memory precursor cell (MPEC) differentiation but promoted short-lived effector cell (SLEC) differentiation. Themis deficiency also enhanced effector cytokine production in memory CD8+ T cells while impairing central memory CD8+ T-cell formation. Mechanistically, we found that Themis mediates PD-1 expression and its signaling in effector CD8+ T cells, which explains the elevated cytokine production in these cells when Themis is disrupted.

Molecular, metabolic and functional CD4 T cell paralysis impedes tumor control

CD4 T cells are important effectors of anti-tumor immunity, yet the regulation of CD4 tumor-specific T (T TS ) cells during cancer development is still unclear. We demonstrate that CD4 T TS cells are initially primed in the tumor draining lymph node and begin to divide following tumor initiation. Distinct from CD8 T TS cells and previously defined exhaustion programs, CD4 T TS cell proliferation is rapidly frozen in place and differentiation stunted by a functional interplay of T regulatory cells and both intrinsic and extrinsic CTLA4 signaling. Together these mechanisms paralyze CD4 T TS cell differentiation, redirecting metabolic and cytokine production circuits, and reducing CD4 T TS cell accumulation in the tumor. Paralysis is actively maintained throughout cancer progression and CD4 T TS cells rapidly resume proliferation and functional differentiation when both suppressive reactions are alleviated. Strikingly, Treg depletion alone reciprocally induced CD4 T TS cells to themselves become tumor-specific Tregs, whereas CTLA4 blockade alone failed to promote T helper differentiation. Overcoming their paralysis established long-term tumor control, demonstrating a novel immune evasion mechanism that specifically cripples CD4 T TS cells to favor tumor progression.

Toward a general model of CD4+ T cell subset specification and memory cell formation

In the past few decades, a number of transformative discoveries have been made regarding memory CD8+ T cell biology; meanwhile, the CD4+ T cell field has lagged behind this progress. This perspective focuses on CD4+ helper T (Th) cell subset specification and memory cell formation. Here, we argue that the sheer number of Th effector and memory cell subsets and a focus on their differences have been a barrier to a general model of CD4+ memory T cell formation that applies to all immune responses. We highlight a bifurcation model that relies on an IL-2 signal-dependent switch as an explanation for the balanced production of diverse Th memory cells that participate in cell-mediated or humoral immunity in most contexts.

Subset-specific and temporal control of effector and memory CD4+ T cell survival

Following their proliferative expansion and differentiation into effector cells like Th1, Tfh, and T central memory precursors (Tcmp), most effector CD4+ T cells die, while some survive and become memory cells. Here, we explored how Bcl-2 family members controlled the survival of CD4+ T cells during distinct phases of mouse acute LCMV infection. During expansion, we found that Th1 cells dominated the response, downregulated expression of Bcl-2, and did not require Bcl-2 for survival. Instead, they relied on the anti-apoptotic protein, A1 for survival. Similarly, Th17 cells in an EAE model also depended on A1 for survival. However, after the peak of the response, CD4+ effector T cells required Bcl-2 to counteract Bim to aid their transition into memory. This Bcl-2 dependence persisted in established memory CD4+ T cells. Combined, these data show a temporal switch in Bcl-2 family-mediated survival of CD4+ T cells over the course of an immune response. This knowledge can help improve T cell survival to boost immunity and conversely, target pathogenic T cells.

Population dynamics and gene regulation of T cells in response to chronic antigen stimulation

T cells are activated by antigen and co-stimulatory receptor signaling and undergo robust proliferation and differentiation into effector cells with protective function. Such quantitatively and qualitatively amplified T cell responses are effective in controlling acute infection and are followed by contraction of the effector population and the formation of resting memory T cells for enhanced protection against previously experienced antigens. However, in the face of persistent antigen during chronic viral infection, in autoimmunity, or in the tumor microenvironment, T cells exhibit distinct responses relative to those in acute insult in several aspects, including reduced clonal expansion and impaired effector function associated with inhibitory receptor expression, a state known as exhaustion. Nevertheless, their responses to chronic infection and tumors are sustained through the establishment of hierarchical heterogeneity, which preserves the duration of the response by generating newly differentiated effector cells. In this review, we highlight recent findings on distinct dynamics of T cell responses under "exhausting" conditions and the roles of the transcription factors that support attenuated yet long-lasting T cell responses as well as the establishment of dysfunctional states.

Clonal lineage tracing reveals mechanisms skewing CD8+ T cell fate decisions in chronic infection

Although recent evidence demonstrates heterogeneity among CD8+ T cells during chronic infection, developmental relationships and mechanisms underlying their fate decisions remain incompletely understood. Using single-cell RNA and TCR sequencing, we traced the clonal expansion and differentiation of CD8+ T cells during chronic LCMV infection. We identified immense clonal and phenotypic diversity, including a subset termed intermediate cells. Trajectory analyses and infection models showed intermediate cells arise from progenitor cells before bifurcating into terminal effector and exhausted subsets. Genetic ablation experiments identified that type I IFN drives exhaustion through an IRF7-dependent mechanism, possibly through an IFN-stimulated subset bridging progenitor and exhausted cells. Conversely, Zeb2 was critical for generating effector cells. Intriguingly, some T cell clones exhibited lineage bias. Mechanistically, we identified that TCR avidity correlates with an exhausted fate, whereas SHP-1 selectively restricts low-avidity effector cell accumulation. Thus, our work elucidates novel mechanisms underlying CD8+ T cell fate determination during persistent infection and suggests two potential pathways leading to exhaustion.

Delineating the transcriptional landscape and clonal diversity of virus-specific CD4+ T cells during chronic viral infection

Although recent evidence indicates that CD4+ T cells responding to chronic viral infection are functionally heterogenous, our understanding of the developmental relationships between these subsets, and a determination of how their transcriptional landscape compares to their acute infection counterparts remains unclear. Additionally, whether cell-intrinsic factors such as TCR usage influence CD4+ T cell fate commitment during persistent infection has not previously been studied. Herein, we perform single-cell RNA sequencing (scRNA-seq) combined with single-cell T cell receptor sequencing (scTCR-seq) on virus-specific CD4+ T cells isolated from mice infected with chronic lymphocytic choriomeningitis virus (LCMV) infection. We identify several transcriptionally distinct states among the Th1, Tfh, and memory-like T cell subsets that form at the peak of infection, including the presence of a previously unrecognized Slamf7+ subset with cytolytic features. We further show that the relative distribution of these populations differs substantially between acute and persistent LCMV infection. Moreover, while the progeny of most T cell clones displays membership within each of these transcriptionally unique populations, overall supporting a one cell-multiple fate model, a small fraction of clones display a biased cell fate decision, suggesting that TCR usage may impact CD4+ T cell development during chronic infection. Importantly, comparative analyses further reveal both subset-specific and core gene expression programs that are differentially regulated between CD4+ T cells responding to acute and chronic LCMV infection. Together, these data may serve as a useful framework and allow for a detailed interrogation into the clonal distribution and transcriptional circuits underlying CD4+ T cell differentiation during chronic viral infection.

Evolution of antigen-specific follicular helper T cell transcription from effector function to memory

Understanding how follicular helper T cells (TFH) regulate the specialization, maturation, and differentiation of adaptive B cell immunity is crucial for developing durable high-affinity immune protection. Using indexed single-cell molecular strategies, we reveal a skewed intraclonal assortment of higher-affinity T cell receptors and the distinct molecular programming of the localized TFH compartment compared with emigrant conventional effector T_H cells. We find a temporal shift in B cell receptor class switch, which permits identification of inflammatory and anti-inflammatory modules of transcriptional programming that subspecialize TFH function before and during the germinal center (GC) reaction. Late collapse of this local primary GC reaction reveals a persistent post-GC TFH population that discloses a putative memory TFH program. These studies define subspecialized antigen-specific TFH transcriptional programs that progressively change with antibody class-specific evolution of high-affinity B cell immunity and a memory TFH transcriptional program that emerges upon local GC resolution.

Autoreactive CD8 T cells in NOD mice exhibit phenotypic heterogeneity but restricted TCR gene usage

Paired scRNA-seq and scTCR-seq reveals that diabetogenic CD8 T cells in the islets and spleens of NOD mice exhibit phenotypic and clonal heterogeneity despite restricted TCR gene usage. Expression of certain TCR genes correlates with clonal proliferation and effector phenotype.

Type 1 diabetes (T1D) is an autoimmune disorder defined by CD8 T cell–mediated destruction of pancreatic β cells. We have previously shown that diabetogenic CD8 T cells in the islets of non-obese diabetic mice are phenotypically heterogeneous, but clonal heterogeneity remains relatively unexplored. Here, we use paired single-cell RNA and T-cell receptor sequencing (scRNA-seq and scTCR-seq) to characterize autoreactive CD8 T cells from the islets and spleens of non-obese diabetic mice. scTCR-seq demonstrates that CD8 T cells targeting the immunodominant β-cell epitope IGRP_206-214 exhibit restricted TCR gene usage. scRNA-seq identifies six clusters of autoreactive CD8 T cells in the islets and six in the spleen, including memory and exhausted cells. Clonal overlap between IGRP_206-214–reactive CD8 T cells in the islets and spleen suggests these cells may circulate between the islets and periphery. Finally, we identify correlations between TCR genes and T-cell clonal expansion and effector fate. Collectively, our work demonstrates that IGRP_206-214–specific CD8 T cells are phenotypically heterogeneous but clonally restricted, raising the possibility of selectively targeting these TCR structures for therapeutic benefit.

Role of cross-reactivity in cellular immune targeting of influenza A M158-66 variant peptide epitopes

The immunologic significance of cross-reactivity of TCR recognition of peptide:MHC complexes is still poorly understood. We have described TCR cross-reactivity in a system involving polyclonal CD8 T cell recognition of the well characterized influenza viral M158-66 epitope. While M158-66 is generally conserved between influenza A isolates, error-prone transcription generates stable variant RNA during infection which could act as novel epitopes. If packaged and viable, variant genomic RNA generates an influenza quasispecies. The stable RNA variants would generate a new transmissible epitope that can select a specific repertoire, which itself should have cross-reactive properties. We tested two candidate peptides in which Thr65 is changed to Ala (A65) or Ser (S65) using recall responses to identify responding T cell clonotypes. Both peptides generated large polyclonal T cell repertoires of their own with repertoire characteristics and cross-reactivity patterns like that observed for the M158-66 repertoire. Both substitutions could be present in viral genomes or mRNA at sufficient frequency during an infection to drive immunity. Peptides from the resulting protein would be a target for CD8 cells irrespective of virus viability or transmissibility. These data support the hypothesis that cross-reactivity is important for immunity against RNA virus infections.

Echidna: integrated simulations of single-cell immune receptor repertoires and transcriptomes

Motivation

Single-cell sequencing now enables the recovery of full-length immune receptor repertoires [B cell receptor (BCR) and T cell receptor (TCR) repertoires], in addition to gene expression information. The feature-rich datasets produced from such experiments require extensive and diverse computational analyses, each of which can significantly influence the downstream immunological interpretations, such as clonal selection and expansion. Simulations produce validated standard datasets, where the underlying generative model can be precisely defined and furthermore perturbed to investigate specific questions of interest. Currently, there is no tool that can be used to simulate single-cell datasets incorporating immune receptor repertoires and gene expression.

Results

We developed Echidna, an R package that simulates immune receptors and transcriptomes at single-cell resolution with user-tunable parameters controlling a wide range of features such as clonal expansion, germline gene usage, somatic hypermutation, transcriptional phenotypes and spatial location. Echidna can additionally simulate time-resolved B cell evolution, producing mutational networks with complex selection histories incorporating class-switching and B cell subtype information. We demonstrated the benchmarking potential of Echidna by simulating clonal lineages and comparing the known simulated networks with those inferred from only the BCR sequences as input. Finally, we simulated immune repertoire information onto existing spatial transcriptomic experiments, thereby generating novel datasets that could be used to develop and integrate methods to profile clonal selection in a spatially resolved manner. Together, Echidna provides a framework that can incorporate experimental data to simulate single-cell immune repertoires to aid software development and bioinformatic benchmarking of clonotyping, phylogenetics, transcriptomics and machine learning strategies.

Availability and implementation

The R package and code used in this manuscript can be found at github.com/alexyermanos/echidna and also in the R package Platypus (Yermanos et al., 2021). Installation instructions and the vignette for Echidna is described in the Platypus Computational Ecosystem (https://alexyermanos.github.io/Platypus/index.html). Publicly available data and corresponding sample accession numbers can be found in Supplementary Tables S2 and S3.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Single-cell immune repertoire sequencing of B and T cells in murine models of infection and autoimmunity

Adaptive immune repertoires are composed by the ensemble of B and T-cell receptors within an individual, reflecting both past and current immune responses. Recent advances in single-cell sequencing enable recovery of the complete adaptive immune receptor sequences in addition to transcriptional information. Here, we recovered transcriptome and immune repertoire information for polyclonal T follicular helper cells following lymphocytic choriomeningitis virus (LCMV) infection, CD8+ T cells with binding specificity restricted to two distinct LCMV peptides, and B and T cells isolated from the nervous system in the context of experimental autoimmune encephalomyelitis. We could relate clonal expansion, germline gene usage, and clonal convergence to cell phenotypes spanning activation, memory, naive, antibody secretion, T-cell inflation, and regulation. Together, this dataset provides a resource for immunologists that can be integrated with future single-cell immune repertoire and transcriptome sequencing datasets.

Divide and Conquer: Phenotypic and Temporal Heterogeneity Within CD8+ T Cell Responses

The advent of technologies that can characterize the phenotypes, functions and fates of individual cells has revealed extensive and often unexpected levels of diversity between cells that are nominally of the same subset. CD8⁺ T cells, also known as cytotoxic T lymphocytes (CTLs), are no exception. Investigations of individual CD8⁺ T cells both in vitro and in vivo have highlighted the heterogeneity of cellular responses at the levels of activation, differentiation and function. This review takes a broad perspective on the topic of heterogeneity, outlining different forms of variation that arise during a CD8⁺ T cell response. Specific attention is paid to the impact of T cell receptor (TCR) stimulation strength on heterogeneity. In particular, this review endeavors to highlight connections between variation at different cellular stages, presenting known mechanisms and key open questions about how variation between cells can arise and propagate.

A CD4+ T cell reference map delineates subtype-specific adaptation during acute and chronic viral infections

CD4+ T cells are critical orchestrators of immune responses against a large variety of pathogens, including viruses. While multiple CD4+ T cell subtypes and their key transcriptional regulators have been identified, there is a lack of consistent definition for CD4+ T cell transcriptional states. In addition, the progressive changes affecting CD4+ T cell subtypes during and after immune responses remain poorly defined. Using single-cell transcriptomics, we characterized the diversity of CD4+ T cells responding to self-resolving and chronic viral infections in mice. We built a comprehensive map of virus-specific CD4+ T cells and their evolution over time, and identified six major cell states consistently observed in acute and chronic infections. During the course of acute infections, T cell composition progressively changed from effector to memory states, with subtype-specific gene modules and kinetics. Conversely, in persistent infections T cells acquired distinct, chronicity-associated programs. By single-cell T cell receptor (TCR) analysis, we characterized the clonal structure of virus-specific CD4+ T cells across individuals. Virus-specific CD4+ T cell responses were essentially private across individuals and most T cells differentiated into both Tfh and Th1 subtypes irrespective of their TCR. Finally, we showed that our CD4+ T cell map can be used as a reference to accurately interpret cell states in external single-cell datasets across tissues and disease models. Overall, this study describes a previously unappreciated level of adaptation of the transcriptional states of CD4+ T cells responding to viruses and provides a new computational resource for CD4+ T cell analysis.

Chemoprevention of Lung Cancer with a Combination of Mitochondria-Targeted Compounds

Simple Summary

Previous reports showed that mitochondria-targeted honokiol and mitochondria-targeted lonidamine potently inhibit complex-I- and complexes-I/II-induced respiration and cancer cell proliferation. In this study, we investigated the efficacy of combining mitochondria-targeted honokiol and mitochondria-targeted lonidamine treatments for lung cancer prevention. We found that their combination exhibited striking tumor inhibition in the benzo[a]pyrene-induced murine lung tumor model without causing detectable side effects. Using single-cell RNA sequencing, we found combined treatment has a clear advantage in that it can significantly inhibit two oncogenic pathways—STAT3 signaling and AKT/mTOR/p70S6K signaling. Such dual inhibition may contribute to the greater efficacy of the combined drug treatment. Therefore, the combination provides a novel option for lung cancer chemoprevention.

Abstract

Combined treatment targeting mitochondria may improve the efficacy of lung cancer chemoprevention. Here, mitochondria-targeted honokiol (Mito-HNK), an inhibitor of mitochondrial complex I and STAT3 phosphorylation, and mitochondria-targeted lonidamine (Mito-LND), an inhibitor of mitochondrial complexes I/II and AKT/mTOR/p70S6K signaling, were evaluated for their combinational chemopreventive efficacy on mouse lung carcinogenesis. All chemopreventive treatments began one-week post-carcinogen treatment and continued daily for 24 weeks. No evidence of toxicity (including liver toxicity) was detected by monitoring serum levels of alanine aminotransferase and aspartate aminotransferase enzymes. Mito-HNK or Mito-LND treatment alone reduced tumor load by 56% and 48%, respectively, whereas the combination of Mito-HNK and Mito-LND reduced tumor load by 83%. To understand the potential mechanism(s) of action for the observed combinatorial effects, single-cell RNA sequencing was performed using mouse tumors treated with Mito-HNK, Mito-LND, and their combination. In lung tumor cells, Mito-HNK treatment blocked the expression of genes involved in mitochondrial complex ǀ, oxidative phosphorylation, glycolysis, and STAT3 signaling. Mito-LND inhibited the expression of genes for mitochondrial complexes I/II, oxidative phosphorylation, and AKT/mTOR/p70S6K signaling in lung tumor cells. In addition to these changes, a combination of Mito-HNK with Mito-LND decreased arginine and proline metabolism, N-glycan biosynthesis, and tryptophan metabolism in lung tumor cells. Our results demonstrate that Mito-LND enhanced the antitumor efficacy of Mito-HNK, where both compounds inhibited common targets (oxidative phosphorylation) as well as unique targets for each agent (STAT3 and mTOR signaling). Therefore, the combination of Mito-HNK with Mito-LND may present an effective strategy for lung cancer chemoprevention.

Clonally Expanded Virus-Specific CD8 T Cells Acquire Diverse Transcriptional Phenotypes During Acute, Chronic, and Latent Infections

CD8+ T cells play a crucial role in the control and resolution of viral infections and can adopt a wide range of phenotypes and effector functions depending on the inflammatory context and the duration and extent of antigen exposure. Similarly, viral infections can exert diverse selective pressures on populations of clonally related T cells. Technical limitations have nevertheless made it challenging to investigate the relationship between clonal selection and transcriptional phenotypes of virus-specific T cells. We therefore performed single-cell T cell receptor (TCR) repertoire and transcriptome sequencing of virus-specific CD8 T cells in murine models of acute, chronic and latent infection. We observed clear infection-specific populations corresponding to memory, effector, exhausted, and inflationary phenotypes. We further uncovered a mouse-specific and polyclonal T cell response, despite all T cells sharing specificity to a single viral epitope, which was accompanied by stereotypic TCR germline gene usage in all three infection types. Persistent antigen exposure during chronic and latent viral infections resulted in a higher proportion of clonally expanded T cells relative to acute infection. We furthermore observed a relationship between transcriptional heterogeneity and clonal expansion for all three infections, with highly expanded clones having distinct transcriptional phenotypes relative to less expanded clones. Together our work relates clonal selection to gene expression in the context of viral infection and further provides a dataset and accompanying software for the immunological community.

Tfh-cell-derived interleukin 21 sustains effector CD8+ T cell responses during chronic viral infection

CD4+ T cell-derived interleukin 21 (IL-21) sustains CD8+ T cell responses during chronic viral infection, but the helper subset that confers this protection remains unclear. Here, we applied scRNA and ATAC-seq approaches to determine the heterogeneity of IL-21+CD4+ T cells during LCMV clone 13 infection. CD4+ T cells were comprised of three transcriptionally and epigenetically distinct populations: Cxcr6+ Th1 cells, Cxcr5+ Tfh cells, and a previously unrecognized Slamf6+ memory-like (Tml) subset. T cell differentiation was specifically redirected toward the Tml subset during chronic, but not acute, LCMV infection. Although this subset displayed an enhanced capacity to accumulate and some developmental plasticity, it remained largely quiescent, which may hinder its helper potential. Conversely, mixed bone marrow chimera experiments revealed that Tfh cell-derived IL-21 was critical to sustain CD8+ T cell responses and viral control. Thus, strategies that bolster IL-21+Tfh cell responses may prove effective in enhancing CD8+ T cell-mediated immunity.

Integrating T cell receptor sequences and transcriptional profiles by clonotype neighbor graph analysis (CoNGA)

Links between T cell clonotypes, as defined by T cell receptor (TCR) sequences, and phenotype, as reflected in gene expression (GEX) profiles, surface protein expression and peptide:major histocompatibility complex binding, can reveal functional relationships beyond the features shared by clonally related cells. Here we present clonotype neighbor graph analysis (CoNGA), a graph theoretic approach that identifies correlations between GEX profile and TCR sequence through statistical analysis of GEX and TCR similarity graphs. Using CoNGA, we uncovered associations between TCR sequence and GEX profiles that include a previously undescribed 'natural lymphocyte' population of human circulating CD8⁺ T cells and a set of TCR sequence determinants of differentiation in thymocytes. These examples show that CoNGA might help elucidate complex relationships between TCR sequence and T cell phenotype in large, heterogeneous, single-cell datasets.

Pre-existing chromatin accessibility and gene expression differences among naive CD4+ T cells influence effector potential

CD4+ T cells have a remarkable potential to differentiate into diverse effector lineages following activation. Here, we probe the heterogeneity present among naive CD4+ T cells before encountering their cognate antigen to ask whether their effector potential is modulated by pre-existing transcriptional and chromatin landscape differences. Single-cell RNA sequencing shows that key drivers of variability are genes involved in T cell receptor (TCR) signaling. Using CD5 expression as a readout of the strength of tonic TCR interactions with self-peptide MHC, and sorting on the ends of this self-reactivity spectrum, we find that pre-existing transcriptional differences among naive CD4+ T cells impact follicular helper T (TFH) cell versus non-TFH effector lineage choice. Moreover, our data implicate TCR signal strength during thymic development in establishing differences in naive CD4+ T cell chromatin landscapes that ultimately shape their effector potential.

Cognate recognition of microbial antigens defines constricted CD4+ T cell receptor repertoires in the inflamed colon

Key aspects of intestinal T cells, including their antigen specificity and their selection by the microbiota and other intestinal antigens, as well as the contribution of individual T cell clones to regulatory and effector functions, remain unresolved. Here we tracked adoptively transferred T cell populations to specify the interrelation of T cell receptor repertoire and the gut antigenic environment. We show that dominant TCRα clonotypes were shared between interferon-γ- and interleukin-17-producing but not regulatory Foxp3⁺ T cells. Identical TCRα clonotypes accumulated in the colon of different individuals, whereas antibiotics or defined colonization correlated with the expansion of distinct expanded T cell clonotypes. Our results demonstrate key aspects of intestinal CD4⁺ T cell activation and suggest that few microbial species exert a dominant effect on the intestinal T cell repertoire during colitis. We speculate that dominant proinflammatory T cell clones might provide a therapeutic target in human inflammatory bowel disease.

Understanding the pathogenesis of infectious diseases by single-cell RNA sequencing

Infections are highly orchestrated and dynamic processes, which involve both pathogen and host. Transcriptional profiling at the single-cell level enables the analysis of cell diversity, heterogeneity of the immune response, and detailed molecular mechanisms underlying infectious diseases caused by bacteria, viruses, fungi, and parasites. Herein, we highlight recent remarkable advances in single-cell RNA sequencing (scRNA-seq) technologies and their applications in the investigation of host-pathogen interactions, current challenges and potential prospects for disease treatment are discussed as well. We propose that with the aid of scRNA-seq, the mechanism of infectious diseases will be further revealed thus inspiring the development of novel interventions and therapies.

Assembly of a spatial circuit of T-bet-expressing T and B lymphocytes is required for antiviral humoral immunity

Effective antiviral immunity requires generation of T and B lymphocytes expressing the transcription factor T-bet, a regulator of type 1 inflammatory responses. Using T-bet expression as an endogenous marker for cells participating in a type 1 response, we report coordinated interactions of T-bet-expressing T and B lymphocytes on the basis of their dynamic colocalization at the T cell zone and B follicle boundary (T-B boundary) and germinal centers (GCs) during lung influenza infection. We demonstrate that the assembly of this circuit takes place in distinct anatomical niches within the draining lymph node, guided by CXCR3 that enables positioning of T_H1 cells at the T-B boundary. The encounter of B and T_H1 cells at the T-B boundary enables IFN-γ produced by the latter to induce IgG2c class switching. Within GCs, T-bet⁺ T_FH cells represent a specialized stable sublineage required for GC growth but dispensable for IgG2c class switching. Our studies show that during respiratory viral infection, T-bet-expressing T and B lymphocytes form a circuit assembled in a spatiotemporally controlled manner that acts as a functional unit enabling a robust and coherent humoral response tailored for optimal antiviral immunity.

Single cell analysis of host response to helminth infection reveals the clonal breadth, heterogeneity, and tissue-specific programming of the responding CD4+ T cell repertoire

The CD4⁺ T cell response is critical to host protection against helminth infection. How this response varies across different hosts and tissues remains an important gap in our understanding. Using IL-4-reporter mice to identify responding CD4⁺ T cells to Nippostrongylus brasiliensis infection, T cell receptor sequencing paired with novel clustering algorithms revealed a broadly reactive and clonally diverse CD4⁺ T cell response. While the most prevalent clones and clonotypes exhibited some tissue selectivity, most were observed to reside in both the lung and lung-draining lymph nodes. Antigen-reactivity of the broader repertoires was predicted to be shared across both tissues and individual mice. Transcriptome, trajectory, and chromatin accessibility analysis of lung and lymph-node repertoires revealed three unique but related populations of responding IL-4⁺ CD4⁺ T cells consistent with T follicular helper, T helper 2, and a transitional population sharing similarity with both populations. The shared antigen reactivity of lymph node and lung repertoires combined with the adoption of tissue-specific gene programs allows for the pairing of cellular and humoral responses critical to the orchestration of anti-helminth immunity.

Using various “omic” approaches, the CD4⁺ T cell receptor (TCR) repertoire was explored after primary helminth infection. Infection generated a broadly reactive and clonally diverse CD4⁺ T cell response with the most prevalent clonotypes and predicted antigen specificities residing in both the lung and lung-draining lymph nodes. Tissue-specific programming of responding CD4⁺ T cells directed the establishment of committed Tfh and Th2 cells, both critical for driving distinct hallmarks of type-2 inflammation. These datasets help to explore the diverse yet tissue-specific nature of anti-helminth immunity.