Gene expression profiles during early differentiation of mouse embryonic stem cells
- Fiona C Mansergh1, 2Email author,
- Carl S Daly†1,
- Anna L Hurley†1,
- Michael A Wride3,
- Susan M Hunter1 and
- Martin J Evans1
© Mansergh et al; licensee BioMed Central Ltd. 2009
Received: 18 July 2008
Accepted: 09 January 2009
Published: 09 January 2009
Understanding the mechanisms controlling stem cell differentiation is the key to future advances in tissue and organ regeneration. Embryonic stem (ES) cell differentiation can be triggered by embryoid body (EB) formation, which involves ES cell aggregation in suspension. EB growth in the absence of leukaemia inhibitory factor (LIF) leads EBs to mimic early embryonic development, giving rise to markers representative of endoderm, mesoderm and ectoderm. Here, we have used microarrays to investigate differences in gene expression between 3 undifferentiated ES cell lines, and also between undifferentiated ES cells and Day 1–4 EBs
An initial array study identified 4 gene expression changes between 3 undifferentiated ES cell lines. Tissue culture conditions for ES differentiation were then optimized to give the maximum range of gene expression and growth. -Undifferentiated ES cells and EBs cultured with and without LIF at each day for 4 days were subjected to microarray analysis. -Differential expression of 23 genes was identified. 13 of these were also differentially regulated in a separate array comparison between undifferentiated ES cells and compartments of very early embryos. A high degree of inter-replicate variability was noted when confirming array results. Using a panel of marker genes, RNA amplification and RT-PCR, we examined expression pattern variation between individual -D4-Lif EBs. We found that individual EBs selected from the same dish were highly variable in gene expression profile.
ES cell lines derived from different mouse strains and carrying different genetic modifications are almost invariant in gene expression profile under conditions used to maintain pluripotency. Tissue culture conditions that give the widest range of gene expression and maximise EB growth involve the use of 20% serum and starting cell numbers of 1000 per EB. 23 genes of importance to early development have been identified; more than half of these are also identified using similar studies, thus validating our results. EBs cultured in the same dish vary widely in terms of their gene expression (and hence, undoubtedly, in their future differentiation potential). This may explain some of the inherent variability in differentiation protocols that use EBs.
ES cells are derived from the inner cell masses of blastocysts and can contribute to all cell types of the embryo proper [1, 2]. The combination of pluripotency and ease of genetic modification has given rise to the revolution in genetic analysis via the use of knockout mice. ES cell pluripotency has also been exploited in vitro; many different cell types can now be generated in culture. Human ES cell lines have been isolated; moreover, ES like cells (iPS cells) can be derived from human skin [3–6]. Therefore, in vitro differentiation protocols for human ES or iPS cells have huge therapeutic potential.
Many in vitro ES differentiation protocols rely on embryoid bodies (EBs); floating aggregates of ES cells which, when grown without LIF, mimic to some extent the early stage embryo, giving rise to precursors of a large number of tissues[7, 8]. Early growth of EBs with LIF favours stem cell renewal and the differentiation of embryonic endoderm , while removal of LIF allows the generation of precursors representative of all three germ layers . To date, protocols for deriving neural, haematopoietic, muscle, bone, pancreatic, hepatic and many other precursor and mature cell types from ES cells in culture [11–13] have been developed, many of which still use EBs as a starting point. EBs can be allowed to differentiate spontaneously from cell suspension or can be formed from a defined cell number using 'hanging drops'. Following aggregation, culture is often allowed to proceed spontaneously for 3–4 days, followed by the addition of factors that promote differentiation of specific precursor types; for example, retinoic acid may be added to promote neuronal specification [14, 15], although, today, more efficient neural differentiation can be achieved in chemically defined medium or via adherent monoculture in the presence of FGF [16–18]. Subsequent growth, followed by disassociation and plating on adherent surfaces, permits the derivation of terminally differentiated cell types.
ES differentiation can provide abundant, partially synchronised sources of transient embryonic precursor types that are present only in very limiting quantities in vivo. Moreover, EBs represent a good model for examining the events of early embryogenesis, as the formation of a pro-amniotic cavity and the expression of markers of early differentiation, for example, are often mimicked by EBs . However, final cell populations are usually heterogeneous, percentages of desired cell types arising often vary from one experimental replicate to another [14, 15, 10]. Reasons for this variation are not hard to identify; differentiation is sensitive to glucose concentration, serum quality, amino acids, growth factors, extracellular matrix proteins, pH, osmolarity, passage number and the identity of the ES cell line used . The existence of ES cell derived chimeras and ensuing mouse lines demonstrates the ability of ES cells to differentiate into all adult cell types. However, we are at present unable to generate the full complement of adult cell types from ES cells in vitro. Greater understanding of gene expression during early differentiation may allow more precise direction of ES cell differentiation and will also widen understanding of early embryonic development.
Dissecting the events of early differentiation has been aided by the development of microarray technology, which allows the examination of global gene expression changes. We have used microarray technology to examine variation between 3 undifferentiated ES cell lines. We then optimised aggregation methods, EB size, and serum concentrations and carried out array analysis using day 1–4 EBs in the presence and absence of LIF. We identified 23 differentially regulated genes, some of which have known roles in early development. However, given lower than expected confirmation rates, and lack of reproducibility in stem cell derived arrays [20–22], we tested the replicability of gene expression patterns arising from individual EBs growing in the same culture, using RNA amplification. We found a high level of variation in gene expression patterns, even between EBs from the same culture dish. This variability may provide explanations for the difficulties involved in obtaining pure cultures of differentiated progeny from EB based protocols.
Tissue Culture and cell lines
This study utilised 3 ES cell lines, IMT11 (derived from 129 mice), HM1 (which is Hprt negative) and SMHBl6 (derived from C57Bl6/J mice). The IMT11 line was selected for all investigations involving differentiation, as it is not genetically modified and is better characterised than SMHBl6. IMT11 cells showed the highest percentage of diploid cells after karyotyping and have been tested for germline transmission.
1) Undifferentiated ES cells were maintained at 37°C in a humidified atmosphere with 5% CO2 on 0.1% gelatin in DMEM, with 2 mM L-glutamine, 50 U/ml penicillin, 50 μg/ml streptomycin (all from Gibco™, Invitrogen Ltd, Paisley, Renfrewshire, UK), 10-4 M β-Mercaptoethanol (Merck KGaA, 64293 Darmstadt, Germany), 10-3 U/ml murine LIF (ESGRO™, Invitrogen, Ltd, Paisley, Renfrewshire, UK), 10% FBS (foetal bovine serum) and 10% NBS (newborn bovine serum) (selected batches, PAA Laboratories GmbH, Linz, A-4020 Austria). All undifferentiated ES cell lines were karyotyped using standard protocols in order to test that the majority of cells showed a normal diploid chromosome number (40XY) prior to differentiation.
2) EB generation: a semi-confluent 100 mM dish of ES cells was trypsinized (0.25% trypsin/EDTA, Invitrogen), followed by trituration in additional ES medium to achieve a single cell suspension. ES medium was prepared as above for + LIF EBs, and without LIF for -LIF differentiations. Cells were counted using a haemocytometer and cell density was adjusted appropriately to the required number of cells per 10 μl. A multichannel pipette (Finnpipette 5–50 μl) was used to deposit approximately 200 10 μl drops on the floor of a 140 mm bacteriological dish (Sterilin). A smaller plate was filled with 1–2 mls PBS and placed in the lid of the bacteriological dish. The plate was inverted and incubated overnight to allow the EBs to aggregate. The following day dishes were righted and flooded with 20 mls of the appropriate differentiation medium, then grown in suspension culture until harvesting. EBs were formed by aggregation of 125, 250, 500, 750 and 1000 cells per 10 μl to determine optimal size. A 50% FCS, 50% NCS mix was prepared from ES batch tested serum samples (PAA), then added to serum free ES medium, at 0%, 5%, 10%, 15% and 20% final volumes. Hanging drops containing either 750 or 1000 cells were generated in order to test serum concentrations.
We also tested EBs that were allowed to aggregate randomly following dissociation. For brevity, we refer to these as "random EBs later in the text". Random aggregation gives rise to greater numbers of EBs per dish, that are of much more variable sizes and shapes than those generated using the "hanging drop" method. For random aggregation of EBs, 1 × 106 cells were suspended in 10 mls media in a 100 mm bacteriological dish and allowed to aggregate spontaneously. Serum quantities were varied as above. Photographs were taken of EBs at day 4 of differentiation without LIF, at varying sizes and serum concentrations and saved in .tif format. Photographs were analysed and EB diameters measured using Scion Image (Scion Corporation). Diameter measurements were used to calculate radius and volume.
RNA extraction and amplification
1) For array analysis, undifferentiated ES cells were washed twice with DPBS-A, treated with trypsin-EDTA (Invitrogen Ltd., Paisley, UK) and counted. Appropriate cell numbers were pelleted by centrifugation at 1000 g. The QIAgen RNeasy™ Midi Kit (Qiagen Ltd., Crawley, Sussex, UK) was used according to the manufacturer's protocol for RNA extraction, followed by OD 260/280 spectrophotometry (Camspec, Cambridge, UK) and gel electrophoresis using dissociating conditions (NorthernMax™ buffers; Ambion, Huntingdon, UK), used according to the manufacturer's protocol, to check RNA concentration and integrity.
2) For RNA extraction from individual embryoid bodies, the mini RNA isolation kit (Zymo) was used, followed by amplification using the RNA amplification kit (Arcturus), both according to the manufacturer's protocol. Following quantification by UV spectrophotometer (Camspec), amplified RNA samples were DNase treated using the Turbo DNAfree kit (Ambion) and reverse transcribed using the random hexamer protocol of the Superscript First Strand Synthesis System for RT-PCR (Invitrogen). (Note, a nuclease step is included in the Arcturus kit to remove residual DNA, but is not entirely sufficient to ensure blank negative control lanes, hence the extra DNase step).
3) For RNA extraction followed by RT-PCR, embryoid bodies were spun down, washed with PBS, then spun down again. EBs were treated with 0.25% trypsin/EDTA and washed with PBS a second time when extracting RNA from larger Day 3–4 EBs 100 ul of EBs+ residual PBS was resuspended in 1 ml TRIzol (Invitrogen); samples were homogenized by pipetting up and down using a 1 ml micropipette. RNA was extracted according to the manufacturer's protocol, quantified, DNase treated and reverse transcribed as above, except that the oligo-dT supplied with the kit was spiked with a 1/10000 dilution of an 18S rRNA gene specific primer, in order to allow 18S rRNA (which doesn't have a poly A tail) to be used as a housekeeping control. Where starting RNA quantity was low, we tripled the amount of RNA used and therefore the reaction volumes per sample. RT reactions were diluted with nuclease free water (Ambion) to 50 ul before PCR analysis. A "no RT" control corresponding to each sample was also produced for all RT-PCR experiments described in this paper; these were treated in exactly the same way as the samples except that reverse transcriptase was not added.
4) For RNA extraction prior to array analysis, EBs were spun down, washed, and trypsinized as above. The QIAgen RNeasy™ Midi Kit (Qiagen Ltd., Crawley, Sussex, UK) was used according to the manufacturer's protocol for RNA extraction. RNA samples were quantified and checked for quality as described in 1) above.
Cy labelling for array analysis
10 μg of total RNA was labelled with either Cy3 or Cy5 dyes using the CyScribe labelling system (GE Healthcare, Chalfont St. Giles, Bucks, UK), according to the manufacturer's protocol. 1 μl of labelled cDNA was combined with 2 ul 50% glycerol/50% TE mixture, run on a microscope slide sized, 1.5% agarose gel (mould manufactured in-house) and scanned using a GeneTac LS IV scanner (Genomic Solutions, Huntingdon, Cambs. UK). Control and experimental samples were then combined and prepared for hybridization.
Array slides were incubated in prehybridization buffer for 1 hour at 42°C (50% formamide, 5 × SSC, 0.1% SDS, 1% BSA). Targets were dried down via vacuum centrifugation then resuspended in 50 μl hybe solution (49.9% de-ionised formamide, 49.9% 20 × SSC, 0.2% SDS) with added 1 μl Cot1 DNA and 1 μl poly A oligo as blocking agents, heated to 95°C for 5 minutes and then added to the face of one slide. The printed face of the second slide of the pair was then placed face to face with the first, using the same probe. Slide pairs were then placed in a humidified container and incubated for 24–48 hours at 42°C. Following hybridization, slides were washed once in Wash solution 1 (1× SSC, 2% SDS, filtered autoclaved ddH2O) for 20 minutes, then twice in Wash solution 2 (0.1× SSC, 0.2% SDS, filtered autoclaved distilled deionised H2O (ddH2O) for 20 minutes each. Slides were dipped in nuclease free filtered water, then spray dried, finally, the backs of the slides were cleaned with filtered autoclaved ddH2O, then wiped with 100% EtOH, then wiped dry and scanned.
Scans were carried out at 12.5 μm, using the averaging setting (GeneTac LSIV scanner, (Genomic Solutions, Cambridgeshire, UK)). It is possible to carry out quick draft scans using this scanner. Gain and black settings, which affect image intensity and background, were varied slightly in order to optimize the signal/noise ratio for each channel and each slide before proper scans were initiated.
In adherence with MIAME standards , all data sets have been submitted to the GEO database http://www.ncbi.nlm.nih.gov/geo/ and are fully MIAME compliant. Undifferentiated ES cell comparison data are described in GSE8625, while EB differentiation comparisons are described in GSE8766. The results from these experiments were compared with an array comparison (described elsewhere) between undifferentiated ES colonies and microdissected inner cell masses from embryonic day 3.5, 4.5, and delayed blastocysts and microdissected day 5.5 and day 6.5 embryonic ectoderm (GSE8881, ).
Undifferentiated ES cell array experiments (GSE8625) were carried out using NIA 15K slides printed in the Cardiff Microarray Facility. These consisted of 17136 spots, printed in 12 × 4 mini grids, each with 17 rows and 21 columns and are described more fully in GEO, platform accession number GPL5530. EB differentiation comparisons (GSE8766) were carried out using NIA 15K slides printed by the HGMP. These consist of 2 slides, 17280 spots per slide, including control spots, empties and landmarks. 2 duplicate spots are printed per slide for each clone. These slides are described more fully in GEO, platform accession number GPL5735. The array slides used for GSE8881 were also printed in Cardiff and contain 16128 spots, printed in 12 × 4 mini grids, each with 16 rows and 21 columns. The GEO platform accession number for these is GPL5771. Differences in spot number relate to differences in the number of landmarks and other controls printed; the core NIA 15K set is present on each of these platforms
Experimental design, Image analysis
Undifferentiated ES cell arrays were carried out using samples from undifferentiated IMT11, SMHBL6.3 and HM1 ES cells. EB arrays were carried out using samples from Day0 (undifferentiated), Day1, Day2, Day3 and Day 4 differentiated EBs, generated in the presence and absence of LIF (9 samples in total). Each experimental sample was hybridized on a slide with a pooled control derived from an equal amount of all experimental samples. ES arrays were repeated 12 times (this print run contained no duplicate spots). EB arrays were repeated twice. Duplicate spots were present on each array, giving four repeats total for each gene. Scanned images were stored and filtered, then analysed using ImaGene™ 5.5 (BioDiscovery). This series of array experiments was carried out using the NIA15K set . ESTs comprising this set were isolated from a variety of embryonic stages and tissues. Arrays were hybridised with fluor switching, in order to counteract any issues of dye bias that may have arisen from direct labelling.
Array analyses for GEO entries GSE8625 and GSE8766 were carried out similarly to work described previously [26, 27]. Output files from ImaGene were saved and analysed in MS Excel spreadsheet format. Each channel from each repetition was normalised via division by the mean intensity value. These data were collated and formatted for Significance Analysis of Microarrays (SAM; http://www-stat.stanford.edu/~tibs/SAM). Genes that showed a fold change of 2 or above and were statistically significant above a delta value of 0.5 (which denotes an error rate of 5%) were selected for further appraisal. SAM, however, has a disadvantage; all replicates have to be in the same order, precluding any data filtration. Hence, we supplemented SAM analysis with a second analysis method as follows:
Following normalisation, we used approximately 700 blank spots per slide to calculate a mean background value + 2 standard deviations of that background value, for each channel. Genes that fell below this cut-off in BOTH control and experimental channels were removed, along with genes with a fold change of < 2. This filtered gene list was compared with that from SAM; genes appearing as differentially regulated using both methods were deemed significant (fold change > 2, above background + 2 SD in at least one channel, delta value of 0.5).
A) Marker/Housekeeping genes (note: Brachyury, Beta-actin, Hprt and Gapdh also came up in the arrays)
PCR fragment size (bp)
B) Primers, Confirmed and interesting genes:
PCR fragment size (bp)
Bioinformatic analysis of confirmed genes was also carried out in order to identify putative function in differentiation. Accession numbers were used to comprehensively search the NCBI databases http://www.ncbi.nlm.nih.gov. Data compiled from various NCBI databases (UniGene, Homologene, OMIM, LocusLink, PubMed etc.) are shown in Figures 1 and 2. Selected marker genes were chosen on the basis of a known role in maintenance of pluripotency, or early development, or because they are well known housekeeping genes. Primers were designed for these from the appropriate reference sequence, using Primer3 as above. Datasets derived from ES cell differentiation subtractive EST library studies [15, 29] or microarrays [30–38] were also searched for the presence of the genes we confirmed in this study (Figures 1 and 2).
RT PCR array confirmations
We used semi-quantitative RT-PCR to confirm differentially regulated genes. As most significant genes appeared to be differentially regulated at various different stages of differentiation, we tested all genes against Day 0, and Day1–4 EBs, + and - LIF. Three biological replicates of 1000 cell Day 1–4 EBs and 3 replicates of 750 cell EBs (+ and - LIF) were used for confirmations, including No RT controls. In order to be deemed confirmed, we asked that a gene demonstrate a consistent expression pattern in at least 4 out of 6 PCR tests. The number of cycles required for minimum visibility was identified and PCRs were optimised at the Tm indicated for each primer pair, such that only single bands appeared. PCR bands obtained for housekeeping controls were of even intensity at minimum visibility cycles before testing other genes. PCRs were carried out in 20 μl volumes using 0.025 μmol of each dNTP (Invitrogen) and 1 unit Taq DNA polymerase (Sigma) per reaction. PCRs were usually carried out using PCR buffer (Sigma) containing 15 mM MgCl2 although magnesium concentrations, along with annealing temperature and cycle number, were varied where necessary in order to optimise results. Final PCR conditions are indicated for each gene in the appropriate figures (see Results). Primer sequences are given for all genes tested in Table 1 and Additional File 1. In order to confirm that semi-quantitative PCR was indeed representative of expression pattern changes, we studied the expression patterns of two confirmed genes, Hspa8 and BG063737, by quantitative PCR. Expression pattern changes achieved by q-PCR were very similar to the methods described above. Methods and results are given in Additional File 2.
Undifferentiated ES cell arrays
Optimisation of embryoid body differentiation
Array results from Day 1–4 EBs, grown in both the presence and the absence of LIF, were compared at all stages with the pooled control (equal amounts of Day 0–4 ES cells + EBs, both with and without LIF), and also with each other, when analysing the arrays. 128 ESTs selected for follow-up were differentially regulated in all 4 repeats and were statistically significant according to SAM. Many ESTs were differentially regulated in more than one comparison. Removal of redundancy (more than one EST mapping to different parts of the same gene) and elimination of those ESTs where sequence quality was insufficient for primer design resulted in a master list of 104 genes. Notably, beta-actin and Gapdh were differentially regulated, so 18S rRNA was used as a housekeeping gene instead.
We then compared this dataset with that from another of our array experiments, which compared undifferentiated IMT11 ES cells with embryonic inner cell masses (the tissue from which ES cells are derived) and day 5.5 and day 6.5 embryonic ectoderm (see Figures 1 and 2, GSE8881, methods described ), Thirteen of the above 26 genes were also up- or down-regulated in this study, indicating their importance in early development and differentiation. Fourteen of the known genes were also noted in other array or subtractive EST library studies using material from different ES cell differentiation experiments (see Figures 1 and 2). Interestingly, two of the four genes that varied between undifferentiated cell lines were also identified by this study.
Amplification from single EBs
We have shown that undifferentiated mouse ES cells, regardless of strain of origin, or absence of Hprt, maintain a very invariant expression pattern under identical culture conditions. Given genetic variation between mouse strains of origin (129Sv/Ev and C57Bl6/J) and absence of the Hprt gene in one line, this result is perhaps surprising, but would imply that the culture conditions that maintain pluripotency may be quite restrictive in terms of gene expression pattern. Hprt was not expressed in HM1 cells as expected. The relevance of Pter differential regulation is unclear, unless this gene is functionally associated with the Hprt pathway. Mki67 is a marker of proliferating cells; its upregulation in IMT11 cells might suggest a higher growth rate, although this has not been tested. Scmh1 may be involved in the meiotic prophase of male germ cells; ES cells are similar to germ cells, but the significance of the differential expression is unclear. We should note that, apart from HM1, expression differences noted were minor. We selected IMT11 cells for further study.
EB growth optimisation shows that EBs of 500–1000 cells reach a size plateau at 20% serum, suggesting that size constraints apply to growth once a certain size is reached. This may be related to the decreased rate of gas and nutrient diffusion to cells in the centres, but EBs are thought to be loosely packed and should therefore not suffer from a buildup of toxic waste products [40, 41] EBs created using cell numbers of 1–1000 were previously found to attain a maximum of approximately 30,000 cells on Day 12 of differentiation, regardless of culture methods used or initial cell number . Alternatively, growth constraints may arise from the absence of some necessary developmental cue, signalling during cavitation or from decreased cell division rates upon initiation of differentiation, which may happen sooner with a larger initial size [41, 42].
Variation of initial size and serum concentration also influences gene expression patterns after 4 days of differentiation. We selected an initial size of 1000 cells and a serum concentration of 20% for our array experiments, on the grounds that these conditions were optimal for growth (and therefore RNA yield), gave the widest range of marker gene expression (implying a greater range of differentiated progeny) and were likely to bear more similarity to early embryonic differentiation.
Microarrays have identified 23 genes, which are differentially regulated during ES cell differentiation; a further 3, while not demonstrating the expression patterns predicted by the array, showed dramatic down-regulation upon induction of differentiation without LIF and are therefore discussed here. Detailed synopses of functional information available is given in Figure 1; a number of functional studies are also of interest [43–59].
13 of these genes are implicated in the early development of pre- and post- implantation embryos, while 14 are highlighted by other studies and/or preferentially expressed in embryonic EST libraries. Given the variation in ES cell lines, differentiation protocols and timepoints studied by others, this overlap is reasonably significant [15, 29–38]. The functions of 7 ESTs and 1 hypothetical protein are presently uncharacterised.
The remaining genes fall into a number of different classes. Hspca (a.k.a Hspc1, Hsp89a, Hsp90a) and Hspa8 are immediately downregulated on induction of differentiation. Both are constitutively expressed heatshock proteins, a class of proteins commonly upregulated in response to cellular stress (OMIM 140571) and also implicated in ES differentiation . Heatshock proteins can function as chaperones; ie, they assist in the correct folding of nascent proteins. Ssr2, another protein downregulated upon onset of differentiation, is also associated with nascent polypeptides, although it forms part of a complex that recognises signal peptides, resulting in transport of the relevant protein across the endoplasmic reticulum. Hspca has been implicated in epigenetic modifications of chromatin (disruption during Drosophila development results in heritable morphogenic alterations) and in increased carcinogenic invasiveness . Pluripotent cells, such as ES cells, demonstrate uniquely dynamic chromatin, "breathing chromatin", which allows availability of a large proportion of the genome for immediate transcriptional activity . Perhaps reductions in Hspca expression during differentiation may result in changes to chromatin structure that result in the tighter histone binding characteristic of differentiated cells and their committed precursors.
Notably, a related class of 5 genes also shows downregulation upon initiation of differentiation; Eif4a1, Rplp1, Rpl8, Rps2 and Rps28. Eif4a1 encodes a translation initiation factor that could be involved in selective regulation of protein expression, while the remaining four genes are structural constituent of the ribosome. Expression of Eif4a1 is associated with increased metastasis in certain cancers , can stimulate cell growth in hepatoma cells  and is enriched in populations of human ES cells . We have previously noted the upregulation of similar genes involved in translation in response to neural differentiation of ES cells ; however, this upregulation was noted after 7 days of EB differentiation, following the addition of retinoic acid on Day 4. It is possible that ribosomal constituents may play a role in complex changes that occur in gene expression in response to changing signals at different stages of differentiation.
A third class of genes comprises those that are involved in energy metabolism and other ubiquitous metabolic processes (Atp5b, Cox4i1, Ftl1, transketolase, Hnrpk). Energy and metabolic requirements may change during differentiation; alternatively, these genes may have alternative functions in these processes. Perhaps surprisingly, only 3 genes, Brachyury, Aurkb and Dppa5 have identified roles in the maintenance of pluripotency, epigenetic remodelling and early development (see Figure 1); we would have expected this class of gene to be better represented. However, larger EBs maintain a core of pluripotent cells that maintain LIF expression and also, therefore, continue to express genes such as Oct4, Nanog and Rex1 for some time. This is obvious from our PCRs of these marker genes; expression can dip a little on Day 1 but normally resumes by Day3–4 of differentiation. Other genes, such as Laminin, are known to display differential expression during early development; laminin was among the 104 genes selected for PCR confirmation. However, expression varied so much between PCR replicates that there was no consistent pattern of variation. This finding was common to a majority of the genes tested; this would not rule out their importance to early developmental processes, however, it would suggest that those genes that did show consistent expression changes may indeed be very important to early development; we were stringent in our selection of differentially regulated genes.
The identification of 26 genes which alter in expression pattern in ES cell differentiation will aid the understanding of early development and in vitro differentiation. However, there are important implications of the low confirmation rates noted in this paper, and of the variability in marker gene expression profiles between individual EBs of homogeneous cell number and size, which originated from two tissue culture dishes. Firstly, ES differentiation protocols using EBs as a starting point may always generate very variable results, despite attempts to standardize them such as using a uniform EB starting size. This would necessarily limit their therapeutic use; development of differentiation protocols that yield more uniform populations of progenitors without the use of EBs may be preferable (such as the majority of current neural differentiation methods)[17, 18]. On the other hand, the ability of EBs to generate a wide variety of precursors that could later be selected for subtypes of choice may be an advantage in certain circumstances. Secondly, this finding would imply that the lack of directional orientation in EBs as compared to early embryos is a source of chaotic variability.
Stem cells provide, potentially, an unparalleled opportunity for treatment of any number of degenerative conditions. A deeper knowledge of stem cell differentiation and of signalling pathways activated therein, will increase our ability to direct the differentiation of stem cells in vitro. This study has advanced our knowledge of early ES cell differentiation in several key respects. Firstly, we have demonstrated that regardless of underlying genetic variation, the constraints of ES cell pluripotency seem to maintain a relatively invariant gene expression profile. Secondly, we have developed optimised tissue culture conditions that allow the widest range of differentiation potential. Thirdly, we have identified genes that are implicated by this study and others in the complex processes of early development; furthermore, given the number of biological repetitions and the EB size variation used when confirming these genes, we would regard our results as robust. Perhaps most importantly, when investigating low confirmation rates, we have shown that gene expression patterns of individual EBs vary markedly from each other, even when grown from the same number of starting cells, in the same culture. This implies that differentiation protocols involving EBs may always yield varying proportions of different cell types, no matter how rigorously conditions are controlled between different experimental replicates, and may suggest that a move away from EB based differentiations may be warranted, where possible.
This work was funded by BBSRC, the Wales Gene Park, the Royal Society and HEFCW. We would like to thank Steve Turner, Steffan Adams and Vicky Workman (Cardiff University Array facility) for their help with microarraying. FCM is currently supported by RP Ireland, Fighting Blindness.
- Evans MJ, Kaufman MH: Establishment in culture of pluripotential cells from mouse embryos. Nature. 1981, 292 (5819): 154-6. 10.1038/292154a0.View ArticlePubMedGoogle Scholar
- Smith AG: Embryo-derived stem cells: of mice and men. Annu Rev Cell Dev Biol. 2001, 17: 435-62. 10.1146/annurev.cellbio.17.1.435.View ArticlePubMedGoogle Scholar
- Thomson JA, Itskovitz-Eldor J, Shapiro SS, Waknitz MA, Swiergiel JJ, Marshall VS, Jones JM: Embryonic stem cell lines derived from human blastocysts. Science. 1998, 282 (5391): 1145-7. 10.1126/science.282.5391.1145.View ArticlePubMedGoogle Scholar
- Wernig M, Meissner A, Foreman R, Brambrink T, Ku M, Hochedlinger K, Bernstein BE, Jaenisch R: In vitro reprogramming of fibroblasts into a pluripotent ES-cell-like state. Nature. 2007, 448 (7151): 318-24. 10.1038/nature05944.View ArticlePubMedGoogle Scholar
- Lowry WE, Richter L, Yachechko R, Pyle AD, Tchieu J, Sridharan R, Clark AT, Plath K: Generation of human induced pluripotent stem cells from dermal fibroblasts. Proc Natl Acad Sci USA. 2008, 105 (8): 2883-8. 10.1073/pnas.0711983105.PubMed CentralView ArticlePubMedGoogle Scholar
- Park IH, Zhao R, West JA, Yabuuchi A, Huo H, Ince TA, Lerou PH, Lensch MW, Daley GQ: Reprogramming of human somatic cells to pluripotency with defined factors. Nature. 2008, 451 (7175): 141-6. 10.1038/nature06534.View ArticlePubMedGoogle Scholar
- Leahy A, Xiong JW, Kuhnert F, Stuhlmann H: Use of developmental marker genes to define temporal and spatial patterns of differentiation during embryoid body formation. J Exp Zool. 1999, 284 (1): 67-81. 10.1002/(SICI)1097-010X(19990615)284:1<67::AID-JEZ10>3.0.CO;2-O.View ArticlePubMedGoogle Scholar
- O'Shea KS: Self-renewal vs. differentiation of mouse embryonic stem cells. Biol Reprod. 2004, 71 (6): 1755-65. 10.1095/biolreprod.104.028100.View ArticlePubMedGoogle Scholar
- Hamazaki T, Oka M, Yamanaka S, Terada N: Aggregation of embryonic stem cells induces Nanog repression and primitive endoderm differentiation. J Cell Sci. 2004, 117 (Pt 23): 5681-6. 10.1242/jcs.01489.View ArticlePubMedGoogle Scholar
- Desbaillets I, Ziegler U, Groscurth P, Gassmann M: Embryoid bodies: an in vitro model of mouse embryogenesis. Exp Physiol. 2000, 85 (6): 645-51. 10.1017/S0958067000021047.View ArticlePubMedGoogle Scholar
- Mansergh FC, Wride MA, Rancourt DE: Neurons from stem cells: Implications for understanding nervous system development and repair. Biochem Cell Biol. 2000, 78: 613-628. 10.1139/bcb-78-5-613.View ArticlePubMedGoogle Scholar
- Mansergh FC, Wride MA, Rancourt DE: Neurons, Stem Cells and Potential Therapies. Stem Cell Handbook. Edited by: Sell S. 2003, Humana Press, Totowa, NewJersey, ISBN: 1-58829-113-8Google Scholar
- Wobus AM, Boheler KR: Embryonic stem cells: prospects for developmental biology and cell therapy. Physiol Rev. 2005, 85: 635-678. 10.1152/physrev.00054.2003.View ArticlePubMedGoogle Scholar
- Bain G, Kitchens D, Yao M, Huettner JE, Gottlieb DI: Embryonic stem cells express neuronal properties in vitro. Dev Biol. 1995, 168 (2): 342-57. 10.1006/dbio.1995.1085.View ArticlePubMedGoogle Scholar
- Bain G, Mansergh FC, Wride MA, Hance JE, Isogawa A, Rancourt SL, Ray WJ, Yoshimura Y, Tsuzuki T, Gottlieb DI, Rancourt DE: ES cell neural differentiation reveals a substantial number of novel ESTs. Funct Integr Genomics. 2000, 1 (2): 127-39. 10.1007/s101420000014.View ArticlePubMedGoogle Scholar
- Okabe S, Forsberg-Nilsson K, Spiro AC, Segal M, McKay RD: Development of neuronal precursor cells and functional postmitotic neurons from embryonic stem cells in vitro. Mech Dev. 1996, 59 (1): 89-102. 10.1016/0925-4773(96)00572-2.View ArticlePubMedGoogle Scholar
- Ying QL, Stavridis M, Griffiths D, Li M, Smith A: Conversion of embryonic stem cells into neuroectodermal precursors in adherent monoculture. Nat Biotechnol. 2003, 21 (2): 183-6. 10.1038/nbt780.View ArticlePubMedGoogle Scholar
- Bouhon IA, Joannides A, Kato H, Chandran S, Allen ND: Embryonic stem cell-derived neural progenitors display temporal restriction to neural patterning. Stem Cells. 2006, 24 (8): 1908-13. 10.1634/stemcells.2006-0031.View ArticlePubMedGoogle Scholar
- Coucouvanis E, Martin GR: BMP signaling plays a role in visceral endoderm differentiation and cavitation in the early mouse embryo. Development. 1999, 126 (3): 535-46.PubMedGoogle Scholar
- Evsikov AV, Solter D: Comment on " 'Stemness':Transcriptional profiling of embryonic and adult stem cells" and " A stem cell molecular signature" (II). Science. 2003, 302: 393-10.1126/science.1082380.View ArticlePubMedGoogle Scholar
- Ivanova NB, Dimos JT, Schaniel C, Hackney JA, Moore KA, Ramalho-Santos M, Yoon S, Matsuzaki Y, Mulligan RC, Melton DA, Lemischka IR: Comment on " 'Stemness':Transcriptional profiling of embryonic and adult stem cells" and " A stem cell molecular signature". Science. 2003, 302: 393d-10.1126/science.1088249.View ArticleGoogle Scholar
- Fortunel NO, Otu HH, Ng HH, Chen J, Mu X, Chevassut T, Li X, Joseph M, Bailey C, Hatzfeld JA, Hatzfeld A, Usta F, Vega VB, Long PM, Libermann TA, Lim B: Comment on " 'Stemness':Transcriptional profiling of embryonic and adult stem cells" and " A stem cell molecular signature" (II). Science. 2003, 302: 393b-10.1126/science.1086384.View ArticleGoogle Scholar
- Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, Gaasterland T, Glenisson P, Holstege FC, Kim IF, Markowitz V, Matese JC, Parkinson H, Robinson A, Sarkans U, Schulze-Kremer S, Stewart J, Taylor R, Vilo J, Vingron M: Minimum information about a microarray experiment (MIAME)-toward standards for microarray data. Nature Genetics. 2001, 29: 365-371. 10.1038/ng1201-365.View ArticlePubMedGoogle Scholar
- Hunter SM, Mansergh FC, Evans MJ: Optimization of minuscule samples for use with cDNA microarrays. J Biochem Biophys Methods. 2008, 70 (6): 1048-1058. 10.1016/j.jprot.2007.11.011.View ArticlePubMedGoogle Scholar
- Kargul GJ, Dudekula DB, Qian Y, Lim MK, Jaradat SA, Tanaka TS, Carter MG, Ko MS: Verification and initial annotation of the NIA mouse 15K cDNA clone set. Nat Genet. 2001, 28 (1): 17-8. 10.1038/88206.PubMedGoogle Scholar
- Mansergh FC, Wride MA, Walker VE, Adams S, Hunter SM, Evans MJ: Gene expression changes during cataract progression in Sparc null mice:differential regulation of mouse globins in the lens. Mol Vis. 2004, 10: 490-511.PubMedGoogle Scholar
- Mansergh FC, Wells T, Elford C, Evans SL, Perry MJ, Evans MJ, Evans BA: Osteopenia in Sparc (osteonectin)-deficient mice: characterisation of phenotypic determinants of femoral strength and changes in gene expression. Physiol Genomics. 2007, 32 (1): 64-73. 10.1152/physiolgenomics.00151.2007.PubMed CentralView ArticlePubMedGoogle Scholar
- Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. Humana Press, Totowa, NJ, 365-386.
- Sharov AA, Piao Y, Matoba R, Dudekula DB, Qian Y, VanBuren V, Falco G, Martin PR, Stagg CA, Bassey UC, Wang Y, Carter MG, Hamatani T, Aiba K, Akatsu H, Sharova L, Tanaka TS, Kimber WL, Yoshikawa T, Jaradat SA, Pantano S, Nagaraja R, Boheler KR, Taub D, Hodes RJ, Longo DL, Schlessinger D, Keller J, Klotz E, Kelsoe G, Umezawa A, Vescovi AL, Rossant J, Kunath T, Hogal BL, Curci A, D'Urso M, Kelso J, Hide W, Ko MS: Transcriptome analysis of mouse stem cells and early embryos. PLOS Biology. 2003, 1 (3): 410-419. 10.1371/journal.pbio.0000074.View ArticleGoogle Scholar
- Kelly DL, Rizzino A: DNA microarray analyses of genes regulated during the differentiation of embryonic stem cells. Mol Reprod Dev. 1999, 56: 113-123. 10.1002/(SICI)1098-2795(200006)56:2<113::AID-MRD1>3.0.CO;2-Q.View ArticleGoogle Scholar
- Gunji W, Kai T, Sameshima E, Iizuka N, Katagi H, Utsugi T, Fujimori F, Murakami Y: Global analysis of the expression patterns of transcriptional regulatory factors in formation of embryoid bodies using sensitive oligonucleotide microarray systems. Biochem Biophys Res Comm. 2004, 325: 265-275. 10.1016/j.bbrc.2004.10.025.View ArticlePubMedGoogle Scholar
- Gissel C, Voolstra C, Doss MX, Koehler C, Winkler J, Hescheler J, Sachinidis A: An optimised embryonic stem cell model for consistent gene expression and developmental studies. Thromb Haemostat. 2005, 94: 719-27.Google Scholar
- Palmqvist L, Glover CH, Hsu L, Lu M, Bossen B, Piret JM, Humhries RK, Helgason CD: Correlation of murine embryonic stem cell expression profiles with functional measures of pluripotency. Stem Cells. 2005, 23: 663-680. 10.1634/stemcells.2004-0157.View ArticlePubMedGoogle Scholar
- Loh Y-H, Wu Q, Chew J-L, Vega VB, Zhang W, Chen X, Bourque G, George J, Leong B, Liu J, Wong K-Y, Sung KW, Lee CWH, Zhao X-D, Chiu K-P, Lipovich L, Kuznetsov VA, Robson P, Stanton LW, Wei C-L, Ruan Y, Lim B, Ng H-H: The Oct4 and Nanog transcription network regulates pluripotency in mouse embryonic stem cells. Nature Genetics. 2006, 38 (4): 431-440. 10.1038/ng1760.View ArticlePubMedGoogle Scholar
- Pritsker M, Ford NR, Jenq HT, Lemischka IR: Genomewide gain-of function genetic screen identifies functionally active genes in mouse embryonic stem cells. PNAS. 2006, 103 (18): 6946-6951. 10.1073/pnas.0509861103.PubMed CentralView ArticlePubMedGoogle Scholar
- Hirst CE, Ng ES, Azzola L, Voss AK, Thomas T, Stanley EG, Elefanty AG: Transcriptional profiling of mouse and human ES cells identifies SLAIN1, a novel stem cell gene. Developmental Biology. 2006, 293: 90-103. 10.1016/j.ydbio.2006.01.023.View ArticlePubMedGoogle Scholar
- Sene KH, Porter CJ, Palidwor G, Perez-Iratxeta C, Muro EM, Campbell PA, Rudnicki MA, Andrade-Navarro MA: Gene function in early mouse embryonic stem cell differentiation. BMC Genomics. 2007, 8: 85-10.1186/1471-2164-8-85.PubMed CentralView ArticleGoogle Scholar
- Sekkai D, Gruel G, Herry M, Moucadel V, Constantinescu SN, Albagli O, Tronik-Le-Roux D, Vainchenker W, Bennacour-Griscelli A: Microarray analysis of Lif/Stat3 transcriptional targets in embryonic stem cells. Stem Cells. 2005, 23: 1634-1642. 10.1634/stemcells.2005-0182.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle Scholar
- Gassmann M, Fandrey J, Bichet S, Wartenburg M, Marti HH, Bauer C, Wenger RH, Acker H: Oxygen supply and oxygen dependent gene expression in differentiating embryonic stem cells. PNAS. 1996, 93: 2867-2872. 10.1073/pnas.93.7.2867.PubMed CentralView ArticlePubMedGoogle Scholar
- Dang SM, Kyba M, Perlingeiro R, Daley GQ, Zandstra PW: Efficiency of embryoid body formation and haematopoietic development from embryonic stem cells in different culture systems. Biotechnology and Bioengineering. 2002, 78 (4): 442-453. 10.1002/bit.10220.View ArticlePubMedGoogle Scholar
- Coucouvanis E, Martin GR: Signals for death and survival: a two step method for cavitation in the vertebrate embryo. Cell. 1995, 83: 279-287. 10.1016/0092-8674(95)90169-8.View ArticlePubMedGoogle Scholar
- Battersby A, Jones RD, Lilley KS, McFarlane RJ, Braig HR, Allen ND, Wakeman JA: Comparative proteomic analysis reveals differential expression of Hsp25 following the directed differentiation of mouse embryonic stem cells. Biochim Biophys Acta. 2007, 1773 (2): 147-56. 10.1016/j.bbamcr.2006.08.030.View ArticlePubMedGoogle Scholar
- Sollars V, Lu X, Xiao L, Wang X, Garfinkel MD, Ruden DM: Evidence for an epigenetic mechanism by which Hsp90 acts as a capacitor for morphological evolution. Nat Genet. 2003, 33 (1): 70-4. 10.1038/ng1067.View ArticlePubMedGoogle Scholar
- Meshorer E, Yellajoshula D, George E, Scambler PJ, Brown DT, Misteli T: Hyperdynamic plasticity of chromatin proteins in pluripotent embryonic stem cells. Dev Cell. 2006, 10 (1): 105-16. 10.1016/j.devcel.2005.10.017.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang WB, Boing S, Zhou XQ, Ji P, Dong Y, Yao Q, Muller-Tidow C: Identification of metastasis-associated genes in early stage non-small cell lung cancer by subtractive hybridization. Sheng Wu Hua Xue Yu Sheng Wu Wu Li Xue Bao (Shanghai). 2002, 34 (3): 273-8.Google Scholar
- Wan D, Gong Y, Qin W, Zhang P, Li J, Wei L, Zhou X, Li H, Qiu X, Zhong F, He L, Yu J, Yao G, Jiang H, Qian L, Yu Y, Shu H, Chen X, Xu H, Guo M, Pan Z, Chen Y, Ge C, Yang S, Gu J: Large-scale cDNA transfection screening for genes related to cancer development and progression. Proc Natl Acad Sci USA. 2004, 101 (44): 15724-9. 10.1073/pnas.0404089101.PubMed CentralView ArticlePubMedGoogle Scholar
- Brandenberger R, Khrebtukova I, Thies RS, Miura T, Jingli C, Puri R, Vasicek T, Lebkowski J, Rao M: MPSS profiling of human embryonic stem cells. BMC Dev Biol. 2004, 10: 4-10.Google Scholar
- Takada Y, Isono K, Shinga J, Turner JM, Kitamura H, Ohara O, Watanabe G, Singh PB, Kamijo T, Jenuwein T, Burgoyne PS, Koseki H: Mammalian Polycomb Scmh1 mediates exclusion of Polycomb complexes from the XY body in the pachytene spermatocytes. Development. 2007, 134 (3): 579-90. 10.1242/dev.02747.View ArticlePubMedGoogle Scholar
- Sabbattini P, Canzonetta C, Sjoberg M, Nikic S, Georgiou A, Kemball-Cook G, Auner HW, Dillon N: A novel role for the Aurora B kinase in epigenetic marking of silent chromatin in differentiated postmitotic cells. EMBO J. 2007, 26 (22): 4657-69. 10.1038/sj.emboj.7601875.PubMed CentralView ArticlePubMedGoogle Scholar
- Fukuda R, Zhang H, Kim JW, Shimoda L, Dang CV, Semenza GL: HIF-1 regulates cytochrome oxidase subunits to optimize efficiency of respiration in hypoxic cells. Cell. 2007, 129 (1): 111-22. 10.1016/j.cell.2007.01.047.View ArticlePubMedGoogle Scholar
- Kim SK, Suh MR, Yoon HS, Lee JB, Oh SK, Moon SY, Moon SH, Lee JY, Hwang JH, Cho WJ, Kim KS: Identification of developmental pluripotency associated 5 expression in human pluripotent stem cells. Stem Cells. 2005, 23 (4): 458-62. 10.1634/stemcells.2004-0245.View ArticlePubMedGoogle Scholar
- Amano H, Itakura K, Maruyama M, Ichisaka T, Nakagawa M, Yamanaka S: Identification and targeted disruption of the mouse gene encoding ESG1 (PH34/ECAT2/DPPA5). Dev Biol. 2006, 28;: 6-11.Google Scholar
- Yano M, Okano HJ, Okano H: Involvement of Hu and heterogeneous nuclear ribonucleoprotein K in neuronal differentiation through p21 mRNA post-transcriptional regulation. J Biol Chem. 2005, 280 (13): 12690-9. 10.1074/jbc.M411119200.View ArticlePubMedGoogle Scholar
- Dzwonek A, Mikula M, Ostrowski J: The diverse involvement of heterogeneous nuclear ribonucleoprotein K in mitochondrial response to insulin. FEBS Lett. 2006, 580 (7): 1839-45. 10.1016/j.febslet.2006.02.043.View ArticlePubMedGoogle Scholar
- Feliers D, Lee MJ, Ghosh-Choudhury G, Bomsztyk K, Kasinath BS: Heterogeneous nuclear ribonucleoprotein K contributes to angiotensin II stimulation of vascular endothelial growth factor mRNA translation. Am J Physiol Renal Physiol. 2007, 293 (2): F607-15. 10.1152/ajprenal.00497.2006.View ArticlePubMedGoogle Scholar
- Nishiyama A, Endo T, Takeda S, Imamura M: Identification and characterization of epsilon-sarcoglycans in the central nervous system. Brain Res Mol Brain Res. 2004, 125 (1–2): 1-12. 10.1016/j.molbrainres.2004.01.012.View ArticlePubMedGoogle Scholar
- Xu ZP, Wawrousek EF, Piatigorsky J: Transketolase haploinsufficiency reduces adipose tissue and female fertility in mice. Mol Cell Biol. 2002, 22 (17): 6142-7. 10.1128/MCB.22.17.6142-6147.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Oh JE, Freilinger A, Gelpi E, Pollak A, Hengstschläger M, Lubec G: Proteins involved in neuronal differentiation of neuroblastoma cell line N1E-115. Electrophoresis. 2007, 28 (12): 2009-17. 10.1002/elps.200600833.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.