Learning the statistics and landscape of somatic mutation-induced insertions and deletions in antibodies

Unveiling inverted D genes and D-D fusions in human antibody repertoires unlocks novel antibody diversity

Antibodies, essential components of adaptive immunity, derive their remarkable diversity primarily from V(D)J gene rearrangements, particularly within the heavy chain complementarity-determining region 3 (CDR-H3) where D genes play a major role. Traditionally, D genes were thought to recombine only in the forward direction, despite having identical recombination signal sequences (12 base pair spacers) at both ends. This observation led us to question whether these symmetrical sequences might enable bidirectional recombination. We identified 25 unique inverted D genes (InvDs) in both naive and memory B cells from antibody repertoires of 13 healthy donors. These InvDs utilize all three reading frames during translation, producing distinct amino acid profiles enriched in histidine, proline, and lysine in CDR-H3s of antibodies with potential functional diversity. Notably, our analysis revealed a broader range of D-D fusions, including D-D, D-InvD, InvD-D, and InvD-InvD configurations, opening new perspectives for antibody engineering and therapeutic development.

Rapid affinity optimization of an anti-TREM2 clinical lead antibody by cross-lineage immune repertoire mining

We describe a process for rapid antibody affinity optimization by repertoire mining to identify clones across B cell clonal lineages based on convergent immune responses where antigen-specific clones with the same heavy (VH) and light chain germline segment pairs, or parallel lineages, bind a single epitope on the antigen. We use this convergence framework to mine unique and distinct VH lineages from rat anti-triggering receptor on myeloid cells 2 (TREM2) antibody repertoire datasets with high diversity in the third complementarity-determining loop region (CDR H3) to further affinity-optimize a high-affinity agonistic anti-TREM2 antibody while retaining critical functional properties. Structural analyses confirm a nearly identical binding mode of anti-TREM2 variants with subtle but significant structural differences in the binding interface. Parallel lineage repertoire mining is uniquely tailored to rationally explore the large CDR H3 sequence space in antibody repertoires and can be easily and generally applied to antibodies discovered in vivo.

Isotype-aware inference of B cell clonal lineage trees from single-cell sequencing data

Single-cell RNA sequencing (scRNA-seq) enables comprehensive characterization of the micro-evolutionary processes of B cells during an adaptive immune response, capturing features of somatic hypermutation (SHM) and class switch recombination (CSR). Existing phylogenetic approaches for reconstructing B cell evolution have primarily focused on the SHM process alone. Here, we present tree inference of B cell clonal lineages (TRIBAL), an algorithm designed to optimally reconstruct the evolutionary history of B cell clonal lineages undergoing both SHM and CSR from scRNA-seq data. Through simulations, we demonstrate that TRIBAL produces more comprehensive and accurate B cell lineage trees compared to existing methods. Using real-world datasets, TRIBAL successfully recapitulates expected biological trends in a model affinity maturation system while reconstructing evolutionary histories with more parsimonious class switching than state-of-the-art methods. Thus, TRIBAL significantly improves B cell lineage tracing, useful for modeling vaccine responses, disease progression, and the identification of therapeutic antibodies.

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.

B cell phylogenetics in the single cell era

The widespread availability of single-cell RNA sequencing (scRNA-seq) has led to the development of new methods for understanding immune responses. Single-cell transcriptome data can now be paired with B cell receptor (BCR) sequences. However, RNA from BCRs cannot be analyzed like most other genes because BCRs are genetically diverse within individuals. In humans, BCRs are shaped through recombination followed by mutation and selection for antigen binding. As these processes co-occur with cell division, B cells can be studied using phylogenetic trees representing the mutations within a clone. B cell trees can link experimental timepoints, tissues, or cellular subtypes. Here, we review the current state and potential of how B cell phylogenetics can be combined with single-cell data to understand immune responses.

Large-scale data mining of four billion human antibody variable regions reveals convergence between therapeutic and natural antibodies that constrains search space for biologics drug discovery

The naïve human antibody repertoire has theoretical access to an estimated > 1015 antibodies. Identifying subsets of this prohibitively large space where therapeutically relevant antibodies may be found is useful for development of these agents. It was previously demonstrated that, despite the immense sequence space, different individuals can produce the same antibodies. It was also shown that therapeutic antibodies, which typically follow seemingly unnatural development processes, can arise independently naturally. To check for biases in how the sequence space is explored, we data mined public repositories to identify 220 bioprojects with a combined seven billion reads. Of these, we created a subset of human bioprojects that we make available as the AbNGS database (https://naturalantibody.com/ngs/). AbNGS contains 135 bioprojects with four billion productive human heavy variable region sequences and 385 million unique complementarity-determining region (CDR)-H3s. We find that 270,000 (0.07% of 385 million) unique CDR-H3s are highly public in that they occur in at least five of 135 bioprojects. Of 700 unique therapeutic CDR-H3, a total of 6% has direct matches in the small set of 270,000. This observation extends to a match between CDR-H3 and V-gene call as well. Thus, the subspace of shared ('public') CDR-H3s shows utility for serving as a starting point for therapeutic antibody design.

TRIBAL: Tree Inference of B cell Clonal Lineages

B cells are a critical component of the adaptive immune system, responsible for producing antibodies that help protect the body from infections and foreign substances. Single cell RNA-sequencing (scRNA-seq) has allowed for both profiling of B cell receptor (BCR) sequences and gene expression. However, understanding the adaptive and evolutionary mechanisms of B cells in response to specific stimuli remains a significant challenge in the field of immunology. We introduce a new method, TRIBAL, which aims to infer the evolutionary history of clonally related B cells from scRNA-seq data. The key insight of TRIBAL is that inclusion of isotype data into the B cell lineage inference problem is valuable for reducing phylogenetic uncertainty that arises when only considering the receptor sequences. Consequently, the TRIBAL inferred B cell lineage trees jointly capture the somatic mutations introduced to the B cell receptor during affinity maturation and isotype transitions during class switch recombination. In addition, TRIBAL infers isotype transition probabilities that are valuable for gaining insight into the dynamics of class switching. Via in silico experiments, we demonstrate that TRIBAL infers isotype transition probabilities with the ability to distinguish between direct versus sequential switching in a B cell population. This results in more accurate B cell lineage trees and corresponding ancestral sequence and class switch reconstruction compared to competing methods. Using real-world scRNA-seq datasets, we show that TRIBAL recapitulates expected biological trends in a model affinity maturation system. Furthermore, the B cell lineage trees inferred by TRIBAL were equally plausible for the BCR sequences as those inferred by competing methods but yielded lower entropic partitions for the isotypes of the sequenced B cell. Thus, our method holds the potential to further advance our understanding of vaccine responses, disease progression, and the identification of therapeutic antibodies.

Fast clonal family inference from large-scale B cell repertoire sequencing data

Advances in high-throughput sequencing technologies have facilitated the large-scale characterization of B cell receptor (BCR) repertoires. However, the vast amount and high diversity of the BCR sequences pose challenges for efficient and biologically meaningful analysis. Here, we introduce fastBCR, an efficient computational approach for inferring B cell clonal families from massive BCR heavy chain sequences. We demonstrate that fastBCR substantially reduces the running time while ensuring high accuracy on simulated datasets with diverse numbers of B cell lineages and varying mutation rates. We apply fastBCR to real BCR sequencing data from peripheral blood samples of COVID-19 patients, showing that the inferred clonal families display disease-associated features, as well as corresponding antigen-binding specificity and affinity. Overall, our results demonstrate the advantages of fastBCR for analyzing BCR repertoire data, which will facilitate the identification of disease-associated antibodies and improve our understanding of the B cell immune response.

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

Adaptive immunity's success relies on the extraordinary diversity of protein receptors on B and T cell membranes. Despite this diversity, the existence of public receptors shared by many individuals gives hope for developing population-wide vaccines and therapeutics. Using probabilistic modeling, we show many of these public receptors are shared by chance in healthy individuals. This predictable overlap is driven not only by biases in the random generation process of receptors, as previously reported, but also by their common functional selection. However, the model underestimates sharing between repertoires of individuals infected with SARS-CoV-2, suggesting strong specific antigen-driven convergent selection. We exploit this discrepancy to identify COVID-associated receptors, which we validate against datasets of receptors with known viral specificity. We study their properties in terms of sequence features and network organization, and use them to design an accurate diagnostic tool for predicting SARS-CoV-2 status from repertoire data.

Lineage Reconstruction of In Vitro Identified Antigen-Specific Autoreactive B Cells from Adaptive Immune Receptor Repertoires

The emergence, survival, growth and maintenance of autoreactive (AR) B-cell clones, the hallmark of humoral autoimmunity, leave their footprints in B-cell receptor repertoires. Collecting IgH sequences related to polyreactive (PR) ones from adaptive immune receptor repertoire (AIRR) datasets make the reconstruction and analysis of PR/AR B-cell lineages possible. We developed a computational approach, named ImmChainTracer, to extract members and to visualize clonal relationships of such B-cell lineages. Our approach was successfully applied on the IgH repertoires of patients suffering from monogenic hypomorphic RAG1 and 2 deficiency (pRD) or polygenic systemic lupus erythematosus (SLE) autoimmune diseases to identify relatives of AR IgH sequences and to track their fate in AIRRs. Signs of clonal expansion, affinity maturation and class-switching events in PR/AR and non-PR/AR B-cell lineages were revealed. An extension of our method towards B-cell expansion caused by any trigger (e.g., infection, vaccination or antibody development) may provide deeper insight into antigen specific B-lymphogenesis.