Transitioning from pluripotency to differentiated cell fates is fundamental to both embryonic development and adult tissue homeostasis. Improving our understanding of this transition would facilitate our ability to manipulate pluripotent cells into tissues for therapeutic use. Here, we show that membrane voltage (Vm) regulates the exit from pluripotency and the onset of germ layer differentiation in the embryo, a process that affects both gastrulation and left-right patterning. By examining candidate genes of congenital heart disease and heterotaxy, we identify KCNH6, a member of the ether-a-go-go class of potassium channels that hyperpolarizes the Vm and thus limits the activation of voltage gated calcium channels, lowering intracellular calcium. In pluripotent embryonic cells, depletion of kcnh6 leads to membrane depolarization, elevation of intracellular calcium levels, and the maintenance of a pluripotent state at the expense of differentiation into ectodermal and myogenic lineages. Using high-resolution temporal transcriptome analysis, we identify the gene regulatory networks downstream of membrane depolarization and calcium signaling and discover that inhibition of the mTOR pathway transitions the pluripotent cell to a differentiated fate. By manipulating Vm using a suite of tools, we establish a bioelectric pathway that regulates pluripotency in vertebrates, including human embryonic stem cells.

Gene expression is tightly regulated, with many genes exhibiting cell-specific silencing when their protein product would disrupt normal cellular function1. This silencing is largely controlled by non-coding elements, and their disruption might cause human disease2. We performed gene-agnostic screening of the non-coding regions to discover new molecular causes of congenital hyperinsulinism. This identified 14 non-coding de novo variants affecting a 42-bp conserved region encompassed by a regulatory element in intron 2 of the hexokinase 1 gene (HK1). HK1 is widely expressed across all tissues except in the liver and pancreatic beta cells and is thus termed a 'disallowed gene' in these specific tissues. We demonstrated that the variants result in a loss of repression of HK1 in pancreatic beta cells, thereby causing insulin secretion and congenital hyperinsulinism. Using epigenomic data accessed from public repositories, we demonstrated that these variants reside within a regulatory region that we determine to be critical for cell-specific silencing. Importantly, this has revealed a disease mechanism for non-coding variants that cause inappropriate expression of a disallowed gene.

ABSTRACT
. The maintenance of pluripotency in mouse embryonic stem cells (ESCs) is governed by the action of an interconnected network of transcription factors. Among them, only Oct4 and Sox2 have been shown to be strictly required for the self-renewal of ESCs and pluripotency, particularly in culture conditions in which differentiation cues are chemically inhibited. Here, we report that the conjunct activity of two orphan nuclear receptors, Esrrb and Nr5a2, parallels the importance of that of Oct4 and Sox2 in naïve mouse ESCs. By occupying a large common set of regulatory elements, these two factors control the binding of Oct4, Sox2 and Nanog to DNA. Consequently, in their absence the pluripotency network collapses and the transcriptome is substantially deregulated, leading to the differentiation of ESCs. Altogether, this work identifies orphan nuclear receptors, previously thought to be performing supportive functions, as a set of core regulators of naïve pluripotency.

Here we consider RNA-Seq, used to measure global gene expression through RNA fragmentation, capture, sequencing, and subsequent computational analysis. Xenopus, with its large number of RNArich, synchronously developing, and accessible embryos, is an excellent model organism for exploiting the power of high-throughput sequencing to understand gene expression during development. Here we present a standard RNA-Seq protocol for performing two-state differential gene expression analysis (between groups of replicates of control and treated embryos) using Illumina sequencing. Samples contain multiple whole embryos, and polyadenylated mRNA is measured under relative normalization. The protocol is divided into two parts: Wet-lab processes to prepare samples for sequencing and downstream computational analysis including quality control, quantification of gene expression, and differential expression.

The access of Transcription Factors (TFs) to their cognate DNA binding motifs requires a precise control over nucleosome positioning. This is especially important following DNA replication and during mitosis, both resulting in profound changes in nucleosome organization over TF binding regions. Using mouse Embryonic Stem (ES) cells, we show that the TF CTCF displaces nucleosomes from its binding site and locally organizes large and phased nucleosomal arrays, not only in interphase steady-state but also immediately after replication and during mitosis. Correlative analyses suggest this is associated with fast gene reactivation following replication and mitosis. While regions bound by other TFs (Oct4/Sox2), display major rearrangement, the post-replication and mitotic nucleosome positioning activity of CTCF is not unique: Esrrb binding regions are also characterized by persistent nucleosome positioning. Therefore, selected TFs such as CTCF and Esrrb act as resilient TFs governing the inheritance of nucleosome positioning at regulatory regions throughout the cell-cycle.

Embryonic development yields many different cell types in response to just a few families of inductive signals. The property of signal-receiving cells that determines how they respond to inductive signals is known as competence, and it differs in different cell types. Here, we explore the ways in which maternal factors modify chromatin to specify initial competence in the frog Xenopus tropicalis. We identify early-engaged regulatory DNA sequences, and infer from them critical activators of the zygotic genome. of these, we show that the pioneering activity of the maternal pluripotency factors Pou5f3 and Sox3 determines competence for germ layer formation by extensively remodelling compacted chromatin before the onset of inductive signalling. This remodelling includes the opening and marking of thousands of regulatory elements, extensive chromatin looping, and the co-recruitment of signal-mediating transcription factors. Our work identifies significant developmental principles that inform our understanding of how pluripotent stem cells interpret inductive signals.

One of the earliest and most significant events in embryonic development is zygotic genome activation (ZGA). In several species, bulk transcription begins at the midblastula transition (MBT) when, after a certain number of cleavages, the embryo attains a particular nuclear-to-cytoplasmic (N/C) ratio, maternal repressors become sufficiently diluted, and the cell cycle slows down. Here we resolve the frog ZGA in time and space by profiling RNA polymerase II (RNAPII) engagement and its transcriptional readout. We detect a gradual increase in both the quantity and the length of RNAPII elongation before the MBT, revealing that >1,000 zygotic genes disregard the N/C timer for their activation and that the sizes of newly transcribed genes are not necessarily constrained by cell cycle duration. We also find that Wnt, Nodal, and BMP signaling together generate most of the spatiotemporal dynamics of regional ZGA, directing the formation of orthogonal body axes and proportionate germ layers.

Transcription factor networks, together with histone modifications and signalling pathways, underlie the establishment and maintenance of gene regulatory architectures associated with the molecular identity of each cell type. However, how master transcription factors individually impact the epigenomic landscape and orchestrate the behaviour of regulatory networks under different environmental constraints is only partially understood. Here, we show that the transcription factor Nanog deploys multiple distinct mechanisms to enhance embryonic stem cell self-renewal. In the presence of LIF, which fosters self-renewal, Nanog rewires the pluripotency network by promoting chromatin accessibility and binding of other pluripotency factors to thousands of enhancers. In the absence of LIF, Nanog blocks differentiation by sustaining H3K27me3, a repressive histone mark, at developmental regulators. Among those, we show that the repression of Otx2 plays a preponderant role. Our results underscore the versatility of master transcription factors, such as Nanog, to globally influence gene regulation during developmental processes.

Mitotic bookmarking transcription factors (BFs) maintain the capacity to bind to their targets during mitosis, despite major rearrangements of the chromatin. While they were thought to propagate gene regulatory information through mitosis by statically occupying their DNA targets, it has recently become clear that BFs are highly dynamic in mitotic cells. This represents both a technical and a conceptual challenge to study and understand the function of BFs: First, formaldehyde has been suggested to be unable to efficiently capture these transient interactions, leading to profound contradictions in the literature; and second, if BFs are not permanently bound to their targets during mitosis, it becomes unclear how they convey regulatory information to daughter cells. Here, comparing formaldehyde to alternative fixatives we clarify the nature of the chromosomal association of previously proposed BFs in embryonic stem cells: While ESRRB can be considered as a canonical BF that binds at selected regulatory regions in mitosis, SOX2 and POU5F1 (also known as OCT4) establish DNA sequence-independent interactions with the mitotic chromosomes, either throughout the chromosomal arms (SOX2) or at pericentromeric regions (POU5F1). Moreover, we show that ordered nucleosomal arrays are retained during mitosis at ESRRB bookmarked sites, whereas regions losing transcription factor binding display a profound loss of order. By maintaining nucleosome positioning during mitosis, ESRRB might ensure the rapid post-mitotic re-establishment of functional regulatory complexes at selected enhancers and promoters. Our results provide a mechanistic framework that reconciles dynamic mitotic binding with the transmission of gene regulatory information across cell division.

Estrogen-related receptor b (Esrrb) is part of a family of three orphan nuclear receptors with broad expression profiles and a generic function in regulating energy metabolism in mammals. However, Esrrb performs specific functions during early mouse development, in pluripotent and multipotent populations of the embryo as well as in primordial germ cells. Moreover, Esrrb also impinges upon the control of self-renewal in embryo-derived stem cells and enhances reprogramming. Here, we review the function of Esrrb with special emphasis on its role in pluripotency. Esrrb activity at crucial regulatory elements of the pluripotency network, coupled with its role as a mitotic bookmarking factor and the ability to reset cellular metabolism, might explain its potent functions in ensuring the stability of pluripotency and driving the late stages of reprogramming. Hence, we argue that Esrrb represents a key addition to the pantheon of transcription factors sustaining pluripotent stem cell identity in mice. Understanding the mechanisms governing the interplay between different estrogen-related receptors (ERRs) and their specificity of action may clarify the role these factors play during preimplantation development and in pluripotent cells in both mouse and humans.

Antisense morpholino oligomers (MOs) have been indispensable tools for developmental biologists to transiently knock down (KD) genes rather than to knock them out (KO). Here we report on the implications of genetic KO versus MO-mediated KD of the mesoderm-specifying Brachyury paralogs in the frog Xenopus tropicalis. While both KO and KD embryos fail to activate the same core gene regulatory network, resulting in virtually identical morphological defects, embryos injected with control or target MOs also show a systemic GC content-dependent immune response and many off-target splicing defects. Optimization of MO dosage and increasing incubation temperatures can mitigate, but not eliminate, these MO side effects, which are consistent with the high affinity measured between MO and off-target sequence in vitro. We conclude that while MOs can be useful to profile loss-of-function phenotypes at a molecular level, careful attention must be paid to their immunogenic and off-target side effects.

Advances in RNA sequencing technologies have led to the surprising discovery that a vast number of transcripts emanate from regions of the genome that are not part of coding genes. Although some of the smaller ncRNAs such as microRNAs have well-characterized functions, the majority of long ncRNA (lncRNA) functions remain poorly understood. Understanding the significance of lncRNAs is an important challenge facing biology today. A powerful approach to uncovering the function of lncRNAs is to explore temporal and spatial expression profiling. This may be particularly useful for classes of lncRNAs that have developmentally important roles as the expression of such lncRNAs will be expected to be both spatially and temporally regulated during development. Here, we take advantage of our ultra-high frequency (temporal) sampling of Xenopus embryos to analyze gene expression trajectories of lncRNA transcripts over the first 3 days of development. We computationally identify 5689 potential single- and multi-exon lncRNAs. These lncRNAs demonstrate clear dynamic expression patterns. A subset of them displays highly correlative temporal expression profiles with respect to those of the neighboring genes. We also identified spatially localized lncRNAs in the gastrula stage embryo. These results suggest that lncRNAs have regulatory roles during early embryonic development.

The changes imposed on the nucleus, chromatin and its regulators during mitosis lead to the dismantlement of most gene regulatory processes. However, an increasing number of transcriptional regulators are being identified as capable of binding their genomic targets during mitosis. These so-called 'mitotic bookmarking factors' encompass transcription factors and chromatin modifiers that are believed to convey gene regulatory information from mother to daughter cells. In this Primer, we review mitotic bookmarking processes in development and stem cells and discuss the interest and potential importance of this concept with regard to epigenetic regulation and cell fate transitions involving cellular proliferation.

BACKGROUND: Because ribosomes are ubiquitously required for protein production, it was long assumed that any inherited defect in ribosome manufacture would be embryonically lethal. However, several human congenital diseases have been found to be associated with mutations in ribosome biogenesis factors. Surprisingly, despite the global requirement for ribosomes, these "ribosomopathies" are characterized by distinct and tissue specific phenotypes. The reasons for such tissue proclivity in ribosomopathies remain mysterious but may include differential expression of ribosome biogenesis factors in distinct tissues. METHODS: Here we use in situ hybridization of labeled antisense mRNA probes and ultra high temporal resolution RNA-Seq data to examine and compare expression of 13 disease associated ribosome biogenesis factors at six key stages in Xenopus tropicalis development. RESULTS: Rather than being ubiquitously expressed during development, mRNAs of all examined ribosome biogenesis factors were highly enriched in specific tissues, including the cranial neural crest and ventral blood islands. Interestingly, expression of ribosome biogenesis factors demonstrates clear differences in timing, transcript number and tissue localization. CONCLUSION: Ribosome biogenesis factor expression is more spatiotemporally regulated during embryonic development than previously expected and correlates closely with many of the common ribosomopathy phenotypes. Our findings provide information on the dynamic use of ribosome production machinery components during development and advance our understanding of their roles in disease.

Transcript regulation is essential for cell function, and misregulation can lead to disease. Despite technologies to survey the transcriptome, we lack a comprehensive understanding of transcript kinetics, which limits quantitative biology. This is an acute challenge in embryonic development, where rapid changes in gene expression dictate cell fate decisions. By ultra-high-frequency sampling of Xenopus embryos and absolute normalization of sequence reads, we present smooth gene expression trajectories in absolute transcript numbers. During a developmental period approximating the first 8 weeks of human gestation, transcript kinetics vary by eight orders of magnitude. Ordering genes by expression dynamics, we find that "temporal synexpression" predicts common gene function. Remarkably, a single parameter, the characteristic timescale, can classify transcript kinetics globally and distinguish genes regulating development from those involved in cellular metabolism. Overall, our analysis provides unprecedented insight into the reorganization of maternal and embryonic transcripts and redefines our ability to perform quantitative biology.

Correct development of the vertebrate body plan requires the early definition of two asymmetric, perpendicular axes. The first axis is established during oocyte maturation, and the second is established by symmetry breaking shortly after fertilization. The physical processes generating the second asymmetric, or dorsal-ventral, axis are well understood, but the specific molecular determinants, presumed to be maternal gene products, are poorly characterized. Whilst enrichment of maternal mRNAs at the animal and vegetal poles in both the oocyte and the early embryo has been studied, little is known about the distribution of maternal mRNAs along either the dorsal-ventral or left-right axes during the early cleavage stages. Here we report an unbiased analysis of the distribution of maternal mRNA on all axes of the Xenopus tropicalis 8-cell stage embryo, based on sequencing of single blastomeres whose positions within the embryo are known. Analysis of pooled data from complete sets of blastomeres from four embryos has identified 908 mRNAs enriched in either the animal or vegetal blastomeres, of which 793 are not previously reported as enriched. In contrast, we find no evidence for asymmetric distribution along either the dorsal-ventral or left-right axes. We confirm that animal pole enrichment is on average distinctly lower than vegetal pole enrichment, and that considerable variation is found between reported enrichment levels in different studies. We use publicly available data to show that there is a significant association between genes with human disease annotation and enrichment at the animal pole. Mutations in the human ortholog of the most animally enriched novel gene, Slc35d1, are causative for Schneckenbecken dysplasia, and we show that a similar phenotype is produced by depletion of the orthologous protein in Xenopus embryos.

Environmental stimuli are known to contribute to psoriasis pathogenesis and that of other autoimmune diseases, but the mechanisms are largely unknown. Here we show that the aryl hydrocarbon receptor (AhR), a transcription factor that senses environmental stimuli, modulates pathology in psoriasis. AhR-activating ligands reduced inflammation in the lesional skin of psoriasis patients, whereas AhR antagonists increased inflammation. Similarly, AhR signaling via the endogenous ligand FICZ reduced the inflammatory response in the imiquimod-induced model of skin inflammation and AhR-deficient mice exhibited a substantial exacerbation of the disease, compared to AhR-sufficient controls. Nonhematopoietic cells, in particular keratinocytes, were responsible for this hyperinflammatory response, which involved upregulation of AP-1 family members of transcription factors. Thus, our data suggest a critical role for AhR in the regulation of inflammatory responses and open the possibility for novel therapeutic strategies in chronic inflammatory disorders. © 2014 the Authors.

The Xenopus mid-blastula transition (MBT) marks the onset of large-scale zygotic transcription, as well as an increase in cell cycle length and a loss of synchronous cell divisions. Little is known about what triggers the activation of transcription or how newly expressed genes interact with each other. Here, we use high-resolution expression profiling to identify three waves of gene activity: a post-fertilisation wave involving polyadenylation of maternal transcripts; a broad wave of zygotic transcription detectable as early as the seventh cleavage and extending beyond the MBT at the twelfth cleavage; and a shorter post-MBT wave of transcription that becomes apparent as development proceeds. Our studies have also allowed us to define a set of maternal mRNAs that are deadenylated shortly after fertilisation, and are likely to be degraded thereafter. Experimental analysis indicates that the polyadenylation of maternal transcripts is necessary for the establishment of proper levels of zygotic transcription at the MBT, and that genes activated in the second wave of expression, including Brachyury and Mixer, contribute to the regulation of genes expressed in the third. Together, our high-resolution time series and experimental studies have yielded a deeper understanding of the temporal organisation of gene regulatory networks in the early Xenopus embryo. © 2014. Published by the Company of Biologists Ltd.

Sex chromosome divergence has been documented across phylogenetically diverse species, with amphibians typically having cytologically nondiverged ("homomorphic") sex chromosomes. With an aim of further characterizing sex chromosome divergence of an amphibian, we used "RAD-tags" and Sanger sequencing to examine sex specificity and heterozygosity in the Western clawed frog Silurana tropicalis (also known as Xenopus tropicalis). Our findings based on approximately 20 million genotype calls and approximately 200 polymerase chain reaction-amplified regions across multiple male and female genomes failed to identify a substantially sized genomic region with genotypic hallmarks of sex chromosome divergence, including in regions known to be tightly linked to the sex-determining region. We also found that expression and molecular evolution of genes linked to the sex-determining region did not differ substantially from genes in other parts of the genome. This suggests that the pseudoautosomal region, where recombination occurs, comprises a large portion of the sex chromosomes of S. tropicalis. These results may in part explain why African clawed frogs have such a high incidence of polyploidization, shed light on why amphibians have a high rate of sex chromosome turnover, and raise questions about why homomorphic sex chromosomes are so prevalent in amphibians. © the Author(s) 2013.

Background: Genomic sequence assemblies are key tools for a broad range of gene function and evolutionary studies. The diploid amphibian Xenopus tropicalis plays a pivotal role in these fields due to its combination of experimental flexibility, diploid genome, and early-branching tetrapod taxonomic position, having diverged from the amniote lineage ~360 million years ago. A genome assembly and a genetic linkage map have recently been made available. Unfortunately, large gaps in the linkage map attenuate long-range integrity of the genome assembly.Results: We laser dissected the short arm of X. tropicalis chromosome 7 for next generation sequencing and computational mapping to the reference genome. This arm is of particular interest as it encodes the sex determination locus, but its genetic map contains large gaps which undermine available genome assemblies. Whole genome amplification of 15 laser-microdissected 7p arms followed by next generation sequencing yielded ~35 million reads, over four million of which uniquely mapped to the X. tropicalis genome. Our analysis placed more than 200 previously unmapped scaffolds on the analyzed chromosome arm, providing valuable low-resolution physical map information for de novo genome assembly.Conclusion: We present a new approach for improving and validating genetic maps and sequence assemblies. Whole genome amplification of 15 microdissected chromosome arms provided sufficient high-quality material for localizing previously unmapped scaffolds and genes as well as recognizing mislocalized scaffolds. © 2013 Seifertova et al.; licensee BioMed Central Ltd.

The design of effective cell replacement therapies requires detailed knowledge of how embryonic stem cells form primary tissues, such as mesoderm or neurectoderm that later become skeletal muscle or nervous system. Members of the T-box transcription factor family are key in the formation of these primary tissues, but their underlying molecular activities are poorly understood. Here, we define invivo genome-wide regulatory inputs of the T-box proteins Brachyury, Eomesodermin, and VegT, which together maintain neuromesodermal stem cells and determine their bipotential fates in frog embryos. These T-box proteins are all recruited to the same genomic recognition sites, from where they activate genes involved in stem cell maintenance and mesoderm formation while repressing neurogenic genes. Consequently, their loss causes embryos to form an oversized neural tube with no mesodermal derivatives. This collaboration between T-box family members thus ensures the continuous formation of correctly proportioned neural and mesodermal tissues in vertebrate embryos during axial elongation

A potential mechanism that allows T cells to reliably discriminate pMHC ligands involves an interplay between kinetic proofreading, negative feedback and a destruction of this negative feedback. We analyse a detailed model of these mechanisms which involves the TCR, SHP1 and ERK. We discover that the behaviour of pSHP1 negative feedback is of primary importance, and particularly the influence of a kinetic proofreading base negative feedback state on pSHP1 dynamics. The CD8 co-receptor is shown to benefit from a kinetic proofreading locking mechanism and is able to overcome pSHP1 negative influences to sensitise a T cell.