Single-Cell RNA Sequencing Reveals mRNA Splice Isoform Switching during Kidney Development

Novel insights into kidney disease: the scRNA-seq and spatial transcriptomics approaches: a literature review

Over the past decade, single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics (ST) have revolutionized biomedical research, particularly in understanding cellular heterogeneity in kidney diseases. This review summarizes the application and development of scRNA-seq combined with ST in the context of kidney disease. By dissecting cellular heterogeneity at an unprecedented resolution, these advanced techniques have identified novel cell subpopulations and their dynamic interactions within the renal microenvironment. The integration of scRNA-seq with ST has been instrumental in elucidating the cellular and molecular mechanisms underlying kidney development, homeostasis, and disease progression. This approach has not only identified key cellular players in renal pathophysiology but also revealed the spatial organization of cells within the kidney, which is crucial for understanding their functional specialization. This paper highlights the transformative impact of these techniques on renal research that have paved the way for targeted therapeutic interventions and personalized medicine in the management of kidney disease.

The crosstalk between alternative splicing and circular RNA in cancer: pathogenic insights and therapeutic implications

RNA splicing is a fundamental step of gene expression. While constitutive splicing removes introns and joins exons unbiasedly, alternative splicing (AS) selectively determines the assembly of exons and introns to generate RNA variants corresponding to the same transcript. The biogenesis of circular RNAs (circRNAs) is inextricably associated with AS. Back-splicing, the biogenic process of circRNA, is a special form of AS. In cancer, both AS and circRNA deviate from the original track. In the present review, we delve into the intricate interplay between AS and circRNAs in the context of cancer. The relationship between AS and circRNAs is intricate, where AS modulates the biogenesis of circRNAs and circRNAs in return regulate AS events. Beyond that, epigenetic and posttranscriptional modifications concurrently regulate AS and circRNAs. On the basis of this modality, we summarize current knowledge on how splicing factors and other RNA binding proteins regulate circRNA biogenesis, and how circRNAs interact with splicing factors to influence AS events. Specifically, the feedback loop regulation between circRNAs and AS events contributes greatly to oncogenesis and cancer progression. In summary, resolving the crosstalk between AS and circRNA will not only provide better insight into cancer biology but also provoke novel strategies to combat cancer.

Renal Endothelial Single-Cell Transcriptomics Reveals Spatiotemporal Regulation and Divergent Roles of Differential Gene Transcription and Alternative Splicing in Murine Diabetic Nephropathy

Endothelial cell (EC) injury is a crucial contributor to the progression of diabetic kidney disease (DKD), but the specific EC populations and mechanisms involved remain elusive. Kidney ECs (n = 5464) were collected at three timepoints from diabetic BTBRob/ob mice and non-diabetic littermates. Their heterogeneity, transcriptional changes, and alternative splicing during DKD progression were mapped using SmartSeq2 single-cell RNA sequencing (scRNAseq) and elucidated through pathway, network, and gene ontology enrichment analyses. We identified 13 distinct transcriptional EC phenotypes corresponding to different kidney vessel subtypes, confirmed through in situ hybridization and immunofluorescence. EC subtypes along nephrons displayed extensive zonation related to their functions. Differential gene expression analyses in peritubular and glomerular ECs in DKD underlined the regulation of DKD-relevant pathways including EIF2 signaling, oxidative phosphorylation, and IGF1 signaling. Importantly, this revealed the differential alteration of these pathways between the two EC subtypes and changes during disease progression. Furthermore, glomerular and peritubular ECs also displayed aberrant and dynamic alterations in alternative splicing (AS), which is strongly associated with DNA repair. Strikingly, genes displaying differential transcription or alternative splicing participate in divergent biological processes. Our study reveals the spatiotemporal regulation of gene transcription and AS linked to DKD progression, providing insight into pathomechanisms and clues to novel therapeutic targets for DKD treatment.

Netrin 1 directs vascular patterning and maturity in the developing kidney

The intricate vascular system of the kidneys supports body fluid and organ homeostasis. However, little is known about how vascular architecture is established during kidney development. More specifically, how signals from the kidney influence vessel maturity and patterning remains poorly understood. Netrin 1 (Ntn1) is a secreted ligand that is crucial for vessel and neuronal guidance. Here, we demonstrate that Ntn1 is expressed by Foxd1+ stromal progenitors in the developing mouse kidney and conditional deletion (Foxd1GC/+;Ntn1fl/fl) results in hypoplastic kidneys with extended nephrogenesis. Wholemount 3D analyses additionally revealed the loss of a predictable vascular pattern in Foxd1GC/+;Ntn1fl/fl kidneys. As vascular patterning has been linked to vessel maturity, we investigated arterialization. Quantification of the CD31+ endothelium at E15.5 revealed no differences in metrics such as the number of branches or branch points, whereas the arterial vascular smooth muscle metrics were significantly reduced at both E15.5 and P0. In support of our observed phenotypes, whole kidney RNA-seq revealed disruptions to genes and programs associated with stromal cells, vasculature and differentiating nephrons. Together, our findings highlight the significance of Ntn1 to proper vascularization and kidney development.

Characterization of Continuous Transcriptional Heterogeneity in High-Risk Blastemal-Type Wilms' Tumors Using Unsupervised Machine Learning

Wilms' tumors are pediatric malignancies that are thought to arise from faulty kidney development. They contain a wide range of poorly differentiated cell states resembling various distorted developmental stages of the fetal kidney, and as a result, differ between patients in a continuous manner that is not well understood. Here, we used three computational approaches to characterize this continuous heterogeneity in high-risk blastemal-type Wilms' tumors. Using Pareto task inference, we show that the tumors form a triangle-shaped continuum in latent space that is bounded by three tumor archetypes with "stromal", "blastemal", and "epithelial" characteristics, which resemble the un-induced mesenchyme, the cap mesenchyme, and early epithelial structures of the fetal kidney. By fitting a generative probabilistic "grade of membership" model, we show that each tumor can be represented as a unique mixture of three hidden "topics" with blastemal, stromal, and epithelial characteristics. Likewise, cellular deconvolution allows us to represent each tumor in the continuum as a unique combination of fetal kidney-like cell states. These results highlight the relationship between Wilms' tumors and kidney development, and we anticipate that they will pave the way for more quantitative strategies for tumor stratification and classification.

Characterization of alternative mRNA splicing in cultured cell populations representing progressive stages of human fetal kidney development

Nephrons are the functional units of the kidney. During kidney development, cells from the cap mesenchyme-a transient kidney-specific progenitor state-undergo a mesenchymal to epithelial transition (MET) and subsequently differentiate into the various epithelial cell types that create the tubular structures of the nephron. Faults in this transition can lead to a pediatric malignancy of the kidney called Wilms' tumor that mimics normal kidney development. While human kidney development has been characterized at the gene expression level, a comprehensive characterization of alternative splicing is lacking. Therefore, in this study, we performed RNA sequencing on cell populations representing early, intermediate, and late developmental stages of the human fetal kidney, as well as three blastemal-predominant Wilms' tumor patient-derived xenografts. Using this newly generated RNAseq data, we identified a set of transcripts that are alternatively spliced between the different developmental stages. Moreover, we found that cells from the earliest developmental stage have a mesenchymal splice-isoform profile that is similar to that of blastemal-predominant Wilms' tumor xenografts. RNA binding motif enrichment analysis suggests that the mRNA binding proteins ESRP1, ESRP2, RBFOX2, and QKI regulate alternative mRNA splicing during human kidney development. These findings illuminate new molecular mechanisms involved in human kidney development and pediatric kidney cancer.

Underlying mechanisms of epithelial splicing regulatory proteins in cancer progression

Cancer is the second-leading disease-related cause of global mortality after cardiovascular disease. Despite significant advances in cancer therapeutic strategies, cancer remains one of the major obstacles to human life extension. Cancer pathogenesis is extremely complicated and not fully understood. Epithelial splicing regulatory proteins (ESRPs), including ESRP1 and ESRP2, belong to the heterogeneous nuclear ribonucleoprotein family of RNA-binding proteins and are crucial regulators of the alternative splicing of messenger RNAs (mRNAs). The expression and activity of ESRPs are modulated by various mechanisms, including post-translational modifications and non-coding RNAs. Although a growing body of evidence suggests that ESRP dysregulation is closely associated with cancer progression, the detailed mechanisms remain inconclusive. In this review, we summarize recent findings on the structures, functions, and regulatory mechanisms of ESRPs and focus on their underlying mechanisms in cancer progression. We also highlight the clinical implications of ESRPs as prognostic biomarkers and therapeutic targets in cancer treatment. The information reviewed herein could be extremely beneficial to the development of individualized therapeutic strategies for cancer patients.

Kidney single-cell transcriptome profile reveals distinct response of proximal tubule cells to SGLT2i and ARB treatment in diabetic mice

Angiotensin receptor blockers (ARBs) and sodium-glucose cotransporter 2 inhibitors (SGLT2i) have been used as the standard therapy for patients with diabetic kidney disease (DKD). However, how these two drugs possess additive renoprotective effects remains unclear. Here, we conducted single-cell RNA sequencing to profile the kidney cell transcriptome of db/db mice treated with vehicle, ARBs, SGLT2i, or ARBs plus SGLT2i, using db/m mice as control. We identified 10 distinct clusters of kidney cells with predominant proximal tubular (PT) cells. We found that ARBs had more anti-inflammatory and anti-fibrotic effects, while SGLT2i affected more mitochondrial function in PT. We also identified a new PT subcluster, was increased in DKD, but reversed by the treatments. This new subcluster was also confirmed by immunostaining of mouse and human kidneys with DKD. Together, our study reveals kidney cell-specific gene signatures in response to ARBs and SGLT2i and identifies a new PT subcluster, which provides new insight into the pathogenesis of DKD.

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.

The epithelial splicing regulator ESRP2 is epigenetically repressed by DNA hypermethylation in Wilms tumour and acts as a tumour suppressor

Wilms tumour (WT), an embryonal kidney cancer, has been extensively characterised for genetic and epigenetic alterations, but a proportion of WTs still lack identifiable abnormalities. To uncover DNA methylation changes critical for WT pathogenesis, we compared the epigenome of foetal kidney with two WT cell lines, filtering our results to remove common cancer-associated epigenetic changes and to enrich for genes involved in early kidney development. This identified four hypermethylated genes, of which ESRP2 (epithelial splicing regulatory protein 2) was the most promising for further study. ESRP2 was commonly repressed by DNA methylation in WT, and this occurred early in WT development (in nephrogenic rests). ESRP2 expression was reactivated by DNA methyltransferase inhibition in WT cell lines. When ESRP2 was overexpressed in WT cell lines, it inhibited cellular proliferation in vitro, and in vivo it suppressed tumour growth of orthotopic xenografts in nude mice. RNA-seq of the ESRP2-expressing WT cell lines identified several novel splicing targets. We propose a model in which epigenetic inactivation of ESRP2 disrupts the mesenchymal to epithelial transition in early kidney development to generate WT.

Effects of alternative splicing events and transcriptome changes on kidney stone formation

During the development of urinary stone disease, the formation of tiny crystals that adhere to the renal tubular epithelium induces epithelial cell damage. This damage and repair of the epithelium is associated with the establishment of more crystal adhesion sites, which in turn stimulates further crystal adhesion and, eventually, stone formation. Deposited crystals typically cause changes in epithelial cell gene expression, such as transcriptome changes and alternative splicing events. Although considered important for regulating gene expression, alternative splicing has not been reported in studies related to kidney stones. To date, whether alternative splicing events are involved in the regulation of stone formation and whether crystallographic cell interactions are regulated by alternative splicing at the transcriptional level have remained unknown. Therefore, we conducted RNA sequencing and alternative splicing-related bioassays by modeling the in vitro stone environment. Many alternative splicing events were associated with crystallographic cell interactions. Moreover, these events regulated transcription and significantly affected the capacity of crystals to adhere to renal tubular epithelial cells and regulate apoptosis.

Probing Isoform Switching Events in Various Cancer Types: Lessons From Pan-Cancer Studies

Alternative splicing is an essential regulatory mechanism for gene expression in mammalian cells contributing to protein, cellular, and species diversity. In cancer, alternative splicing is frequently disturbed, leading to changes in the expression of alternatively spliced protein isoforms. Advances in sequencing technologies and analysis methods led to new insights into the extent and functional impact of disturbed alternative splicing events. In this review, we give a brief overview of the molecular mechanisms driving alternative splicing, highlight the function of alternative splicing in healthy tissues and describe how alternative splicing is disrupted in cancer. We summarize current available computational tools for analyzing differential transcript usage, isoform switching events, and the pathogenic impact of cancer-specific splicing events. Finally, the strategies of three recent pan-cancer studies on isoform switching events are compared. Their methodological similarities and discrepancies are highlighted and lessons learned from the comparison are listed. We hope that our assessment will lead to new and more robust methods for cancer-specific transcript detection and help to produce more accurate functional impact predictions of isoform switching events.

A Comprehensive Map of mRNAs and Their Isoforms across All 14 Renal Tubule Segments of Mouse

BACKGROUND

The repertoire of protein expression along the renal tubule depends both on regulation of transcription and regulation of alternative splicing that can generate multiple proteins from a single gene.

METHODS

A full-length, small-sample RNA-seq protocol profiled transcriptomes for all 14 renal tubule segments microdissected from mouse kidneys.

RESULTS

This study identified >34,000 transcripts, including 3709 that were expressed in a segment-specific manner. All data are provided as an online resource (https://esbl.nhlbi.nih.gov/MRECA/Nephron/). Many of the genes expressed in unique patterns along the renal tubule were solute carriers, transcription factors, or G protein-coupled receptors that account for segment-specific function. Mapping the distribution of transcripts associated with Wnk-SPAK-PKA signaling, renin-angiotensin-aldosterone signaling, and cystic diseases of the kidney illustrated the applications of the online resource. The method allowed full-length mapping of RNA-seq reads, which facilitated comprehensive, unbiased characterization of alternative exon usage along the renal tubule, including known isoforms of Cldn10, Kcnj1 (ROMK), Slc12a1 (NKCC2), Wnk1, Stk39 (SPAK), and Slc14a2 (UT-A urea transporter). It also identified many novel isoforms with segment-specific distribution. These included variants associated with altered protein structure (Slc9a8, Khk, Tsc22d1, and Scoc), and variants that may affect untranslated, regulatory regions of transcripts (Pth1r, Pkar1a, and Dab2).

CONCLUSIONS

Full-length, unbiased sequencing of transcripts identified gene-expression patterns along the mouse renal tubule. The data, provided as an online resource, include both quantitative and qualitative differences in transcripts. Identification of alternative splicing along the renal tubule may prove critical to understanding renal physiology and pathophysiology.