Translational Vision Science & Technology Cover Image for Volume 14, Issue 4
April 2025
Volume 14, Issue 4
Open Access
Cornea & External Disease  |   April 2025
Metabolic Biomarkers Mediate Allergic Conjunctivitis via Circulating Inflammatory Proteins: Evidence From a Mendelian Randomization Study
Author Affiliations & Notes
  • Xiang Cao
    Department of Ophthalmology, Affiliated People's Hospital, Jiangsu University, Zhenjiang, Jiangsu, China
    Zhenjiang Kangfu Eye Hospital, Zhenjiang, Jiangsu, China
  • Zijiao Xu
    School of Medicine, Nankai University, Tianjin, China
    Tianjin Eye Hospital, Nankai University Affiliated Eye Hospital, Clinical College of Ophthalmology, Tianjin Medical University, Tianjin Eye Institute, Tianjin Key Laboratory of Ophthalmology and Visual Science, Tianjin, China
  • Boyang Zhang
    School of Medicine, Nankai University, Tianjin, China
    Tianjin Eye Hospital, Nankai University Affiliated Eye Hospital, Clinical College of Ophthalmology, Tianjin Medical University, Tianjin Eye Institute, Tianjin Key Laboratory of Ophthalmology and Visual Science, Tianjin, China
  • Qingyu Li
    Tianjin Eye Hospital, Nankai University Affiliated Eye Hospital, Clinical College of Ophthalmology, Tianjin Medical University, Tianjin Eye Institute, Tianjin Key Laboratory of Ophthalmology and Visual Science, Tianjin, China
  • Zhixin Jiang
    Tianjin Eye Hospital, Nankai University Affiliated Eye Hospital, Clinical College of Ophthalmology, Tianjin Medical University, Tianjin Eye Institute, Tianjin Key Laboratory of Ophthalmology and Visual Science, Tianjin, China
  • Xiaoyong Yuan
    Tianjin Eye Hospital, Nankai University Affiliated Eye Hospital, Clinical College of Ophthalmology, Tianjin Medical University, Tianjin Eye Institute, Tianjin Key Laboratory of Ophthalmology and Visual Science, Tianjin, China
  • Correspondence: Zhixin Jiang, Tianjin Eye Hospital, Nankai University Affiliated Eye Hospital, Clinical College of Ophthalmology, Tianjin Medical University, Tianjin Eye Institute, Tianjin Key Laboratory of Ophthalmology and Visual Science, No. 4 Gansu Road, Heping District, Tianjin 300020, China. e-mail: [email protected] 
  • Xiaoyong Yuan, Tianjin Eye Hospital, Nankai University Affiliated Eye Hospital, Clinical College of Ophthalmology, Tianjin Medical University, Tianjin Eye Institute, Tianjin Key Laboratory of Ophthalmology and Visual Science, No. 4 Gansu Road, Heping District, Tianjin 300020, China. e-mail: [email protected] 
  • Footnotes
     XC and ZX contributed equally to this article and should be considered co-first authors.
Translational Vision Science & Technology April 2025, Vol.14, 12. doi:https://doi.org/10.1167/tvst.14.4.12
  • Views
  • PDF
  • Share
  • Tools
    • Alerts
      ×
      This feature is available to authenticated users only.
      Sign In or Create an Account ×
    • Get Citation

      Xiang Cao, Zijiao Xu, Boyang Zhang, Qingyu Li, Zhixin Jiang, Xiaoyong Yuan; Metabolic Biomarkers Mediate Allergic Conjunctivitis via Circulating Inflammatory Proteins: Evidence From a Mendelian Randomization Study. Trans. Vis. Sci. Tech. 2025;14(4):12. https://doi.org/10.1167/tvst.14.4.12.

      Download citation file:


      © ARVO (1962-2015); The Authors (2016-present)

      ×
  • Supplements
Abstract

Purpose: This study aimed to investigate the mediating role of metabolic biomarkers (MBs) in the association between circulating inflammatory proteins (CIPs) and allergic conjunctivitis (AC) and identify potential therapeutic targets.

Methods: A Mendelian randomization (MR) study was conducted, leveraging genetic variants as instrumental variables to infer causal relationships. Data were obtained from genome-wide association studies (GWASs), and a two-sample MR was employed to estimate the direct and indirect effects of CIPs on AC through MBs. Inverse-variance weighting (IVW) served as the primary analysis method, supplemented by sensitivity analyses to assess the robustness of findings.

Results: Six CIPs were found to have significant causal effects on AC. Natural killer cell receptor 2B4 (CD244) exhibited a protective effect, and interleukin-18 receptor 1 (IL-18R1), IL-6, leukemia inhibitory factor (LIF), protein S100-A12 (EN-RAGE), and TNF-related activation-induced cytokine (TRANCE) were identified as risk factors. The MR analysis indicated the mediation role of specific MBs in these associations, with MBs such as 4-oxo-retinoic acid, gulonate, 3-(4-hydroxyphenyl) lactate, S-adenosylhomocysteine, and sphingomyelin, significantly influencing the pathway from CIPs to AC.

Conclusions: This study highlights the crucial role of MBs in mediating the association between CIPs and AC. These findings offer new insights into the pathophysiology of AC and suggest potential metabolic targets for novel therapeutic approaches.

Translational Relevance: This study underscores the potential for targeting specific MBs as novel therapeutic approaches to modulate the inflammatory pathways implicated in AC.

Introduction
Allergic conjunctivitis (AC) is a common immune-mediated hypersensitivity disorder characterized by ocular symptoms such as itching, redness, and tearing.13 These symptoms result from a complex immune response, where environmental allergens such as pollen or dust mites trigger an inflammatory cascade in genetically predisposed individuals.4,5 The pathophysiology of AC involves interactions between genetic and environmental factors, leading to inflammation in the conjunctiva and subsequent tissue damage.68 In addition to these ocular manifestations, AC often results in fatigue and difficulty concentrating, which can severely impact an individual's quality of life and daily activities.9 AC remains a significant public health challenge, severely diminishing quality of life and imposing a substantial economic burden on society due to direct medical costs and lost productivity.3,6,10 Although current treatments, such as antihistamines and corticosteroids, are effective in alleviating symptoms, they do not address the underlying mechanisms of the disease.1,2 This gap in therapeutic options underscores the need for a deeper understanding of the molecular processes involved, which could lead to the development of more targeted and effective therapies. 
Recent research has highlighted the crucial roles of circulating inflammatory proteins (CIPs) and metabolic biomarkers (MBs) in allergic diseases.4,8,11 CIPs play a key role in regulating immune and inflammatory responses, with their levels often altered in various diseases, which is highly valuable in clinical medicine.12 MBs serve as the footprints of biological processes and can function as signaling mediators and immunomodulators, significantly impacting the pathophysiology of diseases.7,13 Moreover, serum metabolomics provides insights into the complex effects of diseases on individuals, helping to identify metabolic changes that occur in response to inflammatory processes.14 Studies have demonstrated that IL-33 is a key inflammatory mediator, playing a central role in modulating immune responses by inducing T helper 2 (Th2) cell immune responses in refractory allergic diseases.4,8 Similarly, omega-3 fatty acids and their MBs (e.g., resolvins, protectins) have shown significant anti-inflammatory properties in combating allergic diseases.15 These MBs, rather than the omega-3 fatty acids themselves, mediate the anti-inflammatory effects, and their levels can be measured in the blood. However, the precise mechanisms by which these biomarkers influence the onset and progression of AC remain unclear. Understanding these causal relationships in AC could lead to new therapeutic targets that are more effective and specific than current treatments. 
Traditional observational studies often struggle to establish causal links between CIPs and allergic diseases due to biases and reverse causality.16 Biases arise from confounding factors such as environmental exposures, lifestyle choices, and comorbidities that can influence both CIP levels and disease risk.16,17 Reverse causality occurs when the disease itself alters CIP levels, making it challenging to discern cause from effect. In AC research, these challenges are intensified by diverse environmental and genetic factors.7,18 Environmental influences include exposure to allergens such as pollen, dust mites, and pet dander, as well as air pollution and climate variations.1,6 Genetic predispositions, particularly variations in genes regulating immune responses, further complicate the landscape.5 The interplay of these factors introduces significant variability, hindering the isolation of specific causal relationships. Mendelian randomization (MR) offers a solution by using genetic variants as instrumental variables to infer causality.19 Because genetic variants are randomly assigned at conception, they are generally free from the confounding factors and reverse causation that plague observational studies.17,19 This method enhances the robustness of causal inferences, providing clearer insights into disease etiology.17 For example, a study by Zhou et al.20 utilized MR to establish a causal link between atopic dermatitis and AC, demonstrating the utility of this approach in identifying causal relationships. Understanding these causal pathways is crucial for developing effective prevention and treatment strategies for AC. By minimizing biases inherent in traditional studies, MR enables a more accurate assessment of the role of CIPs in allergic diseases, paving the way for innovative therapeutic interventions.18 In this study, we pioneered the use of a mediated MR approach to explore the roles of 1400 MBs in the association of 91 CIPs with AC. By integrating genome-wide association study (GWAS) data and relevant biomarkers, we aimed to illuminate the intricate connections between genetic susceptibility, inflammatory responses, and allergic diseases. Our findings promise to deepen our understanding of this prevalent allergic condition, ultimately guiding the development of targeted and effective interventions. 
Methods
Study Design
In this study, a two-step MR was employed to investigate the mediating role of MBs in the relationship between CIPs and the risk of AC. Single-nucleotide polymorphisms (SNPs) were defined as instrumental variables (IVs). MR studies are based on three core assumptions: (1) IVs are strongly associated with exposure; (2) IVs are unrelated to confounding factors; and (3) IVs can only affect the outcome through exposure.17 The design of this study is illustrated in Figure 1
Figure 1.
 
Schematic diagram of this study. (A) The total effect c was analyzed using genetically predicted CIPs as the exposure and AC as the outcome. Additionally, the total effect d was also analyzed using genetically predicted AC as the exposure and CIPs as the outcome. The relevance and independence of IVs were considered for both exposure and outcome, and potential confounders were accounted for. (B) The total effect was further broken down as follows: The indirect effect was assessed using a two-step approach involving the effects of CIPs on MBs (a) and the effects of MBs on AC (b), as well as the product method (a × b), and the direct effect was calculated as c′ = ca × b. The proportion mediated was determined by dividing the indirect effect by the total effect.
Figure 1.
 
Schematic diagram of this study. (A) The total effect c was analyzed using genetically predicted CIPs as the exposure and AC as the outcome. Additionally, the total effect d was also analyzed using genetically predicted AC as the exposure and CIPs as the outcome. The relevance and independence of IVs were considered for both exposure and outcome, and potential confounders were accounted for. (B) The total effect was further broken down as follows: The indirect effect was assessed using a two-step approach involving the effects of CIPs on MBs (a) and the effects of MBs on AC (b), as well as the product method (a × b), and the direct effect was calculated as c′ = ca × b. The proportion mediated was determined by dividing the indirect effect by the total effect.
Data Source
The GWAS summary statistics used in this study were publicly available and involved secondary analysis of previously published data, thereby exempting the need for ethical approval. Genetic associations for AC were obtained from the FinnGen consortium, which includes 23,665 cases and 388,516 controls, primarily of European descent. Further details about this dataset can be accessed at https://www.finngen.fi/en/access_results. The FinnGen project was selected due to its large sample size, data diversity, and rigorous quality control standards, making it an ideal resource for GWAS summary data on AC. Ninety-one CIPs were derived from the GWAS summary data of 14,824 individuals of European ancestry.21 This dataset was chosen for its comprehensive coverage of inflammatory proteins and their genetic associations. Complete GWAS summary statistics for these proteins are available at https://www.phpc.cam.ac.uk/ceu/proteins
The plasma metabolome data were sourced from the Canadian Longitudinal Study of Aging cohort, which included 8299 participants.22 This dataset offers a broad analysis of plasma metabolite levels. The summary data for 1400 MBs included 1091 blood metabolites and 309 metabolite ratios, obtained from the GWAS Catalog (https://www.ebi.ac.uk/gwas). The dataset was selected for its comprehensive coverage of metabolic biomarkers relevant to allergic diseases and inflammation.7 
Selection of Genetic Instruments
To meet the correlation criteria and obtain the required number of SNPs for analysis, the included SNPs had to be genome-wide significant (P < 1 × 10–5).11 SNPs with r2 < 0.001 were selected to ensure independence, with a distance greater than 10,000 kb between each pair of SNPs to avoid linkage disequilibrium due to physical proximity.23 In the reverse MR analysis, SNPs in the AC summary data were selected with a threshold of P < 5 × 10–8, and linkage disequilibrium was eliminated using the same method (kb = 10,000, r2 < 0.001).17 Palindromic SNPs with the same allele on the forward and reverse strands were excluded due to ambiguity in determining the strand of the effect allele.24 Weak IVs were excluded by calculating the F statistic, F = R2(nk – 1)/k(1-R2), with weak IVs (F < 10) being removed.25 
Statistical and Sensitivity Analysis
A schematic diagram of the study analysis is presented in Figure 1. Initially, bidirectional MR analysis was conducted to assess the potential causal relationship between SNPs associated with 91 CIPs and AC, determining the total effect (c). Subsequently, the two-step MR analysis was performed to investigate the mediating role of MBs. The first step involved examining the effects of CIPs on MBs, followed by assessing the impact of MBs on AC. The indirect effect (a × b) was calculated as the product of these estimates, and the direct effect (c′) was derived by subtracting the indirect effect from the total effect (ca × b).26 
Five MR analysis methods were employed, including inverse variance weighting (IVW), Mendelian randomization–Egger (MR-Egger), weighted median method, simple model, and weighted model method.17 Results primarily relied on IVW, with statistical significance defined as P < 0.05.26 Heterogeneity was evaluated using the Cochran Q test, with P < 0.05 as the threshold for heterogeneity.27 The IVW random-effects model was utilized in this scenario, and the fixed-effects model was employed otherwise.28 To enhance robustness, the MR-Egger method and Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO) method were used to identify horizontal pleiotropy.29,30 MR-Egger regression assessed horizontal pleiotropy of SNPs acting as IVs through the deviation from zero of the intercept term, with a non-zero intercept indicating bias in the IVW estimate.31 The MR-PRESSO outlier test was employed to eliminate anomalous SNPs and calculate corrected results to further mitigate horizontal pleiotropy.30 Furthermore, leave-one-out analysis was conducted to exclude outlier SNPs and prevent a single SNP from unduly influencing the overall causal estimate. 
These findings were visually depicted through forest plots, funnel plots, scatterplots, and leave-one-out plots. All statistical analyses were carried out using the TwoSampleMR 0.5.6 and MRPRESSO 1.0 packages within R 4.3.2 (R Foundation for Statistical Computing, Vienna, Austria). 
Results
Detailed Information on the Included SNPs of the Analysis
According to the selection criteria for IVs (Fig. 2), a total of 2945 SNPs were selected as IVs for 91 CIPs and 34,843 SNPs were selected as IVs for 1400 MBs (Supplementary Tables S1, S2). 
Figure 2.
 
Flowchart of data source and SNP selection process for exposure (CIPs), mediator (MBs), and outcome (AC). Linkage disequilibrium, LD.
Figure 2.
 
Flowchart of data source and SNP selection process for exposure (CIPs), mediator (MBs), and outcome (AC). Linkage disequilibrium, LD.
Causal Effects Between CIPs and AC
To obtain the total effect (c), the association between CIPs and AC was investigated using a two-sample MR analysis. As shown in Figures 3A and 3B, the IVW method was employed as the primary MR analysis method, which identified six CIPs with a causal relationship with AC. Among them, natural killer cell receptor 2B4 (CD244) (odds ratio [OR] = 0.928; 95% confidence interval [CI], 0.866–0.994; P = 0.0331) had a protective effect, and tumor necrosis factor–related activation-induced cytokine (TRANCE) (OR = 1.045; 95% CI, 1.002–1.089; P = 0.0386), protein S100-A12 (EN-RACE) (OR = 1.076; 95% CI, 1.015–1.089; P = 0.0140), interleukin-18 receptor 1 (IL-18R1) (OR = 1.046; 95% CI, 1.013–1.079; P = 0.0052), leukemia inhibitory factor (LIF) (OR = 1.078; 95% CI, 1.013–1.148; P = 0.0184), and interleukin 6 (IL-6) (OR = 1.129; 95% CI, 1.001–1.273; P = 0.0473) were risk factors. The causal effects were visualized using circular heatmaps (Figs. 3A, 3B). Comprehensive insights into the causal relationship were provided through the scatterplots shown in Figures 3C to 3H and Supplementary Figure S1, with sensitivity analyses showing no substantial heterogeneity and pleiotropy (Supplementary Table S3). Further MR results are detailed in Supplementary Table S4. The reverse MR analysis showed no significant results (Figs. 3I, 3J). 
Figure 3.
 
Bidirectional MR analysis between CIPs and AC. (A) The circular heatmap of P values illustrates the causal relationship between CIPs as exposures and AC as the outcome using various MR methods. The vertical axis represents six CIPs (TRANCE, LIF, IL-6, IL-18R1, EN-RAGE, and CD244), and the horizontal axis, from outermost to innermost, represents five MR methods (weighted mode, simple mode, IVW, weighted median, and MR Egger). Lower P values provide stronger evidence for a causal relationship; the IVW method identified six CIPs with significant causal associations with AC. (B) The circular heatmap of OR visualizes the strength and direction of causal associations between CIPs as exposures and AC as the outcome. Six CIPs showed significant associations, with OR analysis identifying both protective and risk factors. TRANCE, LIF, IL-6, IL-18R1, and EN-RAGE were risk factors, and CD244 was a protective factor. (CH) Scatterplots of causal effects of SNPs associated with TRANCE (C), EN-RAGE (D), IL-18R1 (E), CD244 (F), LIF (G), and IL-6 (H) on AC. These plots demonstrate individual causal relationships and provide visual evidence of consistency across different SNPs. (I) The circular heatmap of P values displays the results of reverse MR analysis with AC as the exposure and CIPs as the outcome. No significant results were observed, indicating the absence of a reverse causal relationship from AC to the six CIPs. (J) Schematic representation of the bidirectional MR analysis results that summarizes the significant causal relationships identified and highlights the lack of significant reverse associations.
Figure 3.
 
Bidirectional MR analysis between CIPs and AC. (A) The circular heatmap of P values illustrates the causal relationship between CIPs as exposures and AC as the outcome using various MR methods. The vertical axis represents six CIPs (TRANCE, LIF, IL-6, IL-18R1, EN-RAGE, and CD244), and the horizontal axis, from outermost to innermost, represents five MR methods (weighted mode, simple mode, IVW, weighted median, and MR Egger). Lower P values provide stronger evidence for a causal relationship; the IVW method identified six CIPs with significant causal associations with AC. (B) The circular heatmap of OR visualizes the strength and direction of causal associations between CIPs as exposures and AC as the outcome. Six CIPs showed significant associations, with OR analysis identifying both protective and risk factors. TRANCE, LIF, IL-6, IL-18R1, and EN-RAGE were risk factors, and CD244 was a protective factor. (CH) Scatterplots of causal effects of SNPs associated with TRANCE (C), EN-RAGE (D), IL-18R1 (E), CD244 (F), LIF (G), and IL-6 (H) on AC. These plots demonstrate individual causal relationships and provide visual evidence of consistency across different SNPs. (I) The circular heatmap of P values displays the results of reverse MR analysis with AC as the exposure and CIPs as the outcome. No significant results were observed, indicating the absence of a reverse causal relationship from AC to the six CIPs. (J) Schematic representation of the bidirectional MR analysis results that summarizes the significant causal relationships identified and highlights the lack of significant reverse associations.
Causal Effects of CIPs on MBs
To obtain the results of the β1 effect (a), the causal effect of CIPs on MBs was analyzed. The two-sample MR analysis demonstrated significant causal effects of CIPs on several MBs. Six CIPs associated with AC were included in the mediation MR analysis. IVW analysis found IL-18R1 protective for 4-oxo-retinoic acid (OR = 0.933; 95% CI, 0.887–0.981; P = 0.0068) but a risk factor for gulonate (OR = 1.068; 95% CI, 1.004–1.136; P = 0.0376). TRANCE increased levels of 3-(4-hydroxyphenyl) lactate (OR = 1.068; 95% CI, 1.004–1.136; P = 0.0376), and protein S100-A12 (EN-RAGE) increased S-adenosylhomocysteine (SAH) levels (OR = 1.130; 95% CI, 1.021–1.251; P = 0.0178). CD244 decreases led to lower sphingomyelin (SM) levels (OR = 0.912; 95% CI, 0.844–0.986; P = 0.0213). These findings are visualized in Supplementary Figure S2 and were confirmed by sensitivity analyses (Supplementary Table S5). Detailed MR results are provided in Supplementary Table S6
Causal Effects of MBs on AC
To obtain the results of the β2 effect (b), the causal effect of MBs on the risk of AC was analyzed. The analysis revealed significant causal effects of certain MBs on the risk of AC. Specifically, increased levels of 3-(4-hydroxyphenyl) lactate (OR = 1.058; 95% CI, 1.003–1.115; P = 0.0375), SM (OR = 1.069; 95% CI, 1.007–1.134; P = 0.0278), SAH (OR = 1.053; 95% CI, 1.007–1.100; P = 0.0231), and gulonate (OR = 1.066; 95% CI, 1.004–1.133; P = 0.0374) were linked to higher AC risk, but 4-oxo-retinoic acid was protective (OR = 0.920; 95% CI, 0.849–0.998; P = 0.0444) (Fig. 4). These relationships were robust in sensitivity analyses (Supplementary Table S7); visualizations are shown in Supplementary Figure S3, and further details are provided in Supplementary Table S8
Figure 4.
 
Forest plot to visualize the causal effects of MBs with CIPs and AC.
Figure 4.
 
Forest plot to visualize the causal effects of MBs with CIPs and AC.
Genetically Predicted MBs Mediate the Association of CIPs With AC
After identifying the effect of exposure on the mediator and the significant mediators affecting AC, the mediated effect proportion and direct effect were quantified. The effects of IL-18R1 on AC were partially mediated by 4-oxo-retinoic acid (mediated proportion = 12.9%; mediated effect = 0.0058) and gulonate (8.36%; –0.0113). The effects of TRANCE were mediated by 3-(4-hydroxyphenyl) lactate (8.36%; 0.0037), EN-RAGE by SAH (8.52%; 0.0063), and CD244 by SM (8.1%; –0.0061) (Supplementary Table S9). Negative mediating roles were noted for gulonate and SM, whereas others had positive mediating roles (Fig. 4). These results indicate that changes in MBs played a significant role in the pathway from CIPs to AC. 
Discussion
The investigation into the roles of CIPs and MBs in AC is crucial for understanding the underlying pathophysiology of this condition.7,32 Allergic diseases are influenced by complex interactions between genetic and environmental factors, with inflammatory proteins and metabolic processes playing significant roles.3,33 In this study, the roles of 91 CIPs and 1400 MBs in the development of AC were initially investigated using the mediation MR approach. Six CIPs were identified to have significant causal effects on AC. CD244 exhibited a protective effect, whereas IL-18R1, LIF, IL-6, EN-RAGE, and TRANCE were identified as risk factors. These findings underscore the critical roles of specific CIPs in the pathogenesis of AC. 
Previous studies have demonstrated that inflammatory cytokines are extensively involved in immune and inflammatory responses.3436 Multiple studies have shown that the application of biologics targeting IL-5, IL-4, and IL-13 antibodies significantly improves the treatment of allergic diseases including asthma, rhinitis, and atopic dermatitis; however, these drugs have not been approved for use in ocular allergic diseases.37 IL-18 and IL-6 function as pro-inflammatory cytokines. IL-18 induces various leukocytes and participates in allergic reactions in conditions such as rhinitis, dermatitis, and asthma.35 IL-18R1, the receptor for IL-18, is essential for mediating IL-18 signaling.38 Inhibiting the binding of IL-18 to IL-18R1 may decrease the risk of inflammatory diseases.38 Research has demonstrated that elevated levels of IL-18R1 are associated with an increased risk of allergies, hay fever, and eczema.39 IL-6 exacerbates allergic reactions by promoting the production of Th2-type cytokines.34 Tao et al.40 found that blocking Th2 signaling can alleviate the inflammatory response and clinical symptoms of allergic conjunctivitis. Moreover, LIF is a pleiotropic cytokine within the IL-6 superfamily.41 Our results support the notion that elevated levels of IL-18R1, IL-6, and LIF are associated with an increased risk of AC. The identification of these inflammatory cytokines as risk factors for AC highlights their potential as therapeutic targets. Targeting IL-18R1 and IL-6 could modulate the inflammatory cascade, potentially leading to new treatments for AC. 
CD244 is a signaling molecule expressed on various immune cells, including natural killer (NK) cells, T cells, monocytes, eosinophils, and mast cells, where it plays a crucial role in immune regulation.42 In CD244–/– mice with mild atopic dermatitis (AD), moderate eosinophil infiltration in the skin suggested that CD244 enhances eosinophil trafficking.43 Furthermore, in a chronic AD model, CD244–/– mice exhibited excessive mast cell degranulation, highlighting the inhibitory role of CD244 in mast cell activation.43 CD244 also participates in eosinophil adhesion and chemotaxis, as its expression increased after the onset of allergic rhinitis (AR).44 This upregulation in AR, a condition closely related to AC, suggests that CD244 may facilitate the recruitment and activation of eosinophils at sites of allergic inflammation. Given the pivotal role of eosinophils in the inflammatory processes of AC, regulation of their trafficking and activation by CD244 could significantly influence the severity of AC. Our findings indicate that CD244 may act as a protective factor against AC. Although its relationship with AC has not yet been studied, it is hypothesized that this protective effect may result from the ability of CD244 to modulate eosinophil function and mast cell activation, which are key drivers of the inflammatory cascade in AC.4,6 
EN-RAGE, a member of the S100 protein family, binds to and activates RAGE signal transduction and is released by various cells under inflammatory conditions, participating in multiple inflammatory responses.45,46 S100A8/9 is highly upregulated in inflammatory bowel disease and can serve as a biomarker of intestinal inflammation due to its high stability in fecal samples.47 Similarly, EN-RAGE is significantly elevated in patients with chronic rhinosinusitis.48 Genetic prediction revealed that elevated EN-RAGE is associated with an increased risk of AC. Although MR analysis suggested that TRANCE was risk factors for allergic conjunctivitis, there is insufficient research supporting their roles in other allergic diseases. The protective effect of CD244 may be related to its role in regulating immune responses, whereas the risk effects of IL-18, IL-6, LIF, EN-RAGE, and TRANCE may be related to their roles in promoting inflammatory responses. In addition, AC itself may influence changes in CIPs. Therefore, we investigated the reverse causality between AC and CIPs and found no significant results. 
MBs, as small-molecule products of metabolic reactions, are influenced by various factors such as genetics and diseases, and they also impact disease risk.22,49 Alterations in the levels of CIPs in the circulatory system may impact metabolic processes, resulting in changes in the levels of MBs.50 Our results highlighted the mediating role of specific metabolites in these associations, with metabolites such as 4-oxo-retinoic acid, gulonate, 3-(4-hydroxyphenyl) lactate, SAH, and SM significantly influencing the pathway from CIPs to AC. Specifically, the effect of IL-18R1 on AC was partially mediated by 4-oxo-retinoic acid and gulonate, TRANCE by 3-(4-hydroxyphenyl) lactate, EN-RAGE by SAH, and CD244 by SM. These findings suggest that metabolites play a crucial role in mediating the effects of CIPs on AC. Therefore, interventions aimed at modulating metabolite levels may represent a novel therapeutic approach. 
Previous studies have demonstrated that these MBs are associated with various diseases and play critical roles in regulating immune and inflammatory responses. For example, 4-oxo-retinoic acid, a bioactive metabolite of vitamin A, primarily exerts its effects through the retinoic acid receptor pathway.51,52 It modulates immune responses by regulating the balance of T helper (Th) cells, specifically promoting regulatory T (Treg) cell differentiation and inhibiting Th2/Th17 responses.52 In the context of allergic inflammation, 4-oxo-retinoic acid may help maintain immune tolerance and reduce inflammation severity by regulating key immune pathways.51,52 Similarly, gulonate, an intermediate product of glucose metabolism, participates in inflammatory processes.53 Recent studies have shown that gulonate increases the risk of age-related macular degeneration, a condition marked by immune dysregulation and inflammation.54 This study suggests that gulonate may contribute to ocular inflammation by modulating immune pathways, potentially influencing conditions such as AC. Furthermore, IL-18, a cytokine involved in immune responses, may modulate AC by regulating the levels of both 4-oxo-retinoic acid and gulonate, highlighting the complex interplay of metabolic pathways in ocular inflammation. 
Another MB, 3-(4-hydroxyphenyl) lactate, a product of tyrosine metabolism, is linked to the interaction between intestinal microbiota and host metabolism.55 This interaction not only influences metabolic disorders such as type 2 diabetes but also promotes systemic inflammation, which could contribute to AC pathogenesis.55,56 Additionally, TRANCE, a cytokine, has been shown to increase the risk of AC by promoting the accumulation of 3-(4-hydroxyphenyl) lactate, potentially exacerbating ocular inflammation. SAH, an intermediate in the methionine cycle, accumulates in response to oxidative stress and inflammation, both of which are critical to AC development.1,57,58 Elevated SAH levels are associated with enhanced inflammatory responses in the eye, with EN-RAGE promoting its accumulation and increasing AC risk. Furthermore, SM, an essential component of cell membranes involved in cell signaling, has been implicated in regulating inflammatory diseases.59 In the context of AC, alterations in SM metabolism may influence immune cell activation and ocular inflammation. CD244 may mitigate AC by reducing SM levels, thereby limiting inflammatory processes within the eye. 
By leveraging MR to infer causal relationships, the influence of confounding factors was reduced.17 The use of GWAS data enhanced the statistical power and robustness of our findings.26 Despite our stringent selection criteria, SNP selection may have inherent limitations. Measurement variability of MBs and CIPs could have introduced bias, and potential unmeasured confounders could have influenced our results. Additionally, potential false-positive results due to multiple hypothesis testing were not corrected by the false discovery rate. Although MR analysis provides evidence of causality, a combination of biological mechanism studies and clinical data is needed to fully understand its role. 
In conclusion, this study reveals significant causal relationships between CIPs and AC, with MBs playing mediating roles. These findings offer new perspectives on the pathophysiology of AC and suggest potential directions for future research and clinical applications. Further exploration of these pathways could inform the development of targeted interventions to improve the management and prevention of AC. 
Acknowledgments
The authors thank the investigators of the original studies for sharing the GWAS summary statistics. 
Supported by grants from the National Natural Science Foundation of China (82371033, 32200684), Tianjin Health Research Project (TJWJ2022QN078), and Tianjin Key Medical Discipline (Specialty) Construction Project (TJYXZDXK-016A). 
Author Contributions: ZJ and XY were responsible for the conception and design of the study. XC, ZX, QL, and BZ were responsible for the acquisition and analysis of data. XC and ZX wrote the paper. ZJ and XY reviewed the paper. All authors read and approved the final manuscript for submission. 
Disclosure: X. Cao, None; Z. Xu, None; B. Zhang, None; Q. Li, None; Z. Jiang, None; X. Yuan, None 
References
Leonardi A, Quintieri L, Presa IJ, et al. Allergic conjunctivitis management: update on ophthalmic solutions. Curr Allergy Asthma Rep. 2024; 24: 347–360. [CrossRef] [PubMed]
Abu SL, Hehar NK, Chigbu DI. Novel therapeutic receptor agonists and antagonists in allergic conjunctivitis. Curr Opin Allergy Clin Immunol. 2024; 24: 380–389. [CrossRef] [PubMed]
Miyazaki D, Fukagawa K, Okamoto S, et al. Epidemiological aspects of allergic conjunctivitis. Allergol Int. 2020; 69: 487–495. [CrossRef] [PubMed]
Chigbu DI, Karbach NJ, Abu SL, Hehar NK. Cytokines in allergic conjunctivitis: unraveling their pathophysiological roles. Life (Basel). 2024; 14: 350. [PubMed]
Falcon RMG, Caoili SEC. Immunologic, genetic, and ecological interplay of factors involved in allergic diseases. Front Allergy. 2023; 4: 1215616. [CrossRef] [PubMed]
Vazirani J, Shukla S, Chhawchharia R, et al. Allergic conjunctivitis in children: current understanding and future perspectives. Curr Opin Allergy Clin Immunol. 2020; 20: 507–515. [CrossRef] [PubMed]
Zou X, Huang H, Tan Y. Genetically determined metabolites in allergic conjunctivitis: a Mendelian randomization study. World Allergy Organ J. 2024; 17: 100894. [CrossRef] [PubMed]
Han X, Krempski JW, Nadeau K. Advances and novel developments in mechanisms of allergic inflammation. Allergy. 2020; 75: 3100–3111. [CrossRef] [PubMed]
Bielory L. Allergic conjunctivitis and the impact of allergic rhinitis. Curr Allergy Asthma Rep. 2010; 10: 122–134. [CrossRef] [PubMed]
Dhami S, Nurmatov U, Arasi S, et al. Allergen immunotherapy for allergic rhinoconjunctivitis: a systematic review and meta-analysis. Allergy. 2017; 72: 1597–1631. [CrossRef] [PubMed]
Chen Z, Suo Y, Du X, Zhao X. Genetically predicted N-methylhydroxyproline levels mediate the association between naive CD8+ T cells and allergic rhinitis: a mediation Mendelian randomization study. Front Immunol. 2024; 15: 1396246. [CrossRef] [PubMed]
Liu C, Chu D, Kalantar-Zadeh K, George J, Young HA, Liu G. Cytokines: from clinical significance to quantification. Adv Sci (Weinh). 2021; 8: e2004433. [CrossRef] [PubMed]
Marchev AS, Vasileva LV, Amirova KM, Savova MS, Balcheve-Sivenova ZP, Georgiev MI. Metabolomics and health: from nutritional crops and plant-based pharmaceuticals to profiling of human biofluids. Cell Mol Life Sci. 2021; 78: 6487–6503. [CrossRef] [PubMed]
James EL, Parkinson EK. Serum metabolomics in animal models and human disease. Curr Opin Clin Nutr Metab Care. 2015; 18: 478–483. [CrossRef] [PubMed]
Miyata J, Arita M. Role of omega-3 fatty acids and their metabolites in asthma and allergic diseases. Allergol Int. 2015; 64: 27–34. [CrossRef] [PubMed]
Li X, Meng X, Timofeeva M, et al. Serum uric acid levels and multiple health outcomes: umbrella review of evidence from observational studies, randomised controlled trials, and Mendelian randomisation studies. BMJ. 2017; 357: j2376. [PubMed]
Emdin CA, Khera AV, Kathiresan S. Mendelian randomization. JAMA. 2017; 318: 1925–1926. [CrossRef] [PubMed]
Zhang X, Yuan W, Xu J, Zhao F. Application of mendelian randomization in ocular diseases: a review. Hum Genomics. 2024; 18: 66. [CrossRef] [PubMed]
Birney E . Mendelian randomization. Cold Spring Harb Perspect Med. 2022; 12: a041302. [PubMed]
Zhou W, Cai J, Li Z, Lin Y. Association of atopic dermatitis with conjunctivitis and other ocular surface diseases: a bidirectional two-sample Mendelian randomization study [published online ahead of print March 13, 2023]. J Eur Acad Dermatol Venereol, https://doi.org/10.1111/jdv.19048.
Zhao JH, Stacey D, Eriksson N, et al. Genetics of circulating inflammatory proteins identifies drivers of immune-mediated disease risk and therapeutic targets. Nat Immunol. 2023; 24: 1540–1551. [CrossRef] [PubMed]
Chen Y, Lu T, Pettersson-Kymmer U, et al. Genomic atlas of the plasma metabolome prioritizes metabolites implicated in human diseases. Nat Genet. 2023; 55: 44–53. [CrossRef] [PubMed]
Slatkin M. Linkage disequilibrium—understanding the evolutionary past and mapping the medical future. Nat Rev Genet. 2008; 9: 477–485. [CrossRef] [PubMed]
Hemani G, Zheng J, Elsworth B, et al. The MR-Base platform supports systematic causal inference across the human phenome. eLife. 2018; 7: e34408. [CrossRef] [PubMed]
Burgess S, Thompson SG. Avoiding bias from weak instruments in Mendelian randomization studies. Int J Epidemiol. 2011; 40: 755–764. [CrossRef] [PubMed]
Carter AR, Sanderson E, Hammerton G, et al. Mendelian randomisation for mediation analysis: current methods and challenges for implementation. Eur J Epidemiol. 2021; 36: 465–478. [CrossRef] [PubMed]
Bowden J, Del Greco M F, Minelli C, et al. Improving the accuracy of two-sample summary-data Mendelian randomization: moving beyond the NOME assumption. Int J Epidemiol. 2019; 48: 728–742. [CrossRef] [PubMed]
Cornish AJ, Law PJ, Timofeeva M, et al. Modifiable pathways for colorectal cancer: a mendelian randomisation analysis. Lancet Gastroenterol Hepatol. 2020; 5: 55–62. [CrossRef] [PubMed]
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–525. [CrossRef] [PubMed]
Verbanck M, Chen C-Y, 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–698. [CrossRef] [PubMed]
Burgess S, Thompson SG. Interpreting findings from Mendelian randomization using the MR-Egger method. Eur J Epidemiol. 2017; 32: 377–389. [CrossRef] [PubMed]
Chai W, Zhang X, Lin M, et al. Allergic rhinitis, allergic contact dermatitis and disease comorbidity belong to separate entities with distinct composition of T-cell subsets, cytokines, immunoglobulins and autoantibodies. Allergy Asthma Clin Immunol. 2022; 18: 10. [CrossRef] [PubMed]
Zhang Y, Lan F, Zhang L. Advances and highlights in allergic rhinitis. Allergy. 2021; 76: 3383–3389. [CrossRef] [PubMed]
Kishimoto T. IL-6: from its discovery to clinical applications. Int Immunol. 2010; 22: 347–352. [CrossRef] [PubMed]
Sanders NL, Mishra A. Role of interleukin-18 in the pathophysiology of allergic diseases. Cytokine Growth Factor Rev. 2016; 32: 31–39. [CrossRef] [PubMed]
Ouyang W, O'Garra A. IL-10 family cytokines IL-10 and IL-22: from basic science to clinical translation. Immunity. 2019; 50: 871–891. [CrossRef] [PubMed]
Fukuda K, Kishimoto T, Sumi T, Yamashiro K, Ebihara N. Biologics for allergy: therapeutic potential for ocular allergic diseases and adverse effects on the eye. Allergol Int. 2023; 72: 234–244. [CrossRef] [PubMed]
Kaur D, Chachi L, Gomez E, Sylvius N, Brightling CE. Interleukin-18, IL-18 binding protein and IL-18 receptor expression in asthma: a hypothesis showing IL-18 promotes epithelial cell differentiation. Clin Transl Immunology. 2021; 10: e1301. [CrossRef] [PubMed]
Ek WE, Karlsson T, Höglund J, Rask-Andersen M, Johansson Å. Causal effects of inflammatory protein biomarkers on inflammatory diseases. Sci Adv. 2021; 7: eabl4359. [CrossRef] [PubMed]
Tao Z, Liu W, Chen Q, et al. Blocking Th2 signaling pathway alleviates the clinical symptoms and inflammation in allergic conjunctivitis. Invest Ophthalmol Vis Sci. 2023; 64: 30. [CrossRef] [PubMed]
Wang J, Chang C-Y, Yang X, et al. Leukemia inhibitory factor, a double-edged sword with therapeutic implications in human diseases. Mol Ther. 2023; 31: 331–343. [CrossRef] [PubMed]
Sun L, Gang X, Li Z, et al. Advances in understanding the roles of CD244 (SLAMF4) in immune regulation and associated diseases. Front Immunol. 2021; 12: 648182. [CrossRef] [PubMed]
Elishmereni M, Fyhrquist N, Singh Gangwar R, Lehtimäki S, Alenius H, Levi-Schaffer F. Complex 2B4 regulation of mast cells and eosinophils in murine allergic inflammation. J Invest Dermatol. 2014; 134: 2928–2937. [CrossRef] [PubMed]
El-Shazly AE, Henket M, Lefebvre PP, Louis R. 2B4 (CD244) is involved in eosinophil adhesion and chemotaxis, and its surface expression is increased in allergic rhinitis after challenge. Int J Immunopathol Pharmacol. 2011; 24: 949–960. [CrossRef] [PubMed]
Hofmann MA, Drury S, Fu C, et al. RAGE mediates a novel proinflammatory axis: a central cell surface receptor for S100/calgranulin polypeptides. Cell. 1999; 97: 889–901. [CrossRef] [PubMed]
Hudson BI, Lippman ME. Targeting RAGE signaling in inflammatory disease. Annu Rev Med. 2018; 69: 349–364. [CrossRef] [PubMed]
Kwapisz L, Gregor J, Chande N, Yan B, Ponich T, Mosli M. The utility of fecal calprotectin in predicting the need for escalation of therapy in inflammatory bowel disease. Scand J Gastroenterol. 2017; 52: 846–850. [CrossRef] [PubMed]
Sumsion JS, Pulsipher A, Alt JA. Differential expression and role of S100 proteins in chronic rhinosinusitis. Curr Opin Allergy Clin Immunol. 2020; 20: 14–22. [CrossRef] [PubMed]
Bar N, Korem T, Weissbrod O, et al. A reference map of potential determinants for the human serum metabolome. Nature. 2020; 588: 135–140. [CrossRef] [PubMed]
Richards JL, Yap YA, McLeod KH, Mackay CR, Mariño E. Dietary metabolites and the gut microbiota: an alternative approach to control inflammatory and autoimmune diseases. Clin Transl Immunology. 2016; 5: e82. [CrossRef] [PubMed]
Larange A, Cheroutre H. Retinoic acid and retinoic acid receptors as pleiotropic modulators of the immune system. Annu Rev Immunol. 2016; 34: 369–394. [CrossRef] [PubMed]
Lotfi R . Retinoic acid (RA): a critical immunoregulatory molecule in asthma and allergies. Immun Inflamm Dis. 2024; 12: e70051. [CrossRef] [PubMed]
Liang L, He M, Zhang Y, et al. Unraveling the 2,3-diketo-l-gulonic acid-dependent and -independent impacts of l-ascorbic acid on somatic cell reprogramming. Cell Biosci. 2023; 13: 218. [CrossRef] [PubMed]
Liu Z-Y, Zhang H, Sun X-L, Liu J-Y. Causal association between metabolites and age-related macular degeneration: a bidirectional two-sample mendelian randomization study. Hereditas. 2024; 161: 51. [CrossRef] [PubMed]
Wikoff WR, Anfora AT, Liu J, et al. Metabolomics analysis reveals large effects of gut microflora on mammalian blood metabolites. Proc Natl Acad Sci U S A. 2009; 106: 3698–3703. [CrossRef] [PubMed]
Vangipurapu J, Fernandes Silva L, Kuulasmaa T, Smith U, Laakso M. Microbiota-related metabolites and the risk of type 2 diabetes. Diabetes Care. 2020; 43: 1319–1325. [CrossRef] [PubMed]
Parkhitko AA, Jouandin P, Mohr SE, Perrimon N. Methionine metabolism and methyltransferases in the regulation of aging and lifespan extension across species. Aging Cell. 2019; 18: e13034. [CrossRef] [PubMed]
Xiao Y, Su X, Huang W, et al. Role of S-adenosylhomocysteine in cardiovascular disease and its potential epigenetic mechanism. Int J Biochem Cell Biol. 2015; 67: 158–166. [CrossRef] [PubMed]
Adada M, Luberto C, Canals D. Inhibitors of the sphingomyelin cycle: sphingomyelin synthases and sphingomyelinases. Chem Phys Lipids. 2016; 197: 45–59. [CrossRef] [PubMed]
Figure 1.
 
Schematic diagram of this study. (A) The total effect c was analyzed using genetically predicted CIPs as the exposure and AC as the outcome. Additionally, the total effect d was also analyzed using genetically predicted AC as the exposure and CIPs as the outcome. The relevance and independence of IVs were considered for both exposure and outcome, and potential confounders were accounted for. (B) The total effect was further broken down as follows: The indirect effect was assessed using a two-step approach involving the effects of CIPs on MBs (a) and the effects of MBs on AC (b), as well as the product method (a × b), and the direct effect was calculated as c′ = ca × b. The proportion mediated was determined by dividing the indirect effect by the total effect.
Figure 1.
 
Schematic diagram of this study. (A) The total effect c was analyzed using genetically predicted CIPs as the exposure and AC as the outcome. Additionally, the total effect d was also analyzed using genetically predicted AC as the exposure and CIPs as the outcome. The relevance and independence of IVs were considered for both exposure and outcome, and potential confounders were accounted for. (B) The total effect was further broken down as follows: The indirect effect was assessed using a two-step approach involving the effects of CIPs on MBs (a) and the effects of MBs on AC (b), as well as the product method (a × b), and the direct effect was calculated as c′ = ca × b. The proportion mediated was determined by dividing the indirect effect by the total effect.
Figure 2.
 
Flowchart of data source and SNP selection process for exposure (CIPs), mediator (MBs), and outcome (AC). Linkage disequilibrium, LD.
Figure 2.
 
Flowchart of data source and SNP selection process for exposure (CIPs), mediator (MBs), and outcome (AC). Linkage disequilibrium, LD.
Figure 3.
 
Bidirectional MR analysis between CIPs and AC. (A) The circular heatmap of P values illustrates the causal relationship between CIPs as exposures and AC as the outcome using various MR methods. The vertical axis represents six CIPs (TRANCE, LIF, IL-6, IL-18R1, EN-RAGE, and CD244), and the horizontal axis, from outermost to innermost, represents five MR methods (weighted mode, simple mode, IVW, weighted median, and MR Egger). Lower P values provide stronger evidence for a causal relationship; the IVW method identified six CIPs with significant causal associations with AC. (B) The circular heatmap of OR visualizes the strength and direction of causal associations between CIPs as exposures and AC as the outcome. Six CIPs showed significant associations, with OR analysis identifying both protective and risk factors. TRANCE, LIF, IL-6, IL-18R1, and EN-RAGE were risk factors, and CD244 was a protective factor. (CH) Scatterplots of causal effects of SNPs associated with TRANCE (C), EN-RAGE (D), IL-18R1 (E), CD244 (F), LIF (G), and IL-6 (H) on AC. These plots demonstrate individual causal relationships and provide visual evidence of consistency across different SNPs. (I) The circular heatmap of P values displays the results of reverse MR analysis with AC as the exposure and CIPs as the outcome. No significant results were observed, indicating the absence of a reverse causal relationship from AC to the six CIPs. (J) Schematic representation of the bidirectional MR analysis results that summarizes the significant causal relationships identified and highlights the lack of significant reverse associations.
Figure 3.
 
Bidirectional MR analysis between CIPs and AC. (A) The circular heatmap of P values illustrates the causal relationship between CIPs as exposures and AC as the outcome using various MR methods. The vertical axis represents six CIPs (TRANCE, LIF, IL-6, IL-18R1, EN-RAGE, and CD244), and the horizontal axis, from outermost to innermost, represents five MR methods (weighted mode, simple mode, IVW, weighted median, and MR Egger). Lower P values provide stronger evidence for a causal relationship; the IVW method identified six CIPs with significant causal associations with AC. (B) The circular heatmap of OR visualizes the strength and direction of causal associations between CIPs as exposures and AC as the outcome. Six CIPs showed significant associations, with OR analysis identifying both protective and risk factors. TRANCE, LIF, IL-6, IL-18R1, and EN-RAGE were risk factors, and CD244 was a protective factor. (CH) Scatterplots of causal effects of SNPs associated with TRANCE (C), EN-RAGE (D), IL-18R1 (E), CD244 (F), LIF (G), and IL-6 (H) on AC. These plots demonstrate individual causal relationships and provide visual evidence of consistency across different SNPs. (I) The circular heatmap of P values displays the results of reverse MR analysis with AC as the exposure and CIPs as the outcome. No significant results were observed, indicating the absence of a reverse causal relationship from AC to the six CIPs. (J) Schematic representation of the bidirectional MR analysis results that summarizes the significant causal relationships identified and highlights the lack of significant reverse associations.
Figure 4.
 
Forest plot to visualize the causal effects of MBs with CIPs and AC.
Figure 4.
 
Forest plot to visualize the causal effects of MBs with CIPs and AC.
×
×

This PDF is available to Subscribers Only

Sign in or purchase a subscription to access this content. ×

You must be signed into an individual account to use this feature.

×