Skip to the main content

Original scientific paper

https://doi.org/10.2478/acph-2026-0022

Adverse event reporting of golimumab: A disproportionality analysis of the FDA Adverse Event Reporting System (FAERS) database

WENJIE MA ; Department of Pharmacy, Southern University of Science and Technology Hospital, Shenzhen, 518055, Guangdong, China
NUAN JIA ; Department of Pharmacy, Southern University of Science and Technology Hospital, Shenzhen, 518055, Guangdong, China
HUI LIAO ; Department of Pharmacy, Southern University of Science and Technology Hospital, Shenzhen, 518055, Guangdong, China
ZHILEI LI ; Department of Pharmacy, Southern University of Science and Technology Hospital, Shenzhen, 518055, Guangdong, China
XIAODAN LI ; Department of Pharmacy, People’s Hospital of Longhua, Shenzhen, 518109, Guangdong, China *

* Corresponding author.


Full text: english pdf 3.465 Kb

downloads: 0

cite

Download JATS file


Abstract

This study conducted a retrospective disproportionality analysis of 50,031 adverse drug event (ADE) reports from the FDA Adverse Event Reporting System (FAERS) database spanning Q2 2009 to Q3 2025, in which golimumab was identified as the primary suspect drug. Four disproportionality algorithms were simultaneously applied: the reporting odds ratio (ROR), proportional reporting ratio (PRR), Bayesian confidence propagation neural network (BCPNN), and multi-item gamma Poisson shrinker (MGPS). A total of 397 positive signals were identified across 25 system organ classes (SOCs). The SOCs with concordant positive signals across all four algorithms were identified for infections and infestations and surgical and medical procedures. High-frequency signals included pneumonia, lower respiratory tract infection, influenza, nasopharyngitis, sinusitis and cellulitis. Comparison with the current prescribing information revealed several previously undocumented ADEs, including urinary tract infection, ear infection, flavivirus infection, Zika virus infection and Dengue fever. Twenty-nine positive malignancy signals absent from the prescribing information were identified, including basal cell carcinoma, colorectal cancer, cervical cancer and breast cancer. Time-to-onset analysis revealed that the risk of ADEs was the highest within the first few months of treatment initiation. Clinicians should rigorously monitor for infectious complications and malignancies when prescribing golimumab to ensure patient safety.

Keywords

golimumab; FAERS; ADE; signal mining

Hrčak ID:

350534

URI

https://hrcak.srce.hr/350534

Publication date:

30.9.2026.

Visits: 0 *




INTRODUCTION

Tumour necrosis factor-alpha (TNF-α) inhibitors are a class of biologic agents used in the treatment of various immune-mediated inflammatory diseases, exerting their therapeutic effects by blocking TNF-α activity and downregulating downstream inflammatory signalling pathways (1). Golimumab is a fully human anti-TNF-α monoclonal antibody that specifically binds to both soluble and transmembrane biologically active forms of TNF-α, forming a complex that blocks the binding of TNF-α to its receptor and thereby inhibits the biological activity of TNF-α (2). It is currently approved primarily for the treatment of rheumatoid arthritis (RA) and ankylosing spondylitis (AS) (3).

Golimumab received its initial marketing authorisation from the United States Food and Drug Administration (FDA) in 2009, followed by approval from the European Medicines Agency (EMA) in the same year, and entered the Chinese market in 2018 (4). Prior to approval, the FDA incorporated black box warnings for serious infections and malignancies into the drug's prescribing information; these warnings have been revised on subsequent occasions, though the scope of updates may not fully reflect all known risks, and the safety profile of golimumab warrants continued attention. Multiple clinical trials conducted both within China and internationally have evaluated the safety of golimumab in RA and AS, identifying numerous adverse drug events (ADEs) that were undetected during pre-approval studies (5–7). Among these, serious adverse events (SAEs) accounted for a certain proportion, and fatal cases have been reported (8). These findings not only reveal the limitations of pre-marketing safety evaluation but also highlight the necessity of conducting post-marketing surveillance based on real-world data. The FDA Adverse Event Reporting System (FAERS), the world's largest spontaneous adverse event reporting database, is widely used in post‑marketing drug safety surveillance and can provide valuable real‑world evidence to support clinical decision‑making.

This study is the first to draw on the FAERS database and jointly apply four pharmacovigilance disproportionality analysis algorithms, namely reporting odds ratio (ROR), proportional reporting ratio (PRR), Bayesian confidence propagation neural network (BCPNN), and multi-item gamma Poisson shrinker (MGPS), to systematically evaluate the real-world safety profile of golimumab, thereby providing clinicians with more reliable evidence‑based data for appropriate drug selection and prevention of adverse drug events (9).

EXPERIMENTAL

Data source and pre-processing

This retrospective analysis draws on quarterly raw data extracted from the FAERS database (https://fis.fda.gov/extensions/FPD-QDE-FAERS/FPD-QDE-FAERS.html). The FAERS dataset comprises seven sections: DEMO (patient demographic and administrative information), DRUG (drug-specific details), REAC (coded adverse events), OUTC (patient outcomes), RPSR (report sources), THER (therapy initiation and cessation dates), and INDI (indications for drug use). Reported drugs in FAERS are categorised into four groups: PS (primary suspect), SS (secondary suspect), C (concomitant), and I (interacting) (9). To ensure comprehensive retrieval, both brand names and generic names for golimumab preparations were used as search terms, including “Simponi," "golimumab," and "Simponi Aria." The query covered reports from the second quarter of 2009 through the third quarter of 2025, and only records in which golimumab was designated as the primary suspect drug were included for analysis.

To address duplicate reports, we adopted the methodology recommended by the FDA. We selected the PRIMARYID, CASEID, and FDA_DT fields from the DEMO table and sorted them by CASEID, FDA_DT, and PRIMARYID. We retained the report with the largest FDA_DT value for cases with the same CASEID. Subsequently, for cases with both the same CASEID and FDA_DT, we kept the report with the largest PRIMARYID value. A list of deleted reports was included in each quarterly data package since Q2 2009, and after data deduplication, reports were excluded based on the CASEID present in the list of deleted reports (10). All preferred terms (PTs) and system organ classes (SOCs) were coded and categorised using the Medical Dictionary for Regulatory Activities (MedDRA), version 28.0 (11).

Statistical analysis

Disproportionality analysis was employed to quantitatively evaluate golimumab-associated ADE signals. This analysis typically comprises two components: frequency‑based statistics and Bayesian statistics. Frequency‑based statistics include the ROR and PRR, whereas Bayesian statistics encompass the MGPS and BCPNN. Although the frequency‑based ROR and PRR methods offer high sensitivity, computational simplicity, and earlier detection of certain signals, they are prone to false‑positive signals when the number of events is small; therefore, their standalone application is susceptible to the masking effect (9). In contrast, BCPNN and MGPS demonstrate superior performance in detecting low‑frequency events and can reduce false‑positive signals compared with the frequency‑based methods. The integrated use of all four methods leverages their complementary strengths and minimises biases inherent to any single approach (10).

In this study, a signal was considered positive only when the positivity thresholds of all four algorithms were simultaneously met, indicating the generation of an ADE risk signal. The positivity criterion for ROR was defined as a 95 % confidence interval lower limit > 1 with at least three reports; for PRR, as PRR ≥ 2 with at least three reports and χ² value ≥ 4; for BCPNN, as IC025 > 0; and for MGPS, as EBGM05 > 2 (9). The 2×2 contingency table used in the descriptive analysis is presented in Table S1, and the specific formulas and positivity thresholds for all four algorithms are provided in Table S2.

Time to onset (TTO) was defined as the interval between the date of ADE occurrence (from the DEMO file) and the date of treatment initiation (from the THER file). Cases in which the ADE occurred before therapy initiation, as well as those with inaccurate or missing key information, were excluded. The Weibull distribution test was applied to evaluate the temporal trend of ADE incidence, allowing for the determination of whether the risk demonstrated an increasing or decreasing tendency (12). Data processing was conducted using MySQL; data visualisation was performed using the "matplotlib" package in Python 3.6.

RESULTS AND DISCUSSION

Basic characteristics of ADE reports

As shown in Fig. 1, the flowchart illustrates the complete process of data extraction and analysis. The initial dataset comprised 21,757,413 records. After removing 3,498,215 duplicate entries, we obtained 18,259,198 unique demographic records (DEMO). Subsequently, after further filtering of drug records (DRUG) and adverse reaction records (REAC) associated with golimumab, a total of 50,031 reports of golimumab-induced adverse events were ultimately identified, involving 123,969 adverse event cases.

image1.jpeg

Fig. 1. Flowchart of the screening process for golimumab-related adverse drug events (ADEs).

From the second quarter of 2009 through the fourth quarter of 2015, the number of adverse event reports increased overall, peaking at 5,098 cases in 2015. After that, the number remained high, with annual totals exceeding 3,000 cases every year from 2016 to 2025 (Fig. S1). This reflects the widespread use of golimumab in real-world clinical practice. Demographic data are shown in Table I. Among reports with complete demographic data, the male-to-female ratio was 0.4:1, which may reflect a higher prevalence of the drug's indications in female patients (13). The most frequent reporter occupations were nurses (21,953 reports, 43.9 %) and physicians (8,489 reports, 17.0 %). Approximately three-quarters of the reports were submitted by healthcare professionals, suggesting good reliability of the reported events. The most common reported outcomes were hospitalisation or prolonged hospitalisation (9,416 reports, 18.8 %) and other serious (important medical event) (21,477 reports, 42.9 %). The majority of reports originated from the United States (17,527 reports, 35.0 %), followed by Canada (7,844 reports, 15.7 %), the United Kingdom (4,886 reports, 9.8 %), and Brazil (4,371 reports, 8.7 %). The geographic distribution characteristics reported herein may be influenced by inter-regional differences in drug utilisation and spontaneous reporting practices.

Table I. Demographic characteristics of ADEs reported in the FAERS database with golimumab as the primary suspect drug

VariableCase numbers (n)Percentage (%)
SexMale13,08726.2
Female31,94863.9
Unknown4,99610.0
Age (years)≤ 172540.5
18–6421,75443.5
65–859,35218.7
> 854580.9
Unknown18,21336.4
ReporterRegistered nurse21,95343.9
Physician8,48917.0
Pharmacist5,45510.9
Lawyer7< 0.1
Other health professional7,38714.8
Unknown6,74013.5
OutcomesDeath1,1362.3
Disability2600.5
Life-threatening3110.6
Hospitalization9,41618.8
Congenital anomaly27< 0.1
Other serious21,47742.9
Unknown17,40434.8
Reported countries (top five)United States17,52735.0
Canada7,84415.7
United Kingdom4,8869.8
Brazil4,3718.7
Japan2,9215.8
Other12,48225.0

Signal detection at the system organ class level

After excluding drug-unrelated MedDRA SOC categories (social circumstances and product issues), 25 SOCs remained, ranked by report frequency. The three most commonly affected SOCs were Infections and infestations (21.9 %), Injury, poisoning and procedural Complications (15.8 %), and General disorders and administration site conditions (14.9 %), as shown in Fig. S2. Positive signals meeting the criteria of all four algorithms were identified for Infections and infestations (ROR: 4.25 [4.19–4.31]; PRR: 3.63 [χ²: 47,114]; EBGM05: 3.56; IC025: 1.83) and Surgical and medical procedures (ROR: 2.39 [2.31–2.46]; PRR: 2.34 [χ²: 3,107]; EBGM05: 2.26; IC025: 1.18). Injury, poisoning and procedural complications (ROR: 1.45 [1.42–1.47]; IC025: 0.45), Musculoskeletal and connective tissue disorders (ROR: 1.51 [1.46–1.53]; IC025: 0.51), and Neoplasms benign, malignant and unspecified (ROR: 1.27 [1.23–1.31]; IC025: 0.29) exceeded the positivity threshold in the ROR and BCPNN methods only, as shown in Table II.

Table II. Signal strength of golimumab ADEs across system organ classes (SOC) in the FAERS database

SOCCase reportsROR (95 % CI)PRR (χ2)EBGM (EBGM05)IC (IC025)
Infections and infestations23,6004.25 (4.19–4.31)3.63 (47,114)3.61 (3.56)1.85 (1.83)
Injury, poisoning and procedural complications17,0621.45 (1.42–1.47)1.39 (2,033)1.39 (1.36)0.47 (0.45)
General disorders and administration site conditions16,1130.71 (0.7–0.72)0.75 (1,669)0.75 (0.74)–0.42 (–0.44)
Musculoskeletal and connective tissue disorders9,3251.51 (1.46–1.53)1.46 (1,412)1.46 (1.43)0.54 (0.51)
Gastrointestinal disorders6,3580.59 (0.58–0.61)0.61 (1,697)0.61 (0.60)–0.71 (–0.74)
Nervous system disorders4,6350.44 (0.43–0.45)0.46 (3,207)0.46 (0.45)–1.12 (–1.16)
Skin and subcutaneous tissue disorders4,1240.60 (0.58–0.62)0.61 (1,063)0.61 (0.59)–0.71 (–0.75)
Neoplasms benign, malignant and unspecified (incl cysts and polyps)4,1051.27 (1.23–1.31)1.26 (226.5)1.26 (1.22)0.33 (0.29)
Surgical and medical procedures4,0032.39 (2.31–2.46)2.34 (3,107)2.34 (2.26)1.22 (1.18)
Respiratory, thoracic and mediastinal disorders3,5980.61 (0.59–0.63)0.62 (858.3)0.62 (0.60)–0.68 (–0.73)
Investigations3,1540.43 (0.41–0.44)0.44 (2,350)0.44 (0.43)–1.17 (–1.23)
Vascular disorders1,6930.66 (0.63–0.69)0.67 (289.8)0.67 (0.63)–0.59 (–0.66)
Psychiatric disorders1,5430.22 (0.21–0.23)0.23 (4,219)0.23 (0.22)–2.12 (–2.19)
Cardiac disorders1,4320.48 (0.45–0.50)0.48 (806.1)0.48 (0.46)–1.05 (–1.12)
Eye disorders1,4150.56 (0.54–0.59)0.57 (471.2)0.57 (0.54)–0.81 (–0.89)
Renal and urinary disorders1,1170.50 (0.47–0.53)0.50 (567.1)0.51 (0.47)–1.00 (–1.08)
Immune system disorders1,0050.73 (0.68–0.77)0.73 (103.6)0.73 (0.68)–0.46 (–0.55)
Metabolism and nutrition disorders9370.36 (0.34–0.38)0.36 (1,073)0.36 (0.34)–1.46 (–1.56)
Blood and lymphatic system disorders7240.35 (0.33–0.38)0.36 (860.8)0.36 (0.33)–1.49 (–1.60)
Hepatobiliary disorders7050.65 (0.60–0.70)0.65 (133.5)0.65 (0.60)–0.62 (–0.73)
Reproductive system and breast disorders4560.46 (0.42–0.5)0.46 (289.8)0.46 (0.42)–1.12 (–1.25)
Ear and labyrinth disorders3740.71 (0.65–0.79)0.71 (42.7)0.72 (0.65)–0.48 (–0.63)
Pregnancy, puerperium and perinatal conditions1760.36 (0.31–0.41)0.36 (205.6)0.36 (0.31)–1.49 (–1.70)
Endocrine disorders1640.52 (0.45–0.61)0.52 (72.4)0.52 (0.45)–0.94 (–1.16)
Congenital, familial and genetic disorders900.26 (0.21–0.32)0.26 (188.9)0.26 (0.21)–1.94 (–2.23)

ADEs – adverse drug events, SOC – system organ class, ROR – reporting odds ratio, 95 % CI – 95 % confidence interval, PRR – proportional reporting ratio, χ2 – chi-squared, EBGM – empirical Bayesian geometric mean, EBGM05 – the lower limit of the 95 % CI of EBGM, IC – information component, IC025 – the lower limit of the 95 % CI of the IC.

The predominance of Infections and infestations is consistent with the known safety profile of TNF-α inhibition regarding infectious risks and supports further evaluation of infection-related events at the preferred term (PT) level. In contrast, the concordant signal for Surgical and medical procedures should be interpreted with caution, as this SOC contains procedure-related terms that may partially reflect underlying disease burden or clinical management rather than a direct pharmacological adverse effect (14).

Signal detection at the preferred term level

A total of 397 PT-level ADE signals simultaneously met the positivity criteria of all four algorithms, distributed across 25 SOCs (Fig. 2). Fig. 3 presents the top 50 adverse drug events associated with golimumab at the preferred term (PT) level, ranked by frequency of reports. The adverse events already documented in the prescribing information include pneumonia (n = 3,047), lower respiratory tract infection (n = 1,743), influenza (n = 1,121), infection (n = 1,030), nasopharyngitis (n = 847), sinusitis (n = 668), cellulitis (n = 632), bronchitis (n = 553), and herpes zoster (n = 483). The most frequently reported positive infection signals not listed on the label were urinary tract infection (n = 1,334), kidney infection (n = 299), ear infection (n = 252), and tooth abscess (n = 187). Serious infections associated with golimumab are subject to a black box warning in the prescribing information. A phase III randomised, placebo-controlled trial in patients with active psoriatic arthritis reported serious infections, including periodontitis, urinary tract infection, and ear infection, during intravenous golimumab treatment through week 60 (15). According to the literature, TNF-α is indispensable for the activation, differentiation, and recruitment of multiple immune cell populations, plays a pivotal role in maintaining granuloma integrity, and host defense against mycobacteria and intracellular pathogens. Consequently, TNF-α blockade substantially elevates the risk of bacterial, viral, and fungal infections (16, 17). In summary, these findings further underscore the importance of systematic infection risk management during golimumab therapy, encompassing comprehensive pre-treatment evaluation of active and latent infections, patient education, and vigilant clinical and laboratory surveillance during and after treatment (18). image2.png

Fig. 2. Venn plot of the number of positive adverse drug event reports for different methods.

image3.png

Fig. 3. The top 50 ADEs of golimumab ranked by case reports at the PTs level. SOC – system organ class, PTs – preferred terms, ROR – reporting odds ratio, 95 % CI – 95 % confidence interval, PRR – proportional reporting ratio, EBGM05 – the lower limit of the 95 % CI for EBGM, IC025 – the lower limit of the 95 % CI for IC, ADEs – adverse drug events.

Malignant events associated with golimumab are also subject to a black-box warning in the prescribing information. In this study, four malignancy events documented in the label with relatively high reporting frequencies were identified, namely malignant neoplasm (n = 416), skin cancer (n = 190), lymphoma (n = 155), and malignant melanoma (n = 148). Several clinical studies and meta-analyses have suggested that TNF-α inhibitors may increase the risk of certain types of tumours (19–22). Specifically, golimumab has been associated with an elevated risk of skin cancer, particularly non-melanoma skin cancer. Furthermore, TNF-α inhibitors may also increase the risk of lymphoma, with higher risks observed in patients with long-standing chronic inflammation or those receiving high-dose therapy (21–22). The potential mechanism underlying malignancy events may be attributable to the physiological role of TNF-α in suppressing tumour growth within the normal immune system; TNF-α inhibitors block the binding of TNF-α to its receptor, thereby inhibiting its biological activity. Therefore, caution should be exercised when using this agent.

Fig. 4 shows the 30 highest-ranked positive signals by strength, of which 11 are listed in the prescribing information, and 19 are not. The three signals with the highest ROR were administration site infection (n = 23; ROR: 143.87 [89.75–230.64]; PRR: 143.85 [χ2: 2,447]; EBGM05: 67.45; IC025: 3.64), flavivirus infection (n = 24; ROR: 100.57 [64.49–156.83]; PRR: 100.55 [χ2: 1,918]; EBGM05: 52.42; IC025: 3.64), and Zika virus infection (n = 13; ROR: 89.06 [49.02–161.80]; PRR: 89.05 [χ2: 938]; EBGM05: 40.72; IC025: 2.73), all of which are absent from the prescribing information. Lower respiratory tract infection (n =1,743; ROR: 20.62 [19.65–21.64]; PRR: 20.34 [χ2: 30,640]; EBGM05: 18.55; IC025: 4.20) and Dengue fever (n = 171; ROR: 37.84 [32.36–44.24]; PRR: 37.78 [χ2: 5,631]; EBGM05: 29.78; IC025: 4.63) were noteworthy for combining both high report frequency and strong signal strength. It is noteworthy that the genus Flavivirus, as well as Zika virus and Dengue virus, all belong to the family Flaviviridae (23). Given that TNF-α plays a critical role in controlling viral infections by recruiting and activating macrophages, natural killer (NK) cells, T cells, and antigen-presenting cells, blockade of TNF-α may increase the risk of infection with viruses of the family Flaviviridae (24, 25). Currently, literature reports on viral infections induced by TNF-α inhibitors primarily involve herpesviruses, cytomegalovirus, and Epstein-Barr virus, with no reported cases of Flaviviridae family infections to date (15–18). Therefore, the adverse events related to the Flaviviridae family identified in this study represent a previously unreported and noteworthy pharmacovigilance finding. These signals warrant prospective clinical investigations to determine whether a causal relationship exists.

image4.png

Fig. 4. The top 30 ADEs of golimumab ranked by ROR at the PTs level. SOC – system organ class, PTs – preferred terms, ROR – reporting odds ratio, 95 % CI – 95 % confidence interval, PRR – proportional reporting ratio, EBGM05 – the lower limit of the 95 % CI for EBGM, IC025 – the lower limit of the 95 % CI for IC, ADEs – adverse drug events.

Fig. 5 presents 29 unlabeled malignancy signals, comprising 562 total reports. The six most frequently reported were basal cell carcinoma (n = 168), colon cancer (n = 81), cervical carcinoma (n = 33), invasive ductal breast carcinoma (n = 31), lung adenocarcinoma (n = 28), and rectal cancer (n = 26). By signal strength (ROR), the top three were ovarian endometrioid carcinoma (ROR: 34.99 [10.79–113.49]), malignant neoplasm of thymus (ROR: 28.77 [8.94–92.58]), and gingival cancer (ROR: 14.71 [4.65–46.50]). These findings are broadly corroborated by existing clinical evidence. Multiple clinical trials have reported malignancies in patients treated with golimumab, including breast cancer, basal cell carcinoma, lung adenocarcinoma, uterine cancer, Hodgkin's disease, and non-small cell lung cancer (16, 22, 26). In a study of golimumab for active RA, 8 patients (1.3 %) developed neoplasms, among whom 2 patients receiving 4 mg kg⁻¹ developed colorectal cancer and gastric cancer, resp. (27). Collectively, these data suggest that golimumab may be associated with a broad spectrum of malignancies, and routine oncological surveillance is therefore advisable in clinical practice. For patients with a history of malignancy, a comprehensive individualized assessment that weighs the tumor type, disease stage, metastatic risk, and patient preferences against the anticipated therapeutic benefits should be performed prior to initiating golimumab therapy.

image5.png

Fig. 5. Positive malignancy events not documented in the prescribing information. PTs – preferred terms, ROR – reporting odds ratio, 95 % CI – 95 % confidence interval, PRR – proportional reporting ratio, EBGM05 – the lower limit of the 95 % CI for EBGM, IC025 – the lower limit of the 95 % CI for IC.

Time-to-onset analysis and Weibull distribution testing

TTO analysis was performed on valid adverse event reports, with results presented in Fig. 6a. ADEs occurred most frequently within the first month of treatment (n = 1,353, 14.5 %), although the frequency of reports declined progressively within the first 180 days following treatment initiation; 4,262 reports (45.7 %) still occurred beyond one-year post-initiation, suggesting that ADEs may exhibit a delayed onset. As shown in Table III, Weibull distribution fitting revealed a median onset time of 292 days (IQR: 71–835.8 days). The upper bound of the 95 % confidence interval for the shape parameter (β) was less than 1.0 (β = 0.72), consistent with an early-failure pattern (β < 1), indicating that the risk of ADEs is highest early in the treatment course and diminishes over time. The Kaplan-Meier curve in Fig. 6b shows the cumulative incidence of ADEs related to golimumab. The curve demonstrates that the cumulative incidence of adverse events increased rapidly within the first 400 days after treatment initiation (reaching 92 % by day 400), after which it gradually plateaued. Collectively, although the risk is predominantly concentrated in the early phase, the occurrence of delayed-onset events nevertheless underscores the need for continued long-term follow-up in clinical practice (28).

image6.jpeg

Fig. 6. Time to onset (TTO) analysis (counted in days) of golimumab-related ADEs: a) distribution of time to onset of golimumab-related ADEs; b) cumulative incidence of golimumab-related ADEs over time.

Table III. Time to onset (TTO) of golimumab-associated adverse drug events and Weibull distribution analysis

TTO (days)Weibull distribution
Case reports (n)Median (d) (IQR)Min-maxScale parameter:α (95 % CI)Shape parameter:β (95 % CI)Type
9,326292.0(71.0–835.8)1.0–5684.0542.58 (533.88–551.28)0.72 (0.71–0.73)Early failure

TTO – time to onset, IQR – interquartile range, α – scale parameter of Weibull distribution, β – shape parameter of Weibull distribution (β < 1: early‑failure pattern; β = 1: constant risk; β > 1: wear‑out pattern), CI – confidence interval

Limitations of the study

The results of this study, which are based on analyses of the FAERS database, still have certain limitations. First, as a spontaneous pharmacovigilance reporting system, the FAERS database is inherently subject to underreporting, reporting bias, and varying data quality and completeness, all of which may affect the accuracy and reliability of signal detection. Second, disproportionality analyses using ROR, PRR, BCPNN and MGPS can only indicate statistical associations between drug exposure and adverse events; the signal strength merely quantifies the strength of association rather than absolute risk, and thus cannot be used to infer causality. Confirmation of causal relationships requires prospective clinical investigation and pharmaco-epidemiological studies. Third, FAERS data are predominantly contributed by reporters from North America and Europe, with a relatively low proportion of reports from Asian populations; this geographic imbalance may limit the generalizability of the findings. Fourth, the FAERS database has inherent limitations in controlling for important confounding factors, including concomitant medication use, disease severity, comorbidities, and a history of malignancy. This may introduce confounding bias into the signal associations identified in this study. To confirm these pharmacovigilance signals, future prospective cohort studies and pharmaco-epidemiological analyses with appropriate covariate adjustment are warranted (29).

CONCLUSIONS

This study provides a comprehensive real-world pharmacovigilance characterisation of golimumab using four simultaneous disproportionality algorithms applied to a 16-year FAERS dataset. The concordant signal approach identified 397 positive ADE signals, including 29 previously unlabeled malignancy signals and novel Flaviviridae-related infection signals, collectively extending the known safety profile of golimumab beyond what is currently described in its prescribing information. Time-to-onset modelling further established a temporal pattern consistent with early failure, offering quantitative benchmarks for post-initiation monitoring. These findings provide a real-world evidence basis to support regulatory review of the existing black-box warning and for evidence - based updates to prescribing guidance. Future prospective pharmaco-epidemiological studies with appropriate covariate adjustment will be essential to establish causal relationships and to further refine the long-term safety profile of golimumab across diverse patient populations.

Acronyms, abbreviations, symbols. – ADEs – adverse drug events, AS – ankylosing spondylitis, BCPNN – Bayesian confidence propagation neural network, C – concomitant, CI – confidence interval, DEMO – patient demographic and administrative information, DRUG – drug-specific details, EBGM – empirical Bayesian geometric mean, EBGM05 – the lower limit of the 95 % CI of EBGM, EMA – European Medicines Agency, FAERS – FDA Adverse Event Reporting System, FDA – Food and Drug Administration, I – interacting, IC – information component, IC025 – the lower limit of the 95 % CI of the IC, INDI – indications for drug use, IQR – interquartile range, MedDRA – medical dictionary for regulatory activities, MGPS – multi-item gamma Poisson shrinker, MySQL – MySQL database management system, NK – natural killer, OUTC – patient outcomes dataset in FAERS, PRR – proportional reporting ratio, PS – primary suspect, PT – preferred term, RA – rheumatoid arthritis, REAC – adverse reaction records, ROR – reporting odds ratio, RPSR – report source, SAEs – serious adverse events, SD – standard deviation, SOC – system organ class, SS – secondary suspect, THER – therapy initiation and cessation dates, TNF-α – tumor necrosis factor-alpha, TTO – time to onset, α – scale parameter of Weibull distribution, β – shape parameter of Weibull distribution, χ2– chi-squared test statistic.

Data availability. – The data supporting the findings of this study are publicly available from the FDA Adverse Event Reporting System (FAERS) database. FAERS data can be accessed at: https://fis.fda.gov/extensions/FPD-QDE-FAERS/FPD-QDE-FAERS.html. The aggregated dataset generated during the current study is available from the corresponding author upon reasonable request.

Acknowledgments. – We appreciate the FDA Adverse Event Reporting System (FAERS) for providing open access.

Funding. – The author(s) received no specific funding for this work.

Competing interests. – The authors declare no competing interests.

Author's contribution. – Conceptualisation, W.M. and X.L.; data curation, W.M., N.J., and H.L.; formal analysis, W.M., N.J., and Z.L.; investigation, W.M., N.J., and H.L.; methodology, W.M., Z.L., and X.L.; supervision, X.L.; writing, original draft preparation, W.M.; writing, review and editing, X.L. All authors have read and agreed to the published version of the manuscript.

REFERENCES

1. Q-Q. Ren, P. Chen, L. Qiao, Y.-K. Du, R.-S. Wang, X.-Q. Huang, H. Guo and J. Yan, TNF-α inhibitor etanercept improves cerebrovascular function in elderly RA patients: Findings from a randomized controlled trial, BMC Rheumatol. 9(1) (2025) Article ID 135 (11 pages);https://doi.org/10.1186/s41927-025-00589-6

2. A. T. Melo, R. Campanilho-Marques and J. E. Fonseca, Golimumab (anti-TNF monoclonal antibody): Where we stand today, Hum. Vaccin. Immunother. 17(6) (2021) 1586–1598;https://doi.org/10.1080/21645515.2020.1836919

3. M. Rossini, S. De Vita, C. Ferri, M. Govoni, G. Paolazzi, C. Salvarani and S. Adami, Golimumab: A novel anti-tumor necrosis factor, Biol. Ther. 3(2) (2013) 83–107; https://doi.org/10.1007/s13554-013-0012-y

4. F. S. Macaluso, S. Renna, A. Orlando and M. Cottone, The biologics of ulcerative colitis, Expert Opin. Biol. Ther. 17(2) (2017) 175–184; https://doi.org/10.1080/14712598.2017.1271871

5. M. E. Husni, A. Deodhar, S. Schwartzman, S. D. Chakravarty, E. C. Hsia, J. H. Leu, Y. Zhou, K. H. Lo and A. Kavanaugh, Pooled safety results across phase 3 randomized trials of intravenous golimumab in rheumatoid arthritis, psoriatic arthritis, and ankylosing spondylitis, Arthritis Res. Ther. 24(1) (2022) Article ID 73 (13 pages);https://doi.org/10.1186/s13075-022-02753-6

6. S. Yokoyama, Y. Ishii and J. Masuda, Persistence and safety of golimumab in elderly patients with rheumatoid arthritis and renal dysfunction in a real-world setting, Drugs Real World Outcomes 10(1) (2023) 51–60;https://doi.org/10.1007/s40801-022-00338-y

7. A. K. Ersbøll, Z. Huang, D. D. Hill, S. M. Hede, V. Andersen, K. Bolin, M. S. Kristensen, S. Esslinger, F. R. Hansen, E. Hertervig, L. Kallio, T. M. Kjærulff, S. Kloster, A. Krumme, J. D. Lewis, L. Mehkri, N. Qvist, L. C. Thygesen, C. Weinstein and A. Green, A longitudinal post-authorization safety study of golimumab in treatment of ulcerative colitis: A cohort study in Denmark and Sweden, 2013–2021, Drug Saf. 48(5) (2025) 541–558;https://doi.org/10.1007/s40264-025-01519-8

8. J. Yu, S. J. Park, H. W. Kim, Y. J. Lim, J. Park, J. M. Cha, B. D. Ye, T. O. Kim, H. S. Kim, H. S. Lee, S. Y. Jung, Y. Kim and C. H. Choi, Effectiveness and safety of golimumab in patients with ulcerative colitis: A multicenter, prospective, postmarketing surveillance study, Gut Liver 16(5) (2022) 764–774; https://doi.org/10.5009/gnl210335

9. Y. Ou, Z. Cui, S. Lou, C. Zhu, J. Chen, L. Zhou, R. Zhao, L. Wang and F. Zou, Analysis of tirzepatide in the US FDA adverse event reporting system (FAERS): A focus on overall patient population and sex-specific subgroups, Front. Pharmacol. 15 (2024) Article ID 1463657 (15 pages); https://doi.org/10.3389/fphar.2024.1463657

10. X. Le, Y. Zhang and J. Ma, Comprehensive analysis of adverse events associated with T-cell engagers using the FAERS database, Expert Opin. Drug Saf. 24(7) (2025) 821–830; https://doi.org/10.1080/14740338.2025.2470875

11. J. Wu, H. Wu, L. Chen, H. Liang, G. Huang, S. Yang, B. Chen, Y. Noguchi and Y. Shen, Safety of daratumumab in the real-world: A pharmacovigilance study based on FAERS database, Expert Opin. Drug Saf. 23(7) (2024) 905–916;https://doi.org/10.1080/14740338.2023.2296966

12. F. Zou, Z. Cui, S. Lou, Y. Ou, C. Zhu, C. Shu, J. Chen, R. Zhao, Z. Wu, L. Wang, Z. Chen, H. Chen and Y. Lan, Adverse drug events associated with linezolid administration: A real-world pharmacovigilance study from 2004 to 2023 using the FAERS database, Front. Pharmacol. 15 (2024) Article ID 1338902; https://doi.org/10.3389/fphar.2024.1338902

13. K. Krüger, G. R. Burmester, S. Wassenberg, M. Bohl-Bühler and M. H. Thomas, Effectiveness and safety of golimumab in patients with rheumatoid arthritis, psoriatic arthritis and ankylosing spondylitis under real-life clinical conditions: Non-interventional GO-NICE study in Germany, BMJ Open 8(6) (2018) Article ID e021082; https://doi.org/10.1136/bmjopen-2017-021082

14. H. Kim, Y. Kim and Y. Lee, Real-world safety and effectiveness of golimumab in rheumatic diseases: Post-marketing surveillance in Korea, Rheumatol. Ther. 8(3) (2021) 1393–1404; https://doi.org/10.1007/s40744-021-00352-6

15. M. E. Husni, A. Kavanaugh, F. Murphy, D. Rekalov, D. D. Harrison, L. Kim, K. H. Lo, J. H. Leu and E. C. Hsia, Efficacy and safety of intravenous golimumab through one year in patients with active psoriatic arthritis, Arthritis Care Res. (Hoboken) 72(6) (2020) 806–813; https://doi.org/10.1002/acr.23905

16. P. Emery, R. M. Fleischmann, I. Strusberg, P. Durez, P. Nash, E. J. B. Amante, M. Churchill, W. Park, B. Pons-Estel, C. Han, T. A. Gathany, S. Xu, Y. Zhou, J. H. Leu and E. C. Hsia, Efficacy and safety of subcutaneous golimumab in methotrexate-naive patients with rheumatoid arthritis: Five-year results of a randomized clinical trial, Arthritis Care Res. (Hoboken) 68(6) (2016) 744–752;https://doi.org/10.1002/acr.22759

17. T. Nemoto, S. Ito, D. Kobayashi, C. Takai, S. Sakai, Y. Kurosawa, E. Hasegawa, R. Okabayashi, A. Abe, H. Otani, H. Lee, A. Murasawa, I. Narita, K. Nakazono, Y. Toyoshima, K. Inagaki and H. Ishikawa, Long-term use of golimumab in daily practice for patients with rheumatoid arthritis, Intern. Med. 60(9) (2021) 1359–1367;https://doi.org/10.2169/internalmedicine.5381-20

18. C. Moriano Morales, J. Graña Gil, N. Brito García, J. L. Martín Varillas, V. Calvo Del Río, P. Moya Alvarado, F. J. Narváez García, G. Espinosa, P. Díaz Del Campo Fontecha, M. Guerra Rodríguez, J. Mateo Arranz, M. López Gómez, F. M. Francisco Hernández, M. M. Trujillo, R. Dos Santos Sobrín, J. I. Martín Sánchez, J. Maese Manzano and J. Suárez Cuba, SER recommendations on treatment of refractory Behçet's syndrome, Reumatol. Clin. (Engl. Ed.) 20(4) (2024) 204–217;https://doi.org/10.1016/j.reumae.2023.12.006

19. R. Arendse, P. Rahman, P. Baer, D. Haaland, L. Bessette, D. Sholter, M. Rachich, E. Rampakakis, A. M. Marrache, A. J. Lehman and O. Asin-Milan, Safety of golimumab in patients with rheumatoid arthritis, psoriatic arthritis, and axial spondyloarthritis: Results from a real-world Canadian setting, Rheumatol. Ther. 12(6) (2025) 1043–1055;https://doi.org/10.1007/s40744-025-00788-0

20. X. Sendaydiego, L. S. Gold, M. Dubreuil, J. S. Andrews, P. Reid, D. F. L. Liew, R. Goulabchand, G. C. Hughes, J. A. Sparks, J. G. Jarvik, S. Singh, J. W. Liew and N. Singh, Comparative safety of biologic and targeted synthetic disease-modifying anti-rheumatic drugs for cardiovascular outcomes in rheumatoid arthritis, Rheumatology (Oxford) 64(6) (2025) 3434–3443;https://doi.org/10.1093/rheumatology/keaf096

21. P. Le Blay, G. Mouterde, T. Barnetche, J. Morel and B. Combe, Short-term risk of total malignancy and nonmelanoma skin cancers with certolizumab and golimumab in patients with rheumatoid arthritis: Meta-analysis of randomized controlled trials, J. Rheumatol. 39(4) (2012) 712–715;https://doi.org/10.3899/jrheum.110982

22. J. S. Smolen, J. Kay, M. K. Doyle, R. Landewé, E. L. Matteson, J. Wollenhaupt, N. Gaylis, F. T. Murphy, J. S. Neal, Y. Zhou, S. Visvanathan, E. C. Hsia and M. U. Rahman, Golimumab in patients with active rheumatoid arthritis after treatment with tumour necrosis factor alpha inhibitors (GO-AFTER study): A multicentre, randomised, double-blind, placebo-controlled, phase III trial, Lancet 374(9685) (2009) 210–221;https://doi.org/10.1016/S0140-6736(09)60506-7

23. N. K. Björkström, B. Strunz and H. G. Ljunggren, Natural killer cells in antiviral immunity, Nat. Rev. Immunol. 22(2) (2022) 112–123; https://doi.org/10.1038/s41577-021-00558-3

24. V. V. Costa, C. T. Fagundes, D. G. Souza and M. M. Teixeira, Inflammatory and innate immune responses in dengue infection: Protection versus disease induction, Am. J. Pathol. 182(6) (2013) 1950–1961; https://doi.org/10.1016/j.ajpath.2013.02.027

25. P. L. Kung, T. W. Chou, M. Lindman, N. P. Chang, I. Estevez, B. D. Buckley, C. Atkins and B. P. Daniels, Zika virus-induced TNF-α signaling dysregulates expression of neurologic genes associated with psychiatric disorders, J. Neuroinflammation 19(1) (2022) Article ID 100; https://doi.org/10.1186/s12974-022-02460-8

26. J. Kay, R. Fleischmann, E. Keystone, E. C. Hsia, B. Hsu, M. Mack, N. Goldstein, J. Braun and A. Kavanaugh, Golimumab 3-year safety update: an analysis of pooled data from the long-term extensions of randomised, double-blind, placebo-controlled trials conducted in patients with rheumatoid arthritis, psoriatic arthritis or ankylosing spondylitis, Ann. Rheum. Dis. 74(3) (2015) 538–546; https://doi.org/10.1136/annrheumdis-2013-204195

27. J. Kremer, C. Ritchlin, A. Mendelsohn, D. Baker, L. Kim, Z. Xu, J. Han and P. Taylor, Golimumab, a new human anti-tumor necrosis factor alpha antibody, administered intravenously in patients with active rheumatoid arthritis: Forty-eight-week efficacy and safety results of a phase III randomized, double-blind, placebo-controlled study, Arthritis Rheum. 62(4) (2010) 917–928;https://doi.org/10.1002/art.27348

28. Q. Liu, Z. Cui, C. Deng, C. Yang and T. Shi, A real-world pharmacovigilance analysis of adverse events associated with irbesartan using the FAERS and JADER databases, Front. Pharmacol. 15 (2024) Article ID 1485190;https://doi.org/10.3389/fphar.2024.1485190

29. Y. C. Zhao, X. Li, C. Q. Wang, Y. Jiao, Y. N. Shen, T. J. Wang and C. H. Zhang, Unveiling the hidden risks: An update decade-long analysis of Abraxane-related adverse events from the FAERS database, Int. J. Nanomedicine 19 (2024) 11847–11858;https://doi.org/10.2147/IJN.S490400

Notes


This display is generated from NISO JATS XML with jats-html.xsl. The XSLT engine is libxslt.