Risk of severe COVID-19 outcomes associated with immune-mediated inflammatory diseases and immune modifying therapies: a nationwide cohort study in the OpenSAFELY platform =========================================================================================================================================================================== * Brian MacKenna * Nicholas A. Kennedy * Amir Mehkar * Anna Rowan * James Galloway * Kathryn E. Mansfield * Katie Bechman * Julian Matthewman * Mark Yates * Jeremy Brown * Anna Schultze * Sam Norton * Alex J. Walker * Caroline E Morton * David Harrison * Krishnan Bhaskaran * Christopher T. Rentsch * Elizabeth Williamson * Richard Croker * Seb Bacon * George Hickman * Tom Ward * Simon Davy * Amelia Green * Louis Fisher * William Hulme * Chris Bates * Helen J. Curtis * John Tazare * Rosalind M. Eggo * David Evans * Peter Inglesby * Jonathan Cockburn * Helen I. McDonald * Laurie A. Tomlinson * Rohini Mathur * Angel YS Wong * Harriet Forbes * John Parry * Frank Hester * Sam Harper * Ian J. Douglas * Liam Smeeth * Charlie W Lees * Stephen JW Evans * Ben Goldacre * Catherine Smith * Sinéad M. Langan ## ABSTRACT **Background** It is unclear if people with immune-mediated inflammatory diseases (IMIDs) (joint, bowel and skin) and on immune modifying therapy have increased risk of serious COVID-19 outcomes. **Methods** With the approval of NHS England we conducted a cohort study, using OpenSAFELY, analysingroutinely-collected primary care data linked to hospital admission, death and previously unavailable hospital prescription data. We used Cox regression (adjusting for confounders) to estimate hazard ratios (HR) comparing risk of COVID-19-death, death/critical care admission, and hospitalisation (March to September 2020) in: 1) people with IMIDs compared to the general population; and 2) people with IMIDs on targeted immune modifying drugs (e.g., biologics) compared to standard systemic treatment (e.g., methotrexate). **Findings** We identified 17,672,065 adults; of 1,163,438 (7%) with IMIDs, 19,119 people received targeted immune modifying drugs, and 200,813 received standard systemics. We saw evidence of increased COVID-19-death (HR 1.23, 95%CI 1.20, 1.27), and COVID-19 hospitalisation (HR 1.32, 95%CI 1.29, 1.35) in individuals with IMIDs overall compared to individuals without IMIDs of the same age, sex, deprivation and smoking status. We saw no evidence of increased COVID-19 deaths with targeted compared to standard systemic treatments (HR 1.03, 95%CI 0.80, 1.33). There was no evidence of increased COVID-19-related death in those prescribed TNF inhibitors, IL-12/23, IL7, IL-6 or JAK inhibitors compared to standard systemics. Rituximab was associated with increased COVID-19 death (HR 1.68, 95%CI 1.11, 2.56); however, this finding may relate to confounding. **Interpretation** COVID-19 death and hospitalisation was higher in people with IMIDs. We saw no increased risk of adverse COVID-19 outcomes in those on most targeted immune modifying drugs for IMIDs compared to standard systemics. **Evidence before this study** We searched PubMed on May 19th, 2021, using the terms “COVID-19”, “SARS-CoV-2” and “rheumatoid arthritis”, “psoriatic arthritis” “ankylosing spondylitis”, “Crohn’s disease” “ulcerative colitis” “hidradenitis suppurativa” and “psoriasis”, to identify primary research articles examining severe COVID-19 outcome risk in individuals with immune-mediated inflammatory diseases (IMIDs) and those on immune modifying therapy. The studies identified (including matched cohort studies and studies in disease-specific registries) were limited by small sample sizes and number of outcomes. Most studies did not show a signal of increased adverse COVID-19 outcomes in those on targeted therapies, with the exception of rituximab. Additionally, disease-specific registries are subject to selection bias and lack denominator populations. **Added value of the study** In our large population-based study of 17 million individuals, including 1 million people with IMIDs and just under 200,000 receiving immune modifying medications, we saw evidence that people with IMIDs had an increased risk of COVID-19-related death compared to the general population after adjusting for potential confounders (age, sex, deprivation, smoking status) (HR 1.23, 95%CI 1.20, 1.27). We saw differences by IMID type, with COVID-19-related death being increased by the most in people with inflammatory joint disease (HR 1.47, 95%CI 1.40, 1.54). We also saw some evidence that those with IMIDs were more likely, compared to the general population, to have COVID-19-related critical care admission/death (HR 1.24, 95%CI 1.21, 1.28) and hospitalisation (HR 1.32, 95%CI 1.29, 1.35). Compared to people with IMIDs taking standard systemics, we saw no evidence of differences in severe COVID-19-related outcomes with TNF inhibitors, IL-17 inhibitors, IL-12/23 inhibitors, IL-6 inhibitors and JAK inhibitors. However, there was some evidence that rituximab was associated with an increased risk of COVID-19-related death (HR 1.68, 95%CI 1.11, 2.56) and death/critical care admission (HR 1.92, 95%CI 1.31, 2.81). We also saw evidence of an increase in COVID-19-related hospital admissions in people prescribed rituximab (HR 1.59, 95%CI 1.16, 2.18) or JAK inhibition (HR 1.81, 95%CI 1.09, 3.01) compared to those on standard systemics, although this could be related to worse underlying health rather than the drugs themselves, and numbers of events were small. This is the first study to our knowledge to use high-cost drug data on medicines supplied by hospitals at a national scale in England (to identify targeted therapies). The availability of these data fills an important gap in the medication record of those with more specialist conditions treated by hospitals creating an important opportunity to generate insights to these conditions and these medications **Implications of all of the available evidence** Our study offers insights into future risk mitigation strategies and SARS-CoV-2 vaccination priorities for individuals with IMIDs, as it highlights that those with IMIDs and those taking rituximab may be at risk of severe COVID-19 outcomes. Critically, our study does not show a link between most targeted immune modifying medications compared to standard systemics and severe COVID-19 outcomes. However, the increased risk of adverse COVID-19 outcomes that we saw in people with IMIDs and those treated with rituximab merits further study. Keywords * immune mediated inflammatory disease * COVID-19 death * COVID-19 critical care admission * COVID-19 hospitalisation * rheumatoid arthritis * psoriatic arthritis * ankylosing spondylitis * inflammatory bowel disease * Crohn’s disease * ulcerative colitis * psoriasis * hidradenitis suppurativa * immune modifying drugs * TNF inhibitor * IL-12/23 inhibitor * I-17 inhibitor * IL-6 inhibitor * JAK inhibitor * rituximab ## BACKGROUND Although most people with coronavirus disease 2019 (COVID-19) experience mild symptoms, 15% develop pneumonia requiring hospitalisation, and 5% progress to severe disease (systemic hyperinflammation, major coagulation problems and acute respiratory distress, with risk of multi-organ failure and death) (of unvaccinated individuals).1 It is unclear to what extent having immune-mediated inflammatory diseases (IMIDs) – including inflammatory joint (rheumatoid arthritis [RA], psoriatic arthritis [PsA], ankylosing spondylitis [AS]), bowel (Crohn’s disease, ulcerative colitis [UC]), and skin (psoriasis, hidradenitis suppurativa [HS]) diseases – or being on targeted immune modifying therapies for IMID treatment, increases severe COVID-19 outcomes.2 Different strategies have been implemented worldwide to limit infection spread. Policies were initially informed by expert opinion and subsequently by risk estimates for COVID-19 death.3 Early in the pandemic, the UK government issued guidance through the National Health Service (NHS), known as “shielding”, advising individuals considered at high risk of severe COVID-19 to stay at home, minimising face-to-face contacts. Shielding advice targeted individuals identified using health records as having specific immunosuppressive diseases or prescriptions for immune modifying drugs.4 Many people with IMIDs taking immune modifying therapies were considered “high risk”, based on prior evidence in other infections.5,6 Disease-specific registries have been used to evaluate COVID-19 risks in individuals with IMIDs and those prescribed immune modifying therapies; results have largely been reassuring.7–9 However, disease-specific registries are limited by small sample sizes, inherent selection bias, and lack of denominators. OpenSAFELY-TPP is a new secure analytics platform for electronic health records (EHRs) to deliver pandemic-related research. We aimed to investigate risks of severe COVID-19 outcomes in people with IMIDs and those on targeted immune modifying therapies using English population-based EHR data linked to a new unique national dataset containing information on targeted immune modifying therapies. ## METHODS ### Study design We conducted a cohort study using electronic health record (EHR) data. We compared risks of COVID-19-related death, critical care admission or death, and hospital admission in: 1) people with IMIDs compared to the general population; and 2) in people with IMIDs prescribed targeted immune modifying drugs compared to those with IMIDs taking standard systemic immune modifying drugs. ### Data source We used primary care records managed by the GP software provider TPP linked to: Office for National Statistics (ONS) death data, COVID-19 testing data, and a unique national hospital medications dataset (including “high-cost drugs”, typically newer, more expensive medications [e.g., adalimumab] supplied by hospitals for ongoing disease management, **Supplementary Text S1**). We accessed all data through OpenSAFELY, a data analytics platform created by our team for NHS England ([https://opensafely.org](https://opensafely.org)). OpenSAFELY provides a secure software interface allowing analysis of pseudonymized primary care records in near real-time within the EHR vendor’s highly secure data centre, avoiding the need for data transfer off-site (minimising re-identification risk). Similarly, pseudonymized datasets from other data providers are securely provided by the EHR vendor and linked to primary care data. The dataset analysed within OpenSAFELY was based on 24 million people currently registered with GP surgeries. ### Study population Our overall study population included adults aged 18–110 years (March 2020), registered with TPP practices with at least 12 months of primary care records prior to March 2020 (**Figure 1**). We followed individuals from March 1st, 2020 (UK SARS-Cov-2 outbreak start) until 30th September 2020 (study end), or until the specific outcome under analysis (i.e., COVID-19-related: death, critical care admission or death, or hospitalisation). ![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/05/06/2021.09.03.21262888/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/F1) Figure 1. Study design We compiled diagnostic and therapeutic code lists (in machine readable languages, e.g., SNOMED-CT, NHS dictionary of medicines and devices) for all study variables (exposures, outcomes, covariates). Detailed information on compilation and sources of code lists are freely available for inspection and re-use ([https://codelists.opensafely.org/](https://codelists.opensafely.org/)). ### Exposures Exposures were: 1) IMIDs: inflammatory joint disease (RA, PsA, AS), inflammatory bowel disease (IBD) (Crohn’s disease, UC, IBD unclassified), and inflammatory skin disease (psoriasis, HS); and 2) prescription of systemic immune modifying medication by general practitioners or supplied by hospitals through high-cost drugs procedures. We identified people with IMIDs using diagnostic morbidity codes in primary care during the three years prior to March 1st, 2020. People with IMID diagnoses in multiple categories (e.g., inflammatory joint and skin disease) contributed to comparisons with the general population for all IMID categories they had records for (e.g., individuals with PsA and psoriasis, contributed to both joint and skin disease). Immune modifying medications were categorised as standard systemic therapy and targeted therapy. Standard systemic therapies included leflunomide, methotrexate, mycophenolate mofetil or mycophenolic acid, ciclosporin, sulphasalazine, mercaptopurine, thioguanine, and azathioprine. Targeted therapies included biologic therapies: TNF inhibitors (etanercept, adalimumab, golimumab, certolizumab, infliximab, and biosimilar versions of these medications), IL-17 inhibitors (secukinumab, ixekizumab, brodalumab), IL-12/23 inhibitors (ustekinumab, guselkumab, risankizumab, tildrakizumab), IL-6 inhibitors (tocilizumab, sarilumab), B-cell depletion therapy (rituximab), and novel small molecules: Janus Kinase (JAK) inhibitors (baricitinib, tofacitinib).10–15 Individuals treated with both systemic therapy and targeted therapies were considered exposed to targeted therapies. We identified standard systemic therapies using primary care prescribing data. We identified targeted immune modifying medications using high-cost drugs invoicing **(Supplementary Text S1)**. Drug exposure was defined by at least one prescription, or delivery of medication to an individual before March 1st, 2020 (date chosen as some of medications were either specifically used, or stopped, due to the pandemic). For each individual, we defined drug exposure based on the closest drug recorded prior to the study start (March 1st 2020) allowing for a maximum of 6 months before the start of the study for all agents apart from rituximab, where we permitted a 12-month exposure window (to allow for the frequency of treatment and the longer duration of response of rituximab, including known prolonged effects on vaccine responses).16,17 We could not evaluate medication switching or adherence during the follow-up period, as data were only available up to March 2020. ### Outcomes Outcomes were COVID-19-related death, critical care admission or death, and hospitalisation. We identified COVID-19-related deaths based on records of COVID-related ICD-10 codes (U071, U072) anywhere on death certificates. We used COVID-19-related critical care admission (Intensive Care National Audit & Research Centre18) or death as a combined endpoint to reflect individuals with severe COVID-19 who died without being admitted to a critical care unit for critical care admission. We identified COVID-19-related hospitalisation as a positive polymerase chain reaction (PCR) test less than 28 days before admission and, up to 5 days post admission to exclude nosocomial infection. ### Covariates We selected potential confounders and mediators *a priori* based on clinical knowledge and prior evidence.2 In the relationship between IMIDs and severe COVID-19 outcomes, we considered age (categorical variable), sex, deprivation (using quintiles of Index of Multiple Deprivation),19,20 and smoking status to be potential confounders; we considered BMI, cardiovascular disease, diabetes and current glucocorticoid use to be potential mediators. In the relationship between (choice of) immune modifying therapy and severe COVID-19 outcomes we considered age, sex, deprivation, smoking status, BMI, specific IMID (inflammatory joint, bowel and skin disease), cardiovascular disease, cancer (excluding non-melanoma skin cancer), stroke, end stage renal failure, chronic liver disease, chronic respiratory disease and diabetes mellitus as potential confounders; we considered current glucocorticoid use as a potential mediator. Details of covariate definitions are included in **Supplementary Text S2**. Figures representing assumed relationships between covariates, primary exposures and outcomes are in **Supplementary Figures S1 and S2.** ### Statistical methods We initially described characteristics of the general population, people with IMIDs, and those with IMIDs prescribed immune modifying therapy. We used Cox regression to estimate hazard ratios (HR) (95%CI), comparing people with IMIDs to the general population, and people with IMIDs on standard systemic drugs to those on targeted therapies. We adjusted models for confounding based on assumptions inherent in our conceptual frameworks (**Supplementary Figures S1, S2**). For comparisons between those with IMIDs to the general population, we adjusted for age, sex, deprivation, and smoking. For comparisons between those on targeted therapies to those on standard systemic drugs, we adjusted for age, sex, deprivation, obesity (BMI>=30), smoking, IMIDs (bowel, joint, skin), and comorbidities (cardiovascular disease, diabetes mellitus, cancer, stroke, end stage renal failure, chronic liver disease, chronic respiratory disease). We also fitted models adjusted for additional covariates considered to be mediators. We tested Cox model assumptions using Schoenfeld residuals. #### Sensitivity analyses We repeated our main analyses in sensitivity analyses assessing robustness of our findings (**Supplementary Table S1, Text S3**). We considered immune-mediated inflammatory disease severity and degree of shielding to be potential unmeasured confounders of associations between specific immune-modifying therapy and COVID-19 outcomes. We therefore conducted quantitative bias analysis using E-values to assess how strongly associated unmeasured confounders would need to be with exposure and outcome to potentially fully explain observed non-null associations (i.e., association adjusted for both measured covariates and the unmeasured confounder would be null).21 ### Software and reproducibility We used Python for data management, and Stata 16/Python for analyses. All code used for data management and analyses is available online ([https://github.com/opensafely/immunosuppressant-meds-research](https://github.com/opensafely/immunosuppressant-meds-research)) including all iterations of the pre-specified study protocol archived with version control. ### Patient and public involvement Patients were not formally involved in developing this study design, which was developed rapidly during a pandemic. We invite patients or members of the public to contact us through our website [https://opensafely.org/](https://opensafely.org/) regarding this study or the broader OpenSAFELY project. Our protocol and paper preprint have been sent to Crohn’s and Colitis UK, Psoriasis Association, National Rheumatoid Arthritis Society, and Versus Arthritis for review from an expert-patient perspective. ### Ethics The study was approved by the Health Research Authority (Research Ethics Committee reference 20/LO/0651) and the London School of Hygiene and Tropical Medicine’s Ethics Board (Reference 21863). ## RESULTS Of 17,672,065 people in the overall study population (**Supplementary Figure S3**), 1,163,438 (7%) had an IMID diagnosis (**Table 1**). Of those with IMIDs: 1) 23% (n=272,452) had inflammatory joint disease (RA 16% n=183,485, PsA 5% n=54,593, AS 3% n=35,138); 2) 17% (n=199,037) had IBD (Crohn’s 35% n=69,788, UC 51% n=100,617, unclassified IBD 16% n=32,093); and 3) 66% (n=769,816) had inflammatory skin disease (psoriasis 90% n=693,178, HS 10% n=76,746). View this table: [Table 1:](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T1) Table 1: Descriptive Characteristics of general population and people with IMIDs. ### IMIDs compared to the general population Compared to the general population, people with IMIDs were older (>=70 years: 25% versus 17% in the general population), more likely to be female (55% versus 50%), of white ethnicity (71% versus 64%), obese (BMI>=30, 30% versus 22%), and with more comorbidities (e.g., cardiovascular disease: 11% versus 7%: diabetes mellitus: 16% versus 11%) (**Table 1**). There were differences between individuals with inflammatory joint, skin and bowel diseases, for example, individuals with inflammatory joint disease were older compared to those with IBD or inflammatory skin disease (>=70 years: 35.3% joint, 23.4% IBD, 21.1% skin). #### COVID-19-related death After adjusting for age and sex, compared to the general population, people with IMIDs had greater risk of COVID-19-related death (HR 1.27, 95%CI 1.23, 1.31). Evidence of association between IMIDs and COVID-19-related death remained after additionally adjusting for confounders (age, sex, deprivation and smoking status, HR 1.23, 95%CI 1.20, 1.27) and after further adjusting for potential mediators (BMI, cardiovascular disease, diabetes and current glucocorticoid use, HR 1.15, 95%CI 1.11, 1.18). (**Figure 2, Supplementary Table S2**). ![](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/05/06/2021.09.03.21262888/F2/graphic-3.medium.gif) [](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/F2/graphic-3) ![](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/05/06/2021.09.03.21262888/F2/graphic-4.medium.gif) [](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/F2/graphic-4) Figure 2. Forest plot of hazard ratios (HRs) for COVID-19-related death, critical care admission/death and hospitalisation for IMID vs general population After adjusting for age and sex, compared to the general population we saw increased COVID-19-related death in people with inflammatory joint (HR 1.51, 95%CI 1.44, 1.58), bowel (HR 1.15,95%CI 1.07, 1.24), and skin (HR 1.16, 95%CI 1.11, 1.20) diseases. After further adjusting for potential confounders, evidence for association between specific IMID types and COVID-19 related death persisted for all IMID types and was greatest for inflammatory joint disease (HR 1.47, 95%CI 1.40, 1.54), with smaller effect estimates for inflammatory skin (HR 1.12, 95%CI 1.08, 1.17) or bowel (HR 1.12, 95%CI 1.04, 1.21) disease, and further attenuation after adjusting for potential mediators. #### COVID-19-related critical care admission or death Individuals with IMIDs had greater risk of COVID-19-related critical care admission or death compared to the general population, which persisted after adjusting for confounders (HR 1.24, 95%CI 1.21, 1.28) and further adjusting for mediators (HR 1.16, 95%CI 1.12, 1.19). Compared to the general population, there was evidence of increased COVID-19-related critical care admission or death in people with inflammatory joint (confounder adjusted HR 1.46, 95%CI 1.39, 1.52), skin (confounder adjusted HR 1.15, 95%CI 1.10, 1.19), and bowel (confounder adjusted HR 1.13, 95%CI 1.06, 1.22) diseases. #### COVID-19-related hospitalisation Compared to the general population, people with IMIDs had greater risk of COVID-19-related hospitalisation, which remained after adjusting for potential confounders (HR 1.32, 95%CI 1.29, 1.35) and mediators (HR 1.20, 95%CI 1.17, 1.23). Compared to the general population, risk of COVID-19 hospitalisation was increased in all IMID categories including inflammatory joint (confounder adjusted HR 1.53, 95%CI 1.47, 1.59), skin (adjusted HR 1.21, 95%CI 1.17, 1.24), and bowel (adjusted HR 1.31, 95%CI 1.24, 1.38) disease. #### Sensitivity analyses Results from sensitivity analyses were broadly similar to the main analysis. (**Supplementary Tables S3-S4**). ### Targeted therapy compared to standard systemic therapy Of those with IMIDs, 200,813 individuals (17%) were prescribed either standard systemic therapy (90%, n=181,694/200,813) or targeted immune modifying therapy (10%, n=19,119/200,813) (**Table 2, Supplementary Table S5**). Compared to those on standard systemic therapy, individuals receiving targeted therapy were younger (<40 years: 22% versus 14%) and less likely to have comorbidities (e.g., cardiovascular disease). The most commonly prescribed targeted therapies were TNF inhibitors (71%, n=13,524/19,119), followed by rituximab (10%, n=1,998/19,119); IL-12/23 inhibitors, (7%, n=1,379/19,119); IL-17 inhibitors (5%, n=1,036/19,119); JAK inhibitors (5%, n=871/19,119) and IL-6 inhibitors (4%, n=758/19,119). View this table: [Table 2:](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T2) Table 2: Descriptive characteristics of IMID population on targeted and standard systemic immune modifying therapy. #### COVID-19-related death After adjusting for potential confounders (age, sex, deprivation, BMI, IMIDs [bowel, joint, skin], cardiovascular disease, cancer, stroke, diabetes mellitus), we saw no evidence of differences in COVID-19-related death in individuals on targeted therapy overall compared to individuals on standard systemic therapy (HR 1.03, 95%CI 0.80,1.33) (**Figure 3, Supplementary Table S6**). ![](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/05/06/2021.09.03.21262888/F3/graphic-7.medium.gif) [](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/F3/graphic-7) ![](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/05/06/2021.09.03.21262888/F3/graphic-8.medium.gif) [](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/F3/graphic-8) Figure 3. Forest plot of hazard ratios (HRs) for COVID-19 death, critical care admission/death and hospitalisation for standard systemic vs targeted immunosuppression Compared to individuals on standard systemic therapy, there was no observed increased risk of COVID-19 related death in individuals on TNF inhibitors, IL-12/23 inhibitors, IL7 inhibitors, JAK inhibitors or IL-6 inhibitors (although confidence intervals were wide). Compared to individuals on standard systemic therapy, people receiving rituximab had an increased risk of COVID-19-related death (HR 1.68, 95%CI 1.11, 2.56), based on 24 deaths in the rituximab group. #### COVID-19-related critical care admission or death After adjusting for confounders, compared to individuals on standard systemic therapy, there was no increase in COVID-19-related critical care admission or death in individuals on TNF inhibitors, IL-12/23 inhibitors, IL-17 inhibitors, JAK inhibitors or IL-6 inhibitors. People receiving rituximab had increased risk of critical care admission or death (HR 1.92, 95%CI 1.31, 2.81), based on 29 events in the rituximab group. #### COVID-19-related hospitalisation After adjusting for confounders, compared to individuals on standard systemic therapy, there was no difference in risk of COVID-19-related hospitalisations in individuals on TNF inhibitors, IL-12/23 inhibitors, IL-17 inhibitors, or IL-6 inhibitors. Compared to standard systemics, we saw increased risk of COVID-19-related hospitalisation in people receiving rituximab (HR 1.59, 95%CI 1.16, 2.18; 40 events), and JAK inhibitors (HR 1.81, 95%CI 1.09, 3.01; 15 events). #### Sensitivity analyses Excluding people with haematological cancers and organ transplants attenuated the effect estimate for rituximab (HR 1.54, 95%CI 0.95, 2.49; 18 events). Otherwise, results from sensitivity analyses were similar to the main analysis (**Supplementary Tables S7-S11**). In quantitative bias analysis comparing individuals with IMIDs taking rituximab or JAK inhibitors to standard systemics we note an unmeasured confounder moderately associated with both exposure and outcome could potentially explain associations of rituximab (by at least a risk ratio of 1.68) and JAK inhibitors (by at least a risk ratio of 1.81) with adverse COVID-19 outcomes. ## DISCUSSION In this large population-based study using data fromOpenSAFELY, we found that people with IMIDs experience higher rates of COVID-19-related death, critical care admission and death, and hospitalisation than people without IMIDs of the same age, sex, deprivation and smoking status. We demonstrated that compared to standard systemic immune modifying therapies for IMIDs, there was no increased risk of COVID-19-related death in those prescribed TNF, IL-12/23, IL7, IL-6 or JAK inhibitors. However, rituximab was associated with an increased risk of death (HR 1.64, 95%CI 1.08, 2.49) and critical care admission, which may relate to confounding. ### Findings in context Our findings — suggesting that people with IMIDs were at increased risk of COVID-19-related death compared to people without IMIDs of the same age, sex, deprivation, and smoking status — have not previously been reported. This study provides granular data on a range of IMIDs but is consistent with previous less detailed analyses.22 Our finding that those with IMIDs were more likely to be hospitalised for COVID-19 compared to the general population is consistent with Canadian and Danish cohort studies 232423 and reports of adverse COVID-19 outcomes for people with specific IMIDs.25 However, factors leading to adverse COVID-19 outcomes are likely multifactorial, with possible explanations including more severe symptoms following COVID-19, unmeasured confounders, better access to care, or a lower physician threshold for admission in those on immune modifying drugs. There are no comparable population-based data sources including information on COVID-19 outcomes in those on targeted therapies for IMIDs. Our observation that those on targeted therapies do not have an increased risk of COVID-19 related death is consistent with data from international registries.7–9 Compared to our findings, a recent meta-analysis using data from registries included 2,766 individuals with autoimmune diseases and COVID-19 diagnoses, reported higher rates of hospitalisation and death in people prescribed combination standard systemic therapy and biologics or JAK inhibitors, but lower rates in those prescribed TNF inhibitor monotherapy.26 The Global Rheumatology Alliance reported no increase in COVID-19-related death with biologic therapies or JAK inhibitors compared to methotrexate monotherapy, but an increase in COVID-19-related death with rituximab.8 The Surveillance Epidemiology of Coronavirus Under Research Exclusion for Inflammatory Bowel Disease (SECURE-IBD) reported no association in COVID-19-related death, critical care or hospitalisation in people on TNF inhibitors compared to those not prescribed TNF inhibitor therapy.7 Our observed increased risk of critical care admission or death in people treated with rituximab is consistent with previous reports of increased mortality in people treated with B cell depleting agents (e.g., rituximab for oncology indications).27 Explanations could include mechanistic roles of B cells in COVID-19 immune responses, unmeasured characteristics of the RA population treated with rituximab or confounding by severity (i.e., rituximab is preferentially prescribed for seropositive disease, with extra-articular manifestations including pulmonary fibrosis) or by cancer.16 We were underpowered to assess effects of regular tocilizumab on COVID-19 outcomes, although trial data has shown benefit in hospitalised and critically ill patients.28,29 ### Strengths and weaknesses The key strengths of this study are the scale and completeness of underlying raw EHR data. The OpenSAFELY platform accesses an unprecedented scale of data; the full dataset of all raw, single-event-level clinical events for all individuals at 40% of all GP practices in England, including all tests, treatments, diagnoses, and clinical and demographic information linked to various sources of hospital data including, for the first time, a comprehensive dataset of medications supplied by hospitals. We recognise some limitations. Although English primary care records are longitudinal and comprehensive, certain confounders were not captured. Shielding, as recommended for groups of clinically vulnerable people by the Chief Medical Officer,30 may have reduced the risk of infection, thus likely biasing results towards the null. In mediator adjusted models, we adjusted for concomitant use of oral glucocorticoids, however this adjustment is likely to be imperfect, leading to residual confounding. We also considered cardiovascular disease and diabetes mellitus to be mediators in the relationship between IMIDs and severe COVID-19 outcomes, however the timing of mediator assessment at index means they could have predated the IMID diagnosis, and hence not be true mediators. Assessment of glucocorticoid doses and usage is challenging due to absent precise dose information, prescribing of reducing dose regimens, low-cost medication administered in hospital alongside high-cost drugs, pandemic stockpiling and patient-led discontinuation due to COVID-19-related concerns (may also apply to immune modifying medication). Finally, as with any large study using routinely-collected data, there is a possibility of misclassification of exposure status. Misclassification of high cost drug exposure is unlikely, as high cost drug information is critical for billing. However, it is possible that we may have systematically excluded some individuals taking standard systemic drugs from analyses comparing targeted to standard systemic therapies. These exclusions may result in underestimation of risks in the standard systemic group due to differential exclusion of “new starters”, where the first prescription was in hospital. In practice, we expect the effects of this misclassification to be minimal due to the short time window. ### Implications and summary We have used one of the largest population-based datasets globally with linked data on immune modifying drugs to describe COVID-19 risks for people with IMIDs. We found that COVID-19 death and hospitalisation were higher in people with IMIDs; however, we saw no increased risk of adverse COVID-19 outcomes in those on most targeted immune modifying drugs for IMIDs compared to standard systemics. The availability of the “high-cost” novel data fills an important gap in available UK EHR data, and routine availability of these data is essential for future research. Our findings provide an evidence-base to inform policy on booster vaccination prioritisation and risk mitigating behaviour advice, but must be interpreted in the context of UK public health policy on shielding. Findings will support health care professionals engaging in shared decision making and communication of risk. ## Data Availability NHS England is the data controller; TPP is the data processor; and the key researchers on OpenSAFELY are acting on behalf of NHS England. This implementation of OpenSAFELY is hosted within the TPP environment which is accredited, the ISO 27001 information security standard and is NHS IG Toolkit compliant; patient data has been pseudonymised for analysis and linkage using industry standard cryptographic hashing techniques; all pseudonymised datasets transmitted for linkage onto OpenSAFELY are encrypted; access, the platform is via a virtual private network (VPN) connection, restricted, a small group of researchers; the researchers hold contracts with NHS England and only access the platform, initiate database queries and statistical models; all database activity is logged; only aggregate statistical outputs leave the platform environment following best practice for anonymisation of results such as statistical disclosure control for low cell counts. The OpenSAFELY research platform adheres, the obligations of the UK General Data Protection Regulation (GDPR) and the Data Protection Act 2018. In March 2020, the Secretary of State for Health and Social Care used powers under the UK Health Service (Control of Patient Information) Regulations 2002 (COPI), require organisations, process confidential patient information for the purposes of protecting public health, providing healthcare services, the public and monitoring and managing the COVID-19 outbreak and incidents of exposure; this sets aside the requirement for patient consent. Taken together, these provide the legal bases, link patient datasets on the OpenSAFELY platform. GP practices, from which the primary care data are obtained, are required, share relevant health information, support the public health response, the pandemic, and have been informed of the OpenSAFELY analytics platform. [https://github.com/opensafely/immunosuppressant-meds-research](https://github.com/opensafely/immunosuppressant-meds-research) ## ADMINISTRATIVE ### Acknowledgements We are very grateful for all the support received from the TPP Technical Operations team throughout this work, and for generous assistance from the information governance and database teams at NHS England / NHSX. Additionally the North East Commissioning Support Unit provided support on behalf of all Commissioning Support Unit to aggregate the high cost drugs data. This study uses electronic health records, data are provided by patients and collected by the National Health Service as part of their care and support. We are very grateful to Professor Joe West from the University of Nottingham and Professor Daniel Prieto-Alhmabra for comments on an early version of the protocol. ### Conflicts of Interest All authors have completed the ICMJE uniform disclosure form at [www.icmje.org/coi\_disclosure.pdf](http://www.icmje.org/coi_disclosure.pdf) and declare the following: BG has received research funding from Health Data Research UK (HDRUK), the Laura and John Arnold Foundation, the Wellcome Trust, the NIHR Oxford Biomedical Research Centre, the NHS National Institute for Health Research School of Primary Care Research, the Mohn-Westlake Foundation, the Good Thinking Foundation, the Health Foundation, and the World Health Organisation; he also receives personal income from speaking and writing for lay audiences on the misuse of science. IJD has received unrestricted research grants and holds shares in GlaxoSmithKline (GSK). CS received departmental research funding from AbbVie, Boehringer Ingelheim, Glaxo SmithKline, Leo, Pfizer, Novartis, Regeneron, SwedishOrphan Biovitrum AB, and Roche and is an investigator within consortia that have industry partners (see biomap.eu and psort.org.uk). JG has received honoraria from Abbvie, Amgen, Celgene, Chugai, Galapagos, Gilead, Janssen, Lilly, Novartis, Pfizer, Roche, Sobi, and UCB, and has research funding from Amgen, Aztra-Zeneca, Gilead, Janssen, Medicago, Novovax and Pfizer. MY has received honoraria from AbbVie and UCB. CWL has received honoraria from Abbvie, Celltrion, Ferring, Galapagos, Gilead, GSK, Iterative Scopes, Janssen, Pfizer and Takeda. ### Funding TPP and NECS CSU provided technical expertise and data infrastructure centre *pro bono* in the context of a national emergency. This work was supported by the Medical Research Council MR/V015737/1. BG’s work on better use of data in healthcare more broadly is currently funded in part by: NIHR Oxford Biomedical Research Centre, NIHR Applied Research Collaboration Oxford and Thames Valley, the Mohn-Westlake Foundation, NHS England, and the Health Foundation; all DataLab staff are supported by BG’s grants on this work. LS reports grants from Wellcome, MRC, NIHR, UKRI, British Council, GSK, British Heart Foundation, and Diabetes UK outside this work. JB is funded by a studentship from GSK. AS is employed by LSHTM on a fellowship sponsored by GSK. KB holds a Sir Henry Dale fellowship jointly funded by Wellcome and the Royal Society. HIM is funded by the National Institute for Health Research (NIHR) Health Protection Research Unit in Immunisation, a partnership between Public Health England and LSHTM. AYSW holds a fellowship from BHF. EW holds grants from MRC. ID holds grants from NIHR and GSK. RM holds a Sir Henry Wellcome Fellowship funded by the Wellcome Trust. HF holds a UKRI fellowship. RME is funded by HDR-UK and the MRC. SML was supported by a Wellcome Trust Senior Research Fellowship in Clinical Science (205039/Z/16/Z). The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the funders. SML was also supported by Health Data Research UK (Grant number: LOND1), which is funded by the UK Medical Research Council, Engineering and Physical Sciences Research Council, Economic and Social Research Council, Department of Health and Social Care (England), Chief Scientist Office of the Scottish Government Health and Social Care Directorates, Health and Social Care Research and Development Division (Welsh Government), Public Health Agency (Northern Ireland), British Heart Foundation and Wellcome Trust. SML is an investigator on the European Union Horizon 2020-funded BIOMAP Consortium ([http://www.biomap-imi.eu/](http://www.biomap-imi.eu/)). CS acknowledges support for this research from the NIHR Biomedical Research Centre (BRC) at King’s College London and Guy’s and St Thomas’ NHS Foundation Trust and the Psoriasis Association. CWL is funded by a UKRI Future Leaders Fellowship. The views expressed are those of the authors and not necessarily those of the NIHR, NHS England, Public Health England or the Department of Health and Social Care. Funders had no role in the study design, collection, analysis, and interpretation of data; in the writing of the report; and in the decision, submit the article for publication. This research was funded in whole or in part by the Wellcome Trust [G205039/Z/16/Z]. For the purpose of Open Access, the author has applied a CC BY public copyright licence to any Author Accepted Manuscript (AAM) version arising from this submission ### Information governance and ethical approval NHS England is the data controller; *TPP is the data processor*; and the key researchers on OpenSAFELY are acting on behalf of NHS England. This implementation of OpenSAFELY is hosted within the *TPP environment which is* accredited, the ISO 27001 information security standard and *is* NHS IG Toolkit compliant;31,32 patient data has been pseudonymised for analysis and linkage using industry standard cryptographic hashing techniques; all pseudonymised datasets transmitted for linkage onto OpenSAFELY are encrypted; access, the platform is via a virtual private network (VPN) connection, restricted, a small group of researchers; the researchers hold contracts with NHS England and only access the platform, initiate database queries and statistical models; all database activity is logged; only aggregate statistical outputs leave the platform environment following best practice for anonymisation of results such as statistical disclosure control for low cell counts.33 The OpenSAFELY research platform adheres, the obligations of the UK General Data Protection Regulation (GDPR) and the Data Protection Act 2018. In March 2020, the Secretary of State for Health and Social Care used powers under the UK Health Service (Control of Patient Information) Regulations 2002 (COPI), require organisations, process confidential patient information for the purposes of protecting public health, providing healthcare services, the public and monitoring and managing the COVID-19 outbreak and incidents of exposure; this sets aside the requirement for patient consent.34 Taken together, these provide the legal bases, link patient datasets on the OpenSAFELY platform. GP practices, from which the primary care data are obtained, are required, share relevant health information, support the public health response, the pandemic, and have been informed of the OpenSAFELY analytics platform. This study was approved by the Health Research Authority (REC Reference 20/LO/0651) and by the LSHTM Ethics Board (Reference 21863). ### Guarantor BG/LS is guarantor. ### Contributions BG conceived the OpenSAFELY platform and the approach. LS and BG led the project overall and are guarantors. SML, CS, NAK, CL, JG, KEM, BMK, SJWE conceptualised the study. SB led on software development. AM led on information governance. CJB, CEM, DH, RC, GH, TW, SCJB, PI, JC, DE, JP, and SH curated data; SML, CS, NAK, CL, JG, KEM, BMK, KB, CTR, BM, CJB, CEM, IJD, AJW, HIM, JC, HJF, HJC, JT, RME, LAT, RME, AYSW and JP conceptualised disease categories and code lists; MY, SN, JG, KSB, JB and JM wrote statistical analysis code; KEM, KSB, JG, SN, MY, JB did data visualisation; EJW, HJC, LS, and BG obtained ethical approvals; CJB, CEM, SCJB, SD, AG, LF, PI, AJW, JC, DE, WJH, and FH contributed to software development; and AS, SML, KEM, DH, CS, KSB, JG, SN, MY, NAK, JM, LAT, RME, AYSW and SJWE reviewed and edited the manuscript. All authors were involved in design and conceptual development and reviewed and approved the final manuscript. ## SUPPLEMENTARY MATERIAL ### Text S1. High-cost drugs dataset To achieve a comprehensive national medicines dataset, we arranged for the unique collation of a single national hospital medication dataset for “high-cost drugs”. High-cost drugs are typically newer and more expensive medications (e.g., adalimumab) that are supplied by hospitals in England for ongoing disease management. High-cost drugs are paid for by commissioners (in England, National Health Service Clinical Commission Groups commission most secondary care services and play a part in the commissioning of GP services) and as part of this process, unique patient identifiers are supplied by hospitals, the responsible commissioner (there are approximately 200 commissioners). The North East Commissioning Support Unit collected high-cost drug data from commissioners and compiled it into a single national dataset (including medication name, date of medication supply and a unique patient identifier – NHS number – allowing linkage with other healthcare data). To our knowledge this is the first time in England that a comprehensive national medicines dataset including high-cost drug data has been available at individual patient level. More information is available in our accompanying short data report.[link will be provided in published paper] Drug exposure was defined by at least one prescription, or delivery of medication to an individual before March 1st, 2020 (date chosen as some of medications were either specifically used, or stopped, due to the pandemic). For each individual, we defined drug exposure based on the closest drug recorded prior to the study start (March 1st 2020) allowing for a maximum of 6 months before the start of the study for all agents apart from rituximab, where we permitted a 12-month exposure window (to allow for the frequency of treatment and the longer duration of response of rituximab, including known prolonged effects on vaccine responses).16,17 We could not evaluate medication switching or adherence during the follow-up period, as data was only available up to March 2020. ### Text S2. Covariates All analytical code for all covariates definitions are openly available online for inspection and re-use. Some brief descriptions are provided below. #### Age Age was defined by the age reached by the 1st March 2020 (study start date). We categorised age groups as: 18–39, 40–49, 50–59, 60–69, 70–79 and 80+ years. #### Smoking and body mass index We categorised smoking status as current-, former- or never-smokers, based on primary care morbidity coding recorded prior to March 1st, 2020. People with missing smoking status were categorised as never smokers. Body mass index (BMI) was calculated from height and weight measures recorded in primary care, based on weight measurements within the last 10 years (before March 1st, 2020), and restricted to records from after age 16 years. BMI was categorised according to World Health Organisation (WHO) BMI classification: no evidence of obesity (BMI < 30); obese class I (BMI 30–34.9); obese class II (BMI 35–39.9); and obese class III (BMI 40+). Individuals with missing BMI data were categorised as having no evidence of obesity (i.e., BMI < 30). We undertook a sensitivity analysis without re-categorising missing BMI and smoking status data (i.e., complete-case analysis of only those with complete smoking and BMI data). Our approach reflects an awareness that BMI and smoking data are likely to be not missing at random, and has precedent in previous OpenSAFELY analyses.1,2 #### Current glucocorticoid use We defined current glucocorticoid use based on a primary care prescription of oral prednisolone at any dose in the three months before study start. #### Diabetes mellitus We defined diabetes mellitus based on the most recent glycosylated haemoglobin (Hba1c) measurement recorded in primary care at any time before study start (March 1st, 2020): controlled diabetes Hba1c<58mmol/mol; uncontrolled diabetes Hba1c≥58mmol/mol; or diabetes (based on morbidity coded diabetes) with no recorded Hba1c measure. #### End stage renal failure We defined end stage renal failure based on estimated glomerular filtration rate calculated from the most recent serum creatinine result recorded 15 months before study start (March 1st, 2020) in primary care data. We only adjusted for end stage renal failure in sensitivity analyses. #### Ethnicity We defined ethnicity based on coding in primary care. We classified ethnicity as: White, Black, South Asian, and mixed or other. We only adjusted for ethnicity in sensitivity analyses due to the large proportion of missing data (**Tables 1 and 2**). #### Other chronic comorbidities We defined cardiovascular disease (chronic heart failure, ischaemic heart disease, and severe valve or congenital heart disease likely to require lifelong follow-up), cancer (excluding non-melanoma skin cancer), and stroke based on any morbidity coding in primary care recorded prior to March 1st, 2020 (study start date). We also adjusted for chronic respiratory disease (chronic obstructive pulmonary disease, fibrosing lung disease, bronchiectasis or cystic fibrosis) and chronic liver disease in a sensitivity analysis, again defined based on any morbidity coding in primary care recorded prior to study start. **References** 1. Bhaskaran K, Bacon S, Evans SJ, et al. Factors associated with deaths due to COVID-19 versus other causes: population- based cohort analysis of UK primary care data and linked national death registrations within the OpenSAFELY platform. Lancet Reg Health Eur 2021;6:100109. 2. Bhaskaran K, Rentsch CT, MacKenna B, et al. HIV infection and COVID-19 death: a population-based cohort analysis of UK primary care data and linked national death registrations within the OpenSAFELY platform. Lancet HIV 2021;8:e24-e32 ![Figure S1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/05/06/2021.09.03.21262888/F4.medium.gif) [Figure S1.](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/F4) Figure S1. Conceptual framework: risk of COVID-19-related death in people with IMIDs compared to the general population The conceptual framework represents the assumed associations between covariates and primary exposure and outcome. Pink circles represent ancestors of the exposure and outcome, blue circles represent ancestors of the outcome, pink lines represent biasing paths (i.e. confounding) and green lines represent causal paths.. The minimally sufficient adjustment set (i.e. the covariates adjusted for in confounder adjusted models) represents covariates such that the adjustment for this set of variables will minimize confounding bias when estimating the association between the exposure and the outcome. ![Figure S2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/05/06/2021.09.03.21262888/F5.medium.gif) [Figure S2.](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/F5) Figure S2. Conceptual framework: risk of COVID-19-related death in people with IMIDS taking targeted immunosuppressive drugs compared to those prescribed standard systemic drugs. The conceptual framework A directed acyclic graph represents the assumed associations between covariates and primary exposure and outcome. Pink circles represent ancestors of the exposure and outcome, blue circles representancestors of the outcome, pink lines represent biasing paths (i.e. confounding) and green lines represent causal paths. The minimally sufficient adjustment set (i.e. the covariates adjusted for in confounder adjusted models) represents covariates such that the adjustment for this set of variables will minimize confounding bias when estimating the association between the exposure and the outcome. View this table: [Table S1.](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T3) Table S1. Sensitivity analyses ![Figure S3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/05/06/2021.09.03.21262888/F6.medium.gif) [Figure S3.](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/F6) Figure S3. Flow diagram for identification of study population View this table: [Table S2](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T4) Table S2 **Main analysis:** Hazard Ratios (HRs) and 95% confidence intervals (CI) for COVID-19-related death, death/critical care admission or hospitalisation, in people with **IMIDs compared to the general population** View this table: [Table S3](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T5) Table S3 **Sensitivity analyses**: **COVID-19-related death** in **IMIDs** compared to the general population. View this table: [Table S4](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T6) Table S4 **Sensitivity analyses**: **COVID-19-related critical care admission/death** in **IMIDs** compared tothe general population. View this table: [Table S5](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T7) Table S5 Number of individuals with each IMID type receiving standard systemic or any targeted immune modifying therapy View this table: [Table S6](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T8) Table S6 **Main analysis:** Hazard Ratios (HRs) and 95% confidence intervals (CI) for COVID-19-related death, death/critical care admission or hospitalisation, in people with **IMIDs on targeted immune modifying therapies** compared to those with **IMIDs on standard systemic immune modifying therapies**. View this table: [Table S7](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T9) Table S7 **Sensitivity analyses:** for **COVID-19-related death** in people with IMIDs on **targeted immune modifying therapies** compared tothose with IMIDs on standard systemic immune modifying therapies. View this table: [Table S8](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T10) Table S8 **Sensitivity analyses**: **COVID-19-related critical care admission/death** in people with IMIDs on **targeted immune modifying therapies** compared to those with IMIDs on standard systemicimmune modifying therapies View this table: [Table S9](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T11) Table S9 **Sensitivity analysis:** Hazard Ratios (95% CI) for COVID-19-related death with **TNF inhibition monotherapy** compared to those receiving combination TNF inhibition and standard systemic therapy. View this table: [Table S10](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T12) Table S10 Descriptive characteristics of population on targeted immunosuppression with **infliximab** View this table: [Table S11](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T13) Table S11 **Sensitivity analysis:** Hazard Ratios (95% CI) for **COVID-19-related death** in people with anIMID on **infliximab therapy** compared to those on standard systemic therapy View this table: [Table S12](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T14) Table S12 **Analysis adjusting for age using a four-knot cubic spline:** Hazard Ratios (HRs) and 95% confidence intervals (CI) for COVID-19-related death, death/critical care admission or hospitalisation, in people with **IMIDs compared to the general population** View this table: [Table S13](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T15) Table S13 **Analysis adjusting for age using a four-knot cubic spline:** Hazard Ratios (HRs) and 95% confidence intervals (CI) for COVID-19-related death, death/critical care admission or hospitalisation, in people with **IMIDs on targeted immune modifying therapies** compared tothose with **IMIDs on standard systemic immune modifying therapies**. ### Text S3. Quantitative bias analysis #### Methods We considered severity of immune-mediated inflammatory disease and degree of shielding to be potential confounders which are not available in the data. As such we conducted quantitative bias analysis using E-values to assess how strongly associated these unmeasured confounders would need to be with exposure and outcome to potentially fully explain observed non-null associations (i.e., the association adjusted for both measured covariates and the unmeasured confounder would be null).21 In order to apply these methods to the hazard ratio, rather than the risk ratio, we make the assumption that the outcomes are rare. #### Results A non-null association was observed between rituximab and COVID-19 death, critical care admission/death, and hospitalisation. For an unmeasured confounder to potentially fully explain the point estimate for the observed confounder-adjusted association the unmeasured confounder would need to be associated with one of exposure or outcome, conditional on measured covariates, by at least risk ratios 2.75, 3.25 and 2.56 for COVID-19 death, critical care admission/death, and hospitalisation respectively (1.46, 1.95, and 1.59 to potentially explain the lower bound of the 95% CI) (**Figures S4-S6, Tables S14-S17**). Furthermore, to potentially fully explain the point estimate the unmeasured confounder would need to be associated with both exposure and outcome, conditional on measured covariates, by at least risk ratio 1.68, 1.92 and 1.59 for COVID-19 death, critical care admission/death and hospitalisation (1.11, 1.31, and 1.16 to potentially explain the lower bound of the 95% CI). Additionally, a non-null association was observed between JAK-inhibitors and COVID-19 hospitalisation. To potentially fully explain the confounder-adjusted point estimate (lower bound of 95% CI) the unmeasured confounder would need to be associated with one of exposure or outcome, conditional on measured covariates, by at least risk ratio 3.02 (1.40) and both of exposure and outcome by 1.81 (1.09) (**Figure S7, Table S17).** #### Discussion An unmeasured confounder moderately associated with both exposure and outcome could potentially explain the associations of rituximab and JAK inhibitors with COVID-19 outcomes. However, whether this is the case in this study and the extent of true residual confounding depends on both the prevalence of and strength of associations with exposure and outcome of the unmeasured confounders, shielding and severity of IMID, which are not known. ![Figure S4.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/05/06/2021.09.03.21262888/F7.medium.gif) [Figure S4.](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/F7) Figure S4. Minimum associations with unmeasured confounder required to potentially explainobserved confounder-adjusted association between rituximab and **COVID-19 death** ![Figure S5.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/05/06/2021.09.03.21262888/F8.medium.gif) [Figure S5.](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/F8) Figure S5. Minimum associations with unmeasured confounder required to potentially explainobserved confounder-adjusted association between rituximab and COVID-19 **critical care admission/death** ![Figure S6.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/05/06/2021.09.03.21262888/F9.medium.gif) [Figure S6.](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/F9) Figure S6. Minimum associations with unmeasured confounder required to potentially explainobserved confounder-adjusted association between rituximab and **COVID-19 hospitalisation** ![Figure S7.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/05/06/2021.09.03.21262888/F10.medium.gif) [Figure S7.](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/F10) Figure S7. Minimum associations with unmeasured confounder required to potentially explainobserved confounder-adjusted association between JAK inhibitors and **COVID-19 hospitalisation** View this table: [Table S14](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T16) Table S14 Minimum possible bias-adjusted association between rituximab and COVID-19 death View this table: [Table S15](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T17) Table S15 Minimum possible bias-adjusted association between rituximab and COVID-19 critical care admission/death View this table: [Table S16](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T18) Table S16 Minimum possible bias-adjusted association between rituximab and COVID-19 hospitalisation View this table: [Table S17](http://medrxiv.org/content/early/2022/05/06/2021.09.03.21262888/T19) Table S17 Minimum possible bias-adjusted association between JAK inhibitors and COVID-19 hospitalisation ## Footnotes * * Joint first authors * ** Joint last authors * This version of the manuscript has been revised to update figures 2 and 3. * Received September 3, 2021. * Revision received May 5, 2022. * Accepted May 6, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## REFERENCES 1. Wu Z, McGoogan JM. Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China: summary of a report of 72 314 cases from the Chinese Center for Disease Control and Prevention. JAMA 2020; 323:1239–42. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2020.2648&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F06%2F2021.09.03.21262888.atom) 2. Williamson EJ, Walker AJ, Bhaskaran K, et al. Factors associated with COVID-19-related death using OpenSAFELY. Nature. 2020; 584: 430–6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2521-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F06%2F2021.09.03.21262888.atom) 3. Han E, Tan MMJ, Turk E, et al. Lessons learnt from easing COVID-19 restrictions: an analysis of countries and regions in Asia Pacific and Europe. Lancet 2020; 396: 1525–34. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F06%2F2021.09.03.21262888.atom) 4. NHS Digital. Data Security and Protection Toolkit. [https://digital.nhs.uk/data-and-information/looking-after-information/data-security-and-information-governance/data-security-and-protection-toolkit](https://digital.nhs.uk/data-and-information/looking-after-information/data-security-and-information-governance/data-security-and-protection-toolkit) (accessed Aug 26, 2021). 5. Singh JA, Wells GA, Christensen R, et al. Adverse effects of biologics: a network meta-analysis and Cochrane overview. Cochrane Database Syst Rev 2011; : CD008794. 6. Galloway JB, Hyrich KL, Mercer LK, et al. Anti-TNF therapy is associated with an increased risk of serious infections in patients with rheumatoid arthritis especially in the first 6 months of treatment: updated results from the British Society for Rheumatology Biologics Register with special emphasis on risks in the elderly. Rheumatology 2011; 50: 124–31. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/rheumatology/keq242&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20675706&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F06%2F2021.09.03.21262888.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000285193500017&link_type=ISI) 7. Brenner EJ, Ungaro RC, Gearry RB, et al. Corticosteroids, But Not TNF Antagonists, Are Associated With Adverse COVID-19 Outcomes in Patients With Inflammatory Bowel Diseases: Results From an International Registry. Gastroenterology 2020; 159:481–91.e3. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/j.gastro.2020.05.032&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F06%2F2021.09.03.21262888.atom) 8. Gianfrancesco M, Hyrich KL, Al-Adely S, et al. Characteristics associated with hospitalisation for COVID-19 in people with rheumatic disease: data from the COVID-19 Global Rheumatology Alliance physician-reported registry. Annals of the Rheumatic Diseases. 2020; 79: 859–66. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6ImFubnJoZXVtZGlzIjtzOjU6InJlc2lkIjtzOjE0OiI3OS9TdXBwbF8xLzg1OSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA1LzA2LzIwMjEuMDkuMDMuMjEyNjI4ODguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 9. Mahil SK, Dand N, Mason KJ, et al. Factors associated with adverse COVID-19 outcomes in patients with psoriasis—insights from a global registry–based study. J Allergy Clin Immunol 2021; 147: 60–71. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaci.2020.10.007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F06%2F2021.09.03.21262888.atom) 10. Psoriasis: assessment and management. Clinical guideline [CG153]. 2012; published online Oct 24. [https://www.nice.org.uk/guidance/cg153/chapter/1-Recommendations](https://www.nice.org.uk/guidance/cg153/chapter/1-Recommendations) (accessed Aug 26, 2021). 11. Rheumatoid arthritis in adults: management. NICE guideline [NG100]. 2018; published online July 11. [https://www.nice.org.uk/guidance/ng100/chapter/Recommendations](https://www.nice.org.uk/guidance/ng100/chapter/Recommendations) (accessed Aug 26, 2021). 12. Spondyloarthritis in over 16s: diagnosis and management. NICE guideline [NG65]. 2017; published online Feb 28. [https://www.nice.org.uk/guidance/NG65/chapter/Recommendations](https://www.nice.org.uk/guidance/NG65/chapter/Recommendations) (accessed Aug 26, 2021). 13. Crohn’s disease: management. NICE guideline [NG129]. 2019; published online May 3. [https://www.nice.org.uk/guidance/ng129/chapter/Recommendations](https://www.nice.org.uk/guidance/ng129/chapter/Recommendations) (accessed Aug 26, 2021). 14. Ulcerative colitis: management. NICE guideline [NG130]. 2019; published online May 3. [https://www.nice.org.uk/guidance/ng130/chapter/Recommendations](https://www.nice.org.uk/guidance/ng130/chapter/Recommendations) (accessed Aug 26, 2021). 15. Adalimumab for treating moderate to severe hidradenitis suppurativa. Technology appraisal guidance [TA392]. 2016; published online June 22. [https://www.nice.org.uk/guidance/TA392/chapter/1-Recommendations](https://www.nice.org.uk/guidance/TA392/chapter/1-Recommendations) (accessed Aug 26, 2021). 16. Buch MH, Smolen JS, Betteridge N, et al. Updated consensus statement on the use of rituximab in patients with rheumatoid arthritis. Ann Rheum Dis 2011; 70: 909–20. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6ImFubnJoZXVtZGlzIjtzOjU6InJlc2lkIjtzOjg6IjcwLzYvOTA5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDUvMDYvMjAyMS4wOS4wMy4yMTI2Mjg4OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 17. van Assen S, Holvast A, Benne CA, et al. Humoral responses after influenza vaccination are severely reduced in patients with rheumatoid arthritis treated with rituximab. Arthritis Rheum 2010; 62: 75–81. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/art.25033&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20039396&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F06%2F2021.09.03.21262888.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000279275700009&link_type=ISI) 18. Richards-Belle A, Orzechowska I, Gould DW, et al. COVID-19 in critical care: epidemiology of the first epidemic wave across England, Wales and Northern Ireland. Intensive Care Med 2020; 46: 2035–47. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00134-020-06267-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F06%2F2021.09.03.21262888.atom) 19. Ministry of Housing, Communities & Local Government. English indices of deprivation 2015. 2015; published online Sept 30. [https://www.gov.uk/government/statistics/english-indices-of-deprivation-2015](https://www.gov.uk/government/statistics/english-indices-of-deprivation-2015) (accessed Aug 26, 2021). 20. Noble M, Wright G, Smith G, Dibben C. Measuring Multiple Deprivation at the Small-Area Level. Environ Plan A 2006; 38: 169–85. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1068/a37168&link_type=DOI) 21. Ding P, VanderWeele TJ. Sensitivity Analysis Without Assumptions. Epidemiology 2016; 27: 368–77. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/EDE.0000000000000457&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26841057&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F06%2F2021.09.03.21262888.atom) 22. Motta F, Sica A, Selmi C. Frailty in Rheumatic Diseases. Front Immunol 2020; 11: 576134. 23. Eder L, Croxford R, Drucker A, et al. OP0285 COVID-19 HOSPITALIZATIONS, ICU ADMISSION, AND DEATH AMONG PATIENTS WITH IMMUNE MEDIATED INFLAMMATORY DISEASES (IMID) – A POPULATION-BASED STUDY. Ann Rheum Dis 2021; 80: 173–173. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6ImFubnJoZXVtZGlzIjtzOjU6InJlc2lkIjtzOjE0OiI4MC9TdXBwbF8xLzE3MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA1LzA2LzIwMjEuMDkuMDMuMjEyNjI4ODguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 24. Attauabi M, Seidelin JB, Felding OK, et al. Coronavirus disease 2019, immune-mediated inflammatory diseases and immunosuppressive therapies - A Danish population-based cohort study. J Autoimmun 2021; 118: 102613. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F06%2F2021.09.03.21262888.atom) 25. Ludvigsson JF, Axelrad J, Halfvarson J, et al. Inflammatory bowel disease and risk of severe COVID-19: A nationwide population-based cohort study in Sweden. United European Gastroenterol J 2021; 9: 177–92. 26. Akiyama S, Hamdeh S, Micic D, Sakuraba A. Prevalence and clinical outcomes of COVID-19 in patients with autoimmune diseases: a systematic review and meta-analysis. Ann Rheum Dis 2020; published online Oct 13. DOI:10.1136/annrheumdis-2020-218946. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6ImFubnJoZXVtZGlzIjtzOjU6InJlc2lkIjtzOjg6IjgwLzMvMzg0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDUvMDYvMjAyMS4wOS4wMy4yMTI2Mjg4OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 27. Lamure S, Dulery R, Delord M, et al. Abstract S09-02: High incidence of persistent COVID-19 among patients with lymphoma treated with B-cell depleting immunotherapy. Clin Cancer Res 2021; 27: S09–02. 28. Investigators R-C. Interleukin-6 Receptor Antagonists in Critically Ill Patients with Covid-19. N Engl J Med 2021; 384: 1491–502. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F06%2F2021.09.03.21262888.atom) 29. Abani O, Abbas A, Abbas F, et al. Tocilizumab in patients admitted to hospital with COVID-19 (RECOVERY): a randomised, controlled, open-label, platform trial. Lancet 2021; 397: 1637–45. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F05%2F06%2F2021.09.03.21262888.atom) 30. Guidance on shielding and protecting people who are clinically extremely vulnerable from COVID-19. [https://www.gov.uk/government/publications/guidance-on-shielding-and-protecting-extremely-vulnerable-persons-from-covid-19/guidance-on-shielding-and-protecting-extremely-vulnerable-persons-from-covid-19](https://www.gov.uk/government/publications/guidance-on-shielding-and-protecting-extremely-vulnerable-persons-from-covid-19/guidance-on-shielding-and-protecting-extremely-vulnerable-persons-from-covid-19) (accessed May 18, 2021). 31. BETA – Data Security Standards - NHS Digital. NHS Digital. [https://digital.nhs.uk/about-nhs-digital/our-work/nhs-digital-data-and-technology-standards/framework/beta\---|data-security-standards](https://digital.nhs.uk/about-nhs-digital/our-work/nhs-digital-data-and-technology-standards/framework/beta\---|data-security-standards) (accessed April 30, 2020). 32. Data Security and Protection Toolkit - NHS Digital. NHS Digital. [https://digital.nhs.uk/data-and-information/looking-after-information/data-security-and-information-governance/data-security-and-protection-toolkit](https://digital.nhs.uk/data-and-information/looking-after-information/data-security-and-information-governance/data-security-and-protection-toolkit) (accessed April 30, 2020). 33. ISB1523: Anonymisation Standard for Publishing Health and Social Care Data - NHS Digital. NHS Digital. [https://digital.nhs.uk/data-and-information/information-standards/information-standards-and-data-collections-including-extractions/publications-and-notifications/standards-and-collections/isb1523-anonymisation-standard-for-publishing-health-and-social-care-data](https://digital.nhs.uk/data-and-information/information-standards/information-standards-and-data-collections-including-extractions/publications-and-notifications/standards-and-collections/isb1523-anonymisation-standard-for-publishing-health-and-social-care-data) (accessed April 30, 2020). 34. Secretary of State for Health and Social Care - UK Government. Coronavirus (COVID-19): notification to organisations to share information. 2020; published online April 1. [https://web.archive.org/web/20200421171727/](https://web.archive.org/web/20200421171727/)[https://www.gov.uk/government/publications/coronavirus-covid-19-notification-of-data-controllers-to-share-information](https://www.gov.uk/government/publications/coronavirus-covid-19-notification-of-data-controllers-to-share-information).