Unbiased kidney-centric molecular categorization of chronic kidney disease as a step towards precision medicine.

Current classification of chronic kidney disease (CKD) into stages using indirect systemic measures (estimated glomerular filtration rate (eGFR) and albuminuria) is agnostic to the heterogeneity of underlying molecular processes in the kidney thereby limiting precision medicine approaches. To generate a novel CKD categorization that directly reflects within kidney disease drivers we analyzed publicly available transcriptomic data from kidney biopsy tissue. A Self-Organizing Maps unsupervised artificial neural network machine-learning algorithm was used to stratify a total of 369 patients with CKD and 46 living kidney donors as healthy controls. Unbiased stratification of the discovery cohort resulted in identification of four novel molecular categories of disease termed CKD-Blue, CKD-Gold, CKD-Olive, CKD-Plum that were replicated in independent CKD and diabetic kidney disease datasets and can be further tested on any external data at kidneyclass.org. Each molecular category spanned across CKD stages and histopathological diagnoses and represented transcriptional activation of distinct biological pathways. Disease progression rates were highly significantly different between the molecular categories. CKD-Gold displayed rapid progression, with significant eGFR-adjusted Cox regression hazard ratio of 5.6 [1.01-31.3] for kidney failure and hazard ratio of 4.7 [1.3-16.5] for composite of kidney failure or a 40% or more eGFR decline. Urine proteomics revealed distinct patterns between the molecular categories, and a 25-protein signature was identified to distinguish CKD-Gold from other molecular categories. Thus, patient stratification based on kidney tissue omics offers a gateway to non-invasive biomarker-driven categorization and the potential for future clinical implementation, as a key step towards precision medicine in CKD.

Interactive Effects of Empagliflozin and Hyperglycemia on Urinary Amino Acids in Individuals With Type 1 Diabetes.

Optimizing energy use in the kidney is critical for normal kidney function. Here, we investigate the effect of hyperglycemia and sodium-glucose cotransporter 2 (SGLT2) inhibition on urinary amino acid excretion in individuals with type 1 diabetes (T1D). The open-label ATIRMA trial assessed the impact of 8 weeks of 25 mg empagliflozin orally once per day in 40 normotensive normoalbuminuric young adults with T1D. A consecutive 2-day assessment of clamped euglycemia and hyperglycemia was evaluated at baseline and posttreatment visits. Principal component analysis was performed on urinary amino acids grouped into representative metabolic pathways using MetaboAnalyst. At baseline, acute hyperglycemia was associated with changes in 25 of the 33 urinary amino acids or their metabolites. The most significant amino acid metabolites affected by acute hyperglycemia were 3-hydroxykynurenine, serotonin, glycyl-histidine, and nicotinic acid. The changes in amino acid metabolites were reflected by the induction of four biosynthetic pathways: aminoacyl-tRNA; valine, leucine, and isoleucine; arginine; and phenylalanine, tyrosine, and tryptophan. In acute hyperglycemia, empagliflozin significantly attenuated the increases in aminoacyl-tRNA biosynthesis and valine, leucine, and isoleucine biosynthesis. Our findings using amino acid metabolomics indicate that hyperglycemia stimulates biosynthetic pathways in T1D. SGLT2 inhibition may attenuate the increase in biosynthetic pathways to optimize kidney energy metabolism.

Tubular dysfunction impairs renal excretion of pseudouridine in diabetic kidney disease.

Plasma nucleosides-pseudouridine (PU) and N2N2-dimethyl guanosine (DMG) predict the progression of type 2 diabetic kidney disease (DKD) to end-stage renal disease, but the mechanisms underlying this relationship are not well understood. We used a well-characterized model of type 2 diabetes (db/db mice) and control nondiabetic mice (db/m mice) to characterize the production and excretion of PU and DMG levels using liquid chromatography-mass spectrometry. The fractional excretion of PU and DMG was decreased in db/db mice compared with control mice at 24 wk before any changes to renal function. We then examined the dynamic changes in nucleoside metabolism using in vivo metabolic flux analysis with the injection of labeled nucleoside precursors. Metabolic flux analysis revealed significant decreases in the ratio of urine-to-plasma labeling of PU and DMG in db/db mice compared with db/m mice, indicating significant tubular dysfunction in diabetic kidney disease. We observed that the gene and protein expression of the renal tubular transporters involved with nucleoside transport in diabetic kidneys in mice and humans was reduced. In conclusion, this study strongly suggests that tubular handling of nucleosides is altered in early DKD, in part explaining the association of PU and DMG with human DKD progression observed in previous studies.NEW & NOTEWORTHY Tubular dysfunction explains the association between the nucleosides pseudouridine and N2N2-dimethyl guanosine and diabetic kidney disease.

Sodium glucose co-transporter 2 inhibition increases epidermal growth factor expression and improves outcomes in patients with type 2 diabetes.

Underlying molecular mechanisms of the kidney protective effects of sodium glucose co-transporter 2 (SGLT2) inhibitors are not fully elucidated. Therefore, we studied the association between urinary epidermal growth factor (uEGF), a mitogenic factor involved in kidney repair, and kidney outcomes in patients with type 2 diabetes (T2D). The underlying molecular mechanisms of the SGLT2 inhibitor canagliflozin on EGF using single-cell RNA sequencing from kidney tissue were examined. Urinary EGF-to-creatinine ratio (uEGF/Cr) was measured in 3521 CANagliflozin cardioVascular Assessment Study (CANVAS) participants at baseline and week 52. Associations of uEGF/Cr with kidney outcome were assessed using multivariable-adjusted Cox regression models. Single-cell RNA sequencing was performed using protocol kidney biopsy tissue from ten young patients with T2D on SGLT2i, six patients with T2D on standard care only, and six healthy controls (HCs). In CANVAS, each doubling in baseline uEGF/Cr was associated with a 12% (95% confidence interval 1-22) decreased risk of kidney outcome. uEGF/Cr decreased after 52 weeks with placebo and remained stable with canagliflozin (between-group difference +7.3% (2.0-12.8). In young persons with T2D, EGF mRNA was primarily expressed in the thick ascending loop of Henle. Expression in biopsies from T2D without SGLT2i was significantly lower compared to HCs, whereas treatment with SGLT2i increased EGF levels closer to the healthy state. In young persons with T2D without SGLT2i, endothelin-1 emerged as a key regulator of the EGF co-expression network. SGLT2i treatment was associated with a shift towards normal EGF expression. Thus, decreased uEGF represents increased risk of kidney disease progression in patients with T2D. Canagliflozin increased kidney tissue expression of EGF and was associated with a downstream signaling cascade linked to tubular repair and reversal of tubular injury.

Defining the molecular correlate of arteriolar hyalinosis in kidney disease progression by integration of single cell transcriptomic analysis and pathology scoring.

Arteriolar hyalinosis in kidneys is an independent predictor of cardiovascular disease, the main cause of mortality in chronic kidney disease (CKD). The underlying molecular mechanisms of protein accumulation in the subendothelial space are not well understood. Using single cell transcriptomic data and whole slide images from kidney biopsies of patients with CKD and acute kidney injury in the Kidney Precision Medicine Project, the molecular signals associated with arteriolar hyalinosis were evaluated. Co-expression network analysis of the endothelial genes yielded three gene set modules as significantly associated with arteriolar hyalinosis. Pathway analysis of these modules showed enrichment of transforming growth factor beta / bone morphogenetic protein (TGFß / BMP) and vascular endothelial growth factor (VEGF) signaling pathways in the endothelial cell signatures. Ligand-receptor analysis identified multiple integrins and cell adhesion receptors as over-expressed in arteriolar hyalinosis, suggesting a potential role of integrin-mediated TGFß signaling. Further analysis of arteriolar hyalinosis associated endothelial module genes identified focal segmental glomerular sclerosis as an enriched term. On validation in gene expression profiles from the Nephrotic Syndrome Study Network cohort, one of the three modules was significantly associated with the composite endpoint (> 40% reduction in estimated glomerular filtration rate (eGFR) or kidney failure) independent of age, sex, race, and baseline eGFR, suggesting poor prognosis with elevated expression of genes in this module. Thus, integration of structural and single cell molecular features yielded biologically relevant gene sets, signaling pathways and ligand-receptor interactions, underlying arteriolar hyalinosis and putative targets for therapeutic intervention.

Intraindividual variations of urinary biomarkers in hospitalized children with glomerular diseases: a prospective observational study.

This study aimed to assess the intraindividual variations of urinary biomarkers in hospitalized children with glomerular diseases. Hospitalized children with glomerular diseases participated in the study. For each patient, an overnight (9:00 p.m.-7:00 a.m.) urine was collected, followed by a 24-h urine (classified into four distinct periods: morning 7:00 a.m.-12:00 p.m., afternoon 12:00 p.m.-4:00 p.m., evening 4:00 p.m.-9:00 p.m., and overnight 9:00 p.m.-7:00 a.m.). The concentrations of protein, albumin, N-acetyl-beta-D-glucosaminidase, and epidermal growth factor (EGF) were measured and normalized by three correction factors (creatinine, osmolality, or specific gravity, respectively). Additionally, the 2nd overnight urine sample was grouped into different aliquots according to centrifugation, additives, storage temperature, or delayed processing. Twenty (14 boys, 6 girls) children were enrolled, with an average age of 11.3 years. Among the three correction factors, creatinine-normalized biomarkers provided the best agreements among different periods over 24 h. There were significant diurnal variations during 24 h in the concentrations of urinary protein, albumin, N-acetyl-beta-D-glucosaminidase, and EGF (p = 0.001, p = 0.003, p = 0.003, and p = 0.003, respectively). Evening urine overestimated 24-h urinary protein and albumin, while overnight urine underestimated 24-h urinary albumin. Urinary EGF showed low variability within a day or between the 2 days (coefficients of variation 10.2% and 10.6%, respectively) and excellent agreements (intraclass correlation coefficients > 0.9) with 24-h urinary concentration. Furthermore, urinary EGF was not affected by centrifugation, additives, storage temperature, or delayed processing of urine samples (all p > 0.05).  Conclusion: Given the diurnal variations of urinary biomarkers, urine samples should be collected during the same time period in clinical practice if possible. The results also extend the evidence for urinary EGF as a relatively stable biomarker applied in the future clinical practice. What is Known: • Urinary biomarkers have been widely used or discussed in making diagnoses and therapy regimens and estimating the prognosis of pediatric glomerular diseases. It remains unclear whether their levels would be affected by the time of sample collection, processing methods, and storage conditions in hospitalized children with glomerular diseases. What is New: • The levels of both commonly used biomarkers and novel biomarkers exhibited diurnal variations in hospitalized children with glomerular diseases. • Our results extend the evidence for urinary EGF as a relatively stable biomarker applied in the future clinical practice.

The deacylase sirtuin 5 reduces malonylation in nonmitochondrial metabolic pathways in diabetic kidney disease.

Early diabetic kidney disease (DKD) is marked by dramatic metabolic reprogramming due to nutrient excess, mitochondrial dysfunction, and increased renal energy requirements from hyperfiltration. We hypothesized that changes in metabolism in DKD may be regulated by Sirtuin 5 (SIRT5), a deacylase that removes posttranslational modifications derived from acyl-coenzyme A and has been demonstrated to regulate numerous metabolic pathways. We found decreased malonylation in the kidney cortex (∼80% proximal tubules) of type 2 diabetic BKS db/db mice, associated with increased SIRT5 expression. We performed a proteomics analysis of malonylated peptides and found that proteins with significantly decreased malonylated lysines in the db/db cortex were enriched in nonmitochondrial metabolic pathways: glycolysis and peroxisomal fatty acid oxidation. To confirm relevance of these findings in human disease, we analyzed diabetic kidney transcriptomic data from a cohort of Southwestern American Indians, which revealed a tubulointerstitial-specific increase in Sirt5 expression. These data were further corroborated by immunofluorescence data of SIRT5 from nondiabetic and DKD cohorts. Furthermore, overexpression of SIRT5 in cultured human proximal tubules demonstrated increased aerobic glycolysis. Conversely, we observed reduced glycolysis with decreased SIRT5 expression. These findings suggest that SIRT5 may lead to differential nutrient partitioning and utilization in DKD. Taken together, our findings highlight a previously unrecognized role for SIRT5 in metabolic reprogramming in DKD.

SGLT2 inhibitors mitigate kidney tubular metabolic and mTORC1 perturbations in youth-onset type 2 diabetes.

The molecular mechanisms of sodium-glucose cotransporter-2 (SGLT2) inhibitors (SGLT2i) remain incompletely understood. Single-cell RNA sequencing and morphometric data were collected from research kidney biopsies donated by young persons with type 2 diabetes (T2D), aged 12 to 21 years, and healthy controls (HCs). Participants with T2D were obese and had higher estimated glomerular filtration rates and mesangial and glomerular volumes than HCs. Ten T2D participants had been prescribed SGLT2i (T2Di[+]) and 6 not (T2Di[-]). Transcriptional profiles showed SGLT2 expression exclusively in the proximal tubular (PT) cluster with highest expression in T2Di(-) patients. However, transcriptional alterations with SGLT2i treatment were seen across nephron segments, particularly in the distal nephron. SGLT2i treatment was associated with suppression of transcripts in the glycolysis, gluconeogenesis, and tricarboxylic acid cycle pathways in PT, but had the opposite effect in thick ascending limb. Transcripts in the energy-sensitive mTORC1-signaling pathway returned toward HC levels in all tubular segments in T2Di(+), consistent with a diabetes mouse model treated with SGLT2i. Decreased levels of phosphorylated S6 protein in proximal and distal tubules in T2Di(+) patients confirmed changes in mTORC1 pathway activity. We propose that SGLT2i treatment benefits the kidneys by mitigating diabetes-induced metabolic perturbations via suppression of mTORC1 signaling in kidney tubules.

Single-cell transcriptomics reveals a mechanosensitive injury signaling pathway in early diabetic nephropathy.

BACKGROUND: Diabetic nephropathy (DN) is the leading cause of end-stage renal disease, and histopathologic glomerular lesions are among the earliest structural alterations of DN. However, the signaling pathways that initiate these glomerular alterations are incompletely understood. METHODS: To delineate the cellular and molecular basis for DN initiation, we performed single-cell and bulk RNA sequencing of renal cells from type 2 diabetes mice (BTBR ob/ob) at the early stage of DN. RESULTS: Analysis of differentially expressed genes revealed glucose-independent responses in glomerular cell types. The gene regulatory network upstream of glomerular cell programs suggested the activation of mechanosensitive transcriptional pathway MRTF-SRF predominantly taking place in mesangial cells. Importantly, activation of MRTF-SRF transcriptional pathway was also identified in DN glomeruli in independent patient cohort datasets. Furthermore, ex vivo kidney perfusion suggested that the regulation of MRTF-SRF is a common mechanism in response to glomerular hyperfiltration. CONCLUSIONS: Overall, our study presents a comprehensive single-cell transcriptomic landscape of early DN, highlighting mechanosensitive signaling pathways as novel targets of diabetic glomerulopathy.

Precision nephrology identified tumor necrosis factor activation variability in minimal change disease and focal segmental glomerulosclerosis.

The diagnosis of nephrotic syndrome relies on clinical presentation and descriptive patterns of injury on kidney biopsies, but not specific to underlying pathobiology. Consequently, there are variable rates of progression and response to therapy within diagnoses. Here, an unbiased transcriptomic-driven approach was used to identify molecular pathways which are shared by subgroups of patients with either minimal change disease (MCD) or focal segmental glomerulosclerosis (FSGS). Kidney tissue transcriptomic profile-based clustering identified three patient subgroups with shared molecular signatures across independent, North American, European, and African cohorts. One subgroup had significantly greater disease progression (Hazard Ratio 5.2) which persisted after adjusting for diagnosis and clinical measures (Hazard Ratio 3.8). Inclusion in this subgroup was retained even when clustering was limited to those with less than 25% interstitial fibrosis. The molecular profile of this subgroup was largely consistent with tumor necrosis factor (TNF) pathway activation. Two TNF pathway urine markers were identified, tissue inhibitor of metalloproteinases-1 (TIMP-1) and monocyte chemoattractant protein-1 (MCP-1), that could be used to predict an individual's TNF pathway activation score. Kidney organoids and single-nucleus RNA-sequencing of participant kidney biopsies, validated TNF-dependent increases in pathway activation score, transcript and protein levels of TIMP-1 and MCP-1, in resident kidney cells. Thus, molecular profiling identified a subgroup of patients with either MCD or FSGS who shared kidney TNF pathway activation and poor outcomes. A clinical trial testing targeted therapies in patients selected using urinary markers of TNF pathway activation is ongoing.

Epigenome-wide meta-analysis identifies DNA methylation biomarkers associated with diabetic kidney disease.

Type 1 diabetes affects over nine million individuals globally, with approximately 40% developing diabetic kidney disease. Emerging evidence suggests that epigenetic alterations, such as DNA methylation, are involved in diabetic kidney disease. Here we assess differences in blood-derived genome-wide DNA methylation associated with diabetic kidney disease in 1304 carefully characterised individuals with type 1 diabetes and known renal status from two cohorts in the United Kingdom-Republic of Ireland and Finland. In the meta-analysis, we identify 32 differentially methylated CpGs in diabetic kidney disease in type 1 diabetes, 18 of which are located within genes differentially expressed in kidneys or correlated with pathological traits in diabetic kidney disease. We show that methylation at 21 of the 32 CpGs predict the development of kidney failure, extending the knowledge and potentially identifying individuals at greater risk for diabetic kidney disease in type 1 diabetes.

Multi-Scalar Data Integration Links Glomerular Angiopoietin-Tie Signaling Pathway Activation With Progression of Diabetic Kidney Disease.

Diabetic kidney disease (DKD) is the leading cause of end-stage kidney disease (ESKD). Prognostic biomarkers reflective of underlying molecular mechanisms are critically needed for effective management of DKD. A three-marker panel was derived from a proteomics analysis of plasma samples by an unbiased machine learning approach from participants (N = 58) in the Clinical Phenotyping and Resource Biobank study. In combination with standard clinical parameters, this panel improved prediction of the composite outcome of ESKD or a 40% decline in glomerular filtration rate. The panel was validated in an independent group (N = 68), who also had kidney transcriptomic profiles. One marker, plasma angiopoietin 2 (ANGPT2), was significantly associated with outcomes in cohorts from the Cardiovascular Health Study (N = 3,183) and the Chinese Cohort Study of Chronic Kidney Disease (N = 210). Glomerular transcriptional angiopoietin/Tie (ANG-TIE) pathway scores, derived from the expression of 154 ANG-TIE signaling mediators, correlated positively with plasma ANGPT2 levels and kidney outcomes. Higher receptor expression in glomeruli and higher ANG-TIE pathway scores in endothelial cells corroborated potential functional effects in the kidney from elevated plasma ANGPT2 levels. Our work suggests that ANGPT2 is a promising prognostic endothelial biomarker with likely functional impact on glomerular pathogenesis in DKD.

Molecular programs associated with glomerular hyperfiltration in early diabetic kidney disease.

Hyperfiltration is a state of high glomerular filtration rate (GFR) observed in early diabetes that damages glomeruli, resulting in an iterative process of increasing filtration load on fewer and fewer remaining functional glomeruli. To delineate underlying cellular mechanisms of damage associated with hyperfiltration, transcriptional profiles of kidney biopsies from Pima Indians with type 2 diabetes with or without early-stage diabetic kidney disease were grouped into two hyperfiltration categories based on annual iothalamate GFR measurements. Twenty-six participants with a peak GFR measurement within two years of biopsy were categorized as the hyperfiltration group, and 26 in whom biopsy preceded peak GFR by over two years were considered pre-hyperfiltration. The hyperfiltration group had higher hemoglobin A1c, higher urine albumin-to-creatinine ratio, increased glomerular basement membrane width and lower podocyte density compared to the pre-hyperfiltration group. A glomerular 1240-gene transcriptional signature identified in the hyperfiltration group was enriched for endothelial stress response signaling genes, including endothelin-1, tec-kinase and transforming growth factor-ß1 pathways, with the majority of the transcripts mapped to endothelial and inflammatory cell clusters in kidney single cell transcriptional data. Thus, our analysis reveals molecular pathomechanisms associated with hyperfiltration in early diabetic kidney disease involving putative ligand-receptor pairs with downstream intracellular targets linked to cellular crosstalk between endothelial and mesangial cells.

Genome-wide meta-analysis and omics integration identifies novel genes associated with diabetic kidney disease.

AIMS/HYPOTHESIS: Diabetic kidney disease (DKD) is the leading cause of kidney failure and has a substantial genetic component. Our aim was to identify novel genetic factors and genes contributing to DKD by performing meta-analysis of previous genome-wide association studies (GWAS) on DKD and by integrating the results with renal transcriptomics datasets. METHODS: We performed GWAS meta-analyses using ten phenotypic definitions of DKD, including nearly 27,000 individuals with diabetes. Meta-analysis results were integrated with estimated quantitative trait locus data from human glomerular (N=119) and tubular (N=121) samples to perform transcriptome-wide association study. We also performed gene aggregate tests to jointly test all available common genetic markers within a gene, and combined the results with various kidney omics datasets. RESULTS: The meta-analysis identified a novel intronic variant (rs72831309) in the TENM2 gene associated with a lower risk of the combined chronic kidney disease (eGFR<60 ml/min per 1.73 m2) and DKD (microalbuminuria or worse) phenotype (p=9.8×10-9; although not withstanding correction for multiple testing, p>9.3×10-9). Gene-level analysis identified ten genes associated with DKD (COL20A1, DCLK1, EIF4E, PTPRN-RESP18, GPR158, INIP-SNX30, LSM14A and MFF; p<2.7×10-6). Integration of GWAS with human glomerular and tubular expression data demonstrated higher tubular AKIRIN2 gene expression in individuals with vs without DKD (p=1.1×10-6). The lead SNPs within six loci significantly altered DNA methylation of a nearby CpG site in kidneys (p<1.5×10-11). Expression of lead genes in kidney tubules or glomeruli correlated with relevant pathological phenotypes (e.g. TENM2 expression correlated positively with eGFR [p=1.6×10-8] and negatively with tubulointerstitial fibrosis [p=2.0×10-9], tubular DCLK1 expression correlated positively with fibrosis [p=7.4×10-16], and SNX30 expression correlated positively with eGFR [p=5.8×10-14] and negatively with fibrosis [p<2.0×10-16]). CONCLUSIONS/INTERPRETATION: Altogether, the results point to novel genes contributing to the pathogenesis of DKD. DATA AVAILABILITY: The GWAS meta-analysis results can be accessed via the type 1 and type 2 diabetes (T1D and T2D, respectively) and Common Metabolic Diseases (CMD) Knowledge Portals, and downloaded on their respective download pages ( https://t1d.hugeamp.org/downloads.html ; https://t2d.hugeamp.org/downloads.html ; https://hugeamp.org/downloads.html ).

Urinary Proteomics Identifies Cathepsin D as a Biomarker of Rapid eGFR Decline in Type 1 Diabetes.

OBJECTIVE: Understanding mechanisms underlying rapid estimated glomerular filtration rate (eGFR) decline is important to predict and treat kidney disease in type 1 diabetes (T1D). RESEARCH DESIGN AND METHODS: We performed a case-control study nested within four T1D cohorts to identify urinary proteins associated with rapid eGFR decline. Case and control subjects were categorized based on eGFR decline ≥3 and <1 mL/min/1.73 m2/year, respectively. We used targeted liquid chromatography-tandem mass spectrometry to measure 38 peptides from 20 proteins implicated in diabetic kidney disease. Significant proteins were investigated in complementary human cohorts and in mouse proximal tubular epithelial cell cultures. RESULTS: The cohort study included 1,270 participants followed a median 8 years. In the discovery set, only cathepsin D peptide and protein were significant on full adjustment for clinical and laboratory variables. In the validation set, associations of cathepsin D with eGFR decline were replicated in minimally adjusted models but lost significance with adjustment for albuminuria. In a meta-analysis with combination of discovery and validation sets, the odds ratio for the association of cathepsin D with rapid eGFR decline was 1.29 per SD (95% CI 1.07-1.55). In complementary human cohorts, urine cathepsin D was associated with tubulointerstitial injury and tubulointerstitial cathepsin D expression was associated with increased cortical interstitial fractional volume. In mouse proximal tubular epithelial cell cultures, advanced glycation end product-BSA increased cathepsin D activity and inflammatory and tubular injury markers, which were further increased with cathepsin D siRNA. CONCLUSIONS: Urine cathepsin D is associated with rapid eGFR decline in T1D and reflects kidney tubulointerstitial injury.

Molecular Characterization of Membranous Nephropathy.

BACKGROUND: Molecular characterization of nephropathies may facilitate pathophysiologic insight, development of targeted therapeutics, and transcriptome-based disease classification. Although membranous nephropathy (MN) is a common cause of adult-onset nephrotic syndrome, the molecular pathways of kidney damage in MN require further definition. METHODS: We applied a machine-learning framework to predict diagnosis on the basis of gene expression from the microdissected kidney tissue of participants in the Nephrotic Syndrome Study Network (NEPTUNE) cohort. We sought to identify differentially expressed genes between participants with MN versus those of other glomerulonephropathies across the NEPTUNE and European Renal cDNA Bank (ERCB) cohorts, to find MN-specific gene modules in a kidney-specific functional network, and to identify cell-type specificity of MN-specific genes using single-cell sequencing data from reference nephrectomy tissue. RESULTS: Glomerular gene expression alone accurately separated participants with MN from those with other nephrotic syndrome etiologies. The top predictive classifier genes from NEPTUNE participants were also differentially expressed in the ERCB participants with MN. We identified a signature of 158 genes that are significantly differentially expressed in MN across both cohorts, finding 120 of these in a validation cohort. This signature is enriched in targets of transcription factor NF-κB. Clustering these MN-specific genes in a kidney-specific functional network uncovered modules with functional enrichments, including in ion transport, cell projection morphogenesis, regulation of adhesion, and wounding response. Expression data from reference nephrectomy tissue indicated 43% of these genes are most highly expressed by podocytes. CONCLUSIONS: These results suggest that, relative to other glomerulonephropathies, MN has a distinctive molecular signature that includes upregulation of many podocyte-expressed genes, provides a molecular snapshot of MN, and facilitates insight into MN's underlying pathophysiology.

Serum Level of Polyubiquitinated PTEN and Loss of Kidney Function in American Indians With Type 2 Diabetes.

RATIONALE & OBJECTIVE: Fibrosis is a major driver of chronic kidney disease, and epithelial-mesenchymal transition (EMT) may contribute to its development. A polyubiquitinated form of phosphatase and tensin homolog (PTENK27polyUb) promotes EMT in vitro. Thus, it is a potentially useful biomarker of progressive kidney fibrosis and may predict loss of kidney function. STUDY DESIGN: Observational cohort study. SETTING & PARTICIPANTS: Southwest United States, American Indians (154 women, 80 men) with or at high risk for diabetic kidney disease (DKD). PREDICTORS: Serum level of PTENK27polyUb. OUTCOME: ≥40% loss of glomerular filtration rate (GFR) or onset of kidney failure. Kidney structural measures in a subset of study participants who underwent research kidney biopsies (n = 77). ANALYTICAL APPROACH: Cox proportional hazards models adjusted for age, sex, diabetes duration, hemoglobin A1c (HbA1c), blood pressure, use of renin angiotensin system (RAS) blockers, measured GFR, and albuminuria. Spearman correlations for associations with structural measures. RESULTS: At baseline, the participants' mean age was 42.8 ± 10.5 (SD) years, diabetes duration 11.5 ± 7.1 years, mean arterial pressure 90.5 ± 9.5 mm Hg, HbA1c 9.3 ± 2.4%, GFR 152 ± 45 mL/min, and median urinary albumin-creatinine ratio 38 (interquartile range, 14-215) mg/g. RAS blockers were being used by 64 participants (27.4%). A higher PTENK27polyUb value was associated with a greater risk of ≥40% loss of GFR during a median follow-up period of 6.3 years (HR for quartile 4 [Q4] vs Q1, 3.95 [95% CI, 2.23-6.98], P < 0.001). Serum PTENK27polyUb was associated with an increased risk of kidney failure over a median follow-up period of 15.8 years (HR for Q4 vs Q1, 5.66 [95% CI, 1.99-16.13], P = 0.001). Baseline serum PTENK27polyUb in the biopsy subset correlated with structural measures including glomerular basement membrane width (ρ = 0.370, P < 0.001) and mesangial fractional volume (ρ = 0.392, P < 0.001). LIMITATIONS: Small study in single population. CONCLUSIONS: Higher serum PTENK27polyUb is associated with increased risk for GFR decline and kidney failure in American Indians with type 2 diabetes.

Glomerular endothelial cell-podocyte stresses and crosstalk in structurally normal kidney transplants.

Increased podocyte detachment begins immediately after kidney transplantation and is associated with long-term allograft failure. We hypothesized that cell-specific transcriptional changes in podocytes and glomerular endothelial cells after transplantation would offer mechanistic insights into the podocyte detachment process. To test this, we evaluated cell-specific transcriptional profiles of glomerular endothelial cells and podocytes from 14 patients of their first-year surveillance biopsies with normal histology from low immune risk recipients with no post-transplant complications and compared these to biopsies of 20 healthy living donor controls. Glomerular endothelial cells from these surveillance biopsies were enriched for genes related to fluid shear stress, angiogenesis, and interferon signaling. In podocytes, pathways were enriched for genes in response to growth factor signaling and actin cytoskeletal reorganization but also showed evidence of podocyte stress as indicated by reduced nephrin (adhesion protein) gene expression. In parallel, transcripts coding for proteins required to maintain podocyte adherence to the underlying glomerular basement membrane were downregulated, including the major glomerular podocyte integrin α3 and the actin cytoskeleton-related gene synaptopodin. The reduction in integrin α3 protein expression in surveillance biopsies was confirmed by immunoperoxidase staining. The combined growth and stress response of patient allografts post-transplantation paralleled similar changes in a rodent model of nephrectomy-induced glomerular hypertrophic stress that progress to develop proteinuria and glomerulosclerosis with shortened kidney life span. Thus, even among patients with apparently healthy allografts with no detectable histologic abnormality including alloimmune injury, transcriptomic changes reflecting cell stresses are already set in motion that could drive hypertrophy-associated glomerular disease progression.

Assessment of Fibrinogen-like 2 (FGL2) in Human Chronic Kidney Disease through Transcriptomics Data Analysis.

Fibrinogen-like 2 (FGL2) was recently found to be associated with fibrosis in a mouse model of kidney damage and was proposed as a potential therapeutic target in chronic kidney disease (CKD). We assessed the association of renal FGL2 mRNA expression with the disease outcome in two independent CKD cohorts (NEPTUNE and Innsbruck CKD cohort) using Kaplan Meier survival analysis. The regulation of FGL2 in kidney biopsies of CKD patients as compared to healthy controls was further assessed in 13 human CKD transcriptomics datasets. The FGL2 protein expression in human renal tissue sections was determined via immunohistochemistry. The regulators of FGL2 mRNA expression in renal tissue were identified in the co-expression and upstream regulator analysis of FGL2-positive renal cells via the use of single-cell RNA sequencing data from the kidney precision medicine project (KPMP). Higher renal FGL2 mRNA expression was positively associated with kidney fibrosis and negatively associated with eGFR. Renal FGL2 mRNA expression was upregulated in CKD as compared with healthy controls and associated with CKD progression in the Innsbruck CKD cohort (p-value = 0.0036) and NEPTUNE cohort (p-value = 0.0048). The highest abundance of FGL2 protein in renal tissue was detected in the thick ascending limb of the loop of Henle and macula densa, proximal tubular cells, as well as in glomerular endothelial cells. The upstream regulator analysis identified TNF, IL1B, IFNG, NFKB1, and SP1 as factors potentially inducing FGL2-co-expressed genes, whereas factors counterbalancing FGL2-co-expressed genes included GLI1, HNF1B, or PPARGC1A. In conclusion, renal FGL2 mRNA expression is elevated in human CKD, and higher FGL2 levels are associated with fibrosis and worse outcomes.

Renin-angiotensin system inhibition reverses the altered triacylglycerol metabolic network in diabetic kidney disease.

OBJECTIVE: Dyslipidemia is a significant risk factor for progression of diabetic kidney disease (DKD). Determining the changes in individual lipids and lipid networks across a spectrum of DKD severity may identify lipids that are pathogenic to DKD progression. METHODS: We performed untargeted lipidomic analysis of kidney cortex tissue from diabetic db/db and db/db eNOS-/- mice along with non-diabetic littermate controls. A subset of mice were treated with the renin-angiotensin system (RAS) inhibitors, lisinopril and losartan, which improves the DKD phenotype in the db/db eNOS-/- mouse model. RESULTS: Of the three independent variables in this study, diabetes had the largest impact on overall lipid levels in the kidney cortex, while eNOS expression and RAS inhibition had smaller impacts on kidney lipid levels. Kidney lipid network architecture, particularly of networks involving glycerolipids such as triacylglycerols, was substantially disrupted by worsening kidney disease in the db/db eNOS-/- mice compared to the db/db mice, a feature that was reversed with RAS inhibition. This was associated with decreased expression of the stearoyl-CoA desaturases, Scd1 and Scd2, with RAS inhibition. CONCLUSIONS: In addition to the known salutary effect of RAS inhibition on DKD progression, our results suggest a previously unrecognized role for RAS inhibition on the kidney triacylglycerol lipid metabolic network.