Computational Strategies for Dissecting the High-Dimensional Complexity of Adaptive Immune Repertoires

A comprehensive evaluation of diversity measures for TCR repertoire profiling

BACKGROUND

T cells play a crucial role in adaptive immunity, as they monitor internal and external immunogenic signals through their specific receptors (TCRs). Using high-throughput sequencing, one can assess TCR repertoire in various clinical settings and describe it quantitatively by calculating a diversity index. Multiple diversity indices that capture the richness of TCRs and the evenness of their distribution have been proposed in the literature; however, there is no consensus on gold-standard measures and interpretation of each index is complex. Our goal was to examine the performance characteristics of 12 commonly used diversity indices in simulated and real-world data.

RESULTS

Simulated data were generated to evaluate how data richness and evenness affect index values using three nonparametric models. Fourteen real-world TCR datasets were obtained to examine differences in indices by analysis protocols and test their robustness to subsampling. Pielou, Basharin, d50, and Gini primarily describe evenness and highly correlate with one another. They are best suited for measuring the representation of TCR clones. Richness is best captured by S index, next Chao1 and ACE which also consider information on evenness. Shannon, Inv.Simspon, D3, D4, and Gini.Simpson measure richness and increasingly more information on evenness. More skewed TCR distributions provided more stable results in subsampling. Gini-Simpson, Pielou, and Basharin were the most robust in both simulated and experimental data.

CONCLUSIONS

Our results could guide investigators to select the best diversity index for their particular experimental question and draw attention to factors that can influence the accuracy and reproducibility of results.

Antigen-driven T cell responses in rheumatic diseases: insights from T cell receptor repertoire studies

Advances in T cell receptor (TCR) profiling techniques have substantially improved our ability to investigate T cell responses to antigens that are presented on HLA class I and class II molecules and associations between autoimmune T cells and rheumatic diseases. Early-stage studies in axial spondyloarthritis (axSpA) identified disease-associated T cell clonotypes, benefiting from the relative genetic homogeneity of the disease. However, both the genetic and the T cell immunological landscape are more complex in other rheumatic diseases. The diversity or redundancy in the TCR repertoire, epitope spreading over disease duration, genetic heterogeneity of HLA genes or other loci, and the diversity of epitopes contributing to disease pathogenesis and persistent inflammation are all likely to contribute to this complexity. TCR profiling holds promise for identifying key antigenic drivers and phenotypic T cell states that sustain autoimmunity in rheumatic diseases. Here, we review key findings from TCR repertoire studies in axSpA and other chronic inflammatory rheumatic diseases including psoriatic arthritis, rheumatoid arthritis, systemic lupus erythematosus and Sjögren syndrome. We explore how TCR profiling technologies, if applied to better controlled studies focused on early disease stages and genetically homogeneous subsets, can facilitate disease monitoring and the development of therapeutics targeting autoimmune T cells, their cognate antigens, or their underlying biology.

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

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

Reading the repertoire: Progress in adaptive immune receptor analysis using machine learning

The adaptive immune system holds invaluable information on past and present immune responses in the form of B and T cell receptor sequences, but we are limited in our ability to decode this information. Machine learning approaches are under active investigation for a range of tasks relevant to understanding and manipulating the adaptive immune receptor repertoire, including matching receptors to the antigens they bind, generating antibodies or T cell receptors for use as therapeutics, and diagnosing disease based on patient repertoires. Progress on these tasks has the potential to substantially improve the development of vaccines, therapeutics, and diagnostics, as well as advance our understanding of fundamental immunological principles. We outline key challenges for the field, highlighting the need for software benchmarking, targeted large-scale data generation, and coordinated research efforts.

Structural trends in antibody-antigen binding interfaces: a computational analysis of 1833 experimentally determined 3D structures

Antibodies are attractive therapeutic candidates due to their ability to bind cognate antigens with high affinity and specificity. Still, the underlying molecular rules governing the antibody-antigen interface remain poorly understood, making in silico antibody design inherently difficult and keeping the discovery and design of novel antibodies a costly and laborious process. This study investigates the characteristics of antibody-antigen binding interfaces through a computational analysis of more than 850,000 atom-atom contacts from the largest reported set of antibody-antigen complexes with 1833 nonredundant, experimentally determined structures. The analysis compares binding characteristics of conventional antibodies and single-domain antibodies (sdAbs) targeting both protein- and peptide antigens. We find clear patterns in the number antibody-antigen contacts and amino acid frequencies in the paratope. The direct comparison of sdAbs and conventional antibodies helps elucidate the mechanisms employed by sdAbs to compensate for their smaller size and the fact that they harbor only half the number of complementarity-determining regions compared to conventional antibodies. Furthermore, we pinpoint antibody interface hotspot residues that are often found at the binding interface and the amino acid frequencies at these positions. These findings have direct potential applications in antibody engineering and the design of improved antibody libraries.

tcrBLOSUM: an amino acid substitution matrix for sensitive alignment of distant epitope-specific TCRs

Deciphering the specificity of T-cell receptor (TCR) repertoires is crucial for monitoring adaptive immune responses and developing targeted immunotherapies and vaccines. To elucidate the specificity of previously unseen TCRs, many methods employ the BLOSUM62 matrix to find TCRs with similar amino acid (AA) sequences. However, while BLOSUM62 reflects the AA substitutions within conserved regions of proteins with similar functions, the remarkable diversity of TCRs means that both TCRs with similar and dissimilar sequences can bind the same epitope. Therefore, reliance on BLOSUM62 may bias detection towards epitope-specific TCRs with similar biochemical properties, overlooking those with more diverse AA compositions. In this study, we introduce tcrBLOSUMa and tcrBLOSUMb, specialized AA substitution matrices for CDR3 alpha and CDR3 beta TCR chains, respectively. The matrices reflect AA frequencies and variations occurring within TCRs that bind the same epitope, revealing that both CDR3 alpha and CDR3 beta display tolerance to a wide range of AA substitutions and differ noticeably from the standard BLOSUM62. By accurately aligning distant TCRs employing tcrBLOSUMb, we were able to improve clustering performance and capture a large number of epitope-specific TCRs with diverse AA compositions and physicochemical profiles overlooked by BLOSUM62. Utilizing both the general BLOSUM62 and specialized tcrBLOSUM matrices in existing computational tools will broaden the range of TCRs that can be associated with their cognate epitopes, thereby enhancing TCR repertoire analysis.

Machine learning for precision diagnostics of autoimmunity

Early and accurate diagnosis is crucial to prevent disease development and define therapeutic strategies. Due to predominantly unspecific symptoms, diagnosis of autoimmune diseases (AID) is notoriously challenging. Clinical decision support systems (CDSS) are a promising method with the potential to enhance and expedite precise diagnostics by physicians. However, due to the difficulties of integrating and encoding multi-omics data with clinical values, as well as a lack of standardization, such systems are often limited to certain data types. Accordingly, even sophisticated data models fall short when making accurate disease diagnoses and presenting data analyses in a user-friendly form. Therefore, the integration of various data types is not only an opportunity but also a competitive advantage for research and industry. We have developed an integration pipeline to enable the use of machine learning for patient classification based on multi-omics data in combination with clinical values and laboratory results. The application of our framework resulted in up to 96% prediction accuracy of autoimmune diseases with machine learning models. Our results deliver insights into autoimmune disease research and have the potential to be adapted for applications across disease conditions.

Unraveling the Immune Web: Advances in SMI Capsular Fibrosis from Molecular Insights to Preclinical Breakthroughs

Breast implant surgery has evolved significantly, yet challenges such as capsular contracture remain a persistent concern. This review presents an in-depth analysis of recent advancements in understanding the immune mechanisms and clinical implications associated with silicone mammary implants (SMIs). The article systematically examines the complex interplay between immune responses and capsular fibrosis, emphasizing the pathophysiological mechanisms of inflammation in the etiology of this fibrotic response. It discusses innovations in biomaterial science, including the development of novel anti-biofilm coatings and immunomodulatory surfaces designed to enhance implant integration and minimize complications. Emphasis is placed on personalized risk assessment strategies, leveraging molecular insights to tailor interventions and improve patient outcomes. Emerging therapeutic targets, advancements in surgical techniques, and the refinement of post-operative care are also explored. Despite notable progress, challenges such as the variability in immune responses, the long-term efficacy of new interventions, and ethical considerations remain. Future research directions are identified, focusing on personalized medicine, advanced biomaterials, and bridging preclinical findings with clinical applications. As we advance from bench to bedside, this review illuminates the path forward, where interdisciplinary collaboration and continued inquiry weave together to enhance the art and science of breast implant surgery, transforming patient care into a realm of precision and excellence.

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

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

mRNA-based influenza vaccine expands breadth of B cell response in humans

Eliciting broad and durable antibody responses against rapidly evolving pathogens like influenza viruses remains a formidable challenge1,2. The germinal center (GC) reaction enables the immune system to generate broad, high-affinity, and durable antibody responses to vaccination3-5. mRNA-based severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccines induce persistent GC B cell responses in humans6-9. Whether an mRNA-based influenza vaccine could induce a superior GC response in humans compared to the conventional inactivated influenza virus vaccine remains unclear. We assessed B cell responses in peripheral blood and draining lymph nodes in cohorts receiving the inactivated or mRNA-based quadrivalent seasonal influenza vaccine. Participants receiving the mRNA-based vaccine produced more robust plasmablast responses and higher antibody titers to H1N1 and H3N2 influenza A viruses and comparable antibody titers against influenza B virus strains. Importantly, mRNA-based vaccination stimulated robust recall B cell responses characterized by sustained GC reactions that lasted at least 26 weeks post-vaccination in three of six participants analyzed. In addition to promoting the maturation of responding B cell clones, these sustained GC reactions resulted in enhanced engagement of low-frequency pre-existing memory B cells, expanding the landscape of vaccine-elicited B cell clones. This translated to expansion of the serological repertoire and increased breadth of serum antibody responses. These findings reveal an important role for the induction of persistent GC responses to influenza vaccination in humans to broaden the repertoire of vaccine-induced antibodies.

B cell tolerance and autoimmunity: Lessons from repertoires

Adaptive immune cell function is regulated by a highly diverse receptor recombined from variable germline-encoded segments that can recognize an almost unlimited array of epitopes. While this diversity enables the recognition of any pathogen, it also poses a risk of self-recognition, leading to autoimmunity. Many layers of regulation are present during both the generation and activation of B cells to prevent this phenomenon, although they are evidently imperfect. In recent years, our ability to analyze immune repertoires at scale has drastically increased, both through advances in sequencing and single-cell analyses. Here, we review the current knowledge on B cell repertoire analyses, focusing on their implication for autoimmunity. These studies demonstrate that a failure of tolerance occurs at multiple independent checkpoints in different autoimmune contexts, particularly during B cell maturation, plasmablast differentiation, and within germinal centers. These failures are marked by distinct repertoire features that may be used to identify disease- or patient-specific therapeutic approaches.

The clinical applications of immunosequencing

Technological advances in high-throughput sequencing have opened the door for the interrogation of adaptive immune responses at unprecedented scale. It is now possible to determine the sequences of antibodies or T-cell receptors produced by individual B and T cells in a sample. This capability, termed immunosequencing, has transformed the study of both infectious and non-infectious diseases by allowing the tracking of dynamic changes in B and T cell clonal populations over time. This has improved our understanding of the pathology of cancers, autoimmune diseases, and infectious diseases. However, to date there has been only limited clinical adoption of the technology. Advances over the last decade and on the horizon that reduce costs and improve interpretability could enable widespread clinical use. Many clinical applications have been proposed and, while most are still undergoing research and development, some methods relying on immunosequencing data have been implemented, the most widespread of which is the detection of measurable residual disease. Here, we review the diagnostic, prognostic, and therapeutic applications of immunosequencing for both infectious and non-infectious diseases.

Benchmarking and integrating human B-cell receptor genomic and antibody proteomic profiling

Immunoglobulins (Ig), which exist either as B-cell receptors (BCR) on the surface of B cells or as antibodies when secreted, play a key role in the recognition and response to antigenic threats. The capability to jointly characterize the BCR and antibody repertoire is crucial for understanding human adaptive immunity. From peripheral blood, bulk BCR sequencing (bulkBCR-seq) currently provides the highest sampling depth, single-cell BCR sequencing (scBCR-seq) allows for paired chain characterization, and antibody peptide sequencing by tandem mass spectrometry (Ab-seq) provides information on the composition of secreted antibodies in the serum. Yet, it has not been benchmarked to what extent the datasets generated by these three technologies overlap and complement each other. To address this question, we isolated peripheral blood B cells from healthy human donors and sequenced BCRs at bulk and single-cell levels, in addition to utilizing publicly available sequencing data. Integrated analysis was performed on these datasets, resolved by replicates and across individuals. Simultaneously, serum antibodies were isolated, digested with multiple proteases, and analyzed with Ab-seq. Systems immunology analysis showed high concordance in repertoire features between bulk and scBCR-seq within individuals, especially when replicates were utilized. In addition, Ab-seq identified clonotype-specific peptides using both bulk and scBCR-seq library references, demonstrating the feasibility of combining scBCR-seq and Ab-seq for reconstructing paired-chain Ig sequences from the serum antibody repertoire. Collectively, our work serves as a proof-of-principle for combining bulk sequencing, single-cell sequencing, and mass spectrometry as complementary methods towards capturing humoral immunity in its entirety.

Shape of the art: TCR-repertoire after allogeneic hematopoietic cell transplantation

The human adaptive immune repertoire is characterized by specificity and diversity to provide immunity against past and future tasks. Such tasks are mainly infections but also malignant transformations of cells. With its multiple lines of defense, the human immune system contains both, rapid reaction forces and the potential to capture, disassemble and analyze strange structures in order to teach the adaptive immune system and mount a specific immune response. Prevention and mitigation of autoimmunity is of equal importance. In the context of allogeneic hematopoietic cell transplantation (HCT) specific challenges exist with the transfer of cells from the adapted donor immune system to the immunosuppressed recipient. Those challenges are immunogenetic disparity between donor and host, reconstitution of immunity early after HCT by expansion of mature immune effector cells, and impaired thymic function, if the recipient is an adult (as it is the case in most HCTs). The possibility to characterize the adaptive immune repertoire by massively parallel sequencing of T-cell receptor gene rearrangements allows for a much more detailed characterization of the T-cell repertoire. In addition, high-dimensional characterization of immune effector cells based on their immunophenotype and single cell RNA sequencing allow for much deeper insights in adaptive immune responses. We here review, existing - still incomplete - information on immune reconstitution after allogeneic HCT. Building on the technological advances much deeper insights into immune recovery after HCT and adaptive immune responses and can be expected in the coming years.

Enhancing comparative T cell receptor repertoire analysis in small biological samples through pooling homologous cell samples from multiple mice

Accurate characterization and comparison of T cell receptor (TCR) repertoires from small biological samples present significant challenges. The main challenge is the low material input, which compromises the quality of bulk sequencing and hinders the recovery of sufficient TCR sequences for robust analyses. We aimed to address this limitation by implementing a strategic approach to pool homologous biological samples. Our findings demonstrate that such pooling indeed enhances the TCR repertoire coverage, particularly for cell subsets of constrained sizes, and enables accurate comparisons of TCR repertoires at different levels of complexity across T cell subsets with different sizes. This methodology holds promise for advancing our understanding of T cell repertoires in scenarios where sample size constraints are a prevailing concern.

Decoding Immuno-Competence: A Novel Analysis of Complete Blood Cell Count Data in COVID-19 Outcomes

BACKGROUND

While 'immuno-competence' is a well-known term, it lacks an operational definition. To address this omission, this study explored whether the temporal and structured data of the complete blood cell count (CBC) can rapidly estimate immuno-competence. To this end, one or more ratios that included data on all monocytes, lymphocytes and neutrophils were investigated.

MATERIALS AND METHODS

Longitudinal CBC data collected from 101 COVID-19 patients (291 observations) were analyzed. Dynamics were estimated with several approaches, which included non-structured (the classic CBC format) and structured data. Structured data were assessed as complex ratios that capture multicellular interactions among leukocytes. In comparing survivors with non-survivors, the hypothesis that immuno-competence may exhibit feedback-like (oscillatory or cyclic) responses was tested.

RESULTS

While non-structured data did not distinguish survivors from non-survivors, structured data revealed immunological and statistical differences between outcomes: while survivors exhibited oscillatory data patterns, non-survivors did not. In survivors, many variables (including IL-6, hemoglobin and several complex indicators) showed values above or below the levels observed on day 1 of the hospitalization period, displaying L-shaped data distributions (positive kurtosis). In contrast, non-survivors did not exhibit kurtosis. Three immunologically defined data subsets included only survivors. Because information was based on visual patterns generated in real time, this method can, potentially, provide information rapidly.

DISCUSSION

The hypothesis that immuno-competence expresses feedback-like loops when immunological data are structured was not rejected. This function seemed to be impaired in immuno-suppressed individuals. While this method rapidly informs, it is only a guide that, to be confirmed, requires additional tests. Despite this limitation, the fact that three protective (survival-associated) immunological data subsets were observed since day 1 supports many clinical decisions, including the early and personalized prognosis and identification of targets that immunomodulatory therapies could pursue. Because it extracts more information from the same data, structured data may replace the century-old format of the CBC.

Prediction of polyspecificity from antibody sequence data by machine learning

Antibodies are generated with great diversity in nature resulting in a set of molecules, each optimized to bind a specific target. Taking advantage of their diversity and specificity, antibodies make up for a large part of recently developed biologic drugs. For therapeutic use antibodies need to fulfill several criteria to be safe and efficient. Polyspecific antibodies can bind structurally unrelated molecules in addition to their main target, which can lead to side effects and decreased efficacy in a therapeutic setting, for example via reduction of effective drug levels. Therefore, we created a neural-network-based model to predict polyspecificity of antibodies using the heavy chain variable region sequence as input. We devised a strategy for enriching antibodies from an immunization campaign either for antigen-specific or polyspecific binding properties, followed by generation of a large sequencing data set for training and cross-validation of the model. We identified important physico-chemical features influencing polyspecificity by investigating the behaviour of this model. This work is a machine-learning-based approach to polyspecificity prediction and, besides increasing our understanding of polyspecificity, it might contribute to therapeutic antibody development.

Disease diagnostics using machine learning of immune receptors

Clinical diagnosis typically incorporates physical examination, patient history, and various laboratory tests and imaging studies, but makes limited use of the human system's own record of antigen exposures encoded by receptors on B cells and T cells. We analyzed immune receptor datasets from 593 individuals to develop MAchine Learning for Immunological Diagnosis (Mal-ID) , an interpretive framework to screen for multiple illnesses simultaneously or precisely test for one condition. This approach detects specific infections, autoimmune disorders, vaccine responses, and disease severity differences. Human-interpretable features of the model recapitulate known immune responses to SARS-CoV-2, Influenza, and HIV, highlight antigen-specific receptors, and reveal distinct characteristics of Systemic Lupus Erythematosus and Type-1 Diabetes autoreactivity. This analysis framework has broad potential for scientific and clinical interpretation of human immune responses.

Dynamics of The Γδtcr Repertoires During The Dedifferentiation Process and Pilot Implications for Immunotherapy of Thyroid Cancer

γδ T cells are evolutionarily conserved T lymphocytes that manifest unique antitumor efficacy independent of tumor mutation burden (TMB) and conventional human leukocyte antigen (HLA) recognition. However, the dynamic changes in their T cell receptor (TCR) repertoire during cancer progression and treatment courses remain unclear. Here, a comprehensive characterization of γδTCR repertoires are performed in thyroid cancers with divergent differentiation states through cross-sectional studies. The findings revealed a significant correlation between the differentiation states and TCR repertoire diversity. Notably, highly expanded clones are prominently enriched in γδ T cell compartment of dedifferentiated patients. Moreover, by longitudinal investigations of the γδ T cell response to various antitumor therapies, it is found that the emergence and expansion of the Vδ2neg subset may be potentially associated with favorable clinical outcomes after post-radiotherapeutic immunotherapy. These findings are further validated at single-cell resolution in both advanced thyroid cancer patients and a murine model, underlining the importance of further investigations into the role of γδTCR in cancer immunity and therapeutic strategies.

Extraction of the CDRH3 sequence of the mouse antibody repertoire selected upon influenza virus infection by subtraction of the background antibody repertoire

Historically, antibody reactivity to pathogens and vaccine antigens has been evaluated using serological measurements of antigen-specific antibodies. However, it is difficult to evaluate all antibodies that contribute to various functions in a single assay, such as the measurement of the neutralizing antibody titer. Bulk antibody repertoire analysis using next-generation sequencing is a comprehensive method for analyzing the overall antibody response; however, it is unreliable for estimating antigen-specific antibodies due to individual variation. To address this issue, we propose a method to subtract the background signal from the repertoire of data of interest. In this study, we analyzed changes in antibody diversity and inferred the heavy-chain complementarity-determining region 3 (CDRH3) sequences of antibody clones that were selected upon influenza virus infection in a mouse model using bulk repertoire analysis. A decrease in the diversity of the antibody repertoire was observed upon viral infection, along with an increase in neutralizing antibody titers. Using kernel density estimation of sequences in a high-dimensional sequence space with background signal subtraction, we identified several clusters of CDRH3 sequences induced upon influenza virus infection. Most of these repertoires were detected more frequently in infected mice than in uninfected control mice, suggesting that infection-specific antibody sequences can be extracted using this method. Such an accurate extraction of antigen- or infection-specific repertoire information will be a useful tool for vaccine evaluation in the future.

IMPORTANCE

As specific interactions between antigens and cell-surface antibodies trigger the proliferation of B-cell clones, the frequency of each antibody sequence in the samples reflects the size of each clonal population. Nevertheless, it is extremely difficult to extract antigen-specific antibody sequences from the comprehensive bulk antibody sequences obtained from blood samples due to repertoire bias influenced by exposure to dietary antigens and other infectious agents. This issue can be addressed by subtracting the background noise from the post-immunization or post-infection repertoire data. In the present study, we propose a method to quantify repertoire data from comprehensive repertoire data. This method allowed subtraction of the background repertoire, resulting in more accurate extraction of expanded antibody repertoires upon influenza virus infection. This accurate extraction of antigen- or infection-specific repertoire information is a useful tool for vaccine evaluation.

Forum on immune digital twins: a meeting report

Medical digital twins are computational models of human biology relevant to a given medical condition, which are tailored to an individual patient, thereby predicting the course of disease and individualized treatments, an important goal of personalized medicine. The immune system, which has a central role in many diseases, is highly heterogeneous between individuals, and thus poses a major challenge for this technology. In February 2023, an international group of experts convened for two days to discuss these challenges related to immune digital twins. The group consisted of clinicians, immunologists, biologists, and mathematical modelers, representative of the interdisciplinary nature of medical digital twin development. A video recording of the entire event is available. This paper presents a synopsis of the discussions, brief descriptions of ongoing digital twin projects at different stages of progress. It also proposes a 5-year action plan for further developing this technology. The main recommendations are to identify and pursue a small number of promising use cases, to develop stimulation-specific assays of immune function in a clinical setting, and to develop a database of existing computational immune models, as well as advanced modeling technology and infrastructure.

Microfluidic-assisted single-cell RNA sequencing facilitates the development of neutralizing monoclonal antibodies against SARS-CoV-2

As a class of antibodies that specifically bind to a virus and block its entry, neutralizing monoclonal antibodies (neutralizing mAbs) have been recognized as a top choice for combating COVID-19 due to their high specificity and efficacy in treating serious infections. Although conventional approaches for neutralizing mAb development have been optimized for decades, there is an urgent need for workflows with higher efficiency due to time-sensitive concerns, including the high mutation rate of SARS-CoV-2. One promising approach is the identification of neutralizing mAb candidates via single-cell RNA sequencing (RNA-seq), as each B cell has a unique transcript sequence corresponding to its secreted antibody. The state-of-the-art high-throughput single-cell sequencing technologies, which have been greatly facilitated by advances in microfluidics, have greatly accelerated the process of neutralizing mAb development. Here, we provide an overview of the general procedures for high-throughput single-cell RNA-seq enabled by breakthroughs in droplet microfluidics, introduce revolutionary approaches that combine single-cell RNA-seq to facilitate the development of neutralizing mAbs against SARS-CoV-2, and outline future steps that need to be taken to further improve development strategies for effective treatments against infectious diseases.

Ovarian cancer is detectable from peripheral blood using machine learning over T-cell receptor repertoires

The extraordinary diversity of T cells and B cells is critical for body maintenance. This diversity has an important role in protecting against tumor formation. In humans, the T-cell receptor (TCR) repertoire is generated through a striking stochastic process called V(D)J recombination, in which different gene segments are assembled and modified, leading to extensive variety. In ovarian cancer (OC), an unfortunate 80% of cases are detected late, leading to poor survival outcomes. However, when detected early, approximately 94% of patients live longer than 5 years after diagnosis. Thus, early detection is critical for patient survival. To determine whether the TCR repertoire obtained from peripheral blood is associated with tumor status, we collected blood samples from 85 women with or without OC and obtained TCR information. We then used machine learning to learn the characteristics of samples and to finally predict, over a set of unseen samples, whether the person is with or without OC. We successfully stratified the two groups, thereby associating the peripheral blood TCR repertoire with the formation of OC tumors. A careful study of the origin of the set of T cells most informative for the signature indicated the involvement of a specific invariant natural killer T (iNKT) clone and a specific mucosal-associated invariant T (MAIT) clone. Our findings here support the proposition that tumor-relevant signal is maintained by the immune system and is coded in the T-cell repertoire available in peripheral blood. It is also possible that the immune system detects tumors early enough for repertoire technologies to inform us near the beginning of tumor formation. Although such detection is made by the immune system, we might be able to identify it, using repertoire data from peripheral blood, to offer a pragmatic way to search for early signs of cancer with minimal patient burden, possibly with enhanced sensitivity.

The dengue-specific immune response and antibody identification with machine learning

Dengue virus poses a serious threat to global health and there is no specific therapeutic for it. Broadly neutralizing antibodies recognizing all serotypes may be an effective treatment. High-throughput adaptive immune receptor repertoire sequencing (AIRR-seq) and bioinformatic analysis enable in-depth understanding of the B-cell immune response. Here, we investigate the dengue antibody response with these technologies and apply machine learning to identify rare and underrepresented broadly neutralizing antibody sequences. Dengue immunization elicited the following signatures on the antibody repertoire: (i) an increase of CDR3 and germline gene diversity; (ii) a change in the antibody repertoire architecture by eliciting power-law network distributions and CDR3 enrichment in polar amino acids; (iii) an increase in the expression of JNK/Fos transcription factors and ribosomal proteins. Furthermore, we demonstrate the applicability of computational methods and machine learning to AIRR-seq datasets for neutralizing antibody candidate sequence identification. Antibody expression and functional assays have validated the obtained results.

Systems immunology spanning tumors, lymph nodes, and periphery

The immune system defines a complex network of tissues and cell types that orchestrate responses across the body in a dynamic manner. The local and systemic interactions between immune and cancer cells contribute to disease progression. Lymphocytes are activated in lymph nodes, traffic through the periphery, and impact cancer progression through their interactions with tumor cells. As a result, therapeutic response and resistance are mediated across tissues, and a comprehensive understanding of lymphocyte dynamics requires a systems-level approach. In this review, we highlight experimental and computational methods that can leverage the study of leukocyte trafficking through an immunomics lens and reveal how adaptive immunity shapes cancer.

Human thymic putative CD8αα precursors exhibit a biased TCR repertoire in single cell AIRR-seq

Thymic T cell development comprises T cell receptor (TCR) recombination and assessment of TCR avidity towards self-peptide-MHC complexes presented by antigen-presenting cells. Self-reactivity may lead to negative selection, or to agonist selection and differentiation into unconventional lineages such as regulatory T cells and CD8[Formula: see text] T cells. To explore the effect of the adaptive immune receptor repertoire on thymocyte developmental decisions, we performed single cell adaptive immune receptor repertoire sequencing (scAIRR-seq) of thymocytes from human young paediatric thymi and blood. Thymic PDCD1+ cells, a putative CD8[Formula: see text] T cell precursor population, exhibited several TCR features previously associated with thymic and peripheral ZNF683+ CD8[Formula: see text] T cells, including enrichment of large and positively charged complementarity-determining region 3 (CDR3) amino acids. Thus, the TCR repertoire may partially explain the decision between conventional vs. agonist selected thymocyte differentiation, an aspect of importance for the development of therapies for patients with immune-mediated diseases.

Identification of TCR repertoire patterns linked with anti-cancer immunotherapy

The highly diverse T cell receptor (TCR) repertoire is a crucial component of the adaptive immune system that aids in the protection against a wide variety of pathogens. This TCR repertoire, comprising the collection of all TCRs in an individual, is a valuable source of information on both recent and ongoing T cell activation. Cancer cells, like pathogens, have the ability to trigger an adaptive immune response. However, because cancer cells use a variety of strategies to escape immune responses, this is often insufficient to completely eradicate them. As a result, immunotherapy is a promising treatment option for cancer patients. This treatment is expected to increase T cell activation and subsequently alter the TCR repertoire composition in these patients. Monitoring TCR repertoires before and after immunotherapy can therefore provide additional insight into T cell responses and might identify cancer-associated TCR sequences. Here we present a computational strategy to identify those changes in the TCR repertoire that occur after treatment with immunotherapy. Since this method allows the identification of TCR patterns that might be treatment-associated, it can help future research by revealing those patterns that are related with response. This TCR analysis workflow is illustrated using public data from three different cancer patients who received anti-PD-1 treatment.

NAIR: Network Analysis of Immune Repertoire

T cells represent a crucial component of the adaptive immune system and mediate anti-tumoral immunity as well as protection against infections, including respiratory viruses such as SARS-CoV-2. Next-generation sequencing of the T-cell receptors (TCRs) can be used to profile the T-cell repertoire. We developed a customized pipeline for Network Analysis of Immune Repertoire (NAIR) with advanced statistical methods to characterize and investigate changes in the landscape of TCR sequences. We first performed network analysis on the TCR sequence data based on sequence similarity. We then quantified the repertoire network by network properties and correlated it with clinical outcomes of interest. In addition, we identified (1) disease-specific/associated clusters and (2) shared clusters across samples based on our customized search algorithms and assessed their relationship with clinical outcomes such as recovery from COVID-19 infection. Furthermore, to identify disease-specific TCRs, we introduced a new metric that incorporates the clonal generation probability and the clonal abundance by using the Bayes factor to filter out the false positives. TCR-seq data from COVID-19 subjects and healthy donors were used to illustrate that the proposed approach to analyzing the network architecture of the immune repertoire can reveal potential disease-specific TCRs responsible for the immune response to infection.

Dietary restriction mitigates the age-associated decline in mouse B cell receptor repertoire diversity

Aging impairs the capacity to respond to novel antigens, reducing immune protection against pathogens and vaccine efficacy. Dietary restriction (DR) extends life- and health span in diverse animals. However, little is known about the capacity of DR to combat the decline in immune function. Here, we study the changes in B cell receptor (BCR) repertoire during aging in DR and control mice. By sequencing the variable region of the BCR heavy chain in the spleen, we show that DR preserves diversity and attenuates the increase in clonal expansions throughout aging. Remarkably, mice starting DR in mid-life have repertoire diversity and clonal expansion rates indistinguishable from chronic DR mice. In contrast, in the intestine, these traits are unaffected by either age or DR. Reduced within-individual B cell repertoire diversity and increased clonal expansions are correlated with higher morbidity, suggesting a potential contribution of B cell repertoire dynamics to health during aging.

Updates on global epidemiology, risk and prognostic factors of gastric cancer

Gastric cancer (GC) is defined as the primary epithelial malignancy derived from the stomach, and it is a complicated and heterogeneous disease with multiple risk factors. Despite its overall declining trend of incidence and mortality in various countries over the past few decades, GC remains the fifth most common malignancy and the fourth leading cause of cancer-related death globally. Although the global burden of GC has shown a significant downward trend, it remains severe in certain areas, such as Asia. GC ranks third in incidence and mortality among all cancer types in China, and it accounts for nearly 44.0% and 48.6% of new GC cases and GC-related deaths in the world, respectively. The regional differences in GC incidence and mortality are obvious, and annual new cases and deaths are increasing rapidly in some developing regions. Therefore, early preventive and screening strategies for GC are urgently needed. The clinical efficacies of conventional treatments for GC are limited, and the developing understanding of GC pathogenesis has increased the demand for new therapeutic regimens, including immune checkpoint inhibitors, cell immunotherapy and cancer vaccines. The present review describes the epidemiology of GC worldwide, especially in China, summarizes its risk and prognostic factors, and focuses on novel immunotherapies to develop therapeutic strategies for the management of GC patients.

Entropic analysis of antigen-specific CDR3 domains identifies essential binding motifs shared by CDR3s with different antigen specificities

Antigen-specific T cell receptor (TCR) sequences can have prognostic, predictive, and therapeutic value, but decoding the specificity of TCR recognition remains challenging. Unlike DNA strands that base pair, TCRs bind to their targets with different orientations and different lengths, which complicates comparisons. We present scanning parametrized by normalized TCR length (SPAN-TCR) to analyze antigen-specific TCR CDR3 sequences and identify patterns driving TCR-pMHC specificity. Using entropic analysis, SPAN-TCR identifies 2-mer motifs that decrease the diversity (entropy) of CDR3s. These motifs are the most common patterns that can predict CDR3 composition, and we identify "essential" motifs that decrease entropy in the same CDR3 α or β chain containing the 2-mer, and "super-essential" motifs that decrease entropy in both chains. Molecular dynamics analysis further suggests that these motifs may play important roles in binding. We then employ SPAN-TCR to resolve similarities in TCR repertoires against different antigens using public databases of TCR sequences.

B cell M-CLL clones retain selection against replacement mutations in their immunoglobulin gene framework regions

Introduction

Chronic lymphocytic leukemia (CLL) is the most common adult leukemia, accounting for 30–40% of all adult leukemias. The dynamics of B-lymphocyte CLL clones with mutated immunoglobulin heavy chain variable region (IgHV) genes in their tumor (M-CLL) can be studied using mutational lineage trees.

Methods

Here, we used lineage tree-based analyses of somatic hypermutation (SHM) and selection in M-CLL clones, comparing the dominant (presumably malignant) clones of 15 CLL patients to their non-dominant (presumably normal) B cell clones, and to those of healthy control repertoires. This type of analysis, which was never previously published in CLL, yielded the following novel insights.

Results

CLL dominant clones undergo – or retain – more replacement mutations that alter amino acid properties such as charge or hydropathy. Although, as expected, CLL dominant clones undergo weaker selection for replacement mutations in the complementarity determining regions (CDRs) and against replacement mutations in the framework regions (FWRs) than non-dominant clones in the same patients or normal B cell clones in healthy controls, they surprisingly retain some of the latter selection in their FWRs. Finally, using machine learning, we show that even the non-dominant clones in CLL patients differ from healthy control clones in various features, most notably their expression of higher fractions of transition mutations.

Discussion

Overall, CLL seems to be characterized by significant loosening – but not a complete loss – of the selection forces operating on B cell clones, and possibly also by changes in SHM mechanisms.

The temporal behavior of the murine T cell receptor repertoire following Immunotherapy

Immunotherapy is now an essential tool for cancer treatment, and the unique features of an individual's T cell receptor repertoire are known to play a key role in its effectiveness. The repertoire, famously vast due to a cascade of cellular mechanisms, can be quantified using repertoire sequencing. In this study, we sampled the repertoire over several time points following treatment with anti-CTLA-4, in a syngeniec mouse model for colorectal cancer, generating a longitudinal dataset of T cell clones and their abundance. The dynamics of the repertoire can be observed in response to treatment over a period of four weeks, as clonal expansion of specific clones ascends and descends. The data made available here can be used to determine treatment and predict its effect, while also providing a unique look at the behavior of the immune system over time.

A Comprehensive View of the Cancer-Immunity Cycle (CIC) in HPV-Mediated Cervical Cancer and Prospects for Emerging Therapeutic Opportunities

Cervical cancer (CC) is the fourth most common cancer in women worldwide, with more than 500,000 new cases each year and a mortality rate of around 55%. Over 80% of these deaths occur in developing countries. The most important risk factor for CC is persistent infection by a sexually transmitted virus, the human papillomavirus (HPV). Conventional treatments to eradicate this type of cancer are accompanied by high rates of resistance and a large number of side effects. Hence, it is crucial to devise novel effective therapeutic strategies. In recent years, an increasing number of studies have aimed to develop immunotherapeutic methods for treating cancer. However, these strategies have not proven to be effective enough to combat CC. This means there is a need to investigate immune molecular targets. An adaptive immune response against cancer has been described in seven key stages or steps defined as the cancer-immunity cycle (CIC). The CIC begins with the release of antigens by tumor cells and ends with their destruction by cytotoxic T-cells. In this paper, we discuss several molecular alterations found in each stage of the CIC of CC. In addition, we analyze the evidence discovered, the molecular mechanisms and their relationship with variables such as histological subtype and HPV infection, as well as their potential impact for adopting novel immunotherapeutic approaches.

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

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

Machine learning identifies T cell receptor repertoire signatures associated with COVID-19 severity

T cell receptor (TCR) repertoires are critical for antiviral immunity. Determining the TCR repertoire composition, diversity, and dynamics and how they change during viral infection can inform the molecular specificity of host responses to viruses such as SARS-CoV-2. To determine signatures associated with COVID-19 disease severity, here we perform a large-scale analysis of over 4.7 billion sequences across 2130 TCR repertoires from COVID-19 patients and healthy donors. TCR repertoire analyses from these data identify and characterize convergent COVID-19-associated CDR3 gene usages, specificity groups, and sequence patterns. Here we show that T cell clonal expansion is associated with the upregulation of T cell effector function, TCR signaling, NF-kB signaling, and interferon-gamma signaling pathways. We also demonstrate that machine learning approaches accurately predict COVID-19 infection based on TCR sequence features, with certain high-power models reaching near-perfect AUROC scores. These analyses provide a systems immunology view of T cell adaptive immune responses to COVID-19.

The evolutionary and functional significance of germline immunoglobulin gene variation

The recombination between immunoglobulin (IG) gene segments determines an individual's naïve antibody repertoire and, consequently, (auto)antigen recognition. Emerging evidence suggests that mammalian IG germline variation impacts humoral immune responses associated with vaccination, infection, and autoimmunity - from the molecular level of epitope specificity, up to profound changes in the architecture of antibody repertoires. These links between IG germline variants and immunophenotype raise the question on the evolutionary causes and consequences of diversity within IG loci. We discuss why the extreme diversity in IG loci remains a mystery, why resolving this is important for the design of more effective vaccines and therapeutics, and how recent evidence from multiple lines of inquiry may help us do so.

simAIRR: simulation of adaptive immune repertoires with realistic receptor sequence sharing for benchmarking of immune state prediction methods

BACKGROUND

Machine learning (ML) has gained significant attention for classifying immune states in adaptive immune receptor repertoires (AIRRs) to support the advancement of immunodiagnostics and therapeutics. Simulated data are crucial for the rigorous benchmarking of AIRR-ML methods. Existing approaches to generating synthetic benchmarking datasets result in the generation of naive repertoires missing the key feature of many shared receptor sequences (selected for common antigens) found in antigen-experienced repertoires.

RESULTS

We demonstrate that a common approach to generating simulated AIRR benchmark datasets can introduce biases, which may be exploited for undesired shortcut learning by certain ML methods. To mitigate undesirable access to true signals in simulated AIRR datasets, we devised a simulation strategy (simAIRR) that constructs antigen-experienced-like repertoires with a realistic overlap of receptor sequences. simAIRR can be used for constructing AIRR-level benchmarks based on a range of assumptions (or experimental data sources) for what constitutes receptor-level immune signals. This includes the possibility of making or not making any prior assumptions regarding the similarity or commonality of immune state-associated sequences that will be used as true signals. We demonstrate the real-world realism of our proposed simulation approach by showing that basic ML strategies perform similarly on simAIRR-generated and real-world experimental AIRR datasets.

CONCLUSIONS

This study sheds light on the potential shortcut learning opportunities for ML methods that can arise with the state-of-the-art way of simulating AIRR datasets. simAIRR is available as a Python package: https://github.com/KanduriC/simAIRR.

Subtle Longitudinal Alterations in Env Sequence Potentiate Differences in Sensitivity to Broadly Neutralizing Antibodies following Acute HIV-1 Subtype C Infection

Broadly neutralizing antibodies (bNAbs) for HIV-1 prevention or cure strategies must inhibit transmitted/founder and reservoir viruses. Establishing sensitivity of circulating viruses to bNAbs and genetic patterns affecting neutralization variability may guide rational bNAbs selection for clinical development. We analyzed 326 single env genomes from nine individuals followed longitudinally following acute HIV-1 infection, with samples collected at ~1 week after the first detection of plasma viremia; 300 to 1,709 days postinfection but prior to initiating antiretroviral therapy (ART) (median = 724 days); and ~1 year post ART initiation. Sequences were assessed for phylogenetic relatedness, potential N- and O-linked glycosylation, and variable loop lengths (V1 to V5). A total of 43 env amplicons (median = 3 per patient per time point) were cloned into an expression vector and the TZM-bl assay was used to assess the neutralization profiles of 15 bNAbs targeting the CD4 binding site, V1/V2 region, V3 supersite, MPER, gp120/gp41 interface, and fusion peptide. At 1 μg/mL, the neutralization breadths were as follows: VRC07-LS and N6.LS (100%), VRC01 (86%), PGT151 (81%), 10-1074 and PGT121 (80%), and less than 70% for 10E8, 3BNC117, CAP256.VRC26, 4E10, PGDM1400, and N123-VRC34.01. Features associated with low sensitivity to V1/V2 and V3 bNAbs were higher potential glycosylation sites and/or relatively longer V1 and V4 domains, including known "signature" mutations. The study shows significant variability in the breadth and potency of bNAbs against circulating HIV-1 subtype C envelopes. VRC07-LS, N6.LS, VRC01, PGT151, 10-1074, and PGT121 display broad activity against subtype C variants, and major determinants of sensitivity to most bNAbs were within the V1/V4 domains. IMPORTANCE Broadly neutralizing antibodies (bNAbs) have potential clinical utility in HIV-1 prevention and cure strategies. However, bNAbs target diverse epitopes on the HIV-1 envelope and the virus may evolve to evade immune responses. It is therefore important to identify antibodies with broad activity in high prevalence settings, as well as the genetic patterns that may lead to neutralization escape. We investigated 15 bNAbs with diverse biophysical properties that target six epitopes of the HIV-1 Env glycoprotein for their ability to inhibit viruses that initiated infection, viruses circulating in plasma at chronic infection before antiretroviral treatment (ART), or viruses that were archived in the reservoir during ART in subtype C infected individuals in South Africa, a high burden country. We identify the antibodies most likely to be effective for clinical use in this setting and describe mutational patterns associated with neutralization escape from these antibodies.

A scalable model for simulating multi-round antibody evolution and benchmarking of clonal tree reconstruction methods

Affinity maturation (AM) of B cells through somatic hypermutations (SHMs) enables the immune system to evolve to recognize diverse pathogens. The accumulation of SHMs leads to the formation of clonal lineages of antibody-secreting b cells that have evolved from a common naïve B cell. Advances in high-throughput sequencing have enabled deep scans of B cell receptor repertoires, paving the way for reconstructing clonal trees. However, it is not clear if clonal trees, which capture microevolutionary time scales, can be reconstructed using traditional phylogenetic reconstruction methods with adequate accuracy. In fact, several clonal tree reconstruction methods have been developed to fix supposed shortcomings of phylogenetic methods. Nevertheless, no consensus has been reached regarding the relative accuracy of these methods, partially because evaluation is challenging. Benchmarking the performance of existing methods and developing better methods would both benefit from realistic models of clonal lineage evolution specifically designed for emulating B cell evolution. In this paper, we propose a model for modeling B cell clonal lineage evolution and use this model to benchmark several existing clonal tree reconstruction methods. Our model, designed to be extensible, has several features: by evolving the clonal tree and sequences simultaneously, it allows modeling selective pressure due to changes in affinity binding; it enables scalable simulations of large numbers of cells; it enables several rounds of infection by an evolving pathogen; and, it models building of memory. In addition, we also suggest a set of metrics for comparing clonal trees and measuring their properties. Our results show that while maximum likelihood phylogenetic reconstruction methods can fail to capture key features of clonal tree expansion if applied naively, a simple post-processing of their results, where short branches are contracted, leads to inferences that are better than alternative methods.

Exploration of shared features of B cell receptor and T cell receptor repertoires reveals distinct clonotype clusters

Although B cells and T cells are integral players of the adaptive immune system and act in co-dependent ways to orchestrate immune responses, existing methods to study the immune repertoire have largely focused on separate analyses of B cell receptor (BCR) and T cell receptor (TCR) repertoires. Based on our hypothesis that the shared history of immune exposures and the shared cellular machinery for recombination result in similarities between BCR and TCR repertoires in an individual, we examine any commonalities and interrelationships between BCR and TCR repertoires. We find that the BCR and TCR repertoires have covarying clonal architecture and diversity, and that the pattern of correlations appears to be altered in immune-mediated diseases. Furthermore, hierarchical clustering of public B and T cell clonotypes in both health and disease based on correlation of clonal proportion revealed distinct clusters of B and T cell clonotypes that exhibit increased sequence similarity, share motifs, and have distinct amino acid characteristics. Our findings point to common principles governing memory formation, recombination, and clonal expansion to antigens in B and T cells within an individual. A significant proportion of public BCR and TCR repertoire can be clustered into nonoverlapping and correlated clusters, suggesting a novel way of grouping B and T cell clonotypes.

Evolution and modulation of antigen-specific T cell responses in melanoma patients

Analyzing antigen-specific T cell responses at scale has been challenging. Here, we analyze three types of T cell receptor (TCR) repertoire data (antigen-specific TCRs, TCR-repertoire, and single-cell RNA + TCRαβ-sequencing data) from 515 patients with primary or metastatic melanoma and compare it to 783 healthy controls. Although melanoma-associated antigen (MAA) -specific TCRs are restricted to individuals, they share sequence similarities that allow us to build classifiers for predicting anti-MAA T cells. The frequency of anti-MAA T cells distinguishes melanoma patients from healthy and predicts metastatic recurrence from primary melanoma. Anti-MAA T cells have stem-like properties and frequent interactions with regulatory T cells and tumor cells via Galectin9-TIM3 and PVR-TIGIT -axes, respectively. In the responding patients, the number of expanded anti-MAA clones are higher after the anti-PD1(+anti-CTLA4) therapy and the exhaustion phenotype is rescued. Our systems immunology approach paves the way for understanding antigen-specific responses in human disorders.

Previous studies have characterized the diversity and dynamics of the T cell receptor (TCR) repertoire in patients with solid cancer. Here, by analyzing TCR repertoire data from multiple datasets, the authors report that melanoma-associated antigen-specific TCRs can be used to separate metastatic melanoma patients from healthy controls and to follow anti-tumor responses in patients treated with immunotherapy.

T cell Repertoire Profiling and the Mechanism by which HLA-B27 Causes Ankylosing Spondylitis

Purpose of Review

Ankylosing spondylitis (AS) is strongly associated with the HLA-B27 gene. The canonical function of HLA-B27 is to present antigenic peptides to CD8 lymphocytes, leading to adaptive immune responses. The ‘arthritogenic peptide’ theory as to the mechanism by which HLA-B27 induces ankylosing spondylitis proposes that HLA-B27 presents peptides derived from exogenous sources such as bacteria to CD8 lymphocytes, which subsequently cross-react with antigens at the site of inflammation of the disease, causing inflammation. This review describes findings of studies in AS involving profiling of T cell expansions and discusses future research opportunities based on these findings.

Recent Findings

Consistent with this theory, there is an expanding body of data showing that expansion of a restricted pool of CD8 lymphocytes is found in most AS patients yet only in a small proportion of healthy HLA-B27 carriers.

Summary

These exciting findings strongly support the theory that AS is driven by presentation of antigenic peptides to the adaptive immune system by HLA-B27. They point to new potential approaches to identify the exogenous and endogenous antigens involved and to potential therapies for the disease.

A guide to systems-level immunomics

The immune system is highly complex and distributed throughout an organism, with hundreds to thousands of cell states existing in parallel with diverse molecular pathways interacting in a highly dynamic and coordinated fashion. Although the characterization of individual genes and molecules is of the utmost importance for understanding immune-system function, high-throughput, high-resolution omics technologies combined with sophisticated computational modeling and machine-learning approaches are creating opportunities to complement standard immunological methods with new insights into immune-system dynamics. Like systems immunology itself, immunology researchers must take advantage of these technologies and form their own diverse networks, connecting with researchers from other disciplines. This Review is an introduction and 'how-to guide' for immunologists with no particular experience in the field of omics but with the intention to learn about and apply these systems-level approaches, and for immunologists who want to make the most of interdisciplinary networks.

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

Motivation

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

Results

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

Availability and implementation

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

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

"Common variable immunodeficiency: Challenges for diagnosis"

Common variable immunodeficiency is a heterogeneous condition characterized by B cell dysfunction with reduced serum immunoglobulin levels and a highly variable spectrum of clinical manifestations ranging from recurrent infections to autoimmune disease. The diagnosis of CVID is often challenging due to the diverse clinical presentation of patients and the existence of multiple diagnostic criteria without a universally adopted consensus. Laboratory evaluation to assist with diagnosis currently includes serum immunoglobulin testing, immunophenotyping, assessment of vaccine response, and genetic testing. Additional emerging techniques include investigation of the B cell repertoire and the use of machine learning algorithms. Advances in our understanding of common variable immunodeficiency will ultimately contribute to earlier diagnosis and novel interventions with the goal of improving prognosis for these patients.

Reference-based comparison of adaptive immune receptor repertoires

B and T cell receptor (immune) repertoires can represent an individual's immune history. While current repertoire analysis methods aim to discriminate between health and disease states, they are typically based on only a limited number of parameters. Here, we introduce immuneREF: a quantitative multidimensional measure of adaptive immune repertoire (and transcriptome) similarity that allows interpretation of immune repertoire variation by relying on both repertoire features and cross-referencing of simulated and experimental datasets. To quantify immune repertoire similarity landscapes across health and disease, we applied immuneREF to >2,400 datasets from individuals with varying immune states (healthy, [autoimmune] disease, and infection). We discovered, in contrast to the current paradigm, that blood-derived immune repertoires of healthy and diseased individuals are highly similar for certain immune states, suggesting that repertoire changes to immune perturbations are less pronounced than previously thought. In conclusion, immuneREF enables the population-wide study of adaptive immune response similarity across immune states.

Machine Learning Approaches to TCR Repertoire Analysis

Sparked by the development of genome sequencing technology, the quantity and quality of data handled in immunological research have been changing dramatically. Various data and database platforms are now driving the rapid progress of machine learning for immunological data analysis. Of various topics in immunology, T cell receptor repertoire analysis is one of the most important targets of machine learning for assessing the state and abnormalities of immune systems. In this paper, we review recent repertoire analysis methods based on machine learning and deep learning and discuss their prospects.