Single-cell analysis identifies the interaction of altered renal tubules with basophils orchestrating kidney fibrosis.

Inflammation is an important component of fibrosis but immune processes that orchestrate kidney fibrosis are not well understood. Here we apply single-cell sequencing to a mouse model of kidney fibrosis. We identify a subset of kidney tubule cells with a profibrotic-inflammatory phenotype characterized by the expression of cytokines and chemokines associated with immune cell recruitment. Receptor-ligand interaction analysis and experimental validation indicate that CXCL1 secreted by profibrotic tubules recruits CXCR2+ basophils. In mice, these basophils are an important source of interleukin-6 and recruitment of the TH17 subset of helper T cells. Genetic deletion or antibody-based depletion of basophils results in reduced renal fibrosis. Human kidney single-cell, bulk gene expression and immunostaining validate a function for basophils in patients with kidney fibrosis. Collectively, these studies identify basophils as contributors to the development of renal fibrosis and suggest that targeting these cells might be a useful clinical strategy to manage chronic kidney disease.

APOL1 risk variants in individuals of African genetic ancestry drive endothelial cell defects that exacerbate sepsis.

The incidence and severity of sepsis is higher among individuals of African versus European ancestry. We found that genetic risk variants (RVs) in the trypanolytic factor apolipoprotein L1 (APOL1), present only in individuals of African ancestry, were associated with increased sepsis incidence and severity. Serum APOL1 levels correlated with sepsis and COVID-19 severity, and single-cell sequencing in human kidneys revealed high expression of APOL1 in endothelial cells. Analysis of mice with endothelial-specific expression of RV APOL1 and in vitro studies demonstrated that RV APOL1 interfered with mitophagy, leading to cytosolic release of mitochondrial DNA and activation of the inflammasome (NLRP3) and the cytosolic nucleotide sensing pathways (STING). Genetic deletion or pharmacological inhibition of NLRP3 and STING protected mice from RV APOL1-induced permeability defects and proinflammatory endothelial changes in sepsis. Our studies identify the inflammasome and STING pathways as potential targets to reduce APOL1-associated health disparities in sepsis and COVID-19.

Stretch regulates alveologenesis and homeostasis via mesenchymal Gαq/11-mediated TGFß2 activation.

Alveolar development and repair require tight spatiotemporal regulation of numerous signalling pathways that are influenced by chemical and mechanical stimuli. Mesenchymal cells play key roles in numerous developmental processes. Transforming growth factor-ß (TGFß) is essential for alveologenesis and lung repair, and the G protein α subunits Gαq and Gα11 (Gαq/11) transmit mechanical and chemical signals to activate TGFß in epithelial cells. To understand the role of mesenchymal Gαq/11 in lung development, we generated constitutive (Pdgfrb-Cre+/-;Gnaqfl/fl;Gna11-/-) and inducible (Pdgfrb-Cre/ERT2+/-;Gnaqfl/fl;Gna11-/-) mesenchymal Gαq/11 deleted mice. Mice with constitutive Gαq/11 gene deletion exhibited abnormal alveolar development, with suppressed myofibroblast differentiation, altered mesenchymal cell synthetic function, and reduced lung TGFß2 deposition, as well as kidney abnormalities. Tamoxifen-induced mesenchymal Gαq/11 gene deletion in adult mice resulted in emphysema associated with reduced TGFß2 and elastin deposition. Cyclical mechanical stretch-induced TGFß activation required Gαq/11 signalling and serine protease activity, but was independent of integrins, suggesting an isoform-specific role for TGFß2 in this model. These data highlight a previously undescribed mechanism of cyclical stretch-induced Gαq/11-dependent TGFß2 signalling in mesenchymal cells, which is imperative for normal alveologenesis and maintenance of lung homeostasis.

How Many Cell Types Are in the Kidney and What Do They Do?

The kidney maintains electrolyte, water, and acid-base balance, eliminates foreign and waste compounds, regulates blood pressure, and secretes hormones. There are at least 16 different highly specialized epithelial cell types in the mammalian kidney. The number of specialized endothelial cells, immune cells, and interstitial cell types might even be larger. The concerted interplay between different cell types is critical for kidney function. Traditionally, cells were defined by their function or microscopical morphological appearance. With the advent of new single-cell modalities such as transcriptomics, epigenetics, metabolomics, and proteomics we are entering into a new era of cell type definition. This new technological revolution provides new opportunities to classify cells in the kidney and understand their functions.

METTL3-Mediated N6-methyladenosine mRNA Modification and cGAS-STING Pathway Activity in Kidney Fibrosis.

BACKGROUND: Chemical modifications on RNA profoundly impact RNA function and regulation. m6A, the most abundant RNA modification in eukaryotes, plays a pivotal role in diverse cellular processes and disease mechanisms. However, its importance is understudied in human chronic kidney disease (CKD) samples regarding its influence on pathological mechanisms. METHODS: LC-MS/MS and Methylated RNA Immunoprecipitation (MeRIP) sequencing were utilized to examine alterations in m6A levels and patterns in CKD samples. Overexpression of the m6A writer METTL3 in cultured kidney tubular cells was performed to confirm the impact of m6A in tubular cells and explore the biological functions of m6A modification on target genes. Additionally, tubule-specific deletion of Mettl3 (Ksp-Cre Mettl3f/f) mice and the use of anti-sense oligonucleotides inhibiting Mettl3 expression were utilized to reduce m6A modification in an animal kidney disease model. RESULTS: By examining 127 human CKD samples, we observed a significant increase in m6A modification and METTL3 expression in diseased kidneys. Epitranscriptomic analysis unveiled an enrichment of m6A modifications in transcripts associated with the activation of inflammatory signaling pathways, particularly the cGAS-STING pathway. m6A hypermethylation increased mRNA stability in cGAS and STING1, as well as elevated the expression of key proteins within the cGAS-STING pathway. Both the tubule-specific deletion of Mettl3 and the use of anti-sense oligonucleotides to inhibit Mettl3 expression protected mice from inflammation, reduced cytokine expression, decreased immune cell recruitment, and attenuated kidney fibrosis. CONCLUSIONS: Our research revealed heightened METTL3-mediated m6A modification in fibrotic kidneys, particularly enriching the cGAS-STING pathway. This hypermethylation increased mRNA stability for cGAS and STING1, leading to sterile inflammation and fibrosis.

The key role of altered tubule cell lipid metabolism in kidney disease development.

Kidney epithelial cells have very high energy requirements, which are largely met by fatty acid oxidation. Complex changes in lipid metabolism are observed in patients with kidney disease. Defects in fatty acid oxidation and increased lipid uptake, especially in the context of hyperlipidemia and proteinuria, contribute to this excess lipid build-up and exacerbate kidney disease development. Recent studies have also highlighted the role of increased de novo lipogenesis in kidney fibrosis. The defect in fatty acid oxidation causes energy starvation. Increased lipid uptake, synthesis, and lower fatty acid oxidation can cause toxic lipid build-up, reactive oxygen species generation, and mitochondrial damage. A better understanding of these metabolic processes may open new treatment avenues for kidney diseases by targeting lipid metabolism.

The 14th International Podocyte Conference 2023: from podocyte biology to glomerular medicine.

The 14th International Podocyte Conference took place in Philadelphia, Pennsylvania, USA from May 23 to 26, 2023. It commenced with an early-career researchers' meeting on May 23, providing young scientists with a platform to present and discuss their research findings. Throughout the main conference, 29 speakers across 9 sessions shared their insights on podocyte biology, glomerular medicine, novel technologic advancements, and translational approaches. Additionally, the event featured 3 keynote lectures addressing engineered chimeric antigen receptor T cell- and mRNA-based therapies and the use of biobanks for enhanced disease comprehension. Furthermore, 4 brief oral abstract sessions allowed scientists to present their findings to a broad audience. The program also included a panel discussion addressing the challenges of conducting human research within the American Black community. Remarkably, after a 5-year hiatus from in-person conferences, the 14th International Podocyte Conference successfully convened scientists from around the globe, fostering the presentation and discussion of crucial research findings, as summarized in this review. Furthermore, to ensure continuous and sustainable education, research, translation, and trial medicine related to podocyte and glomerular diseases for the benefit of patients, the International Society of Glomerular Disease was officially launched during the conference.

The role of complement in kidney disease: conclusions from a Kidney Disease: Improving Global Outcomes (KDIGO) Controversies Conference.

Uncontrolled complement activation can cause or contribute to glomerular injury in multiple kidney diseases. Although complement activation plays a causal role in atypical hemolytic uremic syndrome and C3 glomerulopathy, over the past decade, a rapidly accumulating body of evidence has shown a role for complement activation in multiple other kidney diseases, including diabetic nephropathy and several glomerulonephritides. The number of available complement inhibitor therapies has also increased during the same period. In 2022, Kidney Diseases: Improving Global Outcomes (KDIGO) convened a Controversies Conference, "The Role of Complement in Kidney Disease," to address the expanding role of complement dysregulation in the pathophysiology, diagnosis, and management of various glomerular diseases, diabetic nephropathy, and other forms of hemolytic uremic syndrome. Conference participants reviewed the evidence for complement playing a primary causal or secondary role in progression for several disease states and considered how evidence of complement involvement might inform management. Participating patients with various complement-mediated diseases and caregivers described concerns related to life planning, implications surrounding genetic testing, and the need for inclusive implementation of effective novel therapies into clinical practice. The value of biomarkers in monitoring disease course and the role of the glomerular microenvironment in complement response were examined, and key gaps in knowledge and research priorities were identified.

Genome-wide association study of hospitalized patients and acute kidney injury.

Acute kidney injury (AKI) is a common and devastating complication of hospitalization. Here, we identified genetic loci associated with AKI in patients hospitalized between 2002-2019 in the Million Veteran Program and data from Vanderbilt University Medical Center's BioVU. AKI was defined as meeting a modified KDIGO Stage 1 or more for two or more consecutive days or kidney replacement therapy. Control individuals were required to have one or more qualifying hospitalizations without AKI and no evidence of AKI during any other observed hospitalizations. Genome-wide association studies (GWAS), stratified by race, adjusting for sex, age, baseline estimated glomerular filtration rate (eGFR), and the top ten principal components of ancestry were conducted. Results were meta-analyzed using fixed effects models. In total, there were 54,488 patients with AKI and 138,051 non-AKI individuals included in the study. Two novel loci reached genome-wide significance in the meta-analysis: rs11642015 near the FTO locus on chromosome 16 (obesity traits) (odds ratio 1.07 (95% confidence interval, 1.05-1.09)) and rs4859682 near the SHROOM3 locus on chromosome 4 (glomerular filtration barrier integrity) (odds ratio 0.95 (95% confidence interval, 0.93-0.96)). These loci colocalized with previous studies of kidney function, and genetic correlation indicated significant shared genetic architecture between AKI and eGFR. Notably, the association at the FTO locus was attenuated after adjustment for BMI and diabetes, suggesting that this association may be partially driven by obesity. Both FTO and the SHROOM3 loci showed nominal evidence of replication from diagnostic-code-based summary statistics from UK Biobank, FinnGen, and Biobank Japan. Thus, our large GWA meta-analysis found two loci significantly associated with AKI suggesting genetics may explain some risk for AKI.

Gaining insight into metabolic diseases from human genetic discoveries.

Human large-scale genetic association studies have identified sequence variations at thousands of genetic risk loci that are more common in patients with diverse metabolic disease compared with healthy controls. While these genetic associations have been replicated in multiple large cohorts and sometimes can explain up to 50% of heritability, the molecular and cellular mechanisms affected by common genetic variation associated with metabolic disease remains mostly unknown. A variety of new genome-wide data types, in conjunction with novel biostatistical and computational analytical methodologies and foundational experimental technologies, are paving the way for a principled approach to systematic variant-to-function (V2F) studies for metabolic diseases, turning associated regions into causal variants, cell types and states of action, effector genes, and cellular and physiological mechanisms. Identification of new target genes and cellular programs for metabolic risk loci will improve mechanistic understanding of disease biology and identification of novel therapeutic strategies.

MRI-based quantification of whole-organ renal metabolic rate of oxygen.

During the early stages of diabetes, kidney oxygen utilization increases. The mismatch between oxygen demand and supply contributes to tissue hypoxia, a key driver of chronic kidney disease. Thus, whole-organ renal metabolic rate of oxygen (rMRO2 ) is a potentially valuable biomarker of kidney function. The key parameters required to determine rMRO2 include the renal blood flow rate (RBF) in the feeding artery and oxygen saturation in the draining renal vein (SvO2 ). However, there is currently no noninvasive method to quantify rMRO2 in absolute physiologic units. Here, a new MRI pulse sequence, Kidney Metabolism of Oxygen via T2 and Interleaved Velocity Encoding (K-MOTIVE), is described, along with evaluation of its performance in the human kidney in vivo. K-MOTIVE interleaves a phase-contrast module before a background-suppressed T2 -prepared balanced steady-state-free-precession (bSSFP) readout to measure RBF and SvO2 in a single breath-hold period of 22 s, yielding rMRO2 via Fick's principle. Variants of K-MOTIVE to evaluate alternative bSSFP readout strategies were studied. Kidney mass was manually determined from multislice gradient recalled echo images. Healthy subjects were recruited to quantify rMRO2 of the left kidney at 3-T field strength (N = 15). Assessments of repeat reproducibility and comparisons with individual measurements of RBF and SvO2 were performed, and the method's sensitivity was evaluated with a high-protein meal challenge (N = 8). K-MOTIVE yielded the following metabolic parameters: T2  = 157 ± 19 ms; SvO2  = 92% ± 6%; RBF = 400 ± 110 mL/min; and rMRO2  = 114 ± 117(µmol O2 /min)/100 g tissue. Reproducibility studies of T2 and RBF (parameters directly measured by K-MOTIVE) resulted in coefficients of variation less than 10% and intraclass correlation coefficients more than 0.75. The high-protein meal elicited an increase in rMRO2 , which was corroborated by serum biomarkers. The K-MOTIVE sequence measures SvO2 and RBF, the parameters necessary to quantify whole-organ rMRO2 , in a single breath-hold. The present work demonstrates that rMRO2 quantification is feasible with good reproducibility. rMRO2 is a potentially valuable physiological biomarker.

Tet2- and Tet3- Mediated Cytosine Hydroxymethylation in Six2 Progenitor Cells in Mice Is Critical for Nephron Progenitor Differentiation and Nephron Endowment.

SIGNIFICANCE STATEMENT: Epigenetic changes have been proposed to mediate nephron endowment during development, a critical determinant of future renal disease development. Hydroxymethyl cytosine, an epigenetic modification important for gene regulation, is abundant in the human kidney, but its physiologic role and the role of DNA demethylase enzymes encoded by the Tet1 , Tet2 , or Tet3 , which mediate cytosine hydroxymethylation, are unclear. By genetically deleting Tet1 , Tet2 , or Tet3 in nephron progenitors in mice, the authors showed that combined Tet2 and Tet3 loss in nephron progenitors cause defective kidney development, leading to kidney failure and perinatal death. Tet2 and Tet3 deletion also caused an alteration in demethylation and expression of genes critical for nephron formation. These findings establish that Tet2- and Tet3 -mediated cytosine hydroxymethylation in nephron progenitors plays a critical role in nephron endowment. BACKGROUND: Nephron endowment is a key determinant of hypertension and renal disease in later life. Epigenetic changes have been proposed to mediate fetal programming and nephron number. DNA cytosine methylation, which plays a critical role in gene regulation, is affected by proteins encoded by the ten-eleven translocation (TET) DNA demethylase gene family ( Tet1 , Tet2 , and Tet3 ), but the roles of TET proteins in kidney development and nephron endowment have not been characterized . METHODS: To study whether epigenetic changes-specifically, active DNA hydroxymethylation mediated by Tet1 , Tet2 , and Tet3- are necessary for nephron progenitor differentiation and nephron endowment, we generated mice with deletion of Tet1 , Tet2 , or Tet3 in Six2-positive nephron progenitors cells (NPCs). We then performed unbiased omics profiling, including whole-genome bisulfite sequencing on isolated Six2-positive NPCs and single-cell RNA sequencing on kidneys from newborn mice. RESULTS: We did not observe changes in kidney development or function in mice with NPC-specific deletion of Tet1 , Tet2 , Tet3 or Tet1 / Tet2 , or Tet1 / Tet3 . On the other hand, mice with combined Tet2 and Tet3 loss in Six2-positive NPCs failed to form nephrons, leading to kidney failure and perinatal death. Tet2 and Tet3 loss in Six2 -positive NPCs resulted in defective mesenchymal to epithelial transition and renal vesicle differentiation. Whole-genome bisulfite sequencing, single-cell RNA sequencing, and gene and protein expression analysis identified a defect in expression in multiple genes, including the WNT- ß -catenin signaling pathway, due to a failure in demethylation of these loci in the absence of Tet2 and Tet3 . CONCLUSIONS: These findings suggest that Tet2- and Tet3 -mediated active cytosine hydroxymethylation in NPCs play a key role in kidney development and nephron endowment.

Unified Mouse and Human Kidney Single-Cell Expression Atlas Reveal Commonalities and Differences in Disease States.

SIGNIFICANCE STATEMENT: Mouse models have been widely used to understand kidney disease pathomechanisms and play an important role in drug discovery. However, these models have not been systematically analyzed and compared. The authors characterized 18 different mouse kidney disease models at both bulk and single-cell gene expression levels and compared single-cell gene expression data from diabetic kidney disease (DKD) mice and from patients with DKD. Although single cell-level gene expression changes were mostly model-specific, different disease models showed similar changes when compared at a pathway level. The authors also found that changes in fractions of cell types are major drivers of bulk gene expression differences. Although the authors found only a small overlap of single cell-level gene expression changes between the mouse DKD model and patients, they observed consistent pathway-level changes. BACKGROUND: Mouse models have been widely used to understand kidney disease pathomechanisms and play an important role in drug discovery. However, these models have not been systematically analyzed and compared. METHODS: We analyzed single-cell RNA sequencing data (36 samples) and bulk gene expression data (42 samples) from 18 commonly used mouse kidney disease models. We compared single-nucleus RNA sequencing data from a mouse diabetic kidney disease model with data from patients with diabetic kidney disease and healthy controls. RESULTS: We generated a uniformly processed mouse single-cell atlas containing information for nearly 300,000 cells, identifying all major kidney cell types and states. Our analysis revealed that changes in fractions of cell types are major drivers of differences in bulk gene expression. Although gene expression changes at the single-cell level were mostly model-specific, different disease models showed similar changes when compared at a pathway level. Tensor decomposition analysis highlighted the important changes in proximal tubule cells in disease states. Specifically, we identified important alterations in expression of metabolic and inflammation-associated pathways. The mouse diabetic kidney disease model and patients with diabetic kidney disease shared only a small number of conserved cell type-specific differentially expressed genes, but we observed pathway-level activation patterns conserved between mouse and human diabetic kidney disease samples. CONCLUSIONS: This study provides a comprehensive mouse kidney single-cell atlas and defines gene expression commonalities and differences in disease states in mice. The results highlight the key role of cell heterogeneity in driving changes in bulk gene expression and the limited overlap of single-cell gene expression changes between animal models and patients, but they also reveal consistent pathway-level changes.

Unbiased Human Kidney Tissue Proteomics Identifies Matrix Metalloproteinase 7 as a Kidney Disease Biomarker.

SIGNIFICANCE STATEMENT: Although gene expression changes have been characterized in human diabetic kidney disease (DKD), unbiased tissue proteomics information for this condition is lacking. The authors conducted an unbiased aptamer-based proteomic analysis of samples from patients with DKD and healthy controls, identifying proteins with levels that associate with kidney function (eGFR) or fibrosis, after adjusting for key covariates. Overall, tissue gene expression only modestly correlated with tissue protein levels. Kidney protein and RNA levels of matrix metalloproteinase 7 (MMP7) strongly correlated with fibrosis and with eGFR. Single-cell RNA sequencing indicated that kidney tubule cells are an important source of MMP7. Furthermore, plasma MMP7 levels predicted future kidney function decline. These findings identify kidney tissue MMP7 as a biomarker of fibrosis and blood MMP7 as a biomarker for future kidney function decline. BACKGROUND: Diabetic kidney disease (DKD) is responsible for close to half of all ESKD cases. Although unbiased gene expression changes have been extensively characterized in human kidney tissue samples, unbiased protein-level information is not available. METHODS: We collected human kidney samples from 23 individuals with DKD and ten healthy controls, gathered associated clinical and demographics information, and implemented histologic analysis. We performed unbiased proteomics using the SomaScan platform and quantified the level of 1305 proteins and analyzed gene expression levels by bulk RNA and single-cell RNA sequencing (scRNA-seq). We validated protein levels in a separate cohort of kidney tissue samples as well as in 11,030 blood samples. RESULTS: Globally, human kidney transcript and protein levels showed only modest correlation. Our analysis identified 14 proteins with kidney tissue levels that correlated with eGFR and found that the levels of 152 proteins correlated with interstitial fibrosis. Of the identified proteins, matrix metalloprotease 7 (MMP7) showed the strongest association with both fibrosis and eGFR. The correlation between tissue MMP7 protein expression and kidney function was validated in external datasets. The levels of MMP7 RNA correlated with fibrosis in the primary and validation datasets. Findings from scRNA-seq pointed to proximal tubules, connecting tubules, and principal cells as likely cellular sources of increased tissue MMP7 expression. Furthermore, plasma MMP7 levels correlated not only with kidney function but also associated with prospective kidney function decline. CONCLUSIONS: Our findings, which underscore the value of human kidney tissue proteomics analysis, identify kidney tissue MMP7 as a diagnostic marker of kidney fibrosis and blood MMP7 as a biomarker for future kidney function decline.

Genetic Inhibition of APOL1 Pore-Forming Function Prevents APOL1-Mediated Kidney Disease.

SIGNIFICANCE STATEMENT: African Americans are at increased risk of CKD in part due to high-risk (HR) variants in the apolipoprotein L1 ( APOL1 ) gene, termed G1/G2. A different APOL1 variant, p.N264K , reduced the risk of CKD and ESKD among carriers of APOL1 HR variants to levels comparable with individuals with APOL1 low-risk variants in an analysis of 121,492 participants of African ancestry from the Million Veteran Program (MVP). Functional genetic studies in cell models showed that APOL1 p.N264K blocked APOL1 pore-forming function and ion channel conduction and reduced toxicity of APOL1 HR mutations. Pharmacologic inhibitors that mimic this mutation blocking APOL1 -mediated pore formation may be able to prevent and/or treat APOL1 -associated kidney disease. BACKGROUND: African Americans are at increased risk for nondiabetic CKD in part due to HR variants in the APOL1 gene. METHODS: We tested whether a different APOL1 variant, p.N264K , modified the association between APOL1 HR genotypes (two copies of G1/G2) and CKD in a cross-sectional analysis of 121,492 participants of African ancestry from the MVP. We replicated our findings in the Vanderbilt University Biobank ( n =14,386) and National Institutes of Health All of Us ( n =14,704). Primary outcome was CKD and secondary outcome was ESKD among nondiabetic patients. Primary analysis compared APOL1 HR genotypes with and without p.N264K . Secondary analyses included APOL1 low-risk genotypes and tested for interaction. In MVP, we performed sequential logistic regression models adjusting for demographics, comorbidities, medications, and ten principal components of ancestry. Functional genomic studies expressed APOL1 HR variants with and without APOL1 p.N264K in cell models. RESULTS: In the MVP cohort, 15,604 (12.8%) had two APOL1 HR variants, of which 582 (0.5%) also had APOL1 p.N264K . In MVP, 18,831 (15%) had CKD, 4177 (3%) had ESKD, and 34% had diabetes. MVP APOL1 HR, without p.N264K , was associated with increased odds of CKD (odds ratio [OR], 1.72; 95% confidence interval [CI], 1.60 to 1.85) and ESKD (OR, 3.94; 95% CI, 3.52 to 4.41). In MVP, APOL1 p.N264K mitigated the renal risk of APOL1 HR, in CKD (OR, 0.43; 95% CI, 0.28 to 0.65) and ESKD (OR, 0.19; CI 0.07 to 0.51). In the replication cohorts meta-analysis, APOL1 p.N264K mitigated the renal risk of APOL1 HR in CKD (OR, 0.40; 95% CI, 0.18 to 0.92) and ESKD (OR, 0.19; 95% CI, 0.05 to 0.79). In the mechanistic studies, APOL1 p.N264K blocked APOL1 pore-forming function and ion channel conduction and reduced toxicity of APOL1 HR variants. CONCLUSIONS: APOL1 p.N264K is associated with reduced risk of CKD and ESKD among carriers of APOL1 HR to levels comparable with individuals with APOL1 low-risk genotypes.

Deletion of protein tyrosine phosphatase SHP-1 restores SUMOylation of podocin and reverses the progression of diabetic kidney disease.

Both clinical and experimental data suggest that podocyte injury is involved in the onset and progression of diabetic kidney disease (DKD). Although the mechanisms underlying the development of podocyte loss are not completely understood, critical structural proteins such as podocin play a major role in podocyte survival and function. We have reported that the protein tyrosine phosphatase SHP-1 expression increased in podocytes of diabetic mice and glomeruli of patients with diabetes. However, the in vivo contribution of SHP-1 in podocytes is unknown. Conditional podocyte-specific SHP-1-deficient mice (Podo-SHP-1-/-) were generated to evaluate the impact of SHP-1 deletion at four weeks of age (early) prior to the onset of diabetes and after 20 weeks (late) of diabetes (DM; Ins2+/C96Y) on kidney function (albuminuria and glomerular filtration rate) and kidney pathology. Ablation of the SHP-1 gene specifically in podocytes prevented and even reversed the elevated albumin/creatinine ratio, glomerular filtration rate progression, mesangial cell expansion, glomerular hypertrophy, glomerular basement membrane thickening and podocyte foot process effacement induced by diabetes. Moreover, podocyte-specific deletion of SHP-1 at an early and late stage prevented diabetes-induced expression of collagen IV, fibronectin, transforming growth factor-ß, transforming protein RhoA, and serine/threonine kinase ROCK1, whereas it restored nephrin, podocin and cation channel TRPC6 expression. Mass spectrometry analysis revealed that SHP-1 reduced SUMO2 post-translational modification of podocin while podocyte-specific deletion of SHP-1 preserved slit diaphragm protein complexes in the diabetic context. Thus, our data uncovered a new role of SHP-1 in the regulation of cytoskeleton dynamics and slit diaphragm protein expression/stability, and its inhibition preserved podocyte function preventing DKD progression.

Antisense oligonucleotides ameliorate kidney dysfunction in podocyte-specific APOL1 risk variant mice.

Coding variants (named G1 and G2) in Apolipoprotein L1 (APOL1) can explain most excess risk of kidney disease observed in African American individuals. It has been proposed that risk variant APOL1 dose, such as increased risk variant APOL1 level serves as a trigger (second hit) for disease development. The goal of this study was to determine whether lowering risk variant APOL1 levels protects from disease development in a podocyte-specific transgenic mouse disease model. We administered antisense oligonucleotides (ASO) targeting APOL1 to podocyte-specific G2APOL1 mice and observed efficient reduction of APOL1 levels. APOL1 ASO1, which more efficiently lowered APOL1 transcript levels, protected mice from albuminuria, glomerulosclerosis, tubulointerstitial fibrosis, and renal failure. Administration of APOL1 ASO1 was effective even for established disease in the NEFTA-rtTA/TRE-G2APOL1 (NEFTA/G2APOL1) mice. We observed a strong correlation between APOL1 transcript level and disease severity. We concluded that APOL1 ASO1 may be an effective therapeutic approach for APOL1-associated glomerular disease.

ASEP: Gene-based detection of allele-specific expression across individuals in a population by RNA sequencing.

Allele-specific expression (ASE) analysis, which quantifies the relative expression of two alleles in a diploid individual, is a powerful tool for identifying cis-regulated gene expression variations that underlie phenotypic differences among individuals. Existing methods for gene-level ASE detection analyze one individual at a time, therefore failing to account for shared information across individuals. Failure to accommodate such shared information not only reduces power, but also makes it difficult to interpret results across individuals. However, when only RNA sequencing (RNA-seq) data are available, ASE detection across individuals is challenging because the data often include individuals that are either heterozygous or homozygous for the unobserved cis-regulatory SNP, leading to sample heterogeneity as only those heterozygous individuals are informative for ASE, whereas those homozygous individuals have balanced expression. To simultaneously model multi-individual information and account for such heterogeneity, we developed ASEP, a mixture model with subject-specific random effect to account for multi-SNP correlations within the same gene. ASEP only requires RNA-seq data, and is able to detect gene-level ASE under one condition and differential ASE between two conditions (e.g., pre- versus post-treatment). Extensive simulations demonstrated the convincing performance of ASEP under a wide range of scenarios. We applied ASEP to a human kidney RNA-seq dataset, identified ASE genes and validated our results with two published eQTL studies. We further applied ASEP to a human macrophage RNA-seq dataset, identified genes showing evidence of differential ASE between M0 and M1 macrophages, and confirmed our findings by results from cardiometabolic trait-relevant genome-wide association studies. To the best of our knowledge, ASEP is the first method for gene-level ASE detection at the population level that only requires the use of RNA-seq data. With the growing adoption of RNA-seq, we believe ASEP will be well-suited for various ASE studies for human diseases.

Systematic integrated analysis of genetic and epigenetic variation in diabetic kidney disease.

Poor metabolic control and host genetic predisposition are critical for diabetic kidney disease (DKD) development. The epigenome integrates information from sequence variations and metabolic alterations. Here, we performed a genome-wide methylome association analysis in 500 subjects with DKD from the Chronic Renal Insufficiency Cohort for DKD phenotypes, including glycemic control, albuminuria, kidney function, and kidney function decline. We show distinct methylation patterns associated with each phenotype. We define methylation variations that are associated with underlying nucleotide variations (methylation quantitative trait loci) and show that underlying genetic variations are important drivers of methylation changes. We implemented Bayesian multitrait colocalization analysis (moloc) and summary data-based Mendelian randomization to systematically annotate genomic regions that show association with kidney function, methylation, and gene expression. We prioritized 40 loci, where methylation and gene-expression changes likely mediate the genotype effect on kidney disease development. Functional annotation suggested the role of inflammation, specifically, apoptotic cell clearance and complement activation in kidney disease development. Our study defines methylation changes associated with DKD phenotypes, the key role of underlying genetic variations driving methylation variations, and prioritizes methylome and gene-expression changes that likely mediate the genotype effect on kidney disease pathogenesis.

Genome-wide meta-analysis and omics integration identifies novel genes associated with diabetic kidney disease.

AIMS/HYPOTHESIS: Diabetic kidney disease (DKD) is the leading cause of kidney failure and has a substantial genetic component. Our aim was to identify novel genetic factors and genes contributing to DKD by performing meta-analysis of previous genome-wide association studies (GWAS) on DKD and by integrating the results with renal transcriptomics datasets. METHODS: We performed GWAS meta-analyses using ten phenotypic definitions of DKD, including nearly 27,000 individuals with diabetes. Meta-analysis results were integrated with estimated quantitative trait locus data from human glomerular (N=119) and tubular (N=121) samples to perform transcriptome-wide association study. We also performed gene aggregate tests to jointly test all available common genetic markers within a gene, and combined the results with various kidney omics datasets. RESULTS: The meta-analysis identified a novel intronic variant (rs72831309) in the TENM2 gene associated with a lower risk of the combined chronic kidney disease (eGFR<60 ml/min per 1.73 m2) and DKD (microalbuminuria or worse) phenotype (p=9.8×10-9; although not withstanding correction for multiple testing, p>9.3×10-9). Gene-level analysis identified ten genes associated with DKD (COL20A1, DCLK1, EIF4E, PTPRN-RESP18, GPR158, INIP-SNX30, LSM14A and MFF; p<2.7×10-6). Integration of GWAS with human glomerular and tubular expression data demonstrated higher tubular AKIRIN2 gene expression in individuals with vs without DKD (p=1.1×10-6). The lead SNPs within six loci significantly altered DNA methylation of a nearby CpG site in kidneys (p<1.5×10-11). Expression of lead genes in kidney tubules or glomeruli correlated with relevant pathological phenotypes (e.g. TENM2 expression correlated positively with eGFR [p=1.6×10-8] and negatively with tubulointerstitial fibrosis [p=2.0×10-9], tubular DCLK1 expression correlated positively with fibrosis [p=7.4×10-16], and SNX30 expression correlated positively with eGFR [p=5.8×10-14] and negatively with fibrosis [p<2.0×10-16]). CONCLUSIONS/INTERPRETATION: Altogether, the results point to novel genes contributing to the pathogenesis of DKD. DATA AVAILABILITY: The GWAS meta-analysis results can be accessed via the type 1 and type 2 diabetes (T1D and T2D, respectively) and Common Metabolic Diseases (CMD) Knowledge Portals, and downloaded on their respective download pages ( https://t1d.hugeamp.org/downloads.html ; https://t2d.hugeamp.org/downloads.html ; https://hugeamp.org/downloads.html ).