Uromodulin in sepsis and severe pneumonia: a two-sample Mendelian randomization study.

The outcome for patients with sepsis-associated acute kidney injury in the intensive care unit (ICU) remains poor. Low serum uromodulin (sUMOD) protein levels have been proposed as a causal mediator of this effect. We investigated the effect of different levels of sUMOD on the risk of sepsis and severe pneumonia and outcomes in these conditions. A two-sample Mendelian randomization (MR) study was performed. Single-nucleotide polymorphisms (SNPs) associated with increased levels of sUMOD were identified and used as instrumental variables for association with outcomes. Data from different cohorts were combined based on disease severity and meta-analyzed. Five SNPs associated with increased sUMOD levels were identified and tested in six datasets from two biobanks. There was no protective effect of increased levels of sUMOD on the risk of sepsis [two cohorts, odds ratio (OR) 0.99 (95% confidence interval 0.95-1.03), P = 0.698, and OR 0.95 (0.91-1.00), P = 0.060, respectively], risk of sepsis requiring ICU admission [OR 1.04 (0.93-1.16), P = 0.467], ICU mortality in sepsis [OR 1.00 (0.74-1.37), P = 0.987], risk of pneumonia requiring ICU admission [OR 1.05 (0.98-1.14), P = 0.181], or ICU mortality in pneumonia [OR 1.17 (0.98-1.39), P = 0.079]. Meta-analysis of hospital-admitted and ICU-admitted patients separately yielded similar results [OR 0.98 (0.95-1.01), P = 0.23, and OR 1.05 (0.99-1.12), P = 0.86, respectively]. Among patients with sepsis and severe pneumonia, there was no protective effect of different levels of sUMOD. Results were consistent regardless of geographic origins and not modified by disease severity. NEW & NOTEWORTHY The presence of acute kidney injury in severe infections increases the likelihood of poor outcome severalfold. A decrease in serum uromodulin (sUMOD), synthetized in the kidney, has been proposed as a mediator of this effect. Using the Mendelian randomization technique, we tested the hypothesis that increased sUMOD is protective in severe infections. Analyses, however, showed no evidence of a protective effect of higher levels of sUMOD in sepsis or severe pneumonia.

Whole-genome sequencing in diverse subjects identifies genetic correlates of leukocyte traits: The NHLBI TOPMed program.

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs.

Potential causal role of l-glutamine in sickle cell disease painful crises: A Mendelian randomization analysis.

In a recent clinical trial, the metabolite l-glutamine was shown to reduce painful crises in sickle cell disease (SCD) patients. To support this observation and identify other metabolites implicated in SCD clinical heterogeneity, we profiled 129 metabolites in the plasma of 705 SCD patients. We tested correlations between metabolite levels and six SCD-related complications (painful crises, cholecystectomy, retinopathy, leg ulcer, priapism, aseptic necrosis) or estimated glomerular filtration rate (eGFR), and used Mendelian randomization (MR) to assess causality. We found a potential causal relationship between l-glutamine levels and painful crises (N = 1278, odds ratio (OR) [95% confidence interval] = 0.68 [0.52-0.89], P = 0.0048). In two smaller SCD cohorts (N = 299 and 406), the protective effect of l-glutamine was observed (OR = 0.82 [0.50-1.34]), although the MR result was not significant (P = 0.44). We identified 66 significant correlations between the levels of other metabolites and SCD-related complications or eGFR. We tested these correlations for causality using MR analyses and found no significant causal relationship. The baseline levels of quinolinic acid were associated with prospectively ascertained survival in SCD patients, and this effect was dependent on eGFR. Metabolomics provide a promising approach to prioritize small molecules that may serve as biomarkers or drug targets in SCD.

Genome-wide association study of erythrocyte density in sickle cell disease patients.

Deoxy-hemoglobin S polymerization into rigid fibers is the direct cause of the clinical sequelae observed in sickle cell disease (SCD). The rate of polymerization of sickle hemoglobin is determined primarily by intracellular hemoglobin concentration, itself dependent on the amount of sickle hemoglobin and on red blood cell (RBC) volume. Dense, dehydrated RBC (DRBC) are observed in SCD patients, and their number correlates with hemolytic parameters and complications such as renal dysfunction, leg ulcers and priapism. To identify new genes involved in RBC hydration in SCD, we performed the first genome-wide association study for DRBC in 374 sickle cell anemia (HbSS) patients. We did not find genome-wide significant results, indicating that variants that modulate DRBC have modest-to-weak effects. A secondary analysis demonstrated a nominal association (P=0.003) between DRBC in SCD patients and a variant associated with mean corpuscular hemoglobin concentration (MCHC) in non-anemic individuals. This intronic variant controls the expression of ATP2B4, the main calcium pump in erythrocytes. Our study highlights ATP2B4 as a promising target for modulation of RBC hydration in SCD patients.

Vitamin D, Cognition, and Alzheimer's Disease: Observational and Two-Sample Mendelian Randomization Studies.

Background: Observational studies have found that vitamin D supplementation is associated with improved cognition. Further, recent Mendelian randomization (MR) studies have shown that higher vitamin D levels, 25(OH)D, may protect against Alzheimer's disease. Thus, it is possible that 25(OH)D may protect against Alzheimer's disease by improving cognition. Objective: We assessed this hypothesis, by examining the relationship between 25(OH)D levels and seven cognitive measurements. Methods: To mitigate bias from confounding, we performed two-sample MR analyses. We used instruments from three publications: Manousaki et al. (2020), Sutherland et al. (2022), and the Emerging Risk Factors Collaboration/EPIC-CVD/Vitamin D Studies Collaboration (2021). Results: Our observational studies suggested a protective association between 25(OH)D levels and cognitive measures. An increase in the natural log of 25(OH)D by 1 SD was associated with a higher PACC score (BetaPACC score = 0.06, 95% CI = (0.04-0.08); p = 1.8×10-10). However, in the MR analyses, the estimated effect of 25(OH)D on cognitive measures was null. Specifically, per 1 SD increase in genetically estimated natural log of 25(OH)D, the PACC scores remained unchanged in the overall population, (BetaPACC score = -0.01, 95% CI (-0.06 to 0.03); p = 0.53), and amongst individuals aged over 60 (BetaPACC score = 0.03, 95% CI (-0.028 to 0.08); p = 0.35). Conclusions: In conclusion, our MR study found no clear evidence to support a protective role of increased 25(OH)D concentrations on cognitive performance in European ancestry individuals. However, our study cannot entirely dismiss the potential beneficial effect on PACC for individuals over the age of 60.

Circulating Metabolite Abundances Associated With Risks of Bipolar Disorder, Schizophrenia, and Depression: A Mendelian Randomization Study.

BACKGROUND: Preventive measures and treatments for psychiatric disorders are limited. Circulating metabolites are potential candidates for biomarker and therapeutic target identification, given their measurability and essential roles in biological processes. METHODS: Leveraging large-scale genome-wide association studies, we conducted Mendelian randomization analyses to assess the associations between circulating metabolite abundances and the risks of bipolar disorder, schizophrenia, and depression. Genetic instruments were selected for 94 metabolites measured in the Canadian Longitudinal Study on Aging cohort (N = 8299). We repeated Mendelian randomization analyses based on the UK Biobank, INTERVAL, and EPIC (European Prospective Investigation into Cancer)-Norfolk studies. RESULTS: After validating Mendelian randomization assumptions and colocalization evidence, we found that a 1 SD increase in genetically predicted circulating abundances of eicosapentaenoate and docosapentaenoate was associated with odds ratios of 0.72 (95% CI, 0.65-0.79) and 0.63 (95% CI, 0.55-0.72), respectively, for bipolar disorder. Genetically increased Ω-3 unsaturated fatty acids abundance and Ω-3-to-total fatty acids ratio, as well as genetically decreased Ω-6-to-Ω-3 ratio, were negatively associated with the risk of bipolar disorder in the UK Biobank. Genetically increased circulating abundances of 3 N-acetyl-amino acids were associated with an increased risk of schizophrenia with a maximum odds ratio of 1.31 (95% CI, 1.18-1.44) per 1 SD increase. Furthermore, a 1 SD increase in genetically predicted circulating abundance of hypotaurine was associated with an odds ratio of 0.85 (95% CI, 0.78-0.93) for depression. CONCLUSIONS: The biological mechanisms that underlie Ω-3 unsaturated fatty acids, NAT8-catalyzed N-acetyl-amino acids, and hypotaurine warrant exploration to identify new biomarkers and potential therapeutic targets.

The influence of trinucleotide repeats in the androgen receptor gene on androgen-related traits and diseases.

CONTEXT: Trinucleotide repeats in the androgen receptor have been proposed to influence testosterone signaling in men, but the clinical relevance of these trinucleotide repeats remains controversial. OBJECTIVE: To examine how androgen receptor trinucleotide repeat lengths affect androgen-related traits and disease risks and whether they influence the clinical importance of circulating testosterone levels. METHODS: We quantified CAG and GGC repeat lengths in the androgen receptor (AR) gene of European-ancestry male participants in UK Biobank from whole-genome and whole-exome sequence data using ExpansionHunter, and tested associations with androgen-related traits and diseases. We also examined whether the associations between testosterone levels and these outcomes were affected by adjustment for the repeat lengths. RESULTS: We successfully quantified the repeat lengths from whole-genome and/or whole-exome sequence data in 181,217 males. Both repeat lengths were shown to be positively associated with circulating total testosterone level and bone mineral density, whereas CAG repeat length was negatively associated with male-pattern baldness, but their effects were relatively small and were not associated with most of the other outcomes. Circulating total testosterone level was associated with various outcomes, but this relationship was not affected by adjustment for the repeat lengths. CONCLUSION: In this large-scale study, we found that longer CAG and GGC repeats in the AR gene influence androgen resistance, elevate circulating testosterone level via a feedback loop and play a role in some androgen-targeted tissues. Generally, however, circulating testosterone level is a more important determinant of androgen action in males than repeat lengths.

Gene-metabolite annotation with shortest reactional distance enhances metabolite genome-wide association studies results.

Studies combining metabolomics and genetics, known as metabolite genome-wide association studies (mGWAS), have provided valuable insights into our understanding of the genetic control of metabolite levels. However, the biological interpretation of these associations remains challenging due to a lack of existing tools to annotate mGWAS gene-metabolite pairs beyond the use of conservative statistical significance threshold. Here, we computed the shortest reactional distance (SRD) based on the curated knowledge of the KEGG database to explore its utility in enhancing the biological interpretation of results from three independent mGWAS, including a case study on sickle cell disease patients. Results show that, in reported mGWAS pairs, there is an excess of small SRD values and that SRD values and p-values significantly correlate, even beyond the standard conservative thresholds. The added-value of SRD annotation is shown for identification of potential false negative hits, exemplified by the finding of gene-metabolite associations with SRD ≤1 that did not reach standard genome-wide significance cut-off. The wider use of this statistic as an mGWAS annotation would prevent the exclusion of biologically relevant associations and can also identify errors or gaps in current metabolic pathway databases. Our findings highlight the SRD metric as an objective, quantitative and easy-to-compute annotation for gene-metabolite pairs that can be used to integrate statistical evidence to biological networks.

Gene-metabolite annotation with shortest reactional distance enhances metabolite genome-wide association studies results.

Metabolite genome-wide association studies (mGWAS) have advanced our understanding of the genetic control of metabolite levels. However, interpreting these associations remains challenging due to a lack of tools to annotate gene-metabolite pairs beyond the use of conservative statistical significance threshold. Here, we introduce the shortest reactional distance (SRD) metric, drawing from the comprehensive KEGG database, to enhance the biological interpretation of mGWAS results. We applied this approach to three independent mGWAS, including a case study on sickle cell disease patients. Our analysis reveals an enrichment of small SRD values in reported mGWAS pairs, with SRD values significantly correlating with mGWAS p values, even beyond the standard conservative thresholds. We demonstrate the utility of SRD annotation in identifying potential false negatives and inaccuracies within current metabolic pathway databases. Our findings highlight the SRD metric as an objective, quantitative and easy-to-compute annotation for gene-metabolite pairs, suitable to integrate statistical evidence to biological networks.

HLA allele-calling using multi-ancestry whole-exome sequencing from the UK Biobank identifies 129 novel associations in 11 autoimmune diseases.

The human leukocyte antigen (HLA) region on chromosome 6 is strongly associated with many immune-mediated and infection-related diseases. Due to its highly polymorphic nature and complex linkage disequilibrium patterns, traditional genetic association studies of single nucleotide polymorphisms do not perform well in this region. Instead, the field has adopted the assessment of the association of HLA alleles (i.e., entire HLA gene haplotypes) with disease. Often based on genotyping arrays, these association studies impute HLA alleles, decreasing accuracy and thus statistical power for rare alleles and in non-European ancestries. Here, we use whole-exome sequencing (WES) from 454,824 UK Biobank (UKB) participants to directly call HLA alleles using the HLA-HD algorithm. We show this method is more accurate than imputing HLA alleles and harness the improved statistical power to identify 360 associations for 11 auto-immune phenotypes (at least 129 likely novel), leading to better insights into the specific coding polymorphisms that underlie these diseases. We show that HLA alleles with synonymous variants, often overlooked in HLA studies, can significantly influence these phenotypes. Lastly, we show that HLA sequencing may improve polygenic risk scores accuracy across ancestries. These findings allow better characterization of the role of the HLA region in human disease.

Proteome-wide Mendelian randomization implicates nephronectin as an actionable mediator of the effect of obesity on COVID-19 severity.

Obesity is a major risk factor for Coronavirus disease (COVID-19) severity; however, the mechanisms underlying this relationship are not fully understood. As obesity influences the plasma proteome, we sought to identify circulating proteins mediating the effects of obesity on COVID-19 severity in humans. Here, we screened 4,907 plasma proteins to identify proteins influenced by body mass index using Mendelian randomization. This yielded 1,216 proteins, whose effect on COVID-19 severity was assessed, again using Mendelian randomization. We found that an s.d. increase in nephronectin (NPNT) was associated with increased odds of critically ill COVID-19 (OR = 1.71, P = 1.63 × 10-10). The effect was driven by an NPNT splice isoform. Mediation analyses supported NPNT as a mediator. In single-cell RNA-sequencing, NPNT was expressed in alveolar cells and fibroblasts of the lung in individuals who died of COVID-19. Finally, decreasing body fat mass and increasing fat-free mass were found to lower NPNT levels. These findings provide actionable insights into how obesity influences COVID-19 severity.

Variant-aware saturating mutagenesis using multiple Cas9 nucleases identifies regulatory elements at trait-associated loci.

Cas9-mediated, high-throughput, saturating in situ mutagenesis permits fine-mapping of function across genomic segments. Disease- and trait-associated variants identified in genome-wide association studies largely cluster at regulatory loci. Here we demonstrate the use of multiple designer nucleases and variant-aware library design to interrogate trait-associated regulatory DNA at high resolution. We developed a computational tool for the creation of saturating-mutagenesis libraries with single or multiple nucleases with incorporation of variants. We applied this methodology to the HBS1L-MYB intergenic region, which is associated with red-blood-cell traits, including fetal hemoglobin levels. This approach identified putative regulatory elements that control MYB expression. Analysis of genomic copy number highlighted potential false-positive regions, thus emphasizing the importance of off-target analysis in the design of saturating-mutagenesis experiments. Together, these data establish a widely applicable high-throughput and high-resolution methodology to identify minimal functional sequences within large disease- and trait-associated regions.