Refining transcriptional programs in kidney development by integration of deep RNA-sequencing and array-based spatial profiling

Postnatal persistence of nonhuman primate sex-dependent renal structural and molecular changes programmed by intrauterine growth restriction

BACKGROUND

Poor nutrition during fetal development programs postnatal kidney function. Understanding postnatal consequences in nonhuman primates (NHP) is important for translation to our understanding the impact on human kidney function and disease risk. We hypothesized that intrauterine growth restriction (IUGR) in NHP persists postnatally, with potential molecular mechanisms revealed by Western-type diet challenge.

METHODS

IUGR juvenile baboons were fed a 7-week Western diet, with kidney biopsies, blood, and urine collected before and after challenge. Transcriptomics and metabolomics were used to analyze biosamples.

RESULTS

Pre-challenge IUGR kidney transcriptome and urine metabolome differed from controls. Post-challenge, sex and diet-specific responses in urine metabolite and renal signaling pathways were observed. Dysregulated mTOR signaling persisted postnatally in female pre-challenge. Post-challenge IUGR male response showed uncoordinated signaling suggesting proximal tubule injury.

CONCLUSION

Fetal undernutrition impacts juvenile offspring kidneys at the molecular level suggesting early-onset blood pressure dysregulation.

Identification and characterization of alternative STK39 transcripts within human and mouse kidneys reveals species-specific regulation of blood pressure

STK39 encodes a serine threonine kinase, SPAK, which is part of a multi-kinase network that determines renal Na+ reabsorption and blood pressure (BP) through regulation of sodium-chloride co-transporters in the kidney. Variants within STK39 are associated with susceptibility to essential hypertension, and constitutively active SPAK mice are hypertensive and hyperkalemic, similar to familial hyperkalemic hyperkalemia in humans. SPAK null mice are hypotensive and mimic Gitelman syndrome, a rare monogenic salt wasting human disorder. Mice exhibit nephron segment-specific expression of full length SPAK and N-terminally truncated SPAK isoforms (SPAK2 and KS-SPAK) with impaired kinase function. SPAK2 and KS-SPAK function to inhibit phosphorylation of cation co-transporters by full length SPAK. However, the existence of orthologous SPAK2 or KS-SPAK within the human kidney, and the role of such SPAK isoforms in nephron segment-specific regulation of Na+ reabsorption, still have not been determined. In this study, we examined both human and mouse kidney transcriptomes to uncover novel transcriptional regulation of STK39. We established that humans also express STK39 transcript isoforms similar to those found in mice but differ in abundance and are transcribed from human-specific promoters. In summary, STK39 undergoes species-specific transcriptional regulation, resulting in differentially expressed alternative transcripts that have implications for the design and testing of novel SPAK-targeting antihypertensive medications.

In developing mouse kidneys, orientation of loop of Henle growth is adaptive and guided by long-range cues from medullary collecting ducts

The path taken by the loop of Henle, from renal cortex to medulla and back, is critical to the ability of the kidney to concentrate urine and recover water. Unlike most developing tubules, which navigate as blind‐ended cylinders, the loop of Henle extends as a sharply bent loop, the apex of which leads the double tubes behind it in a ‘V’ shape. Here, we show that, in normal kidney development, loops of Henle extend towards the centroid of the kidney with an accuracy that increases the longer they extend. Using cultured kidney rudiments, and manipulations that rotate or remove portions of the organ, we show that loop orientation depends on long‐range cues from the medulla rather than either the orientation of the parent nephron or local cues in the cortex. The loops appear to be attracted to the most mature branch point of the collecting duct system but, if this is removed, they will head towards the most mature collecting duct branch available to them. Our results demonstrate the adaptive nature of guidance of this unusual example of a growing epithelium, and set the stage for later work devoted to understanding the molecules and mechanisms that underlie it.

The Lhx1-Ldb1 complex interacts with Furry to regulate microRNA expression during pronephric kidney development

The molecular events driving specification of the kidney have been well characterized. However, how the initial kidney field size is established, patterned, and proportioned is not well characterized. Lhx1 is a transcription factor expressed in pronephric progenitors and is required for specification of the kidney, but few Lhx1 interacting proteins or downstream targets have been identified. By tandem-affinity purification, we isolated FRY like transcriptional coactivator (Fryl), one of two paralogous genes, fryl and furry (fry), have been described in vertebrates. Both proteins were found to interact with the Ldb1-Lhx1 complex, but our studies focused on Lhx1/Fry functional roles, as they are expressed in overlapping domains. We found that Xenopus embryos depleted of fry exhibit loss of pronephric mesoderm, phenocopying the Lhx1-depleted animals. In addition, we demonstrated a synergism between Fry and Lhx1, identified candidate microRNAs regulated by the pair, and confirmed these microRNA clusters influence specification of the kidney. Therefore, our data shows that a constitutively-active Ldb1-Lhx1 complex interacts with a broadly expressed microRNA repressor, Fry, to establish the kidney field.

A novel FADS2 isoform identified in human milk fat globule suppresses FADS2 mediated Δ6-desaturation of omega-3 fatty acids

INTRODUCTION

The only known non-pharmacological means to alter long chain polyunsaturated fatty acid (LCPUFA) abundance in mammalian tissue is by altering substrate fatty acid ratios. Alternative mRNA splicing is increasingly recognized as a modulator of protein structure and function. Here we report identification of a novel alternative transcript (AT) of fatty acid desaturase 2 (FADS2) that inhibits production of omega-3 but not omega-6 LCPUFA, discovered during study of ATs in human milk fat globules (MFG).

METHODS

Human breastmilk collected from a single donor was used to isolate MFG. An mRNA-sequencing library was constructed from the total RNA isolated from the MFG. The constructed library was sequenced using an Illumina HiSeq instrument operating in high output mode. Expression levels of evolutionary conserved FADSAT were measured using cDNA from MFG by semi-quantitative RT-PCR assay.

RESULTS

RNA sequencing revealed >15,000 transcripts, including moderate expression of the FADS2 classical transcript (CS). A novel FADS2 alternative transcript (FADS2AT2) with 386 amino acids was discovered. When FADS2AT2 was transiently transfected into MCF7 cells stably expressing FADS2, delta-6 desaturation (D6D) of alpha-linolenic acid 18:3n-3 → 18:4n-3 was suppressed as were downstream products 20:4n-3 and 20:5n-3. In contrast, no significant effect on D6D of linoleic acid 18:2n-6 → 18:3n-6 or downstream products was observed. FADS2, FADS2AT1 and 5 out of 8 known FADS3AT were expressed in MFG. FADS1, FADS3AT3, and FADS3AT5 are undetectable.

CONCLUSION

The novel, noncatalytic FADS2AT2 regulates FADS2CS-mediated Δ6-desaturation of omega-3 but not omega-6 PUFA biosynthesis. This spliced isoform mediated interaction is the first molecular mechanism by which desaturation of one PUFA family but not the other is modulated.

Engineering kidney cells: reprogramming and directed differentiation to renal tissues

Growing knowledge of how cell identity is determined at the molecular level has enabled the generation of diverse tissue types, including renal cells from pluripotent or somatic cells. Recently, several in vitro protocols involving either directed differentiation or transcription-factor-based reprogramming to kidney cells have been established. Embryonic stem cells or induced pluripotent stem cells can be guided towards a kidney fate by exposing them to combinations of growth factors or small molecules. Here, renal development is recapitulated in vitro resulting in kidney cells or organoids that show striking similarities to mammalian embryonic nephrons. In addition, culture conditions are also defined that allow the expansion of renal progenitor cells in vitro. Another route towards the generation of kidney cells is direct reprogramming. Key transcription factors are used to directly impose renal cell identity on somatic cells, thus circumventing the pluripotent stage. This complementary approach to stem-cell-based differentiation has been demonstrated to generate renal tubule cells and nephron progenitors. In-vitro-generated renal cells offer new opportunities for modelling inherited and acquired renal diseases on a patient-specific genetic background. These cells represent a potential source for developing novel models for kidney diseases, drug screening and nephrotoxicity testing and might represent the first steps towards kidney cell replacement therapies. In this review, we summarize current approaches for the generation of renal cells in vitro and discuss the advantages of each approach and their potential applications.

Ablation of the epithelial-specific splicing factor Esrp1 results in ureteric branching defects and reduced nephron number

Background:

Abnormalities in ureteric bud (UB) branching morphogenesis lead to congenital anomalies of the kidney and reduced nephron numbers associated with chronic kidney disease (CKD) and hypertension. Previous studies showed that the epithelial fibroblast growth factor receptor 2 (Fgfr2) IIIb splice variant supports ureteric morphogenesis in response to ligands from the metanephric mesenchyme during renal organogenesis. The epithelial‐specific splicing regulator Esrp1 is required for expression of Fgfr2‐IIIb and other epithelial‐specific splice variants. Our objective was to determine whether Esrp1 is required for normal kidney development.

Results:

Ablation of Esrp1 in mice, alone or together with its paralog Esrp2, was associated with reduced kidney size and increased incidence of renal aplasia. Three‐dimensional imaging showed that embryonic Esrp1 knockout (KO) kidneys had fewer ureteric tips and reduced nephron numbers. Analysis of alternative splicing in Esrp‐null ureteric epithelial cells by RNA‐Seq confirmed a splicing switch in Fgfr2 as well as numerous other transcripts.

Conclusions:

Our findings reveal that Esrp1‐regulated splicing in ureteric epithelial cells plays an important role in renal development. Defects in Esrp1 KO kidneys likely reflect reduced and/or absent ureteric branching, leading to decreased nephron induction secondary to incorrect Fgfr2 splicing and other splicing alterations. Developmental Dynamics 245:991–1000, 2016. © 2016 The Authors. Developmental Dynamics published by Wiley Periodicals, Inc. on behalf of American Association of Anatomists.

Abnormalities in Ureteric bud (UB) branching morphogenesis lead to congenital anomalies of the kidney and reduced nephron numbers associated with chronic kidney disease (CKD) and hypertension.

We investigated the consequences of ablating the epithelial‐specific splicing regulator Esrp1 on renal organogenesis and determined that Esrp1 KO mice have reduced kidney size, fewer ureteric tips, reduced nephron numbers, and increased incidence of renal aplasia.

Analysis of alternative splicing in Esrp null ureteric epithelial cells by RNA‐Seq identified numerous alterations in splicing.

These findings reveal that Esrp1 regulated splicing in ureteric epithelial cells plays an important role in the kidney and illustrate the importance of alternative splicing for normal renal organogenesis.

An integrated miRNA functional screening and target validation method for organ morphogenesis

The relative ease of identifying microRNAs and their increasing recognition as important regulators of organogenesis motivate the development of methods to efficiently assess microRNA function during organ morphogenesis. In this context, embryonic organ explants provide a reliable and reproducible system that recapitulates some of the important early morphogenetic processes during organ development. Here we present a method to target microRNA function in explanted mouse embryonic organs. Our method combines the use of peptide-based nanoparticles to transfect specific microRNA inhibitors or activators into embryonic organ explants, with a microRNA pulldown assay that allows direct identification of microRNA targets. This method provides effective assessment of microRNA function during organ morphogenesis, allows prioritization of multiple microRNAs in parallel for subsequent genetic approaches, and can be applied to a variety of embryonic organs.

Genetic, environmental, and epigenetic factors involved in CAKUT

Congenital anomalies of the kidney and urinary tract (CAKUT) refer to a spectrum of structural renal malformations and are the leading cause of end-stage renal disease in children. The genetic diagnosis of CAKUT has proven to be challenging due to genetic and phenotypic heterogeneity and incomplete genetic penetrance. Monogenic causes of CAKUT have been identified using different approaches, including single gene screening, and gene panel and whole exome sequencing. The majority of the identified mutations, however, lack substantial evidence to support a pathogenic role in CAKUT. Copy number variants or single nucleotide variants that are associated with CAKUT have also been identified. Numerous studies support the influence of epigenetic and environmental factors on kidney development and the natural history of CAKUT, suggesting that the pathogenesis of this syndrome is multifactorial. In this Review we describe the current knowledge regarding the genetic susceptibility underlying CAKUT and the approaches used to investigate the genetic basis of CAKUT. We outline the associated environmental risk factors and epigenetic influences on CAKUT and discuss the challenges and strategies used to fully address the involvement and interplay of these factors in the pathogenesis of the disease.

Improving our resolution of kidney morphogenesis across time and space

As with many mammalian organs, size and cellular complexity represent considerable challenges to the comprehensive analysis of kidney organogenesis. Traditional analyses in the mouse have revealed early patterning events and spatial cellular relationships. However, an understanding of later events is lacking. The generation of a comprehensive temporospatial atlas of gene expression during kidney development has facilitated advances in lineage definition, as well as selective compartment ablation. Advances in quantitative and dynamic imaging have allowed comprehensive analyses at the level of organ, component tissue and cell across kidney organogenesis. Such approaches will enhance our understanding of the links between kidney development and final postnatal organ function. The final frontier will be translating this understanding to outcomes for renal disease in humans.

Cell-cell interactions driving kidney morphogenesis

The mammalian kidney forms via cell-cell interactions between an epithelial outgrowth of the nephric duct and the surrounding nephrogenic mesenchyme. Initial morphogenetic events include ureteric bud branching to form the collecting duct (CD) tree and mesenchymal-to-epithelial transitions to form the nephrons, requiring reciprocal induction between adjacent mesenchyme and epithelial cells. Within the tips of the branching ureteric epithelium, cells respond to mesenchyme-derived trophic factors by proliferation, migration, and mitosis-associated cell dispersal. Self-inhibition signals from one tip to another play a role in branch patterning. The position, survival, and fate of the nephrogenic mesenchyme are regulated by ECM and secreted signals from adjacent tip and stroma. Signals from the ureteric tip promote mesenchyme self-renewal and trigger nephron formation. Subsequent fusion to the CDs, nephron segmentation and maturation, and formation of a patent glomerular basement membrane also require specialized cell-cell interactions. Differential cadherin, laminin, nectin, and integrin expression, as well as intracellular kinesin and actin-mediated regulation of cell shape and adhesion, underlies these cell-cell interactions. Indeed, the capacity for the kidney to form via self-organization has now been established both via the recapitulation of expected morphogenetic interactions after complete dissociation and reassociation of cellular components during development as well as the in vitro formation of 3D kidney organoids from human pluripotent stem cells. As we understand more about how the many cell-cell interactions required for kidney formation operate, this enables the prospect of bioengineering replacement structures based on these self-organizing properties.

Defining nephrotic syndrome from an integrative genomics perspective

Nephrotic syndrome (NS) is a clinical condition with a high degree of morbidity and mortality, caused by failure of the glomerular filtration barrier, resulting in massive proteinuria. Our current diagnostic, prognostic and therapeutic decisions in NS are largely based upon clinical or histological patterns such as "focal segmental glomerulosclerosis" or "steroid sensitive". Yet these descriptive classifications lack the precision to explain the physiologic origins and clinical heterogeneity observed in this syndrome. A more precise definition of NS is required to identify mechanisms of disease and capture various clinical trajectories. An integrative genomics approach to NS applies bioinformatics and computational methods to comprehensive experimental, molecular and clinical data for holistic disease definition. A unique aspect is analysis of data together to discover NS-associated molecules, pathways, and networks. Integrating multidimensional datasets from the outset highlights how molecular lesions impact the entire individual. Data sets integrated range from genetic variation to gene expression, to histologic changes, to progression of chronic kidney disease (CKD). This review will introduce the tenets of integrative genomics and suggest how it can increase our understanding of NS from molecular and pathophysiological perspectives. A diverse group of genome-scale experiments are presented that have sought to define molecular signatures of NS. Finally, the Nephrotic Syndrome Study Network (NEPTUNE) will be introduced as an international, prospective cohort study of patients with NS that utilizes an integrated systems genomics approach from the outset. A major NEPTUNE goal is to achieve comprehensive disease definition from a genomics perspective and identify shared molecular drivers of disease.

MicroRNAs: potential regulators of renal development genes that contribute to CAKUT

Congenital anomalies of the kidney and urinary tract (CAKUT) are the leading cause of childhood chronic kidney disease (CKD). While mutations in several renal development genes have been identified as causes for CAKUT, most cases have not yet been linked to known mutations. Furthermore, the genotype-phenotype correlation is variable, suggesting that there might be additional factors that have an impact on the severity of CAKUT. MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression at the post-transcriptional level, and are involved in many developmental processes. Although little is known about the function of specific miRNAs in kidney development, several have recently been shown to regulate the expression of, and/or are regulated by, crucial renal development genes present in other organ systems. In this review, we discuss how miRNA regulation of common developmental signaling pathways may be applicable to renal development. We focus on genes that are known to contribute to CAKUT in humans, for which miRNA interactions in other contexts have been identified, with miRNAs that are present in the kidney. We hypothesize that miRNA-mediated processes might play a role in kidney development through similar mechanisms, and speculate that genotypic variations in these small RNAs or their targets could be associated with CAKUT.

Imperfect centered miRNA binding sites are common and can mediate repression of target mRNAs

Background

MicroRNAs (miRNAs) bind to mRNAs and target them for translational inhibition or transcriptional degradation. It is thought that most miRNA-mRNA interactions involve the seed region at the 5′ end of the miRNA. The importance of seed sites is supported by experimental evidence, although there is growing interest in interactions mediated by the central region of the miRNA, termed centered sites. To investigate the prevalence of these interactions, we apply a biotin pull-down method to determine the direct targets of ten human miRNAs, including four isomiRs that share centered sites, but not seeds, with their canonical partner miRNAs.

Results

We confirm that miRNAs and their isomiRs can interact with hundreds of mRNAs, and that imperfect centered sites are common mediators of miRNA-mRNA interactions. We experimentally demonstrate that these sites can repress mRNA activity, typically through translational repression, and are enriched in regions of the transcriptome bound by AGO. Finally, we show that the identification of imperfect centered sites is unlikely to be an artifact of our protocol caused by the biotinylation of the miRNA. However, the fact that there was a slight bias against seed sites in our protocol may have inflated the apparent prevalence of centered site-mediated interactions.

Conclusions

Our results suggest that centered site-mediated interactions are much more frequent than previously thought. This may explain the evolutionary conservation of the central region of miRNAs, and has significant implications for decoding miRNA-regulated genetic networks, and for predicting the functional effect of variants that do not alter protein sequence.

Dicer function is required in the metanephric mesenchyme for early kidney development

MicroRNAs (miRNAs) are small, noncoding regulatory RNAs that act as posttranscriptional repressors by binding to the 3'-untranslated region (3'-UTR) of target genes. They require processing by Dicer, an RNase III enzyme, to become mature regulatory RNAs. Previous work from our laboratory revealed critical roles for miRNAs in nephron progenitors at midgestation (Ho J, Pandey P, Schatton T, Sims-Lucas S, Khalid M, Frank MH, Hartwig S, Kreidberg JA. J Am Soc Nephrol 22: 1053-1063, 2011). To interrogate roles for miRNAs in the early metanephric mesenchyme, which gives rise to nephron progenitors as well as the renal stroma during kidney development, we conditionally ablated Dicer function in this lineage. Despite normal ureteric bud outgrowth and condensation of the metanephric mesenchyme to form nephron progenitors, early loss of miRNAs in the metanephric mesenchyme resulted in severe renal dysgenesis. Nephron progenitors are initially correctly specified in the mutant kidneys, with normal expression of several transcription factors known to be critical in progenitors, including Six2, Pax2, Sall1, and Wt1. However, there is premature loss of the nephron progenitor marker Cited1, marked apoptosis, and increased expression of the proapoptotic protein Bim shortly after the initial inductive events in early kidney development. Subsequently, there is a failure in ureteric bud branching and nephron progenitor differentiation. Taken together, our data demonstrate a previously undetermined requirement for miRNAs during early kidney organogenesis and indicate a crucial role for miRNAs in regulating the survival of this lineage.

Divergent transcription is associated with promoters of transcriptional regulators

Background

Divergent transcription is a wide-spread phenomenon in mammals. For instance, short bidirectional transcripts are a hallmark of active promoters, while longer transcripts can be detected antisense from active genes in conditions where the RNA degradation machinery is inhibited. Moreover, many described long non-coding RNAs (lncRNAs) are transcribed antisense from coding gene promoters. However, the general significance of divergent lncRNA/mRNA gene pair transcription is still poorly understood. Here, we used strand-specific RNA-seq with high sequencing depth to thoroughly identify antisense transcripts from coding gene promoters in primary mouse tissues.

Results

We found that a substantial fraction of coding-gene promoters sustain divergent transcription of long non-coding RNA (lncRNA)/mRNA gene pairs. Strikingly, upstream antisense transcription is significantly associated with genes related to transcriptional regulation and development. Their promoters share several characteristics with those of transcriptional developmental genes, including very large CpG islands, high degree of conservation and epigenetic regulation in ES cells. In-depth analysis revealed a unique GC skew profile at these promoter regions, while the associated coding genes were found to have large first exons, two genomic features that might enforce bidirectional transcription. Finally, genes associated with antisense transcription harbor specific H3K79me2 epigenetic marking and RNA polymerase II enrichment profiles linked to an intensified rate of early transcriptional elongation.

Conclusions

We concluded that promoters of a class of transcription regulators are characterized by a specialized transcriptional control mechanism, which is directly coupled to relaxed bidirectional transcription.

MicroRNAs in renal development

The discovery of microRNAs (miRNAs) as novel regulators of gene expression has led to a marked change in how gene regulation is viewed, with important implications for development and disease. MiRNAs are endogenous, small, noncoding RNAs that largely repress their target mRNAs post-transcriptionally. The regulation of gene expression by miRNAs represents an evolutionarily conserved mechanism that is broadly applicable to most biological processes. Recent studies have begun to define the role of miRNAs in different cell lineages during kidney development, and to implicate specific miRNAs in developmental and pathophysiological processes in the kidney. This review will focus on novel insights into the role(s) of miRNAs in kidney development, and discuss the implications for pediatric renal disease.

MicroRNAs in kidney disease: an emerging understanding

MicroRNAs (miRNAs) are short noncoding RNA molecules that function by negatively regulating the expression of their target genes in a tightly controlled manner. Accumulating evidence, based in part on effects seen after miRNA overexpression and/or knockdown, points to the critical involvement of miRNAs in kidney function in health and disease. In this review, we provide a quick overview of the biogenesis of miRNAs and their potential involvement in kidney development and normal function. We also discuss the current literature that has begun to uncover the role of miRNAs in the pathogenesis of kidney diseases, including diabetic nephropathy, hypertension, glomerulonephritis, and cancer. As such, miRNAs have potential utility in the clinical realm as disease biomarkers. Moreover, miRNAs represent an attractive therapeutic target for a number of kidney diseases. We close by discussing a number of potential challenges that face the field of miRNA research and clinical use.

Application of next generation qPCR and sequencing platforms to mRNA biomarker analysis

Recent years have seen the emergence of new high-throughput PCR and sequencing platforms with the potential to bring analysis of transcriptional biomarkers to a broader range of clinical applications and to provide increasing depth to our understanding of the transcriptome. We present an overview of how to process clinical samples for RNA biomarker analysis in terms of RNA extraction and mRNA enrichment, and guidelines for sample analysis by RT-qPCR and digital PCR using nanofluidic real-time PCR platforms. The options for quantitative gene expression profiling and whole transcriptome sequencing by next generation sequencing are reviewed alongside the bioinformatic considerations for these approaches. Considering the diverse technologies now available for transcriptome analysis, methods for standardising measurements between platforms will be paramount if their diagnostic impact is to be maximised. Therefore, the use of RNA standards and other reference materials is also discussed.

RNA-Seq defines novel genes, RNA processing patterns and enhancer maps for the early stages of nephrogenesis: Hox supergenes

During kidney development the cap mesenchyme progenitor cells both self renew and differentiate into nephrons. The balance between renewal and differentiation determines the final nephron count, which is of considerable medical importance. An important goal is to create a precise genetic definition of the early differentiation of cap mesenchyme progenitors. We used RNA-Seq to transcriptional profile the cap mesenchyme progenitors and their first epithelial derivative, the renal vesicles. The results provide a global view of the changing gene expression program during this key period, defining expression levels for all transcription factors, growth factors, and receptors. The RNA-Seq was performed using two different biochemistries, with one examining only polyadenylated RNA and the other total RNA. This allowed the analysis of noncanonical transcripts, which for many genes were more abundant than standard exonic RNAs. Since a large fraction of enhancers are now known to be transcribed the results also provide global maps of potential enhancers. Further, the RNA-Seq data defined hundreds of novel splice patterns and large numbers of new genes. Particularly striking was the extensive sense/antisense transcription and changing RNA processing complexities of the Hox clusters.