Supplementary Table 1 from HIF Induces Human Embryonic Stem Cell Markers in Cancer Cells

Supplementary Table 7 from HIF Induces Human Embryonic Stem Cell Markers in Cancer Cells

Supplementary Methods, Figures 1-11, Table Legends 1-7 from HIF Induces Human Embryonic Stem Cell Markers in Cancer Cells

Supplementary Table 2 from HIF Induces Human Embryonic Stem Cell Markers in Cancer Cells

Supplementary Table 5 from HIF Induces Human Embryonic Stem Cell Markers in Cancer Cells

Supplementary Table 3 from HIF Induces Human Embryonic Stem Cell Markers in Cancer Cells

Supplementary Table 4 from HIF Induces Human Embryonic Stem Cell Markers in Cancer Cells

Supplementary Table 6 from HIF Induces Human Embryonic Stem Cell Markers in Cancer Cells

Embryonic development is a continuum in vivo. Transcriptional analysis can separate established human embryonic stem cells (hESC) into at least four distinct developmental pluripotent stages, two naïve and two primed, early and late relative to the intact epiblast. In this study we primarily show that exposure of frozen human blastocysts to an inhibitor of checkpoint kinase 1 (CHK1) upon thaw greatly enhances establishment of karyotypically normal late naïve hESC cultures. These late naïve cells are plastic and can be toggled back to early naïve and forward to early primed pluripotent stages. The early primed cells are transcriptionally equivalent to the post inner cell mass intermediate (PICMI) stage seen one day following transfer of human blastocysts into in vitro culture and are stable at an earlier stage than conventional primed hESC.

Regulation of embryonic diapause, dormancy that interrupts the tight connection between developmental stage and time, is still poorly understood. Here, we characterize the transcriptional and metabolite profiles of mouse diapause embryos and identify unique gene expression and metabolic signatures with activated lipolysis, glycolysis, and metabolic pathways regulated by AMPK. Lipolysis is increased due to mTORC2 repression, increasing fatty acids to support cell survival. We further show that starvation in pre-implantation ICM-derived mouse ESCs induces a reversible dormant state, transcriptionally mimicking the in vivo diapause stage. During starvation, Lkb1, an upstream kinase of AMPK, represses mTOR, which induces a reversible glycolytic and epigenetically H4K16Ac-negative, diapause-like state. Diapause furthermore activates expression of glutamine transporters SLC38A1/2. We show by genetic and small molecule inhibitors that glutamine transporters are essential for the H4K16Ac-negative, diapause state. These data suggest that mTORC1/2 inhibition, regulated by amino acid levels, is causal for diapause metabolism and epigenetic state.

Single-cell Hi-C (scHi-C) interrogates genome-wide chromatin interaction in individual cells, allowing us to gain insights into 3D genome organization. However, the extremely sparse nature of scHi-C data poses a significant barrier to analysis, limiting our ability to tease out hidden biological information. In this work, we approach this problem by applying topic modeling to scHi-C data. Topic modeling is well-suited for discovering latent topics in a collection of discrete data. For our analysis, we generate nine different single-cell combinatorial indexed Hi-C (sci-Hi-C) libraries from five human cell lines (GM12878, H1Esc, HFF, IMR90, and HAP1), consisting over 19,000 cells. We demonstrate that topic modeling is able to successfully capture cell type differences from sci-Hi-C data in the form of "chromatin topics." We further show enrichment of particular compartment structures associated with locus pairs in these topics.

During mammalian embryogenesis, changes in morphology and gene expression are concurrent with epigenomic reprogramming. Using human embryonic stem cells representing the preimplantation blastocyst (naive) and postimplantation epiblast (primed), our data in 2iL/I/F naive cells demonstrate that a substantial portion of known human enhancers are premarked by H3K4me1, providing an enhanced open chromatin state in naive pluripotency. The 2iL/I/F enhancer repertoire occupies 9% of the genome, three times that of primed cells, and can exist in broad chromatin domains over 50 kb. Enhancer chromatin states are largely poised. Seventy-seven percent of 2iL/I/F enhancers are decommissioned in a stepwise manner as cells become primed. While primed topologically associating domains are largely unaltered upon differentiation, naive 2iL/I/F domains expand across primed boundaries, affecting three-dimensional genome architecture. Differential topologically associating domain edges coincide with 2iL/I/F H3K4me1 enrichment. Our results suggest that naive-derived 2iL/I/F cells have a unique chromatin landscape, which may reflect early embryogenesis.

Chromatin immunoprecipitation (ChIP) is the most widely used approach for identification of genome-associated proteins and their modifications. We have previously introduced a microplate-based ChIP platform, Matrix ChIP, where the entire ChIP procedure is done on the same plate without sample transfers. Compared to conventional ChIP protocols, the Matrix ChIP assay is faster and has increased throughput. However, even with microplate ChIP assays, sample preparation and chromatin fragmentation (which is required to map genomic locations) remains a major bottleneck. We have developed a novel technology (termed 'PIXUL') utilizing an array of ultrasound transducers for simultaneous shearing of samples in standard 96-well microplates. We integrated PIXUL with Matrix ChIP ('PIXUL-ChIP'), that allows for fast, reproducible, low-cost and high-throughput sample preparation and ChIP analysis of 96 samples (cell culture or tissues) in one day. Further, we demonstrated that chromatin prepared using PIXUL can be used in an existing ChIP-seq workflow. Thus, the high-throughput capacity of PIXUL-ChIP provides the means to carry out ChIP-qPCR or ChIP-seq experiments involving dozens of samples. Given the complexity of epigenetic processes, the use of PIXUL-ChIP will advance our understanding of these processes in health and disease, as well as facilitate screening of epigenetic drugs.

Efficient stem cell differentiation into pancreatic islet cells is of critical importance for the development of cell replacement therapies for diabetes. Here, we identify the expression pattern of connexin 43 (Cx43), a gap junction (GJ) channel protein, in human embryonic stem cell (hESC)-derived definitive endoderm (DE) and primitive gut tube cells, representing early lineages for posterior foregut (PF), pancreatic progenitors (PP), pancreatic endocrine progenitors (PE), and islet cells. As the function of GJ channels is dependent on their gating status, we tested the impact of supplementing hESC-derived PP cell cultures with AAP10, a peptide that promotes Cx43 GJ channel opening. We found that this treatment promotes the expression of DE markers FoxA2 and Sox17, leads to a more efficient derivation of DE, and improves the yield of PF, PP, and PE cells. These results demonstrate a functional involvement of GJ channels in the differentiation of embryonic stem cells into pancreatic cell lineages.

To easily edit the genome of naïve human embryonic stem cells (hESC), we introduced a dual cassette encoding an inducible Cas9 into the AAVS1 site of naïve hESC (iCas9). The iCas9 line retained karyotypic stability, expression of pluripotency markers, differentiation potential, and stability in 5iLA and EPS pluripotency conditions. The iCas9 line induced efficient homology-directed repair (HDR) and non-homologous end joining (NHEJ) based mutations through CRISPR-Cas9 system. We utilized the iCas9 line to study the epigenetic regulator, PRC2 in early human pluripotency. The PRC2 requirement distinguishes between early pluripotency stages, however, what regulates PRC2 activity in these stages is not understood. We show reduced H3K27me3 and pluripotency markers in JARID2 2iL-I-F hESC mutants, indicating JARID2 requirement in maintenance of hESC 2iL-I-F state. These data suggest that JARID2 regulates PRC2 in 2iL-I-F state and the lack of PRC2 function in 5iLA state may be due to lack of sufficient JARID2 protein.

The folding and three-dimensional (3D) organization of chromatin in the nucleus critically impacts genome function. The past decade has witnessed rapid advances in genomic tools for delineating 3D genome architecture. Among them, chromosome conformation capture (3C)-based methods such as Hi-C are the most widely used techniques for mapping chromatin interactions. However, traditional Hi-C protocols rely on restriction enzymes (REs) to fragment chromatin and are therefore limited in resolution. We recently developed DNase Hi-C for mapping 3D genome organization, which uses DNase I for chromatin fragmentation. DNase Hi-C overcomes RE-related limitations associated with traditional Hi-C methods, leading to improved methodological resolution. Furthermore, combining this method with DNA capture technology provides a high-throughput approach (targeted DNase Hi-C) that allows for mapping fine-scale chromatin architecture at exceptionally high resolution. Hence, targeted DNase Hi-C will be valuable for delineating the physical landscapes of cis-regulatory networks that control gene expression and for characterizing phenotype-associated chromatin 3D signatures. Here, we provide a detailed description of method design and step-by-step working protocols for these two methods.

In storing and transmitting epigenetic information, organisms must balance the need to maintain information about past conditions with the capacity to respond to information in their current and future environments. Some of this information is encoded by DNA methylation, which can be transmitted with variable fidelity from parent to daughter strand. High fidelity confers strong pattern matching between the strands of individual DNA molecules and thus pattern stability over rounds of DNA replication; lower fidelity confers reduced pattern matching, and thus greater flexibility. Here, we present a new conceptual framework, Ratio of Concordance Preference (RCP), that uses double-stranded methylation data to quantify the flexibility and stability of the system that gave rise to a given set of patterns. We find that differentiated mammalian cells operate with high DNA methylation stability, consistent with earlier reports. Stem cells in culture and in embryos, in contrast, operate with reduced, albeit significant, methylation stability. We conclude that preference for concordant DNA methylation is a consistent mode of information transfer, and thus provides epigenetic stability across cell divisions, even in stem cells and those undergoing developmental transitions. Broader application of our RCP framework will permit comparison of epigenetic-information systems across cells, developmental stages, and organisms whose methylation machineries differ substantially or are not yet well understood.

DNA methylation has been studied extensively in many developmental systems. Little attention, however, has been given to identifying methylation features that distinguish loci whose patterns are in flux in a given cell lineage from those whose patterns are stable. Here, we develop a new metric, the Ratio of Concordance Preference (RCP), to quantify and compare epigenetic flexibility and stability across loci, cell types, and developmental stages, without assuming any specific biochemical mechanisms. We apply RCP to double-stranded DNA methylation data from human and murine cells and conclude that: (i) preference for concordant DNA methylation is reduced but not eliminated in stem relative to differentiated cells; (ii) cellular differentiation is characterized by increasing preference for concordant methylation states; and (iii) while concordance preference remains substantial through embryonic totipotency and early stages of pluripotency, primordial germ cells initially have nearly no preference for concordance, perhaps reflecting the high level of epigenetic flexibility en route to production of gametes. The mechanism-free nature of RCP will enable comparison of DNA methylation systems not only across cell types and developmental stages, but also across organisms whose methylation machineries are not well understood or may differ significantly.