key: cord-1018641-p6sb2g77 authors: Murphy, N. M.; Dite, G. S.; Allman, R. title: Genetic associations with severe COVID-19 date: 2021-03-31 journal: nan DOI: 10.1101/2021.03.29.21254509 sha: ff974e9df762b68bb8fbed5eb770755eeb2a6e71 doc_id: 1018641 cord_uid: p6sb2g77 Identification of host genetic factors that predispose individuals to severe COVID-19 is important, not only for understanding the disease and guiding the development of treatments, but also for risk prediction when combined to form a polygenic risk score (PRS). Using population controls, Pairo-Castineira et al. identified 12 SNPs (a panel of 8 SNPs and a panel of 6 SNPs, with two SNPs in both panels) associated with severe COVID-19. Using controls with asymptomatic or mild COVID-19, we were able to replicate the association with severe COVID-19 for only three of their SNPs and found marginal evidence for an association for one other. When combined as an 8-SNP PRS and a 6-SNP PRS, we found no evidence of association with severe COVID-19. The difference in our results and the results of Pairo-Castineira et al. might be the choice of controls: population controls vs controls with asymptomatic or mild COVID-19. Pairo-Castineira et al. 1 identified 12 single nucleotide polymorphisms (SNPs) that were associated with critical illness in COVID-19, in particular a panel of 8 SNPs (associated with risk of intensive care admission) from independent genome-wide significant regions and a panel of 6 SNPs (associated with hospitalization) from a meta-analysis of overlapping SNPs between GenOMICC, Host Genetics Initiative and 23andMe studies (two of these SNPs were also in the panel of 8 SNPs). These SNP panels were identified and validated using population controls. Using the UK Biobank, [2] [3] [4] we evaluated the ability of the two SNP panels to predict severe COVID-19should a person become infectedwhen combined as PRSs. We downloaded the UK Biobank COVID-19 results file on 8 January 2021 and identified 8672 active participants with a positive SARS-CoV-2 test result, 8374 of whom had SNP data available. As we did previously, 5 we used source of test result as a proxy for disease severity, with 2288 (27.3%) participant results coming from an inpatient setting (cases) and 6086 (72.7%) participant results coming from an outpatient setting (controls). We used the estimates of the odds ratio (OR) per effect allele provided in Pairo-Castineira et al. 1 to construct the 8-SNP and 6-SNP PRS (the meta-analysis ORs for the 6-SNP PRS). For the SNPs in their Table 1 and the two overlapping SNPs from their Table 2, we used the GenoMICC risk allele frequencies provided, and for the remaining SNPs in their Table 2 , we used allele frequencies from SNPnexus. 6 For the construction of the PRS, 7 we assumed independent and additive risks on the log OR scale. For each SNP, we calculated the is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted March 31, 2021. ; https://doi.org/10.1101/2021.03.29.21254509 doi: medRxiv preprint Next, for each SNP, adjusted risks (with a population average risk equal to 1) were calculated as: where N is the number of effect alleles. The PRS was then calculated as the product of the adjusted risk values for each of the SNPs. We first used logistic regression to estimate the unadjusted per allele effects for each SNP (see Table 1 ). Of the 12 SNPs, three (rs71325088, rs73064425, and rs9380142) were associated with severe disease and one (rs13050728) was marginally statistically significant using a nominal statistical significance of P<0.05. These SNPs all had smaller effect sizes in our analysis compared with those reported in Pairo-Castineira et al. 1 Using population controls rather than people who are SARS-CoV-2 positive but are asymptomatic or have mild disease may have done more than bias associations towards the null, as suggested by Pairo-Castineira et al. 1 If the associations seen using population controls were biased towards the null, we would expect the effects of loci identified with population controls to be stronger when using asymptomatic or mild SARS-CoV-2 positive controls. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint Figure 1A and B). Interestingly, three of these regions (3p21.31, 9q34.2, and 12q24.13) overlapped with the Manhattan plot for all COVID-19 (i.e. infection): the COVID-19 vs population controls metaanalysis (see Figure 1B and C). The 3p21.31 locus has been identified previously in studies that used population controls, 10,11 and rs71325088 and rs73064425 SNP identified by Pairo-Castineira et al. 1 are both in 3p21.31. We also see an association for those two SNPs, albeit smaller in magnitude, using SARS-CoV-2 positive controls. The use of population controls might mean that Pairo-Castineira et al. 1 SNPs associated with becoming infected with SARS-CoV-2 rather than SNPs associated with severe COVID-19, a limitation that has been identified by other authors using population controls. 12 Researchers who use the results of large studies with population controls to inform their investigation of severe COVID-19 may therefore be wasting their time and resources. We urge care in the interpretation of results from large-scale studies that use population controls and invite discussion on this issue. is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint is the author/funder, who has granted medRxiv a license to display the preprint in (which was not certified by peer review) preprint The copyright holder for this this version posted March 31, 2021. ; https://doi.org/10.1101/2021.03.29.21254509 doi: medRxiv preprint Genetic mechanisms of critical illness in Covid-19 The UK Biobank resource with deep phenotyping and genomic data UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age UK Biobank makes health data available to tackle COVID-19 An integrated clinical and genetic model for predicting risk of severe COVID-19: A population-based case-control study SNPnexus: a web server for functional annotation of human genome sequence variation (2020 update) Assessment of clinical validity of a breast cancer risk model combining genetic and clinical information The COVID-19 Host Genetics Initiative, a global initiative to elucidate the role of host genetic factors in susceptibility and severity of the SARS-CoV-2 virus pandemic COVID-19 Host Genetics Initiative. COVID19hg GWAS meta-analyses round 5 The UK Biobank has Research Tissue Bank approval (REC #11/NW/0382) that covers analysis of data by approved researchers. All participants provided written informed consent to the UK Biobank before data collection began. This research has been conducted using the UK Biobank resource under Application Number 47401. GSD, NMM, and RA are employees of Genetic Technologies Limited. Genetic Technologies Limited had no role in the conceptualization, design, data analysis, decision to publish or preparation of the manuscript. All authors were involved in the conceptual development of this study. NMM and RA examined the Host Genetics Initiative meta-analysis results. GSD undertook data management and conducted the statistical analyses. All authors contributed to writing the first draft and reviewed and edited the manuscript. All authors have approved the final