Developmental gene regulatory networks in the zebrafish embryo

Evolutionarily conserved Wnt/Sp5 signaling is critical for anterior-posterior axis patterning in sea urchin embryos

Studies across a diverse group of metazoan embryos indicate that Wnt signaling often activates the transcription factor Sp5, forming a signaling 'cassette' that plays critical roles in many developmental processes. This study explores the role of Wnt/Sp5 signaling during the specification and patterning of the primary germ layers during early anterior-posterior axis formation in the deuterostome sea urchin embryo. Our functional analyses show that Sp5 is critical for endomesoderm specification downstream of Wnt/β-catenin in posterior cells as well as anterior neuroectoderm patterning downstream of non-canonical Wnt/JNK signaling in anterior cells. Interestingly, expression and functional data comparisons show that Wnt/Sp5 signaling often plays similar roles in posterior endomesoderm as well as neuroectoderm patterning along the AP axis of several deuterostome embryos, including vertebrates. Thus, our findings provide strong support for the idea that Wnt-Sp5 signaling cassettes were critical for the establishment of early germ layers in the common deuterostome ancestor.

Tetrabromobisphenol a exacerbates the overall radioactive hazard to zebrafish (Danio rerio)

The major health risks of dual exposure to two hazardous factors of plastics and radioactive contamination are obscure. In the present study, we systematically evaluated the combinational toxic effects of tetrabromobisphenol A (TBBPA), one of the most influential plastic ingredients, mainly from electronic wastes, and γ-irradiation in zebrafish for the first time. TBBPA (0.25 μg/mL for embryos and larvae, 300 μg/L for adults) contamination aggravated the radiation (6 Gy for embryos and larvae, 20 Gy for adults)-induced early dysplasia and aberrant angiogenesis of embryos, further impaired the locomotor vitality of irradiated larvae, and worsened the radioactive multiorganic histologic injury, neurobehavioural disturbances and dysgenesis of zebrafish adults as well as the inter-generational neurotoxicity in offspring. TBBPA exaggerated the radiative toxic effects not only by enhancing the inflammatory and apoptotic response but also by further unbalancing the endocrine system and disrupting the underlying gene expression profiles. In conclusion, TBBPA exacerbates radiation-induced injury in zebrafish, including embryos, larvae, adults and even the next generation. Our findings provide new insights into the toxicology of TBBPA and γ-irradiation, shedding light on the severity of cocontamination of MP components and radioactive substances and thereby inspiring novel remediation and rehabilitation strategies for radiation-injured aqueous organisms and radiotherapy patients.

The fate and oxidative stress of different sized SiO2 nanoparticles in zebrafish (Danio rerio) larvae

Nanoparticle such as silicon dioxide nanoparticles (nano-SiO₂) are extensively produced and applied in society. Hence there is an increasing concern about their exposure and toxicity to human and wildlife. To understand the effects of sizes of NPs on their bioavailability and toxicity, zebrafish (Danio rerio) embryos (2 h post-fertilization, hpf) were exposed to 25, 50 and 100 mg/L of 15 or 30 nm nano-SiO₂ for 5 days respectively. The results showed that SiO₂ could be readily uptaken by zebrafish, and the accumulation of SiO₂ was significantly higher in 15 nm treatments groups compared to 30 nm nano-SiO₂ treated groups. Furthermore, exposure to 15 nm nano-SiO₂ at the concentration of 100 mg/L resulted in more significant changes in reactive oxygen species (ROS) levels, perturbation of lipid peroxidative and antioxidant system than the same concentration of 30 nm nano-SiO₂, indicating small sized nano-SiO₂ evoked severer oxidative stress in zebrafish larvae. In addition, the more significant up-regulation of transcription of genes related to oxidative stress (e.g., nrf2 and sod1) in 15 nm nano-SiO₂ at the 100 mg/L treatments groups provided more evidence for this speculation. Given the above, 15 nm nano-SiO₂ were more likely to enter and accumulate in zebrafish larvae, thus causing more serious oxidative stress in vivo. These results may provide additional information on the fate and toxicities of different sizes of NPs.

Developmental regulation of Wnt signaling by Nagk and the UDP-GlcNAc salvage pathway

In a screen for human kinases that regulate Xenopus laevis embryogenesis, we identified Nagk and other components of the UDP-GlcNAc glycosylation salvage pathway as regulators of anteroposterior patterning and Wnt signaling. We find that the salvage pathway does not affect other major embryonic signaling pathways (Fgf, TGFβ, Notch, or Shh), thereby demonstrating specificity for Wnt signaling. We show that the role of the salvage pathway in Wnt signaling is evolutionarily conserved in zebrafish and Drosophila. Finally, we show that GlcNAc is essential for the growth of intestinal enteroids, which are highly dependent on Wnt signaling for growth and maintenance. We propose that the Wnt pathway is sensitive to alterations in the glycosylation state of a cell and acts as a nutritional sensor in order to couple growth/proliferation with its metabolic status. We also propose that the clinical manifestations observed in congenital disorders of glycosylation (CDG) in humans may be due, in part, to their effects on Wnt signaling during development.

Methods for the experimental and computational analysis of gene regulatory networks in sea urchins

The discovery of gene regulatory networks (GRNs) has opened a gate to access the genomic mechanisms controlling development. GRNs are systems of transcriptional regulatory circuits that control the differential specification of cell fates during development by regulating gene expression. The experimental analysis of GRNs involves a collection of methods, each revealing aspects of the overall control process. This review provides an overview of experimental and computational methods that have been successfully applied for solving developmental GRNs in the sea urchin embryo. The key in this approach is to obtain experimental evidence for functional interactions between transcription factors and regulatory DNA. In the second part of this review, a more generally applicable strategy is discussed that shows a path from experimental evidence to annotation of regulatory linkages to the generation of GRN models.

Nodal signaling has dual roles in fate specification and directed migration during germ layer segregation in zebrafish

During gastrulation, endodermal cells actively migrate to the interior of the embryo, but the signals that initiate and coordinate this migration are poorly understood. By transplanting ectopically induced endodermal cells far from the normal location of endoderm specification, we identified the inputs that drive internalization without the confounding influences of fate specification and global morphogenic movements. We find that Nodal signaling triggers an autocrine circuit for initiating endodermal internalization. Activation of the Nodal receptor directs endodermal specification through sox32 and also induces expression of more Nodal ligands. These ligands act in an autocrine fashion to initiate endodermal cell sorting. Our work defines an 'AND' gate consisting of sox32-dependent endodermal specification and Nodal ligand reception controlling endodermal cell sorting to the inner layer of the embryo at the onset of gastrulation.

Analysis of novel caudal hindbrain genes reveals different regulatory logic for gene expression in rhombomere 4 versus 5/6 in embryonic zebrafish

BACKGROUND

Previous work aimed at understanding the gene regulatory networks (GRNs) governing caudal hindbrain formation identified morphogens such as Retinoic Acid (RA) and Fibroblast growth factors (FGFs), as well as transcription factors like hoxb1b, hoxb1a, hnf1ba, and valentino as being required for rhombomere (r) r4-r6 formation in zebrafish. Considering that the caudal hindbrain is relatively complex - for instance, unique sets of neurons are formed in each rhombomere segment - it is likely that additional essential genes remain to be identified and integrated into the caudal hindbrain GRN.

METHODS

By taking advantage of gene expression data available in the Zebrafish Information Network (ZFIN), we identified 84 uncharacterized genes that are expressed in r4-r6. We selected a representative set of 22 genes and assayed their expression patterns in hoxb1b, hoxb1a, hnf1b, and valentino mutants with the goal of positioning them in the caudal hindbrain GRN. We also investigated the effects of RA and FGF on the expression of this gene set. To examine whether these genes are necessary for r4-r6 development, we analyzed germline mutants for six of the genes (gas6, gbx1, sall4, eglf6, celf2, and greb1l) for defects in hindbrain development.

RESULTS

Our results reveal that r4 gene expression is unaffected by the individual loss of hoxb1b, hoxb1a or RA, but is under the combinatorial regulation of RA together with hoxb1b. In contrast, r5/r6 gene expression is dependent on RA, FGF, hnf1ba and valentino - as individual loss of these factors abolishes r5/r6 gene expression. Our analysis of six mutant lines did not reveal rhombomere or neuronal defects, but transcriptome analysis of one line (gas6 mutant) identified expression changes for genes involved in several developmental processes - suggesting that these genes may have subtle roles in hindbrain development.

CONCLUSION

We conclude that r4-r6 formation is relatively robust, such that very few genes are absolutely required for this process. However, there are mechanistic differences in r4 versus r5/r6, such that no single factor is required for r4 development while several genes are individually required for r5/r6 formation.

Tetrabromobisphenol A caused neurodevelopmental toxicity via disrupting thyroid hormones in zebrafish larvae

Tetrabromobisphenol A (TBBPA), one of the most widely used brominated flame retardants (BFRs), has resulted in worldwide environmental contamination. TBBPA has been reported as a thyroid endocrine disruptor and a potential neurotoxicant. However, the underlying mechanism is still not clear. In this study, zebrafish (Danio rerio) embryos (2 h post-fertilization, hpf) were exposed to different concentrations of TBBPA (50, 100, 200 and 400 μg/L) alone or in combination with 3,3',5-triiodo-l-thyronine (T3, 20 μg/L + TBBPA, 200 μg/L). The results confirmed that TBBPA could evoke thyroid disruption by observations of increased T4 contents and decreased T3 contents, accompanied by up-regulated tshβ, tg mRNA and down-regulated ttr and trβ mRNA levels in zebafish larvae. TBBPA-induced neurodevelopmental toxicity was also indicated by down-regulated transcription of genes related to central nervous system (CNS) development (e.g., α1-tubulin, mbp and shha), and decreased locomotor activity and average swimming speed. Our results further demonstrated that treatment with T3 could reverse or eliminate TBBPA-induced effects on thyroidal and neurodevelopmental parameters. Given the above, we hypothesize that the observed neurodevelopmental toxicity in the present study could be attributed to the thyroid hormone disruptions by TBBPA.

The maternal control in the embryonic development of zebrafish

The maternal control directing the very first hours of life is of pivotal importance for ensuring proper development to the growing embryo. Thanks to the finely regulated inheritance of maternal factors including mRNAs and proteins produced during oogenesis and stored into the mature oocyte, the embryo is sustained throughout the so-called maternal-to-zygotic transition, a period in development characterized by a species-specific length in time, during which critical biological changes regarding cell cycle and zygotic transcriptional activation occur. In order not to provoke any kind of persistent damage, the process must be delicately balanced. Surprisingly, our knowledge as to the possible effects of beneficial bacteria regarding the modulation of the quality and/or quantity of both maternally-supplied and zygotically-transcribed mRNAs, is very limited. To date, only one group has investigated the consequences of the parentally-supplied Lactobacillus rhamnosus on the storage of mRNAs into mature oocytes, leading to an altered maternal control process in the F1 generation. Particular attention was called on the monitoring of several biomarkers involved in autophagy, apoptosis and axis patterning, while data on miRNA generation and pluripotency maintenance are herein presented for the first time, and can assist in laying the ground for further investigations in this field. In this review, the reader is supplied with the current knowledge on the above-mentioned biological process, first by drawing the general background and then by emphasizing the most important findings that have highlighted their focal role in normal animal development.

Modulation of N-Methyl-D-Aspartate Receptors (NMDAR), Bcl-2 and C-Fos Gene Expressions on Exposure to Individual and Mixtures of Low Concentration Metals in Zebrafish (Danio rerio)

Currently, there is limited information on the toxicity of low concentration of metal mixtures in the environment. Of particular interest is the effect of low levels of metal mixtures on neurodevelopment of aquatic organisms. This study reports the neurological gene expressions after exposing zebrafish embryos to low concentration toxic heavy metals, 120 h post fertilization (hpf). Embryos were exposed to low concentration individual and mixtures of lead (Pb), mercury (Hg), arsenic (As), and cadmium (Cd). Quantitative real-time PCR was used to assess gene expressions. The findings of this study confirmed that exposure to low concentration heavy metals upregulated N-methyl-D-aspartate (NMDA) receptor subunits NMDAR2A (NR2A), NMDAR2B (NR2B), and NMDAR2D (NR2D) and B cell lymphoma (Bcl-2) genes. NR2A genes were significantly upregulated by 90 and 74%, respectively, on exposure to Pb + As and Pb + Cd. NR2B genes were upregulated by 85.3, 68.6, 62.7, and 62.7% on exposure to As, Pb + Hg, Pb + As, and Pb + Cd, respectively. Exposure to As, Pb + Cd, and Pb + Hg + As significantly upregulated Bcl-2 genes by 2.01-, 1.84-, and 1.80-fold, respectively. NR1A and C-fos gene expressions were not significantly different from control. Upregulation of NMDAR subunits and Bcl-2 genes in this study was largely a counter measure against insults from exposure to low concentration heavy metals. Principal component analysis confirmed the influence of low concentration individual and mixtures of Pb, Hg, As, and Cd on gene expression of NMDAR subunits and Bcl-2. These data suggest that altered expression of NMDA receptor subunits and Bcl-2 genes may explain toxicity of low concentration individual and mixtures of Pb, Hg, As, and Cd.

Deployment of regulatory genes during gastrulation and germ layer specification in a model spiralian mollusc Crepidula

BACKGROUND

During gastrulation, endoderm and mesoderm are specified from a bipotential precursor (endomesoderm) that is argued to be homologous across bilaterians. Spiralians also generate mesoderm from ectodermal precursors (ectomesoderm), which arises near the blastopore. While a conserved gene regulatory network controls specification of endomesoderm in deuterostomes and ecdysozoans, little is known about genes controlling specification or behavior of either source of spiralian mesoderm or the digestive tract.

RESULTS

Using the mollusc Crepidula, we examined conserved regulatory factors and compared their expression to fate maps to score expression in the germ layers, blastopore lip, and digestive tract. Many genes were expressed in both ecto- and endomesoderm, but only five were expressed in ectomesoderm exclusively. The latter may contribute to epithelial-to-mesenchymal transition seen in ectomesoderm.

CONCLUSIONS

We present the first comparison of genes expressed during spiralian gastrulation in the context of high-resolution fate maps. We found variation of genes expressed in the blastopore lip, mouth, and cells that will form the anus. Shared expression of many genes in both mesodermal sources suggests that components of the conserved endomesoderm program were either co-opted for ectomesoderm formation or that ecto- and endomesoderm are derived from a common mesodermal precursor that became subdivided into distinct domains during evolution.

Understanding transcriptional regulatory networks using computational models

Transcriptional regulatory networks (TRNs) encode instructions for animal development and physiological responses. Recent advances in genomic technologies and computational modeling have revolutionized our ability to construct models of TRNs. Here, we survey current computational methods for inferring TRN models using genome-scale data. We discuss their advantages and limitations. We summarize representative TRNs constructed using genome-scale data in both normal and disease development. We discuss lessons learned about the structure/function relationship of TRNs, based on examining various large-scale TRN models. Finally, we outline some open questions regarding TRNs, including how to improve model accuracy by integrating complementary data types, how to infer condition-specific TRNs, and how to compare TRNs across conditions and species in order to understand their structure/function relationship.

BioTapestry now provides a web application and improved drawing and layout tools

Gene regulatory networks (GRNs) control embryonic development, and to understand this process in depth, researchers need to have a detailed understanding of both the network architecture and its dynamic evolution over time and space. Interactive visualization tools better enable researchers to conceptualize, understand, and share GRN models. BioTapestry is an established application designed to fill this role, and recent enhancements released in Versions 6 and 7 have targeted two major facets of the program. First, we introduced significant improvements for network drawing and automatic layout that have now made it much easier for the user to create larger, more organized network drawings. Second, we revised the program architecture so it could continue to support the current Java desktop Editor program, while introducing a new BioTapestry GRN Viewer that runs as a JavaScript web application in a browser. We have deployed a number of GRN models using this new web application. These improvements will ensure that BioTapestry remains viable as a research tool in the face of the continuing evolution of web technologies, and as our understanding of GRN models grows.

Shaped 3D singular spectrum analysis for quantifying gene expression, with application to the early zebrafish embryo

Recent progress in microscopy technologies, biological markers, and automated processing methods is making possible the development of gene expression atlases at cellular-level resolution over whole embryos. Raw data on gene expression is usually very noisy. This noise comes from both experimental (technical/methodological) and true biological sources (from stochastic biochemical processes). In addition, the cells or nuclei being imaged are irregularly arranged in 3D space. This makes the processing, extraction, and study of expression signals and intrinsic biological noise a serious challenge for 3D data, requiring new computational approaches. Here, we present a new approach for studying gene expression in nuclei located in a thick layer around a spherical surface. The method includes depth equalization on the sphere, flattening, interpolation to a regular grid, pattern extraction by Shaped 3D singular spectrum analysis (SSA), and interpolation back to original nuclear positions. The approach is demonstrated on several examples of gene expression in the zebrafish egg (a model system in vertebrate development). The method is tested on several different data geometries (e.g., nuclear positions) and different forms of gene expression patterns. Fully 3D datasets for developmental gene expression are becoming increasingly available; we discuss the prospects of applying 3D-SSA to data processing and analysis in this growing field.

Analysis of a spatial gene expression database for sea anemone Nematostella vectensis during early development

Background

The spatial distribution of many genes has been visualized during the embryonic development in the starlet sea anemone Nematostella vectensis in the last decade. In situ hybridization images are available in the Kahi Kai gene expression database, and a method has been developed to quantify spatial gene expression patterns of N. vectensis. In this paper, gene expression quantification is performed on a wide range of gene expression patterns from this database and descriptions of observed expression domains are stored in a separate database for further analysis.

Methods

Spatial gene expression from suitable in situ hybridization images has been quantified with the GenExp program. A correlation analysis has been performed on the resulting numerical gene expression profiles for each stage. Based on the correlated clusters of spatial gene expression and detailed descriptions of gene expression domains, various mechanisms for developmental gene expression are proposed.

Results

In the blastula and gastrula stages of development in N. vectensis, its continuous sheet of cells is partitioned into correlating gene expression domains. During progressing development, these regions likely correspond to different fates. A statistical analysis shows that genes generally remain expressed during the planula stages in those major regions that they occupy at the end of gastrulation.

Discussion

Observed shifts in gene expression domain boundaries suggest that elongation in the planula stage mainly occurs in the vegetal ring under the influence of the gene Rx. The secondary body axis in N. vectensis is proposed to be determined at the mid blastula transition.

Conclusions

Early gene expression domains in N. vectensis appear to maintain a positional order along the primary body axis. Early determination in N. vectensis occurs in two stages: expression in broad circles and rings in the blastula is consolidated during gastrulation, and more complex expression patterns appear in the planula within these broad regions. Quantification and comparison of gene expression patterns across a database can generate hypotheses about collective cell movements before these movements are measured directly.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0209-4) contains supplementary material, which is available to authorized users.

Dose-dependent nuclear β-catenin response segregates endomesoderm along the sea star primary axis

In many invertebrates, the nuclearization of β-catenin at one pole of the embryo initiates endomesoderm specification. An intriguing possibility is that a gradient of nuclear β-catenin (nβ-catenin), similar to that operating in vertebrate neural tube patterning, functions to distinguish cell fates in invertebrates. To test this hypothesis, we determined the function of nβ-catenin during the early development of the sea star, which undergoes a basal deuterostomal mode of embryogenesis. We show that low levels of nβ-catenin activity initiate bra, which is expressed in the future posterior endoderm-fated territory; intermediate levels are required for expression of foxa and gata4/5/6, which are later restricted to the endoderm; and activation of ets1 and erg in the mesoderm-fated territory requires the highest nβ-catenin activity. Transcription factors acting downstream of high nβ-catenin segregate the endoderm/mesoderm boundary, which is further reinforced by Delta/Notch signaling. Significantly, therefore, in sea stars, endomesoderm segregation arises through transcriptional responses to levels of nβ-catenin activity. Here, we describe the first empirical evidence of a dose-dependent response to a dynamic spatiotemporal nβ-catenin activity that patterns cell fates along the primary axis in an invertebrate.

Current directions and future perspectives from the third Nematostella research conference

The third Nematostella vectensis Research Conference took place in December 2013 in Eilat, Israel, as a satellite to the 8th International Conference on Coelenterate Biology. The starlet sea anemone, N. vectensis, has emerged as a powerful cnidarian model, in large part due to the extensive genomic and transcriptomic resources and molecular approaches that are becoming available for Nematostella, which were the focus of several presentations. In addition, research was presented highlighting the broader utility of this species for studies of development, circadian rhythms, signal transduction, and gene-environment interactions.

A digital framework to build, visualize and analyze a gene expression atlas with cellular resolution in zebrafish early embryogenesis

A gene expression atlas is an essential resource to quantify and understand the multiscale processes of embryogenesis in time and space. The automated reconstruction of a prototypic 4D atlas for vertebrate early embryos, using multicolor fluorescence in situ hybridization with nuclear counterstain, requires dedicated computational strategies. To this goal, we designed an original methodological framework implemented in a software tool called Match-IT. With only minimal human supervision, our system is able to gather gene expression patterns observed in different analyzed embryos with phenotypic variability and map them onto a series of common 3D templates over time, creating a 4D atlas. This framework was used to construct an atlas composed of 6 gene expression templates from a cohort of zebrafish early embryos spanning 6 developmental stages from 4 to 6.3 hpf (hours post fertilization). They included 53 specimens, 181,415 detected cell nuclei and the segmentation of 98 gene expression patterns observed in 3D for 9 different genes. In addition, an interactive visualization software, Atlas-IT, was developed to inspect, supervise and analyze the atlas. Match-IT and Atlas-IT, including user manuals, representative datasets and video tutorials, are publicly and freely available online. We also propose computational methods and tools for the quantitative assessment of the gene expression templates at the cellular scale, with the identification, visualization and analysis of coexpression patterns, synexpression groups and their dynamics through developmental stages.

We propose a workflow to map the expression domains of multiple genes onto a series of 3D templates, or “atlas”, during early embryogenesis. It was applied to the zebrafish at different stages between 4 and 6.3 hpf, generating 6 templates. Our system overcomes the lack of significant morphological landmarks in early development by relying on the expression of a reference gene (goosecoid, gsc) and nuclear staining to guide the registration of the analyzed genes. The proposed method also successfully maps gene domains from partially imaged embryos, thus allowing greater microscope magnification and cellular resolution. By using the workflow to construct a spatiotemporal database of zebrafish, we opened the way to a systematic analysis of vertebrate embryogenesis. The atlas database, together with the mapping software (Match-IT), a custom-made visualization platform (Atlas-IT), and step-by-step user guides are available from the Supplementary Material. We expect that this will encourage other laboratories to generate, map, visualize and analyze new gene expression datasets.

Developmental gene regulatory network evolution: insights from comparative studies in echinoderms

One of the central concerns of Evolutionary Developmental biology is to understand how the specification of cell types can change during evolution. In the last decade, developmental biology has progressed toward a systems level understanding of cell specification processes. In particular, the focus has been on determining the regulatory interactions of the repertoire of genes that make up gene regulatory networks (GRNs). Echinoderms provide an extraordinary model system for determining how GRNs evolve. This review highlights the comparative GRN analyses arising from the echinoderm system. This work shows that certain types of GRN subcircuits or motifs, i.e., those involving positive feedback, tend to be conserved and may provide a constraint on development. This conservation may be due to a required arrangement of transcription factor binding sites in cis regulatory modules. The review will also discuss ways in which novelty may arise, in particular through the co-option of regulatory genes and subcircuits. The development of the sea urchin larval skeleton, a novel feature that arose in echinoderms, has provided a model for study of co-option mechanisms. Finally, the types of GRNs that can permit the great diversity in the patterns of ciliary bands and their associated neurons found among these taxa are discussed. The availability of genomic resources is rapidly expanding for echinoderms, including genome sequences not only for multiple species of sea urchins but also a species of sea star, sea cucumber, and brittle star. This will enable echinoderms to become a particularly powerful system for understanding how developmental GRNs evolve.

Targeting the Wnt pathway in zebrafish as a screening method to identify novel therapeutic compounds

Activating mutations in the Wnt signaling pathway account for the initiation of greater than 90% of all colorectal cancers and this pathway has been implicated in numerous other diseases. Therefore, identifying small molecule inhibitors of this pathway is of critical importance towards identifying clinically relevant drugs. Numerous screens have been employed to identify therapeutic reagents, but none have made it to advanced clinical trials, suggesting that traditional screening methods are ineffective at identifying clinically relevant targets. Here, we describe a novel in vivo screen to identify small molecule inhibitors of the Wnt pathway. Specifically, treatment of zebrafish embryos with LiCl inhibits GSK3 kinase function, resulting in hyperactivation of the signaling pathway and an eyeless phenotype at 1 day post fertilization. Using the small molecule XAV939, a known inhibitor of Wnt signaling, we rescued the LiCl induced eyeless phenotype, confirming efficacy of the screen. We next tested our assay with 400 known small molecule kinase inhibitors, none of which should inhibit Wnt signaling below the level of GSK3 based on their known targets. Accordingly, none of these small molecules rescued the eyeless phenotype, which demonstrates the stringency of the assay. However, several of these small molecule kinase inhibitors did generate a non-Wnt phenotype in accordance with the kinase they targeted. Therefore, combining the efficacy, sensitivity, and stringency of this preliminary screen, this model will provide an alternative to the traditional in vitro screen, generating potentially clinical relevant drugs in a rapid and cost-effective way.

Transcriptome dynamics and diversity in the early zebrafish embryo

Recent years advances in high-throughput sequencing have improved our understanding of how transcripts regulate early vertebrate development. Here, we review the transcriptome dynamics and diversity during early stages of zebrafish embryogenesis. Transcriptome dynamics is characterized by different patterns of mRNA degradation, activation of dormant transcripts and onset of transcription. Several studies have shown a striking diversity of both coding and non-coding transcripts. However, in the aftermath of this immense increase in data, functional studies of both protein-coding and non-coding transcripts are lagging behind. We anticipate that the forthcoming years will see studies relying on different high-throughput sequencing technologies and genomic tools developed for zebrafish embryos to further pin down yet un-annotated transcript-function relationships.

Genetically encoded light-activated transcription for spatiotemporal control of gene expression and gene silencing in mammalian cells

Photocaging provides a method to spatially and temporally control biological function and gene expression with high resolution. Proteins can be photochemically controlled through the site-specific installation of caging groups on amino acid side chains that are essential for protein function. The photocaging of a synthetic gene network using unnatural amino acid mutagenesis in mammalian cells was demonstrated with an engineered bacteriophage RNA polymerase. A caged T7 RNA polymerase was expressed in cells with an expanded genetic code and used in the photochemical activation of genes under control of an orthogonal T7 promoter, demonstrating tight spatial and temporal control. The synthetic gene expression system was validated with two reporter genes (luciferase and EGFP) and applied to the light-triggered transcription of short hairpin RNA constructs for the induction of RNA interference.

Full transcriptome analysis of early dorsoventral patterning in zebrafish

Understanding the molecular interactions that lead to the establishment of the major body axes during embryogenesis is one of the main goals of developmental biology. Although the past two decades have revolutionized our knowledge about the genetic basis of these patterning processes, the list of genes involved in axis formation is unlikely to be complete. In order to identify new genes involved in the establishment of the dorsoventral (DV) axis during early stages of zebrafish embryonic development, we employed next generation sequencing for full transcriptome analysis of normal embryos and embryos lacking overt DV pattern. A combination of different statistical approaches yielded 41 differentially expressed candidate genes and we confirmed by in situ hybridization the early dorsal expression of 32 genes that are transcribed shortly after the onset of zygotic transcription. Although promoter analysis of the validated genes suggests no general enrichment for the binding sites of early acting transcription factors, most of these genes carry “bivalent” epigenetic histone modifications at the time when zygotic transcription is initiated, suggesting a “poised” transcriptional status. Our results reveal some new candidates of the dorsal gene regulatory network and suggest that a plurality of the earliest upregulated genes on the dorsal side have a role in the modulation of the canonical Wnt pathway.

Differential transcript isoform usage pre- and post-zygotic genome activation in zebrafish

Background

Zebrafish embryos are transcriptionally silent until activation of the zygotic genome during the 10^th cell cycle. Onset of transcription is followed by cellular and morphological changes involving cell speciation and gastrulation. Previous genome-wide surveys of transcriptional changes only assessed gene expression levels; however, recent studies have shown the necessity to map isoform-specific transcriptional changes. Here, we perform isoform discovery and quantification on transcriptome sequences from before and after zebrafish zygotic genome activation (ZGA).

Results

We identify novel isoforms and isoform switches during ZGA for genes related to cell adhesion, pluripotency and DNA methylation. Isoform switching events include alternative splicing and changes in transcriptional start sites and in 3’ untranslated regions. New isoforms are identified even for well-characterized genes such as pou5f1, sall4 and dnmt1. Genes involved in cell-cell interactions such as f11r and magi1 display isoform switches with alterations of coding sequences. We also detect over 1000 transcripts that acquire a longer 3’ terminal exon when transcribed by the zygote compared to their maternal transcript counterparts. ChIP-sequencing data mapped onto skipped exon events reveal a correlation between histone H3K36 trimethylation peaks and skipped exons, suggesting epigenetic marks being part of alternative splicing regulation.

Conclusions

The novel isoforms and isoform switches reported here include regulators of transcriptional, cellular and morphological changes taking place around ZGA. Our data display an array of isoform-related functional changes and represent a valuable resource complementary to existing early embryo transcriptomes.

Spatial gene expression quantification: a tool for analysis of in situ hybridizations in sea anemone Nematostella vectensis

BACKGROUND

Spatial gene expression quantification is required for modeling gene regulation in developing organisms. The fruit fly Drosophila melanogaster is the model system most widely applied for spatial gene expression analysis due to its unique embryonic properties: the shape does not change significantly during its early cleavage cycles and most genes are differentially expressed along a straight axis. This system of development is quite exceptional in the animal kingdom.In the sea anemone Nematostella vectensis the embryo changes its shape during early development; there are cell divisions and cell movement, like in most other metazoans. Nematostella is an attractive case study for spatial gene expression since its transparent body wall makes it accessible to various imaging techniques.

FINDINGS

Our new quantification method produces standardized gene expression profiles from raw or annotated Nematostella in situ hybridizations by measuring the expression intensity along its cell layer. The procedure is based on digital morphologies derived from high-resolution fluorescence pictures. Additionally, complete descriptions of nonsymmetric expression patterns have been constructed by transforming the gene expression images into a three-dimensional representation.

CONCLUSIONS

We created a standard format for gene expression data, which enables quantitative analysis of in situ hybridizations from embryos with various shapes in different developmental stages. The obtained expression profiles are suitable as input for optimization of gene regulatory network models, and for correlation analysis of genes from dissimilar Nematostella morphologies. This approach is potentially applicable to many other metazoan model organisms and may also be suitable for processing data from three-dimensional imaging techniques.

BioTapestry: a tool to visualize the dynamic properties of gene regulatory networks

BioTapestry is an open source, freely available software tool that has been developed to handle the -challenges of modeling genetic regulatory networks (GRNs). Using BioTapestry, a researcher can -construct a network model and use it to visualize and understand the dynamic behavior of a complex, spatially and temporally distributed GRN. Here we provide a step-by-step example of a way to use BioTapestry to build a GRN model and discuss some common issues that can arise during this process.

An evolutionarily conserved kernel of gata5, gata6, otx2 and prdm1a operates in the formation of endoderm in zebrafish

An evolutionarily conserved subcircuit (kernel) dedicated to a specific developmental function is found at the top of the gene regulatory networks (GRNs) hierarchy. Here we comprehensively demonstrate that a pan-deuterostome endoderm specification kernel exists in zebrafish. We analyzed interactions among gata5, gata6, otx2 and prdm1a using specific morpholino knockdowns and measured the gene expression profiles by quantitative real-time RT-PCR and in situ hybridization. The mRNA rescue experiment validated the specificity of the morpholino knockdown. We found that the interactions among gata5, gata6, otx2 and prdm1a determine the initial specification of the zebrafish endoderm. Although otx2 can activate both gata5 and gata6, and the prdm1a/krox homologue also activates some endoderm transcription factors, a feedback loop from Gata to otx2 and prdm1a is missing. Furthermore, we found the positive regulation between gata5 and gata6 to further lock-on the mesendoderm specification by the Gata family. Chromatin immunoprecipitation was used to further validate the recruitment of Otx2 to the gata5 and gata6 loci. Functional assays revealed that module B of gata6 and the basal promoter of gata5 drive the gene at the mesendoderm, and mutational analysis demonstrated that Otx2 and Gata5/6 contribute to reporter gene activation. This is the first direct evidence for evolutionarily conserved endoderm specification across echinoderms and vertebrates.

Emerging properties of animal gene regulatory networks

Gene regulatory networks (GRNs) provide system level explanations of developmental and physiological functions in the terms of the genomic regulatory code. Depending on their developmental functions, GRNs differ in their degree of hierarchy, and also in the types of modular sub-circuit of which they are composed, although there is a commonly employed sub-circuit repertoire. Mathematical modelling of some types of GRN sub-circuit has deepened biological understanding of the functions they mediate. The structural organization of various kinds of GRN reflects their roles in the life process, and causally illuminates both developmental and evolutionary process.

Size matters: the contribution of cell proliferation to the progression of the specification Drosophila eye gene regulatory network

Organ development is a complex process in which the activity of scores of interacting transcription factors and signaling pathways need to be integrated with proliferative growth. Developmental gene regulatory networks (GRNs) allow capturing essential regulatory pathways at a systems-level and provide an effective way of approaching such complexity. However typical GRNs studies focus on very early embryonic stages (usually pre-gastrulation) or late stages, when there is little or no cell proliferation, and therefore do not consider how organ growth is integrated in the developmental process. This can be conveniently investigated in the Drosophila melanogaster eye primordium. Here we present a working model meant to illustrate how during a critical period, the second larval stage, changes in cells' proliferative pattern are coordinated with the initiation of the Retinal Determination (RD) gene program. Such changes are regulated in response to two different sources of signal (Wnt1/wg and BMP2/4/dpp) produced by the anterior and posterior ends of the primordium, respectively. The dpp signaling is necessary to trigger the RD program. However in order for it to be effective, cells receiving Dpp have to be out of the wg signaling range. This is obtained thanks to the proliferative growth that precedes the onset of RD expression. With this network model many of the gene regulatory steps previously known to participate in growth and patterning are linked. Analysis of the model highlights a few essential regulatory principles, as well as poses new questions. In addition, these principles might operate during the growth and patterning of other organs.

Complexity of cis-regulatory organization of six3a during forebrain and eye development in zebrafish

BACKGROUND

Six3a belongs to the SIX family of homeodomain proteins and is expressed in the most anterior neural plate at the beginning of neurogenesis in various species. Though the function of Six3a as a crucial regulator of eye and forebrain development has been thoroughly investigated, the transcriptional regulation of six3a is not well understood.

RESULTS

To elucidate the transcriptional regulation of six3a, we performed an in vivo reporter assay. Alignment of the 21-kb region surrounding the zebrafish six3a gene with the analogous region from different species identified several conserved non-coding modules. Transgenesis in zebrafish identified two enhancer elements and one suppressor. The D module drives the GFP reporter in the forebrain and eyes at an early stage, while the A module is responsible for the later expression. The A module also works as a repressor suppressing ectopic expression from the D module. Mutational analysis further minimized the A module to four highly conserved elements and the D module to three elements. Using electrophoresis mobility shift assays, we also provided evidence for the presence of DNA-binding proteins in embryonic nuclear extracts. The transcription factors that may occupy those highly conserved elements were also predicted.

CONCLUSION

This study provides a comprehensive view of six3a transcription regulation during brain and eye development and offers an opportunity to establish the gene regulatory networks underlying neurogenesis in zebrafish.

Transcriptional rewiring of the sex determining dmrt1 gene duplicate by transposable elements

Control and coordination of eukaryotic gene expression rely on transcriptional and posttranscriptional regulatory networks. Evolutionary innovations and adaptations often require rapid changes of such networks. It has long been hypothesized that transposable elements (TE) might contribute to the rewiring of regulatory interactions. More recently it emerged that TEs might bring in ready-to-use transcription factor binding sites to create alterations to the promoters by which they were captured. A process where the gene regulatory architecture is of remarkable plasticity is sex determination. While the more downstream components of the sex determination cascades are evolutionary conserved, the master regulators can switch between groups of organisms even on the interspecies level or between populations. In the medaka fish (Oryzias latipes) a duplicated copy of dmrt1, designated dmrt1bY or DMY, on the Y chromosome was shown to be the master regulator of male development, similar to Sry in mammals. We found that the dmrt1bY gene has acquired a new feedback downregulation of its expression. Additionally, the autosomal dmrt1a gene is also able to regulate transcription of its duplicated paralog by binding to a unique target Dmrt1 site nested within the dmrt1bY proximal promoter region. We could trace back this novel regulatory element to a highly conserved sequence within a new type of TE that inserted into the upstream region of dmrt1bY shortly after the duplication event. Our data provide functional evidence for a role of TEs in transcriptional network rewiring for sub- and/or neo-functionalization of duplicated genes. In the particular case of dmrt1bY, this contributed to create new hierarchies of sex-determining genes.

Evolutionary innovations and adaptations often require rapid changes in gene regulation. Transposable elements constitute the most dynamic part of eukaryotic genomes. Insertions of transposable elements can influence the expression of surrounding genes by donating new regulatory elements. A longstanding hypothesis postulates that the dispersal of transposable elements may rewire regulatory links between genes, thereby changing regulatory networks and shuffling regulatory cascades. A regulatory hierarchy of remarkable plasticity is the sex determination cascade. In the course of animal evolution, new master regulators frequently replace the sex determination gene on top of the hierarchy. In the medaka fish, a duplicate of the dmrt1 transcription factor gene, dmrt1bY, has become the sex master regulator. Its ancestor dmrt1a, in contrast, has a downstream position in the sex determination cascade. We show that after the duplication of the dmrt1 gene, the new hierarchy has been established by the insertion of a transposable element into the regulatory region of the dmrt1bY gene on the sex chromosome. This transposable element, harboring a Dmrt1 binding site, enables the self- and cross-regulation of dmrt1bY expression by Dmrt1 proteins. Our study therefore provides strong evidence for the important role of transposable elements in the rewiring of gene regulatory networks.

Vertebrate endoderm development and organ formation

The endoderm germ layer contributes to the respiratory and gastrointestinal tracts and to all of their associated organs. Over the past decade, studies in vertebrate model organisms, including frog, fish, chick, and mouse, have greatly enhanced our understanding of the molecular basis of endoderm organ development. We review this progress with a focus on early stages of endoderm organogenesis including endoderm formation, gut tube morphogenesis and patterning, and organ specification. Lastly, we discuss how developmental mechanisms that regulate endoderm organogenesis are used to direct differentiation of embryonic stem cells into specific adult cell types, which function to alleviate disease symptoms in animal models.

myGRN: a database and visualisation system for the storage and analysis of developmental genetic regulatory networks

Background

Biological processes are regulated by complex interactions between transcription factors and signalling molecules, collectively described as Genetic Regulatory Networks (GRNs). The characterisation of these networks to reveal regulatory mechanisms is a long-term goal of many laboratories. However compiling, visualising and interacting with such networks is non-trivial. Current tools and databases typically focus on GRNs within simple, single celled organisms. However, data is available within the literature describing regulatory interactions in multi-cellular organisms, although not in any systematic form. This is particularly true within the field of developmental biology, where regulatory interactions should also be tagged with information about the time and anatomical location of development in which they occur.

Description

We have developed myGRN (), a web application for storing and interrogating interaction data, with an emphasis on developmental processes. Users can submit interaction and gene expression data, either curated from published sources or derived from their own unpublished data. All interactions associated with publications are publicly visible, and unpublished interactions can only be shared between collaborating labs prior to publication. Users can group interactions into discrete networks based on specific biological processes. Various filters allow dynamic production of network diagrams based on a range of information including tissue location, developmental stage or basic topology. Individual networks can be viewed using myGRV, a tool focused on displaying developmental networks, or exported in a range of formats compatible with third party tools. Networks can also be analysed for the presence of common network motifs. We demonstrate the capabilities of myGRN using a network of zebrafish interactions integrated with expression data from the zebrafish database, ZFIN.

Conclusion

Here we are launching myGRN as a community-based repository for interaction networks, with a specific focus on developmental networks. We plan to extend its functionality, as well as use it to study networks involved in embryonic development in the future.

Visualization, documentation, analysis, and communication of large-scale gene regulatory networks

Genetic regulatory networks (GRNs) are complex, large-scale, and spatially and temporally distributed. These characteristics impose challenging demands on software tools for building GRN models, and so there is a need for custom tools. In this paper, we report on our ongoing development of BioTapestry, an open source, freely available computational tool designed specifically for building GRN models. We also outline our future development plans, and give some examples of current applications of BioTapestry.