Characterization of Alternative Splicing in High-Risk Wilms' Tumors.

The significant heterogeneity of Wilms' tumors between different patients is thought to arise from genetic and epigenetic distortions that occur during various stages of fetal kidney development in a way that is poorly understood. To address this, we characterized the heterogeneity of alternative mRNA splicing in Wilms' tumors using a publicly available RNAseq dataset of high-risk Wilms' tumors and normal kidney samples. Through Pareto task inference and cell deconvolution, we found that the tumors and normal kidney samples are organized according to progressive stages of kidney development within a triangle-shaped region in latent space, whose vertices, or "archetypes", resemble the cap mesenchyme, the nephrogenic stroma, and epithelial tubular structures of the fetal kidney. We identified a set of genes that are alternatively spliced between tumors located in different regions of latent space and found that many of these genes are associated with the epithelial-to-mesenchymal transition (EMT) and muscle development. Using motif enrichment analysis, we identified putative splicing regulators, some of which are associated with kidney development. Our findings provide new insights into the etiology of Wilms' tumors and suggest that specific splicing mechanisms in early stages of development may contribute to tumor development in different patients.

OCT4 induces long-lived dedifferentiated kidney progenitors poised to redifferentiate in 3D kidney spheroids.

Upscaling of kidney epithelial cells is crucial for renal regenerative medicine. Nonetheless, the adult kidney lacks a distinct stem cell hierarchy, limiting the ability to long-term propagate clonal populations of primary cells that retain renal identity. Toward this goal, we tested the paradigm of shifting the balance between differentiation and stemness in the kidney by introducing a single pluripotency factor, OCT4. Here we show that ectopic expression of OCT4 in human adult kidney epithelial cells (hKEpC) induces the cells to dedifferentiate, stably proliferate, and clonally emerge over many generations. Control hKEpC dedifferentiate, assume fibroblastic morphology, and completely lose clonogenic capacity. Analysis of gene expression and histone methylation patterns revealed that OCT4 represses the HNF1B gene module, which is critical for kidney epithelial differentiation, and concomitantly activates stemness-related pathways. OCT4-hKEpC can be long-term expanded in the dedifferentiated state that is primed for renal differentiation. Thus, when expanded OCT4-hKEpC are grown as kidney spheroids (OCT4-kSPH), they reactivate the HNF1B gene signature, redifferentiate, and efficiently generate renal structures in vivo. Hence, changes occurring in the cellular state of hKEpC following OCT4 induction, long-term propagation, and 3D aggregation afford rapid scale-up technology of primary renal tissue-forming cells.

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.

Applications of Community Detection Algorithms to Large Biological Datasets.

Recent advances in data acquiring technologies in biology have led to major challenges in mining relevant information from large datasets. For example, single-cell RNA sequencing technologies are producing expression and sequence information from tens of thousands of cells in every single experiment. A common task in analyzing biological data is to cluster samples or features (e.g., genes) into groups sharing common characteristics. This is an NP-hard problem for which numerous heuristic algorithms have been developed. However, in many cases, the clusters created by these algorithms do not reflect biological reality. To overcome this, a Networks Based Clustering (NBC) approach was recently proposed, by which the samples or genes in the dataset are first mapped to a network and then community detection (CD) algorithms are used to identify clusters of nodes.Here, we created an open and flexible python-based toolkit for NBC that enables easy and accessible network construction and community detection. We then tested the applicability of NBC for identifying clusters of cells or genes from previously published large-scale single-cell and bulk RNA-seq datasets.We show that NBC can be used to accurately and efficiently analyze large-scale datasets of RNA sequencing experiments.

Human kidney clonal proliferation disclose lineage-restricted precursor characteristics.

In-vivo single cell clonal analysis in the adult mouse kidney has previously shown lineage-restricted clonal proliferation within varying nephron segments as a mechanism responsible for cell replacement and local regeneration. To analyze ex-vivo clonal growth, we now preformed limiting dilution to generate genuine clonal cultures from one single human renal epithelial cell, which can give rise to up to 3.4 * 106 cells, and analyzed their characteristics using transcriptomics. A comparison between clonal cultures revealed restriction to either proximal or distal kidney sub-lineages with distinct cellular and molecular characteristics; rapidly amplifying de-differentiated clones and a stably proliferating cuboidal epithelial-appearing clones, respectively. Furthermore, each showed distinct molecular features including cell-cycle, epithelial-mesenchymal transition, oxidative phosphorylation, BMP signaling pathway and cell surface markers. In addition, analysis of clonal versus bulk cultures show early clones to be more quiescent, with elevated expression of renal developmental genes and overall reduction in renal identity markers, but with an overlapping expression of nephron segment identifiers and multiple identity. Thus, ex-vivo clonal growth mimics the in-vivo situation displaying lineage-restricted precursor characteristics of mature renal cells. These data suggest that for reconstruction of varying renal lineages with human adult kidney based organoid technology and kidney regeneration ex-vivo, use of multiple heterogeneous precursors is warranted.

LEFTY1 Is a Dual-SMAD Inhibitor that Promotes Mammary Progenitor Growth and Tumorigenesis.

SMAD pathways govern epithelial proliferation, and transforming growth factor ß (TGF-ß and BMP signaling through SMAD members has distinct effects on mammary development and homeostasis. Here, we show that LEFTY1, a secreted inhibitor of NODAL/SMAD2 signaling, is produced by mammary progenitor cells and, concomitantly, suppresses SMAD2 and SMAD5 signaling to promote long-term proliferation of normal and malignant mammary epithelial cells. In contrast, BMP7, a NODAL antagonist with context-dependent functions, is produced by basal cells and restrains progenitor cell proliferation. In normal mouse epithelium, LEFTY1 expression in a subset of luminal cells and rare basal cells opposes BMP7 to promote ductal branching. LEFTY1 binds BMPR2 to suppress BMP7-induced activation of SMAD5, and this LEFTY1-BMPR2 interaction is specific to tumor-initiating cells in triple-negative breast cancer xenografts that rely on LEFTY1 for growth. These results suggest that LEFTY1 is an endogenous dual-SMAD inhibitor and that suppressing its function may represent a therapeutic vulnerability in breast cancer.

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

BACKGROUND: During mammalian kidney development, nephron progenitors undergo a mesenchymal-to-epithelial transition and eventually differentiate into the various tubular segments of the nephron. Recently, Drop-seq single-cell RNA sequencing technology for measuring gene expression from thousands of individual cells identified the different cell types in the developing kidney. However, that analysis did not include the additional layer of heterogeneity that alternative mRNA splicing creates. METHODS: Full transcript length single-cell RNA sequencing characterized the transcriptomes of 544 individual cells from mouse embryonic kidneys. RESULTS: Gene expression levels measured with full transcript length single-cell RNA sequencing identified each cell type. Further analysis comprehensively characterized splice isoform switching during the transition between mesenchymal and epithelial cellular states, which is a key transitional process in kidney development. The study also identified several putative splicing regulators, including the genes Esrp1/2 and Rbfox1/2. CONCLUSIONS: Discovery of the sets of genes that are alternatively spliced as the fetal kidney mesenchyme differentiates into tubular epithelium will improve our understanding of the molecular mechanisms that drive kidney development.

Ex Vivo Expanded 3D Human Kidney Spheres Engraft Long Term and Repair Chronic Renal Injury in Mice.

End-stage renal disease is a worldwide epidemic requiring renal replacement therapy. Harvesting tissue from failing kidneys and autotransplantation of tissue progenitors could theoretically delay the need for dialysis. Here we use healthy and end-stage human adult kidneys to robustly expand proliferative kidney epithelial cells and establish 3D kidney epithelial cultures termed "nephrospheres." Formation of nephrospheres reestablishes renal identity and function in primary cultures. Transplantation into NOD/SCID mice shows that nephrospheres restore self-organogenetic properties lost in monolayer cultures, allowing long-term engraftment as tubular structures, potentially adding nephron segments and demonstrating self-organization as critical to survival. Furthermore, long-term tubular engraftment of nephrospheres is functionally beneficial in murine models of chronic kidney disease. Remarkably, nephrospheres inhibit pro-fibrotic collagen production in cultured fibroblasts via paracrine modulation, while transplanted nephrospheres induce transcriptional signatures of proliferation and release from quiescence, suggesting re-activation of endogenous repair. These data support the use of human nephrospheres for renal cell therapy.

Chromatin occupancy and epigenetic analysis reveal new insights into the function of the GATA1 N terminus in erythropoiesis.

Mutations in GATA1, which lead to expression of the GATA1s isoform that lacks the GATA1 N terminus, are seen in patients with Diamond-Blackfan anemia (DBA). In our efforts to better understand the connection between GATA1s and DBA, we comprehensively studied erythropoiesis in Gata1s mice. Defects in yolks sac and fetal liver hematopoiesis included impaired terminal maturation and reduced numbers of erythroid progenitors. RNA-sequencing revealed that both erythroid and megakaryocytic gene expression patterns were altered by the loss of the N terminus, including aberrant upregulation of Gata2 and Runx1. Dysregulation of global H3K27 methylation was found in the erythroid progenitors upon loss of N terminus of GATA1. Chromatin-binding assays revealed that, despite similar occupancy of GATA1 and GATA1s, there was a striking reduction of H3K27me3 at regulatory elements of the Gata2 and Runx1 genes. Consistent with the observation that overexpression of GATA2 has been reported to impair erythropoiesis, we found that haploinsufficiency of Gata2 rescued the erythroid defects of Gata1s fetuses. Together, our integrated genomic analysis of transcriptomic and epigenetic signatures reveals that, Gata1 mice provide novel insights into the role of the N terminus of GATA1 in transcriptional regulation and red blood cell maturation which may potentially be useful for DBA patients.

Macrophages restrict the nephrogenic field and promote endothelial connections during kidney development.

The origins and functions of kidney macrophages in the adult have been explored, but their roles during development remain largely unknown. Here we characterise macrophage arrival, localisation, heterogeneity, and functions during kidney organogenesis. Using genetic approaches to ablate macrophages, we identify a role for macrophages in nephron progenitor cell clearance as mouse kidney development begins. Throughout renal organogenesis, most kidney macrophages are perivascular and express F4/80 and CD206. These macrophages are enriched for mRNAs linked to developmental processes, such as blood vessel morphogenesis. Using antibody-mediated macrophage-depletion, we show macrophages support vascular anastomoses in cultured kidney explants. We also characterise a subpopulation of galectin-3+ (Gal3+) myeloid cells within the developing kidney. Our findings may stimulate research into macrophage-based therapies for renal developmental abnormalities and have implications for the generation of bioengineered kidney tissues.

A cluster robustness score for identifying cell subpopulations in single cell gene expression datasets from heterogeneous tissues and tumors.

MOTIVATION: A major aim of single cell biology is to identify important cell types such as stem cells in heterogeneous tissues and tumors. This is typically done by isolating hundreds of individual cells and measuring expression levels of multiple genes simultaneously from each cell. Then, clustering algorithms are used to group together similar single-cell expression profiles into clusters, each representing a distinct cell type. However, many of these clusters result from overfitting, meaning that rather than representing biologically meaningful cell types, they describe the intrinsic 'noise' in gene expression levels due to limitations in experimental precision or the intrinsic randomness of biochemical cellular processes. Consequentially, these non-meaningful clusters are most sensitive to noise: a slight shift in gene expression levels due to a repeated measurement will rearrange the grouping of data points such that these clusters break up. RESULTS: To identify the biologically meaningful clusters we propose a 'cluster robustness score': We add increasing amounts of noise (zero mean and increasing variance) and check which clusters are most robust in the sense that they do not mix with their neighbors up to high levels of noise. We show that biologically meaningful cell clusters that were manually identified in previously published single cell expression datasets have high robustness scores. These scores are higher than what would be expected in corresponding randomized homogeneous datasets having the same expression level statistics. We believe that this scoring system provides a more automated way to identify cell types in heterogeneous tissues and tumors. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

In Vivo Expansion of Cancer Stemness Affords Novel Cancer Stem Cell Targets: Malignant Rhabdoid Tumor as an Example.

Cancer stem cell (CSC) identification relies on transplantation assays of cell subpopulations sorted from fresh tumor samples. Here, we attempt to bypass limitations of abundant tumor source and predetermined immune selection by in vivo propagating patient-derived xenografts (PDX) from human malignant rhabdoid tumor (MRT), a rare and lethal pediatric neoplasm, to an advanced state in which most cells behave as CSCs. Stemness is then probed by comparative transcriptomics of serial PDXs generating a gene signature of epithelial to mesenchymal transition, invasion/motility, metastasis, and self-renewal, pinpointing putative MRT CSC markers. The relevance of these putative CSC molecules is analyzed by sorting tumorigenic fractions from early-passaged PDX according to one such molecule, deciphering expression in archived primary tumors, and testing the effects of CSC molecule inhibition on MRT growth. Using this platform, we identify ALDH1 and lysyl oxidase (LOX) as relevant targets and provide a larger framework for target and drug discovery in rare pediatric cancers.

Geometry of Gene Expression Space of Wilms' Tumors From Human Patients.

Wilms' tumor is a pediatric malignancy that is thought to originate from faulty kidney development during the embryonic stage. However, there is a large variation between tumors from different patients in both histology and gene expression that is not well characterized. Here we use a meta-analysis of published microarray datasets to show that Favorable Histology Wilms' Tumors (FHWT's) fill a triangle-shaped continuum in gene expression space of which the vertices represent three idealized "archetypes". We show that these archetypes have predominantly renal blastemal, stromal, and epithelial characteristics and that they correlate well with the three major lineages of the developing embryonic kidney. Moreover, we show that advanced stage tumors shift towards the renal blastemal archetype. These results illustrate the potential of this methodology for characterizing the cellular composition of Wilms' tumors and for assessing disease progression.

A brief review of single-cell transcriptomic technologies.

In recent years, there has been an effort to develop new technologies for measuring gene expression and sequence information from thousands of individual cells. Large data sets that were obtained using these 'single cell' technologies have allowed scientists to address fundamental questions in biomedicine ranging from stems cells and development to cancer and immunology. Here, we provide a brief review of recent developments in single-cell technology. Our intention is to provide a quick background for newcomers to the field as well as a deeper description of some of the leading technologies to date.

Role of epithelial to mesenchymal transition associated genes in mammary gland regeneration and breast tumorigenesis.

Previous studies have proposed that epithelial to mesenchymal transition (EMT) in breast cancer cells regulates metastasis, stem cell properties and chemo-resistance; most studies were based on in vitro culture of cell lines and mouse transgenic cancer models. However, the identity and function of cells expressing EMT-associated genes in normal murine mammary gland homeostasis and human breast cancer still remains under debate. Using in vivo lineage tracing and triple negative breast cancer (TNBC) patient derived xenografts we demonstrate that the repopulating capacity in normal mammary epithelial cells and tumorigenic capacity in TNBC is independent of expression of EMT-associated genes. In breast cancer, while a subset of cells with epithelial and mesenchymal phenotypes have stem cell activity, in many cells that have lost epithelial characteristics with increased expression of mesenchymal genes, have decreased tumor-initiating capacity and plasticity. These findings have implications for the development of effective therapeutic agents targeting tumor-initiating cells.

Evidence of In Vitro Preservation of Human Nephrogenesis at the Single-Cell Level.

During nephrogenesis, stem/progenitor cells differentiate and give rise to early nephron structures that segment to proximal and distal nephron cell types. Previously, we prospectively isolated progenitors from human fetal kidney (hFK) utilizing a combination of surface markers. However, upon culture nephron progenitors differentiated and could not be robustly maintained in vitro. Here, by culturing hFK in a modified medium used for in vitro growth of mouse nephron progenitors, and by dissection of NCAM+/CD133- progenitor cells according to EpCAM expression (NCAM+/CD133-/EpCAM-, NCAM+/CD133-/EpCAMdim, NCAM+/CD133-/EpCAMbright), we show at single-cell resolution a preservation of uninduced and induced cap mesenchyme as well as a transitioning mesenchymal-epithelial state. Concomitantly, differentiating and differentiated epithelial lineages are also maintained. In vitro expansion of discrete stages of early human nephrogenesis in nephron stem cell cultures may be used for drug screening on a full repertoire of developing kidney cells and for prospective isolation of mesenchymal or epithelial renal lineages for regenerative medicine.

Peroxisome proliferator-activated receptor gamma (PPARγ) is central to the initiation and propagation of human angiomyolipoma, suggesting its potential as a therapeutic target

Angiomyolipoma (AML), the most common benign renal tumor, can result in severe morbidity from hemorrhage and renal failure. While mTORC1 activation is involved in its growth, mTORC1 inhibitors fail to eradicate AML, highlighting the need for new therapies. Moreover, the identity of the AML cell of origin is obscure. AML research, however, is hampered by the lack of in vivo models. Here, we establish a human AML-xenograft (Xn) model in mice, recapitulating AML at the histological and molecular levels. Microarray analysis demonstrated tumor growth in vivo to involve robust PPARγ-pathway activation. Similarly, immunostaining revealed strong PPARγ expression in human AML specimens. Accordingly, we demonstrate that while PPARγ agonism accelerates AML growth, PPARγ antagonism is inhibitory, strongly suppressing AML proliferation and tumor-initiating capacity, via a TGFB-mediated inhibition of PDGFB and CTGF. Finally, we show striking similarity between AML cell lines and mesenchymal stem cells (MSCs) in terms of antigen and gene expression and differentiation potential. Altogether, we establish the first in vivo human AML model, which provides evidence that AML may originate in a PPARγ-activated renal MSC lineage that is skewed toward adipocytes and smooth muscle and away from osteoblasts, and uncover PPARγ as a regulator of AML growth, which could serve as an attractive therapeutic target.