Analysis of metagene portraits reveals distinct transitions during kidney organogenesis

Early prediction of developing spontaneous activity in cultured neuronal networks

Synchronization and bursting activity are intrinsic electrophysiological properties of in vivo and in vitro neural networks. During early development, cortical cultures exhibit a wide repertoire of synchronous bursting dynamics whose characterization may help to understand the parameters governing the transition from immature to mature networks. Here we used machine learning techniques to characterize and predict the developing spontaneous activity in mouse cortical neurons on microelectrode arrays (MEAs) during the first three weeks in vitro. Network activity at three stages of early development was defined by 18 electrophysiological features of spikes, bursts, synchrony, and connectivity. The variability of neuronal network activity during early development was investigated by applying k-means and self-organizing map (SOM) clustering analysis to features of bursts and synchrony. These electrophysiological features were predicted at the third week in vitro with high accuracy from those at earlier times using three machine learning models: Multivariate Adaptive Regression Splines, Support Vector Machines, and Random Forest. Our results indicate that initial patterns of electrical activity during the first week in vitro may already predetermine the final development of the neuronal network activity. The methodological approach used here may be applied to explore the biological mechanisms underlying the complex dynamics of spontaneous activity in developing neuronal cultures.

Regulation of Solute Carriers OCT2 and OAT1/3 in the Kidney: A Phylogenetic, Ontogenetic and Cell Dynamic Perspective

Over the course of more than 500 million years, the kidneys have undergone a remarkable evolution from primitive nephric tubes to intricate filtration-reabsorption systems that maintain homeostasis and remove metabolic end products from the body. The evolutionarily conserved solute carriers Organic Cation Transporter 2 (OCT2), and Organic Anion Transporters 1 and 3 (OAT1/3) coordinate the active secretion of a broad range of endogenous and exogenous substances, many of which accumulate in the blood of patients with kidney failure despite dialysis. Harnessing OCT2 and OAT1/3 through functional preservation or regeneration could alleviate the progression of kidney disease. Additionally, it would improve current in vitro test models that lose their expression in culture. With this review, we explore OCT2 and OAT1/3 regulation using different perspectives: phylogenetic, ontogenetic and cell dynamic. Our aim is to identify possible molecular targets to both help prevent or compensate for the loss of transport activity in patients with kidney disease, and to enable endogenous OCT2 and OAT1/3 induction in vitro in order to develop better models for drug development.

Beyond Boltzmann-Gibbs-Shannon in Physics and Elsewhere

The pillars of contemporary theoretical physics are classical mechanics, Maxwell electromagnetism, relativity, quantum mechanics, and Boltzmann–Gibbs (BG) statistical mechanics –including its connection with thermodynamics. The BG theory describes amazingly well the thermal equilibrium of a plethora of so-called simple systems. However, BG statistical mechanics and its basic additive entropy SBG started, in recent decades, to exhibit failures or inadequacies in an increasing number of complex systems. The emergence of such intriguing features became apparent in quantum systems as well, such as black holes and other area-law-like scenarios for the von Neumann entropy. In a different arena, the efficiency of the Shannon entropy—as the BG functional is currently called in engineering and communication theory—started to be perceived as not necessarily optimal in the processing of images (e.g., medical ones) and time series (e.g., economic ones). Such is the case in the presence of generic long-range space correlations, long memory, sub-exponential sensitivity to the initial conditions (hence vanishing largest Lyapunov exponents), and similar features. Finally, we witnessed, during the last two decades, an explosion of asymptotically scale-free complex networks. This wide range of important systems eventually gave support, since 1988, to the generalization of the BG theory. Nonadditive entropies generalizing the BG one and their consequences have been introduced and intensively studied worldwide. The present review focuses on these concepts and their predictions, verifications, and applications in physics and elsewhere. Some selected examples (in quantum information, high- and low-energy physics, low-dimensional nonlinear dynamical systems, earthquakes, turbulence, long-range interacting systems, and scale-free networks) illustrate successful applications. The grounding thermodynamical framework is briefly described as well.

Evolution, kidney development, and chronic kidney disease

There is a global epidemic of chronic kidney disease (CKD) characterized by a progressive loss of nephrons, ascribed in large part to a rising incidence of hypertension, metabolic syndrome, and type 2 diabetes mellitus. There is a ten-fold variation in nephron number at birth in the general population, and a 50% overall decrease in nephron number in the last decades of life. The vicious cycle of nephron loss stimulating hypertrophy by remaining nephrons and resulting in glomerulosclerosis has been regarded as maladaptive, and only partially responsive to angiotensin inhibition. Advances over the past century in kidney physiology, genetics, and development have elucidated many aspects of nephron formation, structure and function. Parallel advances have been achieved in evolutionary biology, with the emergence of evolutionary medicine, a discipline that promises to provide new insight into the treatment of chronic disease. This review provides a framework for understanding the origins of contemporary developmental nephrology, and recent progress in evolutionary biology. The establishment of evolutionary developmental biology (evo-devo), ecological developmental biology (eco-devo), and developmental origins of health and disease (DOHaD) followed the discovery of the hox gene family, the recognition of the contribution of cumulative environmental stressors to the changing phenotype over the life cycle, and mechanisms of epigenetic regulation. The maturation of evolutionary medicine has contributed to new investigative approaches to cardiovascular disease, cancer, and infectious disease, and promises the same for CKD. By incorporating these principles, developmental nephrology is ideally positioned to answer important questions regarding the fate of nephrons from embryo through senescence.

Developmental regulation of kidney and liver solute carrier and ATP-binding cassette drug transporters and drug metabolizing enzymes: the role of remote organ communication

INTRODUCTION

The ontogeny of drug transport and metabolism is generally studied independently in tissues, yet in the immediate postnatal period the developmental regulation of SLC and ABC transporters and metabolizing enzymes must be coordinated. Using the Remote Sensing and Signaling Hypothesis as a framework, we describe how a systems physiology view helps to make sense of how inter-organ communication via hepatic, renal, and intestinal transporters and drug metabolizing enzymes (DMEs) is regulated from the immediate postnatal period through adulthood. Areas covered: This review examines patterns of developmental expression and function of transporters and DMEs with a focus on how cross-talk between these proteins in the kidney, liver and other organs (e.g., intestine) may be coordinated postnatally to optimize levels of metabolites and endogenous signaling molecules as well as gut-microbiome products. Expert opinion/commentary: Developmental expression is considered in terms of the Remote Sensing and Signaling Hypothesis, which addresses how transporters and DMEs participate in inter-organ and inter-organism small molecule communication in health, development, and disease. This hypothesis, for which there is growing support, is particularly relevant to the 'birth transition' and post-natal developmental physiology when organs must deal with critical physiological tasks distinct from the fetal period and where remote inter-organ and possibly inter-organismal (e.g. infant-gut microbiome) communication is likely to be critical to maintain homeostasis.

Developmental Programming of Branching Morphogenesis in the Kidney

The kidney developmental program encodes the intricate branching and organization of approximately 1 million functional units (nephrons). Branching regulation is poorly understood, as is the source of a 10-fold variation in nephron number. Notably, low nephron count increases the risk for developing hypertension and renal failure. To better understand the source of this variation, we analyzed the complete gestational trajectory of mouse kidney development. We constructed a computerized architectural map of the branching process throughout fetal life and found that organogenesis is composed of two distinct developmental phases, each with stage-specific rate and morphologic parameters. The early phase is characterized by a rapid acceleration in branching rate and by branching divisions that repeat with relatively reproducible morphology. The latter phase, however, is notable for a significantly decreased yet constant branching rate and the presence of nonstereotyped branching events that generate progressive variability in tree morphology until birth. Our map identifies and quantitates the contribution of four developmental mechanisms that guide organogenesis: growth, patterning, branching rate, and nephron induction. When applied to organs that developed under conditions of malnutrition or in the setting of growth factor mutation, our normative map provided an essential link between kidney architecture and the fundamental morphogenetic mechanisms that guide development. This morphogenetic map is expected to find widespread applications and help identify modifiable targets to prevent developmental programming of common diseases.

The organic anion transporter (OAT) family: a systems biology perspective

The organic anion transporter (OAT) subfamily, which constitutes roughly half of the SLC22 (solute carrier 22) transporter family, has received a great deal of attention because of its role in handling of common drugs (antibiotics, antivirals, diuretics, nonsteroidal anti-inflammatory drugs), toxins (mercury, aristolochic acid), and nutrients (vitamins, flavonoids). Oats are expressed in many tissues, including kidney, liver, choroid plexus, olfactory mucosa, brain, retina, and placenta. Recent metabolomics and microarray data from Oat1 [Slc22a6, originally identified as NKT (novel kidney transporter)] and Oat3 (Slc22a8) knockouts, as well as systems biology studies, indicate that this pathway plays a central role in the metabolism and handling of gut microbiome metabolites as well as putative uremic toxins of kidney disease. Nuclear receptors and other transcription factors, such as Hnf4α and Hnf1α, appear to regulate the expression of certain Oats in conjunction with phase I and phase II drug metabolizing enzymes. Some Oats have a strong selectivity for particular signaling molecules, including cyclic nucleotides, conjugated sex steroids, odorants, uric acid, and prostaglandins and/or their metabolites. According to the "Remote Sensing and Signaling Hypothesis," which is elaborated in detail here, Oats may function in remote interorgan communication by regulating levels of signaling molecules and key metabolites in tissues and body fluids. Oats may also play a major role in interorganismal communication (via movement of small molecules across the intestine, placental barrier, into breast milk, and volatile odorants into the urine). The role of various Oat isoforms in systems physiology appears quite complex, and their ramifications are discussed in the context of remote sensing and signaling.

GDNF and MAPK-ERK pathway signaling is reduced during nephrogenesis following maternal under-nutrition

Maternal under-nutrition (MUN) during gestation results in growth-restricted newborns with reduced glomerular number and subsequent hypertension. We investigated dysregulation of glial derived neurotrophic factor (GDNF) and MAPK-ERK (mitogen-activated protein kinase-extracellular signal-regulated protein kinase) signal pathway gene expression following MUN. MUN rats were 50% food restricted from embryonic day 10 till postnatal day 1. Kidneys were harvested at embryonic day (E)20, and postnatal days (P)1 and 21. Kidney protein expression was determined by Western blot. At E20, protein expression of growth factor receptor alpha 1 (GFRα1) and phosphorylated ERK1/2 and mitogen-activated protein kinase kinase (MEK)1/2 were reduced significantly, and immunohistochemistry confirmed reduction of phosphorylated ERK (pERK) with maintenance of pERK localization. Total MEK and ERK were unchanged. At P1, only GFRα1 and pERK1/2 were reduced significantly while at P21, expression of all growth factors except total MEK was unchanged. Total MEK was increased. Glomerular number was decreased by 19% in P21 kidneys and blood pressure was increased in 12-week-old rats. In conclusion, GDNF and MAPK-ERK signaling are dysregulated during active nephrogenesis in fetal and early newborn offspring kidneys in the MUN model. This may be a key mechanism in reduced offspring nephrogenesis and programmed hypertension.

Hepatocyte nuclear factors 4α and 1α regulate kidney developmental expression of drug-metabolizing enzymes and drug transporters

The transcriptional regulation of drug-metabolizing enzymes and transporters (here collectively referred to as DMEs) in the developing proximal tubule (PT) is not well understood. As in the liver, DME regulation in the PT may be mediated through nuclear receptors, which are thought to "sense" deviations from homeostasis by being activated by ligands, some of which are handled by DMEs, including drug transporters. Systems analysis of transcriptomic data during kidney development predicted a set of upstream transcription factors, including hepatocyte nuclear factor 4α (Hnf4a) and Hnf1a, as well as Nr3c1 (Gr), Nfe2l2 (Nrf2), peroxisome proliferator-activated receptor α (Pparα), and Tp53. Motif analysis of cis-regulatory enhancers further suggested that Hnf4a and Hnf1a are the main transcriptional regulators of DMEs in the PT. Available expression data from tissue-specific Hnf4a knockout tissues revealed that distinct subsets of DMEs were regulated by Hnf4a in a tissue-specific manner. Chromatin immunoprecipitation combined with massively parallel DNA sequencing was performed to characterize the PT-specific binding sites of Hnf4a in rat kidneys at three developmental stages (prenatal, immature, adult), which further supported a major role for Hnf4a in regulating PT gene expression, including DMEs. In ex vivo kidney organ culture, an antagonist of Hnf4a (but not a similar inactive compound) led to predicted changes in DME expression, including among others Fmo1, Cyp2d2, Cyp2d4, Nqo2, as well as organic cation transporters and organic anion transporters Slc22a1 (Oct1), Slc22a2 (Oct2), Slc22a6 (Oat1), Slc22a8 (Oat3), and Slc47a1 (Mate1). Conversely, overexpression of Hnf1a and Hnf4a in primary mouse embryonic fibroblasts, sometimes considered a surrogate for mesenchymal stem cells, induced expression of several of these proximal tubule DMEs, as well as epithelial markers and a PT-enriched brush border marker Ggt1. These cells had organic anion transporter function. Taken together, the data strongly supports a critical role for HNF4a and Hnf1a in the tissue-specific regulation of drug handling and differentiation toward a PT-like cellular identity. We discuss our data in the context of the "remote sensing and signaling hypothesis" (Ahn and Nigam, 2009; Wu et al., 2011).

Concise review: can the intrinsic power of branching morphogenesis be used for engineering epithelial tissues and organs?

Branching morphogenesis is critical to the development of organs such as kidney, lung, mammary gland, prostate, pancreas, and salivary gland. Essentially, an epithelial bud becomes an iterative tip-stalk generator (ITSG) able to form a tree of branching ducts and/or tubules. In different organs, branching morphogenesis is governed by similar sets of genes. Epithelial branching has been recapitulated in vitro (or ex vivo) using three-dimensional cell culture and partial organ culture systems, and several such systems relevant to kidney tissue engineering are discussed here. By adapting systems like these it may be possible to harness the power inherent in the ITSG program to propagate and engineer epithelial tissues and organs. It is also possible to conceive of a universal ITSG capable of propagation that may, by recombination with organ-specific mesenchymal cells, be used for engineering many organ-like tissues similar to the organ from which the mesenchyme cells were derived, or toward which they are differentiated (from stem cells). The three-dimensional (3D) branched epithelial structure could act as a dynamic branching cellular scaffold to establish the architecture for the rest of the tissue. Another strategy-that of recombining propagated organ-specific ITSGs in 3D culture with undifferentiated mesenchymal stem cells-is also worth exploring. If feasible, such engineered tissues may be useful for the ex vivo study of drug toxicity, developmental biology, and physiology in the laboratory. Over the long term, they have potential clinical applications in the general fields of transplantation, regenerative medicine, and bioartificial medical devices to aid in the treatment of chronic kidney disease, diabetes, and other diseases.

Chronic unilateral ureteral obstruction in the neonatal mouse delays maturation of both kidneys and leads to late formation of atubular glomeruli

Unilateral ureteral obstruction (UUO) in the adult mouse is the most widely used model of progressive renal disease: the proximal tubule is the nephron segment most severely affected and atubular glomeruli are formed after only 7 days of UUO. To determine the proximal nephron response to UUO in the maturing kidney, neonatal mice were examined 7 to 28 days following complete UUO under general anesthesia. Proximal tubular mass and maturation were determined by staining with Lotus tetragolonobus lectin. Superoxide was localized by nitroblue tetrazolium and collagen by Sirius red. Cell proliferation, cell death, PAX-2, megalin, α-smooth muscle actin (α-SMA), renin, and fibronectin were identified by immunohistochemistry. During the first 14 days of ipsilateral UUO, despite oxidative stress (4-hydroxynonenal staining), glomerulotubular continuity was maintained and mitochondrial superoxide production persisted. However, from 14 to 28 days, papillary growth was impaired and proximal tubules collapsed with increased apoptosis, autophagy, mitochondrial loss, and formation of atubular glomeruli. Fibronectin, α-SMA, and collagen increased in the obstructed kidney. Oxidative stress was present also in the contralateral kidney: renin was decreased, glomerulotubular maturation and papillary growth were delayed, followed by increased cortical and medullary growth. We conclude that neonatal UUO initially delays renal maturation and results in oxidative stress in both kidneys. In contrast to the adult, proximal tubular injury in the neonatal obstructed kidney is delayed at 14 days, followed only later by the formation of atubular glomeruli. Antioxidant therapies directed at proximal tubular mitochondria during early renal maturation may slow progression of congenital obstructive nephropathy.

Multispecific drug transporter Slc22a8 (Oat3) regulates multiple metabolic and signaling pathways

Multispecific drug transporters of the solute carrier and ATP-binding cassette families are highly conserved through evolution, but their true physiologic role remains unclear. Analyses of the organic anion transporter 3 (OAT3; encoded by Slc22a8/Oat3, originally Roct) knockout mouse have confirmed its critical role in the renal handling of common drugs (e.g., antibiotics, antivirals, diuretics) and toxins. Previous targeted metabolomics of the knockout of the closely related Oat1 have demonstrated a central metabolic role, but the same approach with Oat3 failed to reveal a similar set of endogenous substrates. Nevertheless, the Oat3 knockout is the only Oat described so far with a physiologically significant phenotype, suggesting the disturbance of metabolic or signaling pathways. Here we analyzed global gene expression in Oat3 knockout tissue, which implicated OAT3 in phase I and phase II metabolism (drug metabolizing enzymes or DMEs), as well as signaling pathways. Metabolic reconstruction with the recently developed "mouse Recon1" supported the involvement of Oat3 in the aforementioned pathways. Untargeted metabolomics were used to determine whether the predicted metabolic alterations could be confirmed. Many significant changes were observed; several metabolites were tested for direct interaction with mOAT3, whereas others were supported by published data. Oat3 thus appears critical for the handling of phase I (hydroxylation) and phase II (glucuronidation) metabolites. Oat3 also plays a role in bioenergetic pathways (e.g., the tricarboxylic acid cycle), as well as those involving vitamins (e.g., folate), steroids, prostaglandins, gut microbiome products, uremic toxins, cyclic nucleotides, amino acids, glycans, and possibly hyaluronic acid. The data seemingly consistent with the Remote Sensing and Signaling Hypothesis (Ahn and Nigam, 2009; Wu et al., 2011), also suggests that Oat3 is essential for the handling of dietary flavonoids and antioxidants.

Organic anion and cation SLC22 "drug" transporter (Oat1, Oat3, and Oct1) regulation during development and maturation of the kidney proximal tubule

Proper physiological function in the pre- and post-natal proximal tubule of the kidney depends upon the acquisition of selective permeability, apical-basolateral epithelial polarity and the expression of key transporters, including those involved in metabolite, toxin and drug handling. Particularly important are the SLC22 family of transporters, including the organic anion transporters Oat1 (originally identified as NKT) and Oat3 as well as the organic cation transporter Oct1. In ex vivo cultures of metanephric mesenchyme (MM; the embryonic progenitor tissue of the nephron) Oat function was evident before completion of nephron segmentation and corresponded with the maturation of tight junctions as measured biochemically by detergent extractability of the tight junction protein, ZO-1. Examination of available time series microarray data sets in the context of development and differentiation of the proximal tubule (derived from both in vivo and in vitro/ex vivo developing nephrons) allowed for correlation of gene expression data to biochemically and functionally defined states of development. This bioinformatic analysis yielded a network of genes with connectivity biased toward Hnf4α (but including Hnf1α, hyaluronic acid-CD44, and notch pathways). Intriguingly, the Oat1 and Oat3 genes were found to have strong temporal co-expression with Hnf4α in the cultured MM supporting the notion of some connection between the transporters and this transcription factor. Taken together with the ChIP-qPCR finding that Hnf4α occupies Oat1, Oat3, and Oct1 proximal promoters in the in vivo differentiating rat kidney, the data suggest a network of genes with Hnf4α at its center plays a role in regulating the terminal differentiation and capacity for drug and toxin handling by the nascent proximal tubule of the kidney.

A protein kinase A and Wnt-dependent network regulating an intermediate stage in epithelial tubulogenesis during kidney development

Genetic interactions regulating intermediate stages of tubulogenesis in the developing kidney have been difficult to define. A systems biology strategy using microarray was combined with in vitro/ex vivo and genetic approaches to identify pathways regulating specific stages of tubulogenesis. Analysis of the progression of the metanephric mesenchyme (MM) through four stages of tubule induction and differentiation (i.e., epithelialization, tubular organization and elongation and early differentiation) revealed signaling pathways potentially involved at each stage and suggested key roles for a number of signaling molecules. A screen of the signaling pathways on in vitro/ex vivo nephron formation implicated a unique regulatory role for protein kinase A (PKA), through PKA-2, in a specific post-epithelialization morphogenetic step (conversion of the renal vesicle to the S-shaped body). Microarray analysis not only confirmed this stage-specificity, but also highlighted the upregulation of Wnt genes. Addition of PKA agonists to LIF-induced nephrons (previously shown to be a Wnt/beta-catenin dependent pathway) disrupted normal tubulogenesis in a manner similar to PKA-agonist treated MM/spinal-cord assays, suggesting that PKA regulates a Wnt-dependent tubulogenesis step. PKA induction of canonical Wnt signaling during tubulogenesis was confirmed genetically using MM from Batgal-reporter mice. Addition of a Wnt synthesis inhibitor to activated PKA cultures rescued tubulogenesis. By re-analysis of existing microarray data from the FGF8, Lim1 and Wnt4 knockouts, which arrest in early tubulogenesis, a network of genes involving PKA, Wnt, Lhx1, FGF8, and hyaluronic acid signaling regulating the transition of nascent epithelial cells to tubular epithelium was derived, helping to reconcile in vivo and in vitro/ex vivo data.

Expression cartography of human tissues using self organizing maps

Background

Parallel high-throughput microarray and sequencing experiments produce vast quantities of multidimensional data which must be arranged and analyzed in a concerted way. One approach to addressing this challenge is the machine learning technique known as self organizing maps (SOMs). SOMs enable a parallel sample- and gene-centered view of genomic data combined with strong visualization and second-level analysis capabilities. The paper aims at bridging the gap between the potency of SOM-machine learning to reduce dimension of high-dimensional data on one hand and practical applications with special emphasis on gene expression analysis on the other hand.

Results

The method was applied to generate a SOM characterizing the whole genome expression profiles of 67 healthy human tissues selected from ten tissue categories (adipose, endocrine, homeostasis, digestion, exocrine, epithelium, sexual reproduction, muscle, immune system and nervous tissues). SOM mapping reduces the dimension of expression data from ten of thousands of genes to a few thousand metagenes, each representing a minicluster of co-regulated single genes. Tissue-specific and common properties shared between groups of tissues emerge as a handful of localized spots in the tissue maps collecting groups of co-regulated and co-expressed metagenes. The functional context of the spots was discovered using overrepresentation analysis with respect to pre-defined gene sets of known functional impact. We found that tissue related spots typically contain enriched populations of genes related to specific molecular processes in the respective tissue. Analysis techniques normally used at the gene-level such as two-way hierarchical clustering are better represented and provide better signal-to-noise ratios if applied to the metagenes. Metagene-based clustering analyses aggregate the tissues broadly into three clusters containing nervous, immune system and the remaining tissues.

Conclusions

The SOM technique provides a more intuitive and informative global view of the behavior of a few well-defined modules of correlated and differentially expressed genes than the separate discovery of the expression levels of hundreds or thousands of individual genes. The program is available as R-package 'oposSOM'.

Functional maturation of drug transporters in the developing, neonatal, and postnatal kidney

Because renal function in newborns is immature, the pharmacokinetics of drugs administered to neonates vary significantly from adult patients. The establishment of drug transport systems is a key process in the functional maturation of the nephron. However, a thorough examination of the expression of the main drug transporters in the kidney throughout all stages of development (embryonic, postnatal, and mature) has yet to be carried out, and the functional (physiological) impact is not well understood. Using time-series microarray data, we analyzed the temporal behavior of mRNA levels for a wide range of SLC and ABC transporters in the rodent kidney throughout a developmental time series. We find dynamic increases between the postnatal and mature stages of development for a number of transporters, including the proximal tubule-specific drug and organic anion transporters (OATs) OAT1 (SLC22a6) and OAT3 (SLC22a8). The OATs are the major multispecific basolateral drug, toxin, and metabolite transporters in the proximal tubule responsible for handling of many drugs, as well as the prototypical OAT substrate para-aminohippurate (PAH). We therefore performed specific in vivo pharmacokinetic analysis of the transport of PAH in postnatal and maturing rodent kidney. We show that there is a 4-fold increase in PAH clearance during this period. Clearance studies in Oat1 and Oat3 knockouts confirm that, as in the adult, Oat1 is the principle transporter of PAH in the postnatal kidney. The substantial differences observed supports the need for better understanding of pharmacokinetics in the newborn and juvenile kidney compared with the adult kidney at the basic and clinical level.

Toward a clinically feasible gene expression-based subclassification strategy for septic shock: proof of concept

OBJECTIVE

To develop a clinically feasible stratification strategy for pediatric septic shock, using gene expression mosaics and a 100-gene signature representing the first 24 hrs of admission to the pediatric intensive care unit.

DESIGN

Prospective, observational study involving microarray-based bioinformatics.

SETTING

Multiple pediatric intensive care units in the United States.

PATIENTS

Ninety-eight children with septic shock.

INTERVENTIONS

None other than standard care.

MEASUREMENTS AND MAIN RESULTS

Patients were classified into three previously published, genome-wide, expression-based subclasses (subclasses A, B, and C) having clinically relevant phenotypic differences. The class-defining 100-gene signature was depicted for each individual patient, using mosaics generated by the Gene Expression Dynamics Inspector (GEDI). Composite mosaics were generated representing the average expression patterns for each of the three subclasses. Nine individual clinicians served as blinded evaluators. Each evaluator was shown the 98 individual patient mosaics and asked to classify each patient into one of the three subclasses, using the composite mosaics as the reference point. The respective sensitivities, specificities, positive predictive values, and negative predictive values of the subclassification strategy were ≥ 4% across the three subclasses. The classification strategy also generated positive likelihood ratios of ≥ 6.8 and negative likelihood ratios of ≤ .2 across the three subclasses. The κ coefficient across all possible interevaluator comparisons was 0.81.

CONCLUSIONS

We have provided initial evidence (proof of concept) for a clinically feasible and robust stratification strategy for pediatric septic shock based on a 100-gene signature and gene expression mosaics.

Patterning a complex organ: branching morphogenesis and nephron segmentation in kidney development

The two major components of the kidney, the collecting system and the nephron, have different developmental histories. The collecting system arises by the reiterated branching of a simple epithelial tube, while the nephron forms from a cloud of mesenchymal cells that coalesce into epithelial vesicles. Each develops into a morphologically complex and highly differentiated structure, and together they provide essential filtration and resorption functions. In this review, we will consider their embryological origin and the genes controlling their morphogenesis, patterning, and differentiation, with a focus on recent advances in several areas.

Neuropeptide Y functions as a facilitator of GDNF-induced budding of the Wolffian duct

Ureteric bud (UB) emergence from the Wolffian duct (WD), the initiating step in metanephric kidney morphogenesis, is dependent on GDNF; however, GDNF by itself is generally insufficient to induce robust budding of the isolated WD in culture. Thus, additional factors, presumably peptides or polypeptide growth factors, might be involved. Microarray data from in vivo budding and non-budding conditions were analyzed using non-negative matrix factorization followed by gene ontology filtering and network analysis to identify sets of genes that are highly regulated during budding. These included the GDNF co-receptors GFRalpha1 and RET, as well as neuropeptide Y (NPY). By using ANOVA with pattern matching, NPY was also found to correlate most significantly to the budded condition with a high degree of connectedness to genes with developmental roles. Exogenous NPY [as well as its homolog, peptide YY (PYY)] augmented GDNF-dependent budding in the isolated WD culture; conversely, inhibition of NPY signaling or perturbation of NPY expression inhibited budding, confirming that NPY facilitates this process. NPY was also found to reverse the decreased budding, the downregulation of RET expression, the mislocalization of GFRalpha1, and the inhibition of AKT phosphorylation that resulted from the addition of BMP4 to the isolated WD cultures, suggesting that NPY acts through the budding pathway and is reciprocally regulated by GDNF and BMP4. Thus, the outgrowth of the UB from the WD might result from a combination of the upregulation of the GDNF receptors together with genes that support GDNF signaling in a feed-forward loop and/or counteraction of the inhibitory pathway regulated by BMP4.

Hs2st mediated kidney mesenchyme induction regulates early ureteric bud branching

Heparan sulfate proteoglycans (HSPGs) are central modulators of developmental processes likely through their interaction with growth factors, such as GDNF, members of the FGF and TGFbeta superfamilies, EGF receptor ligands and HGF. Absence of the biosynthetic enzyme, heparan sulfate 2-O-sulfotransferase (Hs2st) leads to kidney agenesis. Using a novel combination of in vivo and in vitro approaches, we have reanalyzed the defect in morphogenesis of the Hs2st(-)(/)(-) kidney. Utilizing assays that separately model distinct stages of kidney branching morphogenesis, we found that the Hs2st(-/-) UB is able to undergo branching and induce mesenchymal-to-epithelial transformation when recombined with control MM, and the isolated Hs2st null UB is able to undergo branching morphogenesis in the presence of exogenous soluble pro-branching growth factors when embedded in an extracellular matrix, indicating that the UB is intrinsically competent. This is in contrast to the prevailing view that the defect underlying the renal agenesis phenotype is due to a primary role for 2-O sulfated HS in UB branching. Unexpectedly, the mutant MM was also fully capable of being induced in recombination experiments with wild-type tissue. Thus, both the mutant UB and mutant MM tissue appear competent in and of themselves, but the combination of mutant tissues fails in vivo and, as we show, in organ culture. We hypothesized a 2OS-dependent defect in the mutual inductive process, which could be on either the UB or MM side, since both progenitor tissues express Hs2st. In light of these observations, we specifically examined the role of the HS 2-O sulfation modification on the morphogenetic capacity of the UB and MM individually. We demonstrate that early UB branching morphogenesis is not primarily modulated by factors that depend on the HS 2-O sulfate modification; however, factors that contribute to MM induction are markedly sensitive to the 2-O sulfation modification. These data suggest that key defect in Hs2st null kidneys is the inability of MM to undergo induction either through a failure of mutual induction or a primary failure of MM morphogenesis. This results in normal UB formation but affects either T-shaped UB formation or iterative branching of the T-shaped UB (possibly two separate stages in collecting system development dependent upon HS). We discuss the possibility that a disruption in the interaction between HS and Wnts (e.g. Wnt 9b) may be an important aspect of the observed phenotype. This appears to be the first example of a defect in the MM preventing advancement of early UB branching past the first bifurcation stage, one of the limiting steps in early kidney development.

The instructive role of metanephric mesenchyme in ureteric bud patterning, sculpting, and maturation and its potential ability to buffer ureteric bud branching defects

Kidney organogenesis depends on reciprocal interactions between the ureteric bud (UB) and the metanephric mesenchyme (MM) to form the UB-derived collecting system and MM-derived nephron. With the advent of in vitro systems, it is clear that UB branching can occur independently of MM contact; however, little has been done to detail the role of MM cellular contact in this process. Here, a model system in which the cultured isolated UB is recombined with uninduced MM is used to isolate the effects of the MM progenitor tissue on the development and maturation of the collecting system. By morphometrics, we demonstrate that cellular contact with the MM is required for vectorial elongation of stalks and tapering of luminal caliber of UB-derived tubules. Expression analysis of developmentally significant genes indicates the cocultured tissue is most similar to an embryonic day 19 (E19) kidney. The likely major contributor to this is the functional maturation of the collecting duct and proximal nephron segments in the UB-induced MM, as measured by quantitative PCR, of the collecting duct-specific arginine vasopressin receptor and the nephron tubule segment-specific organic anion transporter OAT1, Na-P(i) type 2 cotransporter, and Tamm-Horsfall protein gene expressions. However, expression of aquaporin-2 is upregulated similarly in isolated UB and cocultured tissue, suggesting that some aspects of functional maturation can occur independently of MM cellular contact. In addition to its sculpting effects, the MM normalized a "branchless" UB morphology induced by FGF7 or heregulin in isolated UB culture. The morphological changes induced by the MM were accompanied by a reassignment of GFRalpha1 (a receptor for GDNF) to tips. Such "quality control" by the MM of UB morphology may provide resiliency to the branching program. This may help to explain a number of knockout phenotypes in which branching and/or cystic defects are less impressive than expected. A second hit in the MM may thus be necessary to make these defects fully apparent.