A Public BCR Present in a Unique Dual-Receptor-Expressing Lymphocyte from Type 1 Diabetes Patients Encodes a Potent T Cell Autoantigen

Decoding B Cells in Autoimmune Diseases Through ScRNA + BCR-Seq: Current Knowledge and Future Directions

The combined application of single-cell RNA sequencing (scRNA-seq) and single-cell B-cell receptor sequencing (scBCR-seq) offers a multidimensional perspective for dissecting the immunopathological mechanisms of B cells in autoimmune diseases. This review systematically summarizes the principles of these techniques, the analytical framework, and their key applications in diseases such as systemic lupus erythematosus et. al. It reveals the dynamic correlations between the transcriptome of B-cell subsets and B-cell receptor (BCR) clones. Furthermore, we focus on the potential roles of dual BCR B cells and B/T biphenotypic cells in autoimmunity, emphasizing their exacerbation of disease progression through abnormal clonal expansion and autoantibody secretion. By sorting through cutting-edge advancements and bottleneck issues, this article aims to propel the innovation of multi-omics research and precision treatment paradigms for autoimmune diseases.

Deciphering Autoimmune Diseases: Unveiling the Diagnostic, Therapeutic, and Prognostic Potential of Immune Repertoire Sequencing

Autoimmune diseases (AIDs) are immune system disorders where the body exhibits an immune response to its own antigens, causing damage to its own tissues and organs. The pathogenesis of AIDs is incompletely understood. However, recent advances in immune repertoire sequencing (IR-seq) technology have opened-up a new avenue to study the IR. These studies have revealed the prevalence in IR alterations, potentially inducing AIDs by disrupting immune tolerance and thereby contributing to our comprehension of AIDs. IR-seq harbors significant potential for the clinical diagnosis, personalized treatment, and prognosis of AIDs. This article reviews the application and progress of IR-seq in diseases, such as multiple sclerosis, systemic lupus erythematosus, rheumatoid arthritis, and type 1 diabetes, to enhance our understanding of the pathogenesis of AIDs and offer valuable references for the diagnosis and treatment of AIDs.

Mutagenesis-based optimal design of plant peptide phytosulfokine for enhanced biological activity

Recognition of phytosulfokine (PSK), a sulfated pentapeptide, by its receptor PSKRs is crucial in regulating plant growth, development, and reproduction. However, designing highly active PSK remains a formidable challenge due to the lack of understanding of the structure-property relationship, structural dynamics, and the binding characteristics of PSK. Here, with a combined theoretical and experimental approach, we have investigated the binding dynamics of key interactions between PSK and AtPSKR1LRR to reveal the molecular mechanism of PSK recognition. Our molecular dynamics simulations and free energy perturbation calculations demonstrate that the sulfated tyrosines (PSKsY1 and PSKsY3) are indispensable for forming stable PSK-AtPSKR1LRR complex, while the alanine substitution at PSKQ5 site is rather tolerated. Furthermore, two promising PSK peptide analogs (PSKQ5A and PSKQ5K) with enhanced biological activity have been designed through in silico mutagenesis studies and in vivo experiments. They have a strong promoting effect (20 % enhancement) on stimulating root development compared with the wild-type PSK treatment. This work offers an effective strategy to design new peptide-based drugs for facilitating plant growth and consequent crop productivity, potentially benefiting efforts to address the global food crisis.

The landscape of cell lineage tracing

Cell fate changes play a crucial role in the processes of natural development, disease progression, and the efficacy of therapeutic interventions. The definition of the various types of cell fate changes, including cell expansion, differentiation, transdifferentiation, dedifferentiation, reprogramming, and state transitions, represents a complex and evolving field of research known as cell lineage tracing. This review will systematically introduce the research history and progress in this field, which can be broadly divided into two parts: prospective tracing and retrospective tracing. The initial section encompasses an array of methodologies pertaining to isotope labeling, transient fluorescent tracers, non-fluorescent transient tracers, non-fluorescent genetic markers, fluorescent protein, genetic marker delivery, genetic recombination, exogenous DNA barcodes, CRISPR-Cas9 mediated DNA barcodes, and base editor-mediated DNA barcodes. The second part of the review covers genetic mosaicism, genomic DNA alteration, TCR/BCR, DNA methylation, and mitochondrial DNA mutation. In the final section, we will address the principal challenges and prospective avenues of enquiry in the field of cell lineage tracing, with a particular focus on the sequencing techniques and mathematical models pertinent to single-cell genetic lineage tracing, and the value of pursuing a more comprehensive investigation at both the spatial and temporal levels in the study of cell lineage tracing.

Renal Tubular Epithelial Cell-Derived IL-7 Drives Expansion and Protective Effects of Double-Negative T Cells in Acute Kidney Injury

BACKGROUND

Previous studies have demonstrated that double-negative (DN) αβ T cells play an important role in immune responses during ischemic acute kidney injury (AKI). Here, we investigate the role of γ-chain cytokines in driving DN T cell proliferation in steady state and AKI, focusing on IL-2, IL-7, and IL-15.

METHODS

We assessed DN T cell proliferation in vitro in response to IL-2, IL-7, and IL-15, with co-culture experiments using renal tubular epithelial cells (RTECs) and IL-7 blockade. In vivo, wild-type and IL-7r knockout mice were studied to evaluate the impact of IL-7 on DN T cell expansion and kidney function during AKI. Human RTECs confirmed the relevance of these findings.

RESULTS

All three cytokines promoted DN T cell proliferation in vitro, with IL-7 inducing the most robust expansion. Co-culture experiments showed RTECs as a key IL-7 source, and blockade reduced DN T cell expansion. In vivo, IL-7 complexes administered to wild-type mice increased DN T cells and selectively expanded the PD-1+ subset. IL-7 deficiency worsened renal outcomes during AKI. Human RTECs activated peripheral human DN T cells in vitro.

CONCLUSION

These results establish IL-7 as pivotal for DN T cell expansion and highlight RTECs' role in DN T cell homeostasis.

Integrating single-cell RNA and T cell/B cell receptor sequencing with mass cytometry reveals dynamic trajectories of human peripheral immune cells from birth to old age

A comprehensive understanding of the evolution of the immune landscape in humans across the entire lifespan at single-cell transcriptional and protein levels, during development, maturation and senescence is currently lacking. We recruited a total of 220 healthy volunteers from the Shanghai Pudong Cohort (NCT05206643), spanning 13 age groups from 0 to over 90 years, and profiled their peripheral immune cells through single-cell RNA-sequencing coupled with single T cell and B cell receptor sequencing, high-throughput mass cytometry, bulk RNA-sequencing and flow cytometry validation experiments. We revealed that T cells were the most strongly affected by age and experienced the most intensive rewiring in cell-cell interactions during specific age. Different T cell subsets displayed different aging patterns in both transcriptomes and immune repertoires; examples included GNLY+CD8+ effector memory T cells, which exhibited the highest clonal expansion among all T cell subsets and displayed distinct functional signatures in children and the elderly; and CD8+ MAIT cells, which reached their peaks of relative abundance, clonal diversity and antibacterial capability in adolescents and then gradually tapered off. Interestingly, we identified and experimentally verified a previously unrecognized 'cytotoxic' B cell subset that was enriched in children. Finally, an immune age prediction model was developed based on lifecycle-wide single-cell data that can evaluate the immune status of healthy individuals and identify those with disturbed immune functions. Our work provides both valuable insights and resources for further understanding the aging of the immune system across the whole human lifespan.

scRNA + BCR-seq identifies proportions and characteristics of dual BCR B cells in the peritoneal cavity of mice and peripheral blood of healthy human donors across different ages

The increased incidence of inflammatory diseases, infectious diseases, autoimmune disorders, and tumors in elderly individuals is closely associated with several well-established features of immunosenescence, including reduced B cell genesis and dampened immune responses. Recent studies have highlighted the critical role of dual receptor lymphocytes in tumors and autoimmune diseases. This study utilized shared data generated through scRNA-seq + scBCR-seq technology to investigate the presence of dual receptor-expressing B cells in the peritoneum of mouse and peripheral blood of healthy volunteers, and whether there are age-related differences in dual receptor B cell populations. In the peritoneum of mice, a high proportion of B cells expressing dual receptors, predominantly dual κ chains, was observed. Notably, there was an increase in dual BCR B cells in elderly mice. Subsequent analysis revealed that the elevated dual BCR B cells in elderly mice primarily originated from B1 cells.Consistent with the results we observed in healthy volunteers of different ages. Furthermore, these cells exhibited differential expressed genes compared to single BCR B cells, including Vim, Ucp2, and Zcwpw1.These findings support a hypothesis that age-related immune changes encompass not only alterations in B cell numbers but also qualitative changes in BCR diversity. Further exploration of the elevated dual BCR B cells in the elderly population can elucidate their function and their association with immune tolerance, revealing their potential role in maintaining immune surveillance and responding to age-related immune challenges.

Antibody sequence determinants of viral antigen specificity

Throughout life, humans experience repeated exposure to viral antigens through infection and vaccination, resulting in the generation of diverse, antigen-specific antibody repertoires. A paramount feature of antibodies that enables their critical contributions in counteracting recurrent and novel pathogens, and consequently fostering their utility as valuable targets for therapeutic and vaccine development, is the exquisite specificity displayed against their target antigens. Yet, there is still limited understanding of the determinants of antibody-antigen specificity, particularly as a function of antibody sequence. In recent years, experimental characterization of antibody repertoires has led to novel insights into fundamental properties of antibody sequences but has been largely decoupled from at-scale antigen specificity analysis. Here, using the LIBRA-seq technology, we generated a large data set mapping antibody sequence to antigen specificity for thousands of B cells, by screening the repertoires of a set of healthy individuals against 20 viral antigens representing diverse pathogens of biomedical significance. Analysis uncovered virus-specific patterns in variable gene usage, gene pairing, somatic hypermutation, as well as the presence of convergent antiviral signatures across multiple individuals, including the presence of public antibody clonotypes. Notably, our results showed that, for B-cell receptors originating from different individuals but leveraging an identical combination of heavy and light chain variable genes, there is a specific CDRH3 identity threshold above which B cells appear to exclusively share the same antigen specificity. This finding provides a quantifiable measure of the relationship between antibody sequence and antigen specificity and further defines experimentally grounded criteria for defining public antibody clonality.IMPORTANCEThe B-cell compartment of the humoral immune system plays a critical role in the generation of antibodies upon new and repeated pathogen exposure. This study provides an unprecedented level of detail on the molecular characteristics of antibody repertoires that are specific to each of the different target pathogens studied here and provides empirical evidence in support of a 70% CDRH3 amino acid identity threshold in pairs of B cells encoded by identical IGHV:IGL(K)V genes, as a means of defining public clonality and therefore predicting B-cell antigen specificity in different individuals. This is of exceptional importance when leveraging public clonality as a method to annotate B-cell receptor data otherwise lacking antigen specificity information. Understanding the fundamental rules of antibody-antigen interactions can lead to transformative new approaches for the development of antibody therapeutics and vaccines against current and emerging viruses.

scRNA+TCR-seq reveals the pivotal role of dual receptor T lymphocytes in the pathogenesis of Kawasaki disease and during IVIG treatment

Introduction

Kawasaki disease (KD), a common cause of acquired heart disease in children in developed countries, is primarily treated with intravenous immunoglobulin (IVIG), but some children demonstrate IVIG resistance with increased coronary artery injury risk. T cells have been demonstrated to be involved in the pathogenesis of KD and its treatment with IVIG. However, the role and mechanism of dual TCR T lymphocytes in the occurrence of KD and IVIG therapy remain unclear.

Methods

This study, based on scRNA-seq combined with TCR-seq technology, clustered the peripheral blood mononuclear cells of 3 healthy controls and 6 KD patients before and after IVIG treatment. Comparative analysis was conducted to investigate the differences in the proportion of single/dual receptor T cells, the characteristics of CDR3 repertoires, cell types, and the expression of transcription factors among the three groups. The study aimed to explore the correlation between dual TCR T cells and KD as well as IVIG treatment.

Results

In our experimental results, we observed the presence of dual TCR T cells in all three groups. However, compared to the healthy control group and the IVIG-treated group, the KD patients before IVIG treatment exhibited a lower proportion of dual TCR T cells, with variability between samples, ranging from 4% to 15%. Notably, after IVIG treatment, the proportion of dual TCR T cells significantly increased, stabilizing above 12%, and these T cells also exhibited clonal expansion and a preference for V gene usage. In addition we found differences in dual TCR T cell subsets among the three groups, for example, IVIG treatment increases the proportion of dual TCR Treg cells, but it still remains below that of healthy control groups, significantly higher proportions of both dual TCR CD8 central and effector memory T cells in IVIG-treated KD patients, and differences in the expression of transcription factors between single and dual TCR T cells. These results suggest dual TCR T cells correlate with KD and IVIG treatment.

Conclusion

Dual TCR T lymphocytes, especially dual TCR CD8 T cells and Treg cells, play crucial roles in the pathogenesis of KD and during IVIG treatment, providing strong support for further elucidating KD pathogenesis and optimizing treatment strategies.

Invariant VD and DJ Motifs Define a Novel Class of Human Antibodies and TCRs Prototyped by antigen receptors of Dual-Expresser Lymphocytes

INTRODUCTION

Dual-expressing lymphocytes (DEs) are unique immune cells that express both B cell receptors (BCRs, surface antibody) and T cell receptors (TCRs). In type 1 diabetes, DE antibodies are predominated by one antibody (x-mAb), an IgM monoclonal antibody with a germline-encoded CDR3 that recognizes self-reactive TCRs. We explored if x-mAb and its interacting TCRs have distinct structural features.

METHODS

Using bioinformatics, we compared x-mAb and its most common interacting TCRαβ to billions of antigen receptor sequences to determine if they were unique or randomly generated.

RESULTS

X-mAb represents a unique class of human antibodies with a conserved CDR3 sequence (CARx1-4DTAMVYYFYDW), consisting of a fixed DJH motif (DTAMVYYFDYW) paired with various VH genes. A public TCRβ clonotype (CASSPGTEAFF) associated with x-mAb on DEs features two invariant segments, VβD (CASSPGT) and DJβ (PGTEAFF), key to two large families of public TCRβ clonotypes-CASSPGT-Jβx and CASSPGT-Jβx-formed by recombining the VβD motif with Jβ genes and the DJβ motif with Vβ genes. B cells also use CASSPGT as a VHD motif for public IGH clonotypes (CASSPGT-Jβx).

DISCUSSION

DEs, unlike conventional T and B cells, use invariant motifs to create public antibodies and TCRs, a trait previously seen only in cartilaginous fish.

Inhibition of DEF-p65 Interactions as a Potential Avenue to Suppress Tumor Growth in Pancreatic Cancer

The limited success of current targeted therapies for pancreatic cancer underscores an urgent demand for novel treatment modalities. The challenge in mitigating this malignancy can be attributed to the digestive organ expansion factor (DEF), a pivotal yet underexplored factor in pancreatic tumorigenesis. The study uses a blend of in vitro and in vivo approaches, complemented by the theoretical analyses, to propose DEF as a promising anti-tumor target. Analysis of clinical samples reveals that high expression of DEF is correlated with diminished survival in pancreatic cancer patients. Crucially, the depletion of DEF significantly impedes tumor growth. The study further discovers that DEF binds to p65, shielding it from degradation mediated by the ubiquitin-proteasome pathway in cancer cells. Based on these findings and computational approaches, the study formulates a DEF-mimicking peptide, peptide-031, designed to disrupt the DEF-p65 interaction. The effectiveness of peptide-031 in inhibiting tumor proliferation has been demonstrated both in vitro and in vivo. This study unveils the oncogenic role of DEF while highlighting its prognostic value and therapeutic potential in pancreatic cancer. In addition, peptide-031 is a promising therapeutic agent with potent anti-tumor effects.

Autoimmune CD8+ T cells in type 1 diabetes: from single-cell RNA sequencing to T-cell receptor redirection

Type 1 diabetes (T1D) is an organ-specific autoimmune disease caused by pancreatic β cell destruction and mediated primarily by autoreactive CD8+ T cells. It has been shown that only a small number of stem cell-like β cell-specific CD8+ T cells are needed to convert normal mice into T1D mice; thus, it is likely that T1D can be cured or significantly improved by modulating or altering self-reactive CD8+ T cells. However, stem cell-type, effector and exhausted CD8+ T cells play intricate and important roles in T1D. The highly diverse T-cell receptors (TCRs) also make precise and stable targeted therapy more difficult. Therefore, this review will investigate the mechanisms of autoimmune CD8+ T cells and TCRs in T1D, as well as the related single-cell RNA sequencing (ScRNA-Seq), CRISPR/Cas9, chimeric antigen receptor T-cell (CAR-T) and T-cell receptor-gene engineered T cells (TCR-T), for a detailed and clear overview. This review highlights that targeting CD8+ T cells and their TCRs may be a potential strategy for predicting or treating T1D.

In silico design of high-affinity antigenic peptides for HLA-B44

Cancer cell-killing by CD8+ T cells demands effective tumor antigen presentation by human leukocyte antigen class I (HLA-I) molecules. Screening and designing highly immunogenic neoantigens require quantitative computations to reliably predict HLA-peptide binding affinities. Here, with all-atom molecular dynamics (MD) simulations and free energy perturbation (FEP) methods, we design a collection of antigenic peptide candidates through in silico mutagenesis studies on immunogenic neoantigens, yielding enhanced binding affinities to HLA-B*44:02. In-depth structural dissection shows that introducing positively charged residues such as arginine to position 6 or lysine to position 7 of the candidates triggers conformational shifts in both peptides and the antigen-binding groove of the HLA, following the "induced-fit" mechanism. Enhancement in binding affinities compared to the wild-type was found in three out of five mutated candidates. The HLA pocket, capable of accommodating positively charged residues in positions from 5 to 7, is designated as the "dynamic pocket". Taken together, we showcase an effective structure-based binding affinity optimization framework for antigenic peptides of HLA-B*44:02 and underscore the importance of dynamic nature of the antigen-binding groove in concert with the anchoring motifs. This work provides structural insights for rational design of favorable HLA-peptide bindings and future developments in neoantigen-based therapeutics.

Chimeric Antigen Receptor (CAR)-Based Cell Therapy for Type 1 Diabetes Mellitus (T1DM); Current Progress and Future Approaches

Type 1 diabetes mellitus (T1DM) is an autoimmune disease that destroys insulin-producing pancreatic β-cells. Insulin replacement therapy is currently the mainstay of treatment for T1DM; however, treatment with insulin does not ameliorate disease progression, as dysregulated immune response and inflammation continue to cause further pancreatic β-cell degradation. Therefore, shifting therapeutic strategies toward immunomodulating approaches could be effective to prevent and reverse disease progression. Different immune-modulatory therapies could be used, e.g., monoclonal-based immunotherapy, mesenchymal stem cell, and immune cell therapy. Since immune-modulatory approaches could have a systemic effect on the immune system and cause toxicity, more specific treatment options should target the immune response against pancreatic β-cells. In this regard, chimeric antigen receptor (CAR)-based immunotherapy could be a promising candidate for modulation of dysregulated immune function in T1DM. CAR-based therapy has previously been approved for a number of hematologic malignancies. Nevertheless, there is renewed interest in CAR T cells' " off-the-shelf " treatment for T1DM. Several pre-clinical studies demonstrated that redirecting antigen-specific CAR T cells, especially regulatory CAR T cells (CAR Tregs), toward the pancreatic β-cells, could prevent diabetes onset and progression in diabetic mice models. Here, we aim to review the current progress of CAR-based immune-cell therapy for T1DM and the corresponding challenges, with a special focus on designing CAR-based immunomodulatory strategies to improve its efficacy in the treatment of T1DM.

Single-cell RNA sequencing unveils the hidden powers of zebrafish kidney for generating both hematopoiesis and adaptive antiviral immunity

The vertebrate kidneys play two evolutionary conserved roles in waste excretion and osmoregulation. Besides, the kidney of fish is considered as a functional ortholog of mammalian bone marrow that serves as a hematopoietic hub for generating blood cell lineages and immunological responses. However, knowledge about the properties of kidney hematopoietic cells, and the functionality of the kidney in fish immune systems remains to be elucidated. To this end, our present study generated a comprehensive atlas with 59 hematopoietic stem/progenitor cell (HSPC) and immune-cells types from zebrafish kidneys via single-cell transcriptome profiling analysis. These populations included almost all known cells associated with innate and adaptive immunity, and displayed differential responses to viral infection, indicating their diverse functional roles in antiviral immunity. Remarkably, HSPCs were found to have extensive reactivities to viral infection, and the trained immunity can be effectively induced in certain HSPCs. In addition, the antigen-stimulated adaptive immunity can be fully generated in the kidney, suggesting the kidney acts as a secondary lymphoid organ. These results indicated that the fish kidney is a dual-functional entity with functionalities of both primary and secondary lymphoid organs. Our findings illustrated the unique features of fish immune systems, and highlighted the multifaced biology of kidneys in ancient vertebrates.

Shift in the B cell subsets between children with type 1 diabetes and/or celiac disease

Our purpose was to characterize the pattern of B cell subsets in children with a combined diagnosis of type 1 diabetes (T1D) and celiac disease (C) since children with single or double diagnosis of these autoimmune diseases may differ in peripheral B cell subset phenotype patterns. B cells were analyzed with flow cytometry for the expression of differentiation/maturation markers to identify transitional, naive, and memory B cells. Transitional (CD24hiCD38hiCD19+) and memory Bregs (mBregs; CD24hiCD27+CD19+, CD1d+CD27+CD19+, and CD5+CD1d+CD19+) were classified as B cells with regulatory capacity. Children with a combined diagnosis of T1D and C showed a pattern of diminished peripheral B cell subsets. The B cells compartment in children with combined diagnosis had higher percentages of memory B subsets and Bregs, including activated subsets, compared to children with either T1D or C. Children with combined diagnosis had a lower percentage of naive B cells (CD27-CD19+; IgD+CD19+) and an increased percentage of memory B cells (CD27+CD19+; IgD-CD19+). A similar alteration was seen among the CD39+ expressing naive and memory B cells. Memory Bregs (CD1d+CD27+CD19+) were more frequent, contrary to the lower percentage of CD5+ transitional Bregs in children with a combined diagnosis. In children with either T1D or C, the peripheral B cell compartment was dominated by naive cells. Differences in the pattern of heterogeneous peripheral B cell repertoire subsets reflect a shifting in the B cell compartment between children with T1D and/or C. This is an immunological challenge of impact on the pathophysiology of these autoimmune diseases.

The Deinococcus protease PprI senses DNA damage by directly interacting with single-stranded DNA

Bacteria have evolved various response systems to adapt to environmental stress. A protease-based derepression mechanism in response to DNA damage was characterized in Deinococcus, which is controlled by the specific cleavage of repressor DdrO by metallopeptidase PprI (also called IrrE). Despite the efforts to document the biochemical, physiological, and downstream regulation of PprI-DdrO, the upstream regulatory signal activating this system remains unclear. Here, we show that single-stranded DNA physically interacts with PprI protease, which enhances the PprI-DdrO interactions as well as the DdrO cleavage in a length-dependent manner both in vivo and in vitro. Structures of PprI, in its apo and complexed forms with single-stranded DNA, reveal two DNA-binding interfaces shaping the cleavage site. Moreover, we show that the dynamic monomer-dimer equilibrium of PprI is also important for its cleavage activity. Our data provide evidence that single-stranded DNA could serve as the signal for DNA damage sensing in the metalloprotease/repressor system in bacteria. These results also shed light on the survival and acquired drug resistance of certain bacteria under antimicrobial stress through a SOS-independent pathway.

Distinct immune escape and microenvironment between RG-like and pri-OPC-like glioma revealed by single-cell RNA-seq analysis

The association of neurogenesis and gliogenesis with glioma remains unclear. By conducting single-cell RNA-seq analyses on 26 gliomas, we reported their classification into primitive oligodendrocyte precursor cell (pri-OPC)-like and radial glia (RG)-like tumors and validated it in a public cohort and TCGA glioma. The RG-like tumors exhibited wild-type isocitrate dehydrogenase and tended to carry EGFR mutations, and the pri-OPC-like ones were prone to carrying TP53 mutations. Tumor subclones only in pri-OPC-like tumors showed substantially down-regulated MHC-I genes, suggesting their distinct immune evasion programs. Furthermore, the two subgroups appeared to extensively modulate glioma-infiltrating lymphocytes in distinct manners. Some specific genes not expressed in normal immune cells were found in glioma-infiltrating lymphocytes. For example, glial/glioma stem cell markers OLIG1/PTPRZ1 and B cell-specific receptors IGLC2/IGKC were expressed in pri-OPC-like and RG-like glioma-infiltrating lymphocytes, respectively. Their expression was positively correlated with those of immune checkpoint genes (e.g., LGALS33) and poor survivals as validated by the increased expression of LGALS3 upon IGKC overexpression in Jurkat cells. This finding indicated a potential inhibitory role in tumor-infiltrating lymphocytes and could provide a new way of cancer immune evasion.

B-Cell Receptor Repertoire: Recent Advances in Autoimmune Diseases

In the field of contemporary medicine, autoimmune diseases (AIDs) are a prevalent and debilitating group of illnesses. However, they present extensive and profound challenges in terms of etiology, pathogenesis, and treatment. A major reason for this is the elusive pathophysiological mechanisms driving disease onset. Increasing evidence suggests the indispensable role of B cells in the pathogenesis of autoimmune diseases. Interestingly, B-cell receptor (BCR) repertoires in autoimmune diseases display a distinct skewing that can provide insights into disease pathogenesis. Over the past few years, advances in high-throughput sequencing have provided powerful tools for analyzing B-cell repertoire to understand the mechanisms during the period of B-cell immune response. In this paper, we have provided an overview of the mechanisms and analytical methods for generating BCR repertoire diversity and summarize the latest research progress on BCR repertoire in autoimmune diseases, including systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), primary Sjögren's syndrome (pSS), multiple sclerosis (MS), and type 1 diabetes (T1D). Overall, B-cell repertoire analysis is a potent tool to understand the involvement of B cells in autoimmune diseases, facilitating the creation of innovative therapeutic strategies targeting specific B-cell clones or subsets.

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.

A biphenotypic lymphocyte subset displays both T- and B-cell functionalities

T cell/B cell mixed phenotypic lymphocytes have been observed in different disease contexts, yet their presence and function in physiological conditions remain elusive. Here, we provide evidence for the existence of a lymphocyte subset endogenously expressing both T- and B-cell lineage markers in mice. The majority of these T/B phenotypic lymphocytes (CD3+CD19+) show an origin of pro/pre B cells and distribute widely in mouse bone marrow, lymph nodes, spleen, and peripheral blood. Functional assays show that these biphenotypic lymphocytes can be activated through stimulating TCR or BCR signaling pathways. Moreover, we show that these cells actively participate both the humoral and cellular immune responses elicited by vaccination. Compared to conventional T cells, these biphenotypic lymphocytes can secrete a higher level of IL-2 but a lower level of TNF-α upon antigen specific stimulation. An equivalent lymphocyte subset is found in freshly isolated human PBMCs and exhibits similar functionality, albeit at a lower frequency than in mice.

Reproducibility in cytometry: Signals analysis and its connection to uncertainty quantification

Signals analysis for cytometry remains a challenging task that has a significant impact on uncertainty. Conventional cytometers assume that individual measurements are well characterized by simple properties such as the signal area, width, and height. However, these approaches have difficulty distinguishing inherent biological variability from instrument artifacts and operating conditions. As a result, it is challenging to quantify uncertainty in the properties of individual cells and perform tasks such as doublet deconvolution. We address these problems via signals analysis techniques that use scale transformations to: (I) separate variation in biomarker expression from effects due to flow conditions and particle size; (II) quantify reproducibility associated with a given laser interrogation region; (III) estimate uncertainty in measurement values on a per-event basis; and (IV) extract the singlets that make up a multiplet. The key idea behind this approach is to model how variable operating conditions deform the signal shape and then use constrained optimization to "undo" these deformations for measured signals; residuals to this process characterize reproducibility. Using a recently developed microfluidic cytometer, we demonstrate that these techniques can account for instrument and measurand induced variability with a residual uncertainty of less than 2.5% in the signal shape and less than 1% in integrated area.

Tracking DNA-based antigen-specific T cell receptors during progression to type 1 diabetes

T cells targeting self-proteins are important mediators in autoimmune diseases. T cells express unique cell-surface receptors (TCRs) that recognize peptides presented by major histocompatibility molecules. TCRs have been identified from blood and pancreatic islets of individuals with type 1 diabetes (T1D). Here, we tracked ~1700 known antigen-specific TCR sequences, islet antigen or viral reactive, in bulk TCRβ sequencing from longitudinal blood DNA samples in at-risk cases who progressed to T1D, age/sex/human leukocyte antigen-matched controls, and a new-onset T1D cohort. Shared and frequent antigen-specific TCRβ sequences were identified in all three cohorts, and viral sequences were present across all ages. Islet sequences had different patterns of accumulation based upon antigen specificity in the at-risk cases. Furthermore, 73 islet-antigen TCRβ sequences were present in higher frequencies and numbers in T1D samples relative to controls. The total number of these disease-associated TCRβ sequences inversely correlated with age at clinical diagnosis, indicating the potential to use disease-relevant TCR sequences as biomarkers in autoimmune disorders.

Translational immunology: Applying fundamental discoveries to human health and autoimmune diseases

Studying the human immune system is challenging. These challenges stem from the complexity of the immune system itself, the heterogeneity of the immune system between individuals, and the many factors that lead to this heterogeneity including the influence of genetics, environment, and immune experience. Studies of the human immune system in the context of disease are increased in complexity as multiple combinations and variations in immune pathways can lead to a single disease. Thus, although individuals with a disease may share clinical features, the underlying disease mechanisms and resulting pathophysiology can be diverse among individuals with the same disease diagnosis. This has consequences for the treatment of diseases, as no single therapy will work for everyone, therapeutic efficacy varies among patients, and targeting a single immune pathway is rarely 100% effective. This review discusses how to address these challenges by identifying and managing the sources of variation, improving access to high-quality, well-curated biological samples by building cohorts, applying new technologies such as single-cell omics and imaging technologies to interrogate samples, and bringing to bear computational expertise in conjunction with immunologists and clinicians to interpret those results. The review has a focus on autoimmune diseases, including rheumatoid arthritis, MS, systemic lupus erythematosus, and type 1 diabetes, but its recommendations are also applicable to studies of other immune-mediated diseases.

Renal tubular epithelial cells are constitutive non-cognate stimulators of resident T cells

Understanding the roles of different cell types in regulating T cell homeostasis in various tissues is critical for understanding adaptive immunity. Here, we show that RTECs (renal tubular epithelial cells) are intrinsically programmed to polyclonally stimulate proliferation of kidney αβ T cells by a cell-cell contact mechanism that is major histocompatibility complex (MHC) independent and regulated by CD155, αVβ3-integrin, and vitronectin. Peripheral CD4 and CD8 are resistant to RTEC-mediated stimulation, while the minor subset of double-negative (DN) T cells are responsive. This functional property of RTEC is discovered by using a coculture system that recapitulates spontaneous in vivo polyclonal proliferation of kidney T cells, which are mainly comprised of central memory T (TCM) and effector memory T (TEM) cells. This robust cell-intrinsic stimulatory role of RTECs could be underlying the steady-state spontaneous proliferation of kidney T cells. The results have conceptual implications for understanding roles of different cell types in regulating systemic and organ-specific T cell homeostasis.

scRNA-seq revealed the special TCR β & α V(D)J allelic inclusion rearrangement and the high proportion dual (or more) TCR-expressing cells

Allelic exclusion, one lymphocyte expresses one antigen receptor, is a fundamental mechanism of immunological self-tolerance and highly specific immune responses to pathogens. However, the phenomenon of V(D)J allelic inclusion (incomplete allelic exclusion or allelic escape) rearrangement and dual TCR T cells have been discovered by multiple laboratories. Despite continuous new discoveries, the proportion and underlying mechanism of dual TCR has been puzzling immunologists. In this study, we observed the presence of single T cells expressing multiple TCR chains in all samples, with the proportion of 15%, 10%, and 20% in the human thymus, human peripheral blood, and mouse lymphoid organs, respectively. The proportion of T cells possessing multiple T-cell receptors (TCR) varied significantly in different physiological states and developmental stages. By analyzing RSS category, RSS direction, and V(D)J gene position at TR locus of T cells which contain multiple TCR chains, we creatively found that one of TCR β (or TCR α) should originate from the transcription of V(D)J combination in T-cell receptor excision circle (TREC) formed after the twice successful rearrangement in the same chromosome. Moreover, human V30 (or mouse V31) gene may participate in reverse recombination and transcription to prevent allelic exclusion. In general, high proportion of T cells with multiple TCR at the transcriptome level was first made public, and we proposed a novel mechanism of secondary (or more) TCR rearrangement on a single chromosome. Our findings also indicated that the single-cell sequencing data should be classified according to the single, multiple, and abnormal TCR when analyzing the T-cell repertoire.

HLA-II immunopeptidome profiling and deep learning reveal features of antigenicity to inform antigen discovery

CD4+ T cell responses are exquisitely antigen specific and directed toward peptide epitopes displayed by human leukocyte antigen class II (HLA-II) on antigen-presenting cells. Underrepresentation of diverse alleles in ligand databases and an incomplete understanding of factors affecting antigen presentation in vivo have limited progress in defining principles of peptide immunogenicity. Here, we employed monoallelic immunopeptidomics to identify 358,024 HLA-II binders, with a particular focus on HLA-DQ and HLA-DP. We uncovered peptide-binding patterns across a spectrum of binding affinities and enrichment of structural antigen features. These aspects underpinned the development of context-aware predictor of T cell antigens (CAPTAn), a deep learning model that predicts peptide antigens based on their affinity to HLA-II and full sequence of their source proteins. CAPTAn was instrumental in discovering prevalent T cell epitopes from bacteria in the human microbiome and a pan-variant epitope from SARS-CoV-2. Together CAPTAn and associated datasets present a resource for antigen discovery and the unraveling genetic associations of HLA alleles with immunopathologies.

Differences in the Tumor Microenvironment of EBV-Associated Gastric Cancers Revealed Using Single-Cell Transcriptome Analysis

Epstein-Barr virus (EBV) is a gamma-herpesvirus associated with nearly 10% of gastric cancers (GCs). These EBV-associated GCs (EBVaGCs) are molecularly, histopathologically, and clinically distinct from EBV-negative GCs (EBVnGCs). While viral genes in EBVaGCs contribute to the carcinogenesis process, viral proteins also represent foreign antigens that could trigger enhanced immune responses compared to EBVnGCs. Despite prior investigations of the EBVaGC tumor microenvironment (TME), the cellular composition has not been thoroughly explored. In this study, cellular subpopulations overrepresented in EBVaGCs were identified and molecularly characterized. Genes consistently expressed across both bulk tumor and single-cell RNA sequencing data were highlighted, with the expression across the identified cellular subpopulations analyzed. As expected, based on existing histopathological analysis, EBVaGC is characterized by abundant lymphocytic infiltration of the stroma. Our molecular analysis identified three unique immune cell subpopulations in EBVaGC: T and B cells expressing high levels of proliferation markers and B cells expressing T cell features. The proliferating T cell cluster also expressed markers of follicular T helper cells. Overall, EBVaGC also exhibited unique features indicative of a higher inflammatory response. These substantial differences within the TME suggest that further detailed exploration of the cellular composition of EBVaGCs is needed, which may identify cellular subpopulations and phenotypes associated with patient outcomes.

Novel Inhibitory Role of Fenofibric Acid by Targeting Cryptic Site on the RBD of SARS-CoV-2

The emergence of the recent pandemic causing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has created an alarming situation worldwide. It also prompted extensive research on drug repurposing to find a potential treatment for SARS-CoV-2 infection. An active metabolite of the hyperlipidemic drug fenofibrate (also called fenofibric acid or FA) was found to destabilize the receptor-binding domain (RBD) of the viral spike protein and therefore inhibit its binding to human angiotensin-converting enzyme 2 (hACE2) receptor. Despite being considered as a potential drug candidate for SARS-CoV-2, FA's inhibitory mechanism remains to be elucidated. We used molecular dynamics (MD) simulations to investigate the binding of FA to the RBD of the SARS-CoV-2 spike protein and revealed a potential cryptic FA binding site. Free energy calculations were performed for different FA-bound RBD complexes. The results suggest that the interaction of FA with the cryptic binding site of RBD alters the conformation of the binding loop of RBD and effectively reduces its binding affinity towards ACE2. Our study provides new insights for the design of SARS-CoV-2 inhibitors targeting cryptic sites on the RBD of SARS-CoV-2.

Single-cell profiles reveal tumor cell heterogeneity and immunosuppressive microenvironment in Waldenström macroglobulinemia

BACKGROUND

Waldenström macroglobulinemia (WM) is a rare and incurable indolent B-cell malignancy. The molecular pathogenesis and the role of immunosuppressive microenvironment in WM development are still incompletely understood.

METHODS

The multicellular ecosystem in bone marrow (BM) of WM were delineated by single-cell RNA-sequencing (scRNA-seq) and investigated the underlying molecular characteristics.

RESULTS

Our data uncovered the heterogeneity of malignant cells in WM, and investigated the kinetic co-evolution of WM and immune cells, which played pivotal roles in disease development and progression. Two novel subpopulations of malignant cells, CD19+CD3+ and CD138+CD3+, co-expressing T-cell marker genes were identified at single-cell resolution. Pseudotime-ordered analysis elucidated that CD19+CD3+ malignant cells presented at an early stage of WM-B cell differentiation. Colony formation assay further identified that CD19+CD3+ malignant cells acted as potential WM precursors. Based on the findings of T cell marker aberrant expressed on WM tumor cells, we speculate the long-time activation of tumor antigen-induced immunosuppressive microenvironment that is involved in the pathogenesis of WM. Therefore, our study further investigated the possible molecular mechanism of immune cell dysfunction. A precursor exhausted CD8-T cells and functional deletion of NK cells were identified in WM, and CD47 would be a potential therapeutic target to reverse the dysfunction of immune cells.

CONCLUSIONS

Our study facilitates further understanding of the biological heterogeneity of tumor cells and immunosuppressive microenvironment in WM. These data may have implications for the development of novel immunotherapies, such as targeting pre-exhausted CD8-T cells in WM.

Peripheral immune landscape and natural killer-like B cells in human Vogt-Koyanagi-Harada disease

Vogt-Koyanagi-Harada (VKH) disease is a systemic autoimmune disorder threatening the eyesight. The pathogenic mechanisms and biomarkers reflecting disease severity and predicting treatment response require further exploration. Here, we performed a single-cell analysis of peripheral blood mononuclear cells (PBMC) obtained from eight patients with VKH disease and eight healthy controls to comprehensively delineate the changes in VKH disease. We showed a mixture of inflammation, effector, and exhausted states for PBMCs in VKH disease. Notably, our study implicated a newly identified B cell subset, natural killer-like B cells (K-BC) characterized by expressing CD19 and CD56, was correlated with VKH disease. K-BCs expanded in VKH disease, fell back after effective treatment, and promoted the differentiation of pathogenic T cells. Overall, we mapped the peripheral immune cell atlas in VKH disease and indicated the pathogenic role and potential value in predicting treatment response of K-BCs.

Functional antibodies exhibit light chain coherence

The vertebrate adaptive immune system modifies the genome of individual B cells to encode antibodies that bind particular antigens1. In most mammals, antibodies are composed of heavy and light chains that are generated sequentially by recombination of V, D (for heavy chains), J and C gene segments. Each chain contains three complementarity-determining regions (CDR1-CDR3), which contribute to antigen specificity. Certain heavy and light chains are preferred for particular antigens2-22. Here we consider pairs of B cells that share the same heavy chain V gene and CDRH3 amino acid sequence and were isolated from different donors, also known as public clonotypes23,24. We show that for naive antibodies (those not yet adapted to antigens), the probability that they use the same light chain V gene is around 10%, whereas for memory (functional) antibodies, it is around 80%, even if only one cell per clonotype is used. This property of functional antibodies is a phenomenon that we call light chain coherence. We also observe this phenomenon when similar heavy chains recur within a donor. Thus, although naive antibodies seem to recur by chance, the recurrence of functional antibodies reveals surprising constraint and determinism in the processes of V(D)J recombination and immune selection. For most functional antibodies, the heavy chain determines the light chain.

Binding Affinity Calculations of Gluten Peptides to HLA Risk Modifiers: DQ2.5 versus DQ7.5

Free energy perturbation (FEP) calculations can predict relative binding affinities of an antigen and its point mutants to the same human leukocyte antigen (HLA) with high accuracy (e.g., within 1.0 kcal/mol to experiment); however, a more challenging task is to compare binding affinities of wholly different antigens binding to completely different HLAs using FEP. Researchers have used a variety of different FEP schemes to compute and compare absolute binding affinities, with varied success. Here, we propose and assess a unifying scheme to compute the relative binding affinities of different antigens binding to completely different HLAs using absolute binding affinity FEP calculations. We apply our affinity calculation technique to HLA-antigen-T-cell receptor (TCR) systems relevant to celiac disease (CeD) by investigating binding affinity differences between HLA-DQ2.5 (enhanced CeD risk) and HLA-DQ7.5 (CeD protective) in the binary (HLA-gliadin) and ternary (HLA-gliadin-TCR) binding complexes for three gliadin derived epitopes: glia-α1, glia-α2, and glia-ω1. Based on FEP calculations with our carefully designed thermodynamic cycles, we demonstrate that HLA-DQ2.5 has higher binding affinity than HLA-DQ7.5 for gliadin and enhanced binding affinity with a common TCR, agreeing with known results that the HLA-DQ2.5 serotype exhibits increased risk for CeD. Our findings reveal that our proposed absolute binding affinity FEP method is appropriate for predicting HLA binding for disparate antigens with different genotypes. We also discuss atomic-level details of HLA genotypes interacting with gluten peptides and TCRs in regard to the pathogenesis of CeD.

Distinguishing cell-cell complexes from dual lineage cells using single-cell transcriptomics is not trivial

In their recent correspondence, Jie et al. strongly defend that the DE cell population they discovered are always dual lineage co-expressing cells and not complexes of B cells and T cells, which we have previously described as frequently present in single-cell RNA sequencing data. Here, we respond to the specific arguments made in their correspondence. Specifically, we demonstrate that the presence of a gene signature in a given cell population is not enough to ascertain that it does not contain cell-cell complexes, or that it represents a biologically distinct cell type. We also show that the gene signature of DE cells contains several genes from the myeloid lineage, suggesting either that their DE cells are a triple-lineage co-expressing cell, or a three-component cell aggregate. Finally, we identify multiple transcriptomic features of DE cells that correspond to B cell-T cell complexes, namely the presence of lower average expression of B- and T-cell specific genes, and a higher number of detected genes per cell. Taken together, our results demonstrate that solely based on their scRNAseq profile, it is not possible to ascertain that DE cells are dual expressing cells and not cell-cell complexes.

Expression of unique gene signature distinguishes TCRαβ+ /BCR+ dual expressers from CD3+ CD14+ doublets

Increasing evidence shows pathophysiological significance of rare immune cells, necessitating the need for reliable and proper methods for their detection and analysis. We have recently identified a new lymphocyte that coexpresses lineage markers of T- and B-cells including T cell receptor and B cell receptor (called dual expressers, DEs). Because of the peculiar phenotype of DEs, we used multiple techniques to authenticate their identity (fluorescence-activated cell sorting [FACS], scRNAseq, EBV cell lines, and imaging flow cytometry). In an recent article published in this journal, Burel and colleagues successfully detected DEs using FACS and imaging microscopy. Yet they claimed, based on the profile of what they called naturally occurring CD3+ CD14+ T cell/monocyte complexes that the scRNAseq signature of DEs resembles that of cell-cell complexes. Serious flaws in their analysis, however, invalidate their conclusions. Unlike the CD3+ CD14+ complexes, DEs have a distinct identity due to expression of a unique set of signature genes. Without a clear explanation, Burel and colleagues excluded these genes from their analysis, thereby effectively stripped DEs from their identity. Inclusion of these genes as described in this communication restores the identity of DEs. Moreover, contrary to the claim of Burel and colleagues, B- and T-cell specific genes are similarly expressed in DE cells.

Characterization of Pipefish Immune Cell Populations Through Single-Cell Transcriptomics

Teleost adaptive immune systems have evolved with more flexibility than previously assumed. A particularly enigmatic system to address immune system modifications in the evolutionary past is represented by the Syngnathids, the family of pipefishes, seahorses and seadragons. These small fishes with their unique male pregnancy have lost the spleen as an important immune organ as well as a functional major histocompatibility class II (MHC II) pathway. How these evolutionary changes have impacted immune cell population dynamics have up to this point remained unexplored. Here, we present the first immune cell repertoire characterization of a syngnathid fish (Syngnathus typhle) using single-cell transcriptomics. Gene expression profiles of individual cells extracted from blood and head-kidney clustered in twelve putative cell populations with eight belonging to those with immune function. Upregulated cell marker genes identified in humans and teleosts were used to define cell clusters. While the suggested loss of CD4+ T-cells accompanied the loss of the MHC II pathway was supported, the upregulation of specific subtype markers within the T-cell cluster indicates subpopulations of regulatory T-cells (il2rb) and cytotoxic T-cells (gzma). Utilizing single-cell RNA sequencing this report is the first to characterize immune cell populations in syngnathids and provides a valuable foundation for future cellular classification and experimental work within the lineage.

Cytotoxicity-Related Gene Expression and Chromatin Accessibility Define a Subset of CD4+ T Cells That Mark Progression to Type 1 Diabetes

Type 1 diabetes in children is heralded by a preclinical phase defined by circulating autoantibodies to pancreatic islet antigens. How islet autoimmunity is initiated and then progresses to clinical diabetes remains poorly understood. Only one study has reported gene expression in specific immune cells of children at risk associated with progression to islet autoimmunity. We analyzed gene expression with RNA sequencing in CD4+ and CD8+ T cells, natural killer (NK) cells, and B cells, and chromatin accessibility by assay for transposase-accessible chromatin sequencing (ATAC-seq) in CD4+ T cells, in five genetically at risk children with islet autoantibodies who progressed to diabetes over a median of 3 years ("progressors") compared with five children matched for sex, age, and HLA-DR who had not progressed ("nonprogressors"). In progressors, differentially expressed genes (DEGs) were largely confined to CD4+ T cells and enriched for cytotoxicity-related genes/pathways. Several top-ranked DEGs were validated in a semi-independent cohort of 13 progressors and 11 nonprogressors. Flow cytometry confirmed that progression was associated with expansion of CD4+ cells with a cytotoxic phenotype. By ATAC-seq, progression was associated with reconfiguration of regulatory chromatin regions in CD4+ cells, some linked to differentially expressed cytotoxicity-related genes. Our findings suggest that cytotoxic CD4+ T cells play a role in promoting progression to type 1 diabetes.

Single-Cell-Based High-Throughput Ig and TCR Repertoire Sequencing Analysis in Rhesus Macaques

Recent advancements in microfluidics and high-throughput sequencing technologies have enabled recovery of paired H and L chains of Igs and VDJ and VJ chains of TCRs from thousands of single cells simultaneously in humans and mice. Despite rhesus macaques being one of the most well-studied model organisms for the human adaptive immune response, high-throughput single-cell immune repertoire sequencing assays are not yet available due to the complexity of these polyclonal receptors. We used custom primers that capture all known rhesus macaque Ig and TCR isotypes and chains that are fully compatible with a commercial solution for single-cell immune repertoire profiling. Using these rhesus-specific assays, we sequenced Ig and TCR repertoires in >60,000 cells from cryopreserved rhesus PBMCs, splenocytes, and FACS-sorted B and T cells. We were able to recover every Ig isotype and TCR chain, measure clonal expansion in proliferating T cells, and pair Ig and TCR repertoires with gene expression profiles of the same single cells. Our results establish the ability to perform high-throughput immune repertoire analysis in rhesus macaques at the single-cell level.

Detection, sorting, and immortalization of dual expresser lymphocytes from human peripheral blood samples

This protocol describes how to identify Dual Expressers (DEs), a rare type of lymphocyte that co-expresses B-cell receptors and T-cell receptors, by flow cytometry using a cocktail of four antibodies. It also shows the subsequent gating strategy for detecting and sorting DEs and the generation of EBV-immortalized DE lymphoblastoid cell lines and clones for antibody production and cloning antigen receptors. Use of this protocol maximizes detection of DEs and minimizes inclusion of doublets.

For complete details on the use and execution of this protocol, please refer to Ahmed et al. (2019).

•

Protocol to detect DEs and distinguish from T and B cells by flow cytometry

•

Sequential gating strategy to maximize and optimize detection of DEs

•

Strategy for efficient FACS-sorting and checking purity of sorted DEs

•

Protocol for generating EBV-immortalized DE lines and clones

Cellular heterogeneity of circulating CD4+CD8+ double-positive T cells characterized by single-cell RNA sequencing

Circulating CD4+CD8+ double-positive (DP) T cells are associated with a variety of disease states. However, unlike conventional T cells, the composition of this population is poorly understood. Here, we used single-cell RNA sequencing (scRNA-seq) to analyze the composition and characteristics of the DP T cell population circulating in the peripheral blood of cynomolgus monkeys. We found that circulating DP T cells not only contain a large number of naïve cells, but also comprise a heterogeneous population (CD4 CTL-, Eomes+ Tr1-, Th2-, Th17-, Tfh-, Treg-, CD8 CTL-, and innate-like cells) with multiple potential functions. Flow cytometry analysis revealed that a substantial number of the naïve DP T cells expressed CD8αβ, as well as CD8αα, along with high expression of CD31. Moreover, the CD4hiCD8lo and CD4hiCD8hi populations, which express high levels of the CD4 coreceptor, comprised subsets characterized by helper and regulatory functions, some of which also exhibited cytotoxic functions. By contrast, the CD4loCD8hi population with high CD8 coreceptor expression comprised a subset characterized by CD8 CTL- and innate-like properties. Taken together, the data show that scRNA-seq analysis identified a more diverse subset of the circulating DP cells than is currently known, despite this population being very small.

Diabetes type 1: Can it be treated as an autoimmune disorder?

Type 1 Diabetes Mellitus (T1DM) is characterized by progressive autoimmune-mediated destruction of the pancreatic beta-cells leading to insulin deficiency and hyperglycemia. It is associated with significant treatment burden and necessitates life-long insulin therapy. The role of immunotherapy in the prevention and management of T1DM is an evolving area of interest which has the potential to alter the natural history of this disease.In this review, we give insight into recent clinical trials related to the use of immunotherapeutic approaches for T1DM, such as proinflammatory cytokine inhibition, cell-depletion and cell-therapy approaches, autoantigen-specific treatments and stem cell therapies. We highlight the timing of intervention, aspects of therapy including adverse effects and the emergence of a novel lymphocyte crucial in T1DM autoimmunity. We also discuss the role of cardiac autoimmunity and its link to excess CVD risk in T1DM.We conclude that significant advances have been made in development of immunotherapeutic targets and agents for the treatment and prevention of T1DM. These immune-based therapies promise preservation of beta-cells and decreasing insulin dependency. In their current state, immunotherapeutic approaches cannot yet halt the progression from a preclinical state to overt T1DM nor can they replace standard insulin therapy in existing T1DM. It remains to be seen whether immunotherapy will ultimately play a key role in the prevention of progression to overt T1DM and whether it may find a place in our therapeutic armamentarium to improve clinical outcomes and quality of life in established T1DM.

Planar Boronic Graphene and Nitrogenized Graphene Heterostructure for Protein Stretch and Confinement

Single-molecule techniques such as electron tunneling and atomic force microscopy have attracted growing interests in protein sequencing. For these methods, it is critical to refine and stabilize the protein sample to a “suitable mode” before applying a high-fidelity measurement. Here, we show that a planar heterostructure comprising boronic graphene (BC₃) and nitrogenized graphene (C₃N) sandwiched stripe (BC₃/C₃N/BC₃) is capable of the effective stretching and confinement of three types of intrinsically disordered proteins (IDPs), including amyloid-β (1–42), polyglutamine (Q42), and α-Synuclein (61–95). Our molecular dynamics simulations demonstrate that the protein molecules interact more strongly with the C₃N stripe than the BC₃ one, which leads to their capture, elongation, and confinement along the center C₃N stripe of the heterostructure. The conformational fluctuations of IDPs are substantially reduced after being stretched. This design may serve as a platform for single-molecule protein analysis with reduced thermal noise.

Dynamics-Based Peptide-MHC Binding Optimization by a Convolutional Variational Autoencoder: A Use-Case Model for CASTELO

An unsolved challenge in the development of antigen-specific immunotherapies is determining the optimal antigens to target. Comprehension of antigen-major histocompatibility complex (MHC) binding is paramount toward achieving this goal. Here, we apply CASTELO, a combined machine learning-molecular dynamics (ML-MD) approach, to identify per-residue antigen binding contributions and then design novel antigens of increased MHC-II binding affinity for a type 1 diabetes-implicated system. We build upon a small-molecule lead optimization algorithm by training a convolutional variational autoencoder (CVAE) on MD trajectories of 48 different systems across four antigens and four HLA serotypes. We develop several new machine learning metrics including a structure-based anchor residue classification model as well as cluster comparison scores. ML-MD predictions agree well with experimental binding results and free energy perturbation-predicted binding affinities. Moreover, ML-MD metrics are independent of traditional MD stability metrics such as contact area and root-mean-square fluctuations (RMSF), which do not reflect binding affinity data. Our work supports the role of structure-based deep learning techniques in antigen-specific immunotherapy design.

Comprehensive B-Cell Immune Repertoire Analysis of Anti-NMDAR Encephalitis and Anti-LGI1 Encephalitis

Anti-N-methyl-D-aspartate receptor encephalitis (anti-NMDARE) and anti-leucine-rich glioma-inactivated 1 encephalitis (anti-LGI1E) are the two most common types of antibody-mediated autoimmune encephalitis. We performed a comprehensive analysis of the B-cell immune repertoire in patients with anti-NMDARE (n = 7) and anti-LGI1E (n = 10) and healthy controls (n = 4). The results revealed the presence of many common clones between patients with these two types of autoimmune encephalitis, which were mostly class-switched. Additionally, many differences were found among the anti-NMDARE, anti-LGI1E, and healthy control groups, including the diversity of the B-cell immune repertoire and gene usage preference. These findings suggest that the same adaptive immune responses occur in patients with anti-NMDARE and anti-LGI1E, which deserves further exploration.

Insights From Single Cell RNA Sequencing Into the Immunology of Type 1 Diabetes- Cell Phenotypes and Antigen Specificity

In the past few years, huge advances have been made in techniques to analyse cells at an individual level using RNA sequencing, and many of these have precipitated exciting discoveries in the immunology of type 1 diabetes (T1D). This review will cover the first papers to use scRNAseq to characterise human lymphocyte phenotypes in T1D in the peripheral blood, pancreatic lymph nodes and islets. These have revealed specific genes such as IL-32 that are differentially expressed in islet -specific T cells in T1D. scRNAseq has also revealed wider gene expression patterns that are involved in T1D and can predict its development even predating autoantibody production. Single cell sequencing of TCRs has revealed V genes and CDR3 motifs that are commonly used to target islet autoantigens, although truly public TCRs remain elusive. Little is known about BCR repertoires in T1D, but scRNAseq approaches have revealed that insulin binding BCRs commonly use specific J genes, share motifs between donors and frequently demonstrate poly-reactivity. This review will also summarise new developments in scRNAseq technology, the insights they have given into other diseases and how they could be leveraged to advance research in the type 1 diabetes field to identify novel biomarkers and targets for immunotherapy.

A Two-Stage Flow Cytometry Strategy to Distinguish Single Cells from Doublets in Heterogeneous Cell Mixtures and Improve Cell Cluster Identification: Application to Human Monocyte Subpopulations

Flow cytometry is a powerful method, widely used to identify cell types present in tissues, to describe their phenotypes, and to purify cells for functional analyses. As a single cell technique, flow cytometry relies on identifying and excluding cell doublets and aggregates present in samples in the initial gating steps. This identification is based on detection of events generating electrical pulses falling outside of linear variations of pulse height, width, and area in a singlet population with increasing cell sizes. In heterogeneous cell mixtures, however, with cell types varying extensively in size and granularity, exclusion of doublets has the risk of removing single cells that co-localize with doublets of another cell type. This is particularly the case when doublets of a smaller cell type overlap with large cells of a distinct, larger cell type. Here, we describe a gating method to reduce this risk. In this protocol, initial gating steps aim to segregate cells according to physical characteristics (such as size and granularity) and gene expression properties in order to obtain more homogeneous cell clusters. Doublet exclusion is then performed separately in each cluster, minimizing the risk of confusion between single cells and doublets. To illustrate this protocol, human blood monocytes are separated and analyzed. By implementing this protocol, we were able to reveal the existence of a population of large monocytes previously unrecognized using conventional gating strategies. In subsequent functional assays, we have shown that this novel population exhibits unique inflammatory responses, highlighting the need and pertinence of this approach to identify and characterize infrequent-yet functionally relevant-cell populations present in complex cell mixtures. © 2021 Wiley Periodicals LLC. Basic Protocol: Distinguishing single cells from doublets in heterogeneous cell mixtures by flow cytometry.

The methods and advances of adaptive immune receptors repertoire sequencing

The adaptive immune response is a powerful tool, capable of recognizing, binding to, and neutralizing a vast number of internal and external threats via T or B lymphatic receptors with widespread sets of antigen specificities. The emergence of high-throughput sequencing technology and bioinformatics provides opportunities for research in the fields of life sciences and medicine. The analysis and annotation for immune repertoire data can reveal biologically meaningful information, including immune prediction, target antigens, and effective evaluation. Continuous improvements of the immunological repertoire sequencing methods and analysis tools will help to minimize the experimental and calculation errors and realize the immunological information to meet the clinical requirements. That said, the clinical application of adaptive immune repertoire sequencing requires appropriate experimental methods and standard analytical tools. At the population cell level, we can acquire the overview of cell groups, but the information about a single cell is not obtained accurately. The information that is ignored may be crucial for understanding the heterogeneity of each cell, gene expression and drug response. The combination of high-throughput sequencing and single-cell technology allows us to obtain single-cell information with low-cost and high-throughput. In this review, we summarized the current methods and progress in this area.

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.

Toward Guided Mutagenesis: Gaussian Process Regression Predicts MHC Class II Antigen Mutant Binding

Antigen-specific immunotherapies (ASI) require successful loading and presentation of antigen peptides into the major histocompatibility complex (MHC) binding cleft. One route of ASI design is to mutate native antigens for either stronger or weaker binding interaction to MHC. Exploring all possible mutations is costly both experimentally and computationally. To reduce experimental and computational expense, here we investigate the minimal amount of prior data required to accurately predict the relative binding affinity of point mutations for peptide-MHC class II (pMHCII) binding. Using data from different residue subsets, we interpolate pMHCII mutant binding affinities by Gaussian process (GP) regression of residue volume and hydrophobicity. We apply GP regression to an experimental data set from the Immune Epitope Database, and theoretical data sets from NetMHCIIpan and Free Energy Perturbation calculations. We find that GP regression can predict binding affinities of nine neutral residues from a six-residue subset with an average R² coefficient of determination value of 0.62 ± 0.04 (±95% CI), average error of 0.09 ± 0.01 kcal/mol (±95% CI), and with an receiver operating characteristic (ROC) AUC value of 0.92 for binary classification of enhanced or diminished binding affinity. Similarly, metrics increase to an R² value of 0.69 ± 0.04, average error of 0.07 ± 0.01 kcal/mol, and an ROC AUC value of 0.94 for predicting seven neutral residues from an eight-residue subset. Our work finds that prediction is most accurate for neutral residues at anchor residue sites without register shift. This work holds relevance to predicting pMHCII binding and accelerating ASI design.