ORIGINAL RESEARCH article

Front. Drug Saf. Regul., 10 June 2026

Sec. Clinical Drug Safety

Volume 6 - 2026 | https://doi.org/10.3389/fdsfr.2026.1827098

Comprehensive pharmacovigilance assessment of pairwise drug–drug interactions and acute pancreatitis in the FDA adverse event reporting system: focus on incretin-based drugs

  • IK

    Issei Kazume 1,2

  • KO

    Kenji Onda 1*

  • ST

    Sachiko Tanaka 1

  • TY

    Takashi Yokokawa 1

  • KS

    Kenichi Suzuki 1

  • 1. Department of Clinical Pharmacology, School of Pharmacy, Tokyo University of Pharmacy and Life Sciences, Tokyo, Japan

  • 2. P-One Clinic, Keikokai Medical Corp, Tokyo, Japan

Abstract

Background:

Drug-induced pancreatitis remains a clinically important condition owing to its severity and relevance to pharmacotherapeutic risk management. Incretin-based antidiabetic agents, including glucagon-like peptide-1 receptor agonists (GLP-1 RAs) and dipeptidyl peptidase-4 (DPP-4) inhibitors, are associated with an increased risk of acute pancreatitis (AP). However, the extent to which concomitant medication modifies AP has not yet been systematically elucidated. We conducted a comprehensive systematic drug–drug interaction (DDI) assessment involving GLP-1 RA and DPP-4 inhibitors using the U.S. Food and Drug Administration’s Adverse Event Reporting System (FAERS).

Methods:

Disproportionality analyses were conducted using the FAERS data (JAPIC FAERS) from the fourth quarter of 1997 to the first quarter of 2024. Crude reporting odds ratios (cRORs) were calculated to assess the association between individual drugs and AP. DDI signals for 1:1 concomitant use involving GLP-1 RAs or DPP-4 inhibitors were comprehensively evaluated using the Ω shrinkage measure, with Ω025 as the signal detection criterion. Ω025 were analyzed for all pooled cases, and temporal trends were assessed across six multi-year periods using a period-stratified analysis.

Results:

Elevated cRORs for AP were observed for several agents, including GLP-1 RAs, DPP-4 inhibitors, and other antidiabetic medications. Pairwise DDI analyses of GLP-1 RAs or DPP-4 inhibitors with other drugs identified positive interaction signals in 308 drug pairs in pooled analyses. Period-stratified analyses narrowed the findings to 11 drug pairs with consistent signals across ≥3 of 6 periods, including dulaglutide–atorvastatin (four periods). However, no concomitant drug combinations involving GLP-1 RAs or DPP-4 inhibitors were identified as being robustly associated with increased AP reporting.

Conclusion:

Period-stratified analysis helped mitigate false-positive signals identified in the pooled DDI analyses. Drug pairs exhibiting recurrent weak AP signals warrant continued pharmacovigilance and further evaluation.

1 Introduction

Acute pancreatitis (AP) is a leading gastrointestinal disorder that requires hospitalization and is a potentially fatal condition (; ). The global annual incidence is estimated at 34 cases per 100,000, and the associated healthcare costs are reported to reach approximately USD 9.3 billion per year. Notably, the incidence of AP continues to increase worldwide (; ; ), with its etiologies including gallstones, alcohol consumption, hypertriglyceridemia, drug-induced causes, and infections. Among these, drug-induced pancreatitis accounts for approximately 2.8%–5.3% of all cases; although relatively uncommon, it necessitates careful consideration in pharmacotherapy. To date, several classes of medications, such as antidiabetic agents, anticancer drugs, and antiviral agents, have been reported to have established or suspected causal associations with AP, based on adverse event and case reports (; ).

Glucagon-like peptide-1 receptor agonists (GLP-1 RAs) and dipeptidyl peptidase-4 (DPP-4) inhibitors are antidiabetic agents that target incretin hormones, and currently play a central role in the treatment of type 2 diabetes mellitus (; ). In addition to their glucose-lowering effects, GLP-1 RAs promote weight loss and have been increasingly used for the treatment of obesity in recent years, with further expansion of their anticipated use (). However, these agents also include warnings regarding the development of AP in their package insert, and their association with AP remains a subject of ongoing debate ().

Although previous studies have examined the relationship between individual agents and the occurrence of pancreatitis (; ; ; ), information regarding the extent to which the risk of AP is influenced when these drugs are used concomitantly with other medications remains limited. When multiple drugs are administered concurrently, drug–drug interactions (DDI) may increase the risk of AP; therefore, elucidating such interactions and appropriately assessing and mitigating the associated risks are essential for ensuring the safety of pharmacotherapy. Nevertheless, existing DDI research on drug-induced pancreatitis is scarce, and evaluations of the pancreatitis risk under polypharmacy conditions remain insufficient.

The U.S. Food and Drug Administration’s Adverse Event Reporting System (FAERS) is a spontaneous reporting database that collects information on adverse events associated with the use of a wide range of therapeutic agents (). Large-scale databases such as FAERS are particularly valuable for analyzing rare adverse events. Moreover, FAERS is an important resource for detecting potential DDIs and identifying previously unrecognized adverse events (); its application has been proposed for evaluating DDIs and identifying medications that may mitigate adverse drug reactions ().

Disproportionality analysis is commonly used to detect signals of drug–adverse event associations. Given the inherent limitations of spontaneous reporting data, such as the inability to estimate the true incidence rates of adverse drug reactions and the presence of reporting biases, direct comparisons of reporting odds ratios (RORs) derived from disproportionality analyses are considered inappropriate. However, multivariate analyses that incorporate logistic regression have been proposed to enable comparisons using adjusted RORs under specific conditions (; ; ). Furthermore, several methods have been developed to detect DDI signals between two agents, including the Ω shrinkage measure model, additive and multiplicative models, and the combination risk ratio model (). Although no standardized approach has been established, the Ω shrinkage measure model has been reported as the most conservative among these methods ().

In the current study, we conducted a 1:1 pairwise DDI analysis of AP reports involving GLP-1 RAs or DPP-4 inhibitors in combination with other medications, using FAERS. Through this approach, we aimed to identify the relevant DDIs associated with AP reporting.

2 Materials and methods

2.1 Data source and mining

This analysis was conducted using data from the Japan Pharmaceutical Information Center FAERS (JAPIC FAERS), covering the fourth quarter of 1997 to the first quarter of 2024. JAPIC FAERS is a database preprocessed by the Japan Pharmaceutical Information Center based on the original US FAERS data and has undergone data cleaning to ensure quality, including the removal of duplicate reports pertaining to the same patient. Because the FAERS is a publicly accessible and anonymized database, institutional review board approval was deemed unnecessary. This study followed the READUS-PV (REporting of A Disproportionality analysis for druUg Safety signal detection using individual case safety reports in PharmacoVigilance) guidelines (). The FAERS database consists of seven distinct tables, of which the drug (DRUG) and adverse event (REAC) tables were used. These tables were connected through PrimaryID and analyzed using a relational database software (Microsoft Access for Microsoft 365). We confirmed the validity of our aggregation system by demonstrating consistency with results reported in previous studies (; ; ).

In the DRUG table, each medication is classified according to its involvement in adverse events as a Primary Suspected (PS) drug, Secondary Suspected (SS) drug, concomitant (C) drug, or an interacting (I) drug. In the present study, all drug classifications were included in both the disproportionality analysis for individual drugs (for the calculation of ROR) and the DDI analysis. Medications registered under the same PrimaryID were defined as concomitant drugs and analyzed accordingly. In the DDI analysis, we regarded secondary suspected and interacting drugs as equivalent to concomitant drugs and did not apply a hierarchical classification based on drug role.

2.2 Definition of adverse events

All event reports were identified using the Preferred Terms (PTs) outlined in the Medical Dictionary for Regulatory Activities (MedDRA; version 27.0). AP was defined based on the PTs included in the Standard MedDRA Queries 20000022 with a narrow scope. A complete list of the PTs is provided in the Supplementary Table S1.

2.3 Disproportionality analysis of acute pancreatitis reports

The overall analysis was conducted in accordance with the flowchart shown in Figure 1. First, a disproportionality analysis of the drugs associated with AP was performed. A 2 × 2 contingency table was constructed based on the number of reports on drug exposure and AP for all drugs registered in FAERS (8,585 distinct agents). Crude reporting odds ratios, 95% confidence intervals, and chi-squared values were subsequently calculated (Figure 2).

FIGURE 1

FIGURE 2

2.4 Pairwise drug–drug interaction analysis of acute pancreatitis reporting involving GLP-1 receptor agonists or DPP-4 inhibitors

A DDI analysis was performed using the Ω shrinkage measure, a validated method for interaction signal detection in spontaneous reporting databases (; ). We comprehensively evaluated the interaction signals of one-to-one combinations of glucagon-like peptide-1 RAs (seven agents: albiglutide, dulaglutide, exenatide, liraglutide, lixisenatide, semaglutide, and tirzepatide) and dipeptidyl peptidase-4 inhibitors (alogliptin, linagliptin, saxagliptin, sitagliptin, and vildagliptin) with all the other reported medications (8,584 agents). The interaction signal metric (Ω025) for each drug pair was calculated and summarized, as illustrated in Figure 3. To assess the temporal consistency of the detected interaction signals, the study period (26.5 years) was divided into six consecutive intervals of several years. Ω025 values were calculated separately for each interval, and temporal trends in interaction signal strength were evaluated. For the period-stratified analysis, we used six time periods derived from the predefined twelve periods in the JAPIC FAERS dataset by combining adjacent intervals to ensure sufficient sample size and stability of the estimates. The periods were defined as follows: Period 1, 1997 Q4–2013 Q4; Period 2, 2014 Q1–2017 Q2; Period 3, 2017 Q3–2019 Q2; Period 4, 2019 Q3–2021 Q1; Period 5, 2021 Q2–2022 Q3; and Period 6, 2022 Q4–2024 Q1.

FIGURE 3

3 Results

3.1 Disproportionality analysis of AP reports

The cleaned FAERS dataset comprised 15,957,559 reports, of which 61,673 were identified as AP cases. Among the evaluated drugs, 1,247 were associated with significantly elevated crude RORs (cRORs). Table 1 presents the top 30 drugs ranked by χ2 values, as identified through the disproportionality analysis for AP. The highest reporting signal was observed for liraglutide, with a cROR (95% CI) of 19.86 (19.14–20.60). Elevated cRORs were also identified for sitagliptin at 10.90 (10.51–11.30), metformin at 4.50 (4.39–4.62), exenatide at 8.99 (8.61–9.39), and eluxadoline at 40.23 (36.00–44.96). When stratified by pharmacological class, a substantial number of drugs with increased cRORs were associated with glucose-lowering therapy, including GLP-1 RAs (liraglutide, exenatide, dulaglutide, and semaglutide), DPP-4 inhibitors (sitagliptin, linagliptin, and saxagliptin), metformin, glimepiride, and insulin. In addition to antidiabetic agents, elevated cRORs were observed for eluxadoline, which is indicated for irritable bowel syndrome; fenofibrate, which is prescribed for hypertriglyceridemia; and several agents used in oncology, human immunodeficiency virus treatment, and antipsychotic therapy.

TABLE 1

Drugn11n12n21n22cRORΧ2
Liraglutide3,23544,19358,43815,851,69319.86 (19.14–20.6)51155.8
Sitagliptin3,13277,66258,54115,818,22410.9 (10.51–11.3)25692
Metformin7,015440,42254,65815,455,4644.5 (4.39–4.62)16687.2
Exenatide2,20865,35459,46515,830,5328.99 (8.61–9.39)14634.4
Eluxadoline3612,32661,31215,893,56040.23 (36–44.96)11885.6
Fenofibrate1,39149,08060,28215,846,8067.45 (7.06–7.86)7384.23
Dulaglutide1,70772,44759,96615,823,4396.22 (5.92–6.53)7100.16
Semaglutide1,10135,40660,57215,860,4808.14 (7.66–8.65)6570.98
Pegaspargase5269,14261,14715,886,74414.95 (13.69–16.33)6418.72
Didanosine4566,99361,21715,888,89316.92 (15.39–18.61)6367.12
Quetiapine2,363144,80659,31015,751,0804.33 (4.16–4.52)5734.71
Olanzapine1,58685,91160,08715,809,9754.86 (4.62–5.11)4648
Asparaginase4258,97361,24815,886,91312.29 (11.14–13.55)4177.87
Azathioprine1,36177,99760,31215,817,8894.58 (4.34–4.83)3656.39
Pancrelipase65323,34161,02015,872,5457.28 (6.73–7.87)3403.26
Linagliptin60522,53061,06815,873,3566.98 (6.44–7.57)2988.95
Stavudine3488,55761,32515,887,32910.54 (9.46–11.73)2869.91
Mesalazine1,14171,31560,53215,824,5714.18 (3.94–4.44)2669.51
Mercaptopurine46015,52261,21315,880,3647.69 (7–8.44)2580.06
Glimepiride99158,47360,68215,837,4134.42 (4.15–4.71)2540.38
Hydrochlorothiazide2,948312,15758,72515,583,7292.51 (2.41–2.6)2517.33
Pantoprazole2,966320,29258,70715,575,5942.46 (2.37–2.55)2416.93
Insulin3,743451,73857,93015,444,1482.21 (2.14–2.28)2307.55
Lisinopril2,756298,19758,91715,597,6892.45 (2.35–2.54)2231.95
Saxagliptin3018,29961,37215,887,5879.39 (8.37–10.54)2166.65
Ezetimibe1,08080,69660,59315,815,1903.49 (3.29–3.71)1863.33
Simvastatin2,400265,36259,27315,630,5242.39 (2.29–2.48)1838.71
Pancreatin3069,84561,36715,886,0418.05 (7.18–9.02)1822.17
Rufloxacin11761,66215,895,879405.1 (157.03–1045.04)1724.07
Tigecycline1663,34161,50715,892,54512.84 (10.98–15.01)1721.65

Disproportionality analysis of acute pancreatitis (AP) reports.

n11, number of reports including both the target drug and acute pancreatitis; n12, reports including the target drug without acute pancreatitis; n21, reports excluding the target drug but including acute pancreatitis; n22, reports excluding both the target drug and acute pancreatitis; cROR, crude reporting odds ratio; χ2, chi-square value calculated from the corresponding 2 × 2 contingency table.

3.2 Pairwise DDI analysis for AP reporting involving GLP-1 RAs or DPP-4 inhibitors

A comprehensive DDI analysis was conducted using data registered in the FAERS database to evaluate the impact of one-to-one drug combinations on AP reporting. Using the Ω shrinkage measure model described in the Materials and Methods section, interaction signal estimates (Ω025) were systematically calculated for pairwise combinations of GLP-1 RAs or DPP-4 inhibitors with all other reported medications. In the analysis encompassing the entire study period from 1997 to 2024, positive interaction signals based on Ω025 values were identified for 308 drug pairs, corresponding to 0.3% of all evaluated combinations (103,008 pairs). The most frequently reported drug pair was insulin in combination with exenatide, with 532 reports exhibiting a Ω025 value of 0.04. This was followed by acetylsalicylate in combination with exenatide, which accounted for 436 reports and exhibited an Ω025 value of 1.1. Other highly-ranked combinations included exenatide with lisinopril (n111 = 381, Ω025 = 1.15), acetylsalicylate with liraglutide (n111 = 347, Ω025 = 0.03), exenatide with simvastatin (n111 = 331, Ω025 = 1.42), exenatide with hydrochlorothiazide (n111 = 329, Ω025 = 1.28), and atorvastatin with exenatide (n111 = 300, Ω025 = 0.84) (Table 2).

TABLE 2

Drug ADrug Bn111n11+E111Ω (Ω025−Ω975)Judge
ExenatideInsulin53213,773476.40.16 (0.04–0.28)P
AcetylsalicylateExenatide4636,536196.991.23 (1.1–1.36)P
ExenatideLisinopril3814,661154.661.3 (1.15–1.44)P
AcetylsalicylateLiraglutide3474,437305.220.18 (0.03–0.34)P
ExenatideSimvastatin3313,342110.851.57 (1.42–1.73)P
ExenatideHydrochlorothiazide3293,575121.131.44 (1.28–1.59)P
AtorvastatinExenatide3004,469149.21.01 (0.84–1.17)P
LiraglutideLisinopril2732,488174.340.65 (0.47–0.82)P
HydrochlorothiazideLiraglutide2472,434172.470.52 (0.34–0.7)P
ExenatideParacetamol2362,15969.671.75 (1.57–1.94)P
LiraglutideSimvastatin2132,000141.180.59 (0.4–0.79)P
ExenatideFenofibrate1991,43871.81.46 (1.26–1.66)P
ExenatideMetoprolol1872,91995.950.96 (0.75–1.17)P
ExenatideFurosemide1812,60385.671.07 (0.86–1.28)P
LiraglutideSitagliptin1701,466141.370.27 (0.05–0.48)P
ExenatideLevothyroxine1693,527117.50.52 (0.3–0.74)P
ExenatideEzetimibe1681,63062.271.42 (1.21–1.64)P
GlimepirideLiraglutide1651,751134.080.3 (0.08–0.52)P
AmlodipineExenatide1632,64689.250.87 (0.64–1.09)P
ExenatideOmeprazole1492,34481.310.87 (0.64–1.1)P
ExenatideRosuvastatin1371,93567.91.01 (0.77–1.25)P
AtorvastatinDulaglutide1333,22682.970.68 (0.43–0.92)P
LiraglutidePioglitazone1211,20390.660.41 (0.16–0.67)P
ClopidogrelExenatide1171,52049.421.23 (0.97–1.5)P
ExenatideValsartan1161,55349.61.22 (0.95–1.48)P
EsomeprazoleExenatide1131,61758.180.95 (0.69–1.22)P
ExenatidePantoprazole1061,228451.23 (0.95–1.5)P
AcetylsalicylateDulaglutide1022,93671.980.5 (0.22–0.78)P
ExenatideRamipril1001,02437.681.4 (1.11–1.68)P
AtenololExenatide1001,14139.731.32 (1.04–1.6)P

Drug–drug interaction analysis of AP between incretin-based drugs and concomitant medications in all pooled periods.

n111, number of reports including Drug A, Drug B, and AP; n11+, reports including both drugs; E111, expected value of the Ω shrinkage measure model; P, positive signal. Drug pairs with n111 ≥ 100 were listed. Drug pairs (Drugs A and B) are presented in alphabetical order.

3.3 Temporal analysis of interaction signals

To assess the temporal consistency of the interaction signals, we examined changes in the number of positive signals of Ω025 across predefined time periods. No drug pair demonstrated a consistently robust interaction signal across six or five periods. Among the evaluated combinations, the atorvastatin–dulaglutide pair showed positive interaction signals across four periods, representing the highest degree of temporal consistency observed in the analysis. Positive signals were identified for several drug pairs across three periods, including acetylsalicylate–exenatide, liraglutide–lisinopril, exenatide–paracetamol, esomeprazole–exenatide, enalapril–exenatide, exenatide–ondansetron, exenatide–lorazepam, bisoprolol–exenatide, exenatide–morphine, and enoxaparin–exenatide (Figure 4). In addition, 48 drug pairs exhibited positive interaction signals across two periods.

FIGURE 4

4 Discussion

Using FAERS, one of the world’s largest spontaneous reporting databases, we conducted disproportionality analyses of individual drugs associated with AP, as well as DDI analyses examining the pairwise (1:1) concomitant use of GLP-1 RAs or DPP-4 inhibitors with other medications. Several considerations are critical for the appropriate interpretation of the present findings. First, because this study was based on a spontaneous reporting database, causal relationships between drug exposure and AP could not be inferred. Second, our results should not be interpreted as quantitative risk estimates. Rather, the identified associations require external validation, and the primary role of this analysis is to support the prioritization of hypotheses for subsequent investigations. Accordingly, this study is exploratory in nature, providing an initial signal detection framework to guide further pharmacological and epidemiological evaluation.

In the disproportionality analyses, cRORs were elevated for GLP-1 RAs, DPP-4 inhibitors, and oral antidiabetic agents, such as metformin, indicating strong associations with reports of AP. A wide range of additional drugs also yielded high signal values, which were largely concordant with previously reported evidence (), thereby supporting the validity of the analysis.

GLP-1 RAs and DPP-4 inhibitors are incretin-based therapies that were first introduced with the approval of exenatide in 2005 and sitagliptin in 2006. Since their initial market entry, multiple agents within these classes have been approved and now play a central role in the management of type 2 diabetes mellitus. Liraglutide, originally approved for the treatment of type 2 diabetes in 2009, was the first GLP-1 RA approved globally for treating obesity in 2014. Given the expanding indications and increasing use of GLP-1 RAs, their clinical relevance is expected to increase. Although warnings regarding AP are included in product labels for both GLP-1 RAs and DPP-4 inhibitors, the causal relationship between these agents and AP remains a subject of ongoing debate (; ; ; ). Importantly, patients treated with these agents often have underlying conditions such as diabetes, obesity, or metabolic syndrome, each of which independently confers an elevated risk of pancreatitis, necessitating careful clinical interpretation that accounts for the baseline risk (; ). In addition to antidiabetic medications, eluxadoline, a treatment for diarrhea-predominant irritable bowel syndrome, has also been implicated in AP development (). Consistent with previous reports, eluxadoline showed a markedly elevated signal in this study. Additionally, drugs previously associated with AP—including acute lymphoblastic leukemia treatments, such as pegaspargase and asparaginase (; ); antiretroviral agents, such as didanosine and stavudine (; ); and antipsychotics, such as quetiapine and olanzapine (; )—also showed elevated signals, further reflecting established clinical knowledge.

In the individual-drug disproportionality analyses, all drugs recorded in the reports were included regardless of their designated role (PS, SS, concomitant, or I). This inclusive approach was adopted because the primary objective of the present study was to evaluate DDI rather than quantify the contribution of individual suspected drugs alone. However, this strategy introduces potential indications and treatment bias, particularly through the inclusion of medications used for the management of pancreatitis. Accordingly, such drugs should be excluded from the interpretation to avoid distortion of signal estimates driven by reverse causality.

Regarding the pairwise DDI analysis with the Ω shrinkage measure model, in the pooled analysis covering the entire observation period, positive signals based on the Ω025 criterion were identified for 308 drug pairs. Drug combinations with the highest number of reports were insulin and exenatides. However, despite its high reporting frequency in the pooled dataset, this combination did not yield a positive shrinkage signal in any period-specific analysis. Notably, the pooled Ω025 value for this pair was 0.04, lying at the margin of the predefined signal-detection threshold. Insulin is widely used in the management of diabetes mellitus and may be administered to treat hyperglycemia secondary to AP. Therefore, the observed association may reflect confounding by indication or treatment-related bias rather than a true pharmacological interaction. The second most frequently reported combination was acetylsalicylate and exenatide, with 436 reports and an Ω025 value of 1.10. In this case, the number of reports for the observed co-occurrence (n111) exceeded the expected count (E111), resulting in a positive Ω025 signal.

Next, to assess temporal consistency, we examined the number of positive signals observed across individual reporting periods. Across the six periods, no drug combination exhibited positive signals in five or more periods, and no combination demonstrated interaction signals that were consistently observed across all periods. Eleven drug pairs showed positive signals in four or three periods, all of which involved GLP-1 RAs. Among these, the atorvastatin–dulaglutide combination showed the highest temporal recurrence, with positive signals observed consecutively from Periods 2 to 5. Drug pairs with positive signals in two periods totaled 48, of which 46 pairs involved GLP-1 RAs, and only two pairs involved DPP-4 inhibitors. Among these, most drug pairs involved exenatides.

Exenatide was the first GLP-1 RA to receive regulatory approval and was subsequently accompanied by an accumulation of post-marketing reports on AP, which initiated broader debate regarding the potential association between GLP-1 RAs and AP. Therefore, the predominance of positive signals observed during the earlier reporting periods may reflect this historical context, including heightened awareness following market introduction and the possibility of transient stimulated reporting. Accordingly, the potential influence of temporal reporting dynamics should be considered when interpreting these findings.

Among the drug pairs that exhibited positive signals in three or four reporting periods, combinations of statins and angiotensin-converting enzyme (ACE) inhibitors were identified. These drug classes have previously been reported to be associated with AP (; ; ; ; ; ; ), and the present disproportionality analysis likewise demonstrated significantly elevated cRORs for their individual use. This observation prompted further investigation of potential DDI signals between GLP-1 RAs and statins or ACE inhibitors, as detailed in the Supplementary Materials S2, S3. However, no consistent pattern of increased reporting was observed for these drug combinations, and no sustained tendency toward enhanced reporting emerged across the reporting periods. Among other combinations, the exenatide–ondansetron and exenatide–paracetamol (acetaminophen) pairs exhibited positive signals in three reporting periods. Although such associations are infrequent, prior reports have described links between paracetamol and AP (), as well as between ondansetron and AP (), which may have contributed, in part, to the observed signals, warranting further investigation.

Analyses based on the entire observation period, and those stratified by discrete time intervals, have distinct methodological strengths and limitations. Pooled analysis across the entire observation period offers the advantage of increased statistical power owing to larger sample sizes, thereby enhancing the sensitivity of signal detection. However, this approach does not allow for the assessment of the temporal consistency of detected signals. As a result, signals identified in the pooled analyses may reflect transient reporting phenomena rather than stable drug–event associations. By contrast, period-specific analyses enable the identification of signals that might be obscured when the data are averaged across the entire study period. This approach facilitates the evaluation of signal consistency over time and may help mitigate the impact of stimulated reporting associated with regulatory actions, media attention, or changes in clinical practice. A period-stratified approach may facilitate the prioritization of candidate signals by assessing temporal consistency across multiple signals emerging from pooled analyses. However, stratification by time inevitably reduces the number of reports within each period, leading to diminished statistical power and an increased risk of false-negative findings. Given these complementary strengths and weaknesses, pooled and period-specific analyses should not be viewed as competing approaches, but rather as mutually reinforcing strategies. Their combined application allows for a more nuanced and robust interpretation of pharmacovigilance signals, balancing sensitivity with temporal validity.

Several limitations inherent to studies based on spontaneous reporting databases should be acknowledged. First, issues of under-reporting and over-reporting are unavoidable, and causal relationships between drug exposure and adverse events cannot be established from such data. Second, the present analysis was restricted to pairwise DDIs, and the potential effects of concomitant use of three or more drugs were not evaluated. In addition, relevant clinical and demographic factors, including age, sex, underlying comorbidities (such as hepatic or renal dysfunction), dosage information, formulation, dosing frequency, route of administration, and lifestyle factors (e.g., alcohol consumption), were not accounted for in the analysis. Concomitant drug use was defined as medications listed within the same report. Consequently, detailed information on treatment duration, temporal sequencing, or the exact overlap of drug exposure was unavailable, and some reports may have included medications administered before or after treatment switches, rather than truly concurrent use. In real-world clinical practice, interpretation of drug safety signals requires careful consideration of patient-specific contextual factors. Regarding the temporal analysis, drug availability may vary across time periods. Recently approved drugs are not captured in earlier periods, whereas some older drugs may be reported less frequently due to the introduction of newer agents or changes in therapeutic guidelines. This may be reflected in Figure 4, where some missing values of period-specific Ω025 were observed for exenatide combinations, particularly in later periods. Additionally, dividing the data into a larger number of time periods may result in an insufficient number of cases to ensure adequate statistical power, whereas using broader time intervals may limit the ability to assess the robustness of the observed signals. Therefore, the optimal division of time periods remains uncertain and should be considered a limitation of this study. Given these limitations, the findings of this study should be interpreted as hypothesis-generating and do not necessarily reflect the actual clinical risk. Verification of these hypotheses requires further investigation, including pharmacological mechanistic studies and epidemiological validation using other large-scale healthcare databases. Despite these limitations, a comprehensive analysis of the world’s largest spontaneous reporting data source remains a valuable initial step toward risk assessment. Such large-scale approaches are particularly well-suited for detecting subtle safety signals that are difficult to identify at the level of individual case reports, thereby providing an important foundation for subsequent confirmatory research.

5 Conclusion

Our analyses demonstrated that period-stratified analysis effectively reduced false-positive signals identified in the pooled DDI analyses. No concomitant medications were found to clearly increase AP reporting when combined with GLP-1 RA or DPP-4 inhibitors. Nevertheless, specific drug pairs exhibited weak but recurrent signals across multiple reporting periods, including dulaglutide–atorvastatin, suggesting that these combinations warrant continued pharmacovigilance, cautious interpretation, and longitudinal monitoring during post-marketing safety surveillance.

Statements

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: https://fis.fda.gov/extensions/FPD-QDE-FAERS/FPD-QDE-FAERS.html.

Ethics statement

Ethical approval was not required for the studies involving humans because FAERS database is publicly available information. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements because FAERS database contains only anonymized, de-identified data.

Author contributions

IK: Data curation, Investigation, Writing – original draft. KO: Conceptualization, Data curation, Project administration, Validation, Writing – review and editing. ST: Writing – review and editing. TY: Writing – review and editing. KS: Writing – review and editing.

Funding

The author(s) declared that financial support was not received for this work and/or its publication.

Acknowledgments

We thank the Japan Pharmaceutical Information Center (JAPIC: https://www.japic.or.jp) for curating the FAERS data.

Conflict of interest

Author IK was employed by P-One Clinic, Keikokai Medical Corp.

The remaining author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was used in the creation of this manuscript. ChatGPT ver. 5.2 was used only to English proofreading.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fdsfr.2026.1827098/full#supplementary-material

References

Summary

Keywords

acute pancreatitis, DPP-4 inhibitors, drug–drug interaction, FAERS, GLP-1 receptor agonists

Citation

Kazume I, Onda K, Tanaka S, Yokokawa T and Suzuki K (2026) Comprehensive pharmacovigilance assessment of pairwise drug–drug interactions and acute pancreatitis in the FDA adverse event reporting system: focus on incretin-based drugs. Front. Drug Saf. Regul. 6:1827098. doi: 10.3389/fdsfr.2026.1827098

Received

10 March 2026

Revised

04 May 2026

Accepted

25 May 2026

Published

10 June 2026

Volume

6 - 2026

Edited by

Vincenzo Patella, “Santa Maria della Speranza” Hospital, Italy

Reviewed by

Wieslawa Agnieszka Fogel, Polish Academy of Sciences, Poland

Alison Yeomans, Drug Safety Research Unit, United Kingdom

Updates

Copyright

*Correspondence: Kenji Onda,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics