key: cord-1020500-y8of8zo8 authors: Wilcox, Christopher R; Islam, Nazrul; Dambha-Miller, Hajira title: Association between influenza vaccination and hospitalisation or all-cause mortality in people with COVID-19: a retrospective cohort study date: 2021-03-04 journal: BMJ Open Respir Res DOI: 10.1136/bmjresp-2020-000857 sha: ad53274b2651e2881c1766e2d0163637d57620af doc_id: 1020500 cord_uid: y8of8zo8 INTRODUCTION: Recent evidence suggests that influenza vaccination may offer protection against COVID-19 severity. Our aim was to quantify the association between influenza vaccination status and risk of hospitalisation or all-cause mortality in people diagnosed with COVID-19. METHODS: A retrospective cohort study using routinely collected health records from patients registered to a General Practitioner (GP) practice in South West England within the Electronic Care and Health Information Analytics database. The cohort included 6921 people with COVID-19 during the first wave of the pandemic (1 January–31 July 2020). Data on influenza vaccination, hospitalisation and all-cause mortality were ascertained through linked clinical and demographic records. We applied propensity score methods (stabilised inverse probability of treatment weight) to quantify the association between influenza vaccination status and COVID-19 outcomes (hospitalisation or all-cause mortality). RESULTS: 2613 (38%) participants received an influenza vaccination between 1 January 2019 and COVID-19 diagnosis. Receipt of influenza vaccination was associated with a significantly lower odds of hospitalisation or all-cause mortality (adjusted OR: 0.85, 95% CI 0.75 to 0.97, p=0.02), and 24% reduced odds of all-cause mortality (adjusted OR: 0.76, 95% CI 0.64 to 0.90). DISCUSSION: Influenza vaccination was associated with a 15%–24% lower odds of severe COVID-19 outcomes. The current UK influenza vaccination programme needs urgent expansion as an integral component of the ongoing response plans to the COVID-19 pandemic. ABSTRACT Introduction Recent evidence suggests that influenza vaccination may offer protection against COVID-19 severity. Our aim was to quantify the association between influenza vaccination status and risk of hospitalisation or all-cause mortality in people diagnosed with COVID- 19 . Methods A retrospective cohort study using routinely collected health records from patients registered to a General Practitioner (GP) practice in South West England within the Electronic Care and Health Information Analytics database. The cohort included 6921 people with COVID-19 during the first wave of the pandemic (1 January-31 July 2020). Data on influenza vaccination, hospitalisation and all-cause mortality were ascertained through linked clinical and demographic records. We applied propensity score methods (stabilised inverse probability of treatment weight) to quantify the association between influenza vaccination status and COVID-19 outcomes (hospitalisation or all-cause mortality). Results 2613 (38%) participants received an influenza vaccination between 1 January 2019 and COVID-19 diagnosis. Receipt of influenza vaccination was associated with a significantly lower odds of hospitalisation or allcause mortality (adjusted OR: 0.85, 95% CI 0.75 to 0.97, p=0.02), and 24% reduced odds of all-cause mortality (adjusted OR: 0.76, 95% CI 0.64 to 0.90). Discussion Influenza vaccination was associated with a 15%-24% lower odds of severe COVID-19 outcomes. The current UK influenza vaccination programme needs urgent expansion as an integral component of the ongoing response plans to the COVID-19 pandemic. BACKGROUND COVID-19 caused by the SARS-CoV-2 is a rapidly evolving global public health emergency. Since the outbreak began, over 87 million cases have been reported globally, with over 1.8 million attributed deaths (as of 7 January 2020). 1 Of those who develop symptoms, approximately 20% will have severe disease, including pneumonia and other lifethreatening complications. 2 Older adults and those with underlying heath conditions are particularly susceptible to severe outcomes; however, recent evidence suggests that even mild COVID-19 may result in prolonged symptoms and/or long-term adverse health outcomes. 3 The pandemic has carried a substantial economic burden globally, 4 and several countries have experienced unprecedented rise in demand on healthcare systems far exceeding capacity to respond. 5 There remain limited effective interventions to treat COVID-19. A number of vaccines against SARS-CoV-2 have recently been approved for use following promising clinical trials 6 7 ; however, ongoing challenges remain regarding the manufacturing, distribution and accessibility of vaccines internationally, 8 and there is an ongoing drive to identify whether any available treatments could be repurposed to prevent COVID-19 acquisition and/or disease severity. One area of interest is in the role of existing vaccines. Vaccination against specific pathogens with both live-attenuated and inactivated vaccines can induce adaptive and innate immunity in the form of neutralising antibodies and/or antigen-specific cellular immunity. Growing evidence demonstrates that vaccination Key messages ► We sought to determine the association between influenza vaccination status and risk of hospitalisation or all-cause mortality in people diagnosed with COVID-19. ► In people with COVID-19 during the first wave of the pandemic, we found that influenza vaccination was associated with a 15%-24% lower odds of hospitalisation or all-cause mortality. ► Influenza vaccination may have the potential to offer an additional means of mitigating serious adverse complications of COVID-19, and we therefore propose that a rapid expansion of the influenza vaccination programme is urgently needed as an integral component of the ongoing response to the COVID-19 pandemic. Open access may also have a role in wider protection against unrelated pathogens. 9 10 The underlying mechanism for this phenomenon remains poorly understood and may occur by induction of innate immune responses, so-called trained innate immunity. 9 11-14 It is hypothesised that influenza vaccination may potentially have a direct protective effect against COVID-19. 6 9 15-17 This is in addition to its potential indirect benefits, such as reducing the burden of additional respiratory illness among both individual patients and healthcare services. [18] [19] [20] This hypothesis has been supported by a limited number of observational studies, including five retrospective cohort studies 15 17 21-23 (three of which were limited by the use of aggregated data), demonstrating a negative association between influenza vaccination and COVID-19 outcomes including COVID-19 acquisition, severity and/ or mortality. One prospective study using individual-level data demonstrated a negative association between influenza vaccination and COVID-19 acquisition but did not report on subsequent COVID-19 outcomes. 24 Further validation is required from large UK-based cohorts using individual-level data, exploring important outcomes including all-cause mortality and hospitalisation. This may help us better understand the impact of influenza vaccination in COVID-19 disease and inform the current influenza vaccination programme during the second wave of the pandemic. The aim of this study was to quantify the association between influenza vaccination status and COVID-19 outcomes including hospitalisation and all-cause mortality. Retrospective cohort analysis. Data were obtained from the Care and Health Information Analytics (CHIA) database. This is an anonymous National Health Service (NHS) database that includes individual-level data from 1.5 million patients across 160 GP practices in Southern England. Data have been extracted from routine general practice records linked to local secondary care hospitals, including demographic data, primary/secondary care consultations, diagnoses, medications and outcomes of primary/secondary care investigations. A cohort of individuals within the CHIA database diagnosed with COVID-19 during the first wave of the COVID-19 pandemic (1 January-31 July 2020) were identified. We included both confirmed and suspected cases. Confirmed cases had a positive RT-PCR assay for SARS-CoV-2 on nasal/pharyngeal swab. Clinically suspected cases were defined in accordance with Public Health England (PHE) guidance at the time of the first wave of the pandemic in the UK, when testing was not readily available. Further detail on the COVID-19 case definition has been reported in other databases, 25 and our rationale for case definition can also be found in a consensus statement by the COVID-19 Primary Care Database Consortium. 26 Recent evidence demonstrates that clinical suspicion of COVID-19 among GPs closely matches with confirmed cases and subsequent clinical outcomes. 27 Sociodemographic and clinical variables Data for the identified cohort were extracted from the CHIA database including age, sex, ethnicity, Index of Multiple Deprivation (IMD), electronic frailty index, comorbidities, smoking status, body mass index (BMI), medications and influenza vaccination status. Receipt of recent influenza vaccination was defined as an electronic recording of vaccination in the CHIA database between 1 January 2019 and the diagnosis of COVID-19. Since death is a competing risk for hospitalisation, 28 our primary outcome was a composite of hospitalisation or all-cause mortality. Our secondary outcome was allcause mortality alone. Incidence of hospitalisation and all-cause mortality was identified using linked hospital records and the Office for National Statistics. We did not analyse COVID-19 specific mortality (death within 28 days of diagnosis), as the exact dates of diagnosis were not sufficiently robust, and there was uncertainty about accuracy of electronic coding and certification of death nationally with regards to COVID-19 specific mortality in the first few months of the pandemic. Baseline demographics/clinical characteristics of the study cohort were summarised. Categorical variables were summarised as frequencies and percentages, and continuous variables were as means and SD. Since our primary exposure variable, influenza vaccination, is based on observational data (as opposed to allocation by randomisation), the group of patients that received the vaccine is not directly comparable with the group that did not receive the vaccine. The reasons why some people received the vaccination leads to so-called treatment indication bias. We attempted to address this confounding by treatment indication bias within marginal structural regression framework, because conventional regression adjustment may not be sufficient to address this issue. 28 29 Specifically, we estimated the comparative effectiveness of influenza vaccination by applying stabilised inverse probability of treatment weights (IPTW) for the receipt of influenza vaccination. 30 Propensity score was estimated by fitting logistic regression using an a Open access priori list of variables including age, sex, BMI, socioeconomic status (IMD), smoking status, frailty score (electronic frailty index), pre-existing comorbidities (chronic obstructive pulmonary disease, stroke, cancer, depression, peripheral arterial disease, rheumatoid arthritis, atrial fibrillation, dementia, chronic kidney disease, heart failure, learning disability, hypertension, other mental health disorder, cardiovascular disease, epilepsy, asthma, osteoporosis, coronary artery disease, osteoarthritis and diabetes) and the number of prescribed medications. The balance of the propensity score between the two vaccination groups was visually checked using histograms. The numerator for the stabilised IPTW was the marginal probability of vaccination receipt, estimated by fitting intercept-only logistic regression, for the group that received the influenza vaccine. The complement of this probability was the numerator for the group that did not receive the vaccination. Similarly, the propensity score was the denominator for the vaccinated group, while the complement of the propensity score was the denominator for the unvaccinated group. 31 The effect estimates were reported as adjusted ORs and corresponding 95% CIs. Statistical analyses were undertaken using Stata SE V.15.1 (StataCorp, College Station, Texas, USA). To check the robustness of our analysis using stabilised IPTW, we estimated the pooled OR by stratifying the propensity score into deciles (10 equal groups). Mantel-Haenszel weighted pooled OR was estimated using Woolf approximation for the calculation of SE. 32 In further sensitivity analysis, we defined influenza receipt only when it was administered within 6 months of COVID-19 pandemic in the UK (ie, when the vaccination was administered on or after 1 September 2019). Finally, we also performed an additional time-to-event analysis adjusting for month of COVID-19 diagnosis. There was missing data on ethnicity, BMI, smoking, IMD and age. Ethnicity is frequently missing from routinely collected primary care records. We removed this variable from the analyses for multiple reasons: (1) high amount of missing data on ethnicity (32.8%), (2) very low number of outcome events in some subgroups (eg, there was one case of all-cause mortality in the Asian ethnic group) and (3) inadequate reduction of bias by traditional missing data imputation methods 32 33 and (4) this geographical area of the UK has low proportion of ethnic minorities. For other variables with missing data, we used multiple imputation by chained equations assuming missing data were missing at random. We produced five sets of imputed datasets with 100 iterations in the burn-in period. 34 35 The non-missing variables used to impute the missing variables were sex, frailty score and pre-existing comorbidities. Ethical considerations CHIA is an anonymous NHS database, and all individuals have consented for collection of their medical records for inclusion in the database. Findings are reported as per Strengthening the Reporting of Observational Studies in Epidemiology and RECORD guidelines for observational studies of routinely collected data. Anonymised individual level data used in this study was extracted from the CHIA database. Direct requests for data sharing can be made to CHIA. Due to the nature of the study, patients or the public were not involved in the design, conduct, reporting or dissemination of the study A total of 6921 participants were included in the study, of which 774 (11.2%) received a confirmed diagnosis of COVID-19 during the follow-up period, and 6147 (88.8%) were diagnosed clinically. The mean age was 52.4 years (±24.5) and the majority of participants (59.5%) were female. Most participants were white (60.2%); however, ethnicity was not recorded in one-third of cases (32.8%). Approximately half of participants were from the two least deprived socioeconomic status groups (32.8% and 22.9% from IMD groups 5 and 4, respectively). Baseline characteristics of the study cohort are summarised in table 1. A total of 2613 (38%) participants in the cohort had received influenza vaccination between 1 January 2019 and COVID-19 diagnosis (see table 1 ). Compared with the unvaccinated people, those who received influenza vaccination on average were older (mean 61.7 years vs 46.7 years) and frailer (mean frailty score 0.3 vs 0.1). Baselines characteristics of the study cohort, stratified by incidence of the primary study outcome (hospitalisation or all-cause mortality), are displayed in table 2. Of the 2613 people who received an influenza vaccination, 372 (14.2%) died during the follow-up time and 1166 (44.6%) either died or were hospitalised. These proportions were 12.8% (n=553) and 36.8% (n=1584) in the 3755 people without vaccination. Those who died or had a hospital admission were on average older (mean 62.7 years vs 45.5), frailer (frailty score 0.2 vs 0.1) and a lower proportion were reported as white ethnicity (52.0% vs 65.6%), with more having missing ethnicity data (42.7% vs 26.3%). The overall distribution of the propensity score by influenza vaccination status is shown in figure 1 while the mean (SD) of the propensity score by deciles of Open access the propensity score is shown in online supplementary table S1. The distribution shows a fairly strong common support and balanced distribution in the two vaccination groups. In the stabilised IPTW-weighted regression model, influenza vaccination was associated with 15% reduced odds of hospitalisation or all-cause mortality (adjusted OR: 0.85, 95% CI 0.75 to 0.97; p=0.02) and 24% reduced odds of all-cause mortality (adjusted OR: 0.76, 95% CI 0.64 to 0.90) (see table 3 ). Furthermore, on performing additional time-to-event analysis, we found that the results were identical to our primary analysis when we additionally adjusted for the month of COVID-19 diagnosis. In the sensitivity analysis using pooled OR after stratification by deciles of propensity score, the results were almost identical for the primary outcome (adjusted OR: 0.87, 95% CI 0.78 to 0.98, p=0.021) with a slightly stronger effect for all-cause mortality (adjusted OR: 0.59, 95% CI 0.51 to 0.69, p<0.001) (see table 3 ). In the sensitivity analysis with a modified definition of influenza vaccination (administered on or after 1 September 2020), the effects were almost identical to our primary analysis (adjusted OR: 0.85, 95% CI 0.74 to 0.97, p=0.016 for hospitalisation or all-cause mortality, and 0.78, 95% CI 0.66 to 0.92, p=0.004 for all-cause mortality). In this sample of 6921 participants who received a diagnosis of COVID-19 during the first wave of the pandemic (1 January-31 July 2019), we found that having prior influenza vaccination was associated with a 15%-24% lower odds of hospitalisation or all-cause mortality. Our study is consistent with observational data from outside the UK demonstrating a negative association between influenza vaccination and COVID-19 acquisition, severity and mortality. 12 14 18 35 36 An recent Italian retrospective study demonstrated a modest negative correlation between the percentage of adults aged >65 years who had received influenza vaccination and the percentage of COVID-19 deaths from each region of Italy (r: −0.59). 15 Similarly, a US cohort of adults aged ≥65 years demonstrated that a 10% increase in vaccination uptake was associated, on average, with a 28% decrease in the COVID-19 mortality rate. 22 Another study examined influenza vaccination data from 34 countries (using Organisation for Economic Cooperation and Development data) with COVID-19 mortality from worldometer data and reported a negative correlation between influenza vaccination status and COVID-19 mortality (R 2 0.338). 17 These studies were limited by the use of aggregated data, however. Two retrospective studies using patient-level data include a large Italian cohort study in which receipt of influenza vaccination was associated with 11% lower odds of COVID-19 diagnosis, and among a subgroup of those aged ≥65 years vaccinated in the first half of the vaccination programme, a 44% and 30% lower odds of hospitalisation and mortality, respectively. Another large retrospective cohort study in Brazil of 53 752 hospitalised patients with COVID-19 found that receipt of recent influenza vaccination was associated with 7%, 17% and 16% lower odds of intensive care admission, requiring invasive respiratory support and mortality, respectively. Finally, there is one prospective observational study in the USA that includes individual patient level data through a registry of 11 672 people tested for COVID-19 (of whom, 7% tested positive). 24 It reports that influenza vaccination was associated with a lower risk of acquiring COVID-19 (p<0.001, OR not reported) but does not report on subsequent COVID-19 outcomes, including mortality. Possible explanations for these findings include a direct protective effect against COVID-19 from influenza vaccination. Growing evidence, particularly from work studying BCG and measles vaccination, 9 11 12 as well as more recently from influenza vaccination, 13 14 demonstrates that vaccination may have a role in wider protection against unrelated pathogens, in addition to inducing adaptive/innate immunity against specific pathogens. The underlying mechanisms for these non-specific effects remain poorly understood. One suggested mechanism is the induction of innate immune responses (including monocytes and natural killer cells), so-called trained innate immunity, following live attenuated vaccination, in a way that is independent of memory T or B cells. 9 11-14 This trained immunity may confer non-specific protection against different pathogens, inducing upregulation of pattern recognition receptors and the secretion of proinflammatory cytokines through epigenetic and metabolic reprogramming. [11] [12] [13] [14] Another proposed mechanism is a process called 'emergency Open access granulopoiesis'. 37 This is based on the recent observation that following BCG vaccination in newborns; for example, there is induction of the growth factor granulocyte colony-stimulating factor which, in turn, leads to the rapid production of neutrophils, resulting in additional protection from subsequent infection. 38 Finally, influenza vaccination may also offer additional indirect protection against COVID-19 by reducing the burden of additional respiratory illness among both individual patients and healthcare services. [18] [19] [20] Strengths and limitations To our knowledge, this is the first study to examine influenza vaccination in relation to all-cause mortality and hospitalisation based on individual-patient data linked to demographic and clinical records in the UK. We applied robust analytic techniques to reduce the treatment indication bias that is inherent to observational studies. However, this does not completely rule out the possibility of residual and unknown confounding. A strength of the study is that it was drawn from a large database of individuals 160 GP practices across South West England; however, the cohort does include low representation of ethnic minority groups, and a relatively large proportion were from high socioeconomic backgrounds. Another limitation is our inability to analyse the effects by ethnicity due to large amount of missing data. Data being from a centralised single-payer healthcare system may limit the generalisability of these findings to other parts of the world. The case definition of COVID-19 disease has also evolved over time, especially as testing has become more readily available in recent months. We included clinically suspected cases of COVID-19 alongside confirmed cases, in accordance with guidance from PHE during the first wave of the pandemic. Recent studies have shown that clinical suspicion of COVID-19 among GPs does closely match with confirmed cases and subsequent clinical outcomes, but it is still plausible that not all those with clinical cases had COVID-19. 27 The clinical condition of patients at the time of diagnosis is also likely to have varied over the course of the study period (again as testing has become more readily available), as has the clinical management of patients, as clinicians have become more familiar with the COVID-19 over time. An additional time-to-event analysis in which we adjusted for the month of COVID-19 diagnosis showed identical results to our primary analysis, however. Our sample was also restricted to individuals with complete consultation notes over the study period. It is possible that individuals with incomplete records, mild symptoms only (or were asymptomatic) or attended hospital directly without first attending general practice may be under-represented in our sample. There is also a well-recognised 'frailty bias', 39 whereby as people get older and frailer, vaccination rates tends to drop, and it is therefore plausible that this could have influenced our findings. However, we tried to take this into consideration by including frailty index, age and comorbidities within the regression model. Finally, the type of influenza vaccination given was not recorded in the database, and we are therefore unable to explore whether there is a difference between inactivated and live-attenuated vaccination on COVID-19 outcomes. We used logistic regression for our primary analysis because the date of COVID-19 diagnosis, hospitalisation or deaths was not precise in the database. However, as stated above, additional time-to-event analysis adjusted for the month of COVID-19 diagnosis showed identical results. Furthermore, exact date of death was not always available for patients who died (only months and years of death were reported). Since the follow-up time was small, this is unlikely to have changed the estimates Figure 1 Distribution of the propensity score for patients who received seasonal influenza vaccination (n=2613) compared with those who did not receive the vaccination (n=4308) since January 2019. substantially. 40 To verify this, we imputed the date as the 15th of the respective months and ran additionally sensitivity analyses on each of the outcomes using stabilised IPTW-weighted Cox's proportional hazards model to estimate the adjusted HR and 95% CI. This robustly supported our findings; adjusted HR for the primary outcome (hospitalisation or all-cause mortality) and the secondary outcome (all-cause mortality) was 0.88 (95% CI 0.79 to 0.98, p=0.024) and 0.79 (95% CI 0.66 to 0.94, p=0.008), respectively. Our results suggest that influenza vaccination may have the potential to offer an additional means of mitigating serious adverse complications of COVID-19. Further research should aim to validate these findings in larger more diverse cohorts. Additional work is also needed to examine differences between influenza vaccine subtypes and whether the duration from vaccination administration to COVID-19 acquisition is important. In people with COVID-19 during the first wave of the pandemic, we found that influenza vaccination was associated with a 15%-24% lower odds of hospitalisation or all-cause mortality. An expansion of the influenza vaccination programme is needed as an integral component of the ongoing response to the COVID-19 pandemic. *IPTW: inverse probability of treatment weighting based on the propensity score for influenza vaccination stabilised by the marginal probability of vaccination in the study population. †The propensity score was divided into deciles, and the OR for the association between influenza vaccination and the outcome was estimated in each of the deciles. Mantel-Haenszel-weighted pooled OR was estimated using Woolf approximation for the calculation of SE. CHIA, Care and Health Information Analytics. COVID-19 coronavirus pandemic Multiple imputation for nonresponse in surveys Short-and potential longterm adverse health outcomes of COVID-19: a rapid review The COVID-19 Pandemic and the $16 Trillion Virus Indications for healthcare surge capacity in European countries facing an exponential increase in coronavirus disease (COVID-19) cases Potential effects of vaccinations on the prevention of COVID-19: rationale, clinical evidence, risks, and public health considerations Safety and immunogenicity of ChAdOx1 nCoV-19 vaccine administered in a prime-boost regimen in young and old adults (COV002): a Open access single-blind, randomised, controlled, phase 2/3 trial COVID-19 vaccines: no time for complacency Can existing live vaccines prevent COVID-19? Association of BCG, DTP, and measles containing vaccines with childhood mortality: systematic review Bcg vaccination protects against experimental viral infection in humans through the induction of cytokines associated with trained immunity Bcg-Induced cross-protection and development of trained immunity: implication for vaccine design Elevated frequency of gamma interferon-producing NK cells in healthy adults vaccinated against influenza virus The effect of influenza vaccination on trained immunity: impact on COVID-19 Epidemiological evidence for association between higher influenza vaccine uptake in the elderly and lower COVID-19 deaths in Italy Reflection on lower rates of COVID-19 in children: does childhood immunizations offer unexpected protection? Considering interim interventions to control COVID-19 associated morbidity and Mortality-Perspectives COVID-19, seasonal influenza and measles: potential triple burden and the role of flu and MMR vaccines Seasonal influenza vaccination during a pandemic Influenza immunization and COVID-19 Association between exposure to influenza vaccination and COVID-19 diagnosis and outcomes Influenza vaccination and COVID19 mortality in the USA Inactivated trivalent influenza vaccination is associated with lower mortality among patients with COVID-19 in Brazil Individualizing risk prediction for positive coronavirus disease 2019 testing: results from 11,672 patients Emergence of a novel coronavirus (COVID-19): protocol for extending surveillance used by the Royal College of general practitioners research and surveillance centre and public health England The use of primary care big data for COVID-19 research: a consensus statement from the COVID-19 primary care database Consortium Excess mortality from COVID-19 in an English sentinel network population Living risk prediction algorithm (QCOVID) for risk of hospital admission and mortality from coronavirus 19 in adults: national derivation and validation cohort study A most stubborn bias: no adjustment method fully resolves confounding by indication in observational studies Marginal structural models and causal inference in epidemiology Constructing inverse probability weights for marginal structural models A Propensity-scorebased fine stratification approach for confounding adjustment when exposure is infrequent Imputing missing Race/ Ethnicity in pediatric electronic health records: reducing bias with use of U.S. census location and Surname data Multiple imputation for nonresponse in surveys Multiple imputation by chained equations: what is it and how does it work? Joint modelling rationale for chained equations The non-specific and sexdifferential effects of vaccines Bcg vaccination-induced emergency granulopoiesis provides rapid protection from neonatal sepsis The importance of frailty in the assessment of influenza vaccine effectiveness against Influenza-Related hospitalization in elderly people A comparison of the logistic risk function and the proportional hazards model in prospective epidemiologic studies Patient consent for publication Not required.Ethics approval Ethical and governance approval for this specific study was obtained from the University of Southampton (ERGO 61740), and Care and Health Information Exchange Information Governance Group (CHIE IGG).Provenance and peer review Not commissioned; externally peer reviewed.Data availability statement Data are available on reasonable request. Anonymised individual level data used in this study was extracted from the Care and Health Information Analytics (CHIA) database. Direct requests for data sharing can be made to CHIA.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.Open access This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https:// creativecommons. org/ licenses/ by/ 4. 0/.