The Diversity and Molecular Evolution of B-Cell Receptors during Infection

The immune system in cardiovascular diseases: from basic mechanisms to therapeutic implications

Immune system plays a crucial role in the physiological and pathological regulation of the cardiovascular system. The exploration history and milestones of immune system in cardiovascular diseases (CVDs) have evolved from the initial discovery of chronic inflammation in atherosclerosis to large-scale clinical studies confirming the importance of anti-inflammatory therapy in treating CVDs. This progress has been facilitated by advancements in various technological approaches, including multi-omics analysis (single-cell sequencing, spatial transcriptome et al.) and significant improvements in immunotherapy techniques such as chimeric antigen receptor (CAR)-T cell therapy. Both innate and adaptive immunity holds a pivotal role in CVDs, involving Toll-like receptor (TLR) signaling pathway, nucleotide-binding oligomerization domain-containing proteins 1 and 2 (NOD1/2) signaling pathway, inflammasome signaling pathway, RNA and DNA sensing signaling pathway, as well as antibody-mediated and complement-dependent systems. Meanwhile, immune responses are simultaneously regulated by multi-level regulations in CVDs, including epigenetics (DNA, RNA, protein) and other key signaling pathways in CVDs, interactions among immune cells, and interactions between immune and cardiac or vascular cells. Remarkably, based on the progress in basic research on immune responses in the cardiovascular system, significant advancements have also been made in pre-clinical and clinical studies of immunotherapy. This review provides an overview of the role of immune system in the cardiovascular system, providing in-depth insights into the physiological and pathological regulation of immune responses in various CVDs, highlighting the impact of multi-level regulation of immune responses in CVDs. Finally, we also discuss pre-clinical and clinical strategies targeting the immune system and translational implications in CVDs.

Profiling the broad antibody diversity of lymphatic filariasis immune antibody repertoire by deep sequencing

Lymphatic filariasis is caused by infections of thread-like filarial worms, namely Wuchereria bancrofti, Brugia Malayi and Brugia timori. However, in-depth analysis of the antibody repertoire against Lymphatic filariasis is lacking. Using high-throughput sequencing of antibody repertoires, immunome analysis of IgG (LG) and IgM (LM) repertoires were studied. Despite significant differences between LG and LM in V(D)J gene usage, IGHV4-34, IGHV6-1, IGHD3-10 and IGHJ4 were preferred in both repertoires. The CDR3 in the LG repertoire showed a longer length than LM. Higher SHM level were observed in LG sequences and presence of oligoclonal sequences indicates the extent of clonal expansion. The prevalence of rare clonotypes in LM repertoire depicts the high clonal diversity when compared to LG repertoire. Monoclonal antibodies against closely related parasitic infections were present within the LG repertoire suggesting that immune repertoires may not be as exclusive and biased against the target infection as initially thought. The characterization of the immune repertoire can provide critical insight into the antibody response patterns in disease state, antibody generation process during infections and future antibody designs.

Autoimmune effector mechanisms associated with a defective immunosuppressive axis in immune thrombocytopenia (ITP)

Immune thrombocytopenia (ITP) is an autoimmune disease characterized by an isolated thrombocytopenia and variable phenotype as some patients suffer no bleeding whilst others have bleeding from mild to severe, which may be fatal. This variability probably reflects the disease's complex pathophysiology; a dysregulated hyperreactive immune effector cell response involving the entire adaptive immune system (e.g. B and T cell subsets) that leads to platelet and megakaryocyte (MK) destruction. It appears that these effector responses are due to a breakdown in immune tolerance, and this is characterized by defects in several immunosuppressive cell types. These include defective T regulatory cells (Tregs), B regulatory cells (Bregs) and Myeloid-derived suppressor cells (MDSC), all of which are all intimately associated with antigen presenting cells (APC) such as dendritic cells (DC). The loss of this immunosuppressive axis allows for the activation of unchecked autoreactive T cells and B cells, leading to the development of autoantibodies and cytotoxic T cells (CTL), which can directly destroy platelets in the periphery and inhibit MK platelet production in the bone marrow (BM). This review will focus on the effector cell mechanisms in ITP and highlight the defective immunosuppressive axis that appears responsible for this platelet-specific immune hyperreactivity.

New generalized metric based on branch length distance to compare B cell lineage trees

The B cell lineage tree encapsulates the successive phases of B cell differentiation and maturation, transitioning from hematopoietic stem cells to mature, antibody-secreting cells within the immune system. Mathematically, this lineage can be conceptualized as an evolutionary tree, where each node represents a distinct stage in B cell development, and the edges reflect the differentiation pathways. To compare these lineage trees, a rigorous mathematical metric is essential. Analyzing B cell lineage trees mathematically and quantifying changes in lineage attributes over time necessitates a comparison methodology capable of accurately assessing and measuring these changes. Addressing the intricacies of multiple B cell lineage tree comparisons, this study introduces a novel metric that enhances the precision of comparative analysis. This metric is formulated on principles of metric theory and evolutionary biology, quantifying the dissimilarities between lineage trees by measuring branch length distance and weight. By providing a framework for systematically classifying lineage trees, this metric facilitates the development of predictive models that are crucial for the creation of targeted immunotherapy and vaccines. To validate the effectiveness of this new metric, synthetic datasets that mimic the complexity and variability of real B cell lineage structures are employed. We demonstrated the ability of the new metric method to accurately capture the evolutionary nuances of B cell lineages.

B cell tolerance and autoimmunity: Lessons from repertoires

Adaptive immune cell function is regulated by a highly diverse receptor recombined from variable germline-encoded segments that can recognize an almost unlimited array of epitopes. While this diversity enables the recognition of any pathogen, it also poses a risk of self-recognition, leading to autoimmunity. Many layers of regulation are present during both the generation and activation of B cells to prevent this phenomenon, although they are evidently imperfect. In recent years, our ability to analyze immune repertoires at scale has drastically increased, both through advances in sequencing and single-cell analyses. Here, we review the current knowledge on B cell repertoire analyses, focusing on their implication for autoimmunity. These studies demonstrate that a failure of tolerance occurs at multiple independent checkpoints in different autoimmune contexts, particularly during B cell maturation, plasmablast differentiation, and within germinal centers. These failures are marked by distinct repertoire features that may be used to identify disease- or patient-specific therapeutic approaches.

Sequence diversity of kappa light chains from patients with AL amyloidosis and multiple myeloma

BACKGROUND

AL amyloidosis (AL) results from the misfolding of immunoglobulin light chains (IG LCs). Aim of this study was to comprehensively analyse kappa LC sequences from AL patients in comparison with multiple myeloma (MM).

OBJECTIVE

We analysed IGKV/IGKJ usage and associated organ tropism and IGKV1/D-33 in terms of mutational analysis and theoretical biochemical properties.

MATERIAL AND METHODS

cDNA and bulk RNA sequencing of the LCs of AL and MM patients.

RESULTS

We studied 41 AL and 83 MM patients showing that IGKV1 was most expressed among kappa AL and MM, with higher frequency in AL (80% vs. 53%, p = .002). IGKV3 was underrepresented in AL (10% vs. 30%, p = .014). IGKJ2 was more commonly used in AL than in MM (39% vs. 29%). Patients with IGKV1/D-33 were associated with heart involvement (75%, p = .024). IGKV1/D-33-segments of AL had a higher mutation count (AL = 12.0 vs. MM = 10.0). FR3 and CDR3 were most frequently mutated in both, with a median mutation count in FR3 being the highest (AL = 4.0; MM = 3.5) and one mutation hotspot (FR3 (83I)) for IGKV1/D-33/IGKJ2 was associated with cardiac involvement.

CONCLUSION

This study confirmed that germline usage has an influence on AL amyloidosis risk and organ involvement.

Single-cell BCR and transcriptome analysis reveals peripheral immune signatures in patients with thyroid-associated ophthalmopathy

Thyroid-associated ophthalmopathy (TAO) is the most prevalent orbital disease in adults caused by an autoimmune disorder, which can lead to disfigurement and vision impairment. Developing effective treatments for this condition presents challenges due to our limited understanding of its underlying immune aberrations. In this study, we profiled the immune components in the peripheral blood of patients with TAO as well as healthy individuals, utilizing single-cell RNA sequencing and B-cell receptor repertoires (BCR) analysis. We observed a significant reduction in the proportions of regulatory B cells (Bregs) and type 2 conventional dendritic cells (DCs) in patients with TAO during the active phase. Conversely, there was a significant increase in the proportion of type 1 DCs. Further analysis of cell differentiation trajectory revealed potential impairment in the transition of B cells towards Breg phenotype during the active phase of TAO. Besides, the activation process of TAO appeared to involve inflammation and immune dysfunction, as indicated by the dynamic changes in the activities of key regulators. The abnormalities in the peripheral immune system, such as the reduced capacity of Bregs to suppress inflammation, were primarily driven by the enhanced interaction among Breg, DCs, and monocytes (i.e., CD22-PTPRC and BTLA-TNFRSF14). Collectively, our findings offer a comprehensive insight into the molecular regulation and cellular reconfiguration during the active phase of TAO at the single-cell level, in order to explore the pathogenesis of TAO and provide new ideas for the future treatment of TAO.

Development, High-Throughput Profiling, and Biopanning of a Large Phage Display Single-Domain Antibody Library

Immunoglobulin G-based monoclonal antibodies (mAbs) have been effective in treating various diseases, but their large molecular size can limit their penetration of tissue and efficacy in multifactorial diseases, necessitating the exploration of alternative forms. In this study, we constructed a phage display library comprising single-domain antibodies (sdAbs; or "VHHs"), known for their small size and remarkable stability, using a total of 1.6 × 109 lymphocytes collected from 20 different alpacas, resulting in approximately 7.16 × 1010 colonies. To assess the quality of the constructed library, next-generation sequencing-based high-throughput profiling was performed, analyzing approximately 5.65 × 106 full-length VHH sequences, revealing 92% uniqueness and confirming the library's diverse composition. Systematic characterization of the library revealed multiple sdAbs with high affinity for three therapeutically relevant antigens. In conclusion, our alpaca sdAb phage display library provides a versatile resource for diagnostics and therapeutics. Furthermore, the library's vast natural VHH antibody repertoire offers insights for generating humanized synthetic sdAb libraries, further advancing sdAb-based therapeutics.

Destabilisation of T cell-dependent humoral immunity in sepsis

Sepsis is a heterogeneous condition defined as life-threatening organ dysfunction caused by a dysregulated host response to infection. For some, sepsis presents as a predominantly suppressive disorder, whilst others experience a pro-inflammatory condition which can culminate in a 'cytokine storm'. Frequently, patients experience signs of concurrent hyper-inflammation and immunosuppression, underpinning the difficulty in directing effective treatment. Although intensive care unit mortality rates have improved in recent years, one-third of discharged patients die within the following year. Half of post-sepsis deaths are due to exacerbation of pre-existing conditions, whilst half are due to complications arising from a deteriorated immune system. It has been suggested that the intense and dysregulated response to infection may induce irreversible metabolic reprogramming in immune cells. As a critical arm of immune protection in vertebrates, alterations to the adaptive immune system can have devastating repercussions. Indeed, a marked depletion of lymphocytes is observed in sepsis, correlating with increased rates of mortality. Such sepsis-induced lymphopenia has profound consequences on how T cells respond to infection but equally on the humoral immune response that is both elicited by B cells and supported by distinct CD4+ T follicular helper (TFH) cell subsets. The immunosuppressive state is further exacerbated by functional impairments to the remaining lymphocyte population, including the presence of cells expressing dysfunctional or exhausted phenotypes. This review will specifically focus on how sepsis destabilises the adaptive immune system, with a closer examination on how B cells and CD4+ TFH cells are affected by sepsis and the corresponding impact on humoral immunity.

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.

Phage display sequencing reveals that genetic, environmental, and intrinsic factors influence variation of human antibody epitope repertoire

Phage-displayed immunoprecipitation sequencing (PhIP-seq) has enabled high-throughput profiling of human antibody repertoires. However, a comprehensive overview of environmental and genetic determinants shaping human adaptive immunity is lacking. In this study, we investigated the effects of genetic, environmental, and intrinsic factors on the variation in human antibody repertoires. We characterized serological antibody repertoires against 344,000 peptides using PhIP-seq libraries from a wide range of microbial and environmental antigens in 1,443 participants from a population cohort. We detected individual-specificity, temporal consistency, and co-housing similarities in antibody repertoires. Genetic analyses showed the involvement of the HLA, IGHV, and FUT2 gene regions in antibody-bound peptide reactivity. Furthermore, we uncovered associations between phenotypic factors (including age, cell counts, sex, smoking behavior, and allergies, among others) and particular antibody-bound peptides. Our results indicate that human antibody epitope repertoires are shaped by both genetics and environmental exposures and highlight specific signatures of distinct phenotypes and genotypes.

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.

B cell M-CLL clones retain selection against replacement mutations in their immunoglobulin gene framework regions

Introduction

Chronic lymphocytic leukemia (CLL) is the most common adult leukemia, accounting for 30–40% of all adult leukemias. The dynamics of B-lymphocyte CLL clones with mutated immunoglobulin heavy chain variable region (IgHV) genes in their tumor (M-CLL) can be studied using mutational lineage trees.

Methods

Here, we used lineage tree-based analyses of somatic hypermutation (SHM) and selection in M-CLL clones, comparing the dominant (presumably malignant) clones of 15 CLL patients to their non-dominant (presumably normal) B cell clones, and to those of healthy control repertoires. This type of analysis, which was never previously published in CLL, yielded the following novel insights.

Results

CLL dominant clones undergo – or retain – more replacement mutations that alter amino acid properties such as charge or hydropathy. Although, as expected, CLL dominant clones undergo weaker selection for replacement mutations in the complementarity determining regions (CDRs) and against replacement mutations in the framework regions (FWRs) than non-dominant clones in the same patients or normal B cell clones in healthy controls, they surprisingly retain some of the latter selection in their FWRs. Finally, using machine learning, we show that even the non-dominant clones in CLL patients differ from healthy control clones in various features, most notably their expression of higher fractions of transition mutations.

Discussion

Overall, CLL seems to be characterized by significant loosening – but not a complete loss – of the selection forces operating on B cell clones, and possibly also by changes in SHM mechanisms.

Immune repertoire and responses to neoadjuvant TCHP therapy in HER2-positive breast cancer

Background

Despite the introduction of trastuzumab, pathologic complete response (pCR) is not attained in approximately 30-40% of Human epithelial growth factor receptor-2-positive breast cancer. Tumor-infiltrating lymphocytes (TIL) have been suggested as a predictive marker of treatment response, albeit not always effective. We investigated the relationship between trastuzumab, docetaxel, carboplatin, and pertuzumab (TCHP) treatment and immune repertoire as a treatment response predictor.

Design

In all, 35 cases were divided into two experimental groups: 10 and 25 cases in the preliminary and main experiments, respectively. In the preliminary experiment, the biopsy tissues before TCHP treatment and the surgical tissues after TCHP treatment were compared. In the main experiment, the biopsy tissues before TCHP treatment were compared according to the TCHP treatment response.

Methods

The T-cell repertoire for TRA, TRB, TRG, and TRD, and B-cell repertoire for immunoglobulin heavy, immunoglobulin kappa, and immunoglobulin lambda were evaluated. Whole transcriptome sequencing was also performed.

Results

In the preliminary experiment, the density and richness of the T-cell receptor (TCR) and B-cell receptor (BCR) repertoires decreased after treatment, regardless of TCHP response. In the main experiment, the Shannon's entropy index, density, and length of CDR3 of the TCR and BCR repertoires did not differ significantly in patients who did and did not achieve pCR. The pCR and non-pCR subgroups according to the level of TILs revealed that the non-pCR/lowTIL group had a higher proportion of low-frequency clones than the pCR/lowTIL group in TRA (non-pCR/lowTIL versus pCR/lowTIL, 0.01-0.1%, 63% versus 45.3%; <0.01%, 32.9% versus 51.8%, p < 0.001) and TRB (non-pCR/lowTIL versus pCR/lowTIL, 0.01-0.1%, 26.5% versus 14.7%; <0.01%, 72.0% versus 84.1%, p < 0.001).

Conclusions

The role of the diversity, richness, and density of the TCR and BCR repertoires as predictive markers for TCHP response was not identified. Compositions of low-frequency clones could be candidates for predictive factors of TCHP response; however, validation studies and further research are necessary.

Reconstructing B cell lineage trees with minimum spanning tree and genotype abundances

B cell receptor (BCR) genes exposed to an antigen undergo somatic hypermutations and Darwinian antigen selection, generating a large BCR-antibody diversity. This process, known as B cell affinity maturation, increases antibody affinity, forming a specific B cell lineage that includes the unmutated ancestor and mutated variants. In a B cell lineage, cells with a higher antigen affinity will undergo clonal expansion, while those with a lower affinity will not proliferate and probably be eliminated. Therefore, cellular (genotype) abundance provides a valuable perspective on the ongoing evolutionary process. Phylogenetic tree inference is often used to reconstruct B cell lineage trees and represents the evolutionary dynamic of BCR affinity maturation. However, such methods should process B-cell population data derived from experimental sampling that might contain different cellular abundances. There are a few phylogenetic methods for tracing the evolutionary events occurring in B cell lineages; best-performing solutions are time-demanding and restricted to analysing a reduced number of sequences, while time-efficient methods do not consider cellular abundances. We propose ClonalTree, a low-complexity and accurate approach to construct B-cell lineage trees that incorporates genotype abundances into minimum spanning tree (MST) algorithms. Using both simulated and experimental data, we demonstrate that ClonalTree outperforms MST-based algorithms and achieves a comparable performance to a method that explores tree-generating space exhaustively. Furthermore, ClonalTree has a lower running time, being more convenient for building B-cell lineage trees from high-throughput BCR sequencing data, mainly in biomedical applications, where a lower computational time is appreciable. It is hundreds to thousands of times faster than exhaustive approaches, enabling the analysis of a large set of sequences within minutes or seconds and without loss of accuracy. The source code is freely available at github.com/julibinho/ClonalTree.

Antibodies Targeting the Cell Wall Induce Protection against Virulent Mycobacterium bovis Infection

Accumulating evidence indicates that antibodies can protect against some intracellular pathogens. Mycobacterium bovis is an intracellular bacterium, and its cell wall (CW) is essential for its virulence and survival. However, the questions of whether antibodies play a protective role in immunity against M. bovis infection and what effects antibodies specific to the CW of M. bovis have still remain unclear. Here, we report that antibodies targeting the CW of an isolated pathogenic M. bovis strain and that of an attenuated bacillus Calmette-Guérin (BCG) strain could induce protection against virulent M. bovis infection in vitro and in vivo. Further research found that the antibody-induced protection was mainly achieved by promoting Fc gamma receptor (FcγR)-mediated phagocytosis, inhibiting bacterial intracellular growth, and enhancing the fusion of phagosomes and lysosomes, and it also depended on T cells for its efficacy. Additionally, we analyzed and characterized the B-cell receptor (BCR) repertoires of CW-immunized mice via next-generation sequencing. CW immunization stimulated BCR changes in the complementarity determining region 3 (CDR3) isotype distribution, gene usage, and somatic hypermutation. Overall, our study validates the idea that antibodies targeting the CW induce protection against virulent M. bovis infection. This study highlights the importance of antibodies targeting the CW in the defense against tuberculosis. IMPORTANCE M. bovis is the causative agent of animal tuberculosis (TB) and human TB. Research on M. bovis is of great public health significance. Currently, TB vaccines are mainly aimed at eliciting protection by enhancement of cell-mediated immunity, and there are few studies on protective antibodies. This is the first report of protective antibodies against M. bovis infection, and the antibodies had both preventive and even therapeutic effects in an M. bovis infection mouse model. Additionally, we reveal the relationship between CDR3 gene diversity and the immune characteristics of the antibodies. These results will provide valuable advice for the rational development of TB vaccines.

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 multi-objective based clustering for inferring BCR clonal lineages from high-throughput B cell repertoire data

The adaptive B cell response is driven by the expansion, somatic hypermutation, and selection of B cell clonal lineages. A high number of clonal lineages in a B cell population indicates a highly diverse repertoire, while clonal size distribution and sequence diversity reflect antigen selective pressure. Identifying clonal lineages is fundamental to many repertoire studies, including repertoire comparisons, clonal tracking, and statistical analysis. Several methods have been developed to group sequences from high-throughput B cell repertoire data. Current methods use clustering algorithms to group clonally-related sequences based on their similarities or distances. Such approaches create groups by optimizing a single objective that typically minimizes intra-clonal distances. However, optimizing several objective functions can be advantageous and boost the algorithm convergence rate. Here we propose MobiLLe, a new method based on multi-objective clustering. Our approach requires V(D)J annotations to obtain the initial groups and iteratively applies two objective functions that optimize cohesion and separation within clonal lineages simultaneously. We show that our method greatly improves clonal lineage grouping on simulated benchmarks with varied mutation rates compared to other tools. When applied to experimental repertoires generated from high-throughput sequencing, its clustering results are comparable to the most performing tools and can reproduce the results of previous publications. The method based on multi-objective clustering can accurately identify clonally-related antibody sequences and presents the lowest running time among state-of-art tools. All these features constitute an attractive option for repertoire analysis, particularly in the clinical context. MobiLLe can potentially help unravel the mechanisms involved in developing and evolving B cell malignancies.

High-throughput sequencing can produce a large set of sequences and has profoundly changed our ability to study immune repertoires, particularly B cell receptor sequences. An important application is the analysis of the clonal lineage composition of B cell populations; it is the starting point of many immune repertoire studies, for instance, to differentiate between healthy individuals and those with lymphoid malignancies or other diseases. Several computational methods have been developed to identify clonal lineages from a set of B cell receptor sequences. Most of them apply clustering algorithms and optimize a single objective function that typically minimizes intra-clonal distances. However, optimizing several objective functions in parallel can benefit and increase the clustering performance and efficiency. We propose MobiLLe, the first multi-objective clonal lineage grouping method, which simultaneously optimizes two objective functions for minimizing intra-clonal diversity and maximizing inter-clonal differences. Our approach greatly improved clonal grouping on simulated benchmarks and performed comparably to the most powerful and recent methods on experimental samples. MobiLLe is computationally more efficient than existing tools and does not require any training process or hyper-parameter optimization. It can easily manage large-scale experimental repertoires, providing useful plots to help researchers detect clonally-related sequences in high-throughput B cell repertoire data.

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.

Transcriptome analysis reveals sex-specific alterations in gonads of green mussel exposed to organophosphorus insecticide triazophos

Triazophos (TP) is a widespread pollutant in aquatic environments. A sex-specific metabolic response in green-lipped mussel Perna viridis to TP exposure was observed in our previous study, and this led us to investigate the mechanisms associated with its toxicity. P. viridis were subjected to chronic exposure (15 days) to TP at 35 μg/L to compare the sex-biased transcriptomic profiles in the gonads of male and female mussels. We identified 632 differentially expressed genes (DEGs) (348 up-regulated and 284 down-regulated) in TP-exposed males, and only 61 DEGs (9 up-regulated and 52 down-regulated) in TP-exposed females. Many DEGs were found to be involved in the nervous, reproductive endocrine, oxidative stress, and immune systems of P. viridis. Additionally, enzymatic activity analysis indicated TP induced neurotoxic effects and oxidative damage to the mussels. Our results demonstrate that the stress response and molecular mechanisms of TP toxicology are different between female and male mussels.

Characterization of B cell receptor H-CDR3 repertoire of spleen in PRV-infected mice

Pseudorabies virus (PRV), also known as suid Alphaherpesvirus 1 (SuHV-1), which is one of the most devastating infectious pathogen of swine industry worldwide. Vaccination is the safest and most effective PRV prevention and control strategy. B cell receptor (BCR) is membrane-bound immunoglobulin located on the surface of B cells capable of specifically binding foreign antigens, which is one of the most important molecules regulating the proliferation and function of B cells. Here, to assess the molecular diversity of BCR H-CDR3 repertoire after different PRV strains infection, we detected the IGHV, IGHD, IGHJ genes usage and CDR3 sequence changes of mice spleen with PRV vaccine strain (Bartha-K61), variant strain (XJ) and mock infection by high-throughput sequencing. We found that PRV-infected groups shared partial BCR sequences, which are most likely to be PRV-specific BCR candidates. However, there were still differences in the IGHV genes usage as well as the combined usage of IGHV and IGHJ genes between the Bartha-K61 strain and XJ strain infection groups. In addition, the CDR3 sequences exhibited large differences in the types and lengths in PRV infection groups. Our study contributes to a better understanding of the host adaptive immune response to PRV infection and provides a theoretical basis for further research on novel and efficient PRV vaccines in the future.

Evaluation of Spike Protein Epitopes by Assessing the Dynamics of Humoral Immune Responses in Moderate COVID-19

The coronavirus disease 2019 (COVID-19) pandemic is caused by a novel coronavirus called severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The spike protein (S) of SARS-CoV-2 is a major target for diagnosis and vaccine development because of its essential role in viral infection and host immunity. Currently, time-dependent responses of humoral immune system against various S protein epitopes are poorly understood. In this study, enzyme-linked immunosorbent assay (ELISA), peptide microarray, and antibody binding epitope mapping (AbMap) techniques were used to systematically analyze the dynamic changes of humoral immune responses against the S protein in a small cohort of moderate COVID-19 patients who were hospitalized for approximately two months after symptom onset. Recombinant truncated S proteins, target S peptides, and random peptides were used as antigens in the analyses. The assays demonstrated the dynamic IgM- and IgG recognition and reactivity against various S protein epitopes with patient-dependent patterns. Comprehensive analysis of epitope distribution along the spike gene sequence and spatial structure of the homotrimer S protein demonstrated that most IgM- and IgG-reactive peptides were clustered into similar genomic regions and were located at accessible domains. Seven S peptides were generally recognized by IgG antibodies derived from serum samples of all COVID-19 patients. The dynamic immune recognition signals from these seven S peptides were comparable to those of the entire S protein or truncated S1 protein. This suggested that the humoral immune system recognized few conserved S protein epitopes in most COVID-19 patients during the entire duration of humoral immune response after symptom onset. Furthermore, in this cohort, individual patients demonstrated stable immune recognition to certain S protein epitopes throughout their hospitalization period. Therefore, the dynamic characteristics of humoral immune responses to S protein have provided valuable information for accurate diagnosis and immunotherapy of COVID-19 patients.

Global Characterization of Peripheral B Cells in Parkinson's Disease by Single-Cell RNA and BCR Sequencing

Immune system plays important roles in the pathogenesis of Parkinson’s disease (PD). However, the role of B cells in this complex disease are still not fully understood. B cells produce antibodies but can also regulate immune responses. In order to decode the relative contribution of peripheral B cell subtypes to the etiology of PD, we performed single cell RNA and BCR sequencing for 10,466 B cells from 8 PD patients and 6 age-matched healthy controls. We observed significant increased memory B cells and significant decreased naïve B cells in PD patients compared to healthy controls. Notably, we also discovered increased IgG and IgA isotypes and more frequent class switch recombination events in PD patients. Moreover, we identified preferential V and J gene segments of B cell receptors in PD patients as the evidence of convergent selection in PD. Finally, we found a marked clonal expanded memory B cell population in PD patients, up-regulating both MHC II genes (HLA-DRB5, HLA-DQA2 and HLA-DPB1) and transcription factor activator protein 1 (AP-1), suggesting that the antigen presentation capacity of B cells was enhanced and B cells were activated in PD patients. Overall, this study conducted a comprehensive analysis of peripheral B cell characteristics of PD patients, which provided novel insights into the humoral immune response in the pathogenesis of PD.

Rapid humoral immune responses are required for recovery from haemorrhagic fever with renal syndrome patients

Haemorrhagic fever with renal syndrome (HFRS) following Hantaan virus (HTNV) infection displays variable clinical signs. Humoral responses elicited during HTNV infections are considered important, however, this process remains poorly understood. Herein, we have investigated the phenotype, temporal dynamics, and characteristics of B-cell receptor (BCR) repertoire in an HFRS cohort. The serological profiles were characterized by a lowered expression level of nucleoprotein (NP)-specific antibody in severe cases. Importantly, B-cell subsets were activated and proliferated within the first two weeks of symptom onset and moderate cases reacted more rapidly. BCR analysis in the recovery phase revealed a dramatic increase in the immunoglobulin gene diversity which was more significantly progressed in moderate infections. In severe cases, B-cell-related transcription was lower with inflammatory sets overactivated. Taken together, these data suggest the clinical signs and disease recovery in HFRS patients were positively impacted by rapid and efficacious humoral responses.

Global characterization of B cell receptor repertoire in COVID-19 patients by single-cell V(D)J sequencing

The world is facing a pandemic of Corona Virus Disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Adaptive immune responses are essential for SARS-CoV-2 virus clearance. Although a large body of studies have been conducted to investigate the immune mechanism in COVID-19 patients, we still lack a comprehensive understanding of the BCR repertoire in patients. In this study, we used the single-cell V(D)J sequencing to characterize the BCR repertoire across convalescent COVID-19 patients. We observed that the BCR diversity was significantly reduced in disease compared with healthy controls. And BCRs tend to skew toward different V gene segments in COVID-19 and healthy controls. The CDR3 sequences of heavy chain in clonal BCRs in patients were more convergent than that in healthy controls. In addition, we discovered increased IgG and IgA isotypes in the disease, including IgG1, IgG3 and IgA1. In all clonal BCRs, IgG isotypes had the most frequent class switch recombination events and the highest somatic hypermutation rate, especially IgG3. Moreover, we found that an IgG3 cluster from different clonal groups had the same IGHV, IGHJ and CDR3 sequences (IGHV4-4-CARLANTNQFYDSSSYLNAMDVW-IGHJ6). Overall, our study provides a comprehensive characterization of the BCR repertoire in COVID-19 patients, which contributes to the understanding of the mechanism for the immune response to SARS-CoV-2 infection.

Immune repertoire: Revealing the "real-time" adaptive immune response in autoimmune diseases

The diversity of the immune repertoire (IR) enables the human immune system to distinguish multifarious antigens (Ags) that humans may encounter throughout life. At the same time, bias or abnormalities in the IR also pay a contribution to the pathogenesis of autoimmune diseases. Rapid advancements in high-throughput sequencing (HTS) technology have ushered in a new era of immune studies, revealing novel molecules and pathways that might result in autoimmunity. In the field of IR, HTS can monitor the immune response status and identify disease-specific immune repertoires. In this review, we summarize updated progress on the mechanisms of the IR and current related studies on four autoimmune diseases, particularly focusing on systemic lupus erythematosus (SLE). These autoimmune diseases can exhibit slightly or significantly skewed IRs and provide novel insights that inform our comprehending of disease pathogenesis and provide potential targets for diagnosis and treatment.

Phylodynamics for cell biologists

Multicellular organisms are composed of cells connected by ancestry and descent from progenitor cells. The dynamics of cell birth, death, and inheritance within an organism give rise to the fundamental processes of development, differentiation, and cancer. Technical advances in molecular biology now allow us to study cellular composition, ancestry, and evolution at the resolution of individual cells within an organism or tissue. Here, we take a phylogenetic and phylodynamic approach to single-cell biology. We explain how "tree thinking" is important to the interpretation of the growing body of cell-level data and how ecological null models can benefit statistical hypothesis testing. Experimental progress in cell biology should be accompanied by theoretical developments if we are to exploit fully the dynamical information in single-cell data.

Correlations in Somatic Hypermutation Between Sites in IGHV Genes Can Be Explained by Interactions Between AID and/or Polη Hotspots

The somatic hypermutation (SHM) of Immunoglobulin (Ig) genes is a key process during antibody affinity maturation in B cells. The mutagenic enzyme activation induced deaminase (AID) is required for SHM and has a preference for WRC hotspots in DNA. Error-prone repair mechanisms acting downstream of AID introduce further mutations, including DNA polymerase eta (Polη), part of the non-canonical mismatch repair pathway (ncMMR), which preferentially generates mutations at WA hotspots. Previously proposed mechanistic models lead to a variety of predictions concerning interactions between hotspots, for example, how mutations in one hotspot will affect another hotspot. Using a large, high-quality, Ig repertoire sequencing dataset, we evaluated pairwise correlations between mutations site-by-site using an unbiased measure similar to mutual information which we termed “mutational association” (MA). Interactions are dominated by relatively strong correlations between nearby sites (short-range MAs), which can be almost entirely explained by interactions between overlapping hotspots for AID and/or Polη. We also found relatively weak dependencies between almost all sites throughout each gene (longer-range MAs), although these arise mostly as a statistical consequence of high pairwise mutation frequencies. The dominant short-range interactions are also highest within the most highly mutating IGHV sub-regions, such as the complementarity determining regions (CDRs), where there is a high hotspot density. Our results suggest that the hotspot preferences for AID and Polη have themselves evolved to allow for greater interactions between AID and/or Polη induced mutations.

Decrease of Clone Diversity in IgM Repertoires of HBV Chronically Infected Individuals With High Level of Viral Replication

High-throughput antibody sequencing allows in-depth insights into human antibody repertoires. To investigate the characteristics of antibody repertoires in patients with chronic HBV infection, we performed Illumina sequencing and IMGT/HighV-QUEST analysis of B lymphocytes from healthy adults and the HBV carriers with high or low level of viral replication. The comparative study revealed high levels of similarity between the IgM and IgG repertoires of the HBV carriers and the healthy adults, including the somatic mutations in V regions, the average CDR3 length, and the occurrence of junctional modifications. Nevertheless, the diversity of the unique clones decreased and some clusters of unique clones expanded in the IgM repertoire of chronic HBV carriers (CHB) compared with healthy adults (HH) and inactive HBV carriers (IHB). Such difference in clone diversity and expansion was not observed in the IgG repertoires of the three populations. More shared antibody clones were found between the IgM repertoires of IHB and HH than that found between CHB and HH (7079 clones vs. 2304 clones). Besides, the biased used IGHD genes were IGHD2-2 and IGHD3-3 in CHB library but were IGHD3-10 and IGHD3-22 in IHB and HH library. In contrast, for IgG repertories, the preferred used VDJ genes were similar in all the three populations. These results indicated that low level of serum HBV might not induce significant changes in BCR repertoires, and high level of HBV replication could have more impacts on IgM repertories than IgG repertoires. Taken together, our findings provide a better understanding of the antibody repertoires of HBV chronically infected individuals.

Infectious disease antibodies for biomedical applications: A mini review of immune antibody phage library repertoire

Antibody phage display is regarded as a critical tool for the development of monoclonal antibodies for infectious diseases. The different classes of antibody libraries are classified based on the source of repertoire used to generate the libraries. Immune antibody libraries are generated from disease infected host or immunization against an infectious agent. Antibodies derived from immune libraries are distinct from those derived from naïve libraries as the host's in vivo immune mechanisms shape the antibody repertoire to yield high affinity antibodies. As the immune system is constantly evolving in accordance to the health state of an individual, immune libraries can offer more than just infection-specific antibodies but also antibodies derived from the memory B-cells much like naïve libraries. The combinatorial nature of the gene cloning process would give rise to a combination of natural and un-natural antibody gene pairings in the immune library. These factors have a profound impact on the coverage of immune antibody libraries to target both disease-specific and non-disease specific antigens. This review looks at the diverse nature of antibody responses for immune library generation and discusses the extended potential of a disease-specified immune library in the context of phage display.

Beyond bulk single-chain sequencing: Getting at the whole receptor

Recent advancements in paired B-cell receptor sequencing technologies have accelerated the development of simpler, high-throughput pipelines for generating native antibody heavy and light chain pairs used to elucidate novel antibodies and provide insights into antibody response against pathogenic targets. These technologies involve single-cell isolation, using either single wells or emulsified droplets to maintain physical separation of individual cells, followed by sequencing. The development of novel single wells and emulsion-based workflows addresses key challenges by improving throughput of single-cell analyses, reducing method complexity, and integrating functional assays into existing workflows. Enabled by paired B-cell receptor sequencing, functional characterization of pathogen-specific antibodies reveals immunological insights beyond bulk sequencing.

Immune Literacy: Reading, Writing, and Editing Adaptive Immunity

Advances in reading, writing, and editing DNA are providing unprecedented insights into the complexity of immunological systems. This combination of systems and synthetic biology methods is enabling the quantitative and precise understanding of molecular recognition in adaptive immunity, thus providing a framework for reprogramming immune responses for translational medicine. In this review, we will highlight state-of-the-art methods such as immune repertoire sequencing, immunoinformatics, and immunogenomic engineering and their application toward adaptive immunity. We showcase novel and interdisciplinary approaches that have the promise of transforming the design and breadth of molecular and cellular immunotherapies.

AncesTree: An interactive immunoglobulin lineage tree visualizer

High-throughput sequencing of human immunoglobulin genes allows analysis of antibody repertoires and the reconstruction of clonal lineage evolution. The study of antibodies (Abs) affinity maturation is of specific interest to understand the generation of Abs with high affinity or broadly neutralizing activities. Moreover, phylogenic analysis enables the identification of the key somatic mutations required to achieve optimal antigen binding. The Immcantation framework provides a start-to-finish set of analytical methods for high-throughput adaptive immune receptor repertoire sequencing (AIRR-Seq; Rep-Seq) data. Furthermore, Immcantation's Change-O package has developed IgPhyML, an algorithm designed to build specifically immunoglobulin (Ig) phylogenic trees. Meanwhile Phylip, an algorithm that has been originally developed for applications in ecology and macroevolution, can also be used for the phylogenic reconstruction of antibodies maturation pathway. To complement Ig lineages made by IgPhyML or Dnaml (Phylip), we developed AncesTree, a graphic user interface (GUI) that aims to give researchers the opportunity to interactively explore antibodies clonal evolution. AncesTree displays interactive immunoglobulins phylogenic tree, Ig related mutations and sequence alignments using additional information coming from specialized antibody tools. The GUI is a Java standalone application allowing interaction with Ig tree that can run under Windows, Linux and Mac OS.

Computational Model Reveals a Stochastic Mechanism behind Germinal Center Clonal Bursts

Germinal centers (GCs) are specialized compartments within the secondary lymphoid organs where B cells proliferate, differentiate, and mutate their antibody genes in response to the presence of foreign antigens. Through the GC lifespan, interclonal competition between B cells leads to increased affinity of the B cell receptors for antigens accompanied by a loss of clonal diversity, although the mechanisms underlying clonal dynamics are not completely understood. We present here a multi-scale quantitative model of the GC reaction that integrates an intracellular component, accounting for the genetic events that shape B cell differentiation, and an extracellular stochastic component, which accounts for the random cellular interactions within the GC. In addition, B cell receptors are represented as sequences of nucleotides that mature and diversify through somatic hypermutations. We exploit extensive experimental characterizations of the GC dynamics to parameterize our model, and visualize affinity maturation by means of evolutionary phylogenetic trees. Our explicit modeling of B cell maturation enables us to characterise the evolutionary processes and competition at the heart of the GC dynamics, and explains the emergence of clonal dominance as a result of initially small stochastic advantages in the affinity to antigen. Interestingly, a subset of the GC undergoes massive expansion of higher-affinity B cell variants (clonal bursts), leading to a loss of clonal diversity at a significantly faster rate than in GCs that do not exhibit clonal dominance. Our work contributes towards an in silico vaccine design, and has implications for the better understanding of the mechanisms underlying autoimmune disease and GC-derived lymphomas.

Diversity in biology: definitions, quantification and models

Diversity indices are useful single-number metrics for characterizing a complex distribution of a set of attributes across a population of interest. The utility of these different metrics or sets of metrics depends on the context and application, and whether a predictive mechanistic model exists. In this topical review, we first summarize the relevant mathematical principles underlying heterogeneity in a large population, before outlining the various definitions of 'diversity' and providing examples of scientific topics in which its quantification plays an important role. We then review how diversity has been a ubiquitous concept across multiple fields, including ecology, immunology, cellular barcoding experiments, and socioeconomic studies. Since many of these applications involve sampling of populations, we also review how diversity in small samples is related to the diversity in the entire population. Features that arise in each of these applications are highlighted.

Transcriptomic Changes in Liver of Juvenile Cynoglossus semilaevis following Perfluorooctane Sulfonate Exposure

Perfluorooctane sulfonate (PFOS) is an increasingly important environmental pollutant that is pervasive in the environment. A number of studies have focused on the toxicological effects of PFOS on model fish species (zebrafish and medaka), but little is known about the impact of PFOS on commercially important marine fish. Thus, the present study examined transcriptome responses to PFOS exposure in the liver of juvenile Cynoglossus semilaevis, an important farmed flatfish in China. Then, in response to PFOS challenges, 1695 and 5244 genes were identified as significantly increased and depressed, respectively. Significant expression changes were observed in immune-related genes (cytokine-cytokine receptor interaction, T-helper [Th]17 cell differentiation, and the chemokine nuclear factor-kappa B and T-cell receptor signaling pathways), indicating that immunotoxicity is a key aspect of the effects of PFOS on C. semilaevis. Exposure to PFOS also altered the gene expression levels of hormones (inhibin, insulin, somatostatin, and glucagon), which could lead to severe metabolic and endocrine dysfunction. As expected from previous studies, several phase I and phase II detoxification enzymes were significantly up-regulated, which could facilitate the biotransformation and detoxification of PFOS in C. semilaevis. The present study provides new insights into the molecular toxicology of PFOS in a commercially important fish species. Environ Toxicol Chem 2020;39:556-564. © 2019 SETAC.

Inferring B cell specificity for vaccines using a Bayesian mixture model

Background

Vaccines have greatly reduced the burden of infectious disease, ranking in their impact on global health second only after clean water. Most vaccines confer protection by the production of antibodies with binding affinity for the antigen, which is the main effector function of B cells. This results in short term changes in the B cell receptor (BCR) repertoire when an immune response is launched, and long term changes when immunity is conferred. Analysis of antibodies in serum is usually used to evaluate vaccine response, however this is limited and therefore the investigation of the BCR repertoire provides far more detail for the analysis of vaccine response.

Results

Here, we introduce a novel Bayesian model to describe the observed distribution of BCR sequences and the pattern of sharing across time and between individuals, with the goal to identify vaccine-specific BCRs. We use data from two studies to assess the model and estimate that we can identify vaccine-specific BCRs with 69% sensitivity.

Conclusion

Our results demonstrate that statistical modelling can capture patterns associated with vaccine response and identify vaccine specific B cells in a range of different data sets. Additionally, the B cells we identify as vaccine specific show greater levels of sequence similarity than expected, suggesting that there are additional signals of vaccine response, not currently considered, which could improve the identification of vaccine specific B cells.

Broadly neutralizing human antibodies against dengue virus identified by single B cell transcriptomics

Eliciting broadly neutralizing antibodies (bNAbs) against the four dengue virus serotypes (DENV1-4) that are spreading into new territories is an important goal of vaccine design. To define bNAb targets, we characterized 28 antibodies belonging to expanded and hypermutated clonal families identified by transcriptomic analysis of single plasmablasts from DENV-infected individuals. Among these, we identified J9 and J8, two somatically related bNAbs that potently neutralized DENV1-4. Mutagenesis studies showed that the major recognition determinants of these bNAbs are in E protein domain I, distinct from the only known class of human bNAbs against DENV with a well-defined epitope. B cell repertoire analysis from acute-phase peripheral blood suggested that J9 and J8 followed divergent somatic hypermutation pathways, and that a limited number of mutations was sufficient for neutralizing activity. Our study suggests multiple B cell evolutionary pathways leading to DENV bNAbs targeting a new epitope that can be exploited for vaccine design.

Repertoire-wide phylogenetic models of B cell molecular evolution reveal evolutionary signatures of aging and vaccination

In order to produce effective antibodies, B cells undergo rapid somatic hypermutation (SHM) and selection for binding affinity to antigen via a process called affinity maturation. The similarities between this process and evolution by natural selection have led many groups to use phylogenetic methods to characterize the development of immunological memory, vaccination, and other processes that depend on affinity maturation. However, these applications are limited by the fact that most phylogenetic models are designed to be applied to individual lineages comprising genetically diverse sequences, while B cell repertoires often consist of hundreds to thousands of separate low-diversity lineages. Further, several features of affinity maturation violate important assumptions in standard phylogenetic models. Here, we introduce a hierarchical phylogenetic framework that integrates information from all lineages in a repertoire to more precisely estimate model parameters while simultaneously incorporating the unique features of SHM. We demonstrate the power of this repertoire-wide approach by characterizing previously undescribed phenomena in affinity maturation. First, we find evidence consistent with age-related changes in SHM hot-spot targeting. Second, we identify a consistent relationship between increased tree length and signs of increased negative selection, apparent in the repertoires of recently vaccinated subjects and those without any known recent infections or vaccinations. This suggests that B cell lineages shift toward negative selection over time as a general feature of affinity maturation. Our study provides a framework for undertaking repertoire-wide phylogenetic testing of SHM hypotheses and provides a means of characterizing dynamics of mutation and selection during affinity maturation.

Sequencing the peripheral blood B and T cell repertoire - Quantifying robustness and limitations

The adaptive immune response generates a large repertoire of T cells with T-cell receptors (TCR_alpha and TCR_beta) and B cells with immunoglobulins (Ig). The repertoire changes in response to antigen stimulation both through amplification of specific cells (clonal expansion) as well as somatic hypermutation of immunoglobulins. Alterations of the immune repertoire have been observed in response to acute disease, such as external pathogens, or chronic diseases, such as autoimmunity and cancer. Here we establish experimental and analytical protocols for quantifying the peripheral blood of healthy human individuals by profiling the immune repertoire for the Complementarity determining region 3 (CDR3) of the variable regions of TCR_beta (CDRβ3) and the IgG heavy chain (CDRH1, CDRH2, CDRH3). The results demonstrate that 40 ml of blood are sufficient to reliably capture the 10,000 most common TCR_beta and 1000 most common IgG and determine their relative frequency in the circulation. We conclude that by using an accessible sample size of human PBMC one is able to robustly monitor alterations in the immune repertoire.

Using High-Throughput Sequencing to Characterize the Development of the Antibody Repertoire During Infections: A Case Study of HIV-1

High throughput sequencing (HTS) approaches have only recently been applied to describing the antibody/B-cell repertoire in fine detail, but these data sets have already become critical to the design of vaccines and therapeutics, and monitoring of cancer immunotherapy. As a case study, we describe the potential and present limitations of HTS studies of the Ab repertoire during infection with HIV-1. Most of the present studies restrict their analyses to lineages of specific bnAbs. We discuss future initiatives to expand this type of analysis to more complete repertoires and to improve comparing and sharing of these Ab repertoire data across studies and institutions.

Structurally Mapping Antibody Repertoires

Every human possesses millions of distinct antibodies. It is now possible to analyze this diversity via next-generation sequencing of immunoglobulin genes (Ig-seq). This technique produces large volume sequence snapshots of B-cell receptors that are indicative of the antibody repertoire. In this paper, we enrich these large-scale sequence datasets with structural information. Enriching a sequence with its structural data allows better approximation of many vital features, such as its binding site and specificity. Here, we describe the structural annotation of antibodies pipeline that maps the outputs of large Ig-seq experiments to known antibody structures. We demonstrate the viability of our protocol on five separate Ig-seq datasets covering ca. 35 m unique amino acid sequences from ca. 600 individuals. Despite the great theoretical diversity of antibodies, we find that the majority of sequences coming from such studies can be reliably mapped to an existing structure.

Novel Concepts of Altered Immunoglobulin G Galactosylation in Autoimmune Diseases

The composition of the conserved N297 glycan in immunoglobulin G (IgG) has been shown to affect antibody effector functions via C1q of the complement system and Fc gamma receptors (FcγR) on immune cells. Changes in the general levels of IgG-glycoforms, such as lowered total IgG galactosylation observed in many autoimmune diseases have been associated with elevated disease severity. Agalactosyslated IgG has therefore been regarded and classified by many as pro-inflammatory. However, and somewhat counterintuitively, agalactosylation has been shown by several groups to decrease affinity for FcγRIII and decrease C1q binding and downstream activation, which seems at odds with this proposed pro-inflammatory nature. In this review, we discuss these circumstances where altered IgG galactosylation/glycosylation is found. We propose a novel model based on these observations and current biochemical evidence, where the levels of IgG galactosylation found in the total bulk IgG affect the threshold required to achieve immune activation by autoantibodies through either C1q or FcγR. Although this model needs experimental verification, it is supported by several clinical observations and reconciles apparent discrepancies in the literature, and suggests a general mechanism in IgG-mediated autoimmune diseases.

Analyzing Immunoglobulin Repertoires

Somatic assembly of T cell receptor and B cell receptor (BCR) genes produces a vast diversity of lymphocyte antigen recognition capacity. The advent of efficient high-throughput sequencing of lymphocyte antigen receptor genes has recently generated unprecedented opportunities for exploration of adaptive immune responses. With these opportunities have come significant challenges in understanding the analysis techniques that most accurately reflect underlying biological phenomena. In this regard, sample preparation and sequence analysis techniques, which have largely been borrowed and adapted from other fields, continue to evolve. Here, we review current methods and challenges of library preparation, sequencing and statistical analysis of lymphocyte receptor repertoire studies. We discuss the general steps in the process of immune repertoire generation including sample preparation, platforms available for sequencing, processing of sequencing data, measurable features of the immune repertoire, and the statistical tools that can be used for analysis and interpretation of the data. Because BCR analysis harbors additional complexities, such as immunoglobulin (Ig) (i.e., antibody) gene somatic hypermutation and class switch recombination, the emphasis of this review is on Ig/BCR sequence analysis.

Computational Strategies for Dissecting the High-Dimensional Complexity of Adaptive Immune Repertoires

The adaptive immune system recognizes antigens via an immense array of antigen-binding antibodies and T-cell receptors, the immune repertoire. The interrogation of immune repertoires is of high relevance for understanding the adaptive immune response in disease and infection (e.g., autoimmunity, cancer, HIV). Adaptive immune receptor repertoire sequencing (AIRR-seq) has driven the quantitative and molecular-level profiling of immune repertoires, thereby revealing the high-dimensional complexity of the immune receptor sequence landscape. Several methods for the computational and statistical analysis of large-scale AIRR-seq data have been developed to resolve immune repertoire complexity and to understand the dynamics of adaptive immunity. Here, we review the current research on (i) diversity, (ii) clustering and network, (iii) phylogenetic, and (iv) machine learning methods applied to dissect, quantify, and compare the architecture, evolution, and specificity of immune repertoires. We summarize outstanding questions in computational immunology and propose future directions for systems immunology toward coupling AIRR-seq with the computational discovery of immunotherapeutics, vaccines, and immunodiagnostics.

Regulation of B cell fate, survival, and function by mitochondria and autophagy

B cells are responsible for protective antibody production after differentiation into antibody-secreting cells during humoral immune responses. From early B cell development in the bone marrow, to their maturation in the periphery, activation in the germinal center, and differentiation into plasma cells or memory B cells, B cells display ever-changing functions and properties. Autophagy and mitochondria play important roles in B cell development, activation, and differentiation to accommodate the phenotypic and environmental changes encountered over the lifetime of the cell. Among their many functions, mitochondria and autophagy generate energy, mediate cell survival, and produce/eliminate reactive oxygen species that can serve as signal molecules to regulate differentiation. As B cells mature and differentiate into plasma or memory cells, both autophagic and mitochondrial functions undergo significant changes. In this review, we aim to provide an overview of the role of the autophagosome and mitochondria in regulating B cell fate, survival, and function. Moreover, we will discuss the interplay between these two highly metabolic organelles during B cell development, maturation, and differentiation.

Immune Repertoire after Immunization As Seen by Next-Generation Sequencing and Proteomics

The immune system produces a diverse repertoire of immunoglobulins in response to foreign antigens. During B-cell development, VDJ recombination and somatic mutations generate diversity, whereas selection processes remove it. Using both proteomic and NGS approaches, we characterized the immune repertoires in groups of rats after immunization with purified antigens. Proteomics and NGS data on the repertoire are in qualitative agreement, but did show quantitative differences that may relate to differences between the biological niches that were sampled for these approaches. Both methods contributed complementary information in the characterization of the immune repertoire. It was found that the immune repertoires resulting from each antigen had many similarities that allowed samples to cluster together, and that mutated immunoglobulin peptides were shared among animals with a response to the same antigen significantly more than for different antigens. However, the number of shared sequences decreased in a log-linear fashion relative to the number of animals that share them, which may affect future applications. A phylogenetic analysis on the NGS reads showed that reads from different individuals immunized with the same antigen populated distinct branches of the phylogram, an indication that the repertoire had converged. Also, similar mutation patterns were found in branches of the phylogenetic tree that were associated with antigen-specific immunoglobulins through proteomics data. Thus, data from different analysis methods and different experimental platforms show that the immunoglobulin repertoires of immunized animals have overlapping and converging features. With additional research, this may enable interesting applications in biotechnology and clinical diagnostics.

Methods and Applications of CRISPR-Mediated Base Editing in Eukaryotic Genomes

The past several years have seen an explosion in development of applications for the CRISPR-Cas9 system, from efficient genome editing, to high-throughput screening, to recruitment of a range of DNA and chromatin-modifying enzymes. While homology-directed repair (HDR) coupled with Cas9 nuclease cleavage has been used with great success to repair and re-write genomes, recently developed base-editing systems present a useful orthogonal strategy to engineer nucleotide substitutions. Base editing relies on recruitment of cytidine deaminases to introduce changes (rather than double-stranded breaks and donor templates) and offers potential improvements in efficiency while limiting damage and simplifying the delivery of editing machinery. At the same time, these systems enable novel mutagenesis strategies to introduce sequence diversity for engineering and discovery. Here, we review the different base-editing platforms, including their deaminase recruitment strategies and editing outcomes, and compare them to other CRISPR genome-editing technologies. Additionally, we discuss how these systems have been applied in therapeutic, engineering, and research settings. Lastly, we explore future directions of this emerging technology.