Introducing Point Mutations into Human Pluripotent Stem Cells using Seamless Genome Editing.

Custom designed endonucleases, such as RNA-guided Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-Cas9, enable efficient genome editing in mammalian cells. Here we describe detailed procedures to seamlessly genome edit the hepatocyte nuclear factor 4 alpha (HNF4α) locus as an example in human pluripotent stem cells. Combining a piggyBac-based donor plasmid and the CRISPR-Cas9 nickase mutant in a two-step genetic selection, we demonstrate correct and efficient targeting of the HNF4α locus.

Multiomics Analyses of HNF4α Protein Domain Function during Human Pluripotent Stem Cell Differentiation.

During mammalian development, liver differentiation is driven by signals that converge on multiple transcription factor networks. The hepatocyte nuclear factor signaling network is known to be essential for hepatocyte specification and maintenance. In this study, we have generated deletion and point mutants of hepatocyte nuclear factor-4alpha (HNF4α) to precisely evaluate the function of protein domains during hepatocyte specification from human pluripotent stem cells. We demonstrate that nuclear HNF4α is essential for hepatic progenitor specification, and the introduction of point mutations in HNF4α's Small Ubiquitin-like Modifier (SUMO) consensus motif leads to disrupted hepatocyte differentiation. Taking a multiomics approach, we identified key deficiencies in cell biology, which included dysfunctional metabolism, substrate adhesion, tricarboxylic acid cycle flux, microRNA transport, and mRNA processing. In summary, the combination of genome editing and multiomics analyses has provided valuable insight into the diverse functions of HNF4α during pluripotent stem cell entry into the hepatic lineage and during hepatocellular differentiation.

Identification and Single-Cell Functional Characterization of an Endodermally Biased Pluripotent Substate in Human Embryonic Stem Cells.

Human embryonic stem cells (hESCs) display substantial heterogeneity in gene expression, implying the existence of discrete substates within the stem cell compartment. To determine whether these substates impact fate decisions of hESCs we used a GFP reporter line to investigate the properties of fractions of putative undifferentiated cells defined by their differential expression of the endoderm transcription factor, GATA6, together with the hESC surface marker, SSEA3. By single-cell cloning, we confirmed that substates characterized by expression of GATA6 and SSEA3 include pluripotent stem cells capable of long-term self-renewal. When clonal stem cell colonies were formed from GATA6-positive and GATA6-negative cells, more of those derived from GATA6-positive cells contained spontaneously differentiated endoderm cells than similar colonies derived from the GATA6-negative cells. We characterized these discrete cellular states using single-cell transcriptomic analysis, identifying a potential role for SOX17 in the establishment of the endoderm-biased stem cell state.

A Human IPS Model Implicates Embryonic B-Myeloid Fate Restriction as Developmental Susceptibility to B Acute Lymphoblastic Leukemia-Associated ETV6-RUNX1.

ETV6-RUNX1 is associated with childhood acute B-lymphoblastic leukemia (cALL) functioning as a first-hit mutation that initiates a clinically silent pre-leukemia in utero. Because lineage commitment hierarchies differ between embryo and adult, and the impact of oncogenes is cell-context dependent, we hypothesized that the childhood affiliation of ETV6-RUNX1 cALL reflects its origins in a progenitor unique to embryonic life. We characterize the first emerging B cells in first-trimester human embryos, identifying a developmentally restricted CD19-IL-7R+ progenitor compartment, which transitions from a myeloid to lymphoid program during ontogeny. This developmental series is recapitulated in differentiating human pluripotent stem cells (hPSCs), thereby providing a model for the initiation of cALL. Genome-engineered hPSCs expressing ETV6-RUNX1 from the endogenous ETV6 locus show expansion of the CD19-IL-7R+ compartment, show a partial block in B lineage commitment, and produce proB cells with aberrant myeloid gene expression signatures and potential: features (collectively) consistent with a pre-leukemic state.

DNA methylation of intragenic CpG islands depends on their transcriptional activity during differentiation and disease.

The human genome contains ∼30,000 CpG islands (CGIs). While CGIs associated with promoters nearly always remain unmethylated, many of the ∼9,000 CGIs lying within gene bodies become methylated during development and differentiation. Both promoter and intragenic CGIs may also become abnormally methylated as a result of genome rearrangements and in malignancy. The epigenetic mechanisms by which some CGIs become methylated but others, in the same cell, remain unmethylated in these situations are poorly understood. Analyzing specific loci and using a genome-wide analysis, we show that transcription running across CGIs, associated with specific chromatin modifications, is required for DNA methyltransferase 3B (DNMT3B)-mediated DNA methylation of many naturally occurring intragenic CGIs. Importantly, we also show that a subgroup of intragenic CGIs is not sensitive to this process of transcription-mediated methylation and that this correlates with their individual intrinsic capacity to initiate transcription in vivo. We propose a general model of how transcription could act as a primary determinant of the patterns of CGI methylation in normal development and differentiation, and in human disease.

Genetic dissection of the α-globin super-enhancer in vivo.

Many genes determining cell identity are regulated by clusters of Mediator-bound enhancer elements collectively referred to as super-enhancers. These super-enhancers have been proposed to manifest higher-order properties important in development and disease. Here we report a comprehensive functional dissection of one of the strongest putative super-enhancers in erythroid cells. By generating a series of mouse models, deleting each of the five regulatory elements of the α-globin super-enhancer individually and in informative combinations, we demonstrate that each constituent enhancer seems to act independently and in an additive fashion with respect to hematological phenotype, gene expression, chromatin structure and chromosome conformation, without clear evidence of synergistic or higher-order effects. Our study highlights the importance of functional genetic analyses for the identification of new concepts in transcriptional regulation.

Intrinsic factors and the embryonic environment influence the formation of extragonadal teratomas during gestation.

BACKGROUND: Pluripotent cells are present in early embryos until the levels of the pluripotency regulator Oct4 drop at the beginning of somitogenesis. Elevating Oct4 levels in explanted post-pluripotent cells in vitro restores their pluripotency. Cultured pluripotent cells can participate in normal development when introduced into host embryos up to the end of gastrulation. In contrast, pluripotent cells efficiently seed malignant teratocarcinomas in adult animals. In humans, extragonadal teratomas and teratocarcinomas are most frequently found in the sacrococcygeal region of neonates, suggesting that these tumours originate from cells in the posterior of the embryo that either reactivate or fail to switch off their pluripotent status. However, experimental models for the persistence or reactivation of pluripotency during embryonic development are lacking. METHODS: We manually injected embryonic stem cells into conceptuses at E9.5 to test whether the presence of pluripotent cells at this stage correlates with teratocarcinoma formation. We then examined the effects of reactivating embryonic Oct4 expression ubiquitously or in combination with Nanog within the primitive streak (PS)/tail bud (TB) using a transgenic mouse line and embryo chimeras carrying a PS/TB-specific heterologous gene expression cassette respectively. RESULTS: Here, we show that pluripotent cells seed teratomas in post-gastrulation embryos. However, at these stages, induced ubiquitous expression of Oct4 does not lead to restoration of pluripotency (indicated by Nanog expression) and tumour formation in utero, but instead causes a severe phenotype in the extending anteroposterior axis. Use of a more restricted T(Bra) promoter transgenic system enabling inducible ectopic expression of Oct4 and Nanog specifically in the posteriorly-located primitive streak (PS) and tail bud (TB) led to similar axial malformations to those induced by Oct4 alone. These cells underwent induction of pluripotency marker expression in Epiblast Stem Cell (EpiSC) explants derived from somitogenesis-stage embryos, but no teratocarcinoma formation was observed in vivo. CONCLUSIONS: Our findings show that although pluripotent cells with teratocarcinogenic potential can be produced in vitro by the overexpression of pluripotency regulators in explanted somitogenesis-stage somatic cells, the in vivo induction of these genes does not yield tumours. This suggests a restrictive regulatory role of the embryonic microenvironment in the induction of pluripotency.

Hes1 desynchronizes differentiation of pluripotent cells by modulating STAT3 activity.

Robust development of the early embryo may benefit from mechanisms that ensure that not all pluripotent cells differentiate at exactly the same time: such mechanisms would build flexibility into the process of lineage allocation. This idea is supported by the observation that pluripotent stem cells differentiate at different rates in vitro. We use a clonal commitment assay to confirm that pluripotent cells commit to differentiate asynchronously even under uniform differentiation conditions. Stochastic variability in expression of the Notch target gene Hes1 has previously been reported to influence neural versus mesodermal differentiation through modulation of Notch activity. Here we report that Hes1 also has an earlier role to delay exit from the pluripotent state into all lineages. The early function of Hes1 to delay differentiation can be explained by an ability of Hes1 to amplify STAT3 responsiveness in a cell-autonomous manner. Variability in Hes1 expression therefore helps to explain why STAT3 responsiveness varies between individual ES cells, and this in turn helps to explain why pluripotent cells commit to differentiate asynchronously.

Tcf15 primes pluripotent cells for differentiation.

The events that prime pluripotent cells for differentiation are not well understood. Inhibitor of DNA binding/differentiation (Id) proteins, which are inhibitors of basic helix-loop-helix (bHLH) transcription factor activity, contribute to pluripotency by blocking sequential transitions toward differentiation. Using yeast-two-hybrid screens, we have identified Id-regulated transcription factors that are expressed in embryonic stem cells (ESCs). One of these, Tcf15, is also expressed in the embryonic day 4.5 embryo and is specifically associated with a novel subpopulation of primed ESCs. An Id-resistant form of Tcf15 rapidly downregulates Nanog and accelerates somatic lineage commitment. We propose that because Tcf15 can be held in an inactive state through Id activity, it may prime pluripotent cells for entry to somatic lineages upon downregulation of Id. We also find that Tcf15 expression is dependent on fibroblast growth factor (FGF) signaling, providing an explanation for how FGF can prime for differentiation without driving cells out of the pluripotent state.

Transposon-mediated BAC transgenesis in human ES cells.

Transgenesis is a cornerstone of molecular biology. The ability to integrate a specifically engineered piece of DNA into the genome of a living system is fundamental to our efforts to understand life and exploit its implications for medicine, nanotechnology and bioprospecting. However, transgenesis has been hampered by position effects and multi-copy integration problems, which are mainly due to the use of small, plasmid-based transgenes. Large transgenes based on native genomic regions cloned into bacterial artificial chromosomes (BACs) circumvent these problems but are prone to fragmentation. Herein, we report that contrary to widely held notions, large BAC-sized constructs do not prohibit transposition. We also report the first reliable method for BAC transgenesis in human embryonic stem cells (hESCs). The PiggyBac or Sleeping Beauty transposon inverted repeats were integrated into BAC vectors by recombineering, followed by co-lipofection with the corresponding transposase in hESCs to generate robust fluorescent protein reporter lines for OCT4, NANOG, GATA4 and PAX6. BAC transposition delivers several advantages, including increased frequencies of single-copy, full-length integration, which will be useful in all transgenic systems but especially in difficult venues like hESCs.

An interspecies analysis reveals a key role for unmethylated CpG dinucleotides in vertebrate Polycomb complex recruitment.

The role of DNA sequence in determining chromatin state is incompletely understood. We have previously demonstrated that large chromosomal segments from human cells recapitulate their native chromatin state in mouse cells, but the relative contribution of local sequences versus their genomic context remains unknown. In this study, we compare orthologous chromosomal regions for which the human locus establishes prominent sites of Polycomb complex recruitment in pluripotent stem cells, whereas the corresponding mouse locus does not. Using recombination-mediated cassette exchange at the mouse locus, we establish the primacy of local sequences in the encoding of chromatin state. We show that the signal for chromatin bivalency is redundantly encoded across a bivalent domain and that this reflects competition between Polycomb complex recruitment and transcriptional activation. Furthermore, our results suggest that a high density of unmethylated CpG dinucleotides is sufficient for vertebrate Polycomb recruitment. This model is supported by analysis of DNA methyltransferase-deficient embryonic stem cells.

Genomically humanized mice: technologies and promises.

Mouse models have become an invaluable tool for understanding human health and disease owing to our ability to manipulate the mouse genome exquisitely. Recent progress in genomic analysis has led to an increase in the number and type of disease-causing mutations detected and has also highlighted the importance of non-coding regions. As a result, there is increasing interest in creating 'genomically' humanized mouse models, in which entire human genomic loci are transferred into the mouse genome. The technical challenges towards achieving this aim are large but are starting to be tackled with success.

Foxn1 regulates lineage progression in cortical and medullary thymic epithelial cells but is dispensable for medullary sublineage divergence.

The forkhead transcription factor Foxn1 is indispensable for thymus development, but the mechanisms by which it mediates thymic epithelial cell (TEC) development are poorly understood. To examine the cellular and molecular basis of Foxn1 function, we generated a novel and revertible hypomorphic allele of Foxn1. By varying levels of its expression, we identified a number of features of the Foxn1 system. Here we show that Foxn1 is a powerful regulator of TEC differentiation that is required at multiple intermediate stages of TE lineage development in the fetal and adult thymus. We find no evidence for a role for Foxn1 in TEC fate-choice. Rather, we show it is required for stable entry into both the cortical and medullary TEC differentiation programmes and subsequently is needed at increasing dosage for progression through successive differentiation states in both cortical and medullary TEC. We further demonstrate regulation by Foxn1 of a suite of genes with diverse roles in thymus development and/or function, suggesting it acts as a master regulator of the core thymic epithelial programme rather than regulating a particular aspect of TEC biology. Overall, our data establish a genetics-based model of cellular hierarchies in the TE lineage and provide mechanistic insight relating titration of a single transcription factor to control of lineage progression. Our novel revertible hypomorph system may be similarly applied to analyzing other regulators of development.

Polycomb eviction as a new distant enhancer function.

Remote distal enhancers may be located tens or thousands of kilobases away from their promoters. How they control gene expression is still poorly understood. Here, we analyze the influence of a remote enhancer on the balance between repression (Polycomb-PcG) and activation (Trithorax-TrxG) of a developmentally regulated gene associated with a CpG island. We reveal its essential, nonredundant role in clearing the PcG complex and H3K27me3 from the CpG island. In the absence of the enhancer, the H3K27me3 demethylase (JMJD3) is not recruited to the CpG island. We propose a new role of long-range regulatory elements in removing repressive PcG complexes.

Dominant role of the p110beta isoform of PI3K over p110alpha in energy homeostasis regulation by POMC and AgRP neurons.

PI3K signaling is thought to mediate leptin and insulin action in hypothalamic pro-opiomelanocortin (POMC) and agouti-related protein (AgRP) neurons, key regulators of energy homeostasis, through largely unknown mechanisms. We inactivated either p110alpha or p110beta PI3K catalytic subunits in these neurons and demonstrate a dominant role for the latter in energy homeostasis regulation. In POMC neurons, p110beta inactivation prevented insulin- and leptin-stimulated electrophysiological responses. POMCp110beta null mice exhibited central leptin resistance, increased adiposity, and diet-induced obesity. In contrast, the response to leptin was not blocked in p110alpha-deficient POMC neurons. Accordingly, POMCp110alpha null mice displayed minimal energy homeostasis abnormalities. Similarly, in AgRP neurons, p110beta had a more important role than p110alpha. AgRPp110alpha null mice displayed normal energy homeostasis regulation, whereas AgRPp110beta null mice were lean, with increased leptin sensitivity and resistance to diet-induced obesity. These results demonstrate distinct metabolic roles for the p110alpha and p110beta isoforms of PI3K in hypothalamic energy regulation.

Chromosome looping at the human alpha-globin locus is mediated via the major upstream regulatory element (HS -40).

Previous studies in the mouse have shown that high levels of alpha-globin gene expression in late erythropoiesis depend on long-range, physical interactions between remote upstream regulatory elements and the globin promoters. Using quantitative chromosome conformation capture (q3C), we have now analyzed all interactions between 4 such elements lying 10 to 50 kb upstream of the human alpha cluster and their interactions with the alpha-globin promoter. All of these elements interact with the alpha-globin gene in an erythroid-specific manner. These results were confirmed in a mouse model of human alpha globin expression in which the human cluster replaces the mouse cluster in situ (humanized mouse). We have also shown that expression and all of the long-range interactions depend largely on just one of these elements; removal of the previously characterized major regulatory element (called HS -40) results in loss of all the interactions and alpha-globin expression. Reinsertion of this element at an ectopic location restores both expression and the intralocus interactions. In contrast to other more complex systems involving multiple upstream elements and promoters, analysis of the human alpha-globin cluster during erythropoiesis provides a simple and tractable model to understand the mechanisms underlying long-range gene regulation.

Association between active genes occurs at nuclear speckles and is modulated by chromatin environment.

Genes on different chromosomes can be spatially associated in the nucleus in several transcriptional and regulatory situations; however, the functional significance of such associations remains unclear. Using human erythropoiesis as a model, we show that five cotranscribed genes, which are found on four different chromosomes, associate with each other at significant but variable frequencies. Those genes most frequently in association lie in decondensed stretches of chromatin. By replacing the mouse alpha-globin gene cluster in situ with its human counterpart, we demonstrate a direct effect of the regional chromatin environment on the frequency of association, whereas nascent transcription from the human alpha-globin gene appears unaffected. We see no evidence that cotranscribed erythroid genes associate at shared transcription foci, but we do see stochastic clustering of active genes around common nuclear SC35-enriched speckles (hence the apparent nonrandom association between genes). Thus, association between active genes may result from their location on decondensed chromatin that enables clustering around common nuclear speckles.

The p110beta isoform of phosphoinositide 3-kinase signals downstream of G protein-coupled receptors and is functionally redundant with p110gamma.

The p110 isoforms of phosphoinositide 3-kinase (PI3K) are acutely regulated by extracellular stimuli. The class IA PI3K catalytic subunits (p110alpha, p110beta, and p110delta) occur in complex with a Src homology 2 (SH2) domain-containing p85 regulatory subunit, which has been shown to link p110alpha and p110delta to Tyr kinase signaling pathways. The p84/p101 regulatory subunits of the p110gamma class IB PI3K lack SH2 domains and instead couple p110gamma to G protein-coupled receptors (GPCRs). Here, we show, using small-molecule inhibitors with selectivity for p110beta and cells derived from a p110beta-deficient mouse line, that p110beta is not a major effector of Tyr kinase signaling but couples to GPCRs. In macrophages, both p110beta and p110gamma contributed to Akt activation induced by the GPCR agonist complement 5a, but not by the Tyr kinase ligand colony-stimulating factor-1. In fibroblasts, which express p110beta but not p110gamma, p110beta mediated Akt activation by the GPCR ligands stromal cell-derived factor, sphingosine-1-phosphate, and lysophosphatidic acid but not by the Tyr kinase ligands PDGF, insulin, and insulin-like growth factor 1. Introduction of p110gamma in these cells reduced the contribution of p110beta to GPCR signaling. Taken together, these data show that p110beta and p110gamma can couple redundantly to the same GPCR agonists. p110beta, which shows a much broader tissue distribution than the leukocyte-restricted p110gamma, could thus provide a conduit for GPCR-linked PI3K signaling in the many cell types where p110gamma expression is low or absent.

Angiogenesis selectively requires the p110alpha isoform of PI3K to control endothelial cell migration.

Phosphoinositide 3-kinases (PI3Ks) signal downstream of multiple cell-surface receptor types. Class IA PI3K isoforms couple to tyrosine kinases and consist of a p110 catalytic subunit (p110alpha, p110beta or p110delta), constitutively bound to one of five distinct p85 regulatory subunits. PI3Ks have been implicated in angiogenesis, but little is known about potential selectivity among the PI3K isoforms and their mechanism of action in endothelial cells during angiogenesis in vivo. Here we show that only p110alpha activity is essential for vascular development. Ubiquitous or endothelial cell-specific inactivation of p110alpha led to embryonic lethality at mid-gestation because of severe defects in angiogenic sprouting and vascular remodelling. p110alpha exerts this critical endothelial cell-autonomous function by regulating endothelial cell migration through the small GTPase RhoA. p110alpha activity is particularly high in endothelial cells and preferentially induced by tyrosine kinase ligands (such as vascular endothelial growth factor (VEGF)-A). In contrast, p110beta in endothelial cells signals downstream of G-protein-coupled receptor (GPCR) ligands such as SDF-1alpha, whereas p110delta is expressed at low level and contributes only minimally to PI3K activity in endothelial cells. These results provide the first in vivo evidence for p110-isoform selectivity in endothelial PI3K signalling during angiogenesis.

Control of axonal growth and regeneration of sensory neurons by the p110delta PI 3-kinase.

The expression and function of the 8 distinct catalytic isoforms of PI 3-kinase (PI3K) in the nervous system are unknown. Whereas most PI3Ks have a broad tissue distribution, the tyrosine kinase-linked p110delta isoform has previously been shown to be enriched in leukocytes. Here we report that p110delta is also highly expressed in the nervous system. Inactivation of p110delta in mice did not affect gross neuronal development but led to an increased vulnerability of dorsal root ganglia neurons to exhibit growth cone collapse and decreases in axonal extension. Loss of p110delta activity also dampened axonal regeneration following peripheral nerve injury in adult mice and impaired functional recovery of locomotion. p110delta inactivation resulted in reduced neuronal signaling through the Akt protein kinase, and increased activity of the small GTPase RhoA. Pharmacological inhibition of ROCK, a downstream effector of RhoA, restored axonal extension defects in neurons with inactive p110delta, suggesting a key role of RhoA in p110delta signaling in neurons. Our data identify p110delta as an important signaling component for efficient axonal elongation in the developing and regenerating nervous system.