- 1United States Department of Agriculture, Agricultural Research Service, Southern Regional Research Center, New Orleans, LA, United States
- 2United States Department of Agriculture, Agricultural Research Service -Genomics and Bioinformatics Research Unit, Gainesville, FL, United States
- 3Aimmune Nestlé Health Science, Brisbane, CA, United States
- 4Delgado Community College, New Orleans, LA, United States
- 5Department of Pediatrics, University of North Carolina School of Medicine, Chapel Hill, NC, United States
Background: Oral immunotherapy (OIT) with peanut (Arachis hypogaea) allergen powder-dnfp (PTAH; Aimmune Therapeutics) is an FDA-approved treatment to desensitize peanut allergic participants.
Objective: Here we assessed shifts in IgE and IgG4 binding to peanut allergens and their epitopes recognized by United States (US) peanut allergic participants (n = 20) enrolled in phase 3 PTAH OIT clinical trials.
Methods: Pre- and post- trial participant sera were collected approximately 12 months apart and tested for IgE binding to intact peanut proteins via ImmunoCAP ISAC immunoassays. IgE and IgG4 linear epitopes were identified based on binding to synthetic overlapping 15-mer linear peptides of 10 peanut allergens (Ara h 1-11) synthesized on microarray slides.
Results: Statistically significant decreases in IgE binding were identified for intact Ara h 2, 3, and 6, and known and newly identified IgE epitopes were shown to exhibit shifts towards IgG4 binding post-OIT, with most linear peptides having increased IgG4 binding after treatment with PTAH. While PTAH does not seem to alter the actual peptide binding patterns significantly after one year of treatment, the IgE and IgG4 binding ratios and intensity are altered.
Conclusion: At a population level, the linear IgE and IgG4 epitopes of 10 peanut allergens overlap and that increase in IgG4 with OIT results in displacement of IgE binding to both conformational and linear epitopes. Furthermore, it appears as though the increase in IgG4 is more important to achieve desensitization at the 12-month timepoint than the decrease in IgE. This type of knowledge can be useful in the identification of IgE and IgG4-binding allergen and peptide biomarkers that may indicate desensitization or sustained unresponsiveness of allergic individuals to peanut.
1. Introduction
Peanut allergy is an immunologic disease characterized by the production of peanut-protein-targeting immunoglobulin E (IgE) antibodies that mediate reactions ranging from mild cutaneous manifestations to life-threatening systemic anaphylaxis. The risk of life-threatening allergic reactions has substantial negative impacts on quality of life for nearly 10% of adults (1) and 8% of children worldwide that live with food allergies (2), while costing an estimated $25 billion per year in the United States alone (1, 3). Peanut oral immunotherapy (pOIT) is a promising treatment to desensitize participants by increasing tolerance via administration of small, escalating doses of allergen over several months (4, 5) and is shown to be safe and highly successful in suppressing IgE-mediated allergic responses, i.e., desensitization (6). However, pOIT participants can experience gastrointestinal issues and allergic reactions (7), and participants who do not continue ingestion of the allergenic food after pOIT exhibit a high rate of regression to an allergic state (7, 8). Repeated oral administration of allergen during pOIT is shown to decrease peanut-specific IgE and increase peanut-specific IgG4, both of which may play roles in desensitization (9–14). Determining the modulation of IgE and IgG4 epitopes during desensitization might inform safer pOIT regimens, as pOIT has been shown to significantly increase peanut-specific IgG4 based on direct serum antibody measurements and linear peptide binding (13, 15). Here, we utilize linear peptide microarrays and intact allergen ImmunoCAP Immuno Solid-phase Allergen Chip (ISAC) immunoassays to determine IgE and IgG4 epitope shifts over the course of the PTAH PALISADE pOIT phase 3 clinical trial.
2. Materials and methods
2.1. Human sera
Sera were collected from peanut-allergic participants (n = 20) during the Aimmune Therapeutics funded peanut allergy oral immunology clinical studies, ARC001 (16) (NCT01987817) and PALISADE (17) (ARC003, NCT02635776) after informed consent. Participant allergy history was confirmed by clinical history and double-blind placebo controlled oral food challenges (DBPCFC). Phase 3 trials for PTAH (peanut allergen dnfp, PALFORZIA®) were administered with DBPCFC, with the Active group (n = 15) receiving PTAH and the Placebo group (n = 5) receiving oat flour (Avena sativa). Sera were collected prior to pOIT trials and after 12 months of treatment. Participants were removed from trials if they had anaphylactic reactions to peanut. Baseline and exit maximum tolerated dose (MTD) of peanut were administered to participants before and after pOIT and ranged from 3 to 1,000 mg, and all participants undergoing the trial passed this oral food challenge. Baseline and exit skin prick tests (SPTs) were administered to gauge IgE reactivity via mean wheal diameter (MWD) measurements at pre and post timepoints (Table 1).
2.2. Peptide microarrays
Microarrays were utilized for detection of participant sera IgE and IgG4 binding to linear peanut peptides. Synthetic 15-mer linear peptides offset by 5 amino acids that represented the entire amino acid sequences of Ara h 1, 2, 3, 5, 6, 7, 8, 9, 10, and 11 (Supplementary Table 1) were commercially synthesized and spotted in triplicate onto microarray slides by JPT Peptide Technologies (Berlin, Germany). Slide preparation was performed as previously described (18). Slide hybridization was performed on an HS400 Pro Hybridization Station (Tecan, San Jose, CA, USA), with each slide placed in an individual chamber. Slides were blocked in 200 μl filtered SuperBlock (Thermo-Fisher, Waltham, MA, USA) for 30 min at room temperature (RT) under agitation, then washed for 2 min with Tris-buffered saline containing polysorbate-20 (100 mM Tris, 274 mM NaCl, 5.4 mM KCl, and 0.5% Tween-20, TBST). Participant sera were injected into slide chambers and incubated at 4°C overnight (∼16 h) with agitation. Slides were washed as described above before injecting with mouse anti-human IgE (Life Technologies, Grand Island, NY) diluted in SuperBlock (3.3 μg/ml), incubated for 30 min at RT, washed and incubated with diluted Cy3-conjugated goat anti-mouse IgG (0.4 μg/ml; Life Technologies) for 30 min at RT. The wash step was repeated, and slides were dried with N2 gas. Washed and dried slides were scanned with a GenePix-4000B scanner (Software: GenePix Pro 7; Molecular Devices, San Jose, CA, USA) to detect IgE binding measured by Cy3-green (532 nm). Slides were then washed and incubated with 160 ul of anti-human IgG4 (rabbit), incubated for 30 min at RT, washed, and incubated with anti-rabbit Alexaflour-red for 30 min. Slides were scanned for IgG4 signals at 635 nm with a GenePix-4000B scanner.
2.3. Microarray data processing, normalization and epitope identification
GenePix GPR output files containing IgE and IgG4 signal intensities were processed with the Limma package v3.52.4 using R v4.2.1 and RStudio v2022.07.0. Median feature pixel intensities at 532 nm (IgE) and 635 nm (IgG4) were background corrected using the backgroundCorrect.matrix function in Limma (offset = 10, method = “normexp”) and were base 2 logarithm transformed, and the median for each background-corrected peptide triplicate was taken. Cyclic-Loess normalization for each peptide was performed between all possible pairs of arrays for 3 iterations, with a loess smoothing window span of 0.3 using the normalizeBetweenArrays function (method = “cyclicloess”, cyclic.method = “pairs”).
Median peptide fluorescence signal-to-noise ratios (SNRs) across microarray chip lots were used to determine major reactive IgE and IgG4 epitopes. Positive IgE or IgG4 binding was identified as a median SNR >3, which is based on previous chicken sera controls having SNR values of 3. Peptides with positive IgE or IgG4 binding by at least 50% of the sera (i.e., 50% of participant sera had median SNRs >3) were considered major epitopes (19). Peptides were mapped to 3D protein structures using the Protein Data Bank 3D viewer (20).
2.4. Within and between-group comparisons of pre- and post-pOIT peptide intensities
To identify shifts in peptide binding intensity from pre- to post-pOIT within the Active and Placebo treatment groups, paired t-tests of IgE/IgG4 binding intensity ratios (IgE/IgG4) were performed in R v4.2.1 [t-test(paired = TRUE), stats package v4.2.1] for each 15-mer linear peptide from pre- to post-pOIT. Prior to conducting the paired t-tests, the approximate normality of peptide binding intensity distributions were assessed with quantile-quantile plots (Supplementary Figures S1–S10) generated with the geom_qq() and geom_qq_line() functions in R with ggplot v3.4.2. False discovery rate control of paired t-test p values was performed via Benjamini-Hochberg (FDR) correction (21), adjusting for multiple comparisons with a protein (i.e., peptide-level t-tests were penalized based on the number of total peptides within their respective protein) and imposing an adjusted p value threshold <0.05. Peptides with significant adjusted p values were considered to have undergone a significant IgE/IgG4 ratio increase or decrease over the course of the study if their mean difference was >1 or <1, respectively. Cohen's d effect size for paired samples (22, 23) was calculated for each peptide (pre vs. post within treatment groups) and corrected to the Hedges’ g for paired samples (24).
Between-group (i.e., Active vs. Placebo) comparisons of per-peptide IgE/IgG4 were performed with linear regression models (LM) using the lm() function in R: IgE/IgG4 Post = IgE/IgG4 Pre + Treatment. LM fit metrics are available in Supplementary Table 2. Raw p values were adjusted for multiple testing via FDR correction (21). Contrasting least-squares means were computed using the emmeans v1.8.7 package (25). Baseline-adjusted geometric mean post-OIT IgE/IgG4 values were compared per peptide for Active and Placebo groups to determine whether the Active group finished the study with significantly lower or higher IgE/IgG4 ratio compared to Placebo for that peptide, assuming both groups began the study with the same antibody binding intensity. Known epitopes for Ara h 1, 2 and 3 were obtained from the Structural Database of Allergenic Proteins (SDAP) v2.0 (26). Linear B-cell epitopes for Ara h 5-11 were predicted with the BepiPred 2.0 webserver (27), and predicted epitopes with scores >5 were considered in downstream analysis.
2.5. Correlation networks
Correlation matrices were constructed from normalized log2 IgE and IgG4 intensities for pre- and post-pOIT, with pairwise Pearson correlation performed between pre- and post-pOIT matrices via the correlate() function in the corrr R package v0.4.4 (28) (parameters: use = “pairwise.complete.obs”, method = “pearson”). Correlation networks were visualized using the visNetwork v2.1.2 package in R.
2.6. Supervised binary classification
Supervised binary classification models were used to predict desensitization to peanuts and determine which peanut allergen peptide features were most important for those models. IgE/IgG4 of peanut-allergic participant (n = 42) sera from a separate peanut sublingual oral immunotherapy (pSLIT) study utilizing identical microarray-bound peptides and a pre/post experimental design with comparable treatment durations (29) were used for train and test sets (80% train, 20% test), and Active group PALISADE phase 2 (n = 2) and phase 3 (n = 15) participant pOIT IgE/IgG4 were used as the validation (unseen) data set. The pSLIT and PALISADE GPR outputs were processed together with the Limma package v3.52.4 using R v4.2.1, as described in Methods section 2.3. Modified z-scores of per-peptide IgE/IgG4 (features) and pre/post timepoints (labels) were used for supervised binary classifiers using the Python machine learning libraries PyCaret v2.3.10 (30) and SciKit-Learn v0.23.2 (31). Feature selection was based on the importance of weights (Scikit-learn SelectFromModel using LightGBM estimator) using mean importance value thresholds of 0.6, 0.7, and 0.8. For each importance value threshold, 40 sessions comparing classifier performance were run in PyCaret using different random seeds states ranging from 0 to 1,000. Within each session (n = 120), train and test sets were sampled from the pSLIT dataset (80% train, 20% test), and all base classification estimators available in the SciKit-Learn v0.23.2 library were trained and evaluated with Stratified KFold cross validation (splits = 10) on the training set. The top 4 models were selected based on accuracy, and hyperparameter tuning was performed for these models with a Random Grid Search with 50 iterations using Stratified KFold cross validation (splits = 10) with F1 score as the optimization metric. The original models were used in downstream analysis if they outperformed their tuned versions. Ensemble methods were performed via bagging (32), stacking (33), and boosting (34), using 15 base estimators, with F1 score as the optimization metric. For all the generated models, the base estimator or ensemble estimator with the highest F1 score was returned, finalized with the entire train and test datasets, and predictions were performed on the unseen PALISADE dataset. The predictive performance on the unseen dataset for each finalized model was compared and the best model was selected based on F1 score and Matthews correlation coefficient (MCC) (35). A tree-based ensemble classifier [Extremely Randomized Trees (ERT) classifier (36)] was selected [parameters: class_weight='balanced', criterion = 'entropy', max_depth = 2, max_features = 'sqrt', min_impurity_decrease = 0.0005, min_samples_leaf = 5, min_samples_split = 7, n_estimators = 200, n_jobs = −1(36)] was selected as the best performing model (FST = 0.6, random_state = 388). Gini Impurity (37) was used to identify peptide features with the largest effect on this particular ERT model's prediction of desensitization.
2.7. ImmunoCAP ISAC immunoassays
ImmunoCAP Immuno Solid-phase Allergen Chip (ISAC) immunoassays (Thermo-Fisher, Upsala, Sweden) were used to measure pre- and post-pOIT participant sera IgE binding to purified intact Ara h 1, 2, 3, 6, 8, 9 proteins. IgE signal intensity levels in ISAC Standardized Units (ISU-E) with an operating range of 0.3–100 ISU-E were generated by Phadia Microarray Image Analysis software v1.2 (Thermo-Fisher Scientific, Waltham, MA, USA). IgE intensities were log2 transformed [log2(x + 1)]. Wilcoxon signed rank tests were used to compare IgE binding intensity distributions from pre to post for each protein using the wilcox.test() function in the R stats package v.4.2.1 (parameters: paired = TRUE, alternative = “two.sided”, exact = TRUE, conf.level = 0.95). Spearman correlation matrices were constructed using the cor() function in the R stats package v4.2.1, and correlation plots were made using the corrplot v0.92 package in R.
3. Results
3.1. IgE binding to intact allergens
ISAC immunoassays were utilized for IgE binding to whole proteins to determine conformational aspects to binding. PALISADE participant IgE binding was shown via Wilcoxon signed rank tests to have significantly different distributions from pre-to-post pOIT for Ara h 2 (p = 0.018), Ara h 3 (p = 0.0087), and Ara h 6 (p = 0.00073), with IgE binding exhibiting an overall decrease in the study population after pOIT (Figure 1A). The Ara h 1 IgE binding distributions of the Active participant group were not significantly different from pre to post pOIT (p = 0.15). The total IgE binding intensity to Ara h 8 and Ara h 9 was much lower compared to Ara h 1, 2, 3, and 6. Ara h 8 binding increased for two and decreased for two participants, however all other participants did not exhibit notable IgE binding shifts to this protein and pre-to-post distributions were not significant (p = 1) (Figures 1A,B). Similarly, only one participant showed increased binding to Ara h 9 from pre-to-post, while all other participants had no notable binding and pre-to-post distributions were not significantly different (p = 0.86; Figures 1A,B). Ara h 8 IgE binding was not strongly correlated with any other intact protein at pre and post-pOIT timepoints (Figures 1C,D). On a per-patient basis, those with higher IgE binding to Ara h 2 exhibit higher binding to the homologous 2S albumin Ara h 6 (Figure 1B), with a population-level Spearman correlation coefficient of 0.95 at the pre-pOIT timepoint (Figure 1C) and 0.9 at post-pOIT (Figure 1D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279290/falgy-04-1279290-HTML-r1/image_m/falgy-04-1279290-g001.jpg)
Figure 1. ImmunoCAP ISAC immunoassays indicate significantly reduced IgE binding to intact Ara h 2, 3, and 6 post pOIT for active group participants. (A) ISAC Pre vs. post pOIT for participant population. P values are for paired Wilcoxon signed rank tests, with ** indicating p < 0.05 and *** for p < 0.01. Arrows between dots indicate pre-to-post log2 IgE intensity shifts for individual participants. Box plot notches indicate the 95% confidence interval of the median. (B) Log2 transformed IgE intensity values per protein for Active group participants, pre and post pOIT. (C) Spearman correlations of IgE intensity values for Active group participants, pre pOIT. (D) Spearman correlations of IgE intensity values for Active group participants, post pOIT.
3.2. IgE and IgG4 binding to linear peanut allergen peptides
Pre to post-pOIT linear epitope shifts were identified using IgE/IgG4. In general, significant shifts in antibody binding occurred within regions of overlapping peptides. Based on paired t-tests, Ara h 1 had 1.6% of peptides exhibiting significantly increased IgE/IgG4 (increasing IgE and/or decreasing IgG4) from pre- to post-pOIT in the Active group, while 0.8% of peptides had significant shifts towards decreased IgE/IgG4 (indicating increasing IgG4 and/or decreasing IgE binding) (Figure 2A, Table 2). Known epitope regions containing significantly different peptides include those containing peptides 5–7 (IgE/IgG4 decrease), peptides 15–22 (IgE/IgG4 decrease), peptides 28–31 (IgE/IgG4 decrease), and peptide 116 (IgE/IgG4 increase). In the Placebo group, 8.87% of Ara h 1 peptides underwent a significant IgE/IgG4 decrease (Table 2, Figure 2B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279290/falgy-04-1279290-HTML-r1/image_m/falgy-04-1279290-g002.jpg)
Figure 2. Pre-to-post pOIT IgE/IgG4 indicate increased shifts from IgE to IgG4 binding in active vs. Placebo group. (A) Active group, Ara h 1, 2, 3, 5, 6. (B) Placebo group, Ara h 1, 2, 3, 5, 6. (C) Active group, Ara h 7, 8, 9, 10, 11. (D) Placebo group, Ara h 7, 8, 9, 10, 11. Orange and blue dots respectively indicate peptides with a significant increase or decrease in IgE/IgG4 ratio from pre-to-post pOIT based on within-group paired t-tests, as determined by an FDR-adjusted p value <0.05. Blue and orange number labels respectively indicate peptides with a significant decrease or increase in pre-to-post pOIT IgE/IgG4 ratio for Active compared to Placebo based on FDR-adjusted p values <0.15 for an LM. Horizontal lines at y = 0.2, 0.5, and 0.8 respectively indicate boundaries for small, medium, and large effect size benchmarks (22). Points with “^” indicate peptide sequences with regions matching previously identified IgE epitopes listed in the SDAP database, while black circles around points indicate linear peptides predicted with BepiPred v2.0.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279290/falgy-04-1279290-HTML-r1/image_m/falgy-04-1279290-t002.jpg)
Table 2. Percentage of peptides that had significant IgE/IgG4 ratio shifts from pre- to post-pOIT based on study population paired t-tests.
For the major peanut allergen Ara h 2 (Ara h 2.0201), 21.2% of the 15-mer peptides exhibited increased IgE/IgG4% and 3% exhibited decreased IgE/IgG4 over the 12-month pOIT period within the Active group based on paired t-tests, while no significant changes were observed in the Placebo group (Table 2). Peptide regions 6–8 and 28–30 form known IgE epitopes comprising α-helices that span protein segments I, III, and IV (38) (39), and these displayed decreased IgE/IgG4 within the Active group based on paired t-tests (Figure 2A). IgE/IgG4 decreased for peptide 31 in the Active group over the course of pOIT with moderate effect size (Figure 2A) and may be a potential IgE and IgG4 epitope based on median SNRs from part of the population (Supplementary Figure 12). This may indicate a previously undescribed linear epitope in segment IV of Ara h 2 that has undergone a shift from IgE to IgG4 binding in the Active group over the course of pOIT. Linear regression analysis did not identify any peptides with significant within-group IgE/IgG4 shifts where the Active group finished the study with significantly higher or lower IgE/IgG4 vs. Placebo group. Using an adjusted p value cutoff <0.05, only Ara h 2 peptide 17 was identified as being significantly different (Active group finished with significantly lower IgE/IgG4). Using a more liberal adjusted p value cutoff of <0.15 to adjust for small library size, known IgE epitope regions (38) spanning Ara h 2 peptides 13 (p = 0.1), 14 (p = 0.1), 16 (p = 0.059), 17 (p = 0.046), and 27 (p = 0.1) were shown to have lower IgE/IgG4 in the Active group vs. Placebo group at the end of pOIT.
The major peanut allergen Ara h 3 (11S storage protein Ara h 3.0101) had 11% of peptides undergo IgE/IgG4 increases and 13% undergo IgE/IgG4 decreases in the Active group, with no significant changes in the Placebo group (Table 2). Within the N-terminal domain, binding at the peptides 1–5 region suggested an interplay of IgE and IgG4, as IgE/IgG4 decreased within peptide region 1–3 while the neighboring peptide 4–5 region that includes known IgE Epitope 1 (40) had increased IgE/IgG4. Significant increases in IgE/IgG4 were observed in peptides 49 and 50, which overlap IgE Epitope 2 (GNIFSGFTPEFLEQA) (40) with sequence EFLEQA. IgE Epitope 3 (VTVRGGLRILSPDRK) (40) underwent a decrease in IgE/IgG4 at peptides 55 and 56, indicating that this IgE epitope shifted to IgG4 binding over the course of pOIT. At peptide 61, which includes the IgE Epitope 4 sequence (DEDEYEYDEEDRRRG), IgE/IgG4 increased and median SNR was >3 for both IgE and IgG4. This suggests that the flexible region of the N-terminal domain containing Epitopes 3 and 4 exhibited flanking IgG4 and IgE epitopes post-pOIT. The C-terminal domain displayed IgE/IgG4 decreases in peptides 95–97 and median SNR >3 for IgE and IgG4 in peptide 99, all of which fall within a known epitope on the flexible loop region on the periphery of the two adjacent trimers that form Ara h 3 (41).
Peanut profilin Ara h 5 is an allergen linked to pollen-associated peanut allergy (42). No peptides were significantly different in the Placebo group, and peptides 17 and 19 were significantly different in the Active group. Peptides 17 and 19 comprise a beta sheet and coil, and respectively exhibited increased IgG4 and IgE from pre to post pOIT with high effect sizes for both. This may indicate that the Ara h 5 peptide 17–19 region is an IgE epitope that experiences competition with IgG4 binding following pOIT.
Ara h 6, a 2S albumin, is a potent peanut allergen with high similarity to Ara h 2. No significant peptide shifts were identified in the Placebo group (Figure 2B). Ara h 6 peptide 5 is a linear residue on the N-terminus nearly identical to Ara h 2 (39) and was shown to have significantly lower IgE/IgG4 from pre to post-pOIT in the Active group based on paired t-tests (Figure 2A). Peptides 23 and 25 also had significantly higher IgG4 in the Active group from pre to post based on paired t-tests (Figure 2A), and these peptides constitute helix-loop-helix residue that is near-identical with Ara h 2 and forms the transition from the 4th to 5th α-helix (39).
Ara h 7 did not exhibit significant shifts in either the Active or Placebo group based on paired t-tests (Figures 2C,D), however median SNRs indicate potential epitopes. Peptides 30–32 forming the C-terminus had high effect sizes as well as IgG4 SNRs that may indicate these peptides to be IgG4 epitopes. Ara h 8, which is known as an inhaled sensitizing protein (43), is linked to oral allergy syndrome and exhibits cross-reactivity with the Bet v 1 family (44). Ara h 8 was shown to have the peptide with the highest effect size in the Active group (peptide 16), and this peptide had a significant IgE/IgG4 increase in the Active group. While median IgE and IgG4 SNRs do not suggest that peptide 16 is an epitope based on these metrics, the overlapping peptide 17 has median IgG4 SNR values indicating an epitope region (Supplementary Figure 17).
The lipid transfer protein Ara h 9 is implicated in oral allergy syndrome and is indicated as a major peanut allergen in Mediterranean populations (45, 46). Ara h 9 had 23.5% of its peptides undergo increases in IgG4 binding within the Active group (Figure 2C, Table 2), and no significant peptide shifts were identified in the Placebo group (Figure 2D). Peptides 6–7, 9, and 16 had decreased IgE/IgG4 in the Active group, while peptide 13 had increased IgE/IgG4. Peptide 16 was predicted to be a linear B-cell epitope by BepiPred 2.0 (Figure 2C). Median SNRs for Ara h 9 reveals possible IgG4 epitopes within the peptide 2–6 region, as well as peptides 9, 10, 14 and 15, all of which have been previously identified as IgE epitopes (18) and overlap with significantly-different peptides (Figure 2C, Supplementary Table 1) (18).
The oleosins Ara h 10 and 11 exhibited no significant IgE/IgG4 shifts in the Placebo group. Active group Ara h 10 peptide 12 had significantly higher IgE/IgG4 from pre to post pOIT and a large effect size (Figure 2B) decreased IgE/IgG4. Ara h 11 had significantly higher IgE/IgG4 on its N-terminus [peptide 1, (M)AEALYY], and peptides 16–20, 22, and 23 suggest an IgE binding region that is flanked by an IgG4 binding region (Figure 2B). Two regions in Ara h 11 comprised of peptides 5–6 and 16–18 exhibited significant IgE/IgG4 shifts. Peptides 5–6 respectively showed increased and decreased IgE/IgG4 on either side of this overlapping peptide, while IgE/IgG4 increased in the peptide 16–18 region.
3.3. Peptide pre-vs.-post pOIT correlations and supervised binary classification
Given the assumption that pre- and post-pOIT sera epitope binding patterns represent a snapshot of a dynamic system, we wished to identify relationships of immunoglobulin binding intensity among individual allergen epitopes that occurred over the course of the study. Pearson correlation networks were constructed to highlight IgE and IgG4 binding relationships of Ara h peptides exhibiting significant IgE/IgG4 shifts over the course of pOIT. The IgE binding intensity for major Ara h 2 epitopes were shown to be positively correlated with known and newly identified epitopes of other peanut proteins (Figure 3A). IgE binding for Ara h 2 peptide region 28–29 was shown to be positively correlated with Ara h 6 peptide 23 (Figure 3A), all of which shifted towards decreased IgE/IgG4 post pOIT and exhibited high Hedges' g effect sizes (>0.8) (Figure 2A). The Ara h 2 peptide 28–29 region exhibited weaker positive IgE correlations with Ara h 3 peptides 55, 56 (Ara h 3 Epitope 3) and 97 [C-terminal region flexible loop epitope (41)]. The known Ara h 2 epitope peptide 7 had a positive IgE correlation with Ara h 3 peptide 97. Ara h 2 peptides 8 and 9 have positive IgE correlation coefficients with Ara h 9 peptide 2 (Figure 3A). Ara h 2 peptide 9 is a known epitope, and Ara h 9 peptide 2 is a candidate major IgE and IgG4 epitope based on median SNR (Supplementary Figure 18). Ara h 3 had strong positive IgE and IgG4 correlations with multiple Ara h 2 IgE epitopes and significantly different Ara h 6 peptides (Figures 3A,B). Ara h 5 peptide 19 was shown to have positive IgG4 correlation coefficients with Ara h 3, Ara h 6, and Ara h 9 peptides that exhibited decreases in IgE/IgG4 (Figure 3B). This may indicate that this Ara h 5 IgE epitope undergoes antibody binding shifts in concert with these other peptides, and that it may be an important peptide for monitoring tolerance due to its high effect size (Figure 2A). Ara h 10 peptide 11 was found to have positive IgE correlations with Ara h 3, Ara h 5, and Ara h 9 peptides that displayed increased IgE/IgG4 at the population level, along with Ara h 6 and Ara h 2 peptides whose IgE/IgG4 decreased (Figure 3A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279290/falgy-04-1279290-HTML-r1/image_m/falgy-04-1279290-g003.jpg)
Figure 3. Correlation network of Ara h peptides for active group participants. (A) IgE and (B) IgG4 Pearson correlations were performed for peptides found to have significantly different within-treatment IgE/IgG4 from pre-to-post pOIT based on paired t-tests. Orange edges indicate Pearson correlation coefficients >0.7, and purple edges indicate coefficients <-0.7. The width of the edges indicates a greater absolute value of the correlation coefficient. Square nodes indicate a peptide with significantly decreasing IgE/IgG4 ratio (increasing IgG4), and triangular nodes indicate increasing IgE/IgG4 ratio (increasing IgE) based on paired t-test results.
We trained an ERT classifier to predict participant desensitization to peanut (Table 3) and used Gini Impurity feature importance to identify peptides that were most important to this particular model (Figure 4). Peptides from Ara h 1 (peptides 116, 22), Ara h 2 (peptides 16, 14, 17, 20, 13), Ara h 3 (peptides 67) and Ara h 5 (peptides 2, 14)), comprised the top 10 most important features for the ERT classifier (Figure 4). Ara h 2 peptides 16 and 14 were the two most important features for the finalized model, followed by predicted IgE epitope Ara h peptide 2 and known IgE epitopes Ara h 2 peptide 17, Ara h 1 peptide 116 and Ara h 1 peptide 22 (Figure 4). Ara h 2 peptide 20, Ara h 5 peptide 14, Ara h 3 peptide 67, and known IgE epitope Ara h 2 peptide 13 comprised the next tier of important features (Figure 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279290/falgy-04-1279290-HTML-r1/image_m/falgy-04-1279290-t003.jpg)
Table 3. Performance metrics for top 10 finalized binary classifier models trained on pSLIT and predicting desensitization on unseen PALISADE dataset.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1279290/falgy-04-1279290-HTML-r1/image_m/falgy-04-1279290-g004.jpg)
Figure 4. Important peptide features for the extremely randomized trees binary classifier in predicting desensitization. The peptide features that contribute most to the scoring of this particular pSLIT-trained ERT classifier to predict desensitization were identified with Gini Impurity feature importance and are shown from top to bottom on the y-axis in the order of decreasing importance.
4. Discussion
This study utilized sera from peanut-allergic individuals undergoing pOIT with PTAH to identify linear peptide-based IgE and IgG4 epitope shifts of Ara h 1, 2, 3, 5, 6, 7, 8, 9, 10 and 11, as well as shifts in IgE binding to intact Ara h 1, 2, 3, 6, 8, and 9 allergens. A larger set of peanut proteins and peptides were assessed in this study compared to most previous pOIT studies, and specific epitope shifts in participants undergoing pOIT using a standardized, FDA-approved treatment were highlighted. In comparison to the Placebo treated, the Active group population exhibited a greater number of significant shifts in IgE epitopes binding to IgG4 after 12 months of pOIT, as seen in previous pOIT studies (6, 14, 47). In 9 of the 15 (60%) Active group participants, the peanut-specific (ps) IgE and total IgE was increased or relatively unchanged after 12 months of pOIT, while the MWD following SPT decreased in 14 of 15 (93.3%-100%) participants (exit MWD was not measured for participant 004-007) (Table 1). However, 11 (73.3%) of the Active group participants exhibited relative increases in psIgG4 ranging from 33.3% to 9,244% following 12 months of pOIT. This indicates that at the 12-month point, the total and ps-IgE for a majority of the subjects has not started to decrease in response to treatment, and that the shifts in IgE/IgG4 of the peptides are more likely due to the increase in IgG4 than decrease in IgE. Therefore, it appears that the rise in IgG4 levels during pOIT may be more important than the decrease in psIgE to modulating desensitization, as evident by maximum tolerated dose (MTD, Table 1) at 12 months. The post-pOIT decrease in participant IgE binding to intact proteins in ISAC immunoassays may be due in part to inhibition of IgE binding caused by the increased presence of competing IgG4 after immunotherapy (47–49), which may better reflect the in vivo state of a participant's allergen sensitivity (48, 50).
In addition to overall IgE and IgG4 levels in patients, this study offers insight into antibody binding shifts for individual proteins and peptides that occur over the course of pOIT. At a population level, the IgE and IgG4 epitopes overlapped and the IgE epitopes showed limited changes with 12 months of pOIT. The stability of IgE epitopes through allergen-specific immunotherapy or vaccination has previously been shown at an individual level in the case of inhaled allergens (51). Another study using the same grass pollen allergen microarrays found that allergen immunotherapy induced expanded IgG-type and IgE antibodies against novel epitopes after one year of treatment compared to untreated controls (52). Ara h 1, 2, 3, and 6 are currently accepted as the four major peanut allergens (53), however it was surprising that the PALISADE participant group did not exhibit significant population-level IgE binding changes to intact Ara h 1. It appears that ISAC Ara h 1 IgE intensity increases in some participants (e.g., participants 062-014 and 035-013) or does not shift at all (participant 062-017) (Figure 1B). Interestingly, participants 062-014 and 035-013 both had a high psIgE/total IgE ratio and were able to reach a MTD of 600 mgs at the 12-month exit timepoint (Table 1). Ara h 2 IgE binding intensity increased after pOIT in these two participants, and their sera IgE did not bind to Ara h 3 (Figure 1B). Meanwhile, 25 of 124 Ara h 1 peptides showed a decrease in IgE/IgG4 with only 3 peptides having increased IgE/IgG4 with high effect size (Figure 2A). It is possible that 2 of these Ara h 1 peptides (107 and 116) are important in conformational epitopes that drive the increased IgE binding to intact Ara h 1 on the ISAC array for certain participants. Ara h 6 shares 60% sequence identity and multiple IgE epitopes with Ara h 2, resulting in a degree of cross-reactivity and the common occurrence of peanut-allergic individuals producing IgE against Ara h 2 and 6 (54). This likely explains why the Active group had significant reductions and comparable intensity values for ISAC IgE binding to Ara h 2 and Ara h 6 from pre- to post-pOIT. In addition, the Ara h 2 IgE epitope spanning peptides 28 and 29 was shown to be an IgE and IgG4 epitope based on median SNR (Supplementary Figure 12) and has positive IgE correlations with Ara h 6 peptide 23 (Figure 3A), which was identified as an IgE and IgG4 epitope (Supplementary Figure 15) and contains a portion of a highly homologous epitope to Ara h 2 IgE epitope 7 (KRELRMLP) (19).
Aside from shifts in antibody binding to major peanut allergens, this study also highlights shifts in binding to minor allergens Ara h 5, 7, 8, 9, 10, and 11 over the course of pOIT. Ara h 5 peptide 19 was shown to have positive IgG4 correlation coefficients >0.7 with Ara h 3, Ara h 6, and Ara h 9 peptides (Figure 3B). Interestingly, the sensitizing protein Ara h 5 was found to contribute important features (peptides 2 and 14) for the ERT classifier predicting peanut desensitization (Figure 4). These Ara h 5 peptides were not found to be significantly different based on paired t-tests, however they were indicated to be IgE epitopes based on median SNR (Supplementary Figure 14). Given the potential role of Ara h 5 as a sensitizing protein for peanut allergy due to its similarity to cross-reactive proteins such as Bet v 2 (birch pollen) and Hev b 8 (latex profilin) (55), gaining a better understanding of Ara h 5's relationship to other peanut allergens with regards to pOIT may indicate its use as a marker for tolerance. Ara h 7 did not have significantly different IgE or IgG4 binding in Active participant sera following pOIT, as opposed to a previous pOIT study (PRROTECT) where participants received omalizumab (56). It is possible that the concurrent 19-week omalizumab treatment with pOIT may have resulted in increased IgG binding to Ara h 7 in PRROTECT sera (56), vs. non-significant changes in binding in the PALISADE study. However, median SNRs for Ara h 7 peptides 4, 15, 21, and 26 identify them as IgE and IgG4 epitopes in this study. This study also explores IgE and IgG4 binding shifts in the oleosins Ara h 10 and 11, which have not been studied with regards to pOIT. This study reveals the presence of IgE and IgG4 binding regions with high effect sizes, indicating that they are not only statistically significant but the magnitude of this binding difference over the course of pOIT lends to the potential importance of these regions in applied peanut allergy diagnostics. Further studies are needed to determine the roles of these oleosin binding regions in peanut allergy and tolerance.
The notion of major allergen epitope regions shifting antibody binding in concert with sensitizing proteins is supported by Ara h 9 peptide links to major allergens. Positive correlations >0.7 were seen for of Ara h 1 peptide 116 with Ara h 9 peptides 7 and 9 (IgE, Figure 3A), Ara h 2 peptide 29 with Ara h 9 peptide 16, and Ara h 6 peptide 5 with Ara h 9 peptide 9 (IgG4, Figure 3B). The positive IgE correlations of Ara h 9 peptides with known major allergen IgE epitopes and epitopes important in predicting desensitization (Ara h 1 peptide 116, Figure 4), along with these peptides experiencing a shift towards increased IgG4 binding, highlights that Ara h 9 may have greater roles in the development and treatment of allergy given that it is a major allergen in Mediterranean populations (18, 45, 46). Ara h 9 peptide 2 was a major epitope for the PALISADE participant pool based on median SNRs (Supplementary Figure 18), and previous studies have shown US participant sera to be IgE reactive to this peptide as well as other previously described binding regions (18). The US participants described in this study did not exhibit IgE binding to the intact Ara h 9 proteins, in line with ISAC results of US patients described in Kronfel et al. (18). Further research is required to determine whether binding to linear Ara h 9 epitopes influences sensitivity to major allergens in US patients.
The oleosins Ara h 10 and Ara h 11 are understudied with regards to peanut allergy, and this study provides early insight into potential roles that these proteins may play in allergy and development of tolerance to peanuts. Besides identifying epitopes and potential biomarkers for pOIT, correlation networks show a first glance into the links these proteins may have with better-studied allergens. For instance, Ara h 10 and Ara h 11 have positive IgE and IgG4 correlation coefficients with peptides in Ara h 1, 2, 3, 5, 6, 8, and 9 (Figure 3). The presence of peptides in Ara h 10 and Ara h 11 with significantly different IgE/IgG4 shifts and high effect sizes (Figure 2C) suggest that these could be important regions in understudied proteins that have applications towards allergy diagnosis and determination of tolerance after pOIT.
This study utilized an ERT binary classifier to predict participant desensitization to peanut and we used Gini Impurity feature importance to identify the most important peptide features for this particular model's classification task (Figure 4). This method identified peptide features that were significantly different based on between-group (LM) and within-group (paired t-test) IgE/IgG4 comparisons (Figure 2, Figure 4). Ara h 2 peptides 13, 14, 16, and 17 were identified as top 10 important features for the ERT model while being significantly different in the Active vs. the Placebo group after the course of pOIT based on an LM. While Ara h 2 peptides 13, 14 and peptide 16 should be interpreted with some caution given their higher adjusted p values (p = 0.1, 0.1, and 0.059, respectively) versus Ara h 2 peptide 17 (p = 0.046), this does support the finding that these peptides were important for predicting desensitization since they showed significantly lower IgE/IgG4 (increased IgG4 and/or decreased IgE) in the Active treatment group, which underwent pOIT and became desensitized, vs. the Placebo group, which did not. While machine learning algorithms can be very powerful tools, the limited number of participant sera obtained for the PALISADE study and necessity of training the model on a larger pSLIT dataset may be responsible for the identification and some inconsistencies in important peptides compared to paired t-tests. Future studies will require larger amounts of training data to increase model performance.
Increasing our understanding of which epitopes shift from IgE to IgG4 binding over the course of oral immunotherapy is vital towards developing a framework for identifying successful pOIT treatments, as well as improving potential peptide-based diagnostic methods for food allergy. It is also important to note that the binding of IgE and IgG4 to extracts, components, or allergen peptides can diverge and that each level of binding may contribute important information regarding the state of immune modulation and subject response to types and duration of treatments. Future studies with more timepoints and clinical data regarding other allergies will aid in determining long-term shifts in IgE and IgG4 epitope binding once participants reach a maintenance phase. This study contributes to our growing knowledge of peanut allergy dynamics over the course of oral immunotherapy, demonstrates which epitopes shift for participants taking PTAH, and provides an initial framework for identifying critical epitopes that determine if peanut-allergic participants have reached desensitization after treatment.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, and further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving human participants were reviewed and approved by Tulane School of Medicine and Aimmune Nestlé Health Science (NCT01987817 and NCT02635776). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin.
Author contributions
IMR organized, reformatted and analysed the data and meta data, wrote code for and performed the analysis, and wrote the manuscript; CMK collected, organized and analysed data and meta data; ARR assisted with code and directed the selection and use of bioinformatics and machine learning software; LTS helped with data collection and analysis; HC optimized and performed microarray experiments; JKM performed ISAC experiments and provided the data; RJS provided US subject sera and related clinical data and funding and helped oversee the project; RR and SAT helped with writing and submission of the manuscript; JBN helped with data collection and analysis; MDK provided data for the machine learning training set. SJM and BKH supervised the project, BKH provided funding for the project, and SJM helped with writing and submission of the manuscript and provided funding for the project. All authors reviewed and approved the manuscript.
Acknowledgments
We are extremely grateful to Thermo-Fisher who provided the ISAC array reader and tremendous support for this project. Thermo-Fisher was not involved in the study design, analysis, interpretation of data, the writing of this article or the decision to submit it for publication. We are also very appreciative of Dr. Yuzhu Zheng and Dr. Heping Cao for their critical review of the manuscript. We greatly appreciate Dr. Quentin Read for evaluating this study's statistical methods. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the USDA. The USDA is an equal opportunity provider and employer.
Conflict of interest
SAT is an employee of Aimmune Nestle Health Science. This study was funded by The United States Department of Agriculture (USDA), Agricultural Research Service and by Aimmune Nestlé Health Science, via a Cooperative Research and Development Agreement (CRADA agreement #58-6054-1-0011). Aimmune Nestlé was involved in performing the OIT study with Palforzia, which is published information, but was not involved in the design and interpretation of this study.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/falgy.2023.1279290/full#supplementary-material
References
1. Gupta RS, Warren CM, Smith BM, Blumenstock JA, Jiang J, Davis MM, et al. The public health impact of parent-reported childhood food allergies in the United States. Pediatrics. (2018) 142(6):e20181235. doi: 10.1542/peds.2018-1235
2. Sicherer SH, Sampson HA. Food allergy: a review and update on epidemiology, pathogenesis, diagnosis, prevention, and management. J Allergy Clin Immunol. (2018) 141(1):41–58. doi: 10.1016/j.jaci.2017.11.003
3. Bunyavanich S, Rifas-Shiman SL, Platts-Mills TA, Workman L, Sordillo JE, Gillman MW, et al. Peanut allergy prevalence among school-age children in a US cohort not selected for any disease. J Allergy Clin Immunol. (2014) 134(3):753–5. doi: 10.1016/j.jaci.2014.05.050
4. Jones SM, Pons L, Roberts JL, Scurlock AM, Perry TT, Kulis M, et al. Clinical efficacy and immune regulation with peanut oral immunotherapy. J Allergy Clin Immunol. (2009) 124(2):292–300.e1-97. doi: 10.1016/j.jaci.2009.05.022
5. Patrawala M, Shih J, Lee G, Vickery B. Peanut oral immunotherapy: a current perspective. Curr Allergy Asthma Rep. (2020) 20(5):14. doi: 10.1007/s11882-020-00908-6
6. Vickery BP, Berglund JP, Burk CM, Fine JP, Kim EH, Kim JI, et al. Early oral immunotherapy in peanut-allergic preschool children is safe and highly effective. J Allergy Clin Immunol. (2017) 139(1):173–81.e8. doi: 10.1016/j.jaci.2016.05.027
7. Chu DK, Wood RA, French S, Fiocchi A, Jordana M, Waserman S, et al. Oral immunotherapy for peanut allergy (PACE): a systematic review and meta-analysis of efficacy and safety. Lancet. (2019) 393(10187):2222–32. doi: 10.1016/S0140-6736(19)30420-9
8. Chinthrajah RS, Purington N, Andorf S, Long A, O'Laughlin KL, Lyu SC, et al. Sustained outcomes in oral immunotherapy for peanut allergy (POISED study): a large, randomised, double-blind, placebo-controlled, phase 2 study. Lancet. (2019) 394(10207):1437–49. doi: 10.1016/S0140-6736(19)31793-3
9. Kim EH, Bird JA, Kulis M, Laubach S, Pons L, Shreffler W, et al. Sublingual immunotherapy for peanut allergy: clinical and immunologic evidence of desensitization. J Allergy Clin Immunol. (2011) 127(3):640–6.e1. doi: 10.1016/j.jaci.2010.12.1083
10. Kim JD, Kim SY, Kwak EJ, Sol IS, Kim MJ, Kim YH, et al. Reduction rate of specific IgE level as a predictor of persistent egg allergy in children. Allergy Asthma Immunol Res. (2019) 11(4):498–507. doi: 10.4168/aair.2019.11.4.498
11. Blumchen K, Trendelenburg V, Ahrens F, Gruebl A, Hamelmann E, Hansen G, et al. Efficacy, safety, and quality of life in a multicenter, randomized, placebo-controlled trial of low-dose peanut oral immunotherapy in children with peanut allergy. J Allergy Clin Immunol Pract. (2019) 7(2):479–91.e10. doi: 10.1016/j.jaip.2018.10.048
12. Orgel K, Burk C, Smeekens J, Suber J, Hardy L, Guo R, et al. Blocking antibodies induced by peanut oral and sublingual immunotherapy suppress basophil activation and are associated with sustained unresponsiveness. Clin Exp Allergy. (2019) 49(4):461–70. doi: 10.1111/cea.13305
13. Uotila R, Kukkonen AK, Greco D, Pelkonen AS, Mäkelä MJ. Peanut oral immunotherapy increases IgG4 to ara h 1, 2, and 6 but does not affect IgG4 to other allergens. Pediatr Allergy Immunol. (2019) 30(2):248–52. doi: 10.1111/pai.13012
14. Vickery BP, Lin J, Kulis M, Fu Z, Steele PH, Jones SM, et al. Peanut oral immunotherapy modifies IgE and IgG4 responses to major peanut allergens. J Allergy Clin Immunol. (2013) 131(1):128–34.e1-3. doi: 10.1016/j.jaci.2012.10.048
15. Chu DK, Baumert JL, Taylor SL, Nordlee JA, Nham T, Bramson J, et al. Peanut allergen reaction thresholds during controlled food challenges in 2 Canadian randomized studies (Canada-ARM1 and PISCES). J Allergy Clin Immunol Pract. (2021) 9(6):2524–6.e2. doi: 10.1016/j.jaip.2021.02.009
16. Bird JA, Spergel JM, Jones SM, Rachid R, Assa'ad AH, Wang J, et al. Efficacy and safety of AR101 in oral immunotherapy for peanut allergy: results of ARC001, a randomized, double-blind, placebo-controlled phase 2 clinical trial. J Allergy Clin Immunol Pract. (2018) 6(2):476–85.e3. doi: 10.1016/j.jaip.2017.09.016
17. Vickery BP, Vereda A, Casale TB, Beyer K, Du Toit G, Hourihane JO, et al. AR101 oral immunotherapy for peanut allergy. N Engl J Med. (2018) 379(21):1991–2001. doi: 10.1056/NEJMoa1812856
18. Kronfel CM, Cheng H, McBride JK, Nesbit JB, Krouse R, Burns P, et al. Ige epitopes of ara h 9, jug r 3, and pru p 3 in peanut-allergic individuals from Spain and the US. Front Allergy. (2022) 3:1090114. doi: 10.3389/falgy.2022.1090114
19. Otsu K, Guo R, Dreskin SC. Epitope analysis of ara h 2 and ara h 6: characteristic patterns of IgE-binding fingerprints among individuals with similar clinical histories. Clin Exp Allergy. (2015) 45(2):471–84. doi: 10.1111/cea.12407
20. Berman HM, Battistuz T, Bhat TN, Bluhm WF, Bourne PE, Burkhardt K, et al. The protein data bank. Acta Crystallogr. (2002) D58(Pt 6 No 1):899–907. doi: 10.1107/s0907444902003451
21. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. (1995) 57(1):289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
22. Cohen J. Statistical power analysis for the behavioral sciences. 2 ed. Hillside, NJ: Lawrence Erlbaum Associates (1988).
24. Hedges LVO, Olkin I. Statistical Methods for Meta-Analysis. New York, NY: Academic Press (1985).
26. Ivanciuc O, Schein CH, Braun W. SDAP: database and computational tools for allergenic proteins. Nucleic Acids Res. (2003) 31(1):359–62. doi: 10.1093/nar/gkg010
27. Jespersen M, Peters B, Nielsen M, Marcatili P. BepiPred-2.0: improving sequence-based B-cell epitope prediction using conformational epitopes. Nucleic Acids Res. (2017) 3:W24–9. doi: 10.1093/nar/gkx346
29. Kim EH, Keet CA, Virkud YV, Chin S, Ye P, Penumarti A, et al. Open-label study of the efficacy, safety, and durability of peanut sublingual immunotherapy in peanut-allergic children. J Allergy Clin Immunol. (2023) 151(6):1558–65.e6. doi: 10.1016/j.jaci.2023.01.036
31. Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, et al. Scikit-learn: machine learning in python. J Mach Learn Res (2011) 12:2825–30.
33. Wolpert DH. Stacked generalization. Neural Netw. (1992) 5(2):241–59. doi: 10.1016/S0893-6080(05)80023-1
34. Friedman JH. Greedy function approximation: a gradient boosting machine. Ann Stat. (2001) 29(5):1189–232. doi: 10.1214/aos/1013203451
35. Chicco DJ G. The advantages of the matthews correlation coefficient (MCC) over F1 score and accuracy in binary classification evaluation. BMC Genomics. (2020) 21(6):1–13. doi: 10.1186/s12864-019-6413-7
36. Geurts PE D, Wehenkel L. Extremely randomized trees. Mach Learn. (2006) 63:3–42. doi: 10.1007/s10994-006-6226-1
38. Stanley JS, King N, Burks AW, Huang SK, Sampson H, Cockrell G, et al. Identification and mutational analysis of the immunodominant IgE binding epitopes of the major peanut allergen ara h 2. Arch Biochem Biophys. (1997) 342(2):244–53. doi: 10.1006/abbi.1997.9998
39. Hazebrouck S, Patil SU, Guillon B, Lahood N, Dreskin SC, Adel-Patient K, et al. Immunodominant conformational and linear IgE epitopes lie in a single segment of ara h 2. J Allergy Clin Immunol. (2022) 150(1):131–9. doi: 10.1016/j.jaci.2021.12.796
40. Jin T, Guo F, Chen Y-W, Howard A, Zhang Y-Z. Crystal structure of ara h 3, a major allergen in peanut. Mol Immunol. (2009) 46(8):1796–804. doi: 10.1016/j.molimm.2009.01.023
41. Zhao LZ L, Zhang B, Robotham JM, Roux KH, Tang H. Identification of a common ara h 3 epitope recognized by both the capture and the detection monoclonal antibodies in an ELISA detection kit. PLoS One. (2017) 12(8):e0182935. doi: 10.1371/journal.pone.0182935
42. Cabanos C, Urabe H, Masuda T, Tandang-Silvas MR, Utsumi S, Mikami B, et al. Crystallization and preliminary x-ray analysis of the major peanut allergen ara h 1 core region. Acta Crystallogr Sect F Struct Biol Cryst Commun. (2010) 66(Pt 9):1071–3. doi: 10.1107/S1744309110029040
43. Asarnoj A, Nilsson C, Lidholm J, Glaumann S, Ostblom E, Hedlin G, et al. Peanut component ara h 8 sensitization and tolerance to peanut. J Allergy Clin Immunol. (2012) 130(2):468–72. doi: 10.1016/j.jaci.2012.05.019
44. Hurlburt BK, Offermann LR, McBride JK, Majorek KA, Maleki SJ, Chruszcz M. Structure and function of the peanut panallergen ara h 8. J Biol Chem. (2013) 288(52):36890–901. doi: 10.1074/jbc.M113.517797
45. Krause S, Reese G, Randow S, Zennaro D, Quaratino D, Palazzo P, et al. Lipid transfer protein (ara h 9) as a new peanut allergen relevant for a Mediterranean allergic population. J Allergy Clin Immunol. (2009) 124(4):771–8.e5. doi: 10.1016/j.jaci.2009.06.008
46. Lauer I, Dueringer N, Pokoj S, Rehm S, Zoccatelli G, Reese G, et al. The non-specific lipid transfer protein, ara h 9, Is an important allergen in peanut. Clin Exp Allergy. (2009) 39(9):1427–37. doi: 10.1111/j.1365-2222.2009.03312.x
47. van Hage M, Hamsten C, Valenta R. ImmunoCAP assays: pros and cons in allergology. J Allergy Clin Immunol. (2017) 140(4):974–7. doi: 10.1016/j.jaci.2017.05.008
48. Hoffmann HJ, Schmid J, Würtzen PA, Dahl R. Increased allergen-specific IgG4 levels are linked to suppressed IgE titers when analysed by ISAC but not by immunocap. Clin Transl Allergy. (2014) 4:O20. doi: 10.1186/2045-7022-4-S2-O20
49. Hamilton RG. Igg antibody interference diminishes accuracy of chip-based allergen-specific IgE assays. J Allergy Clin Immunol. (2010) 125(2):AB223 doi: 10.1016/j.jaci.2009.12.872
50. Wollmann EL C, Kundi M, Selb R, Niederberger V, Valenta R. Reduction in allergen-specific IgE binding as measured by microarray: a possible surrogate marker for effects of specific immunotherapy. J Allergy Clin Immunol. (2015) 136(3):806–9.e7. doi: 10.1016/j.jaci.2015.02.034
51. Mikus M, Zandian A, Sjöberg R, Hamsten C, Forsström B, Andersson M, et al. Allergome-wide peptide microarrays enable epitope deconvolution in allergen-specific immunotherapy. J Allergy Clin Immunol. (2021) 147(3):1077–86. doi: 10.1016/j.jaci.2020.08.002
52. Thörnqvist L, Sjöberg R, Greiff L, van Hage M, Ohlin M. Linear epitope binding patterns of grass pollen-specific antibodies in allergy and in response to allergen-specific immunotherapy. Front Allergy. (2022) 3:859126. doi: 10.3389/falgy.2022.859126
53. Kukkonen AK, Pelkonen AS, Makinen-Kiljunen S, Voutilainen H, Makela MJ. Ara h 2 and ara 6 are the best predictors of severe peanut allergy: a double-blind placebo-controlled study. Allergy. (2015) 70(10):1239–45. doi: 10.1111/all.12671
54. Hemmings GDT O, Radulovic S, Lack G, Santos AF. Ara h 2 is the dominant peanut allergen despite similarities with ara h 6. J Allergy Clin Immunol. (2020) 146(3):621–30.e5. doi: 10.1016/j.jaci.2020.03.026
55. Wang Y, Fu T-J, Howard A, Kothary MH, McHugh TH, Zhang Y. Crystal structure of peanut (Arachis hypogaea) allergen ara h 5. J Agric Food Chem. (2013) 61(7):1573–8. doi: 10.1021/jf303861p
Keywords: peanut allergy, epitope, immunoglobulin, oral immunotherapy, microarray, ISAC, linear peptide, machine learning
Citation: Rambo IM, Kronfel CM, Rivers AR, Swientoniewski LT, McBride JK, Cheng H, Simon RJ, Ryan R, Tilles SA, Nesbit JB, Kulis MD, Hurlburt BK and Maleki SJ (2023) IgE and IgG4 epitopes of the peanut allergens shift following oral immunotherapy. Front. Allergy 4:1279290. doi: 10.3389/falgy.2023.1279290
Received: 17 August 2023; Accepted: 30 October 2023;
Published: 29 November 2023.
Edited by:
Kazuyuki Nakagome, Saitama Medical University, JapanReviewed by:
Mats Ohlin, Lund University, SwedenMary H. Grace, North Carolina State University, United States
© 2023 Rambo, Kronfel, Rivers, Swientoniewski, McBride, Cheng, Simon, Ryan, Tilles, Nesbit, Kulis, Hurlburt and Maleki. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Soheila J. Maleki soheila.maleki@usda.gov
†These authors have contributed equally to this work and share senior authorship