Skip to main content
Download PDF

Identification of plasma protein markers of allergic disease risk: a mendelian randomization approach to proteomic analysis

BMC Genomics volume 25, Article number: 503 (2024) Cite this article

Abstract

Background

While numerous allergy-related biomarkers and targeted treatment strategies have been developed and employed, there are still signifcant limitations and challenges in the early diagnosis and targeted treatment for allegic diseases. Our study aims to identify circulating proteins causally associated with allergic disease-related traits through Mendelian randomization (MR)-based analytical framework.

Methods

Large-scale cis-MR was employed to estimate the effects of thousands of plasma proteins on five main allergic diseases. Additional analyses including MR Steiger analyzing and Bayesian colocalisation, were performed to test the robustness of the associations; These findings were further validated utilizing meta-analytical methods in the replication analysis. Both proteome- and transcriptome-wide association studies approach was applied, and then, a protein-protein interaction was conducted to examine the interplay between the identified proteins and the targets of existing medications.

Results

Eleven plasma proteins were identified with links to atopic asthma (AA), atopic dermatitis (AD), and allergic rhinitis (AR). Subsequently, these proteins were classified into four distinct target groups, with a focus on tier 1 and 2 targets due to their higher potential to become drug targets. MR analysis and extra validation revealed STAT6 and TNFRSF6B to be Tier 1 and IL1RL2 and IL6R to be Tier 2 proteins with the potential for AA treatment. Two Tier 1 proteins, CRAT and TNFRSF6B, and five Tier 2 proteins, ERBB3, IL6R, MMP12, ICAM1, and IL1RL2, were linked to AD, and three Tier 2 proteins, MANF, STAT6, and TNFSF8, to AR.

Conclusion

Eleven Tier 1 and 2 protein targets that are promising drug target candidates were identified for AA, AD, and AR, which influence the development of allergic diseases and expose new diagnostic and therapeutic targets.

Peer Review reports

Introduction

Allergic diseases arise from inappropriate initiation of type 2 immune responses against innocuous environmental antigens and include disorders such as atopic dermatitis (AD), allergic asthma (AA), allergic rhinitis (AR), allergic conjunctivitis (AC) and allergic urticaria (AU) [1, 2]. The development of AD in early life, followed by other allergies, such as asthma, has been described as the atopic march [3, 4]. Allergic disease incidence has increased over the past 3 decades to affect an estimated 433 million people worldwide and exert a considerable economic and social burden [5, 6]. A combination of genetics, exposure to environmental allergens and irritants, microbial interactions, and abnormal immune responses contribute to inflammation and the atopic march [7, 8]. However, allergic disease pathogenesis is complex and poorly understood with limited drug options to tackle the recurrent and potentially life-long nature. These observations illustrate the need for further mechanistic studies and identification of drug targets.

The human plasma proteome participates in inter-tissue communications via metabolic, signaling, and physiological pathways and includes potential drug targets [9,10,11,12,13]. The link between levels of specific plasma proteins and allergic disease risk has been reported previously but residual confounders and reverse causation mean that observational studies do not demonstrate causality [14, 15]. Mendelian randomization (MR) uses genetic variation as an instrumental variable (IV) to explore the causal effect of exposure on outcomes. The impact of genes or proteins is difficult to interpret from isolated GWAS outcomes due to the occurrence of single nucleotide polymorphisms (SNPs). However, proteome-wide association studies (PWAS) give a clearer picture of the influence of proteins on disease [16]. In addition, functional summary-based imputation software allows transcriptome-wide association studies (TWAS) to assess the impact of whole blood protein-coding gene expression on allergic disease risk [17].

The current study used a PWAS approach, combining GWAS data with protein quantitative trait locus (pQTL) and TWAS with expression quantitative trait loci (eQTL), to investigate potential proteinaceous diagnostic and therapeutic targets in allergic disease. Bayesian colocalization analysis and phenotype scanning were conducted to verify causality between candidate proteins and disease pathogenesis. A protein-protein interaction (PPI) network was constructed and pathway enrichment was performed to illuminate potential mechanisms. The detailed research flowchart is presented in Fig. 1.

Fig. 1
figure 1

Flowchart of evidence-based grading of potential drug targets from plasma protein on allergic diseases. AA: Allergic asthma; AC: Allergic conjunctivitis; AD: Atopic dermatitis; AR: Allergic rhinitis; AU: Allergic urticaria; PPI: Protein–protein interaction; pQTL: Protein quantitative trait loci; PPH4: posterior probability of hypothesis 4; PWAS: proteome-wide association studies; TWAS: transcriptome-wide association studies

Full size image

Materials and methods

Data sources and selection of IVs

Whole blood pQTL data concerning 4,657 proteins derived from 7,213 participants of European descent in the ARIC study was analyzed [18]. SomaScan technology on the v.4.1 platform (SomaLogic) was used for proteomic profiling. pQTLs were filtered according to the following criteria: (1) variants had cis-acting effects, defined as being situated within a 1 Mb region up- or downstream of the gene encoding the plasma protein; (2) pQTLs met the genome-wide significance threshold of p < 5 × 10 − 8; (3) no significant linkage disequilibrium (LD) was present among pQTLs (r2 < 0.001); (4) pQTLs were located outside the major histocompatibility complex (MHC) region (chr6, 26-34 Mb). The analysis was repeated using pQTL data from 35,559 Icelandic individuals, including 4,907 plasma proteins (Icelandic Cancer Project and deCODE genetics, Reykjavík, Iceland) on SomaScan v 4.0 platform [19].

Data for individual allergic diseases, including 38,369 cases and 411,131 controls for AA, 613 cases and 464,657 controls for AC, and 1,057 cases and 482,892 controls for AU were obtained from a UKB cohort GWAS study [20]. A UKB cohort from a trans-biobank meta-analysis study [21] yielded data for 22,474 cases of AD and 774,187 controls. Data on 27,415 cases of AR and 457,183 controls came from a separate GWAS study [22]. Summary statistics for the allergic diseases under scrutiny were obtained from the FinnGen Biobank Analysis Consortium database (https://finngen.gitbook.io/documentation/). Allergic diseases were diagnosed according to ICD-10 (International Classification of Diseases) criteria. Data sources are presented in detail in S. Table 1.

Statistical analysis

Mendelian randomization analysis

MR analysis was performed using R software (version 4.1.3) and the TwoSampleMR package (version 0.5.6) [23]. Effect estimates of plasma proteins with a single SNP were generated using the Wald ratio [24] and of those with two or more SNPs by the Inverse Variance Weighted (IVW) method [25]. A false discovery rate (FDR, p < 0.05) correction was applied to account for multiple comparisons in validating MR results [26]. Results are expressed as odds ratios per standard deviation increase in genetically determined plasma protein levels. Considering the genetic representativeness of the results, the ARIC/UKB combination, which includes a mixed European genetic background, was chose as the discovery analysis.

Sensitivity analysis

Reverse causality between the proteins identified by discovery analysis and allergic diseases was evaluated by bidirectional MR analysis and MR Steiger test [26]. Four additional methods were employed in bidirectional MR analysis to assess the IV validity focusing on the three IV principles: strong association with exposure, direct influence on the outcome through exposure, and no association with outcome confounders. Violations of these principles could compromise the accuracy of the results. MR-Egger regression slope and intercept were used to estimate pleiotropy across IVs to give an adjusted estimate independent of IV validity [27, 28]. MR-PRESSO was used to identify outliers causing significant pleiotropy and heterogeneity to give a corrected causal effect assessment [29]. The weighted-median method enabled consistent inference even when over 50% of IVs were valid [30]. MR-Robust Adjusted Profile Score (MRAPS) increased statistical power to give robust estimates in the presence of weak instrumental bias and horizontal pleiotropy [31]. By default, IVW results are preferred [32], but we turn to MR-Egger when significant pleiotropy is detected by the MR-Egger pleiotropy test. If the MR-PRESSO global test identifies significant outliers, we prioritize results corrected by MR-PRESSO. A Bonferroni correction was applied to address multiple comparison errors and a corrected value of p < 0.01 was considered to indicate significance in reverse MR. Plasma proteins that appeared significant during both MR Steiger and reverse MR analysis, suggestive of reverse causality, were excluded from further analysis.

Potential links between all identified proteins and confounders were investigated via phenotype scanning using the Phenoscanner database [33] with a genome-wide significance threshold of p < 5 × 10−8. pQTLs linked to known allergic disease factors, indicative of pleiotropic effects, were interpreted with caution.

Then, proteins identified by discovery analysis were selected for replication analysis in further multi-center MR studies. Validation alternated between the two sets of pQTL and outcome data from FinnGen and UKB (including three validation sets: ARIC/FinnGen, deCODE/UKB, and deCODE/FinnGen), using genome-wide significant SNPs as genetic instruments. The stability of causal associations was evaluated through meta-analysis with a value of I²>50% indicating significant heterogeneity, necessitating a random-effects model [34].

Bayesian colocalization analysis

Colocalization analysis [35] was employed to determine whether a specific genetic variant influenced both an exposure factor and an outcome by modulating gene expression at common loci. Bayesian analysis to calculate the posterior probability of a shared causal variant influencing two traits was performed using the R package ‘coloc’ (version 5.0, available at https://github.com/chr1swallace/coloc) with a default prior probabilities: a prior probability of 1e-4 for any single SNP being associated with each trait (P1 and P2) and a prior probability of 1e-5 for a SNP being associated with both traits (P12) [36]. Assuming a single causal variant, four hypotheses were considered: H0: no causal variants for either trait; H1: a causal variant for the first trait only; H2: a causal variant for the second trait only; H3: separate causal variants for both traits and H4: a shared causal variant for both traits [37]. Significant colocalization was inferred when the posterior probability of H4 was > 0.8, implying strong evidence of a shared causal influence [38].

Extra validation analysis

Considering that gene expression and protein synthesis are influenced by numerous factors beyond simple genetic processes, we conducted extra validation analyses to verify the results of our discovery analysis at both the tissue and protein levels. This validation was performed using both TWAS and PWAS methods predictive of gene influence on phenotype generated by Functional Summary-Based Imputation (FUSION) software (available at http://gusevlab.org/projects/fusion), based on the utility of GWAS summary statistics to indicate associations between GWAS phenotypes and functional phenotypes. TWAS indicated the association of protein-coding genes with allergic disease risk at the tissue level and was used as external validation analysis which utilized the pre-computed eQTL reference panel for target proteins derived from the GTEx8 (Genotype-Tissue Expression version 8) database. Likewise, the PWAS served as an internal validation analysis that integrated the GWAS summary statistics and the pre-computed plasma proteome reference weight also from the ARIC study [17] to calculate the genetic impact on allergic disease. Thus, the impact of significant SNPs from the GWAS on protein abundance could be evaluated and candidate genes linked to allergic disease that regulate plasma protein levels identified. An similar FDR corrected p value < 0.05 was the threshold of significance in the extra validation analysis.

A PPI network was constructed using the Search Tool for the Retrieval of Interacting Genes (STRING) database (version 11.5) [39] with a minimum required interaction score (IAS) threshold of 0.4 to indicate interactions among identified proteins and pre-existing anti-allergy drug targets [40]. Information on anti-allergy drug targets was sourced from the DrugBank database.

Evidence-based grading of potential drug targets

Proteins were graded according to the criteria of Feihong Ren [41].

  • Tier 1 Targets: substantial evidence (PPH4 > 0.8) for drug targeting, confirmed by replication analysis and extra TWAS or PWAS validation.

  • Tier 2 Targets: direct linkage to known drug targets within the PPI network, validated either by replication analysis or extra TWAS or PWAS.

  • Tier 3 Targets: proteins with a PPH4 > 0.8, validated by either replication analysis or extra TWAS criteria or linked to known drug targets within the PPI network.

  • Tier 4 Targets: Proteins not classified under the first three tiers.

Results

Discovery MR analysis

A total of 1,544 proteins associated with 2,810 SNPs (S. Table 2) were identified during initial discovery analysis. Chromosomal locations are shown by Manhattan plots (Fig. 2). After FDR adjustment, 12 proteins were linked to AA, 6 to AD and 10 to AR but no proteins associated with AC or AU were found (Fig. 3). APOE (OR: 1.1591, 95% CI: 1.0972, 1.2245), MAX (OR: 1.0710, CI: 1.0326, 1.1109), NPNT (OR: 1.1607, CI: 1.1079, 1.2161), PILRA (OR: 1.0231, CI: 1.0110, 1.0354), STAT6 (OR: 1.5608, CI: 1.3622, 1.7884) and VTA1 (OR: 1.1213, CI: 1.0661, 1.1794) were associated with increased AA risk and GALK1 (OR: 0.7228, CI: 0.6363, 0.8211), IL1RL2 (OR: 0.8717, CI: 0.8360, 0.9089), IL6R (OR: 0.9632, CI: 0.9459, 0.9807), LRRC32 (OR: 0.7813, CI: 0.6941, 0.8795), PRSS8 (OR: 0.7769; CI: 0.6830, 0.8837) and TNFRSF6B (OR: 0.8038, CI: 0.7285, 0.8868) with decreased risk. ERBB3 (OR: 1.0123, CI: 1.0075, 1.0171), ICAM1 (OR: 1.0018, CI: 1.0008, 1.0028), IL7R (OR: 1.0127, CI: 1.0061, 1.0195), MANF (OR: 1.0071, CI: 1.0033, 1.0110) and STAT6 (OR: 1.0181; CI: 1.0092, 1.0271) were associated with increased risk of AR and FCRLB (OR: 0.9951, CI: 0.9927 to 0.9976), IL1R1 (OR: 0.9908, CI: 0.9864, 0.9952), IL1RL2 (OR: 0.9924, CI: 0.9896, 0.9953), PILRA (OR: 0.9985, CI: 0.9977, 0.9993) and TNFSF8 (OR: 0.9911, CI: 0.9872, 0.9949) with decreased risk. All proteins, with the exception of ERBB3 (OR: 1.2775, CI: 1.1548, 1.4132), were linked to lower AD risk. Results are visualized in forest plots (Fig. 4).

Fig. 2
figure 2

Manhattan plot illustrating the the chromosomal distribution of identified plasma proteins for allergic diseases. The standard line in the plot represents the threshold of FDR P = 0.05. A Allergic asthma; B Allergic conjunctivitis; C Atopic dermatitis; D Allergic rhinitis; E Allergic urticaria

Full size image
Fig. 3
figure 3

Volcano plots of the MR results from the discovery analysis, displaying the associations between 1,544 proteins from ARIC and the risk of allergic diseases. The increased OR for allergic diseases risk is represented as increments in SD of plasma protein levels. Red dots indicate significant proteins. ‘ln’ refers to the natural logarithm; ‘PVE’ stands for the proportion of variance explained. A Allergic asthma; B Allergic conjunctivitis; C Atopic dermatitis; D Allergic rhinitis; E Allergic urticaria

Full size image
Fig. 4
figure 4

Forest plot of MR results from the discovery analysis. AA: Allergic asthma; AD: Atopic dermatitis; AR: Allergic rhinitis

Full size image

Replicative MR and meta-analysis

A total of 1,394 novel proteins with 4,144 SNPs were identified from deCode data by a similar IV selection process (S. Table 3) and used to validate the findings of discovery analysis from the FinnGen datasets (Table 1). Replication analysis using deCode data showed that GALK1, IL1RL2, and TNFRSF6B failed to replicate for AA, as did VTA1 and TNFRSF6B for AD and FCRLB, IL1RL2 and MANF for AR. IL7R could not be replicated during any iteration and was excluded from further analysis. Meta-analysis showed robust associations for other proteins with AA except APOE (p = 0.8196; OR: 1.0152, 95% CI: 0.8918, 1.1557) and LRRC32 (p = 0.1864; OR: 0.8048, CI: 0.5832, 1.1107, S. Figure 1). Significant associations also remained for proteins with AD, apart from LRRC32 (p = 0.0990; OR: 0.6015, 95% CI: 0.3288, 1.1002, S. Figure 2). Only STAT6 and PILRA retained a significant association with AR following post-replication analysis (S. Figure 3). IL7R in AR was excluded from further analysis due to the failure of replication.

Table 1 Detailes of the discovery analysis and replicative validations
Full size table

Sensitivity analysis

Steiger filtering confirmed causal relationship directionality with only the relationship of APOE to AA failing to pass the Steiger test (R2xy = 3.36 × 10-6; p = 0.878). Further bidirectional MR analysis confirmed the significant reverse causal association between AA and APOE (OR: 2.5883, CI: 1.8731, 3.5766, S. Table 4). Consequently, AA was excluded from further analysis.

Phenotype scanning indicated potential pleiotropic effects with APOE, GALK1, ICAM1, MAX, PRSS8, and VTA1 being associated with body mass and blood lipid levels. APOE and PILRA have been previously linked to diabetes, a condition that may be comorbid with allergic diseases. Furthermore, APOE has been associated with the phenotype of maternal diabetes, a condition that correlates with a higher risk of allergic disease in offspring. In addition, ERBB3, IL1R1, IL1RL2, IL6R, LRRC32, STAT6, and TNFRSF6B have been strongly linked to allergic diseases, such as asthma and rhinitis (S. Table 5).

Colocalization analysis and extra validation

Colocalization analysis was performed to assess any shared genetic signals between the proteins and AA or AD with a predefined threshold of PPH4 ≥ 0.8. The majority of proteins were found to colocalize with AA, except for APOE (PPH4 = 1.54%), IL1RL2 (PPH4 = 0.00%), MAX (PPH4 = 75.70%), and PILRA (PPH4 = 0.19%). Only LRRC32 (PPH4 = 27.4%) failed to show colocalization with AD. However, substantial colocalization with AR was only observed for ERBB3 (PPH4 = 98.40%) and TNFSF8 (PPH4 = 88.00%, S. Table 6 and S. Figures 4–6).

Two different sets of reference panels from the GTEx database were used for TWAS of AD skin tissue and AA or AR respiratory tract tissue. IL1RL2 and NPNT were mismatched and only STAT6 (Ptwas = 0.0011), TNFRSF6B (Ptwas = 0.0058), and VTA1 (Ptwas = 0.0052) showed significant association with AA among the remaining ten proteins. Similarly, ICAM1, IL1RL2, TNFSF8, and FCRLB were absent from the database but significant associations with AR were found for IL7R (Ptwas = 0.0023), ERBB3 (Ptwas = 0.0444), MANF (Ptwas = 0.0122) and PILRA (Ptwas = 0.0486). CRAT (Ptwas = 0.0000), TNFRSF6B (Ptwas = 0.0002), and LRRC32 (Ptwas = 0.0000) levels correlated with AD. PWAS was performed on the reference panel from the ARIC study and 4 proteins, STAT6, TNFRSF6B, GALK1, and CRAT, were excluded due to mismatches. Other proteins maintained stable associations with the exception of APOE (ppwas = 0.4915), NPNT (ppwas = 0.1310), and LRRC32 (ppwas = 0.0831) for AA and IL1R1 (ppwas = 0.1330) for AR (Table 2 and S. Figure 7).

Table 2 Assessment of druggability and evidence grading
Full size table

A PPI network was constructed from DrugBank data to illustrate interactions among anti-allergic drug targets and proteins of interest (S.Table 7). Interactions were found between STAT6, TNFRSF6B, IL1RL2, IL6R, and established drug targets, as evidenced by an IAS greater than 0.4, in the PPI network for AA (S. Figure 8). Similarly, interactions were found for ERBB3, IL6R, and MMP12 in the AD-specific PPI network (S. Figure 9) and for ERBB3, ICAM1, IL1RL2, MANF, STAT6 and TNFSF8 in the AR-specific PPI network (S. Figure 10). Further gene-disease enrichment analysis revealed STAT6 was enriched in various allergic diseases, most notably in AA, with a strength of 1.86 and an FDR P-value of 1.8e-13. IL1RL2, IL6R, ERBB3, ICAM1, IL1R1, and GALK1 were also enriched in categories like immune system disease, autoimmune disease, skin disease, and disease of anatomical entity (S.Table 8). We also conducted KEGG and GO enrichment analyses on genes within the PPI network. With the exception of GALK1, NPNT, PRSS8, VTA1, LRRC32, PILRA, MMP12, MANF, TNFSF8, and FCRLB, genes of other proteins were successfully enriched in multiple pathways, particularly STAT6, IL6R, and IL1R1 proteins, which were involved in several immune and inflammation-related pathways. Except for FCRLB, MANF, and TNFSF8, other proteins were all successfully enriched in specific biological processes (see S.Table 9, and 10 and S. Figure 11).

Potential drug targets

Ultimately, 11 proteins for allergic asthma (AA), 6 for atopic dermatitis (AD), and 9 for allergic rhinitis (AR) were evaluated for their potential as drug targets. Among these, STAT6 was identified as excellent Tier 1 potential drug targets for AA as were CRAT and TNFRSF6B for AD. No Tier 1 proteins were identified for AR. Tier 2 proteins included TNFRSF6B, and IL1RL2 for AA; ERBB3, IL6R, MMP12, ICAM1 and IL1RL2 for AD, and ICAM1, IL1RL2, MANF, STAT6 and TNFSF8 for AR. Other proteins were assigned to Tier 3 or below (Table 2).

Discussion

To the best of our knowledge, this study is the first to scrutinize causal associations between plasma proteins and allergic diseases through the integrated approach of MR, colocalization, Steiger filtering analysis, PAV assessment, eQTLs overlap determination, PPI analysis, pathway enrichment, and drug target evaluation. Eleven plasma proteins were identified with links to AA, AD, and AR. MR analysis and extra validation revealed STAT6 and TNFRSF6B to be Tier 1 and IL1RL2 and IL6R to be Tier 2 proteins with the potential for AA treatment. Two Tier 1 proteins, CRAT and TNFRSF6B, and five Tier 2 proteins, ERBB3, IL6R, MMP12, ICAM1, and IL1RL2, were linked to AD and three Tier 2 proteins, MANF, STAT6, and TNFSF8, to AR.

Many novel biomarkers have been identified by proteomic and metabolomic analyses, although studies on allergic diseases have generally used low throughput methods [42]. Niet-Fontarigo et al. [43] identified 18 potential biomarkers of asthma phenotype and disease severity, including HSPG2 and IGFALS for AA, through a bottom-up/non-targeted proteomics approach. The sample size was modest with 32 healthy controls, 43 AR patients, and 192 asthmatics and failed to distinguish protein biomarkers from pathogenic factors for allergic disease due to the reverse causal characteristics of an observational study.

Several proteins identified by the current study have previously been linked to allergic disease by epidemiological or laboratory studies. Indeed, STAT6 is known to participate in IL-4 signaling and its role in asthma has been extensively studied since both doctor-diagnosed asthma and blood eosinophil counts are known to be linked to STAT6 signaling and the IL-1 receptor family [44]. Baris S et al. [45] have identified a novel inborn error of immunity arising from a STAT6 gain-of-function mutation causing severe allergic dysregulation which is treated by Janus kinase inhibitor therapy. The TNFRSF6 (also called Fas) receptor binds TNFSF6 (FasL) ligands expressed on CD8 + T cells and oligodendrocytes [46,47,48]. Th1 cells secrete IFN-γ to activate the Fas/FasL system and induce keratinocyte apoptosis in the spongiosis area which may influence the progression of AD. INF-γ and Fas ligand are secreted by activated CD4 + T cells, TNFRSF6 expressed on keratinocytes, and tumor necrosis factor (TNF)-α secreted by both the activated CD4 + T cells and keratinocytes, with cell-mediated cytotoxicity induced by perforin and granzyme B released by CD8 + T cells [46,47,48].

The IL-36 receptor (IL-36R, IL-1Rrp2, IL1RL2, or IL-1R6) binds all α, β, and γ members of the IL-36 family. IL-36R is expressed in skin, mammary, and mucosal epithelial cell lines and IL-36 mediates intracellular signaling through the IL-36R and IL-1 receptor accessory protein (IL-1RAcP) [49]. IL-36α, IL-36β, and IL-36γ bind IL-36R, form a signal transduction complex with IL-1RAcP, and recruit myeloid differentiation factor 88 to activate mitogen-activated protein kinases mediated by c-Jun N-terminal kinase, extracellular regulated protein kinases 1/2 and the nuclear factor kappa B pathway. The resulting inflammatory mediators have roles in adaptive immunity. IL-36 cytokines released from keratinocytes increase the immunoglobulin (Ig)E production mediated by IL-4 in B cells from AD patients and treatment with anti-IL-36R antibodies decreases IgE and alleviates the disease phenotype [50, 51]. The RNA helicase, DDX5, which regulates the alternative splicing of IL-36R pre-mRNA was found to be down-regulated in keratinocytes from AD patients which promoted the inflammatory response [52]. TNFSF8 (CD30L) is a ligand for the cell surface antigen and marker for Hodgkin lymphoma and related hematologic malignancies, TNFRSF8/CD30. It is considered that inhibition of CD30L/CD30 signaling may constitute a novel biological therapy for AR, since CD30L was shown to amplify Th2 cell effector response in animal models of AR. In vivo treatment with anti-CD30 antibody suppressed AR development and this may be a sufficient target for the treatment of allergic inflammation [53].

Some novel proteins were suggested to have potential causal effects on allergic diseases by the current work. For example, CRAT is a mitochondrial enzyme that transfers acetyl groups between CoA and carnitine during lipid metabolism and links with dermatitis have not been previously reported. CRAT is a key regulator of mitochondrial dysfunction-induced cellular senescence in dermal fibroblasts [54, 55]. Silencing of CRAT is known to cause mitochondrial dysfunction, inflammation and senescence via activation of the cGAS-STING and NF-ĸB pathways [54]. In addition, functional variants of the IL6R have been linked to increased risk of AA and AD but mechanisms remain unclear, although IL-6/soluble IL-6R trans-signaling may affect AD and AA development [56,57,58]. Genetic variants of ERBB3 have been identified as AD susceptibility factors and serum MMP12 may be an indicator of AD and AR disease pathways [57, 59,60,61]. Adhesion molecules are known to be involved in T cell homing to skin lesions in AD patients, one example being ICAM-1 which is highly expressed and may have a pathogenic role [62, 63]. Lastly, little attention has been paid to MANF, although the protein is measurable in serum and reflective of extracellular biomarkers in AD [64]. However, these putative mechanisms are speculative and experimental mechanistic studies are required to extend the findings of the present study.

We acknowledge several limitations to the current study. First, the current focus was on serum proteins which may differ from those within cells and tissues and should also be explored for disease associations. Second, European participants accounted for the vast majority of the current cohort and findings may not be generalizable to populations with different ethnicities. Thirdly, a cis-pQTL coding variant might change a protein’s amino acid sequence without necessarily impacting its function or level. Equating sequence alterations with functional changes could lead to incorrect conclusions. Caution should be exercised in the interpretation of these results. Last, publicly available datasets were used and represent data resources for target identification which are not new, although novel insights and perspectives may be drawn from them.

Conclusion

A MR analysis was conducted to explore the proteomic pathogenesis of allergic disease. Examples of Tier 1 and 2 protein targets that are promising drug target candidates were STAT6, TNFRSF6B, IL1RL2, and IL6R for AA; CRAT, TNFRSF6B, ERBB3, IL6R, MMP12, ICAM1 and IL1RL2 for AD, and ICAM1, IL1RL2, MANF, STAT6 and TNFSF8 for AR. These proteins may influence the development of allergic diseases and expose new diagnostic and therapeutic targets. Further experiments are required to validate the current findings regarding proteinaceous allergic disease markers.

Availability of data and materials

All analyses were conducted using publicly available data. The data that support this study are openly available in UK Biobank at https://www.ukbiobank.ac.uk/, and FinnGen, at https://www.finngen.fi/en. Code Availability: the analysis code in R is available on request.

References

  1. Julia V, Macia L, Dombrowicz D. The impact of diet on asthma and allergic diseases. Nat Rev Immunol. 2015;15:308–22.

    Article  CAS  PubMed  Google Scholar 

  2. Yamamoto T, Endo Y, Onodera A, Hirahara K, Asou HK, Nakajima T, Kanno T, Ouchi Y, Uematsu S, Nishimasu H, et al. DUSP10 constrains innate IL-33-mediated cytokine production in ST2(hi) memory-type pathogenic Th2 cells. Nat Commun. 2018;9:4231.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Aw M, Penn J, Gauvreau GM, Lima H, Sehmi R. Atopic march: collegium internationale allergologicum update 2020. Int Arch Allergy Immunol. 2020;181:1–10.

    Article  PubMed  Google Scholar 

  4. Toledo MF, Saraiva-Romanholo BM, Oliveira RC, Saldiva PH, Silva LF, Nascimento LF, Solé D. Changes over time in the prevalence of asthma, rhinitis and atopic eczema in adolescents from Taubaté, São Paulo, Brazil (2005–2012): relationship with living near a heavily travelled highway. Allergol Immunopathol (Madr). 2016;44:439–44.

    Article  CAS  PubMed  Google Scholar 

  5. GBD 2019 Diseases and Injuries Collaborators. Global burden of 369 diseases and injuries in 204 countries and territories, 1990-2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet. 2020;396:1204–22. https://doi.org/10.1016/S0140-6736(20)30925-9. Erratum in: Lancet. 2020 Nov 14;396(10262):1562.

  6. GBD 2019 IMID Collaborators. Global, regional, and national incidence of six major immune-mediated inflammatory diseases: findings from the global burden of disease study 2019. EClinicalMedicine. 2023;64:102193. https://doi.org/10.1016/j.eclinm.2023.102193.

  7. Fuchs O, Bahmer T, Rabe KF, von Mutius E. Asthma transition from childhood into adulthood. Lancet Respir Med. 2017;5:224–34.

    Article  PubMed  Google Scholar 

  8. Langan SM, Irvine AD, Weidinger S. Atopic dermatitis. Lancet. 2020;396:345–60.

    Article  CAS  PubMed  Google Scholar 

  9. Gudjonsson A, Gudmundsdottir V, Axelsson GT, Gudmundsson EF, Jonsson BG, Launer LJ, Lamb JR, Jennings LL, Aspelund T, Emilsson V, Gudnason V. A genome-wide association study of serum proteins reveals shared loci with common diseases. Nat Commun. 2022;13:480.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Sun BB, Maranville JC, Peters JE, Stacey D, Staley JR, Blackshaw J, Burgess S, Jiang T, Paige E, Surendran P, et al. Genomic atlas of the human plasma proteome. Nature. 2018;558:73–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Imming P, Sinning C, Meyer A. Drugs, their targets and the nature and number of drug targets. Nat Rev Drug Discov. 2006;5:821–34.

    Article  CAS  PubMed  Google Scholar 

  12. Lam MP, Ping P, Murphy E. Proteomics research in cardiovascular medicine and biomarker discovery. J Am Coll Cardiol. 2016;68:2819–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Gashaw I, Ellinghaus P, Sommer A, Asadullah K. What makes a good drug target? Drug Discov Today. 2011;16:1037–43.

    Article  CAS  PubMed  Google Scholar 

  14. Schnabel CL, Jentsch MC, Lübke S, Kaiser-Thom S, Gerber V, Vrtala S, Huang HJ, Rhyner C, Wagner B, Hoffmann R, Volke D. Immunoproteomics reveal increased serum IgG3/5 binding to Dermatophagoides and yeast protein antigens in severe equine asthma in a preliminary study. Front Immunol. 2023;14:1293684.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Asamoah K, Chung KF, Zounemat Kermani N, Bodinier B, Dahlen SE, Djukanovic R, Bhavsar PK, Adcock IM, Vuckovic D, Chadeau-Hyam M. Proteomic signatures of eosinophilic and neutrophilic asthma from serum and sputum. EBioMedicine. 2024;99:104936.

    Article  CAS  PubMed  Google Scholar 

  16. Yao C, Chen G, Song C, Keefe J, Mendelson M, Huan T, Sun BB, Laser A, Maranville JC, Wu H, et al. Genome-wide mapping of plasma protein QTLs identifies putatively causal genes and pathways for cardiovascular disease. Nat Commun. 2018;9:3268.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Gusev A, Ko A, Shi H, Bhatia G, Chung W, Penninx BW, Jansen R, de Geus EJ, Boomsma DI, Wright FA, et al. Integrative approaches for large-scale transcriptome-wide association studies. Nat Genet. 2016;48:245–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Zhang J, Dutta D, Köttgen A, Tin A, Schlosser P, Grams ME, Harvey B, Yu B, Boerwinkle E, Coresh J, Chatterjee N. Plasma proteome analyses in individuals of European and African ancestry identify cis-pQTLs and models for proteome-wide association studies. Nat Genet. 2022;54:593–602.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Ferkingstad E, Sulem P, Atlason BA, Sveinbjornsson G, Magnusson MI, Styrmisdottir EL, Gunnarsdottir K, Helgason A, Oddsson A, Halldorsson BV, et al. Large-scale integration of the plasma proteome with genetics and disease. Nat Genet. 2021;53:1712–21.

    Article  CAS  PubMed  Google Scholar 

  20. Sakaue S, Kanai M, Tanigawa Y, Karjalainen J, Kurki M, Koshiba S, Narita A, Konuma T, Yamamoto K, Akiyama M, et al. A cross-population atlas of genetic associations for 220 human phenotypes. Nat Genet. 2021;53:1415–24.

    Article  CAS  PubMed  Google Scholar 

  21. Sliz E, Huilaja L, Pasanen A, Laisk T, Reimann E, Mägi R, Hannula-Jouppi K, Peltonen S, Salmi T, Koulu L, et al. Uniting biobank resources reveals novel genetic pathways modulating susceptibility for atopic dermatitis. J Allergy Clin Immunol. 2022;149:1105-1112.e1109.

    Article  CAS  PubMed  Google Scholar 

  22. Dönertaş HM, Fabian DK, Valenzuela MF, Partridge L, Thornton JM. Common genetic associations between age-related diseases. Nat Aging. 2021;1:400–12.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, Laurin C, Burgess S, Bowden J, Langdon R, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7:7.

    Article  Google Scholar 

  24. Burgess S, Small DS, Thompson SG. A review of instrumental variable estimators for mendelian randomization. Stat Methods Med Res. 2017;26:2333–55.

    Article  PubMed  Google Scholar 

  25. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. 2013;37:658–65.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57:289–300.

    Article  Google Scholar 

  27. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44:512–25.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Burgess S, Thompson SG. Interpreting findings from mendelian randomization using the MR-Egger method. Eur J Epidemiol. 2017;32:377–89.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from mendelian randomization between complex traits and diseases. Nat Genet. 2018;50:693–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in mendelian randomization with some Invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40:304–14.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Zhao Q, Wang J, Hemani G, Bowden J, Small DS. Statistical inference in two-sample summary-data Mendelian randomization using robust adjusted profile score. Ann Stat. 2020;48:1742–69. https://doi.org/10.1214/19-AOS1866.

    Article  Google Scholar 

  32. Slob EAW, Burgess S. A comparison of robust mendelian randomization methods using summary data. Genet Epidemiol. 2020;44:313–29.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Kamat MA, Blackshaw JA, Young R, Surendran P, Burgess S, Danesh J, Butterworth AS, Staley JR. PhenoScanner V2: an expanded tool for searching human genotype-phenotype associations. Bioinformatics. 2019;35:4851–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ. 2003;327:557–60.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Folkersen L, Fauman E, Sabater-Lleal M, Strawbridge RJ, Frånberg M, Sennblad B, Baldassarre D, Veglia F, Humphries SE, Rauramaa R, et al. Mapping of 79 loci for 83 plasma protein biomarkers in cardiovascular disease. PLoS Genet. 2017;13:e1006706.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Wang G, Sarkar A, Carbonetto P, Stephens M. A simple new approach to variable selection in regression, with application to genetic fine mapping. J R Stat Soc Ser B Stat Methodol. 2020;82:1273–300.

    Article  Google Scholar 

  37. Giambartolomei C, Vukcevic D, Schadt EE, Franke L, Hingorani AD, Wallace C, Plagnol V. Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet. 2014;10:e1004383.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Chen J, Xu F, Ruan X, Sun J, Zhang Y, Zhang H, Zhao J, Zheng J, Larsson SC, Wang X, et al. Therapeutic targets for inflammatory bowel disease: proteome-wide mendelian randomization and colocalization analyses. EBioMedicine. 2023;89:104494.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Szklarczyk D, Kirsch R, Koutrouli M, Nastou K, Mehryary F, Hachilif R, Gable AL, Fang T, Doncheva NT, Pyysalo S, et al. The STRING database in 2023: protein-protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Res. 2023;51:D638-646.

    Article  CAS  PubMed  Google Scholar 

  40. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47:D607-613.

    Article  CAS  PubMed  Google Scholar 

  41. Ren F, Jin Q, Liu T, Ren X, Zhan Y. Proteome-wide mendelian randomization study implicates therapeutic targets in common cancers. J Transl Med. 2023;21:646.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Breiteneder H, Peng YQ, Agache I, Diamant Z, Eiwegger T, Fokkens WJ, Traidl-Hoffmann C, Nadeau K, O’Hehir RE, O’Mahony L, et al. Biomarkers for diagnosis and prediction of therapy responses in allergic diseases and asthma. Allergy. 2020;75:3039–68.

    Article  PubMed  Google Scholar 

  43. Nieto-Fontarigo JJ, González-Barcala FJ, Andrade-Bulos LJ, San-José ME, Cruz MJ, Valdés-Cuadrado L, Crujeiras RM, Arias P, Salgado FJ. iTRAQ-based proteomic analysis reveals potential serum biomarkers of allergic and nonallergic asthma. Allergy. 2020;75:3171–83.

    Article  CAS  PubMed  Google Scholar 

  44. Foster PS, Maltby S, Rosenberg HF, Tay HL, Hogan SP, Collison AM, Yang M, Kaiko GE, Hansbro PM, Kumar RK, Mattes J. Modeling T(H) 2 responses and airway inflammation to understand fundamental mechanisms regulating the pathogenesis of asthma. Immunol Rev. 2017;278:20–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Baris S, Benamar M, Chen Q, Catak MC, Martínez-Blanco M, Wang M, Fong J, Massaad MJ, Sefer AP, Kara A, et al. Severe allergic dysregulation due to a gain of function mutation in the transcription factor STAT6. J Allergy Clin Immunol. 2023;152:182-194.e187.

    Article  CAS  PubMed  Google Scholar 

  46. Kerstan A, Bröcker EB, Trautmann A. Decisive role of tumor necrosis factor-α for spongiosis formation in acute eczematous dermatitis. Arch Dermatol Res. 2011;303:651–8.

    Article  CAS  PubMed  Google Scholar 

  47. Trautmann A, Akdis M, Kleemann D, Altznauer F, Simon HU, Graeve T, Noll M, Bröcker EB, Blaser K, Akdis CA. T cell-mediated Fas-induced keratinocyte apoptosis plays a key pathogenetic role in eczematous dermatitis. J Clin Invest. 2000;106:25–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Yawalkar N, Schmid S, Braathen LR, Pichler WJ. Perforin and granzyme B may contribute to skin inflammation in atopic dermatitis and psoriasis. Br J Dermatol. 2001;144:1133–9.

    Article  CAS  PubMed  Google Scholar 

  49. Peñaloza HF, van der Geest R, Ybe JA, Standiford TJ, Lee JS. Interleukin-36 cytokines in infectious and non-infectious lung diseases. Front Immunol. 2021;12:754702.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Ahmad F, Alam MA, Ansari AW, Jochebeth A, Leo R, Al-Abdulla MN, Al-Khawaga S, AlHammadi A, Al-Malki A, Al Naama K, et al. Emerging role of the IL-36/IL-36R Axis in multiple inflammatory skin diseases. J Invest Dermatol. 2024;144:206–24.

    Article  CAS  PubMed  Google Scholar 

  51. Fukaura R, Akiyama M. Targeting IL-36 in inflammatory skin diseases. BioDrugs. 2023;37:279–93.

    Article  CAS  PubMed  Google Scholar 

  52. Ni X, Xu Y, Wang W, Kong B, Ouyang J, Chen J, Yan M, Wu Y, Chen Q, Wang X, et al. IL-17D-induced inhibition of DDX5 expression in keratinocytes amplifies IL-36R-mediated skin inflammation. Nat Immunol. 2022;23:1577–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Fuchiwaki T, Sun X, Fujimura K, Yamada H, Shibata K, Muta H, Podack ER, Kawauchi H, Yoshikai Y. The central role of CD30L/CD30 interactions in allergic rhinitis pathogenesis in mice. Eur J Immunol. 2011;41:2947–54.

    Article  CAS  PubMed  Google Scholar 

  54. Song MJ, Park CH, Kim H, Han S, Lee SH, Lee DH, Chung JH. Carnitine acetyltransferase deficiency mediates mitochondrial dysfunction-induced cellular senescence in dermal fibroblasts. Aging Cell. 2023;22:e14000.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Davies MN, Kjalarsdottir L, Thompson JW, Dubois LG, Stevens RD, Ilkayeva OR, Brosnan MJ, Rolph TP, Grimsrud PA, Muoio DM. The acetyl group buffering action of carnitine acetyltransferase offsets macronutrient-induced lysine acetylation of mitochondrial proteins. Cell Rep. 2016;14:243–54.

    Article  CAS  PubMed  Google Scholar 

  56. Zhang M, Bai Y, Wang Y, Cui H, Tang M, Wang L, Wang X, Gu D. Cumulative evidence for associations between genetic variants in interleukin 6 receptor gene and human diseases and phenotypes. Front Immunol. 2022;13:860703.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Esparza-Gordillo J, Schaarschmidt H, Liang L, Cookson W, Bauerfeind A, Lee-Kirsch MA, Nemat K, Henderson J, Paternoster L, Harper JI, et al. A functional IL-6 receptor (IL6R) variant is a risk factor for persistent atopic dermatitis. J Allergy Clin Immunol. 2013;132:371–7.

    Article  CAS  PubMed  Google Scholar 

  58. Wang Y, Hu H, Wu J, Zhao X, Zhen Y, Wang S, Li W, Liang M, Wu B, Ma G. The IL6R gene polymorphisms are associated with sIL-6R, IgE and lung function in Chinese patients with asthma. Gene. 2016;585:51–7.

    Article  CAS  PubMed  Google Scholar 

  59. Martel BC, Litman T, Hald A, Norsgaard H, Lovato P, Dyring-Andersen B, Skov L, Thestrup-Pedersen K, Skov S, Skak K, Poulsen LK. Distinct molecular signatures of mild extrinsic and intrinsic atopic dermatitis. Exp Dermatol. 2016;25:453–9.

    Article  CAS  PubMed  Google Scholar 

  60. Brunner PM, Suárez-Fariñas M, He H, Malik K, Wen HC, Gonzalez J, Chan TC, Estrada Y, Zheng X, Khattri S, et al. The atopic dermatitis blood signature is characterized by increases in inflammatory and cardiovascular risk proteins. Sci Rep. 2017;7:8707.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Zhou Y, Xu M, Gong W, Kang X, Guo R, Wen J, Zhou D, Wang M, Shi D, Jing Q. Circulating MMP-12 as potential biomarker in evaluating disease severity and efficacy of sublingual immunotherapy in allergic rhinitis. Mediators Inflamm. 2022;2022:3378035.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Marinović Kulišić S, Takahashi M, Himelreich Perić M, Mužić Radović V, Jurakić Tončić R. Immunohistochemical analysis of adhesion molecules e-selectin, intercellular adhesion molecule-1, and vascular cell adhesion molecule-1 in inflammatory lesions of atopic dermatitis. Life (Basel). 2023;13:933.

    PubMed  Google Scholar 

  63. Lugović L, Cupić H, Lipozencić J, Jakić-Razumović J. The role of adhesion molecules in atopic dermatitis. Acta Dermatovenerol Croat. 2006;14:2–7.

    PubMed  Google Scholar 

  64. Bang H, Kim JE, Lee HS, Park SM, Park DJ, Lee EJ. Integrated bioinformatic analysis of gene expression profiling data to identify combinatorial biomarkers in inflammatory skin disease. Sci Rep. 2022;12:5889.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors would like to express their gratitude to EditSprings (https://www.editsprings.com/) for the expert linguistic services provided.

Funding

This work was supported by the National Key R&D Program of China [No. 2019YFA0112100]. The study funders/sponsors had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Author information

Authors and Affiliations

  1. Department of Spine Surgery and Orthopaedics, Xiangya Hospital, Central South University, Xiangya Road 87, Changsha, 410000, China

    Ziqin Cao & Jianhuang Wu

  2. National Clinical Research Center for Geriatric Disorders, Central South University, Changsha, China

    Ziqin Cao, Qiangxiang Li, Yajia Li & Jianhuang Wu

  3. Department of Dermatology, Xiangya Hospital, Central South University, Changsha, Hunan, 410011, China

    Yajia Li

Authors
  1. Ziqin Cao
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Qiangxiang Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Yajia Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Jianhuang Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Dr. Yajia Li and Pro. Jianhuang Wu had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Concept and design: Jianhuang Wu and  Yajia Li; Acquisition, analysis, or interpretation of data:Ziqin Cao and  Yajia Li; Drafting of the manuscript: All of the authors; Critical revision of the manuscript for important intellectual content: Ziqin Cao and Jianhuang Wu; Statistical analysis: Ziqin Cao, Yajia Li and Qiangxiang Li; Administrative, technical, or material support: Ziqin Cao and Qiangxiang Li. Dr. Yajia Li and Pro. Jianhuang Wu contributed equally to the work.

Corresponding authors

Correspondence to Yajia Li or Jianhuang Wu.

Ethics declarations

Ethics approval and consent to participate

The manuscript does not contain clinical studies or patient data. Our study is based on the large-scale GWAS datasets, and not the individual-level data. All participants gave informed consent in all the corresponding original studies, and no additional ethical approval was applicable. No animal subjects were used in this work.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cao, Z., Li, Q., Li, Y. et al. Identification of plasma protein markers of allergic disease risk: a mendelian randomization approach to proteomic analysis. BMC Genomics 25, 503 (2024). https://doi.org/10.1186/s12864-024-10412-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-024-10412-0

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us