Zscan4 regulates telomere elongation and genomic stability in ES cells.

Exceptional genomic stability is one of the hallmarks of mouse embryonic stem (ES) cells. However, the genes contributing to this stability remain obscure. We previously identified Zscan4 as a specific marker for two-cell embryo and ES cells. Here we show that Zscan4 is involved in telomere maintenance and long-term genomic stability in ES cells. Only 5% of ES cells express Zscan4 at a given time, but nearly all ES cells activate Zscan4 at least once during nine passages. The transient Zscan4-positive state is associated with rapid telomere extension by telomere recombination and upregulation of meiosis-specific homologous recombination genes, which encode proteins that are colocalized with ZSCAN4 on telomeres. Furthermore, Zscan4 knockdown shortens telomeres, increases karyotype abnormalities and spontaneous sister chromatid exchange, and slows down cell proliferation until reaching crisis by passage eight. Together, our data show a unique mode of genome maintenance in ES cells.

Essential role of chromatin remodeling protein Bptf in early mouse embryos and embryonic stem cells.

We have characterized the biological functions of the chromatin remodeling protein Bptf (Bromodomain PHD-finger Transcription Factor), the largest subunit of NURF (Nucleosome Remodeling Factor) in a mammal. Bptf mutants manifest growth defects at the post-implantation stage and are reabsorbed by E8.5. Histological analyses of lineage markers show that Bptf(-/-) embryos implant but fail to establish a functional distal visceral endoderm. Microarray analysis at early stages of differentiation has identified Bptf-dependent gene targets including homeobox transcriptions factors and genes essential for the development of ectoderm, mesoderm, and both definitive and visceral endoderm. Differentiation of Bptf(-/-) embryonic stem cell lines into embryoid bodies revealed its requirement for development of mesoderm, endoderm, and ectoderm tissue lineages, and uncovered many genes whose activation or repression are Bptf-dependent. We also provide functional and physical links between the Bptf-containing NURF complex and the Smad transcription factors. These results suggest that Bptf may co-regulate some gene targets of this pathway, which is essential for establishment of the visceral endoderm. We conclude that Bptf likely regulates genes and signaling pathways essential for the development of key tissues of the early mouse embryo.

BAF250B-associated SWI/SNF chromatin-remodeling complex is required to maintain undifferentiated mouse embryonic stem cells.

Whether SWI/SNF chromatin remodeling complexes play roles in embryonic stem (ES) cells remains unknown. Here we show that SWI/SNF complexes are present in mouse ES cells, and their composition is dynamically regulated upon induction of ES cell differentiation. For example, the SWI/SNF purified from undifferentiated ES cells contains a high level of BAF155 and a low level of BAF170 (both of which are homologs of yeast SWI3 protein), whereas that from differentiated cells contains nearly equal amounts of both. Moreover, the levels of BAF250A and BAF250B decrease during the differentiation of ES cells, whereas that of BRM increases. The altered expression of SWI/SNF components hinted that these complexes could play roles in ES cell maintenance or differentiation. We therefore generated ES cells with biallelic inactivation of BAF250B and found that these cells display a reduced proliferation rate and an abnormal cell cycle. Importantly, these cells are deficient in the self-renewal capacity of undifferentiated ES cells and exhibit certain phenotypes of differentiated cells, including reduced expression of several pluripotency-related genes and increased expression of some differentiation-related genes. These data suggest that the BAF250B-associated SWI/SNF is essential for mouse ES cells to maintain their normal proliferation and pluripotency. The work presented here underscores the importance of SWI/SNF chromatin remodeling complexes in pluripotent stem cells.

Identification of Pou5f1, Sox2, and Nanog downstream target genes with statistical confidence by applying a novel algorithm to time course microarray and genome-wide chromatin immunoprecipitation data.

BACKGROUND: Target genes of a transcription factor (TF) Pou5f1 (Oct3/4 or Oct4), which is essential for pluripotency maintenance and self-renewal of embryonic stem (ES) cells, have previously been identified based on their response to Pou5f1 manipulation and occurrence of Chromatin-immunoprecipitation (ChIP)-binding sites in promoters. However, many responding genes with binding sites may not be direct targets because response may be mediated by other genes and ChIP-binding site may not be functional in terms of transcription regulation. RESULTS: To reduce the number of false positives, we propose to separate responding genes into groups according to direction, magnitude, and time of response, and to apply the false discovery rate (FDR) criterion to each group individually. Using this novel algorithm with stringent statistical criteria (FDR < 0.2) to a compendium of published and new microarray data (3, 6, 12, and 24 hr after Pou5f1 suppression) and published ChIP data, we identified 420 tentative target genes (TTGs) for Pou5f1. The majority of TTGs (372) were down-regulated after Pou5f1 suppression, indicating that the Pou5f1 functions as an activator of gene expression when it binds to promoters. Interestingly, many activated genes are potent suppressors of transcription, which include polycomb genes, zinc finger TFs, chromatin remodeling factors, and suppressors of signaling. Similar analysis showed that Sox2 and Nanog also function mostly as transcription activators in cooperation with Pou5f1. CONCLUSION: We have identified the most reliable sets of direct target genes for key pluripotency genes - Pou5f1, Sox2, and Nanog, and found that they predominantly function as activators of downstream gene expression. Thus, most genes related to cell differentiation are suppressed indirectly.

An in situ hybridization-based screen for heterogeneously expressed genes in mouse ES cells.

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.

Transcriptome analysis of mouse stem cells and early embryos.

Understanding and harnessing cellular potency are fundamental in biology and are also critical to the future therapeutic use of stem cells. Transcriptome analysis of these pluripotent cells is a first step towards such goals. Starting with sources that include oocytes, blastocysts, and embryonic and adult stem cells, we obtained 249,200 high-quality EST sequences and clustered them with public sequences to produce an index of approximately 30,000 total mouse genes that includes 977 previously unidentified genes. Analysis of gene expression levels by EST frequency identifies genes that characterize preimplantation embryos, embryonic stem cells, and adult stem cells, thus providing potential markers as well as clues to the functional features of these cells. Principal component analysis identified a set of 88 genes whose average expression levels decrease from oocytes to blastocysts, stem cells, postimplantation embryos, and finally to newborn tissues. This can be a first step towards a possible definition of a molecular scale of cellular potency. The sequences and cDNA clones recovered in this work provide a comprehensive resource for genes functioning in early mouse embryos and stem cells. The nonrestricted community access to the resource can accelerate a wide range of research, particularly in reproductive and regenerative medicine.

Induction of specific neuron types by overexpression of single transcription factors.

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.

Emergence of undifferentiated colonies from mouse embryonic stem cells undergoing differentiation by retinoic acid treatment.

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.

Zscan4 restores the developmental potency of embryonic stem cells.

The developmental potency of mouse embryonic stem (ES) cells, which is the ability to contribute to a whole embryo, is known to deteriorate during long-term cell culture. Previously, we have shown that ES cells oscillate between Zscan4(-) and Zscan4(+) states, and the transient activation of Zscan4 is required for the maintenance of telomeres and genome stability of ES cells. Here we show that increasing the frequency of Zscan4 activation in mouse ES cells restores and maintains their developmental potency in long-term cell culture. Injection of a single ES cell with such increased potency into a tetraploid blastocyst gives rise to an entire embryo with a higher success rate. These results not only provide a means to rejuvenate ES cells by manipulating Zscan4 expression, but also indicate the active roles of Zscan4 in the long-term maintenance of ES cell potency.

Systematic repression of transcription factors reveals limited patterns of gene expression changes in ES cells.

Networks of transcription factors (TFs) are thought to determine and maintain the identity of cells. Here we systematically repressed each of 100 TFs with shRNA and carried out global gene expression profiling in mouse embryonic stem (ES) cells. Unexpectedly, only the repression of a handful of TFs significantly affected transcriptomes, which changed in two directions/trajectories: one trajectory by the repression of either Pou5f1 or Sox2; the other trajectory by the repression of either Esrrb, Sall4, Nanog, or Tcfap4. The data suggest that the trajectories of gene expression change are already preconfigured by the gene regulatory network and roughly correspond to extraembryonic and embryonic fates of cell differentiation, respectively. These data also indicate the robustness of the pluripotency gene network, as the transient repression of most TFs did not alter the transcriptomes.

Generation of mouse ES cell lines engineered for the forced induction of transcription factors.

Here we report the generation and characterization of 84 mouse ES cell lines with doxycycline-controllable transcription factors (TFs) which, together with the previous 53 lines, cover 7-10% of all TFs encoded in the mouse genome. Global gene expression profiles of all 137 lines after the induction of TFs for 48 hrs can associate each TF with the direction of ES cell differentiation, regulatory pathways, and mouse phenotypes. These cell lines and microarray data provide building blocks for a variety of future biomedical research applications as a community resource.

Database for mRNA half-life of 19 977 genes obtained by DNA microarray analysis of pluripotent and differentiating mouse embryonic stem cells.

Degradation of mRNA is one of the key processes that control the steady-state level of gene expression. However, the rate of mRNA decay for the majority of genes is not known. We successfully obtained the rate of mRNA decay for 19 977 non-redundant genes by microarray analysis of RNA samples obtained from mouse embryonic stem (ES) cells. Median estimated half-life was 7.1 h and only <100 genes, including Prdm1, Myc, Gadd45 g, Foxa2, Hes5 and Trib1, showed half-life less than 1 h. In general, mRNA species with short half-life were enriched among genes with regulatory functions (transcription factors), whereas mRNA species with long half-life were enriched among genes related to metabolism and structure (extracellular matrix, cytoskeleton). The stability of mRNAs correlated more significantly with the structural features of genes than the function of genes: mRNA stability showed the most significant positive correlation with the number of exon junctions per open reading frame length, and negative correlation with the presence of PUF-binding motifs and AU-rich elements in 3'-untranslated region (UTR) and CpG di-nucleotides in the 5'-UTR. The mRNA decay rates presented in this report are the largest data set for mammals and the first for ES cells.

Defining developmental potency and cell lineage trajectories by expression profiling of differentiating mouse embryonic stem cells.

Biologists rely on morphology, function and specific markers to define the differentiation status of cells. Transcript profiling has expanded the repertoire of these markers by providing the snapshot of cellular status that reflects the activity of all genes. However, such data have been used only to assess relative similarities and differences of these cells. Here we show that principal component analysis of global gene expression profiles map cells in multidimensional transcript profile space and the positions of differentiating cells progress in a stepwise manner along trajectories starting from undifferentiated embryonic stem (ES) cells located in the apex. We present three 'cell lineage trajectories', which represent the differentiation of ES cells into the first three lineages in mammalian development: primitive endoderm, trophoblast and primitive ectoderm/neural ectoderm. The positions of the cells along these trajectories seem to reflect the developmental potency of cells and can be used as a scale for the potential of cells. Indeed, we show that embryonic germ cells and induced pluripotent cells are mapped near the origin of the trajectories, whereas mouse embryo fibroblast and fibroblast cell lines are mapped near the far end of the trajectories. We suggest that this method can be used as the non-operational semi-quantitative definition of cell differentiation status and developmental potency. Furthermore, the global expression profiles of cell lineages provide a framework for the future study of in vitro and in vivo cell differentiation.

Uncovering early response of gene regulatory networks in ESCs by systematic induction of transcription factors.

To examine transcription factor (TF) network(s), we created mouse ESC lines, in each of which 1 of 50 TFs tagged with a FLAG moiety is inserted into a ubiquitously controllable tetracycline-repressible locus. Of the 50 TFs, Cdx2 provoked the most extensive transcriptome perturbation in ESCs, followed by Esx1, Sox9, Tcf3, Klf4, and Gata3. ChIP-Seq revealed that CDX2 binds to promoters of upregulated target genes. By contrast, genes downregulated by CDX2 did not show CDX2 binding but were enriched with binding sites for POU5F1, SOX2, and NANOG. Genes with binding sites for these core TFs were also downregulated by the induction of at least 15 other TFs, suggesting a common initial step for ESC differentiation mediated by interference with the binding of core TFs to their target genes. These ESC lines provide a fundamental resource to study biological networks in ESCs and mice.

Global gene expression profiling reveals similarities and differences among mouse pluripotent stem cells of different origins and strains.

Pluripotent stem cell lines with similar phenotypes can be derived from both blastocysts (embryonic stem cells, ESC) and primordial germ cells (embryonic germ cells, EGC). Here, we present a compendium DNA microarray analysis of multiple mouse ESCs and EGCs from different genetic backgrounds (strains 129 and C57BL/6) cultured under standard conditions and in differentiation-promoting conditions by the withdrawal of Leukemia Inhibitory Factor (LIF) or treatment with retinoic acid (RA). All pluripotent cell lines showed similar gene expression patterns, which separated them clearly from other tissue stem cells with lower developmental potency. Differences between pluripotent lines derived from different sources (ESC vs. EGC) were smaller than differences between lines derived from different mouse strains (129 vs. C57BL/6). Even in the differentiation-promoting conditions, these pluripotent cells showed the same general trends of gene expression changes regardless of their origin and genetic background. These data indicate that ESCs and EGCs are indistinguishable based on global gene expression patterns alone. On the other hand, a detailed comparison between a group of ESC lines and a group of EGC lines identified 20 signature genes whose average expression levels were consistently higher in ESC lines, and 84 signature genes whose average expression levels were consistently higher in EGC lines, irrespective of mouse strains. Similar analysis identified 250 signature genes whose average expression levels were consistently higher in a group of 129 cell lines, and 337 signature genes whose average expression levels were consistently higher in a group of C57BL/6 cell lines. Although none of the genes was exclusively expressed in either ESCs versus EGCs or 129 versus C57BL/6, in combination these signature genes provide a reliable separation and identification of each cell type. Differentiation-promoting conditions also revealed some minor differences between the cell lines. For example, in the presence of RA, EGCs showed a lower expression of muscle- and cardiac-related genes and a higher expression of gonad-related genes than ESCs. Taken together, the results provide a rich source of information about the similarities and differences between ESCs and EGCs as well as 129 lines and C57BL/6 lines. Such information will be crucial to our understanding of pluripotent stem cells. The results also underscore the importance of studying multiple cell lines from different strains when making comparisons based on gene expression analysis.

Esg1, expressed exclusively in preimplantation embryos, germline, and embryonic stem cells, is a putative RNA-binding protein with broad RNA targets.

In our earlier attempt to identify genes involved in the maintenance of cellular pluripotency, we found that KH-domain protein Embryonal stem cell-specific gene 1 (Esg1) showed similar expression patterns to those of Oct3/4 (Pou5f1), whereas the forced repression of Oct3/4 in mouse embryonic stem cells immediately downregulated the expression of Esg1. Here we further confirm this overlap by in situ hybridization and immunohistochemical analyses. Both Esg1 transcript and protein exist in the egg and preimplantation embryos. At embryonic day 3.5, blastocyst stage, however, ESG1 protein was more abundant in the inner cell mass (ICM) than in trophectoderm (TE), whereas Esg1 transcript was detected in both the ICM and the TE, particularly in the polar trophectoderm. The presence of an RNA-binding KH-domain in ESG1 led us to search for and identify 902 target transcripts by microarray analysis of immunoprecipitated ESG1 complex. Interaction of 20 target mRNA with ESG1, including Cdc25a, Cdc42, Ezh2, Nfyc and Nr5a2, was further validated by reverse transcriptase-polymerase chain reaction of the immunoprecipitation material, supporting the notion that ESG1 is an RNA-binding protein which associates with specific target transcripts.

Evolutionary implications of substitution patterns in prolamin genes of Oryza glaberrima (African rice, Poaceae) and related species.

Patterns of sequence variation of nuclear genes encoding 10-kDa and 16-kDa prolamin seed storage proteins were examined in Oryza glaberrima (African rice, Poaceae) and O. barthii and compared to available sequences for the genus to assess potential application of these gene families in evolutionary studies. Sequence variation among species in 10-kDa genes was very low. In contrast, the 16-kDa genes have undergone rapid evolution, displaying a larger number of length and point mutations that in some cases result in frame shift or produce truncated protein or pseudogenes. The proportion of nonsynonymous substitution is high in both genes. Although nonsynonymous mutations did not alter the overall profile of the protein, pronounced shifts in proportions of some amino acids were evident and could have systematic application. The data provide support for a proposed direct evolution of the Asian (O. sativa) and African rice from O. rufipogon and O. barthii, respectively. Patterns of amino acid frequencies of the 10-kDa genes show the distinctness of O. rufipogon and O. longistaminata from the other species. The study underscores the potential application of the prolamin genes as markers from the nuclear genome for evolutionary studies in grasses at different taxonomic levels.