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.

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.

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.

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.

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.

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.

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.

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.

A conditionally immortalized Gli1-positive kidney mesenchymal cell line models myofibroblast transition.

Glioma-associated oncogene homolog-1 (Gli1)-positive resident mesenchymal stem cell-like cells are the predominant source of kidney myofibroblasts in fibrosis, but investigating Gli1-positive myofibroblast progenitor activation is hampered by the difficulty of isolating and propagating primary cultures of these cells. Using a genetic strategy with positive and negative selection, we isolated Kidney-Gli1 (KGli1) cells that maintain expression of appropriate mesenchymal stem cell-like cell markers, respond to hedgehog pathway activation, and display robust myofibroblast differentiation upon treatment with transforming growth factor-ß (TGF-ß). Coculture of KGli1 cells with endothelium stabilizes capillary formation. Single-cell RNA sequencing (scRNA-seq) analysis during differentiation identified autocrine ligand-receptor pair upregulation and a strong focal adhesion pathway signal. This led us to test the serum response factor inhibitor CCG-203971 that potently inhibited TGF-ß-induced pericyte-to-myofibroblast transition. scRNA-seq also identified the unexpected upregulation of nerve growth factor (NGF), which we confirmed in two mouse kidney fibrosis models. The Ngf receptor Ntrk1 is expressed in tubular epithelium in vivo, suggesting a novel interstitial-to-tubule paracrine signaling axis. Thus, KGli1 cells accurately model myofibroblast activation in vitro, and the development of this cell line provides a new tool to study resident mesenchymal stem cell-like progenitors in health and disease.

FBXL12 regulates T-cell differentiation in a cell-autonomous manner.

Aldehyde dehydrogenase (ALDH) activity is a hallmark of stem cells including embryonic, adult tissue and cancer stem cells. The SCF(FBXL) (12) complex is an authentic ubiquitin ligase that targets ALDH3 for degradation. FBXL12 is essential for the differentiation of trophoblast stem cells into specific cell types in the placenta during mouse embryogenesis, but its physiological functions in adult tissues have remained unknown. We have now investigated the role of the FBXL12-ALDH3 axis in the thymus, in which FBXL12 was most abundant among adult mouse tissues examined. During T-cell differentiation, FBXL12 is most abundant in CD4(+) CD8(+) (DP) cells, with its expression declining as these cells differentiate into CD4(+) CD8(-) or CD4(-) CD8(+) (SP) cells. T cells of FBXL12-null mice manifested a differentiation block at the DP-SP transition that was associated with ALDH3 accumulation in DP cells. This differentiation block was also apparent in wild-type mouse recipients of FBXL12-null bone marrow transplants as well as in FBXL12-null fetal thymic organ culture, suggesting that it is a cell-autonomous phenomenon in the thymus rather than an indirect effect of altered systemic conditions. Our results thus indicate that, in addition to its role in placental development, the FBXL12-ALDH3 axis is required for maturation of undifferentiated thymocytes.

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.

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.

Injury-induced FoxM1 expression in mouse kidney drives epithelial proliferation by a Cyclin F dependent mechanism.

Acute kidney injury strongly upregulates the transcription factor Foxm1 in 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 reduces proximal tubule proliferation after ischemic injury in vivo. Foxm1 deletion led to increased AKI-to-CKD transition with enhanced fibrosis and ongoing tubule injury 6 weeks after injury. We report extracellular signal-regulated kinase (ERK) mediates FOXM1 induction downstream of the epidermal growth factor receptor (EGFR) in primary proximal tubule cells. We defined FOXM1 genomic binding sites by Cleavage Under Targets & 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 datasets revealed the cell cycle regulator cyclin F (CCNF) as a putative FOXM1 target. We identify two cis regulatory elements that bind FOXM1 and regulate CCNF expression, demonstrate 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.

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.

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.

Repeated Epstein-Barr virus-associated hemophagocytic syndrome in a lupus nephritis patient.

A 55-year-old Japanese female was diagnosed with systemic lupus erythematosus (SLE) and developed nephrotic syndrome. She was diagnosed with lupus nephritis by a percutaneous renal biopsy. She was treated with intravenous steroid pulse therapy twice, but it proved to be ineffective. She achieved a complete remission after intravenous cyclophosphamide pulse (CPAIV) therapy. Thereafter, her lupus nephritis was well controlled and demonstrated only a low activity. However, she suffered Epstein- Barr virus (EBV)-associated hemophagocytic syndrome (HPS) twice, and in each case she was treated with anticancer drugs and achieved a complete remission. This was a rare case of lupus nephritis who showed repeated EBV-associated HPS.

Development of Wernicke's encephalopathy during initiation of hemodialysis in an elderly non-alcoholic patient.

A 79-year-old man with chronic renal failure developed general fatigue and loss of appetite. He was diagnosed with endstage renal disease and was started on hemodialysis (HD). The symptoms improved immediately, but the mental status deteriorated gradually, reaching Glasgow Coma Scale (GCS) 5. Computed tomography showed no significant intracranial lesion, but magnetic resonance images showed symmetric high-intensity changes in the periaqueductal area, suggestive of Wernicke's encephalopathy (WE). He was immediately treated with intravenous infusion of thiamine. Five days later, the mental status level improved up to GCS 14, and the above MRI findings disappeared. To our knowledge, this is the first report describing the clinical outcome of a non-alcoholic patient who developed WE during initiation of HD. WE should be suspected in patients who are on chronic HD as well as those on initiation of HD with unexplained neurological abnormalities.

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.

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.