Pluripotent stem cells are able to self-renew, and to differentiate into all adult cell types. Many studies report data describing these cells, and characterize them in molecular terms. Machine learning yields classifiers that can accurately identify pluripotent stem cells, but there is a lack of studies yielding minimal sets of best biomarkers (genes/features). We assembled gene expression data of pluripotent stem cells and non-pluripotent cells from the mouse. After normalization and filtering, we applied machine learning, classifying samples into pluripotent and non-pluripotent with high cross-validated accuracy. Furthermore, to identify minimal sets of best biomarkers, we used three methods: information gain, random forests and a wrapper of genetic algorithm and support vector machine (GA/SVM). We demonstrate that the GA/SVM biomarkers work best in combination with each other; pathway and enrichment analyses show that they cover the widest variety of processes implicated in pluripotency. The GA/SVM wrapper yields best biomarkers, no matter which classification method is used. The consensus best biomarker based on the three methods is Tet1, implicated in pluripotency just recently. The best biomarker based on the GA/SVM wrapper approach alone is Fam134b, possibly a missing link between pluripotency and some standard surface markers of unknown function processed by the Golgi apparatus.

Much of embryonic stem cell biology has focused on transcriptional expression and regulation of genes that could mediate its unique potential in self-renewal or pluripotency. In alignment with our present understanding on the genetic, protein, and epigenetic factors that may direct cell fate, we present a short overview of the often overlooked contribution of alternative splice variants to regulatory diversity. Progressing beyond the limitations of a fixed genomic sequence, alternative splicing offers an additional layer of complexity to produce protein variants that may differ in function and localization that can direct embryonic stem cells to specific differentiation pathways. In light of the number of variants that can be produced at key ES cell genes alone, it is challenging to consider how much more multifaceted transcriptional regulation truly is, and if this can be captured more fully in future works.

The role of alternative splicing in self-renewal, pluripotency and tissue lineage specification of human embryonic stem cells (hESCs) is largely unknown. To better define these regulatory cues, we modified the H9 hESC line to allow selection of pluripotent hESCs by neomycin resistance and cardiac progenitors by puromycin resistance. Exon-level microarray expression data from undifferentiated hESCs and cardiac and neural precursors were used to identify splice isoforms with cardiac-restricted or common cardiac/neural differentiation expression patterns. Splice events for these groups corresponded to the pathways of cytoskeletal remodeling, RNA splicing, muscle specification, and cell cycle checkpoint control as well as genes with serine/threonine kinase and helicase activity. Using a new program named AltAnalyze (http://www.AltAnalyze.org), we identified novel changes in protein domain and microRNA binding site architecture that were predicted to affect protein function and expression. These included an enrichment of splice isoforms that oppose cell-cycle arrest in hESCs and that promote calcium signaling and cardiac development in cardiac precursors. By combining genome-wide predictions of alternative splicing with new functional annotations, our data suggest potential mechanisms that may influence lineage commitment and hESC maintenance at the level of specific splice isoforms and microRNA regulation.

Alternative splicing contributes to both gene regulation and protein diversity. To discover broad relationships between regulation of alternative splicing and sequence conservation, we applied a systems approach, using oligonucleotide microarrays designed to capture splicing information across the mouse genome. In a set of 22 adult tissues, we observe differential expression of RNA containing at least two alternative splice junctions for about 40% of the 6,216 alternative events we could detect. Statistical comparisons identify 171 cassette exons whose inclusion or skipping is different in brain relative to other tissues and another 28 exons whose splicing is different in muscle. A subset of these exons is associated with unusual blocks of intron sequence whose conservation in vertebrates rivals that of protein-coding exons. By focusing on sets of exons with similar regulatory patterns, we have identified new sequence motifs implicated in brain and muscle splicing regulation. Of note is a motif that is strikingly similar to the branchpoint consensus but is located downstream of the 5' splice site of exons included in muscle. Analysis of three paralogous membrane-associated guanylate kinase genes reveals that each contains a paralogous tissue-regulated exon with a similar tissue inclusion pattern. While the intron sequences flanking these exons remain highly conserved among mammalian orthologs, the paralogous flanking intron sequences have diverged considerably, suggesting unusually complex evolution of the regulation of alternative splicing in multigene families.

MOTIVATION: Many or most mammalian genes undergo alternative splicing, generating a variety of transcripts from a single gene. New information on splice variation is becoming available through technology for measuring expression levels of several exons or splice junctions per gene. We have developed a statistical method, ANalysis Of Splice VAriation (ANOSVA) to detect alternative splicing from expression data. Since ANOSVA requires no transcript information, it can be applied when the level of annotation is poor. When validated against spiked clone data, it generated no false positives and few false negatives. We demonstrated ANOSVA with data from a prototype mouse alternative splicing array, run against normal adult tissues, yielding a set of genes with evidence of tissue-specific splice variation. AVAILABILITY: The results are available at the supplementary information site. CONTACT: cline@pasteur.fr SUPPLEMENTARY INFORMATION: The results are available at the supplementary information site https://bioinfo.affymetrix.com/Papers/ANOSVA/

Although a substantial number of hormones and drugs increase cellular cAMP levels, the global impact of cAMP and its major effector mechanism, protein kinase A (PKA), on gene expression is not known. Here we show that treatment of murine wild-type S49 lymphoma cells for 24 h with 8-(4-chlorophenylthio)-cAMP (8-CPT-cAMP), a PKA-selective cAMP analog, alters the expression of approximately 4,500 of approximately 13,600 unique genes. By contrast, gene expression was unaltered in Kin- S49 cells (that lack PKA) incubated with 8-CPT-cAMP. Changes in mRNA and protein expression of several cell-cycle regulators accompanied cAMP-induced G1-phase cell-cycle arrest of wild-type S49 cells. Within 2 h, 8-CPT-cAMP altered expression of 152 genes that contain evolutionarily conserved cAMP-response elements within 5 kb of transcriptional start sites, including the circadian clock gene Per1. Thus, cAMP through its activation of PKA produces extensive transcriptional regulation in eukaryotic cells. These transcriptional networks include a primary group of cAMP-response element-containing genes and secondary networks that include the circadian clock.

BACKGROUND: Early transition to labor remains a major cause of infant mortality, yet the causes are largely unknown. Although several marker genes have been identified, little is known about the underlying global gene expression patterns and pathways that orchestrate these striking changes. RESULTS: We performed a detailed time-course study of over 9,000 genes in mouse myometrium at defined physiological states: non-pregnant, mid-gestation, late gestation, and postpartum. This dataset allowed us to identify distinct patterns of gene expression that correspond to phases of myometrial 'quiescence','term activation', and 'postpartum involution'. Using recently developed functional mapping tools (HOPACH (hierarchical ordered partitioning and collapsing hybrid) and GenMAPP 2.0), we have identified new potential transcriptional regulatory gene networks mediating the transition from quiescence to term activation. CONCLUSIONS: These results implicate the myometrium as an essential regulator of endocrine hormone (cortisol and progesterone synthesis) and signaling pathways (cyclic AMP and cyclic GMP stimulation) that direct quiescence via the transcriptional upregulation of both novel and previously associated regulators. With term activation, we observe the upregulation of cytoskeletal remodeling mediators (intermediate filaments), cell junctions, transcriptional regulators, and the coordinate downregulation of negative control checkpoints of smooth muscle contractile signaling. This analysis provides new evidence of multiple parallel mechanisms of uterine contractile regulation and presents new putative targets for regulating myometrial transformation and contraction.

BACKGROUND: Developmental, physiological and tissue engineering studies critical to the development of successful myocardial regeneration therapies require new ways to effectively visualize and isolate large numbers of fluorescently labeled, functional cardiomyocytes. METHODOLOGY/PRINCIPAL FINDINGS: Here we describe methods for the clonal expansion of engineered hESCs and make available a suite of lentiviral vectors for that combine Blasticidin, Neomycin and Puromycin resistance based drug selection of pure populations of stem cells and cardiomyocytes with ubiquitous or lineage-specific promoters that direct expression of fluorescent proteins to visualize and track cardiomyocytes and their progenitors. The phospho-glycerate kinase (PGK) promoter was used to ubiquitously direct expression of histone-2B fused eGFP and mCherry proteins to the nucleus to monitor DNA content and enable tracking of cell migration and lineage. Vectors with T/Brachyury and alpha-myosin heavy chain (alphaMHC) promoters targeted fluorescent or drug-resistance proteins to early mesoderm and cardiomyocytes. The drug selection protocol yielded 96% pure cardiomyocytes that could be cultured for over 4 months. Puromycin-selected cardiomyocytes exhibited a gene expression profile similar to that of adult human cardiomyocytes and generated force and action potentials consistent with normal fetal cardiomyocytes, documenting these parameters in hESC-derived cardiomyocytes and validating that the selected cells retained normal differentiation and function. CONCLUSION/SIGNIFICANCE: The protocols, vectors and gene expression data comprise tools to enhance cardiomyocyte production for large-scale applications.

The second messenger cAMP acts via protein kinase A (PKA) to induce apoptosis by mechanisms that are poorly understood. Here, we assessed a role for mitochondria and analyzed gene expression in cAMP/PKA-promoted apoptosis by comparing wild-type (WT) S49 lymphoma cells and the S49 variant, D-(cAMP-deathless), which lacks cAMP-promoted apoptosis but has wild-type levels of PKA activity and cAMP-promoted G1 growth arrest. Treatment of WT, but not D-, S49 cells with 8-CPT-cAMP for 24 h induced loss of mitochondrial membrane potential, mitochondrial release of cytochrome c and Smac and increase in caspase-3 activity. Gene expression analysis (using Affymetrix 430 2.0 Arrays) revealed that WT and D- cells incubated with 8-CPT-cAMP have similar, but non-identical, extents of cAMP-regulated gene expression at 2h (~800 transcripts) and 6h (~1000 transcripts)(|Fold|>2, P<0.06); by contrast, at 24h ~2500 and ~1100 transcripts were changed in WT and D- cells, respectively. Using an approach that combined regression analysis, clustering and functional annotation to identify transcripts that showed differential expression between WT and D- cells, we found differences in cAMP-mediated regulation of mRNAs involved in transcriptional repression, apoptosis, the cell cycle, RNA splicing, Golgi and lysosomes. The 2 cell lines differed in CREB phosphorylation and expression of the transcriptional inhibitor Icer and in cAMP-regulated expression of genes in the Inhibitor of apoptosis (IAP) and Bcl families. The findings indicate that cAMP/PKA-promoted apoptosis of lymphoid cells occurs via mitochondrial-mediated events and imply that such apoptosis involves gene networks in multiple biochemical pathways.

Cytoscape is a free software package for visualizing, modeling and analyzing molecular and genetic interaction networks. This protocol explains how to use Cytoscape to analyze the results of mRNA expression profiling, and other functional genomics and proteomics experiments, in the context of an interaction network obtained for genes of interest. Five major steps are described:(i) obtaining a gene or protein network,(ii) displaying the network using layout algorithms,(iii) integrating with gene expression and other functional attributes,(iv) identifying putative complexes and functional modules and (v) identifying enriched Gene Ontology annotations in the network. These steps provide a broad sample of the types of analyses performed by Cytoscape.

BACKGROUND: High-throughput mutagenesis of the mammalian genome is a powerful means to facilitate analysis of gene function. Gene trapping in embryonic stem cells (ESCs) is the most widely used form of insertional mutagenesis in mammals. However, the rules governing its efficiency are not fully understood, and the effects of vector design on the likelihood of gene-trapping events have not been tested on a genome-wide scale. METHODOLOGY/PRINCIPAL FINDINGS: In this study, we used public gene-trap data to model gene-trap likelihood. Using the association of gene length and gene expression with gene-trap likelihood, we constructed spline-based regression models that characterize which genes are susceptible and which genes are resistant to gene-trapping techniques. We report results for three classes of gene-trap vectors, showing that both length and expression are significant determinants of trap likelihood for all vectors. Using our models, we also quantitatively identified hotspots of gene-trap activity, which represent loci where the high likelihood of vector insertion is controlled by factors other than length and expression. These formalized statistical models describe a high proportion of the variance in the likelihood of a gene being trapped by expression-dependent vectors and a lower, but still significant, proportion of the variance for vectors that are predicted to be independent of endogenous gene expression. CONCLUSIONS/SIGNIFICANCE: The findings of significant expression and length effects reported here further the understanding of the determinants of vector insertion. Results from this analysis can be applied to help identify other important determinants of this important biological phenomenon and could assist planning of large-scale mutagenesis efforts.

BACKGROUND: Skeletal muscle remodeling is a critical component of an organism's response to environmental changes. Exercise causes structural changes in muscle and can induce phase shifts in circadian rhythms, fluctuations in physiology and behavior with a period of around 24 hours that are maintained by a core clock mechanism. Both exercise-induced remodeling and circadian rhythms rely on the transcriptional regulation of key genes. RESULTS: We used DNA microarrays to determine the effects of resistance exercise (RE) on gene regulation in biopsy samples of human quadriceps muscle obtained 6 and 18 hours after an acute bout of isotonic exercise with one leg. We also profiled diurnal gene regulation at the same time points (2000 and 0800 hours) in the non-exercised leg. Comparison of our results with published circadian gene profiles in mice identified 44 putative genes that were regulated in a circadian fashion. We then used quantitative PCR to validate the circadian expression of selected gene orthologs in mouse skeletal muscle. CONCLUSIONS: The coordinated regulation of the circadian clock genes Cry1, Per2, and Bmal1 6 hours after RE and diurnal genes 18 hours after RE in the exercised leg suggest that RE may directly modulate circadian rhythms in human skeletal muscle.

p53 is critical in regulating the differentiation of ES and induced pluripotent stem (iPS) cells. Here, we report a whole-genome study of p53-mediated DNA damage signaling in mouse ES cells. Systems analyses reveal that binding of p53 at the promoter region significantly correlates with gene activation but not with repression. Unexpectedly, we identify a regulatory mode for p53-mediated repression through interfering with distal enhancer activity. Importantly, many ES cell-enriched core transcription factors are p53-repressed genes. Further analyses demonstrate that p53-repressed genes are functionally associated with ES/iPS cell status while p53-activated genes are linked to differentiation. p53-activated genes and -repressed genes also display distinguishable features of expression levels and epigenetic markers. Upon DNA damage, p53 regulates the self-renewal and pluripotency of ES cells. Together, these results support a model where, in response to DNA damage, p53 affects the status of ES cells through activating differentiation-associated genes and repressing ES cell-enriched genes.

Vascular endothelial cells derived from human pluripotent stem cells have substantial potential for the development of novel vascular therapeutics and cell based therapies for the repair of ischemic damage. To gain maximum benefit from this source of cells a complete understanding of the changes in gene expression and how they are regulated is required. microRNAs have been demonstrated to play a critical roles in controlling stem cell pluripotency and differentiation and are important for mature endothelial cell function. Specific microRNAs which determine stem cell fate have been identified for a number of different cell lineages, however in the case of differentiation and specification of vascular endothelial cells this is yet to be fully elucidated.

The process of differentiation of embryonic stem cells (ESCs) is currently becoming the focus of many systems biologists not only due to mechanistic interest but also since it is expected to play an increasingly important role in regenerative medicine, in particular with the advert to induced pluripotent stem cells. These ESCs give rise to the formation of the three germ layers and therefore to the formation of all tissues and organs. Here, we present a computational method for inferring regulatory interactions between the genes involved in ESC differentiation based on time resolved microarray profiles. Fully quantitative methods are commonly unavailable on such large-scale data; on the other hand, purely qualitative methods may fail to capture some of the more detailed regulations. Our method combines the beneficial aspects of qualitative and quantitative (ODE-based) modeling approaches searching for quantitative interaction coefficients in a discrete and qualitative state space. We further optimize on an ensemble of networks to detect essential properties and compare networks with respect to robustness. Applied to a toy model our method is able to reconstruct the original network and outperforms an entire discrete boolean approach. In particular, we show that including prior knowledge leads to more accurate results. Applied to data from differentiating mouse ESCs reveals new regulatory interactions, in particular we confirm the activation of Foxh1 through Oct4, mediating Nodal signaling.

New evidence has revealed interesting aspects of how the Wnt-β-catenin pathway controls self-renewal and lineage differentiation of pluripotent embryonic stem cells. Although Wnt-β-catenin signaling is dispensable for the self-renewal of naive mouse embryonic stem cells, it facilitates their expansion and resistance to differentiation through an unconventional dual mechanism involving the transcriptional repressor T cell factor (TCF) 3 and the transcriptional activator TCF1.

Embryonic stem (ES) cells are distinguished by their ability to undergo unlimited self-renewal while retaining pluripotency, the capacity to specify cells of all germ layers. Alternative splicing contributes to these biological processes by vastly increasing the protein coding repertoire, enabling genes to code for novel variants that may confer different biological functions. The homeodomain transcription factor Nanog acts collaboratively with the core factors Oct4 and Sox2 to govern the maintenance of pluripotency. We have discovered that Nanog is regulated by alternative splicing. Two novel exons and six sub-exons have been identified that extend the known Nanog gene structure and protein coding capacity. Alternative splicing results in two novel Nanog protein variants with altered capacities for self-renewal and pluripotency in ES cells. Our previous results have implicated the C-terminal domain, including the tryptophan-rich (WR) domain of Nanog, to be important for the function of Nanog. Using point mutation analyses, the serine 2 (Ser-2) residue of Nanog has been identified as critical for ES cell self-renewal and for stabilizing a pluripotent gene signature. An inducible conditional knockout was created to test the ability of new Nanog variants to genetically complement Nanog null ES cells. These results reveal for the first time an expanded Nanog protein coding capacity. We further reveal that a short region of the N-terminal domain and a single phosphorylatable Ser-2 is essential for the maintenance of self-renewal and pluripotency, demonstrating that this region of the protein is highly regulated.

In embryonic stem cells (ESCs), Wnt-responsive development-related genes are silenced to maintain pluripotency and their expression is activated during differentiation. Acetylation of histones by histone acetyltransferases stimulates transcription, whereas deacetylation of histones by HDACs is correlated with transcriptional repression. Although Wnt-mediated gene transcription has been intimately linked to the acetylation or deacetylation of histones, how Wnt signaling regulates this type of histone modification is poorly understood. Here, we report that Smek, a regulatory subunit of protein phosphatase 4 (PP4) complex, plays an important role in histone deacetylation and silencing of the Wnt-responsive gene, brachyury, in ESCs. Smek mediates recruitment of PP4c and HDAC1 to the Tcf/Lef binding site of the brachyury promoter and inhibits brachyury expression in ESCs. Activation of Wnt signaling during differentiation causes disruption of the Smek/PP4c/HDAC1 complex, resulting in an increase in histones H3 and H4 acetylation at the brachyury gene locus. These results suggest that the Smek-containing PP4 complex represses transcription of Wnt-responsive development-related genes through histone deacetylation, and that this complex is essential for ESC pluripotency maintenance.

Mouse Embryonic Stem (ES) cells express a unique set of microRNAs (miRNAs), the miR-290-295 cluster. To elucidate the role of these miRNAs and how they integrate into the ES cell regulatory network requires identification of their direct regulatory targets. The difficulty, however, arises from the limited complementarity of metazoan miRNAs to their targets, with the interaction requiring as few as six nucleotides of the miRNA seed sequence. To identify miR-294 targets, we used Dicer1-null ES cells, which lack all endogenous mature miRNAs, and introduced just miR-294 into these ES cells. We then employed two approaches to discover miR-294 targets in mouse ES cells: transcriptome profiling using microarrays and a biochemical approach to isolate mRNA targets associated with the Argonaute2 (Ago2) protein of the RISC (RNA Induced Silencing Complex) effector, followed by RNA-sequencing. In the absence of Dicer1, the RISC complexes are largely devoid of mature miRNAs and should therefore contain only transfected miR-294 and its base-paired targets. Our data suggest that miR-294 may promote pluripotency by regulating a subset of c-Myc target genes and upregulating pluripotency-associated genes such as Lin28.

The mmu-miR-290-295 cluster codes for a family of microRNAs (miRNAs) that are expressed de novo during early embryogenesis and are specific for mouse embryonic stem cells (ESC) and embryonic carcinoma cells (ECC). Detailed sequence analysis and alignment studies of miR-290-295 precursors demonstrated that the cluster has evolved by repeated duplication events of the ancient miR-290 precursor. We show that under serum starvation, overexpression of miR-290-295 miRNAs withhold ES cells from early differentiation, ensures their high proliferation rate and capacity for forming alkaline phosphate positive colonies. Transcriptome analysis revealed that differentiation related marker genes are underexpressed upon high miR-290-295 level. Importantly, miR-290-295 overexpression prevents ES cells from accumulation in G1 phase at low serum level, and seems to regulate cell cycle in different phases. Our data underline that miR-290-295 miRNAs contribute to the natural absence of G1 checkpoint in embryonic stem cells. We define the cell cycle regulators Wee1 and Fbxl5 as potential direct targets of miR-290-295 miRNAs in vitro. Our results suggest that miR-290-295 miRNAs exhibit their effect predominantly through the regulation of cell cycle phase distribution.

Murine embryonic stem (ES) cells are defined by continuous self-renewal and pluripotency. A diverse repertoire of protein isoforms arising from alternative splicing is expressed in ES cells without defined biological roles. Sall4, a transcription factor essential for pluripotency, exists as two isoforms (Sall4a and Sall4b). Both isoforms can form homodimers and a heterodimer with each other, and each can interact with Nanog. By genomewide location analysis, we determined that Sall4a and Sall4b have overlapping, but not identical binding sites within the ES cell genome. In addition, Sall4b, but not Sall4a, binds preferentially to highly expressed loci in ES cells. Sall4a and Sall4b binding sites are distinguished by both epigenetic marks at target loci and their clustering with binding sites of other pluripotency factors. When ES cells expressing a single isoform of Sall4 are generated, Sall4b alone could maintain the pluripotent state, although it could not completely suppress all differentiation markers. Sall4a and Sall4b collaborate in maintenance of the pluripotent state but play distinct roles. Our work is novel in establishing such isoform-specific differences in ES cells.

Members of the transforming growth factor-beta superfamily play essential roles in both the pluripotency and differentiation of embryonic stem (ES) cells. Although bone morphogenic proteins (BMPs) maintain pluripotency of undifferentiated mouse ES cells, the role of autocrine Nodal signaling is less clear. Pharmacological, molecular, and genetic methods were used to further understand the roles and potential interactions of these pathways. Treatment of undifferentiated ES cells with SB431542, a pharmacological inhibitor of Smad2 signaling, resulted in a rapid reduction of phosphorylated Smad2 and altered the expression of several putative downstream targets. Unexpectedly, inhibition of the Nodal signaling pathway resulted in enhanced BMP signaling, as assessed by Smad1/5 phosphorylation. SB431542-treated cells also demonstrated significant induction of the Id genes, which are known direct targets of BMP signaling and important factors in ES cell pluripotency. Inhibition of BMP signaling decreased the SB431542-mediated phosphorylation of Smad1/5 and induction of Id genes, suggesting that BMP signaling is necessary for some Smad2-mediated activity. Because Smad7, a known inhibitory factor to both Nodal and BMP signaling, was down-regulated following inhibition of Nodal-Smad2 signaling, the contribution of Smad7 to the cross-talk between the transforming growth factor-beta pathways in ES cells was examined. Biochemical manipulation of Smad7 expression, through shRNA knockdown or inducible gene expression, significantly reduced the SB431542-mediated phosphorylation of Smad1/5 and induction of the Id genes. We conclude that autocrine Nodal signaling in undifferentiated mouse ES cells modulates the vital pluripotency pathway of BMP signaling.