A Federated Database for Obesity Research: An IMI-SOPHIA Study

Obesity is considered by many as a lifestyle choice rather than a chronic progressive disease. The Innovative Medicines Initiative (IMI) SOPHIA (Stratification of Obesity Phenotypes to Optimize Future Obesity Therapy) project is part of a momentum shift aiming to provide better tools for the stratification of people with obesity according to disease risk and treatment response. One of the challenges to achieving these goals is that many clinical cohorts are siloed, limiting the potential of combined data for biomarker discovery. In SOPHIA, we have addressed this challenge by setting up a federated database building on open-source DataSHIELD technology. The database currently federates 16 cohorts that are accessible via a central gateway. The database is multi-modal, including research studies, clinical trials, and routine health data, and is accessed using the R statistical programming environment where statistical and machine learning analyses can be performed at a distance without any disclosure of patient-level data. We demonstrate the use of the database by providing a proof-of-concept analysis, performing a federated linear model of BMI and systolic blood pressure, pooling all data from 16 studies virtually without any analyst seeing individual patient-level data. This analysis provided similar point estimates compared to a meta-analysis of the 16 individual studies. Our approach provides a benchmark for reproducible, safe federated analyses across multiple study types provided by multiple stakeholders.

Discovery of drug-omics associations in type 2 diabetes with generative deep-learning models

The application of multiple omics technologies in biomedical cohorts has the potential to reveal patient-level disease characteristics and individualized response to treatment. However, the scale and heterogeneous nature of multi-modal data makes integration and inference a non-trivial task. We developed a deep-learning-based framework, multi-omics variational autoencoders (MOVE), to integrate such data and applied it to a cohort of 789 people with newly diagnosed type 2 diabetes with deep multi-omics phenotyping from the DIRECT consortium. Using in silico perturbations, we identified drug-omics associations across the multi-modal datasets for the 20 most prevalent drugs given to people with type 2 diabetes with substantially higher sensitivity than univariate statistical tests. From these, we among others, identified novel associations between metformin and the gut microbiota as well as opposite molecular responses for the two statins, simvastatin and atorvastatin. We used the associations to quantify drug-drug similarities, assess the degree of polypharmacy and conclude that drug effects are distributed across the multi-omics modalities.

Investigating the causal relationships between excess adiposity and cardiometabolic health in men and women

AIMS/HYPOTHESIS

Excess adiposity is differentially associated with increased risk of cardiometabolic disease in men and women, according to observational studies. Causal inference studies largely assume a linear relationship between BMI and cardiometabolic outcomes, which may not be the case. In this study, we investigated the shapes of the causal relationships between BMI and cardiometabolic diseases and risk factors. We further investigated sex differences within the causal framework.

METHODS

To assess causal relationships between BMI and the outcomes, we used two-stage least-squares Mendelian randomisation (MR), with a polygenic risk score for BMI as the instrumental variable. To elucidate the shapes of the causal relationships, we used a non-linear MR fractional polynomial method, and used piecewise MR to investigate threshold relationships and confirm the shapes.

RESULTS

BMI was associated with type 2 diabetes (OR 3.10; 95% CI 2.73, 3.53), hypertension (OR 1.53; 95% CI 1.44, 1.62) and coronary artery disease (OR 1.20; 95% CI 1.08, 1.33), but not chronic kidney disease (OR 1.08; 95% CI 0.67, 1.72) or stroke (OR 1.08; 95% CI 0.92, 1.28). The data suggest that these relationships are non-linear. For cardiometabolic risk factors, BMI was positively associated with glucose, HbA1c, triacylglycerol levels and both systolic and diastolic BP. BMI had an inverse causal relationship with total cholesterol, LDL-cholesterol and HDL-cholesterol. The data suggest a non-linear causal relationship between BMI and BP and other biomarkers (p<0.001) except lipoprotein A. The piecewise MR results were consistent with the fractional polynomial results. The causal effect of BMI on coronary artery disease, total cholesterol and LDL-cholesterol was different in men and women, but this sex difference was only significant for LDL-cholesterol after controlling for multiple testing (p<0.001). Further, the causal effect of BMI on coronary artery disease varied by menopause status in women.

CONCLUSIONS/INTERPRETATION

We describe the shapes of causal effects of BMI on cardiometabolic diseases and risk factors, and report sex differences in the causal effects of BMI on LDL-cholesterol. We found evidence of non-linearity in the causal effect of BMI on diseases and risk factor biomarkers. Reducing excess adiposity is highly beneficial for health, but there is greater need to consider biological sex in the management of adiposity.

Four groups of type 2 diabetes contribute to the etiological and clinical heterogeneity in newly diagnosed individuals: An IMI DIRECT study

The presentation and underlying pathophysiology of type 2 diabetes (T2D) is complex and heterogeneous. Recent studies attempted to stratify T2D into distinct subgroups using data-driven approaches, but their clinical utility may be limited if categorical representations of complex phenotypes are suboptimal. We apply a soft-clustering (archetype) method to characterize newly diagnosed T2D based on 32 clinical variables. We assign quantitative clustering scores for individuals and investigate the associations with glycemic deterioration, genetic risk scores, circulating omics biomarkers, and phenotypic stability over 36 months. Four archetype profiles represent dysfunction patterns across combinations of T2D etiological processes and correlate with multiple circulating biomarkers. One archetype associated with obesity, insulin resistance, dyslipidemia, and impaired β cell glucose sensitivity corresponds with the fastest disease progression and highest demand for anti-diabetic treatment. We demonstrate that clinical heterogeneity in T2D can be mapped to heterogeneity in individual etiological processes, providing a potential route to personalized treatments.

Refining the accuracy of validated target identification through coding variant fine-mapping in type 2 diabetes

We aggregated coding variant data for 81,412 type 2 diabetes cases and 370,832 controls of diverse ancestry, identifying 40 coding variant association signals (P < 2.2 × 10-7); of these, 16 map outside known risk-associated loci. We make two important observations. First, only five of these signals are driven by low-frequency variants: even for these, effect sizes are modest (odds ratio ≤1.29). Second, when we used large-scale genome-wide association data to fine-map the associated variants in their regional context, accounting for the global enrichment of complex trait associations in coding sequence, compelling evidence for coding variant causality was obtained for only 16 signals. At 13 others, the associated coding variants clearly represent 'false leads' with potential to generate erroneous mechanistic inference. Coding variant associations offer a direct route to biological insight for complex diseases and identification of validated therapeutic targets; however, appropriate mechanistic inference requires careful specification of their causal contribution to disease predisposition.

Erratum: Sequence data and association statistics from 12,940 type 2 diabetes cases and controls

This corrects the article DOI: 10.1038/sdata.2017.179.

Sequence data and association statistics from 12,940 type 2 diabetes cases and controls

To investigate the genetic basis of type 2 diabetes (T2D) to high resolution, the GoT2D and T2D-GENES consortia catalogued variation from whole-genome sequencing of 2,657 European individuals and exome sequencing of 12,940 individuals of multiple ancestries. Over 27M SNPs, indels, and structural variants were identified, including 99% of low-frequency (minor allele frequency [MAF] 0.1-5%) non-coding variants in the whole-genome sequenced individuals and 99.7% of low-frequency coding variants in the whole-exome sequenced individuals. Each variant was tested for association with T2D in the sequenced individuals, and, to increase power, most were tested in larger numbers of individuals (>80% of low-frequency coding variants in ~82 K Europeans via the exome chip, and ~90% of low-frequency non-coding variants in ~44 K Europeans via genotype imputation). The variants, genotypes, and association statistics from these analyses provide the largest reference to date of human genetic information relevant to T2D, for use in activities such as T2D-focused genotype imputation, functional characterization of variants or genes, and other novel analyses to detect associations between sequence variation and T2D.