Aims/hypothesis

Kidney injury molecule 1 (KIM-1), liver fatty acid-binding protein (L-FABP), N-acetyl-β-D-glucosaminidase (NAG) and neutrophil gelatinase-associated lipocalin (NGAL) are urinary biomarkers of renal tubular injury. We examined their association with incident end-stage renal disease (ESRD) and all-cause mortality in American Indians with type 2 diabetes.

Methods

Biomarker concentrations were measured in baseline urine samples in 260 Pima Indians who were followed for a median of 14 years. HRs were reported per SD of creatinine (Cr)-normalised log-transformed KIM-1, NAG and NGAL, and for three categories of L-FABP.

Results

During follow-up, 74 participants developed ESRD and 101 died. Median concentrations of KIM-1/Cr, NAG/Cr and NGAL/Cr and the proportion of detectable L-FABP were highest in those with macroalbuminuria (p < 0.001 for KIM-1/Cr, NAG/Cr and L-FABP; p = 0.006 for NGAL/Cr). After multivariable adjustment, NGAL/Cr was positively associated with ESRD (HR 1.59, 95% CI 1.20, 2.11) and mortality (HR 1.39, 95% CI 1.06, 1.82); L-FABP/Cr was inversely associated with ESRD (HR [for highest vs lowest tertile] 0.40, 95% CI 0.19, 0.83). Addition of NGAL/Cr to models that included albuminuria and glomerular filtration rate increased the c-statistic for predicting ESRD from 0.828 to 0.833 (p = 0.001) and for death from 0.710 to 0.722 (p = 0.018). Addition of L-FABP/Cr increased the c-statistic for ESRD from 0.828 to 0.832 (p = 0.042).

Conclusions/interpretation

In Pima Indians with type 2 diabetes, urinary concentrations of NGAL and L-FABP are associated with important health outcomes, but they are unlikely to add to risk prediction with standard markers in a clinically meaningful way given the small increase in the c-statistic.

Genome-wide significant loci for metformin response in type 2 diabetes reported elsewhere have not been replicated in the Diabetes Prevention Program (DPP). To assess pharmacogenetic interactions in prediabetes, we conducted a genome-wide association study (GWAS) in the DPP. Cox proportional hazards models tested associations with diabetes incidence in the metformin (MET; n = 876) and placebo (PBO; n = 887) arms. Multiple linear regression assessed association with 1-year change in metformin-related quantitative traits, adjusted for baseline trait, age, sex, and 10 ancestry principal components. We tested for gene-by-treatment interaction. No significant associations emerged for diabetes incidence. We identified four genome-wide significant variants after correcting for correlated traits (P < 9 × 10-9). In the MET arm, rs144322333 near ENOSF1 (minor allele frequency [MAF]AFR = 0.07; MAFEUR = 0.002) was associated with an increase in percentage of glycated hemoglobin (per minor allele, β = 0.39 [95% CI 0.28, 0.50]; P = 2.8 × 10-12). rs145591055 near OMSR (MAF = 0.10 in American Indians) was associated with weight loss (kilograms) (per G allele, β = -7.55 [95% CI -9.88, -5.22]; P = 3.2 × 10-10) in the MET arm. Neither variant was significant in PBO; gene-by-treatment interaction was significant for both variants [P(G×T) < 1.0 × 10-4]. Replication in individuals with diabetes did not yield significant findings. A GWAS for metformin response in prediabetes revealed novel ethnic-specific associations that require further investigation but may have implications for tailored therapy.

Platelets play an essential role in hemostasis and thrombosis. We performed a genome-wide association study of platelet count in 12,491 participants of the Hispanic Community Health Study/Study of Latinos by using a mixed-model method that accounts for admixture and family relationships. We discovered and replicated associations with five genes (ACTN1, ETV7, GABBR1-MOG, MEF2C, and ZBTB9-BAK1). Our strongest association was with Amerindian-specific variant rs117672662 (p value = 1.16 × 10(-28)) in ACTN1, a gene implicated in congenital macrothrombocytopenia. rs117672662 exhibited allelic differences in transcriptional activity and protein binding in hematopoietic cells. Our results underscore the value of diverse populations to extend insights into the allelic architecture of complex traits.

Background

The presence of population structure in a sample may confound the search for important genetic loci associated with disease. Our four samples in the Family Investigation of Nephropathy and Diabetes (FIND), European Americans, Mexican Americans, African Americans, and American Indians are part of a genome- wide association study in which population structure might be particularly important. We therefore decided to study in detail one component of this, individual genetic ancestry (IGA). From SNPs present on the Affymetrix 6.0 Human SNP array, we identified 3 sets of ancestry informative markers (AIMs), each maximized for the information in one the three contrasts among ancestral populations: Europeans (HAPMAP, CEU), Africans (HAPMAP, YRI and LWK), and Native Americans (full heritage Pima Indians). We estimate IGA and present an algorithm for their standard errors, compare IGA to principal components, emphasize the importance of balancing information in the ancestry informative markers (AIMs), and test the association of IGA with diabetic nephropathy in the combined sample.

Results

A fixed parental allele maximum likelihood algorithm was applied to the FIND to estimate IGA in four samples: 869 American Indians; 1385 African Americans; 1451 Mexican Americans; and 826 European Americans. When the information in the AIMs is unbalanced, the estimates are incorrect with large error. Individual genetic admixture is highly correlated with principle components for capturing population structure. It takes ~700 SNPs to reduce the average standard error of individual admixture below 0.01. When the samples are combined, the resulting population structure creates associations between IGA and diabetic nephropathy.

Conclusions

The identified set of AIMs, which include American Indian parental allele frequencies, may be particularly useful for estimating genetic admixture in populations from the Americas. Failure to balance information in maximum likelihood, poly-ancestry models creates biased estimates of individual admixture with large error. This also occurs when estimating IGA using the Bayesian clustering method as implemented in the program STRUCTURE. Odds ratios for the associations of IGA with disease are consistent with what is known about the incidence and prevalence of diabetic nephropathy in these populations.

Background

Genetic factors that influence kidney traits have been understudied for low frequency and ancestry-specific variants.

Methods

We combined whole genome sequencing (WGS) data from 23,732 participants from 10 NHLBI Trans-Omics for Precision Medicine (TOPMed) Program multi-ethnic studies to identify novel loci for estimated glomerular filtration rate (eGFR). Participants included European, African, East Asian, and Hispanic ancestries. We applied linear mixed models using a genetic relationship matrix estimated from the WGS data and adjusted for age, sex, study, and ethnicity.

Findings

When testing single variants, we identified three novel loci driven by low frequency variants more commonly observed in non-European ancestry (PRKAA2, rs180996919, minor allele frequency [MAF] 0.04%, P = 6.1 × 10^-11; METTL8, rs116951054, MAF 0.09%, P = 4.5 × 10^-9; and MATK, rs539182790, MAF 0.05%, P = 3.4 × 10^-9). We also replicated two known loci for common variants (rs2461702, MAF=0.49, P = 1.2 × 10^-9, nearest gene GATM, and rs71147340, MAF=0.34, P = 3.3 × 10^-9, CDK12). Testing aggregated variants within a gene identified the MAF gene. A statistical approach based on local ancestry helped to identify replication samples for ancestry-specific variants.

Interpretation

This study highlights challenges in studying variants influencing kidney traits that are low frequency in populations and more common in non-European ancestry.

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.

Disease incidences increase with age, but the molecular characteristics of ageing that lead to increased disease susceptibility remain inadequately understood. Here we perform a whole-blood gene expression meta-analysis in 14,983 individuals of European ancestry (including replication) and identify 1,497 genes that are differentially expressed with chronological age. The age-associated genes do not harbor more age-associated CpG-methylation sites than other genes, but are instead enriched for the presence of potentially functional CpG-methylation sites in enhancer and insulator regions that associate with both chronological age and gene expression levels. We further used the gene expression profiles to calculate the 'transcriptomic age' of an individual, and show that differences between transcriptomic age and chronological age are associated with biological features linked to ageing, such as blood pressure, cholesterol levels, fasting glucose, and body mass index. The transcriptomic prediction model adds biological relevance and complements existing epigenetic prediction models, and can be used by others to calculate transcriptomic age in external cohorts.

Coronary artery disease (CAD) is more frequent among individuals with dysglycemia. Preventive interventions for diabetes can improve cardiometabolic risk factors (CRFs), but it is unclear whether the benefits on CRFs are similar for individuals at different genetic risk for CAD. We built a 201-variant polygenic risk score (PRS) for CAD and tested for interaction with diabetes prevention strategies on 1-year changes in CRFs in 2,658 Diabetes Prevention Program (DPP) participants. We also examined whether separate lifestyle behaviors interact with PRS and affect changes in CRFs in each intervention group. Participants in both the lifestyle and metformin interventions had greater improvement in the majority of recognized CRFs compared with placebo (P < 0.001) irrespective of CAD genetic risk (P _interaction > 0.05). We detected nominal significant interactions between PRS and dietary quality and physical activity on 1-year change in BMI, fasting glucose, triglycerides, and HDL cholesterol in individuals randomized to metformin or placebo, but none of them achieved the multiple-testing correction for significance. This study confirms that diabetes preventive interventions improve CRFs regardless of CAD genetic risk and delivers hypothesis-generating data on the varying benefit of increasing physical activity and improving diet on intermediate cardiovascular risk factors depending on individual CAD genetic risk profile.