During the reprogramming of mouse embryonic fibroblasts (MEFs) to induced pluripotent stem cells, the activation of pluripotency genes such as NANOG occurs after the mesenchymal to epithelial transition. Here we report that both adult stem cells (neural stem cells) and differentiated cells (astrocytes) of the neural lineage can activate NANOG in the absence of cadherin expression during reprogramming. Gene expression analysis revealed that only the NANOG+E-cadherin+ populations expressed stabilization markers, had upregulated several cell cycle genes; and were transgene independent. Inhibition of DOT1L activity enhanced both the numbers of NANOG+ and NANOG+E-cadherin+ colonies in neural stem cells. Expressing SOX2 in MEFs prior to reprogramming did not alter the ratio of NANOG colonies that express E-cadherin. Taken together these results provide a unique pathway for reprogramming taken by cells of the neural lineage.

Following injury, pathologically activated vocal fold fibroblasts (VFFs) can engage in disordered extracellular matrix (ECM) remodeling, leading to VF fibrosis and impaired voice function. Given the importance of scar VFFs to phenotypically appropriate in vitro modeling of VF fibrosis, we pursued detailed characterization of scar VFFs obtained from surgically injured rat VF mucosae, compared with those obtained from experimentally naïve, age-matched tissue. Scar VFFs initially exhibited a myofibroblast phenotype characterized by increased proliferation, increased Col1a1 transcription and collagen, type I synthesis, increased Acta2 transcription and α-smooth muscle actin synthesis, and enhanced contractile function. These features were most distinct at passage 1 (P1); we observed a coalescence of the scar and naïve VFF phenotypes at later passages. An empirical Bayes statistical analysis of the P1 cell transcriptome identified 421 genes that were differentially expressed by scar, compared with naïve, VFFs. These genes were primarily associated with the wound response, ECM regulation, and cell proliferation. Follow-up comparison of P1 scar VFFs and their in vivo tissue source showed substantial transcriptomic differences. Finally, P1 scar VFFs responded to treatment with hepatocyte growth factor and transforming growth factor-β3, two biologics with reported therapeutic value. Despite the practical limitations inherent to working with early passage cells, this experimental model is easily implemented in any suitably equipped laboratory and has the potential to improve the applicability of preclinical VF fibrosis research.

The recent identification of hemogenic endothelium (HE) in human pluripotent stem cell (hPSC) cultures presents opportunities to investigate signaling pathways that are essential for blood development from endothelium and provides an exploratory platform for de novo generation of hematopoietic stem cells (HSCs). However, the use of poorly defined human or animal components limits the utility of the current differentiation systems for studying specific growth factors required for HE induction and manufacturing clinical-grade therapeutic blood cells. Here, we identified chemically defined conditions required to produce HE from hPSCs growing in Essential 8 (E8) medium and showed that Tenascin C (TenC), an extracellular matrix protein associated with HSC niches, strongly promotes HE and definitive hematopoiesis in this system. hPSCs differentiated in chemically defined conditions undergo stages of development similar to those previously described in hPSCs cocultured on OP9 feeders, including the formation of VE-Cadherin(+)CD73(-)CD235a/CD43(-) HE and hematopoietic progenitors with myeloid and T lymphoid potential.

Long range regulatory interactions among distal enhancers and target genes are important for tissue-specific gene expression. Genome-scale identification of these interactions in a cell line-specific manner, especially using the fewest possible datasets, is a significant challenge. We develop a novel computational approach, Regulatory Interaction Prediction for Promoters and Long-range Enhancers (RIPPLE), that integrates published Chromosome Conformation Capture (3C) data sets with a minimal set of regulatory genomic data sets to predict enhancer-promoter interactions in a cell line-specific manner. Our results suggest that CTCF, RAD21, a general transcription factor (TBP) and activating chromatin marks are important determinants of enhancer-promoter interactions. To predict interactions in a new cell line and to generate genome-wide interaction maps, we develop an ensemble version of RIPPLE and apply it to generate interactions in five human cell lines. Computational validation of these predictions using existing ChIA-PET and Hi-C data sets showed that RIPPLE accurately predicts interactions among enhancers and promoters. Enhancer-promoter interactions tend to be organized into subnetworks representing coordinately regulated sets of genes that are enriched for specific biological processes and cis-regulatory elements. Overall, our work provides a systematic approach to predict and interpret enhancer-promoter interactions in a genome-wide cell-type specific manner using a few experimentally tractable measurements.

Previously, we have used clinical and gene expression data from The Cancer Genome Atlas (TCGA) to model a pathway-based index predicting outcomes in ovarian carcinoma. This data were obtained from snap-frozen tissue measured with the Affymetrix U133 platform. In the current study, we correlate the data used to model with data derived from TaqMan qPCR both snap frozen and paraffin embedded (FFPE) samples.

Advancing pluripotent stem cell technologies for modelling haematopoietic stem cell development and blood therapies requires identifying key regulators of haematopoietic commitment from human pluripotent stem cells (hPSCs). Here, by screening the effect of 27 candidate factors, we reveal two groups of transcriptional regulators capable of inducing distinct haematopoietic programs from hPSCs: pan-myeloid (ETV2 and GATA2) and erythro-megakaryocytic (GATA2 and TAL1). In both cases, these transcription factors directly convert hPSCs to endothelium, which subsequently transform into blood cells with pan-myeloid or erythro-megakaryocytic potential. These data demonstrate that two distinct genetic programs regulate the haematopoietic development from hPSCs and that both of these programs specify hPSCs directly to haemogenic endothelial cells. In addition, this study provides a novel method for the efficient induction of blood and endothelial cells from hPSCs via the overexpression of modified mRNA for the selected transcription factors.

Oscillatory gene expression is fundamental to development, but technologies for monitoring expression oscillations are limited. We have developed a statistical approach called Oscope to identify and characterize the transcriptional dynamics of oscillating genes in single-cell RNA-seq data from an unsynchronized cell population. Applying Oscope to a number of data sets, we demonstrated its utility and also identified a potential artifact in the Fluidigm C1 platform.

The salamander has the remarkable ability to regenerate its limb after amputation. Cells at the site of amputation form a blastema and then proliferate and differentiate to regrow the limb. To better understand this process, we performed deep RNA sequencing of the blastema over a time course in the axolotl, a species whose genome has not been sequenced. Using a novel comparative approach to analyzing RNA-seq data, we characterized the transcriptional dynamics of the regenerating axolotl limb with respect to the human gene set. This approach involved de novo assembly of axolotl transcripts, RNA-seq transcript quantification without a reference genome, and transformation of abundances from axolotl contigs to human genes. We found a prominent burst in oncogene expression during the first day and blastemal/limb bud genes peaking at 7 to 14 days. In addition, we found that limb patterning genes, SALL genes, and genes involved in angiogenesis, wound healing, defense/immunity, and bone development are enriched during blastema formation and development. Finally, we identified a category of genes with no prior literature support for limb regeneration that are candidates for further evaluation based on their expression pattern during the regenerative process.

Among the three embryonic germ layers, the mesoderm is a major source of the mesenchymal precursors giving rise to skeletal and connective tissues, but these precursors have not previously been identified and characterized. Using human embryonic stem cells directed toward mesendodermal differentiation, we show that mesenchymal stem/stromal cells (MSCs) originate from a population of mesodermal cells identified by expression of apelin receptor. In semisolid medium, these precursors form FGF2-dependent compact spheroid colonies containing mesenchymal cells with a transcriptional profile representative of mesoderm-derived embryonic mesenchyme. When transferred to adherent cultures, individual colonies give rise to MSC lines with chondro-, osteo-, and adipogenic differentiation potentials. Although the MSC lines lacked endothelial potential, endothelial cells could be derived from the mesenchymal colonies, suggesting that, similar to hematopoietic cells, MSCs arise from precursors with angiogenic potential. Together, these studies identified a common precursor of mesenchymal and endothelial cells, mesenchymoangioblast, as the source of mesoderm-derived MSCs.

An open chromatin largely devoid of heterochromatin is a hallmark of stem cells. It remains unknown whether an open chromatin is necessary for the differentiation potential of stem cells, and which molecules are needed to maintain open chromatin. Here we show that the chromatin remodelling factor Chd1 is required to maintain the open chromatin of pluripotent mouse embryonic stem cells. Chd1 is a euchromatin protein that associates with the promoters of active genes, and downregulation of Chd1 leads to accumulation of heterochromatin. Chd1-deficient embryonic stem cells are no longer pluripotent, because they are incapable of giving rise to primitive endoderm and have a high propensity for neural differentiation. Furthermore, Chd1 is required for efficient reprogramming of fibroblasts to the pluripotent stem cell state. Our results indicate that Chd1 is essential for open chromatin and pluripotency of embryonic stem cells, and for somatic cell reprogramming to the pluripotent state.

Although numerous quantitative trait loci (QTL) influencing disease-related phenotypes have been detected through gene mapping and positional cloning, identification of the individual gene(s) and molecular pathways leading to those phenotypes is often elusive. One way to improve understanding of genetic architecture is to classify phenotypes in greater depth by including transcriptional and metabolic profiling. In the current study, we have generated and analyzed mRNA expression and metabolic profiles in liver samples obtained in an F2 intercross between the diabetes-resistant C57BL/6 leptin(ob/ob) and the diabetes-susceptible BTBR leptin(ob/ob) mouse strains. This cross, which segregates for genotype and physiological traits, was previously used to identify several diabetes-related QTL. Our current investigation includes microarray analysis of over 40,000 probe sets, plus quantitative mass spectrometry-based measurements of sixty-seven intermediary metabolites in three different classes (amino acids, organic acids, and acyl-carnitines). We show that liver metabolites map to distinct genetic regions, thereby indicating that tissue metabolites are heritable. We also demonstrate that genomic analysis can be integrated with liver mRNA expression and metabolite profiling data to construct causal networks for control of specific metabolic processes in liver. As a proof of principle of the practical significance of this integrative approach, we illustrate the construction of a specific causal network that links gene expression and metabolic changes in the context of glutamate metabolism, and demonstrate its validity by showing that genes in the network respond to changes in glutamine and glutamate availability. Thus, the methods described here have the potential to reveal regulatory networks that contribute to chronic, complex, and highly prevalent diseases and conditions such as obesity and diabetes.

From bacteria to humans, individual cells within isogenic populations can show significant variation in stress tolerance, but the nature of this heterogeneity is not clear. To investigate this, we used single-cell RNA sequencing to quantify transcript heterogeneity in single Saccharomyces cerevisiae cells treated with and without salt stress to explore population variation and identify cellular covariates that influence the stress-responsive transcriptome. Leveraging the extensive knowledge of yeast transcriptional regulation, we uncovered significant regulatory variation in individual yeast cells, both before and after stress. We also discovered that a subset of cells appears to decouple expression of ribosomal protein genes from the environmental stress response in a manner partly correlated with the cell cycle but unrelated to the yeast ultradian metabolic cycle. Live-cell imaging of cells expressing pairs of fluorescent regulators, including the transcription factor Msn2 with Dot6, Sfp1, or MAP kinase Hog1, revealed both coordinated and decoupled nucleocytoplasmic shuttling. Together with transcriptomic analysis, our results suggest that cells maintain a cellular filter against decoupled bursts of transcription factor activation but mount a stress response upon coordinated regulation, even in a subset of unstressed cells.

During vertebrate development, progenitor cells give rise to tissues and organs through a complex choreography that commences at gastrulation. A hallmark event of gastrulation is the formation of the primitive streak, a linear assembly of cells along the anterior-posterior (AP) axis of the developing organism. To examine the primitive streak at a single-cell resolution, we measured the transcriptomes of individual chick cells from the streak or the surrounding tissue (the rest of the area pellucida) in Hamburger-Hamilton stage 4 embryos. The single-cell transcriptomes were then ordered by the statistical method Wave-Crest to deduce both the relative position along the AP axis and the prospective lineage of single cells. The ordered transcriptomes reveal intricate patterns of gene expression along the primitive streak.

The ability to quantify cellular heterogeneity is a major advantage of single-cell technologies. However, statistical methods often treat cellular heterogeneity as a nuisance. We present a novel method to characterize differences in expression in the presence of distinct expression states within and among biological conditions. We demonstrate that this framework can detect differential expression patterns under a wide range of settings. Compared to existing approaches, this method has higher power to detect subtle differences in gene expression distributions that are more complex than a mean shift, and can characterize those differences. The freely available R package scDD implements the approach.

Elucidating the mechanism of reprogramming is confounded by heterogeneity due to the low efficiency and differential kinetics of obtaining induced pluripotent stem cells (iPSCs) from somatic cells. Therefore, we increased the efficiency with a combination of epigenomic modifiers and signaling molecules and profiled the transcriptomes of individual reprogramming cells. Contrary to the established temporal order, somatic gene inactivation and upregulation of cell cycle, epithelial, and early pluripotency genes can be triggered independently such that any combination of these events can occur in single cells. Sustained co-expression of Epcam, Nanog, and Sox2 with other genes is required to progress toward iPSCs. Ehf, Phlda2, and translation initiation factor Eif4a1 play functional roles in robust iPSC generation. Using regulatory network analysis, we identify a critical role for signaling inhibition by 2i in repressing somatic expression and synergy between the epigenomic modifiers ascorbic acid and a Dot1L inhibitor for pluripotency gene activation.

Inhibiting the histone 3 lysine 79 (H3K79) methyltransferase, disruptor of telomeric silencing 1-like (DOT1L), increases the efficiency of reprogramming somatic cells to induced pluripotent stem cells (iPSCs). Here, we find that, despite the enrichment of H3K79 methylation on thousands of actively transcribed genes in somatic cells, DOT1L inhibition (DOT1Li) does not immediately cause the shutdown of the somatic transcriptional profile to enable transition to pluripotency. Contrary to the prevalent view, DOT1Li promotes iPSC generation beyond the mesenchymal to epithelial transition and even from already epithelial cell types. DOT1Li is most potent at the midpoint of reprogramming in part by repressing Nfix that persists at late stages of reprogramming. Importantly, regulation of single genes cannot substitute for DOT1Li, demonstrating that H3K79 methylation has pleiotropic effects in maintaining cell identity.

Human pluripotent stem cells hold significant promise for regenerative medicine. However, long differentiation protocols and immature characteristics of stem cell-derived cell types remain challenges to the development of many therapeutic applications. In contrast to the slow differentiation of human stem cells in vitro that mirrors a nine-month gestation period, mouse stem cells develop according to a much faster three-week gestation timeline. Here, we tested if co-differentiation with mouse pluripotent stem cells could accelerate the differentiation speed of human embryonic stem cells. Following a six-week RNA-sequencing time course of neural differentiation, we identified 929 human genes that were upregulated earlier and 535 genes that exhibited earlier peaked expression profiles in chimeric cell cultures than in human cell cultures alone. Genes with accelerated upregulation were significantly enriched in Gene Ontology terms associated with neurogenesis, neuron differentiation and maturation, and synapse signaling. Moreover, chimeric mixed samples correlated with in utero human embryonic samples earlier than human cells alone, and acceleration was dose-dependent on human-mouse co-culture ratios. The altered gene expression patterns and developmental rates described in this report have implications for accelerating human stem cell differentiation and the use of interspecies chimeric embryos in developing human organs for transplantation.

Comparative time series transcriptome analysis is a powerful tool to study development, evolution, aging, disease progression and cancer prognosis. We develop TimeMeter, a statistical method and tool to assess temporal gene expression similarity, and identify differentially progressing genes where one pattern is more temporally advanced than the other. We apply TimeMeter to several datasets, and show that TimeMeter is capable of characterizing complicated temporal gene expression associations. Interestingly, we find: (i) the measurement of differential progression provides a novel feature in addition to pattern similarity that can characterize early developmental divergence between two species; (ii) genes exhibiting similar temporal patterns between human and mouse during neural differentiation are under strong negative (purifying) selection during evolution; (iii) analysis of genes with similar temporal patterns in mouse digit regeneration and axolotl blastema differentiation reveals common gene groups for appendage regeneration with potential implications in regenerative medicine.

Development of treatments for vocal dysphonia has been inhibited by lack of human vocal fold (VF) mucosa models because of difficulty in procuring VF epithelial cells, epithelial cells' limited proliferative capacity and absence of cell lines. Here we report development of engineered VF mucosae from hiPSC, transfected via TALEN constructs for green fluorescent protein, that mimic development of VF epithelial cells in utero. Modulation of FGF signaling achieves stratified squamous epithelium from definitive and anterior foregut derived cultures. Robust culturing of these cells on collagen-fibroblast constructs produces three-dimensional models comparable to in vivo VF mucosa. Furthermore, we demonstrate mucosal inflammation upon exposure of these constructs to 5% cigarette smoke extract. Upregulation of pro-inflammatory genes in epithelium and fibroblasts leads to aberrant VF mucosa remodeling. Collectively, our results demonstrate that hiPSC-derived VF mucosa is a versatile tool for future investigation of genetic and molecular mechanisms underlying epithelium-fibroblasts interactions in health and disease.

Defects in somitogenesis result in vertebral malformations at birth known as spondylocostal dysostosis (SCDO). Somites are formed with a species-specific periodicity controlled by the "segmentation clock," which comprises a group of oscillatory genes in the presomitic mesoderm. Here, we report that a segmentation clock model derived from human embryonic stem cells shows many hallmarks of the mammalian segmentation clock in vivo, including a dependence on the NOTCH and WNT signaling pathways. The gene expression oscillations are highly synchronized, displaying a periodicity specific to the human clock. Introduction of a point of mutation into HES7, a specific mutation previously associated with clinical SCDO, eliminated clock gene oscillations, successfully reproducing the defects in the segmentation clock. Thus, we provide a model for studying the previously inaccessible human segmentation clock to better understand the mechanisms contributing to congenital skeletal defects.