Metastability and dynamics of stem cells: from direct observations to inference at the single cell level

Organismal development, homeostasis, and pathology are rooted in inherently probabilistic events. From gene expression to cellular differentiation, rates and likelihoods shape the form and function of biology. Processes ranging from growth to cancer homeostasis to reprogramming of stem cells all require transitions between distinct phenotypic states, and these occur at defined rates. Therefore, measuring the fidelity and dynamics with which such transitions occur is central to understanding natural biological phenomena and is critical for therapeutic interventions.

While these processes may produce robust population-level behaviors, decisions are made by individual cells. In certain circumstances, these minuscule computing units effectively roll dice to determine their fate. And while the 'omics' era has provided vast amounts of data on what these populations are doing en masse, the behaviors of the underlying units of these processes get washed out in averages.

Therefore, in order to understand the behavior of a sample of cells, it is critical to reveal how its underlying components, or mixture of cells in distinct states, each contribute to the overall phenotype. As such, we must first define what states exist in the population, determine what controls the stability of these states, and measure in high dimensionality the dynamics with which these cells transition between states.

To address a specific example of this general problem, we investigate the heterogeneity and dynamics of mouse embryonic stem cells (mESCs). While a number of reports have identified particular genes in ES cells that switch between 'high' and 'low' metastable expression states in culture, it remains unclear how levels of many of these regulators combine to form states in transcriptional space. Using a method called single molecule mRNA fluorescent in situ hybridization (smFISH), we quantitatively measure and fit distributions of core pluripotency regulators in single cells, identifying a wide range of variabilities between genes, but each explained by a simple model of bursty transcription. From this data, we also observed that strongly bimodal genes appear to be co-expressed, effectively limiting the occupancy of transcriptional space to two primary states across genes studied here. However, these states also appear punctuated by the conditional expression of the most highly variable genes, potentially defining smaller substates of pluripotency.

Having defined the transcriptional states, we next asked what might control their stability or persistence. Surprisingly, we found that DNA methylation, a mark normally associated with irreversible developmental progression, was itself differentially regulated between these two primary states. Furthermore, both acute or chronic inhibition of DNA methyltransferase activity led to reduced heterogeneity among the population, suggesting that metastability can be modulated by this strong epigenetic mark.

Finally, because understanding the dynamics of state transitions is fundamental to a variety of biological problems, we sought to develop a high-throughput method for the identification of cellular trajectories without the need for cell-line engineering. We achieved this by combining cell-lineage information gathered from time-lapse microscopy with endpoint smFISH for measurements of final expression states. Applying a simple mathematical framework to these lineage-tree associated expression states enables the inference of dynamic transitions. We apply our novel approach in order to infer temporal sequences of events, quantitative switching rates, and network topology among a set of ESC states.

Taken together, we identify distinct expression states in ES cells, gain fundamental insight into how a strong epigenetic modifier enforces the stability of these states, and develop and apply a new method for the identification of cellular trajectories using scalable in situ readouts of cellular state.

Electron microprobe investigations of metamorphic reactions and mineral growth histories, Kwoiek area, British Columbia

The Kwoiek Area of British Columbia contains a pendant or screen of metamorphosed sedimentary and volcanic rocks almost entirely surrounded by a portion of the Coast Range Batholith, and intruded by several dozen stocks. The major metamorphic effects were produced by the quartz diorite batholithic rocks, with minor and later effects by the quartz diorite stocks. The sequence of important metamorphic reactions in the metasedimentary and metavolcanic rocks, ranging in grade from chlorite to sillimanite, is:

1. chlorite + carbonate + muscovite → epidote + biotite

2. chlorite + carbonate → actinolite + epidote

3. chlorite + muscovite → garnet + biotite

4. chlorite + epidote → garnet + hornblende

5. chlorite + muscovite → garnet + staurolite + biotite

6. chlorite + muscovite → aluminum silicate + biotite

7. muscovite + staurolite → garnet + aluminum silicate + biotite

8. staurolite → garnet + aluminum silicate

Continuous reactions, occurring between reactions 5 and 7, are:

A. chlorite + (high Ti) biotite + Al₂O₃ (from plagioclase?)→ garnet + staurolite + (low Ti) biotite + O₂

B. muscovite (phengitic) → garnet + staurolite +muscovite (less phengitic) + O₂ (?)

Detailed electron microprobe work on garnet, staurolite, biotite, and chlorite shows that:

(1) The garnet porphyroblasts are zoned according to a depletion model, called the Rayleigh depletion model, which assumes equilibrium between the edge of a growing garnet and the minerals which are unzoned, notably biotite, chlorite, and muscovite, but which assumes disequilibrium within the garnet.

(2) The staurolite porphyroblasts are also zoned, and from their zoning patterns reactions A, B, and 5 are documented. Progressive reduction of iron with increasing grade of metamorphism is also inferred from the staurolite zoning patterns.

(3) During a late period of falling temperature garnet continued to grow and the biotite and chlorite reequilibrated. The biotite, chlorite, and garnet edge compositions can vary from point to point in a given thin section, indicating that the volume of equilibrium at the final stage of metamorphism was only a few cubic microns.

(4) The horizon within the garnet that grew at maximum temperature can be identified. The Mg/Fe ratio of this horizon, if the garnet composition is a limiting composition in the Al₂O₃ - K₂O - FeO - MgO tetrahedron, increases systematically with increasing metamorphic grade. Biotite and chlorite compositions also show a general increase in Mg/Fe ratio with increasing metamorphic grade, but staurolite appears to show the reverse effect.

(5) The Mg/Fe ratio at the maximum temperature horizon of the garnet porphyroblasts is a function of its Mn content as evidenced from the study of five garnet-bearing rocks, collected from one outcrop area, with the same assemblage but with differing proportions of minerals.

An important implication of zoned minerals is that the effective composition of a system in a phase lies on the join between the homogeneous minerals (if there are two) and not within three-or- four-phase fields when a zoned mineral, such as garnet or staurolite, is present in the assemblage.

Study of the three aluminum silicates found in the Kwoiek Area showed that a constant pressure change in polymorphs from andalusite to kyanite to sillimanite took place with increasing temperature. This transition series is best explained by the metastable formation of andalusite.

Photographic materials on pages 15, 121, 160, 162, and 164 are essential and will not reproduce clearly on Xerox copies. Photographic copies should be ordered.