Article Text

Download PDFPDF

Original research
Persistence and safety of anti-TNF biosimilars versus originators in immune-mediated inflammatory diseases: an observational study on the French National Health Data System
  1. Hugo Jourdain1,
  2. Léa Hoisnard2,3,4,
  3. Emilie Sbidian1,2,3,4,5 and
  4. Mahmoud Zureik1,6
  1. 1EPI-PHARE, French National Agency for Medicines and Health Products Safety (ANSM) and French National Health Insurance (CNAM), Saint-Denis, France
  2. 2Fédération Hospitalo-Universitaire TRUE InnovaTive theRapy for immUne disordErs, Assistance Publique-Hôpitaux de Paris (AP-HP), Henri Mondor Hospital, Créteil, France
  3. 3Centre d’Investigation Clinique 1430, INSERM, Créteil, France
  4. 4EpiDermE Epidemiology in Dermatology and Evaluation of Therapeutics, EA7379, Paris Est Créteil University UPEC, Créteil, France
  5. 5Department of Dermatology, Assistance Publique-Hôpitaux de Paris (AP-HP), Henri Mondor Hospital, Créteil, France
  6. 6Anti-Infective Evasion and Pharmacoepidemiology, CESP, University Paris-Saclay - UVSQ, Montigny le Bretonneux, France
  1. Correspondence to Hugo Jourdain; hugo.jourdain{at}assurance-maladie.fr

Abstract

Objectives Biosimilar-originator equivalence has been demonstrated in phase 3 trials in a few indications of infliximab, etanercept and adalimumab. The objective of our study was to compare the persistence and safety of biosimilars versus originators in all the licensed indications of these molecules.

Methods We used data from the French National Health Data System (SNDS), covering 99% of the French population, to identify infliximab, etanercept and adalimumab initiators from biosimilar launch (January 2015, May 2016 and October 2018, respectively) to 30 June 2021. Patients were then followed for 1 year. Treatment persistence (duration without treatment discontinuation or modification) and safety (including severe infections, all-cause hospitalisation and death) were compared between originator and biosimilar users by Cox regressions weighting the populations on the inverse probability of treatment. Analyses were performed by molecule, by disease and by biosimilar product.

Results From January 2015 to June 2021, 86 776 patients were included in the study: 22 670, 24 442 and 39 664 patients had initiated infliximab, etanercept and adalimumab, respectively; 49 752 (53%) were biosimilar initiators. We did not find any risk of discontinuation (HRs were below or around 1, here all pathologies and products together: infliximab 0.88 (0.80–0.97), etanercept 0.85 (0.81–0.90) and adalimumab 0.96 (0.91–1.00)) or safety event (infection: infliximab 0.97 (0.78–1.21), etanercept 1.04 (0.81–1.33) and adalimumab 0.98 (0.83–1.16); hospitalisation: infliximab 1.08 (0.96–1.23), etanercept 0.99 (0.87–1.11) and adalimumab 0.91 (0.83–0.99)) associated with biosimilar versus originator use.

Conclusions Our study shows reassuring results regarding the persistence and safety of biosimilar tumour necrosis factor-alpha inhibitors compared with originators in all licensed indications.

  • Biological Therapy
  • Biosimilar Pharmaceuticals
  • Tumor Necrosis Factor Inhibitors
  • Arthritis, Rheumatoid
  • Spondylitis, Ankylosing

Data availability statement

No data are available. In accordance with data protection legislation and the French regulation, the authors cannot publicly release the data from the French National Health Data System.

http://creativecommons.org/licenses/by-nc/4.0/

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, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

WHAT IS ALREADY KNOWN ON THIS TOPIC

  • Biosimilar-originator equivalence was proven in phase 3 randomised clinical trials in few indications of infliximab, etanercept and adalimumab, and then extrapolated to all the other indications of these molecules.

  • We thus aimed to address the persistence and safety of biosimilar products compared with originator products for each brand name product and each licensed indication of each molecule.

WHAT THIS STUDY ADDS

  • This is the largest cohort of biosimilar or originator infliximab, etanercept and adalimumab initiators from 2015 to 2021.

  • This is the first observational study to address persistence and safety endpoints by comparing marketed biosimilar products and originator products of the tumour necrosis factor (TNF)-alpha inhibitor family in a systematic approach for each licensed indication and each branded product.

  • We did not show any difference in persistence or safety between biosimilar products and their originators in any licensed indication.

HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY

  • This study provides reassuring results regarding the use of anti-TNF alpha biosimilar products in real life.

  • Marketed biosimilars of TNF-alpha inhibitors can thus be used more broadly in order to limit the financial impact of these molecules on health systems.

Introduction

Tumour necrosis factor-alpha (TNF-alpha) inhibitors are a class of biotherapies that have significantly improved the therapeutic care of immune-mediated inflammatory diseases. Infliximab (IFX), etanercept (ETA) and adalimumab (ADA), the first three marketed molecules of this class, were used and reimbursed respectively in 2000, 2003 and 2005 in France and indicated, according to the molecule, in rheumatoid arthritis, ankylosing spondylitis and psoriatic arthritis in rheumatology, psoriasis and hidradenitis suppurativa in dermatology, Crohn’s disease and ulcerative colitis in gastroenterology and uveitis in ophthalmology.1 However, these treatments were expensive, leading to an issue of either access for patients or healthcare system sustainability due to the large number of patients treated. TNF-alpha inhibitors patent expiry and biosimilars progressive market approvals partly answered these issues by lowering the cost of treatment.2 The active substance of biosimilars is highly similar to the originator’s. Biosimilar products are approved provided that the data demonstrating biosimilarity to the originator product include physiochemical, non-clinical (mostly in vitro studies), pharmacokinetics and pharmacodynamics studies and generally one comparative phase 3 clinical trial to demonstrate their bioequivalence.3 The bioequivalence of the biosimilar products is then extrapolated to the other indications of the originator product. Several observational studies assessed the effectiveness of biosimilar products in these remaining indications and showed that there was no difference between originator and biosimilar products. However, they were often small,4 5 non-comparative6–8 or pathology-specific,9–12 as we already did for inflammatory bowel diseases using the same database from 2015 to 2017.13 14

In Europe, four IFX, four ETA and 12 ADA biosimilars have been approved from 2015, 2016 and 2018, respectively,15 based on bioequivalence data demonstrated in 1–3 indications. In France, the use of biosimilar products has experienced a significant increase since their market approval.16

The aim of the present study was to assess and compare the persistence under biosimilar versus originator treatment among initiators for each biosimilar product, taken separately and overall, for each pathology for which the originator products of IFX, ETA and ADA were indicated. In addition, we also compared the safety profiles of biosimilar and originator users in major adverse events.

Methods

Study design

We built cohorts of IFX, ETA and ADA initiators, exposed to the originator molecule or to a biosimilar molecule, from the date of reimbursement of the first biosimilar product for each molecule (2015, 2016 and 2018, respectively) and followed the patients for 1 year.

Data source

We used data from the French National Health Data System (SNDS). The SNDS covers almost the totality (>99%) of the French population—68 million residents. Each person is identified by a unique and anonymous number. The SNDS records comprehensive outpatient (procedures and pharmacy deliveries of reimbursed drugs) and inpatient (pharmacy deliveries of expensive drugs, procedures performed during hospital stays and discharge diagnoses coded according to the International Statistical Classification of Diseases and Related Health Problems, Tenth Revision (ICD-10)) reimbursement information since 2006. The SNDS also contains sociodemographic information on sex, age, place of residence and vital status, among others. Patients’ status for 100% reimbursement of care related to a severe and costly long-term disease (LTD) is recorded and LTD diagnosis is coded according to the ICD-10. The SNDS has been extensively used to conduct pharmacoepidemiological studies, especially on the use, safety and effectiveness of health products.17–21

Study periods and populations

Inclusion periods differed for each originator molecule, according to the date of reimbursement in France of the first biosimilar for each of the three studied TNF-alpha inhibitor originators: IFX (ATC L04AB02), ETA (ATC L04AB01) and ADA (ATC L04AB04). Inclusions started on 27 January 2015 for IFX users, on 10 May 2016 for ETA users and on 9 October 2018 for ADA users. Every patient initiating treatment, that is, having a first delivery with no delivery of IFX, ETA or ADA within the year before initiation, between the inclusion start date and 30 June 2021, was included in the study.

Patient characteristics

Sociodemographic characteristics (sex, age, affiliation to complementary universal health insurance—CSS, deprivation index and region of residence) were collected at the index date. Comorbidities were assessed within 5 years before the index date. The history of visits to a specialist was assessed within 2 years before the index date. The history of dispensing of corticosteroids, non-steroidal anti-inflammatory drugs, non-biological systemic drugs on the one side and biological and targeted drugs on the other side was assessed within the year before the index date (ATC codes and grouping specified in online supplemental table 1). The pathology treated by the TNF-alpha inhibitor was defined at the index date. ATC codes, ICD-10 codes and the procedure used for pathology definition were presented in a previous study.16 Patients under 18 or treated with a molecule not indicated for their pathology were excluded.

Exposure, outcome, censoring and safety

The exposure of interest was, within a molecule, originator use versus biosimilar use. Specific exposure to each marketed biosimilar product (CT-P13 and SB2 for IFX; SB4 and GP2015 for ETA; ABP501, MSB11022, FKB327, GP2017 and SB5 for ADA; the other biosimilar products were excluded due to a lack of patients) was assessed as a subgroup analysis. Exposure was defined at the index date and reassessed at each subsequent delivery or non-delivery. The main outcome was persistence in treatment, defined as a period of continued treatment with the same product. Discontinuation, specified as a treatment gap of more than 60 days after the theoretical coverage period of each molecule (56 days for IFX, 28 days for ETA and ADA), and treatment modification, defined as the delivery of a different biologic molecule, were considered outcome events. Censoring events included intra-molecule switching (at least two consecutive deliveries of a different product from the same molecule), death and the end of the 1-year follow-up period. Safety events included severe infection,22 major adverse cardiovascular events,23 immunological disorders, cancer, anaphylaxis,24 identified through ICD-10-coded hospitalisation discharges (specified in online supplemental table 2), death and all-cause hospitalisations for at least one night, except pathology-related hospitalisations.

Statistical analysis

Patients’ characteristics at inclusion and during the follow-up (follow-up duration and outcome, censoring and safety event rates) were described. To account for the indication bias on exposure groups, a stabilised propensity score of treatment was computed using a logistic regression modelling the probability of initiating a biosimilar compared with the originator adjusted for all the covariates measured at inclusion (year of inclusion, age in categories, sex, affiliation to complementary universal health insurance, deprivation index, region of residence, main comorbidities, number of consultations and hospitalisations linked to the pathology in categories, drug history). Covariate balance was assessed graphically, plotting the standardised mean difference in covariates and the distribution of the propensity score before and after weighting. Inverse probability of treatment weighting (IPTW), truncating the first and 99th percentiles of weights to keep non-extreme weights, was used in Cox proportional hazards regressions comparing the risk of non-persistence for biosimilar users versus originator users. As ponderation was not sufficient to balance the covariates completely (as shown in online supplemental figures 1–3), we further adjusted the model for all the covariates included in the propensity score model. A univariate unweighted Cox proportional hazards regression was also computed. All the analyses were repeated for each pathology and each biosimilar treatment group (all together and each biosimilar taken separately). The proportionality of the hazards hypothesis was confirmed using the test of Schoenfeld residuals. The 0.05 p value threshold was corrected by dividing it by the number of pathologies tested for each molecule (six for IFX, four for ETA and eight for ADA) to follow the Bonferroni correction,25 and two-sided tests were carried out. The p values and p value thresholds were reported in Supplementary tables. The whole method was repeated for three of the safety events, namely severe infections, all-cause hospitalisation and death. Lastly, we only described (effectives and proportions) the other events as they were less frequent and we lacked the power to carry out multivariable analyses.

All extractions from the SNDS were carried out with SAS Enterprise Guide software V.7.15; analyses were done with R26 V.3.5.2, using multiple packages including dplyr,27 ggplot2,28 survival29 and cobalt.30

Sensitivity analyses

Several sensitivity analyses were carried out. First, the biosimilar to originator switch was considered as a treatment discontinuation outcome, as we made the hypothesis that this switch would correspond to a medical switch for lack of efficacy of the biosimilar. Second, the intra-molecule switch was defined as a three-consecutive delivery change of product (and not two), not to misclassify a punctual switch due to stock issues. We also considered the intra-molecule switch at the first delivery to be comprehensive. In addition, we modified the definition of discontinuation, taking a 30-day or 90-day (vs 60-day gap in the main analysis) gap in addition to the theoretical covering period of the molecules. Fourth, we restricted the analyses to at least 6-month persistent patients, to analyse only treatment responders. Fifth, we modified the inclusion criteria, carrying out an additional analysis by changing the definition of TNF initiation to ‘no delivery of any of the three TNF-alpha within the past 5 years’ to prevent including non-initiating patients. We also restricted the analyses to calendar years with at least 10% of patients included in each group to have a minimal balance for this variable, and we repeated the analyses for each year of inclusion as a subgroup analysis. Additionally, the follow-up was extended to 2 years and also restricted to the pre-COVID-19 time period to exclude pandemic-specific discontinuations and the delay in care for new biologics users in the postlockdown period. Finally, the relative risk (RR) of non-persistence at 1 year was computed using the IPTW Kaplan-Meier survival estimate at 1 year; the CI for this estimate was assessed using the percentiles of the RR estimates on 1000 bootstraps.

Results

Patient characteristics

A total of 86 776 patients were included, of whom 22 670 were IFX users, 24 442 were ETA users and 39 664 were ADA users. Among IFX users, 5292 (23%) initiated treatment with the originator product, 14 678 (65%) with CT-P13 and 2700 (12%) with SB2; 2042 (9%) had rheumatoid arthritis, 4059 (18%) ankylosing spondylitis, 611 (3%) psoriatic arthritis, 9531 (42%) Crohn’s disease, 5054 (22%) ulcerative colitis and 1373 (6%) psoriasis. Among ETA users, 13 354 (55%) initiated treatment with the originator product, 8917 (36%) with SB4 and 2171 (9%) with GP2015; 11 265 (46%) had rheumatoid arthritis, 9702 (40%) ankylosing spondylitis, 1556 (6%) psoriatic arthritis and 1919 (8%) psoriasis. Among ADA users, 18 378 (46%) initiated treatment with the originator product, 10 048 (25%) with ABP501, 4189 (11%) with FKB327, 3939 (10%) with SB5, 1870 (5%) with MSB11022 and 1240 (3%) with GP2017; 5111 (13%) had rheumatoid arthritis, 11 091 (28%) ankylosing spondylitis, 1857 (5%) psoriatic arthritis, 11 047 (28%) Crohn’s disease, 5991 (15%) ulcerative colitis, 3651 (9%) psoriasis, 231 (1%) hidradenitis suppurative and 685 (1%) uveitis.

Patients’ characteristics at inclusion were similar between biosimilar and originator product users within each molecule, as shown in table 1 and online supplemental tables 3–5. The only significant difference came from the year of inclusion, as biosimilar product uptake was progressive throughout the years.

Table 1

Main patients’ characteristics at baseline

Persistence

All pathologies taken together, patients treated with a biosimilar of IFX were followed on average for 280 days (SD 112 days), slightly more than originator IFX users (269 (120) days). Both groups had similar discontinuation or molecule switch rates at 1 year (38.9% (n=6761) in the biosimilar group and 38.3% (n=2026) in the originator group).

In ETA users, patient follow-ups in the biosimilar versus originator groups were close (264 (118) days vs 253 (120) days). The biosimilar group experienced fewer events than the originator group (45.7% (n=5064) compared with 52.2% (n=6965)).

Patients treated with a biosimilar of ADA were followed on average for 267 days (SD 117 days), slightly less than originator ADA users (273 (114) days). Event rates at 1 year in both groups were similar (42.8% (n=9104) in the biosimilar group and 43.1% (n=7916) in the originator group).

IPTW aHR (presented in tables 2–4 and online supplemental tables 9–17) of event were close to 1 and CIs covered 1, except for several subgroups: in the IFX cohort, biosimilar use showed a better persistence in Crohn’s disease (aHR 0.83 (95% CI 0.71 to 0.96), driven by the CTP13); in the ETA cohort, biosimilar use showed a better persistence in rheumatology (rheumatoid arthritis: aHR 0.85 (95% CI 0.78 to 0.93); ankylosing spondylitis: aHR 0.88 (95% CI 0.80 to 0.96) and psoriasis: aHR 0.80 (95% CI 0.65 to 1.00); in the ADA cohort, biosimilar use showed a better persistence in ankylosing spondylitis (aHR 0.89 (95% CI 0.81 to 0.97). In addition, each biosimilar product subgroup had better or similar persistence than the originator groups (online supplemental tables 9–17). All the sensitivity analyses (shown in online supplemental figure 4) were consistent with the results of the main analysis.

Table 2

Risk of non-persistence of biosimilar use versus originator use in infliximab-initiating patients, by pathology

Table 3

Risk of non-persistence of biosimilar use versus originator use in etanercept-initiating patients, by pathology

Table 4

Risk of non-persistence of biosimilar use versus originator use in adalimumab-initiating patients, by pathology

Safety

Across molecules, adverse event rates were mostly similar. The most common adverse events were severe infection, all-cause hospitalisation and death, with 6.0% (n=1370), 2.3% (n=566) and 2.5% (n=982) of IFX, ETA and ADA patients having experienced hospitalisation for severe infection, respectively; 19.5% (n=4424), 10.2% (n=2495) and 8.8% (n=3505) of IFX, ETA and ADA patients having experienced all-cause hospitalisation, respectively; and 0.8% (n=173), 0.3% (n=71) and 0.2% (n=74) of IFX, ETA and ADA patients having died before the end of the follow-up, respectively. Weighted Cox proportional hazard models showed no difference between biosimilar and originator users regarding these adverse events (figure 1, online supplemental tables 18–20). The description of other safety outcomes for each group is available in online supplemental tables 6–8. All the sensitivity analyses (shown in online supplemental figure 5) were consistent with the results of the main analysis.

Figure 1

Estimation of the association between product type exposure (biosimilar vs originator, with originator as a reference) and major adverse events in antitumour necrosis factor-alpha-initiating patients, by molecule and adverse event type (crude HRs are presented in grey, whereas IPTW HRs, adjusted on all the covariates, are presented in black). HR values are presented in point estimates (CI) for the IPTW models only. aHR adjusted HR; IPTW, inverse probability of treatment weighting.

Discussion

The marketed biosimilar TNF-alpha inhibitors had been approved based on specific phase 3 randomised controlled trials in 1–3 indications, and their equivalence was then extrapolated in all the other indications of the molecule. In this systematic observational study based on the SNDS gathering all French patients, we showed no difference in terms of persistence or safety between biosimilar TNF-alpha inhibitors and originator products, within each molecule analysed, namely, IFX, ETA and ADA, in all their respective licensed indications.

We have consistent results with the clinical trials that led to the market authorisation of the biosimilar products (as shown in table 5) but also with observational studies in other indications demonstrating either no difference between biosimilars and originators9 11 12 14 or a better persistence in biosimilar-treated patients (IFX CT-P13 had already been shown to be slightly beneficial in Crohn’s disease in a study carried out on the same database,13 and ETA biosimilars had better retention in a French recent study).31 For the first time, our study assessed, in a large unselected population, persistence and safety in all the indications and all the anti-TNF-alpha products altogether.

Table 5

Results overview filling the randomised controlled trial gaps with subgroup analyses by biosimilar product in all the indicated pathologies

While clinical trials mainly focus on efficacy endpoints, our study intends to assess biosimilar versus originator equivalence in its globality: persistence first, as a proxy of real-world effectiveness, but also safety, through a complete description and advanced statistical comparative analyses.

Biosimilar products share the same amino acid sequence as originator products, but their tri-dimensional structure, exact formulation and excipients may vary from one product to another. Manufacturing processes may alter the composition of biosimilar or originator products over time (referred to as molecule drift),32 potentially triggering a divergence between products from the clinical trial lots to the actual packs used in the real world. In addition, the extrapolation paradigm can be questioned by practitioners and patients, even irrationally, and educational efforts had to be made to counter prescription inertia33 and the nocebo effect.34 The aim of our study was then to break down the remaining barriers to the adoption of biosimilar products.

Our study has several strengths. First, as it was based on a comprehensive database, we were able to include almost every patient living in France and having initiated one of the three studied anti-TNF-alpha, which represented almost 90 000 patients, and to follow them through their contact with health services. It allowed us to use a systematic approach to compare biosimilar and originator products in a real-life setting based on the events of interest in all of the licensed indications of the three TNF-alpha inhibitors at stake and to carry out subgroup analyses on each biosimilar product. Second, we tested the robustness of our results with many sensitivity analyses. The selection of the population could have been a source of bias, in particular as the definition of initiation (no delivery of any of the three TNF-alpha within the past year) allowed multiple initiations for different products. The sensitivity analysis we carried out with a modified definition (no delivery of any of the three TNF-alpha within the past 5 years) did not differ from the main analysis. We also varied the outcome definition, the follow-up duration and the statistical methods. None of these modifications changed the main results of our study. Third, the indication bias, often argued in observational studies, should be limited in our case, as the originator and biosimilar products are vastly considered interchangeable (as demonstrated by the well-balanced groups at inclusion and the fast adoption of the biosimilar products).16

Our study has nonetheless several limitations. First, we do not report long-term persistence and safety data, censoring patient follow-ups for 1 year only. We adopted this approach to be comparable with randomised clinical trials on biosimilars, which generally report equivalence data at 52 weeks, and to have similar follow-ups between both exposure groups. As biosimilar uptake was progressive, the biosimilar group would have had a shorter follow-up without our arbitrary right censoring. Second, we included a sufficient number of patients to get enough events for most of our analyses. However, for several pathology subgroups, we did not have enough sample size to show results at the drug level. Moreover, as we took an interest in adverse events in the therapeutic sequence with limited follow-up, safety was reported only for a few event types and results had to be pulled at the molecule level. Third, we were somewhat limited by the nature of the database we used, as it does not capture most patients and disease characteristics, such as baseline disease activity or autoantibody status. We were not able to use the same clinical outcomes as in clinical trials, such as the American College of Rheumatology criteria, the Psoriasis Area Severity Index and the Crohn’s Disease Activity Index. We had thus to identify many covariates (such as pathology) and outcomes (medical discontinuation) with indirect information or proxies, which can be questionable. Specifically regarding intra-molecule switches, we could not distinguish between medical and non-medical switches and thus chose to consider switches as censoring events, which could be an explanation for the beneficial effects found in some biosimilar groups. Also, other factors, such as patient and physician preferences, which are not recorded in the database, could influence drug prescription and continuation. As a consequence, in spite of the fact that our algorithms and methodological choices were validated by a large number of previous studies based on the SNDS and specific to the diseases we studied,13 14 22 35 36 residual bias by indication and misclassification error cannot be precluded. Finally, our systematic approach can be disputable, as we used the same outcome criteria and adjustment covariates for each of our analyses. We think this is in fact one of the interests of the study, as our methodology is highly reproducible in time (for data updates with more patients and more follow-up) and space (in other countries and data sources).

To conclude, our study shows reassuring results regarding the persistence and safety of biosimilars compared with the originator anti-TNF-alpha product in all licensed indications of these molecules in a real-world use context. As originator and biosimilar products seem interchangeable, biosimilar products should be used more broadly.

Data availability statement

No data are available. In accordance with data protection legislation and the French regulation, the authors cannot publicly release the data from the French National Health Data System.

Ethics statements

Patient consent for publication

Ethics approval

EPI-PHARE has permanent regulatory access to the data via its constitutive bodies ANSM and CNAM, thus this present work did not require the approval from the French Data Protection Authority (CNIL). The study was registered on the study register of EPI-PHARE under the number T-2022-11-433.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Contributors Literature search: all; study design: all; data extraction and collection: HJ; statistical analysis: HJ; interpretation of the results: all; drafting of the manuscript: HJ; revision of the manuscript: all; supervision: ES and MZ and guarantor: MZ. HJ and MZ confirm that they had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.

  • Funding The authors are employees of the French National Health Insurance (CNAM), the French National Agency for Medicines and Health Products Safety (ANSM) and the Assistance Publique-Hôpitaux de Paris (AP-HP) and received no funding for this study. The funding source did not intervene at any step of the study.

  • Competing interests None declared.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.