Immunology in the Era of Single-Cell Technologies.

Immune cells are characterized by diversity, specificity, plasticity, and adaptability-properties that enable them to contribute to homeostasis and respond specifically and dynamically to the many threats encountered by the body. Single-cell technologies, including the assessment of transcriptomics, genomics, and proteomics at the level of individual cells, are ideally suited to studying these properties of immune cells. In this review we discuss the benefits of adopting single-cell approaches in studying underappreciated qualities of immune cells and highlight examples where these technologies have been critical to advancing our understanding of the immune system in health and disease.

TCR+/BCR+ dual-expressing cells and their associated public BCR clonotype are not enriched in type 1 diabetes.

Ahmed and colleagues recently described a novel hybrid lymphocyte expressing both a B and T cell receptor, termed double expresser (DE) cells. DE cells in blood of type 1 diabetes (T1D) subjects were present at increased numbers and enriched for a public B cell clonotype. Here, we attempted to reproduce these findings. While we could identify DE cells by flow cytometry, we found no association between DE cell frequency and T1D status. We were unable to identify the reported public B cell clone, or any similar clone, in bulk B cells or sorted DE cells from T1D subjects or controls. We also did not observe increased usage of the public clone VH or DH genes in B cells or in sorted DE cells. Taken together, our findings suggest that DE cells and their alleged public clonotype are not enriched in T1D. This Matters Arising paper is in response to Ahmed et al. (2019), published in Cell. See also the response by Ahmed et al. (2021), published in this issue.

Early developing B cells undergo negative selection by central nervous system-specific antigens in the meninges.

Self-reactive B cell progenitors are eliminated through central tolerance checkpoints, a process thought to be restricted to the bone marrow in mammals. Here, we identified a consecutive trajectory of B cell development in the meninges of mice and non-human primates. The meningeal B cells were located predominantly at the dural sinuses, where endothelial cells expressed essential niche factors to support B cell development. Parabiosis experiments together with lineage tracing showed that meningeal developing B cells were replenished continuously from hematopoietic stem cell (HSC)-derived progenitors via a circulation-independent route. Autoreactive immature B cells that recognized myelin oligodendrocyte glycoprotein (MOG), a central nervous system-specific antigen, were eliminated specifically from the meninges. Furthermore, genetic deletion of the Mog gene restored the self-reactive B cell population in the meninges. These findings identify the meninges as a distinct reservoir for B cell development, allowing in situ negative selection to ensure a locally non-self-reactive immune repertoire.

The Transcription Factor T-bet Resolves Memory B Cell Subsets with Distinct Tissue Distributions and Antibody Specificities in Mice and Humans.

B cell subsets expressing the transcription factor T-bet are associated with humoral immune responses and autoimmunity. Here, we examined the anatomic distribution, clonal relationships, and functional properties of T-bet+ and T-bet- memory B cells (MBCs) in the context of the influenza-specific immune response. In mice, both T-bet- and T-bet+ hemagglutinin (HA)-specific B cells arose in germinal centers, acquired memory B cell markers, and persisted indefinitely. Lineage tracing and IgH repertoire analyses revealed minimal interconversion between T-bet- and T-bet+ MBCs, and parabionts showed differential tissue residency and recirculation properties. T-bet+ MBCs could be subdivided into recirculating T-betlo MBCs and spleen-resident T-bethi MBCs. Human MBCs displayed similar features. Conditional gene deletion studies revealed that T-bet expression in B cells was required for nearly all HA stalk-specific IgG2c antibodies and for durable neutralizing titers to influenza. Thus, T-bet expression distinguishes MBC subsets that have profoundly different homing, residency, and functional properties, and mediate distinct aspects of humoral immune memory.

Genotype-phenotype landscapes for immune-pathogen coevolution.

Our immune systems constantly coevolve with the pathogens that challenge them, as pathogens adapt to evade our defense responses, with our immune repertoires shifting in turn. These coevolutionary dynamics take place across a vast and high-dimensional landscape of potential pathogen and immune receptor sequence variants. Mapping the relationship between these genotypes and the phenotypes that determine immune-pathogen interactions is crucial for understanding, predicting, and controlling disease. Here, we review recent developments applying high-throughput methods to create large libraries of immune receptor and pathogen protein sequence variants and measure relevant phenotypes. We describe several approaches that probe different regions of the high-dimensional sequence space and comment on how combinations of these methods may offer novel insight into immune-pathogen coevolution.

Inferring the T cell repertoire dynamics of healthy individuals.

The adaptive immune system is a diverse ecosystem that responds to pathogens by selecting cells with specific receptors. While clonal expansion in response to particular immune challenges has been extensively studied, we do not know the neutral dynamics that drive the immune system in the absence of strong stimuli. Here, we learn the parameters that underlie the clonal dynamics of the T cell repertoire in healthy individuals of different ages, by applying Bayesian inference to longitudinal immune repertoire sequencing (RepSeq) data. Quantifying the experimental noise accurately for a given RepSeq technique allows us to disentangle real changes in clonal frequencies from noise. We find that the data are consistent with clone sizes following a geometric Brownian motion and show that its predicted steady state is in quantitative agreement with the observed power-law behavior of the clone-size distribution. The inferred turnover time scale of the repertoire increases with patient age and depends on the clone size in some individuals.

IMperm: a fast and comprehensive IMmune Paired-End Reads Merger for sequencing data.

The adaptive immune receptor repertoire (AIRR), consisting of T- and B-cell receptors, is the core component of the immune system. The AIRR sequencing is commonly used in cancer immunotherapy and minimal residual disease (MRD) detection of leukemia and lymphoma. The AIRR is captured by primers and sequenced to yield paired-end (PE) reads. The PE reads could be merged into one sequence by the overlapped region between them. However, the wide range of AIRR data raises the difficulty, so a special tool is required. We developed a software package for IMmune PE reads merger of sequencing data, named IMperm. We used the k-mer-and-vote strategy to pin down the overlapped region rapidly. IMperm could handle all types of PE reads, eliminate adapter contamination and successfully merge low-quality and minor/non-overlapping reads. Compared with existing tools, IMperm performed better in both simulated and sequencing data. Notably, IMperm was well suited to processing the data of MRD detection in leukemia and lymphoma and detected 19 novel MRD clones in 14 patients with leukemia from previously published data. Additionally, IMperm can handle PE reads from other sources, and we demonstrated its effectiveness on two genomic and one cell-free deoxyribonucleic acid datasets. IMperm is implemented in the C programming language and consumes little runtime and memory. It is freely available at https://github.com/zhangwei2015/IMperm.

Anchor Clustering for million-scale immune repertoire sequencing data.

BACKGROUND: The clustering of immune repertoire data is challenging due to the computational cost associated with a very large number of pairwise sequence comparisons. To overcome this limitation, we developed Anchor Clustering, an unsupervised clustering method designed to identify similar sequences from millions of antigen receptor gene sequences. First, a Point Packing algorithm is used to identify a set of maximally spaced anchor sequences. Then, the genetic distance of the remaining sequences to all anchor sequences is calculated and transformed into distance vectors. Finally, distance vectors are clustered using unsupervised clustering. This process is repeated iteratively until the resulting clusters are small enough so that pairwise distance comparisons can be performed. RESULTS: Our results demonstrate that Anchor Clustering is faster than existing pairwise comparison clustering methods while providing similar clustering quality. With its flexible, memory-saving strategy, Anchor Clustering is capable of clustering millions of antigen receptor gene sequences in just a few minutes. CONCLUSIONS: This method enables the meta-analysis of immune-repertoire data from different studies and could contribute to a more comprehensive understanding of the immune repertoire data space.

Computational estimation of clonal diversity in autoimmunity.

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

Identification of the immune repertoire of γδ T-cell receptors in psoriasis.

The γδ T cells are a subpopulation of T cells that are abundantly found in the skin and mucous membranes. Their reactivity to self-antigens and ability to secrete various cytokines make them a key component in psoriasis development. Although the correlation between the immune repertoire (IR) of γδ T-cell receptors and the occurrence and severity of psoriasis remains incompletely explored, high-throughput sequencing of γδ T cells has led to a deeper understanding of IR in psoriasis. This study investigated the differences between γδ T cells in patients with psoriasis and healthy controls. The γδ T cells were identified via immunofluorescence staining and a correlation analysis was performed according to the psoriasis area and severity index (PASI) scores. The IR sequencing method was used to detect IR in the γδ T-cell receptors. The findings demonstrated more skin γδ T cells in patients with psoriasis, which were positively correlated with the PASI score. There were subtle differences in most variable (V), diversity (D) and joining (J) gene segments and VJ/VDJ combination segments between patients with psoriasis and healthy controls. However, a higher diversity of complementarity-determining region 3 (CDR3) was observed in patients with psoriasis. In summary, the IR of skin γδ T cells was significantly altered in patients with psoriasis, and the diversity in the cell's CDR3 population is a promising biomarker for assessment of psoriasis severity.

TCRß rearrangements without a D segment are common, abundant, and public.

T cells play an important role in adaptive immunity. An enormous clonal diversity of T cells with a different specificity, encoded by the T cell receptor (TCR), protect the body against infection. Most TCRß chains are generated from a V, D, and J segment during recombination in the thymus. Although complete absence of the D segment is not easily detectable from sequencing data, we find convincing evidence for a substantial proportion of TCRß rearrangements lacking a D segment. Additionally, sequences without a D segment are more likely to be abundant within individuals and/or shared between individuals. Our analysis indicates that such sequences are preferentially generated during fetal development and persist within the elderly. Summarizing, TCRß rearrangements without a D segment are not uncommon, and tend to allow for TCRß chains with a high abundance in the naive repertoire.

Nucleic acid-triggered tumoral immunity propagates pH-selective therapeutic antibodies through tumor-driven epitope spreading.

Important roles of humoral tumor immunity are often pointed out; however, precise profiles of dominant antigens and developmental mechanisms remain elusive. We systematically investigated the humoral antigens of dominant intratumor immunoglobulin clones found in human cancers. We found that approximately half of the corresponding antigens were restricted to strongly and densely negatively charged polymers, resulting in simultaneous reactivities of the antibodies to both densely sulfated glycosaminoglycans (dsGAGs) and nucleic acids (NAs). These anti-dsGAG/NA antibodies matured and expanded via intratumoral immunological driving force of innate immunity via NAs. These human cancer-derived antibodies exhibited acidic pH-selective affinity across both antigens and showed specific reactivity to diverse spectrums of human tumor cells. The antibody-drug conjugate exerted therapeutic effects against multiple cancers in vivo by targeting cell surface dsGAG antigens. This study reveals that intratumoral immunological reactions propagate tumor-oriented immunoglobulin clones and demonstrates a new therapeutic modality for the universal treatment of human malignancies.

Systemic lupus erythematosus patients contain B-cell receptor repertoires sensitive to immunosuppressive drugs.

Immune repertoire (IR) during treatment may be a surrogate biomarker for disease response. Changes of the IR in systemic lupus erythematosus patients in response to immunosuppressive drugs were identified in ten SLE patients. Patients provided peripheral blood mononuclear cells at two time points for sequencing. They were divided into sensitive and nonsensitive groups by their clinical responses to immunosuppressive drugs. After treatment, the BCR expression significantly decreased in patients from the sensitive group while there was no change in patients from the nonsensitive group. IgM comprised a dominant portion of the BCR repertoire and increased slightly in all patients in the sensitive group but decreased in the nonsensitive group. IgA also exhibited opposing changes between the two groups. Shorter CDR3 of TRB and TRG chains occurred in the sensitive group. CDR3 length of IGK decreased significantly in the sensitive group. CDR3 of TCR δ/γ changed distinctly between time points in the sensitive group. Six immune-related genes showed differential expression levels in sensitive and nonsensitive groups. Our study shows that it is BCR repertoire sensitivity to immunosuppressive drugs in SLE patients and sheds light on personalized therapy for SLE.

Immune repertoire sequencing reveals an abnormal adaptive immune system in COVID-19 survivors.

Accumulating evidence suggests that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) impairs the adaptive immune system during acute infection. Still, it remains largely unclear whether the frequency and functions of T and B cells return to normal after the recovery of Coronavirus Disease 2019 (COVID-19). Here, we analyzed immune repertoires and SARS-CoV-2-specific neutralization antibodies in a prospective cohort of 40 COVID-19 survivors with a 6-month follow-up after hospital discharge. Immune repertoire sequencing revealed abnormal T- and B-cell expression and function with large T cell receptor/B cell receptor clones, decreased diversity, abnormal class-switch recombination, and somatic hypermutation. A decreased number of B cells but an increased proportion of CD19+ CD138+ B cells were found in COVID-19 survivors. The proportion of CD4+ T cells, especially circulating follicular helper T (cTfh) cells, was increased, whereas the frequency of CD3+ CD4- T cells was decreased. SARS-CoV-2-specific neutralization IgG and IgM antibodies were identified in all survivors, especially those recorded with severe COVID-19 who showed a higher inhibition rate of neutralization antibodies. All severe cases complained of more than one COVID-19 sequelae after 6 months of recovery. Overall, our findings indicate that SARS-CoV-2-specific antibodies remain detectable even after 6 months of recovery. Because of their abnormal adaptive immune system with a low number of CD3+ CD4- T cells and high susceptibility to infections, COVID-19 patients might need more time and medical care to fully recover from immune abnormalities and tissue damage.

Immune repertoire profiling for disease pathobiology.

Lymphocytes consist of highly heterogeneous populations, each expressing a specific cell surface receptor corresponding to a particular antigen. Lymphocytes are both the cause and regulator of various diseases, including autoimmune/allergic diseases, lifestyle diseases, neurodegenerative diseases, and cancers. Recently, immune repertoire sequencing has attracted much attention because it helps obtain global profiles of the immune receptor sequences of infiltrating T and B cells in specimens. Immune repertoire sequencing not only helps deepen our understanding of the molecular mechanisms of immune-related pathology but also assists in discovering novel therapeutic modalities for diseases, thereby shedding colorful light on otherwise tiny monotonous cells when observed under a microscope. In this review article, we introduce and detail the background and methodology of immune repertoire sequencing and summarize recent scientific achievements in association with human diseases. Future perspectives on this genetic technique in the field of histopathological research will also be discussed.

A pathogenic and clonally expanded B cell transcriptome in active multiple sclerosis.

Central nervous system B cells have several potential roles in multiple sclerosis (MS): secretors of proinflammatory cytokines and chemokines, presenters of autoantigens to T cells, producers of pathogenic antibodies, and reservoirs for viruses that trigger demyelination. To interrogate these roles, single-cell RNA sequencing (scRNA-Seq) was performed on paired cerebrospinal fluid (CSF) and blood from subjects with relapsing-remitting MS (RRMS; n = 12), other neurologic diseases (ONDs; n = 1), and healthy controls (HCs; n = 3). Single-cell immunoglobulin sequencing (scIg-Seq) was performed on a subset of these subjects and additional RRMS (n = 4), clinically isolated syndrome (n = 2), and OND (n = 2) subjects. Further, paired CSF and blood B cell subsets (RRMS; n = 7) were isolated using fluorescence activated cell sorting for bulk RNA sequencing (RNA-Seq). Independent analyses across technologies demonstrated that nuclear factor kappa B (NF-κB) and cholesterol biosynthesis pathways were activated, and specific cytokine and chemokine receptors were up-regulated in CSF memory B cells. Further, SMAD/TGF-ß1 signaling was down-regulated in CSF plasmablasts/plasma cells. Clonally expanded, somatically hypermutated IgM+ and IgG1+ CSF B cells were associated with inflammation, blood-brain barrier breakdown, and intrathecal Ig synthesis. While we identified memory B cells and plasmablast/plasma cells with highly similar Ig heavy-chain sequences across MS subjects, similarities were also identified with ONDs and HCs. No viral transcripts, including from Epstein-Barr virus, were detected. Our findings support the hypothesis that in MS, CSF B cells are driven to an inflammatory and clonally expanded memory and plasmablast/plasma cell phenotype.

Comprehensive analysis of B and T cell receptor repertoire in patients after kidney transplantation by high-throughput sequencing.

PURPOSE: The dynamic immunity of kidney transplant patients has not been fully elucidated. In this study, we explored the repertoire features of B/T cell receptor (BCR/TCR) of kidney transplant patients. METHODS: Using combined multiplex PCR amplification and high-throughput sequencing technique, we analyzed the uremic patients' BCR H chain and TCR beta chain repertoire which obtained 1 day before kidney transplantation (PRE-1), 1 day and 7 day after kidney transplantation (POST-1 and POST-7). RESULTS: Our analysis results showed the diversity of TCRß CDR3 in POST-7 group was highest. In addition, there were specific skewed usage of TRBV gene subfamilies, and V-J combinations in different time points during kidney transplantation. Moreover, the overlap degrees of BCR-H (TCR-ß) CDR3 repertoire among each group were identified. Notably, the abundance of some TCR-ß CDR3 sequences changed regularly in the time point of kidney transplantation. CONCLUSIONS: The BCR-H (TCR-ß) CDR3 repertoire of kidney transplant patients changed dynamically.

Analysis of T cell receptor clonotypes in tumor microenvironment identifies shared cancer-type-specific signatures.

Despite the conventional view that a truly random V(D)J recombination process should generate a highly diverse immune repertoire, emerging reports suggest that there is a certain bias toward the generation of shared/public immune receptor chains. These studies were performed in viral diseases where public T cell receptors (TCR) appear to confer better protective responses. Selective pressures generating common TCR clonotypes are currently not well understood, but it is believed that they confer a growth advantage. As very little is known about public TCR clonotypes in cancer, here we set out to determine the extent of shared TCR clonotypes in the intra-tumor microenvironments of virus- and non-virus-driven head and neck cancers using TCR sequencing. We report that tumor-infiltrating T cell clonotypes were indeed shared across individuals with the same cancer type, where the majority of shared sequences were specific to the cancer type (i.e., viral versus non-viral). These shared clonotypes were not particularly enriched in EBV-associated nasopharynx cancer but, in both cancers, exhibited distinct characteristics, namely shorter CDR3 lengths, restricted V- and J-gene usages, and also demonstrated convergent V(D)J recombination. Many of these shared TCRs were expressed in patients with a shared HLA background. Pattern recognition of CDR3 amino acid sequences revealed strong convergence to specific pattern motifs, and these motifs were uniquely found to each cancer type. This suggests that they may be enriched for specificity to common antigens found in the tumor microenvironment of different cancers. The identification of shared TCRs in infiltrating tumor T cells not only adds to our understanding of the tumor-adaptive immune recognition but could also serve as disease-specific biomarkers and guide the development of future immunotherapies.

The intrarenal landscape of T cell receptor repertoire in clear cell renal cell cancer.

BACKGROUND: Clear cell renal cell cancer (ccRCC) is accompanied by T-cell infiltration. In this study, we sought to determine the difference in T-cell infiltration and the T-cell receptor (TCR) immune repertoire between ccRCC and peritumour tissue. METHODS: T-cell infiltration was examined using immunohistochemistry (IHC) and haematoxylin and eosin (HE) staining. The chi-squared test and Pearson correlation analysis were applied to evaluate the relationship between clinical traits and CD3, CD4, and CD8 expression. Immune repertoire sequencing (IR-Seq) was used to describe the profile of the TCR repertoire. RESULTS: The adjacent tissue showed increased expression of CD3, CD4 and CD8 compared with ccRCC tissue (PCD3 = 0.033; PCD4 = 0.014; PCD8 = 0.004). Indicated CD3+ T-cell density in ccRCC tissue was positively correlated with that in peritumour tissue (P = 0.010, r = 0.514), which implied the T cells in peritumour tissue directly infect the number of cells infiltrating in ccRCC tissue. Moreover, there was a positive correlation between Vimentin expression and indicated positive T-cell marker in ccRCC tissue (PCD3 = 0.035; PCD4 = 0.020; PCD8 = 0.027). Advanced stage revealed less CD4+ T-cell infiltration in ccRCC tissue (PCD4 = 0.023). The results from IR-Seq revealed an obvious increase in VJ and VDJ segment usage, as well as higher complementarity-determining region 3 (CDR3) amino acid (aa) clonotypes in ccRCC. The matched antigen recognized by the TCR of ccRCC may be potential targets. CONCLUSIONS: The current study collectively demonstrates diminished T-cell infiltration and increased CDR3 aa diversity in ccRCC, which may be associated with immunotherapeutic targets for ccRCC patients.

Heterogeneous origin of IgE in atopic dermatitis and psoriasis revealed by B cell receptor repertoire analysis.

BACKGROUND: Epicutaneous sensitization is an important route for the production of IgE, and skin inflammation-induced IgE has recently been reported having features of natural antibody. Atopic dermatitis (AD) and psoriasis have differentially increased level of serum IgE; however, the production mechanism of IgE in these inflammatory skin diseases remains unknown. OBJECTIVE: To explore the origin of IgE in AD and psoriasis by analyzing the B cell receptor repertoire. METHODS: mRNA was prepared from peripheral blood mononuclear cells of AD and psoriasis patients that had elevated serum levels of IgE, and immunoglobulin heavy chain (IGH) repertoires were sequenced after reverse transcription. Clonal lineages of B cells containing members expressing IgE were identified, and somatic hypermutations in IGH inherited from common ancestors within the clonal lineage were used to infer the relationships between B cells. RESULTS: The proportions of IGHE from AD and psoriasis were higher than that of normal control, which were positively correlated with the levels of serum total IgE. The somatic hypermutation value of IGHE variable region was lower than that of IGHG and IGHA, but higher than IGHM and IGHD, indicating a mixed natural and adaptive origins of IgE; and psoriasis demonstrated lower level of hypermutation than AD. The Shannon indexes of CDR3 in IGHE of AD and psoriasis were higher than that of normal control, also supporting the natural origin. The VH usage of IgE was weakly biased in AD and psoriasis patients with high level of house dust mite-specific IgE. Comparison of the number of shared mutations in multi-isotype lineages containing IgE showed that isotype-switching from IgG-expressing B cells might be the major source of IgE in AD and psoriasis. CONCLUSION: IgE has heterogeneous origin in AD and psoriasis, and skin inflammation may contribute to the increased production of natural IgE.