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.

Asymmetric cell division safeguards memory CD8 T cell development.

The strength of T cell receptor (TCR) stimulation and asymmetric distribution of fate determinants are both implied to affect T cell differentiation. Here, we uncover asymmetric cell division (ACD) as a safeguard mechanism for memory CD8 T cell generation specifically upon strong TCR stimulation. Using live imaging approaches, we find that strong TCR stimulation induces elevated ACD rates, and subsequent single-cell-derived colonies comprise both effector and memory precursor cells. The abundance of memory precursor cells emerging from a single activated T cell positively correlates with first mitosis ACD. Accordingly, preventing ACD by inhibition of protein kinase Cζ (PKCζ) during the first mitosis upon strong TCR stimulation markedly curtails the formation of memory precursor cells. Conversely, no effect of ACD on fate commitment is observed upon weak TCR stimulation. Our data provide relevant mechanistic insights into the role of ACD for CD8 T cell fate regulation upon different activation conditions.

Generation of a single-cell B cell atlas of antibody repertoires and transcriptomes to identify signatures associated with antigen specificity.

Although new genomics-based pipelines have potential to augment antibody discovery, these methods remain in their infancy due to an incomplete understanding of the selection process that governs B cell clonal selection, expansion, and antigen specificity. Furthermore, it remains unknown how factors such as aging and reduction of tolerance influence B cell selection. Here we perform single-cell sequencing of antibody repertoires and transcriptomes of murine B cells following immunizations with a model therapeutic antigen target. We determine the relationship between antibody repertoires, gene expression signatures, and antigen specificity across 100,000 B cells. Recombinant expression and characterization of 227 monoclonal antibodies revealed the existence of clonally expanded and class-switched antigen-specific B cells that were more frequent in young mice. Although integrating multiple repertoire features such as germline gene usage and transcriptional signatures failed to distinguish antigen-specific from nonspecific B cells, other features such as immunoglobulin G (IgG) subtype and sequence composition correlated with antigen specificity.

Persistent virus-specific and clonally expanded antibody-secreting cells respond to induced self-antigen in the CNS.

B cells contribute to the pathogenesis of both cellular- and humoral-mediated central nervous system (CNS) inflammatory diseases through a variety of mechanisms. In such conditions, B cells may enter the CNS parenchyma and contribute to local tissue destruction. It remains unexplored, however, how infection and autoimmunity drive transcriptional phenotypes, repertoire features, and antibody functionality. Here, we profiled B cells from the CNS of murine models of intracranial (i.c.) viral infections and autoimmunity. We identified a population of clonally expanded, antibody-secreting cells (ASCs) that had undergone class-switch recombination and extensive somatic hypermutation following i.c. infection with attenuated lymphocytic choriomeningitis virus (rLCMV). Recombinant expression and characterisation of these antibodies revealed specificity to viral antigens (LCMV glycoprotein GP), correlating with ASC persistence in the brain weeks after resolved infection. Furthermore, these virus-specific ASCs upregulated proliferation and expansion programs in response to the conditional and transient induction of the LCMV GP as a neo-self antigen by astrocytes. This class-switched, clonally expanded, and mutated population persisted and was even more pronounced when peripheral B cells were depleted prior to autoantigen induction in the CNS. In contrast, the most expanded B cell clones in mice with persistent expression of LCMV GP in the CNS did not exhibit neo-self antigen specificity, potentially a consequence of local tolerance induction. Finally, a comparable population of clonally expanded, class-switched, and proliferating ASCs was detected in the cerebrospinal fluid of relapsing multiple sclerosis (RMS) patients. Taken together, our findings support the existence of B cells that populate the CNS and are capable of responding to locally encountered autoantigens.

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.

Tissue-resident CD8+ T cells drive compartmentalized and chronic autoimmune damage against CNS neurons.

The mechanisms underlying the chronicity of autoimmune diseases of the central nervous system (CNS) are largely unknown. In particular, it is unclear whether tissue-resident memory T cells (TRM) contribute to lesion pathogenesis during chronic CNS autoimmunity. Here, we observed that a high frequency of brain-infiltrating CD8+ T cells exhibit a TRM-like phenotype in human autoimmune encephalitis. Using mouse models of neuronal autoimmunity and a combination of T single-cell transcriptomics, high-dimensional flow cytometry, and histopathology, we found that pathogenic CD8+ T cells behind the blood-brain barrier adopt a characteristic TRM differentiation program, and we revealed their phenotypic and functional heterogeneity. In the diseased CNS, autoreactive tissue-resident CD8+ T cells sustained focal neuroinflammation and progressive loss of neurons, independently of recirculating CD8+ T cells. Consistently, a large fraction of autoreactive tissue-resident CD8+ T cells exhibited proliferative potential as well as proinflammatory and cytotoxic properties. Persistence of tissue-resident CD8+ T cells in the CNS and their functional output, but not their initial differentiation, were crucially dependent on CD4+ T cells. Collectively, our results point to tissue-resident CD8+ T cells as essential drivers of chronic CNS autoimmunity and suggest that therapies targeting this compartmentalized autoreactive T cell subset might be effective for treating CNS autoimmune diseases.

Tissue-resident memory CD8+ T cells cooperate with CD4+ T cells to drive compartmentalized immunopathology in the CNS.

In chronic inflammatory diseases of the central nervous system (CNS), immune cells persisting behind the blood-brain barrier are supposed to promulgate local tissue destruction. The drivers of such compartmentalized inflammation remain unclear, but tissue-resident memory T cells (TRM) represent a potentially important cellular player in this process. Here, we investigated whether resting CD8+ TRM persisting after cleared infection with attenuated lymphocytic choriomeningitis virus (LCMV) can initiate immune responses directed against cognate self-antigen in the CNS. We demonstrated that time-delayed conditional expression of the LCMV glycoprotein as neo-self-antigen by glia cells reactivated CD8+ TRM. Subsequently, CD8+ TRM expanded and initiated CNS inflammation and immunopathology in an organ-autonomous manner independently of circulating CD8+ T cells. However, in the absence of CD4+ T cells, TCF-1+ CD8+ TRM failed to expand and differentiate into terminal effectors. Similarly, in human demyelinating CNS autoimmune lesions, we found CD8+ T cells expressing TCF-1 that predominantly exhibited a TRM-like phenotype. Together, our study provides evidence for CD8+ TRM-driven CNS immunopathology and sheds light on why inflammatory processes may evade current immunomodulatory treatments in chronic autoimmune CNS conditions.

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.

Profiling the specificity of clonally expanded plasma cells during chronic viral infection by single-cell analysis.

Plasma cells and their secreted antibodies play a central role in the long-term protection against chronic viral infection. However, due to experimental limitations, a comprehensive description of linked genotypic, phenotypic, and antibody repertoire features of plasma cells (gene expression, clonal frequency, virus specificity, and affinity) has been challenging to obtain. To address this, we performed single-cell transcriptome and antibody repertoire sequencing of the murine BM plasma cell population following chronic lymphocytic choriomeningitis virus infection. Our single-cell sequencing approach recovered full-length and paired heavy- and light-chain sequence information for thousands of plasma cells and enabled us to perform recombinant antibody expression and specificity screening. Antibody repertoire analysis revealed that, relative to protein immunization, chronic infection led to increased levels of clonal expansion, class-switching, and somatic variants. Furthermore, antibodies from the highly expanded and class-switched (IgG) plasma cells were found to be specific for multiple viral antigens and a subset of clones exhibited cross-reactivity to nonviral and autoantigens. Integrating single-cell transcriptome data with antibody specificity suggested that plasma cell transcriptional phenotype was correlated to viral antigen specificity. Our findings demonstrate that chronic viral infection can induce and sustain plasma cell clonal expansion, combined with significant somatic hypermutation, and can generate cross-reactive antibodies.

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.