Rare variants affecting host defense against pathogens may be involved in COVID-19 severity, but most rare variants are not expected to have a major impact on the course of COVID-19. We hypothesized that the accumulation of weak effects of many rare functional variants throughout the exome may contribute to the overall risk in patients with severe disease. This assumption is consistent with the omnigenic model of the relationship between genetic and phenotypic variation in complex traits, according to which association signals tend to spread across most of the genome through gene regulatory networks from genes outside the major pathways to disease-related genes. We performed whole-exome sequencing and compared the burden of rare variants in 57 patients with severe and 29 patients with mild/moderate COVID-19. At the whole-exome level, we observed an excess of rare, predominantly high-impact (HI) variants in the group with severe COVID-19. Restriction to genes intolerant to HI or damaging missense variants increased enrichment for these classes of variants. Among various sets of genes, an increased signal of rare HI variants was demonstrated predominantly for primary immunodeficiency genes and the entire set of genes associated with immune diseases, as well as for genes associated with respiratory diseases. We advocate taking the ideas of the omnigenic model into account in COVID-19 studies.

COVID-19 patients with acute respiratory distress syndrome (ARDS) have an immune imbalance when systemic inflammation and dysfunction of circulating T and B cells lead to a more severe disease. Using TREC/KREC analysis, we studied the level of mature naive T and B cells in peripheral blood of COVID-19 patients and its relationship with clinical and laboratory data. TREC/KREC analysis was performed by multiplex real-time quantitative PCR on a sample of 36 patients aged 45 years or younger. The reduced TREC/KREC level was observed in ARDS patients compared with non-ARDS patients, and similar results were found for the deceased patients. During days 6 to 20 of hospitalization, a higher neutrophil-to-lymphocyte ratio (NLR) was detected in ARDS patients compared with non-ARDS patients. TREC/KREC negatively correlated with NLR; the highest correlation was recorded for TREC per 100,000 cells with the coefficient of determination R2 = 0.527. Thus, TREC/KREC analysis is a potential prognostic marker for assessing the severity and outcome in COVID-19.

Dysregulation in cytokine production has been linked to the pathogenesis of various immune-mediated traits, in which genetic variability contributes to the etiopathogenesis. GWA studies have identified many genetic variants in or near cytokine genes, nonetheless, the translation of these findings into knowledge of functional determinants of complex traits remains a fundamental challenge. In this study we aimed at collection, analysis and interpretation of data on cytokines focused on their tissue-specific expression, eQTLs and GWAS traits. Using GO annotations, we generated a list of 314 cytokines and analyzed them with the GTEx resource. Cytokines were highly tissue-specific, 82.3% of cytokines had Tau expression metrics ≥ 0.8. In total, 3077 associations for 1760 unique SNPs in or near 244 cytokines were mapped in the NHGRI-EBI GWAS Catalog. According to the Experimental Factor Ontology resource, the largest numbers of disease associations were related to ‘Inflammatory disease’, ‘Immune system disease’ and ‘Asthma’. The GTEx-based analysis revealed that among GWAS SNPs, 1142 SNPs had eQTL effects and influenced expression levels of 999 eGenes, among them 178 cytokines. Several types of enrichment analysis showed that it was cytokines expression variability that fundamentally contributed to the molecular origins of considered immune-mediated conditions.