Evolution of non-coding regulatory sequences involved in the developmental process: reflection of differential employment of paralogous genes as highlighted by Sox2 and group B1 Sox genes

The Origin and Regulation of Neuromesodermal Progenitors (NMPs) in Embryos

Neuromesodermal progenitors (NMPs), serving as the common origin of neural and paraxial mesodermal development in a large part of the trunk, have recently gained significant attention because of their critical importance in the understanding of embryonic organogenesis and the design of in vitro models of organogenesis. However, the nature of NMPs at many essential points remains only vaguely understood or even incorrectly assumed. Here, we discuss the nature of NMPs, focusing on their dynamic migratory behavior during embryogenesis and the mechanisms underlying their neural vs. mesodermal fate choice. The discussion points include the following: (1) How the sinus rhomboidals is organized; the tissue where the neural or mesodermal fate choice of NMPs occurs. (2) NMPs originating from the broad posterior epiblast are associated with Sox2 N1 enhancer activity. (3) Tbx6-dependent Sox2 repression occurs during NMP-derived paraxial mesoderm development. (4) The nephric mesenchyme, a component of the intermediate mesoderm, was newly identified as an NMP derivative. (5) The transition of embryonic tissue development from tissue-specific progenitors in the anterior part to that from NMPs occurs at the forelimb bud axial level. (6) The coexpression of Sox2 and Bra in NMPs is conditional and is not a hallmark of NMPs. (7) The ability of the NMP pool to sustain axial embryo growth depends on Wnt3a signaling in the NMP population. Current in vitro models of NMPs are also critically reviewed.

Wnt signal-dependent antero-posterior specification of early-stage CNS primordia modeled in EpiSC-derived neural stem cells

The specification of the embryonic central nervous system (CNS) into future brain (forebrain, midbrain, or hindbrain) and spinal cord (SC) regions is a critical step of CNS development. A previous chicken embryo study indicated that anterior epiblast cells marked by Sox2 N2 enhancer activity are specified to the respective brain regions during the transition phase of the epiblast to the neural plate-forming neural primordium. The present study showed that the SC precursors positioned posterior to the hindbrain precursors in the anterior epiblast migrated posteriorly in contrast to the anterior migration of brain precursors. The anteroposterior specification of the CNS precursors occurs at an analogous time (∼E7.5) in mouse embryos, in which an anterior-to-posterior incremental gradient of Wnt signal strength was observed. To examine the possible Wnt signal contribution to the anteroposterior CNS primordium specification, we utilized mouse epiblast stem cell (EpiSC)-derived neurogenesis in culture. EpiSCs maintained in an activin- and FGF2-containing medium start neural development after the removal of activin, following a day in a transitory state. We placed activin-free EpiSCs in EGF- and FGF2-containing medium to arrest neural development and expand the cells into neural stem cells (NSCs). Simultaneously, a Wnt antagonist or agonist was added to the culture, with the anticipation that different levels of Wnt signals would act on the transitory cells to specify CNS regionality; then, the Wnt-treated cells were expanded as NSCs. Gene expression profiles of six NSC lines were analyzed using microarrays and single-cell RNA-seq. The NSC lines demonstrated anteroposterior regional specification in response to increasing Wnt signal input levels: forebrain-midbrain-, hindbrain-, cervical SC-, and thoracic SC-like lines. The regional coverage of these NSC lines had a range; for instance, the XN1 line expressed Otx2 and En2, indicating midbrain characteristics, but additionally expressed the SC-characteristic Hoxa5. The ranges in the anteroposterior specification of neural primordia may be narrowed as neural development proceeds. The thoracic SC is presumably the posterior limit of the contribution by anterior epiblast-derived neural progenitors, as the characteristics of more posterior SC regions were not displayed.

Enhancer Arrays Regulating Developmental Genes: Sox2 Enhancers as a Paradigm

Enhancers are the primary regulatory DNA sequences in eukaryotes and are mostly located in the non-coding sequences of genes, namely, intergenic regions and introns. The essential characteristic of an enhancer is the ability to activate proximal genes, e.g., a reporter gene in a reporter assay, regardless of orientation, relative position, and distance from the gene. These characteristics are ascribed to the interaction (spatial proximity) of the enhancer sequence and the gene promoter via DNA looping, discussed in the latter part of this chapter.Developmentally regulated genes are associated with multiple enhancers carrying distinct cell and developmental stage specificities, which form arrays on the genome. We discuss the array of enhancers regulating the Sox2 gene as a paradigm. Sox2 enhancers are the best studied enhancers of a single gene in developmental regulation. In addition, the Sox2 gene is located in a genomic region with a very sparse gene distribution (no other protein-coding genes in ~1.6 Mb in the mouse genome), termed a "gene desert," which means that most identified enhancers in the region are associated with Sox2 regulation. Furthermore, the importance of the Sox2 gene in stem cell regulation and neural development justifies focusing on Sox2-associated enhancers.

Expression and significance of SOX B1 genes in glioblastoma multiforme patients

The overall survival of glioblastoma multiforme (GBM) patients remains poor. To improve patient outcomes, effective diagnostic and prognostic biomarkers for GBM are needed. In this study, we first applied bioinformatic analyses to identify biomarkers for GBM, focusing on SOX (sex‐determining region on the Y chromosome (SRY)‐related high mobility group (HMG) box) B1 family members. The ONCOMINE, GEPIA, LinkedOmics and CCLE databases were used to assess mRNA expression levels of the SOX B1 family members in different cancers and normal tissue. Further bioinformatic analysis was performed using the ONCOMINE database in combination with the LinkedOmics data set to identify the prognostic value of SOX B1 family members for GBM. We found mRNA expression levels of all tested SOX B1 genes were significantly increased in GBM. In the LinkedOmics database, increased expression of SOX3 indicated a better overall survival. In GEPIA databases, increased expression of all SOX B1 family members suggested an improved overall survival, but none of them were statistically different. Then, Transwell assays and wound healing were employed to evaluate the motility and invasive captivity of U251 cells when silencing SOX2 and SOX3. We found exogenous inhibition of SOX2 appeared to reduce the migration and invasion of U251 cells in vitro. Collectively, our research suggested that SOX2 might serve as a cancer‐promoting gene to identify high‐risk GBM patients, and SOX3 had the potential to be a prognostic biomarker for GBM patients.

Neural Differentiation of Mouse Embryonic Stem Cells-An in vitro Approach to Profile DNA Methylation of Reprogramming Factor Sox2-SRR2

Sox2 is one of the core transcription factors maintaining the embryonic stem cells (ES) pluripotency and, also indispensable for cellular reprogramming. However, limited data is available about the DNA methylation of pluripotency genes during lineage-specific differentiations. This study investigated the DNA methylation of Sox2 regulatory region 2 (SRR2) during directed differentiation of mouse ES into neural lineage. ES cells were first grown to form embryoid bodies in suspension which were then dissociated, and cultured in defined medium to promote neural differentiation. Typical neuronal morphology together with the up-regulation of Pax6, neuroepithelial stem cell intermediate filament and β-tubulin III and, down-regulation of pluripotency genes Oct4, Nanog and Sox2 showed the existence of neural phenotype in cells undergoing differentiation. Three CpGs in the core enhancer region of neural-specific SRR2 were individually investigated by direct DNA sequencing post-bisulfite treatment and, found to be unmethylated in differentiated cells at time-points chosen for analysis. This analysis does not limit the possibility of methylation at other CpG sites than those profiled here and/or transient methylation. Hence, similar analyses exploring the DNA methylation at other regions of the Sox2 gene could unravel the onset and transitions of epigenetic signatures influencing the outcome of differentiation pathways and neural development. The data presented here shows that in vitro neural differentiation of embryonic stem cells can be employed to study and characterize molecular regulatory mechanisms governing neurogenesis by applying diverse pharmacological and toxicological agents.

ERK Activity Dynamics during Zebrafish Embryonic Development

During vertebrate development, extracellular signal-regulated kinase (ERK) is activated by growth factors such as fibroblast growth factor (FGF), and it regulates the formation of tissues/organs including eyes, brains, somites, limbs, and inner ears. However, an experimental system to monitor ERK activity dynamics in the entire body of the vertebrate embryo is lacking. We recently studied ERK activity dynamics in the pre-somitic mesoderm of living zebrafish embryos injected with mRNAs encoding a Förster resonance energy transfer (FRET)-based ERK biosensor. In this study, transgenic zebrafish stably and ubiquitously expressing the ERK biosensor were generated to monitor ERK activity dynamics throughout embryonic development. The system allowed the identification of ERK activation domains in embryos from the late blastula to the late segmentation stage, consistent with immunostaining patterns obtained using anti-phosphorylated ERK antibody. A spatiotemporal map of ERK activity in the entire body during zebrafish embryogenesis was generated, and previously unidentified activation dynamics and ERK domains were identified. The proposed system is the first reported method to monitor ERK activity dynamics during vertebrate embryogenesis, providing insight into the role of ERK activity in normal and abnormal development in living vertebrate embryos.

SoxB2 in sea urchin development: implications in neurogenesis, ciliogenesis and skeletal patterning

Background

Current studies in evolutionary developmental biology are focused on the reconstruction of gene regulatory networks in target animal species. From decades, the scientific interest on genetic mechanisms orchestrating embryos development has been increasing in consequence to the fact that common features shared by evolutionarily distant phyla are being clarified. In 2011, a study across eumetazoan species showed for the first time the existence of a highly conserved non-coding element controlling the SoxB2 gene, which is involved in the early specification of the nervous system. This discovery raised several questions about SoxB2 function and regulation in deuterostomes from an evolutionary point of view.

Results

Due to the relevant phylogenetic position within deuterostomes, the sea urchin Strongylocentrotus purpuratus represents an advantageous animal model in the field of evolutionary developmental biology. Herein, we show a comprehensive study of SoxB2 functions in sea urchins, in particular its expression pattern in a wide range of developmental stages, and its co-localization with other neurogenic markers, as SoxB1, SoxC and Elav. Moreover, this work provides a detailed description of the phenotype of sea urchin SoxB2 knocked-down embryos, confirming its key function in neurogenesis and revealing, for the first time, its additional roles in oral and aboral ectoderm cilia and skeletal rod morphology.

Conclusions

We concluded that SoxB2 in sea urchins has a neurogenic function; however, this gene could have multiple roles in sea urchin embryogenesis, expanding its expression in non-neurogenic cells. We showed that SoxB2 is functionally conserved among deuterostomes and suggested that in S. purpuratus this gene acquired additional functions, being involved in ciliogenesis and skeletal patterning.

Electronic supplementary material

The online version of this article (10.1186/s13227-018-0094-1) contains supplementary material, which is available to authorized users.

Axial level-dependent molecular and cellular mechanisms underlying the genesis of the embryonic neural plate

The transcription factor gene Sox2, centrally involved in neural primordial regulation, is activated by many enhancers. During the early stages of embryonic development, Sox2 is regulated by the enhancers N2 and N1 in the anterior neural plate (ANP) and posterior neural plate (PNP), respectively. This differential use of the enhancers reflects distinct regulatory mechanisms underlying the genesis of ANP and PNP. The ANP develops directly from the epiblast, triggered by nodal signal inhibition, and via the combined action of TFs SOX2, OTX2, POU3F1, and ZIC2, which promotes the the ANP development and inhibits other cell lineages. In contrast, the PNP is derived from neuromesodermal bipotential axial stem cells that develop into the neural plate when Sox2 is activated by the N1 enhancer, whereas they develop into the paraxial mesoderm when the N1 enhancer is repressed by the action of TBX6. The axial stem cells are maintained by the activity of WNT3a and T (Brachyury). However, at axial levels more anterior to the 8th somites (cervical levels), the development of both the neural plate and somite proceeds in the absence of WNT3a, T, or TBX6. These observations indicate that distinct molecular and cellular mechanisms determine neural plate genesis based on the axial level, and contradict the classical concept of the term "neural induction," which assumes a pan-neural plate mechanism.

Epigenetic activation of Sox2 gene in the developing vertebrate neural plate

The in vivo requirement of the histone demethylase JmjD2A, together with the kinase MSK1, results in a series of epigenetic events necessary for early activation of Sox2 and subsequent neural fate commitment in vertebrates.

One of the earliest manifestations of neural induction is onset of expression of the neural marker Sox2, mediated by the activation of the enhancers N1 and N2. By using loss and gain of function, we find that Sox2 expression requires the activity of JmjD2A and the Msk1 kinase, which can respectively demethylate the repressive H3K9me3 mark and phosphorylate the activating H3S10 (H3S10ph) mark. Bimolecular fluorescence complementation reveals that the adaptor protein 14-3-3, known to bind to H3S10ph, interacts with JMJD2A and may be involved in its recruitment to regulatory regions of the Sox2 gene. Chromatin immunoprecipitation reveals dynamic binding of JMJD2A to the Sox2 promoter and N-1 enhancer at the time of neural plate induction. Finally, we show a clear temporal antagonism on the occupancy of H3K9me3 and H3S10ph modifications at the promoter of the Sox2 locus before and after the neural plate induction. Taken together, our results propose a series of epigenetic events necessary for the early activation of the Sox2 gene in neural progenitor cells.

Transcription, Signaling Receptor Activity, Oxidative Phosphorylation, and Fatty Acid Metabolism Mediate the Presence of Closely Related Species in Distinct Intertidal and Cold-Seep Habitats

Bathyal cold seeps are isolated extreme deep-sea environments characterized by low species diversity while biomass can be high. The Håkon Mosby mud volcano (Barents Sea, 1,280 m) is a rather stable chemosynthetic driven habitat characterized by prominent surface bacterial mats with high sulfide concentrations and low oxygen levels. Here, the nematode Halomonhystera hermesi thrives in high abundances (11,000 individuals 10 cm⁻²). Halomonhystera hermesi is a member of the intertidal Halomonhystera disjuncta species complex that includes five cryptic species (GD1-5). GD1-5’s common habitat is characterized by strong environmental fluctuations. Here, we compared the transcriptomes of H. hermesi and GD1, H. hermesi’s closest relative. Genes encoding proteins involved in oxidative phosphorylation are more strongly expressed in H. hermesi than in GD1, and many genes were only observed in H. hermesi while being completely absent in GD1. Both observations could in part be attributed to high sulfide concentrations and low oxygen levels. Additionally, fatty acid elongation was also prominent in H. hermesi confirming the importance of highly unsaturated fatty acids in this species. Significant higher amounts of transcription factors and genes involved in signaling receptor activity were observed in GD1 (many of which were completely absent in H. hermesi), allowing fast signaling and transcriptional reprogramming which can mediate survival in dynamic intertidal environments. GC content was approximately 8% higher in H. hermesi coding unigenes resulting in differential codon usage between both species and a higher proportion of amino acids with GC-rich codons in H. hermesi. In general our results showed that most pathways were active in both environments and that only three genes are under natural selection. This indicates that also plasticity should be taken in consideration in the evolutionary history of Halomonhystera species. Such plasticity, as well as possible preadaptation to low oxygen and high sulfide levels might have played an important role in the establishment of a cold-seep Halomonhystera population.

First evidence of molecular characterization of rohu carp Sox2 gene being expressed in proliferating spermatogonial cells

Because little is known about the function of Sox2 (Sry-related box-2) in teleosts, the objective of this study was to clone and characterize Sox2 complementary DNA (cDNA) from the testis of Indian major carp, Labeo rohita (rohu). The full-length cDNA contained an open reading frame of 936 nucleotides bearing the typical structural features. Phylogenetically, Sox2 of L rohita was most closely related to freshwater counterparts than marine water. The sequence information of cDNA and genomic DNA together revealed that the Sox2 gene is encoded by an uninterrupted exon. Furthermore, comparative mRNA expression profile in various organs including proliferating spermatogonial stem cells (SSCs) suggested about the participatory role of Sox2 during fish male germ cell development and maintenance of stem cells. In support, we have also provided evidence that Sox2 protein is indeed present in rohu SSCs by Western blot analysis. The evolutionarily conserved high-mobility group box domain indicated its possible involvement in common networking pathways for stem cell maintenance and pluripotency between mammals and nonmammals. Our findings could be the first step toward the use of Sox2 as a potential biomarker for proliferating SSCs and understanding the transcriptional regulatory network involved during male germ cell development and maintenance in fish species.

Wnt/β-catenin and FGF signalling direct the specification and maintenance of a neuromesodermal axial progenitor in ensembles of mouse embryonic stem cells

The development of the central nervous system is known to result from two sequential events. First, an inductive event of the mesoderm on the overlying ectoderm that generates a neural plate that, after rolling into a neural tube, acts as the main source of neural progenitors. Second, the axial regionalization of the neural plate that will result in the specification of neurons with different anteroposterior identities. Although this description of the process applies with ease to amphibians and fish, it is more difficult to confirm in amniote embryos. Here, a specialized population of cells emerges at the end of gastrulation that, under the influence of Wnt and FGF signalling, expands and generates the spinal cord and the paraxial mesoderm. This population is known as the long-term neuromesodermal precursor (NMp). Here, we show that controlled increases of Wnt/β-catenin and FGF signalling during adherent culture differentiation of mouse embryonic stem cells (mESCs) generates a population with many of the properties of the NMp. A single-cell analysis of gene expression within this population reveals signatures that are characteristic of stem cell populations. Furthermore, when this activation is triggered in three-dimensional aggregates of mESCs, the population self-organizes macroscopically and undergoes growth and axial elongation that mimics some of the features of the embryonic spinal cord and paraxial mesoderm. We use both adherent and three-dimensional cultures of mESCs to probe the establishment and maintenance of NMps and their differentiation.

Sox proteins: regulators of cell fate specification and differentiation

Sox transcription factors play widespread roles during development; however, their versatile funtions have a relatively simple basis: the binding of a Sox protein alone to DNA does not elicit transcriptional activation or repression, but requires binding of a partner transcription factor to an adjacent site on the DNA. Thus, the activity of a Sox protein is dependent upon the identity of its partner factor and the context of the DNA sequence to which it binds. In this Primer, we provide an mechanistic overview of how Sox family proteins function, as a paradigm for transcriptional regulation of development involving multi-transcription factor complexes, and we discuss how Sox factors can thus regulate diverse processes during development.

A Systematic Survey and Characterization of Enhancers that Regulate Sox3 in Neuro-Sensory Development in Comparison with Sox2 Enhancers

Development of neural and sensory primordia at the early stages of embryogenesis depends on the activity of two B1 Sox transcription factors, Sox2 and Sox3. The embryonic expression patterns of the Sox2 and Sox3 genes are similar, yet they show gene-unique features. We screened for enhancers of the 231-kb genomic region encompassing Sox3 of chicken, and identified 13 new enhancers that showed activity in different domains of the neuro-sensory primordia. Combined with the three Sox3-proximal enhancers determined previously, at least 16 enhancers were involved in Sox3 regulation. Starting from the NP1 enhancer, more enhancers with different specificities are activated in sequence, resulting in complex overlapping patterns of enhancer activities. NP1 was activated in the caudal lateral epiblast adjacent to the posterior growing end of neural plate, and by the combined action of Wnt and Fgf signaling, similar to the Sox2 N1 enhancer involved in neural/mesodermal dichotomous cell lineage segregation. The Sox3 D5 enhancer and Sox2 N3 enhancer were also activated similarly in the diencephalon, optic vesicle and lens placode, suggesting analogies in their regulation. In general, however, the specificities of the enhancers were not identical between Sox3 and Sox2, including the cases of the NP1 and D5 enhancers.

Do not put all teleosts in one net: focus on the sox2 and pou2 genes

The Pou2 and Sox2 proteins are major transcription factors for development and cell differentiation. In teleosts, the expression patterns of pou2 or sox2 are different between species from distant families, suggesting different regulatory mechanisms of gene expression. In this study, we assessed the divergences among teleosts, including within closely related species. The pou2 and sox2 gene expression patterns were characterised over several developmental stages in a cyprinid model, i.e., the goldfish, and the potential regulation sites of these genes within teleost conserved regions were localised. During embryonic development, differences in the expression patterns between the goldfish and other teleosts, including zebrafish, were observed for both genes. The in silico analysis of the 5' flanking regions of the pou2 gene showed high conservation within teleosts, whereas the sox2 sequence diverged in tetraodontiforms. Certain putative cis regulatory elements were common to all teleosts, whereas others were found only in cyprinids. The analysis of the DNA methylation patterns of the pou2 and sox2 upstream sequences revealed that the studied CpG sites remained hypomethylated at all stages of embryo development in both genes. In contrast, in the adult fin, the studied CpG sites were hypermethylated in pou2 but not in sox2, suggesting the existence of methylation-sensitive regions in pou2. Overall, although most similarities at the level of the gene regulatory sites were found within cyprinids, the expression pattern of pou2 or sox2 during development differs between cyprinids species.

Application of a novel strategy of engineering conditional alleles to a single exon gene, Sox2

BACKGROUND

The Conditional by Inversion (COIN) method for engineering conditional alleles relies on an invertible optimized gene trap-like element, the COIN module, for imparting conditionality. The COIN module contains an optimized 3' splice site-polyadenylation signal pair, but is inserted antisense to the target gene and therefore does not alter transcription, until it is inverted by Cre recombinase. In order to make COIN applicable to all protein-coding genes, the COIN module has been engineered within an artificial intron, enabling insertion into an exon.

METHODOLOGY/PRINCIPAL FINDINGS

Therefore, theoretically, the COIN method should be applicable to single exon genes, and to test this idea we engineered a COIN allele of Sox2. This single exon gene presents additional design challenges, in that its proximal promoter and coding region are entirely contained within a CpG island, and are also spanned by an overlapping transcript, Sox2Ot, which contains mmu-miR1897. Here, we show that despite disruption of the CpG island by the COIN module intron, the COIN allele of Sox2 (Sox2(COIN)) is phenotypically wild type, and also does not interfere with expression of Sox2Ot and miR1897. Furthermore, the inverted COIN allele of Sox2, Sox2(INV) is functionally null, as homozygotes recapitulate the phenotype of Sox2(βgeo/βgeo) mice, a well-characterized Sox2 null. Lastly, the benefit of the eGFP marker embedded in the COIN allele is demonstrated as it mirrors the expression pattern of Sox2.

CONCLUSIONS/SIGNIFICANCE

Our results demonstrate the applicability of the COIN technology as a method of choice for targeting single exon genes.

Visualization and exploration of conserved regulatory modules using ReXSpecies 2

Background

The prediction of transcription factor binding sites is difficult for many reasons. Thus, filtering methods are needed to enrich for biologically relevant (true positive) matches in the large amount of computational predictions that are frequently generated from promoter sequences.

Results

ReXSpecies 2 filters predictions of transcription factor binding sites and generates a set of figures displaying them in evolutionary context. More specifically, it uses position specific scoring matrices to search for motifs that specify transcription factor binding sites. It removes redundant matches and filters the remaining matches by the phylogenetic group that the matrices belong to. It then identifies potential transcriptional modules, and generates figures that highlight such modules, taking evolution into consideration. Module formation, scoring by evolutionary criteria and visual clues reduce the amount of predictions to a manageable scale. Identification of transcription factor binding sites of particular functional importance is left to expert filtering. ReXSpecies 2 interacts with genome browsers to enable scientists to filter predictions together with other sequence-related data.

Conclusions

Based on ReXSpecies 2, we derive plausible hypotheses about the regulation of pluripotency. Our tool is designed to analyze transcription factor binding site predictions considering their common pattern of occurrence, highlighting their evolutionary history.

Antigenic and Genotypic Similarity between Primary Glioblastomas and Their Derived Neurospheres

Formation of neurospheres (NS) in cultures of glioblastomas (GBMs), with self-renewal, clonogenic capacities, and tumorigenicity following transplantation into immunodeficient mice, may denounce the existence of brain tumor stem cells (BTSCs) in vivo. In sixteen cell lines from resected primary glioblastomas, NS showed the same genetic alterations as primary tumors and the expression of stemness antigens. Adherent cells (AC), after adding 10% of fetal bovine serum (FBS) to the culture, were genetically different from NS and prevailingly expressed differentiation antigens. NS developed from a highly malignant tumor phenotype with proliferation, circumscribed necrosis, and high vessel density. Beside originating from transformed neural stem cells (NSCs), BTSCs may be contained within or correspond to dedifferentiated cells after mutation accumulation, which reacquire the expression of stemness antigens.

Learning biomarkers of pluripotent stem cells in mouse

Pluripotent stem cells are able to self-renew, and to differentiate into all adult cell types. Many studies report data describing these cells, and characterize them in molecular terms. Machine learning yields classifiers that can accurately identify pluripotent stem cells, but there is a lack of studies yielding minimal sets of best biomarkers (genes/features). We assembled gene expression data of pluripotent stem cells and non-pluripotent cells from the mouse. After normalization and filtering, we applied machine learning, classifying samples into pluripotent and non-pluripotent with high cross-validated accuracy. Furthermore, to identify minimal sets of best biomarkers, we used three methods: information gain, random forests and a wrapper of genetic algorithm and support vector machine (GA/SVM). We demonstrate that the GA/SVM biomarkers work best in combination with each other; pathway and enrichment analyses show that they cover the widest variety of processes implicated in pluripotency. The GA/SVM wrapper yields best biomarkers, no matter which classification method is used. The consensus best biomarker based on the three methods is Tet1, implicated in pluripotency just recently. The best biomarker based on the GA/SVM wrapper approach alone is Fam134b, possibly a missing link between pluripotency and some standard surface markers of unknown function processed by the Golgi apparatus.

Tbx6-dependent Sox2 regulation determines neural or mesodermal fate in axial stem cells

The classical view of neural plate development held that it arises from the ectoderm, after its separation from the mesodermal and endodermal lineages. However, recent cell lineage tracing experiments indicate that the caudal neural plate and paraxial mesoderm are generated from common bipotential axial stem cells originating from the caudal lateral epiblast (CLE)1,2. Tbx6 null mutant mouse embryos which produce ectopic neural tubes at the expense of paraxial mesoderm3 must provide a clue to the regulatory mechanism underlying this neural versus mesodermal fate choice. Here we demonstrate that Tbx6-dependent regulation of Sox2 determines the fate of axial stem cells. In wild-type embryos, enhancer N1 of the neural primordial gene Sox2 is activated in the CLE, and the cells staying in the superficial layer sustain N1 activity and activate Sox2 expression in the neural plate4-6. In contrast, the cells destined to become mesoderm activate Tbx6 and turn off enhancer N1 before migrating into the paraxial mesoderm compartment. In Tbx6 mutant embryos, however, enhancer N1 activity persists in the paraxial mesoderm compartment, eliciting ectopic Sox2 activation and transforming the paraxial mesoderm into neural tubes. An enhancer N1-specific deletion mutation introduced into Tbx6 mutant embryos prevented this Sox2 activation in the mesodermal compartment and subsequent development of ectopic neural tubes, indicating that Tbx6 regulates Sox2 via enhancer N1. Tbx6-dependent repression of Wnt3a in the paraxial mesodermal compartment is implicated in this regulatory process. Paraxial mesoderm-specific misexpression of a Sox2 transgene in wild type embryos resulted in ectopic neural tube development. Thus, Tbx6 represses Sox2 by inactivating enhancer N1 to inhibit neural development, and this is an essential step for the specification of paraxial mesoderm from the axial stem cells.

Evolution of gene regulation of pluripotency--the case for wiki tracks at genome browsers

Background

Experimentally validated data on gene regulation are hard to obtain. In particular, information about transcription factor binding sites in regulatory regions are scattered around in the literature. This impedes their systematic in-context analysis, e.g. the inference of their conservation in evolutionary history.

Results

We demonstrate the power of integrative bioinformatics by including curated transcription factor binding site information into the UCSC genome browser, using wiki and custom tracks, which enable easy publication of annotation data. Data integration allows to investigate the evolution of gene regulation of the pluripotency-associated genes Oct4, Sox2 and Nanog. For the first time, experimentally validated transcription factor binding sites in the regulatory regions of all three genes were assembled together based on manual curation of data from 39 publications. Using the UCSC genome browser, these data were then visualized in the context of multi-species conservation based on genomic alignment. We confirm previous hypotheses regarding the evolutionary age of specific regulatory patterns, establishing their "deep homology". We also confirm some other principles of Carroll's "Genetic theory of Morphological Evolution", such as "mosaic pleiotropy", exemplified by the dual role of Sox2 reflected in its regulatory region.

Conclusions

We were able to elucidate some aspects of the evolution of gene regulation for three genes associated with pluripotency. Based on the expected return on investment for the community, we encourage other scientists to contribute experimental data on gene regulation (original work as well as data collected for reviews) to the UCSC system, to enable studies of the evolution of gene regulation on a large scale, and to report their findings.

Reviewers

This article was reviewed by Dr. Gustavo Glusman and Dr. Juan Caballero, Institute for Systems Biology, Seattle, USA (nominated by Dr. Doron Lancet, Department of Molecular Genetics, Weizmann Institute of Science, Rehovot, Israel), Dr. Niels Grabe, TIGA Center (BIOQUANT) and Medical Systems Biology Group, Institute of Medical Biometry and Informatics, University Hospital Heidelberg, Germany (nominated by Dr. Mikhail Gelfand, Department of Bioinformatics, Institute of Information Transfer Problems, Russian Academy of Science, Moscow, Russian Federation) and Dr. Franz-Josef Müller, Center for Regenerative Medicine, The Scripps Research Institute, La Jolla, CA, USA and University Hospital for Psychiatry and Psychotherapy (part of ZIP gGmbH), University of Kiel, Germany (nominated by Dr. Trey Ideker, University of California, San Diego, La Jolla CA, United States).

The Pou5f1/Pou3f-dependent but SoxB-independent regulation of conserved enhancer N2 initiates Sox2 expression during epiblast to neural plate stages in vertebrates

The transcription factor Sox2 is a core component of the pluripotency control circuits in the early embryo, and later controls many aspects of neural development. Here, we demonstrate that Sox2 expression in the epiblast (mouse blastoderm) and anterior neural plate (ANP) is determined by the upstream enhancer N2. The mouse enhancer N2 exhibits strong activity in mouse ES cells, epiblast and ANP, and is regulated correctly in chicken and zebrafish embryos. Targeted deletion of this enhancer in mouse embryos caused a large reduction of Sox2 expression to 10% of that of wild-type levels in epiblast and ANP. However, this was tolerated by mouse embryo, probably due to functional compensation by Sox3. The activity of enhancer N2 depends on phylogenetically conserved bipartite POU factor-binding motifs in a 73-bp core sequence that function synergistically, but this activation does not involve Sox2. The major POU factor expressed at the epiblastic stage is Pou5f1 (Oct3/4), while those in the anterior neural plate are Pou3f factors (Oct6, Brn2 etc.). These factors are gradually exchanged during the transition from epiblast to ANP stages in mouse embryos and epiblast stem cells (EpiSC). Consistently, enhancer N2 activity changes from full Pou5f1 dependence to Pou3f dependence during the development of neural plate cells (NPC) from EpiSC, as assessed by specific POU factor knockdown in these cells. Zebrafish mutant embryos completely devoid of Pou5f1 activity failed to activate enhancer N2 and to express Sox2 in the blastoderm and ANP, and these defects were rescued by exogenous supply of pou5f1. Previously, Pou5f1-Sox2 synergism-dependent Sox2 activation through enhancer SRR2 in ES cells has been highlighted, but this mechanism is limited to ES cells and amniotes. In contrast, the enhancer N2-mediated, POU factor-dependent activation of Sox2, without involvement of Sox2, is a phylogenetically conserved core mechanism that functions in gene regulatory networks at early embryonic stages.

Evolutionarily conserved, growth plate zone-specific regulation of the matrilin-1 promoter: L-Sox5/Sox6 and Nfi factors bound near TATA finely tune activation by Sox9

To help uncover the mechanisms underlying the staggered expression of cartilage-specific genes in the growth plate, we dissected the transcriptional mechanisms driving expression of the matrilin-1 gene (Matn1). We show that a unique assembly of evolutionarily conserved cis-acting elements in the Matn1 proximal promoter restricts expression to the proliferative and prehypertrophic zones of the growth plate. These elements functionally interact with distal elements and likewise are capable of restricting the domain of activity of a pancartilaginous Col2a1 enhancer. The proximal elements include a Pe1 element binding the chondrogenic L-Sox5, Sox6, and Sox9 proteins, a SI element binding Nfi proteins, and an initiator Ine element binding the Sox trio and other factors. Sox9 binding to Pe1 is indispensable for functional interaction with the distal promoter. Binding of L-Sox5/Sox6 to Ine and Nfib to SI modulates Sox9 transactivation in a protein dose-dependent manner, possibly to enhance Sox9 activity in early stages of chondrogenesis and repress it at later stages. Hence, our data suggest a novel model whereby Sox and Nfi proteins bind to conserved Matn1 proximal elements and functionally interact with each other to finely tune gene expression in specific zones of the cartilage growth plate.

Common motifs shared by conserved enhancers of Drosophila midline glial genes

Coding sequences are usually the most highly conserved sectors of DNA, but genomic regions controlling the expression pattern of certain genes can also be conserved across diverse species. In this study, we identify five enhancers capable of activating transcription in the midline glia of Drosophila melanogaster and each contains sequences conserved across at least 11 Drosophila species. In addition, the conserved sequences contain reiterated motifs for binding sites of the known midline transcriptional activators, Single-minded, Tango, Dichaete, and Pointed. To understand the molecular basis for the highly conserved genomic subregions within enhancers of the midline genes, we tested the ability of various motifs to affect midline expression, both individually and in combination, within synthetic reporter constructs. Multiple copies of the binding site for the midline regulators Single-minded and Tango can drive expression in midline cells; however, small changes to the sequences flanking this transcription factor binding site can inactivate expression in midline cells and activate expression in tracheal cells instead. For the midline genes described in this study, the highly conserved sequences appear to juxtapose positive and negative regulatory factors in a configuration that activates genes specifically in the midline glia, while maintaining them inactive in other tissues, including midline neurons and tracheal cells.

Evolution of gene regulation--on the road towards computational inferences

If fragments of DNA are transcribed (expressed), they deserve to be called (parts of) a gene. Whether transcription takes place depends on the 'gene regulatory network'. This network is defined as the complex interplay of the sequence, biochemical modifications and structure of the chromosomal DNA with the regulatory proteins/RNA (transcription factors, co-factors, regulating RNA and the transcriptional apparatus itself). Gene regulatory networks play a role in various stages of development as well as in the maintenance of the organism; in this review we will concentrate on the former. Their evolutionary reconstruction is daunting (to say the least), and bioinformatics tools are in their infancy. However, gain of understanding offers a reward beyond itself, since evolutionary considerations can enable discoveries in the first place, e.g. the computational identification of conserved transcription factor binding sites. We discuss the evolution of gene regulation in the context of the 'Genetic Theory of Morphological Evolution' as described by Carroll, identifying those parts of the theory that are relevant for bioinformatics, and their implications. We discuss the important question of how bioinformatics analysis results on the evolution of gene regulation may be validated. Finally, we briefly exemplify use of the UCSC genome browser, exploiting its pre-computed alignments to describe the evolution of gene regulation.

B1 SOX coordinate cell specification with patterning and morphogenesis in the early zebrafish embryo

The B1 SOX transcription factors SOX1/2/3/19 have been implicated in various processes of early embryogenesis. However, their regulatory functions in stages from the blastula to early neurula remain largely unknown, primarily because loss-of-function studies have not been informative to date. In our present study, we systematically knocked down the B1 sox genes in zebrafish. Only the quadruple knockdown of the four B1 sox genes sox2/3/19a/19b resulted in very severe developmental abnormalities, confirming that the B1 sox genes are functionally redundant. We characterized the sox2/3/19a/19b quadruple knockdown embryos in detail by examining the changes in gene expression through in situ hybridization, RT–PCR, and microarray analyses. Importantly, these phenotypic analyses revealed that the B1 SOX proteins regulate the following distinct processes: (1) early dorsoventral patterning by controlling bmp2b/7; (2) gastrulation movements via the regulation of pcdh18a/18b and wnt11, a non-canonical Wnt ligand gene; (3) neural differentiation by regulating the Hes-class bHLH gene her3 and the proneural-class bHLH genes neurog1 (positively) and ascl1a (negatively), and regional transcription factor genes, e.g., hesx1, zic1, and rx3; and (4) neural patterning by regulating signaling pathway genes, cyp26a1 in RA signaling, oep in Nodal signaling, shh, and mdkb. Chromatin immunoprecipitation analysis of the her3, hesx1, neurog1, pcdh18a, and cyp26a1 genes further suggests a direct regulation of these genes by B1 SOX. We also found an interesting overlap between the early phenotypes of the B1 sox quadruple knockdown embryos and the maternal-zygotic spg embryos that are devoid of pou5f1 activity. These findings indicate that the B1 SOX proteins control a wide range of developmental regulators in the early embryo through partnering in part with Pou5f1 and possibly with other factors, and suggest that the B1 sox functions are central to coordinating cell fate specification with patterning and morphogenetic processes occurring in the early embryo.

In the developing embryo, various processes such as cell fate specification, embryo patterning, and morphogenesis take place concurrently. The embryo must control gene expression in order to coordinate these processes and thereby enable the proper organization of its structures. The B1 sox transcription factor genes, exemplified by the “stem cell gene” sox2, are thought to play a key role in these embryonic processes from the blastoderm stage to the neural stage. However, the precise regulatory functions of these genes are largely unknown due to the lack of loss-of-function studies. In our current study, we took advantage of the zebrafish system and successfully depleted B1 sox activity from the early embryo using antisense knockdown technology. This approach enabled us to further uncover the regulatory functions of B1 sox in early embryos. We found that the activity of the B1 sox genes is required for the expression of a wide range of developmental regulators including transcription factors, signaling pathway components, and cell adhesion molecules. These findings suggest that the B1 sox functions are central to coordinating diverse embryonic processes, particularly those that occur during the development of the primordium of the central nervous system.

Involvement of ubiquitous and tale transcription factors, as well as liganded RXRα, in the regulation of human SOX2 gene expression in the NT2/D1 embryonal carcinoma cell line

SOX2 is a key transcription factor in embryonic development representing a universal marker of pluripotent stem cells. Based on the functional redundancy and overlapping expression patterns of SOXB1 subgroup members during development, the goal of this study has been to analyze if some aspects of regulation of expression are preserved between human SOX2 and SOX3 genes. Thus, we have tested several transcription factors previously demonstrated to play roles in controlling SOX3 gene activity for potential participation in the regulation of SOX2 gene expression in NT2/D1 cells. Here we report on the activation of SOX2 expression by ubiquitous transcription factors (NF-Y, Sp1 and MAZ), TALE family members (Pbx1 and Meis1), as well as liganded RXR?. Elucidating components involved in the regulation of SOX gene expression represent a valuable contribution in unraveling the regulatory networks operating in pluripotent embryonic cells.

Role of SoxB1 transcription factors in development

SoxB1 factors, which include Sox1, 2, and 3, share more than 90% amino acid identity in their DNA binding HMG box and participate in diverse developmental events. They are known to exert cell-type-specific functions in concert with other transcription factors on Sox factor-dependent regulatory enhancers. Due to the high degree of sequence similarity both within and outside the HMG box, SoxB1 members show almost identical biological activities. As a result, they exhibit strong functional redundancy in regions where SoxB1 members are coexpressed, such as neural stem/progenitor cells in the developing central nervous system.