Maria Herberg

Nanog Variability and Pluripotency Regulation of Embryonic Stem Cells - Insights from a Mathematical Model Analysis

The expression of the transcription factors Oct4, Sox2, and Nanog is commonly associated with pluripotency of mouse embryonic stem (ES) cells. However, recent observations suggest that ES cell populations are heterogeneous with respect to the expression of Nanog and that individual ES cells reversibly change their Nanog expression level. Furthermore, it has been shown that cells exhibiting a low Nanog level are more likely to undergo differentiation. Applying a novel mathematical transcription factor network model we explore mechanisms and feedback regulations to describe the observed variation of the Nanog levels in mouse ES cells. In particular we show that these variations can occur under different assumptions yielding similar experimental characteristics. Based on model predictions we propose experimental strategies to distinguish between these explanations. Concluding from our results we argue that the heterogeneity with respect to the Nanog concentrations is most likely a functional element to control the differentiation propensity of an ES cell population. Furthermore, we provide a conceptual framework that consistently explains Nanog variability and a potential “gate-keeper” function of Nanog expression with respect to the control of ES cell differentiation.

A Model-Based Analysis of Culture-Dependent Phenotypes of mESCs

Mouse embryonic stem cells (mESCs) can be maintained in a proliferative and undifferentiated state over many passages (self-renewal) while retaining the potential to give rise to every cell type of the organism (pluripotency). Autocrine FGF4/Erk signalling has been identified as a major stimulus for fate decisions and lineage commitment in these cells. Recent findings on serum-free culture conditions with specific inhibitors (known as 2i) demonstrate that the inhibition of this pathway reduces transcription factor heterogeneity and is vital to maintain ground state pluripotency of mESCs. We suggest a novel mathematical model to explicitly integrate FGF4/Erk signalling into an interaction network of key pluripotency factors (namely Oct4, Sox2, Nanog and Rex1). The envisaged model allows to explore whether and how proposed mechanisms and feedback regulations can account for different expression patterns in mESC cultures. We demonstrate that an FGF4/Erk-mediated negative feedback is sufficient to induce molecular heterogeneity with respect to Nanog and Rex1 expression and thus critically regulates the propensity for differentiation and the loss of pluripotency. Furthermore, we compare simulation results on the transcription factor dynamics in different self-renewing states and during differentiation with experimental data on a Rex1GFPd2 reporter cell line using flow cytometry and qRT-PCR measurements. Concluding from our results we argue that interaction between FGF4/Erk signalling and Nanog expression qualifies as a key mechanism to manipulate mESC pluripotency. In particular, we infer that ground state pluripotency under 2i is achieved by shifting stable expression pattern of Nanog from a bistable into a monostable regulation impeding stochastic state transitions. Furthermore, we derive testable predictions on altering the degree of Nanog heterogeneity and on the frequency of state transitions in LIF/serum conditions to challenge our model assumptions.

Imaging, quantification and visualization of spatio-temporal patterning in mESC colonies under different culture conditions

Motivation: Mouse embryonic stem cells (mESCs) have developed into a prime system to study the regulation of pluripotency in stable cell lines. It is well recognized that different, established protocols for the maintenance of mESC pluripotency support morphologically and functionally different cell cultures. However, it is unclear how characteristic properties of cell colonies develop over time and how they are re-established after cell passage depending on the culture conditions. Furthermore, it appears that cell colonies have an internal structure with respect to cell size, marker expression or biomechanical properties, which is not sufficiently understood. The analysis of these phenotypic properties is essential for a comprehensive understanding of mESC development and ultimately requires a bioinformatics approach to guarantee reproducibility and high-throughput data analysis. Results: We developed an automated image analysis and colony tracking framework to obtain an objective and reproducible quantification of structural properties of cell colonies as they evolve in space and time. In particular, we established a method that quantifies changes in colony shape and (internal) motion using fluid image registration and image segmentation. The methodology also allows to robustly track motion, splitting and merging of colonies over a sequence of images. Our results provide a first quantitative assessment of temporal mESC colony formation and estimates of structural differences between colony growth under different culture conditions. Furthermore, we provide a stream-based visualization of structural features of individual colonies over time for the whole experiment, facilitating visual comprehension of differences between experimental conditions. Thus, the presented method establishes the basis for the model-based analysis of mESC colony development. It can be easily extended to integrate further functional information using fluorescence signals and differentiation markers. Availability: The analysis tool is implemented C++ and Mathematica 8.0 (Wolfram Research Inc., Champaign, IL, USA). The tool is freely available from the authors. We will also provide the source code upon request. Contact: nico.scherf@tu-dresden.de

Computational modelling of embryonic stem-cell fate control

The maintenance of pluripotency in embryonic stem cells (ESCs), its loss during lineage specification or its re-induction to generate induced pluripotent stem cells are central topics in stem cell biology. To uncover the molecular basis and the design principles of pluripotency control, a multitude of experimental, but also an increasing number of computational, studies have been published. Here, we consider recent reports that apply computational or mathematical modelling approaches to describe the regulatory processes that underlie cell fate decisions in mouse ESCs. We summarise the principles, the strengths and potentials but also the limitations of different computational strategies. Summary: This Review discusses the principles, applications and limitations of different computational strategies used to study the mechanisms underlying pluripotency regulation in embryonic stem cells.

Genomic and transcriptomic heterogeneity of colorectal tumors arising in Lynch Syndrome

Colorectal cancer (CRC) arising in Lynch syndrome (LS) comprises tumours with constitutional mutations in DNA mismatch repair genes. There is still a lack of whole‐genome and transcriptome studies of LS‐CRC to address questions about similarities and differences in mutation and gene expression characteristics between LS‐CRC and sporadic CRC, about the molecular heterogeneity of LS‐CRC, and about specific mechanisms of LS‐CRC genesis linked to dysfunctional mismatch repair in LS colonic mucosa and the possible role of immune editing. Here, we provide a first molecular characterization of LS tumours and of matched tumour‐distant reference colonic mucosa based on whole‐genome DNA‐sequencing and RNA‐sequencing analyses. Our data support two subgroups of LS‐CRCs, G1 and G2, whereby G1 tumours show a higher number of somatic mutations, a higher amount of microsatellite slippage, and a different mutation spectrum. The gene expression phenotypes support this difference. Reference mucosa of G1 shows a strong immune response associated with the expression of HLA and immune checkpoint genes and the invasion of CD4+ T cells. Such an immune response is not observed in LS tumours, G2 reference and normal (non‐Lynch) mucosa, and sporadic CRC. We hypothesize that G1 tumours are edited for escape from a highly immunogenic microenvironment via loss of HLA presentation and T‐cell exhaustion. In contrast, G2 tumours seem to develop in a less immunogenic microenvironment where tumour‐promoting inflammation parallels tumourigenesis. Larger studies on non‐neoplastic mucosa tissue of mutation carriers are required to better understand the early phases of emerging tumours. Copyright

Epigenetic Nanog regulation and the role of functional heterogeneity

Comment on: Villasante A, et al. Cell Cycle 2011; 10:1488-98.

Linking DNA Damage and Age-Related Promoter DNA Hyper-Methylation in the Intestine

Aberrant DNA methylation in stem cells is a hallmark of aging and tumor development. Here, we explore whether and how DNA damage repair might impact on these time-dependent changes, in particular in proliferative intestinal stem cells. We introduce a 3D multiscale computer model of intestinal crypts enabling simulation of aberrant DNA and histone methylation of gene promoters during aging. We assume histone state-dependent activity of de novo DNA methyltransferases (DNMTs) and methylation-dependent binding of maintenance DNMTs to CpGs. We simulate aging with and without repeated DNA repair. Motivated by recent findings on the histone demethylase KDM2b, we consider that DNA repair is associated with chromatin opening and improved recruitment of de novo DNMTs. Our results suggest that methylation-dependent binding of maintenance DNMTs to CpGs, establishing bistable DNA methylation states, is a prerequisite to promoter hyper-methylation following DNA repair. With this, the transient increase in de novo DNMT activity during repair can induce switches from low to high methylation states. These states remain stable after repair, leading to an epigenetic drift. The switches are most frequent in genes with H3K27me3 modified promoters. Our model provides a mechanistic explanation on how even successful DNA repair might confer long term changes of the epigenome.

The Regulatory Capacity of Bivalent Genes—A Theoretical Approach

Bivalent genes are frequently associated with developmental and lineage specification processes. Resolving their bivalency enables fast changes in their expression, which potentially can trigger cell fate decisions. Here, we provide a theoretical model of bivalency that allows for predictions on the occurrence, stability and regulatory capacity of this prominent modification state. We suggest that bivalency enables balanced gene expression heterogeneity that constitutes a prerequisite of robust lineage priming in somatic stem cells. Moreover, we demonstrate that interactions between the histone and DNA methylation machineries together with the proliferation activity control the stability of the bivalent state and can turn it into an unmodified state. We suggest that deregulation of these interactions underlies cell transformation processes as associated with acute myeloid leukemia (AML) and provide a model of AML blast formation following deregulation of the Ten-eleven Translocation (TET) pathway.

Linking stem cell function and growth pattern of intestinal organoids

Intestinal stem cells (ISCs) require well-defined signals from their environment in order to carry out their specific functions. Most of these signals are provided by neighboring cells that form a stem cell niche, whose shape and cellular composition self-organize. Major features of this self-organization can be studied in ISC-derived organoid culture. In this system, manipulation of essential pathways of stem cell maintenance and differentiation results in well-described growth phenotypes. We here provide an individual cell-based model of intestinal organoids that enables a mechanistic explanation of the observed growth phenotypes. In simulation studies of the 3D structure of expanding organoids, we investigate interdependences between Wnt- and Notch-signaling which control the shape of the stem cell niche and, thus, the growth pattern of the organoids. Similar to in vitro experiments, changes of pathway activities alter the cellular composition of the organoids and, thereby, affect their shape. Exogenous Wnt enforces transitions from branched into a cyst-like growth pattern; known to occur spontaneously during long term organoid expansion. Based on our simulation results, we predict that the cyst-like pattern is associated with biomechanical changes of the cells which assign them a growth advantage. The results suggest ongoing stem cell adaptation to in vitro conditions during long term expansion by stabilizing Wnt-activity. Our study exemplifies the potential of individual cell-based modeling in unraveling links between molecular stem cell regulation and 3D growth of tissues. This kind of modeling combines experimental results in the fields of stem cell biology and cell biomechanics constituting a prerequisite for a better understanding of tissue regeneration as well as developmental processes.