GWAS and polygenic risk score of severe COVID-19 in Eastern Europe

Background: COVID-19 disease has infected more than 772 million people, leading to 7 million deaths. Although the severe course of COVID-19 can be prevented using appropriate treatments, effective interventions require a thorough research of the genetic factors involved in its pathogenesis.

Methods: We conducted a genome-wide association study (GWAS) on 7,124 individuals (comprising 6,400 controls who had mild to moderate COVID-19 and 724 cases with severe COVID-19). The inclusion criteria were acute respiratory distress syndrome (ARDS), acute respiratory failure (ARF) requiring respiratory support, or CT scans indicative of severe COVID-19 infection without any competing diseases. We also developed a polygenic risk score (PRS) model to identify individuals at high risk.

Results: We identified two genome-wide significant loci (P-value <5 × 10⁻⁸) and one locus with approximately genome-wide significance (P-value = 5.92 × 10⁻⁸-6.15 × 10⁻⁸). The most genome-wide significant variants were located in the leucine zipper transcription factor like 1 (LZTFL1) gene, which has been highlighted in several previous GWAS studies. Our PRS model results indicated that individuals in the top 10% group of the PRS had twice the risk of severe course of the disease compared to those at median risk [odds ratio = 2.18 (1.66, 2.86), P-value = 8.9 × 10⁻⁹].

Conclusion: We conducted one of the largest studies to date on the genetics of severe COVID-19 in an Eastern European cohort. Our results are consistent with previous research and will guide further epidemiologic studies on host genetics, as well as for the development of targeted treatments.

Background

Coronavirus infection, or COVID-19, caused by the SARS-CoV-2 virus, has resulted in one of the largest pandemics in human history (1). Epidemiological data on this disease emphasizes the remarkable heterogeneity in the course of the disease, ranging from completely asymptomatic cases to ICU hospitalizations with ventilator support and even fatal outcomes (2–5).

COVID-19 is associated with several comorbidities and non-genetic risk factors that increase the likelihood of developing the disease and experiencing severe progression, including respiratory failure. These risk factors include older age, male gender, and certain medical conditions such as cardiovascular disease, diabetes, obesity, chronic respiratory and kidney diseases, immunodeficiency, and neurological disorders (6–10).

The observed variability in susceptibility to and course of coronavirus infection suggests that both nongenetic risk factors and genetic variants may contribute to the clinical heterogeneity of COVID-19 in the population. Therefore, the COVID-19 Host Genetic Initiative (HGI) conducted a large study on a multiethnic sample of more than 49,500 COVID-19 patients from 46 studies in 19 countries (11). The study identified 13 genome-wide significant loci that are associated with SARS-CoV-2 infection or severe manifestations of COVID-19. Along with other genome-wide association studies (GWAS), several key genetic variants associated with severe COVID-19 outcomes were identified, which included those related to the function of the SARS-CoV-2 receptor (ACE2, ABO, TMPRSS2, and SLC6A20) (12–15) and immune response to the virus (HLA-region) (16, 17).

As epidemiological data and associated genetic variants for COVID-19 disease continue to accumulate, the combined assessment of disease severity for clinical implications remains an area of ongoing research. It is hypothesized that a combination of different gene variants may determine the severity of COVID-19 disease course. In other words, the complex polymorphic pathogenesis of coronavirus infection suggests a polygenic architecture.

Horowitz et al. (13) used the data from the COVID-19 HGI to establish a polygenic risk score (PRS) model. They demonstrated that individuals in the top 10% of the COVID-19 PRS among Europeans (n = 44.958) had a 1.38-fold increased risk of hospitalization and 1.58-fold increased risk of severe disease. Similarly, Farooqi et al. (18) reported a 1.57-fold higher risk of severe COVID-19 in high-risk patients compared to low-risk individuals. Another study by Crossfield et al. (19) involving 9,560 UK Biobank participants showed an adjusted odds ratio (OR) of 1.32, [95% confidence interval (CI): 1.11–1.58] for the highest PRS quintile compared with the lowest one. Nostaeva et al. (20) applied the same HGI statistics to a small cohort of Russian patients (1,085 participants, 347 individuals with severe COVID-19, and 738 with moderate or without disease) low-pass whole genome sequencing (LP-WGS). They found that more than one million genetic variants can be included in calculating polygenic scores to stratify patients by the risk of severe COVID-19. Individuals in the top 10% of the PRS distribution had over a two-fold increased risk of severe COVID-19 (odds ratio, OR: 2.2; 95% CI: 1.3–3.3, P-value = 0.0001).

Since coronavirus infection can develop rapidly in just over a week, identifying individuals at risk of developing severe COVID-19 through genetic variants may help identify targeting agents for investigating appropriate therapeutic interventions. Despite large-scale vaccination programs, optimal treatment selection remains a topical challenge. Therefore, the generation of a PRS model and patient profiling based on the risk of severe COVID-19 course can serve as a valuable tool for the healthcare industry.

In this research, we investigate the genetic factors that contribute to the severity of COVID-19 within the Eastern European population, which is underrepresented in many studies. By utilizing GWAS, we aim to identify genetic variants associated with COVID-19 and explore how these variants differ from those found in other populations. Additionally, we develop a PRS model to identify individuals at a high risk of experiencing severe COVID-19 outcomes.

Materials and methods

Study cohort

We analyzed the genetic data of 787 individuals, with 691 from the N.V. Sklifosovsky Research Institute for Emergency Medicine and 96 from the Ufa Federal Research Center of the Russian Academy of Sciences (UFRC RAS). Biomaterials, specifically blood or saliva, were collected from individuals who had a history of severe COVID-19.

Patients were classified as having an extremely severe course of COVID-19 if they met at least one of the following inclusion criteria: acute respiratory distress syndrome (ARDS); acute respiratory failure (ARF) requiring respiratory support, which could include high-flow non-invasive or invasive ventilation; lung changes on CT scans indicative of viral damage, such as significant or subtotal lesion volume (CT grade 4); or a clinical presentation consistent with ARDS. In addition, two mandatory criteria had to be met: a confirmed COVID-19 infection with the virus identified (ICD-10 code U07.1) and the absence of other diseases that could potentially worsen the patient's condition, such as acute myocardial infarction, exacerbation of bronchial asthma or chronic obstructive pulmonary disease (COPD), or decompensation of chronic heart failure.

This research was approved by the Genotek Ethics Committee (protocol No15 “GWAS of severe COVID-19 in the Russian population”) and performed in accordance with the Declaration of Helsinki. The individuals who were included in our analysis provided informed consent for their data to be used for research purposes and responded to an online questionnaire.

Genotek customer data were added to the study cohort for further analysis as controls (N = 6,400). The participants of the study were selected based on responses to a questionnaire, adhering to the following criteria: consent was provided for the use of anonymized data in scientific research and individuals aged 40 years or above should self-report a COVID-19 diagnosis, confirmed by antibody tests, PCR tests, CT scans, or a physician's diagnosis. These individuals may have experienced mild symptoms such as general fatigue, cough, and loss of smell or taste, along with other non-severe symptoms, without requiring hospitalization.

The number of cases was predetermined as we have obtained their data from hospitals. Although we had the option to adjust the number of cases by modifying the age threshold, we decided to maintain a ratio of ~10 controls per case to enhance the statistical power of our study. The rationale for selecting this sample size has been discussed in a study by Katki et al. (21).

Genotyping

DNA extraction and genotyping were performed on saliva samples that were genotyped on Illumina Infinium Global Screening Array v.3 microarrays [~650,000 single nucleotide polymorphisms (SNPs)]. All samples in the cohort were processed in batches with 192–768 samples per batch using the Genotek microarray data processing pipeline. The pipeline involves variant detection based on iaap-cli 1.1.0 and bcftools +gtc2vcf plugin 1.11, followed by subsequent filtering and analysis. GenomeStudio software (Illumina, San Diego, CA) and manually created cluster files were used to cluster the raw signals and call the genotypes. SNPs with a call rate of < 0.9 within the batch were removed.

Quality control and data preparation

Two-stage sequential filtering was performed based on the number of variants with undetermined genotypes. First, we filtered out the genetic variants with undetermined genotypes in more than 20% of the samples because the quality of detection of these variants was likely to be low, potentially leading to inaccurate conclusions in subsequent stages of the analysis. Subsequently, samples with undetermined genotypes in more than 20% of positions were excluded because the quality of collection, preparation, or analysis of these samples was likely to be low, which may lead to inaccurate conclusions in subsequent stages of the analysis. Following this step, filtering was repeated for positions and samples with a threshold of 2%.

After filtering for variant and sample quality, heterozygosity analysis was performed. The samples with abnormal heterozygosity were filtered using PLINK 1.9. We excluded the samples in which the observed heterozygosity deviated by more than 3 standard deviations from the mean. Heterozygosity estimation was performed on the cohort after filtering for genetic variants in linkage disequilibrium, using a search window of 50 SNPs, with five SNPs to shift the window at the end of each step, and an r² value of SNPs < 0.2.

Subsequently, the human genotype was determined at positions not represented on the microarray by applying linkage disequilibrium (LD). This procedure was performed using the Beagle 5.1 program (22) using two reference panels: 1000 Genomes (23) and the Haplotype Reference Consortium (23, 24). Only those positions that achieved a high-quality metric for defining human genotypes with DR2 < 0.7 were used in further analysis. Multi-allelic substitutions were excluded from further analysis.

In the following step, the positions on sex chromosomes and mitochondrial DNA (mtDNA) were excluded. In addition, positions that violated the Hardy–Weinberg equilibrium were filtered out; specifically, we removed the positions with a significant difference between the observed genotype frequencies and the expected frequencies according to the Hardy–Weinberg test (P-value < 1 × 10⁻⁵). Finally, the positions with a low minor allele frequency (MAF < 0.05) were excluded from further analysis.

Identification of close relatives within the study cohort was performed using the PRIMUS (Rapid Reconstruction of Pedigrees from Genome-wide Estimates of Identity by Descent) program (25). Pairs with a PI_HAT score of >0.15 were considered as related. The cohort was filtered to ensure that it contained no pairs of relatives.

Genome-wide association study (GWAS) and heritability

We performed population stratification and filtered outliers before conducting the GWAS analysis. Initially, the principal component analysis (PCA) algorithm (MultiDimensional Scaling) was applied for dimensionality reduction. Positions filtered by non-equilibrium coupling were used, considering a search window of 50 SNPs, with five SNPs to shift the window at the end of each step and an r² value of SNPs < 0.2). Based on the values of the first and second components, clustering was performed using the DBSCAN algorithm (26). The final PCA plot is shown in Supplementary Figure S1. After selecting the largest cluster, outliers and samples not from this cluster were excluded. The top 20 components were subsequently used as covariates to account for population stratification.

The GWAS analysis was performed using the PLINK 2 program. A logistic regression model was used, and 20 components of PCA and gender were included as covariates.

The statistical package ldsc (https://github.com/bulik/ldsc) was used to estimate SNP heritability.

Polygenic risk score

To construct the PRS model, summary statistics from the International Consortium COVID-19 Host Genetics Initiative (27) were employed (A2 phenotype, ALL_leave_23andme cohort, release 7). Variants with complementary alleles and variants with repeated rs_id were removed during data preprocessing. The PRS was trained using the LDPred2 tool (28) on the data from the entire cohort of 7,124 individuals. The scores were adjusted for sex and the first 20 components of PCA coordinates.

Results

To account for population stratification, we constructed a PCA plot with two principal components (Supplementary Figure S1) and performed clustering to detect minor populations and outliers. The largest cluster (number 1) was selected for further analysis and all other samples were considered outliers.

We performed a GWAS analysis on 7,124 individuals (with 53% being female individuals) from the selected cluster. Among them, 6,400 individuals were controls (i.e., they were aged over 40 years and had recovered from COVID-19 without experiencing a severe outcome), and the remaining 724 patients were cases (i.e., they experienced a severe course of COVID-19). The mean age was 64.77 years for cases and 49.53 years for controls. The characteristics of the final cohort are summarized in Table 1.

Table 1

Table 1. Characteristics of the study population.

The GWAS analysis revealed several genome-significant loci (P-value < 5 × 10⁻⁸), (Table 2). Figure 1 shows the Manhattan plot with the GWAS results.

Table 2

Table 2. Top GWAS results (P-value <10⁻⁸).

Figure 1

Figure 1. Manhattan plot with the GWAS results for the severe COVID-19 phenotype (P-value < 5 × 10⁻⁸).

Particularly, the most genome-wide significant loci were identified in the leucine zipper transcription factor like 1 (LZTFL1) gene: rs35624553 (OR = 1.58, 95% CI: 1.34–1.85, P-value = 3.26 × 10⁻⁸) and rs10490770 (OR = 1.57, 95% CI: 0.082–1.34, P-value = 3.51 × 10⁻⁸). These findings were previously mentioned in the summary statistics of the COVID-19 HGI relevant to our phenotype. Three intergenic variants on chromosome 3 (rs17713054, rs13078854, and rs71325088) had the same association with severe forms of COVID-19 as observed with rs10490770.

Furthermore, we compared the effects of SNPs identified by the COVID-19 HGI (11), which are associated with critical illness (similar to our phenotype) and the current GWAS results. We could not compare the effect of rs912805253 as it is absent in hg19. The remaining SNPs exhibited poor approximation with HGI results; however, the intergenic variants highlighted in the HGI summary statistics (rs1819040 and rs74956615) demonstrated the effects that were comparable to our results (rs1819040: OR_hgi = 0.906, OR_genotek = 0.920, rs74956615: OR_hgi = 1.434, OR_genotek = 1.461). Figure 2 illustrates these results.

Figure 2

Figure 2. Scatter plot of comparison of odds ratios (OR) from the top associated SNPs from the COVID-19 Host Genetics Initiative with our results (OR_hgi indicates results from COVID-19 HGI GWAS and OR_gwas indicates results from the current study). Horizontal and vertical bars represent 95% confidence intervals.

The heritability h2 coefficient was found to be 0.05 ± 0.0432.

Furthermore, we generated the PRS model to identify the group of individuals with a high risk of severe COVID-19 (Figure 3). In this group, the risk of experiencing a severe course of the disease is approximately twice that of the median risk [with OR = 2.25 (1.7, 2.95), P-value = 3.1 × 10⁻⁹]. The area under the curve (AUC) for the developed PRS is equal to 0.6 (0.58–0.62). The PRS was constructed using the grid model in LDPred2 and included 955,503 SNPs.

Figure 3

Figure 3. Quantile plot of the severe course of COVID-19 PRS developed by LDPred2. The odds ratio represents comparison of PRS odds from different quantiles with the reference quantile (40%−60%). The bars represent the standard deviation (SD).

Discussion

In this study, we identified several genome-wide significant loci associated with the severe forms of COVID-19. The loci include LZTFL1 and PRKCH-AS1/PRKCH genes, intergenic variants on chromosome 3, and a novel loci with intergenic variants on chromosome 6 with genome-wide significance (P-value = 5.92 × 10⁻⁸ to 6.15 × 10⁻⁸). In addition, several genetic variants in the LZTFL1 gene showed genome-wide significance.

The intron variant rs1989566, located within the PRKCH-AS1/PRKCH gene, was newly associated with the susceptibility and severity of COVID-19. The PKC family serves as a mediator for diverse signaling pathways and governs numerous crucial cellular functions, including proliferation, differentiation, and apoptosis. One of the PKC family members, protein kinase C eta protein (PKCη), which is encoded by the PRKCH gene, is a serine-threonine kinase. It is predominantly expressed in vascular endothelial cells and plays a role in the progression and exacerbation of atherosclerosis and subsequently stroke, as indicated by previous studies (29, 30).

Furthermore, the identified variant may potentially impact the expression of the PRKCH-AS1/PRKCH gene and subsequently disrupt endothelial function, which can cause an imbalance in hemostasis favoring a procoagulant state. This is characterized by impaired vasodilator release, increased release of vasoconstrictors, heightened microvasculature spastic reactions, increased leukocyte migration across the endothelium, and the initiation of localized inflammation. Prolonged exposure to factors that induce endothelial dysfunction can contribute to a pro-inflammatory and prothrombotic phenotype in endothelial cells (31). This exposure also results in a depletion of the pool of progenitor endothelial cells (32), ultimately limiting the capacity for restoring their normal phenotype and function.

In addition to regulating cell proliferation, differentiation, and cell death, PKCη is also expressed in the lung tissue, immune system, and proliferation pathways, for instance, in activating nuclear factor κB (NF-κB) signaling, which leads to anti-cancer drug resistance (33–35).

We confirmed the effect of rs10490770 on severe course of COVID-19 (OR = 1.57, P-value = 3.51 × 10⁻⁸) in the LZTFL1 gene. Previous studies (11, 36) have shown that rs10490770 increases the risks of all-cause mortality (HR = 1.4), severe respiratory failure (OR, 2.1), venous thromboembolism (OR = 1.7), and hepatic injury (OR = 1.5). For elderly people, it significantly increases the risk of mortality or severe respiratory failure by more than a double (OR, 2.7).

The LZTFL1 gene, expressed in the normal lung epithelium, is involved in protein transport to the cilia of the ciliated epithelium respiratory cells (37). The transcriptome analysis of lung biopsies from patients with COVID-19 showed the presence of signals associated with epithelial-mesenchymal transition of lung cells (EMT) or pulmonary fibrosis, which is regulated by LZTFL1, suggesting that this locus may serve as a potential therapeutic target (38). Other studies have also linked the LZTFL1 gene at the 3p21.31 locus to COVID-19 infection (39, 40).

Several limitations could affect our GWAS results, including relatively small number of controls, low trait heritability (h2 coefficient), and high polygenicity. Previous research has established that age is a significant risk factor for severe COVID-19, with older individuals being at a higher risk. In our study, the patient cases were sourced from hospitals with an average age of 64.77 years, while controls were selected from the database of a genetic testing company, Genotek, with an average age of ~35 years. We did not use a population control, as is commonly done in many studies. Instead, we defined controls as individuals aged 40 years or above who reported having COVID-19 but did not experience a severe course of the disease. The age threshold was set at 40 years to balance the age distribution and the sizes of our case and control cohorts.

Our primary aim was to evaluate the predictive ability of the PRS of COVID-19 severity in the Eastern European cohort. Therefore, we performed the GWAS analysis on a cohort of 7,124 individuals and constructed a polygenic risk model. Stratification of individuals based on PRS quantiles revealed that the high-risk category (top 10%) had twice the risk of severe course of COVID-19 compared to the median risk group. These results support findings from previous European cohort studies demonstrating similar associations of the PRS with severe COVID-19 (13, 18–20). Before implementation, the PRS must be validated in independent cohorts. Additional prospective studies can be beneficial to assess the clinical utility of the PRS in a practical setting. The potential application of the PRS could involve stratifying all tested individuals into high- and low-risk groups. In our study, we demonstrated the predictive power of the developed PRS. The selection of a threshold for the PRS that delineates the high-risk group should be determined based on factors such as the anticipated increase in hospitalization risk, mortality rates, and economic considerations.

Data availability statement

For the Genotek dataset, the user agreement (available at https://www.genotek.ru) states that disclosure of individual-level genetic information and/or self-reported information to third parties for research purposes will not occur without explicit consent, and the consent was not obtained from the individuals. Due to the user agreement, the individual level cannot be made directly available, and the dataset could pose a threat to confidentiality. Data have to be accessed indirectly via Genotek Ltd, https://www.genotek.ru/. Data requests should be sent to Genotek Ltd at aW5mbyYjeDAwMDQwO2dlbm90ZWsucnU=. The summary statistics of the severe course of COVID-19 GWAS based on the 7,124 Eastern-European participants is freely available at the NHGRI-EBI GWAS Catalog (https://www.ebi.ac.uk/gwas/home) with the accession number GCST90444397. PRS weights are available via the PGS catalog (https://www.pgscatalog.org), with the publication ID: PGP000666 and score IDs: PGS004938.

Ethics statement

The studies involving humans were approved by the Ethics Committee of Genotek Ltd. (protocol No15 “GWAS of severe COVID-19 in the Russian population”). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

EKo: Conceptualization, Data curation, Validation, Writing – review & editing. LS: Conceptualization, Formal analysis, Writing – original draft. EV: Writing – original draft, Writing – review & editing. AKa: Conceptualization, Writing – original draft. VR: Formal analysis, Software, Validation, Writing – original draft. DK: Methodology, Software, Visualization, Writing – original draft. AKi: Conceptualization, Methodology, Project administration, Writing – original draft. NP: Conceptualization, Project administration, Supervision, Writing – review & editing. AE: Project administration, Supervision, Validation, Writing – review & editing. NB: Data curation, Methodology, Writing – original draft, Writing – review & editing. MS: Data curation, Writing – original draft. SS: Data curation, Methodology, Validation, Writing – original draft. IG: Data curation, Methodology, Writing – original draft. PM: Methodology, Project administration, Writing – original draft. EKh: Project administration, Writing – original draft, Writing – review & editing. SP: Methodology, Resources, Supervision, Writing – original draft. AI: Project administration, Writing – review & editing. VI: Writing – review & editing, Project administration. AR: Project administration, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. Genotyping of severe COVID-19 samples (cases) was sponsored by Yandex LLC. The funder Yandex LLC was not involved in the study design, collection, analysis, interpretation of data, the writing of this article, or the decision to submit it for publication.

Acknowledgments

We would like to extend our sincere thanks to Yandex LLC for financial support to perform the current study.

Conflict of interest

EKo, LS, EV, AKa, VR, DK, AKi, NP, and AR were employed by Genotek Ltd. AE was employed by Genetic Technologies Ltd. AI and VI were employed by Eligens SIA.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2024.1409714/full#supplementary-material

References

1. WHO coronavirus (COVID-19) dashboard. WHO Coronavirus (COVID-19) Dashboard with Vaccination Data. Available at: https://covid19.who.int/ (accessed November 30, 2023).

Google Scholar

2. Booth A, Reed AB, Ponzo S, Yassaee A, Aral M, Plans D, et al. Population risk factors for severe disease and mortality in COVID-19: a global systematic review and meta-analysis. PLoS ONE. (2021) 16:e0247461. doi: 10.1371/journal.pone.0247461

PubMed Abstract | Crossref Full Text | Google Scholar

3. Xie Y, Wang Z, Liao H, Marley G, Wu D, Tang W, et al. Epidemiologic, clinical, and laboratory findings of the COVID-19 in the current pandemic: systematic review and meta-analysis. BMC Infect Dis. (2020) 20:640. doi: 10.1186/s12879-020-05371-2

PubMed Abstract | Crossref Full Text | Google Scholar

4. Zheng Z, Peng F, Xu B, Zhao J, Liu H, Peng J, et al. Risk factors of critical and mortal COVID-19 cases: a systematic literature review and meta-analysis. J Infect. (2020) 81:e16–25. doi: 10.1016/j.jinf.2020.04.021

PubMed Abstract | Crossref Full Text | Google Scholar

5. Casanova JL, Su HC, COVID Human Genetic Effort. A global effort to define the human genetics of protective immunity to SARS-CoV-2 infection. Cell. (2020) 181:1194–9. doi: 10.1016/j.cell.2020.05.016

PubMed Abstract | Crossref Full Text | Google Scholar

6. Russell CD, Lone NI, Baillie JK. Comorbidities, multimorbidity and COVID-19. Nat Med. (2023) 29:334–43. doi: 10.1038/s41591-022-02156-9

PubMed Abstract | Crossref Full Text | Google Scholar

7. Chatterjee S, Nalla LV, Sharma M, Sharma N, Singh AA, Malim FM, et al. Association of COVID-19 with comorbidities: an update. ACS Pharmacol Transl Sci. (2023) 6:334–54. doi: 10.1021/acsptsci.2c00181

PubMed Abstract | Crossref Full Text | Google Scholar

8. Ng WH. Comorbidities in SARS-CoV-2 patients: a systematic review and meta-analysis. MBio. (2021). doi: 10.1128/mBio.03647-20

PubMed Abstract | Crossref Full Text | Google Scholar

9. Petrilli CM, Jones SA, Yang J, Rajagopalan H, O'Donnell L, Chernyak Y, et al. Factors associated with hospital admission and critical illness among 5279 people with coronavirus disease 2019 in New York City: prospective cohort study. BMJ. (2020) 369:m1966. doi: 10.1136/bmj.m1966

PubMed Abstract | Crossref Full Text | Google Scholar

10. Grasselli G, Zangrillo A, Zanella A, Antonelli M, Cabrini L, Castelli A, et al. Baseline characteristics and outcomes of 1591 patients infected with SARS-CoV-2 admitted to ICUs of the Lombardy Region, Italy. JAMA. (2020) 323:1574–81. doi: 10.1001/jama.2020.5394

PubMed Abstract | Crossref Full Text | Google Scholar

11. COVID-19 Host Genetics Initiative. Mapping the human genetic architecture of COVID-19. Nature. (2021) 600:472–7. doi: 10.1038/s41586-021-03767-x

PubMed Abstract | Crossref Full Text | Google Scholar

12. COVID-19 Host Genetics Initiative. A first update on mapping the human genetic architecture of COVID-19. Nature. (2022) 608:E1–10. doi: 10.1038/s41586-022-04826-7

PubMed Abstract | Crossref Full Text | Google Scholar

13. 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–92. doi: 10.1038/s41588-021-01006-7

PubMed Abstract | Crossref Full Text | Google Scholar

14. Ishak A, Mehendale M, AlRawashdeh MM, Sestacovschi C, Sharath M, Pandav K, et al. The association of COVID-19 severity and susceptibility and genetic risk factors: a systematic review of the literature. Gene. (2022) 836:146674. doi: 10.1016/j.gene.2022.146674

PubMed Abstract | Crossref Full Text | Google Scholar

15. Li Y, Ke Y, Xia X, Wang Y, Cheng F, Liu X, et al. Genome-wide association study of COVID-19 severity among the Chinese population. Cell Discov. (2021) 7:1–16. doi: 10.1038/s41421-021-00318-6

PubMed Abstract | Crossref Full Text | Google Scholar

16. 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. doi: 10.1016/j.ebiom.2021.103629

PubMed Abstract | Crossref Full Text | Google Scholar

17. Weiner J, Suwalski P, Holtgrewe M, Rakitko A, Thibeault C, Müller M, et al. Increased risk of severe clinical course of COVID-19 in carriers of HLA-C^*04:01. EClinicalMedicine. (2021) 40:101099. doi: 10.1016/j.eclinm.2021.101099

PubMed Abstract | Crossref Full Text | Google Scholar

18. 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:1–15. doi: 10.1186/s12920-023-01584-x

PubMed Abstract | Crossref Full Text | Google Scholar

19. Crossfield SSR, Chaddock NJM, Iles MM, Pujades-Rodriguez M, Morgan AW. Interplay between demographic, clinical and polygenic risk factors for severe COVID-19. Int J Epidemiol. (2022) 51:1384–95. doi: 10.1093/ije/dyac137

PubMed Abstract | Crossref Full Text | Google Scholar

20. Nostaeva AV, Shimansky VS, Apalko SV, Kuznetsov IA, Sushentseva NN, Popov OS, et al. Analysis of associations between polygenic risk score and COVID-19 severity in russian population using low-pass genome sequencing. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory. doi: 10.1101/2023.11.20.23298335

Crossref Full Text | Google Scholar

21. Katki HA, Berndt SI, Machiela MJ, Stewart DR, Garcia-Closas M, Kim J, et al. Increase in power by obtaining 10 or more controls per case when type-1 error is small in large-scale association studies. BMC Med Res Methodol. (2023) 23:153. doi: 10.1186/s12874-023-01973-x

PubMed Abstract | Crossref Full Text | Google Scholar

22. Browning SR, Browning BL. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet. (2007) 81:1084–97. doi: 10.1086/521987

PubMed Abstract | Crossref Full Text | Google Scholar

23. 1000 Genomes Project Consortium, Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM. A global reference for human genetic variation. Nature. (2015) 526:68–74. doi: 10.1038/nature15393

PubMed Abstract | Crossref Full Text | Google Scholar

24. McCarthy S, Das S, Kretzschmar W, Delaneau O, Wood AR, Teumer A, et al. A reference panel of 64,976 haplotypes for genotype imputation. Nat Genet. (2016) 48:1279–83. doi: 10.1038/ng.3643

PubMed Abstract | Crossref Full Text | Google Scholar

25. Staples J Qiao D Cho MH Silverman EK University University of Washington Center for Mendelian Genomics Nickerson DA . PRIMUS: rapid reconstruction of pedigrees from genome-wide estimates of identity by descent. Ame J Hum Genet. (2014) 95:553–64. doi: 10.1016/j.ajhg.2014.10.005

PubMed Abstract | Crossref Full Text | Google Scholar

26. Hahsler M, Piekenbrock M, Doran D. dbscan: fast density-based clustering with R. J Stat Softw. (2019) 91:1–30. doi: 10.18637/jss.v091.i01

Crossref Full Text | Google Scholar

27. COVID-19 Host Genetics Initiative. Available at: https://www.covid19hg.org/ (accessed November 30, 2023).

Google Scholar

28. Privé F, Arbel J. LDpred2: better, faster, stronger. Bioinformatics. (2020) 36:5424–31. doi: 10.1093/bioinformatics/btaa1029

PubMed Abstract | Crossref Full Text | Google Scholar

29. Kubo M, Hata J, Ninomiya T, Matsuda K, Yonemoto K, Nakano T, et al. A nonsynonymous SNP in PRKCH (protein kinase C η) increases the risk of cerebral infarction. Nat Genet. (2007) 39:212–7. doi: 10.1038/ng1945

PubMed Abstract | Crossref Full Text | Google Scholar

30. Zhang Z, Xu G, Zhu W, Cao L, Yan B, Liu X. PRKCH 1425G/A polymorphism predicts recurrence of ischemic stroke in a Chinese population. Mol Neurobiol. (2014) 52:1648–53. doi: 10.1007/s12035-014-8964-6

PubMed Abstract | Crossref Full Text | Google Scholar

31. Chiva-Blanch G, Sala-Vila A, Crespo J, Ros E, Estruch R, Badimon L, et al. The Mediterranean diet decreases prothrombotic microvesicle release in asymptomatic individuals at high cardiovascular risk. Clin Nutr. (2020) 39:3377–84. doi: 10.1016/j.clnu.2020.02.027

PubMed Abstract | Crossref Full Text | Google Scholar

32. Schmidt-Lucke C, Rössig L, Fichtlscherer S, Vasa M, Britten M, Kämper U, et al. Reduced number of circulating endothelial progenitor cells predicts future cardiovascular events. Circulation. (2005) 111:2981–7. doi: 10.1161/CIRCULATIONAHA.104.504340

PubMed Abstract | Crossref Full Text | Google Scholar

33. Fu G, Gascoigne NRJ. The role of protein kinase Cη in T cell biology. Front Immunol. (2012) 3. doi: 10.3389/fimmu.2012.00177

PubMed Abstract | Crossref Full Text | Google Scholar

34. Basu A. The enigmatic protein kinase C-eta. Cancers. (2019) 11. doi: 10.3390/cancers11020214

PubMed Abstract | Crossref Full Text | Google Scholar

35. Raveh-Amit H, Hai N, Rotem-Dai N, Shahaf G, Gopas J, Livneh E, et al. Protein kinase Cη activates NF-κB in response to camptothecin-induced DNA damage. Biochem Biophys Res Commun. (2011) 412:313–7. doi: 10.1016/j.bbrc.2011.07.090

PubMed Abstract | Crossref Full Text | Google Scholar

36. Nakanishi T, Pigazzini S, Degenhardt F, Cordioli M, Butler-Laporte G, Maya-Miles D, et al. Age-dependent impact of the major common genetic risk factor for COVID-19 on severity and mortality. J Clin Invest. (2021) 131. doi: 10.1172/JCI152386

PubMed Abstract | Crossref Full Text | Google Scholar

37. Wei Q, Chen ZH, Wang L, Zhang T, Duan L, Behrens C, et al. LZTFL1 suppresses lung tumorigenesis by maintaining differentiation of lung epithelial cells. Oncogene. (2015) 35:2655–63. doi: 10.1038/onc.2015.328

PubMed Abstract | Crossref Full Text | Google Scholar

38. Downes DJ, Cross AR, Hua P, Roberts N, Schwessinger R, Cutler AJ, et al. Identification of LZTFL1 as a candidate effector gene at a COVID-19 risk locus. Nat Genet. (2021) 53:1606–15. doi: 10.1038/s41588-021-00955-3

PubMed Abstract | Crossref Full Text | Google Scholar

39. Ma Y, Huang Y, Zhao S, Yao Y, Zhang Y, Qu J, et al. Integrative genomics analysis reveals a 21q22.11 locus contributing risk to COVID-19. Hum Mol Genet. (2021) 30:1247–58. doi: 10.1093/hmg/ddab125

PubMed Abstract | Crossref Full Text | Google Scholar

40. Ma Y, Qiu F, Deng C, Li J, Huang Y, Wu Z, et al. Integrating single-cell sequencing data with GWAS summary statistics reveals CD16+ monocytes and memory CD8+ T cells involved in severe COVID-19. Genome Med. (2022) 141:16. doi: 10.1101/2022.02.06.21266924

PubMed Abstract | Crossref Full Text | Google Scholar