Introduction The maintenance of a healthy epithelial-endothelial juxtaposition requires crosstalk within glomerular cellular niches. Here, we sought to understand the spatially anchored regulation and transition of endothelial and mesangial cells from health to injury in DKD. Methods From 132 human kidney samples, an integrated multiomics approach was leveraged to identify cellular niches, cell-cell communication, cell injury trajectories, and regulatory transcription factor networks in glomerular capillary endothelial (EC-GC) and mesangial cells. Data were culled from single nucleus RNA, ATAC-sequencing and four orthogonal spatial transcriptomic technologies for correlation with histopathological and clinical trial data. Results We identified a cellular niche in diabetic glomeruli enriched in a proliferative endothelial cell subtype (prEC) and altered vascular smooth muscle cells (VSMCs). Cellular communication within this niche maintained pro-angiogenic signaling with loss of anti-angiogenic factors. We identified a transcription factor network of MEF2C, MEF2A, and TRPS1 which regulated SEMA6A and PLXNA2, a receptor-ligand pair opposing angiogenesis. In silico knockout of the transcription factor network accelerated the transition from healthy EC-GCs toward a degenerative injured endothelial phenotype, with concomitant disruption of EC-GC and prEC expression patterns. Glomeruli enriched in the prEC niche had histologic evidence of neovascularization. MEF2C activity was increased in diabetic glomeruli with nodular mesangial sclerosis. The gene regulatory network (GRN) of MEF2C was dysregulated in EC-GCs of patients with DKD, but sodium glucose transporter-2 inhibitor (SGLT2i) treatment reversed the MEF2C GRN effects of DKD. Conclusions The MEF2C, MEF2A, and TRPS1 transcription factor network carefully balances the fate of the EC-GC in DKD. When the transcription factor network is “on” or over-expressed in DKD, EC-GCs may progress to a prEC state, while transcription factor suppression leads to cell death. SGLT2i therapy may restore the balance of MEF2C activity.
The organizational principles of nephronal segments are based on anatomical and physiological attributes that are linked to the homeostatic functions of the kidney. Recent molecular approaches have uncovered layers of deeper signatures and states in tubular cells that arise at various time points on the disease trajectory. Here, we introduce an analytical pipeline of multiplexed spatial protein imaging integrated with RNA expression to characterize proximal tubular subpopulations and neighborhoods in human kidney tissue. We demonstrate that, in reference tissue, a large proportion of S1 proximal tubular epithelial cells expresses thymus antigen 1 (THY1), a mesenchymal stromal and stem cell marker that regulates differentiation. Kidney disease is associated with loss of THY1 and transition toward expression of prominin 1 (PROM1), another stem cell marker recently linked to failed repair. Our data support a model in which the interplay between THY1 and PROM1 expression in proximal tubules associates with their regenerative potential and marks the timeline of disease progression.
There is a need to define regions of gene activation or repression that control human kidney cells in states of health, injury, and repair to understand the molecular pathogenesis of kidney disease and design therapeutic strategies. Comprehensive integration of gene expression with epigenetic features that define regulatory elements remains a significant challenge. We measure dual single nucleus RNA expression and chromatin accessibility, DNA methylation, and H3K27ac, H3K4me1, H3K4me3, and H3K27me3 histone modifications to decipher the chromatin landscape and gene regulation of the kidney in reference and adaptive injury states. We establish a spatially-anchored epigenomic atlas to define the kidney's active, silent, and regulatory accessible chromatin regions across the genome. Using this atlas, we note distinct control of adaptive injury in different epithelial cell types. A proximal tubule cell transcription factor network of ELF3, KLF6, and KLF10 regulates the transition between health and injury, while in thick ascending limb cells this transition is regulated by NR2F1. Further, combined perturbation of ELF3, KLF6, and KLF10 distinguishes two adaptive proximal tubular cell subtypes, one of which manifested a repair trajectory after knockout. This atlas will serve as a foundation to facilitate targeted cell-specific therapeutics by reprogramming gene regulatory networks.
Long, Anne; Zollman, Amy; Hato, Takashi; Sutton, Timothy A.; Overstreet, Jessica M.; Dagher, Pierre C. Author Information
The Kidney Precision Medicine Project (KPMP) aims to create a kidney tissue atlas, define disease subgroups, and identify critical cells, pathways, and targets for novel therapies through molecular investigation of human kidney biopsies obtained from participants with AKI or CKD. We present the case of a 66-year-old woman with diabetic kidney disease who underwent a protocol KPMP kidney biopsy. Her clinical history included diabetes mellitus complicated by neuropathy and eye disease, increased insulin resistance, hypertension, albuminuria, and relatively preserved glomerular filtration rate (early CKD stage 3a). The patient's histopathology was consistent with diabetic nephropathy and arterial and arteriolar sclerosis. Three-dimensional, immunofluorescence imaging of the kidney biopsy specimen revealed extensive periglomerular neovascularization that was underestimated by standard histopathologic approaches. Spatial transcriptomics was performed to obtain gene expression signatures at discrete areas of the kidney biopsy. Gene expression in the areas of glomerular neovascularization revealed increased expression of genes involved in angiogenic signaling, proliferation, and survival of endothelial cells, as well as new vessel maturation and stability. This molecular correlation provides additional insights into the development of kidney disease in patients with diabetes and spotlights how novel molecular techniques used by the KPMP can supplement and enrich the histopathologic diagnosis obtained from a kidney biopsy.
The human kidney is a complex organ with various cell types that are intricately organized to perform key physiological functions and maintain homeostasis. New imaging modalities, such as mesoscale and highly multiplexed fluorescence microscopy, are increasingly being applied to human kidney tissue to create single-cell resolution data sets that are both spatially large and multidimensional. These single-cell resolution high-content imaging data sets have great potential to uncover the complex spatial organization and cellular makeup of the human kidney. Tissue cytometry is a novel approach used for the quantitative analysis of imaging data; however, the scale and complexity of such data sets pose unique challenges for processing and analysis. We have developed the Volumetric Tissue Exploration and Analysis (VTEA) software, a unique tool that integrates image processing, segmentation, and interactive cytometry analysis into a single framework on desktop computers. Supported by an extensible and open-source framework, VTEA’s integrated pipeline now includes enhanced analytical tools, such as machine learning, data visualization, and neighborhood analyses, for hyperdimensional large-scale imaging data sets. These novel capabilities enable the analysis of mesoscale 2- and 3-dimensional multiplexed human kidney imaging data sets (such as co-detection by indexing and 3-dimensional confocal multiplexed fluorescence imaging). We demonstrate the utility of this approach in identifying cell subtypes in the kidney on the basis of labels, spatial association, and their microenvironment or neighborhood membership. VTEA provides an integrated and intuitive approach to decipher the cellular and spatial complexity of the human kidney and complements other transcriptomics and epigenetic efforts to define the landscape of kidney cell types.
Efforts toward achieving diversity, equity, inclusion, and justice (DEIJ) within graduate medical education (GME) often begin with the formation of a DEIJ committee that steers the work. Little is known about the experiences and the challenges faced by those serving on such committees. We sought to describe the experiences of members of our institutional GME DEIJ committee to gain knowledge that would propel this work forward. An open-ended survey was electronically administered to members of our institutional GME DEIJ committee. Responses were analyzed using a rapid qualitative analytical approach. Eighteen members (58%) responded. Of these, (67%) were women and five (28%) were Black. Six domains emerged: "motivation," "challenges," "emotional response," "highs," "facilitators," and "advice." Black respondents more often cited the need to increase diversity as a motivator to join this work. Women and Black respondents more often identified time constraints as a challenge to participation. Some members found the work emotionally draining; others described it as uplifting. Two themes emerged as high points of participation-pride and achievement around the work completed and the personal benefits of building a community with a shared purpose. Three themes emerged as facilitators: effective leadership, support, and establishing psychological safety during the meetings. Many arrived at the realization that change would take time and advocated for patience and perseverance. Protected time and DEIJ expertise were identified as integral to successful committee work. Our findings provide novel insights into the experience of serving on a GME DEIJ committee and highlights infrastructural and institutional prerequisites for success.
Diabetic kidney disease (DKD) remains the leading cause of end-stage kidney disease despite decades of study. Alterations in the glomerulus and kidney tubules both contribute to the pathogenesis of DKD although the majority of investigative efforts have focused on the glomerulus. We sought to examine the differential expression signature of human DKD in the glomerulus and proximal tubule and corroborate our findings in the db/db mouse model of diabetes. A transcriptogram network analysis of RNAseq data from laser microdissected (LMD) human glomerulus and proximal tubule of DKD and reference nephrectomy samples revealed enriched pathways including rhodopsin-like receptors, olfactory signaling, and ribosome (protein translation) in the proximal tubule of human DKD biopsy samples. The translation pathway was also enriched in the glomerulus. Increased translation in diabetic kidneys was validated using polyribosomal profiling in the db/db mouse model of diabetes. Using single nuclear RNA sequencing (snRNAseq) of kidneys from db/db mice, we prioritized additional pathways identified in human DKD. The top overlapping pathway identified in the murine snRNAseq proximal tubule clusters and the human LMD proximal tubule compartment was carboxylic acid catabolism. Using ultra-performance liquid chromatography-mass spectrometry, the fatty acid catabolism pathway was also found to be dysregulated in the db/db mouse model. The Acetyl-CoA metabolite was down-regulated in db/db mice, aligning with the human differential expression of the genes ACOX1 and ACACB. In summary, our findings demonstrate that proximal tubular alterations in protein translation and carboxylic acid catabolism are key features in both human and murine DKD.
Despite important advances in studying experimental and clinical acute kidney injury (AKI), the pathogenesis of this disease remains incompletely understood. Single cell sequencing studies have closed this knowledge gap by characterizing the transcriptomic signature of different cell types within the kidney. However, the spatial distribution of injury can be regional and affect cells heterogeneously. We first optimized coordination of spatial transcriptomics and single nuclear sequencing datasets, mapping 30 dominant cell types to a human nephrectomy sample. The predicted cell type spots corresponded with the underlying hematoxylin and eosin histopathology. To study the implications of acute kidney injury on the distribution of transcript expression, we then characterized the spatial transcriptomic signature of two murine AKI models: ischemia reperfusion injury (IRI) and cecal ligation puncture (CLP). Localized regions of reduced overall expression were found associated with tissue injury pathways. Using single cell sequencing, we deconvoluted the signature of each spatial transcriptomic spot, identifying patterns of colocalization between immune and epithelial cells. As expected, neutrophils infiltrated the renal medullary outer stripe in the ischemia model. Atf3 was identified as a chemotactic factor in S3 proximal tubule cells. In the CLP model, infiltrating macrophages dominated the outer cortical signature and Mdk was identified as a corresponding chemotactic factor. The regional distribution of these immune cells was validated with multiplexed CO-Detection by inDEXing (CODEX) immunofluorescence. Spatial transcriptomic sequencing can aid in uncovering the mechanisms driving immune cell infiltration and allow detection of relevant subpopulations in single cell sequencing. The complementarity of these technologies facilitates the development of a transcriptomic kidney atlas in health and disease.
Single-cell sequencing studies have characterized the transcriptomic signature of cell types within the kidney. However, the spatial distribution of acute kidney injury (AKI) is regional and affects cells heterogeneously. We first optimized coordination of spatial transcriptomics and single-nuclear sequencing data sets, mapping 30 dominant cell types to a human nephrectomy. The predicted cell-type spots corresponded with the underlying histopathology. To study the implications of AKI on transcript expression, we then characterized the spatial transcriptomic signature of 2 murine AKI models: ischemia/reperfusion injury (IRI) and cecal ligation puncture (CLP). Localized regions of reduced overall expression were associated with injury pathways. Using single-cell sequencing, we deconvoluted the signature of each spatial transcriptomic spot, identifying patterns of colocalization between immune and epithelial cells. Neutrophils infiltrated the renal medulla in the ischemia model. Atf3 was identified as a chemotactic factor in S3 proximal tubules. In the CLP model, infiltrating macrophages dominated the outer cortical signature, and Mdk was identified as a corresponding chemotactic factor. The regional distribution of these immune cells was validated with multiplexed CO-Detection by indEXing (CODEX) immunofluorescence. Spatial transcriptomic sequencing complemented single-cell sequencing by uncovering mechanisms driving immune cell infiltration and detection of relevant cell subpopulations.
The gene expression signature of the human kidney interstitium is incompletely understood. The cortical interstitium (excluding tubules, glomeruli, and vessels) in reference nephrectomies ( N = 9) and diabetic kidney biopsy specimens ( N = 6) was laser microdissected (LMD) and sequenced. Samples underwent RNA sequencing. Gene signatures were deconvolved using single nuclear RNA sequencing (snRNAseq) data derived from overlapping specimens. Interstitial LMD transcriptomics uncovered previously unidentified markers including KISS1 , validated with in situ hybridization. LMD transcriptomics and snRNAseq revealed strong correlation of gene expression within corresponding kidney regions. Relevant enriched interstitial pathways included G-protein coupled receptor. binding and collagen biosynthesis. The diabetic interstitium was enriched for extracellular matrix organization and small-molecule catabolism. Cell type markers with unchanged expression ( NOTCH3 , EGFR , and HEG1 ) and those down-regulated in diabetic nephropathy ( MYH11 , LUM , and CCDC3 ) were identified. LMD transcriptomics complements snRNAseq; together, they facilitate mapping of interstitial marker genes to aid interpretation of pathophysiology in precision medicine studies.
The advent of personalized medicine has driven the development of novel approaches for obtaining detailed cellular and molecular information from clinical tissue samples. Tissue cytometry is a promising new technique that can be used to enumerate and characterize each cell in a tissue and, unlike flow cytometry and other single-cell techniques, does so in the context of the intact tissue, preserving spatial information that is frequently crucial to understanding a cell's physiology, function, and behavior. However, the wide-scale adoption of tissue cytometry as a research tool has been limited by the fact that published examples utilize specialized techniques that are beyond the capabilities of most laboratories. Here we describe a complete and accessible pipeline, including methods of sample preparation, microscopy, image analysis, and data analysis for large-scale three-dimensional tissue cytometry of human kidney tissues. In this workflow, multiphoton microscopy of unlabeled tissue is first conducted to collect autofluorescence and second-harmonic images. The tissue is then labeled with eight fluorescent probes, and imaged using spectral confocal microscopy. The raw 16-channel images are spectrally deconvolved into 8-channel images, and analyzed using the Volumetric Tissue Exploration and Analysis (VTEA) software developed by our group. We applied this workflow to analyze millimeter-scale tissue samples obtained from human nephrectomies and from renal biopsies from individuals diagnosed with diabetic nephropathy, generating a quantitative census of tens of thousands of cells in each. Such analyses can provide useful insights that can be linked to the biology or pathology of kidney disease. The approach utilizes common laboratory techniques, is compatible with most commercially-available confocal microscope systems and all image and data analysis is conducted using the VTEA image analysis software, which is available as a plug-in for ImageJ.