Analysis of associations between polygenic risk score and COVID-19 severity in Russian population using low-pass genome sequencing ================================================================================================================================== * Arina V. Nostaeva * Valentin S. Shimansky * Svetlana V. Apalko * Ivan A. Kuznetsov * Natalya N. Sushentseva * Oleg S. Popov * Yurii S. Aulchenko * Sergey G. Shcherbak ## ABSTRACT The course of COVID-19 is characterized by wide variability, with genetics playing a contributing role. Through large-scale genetic association studies, a significant link between genetic variants and disease severity was established. However, individual genetic variants identified thus far have shown modest effects, indicating a polygenic nature of this trait. To address this, a polygenic risk score (PRS) can be employed to aggregate the effects of multiple single nucleotide polymorphisms (SNPs) into a single number, allowing practical application to individuals within a population. In this work, we investigated the performance of a PRS model in the context of COVID-19 severity in 1085 Russian participants using low-coverage NGS sequencing. By developing a genome-wide PRS model based on summary statistics from the COVID-19 Host Genetics Initiative consortium, we demonstrated that the PRS, which incorporates information from over a million common genetic variants, can effectively identify individuals at significantly higher risk for severe COVID-19. The findings revealed that individuals in the top 10% of the PRS distribution had a markedly elevated risk of severe COVID-19, with an odds ratio (OR) of 2.2 (95% confidence interval (CI): 1.3-3.3, p-value=0.0001). Furthermore, incorporating the PRS into the prediction model significantly improved its accuracy compared to a model that solely relied on demographic information (p-value < 0.0001). This study highlights the potential of PRS as a valuable tool for identifying individuals at increased risk of severe COVID-19 based on their genetic profile. ## INTRODUCTION COVID-19, also known as coronavirus infection, is a contagious illness caused by the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The majority of individuals who contract the virus exhibit mild to moderate respiratory symptoms and can recover without requiring specific medical treatment. However, in certain cases, the disease can manifest in a severe form, requiring medical intervention [1,2]. Apart from external factors like virus characteristics and the effectiveness of public health, certain host-related factors such as older age, male gender, and pre-existing chronic diseases like hypertension and diabetes have been associated with susceptibility and severity of COVID-19 [3,4]. However, these risk factors alone cannot fully explain the wide variation observed in the disease severity. The course of COVID-19 can range from asymptomatic cases to acute respiratory distress and even death [5,6]. Early in the pandemic, it was noted that clinical factors alone were insufficient to account for the variability in disease severity across individuals, as severe cases were observed in young people without apparent predisposing factors, often within families [7]. This suggests that human genetics may play a role in the development of the disease. To gain insights into the aetiology of COVID-19, large-scale genetic association studies incorporating both rare and common genetic variants have employed various study designs. These investigations, along with subsequent follow-up studies, have expanded our understanding of the disease and provided potential avenues for its treatment. The COVID-19 Host Genetics Initiative (HGI) was established to identify genetic loci that impact the severity and susceptibility of COVID-19 [8]. This global effort aims to conduct a meta-analysis of multiple COVID-19 genome-wide association studies (GWAS), and to identify significant single nucleotide polymorphisms (SNPs) associated with infection, hospitalization, and mortality. Through comparisons of genomes of millions of COVID-19 patients and healthy individuals, these studies have implicated genetic variants in 13 loci associated with the severity of the disease [9]. The COVID-19-associated genetic variants could be related to the regulation of processes such as innate antiviral defence signalling, regulation of inflammatory organ damage, and upregulation of cell receptors [10]. Modulation of these pathways can impact susceptibility to infection and subsequent disease manifestation [11]. The effects of individual genetic variants identified so far are generally small, consistent with the polygenic architecture of this trait. An individual who tests negative for a specific risk variant may still have a high genetic risk due to other unmeasured genetic factors. While each single variant only explains a small portion of the risk for severe COVID-19, combining multiple genetic variants into a polygenic risk score (PRS) can offer a better prediction of the risk. PRS allows for the aggregation of the effects of multiple SNPs into a single score, which can be practically applied to individuals within a population [12]. Conventionally, a polygenic score is defined as a weighted linear combination of allele counts for SNPs observed in an individual’s genome. The PRS model consists of the weights of a set of SNPs, with the weights proportional to the estimated effects of the SNPs on the trait being studied [13]. Modern polygenic risk score models for human traits are typically estimated using summary statistics obtained from a genome-wide association meta-analysis (GWAMA) and a reference panel reflecting linkage disequilibrium (LD) in the population [13,14]. Over the past decade, PRS predictive performance has significantly improved due to larger GWAS sample sizes and advancements in methods for variable selection and effect estimation [15–24]. Polygenic scores can be utilized to rank individuals within a group based on their genetic predisposition to a disease [25–27]. This approach considers an individual’s genetic predisposition relative to the genetic predisposition of others in the same group, often expressed as a percentile representing where the individual’s PRS falls within the overall distribution of the group’s PRS. Several studies have explored the development and application of PRS using variants associated with COVID-19, revealing clear associations between PRS and the risk of severe disease. However, most PRS models have been applied to cohorts consisting predominantly of individuals of Western European ancestry [28–31]. Using 1,582 SARS-CoV-2 positive participants from the UK Biobank (1,018 with severe COVID-19 and 564 without severe COVID-19) and 64 SNPs for PRS calculation, Dite et al. developed and validated a clinical and genetic model for predicting the risk of severe COVID-19. Only 13% of participants from this study were non-white, and PRS alone had an area under the receiver operating characteristic curve (AUC) of 68% [31]. While one recent study included African and South Asian groups, the associations with COVID-19 outcomes were limited by applying a PRS based on only six SNPs [32]. Another study that considered non-Western European populations was constrained by its focus on a specific Russian cohort (athletes) and also included only six genetic polymorphisms in the PRS assessment [33]. The multi-ethnic approach implemented in a very recent paper using UK biobank data, allowed the applicability of PRS, based on 17 SNPs, to diverse populations, with the severity model performing well within Black and Asian cohorts [34,35]. Overall, results highlight the potential of PRS as a predictive marker for disease severity and provide further support for its application in risk stratification and personalized healthcare approaches in the context of COVID-19. Our study aimed to investigate the performance of the PRS model in the Russian population. The genomes of study participants (347 individuals with severe COVID-19 and 738 with moderate or without disease) were assessed using low-coverage (with mean depth x3) sequencing. Next, we developed a genome-wide PRS model for COVID-19 severity using the summary statistics from the COVID-19 host genetics initiative consortium. We demonstrated that PRS, incorporating information from more than a million common genetic variants, for COVID-19 severity can identify individuals with markedly elevated risk of severe COVID-19 course: OR=2.2 (95% confidence interval (CI): 1.3-3.3, p-value=0.0001) for individuals in the top 10% of the PRS distribution, and produces a significant improvement in the quality of prediction (p-value < 0.0001) compared to a model including only demographic information. ## RESULTS ### Participant Characteristics The participants of the study were the patients of the infectious disease department of the St. Petersburg State Health Care Institution “City Hospital No. 40, Kurortny District” who were admitted for treatment with coronavirus infection (confirmed by polymerase chain reaction), and healthy individuals. Healthy individuals are defined as people who did not require COVID-19 medical treatment at the time of the study (between April 2020 and March 2022). Table 1 shows the participants’ characteristics. Of the 1085 participants, 479 (44%) were female, with a mean age of 60 years, while for 606 (56%) males the mean age was equal to 56 years. Overall, 895 (82%) of all participants had COVID-19, of which 347 (39%) had severe COVID-19. Separation according to the severity of the disease was carried out according to the following criteria: the case group included 347 patients (214 men and 133 women, 63±15 years) with lung damage more than 50% (computed tomography (CT)-3 and CT-4), the control group included 738 patients (392 men and 346 women, 56±16 years), with lung damage less than 50% or without COVID-19. View this table: [Table 1.](http://medrxiv.org/content/early/2023/11/22/2023.11.20.23298335/T1) Table 1. The demographic and clinical characteristics of the participants. ### Low coverage sequencing and imputation For all samples low coverage sequencing, also called LP-WGS (low-pass whole genome sequencing), was performed with a depth of x3 genome coverage. LP-WGS is the type of WGS with genome coverage from x0.5 to x5 [36,37]. Due to low-coverage data often having poor genotype quality and resulting in high missing genotype rates, the genotype likelihoods (GL) need to be updated using a reference panel for more accurate genotype imputation [38,39]. We used a recent method called GLIMPSE, which performs haplotype phasing and genotype imputation for LP-WGS data through a Gibbs sampling procedure, leading to improved accuracy [38]. As a reference panel, we used the 1000 Genomes data [40]. To evaluate the efficiency of LP-WGS within PRS, we calculated PRS values for a sample (not included in the study population) sequenced 45 times (in each of the batches to control the quality of the sequencing process). The coefficient of variation (CV) for PRS values was equal to 0.5% demonstrating a good method performance. ### Overview of the approach Computation of PRSs requires both genotype data of target individuals and the PRS model. To build the PRS model we used summary statistics from the COVID-19 Host Genetics Initiative consortium (release 7) [8]. These results were obtained by the meta-analysis, which combined the results of 60 individual studies from 25 countries, with a total of 18,000 severe cases of COVID-19 and more than a million controls who either did not have a severe disease course or were not affected by COVID-19 during the study period. From the obtained summary statistics, we generated the PRS model using the Bayesian approach SBayesR with default parameters, implemented in the GCTB software [19,41,42]. Finally, we calculated individual PRS values using the PRS model (Fig. 1, Methods). ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/11/22/2023.11.20.23298335/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2023/11/22/2023.11.20.23298335/F1) Figure 1. Study design and workflow. The PRS model for COVID-19 severity was derived by combining summary association statistics from the COVID-19 Host Genetics Initiative consortium and a linkage disequilibrium reference panel of 50,000 individuals of European ancestry from the UK Biobank data set. As a computational algorithm, SBayesR was used, which is a Bayesian approach to calculate a posterior mean effect for all variants based on a prior (effect size in the previous GWAS) and subsequent shrinkage based on linkage disequilibrium. PRS model was restricted by a list of variants from HapMap3 and included about one million variants. ### Testing associations between PRS and severe COVID-19 We compared the distributions of PRS values between severe cases and the control group combining the milder forms of COVID-19 and healthy individuals (Fig. 2). Comparison of the mean PRS values, performed using Student’s t-test for two independent samples, showed significant difference (p-value=1.2e-06). ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/11/22/2023.11.20.23298335/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2023/11/22/2023.11.20.23298335/F2) Figure 2. Comparison of distributions of PRS values between the groups with and without severe COVID-19. a) Distribution of PRS in the groups with (Ncases=347) and without (Ncontrols=738) severe COVID-19. The x-axis represents PRS, with values scaled to a mean of 0 and a standard deviation of 1 (in the total sample) to facilitate interpretation. b) PRS values among cases versus controls. Within each box plot, the horizontal lines reflect the median, the top, and bottom of each box reflect the interquartile range, and the whiskers reflect the rest of the distribution, except for points that are determined to be “outliers”. Across the study population, PRS was normally distributed with the risk of severe COVID-19 rising in the right tail of the distribution, from 17% in the lowest decile to around 47% in the highest decile (Fig. 3). ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/11/22/2023.11.20.23298335/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2023/11/22/2023.11.20.23298335/F3) Figure 3. Prevalence of the severe COVID-19 according to PRS decile. All participants (N=1,085) were stratified by decile of the PRS distribution. The average prevalence in percent and 95% CI within each decile are displayed. Next, we found that 20% of the population with the highest PRS values had inherited a genetic predisposition that conferred OR=1.8 for severe COVID-19 (95% CI: 1.3-2.4, p-value=0.0003) in comparison with all others. The 10% of the population with the highest PRS values had an OR=2.2 for COVID-19 (95% CI: 1.3-3.3, p-value=0.0001). ### Evaluating the relationship between PRS and COVID-19 outcome The severe form of the disease is associated with an increased risk of death. To assess how much the risk of death is associated with an increased PRS value, next, we calculated the odds ratio (OR) for death between the group with the highest PRS values (10%) and all others. The resulting OR was 1.9 (95% CI: 1.1-3.1) with p-value = 0.018. Thus, in the group with the highest PRS values, the probability of death due to severe disease was almost doubled. We also compared the mean PRS values for groups with different COVID-19 outcomes (death vs no death or no disease). Results showed significant difference in mean PRS (p-value = 0.02, Fig. 4). ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/11/22/2023.11.20.23298335/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2023/11/22/2023.11.20.23298335/F4) Figure 4. Comparison of distributions of PRS values between the groups with and without death outcome. a) Distribution of PRS in the groups with (Ndeath=146) and without (Nno death=939) death outcome of COVID-19. The x-axis represents PRS, with values scaled to a mean of 0 and a standard deviation of 1 (in the total sample) to facilitate interpretation. b) PRS values among cases versus controls. Within each box plot, the horizontal lines reflect the median, the top, and bottom of each box reflect the interquartile range, and the whiskers reflect the rest of the distribution, except for points that are determined to be “outliers”. Next, we hypothesized that PRS for severe COVID-19 would be associated with a higher risk of severe COVID-19 in early age. In Kaplan–Meier analyses, which is a non-parametric statistic used to estimate the survival function from lifetime data, we divided the sample into three groups: 10% of all individuals with the highest PRS values, 10% of all individuals with the lowest PRS values and the rest (Fig. 5). The analysis showed that people from the group of high PRS values start to have increased risk in comparison with other groups already before the age of 40 years (p-value<3.7e-10 for the log rank test). For example, the average risk of a severe course, which is reached at the age of 60 years, in the group with the highest PRS is reached already at 50 years of age. ![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/11/22/2023.11.20.23298335/F5.medium.gif) [Figure 5.](http://medrxiv.org/content/early/2023/11/22/2023.11.20.23298335/F5) Figure 5. Association of PRS with Incident Severe COVID-19. All participants (N=1,085) were stratified, based on their PRS, into three categories: bottom decile, deciles 2–9, and top decile. Incident severe COVID-19 is plotted according to the PRS category. ### Receiver Operating Curve (ROC) analysis Next we analysed the association between PRS and severe COVID-19 using a multivariate logistic regression model adjusted for sex, age, and the first 10 principal components of genetic variation. In the adjusted model, a significant association between PRS and severe COVID-19 was found: OR=1.48 per standard deviation (95% CI: 1.3-1.7 with p-value < 0.0001). High values of PRS (the 10% of PRS distribution) were associated with the adjusted OR=2.7 (95% CI: 1.8-4.2, p-value < 0.0001). Analyses showed significant (p-value < 0.0001) improvements in AUC with the addition of PRS to the base model containing only the demographic predictors. Figure 6 shows that a model predicting the risk of severe COVID-19 had an AUC of 65% (95% CI: 62-69% by the formula given by Hanley and McNeil [43]) for a model excluding PRS, and it increased up to 67% (95% CI: 64-71%) when PRS was included. ![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/11/22/2023.11.20.23298335/F6.medium.gif) [Figure 6.](http://medrxiv.org/content/early/2023/11/22/2023.11.20.23298335/F6) Figure 6. The comparison of receiving operating curves for three logistic regression models. The full model included the demographic predictors (sex and age), PRS, and the first 10 principal components of genetic variation, while the covariates-only model excluded PRS. ## DISCUSSION In this study, we constructed a polygenic risk model for the prediction of the severity of COVID-19 and applied it to a target cohort of 1085 Russian participants. Comparing the distributions of PRS, incorporating information from one million common genetic variants, between the case and control groups revealed significant differences, indicating meaningful associations between PRS and corresponding COVID-19 outcomes. We also demonstrated the potential of LP-WGS with coverage less than x5 to be used for predicting the severity of COVID-19. Our main objective was to evaluate the predictive ability of PRS for COVID-19 severity. To achieve this, we developed a logistic regression model that included only demographic and technical covariates and the full model that also incorporated PRS. Comparison between these models demonstrated that incorporating PRS significantly enhanced the predictive accuracy. These findings align with a previous analysis made by Huang et al., where PRS values for severe COVID-19 were constructed by using 112 SNPs in 430,582 participants from the UK Biobank study [29]. In this work, AUC was calculated for a model including only demographic and clinical parameters, and for the full model, which also included PRS. For the first model, the AUC was 0.789, while in the full mode, the AUC was 0.794 (p-value=0.002 for increment in AUC). Higher overall prediction accuracy of the model could be attributed to utilisation of information on comorbidities (cardiovascular disease, hypertension, diabetes, chronic respiratory infections, asthma, and chronic obstructive pulmonary disease). Our PRS, based on approximately one million SNPs, gave a comparable improvement in AUC (0.5% vs 2%, respectively). The higher contribution of PRS in our case can be explained by the much larger number of genetic variants used but also by the absence of clinical factors in our model. Indeed, it is often observed that adding a predictor to a model having a high AUC improves it by an amount smaller than that that could be achieved by adding the same factor to a poorer model. Furthermore, stratifying individuals by PRS quantiles revealed an association with a distinctive risk of severe COVID-19 in resulting groups. The highest PRS categories generally exhibited higher (up to 2.2 for the top 10% PRS) odds ratios. This genetic basis for differences in disease severity among individuals also extended to the occurrence of fatalities due to COVID-19 (OR=1.9 for the top 10% PRS). These results demonstrate that polygenic risks can be employed to stratify patients and assess their risk of severe disease and mortality related to COVID-19. Additional survival analysis using the non-parametric Kaplan-Meier estimation revealed that the highest risk categories as defined by PRS not only exhibited higher odds ratios for COVID-19 severity but also experienced an earlier onset of increased risk compared to the mean- and low-risk categories. These findings provide insights into both the overall risk for severe COVID-19 and how the risk varies by age. These results can have practical implications for protecting individuals with a greater genetic vulnerability during potential future outbreaks. Targeted public health interventions, such as shielding measures, closer monitoring, protection from high-risk frontline work, and prioritization for vaccination, could help to mitigate the associated risk. Hospital-based applications of PRS could facilitate the screening of COVID-19 patients and aid in the early detection of severe disease [28]. Moreover, informing patients about their increased polygenic risk has shown some evidence of positive behavioural impact [44], potentially leading to a decrease in risk-taking behaviours and promoting better outcomes. A few limitations of the study should be noted. Firstly, despite the multi-ethnic and global nature of the HGI Release 7 meta-analysis, the participants were mostly of Western European descent, which may have affected the accuracy of the predictions in non-Western European populations [9]. Additionally, the lack of detailed clinical data led to the use of CT scans as a criterion for disease severity, which could have introduced some inaccuracy in the classification of the outcome measure for some participants. ## METHODS ### Study population and genetic sequencing As part of the COVID-19 study, biomaterial (blood) and clinical data from COVID-19 patients hospitalized in the infectious disease department of the St. Petersburg State Budgetary Healthcare Institution “City Hospital No. 40 of Kurortny District” were collected. In this work, low-coverage (x2-5) sequencing was performed for 1085 samples divided into 45 batch sizes. Low-coverage sequencing, also called LP-WGS (low-pass whole genome sequencing), is a low-cost, high-throughput DNA sequencing technology used to accurately detect genetic variation in the genomes of multiple species [45]. Using imputation algorithms, this technology provides high variant detection accuracy with very low sequence coverage. LP-WGS and subsequent imputation yield more accurate genotypes than imputation using genotyping data, allowing for increased power in GWAS studies and more accurate results in polygenic risk studies [46]. Prior to sequencing, preliminary analysis and quality control of the case database were performed, and preliminary analysis of samples from each batch was performed to exclude bias for any of the sample characteristics: age, sex, and case/control. Genome DNA isolation was performed with QIAcube, using QIAamp DNA Blood Mini Kit. DNA concentration is measured with Promega QuantiFluor dsDNA System. Library preparation was done using Roche KAPA HyperPlus Kit. Quality control electrophoresis was done on QIAxcel station using QIAxcel High Resolution Kit. Circularization was made with MGIEasy Circularization Kit. Sequencing was done on MGISEQ-2000 sequencing machine with DNBSEQ-G400RS High-throughput Sequencing Set (FCL PE150, 540 G). ### Variant calling, imputation, and quality control Quality analysis (FastQC) [47], alignment (BWA) [48], deduplication (samtools), and variant collation (bcftools) were performed for the reads obtained from sequencing [49]. Imputation of the resulting data was then performed using the GLIMPSE tool [38], which allows imputation of low-coverage sequencing data. To improve imputation quality, only bi-allelic sites were retained from the LP-WGS BAM data and processed with bcftools. Then iterative refinement of GL using the reference panels with segmentation size of 2 Mb with buffer size of 200 kb produced imputed dosages and multiple chunks within each chromosome were ligated. A panel of 1000 Genomes with high coverage [40], including high-quality SNV- and INDELs from over 3,000 samples, was used as a reference sample. Then, we filtered imputed variants by an imputation INFO score, where variants with score < 0.7 and a minor allele frequency < 0.1% were removed from the analysis [9,13]. We focused on the variants and individuals with a call rate of more than 90%. We also removed close relatives from the analysis. We used the KING-robust method to identify relatives [50]. Using a threshold (kinship > 0.125), we found pairs of first- and second-degree relatives. We restricted our analyses to a list of variants from HapMap3 [51], which are included in the PRS models. PLINK 1.90 software [52] was utilised for all genotype extraction and quality control. ### Establishing COVID-19 outcomes The severity of the course was divided according to the following criteria: the case group included samples with lung lesions greater than 50% (computed tomography (CT)-3 and CT-4), while the control group included all other samples. As a result, the case group included 347 patients (214 men and 133 women, 63±15 years) with lung damage more than 50% (computed tomography (CT)-3 and CT-4), the control group included 738 patients (392 men and 346 women, 56±16 years), with lung damage less than 50% or without COVID-19. ### Construction of PRS models The calculation of PRSs relies on both genotype data from the target individuals and a PRS model. To derive a PRS model, GWAS are used to estimate the effect sizes of SNPs [53]. However, the GWAS gives the marginal effect size for each SNP estimated by a regression model that ignores linkage disequilibrium (LD) structure. As a result, to construct a PRS model that incorporates multiple SNPs, the SNP effects must be re-estimated while accounting for LD structure. As the summary statistics, we used summary statistics from the COVID-19 Host Genetics Initiative consortium (release 7). These results were obtained by the meta-analysis, which combined the results of 60 individual studies from 25 countries. To re-weight the effect sizes, we used SBayesR, a software tool that has demonstrated superior performance compared to similar tools [19]. This tool re-weights the effects of each variant based on the marginal estimate of its effect size, statistical strength of association, the degree of correlation between the variant and other variants nearby, and tuning parameters. It also requires a GCTB-compatible LD matrix file based on individual-level data from a reference population, and for this analysis, we used a shrunk sparse GCTB LD matrix from 50,000 individuals of European ancestry in the UK Biobank dataset [41]. PRS values were calculated as a weighted sum of allele counts: ![Formula][1] with β*j* the re-weighted effect size of the *j**th* SNP, *G**ij* the genotype of the *j**th* SNP for *i**th* individual. PLINK 1.90 software [52] was utilised for PRS calculation. ### Statistical analysis and association testing Logistic regression of PRS categories against COVID-19 severity outcomes was then conducted using R [54] and Python3 [55], fully adjusted for covariates, such as sex and age. Data on comorbidities were not available for the majority of patients, as well as other clinical data, so parameters for these were not included in the model to cover as much data as possible. The first 10 principal genetic components (PCs) were also included as covariates to adjust for population genetic structures and avoid bias, as per current recommendations [13]. The discriminative power of models in identifying high-risk individuals was then assessed using receiver operating curve (ROC) analysis. Area under the ROC (AUC) was calculated for full models (consisting of covariates and PRS) and base models (covariates only). The confidence interval for AUC was calculated using the formula given by Hanley and McNeil [43]. Increment in AUC (ΔAUC) was reported based on the difference between the two models, reported as the discriminative or predictive power conferred by PRS. The permutation test for differences between classifiers was used to estimate the significance (p-value) of an increment in AUC. Once PRS was calculated, individuals were separately stratified into quintiles for susceptibility and severity PRS, then categorised into low genetic risk (decile 1, bottom 10% of cohort), intermediate risk (decile 2–9, middle 80%) and high risk (decile 10, top 10%) for each outcome. In each group, we estimated the cumulative hazard curve using the non-parametric method called the Kaplan-Meier estimator [56]. For each pair of groups, the log rank test was applied, which is the statistical test for comparing the survival distributions of two or more groups. ## Data Availability All data produced in the present study are available upon reasonable request to the authors. ## DATA AND CODE AVAILABILITY Personal genetic and clinical data are under restrictions and are available through collaboration with the St. Petersburg State Health Care Institution “City Hospital No. 40, Kurortny District” hospital. ## CONFLICT OF INTEREST YSA is a co-owner of PolyKnomics BV, a private organization providing services, research, and development in the field of computational and statistical genomics. YSA is currently a full-time employee of GSK. The other authors declare that they have no competing interests. ## ACKNOWLEDGEMENTS The authors are grateful to the study participants and the staff from the St. Petersburg State Health Care Institution “City Hospital No. 40, Kurortny District” hospital. The authors would like to thank all authors of the included studies for their valuable contributions to data collection. This work was supported by Saint Petersburg State University, project ID: 94029859. ## Footnotes * Title updated * Received November 20, 2023. * Revision received November 22, 2023. * Accepted November 22, 2023. * © 2023, 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.Bose S, Adapa S, Aeddula NR, Roy S, Nandikanti D, Vupadhyayula PM, et al. Medical Management of COVID-19: Evidence and Experience. J Clin Med Res. 2020;12: 329–343. 2. 2.Bevova MR, Netesov SV, Aulchenko YS. The New Coronavirus COVID-19 Infection. Mol Gen Microbiol Virol. 2020;35: 53–60. 3. 3.COVID-19 National Preparedness Collaborators. Pandemic preparedness and COVID-19: an exploratory analysis of infection and fatality rates, and contextual factors associated with preparedness in 177 countries, from Jan 1, 2020, to Sept 30, 2021. Lancet. 2022;399: 1489–1512. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 4. 4.Biswas M, Rahaman S, Biswas TK, Haque Z, Ibrahim B. Association of Sex, Age, and Comorbidities with Mortality in COVID-19 Patients: A Systematic Review and Meta-Analysis. Intervirology. 2020; 1–12. 5. 5.Wang Y, Wang Y, Chen Y, Qin Q. Unique epidemiological and clinical features of the emerging 2019 novel coronavirus pneumonia (COVID-19) implicate special control measures. J Med Virol. 2020;92: 568–576. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 6. 6.Fricke-Galindo I, Falfán-Valencia R. Genetics Insight for COVID-19 Susceptibility and Severity: A Review. Front Immunol. 2021;12: 622176. 7. 7.Yousefzadegan S, Rezaei N. Case Report: Death due to COVID-19 in Three Brothers. Am J Trop Med Hyg. 2020;102: 1203–1204. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4269/ajtmh.20-0240&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 8. 8.COVID-19 Host Genetics Initiative. 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. Eur J Hum Genet. 2020;28: 715–718. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41431-020&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 9. 9.Mapping the human genetic architecture of COVID-19. Nature. 2021;600: 472–477. 10. 10.Pairo-Castineira E, Clohisey S, Klaric L, Bretherick AD, Rawlik K, Pasko D, et al. Genetic mechanisms of critical illness in COVID-19. Nature. 2021;591: 92–98. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-03065-y&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33307546&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 11. 11.Velavan TP, Pallerla SR, Rüter J, Augustin Y, Kremsner PG, Krishna S, et al. Host genetic factors determining COVID-19 susceptibility and severity. EBioMedicine. 2021;72: 103629. 12. 12.Dudbridge F. Power and predictive accuracy of polygenic risk scores. PLoS Genet. 2013;9: e1003348. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1003348&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23555274&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 13. 13.Choi SW, Mak TS-H, O’Reilly PF. Tutorial: a guide to performing polygenic risk score analyses. Nat Protoc. 2020;15: 2759–2772. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 14. 14.de los Campos G, Gianola D, Allison DB. Predicting genetic predisposition in humans: the promise of whole-genome markers. Nat Rev Genet. 2010;11: 880–886. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrg2898&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21045869&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 15. 15.Vilhjálmsson BJ, Yang J, Finucane HK, Gusev A, Lindström S, Ripke S, et al. Modeling Linkage Disequilibrium Increases Accuracy of Polygenic Risk Scores. Am J Hum Genet. 2015;97: 576–592. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2015.09.001&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26430803&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 16. 16.Zhu X, Stephens M. BAYESIAN LARGE-SCALE MULTIPLE REGRESSION WITH SUMMARY STATISTICS FROM GENOME-WIDE ASSOCIATION STUDIES. Ann Appl Stat. 2017;11: 1561–1592. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1214/17-AOAS1046&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29399241&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 17. 17.Khera AV, Chaffin M, Aragam KG, Haas ME, Roselli C, Choi SH, et al. Genome-wide polygenic scores for common diseases identify individuals with risk equivalent to monogenic mutations. Nat Genet. 2018;50: 1219–1224. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018–0183-z&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30104762&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 18. 18.Choi SW, O’Reilly PF. PRSice-2: Polygenic Risk Score software for biobank-scale data. Gigascience. 2019;8. doi:10.1093/gigascience/giz082 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/gigascience/giz082&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31307061&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 19. 19.Lloyd-Jones LR, Zeng J, Sidorenko J, Yengo L, Moser G, Kemper KE, et al. Improved polygenic prediction by Bayesian multiple regression on summary statistics. Nat Commun. 2019;10: 5086. 20. 20.Li R, Chang C, Tanigawa Y, Narasimhan B, Hastie T, Tibshirani R, et al. Fast numerical optimization for genome sequencing data in population biobanks. Bioinformatics. 2021;37: 4148–4155. 21. 21.Privé F, Arbel J, Vilhjálmsson BJ. LDpred2: better, faster, stronger. Bioinformatics. 2021;36: 5424–5431. 22. 22.Ojavee SE, Kousathanas A, Trejo Banos D, Orliac EJ, Patxot M, Läll K, et al. Genomic architecture and prediction of censored time-to-event phenotypes with a Bayesian genome-wide analysis. Nat Commun. 2021;12: 2337. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-22538-w&link_type=DOI) 23. 23.Wand H, Lambert SA, Tamburro C, Iacocca MA, O’Sullivan JW, Sillari C, et al. Improving reporting standards for polygenic scores in risk prediction studies. Nature. 2021;591: 211–219. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-021-03243-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33692554&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 24. 24.Lambert SA, Gil L, Jupp S, Ritchie SC, Xu Y, Buniello A, et al. The Polygenic Score Catalog as an open database for reproducibility and systematic evaluation. Nat Genet. 2021;53: 420–425. 25. 25.Mavaddat N, Michailidou K, Dennis J, Lush M, Fachal L, Lee A, et al. Polygenic Risk Scores for Prediction of Breast Cancer and Breast Cancer Subtypes. Am J Hum Genet. 2019;104: 21–34. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2018.11.002&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30554720&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 26. 26.Musliner KL, Mortensen PB, McGrath JJ, Suppli NP, Hougaard DM, Bybjerg-Grauholm J, et al. Association of Polygenic Liabilities for Major Depression, Bipolar Disorder, and Schizophrenia With Risk for Depression in the Danish Population. JAMA Psychiatry. 2019;76: 516–525. 27. 27.Cupido AJ, Tromp TR, Hovingh GK. The clinical applicability of polygenic risk scores for LDL-cholesterol: considerations, current evidence and future perspectives. Curr Opin Lipidol. 2021;32: 112–116. 28. 28.Lambert SA, Abraham G, Inouye M. Towards clinical utility of polygenic risk scores. Hum Mol Genet. 2019;28: R133–R142. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddz187&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 29. 29.Huang Q-M, Zhang P-D, Li Z-H, Zhou J-M, Liu D, Zhang X-R, et al. Genetic Risk and Chronic Obstructive Pulmonary Disease Independently Predict the Risk of Incident Severe COVID-19. Ann Am Thorac Soc. 2022;19: 58–65. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 30. 30.Dite GS, Murphy NM, Allman R. Development and validation of a clinical and genetic model for predicting risk of severe COVID-19. Epidemiol Infect. 2021;149: e162. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S095026882100145X&link_type=DOI) 31. 31.Dite GS, Murphy NM, Allman R. An integrated clinical and genetic model for predicting risk of severe COVID-19: A population-based case-control study. PLoS One. 2021;16: e0247205. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0247205&link_type=DOI) 32. 32.Horowitz JE, Kosmicki JA, Damask A, Sharma D, Roberts GHL, Justice AE, et al. Genome-wide analysis provides genetic evidence that ACE2 influences COVID-19 risk and yields risk scores associated with severe disease. Nat Genet. 2022;54: 382–392. 33. 33.Ahmetov II, Borisov OV, Semenova EA, Andryushchenko ON, Andryushchenko LB, Generozov EV, et al. Team sport, power, and combat athletes are at high genetic risk for coronavirus disease-2019 severity. J Sport Health Sci. 2020;9: 430–431. 34. 34.Farooqi R, Kooner JS, Zhang W. Associations between polygenic risk score and covid-19 susceptibility and severity across ethnic groups: UK Biobank analysis. BMC Med Genomics. 2023;16: 150. 35. 35.Sudlow C, Gallacher J, Allen N, Beral V, Burton P, Danesh J, et al. UK biobank: an open access resource for identifying the causes of a wide range of complex diseases of middle and old age. PLoS Med. 2015;12: e1001779. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1001779&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25826379&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 36. 36.Chaubey A, Shenoy S, Mathur A, Ma Z, Valencia CA, Reddy Nallamilli BR, et al. Low-Pass Genome Sequencing: Validation and Diagnostic Utility from 409 Clinical Cases of Low-Pass Genome Sequencing for the Detection of Copy Number Variants to Replace Constitutional Microarray. J Mol Diagn. 2020;22: 823–840. 37. 37.Li JH, Mazur CA, Berisa T, Pickrell JK. Low-pass sequencing increases the power of GWAS and decreases measurement error of polygenic risk scores compared to genotyping arrays. Genome Res. 2021;31: 529–537. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjg6IjMxLzQvNTI5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTEvMjIvMjAyMy4xMS4yMC4yMzI5ODMzNS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 38. 38.Rubinacci S, Ribeiro DM, Hofmeister RJ, Delaneau O. Efficient phasing and imputation of low-coverage sequencing data using large reference panels. Nat Genet. 2021;53: 120–126. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-020-00756-0&link_type=DOI) 39. 39.Hui R, D’Atanasio E, Cassidy LM, Scheib CL, Kivisild T. Evaluating genotype imputation pipeline for ultra-low coverage ancient genomes. Sci Rep. 2020;10: 18542. 40. 40.1000 Genomes Project Consortium, Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, et al. A global reference for human genetic variation. Nature. 2015;526: 68–74. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature15393&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26432245&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 41. 41.Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562: 203–209. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-018-0579-z&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30305743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 42. 42.Jiang L, Zheng Z, Qi T, Kemper KE, Wray NR, Visscher PM, et al. A resource-efficient tool for mixed model association analysis of large-scale data. Nat Genet. 2019;51: 1749–1755. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588–019–0530–8&link_type=DOI) 43. 43.Hanley JA, McNeil BJ. The meaning and use of the area under a receiver operating characteristic (ROC) curve. Radiology. 1982;143: 29–36. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1148/radiology.143.1.7063747&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7063747&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1982NG95400006&link_type=ISI) 44. 44.Frieser MJ, Wilson S, Vrieze S. Behavioral impact of return of genetic test results for complex disease: Systematic review and meta-analysis. Health Psychol. 2018;37: 1134–1144. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1037/hea0000683&link_type=DOI) 45. 45.Alex Buerkle C, Gompert Z. Population genomics based on low coverage sequencing: how low should we go? Mol Ecol. 2013;22: 3028–3035. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/mec.12105&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23174005&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 46. 46.Homburger JR, Neben CL, Mishne G, Zhou AY, Kathiresan S, Khera AV. Low coverage whole genome sequencing enables accurate assessment of common variants and calculation of genome-wide polygenic scores. Genome Med. 2019;11: 74. 47. 47.Andrews S. FastQC: A Quality Control Tool for High Throughput Sequence Data. In: FastQC [Internet]. 2010 [cited 2010]. Available: [http://www.bioinformatics.babraham.ac.uk/projects/fastqc/](http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) 48. 48.Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25: 1754–1760. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btp324&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19451168&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000267665900006&link_type=ISI) 49. 49.Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25: 2078–2079. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btp352&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19505943&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268808600014&link_type=ISI) 50. 50.Manichaikul A, Mychaleckyj JC, Rich SS, Daly K, Sale M, Chen W-M. Robust relationship inference in genome-wide association studies. Bioinformatics. 2010;26: 2867–2873. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq559&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20926424&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000283919800010&link_type=ISI) 51. 51.HapMap Project. HapMap Project. In: HapMap 3 [Internet]. Available: [https://www.sanger.ac.uk/data/hapmap-3/](https://www.sanger.ac.uk/data/hapmap-3/) 52. 52.Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81: 559–575. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/519795&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17701901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) 53. 53.Uffelmann E, Huang QQ, Munung NS, de Vries J, Okada Y, Martin AR, et al. Genome-wide association studies. Nature Reviews Methods Primers. 2021;1: 1–21. 54. 54.Team RDC. R: A language and environment for statistical computing. (No Title). 2010 [cited 20 Oct 2023]. Available: [https://cir.nii.ac.jp/crid/1370294721063650048](https://cir.nii.ac.jp/crid/1370294721063650048) 55. 55.Van Rossum G, Drake FL. Python 3 Reference Manual; CreateSpace: Scotts Valley, CA, USA, 2009. Google Scholar. 56. 56.Goel MK, Khanna P, Kishore J. Understanding survival analysis: Kaplan-Meier estimate. Int J Ayurveda Res. 2010;1: 274–278. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4103/0974-7788.76794&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21455458&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F11%2F22%2F2023.11.20.23298335.atom) [1]: /embed/graphic-8.gif