Abstract
Only a few studies investigated the association between proton pump inhibitor (PPI) use and pancreatic cancer, with inconsistent results. Moreover, these studies had a number of methodologic limitations. Our objective was to assess this association in a nationwide case–control study.
We used the French National Health Data System (SNDS), covering 99% of the French population since 2006. Incident cases of pancreatic cancer, identified between 2014 and 2018, were matched with up to four controls on year of birth, sex, frequency of hospitalization within 8 years prior to index date, and department of residence. Associations between PPIs and pancreatic cancer were estimated using conditional logistic regression models adjusted for sociodemographic characteristics, risk factors of pancreatic cancer (including diabetes mellitus, tobacco-related diseases, and morbid obesity), and other comorbidities.
A total of 23,321 cases of pancreatic cancer (mean age, 69.8 years; 51.7% males) and 75,937 matched controls were included. Overall, 77.8% of cases and 75.5% of controls were PPI ever users. Ever (vs. never) PPI use was associated with an increased risk of pancreatic cancer [adjusted OR (aOR) = 1.05, 95% confidence interval (CI), 1.01–1.09]. A dose–response relationship was observed [1–30 cumulative defined daily dose (cDDD): aOR = 0.92, 95% CI, 0.87–0.97; 31–180 cDDD: aOR = 1.05, 95% CI, 1.00–1.11; 181–1,080 cDDD: aOR = 1.18, 95% CI, 1.12–1.24; >1,080 cDDD: aOR = 1.17, 95% CI, 1.10–1.23].
On the basis of these findings, a slight increase in the risk of pancreatic cancer associated with high cumulative doses of PPIs cannot be excluded.
Given the overuse of PPIs, efforts should be continued to limit treatments to appropriate indications and durations.
Introduction
Since their market introduction in the late 1980s, proton pump inhibitors (PPI) have proven their efficacy and have become a standard treatment for acid-related conditions such as peptic ulcer disease and gastroesophageal reflux disease (1). Their use has steadily increased, and they are currently one of the most commonly prescribed classes of drug worldwide (2, 3). Misuse, such as coprescribing with nonsteroidal anti-inflammatory drugs (NSAID) without gastrointestinal risk, or prescribing without a clear indication (3), is widespread, reaching on average 50% among outpatients (4). Despite a good overall tolerance at short term, prolonged PPI use has raised safety concerns. Risks of long-term PPI therapy have been extensively explored in the literature, although some of them are still debated (5). Potential gastrointestinal health outcomes include infections (6), inflammatory bowel diseases (7), and malignancies (8–10). However, only a few observational studies have examined the association between PPI use and the risk of pancreatic cancer as main outcome, with inconsistent results and methodologic limitations arising from limited numbers of cases or long-term PPI users, concerns of reverse causality, or inability to capture important confounders (11–18). Plausible mechanisms have been suggested for the potential carcinogenic effect of PPIs in pancreatic cancer, related to induced hypergastrinemia (19) and microbiome alterations (20). While risk factors of pancreatic cancer are still insufficiently known, incidence rates are rising in developed countries, with over seven cases per 100,000 person-years in Northern America and Western Europe. After diagnosis, pancreatic cancer is associated with a poor prognosis, with an estimated 5-year survival rate of less than 5% (21–23). Thus, clarifying the impact of PPI exposure on the risk of pancreatic cancer is of major importance.
The aim of this study was to investigate the association between PPI use and the risk of pancreatic cancer in France, based on a large nationwide, population-based case–control study, addressing methodologic limitations of previous studies.
Materials and Methods
Data sources
This study was conducted using the French National Health Data System (Système National des Données de Santé, SNDS), consisting of comprehensive sociodemographic and medical individual information for 99% of the population living in France (about 67 million people) since 2006. The database contains data about all outpatient services reimbursed by the National Health Insurance, including drugs [coded according to the Anatomical Therapeutic Chemical Classification System (ATC; ref. 24)], and physician visits. Patients with costly chronic diseases (long-term diseases), such as cancer, are fully reimbursed for their health expenditures, and the diagnosis are recorded [coded according to the International Classification of Diseases, 10th Revision (ICD-10; ref. 25)]. The database also contains diagnoses related to hospital admissions, and procedures performed during hospital stays. A detailed presentation of the SNDS databases is available in the Supplementary Data S1.
Study population
Cases
We identified all patients ages 40 to 85 years, with an incident primary pancreatic cancer between January 1, 2014 and December 31, 2018. The index date was the date of first mention of pancreatic cancer (ICD-10 code C25), either as a hospital discharge diagnosis, or as a cause of long-term disease if it was further followed by a hospital discharge diagnosis within the next 3 months. We focused on pancreatic adenocarcinoma, which accounts for about 9 of 10 of all pancreatic cancers (21, 23). Thus, patients with a neuroendocrine neoplasm of the pancreas were not included. They were identified by an ICD-10 code C25.4, and/or an outpatient treatment with somatostatin analogs (ATC codes H01CB02, H01CB03) in the year following index date. To note, cases of pancreatic adenocarcinoma receiving a somatostatin analog in the perioperative setting, to reduce the risk of postoperative pancreatic leaks, were not excluded from the analyses.
Controls
Four controls with no diagnosis of pancreatic cancer at index date were randomly selected after matching on year of birth, sex, frequency of hospitalization within 8 years before index date (figured by the number of calendar years with at least one hospital admission, categorized as: 0, 1, 2, or ≥3), and department of residence. The index date of each case was assigned to the matched controls.
Exclusion criteria
Exclusion criteria for cases and controls were death on index date, absence of outpatient claim 7 or 8 years prior to index date (to assign the same length of observation to cases and to their matched controls, ensuring equal time windows to measure exposure and to identify comorbidities), and history of cancer (all causes) or pancreatic abnormality within 8 years before index date (ICD-10 codes available in the Supplementary Table S1).
Exposure to PPIs
Exposure to PPIs was defined as redeeming at least one prescription of a PPI marketed in France, namely omeprazole, esomeprazole, lansoprazole, pantoprazole, or rabeprazole (ATC codes A02BC01 to A02BC05), between January 1, 2006 and the index date. For ever users, we calculated the cumulative defined daily dose (cDDD), classified in quartile categories based on the distribution of use in cases. There is no consensus on optimal treatment duration or agreed definition of long‐term PPI use (26, 27). On the basis of information contained in the Summary of Product Characteristics, individuals were considered long-term users if they had been exposed to a cumulative dose of 181 DDD, equivalent to a 6-month therapy within the study period.
The SNDS does not contain information neither on inpatient nor on over-the-counter (OTC) PPI use. However, this use accounts for a limited proportion in France. In 2015, 92% of PPI boxes were delivered to outpatients and almost 97% of them were obtained from prescriptions (source: French National Agency for Medicines and Health Products Safety, ANSM).
Exclusion of incident PPI users
To allow for latency, and to minimize reverse causality, incident PPI users within 2 years before index date (defined by at least one redeemed prescription of PPI within 24 months prior to index date, and no redeemed prescription of PPI between 24 and 36 months prior to index date) were excluded from the set of cases and matched controls. Strata containing only cases or controls after exclusion of incident PPI users within 2 years before index date were removed.
Covariates
Sociodemographic characteristics were affiliation to complementary universal health insurance (CMUC, free access to healthcare for low-income people under 65 years old), and social deprivation index (levels of disadvantage calculated across small geographic areas). We also identified medical covariates, defined by a diagnosis, or, if appropriate, by at least three redeemed prescriptions, within 8 years prior to index date (codes available in the Supplementary Table S2): (i) potential risk factors of pancreatic cancer, defined according to the current best available evidence (21, 22, 28, 29): diabetes mellitus, tobacco-related diseases (including COPD diagnosis) or drug use, morbid obesity, alcohol-related diseases or drug use, acute pancreatitis, chronic pancreatitis, pancreatic cyst, gallstones, hepatitis B or C (ii); proxies of potential contexts of PPI treatment: gastroesophageal reflux disease, peptic ulcer, Helicobacter pylori eradication (iii); and other comorbidities or drug use: myocardial infarction, congestive heart failure, peripheral vascular disease, cerebrovascular disease, dementia, chronic obstructive pulmonary disease, connective tissue disease, mild liver disease, hemiplegia, moderate to severe liver disease, moderate to severe chronic kidney disease, HIV/AIDS, antihypertensive drug use, NSAID use, statin use.
Statistical analysis
Characteristics of the cases and matched controls were presented using descriptive statistics.
Associations between exposure to PPIs and pancreatic cancer were estimated based on crude and adjusted ORs (aOR) and their 95% confidence intervals (CI) obtained using conditional logistic regression models. In addition to the matching variables (year of birth, sex, history of hospitalizations, and department of residence), and calendar year, accounted for by design, potential risk factors of pancreatic cancer (diabetes mellitus, tobacco-related diseases or drug use, morbid obesity, alcohol-related diseases or drugs use, acute pancreatitis, chronic pancreatitis, pancreatic cyst, gallstones, hepatitis B or C) were forced in the adjusted models. Then, final models were run introducing the remaining covariates through forward selection process.
Complementary analyses stratified by age group, sex, calendar year, and cancer localization were conducted. We also assessed whether the main modifiable risk factors of pancreatic cancer (28, 29), which were reported with a prevalence above 1% among cases and controls, were effect modifier. To this end, we included in the model an interaction term between PPI use and the following covariates: diabetes mellitus, tobacco-related diseases or drug use, or morbid obesity.
The robustness of the main results was assessed in four sensitivity analyses. First, we applied a 2-year and a 4-year lag before the index date, disregarding PPI exposure in these periods, to maximize the control of reverse causality (30). Second, analyses were restricted to new PPI users, excluding patients who received a PPI in 2006. Third, we used inverse probability of treatment weighting (IPTW) to reduce confounding. Stabilized weights were computed from a logistic model using all cases and controls, adjusted for all covariates (except for CMUC, only available among individuals ages under 65 years). Weights were then introduced in a logistic regression model for the outcome, including no other predictor than exposure (31–33). Fourth, a sensitivity analysis considering histamine‐2‐receptor antagonists (H2RA, ATC codes A02BA01 to A02BA04) as an active comparator was conducted (34). These drugs have similar indications to PPIs, but inhibit acid secretion less profoundly than PPIs (35). Because of their well-documented superiority in relieving symptoms and healing mucosal lesions, PPIs have rapidly replaced H2RAs in treating any clinical acid-related condition (4). Most previous studies suggested a lack of association between pancreatic cancer development and H2RAs use (11, 13, 18). Thus, we compared the risk of pancreatic cancer between PPI and H2RA ever users, excluding those with incident PPI or H2RA use within 2 years before index date. Subjects using both PPIs and H2RAs were defined as PPI users.
All analyses were conducted using SAS EG (Copyright 2017 SAS Institute Inc.).
Data availability
The data generated in this study are not publicly available. EPIPHARE (https://www.epi-phare.fr/en/) has a regulatory permanent access to the data from the French National Health Data System (SNDS) via its constitutive bodies ANSM and CNAM. This permanent access is given according the French Decree No. 2016-1871 of December 26, 2016 relating to the processing of personal data called "National Health Data System" [Décret no. 2016–1871 du 26 décembre 2016 relatif au traitement de données à caractère personnel dénommé « système national des données de santé » (Internet). 2016 (cited 2021 Mar 12). Available from: https://www.legifrance.gouv.fr/affichTexte.do?cidTexte=JORFTEXT000033702840&categorieLien=id] and French law articles Art. R. 1461-13 [Article R1461-13 - Code de la santé publique - Légifrance (Internet). (cited 2021 Mar 12). Available from: https://www.legifrance.gouv.fr/codes/article_lc/LEGIARTI000038789574/] and 14 [Article R1461-14 - Code de la santé publique - Légifrance (Internet). (cited 2021 Mar 12). Available from: https://www.legifrance.gouv.fr/codes/article_lc/LEGIARTI000037678676/].
All data were deidentified, thus, informed consent was not necessary.
Results
We identified a total of 64,348 individuals ages 40–85 years, with incident pancreatic cancer between January 1, 2014 and December 31, 2018, including 58,599 (91.1%) patients with pancreatic adenocarcinoma. Among them, 37.1% (N = 19,686) were excluded because of history of cancer within 8 years prior to index date. Seven cases failed to be matched with controls. Almost 30% (N = 9,927) of cases and 18.7% (N = 24,944) of matched controls were excluded because of incident PPI use within 2 years prior to index date. Finally, the study population comprised 23,321 cases of pancreatic cancers, and 75,937 matched population controls, with a mean number of 3.3 controls per case (Fig. 1).
The characteristics of cases and controls are shown in Table 1. Among cases, mean age at diagnosis was 69.8 ± 10.1 years, and 51.7% were males. Compared with controls, cases had a higher prevalence of diabetes mellitus, tobacco- and alcohol-related diseases or drug use, morbid obesity, and history of acute or chronic pancreatitis. Cases were also more likely to present with other comorbidities such as chronic obstructive pulmonary disease. Pancreatic cancers were most often localized to the head of the pancreas (Table 2). More than half of the cases died within 1 year after the index date (52.3%, N = 12,202). Overall, 77.8% (N = 18,141) of cases and 75.5% (N = 57,307) of controls were PPI ever users, and 43.9% and 37.9% redeemed prescriptions for 181 cDDD or more (Supplementary Table S3). Cases and controls were respectively exposed to 658.3 ± 1,079.1 and 560.8 ± 1,009.7 cDDD in mean during the study period. Omeprazole was the most frequently prescribed drug (in 50.8% of cases, and 47.9% of controls), followed by esomeprazole (46.0% of cases and 40.6% of controls; Supplementary Table S3).
The results of the main analysis are shown in Table 3. Ever use of PPIs was associated with a slightly increased risk of pancreatic cancer when compared with never use [crude OR = 1.15, 95% CI, 1.10–1.19; aOR (final model) = 1.05, 95% CI, 1.01–1.09]. A dose–response relationship was observed (1–30 cDDD: aOR = 0.92, 95% CI, 0.87–0.97; 31–180 cDDD: aOR = 1.05, 95% CI, 1.00–1.11; 181–1,080 cDDD: aOR = 1.18, 95% CI, 1.12–1.24; >1,080 cDDD: aOR = 1.17, 95% CI, 1.10–1.23). Analyses by PPI subtype showed a higher risk with esomeprazole (aOR = 1.18, 95% CI, 1.14–1.22). Similar results were found when covariates included in the model were restricted to potential risk factors of pancreatic cancer only (Table 3).
Stratified analyses are shown in Table 4. The magnitude of the association between PPI use and risk of pancreatic cancer remained consistent across all subgroup analyses. Ever use of PPIs was associated with a significantly increased risk of pancreatic cancer among females (aOR = 1.08, 95% CI, 1.02–1.15), subjects without history of diabetes mellitus (aOR = 1.07, 95% CI,1.02–1.12), without history of tobacco-related disease or drug uses (aOR = 1.05, 95% CI, 1.01–1.10), or without morbid obesity (aOR = 1.05, 95% CI, 1.01–1.10). The associations of these risk factors with pancreatic cancer are presented in Supplementary Table S4 (history of diabetes mellitus: aOR = 2.07, 95% CI, 1.99–2.16; tobacco-related diseases or drug use: aOR = 1.35, 95% CI, 1.28–1.42; morbid obesity: aOR = 1.00, 95% CI, 0.95–1.06).
In sensitivity analyses, a dose–response relationship persisted after introduction of a 2-year or a 4-year lag period on PPI exposure, although the associations were of lower magnitude compared with the main analyses. Statistically significant associations were still observed, above 180 cDDD, with the 2-year lag period (2-year lag analyses: 181–1,080 cDDD: aOR = 1.08, 95% CI, 1.03–1.14; >1,080 cDDD: aOR = 1.11, 95% CI, 1.05–1.19; 4-year lag analyses: 181–1,080 cDDD: aOR = 1.02, 95% CI, 0.97–1.07; >1,080 cDDD: aOR = 1.08, 95% CI, 1.01–1.14; Supplementary Table S5). Restriction to new PPI users (Supplementary Table S6), or IPTW approach (Supplementary Table S7) produced results consistent with those of the main analysis. We observed an increased risk of pancreatic cancer associated with PPI use compared with H2RA use, more marked at high PPI cumulative doses (181–1,080 cDDD: crude OR = 1.43, 95% CI, 1.06–1.92; aOR = 1.15, 95% CI, 0.84–1.56; >1,080 cDDD: crude OR = 1.55, 95% CI, 1.15–2.08; aOR = 1.15, 95% CI, 0.84–1.57; Supplementary Table S8).
Discussion
Principal findings
To our knowledge, the current study is the largest investigation on the risk of pancreatic cancer associated with PPI use, with 23,321 cases included. More than 3 of 4 individuals were PPI users over the study period, with a large proportion exposed to high cumulative doses (>180 cDDD). PPI use was associated with a slightly increased risk of pancreatic cancer, especially for cumulative exposure over 180 DDD. Overall, the results remained robust across subgroups, and in sensitivity analyses.
Comparison with the literature
Three previous observational studies found no association between PPI use and pancreatic cancer (11, 13, 17). Limited power for analyses, or low proportion of long-term PPI users may have explained these null findings. In contrast, five studies reported increased risks, three of them conducted in Asian countries (14, 16, 18), while two set in European countries (12, 15). However, these studies also had limitations. First, regional specificities in the distribution of pancreatic cancer risk factors or patterns of PPI use preclude generalization of their findings (14, 16, 18). In a study conducted in Taiwan (18), the prevalence of viral hepatitis was 10-fold higher than those observed in our study. Second, another study (15) found disproportionate numbers of short-term PPI users in the case group compared with the controls, leading to concerns of reverse causality. Finally, some of these studies were unable to capture major confounders such as tobacco smoking, obesity, or pancreatitis (12, 15, 16, 18). In the present study, we sought to address such limitations through careful study design and various sensitivity analyses.
We observed higher risks of pancreatic cancer among long-term PPI users, or with esomeprazole, one of the most potent PPI in decreasing gastric acidity (36). These findings were consistent with the physiopathology of PPIs described in the literature. Hypergastrinemia, produced as a negative feedback of prolonged PPI use, might stimulate the overgrowth of pancreatic cells via CCK-B/gastrin like receptors. However, although exogenous administration of gastrin promotes pancreatic cancer in animal models, in humans, underlying factors are needed to reactivate CCK-B/gastrin like receptors reexpression from their postnatal silenced state to active state in cancer (37). PPI induced hypochlorhydria can also lead to major changes in the gut microbiome, with consequent potential retrograde microbe migration from the gastrointestinal tract, and modulation of the intratumor microbiome. There is strong evidence for the role of the gut and tumor microbiome in pancreatic cancer, that may impact pancreatic carcinogenesis, progression and resistance to therapy (20, 38).
Strengths and limitations
Our study has a number of strengths. First, it was based on a nationwide database, with comprehensive sociodemographic and medical information on both outpatient and inpatient data, recorded since 2006. This allowed the inclusion of more than 23,000 pancreatic cancers over a 5-year period. Second, this database is a valuable tool for detecting cancers, with expected good predictive value and sensitivity (39, 40), which has been used in several studies (41–45). To identify only primary pancreatic cancers, but not pancreatic metastases or secondary pancreatic cancers, we excluded patients with a history of all causes cancers before the index date, which accounted for about one-third of cases. Nevertheless, one-fifth of diagnoses that led to these exclusions were suggestive of misclassified pancreatic cancers (namely, ICD-10 codes D01: Carcinoma in situ of other and unspecified digestive organs; D37: Neoplasm of uncertain or unknown behavior of oral cavity and digestive organs; C24: Malignant neoplasm of other and unspecified parts of biliary tract), half of them identified in the 2 months preceding the index date. Consequently, our selection procedure was very conservative. However, there is no reason to believe that this could have biased the results or prevented their generalization. Moreover, we found that about 9 of 10 of all pancreatic cancers were adenocarcinoma, most often localized in the head of the pancreas, which is consistent with the epidemiology of the disease (21, 23). Third, many covariates were available in the SNDS, and could have been taken into account in the analyses. Among them, smoking is a major risk factor for pancreatic cancer. In this study, prevalence of tobacco use was consistent with the figures of daily smoking reported within the same age groups in a national survey (46). Moreover, the magnitude of the association with pancreatic cancer was comparable with those of a meta-analysis of 82 studies (47). Fourth, PPI exposure could have been measured during a period of up to 13 years (2006–2018), in a time frame compatible with the development of pancreatic cancer. Finally, the careful implementation of the study design, and numerous sensitivity analyses contributed to the robustness of our results. The case–control and the cohort design are two observational designs relevant for studying drug–cancer associations, with similar underlying concepts (48). Here, a cumulative dose–response investigation was needed for establishing plausibility of a causal effect. Thus, the case–control approach was privileged to compute the exposure level of cases and controls. The results were consistent across several sensitivity analyses. Notably, an increased risk of pancreatic cancer was also observed with PPI compared with H2RA use, suggesting that confounding by indication was likely to be limited.
Our study also has some limitations. Given its observational nature, it is prone to bias (48, 49), including residual confounding, time-related bias, and misclassification of exposure. Residual confounding may have occurred, first, because information on genetic, family history, lifestyle, and environmental risk factors for cancers was not available. Analyses, though, were adjusted for other identifiable potential risk factors of pancreatic cancers, including diabetes mellitus, tobacco, and morbid obesity. Second, the indication for PPI treatment was not recorded in the databases, and thus could not be taken into account in the analyses. Third, the lack of an active comparator in the main analyses may also have led to residual confounding (34). The results of the sensitivity analysis considering H2RAs as an active comparator must be interpreted cautiously, because H2RA use is restricted to a small number of users with specific profiles in France (50). However, they support the finding of an excess risk of pancreatic cancer development associated with PPI exposure as compared to H2RAs. Time-related biases were limited by design. Exclusion of cases and controls with an observation period under 7 years resulted in similar duration of exposure opportunity time, minimizing time-window bias. Nevertheless, even studies with similar observation periods between cases and controls can in some instances, introduce differential drug-treated time window, when the duration of treated disease is different (49). Here, information on the nature, and onset date of the condition that led to the initiation of PPI therapy was not available. Thus, time-window bias cannot be fully ruled out. However, given the careful selection of controls and their matched index dates, such a bias is likely to be limited if any. We employed a very conservative method to address latency time bias and reverse causality (or protopathic bias; ref. 30), excluding new PPI users in 2 years before the index date in the main analysis, and applying 2-year and 4-year lag-times in sensitivity analyses. The lagged analyses tended to decrease the magnitude of the associations with the highest cDDD categories. This may either suggest potential residual reverse causality, or reflect excessive caution in the choice of the delay. Using very long lag periods could tend to unjustifiably consider lower level of exposures, which translates to lower ORs in the case of a dose–response association (51). Finally, potential misclassification of exposure status may have occurred. Cumulative exposure to PPIs was estimated on the basis of the quantity redeemed, but there is no guarantee on patient's adherence to the prescription, or even that the patient actually took the drug. This bias is not expected to affect long-term users, with regularly redeemed prescription. Otherwise, the SNDS does not contain information neither on inpatient nor on OTC PPI use. However, these uses are quantitatively much lower than outpatient use. Moreover, in this study, rates of inpatient or OTC PPI uses were not supposed to be different between cases and controls nor to introduce differential bias.
Conclusion
On the basis of these findings, a slight increase in the risk of pancreatic cancers associated with the use of PPIs at high cumulative doses cannot be excluded. Given the massive PPI use, even a relatively modest association would have important public health implications. Therefore, efforts should be continued to limit PPI treatments to appropriate indications and durations. Regular monitoring and reevaluation of treatment are needed.
Authors' Disclosures
M. Camus reports personal fees from Cook Medical, Boston Scientific, Ambu, and Alfasigma outside the submitted work. J. Kirchgesner reports personal fees from Pfizer, Roche, and Gilead outside the submitted work. No disclosures were reported by the other authors.
Authors' Contributions
M. Lassalle: Conceptualization, software, formal analysis, methodology, writing–original draft, writing–review and editing. T. Le Tri: Software, writing–review and editing. P. Afchain: Methodology, writing–review and editing. M. Camus: Methodology, writing–review and editing. J. Kirchgesner: Methodology, writing–review and editing. M. Zureik: Conceptualization, supervision, validation, methodology, project administration, writing–review and editing. R. Dray-Spira: Conceptualization, supervision, validation, methodology, project administration, writing–review and editing.
Acknowledgments
The study was performed at EPI-PHARE in cooperation with clinicians. No external funding was received.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Note: Supplementary data for this article are available at Cancer Epidemiology, Biomarkers & Prevention Online (http://cebp.aacrjournals.org/).