The promise and challenge of high-throughput sequencing of the antibody repertoire

Single-cell 5' RNA sequencing of camelid peripheral B cells provides insights into cellular basis of heavy-chain antibody production

Camelids produce both conventional tetrameric antibodies (Abs) and dimeric heavy-chain antibodies (HCAbs). Although B cells that generate these two types of Abs exhibit distinct B cell receptors (BCRs), whether these two B cell populations differ in their phenotypes and developmental processes remains unclear. Here, we performed single-cell 5' RNA profiling of peripheral blood mononuclear cell samples from Bactrian camels before and after immunization. We characterized the functional subtypes and differentiation trajectories of circulating B cells in camels, and reconstructed single-cell BCR sequences. We found that in contrast to humans, the proportion of T-bet+ B cells was high among camelid peripheral B cells. Several marker genes of human B cell subtypes, including CD27 and IGHD, were expressed at low levels in the corresponding camel B cell subtypes. Camelid B cells expressing variable genes of HACbs (VHH) were widely present in various functional subtypes and showed highly overlapping differentiation trajectories with B cells expressing variable genes of conventional Abs (VH). After immunization, the transcriptional changes in VHH+ and VH+ B cells were largely consistent. Through structure modeling, we identified a variety of scaffold types among the reconstructed VHH sequences. Our study provides insights into the cellular context of HCAb production in camels and lays the foundation for developing single-B cell-based camelid single-domain Ab screening.

Improved detection of tryptic immunoglobulin variable region peptides by chromatographic and gas-phase fractionation techniques

The polyclonal repertoire of circulating antibodies potentially holds valuable information about an individual's humoral immune state. While bottom-up proteomics is well suited for serum proteomics, the vast number of antibodies and dynamic range of serum challenge this analysis. To acquire the serum proteome more comprehensively, we incorporated high-field asymmetric waveform ion-mobility spectrometry (FAIMS) or two-dimensional chromatography into standard trypsin-based bottom-up proteomics. Thereby, the number of variable region (VR)-related spectra increased 1.7-fold with FAIMS and 10-fold with chromatography fractionation. To match antibody VRs to spectra, we combined de novo searching and BLAST alignment. Validation of this approach showed that, as peptide length increased, the de novo accuracy decreased and BLAST performance increased. Through in silico calculations on antibody repository sequences, we determined the uniqueness of tryptic VR peptides and their suitability as antibody surrogate. Approximately one-third of these peptides were unique, and about one-third of all antibodies contained at least one unique peptide.

Immunoglobulin and T cell receptor repertoire changes induced by a prototype vaccine against Chagas disease in naïve rhesus macaques

BACKGROUND

A vaccine against Trypanosoma cruzi, the agent of Chagas disease, would be an excellent additional tool for disease control. A recombinant vaccine based on Tc24 and TSA1 parasite antigens was found to be safe and immunogenic in naïve macaques.

METHODS

We used RNA-sequencing and performed a transcriptomic analysis of PBMC responses to vaccination of naïve macaques after each vaccine dose, to shed light on the immunogenicity of this vaccine and guide the optimization of doses and formulation. We identified differentially expressed genes and pathways and characterized immunoglobulin and T cell receptor repertoires.

RESULTS

RNA-sequencing analysis indicated a clear transcriptomic response of PBMCs after three vaccine doses, with the up-regulation of several immune cell activation pathways and a broad non-polarized immune profile. Analysis of the IgG repertoire showed that it had a rapid turnover with novel IgGs produced following each vaccine dose, while the TCR repertoire presented several persisting clones that were expanded after each vaccine dose.

CONCLUSIONS

These data suggest that three vaccine doses may be needed for optimum immunogenicity and support the further evaluation of the protective efficacy of this vaccine.

The presence of broadly neutralizing anti-SARS-CoV-2 RBD antibodies elicited by primary series and booster dose of COVID-19 vaccine

Antibody-mediated immunity plays a key role in protection against SARS-CoV-2. We characterized B-cell-derived anti-SARS-CoV-2 RBD antibody repertoires from vaccinated and infected individuals and elucidate the mechanism of action of broadly neutralizing antibodies and dissect antibodies at the epitope level. The breadth and clonality of anti-RBD B cell response varies among individuals. The majority of neutralizing antibody clones lose or exhibit reduced activities against Beta, Delta, and Omicron variants. Nevertheless, a portion of anti-RBD antibody clones that develops after a primary series or booster dose of COVID-19 vaccination exhibit broad neutralization against emerging Omicron BA.2, BA.4, BA.5, BQ.1.1, XBB.1.5 and XBB.1.16 variants. These broadly neutralizing antibodies share genetic features including a conserved usage of the IGHV3-53 and 3-9 genes and recognize three clustered epitopes of the RBD, including epitopes that partially overlap the classically defined set identified early in the pandemic. The Fab-RBD crystal and Fab-Spike complex structures corroborate the epitope grouping of antibodies and reveal the detailed binding mode of broadly neutralizing antibodies. Structure-guided mutagenesis improves binding and neutralization potency of antibody with Omicron variants via a single amino-substitution. Together, these results provide an immunological basis for partial protection against severe COVID-19 by the ancestral strain-based vaccine and indicate guidance for next generation monoclonal antibody development and vaccine design.

Accurate prediction of antibody function and structure using bio-inspired antibody language model

In recent decades, antibodies have emerged as indispensable therapeutics for combating diseases, particularly viral infections. However, their development has been hindered by limited structural information and labor-intensive engineering processes. Fortunately, significant advancements in deep learning methods have facilitated the precise prediction of protein structure and function by leveraging co-evolution information from homologous proteins. Despite these advances, predicting the conformation of antibodies remains challenging due to their unique evolution and the high flexibility of their antigen-binding regions. Here, to address this challenge, we present the Bio-inspired Antibody Language Model (BALM). This model is trained on a vast dataset comprising 336 million 40% nonredundant unlabeled antibody sequences, capturing both unique and conserved properties specific to antibodies. Notably, BALM showcases exceptional performance across four antigen-binding prediction tasks. Moreover, we introduce BALMFold, an end-to-end method derived from BALM, capable of swiftly predicting full atomic antibody structures from individual sequences. Remarkably, BALMFold outperforms those well-established methods like AlphaFold2, IgFold, ESMFold and OmegaFold in the antibody benchmark, demonstrating significant potential to advance innovative engineering and streamline therapeutic antibody development by reducing the need for unnecessary trials. The BALMFold structure prediction server is freely available at https://beamlab-sh.com/models/BALMFold.

Rapid generation of human recombinant monoclonal antibodies from antibody-secreting cells using ferrofluid-based technology

Monoclonal antibodies (mAbs) are one of the most important classes of biologics with high therapeutic and diagnostic value, but traditional methods for mAbs generation, such as hybridoma screening and phage display, have limitations, including low efficiency and loss of natural chain pairing. To overcome these challenges, novel single B cell antibody technologies have emerged, but they also have limitations such as in vitro differentiation of memory B cells and expensive cell sorters. In this study, we present a rapid and efficient workflow for obtaining human recombinant monoclonal antibodies directly from single antigen-specific antibody secreting cells (ASCs) in the peripheral blood of convalescent COVID-19 patients using ferrofluid technology. This process allows the identification and expression of recombinant antigen-specific mAbs in less than 10 days, using RT-PCR to generate linear Ig heavy and light chain gene expression cassettes, called "minigenes", for rapid expression of recombinant antibodies without cloning procedures. This approach has several advantages. First, it saves time and resources by eliminating the need for in vitro differentiation. It also allows individual antigen-specific ASCs to be screened for effector function prior to recombinant antibody cloning, enabling the selection of mAbs with desired characteristics and functional activity. In addition, the method allows comprehensive analysis of variable region repertoires in combination with functional assays to evaluate the specificity and function of the generated antigen-specific antibodies. Our approach, which rapidly generates recombinant monoclonal antibodies from single antigen-specific ASCs, could help to identify functional antibodies and deepen our understanding of antibody dynamics in the immune response through combined antibody repertoire sequence analysis and functional reactivity testing.

B-cell characteristics define HCV reinfection outcome

BACKGROUND & AIMS

In individuals highly exposed to HCV, reinfection is common, suggesting that natural development of sterilising immunity is difficult. In those that are reinfected, some will develop a persistent infection, while a small proportion repeatedly clear the virus, suggesting natural protection is possible. The aim of this study was to characterise immune responses associated with rapid natural clearance of HCV reinfection.

METHODS

Broad neutralising antibodies (nAbs) and Envelope 2 (E2)-specific memory B cell (MBC) responses were examined longitudinally in 15 individuals with varied reinfection outcomes.

RESULTS

Broad nAb responses were associated with MBC recall, but not with clearance of reinfection. Strong evidence of antigen imprinting was found, and the B-cell receptor repertoire showed a high level of clonality with ongoing somatic hypermutation of many clones over subsequent reinfection events. Single-cell transcriptomic analyses showed that cleared reinfections featured an activated transcriptomic profile in HCV-specific B cells that rapidly expanded upon reinfection.

CONCLUSIONS

MBC quality, but not necessarily breadth of nAb responses, is important for protection against antigenically diverse variants, which is encouraging for HCV vaccine development.

IMPACT AND IMPLICATIONS

HCV continues to have a major health burden globally. Limitations in the health infrastructure for diagnosis and treatment, as well as high rates of reinfection, indicate that a vaccine that can protect against chronic HCV infection will greatly complement current efforts to eliminate HCV-related disease. With alternative approaches to testing vaccines, such as controlled human inoculation trials under consideration, we desperately need to identify the correlates of immune protection. In this study, in a small but rare cohort of high-risk injecting drug users who were reinfected multiple times, breadth of neutralisation was not associated with ultimate clearance of the reinfection event. Alternatively, characteristics of the HCV-specific B-cell response associated with B-cell proliferation were. This study indicates that humoral responses are important for protection and suggests that for genetically very diverse viruses, such as HCV, it may be beneficial to look beyond just antibodies as correlates of protection.

Innovations and trends in antibody repertoire analysis

Monoclonal antibodies have revolutionized the treatment of human diseases, which has made them the fastest-growing class of therapeutics, with global sales expected to reach $346.6 billion USD by 2028. Advances in antibody engineering and development have led to the creation of increasingly sophisticated antibody-based therapeutics (e.g. bispecific antibodies and chimeric antigen receptor T cells). However, approaches for antibody discovery have remained comparatively grounded in conventional yet reliable in vitro assays. Breakthrough developments in high-throughput single B-cell sequencing and immunoglobulin proteomic serology, however, have enabled the identification of high-affinity antibodies directly from endogenous B cells or circulating immunoglobulin produced in vivo. Moreover, advances in artificial intelligence offer vast potential for antibody discovery and design with large-scale repertoire datasets positioned as the optimal source of training data for such applications. We highlight advances and recent trends in how these technologies are being applied to antibody repertoire analysis.

De novo antibody discovery in human blood from full-length single B cell transcriptomics and matching haplotyped-resolved germline assemblies

Immunoglobulin (IGH, IGK, IGL) loci in the human genome are highly polymorphic regions that encode the building blocks of the light and heavy chain IG proteins that dimerize to form antibodies. The processes of V(D)J recombination and somatic hypermutation in B cells are responsible for creating an enormous reservoir of highly specific antibodies capable of binding a vast array of possible antigens. However, the antibody repertoire is fundamentally limited by the set of variable (V), diversity (D), and joining (J) alleles present in the germline IG loci. To better understand how the germline IG haplotypes contribute to the expressed antibody repertoire, we combined genome sequencing of the germline IG loci with single-cell transcriptome sequencing of B cells from the same donor. Sequencing and assembly of the germline IG loci captured the IGH locus in a single fully-phased contig where the maternal and paternal contributions to the germline V, D, and J repertoire can be fully resolved. The B cells were collected following a measles, mumps, and rubella (MMR) vaccination, resulting in a population of cells that were activated in response to this specific immune challenge. Single-cell, full-length transcriptome sequencing of these B cells resulted in whole transcriptome characterization of each cell, as well as highly-accurate consensus sequences for the somatically rearranged and hypermutated light and heavy chain IG transcripts. A subset of antibodies synthesized based on their consensus heavy and light chain transcript sequences demonstrated binding to measles antigens and neutralization of measles live virus.

The rise of big data: deep sequencing-driven computational methods are transforming the landscape of synthetic antibody design

Synthetic antibodies (Abs) represent a category of artificial proteins capable of closely emulating the functions of natural Abs. Their in vitro production eliminates the need for an immunological response, streamlining the process of Ab discovery, engineering, and development. These artificially engineered Abs offer novel approaches to antigen recognition, paratope site manipulation, and biochemical/biophysical enhancements. As a result, synthetic Abs are fundamentally reshaping conventional methods of Ab production. This mirrors the revolution observed in molecular biology and genomics as a result of deep sequencing, which allows for the swift and cost-effective sequencing of DNA and RNA molecules at scale. Within this framework, deep sequencing has enabled the exploration of whole genomes and transcriptomes, including particular gene segments of interest. Notably, the fusion of synthetic Ab discovery with advanced deep sequencing technologies is redefining the current approaches to Ab design and development. Such combination offers opportunity to exhaustively explore Ab repertoires, fast-tracking the Ab discovery process, and enhancing synthetic Ab engineering. Moreover, advanced computational algorithms have the capacity to effectively mine big data, helping to identify Ab sequence patterns/features hidden within deep sequencing Ab datasets. In this context, these methods can be utilized to predict novel sequence features thereby enabling the successful generation of de novo Ab molecules. Hence, the merging of synthetic Ab design, deep sequencing technologies, and advanced computational models heralds a new chapter in Ab discovery, broadening our comprehension of immunology and streamlining the advancement of biological therapeutics.

Seven-chain adaptive immune receptor repertoire analysis in rheumatoid arthritis reveals novel features associated with disease and clinically relevant phenotypes

BACKGROUND

In rheumatoid arthritis (RA), the activation of T and B cell clones specific for self-antigens leads to the chronic inflammation of the synovium. Here, we perform an in-depth quantitative analysis of the seven chains that comprise the adaptive immune receptor repertoire (AIRR) in RA.

RESULTS

In comparison to controls, we show that RA patients have multiple and strong differences in the B cell receptor repertoire including reduced diversity as well as altered isotype, chain, and segment frequencies. We demonstrate that therapeutic tumor necrosis factor inhibition partially restores this alteration but find a profound difference in the underlying biochemical reactivities between responders and non-responders. Combining the AIRR with HLA typing, we identify the specific T cell receptor repertoire associated with disease risk variants. Integrating these features, we further develop a molecular classifier that shows the utility of the AIRR as a diagnostic tool.

CONCLUSIONS

Simultaneous sequencing of the seven chains of the human AIRR reveals novel features associated with the disease and clinically relevant phenotypes, including response to therapy. These findings show the unique potential of AIRR to address precision medicine in immune-related diseases.

Epitope Mapping of Human Polyclonal Antibodies to the fHbp Antigen of a Neisseria Meningitidis Vaccine by Hydrogen-Deuterium Exchange Mass Spectrometry (HDX-MS)

Antigen-antibody interactions play a key role in the immune response post vaccination and the mechanism of action of antibody-based biopharmaceuticals. 4CMenB is a multicomponent vaccine against Neisseria meningitidis serogroup B in which factor H binding protein (fHbp) is one of the key antigens. In this study, we use hydrogen/deuterium exchange mass spectrometry (HDX-MS) to identify epitopes in fHbp recognized by polyclonal antibodies (pAb) from two human donors (HDs) vaccinated with 4CMenB. Our HDX-MS data reveal several epitopes recognized by the complex mixture of human pAb. Furthermore, we show that the pAb from the two HDs recognize the same epitope regions. Epitope mapping of total pAb and purified fHbp-specific pAb from the same HD reveals that the two antibody samples recognize the same main epitopes, showing that HDX-MS based epitope mapping can, in this case at least, be performed directly using total IgG pAb samples that have not undergone Ab-selective purification. Two monoclonal antibodies (mAb) were previously produced from B-cell repertoire sequences from one of the HDs and used for epitope mapping of fHbp with HDX-MS. The epitopes identified for the pAb from the same HD in this study, overlap with the epitopes recognized by the two individual mAbs. Overall, HDX-MS epitope mapping appears highly suitable for simultaneous identification of epitopes recognized by pAb from human donors and to thus both guide vaccine development and study basic human immunity to pathogens, including viruses.

Building Representation Learning Models for Antibody Comprehension

Antibodies are versatile proteins with both the capacity to bind a broad range of targets and a proven track record as some of the most successful therapeutics. However, the development of novel antibody therapeutics is a lengthy and costly process. It is challenging to predict the functional and biophysical properties of antibodies from their amino acid sequence alone, requiring numerous experiments for full characterization. Machine learning, specifically deep representation learning, has emerged as a family of methods that can complement wet lab approaches and accelerate the overall discovery and engineering process. Here, we review advances in antibody sequence representation learning, and how this has improved antibody structure prediction and facilitated antibody optimization. We discuss challenges in the development and implementation of such models, such as the lack of publicly available, well-curated antibody function data and highlight opportunities for improvement. These and future advances in machine learning for antibody sequences have the potential to increase the success rate in developing new therapeutics, resulting in broader access to transformative medicines and improved patient outcomes.

Antibodies, repertoires and microdevices in antibody discovery and characterization

Therapeutic antibodies are paramount in treating a wide range of diseases, particularly in auto-immunity, inflammation and cancer, and novel antibody candidates recognizing a vast array of novel antigens are needed to expand the usefulness and applications of these powerful molecules. Microdevices play an essential role in this challenging endeavor at various stages since many general requirements of the overall process overlap nicely with the general advantages of microfluidics. Therefore, microfluidic devices are rapidly taking over various steps in the process of new candidate isolation, such as antibody characterization and discovery workflows. Such technologies can allow for vast improvements in time-lines and incorporate conservative antibody stability and characterization assays, but most prominently screenings and functional characterization within integrated workflows due to high throughput and standardized workflows. First, we aim to provide an overview of the challenges of developing new therapeutic candidates, their repertoires and requirements. Afterward, this review focuses on the discovery of antibodies using microfluidic systems, technological aspects of micro devices and small-scale antibody protein characterization and selection, as well as their integration and implementation into antibody discovery workflows. We close with future developments in microfluidic detection and antibody isolation principles and the field in general.

Large-scale integrative analysis of juvenile idiopathic arthritis for new insight into its pathogenesis

BACKGROUND

Juvenile idiopathic arthritis (JIA) is one of the most prevalent rheumatic disorders in children and is classified as an autoimmune disease (AID). While a robust genetic contribution to JIA etiology has been established, the exact pathogenesis remains unclear.

METHODS

To prioritize biologically interpretable susceptibility genes and proteins for JIA, we conducted transcriptome-wide and proteome-wide association studies (TWAS/PWAS). Then, to understand the genetic architecture of JIA, we systematically analyzed single-nucleotide polymorphism (SNP)-based heritability, a signature of natural selection, and polygenicity. Next, we conducted HLA typing using multi-ethnicity RNA sequencing data. Additionally, we examined the T cell receptor (TCR) repertoire at a single-cell level to explore the potential links between immunity and JIA risk.

RESULTS

We have identified 19 TWAS genes and two PWAS proteins associated with JIA risks. Furthermore, we observe that the heritability and cell type enrichment analysis of JIA are enriched in T lymphocytes and HLA regions and that JIA shows higher polygenicity compared to other AIDs. In multi-ancestry HLA typing, B*45:01 is more prevalent in African JIA patients than in European JIA patients, whereas DQA1*01:01, DQA1*03:01, and DRB1*04:01 exhibit a higher frequency in European JIA patients. Using single-cell immune repertoire analysis, we identify clonally expanded T cell subpopulations in JIA patients, including CXCL13+BHLHE40+ TH cells which are significantly associated with JIA risks.

CONCLUSION

Our findings shed new light on the pathogenesis of JIA and provide a strong foundation for future mechanistic studies aimed at uncovering the molecular drivers of JIA.

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.

Structural insights into the unique recognition module between α-synuclein peptide and nanobody

Nanobodies are single-domain fragments of antibodies with comparable specificity and affinity to antibodies. They are emerging as versatile tools in biology due to their relatively small size. Here, we report the crystal structure of a specific nanobody Nbα-syn01, bound to a 14 amino acid long peptide of α-synuclein (αSyn), a 140-residue protein whose aggregation is associated with Parkinson's disease. The complex structure exhibits a unique binding pattern where the αSyn peptide replaces the N-terminal region of nanobody. Recognition is mediated principally by extended main chain interaction of the αSyn peptide and specificity of the interaction lies in the central 48-52 region of αSyn peptide. Structure-guided truncation of Nbα-syn01 shows tighter binding to αSyn peptide and improved inhibition of α-synuclein aggregation. The structure of the truncated complex was subsequently determined and was indistinguishable to full length complex as the full-length form had no visible electron density for the N-terminal end. These findings reveal the molecular basis for a previously unobserved binding mode for nanobody recognition of α-synuclein, providing an explanation for the enhanced binding, and potential for an alternate framework for structure-based protein engineering of nanobodies to develop better diagnostic and therapeutic tools.

B-Cell Receptor Repertoire: Recent Advances in Autoimmune Diseases

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

nf-core/airrflow: an adaptive immune receptor repertoire analysis workflow employing the Immcantation framework

Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) is a valuable experimental tool to study the immune state in health and following immune challenges such as infectious diseases, (auto)immune diseases, and cancer. Several tools have been developed to reconstruct B cell and T cell receptor sequences from AIRR-seq data and infer B and T cell clonal relationships. However, currently available tools offer limited parallelization across samples, scalability or portability to high-performance computing infrastructures. To address this need, we developed nf-core/airrflow, an end-to-end bulk and single-cell AIRR-seq processing workflow which integrates the Immcantation Framework following BCR and TCR sequencing data analysis best practices. The Immcantation Framework is a comprehensive toolset, which allows the processing of bulk and single-cell AIRR-seq data from raw read processing to clonal inference. nf-core/airrflow is written in Nextflow and is part of the nf-core project, which collects community contributed and curated Nextflow workflows for a wide variety of analysis tasks. We assessed the performance of nf-core/airrflow on simulated sequencing data with sequencing errors and show example results with real datasets. To demonstrate the applicability of nf-core/airrflow to the high-throughput processing of large AIRR-seq datasets, we validated and extended previously reported findings of convergent antibody responses to SARS-CoV-2 by analyzing 97 COVID-19 infected individuals and 99 healthy controls, including a mixture of bulk and single-cell sequencing datasets. Using this dataset, we extended the convergence findings to 20 additional subjects, highlighting the applicability of nf-core/airrflow to validate findings in small in-house cohorts with reanalysis of large publicly available AIRR datasets. nf-core/airrflow is available free of charge, under the MIT license on GitHub (https://github.com/nf-core/airrflow). Detailed documentation and example results are available on the nf-core website at (https://nf-co.re/airrflow).

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.

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.

Comprehensive characterization and database construction of immune repertoire in the largest Chinese glioma cohort

Immune receptor repertoire is valuable for developing immunotherapeutic interventions, but remains poorly understood across glioma subtypes including IDH wild type, IDH mutation without 1p/19q codeletion (IDHmut-noncodel) and IDH mutation with 1p/19q codeletion (IDHmut-codel). We assembled over 320,000 TCR/BCR clonotypes from the largest glioma cohort of 913 RNA sequencing samples in the Chinese population, finding that immune repertoire diversity was more prominent in the IDH wild type (the most aggressive glioma). Fewer clonotypes were shared within each glioma subtype, indicating high heterogeneity of the immune repertoire. The TRA-CDR3 was longer in private than in public clonotypes in IDH wild type. CDR3 variable motifs had higher proportions of hydrophobic residues in private than in public clonotypes, suggesting private CDR3 sequences have greater potential for tumor antigen recognition. Finally, we developed GTABdb, a web-based database designed for hosting, exploring, visualizing, and analyzing glioma immune repertoire. Our study will facilitate developing glioma immunotherapy.

Current advancements in B-cell receptor sequencing fast-track the development of synthetic antibodies

Synthetic antibodies (Abs) are a class of engineered proteins designed to mimic the functions of natural Abs. These are produced entirely in vitro, eliminating the need for an immune response. As such, synthetic Abs have transformed the traditional methods of raising Abs. Likewise, deep sequencing technologies have revolutionized genomics and molecular biology. These enable the rapid and cost-effective sequencing of DNA and RNA molecules. They have allowed for accurate and inexpensive analysis of entire genomes and transcriptomes. Notably, via deep sequencing it is now possible to sequence a person's entire B-cell receptor immune repertoire, termed BCR sequencing. This procedure allows for big data explorations of natural Abs associated with an immune response. Importantly, the identified sequences have the ability to improve the design and engineering of synthetic Abs by offering an initial sequence framework for downstream optimizations. Additionally, machine learning algorithms can be introduced to leverage the vast amount of BCR sequencing datasets to rapidly identify patterns hidden in big data to effectively make in silico predictions of antigen selective synthetic Abs. Thus, the convergence of BCR sequencing, machine learning, and synthetic Ab development has effectively promoted a new era in Ab therapeutics. The combination of these technologies is driving rapid advances in precision medicine, diagnostics, and personalized treatments.

Prediction of antigen-responding VHH antibodies by tracking the evolution of antibody along the time course of immunization

Antibody maturation is the central function of the adaptive immune response. This process is driven by the repetitive selection of mutations that increase the affinity toward antigens. We hypothesized that a precise observation of this process by high-throughput sequencing along the time course of immunization will enable us to predict the antibodies reacting to the immunized antigen without any additional in vitro screening. An alpaca was immunized with IgG fragments using multiple antigen injections, and the antibody repertoire development was traced via high-throughput sequencing periodically for months. The sequences were processed into clusters, and the antibodies in the 16 most abundant clusters were generated to determine whether the clusters included antigen-binding antibodies. The sequences of most antigen-responsive clusters resembled those of germline cells in the early stages. These sequences were observed to accumulate significant mutations and also showed a continuous sequence turnover throughout the experimental period. The foregoing characteristics gave us >80% successful prediction of clusters composed of antigen-responding VHHs against IgG fragment. Furthermore, when the prediction method was applied to the data from other alpaca immunized with epidermal growth factor receptor, the success rate exceeded 80% as well, confirming the general applicability of the prediction method. Superior to previous studies, we identified the immune-responsive but very rare clusters or sequences from the immunized alpaca without any empirical screening data.

Direct Comparison of the Hinge-Cleaving Proteases IgdE and BdpK for LC-MS-Based IgG1 Clonal Profiling

Human antibodies are heterogeneous molecules primarily due to clonal sequence variations. Analytical techniques to assess antibody levels quantitatively, such as ELISA, lack the power to resolve abundances at the clonal level. Recently, we introduced an LC-MS-based approach that can distinguish and quantify antibody clones using the mass and retention time of their corresponding Fab-fragments. We used specific hinge-cleaving protease IgdE (FabALACTICA) to release the Fab-fragments from the constant Fc region of the antibody. Here, we explore an alternative IgG1 hinge-cleaving protease, BdpK (FabDELLO), and compare it directly to IgdE for use in IgG1 repertoire profiling. We used IgdE and BdpK in parallel to digest all IgG1s from the same set of plasma samples. Both proteases cleave IgG1 specifically in the hinge, albeit via different mechanisms and at two distinct cleavage sites. Notwithstanding these differences, the Fab fragments generated by IgdE or BdpK produced highly similar clonal repertoires. However, IgdE required ∼16 h of incubation to digest plasma IgG1s, while BdpK required ∼2 h. We authenticated the similarity of the clones by top-down proteomics using electron transfer dissociation. We conclude that BdpK performs very well in digesting polyclonal plasma IgG1s and that neither BdpK nor IgdE displays detectable biases in cleaving IgG1s. We anticipate that BdpK may emerge as the preferred protease for IgG1 hinge-digestion because it offers a shorter digestion time compared to IgdE, an equally specific digestion site, and no bias against any IgG1 present in plasma.

Evaluation of AlphaFold antibody-antigen modeling with implications for improving predictive accuracy

High resolution antibody-antigen structures provide critical insights into immune recognition and can inform therapeutic design. The challenges of experimental structural determination and the diversity of the immune repertoire underscore the necessity of accurate computational tools for modeling antibody-antigen complexes. Initial benchmarking showed that despite overall success in modeling protein-protein complexes, AlphaFold and AlphaFold-Multimer have limited success in modeling antibody-antigen interactions. In this study, we performed a thorough analysis of AlphaFold's antibody-antigen modeling performance on 427 nonredundant antibody-antigen complex structures, identifying useful confidence metrics for predicting model quality, and features of complexes associated with improved modeling success. Notably, we found that the latest version of AlphaFold improves near-native modeling success to over 30%, versus approximately 20% for a previous version, while increased AlphaFold sampling gives approximately 50% success. With this improved success, AlphaFold can generate accurate antibody-antigen models in many cases, while additional training or other optimization may further improve performance.

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.

Immunoglobin and T cell receptor repertoire changes induced by a prototype vaccine against Chagas disease in naïve rhesus macaques

A vaccine against Trypanosoma cruzi, the agent of Chagas disease, would be an excellent additional tool for disease control. A recombinant vaccine based on Tc24 and TSA1 parasite antigens was found to be safe and immunogenic in naïve macaques. Here we performed a transcriptomic analysis of PBMC responses to vaccination, to shed light on the immunogenicity of this vaccine and guide the optimization of doses and formulation. RNA-sequencing analysis indicated a clear transcriptomic response of PBMCs from macaques after three vaccine doses, with the up-regulation of several immune cell activation pathways and a broad non-polarized immune profile. Analysis of the IgG repertoire showed that it had a rapid turnover with novel IgGs produced following each vaccine dose, while the TCR repertoire presented several persisting clones that were expanded after each vaccine dose. These data suggest that three vaccine doses may be needed for optimum immunogenecity and support the further evaluation of the protective efficacy of this vaccine.

Evaluation of AlphaFold Antibody-Antigen Modeling with Implications for Improving Predictive Accuracy

High resolution antibody-antigen structures provide critical insights into immune recognition and can inform therapeutic design. The challenges of experimental structural determination and the diversity of the immune repertoire underscore the necessity of accurate computational tools for modeling antibody-antigen complexes. Initial benchmarking showed that despite overall success in modeling protein-protein complexes, AlphaFold and AlphaFold-Multimer have limited success in modeling antibody-antigen interactions. In this study, we performed a thorough analysis of AlphaFold's antibody-antigen modeling performance on 429 nonredundant antibody-antigen complex structures, identifying useful confidence metrics for predicting model quality, and features of complexes associated with improved modeling success. We show the importance of bound-like component modeling in complex assembly accuracy, and that the current version of AlphaFold improves near-native modeling success to over 30%, versus approximately 20% for a previous version. With this improved success, AlphaFold can generate accurate antibody-antigen models in many cases, while additional training may further improve its performance.

Flow Cytometry: The Next Revolution

Unmasking the subtleties of the immune system requires both a comprehensive knowledge base and the ability to interrogate that system with intimate sensitivity. That task, to a considerable extent, has been handled by an iterative expansion in flow cytometry methods, both in technological capability and also in accompanying advances in informatics. As the field of fluorescence-based cytomics matured, it reached a technological barrier at around 30 parameter analyses, which stalled the field until spectral flow cytometry created a fundamental transformation that will likely lead to the potential of 100 simultaneous parameter analyses within a few years. The simultaneous advance in informatics has now become a watershed moment for the field as it competes with mature systematic approaches such as genomics and proteomics, allowing cytomics to take a seat at the multi-omics table. In addition, recent technological advances try to combine the speed of flow systems with other detection methods, in addition to fluorescence alone, which will make flow-based instruments even more indispensable in any biological laboratory. This paper outlines current approaches in cell analysis and detection methods, discusses traditional and microfluidic sorting approaches as well as next-generation instruments, and provides an early look at future opportunities that are likely to arise.

Immunoglobulin repertoire restriction characterizes the serological responses of patients with predominantly antibody deficiency

BACKGROUND

Predominantly antibody deficiency (PAD) is the most common category of inborn errors of immunity and is underpinned by impaired generation of appropriate antibody diversity and quantity. In the clinic, responses are interrogated by assessment of vaccination responses, which is central to many PAD diagnoses. However, the composition of the generated antibody repertoire is concealed from traditional quantitative measures of serological responses. Leveraging modern mass spectrometry-based proteomics (MS-proteomics), it is possible to elaborate the molecular features of specific antibody repertoires, which may address current limitations of diagnostic vaccinology.

OBJECTIVES

We sought to evaluate serum antibody responses in patients with PAD following vaccination with a neo-antigen (severe acute respiratory syndrome coronavirus-2 vaccination) using MS-proteomics.

METHODS

Following severe acute respiratory syndrome coronavirus-2 vaccination, serological responses in individuals with PAD and healthy controls (HCs) were assessed by anti-S1 subunit ELISA and neutralization assays. Purified anti-S1 subunit IgG and IgM was profiled by MS-proteomics for IGHV subfamily usage and somatic hypermutation analysis.

RESULTS

Twelve patients with PAD who were vaccine-responsive were recruited with 11 matched vaccinated HCs. Neutralization and end point anti-S1 titers were lower in PAD. All subjects with PAD demonstrated restricted anti-S1 IgG antibody repertoires, with usage of <5 IGHV subfamilies (median: 3; range 2-4), compared to ≥5 for the 11 HC subjects (P < .001). IGHV3-7 utilization was far less common in patients with PAD than in HCs (2 of 12 vs 10 of 11; P = .001). Amino acid substitutions due to somatic hypermutation per subfamily did not differ between groups. Anti-S1 IgM was present in 64% and 50% of HC and PAD cohorts, respectively, and did not differ significantly between HCs and patients with PAD.

CONCLUSIONS

This study demonstrates the breadth of anti-S1 antibodies elicited by vaccination at the proteome level and identifies stereotypical restriction of IGHV utilization in the IgG repertoire in patients with PAD compared with HC subjects. Despite uniformly pauci-clonal antibody repertoires some patients with PAD generated potent serological responses, highlighting a possible limitation of traditional serological techniques. These findings suggest that IgG repertoire restriction is a key feature of antibody repertoires in PAD.

Fast, accurate antibody structure prediction from deep learning on massive set of natural antibodies

Antibodies have the capacity to bind a diverse set of antigens, and they have become critical therapeutics and diagnostic molecules. The binding of antibodies is facilitated by a set of six hypervariable loops that are diversified through genetic recombination and mutation. Even with recent advances, accurate structural prediction of these loops remains a challenge. Here, we present IgFold, a fast deep learning method for antibody structure prediction. IgFold consists of a pre-trained language model trained on 558 million natural antibody sequences followed by graph networks that directly predict backbone atom coordinates. IgFold predicts structures of similar or better quality than alternative methods (including AlphaFold) in significantly less time (under 25 s). Accurate structure prediction on this timescale makes possible avenues of investigation that were previously infeasible. As a demonstration of IgFold's capabilities, we predicted structures for 1.4 million paired antibody sequences, providing structural insights to 500-fold more antibodies than have experimentally determined structures.

Complete variable domain sequences of monoclonal antibody light chains identified from untargeted RNA sequencing data

Introduction

Monoclonal antibody light chain proteins secreted by clonal plasma cells cause tissue damage due to amyloid deposition and other mechanisms. The unique protein sequence associated with each case contributes to the diversity of clinical features observed in patients. Extensive work has characterized many light chains associated with multiple myeloma, light chain amyloidosis and other disorders, which we have collected in the publicly accessible database, AL-Base. However, light chain sequence diversity makes it difficult to determine the contribution of specific amino acid changes to pathology. Sequences of light chains associated with multiple myeloma provide a useful comparison to study mechanisms of light chain aggregation, but relatively few monoclonal sequences have been determined. Therefore, we sought to identify complete light chain sequences from existing high throughput sequencing data.

Methods

We developed a computational approach using the MiXCR suite of tools to extract complete rearranged IGVL-IGJL sequences from untargeted RNA sequencing data. This method was applied to whole-transcriptome RNA sequencing data from 766 newly diagnosed patients in the Multiple Myeloma Research Foundation CoMMpass study.

Results

Monoclonal IGVL-IGJL sequences were defined as those where >50% of assigned IGK or IGL reads from each sample mapped to a unique sequence. Clonal light chain sequences were identified in 705/766 samples from the CoMMpass study. Of these, 685 sequences covered the complete IGVL-IGJL region. The identity of the assigned sequences is consistent with their associated clinical data and with partial sequences previously determined from the same cohort of samples. Sequences have been deposited in AL-Base.

Discussion

Our method allows routine identification of clonal antibody sequences from RNA sequencing data collected for gene expression studies. The sequences identified represent, to our knowledge, the largest collection of multiple myeloma-associated light chains reported to date. This work substantially increases the number of monoclonal light chains known to be associated with non-amyloid plasma cell disorders and will facilitate studies of light chain pathology.

Application of germline antibody features to vaccine development, antibody discovery, antibody optimization and disease diagnosis

Although the efficacy and commercial success of vaccines and therapeutic antibodies have been tremendous, designing and discovering new drug candidates remains a labor-, time- and cost-intensive endeavor with high risks. The main challenges of vaccine development are inducing a strong immune response in broad populations and providing effective prevention against a group of highly variable pathogens. Meanwhile, antibody discovery faces several great obstacles, especially the blindness in antibody screening and the unpredictability of the developability and druggability of antibody drugs. These challenges are largely due to poorly understanding of germline antibodies and the antibody responses to pathogen invasions. Thanks to the recent developments in high-throughput sequencing and structural biology, we have gained insight into the germline immunoglobulin (Ig) genes and germline antibodies and then the germline antibody features associated with antigens and disease manifestation. In this review, we firstly outline the broad associations between germline antibodies and antigens. Moreover, we comprehensively review the recent applications of antigen-specific germline antibody features, physicochemical properties-associated germline antibody features, and disease manifestation-associated germline antibody features on vaccine development, antibody discovery, antibody optimization, and disease diagnosis. Lastly, we discuss the bottlenecks and perspectives of current and potential applications of germline antibody features in the biotechnology field.

Artificial intelligence in cancer immunotherapy: Applications in neoantigen recognition, antibody design and immunotherapy response prediction

Cancer immunotherapy is a method of controlling and eliminating tumors by reactivating the body's cancer-immunity cycle and restoring its antitumor immune response. The increased availability of data, combined with advancements in high-performance computing and innovative artificial intelligence (AI) technology, has resulted in a rise in the use of AI in oncology research. State-of-the-art AI models for functional classification and prediction in immunotherapy research are increasingly used to support laboratory-based experiments. This review offers a glimpse of the current AI applications in immunotherapy, including neoantigen recognition, antibody design, and prediction of immunotherapy response. Advancing in this direction will result in more robust predictive models for developing better targets, drugs, and treatments, and these advancements will eventually make their way into the clinical setting, pushing AI forward in the field of precision oncology.

Exploring the sheep (Ovis aries) immunoglobulin repertoire by next generation sequencing

Next-generation sequencing (NGS) has revolutionized the way we determine the antibody repertoires encoded by B cells in the blood or lymphoid organs and transformed our understanding of adaptive immune responses in many species. Sheep (Ovis aries) have been widely used as a host for therapeutic antibody production since the early 1980s, however, little is known about their immune repertoires or immunological processes affecting the antibody generation. The objective of this study was to employ NGS for a comprehensive analysis of immunoglobulin heavy and light chain repertoires in four healthy sheep. We obtained > 90 % complete antibody sequences and nearly 130,000, 48,000 and 218,000 unique CDR3 reads for the heavy chain (IGH), kappa chain (IGK), and lambda chain (IGL) loci, respectively. Consistent with other species, we observed biased usage of germline variable (V), diversity (D) and joining (J) genes in the heavy and kappa loci, but not in the lambda loci. Moreover, the enormous diversity of CDR3 sequences was observed through sequence clustering and convergent recombination. These data will build a foundation for future studies investigating immune repertoires in health and disease as well as contribute to further refinement of ovine-derived therapeutic antibody drugs.

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

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

Counting is almost all you need

The immune memory repertoire encodes the history of present and past infections and immunological attributes of the individual. As such, multiple methods were proposed to use T-cell receptor (TCR) repertoires to detect disease history. We here show that the counting method outperforms two leading algorithms. We then show that the counting can be further improved using a novel attention model to weigh the different TCRs. The attention model is based on the projection of TCRs using a Variational AutoEncoder (VAE). Both counting and attention algorithms predict better than current leading algorithms whether the host had CMV and its HLA alleles. As an intermediate solution between the complex attention model and the very simple counting model, we propose a new Graph Convolutional Network approach that obtains the accuracy of the attention model and the simplicity of the counting model. The code for the models used in the paper is provided at: https://github.com/louzounlab/CountingIsAlmostAllYouNeed.

Bioinformatic Analysis of Natively Paired VH:VL Antibody Repertoires for Antibody Discovery

Next-generation DNA sequencing (NGS) of human antibody repertoires has been extensively implemented to discover novel antibody drugs, to analyze B-cell developmental features, and to investigate antibody responses to infectious diseases and vaccination. Because the antibody repertoire encoded by human B cells is highly diverse, NGS analyses of antibody genes have provided a new window into understanding antibody responses for basic immunology, biopharmaceutical drug discovery, and immunotherapy. However, many antibody discovery protocols analyze the heavy and light chains separately due to the short-read nature of most NGS technologies, whereas paired heavy and light chain data are required for complete antibody characterization. Here, we describe a computational workflow to process millions of paired antibody heavy and light chain DNA sequence reads using the Illumina MiSeq 2x300 NGS platform. In this workflow, we describe raw NGS read processing and initial quality filtering, the annotation and assembly of antibody clonotypes relating to paired heavy and light chain antibody lineages, and the generation of complete heavy+light consensus sequences for the downstream cloning and expression of human antibody proteins.

Analyzing Antibody Repertoire Using Next-Generation Sequencing and Machine Learning

Advances in high-throughput sequencing technologies have enabled comprehensive sequencing of the immune repertoire. Since repertoire analysis can help to explain the relationship between the immune system and diseases, several methods have been developed for repertoire analysis. Here, using simulated and real-world datasets, we describe how to use DeepRC, a method that applies cutting-edge machine learning techniques.

Characterization of the tumor-infiltrating immune repertoire in muscle invasive bladder cancer

Introduction

Muscle-invasive bladder cancer (MIBC) is a heterogeneous disease with several taxonomic molecular subtypes showing different genetic, clinical, and epidemiological profiles. It has been suggested that MIBC-subtypes follow different tumorigenesis pathways playing decisive roles at different stages of tumor development, resulting in distinct tumor microenvironment containing both innate and adaptive immune cells (T and B lymphocytes). We aim to characterize the MIBC tumor microenvironment by analyzing the tumor-infiltrating B and T cell repertoire according to the taxonomic molecular subtypes.

Methods

RNAseq data from 396 MIBC samples included in TCGA were considered. The subtype information was collected from the international consensus taxonomic classification describing six subtypes: Basal/Squamous-like (Ba/Sq), Luminal papillary (LumP), Luminal non-Specify (LumNS), Luminal unstable (LumU), Stroma-rich, and Neuroendocrine-like (NE-like). Using MiXCR, we mapped the RNA read sequences to their respective B-cell receptor (BCR) and T-cell receptor (TCR) clonotypes. To evaluate the BCR and TCR differences among subtypes, we compared diversity measures (richness and diversity) using a Wilcoxon test and we performed a network analysis to characterize the clonal expansion. For the survival analysis stratified by subtypes, Cox regression models adjusted for age, region, and pathological stage were performed.

Results

Overall, we found different patterns of tumor-infiltrating immune repertoire among the different MIBC subtypes. Stroma-rich and Ba/Sq tumors showed the highest BCR and TCR infiltration while LumP showed the lowest. In addition, we observed that the Ba/Sq and Stroma-rich tumors were more clonally expanded than the Luminal subtypes. Moreover, higher TCR richness and diversity were significantly associated with better survival in the Stroma-rich and Ba/Sq subtypes.

Discussion

This study provides evidence that MIBC subtypes present differences in the tumor microenvironment, in particular, the Ba/Sq and the Stroma-rich are related with a higher tumoral-infiltrating immune repertoire, which seems to be translated into better survival. Determining the causes of the different tumoral-infiltrating immune repertoire according to the MIBC molecular subtypes will help to improve our understanding of the disease and the distinct responses to immunotherapy of MIBC.

Applications of High-Throughput DNA Sequencing to Single-Domain Antibody Discovery and Engineering

Next-generation DNA sequencing (NGS) technologies have made it possible to interrogate antibody repertoires to unprecedented depths, typically via sequencing of cDNAs encoding immunoglobulin variable domains. In the absence of heavy-light chain pairing, the variable domains of heavy chain-only antibodies (HCAbs), referred to as single-domain antibodies (sdAbs), are uniquely amenable to NGS analyses. In this chapter, we provide simple and rapid protocols for producing and sequencing multiplexed immunoglobulin variable domain (VHH, VH, or VL) amplicons derived from a variety of sources using the Illumina MiSeq platform. Generation of such amplicon libraries is relatively inexpensive, requiring no specialized equipment and only a limited set of PCR primers. We also present several applications of NGS to sdAb discovery and engineering, including: (1) evaluation of phage-displayed sdAb library sequence diversity and monitoring of panning experiments; (2) identification of sdAbs of predetermined epitope specificity following competitive elution of phage-displayed sdAb libraries; (3) direct selection of B cells expressing antigen-specific, membrane-bound HCAb using antigen-coupled magnetic beads and identification of antigen-specific sdAbs, and (4) affinity maturation of lead sdAbs using tandem phage display selection and NGS. These methods can easily be adapted to other types of proteins and libraries and expand the utility of in vitro display technology.

Immune repertoire profiling for disease pathobiology

Lymphocytes consist of highly heterogeneous populations, each expressing a specific cell surface receptor corresponding to a particular antigen. Lymphocytes are both the cause and regulator of various diseases, including autoimmune/allergic diseases, lifestyle diseases, neurodegenerative diseases, and cancers. Recently, immune repertoire sequencing has attracted much attention because it helps obtain global profiles of the immune receptor sequences of infiltrating T and B cells in specimens. Immune repertoire sequencing not only helps deepen our understanding of the molecular mechanisms of immune-related pathology but also assists in discovering novel therapeutic modalities for diseases, thereby shedding colorful light on otherwise tiny monotonous cells when observed under a microscope. In this review article, we introduce and detail the background and methodology of immune repertoire sequencing and summarize recent scientific achievements in association with human diseases. Future perspectives on this genetic technique in the field of histopathological research will also be discussed.

Characterization of heavy-chain antibody gene repertoires in Bactrian camels

Camelids are the only mammals that can produce functional heavy-chain antibodies (HCAbs). Although HCAbs were discovered over 30 years ago, the antibody gene repertoire of Bactrian camels remains largely underexplored. To characterize the diversity of variable genes of HCAbs (VHHs), germline and rearranged VHH repertoires are constructed. Phylogenetics analysis shows that all camelid VHH genes are derived from a common ancestor and the nucleotide diversity of VHHs is similar across all camelid species. While species-specific hallmark sites are identified, the non-canonical cysteines specific to VHHs are distinct in Bactrian camels and dromedaries compared with alpacas. Though low divergence at the germline repertoire between wild and domestic Bactrian camels, higher expression of VHHs is observed in some wild Bactrian camels than that of domestic ones. This study not only adds our understanding of VHH repertoire diversity across camelids, but also provides useful resources for HCAb engineering.

Homology Modeling of Antibody Variable Regions: Methods and Applications

Adaptive immunity specifically protects us from antigenic challenges. Antibodies are key effector proteins of adaptive immunity, and they are remarkable in their ability to recognize a virtually limitless number of antigens. Fragment variable (F_V), the antigen-binding region of antibodies, can be split into two main components, namely, framework and complementarity determining regions. The framework (FR) consists of light-chain framework (FR_L) and heavy-chain framework (FR_H). Similarly, the complementarity determining regions (CDRs) comprises of light-chain CDRs 1-3 (CDRs L1-3) and heavy-chain CDRs 1-3 (CDRs H1-3). While FRs are relatively constant in sequence and structure across diverse antibodies, sequence variation in CDRs leading to differential conformations of CDR loops accounts for the distinct antigenic specificities of diverse antibodies. The conserved structural features in FRs and conformity of CDRs to a limited set of standard conformations allow for the accurate prediction of F_V models using homology modeling techniques. Antibody structure prediction from its amino acid sequence has numerous important applications including prediction of antibody-antigen interaction interfaces and redesign of therapeutically and biotechnologically useful antibodies with improved affinity. This chapter summarizes the current practices employed in the successful homology modeling of antibody variable regions and the potential applications of the generated homology models.

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.

Multi-state modeling of antibody-antigen complexes with SAXS profiles and deep-learning models

Antibodies are an established class of human therapeutics. Epitope characterization is an important part of therapeutic antibody discovery. However, structural characterization of antibody-antigen complexes remains challenging. On the one hand, X-ray crystallography or cryo-electron microscopy provide atomic resolution characterization of the epitope, but the data collection process is typically long and the success rate is low. On the other hand, computational methods for modeling antibody-antigen structures from the individual components frequently suffer from a high false positive rate, rarely resulting in a unique solution. Recent deep learning models for structure prediction are also successful in predicting protein-protein complexes. However, they do not perform well for antibody-antigen complexes. Small Angle X-ray Scattering (SAXS) is a reliable technique for rapid structural characterization of protein samples in solution albeit at low resolution. Here, we present an integrative approach for modeling antigen-antibody complexes using the antibody sequence, antigen structure, and experimentally determined SAXS profiles of the antibody, antigen, and the complex. The method models antibody structures using a novel deep-learning approach, NanoNet. The structures of the antibodies and antigens are represented using multiple 3D conformations to account for compositional and conformational heterogeneity of the protein samples that are used to collect the SAXS data. The complexes are predicted by integrating the SAXS profiles with scoring functions for protein-protein interfaces that are based on statistical potentials and antibody-specific deep-learning models. We validated the method via application to four Fab:EGFR and one Fab:PCSK9 antibody:antigen complexes with experimentally available SAXS datasets. The integrative approach returns accurate predictions (interface RMSD<4Å) in the top five predictions for four out of five complexes (respective interface RMSD values of 1.95, 2.18, 2.66 and 3.87Å), providing support for the utility of such a computational pipeline for epitope characterization during therapeutic antibody discovery.

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

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

Inference of B cell clonal families using heavy/light chain pairing information

Next generation sequencing of B cell receptor (BCR) repertoires has become a ubiquitous tool for understanding the antibody-mediated immune response: it is now common to have large volumes of sequence data coding for both the heavy and light chain subunits of the BCR. However, until the recent development of high throughput methods of preserving heavy/light chain pairing information, these samples contained no explicit information on which heavy chain sequence pairs with which light chain sequence. One of the first steps in analyzing such BCR repertoire samples is grouping sequences into clonally related families, where each stems from a single rearrangement event. Many methods of accomplishing this have been developed, however, none so far has taken full advantage of the newly-available pairing information. This information can dramatically improve clustering performance, especially for the light chain. The light chain has traditionally been challenging for clonal family inference because of its low diversity and consequent abundance of non-clonal families with indistinguishable naive rearrangements. Here we present a method of incorporating this pairing information into the clustering process in order to arrive at a more accurate partition of the data into clonally related families. We also demonstrate two methods of fixing imperfect pairing information, which may allow for simplified sample preparation and increased sequencing depth. Finally, we describe several other improvements to the partis software package.