Plasma Metabolomics Identifies Markers of Impaired Renal Function: A Meta-analysis of 3089 Persons with Type 2 Diabetes

CONTEXT: There is a need for novel biomarkers and better understanding of the pathophysiology of diabetic kidney disease. OBJECTIVE: To investigate associations between plasma metabolites and kidney function in people with type 2 diabetes (T2D). DESIGN: 3089 samples from individuals with T2D, collected between 1999 and 2015, from 5 independent Dutch cohort studies were included. Up to 7 years follow-up was available in 1100 individuals from 2 of the cohorts. MAIN OUTCOME MEASURES: Plasma metabolites (n = 149) were measured by nuclear magnetic resonance spectroscopy. Associations between metabolites and estimated glomerular filtration rate (eGFR), urinary albumin-to-creatinine ratio (UACR), and eGFR slopes were investigated in each study followed by random effect meta-analysis. Adjustments included traditional cardiovascular risk factors and correction for multiple testing. RESULTS: In total, 125 metabolites were significantly associated (PFDR = 1.5×10-32 - 0.046; β = -11.98-2.17) with eGFR. Inverse associations with eGFR were demonstrated for branched-chain and aromatic amino acids (AAAs), glycoprotein acetyls, triglycerides (TGs), lipids in very low-density lipoproteins (VLDL) subclasses, and fatty acids (PFDR < 0.03). We observed positive associations with cholesterol and phospholipids in high-density lipoproteins (HDL) and apolipoprotein A1 (PFDR < 0.05). Albeit some metabolites were associated with UACR levels (P < 0.05), significance was lost after correction for multiple testing. Tyrosine and HDL-related metabolites were positively associated with eGFR slopes before adjustment for multiple testing (PTyr = 0.003; PHDLrelated < 0.05), but not after. CONCLUSIONS: This study identified metabolites associated with impaired kidney function in T2D, implying involvement of lipid and amino acid metabolism in the pathogenesis. Whether these processes precede or are consequences of renal impairment needs further investigation

Bivariate genome-wide association meta-analysis of pediatric musculoskeletal traits reveals pleiotropic effects at the SREBF1/TOM1L2 locus

Bone mineral density is known to be a heritable, polygenic trait whereas genetic variants contributing to lean mass variation remain largely unknown. We estimated the shared SNP heritability and performed a bivariate GWAS meta-analysis of total-body lean mass (TB-LM) and total-body less head bone mineral density (TBLH-BMD) regions in 10,414 children. The estimated SNP heritability is 43% for TBLH-BMD, and 39% for TB-LM, with a shared genetic component of 43%. We identify variants with pleiotropic effects in eight loci, including seven established bone mineral density loci: _WNT4, GALNT3, MEPE, CPED1/WNT16, TNFSF11, RIN3, and PPP6R3/LRP5_. Variants in the _TOM1L2/SREBF1_ locus exert opposing effects TB-LM and TBLH-BMD, and have a stronger association with the former trait. We show that _SREBF1_ is expressed in murine and human osteoblasts, as well as in human muscle tissue. This is the first bivariate GWAS meta-analysis to demonstrate genetic factors with pleiotropic effects on bone mineral density and lean mass

Genome-wide meta-analysis uncovers novel loci influencing circulating leptin levels.

Leptin is an adipocyte-secreted hormone, the circulating levels of which correlate closely with overall adiposity. Although rare mutations in the leptin (LEP) gene are well known to cause leptin deficiency and severe obesity, no common loci regulating circulating leptin levels have been uncovered. Therefore, we performed a genome-wide association study (GWAS) of circulating leptin levels from 32,161 individuals and followed up loci reaching P&lt;10(-6) in 19,979 additional individuals. We identify five loci robustly associated (P&lt;5 × 10(-8)) with leptin levels in/near LEP, SLC32A1, GCKR, CCNL1 and FTO. Although the association of the FTO obesity locus with leptin levels is abolished by adjustment for BMI, associations of the four other loci are independent of adiposity. The GCKR locus was found associated with multiple metabolic traits in previous GWAS and the CCNL1 locus with birth weight. Knockdown experiments in mouse adipose tissue explants show convincing evidence for adipogenin, a regulator of adipocyte differentiation, as the novel causal gene in the SLC32A1 locus influencing leptin levels. Our findings provide novel insights into the regulation of leptin production by adipose tissue and open new avenues for examining the influence of variation in leptin levels on adiposity and metabolic health

Multiethnic meta-analysis identifies ancestry-specific and cross-ancestry loci for pulmonary function

Nearly 100 loci have been identified for pulmonary function, almost exclusively in studies of European ancestry populations. We extend previous research by meta-analyzing genome-wide association studies of 1000 Genomes imputed variants in relation to pulmonary function in a multiethnic population of 90,715 individuals of European (N = 60,552), African (N = 8429), Asian (N = 9959), and Hispanic/Latino (N = 11,775) ethnicities. We identify over 50 additional loci at genome-wide significance in ancestry-specific or multiethnic meta-analyses. Using recent fine-mapping methods incorporating functional annotation, gene expression, and differences in linkage disequilibrium between ethnicities, we further shed light on potential causal variants and genes at known and newly identified loci. Several of the novel genes encode proteins with predicted or established drug targets, including KCNK2 and CDK12. Our study highlights the utility of multiethnic and integrative genomics approaches to extend existing knowledge of the genetics of l

Trans-ancestry genome-wide association study identifies 12 genetic loci influencing blood pressure and implicates a role for DNA methylation

We carried out a trans-ancestry genome-wide association and replication study of blood pressure phenotypes among up to 320,251 individuals of East Asian, European and South Asian ancestry. We find genetic variants at 12 new loci to be associated with blood pressure (P = 3.9 &times; 10-11 to 5.0 &times; 10-21). The sentinel blood pressure SNPs are enriched for association with DNA methylation at multiple nearby CpG sites, suggesting that, at some of the loci identified, DNA methylation may lie on the regulatory pathway linking sequence variation to blood pressure. The sentinel SNPs at the 12 new loci point to genes involved in vascular smooth muscle (IGFBP3, KCNK3, PDE3A and PRDM6) and renal (ARHGAP24, OSR1, SLC22A7 and TBX2) function. The new and known genetic variants predict increased left ventricular mass, circulating levels of NT-proBNP, and cardiovascular and all-cause mortality (P = 0.04 to 8.6 &times; 10-6). Our results provide new evidence for the role of DNA methylation in blood pressure regulation

Genome-wide analyses identify a role for SLC17A4 and AADAT in thyroid hormone regulation.

Thyroid dysfunction is an important public health problem, which affects 10% of the general population and increases the risk of cardiovascular morbidity and mortality. Many aspects of thyroid hormone regulation have only partly been elucidated, including its transport, metabolism, and genetic determinants. Here we report a large meta-analysis of genome-wide association studies for thyroid function and dysfunction, testing 8 million genetic variants in up to 72,167 individuals. One-hundred-and-nine independent genetic variants are associated with these traits. A genetic risk score, calculated to assess their combined effects on clinical end points, shows significant associations with increased risk of both overt (Graves' disease) and subclinical thyroid disease, as well as clinical complications. By functional follow-up on selected signals, we identify a novel thyroid hormone transporter (SLC17A4) and a metabolizing enzyme (AADAT). Together, these results provide new knowledge about thyroid hormone physiology and disease, opening new possibilities for therapeutic targets

Genome-Wide Association Study of Apparent Treatment-Resistant Hypertension in the CHARGE Consortium: The CHARGE Pharmacogenetics Working Group

BACKGROUND: Only a handful of genetic discovery efforts in apparent treatment-resistant hypertension (aTRH) have been described. METHODS: We conducted a case-control genome-wide association study of aTRH among persons treated for hypertension, using data from 10 cohorts of European ancestry (EA) and 5 cohorts of African ancestry (AA). Cases were treated with 3 different antihypertensive medication classes and had blood pressure (BP) above goal (systolic BP ≥ 140 mm Hg and/or diastolic BP ≥ 90 mm Hg) or 4 or more medication classes regardless of BP control (nEA = 931, nAA = 228). Both a normotensive control group and a treatment-responsive control group were considered in separate analyses. Normotensive controls were untreated (nEA = 14,210, nAA = 2,480) and had systolic BP/diastolic BP &lt; 140/90 mm Hg. Treatment-responsive controls (nEA = 5,266, nAA = 1,817) had BP at goal (&lt;140/90 mm Hg), while treated with one antihypertensive medication class. Individual cohorts used logistic regression with adjustment for age, sex, study site, and principal components for ancestry to examine the association of single-nucleotide polymorphisms with case-control status. Inverse variance-weighted fixed-effects meta-analyses were carried out using METAL. RESULTS: The known hypertension locus, CASZ1, was a top finding among EAs (P = 1.1 × 10-8) and in the race-combined analysis (P = 1.5 × 10-9) using the normotensive control group (rs12046278, odds ratio = 0.71 (95% confidence interval: 0.6-0.8)). Single-nucleotide polymorphisms in this locus were robustly replicated in the Million Veterans Program (MVP) study in consideration of a treatment-responsive control group. There were no statistically significant findings for the discovery analyses including treatment-responsive controls. CONCLUSION: This genomic discovery effort for aTRH identified CASZ1 as an aTRH risk locus

Bone mineral density loci specific to the skull portray potential pleiotropic effects on craniosynostosis

Skull bone mineral density (SK-BMD) provides a suitable trait for the discovery of key genes in bone biology, particularly to intramembranous ossification, not captured at other skeletal sites. We perform a genome-wide association meta-analysis (n ~ 43,800) of SK-BMD, identifying 59 loci, collectively explaining 12.5% of the trait variance. Association signals cluster within gene-sets involved in skeletal development and osteoporosis. Among the four novel loci (ZIC1, PRKAR1A, AZIN1/ATP6V1C1, GLRX3), there are factors implicated in intramembranous ossification and as we show, inherent to craniosynostosis processes. Functional follow-up in zebrafish confirms the importance of ZIC1 on cranial suture patterning. Likewise, we observe abnormal cranial bone initiation that culminates in ectopic sutures and reduced BMD in mosaic atp6v1c1 knockouts. Mosaic prkar1a knockouts present asymmetric bone growth and, conversely, elevated BMD. In light of this evidence linking SK-BMD loci to craniofacial abnormalities, our study provides new insight into the pathophysiology, diagnosis and treatment of skeletal diseases. Prevention, Population and Disease management (PrePoD)Public Health and primary car

FTO genetic variants, dietary intake and body mass index: insights from 177,330 individuals.

FTO is the strongest known genetic susceptibility locus for obesity. Experimental studies in animals suggest the potential roles of FTO in regulating food intake. The interactive relation among FTO variants, dietary intake and body mass index (BMI) is complex and results from previous often small-scale studies in humans are highly inconsistent. We performed large-scale analyses based on data from 177,330 adults (154 439 Whites, 5776 African Americans and 17 115 Asians) from 40 studies to examine: (i) the association between the FTO-rs9939609 variant (or a proxy single-nucleotide polymorphism) and total energy and macronutrient intake and (ii) the interaction between the FTO variant and dietary intake on BMI. The minor allele (A-allele) of the FTO-rs9939609 variant was associated with higher BMI in Whites (effect per allele = 0.34 [0.31, 0.37] kg/m(2), P = 1.9 × 10(-105)), and all participants (0.30 [0.30, 0.35] kg/m(2), P = 3.6 × 10(-107)). The BMI-increasing allele of the FTO variant showed a significant association with higher dietary protein intake (effect per allele = 0.08 [0.06, 0.10] %, P = 2.4 × 10(-16)), and relative weak associations with lower total energy intake (-6.4 [-10.1, -2.6] kcal/day, P = 0.001) and lower dietary carbohydrate intake (-0.07 [-0.11, -0.02] %, P = 0.004). The associations with protein (P = 7.5 × 10(-9)) and total energy (P = 0.002) were attenuated but remained significant after adjustment for BMI. We did not find significant interactions between the FTO variant and dietary intake of total energy, protein, carbohydrate or fat on BMI. Our findings suggest a positive association between the BMI-increasing allele of FTO variant and higher dietary protein intake and offer insight into potential link between FTO, dietary protein intake and adiposity

Mendelian randomization analysis does not support causal associations of birth weight with hypertension risk and blood pressure in adulthood

Epidemiology studies suggested that low birthweight was associated with a higher risk of hypertension in later life. However, little is known about the causality of such associations. In our study, we evaluated the causal association of low birthweight with adulthood hypertension following a standard analytic protocol using the study-level data of 183,433 participants from 60 studies (CHARGE-BIG consortium), as well as that with blood pressure using publicly available summary-level genome-wide association data from EGG consortium of 153,781 participants, ICBP consortium and UK Biobank cohort together of 757,601 participants. We used seven SNPs as the instrumental variable in the study-level analysis and 47 SNPs in the summary-level analysis. In the study-level analyses, decreased birthweight was associated with a higher risk of hypertension in adults (the odds ratio per 1 standard deviation (SD) lower birthweight, 1.22; 95% CI 1.16 to 1.28), while no association was found between genetically instrumented birthweight and hypertension risk (instrumental odds ratio for causal effect per 1 SD lower birthweight, 0.97; 95% CI 0.68 to 1.41). Such results were consistent with that from the summary-level analyses, where the genetically determined low birthweight was not associated with blood pressure measurements either. One SD lower genetically determined birthweight was not associated with systolic blood pressure (β = − 0.76, 95% CI − 2.45 to 1.08 mmHg), 0.06 mmHg lower diastolic blood pressure (β = − 0.06, 95% CI − 0.93 to 0.87 mmHg), or pulse pressure (β = − 0.65, 95% CI − 1.38 to 0.69 mmHg, all p > 0.05). Our findings suggest that the inverse association of birthweight with hypertension risk from observational studies was not supported by large Mendelian randomization analyses