Analysis of the combined effect of rs699 and rs5051 on angiotensinogen expression and hypertension.

Background: Hypertension (HTN) involves genetic variability in the renin-angiotensin system and influences antihypertensive response. We previously reported that angiotensinogen (AGT) messenger RNA (mRNA) is endogenously bound by miR-122-5p and rs699 A > G decreases reporter mRNA in the microRNA functional-assay PASSPORT-seq. The AGT promoter variant rs5051 C > T is in linkage disequilibrium (LD) with rs699 A > G and increases AGT transcription. The independent effect of these variants is understudied due to their LD therefore we aimed to test the hypothesis that increased AGT by rs5051 C > T counterbalances AGT decreased by rs699 A > G, and when these variants occur independently, it translates to HTN-related phenotypes. Methods: We used in silico, in vitro, in vivo, and retrospective models to test this hypothesis. Results: In silico, rs699 A > G is predicted to increase miR-122-5p binding affinity by 3%. Mir-eCLIP results show rs699 is 40-45 nucleotides from the strongest microRNA-binding site in the AGT mRNA. Unexpectedly, rs699 A > G increases AGT mRNA in an AGT-plasmid-cDNA HepG2 expression model. Genotype-Tissue Expression (GTEx) and UK Biobank analyses demonstrate liver AGT expression and HTN phenotypes are not different when rs699 A > G occurs independently from rs5051 C > T. However, GTEx and the in vitro experiments suggest rs699 A > G confers cell-type-specific effects on AGT mRNA abundance, and suggest paracrine renal renin-angiotensin-system perturbations could mediate the rs699 A > G associations with HTN. Conclusions: We found that rs5051 C > T and rs699 A > G significantly associate with systolic blood pressure in Black participants in the UK Biobank, demonstrating a fourfold larger effect than in White participants. Further studies are warranted to determine if altered antihypertensive response in Black individuals might be due to rs5051 C > T or rs699 A > G. Studies like this will help clinicians move beyond the use of race as a surrogate for genotype.

Prognostic implications of unbiased molecular categorization in kidney disease.

In this commentary, a novel approach to the reclassification of chronic kidney disease is reviewed. In the revisited study, the investigators identify 4 distinct subtypes of kidney disease derived from an unbiased self-organizing map of transcriptomic data from kidney biopsy samples. These molecular subtypes then are characterized by biologic cell processes, clinical and histopathologic features, urinary proteomics, and disease progression. The strengths and limitations of the self-organizing map approach are assessed; the prognostic, diagnostic, and therapeutic implications are considered briefly.

Spatial transcriptomics in health and disease.

The ability to localize hundreds of macromolecules to discrete locations, structures and cell types in a tissue is a powerful approach to understand the cellular and spatial organization of an organ. Spatially resolved transcriptomic technologies enable mapping of transcripts at single-cell or near single-cell resolution in a multiplex manner. The rapid development of spatial transcriptomic technologies has accelerated the pace of discovery in several fields, including nephrology. Its application to preclinical models and human samples has provided spatial information about new cell types discovered by single-cell sequencing and new insights into the cell-cell interactions within neighbourhoods, and has improved our understanding of the changes that occur in response to injury. Integration of spatial transcriptomic technologies with other omics methods, such as proteomics and spatial epigenetics, will further facilitate the generation of comprehensive molecular atlases, and provide insights into the dynamic relationships of molecular components in homeostasis and disease. This Review provides an overview of current and emerging spatial transcriptomic methods, their applications and remaining challenges for the field.

Nicotinamide riboside activates renal metabolism and protects the kidney in a model of Alport syndrome.

Chronic kidney disease (CKD) is associated with renal metabolic disturbances, including impaired fatty acid oxidation (FAO). Nicotinamide adenine dinucleotide (NAD + ) is a small molecule that participates in hundreds of metabolism-related reactions. NAD + levels are decreased in CKD, and NAD + supplementation is protective. However, both the mechanism of how NAD + supplementation protects from CKD, as well as the cell types most responsible, are poorly understood. Using a mouse model of Alport syndrome, we show that nicotinamide riboside (NR), an NAD + precursor, stimulates renal peroxisome proliferator-activated receptor α signaling and restores FAO in the proximal tubules, thereby protecting from CKD in both sexes. Bulk RNA-sequencing shows that renal metabolic pathways are impaired in Alport mice and dramatically activated by NR in both sexes. These transcriptional changes are confirmed by orthogonal imaging techniques and biochemical assays. Single nuclei RNA-sequencing and spatial transcriptomics, both the first of their kind from Alport mice, show that NAD + supplementation restores FAO in the proximal tubules with minimal effects on the podocytes. Finally, we also report, for the first time, sex differences at the transcriptional level in this Alport model. Male Alport mice had more severe inflammation and fibrosis than female mice at the transcriptional level. In summary, the data herein identify both the protective mechanism and location of NAD + supplementation in this model of CKD.

The chromatin landscape of healthy and injured cell types in the human kidney.

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.

Precision Medicine in Nephrology: An Integrative Framework of Multidimensional Data in the Kidney Precision Medicine Project.

Chronic kidney disease (CKD) and acute kidney injury (AKI) are heterogeneous syndromes defined clinically by serial measures of kidney function. Each condition possesses strong histopathologic associations, including glomerular obsolescence or acute tubular necrosis, respectively. Despite such characterization, there remains wide variation in patient outcomes and treatment responses. Precision medicine efforts, as exemplified by the Kidney Precision Medicine Project (KPMP), have begun to establish evolving, spatially anchored, cellular and molecular atlases of the cell types, states, and niches of the kidney in health and disease. The KPMP atlas provides molecular context for CKD and AKI disease drivers and will help define subtypes of disease that are not readily apparent from canonical functional or histopathologic characterization but instead are appreciable through advanced clinical phenotyping, pathomic, transcriptomic, proteomic, epigenomic, and metabolomic interrogation of kidney biopsy samples. This perspective outlines the structure of the KPMP, its approach to the integration of these diverse datasets, and its major outputs relevant to future patient care.

Dimension-agnostic and granularity-based spatially variable gene identification using BSP.

Identifying spatially variable genes (SVGs) is critical in linking molecular cell functions with tissue phenotypes. Spatially resolved transcriptomics captures cellular-level gene expression with corresponding spatial coordinates in two or three dimensions and can be used to infer SVGs effectively. However, current computational methods may not achieve reliable results and often cannot handle three-dimensional spatial transcriptomic data. Here we introduce BSP (big-small patch), a non-parametric model by comparing gene expression pattens at two spatial granularities to identify SVGs from two or three-dimensional spatial transcriptomics data in a fast and robust manner. This method has been extensively tested in simulations, demonstrating superior accuracy, robustness, and high efficiency. BSP is further validated by substantiated biological discoveries in cancer, neural science, rheumatoid arthritis, and kidney studies with various types of spatial transcriptomics technologies.

Inflammation primes the kidney for recovery by activating AZIN1 A-to-I editing.

The progression of kidney disease varies among individuals, but a general methodology to quantify disease timelines is lacking. Particularly challenging is the task of determining the potential for recovery from acute kidney injury following various insults. Here, we report that quantitation of post-transcriptional adenosine-to-inosine (A-to-I) RNA editing offers a distinct genome-wide signature, enabling the delineation of disease trajectories in the kidney. A well-defined murine model of endotoxemia permitted the identification of the origin and extent of A-to-I editing, along with temporally discrete signatures of double-stranded RNA stress and Adenosine Deaminase isoform switching. We found that A-to-I editing of Antizyme Inhibitor 1 (AZIN1), a positive regulator of polyamine biosynthesis, serves as a particularly useful temporal landmark during endotoxemia. Our data indicate that AZIN1 A-to-I editing, triggered by preceding inflammation, primes the kidney and activates endogenous recovery mechanisms. By comparing genetically modified human cell lines and mice locked in either A-to-I edited or uneditable states, we uncovered that AZIN1 A-to-I editing not only enhances polyamine biosynthesis but also engages glycolysis and nicotinamide biosynthesis to drive the recovery phenotype. Our findings implicate that quantifying AZIN1 A-to-I editing could potentially identify individuals who have transitioned to an endogenous recovery phase. This phase would reflect their past inflammation and indicate their potential for future recovery.

Clinical and Histopathologic Characteristics of Pediatric Patients With Primary Membranous Nephropathy.

Introduction: Primary membranous nephropathy (PMN) is uncommon in children. Therefore, data on the clinical course of affected children are scarce. In recent years, several novel antigens have been implicated in the pathogenesis of PMN. However, the histopathologic characteristics of pediatric patients with PMN remain poorly represented in the literature. Methods: We have retrospectively analyzed the clinical presentation and outcomes data of 21 children with PMN from 3 centers in the United States. In addition, we have identified novel antigens in biopsy specimens from these patients and correlated their presence or absence to clinical outcomes. Finally, we compared the results of the novel antigen staining from our clinical cohort to a validation cohort of 127 biopsy specimens from children with PMN at Arkana Laboratories. Results: The data from the 2 cohorts demonstrated similar overall antigen positivity rates of 62% to 63%, with phospholipase A2 receptor (PLA2R) and exostosin 1 (EXT1) being the most commonly found antigens. Results from the clinical cohort showed that overall, the kidney prognosis for children with PMN was good, with 17 of 21 patients entering a complete or partial remission. Children who were positive for PLA2R or EXT1 were significantly more likely to enter remission than those in the antigen negative group. Conclusion: Approximately 60% of pediatric membranous cases are positive for a novel antigen on kidney biopsy and the clinical prognosis is generally favorable. More studies are needed to understand the clinical implications of each specific novel antigen.

Proteomic analysis of murine kidney proximal tubule sub-segment derived cell lines reveals preferences in mitochondrial pathway activity.

The proximal tubule (PT) is a nephron segment that is responsible for the majority of solute and water reabsorption in the kidney. Each of its sub-segments have specialized functions; however, little is known about the genes and proteins that determine the oxidative phosphorylation capacity of the PT sub-segments. This information is critical to understanding kidney function and will provide a comprehensive landscape of renal cell adaptations to injury, physiologic stressors, and development. This study analyzed three immortalized murine renal cell lines (PT S1, S2, and S3 segments) for protein content and compared them to a murine fibroblast cell line. All three proximal tubule cell lines generate ATP predominantly by oxidative phosphorylation while the fibroblast cell line is glycolytic. The proteomic data demonstrates that the most significant difference in proteomic signatures between the cell lines are proteins known to be localized in the nucleus followed by mitochondrial proteins. Mitochondrial metabolic substrate utilization assays were performed using the proximal tubule cell lines to determine substrate utilization kinetics thereby providing a physiologic context to the proteomic dataset. This data will allow researchers to study differences in nephron-specific cell lines, between epithelial and fibroblast cells, and between actively respiring cells and glycolytic cells. SIGNIFICANCE: Proteomic analysis of proteins expressed in immortalized murine renal proximal tubule cells was compared to a murine fibroblast cell line proteome. The proximal tubule segment specific cell lines: S1, S2 and S3 are all grown under conditions whereby the cells generate ATP by oxidative phosphorylation while the fibroblast cell line utilizes anaerobic glycolysis for ATP generation. The proteomic studies allow for the following queries: 1) comparisons between the proximal tubule segment specific cell lines, 2) comparisons between polarized epithelia and fibroblasts, 3) comparison between cells employing oxidative phosphorylation versus anaerobic glycolysis and 4) comparisons between cells grown on clear versus opaque membrane supports. The data finds major differences in nuclear protein expression and mitochondrial proteins. This proteomic data set will be an important baseline dataset for investigators who need immortalized renal proximal tubule epithelial cells for their research.

The INGENIOUS trial: Impact of pharmacogenetic testing on adverse events in a pragmatic clinical trial.

Adverse drug events (ADEs) account for a significant mortality, morbidity, and cost burden. Pharmacogenetic testing has the potential to reduce ADEs and inefficacy. The objective of this INGENIOUS trial (NCT02297126) analysis was to determine whether conducting and reporting pharmacogenetic panel testing impacts ADE frequency. The trial was a pragmatic, randomized controlled clinical trial, adapted as a propensity matched analysis in individuals (N = 2612) receiving a new prescription for one or more of 26 pharmacogenetic-actionable drugs across a community safety-net and academic health system. The intervention was a pharmacogenetic testing panel for 26 drugs with dosage and selection recommendations returned to the health record. The primary outcome was occurrence of ADEs within 1 year, according to modified Common Terminology Criteria for Adverse Events (CTCAE). In the propensity-matched analysis, 16.1% of individuals experienced any ADE within 1-year. Serious ADEs (CTCAE level ≥ 3) occurred in 3.2% of individuals. When combining all 26 drugs, no significant difference was observed between the pharmacogenetic testing and control arms for any ADE (Odds ratio 0.96, 95% CI: 0.78-1.18), serious ADEs (OR: 0.91, 95% CI: 0.58-1.40), or mortality (OR: 0.60, 95% CI: 0.28-1.21). However, sub-group analyses revealed a reduction in serious ADEs and death in individuals who underwent pharmacogenotyping for aripiprazole and serotonin or serotonin-norepinephrine reuptake inhibitors (OR 0.34, 95% CI: 0.12-0.85). In conclusion, no change in overall ADEs was observed after pharmacogenetic testing. However, limitations incurred during INGENIOUS likely affected the results. Future studies may consider preemptive, rather than reactive, pharmacogenetic panel testing.

miRNA and mRNA Signatures in Human Acute Kidney Injury Tissue.

Acute kidney injury (AKI) is an important contributor to the development of chronic kidney disease (CKD). There is a need to understand molecular mediators that drive either recovery or progression to CKD. In particular, the role of miRNA and its regulatory role in AKI is poorly understood. We performed miRNA and mRNA sequencing on biobanked human kidney tissues obtained in the routine clinical care of patients with the diagnoses of AKI and minimal change disease (MCD), in addition to nephrectomized (Ref) tissue from individuals without known kidney disease. Transcriptomic analysis of mRNA revealed that Ref tissues exhibited a similar injury signature to AKI, not identified in MCD samples. The transcriptomic signature of human AKI was enriched with genes in pathways involved in cell adhesion and epithelial-to-mesenchymal transition (e.g., CDH6, ITGB6, CDKN1A ). miRNA DE analysis revealed upregulation of miRNA associated with immune cell recruitment and inflammation (e.g., miR-146a, miR-155, miR-142, miR-122). These miRNA (i.e., miR-122, miR-146) are also associated with downregulation of mRNA such as DDR2 and IGFBP6 , respectively. These findings suggest integrated interactions between miRNAs and target mRNAs in AKI-related processes such as inflammation, immune cell activation and epithelial-to-mesenchymal transition. These data contribute several novel findings when describing the epigenetic regulation of AKI by miRNA, and also underscores the importance of utilizing an appropriate reference control tissue to understand canonical pathway alterations in AKI.

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

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

A spatially anchored transcriptomic atlas of the human kidney papilla identifies significant immune injury in patients with stone disease.

Kidney stone disease causes significant morbidity and increases health care utilization. In this work, we decipher the cellular and molecular niche of the human renal papilla in patients with calcium oxalate (CaOx) stone disease and healthy subjects. In addition to identifying cell types important in papillary physiology, we characterize collecting duct cell subtypes and an undifferentiated epithelial cell type that was more prevalent in stone patients. Despite the focal nature of mineral deposition in nephrolithiasis, we uncover a global injury signature characterized by immune activation, oxidative stress and extracellular matrix remodeling. We also identify the association of MMP7 and MMP9 expression with stone disease and mineral deposition, respectively. MMP7 and MMP9 are significantly increased in the urine of patients with CaOx stone disease, and their levels correlate with disease activity. Our results define the spatial molecular landscape and specific pathways contributing to stone-mediated injury in the human papilla and identify associated urinary biomarkers.

An atlas of healthy and injured cell states and niches in the human kidney.

Understanding kidney disease relies on defining the complexity of cell types and states, their associated molecular profiles and interactions within tissue neighbourhoods1. Here we applied multiple single-cell and single-nucleus assays (>400,000 nuclei or cells) and spatial imaging technologies to a broad spectrum of healthy reference kidneys (45 donors) and diseased kidneys (48 patients). This has provided a high-resolution cellular atlas of 51 main cell types, which include rare and previously undescribed cell populations. The multi-omic approach provides detailed transcriptomic profiles, regulatory factors and spatial localizations spanning the entire kidney. We also define 28 cellular states across nephron segments and interstitium that were altered in kidney injury, encompassing cycling, adaptive (successful or maladaptive repair), transitioning and degenerative states. Molecular signatures permitted the localization of these states within injury neighbourhoods using spatial transcriptomics, while large-scale 3D imaging analysis (around 1.2 million neighbourhoods) provided corresponding linkages to active immune responses. These analyses defined biological pathways that are relevant to injury time-course and niches, including signatures underlying epithelial repair that predicted maladaptive states associated with a decline in kidney function. This integrated multimodal spatial cell atlas of healthy and diseased human kidneys represents a comprehensive benchmark of cellular states, neighbourhoods, outcome-associated signatures and publicly available interactive visualizations.

The chromatin landscape of healthy and injured cell types in the human kidney.

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. However, comprehensive integration of gene expression with epigenetic features that define regulatory elements remains a significant challenge. We measured 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 established a comprehensive and spatially-anchored epigenomic atlas to define the kidney's active, silent, and regulatory accessible chromatin regions across the genome. Using this atlas, we noted distinct control of adaptive injury in different epithelial cell types. A proximal tubule cell transcription factor network of ELF3 , KLF6 , and KLF10 regulated the transition between health and injury, while in thick ascending limb cells this transition was regulated by NR2F1 . Further, combined perturbation of ELF3 , KLF6 , and KLF10 distinguished 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.

Analysis of the Combined Effect of rs699 and rs5051 on Angiotensinogen Expression and Hypertension.

Hypertension (HTN) involves genetic variability in the renin-angiotensin system and characterizing this variability will help advance precision antihypertensive treatments. We previously reported that angiotensinogen (AGT) mRNA is endogenously bound by mir-122-5p and that rs699 A>G significantly decreases reporter mRNA in the functional mirSNP assay PASSPORT-seq. The AGT promoter variant rs5051 C>T is in linkage disequilibrium (LD) with rs699 A>G and increases AGT transcription. We hypothesized that the increased AGT by rs5051 C>T counterbalances AGT decrease by rs699 A>G, and when these variants occur independently, would translate to HTN-related phenotypes. The independent effect of each of these variants is understudied due to their LD, therefore, we used in silico, in vitro, in vivo, and retrospective clinical and biobank analyses to assess HTN and AGT expression phenotypes where rs699 A>G occurs independently from rs5051 C>T. In silico, rs699 A>G is predicted to increase mir-122-5p binding strength by 3%. Mir-eCLIP assay results show that rs699 is 40-45 nucleotides from the strongest microRNA binding site in the AGT mRNA. Unexpectedly, rs699 A>G increases AGT mRNA in a plasmid cDNA HepG2 expression model. GTEx and UK Biobank analyses demonstrate that liver AGT expression and HTN phenotypes were not different when rs699 A>G occurs independently from rs5051 C>T, allowing us to reject the original hypothesis. However, both GTEx and our in vitro experiments suggest rs699 A>G confers cell-type specific effects on AGT mRNA abundance. We found that rs5051 C>T and rs699 A>G significantly associate with systolic blood pressure in Black participants in the UK Biobank, demonstrating a 4-fold larger effect than in White participants. Further studies are warranted to determine if the altered antihypertensive response in Black individuals might be due to rs5051 C>T or rs699 A>G. Studies like this will help clinicians move beyond the use of race as a surrogate for genotype.

Implementation of Clinical Cytochrome P450 3A Genotyping for Tacrolimus Dosing in a Large Kidney Transplant Program.

Tacrolimus is a calcineurin inhibitor with a narrow therapeutic range and is metabolized by cytochrome P450 (CYP) isoenzymes CYP3A4 and CYP3A5. The Clinical Pharmacogenetic Implementation Consortium published evidence-based guidelines for CYP3A5 normal/intermediate metabolizers prescribed tacrolimus, yet few transplant centers have implemented routine testing. The objective of this study was to implement preemptive CYP3A genotyping into clinical practice in a large kidney transplant program and to evaluate workflow feasibility, potential clinical benefit, and reimbursement to identify barriers and determine sustainability. Preemptive pharmacogenetic testing for CYP3A5 and CYP3A4 was implemented in all patients listed for a kidney transplant as part of standard clinical care. Genotyping was performed at the listing appointment, results were reported as discrete data in the electronic medical record, and education and clinical decision support alerts were developed to provide pharmacogenetic-recommended tacrolimus dosing. During this initial phase, all patients were administered standard tacrolimus dosing, and clinical and reimbursement outcomes were collected. Greater than 99.5% of genotyping claims were reimbursed by third-party payers. CYP3A5 normal/intermediate metabolizers had significantly fewer tacrolimus trough concentrations within the target range and a significantly longer time to their first therapeutic trough compared to poor metabolizers. The challenge of tacrolimus dosing is magnified in the African American population. The US Food and Drug Administration drug label recommends increased starting doses in African ancestry, yet only ≈66% of African Americans in our cohort were normal/intermediate metabolizers who required higher doses. Routine CYP3A5 genotyping may overcome this issue by using genotype over race as a more accurate predictor of drug response.

Integrated Cytometry With Machine Learning Applied to High-Content Imaging of Human Kidney Tissue for In Situ Cell Classification and Neighborhood Analysis.

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.