Mouse embryonic stem cells (mESCs) cultured in the presence of LIF occupy a ground state with highly active pluripotency-associated transcriptional and epigenetic circuitry. However, ground state pluripotency in some inbred strain backgrounds is unstable in the absence of ERK1/2 and GSK3 inhibition. Using an unbiased genetic approach, we dissect the basis of this divergent response to extracellular cues by profiling gene expression and chromatin accessibility in 170 genetically heterogeneous mESCs. We map thousands of loci affecting chromatin accessibility and/or transcript abundance, including 10 QTL hotspots where genetic variation at a single locus coordinates the regulation of genes throughout the genome. For one hotspot, we identify a single enhancer variant $10 kb upstream of Lifr associated with chromatin accessibility and mediating a cascade of molecular events affecting pluripotency. We validate causation through reciprocal allele swaps, demonstrating the functional consequences of noncoding variation in gene regulatory networks that stabilize pluripotent states in vitro. ll

c Cells respond to environmental stimuli by fine-tuned regulation of gene expression. Here we investigated the dose-dependent modulation of gene expression at high temporal resolution in response to nutrient and stress signals in yeast. The GAL1 activity in cell populations is modulated in a well-defined range of galactose concentrations, correlating with a dynamic change of histone remodeling and RNA polymerase II (RNAPII) association. This behavior is the result of a heterogeneous induction delay caused by decreasing inducer concentrations across the population. Chromatin remodeling appears to be the basis for the dynamic GAL1 expression, because mutants with impaired histone dynamics show severely truncated dose-response profiles. In contrast, the GRE2 promoter operates like a rapid off/on switch in response to increasing osmotic stress, with almost constant expression rates and exclusively temporal regulation of histone remodeling and RNAPII occupancy. The Gal3 inducer and the Hog1 mitogen-activated protein (MAP) kinase seem to determine the different dose-response strategies at the two promoters. Accordingly, GAL1 becomes highly sensitive and dose independent if previously stimulated because of residual Gal3 levels, whereas GRE2 expression diminishes upon repeated stimulation due to acquired stress resistance. Our analysis reveals important differences in the way dynamic signals create dose-sensitive gene expression outputs.C ells continuously adapt their protein composition to changing environmental conditions. The regulation of gene expression is one of the fundamental mechanisms to adjust the global protein repertoire of the cell in order to maintain cell function and integrity in response to environmental challenges. Budding yeast is a powerful model to unravel the modes of transcriptional adaptation at the levels both of specific genes and of the whole organism (1, 2). Additionally, the basic structure of the signaling cascades responding to environmental perturbations is conserved from yeast to humans. It implies the alteration of core kinase activities, which modulate the expression of defense genes through a range of specific transcription factors. Extensive knowledge which precisely describes the molecular machinery and its global impact on gene expression in response to many types of stress has accumulated (3-7). However, the vast majority of these studies are performed with harsh environmental insults and therefore saturating stimulation. As a consequence, only very limited information or approaches are available to understand how cells adapt their gene expression programs to small or gradual changes in their environment.It is assumed that cells have acquired mechanisms that ensure a transcriptional response which is finely adjusted according to the strength of the stress or stimulation. However, the nature of the signaling molecules which confer gradual transcription outputs remains to be determined in most cases. Fine-tuning of gene expression responses can occur with different purposes, and th...

Summary Variability among pluripotent stem cell (PSC) lines is a prevailing issue that hampers not only experimental reproducibility but also large-scale applications and personalized cell-based therapy. This variability could result from epigenetic and genetic factors that influence stem cell behavior. Naive culture conditions minimize epigenetic fluctuation, potentially overcoming differences in PSC line differentiation potential. Here we derived PSCs from distinct mouse strains under naive conditions and show that lines from distinct genetic backgrounds have divergent differentiation capacity, confirming a major role for genetics in PSC phenotypic variability. This is explained in part through inconsistent activity of extra-cellular signaling, including the Wnt pathway, which is modulated by specific genetic variants. Overall, this study shows that genetic background plays a dominant role in driving phenotypic variability of PSCs.

Mouse embryonic stem cells (mESCs) cultured under controlled conditions occupy a stable ground state where pluripotency-associated transcriptional and epigenetic circuitry are highly active. However, mESCs from some genetic backgrounds exhibit metastability, where ground state pluripotency is lost in the absence of ERK1/2 and GSK3 inhibition. We dissected the genetic basis of metastability by profiling gene expression and chromatin accessibility in 185 genetically heterogeneous mESCs. We mapped thousands of loci affecting chromatin accessibility and/or transcript abundance, including eleven instances where distant QTL co-localized in clusters. For one cluster we identified Lifr transcript abundance as the causal intermediate regulating 122 distant genes enriched for roles in maintenance of pluripotency. Joint mediation analysis implicated a single enhancer variant ~10kb upstream of Lifr that alters chromatin accessibility and precipitates a cascade of molecular events affecting maintenance of pluripotency. We validated this hypothesis using reciprocal allele swaps, revealing mechanistic details underlying variability in ground state metastability in mESCs. P < 0.05) and 6,589 eQTL for 5,853 distinct genes (Fig 2C, right; LOD > 7.6, permutation-based genome-wide

No abstract

Genetic background is a major driver of phenotypic variability in pluripotent stem cells (PSCs). Most studies of variation in PSCs have relied on transcript abundance as the primary molecular readout of cell state. However, little is known about how proteins, the primary functional units in the cell, vary across genetically diverse PSCs, how protein abundance relates to variation in other cell characteristics, and how genetic background confers these effects. Here we present a comprehensive genetic study characterizing the pluripotent proteome of 190 unique mouse embryonic stem cell lines (mESCs) derived from genetically heterogeneous Diversity Outbred (DO) mice. The quantitative proteome is highly variable across DO mESCs, and we identified differentially activated pluripotency-associated pathways in the proteomics data that were not evident in transcriptome data from the same cell lines. Comparisons of protein abundance to transcript levels and chromatin accessibility show broad co-variation across molecular layers and variable correlation across samples, with some lines showing high and others low correlation between these multi-omics datasets. Integration of these three molecular data types using multi-omics factor analysis revealed shared and unique drivers of quantitative variation in pluripotency-associated pathways. QTL mapping localized the genetic drivers of this quantitative variation to a number of genomic hotspots, and demonstrated that multi-omics data integration consolidates the influence of genetic signals shared across molecular traits to increase QTL detection power and overcome the limitations inherent in mapping individual molecular features. This study reveals transcriptional and post-transcriptional mechanisms and genetic interactions that underlie quantitative variability in the pluripotent proteome, and in so doing provides a regulatory map for mouse ESCs that can provide a rational basis for future mechanistic studies, including studies of human PSCs.