Oct4, Sox2, Klf4, and cMyc (OSKM) reprogram somatic cells to pluripotency. To gain a mechanistic understanding of their function, we mapped OSKM-binding, stage-specific transcription factors (TFs), and chromatin states in discrete reprogramming stages and performed loss- and gain-of-function experiments. We found that OSK predominantly bind active somatic enhancers early in reprogramming and immediately initiate their inactivation genome-wide by inducing the redistribution of somatic TFs away from somatic enhancers to sites elsewhere engaged by OSK, recruiting Hdac1, and repressing the somatic TF Fra1. Pluripotency enhancer selection is a stepwise process that also begins early in reprogramming through collaborative binding of OSK at sites with high OSK-motif density. Most pluripotency enhancers are selected later in the process and require OS and other pluripotency TFs. Somatic and pluripotency TFs modulate reprogramming efficiency when overexpressed by altering OSK targeting, somatic-enhancer inactivation, and pluripotency enhancer selection. Together, our data indicate that collaborative interactions among OSK and with stage-specific TFs direct both somatic-enhancer inactivation and pluripotency-enhancer selection to drive reprogramming.

The epidermis and its appendage, the hair follicle, represent an elegant developmental system in which cells are replenished with regularity because of controlled proliferation, lineage specification, and terminal differentiation. Although transcriptome data exists for human epidermal and dermal cells, the hair follicle remains poorly characterized. Through single-cell resolution profiling of the epidermis and anagen hair follicle, we characterized the anatomical, transcriptional, functional, and pathological profiles of distinct epidermal, hair follicle, and hair follicle-associated cell subpopulations including melanocytes, endothelial cells, and immune cells. We additionally traced the differentiation trajectory of interfollicular and matrix cell progenitors and explored the association of specific cell subpopulations to known molecular signatures of common skin conditions. These data simultaneously corroborate prior murine and human studies while offering new insights into epidermal and hair follicle differentiation and pathogenesis.

Naive human embryonic stem cells (hESCs) can be derived from primed hESCs or directly from blastocysts, but their X chromosome state has remained unresolved. Here, we show that the inactive X chromosome (X_i) of primed hESCs was reactivated in naive culture conditions. Like cells of the blastocyst, the resulting naive cells contained two active X chromosomes with XIST expression and chromosome-wide transcriptional dampening and initiated XIST-mediated X inactivation upon differentiation. Both establishment of and exit from the naive state (differentiation) happened via an XIST-negative X_aX_a intermediate. Together, these findings identify a cell culture system for functionally exploring the two X chromosome dosage compensation processes in early human development: X dampening and X inactivation. However, remaining differences between naive hESCs and embryonic cells related to mono-allelic XIST expression and non-random X inactivation highlight the need for further culture improvement. As the naive state resets X_i abnormalities seen in primed hESCs, it may provide cells better suited for downstream applications.

The developmental trajectory of human skeletal myogenesis and the transition between progenitor and stem cell states are unclear. We used single-cell RNA sequencing to profile human skeletal muscle tissues from embryonic, fetal, and postnatal stages. In silico, we identified myogenic as well as other cell types and constructed a "roadmap" of human skeletal muscle ontogeny across development. In a similar fashion, we also profiled the heterogeneous cell cultures generated from multiple human pluripotent stem cell (hPSC) myogenic differentiation protocols and mapped hPSC-derived myogenic progenitors to an embryonic-to-fetal transition period. We found differentially enriched biological processes and discovered co-regulated gene networks and transcription factors present at distinct myogenic stages. This work serves as a resource for advancing our knowledge of human myogenesis. It also provides a tool for a better understanding of hPSC-derived myogenic progenitors for translational applications in skeletal muscle-based regenerative medicine.

Human embryonic stem cells (hESCs) can be captured in a primed state in which they resemble the postimplantation epiblast, or in a naive state where they resemble the preimplantation epiblast. Naive-cell-specific culture conditions allow the study of preimplantation development ex vivo but reportedly lead to chromosomal abnormalities, which compromises their utility in research and potential therapeutic applications. Although MEK inhibition is essential for the naive state, here we show that reduced MEK inhibition facilitated the establishment and maintenance of naive hESCs that retained naive-cell-specific features, including global DNA hypomethylation, HERVK expression, and two active X chromosomes. We further show that hESCs cultured under these modified conditions proliferated more rapidly; accrued fewer chromosomal abnormalities; and displayed changes in the phosphorylation levels of MAPK components, regulators of DNA damage/repair, and cell cycle. We thus provide a simple modification to current methods that can enable robust growth and reduced genomic instability in naive hESCs.

(Cell Stem Cell 27, 158–176.e1–e10; July 2, 2020) It came to our attention that during revision we updated the tSNE plot for the left panel of Figure 1D but mistakenly did not update the right panel of Figure 1D or Figure S1D. This error on our part does not affect any of the results or conclusions of our paper. The abovementioned figures have now been updated with the correct tSNE coordinates and reproduced below. We apologize for this oversight and any inconvenience the readers might have encountered. [Figure presented] [Formula presented][Formula presented] [Formula presented]

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.

The specification of inhibitory neurons has been described for the mouse and human brain, and many studies have shown that pluripotent stem cells (PSCs) can be used to create interneurons in vitro. It is unclear whether in vitro methods to produce human interneurons generate all the subtypes found in brain, and how similar in vitro and in vivo interneurons are. We applied single-nuclei and single-cell transcriptomics to model interneuron development from human cortex and interneurons derived from PSCs. We provide a direct comparison of various in vitro interneuron derivation methods to determine the homogeneity achieved. We find that PSC-derived interneurons capture stages of development prior to mid-gestation, and represent a minority of potential subtypes found in brain. Comparison with those found in fetal or adult brain highlighted decreased expression of synapse-related genes. These analyses highlight the potential to tailor the method of generation to drive formation of particular subtypes.

Background: Human embryonic stem cell-derived cardiomyocytes (hESC-CMs) can be used as a source for cell delivery to remuscularize the heart after myocardial infarction. Despite their therapeutic potential, the emergence of ventricular arrhythmias has limited their application. We previously developed a double reporter hESC line to isolate first heart field (FHF: TBX5 ⁺ NKX2-5 ⁺) and second heart field (SHF: TBX5 ^- NKX2-5 ⁺ ) CMs. Herein, we explore the role of TBX5 and its effects on underlying gene regulatory networks driving phenotypical and functional differences between these two populations. Methods: We used a combination of tools and techniques for rapid and unsupervised profiling of FHF and SHF populations at the transcriptional, translational, and functional level including single cell RNA (scRNA) and bulk RNA sequencing, atomic force and quantitative phase microscopy, respirometry, and electrophysiology. Results: Gene ontology analysis revealed three biological processes attributed to TBX5 expression: sarcomeric structure, oxidative phosphorylation, and calcium ion handling. Interestingly, migratory pathways were enriched in SHF population. SHF-like CMs display less sarcomeric organization compared to FHF-like CMs, despite prolonged in vitro culture. Atomic force and quantitative phase microscopy showed increased cellular stiffness and decreased mass distribution over time in FHF compared to SHF populations, respectively. Electrophysiological studies showed longer plateau in action potentials recorded from FHF-like CMs, consistent with their increased expression of calcium handling genes. Interestingly, both populations showed nearly identical respiratory profiles with the only significant functional difference being higher ATP generation-linked oxygen consumption rate in FHF-like CMs. Our findings suggest that FHF-like CMs display more mature features given their enhanced sarcomeric alignment, calcium handling, and decreased migratory characteristics. Finally, pseudotime analyses revealed a closer association of the FHF population to human fetal CMs along the developmental trajectory. Conclusion: Our studies reveal that distinguishing FHF and SHF populations based on TBX5 expression leads to a significant impact on their downstream functional properties. FHF CMs display more mature characteristics such as enhanced sarcomeric organization and improved calcium handling, with closer positioning along the differentiation trajectory to human fetal hearts. These data suggest that the FHF CMs may be a more suitable candidate for cardiac regeneration.