High-Throughput Single-Cell Analysis of B Cell Receptor Usage among Autoantigen-Specific Plasma Cells in Celiac Disease

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

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

New Insights on Genes, Gluten, and Immunopathogenesis of Celiac Disease

Celiac disease (CeD) is a gluten-induced enteropathy that develops in genetically susceptible individuals upon consumption of cereal gluten proteins. It is a unique and complex immune disorder to study as the driving antigen is known and the tissue targeted by the immune reaction can be interrogated. This review integrates findings gained from genetic, biochemical, and immunologic studies, which together have revealed mechanisms of gluten peptide modification and HLA binding, thereby enabling a maladapted anti-gluten immune response. Observations in human samples combined with experimental mouse models have revealed that the gluten-induced immune response involves CD4+ T cells, cytotoxic CD8+ T cells, and B cells; their cross-talks are critical for the tissue-damaging response. The emergence of high-throughput technologies is increasing our understanding of the phenotype, location, and presumably function of the gluten-specific cells, which are all required to identify novel therapeutic targets and strategies for CeD.

Selective activation of naïve B cells with unique epitope specificity shapes autoantibody formation in celiac disease

Many antibody responses induced by infection, vaccination or autoimmunity show signs of convergence across individuals with epitope-dependent selection of particular variable region gene segments and complementarity determining region 3 properties. However, not much is known about the relationship between antigen-specific effector cells and antigen-specific precursors present in the naïve B-cell repertoire. Here, we sought to address this relationship in the context of celiac disease, where there is a stereotyped autoantibody response against the enzyme transglutaminase 2 (TG2). By generating TG2-specific monoclonal antibodies from both duodenal plasma cells and circulating naïve B cells, we demonstrate a discord between the naïve TG2-specific repertoire and the cells that are selected for autoantibody production. Hence, the naïve repertoire does not fully reflect the epitope preference and gene usage observed for memory B cells and plasma cells. Instead, distinct naïve B cells that target particular TG2 epitopes appear to be selectively activated at the expense of TG2-binding B cells targeting other epitopes.

Generation of circulating autoreactive pre-plasma cells fueled by naive B cells in celiac disease

Autoantibodies against the enzyme transglutaminase 2 (TG2) are characteristic of celiac disease (CeD), and TG2-specific immunoglobulin (Ig) A plasma cells are abundant in gut biopsies of patients. Here, we describe the corresponding population of autoreactive B cells in blood. Circulating TG2-specific IgA cells are present in untreated patients on a gluten-containing diet but not in controls. They are clonally related to TG2-specific small intestinal plasma cells, and they express gut-homing molecules, indicating that they are plasma cell precursors. Unlike other IgA-switched cells, the TG2-specific cells are negative for CD27, placing them in the double-negative (IgD-CD27-) category. They have a plasmablast or activated memory B cell phenotype, and they harbor fewer variable region mutations than other IgA cells. Based on their similarity to naive B cells, we propose that autoreactive IgA cells in CeD are generated mainly through chronic recruitment of naive B cells via an extrafollicular response involving gluten-specific CD4+ T cells.

Coeliac disease: the paradox of diagnosing a food hypersensitivity disorder with autoantibodies

Serum antibodies to the autoantigen transglutaminase 2 (TG2) are increasingly harnessed to diagnose coeliac disease. Diagnostic guidelines for children give recommendation for a no-biopsy-based diagnosis through detection of high amounts of IgA anti-TG2 antibodies in serum with confirmation of positivity in a separate blood sample by characteristic autoantibody-staining of tissue. While measurement of IgA anti-TG2 also is important in the diagnostic workup of adults, the adult guidelines still mandate examination of gut biopsies. This requirement might well change in the future, as might the necessity for confirming autoantibody positivity by tissue staining. The key role of autoantibody serology for diagnosis of coeliac disease is paradoxical. Coeliac disease was considered, and still can be considered, a food intolerance disorder where autoantibodies at face value are out of place. The immunological mechanisms underlying the formation of autoantibodies in response to gluten exposure have been dissected. This review presents the current insights demonstrating that the autoantibodies in coeliac disease are intimately integrated in the maladapted immune response to gluten.

Separate Gut Plasma Cell Populations Produce Auto-Antibodies against Transglutaminase 2 and Transglutaminase 3 in Dermatitis Herpetiformis

Dermatitis herpetiformis (DH) is an inflammatory skin disorder often considered as an extra intestinal manifestation of celiac disease (CeD). Hallmarks of CeD and DH are auto-antibodies to transglutaminase 2 (TG2) and transglutaminase 3 (TG3), respectively. DH patients have auto-antibodies reactive with both transglutaminase enzymes. Here it is reported that in DH both gut plasma cells and serum auto-antibodies are specific for either TG2 or TG3 with no TG2-TG3 cross reactivity. By generating monoclonal antibodies from TG3-specific duodenal plasma cells of DH patients, three conformational epitope groups are defined. Both TG2-specific and TG3-specific gut plasma cells have few immunoglobulin (Ig) mutations, and the two transglutaminase-reactive populations show distinct selection of certain heavy and light chain V-genes. Mass spectrometry analysis of TG3-specific serum IgA corroborates preferential usage of IGHV2-5 in combination with IGKV4-1. Collectively, these results demonstrate parallel induction of anti-TG2 and anti-TG3 auto-antibody responses involving separate B-cell populations in DH patients.

Genetic variation in the immunoglobulin heavy chain locus shapes the human antibody repertoire

Variation in the antibody response has been linked to differential outcomes in disease, and suboptimal vaccine and therapeutic responsiveness, the determinants of which have not been fully elucidated. Countering models that presume antibodies are generated largely by stochastic processes, we demonstrate that polymorphisms within the immunoglobulin heavy chain locus (IGH) impact the naive and antigen-experienced antibody repertoire, indicating that genetics predisposes individuals to mount qualitatively and quantitatively different antibody responses. We pair recently developed long-read genomic sequencing methods with antibody repertoire profiling to comprehensively resolve IGH genetic variation, including novel structural variants, single nucleotide variants, and genes and alleles. We show that IGH germline variants determine the presence and frequency of antibody genes in the expressed repertoire, including those enriched in functional elements linked to V(D)J recombination, and overlapping disease-associated variants. These results illuminate the power of leveraging IGH genetics to better understand the regulation, function, and dynamics of the antibody response in disease.

Untangling immunoglobulin genotype-function associations

Immunoglobulin (IG) genes, encoding B cell receptors (BCRs), are fundamental components of the mammalian immune system, which evolved to recognize the diverse antigenic universe present in nature. To handle these myriad inputs, BCRs are generated through combinatorial recombination of a set of highly polymorphic germline genes, resulting in a vast repertoire of antigen receptors that initiate responses to pathogens and regulate commensals. Following antigen recognition and B cell activation, memory B cells and plasma cells form, allowing for the development of anamnestic antibody (Ab) responses. How inherited variation in IG genes impacts host traits, disease susceptibility, and Ab recall responses is a topic of great interest. Here, we consider approaches to translate emerging knowledge about IG genetic diversity and expressed repertoires to inform our understanding of Ab function in health and disease etiology. As our understanding of IG genetics grows, so will our need for tools to decipher preferences for IG gene or allele usage in different contexts, to better understand antibody responses at the population level.

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

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

Tango of B cells with T cells in the making of secretory antibodies to gut bacteria

Polymeric IgA and IgM are transported across the epithelial barrier from plasma cells in the lamina propria to exert a function in the gut lumen as secretory antibodies. Many secretory antibodies are reactive with the gut bacteria, and mounting evidence suggests that these antibodies are important for the host to control gut bacterial communities. However, we have incomplete knowledge of how bacteria-reactive secretory antibodies are formed. Antibodies from gut plasma cells often show bacterial cross-species reactivity, putting the degree of specificity behind anti-bacterial antibody responses into question. Such cross-species reactive antibodies frequently recognize non-genome-encoded membrane glycan structures. On the other hand, the T cell epitopes are peptides encoded in the bacterial genomes, thereby allowing a higher degree of predictable specificity on the T cell side of anti-bacterial immune responses. In this Perspective, we argue that the production of bacteria-reactive secretory antibodies is mainly controlled by the antigen specificity of T cells, which provide help to B cells. To be able to harness this system (for instance, for manipulation with vaccines), we need to obtain insight into the bacterial epitopes recognized by T cells in addition to characterizing the reactivity of the antibodies.

Functional antibodies exhibit light chain coherence

The vertebrate adaptive immune system modifies the genome of individual B cells to encode antibodies that bind particular antigens1. In most mammals, antibodies are composed of heavy and light chains that are generated sequentially by recombination of V, D (for heavy chains), J and C gene segments. Each chain contains three complementarity-determining regions (CDR1-CDR3), which contribute to antigen specificity. Certain heavy and light chains are preferred for particular antigens2-22. Here we consider pairs of B cells that share the same heavy chain V gene and CDRH3 amino acid sequence and were isolated from different donors, also known as public clonotypes23,24. We show that for naive antibodies (those not yet adapted to antigens), the probability that they use the same light chain V gene is around 10%, whereas for memory (functional) antibodies, it is around 80%, even if only one cell per clonotype is used. This property of functional antibodies is a phenomenon that we call light chain coherence. We also observe this phenomenon when similar heavy chains recur within a donor. Thus, although naive antibodies seem to recur by chance, the recurrence of functional antibodies reveals surprising constraint and determinism in the processes of V(D)J recombination and immune selection. For most functional antibodies, the heavy chain determines the light chain.

Pathophysiology and immunogenetics of celiac disease

Celiac disease (CD) is a chronic inflammatory enteropathy caused by gluten (protein from wheat, rye and, barley) in genetically predisposed individuals carrying the HLA-DQ2/HLA-DQ8 genotype. This pathology has a multifactorial etiology in which HLA genes, the microbiome, gluten and, other environmental factors are involved in the development of the disease. Its pathogenesis involves both innate and adaptive immunity as well as upregulation of IL-15. The objective of this review is to examine the results of current studies on genetic and environmental variables to better understand the pathogenesis of this enteropathy. The complex etiology of celiac disease makes our understanding of the pathogenesis of the disease incomplete, and a better knowledge of the many genetic and environmental components would help us better understand the pathophysiology of celiac disease.

Single-cell approaches to dissect adaptive immune responses involved in autoimmunity: the case of celiac disease

Single-cell analysis is a powerful technology that has found widespread use in recent years. For diseases with involvement of adaptive immunity, single-cell analysis of antigen-specific T cells and B cells is particularly informative. In autoimmune diseases, the adaptive immune system is obviously at play, yet the ability to identify the culprit T and B cells recognizing disease-relevant antigen can be difficult. Celiac disease, a widespread disorder with autoimmune components, is unique in that disease-relevant antigens for both T cells and B cells are well defined. Furthermore, the celiac disease gut lesion is readily accessible allowing for sampling of tissue-resident cells. Thus, disease-relevant T cells and B cells from the gut and blood can be studied at the level of single cells. Here we review single-cell studies providing information on such adaptive immune cells and outline some future perspectives in the area of single-cell analysis in autoimmune diseases.

Individualized VDJ recombination predisposes the available Ig sequence space

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

Signatures of B Cell Receptor Repertoire Following Pneumocystis Infection

B cells play vital roles in host defense against Pneumocystis infection. However, the features of the B cell receptor (BCR) repertoire in disease progression remain unclear. Here, we integrated single-cell RNA sequencing and single-cell BCR sequencing of immune cells from mouse lungs in an uninfected state and 1–4 weeks post-infection in order to illustrate the dynamic nature of B cell responses during Pneumocystis infection. We identified continuously increased plasma cells and an elevated ratio of (IgA + IgG) to (IgD + IgM) after infection. Moreover, Pneumocystis infection was associated with an increasing naïve B subset characterized by elevated expression of the transcription factor ATF3. The proportion of clonal expanded cells progressively increased, while BCR diversity decreased. Plasma cells exhibited higher levels of somatic hypermutation than naïve B cells. Biased usage of V(D)J genes was observed, and the usage frequency of IGHV9-3 rose. Overall, these results present a detailed atlas of B cell transcriptional changes and BCR repertoire features in the context of Pneumocystis infection, which provides valuable information for finding diagnostic biomarkers and developing potential immunotherapeutic targets.

Machine Learning Analysis of Naïve B-Cell Receptor Repertoires Stratifies Celiac Disease Patients and Controls

Celiac disease (CeD) is a common autoimmune disorder caused by an abnormal immune response to dietary gluten proteins. The disease has high heritability. HLA is the major susceptibility factor, and the HLA effect is mediated via presentation of deamidated gluten peptides by disease-associated HLA-DQ variants to CD4+ T cells. In addition to gluten-specific CD4+ T cells the patients have antibodies to transglutaminase 2 (autoantigen) and deamidated gluten peptides. These disease-specific antibodies recognize defined epitopes and they display common usage of specific heavy and light chains across patients. Interactions between T cells and B cells are likely central in the pathogenesis, but how the repertoires of naïve T and B cells relate to the pathogenic effector cells is unexplored. To this end, we applied machine learning classification models to naïve B cell receptor (BCR) repertoires from CeD patients and healthy controls. Strikingly, we obtained a promising classification performance with an F1 score of 85%. Clusters of heavy and light chain sequences were inferred and used as features for the model, and signatures associated with the disease were then characterized. These signatures included amino acid (AA) 3-mers with distinct bio-physiochemical characteristics and enriched V and J genes. We found that CeD-associated clusters can be identified and that common motifs can be characterized from naïve BCR repertoires. The results may indicate a genetic influence by BCR encoding genes in CeD. Analysis of naïve BCRs as presented here may become an important part of assessing the risk of individuals to develop CeD. Our model demonstrates the potential of using BCR repertoires and in particular, naïve BCR repertoires, as disease susceptibility markers.

Convergent antibody evolution and clonotype expansion following influenza virus vaccination

Recent advances in high-throughput single cell sequencing have opened up new avenues into the investigation of B cell receptor (BCR) repertoires. In this study, PBMCs were collected from 17 human participants vaccinated with the split-inactivated influenza virus vaccine during the 2016-2017 influenza season. A combination of Immune Repertoire Capture (IRCTM) technology and IgG sequencing was performed on ~7,800 plasmablast (PB) cells and preferential IgG heavy-light chain pairings were investigated. In some participants, a single expanded clonotype accounted for ~22% of their PB BCR repertoire. Approximately 60% (10/17) of participants experienced convergent evolution, possessing public PBs that were elicited independently in multiple participants. Binding profiles of one private and three public PBs confirmed they were all subtype-specific, cross-reactive hemagglutinin (HA) head-directed antibodies. Collectively, this high-resolution antibody repertoire analysis demonstrated the impact evolution can have on BCRs in response to influenza virus vaccination, which can guide future universal influenza prophylactic approaches.

Longevity, clonal relationship, and transcriptional program of celiac disease-specific plasma cells

Disease-specific plasma cells (PCs) reactive with transglutaminase 2 (TG2) or deamidated gluten peptides (DGPs) are abundant in celiac disease (CeD) gut lesions. Their contribution toward CeD pathogenesis is unclear. We assessed expression of markers associated with PC longevity in 15 untreated and 26 treated CeD patients in addition to 13 non-CeD controls and performed RNA sequencing with clonal inference and transcriptomic analysis of 3,251 single PCs. We observed antigen-dependent V-gene selection and stereotypic antibodies. Generation of recombinant DGP-specific antibodies revealed a key role of a heavy chain residue that displays polymorphism, suggesting that immunoglobulin gene polymorphisms may influence CeD-specific antibody responses. We identified transcriptional differences between CeD-specific and non-disease-specific PCs and between short-lived and long-lived PCs. The short-lived CD19+CD45+ phenotype dominated in untreated and short-term-treated CeD, in particular among disease-specific PCs but also in the general PC population. Thus, the disease lesion of untreated CeD is characterized by massive accumulation of short-lived PCs that are not only directed against disease-specific antigens.

A computational method for immune repertoire mining that identifies novel binders from different clonotypes, demonstrated by identifying anti-pertussis toxoid antibodies

Due to their shared genetic history, antibodies from the same clonotype often bind to the same epitope. This knowledge is used in immune repertoire mining, where known binders are used to search bulk sequencing repertoires to identify new binders. However, current computational methods cannot identify epitope convergence between antibodies from different clonotypes, limiting the sequence diversity of antigen-specific antibodies that can be identified. We describe how the antibody binding site, the paratope, can be used to cluster antibodies with common antigen reactivity from different clonotypes. Our method, paratyping, uses the predicted paratope to identify these novel cross clonotype matches. We experimentally validated our predictions on a pertussis toxoid dataset. Our results show that even the simplest abstraction of the antibody binding site, using only the length of the loops involved and predicted binding residues, is sufficient to group antigen-specific antibodies and provide additional information to conventional clonotype analysis.

Abbreviations: BCR: B-cell receptor; CDR: complementarity-determining region; PTx: pertussis toxoid

Immunoglobulin germline gene variation and its impact on human disease

Immunoglobulins (Ig) play an important role in the immune system both when expressed as antigen receptors on the cell surface of B cells and as antibodies secreted into extracellular fluids. The advent of high-throughput sequencing methods has enabled the investigation of human Ig repertoires at unprecedented depth. This has led to the discovery of many previously unreported germline Ig alleles. Moreover, it is becoming clear that convergent and stereotypic antibody responses are common where different individuals recognise defined antigenic epitopes with the use of the same Ig V genes. Thus, germline V gene variation is increasingly being linked to the differential capacity of generating an effective immune response, which might lead to varying disease susceptibility. Here, we review recent evidence of how germline variation in Ig genes impacts the Ig repertoire and its subsequent effects on the adaptive immune response in vaccination, infection, and autoimmunity.

The Usage of Human IGHJ Genes Follows a Particular Non-random Selection: The Recombination Signal Sequence May Affect the Usage of Human IGHJ Genes

The formation of the B cell receptor (BCR) heavy chain variable region is derived from the germline V(D)J gene rearrangement according to the “12/23” rule and the “beyond 12/23” rule. The usage frequency of each V(D)J gene in the peripheral BCR repertoires is related to the initial recombination, self-tolerance selection, and the clonal proliferative response. However, their specific differences and possible mechanisms are still unknown. We analyzed in-frame and out-of-frame BCR-H repertoires from human samples with normal physiological and various pathological conditions by high-throughput sequencing. Our results showed that IGHJ gene frequency follows a similar pattern which is previously known, where IGHJ4 is used at high frequency (>40%), IGHJ6/IGHJ3/IGHJ5 is used at medium frequencies (10∼20%), and IGH2/IGHJ1 is used at low frequency (<4%) under whether normal physiological or various pathological conditions. However, our analysis of the recombination signal sequences suggested that the conserved non-amer and heptamer and certain 23 bp spacer length may affect the initial IGHD-IGHJ recombination, which results in different frequencies of IGHJ genes among the initial BCR-H repertoire. Based on this “initial repertoire,” we recommend that re-evaluation and further investigation are needed when analyzing the significance and mechanism of IGHJ gene frequency in self-tolerance selection and the clonal proliferative response.

Emerging Strategies for Therapeutic Antibody Discovery from Human B Cells

Monoclonal antibodies from human sources are being increasingly recognized as valuable options in many therapeutic areas. These antibodies can show exquisite specificity and high potency while maintaining a desirable safety profile, having been matured and tolerized within human patients. However, the discovery of these antibodies presents important challenges, since the B cells encoding therapeutic antibodies can be rare in a typical blood draw and are short-lived ex vivo. Furthermore, the unique pairing of V_H and V_L domains in each B cell contributes to specificity and function; therefore, maintaining antibody chain pairing presents a throughput limitation. This work will review the various approaches aimed at addressing these challenges with an eye to next-generation methods for high-throughput discovery from the human B-cell repertoire.

B cell tolerance and antibody production to the celiac disease autoantigen transglutaminase 2

Autoantibodies to transglutaminase 2 (TG2) are hallmarks of celiac disease. To address B cell tolerance and autoantibody formation to TG2, we generated immunoglobulin knock-in (Ig KI) mice that express a prototypical celiac patient-derived anti-TG2 B cell receptor equally reactive to human and mouse TG2. We studied B cell development in the presence/absence of autoantigen by crossing the Ig KI mice to Tgm2-/- mice. Autoreactive B cells in Tgm2+/+ mice were indistinguishable from their naive counterparts in Tgm2-/- mice with no signs of clonal deletion, receptor editing, or B cell anergy. The autoreactive B cells appeared ignorant to their antigen, and they produced autoantibodies when provided T cell help. The findings lend credence to a model of celiac disease where gluten-reactive T cells provide help to autoreactive TG2-specific B cells by involvement of gluten-TG2 complexes, and they outline a general mechanism of autoimmunity with autoantibodies being produced by ignorant B cells on provision of T cell help.

Polymorphisms in human immunoglobulin heavy chain variable genes and their upstream regions

Germline variations in immunoglobulin genes influence the repertoire of B cell receptors and antibodies, and such polymorphisms may impact disease susceptibility. However, the knowledge of the genomic variation of the immunoglobulin loci is scarce. Here, we report 25 potential novel germline IGHV alleles as inferred from rearranged naïve B cell cDNA repertoires of 98 individuals. Thirteen novel alleles were selected for validation, out of which ten were successfully confirmed by targeted amplification and Sanger sequencing of non-B cell DNA. Moreover, we detected a high degree of variability upstream of the V-REGION in the 5′UTR, L-PART1 and L-PART2 sequences, and found that identical V-REGION alleles can differ in upstream sequences. Thus, we have identified a large genetic variation not only in the V-REGION but also in the upstream sequences of IGHV genes. Our findings provide a new perspective for annotating immunoglobulin repertoire sequencing data.

Efficient T cell-B cell collaboration guides autoantibody epitope bias and onset of celiac disease

B cells play important roles in autoimmune diseases through autoantibody production, cytokine secretion, or antigen presentation to T cells. In most cases, the contribution of B cells as antigen-presenting cells is not well understood. We have studied the autoantibody response against the enzyme transglutaminase 2 (TG2) in celiac disease patients by generating recombinant antibodies from single gut plasma cells reactive with discrete antigen domains and by undertaking proteomic analysis of anti-TG2 serum antibodies. The majority of the cells recognized epitopes in the N-terminal domain of TG2. Antibodies recognizing C-terminal epitopes interfered with TG2 cross-linking activity, and B cells specific for C-terminal epitopes were inefficient at taking up TG2-gluten complexes for presentation to gluten-specific T cells. The bias toward N-terminal epitopes hence reflects efficient T-B collaboration. Production of antibodies against N-terminal epitopes coincided with clinical onset of disease, suggesting that TG2-reactive B cells with certain epitope specificities could be the main antigen-presenting cells for pathogenic, gluten-specific T cells. The link between B cell epitopes, antigen presentation, and disease onset provides insight into the pathogenic mechanisms of a T cell-mediated autoimmune condition.

Composition and diversity analysis of the B-cell receptor immunoglobulin heavy chain complementarity-determining region 3 repertoire in patients with acute rejection after kidney transplantation using high-throughput sequencing

The aim of the present study was to assess the genetic diversity of the B-cell receptor (BCR) in kidney transplant recipients with acute rejection. A total of three patients with acute rejection after kidney transplantation were examined by performing a composition and diversity analysis of the BCR immunoglobulin heavy chain (IGH) complementarity-determining region 3 (H-CDR3) repertoire. The peripheral blood mononuclear cells of patients were collected at 1 day prior to (Pre1), as well as 1 day (Post1) and 7 days (Post7) after the transplantation, and DNA was extracted. High-throughput sequencing technology was applied to determine the BCR repertoire. Raw sequences in FASTQ format were analyzed with the Basic Local Alignment Search Tool. The diversity of the BCR repertoire was assessed by calculating Shannon entropy, Simpson's diversity index, the Gini coefficient and highly expanded clone distributions. The diversity of the BCR repertoire at Pre1 was greater than that at Post1 or Post7. The diversity of the BCR repertoire was the lowest at Post1 and increased at Post7 but failed to reach the pre-transplantation levels. Patients exhibited the loss of seven IGH variable (IGHV)3 family genes, while five new genes were expressed at a low frequency. Furthermore, five IGHV-IGH joining (IGHJ) gene pairings, including IGHJ6-IGHV3-11, were detected in the patients. Up- and downregulated genes were assessed by calculating the expression frequencies of the IGH diversity and IGHV gene families at Post1 and Post7. The results of the H-CDR3 length distribution and H-CDR3 amino acid (AA) usage analyses indicated that in Case 1 and 2, the AA length was similar at mostly 14–18 AA, while that in Case 3 was relatively stable at 12–16 AA. In conclusion, the present results illustrate the diversity of H-CDR3 in patients with acute rejection after kidney transplantation may provide novel ideas, methods and means of monitoring and analyzing the immune status of patients under physiological and pathological conditions.

A Functional Idiotype/Anti-Idiotype Network Is Active in Genetically Gluten-Intolerant Individuals Negative for Both Celiac Disease-Related Intestinal Damage and Serum Autoantibodies

An unbalance between Abs that recognize an autoantigen (idiotypes; IDs) and Igs that bind such Abs (anti-IDs) is considered a functional event in autoimmune disorders. We investigated the presence of an ID/anti-ID network in celiac disease (CD), a condition in which antitissue transglutaminase 2 (TG2) Abs are suspected to contribute to CD pathogenesis. To characterize the ID side, we reproduced by in vitro yeast display the intestine-resident Abs from CD and control patients. These TG2-specific IDs were used to identify potential anti-IDs in the serum. We observed elevated titers of anti-IDs in asymptomatic patients with predisposition to CD and demonstrated that anti-ID depletion from the serum restores a detectable humoral response against TG2. Our study provides an alternative approach to quantify CD-related autoantibodies in cases that would be defined "negative serology" with current diagnostic applications. Therefore, we suggest that developments of this technology could be designed for perspective routine tests.

Benchmarking Tree and Ancestral Sequence Inference for B Cell Receptor Sequences

B cell receptor sequences evolve during affinity maturation according to a Darwinian process of mutation and selection. Phylogenetic tools are used extensively to reconstruct ancestral sequences and phylogenetic trees from affinity-matured sequences. In addition to using general-purpose phylogenetic methods, researchers have developed new tools to accommodate the special features of B cell sequence evolution. However, the performance of classical phylogenetic techniques in the presence of B cell-specific features is not well understood, nor how much the newer generation of B cell specific tools represent an improvement over classical methods. In this paper we benchmark the performance of classical phylogenetic and new B cell-specific tools when applied to B cell receptor sequences simulated from a forward-time model of B cell receptor affinity maturation toward a mature receptor. We show that the currently used tools vary substantially in terms of tree structure and ancestral sequence inference accuracy. Furthermore, we show that there are still large performance gains to be achieved by modeling the special mutation process of B cell receptors. These conclusions are further strengthened with real data using the rules of isotype switching to count possible violations within each inferred phylogeny.

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.

Strong Clonal Relatedness between Serum and Gut IgA despite Different Plasma Cell Origins

Mucosal antigens induce generation of lamina propria plasma cells (PCs) that secrete dimeric immunoglobulin A (IgA) destined for transport across the epithelium. In addition, blood contains monomeric IgA. To study the relationship between mucosal and systemic antibody responses, we took advantage of celiac disease patient samples for isolation of gut PCs as well as serum IgA and IgG reactive with a gluten-derived peptide or the autoantigen transglutaminase 2. Proteomic analysis of serum IgA revealed antigen-specific V-gene preferences, which matched those found in gut PCs. Further, gut PC CDR-H3 sequences were abundant in serum IgA but also detectable in serum IgG. Our data indicate that the same B cell clones that give rise to gut PCs also contribute to the serum antibody pool. However, serum IgA antibodies had a molecular composition distinct from that of IgA antibodies secreted in the gut, suggesting that individual B cell clones give rise to different PC populations.

Characterization of the Small Intestinal Lesion in Celiac Disease by Label-Free Quantitative Mass Spectrometry

Global characterization of tissue proteomes from small amounts of biopsy material has become feasible because of advances in mass spectrometry and bioinformatics tools. In celiac disease (CD), dietary gluten induces an immune response that is accompanied by pronounced remodeling of the small intestine. Removal of gluten from the diet abrogates the immune response, and the tissue architecture normalizes. In this study, differences in global protein expression of small intestinal biopsy specimens from CD patients were quantified by analyzing formalin-fixed, paraffin-embedded material using liquid chromatography-mass spectrometry and label-free protein quantitation. Protein expression was compared in biopsy specimens collected from the same patients before and after 1-year treatment with gluten-free diet (n = 10) or before and after 3-day gluten provocation (n = 4). Differential expression of proteins in particular from mature enterocytes, neutrophils, and plasma cells could distinguish untreated from treated CD mucosa, and Ig variable region IGHV5-51 expression was found to serve as a CD-specific marker of ongoing immune activation. In patients who had undergone gluten challenge, coordinated up-regulation of wound response proteins, including the CD autoantigen transglutaminase 2, was observed. Our study provides a global and unbiased assessment of antigen-driven changes in protein expression in the celiac intestinal mucosa.

Evaluation of Antigen-Conjugated Fluorescent Beads to Identify Antigen-Specific B Cells

Selection of single antigen-specific B cells to identify their expressed antibodies is of considerable interest for evaluating human immune responses. Here, we present a method to identify single antibody-expressing cells using antigen-conjugated fluorescent beads. To establish this, we selected Folate Receptor alpha (FRα) as a model antigen and a mouse B cell line, expressing both the soluble and the membrane-bound forms of a human/mouse chimeric antibody (MOv18 IgG1) specific for FRα, as test antibody-expressing cells. Beads were conjugated to FRα using streptavidin/avidin-biotin bridges and used to select single cells expressing the membrane-bound form of anti-FRα. Bead-bound cells were single cell-sorted and processed for single cell RNA retrotranscription and PCR to isolate antibody heavy and light chain variable regions. Variable regions were then cloned and expressed as human IgG1/k antibodies. Like the original clone, engineered antibodies from single cells recognized native FRα. To evaluate whether antigen-coated beads could identify specific antibody-expressing cells in mixed immune cell populations, human peripheral blood mononuclear cells (PBMCs) were spiked with test antibody-expressing cells. Antigen-specific cells could comprise up to 75% of cells selected with antigen-conjugated beads when the frequency of the antigen-positive cells was 1:100 or higher. In PBMC pools, beads conjugated to recombinant antigens FRα and HER2 bound antigen-specific anti-FRα MOv18 and anti-HER2 Trastuzumab antibody-expressing cells, respectively. From melanoma patient-derived B cells selected with melanoma cell line-derived protein-coated fluorescent beads, we generated a monoclonal antibody that recognized melanoma antigen-coated beads. This approach may be further developed to facilitate analysis of B cells and their antibody profiles at the single cell level and to help unravel humoral immune repertoires.

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.

Mechanisms underlying B cell immune dysregulation and autoantibody production in MuSK myasthenia gravis

Pathogenic autoantibodies to muscle-specific tyrosine kinase (MuSK) can be found in patients with myasthenia gravis (MG) who do not have detectable antibodies to the acetylcholine receptor. Although the autoantibody-mediated pathology is well understood, much remains to be learned about the cellular immunology that contributes to autoantibody production. To that end, our laboratory has investigated particular components associated with the cellular immunopathology of MuSK MG. First, we found that B cell tolerance defects contribute to the abnormal development of the naive repertoire, which indicates that dysregulation occurs before the production of autoantibodies. Second, both the naive and antigen-experienced memory B cell repertoire, which we examined through the application of high-throughput adaptive immune receptor repertoire sequencing, include abnormalities not found in healthy controls. This highlights a broad immune dysregulation. Third, using complementary approaches, including production of human monoclonal antibodies, we determined that circulating plasmablasts directly contribute to the production of MuSK-specific autoantibodies in patients experiencing relapse following B cell depletion therapy. These collective findings contribute to defining a mechanistic model that describes MuSK MG immunopathogenesis.

New Insights into the Pathogenesis of Celiac Disease

Celiac disease (CD) is an autoimmune and multisystem gluten-related disorder that causes symptoms involving the gastrointestinal tract and other organs. Pathogenesis of CD is only partially known. It had been established that ingestion of gluten proteins present in wheat and other cereals are necessary for the disease and develops in individuals genetically predisposed carrying the DQ2 or DQ8 human leukocyte antigen haplotypes. In this review, we had pay specific attention on the last discoveries regarding the three cellular components mainly involved in the development and maintenance of CD: T-cells, B-cells, and microbioma. All of them had been showed critical for the interaction between inflammatory immune response and gluten peptides. Although the mechanisms of interaction among overall these components are not yet fully understood, recent proteomics and molecular studies had shed some lights in the pathogenic role of tissue transglutaminase 2 in CD and in the alteration of the intestinal barrier function induced by host microbiota.