Network Signatures of IgG Immune Repertoires in Hepatitis B Associated Chronic Infection and Vaccination Responses

Yclon: Ultrafast clustering of B cell clones from high-throughput immunoglobulin repertoire sequencing data

MOTIVATION

The next-generation sequencing technologies have transformed our understanding of immunoglobulin (Ig) profiles in various immune states. Clonotyping, which groups Ig sequences into B cell clones, is crucial in investigating the diversity of repertoires and changes in antigen exposure. Despite its importance, there is no widely accepted method for clonotyping, and existing methods are computationally intensive for large sequencing datasets.

RESULTS

To address this challenge, we introduce YClon, a fast and efficient approach for clonotyping Ig repertoire data. YClon uses a hierarchical clustering approach, similar to other methods, to group Ig sequences into B cell clones in a highly sensitive and specific manner. Notably, our approach outperforms other methods by being more than 30 to 5000 times faster in processing the repertoires analyzed. Astonishingly, YClon can effortlessly handle up to 2 million Ig sequences on a standard laptop computer. This enables in-depth analysis of large and numerous antibody repertoires.

AVAILABILITY AND IMPLEMENTATION

YClon was implemented in Python3 and is freely available on GitHub.

Exploring the impact of clonal definition on B-cell diversity: implications for the analysis of immune repertoires

The adaptive immune system has the extraordinary ability to produce a broad range of immunoglobulins that can bind a wide variety of antigens. During adaptive immune responses, activated B cells duplicate and undergo somatic hypermutation in their B-cell receptor (BCR) genes, resulting in clonal families of diversified B cells that can be related back to a common ancestor. Advances in high-throughput sequencing technologies have enabled the high-throughput characterization of B-cell repertoires, however, the accurate identification of clonally related BCR sequences remains a major challenge. In this study, we compare three different clone identification methods on both simulated and experimental data, and investigate their impact on the characterization of B-cell diversity. We observe that different methods lead to different clonal definitions, which affects the quantification of clonal diversity in repertoire data. Our analyses show that direct comparisons between clonal clusterings and clonal diversity of different repertoires should be avoided if different clone identification methods were used to define the clones. Despite this variability, the diversity indices inferred from the repertoires' clonal characterization across samples show similar patterns of variation regardless of the clonal identification method used. We find the Shannon entropy to be the most robust in terms of the variability of diversity rank across samples. Our analysis also suggests that the traditional germline gene alignment-based method for clonal identification remains the most accurate when the complete information about the sequence is known, but that alignment-free methods may be preferred for shorter sequencing read lengths. We make our implementation freely available as a Python library cdiversity.

Lower Somatic Mutation Levels in the λ Light-Chain Repertoires with Chronic HBV Infection

To investigate the characteristics of the immunoglobulin light-chain repertoires with chronic HBV infection, the high-throughput sequencing and IMGT/HighV-QUEST were adapted to analyze the κ (IgK) and λ (IgL) light-chain repertoires from the inactive HBV carriers (IHB) and the healthy adults (HH). The comparative analysis revealed high similarity between the κ light-chain repertoires of the HBV carriers and the healthy adults. Nevertheless, the proportion of IGLV genes with ≥90% identity as the germline genes was higher in the IgL light-chain repertoire of the IHB library compared with that of HH library (74.6% vs. 69.1%). Besides, the frequency of amino acid mutations in the CDR1 regions was significantly lower in the IgL light-chain repertoire of the IHB library than that of the HH library (65.52% vs. 56.0%). These results suggested the lower somatic mutation level in the IgL repertoire of IHB library, which might indicate the biased selection of IGLV genes in the IgL repertoire with chronic HBV infection. These findings might lead to a better understanding of the characteristics of the light-chain repertoires of HBV chronically infected individuals.

Human immune repertoire in hepatitis B virus infection

Hepatitis B virus (HBV) infection is a public health threat that affects 257 million people worldwide and can progress to liver cirrhosis, liver failure, and hepatocellular carcinoma. The HBV antigen- induced adaptive immune response plays an important role in HBV clearance. Immune repertoire sequencing (IRS) has been used to investigate the molecular mechanisms behind the immune system, find novel ways to treat HBV infection, and evaluate the genetic responses and immune characteristics of individuals infected by HBV or immunized by HBV vaccine. This review summarizes the human immune repertoire analysis methodology, and the application of the IRS in the prediction of HBV infection progression, treatment, and vaccination.

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.

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.

Large-scale network analysis reveals the sequence space architecture of antibody repertoires

The architecture of mouse and human antibody repertoires is defined by the sequence similarity networks of the clones that compose them. The major principles that define the architecture of antibody repertoires have remained largely unknown. Here, we establish a high-performance computing platform to construct large-scale networks from comprehensive human and murine antibody repertoire sequencing datasets (>100,000 unique sequences). Leveraging a network-based statistical framework, we identify three fundamental principles of antibody repertoire architecture: reproducibility, robustness and redundancy. Antibody repertoire networks are highly reproducible across individuals despite high antibody sequence dissimilarity. The architecture of antibody repertoires is robust to the removal of up to 50-90% of randomly selected clones, but fragile to the removal of public clones shared among individuals. Finally, repertoire architecture is intrinsically redundant. Our analysis provides guidelines for the large-scale network analysis of immune repertoires and may be used in the future to define disease-associated and synthetic repertoires.

Next-generation sequencing reveals new insights about gene usage and CDR-H3 composition in the horse antibody repertoire

Horse serum antibodies have been used for greater than a century for the treatment and prophylaxis of infectious diseases and envenomations. Little is known, however, about the immunogenetic diversity that produces horse serum antibodies. Here, we employed next-generation sequencing for a first-in-kind comprehensive analysis of the equine B-cell repertoire. Nearly 45,000 and 30,000 clonotypes were obtained for the heavy-chain (IGH) and lambda light-chain (IGL) loci, respectively. We observed skewed use of the common subgroups IGHV2 (92.49%) and IGLV8 (82.50%), consistent with previous reports, but also novel use of the rare genes IGHV6S1 and IGLV4S2. CDR-H3 amino acid composition revealed different amino acid patterns at positions 106 and 116 compared to human, rabbit, and mouse, suggesting that an extended conformation predominates among horse CDR-H3 loops. Our analysis provides new insights regarding the mechanisms employed to generate antibody diversity in the horse, and could be applicable to the optimized design of synthetic antibodies intended for future therapeutic use.

VDJServer: A Cloud-Based Analysis Portal and Data Commons for Immune Repertoire Sequences and Rearrangements

Background

Recent technological advances in immune repertoire sequencing have created tremendous potential for advancing our understanding of adaptive immune response dynamics in various states of health and disease. Immune repertoire sequencing produces large, highly complex data sets, however, which require specialized methods and software tools for their effective analysis and interpretation.

Results

VDJServer is a cloud-based analysis portal for immune repertoire sequence data that provide access to a suite of tools for a complete analysis workflow, including modules for preprocessing and quality control of sequence reads, V(D)J gene segment assignment, repertoire characterization, and repertoire comparison. VDJServer also provides sophisticated visualizations for exploratory analysis. It is accessible through a standard web browser via a graphical user interface designed for use by immunologists, clinicians, and bioinformatics researchers. VDJServer provides a data commons for public sharing of repertoire sequencing data, as well as private sharing of data between users. We describe the main functionality and architecture of VDJServer and demonstrate its capabilities with use cases from cancer immunology and autoimmunity.

Conclusion

VDJServer provides a complete analysis suite for human and mouse T-cell and B-cell receptor repertoire sequencing data. The combination of its user-friendly interface and high-performance computing allows large immune repertoire sequencing projects to be analyzed with no programming or software installation required. VDJServer is a web-accessible cloud platform that provides access through a graphical user interface to a data management infrastructure, a collection of analysis tools covering all steps in an analysis, and an infrastructure for sharing data along with workflows, results, and computational provenance. VDJServer is a free, publicly available, and open-source licensed resource.

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.

Shared IgG Infection Signatures vs. Hemorrhage-Restricted IgA Clusters in Human Dengue: A Phenotype of Differential Class-Switch via TGFβ1

Phenotypic manifestations of infectious diseases are closely related to individual immune responses. Methods to extract information from patients' own immune reactions would be of great use for both diagnosis and treatment. Dengue fever is one of the diseases that clinical aggravations could occur paradoxically after humoral immunity appears. This property makes dengue fever an excellent disease model to explore. A principal component analyses (PCAs)-based framework derived from a prior vaccination study was developed. The framework was verified by successful demonstrations of known IgG signatures from a Mexico Dengue data set. Afterward the pipeline was tested upon de novo IgG and IgA libraries of Dengue patients from southern Taiwan. We discovered four infection signatures within IgG repertoires, two of which were identical to previous reports. However, it was IgA but not IgG that could differentiate hemorrhagic from non-hemorrhagic patients. IgA repertoires were found more diversified among bleeders, from whom seven signature clusters were characterized. The expressions of transforming growth factor beta 1 (TGFβ1) and accordingly mediated class-switch activity of IgA were distinct only among the PCA-segregated bleeding group. In sum, intercontinental sharing of IgG signatures in dengue fever was demonstrated via a unified working flow. Differential regulation of IgA class-switch with associated diversity expansion plus existences of hemorrhage-restricted clusters were shown. The ability of the framework to find common IgG signatures would implicate applications to infections even from unknown pathogens. The clusters within IgA repertoires could offer perspectives to other IgA-related bleeding disorders such as Henoch-Schönlein purpura or IgA nephropathy. Substantiated grounds for IgA-specific effector function via TGFβ1-mediated class-switch would be a new factor to consider for infectious diseases.

Characteristics Peripheral Blood IgG and IgM Heavy Chain Complementarity Determining Region 3 Repertoire before and after Immunization with Recombinant HBV Vaccine

Immunization with recombinant HBV vaccine induces specific immune responses in human causing B lymphocytes to produce protective HBsAb, and to form memory B lymphocytes, thereby facilitating HBV immunity in the body. However, B lymphocytes heterogeneity and characteristics are not fully elucidated. In this study, we conducted high-throughput sequencing of BCR heavy chain CDR3 repertoires in 3 healthy volunteers before and after the third immunization with recombinant HBV vaccine. We used Roche 454 Genome Sequencer FLX system to perform a comparative analysis of IgM and IgG H chain CDR3 repertoires. First, we found that the diversity of IgG H chain CDR3 repertoires was 1/6 of IgM on average. Moreover, after the third immunization with HBV vaccine, the diversity of IgG H chain CDR3 repertoires was 1/26 of IgM on average. Second, we detected relatively high levels of HBsAbs in all the healthy volunteers after immunization with HBV vaccine. The volunteers shared a small number of CDR3 sequences before and after immunization, and among each other. However, we did not find completely identical BCR H chain CDR3 amino acid sequences in these volunteers. Lastly, after immunization with recombinant HBV vaccine, the volunteers showed high frequency of IgG H chain CDR3 amino acid sequences mostly resulting from rearrangements of IGHV, IGHD and IGHJ, suggesting that the mechanism of high frequency CDR3 generation might be associated with the maturation of IgG affinity (somatic hypermutation) during the recombinant HBV vaccine-induced B lymphocyte responses. This study identified the characteristics and changes of BCR CDR3 repertoires before and after immunization with HBV vaccine, and evaluated the performance of the sequencing technology for this application. Our findings provide a basis for further research in B lymphocyte generated HBsAb heterogeneity and monitoring the maintenance of memory B lymphocytes.