Longitudinal metabolomics of increasing body-mass index and waist-hip ratio reveals two dynamic patterns of obesity pandemic.

BACKGROUND/OBJECTIVE: This observational study dissects the complex temporal associations between body-mass index (BMI), waist-hip ratio (WHR) and circulating metabolomics using a combination of longitudinal and cross-sectional population-based datasets and new systems epidemiology tools. SUBJECTS/METHODS: Firstly, a data-driven subgrouping algorithm was employed to simplify high-dimensional metabolic profiling data into a single categorical variable: a self-organizing map (SOM) was created from 174 metabolic measures from cross-sectional surveys (FINRISK, n = 9708, ages 25-74) and a birth cohort (NFBC1966, n = 3117, age 31 at baseline, age 46 at follow-up) and an expert committee defined four subgroups of individuals based on visual inspection of the SOM. Secondly, the subgroups were compared regarding BMI and WHR trajectories in an independent longitudinal dataset: participants of the Young Finns Study (YFS, n = 1286, ages 24-39 at baseline, 10 years follow-up, three visits) were categorized into the four subgroups and subgroup-specific age-dependent trajectories of BMI, WHR and metabolic measures were modelled by linear regression. RESULTS: The four subgroups were characterised at age 39 by high BMI, WHR and dyslipidemia (designated TG-rich); low BMI, WHR and favourable lipids (TG-poor); low lipids in general (Low lipid) and high low-density-lipoprotein cholesterol (High LDL-C). Trajectory modelling of the YFS dataset revealed a dynamic BMI divergence pattern: despite overlapping starting points at age 24, the subgroups diverged in BMI, fasting insulin (three-fold difference at age 49 between TG-rich and TG-poor) and insulin-associated measures such as triglyceride-cholesterol ratio. Trajectories also revealed a WHR progression pattern: despite different starting points at the age of 24 in WHR, LDL-C and cholesterol-associated measures, all subgroups exhibited similar rates of change in these measures, i.e. WHR progression was uniform regardless of the cross-sectional metabolic profile. CONCLUSIONS: Age-associated weight variation in adults between 24 and 49 manifests as temporal divergence in BMI and uniform progression of WHR across metabolic health strata.

Metabolic profile-based subgroups can identify differences in brain volumes and brain iron deposition.

AIMS: To evaluate associations of metabolic profiles and biomarkers with brain atrophy, lesions, and iron deposition to understand the early risk factors associated with dementia. MATERIALS AND METHODS: Using data from 26 239 UK Biobank participants free from dementia and stroke, we assessed the associations of metabolic subgroups, derived using an artificial neural network approach (self-organizing map), and 39 individual biomarkers with brain MRI measures: total brain volume (TBV), grey matter volume (GMV), white matter volume (WMV), hippocampal volume (HV), white matter hyperintensity (WMH) volume, and caudate iron deposition. RESULTS: In metabolic subgroup analyses, participants characterized by high triglycerides and liver enzymes showed the most adverse brain outcomes compared to the healthy reference subgroup with high-density lipoprotein cholesterol and low body mass index (BMI) including associations with GMV (ßstandardized -0.20, 95% confidence interval [CI] -0.24 to -0.16), HV (ßstandardized -0.09, 95% CI -0.13 to -0.04), WMH volume (ßstandardized 0.22, 95% CI 0.18 to 0.26), and caudate iron deposition (ßstandardized 0.30, 95% CI 0.25 to 0.34), with similar adverse associations for the subgroup with high BMI, C-reactive protein and cystatin C, and the subgroup with high blood pressure (BP) and apolipoprotein B. Among the biomarkers, striking associations were seen between basal metabolic rate (BMR) and caudate iron deposition (ßstandardized 0.23, 95% CI 0.22 to 0.24 per 1 SD increase), GMV (ßstandardized -0.15, 95% CI -0.16 to -0.14) and HV (ßstandardized -0.11, 95% CI -0.12 to -0.10), and between BP and WMH volume (ßstandardized 0.13, 95% CI 0.12 to 0.14 for diastolic BP). CONCLUSIONS: Metabolic profiles were associated differentially with brain neuroimaging characteristics. Associations of BMR, BP and other individual biomarkers may provide insights into actionable mechanisms driving these brain associations.

Multi-Cohort Transcriptomic Subtyping of B-Cell Acute Lymphoblastic Leukemia.

RNA sequencing provides a snapshot of the functional consequences of genomic lesions that drive acute lymphoblastic leukemia (ALL). The aims of this study were to elucidate diagnostic associations (via machine learning) between mRNA-seq profiles, independently verify ALL lesions and develop easy-to-interpret transcriptome-wide biomarkers for ALL subtyping in the clinical setting. A training dataset of 1279 ALL patients from six North American cohorts was used for developing machine learning models. Results were validated in 767 patients from Australia with a quality control dataset across 31 tissues from 1160 non-ALL donors. A novel batch correction method was introduced and applied to adjust for cohort differences. Out of 18,503 genes with usable expression, 11,830 (64%) were confounded by cohort effects and excluded. Six ALL subtypes (ETV6::RUNX1, KMT2A, DUX4, PAX5 P80R, TCF3::PBX1, ZNF384) that covered 32% of patients were robustly detected by mRNA-seq (positive predictive value ≥ 87%). Five other frequent subtypes (CRLF2, hypodiploid, hyperdiploid, PAX5 alterations and Ph-positive) were distinguishable in 40% of patients at lower accuracy (52% ≤ positive predictive value ≤ 73%). Based on these findings, we introduce the Allspice R package to predict ALL subtypes and driver genes from unadjusted mRNA-seq read counts as encountered in real-world settings. Two examples of Allspice applied to previously unseen ALL patient samples with atypical lesions are included.

Gene networks and pathways for plasma lipid traits via multitissue multiomics systems analysis.

Genome-wide association studies (GWASs) have implicated ∼380 genetic loci for plasma lipid regulation. However, these loci only explain 17-27% of the trait variance, and a comprehensive understanding of the molecular mechanisms has not been achieved. In this study, we utilized an integrative genomics approach leveraging diverse genomic data from human populations to investigate whether genetic variants associated with various plasma lipid traits, namely, total cholesterol, high and low density lipoprotein cholesterol (HDL and LDL), and triglycerides, from GWASs were concentrated on specific parts of tissue-specific gene regulatory networks. In addition to the expected lipid metabolism pathways, gene subnetworks involved in "interferon signaling," "autoimmune/immune activation," "visual transduction," and "protein catabolism" were significantly associated with all lipid traits. In addition, we detected trait-specific subnetworks, including cadherin-associated subnetworks for LDL; glutathione metabolism for HDL; valine, leucine, and isoleucine biosynthesis for total cholesterol; and insulin signaling and complement pathways for triglyceride. Finally, by using gene-gene relations revealed by tissue-specific gene regulatory networks, we detected both known (e.g., APOH, APOA4, and ABCA1) and novel (e.g., F2 in adipose tissue) key regulator genes in these lipid-associated subnetworks. Knockdown of the F2 gene (coagulation factor II, thrombin) in 3T3-L1 and C3H10T1/2 adipocytes altered gene expression of Abcb11, Apoa5, Apof, Fabp1, Lipc, and Cd36; reduced intracellular adipocyte lipid content; and increased extracellular lipid content, supporting a link between adipose thrombin and lipid regulation. Our results shed light on the complex mechanisms underlying lipid metabolism and highlight potential novel targets for lipid regulation and lipid-associated diseases.

Adiposity and cancer: a Mendelian randomization analysis in the UK biobank.

BACKGROUND: Observational and Mendelian randomization (MR) studies link obesity and cancer, but it remains unclear whether these depend upon related metabolic abnormalities. METHODS: We used information from 321,472 participants in the UK biobank, including 30,561 cases of obesity-related cancer. We constructed three genetic instruments reflecting higher adiposity together with either "unfavourable" (82 SNPs), "favourable" (24 SNPs) or "neutral" metabolic profile (25 SNPs). We looked at associations with 14 types of cancer, previously suggested to be associated with obesity. RESULTS: All genetic instruments had a strong association with BMI (p < 1 × 10-300 for all). The instrument reflecting unfavourable adiposity was also associated with higher CRP, HbA1c and adverse lipid profile, while instrument reflecting metabolically favourable adiposity was associated with lower HbA1c and a favourable lipid profile. In MR-inverse-variance weighted analysis unfavourable adiposity was associated with an increased risk of non-hormonal cancers (OR = 1.22, 95% confidence interval [CI]:1.08, 1.38), but a lower risk of hormonal cancers (OR = 0.80, 95%CI: 0.72, 0.89). From individual cancers, MR analyses suggested causal increases in the risk of multiple myeloma (OR = 1.36, 95%CI: 1.09, 1.70) and endometrial cancer (OR = 1.77, 95%CI: 1.16, 2.68) by greater genetically instrumented unfavourable adiposity but lower risks of breast and prostate cancer (OR = 0.72, 95%CI: 0.61, 0.83 and OR = 0.81, 95%CI: 0.68, 0.97, respectively). Favourable or neutral adiposity were not associated with the odds of any individual cancer. CONCLUSIONS: Higher adiposity associated with a higher risk of non-hormonal cancer but a lower risk of some hormone related cancers. Presence of metabolic abnormalities might aggravate the adverse effects of higher adiposity on cancer. Further studies are warranted to investigate whether interventions on adverse metabolic health may help to alleviate obesity-related cancer risk.

Insulin resistance and systemic metabolic changes in oral glucose tolerance test in 5340 individuals: an interventional study.

BACKGROUND: Insulin resistance (IR) is predictive for type 2 diabetes and associated with various metabolic abnormalities in fasting conditions. However, limited data are available on how IR affects metabolic responses in a non-fasting setting, yet this is the state people are mostly exposed to during waking hours in the modern society. Here, we aim to comprehensively characterise the metabolic changes in response to an oral glucose test (OGTT) and assess the associations of these changes with IR. METHODS: Blood samples were obtained at 0 (fasting baseline, right before glucose ingestion), 30, 60, and 120 min during the OGTT. Seventy-eight metabolic measures were analysed at each time point for a discovery cohort of 4745 middle-aged Finnish individuals and a replication cohort of 595 senior Finnish participants. We assessed the metabolic changes in response to glucose ingestion (percentage change in relative to fasting baseline) across the four time points and further compared the response profile between five groups with different levels of IR and glucose intolerance. Further, the differences were tested for covariate adjustment, including gender, body mass index, systolic blood pressure, fasting, and 2-h glucose levels. The groups were defined as insulin sensitive with normal glucose (IS-NGT), insulin resistant with normal glucose (IR-NGT), impaired fasting glucose (IFG), impaired glucose tolerance (IGT), and new diabetes (NDM). IS-NGT and IR-NGT were defined as the first and fourth quartile of fasting insulin in NGT individuals. RESULTS: Glucose ingestion induced multiple metabolic responses, including increased glycolysis intermediates and decreased branched-chain amino acids, ketone bodies, glycerol, and triglycerides. The IR-NGT subgroup showed smaller responses for these measures (mean + 23%, interquartile 9-34% at 120 min) compared to IS-NGT (34%, 23-44%, P < 0.0006 for difference, corrected for multiple testing). Notably, the three groups with glucose abnormality (IFG, IGT, and NDM) showed similar metabolic dysregulations as those of IR-NGT. The difference between the IS-NGT and the other subgroups was largely explained by fasting insulin, but not fasting or 2 h glucose. The findings were consistent after covariate adjustment and between the discovery and replication cohort. CONCLUSIONS: Insulin-resistant non-diabetic individuals are exposed to a similar adverse postprandial metabolic milieu, and analogous cardiometabolic risk, as those with type 2 diabetes. The wide range of metabolic abnormalities associated with IR highlights the necessity of diabetes diagnostics and clinical care beyond glucose management.

Statistical reporting of metabolomics data: experience from a high-throughput NMR platform and epidemiological applications.

INTRODUCTION: Meta-analysis is the cornerstone of robust biomedical evidence. OBJECTIVES: We investigated whether statistical reporting practices facilitate metabolomics meta-analyses. METHODS: A literature review of 44 studies that used a comparable platform. RESULTS: Non-numeric formats were used in 31 studies. In half of the studies, less than a third of all measures were reported. Unadjusted P-values were missing from 12 studies and exact P-values from 9 studies. CONCLUSION: Reporting practices can be improved. We recommend (i) publishing all results as numbers, (ii) reporting effect sizes of all measured metabolites and (iii) always reporting unadjusted exact P-values.

Genetic variation within endolysosomal system is associated with late-onset Alzheimer's disease.

Late-onset Alzheimer's disease is the most common dementia type, yet no treatment exists to stop the neurodegeneration. Evidence from monogenic lysosomal diseases, neuronal pathology and experimental models suggest that autophagic and endolysosomal dysfunction may contribute to neurodegeneration by disrupting the degradation of potentially neurotoxic molecules such as amyloid-ß and tau. However, it is uncertain how well the evidence from rare disorders and experimental models capture causal processes in common forms of dementia, including late-onset Alzheimer's disease. For this reason, we set out to investigate if autophagic and endolysosomal genes were enriched for genetic variants that convey increased risk of Alzheimer's disease; such a finding would provide population-based support for the endolysosomal hypothesis of neurodegeneration. We quantified the collective genetic associations between the endolysosomal system and Alzheimer's disease in three genome-wide associations studies (combined n = 62 415). We used the Mergeomics pathway enrichment algorithm that incorporates permutations of the full hierarchical cascade of SNP-gene-pathway to estimate enrichment. We used a previously published collection of 891 autophagic and endolysosomal genes (denoted as AphagEndoLyso, and derived from the Lysoplex sequencing platform) as a proxy for cellular processes related to autophagy, endocytosis and lysosomal function. We also investigated a subset of 142 genes of the 891 that have been implicated in Mendelian diseases (MenDisLyso). We found that both gene sets were enriched for genetic Alzheimer's associations: an enrichment score 3.67 standard deviations from the null model (P = 0.00012) was detected for AphagEndoLyso, and a score 3.36 standard deviations from the null model (P = 0.00039) was detected for MenDisLyso. The high enrichment score was specific to the AphagEndoLyso gene set (stronger than 99.7% of other tested pathways) and to Alzheimer's disease (stronger than all other tested diseases). The APOE locus explained most of the MenDisLyso signal (1.16 standard deviations after APOE removal, P = 0.12), but the AphagEndoLyso signal was less affected (3.35 standard deviations after APOE removal, P = 0.00040). Additional sensitivity analyses further indicated that the AphagEndoLyso Gene Set contained an aggregate genetic association that comprised a combination of subtle genetic signals in multiple genes. We also observed an enrichment of Parkinson's disease signals for MenDisLyso (3.25 standard deviations) and for AphagEndoLyso (3.95 standard deviations from the null model), and a brain-specific pattern of gene expression for AphagEndoLyso in the Gene Tissue Expression Project dataset. These results provide evidence that a diffuse aggregation of genetic perturbations to the autophagy and endolysosomal system may mediate late-onset Alzheimer's risk in human populations.

Data-driven metabolic subtypes predict future adverse events in individuals with type 1 diabetes.

AIMS/HYPOTHESIS: Previously, we proposed that data-driven metabolic subtypes predict mortality in type 1 diabetes. Here, we analysed new clinical endpoints and revisited the subtypes after 7 years of additional follow-up. METHODS: Finnish individuals with type 1 diabetes (2059 men and 1924 women, insulin treatment before 35 years of age) were recruited by the national multicentre FinnDiane Study Group. The participants were assigned one of six metabolic subtypes according to a previously published self-organising map from 2008. Subtype-specific all-cause and cardiovascular mortality rates in the FinnDiane cohort were compared with registry data from the entire Finnish population. The rates of incident diabetic kidney disease and cardiovascular endpoints were estimated based on hospital records. RESULTS: The advanced kidney disease subtype was associated with the highest incidence of kidney disease progression (67.5% per decade, p < 0.001), ischaemic heart disease (26.4% per decade, p < 0.001) and all-cause mortality (41.5% per decade, p < 0.001). Across all subtypes, mortality rates were lower in women compared with men, but standardised mortality ratios (SMRs) were higher in women. SMRs were indistinguishable between the original study period (1994-2007) and the new period (2008-2014). The metabolic syndrome subtype predicted cardiovascular deaths (SMR 11.0 for men, SMR 23.4 for women, p < 0.001), and women with the high HDL-cholesterol subtype were also at high cardiovascular risk (SMR 16.3, p < 0.001). Men with the low-cholesterol or good glycaemic control subtype showed no excess mortality. CONCLUSIONS/INTERPRETATION: Data-driven multivariable metabolic subtypes predicted the divergence of complication burden across multiple clinical endpoints simultaneously. In particular, men with the metabolic syndrome and women with high HDL-cholesterol should be recognised as important subgroups in interventional studies and public health guidelines on type 1 diabetes.

Integrative genomics reveals novel molecular pathways and gene networks for coronary artery disease.

The majority of the heritability of coronary artery disease (CAD) remains unexplained, despite recent successes of genome-wide association studies (GWAS) in identifying novel susceptibility loci. Integrating functional genomic data from a variety of sources with a large-scale meta-analysis of CAD GWAS may facilitate the identification of novel biological processes and genes involved in CAD, as well as clarify the causal relationships of established processes. Towards this end, we integrated 14 GWAS from the CARDIoGRAM Consortium and two additional GWAS from the Ottawa Heart Institute (25,491 cases and 66,819 controls) with 1) genetics of gene expression studies of CAD-relevant tissues in humans, 2) metabolic and signaling pathways from public databases, and 3) data-driven, tissue-specific gene networks from a multitude of human and mouse experiments. We not only detected CAD-associated gene networks of lipid metabolism, coagulation, immunity, and additional networks with no clear functional annotation, but also revealed key driver genes for each CAD network based on the topology of the gene regulatory networks. In particular, we found a gene network involved in antigen processing to be strongly associated with CAD. The key driver genes of this network included glyoxalase I (GLO1) and peptidylprolyl isomerase I (PPIL1), which we verified as regulatory by siRNA experiments in human aortic endothelial cells. Our results suggest genetic influences on a diverse set of both known and novel biological processes that contribute to CAD risk. The key driver genes for these networks highlight potential novel targets for further mechanistic studies and therapeutic interventions.

Mergeomics: a web server for identifying pathological pathways, networks, and key regulators via multidimensional data integration.

BACKGROUND: Human diseases are commonly the result of multidimensional changes at molecular, cellular, and systemic levels. Recent advances in genomic technologies have enabled an outpour of omics datasets that capture these changes. However, separate analyses of these various data only provide fragmented understanding and do not capture the holistic view of disease mechanisms. To meet the urgent needs for tools that effectively integrate multiple types of omics data to derive biological insights, we have developed Mergeomics, a computational pipeline that integrates multidimensional disease association data with functional genomics and molecular networks to retrieve biological pathways, gene networks, and central regulators critical for disease development. RESULTS: To make the Mergeomics pipeline available to a wider research community, we have implemented an online, user-friendly web server ( http://mergeomics. RESEARCH: idre.ucla.edu/ ). The web server features a modular implementation of the Mergeomics pipeline with detailed tutorials. Additionally, it provides curated genomic resources including tissue-specific expression quantitative trait loci, ENCODE functional annotations, biological pathways, and molecular networks, and offers interactive visualization of analytical results. Multiple computational tools including Marker Dependency Filtering (MDF), Marker Set Enrichment Analysis (MSEA), Meta-MSEA, and Weighted Key Driver Analysis (wKDA) can be used separately or in flexible combinations. User-defined summary-level genomic association datasets (e.g., genetic, transcriptomic, epigenomic) related to a particular disease or phenotype can be uploaded and computed real-time to yield biologically interpretable results, which can be viewed online and downloaded for later use. CONCLUSIONS: Our Mergeomics web server offers researchers flexible and user-friendly tools to facilitate integration of multidimensional data into holistic views of disease mechanisms in the form of tissue-specific key regulators, biological pathways, and gene networks.

Mergeomics: multidimensional data integration to identify pathogenic perturbations to biological systems.

BACKGROUND: Complex diseases are characterized by multiple subtle perturbations to biological processes. New omics platforms can detect these perturbations, but translating the diverse molecular and statistical information into testable mechanistic hypotheses is challenging. Therefore, we set out to create a public tool that integrates these data across multiple datasets, platforms, study designs and species in order to detect the most promising targets for further mechanistic studies. RESULTS: We developed Mergeomics, a computational pipeline consisting of independent modules that 1) leverage multi-omics association data to identify biological processes that are perturbed in disease, and 2) overlay the disease-associated processes onto molecular interaction networks to pinpoint hubs as potential key regulators. Unlike existing tools that are mostly dedicated to specific data type or settings, the Mergeomics pipeline accepts and integrates datasets across platforms, data types and species. We optimized and evaluated the performance of Mergeomics using simulation and multiple independent datasets, and benchmarked the results against alternative methods. We also demonstrate the versatility of Mergeomics in two case studies that include genome-wide, epigenome-wide and transcriptome-wide datasets from human and mouse studies of total cholesterol and fasting glucose. In both cases, the Mergeomics pipeline provided statistical and contextual evidence to prioritize further investigations in the wet lab. The software implementation of Mergeomics is freely available as a Bioconductor R package. CONCLUSION: Mergeomics is a flexible and robust computational pipeline for multidimensional data integration. It outperforms existing tools, and is easily applicable to datasets from different studies, species and omics data types for the study of complex traits.

Metabolic profiling of pregnancy: cross-sectional and longitudinal evidence.

BACKGROUND: Pregnancy triggers well-known alterations in maternal glucose and lipid balance but its overall effects on systemic metabolism remain incompletely understood. METHODS: Detailed molecular profiles (87 metabolic measures and 37 cytokines) were measured for up to 4260 women (24-49 years, 322 pregnant) from three population-based cohorts in Finland. Circulating molecular concentrations in pregnant women were compared to those in non-pregnant women. Metabolic profiles were also reassessed for 583 women 6 years later to uncover the longitudinal metabolic changes in response to change in the pregnancy status. RESULTS: Compared to non-pregnant women, all lipoprotein subclasses and lipids were markedly increased in pregnant women. The most pronounced differences were observed for the intermediate-density, low-density and high-density lipoprotein triglyceride concentrations. Large differences were also seen for many fatty acids and amino acids. Pregnant women also had higher concentrations of low-grade inflammatory marker glycoprotein acetyls, higher concentrations of interleukin-18 and lower concentrations of interleukin-12p70. The changes in metabolic concentrations for women who were not pregnant at baseline but pregnant 6 years later (or vice versa) matched (or were mirror-images of) the cross-sectional association pattern. Cross-sectional results were consistent across the three cohorts and similar longitudinal changes were seen for 653 women in 4-year and 497 women in 10-year follow-up. For multiple metabolic measures, the changes increased in magnitude across the three trimesters. CONCLUSIONS: Pregnancy initiates substantial metabolic and inflammatory changes in the mothers. Comprehensive characterisation of normal pregnancy is important for gaining understanding of the key nutrients for fetal growth and development. These findings also provide a valuable molecular reference in relation to studies of adverse pregnancy outcomes.

Systems Genetics Analysis of Genome-Wide Association Study Reveals Novel Associations Between Key Biological Processes and Coronary Artery Disease.

OBJECTIVE: Genome-wide association studies have identified multiple genetic variants affecting the risk of coronary artery disease (CAD). However, individually these explain only a small fraction of the heritability of CAD and for most, the causal biological mechanisms remain unclear. We sought to obtain further insights into potential causal processes of CAD by integrating large-scale GWA data with expertly curated databases of core human pathways and functional networks. APPROACHES AND RESULTS: Using pathways (gene sets) from Reactome, we carried out a 2-stage gene set enrichment analysis strategy. From a meta-analyzed discovery cohort of 7 CAD genome-wide association study data sets (9889 cases/11 089 controls), nominally significant gene sets were tested for replication in a meta-analysis of 9 additional studies (15 502 cases/55 730 controls) from the Coronary ARtery DIsease Genome wide Replication and Meta-analysis (CARDIoGRAM) Consortium. A total of 32 of 639 Reactome pathways tested showed convincing association with CAD (replication P<0.05). These pathways resided in 9 of 21 core biological processes represented in Reactome, and included pathways relevant to extracellular matrix (ECM) integrity, innate immunity, axon guidance, and signaling by PDRF (platelet-derived growth factor), NOTCH, and the transforming growth factor-ß/SMAD receptor complex. Many of these pathways had strengths of association comparable to those observed in lipid transport pathways. Network analysis of unique genes within the replicated pathways further revealed several interconnected functional and topologically interacting modules representing novel associations (eg, semaphoring-regulated axonal guidance pathway) besides confirming known processes (lipid metabolism). The connectivity in the observed networks was statistically significant compared with random networks (P<0.001). Network centrality analysis (degree and betweenness) further identified genes (eg, NCAM1, FYN, FURIN, etc) likely to play critical roles in the maintenance and functioning of several of the replicated pathways. CONCLUSIONS: These findings provide novel insights into how genetic variation, interpreted in the context of biological processes and functional interactions among genes, may help define the genetic architecture of CAD.

New susceptibility loci associated with kidney disease in type 1 diabetes.

Diabetic kidney disease, or diabetic nephropathy (DN), is a major complication of diabetes and the leading cause of end-stage renal disease (ESRD) that requires dialysis treatment or kidney transplantation. In addition to the decrease in the quality of life, DN accounts for a large proportion of the excess mortality associated with type 1 diabetes (T1D). Whereas the degree of glycemia plays a pivotal role in DN, a subset of individuals with poorly controlled T1D do not develop DN. Furthermore, strong familial aggregation supports genetic susceptibility to DN. However, the genes and the molecular mechanisms behind the disease remain poorly understood, and current therapeutic strategies rarely result in reversal of DN. In the GEnetics of Nephropathy: an International Effort (GENIE) consortium, we have undertaken a meta-analysis of genome-wide association studies (GWAS) of T1D DN comprising ~2.4 million single nucleotide polymorphisms (SNPs) imputed in 6,691 individuals. After additional genotyping of 41 top ranked SNPs representing 24 independent signals in 5,873 individuals, combined meta-analysis revealed association of two SNPs with ESRD: rs7583877 in the AFF3 gene (P = 1.2 × 10(-8)) and an intergenic SNP on chromosome 15q26 between the genes RGMA and MCTP2, rs12437854 (P = 2.0 × 10(-9)). Functional data suggest that AFF3 influences renal tubule fibrosis via the transforming growth factor-beta (TGF-ß1) pathway. The strongest association with DN as a primary phenotype was seen for an intronic SNP in the ERBB4 gene (rs7588550, P = 2.1 × 10(-7)), a gene with type 2 diabetes DN differential expression and in the same intron as a variant with cis-eQTL expression of ERBB4. All these detected associations represent new signals in the pathogenesis of DN.

Novel genetic susceptibility loci for diabetic end-stage renal disease identified through robust naive Bayes classification.

AIMS/HYPOTHESIS: Diabetic nephropathy is a major diabetic complication, and diabetes is the leading cause of end-stage renal disease (ESRD). Family studies suggest a hereditary component for diabetic nephropathy. However, only a few genes have been associated with diabetic nephropathy or ESRD in diabetic patients. Our aim was to detect novel genetic variants associated with diabetic nephropathy and ESRD. METHODS: We exploited a novel algorithm, 'Bag of Naive Bayes', whose marker selection strategy is complementary to that of conventional genome-wide association models based on univariate association tests. The analysis was performed on a genome-wide association study of 3,464 patients with type 1 diabetes from the Finnish Diabetic Nephropathy (FinnDiane) Study and subsequently replicated with 4,263 type 1 diabetes patients from the Steno Diabetes Centre, the All Ireland-Warren 3-Genetics of Kidneys in Diabetes UK collection (UK-Republic of Ireland) and the Genetics of Kidneys in Diabetes US Study (GoKinD US). RESULTS: Five genetic loci (WNT4/ZBTB40-rs12137135, RGMA/MCTP2-rs17709344, MAPRE1P2-rs1670754, SEMA6D/SLC24A5-rs12917114 and SIK1-rs2838302) were associated with ESRD in the FinnDiane study. An association between ESRD and rs17709344, tagging the previously identified rs12437854 and located between the RGMA and MCTP2 genes, was replicated in independent case-control cohorts. rs12917114 near SEMA6D was associated with ESRD in the replication cohorts under the genotypic model (p < 0.05), and rs12137135 upstream of WNT4 was associated with ESRD in Steno. CONCLUSIONS/INTERPRETATION: This study supports the previously identified findings on the RGMA/MCTP2 region and suggests novel susceptibility loci for ESRD. This highlights the importance of applying complementary statistical methods to detect novel genetic variants in diabetic nephropathy and, in general, in complex diseases.

Genome-wide association study of urinary albumin excretion rate in patients with type 1 diabetes.

AIMS/HYPOTHESIS: An abnormal urinary albumin excretion rate (AER) is often the first clinically detectable manifestation of diabetic nephropathy. Our aim was to estimate the heritability and to detect genetic variation associated with elevated AER in patients with type 1 diabetes. METHODS: The discovery phase genome-wide association study (GWAS) included 1,925 patients with type 1 diabetes and with data on 24 h AER. AER was analysed as a continuous trait and the analysis was stratified by the use of antihypertensive medication. Signals with a p value <10(-4) were followed up in 3,750 additional patients with type 1 diabetes from seven studies. RESULTS: The narrow-sense heritability, captured with our genotyping platform, was estimated to explain 27.3% of the total AER variability, and 37.6% after adjustment for covariates. In the discovery stage, five single nucleotide polymorphisms in the GLRA3 gene were strongly associated with albuminuria (p < 5 × 10(-8)). In the replication group, a nominally significant association (p = 0.035) was observed between albuminuria and rs1564939 in GLRA3, but this was in the opposite direction. Sequencing of the surrounding genetic region in 48 Finnish and 48 UK individuals supported the possibility that population-specific rare variants contribute to the synthetic association observed at the common variants in GLRA3. The strongest replication (p = 0.026) was obtained for rs2410601 between the PSD3 and SH2D4A genes. Pathway analysis highlighted natural killer cell mediated immunity processes. CONCLUSIONS/INTERPRETATION: This study suggests novel pathways and molecular mechanisms for the pathogenesis of albuminuria in type 1 diabetes.

Long-term leisure-time physical activity and serum metabolome.

BACKGROUND: Long-term physical inactivity seems to cause many health problems. We studied whether persistent physical activity compared with inactivity has a global effect on serum metabolome toward reduced cardiometabolic disease risk. METHODS AND RESULTS: Sixteen same-sex twin pairs (mean age, 60 years) were selected from a cohort of twin pairs on the basis of their >30-year discordance for physical activity. Persistently (≥5 years) active and inactive groups in 3 population-based cohorts (mean ages, 31-52 years) were also studied (1037 age- and sex-matched pairs). Serum metabolome was quantified by nuclear magnetic resonance spectroscopy. We used permutation analysis to estimate the significance of the multivariate effect combined across all metabolic measures; univariate effects were estimated by paired testing in twins and in matched pairs in the cohorts, and by meta-analysis over all substudies. Persistent physical activity was associated with the multivariate metabolic profile in the twins (P=0.003), and a similar pattern was observed in all 3 population cohorts with differing mean ages. Isoleucine, α1-acid glycoprotein, and glucose were lower in the physically active than in the inactive individuals (P<0.001 in meta-analysis); serum fatty acid composition was shifted toward a less saturated profile; and lipoprotein subclasses were shifted toward lower very-low-density lipoprotein (P<0.001) and higher large and very large high-density lipoprotein (P<0.001) particle concentrations. The findings persisted after adjustment for body mass index. CONCLUSIONS: The numerous differences found between persistently physically active and inactive individuals in the circulating metabolome together indicate better metabolic health in the physically active than in inactive individuals.

Chromosome 2q31.1 associates with ESRD in women with type 1 diabetes.

Sex and genetic variation influence the risk of developing diabetic nephropathy and ESRD in patients with type 1 diabetes. We performed a genome-wide association study in a cohort of 3652 patients from the Finnish Diabetic Nephropathy (FinnDiane) Study with type 1 diabetes to determine whether sex-specific genetic risk factors for ESRD exist. A common variant, rs4972593 on chromosome 2q31.1, was associated with ESRD in women (P<5×10(-8)) but not in men (P=0.77). This association was replicated in the meta-analysis of three independent type 1 diabetes cohorts (P=0.02) and remained significant for women (P<5×10(-8); odds ratio, 1.81 [95% confidence interval, 1.47 to 2.24]) upon combined meta-analysis of the discovery and replication cohorts. rs4972593 is located between the genes that code for the Sp3 transcription factor, which interacts directly with estrogen receptor α and regulates the expression of genes linked to glomerular function and the pathogenesis of nephropathy, and the CDCA7 transcription factor, which regulates cell proliferation. Further examination revealed potential transcription factor-binding sites within rs4972593 and predicted eight estrogen-responsive elements within 5 kb of this locus. Moreover, we found sex-specific differences in the glomerular expression levels of SP3 (P=0.004). Overall, these results suggest that rs4972593 is a sex-specific genetic variant associated with ESRD in patients with type 1 diabetes and may underlie the sex-specific protection against ESRD.