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.

C-peptide and metabolic outcomes in trials of disease modifying therapy in new-onset type 1 diabetes: an individual participant meta-analysis.

BACKGROUND: Metabolic outcomes in type 1 diabetes remain suboptimal. Disease modifying therapy to prevent ß-cell loss presents an alternative treatment framework but the effect on metabolic outcomes is unclear. We, therefore, aimed to define the relationship between insulin C-peptide as a marker of ß-cell function and metabolic outcomes in new-onset type 1 diabetes. METHODS: 21 trials of disease-modifying interventions within 100 days of type 1 diabetes diagnosis comprising 1315 adults (ie, those 18 years and older) and 1396 children (ie, those younger than 18 years) were combined. Endpoints assessed were stimulated area under the curve C-peptide, HbA1c, insulin use, hypoglycaemic events, and composite scores (such as insulin dose adjusted A1c, total daily insulin, U/kg per day, and BETA-2 score). Positive studies were defined as those meeting their primary endpoint. Differences in outcomes between active and control groups were assessed using the Wilcoxon rank test. FINDINGS: 6 months after treatment, a 24·8% greater C-peptide preservation in positive studies was associated with a 0·55% lower HbA1c (p<0·0001), with differences being detectable as early as 3 months. Cross-sectional analysis, combining positive and negative studies, was consistent with this proportionality: a 55% improvement in C-peptide preservation was associated with 0·64% lower HbA1c (p<0·0001). Higher initial C-peptide levels and greater preservation were associated with greater improvement in HbA1c. For HbA1c, IDAAC, and BETA-2 score, sample size predictions indicated that 2-3 times as many participants per group would be required to show a difference at 6 months as compared with C-peptide. Detecting a reduction in hypoglycaemia was affected by reporting methods. INTERPRETATION: Interventions that preserve ß-cell function are effective at improving metabolic outcomes in new-onset type 1 diabetes, confirming their potential as adjuncts to insulin. We have shown that improvements in HbA1c are directly proportional to the degree of C-peptide preservation, quantifying this relationship, and supporting the use of C-peptides as a surrogate endpoint in clinical trials. FUNDING: JDRF and Diabetes UK.

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.

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.

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.

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.

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.

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.

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.

Influence of Uridine Diphosphate Glucuronosyltransferase Family 1 Member A1 and Solute Carrier Organic Anion Transporter Family 1 Member B1 Polymorphisms and Efavirenz on Bilirubin Disposition in Healthy Volunteers.

Chronic administration of efavirenz is associated with decreased serum bilirubin levels, probably through induction of UGT1A1 We assessed the impact of efavirenz monotherapy and UGT1A1 phenotypes on total, conjugated, and unconjugated serum bilirubin levels in healthy volunteers. Healthy volunteers were enrolled into a clinical study designed to address efavirenz pharmacokinetics, drug interactions, and pharmacogenetics. Volunteers received multiple oral doses (600 mg/day for 17 days) of efavirenz. Serum bilirubin levels were obtained at study entry and 1 week after completion of the study. DNA genotyping was performed for UGT1A1 [*80 (C>T), *6 (G>A), *28 (TA7), *36 (TA5), and *37 (TA8)] and for SLCO1B1 [*5 (521T>C) and *1b (388A>G] variants. Diplotype predicted phenotypes were classified as normal, intermediate, and slow metabolizers. Compared with bilirubin levels at screening, treatment with efavirenz significantly reduced total, conjugated, and unconjugated bilirubin. After stratification by UGT1A1 phenotypes, there was a significant decrease in total bilirubin among all phenotypes, conjugated bilirubin among intermediate metabolizers, and unconjugated bilirubin among normal and intermediate metabolizers. The data also show that UGT1A1 genotype predicts serum bilirubin levels at baseline, but this relationship is lost after efavirenz treatment. SLCO1B1 genotypes did not predict bilirubin levels at baseline or after efavirenz treatment. Our data suggest that efavirenz may alter bilirubin disposition mainly through induction of UGT1A1 metabolism and efflux through multidrug resistance-associated protein 2. SIGNIFICANCE STATEMENT: Efavirenz likely alters the pharmacokinetics of coadministered drugs, potentially causing lack of efficacy or increased adverse effects, as well as the disposition of endogenous compounds relevant in homeostasis through upregulation of UGT1A1 and multidrug resistance-associated protein 2. Measurement of unconjugated and conjugated bilirubin during new drug development may provide mechanistic understanding regarding enzyme and transporters modulated by the new drug.

Analytical validity of a genotyping assay for use with personalized antihypertensive and chronic kidney disease therapy.

Hypertension and chronic kidney disease are inextricably linked. Hypertension is a well-recognized contributor to chronic kidney disease progression and, in turn, renal disease potentiates hypertension. A generalized approach to drug selection and dosage has not proven effective in managing these conditions, in part, because patients with heterogeneous kidney disease and hypertension etiologies are frequently grouped according to functional or severity classifications. Genetic testing may serve as an important tool in the armamentarium of clinicians who embrace precision medicine. Increasing scientific evidence has supported the utilization of genomic information to select efficacious antihypertensive therapy and understand hereditary contributors to chronic kidney disease progression. Given the wide array of antihypertensive agents available and diversity of genetic renal disease predictors, a panel-based approach to genotyping may be an efficient and economic means of establishing an individualized blood pressure response profile for patients with various forms of chronic kidney disease and hypertension. In this manuscript, we discuss the validation process of a Clinical Laboratory Improvement Amendments-approved genetic test to relay information on 72 genetic variants associated with kidney disease progression and hypertension therapy. These genomic-based interventions, in addition to routine clinical data, may help inform physicians to provide personalized therapy.

Improving the Patency of Jugular Vein Catheters in Sprague-Dawley Rats by Using an Antiseptic Nitrocellulose Coating.

Preclinical studies in animals often require frequent blood sampling over prolonged periods. A preferred method in rats is the implantation of a polyurethane catheter into the jugular vein, with heparinized glycerol as a lock solution. However, analysis of various biologic compounds (for example, microRNA) precludes the use of heparin. We used sodium citrate as an alternative to heparin but observed more frequent loss of catheter patency. We hypothesized that this effect was due to evaporation of lock solution at the exteriorized portion of the catheter, subsequent blood infiltration into the catheter, and ultimately clot formation within the catheter. We therefore tested evaporation and its variables in vitro by using 5 common catheter materials. We used the migration of dye into vertically anchored catheters as a measure of lock displacement due to evaporation. Exposure to dry room-temperature air was sufficient to cause dye migration against gravity, whereas a humid environment and adding glycerol to the lock solution mitigated this effect, thus confirming loss of the lock solution from the catheter by evaporation. We tested 4 catheter treatments for the ability to reduce lock evaporation. Results were validated in vivo by using male Sprague-Dawley rats (n = 12) implanted with polyurethane jugular vein catheters and randomized to receive a nitrocellulose-based coating on the exteriorized portion of the catheter. Coating the catheters significantly improved patency, as indicated by a Kaplan-Meier log-rank hazard ratio greater than 5 in untreated catheters. We here demonstrate that a simple nitrocellulose coating reduces evaporation from and thus prolongs the patency of polyurethane catheters in rats.

PASSPORT-seq: A Novel High-Throughput Bioassay to Functionally Test Polymorphisms in Micro-RNA Target Sites.

Next-generation sequencing (NGS) studies have identified large numbers of genetic variants that are predicted to alter miRNA-mRNA interactions. We developed a novel high-throughput bioassay, PASSPORT-seq, that can functionally test in parallel 100s of these variants in miRNA binding sites (mirSNPs). The results are highly reproducible across both technical and biological replicates. The utility of the bioassay was demonstrated by testing 100 mirSNPs in HEK293, HepG2, and HeLa cells. The results of several of the variants were validated in all three cell lines using traditional individual luciferase assays. Fifty-five mirSNPs were functional in at least one of three cell lines (FDR ≤ 0.05); 11, 36, and 27 of them were functional in HEK293, HepG2, and HeLa cells, respectively. Only four of the variants were functional in all three cell lines, which demonstrates the cell-type specific effects of mirSNPs and the importance of testing the mirSNPs in multiple cell lines. Using PASSPORT-seq, we functionally tested 111 variants in the 3' UTR of 17 pharmacogenes that are predicted to alter miRNA regulation. Thirty-three of the variants tested were functional in at least one cell line.

In Vivo siRNA Delivery and Rebound of Renal LRP2 in Mice.

siRNA stabilized for in vivo applications is filtered and reabsorbed in the renal proximal tubule (PT), reducing mRNA expression transiently. Prior siRNA efforts have successfully prevented upregulation of mRNA in response to injury. We proposed reducing constitutive gene and protein expression of LRP2 (megalin) in order to understand its molecular regulation in mice. Using siRNA targeting mouse LRP2 (siLRP2), reduction of LRP2 mRNA expression was compared to scrambled siRNA (siSCR) in mouse PT cells. Mice received siLRP2 administration optimized for dose, administration site, carrier solution, administration frequency, and administration duration. Kidney cortex was collected upon sacrifice. Renal gene and protein expression were compared by qRT-PCR, immunoblot, and immunohistochemistry (IHC). Compared to siSCR, siLRP2 reduced mRNA expression in PT cells to 16.6% ± 0.6%. In mouse kidney cortex, siLRP2 reduced mRNA expression to 74.8 ± 6.3% 3 h and 70.1 ± 6.3% 6 h after administration. mRNA expression rebounded at 12 h (160.6 ± 11.2%). No megalin renal protein expression reduction was observed by immunoblot or IHC, even after serial twice daily dosing for 3.5 days. Megalin is a constitutively expressed protein. Although LRP2 renal mRNA expression reduction was achieved, siRNA remains a costly and inefficient intervention to reduce in vivo megalin protein expression.