Plasma Proteins associated with Chronic Histopathologic Lesions on Kidney Biopsy

BACKGROUND

The severity of chronic histopathologic lesions on kidney biopsy is independently associated with higher risk of progressive chronic kidney disease (CKD). Because kidney biopsies are invasive, identification of blood markers that report on underlying kidney histopathology has the potential to enhance CKD care.

METHODS

We examined the association between 6592 plasma protein levels measured by aptamers and the severity of interstitial fibrosis and tubular atrophy (IFTA), glomerulosclerosis, arteriolar sclerosis, and arterial sclerosis among 434 participants of the Boston Kidney Biopsy Cohort. For proteins significantly associated with at least one histologic lesion, we assessed renal arteriovenous protein gradients among 21 individuals who had undergone invasive catheterization and assessed the expression of the cognate gene among 47 individuals with single cell RNA sequencing data in the Kidney Precision Medicine Project.

RESULTS

In models adjusted for estimated glomerular filtration rate (eGFR), proteinuria, and demographic factors, we identified 35 proteins associated with one or more chronic histologic lesions, including 20 specific for IFTA, 8 specific for glomerulosclerosis, and 1 specific for arteriolar sclerosis. In general, higher levels of these proteins were associated with more severe histologic score and lower eGFR. Exceptions included testican-2 and NELL1, which were associated with less glomerulosclerosis and IFTA, respectively, and higher eGFR; notably, both of these proteins demonstrated significantly higher levels from artery to renal vein, demonstrating net kidney release. In the Kidney Precision Medicine Project, 13 of the 35 protein hits had cognate gene expression enriched in one or more cell types in the kidney, including podocyte expression of select glomerulosclerosis markers (including testican-2) and tubular expression of several IFTA markers (including NELL1).

CONCLUSIONS

Proteomic analysis identified circulating proteins associated with chronic histopathologic lesions, some of which have concordant site-specific expression within the kidney.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Development and external validation of a diagnostic model for biopsy-proven acute interstitial nephritis using electronic health record data

BACKGROUND

Patients with acute interstitial nephritis (AIN) can present without typical clinical features, leading to a delay in diagnosis and treatment. We therefore developed and validated a diagnostic model to identify patients at risk of AIN using variables from the electronic health record.

METHODS

In patients who underwent a kidney biopsy at Yale University between 2013 and 2018, we tested the association of >150 variables with AIN, including demographics, comorbidities, vital signs and laboratory tests (training set 70%). We used least absolute shrinkage and selection operator methodology to select prebiopsy features associated with AIN. We performed area under the receiver operating characteristics curve (AUC) analysis with internal (held-out test set 30%) and external validation (Biopsy Biobank Cohort of Indiana). We tested the change in model performance after the addition of urine biomarkers in the Yale AIN study.

RESULTS

We included 393 patients (AIN 22%) in the training set, 158 patients (AIN 27%) in the test set, 1118 patients (AIN 11%) in the validation set and 265 patients (AIN 11%) in the Yale AIN study. Variables in the selected model included serum creatinine {adjusted odds ratio [aOR] 2.31 [95% confidence interval (CI) 1.42-3.76]}, blood urea nitrogen:creatinine ratio [aOR 0.40 (95% CI 0.20-0.78)] and urine dipstick specific gravity [aOR 0.95 (95% CI 0.91-0.99)] and protein [aOR 0.39 (95% CI 0.23-0.68)]. This model showed an AUC of 0.73 (95% CI 0.64-0.81) in the test set, which was similar to the AUC in the external validation cohort [0.74 (95% CI 0.69-0.79)]. The AUC improved to 0.84 (95% CI 0.76-0.91) upon the addition of urine interleukin-9 and tumor necrosis factor-α.

CONCLUSIONS

We developed and validated a statistical model that showed a modest AUC for AIN diagnosis, which improved upon the addition of urine biomarkers. Future studies could evaluate this model and biomarkers to identify unrecognized cases of AIN.

Design and rationale of GUARDD-US: A pragmatic, randomized trial of genetic testing for APOL1 and pharmacogenomic predictors of antihypertensive efficacy in patients with hypertension

RATIONALE AND OBJECTIVE

APOL1 risk alleles are associated with increased cardiovascular and chronic kidney disease (CKD) risk. It is unknown whether knowledge of APOL1 risk status motivates patients and providers to attain recommended blood pressure (BP) targets to reduce cardiovascular disease.

STUDY DESIGN

Multicenter, pragmatic, randomized controlled clinical trial.

SETTING AND PARTICIPANTS

6650 individuals with African ancestry and hypertension from 13 health systems.

INTERVENTION

APOL1 genotyping with clinical decision support (CDS) results are returned to participants and providers immediately (intervention) or at 6 months (control). A subset of participants are re-randomized to pharmacogenomic testing for relevant antihypertensive medications (pharmacogenomic sub-study). CDS alerts encourage appropriate CKD screening and antihypertensive agent use.

OUTCOMES

Blood pressure and surveys are assessed at baseline, 3 and 6 months. The primary outcome is change in systolic BP from enrollment to 3 months in individuals with two APOL1 risk alleles. Secondary outcomes include new diagnoses of CKD, systolic blood pressure at 6 months, diastolic BP, and survey results. The pharmacogenomic sub-study will evaluate the relationship of pharmacogenomic genotype and change in systolic BP between baseline and 3 months.

RESULTS

To date, the trial has enrolled 3423 participants.

CONCLUSIONS

The effect of patient and provider knowledge of APOL1 genotype on systolic blood pressure has not been well-studied. GUARDD-US addresses whether blood pressure improves when patients and providers have this information. GUARDD-US provides a CDS framework for primary care and specialty clinics to incorporate APOL1 genetic risk and pharmacogenomic prescribing in the electronic health record.

TRIAL REGISTRATION

ClinicalTrials.govNCT04191824.

Incidence and Importance of Calcium Deposition in Kidney Biopsy Specimens

INTRODUCTION

Calcification on native kidney biopsy specimens is often noted by pathologists, but the consequence is unknown.

METHODS

We searched the pathology reports in the Biopsy Biobank Cohort of Indiana for native biopsy specimens with calcification.

RESULTS

Of the 4,364 specimens, 416 (9.8%) had calcification. We compared clinical and histopathology findings in those with calcification (n = 429) compared to those without calcification (n = 3,936). Patients with calcification were older, had more comorbidities, lower estimated glomerular filtration rates (eGFR), were more likely to have hyaline arteriosclerosis, interstitial fibrosis/tubular atrophy, and a primary pathologic diagnosis of acute tubular injury or acute tubular necrosis when compared to patients without calcification. Patients with calcium oxalate deposition alone, compared to calcium phosphate or mixed calcifications, had fewer comorbidities but were more likely to have a history of gastric bypass surgery or malabsorption and take vitamin D. In patients with two or more years of follow-up, multivariate analyses showed the presence of calcification (HR 0.59, 0.38-0.92, p = 0.02) and higher eGFR (HR 0.76, 0.73-0.79, p < 0.001), was associated with decreased likelihood of progressing to end-stage renal disease. The presence of calcification was also associated with a reduced slope/decline in eGFR compared to known biopsy and clinical risk factors for decline in kidney function. We hypothesized this was due to more recoverable acute kidney injury (AKI) and found more severe acute kidney injury network stage in patients with kidney calcification but also greater improvement over time.

DISCUSSION/CONCLUSION

In summary, we demonstrated that calcification on kidney biopsy specimens was associated with a better prognosis than those without calcification due to the association with recoverable AKI.

A reference tissue atlas for the human kidney

Kidney Precision Medicine Project (KPMP) is building a spatially specified human kidney tissue atlas in health and disease with single-cell resolution. Here, we describe the construction of an integrated reference map of cells, pathways, and genes using unaffected regions of nephrectomy tissues and undiseased human biopsies from 56 adult subjects. We use single-cell/nucleus transcriptomics, subsegmental laser microdissection transcriptomics and proteomics, near-single-cell proteomics, 3D and CODEX imaging, and spatial metabolomics to hierarchically identify genes, pathways, and cells. Integrated data from these different technologies coherently identify cell types/subtypes within different nephron segments and the interstitium. These profiles describe cell-level functional organization of the kidney following its physiological functions and link cell subtypes to genes, proteins, metabolites, and pathways. They further show that messenger RNA levels along the nephron are congruent with the subsegmental physiological activity. This reference atlas provides a framework for the classification of kidney disease when multiple molecular mechanisms underlie convergent clinical phenotypes.

Spatial transcriptomics and the kidney

PURPOSE OF REVIEW

The application of spatial transcriptomics technologies to the interrogation of kidney tissue is a burgeoning effort. These technologies share a common purpose in mapping both the expression of individual molecules and entire transcriptomic signatures of kidney cell types and structures. Such information is often superimposed upon a histologic image. The resulting datasets are readily merged with other imaging and transcriptomic techniques to establish a spatially anchored atlas of the kidney. This review provides an overview of the various spatial transcriptomic technologies and recent studies in kidney disease. Potential applications gleaned from the interrogation of other organ systems, but relative to the kidney, are also discussed.

RECENT FINDINGS

Spatial transcriptomic technologies have enabled localization of whole transcriptome mRNA expression, correlation of mRNA to histology, measurement of in situ changes in expression across time, and even subcellular localization of transcripts within the kidney. These innovations continue to aid in the development of human cellular atlases of the kidney, the reclassification of disease, and the identification of important therapeutic targets.

SUMMARY

Spatial localization of gene expression will complement our current understanding of disease derived from single cell RNA sequencing, histopathology, protein immunofluorescence, and electron microscopy. Although spatial technologies continue to evolve rapidly, their importance in the localization of disease signatures is already apparent. Further efforts are required to integrate whole transcriptome and subcellular expression signatures into the individualized assessment of human kidney disease.

Molecular Signatures of Diabetic Kidney Disease Hiding in a Patient with Hypertension-Related Kidney Disease: A Clinical Pathologic Molecular Correlation

The Kidney Precision Medicine Project (KPMP) seeks to establish a molecular atlas of the kidney in health and disease and improve our understanding of the molecular drivers of CKD and AKI. Herein, we describe the case of a 66-year-old woman with CKD who underwent a protocol KPMP kidney biopsy. Her clinical history included well-controlled diabetes mellitus, hypertension, and proteinuria. The patient's histopathology was consistent with modest hypertension-related kidney injury, without overt diabetic kidney disease. Transcriptomic signatures of the glomerulus, interstitium, and tubular subsegments were obtained from laser microdissected tissue. The molecular signatures that were uncovered revealed evidence of early diabetic kidney disease adaptation and ongoing active tubular injury with enriched pathways related to mesangial cell hypertrophy, glycosaminoglycan biosynthesis, and apoptosis. Molecular evidence of diabetic kidney disease was found across the nephron. Novel molecular assays can supplement and enrich the histopathologic diagnosis obtained from a kidney biopsy.

Alterations in Protein Translation and Carboxylic Acid Catabolic Processes in Diabetic Kidney Disease

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.

Cellular and molecular interrogation of kidney biopsy specimens

PURPOSE OF REVIEW

Traditional histopathology of the kidney biopsy specimen has been an essential and successful tool for the diagnosis and staging of kidney diseases. However, it is likely that the full potential of the kidney biopsy has not been tapped so far. Indeed, there is now a concerted worldwide effort to interrogate kidney biopsy samples at the cellular and molecular levels with unprecedented rigor and depth. This review examines these novel approaches to study kidney biopsy specimens and highlights their potential to refine our understanding of the pathophysiology of kidney disease and lead to precision-based diagnosis and therapy.

RECENT FINDINGS

Several consortia are now active at studying kidney biopsy samples from various patient cohorts with state-of-the art cellular and molecular techniques. These include advanced imaging approaches as well as deep molecular interrogation with tools such as epigenetics, transcriptomics, proteomics and metabolomics. The emphasis throughout is on rigor, reproducibility and quality control.

SUMMARY

Although these techniques to study kidney biopsies are complementary, each on its own can yield novel ways to define and classify kidney disease. Therefore, great efforts are needed in order to generate an integrated output that can propel the diagnosis and treatment of kidney disease into the realm of precision medicine.

Label-free imaging of non-deparaffinized sections of the human kidney to determine tissue quality and signatures of disease

Label-free fluorescence imaging of kidney sections can provide important morphological information, but its utility has not been tested in a histology processing workflow. We tested the feasibility of label-free imaging of paraffin-embedded sections without deparaffinization and its potential usefulness in generating actionable data. Kidney tissue specimens were obtained during percutaneous nephrolithotomy or via diagnostic needle biopsy. Unstained non-deparaffinized sections were imaged using widefield fluorescence microscopy to capture endogenous fluorescence. Some samples were also imaged with confocal microscopy and multiphoton excitation to collect second harmonic generation (SHG) signal to obtain high-quality autofluorescence images with optical sectioning. To adjudicate the label-free signal, the samples or corresponding contiguous sections were subsequently deparaffinized and stained with Lillie's allochrome. Label-free imaging allowed the recognition of various kidney structures and enabled morphological qualification for adequacy. SHG and confocal imaging yielded quantifiable high-quality images for tissue collagens and revealed specific patterns in glomeruli and various tubules. Disease specimens from patients with diabetic kidney disease and focal segmental glomerulosclerosis showed distinctive signatures compared to specimens from healthy controls with normal kidney function. Quantitative cytometry could also be performed when DAPI is added in situ before imaging. These results show that label-free imaging of non-deparaffinized sections provides useful information about tissue quality that could be beneficial to nephropathologists by maximizing the use of scarce kidney tissue. This approach also provides quantifiable features that could inform on the biology of health and disease.

Deconvolution Tactics and Normalization in Renal Spatial Transcriptomics

The kidney is composed of heterogeneous groups of epithelial, endothelial, immune, and stromal cells, all in close anatomic proximity. Spatial transcriptomic technologies allow the interrogation of in situ expression signatures in health and disease, overlaid upon a histologic image. However, some spatial gene expression platforms have not yet reached single-cell resolution. As such, deconvolution of spatial transcriptomic spots is important to understand the proportion of cell signature arising from these varied cell types in each spot. This article reviews the various deconvolution strategies discussed in the 2021 Indiana O'Brien Center for Microscopy workshop. The unique features of Seurat transfer score methodology, SPOTlight, Robust Cell Type Decomposition, and BayesSpace are reviewed. The application of normalization and batch effect correction across spatial transcriptomic samples is also discussed.

Pharmacogenomics of Hypertension in CKD: The CKD-PGX Study

Background

Patients with CKD often have uncontrolled hypertension despite polypharmacy. Pharmacogenomic drug-gene interactions (DGIs) may affect the metabolism or efficacy of antihypertensive agents. We report changes in hypertension control after providing a panel of 11 pharmacogenomic predictors of antihypertensive response.

Methods

A prospective cohort with CKD and hypertension was followed to assess feasibility of pharmacogenomic testing implementation, self-reported provider utilization, and BP control. The analysis population included 382 subjects with hypertension who were genotyped for cross-sectional assessment of DGIs, and 335 subjects followed for 1 year to assess systolic BP (SBP) and diastolic BP (DBP).

Results

Most participants (58%) with uncontrolled hypertension had a DGI reducing the efficacy of one or more antihypertensive agents. Subjects with a DGI had 1.85-fold (95% CI, 1.2- to 2.8-fold) higher odds of uncontrolled hypertension, as compared with those without a DGI, adjusted for race, health system (safety-net hospital versus other locations), and advanced CKD (eGFR <30 ml/min). CYP2C9-reduced metabolism genotypes were associated with losartan response and uncontrolled hypertension (odds ratio [OR], 5.2; 95% CI, 1.9 to 14.7). CYP2D6-intermediate or -poor metabolizers had less frequent uncontrolled hypertension compared with normal metabolizers taking metoprolol or carvedilol (OR, 0.55; 95% CI, 0.3 to 0.95). In 335 subjects completing 1-year follow-up, SBP (-4.0 mm Hg; 95% CI, 1.6 to 6.5 mm Hg) and DBP (-3.3 mm Hg; 95% CI, 2.0 to 4.6 mm Hg) were improved. No significant difference in SBP or DBP change were found between individuals with and without a DGI.

Conclusions

There is a potential role for the addition of pharmacogenomic testing to optimize antihypertensive regimens in patients with CKD.

Integration of spatial and single-cell transcriptomics localizes epithelial cell-immune cross-talk in kidney injury

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.

Rationale and design of the Kidney Precision Medicine Project

Chronic kidney disease (CKD) and acute kidney injury (AKI) are common, heterogeneous, and morbid diseases. Mechanistic characterization of CKD and AKI in patients may facilitate a precision-medicine approach to prevention, diagnosis, and treatment. The Kidney Precision Medicine Project aims to ethically and safely obtain kidney biopsies from participants with CKD or AKI, create a reference kidney atlas, and characterize disease subgroups to stratify patients based on molecular features of disease, clinical characteristics, and associated outcomes. An additional aim is to identify critical cells, pathways, and targets for novel therapies and preventive strategies. This project is a multicenter prospective cohort study of adults with CKD or AKI who undergo a protocol kidney biopsy for research purposes. This investigation focuses on kidney diseases that are most prevalent and therefore substantially burden the public health, including CKD attributed to diabetes or hypertension and AKI attributed to ischemic and toxic injuries. Reference kidney tissues (for example, living-donor kidney biopsies) will also be evaluated. Traditional and digital pathology will be combined with transcriptomic, proteomic, and metabolomic analysis of the kidney tissue as well as deep clinical phenotyping for supervised and unsupervised subgroup analysis and systems biology analysis. Participants will be followed prospectively for 10 years to ascertain clinical outcomes. Cell types, locations, and functions will be characterized in health and disease in an open, searchable, online kidney tissue atlas. All data from the Kidney Precision Medicine Project will be made readily available for broad use by scientists, clinicians, and patients.

Multi-Institutional Implementation of Clinical Decision Support for APOL1, NAT2, and YEATS4 Genotyping in Antihypertensive Management

(1) Background: Clinical decision support (CDS) is a vitally important adjunct to the implementation of pharmacogenomic-guided prescribing in clinical practice. A novel CDS was sought for the APOL1, NAT2, and YEATS4 genes to guide optimal selection of antihypertensive medications among the African American population cared for at multiple participating institutions in a clinical trial. (2) Methods: The CDS committee, made up of clinical content and CDS experts, developed a framework and contributed to the creation of the CDS using the following guiding principles: 1. medical algorithm consensus; 2. actionability; 3. context-sensitive triggers; 4. workflow integration; 5. feasibility; 6. interpretability; 7. portability; and 8. discrete reporting of lab results. (3) Results: Utilizing the principle of discrete patient laboratory and vital information, a novel CDS for APOL1, NAT2, and YEATS4 was created for use in a multi-institutional trial based on a medical algorithm consensus. The alerts are actionable and easily interpretable, clearly displaying the purpose and recommendations with pertinent laboratory results, vitals and links to ordersets with suggested antihypertensive dosages. Alerts were either triggered immediately once a provider starts to order relevant antihypertensive agents or strategically placed in workflow-appropriate general CDS sections in the electronic health record (EHR). Detailed implementation instructions were shared across institutions to achieve maximum portability. (4) Conclusions: Using sound principles, the created genetic algorithms were applied across multiple institutions. The framework outlined in this study should apply to other disease-gene and pharmacogenomic projects employing CDS.

Large-scale, three-dimensional tissue cytometry of the human kidney: a complete and accessible pipeline

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.

Establishing the value of genomics in medicine: the IGNITE Pragmatic Trials Network

PURPOSE

A critical gap in the adoption of genomic medicine into medical practice is the need for the rigorous evaluation of the utility of genomic medicine interventions.

METHODS

The Implementing Genomics in Practice Pragmatic Trials Network (IGNITE PTN) was formed in 2018 to measure the clinical utility and cost-effectiveness of genomic medicine interventions, to assess approaches for real-world application of genomic medicine in diverse clinical settings, and to produce generalizable knowledge on clinical trials using genomic interventions. Five clinical sites and a coordinating center evaluated trial proposals and developed working groups to enable their implementation.

RESULTS

Two pragmatic clinical trials (PCTs) have been initiated, one evaluating genetic risk APOL1 variants in African Americans in the management of their hypertension, and the other to evaluate the use of pharmacogenetic testing for medications to manage acute and chronic pain as well as depression.

CONCLUSION

IGNITE PTN is a network that carries out PCTs in genomic medicine; it is focused on diversity and inclusion of underrepresented minority trial participants; it uses electronic health records and clinical decision support to deliver the interventions. IGNITE PTN will develop the evidence to support (or oppose) the adoption of genomic medicine interventions by patients, providers, and payers.

Clinical, histopathologic and molecular features of idiopathic and diabetic nodular mesangial sclerosis in humans

BACKGROUND

Idiopathic nodular mesangial sclerosis, also called idiopathic nodular glomerulosclerosis (ING), is a rare clinical entity with an unclear pathogenesis. The hallmark of this disease is the presence of nodular mesangial sclerosis on histology without clinical evidence of diabetes mellitus or other predisposing diagnoses. To achieve insights into its pathogenesis, we queried the clinical, histopathologic and transcriptomic features of ING and nodular diabetic nephropathy (DN).

METHODS

All renal biopsy reports accessioned at Indiana University Health from 2001 to 2016 were reviewed to identify 48 ING cases. Clinical and histopathologic features were compared between individuals with ING and DN (n = 751). Glomeruli of ING (n = 5), DN (n = 18) and reference (REF) nephrectomy (n = 9) samples were isolated by laser microdissection and RNA was sequenced. Immunohistochemistry of proline-rich 36 (PRR36) protein was performed.

RESULTS

ING subjects were frequently hypertensive (95.8%) with a smoking history (66.7%). ING subjects were older, had lower proteinuria and had less hyaline arteriolosclerosis than DN subjects. Butanoate metabolism was an enriched pathway in ING samples compared with either REF or DN samples. The top differentially expressed gene, PRR36, had increased expression in glomeruli 248-fold [false discovery rate (FDR) P = 5.93 × 10-6] compared with the REF and increased 109-fold (FDR P = 1.85 × 10-6) compared with DN samples. Immunohistochemistry revealed a reduced proportion of cells with perinuclear reaction in ING samples as compared to DN.

CONCLUSIONS

Despite similar clinical and histopathologic characteristics in ING and DN, the uncovered transcriptomic signature suggests that ING has distinct molecular features from nodular DN. Further study is warranted to understand these relationships.

Molecular characterization of the human kidney interstitium in health and disease

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 orchestrated cellular and molecular responses of the kidney to endotoxin define a precise sepsis timeline

Sepsis is a dynamic state that progresses at variable rates and has life-threatening consequences. Staging patients along the sepsis timeline requires a thorough knowledge of the evolution of cellular and molecular events at the tissue level. Here, we investigated the kidney, an organ central to the pathophysiology of sepsis. Single-cell RNA-sequencing in a murine endotoxemia model revealed the involvement of various cell populations to be temporally organized and highly orchestrated. Endothelial and stromal cells were the first responders. At later time points, epithelial cells upregulated immune-related pathways while concomitantly downregulating physiological functions such as solute homeostasis. Sixteen hours after endotoxin, there was global cell–cell communication failure and organ shutdown. Despite this apparent organ paralysis, upstream regulatory analysis showed significant activity in pathways involved in healing and recovery. This rigorous spatial and temporal definition of murine endotoxemia will uncover precise biomarkers and targets that can help stage and treat human sepsis.

Multiple reaction monitoring profiling as an analytical strategy to investigate lipids in extracellular vesicles

Extracellular vesicles (EVs) convey information used in cell-to-cell interactions. Lipid analysis of EVs remains challenging because of small sample amounts available. Lipid discovery using traditional mass spectrometry platforms based on liquid chromatography and high mass resolution typically employs milligram sample amounts. We report a simple workflow for lipid profiling of EVs based on multiple reaction monitoring (MRM) profiling that uses microgram amounts of sample. After liquid-liquid extraction, individual EV samples were injected directly into the electrospray ionization (ESI) ion source at low flow rates (10 μl/min) and screened for 197 MRM transitions chosen to be a characteristic of several classes of lipids. This choice was based on a discovery experiment, which applied 1,419 MRMs associated with multiple lipid classes to a representative pooled sample. EVs isolated from 12 samples of human lymphocytes and 16 replicates from six different rat cells lines contained an estimated amount of total lipids of 326 to 805 μg. Samples showed profiles that included phosphatidylcholine (PC), sphingomyelin (SM), cholesteryl ester (CE), and ceramide (Cer) lipids, as well as acylcarnitines. The lipid profiles of human lymphocyte EVs were distinguishable using principal component and cluster analysis in terms of prior antibody and drug exposure. Lipid profiles of rat cell lines EV's were distinguishable by their tissue of origin.

A multimodal and integrated approach to interrogate human kidney biopsies with rigor and reproducibility: guidelines from the Kidney Precision Medicine Project

Comprehensive and spatially mapped molecular atlases of organs at a cellular level are a critical resource to gain insights into pathogenic mechanisms and personalized therapies for diseases. The Kidney Precision Medicine Project (KPMP) is an endeavor to generate three-dimensional (3-D) molecular atlases of healthy and diseased kidney biopsies by using multiple state-of-the-art omics and imaging technologies across several institutions. Obtaining rigorous and reproducible results from disparate methods and at different sites to interrogate biomolecules at a single-cell level or in 3-D space is a significant challenge that can be a futile exercise if not well controlled. We describe a "follow the tissue" pipeline for generating a reliable and authentic single-cell/region 3-D molecular atlas of human adult kidney. Our approach emphasizes quality assurance, quality control, validation, and harmonization across different omics and imaging technologies from sample procurement, processing, storage, shipping to data generation, analysis, and sharing. We established benchmarks for quality control, rigor, reproducibility, and feasibility across multiple technologies through a pilot experiment using common source tissue that was processed and analyzed at different institutions and different technologies. A peer review system was established to critically review quality control measures and the reproducibility of data generated by each technology before their being approved to interrogate clinical biopsy specimens. The process established economizes the use of valuable biopsy tissue for multiomics and imaging analysis with stringent quality control to ensure rigor and reproducibility of results and serves as a model for precision medicine projects across laboratories, institutions and consortia.

In Situ Classification of Cell Types in Human Kidney Tissue Using 3D Nuclear Staining

To understand the physiology and pathology of disease, capturing the heterogeneity of cell types within their tissue environment is fundamental. In such an endeavor, the human kidney presents a formidable challenge because its complex organizational structure is tightly linked to key physiological functions. Advances in imaging-based cell classification may be limited by the need to incorporate specific markers that can link classification to function. Multiplex imaging can mitigate these limitations, but requires cumulative incorporation of markers, which may lead to tissue exhaustion. Furthermore, the application of such strategies in large scale 3-dimensional (3D) imaging is challenging. Here, we propose that 3D nuclear signatures from a DNA stain, DAPI, which could be incorporated in most experimental imaging, can be used for classifying cells in intact human kidney tissue. We developed an unsupervised approach that uses 3D tissue cytometry to generate a large training dataset of nuclei images (NephNuc), where each nucleus is associated with a cell type label. We then devised various supervised machine learning approaches for kidney cell classification and demonstrated that a deep learning approach outperforms classical machine learning or shape-based classifiers. Specifically, a custom 3D convolutional neural network (NephNet3D) trained on nuclei image volumes achieved a balanced accuracy of 80.26%. Importantly, integrating NephNet3D classification with tissue cytometry allowed in situ visualization of cell type classifications in kidney tissue. In conclusion, we present a tissue cytometry and deep learning approach for in situ classification of cell types in human kidney tissue using only a DNA stain. This methodology is generalizable to other tissues and has potential advantages on tissue economy and non-exhaustive classification of different cell types.

Risk Factors for Bleeding and Clinical Ineffectiveness Associated With Clopidogrel Therapy: A Comprehensive Meta-Analysis

Although clopidogrel is a frequently used antiplatelet medication to treat and prevent atherothrombotic disease, clinicians must balance its clinical effectiveness with the potential side effect of bleeding. However, many previous studies have evaluated beneficial and adverse factors separately. The objective of our study was to perform a comprehensive meta‐analysis of studies of clopidogrel’s clinical effectiveness and/or risk of bleeding in order to identify and assess all reported risk factors, thus helping clinicians to balance patient safety with drug efficacy. We analyzed randomized controlled trials (RCTs) of maintenance use in four stages: search for relevant primary articles; abstract and full article screening; quality assessment and data extraction; and synthesis and data analysis. Screening of 7,109 articles yielded 52 RCTs that met the inclusion criteria. Twenty‐seven risk factors were identified. “Definite risk factors” were defined as those with aggregated odds ratios (ORs) > 1 and confidence intervals (CIs) > 1 if analyzed in more than one study. Definite risk factors for major bleeding were concomitant aspirin use (OR 2.83, 95% CI 2.04–3.94) and long duration of clopidogrel therapy (> 6 months) (OR 1.74, 95% CI 1.21–2.50). Dual antiplatelet therapy, extended clopidogrel therapy, and high maintenance dose (150 mg/day) of clopidogrel were definite risk factors for any bleeding. Reduced renal function, both mild and severe, was the only definite risk factor for clinical ineffectiveness. These findings can help clinicians predict the risks and effectiveness of clopidogrel use for their patients and be used in clinical decision support tools.

Targeting fibroblast growth factor 23-responsive pathways uncovers controlling genes in kidney mineral metabolism

Fibroblast Growth Factor 23 (FGF23) is a bone-derived hormone that reduces kidney phosphate reabsorption and 1,25(OH)2 vitamin D synthesis via its required co-receptor alpha-Klotho. To identify novel genes that could serve as targets to control FGF23-mediated mineral metabolism, gene array and single-cell RNA sequencing were performed in wild type mouse kidneys. Gene array demonstrated that heparin-binding EGF-like growth factor (HBEGF) was significantly up-regulated following one-hour FGF23 treatment of wild type mice. Mice injected with HBEGF had phenotypes consistent with partial FGF23-mimetic activity including robust induction of Egr1, and increased Cyp24a1 mRNAs. Single cell RNA sequencing showed overlapping HBEGF and EGF-receptor expression mostly in the proximal tubule, and alpha-Klotho expression in proximal and distal tubule segments. In alpha-Klotho-null mice devoid of canonical FGF23 signaling, HBEGF injections significantly increased Egr1 and Cyp24a1 with correction of basally elevated Cyp27b1. Additionally, mice placed on a phosphate deficient diet to suppress FGF23 had endogenously increased Cyp27b1 mRNA, which was rescued in mice receiving HBEGF. In HEK293 cells with stable alpha-Klotho expression, FGF23 and HBEGF increased CYP24A1 mRNA expression. HBEGF, but not FGF23 bioactivity was blocked with EGF-receptor inhibition. Thus, our findings support that the paracrine/autocrine factor HBEGF could play novel roles in controlling genes downstream of FGF23 via targeting common signaling pathways.

Genotype-Guided Hydralazine Therapy

BACKGROUND

Despite its approval in 1953, hydralazine hydrochloride continues to be used in the management of resistant hypertension, a condition frequently managed by nephrologists and other clinicians. Hydralazine hydrochloride undergoes metabolism by the N-acetyltransferase 2 (NAT2) enzyme. NAT2 is highly polymorphic as approximately 50% of the general population are slow acetylators. In this review, we first evaluate the link between NAT2 genotype and phenotype. We then assess the evidence available for genotype-guided therapy of hydralazine, specifically addressing associations of NAT2 acetylator status with hydralazine pharmacokinetics, antihypertensive efficacy, and toxicity.

SUMMARY

There is a critical need to use hydralazine in some patients with resistant hypertension. Available evidence supports a significant link between genotype and NAT2 enzyme activity as 29 studies were identified with an overall concordance between genotype and phenotype of 92%. The literature also supports an association between acetylator status and hydralazine concentration, as fourteen of fifteen identified studies revealed significant relationships with a consistent direction of effect. Although fewer studies are available to directly link acetylator status with hydralazine antihypertensive efficacy, the evidence from this smaller set of studies is significant in 7 of 9 studies identified. Finally, 5 studies were identified which support the association of acetylator status with hydralazine-induced lupus. Clinicians should maintain vigilance when prescribing maximum doses of hydralazine. Key Messages: NAT2 slow acetylator status predicts increased hydralazine levels, which may lead to increased efficacy and adverse effects. Caution should be exercised in slow acetylators with total daily hydralazine doses of 200 mg or more. Fast acetylators are at risk for inefficacy at lower doses of hydralazine. With appropriate guidance on the usage of NAT2 genotype, clinicians can adopt a personalized approach to hydralazine dosing and prescription, enabling more efficient and safe treatment of resistant hypertension.

Application of Laser Microdissection to Uncover Regional Transcriptomics in Human Kidney Tissue

Gene expression analysis of human kidney tissue is an important tool to understand homeostasis and disease pathophysiology. Increasing the resolution and depth of this technology and extending it to the level of cells within the tissue is needed. Although the use of single nuclear and single cell RNA sequencing has become widespread, the expression signatures of cells obtained from tissue dissociation do not maintain spatial context. Laser microdissection (LMD) based on specific fluorescent markers would allow the isolation of specific structures and cell groups of interest with known localization, thereby enabling the acquisition of spatially-anchored transcriptomic signatures in kidney tissue. We have optimized an LMD methodology, guided by a rapid fluorescence-based stain, to isolate five distinct compartments within the human kidney and conduct subsequent RNA sequencing from valuable human kidney tissue specimens. We also present quality control parameters to enable the assessment of adequacy of the collected specimens. The workflow outlined in this manuscript shows the feasibility of this approach to isolate sub-segmental transcriptomic signatures with high confidence. The methodological approach presented here may also be applied to other tissue types with substitution of relevant antibody markers.

Interindividual Variability in Lymphocyte Stimulation and Transcriptomic Response Predicts Mycophenolic Acid Sensitivity in Healthy Volunteers

Mycophenolic acid (MPA) is an immunosuppressant commonly used to prevent renal transplant rejection and treat glomerulonephritis. MPA inhibits IMPDH2 within stimulated lymphocytes, reducing guanosine synthesis. Despite the widespread use of MPA, interindividual variability in response remains with rates of allograft rejection up to 15% and approximately half of individuals fail to achieve complete remission to lupus nephritis. We sought to identify contributors to interindividual variability in MPA response, hypothesizing that the HPRT1 salvage guanosine synthesis contributes to variability. MPA sensitivity was measured in 40 healthy individuals using an ex vivo lymphocyte viability assay. Measurement of candidate gene expression (n ± 40) and single‐cell RNA‐sequencing (n ± 6) in lymphocytes was performed at baseline, poststimulation, and post‐MPA treatment. After stimulation, HPRT1 expression was 2.1‐fold higher in resistant individuals compared with sensitive individuals (P ± 0.049). Knockdown of HPRT1 increased MPA sensitivity (12%; P ± 0.003), consistent with higher expression levels in resistant individuals. Sensitive individuals had higher IMPDH2 expression and 132% greater stimulation. In lymphocyte subpopulations, differentially expressed genes between sensitive and resistant individuals included KLF2 and LTB. Knockdown of KLF2 and LTB aligned with the predicted direction of effect on proliferation. In sensitive individuals, more frequent receptor‐ligand interactions were observed after stimulation (P ± 0.0004), but fewer interactions remained after MPA treatment (P ± 0.0014). These data identify a polygenic transcriptomic signature in lymphocyte subpopulations predictive of MPA response. The degree of lymphocyte stimulation, HPRT1, KLF2, and LTB expression may serve as markers of MPA efficacy.

Kidney Histopathology and Prediction of Kidney Failure: A Retrospective Cohort Study

RATIONALE & OBJECTIVE

The use of kidney histopathology for predicting kidney failure is not established. We hypothesized that the use of histopathologic features of kidney biopsy specimens would improve prediction of clinical outcomes made using demographic and clinical variables alone.

STUDY DESIGN

Retrospective cohort study and development of a clinical prediction model.

SETTING & PARTICIPANTS

All 2,720 individuals from the Biopsy Biobank Cohort of Indiana who underwent kidney biopsy between 2002 and 2015 and had at least 2 years of follow-up.

NEW PREDICTORS & ESTABLISHED PREDICTORS

Demographic variables, comorbid conditions, baseline clinical characteristics, and histopathologic features.

OUTCOMES

Time to kidney failure, defined as sustained estimated glomerular filtration rate ≤ 10mL/min/1.73m².

ANALYTICAL APPROACH

Multivariable Cox regression model with internal validation by bootstrapping. Models including clinical and demographic variables were fit with the addition of histopathologic features. To assess the impact of adding a histopathology variable, the amount of variance explained (r²) and the C index were calculated. The impact on prediction was assessed by calculating the net reclassification index for each histopathologic variable and for all combined.

RESULTS

Median follow-up was 3.1 years. Within 5 years of biopsy, 411 (15.1%) patients developed kidney failure. Multivariable analyses including demographic and clinical variables revealed that severe glomerular obsolescence (adjusted HR, 2.03; 95% CI, 1.51-2.03), severe interstitial fibrosis and tubular atrophy (adjusted HR, 1.99; 95% CI, 1.52-2.59), and severe arteriolar hyalinosis (adjusted HR, 1.53; 95% CI, 1.14-2.05) were independently associated with the primary outcome. The addition of all histopathologic variables to the clinical model yielded a net reclassification index for kidney failure of 5.1% (P < 0.001) with a full model C statistic of 0.915. Analyses addressing the competing risk for death, optimism, or shrinkage did not significantly change the results.

LIMITATIONS

Selection bias from the use of clinically indicated biopsies and exclusion of patients with less than 2 years of follow-up, as well as reliance on surrogate indicators of kidney failure onset.

CONCLUSIONS

A model incorporating histopathologic features from kidney biopsy specimens improved prediction of kidney failure and may be valuable clinically. Future studies will be needed to understand whether even more detailed characterization of kidney tissue may further improve prognostication about the future trajectory of estimated glomerular filtration rate.

Implementation of a Renal Precision Medicine Program: Clinician Attitudes and Acceptance

A precision health initiative was implemented across a multi-hospital health system, wherein a panel of genetic variants was tested and utilized in the clinical care of chronic kidney disease (CKD) patients. Pharmacogenomic predictors of antihypertensive response and genomic predictors of CKD were provided to clinicians caring for nephrology patients. To assess clinician knowledge, attitudes, and willingness to act on genetic testing results, a Likert-scale survey was sent to and self-administered by these nephrology providers (N = 76). Most respondents agreed that utilizing pharmacogenomic-guided antihypertensive prescribing is valuable (4.0 ± 0.7 on a scale of 1 to 5, where 5 indicates strong agreement). However, the respondents also expressed reluctance to use genetic testing for CKD risk stratification due to a perceived lack of supporting evidence (3.2 ± 0.9). Exploratory sub-group analyses associated this reluctance with negative responses to both knowledge and attitude discipline questions, thus suggesting reduced exposure to and comfort with genetic information. Given the evolving nature of genomic implementation in clinical care, further education is warranted to help overcome these perception barriers.

Nephrotoxicity in a Patient With Inadequate Pain Control: Potential Role of Pharmacogenetic Testing for Cytochrome P450 2D6 and Apolipoprotein L1

A case is presented which demonstrates the perils of opioid inefficacy and how pharmacogenomic testing may have prevented nonsteroidal anti-inflammatory drug (NSAID)-induced nephrotoxicity and progression to chronic kidney disease (CKD). A 62 year-old female with back pain was treated with tramadol and hydrocodone; however, neither proved effective. Consequently, to control her pain, she resorted to cocaine, marijuana, and high dose nonsteroidal anti-inflammatory drugs (NSAIDs). She eventually developed CKD. To identify CKD contributors, she underwent genotyping for Apolipoprotein L1 (APOL1), a known risk factor of CKD, as well as relevant pharmacogenomic genes. Her APOL1 genotype was *G1(GM)/*G1(GM), placing her at increased risk of CKD progression. Her CYP2D6 genotype was *5/*17, consistent with intermediate metabolism, making opioid drugs reliant on CYP2D6 activation, such as tramadol and hydrocodone, relatively ineffective in this patient. Thus, this patient was at genetic risk for CKD and reduced opioid efficacy. We conclude that this genetic combination likely contributed to opioid inefficacy and the eventual progression to CKD.

RegSNPs-intron: a computational framework for predicting pathogenic impact of intronic single nucleotide variants

Single nucleotide variants (SNVs) in intronic regions have yet to be systematically investigated for their disease-causing potential. Using known pathogenic and neutral intronic SNVs (iSNVs) as training data, we develop the RegSNPs-intron algorithm based on a random forest classifier that integrates RNA splicing, protein structure, and evolutionary conservation features. RegSNPs-intron showed excellent performance in evaluating the pathogenic impacts of iSNVs. Using a high-throughput functional reporter assay called ASSET-seq (ASsay for Splicing using ExonTrap and sequencing), we evaluate the impact of RegSNPs-intron predictions on splicing outcome. Together, RegSNPs-intron and ASSET-seq enable effective prioritization of iSNVs for disease pathogenesis.

Increasing Nephrologist Awareness of Symptom Burden in Older Hospitalized End-Stage Renal Disease Patients

BACKGROUND

End-stage renal disease (ESRD) patients have significant symptom burden. Reduced provider awareness of symptoms contributes to underutilization of symptom management resources.

METHOD

We hypothesized that improved nephrologist awareness of symptoms leads to symptom improvement. In this prospective, multicenter interventional study, 53 (age >65) ESRD inpatients underwent symptom assessment using the modified Edmonton Symptom Assessment System (ESAS) at admission and 1-week post-discharge. Physicians caring for the enrollees were asked if they felt their patients would die within the year, and then sequentially randomized to receive the results of the baseline survey (group 1) or to not receive the results (group 2).

RESULTS

Fifty-two patients completed the study; 1 died. Baseline characteristics were compared. For 70% of the total cohort, physicians reported that they would not be surprised if their patient died within a year. There was no difference in baseline scores of the patients between the 2 physician groups. Severity ratings were compared between in-hospital and post discharge scores and between physicians who received the results versus those that did not. Total ESAS scores improved more in group 1 (12.9) than in group 2 (9.2; p = 0.04). Among individual symptoms, there was greater improvement in pain control (p = 0.02), and nominal improvement in itching (p = 0.03) in group 1 as compared to group 2. There were 3 palliative care consults.

CONCLUSIONS

Our findings reinforce the high symptom burden prevalent in older ESRD patients. The improvement in total scores, and individual symptoms of pain and itching in group 1 indicates better symptom control when physician awareness is increased. Residual symptoms post hospitalization and low utilization of palliative care resources are suggestive of a missed opportunity by nephrologists to address the high symptom burden at the inpatient encounter, which is selective for sick patients and/or indication of inadequacy of dialysis to control these symptoms.