Inhaled corticosteroid use and risk COVID-19 related death among 966,461 patients with COPD or asthma: an OpenSAFELY analysis ============================================================================================================================= * The OpenSAFELY Collaborative * Anna Schultze * Alex J Walker * Brian MacKenna * Caroline E Morton * Krishnan Bhaskaran * Jeremy P Brown * Christopher T Rentsch * Elizabeth Williamson * Henry Drysdale * Richard Croker * Seb Bacon * William Hulme * Chris Bates * Helen J Curtis * Amir Mehrkar * David Evans * Peter Inglesby * Jonathan Cockburn * Helen I McDonald * Laurie Tomlinson * Rohini Mathur * Kevin Wing * Angel YS Wong * Harriet Forbes * John Parry * Frank Hester * Sam Harper * Stephen JW Evans * Jennifer Quint * Liam Smeeth * Ian J Douglas * Ben Goldacre ## Abstract **Background** Early descriptions of the coronavirus outbreak showed a lower prevalence of asthma and COPD than was expected for people diagnosed with COVID-19, leading to speculation that inhaled corticosteroids (ICS) may protect against infection with SARS-CoV-2, and development of serious sequelae. We evaluated the association between ICS and COVID-19 related death using linked electronic health records in the UK. **Methods** We conducted cohort studies on two groups of people (COPD and asthma) using the OpenSAFELY platform to analyse data from primary care practices linked to national death registrations. People receiving an ICS were compared to those receiving alternative respiratory medications. Our primary outcome was COVID-19 related death. **Findings** We identified 148,588 people with COPD and 817,973 people with asthma receiving relevant respiratory medications in the four months prior to 01 March 2020. People with COPD receiving ICS were at a greater risk of COVID-19 related death compared to those receiving a long-acting beta agonist (LABA) and a long-acting muscarinic antagonist (LAMA) (adjusted HR = 1.38, 95% CI = 1.08 – 1.75). People with asthma receiving high dose ICS were at an increased risk of death compared to those receiving a short-acting beta agonist (SABA) only (adjusted HR = 1.52, 95%CI = 1.08 – 2.14); the adjusted HR for those receiving low-medium dose ICS was 1.10 (95% CI = 0.82 – 1.49). Quantitative bias analyses indicated that an unmeasured confounder of only moderate strength of association with exposure and outcome could explain the observed associations in both populations. **Interpretation** These results do not support a major role of ICS in protecting against COVID-19 related deaths. Observed increased risks of COVID-19 related death among people with COPD and asthma receiving ICS can be plausibly explained by unmeasured confounding due to disease severity. **Funding** This work was supported by the Medical Research Council MR/V015737/1. Keywords * COVID-19 * Pulmonary Disease * Chronic Obstructive * Asthma * Steroids * Pharmacoepidemiology ## Introduction The ongoing pandemic due to the novel coronavirus, SARS-CoV-2, has now affected over 7 million people worldwide with at least 400,000 people having died with COVID-191. People with more severe COVID-19 outcomes, including hospitalisation or death, tend to be older and have pre-existing comorbidities 2,3,4,5,6,7,8. Severe outcomes are often a result of lung complications, such as acute respiratory distress syndrome (ARDS) and respiratory failure. However, early reports of COVID-19 patients described an unexpectedly low prevalence of chronic respiratory conditions among hospitalised patients9. Although other studies suggest that chronic lung diseases, including chronic obstructive pulmonary disease (COPD), increase the risk of severe outcomes6–8, reported effect sizes for asthma have been modest 6,7. This has led to speculation that treatments for respiratory disease, specifically inhaled corticosteroids (ICS), may have a protective effect against SARS-CoV-2.9–11. ICS are used to reduce airway inflammation, oedema, and mucus secretions. In-vitro evidence indicates that the ICS ciclesonide can suppress SARS-CoV-2 replication12, and budesonide combined with glycopyrronium and formoterol inhibits the production of cytokines in cells exposed to HCoV-229E, another human coronavirus13. The oral/IV steroid dexamethasone has recently been shown to reduce the risk of death in severe COVID-1914. Conversely, although ICS have low systemic absorption, they have been associated with increased risk of developing pneumonia in people with COPD 15–17, as well as other systemic steroid-related adverse effects18. A recent systematic review of the role of ICS in SARS-CoV-2, SARS-CoV-1, and MERS found no studies investigating the impact of prior ICS use on outcomes in any of these infections 10. We therefore set out to explore the association between current ICS use and outcomes in COVID-19, using the OpenSAFELY platform which contains linked primary care electronic health record data for approximately 40% of the population in England. ## Methods ### Study Design We conducted two cohort studies using primary care electronic health record (EHR) data linked to death data from the Office for National Statistics (ONS). The index date (start of follow up) for both cohorts was 01 Mar 2020; follow-up lasted until 06 May 2020. ### Data Source Primary care records managed by the GP software provider The Phoenix Partnership (TPP) were linked to ONS death data through OpenSAFELY, a data analytics platform created by our team on behalf of NHS England19 to address urgent COVID-19 research questions 7 ([https://opensafely.org](https://opensafely.org)). OpenSAFELY provides a secure software interface allowing the analysis of pseudonymised primary care patient records from England in near real-time within the EHR vendor’s highly secure data centre, avoiding the need for large volumes of potentially disclosive pseudonymised patient data to be transferred off-site. This, in addition to other technical and organisational controls, minimises the risk of re-identification. Similarly pseudonymised datasets from other data providers are securely provided to the EHR vendor and linked to the primary care data. The dataset analysed within OpenSAFELY is based on 24 million people currently registered with GP surgeries using TPP SystmOne software. It includes pseudonymised data such as coded diagnoses, medications and physiological parameters. No free text data is included. ### Study Populations The COPD cohort included adults older than 35 years with COPD and current or former smoking recorded any time before the index date20. We excluded people with prior diagnoses of other chronic respiratory conditions, or with asthma in the three years before the index date21, and those receiving nebulised COPD medications in the twelve months before the index date or a leukotriene receptor antagonist (indicating potential asthma) in the four months before the index date. The asthma cohort included adults older than 18 years with asthma recorded within three years prior to the index date. People with COPD or other chronic respiratory conditions prior to the index date were excluded, as were those receiving a LAMA without an ICS, as this indicates possible COPD22. People with missing data for gender, index of multiple deprivation (IMD), or less than one year of primary care records were excluded (supplemental figure 1). ### Exposures In the COPD population, people issued at least one ICS prescription within four months prior to the index date either in combination with LABA or LAMA/LABA, or as single therapy provided there was also at least one prescription record of a LABA, were compared with those with a prescription for a LABA/LAMA (combined or as separate single therapy prescriptions) only22. We did not include patients receiving LAMA monotherapy, as we were expecting greater clinical comparability between the LAMA/LABA and ICS-based therapy groups. In the asthma population, people prescribed high dose ICS and low/medium-dose ICS during the four months before index date were compared with those prescribed SABA only. Exposure for people prescribed both high and low/medium dose ICS was assigned according to their most recent prescription. Inhalers were assigned to low/medium or high dose based the [OpenPrescribing.net](http://OpenPrescribing.net) prescribing explorer which was developed based on BTS/SIGN guidance23. Studies have shown that a significant percentage of people with asthma receiving SABA only are eligible for ICS treatment24, suggesting they have similar disease severity to those receiving ICS and therefore represent a reasonable comparator group. The characteristics of all other people are described in supplementary material (supplemental table 1-2), however they are excluded from regression models to avoid comparisons to individuals not prescribed drugs of interest 25. ### Outcomes The outcome was COVID-19 related death as registered in ONS data using ICD-10 codes U07.1 (“COVID-19, virus identified”) and U07.2 (“COVID-19, virus not identified”) listed either as the underlying or any contributing cause of death. The latter ICD-10 code is used when laboratory testing is inconclusive or unavailable26. ### Covariates Potential determinants of exposures and outcomes were identified by reviewing literature and through discussions with practising clinicians. As this is a study of current users, determinants of exposures include both factors that may affect the initial choice of treatment as well as those that influence whether patients remain on a certain treatment. The final list of potential confounders can be seen in box 1. Our methodology for creating codelists associated with these confounders has been previously described7: this included clinical and epidemiological review and sign-off by at least two authors. Detailed information on every codelist is shared at [https://codelists.opensafely.org/](https://codelists.opensafely.org/). ### Statistical Methods Individuals characteristics were summarised using descriptive statistics, stratified by exposure status. Time to the primary outcome is displayed in Kaplan-Meier (KM) plots with time in study as the timescale. The competing risk of death from non-COVID-19 causes was dealt with by analysing the cause-specific hazard, with people dying from other causes censored at their date of death27. We used cause-specific Cox regression models to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for the association between exposure categories and the outcome in each population. Univariable models, models adjusted for age (using restricted cubic splines) and sex as well as fully adjusted models including all covariates were fitted. Region was included as a stratification variable in fully adjusted models. We evaluated an *a priori* specified interaction between ICS exposure and age, to see if we could distinguish a differential effect in groups known to be at higher risk. ### Sensitivity Analyses In sensitivity analyses, first, we split the exposure categories in the COPD population to examine the effect of ICS with LABA/LAMA (triple combination) and ICS with LABA (dual combination) separately, anticipating greater underlying disease severity in people receiving triple therapy. Second, we restricted analyses to the largest ethnicity group (i.e., white British) to exclude any substantial confounding by ethnicity. We did not adjust for ethnicity in the main models as this was not anticipated to be a strong confounder and due to a sizable proportion of individuals with missing ethnicity (∼25%). In the asthma population, we varied the sample definition to include people with asthma diagnosed at any time, and a recent prescription for any asthma medication. ### Negative Control Outcomes We hypothesised that disease severity, but not ICS use, may influence the risk of non-COVID-19 related death. Analyses were therefore conducted using non-COVID-19 death as a negative control outcome censoring people at time of COVID-19 related death. If any potentially harmful association observed in primary analyses was due to confounding (i.e. people prescribed ICS had more severe underlying respiratory disease than those who did not) we expected to observe a similar association with non-COVID-19 related death in people prescribed ICS. ### Quantitative Bias Analysis We used e-value formulae to calculate the minimum strengths of association between an unmeasured confounder and exposure or outcome, conditional on measured covariates, necessary to fully explain observed associations28. ### Software and Reproducibility Data management was performed using Python 3.8 and SQL, with analysis carried out using Stata 16.1. All of the code used for data management and analyses is openly shared online for review and re-use ([https://github.com/opensafely/ics-research](https://github.com/opensafely/ics-research)). All iterations of the pre-specified study protocol are archived with version control ([https://github.com/opensafely/ics-research/tree/master/protocol](https://github.com/opensafely/ics-research/tree/master/protocol)). ### Patient and Public Involvement Patients were not formally involved in developing this specific study design that was developed rapidly in the context of a global health emergency. We have developed a publicly available website [https://opensafely.org/](https://opensafely.org/) through which we invite any patient or member of the public to contact us regarding this study or the broader OpenSAFELY project. The protocol and draft paper have been sent to the Asthma UK and British Lung Foundation Partnership for review from an expert-patient perspective. ## Results ### Patient Characteristics 148,488 people with COPD and 817,973 people with asthma and a relevant prescription within the four months before 01 Mar 2020 were included (supplemental figure 2-3). #### COPD Population Table 1 (see end) shows the characteristics of the COPD population. 43,278 (29%) received a LABA/LAMA prescription in the 4 months before index, and 105,210 (71%) prescriptions for ICS+LABA or ICS+LABA/LAMA. Demographic characteristics of treatment groups were similar; median age was 71 (IQR = 63 −77) in the LABA/LAMA group and 72 (IQR = 64 − 78) in the ICS group, and just over half were men (54.5% in the LABA/LAMA group and 53.7% in the ICS group). The presence of comorbidities was similar across the two treatment groups, except prior asthma diagnosis, which was more common among ICS recipients (12.9% in the LABA/LAMA group compared to 27.7% in the ICS group). The percentage of people with an exacerbation in the last year was lower in the LABA/LAMA group (19.7% vs 26.0%). View this table: [Table 1.](http://medrxiv.org/content/early/2020/06/20/2020.06.19.20135491/T1) Table 1. Demographic and Clinical Characteristics - COPD #### Asthma Population Table 2 (see end) shows the characteristics of the asthma population. Most (608,583, 74%) had a prescription for low/medium dose ICS in the four months prior to index date; fewer received high dose ICS (101,010, 12%) and SABA only (108,380, 13%). View this table: [Table 2.](http://medrxiv.org/content/early/2020/06/20/2020.06.19.20135491/T2) Table 2. Demographic and Clinical Characteristics - Asthma The asthma treatment groups differed in terms of demographic and clinical characteristics. The median age was 48 (IQR = 35 - 60), 53 (IQR = 40 - 66), and 55 (IQR = 44 - 67) for SABA only, low/medium dose ICS and high dose ICS respectively. The proportion of men was highest in the SABA group (43.0%), followed by the low/medium dose ICS (40.3%) and high dose ICS (38.2%). The prevalence of most comorbidities was lowest among people in the SABA only group and highest in the high dose ICS group. The proportion with an asthma exacerbation in the past year increased markedly across treatment groups (14.0%, 20.0% and 36.3% among SABA, low/medium dose and high dose ICS respectively). ### Univariable and Multivariable Results #### COPD 421 COVID-19 related deaths occurred in the treated COPD population; time to COVID-19 related death by treatment group can be seen in Figure 1a. In univariable models, ICS use was associated with an increased risk of COVID-19 related death (HR = 1.52, 95% CI = 1.20 - 1.91; Figure 2). This association decreased on adjustment for age and gender, and the remaining comorbidities (aHR = 1.38; 95% CI = 1.08 - 1.75; Figure 2). Post-hoc analyses revealed that adjustment for prior exacerbations as a binary variable and smoking status (current/former) had the largest impact on reducing the hazard ratio (supplemental table 6). There was no evidence of a (pre-specified) interaction with age (supplemental table 4). We detected significant deviations from the proportional hazards assumptions by testing for a zero slope in the scaled Schoenfeld residuals as well as by graphical inspection of plots of the Schoenfeld residuals against time (supplemental table 5, Figure 4-6). The KM curve indicated that the hazard ratio was likely above one throughout the follow-up, with the effect size growing over time. ![Figure 1a:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/20/2020.06.19.20135491/F1.medium.gif) [Figure 1a:](http://medrxiv.org/content/early/2020/06/20/2020.06.19.20135491/F1) Figure 1a: Kaplan-Meier plots in in the COPD population ![Figure 1b:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/20/2020.06.19.20135491/F2.medium.gif) [Figure 1b:](http://medrxiv.org/content/early/2020/06/20/2020.06.19.20135491/F2) Figure 1b: Kaplan-Meier plots in in the asthma population ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/20/2020.06.19.20135491/F3.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/06/20/2020.06.19.20135491/F3) Figure 2: Univariable and Multivariable Models, COPD population #### Asthma There were 515 COVID-19 related deaths in the treated asthma population; time to COVID-19 related death by treatment group can be seen in Figure 1b. In univariable models, receipt of both low/medium dose and high dose ICS was associated with an increased risk of COVID-19 related death (HR = 1.32, 95% CI = 0.98 - 1.77 and HR = 2.28, 1.62 - 3.20 respectively; Figure 3). These associations reduced markedly upon adjustment for age and gender, and the remaining pre-specified comorbidities (aHR = 1.10, 95% CI = 0.82 - 1.49 and aHR = 1.52, 95% CI = 1.08 - 2.14 for low/medium dose and high dose ICS respectively; Figure 3). Post-hoc analyses revealed that the greatest reduction in the strength of the association after the age and sex adjustment was adjustment for previous exacerbations (supplemental table 10). There was no evidence for a (pre-specified) interaction with age, and no deviations from the proportional hazards assumption (supplemental table 8 and 9, figure 7-12). ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/20/2020.06.19.20135491/F4.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/06/20/2020.06.19.20135491/F4) Figure 3: Univariable and Multivariable Models, Asthma population ### Sensitivity Analyses Figure 2 and 3 show results from sensitivity analyses. When considering receipt of ICS + LABA and ICS + LABA/LAMA separately in the COPD population, the risk of death was higher among those receiving ICS + LABA/LAMA (aHR = 1.41 95% CI = 1.10 - 1.81) but less markedly amongst those receiving ICS + LABA (aHR = 1.29, 95% CI = 0.95 - 1.74). Restricting analyses to people of white ethnicity to attempt to control for potential confounding by ethnicity led to a reduction in the hazard ratios in the COPD population (Figure 2), but not in the asthma population (Figure 3). Changing the population definition for the asthma population had a negligible impact on the results (Figure 3). ### Negative Control Analysis In the COPD population, the risk of non-COVID-19 related death was higher among individuals in the ICS combination group with an adjusted HR = 1.23 (95% CI = 1.07 - 1.40). In the asthma population, there was no evidence of an association between receipt of high or low/medium dose ICS and non-COVID-19 related death with adjusted HRs of 0.85 (95% CI = 0.71 - 1.00) and 0.93 (95% CI = 0.75 - 1.15), respectively. ### Quantitative Bias Analysis To fully explain the lower bound of the 95% CI in the COPD cohort (1.08), or for the high-dose ICS association in the asthma cohort (1.08), an unmeasured confounder would need to be associated (conditional on measured covariates) with either exposure or outcome by at least risk ratio (RR) 1.37 (e-value/high threshold) and with both exposure and outcome by at least RR 1.08 (low threshold; Figure 4a-b). An unmeasured confounder would need to have a stronger association with either exposure or outcome to move the observed HRs to a clinically meaningful protective effect of 0.8 (high threshold COPD: 2.84, asthma: 3.21, supplemental figure 14-15; low threshold COPD: 1.73, asthma: 1.90). ![Figure 4a:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/20/2020.06.19.20135491/F5.medium.gif) [Figure 4a:](http://medrxiv.org/content/early/2020/06/20/2020.06.19.20135491/F5) Figure 4a: E-value for the lower 95% CI and point estimate in the COPD ![Figure 4b:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/06/20/2020.06.19.20135491/F6.medium.gif) [Figure 4b:](http://medrxiv.org/content/early/2020/06/20/2020.06.19.20135491/F6) Figure 4b: E-value for the lower 95% CI and point estimate in the Asthma population ## Discussion ### Summary This is the first study to investigate the association between regular ICS use and COVID-19 related death. Compared with non-ICS based treatments, ICS in people with COPD was associated with a ∼40% increased risk of COVID-19 related death. In people with asthma, high dose ICS use was associated with an ∼50% increased risk of COVID-19 related death, with little evidence of any association for low/medium dose ICS. Our findings do not provide any strong support for a protective effect from ICS use in these populations, as has been previously hypothesised. Our analyses overall indicate that the observed harmful associations could readily be explained by confounding due to underlying health differences between people prescribed ICS and those using other medications for asthma and COPD, rather than a causally harmful effect of ICS. Specifically, we observed a stronger association with COVID-19 related death between ICS triple therapy than ICS dual therapy in the COPD cohort; the ICS content of these two regimens is similar, and a causal effect of ICS would be expected to be comparable in these two groups. If we had successfully controlled for disease severity differences between treatment groups, we would also expect to see no association between ICS and the negative control outcome of non-COVID-19 death. The harmful association we observed suggests we had not perfectly captured all markers of disease severity, resulting in an association that is unlikely to be causal. The null finding between ICS and non-COVID-19 related death in the asthma cohort is less surprising, since asthma is less markedly associated with overall mortality compared to COPD29. Finally, quantitative bias analysis confirms that a hypothetical unmeasured confounder of modest strength could fully explain the observed results. ### Findings in Context A literature review (box 2) found no other epidemiological studies or randomised controlled trials (RCTs) assessing the role of ICS in COVID-19. Two ongoing RCTs investigating the role of ICS in people hospitalised with laboratory confirmed COVID-19 and mild COVID-19 respectively ([NCT04331054](http://medrxiv.org/lookup/external-ref?link\_type=CLINTRIALGOV&access_num=NCT04331054&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F20%2F2020.06.19.20135491.atom), [NCT04330586](http://medrxiv.org/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT04330586&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F20%2F2020.06.19.20135491.atom)) are due to complete later this year. The hypothesis of a protective effect of ICS in COVID-19 was based partly on the prevalence of chronic lung disease among outpatient and inpatient COVID-19 in China. However, more recent studies do not support initial assertions that people with chronic lung diseases (including COPD) are significantly underrepresented among COVID-19 patients 6,7,30. In addition, a number of studies have found that people with COPD are at greater risk of severe COVID-19 and death from COVID-19 once infected 6,7. The evidence for asthma in COVID-19 is more varied, with studies reporting both null and moderately harmful associations6,7. It may be that features other than ICS use, such as shielding, influence the risk of acquiring SARS-CoV-2 among asthmatics. Studies investigating the causal effect of chronic respiratory disease, including COPD, and asthma on SARS-CoV-2 infection risk and COVID-19 disease, ideally taking into account the relatively large degree of heterogeneity that exists within each of these diagnostic categories, are urgently needed to help inform about levels of risk for these patients. ### Strengths and weaknesses The greatest strength of this study was the power we had to look at multiple drug treatments as our dataset included medical records from almost 24 million individuals. Our study is further strengthened by the use of two different study populations and active comparators, as well as sensitivity analyses to quantify the potential impact of unmeasured confounding on results. Another strength is our use of open methods: we pre-specified our analysis plan and have shared all analytical code. We also recognise possible limitations. The primary limitation is the risk of confounding by indication due to unmeasured or imperfectly defined potential confounding variables. Decisions regarding treatment choices involve factors that may not be well recorded in electronic health records including measures such as spirometry, and likely steroid responsiveness. As we did not have secondary care data our assessment of exacerbation history was incomplete, limiting our ability to adjust for this. Our sensitivity analyses confirm that unmeasured confounding is a plausible reason for the harmful associations we observed. The proportional hazards assumption was not met for the COPD models, with KM plots indicating that the hazard ratio for this exposure increased over time. This is perhaps not surprising, as the risk of acquiring COVID-19 was lower in the early stages of the pandemic. The HR for the COPD population should be interpreted as an average over the entire follow-up period. Finally, it is important to note that the outcome of COVID-19 related death will reflect the risk both of becoming infected as well as the risk of developing severe disease and dying. It is possible that ICS use has a different effect on the risk of infection and on disease severity. ### Policy Implications and Future Research We find no evidence that ICS has a strongly beneficial effect on COVID-19 related mortality, and therefore we cannot recommend that they are used to treat people with COVID-19 outside of the context of RCTs. Importantly, from the totality of the evidence provided here, including our sensitivity analyses, our results do not support the interpretation that regular ICS therapy for asthma or COPD increases risk of death from COVID-19, and do not provide evidence to support adjustments in ICS therapy among COVID-19 patients. Future observational studies of this clinical question are likely to face similar challenges around unmeasured confounding. The UK has an unusually large volume of detailed longitudinal patient data. We have demonstrated that it is feasible to rapidly address specific hypotheses about medicines in a transparent manner inside the secure environment of an EHR vendor in order to minimise the large volumes of potentially disclosive data that would otherwise have to move into separate systems. We will use the OpenSAFELY platform to further inform the global response about drug treatments during the COVID-19 emergency. ### Summary We found no evidence of a beneficial effect of regular ICS use on COVID-19 related mortality. Although we report a small harmful association, the pattern of results we observed suggests this could readily be explained by differences in underlying health between people receiving ICS and those receiving other respiratory medications. People currently taking ICS should continue taking them if recommended as part of routine care. **Box 1: Prespecified hypothetical confounders** * Age * Sex * Body Mass Index (BMI): Measurement of weight in the last decade * Indices of multiple deprivation (IMD): Quintiles from IMD 2019 * Diagnosed hypertension * Heart disease: categorised as heart failure and other heart disease * Diabetes: categorised as controlled (HbA1c < 58 mmol/mol), uncontrolled (HbA1c * ≥ 58 mmol/mol) or HbA1c not measured measured within the last 12 months * Cancer * Immunosuppressive conditions: organ transplant, sickle cell anaemia and splenectomy * Chronic kidney disease: based on creatinine measurements within the last 12 months or ever having a code for renal dialysis * Influenza vaccination status: recorded between 01Sep 2019 and 01 Mar 2020 * Pneumococcal vaccination status: record in the five years prior to 01 Mar 2020 * Statin use: recorded within four months prior to 1st of March 2020, * Exacerbation history: different methods used for asthma 31 and COPD population 32with COPD models additionally adjusted for a history of asthma. **Box 2: Research in Context** **Evidence before this study** At the start of the global coronavirus outbreak, Inhaled Corticosteroids (ICS) were hypothesised to offer some protection against either infection with SARS-CoV-2 or against severe outcomes from COVID-19, such as acute respiratory distress syndrome (ARDS) and respiratory failure9, despite these medications being known to increase the risk of pneumonia and other respiratory tract infections17,33. The hypothesis was based at least in part on epidemiological data showing a low prevalence of chronic respiratory disease among Chinese COVID-19 patients2, although there was also some support of a potential protective effect from in-vitro studies12,13. Most recently, ICS exposure has been found to correlate with a lower expression of ACE2 and TMPRSS2, the entry receptors used by SARS-CoV-2, in sputum cells34. A recent systematic review evaluating whether administration of ICS was associated with clinical outcomes in COVID-19, SARS or MERS identified no relevant studies10. **Added value of this study** Our study was specifically designed to assess the role of pre-morbid ICS use in COVID-19. We included two cohorts of participants: people with asthma, and people with COPD, both of whom have a possible indication for ICS. Neither analysis was strongly suggestive that regular ICS therapy for asthma or COPD has a clinically important causal effect on COVID-19 mortality in either direction. Our study has several key strengths: Firstly, it includes almost a million participants making it the largest contemporary study of ICS use to date. Secondly, we used active comparators and multiple sensitivity analyses to reduce and quantify the impact of possible unmeasured confounding. Finally, our analyses were pre-specified and we used open methods throughout the study with code and codelists available for examination and reuse. **Implications of all the available evidence** Evidence from our study and other research suggests there is neither a demonstrable benefit nor clear harm from ICS use against COVID-19 related mortality and therefore at present there is no evidence people should alter their ICS therapies during the pandemic. We also cannot recommend that ICS be used specifically to treat people with COVID-19 outside of the context of clinical trials. Future observational research is likely to be subject to similar issues around unmeasured confounding, and evidence from RCTs will provide more definitive answers regarding the role of ICS in COVID-19 later this year. ## Data Availability Data management was performed using Python 3.8 and SQL, with analysis carried out using Stata 16.1. All of the code used for data management and analyses is openly shared online for review and re-use ([https://github.com/opensafely/ics-research](https://github.com/opensafely/ics-research)). All iterations of the pre-specified study protocol are archived with version control ([https://github.com/opensafely/ics-research/tree/master/protocol](https://github.com/opensafely/ics-research/tree/master/protocol)) [https://github.com/opensafely/ics-research](https://github.com/opensafely/ics-research) ## Administrative ## 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). ## Funding This work was supported by the Medical Research Council MR/V015737/1. TPP provided technical expertise and infrastructure within their data centre *pro bono* in the context of a national emergency. 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. JPB 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. RM holds a Sir Henry Wellcome fellowship. EW holds grants from MRC. RG holds grants from NIHR and 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. 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 to submit the article for publication. ## Ethical approval This study was approved by the Health Research Authority (REC reference 20/LO/0651) and by the LSHTM Ethics Board (ref 21863). No further ethical or research governance approval was required by the University of Oxford but copies of the approval documents were reviewed and held on record. ## Guarantor BG ## Contributorship Contributions are as follows: Conceptualization LS BG ID; Data curation CB JP JC SH SB DE PI CM; Formal Analysis AS AJW BM CM JB; Funding acquisition BG LS; Information governance AM BG CB JP; Methodology ID AS LT CTR AW KB EW SJWE JQ LS JB CM AJW BM SB BG; Disease category conceptualisation and codelists CM AJW AS CTR PI SB DE CB JC JP SH HD HC KB SB AM BM LT ID HM RM HF JQ; Ethics approval HC EW LS BG; Project administration AS CM HC CB SB AM LS BG; Resources BG LS; Software SB DE PI AJW CM CB FH JC SH; Supervision ID LS BG; Visualisation CTR AJW; Writing (original draft) AS ID CM; Writing (review & editing) AS AJW BM CM CTR KB EW HJC HD SB CB AM DE PI HM LT RH KW AYSW HF RC SJWE JQ LS ID BG; All authors were involved in design and conceptual development and reviewed and approved the final manuscript. ## Acknowledgements We are very grateful for all the support received from the TPP Technical Operations team throughout this work; for generous assistance from the information governance and database teams at NHS England / NHSX. ## Footnotes * † Joint principal investigators * Received June 19, 2020. * Revision received June 19, 2020. * Accepted June 20, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.COVID-19 Map. Johns Hopkins Coronavirus Resource Center [https://coronavirus.jhu.edu/map.html](https://coronavirus.jhu.edu/map.html). 2. 2.Guan, W.-J. et al. Clinical Characteristics of Coronavirus Disease 2019 in China. N. Engl. J. Med. 382, 1708–1720 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2002032&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F20%2F2020.06.19.20135491.atom) 3. 3.Chen, N. et al. Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study. Lancet 395, 507–513 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30211-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F20%2F2020.06.19.20135491.atom) 4. 4.Garg, S. et al. Hospitalization Rates and Characteristics of Patients Hospitalized with Laboratory-Confirmed Coronavirus Disease 2019 - COVID-NET, 14 States, March 1-30, 2020. MMWR Morb. Mortal. Wkly. Rep. 69, 458–464 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.15585/mmwr.mm6915e3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F20%2F2020.06.19.20135491.atom) 5. 5.Richardson, S. et al. Presenting Characteristics, Comorbidities, and Outcomes Among 5700 Patients Hospitalized With COVID-19 in the New York City Area. JAMA (2020) doi:10.1001/jama.2020.6775. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2020.6775&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32320003&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F20%2F2020.06.19.20135491.atom) 6. 6.Docherty, A. B. et al. Features of 20 133 UK patients in hospital with covid-19 using the ISARIC WHO Clinical Characterisation Protocol: prospective observational cohort study. BMJ 369, m1985 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzNjkvbWF5MjJfMS9tMTk4NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA2LzIwLzIwMjAuMDYuMTkuMjAxMzU0OTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 7. 7.The OpenSAFELY Collaborative et al. OpenSAFELY: factors associated with COVID- 19-related hospital death in the linked electronic health records of 17 million adult NHS patients. Epidemiology (2020) doi:10.1101/2020.05.06.20092999. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNS4wNi4yMDA5Mjk5OXYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDYvMjAvMjAyMC4wNi4xOS4yMDEzNTQ5MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 8. 8.Reilev, M. et al. Characteristics and predictors of hospitalization and death in the first 9,519 cases with a positive RT-PCR test for SARS-CoV-2 in Denmark: A nationwide cohort. Epidemiology (2020) doi:10.1101/2020.05.24.20111823. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNS4yNC4yMDExMTgyM3YxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDYvMjAvMjAyMC4wNi4xOS4yMDEzNTQ5MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 9. 9.Halpin, D. M. G., Faner, R., Sibila, O., Badia, J. R. & Agusti, A. Do chronic respiratory diseases or their treatment affect the risk of SARS-CoV-2 infection? Lancet Respir Med 8, 436–438 (2020). 10. 10.Halpin, D. M. G., Singh, D. & Hadfield, R. M. Inhaled corticosteroids and COVID-19: a systematic review and clinical perspective. Eur. Respir. J. 55, (2020). 11. 11.Inhaled corticosteroids and COVID-19 outcomes in airways disease. (2020). 12. 12.Matsuyama, S. et al. The inhaled corticosteroid ciclesonide blocks coronavirus RNA replication by targeting viral NSP15. bioRxiv 2020.03.11.987016 (2020) doi:10.1101/2020.03.11.987016. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4wMy4xMS45ODcwMTZ2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA2LzIwLzIwMjAuMDYuMTkuMjAxMzU0OTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 13. 13.Yamaya, M. et al. Inhibitory effects of glycopyrronium, formoterol, and budesonide on coronavirus HCoV-229E replication and cytokine production by primary cultures of human nasal and tracheal epithelial cells. Respir. Investig. 58, 155–168 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.resinv.2019.12.005&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F20%2F2020.06.19.20135491.atom) 14. 14.Low-cost dexamethasone reduces death by up to one third in hospitalised patients with severe respiratory complications of COVID-19. [https://www.recoverytrial.net/news/low-cost-dexamethasone-reduces-death-by-up-to-one-third-in-hospitalised-patients-with-severe-respiratory-complications-of-covid-19](https://www.recoverytrial.net/news/low-cost-dexamethasone-reduces-death-by-up-to-one-third-in-hospitalised-patients-with-severe-respiratory-complications-of-covid-19). 15. 15.Ernst, P., Gonzalez, A. V., Brassard, P. & Suissa, S. Inhaled corticosteroid use in chronic obstructive pulmonary disease and the risk of hospitalization for pneumonia. Am. J. Respir. Crit. Care Med. 176, 162–166 (2007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/rccm.200611-1630OC&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17400730&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F20%2F2020.06.19.20135491.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248138300011&link_type=ISI) 16. 16.Suissa, S., Patenaude, V., Lapi, F. & Ernst, P. Inhaled corticosteroids in COPD and the risk of serious pneumonia. Thorax 68, 1029–1036 (2013). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToidGhvcmF4am5sIjtzOjU6InJlc2lkIjtzOjEwOiI2OC8xMS8xMDI5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDYvMjAvMjAyMC4wNi4xOS4yMDEzNTQ5MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 17. 17.Crim, C. et al. Pneumonia risk in COPD patients receiving inhaled corticosteroids alone or in combination: TORCH study results. Eur. Respir. J. 34, 641–647 (2009). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiZXJqIjtzOjU6InJlc2lkIjtzOjg6IjM0LzMvNjQxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDYvMjAvMjAyMC4wNi4xOS4yMDEzNTQ5MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 18. 18.Corticosteroids - inhaled - NICE CKS. [https://cks.nice.org.uk/corticosteroids-inhaled](https://cks.nice.org.uk/corticosteroids-inhaled). 19. 19.England, N. H. S. NHS England » Coronavirus (COVID-19) Research Platform. [https://www.england.nhs.uk/contact-us/privacy-notice/how-we-use-your-information/covid-19-response/coronavirus-covid-19-research-platform/](https://www.england.nhs.uk/contact-us/privacy-notice/how-we-use-your-information/covid-19-response/coronavirus-covid-19-research-platform/). 20. 20.Quint, J. K. et al. Validation of chronic obstructive pulmonary disease recording in the Clinical Practice Research Datalink (CPRD-GOLD). BMJ Open 4, e005540 (2014). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYm1qb3BlbiI7czo1OiJyZXNpZCI7czoxMToiNC83L2UwMDU1NDAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNi8yMC8yMDIwLjA2LjE5LjIwMTM1NDkxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 21. 21.Nissen, F. et al. Concomitant diagnosis of asthma and COPD: a quantitative study in UK primary care. Br. J. Gen. Pract. 68, e775–e782 (2018). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYmpncCI7czo1OiJyZXNpZCI7czoxMToiNjgvNjc2L2U3NzUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNi8yMC8yMDIwLjA2LjE5LjIwMTM1NDkxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 22. 22.Overview | Chronic obstructive pulmonary disease in over 16s: diagnosis and management | Guidance | NICE. 23. 23.OpenSAFELY Codelists: High dose ICS Inhalers. [https://codelists.opensafely.org/codelist/opensafely/high-dose-ics-inhalers/](https://codelists.opensafely.org/codelist/opensafely/high-dose-ics-inhalers/). 24. 24.Janson, C. et al. SABINA: An Overview of Short-Acting β2-Agonist Use in Asthma in European Countries. Adv. Ther. 37, 1124–1135 (2020). 25. 25.Lund, J. L., Richardson, D. B. & Stürmer, T. The active comparator, new user study design in pharmacoepidemiology: historical foundations and contemporary application. Curr Epidemiol Rep 2, 221–228 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s40471-015-0053-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26954351&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F20%2F2020.06.19.20135491.atom) 26. 26.WHO | Emergency use ICD codes for COVID-19 disease outbreak. (2020). 27. 27.Austin, P. C., Lee, D. S. & Fine, J. P. Introduction to the Analysis of Survival Data in the Presence of Competing Risks. Circulation 133, 601–609 (2016). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjk6IjEzMy82LzYwMSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA2LzIwLzIwMjAuMDYuMTkuMjAxMzU0OTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 28. 28.VanderWeele, T. J. & Ding, P. Sensitivity Analysis in Observational Research: Introducing the E-Value. Ann. Intern. Med. 167, 268–274 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/M16-2607&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28693043&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F20%2F2020.06.19.20135491.atom) 29. 29.GBD 2015 Chronic Respiratory Disease Collaborators. Global, regional, and national deaths, prevalence, disability-adjusted life years, and years lived with disability for chronic obstructive pulmonary disease and asthma, 1990-2015: a systematic analysis for the Global Burden of Disease Study 2015. Lancet Respir Med 5, 691–706 (2017). 30. 30.MMWR - Preliminary Estimates of the Prevalence of Selected Underlying Health Conditions Among Patients with Coronavirus Disease 2019 — United States, February 12–March 28, 2020. 31. 31.Bloom, C. I., Palmer, T., Feary, J., Quint, J. K. & Cullinan, P. Exacerbation Patterns in Adults with Asthma in England. A Population-based Study. Am. J. Respir. Crit. Care Med. 199, 446–453 (2019). 32. 32.Rothnie, K. J. et al. Validation of the Recording of Acute Exacerbations of COPD in UK Primary Care Electronic Healthcare Records. PLoS One 11, e0151357 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0151357&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26959820&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F20%2F2020.06.19.20135491.atom) 33. 33.Chalmers, J. D. & Keir, H. R. 10 years since TORCH: shining a new light on the risks of inhaled corticosteroids in COPD. The European respiratory journal: official journal of the European Society for Clinical Respiratory Physiology vol. 50 (2017). 34. 34.Peters, M. C. et al. COVID-19 Related Genes in Sputum Cells in Asthma: Relationship to Demographic Features and Corticosteroids. Am. J. Respir. Crit. Care Med. (2020) doi:10.1164/rccm.202003-0821OC. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1164/rccm.202003-0821OC&link_type=DOI)