A cross-disorder dosage sensitivity map of the human genome.

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the properties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

Multilayered genetic and omics dissection of mitochondrial activity in a mouse reference population.

The manner by which genotype and environment affect complex phenotypes is one of the fundamental questions in biology. In this study, we quantified the transcriptome--a subset of the metabolome--and, using targeted proteomics, quantified a subset of the liver proteome from 40 strains of the BXD mouse genetic reference population on two diverse diets. We discovered dozens of transcript, protein, and metabolite QTLs, several of which linked to metabolic phenotypes. Most prominently, Dhtkd1 was identified as a primary regulator of 2-aminoadipate, explaining variance in fasted glucose and diabetes status in both mice and humans. These integrated molecular profiles also allowed further characterization of complex pathways, particularly the mitochondrial unfolded protein response (UPR(mt)). UPR(mt) shows strikingly variant responses at the transcript and protein level that are remarkably conserved among C. elegans, mice, and humans. Overall, these examples demonstrate the value of an integrated multilayered omics approach to characterize complex metabolic phenotypes.

Genetic insights into the age-specific biological mechanisms governing human ovarian aging.

There is currently little evidence that the genetic basis of human phenotype varies significantly across the lifespan. However, time-to-event phenotypes are understudied and can be thought of as reflecting an underlying hazard, which is unlikely to be constant through life when values take a broad range. Here, we find that 74% of 245 genome-wide significant genetic associations with age at natural menopause (ANM) in the UK Biobank show a form of age-specific effect. Nineteen of these replicated discoveries are identified only by our modeling framework, which determines the time dependency of DNA-variant age-at-onset associations without a significant multiple-testing burden. Across the range of early to late menopause, we find evidence for significantly different underlying biological pathways, changes in the signs of genetic correlations of ANM to health indicators and outcomes, and differences in inferred causal relationships. We find that DNA damage response processes only act to shape ovarian reserve and depletion for women of early ANM. Genetically mediated delays in ANM were associated with increased relative risk of breast cancer and leiomyoma at all ages and with high cholesterol and heart failure for late-ANM women. These findings suggest that a better understanding of the age dependency of genetic risk factor relationships among health indicators and outcomes is achievable through appropriate statistical modeling of large-scale biobank data.

The impact of 22q11.2 copy-number variants on human traits in the general population.

While extensively studied in clinical cohorts, the phenotypic consequences of 22q11.2 copy-number variants (CNVs) in the general population remain understudied. To address this gap, we performed a phenome-wide association scan in 405,324 unrelated UK Biobank (UKBB) participants by using CNV calls from genotyping array. We mapped 236 Human Phenotype Ontology terms linked to any of the 90 genes encompassed by the region to 170 UKBB traits and assessed the association between these traits and the copy-number state of 504 genotyping array probes in the region. We found significant associations for eight continuous and nine binary traits associated under different models (duplication-only, deletion-only, U-shape, and mirror models). The causal effect of the expression level of 22q11.2 genes on associated traits was assessed through transcriptome-wide Mendelian randomization (TWMR), revealing that increased expression of ARVCF increased BMI. Similarly, increased DGCR6 expression causally reduced mean platelet volume, in line with the corresponding CNV effect. Furthermore, cross-trait multivariable Mendelian randomization (MVMR) suggested a predominant role of genuine (horizontal) pleiotropy in the CNV region. Our findings show that within the general population, 22q11.2 CNVs are associated with traits previously linked to genes in the region, and duplications and deletions act upon traits in different fashions. We also showed that gain or loss of distinct segments within 22q11.2 may impact a trait under different association models. Our results have provided new insights to help further the understanding of the complex 22q11.2 region.

Liability-scale heritability estimation for biobank studies of low-prevalence disease.

Theory for liability-scale models of the underlying genetic basis of complex disease provides an important way to interpret, compare, and understand results generated from biological studies. In particular, through estimation of the liability-scale heritability (LSH), liability models facilitate an understanding and comparison of the relative importance of genetic and environmental risk factors that shape different clinically important disease outcomes. Increasingly, large-scale biobank studies that link genetic information to electronic health records, containing hundreds of disease diagnosis indicators that mostly occur infrequently within the sample, are becoming available. Here, we propose an extension of the existing liability-scale model theory suitable for estimating LSH in biobank studies of low-prevalence disease. In a simulation study, we find that our derived expression yields lower mean square error (MSE) and is less sensitive to prevalence misspecification as compared to previous transformations for diseases with ≤2% population prevalence and LSH of ≤0.45, especially if the biobank sample prevalence is less than that of the wider population. Applying our expression to 13 diagnostic outcomes of ≤3% prevalence in the UK Biobank study revealed important differences in LSH obtained from the different theoretical expressions that impact the conclusions made when comparing LSH across disease outcomes. This demonstrates the importance of careful consideration for estimation and prediction of low-prevalence disease outcomes and facilitates improved inference of the underlying genetic basis of ≤2% population prevalence diseases, especially where biobank sample ascertainment results in a healthier sample population.

The individual and global impact of copy-number variants on complex human traits.

The impact of copy-number variations (CNVs) on complex human traits remains understudied. We called CNVs in 331,522 UK Biobank participants and performed genome-wide association studies (GWASs) between the copy number of CNV-proxy probes and 57 continuous traits, revealing 131 signals spanning 47 phenotypes. Our analysis recapitulated well-known associations (e.g., 1q21 and height), revealed the pleiotropy of recurrent CNVs (e.g., 26 and 16 traits for 16p11.2-BP4-BP5 and 22q11.21, respectively), and suggested gene functionalities (e.g., MARF1 in female reproduction). Forty-eight CNV signals (38%) overlapped with single-nucleotide polymorphism (SNP)-GWASs signals for the same trait. For instance, deletion of PDZK1, which encodes a urate transporter scaffold protein, decreased serum urate levels, while deletion of RHD, which encodes the Rhesus blood group D antigen, associated with hematological traits. Other signals overlapped Mendelian disorder regions, suggesting variable expressivity and broad impact of these loci, as illustrated by signals mapping to Rotor syndrome (SLCO1B1/3), renal cysts and diabetes syndrome (HNF1B), or Charcot-Marie-Tooth (PMP22) loci. Total CNV burden negatively impacted 35 traits, leading to increased adiposity, liver/kidney damage, and decreased intelligence and physical capacity. Thirty traits remained burden associated after correcting for CNV-GWAS signals, pointing to a polygenic CNV architecture. The burden negatively correlated with socio-economic indicators, parental lifespan, and age (survivorship proxy), suggesting a contribution to decreased longevity. Together, our results showcase how studying CNVs can expand biological insights, emphasizing the critical role of this mutational class in shaping human traits and arguing in favor of a continuum between Mendelian and complex diseases.

Using genetic variation to disentangle the complex relationship between food intake and health outcomes.

Diet is considered as one of the most important modifiable factors influencing human health, but efforts to identify foods or dietary patterns associated with health outcomes often suffer from biases, confounding, and reverse causation. Applying Mendelian randomization in this context may provide evidence to strengthen causality in nutrition research. To this end, we first identified 283 genetic markers associated with dietary intake in 445,779 UK Biobank participants. We then converted these associations into direct genetic effects on food exposures by adjusting them for effects mediated via other traits. The SNPs which did not show evidence of mediation were then used for MR, assessing the association between genetically predicted food choices and other risk factors, health outcomes. We show that using all associated SNPs without omitting those which show evidence of mediation, leads to biases in downstream analyses (genetic correlations, causal inference), similar to those present in observational studies. However, MR analyses using SNPs which have only a direct effect on the exposure on food exposures provided unequivocal evidence of causal associations between specific eating patterns and obesity, blood lipid status, and several other risk factors and health outcomes.

Bias correction for inverse variance weighting Mendelian randomization.

Inverse-variance weighted two-sample Mendelian randomization (IVW-MR) is the most widely used approach that utilizes genome-wide association studies (GWAS) summary statistics to infer the existence and the strength of the causal effect between an exposure and an outcome. Estimates from this approach can be subject to different biases due to the use of weak instruments and winner's curse, which can change as a function of the overlap between the exposure and outcome samples. We developed a method (MRlap) that simultaneously considers weak instrument bias and winner's curse while accounting for potential sample overlap. Assuming spike-and-slab genomic architecture and leveraging linkage disequilibrium score regression and other techniques, we could analytically derive, reliably estimate, and hence correct for the bias of IVW-MR using association summary statistics only. We tested our approach using simulated data for a wide range of realistic settings. In all the explored scenarios, our correction reduced the bias, in some situations by as much as 30-fold. In addition, our results are consistent with the fact that the strength of the biases will decrease as the sample size increases and we also showed that the overall bias is also dependent on the genetic architecture of the exposure, and traits with low heritability and/or high polygenicity are more strongly affected. Applying MRlap to obesity-related exposures revealed statistically significant differences between IVW-based and corrected effects, both for nonoverlapping and fully overlapping samples. Our method not only reduces bias in causal effect estimation but also enables the use of much larger GWAS sample sizes, by allowing for potentially overlapping samples.

Identification of four novel loci associated with psychotropic drug-induced weight gain in a Swiss psychiatric longitudinal study: A GWAS analysis.

Patients suffering from mental disorders are at high risk of developing cardiovascular diseases, leading to a reduction in life expectancy. Genetic variants can display greater influence on cardiometabolic features in psychiatric cohorts compared to the general population. The difference is possibly due to an intricate interaction between the mental disorder or the medications used to treat it and metabolic regulations. Previous genome wide association studies (GWAS) on antipsychotic-induced weight gain included a low number of participants and/or were restricted to patients taking one specific antipsychotic. We conducted a GWAS of the evolution of body mass index (BMI) during early (i.e., ≤ 6) months of treatment with psychotropic medications inducing metabolic disturbances (i.e., antipsychotics, mood stabilizers and some antidepressants) in 1135 patients from the PsyMetab cohort. Six highly correlated BMI phenotypes (i.e., BMI change and BMI slope after distinct durations of psychotropic treatment) were considered in the analyses. Our results showed that four novel loci were associated with altered BMI upon treatment at genome-wide significance (p < 5 × 10-8): rs7736552 (near MAN2A1), rs11074029 (in SLCO3A1), rs117496040 (near DEFB1) and rs7647863 (in IQSEC1). Associations between the four loci and alternative BMI-change phenotypes showed consistent effects. Replication analyses in 1622 UK Biobank participants under psychotropic treatment showed a consistent association between rs7736552 and BMI slope (p = 0.017). These findings provide new insights into metabolic side effects induced by psychotropic drugs and underline the need for future studies to replicate these associations in larger cohorts.

Influence of Genetic Ancestry on Human Serum Proteome.

Disease risk varies significantly between ethnic groups, however, the clinical significance and implications of these observations are poorly understood. Investigating ethnic differences within the human proteome may shed light on the impact of ancestry on disease risk. We used admixture mapping to explore the impact of genetic ancestry on 237 cardiometabolic biomarkers in 2,216 Latin Americans within the Outcomes Reduction with an Initial Glargine Intervention (ORIGIN) study. We developed a variance component model in order to determine the proportion of variance explained by inter-ancestry differences, and we applied it to the biomarker panel. Multivariable linear regression was used to identify and localize genetic loci affecting biomarker variability between ethnicities. Variance component analysis revealed that 5% of biomarkers were significantly impacted by genetic admixture (p < 0.05/237), including C-peptide, apolipoprotein-E, and intercellular adhesion molecule 1. We also identified 46 regional associations across 40 different biomarkers (p < 1.13 × 10-6). An independent analysis revealed that 34 of these 46 regions were associated at genome-wide significance (p < 5 × 10-8) with their respective biomarker in either Europeans or Latin populations. Additional analyses revealed that an admixture mapping signal associated with increased C-peptide levels was also associated with an increase in diabetes risk (odds ratio [OR] = 6.07 per SD, 95% confidence interval [CI] 1.44 to 25.56, p = 0.01) and surrogate measures of insulin resistance. Our results demonstrate the impact of ancestry on biomarker levels, suggesting that some of the observed differences in disease prevalence have a biological basis, and that reference intervals for those biomarkers should be tailored to ancestry. Specifically, our results point to a strong role of ancestry in insulin resistance and diabetes risk.

Gene set enrichment analysis of pathophysiological pathways highlights oxidative stress in psychosis.

Polygenic risk prediction remains an important aim of genetic association studies. Currently, the predictive power of schizophrenia polygenic risk scores (PRSs) is not large enough to allow highly accurate discrimination between cases and controls and thus is not adequate for clinical integration. Since PRSs are rarely used to reveal biological functions or to validate candidate pathways, to fill this gap, we investigated whether their predictive ability could be improved by building genome-wide (GW-PRSs) and pathway-specific PRSs, using distance- or expression quantitative trait loci (eQTLs)- based mapping between genetic variants and genes. We focused on five pathways (glutamate, oxidative stress, GABA/interneurons, neuroimmune/neuroinflammation and myelin) which belong to a critical hub of schizophrenia pathophysiology, centred on redox dysregulation/oxidative stress. Analyses were first performed in the Lausanne Treatment and Early Intervention in Psychosis Program (TIPP) study (n = 340, cases/controls: 208/132), a sample of first-episode of psychosis patients and matched controls, and then validated in an independent study, the epidemiological and longitudinal intervention program of First-Episode Psychosis in Cantabria (PAFIP) (n = 352, 224/128). Our results highlighted two main findings. First, GW-PRSs for schizophrenia were significantly associated with early psychosis status. Second, oxidative stress was the only significantly associated pathway that showed an enrichment in both the TIPP (p = 0.03) and PAFIP samples (p = 0.002), and exclusively when gene-variant linking was done using eQTLs. The results suggest that the predictive accuracy of polygenic risk scores could be improved with the inclusion of information from functional annotations, and through a focus on specific pathways, emphasizing the need to build and study functionally informed risk scores.

Cannabis use and atherosclerotic cardiovascular disease: a Mendelian randomization study.

BACKGROUND: Association between cannabis use and development of atherosclerotic cardiovascular disease (ASCVD) is inconsistent and challenging to interpret, given existing study limitations. METHODS: Sixty five independent single-nucleotide polymorphisms (SNPs), obtained from a genome-wide association study on lifetime cannabis use, were employed as genetic instruments to estimate the effects of genetically indexed cannabis use on risk of coronary artery disease (CAD) and acute ischemic stroke (IS) using a two-sample Mendelian randomization (MR) approach. Summary statistics on CAD (CARDIoGRAMplusC4D; 60,801 cases and 123,504 controls) and IS (MEGASTROKE; 34,217 cases and 406,111 controls) were obtained separately. A comprehensive review of the observational literature on cannabis use and CAD or IS was also performed and contrasted with MR results. RESULTS: There was no causal effect of cannabis use on the risk of CAD (odds ratio (OR) per ever-users vs. never-users 0.93; 95% confidence interval (CI), 0.83 to 1.03) or IS (OR 1.05; 95%CI, 0.93 to 1.19). Sensitivity analyses yielded similar results, and no heterogeneity and directional pleiotropy was observed. Our meta-analysis of observational studies showed no significant association between ever use of cannabis with risk of CAD (k = 6 studies; ORpooled = 1.23, 95%CI 0.78 to 1.69), nor with IS (k = 6 studies; ORpooled = 1.22, 95%CI 0.95 to 1.50). CONCLUSION: Using a genetic approach approximating a clinical trial does not provide evidence consistent with a causal effect of genetic predisposition to cannabis use on CAD or IS development. Further studies are needed to replicate our findinds, an to investigate more precisely the risk of ASCVD in relation to the quantity, type, route of administration, or the age at exposure to cannabis.

The Human-Specific BOLA2 Duplication Modifies Iron Homeostasis and Anemia Predisposition in Chromosome 16p11.2 Autism Individuals.

Human-specific duplications at chromosome 16p11.2 mediate recurrent pathogenic 600 kbp BP4-BP5 copy-number variations, which are among the most common genetic causes of autism. These copy-number polymorphic duplications are under positive selection and include three to eight copies of BOLA2, a gene involved in the maturation of cytosolic iron-sulfur proteins. To investigate the potential advantage provided by the rapid expansion of BOLA2, we assessed hematological traits and anemia prevalence in 379,385 controls and individuals who have lost or gained copies of BOLA2: 89 chromosome 16p11.2 BP4-BP5 deletion carriers and 56 reciprocal duplication carriers in the UK Biobank. We found that the 16p11.2 deletion is associated with anemia (18/89 carriers, 20%, p = 4e-7, OR = 5), particularly iron-deficiency anemia. We observed similar enrichments in two clinical 16p11.2 deletion cohorts, which included 6/63 (10%) and 7/20 (35%) unrelated individuals with anemia, microcytosis, low serum iron, or low blood hemoglobin. Upon stratification by BOLA2 copy number, our data showed an association between low BOLA2 dosage and the above phenotypes (8/15 individuals with three copies, 53%, p = 1e-4). In parallel, we analyzed hematological traits in mice carrying the 16p11.2 orthologous deletion or duplication, as well as Bola2+/- and Bola2-/- animals. The Bola2-deficient mice and the mice carrying the deletion showed early evidence of iron deficiency, including a mild decrease in hemoglobin, lower plasma iron, microcytosis, and an increased red blood cell zinc-protoporphyrin-to-heme ratio. Our results indicate that BOLA2 participates in iron homeostasis in vivo, and its expansion has a potential adaptive role in protecting against iron deficiency.

Exome-Derived Adiponectin-Associated Variants Implicate Obesity and Lipid Biology.

Circulating levels of adiponectin, an adipocyte-secreted protein associated with cardiovascular and metabolic risk, are highly heritable. To gain insights into the biology that regulates adiponectin levels, we performed an exome array meta-analysis of 265,780 genetic variants in 67,739 individuals of European, Hispanic, African American, and East Asian ancestry. We identified 20 loci associated with adiponectin, including 11 that had been reported previously (p < 2 × 10-7). Comparison of exome array variants to regional linkage disequilibrium (LD) patterns and prior genome-wide association study (GWAS) results detected candidate variants (r2 > .60) spanning as much as 900 kb. To identify potential genes and mechanisms through which the previously unreported association signals act to affect adiponectin levels, we assessed cross-trait associations, expression quantitative trait loci in subcutaneous adipose, and biological pathways of nearby genes. Eight of the nine loci were also associated (p < 1 × 10-4) with at least one obesity or lipid trait. Candidate genes include PRKAR2A, PTH1R, and HDAC9, which have been suggested to play roles in adipocyte differentiation or bone marrow adipose tissue. Taken together, these findings provide further insights into the processes that influence circulating adiponectin levels.

Assessment of network module identification across complex diseases.

Many bioinformatics methods have been proposed for reducing the complexity of large gene or protein networks into relevant subnetworks or modules. Yet, how such methods compare to each other in terms of their ability to identify disease-relevant modules in different types of network remains poorly understood. We launched the 'Disease Module Identification DREAM Challenge', an open competition to comprehensively assess module identification methods across diverse protein-protein interaction, signaling, gene co-expression, homology and cancer-gene networks. Predicted network modules were tested for association with complex traits and diseases using a unique collection of 180 genome-wide association studies. Our robust assessment of 75 module identification methods reveals top-performing algorithms, which recover complementary trait-associated modules. We find that most of these modules correspond to core disease-relevant pathways, which often comprise therapeutic targets. This community challenge establishes biologically interpretable benchmarks, tools and guidelines for molecular network analysis to study human disease biology.

Haematological changes from conception to childbirth: An indicator of major pregnancy complications.

BACKGROUND: About 800 women die every day worldwide from pregnancy-related complications, including excessive blood loss, infections and high-blood pressure (World Health Organization, 2019). To improve screening for high-risk pregnancies, we set out to identify patterns of maternal hematological changes associated with future pregnancy complications. METHODS: Using mixed effects models, we established changes in 14 complete blood count (CBC) parameters for 1710 healthy pregnancies and compared them to measurements from 98 pregnancy-induced hypertension, 106 gestational diabetes and 339 postpartum hemorrhage cases. RESULTS: Results show interindividual variations, but good individual repeatability in CBC values during physiological pregnancies, allowing the identification of specific alterations in women with obstetric complications. For example, in women with uncomplicated pregnancies, haemoglobin count decreases of 0.12 g/L (95% CI -0.16, -0.09) significantly per gestation week (p value <.001). Interestingly, this decrease is three times more pronounced in women who will develop pregnancy-induced hypertension, with an additional decrease of 0.39 g/L (95% CI -0.51, -0.26). We also confirm that obstetric complications and white CBC predict the likelihood of giving birth earlier during pregnancy. CONCLUSION: We provide a comprehensive description of the associations between haematological changes through pregnancy and three major obstetric complications to support strategies for prevention, early-diagnosis and maternal care.

Untargeted Metabolome- and Transcriptome-Wide Association Study Suggests Causal Genes Modulating Metabolite Concentrations in Urine.

Gene products can affect the concentrations of small molecules (aka "metabolites"), and conversely, some metabolites can modulate the concentrations of gene transcripts. While many specific instances of this interplay have been revealed, a global approach to systematically uncover human gene-metabolite interactions is still lacking. We performed a metabolome- and transcriptome-wide association study to identify genes influencing the human metabolome using untargeted metabolome features, extracted from 1H nuclear magnetic resonance spectroscopy (NMR) of urine samples, and gene expression levels, quantified from RNA-Seq of lymphoblastoid cell lines (LCL) from 555 healthy individuals. We identified 20 study-wide significant associations corresponding to 15 genes, of which 5 associations (with 2 genes) were confirmed with follow-up NMR data. Using metabomatching, we identified the metabolites corresponding to metabolome features associated with the genes, namely, N-acetylated compounds with ALMS1 and trimethylamine (TMA) with HPS1. Finally, Mendelian randomization analysis supported a potential causal link between the expression of genes in both the ALMS1- and HPS1-loci and their associated metabolite concentrations. In the case of HPS1, we additionally observed that TMA concentration likely exhibits a reverse causal effect on HPS1 expression levels, indicating a negative feedback loop. Our study highlights how the integration of metabolomics, gene expression, and genetic data can pinpoint causal genes modulating metabolite concentrations.

FADS3 is a Δ14Z sphingoid base desaturase that contributes to gender differences in the human plasma sphingolipidome.

Sphingolipids (SLs) are structurally diverse lipids that are defined by the presence of a long-chain base (LCB) backbone. Typically, LCBs contain a single Δ4E double bond (DB) (mostly d18:1), whereas the dienic LCB sphingadienine (d18:2) contains a second DB at the Δ14Z position. The enzyme introducing the Δ14Z DB is unknown. We analyzed the LCB plasma profile in a gender-, age-, and BMI-matched subgroup of the CoLaus cohort (n = 658). Sphingadienine levels showed a significant association with gender, being on average ∼30% higher in females. A genome-wide association study (GWAS) revealed variants in the fatty acid desaturase 3 (FADS3) gene to be significantly associated with the plasma d18:2/d18:1 ratio (p = -log 7.9). Metabolic labeling assays, FADS3 overexpression and knockdown approaches, and plasma LCB profiling in FADS3-deficient mice confirmed that FADS3 is a bona fide LCB desaturase and required for the introduction of the Δ14Z double bond. Moreover, we showed that FADS3 is required for the conversion of the atypical cytotoxic 1-deoxysphinganine (1-deoxySA, m18:0) to 1-deoxysphingosine (1-deoxySO, m18:1). HEK293 cells overexpressing FADS3 were more resistant to m18:0 toxicity than WT cells. In summary, using a combination of metabolic profiling and GWAS, we identified FADS3 to be essential for forming Δ14Z DB containing LCBs, such as d18:2 and m18:1. Our results unravel FADS3 as a Δ14Z LCB desaturase, thereby disclosing the last missing enzyme of the SL de novo synthesis pathway.

Mendelian randomization to assess causality between uromodulin, blood pressure and chronic kidney disease.

UMOD variants associated with higher levels of urinary uromodulin (uUMOD) increase the risk of chronic kidney disease (CKD) and hypertension. However, uUMOD levels also reflect functional kidney tubular mass in observational studies, questioning the causal link between uromodulin production and kidney damage. We used Mendelian randomization to clarify causality between uUMOD levels, kidney function and blood pressure in individuals of European descent. The link between uUMOD and estimated glomerular filtration rate (eGFR) was first investigated in a population-based cohort of 3851 individuals. In observational data, higher uUMOD associated with higher eGFR. Conversely, when using rs12917707 (an UMOD polymorphism) as an instrumental variable in one-sample Mendelian randomization, higher uUMOD strongly associated with eGFR decline. We next applied two-sample Mendelian randomization on four genome wide association study consortia to explore causal links between uUMOD and eGFR, CKD risk (567,460 individuals) and blood pressure (757,461 individuals). Higher uUMOD levels significantly associated with lower eGFR, higher odds for eGFR decline or CKD, and higher systolic or diastolic blood pressure. Each one standard deviation (SD) increase of uUMOD decreased log-transformed eGFR by -0.15 SD (95% confidence interval -0.17 to -0.13) and increased log-odds CKD by 0.13 SD (0.12 to 0.15). One SD increase of uUMOD increased systolic blood pressure by 0.06 SD (0.03 to 0.09) and diastolic blood pressure by 0.08 SD (0.05 to 0.12). The effect of uUMOD on blood pressure was mediated by eGFR, whereas the effect on eGFR was not mediated by blood pressure. Thus, our data support that genetically driven levels of uromodulin have a direct, causal and adverse effect on kidney function outcome in the general population, not mediated by blood pressure.