Published papers

Pluripotent human embryonic stem cells (hESCs) can differentiate into multiple cell lineages, thus, providing one of the best platforms to study molecular mechanisms during cell differentiation. Recently, we have reported rapid and efficient differentiation of hESCs into functional neurons by introducing a cocktail of synthetic mRNAs encoding five transcription factors (TFs): NEUROG1, NEUROG2, NEUROG3, NEUROD1, and NEUROD2. Here we further tested a possibility that even single transcription factors, when expressed ectopically, can differentiate hESCs into neurons. To this end, we established hESC lines in which each of these TFs can be overexpressed by the doxycycline-inducible piggyBac vector. The overexpression of any of these five TFs indeed caused a rapid and rather uniform differentiation of hESCs, which were identified as neurons based on their morphologies, qRT-PCR, and immunohistochemistry. Furthermore, calcium-imaging analyses and patch clamp recordings demonstrated that these differentiated cells are electrophysiologically functional. Interestingly, neural differentiations occurred despite the cell culture conditions that rather promote the maintenance of the undifferentiated state. These results indicate that over-expression of each of these five TFs can override the pluripotency-specific gene network and force hESCs to differentiate into neurons.

Efficient differentiation of human pluripotent stem cells (hPSCs) into neurons is paramount for disease modeling, drug screening, and cell transplantation therapy in regenerative medicine. In this manuscript, we report the capability of five transcription factors (TFs) toward this aim: NEUROG1, NEUROG2, NEUROG3, NEUROD1, and NEUROD2. In contrast to previous methods that have shortcomings in their speed and efficiency, a cocktail of these TFs as synthetic mRNAs can differentiate hPSCs into neurons in 7 days, judged by calcium imaging and electrophysiology. They exhibit motor neuron phenotypes based on immunostaining. These results indicate the establishment of a novel method for rapid, efficient, and footprint-free differentiation of functional neurons from hPSCs.

Zygotic genome activation (ZGA, a.k.a. zygotic gene activation) is a critical event in development, when the paternally derived genome and maternally derived genome begin to be activated and transcribed after fertilization. Major ZGA occurs at the two-cell stage in mice and the four- to eight-cell stage in human preimplantation embryos. It has been thought that ZGA exists to provide RNAs and proteins supporting embryonic development after supplies stored in oocytes are used up; however, this paradigm does not seem to explain recent findings. For example, many ZGA genes-once activated-are quickly turned off, and thus ZGA forms a transient wave of transcriptional activation. In addition, ZGA genes are not evolutionarily conserved. In this review, we address these issues by focusing on Zscan4 (zinc finger and SCAN domain-containing 4), which was identified for its specific expression in preimplantation embryos during ZGA. Detailed molecular analyses of Zscan4 expression and function have revealed common features of Zscan4-associated events (Z4 events) in mouse embryonic stem cells and ZGA in preimplantation embryos. One feature is a rapid derepression and rerepression of constitutive heterochromatin, which includes pericentromeric major satellites and telomeres, and facultative heterochromatin, which includes retrotransposons and Z4 event-associated genes. We propose that the Z4 event superimposed on ZGA plays a critical role in the maintenance of genome and chromosome integrity in preimplantation embryos by promoting correction of DNA damage and chromosome abnormalities.

Harnessing epigenetic regulation is crucial for the efficient and proper differentiation of pluripotent stem cells (PSCs) into desired cell types. Histone H3 lysine 27 trimethylation (H3K27me3) functions as a barrier against cell differentiation through the suppression of developmental gene expression in PSCs. Here, we have generated human PSC (hPSC) lines in which genome-wide reduction of H3K27me3 can be induced by ectopic expression of the catalytic domain of the histone demethylase JMJD3 (called JMJD3c). We found that transient, forced demethylation of H3K27me3 alone triggers the upregulation of mesoendodermal genes, even when the culture conditions for the hPSCs are not changed. Furthermore, transient and forced expression of JMJD3c followed by the forced expression of lineage-defining transcription factors enabled the hPSCs to activate tissue-specific genes directly. We have also shown that the introduction of JMJD3c facilitates the differentiation of hPSCs into functional hepatic cells and skeletal muscle cells. These results suggest the utility of the direct manipulation of epigenomes for generating desired cell types from hPSCs for cell transplantation therapy and platforms for drug screenings.

Mouse zinc finger and SCAN domain containing 4 (Zscan4) proteins, which are encoded by multiple copies of Zscan4 genes, are expressed specifically in preimplantation embryos in vivo and embryonic stem (ES) cells in vitro. However, the expression patterns of mouse Zscan4 in vivo have been largely elusive. Here, we show that Zscan4 proteins are expressed in adult ovaries and testes. In ovaries, Zscan4 proteins were detected in germinal vesicle (GV) stage oocytes in antral follicles, indicating that Zscan4 genes are activated during the diplotene/dictyate stage in meiotic prophase I. Remarkably, Zscan4 showed different spatial localization patterns between two distinct GV oocytes, which can be distinguished by global chromatin organization—surrounded nucleolus (SN) and non-surrounded nucleolus (NSN). These spatiotemporal differences in Zscan4 localizations correlated with the transition of RNA polymerase II-mediated transcriptional status during GV oocyte maturation. In testes, Zscan4 proteins were detected in spermatocytes at late pachytene/diplotene stages and in Sertoli cells. These results suggest that Zscan4 may play critical roles during late meiotic prophase in both males and females.

Mouse Zinc finger and SCAN domain containing 4 (Zscan4) is encoded in multiple copies of Zscan4 genes, which are expressed in late two-cell stage preimplantation embryos and in 1–5% of the embryonic stem (ES) cell population at a given time. Due to the highly identical nucleotide sequences of multiple copies of Zscan4 paralogs and pseudogenes in the mouse Zscan4 genomic cluster, previous analyses have been done using exogenous transgenes under the regulation of Zscan4c promoter. In this manuscript, we generated knock-in mouse ES cell lines and mouse lines, in which the expression of endogenous Zscan4c, one of the Zscan4 genes, can be specifically monitored with a green fluorescent protein variant, Emerald. Interestingly, we found that only ∼30% of Zscan4-immunopositive ES cells were Emerald positive, suggesting that even when the Zscan4 locus is active, not all Zscan4 genes are expressed synchronously. We also carried out mass spectrometry of protein complexes associated with endogenous Zscan4 proteins. Taken together, our genetic engineering at an endogenous Zscan4c gene provides the first clue for the expression and function of each gene copy of Zscan4 locus in a physiological context.

Summary: Specific neuronal types derived from embryonic stem cells (ESCs) can facilitate mechanistic studies and potentially aid in regenerative medicine. Existing induction methods, however, mostly rely on the effects of the combined action of multiple added growth factors, which generally tend to result in mixed populations of neurons. Here, we report that overexpression of specific transcription factors (TFs) in ESCs can rather guide the differentiation of ESCs towards specific neuron lineages. Analysis of data on gene expression changes 2 d after induction of each of 185 TFs implicated candidate TFs for further ESC differentiation studies. Induction of 23 TFs (out of 49 TFs tested) for 6 d facilitated neural differentiation of ESCs as inferred from increased proportion of cells with neural progenitor marker PSA-NCAM. We identified early activation of the Notch signaling pathway as a common feature of most potent inducers of neural differentiation. The majority of neuron-like cells generated by induction of Ascl1, Smad7, Nr2f1, Dlx2, Dlx4, Nr2f2, Barhl2, and Lhx1 were GABA-positive and expressed other markers of GABAergic neurons. In the same way, we identified Lmx1a and Nr4a2 as inducers for neurons bearing dopaminergic markers and Isl1, Fezf2, and St18 for cholinergic motor neurons. A time-course experiment with induction of Ascl1 showed early upregulation of most neural-specific messenger RNA (mRNA) and microRNAs (miRNAs). Sets of Ascl1-induced mRNAs and miRNAs were enriched in Ascl1 targets. In further studies, enrichment of cells obtained with the induction of Ascl1, Smad7, and Nr2f1 using microbeads resulted in essentially pure population of neuron-like cells with expression profiles similar to neural tissues and expressed markers of GABAergic neurons. In summary, this study indicates that induction of transcription factors is a promising approach to generate cultures that show the transcription profiles characteristic of specific neural cell types.

Mouse embryonic stem cells (ESCs) can differentiate into a wide range-and possibly all cell types in vitro, and thus provide an ideal platform to study systematically the action of transcription factors (TFs) in cell differentiation. Previously, we have generated and analyzed 137 TF-inducible mouse ESC lines. As an extension of this "NIA Mouse ESC Bank" we generated and characterized 48 additional mouse ESC lines, in which single TFs in each line could be induced in a doxycycline-controllable manner. Together, with the previous ESC lines, the bank now comprises 185 TF-manipulable ESC lines (>10% of all mouse TFs). Global gene expression (transcriptome) profiling revealed that the induction of individual TFs in mouse ESCs for 48 hours shifts their transcriptomes toward specific differentiation fates (e.g., neural lineages by Myt1 Isl1, and St18; mesodermal lineages by Pitx1, Pitx2, Barhl2, and Lmx1a; white blood cells by Myb, Etv2, and Tbx6, and ovary by Pitx1, Pitx2, and Dmrtc2). These data also provide and lists of inferred target genes of each TF and possible functions of these TFs. The results demonstrate the utility of mouse ESC lines and their transcriptome data for understanding the mechanism of cell differentiation and the function of TFs.

Retinoic acid (RA) is one of the most potent inducers of differentiation of mouse embryonic stem cells (ESCs). However, previous studies show that RA treatment of cells cultured in the presence of a leukemia inhibitory factor (LIF) also result in the upregulation of a gene called Zscan4, whose transient expression is a marker for undifferentiated ESCs. We explored the balance between these two seemingly antagonistic effects of RA. ESCs indeed differentiated in the presence of LIF after RA treatment, but colonies of undifferentiated ESCs eventually emerged from these differentiated cells — even in the presence of RA. These colonies, named secondary colonies, consist of three cell types: typical undifferentiated ESCs expressing pluripotency genes such as Pou5f1, Sox2, and Nanog; cells expressing Zscan4; and endodermal-like cells located at the periphery of the colony. The capacity to form secondary colonies was confirmed for all eight tested ESC lines. Cells from the secondary colonies — after transfer to the standard ESC medium — retained pluripotency, judged by their strong alkaline phosphatase (ALP) staining, typical colony morphology, gene expression profile, stable karyotype, capacity to differentiate into all three germ layers in embryoid body formation assays, and successful contribution to chimeras after injection into blastocysts. Based on flow cytometry analysis (FACS), the proportion of Zscan4-positive cells in secondary colonies was higher than in standard ESC colonies, which may explain the capacity of ESCs to resist the differentiating effects of RA and instead form secondary colonies of undifferentiated ESCs. This hypothesis is supported by cell-lineage tracing analysis, which showed that most cells in the secondary colonies were descendents of cells transiently expressing Zscan4.

Neural crest cells (NC cells) are multipotent cells that emerge from the edge of the neural folds and migrate throughout the developing embryo. Although the gene regulatory network for generation of NC cells has been elucidated in detail, it has not been revealed which of the factors in the network are pivotal to directing NC identity. In this study we analyzed the gene expression profile of a pure NC subpopulation isolated from Sox10-IRES-Venus mice and investigated whether these genes played a key role in the direct conversion of Sox10-IRES-Venus mouse embryonic fibroblasts (MEFs) into NC cells. The comparative molecular profiles of NC cells and neural tube cells in 9.5-day embryos revealed genes including transcription factors selectively expressed in developing trunk NC cells. Among 25 NC cell-specific transcription factor genes tested, SOX10 and SOX9 were capable of converting MEFs into SOX10-positive (SOX10+) cells. The SOX10+ cells were then shown to differentiate into neurons, glial cells, smooth muscle cells, adipocytes and osteoblasts. These SOX10+ cells also showed limited self-renewal ability, suggesting that SOX10 and SOX9 directly converted MEFs into NC cells. Conversely, the remaining transcription factors, including well-known NC cell specifiers, were unable to convert MEFs into SOX10+ NC cells. These results suggest that SOX10 and SOX9 are the key factors necessary for the direct conversion of MEFs into NC cells.

We have developed ExAtlas, an on-line software tool for meta-analysis and visualization of gene expression data. In contrast to existing software tools, ExAtlas compares multi-component data sets and generates results for all combinations (e.g. all gene expression profiles versus all Gene Ontology annotations). ExAtlas handles both users' own data and data extracted semi-automatically from the public repository (GEO/NCBI database). ExAtlas provides a variety of tools for meta-analyses: (1) standard meta-analysis (fixed effects, random effects, z-score, and Fisher's methods); (2) analyses of global correlations between gene expression data sets; (3) gene set enrichment; (4) gene set overlap; (5) gene association by expression profile; (6) gene specificity; and (7) statistical analysis (ANOVA, pairwise comparison, and PCA). ExAtlas produces graphical outputs, including heatmaps, scatter-plots, bar-charts, and three-dimensional images. Some of the most widely used public data sets (e.g. GNF/BioGPS, Gene Ontology, KEGG, GAD phenotypes, BrainScan, ENCODE ChIP-seq, and protein-protein interaction) are pre-loaded and can be used for functional annotations.

Human pluripotent stem cells (hPSCs) have the capacity to differentiate into essentially all cell types in the body. Such differentiation can be directed to specific cell types by appropriate cell culture conditions or overexpressing lineage-defining transcription factors (TFs). Especially, for the activation of myogenic program, early studies have shown the effectiveness of enforced expression of TFs associated with myogenic differentiation, such as PAX7 and MYOD1. However, the efficiency of direct differentiation was rather low, most likely due to chromatin features unique to hPSCs, which hinder the access of TFs to genes involved in muscle differentiation. Indeed, recent studies have demonstrated that ectopic expression of epigenetic-modifying factors such as a histone demethylase and an ATP-dependent remodeling factor significantly enhances myogenic differentiation from hPSCs. In this article, we review the recent progress for in vitro generation of skeletal muscles from hPSCs through forced epigenetic and transcriptional manipulation.

Motivation: The analysis of RNA-Seq data from individual differentiating cells enables us to reconstruct the differentiation process and the degree of differentiation (in pseudo-time) of each cell. Such analyses can reveal detailed expression dynamics and functional relationships for differentiation. To further elucidate differentiation processes, more insight into gene regulatory networks is required. The pseudo-time can be regarded as time information and, therefore, single-cell RNASeq data are time-course data with high time resolution. Although time-course data are useful for inferring networks, conventional inference algorithms for such data suffer from high time complexity when the number of samples and genes is large. Therefore, a novel algorithm is necessary to infer networks from single-cell RNA-Seq during differentiation. Results: In this study, we developed the novel and efficient algorithm SCODE to infer regulatory networks, based on ordinary differential equations. We applied SCODE to three single-cell RNASeq datasets and confirmed that SCODE can reconstruct observed expression dynamics. We evaluated SCODE by comparing its inferred networks with use of a DNaseI-footprint based network. The performance of SCODE was best for two of the datasets and nearly best for the remaining dataset. We also compared the runtimes and showed that the runtimes for SCODE are significantly shorter than for alternatives. Thus, our algorithm provides a promising approach for further singlecell differentiation analyses.

We previously reported that Zscan4 showed heterogeneous expression patterns in mouse embryonic stem (ES) cells. To identify genes that show similar expression patterns, we carried out high-throughput in situ hybridization assays on ES cell cultures for 244 genes. Most of the genes are involved in transcriptional regulation, and were selected using microarray-based comparisons of gene expression profiles in ES and embryonal carcinoma (EC) cells versus differentiated cell types. Pou5f1 (Oct4, Oct3/4) and Krt8 (EndoA) were used as controls. Hybridization signals were detected on ES cell colonies for 147 genes (60%). The majority (136 genes) of them showed relatively homogeneous expression in ES cell colonies. However, we found that two genes unequivocally showed Zscan4-like spotted expression pattern (spot-in-colony pattern; Whsc2 and Rhox9). We also found that nine genes showed relatively heterogeneous expression pattern (mosaic-in-colony pattern: Zfp42/Rex1, Rest, Atf4, Pa2g4, E2f2, Nanog, Dppa3/Pgc7/Stella, Esrrb, and Fscn1). Among these genes, Zfp42/Rex1 showed unequivocally heterogeneous expression in individual ES cells prepared by the CytoSpin. These results show the presence of different types or states of cells within ES cell cultures otherwise thought to be undifferentiated and homogeneous, suggesting a previously unappreciated complexity in ES cell cultures.

Background: Recently, several populations of postnatal stem cells, such as multipotent adult progenitor cells (MAPCs), have been described that have broader differentiation ability than classical adult stem cells. Here we compare the transcriptome of pluripotent embryonic stem cells (ESCs), MAPCs, and lineage-restricted mesenchymal stem cells (MSCs) to determine their relationship. Results: Applying principal component analysis, non-negative matrix factorization and k-means clustering algorithms to the gene-expression data, we identified a unique gene-expression profile for MAPCs. Apart from the ESC-specific transcription factor Oct4 and other ESC transcripts, some of them associated with maintaining ESC pluripotency, MAPCs also express transcripts characteristic of early endoderm and mesoderm. MAPCs do not, however, express Nanog or Sox2, two other key transcription factors involved in maintaining ESC properties. This unique molecular signature was seen irrespective of the microarray platform used and was very similar for both mouse and rat MAPCs. As MSC-like cells isolated under MAPC conditions are virtually identical to MSCs, and MSCs cultured in MAPC conditions do not upregulate MAPC-expressed transcripts, the MAPC signature is cell-type specific and not merely the result of differing culture conditions. Conclusion: Multivariate analysis techniques clustered stem cells on the basis of their expressed gene profile, and the genes determining this clustering reflected the stem cells' differentiation potential in vitro. This comparative transcriptome analysis should significantly aid the isolation and culture of MAPCs and MAPC-like cells, and form the basis for studies to gain insights into genes that confer on these cells their greater developmental potency. © 2007 Ulloa-Montoya et al.; licensee BioMed Central Ltd.

Loss of imprinting (LOI) of the insulin-like growth factor-II gene (IGF2), leading to abnormal activation of the normally silent maternal allele, is a common human epigenetic population variant associated with a 5-fold increased frequency of colorectal neoplasia. Here, we show first that LOI leads specifically to increased expression of proliferation-related genes in mouse intestinal crypts. Surprisingly, LOI(+) mice also have enhanced sensitivity to IGF-II signaling, not simply increased IGF-II levels, because in vivo blockade with NVP-AEW541, a specific inhibitor of the IGF-II signaling receptor, showed reduction of proliferation-related gene expression to levels half that seen in LOI(-) mice. Signal transduction assays in microfluidic chips confirmed this enhanced sensitivity with marked augmentation of Akt/PKB signaling in LOI(+) cells at low doses of IGF-II, which was reduced in the presence of the inhibitor to levels below those found in LOI(-) cells, and was associated with increased expression of the IGF1 and insulin receptor genes. We exploited this increased IGF-II sensitivity to develop an in vivo chemopreventive strategy using the azoxymethane (AOM) mutagenesis model. LOI(+) mice treated with AOM showed a 60% increase in premalignant aberrant crypt foci (ACF) formation over LOI(-) mice. In vivo IGF-II blockade with NVP-AEW541 abrogated this effect, reducing ACF to a level 30% lower even than found in exposed LOI(-) mice. Thus, LOI increases cancer risk in a counterintuitive way, by increasing the sensitivity of the IGF-II signaling pathway itself, providing a previously undescribed epigenetic chemoprevention strategy in which cells with LOI are "IGF-II addicted" and undergo reduced tumorigenesis in the colon upon IGF-II pathway blockade. © 2007 by The National Academy of Sciences of the USA.

We present the AGEMAP (Atlas of Gene Expression in Mouse Aging Project) gene expression database, which is a resource that catalogs changes in gene expression as a function of age in mice. The AGEMAP database includes expression changes for 8,932 genes in 16 tissues as a function of age. We found great heterogeneity in the amount of transcriptional changes with age in different tissues. Some tissues displayed large transcriptional differences in old mice, suggesting that these tissues may contribute strongly to organismal decline. Other tissues showed few or no changes in expression with age, indicating strong levels of homeostasis throughout life. Based on the pattern of agerelated transcriptional changes, we found that tissues could be classified into one of three aging processes: (1) a pattern common to neural tissues, (2) a pattern for vascular tissues, and (3) a pattern for steroid-responsive tissues. We observed that different tissues age in a coordinated fashion in individual mice, such that certain mice exhibit rapid aging, whereas others exhibit slow aging for multiple tissues. Finally, we compared the transcriptional profiles for aging in mice to those from humans, flies, and worms. We found that genes involved in the electron transport chain show common age regulation in all four species, indicating that these genes may be exceptionally good markers of aging. However, we saw no overall correlation of age regulation between mice and humans, suggesting that aging processes in mice and humans may be fundamentally different.