Atlas of regulated target genes of transcription factors (ART-TF) in human ES cells.

BACKGROUND: Transcription factors (TFs) play central roles in maintaining "stemness" of embryonic stem (ES) cells and their differentiation into several hundreds of adult cell types. The regulatory competence of TFs is routinely assessed by detecting target genes to which they bind. However, these data do not indicate which target genes are activated, repressed, or not affected by the change of TF abundance. There is a lack of large-scale studies that compare the genome binding of TFs with the expression change of target genes after manipulation of each TF. RESULTS: In this paper we associated human TFs with their target genes by two criteria: binding to genes, evaluated from published ChIP-seq data (n = 1868); and change of target gene expression shortly after induction of each TF in human ES cells. Lists of direction- and strength-specific regulated target genes are generated for 311 TFs (out of 351 TFs tested) with expected proportion of false positives less than or equal to 0.30, including 63 new TFs not present in four existing databases of target genes. Our lists of direction-specific targets for 152 TFs (80.0%) are larger that in the TRRUST database. In average, 30.9% of genes that respond greater than or equal to twofold to the induction of TFs are regulated targets. Regulated target genes indicate that the majority of TFs are either strong activators or strong repressors, whereas sets of genes that responded greater than or equal to twofold to the induction of TFs did not show strong asymmetry in the direction of expression change. The majority of human TFs (82.1%) regulated their target genes primarily via binding to enhancers. Repression of target genes is more often mediated by promoter-binding than activation of target genes. Enhancer-promoter loops are more abundant among strong activator and repressor TFs. CONCLUSIONS: We developed an atlas of regulated targets of TFs (ART-TF) in human ES cells by combining data on TF binding with data on gene expression change after manipulation of individual TFs. Sets of regulated gene targets were identified with a controlled rate of false positives. This approach contributes to the understanding of biological functions of TFs and organization of gene regulatory networks. This atlas should be a valuable resource for ES cell-based regenerative medicine studies.

Generation and Profiling of 2,135 Human ESC Lines for the Systematic Analyses of Cell States Perturbed by Inducing Single Transcription Factors.

Transcription factors (TFs) play a pivotal role in determining cell states, yet our understanding of the causative relationship between TFs and cell states is limited. Here, we systematically examine the state changes of human pluripotent embryonic stem cells (hESCs) by the large-scale manipulation of single TFs. We establish 2,135 hESC lines, representing three clones each of 714 doxycycline (Dox)-inducible genes including 481 TFs, and obtain 26,998 microscopic cell images and 2,174 transcriptome datasets-RNA sequencing (RNA-seq) or microarrays-48 h after the presence or absence of Dox. Interestingly, the expression of essentially all the genes, including genes located in heterochromatin regions, are perturbed by these TFs. TFs are also characterized by their ability to induce differentiation of hESCs into specific cell lineages. These analyses help to provide a way of classifying TFs and identifying specific sets of TFs for directing hESC differentiation into desired cell types.

A single-cell transcriptome atlas of the adult human retina.

The retina is a specialized neural tissue that senses light and initiates image processing. Although the functional organization of specific retina cells has been well studied, the molecular profile of many cell types remains unclear in humans. To comprehensively profile the human retina, we performed single-cell RNA sequencing on 20,009 cells from three donors and compiled a reference transcriptome atlas. Using unsupervised clustering analysis, we identified 18 transcriptionally distinct cell populations representing all known neural retinal cells: rod photoreceptors, cone photoreceptors, Müller glia, bipolar cells, amacrine cells, retinal ganglion cells, horizontal cells, astrocytes, and microglia. Our data captured molecular profiles for healthy and putative early degenerating rod photoreceptors, and revealed the loss of MALAT1 expression with longer post-mortem time, which potentially suggested a novel role of MALAT1 in rod photoreceptor degeneration. We have demonstrated the use of this retina transcriptome atlas to benchmark pluripotent stem cell-derived cone photoreceptors and an adult Müller glia cell line. This work provides an important reference with unprecedented insights into the transcriptional landscape of human retinal cells, which is fundamental to understanding retinal biology and disease.

Composite Agency: Semiotics of Modularity and Guiding Interactions.

Principles of constructivism are used here to explore how organisms develop tools, subagents, scaffolds, signs, and adaptations. Here I discuss reasons why organisms have composite nature and include diverse subagents that interact in partially cooperating and partially conflicting ways. Such modularity is necessary for efficient and robust functionality, including mutual construction and adaptability at various time scales. Subagents interact via material and semiotic relations, some of which force or prescribe actions of partners. Other interactions, which I call "guiding", do not have immediate effects and do not disrupt the evolution and learning capacity of partner agents. However, they modify the extent of learning and evolutionary possibilities of partners via establishment of scaffolds and constraints. As a result, subagents construct reciprocal scaffolding for each other to rebalance their communal evolution and learning. As an example, I discuss guiding interactions between the body and mind of animals, where the pain system adjusts mind-based learning to the physical and physiological constraints of the body. Reciprocal effects of mind and behaviors on the development and evolution of the body includes the effects of Lamarck and Baldwin.

Evolution of natural agents: preservation, advance, and emergence of functional information.

Biological evolution is often viewed narrowly as a change of morphology or allele frequency in a sequence of generations. Here I pursue an alternative informational concept of evolution, as preservation, advance, and emergence of functional information in natural agents. Functional information is a network of signs (e.g., memory, transient messengers, and external signs) that are used by agents to preserve and regulate their functions. Functional information is preserved in evolution via complex interplay of copying and construction processes: the digital components are copied, whereas interpreting subagents together with scaffolds, tools, and resources, are constructed. Some of these processes are simple and invariant, whereas others are complex and contextual. Advance of functional information includes improvement and modification of already existing functions. Although the genome information may change passively and randomly, the interpretation is active and guided by the logic of agent behavior and embryonic development. Emergence of new functions is based on the reinterpretation of already existing information, when old tools, resources, and control algorithms are adopted for novel functions. Evolution of functional information progressed from protosemiosis, where signs correspond directly to actions, to eusemiosis, where agents associate signs with objects. Language is the most advanced form of eusemiosis, where the knowledge of objects and models is communicated between agents.

Induction of specific neuron types by overexpression of single transcription factors.

Specific neuronal types derived from embryonic stem cells (ESCs) can facilitate mechanistic studies and potentially aid in regenerative medicine. Existing induction methods, however, mostly rely on the effects of the combined action of multiple added growth factors, which generally tend to result in mixed populations of neurons. Here, we report that overexpression of specific transcription factors (TFs) in ESCs can rather guide the differentiation of ESCs towards specific neuron lineages. Analysis of data on gene expression changes 2 d after induction of each of 185 TFs implicated candidate TFs for further ESC differentiation studies. Induction of 23 TFs (out of 49 TFs tested) for 6 d facilitated neural differentiation of ESCs as inferred from increased proportion of cells with neural progenitor marker PSA-NCAM. We identified early activation of the Notch signaling pathway as a common feature of most potent inducers of neural differentiation. The majority of neuron-like cells generated by induction of Ascl1, Smad7, Nr2f1, Dlx2, Dlx4, Nr2f2, Barhl2, and Lhx1 were GABA-positive and expressed other markers of GABAergic neurons. In the same way, we identified Lmx1a and Nr4a2 as inducers for neurons bearing dopaminergic markers and Isl1, Fezf2, and St18 for cholinergic motor neurons. A time-course experiment with induction of Ascl1 showed early upregulation of most neural-specific messenger RNA (mRNA) and microRNAs (miRNAs). Sets of Ascl1-induced mRNAs and miRNAs were enriched in Ascl1 targets. In further studies, enrichment of cells obtained with the induction of Ascl1, Smad7, and Nr2f1 using microbeads resulted in essentially pure population of neuron-like cells with expression profiles similar to neural tissues and expressed markers of GABAergic neurons. In summary, this study indicates that induction of transcription factors is a promising approach to generate cultures that show the transcription profiles characteristic of specific neural cell types.

Emergence of undifferentiated colonies from mouse embryonic stem cells undergoing differentiation by retinoic acid treatment.

Retinoic acid (RA) is one of the most potent inducers of differentiation of mouse embryonic stem cells (ESCs). However, previous studies show that RA treatment of cells cultured in the presence of a leukemia inhibitory factor (LIF) also result in the upregulation of a gene called Zscan4, whose transient expression is a marker for undifferentiated ESCs. We explored the balance between these two seemingly antagonistic effects of RA. ESCs indeed differentiated in the presence of LIF after RA treatment, but colonies of undifferentiated ESCs eventually emerged from these differentiated cells - even in the presence of RA. These colonies, named secondary colonies, consist of three cell types: typical undifferentiated ESCs expressing pluripotency genes such as Pou5f1, Sox2, and Nanog; cells expressing Zscan4; and endodermal-like cells located at the periphery of the colony. The capacity to form secondary colonies was confirmed for all eight tested ESC lines. Cells from the secondary colonies - after transfer to the standard ESC medium - retained pluripotency, judged by their strong alkaline phosphatase (ALP) staining, typical colony morphology, gene expression profile, stable karyotype, capacity to differentiate into all three germ layers in embryoid body formation assays, and successful contribution to chimeras after injection into blastocysts. Based on flow cytometry analysis (FACS), the proportion of Zscan4-positive cells in secondary colonies was higher than in standard ESC colonies, which may explain the capacity of ESCs to resist the differentiating effects of RA and instead form secondary colonies of undifferentiated ESCs. This hypothesis is supported by cell-lineage tracing analysis, which showed that most cells in the secondary colonies were descendents of cells transiently expressing Zscan4.

Generation and gene expression profiling of 48 transcription-factor-inducible mouse embryonic stem cell lines.

Mouse embryonic stem cells (ESCs) can differentiate into a wide range - and possibly all cell types in vitro, and thus provide an ideal platform to study systematically the action of transcription factors (TFs) in cell differentiation. Previously, we have generated and analyzed 137 TF-inducible mouse ESC lines. As an extension of this "NIA Mouse ESC Bank," we generated and characterized 48 additional mouse ESC lines, in which single TFs in each line could be induced in a doxycycline-controllable manner. Together, with the previous ESC lines, the bank now comprises 185 TF-manipulable ESC lines (>10% of all mouse TFs). Global gene expression (transcriptome) profiling revealed that the induction of individual TFs in mouse ESCs for 48 hours shifts their transcriptomes toward specific differentiation fates (e.g., neural lineages by Myt1 Isl1, and St18; mesodermal lineages by Pitx1, Pitx2, Barhl2, and Lmx1a; white blood cells by Myb, Etv2, and Tbx6, and ovary by Pitx1, Pitx2, and Dmrtc2). These data also provide and lists of inferred target genes of each TF and possible functions of these TFs. The results demonstrate the utility of mouse ESC lines and their transcriptome data for understanding the mechanism of cell differentiation and the function of TFs.

Coenzyme world model of the origin of life.

The origin of life means the emergence of heritable and evolvable self-reproduction. However the mechanisms of primordial heredity were different from those in contemporary cells. Here I argue that primordial life had no nucleic acids; instead heritable signs were represented by isolated catalytically active self-reproducing molecules, similar to extant coenzymes, which presumably colonized surfaces of oil droplets in water. The model further assumes that coenzyme-like molecules (CLMs) changed surface properties of oil droplets (e.g., by oxidizing terminal carbons), and in this way created and sustained favorable conditions for their own self-reproduction. Such niche-dependent self-reproduction is a necessary condition for cooperation between different kinds of CLMs because they have to coexist in the same oil droplets and either succeed or perish together. Additional kinds of hereditary molecules were acquired via coalescence of oil droplets carrying different kinds of CLMs or via modification of already existing CLMs. Eventually, polymerization of CLMs became controlled by other polymers used as templates; and this kind of template-based synthesis eventually resulted in the emergence of RNA-like replicons. Apparently, oil droplets transformed into the outer membrane of cells via engulfing water, stabilization of the surface, and osmoregulation. In result, the metabolism was internalized allowing cells to accumulate free-floating resources (e.g., animoacids, ATP), which was a necessary condition for the development of protein synthesis. Thus, life originated from simple but already functional molecules, and its gradual evolution towards higher complexity was driven by cooperation and natural selection.

Evolutionary biosemiotics and multilevel construction networks.

In contrast to the traditional relational semiotics, biosemiotics decisively deviates towards dynamical aspects of signs at the evolutionary and developmental time scales. The analysis of sign dynamics requires constructivism (in a broad sense) to explain how new components such as subagents, sensors, effectors, and interpretation networks are produced by developing and evolving organisms. Semiotic networks that include signs, tools, and subagents are multilevel, and this feature supports the plasticity, robustness, and evolvability of organisms. The origin of life is described here as the emergence of simple self-constructing semiotic networks that progressively increased the diversity of their components and relations. Primitive organisms have no capacity to classify and track objects; thus, we need to admit the existence of proto-signs that directly regulate activities of agents without being associated with objects. However, object recognition and handling became possible in eukaryotic species with the development of extensive rewritable epigenetic memory as well as sensorial and effector capacities. Semiotic networks are based on sequential and recursive construction, where each step produces components (i.e., agents, scaffolds, signs, and resources) that are needed for the following steps of construction. Construction is not limited to repair and reproduction of what already exists or is unambiguously encoded, it also includes production of new components and behaviors via learning and evolution. A special case is the emergence of new levels of organization known as metasystem transition. Multilevel semiotic networks reshape the phenotype of organisms by combining a mosaic of features developed via learning and evolution of cooperating and/or conflicting subagents.

Transient bursts of Zscan4 expression are accompanied by the rapid derepression of heterochromatin in mouse embryonic stem cells.

Mouse embryonic stem cells (mESCs) have a remarkable capacity to maintain normal genome stability and karyotype in culture. We previously showed that infrequent bursts of Zscan4 expression (Z4 events) are important for the maintenance of telomere length and genome stability in mESCs. However, the molecular details of Z4 events remain unclear. Here we show that Z4 events involve unexpected transcriptional derepression in heterochromatin regions that usually remain silent. During a Z4 event, we see rapid derepression and rerepression of heterochromatin leading to a burst of transcription that coincides with transient histone hyperacetylation and DNA demethylation, clustering of pericentromeric heterochromatin around the nucleolus, and accumulation of activating and repressive chromatin remodelling complexes. This heterochromatin-based transcriptional activity suggests that mESCs may maintain their extraordinary genome stability at least in part by transiently resetting their heterochromatin.

ExAtlas: An interactive online tool for meta-analysis of gene expression data.

We have developed ExAtlas, an on-line software tool for meta-analysis and visualization of gene expression data. In contrast to existing software tools, ExAtlas compares multi-component data sets and generates results for all combinations (e.g. all gene expression profiles versus all Gene Ontology annotations). ExAtlas handles both users' own data and data extracted semi-automatically from the public repository (GEO/NCBI database). ExAtlas provides a variety of tools for meta-analyses: (1) standard meta-analysis (fixed effects, random effects, z-score, and Fisher's methods); (2) analyses of global correlations between gene expression data sets; (3) gene set enrichment; (4) gene set overlap; (5) gene association by expression profile; (6) gene specificity; and (7) statistical analysis (ANOVA, pairwise comparison, and PCA). ExAtlas produces graphical outputs, including heatmaps, scatter-plots, bar-charts, and three-dimensional images. Some of the most widely used public data sets (e.g. GNF/BioGPS, Gene Ontology, KEGG, GAD phenotypes, BrainScan, ENCODE ChIP-seq, and protein-protein interaction) are pre-loaded and can be used for functional annotations.

Protosemiosis: agency with reduced representation capacity.

Life has semiotic nature; and as life forms differ in their complexity, functionality, and adaptability, we assume that forms of semiosis also vary accordingly. Here we propose a criterion to distinguish between the primitive kind of semiosis, which we call "protosemiosis" (following Prodi) from the advanced kind of semiosis, or "eusemiosis". In protosemiosis, agents associate signs directly with actions without considering objects, whereas in eusemiosis, agents associate signs with objects and only then possibly with actions. Protosemiosis started from the origin of life, and eusemiosis started when evolving agents acquired the ability to track and classify objects. Eusemiosis is qualitatively different from protosemiosis because it can not be reduced to a small number of specific signaling pathways. Proto-signs can be classified into proto-icons that signal via single specific interaction, proto-indexes that combine several functions, and proto-symbols that are processed by a universal subagent equipped with a set of heritable adapters. Prefix "proto" is used here to characterize signs at the protosemiotic level. Although objects are not recognized by protosemiotic agents, they can be reliably reconstructed by human observers. In summary, protosemiosis is a primitive kind of semiosis that supports "know-how" without "know-what". Without studying protosemiosis, the biosemiotics theory would be incomplete.

Evolutionary constraints or opportunities?

Natural selection is traditionally viewed as a leading factor of evolution, whereas variation is assumed to be random and non-directional. Any order in variation is attributed to epigenetic or developmental constraints that can hinder the action of natural selection. In contrast I consider the positive role of epigenetic mechanisms in evolution because they provide organisms with opportunities for rapid adaptive change. Because the term "constraint" has negative connotations, I use the term "regulated variation" to emphasize the adaptive nature of phenotypic variation, which helps populations and species to survive and evolve in changing environments. The capacity to produce regulated variation is a phenotypic property, which is not described in the genome. Instead, the genome acts as a switchboard, where mostly random mutations switch "on" or "off" preexisting functional capacities of organism components. Thus, there are two channels of heredity: informational (genomic) and structure-functional (phenotypic). Functional capacities of organisms most likely emerged in a chain of modifications and combinations of more simple ancestral functions. The role of DNA has been to keep records of these changes (without describing the result) so that they can be reproduced in the following generations. Evolutionary opportunities include adjustments of individual functions, multitasking, connection between various components of an organism, and interaction between organisms. The adaptive nature of regulated variation can be explained by the differential success of lineages in macro-evolution. Lineages with more advantageous patterns of regulated variation are likely to produce more species and secure more resources (i.e., long-term lineage selection).

Inferring directions of evolution from patterns of variation: the legacy of Sergei Meyen.

In the era of the extended evolutionary synthesis, which no longer considers natural selection as the only leading factor of evolution, it is meaningful to revisit the legacy of biologists who discussed the role of alternative factors. Here we analyze the evolutionary views of Sergei Meyen (1935-1987), a paleobotanist who argued that the theory of evolution should incorporate a "nomothetical" approach which infers the laws of morphogenesis (i.e., form generation) from the observed patterns of variation in living organisms and in the fossil records. Meyen developed a theory of "repeated polymorphic sets" (RPSs), which he applied consistently to describe inter-organism variation in populations, intra-organism variation of metameric organs, variation of abnormalities, heterotopy, changes during embryo development, and inter-species variation within evolutionary lineages. The notion of RPS assumes the active nature of organisms that possess hidden morphogenic and behavioral capacities. Meyen's theory is compatible with Darwin's natural selection; however, Meyen emphasized the importance of other forms of selection (e.g., selection of developmental trajectories, habitats, and behaviors) in choosing specific elements from the RPS. Finally, Meyen developed a new typological concept of time, where time represents variability (i.e., change) of real objects such as living organisms or geological formations.

Chromatin properties of regulatory DNA probed by manipulation of transcription factors.

Transcription factors (TFs) bind to DNA and regulate the transcription of nearby genes. However, only a small fraction of TF binding sites have such regulatory effects. Here we search for the predictors of functional binding sites by carrying out a systematic computational screening of a variety of contextual factors (histone modifications, nuclear lamin-bindings, and cofactor bindings). We used regression analysis to test if contextual factors are associated with upregulation or downregulation of neighboring genes following the induction or knockdown of the 9 TFs in mouse embryonic stem (ES) cells. Functional TF binding sites appeared to be either active (i.e., bound by P300, CHD7, mediator, cohesin, and SWI/SNF) or repressed (i.e., with H3K27me3 histone marks and bound by Polycomb factors). Active binding sites mediated the downregulation of nearby genes upon knocking down the activating TFs or inducing repressors. Repressed TF binding sites mediated the upregulation of nearby genes (e.g., poised developmental regulators) upon inducing TFs. In addition, repressed binding sites mediated repressive effects of TFs, identified by the downregulation of target genes after the induction of TFs or by the upregulation of target genes after the knockdown of TFs. The contextual factors associated with functions of DNA-bound TFs were used to improve the identification of candidate target genes regulated by TFs.

Evolutionary constraints or opportunities?

Natural selection is traditionally viewed as a leading factor of evolution, whereas variation is assumed to be random and non-directional. Any order in variation is attributed to epigenetic or developmental constraints that can hinder the action of natural selection. In contrast I consider the positive role of epigenetic mechanisms in evolution because they provide organisms with opportunities for rapid adaptive change. Because the term "constraint" has negative connotations, I use the term "regulated variation" to emphasize the adaptive nature of phenotypic variation, which helps populations and species to survive and evolve in changing environments. The capacity to produce regulated variation is a phenotypic property, which is not described in the genome. Instead, the genome acts as a switchboard, where mostly random mutations switch "on" or "off" preexisting functional capacities of organism components. Thus, there are two channels of heredity: informational (genomic) and structure-functional (phenotypic). Functional capacities of organisms most likely emerged in a chain of modifications and combinations of more simple ancestral functions. The role of DNA has been to keep records of these changes (without describing the result) so that they can be reproduced in the following generations. Evolutionary opportunities include adjustments of individual functions, multitasking, connection between various components of an organism, and interaction between organisms. The adaptive nature of regulated variation can be explained by the differential success of lineages in macro-evolution. Lineages with more advantageous patterns of regulated variation are likely to produce more species and secure more resources (i.e., long-term lineage selection).

Identification of transcription factors for lineage-specific ESC differentiation.

A network of transcription factors (TFs) determines cell identity, but identity can be altered by overexpressing a combination of TFs. However, choosing and verifying combinations of TFs for specific cell differentiation have been daunting due to the large number of possible combinations of ∼2,000 TFs. Here, we report the identification of individual TFs for lineage-specific cell differentiation based on the correlation matrix of global gene expression profiles. The overexpression of identified TFs-Myod1, Mef2c, Esx1, Foxa1, Hnf4a, Gata2, Gata3, Myc, Elf5, Irf2, Elf1, Sfpi1, Ets1, Smad7, Nr2f1, Sox11, Dmrt1, Sox9, Foxg1, Sox2, or Ascl1-can direct efficient, specific, and rapid differentiation into myocytes, hepatocytes, blood cells, and neurons. Furthermore, transfection of synthetic mRNAs of TFs generates their appropriate target cells. These results demonstrate both the utility of this approach to identify potent TFs for cell differentiation, and the unanticipated capacity of single TFs directly guides differentiation to specific lineage fates.

A conserved Oct4/POUV-dependent network links adhesion and migration to progenitor maintenance.

BACKGROUND: The class V POU domain transcription factor Oct4 (Pou5f1) is a pivotal regulator of embryonic stem cell (ESC) self-renewal and reprogramming of somatic cells to induced pluripotent stem (iPS) cells. Oct4 is also an important evolutionarily conserved regulator of progenitor cell differentiation during embryonic development. RESULTS: Here we examine the function of Oct4 homologs in Xenopus embryos and compare this to the role of Oct4 in maintaining mammalian embryo-derived stem cells. Based on a combination of expression profiling of Oct4/POUV-depleted Xenopus embryos and in silico analysis of existing mammalian Oct4 target data sets, we defined a set of evolutionary-conserved Oct4/POUV targets. Most of these targets were regulators of cell adhesion. This is consistent with Oct4/POUV phenotypes observed in the adherens junctions in Xenopus ectoderm, mouse embryonic, and epiblast stem cells. A number of these targets could rescue both Oct4/POUV phenotypes in cellular adhesion and multipotent progenitor cell maintenance, whereas expression of cadherins on their own could only transiently support adhesion and block differentiation in both ESC and Xenopus embryos. CONCLUSIONS: Currently, the list of Oct4 transcriptional targets contains thousands of genes. Using evolutionary conservation, we identified a core set of functionally relevant factors that linked the maintenance of adhesion to Oct4/POUV. We found that the regulation of adhesion by the Oct4/POUV network occurred at both transcriptional and posttranslational levels and was required for pluripotency.