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.

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.

Wnt signaling is intrinsic to mouse embryonic stem cell self-renewal. Therefore, it is surprising that reprogramming of somatic cells to induced pluripotent stem cells (iPSCs) is not strongly enhanced by Wnt signaling. Here, we demonstrate that active Wnt signaling inhibits the early stage of reprogramming to iPSCs, whereas it is required and even stimulating during the late stage. Mechanistically, this biphasic effect of Wnt signaling is accompanied by a change in the requirement of all four of its transcriptional effectors: T cell factor 1 (Tcf1), Lef1, Tcf3, and Tcf4. For example, Tcf3 and Tcf4 are stimulatory early but inhibitory late in the reprogramming process. Accordingly, ectopic expression of Tcf3 early in reprogramming combined with its loss of function late enables efficient reprogramming in the absence of ectopic Sox2. Together, our data indicate that the stepwise process of reprogramming to iPSCs is critically dependent on the stage-specific control and action of all four Tcfs and Wnt signaling.

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.

Oct4 is a key component of the molecular circuitry which regulates embryonic stem cell proliferation and differentiation. It is essential for maintenance of undifferentiated, pluripotent cell populations, and accomplishes these tasks by binding DNA in multiple heterodimer and homodimer configurations. Very little is known about how formation of these complexes is regulated, or the mechanisms through which Oct4 proteins respond to complex extracellular stimuli which regulate pluripotency. Here, we provide evidence for a phosphorylation-based mechanism which regulates specific Oct4 homodimer conformations. Point mutations of a putative phosphorylation site can specifically abrogate transcriptional activity of a specific homodimer assembly, with little effect on other configurations. Moreover, we performed bioinformatic predictions to identify a subset of Oct4 target genes which may be regulated by this specific assembly, and show that altering Oct4 protein levels affects transcription of Oct4 target genes which are regulated by this assembly but not others. Finally, we identified several signaling pathways which may mediate this phosphorylation and act in combination to regulate Oct4 transcriptional activity and protein stability. These results provide a mechanism for rapid and reversible alteration of Oct4 transactivation potential in response to extracellular signals.

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.

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.

Nuclear compartments have diverse roles in regulating gene expression, yet the molecular forces and components that drive compartment formation remain largely unclear1. The long non-coding RNA Xist establishes an intra-chromosomal compartment by localizing at a high concentration in a territory spatially close to its transcription locus2 and binding diverse proteins3-5 to achieve X-chromosome inactivation (XCI)6,7. The XCI process therefore serves as a paradigm for understanding how RNA-mediated recruitment of various proteins induces a functional compartment. The properties of the inactive X (Xi)-compartment are known to change over time, because after initial Xist spreading and transcriptional shutoff a state is reached in which gene silencing remains stable even if Xist is turned off8. Here we show that the Xist RNA-binding proteins PTBP19, MATR310, TDP-4311 and CELF112 assemble on the multivalent E-repeat element of Xist7 and, via self-aggregation and heterotypic protein-protein interactions, form a condensate1 in the Xi. This condensate is required for gene silencing and for the anchoring of Xist to the Xi territory, and can be sustained in the absence of Xist. Notably, these E-repeat-binding proteins become essential coincident with transition to the Xist-independent XCI phase8, indicating that the condensate seeded by the E-repeat underlies the developmental switch from Xist-dependence to Xist-independence. Taken together, our data show that Xist forms the Xi compartment by seeding a heteromeric condensate that consists of ubiquitous RNA-binding proteins, revealing an unanticipated mechanism for heritable gene silencing.

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.

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 lncRNA Xist forms ∼50 diffraction-limited foci to transcriptionally silence one X chromosome. How this small number of RNA foci and interacting proteins regulate a much larger number of X-linked genes is unknown. We show that Xist foci are locally confined, contain ∼2 RNA molecules, and nucleate supramolecular complexes (SMACs) that include many copies of the critical silencing protein SPEN. Aggregation and exchange of SMAC proteins generate local protein gradients that regulate broad, proximal chromatin regions. Partitioning of numerous SPEN molecules into SMACs is mediated by their intrinsically disordered regions and essential for transcriptional repression. Polycomb deposition via SMACs induces chromatin compaction and the increase in SMACs density around genes, which propagates silencing across the X chromosome. Our findings introduce a mechanism for functional nuclear compartmentalization whereby crowding of transcriptional and architectural regulators enables the silencing of many target genes by few RNA molecules.

High-field asymmetric waveform ion mobility spectrometry (FAIMS) enables gas-phase separations on a chromatographic time scale and has become a useful tool for proteomic applications. Despite its emerging utility, however, the molecular determinants underlying peptide separation by FAIMS have not been systematically investigated. Here, we characterize peptide transmission in a FAIMS device across a broad range of compensation voltages (CVs) and used machine learning to identify charge state and three-dimensional (3D) electrostatic peptide potential as major contributors to peptide intensity at a given CV. We also demonstrate that the machine learning model can be used to predict optimized CV values for peptides, which significantly improves parallel reaction monitoring workflows. Together, these data provide insight into peptide separation by FAIMS and highlight its utility in targeted proteomic applications.

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.

Cells secrete numerous bioactive molecules essential for the function of healthy organisms. However, there are no scalable methods to link individual cell secretions to their transcriptional state. By developing and using secretion encoded single-cell sequencing (SEC-seq), which exploits hydrogel nanovials to capture individual cells and their secretions, we simultaneously measured the secretion of vascular endothelial growth factor A (VEGF-A) and the transcriptome for thousands of individual mesenchymal stromal cells (MSCs). We found that VEGF-A secretion is heterogeneous across the cell population and lowly correlated with the VEGFA transcript level. While there is a modest population-wide increase in VEGF-A secretion by hypoxic induction, highest VEGF-A secretion across normoxic and hypoxic culture conditions occurs in a subpopulation of MSCs characterized by a unique gene expression signature. Taken together, SEC-seq enables the identification of specific genes involved in the control of secretory states, which may be exploited for developing means to modulate cellular secretion for disease treatment.

Dominantly inherited disorders are not typically considered to be therapeutic candidates for gene augmentation. Here, we utilized induced pluripotent stem cell-derived retinal pigment epithelium (iPSC-RPE) to test the potential of gene augmentation to treat Best disease, a dominant macular dystrophy caused by over 200 missense mutations in BEST1. Gene augmentation in iPSC-RPE fully restored BEST1 calcium-activated chloride channel activity and improved rhodopsin degradation in an iPSC-RPE model of recessive bestrophinopathy as well as in two models of dominant Best disease caused by different mutations in regions encoding ion-binding domains. A third dominant Best disease iPSC-RPE model did not respond to gene augmentation, but showed normalization of BEST1 channel activity following CRISPR-Cas9 editing of the mutant allele. We then subjected all three dominant Best disease iPSC-RPE models to gene editing, which produced premature stop codons specifically within the mutant BEST1 alleles. Single-cell profiling demonstrated no adverse perturbation of retinal pigment epithelium (RPE) transcriptional programs in any model, although off-target analysis detected a silent genomic alteration in one model. These results suggest that gene augmentation is a viable first-line approach for some individuals with dominant Best disease and that non-responders are candidates for alternate approaches such as gene editing. However, testing gene editing strategies for on-target efficiency and off-target events using personalized iPSC-RPE model systems is warranted. In summary, personalized iPSC-RPE models can be used to select among a growing list of gene therapy options to maximize safety and efficacy while minimizing time and cost. Similar scenarios likely exist for other genotypically diverse channelopathies, expanding the therapeutic landscape for affected individuals.

Changes in transcriptional regulatory networks can significantly alter cell fate. To gain insight into transcriptional dynamics, several studies have profiled bulk multi-omic data sets with parallel transcriptomic and epigenomic measurements at different stages of a developmental process. However, integrating these data to infer cell type-specific regulatory networks is a major challenge. We present dynamic regulatory module networks (DRMNs), a novel approach to infer cell type-specific cis-regulatory networks and their dynamics. DRMN integrates expression, chromatin state, and accessibility to predict cis-regulators of context-specific expression, where context can be cell type, developmental stage, or time point, and uses multitask learning to capture network dynamics across linearly and hierarchically related contexts. We applied DRMNs to study regulatory network dynamics in three developmental processes, each showing different temporal relationships and measuring a different combination of regulatory genomic data sets: cellular reprogramming, liver dedifferentiation, and forward differentiation. DRMN identified known and novel regulators driving cell type-specific expression patterns, showing its broad applicability to examine dynamics of gene regulatory networks from linearly and hierarchically related multi-omic data sets.

We performed RNA sequencing on 40,000 cells to create a high-resolution single-cell gene expression atlas of developing human cortex, providing the first single-cell characterization of previously uncharacterized cell types, including human subplate neurons, comparisons with bulk tissue, and systematic analyses of technical factors. These data permit deconvolution of regulatory networks connecting regulatory elements and transcriptional drivers to single-cell gene expression programs, significantly extending our understanding of human neurogenesis, cortical evolution, and the cellular basis of neuropsychiatric disease. We tie cell-cycle progression with early cell fate decisions during neurogenesis, demonstrating that differentiation occurs on a transcriptomic continuum; rather than only expressing a few transcription factors that drive cell fates, differentiating cells express broad, mixed cell-type transcriptomes before telophase. By mapping neuropsychiatric disease genes to cell types, we implicate dysregulation of specific cell types in ASD, ID, and epilepsy. We developed CoDEx, an online portal to facilitate data access and browsing.