Practical guidelines for B-cell receptor repertoire sequencing analysis

Persistent IgG1 clones dominate and personalize the plasma antibody repertoire

Antibodies play a pivotal role in the immune defense and long-term immunity. Yet, while several studies have highlighted the persistence of antigen-specific antibody responses, it is unclear whether this stems from the continuous production of the same clones or recurrent activation of B cells generating new clones. To examine the stability of the human antibody repertoire, we monitored the concentrations of the most abundant IgG1 clones in plasma samples of 11 healthy donors at nine sampling points over a year. During this year, each donor received three doses of a COVID-19 vaccine. Notwithstanding these vaccinations, the concentrations of the most abundant IgG1 clones remained constant. Given the 2- to 3-week half-life of IgG1 molecules in blood, our data suggest that these clones are associated with long-term immunity and do not undergo somatic hypermutation which would imply short-lived plasma cells. Overall, our data suggest that most of the abundant IgG1 clones in plasma are persistently produced by long-lived plasma cells.

Pan-cancer secreted proteome and skeletal muscle regulation: insight from a proteogenomic data-driven knowledge base

Large-scale, pan-cancer analysis is enabled by data driven knowledge bases that link tumor molecular profiles with phenotypes. A debilitating cancer-related phenotype is skeletal muscle loss, or cachexia, which occurs partly from tumor products secreted into circulation. Using the LinkedOmicsKB knowledge base assembled from the Clinical Proteomics Tumor Analysis Consortium proteogenomic analysis, along with catalogs of human secretome proteins, ligand-receptor pairs and molecular signatures, we sought to identify candidate pan-cancer proteins secreted to blood that could regulate skeletal muscle phenotypes in multiple solid cancers. Tumor proteins having significant pan-cancer associations with muscle were referenced against secretome proteins secreted to blood from the Human Protein Atlas, then verified as increased in paired tumor vs. normal tissues in pan-cancer manner. This workflow revealed seven secreted proteins from cancers afflicting kidneys, head and neck, lungs and pancreas that classified as protein-binding activity modulator, extracellular matrix protein or intercellular signaling molecule. Concordance of these biomarkers with validated molecular signatures of cachexia and senescence supported relevance to muscle and cachexia disease biology, and high tumor expression of the biomarker set associated with lower overall survival. In this article, we discuss avenues by which skeletal muscle and cachexia may be regulated by these candidate pan-cancer proteins secreted to blood, and conceptualize a strategy that considers them collectively as a biomarker signature with potential for refinement by data analytics and radiogenomics for predictive testing of future risk in a non-invasive, blood-based panel amenable to broad uptake and early management.

Applying phylogenetic methods for species delimitation to distinguish B-cell clonal families

The adaptive immune system generates a diverse array of B-cell receptors through the processes of V(D)J recombination and somatic hypermutation. B-cell receptors that bind to an antigen will undergo clonal expansion, creating a Darwinian evolutionary dynamic within individuals. A key step in studying these dynamics is to identify sequences derived from the same ancestral V(D)J recombination event (i.e. a clonal family). There are a number of widely used methods for accomplishing this task but a major limitation of all of them is that they rely, at least in part, on the ability to map sequences to a germline reference set. This requirement is particularly problematic in non-model systems where we often know little about the germline allelic diversity in the study population. Recognizing that delimiting B-cell clonal families is analogous to delimiting species from single locus data, we propose a novel strategy of reconstructing the phylogenetic tree of all B-cell sequences in a sample and using a popular species delimitation method, multi-rate Poisson Tree Processes (mPTP), to delimit clonal families. Using extensive simulations, we show that not only does this phylogenetically explicit approach perform well for the purpose of delimiting clonal families when no reference allele set is available, it performs similarly to state-of-the-art techniques developed specifically for B-cell data even when we have a complete reference allele set. Additionally, our analysis of an empirical dataset shows that mPTP performs similarly to leading methods in the field. These findings demonstrate the utility of using off-the-shelf phylogenetic techniques for analyzing B-cell clonal dynamics in non-model systems, and suggests that phylogenetic inference techniques may be potentially combined with mapping based approaches for even more robust inferences, even in model systems.

Conserved B cell signaling, activation, and differentiation in porcine jejunal and ileal Peyer's patches despite distinct immune landscapes

Peyer's patches (PPs) are B cell-rich sites of intestinal immune induction, yet PP-associated B cell signaling, activation, and differentiation are poorly defined. Single-cell and spatial transcriptomics were completed to study B cells from porcine jejunum and ileum containing PPs. Intestinal locations had distinct immune landscapes, including more follicular B cells in ileum and increased MHC-II-encoding gene expression in jejunal B cells. Despite distinct landscapes, conserved B cell dynamics were detected across intestinal locations, including B cell signaling to CD4+ macrophages that are putative phagocytic, cytotoxic, effector cells and deduced routes of B cell activation/differentiation, including resting B cells migrating into follicles to replicate/divide or differentiate into antibody-secreting cells residing in intestinal crypts. A six-biomarker panel recapitulated transcriptomics findings of B cell phenotypes, frequencies, and spatial locations via ex vivo and in situ staining. Findings convey conserved B cell dynamics across intestinal locations containing PPs, despite location-specific immune environments. Results establish a benchmark of B cell dynamics for understanding intestinal immune induction important to promoting gut/overall health.

Minimal framework for optimizing vaccination protocols targeting highly mutable pathogens

A persistent public health challenge is identifying immunization schemes that are effective against highly mutable pathogens such as HIV and influenza viruses. To address this, we analyze a simplified model of affinity maturation, the Darwinian evolutionary process B cells undergo during immunization. The vaccination protocol determines the selection forces that steer affinity maturation to generate antibodies. We focus on identifying the optimal selection forces exerted by a generic time-dependent vaccination protocol to maximize the production of broadly neutralizing antibodies (bnAbs) that can protect against a broad spectrum of pathogen strains. The model utilizes a path integral representation and operator approximations within a mean-field limit and provides guiding principles for optimizing time-dependent vaccine-induced selection forces to enhance bnAb generation. We compare our analytical mean-field results with the outcomes of stochastic simulations, and we discuss their similarities and differences.

An unbiased comparison of immunoglobulin sequence aligners

Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) is critical for our understanding of the adaptive immune system's dynamics in health and disease. Reliable analysis of AIRR-seq data depends on accurate rearranged immunoglobulin (Ig) sequence alignment. Various Ig sequence aligners exist, but there is no unified benchmarking standard representing the complexities of AIRR-seq data, obscuring objective comparisons of aligners across tasks. Here, we introduce GenAIRR, a modular simulation framework for generating Ig sequences alongside their ground truths. GenAIRR realistically simulates the intricacies of V(D)J recombination, somatic hypermutation, and an array of sequence corruptions. We comprehensively assessed prominent Ig sequence aligners across various metrics, unveiling unique performance characteristics for each aligner. The GenAIRR-produced datasets, combined with the proposed rigorous evaluation criteria, establish a solid basis for unbiased benchmarking of immunogenetics computational tools. It sets up the ground for further improving the crucial task of Ig sequence alignment, ultimately enhancing our understanding of adaptive immunity.

bcRflow: a Nextflow pipeline for characterizing B cell receptor repertoires from non-targeted transcriptomic data

B cells play a critical role in the adaptive recognition of foreign antigens through diverse receptor generation. While targeted immune sequencing methods are commonly used to profile B cell receptors (BCRs), they have limitations in cost and tissue availability. Analyzing B cell receptor profiling from non-targeted transcriptomics data is a promising alternative, but a systematic pipeline integrating tools for accurate immune repertoire extraction is lacking. Here, we present bcRflow, a Nextflow pipeline designed to characterize BCR repertoires from non-targeted transcriptomics data, with functional modules for alignment, processing, and visualization. bcRflow is a comprehensive, reproducible, and scalable pipeline that can run on high-performance computing clusters, cloud-based computing resources like Amazon Web Services (AWS), the Open OnDemand framework, or even local desktops. bcRflow utilizes institutional configurations provided by nf-core to ensure maximum portability and accessibility. To demonstrate the functionality of the bcRflow pipeline, we analyzed a public dataset of bulk transcriptomic samples from COVID-19 patients and healthy controls. We have shown that bcRflow streamlines the analysis of BCR repertoires from non-targeted transcriptomics data, providing valuable insights into the B cell immune response for biological and clinical research. bcRflow is available at https://github.com/Bioinformatics-Core-at-Childrens/bcRflow.

Antibody gene features associated with binding and functional activity in malaria vaccine-derived human mAbs

The impact of adjuvants on malaria vaccine-induced antibody repertoire is poorly understood. Here, we characterize the impact of two adjuvants, Alhydrogel® and AS01, on antibody clonotype diversity, binding and function, post malaria vaccination. We expressed 132 recombinant anti-Pfs230D1 human monoclonal antibodies (mAbs) from participants immunized with malaria transmission-blocking vaccine Pfs230D1, formulated with either Alhydrogel® or AS01. Anti-Pfs230D1 mAbs generated by Alhydrogel® formulation showed higher binding frequency to Pfs230D1 compared to AS01 formulation, although the frequency of functional mAbs was similar between adjuvant groups. Overall, the AS01 formulation induced anti-Pfs230D1 functional antibodies from a broader array of germline sequences versus the Alhydrogel® formulation. All mAbs using IGHV1-69 gene from the Alhydrogel® cohort bound to recombinant Pfs230D1, but did not block parasite transmission to mosquitoes, similar to the IGHV1-69 mAbs isolated from the AS01 cohort. These findings may help inform vaccine design and adjuvant selection for immunization with Plasmodium antigens.

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

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

Resolving haplotype variation and complex genetic architecture in the human immunoglobulin kappa chain locus in individuals of diverse ancestry

Immunoglobulins (IGs), critical components of the human immune system, are composed of heavy and light protein chains encoded at three genomic loci. The IG Kappa (IGK) chain locus consists of two large, inverted segmental duplications. The complexity of the IG loci has hindered use of standard high-throughput methods for characterizing genetic variation within these regions. To overcome these limitations, we use long-read sequencing to create haplotype-resolved IGK assemblies in an ancestrally diverse cohort (n = 36), representing the first comprehensive description of IGK haplotype variation. We identify extensive locus polymorphism, including novel single nucleotide variants (SNVs) and novel structural variants harboring functional IGKV genes. Among 47 functional IGKV genes, we identify 145 alleles, 67 of which were not previously curated. We report inter-population differences in allele frequencies for 10 IGKV genes, including alleles unique to specific populations within this dataset. We identify haplotypes carrying signatures of gene conversion that associate with SNV enrichment in the IGK distal region, and a haplotype with an inversion spanning the proximal and distal regions. These data provide a critical resource of curated genomic reference information from diverse ancestries, laying a foundation for advancing our understanding of population-level genetic variation in the IGK locus.

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.

Single-cell genomics-based immune and disease monitoring in blood malignancies

Achieving long-term disease control using therapeutic immunomodulation is a long-standing concept with a strong tradition in blood malignancies. Besides allogeneic hematopoietic stem cell transplantation that continues to provide potentially curative treatment for otherwise challenging diagnoses, recent years have seen impressive progress in immunotherapies for leukemias and lymphomas with immune checkpoint blockade, bispecific monoclonal antibodies, and CAR T cell therapies. Despite their success, non-response, relapse, and immune toxicities remain frequent, thus prioritizing the elucidation of the underlying mechanisms and identifying predictive biomarkers. The increasing availability of single-cell genomic tools now provides a system's immunology view to resolve the molecular and cellular mechanisms of immunotherapies at unprecedented resolution. Here, we review recent studies that leverage these technological advancements for tracking immune responses, the emergence of immune resistance, and toxicities. As single-cell immune monitoring tools evolve and become more accessible, we expect their wide adoption for routine clinical applications to catalyze more precise therapeutic steering of personal immune responses.

Evaluating methods for B-cell clonal family assignment

The adaptive immune response relies on a diverse repertoire of B-cell receptors, each of which is characterized by a distinct sequence resulting from VDJ-recombination. Upon binding to an antigen, B-cells undergo clonal expansion and in a process unique to B-cells the overall binding affinity of the repertoire is further enhanced by somatic hypermutations in the receptor sequence. For B-cell repertoires it is therefore particularly important to analyze the dynamics of clonal expansion and patterns of somatic hypermutations and thus it is necessary to group the sequences into distinct clones to determine the number and identity of expanding clonal families responding to an antigen. Multiple methods are currently used to identify clones from sequences, employing distinct approaches to the problem. Until now there has not been an extensive comparison of how well these methods perform under the same conditions. Furthermore, since this is fundamentally a phylogenetics problem, we speculated that the mPTP method, which delimits species based on an analysis of changes in the underlying process of diversification, might perform as well as or better than existing methods. Here we conducted extensive simulations of B-cell repertoires under a diverse set of conditions and studied errors in clonal assignment and in downstream ancestral state reconstruction. We demonstrated that SCOPer-H consistently yielded superior results across parameters. However, this approach relies on a good reference assembly for the germline immunoglobulin genes which is lacking for many species. Using mPTP had lower error rates than tailor-made immunogenetic methods and should therefore be considered by researchers studying antibody evolution in non-model organisms without a reference genome.

Inferring B Cell Phylogenies from Paired H and L Chain BCR Sequences with Dowser

Abs are vital to human immune responses and are composed of genetically variable H and L chains. These structures are initially expressed as BCRs. BCR diversity is shaped through somatic hypermutation and selection during immune responses. This evolutionary process produces B cell clones, cells that descend from a common ancestor but differ by mutations. Phylogenetic trees inferred from BCR sequences can reconstruct the history of mutations within a clone. Until recently, BCR sequencing technologies separated H and L chains, but advancements in single-cell sequencing now pair H and L chains from individual cells. However, it is unclear how these separate genes should be combined to infer B cell phylogenies. In this study, we investigated strategies for using paired H and L chain sequences to build phylogenetic trees. We found that incorporating L chains significantly improved tree accuracy and reproducibility across all methods tested. This improvement was greater than the difference between tree-building methods and persisted even when mixing bulk and single-cell sequencing data. However, we also found that many phylogenetic methods estimated significantly biased branch lengths when some L chains were missing, such as when mixing single-cell and bulk BCR data. This bias was eliminated using maximum likelihood methods with separate branch lengths for H and L chain gene partitions. Thus, we recommend using maximum likelihood methods with separate H and L chain partitions, especially when mixing data types. We implemented these methods in the R package Dowser: https://dowser.readthedocs.io.

Single-cell immune repertoire analysis

Single-cell T cell and B cell antigen receptor-sequencing data analysis can potentially perform in-depth assessments of adaptive immune cells that inform on understanding immune cell development to tracking clonal expansion in disease and therapy. However, it has been extremely challenging to analyze and interpret T cells and B cells and their adaptive immune receptor repertoires at the single-cell level due to not only the complexity of the data but also the underlying biology. In this Review, we delve into the computational breakthroughs that have transformed the analysis of single-cell T cell and B cell antigen receptor-sequencing data.

Guidelines for reproducible analysis of adaptive immune receptor repertoire sequencing data

Enhancing the reproducibility and comprehension of adaptive immune receptor repertoire sequencing (AIRR-seq) data analysis is critical for scientific progress. This study presents guidelines for reproducible AIRR-seq data analysis, and a collection of ready-to-use pipelines with comprehensive documentation. To this end, ten common pipelines were implemented using ViaFoundry, a user-friendly interface for pipeline management and automation. This is accompanied by versioned containers, documentation and archiving capabilities. The automation of pre-processing analysis steps and the ability to modify pipeline parameters according to specific research needs are emphasized. AIRR-seq data analysis is highly sensitive to varying parameters and setups; using the guidelines presented here, the ability to reproduce previously published results is demonstrated. This work promotes transparency, reproducibility, and collaboration in AIRR-seq data analysis, serving as a model for handling and documenting bioinformatics pipelines in other research domains.

Blood transcriptomics analysis offers insights into variant-specific immune response to SARS-CoV-2

Bulk RNA sequencing (RNA-seq) of blood is typically used for gene expression analysis in biomedical research but is still rarely used in clinical practice. In this study, we propose that RNA-seq should be considered a diagnostic tool, as it offers not only insights into aberrant gene expression and splicing but also delivers additional readouts on immune cell type composition as well as B-cell and T-cell receptor (BCR/TCR) repertoires. We demonstrate that RNA-seq offers insights into a patient's immune status via integrative analysis of RNA-seq data from patients infected with various SARS-CoV-2 variants (in total 196 samples with up to 200 million reads sequencing depth). We compare the results of computational cell-type deconvolution methods (e.g., MCP-counter, xCell, EPIC, quanTIseq) to complete blood count data, the current gold standard in clinical practice. We observe varying levels of lymphocyte depletion and significant differences in neutrophil levels between SARS-CoV-2 variants. Additionally, we identify B and T cell receptor (BCR/TCR) sequences using the tools MiXCR and TRUST4 to show that-combined with sequence alignments and BLASTp-they could be used to classify a patient's disease. Finally, we investigated the sequencing depth required for such analyses and concluded that 10 million reads per sample is sufficient. In conclusion, our study reveals that computational cell-type deconvolution and BCR/TCR methods using bulk RNA-seq analyses can supplement missing CBC data and offer insights into immune responses, disease severity, and pathogen-specific immunity, all achievable with a sequencing depth of 10 million reads per sample.

Anchor Clustering for million-scale immune repertoire sequencing data

BACKGROUND

The clustering of immune repertoire data is challenging due to the computational cost associated with a very large number of pairwise sequence comparisons. To overcome this limitation, we developed Anchor Clustering, an unsupervised clustering method designed to identify similar sequences from millions of antigen receptor gene sequences. First, a Point Packing algorithm is used to identify a set of maximally spaced anchor sequences. Then, the genetic distance of the remaining sequences to all anchor sequences is calculated and transformed into distance vectors. Finally, distance vectors are clustered using unsupervised clustering. This process is repeated iteratively until the resulting clusters are small enough so that pairwise distance comparisons can be performed.

RESULTS

Our results demonstrate that Anchor Clustering is faster than existing pairwise comparison clustering methods while providing similar clustering quality. With its flexible, memory-saving strategy, Anchor Clustering is capable of clustering millions of antigen receptor gene sequences in just a few minutes.

CONCLUSIONS

This method enables the meta-analysis of immune-repertoire data from different studies and could contribute to a more comprehensive understanding of the immune repertoire data space.

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.

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.

Historical perspective and future directions: computational science in immuno-oncology

Immuno-oncology holds promise for transforming patient care having achieved durable clinical response rates across a variety of advanced and metastatic cancers. Despite these achievements, only a minority of patients respond to immunotherapy, underscoring the importance of elucidating molecular mechanisms responsible for response and resistance to inform the development and selection of treatments. Breakthroughs in molecular sequencing technologies have led to the generation of an immense amount of genomic and transcriptomic sequencing data that can be mined to uncover complex tumor-immune interactions using computational tools. In this review, we discuss existing and emerging computational methods that contextualize the composition and functional state of the tumor microenvironment, infer the reactivity and clonal dynamics from reconstructed immune cell receptor repertoires, and predict the antigenic landscape for immune cell recognition. We further describe the advantage of multi-omics analyses for capturing multidimensional relationships and artificial intelligence techniques for integrating omics data with histopathological and radiological images to encapsulate patterns of treatment response and tumor-immune biology. Finally, we discuss key challenges impeding their widespread use and clinical application and conclude with future perspectives. We are hopeful that this review will both serve as a guide for prospective researchers seeking to use existing tools for scientific discoveries and inspire the optimization or development of novel tools to enhance precision, ultimately expediting advancements in immunotherapy that improve patient survival and quality of life.

B cell phylogenetics in the single cell era

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

Construction of multiplexed transcriptome NGS libraries of microdissected tissue samples based on combinational DNA barcode microbeads

The combination of single-cell RNA sequencing and microdissection techniques that preserves positional information has become a major tool for spatial transcriptome analyses. However, high costs and time requirements, especially for experiments at the single cell scale, make it challenging for this approach to meet the demand for increased throughput. Therefore, we proposed combinational DNA barcode (CDB)-seq as a medium-throughput, multiplexed approach combining Smart-3SEQ and CDB magnetic microbeads for transcriptome analyses of microdissected tissue samples. We conducted a comprehensive comparison of conditions for CDB microbead preparation and related factors and then applied CDB-seq to RNA extracts, fresh frozen (FF) and formalin-fixed paraffin-embedded (FFPE) mouse brain tissue samples. CDB-seq transcriptomic profiles of tens of microdissected samples could be obtained in a simple, cost-effective way, providing a promising method for future spatial transcriptomics.

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

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

Distinct features of B cell receptors in neuromyelitis optica spectrum disorder among CNS inflammatory demyelinating diseases

BACKGROUND

Neuromyelitis optica spectrum disorder (NMOSD) stands out among CNS inflammatory demyelinating diseases (CIDDs) due to its unique disease characteristics, including severe clinical attacks with extensive lesions and its association with systemic autoimmune diseases. We aimed to investigate whether characteristics of B cell receptors (BCRs) differ between NMOSD and other CIDDs using high-throughput sequencing.

METHODS

From a prospective cohort, we recruited patients with CIDDs and categorized them based on the presence and type of autoantibodies: NMOSD with anti-aquaporin-4 antibodies, myelin oligodendrocyte glycoprotein antibody-associated disease (MOGAD) with anti-myelin oligodendrocyte glycoprotein antibodies, double-seronegative demyelinating disease (DSN), and healthy controls (HCs). The BCR features, including isotype class, clonality, somatic hypermutation (SHM), and the third complementarity-determining region (CDR3) length, were analyzed and compared among the different disease groups.

RESULTS

Blood samples from 33 patients with CIDDs (13 NMOSD, 12 MOGAD, and 8 DSN) and 34 HCs were investigated for BCR sequencing. Patients with NMOSD tended to have more activated BCR features compare to the other disease groups. They showed a lower proportion of unswitched isotypes (IgM and IgD) and a higher proportion of switched isotypes (IgG), increased clonality of BCRs, higher rates of SHM, and shorter lengths of CDR3. Notably, advanced age was identified as a clinical factor associated with these activated BCR features, including increased levels of clonality and SHM rates in the NMOSD group. Conversely, no such clinical factors were found to be associated with activated BCR features in the other CIDD groups.

CONCLUSIONS

NMOSD patients, among those with CIDDs, displayed the most pronounced B cell activation, characterized by higher levels of isotype class switching, clonality, SHM rates, and shorter CDR3 lengths. These findings suggest that B cell-mediated humoral immune responses and characteristics in NMOSD patients are distinct from those observed in the other CIDDs, including MOGAD. Age was identified as a clinical factor associated with BCR activation specifically in NMOSD, implying the significance of persistent B cell activation attributed to anti-aquaporin-4 antibodies, even in the absence of clinical relapses throughout an individual's lifetime.

Inferring B cell phylogenies from paired heavy and light chain BCR sequences with Dowser

Antibodies are vital to human immune responses and are composed of genetically variable heavy and light chains. These structures are initially expressed as B cell receptors (BCRs). BCR diversity is shaped through somatic hypermutation and selection during immune responses. This evolutionary process produces B cell clones, cells that descend from a common ancestor but differ by mutations. Phylogenetic trees inferred from BCR sequences can reconstruct the history of mutations within a clone. Until recently, BCR sequencing technologies separated heavy and light chains, but advancements in single cell sequencing now pair heavy and light chains from individual cells. However, it is unclear how these separate genes should be combined to infer B cell phylogenies. In this study, we investigated strategies for using paired heavy and light chain sequences to build phylogenetic trees. We found incorporating light chains significantly improved tree accuracy and reproducibility across all methods tested. This improvement was greater than the difference between tree building methods and persisted even when mixing bulk and single cell sequencing data. However, we also found that many phylogenetic methods estimated significantly biased branch lengths when some light chains were missing, such as when mixing single cell and bulk BCR data. This bias was eliminated using maximum likelihood methods with separate branch lengths for heavy and light chain gene partitions. Thus, we recommend using maximum likelihood methods with separate heavy and light chain partitions, especially when mixing data types. We implemented these methods in the R package Dowser: https://dowser.readthedocs.io.

IGHV allele similarity clustering improves genotype inference from adaptive immune receptor repertoire sequencing data

In adaptive immune receptor repertoire analysis, determining the germline variable (V) allele associated with each T- and B-cell receptor sequence is a crucial step. This process is highly impacted by allele annotations. Aligning sequences, assigning them to specific germline alleles, and inferring individual genotypes are challenging when the repertoire is highly mutated, or sequence reads do not cover the whole V region. Here, we propose an alternative naming scheme for the V alleles, as well as a novel method to infer individual genotypes. We demonstrate the strengths of the two by comparing their outcomes to other genotype inference methods. We validate the genotype approach with independent genomic long-read data. The naming scheme is compatible with current annotation tools and pipelines. Analysis results can be converted from the proposed naming scheme to the nomenclature determined by the International Union of Immunological Societies (IUIS). Both the naming scheme and the genotype procedure are implemented in a freely available R package (PIgLET https://bitbucket.org/yaarilab/piglet). To allow researchers to further explore the approach on real data and to adapt it for their uses, we also created an interactive website (https://yaarilab.github.io/IGHV_reference_book).

ePlatypus: an ecosystem for computational analysis of immunogenomics data

MOTIVATION

The maturation of systems immunology methodologies requires novel and transparent computational frameworks capable of integrating diverse data modalities in a reproducible manner.

RESULTS

Here, we present the ePlatypus computational immunology ecosystem for immunogenomics data analysis, with a focus on adaptive immune repertoires and single-cell sequencing. ePlatypus is an open-source web-based platform and provides programming tutorials and an integrative database that helps elucidate signatures of B and T cell clonal selection. Furthermore, the ecosystem links novel and established bioinformatics pipelines relevant for single-cell immune repertoires and other aspects of computational immunology such as predicting ligand-receptor interactions, structural modeling, simulations, machine learning, graph theory, pseudotime, spatial transcriptomics, and phylogenetics. The ePlatypus ecosystem helps extract deeper insight in computational immunology and immunogenomics and promote open science.

AVAILABILITY AND IMPLEMENTATION

Platypus code used in this manuscript can be found at github.com/alexyermanos/Platypus.

Antibody gene features associated with binding and functional activity in vaccine-derived human mAbs targeting malaria parasites

Adjuvants have been essential to malaria vaccine development, but their impact on the vaccine-induced antibody repertoire is poorly understood. Here, we used cDNA sequences from antigen-specific single memory B cells to express 132 recombinant human anti-Pfs230 monoclonal antibodies (mAbs). Alhydrogel®-induced mAbs demonstrated higher binding to Pfs230D1, although functional activity was similar between adjuvants. All Alhydrogel® mAbs using IGHV1-69 gene bound to recombinant Pfs230D1, but none blocked parasite transmission to mosquitoes; similarly, no AS01 mAb using IGHV1-69 blocked transmission. Functional mAbs from both Alhydrogel® and AS01 vaccines used IGHV3-21 and IGHV3-30 genes. Antibodies with the longest CDR3 sequences were associated with binding but not functional activity. This study assesses adjuvant effects on antibody clonotype diversity during malaria vaccination.

AIRRSHIP: simulating human B cell receptor repertoire sequences

SUMMARY

Adaptive Immune Receptor Repertoire Sequencing is a rapidly developing field that has advanced understanding of the role of the adaptive immune system in health and disease. Numerous tools have been developed to analyse the complex data produced by this technique but work to compare their accuracy and reliability has been limited. Thorough, systematic assessment of their performance is dependent on the ability to produce high quality simulated datasets with known ground truth. We have developed AIRRSHIP, a flexible and fast Python package that produces synthetic human B cell receptor sequences. AIRRSHIP uses a comprehensive set of reference data to replicate key mechanisms in the immunoglobulin recombination process, with a particular focus on junctional complexity. Repertoires generated by AIRRSHIP are highly similar to published data and all steps in the sequence generation process are recorded. These data can be used to not only determine the accuracy of repertoire analysis tools but can also, by tuning of the large number of user-controllable parameters, give insight into factors that contribute to inaccuracies in results.

AVAILABILITY AND IMPLEMENTATION

AIRRSHIP is implemented in Python. It is available via https://github.com/Cowanlab/airrship and on PyPI at https://pypi.org/project/airrship/. Documentation can be found at https://airrship.readthedocs.io/.

ViCloD, an interactive web tool for visualizing B cell repertoires and analyzing intraclonal diversities: application to human B-cell tumors

High throughput sequencing of adaptive immune receptor repertoire (AIRR-seq) has provided numerous human immunoglobulin (IG) sequences allowing specific B cell receptor (BCR) studies such as the antigen-driven evolution of antibodies (soluble forms of the membrane-bound IG part of the BCR). AIRR-seq data allows researchers to examine intraclonal differences caused primarily by somatic hypermutations in IG genes and affinity maturation. Exploring this essential adaptive immunity process could help elucidate the generation of antibodies with high affinity or broadly neutralizing activities. Retracing their evolutionary history could also clarify how vaccines or pathogen exposition drive the humoral immune response, and unravel the clonal architecture of B cell tumors. Computational methods are necessary for large-scale analysis of AIRR-seq properties. However, there is no efficient and interactive tool for analyzing intraclonal diversity, permitting users to explore adaptive immune receptor repertoires in biological and clinical applications. Here we present ViCloD, a web server for large-scale visual analysis of repertoire clonality and intraclonal diversity. ViCloD uses preprocessed data in the format defined by the Adaptive Immune Receptor Repertoire (AIRR) Community. Then, it performs clonal grouping and evolutionary analyses, producing a collection of useful plots for clonal lineage inspection. The web server presents diverse functionalities, including repertoire navigation, clonal abundance analysis, and intraclonal evolutionary tree reconstruction. Users can download the analyzed data in different table formats and save the generated plots as images. ViCloD is a simple, versatile, and user-friendly tool that can help researchers and clinicians to analyze B cell intraclonal diversity. Moreover, its pipeline is optimized to process hundreds of thousands of sequences within a few minutes, allowing an efficient investigation of large and complex repertoires.

Abalign: a comprehensive multiple sequence alignment platform for B-cell receptor immune repertoires

The utilization of high-throughput sequencing (HTS) for B-cell receptor (BCR) immune repertoire analysis has become widespread in the fields of adaptive immunity and antibody drug development. However, the sheer volume of sequences generated by these experiments presents a challenge in data processing. Specifically, multiple sequence alignment (MSA), a critical aspect of BCR analysis, remains inadequate for handling massive BCR sequencing data and lacks the ability to provide immunoglobulin-specific information. To address this gap, we introduce Abalign, a standalone program specifically designed for ultrafast MSA of BCR/antibody sequences. Benchmark tests demonstrate that Abalign achieves comparable or even better accuracy than state-of-the-art MSA tools, and shows remarkable advantages in terms of speed and memory consumption, reducing the time required for high-throughput analysis from weeks to hours. In addition to its alignment capabilities, Abalign offers a broad range of BCR analysis features, including extracting BCRs, constructing lineage trees, assigning VJ genes, analyzing clonotypes, profiling mutations, and comparing BCR immune repertoires. With its user-friendly graphic interface, Abalign can be easily run on personal computers instead of computing clusters. Overall, Abalign is an easy-to-use and effective tool that enables researchers to analyze massive BCR/antibody sequences, leading to new discoveries in the field of immunoinformatics. The software is freely available at http://cao.labshare.cn/abalign/.

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.

Protocol to determine antibody affinity and concentration in complex solutions using microfluidic antibody affinity profiling

Conventional methods of measuring affinity are limited by artificial immobilization, large sample volumes, and homogeneous solutions. This protocol describes microfluidic antibody affinity profiling on complex human samples in solution to obtain a fingerprint reflecting both affinity and active concentration of the target protein. To illustrate the protocol, we analyze the antibody response in SARS-CoV-2 omicron-naïve samples against different SARS-CoV-2 variants of concern. However, the protocol and the technology are amenable to a broad spectrum of biomedical questions. For complete details on the use and execution of this protocol, please refer to Emmenegger et al. (2022),1 Schneider et al. (2022),2 and Fiedler et al. (2022).3.

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.

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.

A somatic hypermutation-based machine learning model stratifies individuals with Crohn's disease and controls

Crohn's disease (CD) is a chronic relapsing-remitting inflammatory disorder of the gastrointestinal tract that is characterized by altered innate and adaptive immune function. Although massively parallel sequencing studies of the T cell receptor repertoire identified oligoclonal expansion of unique clones, much less is known about the B cell receptor (BCR) repertoire in CD. Here, we present a novel BCR repertoire sequencing data set from ileal biopsies from pediatric patients with CD and controls, and identify CD-specific somatic hypermutation (SHM) patterns, revealed by a machine learning (ML) algorithm trained on BCR repertoire sequences. Moreover, ML classification of a different data set from blood samples of adults with CD versus controls identified that V gene usage, clusters, or mutation frequencies yielded excellent results in classifying the disease (F1 > 90%). In summary, we show that an ML algorithm enables the classification of CD based on unique BCR repertoire features with high accuracy.

Clinical application of immune repertoire sequencing in solid organ transplant

Background

Measurement of T cell receptor (TCR) or B cell receptor (BCR) gene utilization may be valuable in monitoring the dynamic changes in donor-reactive clonal populations following transplantation and enabling adjustment in therapy to avoid the consequences of excess immune suppression or to prevent rejection with contingent graft damage and to indicate the development of tolerance.

Objective

We performed a review of current literature to examine research in immune repertoire sequencing in organ transplantation and to assess the feasibility of this technology for clinical application in immune monitoring.

Methods

We searched MEDLINE and PubMed Central for English-language studies published between 2010 and 2021 that examined T cell/B cell repertoire dynamics upon immune activation. Manual filtering of the search results was performed based on relevancy and predefined inclusion criteria. Data were extracted based on study and methodology characteristics.

Results

Our initial search yielded 1933 articles of which 37 met the inclusion criteria; 16 of these were kidney transplant studies (43%) and 21 were other or general transplantation studies (57%). The predominant method for repertoire characterization was sequencing the CDR3 region of the TCR β chain. Repertoires of transplant recipients were found to have decreased diversity in both rejectors and non-rejectors when compared to healthy controls. Rejectors and those with opportunistic infections were more likely to have clonal expansion in T or B cell populations. Mixed lymphocyte culture followed by TCR sequencing was used in 6 studies to define an alloreactive repertoire and in specialized transplant settings to track tolerance.

Conclusion

Methodological approaches to immune repertoire sequencing are becoming established and offer considerable potential as a novel clinical tool for pre- and post-transplant immune monitoring.

Profiling the B cell immune response elicited by vaccination against the respiratory virus SARS-CoV-2

B cells play a fundamental role in host defenses against viral infections. Profiling the B cell response elicited by SARS-CoV-2 vaccination, including the generation and persistence of antigen-specific memory B cells, is essential for improving the knowledge of vaccine immune responsiveness, beyond the antibody response. mRNA-based vaccines have shown to induce a robust class-switched memory B cell response that persists overtime and is boosted by further vaccine administration, suggesting that memory B cells are critical in driving a recall response upon re-exposure to SARS-CoV-2 antigens. Here, we focus on the role of the B cell response in the context of SARS-CoV-2 vaccination, offering an overview of the different technologies that can be used to identify spike-specific B cells, characterize their phenotype using machine learning approaches, measure their capacity to reactivate following antigen encounter, and tracking the maturation of the B cell receptor antigenic affinity.

Advances in antibody discovery from human BCR repertoires

Antibodies make up an important and growing class of compounds used for the diagnosis or treatment of disease. While traditional antibody discovery utilized immunization of animals to generate lead compounds, technological innovations have made it possible to search for antibodies targeting a given antigen within the repertoires of B cells in humans. Here we group these innovations into four broad categories: cell sorting allows the collection of cells enriched in specificity to one or more antigens; BCR sequencing can be performed on bulk mRNA, genomic DNA or on paired (heavy-light) mRNA; BCR repertoire analysis generally involves clustering BCRs into specificity groups or more in-depth modeling of antibody-antigen interactions, such as antibody-specific epitope predictions; validation of antibody-antigen interactions requires expression of antibodies, followed by antigen binding assays or epitope mapping. Together with innovations in Deep learning these technologies will contribute to the future discovery of diagnostic and therapeutic antibodies directly from humans.

A guide to systems-level immunomics

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

Effects of Prior Infection with SARS-CoV-2 on B Cell Receptor Repertoire Response during Vaccination

Understanding the B cell response to SARS-CoV-2 vaccines is a high priority. High-throughput sequencing of the B cell receptor (BCR) repertoire allows for dynamic characterization of B cell response. Here, we sequenced the BCR repertoire of individuals vaccinated by the Pfizer SARS-CoV-2 mRNA vaccine. We compared BCR repertoires of individuals with previous COVID-19 infection (seropositive) to individuals without previous infection (seronegative). We discovered that vaccine-induced expanded IgG clonotypes had shorter heavy-chain complementarity determining region 3 (HCDR3), and for seropositive individuals, these expanded clonotypes had higher somatic hypermutation (SHM) than seronegative individuals. We uncovered shared clonotypes present in multiple individuals, including 28 clonotypes present across all individuals. These 28 shared clonotypes had higher SHM and shorter HCDR3 lengths compared to the rest of the BCR repertoire. Shared clonotypes were present across both serotypes, indicating convergent evolution due to SARS-CoV-2 vaccination independent of prior viral exposure.

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.

B-Cell Receptor Repertoire Sequencing: Deeper Digging into the Mechanisms and Clinical Aspects of Immune-mediated Diseases

B cells play an essential role in adaptive immunity and are intimately correlated with pleiotropic immune-mediated diseases. Each B cell occupies a unique B cell receptor (BCR), and all BCRs throughout our body form “BCR repertoire.” With the development of sequencing technology and coupled bioinformatics, accumulating evidence indicates that BCR repertoire largely varies under physiological and pathological conditions. Therefore, comprehensive grasp of BCR repertoire will provide new insights into the pathogenesis of immune-mediated diseases and help exploit efficient diagnostic and treatment strategies. In this review, we start with an overview of BCR repertoire and related sequencing technologies and summarize their current applications in immune-mediated diseases. We also underscore the challenges of this emerging field and propose promising future directions in advancing BCR repertoire exploration.

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.

Antibody repertoire sequencing analysis

High-throughput sequencing for B cell receptor (BCR) repertoire provides useful insights for the adaptive immune system. With the continuous development of the BCR-seq technology, many efforts have been made to develop methods for analyzing the ever-increasing BCR repertoire data. In this review, we comprehensively outline different BCR repertoire library preparation protocols and summarize three major steps of BCR-seq data analysis, i. e., V(D)J sequence annotation, clonal phylogenetic inference, and BCR repertoire profiling and mining. Different from other reviews in this field, we emphasize background intuition and the statistical principle of each method to help biologists better understand it. Finally, we discuss data mining problems for BCR-seq data and with a highlight on recently emerging multiple-sample analysis.

Single-cell sequencing: a promising approach for uncovering the mechanisms of tumor metastasis

Single-cell sequencing (SCS) is an emerging high-throughput technology that can be used to study the genomics, transcriptomics, and epigenetics at a single cell level. SCS is widely used in the diagnosis and treatment of various diseases, including cancer. Over the years, SCS has gradually become an effective clinical tool for the exploration of tumor metastasis mechanisms and the development of treatment strategies. Currently, SCS can be used not only to analyze metastasis-related malignant biological characteristics, such as tumor heterogeneity, drug resistance, and microenvironment, but also to construct metastasis-related cell maps for predicting and monitoring the dynamics of metastasis. SCS is also used to identify therapeutic targets related to metastasis as it provides insights into the distribution of tumor cell subsets and gene expression differences between primary and metastatic tumors. Additionally, SCS techniques in combination with artificial intelligence (AI) are used in liquid biopsy to identify circulating tumor cells (CTCs), thereby providing a novel strategy for treating tumor metastasis. In this review, we summarize the potential applications of SCS in the field of tumor metastasis and discuss the prospects and limitations of SCS to provide a theoretical basis for finding therapeutic targets and mechanisms of metastasis.

Moving beyond Titers

Vaccination to prevent and even eliminate disease is amongst the greatest achievements of modern medicine. Opportunities remain in vaccine development to improve protection across the whole population. A next step in vaccine development is the detailed molecular characterization of individual humoral immune responses against a pathogen, especially the rapidly evolving pathogens. New technologies such as sequencing the immune repertoire in response to disease, immunogenomics/vaccinomics, particularly the individual HLA variants, and high-throughput epitope characterization offer new insights into disease protection. Here, we highlight the emerging technologies that could be used to identify variation within the human population, facilitate vaccine discovery, improve vaccine safety and efficacy, and identify mechanisms of generating immunological memory. In today’s vaccine-hesitant climate, these techniques used individually or especially together have the potential to improve vaccine effectiveness and safety and thus vaccine uptake rates. We highlight the importance of using these techniques in combination to understand the humoral immune response as a whole after vaccination to move beyond neutralizing titers as the standard for immunogenicity and vaccine efficacy, especially in clinical trials.

Affinity maturation for an optimal balance between long-term immune coverage and short-term resource constraints

Humoral immunity relies on the mutation and selection of B cells to better recognize pathogens. This affinity maturation process produces cells with diverse recognition capabilities. Examining optimal immune strategies that maximize the long-term immune coverage at a minimal metabolic cost, we show when the immune system should mount a de novo response rather than rely on existing memory cells. Our theory recapitulates known modes of the B cell response, predicts the empirical form of the distribution of clone sizes, and rationalizes as a trade-off between metabolic and immune costs the antigenic imprinting effects that limit the efficacy of vaccines (original antigenic sin). Our predictions provide a framework to interpret experimental results that could be used to inform vaccination strategies.

In order to target threatening pathogens, the adaptive immune system performs a continuous reorganization of its lymphocyte repertoire. Following an immune challenge, the B cell repertoire can evolve cells of increased specificity for the encountered strain. This process of affinity maturation generates a memory pool whose diversity and size remain difficult to predict. We assume that the immune system follows a strategy that maximizes the long-term immune coverage and minimizes the short-term metabolic costs associated with affinity maturation. This strategy is defined as an optimal decision process on a finite dimensional phenotypic space, where a preexisting population of cells is sequentially challenged with a neutrally evolving strain. We show that the low specificity and high diversity of memory B cells—a key experimental result—can be explained as a strategy to protect against pathogens that evolve fast enough to escape highly potent but narrow memory. This plasticity of the repertoire drives the emergence of distinct regimes for the size and diversity of the memory pool, depending on the density of de novo responding cells and on the mutation rate of the strain. The model predicts power-law distributions of clonotype sizes observed in data and rationalizes antigenic imprinting as a strategy to minimize metabolic costs while keeping good immune protection against future strains.