Retinoic Acid Fluctuation Activates an Uneven, Direction-Dependent Network-Wide Robustness Response in Early Embryogenesis.

Robustness is a feature of regulatory pathways to ensure signal consistency in light of environmental changes or genetic polymorphisms. The retinoic acid (RA) pathway, is a central developmental and tissue homeostasis regulatory signal, strongly dependent on nutritional sources of retinoids and affected by environmental chemicals. This pathway is characterized by multiple proteins or enzymes capable of performing each step and their integration into a self-regulating network. We studied RA network robustness by transient physiological RA signaling disturbances followed by kinetic transcriptomic analysis of the recovery during embryogenesis. The RA metabolic network was identified as the main regulated module to achieve signaling robustness using an unbiased pattern analysis. We describe the network-wide responses to RA signal manipulation and found the feedback autoregulation to be sensitive to the direction of the RA perturbation: RA knockdown exhibited an upper response limit, whereas RA addition had a minimal feedback-activation threshold. Surprisingly, our robustness response analysis suggests that the RA metabolic network regulation exhibits a multi-objective optimization, known as Pareto optimization, characterized by trade-offs between competing functionalities. We observe that efficient robustness to increasing RA is accompanied by worsening robustness to reduced RA levels and vice versa. This direction-dependent trade-off in the network-wide feedback response, results in an uneven robustness capacity of the RA network during early embryogenesis, likely a significant contributor to the manifestation of developmental defects.

A single cell transcriptomics map of paracrine networks in the intrinsic cardiac nervous system.

We developed a spatially-tracked single neuron transcriptomics map of an intrinsic cardiac ganglion, the right atrial ganglionic plexus (RAGP) that is a critical mediator of sinoatrial node (SAN) activity. This 3D representation of RAGP used neuronal tracing to extensively map the spatial distribution of the subset of neurons that project to the SAN. RNA-seq of laser capture microdissected neurons revealed a distinct composition of RAGP neurons compared to the central nervous system and a surprising finding that cholinergic and catecholaminergic markers are coexpressed, suggesting multipotential phenotypes that can drive neuroplasticity within RAGP. High-throughput qPCR of hundreds of laser capture microdissected single neurons confirmed these findings and revealed a high dimensionality of neuromodulatory factors that contribute to dynamic control of the heart. Neuropeptide-receptor coexpression analysis revealed a combinatorial paracrine neuromodulatory network within RAGP informing follow-on studies on the vagal control of RAGP to regulate cardiac function in health and disease.

Innervation and Neuronal Control of the Mammalian Sinoatrial Node a Comprehensive Atlas.

[Figure: see text].

A Comprehensive Integrated Anatomical and Molecular Atlas of Rat Intrinsic Cardiac Nervous System.

We have developed and integrated several technologies including whole-organ imaging and software development to support an initial precise 3D neuroanatomical mapping and molecular phenotyping of the intracardiac nervous system (ICN). While qualitative and gross anatomical descriptions of the anatomy of the ICN have each been pursued, we here bring forth a comprehensive atlas of the entire rat ICN at single-cell resolution. Our work precisely integrates anatomical and molecular data in the 3D digitally reconstructed whole heart with resolution at the micron scale. We now display the full extent and the position of neuronal clusters on the base and posterior left atrium of the rat heart, and the distribution of molecular phenotypes that are defined along the base-to-apex axis, which had not been previously described. The development of these approaches needed for this work has produced method pipelines that provide the means for mapping other organs.

Single-Cell Gene Expression Analysis Identifies Chronic Alcohol-Mediated Shift in Hepatocyte Molecular States After Partial Hepatectomy.

The analysis of molecular states of individual cells, as defined by their mRNA expression profiles and protein composition, has gained widespread interest in studying biological phenomena ranging from embryonic development to homeostatic tissue function and genesis and evolution of cancers. Although the molecular content of individual cells in a tissue can vary widely, their molecular states tend to be constrained within a transcriptional landscape partly described by the canonical archetypes of a population of cells. In this study, we sought to characterize the effects of an acute (partial hepatectomy) and chronic (alcohol consumption) perturbation on the molecular states of individual hepatocytes during the onset and progression of liver regeneration. We analyzed the expression of 84 genes across 233 individual hepatocytes acquired using laser capture microdissection. Analysis of the single-cell data revealed that hepatocyte molecular states can be considered as distributed across a set of four states irrespective of perturbation, with the proportions of hepatocytes in these states being dependent on the perturbation. In addition to the quiescent, primed, and replicating hepatocytes, we identified a fourth molecular state lying between the primed and replicating subpopulations. Comparison of the proportions of hepatocytes from each experimental condition in these four molecular states suggested that, in addition to aberrant priming, a slower transition from primed to replication state could contribute toward ethanol-mediated suppression of liver regenerative response to partial hepatectomy.

Cellular network modeling and single cell gene expression analysis reveals novel hepatic stellate cell phenotypes controlling liver regeneration dynamics.

BACKGROUND: Recent results from single cell gene and protein regulation studies are starting to uncover the previously underappreciated fact that individual cells within a population exhibit high variability in the expression of mRNA and proteins (i.e., molecular variability). By combining cellular network modeling, and high-throughput gene expression measurements in single cells, we seek to reconcile the high molecular variability in single cells with the relatively low variability in tissue-scale gene and protein expression and the highly coordinated functional responses of tissues to physiological challenges. In this study, we focus on relating the dynamic changes in distributions of hepatic stellate cell (HSC) functional phenotypes to the tightly regulated physiological response of liver regeneration. RESULTS: We develop a mathematical model describing contributions of HSC functional phenotype populations to liver regeneration and test model predictions through isolation and transcriptional characterization of single HSCs. We identify and characterize four HSC transcriptional states contributing to liver regeneration, two of which are described for the first time in this work. We show that HSC state populations change in vivo in response to acute challenges (in this case, 70% partial hepatectomy) and chronic challenges (chronic ethanol consumption). Our results indicate that HSCs influence the dynamics of liver regeneration through steady-state tissue preconditioning prior to an acute insult and through dynamic control of cell state balances. Furthermore, our modeling approach provides a framework to understand how balances among cell states influence tissue dynamics. CONCLUSIONS: Taken together, our combined modeling and experimental studies reveal novel HSC transcriptional states and indicate that baseline differences in HSC phenotypes as well as a dynamic balance of transitions between these phenotypes control liver regeneration responses.

Integrated live imaging and molecular profiling of embryoid bodies reveals a synchronized progression of early differentiation.

Embryonic stem cells can spontaneously differentiate into cell types of all germ layers within embryoid bodies (EBs) in a highly variable manner. Whether there exists an intrinsic differentiation program common to all EBs is unknown. Here, we present a novel combination of high-throughput live two-photon imaging and gene expression profiling to study early differentiation dynamics spontaneously occurring within developing EBs. Onset timing of Brachyury-GFP was highly variable across EBs, while the spatial patterns as well as the dynamics of mesendodermal progression following onset were remarkably similar. We therefore defined a 'developmental clock' using the Brachyury-GFP signal onset timing. Mapping snapshot gene expression measurements to this clock revealed their temporal trends, indicating that loss of pluripotency, formation of primitive streak and mesodermal lineage progression are synchronized in EBs. Exogenous activation of Wnt or BMP signaling accelerated the intrinsic clock. CHIR down-regulated Wnt3, allowing insights into dependency mechanisms between canonical Wnt signaling and multiple genes. Our findings reveal a developmental clock characteristic of an early differentiation program common to all EBs, further establishing them as an in vitro developmental model.

Prolactin-stimulated activation of ERK1/2 mitogen-activated protein kinases is controlled by PI3-kinase/Rac/PAK signaling pathway in breast cancer cells.

There is strong evidence that deregulation of prolactin (PRL) signaling contributes to pathogenesis and chemoresistance of breast cancer. Therefore, understanding cross-talk between distinct signal transduction pathways triggered by activation of the prolactin receptor (PRL-R), is essential for elucidating the pathogenesis of metastatic breast cancer. In this study, we applied a sequential inhibitory analysis of various signaling intermediates to examine the hierarchy of protein interactions within the PRL signaling network and to evaluate the relative contributions of multiple signaling branches downstream of PRL-R to the activation of the extracellular signal-regulated kinases ERK1 and ERK2 in T47D and MCF-7 human breast cancer cells. Quantitative measurements of the phosphorylation/activation patterns of proteins showed that PRL simultaneously activated Src family kinases (SFKs) and the JAK/STAT, phosphoinositide-3 (PI3)-kinase/Akt and MAPK signaling pathways. The specific blockade or siRNA-mediated suppression of SFK/FAK, JAK2/STAT5, PI3-kinase/PDK1/Akt, Rac/PAK or Ras regulatory circuits revealed that (1) the PI3-kinase/Akt pathway is required for activation of the MAPK/ERK signaling cascade upon PRL stimulation; (2) PI3-kinase-mediated activation of the c-Raf-MEK1/2-ERK1/2 cascade occurs independent of signaling dowstream of STATs, Akt and PKC, but requires JAK2, SFKs and FAK activities; (3) activated PRL-R mainly utilizes the PI3-kinase-dependent Rac/PAK pathway rather than the canonical Shc/Grb2/SOS/Ras route to initiate and sustain ERK1/2 signaling. By interconnecting diverse signaling pathways PLR may enhance proliferation, survival, migration and invasiveness of breast cancer cells.