ERK signalling eliminates Nanog and maintains Oct4 to drive the formative pluripotency transition.

Naïve epiblast cells in the embryo and pluripotent stem cells in vitro undergo developmental progression to a formative state competent for lineage specification. During this transition, transcription factors and chromatin are rewired to encode new functional features. Here, we examine the role of mitogen-activated protein kinase (ERK1/2) signalling in pluripotent state transition. We show that a primary consequence of ERK activation in mouse embryonic stem cells is elimination of Nanog, which precipitates breakdown of the naïve state gene regulatory network. Variability in pERK dynamics results in heterogeneous loss of Nanog and metachronous state transition. Knockdown of Nanog allows exit without ERK activation. However, transition to formative pluripotency does not proceed and cells collapse to an indeterminate identity. This outcome is due to failure to maintain expression of the central pluripotency factor Oct4. Thus, during formative transition ERK signalling both dismantles the naïve state and preserves pluripotency. These results illustrate how a single signalling pathway can both initiate and secure transition between cell states.

Individual transcription factors modulate both the micromovement of chromatin and its long-range structure.

The control of eukaryotic gene expression is intimately connected to highly dynamic chromatin structures. Gene regulation relies on activator and repressor transcription factors (TFs) that induce local chromatin opening and closing. However, it is unclear how nucleus-wide chromatin organization responds dynamically to the activity of specific TFs. Here, we examined how two TFs with opposite effects on local chromatin accessibility modulate chromatin dynamics nucleus-wide. We combine high-resolution diffusion mapping and dense flow reconstruction and correlation in living cells to obtain an imaging-based, nanometer-scale analysis of local diffusion processes and long-range coordinated movements of both chromatin and TFs. We show that the expression of either an individual transcriptional activator (CDX2) or repressor (SIX6) with large numbers of binding sites increases chromatin mobility nucleus-wide, yet they induce opposite coherent chromatin motions at the micron scale. Hi-C analysis of higher-order chromatin structures shows that induction of the pioneer factor CDX2 leads both to changes in local chromatin interactions and the distribution of A and B compartments, thus relating the micromovement of chromatin with changes in compartmental structures. Given that inhibition of transcription initiation and elongation by RNA Pol II has a partial impact on the global chromatin dynamics induced by CDX2, we suggest that CDX2 overexpression alters chromatin structure dynamics both dependently and independently of transcription. Our biophysical analysis shows that sequence-specific TFs can influence chromatin structure on multiple architectural levels, arguing that local chromatin changes brought by TFs alter long-range chromatin mobility and its organization.

An out-of-equilibrium definition of protein turnover.

Protein turnover (PT) has been formally defined only in equilibrium conditions, which is ill-suited to quantify PT during dynamic processes that occur during embryogenesis or (extra) cellular signaling. In this Hypothesis, we propose a definition of PT in an out-of-equilibrium regime that allows the quantification of PT in virtually any biological context. We propose a simple mathematical and conceptual framework applicable to a broad range of available data, such as RNA sequencing coupled with pulsed-SILAC datasets. We apply our framework to a published dataset and show that stimulation of mouse dendritic cells with LPS leads to a proteome-wide change in PT. This is the first quantification of PT out-of-equilibrium, paving the way for the analysis of biological systems in other contexts.

A role for mitotic bookmarking of SOX2 in pluripotency and differentiation.

Mitotic bookmarking transcription factors remain bound to chromosomes during mitosis and were proposed to regulate phenotypic maintenance of stem and progenitor cells at the mitosis-to-G1 (M-G1) transition. However, mitotic bookmarking remains largely unexplored in most stem cell types, and its functional relevance for cell fate decisions remains unclear. Here we screened for mitotic chromosome binding within the pluripotency network of embryonic stem (ES) cells and show that SOX2 and OCT4 remain bound to mitotic chromatin through their respective DNA-binding domains. Dynamic characterization using photobleaching-based methods and single-molecule imaging revealed quantitatively similar specific DNA interactions, but different nonspecific DNA interactions, of SOX2 and OCT4 with mitotic chromatin. Using ChIP-seq (chromatin immunoprecipitation [ChIP] combined with high-throughput sequencing) to assess the genome-wide distribution of SOX2 on mitotic chromatin, we demonstrate the bookmarking activity of SOX2 on a small set of genes. Finally, we investigated the function of SOX2 mitotic bookmarking in cell fate decisions and show that its absence at the M-G1 transition impairs pluripotency maintenance and abrogates its ability to induce neuroectodermal differentiation but does not affect reprogramming efficiency toward induced pluripotent stem cells. Our study demonstrates the mitotic bookmarking property of SOX2 and reveals its functional importance in pluripotency maintenance and ES cell differentiation.

Single-molecule dynamics and genome-wide transcriptomics reveal that NF-kB (p65)-DNA binding times can be decoupled from transcriptional activation.

Transcription factors (TFs) regulate gene expression in both prokaryotes and eukaryotes by recognizing and binding to specific DNA promoter sequences. In higher eukaryotes, it remains unclear how the duration of TF binding to DNA relates to downstream transcriptional output. Here, we address this question for the transcriptional activator NF-κB (p65), by live-cell single molecule imaging of TF-DNA binding kinetics and genome-wide quantification of p65-mediated transcription. We used mutants of p65, perturbing either the DNA binding domain (DBD) or the protein-protein transactivation domain (TAD). We found that p65-DNA binding time was predominantly determined by its DBD and directly correlated with its transcriptional output as long as the TAD is intact. Surprisingly, mutation or deletion of the TAD did not modify p65-DNA binding stability, suggesting that the p65 TAD generally contributes neither to the assembly of an "enhanceosome," nor to the active removal of p65 from putative specific binding sites. However, TAD removal did reduce p65-mediated transcriptional activation, indicating that protein-protein interactions act to translate the long-lived p65-DNA binding into productive transcription.

Modulation of transcriptional burst frequency by histone acetylation.

Many mammalian genes are transcribed during short bursts of variable frequencies and sizes that substantially contribute to cell-to-cell variability. However, which molecular mechanisms determine bursting properties remains unclear. To probe putative mechanisms, we combined temporal analysis of transcription along the circadian cycle with multiple genomic reporter integrations, using both short-lived luciferase live microscopy and single-molecule RNA-FISH. Using the Bmal1 circadian promoter as our model, we observed that rhythmic transcription resulted predominantly from variations in burst frequency, while the genomic position changed the burst size. Thus, burst frequency and size independently modulated Bmal1 transcription. We then found that promoter histone-acetylation level covaried with burst frequency, being greatest at peak expression and lowest at trough expression, while remaining unaffected by the genomic location. In addition, specific deletions of ROR-responsive elements led to constitutively elevated histone acetylation and burst frequency. We then investigated the suggested link between histone acetylation and burst frequency by dCas9p300-targeted modulation of histone acetylation, revealing that acetylation levels influence burst frequency more than burst size. The correlation between acetylation levels at the promoter and burst frequency was also observed in endogenous circadian genes and in embryonic stem cell fate genes. Thus, our data suggest that histone acetylation-mediated control of transcription burst frequency is a common mechanism to control mammalian gene expression.

Endogenous fluctuations of OCT4 and SOX2 bias pluripotent cell fate decisions.

SOX2 and OCT4 are pioneer transcription factors playing a key role in embryonic stem (ES) cell self-renewal and differentiation. How temporal fluctuations in their expression levels bias lineage commitment is unknown. Here, we generated knock-in reporter fusion ES cell lines allowing to monitor endogenous SOX2 and OCT4 protein fluctuations in living cells and to determine their impact on mesendodermal and neuroectodermal commitment. We found that small differences in SOX2 and OCT4 levels impact cell fate commitment in G1 but not in S phase. Elevated SOX2 levels modestly increased neuroectodermal commitment and decreased mesendodermal commitment upon directed differentiation. In contrast, elevated OCT4 levels strongly biased ES cells towards both neuroectodermal and mesendodermal fates in undirected differentiation. Using ATAC-seq on ES cells gated for different endogenous SOX2 and OCT4 levels, we found that high OCT4 levels increased chromatin accessibility at differentiation-associated enhancers. This suggests that small endogenous fluctuations of pioneer transcription factors can bias cell fate decisions by concentration-dependent priming of differentiation-associated enhancers.

Spatial organization of RNA polymerase II inside a mammalian cell nucleus revealed by reflected light-sheet superresolution microscopy.

Superresolution microscopy based on single-molecule centroid determination has been widely applied to cellular imaging in recent years. However, quantitative imaging of the mammalian nucleus has been challenging due to the lack of 3D optical sectioning methods for normal-sized cells, as well as the inability to accurately count the absolute copy numbers of biomolecules in highly dense structures. Here we report a reflected light-sheet superresolution microscopy method capable of imaging inside the mammalian nucleus with superior signal-to-background ratio as well as molecular counting with single-copy accuracy. Using reflected light-sheet superresolution microscopy, we probed the spatial organization of transcription by RNA polymerase II (RNAP II) molecules and quantified their global extent of clustering inside the mammalian nucleus. Spatiotemporal clustering analysis that leverages on the blinking photophysics of specific organic dyes showed that the majority (>70%) of the transcription foci originate from single RNAP II molecules, and no significant clustering between RNAP II molecules was detected within the length scale of the reported diameter of "transcription factories." Colocalization measurements of RNAP II molecules equally labeled by two spectrally distinct dyes confirmed the primarily unclustered distribution, arguing against a prevalent existence of transcription factories in the mammalian nucleus as previously proposed. The methods developed in our study pave the way for quantitative mapping and stoichiometric characterization of key biomolecular species deep inside mammalian cells.

Single-molecule imaging of transcription factor binding to DNA in live mammalian cells.

Imaging single fluorescent proteins in living mammalian cells is challenged by out-of-focus fluorescence excitation. To reduce out-of-focus fluorescence we developed reflected light-sheet microscopy (RLSM), a fluorescence microscopy method allowing selective plane illumination throughout the nuclei of living mammalian cells. A thin light sheet parallel to the imaging plane and close to the sample surface is generated by reflecting an elliptical laser beam incident from the top by 90° with a small mirror. The thin light sheet allows for an increased signal-to-background ratio superior to that in previous illumination schemes and enables imaging of single fluorescent proteins with up to 100-Hz time resolution. We demonstrated the single-molecule sensitivity of RLSM by measuring the DNA-bound fraction of glucocorticoid receptor (GR) and determining the residence times on DNA of various oligomerization states and mutants of GR and estrogen receptor-α (ER), which permitted us to resolve different modes of DNA binding of GR. We demonstrated two-color single-molecule imaging by observing the spatiotemporal colocalization of two different protein pairs. Our single-molecule measurements and statistical analysis revealed dynamic properties of transcription factors.

Stimulus-induced modulation of transcriptional bursting in a single mammalian gene.

Mammalian genes are often transcribed discontinuously as short bursts of RNA synthesis followed by longer silent periods. However, how these "on" and "off" transitions, together with the burst sizes, are modulated in single cells to increase gene expression upon stimulation is poorly characterized. By combining single-cell time-lapse luminescence imaging with stochastic modeling of the time traces, we quantified the transcriptional responses of the endogenous connective tissue growth factor gene to different physiological stimuli: serum and TGF-ß1. Both stimuli caused a rapid and acute increase in burst sizes. Whereas TGF-ß1 showed prolonged transcriptional activation mediated by an increase of transcription rate, serum stimulation resulted in a large and temporally tight first transcriptional burst, followed by a refractory period in the range of hours. Our study thus reveals how different physiological stimuli can trigger kinetically distinct transcriptional responses of the same gene.

Gene expression flux analysis reveals specific regulatory modalities of gene expression.

The level of a given protein is determined by the synthesis and degradation rates of its mRNA and protein. While several studies have quantified the contribution of different gene expression steps in regulating protein levels, these are limited by using equilibrium approximations in out-of-equilibrium biological systems. Here, we introduce gene expression flux analysis to quantitatively dissect the dynamics of the expression level for specific proteins and use it to analyze published transcriptomics and proteomics datasets. Our analysis reveals distinct regulatory modalities shared by sets of genes with clear functional signatures. We also find that protein degradation plays a stronger role than expected in the adaptation of protein levels. These findings suggest that shared regulatory strategies can lead to versatile responses at the protein level and highlight the importance of going beyond equilibrium approximations to dissect the quantitative contribution of different steps of gene expression to protein dynamics.

Isoform- and dose-sensitive feedback interactions between paired box 6 gene and delta-catenin in cell differentiation and death.

Pax6, a mammalian homolog of the Drosophila paired box gene family member expressed in stem and progenitor cells, resides at the top of the genetic hierarchy in controlling cell fates and morphogenesis. While Pax6 activation can lead to mitotic arrest, premature neurogenesis, and apoptosis, the underlying molecular mechanisms have not been resolved. Here we report that either Pax6(+5a) or Pax6(-5a) was sufficient to promote, whereas their knockdown reduced the expression of delta-catenin (CTNND2), a neural specific member of the armadillo/beta-catenin superfamily. Pax6(+5a) elicited stronger effects on delta-catenin than Pax6(-5a). Inducible Pax6(+5a) expression demonstrated a biphasic and dose-dependent regulation of delta-catenin expression and cell fates. A moderate upregulation of Pax6(+5a) promoted delta-catenin expression and induced neurite-like cellular protrusions, but increasing expression of Pax6(+5a) reversed these processes. Furthermore, sustained high expression of Pax6(+5a) triggered apoptosis as determined by the reduction of phospho-Bad, Bcl-2, survivin and procaspases, as well as the increases in Bax and cleaved poly(ADP-ribose) polymerase. Importantly, re-introducing delta-catenin by ectopic expression elicited a feedback suppression on Pax6(+5a) expression and reduced Pax6(+5a) induced apoptosis. Therefore, delta-catenin expression is not only controlled by Pax6, but it also provides a feedback suppression mechanism for their functional interactions with important implications in cellular morphogenesis, apoptosis, and cancer.

A Sox1 to Pax6 switch drives neuroectoderm to radial glia progression during differentiation of mouse embryonic stem cells.

The transcription factors Sox1 and Pax6 are expressed sequentially during early mouse embryonic neurogenesis. Sox1 expression starts upon formation of neuroectoderm, whereas Pax6 is subsequently expressed in radial glial cells, the latter giving rise to most neurons of the cerebral cortex. Here we used mouse embryonic stem (ES) cells to study the role of Sox1 and Pax6 in regulating differentiation of neural progenitors. For this purpose, we investigated the effect of overexpression and knockdown of Sox1 and Pax6, using three differentiation protocols. We show that (a) expression of Sox1 or Pax6 in uncommitted ES cells favored neuroectodermal lineage choice; (b) continuous Sox1 expression maintained cells at the neuroepithelial stage and prevented expression of Pax6 and other radial glial cell markers; (c) Sox1 knockdown facilitated exit from the progenitor stage, whereas Pax6 knockdown decreased formation of radial glia; (d) forced Pax6 expression in neuroepithelial cells triggered their differentiation into radial glia and neurons; and (e) Pax6 expression induced cell migration, a feature typical of radial glia-derived early neurons. We conclude that Sox1 enhances neuroectodermal commitment and maintenance but blocks further differentiation. In contrast, Pax6 is involved in the progression of neuroectoderm toward radial glia.

Development of human nervous tissue upon differentiation of embryonic stem cells in three-dimensional culture.

Researches on neural differentiation using embryonic stem cells (ESC) require analysis of neurogenesis in conditions mimicking physiological cellular interactions as closely as possible. In this study, we report an air-liquid interface-based culture of human ESC. This culture system allows three-dimensional cell expansion and neural differentiation in the absence of added growth factors. Over a 3-month period, a macroscopically visible, compact tissue developed. Histological coloration revealed a dense neural-like neural tissue including immature tubular structures. Electron microscopy, immunochemistry, and electrophysiological recordings demonstrated a dense network of neurons, astrocytes, and oligodendrocytes able to propagate signals. Within this tissue, tubular structures were niches of cells resembling germinal layers of human fetal brain. Indeed, the tissue contained abundant proliferating cells expressing markers of neural progenitors. Finally, the capacity to generate neural tissues on air-liquid interface differed for different ESC lines, confirming variations of their neurogenic potential. In conclusion, this study demonstrates in vitro engineering of a human neural-like tissue with an organization that bears resemblance to early developing brain. As opposed to previously described methods, this differentiation (a) allows three-dimensional organization, (b) yields dense interconnected neural tissue with structurally and functionally distinct areas, and (c) is spontaneously guided by endogenous developmental cues.

Transcription Factors and DNA Play Hide and Seek.

Transcription factors (TFs) bind to specific DNA motifs to regulate the expression of target genes. To reach their binding sites, TFs diffuse in 3D and perform local motions such as 1D sliding, hopping, or intersegmental transfer. TF-DNA interactions depend on multiple parameters, such as the chromatin environment, TF partitioning into distinct subcellular regions, and cooperativity with other DNA-binding proteins. In this review, how current understanding of the search process has initially been shaped by prokaryotic studies is discussed, as well as what is known about the parameters regulating TF search efficiency in the context of the complex eukaryotic chromatin landscape.

Phenazopyridine induces and synchronizes neuronal differentiation of embryonic stem cells.

Embryonic stem (ES) cells are powerful tools to understand mechanisms of neuronal differentiation and to engineer neurons for in vitro studies and cell therapy. We developed a screening approach to identify small organic molecules driving neuronal differentiation of ES cells. For this purpose, we used a lentivector carrying a dual luciferase reporter system to engineer an ES cell line which allowed us to screen for small organic molecules enhancing neuronal differentiation. One of them, phenazopyridine, was further analysed in human ES cells. Phenazopyridine: (i) enhanced neuronal differentiation, (ii) increased cell survival, (iii) decreased the amount of non-neuronal and undifferentiated cells and (iv) synchronized the cellular differentiation state. Phenazopyridine allowed the development of a differentiation protocol compatible with the generation of clinical grade neural precursors, which were able differentiate into different neuronal subtypes, astrocytes and oligodendrocytes. In summary, we describe a powerful approach to identify small molecules directing stem cell differentiation. This led to the establishment of a new application for an old drug and the development of a novel clinical grade protocol for neuronal differentiation of ES cells.

A pure population of ectodermal cells derived from human embryonic stem cells.

Embryonic stem (ES) cells represent a unique cellular model to recapitulate in vitro early steps of embryonic development and an unlimited cellular source in therapy for many diseases, as well as targets for drug discovery and toxicology screens. Although previous studies have reported epidermal differentiation of mouse and human embryonic stem (huES) cells, the heterogeneity of the resulting cell culture impairs the evaluation of differentiated cells for cell therapy. We report here the reproducible isolation of a homogenous ectodermal cell population, IT1, from human ES cells. Like primary cells, IT1 cells remain homogenous over 15 passages, expand up to 60 population doublings, and then die through senescence. Accordingly, IT1 cells display a normal karyotype and a somatic cell cycle kinetics and do not produce teratoma in nude mice. The production of K14-expressing epithelial cells driven by p63 expression strengthens the ectodermal nature of IT1 cells. Since IT1 can be isolated from different huES cell lines, it may provide a ready source of ectodermal progenitors for the development of a toxicology cell model, new-drug-screening strategies, and cell therapy transplantation.

Quantitative relationships between SMAD dynamics and target gene activation kinetics in single live cells.

The transduction of extracellular signals through signaling pathways that culminate in a transcriptional response is central to many biological processes. However, quantitative relationships between activities of signaling pathway components and transcriptional output of target genes remain poorly explored. Here we developed a dual bioluminescence imaging strategy allowing simultaneous monitoring of nuclear translocation of the SMAD4 and SMAD2 transcriptional activators upon TGF-ß stimulation, and the transcriptional response of the endogenous connective tissue growth factor (ctgf) gene. Using cell lines allowing to vary exogenous SMAD4/2 expression levels, we performed quantitative measurements of the temporal profiles of SMAD4/2 translocation and ctgf transcription kinetics in hundreds of individual cells at high temporal resolution. We found that while nuclear translocation efficiency had little impact on initial ctgf transcriptional activation, high total cellular SMAD4 but not SMAD2 levels increased the probability of cells to exhibit a sustained ctgf transcriptional response. The approach we present here allows time-resolved single cell quantification of transcription factor dynamics and transcriptional responses and thereby sheds light on the quantitative relationship between SMADs and target gene responses.

Memory and relatedness of transcriptional activity in mammalian cell lineages.

Phenotypically identical mammalian cells often display considerable variability in transcript levels of individual genes. How transcriptional activity propagates in cell lineages, and how this varies across genes is poorly understood. Here we combine live-cell imaging of short-lived transcriptional reporters in mouse embryonic stem cells with mathematical modelling to quantify the propagation of transcriptional activity over time and across cell generations in phenotypically homogenous cells. In sister cells we find mean transcriptional activity to be strongly correlated and transcriptional dynamics tend to be synchronous; both features control how quickly transcriptional levels in sister cells diverge in a gene-specific manner. Moreover, mean transcriptional activity is transmitted from mother to daughter cells, leading to multi-generational transcriptional memory and causing inter-family heterogeneity in gene expression.