Integrated analysis of blood DNA methylation, genetic variants, circulating proteins, microRNAs, and kidney failure in type 1 diabetes.

Variation in DNA methylation (DNAmet) in white blood cells and other cells/tissues has been implicated in the etiology of progressive diabetic kidney disease (DKD). However, the specific mechanisms linking DNAmet variation in blood cells with risk of kidney failure (KF) and utility of measuring blood cell DNAmet in personalized medicine are not clear. We measured blood cell DNAmet in 277 individuals with type 1 diabetes and DKD using Illumina EPIC arrays; 51% of the cohort developed KF during 7 to 20 years of follow-up. Our epigenome-wide analysis identified DNAmet at 17 CpGs (5'-cytosine-phosphate-guanine-3' loci) associated with risk of KF independent of major clinical risk factors. DNAmet at these KF-associated CpGs remained stable over a median period of 4.7 years. Furthermore, DNAmet variations at seven KF-associated CpGs were strongly associated with multiple genetic variants at seven genomic regions, suggesting a strong genetic influence on DNAmet. The effects of DNAmet variations at the KF-associated CpGs on risk of KF were partially mediated by multiple KF-associated circulating proteins and KF-associated circulating miRNAs. A prediction model for risk of KF was developed by adding blood cell DNAmet at eight selected KF-associated CpGs to the clinical model. This updated model significantly improved prediction performance (c-statistic = 0.93) versus the clinical model (c-statistic = 0.85) at P = 6.62 × 10-14. In conclusion, our multiomics study provides insights into mechanisms through which variation of DNAmet may affect KF development and shows that blood cell DNAmet at certain CpGs can improve risk prediction for KF in T1D.

Neuroblastoma suppressor of tumorigenicity 1 is a circulating protein associated with progression to end-stage kidney disease in diabetes.

Circulating proteins associated with transforming growth factor-ß (TGF-ß) signaling are implicated in the development of diabetic kidney disease (DKD). It remains to be comprehensively examined which of these proteins are involved in the pathogenesis of DKD and its progression to end-stage kidney disease (ESKD) in humans. Using the SOMAscan proteomic platform, we measured concentrations of 25 TGF-ß signaling family proteins in four different cohorts composed in total of 754 Caucasian or Pima Indian individuals with type 1 or type 2 diabetes. Of these 25 circulating proteins, we identified neuroblastoma suppressor of tumorigenicity 1 (NBL1, aliases DAN and DAND1), a small secreted protein known to inhibit members of the bone morphogenic protein family, to be most strongly and independently associated with progression to ESKD during 10-year follow-up in all cohorts. The extent of damage to podocytes and other glomerular structures measured morphometrically in 105 research kidney biopsies correlated strongly with circulating NBL1 concentrations. Also, in vitro exposure to NBL1 induced apoptosis in podocytes. In conclusion, circulating NBL1 may be involved in the disease process underlying progression to ESKD, and its concentration in circulation may identify subjects with diabetes at increased risk of progression to ESKD.

Results of untargeted analysis using the SOMAscan proteomics platform indicates novel associations of circulating proteins with risk of progression to kidney failure in diabetes.

This study applies a large proteomics panel to search for new circulating biomarkers associated with progression to kidney failure in individuals with diabetic kidney disease. Four independent cohorts encompassing 754 individuals with type 1 and type 2 diabetes and early and late diabetic kidney disease were followed to ascertain progression to kidney failure. During ten years of follow-up, 227 of 754 individuals progressed to kidney failure. Using the SOMAscan proteomics platform, we measured baseline concentration of 1129 circulating proteins. In our previous publications, we analyzed 334 of these proteins that were members of specific candidate pathways involved in diabetic kidney disease and found 35 proteins strongly associated with risk of progression to kidney failure. Here, we examined the remaining 795 proteins using an untargeted approach. Of these remaining proteins, 11 were significantly associated with progression to kidney failure. Biological processes previously reported for these proteins were related to neuron development (DLL1, MATN2, NRX1B, KLK8, RTN4R and ROR1) and were implicated in the development of kidney fibrosis (LAYN, DLL1, MAPK11, MATN2, endostatin, and ROR1) in cellular and animal studies. Specific mechanisms that underlie involvement of these proteins in progression of diabetic kidney disease must be further investigated to assess their value as targets for kidney-protective therapies. Using multivariable LASSO regression analysis, five proteins (LAYN, ESAM, DLL1, MAPK11 and endostatin) were found independently associated with risk of progression to kidney failure. Thus, our study identified proteins that may be considered as new candidate prognostic biomarkers to predict risk of progression to kidney failure in diabetic kidney disease. Furthermore, three of these proteins (DLL1, ESAM, and MAPK11) were selected as candidate biomarkers when all SOMAscan results were evaluated.

Circulating proteins protect against renal decline and progression to end-stage renal disease in patients with diabetes.

Diabetic kidney disease (DKD) and its major clinical manifestation, progressive renal decline that leads to end-stage renal disease (ESRD), are a major health burden for individuals with diabetes. The disease process that underlies progressive renal decline comprises factors that increase risk as well as factors that protect against this outcome. Using untargeted proteomic profiling of circulating proteins from individuals in two independent cohorts with type 1 and type 2 diabetes and varying stages of DKD followed for 7 to 15 years, we identified three elevated plasma proteins-fibroblast growth factor 20 (OR, 0.69; 95% CI, 0.54 to 0.88), angiopoietin-1 (OR, 0.72; 95% CI, 0.57 to 0.91), and tumor necrosis factor ligand superfamily member 12 (OR, 0.75; 95% CI, 0.59 to 0.95)-that were associated with protection against progressive renal decline and progression to ESRD. The combined effect of these three protective proteins was demonstrated by very low cumulative risk of ESRD in those who had baseline concentrations above median for all three proteins, whereas the cumulative risk of ESRD was high in those with concentrations below median for these proteins at the beginning of follow-up. This protective effect was shown to be independent from circulating inflammatory proteins and clinical covariates and was confirmed in a third cohort of diabetic individuals with normal renal function. These three protective proteins may serve as biomarkers to stratify diabetic individuals according to risk of progression to ESRD and might also be investigated as potential therapeutics to delay or prevent the onset of ESRD.

Comprehensive Search for Novel Circulating miRNAs and Axon Guidance Pathway Proteins Associated with Risk of ESKD in Diabetes.

BACKGROUND: Mechanisms underlying the pro gression of diabetic kidney disease to ESKD are not fully understood. METHODS: We performed global microRNA (miRNA) analysis on plasma from two cohorts consisting of 375 individuals with type 1 and type 2 diabetes with late diabetic kidney disease, and targeted proteomics analysis on plasma from four cohorts consisting of 746 individuals with late and early diabetic kidney disease. We examined structural lesions in kidney biopsy specimens from the 105 individuals with early diabetic kidney disease. Human umbilical vein endothelial cells were used to assess the effects of miRNA mimics or inhibitors on regulation of candidate proteins. RESULTS: In the late diabetic kidney disease cohorts, we identified 17 circulating miRNAs, represented by four exemplars (miR-1287-5p, miR-197-5p, miR-339-5p, and miR-328-3p), that were strongly associated with 10-year risk of ESKD. These miRNAs targeted proteins in the axon guidance pathway. Circulating levels of six of these proteins-most notably, EFNA4 and EPHA2-were strongly associated with 10-year risk of ESKD in all cohorts. Furthermore, circulating levels of these proteins correlated with severity of structural lesions in kidney biopsy specimens. In contrast, expression levels of genes encoding these proteins had no apparent effects on the lesions. In in vitro experiments, mimics of miR-1287-5p and miR-197-5p and inhibitors of miR-339-5p and miR-328-3p upregulated concentrations of EPHA2 in either cell lysate, supernatant, or both. CONCLUSIONS: This study reveals novel mechanisms involved in progression to ESKD and points to the importance of systemic factors in the development of diabetic kidney disease. Some circulating miRNAs and axon guidance pathway proteins represent potential targets for new therapies to prevent and treat this condition.

Effect of TNFα stimulation on expression of kidney risk inflammatory proteins in human umbilical vein endothelial cells cultured in hyperglycemia.

We recently identified a kidney risk inflammatory signature (KRIS), comprising 6 TNF receptors (including TNFR1 and TNFR2) and 11 inflammatory proteins. Elevated levels of these proteins in circulation were strongly associated with risk of the development of end-stage kidney disease (ESKD) during 10-year follow-up. It has been hypothesized that elevated levels of these proteins in circulation might reflect (be markers of) systemic exposure to TNFα. In this in vitro study, we examined intracellular and extracellular levels of these proteins in human umbilical vein endothelial cells (HUVECs) exposed to TNFα in the presence of hyperglycemia. KRIS proteins as well as 1300 other proteins were measured using the SOMAscan proteomics platform. Four KRIS proteins (including TNFR1) were down-regulated and only 1 protein (IL18R1) was up-regulated in the extracellular fraction of TNFα-stimulated HUVECs. In the intracellular fraction, one KRIS protein was down-regulated (CCL14) and 1 protein was up-regulated (IL18R1). The levels of other KRIS proteins were not affected by exposure to TNFα. HUVECs exposed to a hyperglycemic and inflammatory environment also showed significant up-regulation of a distinct set of 53 proteins (mainly in extracellular fraction). In our previous study, circulating levels of these proteins were not associated with progression to ESKD in diabetes.

Circulating short and medium chain fatty acids are associated with normoalbuminuria in type 1 diabetes of long duration.

A substantial number of subjects with Type 1 Diabetes (T1D) of long duration never develop albuminuria or renal function impairment, yet the underlying protective mechanisms remain unknown. Therefore, our study included 308 Joslin Kidney Study subjects who had T1D of long duration (median: 24 years), maintained normal renal function and had either normoalbuminuria or a broad range of albuminuria within the 2 years preceding the metabolomic determinations. Serum samples were subjected to global metabolomic profiling. 352 metabolites were detected in at least 80% of the study population. In the logistic analyses adjusted for multiple testing (Bonferroni corrected α = 0.000028), we identified 38 metabolites associated with persistent normoalbuminuria independently from clinical covariates. Protective metabolites were enriched in Medium Chain Fatty Acids (MCFAs) and in Short Chain Fatty Acids (SCFAs) and particularly involved odd-numbered and dicarboxylate Fatty Acids. One quartile change of nonanoate, the top protective MCFA, was associated with high odds of having persistent normoalbuminuria (OR (95% CI) 0.14 (0.09, 0.23); p < 10-12). Multivariable Random Forest analysis concordantly indicated to MCFAs as effective classifiers. Associations of the relevant Fatty Acids with albuminuria seemed to parallel associations with tubular biomarkers. Our findings suggest that MCFAs and SCFAs contribute to the metabolic processes underlying protection against albuminuria development in T1D that are independent from mechanisms associated with changes in renal function.

A profile of multiple circulating tumor necrosis factor receptors associated with early progressive kidney decline in Type 1 Diabetes is similar to profiles in autoimmune disorders.

This study comprehensively evaluated the association between known circulating tumor necrosis factor (TNF) superfamily ligands and receptors and the development of early progressive kidney decline (PKD) leading to end-stage kidney disease (ESKD) in Type 1 diabetes. Participants for the study were from the Macro-Albuminuria Study (198 individuals), and the Micro-Albuminuria Study (148 individuals) of the Joslin Kidney Study. All individuals initially had normal kidney function and were followed for seven-fifteen years to determine the slope of the estimate glomerular filtration rate and to ascertain onset of ESKD. Plasma concentrations of 25 TNF superfamily proteins were measured using proximity extension assay applied in the OLINK proteomics platform. In the both studies risk of early PKD, determined as estimated glomerular filtration rate loss greater than or equal to three ml/min/1.73m2/year, was associated with elevated circulating levels of 13 of 19 TNF receptors examined. In the Macro-Albuminuria Study, we obtained similar findings for risk of progression to ESKD. These receptors comprised: TNF-R1A, -R1B, -R3, -R4, -R6, -R6B, -R7, -R10A, -R10B, -R11A, -R14, -R21, and -R27. Serial measurements showed that circulating levels of these TNF receptors had increased before the onset of PKD. In contrast, none of the six measured TNF ligands showed association with risk of early PKD. Of significance, the disease process that underlies PKD leading to ESKD in Type 1 diabetes has a profile also seen in autoimmune disorders. The mechanisms of this enrichment may be causally related to the development of PKD in Type 1 diabetes and must be investigated further. Thus, some of these receptors may be used as new risk predictors of ESKD.

Profibrotic Circulating Proteins and Risk of Early Progressive Renal Decline in Patients With Type 2 Diabetes With and Without Albuminuria.

OBJECTIVE: The role of fibrosis in early progressive renal decline in type 2 diabetes is unknown. Circulating WFDC2 (WAP four-disulfide core domain protein 2) and matrix metalloproteinase 7 (MMP-7; Matrilysin) are postulated to be biomarkers of renal fibrosis. This study examined an association of circulating levels of these proteins with early progressive renal decline. RESEARCH DESIGN AND METHODS: Individuals with type 2 diabetes enrolled in the Joslin Kidney Study with an estimated glomerular filtration rate (eGFR) ≥60 mL/min/1.73 m2 were monitored for 6-12 years to ascertain fast early progressive renal decline, defined as eGFR loss ≥5 mL/min/1.73 m2/year. RESULTS: A total of 1,181 individuals were studied: 681 without and 500 with albuminuria. Median eGFR and albumin-to-creatinine ratio (ACR) at baseline were 97 mL/min/1.73 m2 and 24 mg/g, respectively. During follow-up, 152 individuals experienced fast early progressive renal decline: 6.9% in those with normoalbuminuria and 21% with albuminuria. In both subgroups, the risk of renal decline increased with increasing baseline levels of WFDC2 (P < 0.0001) and MMP-7 (P < 0.0001). After adjustment for relevant clinical characteristics and known biomarkers, an increase by one quartile in the fibrosis index (combination of levels of WFDC2 and MMP-7) was associated with higher risk of renal decline (odds ratio 1.63; 95% CI 1.30-2.04). The association was similar and statistically significant among patients with and without albuminuria. CONCLUSIONS: Elevation of circulating profibrotic proteins is associated with the development of early progressive renal decline in type 2 diabetes. This association is independent from albuminuria status and points to the importance of the fibrotic process in the development of early renal decline.

A signature of circulating inflammatory proteins and development of end-stage renal disease in diabetes.

Chronic inflammation is postulated to be involved in the development of end-stage renal disease in diabetes, but which specific circulating inflammatory proteins contribute to this risk remain unknown. To study this, we examined 194 circulating inflammatory proteins in subjects from three independent cohorts with type 1 and type 2 diabetes. In each cohort, we identified an extremely robust kidney risk inflammatory signature (KRIS), consisting of 17 proteins enriched in tumor necrosis factor-receptor superfamily members, that was associated with a 10-year risk of end-stage renal disease. All these proteins had a systemic, non-kidney source. Our prospective study findings provide strong evidence that KRIS proteins contribute to the inflammatory process underlying end-stage renal disease development in both types of diabetes. These proteins point to new therapeutic targets and new prognostic tests to identify subjects at risk of end-stage renal disease, as well as biomarkers to measure responses to treatment of diabetic kidney disease.

Mycophenolic acid concentrations in peripheral blood mononuclear cells are associated with the incidence of rejection in renal transplant recipients.

AIMS: Although therapeutic drug monitoring of plasma mycophenolic acid (MPA) concentrations has been recommended to individualize dosage in transplant recipients, little is known regarding lymphocyte concentrations of MPA, where MPA inhibits inosine monophosphate dehydrogenase (IMPDH). This study investigated the utility of measuring predose MPA concentrations in peripheral blood mononuclear cells (C0C ) and predose IMPDH activity, as predictors of graft rejection in renal transplant recipients. METHODS: Forty-eight patients commencing mycophenolate mofetil (1 g twice daily) in combination with tacrolimus and prednisolone were recruited. Blood was collected for determination of trough total (C0P ) and unbound (C0u ) plasma MPA concentrations. Peripheral blood mononuclear cells were isolated for determination of C0C and IMPDH activity. The incidence of rejection within 2 days of sample collection was determined histologically and classified according to the Banff 2007 criteria. RESULTS: There was no association between MPA C0C and C0P (rs  = 0.28, P = 0.06), however, MPA C0C were weakly correlated with MPA C0u (rs  = 0.42, P = 0.013). Multivariate analysis indicated that MPA C0C was the only covariate independently associated with rejection (FDR-adjusted P = 0.033). The receiver operating characteristic area under the curve (AUC) for the prediction of severe rejection using MPA C0C was 0.75 (P = 0.013), with 73% sensitivity and specificity at a C0C threshold of 0.5 ng 10-7 cells. However, predose IMPDH activity was not a predictor of rejection (P > 0.15). CONCLUSIONS: MPA C0C measurement within the early post-transplant period may be useful to facilitate early titration of MPA dosing to significantly reduce rejection.

Circulating miRNA Profiles Associated With Hyperglycemia in Patients With Type 1 Diabetes.

We investigated plasma microRNA (miRNA) profiles associated with variation of hyperglycemia, measured as hemoglobin A1c (HbA1c), in two panels of patients with type 1 diabetes (T1D). Using the HTG Molecular Diagnostics EdgeSeq platform, 2,083 miRNAs were measured in plasma from 71 patients included in a screening panel. Quantitative real-time PCR was used to measure the candidate miRNAs in plasma from 95 patients included in an independent replication panel. We found 10 miRNAs replicated in both panels and 4 with high statistical significance. The strongest positive correlations with HbA1c were found with miR-125b-5p (rs = 0.40, P = 6.0 × 10-5) and miR-365a-3p (rs = 0.35, P = 5.9 × 10-4). The strongest negative correlations were found with miR-5190 (rs = -0.30, P = 0.003) and miR-770-5p (rs = -0.27, P = 0.008). Pathway analysis revealed that 50 Kyoto Encyclopedia of Genes and Genomes pathways were significantly enriched by genes targeted by these four miRNAs. The axon guidance signaling pathway was enriched (P < 1 × 10-7) by genes targeted by all four miRNAs. In addition, three other pathways (Rap1 signaling, focal adhesion, and neurotrophin signaling) were also significantly enriched but with genes targeted by only by three of the identified miRNAs. In conclusion, our study identified four circulating miRNAs that were influenced by variation in hyperglycemia. Dysregulation of these miRNAs, which are associated with hyperglycemia in patients with T1D, may contribute to the development of diabetes complications. However, there are multitudes of possible mechanisms/pathways through which dysregulation of these miRNAs may impact risk of diabetes complications.

Validation of an LC-MS/MS method for the quantification of mycophenolic acid in human kidney transplant biopsies.

Mycophenolic acid (MPA) has a low therapeutic index and large inter-individual pharmacokinetic variability necessitating therapeutic drug monitoring to individualise dosing after transplantation. There is an ongoing discrepancy as to whether plasma MPA concentrations sufficiently predict kidney rejection or toxicity and whether immunosuppressant concentrations within the graft tissue may better predict transplant outcomes. The aim of the study was to develop an LC-MS/MS method for the quantification of MPA concentrations in human kidney biopsies taken as part of routine clinical procedures. A total of 4 surplus human kidney biopsies obtained from 4 different kidney transplant recipients were available to use for this study. MPA was also quantified in 2 kidney samples from rats administered MPA to assess tissue extraction reproducibility. Human kidney biopsies and rat kidneys were homogenised mechanically and underwent liquid-liquid extraction before analysis by LC-MS/MS. MPA-free human kidney tissue was used in calibrators and quality control samples. Analyte detection was achieved from multiple reaction monitoring of the ammonium adducts of both MPA (m/z 321.1→207.3) and N-phthaloyl-l-phenylalanine (PPA, internal standard, m/z 296.2→250.2) using positive electrospray ionisation. The method was linear (calibration curves R(2)>0.99, n=10), precise, and accurate with coefficients of variation and bias less than 15%. Extraction efficiencies for MPA and PPA were approximately 97% and 86%, respectively, and matrix effects were minimal. In 4 kidney transplant recipients, tissue MPA concentrations ranged from 1.3 to 7.7ng/mg of tissue, however, the correlation between blood (C0) and tissue MPA concentrations could not be established. The method was successfully applied to the quantification of MPA in human kidney biopsies without the need to alter current clinical protocols.