Single-cell analysis of progenitor cell dynamics and lineage specification in the human fetal kidney

Polysialic acid regulates glomerular microvasculature formation by interaction with VEGF-A188 in mice

Vascular endothelial growth factor A (VEGF-A) is a key signalling protein that stimulates blood vessel development and repair. Its tight control is essential for organ development and tissue homeostasis. However, the complex regulatory network for balanced bioavailability of VEGF-A is not fully understood. Here, we assessed the role of the glycocalyx component polysialic acid (polySia) for kidney development and its potential interactions with VEGF-A isoforms, in vitro and in vivo, using mouse models of polySia deficiency. PolySia acts as negative regulator of cell adhesion, but also may interact with extracellular components. In murine kidney, polySia was identified on nephron progenitor and endothelial cell subsets in developing nephrons with declining expression during maturation. Loss of polySia in Ncam-/- mice revealed the neural cell adhesion molecule NCAM as major protein carrier. Both polysialyltransferase-negative and Ncam-/- mice displayed impaired glomerular microvasculature development with reduced endothelial cell numbers, reminiscent to the phenotype of mice with impaired VEGF-A signalling. In vitro, immobilized polySia specifically interacted with the VEGF-A188 isoform demonstrating an isoform-specific direct interaction. Single cell RNA sequencing data analysis of newborn mouse kidneys implicated activation of VEGF-A-signalling in polysialyltransferase-positive endothelial cells. Consistently, loss of polySia resulted in diminished VEGFR2 activation in perinatal kidney and human endothelial cells. At transcriptional level, the expression of polysialyltransferases and known polySia carrier proteins is conserved in human developing kidney. Together, these data demonstrate a direct impact of polySia on VEGF-A signalling with the perspective that polysialylation could be a therapeutic target to ameliorate microvasculature repair after renal injury.

Single cell transcription revealing key transcription factors in embryonic kidney development

The development of the kidney is a complex process involving the differentiation of renal progenitor cells into various specialized cell types, including renal tubular epithelial cells and podocytes. Understanding the molecular mechanisms that govern this differentiation process is crucial for unraveling the intricacies of kidney development. Recent advances in single-cell RNA sequencing technology have enabled researchers to explore the heterogeneity of renal progenitor cells at an unprecedented resolution, offering new insights into the distinct subpopulations and transcriptional landscapes within these cells. Our analysis revealed the presence of distinct subpopulations within renal progenitor cells. Furthermore, we identified key transcription factors that are crucial for the differentiation of these progenitors into renal tubular epithelial cells and podocytes. These findings provide new insights into the molecular mechanisms of kidney development, highlighting the role of specific transcription factors in the differentiation of renal progenitor cells. Although further research is needed to fully validate these findings and their therapeutic potential, this study provides a valuable foundation for future exploration of kidney regeneration strategies.

Parallel multiOMIC analysis reveals glutamine deprivation enhances directed differentiation of renal organoids

Metabolic pathways play a critical role in driving differentiation but remain poorly understood in the development of kidney organoids. In this study, parallel metabolite and transcriptome profiling of differentiating human pluripotent stem cells (hPSCs) to multicellular renal organoids revealed key metabolic drivers of the differentiation process. In the early stage, transitioning from hPSCs to nephron progenitor cells (NPCs), both the glutamine and the alanine-aspartate-glutamate pathways changed significantly, as detected by enrichment and pathway impact analyses. Intriguingly, hPSCs maintained their ability to generate NPCs, even when deprived of both glutamine and glutamate. Surprisingly, single cell RNA-Seq analysis detected enhanced maturation and enrichment for podocytes under glutamine-deprived conditions. Together, these findings illustrate a novel role of glutamine metabolism in regulating podocyte development.

Integrative transcriptome analysis reveals Serpine2 promotes glomerular mesangial cell proliferation and extracellular matrix accumulation via activating ERK1/2 signalling pathway in diabetic nephropathy

BACKGROUND

Diabetic nephropathy (DN) is one of the main causes of end-stage renal disease (ESRD), but its mechanism has not been clearly studied. We utilized integrative transcriptome analysis to explore the pathogenesis of DN.

METHODS

We conducted an analysis by combining bulk dataset and single-cell transcriptome dataset. Through this approach, we identified that Serpine2 may regulate the 'collagen-containing extracellular matrix' pathway involved in DN. Subsequently, we established DN animal and cell models using db/db mice and mesangial cells (MCs) to validate the role of Serpine2 in DN. In the animal model, we detected the expression level of Serpine2 in DN using western blotting (WB) and immunofluorescence (IF) assays. To further clarify the molecular mechanism of Serpine2 in DN, we knocked down Serpine2 and observed its effects on MCs proliferation and extracellular matrix (ECM) accumulation.

RESULTS

Our single-cell analysis of DN models highlighted a pivotal role for MCs in the disease's initiation. Next, through Cytoscape analysis of differentially expressed genes (DEGs) in MCs, we identified the following 10 hub genes: Acta2, Angpt2, Ccn1, Col4a1, Col4a2, Col8a1, Kdr, Thbs1, Tpm4 and Serpine2. Subsequently, we identified that Serpine2 and Kdr were also significantly DEGs in the bulk analysis of glomeruli. Additionally, our integrated gene set enrichment analysis of bulk dataset and single-cell RNA dataset revealed that the 'collagen-containing extracellular matrix' was a key pathway in DN progression. Serpine2 was one of the crucial genes involved in regulating this pathway. Therefore, we speculated that the regulation of the 'collagen-containing extracellular matrix' pathway by Serpine2 was an important mechanism. Importantly, WB and IF staining confirmed that Serpine2 expression was upregulated in the MCs of diabetic mice. Knockdown of Serpine2 in cultured MCs alleviated high-glucose-induced excessive MCs proliferation and ECM accumulation. Finally, we found that ERK agonist Ro 67-7476 eliminated the effect of Serpine2 siRNA.

CONCLUSIONS

In summary, Serpine2 regulates MCs proliferation and ECM synthesis through activation of the ERK1/2 pathway, which is an important pathogenesis mechanism of DN. These findings offer fresh perspectives on the mechanisms of glomerulosclerosis in DN pathogenesis and may provide new targets for treating DN.

Single-cell RNA sequencing in diabetic kidney disease: a literature review

Diabetic kidney disease (DKD) is the leading cause of end-stage renal disease (ESRD), and its pathogenesis has not been clarified. Current research suggests that DKD involves multiple cell types and extra-renal factors, and it is particularly important to clarify the pathogenesis and identify new therapeutic targets. Single-cell RNA sequencing (scRNA-seq) technology is high-throughput sequencing of the transcriptomes of individual cells at the single-cell level, which is an effective technology for exploring the development of diseases by comparing genetic information, reflecting the differences in genetic information between cells, and identifying different cell subpopulations. Accumulating evidence supports the role of scRNA-seq in revealing the pathogenesis of diabetes and strengthening our understanding of the molecular mechanisms of DKD. We reviewed the scRNA-seq data this time. Then, we analyzed and discussed the applications of scRNA-seq technology in DKD research, including annotation of cell types, identification of novel cell types (or subtypes), identification of intercellular communication, analysis of cell differentiation trajectories, gene expression detection, and analysis of gene regulatory networks, and lastly, we explored the future perspectives of scRNA-seq technology in DKD research.

Proximal tubule cell maturation rate and function are controlled by PPARα signaling in kidney organoids

Human pluripotent stem cell-derived kidney organoids are expected to be a useful tool for new drug discoveries, however, the immaturation of kidney organoids causes difficulties in recapitulating renal pharmacokinetics using organoids. Here, we performed time-course single-cell RNA sequencing of kidney organoids and revealed cell heterogeneity in the maturation rate of the proximal tubule. An unbiased analysis to identify upstream targets of genes that are expressed differentially between cells with low and high maturation rates revealed a higher activation of PPARα signaling in rapidly maturing cells. Treatment with a combination of a PPARα agonist and an RXRα agonist induced genes related to proximal tubule maturation and increased the capacity for protein uptake as well as the sensitivity to nephrotoxicity by cisplatin. This method to promote the maturation rate of proximal tubule cells has the potential to be utilized in microphysiological systems to recapitulate proximal tubule functions and to screen nephrotoxic drugs.

Comparative single-cell analyses identify shared and divergent features of human and mouse kidney development

The mammalian kidney maintains fluid homeostasis through diverse epithelial cell types generated from nephron and ureteric progenitor cells. To extend a developmental understanding of the kidney's epithelial networks, we compared chromatin organization (single-nuclear assay for transposase-accessible chromatin sequencing [ATAC-seq]; 112,864 nuclei) and gene expression (single-cell/nuclear RNA sequencing [RNA-seq]; 109,477 cells/nuclei) in the developing human (10.6-17.6 weeks; n = 10) and mouse (post-natal day [P]0; n = 10) kidney, supplementing analysis with published mouse datasets from earlier stages. Single-cell/nuclear datasets were analyzed at a species level, and then nephron and ureteric cellular lineages were extracted and integrated into a common, cross-species, multimodal dataset. Comparative computational analyses identified conserved and divergent features of chromatin organization and linked gene activity, identifying species-specific and cell-type-specific regulatory programs. In situ validation of human-enriched gene activity points to human-specific signaling interactions in kidney development. Further, human-specific enhancer regions were linked to kidney diseases through genome-wide association studies (GWASs), highlighting the potential for clinical insight from developmental modeling.

Attenuated kidney oxidative metabolism in young adults with type 1 diabetes

BACKGROUNDIn type 1 diabetes (T1D), impaired insulin sensitivity may contribute to the development of diabetic kidney disease (DKD) through alterations in kidney oxidative metabolism.METHODSYoung adults with T1D (n = 30) and healthy controls (HCs) (n = 20) underwent hyperinsulinemic-euglycemic clamp studies, MRI, 11C-acetate PET, kidney biopsies, single-cell RNA-Seq, and spatial metabolomics to assess this relationship.RESULTSParticipants with T1D had significantly higher glomerular basement membrane (GBM) thickness compared with HCs. T1D participants exhibited lower insulin sensitivity and cortical oxidative metabolism, correlating with higher insulin sensitivity. Proximal tubular transcripts of TCA cycle and oxidative phosphorylation enzymes were lower in T1D. Spatial metabolomics showed reductions in tubular TCA cycle intermediates, indicating mitochondrial dysfunction. The Slingshot algorithm identified a lineage of proximal tubular cells progressing from stable to adaptive/maladaptive subtypes, using pseudotime trajectory analysis, which computationally orders cells along a continuum of states. This analysis revealed distinct distribution patterns between T1D and HCs, with attenuated oxidative metabolism in T1D attributed to a greater proportion of adaptive/maladaptive subtypes with low expression of TCA cycle and oxidative phosphorylation transcripts. Pseudotime progression associated with higher HbA1c, BMI, and GBM, and lower insulin sensitivity and cortical oxidative metabolism.CONCLUSIONThese early structural and metabolic changes in T1D kidneys may precede clinical DKD.TRIAL REGISTRATIONClinicalTrials.gov NCT04074668.FUNDINGUniversity of Michigan O'Brien Kidney Translational Core Center grant (P30 DK081943); CROCODILE studies by National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) (P30 DK116073), Juvenile Diabetes Research Foundation (JDRF) (2-SRA-2019-845-S-B), Boettcher Foundation, Intramural Research Program at NIDDK and Centers for Disease Control and Prevention (CKD Initiative) under Inter-Agency Agreement #21FED2100157DPG.

INF2 formin variants linked to human inherited kidney disease reprogram the transcriptome, causing mitotic chaos and cell death

Mutations in the human INF2 gene cause autosomal dominant focal segmental glomerulosclerosis (FSGS)-a condition characterized by podocyte loss, scarring, and subsequent kidney degeneration. To understand INF2-linked pathogenicity, we examined the effect of pathogenic INF2 on renal epithelial cell lines and human primary podocytes. Our study revealed an increased incidence of mitotic cells with surplus microtubule-organizing centers fostering multipolar spindle assembly, leading to nuclear abnormalities, particularly multi-micronucleation. The levels of expression of exogenous pathogenic INF2 were similar to those of endogenous INF2. The aberrant nuclear phenotypes were observed regardless of the expression method used (retrovirus infection or plasmid transfection) or the promoter (LTR or CMV) used, and were absent with exogenous wild type INF2 expression. This indicates that the effect of pathogenic INF2 is not due to overexpression or experimental cell manipulation, but instead to the intrinsic properties of pathogenic INF2. Inactivation of the INF2 catalytic domain prevented aberrant nuclei formation. Pathogenic INF2 triggered the translocation of the transcriptional cofactor MRTF into the nucleus. RNA sequencing revealed a profound alteration in the transcriptome that could be primarily attributed to the sustained activation of the MRTF-SRF transcriptional complex. Cells eventually underwent mitotic catastrophe and death. Reducing MRTF-SRF activation mitigated multi-micronucleation, reducing the extent of cell death. Our results, if validated in animal models, could provide insights into the mechanism driving glomerular degeneration in INF2-linked FSGS and may suggest potential therapeutic strategies for impeding FSGS progression.

Selective induction of human renal interstitial progenitor-like cell lineages from iPSCs reveals development of mesangial and EPO-producing cells

Recent regenerative studies using human pluripotent stem cells (hPSCs) have developed multiple kidney-lineage cells and organoids. However, to further form functional segments of the kidney, interactions of epithelial and interstitial cells are required. Here we describe a selective differentiation of renal interstitial progenitor-like cells (IPLCs) from human induced pluripotent stem cells (hiPSCs) by modifying our previous induction method for nephron progenitor cells (NPCs) and analyzing mouse embryonic interstitial progenitor cell (IPC) development. Our IPLCs combined with hiPSC-derived NPCs and nephric duct cells form nephrogenic niche- and mesangium-like structures in vitro. Furthermore, we successfully induce hiPSC-derived IPLCs to differentiate into mesangial and erythropoietin-producing cell lineages in vitro by screening differentiation-inducing factors and confirm that p38 MAPK, hypoxia, and VEGF signaling pathways are involved in the differentiation of mesangial-lineage cells. These findings indicate that our IPC-lineage induction method contributes to kidney regeneration and developmental research.

Glucocorticoid- and pioglitazone-induced proteinuria reduction in experimental NS both correlate with glomerular ECM modulation

Idiopathic nephrotic syndrome (NS) is a common glomerular disease. Although glucocorticoids (GC) are the primary treatment, the PPARγ agonist pioglitazone (Pio) also reduces proteinuria in patients with NS and directly protects podocytes from injury. Because both drugs reduce proteinuria, we hypothesized these effects result from overlapping transcriptional patterns. Systems biology approaches compared glomerular transcriptomes from rats with PAN-induced NS treated with GC vs. Pio and identified 29 commonly regulated genes-of-interest, primarily involved in extracellular matrix (ECM) remodeling. Correlation with clinical idiopathic NS patient datasets confirmed glomerular ECM dysregulation as a potential mechanism of injury. Cellular deconvolution in silico revealed GC- and Pio-induced amelioration of altered genes primarily within podocytes and mesangial cells. While validation studies are indicated, these analyses identified molecular pathways involved in the early stages of NS (prior to scarring), suggesting that targeting glomerular ECM dysregulation may enable a future non-immunosuppressive approach for proteinuria reduction in idiopathic NS.

Nephron progenitors rhythmically alternate between renewal and differentiation phases that synchronize with kidney branching morphogenesis

The mammalian kidney achieves massive parallelization of function by exponentially duplicating nephron-forming niches during development. Each niche caps a tip of the ureteric bud epithelium (the future urinary collecting duct tree) as it undergoes branching morphogenesis, while nephron progenitors within niches balance self-renewal and differentiation to early nephron cells. Nephron formation rate approximately matches branching rate over a large fraction of mouse gestation, yet the nature of this apparent pace-maker is unknown. Here we correlate spatial transcriptomics data with branching 'life-cycle' to discover rhythmically alternating signatures of nephron progenitor differentiation and renewal across Wnt, Hippo-Yap, retinoic acid (RA), and other pathways. We then find in human stem-cell derived nephron progenitor organoids that Wnt/β-catenin-induced differentiation is converted to a renewal signal when it temporally overlaps with YAP activation. Similar experiments using RA activation indicate a role in setting nephron progenitor exit from the naive state, the spatial extent of differentiation, and nephron segment bias. Together the data suggest that nephron progenitor interpretation of consistent Wnt/β-catenin differentiation signaling in the niche may be modified by rhythmic activity in ancillary pathways to set the pace of nephron formation. This would synchronize nephron formation with ureteric bud branching, which creates new sites for nephron condensation. Our data bring temporal resolution to the renewal vs. differentiation balance in the nephrogenic niche and inform new strategies to achieve self-sustaining nephron formation in synthetic human kidney tissues.

Lung patterning: Is a distal-to-proximal gradient of cell allocation and fate decision a general paradigm?: A gradient of distal-to-proximal distribution and differentiation of tip progenitors produces distinct compartments in the lung

Recent studies support a model in which the progeny of SOX9+ epithelial progenitors at the distal tip of lung branches undergo cell allocation and differentiation sequentially along the distal-to-proximal axis. Concomitant with the elongation and ramification of lung branches, the descendants of the distal SOX9+ progenitors are distributed proximally, express SOX2, and differentiate into cell types in the conducting airways. Amid subsequent sacculation, the distal SOX9+ progenitors generate alveolar epithelial cells to form alveoli. Sequential cell allocation and differentiation are integrated with the branching process to generate a functional branching organ. This review focuses on the roles of SOX9+ cells as precursors for new branches, as the source of various cell types in the conducting airways, and as progenitors of the alveolar epithelium. All of these processes are controlled by multiple signaling pathways. Many mouse mutants with defective lung branching contain underlying defects in one or more steps of cell allocation and differentiation of SOX9+ progenitors. This model provides a framework to understand the molecular basis of lung phenotypes and to elucidate the molecular mechanisms of lung patterning. It builds a foundation on which comparing and contrasting the mechanisms employed by different branching organs in diverse species can be made.

Generation of renal tubular organoids from adult SOX9+ kidney progenitor cells

The pathogenesis of several kidney diseases results in the eventual destruction of the renal tubular system, which can progress to end-stage renal disease. Previous studies have demonstrated the involvement of a population of SOX9-positive cells in kidney regeneration and repair process following kidney injury. However, the ability of these cells to autonomously generate kidney organoids has never been investigated. Here, we isolated SOX9+ kidney progenitor cells (KPCs) from both mice and humans and tested their differentiation potential in vitro. The data showed that the human SOX9+ KPC could self-assemble into organoids with kidney-like morphology. We also used single-cell RNA sequencing to characterize the organoid cell populations and identified four distinct types of renal tubular cells. Compared to the induced pluripotent stem cell-derived kidney organoids, KPC demonstrated more tubular differentiation potential but failed to differentiate into glomerular cells. KPC-derived organoid formation involved the expression of genes related to metanephric development and followed a similar mechanism to renal injury repair in acute kidney injury patients. Altogether, our study provided a potentially useful approach to generating kidney tubular organoids for future application.

Sctensor detects many-to-many cell-cell interactions from single cell RNA-sequencing data

BACKGROUND

Complex biological systems are described as a multitude of cell-cell interactions (CCIs). Recent single-cell RNA-sequencing studies focus on CCIs based on ligand-receptor (L-R) gene co-expression but the analytical methods are not appropriate to detect many-to-many CCIs.

RESULTS

In this work, we propose scTensor, a novel method for extracting representative triadic relationships (or hypergraphs), which include ligand-expression, receptor-expression, and related L-R pairs.

CONCLUSIONS

Through extensive studies with simulated and empirical datasets, we have shown that scTensor can detect some hypergraphs that cannot be detected using conventional CCI detection methods, especially when they include many-to-many relationships. scTensor is implemented as a freely available R/Bioconductor package.

Sodium glucose co-transporter 2 inhibition increases epidermal growth factor expression and improves outcomes in patients with type 2 diabetes

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.

Prioritization of Monogenic Congenital Anomalies of the Kidney and Urinary Tract Candidate Genes with Existing Single-Cell Transcriptomics Data of the Human Fetal Kidney

INTRODUCTION

Congenital anomalies of the kidney and urinary tract (CAKUT) are the most common cause of chronic kidney disease in the first 3 decades of life. Over 40 genes have been identified as causative for isolated human CAKUT. However, many genes remain unknown, and the prioritization of potential CAKUT candidate genes is challenging. To develop an independent approach to prioritize CAKUT candidate genes, we hypothesized that monogenic CAKUT genes are most likely co-expressed along a temporal axis during kidney development and that genes with coinciding high expression may represent strong novel CAKUT candidate genes.

METHODS

We analyzed single-cell mRNA (sc-mRNA) transcriptomics data of human fetal kidney for temporal sc-mRNA co-expression of 40 known CAKUT genes. A maximum of high expression in consecutive timepoints of kidney development was found for four of the 40 genes (EYA1, SIX1, SIX2, and ITGA8) in nephron progenitor cells a, b, c, d (NPCa-d). We concluded that NPCa-d are relevant for CAKUT pathogenesis and intersected two lists of CAKUT candidate genes resulting from unbiased whole-exome sequencing (WES) with the 100 highest expressed genes in NPCa-d.

RESULTS

Intersection of the 100 highest expressed genes in NPCa-d with WES-derived CAKUT candidate genes identified an overlap with the candidate genes KIF19, TRIM36, USP35, CHTF18, in each of which a biallelic variant was detected in different families with CAKUT.

CONCLUSION

Sc-mRNA expression data of human fetal kidney can be utilized to prioritize WES-derived CAKUT candidate genes. KIF19, TRIM36, USP35, and CHTF18 may represent strong novel candidate genes for CAKUT.

scRNA-Seq of Cultured Human Amniotic Fluid from Fetuses with Spina Bifida Reveals the Origin and Heterogeneity of the Cellular Content

Amniotic fluid has been proposed as an easily available source of cells for numerous applications in regenerative medicine and tissue engineering. The use of amniotic fluid cells in biomedical applications necessitates their unequivocal characterization; however, the exact cellular composition of amniotic fluid and the precise tissue origins of these cells remain largely unclear. Using cells cultured from the human amniotic fluid of fetuses with spina bifida aperta and of a healthy fetus, we performed single-cell RNA sequencing to characterize the tissue origin and marker expression of cultured amniotic fluid cells at the single-cell level. Our analysis revealed nine different cell types of stromal, epithelial and immune cell phenotypes, and from various fetal tissue origins, demonstrating the heterogeneity of the cultured amniotic fluid cell population at a single-cell resolution. It also identified cell types of neural origin in amniotic fluid from fetuses with spina bifida aperta. Our data provide a comprehensive list of markers for the characterization of the various progenitor and terminally differentiated cell types in cultured amniotic fluid. This study highlights the relevance of single-cell analysis approaches for the characterization of amniotic fluid cells in order to harness their full potential in biomedical research and clinical applications.

Transfer learning enables predictions in network biology

Mapping gene networks requires large amounts of transcriptomic data to learn the connections between genes, which impedes discoveries in settings with limited data, including rare diseases and diseases affecting clinically inaccessible tissues. Recently, transfer learning has revolutionized fields such as natural language understanding1,2 and computer vision3 by leveraging deep learning models pretrained on large-scale general datasets that can then be fine-tuned towards a vast array of downstream tasks with limited task-specific data. Here, we developed a context-aware, attention-based deep learning model, Geneformer, pretrained on a large-scale corpus of about 30 million single-cell transcriptomes to enable context-specific predictions in settings with limited data in network biology. During pretraining, Geneformer gained a fundamental understanding of network dynamics, encoding network hierarchy in the attention weights of the model in a completely self-supervised manner. Fine-tuning towards a diverse panel of downstream tasks relevant to chromatin and network dynamics using limited task-specific data demonstrated that Geneformer consistently boosted predictive accuracy. Applied to disease modelling with limited patient data, Geneformer identified candidate therapeutic targets for cardiomyopathy. Overall, Geneformer represents a pretrained deep learning model from which fine-tuning towards a broad range of downstream applications can be pursued to accelerate discovery of key network regulators and candidate therapeutic targets.

Comparative single-cell analyses identify shared and divergent features of human and mouse kidney development

Mammalian kidneys maintain fluid homeostasis through the cellular activity of nephrons and the conjoined collecting system. Each epithelial network originates from distinct progenitor cell populations that reciprocally interact during development. To extend our understanding of human and mouse kidney development, we profiled chromatin organization (ATAC-seq) and gene expression (RNA-seq) in developing human and mouse kidneys. Data were analyzed at a species level and then integrated into a common, cross-species multimodal data set. Comparative analysis of cell types and developmental trajectories identified conserved and divergent features of chromatin organization and linked gene activity, revealing species- and cell-type specific regulatory programs. Identification of human-specific enhancer regions linked through GWAS studies to kidney disease highlights the potential of developmental modeling to provide clinical insight.

Podocyte-Parietal Epithelial Cell Interdependence in Glomerular Development and Disease

Podocytes and parietal epithelial cells (PECs) are among the few principal cell types within the kidney glomerulus, the former serving as a crucial constituent of the kidney filtration barrier and the latter representing a supporting epithelial layer that adorns the inner wall of Bowman's capsule. Podocytes and PECs share a circumscript developmental lineage that only begins to diverge during the S-shaped body stage of nephron formation-occurring immediately before the emergence of the fully mature nephron. These two cell types, therefore, share a highly conserved gene expression program, evidenced by recently discovered intermediate cell types occupying a distinct spatiotemporal gene expression zone between podocytes and PECs. In addition to their homeostatic functions, podocytes and PECs also have roles in kidney pathogenesis. Rapid podocyte loss in diseases, such as rapidly progressive GN and collapsing and cellular subtypes of FSGS, is closely allied with PEC proliferation and migration toward the capillary tuft, resulting in the formation of crescents and pseudocrescents. PECs are thought to contribute to disease progression and severity, and the interdependence between these two cell types during development and in various manifestations of kidney pathology is the primary focus of this review.

Growth and differentiation of human induced pluripotent stem cell (hiPSC)-derived kidney organoids using fully synthetic peptide hydrogels

Human induced pluripotent stem cell (hiPSC)-derived kidney organoids have prospective applications ranging from basic disease modelling to personalised medicine. However, there remains a necessity to refine the biophysical and biochemical parameters that govern kidney organoid formation. Differentiation within fully-controllable and physiologically relevant 3D growth environments will be critical to improving organoid reproducibility and maturation. Here, we matured hiPSC-derived kidney organoids within fully synthetic self-assembling peptide hydrogels (SAPHs) of variable stiffness (storage modulus, G'). The resulting organoids contained complex structures comparable to those differentiated within the animal-derived matrix, Matrigel. Single-cell RNA sequencing (scRNA-seq) was then used to compare organoids matured within SAPHs to those grown within Matrigel or at the air-liquid interface. A total of 13,179 cells were analysed, revealing 14 distinct clusters. Organoid compositional analysis revealed a larger proportion of nephron cell types within Transwell-derived organoids, while SAPH-derived organoids were enriched for stromal-associated cell populations. Notably, differentiation within a higher G' SAPH generated podocytes with more mature gene expression profiles. Additionally, maturation within a 3D microenvironment significantly reduced the derivation of off-target cell types, which are a known limitation of current kidney organoid protocols. This work demonstrates the utility of synthetic peptide-based hydrogels with a defined stiffness, as a minimally complex microenvironment for the selected differentiation of kidney organoids.

Vascular endothelial cell development and diversity

Vascular endothelial cells form the inner layer of blood vessels where they have a key role in the development and maintenance of the functional circulatory system and provide paracrine support to surrounding non-vascular cells. Technical advances in the past 5 years in single-cell genomics and in in vivo genetic labelling have facilitated greater insights into endothelial cell development, plasticity and heterogeneity. These advances have also contributed to a new understanding of the timing of endothelial cell subtype differentiation and its relationship to the cell cycle. Identification of novel tissue-specific gene expression patterns in endothelial cells has led to the discovery of crucial signalling pathways and new interactions with other cell types that have key roles in both tissue maintenance and disease pathology. In this Review, we describe the latest findings in vascular endothelial cell development and diversity, which are often supported by large-scale, single-cell studies, and discuss the implications of these findings for vascular medicine. In addition, we highlight how techniques such as single-cell multimodal omics, which have become increasingly sophisticated over the past 2 years, are being utilized to study normal vascular physiology as well as functional perturbations in disease.

SGLT2 inhibitors mitigate kidney tubular metabolic and mTORC1 perturbations in youth-onset type 2 diabetes

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.

Modeling Podocyte Ontogeny and Podocytopathies with the Zebrafish

Podocytes are exquisitely fashioned kidney cells that serve an essential role in the process of blood filtration. Congenital malformation or damage to podocytes has dire consequences and initiates a cascade of pathological changes leading to renal disease states known as podocytopathies. In addition, animal models have been integral to discovering the molecular pathways that direct the development of podocytes. In this review, we explore how researchers have used the zebrafish to illuminate new insights about the processes of podocyte ontogeny, model podocytopathies, and create opportunities to discover future therapies.

Genetic Variants in ARHGEF6 Cause Congenital Anomalies of the Kidneys and Urinary Tract in Humans, Mice, and Frogs

BACKGROUND

About 40 disease genes have been described to date for isolated CAKUT, the most common cause of childhood CKD. However, these genes account for only 20% of cases. ARHGEF6, a guanine nucleotide exchange factor that is implicated in biologic processes such as cell migration and focal adhesion, acts downstream of integrin-linked kinase (ILK) and parvin proteins. A genetic variant of ILK that causes murine renal agenesis abrogates the interaction of ILK with a murine focal adhesion protein encoded by Parva , leading to CAKUT in mice with this variant.

METHODS

To identify novel genes that, when mutated, result in CAKUT, we performed exome sequencing in an international cohort of 1265 families with CAKUT. We also assessed the effects in vitro of wild-type and mutant ARHGEF6 proteins, and the effects of Arhgef6 deficiency in mouse and frog models.

RESULTS

We detected six different hemizygous variants in the gene ARHGEF6 (which is located on the X chromosome in humans) in eight individuals from six families with CAKUT. In kidney cells, overexpression of wild-type ARHGEF6 -but not proband-derived mutant ARHGEF6 -increased active levels of CDC42/RAC1, induced lamellipodia formation, and stimulated PARVA-dependent cell spreading. ARHGEF6-mutant proteins showed loss of interaction with PARVA. Three-dimensional Madin-Darby canine kidney cell cultures expressing ARHGEF6-mutant proteins exhibited reduced lumen formation and polarity defects. Arhgef6 deficiency in mouse and frog models recapitulated features of human CAKUT.

CONCLUSIONS

Deleterious variants in ARHGEF6 may cause dysregulation of integrin-parvin-RAC1/CDC42 signaling, thereby leading to X-linked CAKUT.

Single-cell transcriptomics reveals a mechanosensitive injury signaling pathway in early diabetic nephropathy

BACKGROUND

Diabetic nephropathy (DN) is the leading cause of end-stage renal disease, and histopathologic glomerular lesions are among the earliest structural alterations of DN. However, the signaling pathways that initiate these glomerular alterations are incompletely understood.

METHODS

To delineate the cellular and molecular basis for DN initiation, we performed single-cell and bulk RNA sequencing of renal cells from type 2 diabetes mice (BTBR ob/ob) at the early stage of DN.

RESULTS

Analysis of differentially expressed genes revealed glucose-independent responses in glomerular cell types. The gene regulatory network upstream of glomerular cell programs suggested the activation of mechanosensitive transcriptional pathway MRTF-SRF predominantly taking place in mesangial cells. Importantly, activation of MRTF-SRF transcriptional pathway was also identified in DN glomeruli in independent patient cohort datasets. Furthermore, ex vivo kidney perfusion suggested that the regulation of MRTF-SRF is a common mechanism in response to glomerular hyperfiltration.

CONCLUSIONS

Overall, our study presents a comprehensive single-cell transcriptomic landscape of early DN, highlighting mechanosensitive signaling pathways as novel targets of diabetic glomerulopathy.

A comparative study of cellular diversity between the Xenopus pronephric and mouse metanephric nephron

The kidney is an essential organ that ensures bodily fluid homeostasis and removes soluble waste products from the organism. Nephrons, the functional units of the kidney, comprise a blood filter, the glomerulus or glomus, and an epithelial tubule that processes the filtrate from the blood or coelom and selectively reabsorbs solutes, such as sugars, proteins, ions, and water, leaving waste products to be eliminated in the urine. Genes coding for transporters are segmentally expressed, enabling the nephron to sequentially process the filtrate. The Xenopus embryonic kidney, the pronephros, which consists of a single large nephron, has served as a valuable model to identify genes involved in nephron formation and patterning. Therefore, the developmental patterning program that generates these segments is of great interest. Prior work has defined the gene expression profiles of Xenopus nephron segments via in situ hybridization strategies, but a comprehensive understanding of the cellular makeup of the pronephric kidney remains incomplete. Here, we carried out single-cell mRNA sequencing of the functional Xenopus pronephric nephron and evaluated its cellular composition through comparative analyses with previous Xenopus studies and single-cell mRNA sequencing of the adult mouse kidney. This study reconstructs the cellular makeup of the pronephric kidney and identifies conserved cells, segments, and associated gene expression profiles. Thus, our data highlight significant conservation in podocytes, proximal and distal tubule cells, and divergence in cellular composition underlying the capacity of each nephron to remove wastes in the form of urine, while emphasizing the Xenopus pronephros as a model for physiology and disease.

A single-cell comparison of adult and fetal human epicardium defines the age-associated changes in epicardial activity

Re-activating quiescent adult epicardium represents a potential therapeutic approach for human cardiac regeneration. However, the exact molecular differences between inactive adult and active fetal epicardium are not known. In this study, we combined fetal and adult human hearts using single-cell and single-nuclei RNA sequencing and compared epicardial cells from both stages. We found that a migratory fibroblast-like epicardial population only in the fetal heart and fetal epicardium expressed angiogenic gene programs, whereas the adult epicardium was solely mesothelial and immune responsive. Furthermore, we predicted that adult hearts may still receive fetal epicardial paracrine communication, including WNT signaling with endocardium, reinforcing the validity of regenerative strategies that administer or reactivate epicardial cells in situ. Finally, we explained graft efficacy of our human embryonic stem-cell-derived epicardium model by noting its similarity to human fetal epicardium. Overall, our study defines epicardial programs of regenerative angiogenesis absent in adult hearts, contextualizes animal studies and defines epicardial states required for effective human heart regeneration.

Single-cell multiomics reveals the complexity of TGFβ signalling to chromatin in iPSC-derived kidney organoids

TGFβ1 plays a regulatory role in the determination of renal cell fate and the progression of renal fibrosis. Here we show an association between SMAD3 and the histone methyltransferase, EZH2, during cell differentiation; ChIP-seq revealed that SMAD3 and EZH2 co-occupy the genome in iPSCs and in iPSC-derived nephron progenitors. Through integration of single cell gene expression and epigenome profiling, we identified de novo ACTA2^+ve/POSTN^+ve myofibroblasts in kidney organoids treated with TGFβ1, characterised by increased SMAD3-dependent cis chromatin accessibility and gene expression associated with fibroblast activation. We have identified fibrosis-associated regulons characterised by enrichment of SMAD3, AP1, the ETS family of transcription factors, and NUAK1, CREB3L1, and RARG, corresponding to enriched motifs at accessible loci identified by scATACseq. Treatment with the EZH2 specific inhibitor GSK343, blocked SMAD3-dependent cis co-accessibility and inhibited myofibroblast activation. This mechanism, through which TGFβ signals directly to chromatin, represents a critical determinant of fibrotic, differentiated states.

Principles of human and mouse nephron development

The mechanisms underlying kidney development in mice and humans is an area of intense study. Insights into kidney organogenesis have the potential to guide our understanding of the origin of congenital anomalies and enable the assembly of genetic diagnostic tools. A number of studies have delineated signalling nodes that regulate positional identities and cell fates of nephron progenitor and precursor cells, whereas cross-species comparisons have markedly enhanced our understanding of conserved and divergent features of mammalian kidney organogenesis. Greater insights into the complex cellular movements that occur as the proximal-distal axis is established have challenged our understanding of nephron patterning and provided important clues to the elaborate developmental context in which human kidney diseases can arise. Studies of kidney development in vivo have also facilitated efforts to recapitulate nephrogenesis in kidney organoids in vitro, by providing a detailed blueprint of signalling events, cell movements and patterning mechanisms that are required for the formation of correctly patterned nephrons and maturation of physiologically functional apparatus that are responsible for maintaining human health.

Assessing kidney development and disease using kidney organoids and CRISPR engineering

The differentiation of human pluripotent stem cells (hPSCs) towards organoids is one of the biggest scientific advances in regenerative medicine. Kidney organoids have not only laid the groundwork for various organ-like tissue systems but also provided insights into kidney embryonic development. Thus, several protocols for the differentiation of renal progenitors or mature cell types have been established. Insights into the interplay of developmental pathways in nephrogenesis and determination of different cell fates have enabled the in vitro recapitulation of nephrogenesis. Here we first provide an overview of kidney morphogenesis and patterning in the mouse model in order to dissect signalling pathways that are key to define culture conditions sustaining renal differentiation from hPSCs. Secondly, we also highlight how genome editing approaches have provided insights on the specific role of different genes and molecular pathways during renal differentiation from hPSCs. Based on this knowledge we further review how CRISPR/Cas9 technology has enabled the recapitulation and correction of cellular phenotypes associated with human renal disease. Last, we also revise how the field has positively benefited from emerging technologies as single cell RNA sequencing and discuss current limitations on kidney organoid technology that will take advantage from bioengineering solutions to help standardizing the use of this model systems to study kidney development and disease.

Single-cell and spatial mapping Identify cell types and signaling Networks in the human ureter

Tissue engineering offers a promising treatment strategy for ureteral strictures, but its success requires an in-depth understanding of the architecture, cellular heterogeneity, and signaling pathways underlying tissue regeneration. Here, we define and spatially map cell populations within the human ureter using single-cell RNA sequencing, spatial gene expression, and immunofluorescence approaches. We focus on the stromal and urothelial cell populations to enumerate the distinct cell types composing the human ureter and infer potential cell-cell communication networks underpinning the bi-directional crosstalk between these compartments. Furthermore, we analyze and experimentally validate the importance of the sonic hedgehog (SHH) signaling pathway in adult progenitor cell maintenance. The SHH-expressing basal cells support organoid generation in vitro and accurately predict the differentiation trajectory from basal progenitor cells to terminally differentiated umbrella cells. Our results highlight the essential processes involved in adult ureter tissue homeostasis and provide a blueprint for guiding ureter tissue engineering.

Transplanted organoids empower human preclinical assessment of drug candidate for the clinic

Pharmacodynamic (PD) studies are an essential component of preclinical drug discovery. Current approaches for PD studies, including the analysis of novel kidney disease targeting therapeutic agents, are limited to animal models with unclear translatability to the human condition. To address this challenge, we developed a novel approach for PD studies using transplanted, perfused human kidney organoids. We performed pharmacokinetic (PK) studies with GFB-887, an investigational new drug now in phase 2 trials. Orally dosed GFB-887 to athymic rats that had undergone organoid transplantation resulted in measurable drug exposure in transplanted organoids. We established the efficacy of orally dosed GFB-887 in PD studies, where quantitative analysis showed significant protection of kidney filter cells in human organoids and endogenous rat host kidneys. This widely applicable approach demonstrates feasibility of using transplanted human organoids in preclinical PD studies with an investigational new drug, empowering organoids to revolutionize drug discovery.

A scalable organoid model of human autosomal dominant polycystic kidney disease for disease mechanism and drug discovery

Human pluripotent stem-cell-derived organoids are models for human development and disease. We report a modified human kidney organoid system that generates thousands of similar organoids, each consisting of 1-2 nephron-like structures. Single-cell transcriptomic profiling and immunofluorescence validation highlighted patterned nephron-like structures utilizing similar pathways, with distinct morphogenesis, to human nephrogenesis. To examine this platform for therapeutic screening, the polycystic kidney disease genes PKD1 and PKD2 were inactivated by gene editing. PKD1 and PKD2 mutant models exhibited efficient and reproducible cyst formation. Cystic outgrowths could be propagated for months to centimeter-sized cysts. To shed new light on cystogenesis, 247 protein kinase inhibitors (PKIs) were screened in a live imaging assay identifying compounds blocking cyst formation but not overall organoid growth. Scaling and further development of the organoid platform will enable a broader capability for kidney disease modeling and high-throughput drug screens.

Single-Cell Chromatin and Gene-Regulatory Dynamics of Mouse Nephron Progenitors

BACKGROUND

We reasoned that unraveling the dynamic changes in accessibility of genomic regulatory elements and gene expression at single-cell resolution will inform the basic mechanisms of nephrogenesis.

METHODS

We performed single-cell ATAC-seq and RNA-seq both individually (singleomes; Six2GFP cells) and jointly in the same cells (multiomes; kidneys) to generate integrated chromatin and transcriptional maps in mouse embryonic and neonatal nephron progenitor cells.

RESULTS

We demonstrate that singleomes and multiomes are comparable in assigning most cell states, identification of new cell type markers, and defining the transcription factors driving cell identity. However, multiomes are more precise in defining the progenitor population. Multiomes identified a "pioneer" bHLH/Fox motif signature in nephron progenitor cells. Moreover, we identified a subset of Fox factors exhibiting high chromatin activity in podocytes. One of these Fox factors, Foxp1, is important for nephrogenesis. Key nephrogenic factors are distinguished by strong correlation between linked gene regulatory elements and gene expression.

CONCLUSION

Mapping the regulatory landscape at single-cell resolution informs the regulatory hierarchy of nephrogenesis. Paired single-cell epigenomes and transcriptomes of nephron progenitors should provide a foundation to understand prenatal programming, regeneration after injury, and ex vivo nephrogenesis.

Kidney Development: Recent Insights from Technological Advances

The kidney is a complex organ, and how it forms is a fascinating process. New technologies, such as single-cell transcriptomics, and enhanced imaging modalities are offering new approaches to understand the complex and intertwined processes during embryonic kidney development.

Identification, discrimination and heterogeneity of fibroblasts

Fibroblasts, the principal cell type of connective tissue, secrete extracellular matrix components during tissue development, homeostasis, repair and disease. Despite this crucial role, the identification and distinction of fibroblasts from other cell types are challenging and laden with caveats. Rapid progress in single-cell transcriptomics now yields detailed molecular portraits of fibroblasts and other cell types in our bodies, which complement and enrich classical histological and immunological descriptions, improve cell class definitions and guide further studies on the functional heterogeneity of cell subtypes and states, origins and fates in physiological and pathological processes. In this review, we summarize and discuss recent advances in the understanding of fibroblast identification and heterogeneity and how they discriminate from other cell types.

In this review, the authors look at how recent progress in single-cell transcriptomics complement and enrich the classical, largely morphological, portraits of fibroblasts. The detailed molecular information now available provides new insights into fibroblast identity, heterogeneity and function.

Integrating spatial transcriptomics with single-cell transcriptomics reveals a spatiotemporal gene landscape of the human developing kidney

BACKGROUND

Research on spatiotemporal gene landscape can provide insights into the spatial characteristics of human kidney development and facilitate kidney organoid cultivation. Here, we profiled the spatiotemporal gene programs of the human embryonic kidneys at 9 and 18 post-conception weeks (PCW) by integrating the application of microarray-based spatial transcriptomics and single-cell transcriptomics.

RESULTS

We mapped transcriptomic signatures of scRNA-seq cell types upon the 9 and 18 PCW kidney sections based on cell-type deconvolution and multimodal intersection analyses, depicting a spatial landscape of developing cell subpopulations. We established the gene characteristics in the medullary regions and revealed a strong mitochondrial oxidative phosphorylation and glycolysis activity in the deeper medullary region. We also built a regulatory network centered on GDNF-ETV4 for nephrogenic niche development based on the weighted gene co-expression network analysis and highlighted the key roles of Wnt, FGF, and JAG1-Notch2 signaling in maintaining renal branching morphogenesis.

CONCLUSIONS

Our findings obtained by this spatiotemporal gene program are expected to improve the current understanding of kidney development.

Enhanced metanephric specification to functional proximal tubule enables toxicity screening and infectious disease modelling in kidney organoids

While pluripotent stem cell-derived kidney organoids are now being used to model renal disease, the proximal nephron remains immature with limited evidence for key functional solute channels. This may reflect early mispatterning of the nephrogenic mesenchyme and/or insufficient maturation. Here we show that enhanced specification to metanephric nephron progenitors results in elongated and radially aligned proximalised nephrons with distinct S1 - S3 proximal tubule cell types. Such PT-enhanced organoids possess improved albumin and organic cation uptake, appropriate KIM-1 upregulation in response to cisplatin, and improved expression of SARS-CoV-2 entry factors resulting in increased viral replication. The striking proximo-distal orientation of nephrons resulted from localized WNT antagonism originating from the organoid stromal core. PT-enhanced organoids represent an improved model to study inherited and acquired proximal tubular disease as well as drug and viral responses.

Comparative whole-genome transcriptome analysis in renal cell populations reveals high tissue specificity of MAPK/ERK targets in embryonic kidney

Background

MAPK/ERK signaling is a well-known mediator of extracellular stimuli controlling intracellular responses to growth factors and mechanical cues. The critical requirement of MAPK/ERK signaling for embryonic stem cell maintenance is demonstrated, but specific functions in progenitor regulation during embryonic development, and in particular kidney development remain largely unexplored. We previously demonstrated MAPK/ERK signaling as a key regulator of kidney growth through branching morphogenesis and normal nephrogenesis where it also regulates progenitor expansion. Here, we performed RNA sequencing-based whole-genome expression analysis to identify transcriptional MAPK/ERK targets in two distinct renal populations: the ureteric bud epithelium and the nephron progenitors.

Results

Our analysis revealed a large number (5053) of differentially expressed genes (DEGs) in nephron progenitors and significantly less (1004) in ureteric bud epithelium, reflecting likely heterogenicity of cell types. The data analysis identified high tissue-specificity, as only a fraction (362) of MAPK/ERK targets are shared between the two tissues. Tissue-specific MAPK/ERK targets participate in the regulation of mitochondrial energy metabolism in nephron progenitors, which fail to maintain normal mitochondria numbers in the MAPK/ERK-deficient tissue. In the ureteric bud epithelium, a dramatic decline in progenitor-specific gene expression was detected with a simultaneous increase in differentiation-associated genes, which was not observed in nephron progenitors. Our experiments in the genetic model of MAPK/ERK deficiency provide evidence that MAPK/ERK signaling in the ureteric bud maintains epithelial cells in an undifferentiated state. Interestingly, the transcriptional targets shared between the two tissues studied are over-represented by histone genes, suggesting that MAPK/ERK signaling regulates cell cycle progression and stem cell maintenance through chromosome condensation and nucleosome assembly.

Conclusions

Using tissue-specific MAPK/ERK inactivation and RNA sequencing in combination with experimentation in embryonic kidneys, we demonstrate here that MAPK/ERK signaling maintains ureteric bud tip cells, suggesting a regulatory role in collecting duct progenitors. We additionally deliver new mechanistic information on how MAPK/ERK signaling regulates progenitor maintenance through its effects on chromatin accessibility and energy metabolism.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01309-z.

Single-Nucleus Transcriptional Profiling of Chronic Kidney Disease after Cisplatin Nephrotoxicity

Cisplatin induces both acute and chronic nephrotoxicity during chemotherapy in patients with cancer. Presented here is the first study of single-nucleus RNA sequencing (snRNA-seq) of cisplatin-induced nephrotoxicity. Repeated low-dose cisplatin treatment (RLDC) led to decreases in renal function and kidney weight in mice at 9 weeks. The kidneys of these mice showed tubular degeneration and dilation. snRNA-seq identified 16 cell types and 17 cell clusters in these kidneys. Cluster-by-cluster comparison demonstrated cell type-specific changes in gene expression and identified a unique proximal tubule (PT) injury/repair cluster that co-expressed the injury marker kidney injury molecule-1 (Kim1) and the proliferation marker Ki-67. Compared with control, post-RLDC kidneys had 424 differentially expressed genes in PT cells, including tubular transporters and cytochrome P450 enzymes involved in lipid metabolism. snRNA-seq also revealed transcriptional changes in potential PT injury markers (Krt222, Eda2r, Ltbp2, and Masp1) and repair marker (Bex4). RLDC induced inflammation and proinflammatory cytokines (RelB, TNF-α, Il7, Ccl2, and Cxcl2) and the expression of fibrosis markers (fibronectin, collagen I, connective tissue growth factor, vimentin, and α-smooth muscle actin). Together, these results provide new insights into RLDC-induced transcriptional changes at the single-cell level that may contribute to the development of chronic kidney problems in patients with cancer after cisplatin chemotherapy.

Mapping the human kidney using single-cell genomics

The field of single-cell genomics and spatial technologies is rapidly evolving and has already provided unprecedented insights into complex tissues. Major advances have been made in dissecting the cellular composition and spatiotemporal interactions that mediate developmental processes in the fetal kidney. Single-cell technologies have also provided detailed insights into the heterogeneity of cell types within the healthy adult and shed light on the complex cellular mechanisms that contribute to kidney disease. The in-depth characterization of specific cell types associated with acute kidney injury and glomerular diseases has potential for the development of prognostic biomarkers and new therapeutics. Analyses of pathway activity in clear-cell renal cell carcinoma can predict the sensitivity of tumour cells to specific inhibitors. The identification of the cell of origin of renal cell carcinoma and of new cell types within the tumour microenvironment also has implications for the development of targeted therapeutics. Similarly, single-cell sequencing has provided new insights into the mechanisms underlying kidney fibrosis, specifically our understanding of myofibroblast origins and the contribution of cell crosstalk within the fibrotic niche to disease progression. These and future studies will enable the creation of a map to aid our understanding of the cellular processes and interactions in the developing, healthy and diseased kidney.

DevKidCC allows for robust classification and direct comparisons of kidney organoid datasets

BACKGROUND

While single-cell transcriptional profiling has greatly increased our capacity to interrogate biology, accurate cell classification within and between datasets is a key challenge. This is particularly so in pluripotent stem cell-derived organoids which represent a model of a developmental system. Here, clustering algorithms and selected marker genes can fail to accurately classify cellular identity while variation in analyses makes it difficult to meaningfully compare datasets. Kidney organoids provide a valuable resource to understand kidney development and disease. However, direct comparison of relative cellular composition between protocols has proved challenging. Hence, an unbiased approach for classifying cell identity is required.

METHODS

The R package, scPred, was trained on multiple single cell RNA-seq datasets of human fetal kidney. A hierarchical model classified cellular subtypes into nephron, stroma and ureteric epithelial elements. This model, provided in the R package DevKidCC ( github.com/KidneyRegeneration/DevKidCC ), was then used to predict relative cell identity within published kidney organoid datasets generated using distinct cell lines and differentiation protocols, interrogating the impact of such variations. The package contains custom functions for the display of differential gene expression within cellular subtypes.

RESULTS

DevKidCC was used to directly compare between distinct kidney organoid protocols, identifying differences in relative proportions of cell types at all hierarchical levels of the model and highlighting variations in stromal and unassigned cell types, nephron progenitor prevalence and relative maturation of individual epithelial segments. Of note, DevKidCC was able to distinguish distal nephron from ureteric epithelium, cell types with overlapping profiles that have previously confounded analyses. When applied to a variation in protocol via the addition of retinoic acid, DevKidCC identified a consequential depletion of nephron progenitors.

CONCLUSIONS

The application of DevKidCC to kidney organoids reproducibly classifies component cellular identity within distinct single-cell datasets. The application of the tool is summarised in an interactive Shiny application, as are examples of the utility of in-built functions for data presentation. This tool will enable the consistent and rapid comparison of kidney organoid protocols, driving improvements in patterning to kidney endpoints and validating new approaches.

DevKidCC allows for robust classification and direct comparisons of kidney organoid datasets

Background

While single-cell transcriptional profiling has greatly increased our capacity to interrogate biology, accurate cell classification within and between datasets is a key challenge. This is particularly so in pluripotent stem cell-derived organoids which represent a model of a developmental system. Here, clustering algorithms and selected marker genes can fail to accurately classify cellular identity while variation in analyses makes it difficult to meaningfully compare datasets. Kidney organoids provide a valuable resource to understand kidney development and disease. However, direct comparison of relative cellular composition between protocols has proved challenging. Hence, an unbiased approach for classifying cell identity is required.

Methods

The R package, scPred, was trained on multiple single cell RNA-seq datasets of human fetal kidney. A hierarchical model classified cellular subtypes into nephron, stroma and ureteric epithelial elements. This model, provided in the R package DevKidCC (github.com/KidneyRegeneration/DevKidCC), was then used to predict relative cell identity within published kidney organoid datasets generated using distinct cell lines and differentiation protocols, interrogating the impact of such variations. The package contains custom functions for the display of differential gene expression within cellular subtypes.

Results

DevKidCC was used to directly compare between distinct kidney organoid protocols, identifying differences in relative proportions of cell types at all hierarchical levels of the model and highlighting variations in stromal and unassigned cell types, nephron progenitor prevalence and relative maturation of individual epithelial segments. Of note, DevKidCC was able to distinguish distal nephron from ureteric epithelium, cell types with overlapping profiles that have previously confounded analyses. When applied to a variation in protocol via the addition of retinoic acid, DevKidCC identified a consequential depletion of nephron progenitors.

Conclusions

The application of DevKidCC to kidney organoids reproducibly classifies component cellular identity within distinct single-cell datasets. The application of the tool is summarised in an interactive Shiny application, as are examples of the utility of in-built functions for data presentation. This tool will enable the consistent and rapid comparison of kidney organoid protocols, driving improvements in patterning to kidney endpoints and validating new approaches.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01023-z.

Alignment of single-cell trajectories by tuMap enables high-resolution quantitative comparison of cancer samples

Single-cell technologies allow characterization of cancer samples as continuous developmental trajectories. Yet, the obtained temporal resolution cannot be leveraged for a comparative analysis due to the large phenotypic heterogeneity existing between patients. Here, we present the tuMap algorithm that exploits high-dimensional single-cell data of cancer samples exhibiting an underlying developmental structure to align them with the healthy development, yielding the tuMap pseudotime axis that allows their systematic, meaningful comparison. We applied tuMap on single-cell mass cytometry data of acute lymphoblastic and myeloid leukemia to reveal associations between the tuMap pseudotime axis and clinics that outperform cellular assignment into developmental populations. Application of the tuMap algorithm on single-cell RNA sequencing data further identified gene signatures of stem cells residing at the very-early parts of the cancer trajectories. The quantitative framework provided by tuMap allows generation of metrics for cancer patients evaluation.

RNA profiles reveal signatures of future health and disease in pregnancy

Maternal morbidity and mortality continue to rise, and pre-eclampsia is a major driver of this burden1. Yet the ability to assess underlying pathophysiology before clinical presentation to enable identification of pregnancies at risk remains elusive. Here we demonstrate the ability of plasma cell-free RNA (cfRNA) to reveal patterns of normal pregnancy progression and determine the risk of developing pre-eclampsia months before clinical presentation. Our results centre on comprehensive transcriptome data from eight independent prospectively collected cohorts comprising 1,840 racially diverse pregnancies and retrospective analysis of 2,539 banked plasma samples. The pre-eclampsia data include 524 samples (72 cases and 452 non-cases) from two diverse independent cohorts collected 14.5 weeks (s.d., 4.5 weeks) before delivery. We show that cfRNA signatures from a single blood draw can track pregnancy progression at the placental, maternal and fetal levels and can robustly predict pre-eclampsia, with a sensitivity of 75% and a positive predictive value of 32.3% (s.d., 3%), which is superior to the state-of-the-art method2. cfRNA signatures of normal pregnancy progression and pre-eclampsia are independent of clinical factors, such as maternal age, body mass index and race, which cumulatively account for less than 1% of model variance. Further, the cfRNA signature for pre-eclampsia contains gene features linked to biological processes implicated in the underlying pathophysiology of pre-eclampsia.

How Times Have Changed! A Cornucopia of Antigens for Membranous Nephropathy

The identification of the major target antigen phospholipase A2 receptor (PLA2R) in the majority of primary (idiopathic) cases of membranous nephropathy (MN) has been followed by the rapid identification of numerous minor antigens that appear to define phenotypically distinct forms of disease. This article serves to review all the known antigens that have been shown to localize to subepithelial deposits in MN, as well as the distinctive characteristics associated with each subtype of MN. We will also shed light on the novel proteomic approaches that have allowed identification of the most recent antigens. The paradigm of an antigen normally expressed on the podocyte cell surface leading to in-situ immune complex formation, complement activation, and subsequent podocyte injury will be discussed and challenged in light of the current repertoire of multiple MN antigens. Since disease phenotypes associated with each individual target antigens can often blur the distinction between primary and secondary disease, we encourage the use of antigen-based classification of membranous nephropathy.

Single-cell RNA sequencing reveals differential cell cycle activity in key cell populations during nephrogenesis

The kidney is a complex organ composed of more than 30 terminally differentiated cell types that all are required to perform its numerous homeostatic functions. Defects in kidney development are a significant cause of chronic kidney disease in children, which can lead to kidney failure that can only be treated by transplant or dialysis. A better understanding of molecular mechanisms that drive kidney development is important for designing strategies to enhance renal repair and regeneration. In this study, we profiled gene expression in the developing mouse kidney at embryonic day 14.5 at single-cell resolution. Consistent with previous studies, clusters with distinct transcriptional signatures clearly identify major compartments and cell types of the developing kidney. Cell cycle activity distinguishes between the "primed" and "self-renewing" sub-populations of nephron progenitors, with increased expression of the cell cycle-related genes Birc5, Cdca3, Smc2 and Smc4 in "primed" nephron progenitors. In addition, augmented expression of cell cycle related genes Birc5, Cks2, Ccnb1, Ccnd1 and Tuba1a/b was detected in immature distal tubules, suggesting cell cycle regulation may be required for early events of nephron patterning and tubular fusion between the distal nephron and collecting duct epithelia.

The origin and role of the renal stroma

The postnatal kidney is predominantly composed of nephron epithelia with the interstitial components representing a small proportion of the final organ, except in the diseased state. This is in stark contrast to the developing organ, which arises from the mesoderm and comprises an expansive stromal population with distinct regional gene expression. In many organs, the identity and ultimate function of an epithelium is tightly regulated by the surrounding stroma during development. However, although the presence of a renal stromal stem cell population has been demonstrated, the focus has been on understanding the process of nephrogenesis whereas the role of distinct stromal components during kidney morphogenesis is less clear. In this Review, we consider what is known about the role of the stroma of the developing kidney in nephrogenesis, where these cells come from as well as their heterogeneity, and reflect on how this information may improve human kidney organoid models.