ABSTRACT
The impact of common and rare variants in COVID-19 host genetics has been widely studied. In particular, in Fallerini et al. (Human genetics, 2022, 141, 147-173), common and rare variants were used to define an interpretable machine learning model for predicting COVID-19 severity. First, variants were converted into sets of Boolean features, depending on the absence or the presence of variants in each gene. An ensemble of LASSO logistic regression models was used to identify the most informative Boolean features with respect to the genetic bases of severity. After that, the Boolean features, selected by these logistic models, were combined into an Integrated PolyGenic Score (IPGS), which offers a very simple description of the contribution of host genetics in COVID-19 severity.. IPGS leads to an accuracy of 55%-60% on different cohorts, and, after a logistic regression with both IPGS and age as inputs, it leads to an accuracy of 75%. The goal of this paper is to improve the previous results, using not only the most informative Boolean features with respect to the genetic bases of severity but also the information on host organs involved in the disease. In this study, we generalize the IPGS adding a statistical weight for each organ, through the transformation of Boolean features into "Boolean quantum features," inspired by quantum mechanics. The organ coefficients were set via the application of the genetic algorithm PyGAD, and, after that, we defined two new integrated polygenic scores (IPGSph1 and IPGSph2). By applying a logistic regression with both IPGS, (IPGSph2 (or indifferently IPGSph1) and age as inputs, we reached an accuracy of 84%-86%, thus improving the results previously shown in Fallerini et al. (Human genetics, 2022, 141, 147-173) by a factor of 10%.
ABSTRACT
Polygenic scores (PGSs) offer the ability to predict genetic risk for complex diseases across the life course; a key benefit over short-term prediction models. To produce risk estimates relevant to clinical and public health decision-making, it is important to account for varying effects due to age and sex. Here, we develop a novel framework to estimate country-, age-, and sex-specific estimates of cumulative incidence stratified by PGS for 18 high-burden diseases. We integrate PGS associations from seven studies in four countries (N = 1,197,129) with disease incidences from the Global Burden of Disease. PGS has a significant sex-specific effect for asthma, hip osteoarthritis, gout, coronary heart disease and type 2 diabetes (T2D), with all but T2D exhibiting a larger effect in men. PGS has a larger effect in younger individuals for 13 diseases, with effects decreasing linearly with age. We show for breast cancer that, relative to individuals in the bottom 20% of polygenic risk, the top 5% attain an absolute risk for screening eligibility 16.3 years earlier. Our framework increases the generalizability of results from biobank studies and the accuracy of absolute risk estimates by appropriately accounting for age- and sex-specific PGS effects. Our results highlight the potential of PGS as a screening tool which may assist in the early prevention of common diseases.
Subject(s)
Genetic Predisposition to Disease , Multifactorial Inheritance , Humans , Male , Female , Multifactorial Inheritance/genetics , Incidence , Middle Aged , Adult , Aged , Diabetes Mellitus, Type 2/genetics , Diabetes Mellitus, Type 2/epidemiology , Risk Factors , Risk Assessment/methods , Global Burden of Disease , Sex Factors , Age FactorsABSTRACT
The clinical manifestations of SARS-CoV-2 infection vary widely among patients, from asymptomatic to life-threatening. Host genetics is one of the factors that contributes to this variability as previously reported by the COVID-19 Host Genetics Initiative (HGI), which identified sixteen loci associated with COVID-19 severity. Herein, we investigated the genetic determinants of COVID-19 mortality, by performing a case-only genome-wide survival analysis, 60 days after infection, of 3904 COVID-19 patients from the GEN-COVID and other European series (EGAS00001005304 study of the COVID-19 HGI). Using imputed genotype data, we carried out a survival analysis using the Cox model adjusted for age, age2, sex, series, time of infection, and the first ten principal components. We observed a genome-wide significant (P-value < 5.0 × 10-8) association of the rs117011822 variant, on chromosome 11, of rs7208524 on chromosome 17, approaching the genome-wide threshold (P-value = 5.19 × 10-8). A total of 113 variants were associated with survival at P-value < 1.0 × 10-5 and most of them regulated the expression of genes involved in immune response (e.g., CD300 and KLR genes), or in lung repair and function (e.g., FGF19 and CDH13). Overall, our results suggest that germline variants may modulate COVID-19 risk of death, possibly through the regulation of gene expression in immune response and lung function pathways.
Subject(s)
COVID-19 , Humans , Genome-Wide Association Study/methods , Genetic Predisposition to Disease , SARS-CoV-2 , GenotypeABSTRACT
Extreme polymorphism of HLA and killer-cell immunoglobulin-like receptors (KIR) differentiates immune responses across individuals. Additional to T cell receptor interactions, subsets of HLA class I act as ligands for inhibitory and activating KIR, allowing natural killer (NK) cells to detect and kill infected cells. We investigated the impact of HLA and KIR polymorphism on the severity of COVID-19. High resolution HLA class I and II and KIR genotypes were determined from 403 non-hospitalized and 1575 hospitalized SARS-CoV-2 infected patients from Italy collected in 2020. We observed that possession of the activating KIR2DS4*001 allotype is associated with severe disease, requiring hospitalization (OR = 1.48, 95% CI 1.20-1.85, pc = 0.017), and this effect is greater in individuals homozygous for KIR2DS4*001 (OR = 3.74, 95% CI 1.75-9.29, pc = 0.003). We also observed the HLA class II allotype, HLA-DPB1*13:01 protects SARS-CoV-2 infected patients from severe disease (OR = 0.49, 95% CI 0.33-0.74, pc = 0.019). These association analyses were replicated using logistic regression with sex and age as covariates. Autoantibodies against IFN-α associated with COVID-19 severity were detected in 26% of 156 hospitalized patients tested. HLA-C*08:02 was more frequent in patients with IFN-α autoantibodies than those without, and KIR3DL1*01502 was only present in patients lacking IFN-α antibodies. These findings suggest that KIR and HLA polymorphism is integral in determining the clinical outcome following SARS-CoV-2 infection, by influencing the course both of innate and adaptive immunity.
Subject(s)
COVID-19 , HLA-DP beta-Chains , Humans , COVID-19/genetics , SARS-CoV-2/genetics , Alleles , Receptors, KIR/genetics , Genotype , Autoantibodies/geneticsABSTRACT
Flavonoids, ubiquitously distributed in the plant world, are regularly ingested with diets rich in fruit, vegetables, wine, and tea. During digestion, they are partially absorbed in the stomach. The present work aimed to assess the in vitro effects of quercetin and ten structurally related flavonoids on the rat gastric fundus smooth muscle, focussing on ATP-dependent K+ (Kir6.1) channels, which play a central role in the regulation of resting membrane potential, membrane excitability and, consequently, of gastric motility. Whole-cell currents through Kir6.1 channels (IKir6.1) were recorded with the patch-clamp technique and the mechanical activity of gastric fundus smooth muscle strips was studied under isometric conditions. Galangin ≈ tamarixetin > quercetin > kaempferol > isorhamnetin ≈ luteolin ≈ fisetin > (±)-taxifolin inhibited pinacidil-evoked, glibenclamide-sensitive IKir6.1 in a concentration-dependent manner. Morin, rutin, and myricetin were ineffective. The steric hindrance of the molecule and the number and position of hydroxyl groups on the B ring played an important role in the activity of the molecule. Molecular docking simulations revealed a possible binding site for flavonoids in the C-terminal domain of the Kir6.1 channel subunit SUR2B, in a flexible loop formed by residues 251 to 254 of chains C and D. Galangin and tamarixetin, but not rutin relaxed both high K+- and carbachol-induced contraction of fundus strips in a concentration-dependent manner. Furthermore, both flavonoids shifted to the right the concentration-relaxation curves to either pinacidil or L-cysteine constructed in strips pre-contracted by high K+, rutin being ineffective. In conclusion, IKir6.1 inhibition exerted by dietary flavonoids might counterbalance their myorelaxant activity, affect gastric accommodation or, at least, some stages of digestion.
Subject(s)
Gastric Fundus , Vasodilator Agents , Rats , Animals , Pinacidil/pharmacology , Vasodilator Agents/pharmacology , Gastric Fundus/metabolism , Quercetin/pharmacology , Molecular Docking Simulation , Potassium Channels/metabolism , Muscle, Smooth/metabolism , Electrophysiology , Rutin , Diet , Sulfonylurea Receptors/metabolismABSTRACT
Male sex represents one of the major risk factors for severe COVID-19 outcome. However, underlying mechanisms that mediate sex-dependent disease outcome are as yet unknown. Here, we identify the CYP19A1 gene encoding for the testosterone-to-estradiol metabolizing enzyme CYP19A1 (also known as aromatase) as a host factor that contributes to worsened disease outcome in SARS-CoV-2-infected males. We analyzed exome sequencing data obtained from a human COVID-19 cohort (n = 2,866) using a machine-learning approach and identify a CYP19A1-activity-increasing mutation to be associated with the development of severe disease in men but not women. We further analyzed human autopsy-derived lungs (n = 86) and detect increased pulmonary CYP19A1 expression at the time point of death in men compared with women. In the golden hamster model, we show that SARS-CoV-2 infection causes increased CYP19A1 expression in the lung that is associated with dysregulated plasma sex hormone levels and reduced long-term pulmonary function in males but not females. Treatment of SARS-CoV-2-infected hamsters with a clinically approved CYP19A1 inhibitor (letrozole) improves impaired lung function and supports recovery of imbalanced sex hormones specifically in males. Our study identifies CYP19A1 as a contributor to sex-specific SARS-CoV-2 disease outcome in males. Furthermore, inhibition of CYP19A1 by the clinically approved drug letrozole may furnish a new therapeutic strategy for individualized patient management and treatment.
Subject(s)
Aromatase , COVID-19 , Female , Humans , Male , Aromatase/genetics , Letrozole , SARS-CoV-2 , COVID-19/genetics , Estradiol , TestosteroneABSTRACT
Resource competition can be the cause of unintended coupling between co-expressed genetic constructs. Here we report the quantification of the resource load imposed by different mammalian genetic components and identify construct designs with increased performance and reduced resource footprint. We use these to generate improved synthetic circuits and optimise the co-expression of transfected cassettes, shedding light on how this can be useful for bioproduction and biotherapeutic applications. This work provides the scientific community with a framework to consider resource demand when designing mammalian constructs to achieve robust and optimised gene expression.
Subject(s)
Mammals , Animals , Mammals/geneticsABSTRACT
Objective: Bipolar disorder is a heritable chronic mental disorder that causes psychosocial impairment through depressive/manic episodes. Familial transmission of bipolar disorder does not follow simple Mendelian patterns of inheritance. The aim of this study was to describe a large family with 12 members affected by bipolar disorder. Whole-exome sequencing was performed for eight members, three of whom were diagnosed with bipolar disorder, and another reported as "borderline." Methods: Whole-exome sequencing data allowed us to select variants that the affected members had in common, including and excluding the "borderline" individual with moderate anxiety and obsessive-compulsive traits. Results: The results favored designating certain genes as predispositional to bipolar disorder: a heterozygous missense variant in CLN6 resulted in a "borderline" phenotype that, if combined with a heterozygous missense variant in ZNF92, is responsible for the more severe bipolar disorder phenotype. Both rare missense changes are predicted to disrupt protein function. Conclusions: Loss of both alleles in CLN6 causes neuronal ceroid lipofuscinosis, a severe progressive childhood neurological disorder. Our results indicate that heterozygous CLN6 carriers, previously reported as healthy, may be susceptible to bipolar disorder later in life if associated with additional variants in ZNF92.
ABSTRACT
BACKGROUND: Neuroblastoma (NB) is the most common solid extracranial paediatric tumour. Genome-wide association studies have driven the discovery of common risk variants, but no large study has investigated the contribution of rare variants to NB susceptibility. Here, we conducted a whole-exome sequencing (WES) of 664 NB cases and 822 controls and used independent validation datasets to identify genes with rare risk variants and involved pathways. METHODS: WES was performed at 50× depth and variants were jointly called in cases and controls. We developed two models to identify mutations with high clinical impact (P/LP model) and to discover less penetrant risk mutations affecting non-canonical cancer pathways (RPV model). We performed a gene-level collapsing test using Firth's logistic regression in 242 selected cancer predisposition genes (CPGs) and a gene-sets burden analysis of biologically-informed pathways. FINDINGS: Twelve percent of patients carried P/LP variants in CPGs and showed a significant enrichment (P = 2.3 × 10-4) compared to controls (6%). We identified P/LP variants in 45 CPGs enriched in homologous recombination (HR) pathway. The most P/LP enriched genes in NB were BRCA1, ALK and RAD51C. Additionally, we found higher RPV burden in gene-sets of neuron differentiation, neural tube development and synapse assembly, and in gene-sets associated with neurodevelopmental disorders (NDD). INTERPRETATION: The high fraction of NB patients with P/LP variants indicates the need of genetic counselling. Furthermore, inherited rare variants predispose to NB development by affecting mechanisms related to HR and neurodevelopmental processes, and demonstrate that NDD genes are altered in NB at the germline level. FUNDING: Associazione Italiana per la Ricerca sul Cancro, Fondazione Italiana per la Lotta al Neuroblastoma, Associazione Oncologia Pediatrica e Neuroblastoma, Regione Campania, Associazione Giulio Adelfio onlus, and Italian Health Ministry.
Subject(s)
Genetic Predisposition to Disease , Neuroblastoma , Humans , Child , Genome-Wide Association Study , Mutation , Neuroblastoma/genetics , Homologous RecombinationABSTRACT
Fast C-type inactivation confers distinctive functional properties to the hERG potassium channel, and its association to inherited and acquired cardiac arrythmias makes the study of the inactivation mechanism of hERG at the atomic detail of paramount importance. At present, two models have been proposed to describe C-type inactivation in K+-channels. Experimental data and computational work on the bacterial KcsA channel support the hypothesis that C-type inactivation results from a closure of the selectivity filter that sterically impedes ion conduction. Alternatively, recent experimental structures of a mutated Shaker channel revealed a widening of the extracellular portion of the selectivity filter, which might diminish conductance by interfering with the mechanism of ion permeation. Here, we performed molecular dynamics simulations of the wild-type hERG, a non-inactivating mutant (hERG-N629D), and a mutant that inactivates faster than the wild-type channel (hERG-F627Y) to find out which and if any of the two reported C-type inactivation mechanisms applies to hERG. Closure events of the selectivity filter were not observed in any of the simulated trajectories but instead, the extracellular section of the selectivity filter deviated from the canonical conductive structure of potassium channels. The degree of widening of the potassium binding sites at the extracellular entrance of the channel was directly related to the degree of inactivation with hERG-F627Y > wild-type hERG > hERG-N629D. These findings support the hypothesis that C-type inactivation in hERG entails a widening of the extracellular entrance of the channel rather than a closure of the selectivity filter.
Subject(s)
Ether-A-Go-Go Potassium Channels , Molecular Dynamics Simulation , Ether-A-Go-Go Potassium Channels/chemistry , Ether-A-Go-Go Potassium Channels/genetics , Potassium/chemistryABSTRACT
OBJECTIVE: Bipolar disorder is a heritable chronic mental disorder that causes psychosocial impairment through depressive/manic episodes. Familial transmission of bipolar disorder does not follow simple Mendelian patterns of inheritance. The aim of this study was to describe a large family with 12 members affected by bipolar disorder. Whole-exome sequencing was performed for eight members, three of whom were diagnosed with bipolar disorder, and another reported as "borderline." METHODS: Whole-exome sequencing data allowed us to select variants that the affected members had in common, including and excluding the "borderline" individual with moderate anxiety and obsessive-compulsive traits. RESULTS: The results favored designating certain genes as predispositional to bipolar disorder: a heterozygous missense variant in CLN6 resulted in a "borderline" phenotype that, if combined with a heterozygous missense variant in ZNF92, is responsible for the more severe bipolar disorder phenotype. Both rare missense changes are predicted to disrupt protein function. CONCLUSIONS: Loss of both alleles in CLN6 causes neuronal ceroid lipofuscinosis, a severe progressive childhood neurological disorder. Our results indicate that heterozygous CLN6 carriers, previously reported as healthy, may be susceptible to bipolar disorder later in life if associated with additional variants in ZNF92.
Subject(s)
Bipolar Disorder , Neuronal Ceroid-Lipofuscinoses , Humans , Membrane Proteins/genetics , Neuronal Ceroid-Lipofuscinoses/diagnosis , Neuronal Ceroid-Lipofuscinoses/geneticsABSTRACT
Carriers of single pathogenic variants of the CFTR (cystic fibrosis transmembrane conductance regulator) gene have a higher risk of severe COVID-19 and 14-day death. The machine learning post-Mendelian model pinpointed CFTR as a bidirectional modulator of COVID-19 outcomes. Here, we demonstrate that the rare complex allele [G576V;R668C] is associated with a milder disease via a gain-of-function mechanism. Conversely, CFTR ultra-rare alleles with reduced function are associated with disease severity either alone (dominant disorder) or with another hypomorphic allele in the second chromosome (recessive disorder) with a global residual CFTR activity between 50 to 91%. Furthermore, we characterized novel CFTR complex alleles, including [A238V;F508del], [R74W;D1270N;V201M], [I1027T;F508del], [I506V;D1168G], and simple alleles, including R347C, F1052V, Y625N, I328V, K68E, A309D, A252T, G542*, V562I, R1066H, I506V, I807M, which lead to a reduced CFTR function and thus, to more severe COVID-19. In conclusion, CFTR genetic analysis is an important tool in identifying patients at risk of severe COVID-19.
Subject(s)
COVID-19 , Cystic Fibrosis , Humans , Cystic Fibrosis Transmembrane Conductance Regulator/genetics , Alleles , Cystic Fibrosis/pathology , COVID-19/genetics , HeterozygoteABSTRACT
Host genetics is a key determinant of COVID-19 outcomes. Previously, the COVID-19 Host Genetics Initiative genome-wide association study used common variants to identify multiple loci associated with COVID-19 outcomes. However, variants with the largest impact on COVID-19 outcomes are expected to be rare in the population. Hence, studying rare variants may provide additional insights into disease susceptibility and pathogenesis, thereby informing therapeutics development. Here, we combined whole-exome and whole-genome sequencing from 21 cohorts across 12 countries and performed rare variant exome-wide burden analyses for COVID-19 outcomes. In an analysis of 5,085 severe disease cases and 571,737 controls, we observed that carrying a rare deleterious variant in the SARS-CoV-2 sensor toll-like receptor TLR7 (on chromosome X) was associated with a 5.3-fold increase in severe disease (95% CI: 2.75-10.05, p = 5.41x10-7). This association was consistent across sexes. These results further support TLR7 as a genetic determinant of severe disease and suggest that larger studies on rare variants influencing COVID-19 outcomes could provide additional insights.
Subject(s)
COVID-19 , Exome , Humans , Exome/genetics , Genome-Wide Association Study , COVID-19/genetics , Genetic Predisposition to Disease , Toll-Like Receptor 7/genetics , SARS-CoV-2/geneticsABSTRACT
We employed a multifaceted computational strategy to identify the genetic factors contributing to increased risk of severe COVID-19 infection from a Whole Exome Sequencing (WES) dataset of a cohort of 2000 Italian patients. We coupled a stratified k-fold screening, to rank variants more associated with severity, with the training of multiple supervised classifiers, to predict severity based on screened features. Feature importance analysis from tree-based models allowed us to identify 16 variants with the highest support which, together with age and gender covariates, were found to be most predictive of COVID-19 severity. When tested on a follow-up cohort, our ensemble of models predicted severity with high accuracy (ACC = 81.88%; AUCROC = 96%; MCC = 61.55%). Our model recapitulated a vast literature of emerging molecular mechanisms and genetic factors linked to COVID-19 response and extends previous landmark Genome-Wide Association Studies (GWAS). It revealed a network of interplaying genetic signatures converging on established immune system and inflammatory processes linked to viral infection response. It also identified additional processes cross-talking with immune pathways, such as GPCR signaling, which might offer additional opportunities for therapeutic intervention and patient stratification. Publicly available PheWAS datasets revealed that several variants were significantly associated with phenotypic traits such as "Respiratory or thoracic disease", supporting their link with COVID-19 severity outcome.
Subject(s)
COVID-19 , Humans , COVID-19/genetics , Genome-Wide Association Study , Genetic Predisposition to Disease , Exome Sequencing , PhenotypeABSTRACT
Whole exome sequencing has provided significant opportunities to discover novel candidate genes for intellectual disability and autism spectrum disorders. Variants in the spectrin genes SPTAN1, SPTBN1, SPTBN2, and SPTBN4 have been associated with neurological disorders; however, SPTBN5 gene-variants have not been associated with any human disorder. This is the first report that associates SPTBN5 gene variants (ENSG00000137877: c.266A>C; p.His89Pro, c.9784G>A; p.Glu3262Lys, c.933C>G; p.Tyr311Ter, and c.8809A>T; p.Asn2937Tyr) causing neurodevelopmental phenotypes in four different families. The SPTBN5-associated clinical traits in our patients include intellectual disability (mild to severe), aggressive tendencies, accompanied by variable features such as craniofacial and physical dysmorphisms, autistic behavior, and gastroesophageal reflux. We also provide a review of the existing literature related to other spectrin genes, which highlights clinical features partially overlapping with SPTBN5.