Systems Analysis Reveals High Genetic and Antigen-Driven Predetermination of Antibody Repertoires throughout B Cell Development

Computational Analysis of T-Cell Receptor Repertoire Workflow: From T-Cell Isolation to Bioinformatics Analysis

The T-cell receptor (TCR) is the key molecule involved in the adaptive immune response. It is generated by the V(D)J recombination, responsible of the enormous diversity of the TCR repertoire, a crucial feature determining the individual capability to response to antigens and to build immunological memory. A pivotal role in the recognition of antigen is played by the hypervariable complementarity-determining region 3 (CDR3) of the V-beta chain of TCR. Investigating the CDR3 supports the understanding of the adaptive immune system dynamics in physiological processes, such as immune aging, and in disease, especially autoimmune disorders in which T cells are main actors. High-throughput sequencing (HTS) paved the way for a great progress in the investigation of TCR repertoire, enhancing the read depth in the process of library generation of sequencing and the number of samples that can be analyzed simultaneously. Therefore, the leverage of big datasets stressed the need to develop computational approach, by bioinformatics, to unravel the characteristics of the TCR repertoire.

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.

The physiological landscape and specificity of antibody repertoires are consolidated by multiple immunizations

Diverse antibody repertoires spanning multiple lymphoid organs (i.e., bone marrow, spleen, lymph nodes) form the foundation of protective humoral immunity. Changes in their composition across lymphoid organs are a consequence of B-cell selection and migration events leading to a highly dynamic and unique physiological landscape of antibody repertoires upon antigenic challenge (e.g., vaccination). However, to what extent B cells encoding identical or similar antibody sequences (clones) are distributed across multiple lymphoid organs and how this is shaped by the strength of a humoral response remains largely unexplored. Here, we performed an in-depth systems analysis of antibody repertoires across multiple distinct lymphoid organs of immunized mice and discovered that organ-specific antibody repertoire features (i.e., germline V-gene usage and clonal expansion profiles) equilibrated upon a strong humoral response (multiple immunizations and high serum titers). This resulted in a surprisingly high degree of repertoire consolidation, characterized by highly connected and overlapping B-cell clones across multiple lymphoid organs. Finally, we revealed distinct physiological axes indicating clonal migrations and showed that antibody repertoire consolidation directly correlated with antigen specificity. Our study uncovered how a strong humoral response resulted in a more uniform but redundant physiological landscape of antibody repertoires, indicating that increases in antibody serum titers were a result of synergistic contributions from antigen-specific B-cell clones distributed across multiple lymphoid organs. Our findings provide valuable insights for the assessment and design of vaccine strategies.

Bulk T cell repertoire sequencing (TCR-Seq) is a powerful technology for understanding inflammation-mediated diseases

Multiple alterations in the T cell repertoire were identified in many chronic inflammatory diseases such as inflammatory bowel disease, multiple sclerosis, and rheumatoid arthritis, suggesting that T cells might, directly or indirectly, be implicated in these pathologies. This has sparked a deep interest in characterizing disease-associated T cell clonotypes as well as in identifying and quantifying their contribution to the pathophysiology of different autoimmune and inflammation-mediated diseases. Bulk T cell repertoire sequencing (TCR-Seq) has emerged as a powerful method to profile the T cell repertoire of a sample in a high throughput fashion. Given the increasing utilization of TCR-Seq, we aimed here to provide a comprehensive, up-to-date review of the method, its extensions, and its ability to investigate chronic and autoimmune diseases. Specifically, we started by introducing the immunological basis of TCR repertoire generation and features, followed by discussing different experimental approach to perform TCR-Seq, then we describe different methods and frameworks for analyzing the generated datasets. Subsequently, different experimental techniques for investigating the antigenicity of T cell clonotypes are described. Lastly, we discuss recent studies that utilized TCR-Seq to understand different inflammation-mediated diseases, discuss fallbacks of the technology and potential future directions to overcome current limitations.

Learning antibody sequence constraints from allelic inclusion

Antibodies and B-cell receptors (BCRs) are produced by B cells, and are built of a heavy chain and a light chain. Although each B cell could express two different heavy chains and four different light chains, usually only a unique pair of heavy chain and light chain is expressed-a phenomenon known as allelic exclusion. However, a small fraction of naive-B cells violate allelic exclusion by expressing two productive light chains, one of which has impaired function; this has been called allelic inclusion. We demonstrate that these B cells can be used to learn constraints on antibody sequence. Using large-scale single-cell sequencing data from humans, we find examples of light chain allelic inclusion in thousands of naive-B cells, which is an order of magnitude larger than existing datasets. We train machine learning models to identify the abnormal sequences in these cells. The resulting models correlate with antibody properties that they were not trained on, including polyreactivity, surface expression, and mutation usage in affinity maturation. These correlations are larger than what is achieved by existing antibody modeling approaches, indicating that allelic inclusion data contains useful new information. We also investigate the impact of similar selection forces on the heavy chain in mouse, and observe that pairing with the surrogate light chain significantly restricts heavy chain diversity.

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.

Linguistics-based formalization of the antibody language as a basis for antibody language models

Apparent parallels between natural language and antibody sequences have led to a surge in deep language models applied to antibody sequences for predicting cognate antigen recognition. However, a linguistic formal definition of antibody language does not exist, and insight into how antibody language models capture antibody-specific binding features remains largely uninterpretable. Here we describe how a linguistic formalization of the antibody language, by characterizing its tokens and grammar, could address current challenges in antibody language model rule mining.

Neutralizing antibody responses against contemporary and future influenza A(H3N2) viruses in paradoxical clades elicited by repeated and single vaccinations

As one of the most effective measures to prevent seasonal influenza viruses, annual influenza vaccination is globally recommended. Nevertheless, evidence regarding the impact of repeated vaccination to contemporary and future influenza has been inconclusive. A total of 100 subjects singly or repeatedly immunized with influenza vaccines including 3C.2a1 or 3C.3a1 A(H3N2) during 2018-2019 and 2019-2020 influenza season were recruited. We investigated neutralization antibody by microneutralization assay using four antigenically distinct A(H3N2) viruses circulating from 2018 to 2023, and tracked the dynamics of B cell receptor (BCR) repertoire for consecutive vaccinations. We found that vaccination elicited cross-reactive antibody responses against future emerging strains. Broader neutralizing antibodies to A(H3N2) viruses and more diverse BCR repertoires were observed in the repeated vaccination. Meanwhile, a higher frequency of BCR sequences shared among the repeated-vaccinated individuals with consistently boosting antibody response was found than those with a reduced antibody response. Our findings suggest that repeated seasonal vaccination could broaden the breadth of antibody responses, which may improve vaccine protection against future emerging viruses.

Re-visiting humoral constitutive antibacterial heterogeneity in bloodstream infections

Although cellular immunity has garnered much attention in the era of single-cell technologies, humoral innate immunity has receded in priority due to its presumed limited roles. Hence, despite the long-recognised bactericidal activity of serum-a functional characteristic of constitutive humoral immunity-much remains unclear regarding mechanisms underlying its inter-individual heterogeneity and clinical implications in bloodstream infections. Recent work suggests that the immediate antimicrobial effect of humoral innate immunity contributes to suppression of the excessive inflammatory responses to infection by reducing the amount of pathogen-associated molecular patterns. In this Personal View, we propose the need to re-explore factors underlying the inter-individual heterogeneity in serum antibacterial competence as a new approach to better understand humoral innate immunity and revisit the clinical use of measuring serum antibacterial activity in the management of bacterial bloodstream infections. Given the current emphasis on subtyping sepsis, a serum bactericidal assay might prove useful in defining a distinct sepsis endotype, to enable more personalised management.

Leveraging Artificial Intelligence to Expedite Antibody Design and Enhance Antibody-Antigen Interactions

This perspective sheds light on the transformative impact of recent computational advancements in the field of protein therapeutics, with a particular focus on the design and development of antibodies. Cutting-edge computational methods have revolutionized our understanding of protein-protein interactions (PPIs), enhancing the efficacy of protein therapeutics in preclinical and clinical settings. Central to these advancements is the application of machine learning and deep learning, which offers unprecedented insights into the intricate mechanisms of PPIs and facilitates precise control over protein functions. Despite these advancements, the complex structural nuances of antibodies pose ongoing challenges in their design and optimization. Our review provides a comprehensive exploration of the latest deep learning approaches, including language models and diffusion techniques, and their role in surmounting these challenges. We also present a critical analysis of these methods, offering insights to drive further progress in this rapidly evolving field. The paper includes practical recommendations for the application of these computational techniques, supplemented with independent benchmark studies. These studies focus on key performance metrics such as accuracy and the ease of program execution, providing a valuable resource for researchers engaged in antibody design and development. Through this detailed perspective, we aim to contribute to the advancement of antibody design, equipping researchers with the tools and knowledge to navigate the complexities of this field.

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.

Revolutionizing Synthetic Antibody Design: Harnessing Artificial Intelligence and Deep Sequencing Big Data for Unprecedented Advances

Synthetic antibodies (Abs) represent a category of engineered proteins meticulously crafted to replicate the functions of their natural counterparts. Such Abs are generated in vitro, enabling advanced molecular alterations associated with antigen recognition, paratope site engineering, and biochemical refinements. In a parallel realm, deep sequencing has brought about a paradigm shift in molecular biology. It facilitates the prompt and cost-effective high-throughput sequencing of DNA and RNA molecules, enabling the comprehensive big data analysis of Ab transcriptomes, including specific regions of interest. Significantly, the integration of artificial intelligence (AI), based on machine- and deep- learning approaches, has fundamentally transformed our capacity to discern patterns hidden within deep sequencing big data, including distinctive Ab features and protein folding free energy landscapes. Ultimately, current AI advances can generate approximations of the most stable Ab structural configurations, enabling the prediction of de novo synthetic Abs. As a result, this manuscript comprehensively examines the latest and relevant literature concerning the intersection of deep sequencing big data and AI methodologies for the design and development of synthetic Abs. Together, these advancements have accelerated the exploration of antibody repertoires, contributing to the refinement of synthetic Ab engineering and optimizations, and facilitating advancements in the lead identification process.

The multilevel extensive diversity across the cynomolgus macaque captured by ultra-deep adaptive immune receptor repertoire sequencing

The extent to which AIRRs differ among and within individuals remains elusive. Via ultra-deep repertoire sequencing of 22 and 25 tissues in three cynomolgus macaques, respectively, we identified 84 and 114 novel IGHV and TRBV alleles, confirming 72 (85.71%) and 100 (87.72%) of them. The heterogeneous V gene usage patterns were influenced, in turn, by genetics, isotype (for BCRs only), tissue group, and tissue. A higher proportion of intragroup shared clones in the intestinal tissues than those in other tissues suggests a close intra-intestinal adaptive immunity network. Significantly higher mutation burdens in the public clones and the inter-tissue shared IgM and IgD clones indicate that they might target the shared antigens. This study reveals the extensive heterogeneity of the AIRRs at various levels and has broad fundamental and clinical implications. The data generated here will serve as an invaluable resource for future studies on adaptive immunity in health and diseases.

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.

Characterizing adjuvants' effects at murine immunoglobulin repertoire level

Generating large-scale, high-fidelity sequencing data is challenging and, furthermore, not much has been done to characterize adjuvants' effects at the repertoire level. Thus, we introduced an IgSeq pipeline that standardized library prep protocols and data analysis functions for accurate repertoire profiling. We then studied systemically effects of CpG and Alum on the Ig heavy chain repertoire using the ovalbumin (OVA) murine model. Ig repertoires of different tissues (spleen and bone marrow) and isotypes (IgG and IgM) were examined and compared in IGHV mutation, gene usage, CDR3 length, clonal diversity, and clonal selection. We found Ig repertoires of different compartments exhibited distinguishable profiles at the non-immunized steady state, and distinctions became more pronounced upon adjuvanted immunizations. Notably, Alum and CpG effects exhibited different tissue- and isotype-preferences. The former led to increased diversity of abundant clones in bone marrow, and the latter promoted the selection of IgG clones in both tissues.

Lupus-associated innate receptors drive extrafollicular evolution of autoreactive B cells

In systemic lupus erythematosus, recent findings highlight the extrafollicular (EF) pathway as prominent origin of autoantibody-secreting cells (ASCs). CD21loCD11c+ B cells, associated with aging, infection, and autoimmunity, are contributors to autoreactive EF ASCs but have an obscure developmental trajectory. To study EF kinetics of autoreactive B cell in tissue, we adoptively transferred WT and gene knockout B cell populations into the 564Igi mice - an autoreactive host enriched with autoantigens and T cell help. Time-stamped analyses revealed TLR7 dependence in early escape of peripheral B cell tolerance and establishment of a pre-ASC division program. We propose CD21lo cells as precursors to EF ASCs due to their elevated TLR7 sensitivity and proliferative nature. Blocking receptor function reversed CD21 loss and reduced effector cell generation, portraying CD21 as a differentiation initiator and a possible target for autoreactive B cell suppression. Repertoire analysis further delineated proto-autoreactive B cell selection and receptor evolution toward self-reactivity. This work elucidates receptor and clonal dynamics in EF development of autoreactive B cells, and establishes modular, native systems to probe mechanisms of autoreactivity.

Assessing the feasibility of statistical inference using synthetic antibody-antigen datasets

Simulation frameworks are useful to stress-test predictive models when data is scarce, or to assert model sensitivity to specific data distributions. Such frameworks often need to recapitulate several layers of data complexity, including emergent properties that arise implicitly from the interaction between simulation components. Antibody-antigen binding is a complex mechanism by which an antibody sequence wraps itself around an antigen with high affinity. In this study, we use a synthetic simulation framework for antibody-antigen folding and binding on a 3D lattice that include full details on the spatial conformation of both molecules. We investigate how emergent properties arise in this framework, in particular the physical proximity of amino acids, their presence on the binding interface, or the binding status of a sequence, and relate that to the individual and pairwise contributions of amino acids in statistical models for binding prediction. We show that weights learnt from a simple logistic regression model align with some but not all features of amino acids involved in the binding, and that predictive sequence binding patterns can be enriched. In particular, main effects correlated with the capacity of a sequence to bind any antigen, while statistical interactions were related to sequence specificity.

Characterizing the dynamics of BCR repertoire from repeated influenza vaccination

Yearly epidemics of seasonal influenza cause an enormous disease burden around the globe. An understanding of the rules behind the immune response with repeated vaccination still presents a significant challenge, which would be helpful for optimizing the vaccination strategy. In this study, 34 healthy volunteers with 16 vaccinated were recruited, and the dynamics of the BCR repertoire for consecutive vaccinations in two seasons were tracked. In terms of diversity, length, network, V and J gene segments usage, somatic hypermutation (SHM) rate and isotype, it was found that the overall changes were stronger in the acute phase of the first vaccination than the second vaccination. However, the V gene segments of IGHV4-39, IGHV3-9, IGHV3-7 and IGHV1-69 were amplified in the acute phase of the first vaccination, with IGHV3-7 dominant. On the other hand, for the second vaccination, the changes were dominated by IGHV1-69, with potential for coding broad neutralizing antibody. Additional analysis indicates that the application of V gene segment for IGHV3-7 in the acute phase of the first vaccination was due to the elevated usage of isotypes IgM and IgG3. While for IGHV1-69 in the second vaccination, it was contributed by isotypes IgG1 and IgG2. Finally, 41 public BCR clusters were identified in the vaccine group, with both IGHV3-7 and IGHV1-69 were involved and representative complementarity determining region 3 (CDR3) motifs were characterized. This study provides insights into the immune response dynamics following repeated influenza vaccination in humans and can inform universal vaccine design and vaccine strategies in the future.

Application and challenges of TCR and BCR sequencing to investigate T- and B-cell clonality in elastase-induced experimental murine abdominal aortic aneurysm

Background

An abdominal aortic aneurysm (AAA) is a life-threatening cardiovascular disease. Although its pathogenesis is still poorly understood, recent evidence suggests that AAA displays autoimmune disease characteristics. Particularly, T cells responding to AAA-related antigens in the aortic wall may contribute to an initial immune response. Single-cell RNA (scRNA) T cell receptor (TCR) and B cell receptor (BCR) sequencing is a powerful tool for investigating clonality. However, difficulties such as limited numbers of isolated cells must be considered during implementation and data analysis, making biological interpretation challenging. Here, we perform a representative single-cell immune repertoire analysis in experimental murine AAA and show a reliable bioinformatic processing pipeline highlighting opportunities and limitations of this approach.

Methods

We performed scRNA TCR and BCR sequencing of isolated lymphocytes from the infrarenal aorta of male C57BL/6J mice 3, 7, 14, and 28 days after AAA induction via elastase perfusion of the aorta. Sham-operated mice at days 3 and 28 and non-operated mice served as controls.

Results

Comparison of complementarity-determining region (CDR3) length distribution of 179 B cells and 796 T cells revealed neither differences between AAA and control nor between the disease stages. We found no clonal expansion of B cells in AAA. For T cells, we identified several clones in 11 of 16 AAA samples and one of eight control samples. Immune receptor repertoire comparison indicated that only a few clones were shared between the individual AAA samples. The most frequently used V-genes in the TCR beta chain in AAA were TRBV3, TRBV19, and the splicing variant TRBV12-2 + TRBV13-2.

Conclusion

We found no clonal expansion of B cells but evidence for clonal expansion of T cells in elastase-induced AAA in mice. Our findings imply that a more precise characterization of TCR and BCR distribution requires a more extensive number of lymphocytes to prevent undersampling and potentially detect rare clones. Thus, further experiments are necessary to confirm our findings. In summary, this paper examines TCR and BCR sequencing results, identifies limitations and pitfalls, and offers guidance for future studies.

Germline-encoded specificities and the predictability of the B cell response

Antibodies result from the competition of B cell lineages evolving under selection for improved antigen recognition, a process known as affinity maturation. High-affinity antibodies to pathogens such as HIV, influenza, and SARS-CoV-2 are frequently reported to arise from B cells whose receptors, the precursors to antibodies, are encoded by particular immunoglobulin alleles. This raises the possibility that the presence of particular germline alleles in the B cell repertoire is a major determinant of the quality of the antibody response. Alternatively, initial differences in germline alleles' propensities to form high-affinity receptors might be overcome by chance events during affinity maturation. We first investigate these scenarios in simulations: when germline-encoded fitness differences are large relative to the rate and effect size variation of somatic mutations, the same germline alleles persistently dominate the response of different individuals. In contrast, if germline-encoded advantages can be easily overcome by subsequent mutations, allele usage becomes increasingly divergent over time, a pattern we then observe in mice experimentally infected with influenza virus. We investigated whether affinity maturation might nonetheless strongly select for particular amino acid motifs across diverse genetic backgrounds, but we found no evidence of convergence to similar CDR3 sequences or amino acid substitutions. These results suggest that although germline-encoded specificities can lead to similar immune responses between individuals, diverse evolutionary routes to high affinity limit the genetic predictability of responses to infection and vaccination.

Characterisation of the immune repertoire of a humanised transgenic mouse through immunophenotyping and high-throughput sequencing

Immunoglobulin loci-transgenic animals are widely used in antibody discovery and increasingly in vaccine response modelling. In this study, we phenotypically characterised B-cell populations from the Intelliselect Transgenic mouse (Kymouse) demonstrating full B-cell development competence. Comparison of the naïve B-cell receptor (BCR) repertoires of Kymice BCRs, naïve human, and murine BCR repertoires revealed key differences in germline gene usage and junctional diversification. These differences result in Kymice having CDRH3 length and diversity intermediate between mice and humans. To compare the structural space explored by CDRH3s in each species' repertoire, we used computational structure prediction to show that Kymouse naïve BCR repertoires are more human-like than mouse-like in their predicted distribution of CDRH3 shape. Our combined sequence and structural analysis indicates that the naïve Kymouse BCR repertoire is diverse with key similarities to human repertoires, while immunophenotyping confirms that selected naïve B cells are able to go through complete development.

Longitudinal single-cell analysis of SARS-CoV-2-reactive B cells uncovers persistence of early-formed, antigen-specific clones

Understanding persistence and evolution of B cell clones after COVID-19 infection and vaccination is crucial for predicting responses against emerging viral variants and optimizing vaccines. Here, we collected longitudinal samples from patients with severe COVID-19 every third to seventh day during hospitalization and every third month after recovery. We profiled their antigen-specific immune cell dynamics by combining single-cell RNA-Seq, Cellular Indexing of Transcriptomes and Epitopes by Sequencing (CITE-Seq), and B cell receptor-Seq (BCR-Seq) with oligo-tagged antigen baits. While the proportion of Spike receptor binding domain-specific memory B cells (MBC) increased from 3 months after infection, the other Spike- and Nucleocapsid-specific B cells remained constant. All patients showed ongoing class switching and sustained affinity maturation of antigen-specific cells, and affinity maturation was not significantly increased early after vaccine. B cell analysis revealed a polyclonal response with limited clonal expansion; nevertheless, some clones detected during hospitalization, as plasmablasts, persisted for up to 1 year, as MBC. Monoclonal antibodies derived from persistent B cell families increased their binding and neutralization breadth and started recognizing viral variants by 3 months after infection. Overall, our findings provide important insights into the clonal evolution and dynamics of antigen-specific B cell responses in longitudinally sampled patients infected with 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.

Adaptive Immunity in Reptiles: Conventional Components but Unconventional Strategies

Recent studies have established that the innate immune system of reptiles is broad and robust, but the question remains: What role does the reptilian adaptive immune system play? Conventionally, adaptive immunity is described as involving T and B lymphocytes that display variable receptors, is highly specific, improves over the course of the response, and produces a memory response. While reptiles do have B and T lymphocytes that utilize variable receptors, their adaptive response is relatively non-specific, generates a prolonged antibody response, and does not produce a typical memory response. This alternative adaptive strategy may allow reptiles to produce a broad adaptive response that complements a strong innate system. Further studies into reptile adaptive immunity cannot only clarify outstanding questions on the reptilian immune system but can shed light on a number of important immunological concepts, including the evolution of the immune system and adaptive immune responses that take place outside of germinal centers.

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.

Unconstrained generation of synthetic antibody-antigen structures to guide machine learning methodology for antibody specificity prediction

Machine learning (ML) is a key technology for accurate prediction of antibody-antigen binding. Two orthogonal problems hinder the application of ML to antibody-specificity prediction and the benchmarking thereof: the lack of a unified ML formalization of immunological antibody-specificity prediction problems and the unavailability of large-scale synthetic datasets to benchmark real-world relevant ML methods and dataset design. Here we developed the Absolut! software suite that enables parameter-based unconstrained generation of synthetic lattice-based three-dimensional antibody-antigen-binding structures with ground-truth access to conformational paratope, epitope and affinity. We formalized common immunological antibody-specificity prediction problems as ML tasks and confirmed that for both sequence- and structure-based tasks, accuracy-based rankings of ML methods trained on experimental data hold for ML methods trained on Absolut!-generated data. The Absolut! framework has the potential to enable real-world relevant development and benchmarking of ML strategies for biotherapeutics design.

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.

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.

Apyrase-mediated amplification of secretory IgA promotes intestinal homeostasis

Secretory immunoglobulin A (SIgA) interaction with commensal bacteria conditions microbiota composition and function. However, mechanisms regulating reciprocal control of microbiota and SIgA are not defined. Bacteria-derived adenosine triphosphate (ATP) limits T follicular helper (Tfh) cells in the Peyer's patches (PPs) via P2X7 receptor (P2X7R) and thereby SIgA generation. Here we show that hydrolysis of extracellular ATP (eATP) by apyrase results in amplification of the SIgA repertoire. The enhanced breadth of SIgA in mice colonized with apyrase-releasing Escherichia coli influences topographical distribution of bacteria and expression of genes involved in metabolic versus immune functions in the intestinal epithelium. SIgA-mediated conditioning of bacteria and enterocyte function is reflected by differences in nutrient absorption in mice colonized with apyrase-expressing bacteria. Apyrase-induced SIgA improves intestinal homeostasis and attenuates barrier impairment and susceptibility to infection by enteric pathogens in antibiotic-induced dysbiosis. Therefore, amplification of SIgA by apyrase can be leveraged to restore intestinal fitness in dysbiotic conditions.

Benchmarking computational methods for B-cell receptor reconstruction from single-cell RNA-seq data

Multiple methods have recently been developed to reconstruct full-length B-cell receptors (BCRs) from single-cell RNA sequencing (scRNA-seq) data. This need emerged from the expansion of scRNA-seq techniques, the increasing interest in antibody-based drug development and the importance of BCR repertoire changes in cancer and autoimmune disease progression. However, a comprehensive assessment of performance-influencing factors such as the sequencing depth, read length or number of somatic hypermutations (SHMs) as well as guidance regarding the choice of methodology is still lacking. In this work, we evaluated the ability of six available methods to reconstruct full-length BCRs using one simulated and three experimental SMART-seq datasets. In addition, we validated that the BCRs assembled in silico recognize their intended targets when expressed as monoclonal antibodies. We observed that methods such as BALDR, BASIC and BRACER showed the best overall performance across the tested datasets and conditions, whereas only BASIC demonstrated acceptable results on very short read libraries. Furthermore, the de novo assembly-based methods BRACER and BALDR were the most accurate in reconstructing BCRs harboring different degrees of SHMs in the variable domain, while TRUST4, MiXCR and BASIC were the fastest. Finally, we propose guidelines to select the best method based on the given data characteristics.

Development and application of a next-generation sequencing protocol and bioinformatics pipeline for the comprehensive analysis of the canine immunoglobulin repertoire

Profiling the adaptive immune repertoire using next generation sequencing (NGS) has become common in human medicine, showing promise in characterizing clonal expansion of B cell clones through analysis of B cell receptors (BCRs) in patients with lymphoid malignancies. In contrast, most work evaluating BCR repertoires in dogs has employed traditional PCR-based approaches analyzing the IGH locus only. The objectives of this study were to: (1) describe a novel NGS protocol to evaluate canine BCRs; (2) develop a bioinformatics pipeline for processing canine BCR sequencing data; and (3) apply these methods to derive insights into BCR repertoires of healthy dogs and dogs undergoing treatment for B-cell lymphoma. RNA from peripheral blood mononuclear cells of healthy dogs (n = 25) and dogs newly diagnosed with intermediate-to-large B-cell lymphoma (n = 18) with intent to pursue chemotherapy was isolated, converted into cDNA and sequenced by NGS. The BCR repertoires were identified and quantified using a novel analysis pipeline. The IGK repertoires of the healthy dogs were far less diverse compared to IGL which, as with IGH, was highly diverse. Strong biases at key positions within the CDR3 sequence were identified within the healthy dog BCR repertoire. For a subset of the dogs with B-cell lymphoma, clonal expansion of specific IGH sequences pre-treatment and reduction post-treatment was observed. The degree of expansion and reduction correlated with the clinical outcome in this subset. Future studies employing these techniques may improve disease monitoring, provide earlier recognition of disease progression, and ultimately lead to more targeted therapeutics.

A BALB/c IGHV Reference Set, Defined by Haplotype Analysis of Long-Read VDJ-C Sequences From F1 (BALB/c x C57BL/6) Mice

The immunoglobulin genes of inbred mouse strains that are commonly used in models of antibody-mediated human diseases are poorly characterized. This compromises data analysis. To infer the immunoglobulin genes of BALB/c mice, we used long-read SMRT sequencing to amplify VDJ-C sequences from F1 (BALB/c x C57BL/6) hybrid animals. Strain variations were identified in the Ighm and Ighg2b genes, and analysis of VDJ rearrangements led to the inference of 278 germline IGHV alleles. 169 alleles are not present in the C57BL/6 genome reference sequence. To establish a set of expressed BALB/c IGHV germline gene sequences, we computationally retrieved IGHV haplotypes from the IgM dataset. Haplotyping led to the confirmation of 162 BALB/c IGHV gene sequences. A musIGHV398 pseudogene variant also appears to be present in the BALB/cByJ substrain, while a functional musIGHV398 gene is highly expressed in the BALB/cJ substrain. Only four of the BALB/c alleles were also observed in the C57BL/6 haplotype. The full set of inferred BALB/c sequences has been used to establish a BALB/c IGHV reference set, hosted at https://ogrdb.airr-community.org. We assessed whether assemblies from the Mouse Genome Project (MGP) are suitable for the determination of the genes of the IGH loci. Only 37 (43.5%) of the 85 confirmed IMGT-named BALB/c IGHV and 33 (42.9%) of the 77 confirmed non-IMGT IGHV were found in a search of the MGP BALB/cJ genome assembly. This suggests that current MGP assemblies are unsuitable for the comprehensive documentation of germline IGHVs and more efforts will be needed to establish strain-specific reference sets.

Profiling the baseline performance and limits of machine learning models for adaptive immune receptor repertoire classification

BACKGROUND

Machine learning (ML) methodology development for the classification of immune states in adaptive immune receptor repertoires (AIRRs) has seen a recent surge of interest. However, so far, there does not exist a systematic evaluation of scenarios where classical ML methods (such as penalized logistic regression) already perform adequately for AIRR classification. This hinders investigative reorientation to those scenarios where method development of more sophisticated ML approaches may be required.

RESULTS

To identify those scenarios where a baseline ML method is able to perform well for AIRR classification, we generated a collection of synthetic AIRR benchmark data sets encompassing a wide range of data set architecture-associated and immune state-associated sequence patterns (signal) complexity. We trained ≈1,700 ML models with varying assumptions regarding immune signal on ≈1,000 data sets with a total of ≈250,000 AIRRs containing ≈46 billion TCRβ CDR3 amino acid sequences, thereby surpassing the sample sizes of current state-of-the-art AIRR-ML setups by two orders of magnitude. We found that L1-penalized logistic regression achieved high prediction accuracy even when the immune signal occurs only in 1 out of 50,000 AIR sequences.

CONCLUSIONS

We provide a reference benchmark to guide new AIRR-ML classification methodology by (i) identifying those scenarios characterized by immune signal and data set complexity, where baseline methods already achieve high prediction accuracy, and (ii) facilitating realistic expectations of the performance of AIRR-ML models given training data set properties and assumptions. Our study serves as a template for defining specialized AIRR benchmark data sets for comprehensive benchmarking of AIRR-ML methods.

Phenotypic determinism and stochasticity in antibody repertoires of clonally expanded plasma cells

B cell clonal selection and expansion from a genetically diverse antibody repertoire guides the immune response to a target antigen. It remains unclear if clonal selection and expansion follow any deterministic rules or are stochastic with regards to phenotypic antibody properties such as antigen-binding, affinity, and epitope specificity. We perform the in-depth genotypic and phenotypic characterization of antibody repertoires following immunization in mice. We identify the degree to which clonal expansion is driven by antibody binding, affinity, and epitope specificity and as such may provide greater insight into vaccine-induced immunity.

The capacity of humoral B cell-mediated immunity to effectively respond to and protect against pathogenic infections is largely driven by the presence of a diverse repertoire of polyclonal antibodies in the serum, which are produced by plasma cells (PCs). Recent studies have started to reveal the balance between deterministic mechanisms and stochasticity of antibody repertoires on a genotypic level (i.e., clonal diversity, somatic hypermutation, and germline gene usage). However, it remains unclear if clonal selection and expansion of PCs follow any deterministic rules or are stochastic with regards to phenotypic antibody properties (i.e., antigen-binding, affinity, and epitope specificity). Here, we report on the in-depth genotypic and phenotypic characterization of clonally expanded PC antibody repertoires following protein immunization. We find that clonal expansion drives antigen specificity of the most expanded clones (top ∼10), whereas among the rest of the clonal repertoire antigen specificity is stochastic. Furthermore, we report both on a polyclonal repertoire and clonal lineage level that antibody-antigen binding affinity does not correlate with clonal expansion or somatic hypermutation. Last, we provide evidence for convergence toward targeting dominant epitopes despite clonal sequence diversity among the most expanded clones. Our results highlight the extent to which clonal expansion can be ascribed to antigen binding, affinity, and epitope specificity, and they have implications for the assessment of effective vaccines.

Clonally Expanded Virus-Specific CD8 T Cells Acquire Diverse Transcriptional Phenotypes During Acute, Chronic, and Latent Infections

CD8+ T cells play a crucial role in the control and resolution of viral infections and can adopt a wide range of phenotypes and effector functions depending on the inflammatory context and the duration and extent of antigen exposure. Similarly, viral infections can exert diverse selective pressures on populations of clonally related T cells. Technical limitations have nevertheless made it challenging to investigate the relationship between clonal selection and transcriptional phenotypes of virus-specific T cells. We therefore performed single-cell T cell receptor (TCR) repertoire and transcriptome sequencing of virus-specific CD8 T cells in murine models of acute, chronic and latent infection. We observed clear infection-specific populations corresponding to memory, effector, exhausted, and inflationary phenotypes. We further uncovered a mouse-specific and polyclonal T cell response, despite all T cells sharing specificity to a single viral epitope, which was accompanied by stereotypic TCR germline gene usage in all three infection types. Persistent antigen exposure during chronic and latent viral infections resulted in a higher proportion of clonally expanded T cells relative to acute infection. We furthermore observed a relationship between transcriptional heterogeneity and clonal expansion for all three infections, with highly expanded clones having distinct transcriptional phenotypes relative to less expanded clones. Together our work relates clonal selection to gene expression in the context of viral infection and further provides a dataset and accompanying software for the immunological community.

Similarity measurements of B cell receptor repertoire in baseline mice showed spectrum convergence of IgM

BACKGROUND

The B cell receptor (BCR) repertoire is highly diverse among individuals. Poor similarity of the spectrum among inbred baseline mice may limit the ability to discriminate true signals from those involving specific experimental factors. The repertoire similarity of the baseline status lacks intensive measurements.

RESULTS

We measured the repertoire similarity of IgH in blood and spleen samples from untreated BALB/c and C57BL/6J mice to investigate the baseline status of the two inbred strains. The antibody pool was stratified by the isotype of IgA, IgG and IgM. Between individuals, the results showed better convergence of CDR3 and clonal lineage profiles in IgM than in IgA and IgG, and better robustness of somatic mutation networks in IgM than in IgA and IgG. It also showed that the CDR3 clonotypes and clonal lineages shared better in the spleen samples than in the blood samples. The animal batch differences were detected in CDR3 evenness, mutated clonotype proportions, and maximal network degrees. A cut-off of 95% identity in the CDR3 nucleotide sequences was suitable for clonal lineage establishment.

CONCLUSIONS

Our findings reveal a natural landscape of BCR repertoire similarities between baseline mice and provide a solid reference for designing studies of mouse BCR repertoires.

SARS-CoV-2 reactive and neutralizing antibodies discovered by single-cell sequencing of plasma cells and mammalian display

Characterization of COVID-19 antibodies has largely focused on memory B cells; however, it is the antibody-secreting plasma cells that are directly responsible for the production of serum antibodies, which play a critical role in resolving SARS-CoV-2 infection. Little is known about the specificity of plasma cells, largely because plasma cells lack surface antibody expression, thereby complicating their screening. Here, we describe a technology pipeline that integrates single-cell antibody repertoire sequencing and mammalian display to interrogate the specificity of plasma cells from 16 convalescent patients. Single-cell sequencing allows us to profile antibody repertoire features and identify expanded clonal lineages. Mammalian display screening is used to reveal that 43 antibodies (of 132 candidates) derived from expanded plasma cell lineages are specific to SARS-CoV-2 antigens, including antibodies with high affinity to the SARS-CoV-2 receptor-binding domain (RBD) that exhibit potent neutralization and broad binding to the RBD of SARS-CoV-2 variants (of concern/interest).

Progress and challenges for the machine learning-based design of fit-for-purpose monoclonal antibodies

Although the therapeutic efficacy and commercial success of monoclonal antibodies (mAbs) are tremendous, the design and discovery of new candidates remain a time and cost-intensive endeavor. In this regard, progress in the generation of data describing antigen binding and developability, computational methodology, and artificial intelligence may pave the way for a new era of in silico on-demand immunotherapeutics design and discovery. Here, we argue that the main necessary machine learning (ML) components for an in silico mAb sequence generator are: understanding of the rules of mAb-antigen binding, capacity to modularly combine mAb design parameters, and algorithms for unconstrained parameter-driven in silico mAb sequence synthesis. We review the current progress toward the realization of these necessary components and discuss the challenges that must be overcome to allow the on-demand ML-based discovery and design of fit-for-purpose mAb therapeutic candidates.

Adaptive Immune Receptor Repertoire (AIRR) Community Guide to Repertoire Analysis

Adaptive immune receptor repertoires (AIRRs) are rich with information that can be mined for insights into the workings of the immune system. Gene usage, CDR3 properties, clonal lineage structure, and sequence diversity are all capable of revealing the dynamic immune response to perturbation by disease, vaccination, or other interventions. Here we focus on a conceptual introduction to the many aspects of repertoire analysis and orient the reader toward the uses and advantages of each. Along the way, we note some of the many software tools that have been developed for these investigations and link the ideas discussed to chapters on methods provided elsewhere in this volume.

In silico proof of principle of machine learning-based antibody design at unconstrained scale

Generative machine learning (ML) has been postulated to become a major driver in the computational design of antigen-specific monoclonal antibodies (mAb). However, efforts to confirm this hypothesis have been hindered by the infeasibility of testing arbitrarily large numbers of antibody sequences for their most critical design parameters: paratope, epitope, affinity, and developability. To address this challenge, we leveraged a lattice-based antibody-antigen binding simulation framework, which incorporates a wide range of physiological antibody-binding parameters. The simulation framework enables the computation of synthetic antibody-antigen 3D-structures, and it functions as an oracle for unrestricted prospective evaluation and benchmarking of antibody design parameters of ML-generated antibody sequences. We found that a deep generative model, trained exclusively on antibody sequence (one dimensional: 1D) data can be used to design conformational (three dimensional: 3D) epitope-specific antibodies, matching, or exceeding the training dataset in affinity and developability parameter value variety. Furthermore, we established a lower threshold of sequence diversity necessary for high-accuracy generative antibody ML and demonstrated that this lower threshold also holds on experimental real-world data. Finally, we show that transfer learning enables the generation of high-affinity antibody sequences from low-N training data. Our work establishes a priori feasibility and the theoretical foundation of high-throughput ML-based mAb design.

Individualized VDJ recombination predisposes the available Ig sequence space

The process of recombination between variable (V), diversity (D), and joining (J) immunoglobulin (Ig) gene segments determines an individual's naive Ig repertoire and, consequently, (auto)antigen recognition. VDJ recombination follows probabilistic rules that can be modeled statistically. So far, it remains unknown whether VDJ recombination rules differ between individuals. If these rules differed, identical (auto)antigen-specific Ig sequences would be generated with individual-specific probabilities, signifying that the available Ig sequence space is individual specific. We devised a sensitivity-tested distance measure that enables inter-individual comparison of VDJ recombination models. We discovered, accounting for several sources of noise as well as allelic variation in Ig sequencing data, that not only unrelated individuals but also human monozygotic twins and even inbred mice possess statistically distinguishable immunoglobulin recombination models. This suggests that, in addition to genetic, there is also nongenetic modulation of VDJ recombination. We demonstrate that population-wide individualized VDJ recombination can result in orders of magnitude of difference in the probability to generate (auto)antigen-specific Ig sequences. Our findings have implications for immune receptor-based individualized medicine approaches relevant to vaccination, infection, and autoimmunity.

Observed Antibody Space: A diverse database of cleaned, annotated, and translated unpaired and paired antibody sequences

The antibody repertoires of individuals and groups have been used to explore disease states, understand vaccine responses, and drive therapeutic development. The arrival of B‐cell receptor repertoire sequencing has enabled researchers to get a snapshot of these antibody repertoires, and as more data are generated, increasingly in‐depth studies are possible. However, most publicly available data only exist as raw FASTQ files, making the data hard to access, process, and compare. The Observed Antibody Space (OAS) database was created in 2018 to offer clean, annotated, and translated repertoire data. In this paper, we describe an update to OAS that has been driven by the increasing volume of data and the appearance of paired (VH/VL) sequence data. OAS is now accessible via a new web server, with standardized search parameters and a new sequence‐based search option. The new database provides both nucleotides and amino acids for every sequence, with additional sequence annotations to make the data Minimal Information about Adaptive Immune Receptor Repertoire compliant, and comments on potential problems with the sequence. OAS now contains 25 new studies, including severe acute respiratory syndrome coronavirus 2 data and paired sequencing data. The new database is accessible at http://opig.stats.ox.ac.uk/webapps/oas/, and all data are freely available for download.

Accelerated antibody discovery targeting the SARS-CoV-2 spike protein for COVID-19 therapeutic potential

Background

Rapid deployment of technologies capable of high-throughput and high-resolution screening is imperative for timely response to viral outbreaks. Risk mitigation in the form of leveraging multiple advanced technologies further increases the likelihood of identifying efficacious treatments in aggressive timelines.

Methods

In this study, we describe two parallel, yet distinct, in vivo approaches for accelerated discovery of antibodies targeting the severe acute respiratory syndrome coronavirus-2 spike protein. Working with human transgenic Alloy-GK mice, we detail a single B-cell discovery workflow to directly interrogate antibodies secreted from plasma cells for binding specificity and ACE2 receptor blocking activity. Additionally, we describe a concurrent accelerated hybridoma-based workflow utilizing a DiversimAb™ mouse model for increased diversity.

Results

The panel of antibodies isolated from both workflows revealed binding to distinct epitopes with both blocking and non-blocking profiles. Sequence analysis of the resulting lead candidates uncovered additional diversity with the opportunity for straightforward engineering and affinity maturation.

Conclusions

By combining in vivo models with advanced integration of screening and selection platforms, lead antibody candidates can be sequenced and fully characterized within one to three months.

An Analysis of the Effects of Spaceflight and Vaccination on Antibody Repertoire Diversity

Ab repertoire diversity plays a critical role in the host's ability to fight pathogens. CDR3 is partially responsible for Ab-Ag binding and is a significant source of diversity in the repertoire. CDR3 diversity is generated during VDJ rearrangement because of gene segment selection, gene segment trimming and splicing, and the addition of nucleotides. We analyzed the Ab repertoire diversity across multiple experiments examining the effects of spaceflight on the Ab repertoire after vaccination. Five datasets from four experiments were analyzed using rank-abundance curves and Shannon indices as measures of diversity. We discovered a trend toward lower diversity as a result of spaceflight but did not find the same decrease in our physiological model of microgravity in either the spleen or bone marrow. However, the bone marrow repertoire showed a reduction in diversity after vaccination. We also detected differences in Shannon indices between experiments and tissues. We did not detect a pattern of CDR3 usage across the experiments. Overall, we were able to find differences in the Ab repertoire diversity across experimental groups and tissues.

Antibody Repertoire Analysis of Tumor-Infiltrating B Cells Reveals Distinct Signatures and Distributions Across Tissues

The role of B cells in the tumor microenvironment (TME) has largely been under investigated, and data regarding the antibody repertoire encoded by B cells in the TME and the adjacent lymphoid organs are scarce. Here, we utilized B cell receptor high-throughput sequencing (BCR-Seq) to profile the antibody repertoire signature of tumor-infiltrating lymphocyte B cells (TIL−Bs) in comparison to B cells from three anatomic compartments in a mouse model of triple-negative breast cancer. We found that TIL-Bs exhibit distinct antibody repertoire measures, including high clonal polarization and elevated somatic hypermutation rates, suggesting a local antigen-driven B-cell response. Importantly, TIL-Bs were highly mutated but non-class switched, suggesting that class-switch recombination may be inhibited in the TME. Tracing the distribution of TIL-B clones across various compartments indicated that they migrate to and from the TME. The data thus suggests that antibody repertoire signatures can serve as indicators for identifying tumor-reactive B cells.

A Single-Cell Atlas of Lymphocyte Adaptive Immune Repertoires and Transcriptomes Reveals Age-Related Differences in Convalescent COVID-19 Patients

COVID-19 disease outcome is highly dependent on adaptive immunity from T and B lymphocytes, which play a critical role in the control, clearance and long-term protection against SARS-CoV-2. To date, there is limited knowledge on the composition of the T and B cell immune receptor repertoires [T cell receptors (TCRs) and B cell receptors (BCRs)] and transcriptomes in convalescent COVID-19 patients of different age groups. Here, we utilize single-cell sequencing (scSeq) of lymphocyte immune repertoires and transcriptomes to quantitatively profile the adaptive immune response in COVID-19 patients of varying age. We discovered highly expanded T and B cells in multiple patients, with the most expanded clonotypes coming from the effector CD8⁺ T cell population. Highly expanded CD8⁺ and CD4⁺ T cell clones show elevated markers of cytotoxicity (CD8: PRF1, GZMH, GNLY; CD4: GZMA), whereas clonally expanded B cells show markers of transition into the plasma cell state and activation across patients. By comparing young and old convalescent COVID-19 patients (mean ages = 31 and 66.8 years, respectively), we found that clonally expanded B cells in young patients were predominantly of the IgA isotype and their BCRs had incurred higher levels of somatic hypermutation than elderly patients. In conclusion, our scSeq analysis defines the adaptive immune repertoire and transcriptome in convalescent COVID-19 patients and shows important age-related differences implicated in immunity against SARS-CoV-2.

Single-cell BCR and transcriptome analysis after influenza infection reveals spatiotemporal dynamics of antigen-specific B cells

B cell responses are critical for antiviral immunity. However, a comprehensive picture of antigen-specific B cell differentiation, clonal proliferation, and dynamics in different organs after infection is lacking. Here, by combining single-cell RNA and B cell receptor (BCR) sequencing of antigen-specific cells in lymph nodes, spleen, and lungs after influenza infection in mice, we identify several germinal center (GC) B cell subpopulations and organ-specific differences that persist over the course of the response. We discover transcriptional differences between memory cells in lungs and lymphoid organs and organ-restricted clonal expansion. Remarkably, we find significant clonal overlap between GC-derived memory and plasma cells. By combining BCR-mutational analyses with monoclonal antibody (mAb) expression and affinity measurements, we find that memory B cells are highly diverse and can be selected from both low- and high-affinity precursors. By linking antigen recognition with transcriptional programming, clonal proliferation, and differentiation, these finding provide important advances in our understanding of antiviral immunity.

Prospective Artificial Intelligence to Dissect the Dengue Immune Response and Discover Therapeutics

Dengue virus (DENV) poses a serious threat to global health as the causative agent of dengue fever. The virus is endemic in more than 128 countries resulting in approximately 390 million infection cases each year. Currently, there is no approved therapeutic for treatment nor a fully efficacious vaccine. The development of therapeutics is confounded and hampered by the complexity of the immune response to DENV, in particular to sequential infection with different DENV serotypes (DENV1-5). Researchers have shown that the DENV envelope (E) antigen is primarily responsible for the interaction and subsequent invasion of host cells for all serotypes and can elicit neutralizing antibodies in humans. The advent of high-throughput sequencing and the rapid advancements in computational analysis of complex data, has provided tools for the deconvolution of the DENV immune response. Several types of complex statistical analyses, machine learning models and complex visualizations can be applied to begin answering questions about the B- and T-cell immune responses to multiple infections, antibody-dependent enhancement, identification of novel therapeutics and advance vaccine research.

The respiratory syncytial virus (RSV) prefusion F-protein functional antibody repertoire in adult healthy donors

Respiratory syncytial virus (RSV) is the leading cause of death from lower respiratory tract infection in infants and children, and is responsible for considerable morbidity and mortality in older adults. Vaccines for pregnant women and elderly which are in phase III clinical studies target people with pre-existing natural immunity against RSV. To investigate the background immunity which will be impacted by vaccination, we single cell-sorted human memory B cells and dissected functional and genetic features of neutralizing antibodies (nAbs) induced by natural infection. Most nAbs recognized both the prefusion and postfusion conformations of the RSV F-protein (cross-binders) while a smaller fraction bound exclusively to the prefusion conformation. Cross-binder nAbs used a wide array of gene rearrangements, while preF-binder nAbs derived mostly from the expansion of B-cell clonotypes from the IGHV1 germline. This latter class of nAbs recognizes an epitope located between Site Ø, Site II, and Site V on the F-protein, identifying an important site of pathogen vulnerability.

Platypus: an open-access software for integrating lymphocyte single-cell immune repertoires with transcriptomes

High-throughput single-cell sequencing (scSeq) technologies are revolutionizing the ability to molecularly profile B and T lymphocytes by offering the opportunity to simultaneously obtain information on adaptive immune receptor repertoires (VDJ repertoires) and transcriptomes. An integrated quantification of immune repertoire parameters, such as germline gene usage, clonal expansion, somatic hypermutation and transcriptional states opens up new possibilities for the high-resolution analysis of lymphocytes and the inference of antigen-specificity. While multiple tools now exist to investigate gene expression profiles from scSeq of transcriptomes, there is a lack of software dedicated to single-cell immune repertoires. Here, we present Platypus, an open-source software platform providing a user-friendly interface to investigate B-cell receptor and T-cell receptor repertoires from scSeq experiments. Platypus provides a framework to automate and ease the analysis of single-cell immune repertoires while also incorporating transcriptional information involving unsupervised clustering, gene expression and gene ontology. To showcase the capabilities of Platypus, we use it to analyze and visualize single-cell immune repertoires and transcriptomes from B and T cells from convalescent COVID-19 patients, revealing unique insight into the repertoire features and transcriptional profiles of clonally expanded lymphocytes. Platypus will expedite progress by facilitating the analysis of single-cell immune repertoire and transcriptome sequencing.