Associations of severe COVID-19 with polypharmacy in the REACT-SCOT case-control study ====================================================================================== * Paul M McKeigue * Sharon Kennedy * Amanda Weir * Jen Bishop * Stuart J McGurnaghan * David McAllister * Chris Robertson * Rachael Wood * Nazir Lone * Janet Murray * Thomas M Caparrotta * Alison Smith-Palmer * David Goldberg * Jim McMenamin * Colin Ramsay * Bruce Guthrie * Sharon Hutchinson * Helen M Colhoun ## Abstract **Objectives** To investigate the relation of severe COVID-19 to prior drug prescribing. **Design** Matched case-control study (REACT-SCOT) based on record linkage to hospital discharges since June 2015 and dispensed prescriptions issued in primary care during the last 240 days. **Setting** Scottish population. **Main outcome measure** Severe COVID-19, defined by entry to critical care or fatal outcome. **Participants** All 4272 cases of severe COVID-19 in Scotland since the start of the epidemic, with 36948 controls matched for age, sex and primary care practice. **Results** Severe COVID-19 was strongly associated with the number of non-cardiovascular drug classes dispensed. This association was strongest in those not resident in care homes, in whom the rate ratio (95% CI) associated with dispensing of 12 or more drug classes versus none was 10.8 (8.7, 13.2), and was not accounted for by treatment of conditions designated as conferring increased risk. Of 17 drug classes postulated at the start of the epidemic to be “medications compromising COVID”, all were associated with increased risk of severe COVID-19. The largest effect was for antipsychotic agents: rate ratio 4.14 (3.39, 5.07). Other drug classes with large effects included proton pump inhibitors (rate rato 2.19 (1.70, 2.80) for >= 2 defined daily doses/day), opioids (3.62 (2.65, 4.94) for >= 50 mg morphine equivalent/day) and gabapentinoids. These associations persisted after adjusting for covariates, and were stronger with recent than with non-recent exposure. **Conclusions** Severe COVID-19 is associated with polypharmacy and with drugs that cause sedation, respiratory depression or dyskinesia, have anticholinergic effects or affect the gastrointestinal system. These associations are not easily explained by co-morbidity. Although the evidence for causality is not conclusive, these results support existing guidance on reducing overprescribing of these drug classes and limiting inappropriate polypharmacy as a potential means of reducing COVID-19 risk. **Registration** ENCEPP number EUPAS35558 **What is already known on this topic** Two previous studies have examined the relationship of severe COVID-19 to drugs for the cardiovascular system. This is the first systematic study of the relationship of severe COVID-19 to prior drug prescribing. **What this study adds** Severe COVID-19 is associated with polypharmacy and with drugs that cause sedation, respiratory depression or dyskinesia, have anticholinergic effects or affect the gastrointestinal system. These associations are not easily explained by co-morbidity. These results support earlier warnings that these drug classes that these drugs might increase susceptibility to COVID-19, and reinforce existing guidance on reducing overprescribing of these drug classes. ## Background In the initial analysis of REACT-SCOT, a matched case control study of risk factors for severe COVID-19 in Scotland, we reported a strong association of severe COVID-19 with having had at least one prescription dispensed in the past year1. The univariate rate ratios associated with at least one prescription varied from 3.8 in those aged under 60 years to 2.3 in those aged 75 years and over. This association persisted after adjusting for care home residence, hospital admission in the last five years, and diagnoses of conditions designated by public health agencies as conferring vulnerability to COVID-19 (hereafter listed conditions). The objective of this study was to investigate this association. ## Methods The design of the REACT-SCOT case-control study has been described in detail elsewhere1. As this study was set up in response to a national emergency no arrangements were made for patient and public involvement. All individuals testing positive for nucleic acid for SARS-CoV-2 in Scotland were ascertained through the Electronic Communication of Surveillance in Scotland (ECOSS) database. Admissions to critical care were obtained from the Scottish Intensive Care Society and Audit Group (SICSAG) database that captures admission to all critical care (intensive care or high dependency) units. Death registrations were obtained from linkage to the National Register of Scotland. Severe or fatal COVID-19 was defined by either (1) a positive nucleic acid test followed by entry to critical care or death within 28 days; or (2) a death certificate with COVID-19 as underlying cause. For each case, the Community Health Index database was used to select up to ten controls matched for sex, one-year age band and registered with the same primary care practice, who were alive on the same day as the first date that the case tested positive. For fatal cases who had not tested positive, the incident date was assigned as 14 days before death. For this analysis based on ascertainment of positive test results up to 6 June 2020, entry to critical care up to 14 June 2020 and deaths up to 12 June 2020 there were 4272 cases and 36948 controls. ### Morbidity and drug prescribing For all cases and controls, ICD-10 diagnostic codes were extracted from the last five years of hospital discharge records in the Scottish Morbidity Record (SMR01), excluding records of discharges less than 25 days before testing positive for SARS-CoV-2, and from the national cancer registry. Diagnoses of diabetes were extracted from linkage to the national diabetes register. British National Formulary (BNF) drug codes for dispensed prescriptions issued in primary care were extracted from the Scottish Prescribing Information System2. A cutoff date of 15 days before the incident date (date of testing positive for SARS-CoV-2, or 14 days before death for fatal cases without a positive test) was set, and prescriptions dispensed in a 240-day interval before this cutoff date were included. For this analysis prescription codes from BNF chapters 14 and above, comprising dressings, appliances, vaccines, anaesthesia and other preparations were grouped as “Other”. We began by testing for association of severe COVID-19 with the number of drug classes (BNF subparagraph codes) for which at least one prescription had been dispensed during the period of observation. In accordance with earlier suggestions that prescribing of multiple cardiovascular drugs, which is supported by evidence-based guidelines, should be considered separately from putatively inappropriate polypharmacy3, we partitioned the number of drug classes dispensed into cardiovascular and other drugs. To define a prespecified hypothesis about drug classes postulated to increase the risk of severe COVID-19, we used a review by Laporte and Healy published on 2 April 2020 that listed 17 widely prescribed drug classes associated with increased pneumonia risk and therefore of concern as “medications compromising COVID”4. Excluding immunosuppressive drugs that are criteria for shielding, the Laporte-Healy list comprised proton pump inhibitors, gastrointestinal antispasmodics, H1 antihistamines, hypnotics and sedatives, antipsychotic drugs, antidepressants, drugs used in nausea and vertigo, opioid analgesics, gabapentinoids, anti-epileptic drugs, antimuscarinic drugs used in parkinsonism, urinary antispasmodics, and non-steroidal anti-inflammatory drugs. For proton pump inhibitors and opioids, it was possible to use equivalent doses to calculate a total dose for all drugs in the class. Defined daily doses (DDDs) of each proton pump inhibitor were obtained from the [DDD/ATC Index](http://WHO%20Collaborating%20Centre%20for%20Drug%20Statistics) of the WHO Collaborating Centre for Drug Statistics Methodology. Morphine milligram equivalents (MMEs) for opioids were obtained from the website of the [Faculty of Pain Medicine](http://WHO%20Collaborating%20Centre%20for%20Drug%20Statistics) of the Royal College of Anaesthestists. Average daily doses of proton pump inhibitors (as DDDs) and opioids (as MME) were calculated for each individual as the sum of (conversion factor *×* strength *×* quantity dispensed) divided by observation period. The dose of opioid was calculated as the sum over opioid-containing items in subparagraph 0407010 (nonopioid and compound preparations) and 0407020 (opioid analgesics). As described previously, we derived indicator variables for a list of conditions that have been designated as risk conditions for COVID-19 by public health agencies5: diabetes, heart disease, asthma or chronic obstructive airway disease, chronic kidney disease, disabling neurological disease, liver disease and immunodeficiency or immunosuppression. ICD-10 diagnostic and BNF drug codes used to derive these conditions are available with the [ENCEPP registration](http://www.encepp.eu/encepp/openAttachment/documents.otherDocument-0/35565). Socioeconomic status was encoded as the quintile of the postcode-based Scottish Index of Multiple Deprivation (SIMD). ## Statistical analysis Rate ratios for severe COVID-19 (hereafter COVID-19) were estimated from conditional logistic regression models, implemented as Cox regression in the R function survival::clogit. To distinguish between causality and confounding, we used several approaches: * Testing for consistency of association with drugs that have a similar mode of action across different indication groups * Testing for a dose-response relationship and stratifying by age group * Adjusting for prespecified covariates. * Comparison between recent (last 120 days of observation) and less recent (1 June 2019 to 121 days before exposure) time windows of dispensing. ## Results ### Relation of severe COVID-19 to polypharmacy Figure 1 shows that the rate of of severe COVID-19 increased steeply with the number of non-cardiovascular drug classes and decreased with the number of noncardiovascular drug classes dispensed. Table 1 shows that this association was restricted to those not resident in care homes. Supplementary Table S1 shows that among those not resident in care homes, the association of severed COVID-19 with number of non-cardiovascular drug classes dispensed was present both in those without any diagnosed listed condition and in those with at least one diagnosed listed condition. ![Fig 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/27/2020.07.23.20160747/F1.medium.gif) [Fig 1.](http://medrxiv.org/content/early/2020/07/27/2020.07.23.20160747/F1) Fig 1. Rate ratios (with standard errors) in a conditional logistic regression of severe COVID-19 on number of cardiovascular (BNF chapter 4) and non-cardiovascular drug classes dispensed View this table: [Table 1.](http://medrxiv.org/content/early/2020/07/27/2020.07.23.20160747/T1) Table 1. Association of severe COVID-19 with number of non-cardiovascular drug classes dispensed, by care home residence ### Associations with specific drug classes In order to identify specific drug classes associated with severe COVID-19, we examined thse associations among those not resident in a care home and without any listed condition, as most people diagnosed with listed conditions such as asthma and diabetes are on medication for these conditions. Table 2 shows univariate associations of dispensing of at least one drug in each BNF subparagraph, filtered to show only drug classes with at least 50 exposed individuals and *p <* 0.001. The drug classes associated with severe COVID-19 include proton pump inhibitors, laxatives, multiple classes of drugs acting on the central nervous system, nutritional supplements, and non-steroidal anti-inflammatory drugs. View this table: [Table 2.](http://medrxiv.org/content/early/2020/07/27/2020.07.23.20160747/T2) Table 2. Univariate rate ratios for severe COVID-19 associated with dispensing of each BNF subparagraph code, restricted to those without a listed condition and not resident in a care home, filtered to show only drug classes with at least 50 exposed individuals and *p <* 0.001: Table 3 shows associations of severe COVID-19 with dispensing of each of the drug classes on the Laporte-Healy list in those not resident in a care home. All drug classes on this list were associated with increased risk in univariate analyses (though the association with NSAIDS was stronger in Table 2 where listed conditions were excluded. This included all the drug classes listed by Laporte and Healy as having anticholinergic effects: H1 antihistamines, antidepressants, urinary antispasmodics, gastrointestinal antispasmodics, drugs for vertigo, antimuscarinic drugs used in the treatment of parkinsonism, and antiepileptic drugs. View this table: [Table 3.](http://medrxiv.org/content/early/2020/07/27/2020.07.23.20160747/T3) Table 3. Associations of severe COVID-19 with drug classes listed by Laporte and Healy (2020), restricted to those not resident in a care home. Adjusted rate ratios are based on a joint model with all drug classes in the table and SIMD quintile as covariates In a multivariable regression, the strongest independent associations were with proton pump inhibitors, antihistamines, antipsychotic drugs, and opioid analgesics. In both univariate and multivariable analyses, the highest rate ratio was that associated with antipsychotic drugs. As others have noted4, the chemical structures and modes of action of drugs used in the treatment of nausea and vertigo overlap with those of antipsychotic drugs. Supplementary Table S2 shows the univariate associations of severe COVID-19 with specific drugs classified in these two subparagraphs of the BNF. Across both groups of indications, phenothiazines and other drugs that are dopamine antagonists were strongly associated with increased rates of severe disease. Rate ratios were elevated both for phenothiazines and for second-generation antipsychotics: aripiprazole, olanzapine, risperidone and amrisulpride. ### Time window analyses Supplementary Table S3 shows for each drug class on the Laporte-Healy list the rate ratio associated with dispensing in the most recent 120-day time window, with dispensing only in the previous time window as reference category. Because most users of these drugs had dispensed prescriptions in both time windows, these analyses are based on relatively small numbers: the drug classes shown are restricted to those with at least 500 cases and controls exposed only in the earlier time window. For most of these drug classes shown, the rate ratio associated with recent exposure only is above 1, but only for opioid analgesics, and to some extent proton pump inhibitors and anxiolytics, does this association reach conventional levels of statistical significance. ### Dose-response analyses Supplementary Table S4 shows the relationship of severe COVID-19 to average daily doses of opioids (as MME) and proton pump inhibitors (as DDDs) over the 240-day observation period. For both these drug classes there were dose-response relationships, and for proton pump inhibitors this relationship was strongest in those aged less than 75 years. With unexposed as baseline, the univariate rate ratio associated with opioid use in this age group was 3.62 (2.65, 4.94) in those with average daily dose of more than 50 mg morphine equivalent (MME), reduced to 2.90 (2.04, 4.11) on adjusting for care home residence, SIMD quintile and any history of neoplasm. For proton pump inhibitors, the univariate rate ratio associated with average dose of 2 or more DDDs / day was 2.19 (1.70, 2.80), reduced to 1.98 (1.50, 2.60) by adjusting for care home residence, SIMD quintile, any diagnosis of ICD-10 codes K20-K31 (diseases of esophagus, stomach and duodenum), non-steroidal anti-inflammatory drugs, anti-platelet agents and anticoagulants. ### Associations with other drug classes Supplementary Tables S5 and S6 show associations with drug classes in BNF chapters 2 (cardiovascular) and 10 (musculoskeletal). These chapters were selected as of interest because because specific hypotheses about possible effects of drugs in these chapters – ACE inhibitors6, anticoagulants7, and hydroxychloroquine8– have been proposed or discussed. Table S5 shows associations with drugs for the cardiovascular system. Prescriptions of loop diuretics and anticoagulants were associated with elevated rate ratios for severe COVID-19 in univariate and multivariate analyses. Over all age groups combined, the univariate rate ratio associated with oral anticoagulants was reduced from 1.88 (1.66, 2.12) to 1.28 (1.12, 1.46) by adjustment for ischaemic heart disese, other heart disease and ever-use of a proton pump inhibitor. Of drug classes commonly used to treat hypertension, ACE inhibitors and angiotensin II receptor antagonists were associated with reduced risk of COVID-19. Table S6 shows associations with drugs for the musculoskeletal system disaggregated by generic name, as the BNF groups all disease-modifying antirheumatic drugs under a single subparagraph code. All antirheumatic drugs were associated with elevated rate ratios for COVID-19. The univariate rate ratio associated with hydroxychloroquine sulfate was 2.10 (1.34, 3.28); adjustment for ever-use of a proton pump inhibitor reduced this to 1.80 (1.15, 2.82). ## Discussion We have shown that severe COVID-19 is associated with polypharmacy, defined by the number of drugs classes dispensed during the period of observation, in individuals without conditions designated as conferring high risk. The rate ratios of 5 to 7 associated with dispensing of more than 10 drug classes are larger than the ratio of about 2 for all-cause mortality associated with this level of polypharmacy in a systematic review9. Attempting to investigate associations with specific drugs with a hypothesis-free approach is difficult because many of the drug classes that are strongly associated with severe COVID-19, such as proton pump inhibitors, opioids and gabapentinoids, are indicators of overprescribing, recognized as such in the Scottish National Therapeutic Index of prescribing quality10. To narrow the hypothesis space we tested a pre-specified list of drugs postulated at the start of the epidemic to increase risk of severe COVID-19, based on previously described associations with pneumonia or activity on relevant pathways, especially anticholinergic agents4. We have shown that all the drugs originally listed are associated with increased risk. Distinguishing between causality and confounding as possible explanations for the association of COVID-19 with these drug classes is difficult because drug exposures are not measured accurately, and the associations are likely to be confounded by other drugs, by co-morbidities and more generally by frailty and socioeconomic deprivation. The strongest association was with antipsychotic drugs. As the recognized indications for these drugs are not known to be associated with susceptibility to severe COVID-19, there is no obvious explanation for this association other than causality. Further evidence favouring causality is that drugs prescribed for nausea that have similar chemical structures and mode of action to antipsychotic drugs show similar associations with severe COVID-19. For both opioids and proton pump inhibitors, four criteria provide moderate but not conclusive evidence favouring causal explanations over confounding. There are strong dose-response relationships of COVID-19 to dispensed average daily dose. Adjustment for covariates pre-specified as likely to confound these associations reduces the effect size only slightly. The effect sizes are larger in younger individuals and when the analysis is restricted to those without any of the designated risk conditions for COVID-19. Among ever-exposed individuals, the rate ratios associated with dispensing only in the most recent time window were higher than the rate ratios associated with dispensing only in an earlier time window. The dose-response relationship of opioid use to COVID-19 is similar in magnitude to that reported for community-acquired pneumonia in a study of people receiving medical care through the Veterans Administration from 2000-201211. Although gabapentinoids are classified in the BNF under subparagraph 0408010 (“Control of epilepsy”), in Scotland they are widely used in combination with or as substitutes for opioid analgesics. A limitation of this study is that we do not have morbidity data from primary care, which would include risk factors such as smoking and coding of presenting complaints. Another limitation is that it is not possible to capture hospital prescribing data which includes biologic agents that have immunosuppressive effects. Strengths of this study are that diagnoses are based on hospital discharge records coded to ICD-10 (rather than the SNOMED-CT codes used in primary care databases), and that drug exposure is based on dispensed rather than issued prescriptions. The mechanisms postulated by Laporte and Healy for drugs to increase risk of severe COVID-19 include sedation, respiratory depression, respiratory dyskinesia and anticholinergic effects4. We note that as SARS-CoV-2 is at least partly an enteric infection12 and the ACE2 receptor is expressed in the intestine, it is plausible that proton pump inhibitors and other drugs acting on the gastrointestinal tract could increase susceptibility to severe infection. It may be relevant to investigate these associations in other countries where COVID-19 epidemics have been especially severe and overprescribing of drug classes such as proton pump inhibitors13–15 or opioids16 has been reported previously. We emphasize that because of the relationship of COVID-19 to polypharmacy, associations with specific drug classes cannot be studied without taking into account how those drug classes are related to the profile of drug prescribing. Of those drug classes that are not on the Laporte-Healy list, there are strong univariate associations of severe COVID-10 with dispensing of antibiotics, laxatives and nutritional supplements. The associations with nutritional supplements is likely to be confounded by overprescribing. For ACE inhibitors and angiotensin-II receptor blockers, our results are consistent with other studies that have found no increased risk associated with these drugs6,17 and indeed suggest some protective effect may be possible. To explore this more fully will require access to other datasets where measurements of blood pressure and other covariates are available. The relation of anticoagulant use to COVID-19 is of interest because coagulopathy is a feature of severe disease7. Although in this study anticoagulants were associated with increased risk in univariate analysis, this association was reduced by adjusting for covariates including diagnosed heart disease and co-prescribing of proton pump inhibitors. Associations with disease-modifying antirheumatic drugs such as hydroxychloroquine are likely to be confounded by hospital-based prescribing of biologic anti-rheumatic drugs, which is not captured by record linkage in Scotland. ### Conclusion Severe COVID-19 is strongly associated with polypharmacy. This association is not easily explained by co-morbidity, and it is strongest in those without hospital diagnoses of conditions that confer high risk of disease. As a prediction of which drug classes would be associated with increased susceptibility to severe COVID-19, the Laporte-Healy list prepared at the start of the epidemic appears to be remarkably accurate. Many of the drug classes on this list are recognized indicators of overprescribing. The consistency of associations with drugs that have similar modes of action across different groups of indication, the dose-response and time window effects support causal explanations for at least some these associations. We recommend that public health agencies should reinforce existing guidelines on avoidance of overprescribing of these drug classes and more generally on inappropriate polypharmacy. ## Declarations ### Ethics approval and information governance This study was conducted under approvals from the Public Benefit and Privacy Panel for Health and Social Care that allow Public Health Scotland staff to link datasets. Datasets were de-identified before analysis. ### Data sharing Requests for acccess to the data may be submitted to the Public Benefit and Privacy Panel. ### Role of the funding source This study was undertaken without any external funding. ### Transparency declaration PM, as the manuscript’s guarantor, affirms: that the manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as originally planned and registered have been explained. This manuscript has been generated directly from the source data by a reproducible research pipeline. All source code used for derivation of variables, statistical analysis and generation of this manuscript is freely available on [GitHub](http://github.com/pmckeigue/covid-scotland_public). ### Registration The study protocol was registered with the European Network of Centres for Pharmacoepidemiology and Pharmacovigilance (ENCEPP number [EUPAS35558](http://www.encepp.eu/encepp/viewResource.htm?id=35559)). ## Data Availability Applications for access to the data may be made to the Public Benefit and Privacy Panel for Health and Social Care, Scotland ## Conflicts of interest HC receives research support and honoraria and is a member of advisory panels or speaker bureaus for Sanofi Aventis, Regeneron, Novartis, Novo-Nordisk and Eli Lilly. HC receives or has recently received non-binding research support from AstraZeneca and Novo-Nordisk. SH received honoraria from Gilead. ## Public Health Scotland COVID-19 Health Protection Study Group Alice Whettlock1, Allan McLeod1, Andrew Gasiorowski1, Andrew Merrick1, Andy McAuley1, April Went1, Calum Purdie1, Colin Fischbacher1, Colin Ramsay1, David Bailey1, David Henderson1, Diogo Marques1, Eisin McDonald1, Genna Drennan1, Graeme Gowans1, Graeme Reid1, Heather Murdoch1, Jade Carruthers1, Janet Fleming1, Jade Carruthers1, Joseph Jasperse1, Josie Murray1, Karen Heatlie1, Lindsay Mathie1, Lorraine Donaldson1, Martin Paton1, Martin Reid1, Melissa Llano1, Michelle Murphy-Hall1, Paul Smith1, Ros Hall1, Ross Cameron1, Susan Brownlie1, Adam Gaffney2, Aynsley Milne2, Christopher Sullivan2, Edward McArdle2, Elaine Glass2, Johanna Young2, William Malcolm2, Jodie McCoubrey 2 1 Health Protection Scotland (Public Health Scotland), Meridian Court, 5 Cadogan Street, Glasgow G2 6QE. 2 NHS National Services Scotland, Meridian Court, 5 Cadogan Street, Glasgow G2 6QE. ### Supplementary material The R and Rmarkdown scripts used to generate this article, which include the code used to derive variables, will be made available when the paper is published. ## Supplementary Tables View this table: [Table S1.](http://medrxiv.org/content/early/2020/07/27/2020.07.23.20160747/T4) Table S1. Association of severe COVID-19 with number of non-cardiovascular drug classes dispensed, in those not resident in a care home by diagnosis of at least one listed condition View this table: [Table S2.](http://medrxiv.org/content/early/2020/07/27/2020.07.23.20160747/T5) Table S2. Univariate associations of severe COVID-19 with drugs in four subparagraphs of BNF chapter 4, in those not resident in a care home and not treated for cancer in last year, filtered to retain rows with at least 20 exposed individuals View this table: [Table S3.](http://medrxiv.org/content/early/2020/07/27/2020.07.23.20160747/T6) Table S3. Comparison of associations with exposure only in the most recent 120-day time window with exposure only in the previous 120-day time window (baseline category for rate ratios) for drug classes on the Laporte-Healy list with at least 500 cases and controls exposed only in the earlier time window View this table: [Table S4.](http://medrxiv.org/content/early/2020/07/27/2020.07.23.20160747/T7) Table S4. Associations of severe COVID-19 with dispensed average daily dose of opioids and proton pump inhibitors, with and without adjustment for prespecified covariates View this table: [Table S5.](http://medrxiv.org/content/early/2020/07/27/2020.07.23.20160747/T8) Table S5. Associations of severe COVID-19 with prescribed drugs in BNF chapter 2 by subparagraph code, restricted to those not resident in a care home and filtered to retain rows with at least 50 exposed individuals View this table: [Table S6.](http://medrxiv.org/content/early/2020/07/27/2020.07.23.20160747/T9) Table S6. Associations of severe COVID-19 in with prescribed drugs in BNF chapter 10 by approved name, restricted to those not resident in a care home and filtered to retain rows with at least 50 exposed individuals ## Acknowledgements We thank all staff in critical care units who submitted data to the SICSAG database, the Scottish Morbidity Record Data Team, the staff of the National Register of Scotland, the Public Health Scotland Terminology Services, the HPS COVID-19 Laboratory & Testing cell and the NHS Scotland Diagnostic Virology Laboratories, and Nicola Rowan (HPS) for coordinating this collaboration. * Received July 23, 2020. * Revision received July 23, 2020. * Accepted July 27, 2020. * © 2020, 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. 1.McKeigue PM, Weir A, Bishop J, et al. Rapid Epidemiological Analysis of Comorbidities and Treatments as risk factors for COVID-19 in Scotland (REACT-SCOT): A population-based case-control study. medRxiv 2020;2020.05.28.20115394. 2. 2.Alvarez-Madrazo S, McTaggart S, Nangle C, Nicholson E, Bennie M. Data Resource Profile: The Scottish National Prescribing Information System (PIS). International Journal of Epidemiology 2016;45(3):714–715f. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyw060&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27165758&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F27%2F2020.07.23.20160747.atom) 3. 3.Appleton SC, Abel GA, Payne RA. Cardiovascular polypharmacy is not associated with unplanned hospitalisation: Evidence from a retrospective cohort study. BMC family practice 2014;15:58. 4. 4.Laporte J-R, Healy D. In the midst of the SARS-CoV-2 pandemia, caution is needed with drugs that increase the risk of pneumonia. Data-based Medicine Americas Ltd; 2020. 5. 5.NHS. Who’s at higher risk from coronavirus - Coronavirus (COVID-19). nhs.uk. 6. 6.Mancia G, Rea F, Ludergnani M, Apolone G, Corrao G. ReninAngiotensinAldosterone System Blockers and the Risk of Covid-19. New England Journal of Medicine 2020;():ull. 7. 7.Connors JM, Levy JH. COVID-19 and its implications for thrombosis and anticoagulation. Blood 2020; 8. 8.Alexander PE, Debono VB, Mammen MJ, et al. COVID-19 research has overall low methodological quality thus far: Case in point for chloroquine/hydroxychloroquine. Journal of Clinical Epidemiology 2020; 9. 9.Leelakanok N, Holcombe AL, Lund BC, Gu X, Schweizer ML. Association between polypharmacy and death: A systematic review and meta-analysis. Journal of the American Pharmacists Association: JAPhA 57(6):729–738.e10. 10. 10.MacBride-Stewart S, Guthrie B, Marwick C, Hurding S. National Therapeutic Indicators in Scotland and Financial Incentives. International Journal of Population Data Science 2019;4(3). 11. 11.Edelman EJ, Gordon KS, Crothers K, et al. Association of Prescribed Opioids With Increased Risk of Community-Acquired Pneumonia Among Patients With and Without HIV. JAMA Internal Medicine 2019;179(3):297–304. 12. 12.Wan Y, Li J, Shen L, et al. Enteric involvement in hospitalised patients with COVID-19 outside Wuhan. The Lancet Gastroenterology & Hepatology 2020; 13. 13.Luo H, Fan Q, Xiao S, Chen K. Changes in proton pump inhibitor prescribing trend over the past decade and pharmacists’ effect on prescribing practice at a tertiary hospital. BMC health services research 2018;18(1):537. 14. 14.Pasina L, Novella A, Elli C, Nobili A, Ianes A. Overuse of proton pump inhibitors in nursing homes: An Italian multicenter observational study. Pharmacoepidemiology and Drug Safety 2020;29(4):461–6. 15. 15.Ramirez E, Lei SH, Borobia AM, et al. Overuse of PPIs in patients at admission, during treatment, and at discharge in a tertiary Spanish hospital. Current Clinical Pharmacology 2010;5(4):288–97. 16. 16.Weiner SG, Malek SK, Price CN. The Opioid Crisis and Its Consequences. Transplantation 2017;101(4):678–81. 17. 17.de Abajo FJ, Rodríguez-Martín S, Lerma V, et al. Use of renin-angiotensin-aldosterone system inhibitors and risk of COVID-19 requiring admission to hospital: A case-population study. Lancet (London, England) 2020;395(10238):1705–14.