Rare Variants in Primary Immunodeficiency Genes and Their Functional Partners in Severe COVID-19.

The development of severe COVID-19, which is a complex multisystem disease, is thought to be associated with many genes whose action is modulated by numerous environmental and genetic factors. In this study, we focused on the ideas of the omnigenic model of heritability of complex traits, which assumes that a small number of core genes and a large pool of peripheral genes expressed in disease-relevant tissues contribute to the genetics of complex traits through interconnected networks. We hypothesized that primary immunodeficiency disease (PID) genes may be considered as core genes in severe COVID-19, and their functional partners (FPs) from protein-protein interaction networks may be considered as peripheral near-core genes. We used whole-exome sequencing data from patients aged ≤ 45 years with severe (n = 9) and non-severe COVID-19 (n = 11), and assessed the cumulative contribution of rare high-impact variants to disease severity. In patients with severe COVID-19, an excess of rare high-impact variants was observed at the whole-exome level, but maximal association signals were detected for PID + FP gene subsets among the genes intolerant to LoF variants, haploinsufficient and essential. Our exploratory study may serve as a model for new directions in the research of host genetics in severe COVID-19.

COVID-19 severity: does the genetic landscape of rare variants matter?

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.

Pneumonia and Related Conditions in Critically Ill Patients-Insights from Basic and Experimental Studies.

Pneumonia is an acute infectious disease with high morbidity and mortality rates. Pneumonia's development, severity and outcome depend on age, comorbidities and the host immune response. In this study, we combined theoretical and experimental investigations to characterize pneumonia and its comorbidities as well as to assess the host immune response measured by TREC/KREC levels in patients with pneumonia. The theoretical study was carried out using the Columbia Open Health Data (COHD) resource, which provides access to clinical concept prevalence and co-occurrence from electronic health records. The experimental study included TREC/KREC assays in young adults (18-40 years) with community-acquired (CAP) (n = 164) or nosocomial (NP) (n = 99) pneumonia and healthy controls (n = 170). Co-occurring rates between pneumonia, sepsis, acute respiratory distress syndrome (ARDS) and some other related conditions common in intensive care units were the top among 4170, 3382 and 963 comorbidities in pneumonia, sepsis and ARDS, respectively. CAP patients had higher TREC levels, while NP patients had lower TREC/KREC levels compared to controls. Low TREC and KREC levels were predictive for the development of NP, ARDS, sepsis and lethal outcome (AUCTREC in the range 0.71-0.82, AUCKREC in the range 0.67-0.74). TREC/KREC analysis can be considered as a potential prognostic test in patients with pneumonia.

Genes and Diseases: Insights from Transcriptomics Studies.

Results of expression studies can be useful to clarify the genotype-phenotype relationship. However, according to data from recent literature, there is a large group of genes that are revealed as differentially expressed (DE) in many studies, regardless of the biological context. Additional analyses could shed more light on the relationships between genes, their differential expression, and diseases. We generated a set of 9972 disease genes from five gene-phenotype databases (OMIM, ORPHANET, DDG2P, DisGeNet and MalaCards) and a report of the International Union of Immunological Societies. To study transcriptomics of disease and non-disease genes in healthy tissues, we obtained data from the Human Protein Atlas (HPA) website. We analyzed the dependency between expression in healthy tissues and gene occurrence in Gene Expression Omnibus series using tools within the Enrichr libraries. The results of expression studies were annotated with Gene Ontology (GO) and Human Phenotype Ontology (HPO) terms. Using transcriptomics analysis of healthy tissues, we validated the previous findings of higher expression levels of disease genes in pathologically linked tissues compared to other tissues. Preferentially DE genes were generally highly expressed in one or multiple tissues and were enriched for disease genes. According to the results of GO enrichment analyses, both down- and up-regulated DE genes most often took part in immune response, translation and tissue-specific processes. A connection between DE-related pathology and the diversity of HPO terms was found. Investigating a link between expression and phenotype contributes to understanding the mode of development and progression of human diseases.

An overview of germline variations in genes of primary immunodeficiences through integrative analysis of ClinVar, HGMD® and dbSNP databases.

Primary immunodeficiencies (PID) are a diverse group of genetic disorders caused by inadequate development and function of immune system. Identifying genetic etiology is important for genetic counselling and treatment decisions. Clinical relevance of genetic variants is a complex problem depending on gene-specific and variant specific genotype-phenotype interactions. To address this challenge, we aimed to characterize the pathogenic landscape of PID genes by combining the analysis of germline variations reported in ClinVar and HGMD® and identification of damaging variations available in dbSNP. We generated a joint ClinVar/HGMD database, which included 111,940 variants, among them 32,452 were classified as pathogenic/likely pathogenic. From a total of 5,415,794 bi- or multiallelic variants in PID genes recorded in dbSNP, we retrieved 38,291 high impact (HI) biallelic variants with presumably disruptive impact in the protein, of them 25,500 variants were not present in ClinVar/HGMD. Using a functional prediction algorithm, we additionally identified 28,507 deleterious and 56,016 neutral missense variants among dbSNP variants and created a collection of damaging and neutral variations in PID genes, not currently present in ClinVar/HGMD, with their allele frequencies and mappings to protein domains. The distribution of pathogenic variants from ClinVar/HGMD, HI variants and deleterious missense variants from dbSNP was analyzed in the context of hereditary pattern and gene specific metrics, such as pLI and haploinsufficiency. Our report summarized data on complex gene-specific variability in PID genes and might be useful for the identification of the most promising variants and gene regions for further study.

Cytokines mapping for tissue-specific expression, eQTLs and GWAS traits.

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.

Expression changes in pelvic organ prolapse: a systematic review and in silico study.

Pelvic organ prolapse (POP) is a highly disabling condition common for a vast number of women worldwide. To contribute to existing knowledge in POP pathogenesis, we performed a systematic review of expression studies on both specific gene and whole-genome/proteome levels and an in silico analysis of publicly available datasets related to POP development. The most extensively investigated genes in individual studies were related to extracellular matrix (ECM) organization. Three premenopausal and two postmenopausal sets from two Gene Expression Omnibus (GEO) studies (GSE53868 and GSE12852) were analyzed; Gene Ontology (GO) terms related to tissue repair (locomotion, biological adhesion, immune processes and other) were enriched in all five datasets. Co-expression was higher in cases than in controls in three premenopausal sets. The shared between two or more datasets up-regulated genes were enriched with those related to inflammatory bowel disease (IBD) in the NHGRI GWAS Catalog. ECM-related genes were not over-represented among differently expressed genes. Up-regulation of genes related to tissue renewal probably reflects compensatory mechanisms aimed at repair of damaged tissue. Inefficiency of this process may have different origins including age-related deregulation of gene expression.

Biological findings from the PheWAS catalog: focus on connective tissue-related disorders (pelvic floor dysfunction, abdominal hernia, varicose veins and hemorrhoids).

Pelvic floor dysfunction, specifically genital prolapse (GP) and stress urinary inconsistency (SUI) presumably co-occur with other connective tissue disorders such as hernia, hemorrhoids, and varicose veins. Observations on non-random coexistence of these disorders have never been summarized in a meta-analysis. The performed meta-analysis demonstrated that varicose veins and hernia are associated with GP. Disease connections on the molecular level may be partially based on shared genetic susceptibility. A unique opportunity to estimate shared genetic susceptibility to disorders is provided by a PheWAS (phenome-wide association study) designed to utilize GWAS data concurrently to many phenotypes. We searched the PheWAS Catalog, which includes the results of the PheWAS study with P value < 0.05, for genes associated with GP, SUI, abdominal hernia, varicose veins and hemorrhoids. We found pronounced signals for the associations of the SLC2A9 gene with SUI (P = 6.0e-05) and the MYH9 gene with varicose veins of lower extremity (P = 0.0001) and hemorrhoids (P = 0.0007). The comparison of the PheWAS Catalog and the NHGRI Catalog data revealed enrichment of genes associated with bone mineral density in GP and with activated partial thromboplastin time in varicose veins of lower extremity. In cross-phenotype associations, genes responsible for peripheral nerve functions seem to predominate. This study not only established novel biologically plausible associations that may warrant further studies but also exemplified an effective use of the PheWAS Catalog data.

Verification of the Chromosome Region 9q21 Association with Pelvic Organ Prolapse Using RegulomeDB Annotations.

Pelvic organ prolapse (POP) is a common highly disabling disorder with a large hereditary component. It is characterized by a loss of pelvic floor support that leads to the herniation of the uterus in or outside the vagina. Genome-wide linkage studies have shown an evidence of POP association with the region 9q21 and six other loci in European pedigrees. The aim of our study was to test the above associations in a case-control study in Russian population. Twelve SNPs including SNPs cited in the above studies and those selected using the RegulomeDB annotations for the region 9q21 were genotyped in 210 patients with POP (stages III-IV) and 292 controls with no even minimal POP. Genotyping was performed using the polymerase chain reaction with confronting two-pair primers (PCR-CTPP). Association analyses were conducted for individual SNPs, 9q21 haplotypes, and SNP-SNP interactions. SNP rs12237222 with the highest RegulomeDB score 1a appeared to be the key SNP in haplotypes associated with POP. Other RegulomeDB Category 1 SNPs, rs12551710 and rs2236479 (scores 1d and 1f, resp.), exhibited epistatic effects. In this study, we verified the region 9q21 association with POP in Russians, using RegulomeDB annotations.