# Leukotriene receptor antagonists and eosinophilic granulomatosis with polyangiitis: a disproportionality analysis from FAERS, JADER, CVAR databases integrated with network pharmacology

**Authors:** Cuilv Liang, Pingping Zhuo, Yin Zhang, Vipula Bataduwaarachchi, Vipula Bataduwaarachchi, Vipula Bataduwaarachchi

PMC · DOI: 10.1371/journal.pone.0343084 · PLOS One · 2026-03-09

## TL;DR

This study finds that leukotriene receptor antagonists may increase the risk of a rare autoimmune disease called EGPA and identifies potential genetic mechanisms behind this association.

## Contribution

The study is the first to integrate pharmacovigilance data with network pharmacology to explore the risk and mechanisms of LTRAs-induced EGPA.

## Key findings

- LTRAs like montelukast are strongly associated with EGPA across multiple databases.
- Montelukast remains significantly linked to EGPA even when corticosteroids are excluded.
- Key genes such as SRC, PTGS2, and EDN1 are identified as central to LTRAs-induced EGPA.

## Abstract

The relationship between leukotriene receptor antagonists (LTRAs) usage and the subsequent occurrence of eosinophilic granulomatosis with polyangiitis (EGPA) remained highly polarizing and controversial in previous studies. We aimed to investigate the risk of EGPA caused by LTRAs and the potential toxicological mechanisms of LTRAs-related EGPA.

In this real-world pharmacovigilance study, we collected adverse event (AE) reports of EGPA associated with LTRAs use from the U.S. FDA Adverse Event Reporting System (FAERS), Japanese Adverse Drug Event Reporting (JADER), and Canadian Vigilance Adverse Reaction (CVAR) databases. The reporting odds ratio (ROR), proportional reporting ratio (PRR), information component (IC), and empirical Bayesian geometric mean (EBGM) were calculated to quantify the strength of the association between LTRAs and EGPA. The Weibull shape parameter (WSP) test was applied to analyze time-to-onset profiles of EGPA toxicity. Network pharmacology analysis was subsequently performed to identify the central genes to determine the potential mechanisms underlying LTRAs-induced EGPA.

LTRAs, including montelukast, zafirlukast, and pranlukast, exhibited a strong association with EGPA in three databases (the lower limit of 95% confidence interval (CI) for ROR > 1, PRR > 2 with χ2 values ≥4, EBGM05 > 2, and IC025 > 0). After excluding corticosteroids as concomitant medication, montelukast remained significantly associated with EGPA in the FAERS database. The median time-to-onset of EGPA associated with LTRAs was 233 (range: 76–660) days, and the WSP test indicated LTRAs had early failure-type profiles. We isolated 81 interactive target genes linking LTRAs to EGPA. Several central genes, including SRC, PTGS2, EDN1, HMOX1, KDR, and OCLN, were revealed via protein-protein interactions analysis and molecular complex detection (MCODE) algorithm.

Our study revealed LTRAs could increase the risk of EGPA, and initially explored potential genes and mechanisms of LTRAs-induced EGPA. It is helpful for clinicians to be alerted to the risk of EGPA during LTRAs administration.

## Linked entities

- **Genes:** SRC (SRC proto-oncogene, non-receptor tyrosine kinase) [NCBI Gene 6714], PTGS2 (prostaglandin-endoperoxide synthase 2) [NCBI Gene 5743], EDN1 (endothelin 1) [NCBI Gene 1906], HMOX1 (heme oxygenase 1) [NCBI Gene 3162], KDR (kinase insert domain receptor) [NCBI Gene 3791], OCLN (occludin) [NCBI Gene 100506658]
- **Chemicals:** montelukast (PubChem CID 5281040), zafirlukast (PubChem CID 5717), pranlukast (PubChem CID 4887)
- **Diseases:** eosinophilic granulomatosis with polyangiitis (MONDO:0015943), EGPA (MONDO:0015943)

## Full-text entities

- **Genes:** EDN1 (endothelin 1) [NCBI Gene 1906] {aka ARCND3, ET1, HDLCQ7, PPET1, QME}, CYSLTR2 (cysteinyl leukotriene receptor 2) [NCBI Gene 57105] {aka CYSLT2, CYSLT2R, GPCR21, HG57, HPN321, KPG_011}, CYSLTR1 (cysteinyl leukotriene receptor 1) [NCBI Gene 10800] {aka CYSLT1, CYSLT1R, CYSLTR, HMTMF81}, HSD17B6 (hydroxysteroid 17-beta dehydrogenase 6) [NCBI Gene 8630] {aka HSE, RODH, SDR9C6}, REN (renin) [NCBI Gene 5972] {aka ADTKD4, HNFJ2, RTD}, VEGFA (vascular endothelial growth factor A) [NCBI Gene 7422] {aka L-VEGF, MVCD1, VEGF, VPF}, OCLN (occludin) [NCBI Gene 100506658] {aka BLCPMG, PPP1R115, PTORCH1}, PTGS2 (prostaglandin-endoperoxide synthase 2) [NCBI Gene 5743] {aka COX-2, COX2, GRIPGHS, PGG/HS, PGHS-2, PHS-2}, GATA3 (GATA binding protein 3) [NCBI Gene 2625] {aka HDR, HDRS}, SRC (SRC proto-oncogene, non-receptor tyrosine kinase) [NCBI Gene 6714] {aka ASV, SRC1, THC6, c-SRC, p60-Src}, LPP (LIM domain containing preferred translocation partner in lipoma) [NCBI Gene 4026], KDR (kinase insert domain receptor) [NCBI Gene 3791] {aka CD309, FLK1, VEGFR, VEGFR2}, BACH2 (BACH transcriptional regulator 2) [NCBI Gene 60468] {aka BTBD25, IMD60}, HLA-DRB4 (major histocompatibility complex, class II, DR beta 4) [NCBI Gene 3126] {aka DR4, DRB4, HLA-DR4B, HLA-DRB, HLA-DRB4*}, HMOX1 (heme oxygenase 1) [NCBI Gene 3162] {aka HMOX1D, HO-1, HSP32, bK286B10}, TSLP (thymic stromal lymphopoietin) [NCBI Gene 85480], GPR17 (G protein-coupled receptor 17) [NCBI Gene 2840], LINC-ROR (long intergenic non-protein coding RNA, regulator of reprogramming) [NCBI Gene 100885779] {aka ROR, lincRNA-RoR, lincRNA-ST8SIA3}
- **Diseases:** immune-mediated disorder (MESH:C567355), chronic obstructive pulmonary disease (MESH:D029424), Churg-Strauss syndrome (MESH:D015267), vasculitis (MESH:D014657), inflammation (MESH:D007249), Asthma (MESH:D001249), eosinophilia (MESH:D004802), diabetes (MESH:D003920), EBGM (MESH:D009800), allergic conditions (MESH:D004342), SRSs (MESH:D005598), LTRAs (MESH:C565439), rhinitis (MESH:D012220), eosinophilic granulomatosis (MESH:D017681), EGPA (MESH:D014890), hypertension (MESH:D006973), WSP (MESH:C562399), melanoma metastasis (MESH:D009362), underlying (MESH:C566142), AE (MESH:D064420), Reaction (MESH:D006967), infections (MESH:D007239)
- **Chemicals:** prostaglandin E2 (MESH:D015232), thiouracil (MESH:D013889), arachidonic acid (MESH:D016718), silicone (MESH:D012828), ibudilast (MESH:C038366), torsemide (MESH:D000077786), prostaglandins (MESH:D011453), zafirlukast (MESH:C062735), pranlukast (MESH:C047681), minocycline (MESH:D008911), steroid (MESH:D013256), prednisone (MESH:D011241), omalizumab (MESH:D000069444), budesonide (MESH:D019819), NAD (MESH:D009243), Montelukast (MESH:C093875), EGFA (-), eicosanoid (MESH:D015777), fluticasone (MESH:D000068298), NADP (MESH:D009249), salmeterol (MESH:D000068299), macrolides (MESH:D018942), hydralazine (MESH:D006830), CysLTs (MESH:C112381), Leukotrienes (MESH:D015289)
- **Species:** Homo sapiens (human, species) [taxon 9606]
- **Cell lines:** Line 419 — Homo sapiens (Human), Finite cell line (CVCL_X131)

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12970897/full.md

## Figures

5 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12970897/full.md

## References

53 references — full list in the complete paper: https://tomesphere.com/paper/PMC12970897/full.md

---
Source: https://tomesphere.com/paper/PMC12970897