Genetic drivers of heterogeneity in type 2 diabetes pathophysiology.

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes1,2 and molecular mechanisms that are often specific to cell type3,4. Here, to characterize the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study data from 2,535,601 individuals (39.7% not of European ancestry), including 428,452 cases of T2D. We identify 1,289 independent association signals at genome-wide significance (P < 5 × 10-8) that map to 611 loci, of which 145 loci are, to our knowledge, previously unreported. We define eight non-overlapping clusters of T2D signals that are characterized by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type-specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial cells and enteroendocrine cells. We build cluster-specific partitioned polygenic scores5 in a further 279,552 individuals of diverse ancestry, including 30,288 cases of T2D, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned polygenic scores are associated with coronary artery disease, peripheral artery disease and end-stage diabetic nephropathy across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings show the value of integrating multi-ancestry genome-wide association study data with single-cell epigenomics to disentangle the aetiological heterogeneity that drives the development and progression of T2D. This might offer a route to optimize global access to genetically informed diabetes care.

Lifeday coverage of oral anticoagulants and one-year relative survival in patients with atrial fibrillation: a population-based study in Estonia.

BACKGROUND: Routine oral anticoagulation (OAC) is recommended for almost all high-risk patients with atrial fibrillation, yet registries show that OACs are still underused. Our aim was to study the lifeday coverage (LDC) of OAC prescriptions and its relationship with one-year mortality rates of AF patients aged ≥ 65 in Estonia for the years 2019 and 2020. METHODS: Medical data for AF patients aged ≥ 65 years from 2018 and alive as of 01.01.2019 (cohort I) and new AF documentation from 2019 and alive as of 01.01.2020 (cohort II) was obtained from the Health Insurance Fund's electronic database. The data was linked to the nationwide Estonian Medical Prescription Centre's database of prescribed OACs. For LDC analysis, daily doses of guideline-recommended OACs were used. The patients were categorized into three LDC groups: 0%, 1-79%, and ≥ 80%. The data was linked to the Estonian Causes of Death Registry to establish the date of death and mortality rate for the whole Estonian population aged ≥ 65. RESULTS: There were 34,018 patients in cohort I and 9,175 patients with new AF documentation (cohort II), previously not included in cohort I. Of the patients, 77.7% and 68.6% had at least one prescription of OAC in cohorts I and II respectively. 57.4% in cohort I and 44.5% in cohort II had an LDC of ≥ 80%. The relative survival estimates at 1 year for LDC lifeday coverage groups 0%, 1-79%, and ≥ 80% were 91.2%, 98.2%, and 98.5% (cohort I), and 91.9%, 95.2%, and 97.6% (cohort II), respectively. CONCLUSIONS: Despite clear indications for OAC use, LDC is still insufficient and anticoagulation is underused for stroke prevention in Estonia. Further education of the medical community and patients is needed to achieve higher lifeday coverage of prescribed OACs.

Comprehensive population-based genome sequencing provides insight into hematopoietic regulatory mechanisms.

Genetic variants affecting hematopoiesis can influence commonly measured blood cell traits. To identify factors that affect hematopoiesis, we performed association studies for blood cell traits in the population-based Estonian Biobank using high-coverage whole-genome sequencing (WGS) in 2,284 samples and SNP genotyping in an additional 14,904 samples. Using up to 7,134 samples with available phenotype data, our analyses identified 17 associations across 14 blood cell traits. Integration of WGS-based fine-mapping and complementary epigenomic datasets provided evidence for causal mechanisms at several loci, including at a previously undiscovered basophil count-associated locus near the master hematopoietic transcription factor CEBPA The fine-mapped variant at this basophil count association near CEBPA overlapped an enhancer active in common myeloid progenitors and influenced its activity. In situ perturbation of this enhancer by CRISPR/Cas9 mutagenesis in hematopoietic stem and progenitor cells demonstrated that it is necessary for and specifically regulates CEBPA expression during basophil differentiation. We additionally identified basophil count-associated variation at another more pleiotropic myeloid enhancer near GATA2, highlighting regulatory mechanisms for ordered expression of master hematopoietic regulators during lineage specification. Our study illustrates how population-based genetic studies can provide key insights into poorly understood cell differentiation processes of considerable physiologic relevance.

Whole-Exome Sequencing Identifies Loci Associated with Blood Cell Traits and Reveals a Role for Alternative GFI1B Splice Variants in Human Hematopoiesis.

Circulating blood cell counts and indices are important indicators of hematopoietic function and a number of clinical parameters, such as blood oxygen-carrying capacity, inflammation, and hemostasis. By performing whole-exome sequence association analyses of hematologic quantitative traits in 15,459 community-dwelling individuals, followed by in silico replication in up to 52,024 independent samples, we identified two previously undescribed coding variants associated with lower platelet count: a common missense variant in CPS1 (rs1047891, MAF = 0.33, discovery + replication p = 6.38 × 10(-10)) and a rare synonymous variant in GFI1B (rs150813342, MAF = 0.009, discovery + replication p = 1.79 × 10(-27)). By performing CRISPR/Cas9 genome editing in hematopoietic cell lines and follow-up targeted knockdown experiments in primary human hematopoietic stem and progenitor cells, we demonstrate an alternative splicing mechanism by which the GFI1B rs150813342 variant suppresses formation of a GFI1B isoform that preferentially promotes megakaryocyte differentiation and platelet production. These results demonstrate how unbiased studies of natural variation in blood cell traits can provide insight into the regulation of human hematopoiesis.

Translating genotype data of 44,000 biobank participants into clinical pharmacogenetic recommendations: challenges and solutions.

PURPOSE: Biomedical databases combining electronic medical records and phenotypic and genomic data constitute a powerful resource for the personalization of treatment. To leverage the wealth of information provided, algorithms are required that systematically translate the contained information into treatment recommendations based on existing genotype-phenotype associations. METHODS: We developed and tested algorithms for translation of preexisting genotype data of over 44,000 participants of the Estonian biobank into pharmacogenetic recommendations. We compared the results obtained by genome sequencing, exome sequencing, and genotyping using microarrays, and evaluated the impact of pharmacogenetic reporting based on drug prescription statistics in the Nordic countries and Estonia. RESULTS: Our most striking result was that the performance of genotyping arrays is similar to that of genome sequencing, whereas exome sequencing is not suitable for pharmacogenetic predictions. Interestingly, 99.8% of all assessed individuals had a genotype associated with increased risks to at least one medication, and thereby the implementation of pharmacogenetic recommendations based on genotyping affects at least 50 daily drug doses per 1000 inhabitants. CONCLUSION: We find that microarrays are a cost-effective solution for creating preemptive pharmacogenetic reports, and with slight modifications, existing databases can be applied for automated pharmacogenetic decision support for clinicians.

Recall by genotype and cascade screening for familial hypercholesterolemia in a population-based biobank from Estonia.

PURPOSE: Large-scale, population-based biobanks integrating health records and genomic profiles may provide a platform to identify individuals with disease-predisposing genetic variants. Here, we recall probands carrying familial hypercholesterolemia (FH)-associated variants, perform cascade screening of family members, and describe health outcomes affected by such a strategy. METHODS: The Estonian Biobank of Estonian Genome Center, University of Tartu, comprises 52,274 individuals. Among 4776 participants with exome or genome sequences, we identified 27 individuals who carried FH-associated variants in the LDLR, APOB, or PCSK9 genes. Cascade screening of 64 family members identified an additional 20 carriers of FH-associated variants. RESULTS: Via genetic counseling and clinical management of carriers, we were able to reclassify 51% of the study participants from having previously established nonspecific hypercholesterolemia to having FH and identify 32% who were completely unaware of harboring a high-risk disease-associated genetic variant. Imaging-based risk stratification targeted 86% of the variant carriers for statin treatment recommendations. CONCLUSION: Genotype-guided recall of probands and subsequent cascade screening for familial hypercholesterolemia is feasible within a population-based biobank and may facilitate more appropriate clinical management.

Single base resolution analysis of 5-hydroxymethylcytosine in 188 human genes: implications for hepatic gene expression.

To improve the epigenomic analysis of tissues rich in 5-hydroxymethylcytosine (hmC), we developed a novel protocol called TAB-Methyl-SEQ, which allows for single base resolution profiling of both hmC and 5-methylcytosine by targeted next-generation sequencing. TAB-Methyl-SEQ data were extensively validated by a set of five methodologically different protocols. Importantly, these extensive cross-comparisons revealed that protocols based on Tet1-assisted bisulfite conversion provided more precise hmC values than TrueMethyl-based methods. A total of 109 454 CpG sites were analyzed by TAB-Methyl-SEQ for mC and hmC in 188 genes from 20 different adult human livers. We describe three types of variability of hepatic hmC profiles: (i) sample-specific variability at 40.8% of CpG sites analyzed, where the local hmC values correlate to the global hmC content of livers (measured by LC-MS), (ii) gene-specific variability, where hmC levels in the coding regions positively correlate to expression of the respective gene and (iii) site-specific variability, where prominent hmC peaks span only 1 to 3 neighboring CpG sites. Our data suggest that both the gene- and site-specific components of hmC variability might contribute to the epigenetic control of hepatic genes. The protocol described here should be useful for targeted DNA analysis in a variety of applications.

RegScan: a GWAS tool for quick estimation of allele effects on continuous traits and their combinations.

UNLABELLED: Genome-wide association studies are becoming computationally more demanding with the growing amounts of data. Combinatorial traits can increase the data dimensions beyond the computational capabilities of the current tools. We addressed this issue by creating an application for quick association analysis that is ten to hundreds of times faster than the leading fast methods. Our tool (RegScan) is designed for performing basic linear regression analysis with continuous traits maximally fast on large data sets. RegScan specifically targets association analysis of combinatorial traits in metabolomics. It can both generate and analyze the combinatorial traits efficiently. RegScan is capable of analyzing any number of traits together without the need to specify each trait individually. The main goal of the article is to show that RegScan can be the preferred analytical tool when large amounts of data need to be analyzed quickly using the allele frequency test. AVAILABILITY: Precompiled RegScan (all major platforms), source code, user guide and examples are freely available at www.biobank.ee/regscan. REQUIREMENTS: Qt 4.4.3 or newer for dynamic compilations.

Exome analysis in an Estonian multiplex family with neural tube defects-a case report.

INTRODUCTION: Neural tube defects (NTDs) are a group of common and severe congenital birth defects that occur during early embryonic development due to incomplete closure of the neural tube. The genetic architecture of human NTDs, including spina bifida and hydrocephalus, is highly heterogeneous, with multiple genes/loci and both gene-gene and gene-environment interactions involved. Hence, the variation in outcomes also most likely relates to a combination of the severity of different variants in multiple genes and genetic modifiers affecting the biochemical traits. METHODS: Here, we present a multiple-spouse family with one pedigree lineage where three brothers are affected with NTDs-two lumbar spina bifidas without hydrocephalus and one obstructive hydrocephalus. We sequenced the exomes of three NTD patients and their parents. RESULTS: The analysis revealed a heterozygous c.844ins68 variant in CBS, which was carried by all affected individuals and inherited from their mother. All affected individuals had a variable set of additional low frequency deleterious variants in PTK7, PLCD4, IL4I1 or RASSF4 as likely causal loci contributing to the disease development. CONCLUSION: This report extends the current knowledge of the genetic background of NTDs and proposes that common and low frequency variants in genes involved mostly in one-carbon metabolism or planar cell polarity (PCP) pathways can act in an additive manner to increase the genetic risk of the disease.

Whole-exome sequencing identifies de novo mutation in the COL1A1 gene to underlie the severe osteogenesis imperfecta.

BACKGROUND: Osteogenesis imperfecta (OI) comprises a clinically and genetically heterogeneous group of connective tissue disorders, characterized by low bone mass, increased bone fragility, and blue-gray eye sclera. OI often results from missense mutations in one of the conserved glycine residues present in the Gly-X-Y sequence repeats of the triple helical region of the collagen type I α chain, which is encoded by the COL1A1 gene. The aim of the present study is to describe the phenotype of OI II patient and a novel mutation, causing current phenotype. RESULTS: We report an undescribed de novo COL1A1 mutation in a patient affected by severe OI. After performing the whole-exome sequencing in a case parent-child trio, we identified a novel heterozygous c.2317G > T missense mutation in the COL1A1 gene, which leads to p.Gly773Cys transversion in the triple helical domain of the collagen type I α chain. The presence of the missense mutation was confirmed with the Sanger sequencing. CONCLUSIONS: Hereby, we report a novel mutation in the COL1A1 gene causing severe, life threatening OI and indicate the role of de novo mutation in the pathogenesis of rare familial diseases. Our study underlines the importance of exome sequencing in disease gene discovery for families where conventional genetic testing does not give conclusive evidence.

Structured education can improve primary-care management of headache: the first empirical evidence, from a controlled interventional study.

BACKGROUND: Headache disorders are under-recognized and under-diagnosed. A principal factor in their suboptimal management is lack of headache-related training among health-care providers, especially in primary care. In Estonia, general practitioners (GPs) refer many headache patients to neurological specialist services, mostly unnecessarily. GPs request "diagnostic" investigations, which are usually unhelpful and therefore wasteful. GP-made headache diagnoses are often arcane and non-specific, and treatments based on these are inappropriate. The aim of this study was to develop, implement and test an educational model intended to improve headache-related primary health care in Estonia. METHODS: This was a controlled study consisting of baseline observation, intervention and follow-up observation using the same measures of effect. It involved six GPs in Põlva and the surrounding region in Southern Estonia, together with their future patients presenting consecutively with headache as their main complaint, all with their consent. The primary outcome measure was referral rate (RR) to neurological specialist services. Secondary measures included number of GP-requested investigations, GP-made headache diagnoses and how these conformed to standard terminology (ICD-10), and GP-recommended or initiated treatments. RESULTS: RR at baseline (n = 490) was 39.5 %, falling to 34.7 % in the post-intervention group (n = 295) (overall reduction 4.8 %; p = 0.21). In the large subgroup of patients (88 %) for whom GPs made clearly headache-related ICD-10 diagnoses, RR fell by one fifth (from 40 to 32 %; p = 0.08), but the only diagnosis-related RR that showed a statistically significant reduction was (pericranial) myalgia (19 to 3 %; p = 0.03). There was a significant increase towards use of more specific diagnoses. Use of investigations in diagnosing headache reduced from 26 to 4 % (p < 0.0001). Initiation of treatment by GPs increased from 58 to 81 % (p < 0.0001). CONCLUSIONS: These were modest changes in GPs' entrenched behaviour. Nevertheless they were empirical evidence that GPs' practice in the field of headache could be improved by structured education. Furthermore, the changes were likely to be cost-saving. To our knowledge this study is the first to produce such evidence.

CDH13 promoter SNPs with pleiotropic effect on cardiometabolic parameters represent methylation QTLs.

CDH13 encodes T-cadherin, a receptor for high molecular weight (HMW) adiponectin and low-density lipoprotein, promoting proliferation and migration of endothelial cells. Genome-wide association studies have mapped multiple variants in CDH13 associated with cardiometabolic traits (CMT) with variable effects across studies. We hypothesized that this heterogeneity might reflect interplay with DNA methylation within the region. Resequencing and EpiTYPER™ assay were applied for the HYPertension in ESTonia/Coronary Artery Disease in Czech (HYPEST/CADCZ; n = 358) samples to identify CDH13 promoter SNPs acting as methylation Quantitative Trait Loci (meQTLs) and to investigate their associations with CMT. In silico data were extracted from genome-wide DNA methylation and genotype datasets of the population-based sample Estonian Genome Center of the University of Tartu (EGCUT; n = 165). HYPEST-CADCZ meta-analysis identified a rare variant rs113460564 as highly significant meQTL for a 134-bp distant CpG site (P = 5.90 × 10(-6); ß = 3.19%). Four common SNPs (rs12443878, rs12444338, rs62040565, rs8060301) exhibited effect on methylation level of up to 3 neighboring CpG sites in both datasets. The strongest association was detected in EGCUT between rs8060301 and cg09415485 (false discovery rate corrected P value = 1.89 × 10(-30)). Simultaneously, rs8060301 showed association with diastolic blood pressure, serum high-density lipoprotein and HMW adiponectin (P < 0.005). Novel strong associations were identified between rare CDH13 promoter meQTLs (minor allele frequency <5%) and HMW adiponectin: rs2239857 (P = 5.50 × 10(-5), ß = -1,841.9 ng/mL) and rs77068073 (P = 2.67 × 10(-4), ß = -2,484.4 ng/mL). Our study shows conclusively that CDH13 promoter harbors meQTLs associated with CMTs. It paves the way to deeper understanding of the interplay between DNA variation and methylation in susceptibility to common diseases.

In-solution hybrid capture of bisulfite-converted DNA for targeted bisulfite sequencing of 174 ADME genes.

DNA methylation is one of the most important epigenetic alterations involved in the control of gene expression. Bisulfite sequencing of genomic DNA is currently the only method to study DNA methylation patterns at single-nucleotide resolution. Hence, next-generation sequencing of bisulfite-converted DNA is the method of choice to investigate DNA methylation profiles at the genome-wide scale. Nevertheless, whole genome sequencing for analysis of human methylomes is expensive, and a method for targeted gene analysis would provide a good alternative in many cases where the primary interest is restricted to a set of genes. Here, we report the successful use of a custom Agilent SureSelect Target Enrichment system for the hybrid capture of bisulfite-converted DNA. We prepared bisulfite-converted next-generation sequencing libraries, which are enriched for the coding and regulatory regions of 174 ADME genes (i.e. genes involved in the metabolism and distribution of drugs). Sequencing of these libraries on Illumina's HiSeq2000 revealed that the method allows a reliable quantification of methylation levels of CpG sites in the selected genes, and validation of the method using pyrosequencing and the Illumina 450K methylation BeadChips revealed good concordance.

Evidence of inbreeding depression on human height.

Stature is a classical and highly heritable complex trait, with 80%-90% of variation explained by genetic factors. In recent years, genome-wide association studies (GWAS) have successfully identified many common additive variants influencing human height; however, little attention has been given to the potential role of recessive genetic effects. Here, we investigated genome-wide recessive effects by an analysis of inbreeding depression on adult height in over 35,000 people from 21 different population samples. We found a highly significant inverse association between height and genome-wide homozygosity, equivalent to a height reduction of up to 3 cm in the offspring of first cousins compared with the offspring of unrelated individuals, an effect which remained after controlling for the effects of socio-economic status, an important confounder (χ(2) = 83.89, df = 1; p = 5.2 × 10(-20)). There was, however, a high degree of heterogeneity among populations: whereas the direction of the effect was consistent across most population samples, the effect size differed significantly among populations. It is likely that this reflects true biological heterogeneity: whether or not an effect can be observed will depend on both the variance in homozygosity in the population and the chance inheritance of individual recessive genotypes. These results predict that multiple, rare, recessive variants influence human height. Although this exploratory work focuses on height alone, the methodology developed is generally applicable to heritable quantitative traits (QT), paving the way for an investigation into inbreeding effects, and therefore genetic architecture, on a range of QT of biomedical importance.

Genetic and epigenetic regulation of gene expression in fetal and adult human livers.

BACKGROUND: The liver plays a central role in the maintenance of homeostasis and health in general. However, there is substantial inter-individual variation in hepatic gene expression, and although numerous genetic factors have been identified, less is known about the epigenetic factors. RESULTS: By analyzing the methylomes and transcriptomes of 14 fetal and 181 adult livers, we identified 657 differentially methylated genes with adult-specific expression, these genes were enriched for transcription factor binding sites of HNF1A and HNF4A. We also identified 1,000 genes specific to fetal liver, which were enriched for GATA1, STAT5A, STAT5B and YY1 binding sites. We saw strong liver-specific effects of single nucleotide polymorphisms on both methylation levels (28,447 unique CpG sites (meQTL)) and gene expression levels (526 unique genes (eQTL)), at a false discovery rate (FDR) < 0.05. Of the 526 unique eQTL associated genes, 293 correlated significantly not only with genetic variation but also with methylation levels. The tissue-specificities of these associations were analyzed in muscle, subcutaneous adipose tissue and visceral adipose tissue. We observed that meQTL were more stable between tissues than eQTL and a very strong tissue-specificity for the identified associations between CpG methylation and gene expression. CONCLUSIONS: Our analyses generated a comprehensive resource of factors involved in the regulation of hepatic gene expression, and allowed us to estimate the proportion of variation in gene expression that could be attributed to genetic and epigenetic variation, both crucial to understanding differences in drug response and the etiology of liver diseases.

Preoperative corticosteroid injections are associated with worse long-term outcome of surgical carpal tunnel release.

BACKGROUND AND PURPOSE: Failed closed treatment of carpal tunnel syndrome (CTS) is often followed by surgery. We investigated whether preoperative steroid injections could have a negative effect on the long-term outcome of the operation. PATIENTS AND METHODS: 174 hands (164 patients) were operated on by a single surgeon at Tartu University Hospital in 2005. The patients were interviewed by telephone 5-6 years after the operation. Self-reported data were gathered retrospectively concerning the number of steroid injections received before the surgery and the perceived regression of symptoms (on a 100-point numeric rating scale) at the time of interview. The patients were also asked about the presence of specific symptoms of CTS if regression of their symptoms had not been complete. RESULTS: 93 of the 174 hands had complete regression of symptoms. Each additional injection was associated with an increased risk of occurrence of pain (RR = 1.1, 95% CI: 1.02-1-2), paresthesiae (RR = 1.1, CI: 1.1-1.2), and nocturnal awakenings (RR = 1.2, CI: 1.1-1.3). There was a weak association between the number of injections and the score given to regression of symptoms. INTERPRETATION: This is the first study to indicate that patients who received a greater number of local steroid injections preoperatively were more likely to have postoperative complaints associated with CTS.

Genetic determinants of plasma protein levels in the Estonian population.

The proteome holds great potential as an intermediate layer between the genome and phenome. Previous protein quantitative trait locus studies have focused mainly on describing the effects of common genetic variations on the proteome. Here, we assessed the impact of the common and rare genetic variations as well as the copy number variants (CNVs) on 326 plasma proteins measured in up to 500 individuals. We identified 184 cis and 94 trans signals for 157 protein traits, which were further fine-mapped to credible sets for 101 cis and 87 trans signals for 151 proteins. Rare genetic variation contributed to the levels of 7 proteins, with 5 cis and 14 trans associations. CNVs were associated with the levels of 11 proteins (7 cis and 5 trans), examples including a 3q12.1 deletion acting as a hub for multiple trans associations; and a CNV overlapping NAIP, a sensor component of the NAIP-NLRC4 inflammasome which is affecting pro-inflammatory cytokine interleukin 18 levels. In summary, this work presents a comprehensive resource of genetic variation affecting the plasma protein levels and provides the interpretation of identified effects.

A Large-Scale Genome-Wide Study of Gene-Sleep Duration Interactions for Blood Pressure in 811,405 Individuals from Diverse Populations.

Although both short and long sleep duration are associated with elevated hypertension risk, our understanding of their interplay with biological pathways governing blood pressure remains limited. To address this, we carried out genome-wide cross-population gene-by-short-sleep and long-sleep duration interaction analyses for three blood pressure traits (systolic, diastolic, and pulse pressure) in 811,405 individuals from diverse population groups. We discover 22 novel gene-sleep duration interaction loci for blood pressure, mapped to genes involved in neurological, thyroidal, bone metabolism, and hematopoietic pathways. Non-overlap between short sleep (12) and long sleep (10) interactions underscores the plausibility of distinct influences of both sleep duration extremes in cardiovascular health. With several of our loci reflecting specificity towards population background or sex, our discovery sheds light on the importance of embracing granularity when addressing heterogeneity entangled in gene-environment interactions, and in therapeutic design approaches for blood pressure management.

The role of adiposity in cardiometabolic traits: a Mendelian randomization analysis.

BACKGROUND: The association between adiposity and cardiometabolic traits is well known from epidemiological studies. Whilst the causal relationship is clear for some of these traits, for others it is not. We aimed to determine whether adiposity is causally related to various cardiometabolic traits using the Mendelian randomization approach. METHODS AND FINDINGS: We used the adiposity-associated variant rs9939609 at the FTO locus as an instrumental variable (IV) for body mass index (BMI) in a Mendelian randomization design. Thirty-six population-based studies of individuals of European descent contributed to the analyses. Age- and sex-adjusted regression models were fitted to test for association between (i) rs9939609 and BMI (n  =  198,502), (ii) rs9939609 and 24 traits, and (iii) BMI and 24 traits. The causal effect of BMI on the outcome measures was quantified by IV estimators. The estimators were compared to the BMI-trait associations derived from the same individuals. In the IV analysis, we demonstrated novel evidence for a causal relationship between adiposity and incident heart failure (hazard ratio, 1.19 per BMI-unit increase; 95% CI, 1.03-1.39) and replicated earlier reports of a causal association with type 2 diabetes, metabolic syndrome, dyslipidemia, and hypertension (odds ratio for IV estimator, 1.1-1.4; all p < 0.05). For quantitative traits, our results provide novel evidence for a causal effect of adiposity on the liver enzymes alanine aminotransferase and gamma-glutamyl transferase and confirm previous reports of a causal effect of adiposity on systolic and diastolic blood pressure, fasting insulin, 2-h post-load glucose from the oral glucose tolerance test, C-reactive protein, triglycerides, and high-density lipoprotein cholesterol levels (all p < 0.05). The estimated causal effects were in agreement with traditional observational measures in all instances except for type 2 diabetes, where the causal estimate was larger than the observational estimate (p  =  0.001). CONCLUSIONS: We provide novel evidence for a causal relationship between adiposity and heart failure as well as between adiposity and increased liver enzymes.