Objectives Hepatic enzyme-inducing antiepileptic drugs (AEDs) increase serum lipid levels and other atherogenic markers via the induction of cytochrome P450 and may therefore increase the risk of vascular events. We sought to assess the risk of ischaemic stroke and myocardial infarction (MI) according to AED enzymatic properties.
Design Population-based cohort study with nested case–control analysis.
Setting 650 general practices in the UK contributing to the Clinical Practice Research Datalink.
Participants A cohort of 252 407 incident AED users aged 18 or older between January 1990 and April 2013. For each case of ischaemic stroke or MI, up to 10 controls were randomly selected among the cohort members in the risk sets defined by the case and matched on age, sex, indication for AED, calendar time and duration of follow-up.
Interventions Current use of enzyme-inducing and enzyme-inhibiting AEDs compared with non-inducing AEDs.
Primary outcome measures Incidence rate ratios (RRs) of ischaemic stroke and MI.
Results 5069 strokes and 3636 MIs were identified during follow-up. Inducing AEDs use was associated with a small increased risk of ischaemic stroke (RR=1.16, 95% CI 1.02 to 1.33) relative to non-inducing AEDs, most likely due to residual confounding. However, current use of inducing AEDs for ≥24 months was associated with a 46% increased risk of MI (RR=1.46, 95% CI 1.15 to 1.85) compared with the same duration of non-inducing AED, corresponding to a risk difference of 1.39/1000 (95% CI 0.33 to 2.45) persons per year. Current use of inhibiting AED was associated with a decreased risk of MI (RR=0.81, 95% CI 0.66 to 1.00).
Conclusions The use of enzyme-inducing AEDs was not associated with an increased risk of ischaemic stroke; a small increase of MI with prolonged use was observed. In contrast, use of inhibiting AEDs was associated with a decreased risk of MI.
- Antiepileptic Drugs
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Statistics from Altmetric.com
Strengths and limitations of this study
We conducted a large population-based cohort study using the UK Clinical Practice Research Datalink (CPRD), allowing precise estimates and generalisability of our results.
We estimated rates of vascular events within groups of patients with the same indication for antiepileptic drugs, minimising the potential for confounding by indication.
Results were consistent when using alternative methods to control for confounding.
Exposure to antiepileptic drugs was based on prescriptions issued by physicians and not on prescriptions actually filled or taken by patients.
Antiepileptic drugs (AEDs) are increasingly used to treat conditions other than epilepsy such as migraine, pain originating from the nervous system or bipolar disorders. As such, it is estimated that more than 1% of the general population is exposed to AEDs.1 AEDs can be distinguished according to their action on the liver enzymatic system as inducing AEDs, non-inducing AEDs and inhibiting AEDs.2 Older AEDs, that is, carbamazepine, phenobarbital, phenytoin and primidone are all inducing AEDs, except sodium valproate, which is the only inhibiting AED. Second-generation AEDs have weak or non-inducing properties.2
Through their potent hepatic enzyme-inducing properties, predominantly on the cytochrome P450 system, inducing AEDs may lead to drug interactions and alter various metabolisms, including lipid metabolism.3 Indeed, several studies have shown that adults with epilepsy treated with inducing AEDs have increased serum levels of total cholesterol, low-density lipoprotein (LDL) cholesterol, triglycerides, lipoprotein (a), as well as C reactive protein and homocysteine.4–7 Similar findings were observed in healthy users exposed to carbamazepine4 and these changes are detectable 2–3 months after onset of treatment.4–8 Moreover, switching from inducing AEDs to non-inducing AEDs led to a significant decrease in these markers of atherothrombotic risk.6–9 Since enzyme-inducing AEDs promote proatherogenic factors, concerns have been raised regarding their use as first-line agents in the absence of studies properly assessing their effect on the risk of cardiovascular and cerebrovascular events.10 Conversely, sodium valproate, which possesses inhibiting enzymatic properties, has not been associated with such metabolic changes. However, the net vascular effect of sodium valproate is unclear as valproate could have a proatherogenic effect through induction of insulin resistance, body weight gain, metabolic syndrome and increased oxidative stress.11 ,12
Several cross-sectional studies have shown that patients with epilepsy have a higher prevalence of vascular risk factors and vascular diseases than the general population.13 Considering the increasing use of AEDs, and the long-term exposure of patients in the context of chronic diseases, a thorough investigation of their vascular risk is warranted. The objective of this study was therefore to assess separately the risk of ischaemic stroke and myocardial infarction (MI) associated with the use of inducing and inhibiting AEDs compared with non-inducing AEDs using a population-based cohort of patients newly treated with AEDs.
The source population was identified using the UK's Clinical Practice Research Datalink (CPRD; previously known as General Practice Research Database (GPRD)),14 ,15 the world's largest computerised database of anonymised primary care medical records on a population of over 12 million patients enrolled with more than 650 general practices. The geographic distribution of the practices, as well as age and sex distributions of patients, broadly reflects those of the UK population.14 Information collected includes demographic characteristics, lifestyle factors, medical symptoms and diagnoses, laboratory tests, prescriptions and referrals to specialists and hospitals. Prescriptions dispensed by the general practitioner are automatically recorded on the computerised file. Read codes are used to enter medical diagnoses,16 and a coded drug dictionary based on the UK Prescription Pricing Authority Dictionary is used for recording prescriptions. General practitioners are specifically trained for the recording of data, and practices fulfilling predefined quality criteria are labelled ‘up-to-standard’ (UTS).17 Completeness and quality of data as well as consistency with medical files are checked regularly, and numerous studies have shown the high quality of recorded data,18 ,19 including diagnostic codes for stroke and MI.20
Study design and cohort definition
We conducted a retrospective cohort study that was analysed using a nested case–control approach,21 ,22 separately for the two outcomes. This nested case–control approach is computationally more efficient than a time-dependent survival analysis while producing equivalent estimates.23 ,24 The cohort was formed of all patients in the CPRD aged 18 years or older who were members of a UTS practice, and who received at least one prescription for one of the following AEDs between 1 January 1990 and 30 April 2013: beclamide, carbamazepine, ethosuximide, gabapentin, lacosamide, lamotrigine, levetiracetam, oxcarbazepine, phenobarbital, phenytoin, piracetam, pregabalin, primidone, rufinamide, sodium valproate, tiagabine, topiramate, vigabatrin and zonisamide. Cohort entry (time zero) was taken as the date of the first AED prescription within the study period. Patients with less than 1 year of information in the database prior to the date of cohort entry were excluded, as were patients who had been dispensed an AED before cohort entry in order to include only new users of AEDs.25 Patients with a history of stroke at any time before cohort entry were excluded when studying stroke, whereas patients with a history of MI were excluded when studying MI. Patients were followed until the date of their first stroke or MI, depending on the outcome under study, the date they transferred out of the UTS practice, death or end of the study period (30 April 2013), whichever occurred first.
For both outcomes, we identified all participants within our cohort with a first-time diagnostic code recorded for stroke or MI after cohort entry. The index date was defined as the calendar date of the first recorded ischaemic stroke, or MI, depending on the outcome under study. Potential cases of stroke were identified using codes referring to ischaemic, haemorrhagic and unspecified types of strokes (eg, cerebrovascular accident). An algorithm integrating examinations, procedures and treatments was developed to classify strokes not further specified as ischaemic or haemorrhagic. The computerised medical records of a random sample (10%) of the potential cases of stroke were reviewed by a neurologist (CR), blinded to exposure status, and the corresponding algorithm was modified until complete agreement between the computer and the reviewer's diagnostic classification was achieved. Following this process, strokes whose type could not be determined (30.6%) were classified as ischaemic strokes. Patients with a first code for haemorrhagic stroke were censored at the date of diagnosis and only ischaemic strokes and strokes not otherwise specified were considered as cases.
For each case, up to 10 controls were randomly selected among the cohort members in the risk sets defined by the case, after matching on sex, age (±1 year), indication for AED prescription (epilepsy, pain, psychiatric disease, other/unknown), date of cohort entry (±1 year) and duration of follow-up. When no controls were found for a case, matching criteria were relaxed for cohort entry (±2 years) and for age (±5 years). The matching on calendar time permits controlling for trends over time in the use of AEDs and changes in incidence of the outcome, whereas matching on indication of AED is essential to control for the baseline vascular risk that is likely to differ between indications. Controls being selected from the risk set defined by each case, they were necessarily alive, active in the practice and event-free when matched to their corresponding case. The date resulting in the same duration of follow-up for the case and controls defined the index date for the controls.
Definition of AED exposure
For each matched set of cases and controls, all AED prescriptions issued between cohort entry and the index date were identified. The duration of each prescription was calculated from the number of prescribed tablets combined with dosage instructions. In order to assess if the risk of ischaemic vascular event varies according to enzyme-inducing properties, AEDs were classified into the three following classes: inducing AED (carbamazepine, phenytoin, phenobarbital and primidone), inhibiting AED (sodium valproate) and non-inducing AED (all other AEDs). Patients were considered current users if the duration of their last prescription included the index date or ended within 90 days prior to the index date. We further distinguished between those who had no AED prescription in the year before the index date except on the 90 days before index date (initiators) and those whose current AED exposure started before the past 90 days (current users). This distinction ensures a relevant exposure time window as it is unlikely that an AED initiated shortly before the event would be responsible for the event based on the biological hypothesis. Moreover, it avoids a potential protopathic bias, especially when studying stroke, where starting an AED could be linked to early manifestations of cerebrovascular disease. Past users were patients whose last AED prescription ended between 365 and 90 days prior to the index date. Non-users were patients with no AED prescription in the year preceding the index date. Patients currently exposed to more than one AED class, for instance with prescriptions for both inhibiting and inducing AEDs, were classified into a separate category (multiple users). Consequently, current users, initiators, past users, multiple and non-users represented mutually exclusive exposure categories. In a subsequent analysis, patients currently exposed (based on the exposure definition above) were categorised according to duration of use during the entire follow-up from cohort entry to index date. The duration of continuous use of an AED was defined by successive prescriptions for the AED class of interest with a grace period of 90 days added to the end of a prescription to allow for refill time. If within the calculated duration a new prescription of the same AED class was issued, exposure was considered continuous, lasting until the last day of the last prescription. Duration of use was stratified into three groups according to the tertiles of the distribution of use in the controls.
In addition to the inherent adjustment by the matching variables, the following potential confounders were included in the multivariable analysis: obesity (body mass index≥30), smoking status, alcohol abuse, transient ischaemic attack, diabetes, hypertension, hyperlipidaemia, atrial fibrillation, coronary artery disease, congestive heart failure, peripheral vascular disease, chronic obstructive pulmonary disease, chronic kidney disease, cancer, blood disorders and number of medical consultations. In addition, history of stroke was included in the analysis when studying MI and history of MI when studying stroke. All variables were measured 366–730 days prior to the index date and identified with relevant diagnostic codes except hypertension, diabetes and hyperlipidaemia which were defined using diagnostic codes or related medications. Medications prescribed within 1 year before the index date, such as oral anticoagulants, aspirin and other antiplatelet drugs, non-steroidal anti-inflammatory agents, antidepressants, oral contraceptives and hormone replacement therapy, were also included as covariates. Missing data were expected for body mass index and smoking, so a separate category was created to classify this missing information; in a sensitivity analysis, we used multiple imputations for these missing values. We also performed a supplementary analysis excluding hyperlipidaemia as it is likely in the causal pathway between AED exposure and outcome. We also removed transient ischaemic attack or coronary artery disease from the potential confounders adjusted for in the multivariate model, depending on the outcome under study, as they may represent minor forms of the outcome of interest.
We calculated crude incidence rates for ischaemic stroke and MI, along with 95% CIs based on the Poisson distribution. Incidence rate ratios (RRs) of ischaemic stroke and MI associated with current use of enzyme-inducing AEDs and enzyme-inhibiting AEDs relative to non-enzyme-inducing AEDs and their corresponding 95% CIs were estimated from the ORs computed using conditional logistic regression. This approach provides unbiased estimates of the RRs that would be obtained from a Cox regression analysis, with little or no loss in precision, and is computationally more efficient, particularly when analysing time-varying exposures in large cohorts.24 In secondary analyses, we stratified cases and their matched controls according to the indication for AED to assess effect modification by the underlying reason for AED prescription. We also determined separately the risk of stroke and MI associated with current exposure to the different inducing AEDs as compared with non-inducing AEDs. Finally, we quantified the risk of stroke and MI according to the duration of AED use in an analysis restricted to current users ever exposed to one AED class only to exclude potential persistent effects of other AED classes.
We also calculated an approximation of absolute adjusted rate differences (RDs) and corresponding 95% CIs as measures of the excess risk associated with the inducing AED. Adjusted absolute RDs were calculated as R0×(RR−1), where R0 was the incidence rate of outcome among the non-inducing AED person time in the cohort (which was estimated to be 4.68/1000 persons per year for ischaemic stroke and 3.61/1000 persons per year for MI), and RR was the adjusted RR for the inducing AED versus non-inducing AED.
Sensitivity analyses were performed to assess the robustness of our results. To explore the effect of potential misclassification of current exposure on the estimated RRs, the duration of the current exposure time window was shortened to 30 days prior to the index date for all cases and controls. Thus, some participants classified as currently exposed to AEDs at the index date based on the previous definition of current use in our main analysis became past or non-users. To assess the potential adjustment for covariates in the causal pathway, we repeated the primary analysis while adjusting for covariates measured at cohort entry. Finally, to assess the potential for residual confounding, we repeated the primary analyses (enzyme-inducing vs non-enzyme-inducing AED and enzyme-inhibiting vs non-enzyme-inducing AED) using propensity scores26 (PS) and high-dimensional propensity scores27 (hdPS). The PS were estimated in controls and applied to cases using covariates measured between 1 and 2 years before the index date. To ensure overlap in the distribution of PS for each of the exposure categories, 5% of cases and controls were excluded at either end of the distribution of the PS. HdPS were estimated in controls using an algorithm described in detail elsewhere.27 Briefly, after prioritising variables according to their potential to bias the exposure–outcome association, the 300 variables most likely to cause bias were included in the hdPS model. Estimated PS (or hdPS) were then included as quintiles in the conditional logistic regression model to assess the risk of stroke and MI associated with the different AED classes. CIs were calculated using a significance level of 5%. All computations were performed using the SAS software V.9.3 (SAS Institute Inc, Cary, North Carolina, USA).
Of the 440 191 patients with at least one AED prescription, 237 050 were eligible for the stroke cohort and 241 405 for the MI cohort (figure 1). The mean age at cohort entry in the stroke cohort was 56.1 (SD=17.5) years; 39.9% were males; and indications for AED prescription were epilepsy (13.5%), pain (49.8%), psychiatric disorders (5.3%) and others/unknown (31.4%). Characteristics were virtually the same in the cohort for studying MI. The mean duration of follow-up was 4.1 years (SD=3.9) in both cohorts, during which 4535 patients were diagnosed with ischaemic stroke and 3636 with MI, yielding overall incidence rates of 4.62 (95% CI 4.49 to 4.75) and 3.63 (95% CI 3.51 to 3.75) per 1000 persons per year, respectively. The characteristics of the cases of ischaemic stroke and MI and their matched controls are presented in table 1. As expected, cases had a higher prevalence of vascular risk factors and comorbidities than controls.
In the stroke cohort, a total of 2743 (60.5%) cases and 24 971 (57.4%) controls had received at least one AED prescription in the year before the index date. Among controls, 32% of these prescriptions were inducing AEDs (carbamazepine in 65%), 18% were inhibiting AEDs and 49% were non-inducing AEDs (gabapentin in 58%). The pattern of prescriptions was similar for the MI cohort. Among current users at index date, in the stroke outcome analysis, 85% of the cases and 87% of the controls remained exposed to the same class of AED during the entire follow-up, that is, they did not switch during follow-up. For MI, the numbers were 84% and 85% for cases and controls, respectively.
Current use of inducing AED was associated with a small increased risk of ischaemic stroke (RR 1.16 (95% CI 1.02 to 1.33)) relative to current use of non-inducing AED and no increased risk of MI (RR 1.12 (95% CI 0.97 to 1.30); table 2). When stratified by indication, the increased risk of stroke was mainly driven by a 47% increased risk in the other/missing category (figure 2). The rate of stroke and MI did not vary according to the different inducing AED molecules (data not shown). When stratified by duration of use, there was no clear increased risk of ischaemic stroke with longer duration of current use of inducing AED (figure 2). However, there was an increased risk of MI with longer duration of use. Current inducing AED use of more than 24 months was associated with a 46% increased risk of MI (RR 1.46 (95% CI 1.15 to 1.85)) compared with the same duration of non-inducing AED use, corresponding to an absolute adjusted RD of 1.39/1000 (95% CI 0.33 to 2.45) persons per year.
Current use of inhibiting AED was not associated with an increased risk of ischaemic stroke (RR 1.13 (95% CI 0.96 to 1.33)) relative to current use of non-inducing AED and a non-statistically significant decreased risk of MI (RR 0.81 (95% CI 0.66 to 1.00); table 2). When stratified by duration of use, there was no clear pattern of association between the risk of ischaemic stroke or MI with longer duration of use of inhibiting AED (figure 3).
Restricting the exposure time window to 30 days did not affect the adjusted RRs of ischaemic stroke and MI estimated in the primary analysis (data not shown). Similarly, measuring covariates at cohort entry rather than 1–2 years before the index date, or controlling for confounding with PS or hdPS, did not materially change our results (figure 4). Restricting the analysis to ischaemic stroke only, that is, excluding strokes not otherwise specified, or using multiple imputations for missing data on smoking and body mass index, yielded similar estimates. Finally, adding dementia to the list of adjustment covariates and removing adjustment for transient ischaemic attack or coronary heart disease during follow-up resulted in very similar estimates.
Using a large population-based cohort of new users of AEDs, we showed no clear increased risk of ischaemic stroke with current use of inducing or inhibiting AEDs compared with non-inducing AEDs, except for a slight increased risk of ischaemic stroke associated with inducing AED, most likely due to residual confounding by indication. We found a small increased risk of MI with longer duration of inducing AED use. However, this relative increase would translate into few additional events and therefore a small clinical impact. Conversely, use of inhibiting AEDs was associated with a decreased risk of MI.
In the light of our results, the metabolic changes produced by highly inducing AEDs through the activation of the cytochrome P450 system do not lead to a clinically detectable increased risk of ischaemic stroke. The small increased risk of MI with longer duration of use but not ischaemic stroke may be related to the heterogeneity of ischaemic stroke including various aetiological subtypes. Conversely, MI is almost always caused by atherothrombosis and the association between high serum lipid levels and cardiovascular risk has been clearly established.28 ,29 The vascular effect of sodium valproate, which inhibits the cytochrome P450 pathway, is unclear. Indeed, valproate could have a proatherogenic effect through induction of insulin resistance, body weight gain, metabolic syndrome and increased oxidative stress.11 ,12 ,30–33 On the other hand, several studies suggest that valproate use is associated with lower LDL and total cholesterol serum levels.34 Moreover, studies in rodents have shown that sodium valproate attenuates the development of atherogenesis in hyperglycaemic mice via inhibition of glycogen synthase kinase-3, a serine/threonine kinase which plays a central role in many cellular processes.35 ,36 Inhibition of glycogen synthase kinase-3 is also a promising target for cardioprotection.37 The capacity of valproate to affect other intracellular enzymes and pathways such as histone deacetylase, which plays a critical role in cardiac pathology, could also lead to cardioprotective effects. Indeed, administration of valproate had a beneficial effect on left ventricular remodelling after MI in rats via inhibition of histone deacetylase activity.38 Valproate may also exert a protective effect on large vessel ischaemic stroke via inhibition of histone deacetylase. Indeed, a genetic association between a variant of HDAC9 (encoding histone deacetylase 9) and this stroke subtype was recently identified.39 However, we did not have information on ischaemic stroke subtype and this association may have been masked in our analyses.
Only one previous study addressed the risk of ischaemic arterial events according to AED enzymatic properties using a US Medicare database and found no clear increased risk. However, AEDs were classified into two categories only, with the enzyme-inhibiting AED (valproate) being classified with the non-inducing AEDs, and the mean duration of follow-up was only 4 months.40 In a Danish cohort study of patients with epilepsy, those treated with AEDs were at increased risk of vascular events as compared with non-treated patients except for MI.41 In a secondary analysis, the authors examined the risk associated with few individual AEDs. Compared with carbamazepine, sodium valproate was associated with a decreased risk of stroke and MI, whereas oxcarbazepine and phenobarbital were associated with an increased risk of death, and oxcarbazepine was associated with an increased risk of stroke. Several methodological differences with our work are noticeable such as the inclusion of incident and prevalent AED users at cohort entry, the assessment of few individual AEDs rather than AEDs according to enzyme properties, and the lack of adjustment for some potential confounding factors such as smoking or obesity. Although no consistent trends were found with AED individually, their finding of an attenuated risk of MI with valproate use is compatible with our results. However, we only found a trend towards a decreased risk for MI, which may be due in particular to the tighter adjustment for confounding in our study. One recent study using the CRPD database examined the risk of ischaemic stroke and MI associated with sodium valproate in a cohort of patients with epilepsy who were incident and prevalent users of one of the four following AEDs: sodium valproate, carbamazepine, phenobarbital and phenytoin.42 The authors reported a decreased risk of MI with sodium valproate compared with the use of other AEDs and a decreased risk of stroke with longer duration of use. However, the same trends were observed with the three other AEDs regarding the risk of stroke, and the results are most likely explained by several methodological issues as acknowledged by the authors. To the best of our knowledge, no previous study has assessed the risk of ischaemic events separately for the three AED categories while taking the indication for AED into account.
Our observational study was conducted in a population-based setting, using the CPRD. The organisation of the UK healthcare system where over 99% of people are registered with a general practitioner ensures representativeness of the study population. The proposed design, consisting of a nested case–control approach within a well-defined cohort, virtually rules out the potential for selection bias. Also, the two outcomes, stroke and MI, have been extensively studied and well validated in the CPRD. Finally, since drug prescriptions are automatically recorded in the computerised file when issued by the physician, the definition of exposure is exempt of recall bias. However, it was based solely on prescriptions issued and not on prescriptions actually filled or taken by the patient, which could result in misclassification of exposure that is likely to be non-differential between cases and controls, therefore potentially biasing the results towards the null. The effect of potential misclassification was examined in sensitivity analyses which produced similar results. One of the potential threats to validity may be the presence of residual confounding due to differential drug prescription according to disease type or disease severity, that is, confounding by indication. By design, this bias was minimised by estimating the risk within groups of patients with the same indication and by choosing a treated comparator group. Indeed, treated patients are fundamentally different from non-treated patients in terms of disease severity, and using the latter group as a comparator would preclude disentangling the respective role of treatment and severity of the disease if the severity itself is associated with vascular events. Moreover, since the hypothesis is that the vascular risk varies with enzyme-inducing properties of AEDs, rather than with AEDs in general, a treated comparison group is more relevant. However, some residual confounding may still be present in the stratified analyses after matching by indication, but to various degrees within each indication. Still, our results were consistent within each matching category, except for the category with missing indication (for the stroke outcome) where we may expect a greater amount of residual confounding. This suggests that matching by indication was an efficient strategy to minimise confounding and that simply adjusting for numerous covariates in a multivariate model as carried out in some previous studies would not be sufficient. Of note, indication for a particular medication is not readily available in administrative databases either. Therefore, previous studies restricted their cohort to patients with epilepsy or simply adjusted for several comorbidities in the analyses. In theory, restriction is a good strategy, but given the incomplete information on indication in databases, this may lead to other biases such as selection bias. We were able to attribute AED prescription to a particular indication in up to 70% of our cohort, which permitted matching and possibly tighter adjustment. Moreover, we could assess that the risk was consistent in various indications, which was not possible in other studies. To further decrease the potential for residual confounding, all models were adjusted for several important potential confounders, including lifestyle variables such as smoking and alcohol abuse, which are usually not recorded in administrative databases and were not adjusted for in previous studies. Finally, we undertook several additional analyses using alternative methods to control for confounding such as PS and hdPS. We showed that, although the estimates varied slightly, particularly when adjusting using hdPS, the effect size remained consistent with a minor clinical impact.
In summary, we did not find a clear and consistent increased risk of ischaemic stroke with the use of inducing or inhibiting AEDs when compared with non-inducing AEDs. We found an increased risk of MI with longer duration of use of inducing AEDs, which is most likely of small clinical importance. Conversely, use of inhibiting AEDs was associated with a decreased risk of MI, consistent with the underlying biological hypothesis.
Contributors CR designed the study. SD and J-FB critically reviewed the study protocol. SD analysed the data. OS participated in the statistical analysis. CR wrote the manuscript. All authors revised the manuscript. CR is the guarantor of the manuscript. All authors had full access to all of the data (including statistical reports and tables) in the study and can take responsibility for the integrity of the data and the accuracy of the data analysis.
Funding This study was funded by an operating grant (DRB-126643) from the Drug Safety and Effectiveness Network of the Canadian Institutes of Health Research (CIHR). CR is the recipient of a Chercheur-Boursier Award from the Fonds de la recherche en santé du Québec (FRSQ). KBF holds a New Investigator Award from the CIHR.
Competing interests CR, KBF and J-FB report grants from the Canadian Institutes of Health Research (CIHR) during the conduct of the study.
Ethics approval The study protocol was approved by the Independent Scientific Advisory Committee of the CPRD (No. 14_181R) and the Research Ethics Committee of the Jewish General Hospital, Montreal, Canada; this protocol was made available to the journal reviewers.
Provenance and peer review Not commissioned; externally peer reviewed.
Data sharing statement No additional data are available.