Modeling and predicting the overlap of B- and T-cell receptor repertoires in healthy and SARS-CoV-2 infected individuals

Application of Immune Repertoire Analysis in Differentiating Thyroid Cancer and Large Benign Thyroid Nodules

This study compares the peripheral T-cell receptor (TCR) and B-cell receptor (BCR) immune repertoires among early-stage papillary thyroid carcinoma (PTC) patients, patients with benign thyroid nodules larger than 4 cm, and healthy controls. Adaptive immune repertoire sequencing is used to analyze peripheral immune profile differences among these groups. Results indicates that early PTC and large benign nodules show significantly higher proportions of expanded clones than healthy controls, reflecting antigen-driven clonal expansion. By introducing the concept of "publicness," disease-specific high-publicness clonotypes is identified. These clonotypes exhibits distinct V-J rearrangement characteristics and strong immune heterogeneity. This study further reveals that this immune heterogeneity may be associated with patients' thyroid hormone levels and autoimmune antibody levels. These findings provides new insights into the immunopathological mechanisms of thyroid disorders.

Learning predictive signatures of HLA type from T-cell repertoires

T cells recognize a wide range of pathogens using surface receptors that interact directly with peptides presented on major histocompatibility complexes (MHC) encoded by the HLA loci in humans. Understanding the association between T cell receptors (TCR) and HLA alleles is an important step towards predicting TCR-antigen specificity from sequences. Here we analyze the TCR alpha and beta repertoires of large cohorts of HLA-typed donors to systematically infer such associations, by looking for overrepresentation of TCRs in individuals with a common allele.TCRs, associated with a specific HLA allele, exhibit sequence similarities that suggest prior antigen exposure. Immune repertoire sequencing has produced large numbers of datasets, however the HLA type of the corresponding donors is rarely available. Using our TCR-HLA associations, we trained a computational model to predict the HLA type of individuals from their TCR repertoire alone. We propose an iterative procedure to refine this model by using data from large cohorts of untyped individuals, by recursively typing them using the model itself. The resulting model shows good predictive performance, even for relatively rare HLA alleles.

The characteristics of TCR CDR3 repertoire in COVID-19 patients and SARS-CoV-2 vaccine recipients

The COVID-19 pandemic and large-scale administration of multiple SARS-CoV-2 vaccines have attracted global attention to the short-term and long-term effects on the human immune system. An analysis of the "traces" left by the body's T-cell immune response is needed, especially for the prevention and treatment of breakthrough infections and long COVID-19 and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant infections. T-cell receptor complementarity determining region 3 (TCR CDR3) repertoire serves as a target molecule for monitoring the effects, mechanisms, and memory of the T-cell response. Furthermore, it has been extensively applied in the elucidation of the infectious mechanism and vaccine refinement of hepatitis B virus (HBV), influenza virus, human immunodeficiency virus (HIV), and SARS-CoV. Laboratories worldwide have utilized high-throughput sequencing (HTS) and scTCR-seq to characterize, share, and apply the TCR CDR3 repertoire in COVID-19 patients and SARS-CoV-2 vaccine recipients. This article focuses on the comparative analysis of the diversity, clonality, V&J gene usage and pairing, CDR3 length, shared CDR3 sequences or motifs, and other characteristics of TCR CDR3 repertoire. These findings provide molecular targets for evaluating T-cell response effects and short-term and long-term impacts on the adaptive immune system following SARS-CoV-2 infection or vaccination and establish a comparative archive of T-cell response "traces."

Chemical-Specific T Cell Tests Aim to Bridge a Gap in Skin Sensitization Evaluation

T cell activation is the final key event (KE4) in the adverse outcome pathway (AOP) of skin sensitization. However, validated new approach methodologies (NAMs) for evaluating this step are missing. Accordingly, chemicals that activate an unusually high frequency of T cells, as does the most prevalent metal allergen nickel, are not yet identified in a regulatory context. T cell reactivity to chemical sensitizers might be especially relevant in real-life scenarios, where skin injury, co-exposure to irritants in chemical mixtures, or infections may trigger the heterologous innate immune stimulation necessary to induce adaptive T cell responses. Additionally, cross-reactivity, which underlies cross-allergies, can only be assessed by T cell tests. To date, several experimental T cell tests are available that use primary naïve and memory CD4+ and CD8+ T cells from human blood. These include priming and lymphocyte proliferation tests and, most recently, activation-induced marker (AIM) assays. All approaches are challenged by chemical-mediated toxicity, inefficient or unknown generation of T cell epitopes, and a low throughput. Here, we summarize solutions and strategies to confirm in vitro T cell signals. Broader application and standardization are necessary to possibly define chemical applicability domains and to strengthen the role of T cell tests in regulatory risk assessment.

Computational detection of antigen-specific B cell receptors following immunization

B cell receptors (BCRs) play a crucial role in recognizing and fighting foreign antigens. High-throughput sequencing enables in-depth sampling of the BCRs repertoire after immunization. However, only a minor fraction of BCRs actively participate in any given infection. To what extent can we accurately identify antigen-specific sequences directly from BCRs repertoires? We present a computational method grounded on sequence similarity, aimed at identifying statistically significant responsive BCRs. This method leverages well-known characteristics of affinity maturation and expected diversity. We validate its effectiveness using longitudinally sampled human immune repertoire data following influenza vaccination and SARS-CoV-2 infections. We show that different lineages converge to the same responding Complementarity Determining Region 3, demonstrating convergent selection within an individual. The outcomes of this method hold promise for application in vaccine development, personalized medicine, and antibody-derived therapeutics.

Combining mutation and recombination statistics to infer clonal families in antibody repertoires

B-cell repertoires are characterized by a diverse set of receptors of distinct specificities generated through two processes of somatic diversification: V(D)J recombination and somatic hypermutations. B-cell clonal families stem from the same V(D)J recombination event, but differ in their hypermutations. Clonal families identification is key to understanding B-cell repertoire function, evolution, and dynamics. We present HILARy (high-precision inference of lineages in antibody repertoires), an efficient, fast, and precise method to identify clonal families from single- or paired-chain repertoire sequencing datasets. HILARy combines probabilistic models that capture the receptor generation and selection statistics with adapted clustering methods to achieve consistently high inference accuracy. It automatically leverages the phylogenetic signal of shared mutations in difficult repertoire subsets. Exploiting the high sensitivity of the method, we find the statistics of evolutionary properties such as the site frequency spectrum and dN/dS ratio do not depend on the junction length. We also identify a broad range of selection pressures spanning two orders of magnitude.

Robust detection of infectious disease, autoimmunity, and cancer from the paratope networks of adaptive immune receptors

Liquid biopsies based on peripheral blood offer a minimally invasive alternative to solid tissue biopsies for the detection of diseases, primarily cancers. However, such tests currently consider only the serum component of blood, overlooking a potentially rich source of biomarkers: adaptive immune receptors (AIRs) expressed on circulating B and T cells. Machine learning-based classifiers trained on AIRs have been reported to accurately identify not only cancers but also autoimmune and infectious diseases as well. However, when using the conventional "clonotype cluster" representation of AIRs, individuals within a disease or healthy cohort exhibit vastly different features, limiting the generalizability of these classifiers. This study aimed to address the challenge of classifying specific diseases from circulating B or T cells by developing a novel representation of AIRs based on similarity networks constructed from their antigen-binding regions (paratopes). Features based on this novel representation, paratope cluster occupancies (PCOs), significantly improved disease classification performance for infectious disease, autoimmune disease, and cancer. Under identical methodological conditions, classifiers trained on PCOs achieved a mean AUC of 0.893 when applied to new individuals, outperforming clonotype cluster-based classifiers (AUC 0.714) and the best-performing published classifier (AUC 0.777). Surprisingly, for cancer patients, we observed that "healthy-biased" AIRs were predicted to target known cancer-associated antigens at dramatically higher rates than healthy AIRs as a whole (Z scores >75), suggesting an overlooked reservoir of cancer-targeting immune cells that could be identified by PCOs.

Deep learning predictions of TCR-epitope interactions reveal epitope-specific chains in dual alpha T cells

T cells have the ability to eliminate infected and cancer cells and play an essential role in cancer immunotherapy. T cell activation is elicited by the binding of the T cell receptor (TCR) to epitopes displayed on MHC molecules, and the TCR specificity is determined by the sequence of its α and β chains. Here, we collect and curate a dataset of 17,715 αβTCRs interacting with dozens of class I and class II epitopes. We use this curated data to develop MixTCRpred, an epitope-specific TCR-epitope interaction predictor. MixTCRpred accurately predicts TCRs recognizing several viral and cancer epitopes. MixTCRpred further provides a useful quality control tool for multiplexed single-cell TCR sequencing assays of epitope-specific T cells and pinpoints a substantial fraction of putative contaminants in public databases. Analysis of epitope-specific dual α T cells demonstrates that MixTCRpred can identify α chains mediating epitope recognition. Applying MixTCRpred to TCR repertoires from COVID-19 patients reveals enrichment of clonotypes predicted to bind an immunodominant SARS-CoV-2 epitope. Overall, MixTCRpred provides a robust tool to predict TCRs interacting with specific epitopes and interpret TCR-sequencing data from both bulk and epitope-specific T cells.

Systematic evaluation of B-cell clonal family inference approaches

The reconstruction of clonal families (CFs) in B-cell receptor (BCR) repertoire analysis is a crucial step to understand the adaptive immune system and how it responds to antigens. The BCR repertoire of an individual is formed throughout life and is diverse due to several factors such as gene recombination and somatic hypermutation. The use of Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) using next generation sequencing enabled the generation of full BCR repertoires that also include rare CFs. The reconstruction of CFs from AIRR-seq data is challenging and several approaches have been developed to solve this problem. Currently, most methods use the heavy chain (HC) only, as it is more variable than the light chain (LC). CF reconstruction options include the definition of appropriate sequence similarity measures, the use of shared mutations among sequences, and the possibility of reconstruction without preliminary clustering based on V- and J-gene annotation. In this study, we aimed to systematically evaluate different approaches for CF reconstruction and to determine their impact on various outcome measures such as the number of CFs derived, the size of the CFs, and the accuracy of the reconstruction. The methods were compared to each other and to a method that groups sequences based on identical junction sequences and another method that only determines subclones. We found that after accounting for data set variability, in particular sequencing depth and mutation load, the reconstruction approach has an impact on part of the outcome measures, including the number of CFs. Simulations indicate that unique junctions and subclones should not be used as substitutes for CF and that more complex methods do not outperform simpler methods. Also, we conclude that different approaches differ in their ability to correctly reconstruct CFs when not considering the LC and to identify shared CFs. The results showed the effect of different approaches on the reconstruction of CFs and highlighted the importance of choosing an appropriate method.

Deciphering the antigen specificities of antibodies by clustering their complementarity determining region sequences

IMPORTANCE

Determining antigen and epitope specificity is an essential step in the discovery of therapeutic antibodies as well as in the analysis adaptive immune responses to disease or vaccination. Despite extensive efforts, deciphering antigen specificity solely from BCR amino acid sequence remains a challenging task, requiring a combination of experimental and computational approaches. Here, we describe and experimentally validate a simple and straightforward approach for grouping antibodies that share antigen and epitope specificities based on their CDR sequence similarity. This approach allows us to identify the specificities of a large number of antibodies whose antigen targets are unknown, using a small fraction of antibodies with well-annotated binding specificities.

Mathematical Characterization of Private and Public Immune Receptor Sequences

Diverse T and B cell repertoires play an important role in mounting effective immune responses against a wide range of pathogens and malignant cells. The number of unique T and B cell clones is characterized by T and B cell receptors (TCRs and BCRs), respectively. Although receptor sequences are generated probabilistically by recombination processes, clinical studies found a high degree of sharing of TCRs and BCRs among different individuals. In this work, we use a general probabilistic model for T/B cell receptor clone abundances to define "publicness" or "privateness" and information-theoretic measures for comparing the frequency of sampled sequences observed across different individuals. We derive mathematical formulae to quantify the mean and the variances of clone richness and overlap. Our results can be used to evaluate the effect of different sampling protocols on abundances of clones within an individual as well as the commonality of clones across individuals. Using synthetic and empirical TCR amino acid sequence data, we perform simulations to study expected clonal commonalities across multiple individuals. Based on our formulae, we compare these simulated results with the analytically predicted mean and variances of the repertoire overlap. Complementing the results on simulated repertoires, we derive explicit expressions for the richness and its uncertainty for specific, single-parameter truncated power-law probability distributions. Finally, the information loss associated with grouping together certain receptor sequences, as is done in spectratyping, is also evaluated. Our approach can be, in principle, applied under more general and mechanistically realistic clone generation models.

The Longitudinal Analysis of Convergent Antibody VDJ Regions in SARS-CoV-2-Positive Patients Using RNA-Seq

Severe acute respiratory syndrome-related coronavirus-2 (SARS-CoV-2) is an ongoing pandemic that continues to evolve and reinfect individuals. To understand the convergent antibody responses that evolved over the course of the pandemic, we evaluated the immunoglobulin repertoire of individuals infected by different SARS-CoV-2 variants for similarity between patients. We utilized four public RNA-seq data sets collected between March 2020 and March 2022 from the Gene Expression Omnibus (GEO) in our longitudinal analysis. This covered individuals infected with Alpha and Omicron variants. In total, from 269 SARS-CoV-2-positive patients and 26 negative patients, 629,133 immunoglobulin heavy-chain variable region V(D)J sequences were reconstructed from sequencing data. We grouped samples based on the SARS-CoV-2 variant type and/or the time they were collected from patients. Our comparison of patients within each SARS-CoV-2-positive group found 1011 common V(D)Js (same V gene, J gene and CDR3 amino acid sequence) shared by more than one patient and no common V(D)Js in the noninfected group. Taking convergence into account, we clustered based on similar CDR3 sequence and identified 129 convergent clusters from the SARS-CoV-2-positive groups. Within the top 15 clusters, 4 contain known anti-SARS-CoV-2 immunoglobulin sequences with 1 cluster confirmed to cross-neutralize variants from Alpha to Omicron. In our analysis of longitudinal groups that include Alpha and Omicron variants, we find that 2.7% of the common CDR3s found within groups were also present in more than one group. Our analysis reveals common and convergent antibodies, which include anti-SARS-CoV-2 antibodies, in patient groups over various stages of the pandemic.