Kazuhiro Aiba

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.

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

BackgroundTarget 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.ResultsTo 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.ConclusionWe 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.

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.

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.

Defining a Developmental Path to Neural Fate by Global Expression Profiling of Mouse Embryonic Stem Cells and Adult Neural Stem/Progenitor Cells

To understand global features of gene expression changes during in vitro neural differentiation, we carried out the microarray analysis of embryonic stem cells (ESCs), embryonal carcinoma cells, and adult neural stem/progenitor (NS) cells. Expression profiling of ESCs during differentiation in monolayer culture revealed three distinct phases: undifferentiated ESCs, primitive ectoderm‐like cells, and neural progenitor cells. Principal component (PC) analysis revealed that these cells were aligned on PC1 over the course of 6 days. This PC1 represents approximately 4,000 genes, the expression of which increased with neural commitment/differentiation. Furthermore, NS cells derived from adult brain and their differentiated cells were positioned along this PC axis further away from undifferentiated ESCs than embryonic stem–derived neural progenitors. We suggest that this PC1 defines a path to neural fate, providing a scale for the degree of commitment/differentiation.

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.

Comparative analysis of oocyte transcript profiles reveals a high degree of conservation among species

Cross-species comparison of gene expression is a powerful approach for discovering genes that have been conserved throughout evolution. Conserved genes are presumably very important in the mechanisms related to the unique molecular functions in oocytes. The objective of this study was to identify genes expressed in the oocyte and conserved across three diverse vertebrate species. We report the global gene expression profiles of Bos taurus and Xenopus laevis oocytes on an NIA mouse development microarray that consists of 60-mer oligonucleotide probes representing more than 20,000 mouse transcripts derived from stem cell, oocyte, and early embryo cDNA libraries. Analysis based on intensity values revealed that 9853 and 10,046 genes are expressed in bovine and Xenopus oocytes respectively. Furthermore, previously published microarray data on preimplantation development in the mouse were used for a comparative analysis of global oocyte gene expression profiles. Interestingly, a substantial proportion of the genes expressed in mouse oocytes is conserved between the three species (74%, 7275 genes). Moreover, functional annotation of these conserved oocyte-expressed genes confirmed that certain functions are conserved among the three species. RNA metabolism and cell cycle were among the over-represented Gene Ontology terms in the biological process category. Finally, a pattern-matching analysis identified 208 conserved maternally expressed genes. Results from these cross-species hybridizations allowed numerous genes expressed in oocytes and conserved between Mus musculus, B. taurus, and X. laevis to be identified. This comparative analysis of oocyte transcript profiles revealed a high degree of conservation among species.

Changes in global gene expression during in vitro decidualization of rat endometrial stromal cells

During the preimplantation phase of pregnancy the endometrial stroma differentiates into decidua, a process that implies numerous morphological changes and is an example of physiological transdifferentiation. Here we show that UIII rat endometrial stromal cells cultured in the presence of calf serum acquired morphological features of decidual cells and expressed decidual markers. To identify genes involved in decidualization we compared gene expression patterns of control and decidualized UIII cells using cDNA microarray. We found 322 annotated genes exhibiting significant differences in expression (>3‐fold, fold discovery rate (FDR) >0.005), of which 312 have not been previously related to decidualization. Analysis of overrepresented functions revealed that protein synthesis, gene expression, and chromatin architecture and remodeling are the most relevant modified functions during decidualization. Relevant genes are also found in the functional terms differentiation, cell proliferation, signal transduction, and matrix/structural proteins. Several of these new genes involved in decidualization (Csdc2, Trim27, Eef1a1, Bmp1, Wt1, Aes, Gna12, and Men1) are shown to be also regulated in uterine decidua during normal pregnancy. Thus, the UIII cell culture model will allow future mechanistic studies to define the transcriptional network regulating reprogramming of stromal cells into decidual cells. J. Cell. Physiol. 222:127–137, 2010.