Whole-exome sequencing study identifies four novel gene loci associated with diabetic kidney disease.

Diabetic kidney disease (DKD) is recognized as an important public health challenge. However, its genomic mechanisms are poorly understood. To identify rare variants for DKD, we conducted a whole-exome sequencing (WES) study leveraging large cohorts well-phenotyped for chronic kidney disease and diabetes. Our two-stage WES study included 4372 European and African ancestry participants from the Chronic Renal Insufficiency Cohort and Atherosclerosis Risk in Communities studies (stage 1) and 11 487 multi-ancestry Trans-Omics for Precision Medicine participants (stage 2). Generalized linear mixed models, which accounted for genetic relatedness and adjusted for age, sex and ancestry, were used to test associations between single variants and DKD. Gene-based aggregate rare variant analyses were conducted using an optimized sequence kernel association test implemented within our mixed model framework. We identified four novel exome-wide significant DKD-related loci through initiating diabetes. In single-variant analyses, participants carrying a rare, in-frame insertion in the DIS3L2 gene (rs141560952) exhibited a 193-fold increased odds [95% confidence interval (CI): 33.6, 1105] of DKD compared with noncarriers (P = 3.59 × 10-9). Likewise, each copy of a low-frequency KRT6B splice-site variant (rs425827) conferred a 5.31-fold higher odds (95% CI: 3.06, 9.21) of DKD (P = 2.72 × 10-9). Aggregate gene-based analyses further identified ERAP2 (P = 4.03 × 10-8) and NPEPPS (P = 1.51 × 10-7), which are both expressed in the kidney and implicated in renin-angiotensin-aldosterone system modulated immune response. In the largest WES study of DKD, we identified novel rare variant loci attaining exome-wide significance. These findings provide new insights into the molecular mechanisms underlying DKD.

Proteomics Analysis of Genetic Liability of Abdominal Aortic Aneurysm Identifies Plasma Neogenin and Kit Ligand: The ARIC Study.

BACKGROUND: Genome-wide association studies have reported 23 gene loci related to abdominal aortic aneurysm (AAA)-a potentially lethal condition characterized by a weakened dilated vessel wall. This study aimed to identify proteomic signatures and pathways related to these risk loci to better characterize AAA genetic susceptibility. METHODS: Plasma concentrations of 4870 proteins were determined using a DNA aptamer-based array. Linear regression analysis estimated the associations between the 23 risk alleles and plasma protein levels with adjustments for potential confounders in a race-stratified analysis of 1671 Black and 7241 White participants. Significant proteins were then evaluated for their prediction of clinical AAA (454 AAA events in 11 064 individuals), and those significantly associated with AAA were further interrogated using Mendelian randomization analysis. RESULTS: Risk variants proximal to PSRC1-CELSR2-SORT1, PCIF1-ZNF335-MMP9, RP11-136O12.2/TRIB1, ZNF259/APOA5, IL6R, PCSK9, LPA, and APOE were associated with 118 plasma proteins in Whites and 59 were replicated in Black participants. Novel associations with clinical AAA incidence were observed for kit ligand (HR, 0.59 [95% CI, 0.42-0.82] for top versus first quintiles) and neogenin (HR, 0.64 [95% CI, 0.46-0.88]) over a median 21.2-year follow-up; neogenin was also associated with ultrasound-detected asymptomatic AAA (N=4295; 57 asymptomatic AAA cases). Mendelian randomization inverse variance weighted estimates suggested that AAA risk is promoted by lower levels of kit ligand (OR per SD=0.67; P=1.4×10-5) and neogenin (OR per SD=0.50; P=0.03). CONCLUSIONS: Low levels of neogenin and kit ligand may be novel risk factors for AAA development in potentially causal pathways. These findings provide insights and potential targets to reduce AAA susceptibility.

Comparison of Proteomic Measurements Across Platforms in the Atherosclerosis Risk in Communities (ARIC) Study.

BACKGROUND: The plasma proteome can be quantified using different types of highly multiplexed technologies, including aptamer-based and proximity-extension immunoassay methods. There has been limited characterization of how these protein measurements correlate across platforms and with absolute measures from targeted immunoassays. METHODS: We assessed the comparability of (a) highly multiplexed aptamer-based (SomaScan v4; Somalogic) and proximity-extension immunoassay (OLINK Proseek® v5003; Olink) methods in 427 Atherosclerosis Risk in Communities (ARIC) Study participants (Visit 5, 2011-2013), and (b) 18 of the SomaScan protein measurements against targeted immunoassays in 110 participants (55 cardiovascular disease cases, 55 controls). We calculated Spearman correlations (r) between the different measurements and compared associations with case-control status. RESULTS: There were 417 protein comparisons (366 unique proteins) between the SomaScan and Olink platforms. The average correlation was r = 0.46 (range: -0.21 to 0.97; 79 [19%] with r ≥ 0.8). For the comparison of SomaScan and targeted immunoassays, 6 of 18 assays (growth differentiation factor 15 [GDF15], interleukin-1 receptor-like 1 [ST2], interstitial collagenase [MMP1], adiponectin, leptin, and resistin) had good correlations (r ≥ 0.8), 2 had modest correlations (0.5 ≤ r < 0.8; osteopontin and interleukin-6 [IL6]), and 10 were poorly correlated (r < 0.5; metalloproteinase inhibitor 1 [TIMP1], stromelysin-1 [MMP3], matrilysin [MMP7], C-C motif chemokine 2 [MCP1], interleukin-10 [IL10], vascular cell adhesion protein 1 [VCAM1], intercellular adhesion molecule 1 [ICAM1], interleukin-18 [IL18], tumor necrosis factor [TNFα], and visfatin) overall. Correlations for SomaScan and targeted immunoassays were similar according to case status. CONCLUSIONS: There is variation in the quantitative measurements for many proteins across aptamer-based and proximity-extension immunoassays (approximately 1/2 showing good or modest correlation and approximately 1/2 poor correlation) and also for correlations of these highly multiplexed technologies with targeted immunoassays. Design and interpretation of protein quantification studies should be informed by the variation across measurement techniques for each protein.

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

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

Mendelian Randomization Analysis as a Tool to Gain Insights into Causes of Diseases: A Primer.

Many Mendelian randomization (MR) studies have been published recently, with inferences on the causal relationships between risk factors and diseases that have potential implications for clinical research. In nephrology, MR methods have been applied to investigate potential causal relationships of traditional risk factors, lifestyle factors, and biomarkers from omics technologies with kidney function or CKD. This primer summarizes the basic concepts of MR studies, highlighting methods used in recent applications, and emphasizes key elements in conducting and reporting of MR studies that are important for interpreting the results.

Polygenic Risk Scores for Kidney Function and Their Associations with Circulating Proteome, and Incident Kidney Diseases.

BACKGROUND: Genome-wide association studies (GWAS) have revealed numerous loci for kidney function (eGFR). The relationship between polygenic predictors of eGFR, risk of incident adverse kidney outcomes, and the plasma proteome is not known. METHODS: We developed a genome-wide polygenic risk score (PRS) for eGFR by applying the LDpred algorithm to summary statistics generated from a multiethnic meta-analysis of CKDGen Consortium GWAS ( n =765,348) and UK Biobank GWAS (90% of the cohort; n =451,508), followed by best-parameter selection using the remaining 10% of UK Biobank data ( n =45,158). We then tested the association of the PRS in the Atherosclerosis Risk in Communities (ARIC) study ( n =8866) with incident CKD, ESKD, kidney failure, and AKI. We also examined associations between the PRS and 4877 plasma proteins measured at middle age and older adulthood and evaluated mediation of PRS associations by eGFR. RESULTS: The developed PRS showed a significant association with all outcomes. Hazard ratios per 1 SD lower PRS ranged from 1.06 (95% CI, 1.01 to 1.11) to 1.33 (95% CI, 1.28 to 1.37). The PRS was significantly associated with 132 proteins at both time points. The strongest associations were with cystatin C, collagen α -1(XV) chain, and desmocollin-2. Most proteins were higher at lower kidney function, except for five proteins, including testican-2. Most correlations of the genetic PRS with proteins were mediated by eGFR. CONCLUSIONS: A PRS for eGFR is now sufficiently strong to capture risk for a spectrum of incident kidney diseases and broadly influences the plasma proteome, primarily mediated by eGFR.

Genome-wide association study of serum metabolites in the African American Study of Kidney Disease and Hypertension.

The genome-wide association study (GWAS) is a powerful means to study genetic determinants of disease traits and generate insights into disease pathophysiology. To date, few GWAS of circulating metabolite levels have been performed in African Americans with chronic kidney disease. Hypothesizing that novel genetic-metabolite associations may be identified in a unique population of African Americans with a lower glomerular filtration rate (GFR), we conducted a GWAS of 652 serum metabolites in 619 participants (mean measured glomerular filtration rate 45 mL/min/1.73m2) in the African American Study of Kidney Disease and Hypertension, a clinical trial of blood pressure lowering and antihypertensive medication in African Americans with chronic kidney disease. We identified 42 significant variant metabolite associations. Twenty associations had been previously identified in published GWAS, and eleven novel associations were replicated in a separate cohort of 818 African Americans with genetic and metabolomic data from the Atherosclerosis Risk in Communities Study. The replicated novel variant-metabolite associations comprised eight metabolites and eleven distinct genomic loci. Nine of the replicated associations represented clear enzyme-metabolite interactions, with high expression in the kidneys as well as the liver. Three loci (ACY1, ACY3, and NAT8) were associated with a common pool of metabolites, acetylated amino acids, but with different individual affinities. Thus, extensive metabolite profiling in an African American population with chronic kidney disease aided identification of novel genome-wide metabolite associations, providing clues about substrate specificity and the key roles of enzymes in modulating systemic levels of metabolites.

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

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

Association between Circulating Protein C Levels and Incident Dementia: The Atherosclerosis Risk in Communities Study.

INTRODUCTION: Hemostasis depends on the delicate balance between coagulants and anticoagulants. Higher levels of circulating coagulants have been associated with higher risk of cerebral infarctions and dementia. In contrast, higher levels of circulating protein C, an endogenous anticoagulant, have been associated with lower risk of cerebral infarctions, and the association between protein C levels and the risk of dementia is unknown. The goal of this study was to evaluate the association of circulating protein C levels in midlife and late life with incident dementia. METHODS: Circulating protein C levels were measured using blood samples collected at the midlife baseline (1987-1989) and the late-life baseline (2011-2013) among 14,462 and 3,614 participants, respectively, in the Atherosclerosis Risk in Communities study. Protein C levels were measured using enzyme-linked immunosorbent assay at midlife and a modified aptamer-based assay at late life. Participants were followed up to 2013 from midlife and up to 2017 from late life. Incident dementia was ascertained during the follow-up periods using in-person cognitive and functional assessment, informant interviews, and International Classification of Diseases codes at hospitalization discharge and on death certificates. Cause-specific Cox regression models were used to evaluate the association between quintiles of circulating protein C and incident dementia. RESULTS: From midlife (mean age of 54), 1,389 incident dementia events were observed over a median follow-up of 23 years. From late life (mean age of 75), 353 incident dementia events were observed over a median follow-up of 4.9 years. At both midlife and late life, circulating protein C had an inverse association with incident dementia after adjusting for demographic, vascular, and hemostatic risk factors, incident stroke as time-dependent covariate, and incorporating stabilized weights based on propensity scores (quintile 5 vs. quintile 1 as the reference, midlife hazard ratio 0.80, 95% confidence interval 0.66-0.96, p value for trend 0.04; late-life hazard ratio 0.84, 95% confidence interval: 0.55-1.28, p value for trend 0.04). DISCUSSION/CONCLUSION: Circulating protein C has an inverse association with incident dementia independent of established risk factors, including stroke. Our results suggest studying anticoagulants in addition to coagulants can increase our understanding on the relationship between hemostasis and dementia.

Integration of GWAS Summary Statistics and Gene Expression Reveals Target Cell Types Underlying Kidney Function Traits.

BACKGROUND: Genetic variants identified in genome-wide association studies (GWAS) are often not specific enough to reveal complex underlying physiology. By integrating RNA-seq data and GWAS summary statistics, novel computational methods allow unbiased identification of trait-relevant tissues and cell types. METHODS: The CKDGen consortium provided GWAS summary data for eGFR, urinary albumin-creatinine ratio (UACR), BUN, and serum urate. Genotype-Tissue Expression Project (GTEx) RNA-seq data were used to construct the top 10% specifically expressed genes for each of 53 tissues followed by linkage disequilibrium (LD) score-based enrichment testing for each trait. Similar procedures were performed for five kidney single-cell RNA-seq datasets from humans and mice and for a microdissected tubule RNA-seq dataset from rat. Gene set enrichment analyses were also conducted for genes implicated in Mendelian kidney diseases. RESULTS: Across 53 tissues, genes in kidney function-associated GWAS loci were enriched in kidney (P=9.1E-8 for eGFR; P=1.2E-5 for urate) and liver (P=6.8·10-5 for eGFR). In the kidney, proximal tubule was enriched in humans (P=8.5E-5 for eGFR; P=7.8E-6 for urate) and mice (P=0.0003 for eGFR; P=0.0002 for urate) and confirmed as the primary cell type in microdissected tubules and organoids. Gene set enrichment analysis supported this and showed enrichment of genes implicated in monogenic glomerular diseases in podocytes. A systematic approach generated a comprehensive list of GWAS genes prioritized by cell type-specific expression. CONCLUSIONS: Integration of GWAS statistics of kidney function traits and gene expression data identified relevant tissues and cell types, as a basis for further mechanistic studies to understand GWAS loci.

GSTM1 Deletion Exaggerates Kidney Injury in Experimental Mouse Models and Confers the Protective Effect of Cruciferous Vegetables in Mice and Humans.

BACKGROUND: GSTM1 encodes glutathione S-transferase µ-1 (GSTM1), which belongs to a superfamily of phase 2 antioxidant enzymes. The highly prevalent GSTM1 deletion variant is associated with kidney disease progression in human cohorts: the African American Study of Kidney Disease and Hypertension and the Atherosclerosis Risk in Communities (ARIC) Study. METHODS: We generated a Gstm1 knockout mouse line to study its role in a CKD model (involving subtotal nephrectomy) and a hypertension model (induced by angiotensin II). We examined the effect of intake of cruciferous vegetables and GSTM1 genotypes on kidney disease in mice as well as in human ARIC study participants. We also examined the importance of superoxide in the mediating pathways and of hematopoietic GSTM1 on renal inflammation. RESULTS: Gstm1 knockout mice displayed increased oxidative stress, kidney injury, and inflammation in both models. The central mechanism for kidney injury is likely mediated by oxidative stress, because treatment with Tempol, an superoxide dismutase mimetic, rescued kidney injury in knockout mice without lowering BP. Bone marrow crosstransplantation revealed that Gstm1 deletion in the parenchyma, and not in bone marrow-derived cells, drives renal inflammation. Furthermore, supplementation with cruciferous broccoli powder rich in the precursor to antioxidant-activating sulforaphane significantly ameliorated kidney injury in Gstm1 knockout, but not wild-type mice. Similarly, among humans (ARIC study participants), high consumption of cruciferous vegetables was associated with fewer kidney failure events compared with low consumption, but this association was observed primarily in participants homozygous for the GSTM1 deletion variant. CONCLUSIONS: Our data support a role for the GSTM1 enzyme in the modulation of oxidative stress, inflammation, and protective metabolites in CKD.

Heritability analysis of nontraditional glycemic biomarkers in the Atherosclerosis Risk in Communities Study.

Nontraditional glycemic biomarkers, including fructosamine, glycated albumin, and 1,5-anhydroglucitol (1,5-AG) are potential alternatives or complement to traditional measures of hyperglycemia. Genetic variants are associated with these biomarkers, but the heritability, or extent to which genetics control their variation, is not known. We estimated pedigree-based, SNP-based, and bivariate heritabilities for traditional glycemic biomarkers (fasting glucose, HbA1c), and nontraditional biomarkers (fructosamine, glycated albumin, 1,5-AG) among white participants in the Atherosclerosis Risk in Communities (ARIC) Study (N = 400 first-degree relatives from sibships, N = 5,575 unrelated individuals). Pedigree-based heritabilities (representing heritability from the entire genome) for nontraditional biomarkers were substantial (0.44-0.55) and comparable to HbA1c (0.34); the fasting glucose estimate was nonsignificant. SNP-based heritabilities (representing heritability from common variants) were lower than pedigree-based heritabilities for all biomarkers. Bivariate heritabilities showed shared genetics between fructosamine and glycated albumin (0.46 pedigree-based, 1.00 SNP-based) and glycated albumin and 1,5-AG (0.50 pedigree-based, 0.47 SNP-based). Genetic factors contribute to a considerable proportion of the variance of fructosamine, glycated albumin, and 1,5-AG and a portion of this heritability likely comes from common variants.

A bidirectional Mendelian randomization study supports causal effects of kidney function on blood pressure.

Blood pressure and kidney function have a bidirectional relation. Hypertension has long been considered as a risk factor for kidney function decline. However, whether intensive blood pressure control could promote kidney health has been uncertain. The kidney is known to have a major role in affecting blood pressure through sodium extraction and regulating electrolyte balance. This bidirectional relation makes causal inference between these two traits difficult. Therefore, to examine the causal relations between these two traits, we performed two-sample Mendelian randomization analyses using summary statistics of large-scale genome-wide association studies. We selected genetic instruments more likely to be specific for kidney function using meta-analyses of complementary kidney function biomarkers (glomerular filtration rate estimated from serum creatinine [eGFRcr], and blood urea nitrogen from the CKDGen Consortium). Systolic and diastolic blood pressure summary statistics were from the International Consortium for Blood Pressure and UK Biobank. Significant evidence supported the causal effects of higher kidney function on lower blood pressure. Based on the mode-based Mendelian randomization method, the effect estimates for one standard deviation (SD) higher in log-transformed eGFRcr was -0.17 SD unit (95 % confidence interval: -0.09 to -0.24) in systolic blood pressure and -0.15 SD unit (95% confidence interval: -0.07 to -0.22) in diastolic blood pressure. In contrast, the causal effects of blood pressure on kidney function were not statistically significant. Thus, our results support causal effects of higher kidney function on lower blood pressure and suggest preventing kidney function decline can reduce the public health burden of hypertension.

Genome-wide association study identifies novel loci for type 2 diabetes-attributed end-stage kidney disease in African Americans.

BACKGROUND: End-stage kidney disease (ESKD) is a significant public health concern disproportionately affecting African Americans (AAs). Type 2 diabetes (T2D) is the leading cause of ESKD in the USA, and efforts to uncover genetic susceptibility to diabetic kidney disease (DKD) have had limited success. A prior genome-wide association study (GWAS) in AAs with T2D-ESKD was expanded with additional AA cases and controls and genotypes imputed to the higher density 1000 Genomes reference panel. The discovery analysis included 3432 T2D-ESKD cases and 6977 non-diabetic non-nephropathy controls (N = 10,409), followed by a discrimination analysis in 2756 T2D non-nephropathy controls to exclude T2D-associated variants. RESULTS: Six independent variants located in or near RND3/RBM43, SLITRK3, ENPP7, GNG7, and APOL1 achieved genome-wide significant association (P < 5 × 10-8) with T2D-ESKD. Following extension analyses in 1910 non-diabetic ESKD cases and 908 non-diabetic non-nephropathy controls, a meta-analysis of 5342 AA all-cause ESKD cases and 6977 AA non-diabetic non-nephropathy controls revealed an additional novel all-cause ESKD locus at EFNB2 (rs77113398; P = 9.84 × 10-9; OR = 1.94). Exclusion of APOL1 renal-risk genotype carriers identified two additional genome-wide significant T2D-ESKD-associated loci at GRAMD3 and MGAT4C. A second variant at GNG7 (rs373971520; P = 2.17 × 10-8, OR = 1.46) remained associated with all-cause ESKD in the APOL1-negative analysis. CONCLUSIONS: Findings provide further evidence for genetic factors associated with advanced kidney disease in AAs with T2D.

APOL1 Kidney Risk Variants and Cardiovascular Disease: An Individual Participant Data Meta-Analysis.

BACKGROUND: Two coding variants in the apo L1 gene (APOL1) are strongly associated with kidney disease in blacks. Kidney disease itself increases the risk of cardiovascular disease, but whether these variants have an independent direct effect on the risk of cardiovascular disease is unclear. Previous studies have had inconsistent results. METHODS: We conducted a two-stage individual participant data meta-analysis to assess the association of APOL1 kidney-risk variants with adjudicated cardiovascular disease events and death, independent of kidney measures. The analysis included 21,305 blacks from eight large cohorts. RESULTS: Over 8.9±5.0 years of follow-up, 2076 incident cardiovascular disease events occurred in the 16,216 participants who did not have cardiovascular disease at study enrollment. In fully-adjusted analyses, individuals possessing two APOL1 kidney-risk variants had similar risk of incident cardiovascular disease (coronary heart disease, myocardial infarction, stroke and heart failure; hazard ratio 1.11, 95% confidence interval, 0.96 to 1.28) compared to individuals with zero or one kidney-risk variant. The risk of coronary heart disease, myocardial infarction, stroke and heart failure considered individually was also comparable by APOL1 genotype. APOL1 genotype was also not associated with death. There was no difference in adjusted associations by level of kidney function, age, diabetes status, or body-mass index. CONCLUSIONS: In this large, two-stage individual participant data meta-analysis, APOL1 kidney-risk variants were not associated with incident cardiovascular disease or death independent of kidney measures.

Serum 6-Bromotryptophan Levels Identified as a Risk Factor for CKD Progression.

Background Metabolite levels reflect physiologic homeostasis and may serve as biomarkers of disease progression. Identifying metabolites associated with APOL1 risk alleles-genetic variants associated with CKD risk commonly present in persons of African descent-may reveal novel markers of CKD progression relevant to other populations.Methods We evaluated associations between the number of APOL1 risk alleles and 760 serum metabolites identified via untargeted profiling in participants of the African American Study of Kidney Disease and Hypertension (AASK) (n=588; Bonferroni significance threshold P<6.5×10-5) and replicated findings in 678 black participants with CKD in BioMe, an electronic medical record-linked biobank. We tested the metabolite association with CKD progression in AASK, BioMe, and the Modification of Diet in Renal Disease (MDRD) Study.Results One metabolite, 6-bromotryptophan, was significant in AASK (P=4.7×10-5) and replicated in BioMe (P=5.7×10-3) participants, with lower levels associated with more APOL1 risk alleles. Lower levels of 6-bromotryptophan were associated with CKD progression in AASK and BioMe participants and in white participants in the MDRD Study, independent of demographics and clinical characteristics, including baseline GFR (adjusted hazard ratio per two-fold higher 6-bromotryptophan level, AASK, 0.76; 95% confidence interval [95% CI], 0.64 to 0.91; BioMe, 0.61; 95% CI, 0.43 to 0.85; MDRD, 0.52; 95% CI, 0.34 to 0.79). The interaction between the APOL1 risk alleles and 6-bromotryptophan was not significant. The identity of 6-bromotryptophan was confirmed in experiments comparing its molecular signature with that of authentic standards of other bromotryptophan isomers.Conclusions Serum 6-bromotryptophan is a consistent and novel risk factor for CKD progression.

Serum metabolomic profile of incident diabetes.

AIMS/HYPOTHESIS: Metabolomic profiling offers the potential to reveal metabolic pathways relevant to the pathophysiology of diabetes and improve diabetes risk prediction. METHODS: We prospectively analysed known metabolites using an untargeted approach in serum specimens from baseline (1987-1989) and incident diabetes through to 31 December 2015 in a subset of 2939 Atherosclerosis Risk in Communities (ARIC) study participants with metabolomics data and without prevalent diabetes. RESULTS: Among the 245 named compounds identified, seven metabolites were significantly associated with incident diabetes after Bonferroni correction and covariate adjustment; these included a food additive (erythritol) and compounds involved in amino acid metabolism [isoleucine, leucine, valine, asparagine, 3-(4-hydoxyphenyl)lactate] and glucose metabolism (trehalose). Higher levels of metabolites were associated with increased risk of incident diabetes (HR per 1 SD increase in isoleucine 2.96, 95% CI 2.02, 4.35, p = 3.18 × 10-8; HR per 1 SD increase in trehalose 1.16, 95% CI 1.09, 1.25, p = 1.87 × 10-5), with the exception of asparagine, which was associated with a lower risk of diabetes (HR per 1 SD increase in asparagine 0.78, 95% CI 0.71, 0.85, p = 4.19 × 10-8). The seven metabolites modestly improved prediction of incident diabetes beyond fasting glucose and established risk factors (C statistics 0.744 vs 0.735, p = 0.001 for the difference in C statistics). CONCLUSIONS/INTERPRETATION: Branched chain amino acids may play a role in diabetes development. Our study is the first to report asparagine as a protective biomarker of diabetes risk. The serum metabolome reflects known and novel metabolic disturbances that improve prediction of diabetes.

Serum metabolites are associated with all-cause mortality in chronic kidney disease.

Chronic kidney disease (CKD) involves significant metabolic abnormalities and has a high mortality rate. Because the levels of serum metabolites in patients with CKD might provide insight into subclinical disease states and risk for future mortality, we determined which serum metabolites reproducibly associate with mortality in CKD using a discovery and replication design. Metabolite levels were quantified via untargeted liquid chromatography and mass spectroscopy from serum samples of 299 patients with CKD in the Modification of Diet in Renal Disease (MDRD) study as a discovery cohort. Six among 622 metabolites were significantly associated with mortality over a median follow-up of 17 years after adjustment for demographic and clinical covariates, including urine protein and measured glomerular filtration rate. We then replicated associations with mortality in 963 patients with CKD from the African American Study of Kidney Disease and Hypertension (AASK) cohort over a median follow-up of ten years. Three of the six metabolites identified in the MDRD cohort replicated in the AASK cohort: fumarate, allantoin, and ribonate, belonging to energy, nucleotide, and carbohydrate pathways, respectively. Point estimates were similar in both studies and in meta-analysis (adjusted hazard ratios 1.63, 1.59, and 1.61, respectively, per doubling of the metabolite). Thus, selected serum metabolites were reproducibly associated with long-term mortality in CKD beyond markers of kidney function in two well characterized cohorts, providing targets for investigation.

APOL1 Risk Variants and Cardiovascular Disease: Results From the AASK (African American Study of Kidney Disease and Hypertension).

OBJECTIVE: Among African Americans, the apolipoprotein L1 (APOL1) risk variants have been associated with various types of kidney disease and chronic kidney disease progression. We aimed to determine whether these same risk variants also confer an increased risk for cardiovascular disease. APPROACH AND RESULTS: In a cohort of African Americans with hypertension-attributed chronic kidney disease followed for up to 12 years, we used Cox proportional hazards models to estimate the relative hazard of a composite cardiovascular disease outcome (cardiovascular death or hospitalization for myocardial infarction, cardiac revascularization procedure, heart failure, or stroke) for the APOL1 high- (2 risk variants) versus low-risk (0-1 risk variant) genotypes. We adjusted for age, sex, ancestry, smoking, heart disease history, body mass index, cholesterol, randomized treatment groups, and baseline and longitudinal estimated glomerular filtration rate, systolic blood pressure, and proteinuria. Among 693 participants with APOL1 genotyping available (23% high risk), the high-risk group had lower mean estimated glomerular filtration rate (44.7 versus 50.1 mL/min per 1.73 m2) and greater proteinuria (median 0.19 versus 0.06) compared with the low-risk group at baseline. There was no significant association between APOL1 genotypes and the composite cardiovascular disease outcome in both unadjusted (hazard ratio=1.23; 95% confidence interval: 0.83-1.81) and fully adjusted (hazard ratio=1.16; 95% confidence interval: 0.77-1.76) models; however, in using an additive model, APOL1 high-risk variants were associated with increased cardiovascular mortality. CONCLUSIONS: Among African Americans with hypertension-attributed chronic kidney disease, APOL1 risk variants were not associated with an overall risk for cardiovascular disease although some signals for cardiovascular mortality were noted.

The Loss of GSTM1 Associates with Kidney Failure and Heart Failure.

Glutathione S-transferase mu 1 (GSTM1) encodes an enzyme that catalyzes the conjugation of electrophilic compounds with glutathione to facilitate their degradation or excretion. The loss of one or both copies of GSTM1 is common in many populations and has been associated with CKD progression. With the hypothesis that the loss of GSTM1 is also associated with incident kidney failure and heart failure, we estimated GSTM1 copy number using exome sequencing reads in the Atherosclerosis Risk in Communities (ARIC) Study, a community-based prospective cohort of white and black participants. Overall, 51.2% and 39.8% of white participants and 25.6% and 48.5% of black participants had zero or one copy of GSTM1, respectively. Over a median follow-up of 24.6 years, 256 kidney failure events occurred in 5715 participants without prevalent kidney failure, and 1028 heart failure events occurred in 5368 participants without prevalent heart failure. In analysis adjusted for demographics, diabetes, and hypertension, having zero or one copy of GSTM1 associated with higher risk of kidney failure and heart failure (adjusted hazard ratio [95% confidence interval] for zero or one versus two copies of GSTM1: kidney failure, 1.66 [1.27 to 2.17]; heart failure, 1.16 [1.04 to 1.29]). Risk did not differ significantly between participants with zero and one copy of GSTM1 (P>0.10). In summary, the loss of GSTM1 was significantly associated with incident kidney and heart failure, independent of traditional risk factors. These results suggest GSTM1 function is a potential treatment target for the prevention of kidney and heart failure.