Search

Diabetic kidney disease (DKD) is the leading cause of end stage kidney failure worldwide, of which cellular insulin resistance is a major driver. Here, we study key human kidney cell types implicated in DKD (podocytes, glomerular endothelial, mesangial and proximal tubular cells) in insulin sensitive and resistant conditions, and perform simultaneous transcriptomics and proteomics for integrated analysis. Our data is further compared with bulk- and single-cell transcriptomic kidney biopsy data from early- and advanced-stage DKD patient cohorts. We identify several consistent changes (individual genes, proteins, and molecular pathways) occurring across all insulin-resistant kidney cell types, together with cell-line-specific changes occurring in response to insulin resistance, which are replicated in DKD biopsies. This study provides a rich data resource to direct future studies in elucidating underlying kidney signalling pathways and potential therapeutic targets in DKD., Competing Interests: Competing interests The authors declare no competing interests., (© 2024. The Author(s).) more...

Introduction: Single-cell RNA sequencing (scRNA-seq) and single-nucleus RNA sequencing (snRNA-seq) provide valuable insights into the cellular states of kidney cells. However, the annotation of cell types often requires extensive domain expertise and time-consuming manual curation, limiting scalability and generalizability. To facilitate this process, we tested the performance of five supervised classification methods for automatic cell type annotation., Results: We analyzed publicly available sc/snRNA-seq datasets from five expert-annotated studies, comprising 62,120 cells from 79 kidney biopsy samples. Datasets were integrated by harmonizing cell type annotations across studies. Five different supervised machine learning algorithms (support vector machines, random forests, multilayer perceptrons, k-nearest neighbors, and extreme gradient boosting) were applied to automatically annotate cell types using four training datasets and one testing dataset. Performance metrics, including accuracy (F1 score) and rejection rates, were evaluated. All five machine learning algorithms demonstrated high accuracies, with a median F1 score of 0.94 and a median rejection rate of 1.8 %. The algorithms performed equally well across different datasets and successfully rejected cell types that were not present in the training data. However, F1 scores were lower when models trained primarily on scRNA-seq data were tested on snRNA-seq data., Conclusions: Despite limitations including the number of biopsy samples, our findings demonstrate that machine learning algorithms can accurately annotate a wide range of adult kidney cell types in scRNA-seq/snRNA-seq data. This approach has the potential to standardize cell type annotation and facilitate further research on cellular mechanisms underlying kidney disease., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2024 The Authors. Published by Elsevier Ltd.) more...

COVID-19 has been a significant public health concern for the last four years; however, little is known about the mechanisms that lead to severe COVID-associated kidney injury. In this multicenter study, we combined quantitative deep urinary proteomics and machine learning to predict severe acute outcomes in hospitalized COVID-19 patients. Using a 10-fold cross-validated random forest algorithm, we identified a set of urinary proteins that demonstrated predictive power for both discovery and validation set with 87% and 79% accuracy, respectively. These predictive urinary biomarkers were recapitulated in non-COVID acute kidney injury revealing overlapping injury mechanisms. We further combined orthogonal multiomics datasets to understand the mechanisms that drive severe COVID-associated kidney injury. Functional overlap and network analysis of urinary proteomics, plasma proteomics and urine sediment single-cell RNA sequencing showed that extracellular matrix and autophagy-associated pathways were uniquely impacted in severe COVID-19. Differentially abundant proteins associated with these pathways exhibited high expression in cells in the juxtamedullary nephron, endothelial cells, and podocytes, indicating that these kidney cell types could be potential targets. Further, single-cell transcriptomic analysis of kidney organoids infected with SARS-CoV-2 revealed dysregulation of extracellular matrix organization in multiple nephron segments, recapitulating the clinically observed fibrotic response across multiomics datasets. Ligand-receptor interaction analysis of the podocyte and tubule organoid clusters showed significant reduction and loss of interaction between integrins and basement membrane receptors in the infected kidney organoids. Collectively, these data suggest that extracellular matrix degradation and adhesion-associated mechanisms could be a main driver of COVID-associated kidney injury and severe outcomes. more...

Acute kidney injury (AKI) is a common condition that lacks effective treatments. In part, this shortcoming is due to an incomplete understanding of the genetic mechanisms that control pathogenesis and recovery. Identifying the molecular and genetic regulators unique to nephron segments that dictate vulnerability to injury and regenerative potential could lead to new therapeutic targets to treat ischemic kidney injury. Pax2 and Pax8 are homologous transcription factors with overlapping functions that are critical for kidney development and are re-activated in AKI. Here, we examined the role of Pax2 and Pax8 in recovery from ischemic AKI and found them upregulated after severe AKI and correlated with chronic injury. Surprisingly, proximal-tubule-selective deletion of Pax2 and Pax8 resulted in a less severe chronic injury phenotype. This effect was mediated by protection against the acute insult, similar to pre-conditioning. Prior to injury, Pax2 and Pax8 mutant mice develop a unique subpopulation of proximal tubule cells in the S3 segment that displayed features usually seen only in acute or chronic injury. The expression signature of these cells was strongly enriched with genes associated with other mechanisms of protection against ischemic AKI including caloric restriction, hypoxic pre-conditioning, and female sex. Thus, our results identified a novel role for Pax2 and Pax8 in mature proximal tubules that regulates critical genes and pathways involved in both the injury response and protection from ischemic AKI., (Copyright © 2023 International Society of Nephrology. Published by Elsevier Inc. All rights reserved.) more...

Introduction: Diabetic kidney disease (DKD) is a common cause of chronic kidney disease with around 25-40% of patients with diabetes being affected. The course of DKD is variable, and estimated glomerular filtration rate (eGFR) and albuminuria, the currently used clinical markers, are not able to accurately predict the individual disease trajectory, in particular in early stages of the disease. The aim of this study was to assess the association of urine levels of selected protein biomarkers with the progression of DKD at an early stage of disease., Methods: We measured 22 protein biomarkers using the Mesoscale Discovery platform in 461 urine samples of the PROVALID cohort, an observational study of patients with type 2 diabetes mellitus followed at the primary health care level for a minimum of 4 years. Odds ratios (ORs) were estimated for the effect of marker values above median on fast progression using unadjusted and adjusted logistic regression models. RNA expression at the single-cell level in kidney biopsy samples obtained from a cohort of young persons with type 2 diabetes mellitus was in addition determined for markers showing significant associations with disease progression., Results: Increased urinary levels of epidermal growth factor (EGF) were linked to lower odds of fast progression (defined as annual eGFR decline greater than 2.58 mL/min per 1.73 m2) with an OR of 0.60 (95% CI: 0.46, 0.78). The association with outcome was even stronger when adjusting for a set of 14 baseline clinical parameters including age, biological sex, eGFR, body mass index, albuminuria, and HbA1c. Elevated urinary levels of fatty acid-binding protein 3 (FABP3) and vascular cell adhesion molecule 1 (VCAM1) were each significantly associated with fast progression with an OR of 1.44 (95% CI: 1.11, 1.87) and an OR of 1.41 (95% CI: 1.08, 1.83), respectively. Enriched expression of EGF and FABP3 was observed in distal convoluted tubular cells and VCAM1 in parietal epithelial cells at single-cell level from biopsies of patients with early DKD., Conclusion: In summary, we show that lower urinary levels of EGF and higher urinary levels of FABP3 and VCAM1 are significantly associated with DKD progression in early-stage disease., (© 2024 The Author(s). Published by S. Karger AG, Basel.) more...

Diabetic kidney disease (DKD) can lead to end-stage kidney disease (ESKD) and mortality; however, few mechanistic biomarkers are available for high-risk patients, especially those without macroalbuminuria. Urine from participants with diabetes from the Chronic Renal Insufficiency Cohort (CRIC) study, the Singapore Study of Macro-angiopathy and Micro-vascular Reactivity in Type 2 Diabetes (SMART2D), and the American Indian Study determined whether urine adenine/creatinine ratio (UAdCR) could be a mechanistic biomarker for ESKD. ESKD and mortality were associated with the highest UAdCR tertile in the CRIC study and SMART2D. ESKD was associated with the highest UAdCR tertile in patients without macroalbuminuria in the CRIC study, SMART2D, and the American Indian study. Empagliflozin lowered UAdCR in nonmacroalbuminuric participants. Spatial metabolomics localized adenine to kidney pathology, and single-cell transcriptomics identified ribonucleoprotein biogenesis as a top pathway in proximal tubules of patients without macroalbuminuria, implicating mTOR. Adenine stimulated matrix in tubular cells via mTOR and stimulated mTOR in mouse kidneys. A specific inhibitor of adenine production was found to reduce kidney hypertrophy and kidney injury in diabetic mice. We propose that endogenous adenine may be a causative factor in DKD. more...

Acute kidney injury (AKI) is a common condition that lacks effective treatments. In part this shortcoming is due to an incomplete understanding of the genetic mechanisms that control pathogenesis and recovery. Pax2 and Pax8 are homologous transcription factors with overlapping functions that are critical for kidney development and are re-activated in AKI. In this report, we examined the role of Pax2 and Pax8 in recovery from ischemic AKI. We found that Pax2 and Pax8 are upregulated after severe AKI and correlate with chronic injury. Surprisingly, we then discovered that proximal-tubule-selective deletion of Pax2 and Pax8 resulted in a less severe chronic injury phenotype. This effect was mediated by protection against the acute insult, similar to preconditioning. Prior to injury, Pax2 and Pax8 mutant mice develop a unique subpopulation of S3 proximal tubule cells that display features usually seen only in acute or chronic injury. The expression signature of these cells was strongly enriched with genes associated with other mechanisms of protection against ischemic AKI including caloric restriction, hypoxic preconditioning, and female sex. Taken together, our results identify a novel role for Pax2 and Pax8 in mature proximal tubules that regulates critical genes and pathways involved in both injury response and protection from ischemic AKI., Translational Statement: Identifying the molecular and genetic regulators unique to the nephron that dictate vulnerability to injury and regenerative potential could lead to new therapeutic targets to treat ischemic kidney injury. Pax2 and Pax8 are two homologous nephron-specific transcription factors that are critical for kidney development and physiology. Here we report that proximal-tubule-selective depletion of Pax2 and Pax8 protects against both acute and chronic injury and induces an expression profile in the S3 proximal tubule with common features shared among diverse conditions that protect against ischemia. These findings highlight a new role for Pax proteins as potential therapeutic targets to treat AKI. more...

Underlying molecular mechanisms of the kidney protective effects of sodium glucose co-transporter 2 (SGLT2) inhibitors are not fully elucidated. Therefore, we studied the association between urinary epidermal growth factor (uEGF), a mitogenic factor involved in kidney repair, and kidney outcomes in patients with type 2 diabetes (T2D). The underlying molecular mechanisms of the SGLT2 inhibitor canagliflozin on EGF using single-cell RNA sequencing from kidney tissue were examined. Urinary EGF-to-creatinine ratio (uEGF/Cr) was measured in 3521 CANagliflozin cardioVascular Assessment Study (CANVAS) participants at baseline and week 52. Associations of uEGF/Cr with kidney outcome were assessed using multivariable-adjusted Cox regression models. Single-cell RNA sequencing was performed using protocol kidney biopsy tissue from ten young patients with T2D on SGLT2i, six patients with T2D on standard care only, and six healthy controls (HCs). In CANVAS, each doubling in baseline uEGF/Cr was associated with a 12% (95% confidence interval 1-22) decreased risk of kidney outcome. uEGF/Cr decreased after 52 weeks with placebo and remained stable with canagliflozin (between-group difference +7.3% (2.0-12.8). In young persons with T2D, EGF mRNA was primarily expressed in the thick ascending loop of Henle. Expression in biopsies from T2D without SGLT2i was significantly lower compared to HCs, whereas treatment with SGLT2i increased EGF levels closer to the healthy state. In young persons with T2D without SGLT2i, endothelin-1 emerged as a key regulator of the EGF co-expression network. SGLT2i treatment was associated with a shift towards normal EGF expression. Thus, decreased uEGF represents increased risk of kidney disease progression in patients with T2D. Canagliflozin increased kidney tissue expression of EGF and was associated with a downstream signaling cascade linked to tubular repair and reversal of tubular injury., (Copyright © 2023 International Society of Nephrology. Published by Elsevier Inc. All rights reserved.) more...

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. more...

Kidney organoids are a promising model to study kidney disease, but their use is constrained by limited knowledge of their functional protein expression profile. Here, we define the organoid proteome and transcriptome trajectories over culture duration and upon exposure to TNFα, a cytokine stressor. Older organoids increase deposition of extracellular matrix but decrease expression of glomerular proteins. Single cell transcriptome integration reveals that most proteome changes localize to podocytes, tubular and stromal cells. TNFα treatment of organoids results in 322 differentially expressed proteins, including cytokines and complement components. Transcript expression of these 322 proteins is significantly higher in individuals with poorer clinical outcomes in proteinuric kidney disease. Key TNFα-associated protein (C3 and VCAM1) expression is increased in both human tubular and organoid kidney cell populations, highlighting the potential for organoids to advance biomarker development. By integrating kidney organoid omic layers, incorporating a disease-relevant cytokine stressor and comparing with human data, we provide crucial evidence for the functional relevance of the kidney organoid model to human kidney disease., (© 2023. Springer Nature Limited.) more...

Understanding kidney disease relies on defining the complexity of cell types and states, their associated molecular profiles and interactions within tissue neighbourhoods 1> . 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., (© 2023. The Author(s).) more...

Diabetic kidney disease (DKD) can lead to end-stage kidney disease (ESKD) and mortality, however, few mechanistic biomarkers are available for high risk patients, especially those without macroalbuminuria. Urine from participants with diabetes from Chronic Renal Insufficiency Cohort (CRIC), Singapore Study of Macro-Angiopathy and Reactivity in Type 2 Diabetes (SMART2D), and the Pima Indian Study determined if urine adenine/creatinine ratio (UAdCR) could be a mechanistic biomarker for ESKD. ESKD and mortality were associated with the highest UAdCR tertile in CRIC (HR 1.57, 1.18, 2.10) and SMART2D (HR 1.77, 1.00, 3.12). ESKD was associated with the highest UAdCR tertile in patients without macroalbuminuria in CRIC (HR 2.36, 1.26, 4.39), SMART2D (HR 2.39, 1.08, 5.29), and Pima Indian study (HR 4.57, CI 1.37-13.34). Empagliflozin lowered UAdCR in non-macroalbuminuric participants. Spatial metabolomics localized adenine to kidney pathology and transcriptomics identified ribonucleoprotein biogenesis as a top pathway in proximal tubules of patients without macroalbuminuria, implicating mammalian target of rapamycin (mTOR). Adenine stimulated matrix in tubular cells via mTOR and stimulated mTOR in mouse kidneys. A specific inhibitor of adenine production was found to reduce kidney hypertrophy and kidney injury in diabetic mice. We propose that endogenous adenine may be a causative factor in DKD. more...

The molecular mechanisms of sodium-glucose cotransporter-2 (SGLT2) inhibitors (SGLT2i) remain incompletely understood. Single-cell RNA sequencing and morphometric data were collected from research kidney biopsies donated by young persons with type 2 diabetes (T2D), aged 12 to 21 years, and healthy controls (HCs). Participants with T2D were obese and had higher estimated glomerular filtration rates and mesangial and glomerular volumes than HCs. Ten T2D participants had been prescribed SGLT2i (T2Di[+]) and 6 not (T2Di[-]). Transcriptional profiles showed SGLT2 expression exclusively in the proximal tubular (PT) cluster with highest expression in T2Di(-) patients. However, transcriptional alterations with SGLT2i treatment were seen across nephron segments, particularly in the distal nephron. SGLT2i treatment was associated with suppression of transcripts in the glycolysis, gluconeogenesis, and tricarboxylic acid cycle pathways in PT, but had the opposite effect in thick ascending limb. Transcripts in the energy-sensitive mTORC1-signaling pathway returned toward HC levels in all tubular segments in T2Di(+), consistent with a diabetes mouse model treated with SGLT2i. Decreased levels of phosphorylated S6 protein in proximal and distal tubules in T2Di(+) patients confirmed changes in mTORC1 pathway activity. We propose that SGLT2i treatment benefits the kidneys by mitigating diabetes-induced metabolic perturbations via suppression of mTORC1 signaling in kidney tubules. more...

Hyperfiltration is a state of high glomerular filtration rate (GFR) observed in early diabetes that damages glomeruli, resulting in an iterative process of increasing filtration load on fewer and fewer remaining functional glomeruli. To delineate underlying cellular mechanisms of damage associated with hyperfiltration, transcriptional profiles of kidney biopsies from Pima Indians with type 2 diabetes with or without early-stage diabetic kidney disease were grouped into two hyperfiltration categories based on annual iothalamate GFR measurements. Twenty-six participants with a peak GFR measurement within two years of biopsy were categorized as the hyperfiltration group, and 26 in whom biopsy preceded peak GFR by over two years were considered pre-hyperfiltration. The hyperfiltration group had higher hemoglobin A1c, higher urine albumin-to-creatinine ratio, increased glomerular basement membrane width and lower podocyte density compared to the pre-hyperfiltration group. A glomerular 1240-gene transcriptional signature identified in the hyperfiltration group was enriched for endothelial stress response signaling genes, including endothelin-1, tec-kinase and transforming growth factor-β1 pathways, with the majority of the transcripts mapped to endothelial and inflammatory cell clusters in kidney single cell transcriptional data. Thus, our analysis reveals molecular pathomechanisms associated with hyperfiltration in early diabetic kidney disease involving putative ligand-receptor pairs with downstream intracellular targets linked to cellular crosstalk between endothelial and mesangial cells., (Copyright © 2022 International Society of Nephrology. All rights reserved.) more...

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. more...

Background: Molecular characterization of nephropathies may facilitate pathophysiologic insight, development of targeted therapeutics, and transcriptome-based disease classification. Although membranous nephropathy (MN) is a common cause of adult-onset nephrotic syndrome, the molecular pathways of kidney damage in MN require further definition., Methods: We applied a machine-learning framework to predict diagnosis on the basis of gene expression from the microdissected kidney tissue of participants in the Nephrotic Syndrome Study Network (NEPTUNE) cohort. We sought to identify differentially expressed genes between participants with MN versus those of other glomerulonephropathies across the NEPTUNE and European Renal cDNA Bank (ERCB) cohorts, to find MN-specific gene modules in a kidney-specific functional network, and to identify cell-type specificity of MN-specific genes using single-cell sequencing data from reference nephrectomy tissue., Results: Glomerular gene expression alone accurately separated participants with MN from those with other nephrotic syndrome etiologies. The top predictive classifier genes from NEPTUNE participants were also differentially expressed in the ERCB participants with MN. We identified a signature of 158 genes that are significantly differentially expressed in MN across both cohorts, finding 120 of these in a validation cohort. This signature is enriched in targets of transcription factor NF-κB. Clustering these MN-specific genes in a kidney-specific functional network uncovered modules with functional enrichments, including in ion transport, cell projection morphogenesis, regulation of adhesion, and wounding response. Expression data from reference nephrectomy tissue indicated 43% of these genes are most highly expressed by podocytes., Conclusions: These results suggest that, relative to other glomerulonephropathies, MN has a distinctive molecular signature that includes upregulation of many podocyte-expressed genes, provides a molecular snapshot of MN, and facilitates insight into MN's underlying pathophysiology., (Copyright © 2022 by the American Society of Nephrology.) more...

Studying isoform expression at the microscopic level has always been a challenging task. A classical example is kidney, where glomerular and tubulo-interstitial compartments carry out drastically different physiological functions and thus presumably their isoform expression also differs. We aim at developing an experimental and computational pipeline for identifying isoforms at microscopic structure-level. We microdissected glomerular and tubulo-interstitial compartments from healthy human kidney tissues from two cohorts. The two compartments were separately sequenced with the PacBio RS II platform. These transcripts were then validated using transcripts of the same samples by the traditional Illumina RNA-Seq protocol, distinct Illumina RNA-Seq short reads from European Renal cDNA Bank (ERCB) samples, and annotated GENCODE transcript list, thus identifying novel transcripts. We identified 14,739 and 14,259 annotated transcripts, and 17,268 and 13,118 potentially novel transcripts in the glomerular and tubulo-interstitial compartments, respectively. Of note, relying solely on either short or long reads would have resulted in many erroneous identifications. We identified distinct pathways involved in glomerular and tubulo-interstitial compartments at the isoform level, creating an important experimental and computational resource for the kidney research community., Competing Interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: YG receives personal payment from Eli Lilly and Company, which may have sponsored certain aspects of the study to MK. more...

Increased podocyte detachment begins immediately after kidney transplantation and is associated with long-term allograft failure. We hypothesized that cell-specific transcriptional changes in podocytes and glomerular endothelial cells after transplantation would offer mechanistic insights into the podocyte detachment process. To test this, we evaluated cell-specific transcriptional profiles of glomerular endothelial cells and podocytes from 14 patients of their first-year surveillance biopsies with normal histology from low immune risk recipients with no post-transplant complications and compared these to biopsies of 20 healthy living donor controls. Glomerular endothelial cells from these surveillance biopsies were enriched for genes related to fluid shear stress, angiogenesis, and interferon signaling. In podocytes, pathways were enriched for genes in response to growth factor signaling and actin cytoskeletal reorganization but also showed evidence of podocyte stress as indicated by reduced nephrin (adhesion protein) gene expression. In parallel, transcripts coding for proteins required to maintain podocyte adherence to the underlying glomerular basement membrane were downregulated, including the major glomerular podocyte integrin α3 and the actin cytoskeleton-related gene synaptopodin. The reduction in integrin α3 protein expression in surveillance biopsies was confirmed by immunoperoxidase staining. The combined growth and stress response of patient allografts post-transplantation paralleled similar changes in a rodent model of nephrectomy-induced glomerular hypertrophic stress that progress to develop proteinuria and glomerulosclerosis with shortened kidney life span. Thus, even among patients with apparently healthy allografts with no detectable histologic abnormality including alloimmune injury, transcriptomic changes reflecting cell stresses are already set in motion that could drive hypertrophy-associated glomerular disease progression., (Copyright © 2021 International Society of Nephrology. Published by Elsevier Inc. All rights reserved.) more...

Kidneys are critical target organs of COVID-19, but susceptibility and responses to infection remain poorly understood. Here, we combine SARS-CoV-2 variants with genome-edited kidney organoids and clinical data to investigate tropism, mechanism, and therapeutics. SARS-CoV-2 specifically infects organoid proximal tubules among diverse cell types. Infections produce replicating virus, apoptosis, and disrupted cell morphology, features of which are revealed in the context of polycystic kidney disease. Cross-validation of gene expression patterns in organoids reflects proteomic signatures of COVID-19 in the urine of critically ill patients indicating interferon pathway upregulation. SARS-CoV-2 viral variants alpha, beta, gamma, kappa, and delta exhibit comparable levels of infection in organoids. Infection is ameliorated in ACE2-/- organoids and blocked via treatment with de novo-designed spike binder peptides. Collectively, these studies clarify the impact of kidney infection in COVID-19 as reflected in organoids and clinical populations, enabling assessment of viral fitness and emerging therapies. more...

BACKGROUNDAssessment of chronic kidney disease (CKD) risk after acute kidney injury (AKI) is based on limited markers primarily reflecting glomerular function. We evaluated markers of cell integrity (EGF) and inflammation (monocyte chemoattractant protein-1, MCP-1) for predicting long-term kidney outcomes after cardiac surgery.METHODSWe measured EGF and MCP-1 in postoperative urine samples from 865 adults who underwent cardiac surgery at 2 sites in Canada and the United States and assessed EGF and MCP-1's associations with the composite outcome of CKD incidence or progression. We used single-cell RNA-Seq (scRNA-Seq) of AKI patient biopsies to perform transcriptomic analysis of programs corregulated with the associated genes.RESULTSOver a median (IQR) follow-up of 5.8 (4.2-7.1) years, 266 (30.8%) patients developed the composite CKD outcome. Postoperatively, higher levels of urinary EGF were protective and higher levels of MCP-1 were associated with the composite CKD outcome (adjusted HR 0.83, 95% CI 0.73-0.95 and 1.10, 95% CI 1.00-1.21, respectively). Intrarenal scRNA-Seq transcriptomes in patients with AKI-defined cell populations revealed concordant changes in EGF and MCP-1 levels and underlying molecular processes associated with loss of EGF expression and gain of CCL2 (encoding MCP-1) expression.CONCLUSIONUrinary EGF and MCP-1 were each independently associated with CKD after cardiac surgery. These markers may serve as noninvasive indicators of tubular damage, supported by tissue transcriptomes, and provide an opportunity for novel interventions in cardiac surgery.TRIAL REGISTRATIONClinicalTrials.gov NCT00774137.FUNDINGThe NIH funded the TRIBE-AKI Consortium and Kidney Precision Medicine Project. Yale O'Brien Kidney Center, American Heart Association, Patterson Trust Fund, Dr. Adam Linton Chair in Kidney Health Analytics, Canadian Institutes of Health Research, ICES, Ontario Ministry of Health and Long-Term Care, Academic Medical Organization of Southwestern Ontario, Schulich School of Medicine & Dentistry, Western University, Lawson Health Research Institute, Chan Zuckerberg Initiative Human Cell Atlas Kidney Seed Network. more...

Skeletal muscle wasting is commonly associated with chronic kidney disease (CKD), resulting in increased morbidity and mortality. However, the link between kidney and muscle function remains poorly understood. Here, we took a complementary interorgan approach to investigate skeletal muscle wasting in CKD. We identified increased production and elevated blood levels of soluble pro-cachectic factors, including activin A, directly linking experimental and human CKD to skeletal muscle wasting programs. Single-cell sequencing data identified the expression of activin A in specific kidney cell populations of fibroblasts and cells of the juxtaglomerular apparatus. We propose that persistent and increased kidney production of pro-cachectic factors, combined with a lack of kidney clearance, facilitates a vicious kidney/muscle signaling cycle, leading to exacerbated blood accumulation and, thereby, skeletal muscle wasting. Systemic pharmacological blockade of activin A using soluble activin receptor type IIB ligand trap as well as muscle-specific adeno-associated virus-mediated downregulation of its receptor ACVR2A/B prevented muscle wasting in different mouse models of experimental CKD, suggesting that activin A is a key factor in CKD-induced cachexia. In summary, we uncovered a crosstalk between kidney and muscle and propose modulation of activin signaling as a potential therapeutic strategy for skeletal muscle wasting in CKD. more...

Objective: Skin inflammation heralds systemic disease in juvenile myositis, yet we lack an understanding of pathogenic mechanisms driving skin inflammation in this disease. We undertook this study to define cutaneous gene expression signatures in juvenile myositis and identify key genes and pathways that differentiate skin disease in juvenile myositis from childhood-onset systemic lupus erythematosus (SLE)., Methods: We used formalin-fixed paraffin-embedded skin biopsy samples from 15 patients with juvenile myositis (9 lesional, 6 nonlesional), 5 patients with childhood-onset SLE, and 8 controls to perform transcriptomic analysis and identify significantly differentially expressed genes (DEGs; q ≤ 5%) between patient groups. We used Ingenuity Pathway Analysis (IPA) to highlight enriched biologic pathways and validated DEGs by immunohistochemistry and quantitative real-time polymerase chain reaction., Results: Comparison of lesional juvenile myositis to control samples revealed 221 DEGs, with the majority of up-regulated genes representing interferon (IFN)-stimulated genes. CXCL10, CXCL9, and IFI44L represented the top 3 DEGs (fold change 23.2, 13.3, and 13.0, respectively; q < 0.0001). IPA revealed IFN signaling as the top canonical pathway. When compared to childhood-onset SLE, lesional juvenile myositis skin shared a similar gene expression pattern, with only 28 unique DEGs, including FBLN2, CHKA, and SLURP1. Notably, patients with juvenile myositis who were positive for nuclear matrix protein 2 (NXP-2) autoantibodies exhibited the strongest IFN signature and also demonstrated the most extensive Mx-1 immunostaining, both in keratinocytes and perivascular regions., Conclusion: Lesional juvenile myositis skin demonstrates a striking IFN signature similar to that previously reported in juvenile myositis muscle and peripheral blood. Further investigation into the association of a higher IFN score with NXP-2 autoantibodies may provide insight into disease endotypes and pathogenesis., (© 2020, American College of Rheumatology.) more...

Podocytes are key to the glomerular filtration barrier by forming a slit diaphragm between interdigitating foot processes; however, the molecular details and functional importance of protein folding and degradation in the ER remain unknown. Here, we show that the SEL1L-HRD1 protein complex of ER-associated degradation (ERAD) is required for slit diaphragm formation and glomerular filtration function. SEL1L-HRD1 ERAD is highly expressed in podocytes of both mouse and human kidneys. Mice with podocyte-specific Sel1L deficiency develop podocytopathy and severe congenital nephrotic syndrome with an impaired slit diaphragm shortly after weaning and die prematurely, with a median lifespan of approximately 3 months. We show mechanistically that nephrin, a type 1 membrane protein causally linked to congenital nephrotic syndrome, is an endogenous ERAD substrate. ERAD deficiency attenuated the maturation of nascent nephrin, leading to its retention in the ER. We also show that various autosomal-recessive nephrin disease mutants were highly unstable and broken down by SEL1L-HRD1 ERAD, which attenuated the pathogenicity of the mutants toward the WT allele. This study uncovers a critical role of SEL1L-HRD1 ERAD in glomerular filtration barrier function and provides insights into the pathogenesis associated with autosomal-recessive disease mutants. more...

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. more...

COVID-19 morbidity and mortality are increased via unknown mechanisms in patients with diabetes and kidney disease. SARS-CoV-2 uses angiotensin-converting enzyme 2 (ACE2) for entry into host cells. Because ACE2 is a susceptibility factor for infection, we investigated how diabetic kidney disease and medications alter ACE2 receptor expression in kidneys. Single cell RNA profiling of kidney biopsies from healthy living donors and patients with diabetic kidney disease revealed ACE2 expression primarily in proximal tubular epithelial cells. This cell-specific localization was confirmed by in situ hybridization. ACE2 expression levels were unaltered by exposures to renin-angiotensin-aldosterone system inhibitors in diabetic kidney disease. Bayesian integrative analysis of a large compendium of public -omics datasets identified molecular network modules induced in ACE2-expressing proximal tubular epithelial cells in diabetic kidney disease (searchable at hb.flatironinstitute.org/covid-kidney) that were linked to viral entry, immune activation, endomembrane reorganization, and RNA processing. The diabetic kidney disease ACE2-positive proximal tubular epithelial cell module overlapped with expression patterns seen in SARS-CoV-2-infected cells. Similar cellular programs were seen in ACE2-positive proximal tubular epithelial cells obtained from urine samples of 13 hospitalized patients with COVID-19, suggesting a consistent ACE2-coregulated proximal tubular epithelial cell expression program that may interact with the SARS-CoV-2 infection processes. Thus SARS-CoV-2 receptor networks can seed further research into risk stratification and therapeutic strategies for COVID-19-related kidney damage., (Copyright © 2020 International Society of Nephrology. All rights reserved.) more...

Transected axons typically fail to regenerate in the central nervous system (CNS), resulting in chronic neurological disability in individuals with traumatic brain or spinal cord injury, glaucoma and ischemia-reperfusion injury of the eye. Although neuroinflammation is often depicted as detrimental, there is growing evidence that alternatively activated, reparative leukocyte subsets and their products can be deployed to improve neurological outcomes. In the current study, we identify a unique granulocyte subset, with characteristics of an immature neutrophil, that had neuroprotective properties and drove CNS axon regeneration in vivo, in part via secretion of a cocktail of growth factors. This pro-regenerative neutrophil promoted repair in the optic nerve and spinal cord, demonstrating its relevance across CNS compartments and neuronal populations. Our findings could ultimately lead to the development of new immunotherapies that reverse CNS damage and restore lost neurological function across a spectrum of diseases. more...

COVID-19 morbidity and mortality is increased in patients with diabetes and kidney disease via unknown mechanisms. SARS-CoV-2 uses angiotensin-converting enzyme 2 (ACE2) for entry into host cells. Since ACE2 is a susceptibility factor for infection, we investigated how diabetic kidney disease (DKD) and medications alter ACE2 receptor expression in kidneys. Single cell RNA profiling of healthy living donor (LD) and DKD kidney biopsies revealed ACE2 expression primarily in proximal tubular epithelial cells (PTEC). This cell specific localization was confirmed by in situ hybridization. ACE2 expression levels were unaltered by exposures to renin angiotensin aldosterone system inhibitors in DKD. Bayesian integrative analysis of a large compendium of public -omics datasets identified molecular network modules induced in ACE2-expressing PTEC in DKD (searchable at hb.flatironinstitute.org/covid-kidney) that were linked to viral entry, immune activation, endomembrane reorganization, and RNA processing. The DKD ACE2-positive PTEC module overlapped with expression patterns seen in SARS-CoV-2 infected cells. Similar cellular programs were seen in ACE2-positive PTEC obtained from urine samples of 13 COVID-19 patients who were hospitalized, suggesting a consistent ACE2-coregulated PTEC expression program that may interact with the SARS-CoV-2 infection processes. Thus SARS-CoV-2 receptor networks can seed further research into risk stratification and therapeutic strategies for COVID-19 related kidney damage. more...

Loss-of-function mutations in phospholipase C-ε1 (PLCE1) have been detected in patients with nephrotic syndrome, but other family members with the same mutation were asymptomatic, suggesting additional stressor are required to cause the full phenotype. Consistent with these observations, we determined that global Plce1 -deficient mice have histologically normal glomeruli and no albuminuria at baseline. Angiotensin II (ANG II) is known to induce glomerular damage in genetically susceptible individuals. Therefore, we tested whether ANG II enhances glomerular damage in Plce1 -deficient mice. ANG II increased blood pressure equally in Plce1 -deficient and wild-type littermates. Additionally, it led to 20-fold increased albuminuria and significantly more sclerotic glomeruli in Plce1 -deficient mice compared with wild-type littermates. Furthermore, Plce1- deficient mice demonstrated diffuse mesangial expansion, podocyte loss, and focal podocyte foot process effacement. To determine whether these effects are mediated by hypertension and hyperfiltration, rather than directly through ANG II, we raised blood pressure to a similar level using DOCA + salt + uninephrectomy and norepinephrine. This caused a fivefold increase in albuminuria in Plce1 -deficient mice and a significant increase in the number of sclerotic glomeruli. Consistent with previous findings in mice, we detected strong PLCE1 transcript expression in podocytes using single cell sequencing of human kidney tissue. In hemagglutinin-tagged Plce1 transgenic mice, Plce1 was detected in podocytes and also in glomerular arterioles using immunohistochemistry. Our data demonstrate that Plce1 deficiency in mice predisposes to glomerular damage secondary to hypertensive insults. more...

To define cellular mechanisms underlying kidney function and failure, the KPMP analyzes biopsy tissue in a multicenter research network to build cell-level process maps of the kidney. This study aimed to establish a single cell RNA sequencing strategy to use cell-level transcriptional profiles from kidney biopsies in KPMP to define molecular subtypes in glomerular diseases. Using multiple sources of adult human kidney reference tissue samples, 22,268 single cell profiles passed KPMP quality control parameters. Unbiased clustering resulted in 31 distinct cell clusters that were linked to kidney and immune cell types using specific cell markers. Focusing on endothelial cell phenotypes, in silico and in situ hybridization methods assigned 3 discrete endothelial cell clusters to distinct renal vascular beds. Transcripts defining glomerular endothelial cells (GEC) were evaluated in biopsies from patients with 10 different glomerular diseases in the NEPTUNE and European Renal cDNA Bank (ERCB) cohort studies. Highest GEC scores were observed in patients with focal segmental glomerulosclerosis (FSGS). Molecular endothelial signatures suggested 2 distinct FSGS patient subgroups with α-2 macroglobulin (A2M) as a key downstream mediator of the endothelial cell phenotype. Finally, glomerular A2M transcript levels associated with lower proteinuria remission rates, linking endothelial function with long-term outcome in FSGS. more...

Heterotopic ossification (HO) is an aberrant regenerative process with ectopic bone induction in response to musculoskeletal trauma, in which mesenchymal stem cells (MSC) differentiate into osteochondrogenic cells instead of myocytes or tenocytes. Despite frequent cases of hospitalized musculoskeletal trauma, the inflammatory responses and cell population dynamics that regulate subsequent wound healing and tissue regeneration are still unclear. Here we examine, using a mouse model of trauma-induced HO, the local microenvironment of the initial post-injury inflammatory response. Single cell transcriptome analyses identify distinct monocyte/macrophage populations at the injury site, with their dynamic changes over time elucidated using trajectory analyses. Mechanistically, transforming growth factor beta-1 (TGFβ1)-producing monocytes/macrophages are associated with HO and aberrant chondrogenic progenitor cell differentiation, while CD47-activating peptides that reduce systemic macrophage TGFβ levels and help ameliorate HO. Our data thus implicate CD47 activation as a therapeutic approach for modulating monocyte/macrophage phenotypes, MSC differentiation and HO formation during wound healing. more...

Background: Although diabetic kidney disease demonstrates both familial clustering and single nucleotide polymorphism heritability, the specific genetic factors influencing risk remain largely unknown., Methods: To identify genetic variants predisposing to diabetic kidney disease, we performed genome-wide association study (GWAS) analyses. Through collaboration with the Diabetes Nephropathy Collaborative Research Initiative, we assembled a large collection of type 1 diabetes cohorts with harmonized diabetic kidney disease phenotypes. We used a spectrum of ten diabetic kidney disease definitions based on albuminuria and renal function., Results: Our GWAS meta-analysis included association results for up to 19,406 individuals of European descent with type 1 diabetes. We identified 16 genome-wide significant risk loci. The variant with the strongest association (rs55703767) is a common missense mutation in the collagen type IV alpha 3 chain ( COL4A3) gene, which encodes a major structural component of the glomerular basement membrane (GBM). Mutations in COL4A3 are implicated in heritable nephropathies, including the progressive inherited nephropathy Alport syndrome. The rs55703767 minor allele (Asp326Tyr) is protective against several definitions of diabetic kidney disease, including albuminuria and ESKD, and demonstrated a significant association with GBM width; protective allele carriers had thinner GBM before any signs of kidney disease, and its effect was dependent on glycemia. Three other loci are in or near genes with known or suggestive involvement in this condition ( BMP7) or renal biology ( COLEC11 and DDR1 )., Conclusions: The 16 diabetic kidney disease-associated loci may provide novel insights into the pathogenesis of this condition and help identify potential biologic targets for prevention and treatment., (Copyright © 2019 by the American Society of Nephrology.) more...

An imbalance of nephroprotective factors and renal damaging molecules contributes to development and progression of chronic kidney disease (CKD). We investigated associations of renoprotective factor gene expression patterns with CKD severity and outcome. Gene expression profiles of 197 previously reported renoprotective factors were analyzed in a discovery cohort in renal biopsies of 63 CKD patients. Downregulation of dicarbonyl and L-xylulose reductase (DCXR) showed the strongest association with disease progression. This significant association was validated in an independent set of 225 patients with nephrotic syndrome from the multicenter NEPTUNE cohort. Reduced expression of DCXR was significantly associated with degree of histological damage as well as with lower estimated glomerular filtration rate and increased urinary protein levels. DCXR downregulation in CKD was confirmed in 3 publicly available transcriptomics data sets in the context of CKD. Expression of DCXR showed positive correlations to enzymes that are involved in dicarbonyl stress detoxification based on transcriptomics profiles. The sodium glucose cotransporter-2 (SGLT2) inhibitors canagliflozin and empagliflozin showed a beneficial effect on renal proximal tubular cells under diabetic stimuli-enhanced DCXR gene expression. In summary, lower expression of the renoprotective factor DCXR in renal tissue is associated with more severe disease and worse outcome in human CKD. more...

Podocyte injury is central to many forms of kidney disease, but transcriptional signatures reflecting podocyte injury and compensation mechanisms are challenging to analyze in vivo. Human kidney organoids derived from pluripotent stem cells (PSCs), a potentially new model for disease and regeneration, present an opportunity to explore the transcriptional plasticity of podocytes. Here, transcriptional profiling of more than 12,000 single cells from human PSC-derived kidney organoid cultures was used to identify robust and reproducible cell lineage gene expression signatures shared with developing human kidneys based on trajectory analysis. Surprisingly, the gene expression signature characteristic of developing glomerular epithelial cells was also observed in glomerular tissue from a kidney disease cohort. This signature correlated with proteinuria and inverse eGFR, and it was confirmed in an independent podocytopathy cohort. Three genes in particular were further characterized as potentially novel components of the glomerular disease signature. We conclude that cells in human PSC-derived kidney organoids reliably recapitulate the developmental transcriptional program of podocytes and other cell lineages in the human kidney and that transcriptional profiles seen in developing podocytes are reactivated in glomerular disease. Our findings demonstrate an approach to identifying potentially novel molecular programs involved in the pathogenesis of glomerulopathies. more...

The mammalian kidney develops through reciprocal interactions between the ureteric bud and the metanephric mesenchyme to give rise to the entire collecting system and the nephrons. Most of our knowledge of the developmental regulators driving this process arises from the study of gene expression and functional genetics in mice and other animal models. In order to shed light on human kidney development, we have used single-cell transcriptomics to characterize gene expression in different cell populations, and to study individual cell dynamics and lineage trajectories during development. Single-cell transcriptome analyses of 6414 cells from five individual specimens identified 11 initial clusters of specific renal cell types as defined by their gene expression profile. Further subclustering identifies progenitors, and mature and intermediate stages of differentiation for several renal lineages. Other lineages identified include mesangium, stroma, endothelial and immune cells. Novel markers for these cell types were revealed in the analysis, as were components of key signaling pathways driving renal development in animal models. Altogether, we provide a comprehensive and dynamic gene expression profile of the developing human kidney at the single-cell level., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2018. Published by The Company of Biologists Ltd.) more...

Expression quantitative trait loci (eQTL) studies illuminate the genetics of gene expression and, in disease research, can be particularly illuminating when using the tissues directly impacted by the condition. In nephrology, there is a paucity of eQTL studies of human kidney. Here, we used whole-genome sequencing (WGS) and microdissected glomerular (GLOM) and tubulointerstitial (TI) transcriptomes from 187 individuals with nephrotic syndrome (NS) to describe the eQTL landscape in these functionally distinct kidney structures. Using MatrixEQTL, we performed cis-eQTL analysis on GLOM (n = 136) and TI (n = 166). We used the Bayesian "Deterministic Approximation of Posteriors" (DAP) to fine-map these signals, eQTLBMA to discover GLOM- or TI-specific eQTLs, and single-cell RNA-seq data of control kidney tissue to identify the cell type specificity of significant eQTLs. We integrated eQTL data with an IgA Nephropathy (IgAN) GWAS to perform a transcriptome-wide association study (TWAS). We discovered 894 GLOM eQTLs and 1,767 TI eQTLs at FDR < 0.05. 14% and 19% of GLOM and TI eQTLs, respectively, had >1 independent signal associated with its expression. 12% and 26% of eQTLs were GLOM specific and TI specific, respectively. GLOM eQTLs were most significantly enriched in podocyte transcripts and TI eQTLs in proximal tubules. The IgAN TWAS identified significant GLOM and TI genes, primarily at the HLA region. In this study, we discovered GLOM and TI eQTLs, identified those that were tissue specific, deconvoluted them into cell-specific signals, and used them to characterize known GWAS alleles. These data are available for browsing and download via our eQTL browser, "nephQTL.", (Copyright © 2018 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.) more...

Organoids derived from human pluripotent stem cells are a potentially powerful tool for high-throughput screening (HTS), but the complexity of organoid cultures poses a significant challenge for miniaturization and automation. Here, we present a fully automated, HTS-compatible platform for enhanced differentiation and phenotyping of human kidney organoids. The entire 21-day protocol, from plating to differentiation to analysis, can be performed automatically by liquid-handling robots, or alternatively by manual pipetting. High-content imaging analysis reveals both dose-dependent and threshold effects during organoid differentiation. Immunofluorescence and single-cell RNA sequencing identify previously undetected parietal, interstitial, and partially differentiated compartments within organoids and define conditions that greatly expand the vascular endothelium. Chemical modulation of toxicity and disease phenotypes can be quantified for safety and efficacy prediction. Screening in gene-edited organoids in this system reveals an unexpected role for myosin in polycystic kidney disease. Organoids in HTS formats thus establish an attractive platform for multidimensional phenotypic screening., (Copyright © 2018 Elsevier Inc. All rights reserved.) more...

Previously, we showed that cortical mineralization is coordinately adjusted to mechanically offset external bone size differences between A/J (narrow) and C57BL/6J (wide) mouse femora to achieve whole bone strength equivalence at adulthood. The identity of the genes and their interactions that are responsible for establishing this homeostatic state (ie, canalization) remain unknown. We hypothesize that these inbred strains, whose interindividual differences in bone structure and material properties mimic that observed among humans, achieve functional homeostasis by differentially adjusting key molecular pathways regulating external bone size and mineralization throughout growth. The cortices of A/J and C57BL/6J male mouse femora were phenotyped and gene expression levels were assessed across growth (ie, ages 2, 4, 6, 8, 12, 16 weeks). A difference in total cross-sectional area (p < 0.01) and cortical tissue mineral density were apparent between mouse strains by age 2 weeks and maintained at adulthood (p < 0.01). These phenotypic dissimilarities corresponded to gene expression level differences among key regulatory pathways throughout growth. A/J mice had a 1.55- to 7.65-fold greater expression among genes inhibitory to Wnt pathway induction, whereas genes involved in cortical mineralization were largely upregulated 1.50- to 3.77-fold to compensate for their narrow diaphysis. Additionally, both mouse strains showed an upregulation among Wnt pathway antagonists corresponding to the onset of adult ambulation (ie, increased physiological loads). This contrasts with other studies showing an increase in Wnt pathway activation after functionally isolated, experimental in vivo loading regimens. A/J and C57BL/6J long bones provide a model to develop a systems-based approach to identify individual genes and the gene-gene interactions that contribute to trait differences between the strains while being involved in the process by which these traits are coordinately adjusted to establish similar levels of mechanical function, thus providing insight into the process of canalization. © 2017 American Society for Bone and Mineral Research., (© 2017 American Society for Bone and Mineral Research.) more...

Tens of thousands of splice isoforms of proteins have been catalogued as predicted sequences from transcripts in humans and other species. Relatively few have been characterized biochemically or structurally. With the extensive development of protein bioinformatics, the characterization and modeling of isoform features, isoform functions, and isoform-level networks have advanced notably. Here we present applications of the I-TASSER family of algorithms for folding and functional predictions and the IsoFunc, MIsoMine, and Hisonet data resources for isoform-level analyses of network and pathway-based functional predictions and protein-protein interactions. Hopefully, predictions and insights from protein bioinformatics will stimulate many experimental validation studies. more...

Background: High doses of the histone deacetylase inhibitor valproic acid (VPA, 150-400 mg/kg) improve outcomes in animal models of lethal insults. We are conducting a US Food and Drug Administration-approved Phase I, double-blind, placebo-controlled trial to evaluate the safety and tolerability of ascending doses of VPA in human volunteers. We hypothesized that VPA would induce significant changes in the proteome of healthy humans when given at doses lower than those used in prior animal studies., Methods: Peripheral blood mononuclear cells were obtained from three healthy subjects randomized to receive VPA (120 mg/kg over 1 hour) at baseline and at 4 and 8 hours following infusion. Detailed proteomic analysis was performed using 1D gel electrophoresis, liquid chromatography, and mass spectrometry. Proteins with differential expression were chosen for functional annotation and pathway analysis using Ingenuity Pathway Analysis (Qiagen GmbH, Hilden, Germany) and Panther Gene Ontology., Results: A total of 3,074 unique proteins were identified. The average number of proteins identified per sample was 1,716 ± 459. There were a total of 140 unique differentially expressed proteins (p < 0.05). There was a minor and inconsistent increase in histone and nonhistone protein acetylation. Functional annotation showed significant enrichment of apoptosis (p = 3.5E-43), cell death (p = 9.9E-72), proliferation of cells (p = 1.6E-40), dementia (p = 9.6E-40), amyloidosis (p = 6.3E-38), fatty acid metabolism (p = 4.6E-76), quantity of steroid (p = 4.2E-75), and cell movement (p = 1.9E-64)., Conclusions: Valproic acid induces significant changes to the proteome of healthy humans when given at a dose of 120 mg/kg. It alters the expression of key proteins and pathways, including those related to cell survival, without significant modification of protein acetylation. In the next part of the ongoing Phase I trial, we will study the effects of VPA on trauma patients in hemorrhagic shock., Level of Evidence: Therapeutic study, level V. more...

Ewing sarcomas are characterized by the presence of EWS/ETS fusion genes in the absence of other recurrent genetic alterations and mechanisms of tumor heterogeneity that contribute to disease progression remain unclear. Mutations in the Wnt/β-catenin pathway are rare in Ewing sarcoma but the Wnt pathway modulator LGR5 is often highly expressed, suggesting a potential role for the axis in tumor pathogenesis. We evaluated β-catenin and LGR5 expression in Ewing sarcoma cell lines and tumors and noted marked intra- and inter-tumor heterogeneity. Tumors with evidence of active Wnt/β-catenin signaling were associated with increased incidence of tumor relapse and worse overall survival. Paradoxically, RNA sequencing revealed a marked antagonism of EWS/ETS transcriptional activity in Wnt/β-catenin-activated tumor cells. Consistent with this, Wnt/β-catenin-activated cells displayed a phenotype that was reminiscent of Ewing sarcoma cells with partial EWS/ETS loss of function. Specifically, activation of Wnt/β-catenin induced alterations to the actin cytoskeleton, acquisition of a migratory phenotype, and upregulation of EWS/ETS-repressed genes. Notably, activation of Wnt/β-catenin signaling led to marked induction of tenascin C (TNC), an established promoter of cancer metastasis, and an EWS/ETS-repressed target gene. Loss of TNC function in Ewing sarcoma cells profoundly inhibited their migratory and metastatic potential. Our studies reveal that heterogeneous activation of Wnt/β-catenin signaling in subpopulations of tumor cells contributes to phenotypic heterogeneity and disease progression in Ewing sarcoma. Significantly, this is mediated, at least in part, by inhibition of EWS/ETS fusion protein function that results in derepression of metastasis-associated gene programs. Cancer Res; 76(17); 5040-53. ©2016 AACR., (©2016 American Association for Cancer Research.) more...

The vast majority of human multiexon genes undergo alternative splicing and produce a variety of splice variant transcripts and proteins, which can perform different functions. These protein-coding splice variants (PCSVs) greatly increase the functional diversity of proteins. Most functional annotation algorithms have been developed at the gene level; the lack of isoform-level gold standards is an important intellectual limitation for currently available machine learning algorithms. The accumulation of a large amount of RNA-seq data in the public domain greatly increases our ability to examine the functional annotation of genes at isoform level. In the present study, we used a multiple instance learning (MIL)-based approach for predicting the function of PCSVs. We used transcript-level expression values and gene-level functional associations from the Gene Ontology database. A support vector machine (SVM)-based 5-fold cross-validation technique was applied. Comparatively, genes with multiple PCSVs performed better than single PCSV genes, and performance also improved when more examples were available to train the models. We demonstrated our predictions using literature evidence of ADAM15, LMNA/C, and DMXL2 genes. All predictions have been implemented in a web resource called "IsoFunc", which is freely available for the global scientific community through http://guanlab.ccmb.med.umich.edu/isofunc . more...

The laboratory mouse is the primary mammalian species used for studying alternative splicing events. Recent studies have generated computational models to predict functions for splice isoforms in the mouse. However, the functional relationship network, describing the probability of splice isoforms participating in the same biological process or pathway, has not yet been studied in the mouse. Here we describe a rich genome-wide resource of mouse networks at the isoform level, which was generated using a unique framework that was originally developed to infer isoform functions. This network was built through integrating heterogeneous genomic and protein data, including RNA-seq, exon array, protein docking and pseudo-amino acid composition. Through simulation and cross-validation studies, we demonstrated the accuracy of the algorithm in predicting isoform-level functional relationships. We showed that this network enables the users to reveal functional differences of the isoforms of the same gene, as illustrated by literature evidence with Anxa6 (annexin a6) as an example. We expect this work will become a useful resource for the mouse genetics community to understand gene functions. The network is publicly available at: http://guanlab.ccmb.med.umich.edu/isoformnetwork. more...

The notion that dietary antioxidants can help fight cancer is popular. However, the mechanism(s) behind the effect of antioxidants in cancer is still unclear. Previous studies indicate that supplements can influence gene expression; however, all of these studies were focused on the coding/exonic gene expression. Studies are now emerging to highlight critical functional roles for RNAs expressed from the non-coding regions. This project was designed to study the effect of antioxidant supplements on non-coding intronic RNA expression in human cancers. Vitamin E, N-Acetyl cysteine (NAC) and Sulforaphane are commonly used supplements to prevent diseases including cancers. We studied the effect of these antioxidant supplements on the non-coding intronic RNA expression using publicly available datasets from a mouse model for lung cancer and prostate cancer cell lines. Although high throughput polyA-enriched RNA-Seq data characterize spliced coding mRNA regions, recent studies reveal the expression of reads from the non-coding intronic regions. Our analyses indicate that cancer cells have higher expression of introns compared to that of normal cells and that treatment with antioxidant supplements reduces the increased expression of introns of several genes. However, we did find high expression of introns of multiple genes including many oncogenes in the supplement treated groups compared to that of the control; this effect was distinct depending on the cell type and the supplement studied. Using RT-PCRs, we validated the expression of introns of two oncogenes, DLK1 and LRG1, known to be key players in lung cancer progression, and demonstrate changed intronic expression with supplement treatment in cancer cells. With regard to the antioxidant system, supplements did not change the intronic RNAs for endogenous antioxidant enzymes except for a significant decrease in the expression of superoxide dismutase (SOD) intronic RNA. Concurrently, we also found that a prolonged (48 h) exposure to Vitamin C, Vitamin E and Green tea extract reduced the enzymatic activity of SOD in lung cancer cells. The results from this study reveal that the antioxidant supplements have a significant effect on the intronic RNA expression of many genes including cancer genes that are not directly linked to the body's antioxidant system. It is important to study this novel effect of antioxidant supplements in detail as it may have a significant role in disease progression. more...

This paper summarizes the recent activities of the Chromosome-Centric Human Proteome Project (C-HPP) consortium, which develops new technologies to identify yet-to-be annotated proteins (termed "missing proteins") in biological samples that lack sufficient experimental evidence at the protein level for confident protein identification. The C-HPP also aims to identify new protein forms that may be caused by genetic variability, post-translational modifications, and alternative splicing. Proteogenomic data integration forms the basis of the C-HPP's activities; therefore, we have summarized some of the key approaches and their roles in the project. We present new analytical technologies that improve the chemical space and lower detection limits coupled to bioinformatics tools and some publicly available resources that can be used to improve data analysis or support the development of analytical assays. Most of this paper's content has been compiled from posters, slides, and discussions presented in the series of C-HPP workshops held during 2014. All data (posters, presentations) used are available at the C-HPP Wiki (http://c-hpp.webhosting.rug.nl/) and in the Supporting Information. more...

One goal of the Human Proteome Project is to identify at least one protein product for each of the ∼20,000 human protein-coding genes. As of October 2014, however, there are 3564 genes (18%) that have no or insufficient evidence of protein existence (PE), as curated by neXtProt; these comprise 2647 PE2-4 missing proteins and 616 PE5 dubious protein entries. We conducted a systematic examination of the 616 PE5 protein entries using cutting-edge protein structure and function modeling methods. Compared to a random sample of high-confidence PE1 proteins, the putative PE5 proteins were found to be over-represented in the membrane and cell surface proteins and peptides fold families. Detailed functional analyses show that most PE5 proteins, if expressed, would belong to transporters and receptors localized in the plasma membrane compartment. The results suggest that experimental difficulty in identifying membrane-bound proteins and peptides could have precluded their detection in mass spectrometry and that special enrichment techniques with improved sensitivity for membrane proteins could be important for the characterization of the PE5 "dark matter" of the human proteome. Finally, we identify 66 high scoring PE5 protein entries and find that six of them were reported in recent mass spectrometry databases; an illustrative annotation of these six is provided. This work illustrates a new approach to examine the potential folding and function of the dubious proteins comprising PE5, which we will next apply to the far larger group of missing proteins comprising PE2-4. more...

We have developed the web-based Michigan Proteome Visualization Tool (MI-PVT) to visualize and compare protein expression and isoform-level function across human chromosomes and tissues (http://guanlab.ccmb.med.umich.edu/mipvt). As proof of principle, we have populated the tool with Human Proteome Map (HPM) data. We were able to observe many biologically interesting features. From the vantage point of our chromosome 17 team, for example, we found more than 300 proteins from chromosome 17 expressed in each of the 30 tissues and cell types studied, with the highest number of expressed proteins being 685 in testis. Comparisons of expression levels across tissues showed low numbers of proteins expressed in esophagus, but esophagus had 12 cytoskeletal proteins coded on chromosome 17 with very high expression (>1000 spectral counts). This customized MI-PVT should be helpful for biologists to browse and study specific proteins and protein data sets across tissues and chromosomes. Users can upload any data of interest in MI-PVT for visualization. Our aim is to integrate extensive mass-spectrometric proteomic data into the tool to facilitate finding chromosome-centric protein expression and correlation across tissues. more...