Single-cell analysis identifies the interaction of altered renal tubules with basophils orchestrating kidney fibrosis.

Inflammation is an important component of fibrosis but immune processes that orchestrate kidney fibrosis are not well understood. Here we apply single-cell sequencing to a mouse model of kidney fibrosis. We identify a subset of kidney tubule cells with a profibrotic-inflammatory phenotype characterized by the expression of cytokines and chemokines associated with immune cell recruitment. Receptor-ligand interaction analysis and experimental validation indicate that CXCL1 secreted by profibrotic tubules recruits CXCR2+ basophils. In mice, these basophils are an important source of interleukin-6 and recruitment of the TH17 subset of helper T cells. Genetic deletion or antibody-based depletion of basophils results in reduced renal fibrosis. Human kidney single-cell, bulk gene expression and immunostaining validate a function for basophils in patients with kidney fibrosis. Collectively, these studies identify basophils as contributors to the development of renal fibrosis and suggest that targeting these cells might be a useful clinical strategy to manage chronic kidney disease.

The key role of altered tubule cell lipid metabolism in kidney disease development.

Kidney epithelial cells have very high energy requirements, which are largely met by fatty acid oxidation. Complex changes in lipid metabolism are observed in patients with kidney disease. Defects in fatty acid oxidation and increased lipid uptake, especially in the context of hyperlipidemia and proteinuria, contribute to this excess lipid build-up and exacerbate kidney disease development. Recent studies have also highlighted the role of increased de novo lipogenesis in kidney fibrosis. The defect in fatty acid oxidation causes energy starvation. Increased lipid uptake, synthesis, and lower fatty acid oxidation can cause toxic lipid build-up, reactive oxygen species generation, and mitochondrial damage. A better understanding of these metabolic processes may open new treatment avenues for kidney diseases by targeting lipid metabolism.

Unbiased Human Kidney Tissue Proteomics Identifies Matrix Metalloproteinase 7 as a Kidney Disease Biomarker.

SIGNIFICANCE STATEMENT: Although gene expression changes have been characterized in human diabetic kidney disease (DKD), unbiased tissue proteomics information for this condition is lacking. The authors conducted an unbiased aptamer-based proteomic analysis of samples from patients with DKD and healthy controls, identifying proteins with levels that associate with kidney function (eGFR) or fibrosis, after adjusting for key covariates. Overall, tissue gene expression only modestly correlated with tissue protein levels. Kidney protein and RNA levels of matrix metalloproteinase 7 (MMP7) strongly correlated with fibrosis and with eGFR. Single-cell RNA sequencing indicated that kidney tubule cells are an important source of MMP7. Furthermore, plasma MMP7 levels predicted future kidney function decline. These findings identify kidney tissue MMP7 as a biomarker of fibrosis and blood MMP7 as a biomarker for future kidney function decline. BACKGROUND: Diabetic kidney disease (DKD) is responsible for close to half of all ESKD cases. Although unbiased gene expression changes have been extensively characterized in human kidney tissue samples, unbiased protein-level information is not available. METHODS: We collected human kidney samples from 23 individuals with DKD and ten healthy controls, gathered associated clinical and demographics information, and implemented histologic analysis. We performed unbiased proteomics using the SomaScan platform and quantified the level of 1305 proteins and analyzed gene expression levels by bulk RNA and single-cell RNA sequencing (scRNA-seq). We validated protein levels in a separate cohort of kidney tissue samples as well as in 11,030 blood samples. RESULTS: Globally, human kidney transcript and protein levels showed only modest correlation. Our analysis identified 14 proteins with kidney tissue levels that correlated with eGFR and found that the levels of 152 proteins correlated with interstitial fibrosis. Of the identified proteins, matrix metalloprotease 7 (MMP7) showed the strongest association with both fibrosis and eGFR. The correlation between tissue MMP7 protein expression and kidney function was validated in external datasets. The levels of MMP7 RNA correlated with fibrosis in the primary and validation datasets. Findings from scRNA-seq pointed to proximal tubules, connecting tubules, and principal cells as likely cellular sources of increased tissue MMP7 expression. Furthermore, plasma MMP7 levels correlated not only with kidney function but also associated with prospective kidney function decline. CONCLUSIONS: Our findings, which underscore the value of human kidney tissue proteomics analysis, identify kidney tissue MMP7 as a diagnostic marker of kidney fibrosis and blood MMP7 as a biomarker for future kidney function decline.

Comprehensive Dipeptide Profiling and Quantitation by Capillary Electrophoresis and Liquid Chromatography Coupled with Tandem Mass Spectrometry.

Dipeptides have attracted much attention as post-amino acids with physical properties and functions different from those of amino acids. However, a given dipeptide cannot be distinguished by mass spectrometry from its structural isomer with an opposite amino acid binding order unless these isomers are separated before introduction, which complicates the comprehensive analysis of dipeptides. Herein, a novel analytical platform for dipeptide analysis by capillary electrophoresis tandem mass spectrometry and liquid chromatography tandem mass spectrometry is developed. This method is used to quantitate 335 dipeptides and achieves excellent separation of structural isomers with opposite binding orders, high correlation coefficients, and low instrumental detection limits (0.088-83.1 nM). Moreover, acceptable recoveries (70-135%) are observed for most tested dipeptides in chicken liver samples spiked both before and after preparation. The developed method is also applied to the quantitation of dipeptides in the livers of mice fed different diets to detect 236 dipeptides, and the shift from a normal diet to a high-fat diet is shown to increase/decrease (p < 0.05, fold-change < 0.5) the contents of 0/29 dipeptides, respectively. The developed method is expected to facilitate the search for new dipeptide applications such as novel functional components of foods and biomarkers of diseases.

Aging-associated renal disease in mice is fructokinase dependent.

Aging-associated kidney disease is usually considered a degenerative process associated with aging. Recently, it has been shown that animals can produce fructose endogenously, and that this can be a mechanism for causing kidney damage in diabetic nephropathy and in association with recurrent dehydration. We therefore hypothesized that low-level metabolism of endogenous fructose might play a role in aging-associated kidney disease. Wild-type and fructokinase knockout mice were fed a normal diet for 2 yr that had minimal (<5%) fructose content. At the end of 2 yr, wild-type mice showed elevations in systolic blood pressure, mild albuminuria, and glomerular changes with mesangial matrix expansion, variable mesangiolysis, and segmental thrombi. The renal injury was amplified by provision of high-salt diet for 3 wk, as noted by the presence of glomerular hypertrophy, mesangial matrix expansion, and alpha smooth muscle actin expression, and with segmental thrombi. Fructokinase knockout mice were protected from renal injury both at baseline and after high salt intake (3 wk) compared with wild-type mice. This was associated with higher levels of active (phosphorylated serine 1177) endothelial nitric oxide synthase in their kidneys. These studies suggest that aging-associated renal disease might be due to activation of specific metabolic pathways that could theoretically be targeted therapeutically, and raise the hypothesis that aging-associated renal injury may represent a disease process as opposed to normal age-related degeneration.

Single-cell transcriptomics and chromatin accessibility profiling elucidate the kidney-protective mechanism of mineralocorticoid receptor antagonists.

Mineralocorticoid excess commonly leads to hypertension (HTN) and kidney disease. In our study, we used single-cell expression and chromatin accessibility tools to characterize the mineralocorticoid target genes and cell types. We demonstrated that mineralocorticoid effects were established through open chromatin and target gene expression, primarily in principal and connecting tubule cells and, to a lesser extent, in segments of the distal convoluted tubule cells. We examined the kidney-protective effects of steroidal and nonsteroidal mineralocorticoid antagonists (MRAs), as well as of amiloride, an epithelial sodium channel inhibitor, in a rat model of deoxycorticosterone acetate, unilateral nephrectomy, and high-salt consumption-induced HTN and cardiorenal damage. All antihypertensive therapies protected against cardiorenal damage. However, finerenone was particularly effective in reducing albuminuria and improving gene expression changes in podocytes and proximal tubule cells, even with an equivalent reduction in blood pressure. We noted a strong correlation between the accumulation of injured/profibrotic tubule cells expressing secreted posphoprotein 1 (Spp1), Il34, and platelet-derived growth factor subunit b (Pdgfb) and the degree of fibrosis in rat kidneys. This gene signature also showed a potential for classifying human kidney samples. Our multiomics approach provides fresh insights into the possible mechanisms underlying HTN-associated kidney disease, the target cell types, the protective effects of steroidal and nonsteroidal MRAs, and amiloride.

Molecular pathways that drive diabetic kidney disease.

Kidney disease is a major driver of mortality among patients with diabetes and diabetic kidney disease (DKD) is responsible for close to half of all chronic kidney disease cases. DKD usually develops in a genetically susceptible individual as a result of poor metabolic (glycemic) control. Molecular and genetic studies indicate the key role of podocytes and endothelial cells in driving albuminuria and early kidney disease in diabetes. Proximal tubule changes show a strong association with the glomerular filtration rate. Hyperglycemia represents a key cellular stress in the kidney by altering cellular metabolism in endothelial cells and podocytes and by imposing an excess workload requiring energy and oxygen for proximal tubule cells. Changes in metabolism induce early adaptive cellular hypertrophy and reorganization of the actin cytoskeleton. Later, mitochondrial defects contribute to increased oxidative stress and activation of inflammatory pathways, causing progressive kidney function decline and fibrosis. Blockade of the renin-angiotensin system or the sodium-glucose cotransporter is associated with cellular protection and slowing kidney function decline. Newly identified molecular pathways could provide the basis for the development of much-needed novel therapeutics.

NAD+ precursor supplementation prevents mtRNA/RIG-I-dependent inflammation during kidney injury.

Our understanding of how global changes in cellular metabolism contribute to human kidney disease remains incompletely understood. Here we show that nicotinamide adenine dinucleotide (NAD+) deficiency drives mitochondrial dysfunction causing inflammation and kidney disease development. Using unbiased global metabolomics in healthy and diseased human kidneys, we identify NAD+ deficiency as a disease signature. Furthermore using models of cisplatin- or ischaemia-reperfusion induced kidney injury in male mice we observed NAD+ depletion Supplemental nicotinamide riboside or nicotinamide mononucleotide restores NAD+ levels and improved kidney function. We find that cisplatin exposure causes cytosolic leakage of mitochondrial RNA (mtRNA) and activation of the cytosolic pattern recognition receptor retinoic acid-inducible gene I (RIG-I), both of which can be ameliorated by restoring NAD+. Male mice with RIG-I knock-out (KO) are protected from cisplatin-induced kidney disease. In summary, we demonstrate that the cytosolic release of mtRNA and RIG-I activation is an NAD+-sensitive mechanism contributing to kidney disease.

Long non-coding RNA lnc-CHAF1B-3 promotes renal interstitial fibrosis by regulating EMT-related genes in renal proximal tubular cells.

Renal interstitial fibrosis (RIF) is a common pathological manifestation of chronic kidney diseases. Epithelial-mesenchymal transition (EMT) of tubular epithelial cells is considered a major cause of RIF. Although long non-coding RNAs (lncRNAs) are reportedly involved in various pathophysiological processes, the roles and underlying molecular mechanisms of lncRNAs in the progression of RIF are poorly understood. In this study, we investigated the function of lncRNAs in RIF. Microarray assays showed that expression of the lncRNA lnc-CHAF1B-3 (also called claudin 14 antisense RNA 1) was significantly upregulated in human renal proximal tubular cells by both transforming growth factor-ß1 (TGF-ß1) and hypoxic stimulation, accompanied with increased expression of EMT-related genes. Knockdown of lnc-CHAF1B-3 significantly suppressed TGF-ß1-induced upregulated expression of collagen type I alpha 1, cadherin-2, plasminogen activator inhibitor-1, snail family transcriptional repressor I (SNAI1) and SNAI2. Quantitative reverse transcriptase PCR analyses of paraffin-embedded kidney biopsy samples from IgA nephropathy patients revealed lnc-CHAF1B-3 expression was correlated positively with urinary protein levels and correlated negatively with estimated glomerular filtration rate. In situ hybridization demonstrated that lnc-CHAF1B-3 is expressed only in proximal tubules. These findings suggest lnc-CHAF1B-3 affects the progression of RIF by regulating EMT-related signaling. Thus, lnc-CHAF1B-3 is a potential target in the treatment of RIF.

ACSS2 gene variants determine kidney disease risk by controlling de novo lipogenesis in kidney tubules.

Worldwide, over 800 million people are affected by kidney disease, yet its pathogenesis remains elusive, hindering the development of novel therapeutics. In this study, we used kidney-specific expression of quantitative traits and single-nucleus open chromatin analysis to show that genetic variants linked to kidney dysfunction on chromosome 20 target the acyl-CoA synthetase short-chain family 2 (ACSS2). By generating ACSS2-KO mice, we demonstrated their protection from kidney fibrosis in multiple disease models. Our analysis of primary tubular cells revealed that ACSS2 regulated de novo lipogenesis (DNL), causing NADPH depletion and increasing ROS levels, ultimately leading to NLRP3-dependent pyroptosis. Additionally, we discovered that pharmacological inhibition or genetic ablation of fatty acid synthase safeguarded kidney cells against profibrotic gene expression and prevented kidney disease in mice. Lipid accumulation and the expression of genes related to DNL were elevated in the kidneys of patients with fibrosis. Our findings pinpoint ACSS2 as a critical kidney disease gene and reveal the role of DNL in kidney disease.

The multifaceted role of kidney tubule mitochondrial dysfunction in kidney disease development.

More than 800 million people suffer from kidney disease. Genetic studies and follow-up animal models and cell biological experiments indicate the key role of proximal tubule metabolism. Kidneys have one of the highest mitochondrial densities. Mitochondrial biogenesis, mitochondrial fusion and fission, and mitochondrial recycling, such as mitophagy are critical for proper mitochondrial function. Mitochondrial dysfunction can lead to an energetic crisis, orchestrate different types of cell death (apoptosis, necroptosis, pyroptosis, and ferroptosis), and influence cellular calcium levels and redox status. Collectively, mitochondrial defects in renal tubules contribute to epithelial atrophy, inflammation, or cell death, orchestrating kidney disease development.

Single-cell analysis highlights differences in druggable pathways underlying adaptive or fibrotic kidney regeneration.

The kidney has tremendous capacity to repair after acute injury, however, pathways guiding adaptive and fibrotic repair are poorly understood. We developed a model of adaptive and fibrotic kidney regeneration by titrating ischemic injury dose. We performed detailed biochemical and histological analysis and profiled transcriptomic changes at bulk and single-cell level (> 110,000 cells) over time. Our analysis highlights kidney proximal tubule cells as key susceptible cells to injury. Adaptive proximal tubule repair correlated with fatty acid oxidation and oxidative phosphorylation. We identify a specific maladaptive/profibrotic proximal tubule cluster after long ischemia, which expresses proinflammatory and profibrotic cytokines and myeloid cell chemotactic factors. Druggability analysis highlights pyroptosis/ferroptosis as vulnerable pathways in these profibrotic cells. Pharmacological targeting of pyroptosis/ferroptosis in vivo pushed cells towards adaptive repair and ameliorates fibrosis. In summary, our single-cell analysis defines key differences in adaptive and fibrotic repair and identifies druggable pathways for pharmacological intervention to prevent kidney fibrosis.

Epigenomic and transcriptomic analyses define core cell types, genes and targetable mechanisms for kidney disease.

More than 800 million people suffer from kidney disease, yet the mechanism of kidney dysfunction is poorly understood. In the present study, we define the genetic association with kidney function in 1.5 million individuals and identify 878 (126 new) loci. We map the genotype effect on the methylome in 443 kidneys, transcriptome in 686 samples and single-cell open chromatin in 57,229 kidney cells. Heritability analysis reveals that methylation variation explains a larger fraction of heritability than gene expression. We present a multi-stage prioritization strategy and prioritize target genes for 87% of kidney function loci. We highlight key roles of proximal tubules and metabolism in kidney function regulation. Furthermore, the causal role of SLC47A1 in kidney disease is defined in mice with genetic loss of Slc47a1 and in human individuals carrying loss-of-function variants. Our findings emphasize the key role of bulk and single-cell epigenomic information in translating genome-wide association studies into identifying causal genes, cellular origins and mechanisms of complex traits.

Genome-wide association studies identify the role of caspase-9 in kidney disease.

Genome-wide association studies (GWAS) have identified hundreds of genetic risk regions for kidney dysfunction [estimated glomerular filtration rate (eGFR)]; however, the causal genes, cell types, and pathways are poorly understood. Integration of GWAS and human kidney expression of quantitative trait analysis using Bayesian colocations, transcriptome-wide association studies, and summary-based Mendelian randomization studies prioritized caspase-9 (CASP9) as a kidney disease risk gene. Human kidney single-cell epigenetic and immunostaining studies indicated kidney tubule cells as a disease-causing cell type. Mice with genetic deletion or pharmacological inhibition of CASP9 showed lower apoptosis while having improved mitophagy, resulting in dampened activation of cytosolic nucleotide sensing pathways (cGAS-STING), reduction of inflammation, and protection from acute kidney disease or renal fibrosis. In summary, here, we prioritized CASP9 as an eGFR GWAS target gene and demonstrated the causal role of CASP9 in kidney disease development via improving mitophagy and lowering inflammation and apoptosis.

Transcriptome-wide association analysis identifies DACH1 as a kidney disease risk gene that contributes to fibrosis.

Genome-wide association studies (GWAS) for kidney function identified hundreds of risk regions; however, the causal variants, target genes, cell types, and disease mechanisms remain poorly understood. Here, we performed transcriptome-wide association studies (TWAS), summary Mendelian randomization, and MetaXcan to identify genes whose expression mediates the genotype effect on the phenotype. Our analyses identified Dachshund homolog 1 (DACH1), a cell-fate determination factor. GWAS risk variant was associated with lower DACH1 expression in human kidney tubules. Human and mouse kidney single-cell open chromatin data (snATAC-Seq) prioritized estimated glomerular filtration rate (eGFR) GWAS variants located on an intronic regulatory region in distal convoluted tubule cells. CRISPR-Cas9-mediated gene editing confirmed the role of risk variants in regulating DACH1 expression. Mice with tubule-specific Dach1 deletion developed more severe renal fibrosis both in folic acid and diabetic kidney injury models. Mice with tubule-specific Dach1 overexpression were protected from folic acid nephropathy. Single-cell RNA sequencing, chromatin immunoprecipitation, and functional analysis indicated that DACH1 controls the expression of cell cycle and myeloid chemotactic factors, contributing to macrophage infiltration and fibrosis development. In summary, integration of GWAS, TWAS, single-cell epigenome, expression analyses, gene editing, and functional validation in different mouse kidney disease models identified DACH1 as a kidney disease risk gene.

DACH1 protects podocytes from experimental diabetic injury and modulates PTIP-H3K4Me3 activity.

Dachshund homolog 1 (DACH1), a key cell-fate determinant, regulates transcription by DNA sequence-specific binding. We identified diminished Dach1 expression in a large-scale screen for mutations that convert injury-resistant podocytes into injury-susceptible podocytes. In diabetic kidney disease (DKD) patients, podocyte DACH1 expression levels are diminished, a condition that strongly correlates with poor clinical outcomes. Global Dach1 KO mice manifest renal hypoplasia and die perinatally. Podocyte-specific Dach1 KO mice, however, maintain normal glomerular architecture at baseline, but rapidly exhibit podocyte injury after diabetes onset. Furthermore, podocyte-specific augmentation of DACH1 expression in mice protects from DKD. Combined RNA sequencing and in silico promoter analysis reveal conversely overlapping glomerular transcriptomic signatures between podocyte-specific Dach1 and Pax transactivation-domain interacting protein (Ptip) KO mice, with upregulated genes possessing higher-than-expected numbers of promoter Dach1-binding sites. PTIP, an essential component of the activating histone H3 lysine 4 trimethylation (H3K4Me3) complex, interacts with DACH1 and is recruited by DACH1 to its promoter-binding sites. DACH1-PTIP recruitment represses transcription and reduces promoter H3K4Me3 levels. DACH1 knockdown in podocytes combined with hyperglycemia triggers target gene upregulation and increases promoter H3K4Me3. These findings reveal that in DKD, diminished DACH1 expression enhances podocyte injury vulnerability via epigenetic derepression of its target genes.

NAD+ flux is maintained in aged mice despite lower tissue concentrations.

NAD+ is an essential coenzyme for all living cells. NAD+ concentrations decline with age, but whether this reflects impaired production or accelerated consumption remains unclear. We employed isotope tracing and mass spectrometry to probe age-related changes in NAD+ metabolism across tissues. In aged mice, we observed modest tissue NAD+ depletion (median decrease ∼30%). Circulating NAD+ precursors were not significantly changed, and isotope tracing showed the unimpaired synthesis of nicotinamide from tryptophan. In most tissues of aged mice, turnover of the smaller tissue NAD+ pool was modestly faster such that absolute NAD+ biosynthetic flux was maintained, consistent with more active NAD+-consuming enzymes. Calorie restriction partially mitigated age-associated NAD+ decline by decreasing consumption. Acute inflammatory stress induced by LPS decreased NAD+ by impairing synthesis in both young and aged mice. Thus, the decline in NAD+ with normal aging is relatively subtle and occurs despite maintained NAD+ production, likely due to increased consumption.

The Nuclear Receptor ESRRA Protects from Kidney Disease by Coupling Metabolism and Differentiation.

Kidney disease is poorly understood because of the organ's cellular diversity. We used single-cell RNA sequencing not only in resolving differences in injured kidney tissue cellular composition but also in cell-type-specific gene expression in mouse models of kidney disease. This analysis highlighted major changes in cellular diversity in kidney disease, which markedly impacted whole-kidney transcriptomics outputs. Cell-type-specific differential expression analysis identified proximal tubule (PT) cells as the key vulnerable cell type. Through unbiased cell trajectory analyses, we show that PT cell differentiation is altered in kidney disease. Metabolism (fatty acid oxidation and oxidative phosphorylation) in PT cells showed the strongest and most reproducible association with PT cell differentiation and disease. Coupling of cell differentiation and the metabolism was established by nuclear receptors (estrogen-related receptor alpha [ESRRA] and peroxisomal proliferation-activated receptor alpha [PPARA]) that directly control metabolic and PT-cell-specific gene expression in mice and patient samples while protecting from kidney disease in the mouse model.

Fructose increases the activity of sodium hydrogen exchanger in renal proximal tubules that is dependent on ketohexokinase.

High fructose intake has been known to induce metabolic syndrome in laboratory animals and humans. Although fructose intake enhances sodium reabsorption and elevates blood pressure, role of fructose metabolism in this process has not been studied. Here we show that by ketohexokinase - the primary enzyme of fructose - is involved in regulation of renal sodium reabsorption and blood pressure via activation of the sodium hydrogen exchanger in renal proximal tubular cells. First, wild-type and ketohexokinase knockout mice (Male, C57BL/6) were fed fructose water or tap water with or without a high salt diet. Only wild type mice fed the combination of fructose water and high salt diet displayed increased systolic blood pressure and decreased urinary sodium excretion. In contrast, ketohexokinase knockout mice were protected. Second, urinary sodium excretion after intraperitoneal saline administration was reduced with the decreased phosphorylation of sodium hydrogen exchanger 3 in fructose-fed WT; these changes were not observed in the ketohexokinase knockout mice, however. Third, knockdown of ketohexokinase attenuated fructose-mediated increases of NHE activity with decreased cAMP levels in porcine renal proximal tubular cells (LLC-PK1). In conclusion, fructose metabolism by ketohexokinase increases sodium hydrogen exchanger activity in renal proximal tubular cells via decreased intracellular cAMP level, resulting in increased renal sodium reabsorption and blood pressure in mice.