Genetically proxied therapeutic prolyl-hydroxylase inhibition and cardiovascular risk.

Prolyl hydroxylase (PHD) inhibitors are in clinical development for anaemia in chronic kidney disease. Epidemiological studies have reported conflicting results regarding safety of long-term therapeutic haemoglobin (Hgb) rises through PHD inhibition on risk of cardiovascular disease. Genetic variation in genes encoding PHDs can be used as partial proxies to investigate the potential effects of long-term Hgb rises. We used Mendelian randomization to investigate the effect of long-term Hgb level rises through genetically proxied PHD inhibition on coronary artery disease (CAD: 60 801 cases; 123 504 controls), myocardial infarction (MI: 42 561 cases; 123 504 controls) or stroke (40 585 cases; 406 111 controls). To further characterize long-term effects of Hgb level rises, we performed a phenome-wide association study (PheWAS) in up to 451 099 UK Biobank individuals. Genetically proxied therapeutic PHD inhibition, equivalent to a 1.00 g/dl increase in Hgb levels, was not associated (at P < 0.05) with increased odds of CAD; odd ratio (OR) [95% confidence intervals (CI)] = 1.06 (0.84, 1.35), MI [OR (95% CI) = 1.02 (0.79, 1.33)] or stroke [OR (95% CI) = 0.91 (0.66, 1.24)]. PheWAS revealed associations with blood related phenotypes consistent with EGLN's role, relevant kidney- and liver-related biomarkers like estimated glomerular filtration rate and microalbuminuria, and non-alcoholic fatty liver disease (Bonferroni-adjusted P < 5.42E-05) but these were not clinically meaningful. These findings suggest that long-term alterations in Hgb through PHD inhibition are unlikely to substantially increase cardiovascular disease risk; using large disease genome-wide association study data, we could exclude ORs of 1.35 for cardiovascular risk with a 1.00 g/dl increase in Hgb.

Identification and single-base gene-editing functional validation of a cis-EPO variant as a genetic predictor for EPO-increasing therapies.

Hypoxia-inducible factor prolyl hydroxylase inhibitors (HIF-PHIs) are currently under clinical development for treating anemia in chronic kidney disease (CKD), but it is important to monitor their cardiovascular safety. Genetic variants can be used as predictors to help inform the potential risk of adverse effects associated with drug treatments. We therefore aimed to use human genetics to help assess the risk of adverse cardiovascular events associated with therapeutically altered EPO levels to help inform clinical trials studying the safety of HIF-PHIs. By performing a genome-wide association meta-analysis of EPO (n = 6,127), we identified a cis-EPO variant (rs1617640) lying in the EPO promoter region. We validated this variant as most likely causal in controlling EPO levels by using genetic and functional approaches, including single-base gene editing. Using this variant as a partial predictor for therapeutic modulation of EPO and large genome-wide association data in Mendelian randomization tests, we found no evidence (at p < 0.05) that genetically predicted long-term rises in endogenous EPO, equivalent to a 2.2-unit increase, increased risk of coronary artery disease (CAD, OR [95% CI] = 1.01 [0.93, 1.07]), myocardial infarction (MI, OR [95% CI] = 0.99 [0.87, 1.15]), or stroke (OR [95% CI] = 0.97 [0.87, 1.07]). We could exclude increased odds of 1.15 for cardiovascular disease for a 2.2-unit EPO increase. A combination of genetic and functional studies provides a powerful approach to investigate the potential therapeutic profile of EPO-increasing therapies for treating anemia in CKD.

Genetic loci and prioritization of genes for kidney function decline derived from a meta-analysis of 62 longitudinal genome-wide association studies.

Estimated glomerular filtration rate (eGFR) reflects kidney function. Progressive eGFR-decline can lead to kidney failure, necessitating dialysis or transplantation. Hundreds of loci from genome-wide association studies (GWAS) for eGFR help explain population cross section variability. Since the contribution of these or other loci to eGFR-decline remains largely unknown, we derived GWAS for annual eGFR-decline and meta-analyzed 62 longitudinal studies with eGFR assessed twice over time in all 343,339 individuals and in high-risk groups. We also explored different covariate adjustment. Twelve genome-wide significant independent variants for eGFR-decline unadjusted or adjusted for eGFR-baseline (11 novel, one known for this phenotype), including nine variants robustly associated across models were identified. All loci for eGFR-decline were known for cross-sectional eGFR and thus distinguished a subgroup of eGFR loci. Seven of the nine variants showed variant-by-age interaction on eGFR cross section (further about 350,000 individuals), which linked genetic associations for eGFR-decline with age-dependency of genetic cross-section associations. Clinically important were two to four-fold greater genetic effects on eGFR-decline in high-risk subgroups. Five variants associated also with chronic kidney disease progression mapped to genes with functional in-silico evidence (UMOD, SPATA7, GALNTL5, TPPP). An unfavorable versus favorable nine-variant genetic profile showed increased risk odds ratios of 1.35 for kidney failure (95% confidence intervals 1.03-1.77) and 1.27 for acute kidney injury (95% confidence intervals 1.08-1.50) in over 2000 cases each, with matched controls). Thus, we provide a large data resource, genetic loci, and prioritized genes for kidney function decline, which help inform drug development pipelines revealing important insights into the age-dependency of kidney function genetics.

Associations of Mitochondrial and Nuclear Mitochondrial Variants and Genes with Seven Metabolic Traits.

Mitochondria (MT), the major site of cellular energy production, are under dual genetic control by 37 mitochondrial DNA (mtDNA) genes and numerous nuclear genes (MT-nDNA). In the CHARGEmtDNA+ Consortium, we studied genetic associations of mtDNA and MT-nDNA associations with body mass index (BMI), waist-hip-ratio (WHR), glucose, insulin, HOMA-B, HOMA-IR, and HbA1c. This 45-cohort collaboration comprised 70,775 (insulin) to 170,202 (BMI) pan-ancestry individuals. Validation and imputation of mtDNA variants was followed by single-variant and gene-based association testing. We report two significant common variants, one in MT-ATP6 associated (p ≤ 5E-04) with WHR and one in the D-loop with glucose. Five rare variants in MT-ATP6, MT-ND5, and MT-ND6 associated with BMI, WHR, or insulin. Gene-based meta-analysis identified MT-ND3 associated with BMI (p ≤ 1E-03). We considered 2,282 MT-nDNA candidate gene associations compiled from online summary results for our traits (20 unique studies with 31 dataset consortia's genome-wide associations [GWASs]). Of these, 109 genes associated (p ≤ 1E-06) with at least 1 of our 7 traits. We assessed regulatory features of variants in the 109 genes, cis- and trans-gene expression regulation, and performed enrichment and protein-protein interactions analyses. Of the identified mtDNA and MT-nDNA genes, 79 associated with adipose measures, 49 with glucose/insulin, 13 with risk for type 2 diabetes, and 18 with cardiovascular disease, indicating for pleiotropic effects with health implications. Additionally, 21 genes related to cholesterol, suggesting additional important roles for the genes identified. Our results suggest that mtDNA and MT-nDNA genes and variants reported make important contributions to glucose and insulin metabolism, adipocyte regulation, diabetes, and cardiovascular disease.

Meta-analysis uncovers genome-wide significant variants for rapid kidney function decline.

Rapid decline of glomerular filtration rate estimated from creatinine (eGFRcrea) is associated with severe clinical endpoints. In contrast to cross-sectionally assessed eGFRcrea, the genetic basis for rapid eGFRcrea decline is largely unknown. To help define this, we meta-analyzed 42 genome-wide association studies from the Chronic Kidney Diseases Genetics Consortium and United Kingdom Biobank to identify genetic loci for rapid eGFRcrea decline. Two definitions of eGFRcrea decline were used: 3 mL/min/1.73m2/year or more ("Rapid3"; encompassing 34,874 cases, 107,090 controls) and eGFRcrea decline 25% or more and eGFRcrea under 60 mL/min/1.73m2 at follow-up among those with eGFRcrea 60 mL/min/1.73m2 or more at baseline ("CKDi25"; encompassing 19,901 cases, 175,244 controls). Seven independent variants were identified across six loci for Rapid3 and/or CKDi25: consisting of five variants at four loci with genome-wide significance (near UMOD-PDILT (2), PRKAG2, WDR72, OR2S2) and two variants among 265 known eGFRcrea variants (near GATM, LARP4B). All these loci were novel for Rapid3 and/or CKDi25 and our bioinformatic follow-up prioritized variants and genes underneath these loci. The OR2S2 locus is novel for any eGFRcrea trait including interesting candidates. For the five genome-wide significant lead variants, we found supporting effects for annual change in blood urea nitrogen or cystatin-based eGFR, but not for GATM or LARP4B. Individuals at high compared to those at low genetic risk (8-14 vs. 0-5 adverse alleles) had a 1.20-fold increased risk of acute kidney injury (95% confidence interval 1.08-1.33). Thus, our identified loci for rapid kidney function decline may help prioritize therapeutic targets and identify mechanisms and individuals at risk for sustained deterioration of kidney function.

Genome-wide Trans-ethnic Meta-analysis Identifies Seven Genetic Loci Influencing Erythrocyte Traits and a Role for RBPMS in Erythropoiesis.

Genome-wide association studies (GWASs) have identified loci for erythrocyte traits in primarily European ancestry populations. We conducted GWAS meta-analyses of six erythrocyte traits in 71,638 individuals from European, East Asian, and African ancestries using a Bayesian approach to account for heterogeneity in allelic effects and variation in the structure of linkage disequilibrium between ethnicities. We identified seven loci for erythrocyte traits including a locus (RBPMS/GTF2E2) associated with mean corpuscular hemoglobin and mean corpuscular volume. Statistical fine-mapping at this locus pointed to RBPMS at this locus and excluded nearby GTF2E2. Using zebrafish morpholino to evaluate loss of function, we observed a strong in vivo erythropoietic effect for RBPMS but not for GTF2E2, supporting the statistical fine-mapping at this locus and demonstrating that RBPMS is a regulator of erythropoiesis. Our findings show the utility of trans-ethnic GWASs for discovery and characterization of genetic loci influencing hematologic traits.

Group IIA Secretory Phospholipase A2, Vascular Inflammation, and Incident Cardiovascular Disease.

Objective- Inflammation is a causal risk factor for cardiovascular disease (CVD). sPLA2-IIA (group IIA secretory phospholipase A2) plays an integral role in regulating vascular inflammation. Although studies investigated sPLA2-IIA in secondary prevention, we prospectively evaluated sPLA2-IIA mass and genetic variants with CVD events in a primary prevention population with chronic inflammation. Approach and Results- The JUPITER trial (Justification for the Use of Statins in Prevention: An Intervention Trial Evaluating Rosuvastatin) randomized participants with LDL (low-density lipoprotein) <130 mg/dL and hsCRP (high-sensitivity C-reactive protein) ≥2 mg/L to high-intensity rosuvastatin versus placebo. Baseline and 1-year plasma sPLA2-IIA mass was measured (N=11 269 baseline; N=9620 1 year). We also identified genetic variants influencing sPLA2-IIA using genome-wide association and examined them with CVD. Three hundred thirteen incident CVD events occurred during follow-up. Baseline sPLA2-IIA mass (median, 25th-75th percentile: 3.81, 2.49-6.03 ng/mL) was associated with increased risk of CVD: risk factor-adjusted hazard ratio (95% CI; P) per SD increment: 1.22 (1.08-1.38; P=0.002). This remained significant (1.18; 1.04-1.35; P=0.01) after incrementally adjusting for hsCRP. Similar estimates were observed in rosuvastatin and placebo groups ( P treatment interaction>0.05). The rs11573156C variant in PLA2G2A (encoding sPLA2-IIA) had the strongest effect on sPLA2-II: median (25th-75th percentile, ng/mL) for CC and GG genotypes: 2.79 (1.97-4.01) and 7.38 (5.38-10.19), respectively; and had nonsignificant trend for higher CVD risk (hazard ratio, 1.11; 95% CI, 0.89-1.38; P=0.34). Conclusions- In the JUPITER population recruited on chronic inflammation, sPLA2-IIA mass was associated with CVD risk relating to vascular inflammation not fully reflected by hsCRP. Additional studies, including larger functional genetic and clinical studies, are needed to determine whether sPLA2-IIA may be a potential pharmacological target for primary prevention of CVD. Clinical Trial Registration- URL: http://www.clinicaltrials.gov . Unique identifier: NCT00239681.

Ranking and characterization of established BMI and lipid associated loci as candidates for gene-environment interactions.

Phenotypic variance heterogeneity across genotypes at a single nucleotide polymorphism (SNP) may reflect underlying gene-environment (G×E) or gene-gene interactions. We modeled variance heterogeneity for blood lipids and BMI in up to 44,211 participants and investigated relationships between variance effects (Pv), G×E interaction effects (with smoking and physical activity), and marginal genetic effects (Pm). Correlations between Pv and Pm were stronger for SNPs with established marginal effects (Spearman's ρ = 0.401 for triglycerides, and ρ = 0.236 for BMI) compared to all SNPs. When Pv and Pm were compared for all pruned SNPs, only BMI was statistically significant (Spearman's ρ = 0.010). Overall, SNPs with established marginal effects were overrepresented in the nominally significant part of the Pv distribution (Pbinomial <0.05). SNPs from the top 1% of the Pm distribution for BMI had more significant Pv values (PMann-Whitney = 1.46×10-5), and the odds ratio of SNPs with nominally significant (<0.05) Pm and Pv was 1.33 (95% CI: 1.12, 1.57) for BMI. Moreover, BMI SNPs with nominally significant G×E interaction P-values (Pint<0.05) were enriched with nominally significant Pv values (Pbinomial = 8.63×10-9 and 8.52×10-7 for SNP × smoking and SNP × physical activity, respectively). We conclude that some loci with strong marginal effects may be good candidates for G×E, and variance-based prioritization can be used to identify them.

Correction: Genome-wide physical activity interactions in adiposity - A meta-analysis of 200,452 adults.

[This corrects the article DOI: 10.1371/journal.pgen.1006528.].

Genome-wide physical activity interactions in adiposity - A meta-analysis of 200,452 adults.

Physical activity (PA) may modify the genetic effects that give rise to increased risk of obesity. To identify adiposity loci whose effects are modified by PA, we performed genome-wide interaction meta-analyses of BMI and BMI-adjusted waist circumference and waist-hip ratio from up to 200,452 adults of European (n = 180,423) or other ancestry (n = 20,029). We standardized PA by categorizing it into a dichotomous variable where, on average, 23% of participants were categorized as inactive and 77% as physically active. While we replicate the interaction with PA for the strongest known obesity-risk locus in the FTO gene, of which the effect is attenuated by ~30% in physically active individuals compared to inactive individuals, we do not identify additional loci that are sensitive to PA. In additional genome-wide meta-analyses adjusting for PA and interaction with PA, we identify 11 novel adiposity loci, suggesting that accounting for PA or other environmental factors that contribute to variation in adiposity may facilitate gene discovery.

Integrating genetic, transcriptional, and biological information provides insights into obesity.

OBJECTIVE: Indices of body fat distribution are heritable, but few genetic signals have been reported from genome-wide association studies (GWAS) of computed tomography (CT) imaging measurements of body fat distribution. We aimed to identify genes associated with adiposity traits and the key drivers that are central to adipose regulatory networks. SUBJECTS: We analyzed gene transcript expression data in blood from participants in the Framingham Heart Study, a large community-based cohort (n up to 4303), as well as implemented an integrative analysis of these data and existing biological information. RESULTS: Our association analyses identified unique and common gene expression signatures across several adiposity traits, including body mass index, waist-hip ratio, waist circumference, and CT-measured indices, including volume and quality of visceral and subcutaneous adipose tissues. We identified six enriched KEGG pathways and two co-expression modules for further exploration of adipose regulatory networks. The integrative analysis revealed four gene sets (Apoptosis, p53 signaling pathway, Proteasome, Ubiquitin-mediated proteolysis) and two co-expression modules with significant genetic variants and 94 key drivers/genes whose local networks were enriched with adiposity-associated genes, suggesting that these enriched pathways or modules have genetic effects on adiposity. Most identified key driver genes are involved in essential biological processes such as controlling cell cycle, DNA repair, and degradation of regulatory proteins are cancer related. CONCLUSIONS: Our integrative analysis of genetic, transcriptional, and biological information provides a list of compelling candidates for further follow-up functional studies to uncover the biological mechanisms underlying obesity. These candidates highlight the value of examining CT-derived and central adiposity traits.

Impact of common genetic determinants of Hemoglobin A1c on type 2 diabetes risk and diagnosis in ancestrally diverse populations: A transethnic genome-wide meta-analysis.

BACKGROUND: Glycated hemoglobin (HbA1c) is used to diagnose type 2 diabetes (T2D) and assess glycemic control in patients with diabetes. Previous genome-wide association studies (GWAS) have identified 18 HbA1c-associated genetic variants. These variants proved to be classifiable by their likely biological action as erythrocytic (also associated with erythrocyte traits) or glycemic (associated with other glucose-related traits). In this study, we tested the hypotheses that, in a very large scale GWAS, we would identify more genetic variants associated with HbA1c and that HbA1c variants implicated in erythrocytic biology would affect the diagnostic accuracy of HbA1c. We therefore expanded the number of HbA1c-associated loci and tested the effect of genetic risk-scores comprised of erythrocytic or glycemic variants on incident diabetes prediction and on prevalent diabetes screening performance. Throughout this multiancestry study, we kept a focus on interancestry differences in HbA1c genetics performance that might influence race-ancestry differences in health outcomes. METHODS & FINDINGS: Using genome-wide association meta-analyses in up to 159,940 individuals from 82 cohorts of European, African, East Asian, and South Asian ancestry, we identified 60 common genetic variants associated with HbA1c. We classified variants as implicated in glycemic, erythrocytic, or unclassified biology and tested whether additive genetic scores of erythrocytic variants (GS-E) or glycemic variants (GS-G) were associated with higher T2D incidence in multiethnic longitudinal cohorts (N = 33,241). Nineteen glycemic and 22 erythrocytic variants were associated with HbA1c at genome-wide significance. GS-G was associated with higher T2D risk (incidence OR = 1.05, 95% CI 1.04-1.06, per HbA1c-raising allele, p = 3 × 10-29); whereas GS-E was not (OR = 1.00, 95% CI 0.99-1.01, p = 0.60). In Europeans and Asians, erythrocytic variants in aggregate had only modest effects on the diagnostic accuracy of HbA1c. Yet, in African Americans, the X-linked G6PD G202A variant (T-allele frequency 11%) was associated with an absolute decrease in HbA1c of 0.81%-units (95% CI 0.66-0.96) per allele in hemizygous men, and 0.68%-units (95% CI 0.38-0.97) in homozygous women. The G6PD variant may cause approximately 2% (N = 0.65 million, 95% CI 0.55-0.74) of African American adults with T2D to remain undiagnosed when screened with HbA1c. Limitations include the smaller sample sizes for non-European ancestries and the inability to classify approximately one-third of the variants. Further studies in large multiethnic cohorts with HbA1c, glycemic, and erythrocytic traits are required to better determine the biological action of the unclassified variants. CONCLUSIONS: As G6PD deficiency can be clinically silent until illness strikes, we recommend investigation of the possible benefits of screening for the G6PD genotype along with using HbA1c to diagnose T2D in populations of African ancestry or groups where G6PD deficiency is common. Screening with direct glucose measurements, or genetically-informed HbA1c diagnostic thresholds in people with G6PD deficiency, may be required to avoid missed or delayed diagnoses.

Genetic risk score and risk of stage 3 chronic kidney disease.

BACKGROUND: We developed a genetic risk score (GRS) and examined whether the GRS may predict incident stage 3 chronic kidney disease (CKD) independent of common clinical risk factors. METHOD: The present study included 2,698 individuals who attended the 15th (1977 to 1979) and the 24th exams (1995 to 1998) in the Framingham Original cohort or the 6th (1995 to 1998) and the 8th exams (2005 to 2008) in the Framingham Offspring cohort. A weighted GRS was constructed combining 53 single nucleotide polymorphisms (SNPs) associated with lower creatinine-based estimated glomerular filtration rate (eGFR). Stage 3 CKD was defined as eGFR <60 mL/min/1.73 m2, and incident cases were identified at follow-up after excluding prevalent cases at baseline. RESULTS: A total of 292 incident cases and 2,406 non-cases were identified over, on average, 11 years of follow-up. After adjustment for sex, age, cohort, baseline eGFR, hypertension, diabetes, and dipstick proteinuria, the odds ratio of incident stage 3 CKD was 1.37 (95%CI: 1.02-1.83) per 10 alleles of the GRS (P = 0.04). There was no statistically significant difference between the C-statistic without and with inclusion of the GRS (0.783 and 0.785, respectively; P = 0.39). CONCLUSIONS: A GRS developed based on 53 SNPs associated with reduced eGFR was prospectively associated with incident stage 3 CKD. However, this score did not substantially improve discrimination of stage 3 CKD beyond the common clinical risk factors.

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.

Dietary Guideline Adherence Index and Kidney Measures in the Framingham Heart Study.

BACKGROUND: No observational studies have directly considered dietary guidelines when examining the prospective association between dietary intake and kidney measures. STUDY DESIGN: Prospective cohort study. SETTING & PARTICIPANTS: We examined participants who attended examinations 7 (1998-2001) and 8 (2005-2008) in the Framingham Offspring Cohort. PREDICTORS: Individual components of Dietary Guidelines for Americans Adherence Index (DGAI) that reflect adherence to key dietary recommendations based on the 2005 guideline. OUTCOMES & MEASURES: The primary outcome was incident low estimated glomerular filtration rate (eGFR) at follow-up after exclusion of prevalent low eGFR at baseline. Low eGFR was defined as serum creatinine-based eGFR<60mL/min/1.73m2. RESULTS: Among 1,822 participants (mean age, 59.4 years; 54.6% women), 181 incident cases of low eGFR were identified. After adjustment for potential confounders, compared to optimal adherence to meat and legume recommendations, low adherence was associated with higher odds of incident low eGFR (P for trend = 0.01); ORs in the lowest and intermediate adherence categories were 2.98 (95% CI, 1.13-7.92) and 1.65 (95% CI, 1.02-2.66), respectively. Low adherence to dairy product recommendations was also associated with higher odds of incident low eGFR compared to optimal adherence (P for trend = 0.03); ORs in the lowest and intermediate adherence categories were 1.98 (95% CI, 1.03-3.82) and 1.59 (95% CI, 0.81-3.11), respectively. In addition, low adherence to meat and legume recommendations was associated with rapid eGFR decline (P for trend = 0.01), and low adherence to dairy product recommendations was associated with rapid eGFR decline (P for trend = 0.01) and incident albuminuria (P for trend = 0.03). LIMITATIONS: The DGAI was developed based on the 2005 Dietary Guidelines for Americans. CONCLUSIONS: Better adherence to dietary recommendations for both meat and legumes and dairy products was associated with lower risk for developing adverse kidney measures.

A null mutation in ANGPTL8 does not associate with either plasma glucose or type 2 diabetes in humans.

BACKGROUND: Experiments in mice initially suggested a role for the protein angiopoietin-like 8 (ANGPTL8) in glucose homeostasis. However, subsequent experiments in model systems have challenged this proposed role. We sought to better understand the importance of ANGPTL8 in human glucose homeostasis by examining the association of a null mutation in ANGPTL8 with fasting glucose levels and risk for type 2 diabetes. METHODS: A naturally-occurring null mutation in human ANGPTL8 (rs145464906; c.361C > T; p.Q121X) is carried by ~1 in 1000 individuals of European ancestry and is associated with higher levels of plasma high-density lipoprotein cholesterol, suggesting that this mutation has functional significance. We examined the association of p.Q121X with fasting glucose levels and risk for type 2 diabetes in up to 95,558 individuals (14,824 type 2 diabetics and 80,734 controls). RESULTS: We found no significant association of p.Q121X with either fasting glucose or type 2 diabetes (p-value = 0.90 and 0.65, respectively). Given our sample sizes, we had >98 % power to detect at least a 0.23 mmol/L effect on plasma glucose and >95 % power to detect a 70 % increase in risk for type 2 diabetes. CONCLUSION: Disruption of ANGPTL8 function in humans does not seem to have a large effect on measures of glucose tolerance.

Gene × physical activity interactions in obesity: combined analysis of 111,421 individuals of European ancestry.

Numerous obesity loci have been identified using genome-wide association studies. A UK study indicated that physical activity may attenuate the cumulative effect of 12 of these loci, but replication studies are lacking. Therefore, we tested whether the aggregate effect of these loci is diminished in adults of European ancestry reporting high levels of physical activity. Twelve obesity-susceptibility loci were genotyped or imputed in 111,421 participants. A genetic risk score (GRS) was calculated by summing the BMI-associated alleles of each genetic variant. Physical activity was assessed using self-administered questionnaires. Multiplicative interactions between the GRS and physical activity on BMI were tested in linear and logistic regression models in each cohort, with adjustment for age, age(2), sex, study center (for multicenter studies), and the marginal terms for physical activity and the GRS. These results were combined using meta-analysis weighted by cohort sample size. The meta-analysis yielded a statistically significant GRS × physical activity interaction effect estimate (Pinteraction  = 0.015). However, a statistically significant interaction effect was only apparent in North American cohorts (n = 39,810, Pinteraction  = 0.014 vs. n = 71,611, Pinteraction  = 0.275 for Europeans). In secondary analyses, both the FTO rs1121980 (Pinteraction  = 0.003) and the SEC16B rs10913469 (Pinteraction  = 0.025) variants showed evidence of SNP × physical activity interactions. This meta-analysis of 111,421 individuals provides further support for an interaction between physical activity and a GRS in obesity disposition, although these findings hinge on the inclusion of cohorts from North America, indicating that these results are either population-specific or non-causal.

Genetic loci associated with circulating levels of very long-chain saturated fatty acids.

Very long-chain saturated fatty acids (VLSFAs) are saturated fatty acids with 20 or more carbons. In contrast to the more abundant saturated fatty acids, such as palmitic acid, there is growing evidence that circulating VLSFAs may have beneficial biological properties. Whether genetic factors influence circulating levels of VLSFAs is not known. We investigated the association of common genetic variation with plasma phospholipid/erythrocyte levels of three VLSFAs by performing genome-wide association studies in seven population-based cohorts comprising 10,129 subjects of European ancestry. We observed associations of circulating VLSFA concentrations with common variants in two genes, serine palmitoyl-transferase long-chain base subunit 3 (SPTLC3), a gene involved in the rate-limiting step of de novo sphingolipid synthesis, and ceramide synthase 4 (CERS4). The SPTLC3 variant at rs680379 was associated with higher arachidic acid (20:0 , P = 5.81 × 10(-13)). The CERS4 variant at rs2100944 was associated with higher levels of 20:0 (P = 2.65 × 10(-40)) and in analyses that adjusted for 20:0, with lower levels of behenic acid (P = 4.22 × 10(-26)) and lignoceric acid (P = 3.20 × 10(-21)). These novel associations suggest an inter-relationship of circulating VLSFAs and sphingolipid synthesis.

Genome-wide association study of kidney function decline in individuals of European descent.

Genome-wide association studies (GWASs) have identified multiple loci associated with cross-sectional eGFR, but a systematic genetic analysis of kidney function decline over time is missing. Here we conducted a GWAS meta-analysis among 63,558 participants of European descent, initially from 16 cohorts with serial kidney function measurements within the CKDGen Consortium, followed by independent replication among additional participants from 13 cohorts. In stage 1 GWAS meta-analysis, single-nucleotide polymorphisms (SNPs) at MEOX2, GALNT11, IL1RAP, NPPA, HPCAL1, and CDH23 showed the strongest associations for at least one trait, in addition to the known UMOD locus, which showed genome-wide significance with an annual change in eGFR. In stage 2 meta-analysis, the significant association at UMOD was replicated. Associations at GALNT11 with Rapid Decline (annual eGFR decline of 3 ml/min per 1.73 m(2) or more), and CDH23 with eGFR change among those with CKD showed significant suggestive evidence of replication. Combined stage 1 and 2 meta-analyses showed significance for UMOD, GALNT11, and CDH23. Morpholino knockdowns of galnt11 and cdh23 in zebrafish embryos each had signs of severe edema 72 h after gentamicin treatment compared with controls, but no gross morphological renal abnormalities before gentamicin administration. Thus, our results suggest a role in the deterioration of kidney function for the loci GALNT11 and CDH23, and show that the UMOD locus is significantly associated with kidney function decline.

Novel locus including FGF21 is associated with dietary macronutrient intake.

Dietary intake of macronutrients (carbohydrate, protein, and fat) has been associated with risk of chronic conditions such as obesity and diabetes. Family studies have reported a moderate contribution of genetics to variation in macronutrient intake. In a genome-wide meta-analysis of a population-based discovery cohort (n = 33 533), rs838133 in FGF21 (19q13.33), rs197273 near TRAF family member-associated NF-kappa-B activator (TANK) (2p24.2), and rs10163409 in FTO (16q12.2) were among the top associations (P < 10(-5)) for percentage of total caloric intake from protein and carbohydrate. rs838133 was replicated in silico in an independent sample from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (CHARGE) Nutrition Working Group (n = 38 360) and attained genome-wide significance in combined analysis (Pjoint = 7.9 × 10(-9)). A cytokine involved in cellular metabolism, FGF21 is a potential susceptibility gene for obesity and type 2 diabetes. Our results highlight the potential of genetic variation for determining dietary macronutrient intake.