Automated workflow-based exploitation of pathway databases provides new insights into genetic associations of metabolite profiles

Background: Genome-wide association studies (GWAS) have identified many common single nucleotide polymorphisms (SNPs) that associate with clinical phenotypes, but these SNPs usually explain just a small part of the heritability and have relatively modest effect sizes. In contrast, SNPs that associate with metabolite levels generally explain a higher percentage of the genetic variation and demonstrate larger effect sizes. Still, the discovery of SNPs associated with metabolite levels is challenging since testing all metabolites measured in typical metabolomics studies with all SNPs comes with a severe multiple testing penalty. We have developed an automated workflow approach that utilizes prior knowledge of biochemical pathways present in databases like KEGG and BioCyc to generate a smaller SNP set relevant to the metabolite. This paper explores the opportunities and challenges in the analysis of GWAS of metabolomic phenotypes and provides novel insights into the genetic basis of metabolic variation through the re-analysis of published GWAS datasets. Results: Re-analysis of the published GWAS dataset from Illig et al. (Nature Genetics, 2010) using a pathway-based workflow (http://www.myexperiment.org/packs/319.html), confirmed previously identified hits and identified a new locus of human metabolic individuality, associating Aldehyde dehydrogenase family1 L1 (ALDH1L1) with serine/glycine ratios in blood. Replication in an independent GWAS dataset of phospholipids (Demirkan et al., PLoS Genetics, 2012) identified two novel loci supported by additional literature evidence: GPAM (Glycerol-3 phosphate acyltransferase) and CBS (Cystathionine beta-synthase). In addition, the workflow approach provided novel insight into the affected pathways and relevance of some of these gene-metabolite pairs in disease development and progression. Conclusions: We demonstrate the utility of automated exploitation of background knowledge present in pathway databases for the analysis of GWAS datasets of metabolomic phenotypes. We report novel loci and potential biochemical mechanisms that contribute to our understanding of the genetic basis of metabolic variation and its relationship to disease development and progression

Genomewide meta-analysis identifies loci associated with IGF-I and IGFBP-3 levels with impact on age-related traits

The growth hormone/insulin-like growth factor (IGF) axis can be manipulated in animal models to promote longevity, and IGF-related proteins including IGF-I and IGF-binding protein-3 (IGFBP-3) have also been implicated in risk of human diseases including cardiovascular diseases, diabetes, and cancer. Throug

Genome-wide analyses identify a role for SLC17A4 and AADAT in thyroid hormone regulation.

Thyroid dysfunction is an important public health problem, which affects 10% of the general population and increases the risk of cardiovascular morbidity and mortality. Many aspects of thyroid hormone regulation have only partly been elucidated, including its transport, metabolism, and genetic determinants. Here we report a large meta-analysis of genome-wide association studies for thyroid function and dysfunction, testing 8 million genetic variants in up to 72,167 individuals. One-hundred-and-nine independent genetic variants are associated with these traits. A genetic risk score, calculated to assess their combined effects on clinical end points, shows significant associations with increased risk of both overt (Graves' disease) and subclinical thyroid disease, as well as clinical complications. By functional follow-up on selected signals, we identify a novel thyroid hormone transporter (SLC17A4) and a metabolizing enzyme (AADAT). Together, these results provide new knowledge about thyroid hormone physiology and disease, opening new possibilities for therapeutic targets

Dipeptidyl peptidase-1 inhibition in patients hospitalised with COVID-19: a multicentre, double-blind, randomised, parallel-group, placebo-controlled trial

Background Neutrophil serine proteases are involved in the pathogenesis of COVID-19 and increased serine protease activity has been reported in severe and fatal infection. We investigated whether brensocatib, an inhibitor of dipeptidyl peptidase-1 (DPP-1; an enzyme responsible for the activation of neutrophil serine proteases), would improve outcomes in patients hospitalised with COVID-19. Methods In a multicentre, double-blind, randomised, parallel-group, placebo-controlled trial, across 14 hospitals in the UK, patients aged 16 years and older who were hospitalised with COVID-19 and had at least one risk factor for severe disease were randomly assigned 1:1, within 96 h of hospital admission, to once-daily brensocatib 25 mg or placebo orally for 28 days. Patients were randomly assigned via a central web-based randomisation system (TruST). Randomisation was stratified by site and age (65 years or ≥65 years), and within each stratum, blocks were of random sizes of two, four, or six patients. Participants in both groups continued to receive other therapies required to manage their condition. Participants, study staff, and investigators were masked to the study assignment. The primary outcome was the 7-point WHO ordinal scale for clinical status at day 29 after random assignment. The intention-to-treat population included all patients who were randomly assigned and met the enrolment criteria. The safety population included all participants who received at least one dose of study medication. This study was registered with the ISRCTN registry, ISRCTN30564012. Findings Between June 5, 2020, and Jan 25, 2021, 406 patients were randomly assigned to brensocatib or placebo; 192 (47·3%) to the brensocatib group and 214 (52·7%) to the placebo group. Two participants were excluded after being randomly assigned in the brensocatib group (214 patients included in the placebo group and 190 included in the brensocatib group in the intention-to-treat population). Primary outcome data was unavailable for six patients (three in the brensocatib group and three in the placebo group). Patients in the brensocatib group had worse clinical status at day 29 after being randomly assigned than those in the placebo group (adjusted odds ratio 0·72 [95% CI 0·57–0·92]). Prespecified subgroup analyses of the primary outcome supported the primary results. 185 participants reported at least one adverse event; 99 (46%) in the placebo group and 86 (45%) in the brensocatib group. The most common adverse events were gastrointestinal disorders and infections. One death in the placebo group was judged as possibly related to study drug. Interpretation Brensocatib treatment did not improve clinical status at day 29 in patients hospitalised with COVID-19

Expanding the clinical phenotype of individuals with a 3-bp in-frame deletion of the NF1 gene (c.2970_2972del): an update of genotype–phenotype correlation

Purpose: Neurofibromatosis type 1 (NF1) is characterized by a highly variable clinical presentation, but almost all NF1-affected adults present with cutaneous and/or subcutaneous neurofibromas. Exceptions are individuals heterozygous for the NF1 in-frame deletion, c.2970_2972del (p.Met992del), associated with a mild phenotype without any externally visible tumors. Methods: A total of 135 individuals from 103 unrelated families, all carrying the constitutional NF1 p.Met992del pathogenic variant and clinically assessed using the same standardized phenotypic checklist form, were included in this study. Results: None of the individuals had externally visible plexiform or histopathologically confirmed cutaneous or subcutaneous neurofibromas. We did not identify any complications, such as symptomatic optic pathway gliomas (OPGs) or symptomatic spinal neurofibromas; however, 4.8% of individuals had nonoptic brain tumors, mostly low-grade and asymptomatic, and 38.8% had cognitive impairment/learning disabilities. In an individual with the NF1 constitutional c.2970_2972del and three astrocytomas, we provided proof that all were NF1-associated tumors given loss of heterozygosity at three intragenic NF1 microsatellite markers and c.2970_297

Pleiotropy among common genetic loci identified for cardiometabolic disorders and C-reactive protein.

Pleiotropic genetic variants have independent effects on different phenotypes. C-reactive protein (CRP) is associated with several cardiometabolic phenotypes. Shared genetic backgrounds may partially underlie these associations. We conducted a genome-wide analysis to identify the shared genetic background of inflammation and cardiometabolic phenotypes using published genome-wide association studies (GWAS). We also evaluated whether the pleiotropic effects of such loci were biological or mediated in nature. First, we examined whether 283 common variants identified for 10 cardiometabolic phenotypes in GWAS are associated with CRP level. Second, we tested whether 18 variants identified for serum CRP are associated with 10 cardiometabolic phenotypes. We used a Bonferroni corrected p-value of 1.1×10-04 (0.05/463) as a threshold of significance. We evaluated the independent pleiotropic effect on both phenotypes using individual level data from the Women Genome Health Study. Evaluating the genetic overlap between inflammation and cardiometabolic phenotypes, we found 13 pleiotropic regions. Additional analyses showed that 6 regions (APOC1, HNF1A, IL6R, PPP1R3B, HNF4A and IL1F10) appeared to have a pleiotropic effect on CRP independent of the effects on the cardiometabolic phenotypes. These included loci where individuals carrying the risk allele for CRP encounter higher lipid levels and risk of type 2 diabetes. In addition, 5 regions (GCKR, PABPC4, BCL7B, FTO and TMEM18) had an effect on CRP largely mediated through the cardiometabolic phenotypes. In conclusion, our results show genetic pleiotropy among inflammation and cardiometabolic phenotypes. In addition to reverse causation, our data suggests that pleiotropic genetic variants partially underlie the association between CRP and cardiometabolic phenotypes

Genome-wide association meta-analyses and fine-mapping elucidate pathways influencing albuminuria

Increased levels of the urinary albumin-to-creatinine ratio (UACR) are associated with higher risk of kidney disease progression and cardiovascular events, but underlying mechanisms are incompletely understood. Here, we conduct trans-ethnic (n = 564,257) and European-ancestry specific meta-analyses of genome-wide association studies of UACR, including ancestry- and diabetes-specific analyses, and identify 68 UACR-associated loci. Genetic correlation analyses and risk score associations in an independent electronic medical records database (n = 192,868) reveal connections with proteinuria, hyperlipidemia, gout, and hypertension. Fine-mapping and trans-Omics analyses with gene expression in 47 tissues and plasma protein levels implicate genes potentially operating through differential expression in kidney (including TGFB1, MUC1, PRKCI, and OAF), and allow coupling of UACR associations to altered plasma OAF concentrations. Knockdown of OAF and PRKCI orthologs in Drosophila nephrocytes reduces albumin endocytosis. Silencing fly PRKCI further impairs slit diaphragm formation. These results generate a priority list of genes and pathways for translational research to reduce albuminuria

Novel loci affecting iron homeostasis and their effects in individuals at risk for hemochromatosis

Variation in body iron is associated with or causes diseases, including anaemia and iron overload. Here, we analyse genetic association data on biochemical markers of iron status from 11 European-population studies, with replication in eight additional cohorts (total up to 48,972 subjects). We find 11 genome-wide-significant (P<5 × 10-8) loci, some including known iron-related genes (HFE, SLC40A1, TF, TFR2, TFRC, TMPRSS6) and others novel (ABO, ARNTL, FADS2, NAT2, TEX14). SNPs at ARNTL, TF, and TFR2 affect iron markers in HFE C282Y homozygotes at risk for hemochromatosis. There is substantial overlap between our iron loci and loci affecting erythrocyte and lipid phenotypes. These results will facilitate investigation of the roles of iron in disease

Genetic variants in RBFOX3 are associated with sleep latency

Time to fall asleep (sleep latency) is a major determinant of sleep quality. Chronic, long sleep latency is a major characteristic of sleep-onset insomnia and/or delayed sleep phase syndrome. In this study we aimed to discover common polymorphisms that contribute to the genetics of sleep latency. We performed a meta-analysis of genome-wide association studies (GWAS) including 2 572 737 single nucleotide polymorphisms (SNPs) established in seven European cohorts including 4242 individuals. We found a cluster of three highly correlated variants (rs9900428, rs9907432 and rs7211029) in the RNA-binding protein fox-1 homolog 3 gene (RBFOX3) associated with sleep latency (P-values=5.77 × 10-08, 6.59 × 10- 08 and 9.17 × 10- 08). These SNPs were replicated in up to 12 independent populations including 30 377 individuals (P-values=1.5 × 10- 02, 7.0 × 10- 03 and 2.5 × 10- 03; combined meta-analysis P-values=5.5 × 10-07, 5.4 × 10-07 and 1.0 × 10-07). A functional prediction of RBFOX3 base

Genome-wide association meta-analyses and fine-mapping elucidate pathways influencing albuminuria

Increased levels of the urinary albumin-to-creatinine ratio (UACR) are associated with higher risk of kidney disease progression and cardiovascular events, but underlying mechanisms are incompletely understood. Here, we conduct trans-ethnic (n = 564,257) and European-ancestry specific meta-analyses of genome-wide association studies of UACR, including ancestry- and diabetes-specific analyses, and identify 68 UACR-associated loci. Genetic correlation analyses and risk score associations in an independent electronic medical records database (n = 192,868) reveal connections with proteinuria, hyperlipidemia, gout, and hypertension. Fine-mapping and trans-Omics analyses with gene expression in 47 tissues and plasma protein levels implicate genes potentially operating through differential expression in kidney (including TGFB1, MUC1, PRKCI, and OAF), and allow coupling of UACR associations to altered plasma OAF concentrations. Knockdown of OAF and PRKCI orthologs in Drosophila nephrocytes reduces albumin endocytosis. Silencing fly PRKCI further impairs slit diaphragm formation. These results generate a priority list of genes and pathways for translational research to reduce albuminuria