Epigenetic reprogramming driving successful and failed repair in acute kidney injury.

Acute kidney injury (AKI) causes epithelial damage followed by subsequent repair. While successful repair restores kidney function, this process is often incomplete and can lead to chronic kidney disease (CKD) in a process called failed repair. To better understand the epigenetic reprogramming driving this AKI-to-CKD transition, we generated a single-nucleus multiomic atlas for the full mouse AKI time course, consisting of ~280,000 single-nucleus transcriptomes and epigenomes. We reveal cell-specific dynamic alterations in gene regulatory landscapes reflecting, especially, activation of proinflammatory pathways. We further generated single-nucleus multiomic data from four human AKI samples including validation by genome-wide identification of nuclear factor κB binding sites. A regularized regression analysis identifies key regulators involved in both successful and failed repair cell fate, identifying the transcription factor CREB5 as a regulator of both successful and failed tubular repair that also drives proximal tubular cell proliferation after injury. Our interspecies multiomic approach provides a foundation to comprehensively understand cell states in AKI.

Multi-omics profiling of mouse polycystic kidney disease progression at a single cell resolution.

Autosomal dominant polycystic kidney disease (ADPKD) is the most common hereditary kidney disease and causes significant morbidity, ultimately leading to end-stage kidney disease. PKD pathogenesis is characterized by complex and dynamic alterations in multiple cell types during disease progression, hampering a deeper understanding of disease mechanism and the development of therapeutic approaches. Here, we generate a single nucleus multimodal atlas of an orthologous mouse PKD model at early, mid and late timepoints, consisting of 125,434 single-nucleus transcriptomic and epigenetic multiomes. We catalogue differentially expressed genes and activated epigenetic regions in each cell type during PKD progression, characterizing cell-type-specific responses to Pkd1 deletion. We describe heterogeneous, atypical collecting duct cells as well as proximal tubular cells that constitute cyst epithelia in PKD. The transcriptional regulation of the cyst lining cell marker GPRC5A is conserved between mouse and human PKD cystic epithelia, suggesting shared gene regulatory pathways. Our single nucleus multiomic analysis of mouse PKD provides a foundation to understand the earliest changes molecular deregulation in a mouse model of PKD at a single-cell resolution.

Injury-induced Foxm1 expression in the mouse kidney drives epithelial proliferation by a cyclin F-dependent mechanism.

Acute kidney injury (AKI) strongly upregulates the transcription factor Foxm1 in the proximal tubule in vivo, and Foxm1 drives epithelial proliferation in vitro. Here, we report that deletion of Foxm1 either with a nephron-specific Cre driver or by inducible global deletion reduced proximal tubule proliferation after ischemic injury in vivo. Foxm1 deletion led to increased AKI to chronic kidney disease transition, with enhanced fibrosis and ongoing tubule injury 6 weeks after injury. We report ERK mediated FOXM1 induction downstream of the EGFR in primary proximal tubule cells. We defined FOXM1 genomic binding sites by cleavage under targets and release using nuclease (CUT&RUN) and compared the genes located near FOXM1 binding sites with genes downregulated in primary proximal tubule cells after FOXM1 knockdown. The aligned data sets revealed the cell cycle regulator cyclin F (CCNF) as a putative FOXM1 target. We identified 2 cis regulatory elements that bound FOXM1 and regulated CCNF expression, demonstrating that Ccnf is strongly induced after kidney injury and that Foxm1 deletion abrogates Ccnf expression in vivo and in vitro. Knockdown of CCNF also reduced proximal tubule proliferation in vitro. These studies identify an ERK/FOXM1/CCNF signaling pathway that regulates injury-induced proximal tubule cell proliferation.

TNIK depletion induces inflammation and apoptosis in injured renal proximal tubule epithelial cells.

In the aftermath of acute kidney injury (AKI), surviving proximal tubule epithelia repopulate injured tubules to promote repair. However, a portion of cells fail to repair [termed failed-repair proximal tubule cells (FR-PTCs)] and exert ongoing proinflammatory and profibrotic effects. To better understand the molecular drivers of the FR-PTC state, we reanalyzed a mouse ischemia-reperfusion injury single-nucleus RNA-sequencing (snRNA-seq) atlas to identify Traf2 and Nck interacting kinase (Tnik) to be exclusively expressed in FR-PTCs but not in healthy or acutely injured proximal tubules after AKI (2 and 6 wk) in mice. We confirmed expression of Tnik protein in injured mouse and human tissues by immunofluorescence. Then, to determine the functional role of Tnik in FR-PTCs, we depleted TNIK with siRNA in two human renal proximal tubule epithelial cell lines (primary and immortalized hRPTECs) and analyzed each by bulk RNA-sequencing. Pathway analysis revealed significant upregulation of inflammatory signaling pathways, whereas pathways associated with differentiated proximal tubules such as organic acid transport were significantly downregulated. TNIK gene knockdown drove reduced cell viability and increased apoptosis, including differentially expressed poly(ADP-ribose) polymerase (PARP) family members, cleaved PARP-1 fragments, and increased annexin V binding to phosphatidylserine. Together, these results indicate that Tnik upregulation in FR-PTCs acts in a compensatory fashion to suppress inflammation and promote proximal tubule epithelial cell survival after injury. Modulating TNIK activity may represent a prorepair therapeutic strategy after AKI.NEW & NOTEWORTHY The molecular drivers of successful and failed repair in the proximal tubule after acute kidney injury (AKI) are incompletely understood. We identified Traf2 and Nck interacting kinase (Tnik) to be exclusively expressed in failed-repair proximal tubule cells after AKI. We tested the effect of siTNIK depletion in two proximal tubule cell lines followed by bulk RNA-sequencing analysis. Our results indicate that TNIK acts to suppress inflammatory signaling and apoptosis in injured renal proximal tubule epithelial cells to promote cell survival.

Epigenetic reprogramming driving successful and failed repair in acute kidney injury.

Acute kidney injury (AKI) causes epithelial damage followed by subsequent repair. While successful repair restores kidney function, this process is often incomplete and can lead to chronic kidney disease (CKD) in a process called failed repair. To better understand the epigenetic reprogramming driving this AKI-to-CKD transition we generated a single nucleus multiomic atlas for the full mouse AKI time course, consisting of ~280,000 single nucleus transcriptomes and epigenomes. We reveal cell-specific dynamic alterations in gene regulatory landscapes reflecting especially activation of proinflammatory pathways. We further generated single nucleus multiomic data from four human AKI samples including validation by genome-wide identification of NF-kB binding sites. A regularized regression analysis identifies key regulators involved in both successful and failed repair cell fate, identifying the transcription factor CREB5 as a regulator of both successful and failed tubular repair that also drives proximal tubule cell proliferation after injury. Our interspecies multiomic approach provides a foundation to comprehensively understand cell states in AKI.

Predicting proximal tubule failed repair drivers through regularized regression analysis of single cell multiomic sequencing.

Renal proximal tubule epithelial cells have considerable intrinsic repair capacity following injury. However, a fraction of injured proximal tubule cells fails to undergo normal repair and assumes a proinflammatory and profibrotic phenotype that may promote fibrosis and chronic kidney disease. The healthy to failed repair change is marked by cell state-specific transcriptomic and epigenomic changes. Single nucleus joint RNA- and ATAC-seq sequencing offers an opportunity to study the gene regulatory networks underpinning these changes in order to identify key regulatory drivers. We develop a regularized regression approach to construct genome-wide parametric gene regulatory networks using multiomic datasets. We generate a single nucleus multiomic dataset from seven adult human kidney samples and apply our method to study drivers of a failed injury response associated with kidney disease. We demonstrate that our approach is a highly effective tool for predicting key cis- and trans-regulatory elements underpinning the healthy to failed repair transition and use it to identify NFAT5 as a driver of the maladaptive proximal tubule state.

Leveraging multimodal chromatin profiling to identify a new potential driver of diabetic kidney disease.

The 3-dimensional nature of chromatin architecture plays crucial roles in regulating gene expression in development, homeostasis, and disease. Until recently, however, comprehensive chromatin profiling in human kidneys has been lacking. In this issue, Eun and Kim et al. employed a multimodal approach by integrating a single-nucleus assay for transposase-accessible chromatin sequencing, chromatin immunoprecipitation sequencing, and Hi-C (a method to comprehensively detect chromatin interactions) to investigate how the epigenetic landscape is altered in diabetic kidney disease.

Mosaic loss of Y chromosome is associated with aging and epithelial injury in chronic kidney disease.

BACKGROUND: Mosaic loss of Y chromosome (LOY) is the most common chromosomal alteration in aging men. Here, we use single-cell RNA and ATAC sequencing to show that LOY is present in the kidney and increases with age and chronic kidney disease. RESULTS: The likelihood of a cell having LOY varies depending on its location in the nephron. Cortical epithelial cell types have a greater proportion of LOY than medullary or glomerular cell types, which may reflect their proliferative history. Proximal tubule cells are the most abundant cell type in the cortex and are susceptible to hypoxic injury. A subset of these cells acquires a pro-inflammatory transcription and chromatin accessibility profile associated with expression of HAVCR1, VCAM1, and PROM1. These injured epithelial cells have the greatest proportion of LOY and their presence predicts future kidney function decline. Moreover, proximal tubule cells with LOY are more likely to harbor additional large chromosomal gains and express pro-survival pathways. Spatial transcriptomics localizes injured proximal tubule cells to a pro-fibrotic microenvironment where they adopt a secretory phenotype and likely communicate with infiltrating immune cells. CONCLUSIONS: We hypothesize that LOY is an indicator of increased DNA damage and potential marker of cellular senescence that can be applied to single-cell datasets in other tissues.

Elucidating the Proximal Tubule HNF4A Gene Regulatory Network in Human Kidney Organoids.

SIGNIFICANCE STATEMENT: HNF4 genes promote proximal tubule differentiation in mice, but their function in human nephrogenesis is not fully defined. This study uses human pluripotent stem cell (PSC)-derived kidney organoids as a model to investigate HNF4A and HNF4G functions. The loss of HNF4A , but not HNF4G , impaired reabsorption-related molecule expression and microvilli formation in human proximal tubules. Cleavage under targets and release using nuclease (CUT&RUN) sequencing and CRISPR-mediated transcriptional activation (CRISPRa) further confirm that HNF4A directly regulates its target genes. Human kidney organoids provide a good model for studying transcriptional regulation in human kidney development. BACKGROUND: The proximal tubule plays a major role in electrolyte homeostasis. Previous studies have shown that HNF4A regulates reabsorption-related genes and promotes proximal tubule differentiation during murine kidney development. However, the functions and gene regulatory mechanisms of HNF4 family genes in human nephrogenesis have not yet been investigated. METHODS: We generated HNF4A -knock out (KO), HNF4G -KO, and HNF4A/4G -double KO human pluripotent stem cell lines, differentiated each into kidney organoids, and used immunofluorescence analysis, electron microscopy, and RNA-seq to analyze them. We probed HNF4A-binding sites genome-wide by cleavage under targets and release using nuclease sequencing in both human adult kidneys and kidney organoid-derived proximal tubular cells. Clustered Regularly Interspaced Short Palindromic Repeats-mediated transcriptional activation validated HNF4A and HNF4G function in proximal tubules during kidney organoid differentiation. RESULTS: Organoids lacking HNF4A , but not HNF4G , showed reduced expression of transport-related, endocytosis-related, and brush border-related genes, as well as disorganized brush border structure in the apical lumen of the organoid proximal tubule. Cleavage under targets and release using nuclease revealed that HNF4A primarily bound promoters and enhancers of genes that were downregulated in HNF4A -KO, suggesting direct regulation. Induced expression of HNF4A or HNF4G by CRISPR-mediated transcriptional activation drove increased expression of selected target genes during kidney organoid differentiation. CONCLUSIONS: This study reveals regulatory mechanisms of HNF4A and HNF4G during human proximal tubule differentiation. The experimental strategy can be applied more broadly to investigate transcriptional regulation in human kidney development.

A single-cell multiomic analysis of kidney organoid differentiation.

Kidney organoids differentiated from pluripotent stem cells are powerful models of kidney development and disease but are characterized by cell immaturity and off-target cell fates. Comparing the cell-specific gene regulatory landscape during organoid differentiation with human adult kidney can serve to benchmark progress in differentiation at the epigenome and transcriptome level for individual organoid cell types. Using single-cell multiome and histone modification analysis, we report more broadly open chromatin in organoid cell types compared to the human adult kidney. We infer enhancer dynamics by cis-coaccessibility analysis and validate an enhancer driving transcription of HNF1B by CRISPR interference both in cultured proximal tubule cells and also during organoid differentiation. Our approach provides an experimental framework to judge the cell-specific maturation state of human kidney organoids and shows that kidney organoids can be used to validate individual gene regulatory networks that regulate differentiation.

Tubular Epithelial Cell HMGB1 Promotes AKI-CKD Transition by Sensitizing Cycling Tubular Cells to Oxidative Stress: A Rationale for Targeting HMGB1 during AKI Recovery.

SIGNIFICANCE STATEMENT: Cells undergoing necrosis release extracellular high mobility group box (HMGB)-1, which triggers sterile inflammation upon AKI in mice. Neither deletion of HMGB1 from tubular epithelial cells, nor HMGB1 antagonism with small molecules, affects initial ischemic tubular necrosis and immediate GFR loss upon unilateral ischemia/reperfusion injury (IRI). On the contrary, tubular cell-specific HMGB1 deficiency, and even late-onset pharmacological HMGB1 inhibition, increased functional and structural recovery from AKI, indicating that intracellular HMGB1 partially counters the effects of extracellular HMGB1. In vitro studies indicate that intracellular HMGB1 decreases resilience of tubular cells from prolonged ischemic stress, as in unilateral IRI. Intracellular HMGB1 is a potential target to enhance kidney regeneration and to improve long-term prognosis in AKI. BACKGROUND: Late diagnosis is a hurdle for treatment of AKI, but targeting AKI-CKD transition may improve outcomes. High mobility group box-1 (HMGB1) is a nuclear regulator of transcription and a driver of necroinflammation in AKI. We hypothesized that HMGB1 would also modulate AKI-CKD transition in other ways. METHODS: We conducted single-cell transcriptome analysis of human and mouse AKI and mouse in vivo and in vitro studies with tubular cell-specific depletion of Hmgb1 and HMGB1 antagonists. RESULTS: HMGB1 was ubiquitously expressed in kidney cells. Preemptive HMGB1 antagonism with glycyrrhizic acid (Gly) and ethyl pyruvate (EP) did not affect postischemic AKI but attenuated AKI-CKD transition in a model of persistent kidney hypoxia. Consistently, tubular Hmgb1 depletion in Pax8 rtTA, TetO Cre, Hmgb1fl/fl mice did not protect from AKI, but from AKI-CKD transition. In vitro studies confirmed that absence of HMGB1 or HMGB1 inhibition with Gly and EP does not affect ischemic necrosis of growth-arrested differentiated tubular cells but increased the resilience of cycling tubular cells that survived the acute injury to oxidative stress. This effect persisted when neutralizing extracellular HMGB1 with 2G7. Consistently, late-onset HMGB1 blockade with EP started after the peak of ischemic AKI in mice prevented AKI-CKD transition, even when 2G7 blocked extracellular HMGB1. CONCLUSION: Treatment of AKI could become feasible when ( 1 ) focusing on long-term outcomes of AKI; ( 2 ) targeting AKI-CKD transition with drugs initiated after the AKI peak; and ( 3 ) targeting with drugs that block HMGB1 in intracellular and extracellular compartments.

Quality assessment and refinement of chromatin accessibility data using a sequence-based predictive model.

Chromatin accessibility assays are central to the genome-wide identification of gene regulatory elements associated with transcriptional regulation. However, the data have highly variable quality arising from several biological and technical factors. To surmount this problem, we developed a sequence-based machine learning method to evaluate and refine chromatin accessibility data. Our framework, gapped k-mer SVM quality check (gkmQC), provides the quality metrics for a sample based on the prediction accuracy of the trained models. We tested 886 DNase-seq samples from the ENCODE/Roadmap projects to demonstrate that gkmQC can effectively identify "high-quality" (HQ) samples with low conventional quality scores owing to marginal read depths. Peaks identified in HQ samples are more accurately aligned at functional regulatory elements, show greater enrichment of regulatory elements harboring functional variants, and explain greater heritability of phenotypes from their relevant tissues. Moreover, gkmQC can optimize the peak-calling threshold to identify additional peaks, especially for rare cell types in single-cell chromatin accessibility data.

Defining cellular complexity in human autosomal dominant polycystic kidney disease by multimodal single cell analysis.

Autosomal dominant polycystic kidney disease (ADPKD) is the leading genetic cause of end stage renal disease characterized by progressive expansion of kidney cysts. To better understand the cell types and states driving ADPKD progression, we analyze eight ADPKD and five healthy human kidney samples, generating single cell multiomic atlas consisting of ~100,000 single nucleus transcriptomes and ~50,000 single nucleus epigenomes. Activation of proinflammatory, profibrotic signaling pathways are driven by proximal tubular cells with a failed repair transcriptomic signature, proinflammatory fibroblasts and collecting duct cells. We identify GPRC5A as a marker for cyst-lining collecting duct cells that exhibits increased transcription factor binding motif availability for NF-κB, TEAD, CREB and retinoic acid receptors. We identify and validate a distal enhancer regulating GPRC5A expression containing these motifs. This single cell multiomic analysis of human ADPKD reveals previously unrecognized cellular heterogeneity and provides a foundation to develop better diagnostic and therapeutic approaches.

Multimodal single cell sequencing implicates chromatin accessibility and genetic background in diabetic kidney disease progression.

The proximal tubule is a key regulator of kidney function and glucose metabolism. Diabetic kidney disease leads to proximal tubule injury and changes in chromatin accessibility that modify the activity of transcription factors involved in glucose metabolism and inflammation. Here we use single nucleus RNA and ATAC sequencing to show that diabetic kidney disease leads to reduced accessibility of glucocorticoid receptor binding sites and an injury-associated expression signature in the proximal tubule. We hypothesize that chromatin accessibility is regulated by genetic background and closely-intertwined with metabolic memory, which pre-programs the proximal tubule to respond differently to external stimuli. Glucocorticoid excess has long been known to increase risk for type 2 diabetes, which raises the possibility that glucocorticoid receptor inhibition may mitigate the adverse metabolic effects of diabetic kidney disease.

Spatially Resolved Transcriptomic Analysis of Acute Kidney Injury in a Female Murine Model.

BACKGROUND: Single-cell sequencing technologies have advanced our understanding of kidney biology and disease, but the loss of spatial information in these datasets hinders our interpretation of intercellular communication networks and regional gene expression patterns. New spatial transcriptomic sequencing platforms make it possible to measure the topography of gene expression at genome depth. METHODS: We optimized and validated a female bilateral ischemia-reperfusion injury model. Using the 10× Genomics Visium Spatial Gene Expression solution, we generated spatial maps of gene expression across the injury and repair time course, and applied two open-source computational tools, Giotto and SPOTlight, to increase resolution and measure cell-cell interaction dynamics. RESULTS: An ischemia time of 34 minutes in a female murine model resulted in comparable injury to 22 minutes for males. We report a total of 16,856 unique genes mapped across our injury and repair time course. Giotto, a computational toolbox for spatial data analysis, enabled increased resolution mapping of genes and cell types. Using a seeded nonnegative matrix regression (SPOTlight) to deconvolute the dynamic landscape of cell-cell interactions, we found that injured proximal tubule cells were characterized by increasing macrophage and lymphocyte interactions even 6 weeks after injury, potentially reflecting the AKI to CKD transition. CONCLUSIONS: In this transcriptomic atlas, we defined region-specific and injury-induced loss of differentiation markers and their re-expression during repair, as well as region-specific injury and repair transcriptional responses. Lastly, we created an interactive data visualization application for the scientific community to explore these results (http://humphreyslab.com/SingleCell/).

Recent advances in lineage tracing for the kidney.

Lineage tracing was originally developed by developmental biologists to identify all progeny of a single cell during morphogenesis. More recently this approach has been applied to other fields, including organ homeostasis and recovery from injury. Modern lineage tracing techniques typically rely on reporter gene expression induced by cell-specific DNA recombination. There have been important scientific advances in the last 10 years that have impacted lineage tracing approaches, including intersectional genetics, optical clearing techniques, and the use of sequencing-based genomic lineage tracing. The latter combines CRISPR-Cas9-based genetic scarring with single-cell RNA-sequencing that, in theory, could allow comprehensive reconstruction of a lineage tree for an entire organism. This review summarizes recent advances in lineage tracing technologies and outlines potential applications for kidney research.

Single cell transcriptional and chromatin accessibility profiling redefine cellular heterogeneity in the adult human kidney.

The integration of single cell transcriptome and chromatin accessibility datasets enables a deeper understanding of cell heterogeneity. We performed single nucleus ATAC (snATAC-seq) and RNA (snRNA-seq) sequencing to generate paired, cell-type-specific chromatin accessibility and transcriptional profiles of the adult human kidney. We demonstrate that snATAC-seq is comparable to snRNA-seq in the assignment of cell identity and can further refine our understanding of functional heterogeneity in the nephron. The majority of differentially accessible chromatin regions are localized to promoters and a significant proportion are closely associated with differentially expressed genes. Cell-type-specific enrichment of transcription factor binding motifs implicates the activation of NF-κB that promotes VCAM1 expression and drives transition between a subpopulation of proximal tubule epithelial cells. Our multi-omics approach improves the ability to detect unique cell states within the kidney and redefines cellular heterogeneity in the proximal tubule and thick ascending limb.

The autism-related protein CHD8 contributes to the stemness and differentiation of mouse hematopoietic stem cells.

Chromodomain helicase DNA-binding protein 8 (CHD8) is an ATP-dependent chromatin-remodeling factor that is encoded by the most frequently mutated gene in individuals with autism spectrum disorder. CHD8 is expressed not only in neural tissues but also in many other organs; however, its functions are largely unknown. Here, we show that CHD8 is highly expressed in and maintains the stemness of hematopoietic stem cells (HSCs). Conditional deletion of Chd8 specifically in mouse bone marrow induces cell cycle arrest, apoptosis, and a differentiation block in HSCs in association with upregulation of the expression of p53 target genes. A colony formation assay and bone marrow transplantation reveal that CHD8 deficiency also compromises the stemness of HSCs. Furthermore, additional ablation of p53 rescues the impaired stem cell function and differentiation block of CHD8-deficient HSCs. Our results thus suggest that the CHD8-p53 axis plays a key role in regulation of the stemness and differentiation of HSCs.

Skp2 contributes to cell cycle progression in trophoblast stem cells and to placental development.

All trophoblast subtypes of the placenta are derived from trophoblast stem cells (TSCs). TSCs have the capacity to self-renew, but how the proliferation of these cells is regulated in the undifferentiated state has been largely unclear. We now show that the F-box protein Skp2 regulates the proliferation of TSCs and thereby plays a pivotal role in placental development in mice on the C57BL/6 background. The placenta of Skp2-/- mouse embryos on the C57BL/6 background was smaller than that of their Skp2+/+ littermates, with the mutant embryos also manifesting intrauterine growth retardation. Although the Skp2-/- mice were born alive, most of them died before postnatal day 21, presumably as a result of placental defects. Depletion of Skp2 in TSCs cultured in the undifferentiated state resulted in a reduced rate of proliferation and arrest of the cell cycle in G1 phase, indicative of a defect in self-renewal capacity. The cell cycle arrest apparent in Skp2-deficient TSCs was reversed by additional ablation of the cyclin-dependent kinase inhibitor (CKI) p57 but not by that of the CKI p27. Our results thus suggest that Skp2-mediated degradation of p57 is an important determinant of the self-renewal capacity of TSCs during placental development, at least in mice of certain genetic backgrounds.

Disruption of FBXL5-mediated cellular iron homeostasis promotes liver carcinogenesis.

Hepatic iron overload is a risk factor for progression of hepatocellular carcinoma (HCC), although the molecular mechanisms underlying this association have remained unclear. We now show that the iron-sensing ubiquitin ligase FBXL5 is a previously unrecognized oncosuppressor in liver carcinogenesis in mice. Hepatocellular iron overload elicited by FBXL5 ablation gave rise to oxidative stress, tissue damage, inflammation, and compensatory proliferation of hepatocytes and to consequent promotion of liver carcinogenesis induced by exposure to a chemical carcinogen. The tumor-promoting outcome of FBXL5 deficiency in the liver was also found to be effective in a model of virus-induced HCC. FBXL5-deficient mice thus constitute the first genetically engineered mouse model of liver carcinogenesis promoted by iron overload. In addition, dysregulation of FBXL5-mediated cellular iron homeostasis was found to be associated with poor prognosis in human HCC, suggesting that FBXL5 plays a key role in defense against hepatocarcinogenesis.