Peripheral IgE Repertoires of Healthy Donors Carry Moderate Mutation Loads and Do Not Overlap With Other Isotypes

IgE-mediated allergic disease represents an increasing health problem. Although numerous studies have investigated IgE sequences in allergic patients, little information is available on the healthy IgE repertoire. IgM, IgG, IgA, and IgE transcripts from peripheral blood B cells of five healthy, non-atopic individuals were amplified by unbiased, template-switching, isotype-specific PCR. Complete VDJ regions were sequenced to near-exhaustion on the PacBio platform. Sequences were analyzed for clonal relationships, degree of somatic hypermutation, IGHV gene usage, evidence of antigenic selection, and N-linked glycosylation motifs. IgE repertoires appeared to be highly oligoclonal with preferential usage of certain IGHV genes compared to the other isotypes. IgE sequences carried more somatic mutations than IgM, yet fewer than IgG and IgA. Many IgE sequences contained N-linked glycosylation motifs. IgE sequences had no clonal relationship with the other isotypes. The IgE repertoire in healthy individuals is derived from relatively few clonal expansions without apparent relations to immune reactions that give rise to IgG or IgA. The mutational burden of normal IgE suggests an origin through direct class-switching from the IgM repertoire with little evidence of antigenic drive, and hence presumably low affinity for specific antigens. These findings are compatible with a primary function of the healthy IgE repertoire to occupy Fcε receptors for competitive protection against mast cell degranulation induced by allergen-specific, high-affinity IgE. This background knowledge may help to elucidate pathogenic mechanisms in allergic disease and to design improved desensitization strategies.

Introduction

Type I hypersensitivity is an immune response triggered by allergen-specific IgE. Binding of allergen-bound IgE to Fcε receptors on mast cells leads to prompt degranulation, which provokes a range of clinical symptoms including atopic dermatitis, asthma, allergic rhinoconjunctivitis, urticarial, and anaphylaxis (1, 2). Despite the recognition that IgE-mediated diseases are becoming an ever increasing health burden, especially in urban societies (3), the mechanisms leading to IgE-mediated disease, as well as the underlying principles causing clinical heterogeneity, remain incompletely understood (4, 5). An improved understanding of the development of IgE repertoires in healthy individuals may aid in the identification of these disease mechanisms and facilitates efficient design of anti-allergic strategies (6, 7).

Two distinct pathways to generate IgE have been identified in mouse models. In the direct pathway, functional B-cell receptor (BCR) genes undergo direct class-switch recombination (CSR) from IgM to IgE isotype and yield predominantly IgE with low affinity to antigen. In the indirect pathway, high-affinity IgE results from secondary CSR in B cells that express IgG or potentially IgA (8, 9). Along with increased affinity, IgE derived from IgG1 B cells also carry significantly more mutations than IgE in IgG1-deficient mice (10). The existence of multiple pathways has been further supported by flow cytometric analyses in humans showing subsets both dependent and independent of germinal centers (11).

Since low and high affinity IgE compete for occupation of the Fcε receptor, higher concentrations of low-affinity IgE than high-affinity IgE can theoretically provide protection against anaphylaxis (10). In non-allergic humans, very limited information on the extent and origin of the IgE repertoire is currently available and consists of only 60 near-full-length IgE sequences from two individuals (12). Although recently a number of studies have applied massive parallel sequencing to gain new insights in the IgE repertoire, none of them generated the full-length sequences essential for comprehensive analysis and/or used primer binding bias-free methodology (13–18).

This study aimed to provide in-depth characterization of the IgE repertoire in healthy, non-allergic individuals as an essential reference for comparative studies in allergic and desensitized individuals. We determined the peripheral blood IgE BCR repertoire of five non-allergic donors to near completion by unbiased, full-length massive parallel sequencing. IgE BCR repertoires were compared to IgM, IgG, and IgA repertoires to support the hypothesis that the direct pathway of IgE generation would be dominant in non-allergic human individuals.

Materials and Methods

Material Collection and Storage

Cryopreserved aliquots of Ficoll-separated mononuclear cells (PBMC) from five healthy, asymptomatic stem cell donors were obtained from the biobank of the Leiden University Medical Center biobank in accordance with local guidelines. Absence of atopic constitution was confirmed by measurement of total IgM, IgG, IgA, and IgE immunoglobulin levels in time-matched serum samples.

Flow Cytometry Analysis for Expression of CD19 and IgE

IgE-expressing B cells were isolated from aliquots of 1 × 10⁵ thawed PBMC by flow cytometry. To avoid artifacts from binding of IgE to Fcε receptor-expressing cells (19), cells were fixed in 200 μl 1% paraformaldehyde in phosphate-buffered saline (B. Braun, Melsungen, Germany) with 5% fetal bovine serum (Bodinco, Alkmaar, The Netherlands) and 0.2% saponin (Sigma Aldrich, St. Louis, MO, USA) for 8 min at 4°C. After washing, fixed cells were permeabilized for 30 min at 4°C in 200 μl of the same buffer without paraformaldehyde. Cells were washed once more and stained with anti-CD19-FITC (BD Biosciences, Franklin Lakes, NJ, USA) and anti-IgE-APC (Miltenyi Biotec, Leiden, The Netherlands) for 30 min at 4°C in the dark, followed by another wash and resuspension. The abundance of CD19⁺ and IgE⁺ cells was determined by flow cytometry in live gated cells.

BCR Repertoire Sequencing

The ARTISAN PCR protocol for unbiased amplification of BCR repertoires (20) was adapted for IgE by designing a series of IgE constant region-specific reverse primers (ε.rt 5′-GGCATAGTGACCAGAGAGCG-3′ for reverse transcription; ε.pcr1 5′-GGTCACCATCACCGGCTCCG-3′ for initial PCR amplification; ε.pcr2 5′-GGCAGCCCAGAGTCACGG-3′ for semi-nested PCR amplification; ε.bc 5′-[barcode]-CGGATGGGCTCTGTGTGG-3′ for barcoding; barcodes: 5′-CCATCTCATATGTAGTACTCT-3′, 5′-CGGATGGGCTCTGTGTGG-3′, 5′-CGGATGGGCTCTGTGTGG-3′, 5′-CGGATGGGCTCTGTGTGG-3′ and 5′-CGGATGGGCTCTGTGTGG-3′). B cells were isolated from aliquots of thawed cells by removal of non-B cells with magnetic beads (B cell isolation kit II; Miltenyi Biotec, Leiden, The Netherlands) and routinely yielded a purity of >99% CD19⁺ B cells as assessed by flow cytometry. For each healthy donor, 2 × 10⁶ B cells were divided into five aliquots.

Messenger RNA isolation and cDNA synthesis were performed separately for each aliquot with addition of the ε.rt primer to the ARTISAN cDNA synthesis mix (20). IgM, IgG, and IgA cDNA was amplified according to the original ARTISAN PCR protocol. Due to the low abundance of IgE⁺ B cells in healthy donor peripheral blood, amplification of functional IgE transcripts was extensively optimized with different primer combinations on serial dilutions (100–0.01%) of the IgE-expressing multiple myeloma cell line U266 (DSMZ, Braunschweig, Germany) in a background of healthy donor PBMCs. For IgE repertoires, cDNA was first amplified for 15 cycles with the ε.pcr1 primer. First-round IgE amplicons were purified by silica spin columns (Promega, Madison, WI, USA) and re-amplified for 15 cycles of semi-nested PCR using the ε.pcr2 primer. Libraries were barcoded at the 5′ terminus to identify the donor and at the 3′ terminus to identify individual aliquots. Pooled libraries were amplified as single molecules in rolling circles on a total of five SMRT cells on the RSII system (Pacific Biosciences, Menlo Park, CA, USA). IgE sequences were sequenced to high depth to achieve near-complete representation of all IgE⁺ B cells present in the sample.

Sequence and Statistical Analysis

Output sequence files were filtered with SMRT portal software for a minimum of eight sequencing passes. All sequences were annotated by IMGT HighV-QUEST (21). For statistical analyses, sequences with identical IGHV genes and amino acid CDR3 sequences within one aliquot were counted as a single sequence.

IGHV gene usage was compared between isotypes by Fisher's exact test and corrected for multiple testing by Bonferroni correction. Cumulative differences in IGHV gene usage were calculated by determining the absolute differences in fractional IGHV usage between two isotypes for every IGHV gene and adding these to a cumulative difference with a theoretical maximum of 200. BCR mutational status and CDR3 length were compared between isotypes by unpaired t-test.

N-linked glycosylation motifs were identified as Asn-X-Ser/Thr motifs (where X may be any amino acid except proline) and their abundance, location, and the mutational status of the corresponding sequence were compared between isotypes using unpaired t and Fisher's exact tests.

Clonal B-cell expansions were defined as the presence of BCR sequences with identical V, D and J gene usage, identical CDR3 length and ≥95% nucleotide overlap in CDR3 in either multiple aliquots or in more than one isotype of any individual donor. Intraclonal sequence variation was defined as the distance from the clonal consensus and determined for the largest 10 IgE clones, as well as the largest 5 IgM, IgG, and IgA clones.

Calculation of Sampling Depth

To assess how many cells from the sample were represented in the VDJ sequence libraries, we performed an in silico simulation to estimate the fraction of observed unique sequences sampled out of a large pool of sequences. A unique numerical identifier was assigned to each unique observed BCR sequence for each donor, and 5 × 10⁵ copies of every identifier were pooled in the simulation. Consequently, random sequences were sampled from the pool up to the number of sequences obtained from massive parallel sequencing, and the number of unique identifiers was counted. This simulation was performed 100 times per donor, and the median number of unique samples identifiers per donor was calculated.

Results

According to applicable stem cell donor regulations, all donors had no atopic constitution, denied any allergic symptoms, and had normal serum immunoglobulin and IgE concentrations (Table 1).

TABLE 1

Table 1. Serum immunoglobulin concentrations in five healthy donors.

Sequence Acquisition

A median of 0.09% of peripheral blood B cells of healthy donors expressed IgE (Table 2), corresponding to a median of 2,000 cells per 2 × 10⁶ B cells. A median of 3,254 total full-length, potentially functional VDJ IgE sequences were obtained per donor by ARTISAN PCR and deep PacBio sequencing, representing a median of 1.52 IgE sequences per IgE⁺ B cell present in the sample. In silico simulations indicated that this amount of oversampling would theoretically return a median of 78% of unique sequences (Supplementary Data). In addition, medians of 2,088 IgM (range: 1,287–2,437), 1,221 IgG (range: 735–2,144), and 2,770 IgA (range: 1,351–3,012) unique, full-length and potentially functional VDJ sequences were obtained per donor.

TABLE 2

Table 2. Quantification of IgE⁺ B cells and sequencing depth.

Clonal B-Cell Expansions

To identify and quantify clonal B cell expansions, we first assigned all individual VDJ sequences within an isotype of each donor to clonotypes according to stringent criteria developed rationally to minimize the calling of false-positive clonal relationships (Supplementary Data). This analysis identified a median of 39 (range: 23–46) putative unique IgE clonotypes per donor and indicated representation of each clonotype by an average of 85 closely related VDJ sequences. Additional manual inspection revealed that some of these putative clonotypes shared substantial numbers of individual mutations in their IGHV, strongly suggesting a common clonal origin despite less than 95% CDR3 identity and despite occasional variation in CD3 length. Therefore, we combined putative clonotypes from an individual that shared at least 70% of their IGHV mutations regardless of CDR3 similarities into definitive clonotypes for all further analyses. We identified a total of 146 definitive IgE clonotypes (median per donor 31; range: 18–38). In a single instance, these criteria indicated the presence of a single clonotype in two donors. Since the libraries of these two donors were sequenced on the same SMRT cell, this phenomenon may originate from barcode contamination and does not permit to conclude the presence of canonical IgE clonotypes across individuals.

Expansion of an IgE-expressing B-cell clone was unequivocally demonstrated by the presence of 69 of the total of 146 IgE clonotypes (47%) in at least two B-cell aliquots (Figure 1), corresponding to a median of 12 (range: 11–19) clonal expansions of IgE-expressing B-cells per donor. In comparison, only 1.3 of IgM, 7.9 of IgG, and 9.1% of IgA VDJ clonotypes were found in multiple aliquots and therefore derived from clonally expanded B cells, corresponding to a median of 14 (range: 6–40; total 97), 46 (range: 28–128; total 301) and 115 (range: 54–192; total 625) clonal expansions per donor and isotype, respectively (Figure 1). IgE clonal expansions were fewer than IgG (p = 0.034) and IgA (p = 0.0017), but similar to IgM expansions (p = 0.39). Intraclonal sequence diversity was lower in IgE than in all other isotypes (IgM: p = 0.0002; IgG: p = 0.0184; IgA: p = 0.0039; Figure 2). Twelve VDJ clonotypes were present in both IgM and IgG BCR repertoires, 24 in IgM and IgA, and 45 in IgG and IgA. Seven clonotypes comprised IgM, IgG, and IgA isotypes. In marked contrast, no relationship was found between any IgE and non-IgE VDJ. Since this finding is in contrast with previous reports (22), we performed extensive in silico simulations that demonstrated frequent false-positive conclusions of clonal overlap when inappropriately relaxed criteria for clonal relationship are applied (Supplementary Data).

FIGURE 1

Figure 1. Clonality in the healthy BCR repertoire. Distribution of unique VDJ sequences present in one or more cellular aliquots per isotype and donor.

FIGURE 2

Figure 2. Intraclonal sequence diversity in the healthy BCR repertoire. Intraclonal variability per isotype. The number of nucleotide differences compared to the respective consensus VDJ sequence is displayed for every VDJ sequence of the 5 largest B-cell clones expressing IgM, IgG, or IgA, respectively, and for the 10 largest IgE-expressing B-cell clones.

BCR Characteristics per Isotype

IgE VDJ utilized more frequently IGHV3-9 and IGHV3-11 (p < 0.0001), and less frequently IGHV4-39 than the remaining isotypes (p = 0.0018; Figure 3A). Concordantly, cumulative gene fraction distances between IgE and other isotypes were significantly larger than distances between IgM and IgG, IgM, and IgA, and IgG, and IgA (Figure 3B).

FIGURE 3

Figure 3. Differential IGHV usage between isotypes. (A) Heat map indicating the relative frequency of all IGHV genes among unique VDJ sequences per isotype. For each isotype, the frequencies shown add up to 100%. (B) Pairwise comparison of cumulative IGHV gene distance between two isotypes. For each IGHV gene, the fractions of VDJ sequences containing that IGHV gene within an isotype were calculated, and the numerical difference of these fractions between two isotypes was determined per IGHV gene. For each pairwise comparison between isotypes, all differences of IGHV gene fractions were added to obtain their cumulative distance.

All isotypes had median CDR3 lengths of 15 amino acids (Figure 4A). Unique IgE VDJ carried fewer mutations (4.9%) than the other class-switched isotypes (p < 0.001) (Figure 4B).

FIGURE 4

Figure 4. Parameters of BCR hypervariability. (A) Comparison of CDR3 length between isotypes. The number of unique VDJ sequences with a given CDR3 amino acid (AA) length as defined by IGMT HighV-Quest is shown per isotype. (B) Comparison of somatic hypermutation between isotypes. The number of unique VDJ sequences with a given number of nucleotide changes compared to their most closely related germ-line IGHV sequence is displayed per isotype. Sequences longer than 75 bp are indicated by “>”. (C) Fractions of VDJ sequences with N-glycosylation motifs per isotype. Fractions indicate the origin of the N-glycosylation motifs and whether a germ-line (GL)-encoded motif was lost by somatic hypermutation (SHM).

Despite their moderate SHM burden, more IgE VDJ (23%) had acquired N-linked glycosylation motifs (NLGM) through SHM than IgG (16%; p < 0.001) and IgA (17%; p < 0.001). When considering all NLGM, i.e., both germline-encoded (from IGHV1-8, IGHV4-34, and IGHV5-10-1) and SHM-derived motifs, their prevalence of 27% in all IgE VDJ was significantly higher (p < 0.0001) than in IgM, IgG, and IgA (Figure 4C). This difference was largely attributable to preservation of germ-line-encoded NLGM in the IgE compartment (Table 3).

TABLE 3

Table 3. Preservation of germline-encoded N-linked glycosylation sites.

Discussion

In conclusion, we provide the first comprehensive inventory of near-complete peripheral blood IgE repertoires from healthy individuals through tailored methodology that lacks primer binding bias, yields full-length VDJ sequences, and detects clonal expansions by standardized parallel analysis of several aliquots (20). The applicable regulations for volunteer stem cell donors precluded the acquisition of additional epidemiological information such as dwelling and other living conditions that are associated with allergy. The observed frequency of 0.09% IgE-expressing cells among peripheral B cells is on the high end of the reported spectrum (23, 24). The consequences of a possible overestimation of the true prevalence of IgE+ B cells would predominantly imply that the sequencing of the IgE repertoires would have been even more exhaustive than indicated by our simulations. In this context, presence of clonally related IgE sequences in several cellular aliquots is a much more reliable indicator of clonal expansions than BCR sequence read counts in massive parallel sequencing experiments. Nevertheless, it is a striking observation that IgE sequences found in all 5 aliquots dominated the sequence libraries (not shown), suggesting that few IgE+ clones actually had expanded strongly.

The presence of SHM indicates that IgE⁺ B cells have passed through germinal center reactions. The intermediate SHM load of IgE between IgM and IgG/IgA could result from direct CSR of the majority of healthy donor PBMC from IgM to IgE (13). Alternatively, non-IgE B cells could have acquired higher SHM loads by repeated GC passages, whereas IgE⁺ B cells have only a limited presence in germinal centers, (25–27). The difference in IgE mutation rate compared to other isotypes was more striking than in another recent study of IgE in non-allergic subjects, yet not as low as in children with atopic dermatitis (15, 16). These findings call for further verification.

The striking lack of clonal relationships between IgE-expressing B cells and B cells expressing the other isotypes supports the important conclusion that IgE-expressing B cells emerge from qualitatively different immune responses. While the IgE repertoire has been sequenced to apparent exhaustion in our study, incomplete sampling of the other isotype compartments cannot completely exclude a low degree of clonal overlap. However, lack of clonal overlap is corroborated by a recent massive parallel sequencing study that also found virtually no clonal relationship between IgE and other isotypes in healthy donors (16). In contrast, marked overlap of allergen-specific IgE clones with IgG and IgA has been described in allergic individuals (14, 17, 28) and in immunized mice (13).

The distinct characteristics of the normal IgE repertoire, i.e., relatively low SHM burden, lack of intraclonal sequence variation, skewed usage of IGHV genes, marked oligoclonality, retention of germline-encoded NLGM, frequent acquisition of additional such motifs, and striking absence of clonal relatedness to IgM-expressing and non-IgE-class-switched B cells cumulatively indicates its origin from qualitatively different immune responses than IgG- and IgA-expressing B cells. Like allergic individuals (15, 17), healthy donors appear to have IgE repertoires composed of a small pool of highly expanded clones. However, in healthy donors, these likely represent rearrangements in a low affinity, non-antigen-specific, “static” state. As previously suggested, such low-affinity IgE could have a protective role against allergy by competition with high-affinity, type I hypersensitivity-inducing IgE for occupation of Fcε receptors (13, 29).

The novel hypothesis of a specific role of N-glycosylation in expansion and maintenance of IgE-expressing B cells in non-allergic individuals warrants further studies. Ubiquitous NLGM acquisition could have a disease-specific role in selection of IgE-expressing B-cells in similarity to follicular lymphoma (30), primary cutaneous follicle center lymphoma (31) and rheumatoid arthritis (32). On the other hand, glycosylation at NLGM positions could non-specifically obstruct antigen recognition and effectively inhibit BCR affinity maturation (33).

Our findings also support direct class switch recombination from IgM to IgE as the origin of such low-affinity IgE antibodies as demonstrated in mouse studies (10, 13, 19).

Some, but not all (17, 18, 28), previous studies found preferential IGHV usage in allergic individuals, most notably of IGHV2, IGHV4, and IGHV5 family genes (15, 34–36). Our data do not indicate this particular IGHV bias in healthy individuals. Although theoretically attributable to study population differences, primer binding bias in multiplexed primer strategies creates an inherent risk of skewing repertoire analyses. Use of a forward primer binding to an artificial uniform sequence at the 5′ cDNA termini effectively alleviates this risk (7, 20, 31, 37). Differential amplification efficiency of different IGHV genes may be another theoretical source of bias for observed IGHV usage. However, IgE sequences underwent only two additional thermocycles than the other isotypes, effectively limiting this particular risk. In addition, ARTISAN PCR employs long extension times to avoid such bias. In previous applications of this method, no preferential V allele amplification was actually observed (20).

In comparison to allergic individuals and derivatives of high-affinity murine clones, healthy donors appear to carry IgE sequences with fewer BCR mutations (12–15, 18, 28) and less restriction to few individual clonotypes (7, 12, 38, 39). Although IgE VDJ from allergic and parasite-exposed patients also lacked evidence of substantial antigenic selection (7, 15), these characteristics indicate that high-affinity (allergen-specific) IgE clones in allergic individuals are probably generated through indirect class switch recombination from IgG and IgA clonal expansions, preferentially utilizing a restricted repertoire of IGHV genes (13).

Overall, our results add various new aspects to the current knowledge on the IgE repertoire (7). Future projects investigating IgE repertoires in allergic disease should be generated by the same high-standard unbiased approach used here to allow for side-to-side comparison with healthy donors.

Ethics Statement

This study was carried out in accordance with the guidelines as outlined by the LUMC Biobank Committee with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the LUMC Medical Ethical Testing Committee.

Author Contributions

CvB obtained samples. MK and IT obtained the data. MK analyzed the data. DA, MN, and SK provided bioinformatics support. MK and HV wrote the manuscript. All authors read and approved the manuscript.

Funding

MK was supported by a grant from the Professor Steenhuis Fund. MN and DA were supported by Fondecyt 11140542 and Fondecyt 1180882.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2019.01543/full#supplementary-material

References

1. Akdis CA. Allergy and hypersensitivity: mechanisms of allergic disease. Curr Opin Immunol. (2006) 18:718–26. doi: 10.1016/j.coi.2006.09.016

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Galli SJ, Tsai M, Piliponsky AM. The development of allergic inflammation. Nature. (2008) 454:445–54. doi: 10.1038/nature07204

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Ring J, Akdis C, Behrendt H, Lauener RP, Schappi G, Akdis M, et al. Davos declaration: allergy as a global problem. Allergy. (2012) 67:141–3. doi: 10.1111/j.1398-9995.2011.02770.x

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Akdis CA. Therapies for allergic inflammation: refining strategies to induce tolerance. Nat Med. (2012) 18:736–49. doi: 10.1038/nm.2754

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Lotvall J, Akdis CA, Bacharier LB, Bjermer L, Casale TB, Custovic A, et al. Asthma endotypes: a new approach to classification of disease entities within the asthma syndrome. J Allergy Clin Immunol. (2011) 127:355–60. doi: 10.1016/j.jaci.2010.11.037

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Akdis M, Akdis CA. Mechanisms of allergen-specific immunotherapy: multiple suppressor factors at work in immune tolerance to allergens. J Allergy Clin Immunol. (2014) 133:621–31. doi: 10.1016/j.jaci.2013.12.1088

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Gadermaier E, Levin M, Flicker S, Ohlin M. The human IgE repertoire. Int Arch Allergy Immunol. (2014) 163:77–91. doi: 10.1159/000355947

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Erazo A, Kutchukhidze N, Leung M, Christ AP, Urban JF Jr, Curotto de. Lafaille MA, et al. Unique maturation program of the IgE response in vivo. Immunity. (2007) 26:191–203. doi: 10.1016/j.immuni.2006.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Mills FC, Thyphronitis G, Finkelman FD, Max EE. Ig mu-epsilon isotype switch in IL-4-treated human B lymphoblastoid cells. evidence for a sequential switch. J Immunol. (1992) 149:1075–85.

PubMed Abstract | Google Scholar

10. Xiong H, Dolpady J, Wabl M, Curotto de. Lafaille MA, Lafaille JJ. Sequential class switching is required for the generation of high affinity IgE antibodies. J Exp Med. (2012) 209:353–64. doi: 10.1084/jem.20111941

CrossRef Full Text | Google Scholar

11. Berkowska MA, Heeringa JJ, Hajdarbegovic E, van der Burg M, Thio HB, van Hagen PM, et al. Human IgE(+) B cells are derived from T cell-dependent and T cell-independent pathways. J Allergy Clin Immunol. (2014) 134:688–97 e6. doi: 10.1016/j.jaci.2014.03.036

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Dahlke I, Nott DJ, Ruhno J, Sewell WA, Collins AM. Antigen selection in the IgE response of allergic and nonallergic individuals. J Allergy Clin Immunol. (2006) 117:1477–83. doi: 10.1016/j.jaci.2005.12.1359

PubMed Abstract | CrossRef Full Text | Google Scholar

13. He JS, Subramaniam S, Narang V, Srinivasan K, Saunders SP, Carbajo D, et al. IgG1 memory B cells keep the memory of IgE responses. Nat Commun. (2017) 8:641. doi: 10.1038/s41467-017-00723-0

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Hoh RA, Joshi SA, Liu Y, Wang C, Roskin KM, Lee JY, et al. Single B-cell deconvolution of peanut-specific antibody responses in allergic patients. J Allergy Clin Immunol. (2016) 137:157–67. doi: 10.1016/j.jaci.2015.05.029

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Kerzel S, Rogosch T, Struecker B, Maier RF, Kabesch M, Zemlin M. Unlike in children with allergic asthma, IgE transcripts from preschool children with atopic dermatitis display signs of superantigen-driven activation. J Immunol. (2016) 196:4885–92. doi: 10.4049/jimmunol.1402889

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Kitaura K, Yamashita H, Ayabe H, Shini T, Matsutani T, Suzuki R. Different somatic hypermutation levels among antibody subclasses disclosed by a new next-generation sequencing-based antibody repertoire analysis. Front Immunol. (2017) 8:389. doi: 10.3389/fimmu.2017.00389

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Levin M, King JJ, Glanville J, Jackson KJ, Looney TJ, Hoh RA, et al. Persistence and evolution of allergen-specific IgE repertoires during subcutaneous specific immunotherapy. J Allergy Clin Immunol. (2016) 137:1535–44. doi: 10.1016/j.jaci.2015.09.027

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Levin M, Levander F, Palmason R, Greiff L, Ohlin M. Antibody-encoding repertoires of bone marrow and peripheral blood-a focus on IgE. J Allergy Clin Immunol. (2017) 139:1026–30. doi: 10.1016/j.jaci.2016.06.040

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Wesemann DR, Magee JM, Boboila C, Calado DP, Gallagher MP, Portuguese AJ, et al. Immature B cells preferentially switch to IgE with increased direct Smu to Sepsilon recombination. J Exp Med. (2011) 208:2733–46. doi: 10.1084/jem.20111155

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Koning MT, Kielbasa SM, Boersma V, Buermans HPJ, van der Zeeuw SAJ, van Bergen CAM, et al. ARTISAN PCR: rapid identification of full-length immunoglobulin rearrangements without primer binding bias. Br J Haematol. (2017) 178:983–6. doi: 10.1111/bjh.14180

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Alamyar E, Duroux P, Lefranc MP, Giudicelli V. IMGT((R)) tools for the nucleotide analysis of immunoglobulin (IG) and T cell receptor (TR) V-(D)-J repertoires, polymorphisms, and IG mutations: IMGT/V-QUEST and IMGT/HighV-QUEST for NGS. Methods Mol Biol. (2012) 882:569–604. doi: 10.1007/978-1-61779-842-9_32

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Looney TJ, Lee JY, Roskin KM, Hoh RA, King J, Glanville J, et al. Human B-cell isotype switching origins of IgE. J Allergy Clin Immunol. (2016) 137:579–86 e7. doi: 10.1016/j.jaci.2015.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

23. King CL, Nutman TB. IgE and IgG subclass regulation by IL-4 and IFN-gamma in human helminth infections. assessment by B cell precursor frequencies. J Immunol. (1993) 151:458–65.

PubMed Abstract | Google Scholar

24. King CL, Poindexter RW, Ragunathan J, Fleisher TA, Ottesen EA, Nutman TB. Frequency analysis of IgE-secreting B lymphocytes in persons with normal or elevated serum IgE levels. J Immunol. (1991) 146:1478–83.

Google Scholar

25. He JS, Meyer-Hermann M, Xiangying D, Zuan LY, Jones LA, Ramakrishna L, et al. The distinctive germinal center phase of IgE+ B lymphocytes limits their contribution to the classical memory response. J Exp Med. (2013) 210:2755–71. doi: 10.1084/jem.20131539

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Talay O, Yan D, Brightbill HD, Straney EE, Zhou M, Ladi E, et al. IgE(+) memory B cells and plasma cells generated through a germinal-center pathway. Nat Immunol. (2012) 13:396–404. doi: 10.1038/ni.2256

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Yang Z, Sullivan BM, Allen CD. Fluorescent in vivo detection reveals that IgE(+) B cells are restrained by an intrinsic cell fate predisposition. Immunity. (2012) 36:857–72. doi: 10.1016/j.immuni.2012.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Croote D, Darmanis S, Nadeau KC, Quake SR. High-affinity allergen-specific human antibodies cloned from single IgE B cell transcriptomes. Science. (2018) 362:1306–9. doi: 10.1126/science.aau2599

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Nopp A, Johansson SG, Ankerst J, Bylin G, Cardell LO, Gronneberg R, et al. Basophil allergen threshold sensitivity: a useful approach to anti-IgE treatment efficacy evaluation. Allergy. (2006) 61:298–302. doi: 10.1111/j.1398-9995.2006.00987.x

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Schneider D, Duhren-von Minden M, Alkhatib A, Setz C, van Bergen CA, Benkisser-Petersen M, et al. Lectins from opportunistic bacteria interact with acquired variable-region glycans of surface immunoglobulin in follicular lymphoma. Blood. (2015) 125:3287–96. doi: 10.1182/blood-2014-11-609404

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Koning MT, Nteleah V, Veelken H, Navarrete MA. Template-switching anchored polymerase chain reaction reliably amplifies functional lambda light chain transcripts of malignant lymphoma. Leuk Lymphoma. (2014) 55:1212–4. doi: 10.3109/10428194.2013.828353

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Vergroesen RD, Slot LM, Hafkenscheid L, Koning MT, van der Voort EIH, Grooff CA, et al. B-cell receptor sequencing of anti-citrullinated protein antibody (ACPA) IgG-expressing B cells indicates a selective advantage for the introduction of N-glycosylation sites during somatic hypermutation. Ann Rheum Dis. (2017) 77:956–8. doi: 10.1136/annrheumdis-2017-212052

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Wright A, Tao MH, Kabat EA, Morrison SL. Antibody variable region glycosylation: position effects on antigen binding and carbohydrate structure. EMBO J. (1991) 10:2717–23. doi: 10.1002/j.1460-2075.1991.tb07819.x

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Coker HA, Harries HE, Banfield GK, Carr VA, Durham SR, Chevretton E, et al. Biased use of VH5 IgE-positive B cells in the nasal mucosa in allergic rhinitis. J Allergy Clin Immunol. (2005) 116:445–52. doi: 10.1016/j.jaci.2005.04.032

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Snow RE, Djukanovic R, Stevenson FK. Analysis of immunoglobulin E VH transcripts in a bronchial biopsy of an asthmatic patient confirms bias towards VH5, and indicates local clonal expansion, somatic mutation and isotype switch events. Immunology. (1999) 98:646–51. doi: 10.1046/j.1365-2567.1999.00910.x

PubMed Abstract | CrossRef Full Text | Google Scholar

36. van der Stoep N, van der Linden J, Logtenberg T. Molecular evolution of the human immunoglobulin E response: high incidence of shared mutations and clonal relatedness among epsilon VH5 transcripts from three unrelated patients with atopic dermatitis. J Exp Med. (1993) 177:99–107. doi: 10.1084/jem.177.1.99

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Wang Y, Jackson KJ, Chen Z, Gaeta BA, Siba PM, Pomat W, et al. IgE sequences in individuals living in an area of endemic parasitism show little mutational evidence of antigen selection. Scand J Immunol. (2011) 73:496–504. doi: 10.1111/j.1365-3083.2011.02525.x

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Edwards MR, Brouwer W, Choi CH, Ruhno J, Ward RL, Collins AM. Analysis of IgE antibodies from a patient with atopic dermatitis: biased V gene usage and evidence for polyreactive IgE heavy chain complementarity-determining region 3. J Immunol. (2002) 168:6305–13. doi: 10.4049/jimmunol.168.12.6305

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Janezic A, Chapman CJ, Snow RE, Hourihane JO, Warner JO, Stevenson FK. Immunogenetic analysis of the heavy chain variable regions of IgE from patients allergic to peanuts. J Allergy Clin Immunol. (1998) 101:391–6. doi: 10.1016/S0091-6749(98)70253-2

PubMed Abstract | CrossRef Full Text | Google Scholar