Measuring Absolute RNA Copy Numbers at High Temporal Resolution Reveals Transcriptome Kinetics in Development

Transcriptomics of dorso-ventral axis determination in Xenopus tropicalis

Amphibian embryos provide a powerful system to study early cell fate determination because their eggs are externally fertilised, large, and easy to manipulate. Ultraviolet (UV) or lithium chloride (LiCl) treatment are classic embryonic manipulations frequently used to perturb specification of the dorso-ventral (DV) axis by affecting the stability of the maternal Wnt mediator β-catenin. Such treatments result in the formation of so-called ventralised or dorsalised embryos. Although these phenotypes have been well described with respect to their morphology and some aspects of gene expression, their whole transcriptomes have never been systematically characterised and compared. Here we show that at the early gastrula stage UV-treated embryos are transcriptionally more closely related to untreated embryos than to LiCl-treated embryos. Transcriptional comparisons with dissected ventral and dorsal regions of unperturbed gastrula embryos indicate that UV and LiCl treatments indeed enrich for ventral and dorsal cells, respectively. However, these treatments also affect the balance of neural induction in the ectodermal germ layer, with LiCl stimulating pro-neural BMP inhibition and UV preferentially generating epidermis because of elevated BMP levels. Thus the transcriptomes of UV- and LiCl-treated embryos can best be described as ventro-epidermalised and dorso-neuralised. These descriptions notwithstanding, our profiling reveals several hitherto uncharacterized genes with differential expression along the DV axis. At least one of these genes, a RNF220-like ubiquitin ligase, is activated dorsally by β-catenin. Our analysis of UV/LiCl-mediated axis perturbation will enhance the mechanistic understanding of DV axis determination in vertebrates.

ADAMTS9, a member of the ADAMTS family, in Xenopus development

Extracellular matrix (ECM) remodeling by metalloproteinases is crucial during development. The ADAMTS (A Disintegrin and Metalloproteinase with Thrombospondin type I motifs) enzymes are secreted, multi-domain matrix-associated zinc metalloendopeptidases that have diverse roles in tissue morphogenesis and patho-physiological remodeling. The human family includes 19 members. In this study we identified the 19 members of the ADAMTS family in Xenopus laevis and Xenopus tropicalis. Gene identification and a phylogenetic study revealed strong conservation of the ADAMTS family and contributed to a better annotation of the Xenopus genomes. Expression of the entire ADAMTS family was studied from early stages to tadpole stages of Xenopus, and detailed analysis of ADAMTS9 revealed expression in many structures during organogenesis such as neural crest (NC) derivative tissues, the pronephros and the pancreas. Versican, a matrix component substrate of ADAMTS9 shows a similar expression pattern suggesting a role of ADAMTS9 in the remodeling of the ECM in these structures by degradation of versican.

SpatialDE: identification of spatially variable genes

Technological advances have made it possible to measure spatially resolved gene expression at high throughput. However, methods to analyze these data are not established. Here, we develop SpatialDE, a statistical test to identify genes with spatial patterns of expression variation from multiplexed imaging or spatial RNA sequencing data. SpatialDE also implements “automatic expression histology”, a spatial gene clustering approach that enables expression-based tissue histology.

Phosphorylation states change Otx2 activity for cell proliferation and patterning in the Xenopus embryo

The homeodomain transcription factor Otx2 has essential roles in head and eye formation via the negative and positive regulation of its target genes, but it remains elusive how this dual activity of Otx2 affects cellular functions. In the current study, we first demonstrated that both exogenous and endogenous Otx2 are phosphorylated at multiple sites. Using Xenopus embryos, we identified three possible cyclin-dependent kinase (Cdk) sites and one Akt site, and analyzed the biological activities of phosphomimetic (4E) and nonphosphorylatable (4A) mutants for those sites. In the neuroectoderm, the 4E but not the 4A mutant downregulated the Cdk inhibitor gene p27xic1 (cdknx) and posterior genes, and promoted cell proliferation, possibly forming a positive-feedback loop consisting of Cdk, Otx2 and p27xic1 for cell proliferation, together with anteriorization. Conversely, the 4A mutant functioned as an activator on its own and upregulated the expression of eye marker genes, resulting in enlarged eyes. Consistent with these results, the interaction of Otx2 with the corepressor Tle1 is suggested to be phosphorylation dependent. These data suggest that Otx2 orchestrates cell proliferation, anteroposterior patterning and eye formation via its phosphorylation state.

Innate Immune Response and Off-Target Mis-splicing Are Common Morpholino-Induced Side Effects in Xenopus

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.

Primordial Germ Cell Transplantation for CRISPR/Cas9-based Leapfrogging in Xenopus

The creation of mutant lines by genome editing is accelerating genetic analysis in many organisms. CRISPR/Cas9 methods have been adapted for use in the African clawed frog, Xenopus, a longstanding model organism for biomedical research. Traditional breeding schemes for creating homozygous mutant lines with CRISPR/Cas9-targeted mutagenesis have several time-consuming and laborious steps. To facilitate the creation of mutant embryos, particularly to overcome the obstacles associated with knocking out genes that are essential for embryogenesis, a new method called leapfrogging was developed. This technique leverages the robustness of Xenopus embryos to "cut and paste" embryological methods. Leapfrogging utilizes the transfer of primordial germ cells (PGCs) from efficiently-mutagenized donor embryos into PGC-ablated wildtype siblings. This method allows for the efficient mutation of essential genes by creating chimeric animals with wildtype somatic cells that carry a mutant germline. When two F0 animals carrying "leapfrog transplants" (i.e., mutant germ cells) are intercrossed, they produce homozygous, or compound heterozygous, null F1 embryos, thus saving a full generation time to obtain phenotypic data. Leapfrogging also provides a new approach for analyzing maternal effect genes, which are refractory to F0 phenotypic analysis following CRISPR/Cas9 mutagenesis. This manuscript details the method of leapfrogging, with special emphasis on how to successfully perform PGC transplantation.

A high-resolution mRNA expression time course of embryonic development in zebrafish

We have produced an mRNA expression time course of zebrafish development across 18 time points from 1 cell to 5 days post-fertilisation sampling individual and pools of embryos. Using poly(A) pulldown stranded RNA-seq and a 3′ end transcript counting method we characterise temporal expression profiles of 23,642 genes. We identify temporal and functional transcript co-variance that associates 5024 unnamed genes with distinct developmental time points. Specifically, a class of over 100 previously uncharacterised zinc finger domain containing genes, located on the long arm of chromosome 4, is expressed in a sharp peak during zygotic genome activation. In addition, the data reveal new genes and transcripts, differential use of exons and previously unidentified 3′ ends across development, new primary microRNAs and temporal divergence of gene paralogues generated in the teleost genome duplication. To make this dataset a useful baseline reference, the data can be browsed and downloaded at Expression Atlas and Ensembl.

Stochastic models of gene transcription with upstream drives: exact solution and sample path characterization

Gene transcription is a highly stochastic and dynamic process. As a result, the mRNA copy number of a given gene is heterogeneous both between cells and across time. We present a framework to model gene transcription in populations of cells with time-varying (stochastic or deterministic) transcription and degradation rates. Such rates can be understood as upstream cellular drives representing the effect of different aspects of the cellular environment. We show that the full solution of the master equation contains two components: a model-specific, upstream effective drive, which encapsulates the effect of cellular drives (e.g. entrainment, periodicity or promoter randomness) and a downstream transcriptional Poissonian part, which is common to all models. Our analytical framework treats cell-to-cell and dynamic variability consistently, unifying several approaches in the literature. We apply the obtained solution to characterize different models of experimental relevance, and to explain the influence on gene transcription of synchrony, stationarity, ergodicity, as well as the effect of time scales and other dynamic characteristics of drives. We also show how the solution can be applied to the analysis of noise sources in single-cell data, and to reduce the computational cost of stochastic simulations.

Molecular dissection of transcriptional reprogramming of steviol glycosides synthesis in leaf tissue during developmental phase transitions in Stevia rebaudiana Bert

Stevia is a natural source of commercially important steviol glycosides (SGs), which share biosynthesis route with gibberellic acids (GAs) through plastidal MEP and cytosolic MVA pathways. Ontogeny-dependent deviation in SGs biosynthesis is one of the key factor for global cultivation of Stevia, has not been studied at transcriptional level. To dissect underlying molecular mechanism, we followed a global transcriptome sequencing approach and generated more than 100 million reads. Annotation of 41,262 de novo assembled transcripts identified all the genes required for SGs and GAs biosynthesis. Differential gene expression and quantitative analysis of important pathway genes (DXS, HMGR, KA13H) and gene regulators (WRKY, MYB, NAC TFs) indicated developmental phase dependent utilization of metabolic flux between SGs and GAs synthesis. Further, identification of 124 CYPs and 45 UGTs enrich the genomic resources, and their PPI network analysis with SGs/GAs biosynthesis proteins identifies putative candidates involved in metabolic changes, as supported by their developmental phase-dependent expression. These putative targets can expedite molecular breeding and genetic engineering efforts to enhance SGs content, biomass and yield. Futuristically, the generated dataset will be a useful resource for development of functional molecular markers for diversity characterization, genome mapping and evolutionary studies in Stevia.

Bayesian Unidimensional Scaling for visualizing uncertainty in high dimensional datasets with latent ordering of observations

BACKGROUND

Detecting patterns in high-dimensional multivariate datasets is non-trivial. Clustering and dimensionality reduction techniques often help in discerning inherent structures. In biological datasets such as microbial community composition or gene expression data, observations can be generated from a continuous process, often unknown. Estimating data points' 'natural ordering' and their corresponding uncertainties can help researchers draw insights about the mechanisms involved.

RESULTS

We introduce a Bayesian Unidimensional Scaling (BUDS) technique which extracts dominant sources of variation in high dimensional datasets and produces their visual data summaries, facilitating the exploration of a hidden continuum. The method maps multivariate data points to latent one-dimensional coordinates along their underlying trajectory, and provides estimated uncertainty bounds. By statistically modeling dissimilarities and applying a DiSTATIS registration method to their posterior samples, we are able to incorporate visualizations of uncertainties in the estimated data trajectory across different regions using confidence contours for individual data points. We also illustrate the estimated overall data density across different areas by including density clouds. One-dimensional coordinates recovered by BUDS help researchers discover sample attributes or covariates that are factors driving the main variability in a dataset. We demonstrated usefulness and accuracy of BUDS on a set of published microbiome 16S and RNA-seq and roll call data.

CONCLUSIONS

Our method effectively recovers and visualizes natural orderings present in datasets. Automatic visualization tools for data exploration and analysis are available at: https://nlhuong.shinyapps.io/visTrajectory/ .

Dual control of pcdh8l/PCNS expression and function in Xenopus laevis neural crest cells by adam13/33 via the transcription factors tfap2α and arid3a

Adam13/33 is a cell surface metalloprotease critical for cranial neural crest (CNC) cell migration. It can cleave multiple substrates including itself, fibronectin, ephrinB, cadherin-11, pcdh8 and pcdh8l (this work). Cleavage of cadherin-11 produces an extracellular fragment that promotes CNC migration. In addition, the adam13 cytoplasmic domain is cleaved by gamma secretase, translocates into the nucleus and regulates multiple genes. Here, we show that adam13 interacts with the arid3a/dril1/Bright transcription factor. This interaction promotes a proteolytic cleavage of arid3a and its translocation to the nucleus where it regulates another transcription factor: tfap2α. Tfap2α in turn activates multiple genes including the protocadherin pcdh8l (PCNS). The proteolytic activity of adam13 is critical for the release of arid3a from the plasma membrane while the cytoplasmic domain appears critical for the cleavage of arid3a. In addition to this transcriptional control of pcdh8l, adam13 cleaves pcdh8l generating an extracellular fragment that also regulates cell migration.

Zygotic Genome Activation in Vertebrates

The first major developmental transition in vertebrate embryos is the maternal-to-zygotic transition (MZT) when maternal mRNAs are degraded and zygotic transcription begins. During the MZT, the embryo takes charge of gene expression to control cell differentiation and further development. This spectacular organismal transition requires nuclear reprogramming and the initiation of RNAPII at thousands of promoters. Zygotic genome activation (ZGA) is mechanistically coordinated with other embryonic events, including changes in the cell cycle, chromatin state, and nuclear-to-cytoplasmic component ratios. Here, we review progress in understanding vertebrate ZGA dynamics in frogs, fish, mice, and humans to explore differences and emphasize common features.

Identifying cell populations with scRNASeq

Single-cell RNASeq (scRNASeq) has emerged as a powerful method for quantifying the transcriptome of individual cells. However, the data from scRNASeq experiments is often both noisy and high dimensional, making the computational analysis non-trivial. Here we provide an overview of different experimental protocols and the most popular methods for facilitating the computational analysis. We focus on approaches for identifying biologically important genes, projecting data into lower dimensions and clustering data into putative cell-populations. Finally we discuss approaches to validation and biological interpretation of the identified cell-types or cell-states.

Chromatin Control of Developmental Dynamics and Plasticity

Chromatin structure is intimately connected with gene expression and cell identity. Here we review recent advances in the field and discuss how establishment of cell identity during development is accompanied by large-scale remodeling of the epigenetic landscape and how this remodeling drives and supports lineage specification and maintenance. We discuss maternal control of the early embryonic epigenetic landscape, selective usage of enhancer clusters via 3D chromatin contacts leading to activation of transcription factor networks, and conserved regulation of developmental pathways by specific DNA demethylation of key regulatory regions. Together, these processes establish an epigenetic framework regulating different phases of embryonic development.

Using Xenopus to understand human disease and developmental disorders

Model animals are crucial to biomedical research. Among the commonly used model animals, the amphibian, Xenopus, has had tremendous impact because of its unique experimental advantages, cost effectiveness, and close evolutionary relationship with mammals as a tetrapod. Over the past 50 years, the use of Xenopus has made possible many fundamental contributions to biomedicine, and it is a cornerstone of research in cell biology, developmental biology, evolutionary biology, immunology, molecular biology, neurobiology, and physiology. The prospects for Xenopus as an experimental system are excellent: Xenopus is uniquely well-suited for many contemporary approaches used to study fundamental biological and disease mechanisms. Moreover, recent advances in high throughput DNA sequencing, genome editing, proteomics, and pharmacological screening are easily applicable in Xenopus, enabling rapid functional genomics and human disease modeling at a systems level.

Competition between histone and transcription factor binding regulates the onset of transcription in zebrafish embryos

Upon fertilization, the genome of animal embryos remains transcriptionally inactive until the maternal-to-zygotic transition. At this time, the embryo takes control of its development and transcription begins. How the onset of zygotic transcription is regulated remains unclear. Here, we show that a dynamic competition for DNA binding between nucleosome-forming histones and transcription factors regulates zebrafish genome activation. Taking a quantitative approach, we found that the concentration of non-DNA-bound core histones sets the time for the onset of transcription. The reduction in nuclear histone concentration that coincides with genome activation does not affect nucleosome density on DNA, but allows transcription factors to compete successfully for DNA binding. In agreement with this, transcription factor binding is sensitive to histone levels and the concentration of transcription factors also affects the time of transcription. Our results demonstrate that the relative levels of histones and transcription factors regulate the onset of transcription in the embryo.

DOI:http://dx.doi.org/10.7554/eLife.23326.001

The DNA in a fertilized egg contains all the information required to form an animal’s body. In order for the animal to develop properly, particular genes encoded in the DNA are only active at specific times. The DNA is wrapped around proteins called histones, which allows the DNA to be tightly packed inside the cell. However, histones can block other proteins called transcription factors from binding to the DNA to activate the genes. Young embryos initially develop with all of their genes switched off, relying on the nutrients and other molecules provided by their mother. After some time, the embryo starts to switch on its own genes to take control of its own development, but it was not clear how this happens.

Joseph et al. investigated how genes are activated in zebrafish embryos, which are often used as models to study how animals develop. The experiments show that competition between histones and transcription factors for binding to DNA controls when genes are switched on. In young fish embryos, there are so many histones present that transcription factors have no opportunity to bind to DNA. Over time, however, the numbers of histones decrease, allowing transcription factors to bind to DNA and switch on genes.

Histones and transcription factors regulate the activity of genes throughout the life of the animal. Therefore, competition between these two types of protein may also control gene activity in other situations. A better understanding of how gene activity is controlled could allow researchers to more easily grow different types of cell in the laboratory or to reprogram specific cells in the body. As such, these new findings may aid the development of therapies to regenerate organs or tissues that have been damaged by injury or disease.

DOI:http://dx.doi.org/10.7554/eLife.23326.002

Power analysis of single-cell RNA-sequencing experiments

Single-cell RNA sequencing (scRNA-seq) has become an established and powerful method to investigate transcriptomic cell-to-cell variation, thereby revealing new cell types and providing insights into developmental processes and transcriptional stochasticity. A key question is how the variety of available protocols compare in terms of their ability to detect and accurately quantify gene expression. Here, we assessed the protocol sensitivity and accuracy of many published data sets, on the basis of spike-in standards and uniform data processing. For our workflow, we developed a flexible tool for counting the number of unique molecular identifiers (https://github.com/vals/umis/). We compared 15 protocols computationally and 4 protocols experimentally for batch-matched cell populations, in addition to investigating the effects of spike-in molecular degradation. Our analysis provides an integrated framework for comparing scRNA-seq protocols.

Foxh1 Occupies cis-Regulatory Modules Prior to Dynamic Transcription Factor Interactions Controlling the Mesendoderm Gene Program

The interplay between transcription factors and chromatin dictates gene regulatory network activity. Germ layer specification is tightly coupled with zygotic gene activation and, in most metazoans, is dependent upon maternal factors. We explore the dynamic genome-wide interactions of Foxh1, a maternal transcription factor that mediates Nodal/TGF-β signaling, with cis-regulatory modules (CRMs) during mesendodermal specification. Foxh1 marks CRMs during cleavage stages and recruits the co-repressor Tle/Groucho in the early blastula. We highlight a population of CRMs that are continuously occupied by Foxh1 and show that they are marked by H3K4me1, Ep300, and Fox/Sox/Smad motifs, suggesting interplay between these factors in gene regulation. We also propose a molecular "hand-off" between maternal Foxh1 and zygotic Foxa at these CRMs to maintain enhancer activation. Our findings suggest that Foxh1 functions at the top of a hierarchy of interactions by marking developmental genes for activation, beginning with the onset of zygotic gene expression.

A gene regulatory program controlling early Xenopus mesendoderm formation: Network conservation and motifs

Germ layer formation is among the earliest differentiation events in metazoan embryos. In triploblasts, three germ layers are formed, among which the endoderm gives rise to the epithelial lining of the gut tube and associated organs including the liver, pancreas and lungs. In frogs (Xenopus), where early germ layer formation has been studied extensively, the process of endoderm specification involves the interplay of dozens of transcription factors. Here, we review the interactions between these factors, summarized in a transcriptional gene regulatory network (GRN). We highlight regulatory connections conserved between frog, fish, mouse, and human endodermal lineages. Especially prominent is the conserved role and regulatory targets of the Nodal signaling pathway and the T-box transcription factors, Vegt and Eomes. Additionally, we highlight network topologies and motifs, and speculate on their possible roles in development.

Single-cell RNA-seq and computational analysis using temporal mixture modelling resolves Th1/Tfh fate bifurcation in malaria

Differentiation of naïve CD4+ T cells into functionally distinct T helper subsets is crucial for the orchestration of immune responses. Due to extensive heterogeneity and multiple overlapping transcriptional programs in differentiating T cell populations, this process has remained a challenge for systematic dissection in vivo. By using single-cell transcriptomics and computational analysis using a temporal mixtures of Gaussian processes model, termed GPfates, we reconstructed the developmental trajectories of Th1 and Tfh cells during blood-stage Plasmodium infection in mice. By tracking clonality using endogenous TCR sequences, we first demonstrated that Th1/Tfh bifurcation had occurred at both population and single-clone levels. Next, we identified genes whose expression was associated with Th1 or Tfh fates, and demonstrated a T-cell intrinsic role for Galectin-1 in supporting a Th1 differentiation. We also revealed the close molecular relationship between Th1 and IL-10-producing Tr1 cells in this infection. Th1 and Tfh fates emerged from a highly proliferative precursor that upregulated aerobic glycolysis and accelerated cell cycling as cytokine expression began. Dynamic gene expression of chemokine receptors around bifurcation predicted roles for cell-cell in driving Th1/Tfh fates. In particular, we found that precursor Th cells were coached towards a Th1 but not a Tfh fate by inflammatory monocytes. Thus, by integrating genomic and computational approaches, our study has provided two unique resources, a database www.PlasmoTH.org, which facilitates discovery of novel factors controlling Th1/Tfh fate commitment, and more generally, GPfates, a modelling framework for characterizing cell differentiation towards multiple fates.

An Enumerative Combinatorics Model for Fragmentation Patterns in RNA Sequencing Provides Insights into Nonuniformity of the Expected Fragment Starting-Point and Coverage Profile

RNA sequencing (RNA-seq) has emerged as the method of choice for measuring the expression of RNAs in a given cell population. In most RNA-seq technologies, sequencing the full length of RNA molecules requires fragmentation into smaller pieces. Unfortunately, the issue of nonuniform sequencing coverage across a genomic feature has been a concern in RNA-seq and is attributed to biases for certain fragments in RNA-seq library preparation and sequencing. To investigate the expected coverage obtained from fragmentation, we develop a simple fragmentation model that is independent of bias from the experimental method and is not specific to the transcript sequence. Essentially, we enumerate all configurations for maximal placement of a given fragment length, F, on transcript length, T, to represent every possible fragmentation pattern, from which we compute the expected coverage profile across a transcript. We extend this model to incorporate general empirical attributes such as read length, fragment length distribution, and number of molecules of the transcript. We further introduce the fragment starting-point, fragment coverage, and read coverage profiles. We find that the expected profiles are not uniform and that factors such as fragment length to transcript length ratio, read length to fragment length ratio, fragment length distribution, and number of molecules influence the variability of coverage across a transcript. Finally, we explore a potential application of the model where, with simulations, we show that it is possible to correctly estimate the transcript copy number for any transcript in the RNA-seq experiment.

Heat Shock Proteins and Maternal Contribution to Oogenesis and Early Embryogenesis

Early embryos develop from fertilized eggs using materials that are stored during oocyte growth and which can be defined as maternal contribution (molecules, factors, or determinants). Several heat shock proteins (HSPs) and the heat shock transcriptional factor (HSF) are part of the maternal contribution that is critical for successful embryogenesis and reproduction. A maternal role for heat shock-related genes was mainly demonstrated in genetic experimental organisms (e.g., fly, nematode, mouse). Nowadays, an increasing number of "omics" data are produced from a large panel of organisms implementing a catalog of maternal and/or embryonic HSPs and HSFs. However, for most of them, it remains to better understand their potential roles in this context. Existing and future genome-wide screens mainly set up to create loss-of-function are likely to improve this situation. This chapter will discuss available data from various experimental organisms following the developmental steps from egg production to early embryogenesis.

Expression of ribosomopathy genes during Xenopus tropicalis embryogenesis

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12861-016-0138-5) contains supplementary material, which is available to authorized users.

MicroRNAs and ectodermal specification I. Identification of miRs and miR-targeted mRNAs in early anterior neural and epidermal ectoderm

The establishment of cell lineages occurs via a dynamic progression of gene regulatory networks (GRNs) that underlie developmental commitment and differentiation. To investigate how microRNAs (miRs) function in this process, we compared miRs and miR targets at the initiation of the two major ectodermal lineages in Xenopus. We used next-generation sequencing to identify over 170 miRs expressed in midgastrula ectoderm expressing either noggin or a constitutively active BMP receptor, reflecting anterior neural or epidermal ectoderm, respectively; 125 had not previously been identified in Xenopus. We identified the locations of the pre-miR sequences in the X. laevis genome. Neural and epidermal ectoderm express broadly similar sets of miRs. To identify targets of miR-dependent translational control, we co-immunoprecipitated Argonaute-Ribonucleoprotein (Ago-RNP) complexes from early neural and epidermal ectoderm and sequenced the associated RNA. The Ago-RNP RNAs from these tissues represent overlapping, yet distinct, subsets of genes. Moreover, the profile of Ago-RNP associated genes differs substantially from the profile of total RNAs in these tissues. We generated target predictions for the "high confidence" Ago-RNP RNAs using the identified ectodermal miRs; These RNAs generally had target sites for multiple miRs. Oct4 orthologues, as well as many of their previously identified transcriptional targets, are represented in the Ago-RNP pool in both tissues, suggesting that miR-dependent regulation contributes to the downregulation of the oct4 gene regulatory network and the reduction in ectodermal pluripotency.

Leapfrogging: primordial germ cell transplantation permits recovery of CRISPR/Cas9-induced mutations in essential genes

CRISPR/Cas9 genome editing is revolutionizing genetic loss-of-function analysis but technical limitations remain that slow progress when creating mutant lines. First, in conventional genetic breeding schemes, mosaic founder animals carrying mutant alleles are outcrossed to produce F1 heterozygotes. Phenotypic analysis occurs in the F2 generation following F1 intercrosses. Thus, mutant analyses will require multi-generational studies. Second, when targeting essential genes, efficient mutagenesis of founders is often lethal, preventing the acquisition of mature animals. Reducing mutagenesis levels may improve founder survival, but results in lower, more variable rates of germline transmission. Therefore, an efficient approach to study lethal mutations would be useful. To overcome these shortfalls, we introduce 'leapfrogging', a method combining efficient CRISPR mutagenesis with transplantation of mutated primordial germ cells into a wild-type host. Tested using Xenopus tropicalis, we show that founders containing transplants transmit mutant alleles with high efficiency. F1 offspring from intercrosses between F0 animals that carry embryonic lethal alleles recapitulate loss-of-function phenotypes, circumventing an entire generation of breeding. We anticipate that leapfrogging will be transferable to other species.

A catalog of Xenopus tropicalis transcription factors and their regional expression in the early gastrula stage embryo

Gene regulatory networks (GRNs) involve highly combinatorial interactions between transcription factors and short sequence motifs in cis-regulatory modules of target genes to control cellular phenotypes. The GRNs specifying most cell types are largely unknown and are the subject of wide interest. A catalog of transcription factors is a valuable tool toward obtaining a deeper understanding of the role of these critical effectors in any biological setting. Here we present a comprehensive catalog of the transcription factors for the diploid frog Xenopus tropicalis. We identify 1235 genes encoding DNA-binding transcription factors, comparable to the numbers found in typical mammalian species. In detail, the repertoire of X. tropicalis transcription factor genes is nearly identical to human and mouse, with the exception of zinc finger family members, and a small number of species/lineage-specific gene duplications and losses relative to the mammalian repertoires. We applied this resource to the identification of transcription factors differentially expressed in the early gastrula stage embryo. We find transcription factor enrichment in Spemann's organizer, the ventral mesoderm, ectoderm and endoderm, and report 218 TFs that show regionalized expression patterns at this stage. Many of these have not been previously reported as expressed in the early embryo, suggesting thus far unappreciated roles for many transcription factors in the GRNs regulating early development. We expect our transcription factor catalog will facilitate myriad studies using Xenopus as a model system to understand basic biology and human disease.

Analysis of embryonic development in the unsequenced axolotl: Waves of transcriptomic upheaval and stability

The axolotl (Ambystoma mexicanum) has long been the subject of biological research, primarily owing to its outstanding regenerative capabilities. However, the gene expression programs governing its embryonic development are particularly underexplored, especially when compared to other amphibian model species. Therefore, we performed whole transcriptome polyA+ RNA sequencing experiments on 17 stages of embryonic development. As the axolotl genome is unsequenced and its gene annotation is incomplete, we built de novo transcriptome assemblies for each stage and garnered functional annotation by comparing expressed contigs with known genes in other organisms. In evaluating the number of differentially expressed genes over time, we identify three waves of substantial transcriptome upheaval each followed by a period of relative transcriptome stability. The first wave of upheaval is between the one and two cell stage. We show that the number of differentially expressed genes per unit time is higher between the one and two cell stage than it is across the mid-blastula transition (MBT), the period of zygotic genome activation. We use total RNA sequencing to demonstrate that the vast majority of genes with increasing polyA+ signal between the one and two cell stage result from polyadenylation rather than de novo transcription. The first stable phase begins after the two cell stage and continues until the mid-blastula transition, corresponding with the pre-MBT phase of transcriptional quiescence in amphibian development. Following this is a peak of differential gene expression corresponding with the activation of the zygotic genome and a phase of transcriptomic stability from stages 9-11. We observe a third wave of transcriptomic change between stages 11 and 14, followed by a final stable period. The last two stable phases have not been documented in amphibians previously and correspond to times of major morphogenic change in the axolotl embryo: gastrulation and neurulation. These results yield new insights into global gene expression during early stages of amphibian embryogenesis and will help to further develop the axolotl as a model species for developmental and regenerative biology.

Developmentally regulated long non-coding RNAs in Xenopus tropicalis

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.

Xenopus genomic data and browser resources

The two species of Xenopus most commonly used in biomedical research are the diploid Xenopus (Silurana) tropicalis and the tetraploid Xenopus laevis. The X. tropicalis genome sequence has been available since 2010 and this year the X. laevis, genome from two distinct genetic backgrounds has been published. Multiple genome assemblies available for both species and transcriptomic and epigenetic data sets are growing rapidly, all of which are available from a variety of web resources. This review describes the contents of these resources, how to locate and download genomic data, and also how to view and manipulate these data on various public genome browsers, with an emphasis on Xenbase, the Xenopus model organism database.