A pseudolikelihood approach for assessing genetic association in case-control studies with unmeasured population structure

The case-control study design is one of the main tools for detecting associations between genetic markers and diseases. It is well known that population substructure can lead to spurious association between disease status and a genetic marker if the prevalence of disease and the marker allele frequency vary across subpopulations. In this paper, we propose a novel statistical method to estimate the association in case-control studies with unmeasured population substructure. The proposed method takes two steps. First, the information on genomic markers and disease status is used to infer the population substructure; second, the association between the disease and the test marker adjusting for the population substructure is modeled and estimated parametrically through polytomous logistic regression. The performance of the proposed method, relative to the existing methods, on bias, coverage probability and computational time, is assessed through simulations. The method is applied to an end-stage renal disease study in African Americans population.

Associations of serum and dialysate electrolytes with QT interval and prolongation in incident hemodialysis: the Predictors of Arrhythmic and Cardiovascular Risk in End-Stage Renal Disease (PACE) study

Background

Prolonged QT interval in hemodialysis patients may be associated with sudden cardiac death, however, few studies examined the longitudinal associations of modifiable factors such as serum and dialysate concentrations of calcium, potassium, and magnesium with corrected QT (QTc) prolongation in incident hemodialysis patients.

Methods

In 330 in-center hemodialysis participants from the PACE study who were followed up for one year, we examined the associations of predialysis serum electrolytes (total calcium [Ca], corrected Ca [cCa], ionized Ca [iCa], potassium [K], magnesium [Mg]), dialysate (dCa and dK), and serum-to-dialysate gradient measures with QTc interval and prolongation (≥460 ms in women and ≥ 450 ms in men).

Results

At the first study visit, 47% had QTc prolongation. Lower iCa and K were associated with longer QTc interval independent of potential confounders (QTc difference = 8.55[95% CI: 2.13, 14.97] ms for iCa; QTc difference = 9.89[1.58, 18.20] ms for K). Lower iCa was also associated with a higher risk of QTc prolongation. At 1 year of follow-up, 31% had persistent QTc prolongation. In longitudinal analyses, the associations of iCa and K with QTc interval remained significant, and lower K was associated with a higher risk of QTc prolongation while the association of iCa with QTc prolongation was borderline statistically significant. Serum Mg, dCa or dK, and respective gradients were not associated with QTc interval or prolongation.

Conclusion

Prolonged QTc is very common in incident hemodialysis participants and persists over follow-up. Ionized Ca and K are consistently inversely associated with QTc prolongation, which suggests closer monitoring for a low calcium or potassium level to mitigate risk.

Electronic supplementary material

The online version of this article (10.1186/s12882-019-1282-5) contains supplementary material, which is available to authorized users.

Trans-ethnic Meta-analysis and Functional Annotation Illuminates the Genetic Architecture of Fasting Glucose and Insulin

Knowledge of the genetic basis of the type 2 diabetes (T2D)-related quantitative traits fasting glucose (FG) and insulin (FI) in African ancestry (AA) individuals has been limited. In non-diabetic subjects of AA (n = 20,209) and European ancestry (EA; n = 57,292), we performed trans-ethnic (AA+EA) fine-mapping of 54 established EA FG or FI loci with detailed functional annotation, assessed their relevance in AA individuals, and sought previously undescribed loci through trans-ethnic (AA+EA) meta-analysis. We narrowed credible sets of variants driving association signals for 22/54 EA-associated loci; 18/22 credible sets overlapped with active islet-specific enhancers or transcription factor (TF) binding sites, and 21/22 contained at least one TF motif. Of the 54 EA-associated loci, 23 were shared between EA and AA. Replication with an additional 10,096 AA individuals identified two previously undescribed FI loci, chrX FAM133A (rs213676) and chr5 PELO (rs6450057). Trans-ethnic analyses with regulatory annotation illuminate the genetic architecture of glycemic traits and suggest gene regulation as a target to advance precision medicine for T2D. Our approach to utilize state-of-the-art functional annotation and implement trans-ethnic association analysis for discovery and fine-mapping offers a framework for further follow-up and characterization of GWAS signals of complex trait loci.

GCKR and PPP1R3B identified as genome-wide significant loci for plasma lactate: the Atherosclerosis Risk in Communities (ARIC) study

AIM

To investigate the genetic influence of circulating lactate level, a marker of oxidative capacity associated with diabetes.

METHODS

We conducted a genome-wide association study of log-transformed plasma lactate levels in 6901 European-American participants in the Atherosclerosis Risk in Communities study. For regions that achieved genome-wide significance in European-American participants, we conducted candidate region analysis in African-American subjects and tested for interaction between metformin use and the index single nucleotide polymorphisms for plasma lactate in European-American subjects.

RESULTS

The genome-wide association study in European-American subjects identified two genome-wide significant loci, GCKR (rs1260326, T allele β=0.08; P=1.8×10(-47) ) and PPP1R3B/LOC157273 (rs9987289, A allele β=0.06; P=1.6×10(-9) ). The index single nucleotide polymorphisms in these two loci explain 3.3% of the variance in log-transformed plasma lactate levels among the European-American subjects. In the African-American subjects, based on a region-significant threshold, the index single nucleotide polymorphism at GCKR was associated with plasma lactate but that at PPP1R3B/LOC157273 was not. Metformin use appeared to strengthen the association between the index single nucleotide polymorphism at PPP1R3B/LOC157273 and plasma lactate in European-American subjects (P for interaction=0.01).

CONCLUSIONS

We identified GCKR and PPP1R3B/LOC157273 as two genome-wide significant loci of plasma lactate. Both loci are associated with other diabetes-related phenotypes. These findings increase our understanding of the genetic control of lactate metabolism.

Electrophysiologic Substrate and Risk of Mortality in Incident Hemodialysis

The single leading cause of mortality on hemodialysis is sudden cardiac death. Whether measures of electrophysiologic substrate independently associate with mortality is unknown. We examined measures of electrophysiologic substrate in a prospective cohort of 571 patients on incident hemodialysis enrolled in the Predictors of Arrhythmic and Cardiovascular Risk in End Stage Renal Disease Study. A total of 358 participants completed both baseline 5-minute and 12-lead electrocardiogram recordings on a nondialysis day. Measures of electrophysiologic substrate included ventricular late potentials by the signal-averaged electrocardiogram and spatial mean QRS-T angle measured on the averaged beat recorded within a median of 106 days (interquartile range, 78-151 days) from dialysis initiation. The cohort was 59% men, and 73% were black, with a mean±SD age of 55±13 years. Transthoracic echocardiography revealed a mean±SD ejection fraction of 65.5%±12.0% and a mean±SD left ventricular mass index of 66.6±22.3 g/m^2.7 During 864.6 person-years of follow-up, 77 patients died; 35 died from cardiovascular causes, of which 15 were sudden cardiac deaths. By Cox regression analysis, QRS-T angle ≥75° significantly associated with increased risk of cardiovascular mortality (hazard ratio, 2.99; 95% confidence interval, 1.31 to 6.82) and sudden cardiac death (hazard ratio, 4.52; 95% confidence interval, 1.17 to 17.40) after multivariable adjustment for demographic, cardiovascular, and dialysis factors. Abnormal signal-averaged electrocardiogram measures did not associate with mortality. In conclusion, spatial QRS-T angle but not abnormal signal-averaged electrocardiogram significantly associates with cardiovascular mortality and sudden cardiac death independent of traditional risk factors in patients starting hemodialysis.

Examination of Potential Modifiers of the Association of APOL1 Alleles with CKD Progression

BACKGROUND AND OBJECTIVES

Common apolipoprotein L1 (APOL1) variants are associated with increased risk of progressive CKD; however, not all individuals with high-risk APOL1 variants experience CKD progression. Identification of factors contributing to heterogeneity has important scientific and clinical implications.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS

Using multivariable Cox models, we analyzed data from 693 participants in the African American Study of Kidney Disease and Hypertension to identify factors that modify the association between APOL1 genotypes and CKD progression (doubling of serum creatinine or incident ESRD).

RESULTS

Participant mean age was 54 years old, median GFR was 49 ml/min per 1.73 m(2), and 23% had the APOL1 high-risk genotype (two copies of the high-risk allele). Over a mean follow-up of 7.8 years, 288 (42%) participants experienced CKD progression. As previously reported, the high-risk genotype was associated with higher risk of CKD progression compared with the low-risk genotype (hazard ratio [HR], 1.88; 95% confidence interval [95% CI], 1.46 to 2.41). Although we found some suggestion that obesity (HR, 1.48; 95% CI, 1.05 to 2.08 and HR, 2.44; 95% CI, 1.66 to 3.57 for body mass index ≥ 30 versus <30 kg/m(2); P interaction =0.04) and increased urinary excretion of urea nitrogen (HR, 1.43; 95% CI, 0.98 to 2.09 versus HR, 2.33; 95% CI, 1.65 to 3.30 for urine urea nitrogen ≥ 8 versus <8 g/d; P interaction =0.04) were associated with lower APOL1-associated risk for CKD progression, these findings were not robust in sensitivity analyses with alternative cut points. No other sociodemographic (e.g., education and income), clinical (e.g., systolic BP and smoking), or laboratory (e.g., net endogenous acid production, urinary sodium and potassium excretions, 25-hydroxy vitamin D, intact parathyroid hormone, or fibroblast growth factor 23) variables modified the association between APOL1 and CKD progression (P interaction >0.05 for each).

CONCLUSIONS

Sociodemographic factors and common risk factors for CKD progression do not seem to alter APOL1-related CKD progression. Additional investigation is needed to identify nontraditional factors that may affect the association between APOL1 and progressive CKD.

Genome-Wide Association and Trans-ethnic Meta-Analysis for Advanced Diabetic Kidney Disease: Family Investigation of Nephropathy and Diabetes (FIND)

Diabetic kidney disease (DKD) is the most common etiology of chronic kidney disease (CKD) in the industrialized world and accounts for much of the excess mortality in patients with diabetes mellitus. Approximately 45% of U.S. patients with incident end-stage kidney disease (ESKD) have DKD. Independent of glycemic control, DKD aggregates in families and has higher incidence rates in African, Mexican, and American Indian ancestral groups relative to European populations. The Family Investigation of Nephropathy and Diabetes (FIND) performed a genome-wide association study (GWAS) contrasting 6,197 unrelated individuals with advanced DKD with healthy and diabetic individuals lacking nephropathy of European American, African American, Mexican American, or American Indian ancestry. A large-scale replication and trans-ethnic meta-analysis included 7,539 additional European American, African American and American Indian DKD cases and non-nephropathy controls. Within ethnic group meta-analysis of discovery GWAS and replication set results identified genome-wide significant evidence for association between DKD and rs12523822 on chromosome 6q25.2 in American Indians (P = 5.74x10^-9). The strongest signal of association in the trans-ethnic meta-analysis was with a SNP in strong linkage disequilibrium with rs12523822 (rs955333; P = 1.31x10^-8), with directionally consistent results across ethnic groups. These 6q25.2 SNPs are located between the SCAF8 and CNKSR3 genes, a region with DKD relevant changes in gene expression and an eQTL with IPCEF1, a gene co-translated with CNKSR3. Several other SNPs demonstrated suggestive evidence of association with DKD, within and across populations. These data identify a novel DKD susceptibility locus with consistent directions of effect across diverse ancestral groups and provide insight into the genetic architecture of DKD.

Type 2 diabetes is the most common cause of severe kidney disease worldwide and diabetic kidney disease (DKD) associates with premature death. Individuals of non-European ancestry have the highest burden of type 2 DKD; hence understanding the causes of DKD remains critical to reducing health disparities. Family studies demonstrate that genes regulate the onset and progression of DKD; however, identifying these genes has proven to be challenging. The Family Investigation of Diabetes and Nephropathy consortium (FIND) recruited a large multi-ethnic collection of individuals with type 2 diabetes with and without kidney disease in order to detect genes associated with DKD. FIND discovered and replicated a DKD-associated genetic locus on human chromosome 6q25.2 (rs955333) between the SCAF8 and CNKSR genes. Findings were supported by significantly different expression of genes in this region from kidney tissue of subjects with, versus without DKD. The present findings identify a novel kidney disease susceptibility locus in individuals with type 2 diabetes which is consistent across subjects of differing ancestries. In addition, FIND results provide a rich catalogue of genetic variation in DKD patients for future research on the genetic architecture regulating this common and devastating disease.

APOL1 Genotype and Glomerular and Tubular Kidney Injury in Women With HIV

BACKGROUND

APOL1 genotype is associated with advanced kidney disease in African Americans, but the pathogenic mechanisms are unclear. Here, associations of APOL1 genotype with urine biomarkers of glomerular and tubular injury and kidney function decline were evaluated.

STUDY DESIGN

Observational study.

SETTING & PARTICIPANTS

431 human immunodeficiency virus (HIV)-infected African American women enrolled in Women's Interagency HIV Study (WIHS).

PREDICTOR

APOL1 genotype.

OUTCOMES

Albumin-creatinine ratio (ACR), 4 tubular injury biomarkers (interleukin 18 [IL-18], kidney injury molecule 1 [KIM-1], neutrophil gelatinase-associated lipocalin [NGAL], and α1-microglobulin [A1M]), and kidney function estimated using the CKD-EPI cystatin C equation.

MEASUREMENTS

Participants were genotyped for APOL1 single-nucleotide polymorphisms rs73885319 (G1 allele) and rs71785313 (G2 allele). Urine biomarkers were measured using stored samples from 1999-2000. Cystatin C was measured using serum collected at baseline and 4- and 8-year follow-ups.

RESULTS

At baseline, ACRs were higher among 47 women with 2 APOL1 risk alleles versus 384 women with 0/1 risk allele (median, 24 vs 11mg/g; P<0.001). Compared with women with 0/1 risk allele, women with 2 risk alleles had 104% higher ACRs (95% CI, 29-223mg/g) and 2-fold greater risk of ACR>30 (95% CI, 1.17-3.44) mg/g after multivariable adjustment. APOL1 genotype showed little association with urine IL-18:Cr ratio, KIM-1:Cr ratio, and NGAL:Cr ratio (estimates of -5% [95% CI, -24% to 18%], -20% [95% CI, -36% to -1%], and 10% [95% CI, -26% to 64%], respectively) or detectable urine A1M (prevalence ratio, 1.13; 95% CI, 0.65-1.97) in adjusted analyses. Compared with women with 0/1 allele, women with 2 risk alleles had faster eGFR decline, by 1.2 (95% CI, 0.2 to 2.2) mL/min/1.73m(2) per year, and 1.7- and 3.4-fold greater rates of incident chronic kidney disease (95% CI, 1.1 to 2.5) and 10% annual eGFR decline (95% CI, 1.7 to 6.7), respectively, with minimal attenuation after adjustment for glomerular and tubular injury biomarker levels.

LIMITATIONS

Results may not be generalizable to men.

CONCLUSIONS

Among HIV-infected African American women, APOL1-associated kidney injury appears to localize to the glomerulus, rather than the tubules.

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.

Rationale and design for the Predictors of Arrhythmic and Cardiovascular Risk in End Stage Renal Disease (PACE) study

Background

Sudden cardiac death occurs commonly in the end-stage renal disease population receiving dialysis, with 25% dying of sudden cardiac death over 5 years. Despite this high risk, surprisingly few prospective studies have studied clinical- and dialysis-related risk factors for sudden cardiac death and arrhythmic precursors of sudden cardiac death in end-stage renal disease.

Methods/Design

We present a brief summary of the risk factors for arrhythmias and sudden cardiac death in persons with end-stage renal disease as the rationale for the Predictors of Arrhythmic and Cardiovascular Risk in End Stage Renal Disease (PACE) study, a prospective cohort study of patients recently initiated on chronic hemodialysis, with the overall goal to understand arrhythmic and sudden cardiac death risk. Participants were screened for eligibility and excluded if they already had a pacemaker or an automatic implantable cardioverter defibrillator. We describe the study aims, design, and data collection of 574 incident hemodialysis participants from the Baltimore region in Maryland, U.S.A.. Participants were recruited from 27 hemodialysis units and underwent detailed clinical, dialysis and cardiovascular evaluation at baseline and follow-up. Cardiovascular phenotyping was conducted on nondialysis days with signal averaged electrocardiogram, echocardiogram, pulse wave velocity, ankle, brachial index, and cardiac computed tomography and angiography conducted at baseline. Participants were followed annually with study visits including electrocardiogram, pulse wave velocity, and ankle brachial index up to 4 years. A biorepository of serum, plasma, DNA, RNA, and nails were collected to study genetic and serologic factors associated with disease.

Discussion

Studies of modifiable risk factors for sudden cardiac death will help set the stage for clinical trials to test therapies to prevent sudden cardiac death in this high-risk population.

Relative performance of gene- and pathway-level methods as secondary analyses for genome-wide association studies

BACKGROUND

Despite the success of genome-wide association studies (GWAS), there still remains "missing heritability" for many traits. One contributing factor may be the result of examining one marker at a time as opposed to a group of markers that are biologically meaningful in aggregate. To address this problem, a variety of gene- and pathway-level methods have been developed to identify putative biologically relevant associations. A simulation was conducted to systematically assess the performance of these methods. Using genetic data from 4,500 individuals in the Wellcome Trust Case Control Consortium (WTCCC), case-control status was simulated based on an additive polygenic model. We evaluated gene-level methods based on their sensitivity, specificity, and proportion of false positives. Pathway-level methods were evaluated on the relationship between proportion of causal genes within the pathway and the strength of association.

RESULTS

The gene-level methods had low sensitivity (20-63%), high specificity (89-100%), and low proportion of false positives (0.1-6%). The gene-level program VEGAS using only the top 10% of associated single nucleotide polymorphisms (SNPs) within the gene had the highest sensitivity (28.6%) with less than 1% false positives. The performance of the pathway-level methods depended on their reliance upon asymptotic distributions or if significance was estimated in a competitive manner. The pathway-level programs GenGen, GSA-SNP and MAGENTA had the best performance while accounting for potential confounders.

CONCLUSIONS

Novel genes and pathways can be identified using the gene and pathway-level methods. These methods may provide valuable insight into the "missing heritability" of traits and provide biological interpretations to GWAS findings.

Association of levels of fasting glucose and insulin with rare variants at the chromosome 11p11.2-MADD locus: Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium Targeted Sequencing Study

BACKGROUND

Common variation at the 11p11.2 locus, encompassing MADD, ACP2, NR1H3, MYBPC3, and SPI1, has been associated in genome-wide association studies with fasting glucose and insulin (FI). In the Cohorts for Heart and Aging Research in Genomic Epidemiology Targeted Sequencing Study, we sequenced 5 gene regions at 11p11.2 to identify rare, potentially functional variants influencing fasting glucose or FI levels.

METHODS AND RESULTS

Sequencing (mean depth, 38×) across 16.1 kb in 3566 individuals without diabetes mellitus identified 653 variants, 79.9% of which were rare (minor allele frequency <1%) and novel. We analyzed rare variants in 5 gene regions with FI or fasting glucose using the sequence kernel association test. At NR1H3, 53 rare variants were jointly associated with FI (P=2.73×10(-3)); of these, 7 were predicted to have regulatory function and showed association with FI (P=1.28×10(-3)). Conditioning on 2 previously associated variants at MADD (rs7944584, rs10838687) did not attenuate this association, suggesting that there are >2 independent signals at 11p11.2. One predicted regulatory variant, chr11:47227430 (hg18; minor allele frequency=0.00068), contributed 20.6% to the overall sequence kernel association test score at NR1H3, lies in intron 2 of NR1H3, and is a predicted binding site for forkhead box A1 (FOXA1), a transcription factor associated with insulin regulation. In human HepG2 hepatoma cells, the rare chr11:47227430 A allele disrupted FOXA1 binding and reduced FOXA1-dependent transcriptional activity.

CONCLUSIONS

Sequencing at 11p11.2-NR1H3 identified rare variation associated with FI. One variant, chr11:47227430, seems to be functional, with the rare A allele reducing transcription factor FOXA1 binding and FOXA1-dependent transcriptional activity.

Modulation of genetic associations with serum urate levels by body-mass-index in humans

We tested for interactions between body mass index (BMI) and common genetic variants affecting serum urate levels, genome-wide, in up to 42569 participants. Both stratified genome-wide association (GWAS) analyses, in lean, overweight and obese individuals, and regression-type analyses in a non BMI-stratified overall sample were performed. The former did not uncover any novel locus with a major main effect, but supported modulation of effects for some known and potentially new urate loci. The latter highlighted a SNP at RBFOX3 reaching genome-wide significant level (effect size 0.014, 95% CI 0.008-0.02, Pinter= 2.6 x 10-8). Two top loci in interaction term analyses, RBFOX3 and ERO1LB-EDARADD, also displayed suggestive differences in main effect size between the lean and obese strata. All top ranking loci for urate effect differences between BMI categories were novel and most had small magnitude but opposite direction effects between strata. They include the locus RBMS1-TANK (men, Pdifflean-overweight= 4.7 x 10-8), a region that has been associated with several obesity related traits, and TSPYL5 (men, Pdifflean-overweight= 9.1 x 10-8), regulating adipocytes-produced estradiol. The top-ranking known urate loci was ABCG2, the strongest known gout risk locus, with an effect halved in obese compared to lean men (Pdifflean-obese= 2 x 10-4). Finally, pathway analysis suggested a role for N-glycan biosynthesis as a prominent urate-associated pathway in the lean stratum. These results illustrate a potentially powerful way to monitor changes occurring in obesogenic environment.

The association between APOL1 risk alleles and longitudinal kidney function differs by HIV viral suppression status

BACKGROUND

Existing data suggest that human immunodeficiency virus (HIV)-infected African Americans carrying 2 copies of the APOL1 risk alleles have greater risk of kidney disease than noncarriers. We sought to determine whether HIV RNA suppression mitigates APOL1-related kidney function decline among African Americans enrolled in the Multicenter AIDS Cohort Study.

METHODS

We genotyped HIV-infected men for the G1 and G2 risk alleles and ancestry informative markers. Mixed-effects models were used to estimate the annual rate of estimated glomerular filtration rate (eGFR) decline, comparing men carrying 2 (high-risk) vs 0-1 risk allele (low-risk). Effect modification by HIV suppression status (defined as HIV type 1 RNA level <400 copies/mL for >90% of follow-up time) was evaluated using interaction terms and stratified analyses.

RESULTS

Of the 333 African American men included in this study, 54 (16%) carried the APOL1 high-risk genotype. Among HIV-infected men with unsuppressed viral loads, those with the high-risk genotype had a 2.42 mL/minute/1.73 m(2) (95% confidence interval [CI], -3.52 to -1.32) faster annual eGFR decline than men with the low-risk genotype. This association was independent of age, comorbid conditions, baseline eGFR, ancestry, and HIV-related factors. In contrast, the rate of decline was similar by APOL1 genotype among men with sustained viral suppression (-0.16 mL/minute/1.73 m(2)/year; 95% CI, -.59 to .27; P for interaction <.001).

CONCLUSIONS

Unsuppressed HIV-infected African Americans with the APOL1 high-risk genotype experience an accelerated rate of kidney function decline; HIV suppression with antiretroviral therapy may reduce these deleterious renal effects.

Both rare and common variants in PCSK9 influence plasma low-density lipoprotein cholesterol level in American Indians

CONTEXT

Elevated LDL cholesterol (LDL-C) is an important risk factor for atherosclerosis and cardiovascular disease. Variants in the proprotein convertase subtilisin/kexin type 9 (PCSK9) gene have been associated not only with plasma LDL-C concentration, but also with ischemic heart disease. Little is known about the genetic architecture of PCSK9 and its influence on LDL-C in American Indians.

OBJECTIVE

We aimed to investigate the genetic architecture in the 1p32 region encompassing PCSK9 and its influence on LDL-C in American Indians.

DESIGN

The Strong Heart Family Study (SHFS) is a family-based genetic study.

PARTICIPANTS

Two thousand four hundred fifty eight American Indians from Arizona, Oklahoma, North Dakota, and South Dakota, who were genotyped by Illumina MetaboChip.

RESULTS

We genotyped 486 SNPs in a 3.9 Mb region at chromosome 1p32 encompassing PCSK9 in 2458 American Indians. We examined the association between these SNPs and LDL-C. For common variants (MAF ≥ 1%), meta-analysis across the three geographic regions showed common variants in PCSK9 were significantly associated with higher LDL-C. The most significant SNP rs12067569 (MAF = 1.7 %, β = 16.9 ± 3.7, P = 5.9 × 10(-6)) was in complete LD (r(2) = 1) with a nearby missense SNP, rs505151 (E670G) (β = 15.0 ± 3.6, P = 3.6 × 10(-5)). For rare variants (MAF < 1%), rs11591147 (R46L, MAF = 0.9%) was associated with lower LDL-C (β = - 31.1 ± 7.1, P = 1.4 × 10(-5)). The mean (SD) of LDL-C was 76.9 (7.8) and 107.4 (1.0) mg/dL for those with and without the R46L mutation, respectively. One person who was homozygous for R46L had LDL-C levels of 11 mg/dL. In one family, 6 out of 8 members carrying the R46L mutation had LDL-C levels below the lower 10% percentile of LDL-C among all study participants.

CONCLUSIONS

Both rare and common variants in PCSK9 influence plasma LDL-C levels in American Indians. Follow-up studies may disclose the influence of these mutations on the risk of CVD and responses to cholesterol-lowering medications.

Effect of zinc supplementation on insulin secretion: interaction between zinc and SLC30A8 genotype in Old Order Amish

AIMS/HYPOTHESIS

SLC30A8 encodes a zinc transporter in the beta cell; individuals with a common missense variant (rs13266634; R325W) in SLC30A8 demonstrate a lower early insulin response to glucose and an increased risk of type 2 diabetes. We hypothesised that zinc supplementation may improve insulin secretion in a genotype-dependent manner.

METHODS

We evaluated the early insulin response to glucose (using frequently sampled intravenous glucose tolerance testing) by R325W genotype before and after 14 days of supplementation with oral zinc acetate (50 mg elemental zinc) twice daily in healthy non-diabetic Amish individuals (N = 55).

RESULTS

Individuals with RW/WW genotypes (n = 32) had the lowest insulin response to glucose at 5 and 10 min at baseline (vs RR homozygotes [n = 23]). After zinc supplementation, the RW/WW group experienced 15% and 14% increases in the insulin response to glucose at 5 and 10 min, respectively (p ≤ 0.04), and, compared with RR homozygotes, experienced a 26% (p = 0.04) increase in insulin at 5 min. We observed reciprocal decreases in proinsulin:insulin in the RW/WW (p = 0.002) vs RR group (p = 0.048), suggesting a genotype-specific improvement in insulin processing.

CONCLUSIONS/INTERPRETATION

Zinc supplementation appears to affect the early insulin response to glucose differentially by rs13266634 genotype and could be beneficial for diabetes prevention and/or treatment for some individuals based on SLC30A8 variation.

TRIAL REGISTRATION

ClinicalTrials.gov NCT00981448.

Low-frequency and rare exome chip variants associate with fasting glucose and type 2 diabetes susceptibility

Fasting glucose and insulin are intermediate traits for type 2 diabetes. Here we explore the role of coding variation on these traits by analysis of variants on the HumanExome BeadChip in 60,564 non-diabetic individuals and in 16,491 T2D cases and 81,877 controls. We identify a novel association of a low-frequency nonsynonymous SNV in GLP1R (A316T; rs10305492; MAF=1.4%) with lower FG (β=-0.09±0.01 mmol l(-1), P=3.4 × 10(-12)), T2D risk (OR[95%CI]=0.86[0.76-0.96], P=0.010), early insulin secretion (β=-0.07±0.035 pmolinsulin mmolglucose(-1), P=0.048), but higher 2-h glucose (β=0.16±0.05 mmol l(-1), P=4.3 × 10(-4)). We identify a gene-based association with FG at G6PC2 (pSKAT=6.8 × 10(-6)) driven by four rare protein-coding SNVs (H177Y, Y207S, R283X and S324P). We identify rs651007 (MAF=20%) in the first intron of ABO at the putative promoter of an antisense lncRNA, associating with higher FG (β=0.02±0.004 mmol l(-1), P=1.3 × 10(-8)). Our approach identifies novel coding variant associations and extends the allelic spectrum of variation underlying diabetes-related quantitative traits and T2D susceptibility.

Genome-wide association analysis of plasma B-type natriuretic peptide in blacks: the Jackson Heart Study

BACKGROUND

Numerous experimental studies suggest that B-type natriuretic peptide (BNP) is cardioprotective; however, in clinical studies, higher plasma BNP concentrations have been associated with incident cardiovascular disease and higher left ventricular mass. Genetic association studies may allow us to determine the true causal directions without confounding by compensatory mechanisms.

METHODS AND RESULTS

We performed a meta-analysis of 2 genome-wide association results from a total of 2790 blacks. We assumed an additive genetic model in an association analysis of imputed 2.5 million single-nucleotide polymorphism dosages with residuals generated from multivariable-adjusted logarithmically transformed BNP controlling for relevant covariates and population stratification. Two loci were genome-wide significant, a candidate gene locus NPPB (rs198389, P=1.18×10(-09)) and a novel missense variant in the KLKB1 locus (rs3733402, P=1.75×10(-11)) that explained 0.4% and 1.9% of variation in log BNP concentration, respectively. The observed increase in BNP concentration was proportional to the number of effect allele copies, and an average of 8.1 pg/mL increase was associated with 2 allele copies. In a companion study, single-nucleotide polymorphisms in this loci were cross-checked with genome-wide association results for the aldosterone/renin ratio in individuals of European ancestry, and rs3733402 was genome-wide significant (P<5.0×10(-8)), suggesting possible shared genetic architecture for these 2 pathways. Other statistically significant relations for these single-nucleotide polymorphisms included the following: rs198389 with systolic blood pressure in blacks (COGENT consortium) and rs198389 and rs3733402 with left ventricular mass in whites (EchoGEN consortium).

CONCLUSIONS

These findings improve our knowledge of the genetic basis of BNP variation in blacks, demonstrate a possible shared allelic architecture for BNP with aldosterone-renin ratio, and motivate further studies of underlying mechanisms.

Polygenic overlap between kidney function and large artery atherosclerotic stroke

BACKGROUND AND PURPOSE

Epidemiological studies show strong associations between kidney dysfunction and risk of ischemic stroke (IS), the mechanisms of which are incompletely understood. We investigated whether these associations may reflect shared heritability because of a common polygenic basis and whether this differed for IS subtypes.

METHODS

Polygenic models were derived using genome-wide association studies meta-analysis results for 3 kidney traits: estimated glomerular filtration rate using serum creatinine (eGFRcrea: n=73 998), eGFR using cystatin C (eGFRcys: n=22 937), and urinary albumin to creatinine ratio (n=31 580). For each, single nucleotide polymorphisms passing 10 P value thresholds were used to form profile scores in 4561 IS cases and 7094 controls from the United Kingdom, Germany, and Australia. Scores were tested for association with IS and its 3 aetiological subtypes: large artery atherosclerosis, cardioembolism, and small vessel disease.

RESULTS

Polygenic scores correlating with higher eGFRcrea were associated with reduced risk of large artery atherosclerosis, with 5 scores reaching P<0.05 (peak P=0.004) and all showing the epidemiologically expected direction of effect. A similar pattern was observed for polygenic scores reflecting higher urinary albumin to creatinine ratio, of which 3 associated with large artery atherosclerosis (peak P=0.01) and all showed the expected directional association. One urinary albumin to creatinine ratio-based score also associated with small vessel disease (P=0.03). The global pattern of results was unlikely to have occurred by chance (P=0.02).

CONCLUSIONS

This study suggests possible polygenic correlation between renal dysfunction and IS. The shared genetic components may be specific to stroke subtypes, particularly large artery atherosclerotic stroke. Further study of the genetic relationships between these disorders seems merited.

Estimating time to ESRD using kidney failure risk equations: results from the African American Study of Kidney Disease and Hypertension (AASK)

BACKGROUND

Planning for renal replacement therapy, such as referral for arteriovenous fistula placement and transplantation, often is guided by level of estimated glomerular filtration rate (eGFR). The use of risk equations might enable more accurate estimation of time to end-stage renal disease (ESRD), thus improving patient care.

STUDY DESIGN

Prospective observational study.

SETTING & PARTICIPANTS

1,094 participants in the African American Study of Kidney Disease and Hypertension (AASK) cohort.

PREDICTOR

Age, sex, urine protein-creatinine ratio ≥ 1g/g, APOL1 high-risk status, and 3-year antecedent eGFR decline.

OUTCOME

Cumulative incidence of ESRD from 5 different starting points: eGFR of 30 and 15mL/min/1.73m(2) and 1-year ESRD risk of 5%, 10%, and 20%, estimated by a published 4-variable kidney failure risk equation.

RESULTS

566 participants developed eGFR of 30mL/min/1.73m(2), 244 developed eGFR of 15mL/min/1.73m(2), and 437, 336, and 259 developed 1-year ESRD risks of 5%, 10%, and 20%, respectively. The 1-year cumulative incidence of ESRD was 4.3% from eGFR of 30mL/min/1.73m(2), 49.0% from eGFR of 15mL/min/1.73m(2), 6.7% from 5% ESRD risk, 15.0% from 10% ESRD risk, and 29% from 20% ESRD risk. From eGFR of 30mL/min/1.73m(2), there were several risk factors that predicted ESRD risk. From eGFR of 15mL/min/1.73m(2), only level of proteinuria did; median time to ESRD was 9 and 19 months in those with higher and lower proteinuria, respectively. Median times were less variable from corresponding ESRD risk thresholds. For example, median times to ESRD from 20% ESRD risk were 22 and 25 months among those with higher and lower proteinuria, respectively.

LIMITATIONS

Relatively homogeneous population of African Americans with hypertensive kidney disease.

CONCLUSIONS

Results of the present study suggest the potential benefit of incorporating kidney failure risk equations into clinical care, with selection of a specific threshold guided by its intended use.

Meta-analysis of genome-wide association studies in African Americans provides insights into the genetic architecture of type 2 diabetes

Type 2 diabetes (T2D) is more prevalent in African Americans than in Europeans. However, little is known about the genetic risk in African Americans despite the recent identification of more than 70 T2D loci primarily by genome-wide association studies (GWAS) in individuals of European ancestry. In order to investigate the genetic architecture of T2D in African Americans, the MEta-analysis of type 2 DIabetes in African Americans (MEDIA) Consortium examined 17 GWAS on T2D comprising 8,284 cases and 15,543 controls in African Americans in stage 1 analysis. Single nucleotide polymorphisms (SNPs) association analysis was conducted in each study under the additive model after adjustment for age, sex, study site, and principal components. Meta-analysis of approximately 2.6 million genotyped and imputed SNPs in all studies was conducted using an inverse variance-weighted fixed effect model. Replications were performed to follow up 21 loci in up to 6,061 cases and 5,483 controls in African Americans, and 8,130 cases and 38,987 controls of European ancestry. We identified three known loci (TCF7L2, HMGA2 and KCNQ1) and two novel loci (HLA-B and INS-IGF2) at genome-wide significance (4.15×10⁻⁹⁴<P<5×10⁻⁸, odds ratio (OR) = 1.09 to 1.36). Fine-mapping revealed that 88 of 158 previously identified T2D or glucose homeostasis loci demonstrated nominal to highly significant association (2.2×10⁻²³ < locus-wide P<0.05). These novel and previously identified loci yielded a sibling relative risk of 1.19, explaining 17.5% of the phenotypic variance of T2D on the liability scale in African Americans. Overall, this study identified two novel susceptibility loci for T2D in African Americans. A substantial number of previously reported loci are transferable to African Americans after accounting for linkage disequilibrium, enabling fine mapping of causal variants in trans-ethnic meta-analysis studies.

Despite the higher prevalence of type 2 diabetes (T2D) in African Americans than in Europeans, recent genome-wide association studies (GWAS) were examined primarily in individuals of European ancestry. In this study, we performed meta-analysis of 17 GWAS in 8,284 cases and 15,543 controls to explore the genetic architecture of T2D in African Americans. Following replication in additional 6,061 cases and 5,483 controls in African Americans, and 8,130 cases and 38,987 controls of European ancestry, we identified two novel and three previous reported T2D loci reaching genome-wide significance. We also examined 158 loci previously reported to be associated with T2D or regulating glucose homeostasis. While 56% of these loci were shared between African Americans and the other populations, the strongest associations in African Americans are often found in nearby single nucleotide polymorphisms (SNPs) instead of the original SNPs reported in other populations due to differential genetic architecture across populations. Our results highlight the importance of performing genetic studies in non-European populations to fine map the causal genetic variants.

Explaining the racial difference in AKI incidence

African Americans face higher risk of AKI than Caucasians. The extent to which this increased risk is because of differences in clinical, socioeconomic, or genetic risk factors is unknown. We evaluated 10,588 African-American and Caucasian participants in the Atherosclerosis Risk in Communities study, a community-based prospective cohort of middle-aged individuals. Participants were followed from baseline study visit (1996-1999) to first hospitalization for AKI (defined by billing code), ESRD, death, or December 31, 2010. African-American participants were slightly younger (61.7 versus 63.1 years, P<0.001), were more often women (64.5% versus 53.2%, P<0.001), and had higher baseline eGFR compared with Caucasians. Annual family income, education level, and prevalence of health insurance were lower among African Americans than Caucasians. The unadjusted incidence of hospitalized AKI was 7.4 cases per 1000 person-years among African Americans and 5.8 cases per 1000 person-years among Caucasians (P=0.002). The elevated risk of AKI among African Americans persisted after adjustment for demographics, cardiovascular risk factors, kidney markers, and time-varying number of hospitalizations (adjusted hazard ratio, 1.20; 95% confidence interval [95% CI], 1.01 to 1.43; P=0.04); however, accounting for differences in income and/or insurance by race attenuated the association (P>0.05). High-risk APOL1 variants did not associate with AKI among African Americans (demographic-adjusted hazard ratio, 1.07; 95% CI, 0.69 to 1.65; P=0.77). In summary, the higher risk of AKI among African Americans may be related to disparities in socioeconomic status.

Pleiotropic genes for metabolic syndrome and inflammation

Metabolic syndrome (MetS) has become a health and financial burden worldwide. The MetS definition captures clustering of risk factors that predict higher risk for diabetes mellitus and cardiovascular disease. Our study hypothesis is that additional to genes influencing individual MetS risk factors, genetic variants exist that influence MetS and inflammatory markers forming a predisposing MetS genetic network. To test this hypothesis a staged approach was undertaken. (a) We analyzed 17 metabolic and inflammatory traits in more than 85,500 participants from 14 large epidemiological studies within the Cross Consortia Pleiotropy Group. Individuals classified with MetS (NCEP definition), versus those without, showed on average significantly different levels for most inflammatory markers studied. (b) Paired average correlations between 8 metabolic traits and 9 inflammatory markers from the same studies as above, estimated with two methods, and factor analyses on large simulated data, helped in identifying 8 combinations of traits for follow-up in meta-analyses, out of 130,305 possible combinations between metabolic traits and inflammatory markers studied. (c) We performed correlated meta-analyses for 8 metabolic traits and 6 inflammatory markers by using existing GWAS published genetic summary results, with about 2.5 million SNPs from twelve predominantly largest GWAS consortia. These analyses yielded 130 unique SNPs/genes with pleiotropic associations (a SNP/gene associating at least one metabolic trait and one inflammatory marker). Of them twenty-five variants (seven loci newly reported) are proposed as MetS candidates. They map to genes MACF1, KIAA0754, GCKR, GRB14, COBLL1, LOC646736-IRS1, SLC39A8, NELFE, SKIV2L, STK19, TFAP2B, BAZ1B, BCL7B, TBL2, MLXIPL, LPL, TRIB1, ATXN2, HECTD4, PTPN11, ZNF664, PDXDC1, FTO, MC4R and TOMM40. Based on large data evidence, we conclude that inflammation is a feature of MetS and several gene variants show pleiotropic genetic associations across phenotypes and might explain a part of MetS correlated genetic architecture. These findings warrant further functional investigation.

Genetic association study of QT interval highlights role for calcium signaling pathways in myocardial repolarization

The QT interval, an electrocardiographic measure reflecting myocardial repolarization, is a heritable trait. QT prolongation is a risk factor for ventricular arrhythmias and sudden cardiac death (SCD) and could indicate the presence of the potentially lethal Mendelian Long QT Syndrome (LQTS). Using a genome-wide association and replication study in up to 100,000 individuals we identified 35 common variant QT interval loci, that collectively explain ∼8-10% of QT variation and highlight the importance of calcium regulation in myocardial repolarization. Rare variant analysis of 6 novel QT loci in 298 unrelated LQTS probands identified coding variants not found in controls but of uncertain causality and therefore requiring validation. Several newly identified loci encode for proteins that physically interact with other recognized repolarization proteins. Our integration of common variant association, expression and orthogonal protein-protein interaction screens provides new insights into cardiac electrophysiology and identifies novel candidate genes for ventricular arrhythmias, LQTS,and SCD.

Genome-wide association study identifies common loci influencing circulating glycated hemoglobin (HbA1c) levels in non-diabetic subjects: the Long Life Family Study (LLFS)

OBJECTIVE

Glycated hemoglobin (HbA1c) is a stable index of chronic glycemic status and hyperglycemia associated with progressive development of insulin resistance and frank diabetes. It is also associated with premature aging and increased mortality. To uncover novel loci for HbA1c that are associated with healthy aging, we conducted a genome-wide association study (GWAS) using non-diabetic participants in the Long Life Family Study (LLFS), a study with familial clustering of exceptional longevity in the US and Denmark.

METHODS

A total of 4088 non-diabetic subjects from the LLFS were used for GWAS discoveries, and a total of 8231 non-diabetic subjects from the Atherosclerosis Risk in Communities Study (ARIC, in the MAGIC Consortium) and the Health, Aging, and Body Composition Study (HABC) were used for GWAS replications. HbA1c was adjusted for age, sex, centers, 20 principal components, without and with BMI. A linear mixed effects model was used for association testing.

RESULTS

Two known loci at GCK rs730497 (or rs2908282) and HK1 rs17476364 were confirmed (p<5e-8). Of 25 suggestive (5e-8<p<1e-5) loci, one known (G6PC2 rs560887, replication p=5e-5) and one novel (OR10R3P/SPTA1- rs12041363, replication p=1e-17) loci were replicated (p<0.0019). Similar findings resulted when HbA1c was further adjusted for BMI. Further validations are crucial for the remaining suggestive loci including the emerged variant near OR10R3P/SPTA1.

CONCLUSIONS

The analysis reconfirmed two known GWAS loci (GCK, HK1) and identified 25 suggestive loci including one reconfirmed variant in G6PC2 and one replicated variant near OR10R3P/SPTA1. Future focused survey of sequence elements containing mainly functional and regulatory variants may yield additional findings.

Biomarkers of vascular calcification and mortality in patients with ESRD

BACKGROUND

Vascular calcification is common among patients undergoing dialysis and is associated with mortality. Factors such as osteoprotegerin (OPG), osteopontin (OPN), bone morphogenic protein-7 (BMP-7), and fetuin-A are involved in vascular calcification.

DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS

OPG, OPN, BMP-7, and fetuin-A were measured in blood samples from 602 incident dialysis patients recruited from United States dialysis centers between 1995 and 1998 as part of the Choices for Healthy Outcomes In Caring for ESRD Study. Their association with all-cause and cardiovascular mortality were assessed using Cox proportional hazards models adjusted for demographic characteristics, comorbidity, serum phosphate, and calcium. An interaction with diabetes was tested because of its known association with vascular calcification. Predictive accuracy of selected biomarkers was explored by C-statistics in nested models with training and validation subcohorts.

RESULTS

Higher OPG and lower fetuin-A levels were associated with higher mortality over up to 13 years of follow-up (median, 3.4 years). The adjusted hazard ratios (HR) for highest versus lowest tertile were 1.49 (95% confidence interval [95% CI], 1.08 to 2.06) for OPG and 0.69 (95% CI, 0.52 to 0.92) for fetuin-A. In stratified models, the highest tertile of OPG was associated with higher mortality among patients without diabetes (HR, 2.42; 95% CI, 1.35 to 4.34), but not patients with diabetes (HR, 1.26; 95% CI, 0.82 to 1.93; P for interaction=0.001). In terms of cardiovascular mortality, higher fetuin-A was associated with lower risk (HR, 0.85 per 0.1 g/L: 95% CI, 0.75 to 0.96). In patients without diabetes, higher OPG was associated with greater risk (HR for highest versus lowest tertile, 2.91; 95% CI, 1.06 to 7.99), but not in patients with diabetes or overall. OPN and BMP-7 were not independently associated with outcomes overall. The addition of OPG and fetuin-A did not significantly improve predictive accuracy of mortality.

CONCLUSIONS

OPG and fetuin-A may be risk factors for all-cause and cardiovascular mortality in patients undergoing dialysis, but do not improve risk prediction.

The pharmacogenetics of type 2 diabetes: a systematic review

OBJECTIVE

We performed a systematic review to identify which genetic variants predict response to diabetes medications.

RESEARCH DESIGN AND METHODS

We performed a search of electronic databases (PubMed, EMBASE, and Cochrane Database) and a manual search to identify original, longitudinal studies of the effect of diabetes medications on incident diabetes, HbA1c, fasting glucose, and postprandial glucose in prediabetes or type 2 diabetes by genetic variation. Two investigators reviewed titles, abstracts, and articles independently. Two investigators abstracted data sequentially and evaluated study quality independently. Quality evaluations were based on the Strengthening the Reporting of Genetic Association Studies guidelines and Human Genome Epidemiology Network guidance.

RESULTS

Of 7,279 citations, we included 34 articles (N = 10,407) evaluating metformin (n = 14), sulfonylureas (n = 4), repaglinide (n = 8), pioglitazone (n = 3), rosiglitazone (n = 4), and acarbose (n = 4). Studies were not standalone randomized controlled trials, and most evaluated patients with diabetes. Significant medication-gene interactions for glycemic outcomes included 1) metformin and the SLC22A1, SLC22A2, SLC47A1, PRKAB2, PRKAA2, PRKAA1, and STK11 loci; 2) sulfonylureas and the CYP2C9 and TCF7L2 loci; 3) repaglinide and the KCNJ11, SLC30A8, NEUROD1/BETA2, UCP2, and PAX4 loci; 4) pioglitazone and the PPARG2 and PTPRD loci; 5) rosiglitazone and the KCNQ1 and RBP4 loci; and 5) acarbose and the PPARA, HNF4A, LIPC, and PPARGC1A loci. Data were insufficient for meta-analysis.

CONCLUSIONS

We found evidence of pharmacogenetic interactions for metformin, sulfonylureas, repaglinide, thiazolidinediones, and acarbose consistent with their pharmacokinetics and pharmacodynamics. While high-quality controlled studies with prespecified analyses are still lacking, our results bring the promise of personalized medicine in diabetes one step closer to fruition.

Association of a cystatin C gene variant with cystatin C levels, CKD, and risk of incident cardiovascular disease and mortality

BACKGROUND

Carriers of the T allele of the single-nucleotide polymorphism rs13038305 tend to have lower cystatin C levels and higher cystatin C-based estimated glomerular filtration rate (eGFRcys). Adjusting for this genetic effect on cystatin C concentrations may improve GFR estimation, reclassify cases of chronic kidney disease (CKD), and strengthen risk estimates for cardiovascular disease (CVD) and mortality.

STUDY DESIGN

Observational.

SETTING & POPULATION

4 population-based cohorts: Atherosclerosis Risk in Communities (ARIC), Cardiovascular Health (CHS), Framingham Heart (FHS), and Health, Aging, and Body Composition (Health ABC) studies.

PREDICTORS

We estimated the association of rs13038305 with eGFRcys and serum creatinine-based eGFR (eGFRcr) and performed longitudinal analyses of the associations of eGFRcys with mortality and cardiovascular events following adjustment for rs13038305.

OUTCOMES

We assessed reclassification by genotype-adjusted eGFRcys across CKD categories: <45, 45-59, 60-89, and ≥ 90 mL/min/1.73 m(2). We compared mortality and CVD outcomes in those reclassified to a worse eGFRcys category with those unaffected. Results were combined using fixed-effect inverse-variance meta-analysis.

RESULTS

In 14,645 participants, each copy of the T allele of rs13038305 (frequency, 21%) was associated with a 6.4% lower cystatin C concentration, 5.5-mL/min/1.73 m(2) higher eGFRcys, and 36% [95% CI, 29%-41%] lower odds of CKD. Associations with CVD (HR, 1.17; 95% CI, 1.14-1.20) and mortality (HR, 1.22; 95% CI, 1.19-1.24) per 10-mL/min/1.73 m(2) lower eGFRcys were similar with or without rs13038305 adjustment. 1,134 (7.7%) participants were reclassified to a worse CKD category following rs13038305 adjustment, and rates of CVD and mortality were higher in individuals who were reclassified. However, the overall net reclassification index was not significant for either outcome, at 0.009 (95% CI, -0.003 to 0.022) for mortality and 0.014 (95% CI, 0.0 to 0.028) for CVD.

LIMITATIONS

rs13038305 explains only a small proportion of cystatin C variation.

CONCLUSIONS

Statistical adjustment can correct a genetic bias in GFR estimates based on cystatin C in carriers of the T allele of rs13038305 and result in changes in disease classification. However, on a population level, the effects on overall reclassification of CKD status are modest.

A genome-wide search for linkage of estimated glomerular filtration rate (eGFR) in the Family Investigation of Nephropathy and Diabetes (FIND)

Objective

Estimated glomerular filtration rate (eGFR), a measure of kidney function, is heritable, suggesting that genes influence renal function. Genes that influence eGFR have been identified through genome-wide association studies. However, family-based linkage approaches may identify loci that explain a larger proportion of the heritability. This study used genome-wide linkage and association scans to identify quantitative trait loci (QTL) that influence eGFR.

Methods

Genome-wide linkage and sparse association scans of eGFR were performed in families ascertained by probands with advanced diabetic nephropathy (DN) from the multi-ethnic Family Investigation of Nephropathy and Diabetes (FIND) study. This study included 954 African Americans (AA), 781 American Indians (AI), 614 European Americans (EA) and 1,611 Mexican Americans (MA). A total of 3,960 FIND participants were genotyped for 6,000 single nucleotide polymorphisms (SNPs) using the Illumina Linkage IVb panel. GFR was estimated by the Modification of Diet in Renal Disease (MDRD) formula.

Results

The non-parametric linkage analysis, accounting for the effects of diabetes duration and BMI, identified the strongest evidence for linkage of eGFR on chromosome 20q11 (log of the odds [LOD] = 3.34; P = 4.4×10⁻⁵) in MA and chromosome 15q12 (LOD = 2.84; P = 1.5×10⁻⁴) in EA. In all subjects, the strongest linkage signal for eGFR was detected on chromosome 10p12 (P = 5.5×10⁻⁴) at 44 cM near marker rs1339048. A subsequent association scan in both ancestry-specific groups and the entire population identified several SNPs significantly associated with eGFR across the genome.

Conclusion

The present study describes the localization of QTL influencing eGFR on 20q11 in MA, 15q21 in EA and 10p12 in the combined ethnic groups participating in the FIND study. Identification of causal genes/variants influencing eGFR, within these linkage and association loci, will open new avenues for functional analyses and development of novel diagnostic markers for DN.

APOL1 risk variants, race, and progression of chronic kidney disease

BACKGROUND

Among patients in the United States with chronic kidney disease, black patients are at increased risk for end-stage renal disease, as compared with white patients.

METHODS

In two studies, we examined the effects of variants in the gene encoding apolipoprotein L1 (APOL1) on the progression of chronic kidney disease. In the African American Study of Kidney Disease and Hypertension (AASK), we evaluated 693 black patients with chronic kidney disease attributed to hypertension. In the Chronic Renal Insufficiency Cohort (CRIC) study, we evaluated 2955 white patients and black patients with chronic kidney disease (46% of whom had diabetes) according to whether they had 2 copies of high-risk APOL1 variants (APOL1 high-risk group) or 0 or 1 copy (APOL1 low-risk group). In the AASK study, the primary outcome was a composite of end-stage renal disease or a doubling of the serum creatinine level. In the CRIC study, the primary outcomes were the slope in the estimated glomerular filtration rate (eGFR) and the composite of end-stage renal disease or a reduction of 50% in the eGFR from baseline.

RESULTS

In the AASK study, the primary outcome occurred in 58.1% of the patients in the APOL1 high-risk group and in 36.6% of those in the APOL1 low-risk group (hazard ratio in the high-risk group, 1.88; P<0.001). There was no interaction between APOL1 status and trial interventions or the presence of baseline proteinuria. In the CRIC study, black patients in the APOL1 high-risk group had a more rapid decline in the eGFR and a higher risk of the composite renal outcome than did white patients, among those with diabetes and those without diabetes (P<0.001 for all comparisons).

CONCLUSIONS

Renal risk variants in APOL1 were associated with the higher rates of end-stage renal disease and progression of chronic kidney disease that were observed in black patients as compared with white patients, regardless of diabetes status. (Funded by the National Institute of Diabetes and Digestive and Kidney Diseases and others.).

Association of SERPINA9 gene variants with carotid artery atherosclerosis: the Atherosclerosis Risk in Communities (ARIC) Carotid MRI Study

The SNP rs11628722 in the SERPINA9 gene was previously associated with incident ischemic stroke in the Atherosclerosis Risk in Communities (ARIC) study. Centerin, the protein encoded by SERPINA9, is involved in maturation and maintenance of naïve B cells, which play a role in atherogenesis. We investigated whether 21 tag SNPs in the SERPINA9 gene are associated with features of carotid artery atherosclerotic plaque measured by magnetic resonance imaging (MRI). Carotid MRI data were obtained from 1,282 European Americans and 341 African Americans of the ARIC Carotid MRI study, which recruited participants from ARIC by a stratified sampling plan that over-sampled participants with carotid intima-media thickening. Five MRI measures, focused on carotid wall volume, wall thickness, and lipid core, were analyzed. Genetic associations between the MRI measurements and each of the 21 SNPs were analyzed in linear regression models with adjustment for sample weights and traditional risk factors. Rs11628722 was tested a priori. In African Americans, rs11628722 was significantly associated with carotid wall volume (p < 0.05). Among the other 20 SNPs, adjusted for multiple testing, rs4905204, which encodes an Ala to Val amino acid change, was significantly associated with maximum wall thickness (p < 0.000625) and suggestively associated with total wall volume (p < 0.0026) in European Americans. In conclusion, SNPs in the SERPINA9 gene showed race-specific associations with characteristics of carotid atherosclerotic plaques. Replications in other populations are needed to validate findings of this study and to establish the SERPINA9 gene as a candidate in the etiology of carotid atherosclerosis.

Oral contraceptive pill use is associated with reduced odds of nonalcoholic fatty liver disease in menstruating women: results from NHANES III

BACKGROUND

Higher prevalence of nonalcoholic fatty liver disease (NAFLD) in men and postmenopausal women than in premenopausal women has suggested a potential role of sex hormones in the pathogenesis of the disease. We sought to evaluate the association between oral contraceptive pills (OCP) and NAFLD and to determine whether adiposity mediates any effect.

METHODS

We included 4338 women aged 20-60 years who were enrolled in the Third National Health and Nutrition Examination Survey from 1988 to 1994 in a population-based cross-sectional study. We defined NAFLD as moderate-severe steatosis on ultrasonography in women without excessive alcohol use or other identifiable causes. OCP use was based on self-report and was categorized as never, former or current use.

RESULTS

The overall weighted prevalence of NAFLD was 11.6 % but lower in current (6.7 %) than in former (12.0 %) or never users (15.6 %, P = 0.016). In the multivariable model, current OCP users experienced a 50 % lower odds of NAFLD than never users (adjusted odds ratio 0.50; 95 % confidence interval 0.26, 0.98) after adjusting for age, race/ethnicity, smoking status, history of diabetes or hypertension and education. Further adjustment for body mass index or waist circumference significantly attenuated the OCP-NAFLD relationship.

CONCLUSIONS

In this large US-representative population, OCP use was associated with reduced odds of NAFLD. However, this association could be mediated or confounded by adiposity. Prospective studies are needed to further clarify the causal role of sex hormone.

Meta-analysis of genome-wide association studies identifies six new Loci for serum calcium concentrations

Calcium is vital to the normal functioning of multiple organ systems and its serum concentration is tightly regulated. Apart from CASR, the genes associated with serum calcium are largely unknown. We conducted a genome-wide association meta-analysis of 39,400 individuals from 17 population-based cohorts and investigated the 14 most strongly associated loci in ≤21,679 additional individuals. Seven loci (six new regions) in association with serum calcium were identified and replicated. Rs1570669 near CYP24A1 (P = 9.1E-12), rs10491003 upstream of GATA3 (P = 4.8E-09) and rs7481584 in CARS (P = 1.2E-10) implicate regions involved in Mendelian calcemic disorders: Rs1550532 in DGKD (P = 8.2E-11), also associated with bone density, and rs7336933 near DGKH/KIAA0564 (P = 9.1E-10) are near genes that encode distinct isoforms of diacylglycerol kinase. Rs780094 is in GCKR. We characterized the expression of these genes in gut, kidney, and bone, and demonstrate modulation of gene expression in bone in response to dietary calcium in mice. Our results shed new light on the genetics of calcium homeostasis.

Calcium is vital to many biological processes and its serum concentration is tightly regulated. Family studies have shown that serum calcium is under strong genetic control. Apart from CASR, the genes associated with serum calcium are largely unknown. We conducted a genome-wide association meta-analysis of 39,400 individuals from 17 population-based cohorts and investigated the 14 most strongly associated loci in ≤21,679 additional individuals. We identified seven loci (six new regions) as being robustly associated with serum calcium. Three loci implicate regions involved in rare monogenic diseases including disturbances of serum calcium levels. Several of the newly identified loci harbor genes linked to the hormonal control of serum calcium. In mice experiments, we characterized the expression of these genes in gut, kidney, and bone, and explored the influence of dietary calcium intake on the expression of these genes in these organs. Our results shed new light on the genetics of calcium homeostasis and suggest a role for dietary calcium intake in bone-specific gene expression.

APOL1 variants associate with increased risk of CKD among African Americans

Although case-control studies suggest that African Americans with common coding variants in the APOL1 gene are 5-29 times more likely than those individuals without such variants to have focal segmental glomerulosclerosis, HIV-associated nephropathy, or ESRD, prospective studies have not yet evaluated the impact of these variants on CKD in a community-based sample of African Americans. Here, we studied whether the APOL1 G1 and G2 risk alleles associate with the development of CKD and progression to ESRD by analyzing data from 3067 African Americans in the Atherosclerosis Risk in Communities Study who did not have CKD at baseline. Carrying two risk alleles associated with a 1.49-fold increased risk of CKD (95% CI=1.02 to 2.17) and a 1.88-fold increased risk of ESRD (95% CI=1.20 to 2.93) compared with zero or one risk allele; associations persisted after adjusting for European ancestry. Among participants who developed CKD, those participants with two risk alleles were more likely to progress to ESRD than their counterparts with zero or one risk allele (HR=2.22, 95% CI=1.01 to 4.84). In conclusion, APOL1 risk variants are risk factors for the development of CKD and progression from CKD to ESRD among African Americans in the general population.

Host APOL1 genotype is independently associated with proteinuria in HIV infection

Proteinuria is associated with adverse clinical outcomes in HIV infection. Here we evaluated whether APOL1 risk alleles, previously associated with advanced kidney disease, are independently associated with proteinuria in HIV infection in a cross-sectional study of HIV-infected women in the Women's Interagency HIV Study. We estimated the percent difference in urine protein excretion and odds of proteinuria (≥200 mg/g) associated with two versus one or no APOL1 risk allele using linear and logistic regression, respectively. Of 1285 women successfully genotyped, 379 carried one and 80 carried two risk alleles. Proteinuria was present in 124 women, 78 of whom had proteinuria confirmed on a second sample. In women without prior AIDS, two risk alleles were independently associated with a 69% higher urine protein excretion (95% confidence interval (CI): 36, 108) and five-fold higher odds of proteinuria (95% CI: 2.45, 10.37) as compared with one or no risk allele. No association was found in women with prior AIDS. Analyses in which women with impaired kidney function were excluded and proteinuria was confirmed by a second urine sample yielded similar estimates. Thus, APOL1 risk alleles are associated with significant proteinuria in HIV-infected persons without prior clinical AIDS, independent of clinical factors traditionally associated with proteinuria. Trials are needed to determine whether APOL1 genotyping identifies individuals who could benefit from earlier intervention to prevent overt renal disease.

Association of functional polymorphism rs2231142 (Q141K) in the ABCG2 gene with serum uric acid and gout in 4 US populations: the PAGE Study

A loss-of-function mutation (Q141K, rs2231142) in the ATP-binding cassette, subfamily G, member 2 gene (ABCG2) has been shown to be associated with serum uric acid levels and gout in Asians, Europeans, and European and African Americans; however, less is known about these associations in other populations. Rs2231142 was genotyped in 22,734 European Americans, 9,720 African Americans, 3,849 Mexican Americans, and 3,550 American Indians in the Population Architecture using Genomics and Epidemiology (PAGE) Study (2008-2012). Rs2231142 was significantly associated with serum uric acid levels (P = 2.37 × 10(-67), P = 3.98 × 10(-5), P = 6.97 × 10(-9), and P = 5.33 × 10(-4) in European Americans, African Americans, Mexican Americans, and American Indians, respectively) and gout (P = 2.83 × 10(-10), P = 0.01, and P = 0.01 in European Americans, African Americans, and Mexican Americans, respectively). Overall, the T allele was associated with a 0.24-mg/dL increase in serum uric acid level (P = 1.37 × 10(-80)) and a 1.75-fold increase in the odds of gout (P = 1.09 × 10(-12)). The association between rs2231142 and serum uric acid was significantly stronger in men, postmenopausal women, and hormone therapy users compared with their counterparts. The association with gout was also significantly stronger in men than in women. These results highlight a possible role of sex hormones in the regulation of ABCG2 urate transporter and its potential implications for the prevention, diagnosis, and treatment of hyperuricemia and gout.

Lack of associations of ten candidate coronary heart disease risk genetic variants and subclinical atherosclerosis in four US populations: the Population Architecture using Genomics and Epidemiology (PAGE) study

BACKGROUND

A number of genetic variants have been discovered by recent genome-wide association studies for their associations with clinical coronary heart disease (CHD). However, it is unclear whether these variants are also associated with the development of CHD as measured by subclinical atherosclerosis phenotypes, ankle brachial index (ABI), carotid artery intima-media thickness (cIMT) and carotid plaque.

METHODS

Ten CHD risk single nucleotide polymorphisms (SNPs) were genotyped in individuals of European American (EA), African American (AA), American Indian (AI), and Mexican American (MA) ancestry in the Population Architecture using Genomics and Epidemiology (PAGE) study. In each individual study, we performed linear or logistic regression to examine population-specific associations between SNPs and ABI, common and internal cIMT, and plaque. The results from individual studies were meta-analyzed using a fixed effect inverse variance weighted model.

RESULTS

None of the ten SNPs was significantly associated with ABI and common or internal cIMT, after Bonferroni correction. In the sample of 13,337 EA, 3809 AA, and 5353 AI individuals with carotid plaque measurement, the GCKR SNP rs780094 was significantly associated with the presence of plaque in AI only (OR = 1.32, 95% confidence interval: 1.17, 1.49, P = 1.08 × 10(-5)), but not in the other populations (P = 0.90 in EA and P = 0.99 in AA). A 9p21 region SNP, rs1333049, was nominally associated with plaque in EA (OR = 1.07, P = 0.02) and in AI (OR = 1.10, P = 0.05).

CONCLUSIONS

We identified a significant association between rs780094 and plaque in AI populations, which needs to be replicated in future studies. There was little evidence that the index CHD risk variants identified through genome-wide association studies in EA influence the development of CHD through subclinical atherosclerosis as assessed by cIMT and ABI across ancestries.

Overlap between common genetic polymorphisms underpinning kidney traits and cardiovascular disease phenotypes: the CKDGen consortium

BACKGROUND

Chronic kidney disease is associated with cardiovascular disease. We tested for evidence of a shared genetic basis to these traits.

STUDY DESIGN

We conducted 2 targeted analyses. First, we examined whether known single-nucleotide polymorphisms (SNPs) underpinning kidney traits were associated with a series of vascular phenotypes. Additionally, we tested whether vascular SNPs were associated with markers of kidney damage. Significance was set to 1.5×10(-4) (0.05/325 tests).

SETTING & PARTICIPANTS

Vascular outcomes were analyzed in participants from the AortaGen (20,634), CARDIoGRAM (86,995), CHARGE Eye (15,358), CHARGE IMT (31,181), ICBP (69,395), and NeuroCHARGE (12,385) consortia. Tests for kidney outcomes were conducted in up to 67,093 participants from the CKDGen consortium.

PREDICTOR

We used 19 kidney SNPs and 64 vascular SNPs.

OUTCOMES & MEASUREMENTS

Vascular outcomes tested were blood pressure, coronary artery disease, carotid intima-media thickness, pulse wave velocity, retinal venular caliber, and brain white matter lesions. Kidney outcomes were estimated glomerular filtration rate and albuminuria.

RESULTS

In general, we found that kidney disease variants were not associated with vascular phenotypes (127 of 133 tests were nonsignificant). The one exception was rs653178 near SH2B3 (SH2B adaptor protein 3), which showed direction-consistent association with systolic (P = 9.3 ×10(-10)) and diastolic (P = 1.6 ×10(-14)) blood pressure and coronary artery disease (P = 2.2 ×10(-6)), all previously reported. Similarly, the 64 SNPs associated with vascular phenotypes were not associated with kidney phenotypes (187 of 192 tests were nonsignificant), with the exception of 2 high-correlated SNPs at the SH2B3 locus (P = 1.06 ×10(-07) and P = 7.05 ×10(-08)).

LIMITATIONS

The combined effect size of the SNPs for kidney and vascular outcomes may be too low to detect shared genetic associations.

CONCLUSIONS

Overall, although we confirmed one locus (SH2B3) as associated with both kidney and cardiovascular disease, our primary findings suggest that there is little overlap between kidney and cardiovascular disease risk variants in the overall population. The reciprocal risks of kidney and cardiovascular disease may not be genetically mediated, but rather a function of the disease milieu itself.

Higher magnesium intake is associated with lower fasting glucose and insulin, with no evidence of interaction with select genetic loci, in a meta-analysis of 15 CHARGE Consortium Studies

Favorable associations between magnesium intake and glycemic traits, such as fasting glucose and insulin, are observed in observational and clinical studies, but whether genetic variation affects these associations is largely unknown. We hypothesized that single nucleotide polymorphisms (SNPs) associated with either glycemic traits or magnesium metabolism affect the association between magnesium intake and fasting glucose and insulin. Fifteen studies from the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) Consortium provided data from up to 52,684 participants of European descent without known diabetes. In fixed-effects meta-analyses, we quantified 1) cross-sectional associations of dietary magnesium intake with fasting glucose (mmol/L) and insulin (ln-pmol/L) and 2) interactions between magnesium intake and SNPs related to fasting glucose (16 SNPs), insulin (2 SNPs), or magnesium (8 SNPs) on fasting glucose and insulin. After adjustment for age, sex, energy intake, BMI, and behavioral risk factors, magnesium (per 50-mg/d increment) was inversely associated with fasting glucose [β = -0.009 mmol/L (95% CI: -0.013, -0.005), P < 0.0001] and insulin [-0.020 ln-pmol/L (95% CI: -0.024, -0.017), P < 0.0001]. No magnesium-related SNP or interaction between any SNP and magnesium reached significance after correction for multiple testing. However, rs2274924 in magnesium transporter-encoding TRPM6 showed a nominal association (uncorrected P = 0.03) with glucose, and rs11558471 in SLC30A8 and rs3740393 near CNNM2 showed a nominal interaction (uncorrected, both P = 0.02) with magnesium on glucose. Consistent with other studies, a higher magnesium intake was associated with lower fasting glucose and insulin. Nominal evidence of TRPM6 influence and magnesium interaction with select loci suggests that further investigation is warranted.

Genome-wide significant locus of beta-trace protein, a novel kidney function biomarker, identified in European and African Americans

BACKGROUND

Beta-trace protein (BTP), measured in serum or plasma, has potential as a novel biomarker for kidney function. Little is known about the genes influencing BTP levels.

METHODS

We conducted a genome-wide association study of log-transformed plasma BTP levels in 6720 European Americans (EAs) and replicated the significant associations in 1734 African Americans (AAs) from the Atherosclerosis Risk in Communities (ARIC) study.

RESULTS

We identified a genome-wide significant locus in EA upstream of Prostaglandin D2 synthase (PTGDS), the gene encoding BTP. Each copy of the A allele at rs57024841 was associated with 5% higher BTP levels (P = 1.2 × 10(-23)). The association at PTGDS was confirmed in AAs (6% higher BTP for each A allele at rs57024841, P = 1.9 × 10(-7)). The index single nucleotide polymorphisms (SNPs) in EAs and AAs explained ∼1.1% of the log(BTP) variance within each population and explained over 30% of the difference in log(BTP) levels between EAs and AAs. The index SNPs at the PTGDS locus in the two populations were not associated with the estimated glomerular filtration rate (eGFR) or the urine albumin creatinine ratio (P > 0.05). We further tested for the associations of BTP with 16 known loci of the eGFR in EA, and BTP was associated with 3 of 16 tested.

CONCLUSIONS

The identification of a novel BTP-specific (non-renal related) locus and the confirmation of several genetic loci of the eGFR with BTP extend our understanding of the metabolism of BTP and inform its use as a kidney filtration biomarker.

Genome-wide association analyses identify 18 new loci associated with serum urate concentrations

Elevated serum urate concentrations can cause gout, a prevalent and painful inflammatory arthritis. By combining data from >140,000 individuals of European ancestry within the Global Urate Genetics Consortium (GUGC), we identified and replicated 28 genome-wide significant loci in association with serum urate concentrations (18 new regions in or near TRIM46, INHBB, SFMBT1, TMEM171, VEGFA, BAZ1B, PRKAG2, STC1, HNF4G, A1CF, ATXN2, UBE2Q2, IGF1R, NFAT5, MAF, HLF, ACVR1B-ACVRL1 and B3GNT4). Associations for many of the loci were of similar magnitude in individuals of non-European ancestry. We further characterized these loci for associations with gout, transcript expression and the fractional excretion of urate. Network analyses implicate the inhibins-activins signaling pathways and glucose metabolism in systemic urate control. New candidate genes for serum urate concentration highlight the importance of metabolic control of urate production and excretion, which may have implications for the treatment and prevention of gout.

DNA repair gene variants in relation to overall cancer risk: a population-based study

The hypothesis that germ-line polymorphisms in DNA repair genes influence cancer risk has previously been tested primarily on a cancer site-specific basis. The purpose of this study was to test the hypothesis that DNA repair gene allelic variants contribute to globally elevated cancer risk by measuring associations with risk of all cancers that occurred within a population-based cohort. In the CLUE II cohort study established in 1989 in Washington County, MD, this study was comprised of all 3619 cancer cases ascertained through 2007 compared with a sample of 2296 with no cancer. Associations were measured between 759 DNA repair gene single nucleotide polymorphisms (SNPs) and risk of all cancers. A SNP in O(6)-methylguanine-DNA methyltransferase, MGMT, (rs2296675) was significantly associated with overall cancer risk [per minor allele odds ratio (OR) 1.30, 95% confidence interval (CI) 1.19-1.43 and P-value: 4.1 × 10(-8)]. The association between rs2296675 and cancer risk was stronger among those aged ≤54 years old than those who were ≥55 years at baseline (P-for-(interaction) = 0.021). OR were in the direction of increased risk for all 15 categories of malignancies studied (P < 0.0001), ranging from 1.22 (P = 0.42) for ovarian cancer to 2.01 (P = 0.008) for urinary tract cancers; the smallest P-value was for breast cancer (OR 1.45, P = 0.0002). The results indicate that the minor allele of MGMT SNP rs2296675, a common genetic marker with 37% carriers, was significantly associated with increased risk of cancer across multiple tissues. Replication is needed to more definitively determine the scientific and public health significance of this observed association.

Integration of genome-wide association studies with biological knowledge identifies six novel genes related to kidney function

In conducting genome-wide association studies (GWAS), analytical approaches leveraging biological information may further understanding of the pathophysiology of clinical traits. To discover novel associations with estimated glomerular filtration rate (eGFR), a measure of kidney function, we developed a strategy for integrating prior biological knowledge into the existing GWAS data for eGFR from the CKDGen Consortium. Our strategy focuses on single nucleotide polymorphism (SNPs) in genes that are connected by functional evidence, determined by literature mining and gene ontology (GO) hierarchies, to genes near previously validated eGFR associations. It then requires association thresholds consistent with multiple testing, and finally evaluates novel candidates by independent replication. Among the samples of European ancestry, we identified a genome-wide significant SNP in FBXL20 (P = 5.6 × 10(-9)) in meta-analysis of all available data, and additional SNPs at the INHBC, LRP2, PLEKHA1, SLC3A2 and SLC7A6 genes meeting multiple-testing corrected significance for replication and overall P-values of 4.5 × 10(-4)-2.2 × 10(-7). Neither the novel PLEKHA1 nor FBXL20 associations, both further supported by association with eGFR among African Americans and with transcript abundance, would have been implicated by eGFR candidate gene approaches. LRP2, encoding the megalin receptor, was identified through connection with the previously known eGFR gene DAB2 and extends understanding of the megalin system in kidney function. These findings highlight integration of existing genome-wide association data with independent biological knowledge to uncover novel candidate eGFR associations, including candidates lacking known connections to kidney-specific pathways. The strategy may also be applicable to other clinical phenotypes, although more testing will be needed to assess its potential for discovery in general.

A population-based study of hedgehog pathway gene variants in relation to the dual risk of basal cell carcinoma plus another cancer

INTRODUCTION

A personal history of basal cell carcinoma (BCC) is associated with increased risk of other malignancies, but the reason is unknown. The hedgehog pathway is critical to the etiology of BCC, and is also believed to contribute to susceptibility to other cancers. This study tested the hypothesis that hedgehog pathway and pathway-related gene variants contribute to the increased risk of subsequent cancers among those with a history of BCC.

METHODS

The study was nested within the ongoing CLUE II cohort study, established in 1989 in Washington County, Maryland, USA. The study consisted of a cancer-free control group (n=2296) compared to three different groups of cancer cases ascertained through 2007, those diagnosed with: (1) Other (non-BCC) cancer only (n=2349); (2) BCC only (n=534); and (3) BCC plus other cancer (n=446). The frequencies of variant alleles were compared among these four groups for 20 single nucleotide polymorphisms (SNPs) in 6 hedgehog pathway genes (SHH, IHH, PTCH2, SMO, GLI1, SUFU), and also 22 SNPs in VDR and 8 SNPs in FAS, which have cross-talk with the hedgehog pathway.

RESULTS

Comparing those with both BCC and other cancer versus those with no cancer, no significant associations were observed for any of the hedgehog pathway SNPs, or for the FAS SNPs. One VDR SNP was nominally significantly associated with the BCC cancer-prone phenotype, rs11574085 [per minor allele odds ratio (OR) 1.38, 95% confidence interval (CI) 1.05-1.82; p-value=0.02].

CONCLUSION

The hedgehog pathway gene SNPs studied, along with the VDR and FAS SNPs studied, are not strongly associated with the BCC cancer-prone phenotype.

Association of estimated glomerular filtration rate and urinary uromodulin concentrations with rare variants identified by UMOD gene region sequencing

BACKGROUND

Recent genome-wide association studies (GWAS) have identified common variants in the UMOD region associated with kidney function and disease in the general population. To identify novel rare variants as well as common variants that may account for this GWAS signal, the exons and 4 kb upstream region of UMOD were sequenced.

METHODOLOGY/PRINCIPAL FINDINGS

Individuals (n = 485) were selected based on presence of the GWAS risk haplotype and chronic kidney disease (CKD) in the ARIC Study and on the extremes of of the UMOD gene product, uromodulin, in urine (Tamm Horsfall protein, THP) in the Framingham Heart Study (FHS). Targeted sequencing was conducted using capillary based Sanger sequencing (3730 DNA Analyzer). Variants were tested for association with THP concentrations and estimated glomerular filtration rate (eGFR), and identified non-synonymous coding variants were genotyped in up to 22,546 follow-up samples. Twenty-four and 63 variants were identified in the 285 ARIC and 200 FHS participants, respectively. In both studies combined, there were 33 common and 54 rare (MAF<0.05) variants. Five non-synonymous rare variants were identified in FHS; borderline enrichment of rare variants was found in the extremes of THP (SKAT p-value = 0.08). Only V458L was associated with THP in the FHS general-population validation sample (p = 9*10(-3), n = 2,522), but did not show direction-consistent and significant association with eGFR in both the ARIC (n = 14,635) and FHS (n = 7,520) validation samples. Pooling all non-synonymous rare variants except V458L together showed non-significant associations with THP and eGFR in the FHS validation sample. Functional studies of V458L revealed no alternations in protein trafficking.

CONCLUSIONS/SIGNIFICANCE

Multiple novel rare variants in the UMOD region were identified, but none were consistently associated with eGFR in two independent study samples. Only V458L had modest association with THP levels in the general population and thus could not account for the observed GWAS signal.

A population-based study of DNA repair gene variants in relation to non-melanoma skin cancer as a marker of a cancer-prone phenotype

For unknown reasons, non-melanoma skin cancer (NMSC) is associated with increased risk of other malignancies. Focusing solely on DNA repair or DNA repair-related genes, this study tested the hypothesis that DNA repair gene variants contribute to the increased cancer risk associated with a personal history of NMSC. From the parent CLUE II cohort study, established in 1989 in Washington County, MD, the study consisted of a cancer-free control group (n 5 2296) compared with three mutually exclusive groups of cancer cases ascertained through 2007: (i) Other (non-NMSC) cancer only (n 5 2349); (ii) NMSC only (n 5 694) and (iii) NMSC plus other cancer (n 5 577). The frequency of minor alleles in 759 DNA repair gene single nucleotide polymorphisms (SNPs) was compared in these four groups. Comparing those with both NMSC and other cancer versus those with no cancer, 10 SNPs had allelic trend P-values <0.01. The two top-ranked SNPs were both within the thymine DNA glycosylase gene (TDG). One was a non-synonymous coding SNP (rs2888805) [per allele odds ratio (OR) 1.40, 95% confidence interval (CI) 1.16-1.70; P-value 5 0.0006] and the other was an intronic SNP in high linkage disequilibrium with rs2888805 (rs4135150). None of the associations had a P-value <6.6310(-5), the threshold for statistical significance after correcting for multiple comparisons. The results pinpoint DNA repair genes most likely to contribute to the NMSC cancer-prone phenotype. A promising lead is genetic variants in TDG, important not only in base excision repair but also in regulating the epigenome and gene expression, which may contribute to the NMSC-associated increase in overall cancer risk.

No interactions between previously associated 2-hour glucose gene variants and physical activity or BMI on 2-hour glucose levels

Gene-lifestyle interactions have been suggested to contribute to the development of type 2 diabetes. Glucose levels 2 h after a standard 75-g glucose challenge are used to diagnose diabetes and are associated with both genetic and lifestyle factors. However, whether these factors interact to determine 2-h glucose levels is unknown. We meta-analyzed single nucleotide polymorphism (SNP) × BMI and SNP × physical activity (PA) interaction regression models for five SNPs previously associated with 2-h glucose levels from up to 22 studies comprising 54,884 individuals without diabetes. PA levels were dichotomized, with individuals below the first quintile classified as inactive (20%) and the remainder as active (80%). BMI was considered a continuous trait. Inactive individuals had higher 2-h glucose levels than active individuals (β = 0.22 mmol/L [95% CI 0.13-0.31], P = 1.63 × 10(-6)). All SNPs were associated with 2-h glucose (β = 0.06-0.12 mmol/allele, P ≤ 1.53 × 10(-7)), but no significant interactions were found with PA (P > 0.18) or BMI (P ≥ 0.04). In this large study of gene-lifestyle interaction, we observed no interactions between genetic and lifestyle factors, both of which were associated with 2-h glucose. It is perhaps unlikely that top loci from genome-wide association studies will exhibit strong subgroup-specific effects, and may not, therefore, make the best candidates for the study of interactions.

Meta-analysis identifies six new susceptibility loci for atrial fibrillation

Atrial fibrillation is a highly prevalent arrhythmia and a major risk factor for stroke, heart failure and death. We conducted a genome-wide association study (GWAS) in individuals of European ancestry, including 6,707 with and 52,426 without atrial fibrillation. Six new atrial fibrillation susceptibility loci were identified and replicated in an additional sample of individuals of European ancestry, including 5,381 subjects with and 10,030 subjects without atrial fibrillation (P < 5 × 10(-8)). Four of the loci identified in Europeans were further replicated in silico in a GWAS of Japanese individuals, including 843 individuals with and 3,350 individuals without atrial fibrillation. The identified loci implicate candidate genes that encode transcription factors related to cardiopulmonary development, cardiac-expressed ion channels and cell signaling molecules.

Genome-wide association and functional follow-up reveals new loci for kidney function

Chronic kidney disease (CKD) is an important public health problem with a genetic component. We performed genome-wide association studies in up to 130,600 European ancestry participants overall, and stratified for key CKD risk factors. We uncovered 6 new loci in association with estimated glomerular filtration rate (eGFR), the primary clinical measure of CKD, in or near MPPED2, DDX1, SLC47A1, CDK12, CASP9, and INO80. Morpholino knockdown of mpped2 and casp9 in zebrafish embryos revealed podocyte and tubular abnormalities with altered dextran clearance, suggesting a role for these genes in renal function. By providing new insights into genes that regulate renal function, these results could further our understanding of the pathogenesis of CKD.

African ancestry and its correlation to type 2 diabetes in African Americans: a genetic admixture analysis in three U.S. population cohorts

The risk of type 2 diabetes is approximately 2-fold higher in African Americans than in European Americans even after adjusting for known environmental risk factors, including socioeconomic status (SES), suggesting that genetic factors may explain some of this population difference in disease risk. However, relatively few genetic studies have examined this hypothesis in a large sample of African Americans with and without diabetes. Therefore, we performed an admixture analysis using 2,189 ancestry-informative markers in 7,021 African Americans (2,373 with type 2 diabetes and 4,648 without) from the Atherosclerosis Risk in Communities Study, the Jackson Heart Study, and the Multiethnic Cohort to 1) determine the association of type 2 diabetes and its related quantitative traits with African ancestry controlling for measures of SES and 2) identify genetic loci for type 2 diabetes through a genome-wide admixture mapping scan. The median percentage of African ancestry of diabetic participants was slightly greater than that of non-diabetic participants (study-adjusted difference = 1.6%, P<0.001). The odds ratio for diabetes comparing participants in the highest vs. lowest tertile of African ancestry was 1.33 (95% confidence interval 1.13–1.55), after adjustment for age, sex, study, body mass index (BMI), and SES. Admixture scans identified two potential loci for diabetes at 12p13.31 (LOD = 4.0) and 13q14.3 (Z score = 4.5, P = 6.6×10⁻⁶). In conclusion, genetic ancestry has a significant association with type 2 diabetes above and beyond its association with non-genetic risk factors for type 2 diabetes in African Americans, but no single gene with a major effect is sufficient to explain a large portion of the observed population difference in risk of diabetes. There undoubtedly is a complex interplay among specific genetic loci and non-genetic factors, which may both be associated with overall admixture, leading to the observed ethnic differences in diabetes risk.