X-chromosome and kidney function: evidence from a multi-trait genetic analysis of 908,697 individuals reveals sex-specific and sex-differential findings in genes regulated by androgen response elements

X-chromosomal genetic variants are understudied but can yield valuable insights into sexually dimorphic human traits and diseases. We performed a sex-stratified cross-ancestry X-chromosome-wide association meta-analysis of seven kidney-related traits (n = 908,697), identifying 23 loci genome-wide significantly associated with two of the traits: 7 for uric acid and 16 for estimated glomerular filtration rate (eGFR), including four novel eGFR loci containing the functionally plausible prioritized genes ACSL4, CLDN2, TSPAN6 and the female-specific DRP2. Further, we identified five novel sex-interactions, comprising male-specific effects at FAM9B and AR/EDA2R, and three sex-differential findings with larger genetic effect sizes in males at DCAF12L1 and MST4 and larger effect sizes in females at HPRT1. All prioritized genes in loci showing significant sex-interactions were located next to androgen response elements (ARE). Five ARE genes showed sex-differential expressions. This study contributes new insights into sex-dimorphisms of kidney traits along with new prioritized gene targets for further molecular research.

The Type 2 Diabetes Knowledge Portal: An open access genetic resource dedicated to type 2 diabetes and related traits

Associations between human genetic variation and clinical phenotypes have become a foundation of biomedical research. Most repositories of these data seek to be disease-agnostic and therefore lack disease-focused views. The Type 2 Diabetes Knowledge Portal (T2DKP) is a public resource of genetic datasets and genomic annotations dedicated to type 2 diabetes (T2D) and related traits. Here, we seek to make the T2DKP more accessible to prospective users and more useful to existing users. First, we evaluate the T2DKP's comprehensiveness by comparing its datasets with those of other repositories. Second, we describe how researchers unfamiliar with human genetic data can begin using and correctly interpreting them via the T2DKP. Third, we describe how existing users can extend their current workflows to use the full suite of tools offered by the T2DKP. We finally discuss the lessons offered by the T2DKP toward the goal of democratizing access to complex disease genetic results.

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.

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.

Genome-wide association analyses of physical activity and sedentary behavior provide insights into underlying mechanisms and roles in disease prevention

Although physical activity and sedentary behavior are moderately heritable, little is known about the mechanisms that influence these traits. Combining data for up to 703,901 individuals from 51 studies in a multi-ancestry meta-analysis of genome-wide association studies yields 99 loci that associate with self-reported moderate-to-vigorous intensity physical activity during leisure time (MVPA), leisure screen time (LST) and/or sedentary behavior at work. Loci associated with LST are enriched for genes whose expression in skeletal muscle is altered by resistance training. A missense variant in ACTN3 makes the alpha-actinin-3 filaments more flexible, resulting in lower maximal force in isolated type IIA muscle fibers, and possibly protection from exercise-induced muscle damage. Finally, Mendelian randomization analyses show that beneficial effects of lower LST and higher MVPA on several risk factors and diseases are mediated or confounded by body mass index (BMI). Our results provide insights into physical activity mechanisms and its role in disease prevention.

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.

Differential and shared genetic effects on kidney function between diabetic and non-diabetic individuals

Reduced glomerular filtration rate (GFR) can progress to kidney failure. Risk factors include genetics and diabetes mellitus (DM), but little is known about their interaction. We conducted genome-wide association meta-analyses for estimated GFR based on serum creatinine (eGFR), separately for individuals with or without DM (nDM = 178,691, nnoDM = 1,296,113). Our genome-wide searches identified (i) seven eGFR loci with significant DM/noDM-difference, (ii) four additional novel loci with suggestive difference and (iii) 28 further novel loci (including CUBN) by allowing for potential difference. GWAS on eGFR among DM individuals identified 2 known and 27 potentially responsible loci for diabetic kidney disease. Gene prioritization highlighted 18 genes that may inform reno-protective drug development. We highlight the existence of DM-only and noDM-only effects, which can inform about the target group, if respective genes are advanced as drug targets. Largely shared effects suggest that most drug interventions to alter eGFR should be effective in DM and noDM.

MO802: Genetics Plays a Limited Role in ESA-Hyporesponsiveness and Haemoglobin Outcomes in End-Stage Renal Disease Patients on Haemodialysis

BACKGROUND AND AIMS

Patients with end-stage renal disease on haemodialysis (HD) frequently have resistance or hyporesponsiveness (hypoR) to erythropoiesis-stimulating agents (ESAs) and are also often iron deficient. Despite evidence of increased comorbidities impacting outcomes such as treatment response and achieving targeted haemoglobin (Hb) and iron levels, there is substantial residual variability in these outcomes that remains unexplained. The ability to identify higher-risk sub-populations may provide additional opportunities to intervene. Genetic loci associated with Hb concentration have been identified in genome-wide association studies. Genetic risk scores (GRS) that provide a means to aggregate the health-related risk of multiple genetic alleles into a single number for risk assessment have been validated using these variants. We evaluated if GRS associated with Hb concentration, in addition to clinical factors, improves prediction of ESA hypoR, Hb and iron outcomes better than either one used alone.

METHOD

A collaborative, non-interventional study between GlaxoSmithKline and DaVita enrolled 3192 African American (AA) and European American (EA) HD patients. Patient demographics, laboratory and vital measurements, concomitant medications, and comorbidity records were drawn from Electronic Health Records, with patient blood samples obtained for genetic analysis. We created GRS for transethnic and euro-centric populations for Hb concentrations from the literature, and a GRS specific to EGLN1/2/3 genes. ESA hypoR status, iron replete status, and time weighted Hb concentration were evaluated as outcomes. ESA hypoR for epoetin alfa was defined by an erythropoietin resistance index ≥ 2.0 U/kg/wk/g/L [intravenous (IV)] or ≥ 1.33 U/kg/wk/g/L [subcutaneous (SC)], or dose ≥ 450 U/kg/wk (IV) or ≥ 300 U/kg/wk (SC). Iron replete was defined as a ferritin ≥ 100 ng/mL and transferrin saturation ≥ 20%. Heritability was estimated using genome complex trait analysis (GCTA). All analyses were performed individually for each population (AA and EA). Significant P-value threshold was 0.0167 (0.05/3 main effects).

RESULTS

None of the genetic risk scores tested had any clinically meaningful association with outcomes. Though the GRS specific to EGLN1/2/3 was significantly associated with Hb concentration (P &lt; .01) in the AA population, the Beta squared effect size (ES) was very low (ES = 0.005). GCTA analyses were inconclusive with moderate heritability (h2 &gt;8%) and high standard error (&gt;30%), suggesting high variability in heritability estimates for outcomes. Ignoring GRS, several baseline measurements were significantly associated with outcomes. Baseline Hb was consistently and significantly associated with ESA hypoR and Hb concentration in both AA and EA populations, while significantly associated with iron replete status in only the EA population (P &lt; .0001). Other variables such as age and baseline albumin were significantly associated with ESA hypoR and iron replete status, but only in either the EA or AA population (Table). Significant differences were also observed in baseline characteristics between the EA and AA patients. AA patients were 8 years younger and on HD for 2 years longer (P &lt; .0001) at enrolment compared with EA patients. Baseline ferritin, creatinine, albumin and hypertensive status were also significantly different at baseline (P &lt; .01).

CONCLUSION

There was limited utility for genetics in predicting ESA hypoR, haemoglobin and iron outcomes in HD patients in AA and EA populations. GCTA analyses were inconclusive, suggesting high variability in heritability estimates. Significant differences in baseline variables between the AA and EA populations that were also significantly associated with ESA hypoR, haemoglobin and iron outcomes suggested that population-level baseline differences must be considered when identifying at-risk sub-populations.

Genetic analysis of dietary intake identifies new loci and functional links with metabolic traits

Dietary intake is a major contributor to the global obesity epidemic and represents a complex behavioural phenotype that is partially affected by innate biological differences. Here, we present a multivariate genome-wide association analysis of overall variation in dietary intake to account for the correlation between dietary carbohydrate, fat and protein in 282,271 participants of European ancestry from the UK Biobank (n = 191,157) and Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (n = 91,114), and identify 26 distinct genome-wide significant loci. Dietary intake signals map exclusively to specific brain regions and are enriched for genes expressed in specialized subtypes of GABAergic, dopaminergic and glutamatergic neurons. We identified two main clusters of genetic variants for overall variation in dietary intake that were differently associated with obesity and coronary artery disease. These results enhance the biological understanding of interindividual differences in dietary intake by highlighting neural mechanisms, supporting functional follow-up experiments and possibly providing new avenues for the prevention and treatment of prevalent complex metabolic diseases.

The trans-ancestral genomic architecture of glycemic traits

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

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.

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

Objective- Inflammation is a causal risk factor for cardiovascular disease (CVD). sPLA₂-IIA (group IIA secretory phospholipase A₂) plays an integral role in regulating vascular inflammation. Although studies investigated sPLA₂-IIA in secondary prevention, we prospectively evaluated sPLA₂-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 sPLA₂-IIA mass was measured (N=11 269 baseline; N=9620 1 year). We also identified genetic variants influencing sPLA₂-IIA using genome-wide association and examined them with CVD. Three hundred thirteen incident CVD events occurred during follow-up. Baseline sPLA₂-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 sPLA₂-IIA) had the strongest effect on sPLA₂-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, sPLA₂-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 sPLA₂-IIA may be a potential pharmacological target for primary prevention of CVD. Clinical Trial Registration- URL: http://www.clinicaltrials.gov . Unique identifier: NCT00239681.

Epigenome-Wide Association Study of Incident Type 2 Diabetes in a British Population: EPIC-Norfolk Study

Epigenetic changes may contribute substantially to risks of diseases of aging. Previous studies reported seven methylation variable positions (MVPs) robustly associated with incident type 2 diabetes mellitus (T2DM). However, their causal roles in T2DM are unclear. In an incident T2DM case-cohort study nested within the population-based European Prospective Investigation into Cancer and Nutrition (EPIC)-Norfolk cohort, we used whole blood DNA collected at baseline, up to 11 years before T2DM onset, to investigate the role of methylation in the etiology of T2DM. We identified 15 novel MVPs with robust associations with incident T2DM and robustly confirmed three MVPs identified previously (near to TXNIP, ABCG1, and SREBF1). All 18 MVPs showed directionally consistent associations with incident and prevalent T2DM in independent studies. Further conditional analyses suggested that the identified epigenetic signals appear related to T2DM via glucose and obesity-related pathways acting before the collection of baseline samples. We integrated genome-wide genetic data to identify methylation-associated quantitative trait loci robustly associated with 16 of the 18 MVPs and found one MVP, cg00574958 at CPT1A, with a possible direct causal role in T2DM. None of the implicated genes were previously highlighted by genetic association studies, suggesting that DNA methylation studies may reveal novel biological mechanisms involved in tissue responses to glycemia.

Target genes, variants, tissues and transcriptional pathways influencing human serum urate levels

Elevated serum urate levels cause gout and correlate with cardiometabolic diseases via poorly understood mechanisms. We performed a trans-ancestry genome-wide association study of serum urate in 457,690 individuals, identifying 183 loci (147 previously unknown) that improve the prediction of gout in an independent cohort of 334,880 individuals. Serum urate showed significant genetic correlations with many cardiometabolic traits, with genetic causality analyses supporting a substantial role for pleiotropy. Enrichment analysis, fine-mapping of urate-associated loci and colocalization with gene expression in 47 tissues implicated the kidney and liver as the main target organs and prioritized potentially causal genes and variants, including the transcriptional master regulators in the liver and kidney, HNF1A and HNF4A. Experimental validation showed that HNF4A transactivated the promoter of ABCG2, encoding a major urate transporter, in kidney cells, and that HNF4A p.Thr139Ile is a functional variant. Transcriptional coregulation within and across organs may be a general mechanism underlying the observed pleiotropy between urate and cardiometabolic traits.

Genome-wide association meta-analyses and fine-mapping elucidate pathways influencing albuminuria

Increased levels of the urinary albumin-to-creatinine ratio (UACR) are associated with higher risk of kidney disease progression and cardiovascular events, but underlying mechanisms are incompletely understood. Here, we conduct trans-ethnic (n = 564,257) and European-ancestry specific meta-analyses of genome-wide association studies of UACR, including ancestry- and diabetes-specific analyses, and identify 68 UACR-associated loci. Genetic correlation analyses and risk score associations in an independent electronic medical records database (n = 192,868) reveal connections with proteinuria, hyperlipidemia, gout, and hypertension. Fine-mapping and trans-Omics analyses with gene expression in 47 tissues and plasma protein levels implicate genes potentially operating through differential expression in kidney (including TGFB1, MUC1, PRKCI, and OAF), and allow coupling of UACR associations to altered plasma OAF concentrations. Knockdown of OAF and PRKCI orthologs in Drosophila nephrocytes reduces albumin endocytosis. Silencing fly PRKCI further impairs slit diaphragm formation. These results generate a priority list of genes and pathways for translational research to reduce albuminuria.

Quality of dietary fat and genetic risk of type 2 diabetes: individual participant data meta-analysis

OBJECTIVE

To investigate whether the genetic burden of type 2 diabetes modifies the association between the quality of dietary fat and the incidence of type 2 diabetes.

DESIGN

Individual participant data meta-analysis.

DATA SOURCES

Eligible prospective cohort studies were systematically sourced from studies published between January 1970 and February 2017 through electronic searches in major medical databases (Medline, Embase, and Scopus) and discussion with investigators.

REVIEW METHODS

Data from cohort studies or multicohort consortia with available genome-wide genetic data and information about the quality of dietary fat and the incidence of type 2 diabetes in participants of European descent was sought. Prospective cohorts that had accrued five or more years of follow-up were included. The type 2 diabetes genetic risk profile was characterized by a 68-variant polygenic risk score weighted by published effect sizes. Diet was recorded by using validated cohort-specific dietary assessment tools. Outcome measures were summary adjusted hazard ratios of incident type 2 diabetes for polygenic risk score, isocaloric replacement of carbohydrate (refined starch and sugars) with types of fat, and the interaction of types of fat with polygenic risk score.

RESULTS

Of 102 305 participants from 15 prospective cohort studies, 20 015 type 2 diabetes cases were documented after a median follow-up of 12 years (interquartile range 9.4-14.2). The hazard ratio of type 2 diabetes per increment of 10 risk alleles in the polygenic risk score was 1.64 (95% confidence interval 1.54 to 1.75, I²=7.1%, τ²=0.003). The increase of polyunsaturated fat and total omega 6 polyunsaturated fat intake in place of carbohydrate was associated with a lower risk of type 2 diabetes, with hazard ratios of 0.90 (0.82 to 0.98, I²=18.0%, τ²=0.006; per 5% of energy) and 0.99 (0.97 to 1.00, I²=58.8%, τ²=0.001; per increment of 1 g/d), respectively. Increasing monounsaturated fat in place of carbohydrate was associated with a higher risk of type 2 diabetes (hazard ratio 1.10, 95% confidence interval 1.01 to 1.19, I²=25.9%, τ²=0.006; per 5% of energy). Evidence of small study effects was detected for the overall association of polyunsaturated fat with the risk of type 2 diabetes, but not for the omega 6 polyunsaturated fat and monounsaturated fat associations. Significant interactions between dietary fat and polygenic risk score on the risk of type 2 diabetes (P>0.05 for interaction) were not observed.

CONCLUSIONS

These data indicate that genetic burden and the quality of dietary fat are each associated with the incidence of type 2 diabetes. The findings do not support tailoring recommendations on the quality of dietary fat to individual type 2 diabetes genetic risk profiles for the primary prevention of type 2 diabetes, and suggest that dietary fat is associated with the risk of type 2 diabetes across the spectrum of type 2 diabetes genetic risk.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

An integrative cross-omics analysis of DNA methylation sites of glucose and insulin homeostasis

Despite existing reports on differential DNA methylation in type 2 diabetes (T2D) and obesity, our understanding of its functional relevance remains limited. Here we show the effect of differential methylation in the early phases of T2D pathology by a blood-based epigenome-wide association study of 4808 non-diabetic Europeans in the discovery phase and 11,750 individuals in the replication. We identify CpGs in LETM1, RBM20, IRS2, MAN2A2 and the 1q25.3 region associated with fasting insulin, and in FCRL6, SLAMF1, APOBEC3H and the 15q26.1 region with fasting glucose. In silico cross-omics analyses highlight the role of differential methylation in the crosstalk between the adaptive immune system and glucose homeostasis. The differential methylation explains at least 16.9% of the association between obesity and insulin. Our study sheds light on the biological interactions between genetic variants driving differential methylation and gene expression in the early pathogenesis of T2D.

A catalog of genetic loci associated with kidney function from analyses of a million individuals

Chronic kidney disease (CKD) is responsible for a public health burden with multi-systemic complications. Through trans-ancestry meta-analysis of genome-wide association studies of estimated glomerular filtration rate (eGFR) and independent replication (n = 1,046,070), we identified 264 associated loci (166 new). Of these, 147 were likely to be relevant for kidney function on the basis of associations with the alternative kidney function marker blood urea nitrogen (n = 416,178). Pathway and enrichment analyses, including mouse models with renal phenotypes, support the kidney as the main target organ. A genetic risk score for lower eGFR was associated with clinically diagnosed CKD in 452,264 independent individuals. Colocalization analyses of associations with eGFR among 783,978 European-ancestry individuals and gene expression across 46 human tissues, including tubulo-interstitial and glomerular kidney compartments, identified 17 genes differentially expressed in kidney. Fine-mapping highlighted missense driver variants in 11 genes and kidney-specific regulatory variants. These results provide a comprehensive priority list of molecular targets for translational research.

Protein-coding variants implicate novel genes related to lipid homeostasis contributing to body-fat distribution

Body-fat distribution is a risk factor for adverse cardiovascular health consequences. We analyzed the association of body-fat distribution, assessed by waist-to-hip ratio adjusted for body mass index, with 228,985 predicted coding and splice site variants available on exome arrays in up to 344,369 individuals from five major ancestries (discovery) and 132,177 European-ancestry individuals (validation). We identified 15 common (minor allele frequency, MAF ≥5%) and nine low-frequency or rare (MAF <5%) coding novel variants. Pathway/gene set enrichment analyses identified lipid particle, adiponectin, abnormal white adipose tissue physiology and bone development and morphology as important contributors to fat distribution, while cross-trait associations highlight cardiometabolic traits. In functional follow-up analyses, specifically in Drosophila RNAi-knockdowns, we observed a significant increase in the total body triglyceride levels for two genes (DNAH10 and PLXND1). We implicate novel genes in fat distribution, stressing the importance of interrogating low-frequency and protein-coding variants.

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.

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 4,303), 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 and are cancer related.

Conclusion:

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.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the version of this article originally published, one of the two authors with the name Wei Zhao was omitted from the author list and the affiliations for both authors were assigned to the single Wei Zhao in the author list. In addition, the ORCID for Wei Zhao (Department of Biostatistics and Epidemiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA) was incorrectly assigned to author Wei Zhou. The errors have been corrected in the HTML and PDF versions of the article.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the published version of this paper, the name of author Emanuele Di Angelantonio was misspelled. This error has now been corrected in the HTML and PDF versions of the article.

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.

Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity.

Genome-wide association meta-analysis of fish and EPA+DHA consumption in 17 US and European cohorts

BACKGROUND

Regular fish and omega-3 consumption may have several health benefits and are recommended by major dietary guidelines. Yet, their intakes remain remarkably variable both within and across populations, which could partly owe to genetic influences.

OBJECTIVE

To identify common genetic variants that influence fish and dietary eicosapentaenoic acid plus docosahexaenoic acid (EPA+DHA) consumption.

DESIGN

We conducted genome-wide association (GWA) meta-analysis of fish (n = 86,467) and EPA+DHA (n = 62,265) consumption in 17 cohorts of European descent from the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) Consortium Nutrition Working Group. Results from cohort-specific GWA analyses (additive model) for fish and EPA+DHA consumption were adjusted for age, sex, energy intake, and population stratification, and meta-analyzed separately using fixed-effect meta-analysis with inverse variance weights (METAL software). Additionally, heritability was estimated in 2 cohorts.

RESULTS

Heritability estimates for fish and EPA+DHA consumption ranged from 0.13-0.24 and 0.12-0.22, respectively. A significant GWA for fish intake was observed for rs9502823 on chromosome 6: each copy of the minor allele (FreqA = 0.015) was associated with 0.029 servings/day (~1 serving/month) lower fish consumption (P = 1.96x10-8). No significant association was observed for EPA+DHA, although rs7206790 in the obesity-associated FTO gene was among top hits (P = 8.18x10-7). Post-hoc calculations demonstrated 95% statistical power to detect a genetic variant associated with effect size of 0.05% for fish and 0.08% for EPA+DHA.

CONCLUSIONS

These novel findings suggest that non-genetic personal and environmental factors are principal determinants of the remarkable variation in fish consumption, representing modifiable targets for increasing intakes among all individuals. Genes underlying the signal at rs72838923 and mechanisms for the association warrant further investigation.

Exome-wide association study of plasma lipids in >300,000 individuals

We screened variants on an exome-focused genotyping array in >300,000 participants (replication in >280,000 participants) and identified 444 independent variants in 250 loci significantly associated with total cholesterol (TC), high-density-lipoprotein cholesterol (HDL-C), low-density-lipoprotein cholesterol (LDL-C), and/or triglycerides (TG). At two loci (JAK2 and A1CF), experimental analysis in mice showed lipid changes consistent with the human data. We also found that: (i) beta-thalassemia trait carriers displayed lower TC and were protected from coronary artery disease (CAD); (ii) excluding the CETP locus, there was not a predictable relationship between plasma HDL-C and risk for age-related macular degeneration; (iii) only some mechanisms of lowering LDL-C appeared to increase risk for type 2 diabetes (T2D); and (iv) TG-lowering alleles involved in hepatic production of TG-rich lipoproteins (TM6SF2 and PNPLA3) tracked with higher liver fat, higher risk for T2D, and lower risk for CAD, whereas TG-lowering alleles involved in peripheral lipolysis (LPL and ANGPTL4) had no effect on liver fat but decreased risks for both T2D and CAD.

Epigenome-wide association studies identify DNA methylation associated with kidney function

Chronic kidney disease (CKD) is defined by reduced estimated glomerular filtration rate (eGFR). Previous genetic studies have implicated regulatory mechanisms contributing to CKD. Here we present epigenome-wide association studies of eGFR and CKD using whole-blood DNA methylation of 2264 ARIC Study and 2595 Framingham Heart Study participants to identify epigenetic signatures of kidney function. Of 19 CpG sites significantly associated (P < 1e-07) with eGFR/CKD and replicated, five also associate with renal fibrosis in biopsies from CKD patients and show concordant DNA methylation changes in kidney cortex. Lead CpGs at PTPN6/PHB2, ANKRD11, and TNRC18 map to active enhancers in kidney cortex. At PTPN6/PHB2 cg19942083, methylation in kidney cortex associates with lower renal PTPN6 expression, higher eGFR, and less renal fibrosis. The regions containing the 243 eGFR-associated (P < 1e-05) CpGs are significantly enriched for transcription factor binding sites of EBF1, EP300, and CEBPB (P < 5e-6). Our findings highlight kidney function associated epigenetic variation.

Genome-wide association studies of kidney function show enrichment of associated genetic variants in regulatory regions. Here, the authors perform epigenome-wide association studies of kidney function and disease, identifying 19 CpG sites significantly associated with these.

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.

Atherogenic Lipoprotein Determinants of Cardiovascular Disease and Residual Risk Among Individuals With Low Low-Density Lipoprotein Cholesterol

Background

Levels of LDL (low‐density lipoprotein) cholesterol in the population are declining, and increasing attention is being focused on residual lipid‐related pathways of atherosclerotic cardiovascular disease risk beyond LDL cholesterol. Among individuals with low (<130 mg/dL) LDL cholesterol, we undertook detailed profiling of circulating atherogenic lipoproteins in relation to incident cardiovascular disease in 2 populations.

Methods and Results

We performed proton nuclear magnetic resonance spectroscopy to quantify concentrations of LDL and VLDL (very low‐density lipoprotein) particle subclasses in 11 984 JUPITER trial participants (NCT00239681). Adjusted Cox models examined cardiovascular disease risk associated with lipoprotein measures according to treatment allocation. Risk (adjusted hazard ratio [95%CI] per SD increment) among placebo‐allocated participants was associated with total LDL particles (1.19 [1.02, 1.38]) and total VLDL particles (1.21 [1.04, 1.41]), as well as apolipoprotein B, non–high‐density lipoprotein cholesterol, and triglycerides, but not LDL‐c. Rosuvastatin reduced LDL measures but had variable effects on triglyceride and VLDL measures. On‐statin levels of the smallest VLDL particle subclass were associated with a 68% per‐SD (adjusted hazard ratio 1.68 [1.28, 2.22]) increase in residual risk—this risk was related to VLDL cholesterol and not triglyceride or larger VLDL particles. There was evidence that residual risk prediction during statin therapy could be significantly improved through the inclusion of key VLDL measures (Harrell C‐index 0.780 versus 0.712; P<0.0001). In an independent, prospective cohort of 4721 individuals referred for cardiac catheterization (CATHGEN), similar patterns of lipoprotein‐related risk were observed.

Conclusions

Atherogenic lipoprotein particle concentrations were associated with cardiovascular disease risk when LDL cholesterol was low. VLDL lipoproteins, particularly the smallest remnant subclass, may represent unused targets for risk prediction and potential therapeutic intervention for reducing residual risk.

Clinical Trial Registration

URL: http://www.clinicaltrials.gov. Unique identifier: NCT00239681.

A Low-Frequency Inactivating AKT2 Variant Enriched in the Finnish Population Is Associated With Fasting Insulin Levels and Type 2 Diabetes Risk

To identify novel coding association signals and facilitate characterization of mechanisms influencing glycemic traits and type 2 diabetes risk, we analyzed 109,215 variants derived from exome array genotyping together with an additional 390,225 variants from exome sequence in up to 39,339 normoglycemic individuals from five ancestry groups. We identified a novel association between the coding variant (p.Pro50Thr) in AKT2 and fasting plasma insulin (FI), a gene in which rare fully penetrant mutations are causal for monogenic glycemic disorders. The low-frequency allele is associated with a 12% increase in FI levels. This variant is present at 1.1% frequency in Finns but virtually absent in individuals from other ancestries. Carriers of the FI-increasing allele had increased 2-h insulin values, decreased insulin sensitivity, and increased risk of type 2 diabetes (odds ratio 1.05). In cellular studies, the AKT2-Thr50 protein exhibited a partial loss of function. We extend the allelic spectrum for coding variants in AKT2 associated with disorders of glucose homeostasis and demonstrate bidirectional effects of variants within the pleckstrin homology domain of AKT2.

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 (P_v), G×E interaction effects (with smoking and physical activity), and marginal genetic effects (P_m). Correlations between P_v and P_m were stronger for SNPs with established marginal effects (Spearman’s ρ = 0.401 for triglycerides, and ρ = 0.236 for BMI) compared to all SNPs. When P_v and P_m 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 P_v distribution (P_binomial <0.05). SNPs from the top 1% of the P_m distribution for BMI had more significant P_v values (P_{Mann–Whitney}= 1.46×10⁻⁵), and the odds ratio of SNPs with nominally significant (<0.05) P_m and P_v was 1.33 (95% CI: 1.12, 1.57) for BMI. Moreover, BMI SNPs with nominally significant G×E interaction P-values (P_int<0.05) were enriched with nominally significant P_v values (P_binomial = 8.63×10⁻⁹ and 8.52×10⁻⁷ 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.

Most contemporary studies of gene-environment interactions focus on gene variants that are known to bear strong and reliable associations with the traits of interest. The strategy is intuitive because it helps limit the number of tests performed by focusing on a relatively small number of gene variants. However, this approach is predicated on an implicit assumption that these loci are strong candidates for interactions owing to their established relationships with the index traits. The counter-argument is that, because these loci have highly consistent signals within and between populations that vary by environmental characteristics, the probability that these variants interact with other factors is low. The current analysis tests whether variants with strong marginal effects signals (i.e., those prioritized through conventional genome-wide association analyses) are strong or weak candidates for gene-environment interactions. Here we describe analyses focused on lipids and BMI that test this hypothesis by comparing marginal effect signals with variance effect signals and those derived from explicit genome-wide, gene-environment interaction analyses. We conclude that for BMI, there are features of the top-ranking marginal effect loci that render them stronger candidates for interactions than is true of variants with weaker marginal effects signals. These findings are likely to help optimize the efficiency of future gene-environment interaction analyses by providing evidence-based rankings for strong candidate loci.

Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution.

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.

Rare and low-frequency coding variants alter human adult height

Height is a highly heritable, classic polygenic trait with approximately 700 common associated variants identified through genome-wide association studies so far. Here, we report 83 height-associated coding variants with lower minor-allele frequencies (in the range of 0.1-4.8%) and effects of up to 2 centimetres per allele (such as those in IHH, STC2, AR and CRISPLD2), greater than ten times the average effect of common variants. In functional follow-up studies, rare height-increasing alleles of STC2 (giving an increase of 1-2 centimetres per allele) compromised proteolytic inhibition of PAPP-A and increased cleavage of IGFBP-4 in vitro, resulting in higher bioavailability of insulin-like growth factors. These 83 height-associated variants overlap genes that are mutated in monogenic growth disorders and highlight new biological candidates (such as ADAMTS3, IL11RA and NOX4) and pathways (such as proteoglycan and glycosaminoglycan synthesis) involved in growth. Our results demonstrate that sufficiently large sample sizes can uncover rare and low-frequency variants of moderate-to-large effect associated with polygenic human phenotypes, and that these variants implicate relevant genes and pathways.

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.

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.

DNA methylation signatures of chronic low-grade inflammation are associated with complex diseases

BACKGROUND

Chronic low-grade inflammation reflects a subclinical immune response implicated in the pathogenesis of complex diseases. Identifying genetic loci where DNA methylation is associated with chronic low-grade inflammation may reveal novel pathways or therapeutic targets for inflammation.

RESULTS

We performed a meta-analysis of epigenome-wide association studies (EWAS) of serum C-reactive protein (CRP), which is a sensitive marker of low-grade inflammation, in a large European population (n = 8863) and trans-ethnic replication in African Americans (n = 4111). We found differential methylation at 218 CpG sites to be associated with CRP (P < 1.15 × 10^-7) in the discovery panel of European ancestry and replicated (P < 2.29 × 10^-4) 58 CpG sites (45 unique loci) among African Americans. To further characterize the molecular and clinical relevance of the findings, we examined the association with gene expression, genetic sequence variants, and clinical outcomes. DNA methylation at nine (16%) CpG sites was associated with whole blood gene expression in cis (P < 8.47 × 10^-5), ten (17%) CpG sites were associated with a nearby genetic variant (P < 2.50 × 10^-3), and 51 (88%) were also associated with at least one related cardiometabolic entity (P < 9.58 × 10^-5). An additive weighted score of replicated CpG sites accounted for up to 6% inter-individual variation (R2) of age-adjusted and sex-adjusted CRP, independent of known CRP-related genetic variants.

CONCLUSION

We have completed an EWAS of chronic low-grade inflammation and identified many novel genetic loci underlying inflammation that may serve as targets for the development of novel therapeutic interventions for inflammation.

SOS2 and ACP1 Loci Identified through Large-Scale Exome Chip Analysis Regulate Kidney Development and Function

Genome-wide association studies have identified >50 common variants associated with kidney function, but these variants do not fully explain the variation in eGFR. We performed a two-stage meta-analysis of associations between genotypes from the Illumina exome array and eGFR on the basis of serum creatinine (eGFRcrea) among participants of European ancestry from the CKDGen Consortium (n_Stage1: 111,666; n_Stage2: 48,343). In single-variant analyses, we identified single nucleotide polymorphisms at seven new loci associated with eGFRcrea (PPM1J, EDEM3, ACP1, SPEG, EYA4, CYP1A1, and ATXN2L; P_Stage1<3.7×10^-7), of which most were common and annotated as nonsynonymous variants. Gene-based analysis identified associations of functional rare variants in three genes with eGFRcrea, including a novel association with the SOS Ras/Rho guanine nucleotide exchange factor 2 gene, SOS2 (P=5.4×10^-8 by sequence kernel association test). Experimental follow-up in zebrafish embryos revealed changes in glomerular gene expression and renal tubule morphology in the embryonic kidney of acp1- and sos2-knockdowns. These developmental abnormalities associated with altered blood clearance rate and heightened prevalence of edema. This study expands the number of loci associated with kidney function and identifies novel genes with potential roles in kidney formation.

A principal component meta-analysis on multiple anthropometric traits identifies novel loci for body shape

Large consortia have revealed hundreds of genetic loci associated with anthropometric traits, one trait at a time. We examined whether genetic variants affect body shape as a composite phenotype that is represented by a combination of anthropometric traits. We developed an approach that calculates averaged PCs (AvPCs) representing body shape derived from six anthropometric traits (body mass index, height, weight, waist and hip circumference, waist-to-hip ratio). The first four AvPCs explain >99% of the variability, are heritable, and associate with cardiometabolic outcomes. We performed genome-wide association analyses for each body shape composite phenotype across 65 studies and meta-analysed summary statistics. We identify six novel loci: LEMD2 and CD47 for AvPC1, RPS6KA5/C14orf159 and GANAB for AvPC3, and ARL15 and ANP32 for AvPC4. Our findings highlight the value of using multiple traits to define complex phenotypes for discovery, which are not captured by single-trait analyses, and may shed light onto new pathways.

Meta-analysis identifies common and rare variants influencing blood pressure and overlapping with metabolic trait loci

Meta-analyses of association results for blood pressure using exome-centric single-variant and gene-based tests identified 31 new loci in a discovery stage among 146,562 individuals, with follow-up and meta-analysis in 180,726 additional individuals (total n = 327,288). These blood pressure-associated loci are enriched for known variants for cardiometabolic traits. Associations were also observed for the aggregation of rare and low-frequency missense variants in three genes, NPR1, DBH, and PTPMT1. In addition, blood pressure associations at 39 previously reported loci were confirmed. The identified variants implicate biological pathways related to cardiometabolic traits, vascular function, and development. Several new variants are inferred to have roles in transcription or as hubs in protein-protein interaction networks. Genetic risk scores constructed from the identified variants were strongly associated with coronary disease and myocardial infarction. This large collection of blood pressure-associated loci suggests new therapeutic strategies for hypertension, emphasizing a link with cardiometabolic risk.

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.73m².

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 genomic approach to therapeutic target validation identifies a glucose-lowering GLP1R variant protective for coronary heart disease

Regulatory authorities have indicated that new drugs to treat type 2 diabetes (T2D) should not be associated with an unacceptable increase in cardiovascular risk. Human genetics may be able to guide development of antidiabetic therapies by predicting cardiovascular and other health endpoints. We therefore investigated the association of variants in six genes that encode drug targets for obesity or T2D with a range of metabolic traits in up to 11,806 individuals by targeted exome sequencing and follow-up in 39,979 individuals by targeted genotyping, with additional in silico follow-up in consortia. We used these data to first compare associations of variants in genes encoding drug targets with the effects of pharmacological manipulation of those targets in clinical trials. We then tested the association of those variants with disease outcomes, including coronary heart disease, to predict cardiovascular safety of these agents. A low-frequency missense variant (Ala316Thr; rs10305492) in the gene encoding glucagon-like peptide-1 receptor (GLP1R), the target of GLP1R agonists, was associated with lower fasting glucose and T2D risk, consistent with GLP1R agonist therapies. The minor allele was also associated with protection against heart disease, thus providing evidence that GLP1R agonists are not likely to be associated with an unacceptable increase in cardiovascular risk. Our results provide an encouraging signal that these agents may be associated with benefit, a question currently being addressed in randomized controlled trials. Genetic variants associated with metabolic traits and multiple disease outcomes can be used to validate therapeutic targets at an early stage in the drug development process.

Genome-wide Association Studies Identify Genetic Loci Associated With Albuminuria in Diabetes

Elevated concentrations of albumin in the urine, albuminuria, are a hallmark of diabetic kidney disease and are associated with an increased risk for end-stage renal disease and cardiovascular events. To gain insight into the pathophysiological mechanisms underlying albuminuria, we conducted meta-analyses of genome-wide association studies and independent replication in up to 5,825 individuals of European ancestry with diabetes and up to 46,061 without diabetes, followed by functional studies. Known associations of variants in CUBN, encoding cubilin, with the urinary albumin-to-creatinine ratio (UACR) were confirmed in the overall sample (P = 2.4 × 10(-10)). Gene-by-diabetes interactions were detected and confirmed for variants in HS6ST1 and near RAB38/CTSC. Single nucleotide polymorphisms at these loci demonstrated a genetic effect on UACR in individuals with but not without diabetes. The change in the average UACR per minor allele was 21% for HS6ST1 (P = 6.3 × 10(-7)) and 13% for RAB38/CTSC (P = 5.8 × 10(-7)). Experiments using streptozotocin-induced diabetic Rab38 knockout and control rats showed higher urinary albumin concentrations and reduced amounts of megalin and cubilin at the proximal tubule cell surface in Rab38 knockout versus control rats. Relative expression of RAB38 was higher in tubuli of patients with diabetic kidney disease compared with control subjects. The loci identified here confirm known pathways and highlight novel pathways influencing albuminuria.

Interaction of methylation-related genetic variants with circulating fatty acids on plasma lipids: a meta-analysis of 7 studies and methylation analysis of 3 studies in the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium

BACKGROUND

DNA methylation is influenced by diet and single nucleotide polymorphisms (SNPs), and methylation modulates gene expression.

OBJECTIVE

We aimed to explore whether the gene-by-diet interactions on blood lipids act through DNA methylation.

DESIGN

We selected 7 SNPs on the basis of predicted relations in fatty acids, methylation, and lipids. We conducted a meta-analysis and a methylation and mediation analysis with the use of data from the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) consortium and the ENCODE (Encyclopedia of DNA Elements) consortium.

RESULTS

On the basis of the meta-analysis of 7 cohorts in the CHARGE consortium, higher plasma HDL cholesterol was associated with fewer C alleles at ATP-binding cassette subfamily A member 1 (ABCA1) rs2246293 (β = -0.6 mg/dL, P = 0.015) and higher circulating eicosapentaenoic acid (EPA) (β = 3.87 mg/dL, P = 5.62 × 10(21)). The difference in HDL cholesterol associated with higher circulating EPA was dependent on genotypes at rs2246293, and it was greater for each additional C allele (β = 1.69 mg/dL, P = 0.006). In the GOLDN (Genetics of Lipid Lowering Drugs and Diet Network) study, higher ABCA1 promoter cg14019050 methylation was associated with more C alleles at rs2246293 (β = 8.84%, P = 3.51 × 10(18)) and lower circulating EPA (β = -1.46%, P = 0.009), and the mean difference in methylation of cg14019050 that was associated with higher EPA was smaller with each additional C allele of rs2246293 (β = -2.83%, P = 0.007). Higher ABCA1 cg14019050 methylation was correlated with lower ABCA1 expression (r = -0.61, P = 0.009) in the ENCODE consortium and lower plasma HDL cholesterol in the GOLDN study (r = -0.12, P = 0.0002). An additional mediation analysis was meta-analyzed across the GOLDN study, Cardiovascular Health Study, and the Multi-Ethnic Study of Atherosclerosis. Compared with the model without the adjustment of cg14019050 methylation, the model with such adjustment provided smaller estimates of the mean plasma HDL cholesterol concentration in association with both the rs2246293 C allele and EPA and a smaller difference by rs2246293 genotypes in the EPA-associated HDL cholesterol. However, the differences between 2 nested models were NS (P > 0.05).

CONCLUSION

We obtained little evidence that the gene-by-fatty acid interactions on blood lipids act through DNA methylation.

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.

Recent Positive Selection Drives the Expansion of a Schizophrenia Risk Nonsynonymous Variant at SLC39A8 in Europeans

Natural selection has played important roles in optimizing complex human adaptations. However, schizophrenia poses an evolutionary paradox during human evolution, as the illness has strongly negative effects on fitness, but persists with a prevalence of ~0.5% across global populations. Recent studies have identified numerous risk variations in diverse populations, which might be able to explain the stable and high rate of schizophrenia morbidity in different cultures and regions, but the questions about why the risk alleles derived and maintained in human gene pool still remain unsolved. Here, we studied the evolutionary pattern of a schizophrenia risk variant rs13107325 (P < 5.0 × 10(-8) in Europeans) in the SLC39A8 gene. We found the SNP is monomorphic in Asians and Africans with risk (derived) T-allele totally absent, and further evolutionary analyses showed the T-allele has experienced recent positive selection in Europeans. Subsequent exploratory analyses implicated that the colder environment in Europe was the likely selective pressures, ie, when modern humans migrated "out of Africa" and moved to Europe mainland (a colder and cooler continent than Africa), new alleles derived due to positive selection and protected humans from risk of hypertension and also helped them adapt to the cold environment. The hypothesis was supported by our pleiotropic analyses with hypertension and energy intake as well as obesity in Europeans. Our data thus provides an intriguing example to illustrate a possible mechanism for maintaining schizophrenia risk alleles in the human gene pool, and further supported that schizophrenia is likely a product caused by pleiotropic effect during human evolution.

Differential Genetic Effects on Statin-Induced Changes Across Low-Density Lipoprotein-Related Measures

BACKGROUND

Statin therapy influences not only low-density lipoprotein (LDL) cholesterol levels but also LDL-related biomarkers, including non-high-density lipoprotein cholesterol (non-HDL-C), apolipoprotein B, total number of LDL particles, and mean LDL particle size. Recent studies have identified many genetic loci influencing circulating lipid levels and statin-induced LDL cholesterol reduction. However, it is unknown how these genetic variants influence statin-induced changes in LDL subfractions and non-HDL-C.

METHODS AND RESULTS

One hundred sixty candidate single-nucleotide polymorphisms for effects on circulating lipid levels or statin-induced LDL-cholesterol lowering were tested for association with response of LDL subfractions and non-HDL-C to rosuvastatin or placebo for 1 year among 7046 participants from the Justification for Use of Statins in Prevention: an Intervention Trial Evaluating Rosuvastatin (JUPITER) trial. Of the 51 single-nucleotide polymorphisms associated with statin response for ≥ 1 of the LDL subfractions or non-HDL-C, 20 single-nucleotide polymorphisms could be clustered according to effects predominantly on LDL particle size, predominantly on LDL particle number, and on apolipoprotein B but not on LDL cholesterol or non-HDL-C.

CONCLUSIONS

These differential associations point to pathways of LDL response to statin therapy and possibly to mechanisms of statin-dependent cardiovascular disease risk reduction.

CLINICAL TRIAL REGISTRATION

URL: http://www.clinicaltrials.gov. Unique identifier: NCT00239681.

Dietary fatty acids modulate associations between genetic variants and circulating fatty acids in plasma and erythrocyte membranes: Meta-analysis of nine studies in the CHARGE consortium

SCOPE

Tissue concentrations of omega-3 fatty acids may reduce cardiovascular disease risk, and genetic variants are associated with circulating fatty acids concentrations. Whether dietary fatty acids interact with genetic variants to modify circulating omega-3 fatty acids is unclear. We evaluated interactions between genetic variants and fatty acid intakes for circulating alpha-linoleic acid, eicosapentaenoic acid, docosahexaenoic acid, and docosapentaenoic acid.

METHODS AND RESULTS

We conducted meta-analyses (N = 11 668) evaluating interactions between dietary fatty acids and genetic variants (rs174538 and rs174548 in FADS1 (fatty acid desaturase 1), rs7435 in AGPAT3 (1-acyl-sn-glycerol-3-phosphate), rs4985167 in PDXDC1 (pyridoxal-dependent decarboxylase domain-containing 1), rs780094 in GCKR (glucokinase regulatory protein), and rs3734398 in ELOVL2 (fatty acid elongase 2)). Stratification by measurement compartment (plasma versus erthyrocyte) revealed compartment-specific interactions between FADS1 rs174538 and rs174548 and dietary alpha-linolenic acid and linoleic acid for docosahexaenoic acid and docosapentaenoic acid.

CONCLUSION

Our findings reinforce earlier reports that genetically based differences in circulating fatty acids may be partially due to differences in the conversion of fatty acid precursors. Further, fatty acids measurement compartment may modify gene-diet relationships, and considering compartment may improve the detection of gene-fatty acids interactions for circulating fatty acid outcomes.