Objectives In March 2020, the World Health Organization declared that an infectious respiratory disease caused by a new severe acute respiratory syndrome coronavirus 2 [SARS-CoV-2, causing coronavirus disease 2019 (COVID-19)] became a pandemic. In our study, we have analyzed a large publicly available dataset, the Genome Aggregation Database (gnomAD), as well as a cohort of 37 Russian patients with COVID-19 to assess the influence of different classes of genetic variants in the angiotensin-converting enzyme-2 ( ACE2 ) gene on the susceptibility to COVID-19 and the severity of disease outcome. Results We demonstrate that the European populations slightly differ in alternative allele frequencies at the 2,754 variant sites in ACE2 identified in the gnomAD database. We find that the Southern European population has a lower frequency of missense variants and slightly higher frequency of regulatory variants. However, we found no statistical support for the significance of these differences. We also show that the Russian population is similar to other European populations when comparing the frequencies of the ACE2 variants. Evaluation of the effect of various classes of ACE2 variants on COVID-19 outcome in a cohort of Russian patients showed that common missense and regulatory variants do not explain the differences in disease severity. At the same time, we find several rare ACE2 variants (including rs146598386, rs73195521, rs755766792, and others) that are likely to affect the outcome of COVID-19. Our results demonstrate that the spectrum of genetic variants in ACE2 may partially explain the differences in severity of the COVID-19 outcome.

ObjectiveA critical role in coronavirus disease 2019 (COVID-19) pathogenesis is played by immune dysregulation that leads to a generalized uncontrolled multisystem inflammatory response, caused by overproduction of proinflammatory cytokines, known as “a cytokine storm” (CS), strongly associated with a severe course of disease. The aim of this study is to identify prognostic biomarkers for CS development in COVID-19 patients and integrate them into a prognostic score for CS-associated risk applicable to routine clinical practice.Materials and MethodsThe authors performed a review of 458 medical records from COVID-19 patients (241 men and 217 women aged 60.0 ± 10.0) who received treatment in the St. Petersburg State Budgetary Institution of Healthcare City Hospital 40 (City Hospital 40, St. Petersburg), from Apr. 18, 2020 to Nov. 21, 2020. The patients were split in two groups: one group included 100 patients with moderate disease symptoms; the other group included 358 patients with progressive moderately severe, severe, and extremely severe disease. The National Early Warning Score (NEWS) score was used alongside with clinical assessment, chest computed tomographic (CT) scans, electrocardiography (ECG), and lab tests, like ferritin, C-reactive protein (CRP), interleukin (IL)-6, lactate dehydrogenase (LDH), and D-dimer.ResultsThe basic risk factors for cytokine storms in COVID-19 patients are male gender, age over 40 years, positive test result for replicative severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA, absolute lymphocyte count, dynamics in the NEWS score, as well as LDH, D-dimer, ferritin, and IL-6 levels. These clinical and instrumental findings can be also used as laboratory biomarkers for diagnosis and dynamic monitoring of cytokine storms. The suggested prognostic scale (including the NEWS score dynamics; serum IL-6 greater than 23 pg/ml; serum CRP 50 mg/L or greater; absolute lymphocyte count less than 0.72 × 109/L; positive test result for replicative coronavirus (SARS-CoV-2) RNA; age 40 years and over) is a useful tool to identify patients at a high risk for cytokine storm, requiring an early onset of anti-inflammatory therapy.

The COVID-19 pandemic has drawn the attention of many researchers to the interaction between pathogen and host genomes. Over the last two years, numerous studies have been conducted to identify the genetic risk factors that predict COVID-19 severity and outcome. However, such an analysis might be complicated in cohorts of limited size and/or in case of limited breadth of genome coverage. In this work, we tried to circumvent these challenges by searching for candidate genes and genetic variants associated with a variety of quantitative and binary traits in a cohort of 840 COVID-19 patients from Russia. While we found no gene- or pathway-level associations with the disease severity and outcome, we discovered eleven independent candidate loci associated with quantitative traits in COVID-19 patients. Out of these, the most significant associations correspond to rs1651553 in MYH14p = 1.4 × 10−7), rs11243705 in SETX (p = 8.2 × 10−6), and rs16885 in ATXN1 (p = 1.3 × 10−5). One of the identified variants, rs33985936 in SCN11A, was successfully replicated in an independent study, and three of the variants were found to be associated with blood-related quantitative traits according to the UK Biobank data (rs33985936 in SCN11A, rs16885 in ATXN1, and rs4747194 in CDH23). Moreover, we show that a risk score based on these variants can predict the severity and outcome of hospitalization in our cohort of patients. Given these findings, we believe that our work may serve as proof-of-concept study demonstrating the utility of quantitative traits and extensive phenotyping for identification of genetic risk factors of severe COVID-19.