Tracking the embryonic stem cell transition from ground state pluripotency

Mouse embryonic stem (ES) cells are locked into self-renewal by shielding from inductive cues. Release from this ground state in minimal conditions offers a system for delineating developmental progression from naive pluripotency. Here we examined the initial transition process. The ES cell population behaves asynchronously. We therefore exploited a short-half-life $\textit{Rex1::GFP}$ reporter to isolate cells either side of exit from naive status. Extinction of ES cell identity in single cells is acute. It occurs only after near-complete elimination of naïve pluripotency factors, but precedes appearance of lineage specification markers. Cells newly departed from the ES cell state display features of early post-implantation epiblast and are distinct from primed epiblast. They also exhibit a genome-wide increase in DNA methylation, intermediate between early and late epiblast. These findings are consistent with the proposition that naive cells transition to a distinct formative phase of pluripotency preparatory to lineage priming.This research was funded by the Wellcome Trust (091484/Z/10/Z and 095645/Z/11/Z), the Biotechnology and Biological Sciences Research Council (BB/M004023/1 and BB/K010867/1), a European Commission Framework 7 project EuroSyStem (HEALTH-F4-2007-200720 EUROSYSTEM), SysStemCell (ERC-2013-AdG 339431), the Medical Research Council (MRC) (G1100526/1) the Louis-Jeantet Foundation and the Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO-VIDI 864.12.007). The Cambridge Stem Cell Institute receives core funding from the Wellcome Trust and Medical Research Council (MRC). A.S. is an MRC Professor. Deposited in PMC for immediate release

A Genome-Wide RNAi Screen Reveals MAP Kinase Phosphatases as Key ERK Pathway Regulators during Embryonic Stem Cell Differentiation

Embryonic stem cells and induced pluripotent stem cells represent potentially important therapeutic agents in regenerative medicine. Complex interlinked transcriptional and signaling networks control the fate of these cells towards maintenance of pluripotency or differentiation. In this study we have focused on how mouse embryonic stem cells begin to differentiate and lose pluripotency and, in particular, the role that the ERK MAP kinase and GSK3 signaling pathways play in this process. Through a genome-wide siRNA screen we have identified more than 400 genes involved in loss of pluripotency and promoting the onset of differentiation. These genes were functionally associated with the ERK and/or GSK3 pathways, providing an important resource for studying the roles of these pathways in controlling escape from the pluripotent ground state. More detailed analysis identified MAP kinase phosphatases as a focal point of regulation and demonstrated an important role for these enzymes in controlling ERK activation kinetics and subsequently determining early embryonic stem cell fate decisions

ESC chimaera dynamics

The process by which pluripotent cells incorporate into host embryos is of interest to investigate cell potency and cell fate decisions. Previous studies suggest that only a minority of the embryonic stem cell (ESC) inoculum contributes to the adult chimaera. How incoming cells are chosen for integration or elimination remains unclear. By comparing a heterogeneous mix of undifferentiated and differentiating ESCs (serum/LIF) with more homogeneous undifferentiated culture (2i/LIF), we examine the role of cellular heterogeneity in this process. Time-lapse ex vivo imaging revealed a drastic elimination of serum/LIF ESCs during early development in comparison with 2i/LIF ESCs. Using a fluorescent reporter for naive pluripotency (Rex1-GFP), we established that the acutely eliminated serum/LIF ESCs had started to differentiate. The rejected cells were apparently killed by apoptosis. We conclude that a selection process exists by which unwanted differentiating cells are eliminated from the embryo. However, occasional Rex1(-) cells were able to integrate. Upregulation of Rex1 occurred in a proportion of these cells, reflecting the potential of the embryonic environment to expedite diversion from differentiation priming to enhance the developing embryonic epiblast.We thank Samuel Jameson and the Stem Cell and Gurdon Institute animal facilities for mouse husbandry, Carla Mulas for providing mKO-Rex1-GFP ESC line and the homozygous Rex1GFPd2 ESC line, Kenneth Jones for expert tissue culture assistance, Bernhard Strauss for instruction in embryo immobilisation, Charles Dumeau and Bill Mansfield for embryo transfers and animal photography. We are also indebted to Tristan Rodriguez and Austin Smith for valuable discussions, and Thorsten Boroviak and Sarra Achourri for helpful comments on the manuscript. This work was supported by the Wellcome Trust, Medical Research Council, University of Cambridge (UK), and by SFB873 funded by the Deutsche Forschungsgemeinschaft (DFG) and the Dietmar Hopp Foundation (to A.T.)This is the final version of the article. It was first available from the Company of Biologists via http://dx.doi.org/10.1242/dev.12460

A genome-wide RNAi library screen identifies factors involved in signal-dependent embryonic stem cell differentiation.

<p>(A) Schematic representation of primary screen and secondary validation screens using mouse embryonic stem cell lines containing <i>GFP</i> reporter constructs under the control of either the endogenous <i>rex1</i> (<i>Rex1</i>GFPd2) and/or <i>oct4</i> (<i>Oct4</i>GFP) promoters. “−2i” indicates that the two kinase inhibitors (CHIR99021 and PD0325901) were removed for 28 hrs (<i>Rex1</i>GFPd2 cells) and 72 hrs (<i>Oct4</i>GFP cells) before quantifying the GFP-positive population of cells. The Venn diagram shows 316 high confident hits resulting from the overlap of both validation screens. (B) The z-score of each of two biological replicates from the primary screens are plotted against each other (left panel). The average of ranked z-scores from the knockdown of individual siRNA pools is shown (right panel). The arrows indicate the z-score threshold (2 and −2). The green and red boxes mark the hits with z-score >2 and <−2, respectively. (C) Representative FACS profiles from control non-targeting and positive siRNA hits from the primary <i>rex1-GFP(d2)</i> screen. The high GFP expressing population is depicted in dark green and numbers above each graph are the corresponding GFP high/GFP low ratios. (D and E) The log₂ ratio of GFP(+)/GFP(−) of each of the two biological replicates are plotted against each other in either the <i>rex1-GFP(d2)</i> (D) or <i>oct4-GFP</i> (E) validation screens (left panels) and graphical representations of the ranked log₂ values of these ratios from the average of two independent experiments upon knockdown of individual genes are shown (right panels). The yellow shaded boxes indicate the positive hits which scored as a GFP(+)/GFP(−) ratio above 1.25× standard deviation of the controls in each of the sub-screens. (F and G) Enriched KEGG (F) and “molecular function” level GO terms (G) amongst the high confidence hits identified in both of the secondary validation screens.</p

The role of ERK pathway-specific hits in the expression of pluripotency and early differentiation marker genes.

<p>(A) Venn diagram (top) and heatmap summary (bottom; see <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003112#pgen-1003112-g002" target="_blank">Figure 2B</a> for details) illustrating the number and a list of selected screen hits used in the subsequent studies. These selected hits are distributed within three categories as indicated on the heatmap. The colour scale represents the ratio of high to low GFP expressing cells for each siRNA pool in each of the “1i” screens. (B) <i>rex1</i> (x-axis) and <i>nanog</i> (y-axis) mRNA expression levels following 2i withdrawal for 36 hrs are plotted upon knockdown of individual genes (see <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003112#pgen.1003112.s009" target="_blank">Figure S9B</a> for details). Data are shown for each siRNA duplex relative to the maximal expression exhibited in the presence of an siRNA pool (taken as 100). Dotted lines represent the expression values >2 standard deviations above the mean of the negative control siRNAs. Red dots represent siRNA duplexes which promote elevated expression of both genes (quadrant 1), whereas green (quadrant 3) and black (quadrant 4) dots represent siRNAs that cause changes at or below this threshold cut-off value for only one gene. The brown dot represents the negative control siRNAs. (C) <i>rex1</i> (y-axis) and the reciprocal of <i>fgf5</i> (x-axis) mRNA expression levels upon 2i withdrawal for 36 hrs and 48 hrs, respectively, are plotted following knockdown of individual genes (see <a href="http://www.plosgenetics.org/article/info:doi/10.1371/journal.pgen.1003112#pgen.1003112.s009" target="_blank">Figure S9B and C f</a>or details). The labeling is as indicated in (B), except that red dots represent siRNA duplexes which promote elevated expression of <i>rex1</i> and lower levels of <i>fgf5</i> (quadrant 1). Blue dots (quadrant 2) represent siRNAs that cause elevated <i>rex1</i> expression but fail to show reductions in <i>fgf5</i> expression. (D) Alkaline phosphatase staining of <i>Rex1</i>GFPd2 ES cells following treatment of cells with siRNAs against <i>gmnn</i> or <i>3830406c13rik</i> or a non-targeting control (ctrl) and release from “2i” for 5 days. Data are means ± SEM (n = 2) (E) RT-PCR analysis of the expression of the indicated lineage marker genes following treatment of cells with siRNAs against <i>gmnn</i> or <i>3830406c13rik</i> or a non-targeting control (ctrl) and release from “2i” for 3 (top) or 5 days (bottom). Data are presented as means ± SEM (n = 2).</p

Secondary screening and association of genes to the ERK and/or GSK3 pathways.

<p>(A) Schematic representation of the strategy used to stratify the high confidence positive hits via either the ERK and/or GSK3 pathways (1i-screens) using the <i>rex1-GFP</i> reporter system. A counter screen was performed in the presence of the two kinase inhibitors (CHIR99021 and PD0325901)(“+2i”) and remaining hits were tested when either the GSK3 inhibitor (CHIR99021; -GSKi) or the MEK inhibitor (PD0325901; -MEKi) was withdrawn. The Venn diagram illustrates the four different hit categories. (B) Heatmap summary depicts the stratification of the hits according to their effects on the GFP(+)/GFP(−) ratio upon withdrawal of the GSK3 inhibitor (-GSK3i) or MEK inhibitor (-MEKi). The numbers of genes in each category are indicated. The colour scale represents the ratio of high to low GFP expressing cells for each siRNA pool in each of the “1i” screens. (C and D) Heatmaps of the enriched GO terms identified for genes corresponding to hits specific to the total dataset from the “2i” screen (274), or hits from the “1i” screens; ERK (133), or GSK3 (168) pathways. Each GO term is scored by −log₁₀(P-value). The associated GO term descriptions are indicated (the GO terms enriched in only the ERK or GSK3 categories are indicated in red or green font respectively). Distinct functional groups corresponding to terms associated with cell signaling (C) and gene expression (D) are manually clustered. (E) STRING network analysis of the core network formed by the 274 genes associated with signal-dependent loss of pluripotency and promoting early differentiation processes in the mouse embryonic stem cells. Genes are grouped according to common biological processes. The coloured lines of edges represent confidence scores of interconnectivity. Dark blue lines represent 0.8–1, light blue lines represent 0.6–0.8, and light grey lines represent 0.4–0.6 confidence levels, respectively.</p

Association of siRNA screen hits with the ERK signaling pathway.

<p>(A) ERK activation levels following 2i withdrawal (−2i) for 20 mins are plotted as the ratio of phospho-ERK2 and ERK2 signals upon depletion of selected genes as indicated. The blue dashed line indicates the threshold level (1.5× SD above the mean of the negative controls) and levels below this are indicated by red bars. The average activity in the presence of control siRNA (ctrl) is shown by the solid grey line and data are plotted relative to the siRNA giving the highest levels of phospho-Erk (taken as 100). The data are presented as means ± SEM and are the average of three biological replicates (n = 3). A heatmap summary of the effect of each siRNA duplex in the “1i” screens is shown on the left. (B) Summary of the points of action of the siRNA screen hits with respect to the ERK pathway. Genes are partitioned according to which class of siRNA hits they belong. (C) Ras activity levels upon depletion of the indicated genes upon 2i withdrawal (−2i) for 2 mins. The blue dashed line indicates the threshold level (2× SD above the mean of the negative controls) and levels below this are indicated by red bars. The average activity in the presence of control (ctrl) siRNA is shown by the solid grey line and data are plotted relative to the siRNA giving the highest levels of Ras activity (taken as 100). Data are presented as means ± SEM and are the average of three biological replicates (n = 3). A heatmap summary of the effect of each siRNA duplex in the “1i” screens is shown on the left. (D) Summary diagram illustrating the point of action of upstream ERK effectors in the ERK pathway as either upstream of Ras or between Ras and Erk. The hit lists are grouped into the ERK-unique or ERK/GSK-shared hit categories. (E) Summary of the points of action of genes encoding transcriptional regulators respect to the ERK pathway.</p