Next Article in Journal
Exploring Inflammatory Asthma Phenotypes: Proteomic Signatures in Serum and Induced Sputum
Next Article in Special Issue
A Rare Case of TP63-Associated Lymphopenia Revealed by Newborn Screening Using TREC
Previous Article in Journal
Use of Enzymatically Converted Cell-Free DNA (cfDNA) Data for Copy Number Variation-Linked Fragmentation Analysis Allows for Early Colorectal Cancer Detection
Previous Article in Special Issue
Effects of Human Leukocyte Antigen DRB1 Genetic Polymorphism on Anti-Cyclic Citrullinated Peptide (ANTI-CCP) and Rheumatoid Factor (RF) Expression in Rheumatoid Arthritis (RA) Patients
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 25
Issue 6
10.3390/ijms25063499
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Blood Transcriptomics Identifies Multiple Gene Expression Pathways Associated with the Clinical Efficacy of Hymenoptera Venom Immunotherapy

by
Ajda Demšar Luzar
1,2,
Peter Korošec
1,3,4,
Mitja Košnik
1,5,
Mihaela Zidarn
1,5 and
Matija Rijavec
1,2,*
1
Laboratory for Clinical Immunology and Molecular Genetics, University Clinic of Respiratory and Allergic Diseases Golnik, 4204 Golnik, Slovenia
2
Biotechnical Faculty, University of Ljubljana, 1000 Ljubljana, Slovenia
3
Faculty of Pharmacy, University of Ljubljana, 1000 Ljubljana, Slovenia
4
Faculty of Medicine, University of Maribor, 2000 Maribor, Slovenia
5
Faculty of Medicine, University of Ljubljana, 1000 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(6), 3499; https://doi.org/10.3390/ijms25063499
Submission received: 29 January 2024 / Revised: 13 March 2024 / Accepted: 18 March 2024 / Published: 20 March 2024
(This article belongs to the Special Issue Genetic Studies of Immune-Related Diseases)
Browse Figures
Review Reports Versions Notes

Abstract

:
Allergen-specific venom immunotherapy (VIT) is a well-established therapy for Hymenoptera venom allergy (HVA). However, the precise mechanism underlying its clinical effect remains uncertain. Our study aimed to identify the molecular mechanisms associated with VIT efficiency. We prospectively included 19 patients with HVA undergoing VIT (sampled before the beginning of VIT, after reaching the maintenance dose, one year after finishing VIT, and after a sting challenge) and 9 healthy controls. RNA sequencing of whole blood was performed on an Illumina sequencing platform. Longitudinal transcriptomic profiling revealed the importance of the inhibition of the NFκB pathway and the downregulation of DUX4 transcripts for the early protection and induction of tolerance after finishing VIT. Furthermore, successful treatment was associated with inhibiting Th2, Th17, and macrophage alternative signalling pathways in synergy with the inhibition of the PPAR pathway and further silencing of the Th2 response. The immune system became activated when reaching the maintenance dose and was suppressed after finishing VIT. Finally, successful VIT restores the immune system’s balance to a state similar to that of healthy individuals. Our results underline the important role of the inhibition of four pathways in the clinical effect of VIT: Th2, Th17, NFκB, and macrophage signalling. Two biomarkers specific for successful VIT, regardless of the time of sampling, were C4BPA and RPS10-NUDT3 and should be further tested as potential biomarkers.

Graphical Abstract

1. Introduction

Hymenoptera venom allergy (HVA) is an IgE-mediated hypersensitivity reaction following a honeybee, vespid, or ant sting [1,2]. An allergic reaction to Hymenoptera species presents with severe manifestations and is potentially life-threatening [3]. HVA varies from large local reactions at the sting site to systemic reactions [4]. Systemic reactions vary significantly in severity, from moderate reactions consisting of generalised skin symptoms to severe life-threatening anaphylactic reactions affecting the cardiac and respiratory systems. The prevalence of systemic reactions is estimated to be up to 7.5% in adults and up to 3.4% in children, with anaphylaxis in 0.3–42.8% of cases [2,5,6]. The incidence of anaphylaxis has been increasing over the years, and its mortality rate is estimated to be between 0.03 and 0.48 per million population per year [7]. However, mortality may have been underestimated due to unrecognised stings in unexplained causes of death [8,9].
The only treatment is allergen-specific venom immunotherapy (VIT), which is reported to be effective in 77–84% of patients treated with honeybee venom and in 91–96% of patients treated with vespid venom [2,10]. Immunotherapy aims to restore immune tolerance and thus eliminate systemic allergic reactions after insect stings [2]. In contrast to aero and food immunotherapies, the immune tolerance established during VIT is considered to be lifelong, even after the discontinuation of treatment [11]. However, VIT is a long-lasting treatment, and in some cases, it should not be discontinued. Patients with known risk factors, namely, an initial severe sting reaction, clonal mast cell disorders, and/or elevated basal serum tryptase (>11.4 µg/L) [2,12,13,14], should receive prolonged or even lifelong treatment.
Many studies have focused on detecting biomarkers for VIT efficiency evaluation. Currently, laboratory parameters used as diagnostic methods to detect sensitisation are specific IgE (sIgE) levels and skin prick tests. However, there is no evidence that monitoring sIgE levels can predict the success of therapy [2]. Furthermore, IgG4 levels were studied. Its levels can be used for an evaluation of immunological responsiveness to VIT. It was suggested that not levels but rather functional activity might correlate with VIT effectiveness [15,16]. Previous studies evaluated the activation of effector cells using the basophil activation test and concluded that successful VIT decreases test sensitivity without changing the reactivity [17,18,19]. Other studies have focused on using serum CD30 levels to indicate a Th2 immune response. A decrease in serum CD30 levels may be considered an effective prognostic factor [20,21]. Most recent studies have highlighted the importance of the type of allergen the patient is sensitised to. Predominant sensitisation to the allergen Api m 1 showed no increased risk of treatment failure, whereas for patients predominantly sensitised to the allergen Api m 10, the risk of VIT failure increased. Therefore, molecular diagnostics has promising potential for the risk stratification of VIT [22,23,24]. But to date, no biomarkers have been established in clinical practice, making a controlled sting challenge the gold standard for evaluating venom tolerance, indicating clinical protection [2,25].
With the progress of science, the rapidly evolving field of genetics could be used to explore mechanistic changes in patients undergoing VIT and further propose potential biomarkers. Only a few transcriptome studies have been conducted in the research field of VIT. One of the first ones was performed by Tsicopoulus A., who investigated the mRNA profiles of IL-4 and IFNγ in enriched T lymphocytes. A significant decrease in IL-4 mRNA expression was observed after 90 days of VIT. Similarly, an increase in IFNγ expression was detected during VIT [26]. Later, in 2006, a microarray analysis of differentially expressed genes (DEGs) in peripheral blood mononuclear cells was performed. Three independent groups of patients were compared: the control group, patients after less than a year of VIT, and patients after five to six years of VIT. Osteopontin, a cytokine known to modulate Th1 activity, was significantly upregulated in patients after VIT [27,28]. In the following years, Niedoszytko et al. studied gene expression reflecting long-term protection in VIT. After comparing a group of patients with successful VIT to a group with unsuccessful VIT, a list of 18 genes forming part of a predictive model of successful VIT was proposed [29]. Furthermore, the same authors compared a group of patients before VIT with patients in the maintenance phase and with patients after VIT. Only 6 of 18 previously described genes were significantly expressed between the first and last aforementioned groups. Signature genes differentially expressed after reaching the maintenance phase of VIT, and signature DEGs still observed in patients after VIT were described [30]. Later, in 2018, a genome-wide analysis of expression before VIT and one year after VIT was conducted. It was the first study to include paired samples from the same patients during VIT. Although the findings did not align with previously described gene expression patterns, they were important, as they described the upregulation of leukocyte-specific genes expressed in antigen-presenting cells: macrophages and monocytes. The authors subsequently suggested that changes in T-cell subsets and the simultaneous restoration of the immune system’s innate branch may play a key role in establishing venom tolerance [31]. Over the last few years, an evolving method of single-cell RNA sequencing has been used for biomarker research [32,33,34], but not yet in the field of HVA or VIT.
Even though it is well documented that VIT provokes venom tolerance in most treated patients, the exact underlying mechanism remains unclear. We aimed to establish a less invasive method than the sting challenge to assess the effect of VIT. With follow-up transcriptome analysis of paired blood samples, we tried to contribute to the multi-level fingerprint for VIT monitoring and the evaluation of its effectiveness. By gaining new insights into underlying mechanisms, novel possibilities for biomarker research are uncovered.

2. Results

2.1. Characterisation of Transcriptome Fingerprint Characteristic of Successful VIT

Of 19 patients enrolled, 14 tolerated exposure to the allergen at the sting challenge and were classified as successfully completing the treatment, while patients with recurrent systemic reaction to sting, evaluated as Ring and Messmer grading I-II, were classified as treatment failure. Interestingly, only one patient with successful VIT had complications during the build-up phase of VIT, whereas all five patients with subsequent treatment failure experienced adverse reactions during the beginning of treatment.
We aimed to characterise the transcriptome fingerprint and its differences between patients who tolerated a sting after VIT and those with treatment failure, enabling a further understanding of the mechanisms underlying VIT and thus contributing to the search for potential biomarkers of VIT efficiency. To our knowledge, this is the first transcriptome analysis designed as a follow-up study with paired samples before VIT, after reaching the maintenance dose, one year after finishing VIT, and right after the following sting challenge. The outcome of transcriptome analysis is a large amount of data, which need to be processed thoroughly. Thus, biomarker selection filtration criteria and a comparison of selected DEGs with previously published data were necessary. Accordingly, we first performed a follow-up analysis at different time points of VIT, followed by a detailed examination of overlapping DEGs (Figure 1) and a comparison to previously proposed biomarkers. Since immune cells represent only 1% of all blood cells [35], the transcriptome changes regarding the immune system were moderate; therefore, a p-value < 0.05 was used as a filtration criterion instead of a stricter false discovery rate p-value. Combined with the stringent criterion of a fold change higher than five, we gained a strong forecast of the selected DEGs’ significance. A detailed list of DEGs characteristic of successful VIT after reaching the maintenance dose (maintenance dose vs. before; maintenance dose vs. after; 76 DEGs), DEGs specific for the maintenance dose and after VIT (maintenance dose vs. before; after vs. before; 91 DEGs), and DEGs specific only for the time point after the treatment (after vs. before; maintenance dose vs. after; 254 + 49 DEGs) are presented in Supplementary S1, Supplementary S2, and Supplementary S3, respectively. DEGs overlapping with treatment failure (14) are unspecific and were excluded from further analysis. Further results are focused on patients with successful VIT.

2.2. The Immune System Seems to Be Activated after Reaching the Maintenance Dose and, Conversely, Significantly Suppressed after Finishing VIT

After the maintenance dose of 100 µg is reached, the patient is believed to be protected from subsequent Hymenoptera stings. The mechanism behind gained protection differs from the mechanisms acquired through the following years of VIT. DEGs (76) specific for the maintenance dose that are different from those at time points before and after immunotherapy are presented in Supplementary S1. We would like to highlight three upregulated immune-related DEGs: SERPINB4, HLA-U, and IGLV2–5 (Table 1). Interestingly, the serin protease inhibitor encoded by the SERPINB4 gene seems to be upregulated only at the time point of the maintenance dose. SERPINB4 has been previously associated with early inflammation in atopic dermatitis and psoriasis [36]. Similarly, the upregulation of the human leukocyte antigen can be observed, which has been directly associated with immune stimulation against external antigens [37]. DEGs that describe the immune system characteristics established after reaching the maintenance dose and that remain the same even after finishing immunotherapy are mainly downregulated (75/91; 83%). DEGs closely related to the immune system (referred to as hub genes) are presented in Table 2. Overall, not only the proportions of down-/upregulated genes but also their functional characterisation suggest that the immune system is activated after reaching the maintenance dose of the treatment and suppressed later after finishing VIT.

2.3. NFκB Pathway Inhibition and Downregulation of DUX4 Transcripts Are Important for Early Protection and Tolerance Induction

The two DEGs specifically downregulated at both time points (maintenance dose and after VIT) are the pseudogenes DUX4L10 and DUX4L12. Furthermore, the DUX4, DUX4L2, and DUX4L14 transcripts were significantly downregulated after finishing VIT. The DUX4 gene has been previously investigated in the context of other diseases [39]. Its transcriptomic profiling revealed its involvement in inflammation, cellular migration, and chemotaxis [40]. The pseudogenes DUX4L10 and DUX4L12 were significantly downregulated after reaching the maintenance dose (FC = −240, p-value = 6.46 × 10−7; FC = −90, p-value = 5.48 × 10−5, respectively), suggesting the limitation of the inflammatory response. Interestingly, the DUX4 gene was significantly downregulated only after finishing VIT (after vs. before; fold change = −18.5, p-value = 0.0007).
The most significant DEG expressed as early as reaching the maintenance dose, which persisted even one year after finishing VIT, is the downregulated IKBKGP1 transcript, influencing the NFкB pathway. NFкB represents a family of transcription factors that regulate the expression of hundreds of genes, most of which are involved in immune and inflammatory responses [41]. Its signalling is essential to effective innate and adaptive immune responses, playing critical roles in the development, survival, and activation of B lymphocytes. There are two pathways of NFкB signalling, namely, canonical and non-canonical [42]. In the canonical pathway, upon the stimulation of B-cell receptors, CD40, or Toll-like receptors, the IкB kinase (IKK) complex liberates NFкB dimers, which translocate to the nucleus, where they control the transcription of target genes. The NF-kappa-B essential modulator (also known as NEMO) is the regulatory subunit of the IKK complex. NEMO is encoded by the Inhibitor of the Nuclear Factor Kappa B Kinase Subunit Gamma (IKBKG) gene. Its pseudogene IKBKGP1 is highly significant for patients who tolerate stings. It is important to note that the significance of pseudogenes is being discovered and has been highlighted in various studies [43,44,45]. The transcript IKBKGP1 was 5 times less expressed after reaching the maintenance dose (p-value = 0.01) and 25 times less expressed after immunotherapy (p-value = 0.00001). Its significant downregulation after VIT suggests lower immune response activation. Comparing successful and treatment failure patients, no difference in the transcript level can be observed before the initiation of VIT. In contrast, its significant downregulation is characteristic of the group with successful VIT after treatment (Figure 2).
Interestingly, the connection between VIT and the NFкB pathway was highlighted in a study by Kempinski et al. [46]. According to their research, the only gene that changed significantly during the build-up phase of VIT was COMMD8, which influences IкB degradation and, thus, the release of NFкB for nuclear translocation [46,47]. Moreover, the downregulation of the COMMD8 transcript was also a part of the prediction model of successful VIT in a study performed by Niedoszytko et al. [29], further highlighting the vital role of the NFкB pathway in mechanisms of allergen tolerance established during VIT.

2.4. Th2- and Th17-Related DEGs Are Significantly Downregulated after Successful Immunotherapy

Of 410 DEGs in the group with successful VIT, those closely related to the immune system are described in Table 3. Of these, 96 overlapped with the time point of the maintenance dose. Our analysis of follow-up samples’ differential expression patterns before VIT, after reaching the maintenance dose, and after finishing VIT in the group with successful VIT shows the apparent downregulation of Th2- and Th17-related immune response transcripts. Furthermore, only one-quarter of DEGs after immunotherapy were upregulated and are primarily involved in negative immune system regulation. Hence, we assume that immunotherapy does not promote specific immune system pathways but, conversely, downregulates excessive immune processes. After successfully finishing immunotherapy, DEGs involved in the Th2 immune pathway (EBI3, UMOD), the Th17 immune pathway (PHBP6), tryptophan metabolism (TDO2), the production of inflammatory mediators (PLA2G10), and the innate branch of the immune system (CPN1, MMD2, NLRP8) were significantly downregulated. Notably, the PHBP6 transcript, a pseudogene of prohibitin (PHB1), was significantly downregulated after immunotherapy (after vs. before; fold change = −24, p-value = 0.001). Prohibitin is known to interact with STAT3 to affect IL17 secretion in Th17 cells. It is also required to activate the NFκB signalling pathway [48]. The downregulation of the Th17 pathway seems to be one of the prevailing mechanisms important for establishing allergen tolerance.
Furthermore, the important involvement of B lymphocytes can be observed through possible isotype switching by the upregulation of IGKV2D-28, which is involved in antigen recognition.
On the other hand, fewer changes are seen through transcriptome analysis in the group of patients with treatment failure. One of the downregulated genes is the IFI44L transcript, which is known to promote macrophage differentiation and facilitate inflammatory cytokine secretion. Furthermore, the TRIL transcript, a component of the TLR4 signalling complex that mediates the innate immune response, was significantly upregulated after immunotherapy. Taken together, VIT seems to impact both humoral and innate immune system pathways in both groups of patients, namely, successful treatment and treatment failure, although seemingly to a much lower extent in the group of patients with treatment failure.

2.5. Successful Venom Immunotherapy Seems to Restore the Balance of the Immune System to a State Similar to That of Healthy Individuals

Comparing the transcriptome signatures of patients with successful treatment or treatment failure with healthy controls, a higher proportion of similarity can be observed between healthy controls and patients with successful VIT after immunotherapy (Figure 3a,b). In the group with successful VIT, there were 1768 DEGs compared to healthy controls before the treatment, whereas there were only 476 DEGs after VIT. In contrast, in the group with treatment failure, there were 1745 DEGs before VIT and 1012 DEGs after VIT compared to healthy controls. Hence, patients tolerating stings have a more similar expression pattern to healthy controls than patients with treatment failure (Figure 3c). This observation was previously described by Niedoszytko et al. [30], who discovered that immunotherapy seems to re-establish a balance similar to that in non-allergic individuals by silencing the immune system. Also noteworthy is the comparison of transcriptome signatures between successful VIT and treatment failure before and after the treatment, which are relatively similar, with only slight significant differences. Nevertheless, the dissimilarity is more notable after the treatment than before VIT (Figure 3a,b).

2.6. Two Biomarkers Specific for Successful VIT, Regardless of the Time of Sampling, Are C4BPA and RPS10-NUDT3

The most important outcome of differential gene expression between groups at different time points of VIT is the identification of two DEGs referred to as biomarkers of successful VIT. The two transcripts significantly expressed in the group with successful VIT were C4b-binding protein alpha chain transcript (C4BPA) and RPS10-NUDT3. The expression of the C4BPA transcript was 10 times higher before and after immunotherapy in comparison to healthy controls, and its expression was 14 times and 9 times higher before and after immunotherapy in comparison to treatment failure, respectively. Conversely, the RPS10-NUDT3 transcript was significantly downregulated in comparison to healthy controls and those with treatment failure (Figure 3a). C4b-binding protein is a complement inhibitor of the uncontrolled activation of the classical and lectin complement pathways. Besides complement inhibition, it has been studied as a regulator of excessive inflammation in chronic diseases, like systemic lupus erythematosus, and in cancer [49]. To our knowledge, C4BPA and RPS10-NUDT3 have not been connected to HVA or VIT before.

2.7. Inhibition of Macrophage Alternative Signalling Pathway in Synergy with Inhibition of PPAR Pathway Results in Silencing of Th2 Immune Response Characteristic of Successful VIT

Many differences in active pathways between patients with successful VIT and patients with treatment failure can be seen when characterising DEGs. The major blood-related pathways before the initiation of VIT demonstrated the activation of IL-6 signalling, TREM-1 signalling, and macrophage alternative activation signalling pathways in successful VIT patients. In contrast, IL-10 signalling, CLEAR signalling, and PPAR signalling pathways seem to be inhibited compared to patients with treatment failure. After immunotherapy, pathways related to B lymphocytes appear to be the most enriched, namely, B-cell receptor signalling, PI3K signalling, and FcµRIIB signalling. Outstanding is the pathway activation of the FcµRIIB receptor expressed on B lymphocytes, the only inhibitory Fc receptor [50]. Its crosslinking with the B-cell receptor (BCR) results in the inhibition of activation, proliferation, antigen internalisation, and antibody secretion [51]. Furthermore, the immune response is silenced by PI3K signalling, which has been shown to inhibit mast cell degranulation [52].
The importance of macrophages and the subsequent involvement of the immune system’s innate branch in establishing venom tolerance was previously observed by Karpinski et al. [31]. Our results further support their theory. The macrophage alternative activation signalling pathway is shown to be inhibited in successful VIT after the treatment, alongside the peroxisome proliferator-activated receptor (PPAR) and retinoid X receptor signalling pathways. The canonical pathway of alternative macrophage activation signalling involves the activation of STAT6 by the Th2 cytokines IL-4 and IL-13. The transcriptional synergy between PPAR/retinoid X receptor regulators and STAT6 is known to sustain the immune effector response [53]. Several studies have highlighted the important role of the activated PPAR pathway in promoting the Th2 immune response associated with allergic disease [54,55,56]. The downregulation of macrophage activation and the downregulation of the PPAR pathway reflect a lower Th2 immune response. These results again emphasise the importance of downregulating the immune system and establishing allergen tolerance gained over the years of VIT.

2.8. Protein–Protein Interaction before and after Venom Immunotherapy Identifies Different Master Regulators

The protein–protein interaction network analysis of selected DEGs in successful VIT vs. treatment failure was performed using Ingenuity Pathway Analysis to identify master regulators and their concomitant causal networks. Before VIT, 39 master regulators were identified. Of these, 27 predicted activation, and 12 predicted inhibition. The causal network with the most target-connected regulators is the EGFR network (Supplementary S4). The EGFR regulator impacts the NFκB and STAT3 pathways, promoting an inflammatory response. After VIT, nine master regulators were identified. Of these, five were related to activation, and four predicted inhibition. The causal network with the most target-connected regulators is RASSF1 (Supplementary S5). The RASSF1 master regulator significantly inhibits the NFκB subunit (RELA), further supporting the follow-up analysis suggesting that the significant inhibition of the NFκB pathway after VIT is characteristic of patients with successful VIT.

2.9. Different Mechanisms after Sting Challenge in Groups with Successful VIT and Treatment Failure

Differences in the immune response after a sting challenge between the group of patients tolerating the sting and the group of patients with treatment failure are evident. After the sting challenge, successful VIT establishes a tolerance to the allergen using different mechanisms, which can be seen in transcriptome signatures, as highly expressed DEGs in the group with treatment failure are completely downregulated in successful VIT (Figure 4). The major blood-related pathways in the group with successful VIT after the sting showed activated pathways of phagosome formation, the PI3K cascade, and potassium channels. The PI3K signalling cascade, inhibiting mast cell degranulation, is characteristic of patients with successful VIT tolerating a sting. Conversely, functional and pathway analysis in the group with treatment failure after the sting resulted in enriched pathways in B-cell receptor signalling, FCƐRI signalling, FCγR-dependent phagocytosis, and the binding and uptake of ligands by scavenger receptors.

3. Discussion

We used transcriptome signatures to characterise differences between patients with successful VIT, who tolerated a sting challenge, and patients with treatment failure. Transcriptome analysis offers a bridge to the gap of known VIT mechanisms of gained allergen tolerance and opens new biomarker possibilities. Our results underline the important role of four immunological mechanisms: inhibition of Th2, Th17, and NFκB pathways and the involvement of macrophages in establishing venom tolerance; they also identify C4BPA and RPS10-NUDT3 as possible biomarkers of successful VIT.
To our knowledge, this is the first study to use follow-up samples before VIT, after reaching the maintenance dose, after finishing VIT, and after the sting challenge. Investigating patients’ transcriptome signatures, we outlined DEGs that were characteristic of patients who tolerated a sting in the controlled sting challenge. Since allergy and immunotherapy both consist of various complicated, interacting immune system mechanisms, we assumed that the specific dynamics of differentially expressed genes indicating different activated or inhibited pathways could be specific for different outcomes of VIT rather than a single measurement of a specific DEG. With this study being the first follow-up study on Hymenoptera venom immunotherapy, various previously described mechanisms have been confirmed through differential expression patterns.
After reaching the maintenance dose, three transcripts, namely, SERPINB4, HLA-U, and IGLV-2, seem to describe the activated state of the immune system. SERPINB4 has been previously associated with early inflammation in atopic dermatitis and psoriasis [36]. Similarly, human leukocyte antigen has been directly associated with immune stimulation against external antigens [37]. In contrast, DEGs that are significantly expressed after reaching the maintenance dose and maintained even after finishing VIT show significant inhibition of the immune system. Overall, the immune response seems to be induced after reaching the maintenance dose and, conversely, suppressed after finishing VIT. The DUX4 gene is significantly downregulated after reaching the maintenance dose and after finishing VIT. Its transcriptomic profiling in previous studies revealed its involvement in inflammation, cellular migration, and chemotaxis [40]. The downregulation of the DUX4 gene suggests a lower inflammatory state, which enables early protection and, later, allergen tolerance. Furthermore, the NFκB pathway is significantly inhibited. NFκB signalling is essential to effective innate and adaptive immune responses. It plays important roles in the development, survival, and activation of B lymphocytes [42]. The inhibition of the NFκB pathway through the downregulation of the IKBKGP1 transcript suggests overall lower immune activation. The set of genes proposed as a prediction model for VIT evaluation, as described by Niedoszytko et al. [30], was analysed. However, of 18 proposed genes, only the COMMD8 transcript indirectly coincides with our observation of the downregulation of the NFκB pathway, which furthermore highlights the vital role of the NFκB pathway in gaining early protection after reaching the maintenance dose and allergen tolerance after finishing VIT.
In the group of patients with successful VIT, a previously suggested shift from the Th2 to Th1 immune response [57] can be observed through a trend of primarily downregulated DEGs. In addition, the downregulation of DEGs related to the Th17 immune pathway (PHBP6), tryptophan metabolism (TDO2), the production of inflammatory mediators (PLA2G10), and innate immune system branches (CPN1, MMD2, NLRP8) is highly significant for successful VIT. According to our transcriptome analysis, the silencing of the Th2 immune response seems to be a consequence of the inhibition of the macrophage alternative signalling pathway in synergy with the inhibition of the PPAR pathway. The importance of the involvement of the innate immune system branch, specifically macrophage involvement, was previously highlighted by Karpinski et al. in the analysis of gene expression one year after VIT [31]. Furthermore, the inhibition of the Th17 pathway seems to be one of the prevailing mechanisms important for establishing allergen tolerance. PHBP6, a pseudogene of prohibitin, is highly downregulated. It is known to interact with STAT3 to affect IL17 secretion in Th17 cells. Also important is the involvement of B lymphocytes, which can be observed through possible isotype switching by the upregulation of the IGKV2D-28 transcript. IGKV2D-28 is involved in antigen recognition.
It is important to state that most DEGs were downregulated after immunotherapy, as was similarly proposed by Konno et al. [27,28] and Niedoszytko et al. [29]. In contrast, Karpinski et al. [31], in a follow-up study, discovered a bigger proportion of upregulated genes one year after VIT. In our study, most of the DEGs were downregulated at both time points of VIT—after reaching the maintenance dose and after finishing VIT (maintenance dose: 385/486; 80%; after VIT: 317/410; 77%). Only roughly one-quarter of all DEGs after immunotherapy were upregulated and are primarily involved in negative immune system regulation. Hence, we assume that immunotherapy does not promote specific immune system pathways but, conversely, downregulates excessive immune processes.
Taken together, VIT seems to impact both humoral and innate immune systems in both groups of patients but seemingly does not result in immune tolerance in the group of patients with treatment failure.
Comparing transcriptome signatures of the three groups—successful VIT, treatment failure, and healthy controls—we identified similar expression patterns between healthy controls and patients with successful VIT after immunotherapy. Immunotherapy seems to restore the balance of the immune system to a state similar to that of non-allergic individuals. Many differences between patients with successful VIT and those with treatment failure can be observed. The major blood-related pathways in successful VIT patients before the initiation of VIT demonstrated the activation of IL-6 signalling, TREM-1 signalling, and macrophage alternative activation signalling pathways and the inhibition of IL-10 signalling, CLEAR signalling, and PPAR signalling pathways. In contrast, after immunotherapy, pathways regarding B lymphocytes appear to be the most enriched, namely, B-cell receptor signalling, PI3K signalling, and FcµRIIB signalling. The most outstanding pathway after immunotherapy in patients with successful VIT is the activation of the FcµRIIB receptor expressed on B lymphocytes, which is the only inhibitory receptor. Its crosslinking with the BCR results in the inhibition of activation, proliferation, antigen internalisation, and antibody secretion [51]. Furthermore, the macrophage alternative activation signalling pathway is shown to be activated before VIT and, conversely, inhibited after VIT in the group with successful VIT. The canonical pathway of alternative macrophage activation signalling involves the activation of STAT6 by the Th2 cytokines IL-4 and IL-13. Additionally, after VIT, the PPAR and retinoid X receptor signalling pathways seem to be inhibited. The transcriptional synergy between PPAR/retinoid X receptor regulators and STAT6 is known to sustain the immune effector response. The inhibition of the macrophage activation signalling pathway and the inhibition of the PPAR pathway reflect a lower Th2 immune response.
The important result of the transcriptome analysis of the follow-up samples is the identification of two biomarkers of successful VIT, the C4BPA and RPS10-NUDT3 transcripts, that differentiate patients from those with treatment failure and from healthy controls before and after VIT. C4BPA acts as a complement inhibitor and has been studied as a regulator of excessive inflammation in chronic diseases [49]. To our knowledge, neither of the two transcripts has ever been associated with HVA or VIT before. Future confirmational studies on larger cohorts of patients are required to evaluate the potential of the two proposed biomarkers.
The previously proposed biomarkers of VIT efficiency (osteopontin, CD30, IL-4, IFNγ) were studied in detail and were not significantly expressed after reaching the maintenance dose or after finishing immunotherapy.
The transcriptome signature of the sting challenge showed clear differences in expression patterns between patients with successful VIT and those with treatment failure. Many more DEGs were upregulated in the group with successful VIT, but those silenced were shown to be significantly upregulated in the group of patients with treatment failure. Even though VIT suppresses the overall immune system in patients who tolerate the sting, it seems that upregulated DEGs also suppress inappropriate responses of the immune system, and those downregulated are responsible for allergic reactions. The major blood-related pathways in the group with successful VIT after the sting showed activated pathways of phagosome formation, PI3K cascade, and potassium channels. The PI3K cascade is known to inhibit mast cell degranulation, contributing to allergen tolerance.
The results of our analysis coincide with other studies performed in the field of immunotherapy to aero and food allergens (AIT), such as house dust mites, grass pollen, and eggs. Wide differences have been pointed out between VIT and AIT. To date, the mechanisms of hypersensitivity reactions that differentiate between respiratory and oral allergies from HVA have not been clarified. The most recent remarkable study by Huang et al. performed on house dust mite-allergic patients treated with AIT, similarly to our study, revealed the significant downregulation of the innate immune response, the downregulation of NFκB and PPAR pathways, and a lower involvement of Th17 signalling with subsequent decreases in Th17 and Th2 cytokines [58]. Furthermore, a single-cell RNA analysis on grass pollen-allergic patients demonstrated the downregulation of genes associated with Th2 signalling and the downregulation of antigen-presenting pathways in B lymphocytes [34]. Notably, a transcriptome study on children receiving oral immunotherapy for egg allergy also highlighted that the majority of DEGs after 8 months of AIT were downregulated, and some of them played an important role in IL-17 signalling [59]. Even though the differences between VIT and AIT for aero and food allergens have been pointed out, our results of the underlying mechanisms of gained allergen tolerance coincide and are very important for further research directions.
We want to point out that the potential limitation of our study is the fact that we used whole-blood mRNA for whole-transcriptome analysis. Since significant changes in gene expression in VIT are likely related to a small subset of cells, single-cell analysis would be reliable for the further evaluation of the presented results, as would polymerase chain reaction studies for the confirmation of identified gene fingerprints. Further, since the sting challenge is highly unpleasant and, in the case of treatment failure, also life-threatening for the patient, a possible limitation of the study is the small number of patients with treatment failure. An advantage of our study is the follow-up samples from the same patients monitored from the beginning of VIT until the end of treatment. Furthermore, the evaluation of VIT success was determined using the controlled sting challenge.
Taken together, our results underline the important role of four main known immunological mechanisms, namely, the inhibition of Th2, Th17, and NFκB pathways and the involvement of macrophages in establishing venom tolerance. We suggest further evaluating two possible biomarkers of successful VIT prognosis: C4BPA and RPS10-NUDT3.
Since the controlled sting challenge is very unpleasant for HVA patients, future studies should strive to find a biomarker of VIT efficacy evaluated from patients’ blood. We suggest that a biomarker fingerprint rather than a single measurement could help predict VIT outcomes.

4. Materials and Methods

4.1. Study Group

In this prospective study, we enrolled 19 patients with Hymenoptera venom allergy undergoing VIT at the University Clinic of Respiratory and Allergic Diseases Golnik (Slovenia) and 9 healthy controls with no diagnosed allergies. The study group characteristics are described in Table 4. According to the clinical history and skin and/or blood testing, honeybee VIT was indicated. Patients underwent honeybee ultra-rush VIT as previously described [2,25,60,61]. Whole-blood samples were collected before the beginning of VIT, after reaching the maintenance dose, one year after finishing VIT, and after a controlled sting challenge [2] (Figure 5). A sting challenge is still the most reliable method for evaluating VIT efficiency and the gold standard for clinicians. The sting challenge confirmed successful treatment in 14 patients, while 5 patients were considered as treatment failure. Treatment failure was characterised as Ring and Messmer grading I-II, and the reaction was immediately treated with epinephrine. The study was approved by the Slovenian National Medical Ethics Committee (KME 0120-443/2020/3). All patients gave their written informed consent to participate in the study.

4.2. Allergy Diagnostic Tests

In the laboratory evaluation, serum concentrations of specific honeybee IgE and total IgE were measured before the beginning of VIT on an Immulite (Siemens Healthcare GmbH, Erlangen, Germany). Sensitisation was defined as a sIgE level of 0.35 kIU/L or higher. The descriptive statistics are presented as medians and min–max ranges for measurement data and percentages for categorical data. Fisher Exact, Chi-square, or Mann–Whitney tests were used as appropriate to determine statistically significant differences using GraphPad Prism 10 software (version 10.2.1 for Windows; GraphPad Software, San Diego, CA, USA).

4.3. RNA Samples

Whole-blood samples containing an RNA-stabilising reagent were collected in PAXgene blood tubes (Qiagen, Hilden, Germany), which contain a proprietary reagent composition based on patented RNA stabilisation technology. Intracellular RNA was then isolated from human whole-blood leukocytes using the PAXgene Blood miRNA Kit on a fully automated QIAcube system (Qiagen) according to the manufacturer’s instructions. The integrity of RNA was determined using an RNA 6000 Nano LabChip kit (Agilent Technologies, Santa Clara, CA, USA) and an Agilent2100 Bioanalyser (Agilent Technologies). All samples with RIN integrity values greater than six were suitable for further sequencing.

4.4. RNA Sequencing and Bioinformatics Analysis

Sequencing libraries were prepared using the TruSeq Stranded TotalRNA Library workflow with Ribo-Zero-Globin (Illumina, San Diego, CA, USA). RNA sequencing (RNAseq) was performed on the Illumina NovaSeq6000 platform with 100 bp paired-end reads by Macrogen Inc. (Seoul, Republic of Korea). RNA sequencing data were then uploaded to the CLC Genomics Workbench (v.23.0.2, Qiagen Bioinformatics) [62]. Briefly, quality trimming to exclude low-quality reads, adapter trimming, alignment with Homo Sapiens (reference sequence_hg19, gene track_ensembl v87, mRNA track_ensembl v87), and global TMM normalisation of the reads were performed. The ENSEMBL gene annotation reference was used, which includes annotations of coding and non-coding genes (long non-coding RNAs, pseudogenes, small RNAs, and IG/TR transcripts) [63]. Genes with low expression (mean expression < 0.01) were filtered out. Differential gene expression was analysed within groups at different time points and between groups using the Wald test against the control group. Filtration criteria for differentially expressed genes were set to a p-value of less than 0.05 and a fold change greater than five. Similarly, filtration criteria to identify differentially expressed genes after the sting challenge were set to a p-value less than 0.05 and a fold change of more than two. Furthermore, a network and pathway analysis was performed using Ingenuity Pathway Analysis (IPA, Qiagen).

5. Conclusions

The transcriptome analysis of follow-up samples from the same patients before VIT, after reaching the maintenance dose, one year after finishing VIT, and after a controlled sting challenge reveals the significant downregulation of DEGs related to the immune system. The significant inhibition of the Th2, Th17, and NFκB pathways is characteristic of successful VIT. Further evaluation of two possible biomarkers of successful VIT prognosis, namely, C4BPA and RPS10-NUDT3, is needed. Overall, the immune system seems to be activated after reaching the maintenance dose of the treatment and later suppressed after finishing VIT.
Since allergy is a complex disease, we suggest a biomarker fingerprint rather than a single biomarker, which could help predict VIT outcomes. Downregulated transcripts indicating lower immune activity and previously proposed biomarkers combined could provide a stronger forecast of VIT efficiency.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms25063499/s1.

Author Contributions

Conceptualisation, M.R. and P.K.; data curation, M.R., M.Z. and A.D.L.; methodology, M.R., A.D.L., M.K. and M.Z.; formal analysis, A.D.L.; investigation, A.D.L., resources, A.D.L., M.R., P.K., M.Z. and M.K.; writing—original draft preparation, A.D.L.; writing—review and editing, M.R. and P.K.; visualisation, A.D.L., M.R. and P.K.; supervision, M.R., M.K. and P.K.; project administration, M.R.; funding acquisition, M.R., M.K. and P.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Slovenian Research and Innovation Agency (grant nos. P3-0360 and 53537).

Institutional Review Board Statement

This study was approved by the Slovenian National Medical Ethics Committee (KME 0120-443/2020/3).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data are contained within the article and Supplementary Files. Any additional supporting data are available from the corresponding author upon reasonable request.

Acknowledgments

The graphical abstract and Figure 5 were created using Biorender.com (accessed on 23 January 2024).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Lydyard, P.; Whelan, A.; Fanger, M.W. Immunology, 2nd ed.; Garland Science/Bios Scientific Publishers: New York, NY, USA, 2004. [Google Scholar]
  2. Sturm, G.J.; Varga, E.; Roberts, G.; Mosbech, H.; Bilò, M.B.; Akdis, C.A.; Antolín-Amérigo, D.; Cichocka-Jarosz, E.; Gawlik, R.; Jakob, T.; et al. EAACI guidelines on allergen immunotherapy: Hymenoptera venom allergy. Allergy Eur. J. Allergy Clin. Immunol. 2018, 73, 744–764. [Google Scholar] [CrossRef]
  3. De Graaf, D.C.; Aerts, M.; Danneels, E.; Devreese, B. Bee, wasp and ant venomics pave the way for a component-resolved diagnosis of sting allergy. J. Proteom. 2009, 72, 145–154. [Google Scholar] [CrossRef]
  4. Biló, B.M.; Rueff, F.; Mosbech, H.; Bonifazi, F.; Oude-Elberink, J.N.G.; The EAACI Interest Group on Insect Venom Hypersensitivity. Diagnosis of Hymenoptera venom allergy. Allergy Eur. J. Allergy Clin. Immunol. 2005, 60, 1339–1349. [Google Scholar] [CrossRef]
  5. Bilo, B.; Bonifazi, F. Epidemiology of insect-venom anaphylaxis. Curr. Opin. Allergy Clin. Immunol. 2008, 8, 330–337. [Google Scholar] [CrossRef] [PubMed]
  6. Bilò, M.B. Anaphylaxis caused by Hymenoptera stings: From epidemiology to treatment. Allergy 2011, 66, 35–37. [Google Scholar] [CrossRef]
  7. Pühringer, V.; Jilma, B.; Herkner, H. Population-based incidence of all-cause anaphylaxis and its development over time: A systematic review and meta-analysis. Front. Allergy 2023, 4, 1249280. [Google Scholar] [CrossRef] [PubMed]
  8. Worm, M.; Monaret-Vautrin, A.; Scherer, K.; Lang, R. First European data from the network of severe allergic reactions (NORA). Eur. J. Allergy Clin. Immunol. 2014, 69, 1397–1404. [Google Scholar] [CrossRef] [PubMed]
  9. Rijavec, M.; Inkret, J.; Bidovec-stojkovi, U. Fatal Hymenoptera Venom—Triggered Anaphylaxis in Patients with Unrecognized Clonal Mast Cell Disorder—Is Mastocytosis to Blame? Int. J. Mol. Sci. 2023, 24, 16368. [Google Scholar] [CrossRef] [PubMed]
  10. Ruëff, F.; Vos, B.; Elberink, J.O.; Bender, A.; Chatelain, R.; Dugas-Breit, S.; Horny, H.; Küchenhoff, H.; Linhardt, A.; Mastnik, S.; et al. Predictors of clinical effectiveness of Hymenoptera venom immunotherapy. Clin. Exp. Allergy 2014, 44, 736–746. [Google Scholar] [CrossRef]
  11. Lang, R.; Hawranek, T. Hymenoptera Venom Immunotherapy and Field Stings. J. Investig. Allergol. Clin. Immunol. 2006, 16, 224. [Google Scholar]
  12. Jarkvist, J.; Salehi, C.; Akin, C.; Gülen, T. Venom immunotherapy in patients with clonal mast cell disorders: IgG4 correlates with protection. Allergy Eur. J. Allergy Clin. Immunol. 2020, 75, 169–177. [Google Scholar] [CrossRef]
  13. Kačar, M.; Rijavec, M.; Šelb, J.; Korošec, P. Clonal mast cell disorders and hereditary α-tryptasemia as risk factors for anaphylaxis. Clin. Exp. Allergy 2023, 53, 392–404. [Google Scholar] [CrossRef]
  14. Bonadonna, P.; Zanotti, R.; Pagani, M.; Bonifacio, M.; Scaffidi, L.; Olivieri, E.; Franchini, M.; Reccardini, F.; Costantino, M.T.; Roncallo, C.; et al. Anaphylactic Reactions After Discontinuation of Hymenoptera Venom Immunotherapy: A Clonal Mast Cell Disorder Should Be Suspected. J. Allergy Clin. Immunol. Pract. 2018, 6, 1368–1372. [Google Scholar] [CrossRef] [PubMed]
  15. Shamji, M.H.; Durham, S.R. Mechanisms of allergen immunotherapy for inhaled allergens and predictive biomarkers. J. Allergy Clin. Immunol. 2017, 140, 1485–1498. [Google Scholar] [CrossRef] [PubMed]
  16. Möbs, C.; Müller, J.; Rudzio, A.; Pickert, J.; Blank, S.; Jakob, T.; Spillner, E.P.W. Decline of Ves v 5-specific blocking capacity in wasp venom-allergic patients after stopping allergen immunotherapy. Allergy 2015, 70, 715–719. [Google Scholar] [CrossRef] [PubMed]
  17. Eržen, R.; Košnik, M.; Silar, M.K.P. Basophil response and the induction of a tolerance in venom immunotherapy: A long-term sting challenge study. Allergy 2012, 67, 822–830. [Google Scholar] [CrossRef] [PubMed]
  18. Žitnik, S.E.; Vesel, T.; Avčin, T.; Šilar, M.; Košnik, M.; Korošec, P. Monitoring honeybee venom immunotherapy in children with the basophil activation test. Pediatr. Allergy Immunol. 2012, 23, 166–172. [Google Scholar] [CrossRef]
  19. Čelesnik, N.; Vesel, T.; Rijavec, M.; Šilar, M.; Eržen, R.; Košnik, M.; Žitnik, S.E.K.; Avčin, T.; Korošec, P. Short-term venom immunotherapy induces desensitization of FcεRI-mediated basophil response. Allergy Eur. J. Allergy Clin. Immunol. 2012, 67, 1594–1600. [Google Scholar] [CrossRef]
  20. Foschi, F.G.; Emiliani, F.; Savini, S.; Quercia, O.; Stefanini, G.F. CD30 serum levels and response to Hymenoptera venom immunotherapy. J. Investig. Allergol. Clin. Immunol. 2008, 18, 279–283. [Google Scholar]
  21. Packi, K.; Matysiak, J.; Matuszewska, E.; Bręborowicz, A.; Kycler, Z.; Matysiak, J. New biomarkers of hymenoptera venom allergy in a group of inflammation factors. Int. J. Environ. Res. Public Health 2021, 18, 4011. [Google Scholar] [CrossRef]
  22. Frick, M.; Fischer, J.; Helbling, A.; Ruëff, F.; Wieczorek, D.; Ollert, M.; Pfützner, W.; Müller, S.; Huss-Marp, J.; Dorn, B.; et al. Predominant Api m 10 sensitization as risk factor for treatment failure in honey bee venom immunotherapy. J. Allergy Clin. Immunol. 2016, 138, 1663–1671.e9. [Google Scholar] [CrossRef]
  23. Jakob, T.; Rauber, M.M.; Perez-Riverol, A.; Spillner, E.; Blank, S. The Honeybee Venom Major Allergen Api m 10 (Icarapin) and Its Role in Diagnostics and Treatment of Hymenoptera Venom Allergy. Curr. Allergy Asthma Rep. 2020, 20, 48. [Google Scholar] [CrossRef] [PubMed]
  24. Blank, S.; Jakwerth, C.A.; Zissler, U.M.; Schmidt-Weber, C.B. Molecular determination of insect venom allergies. Expert Rev. Mol. Diagn. 2022, 22, 983–986. [Google Scholar] [CrossRef] [PubMed]
  25. Demšar Luzar, A.; Korošec, P.; Košnik, M.; Zidarn, M.; Rijavec, M. Hymenoptera venom immunotherapy: Immune mechanisms of induced protection and tolerance. Cells 2021, 10, 1575. [Google Scholar] [CrossRef] [PubMed]
  26. Tsicopoulos, A. Venom immunotherapy modulates interleukin-4 and interferon-γ messenger RNA expression of peripheral T lymphocytes. Immunology 1996, 87, 593–598. [Google Scholar] [CrossRef]
  27. Konno, S.; Hizawa, N.; Nishimura, M.; Huang, S.K. Osteopontin: A potential biomarker for successful bee venom immunotherapy and a potential molecule for inhibiting IgE-mediated allergic responses. Allergol. Int. 2006, 55, 355–359. [Google Scholar] [CrossRef]
  28. Konno, S.; Golden, D.B.K.; Schroeder, J.; Hamilton, R.G.; Lichtenstein, L.M.; Huang, S.K. Increased expression of osteopontin is associated with long-term bee venom immunotherapy. J. Allergy Clin. Immunol. 2005, 115, 1063–1067. [Google Scholar] [CrossRef]
  29. Niedoszytko, M.; Bruinenberg, M.; de Monchy, J.; Wijmenga, C.; Platteel, M.; Jassem, E.; Elberink, J.N.O. Gene expression analysis in predicting the effectiveness of insect venom immunotherapy. J. Allergy Clin. Immunol. 2010, 125, 1092–1097. [Google Scholar] [CrossRef]
  30. Niedoszytko, M.; Bruinenberg, M.; de Monchy, J.; Weersma, R.K.; Wijmenga, C.; Jassem, E.; Elberink, J.N.O. Changes in gene expression caused by insect venom immunotherapy responsible for the long-term protection of insect venom-allergic patients. Ann. Allergy Asthma Immunol. 2011, 106, 502–510. [Google Scholar] [CrossRef] [PubMed]
  31. Karpinski, P.; Skiba, P.; Kosinska, M.; Rosiek-Biegus, M.; Królewicz, E.; Blin, N.; Meese, E.; Panaszek, B.; Nittner-Marszalska, M.; Sasiadek, M.M. Genome-wide analysis of gene expression after one year of venom immunotherapy. Immunol. Lett. 2018, 204, 23–28. [Google Scholar] [CrossRef]
  32. McKenzie, C.I.; Varese, N.; Aui, P.M.; Reinwald, S.; Wines, B.D.; Hogarth, P.M.; Thien, F.; Hew, M.; Rolland, J.M.; O’Hehir, R.E.; et al. RNA sequencing of single allergen-specific memory B cells after grass pollen immunotherapy: Two unique cell fates and CD29 as a biomarker for treatment effect. Allergy Eur. J. Allergy Clin. Immunol. 2023, 78, 822–835. [Google Scholar] [CrossRef]
  33. Chao, W.; Liao, W.; Wang, J.; Sung, H.; Chen, H.; Lin, F.; Huang, J.; Liu, K.; Ko, T. Single-cell profiling reveals novel cellular heterogeneity of monocytes during Hymenoptera venom allergy. Clin. Transl. Allergy 2022, 12, e12151. [Google Scholar] [CrossRef]
  34. Layhadi, J.A.; Moya, R.; Tan, T.J.; Lenormand, M.M.; Sharif, H.; Parkin, R.V.; Vila-Nadal, G.; Fedina, O.; Zhu, R.; Laisuan, W.; et al. Single-cell RNA sequencing identifies precise tolerogenic cellular and molecular pathways induced by depigmented-polymerized grass pollen allergen extract. J. Allergy Clin. Immunol. 2023, 151, 1357–1370.e9. [Google Scholar] [CrossRef] [PubMed]
  35. Alberts, B.; Johnson, A.; Lewis, J.; Raff, M.; Roberts, K.; Walter, P. Table 22-1 Blood Cells. In Molecular Biology of the Cell, 4th ed.; Garland Science: New York, NY, USA, 2002. [Google Scholar]
  36. Sivaprasad, U.; Kinker, K.G.; Ericksen, M.B.; Lindsey, M.; Gibson, A.M.; Bass, S.A.; Hershey, N.S.; Deng, J.; Medvedovic, M.; Hershey, G.K.K. SERPINB3/B4 Contributes to Early Inflammation and Barrier Dysfunction in an Experimental Murine Model of Atopic Dermatitis. J. Investig. Dermatol. 2015, 135, 160–169. [Google Scholar] [CrossRef] [PubMed]
  37. Howell, W.M.; Carter, V.; Clark, B. The HLA system: Immunobiology, HLA typing, antibody screening and crossmatching techniques. J. Clin. Pathol. 2010, 63, 387–390. [Google Scholar] [CrossRef]
  38. Gene Cards®. The Human Gene Database n.d. Available online: https://www.genecards.org/ (accessed on 23 January 2024).
  39. Mocciaro, E.; Runfola, V.; Ghezzi, P.; Pannese, M.; Gabellini, D. DUX4 role in normal physiology and in FSHD muscular dystrophy. Cells 2021, 10, 3322. [Google Scholar] [CrossRef] [PubMed]
  40. Dmitriev, P.; Kiseleva, E.; Kharchenko, O.; Ivashkin, E.; Pichugin, A.; Dessen, P.; Robert, T.; Coppée, F.; Belayew, A.; Carnac, G.; et al. Dux4 controls migration of mesenchymal stem cells through the Cxcr4-Sdf1 axis. Oncotarget 2016, 7, 65090–65108. [Google Scholar] [CrossRef]
  41. Liu, T.; Zhang, L.; Joo, D.; Sun, S.C. NF-κB signaling in inflammation. Signal Transduct. Target. Ther. 2017, 2, 17023. [Google Scholar] [CrossRef]
  42. Hayden, M.S.; Ghosh, S. NF-κB, the first quarter-century: Remarkable progress and outstanding questions. Genes Dev. 2012, 26, 203–234. [Google Scholar] [CrossRef]
  43. Gustavsson, E.K.; Sethi, S.; Gao, Y.; Brenton, J.W.; García-Ruiz, S.; Zhang, D.; Garza, R.; Reynolds, R.H.; Evans, J.R.; Chen, Z.; et al. The annotation and function of the Parkinson’s and Gaucher disease-linked gene GBA1 has been concealed by its protein-coding pseudogene GBAP1. BioRxiv 2023. [Google Scholar] [CrossRef]
  44. Pink, R.C.; Wicks, K.; Caley, D.P.; Punch, E.K.; Jacobs, L.; Carter, D.R.F. Pseudogenes: Pseudo-functional or key regulators in health and disease? RNA 2011, 17, 792–798. [Google Scholar] [CrossRef]
  45. Tan, L.; Qin, Y.; Xie, R.; Xia, T.; Duan, X.; Peng, L.; You, R.; Liu, Y.; Zou, X.; Zhang, M.; et al. Theranostics N6-methyladenosine-associated prognostic pseudogenes contribute to predicting immunotherapy benefits and therapeutic agents in head and neck squamous cell carcinoma. Theranostics 2022, 12, 7267–7288. [Google Scholar] [CrossRef] [PubMed]
  46. Kempiński, K.; Romantowski, J.; Maciejewska, A.; Pawłowski, R.; Chełmińska, M.; Jassem, E.; Niedoszytko, M. COMMD8 changes expression during initial phase of wasp venom immunotherapy. J. Gene Med. 2020, 22, e3243. [Google Scholar] [CrossRef] [PubMed]
  47. Starokadomskyy, P.; Gluck, N.; Li, H.; Chen, B.; Wallis, M.; Maine, G.N.; Mao, X.; Zaidi, I.W.; Hein, M.Y.; McDonald, F.J.; et al. CCDC22 deficiency in humans blunts activation of proinfammatory NF-κB signaling. J. Clin. Investig. 2013, 123, 2244–2256. [Google Scholar] [CrossRef] [PubMed]
  48. Ande, S.R.; Xu, Y.X.Z.; Mishra, S. Prohibitin: A potential therapeutic target in tyrosine kinase signaling. Signal Transduct. Target. Ther. 2017, 2, 17059. [Google Scholar] [CrossRef] [PubMed]
  49. Werner, L.M.; Criss, A.K. Diverse Functions of C4b-Binding Protein in Health and Disease. J. Immunol. 2023, 211, 1443–1449. [Google Scholar] [CrossRef] [PubMed]
  50. Sjef Verbeek, J.; Hirose, S.; Nishimura, H. The complex association of fcγriib with autoimmune susceptibility. Front. Immunol. 2019, 10, 2061. [Google Scholar] [CrossRef]
  51. Phillips, N.E.; Parker, D.C. Cross-linking of B lymphocyte Fc gamma receptors and membrane immunoglobulin inhibits anti-immunoglobulin-induced blastogenesis. J. Immunol. 1984, 132, 627–632. [Google Scholar] [CrossRef]
  52. Varma, R.; Puri, N. Dampening of mast cell secondary responses to allergen involves specific signalling and epigenetic changes. Cell Immunol. 2019, 344, 103944. [Google Scholar] [CrossRef]
  53. Odegaard, J.I.; Chawla, A. Alternative Macrophage Activation and Metabolism. Annu. Rev. Pathol. 2011, 6, 275–297. [Google Scholar] [CrossRef]
  54. Nobs, S.P.; Natali, S.; Pohlmeier, L.; Okreglicka, K.; Schneider, C.; Kurrer, M.; Sallusto, F.; Kopf, M. PPARγ in dendritic cells and T cells drives pathogenic type-2 effector responses in lung inflammation. J. Exp. Med. 2017, 214, 3015–3035. [Google Scholar] [CrossRef]
  55. Chen, T.; Tibbitt, C.A.; Feng, X.; Stark, J.M.; Rohrbeck, L.; Rausch, L.; Sedimbi, S.K.; Karlsson, M.C.I.; Lambrecht, B.N.; Hedestam, G.B.K.; et al. PPAR-promotes type 2 immune responses in allergy and nematode infection. Sci. Immunol. 2017, 2, 1–11. [Google Scholar] [CrossRef] [PubMed]
  56. Micossé, C.; von Meyenn, L.; Steck, O.; Kipfer, E.; Adam, C.; Simillion, C.; Jafari, S.M.S.; Olah, P.; Yawlkar, N.; Simon, D.; et al. Human “TH9” cells are a subpopulation of PPAR-+ TH2 cells. Sci. Immunol. 2019, 4, 1–14. [Google Scholar] [CrossRef]
  57. McHugh, S.M.; Deighton, J.; Stewart, A.G.; Lachmann, P.J.; Ewan, P.W. Bee venom immunotherapy induces a shift in cytokine responses from a TH-2 to a TH-1 dominant pattern: Comparison of rush and conventional immunotherapy. Clin. Exp. Allergy 1995, 25, 828–838. [Google Scholar] [CrossRef] [PubMed]
  58. Huang, J.; Zhang, W.; Xiang, R.; Tan, L.; Liu, P.; Tao, Z.; Deng, Y.; Tong, H.; Xu, Y. The early-phase transcriptome and the clinical efficacy analysis in three modes of subcutaneous immunotherapy for allergic rhinitis. World Allergy Organ. J. 2023, 16, 100811. [Google Scholar] [CrossRef] [PubMed]
  59. Karisola, P.; Palosuo, K.; Hinkkanen, V.; Wisgrill, L.; Savinko, T.; Fyhrquist, N.; Alenius, H.; Mäkelä, M.J. Integrative Transcriptomics Reveals Activation of Innate Immune Responses and Inhibition of Inflammation during Oral Immunotherapy for Egg Allergy in Children. Front. Immunol. 2021, 12, 704633. [Google Scholar] [CrossRef]
  60. Gruzelle, V.; Mailhol, C.; Waters, D.W.; Guilleminault, L. Clinical utility of rush venom immunotherapy: Current status. J. Asthma Allergy 2020, 13, 1–10. [Google Scholar] [CrossRef]
  61. Pospischil, I.M.; Kagerer, M.; Cozzio, A.; Angelova-Fischer, I.; Guenova, E.; Ballmer-Weber, B.; Hoetzenecker, W. Comparison of the Safety Profiles of 3 Different Hymenoptera Venom Immunotherapy Protocols: A Retrospective 2-Center Study of 143 Patients. Int. Arch. Allergy Immunol. 2020, 181, 783–789. [Google Scholar] [CrossRef]
  62. CLC Genomics Workbench 23.0.2 Online Manual. Available online: https://resources.qiagenbioinformatics.com/manuals/clcgenomicsworkbench/2302/index.php?manual=Introduction_CLC_Genomics_Workbench.html (accessed on 14 December 2023).
  63. Aken, B.L.; Ayling, S.; Barrell, D.; Clarke, L.; Curwen, V.; Fairley, S.; Banet, J.F.; Billis, K.; Girón, C.G.; Hourlier, T.; et al. The Ensembl gene annotation system. Database 2016, 2016, baw093. [Google Scholar] [CrossRef]
Ijms 25 03499 g001
Figure 1. Follow-up analysis and overlapping of differentially expressed genes at three time points (before, after reaching maintenance dose, and after) for successful venom immunotherapy (blue) and at two time points (before and after) for treatment failure (red). (Cut-off: fold change > 5, p-value < 0.05.) VIT: venom immunotherapy; DEGs: differentially expressed genes.
Figure 1. Follow-up analysis and overlapping of differentially expressed genes at three time points (before, after reaching maintenance dose, and after) for successful venom immunotherapy (blue) and at two time points (before and after) for treatment failure (red). (Cut-off: fold change > 5, p-value < 0.05.) VIT: venom immunotherapy; DEGs: differentially expressed genes.
Ijms 25 03499 g001
Ijms 25 03499 g002
Figure 2. A volcano plot of differentially expressed genes between successful venom immunotherapy and treatment failure groups. Highlighted and with red arrows marked is the transcript of interest, IKBKGP1. No difference between the two groups before the initiation of venom immunotherapy is observed (a), whereas a significant difference between the two groups after venom immunotherapy is pointed out (b). The IKBKGP1 transcript is significantly downregulated after venom immunotherapy in the group with successful treatment. The kinetics of the IKBKGP1 transcripts of individual samples from successful venom immunotherapy are shown, with significant downregulation after venom immunotherapy (p-value = 0.0156) (c). RPKM: reads per kilobase per million mapped reads; *: statistical significance.
Figure 2. A volcano plot of differentially expressed genes between successful venom immunotherapy and treatment failure groups. Highlighted and with red arrows marked is the transcript of interest, IKBKGP1. No difference between the two groups before the initiation of venom immunotherapy is observed (a), whereas a significant difference between the two groups after venom immunotherapy is pointed out (b). The IKBKGP1 transcript is significantly downregulated after venom immunotherapy in the group with successful treatment. The kinetics of the IKBKGP1 transcripts of individual samples from successful venom immunotherapy are shown, with significant downregulation after venom immunotherapy (p-value = 0.0156) (c). RPKM: reads per kilobase per million mapped reads; *: statistical significance.
Ijms 25 03499 g002
Ijms 25 03499 g003aIjms 25 03499 g003b
Figure 3. Comparison of differentially expressed genes between the three groups—successful VIT, treatment failure, and healthy controls. Overlapping genes specific for successful VIT with the two transcripts characteristic of successful VIT, regardless of the time point of sampling (From: Gene Cards®. The Human Gene Database n.d. Available online: https://www.genecards.org/, accessed on 23 January 2024) (a). Overlapping genes specific for treatment failure (b). Number of DEGs after VIT in the analysis of successful VIT and healthy controls compared to treatment failure and healthy controls (c). (Cut-off: fold change > 5, p-value < 0.05.) VIT: venom immunotherapy; DEGs: differentially expressed genes.
Figure 3. Comparison of differentially expressed genes between the three groups—successful VIT, treatment failure, and healthy controls. Overlapping genes specific for successful VIT with the two transcripts characteristic of successful VIT, regardless of the time point of sampling (From: Gene Cards®. The Human Gene Database n.d. Available online: https://www.genecards.org/, accessed on 23 January 2024) (a). Overlapping genes specific for treatment failure (b). Number of DEGs after VIT in the analysis of successful VIT and healthy controls compared to treatment failure and healthy controls (c). (Cut-off: fold change > 5, p-value < 0.05.) VIT: venom immunotherapy; DEGs: differentially expressed genes.
Ijms 25 03499 g003aIjms 25 03499 g003b
Ijms 25 03499 g004
Figure 4. A heat map of normalised expression values of differentially expressed genes after the sting challenge in the groups with successful venom immunotherapy (dark blue) and treatment failure (dark red). Each column represents a patient’s sample; each row represents an individual gene. For each gene, bright blue represents under-expression, and bright red over-expression. (Cut-off: fold change > 5, p-value < 0.05.) VIT: venom immunotherapy.
Figure 4. A heat map of normalised expression values of differentially expressed genes after the sting challenge in the groups with successful venom immunotherapy (dark blue) and treatment failure (dark red). Each column represents a patient’s sample; each row represents an individual gene. For each gene, bright blue represents under-expression, and bright red over-expression. (Cut-off: fold change > 5, p-value < 0.05.) VIT: venom immunotherapy.
Ijms 25 03499 g004
Ijms 25 03499 g005
Figure 5. The follow-up study design. Whole-blood samples were collected before the beginning of VIT, after reaching the maintenance dose, one year after finishing VIT, and after a controlled sting challenge. VIT: venom immunotherapy.
Figure 5. The follow-up study design. Whole-blood samples were collected before the beginning of VIT, after reaching the maintenance dose, one year after finishing VIT, and after a controlled sting challenge. VIT: venom immunotherapy.
Ijms 25 03499 g005
Table 1. Hub genes differentially expressed after reaching maintenance dose. Positive fold change represents genes upregulated after reaching maintenance dose. Negative fold change represents genes downregulated after reaching maintenance dose.
Table 1. Hub genes differentially expressed after reaching maintenance dose. Positive fold change represents genes upregulated after reaching maintenance dose. Negative fold change represents genes downregulated after reaching maintenance dose.
GeneOntology [38]Maintenance Dose vs. BeforeMaintenance Dose vs. After
Fold Changep-ValueFold Changep-Value
IGLV2–5Immunoglobulin Lambda Variable 2–56.40.046.50.03
HLA-UMajor Histocompatibility Complex, Class I, U 6.50.046.60.03
SERPINB4May act as a protease inhibitor to modulate the host immune response against tumour cells7.20.027.70.01
Table 2. Hub genes expressed after reaching maintenance dose and one year after finishing venom immunotherapy. Positive fold change represents genes upregulated after reaching maintenance dose and after venom immunotherapy. Negative fold change represents genes downregulated after reaching maintenance dose and after venom immunotherapy.
Table 2. Hub genes expressed after reaching maintenance dose and one year after finishing venom immunotherapy. Positive fold change represents genes upregulated after reaching maintenance dose and after venom immunotherapy. Negative fold change represents genes downregulated after reaching maintenance dose and after venom immunotherapy.
GeneOntology [38]Maintenance Dose vs. BeforeAfter vs. Before
Fold Changep-ValueFold Changep-Value
DUX4L10Involvement in inflammation, cellular migration, and chemotaxis−240.16.46 × 10−7−7.30.0009
DUX4L12Involvement in inflammation, cellular migration, and chemotaxis−90.15.48 × 10−5−6.50.002
CFL1P2Required for the upregulation of the atypical chemokine receptor ACKR2 −9.90.02−8.30.02
HTN3They function as antimicrobial peptides and are important components of the innate immune system−7.60.009−11.70.002
CDH9Cadherins are calcium-dependent cell adhesion proteins−7.00.008−5.60.01
PCDHB8Calcium-dependent cell adhesion protein involved in cell self-recognition and non-self discrimination−6.90.01−8.80.002
USP17L22Regulates different cellular processes, which may include cell proliferation, progression through the cell cycle, apoptosis, cell migration, and the cellular response to viral infection−6.50.01−130.004
IKBKGP1The NF-kappa-B essential modulator (also known as NEMO) regulating the NFкB pathway, which is involved in many immune and inflammatory responses −5.50.01−24.90.00001
MAPK15Regulates several processes, such as autophagy, protein trafficking, and genome integrity, in a kinase activity-dependent manner−50.03−5.60.01
Table 3. Hub genes in the group of patients who tolerated a sting (successful venom immunotherapy; after vs. before). A positive fold change represents genes upregulated after venom immunotherapy. A negative fold change represents genes downregulated after venom immunotherapy.
Table 3. Hub genes in the group of patients who tolerated a sting (successful venom immunotherapy; after vs. before). A positive fold change represents genes upregulated after venom immunotherapy. A negative fold change represents genes downregulated after venom immunotherapy.
GeneOntology [38]Fold Changep-Value
PHBP6Interacts with STAT3 to affect IL17 secretion in T-helper Th17 cells−240.001
DUX4Involvement in inflammation, cellular migration, and chemotaxis−18.50.0007
ZG16Predicted to act upstream of or in defence response to Gram-positive bacteria−120.007
CACNG5Involved in TCR signalling−10.40.01
CPN1Among its related pathways are the complement cascade and innate immune system−100.011
TMEM196Predicted to be an integral component of the membrane. Diseases associated with TMEM196 include Adult-Onset Severe Asthma−90.002
IGSF21Proteins in this superfamily are usually found on or in cell membranes and act as receptors in immune response pathways−90.007
DUX4L14Involvement in inflammation, cellular migration, and chemotaxis−8.30.03
TDO2This gene encodes a heme enzyme that plays a critical role in tryptophan metabolism−80.004
MMD2Monocyte To Macrophage Differentiation Associated 2; annotations related to this gene include protein kinase activity−80.026
TWIST1Represses expression of pro-inflammatory cytokines such as TNFA and IL1B−70.014
DUX4L2Involvement in inflammation, cellular migration, and chemotaxis−6.50.01
UMODMay serve as a receptor for the binding and endocytosis of cytokines (IL-1, IL-2) and TNF−60.043
DPTDermatopontin is postulated to modify the behaviour of TGF-beta through interaction with decorin. Enhances TGFB1 activity.−60.014
SLC17A2Predicted to be located in the membrane. Predicted to be active in the lysosome. −60.018
GNB3Related to chronic rhinosinusitis.−60.006
PLA2G10This gene encodes a member of the phospholipase A2 family of proteins. This enzyme plays a role in the production of various inflammatory lipid mediators, such as prostaglandins.−50.016
NLRP8NLRP genes play roles in the mammalian innate immune system through inflammasome formation and the activation of caspases−50.042
EBI3It encodes a secreted glycoprotein and forms interleukin 27 (IL-27). IL-27 partially regulates T-cell and inflammatory responses by activating the Jak/STAT pathway of CD4+ T cells. Among its related pathways are IL27-mediated signalling events and interleukin-12 family signalling−50.040
PCDHA9Potential calcium-dependent cell adhesion protein. It may be involved in establishing and maintaining specific neuronal connections in the brain−50.013
IGKV2D-28The V region of the variable domain of immunoglobulin light chains that participate in the antigen recognition60.002
E2F4P1It plays a vital role in suppressing proliferation-associated genes, and its gene mutation and increased expression may be associated with human cancer70.042
Table 4. Study group characteristics.
Table 4. Study group characteristics.
Patients on VITHealthy Controlsp-Value
Successful VITUnsuccessful VIT
No. of patients1459/
Age [years], median (range)50 (27–63)49 (26–59)47 (24–58)ns 1
Sex, no. of patients (%)
Female
Male
3 (21%)
11 (79%)
0
5 (100%)
2 (22%)
7 (78%)		ns
Specific honeybee IgE [kIU/L], median (range) 2 5.48 (0.45–52.40)3.46 (0.82–11.4)/ns
Double sensitisation; honeybee and wasp sIgE, no. of patients (%) 210 (71%)1 (20%)/ns
Total IgE [IU/mL], median (range) 2113 (3–661)32 (8–297)/ns
Years of VIT, median (range)6 (6–7)7 (6–7)/ns
Complications during VIT, n (%)1 (7%)5 (100%)/0.0005
1 ns, not significant; 2 measurements of specific IgE and total IgE were made at the beginning of venom immunotherapy.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Demšar Luzar, A.; Korošec, P.; Košnik, M.; Zidarn, M.; Rijavec, M. Blood Transcriptomics Identifies Multiple Gene Expression Pathways Associated with the Clinical Efficacy of Hymenoptera Venom Immunotherapy. Int. J. Mol. Sci. 2024, 25, 3499. https://doi.org/10.3390/ijms25063499

AMA Style

Demšar Luzar A, Korošec P, Košnik M, Zidarn M, Rijavec M. Blood Transcriptomics Identifies Multiple Gene Expression Pathways Associated with the Clinical Efficacy of Hymenoptera Venom Immunotherapy. International Journal of Molecular Sciences. 2024; 25(6):3499. https://doi.org/10.3390/ijms25063499

Chicago/Turabian Style

Demšar Luzar, Ajda, Peter Korošec, Mitja Košnik, Mihaela Zidarn, and Matija Rijavec. 2024. "Blood Transcriptomics Identifies Multiple Gene Expression Pathways Associated with the Clinical Efficacy of Hymenoptera Venom Immunotherapy" International Journal of Molecular Sciences 25, no. 6: 3499. https://doi.org/10.3390/ijms25063499

APA Style

Demšar Luzar, A., Korošec, P., Košnik, M., Zidarn, M., & Rijavec, M. (2024). Blood Transcriptomics Identifies Multiple Gene Expression Pathways Associated with the Clinical Efficacy of Hymenoptera Venom Immunotherapy. International Journal of Molecular Sciences, 25(6), 3499. https://doi.org/10.3390/ijms25063499

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop