Aging of organ and abnormal tissue regeneration are recurrent problems in physiological and pathophysiological conditions. This is most crucial in case of high-turnover tissues, like bone marrow (BM). Using reciprocal transplantation experiments in mouse, we have shown that self-renewal potential of hematopoietic stem and progenitor cells (HSPCs) and BM cellularity are markedly influenced with the age of the recipient mice rather than donor mice. Moreover, accumulation of excessive reactive oxygen species (ROS) in BM stromal cells compared to HSPC compartment, in time-dependent manner, suggests that oxidative stress is involved in suppression of BM cellularity by affecting microenvironment in aged mice. Treatment of these mice with a polyphenolic antioxidant curcumin is found to partially quench ROS, thereby rescues stromal cells from oxidative stress-dependent cellular injury. This rejuvenation of stromal cells significantly improves hematopoietic reconstitution in 18-month-old mice compared to age control mice. In conclusion, this study implicates the role of ROS in perturbation of stromal cell function upon aging, which in turn affects BM's reconstitution ability in aged mice. Thus, a rejuvenation therapy using curcumin, before HSPC transplantation, is found to be an efficient strategy for successful marrow reconstitution in older mice.

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.

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.

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.

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.

Computational prediction of protein interactions typically use protein domains as classifier features because they capture conserved information of interaction surfaces. However, approaches relying on domains as features cannot be applied to proteins without any domain information. In this paper, we explore the contribution of pure amino acid composition (AAC) for protein interaction prediction. This simple feature, which is based on normalized counts of single or pairs of amino acids, is applicable to proteins from any sequenced organism and can be used to compensate for the lack of domain information.

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.

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.

Seminal studies of vertebrate protein evolution speculated that gene regulatory changes can drive anatomical innovations. However, very little is known about gene regulatory network (GRN) evolution associated with phenotypic effect across ecologically diverse species. Here we use a novel approach for comparative GRN analysis in vertebrate species to study GRN evolution in representative species of the most striking examples of adaptive radiations, the East African cichlids. We previously demonstrated how the explosive phenotypic diversification of East African cichlids can be attributed to diverse molecular mechanisms, including accelerated regulatory sequence evolution and gene expression divergence.

Leading CRISPR-Cas technologies employ Cas9 and Cas12 enzymes that generate RNA-guided dsDNA breaks. Yet, the most abundant microbial adaptive immune systems, Type I CRISPRs, are under-exploited for eukaryotic applications. Here, we report the adoption of a minimal CRISPR-Cas3 from Neisseria lactamica (Nla) type I-C system to create targeted large deletions in the human genome. RNP delivery of its processive Cas3 nuclease and target recognition complex Cascade can confer ∼95% editing efficiency. Unexpectedly, NlaCascade assembly in bacteria requires internal translation of a hidden component Cas11 from within the cas8 gene. Furthermore, expressing a separately encoded NlaCas11 is the key to enable plasmid- and mRNA-based editing in human cells. Finally, we demonstrate that supplying cas11 is a universal strategy to systematically implement divergent I-C, I-D, and I-B CRISPR-Cas3 editors with compact sizes, distinct PAM preferences, and guide orthogonality. These findings greatly expand our ability to engineer long-range genome edits.

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.

Cancer risk by environmental exposure is modulated by an individual's genetics and age at exposure. This age-specific period of susceptibility is referred to as the "Window of Susceptibility" (WOS). Rats have a similar WOS for developing breast cancer. A previous study in rat identified an age-specific long-range regulatory interaction for the cancer gene, Pappa, that is associated with breast cancer susceptibility. However, the global role of three-dimensional genome organization and downstream gene expression programs in the WOS is not known. Therefore, we generated Hi-C and RNA-seq data in rat mammary epithelial cells within and outside the WOS. To systematically identify higher-order changes in 3D genome organization, we developed NE-MVNMF that combines network enhancement followed by multitask non-negative matrix factorization. We examined three-dimensional genome organization dynamics at the level of individual loops as well as higher-order domains. Differential chromatin interactions tend to be associated with differentially up-regulated genes with the WOS and recapitulate several human SNP-gene interactions associated with breast cancer susceptibility. Our approach identified genomic blocks of regions with greater overall differences in contact count between the two time points when the cluster assignments change and identified genes and pathways implicated in early carcinogenesis and cancer treatment. Our results suggest that WOS-specific changes in 3D genome organization are linked to transcriptional changes that may influence susceptibility to breast cancer.

The divergence of regulatory regions and gene regulatory network (GRN) rewiring is a key driver of cichlid phenotypic diversity. However, the contribution of miRNA-binding site turnover has yet to be linked to GRN evolution across cichlids. Here, we extend our previous studies by analyzing the selective constraints driving evolution of miRNA and transcription factor (TF)-binding sites of target genes, to infer instances of cichlid GRN rewiring associated with regulatory binding site turnover. Comparative analyses identified increased species-specific networks that are functionally associated to traits of cichlid phenotypic diversity. The evolutionary rewiring is associated with differential models of miRNA- and TF-binding site turnover, driven by a high proportion of fast-evolving polymorphic sites in adaptive trait genes compared with subsets of random genes. Positive selection acting upon discrete mutations in these regulatory regions is likely to be an important mechanism in rewiring GRNs in rapidly radiating cichlids. Regulatory variants of functionally associated miRNA- and TF-binding sites of visual opsin genes differentially segregate according to phylogeny and ecology of Lake Malawi species, identifying both rewired, for example, clade-specific and conserved network motifs of adaptive trait associated GRNs. Our approach revealed several novel candidate regulators, regulatory regions, and three-node motifs across cichlid genomes with previously reported associations to known adaptive evolutionary traits.

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.

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.