Kinetic modeling of stem cell transcriptome dynamics to identify regulatory modules of normal and disturbed neuroectodermal differentiation

Thousands of transcriptome data sets are available, but approaches for their use in dynamic cell response modelling are few, especially for processes affected simultaneously by two orthogonal influencing variables. We approached this problem for neuroepithelial development of human pluripotent stem cells (differentiation variable), in the presence or absence of valproic acid (signaling variable). Using few basic assumptions (sequential differentiation states of cells; discrete on/off states for individual genes in these states), and time-resolved transcriptome data, a comprehensive model of spontaneous and perturbed gene expression dynamics was developed. The model made reliable predictions (average correlation of 0.85 between predicted and subsequently tested expression values). Even regulations predicted to be non-monotonic were successfully validated by PCR in new sets of experiments. Transient patterns of gene regulation were identified from model predictions. They pointed towards activation of Wnt signaling as a candidate pathway leading to a redirection of differentiation away from neuroepithelial cells towards neural crest. Intervention experiments, using a Wnt/beta-catenin antagonist, led to a phenotypic rescue of this disturbed differentiation. Thus, our broadly applicable model allows the analysis of transcriptome changes in complex time/perturbation matrices

Radiation Response of Murine Embryonic Stem Cells

To understand the mechanisms of disturbed differentiation and development by radiation, murine CGR8 embryonic stem cells (mESCs) were exposed to ionizing radiation and differentiated by forming embryoid bodies (EBs). The colony forming ability test was applied for survival and the MTT test for viability determination after X-irradiation. Cell cycle progression was determined by flow cytometry of propidium iodide-stained cells, and DNA double strand break (DSB) induction and repair by γH2AX immunofluorescence. The radiosensitivity of mESCs was slightly higher compared to the murine osteoblast cell line OCT-1. The viability 72 h after X-irradiation decreased dose-dependently and was higher in the presence of leukemia inhibitory factor (LIF). Cells exposed to 2 or 7 Gy underwent a transient G2 arrest. X-irradiation induced γH2AX foci and they disappeared within 72 h. After 72 h of X-ray exposure, RNA was isolated and analyzed using genome-wide microarrays. The gene expression analysis revealed amongst others a regulation of developmental genes (Ada, Baz1a, Calcoco2, Htra1, Nefh, S100a6 and Rassf6), downregulation of genes involved in glycolysis and pyruvate metabolism whereas upregulation of genes related to the p53 signaling pathway. X-irradiated mESCs formed EBs and differentiated toward cardiomyocytes but their beating frequencies were lower compared to EBs from unirradiated cells. These results suggest that X-irradiation of mESCs deregulate genes related to the developmental process. The most significant biological processes found to be altered by X-irradiation in mESCs were the development of cardiovascular, nervous, circulatory and renal system. These results may explain the X-irradiation induced-embryonic lethality and malformations observed in animal studies

Transcriptomic and Functional Evidence That miRNA193a-3p Inhibits Lymphatic Endothelial Cell (LEC) and LEC + MCF-7 Spheroid Growth Directly and by Altering MCF-7 Secretome

MicroRNA 193a-3p (miR193a-3p) is a short non-coding RNA with tumor suppressor properties. Breast cancer (BC) progression is governed by active interaction between breast cancer cells, vascular (V)/lymphatic (L) endothelial cells (ECs), and BC secretome. We have recently shown that miR193a-3p, a tumor suppressor miRNA, inhibits MCF-7 BC cell-driven growth of VECs via direct antimitogenic actions and alters MCF-7 secretome. Since LEC-BC cross-talk plays a key role in BC progression, we investigated the effects of miR193a-3p on MCF-7 secretome and estradiol-mediated growth effects in LECs and LEC + MCF-7 spheroids, and delineated the underlying mechanisms. Transfection of LECs with miR193a-3p, as well as secretome from MCF-7 transfected cells, inhibited LEC growth, and these effects were mimicked in LEC + MCF-7 spheroids. Moreover, miR193a-3p inhibited ERK1/2 and Akt phosphorylation in LECs and LEC + MCF-7 spheroids, which are importantly involved in promoting cancer development and metastasis. Treatment of LECs and LEC + MCF-7 spheroids with estradiol (E2)-induced growth, as well as ERK1/2 and Akt phosphorylation, and was abrogated by miR193a-3p and secretome from MCF-7 transfected cells. Gene expression analysis (GEA) in LEC + MCF-7 spheroids transfected with miR193a-3p showed significant upregulation of 54 genes and downregulation of 73 genes. Pathway enrichment analysis of regulated genes showed significant modulation of several pathways, including interferon, interleukin/cytokine-mediated signaling, innate immune system, ERK1/2 cascade, apoptosis, and estrogen receptor signaling. Transcriptomic analysis showed downregulation in interferon and anti-apoptotic and pro-growth molecules, such as IFI6, IFIT1, OSA1/2, IFITM1, HLA-A/B, PSMB8/9, and PARP9, which are known to regulate BC progression. The cytokine proteome array of miR193a-3p transfected MCF secretome and confirmed the upregulation of several growth inhibitory cytokines, including IFNγ, Il-1a, IL-1ra, IL-32, IL-33, IL-24, IL-27, cystatin, C-reactive protein, Fas ligand, MIG, and sTIM3. Moreover, miR193a-3p alters factors in MCF-7 secretome, which represses ERK1/2 and Akt phosphorylation, induces pro-apoptotic protein and apoptosis in LECs, and downregulates interferon-associated proteins known to promote cancer growth and metastasis. In conclusion, miR193a-3p can potentially modify the tumor microenvironment by altering pro-growth BC secretome and inhibiting LEC growth, and may represent a therapeutic molecule to target breast tumors/cancer

Transcriptomic and Functional Evidence That miRNA193a-3p Inhibits Lymphatic Endothelial Cell (LEC) and LEC + MCF-7 Spheroid Growth Directly and by Altering MCF-7 Secretome

MicroRNA 193a-3p (miR193a-3p) is a short non-coding RNA with tumor suppressor properties. Breast cancer (BC) progression is governed by active interaction between breast cancer cells, vascular (V)/lymphatic (L) endothelial cells (ECs), and BC secretome. We have recently shown that miR193a-3p, a tumor suppressor miRNA, inhibits MCF-7 BC cell-driven growth of VECs via direct antimitogenic actions and alters MCF-7 secretome. Since LEC-BC cross-talk plays a key role in BC progression, we investigated the effects of miR193a-3p on MCF-7 secretome and estradiol-mediated growth effects in LECs and LEC + MCF-7 spheroids, and delineated the underlying mechanisms. Transfection of LECs with miR193a-3p, as well as secretome from MCF-7 transfected cells, inhibited LEC growth, and these effects were mimicked in LEC + MCF-7 spheroids. Moreover, miR193a-3p inhibited ERK1/2 and Akt phosphorylation in LECs and LEC + MCF-7 spheroids, which are importantly involved in promoting cancer development and metastasis. Treatment of LECs and LEC + MCF-7 spheroids with estradiol (E2)-induced growth, as well as ERK1/2 and Akt phosphorylation, and was abrogated by miR193a-3p and secretome from MCF-7 transfected cells. Gene expression analysis (GEA) in LEC + MCF-7 spheroids transfected with miR193a-3p showed significant upregulation of 54 genes and downregulation of 73 genes. Pathway enrichment analysis of regulated genes showed significant modulation of several pathways, including interferon, interleukin/cytokine-mediated signaling, innate immune system, ERK1/2 cascade, apoptosis, and estrogen receptor signaling. Transcriptomic analysis showed downregulation in interferon and anti-apoptotic and pro-growth molecules, such as IFI6, IFIT1, OSA1/2, IFITM1, HLA-A/B, PSMB8/9, and PARP9, which are known to regulate BC progression. The cytokine proteome array of miR193a-3p transfected MCF secretome and confirmed the upregulation of several growth inhibitory cytokines, including IFNγ, Il-1a, IL-1ra, IL-32, IL-33, IL-24, IL-27, cystatin, C-reactive protein, Fas ligand, MIG, and sTIM3. Moreover, miR193a-3p alters factors in MCF-7 secretome, which represses ERK1/2 and Akt phosphorylation, induces pro-apoptotic protein and apoptosis in LECs, and downregulates interferon-associated proteins known to promote cancer growth and metastasis. In conclusion, miR193a-3p can potentially modify the tumor microenvironment by altering pro-growth BC secretome and inhibiting LEC growth, and may represent a therapeutic molecule to target breast tumors/cancer

Transcriptomic and Functional Evidence That miRNA193a-3p Inhibits Lymphatic Endothelial Cell (LEC) and LEC + MCF-7 Spheroid Growth Directly and by Altering MCF-7 Secretome

MicroRNA 193a-3p (miR193a-3p) is a short non-coding RNA with tumor suppressor properties. Breast cancer (BC) progression is governed by active interaction between breast cancer cells, vascular (V)/lymphatic (L) endothelial cells (ECs), and BC secretome. We have recently shown that miR193a-3p, a tumor suppressor miRNA, inhibits MCF-7 BC cell-driven growth of VECs via direct antimitogenic actions and alters MCF-7 secretome. Since LEC-BC cross-talk plays a key role in BC progression, we investigated the effects of miR193a-3p on MCF-7 secretome and estradiol-mediated growth effects in LECs and LEC + MCF-7 spheroids, and delineated the underlying mechanisms. Transfection of LECs with miR193a-3p, as well as secretome from MCF-7 transfected cells, inhibited LEC growth, and these effects were mimicked in LEC + MCF-7 spheroids. Moreover, miR193a-3p inhibited ERK1/2 and Akt phosphorylation in LECs and LEC + MCF-7 spheroids, which are importantly involved in promoting cancer development and metastasis. Treatment of LECs and LEC + MCF-7 spheroids with estradiol (E2)-induced growth, as well as ERK1/2 and Akt phosphorylation, and was abrogated by miR193a-3p and secretome from MCF-7 transfected cells. Gene expression analysis (GEA) in LEC + MCF-7 spheroids transfected with miR193a-3p showed significant upregulation of 54 genes and downregulation of 73 genes. Pathway enrichment analysis of regulated genes showed significant modulation of several pathways, including interferon, interleukin/cytokine-mediated signaling, innate immune system, ERK1/2 cascade, apoptosis, and estrogen receptor signaling. Transcriptomic analysis showed downregulation in interferon and anti-apoptotic and pro-growth molecules, such as IFI6, IFIT1, OSA1/2, IFITM1, HLA-A/B, PSMB8/9, and PARP9, which are known to regulate BC progression. The cytokine proteome array of miR193a-3p transfected MCF secretome and confirmed the upregulation of several growth inhibitory cytokines, including IFNγ, Il-1a, IL-1ra, IL-32, IL-33, IL-24, IL-27, cystatin, C-reactive protein, Fas ligand, MIG, and sTIM3. Moreover, miR193a-3p alters factors in MCF-7 secretome, which represses ERK1/2 and Akt phosphorylation, induces pro-apoptotic protein and apoptosis in LECs, and downregulates interferon-associated proteins known to promote cancer growth and metastasis. In conclusion, miR193a-3p can potentially modify the tumor microenvironment by altering pro-growth BC secretome and inhibiting LEC growth, and may represent a therapeutic molecule to target breast tumors/cancer

Differentiation and Selection of Hepatocyte Precursors in Suspension Spheroid Culture of Transgenic Murine Embryonic Stem Cells

Embryonic stem cell-derived hepatocyte precursor cells represent a promising model for clinical transplantations to diseased livers, as well as for establishment of in vitro systems for drug metabolism and toxicology investigations. This study aimed to establish an in vitro culture system for scalable generation of hepatic progenitor cells. We used stable transgenic clones of murine embryonic stem cells possessing a reporter/selection vector, in which the enhanced green fluorescent protein-and puromycin N-acetyltransferase-coding genes are driven by a common alpha-fetoprotein gene promoter. This allowed for live monitoring and puromycin selection of the desired differentiating cell type possessing the activated alpha-fetoprotein gene. A rotary culture system was established, sequentially yielding initially partially selected hepatocyte lineage-committed cells, and finally, a highly purified cell population maintained as a dynamic suspension spheroid culture, which progressively developed the hepatic gene expression phenotype. The latter was confirmed by quantitative RT-PCR analysis, which showed a progressive up-regulation of hepatic genes during spheroid culture, indicating development of a mixed hepatocyte precursor-/fetal hepatocyte-like cell population. Adherent spheroids gave rise to advanced differentiated hepatocyte-like cells expressing hepatic proteins such as albumin, alpha-1-antitrypsin, cytokeratin 18, E-cadherin, and liver-specific organic anion transporter 1, as demonstrated by fluorescent immunostaining. A fraction of adherent cells was capable of glycogen storage and of reversible up-take of indocyanine green, demonstrating their hepatocyte-like functionality. Moreover, after transplantation of spheroids into the mouse liver, the spheroid-derived cells integrated into recipient. These results demonstrate that large-scale hepatocyte precursor-/hepatocyte-like cultures can be established for use in clinical trials, as well as in in vitro screening assays

Distinct transcriptomic changes in E14.5 mouse skeletal muscle lacking RYR1 or Ca(v)1.1 converge at E18.5

In skeletal muscle the coordinated actions of two mechanically coupled Ca2+ channels-the 1,4-dihydropyridine receptor (Ca(v)1.1) and the type 1 ryanodine receptor (RYR1)-underlie the molecular mechanism of rapid cytosolic [Ca2+] increase leading to contraction. While both [Ca2+](i) and contractile activity have been implicated in the regulation of myogenesis, less is known about potential specific roles of Ca(v)1.1 and RYR1 in skeletal muscle development. In this study, we analyzed the histology and the transcriptomic changes occurring at E14.5 -the end of primary myogenesis and around the onset of intrauterine limb movement, and at E18.5 -the end of secondary myogenesis, in WT, RYR1(-/-), and Ca(v)1.1(-/-) murine limb skeletal muscle. At E14.5 the muscle histology of both mutants exhibited initial alterations, which became much more severe at E18.5. Immunohistological analysis also revealed higher levels of activated caspase-3 in the Cav1.1(-/-) muscles at E14.5, indicating an increase in apoptosis. With WT littermates as controls, microarray analyses identified 61 and 97 differentially regulated genes (DEGs) at E14.5, and 493 and 1047 DEGs at E18.5, in RYR1(-/-) and Cav1.1(-/-) samples, respectively. Gene enrichment analysis detected no overlap in the affected biological processes and pathways in the two mutants at E14.5, whereas at E18.5 there was a significant overlap of DEGs in both mutants, affecting predominantly processes linked to muscle contraction. Moreover, the E18.5 vs. E14.5 comparison revealed multiple genotype-specific DEGs involved in contraction, cell cycle and miRNA-mediated signaling in WT, neuronal and bone development in RYR1(-/-), and lipid metabolism in Ca(v)1.1(-/-) samples. Taken together, our study reveals discrete changes in the global transcriptome occurring in limb skeletal muscle from E14.5 to E18.5 in WT, RYR1(-/-) and Ca(v)1.1(-/-) mice. Our results suggest distinct functional roles for RYR1 and Cav1.1 in skeletal primary and secondary myogenesis

Modulation of Differentiation Processes in Murine Embryonic Stem Cells Exposed to Parabolic Flight-Induced Acute Hypergravity and Microgravity

Embryonic developmental studies under microgravity conditions in space are very limited. To study the effects of short-term altered gravity on embryonic development processes, we exposed mouse embryonic stem cells (mESCs) to phases of hypergravity and microgravity and studied the differentiation potential of the cells using wide-genome microarray analysis. During the 64th European Space Agency's parabolic flight campaign, mESCs were exposed to 31 parabolas. Each parabola comprised phases lasting 22 s of hypergravity, microgravity, and a repeat of hypergravity. On different parabolas, RNA was isolated for microarray analysis. After exposure to 31 parabolas, mESCs (P31 mESCs) were further differentiated under normal gravity (1 g) conditions for 12 days, producing P31 12-day embryoid bodies (EBs). After analysis of the microarrays, the differentially expressed genes were analyzed using different bioinformatic tools to identify developmental and nondevelopmental biological processes affected by conditions on the parabolic flight experiment. Our results demonstrated that several genes belonging to GOs associated with cell cycle and proliferation were downregulated in undifferentiated mESCs exposed to gravity changes. However, several genes belonging to developmental processes, such as vasculature development, kidney development, skin development, and to the TGF-β signaling pathway, were upregulated. Interestingly, similar enriched and suppressed GOs were obtained in P31 12-day EBs compared with ground control 12-day EBs. Our results show that undifferentiated mESCs exposed to alternate hypergravity and microgravity phases expressed several genes associated with developmental/differentiation and cell cycle processes, suggesting a transition from the undifferentiated pluripotent to a more differentiated stage of mESCs

Stem Cell-Derived Immature Human Dorsal Root Ganglia Neurons to Identify Peripheral Neurotoxicants

Safety sciences and the identification of chemical hazards have been seen as one of the most immediate practical applications of human pluripotent stem cell technology. Protocols for the generation of many desirable human cell types have been developed, but optimization of neuronal models for toxicological use has been astonishingly slow, and the wide, clinically important field of peripheral neurotoxicity is still largely unexplored. A two-step protocol to generate large lots of identical peripheral human neuronal precursors was characterized and adapted to the measurement of peripheral neurotoxicity. High content imaging allowed an unbiased assessment of cell morphology and viability. The computational quantification of neurite growth as a functional parameter highly sensitive to disturbances by toxicants was used as an endpoint reflecting specific neurotoxicity. The differentiation of cells toward dorsal root ganglia neurons was tracked in relation to a large background data set based on gene expression microarrays. On this basis, a peripheral neurotoxicity (PeriTox) test was developed as a first toxicological assay that harnesses the potential of human pluripotent stem cells to generate cell types/tissues that are not otherwise available for the prediction of human systemic organ toxicity. Testing of more than 30 chemicals showed that human neurotoxicants and neurite growth enhancers were correctly identified. Various classes of chemotherapeutic agents causing human peripheral neuropathies were identified, and they were missed when tested on human central neurons. The PeriTox test we established shows the potential of human stem cells for clinically relevant safety testing of drugs in use and of new emerging candidates