- 1Leibniz Institute for Prevention Research and Epidemiology – BIPS, Bremen, Germany
- 2Techniker Krankenkasse – TK, Hamburg, Germany
- 3Institute for Experimental and Clinical Pharmacology and Toxicology, University of Lübeck, Lübeck, Germany
- 4Faculty of Mathematics and Computer Science, University of Bremen, Bremen, Germany
- 5Section of Biostatistics, Department of Public Health, University of Copenhagen, Copenhagen, Denmark
Introduction: Pharmacovigilance is vital for drug safety. The process typically involves two key steps: initial signal generation from spontaneous reporting systems (SRSs) and subsequent expert review to assess the signals’ (potential) causality and decide on the appropriate action.
Methods: We propose a novel discovery and verification approach to pharmacovigilance based on electronic healthcare data. We enhance the signal detection phase by introducing an ensemble of methods which generated signals are combined using Borda count ranking; a method designed to emphasize consensus. Ensemble methods tend to perform better when data is noisy and leverage the strengths of individual classifiers, while trying to mitigate some of their limitations. Additionally, we offer the committee of medical experts with the option to perform an in-depth investigation of selected signals through tailored pharmacoepidemiological studies to evaluate their plausibility or spuriousness. To illustrate our approach, we utilize data from the German Pharmacoepidemiological Research Database, focusing on drug reactions to the direct oral anticoagulant rivaroxaban.
Results: In this example, the ensemble method is built upon the Bayesian confidence propagation neural network, longitudinal Gamma Poisson shrinker, penalized regression and random forests. We also conduct a pharmacoepidemiological verification study in the form of a nested active comparator case-control study, involving patients diagnosed with atrial fibrillation who initiated anticoagulant treatment between 2011 and 2017.
Discussion: The case study reveals our ability to detect known adverse drug reactions and discover new signals. Importantly, the ensemble method is computationally efficient. Hasty false conclusions can be avoided by a verification study, which is, however, time-consuming to carry out. We provide an online tool for easy application: https://borda.bips.eu.
1 Introduction
Evidence of the safety of newly approved drugs is often limited. The pivotal randomized clinical trials (RCTs) are powered to assess efficacy so that the sample sizes are too small to examine rare safety outcomes. Patients in RCTs usually have to fulfill several inclusion criteria and especially vulnerable groups such as pregnant women, elderly, and multi-morbid persons are often excluded or underrepresented. Moreover, patients are followed up very closely under controlled conditions over a limited period of time. Post-market surveillance or pharmacovigilance is, therefore, essential to guarantee the safety of drugs in routine care (Routledge, 1998; World Health Organization, 2002; Härmark and van, 2008; Suling and Pigeot, 2012; Aronson, 2023). Its primary goal is to promptly identify any previously unknown adverse drug reactions (ADRs) associated with drugs already on the market, allowing necessary precautions to be taken to safeguard the population. The post-market surveillance process typically comprises of three key elements or stages, as illustrated in Figure 1A: 1) the use of spontaneous reporting system (SRS) data, 2) a signal detection phase, and 3) the evaluation of the results by a committee of medical experts. In this context, we introduce a novel discovery and verification approach to pharmacovigilance, shown in Figure 1B. We start by considering each of the three components of the conventional approach individually and then propose potential improvements for each step.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1426323/fphar-15-1426323-HTML/image_m/fphar-15-1426323-g001.jpg)
Figure 1. (A) The conventional pharmacovigilance approach involves applying a single detection method to SRS data to generate a single ranking of ADRs pairs. This ranking is subsequently submitted to a committee of medical experts. (B) The approach proposed in this paper entails the application of
Pharmacovigilance traditionally relies on SRS data as the primary data source (Routledge, 1998; World Health Organization, 2002; Suling and Pigeot, 2012; Huang et al., 2014; Bailey et al., 2016; Alomar et al., 2020; Noguchi et al., 2021; Costa et al., 2023; Cutroneo et al., 2024). In the last decades, research has shifted its focus to electronic healthcare (EHC) data which encompass individual patients’ drug prescriptions and medical events over time, in addition to personal details, such as age, sex, region of residence and other factors (Suling and Pigeot, 2012; Coloma et al., 2012; Zorych et al., 2011; Patadia et al., 2015; Li et al., 2015; Pacurariu et al., 2015). EHC data offer several advantages over SRSs: 1) the total count of patients prescribed a drug and/or experiencing an ADR is available, whereas spontaneous reports are only filed when the drug was prescribed and an ADR occurred; 2) issues of under- and over-reporting are less prominent (Schroeder, 1998; van der Heijden et al., 2002); and 3) some EHC databases include follow-up data for several million individuals (Haug and Schink, 2021). There is a strong push to make electronic healthcare data more accessible. One such initiative is the European Health Data Space (EHDS; (European Commission, 2024)), which aims to standardize data formats, enhance interoperability, and ensure privacy protections across the European Union.
A wide variety of statistical signal detection methods have been proposed for the detection phase in which one aims to identify associations between drugs and ADRs, often referred to as signals (Pacurariu et al., 2015; Simpson, 2011; Arnaud et al., 2017; Dijkstra et al., 2020). These methods encompass a wide spectrum, including disproportionality measures (e.g., reporting odds ratio (Stricker and Tijssen, 1992)), hypothesis tests (e.g., Poisson test (DuMouchel, 1999)), Bayesian shrinkage estimates (e.g., Bayesian confidence propagation neural network or BCPNN (Bate et al., 1998; Norén et al., 2006)), random forests (RFs (Breiman, 2001)), and sparse regression methods like LASSO (Caster et al., 2010; Tibshirani, 1996). Each method assigns a score to each of the drug-ADR pairs in the data. Following signal generation, these scores are utilized to create a ranking of the drug-ADR pairs, reflecting the strength of associations between them.
Each of these signal detection methods possesses distinct strengths and weaknesses. For instance, the widely-used BCPNN excels in cases involving innocent bystanders, in which a drug is mistakenly linked to an ADR since it is frequently co-prescribed with the actual causative drug (Dijkstra et al., 2020). However, it has the limitation of evaluating each drug-ADR pair independently. In contrast, LASSO, while less effective in scenarios with innocent bystanders (Dijkstra et al., 2020), does not have this constraint and can assess multiple ADRs simultaneously, even accommodating potential confounders (Dijkstra et al., 2020; Caster et al., 2010; Tibshirani, 1996). RFs stand out due to their capability to handle non-linear relationships (Breiman, 2001); a feature that sets them firmly apart from the other methods. The choice of signal detection method is, therefore, non-trivial and depends on many factors, as highlighted in previous research (Dijkstra et al., 2020; Candore et al., 2015). In such scenarios, ensemble methods have shown their effectiveness. They have the capacity to combine the strengths of individual classifiers and mitigate some of their limitations (Rokach, 2009; Antczak, 2016; Tabassum and Ahmed, 2016). Furthermore, ensemble methods are particularly adept at handling noisy data (Ho et al., 1994).
After the signal detection phase, the resulting ranking of drug-ADR pairs is presented to a committee of medical experts, which triages the signals and decides on any possible actions, e.g., issue warnings, request label changes or, in extreme cases, recommend to withdraw the license (Alomar et al., 2019).
We propose the alternative approach to pharmacovigilance shown in Figure 1B, with the goal of addressing the challenges outlined earlier. We choose to utilize EHC data which not only offers the previously mentioned advantages but also provides the possibility to execute a subsequent verification phase, as elaborated on later. Second, to leverage the strengths and address the limitations of various signal detection methods, we advocate for employing an ensemble of more than one method, see Figure 1B. The scores generated by these methods result in varying rankings, which are then aggregated into a unified ranking using the Borda count, a method originally developed in the field of electoral mathematics (Emerson, 2011). The Borda count ranking is advantageous as it emphasizes consensus (Saari, 1985; 2022). Furthermore, we extend the options available to the committee by introducing the option to conduct a verification study using pharmacoepidemiological methods. This verification phase allows the committee to assess whether signals persist when standard methods, such as controlling for confounding, are applied.
In conclusion, we define a new approach to pharmacovigilance that utilizes an ensemble of signal detection methods and offers the committee of medical experts the opportunity to verify select signals through a tailored pharmacoepidemiological study. We assess the approach’s applicability by conducting a case study using real-world EHC data, where we apply an ensemble of four established signal detection methods to uncover potential unknown ADRs associated with the direct oral anticoagulant rivaroxaban. We perform a pharmacoepidemiological signal verification study on a subset of interesting signals selected by the committee.
2 Methods
This section outlines the proposed discovery and verification approach in general terms, structured around three phases as shown in Figure 1B: 1) signal detection, 2) signal triage, and 3) signal verification. To prevent any potential confusion, we define certain terms used throughout the paper here, as their interpretations may vary across different fields. Here, we define an adverse drug reaction (ADR) as any individual event that can be reported to a SRS or documented in electronic healthcare records, even if the event is not commonly thought of as a reaction to a drug. A drug-ADR pair denotes a combination of a specific drug and a particular ADR. A signal is a drug-ADR pair that is deemed associated after applying a signal detection method. The identification of a signal is contingent upon the chosen detection method and the threshold used to determine whether the drug-ADR pair should be submitted to the committee. Additionally, we use the term health outcome to denote a collection of various drug-ADR pairs, typically defined by the committee of medical experts.
2.1 Signal detection
Typically, the detection of potential previously-unknown ADRs is carried out using a single signal detection method. This method assigns a score to each drug-ADR pair within the dataset, reflecting the strength of the association between the drug and the ADR in question. These scores may take various forms, such as odds ratios (van Puijenbroek et al., 2002),
We propose the use of an ensemble of more than one signal detection methods. Given the diverse nature of the scores produced by signal detection methods, it is challenging to directly aggregating them into a single score. Instead, we opt for aggregating the rankings generated by the various signal detection methods into a unified ranking using the Borda count method (Emerson, 2011; Saari, 2022; de Borda, 1781). Drawing inspiration from electoral mathematics, the Borda count treats each signal detection method as a “voter” that ranks drug-ADR pairs in order of preference. The Borda count for each drug-ADR pair is calculated as the sum of the individual ranks assigned by the voters/methods. The drug-ADR pair with the lowest Borda count is deemed to exhibit the strongest signal and is positioned at the top, while the pair with the highest Borda count is placed at the bottom, indicating no signal. Tournament-style counting is employed to handle ties (Narodytska and Walsh, 2014). See the Supplementary Material for a formal definition of Borda count ranking.
The essence of the Borda count ranking lies in capturing the “consensus” among the voters/methods, rather than relying on a simple majority (Emerson, 2011; de Borda, 1781). To see this, let us consider an example where the ensemble consists of three methods and the EHC dataset contains five drug-ADR pairs: A, B, C, D and E. Table 1 shows the ranking of each method. For instance, Method 1 assigned D rank 1, corresponding to the strongest signal, B second etc., and ending with A which is deemed to weakest signal. The right column shows the ranking using the Borda method, with the corresponding Borda count shown in brackets. For example, B is ranked second by all methods, resulting in a Borda count of 6. D ranked first by Method 1 and Method 2, but last by Method 3, giving it a Borda count of
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1426323/fphar-15-1426323-HTML/image_m/fphar-15-1426323-t001.jpg)
Table 1. Example of Borda count ranking with three signal detection methods and five drug-ADR pairs (A, B, C, D,and E).
2.2 Signal triage
Following the signal detection phase, the committee of medical experts is provided with the drug-ADR pairs ranking derived from the Borda count. This committee is entrusted with the task of filtering out signals, aiming to identify severe and previously unknown ADRs and taking appropriate actions.
It is important to note that expert knowledge in this part of the process is paramount. Subject-matter knowledge is required to assess the biological plausibility that an ADR is indeed caused by the drug currently under suspicion. Relying solely on the top-ranked drug-ADR pairs proves insufficient, as these are likely to encompass known ADRs as well, whereas previously unknown ADRs are more likely to be further down the list, i.e., they tend to generate weaker signals. The possibility of the drug in question to be an innocent bystander (Dijkstra et al., 2020), or whether the signal is due to indication bias, has to be ruled out.
Based on medical or pharmacological considerations, the committee has the option to merge specific ADRs. This involves aggregating ICD codes from multiple signals into a single health outcome. See Section 3.3 for an example. A formal definition of a health outcome can be found in the Supplementary Material.
To facilitate the signal triage process, we have developed an online tool based on Borda count rankings, accessible at https://borda.bips.eu. Our R-Shiny app computes the Borda count ranking using uploaded data and offers graphical visualizations to aid the selection of drug-ADR pairs for further investigation.
2.3 Signal verification
The presented approach includes a verification phase (see Figure 1B), providing the committee with the opportunity to assess the reliability of filtered signals through an observational study focused on the specific drug-ADR pair or health outcome of interest. This study should be designed to minimize bias and control for potential confounding, for instance, by employing a new user active comparator approach (Lund et al., 2015), evaluating all potential confounders, and utilizing appropriate statistical methods. It is important to note that due to the substantial number of drug-ADR pairs, systematically verifying all signals through traditional pharmacoepidemiological studies is not feasible. As a result, verification is limited to a small subset of drug-ADR pairs/health outcomes.
3 Case study
We demonstrate the applicability of the approach outlined in the previous section using a case study. We begin by introducing the data source and the drug of interest, rivaroxaban (RVX, ATC: B01AF01). We are focusing on this drug in particular because it is central to the PV-Monitor project from which this work stems. For further details, please refer to the Funding section. The subsequent part of the section adheres to the same structure as the preceding one, starting with signal detection, followed by signal triage, and concluding with signal verification. The results of the case study are available at https://borda.bips.eu as well. The website also provides the option to conduct a similar analysis with one’s own dataset.
3.1 Data source and exposure of interest
We used the German Pharmacoepidemiological Research Database (GePaRD) that is based on claims data from four statutory health insurance (SHI) providers in Germany. It currently includes information on approximately 25 million persons who have been insured with one of the participating providers since 2004 or later (Haug and Schink, 2021).
To illustrate our methodology, we focus on a single drug: RVX, a direct oral anticoagulant (DOAC) that has been approved for, e.g., the prevention of stroke and systemic embolism in patients with atrial fibrillation, the treatment of deep vein thrombosis and pulmonary embolism, and prevention of recurrent deep vein thrombosis and pulmonary embolism in adult patients, and for the prevention of venous thromboembolism in adult patients undergoing elective hip or knee replacement surgery. As is the case with other anticoagulants, clinical studies of rivaroxaban identified hemorrhage as an important safety outcome (Ruff et al., 2014). Since the approval of RVX on the German market in 2008, numerous ADRs have been reported to spontaneous reporting systems (Morgovan et al., 2023). In this case study, however, we are interested in as of yet unknown ADRs. GePaRD and the cohort of the signal detection study are described in detail in the Supplementary Material.
3.2 Signal detection
3.2.1 General principle
We employ four signal detection methods as the foundation for the ensemble. These methods consist of the BCPNN (Norén et al., 2006)), the LGPS (Schuemie, 2011), LASSO (Caster et al., 2010; Tibshirani, 1996), and RFs (Breiman, 2001). See the Supplementary Material for a full description of all four methods. We opt for the BCPNN method, since it is the method of choice for SRSs. The LGPS is a disproportionality method tailored specifically to longitudinal data, demonstrating superior performance in the Observational Medical Outcomes Partnership (OMOP) competition (Schuemie, 2011; The Observational Medical Outcomes Partnership, 2010). The inclusion of LASSO is motivated by its ability to regress all ADRs against the drug of interest while considering patient information and accounting for potential confounding variables. RFs are included in our analysis since they can handle non-linear effects (Breiman, 2001). We use the corrected version of the impurity measure for the RFs (Nembrini et al., 2018). It is important to emphasize that this serves as an illustrative example, and alternative signal detection methods could be employed as well.
We included persons in the signal detection study who.
1. Had been insured with the health insurance Techniker Krankenkasse (TK), which is one of the largest nationwide SHI providers in Germany, for at least 12 consecutive months between January 2015 and December 2016;
2. Who had at least one dispensation of RVX at least 90 days after cohort entry and before cohort exit, respectively, and
3. Who had no DOAC use within 12 months preceding cohort entry.
In our analysis, we incorporate covariates such as age, sex, other diagnoses, and prescribed medications to address potential interactions. See the Supplementary Material for details.
3.2.2 Results
The dataset included 3,795 ICD-codes that were treated as potential ADRs. Figure 2 illustrates the Kendall’s
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1426323/fphar-15-1426323-HTML/image_m/fphar-15-1426323-g002.jpg)
Figure 2. Kendall’s
For instance, consider the ADR acute cystitis with the ICD-code N30.0, which is ranked quite disparately by these methods. BCPNN, LGPS, and RF position this ICD-code among the top 2% of the strongest signals, whereas LASSO does not even include it in the top 10%. Depending solely on the LASSO output could result in overlooking this ADR as a potential signal.
Table 2 presents the top 10 ICD-codes with the highest rankings according to the Borda count. It is worth noting that K92.2 (gastrointestinal hemorrhage, unspecified) is a recognized ADR associated with DOACs.
3.3 Signal triage
3.3.1 General principle
A committee of pharmacoepidemiologists, pharmacologists, physicians and statisticians reviewed the rankings and discussed plausibility, severity and the novelty of the signals.
3.3.2 Results
Based on the ranking obtained during the signal detection phase, the committee identified four health outcomes of interest: acute liver injury (ALI), acute cystitis (CYS), epilepsy and seizures (EPI), and sepsis. Additionally, two health outcomes were chosen as positive controls – gastrointestinal bleeding (GB) and intracranial bleeding (ICB) – as they are known adverse outcomes of RVX. See the Supplementary Material for the ICD codes corresponding to each health outcome.
This selection was made based on medical plausibility, whether the event was already known before, and with the help of plots as shown in Figure 3. This figure shows the relative Borda ranks for the ICD-codes associated with each of the six health outcomes and a negative control: Fracture of head and neck of femur (S72.0). Each dot represents an individual ICD-code. A relative rank of 1 indicates the ADR is ranked last (indicating no signal), while a relative rank of 0 signifies the item is at the top of the ranking (indicating a strong signal). We prefer to use the relative rank for its clarity.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1426323/fphar-15-1426323-HTML/image_m/fphar-15-1426323-g003.jpg)
Figure 3. The relative ranks represent the positions of seven adverse outcomes, with five consisting of multiple individual ADRs. Each dot represents an individual ICD-code. A relative rank of 1 signifies the highest rank, indicating no signal, while a relative rank of 0 indicates the top rank, suggesting a very strong signal.
GB and ICB are established adverse outcomes of RVX, with many associated ADRs exhibiting strong signals. Acute cystitis, corresponding to a single ADR (N30.0), produces a very strong signal. The ADRs linked to epilepsy and sepsis exhibit significant variation; there are both groups of ADRs that display strong associations and those that suggest there is no connection with RVX. Especially the latter two can benefit from performing a verification study. It is evident from the bottom row that there is no signal for the negative control. The negative control is, therefore, not considered in the verification study.
3.4 Signal verification
3.4.1 General principle
The verification study was designed as an active comparator case-control study nested in a cohort of new users of RVX and phenprocoumon (PPC, ATC: B01AA04, the most frequently used vitamin K antagonist in Germany) with atrial fibrillation and included data from 2011 to 2017. Cases were defined as patients with a diagnosis of the respective outcome of interest. Each case was matched with up to 10 controls by sex, age at index day (
3.4.2 Results
The study cohort of the verification study was based on 97,400 new users of PPC and 71,917 users of RVX. The median age at cohort entry was 73 (IQR: 65, 79) for RVX and 75 (IQR: 69, 81) for PPC new users. The proportion of women was 47% (RVX) and 48% (PPC). We observed 5,053 cases of acute cystitis, 322 cases of acute liver injury, 3,070 cases of epilepsy, 3,504 cases of sepsis, 6,705 cases of gastrointestinal bleeding, and 2,974 cases of intracranial bleeding. A full description of the study cohort, including baseline characteristics, relevant medical history and dispensed medication at cohort entry, can be found in the Supplementary Material.
Table 3 shows that the adjusted risk of gastrointestinal bleeding was increased in current single use of RVX compared to current single use of PPC (OR = 1.37, 95% CI: (1.27; 1.47)) but decreased for intracranial bleeding (OR = 0.74, 95% CI: (0.66; 0.83)). This is what we expected from the literature (Hohnloser et al., 2017). The risk of cystitis was slightly increased (OR = 1.11, 95% CI: (1.02; 1.21)) whereas the risk of acute liver injury was much lower in patients using RVX compared to PPC (OR = 0.50, 95% CI: (0.35; 0.72)). Patients receiving current treatment with RVX showed a considerably higher risk of epilepsy and seizures compared to PPC users (OR = 1.26, 95% CI: (1.11; 1.43)). The risk of sepsis was similar in current users of RVX compared to current use of PPC (OR = 1.10, 95% CI: (0.99; 1.23)).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1426323/fphar-15-1426323-HTML/image_m/fphar-15-1426323-t003.jpg)
Table 3. Results of the verification study: matched crude and adjusted odds ratios with corresponding 95% confidence intervals.
4 Discussion and conclusion
In this paper, we introduced a novel approach to pharmacovigilance. Conventionally, a single signal detection method is applied to SRS data, and the resulting ranking of drug-ADR pairs is presented to a committee of medical experts. Here, we enhance this process in three significant ways: First, we opt for EHC data due to its richness and its capacity to facilitate a pharmacoepidemiological study; something that cannot be accomplished with SRS data. Second, we simultaneously apply multiple signal detection methods, aggregating the obtained rankings through the Borda count and we provide the committee with a visual representation for defining health outcomes. See, for example, Figure 3. Third, we recommend verifying a selection of health outcomes in a pharmacoepidemiological verification study. Such a study can be immensely beneficial, as it may offer support for either confirming the existence of the association or revealing its spurious nature (Platzbecker et al., 2023).
To illustrate this approach, we conducted a case study with the direct oral anticoagulant rivaroxaban (RVX) based on the healthcare claims database GePaRD. In our case study, the ensemble method comprised of four signal detection methods: BCPNN, LGPS, LASSO, and RF. After applying the ensemble method to the data, the resulting ranking was sent to the committee. During signal triage, a total of four health outcomes were selected to be further investigated in a pharmacoepidemiogical verification study: acute liver injury (ALI), acute cystitis (CYS), epilepsy and seizures (EPI), and sepsis. In addition, gastrointestinal bleeding (GB) and intracranial bleeding (ICB) served as positive controls, being recognized adverse health outcomes of RVX. We performed an active comparator, case-control study nested in GePaRD in which we compared patients with atrial fibrillation who initiated treatment with RVX or phenprocomoun (PPC). For each health outcome, possible confounding variables were collected from the literature and adjusted for. Our approach successfully identified both the positive controls GB and ICB. Additionally, the results indicated no significant safety issues in RVX users compared to PPC concerning the occurrences of ALI, CYS, and sepsis, which were ranked in the top third of the Borda count ranking, as depicted in Figure 3. However, the verification study confirmed an elevated risk for epilepsy and seizures in RVX users compared to PPC users. This finding was further investigated and published by Platzbecker et al. (2023). Within this paper, also the limitations of the pharmacoepidemiological studies are discussed in more detail. These examples illustrate that signal detection using EHC data is feasible, but also underscore the importance of verification to address biases such as confounding, selection, or measurement bias (Prada-Ramallal et al., 2019).
The ensemble method proposed for the signal detection phase is characterized by its flexibility, allowing for easy modification of the signal detection methods incorporated in the ensemble, and it is computationally fast. By employing an ensemble, the goal is to leverage the individual advantages of different methods while potentially mitigating some of their limitations (Rokach, 2009; Antczak, 2016; Tabassum and Ahmed, 2016). Moreover, ensemble techniques are adept at handling noisy data (Ho et al., 1994). The use of Borda count ranking for aggregation offers the advantage of independence from the interpretation or characteristics of the scores generated by each individual signal detection method. The method aims to promote consensus (Emerson, 2011; Saari, 2022), as discussed in Section 2. The case study highlights that individual signal detection methods may produce significantly different rankings, as shown in Figure 2. The Borda method also allows for assigning weights to the rankings produced by individual signal detection methods. This adaptability facilitates the prioritization of output from specific methods over others, proving particularly valuable when there is an expectation that a particular method will excel in the given scenario.
Conducting a pharmacoepidemiological verification study is generally laborious and if the study confirms a safety issue, a follow-up with sensitivity analyses to investigate different sources of potential bias is required. Hence, we considered only a limited number of health outcomes in our case study. It is important to recognize that the decision for the additional verification step in the presented approach leads to a time-consuming process. While this may be acceptable in certain situations, the committee tasked with evaluating signals must factor in this aspect when deciding on the course of action.
In the case study, we focused on a single medication rather than multiple drugs simultaneously, which is more typical in pharmacovigilance. However, this does not prohibit the application of the proposed approach to scenarios involving multiple drugs. The resulting Borda ranking presented to the committee would be longer; however, this is standard in the pharmacovigilance context. Most signal detection methods are computationally efficient and capable of handling large datasets with numerous drug-ADR pairs. Consequently, running multiple methods simultaneously should not present significant challenges.
It is standard practice to partition the dataset into two segments for detection and verification to avoid overfitting. Nevertheless, when the available data are limited, such as when the drug of interest is seldom prescribed or the ADR is rare, this approach might significantly diminish the power to detect relevant signals, possibly leaving crucial associations between drugs and ADRs unnoticed. In such scenarios, the option of not splitting the data could be considered; however, we caution that this decision should be taken into account when interpreting the results of the verification study. In our RVX case study, the datasets used for signal detection and the verification study overlap slightly (2 years of data from one health insurance company compared to data over 7 years from four health insurance companies). In addition, the inclusion and exclusion criteria for signal generation and the verification studies were very different. For more details, please see Section 3.2 and Section 3.4. In Platzbecker et al. (2023), we conducted a sensitivity analysis by repeating our main analysis without the data used for signal detection and were able to confirm our results.
The availability of EHC data is often delayed, meaning that there can be a significant gap between the occurrence of the event itself and the moment the observation becomes available for analysis. This delay must be taken into account when utilizing EHC data for detection, as it can be challenging, if not impossible, to detect ADRs for drugs that very recently entered the market. In those cases, SRS data might be a more appropriate choice, since reports tend to be processed quickly. However, for drugs that have been on the market for a longer time, EHC data may be preferable due to its richer nature, allowing for easier detection of rare ADRs and long-term effects. Therefore, we recommend users to carefully consider which data source to use for signal detection, as it heavily depends on their specific goals.
In efforts to improve post-market surveillance, research typically focuses on one or more of three key aspects: 1) the data source being used, 2) the statistical methods employed for signal detection, and 3) the signal triage process. Alongside traditional data sources such as SRSs or EHC data, alternative sources like social media posts, journal publications, sensor data, and reporting apps are also being explored (Lavertu et al., 2021; Satwika et al., 2021; Shin et al., 2022; Pilipiec et al., 2022). Instead of relying solely on a single source, some studies delve into multimodal approaches, integrating multiple sources simultaneously (Dimitriadis et al., 2023). In this particular study, our choice for the signal detection phase was EHC data since it can be used for the verification study as well. However, other data sources could also be used. However, alternative data forms such as SRS and social media are unsuitable for performing a pharmacoepidemiological verification study.
A considerable amount of research is directed towards enhancing the signal detection phase through the development of novel methodologies or the refinement of existing ones. Particularly in recent years, there has been a strong emphasis on artificial intelligence methods within the field (Bate and Hobbiger, 2021; Huysentruyt et al., 2021; Ball and Dal Pan, 2022; Kompa et al., 2022; Painter et al., 2023), although classical methods continue to be considered as well. For example, Courtois et al. (2021) proposed employing the adaptive LASSO to simplify the selection of tuning parameters. Note that while it is feasible to integrate multiple methods simultaneously within our framework, determining the optimal set of methods to use requires further investigation.
Hsieh et al. (2023) conducted a study that employs two signal detection methods simultaneously (Hsieh et al., 2023): sequence symmetry analysis (SSA; (Lai et al., 2017)) and tree-based scan statistics (TreeScan; (Kulldorff et al., 2013)). In their research, they examine two drugs using EHC data. Any ADR flagged by either SSA or TreeScan is treated as a signal. Consequently, their approach requires one to select appropriate thresholds for both methods, which can be challenging. In contrast, our approach simplifies this process by requiring only the selection of a threshold for the Borda count. How to extend their approach to accommodate an arbitrary number of methods and drugs is unclear.
Furthermore, research addresses the third phase, signal triage, by aiming to further refine the list of signals by filtering out known signals. In the same study by Hsieh et al. (2023), signals are automatically categorized into various groups, e.g., known ADRs, ADRs associated with indications and ADRs linked to patient characteristics. However, this study only considers two drugs. Subsequent research is need to explore how to extend this automated triage process to encompass more drugs. Thus far, we have not encountered any research proposing to enhance the post-market surveillance process with a verification step.
To facilitate the implementation of the approach described in this paper, we have developed an online tool accessible at https://borda.bips.eu.
Data availability statement
The data analyzed in this study is subject to the following licenses/restrictions: As we are not the owners of the data we are not legally entitled to grant access to the data of the German Pharmacoepidemiological Research Database. In accordance with German data protection regulations, access to the data is granted only to BIPS employees on the BIPS premises and in the context of approved research projects. Third parties may only access the data in cooperation with BIPS and after signing an agreement for guest researchers at BIPS. Requests to access these datasets should be directed to gepard@leibniz-bips.de.
Ethics statement
In Germany, the utilization of health insurance data for scientific research is regulated by the Code of Social Law. All involved health insurance providers as well as the German Federal (Social) Insurance Office and the Senator for Science, Health, and Consumer Protection in Bremen as their responsible authorities approved the use of GePaRD data for this study. Informed consent for studies based on GePaRD is not required by law, and according to the Ethics Committee of the University of Bremen, these studies are exempt from institutional review board review.
Author contributions
LD: Conceptualization, Formal Analysis, Investigation, Methodology, Software, Visualization, Writing–original draft, Writing–review and editing. TS: Conceptualization, Investigation, Methodology, Writing–review and editing. RL: Funding acquisition, Investigation, Resources, Writing–review and editing. MS: Funding acquisition, Investigation, Writing–review and editing. IP: Funding acquisition, Investigation, Resources, Supervision, Writing–review and editing. MW: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Writing–original draft, Writing–review and editing. RF: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This study was supported by the innovation fund (“Innovationsfonds”) of the Federal Joint Committee in Germany (grant number: 01VSF16020). MNW is supported by the German Research Foundation (DFG) – Emmy Noether Grant 437611051.
Acknowledgments
The authors would like to thank all statutory health insurances which provided data for this study, namely AOK Bremen/Bremerhaven, DAK-Gesundheit, Die Techniker (TK), and hkk Krankenkasse. The authors would also like to thank Arne Godt for the implementation of the online tool.
Conflict of interest
Author RL was employed by Techniker Krankenkasse – TK.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
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/fphar.2024.1426323/full#supplementary-material
References
Alomar, M., Palaian, S., and Al-tabakha, M. M. (2019). Pharmacovigilance in perspective: drug withdrawals, data mining and policy implications. F1000Research 8, 2109. doi:10.12688/f1000research.21402.1
Alomar, M., Tawfiq, A. M., Hassan, N., and Palaian, S. (2020). Post marketing surveillance of suspected adverse drug reactions through spontaneous reporting: current status, challenges and the future. Ther. Adv. Drug Saf. 11, 2042098620938595. doi:10.1177/2042098620938595
Antczak, K. (2016). Rank thresholds in classifier ensembles in medical diagnosis. Comput. Sci. Math. Model. 0, 5–12. doi:10.5604/01.3001.0009.4412
Arnaud, M., Bégaud, B., Thurin, N., Moore, N., Pariente, A., and Salvo, F. (2017). Methods for safety signal detection in healthcare databases: a literature review. Expert Opin. Drug Saf. 16, 721–732. doi:10.1080/14740338.2017.1325463
Aronson, J. K. (2023). When I use a word. . . Medical definitions: pharmacovigilance signals. BMJ 381, 1032. doi:10.1136/bmj.p1032
Bailey, C., Peddie, D., Wickham, M. E., Badke, K., Small, S. S., Doyle-Waters, M. M., et al. (2016). Adverse drug event reporting systems: a systematic review. Br. J. Clin. Pharmacol. 82, 17–29. doi:10.1111/bcp.12944
Ball, R., and Dal Pan, G. (2022). Artificial intelligence for pharmacovigilance: ready for prime time? Drug Saf. 45, 429–438. doi:10.1007/s40264-022-01157-4
Bate, A., and Hobbiger, S. F. (2021). Artificial intelligence, real-world automation and the safety of medicines. Drug Saf. 44, 125–132. doi:10.1007/s40264-020-01001-7
Bate, A., Lindquist, M., Edwards, I. R., Olsson, S., Orre, R., Lansner, A., et al. (1998). A Bayesian neural network method for adverse drug reaction signal generation. Eur. J. Clin. Pharmacol. 54, 315–321. doi:10.1007/s002280050466
Candore, G., Juhlin, K., Manlik, K., Thakrar, B., Quarcoo, N., Seabroke, S., et al. (2015). Comparison of statistical signal detection methods within and across spontaneous reporting databases. Drug Saf. 38, 577–587. doi:10.1007/s40264-015-0289-5
Caster, O., Norén, G. N., Madigan, D., and Bate, A. (2010). Large-scale regression-based pattern discovery: the example of screening the WHO global drug safety database. Stat. Analysis Data Min. ASA Data Sci. J. 3, 197–208. doi:10.1002/sam.10078
Coloma, P. M., Avillach, P., Salvo, F., Schuemie, M. J., Ferrajolo, C., Pariente, A., et al. (2012). A reference standard for evaluation of methods for drug safety signal detection using electronic healthcare record databases. Drug Saf. 36, 13–23. doi:10.1007/s40264-012-0002-x
Costa, C., Abeijon, P., Rodrigues, D. A., Figueiras, A., Herdeiro, M. T., and Torre, C. (2023). Factors associated with underreporting of adverse drug reactions by patients: a systematic review. Int. J. Clin. Pharm. 45, 1349–1358. doi:10.1007/s11096-023-01592-y
Courtois, É., Tubert-Bitter, P., and Ahmed, I. (2021). New adaptive Lasso approaches for variable selection in automated pharmacovigilance signal detection. BMC Med. Res. Methodol. 21, 271. doi:10.1186/s12874-021-01450-3
Cutroneo, P. M., Sartori, D., Tuccori, M., Crisafulli, S., Battini, V., Carnovale, C., et al. (2024). Conducting and interpreting disproportionality analyses derived from spontaneous reporting systems. Front. Drug Saf. Regul. 3, 1323057. doi:10.3389/fdsfr.2023.1323057
Dijkstra, L., Garling, M., Foraita, R., and Pigeot, I. (2020). Adverse drug reaction or innocent bystander? A systematic comparison of statistical discovery methods for spontaneous reporting systems. Pharmacoepidemiol. Drug Saf. 29, 396–403. doi:10.1002/pds.4970
Dimitriadis, V. K., Dimitsaki, S., Chytas, A., Gavriilidis, G. I., Kakalou, C., Bonotis, P., et al. (2023). An open-source platform integrating emerging data sources to support multi-modal active pharmacovigilance. Front. Drug Saf. Regul. 2, 1016042. doi:10.3389/fdsfr.2022.1016042
DuMouchel, W. (1999). Bayesian data mining in large frequency tables, with an application to the FDA spontaneous reporting system. Am. Statistician 53, 177–190. doi:10.1080/00031305.1999.10474456
Emerson, P. (2011). The original Borda count and partial voting. Soc. Choice Welf. 40, 353–358. doi:10.1007/s00355-011-0603-9
European Commission (2024). European health data Space. Available at: https://health.ec.europa.eu/ehealth-digital-health-and-care/european-health-data-space_en (Accessed July 22, 2024).
Goldstein, L., and Langholz, B. (1996). Risk set sampling in epidemiologic cohort studies. Stat. Sci. 11. doi:10.1214/ss/1032209663
Härmark, L., and van, G. (2008). Pharmacovigilance: methods, recent developments and future perspectives. Eur. J. Clin. Pharmacol. 64, 743–752. doi:10.1007/s00228-008-0475-9
Haug, U., and Schink, T. (2021). German pharmacoepidemiological research database (GePaRD). Databases Pharmacoepidemiol. Res., 119–124. doi:10.1007/978-3-030-51455-6_8
Ho, T. K., Hull, J. J., and Srihari, S. N. (1994). Decision combination in multiple classifier systems. IEEE Trans. pattern analysis Mach. Intell. 16, 66–75. doi:10.1109/34.273716
Hohnloser, S., Basic, E., and Nabauer, M. (2017). Comparative risk of major bleeding with new oral anticoagulants (NOACs) and phenprocoumon in patients with atrial fibrillation: a post-marketing surveillance study. Clin. Res. Cardiol. 106, 618–628. doi:10.1007/s00392-017-1098-x
Hsieh, M. H.-C., Liang, H.-Y., Tsai, C.-Y., Tseng, Y.-T., Chao, P.-H., Huang, W.-I., et al. (2023). A new drug safety signal detection and triage system integrating sequence symmetry analysis and tree-based scan statistics with longitudinal data. Clin. Epidemiol. 15, 91–107. doi:10.2147/CLEP.S395922
Huang, Y.-L., Moon, J., and Segal, J. B. (2014). A comparison of active adverse event surveillance systems worldwide. Drug Saf. 37, 581–596. doi:10.1007/s40264-014-0194-3
Huysentruyt, K., Kjoersvik, O., Dobracki, P., Savage, E., Mishalov, E., Cherry, M., et al. (2021). Validating intelligent automation systems in pharmacovigilance: insights from good manufacturing practices. Drug Saf. 44, 261–272. doi:10.1007/s40264-020-01030-2
Kompa, B., Hakim, J. B., Palepu, A., Kompa, K. G., Smith, M., Bain, P. A., et al. (2022). Artificial intelligence based on machine learning in pharmacovigilance: a scoping review. Drug Saf. 45, 477–491. doi:10.1007/s40264-022-01176-1
Kulldorff, M., Dashevsky, I., Avery, T. R., Chan, A. K., Davis, R. L., Graham, D., et al. (2013). Drug safety data mining with a tree-based scan statistic. Pharmacoepidemiol. Drug Saf. 22, 517–523. doi:10.1002/pds.3423
Lai, E. C.-C., Pratt, N., Hsieh, C.-Y., Lin, S.-J., Pottegård, A., Roughead, E. E., et al. (2017). Sequence symmetry analysis in pharmacovigilance and pharmacoepidemiologic studies. Eur. J. Epidemiol. 32, 567–582. doi:10.1007/s10654-017-0281-8
Lavertu, A., Vora, B., Giacomini, K. M., Altman, R., and Rensi, S. (2021). A new era in pharmacovigilance: toward real-world data and digital monitoring. Clin. Pharmacol. and Ther. 109, 1197–1202. doi:10.1002/cpt.2172
Li, Y., Ryan, P., Wei, Y., and Friedman, C. (2015). A method to combine signals from spontaneous reporting systems and observational healthcare data to detect adverse drug reactions. Drug Saf. 38, 895–908. doi:10.1007/s40264-015-0314-8
Lund, J. L., Richardson, D. B., and Stürmer, T. (2015). The active comparator, new user study design in pharmacoepidemiology: historical foundations and contemporary application. Curr. Epidemiol. Rep. 2, 221–228. doi:10.1007/s40471-015-0053-5
Morgovan, C., Dobrea, C. M., Chis, A. A., Juncan, A. M., Arseniu, A. M., Rus, L. L., et al. (2023). A descriptive analysis of direct oral anticoagulant drugs dosing errors based on spontaneous reports from the EudraVigilance database. Pharmaceuticals 16, 455. doi:10.3390/ph16030455
Narodytska, N., and Walsh, T. (2014). The computational impact of partial votes on strategic voting. Available at: http://arxiv.org/abs/1405.7714 (Accessed May 28, 2014).
Nembrini, S., König, I. R., and Wright, M. N. (2018). The revival of the Gini importance? Bioinformatics 34, 3711–3718. doi:10.1093/bioinformatics/bty373
Noguchi, Y., Tachi, T., and Teramachi, H. (2021). Detection algorithms and attentive points of safety signal using spontaneous reporting systems as a clinical data source. Briefings Bioinforma. 22, bbab347. doi:10.1093/bib/bbab347
Norén, G. N., Bate, A., Orre, R., and Edwards, I. R. (2006). Extending the methods used to screen the WHO drug safety database towards analysis of complex associations and improved accuracy for rare events. Statistics Med. 25, 3740–3757. doi:10.1002/sim.2473
Pacurariu, A. C., Straus, S. M., Trifirò, G., Schuemie, M. J., Gini, R., Herings, R., et al. (2015). Useful interplay between spontaneous ADR reports and electronic healthcare records in signal detection. Drug Saf. 38, 1201–1210. doi:10.1007/s40264-015-0341-5
Painter, J. L., Kassekert, R., and Bate, A. (2023). An industry perspective on the use of machine learning in drug and vaccine safety. Front. Drug Saf. Regul. 3, 1110498. doi:10.3389/fdsfr.2023.1110498
Patadia, V. K., Schuemie, M., Coloma, P., Herings, R., van der Lei, J., Straus, S., et al. (2015). Evaluating performance of electronic healthcare records and spontaneous reporting data in drug safety signal detection. Int. J. Clin. Pharm. 37, 94–104. doi:10.1007/s11096-014-0044-5
Pilipiec, P., Liwicki, M., and Bota, A. (2022). Using machine learning for pharmacovigilance: a systematic review. Pharmaceutics 14, 266. doi:10.3390/pharmaceutics14020266
Platzbecker, K., Müller-Fielitz, H., Foraita, R., Koepp, M. J., Voss, A., Pflock, R., et al. (2023). In atrial fibrillation epilepsy risk differs between oral anticoagulants: active comparator, nested case-control study. Europace 25, euad087. doi:10.1093/europace/euad087
Prada-Ramallal, G., Takkouche, B., and Figueiras, A. (2019). Bias in pharmacoepidemiologic studies using secondary health care databases: a scoping review. BMC Med. Res. Methodol. 19, 53. doi:10.1186/s12874-019-0695-y
Rokach, L. (2009). Ensemble-based classifiers. Artif. Intell. Rev. 33, 1–39. doi:10.1007/s10462-009-9124-7
Rothman, K. J., Greenland, S., and Lash, T. L. (2008). “Case-control studies,” in Modern Epidemiology. 3rd Edn, Editors K. J. Rothman, S. Greenland, and T. L. Lash (Philadelphia: Lippincott Williams and Wilkins), 111–127.
Routledge, P. (1998). 150 Years of pharmacovigilance. Lancet 351, 1200–1201. doi:10.1016/S0140-6736(98)03148-1
Ruff, C. T., Giugliano, R. P., Braunwald, E., Hoffman, E. B., Deenadayalu, N., Ezekowitz, M. D., et al. (2014). Comparison of the efficacy and safety of new oral anticoagulants with warfarin in patients with atrial fibrillation: a meta-analysis of randomised trials. Lancet 383, 955–962. doi:10.1016/S0140-6736(13)62343-0
Saari, D. G. (2022). Selecting a voting method: the case for the Borda count. Const. Polit. Econ. 34, 357–366. doi:10.1007/s10602-022-09380-y
Satwika, M. V., Sushma, D. S., Jaiswal, V., Asha, S., and Pal, T. (2021). The role of advanced technologies supplemented with traditional methods in pharmacovigilance sciences. Recent Pat. Biotechnol. 15, 34–50. doi:10.2174/1872208314666201021162704
Schroeder, D. (1998). Statistics: detecting a rare adverse drug reaction using spontaneous reports. Regional Anesth. Pain Med. 23, 183–189. doi:10.1016/s1098-7339(98)90145-6
Schuemie, M. J. (2011). Methods for drug safety signal detection in longitudinal observational databases: LGPS and LEOPARD. Pharmacoepidemiol. Drug Saf. 20, 292–299. doi:10.1002/pds.2051
Shin, H., Park, S., and Lee, S. (2022). Adverse drug reaction analysis methods and research trends by data sources for post-marketing surveillance. J. Health Inf. Statistics 47, S41–S50. doi:10.21032/jhis.2022.47.s3.s41
Simpson, S. E. (2011). Self-controlled methods for postmarketing drug safety surveillance in large-scale longitudinal data. PhD thesis. New York City (United States): Columbia University.
Stricker, B. H., and Tijssen, J. G. (1992). Serum sickness-like reactions to cefaclor. J. Clin. Epidemiol. 45, 1177–1184. doi:10.1016/0895-4356(92)90158-j
Suling, M., and Pigeot, I. (2012). Signal detection and monitoring based on longitudinal healthcare data. Pharmaceutics 4, 607–640. doi:10.3390/pharmaceutics4040607
Tabassum, N., and Ahmed, T. (2016). “A theoretical study on classifier ensemble methods and its applications,” in 2016 3rd International Conference on Computing for Sustainable Global Development (INDIACom), New Delhi, India, March 16-18, 2016 (IEEE), 374–378.
The Observational Medical Outcomes Partnership (OMOP) (2010). Research partnership announces competition winners: OMOP cup challenged contestants to develop algorithms to improve drug safety. Available at: https://fnih.org/press-release/research-partnership-announces-competition-winners/ (Accesssed May 24, 2010).
Tibshirani, R. (1996). Regression shrinkage and selection via the Lasso. J. R. Stat. Soc. Ser. B Methodol. 58, 267–288. doi:10.1111/j.2517-6161.1996.tb02080.x
van der Heijden, P. G. M., van Puijenbroek, E. P., van Buuren, S., and van der Hofstede, J. W. (2002). On the assessment of adverse drug reactions from spontaneous reporting systems: the influence of under-reporting on odds ratios. Statistics Med. 21, 2027–2044. doi:10.1002/sim.1157
van Puijenbroek, E. P., Bate, A., Leufkens, H. G. M., Lindquist, M., Orre, R., and Egberts, A. C. G. (2002). A comparison of measures of disproportionality for signal detection in spontaneous reporting systems for adverse drug reactions. Pharmacoepidemiol. Drug Saf. 11, 3–10. doi:10.1002/pds.668
Keywords: Borda count, DOAC, drug safety, GEPARD, machine learning, medical records
Citation: Dijkstra L, Schink T, Linder R, Schwaninger M, Pigeot I, Wright MN and Foraita R (2024) A discovery and verification approach to pharmacovigilance using electronic healthcare data. Front. Pharmacol. 15:1426323. doi: 10.3389/fphar.2024.1426323
Received: 01 May 2024; Accepted: 19 August 2024;
Published: 04 September 2024.
Edited by:
Marco Tuccori, University of Pisa, ItalyReviewed by:
Pantelis Natsiavas, Centre for Research and Technology Hellas (INAB|CERTH), GreeceAditi Shendre, The Ohio State University, United States
Copyright © 2024 Dijkstra, Schink, Linder, Schwaninger, Pigeot, Wright and Foraita. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ronja Foraita, foraita@leibniz-bips.de