Developmental Gene Regulatory Networks Research Articles

Network architecture and regulatory logic in neural crest development.

The neural crest is an ectodermal cell population that gives rise to over 30 cell types during vertebrate embryogenesis. These stem cells are formed at the border of the developing central nervous system and undergo extensive migration before differentiating into components of multiple tissues and organs. Neural crest formation and differentiation is a multistep process, as these cells transition through sequential regulatory states before adopting their adult phenotype. Such changes are governed by a complex gene regulatory network (GRN) that integrates environmental and cell-intrinsic inputs to regulate cell identity. Studies of neural crest cells in a variety of vertebrate models have elucidated the function and regulation of dozens of the molecular players that are part of this network. The neural crest GRN has served as a platform to explore the molecular control of multipotency, cell differentiation, and the evolution of vertebrates. In this review, we employ this genetic program as a stepping-stone to explore the architecture and the regulatory principles of developmental GRNs. We also discuss how modern genomic approaches can further expand our understanding of genetic networks in this system and others. This article is categorized under: Physiology > Mammalian Physiology in Health and Disease Biological Mechanisms > Cell Fates Developmental Biology > Lineages Models of Systems Properties and Processes > Cellular Models.

Immune activity at the gut epithelium in the larval sea urchin.

The embryo of the purple sea urchin has been a fruitful model for the study of developmental gene regulatory networks. For similar reasons, the feeding sea urchin larva provides a gene regulatory model to investigate immune interactions at the gut epithelium. Here we describe what is known of the gut structure and immune cells of the sea urchin larva, and the cellular and gene expression response of the larva to gut-associated immune challenge. As a focused example of how the sea urchin larva can be compared with vertebrate systems, we discuss the expression and function of the IL-17 signalling system in the course of thelarval immune response.

Establishment of knockout adult sea urchins by using a CRISPR-Cas9 system.

Sea urchins are used as a model organism for research on developmental biology and gene regulatory networks during early development. Gene knockdown by microinjection of morpholino antisense oligonucleotide (MASO) has been used to analyze gene function in early sea urchin embryos. However, as the effect of MASO is not long lasting, it is impossible to perturb genes expressed during late development by MASO. Recent advances in genome editing technologies have enabled gene modification in various organisms. We previously reported genome editing in the sea urchin Hemicentrotus pulcherrimus using zinc-finger nuclease (ZFN) and transcription activator-like effector nuclease (TALEN); however, the efficiencies of these technologies were not satisfactory. Here, we applied clustered regularly interspaced short palindromic repeat (CRISPR)-CRISPR-associated nuclease 9 (Cas9) technology to knock out the Pks1 gene in H.pulcherrimus. When sgRNAs targeting Pks1, which is required for the biosynthesis of larval pigment, were microinjected into fertilized eggs with SpCas9 mRNA, high-efficiency mutagenesis was achieved within 24hr post fertilization and SpCas9/sgRNA-injected pluteus larvae had an albino phenotype. One of the sgRNAs yielded 100% mutagenesis efficiency, and no off-target effect was detected. In addition, the albino phenotype was maintained in juvenile sea urchins after metamorphosis, and the knockout sea urchins survived for at least oneyear and grew to albino adult sea urchins. These findings suggest that knockout adult sea urchins were successfully established and the CRISPR-Cas9 system is a feasible method for analyzing gene functions from late developmental to adult stage.

Repressive Gene Regulation Synchronizes Development with Cellular Metabolism

Metabolic conditions affect the developmental tempo of animals. Developmental gene regulatory networks (GRNs) must therefore synchronize their dynamics with a variable timescale. We find that layered repression of genes couples GRN output with variable metabolism. When repressors of transcription or mRNA and protein stability are lost, fewer errors in Drosophila development occur when metabolism is lowered. We demonstrate the universality of this phenomenon by eliminating the entire microRNA family of repressors and find that development to maturity can be largely rescued when metabolism is reduced. Using a mathematical model that replicates GRN dynamics, we find that lowering metabolism suppresses the emergence of developmental errors by curtailing the influence of auxiliary repressors on GRN output. We experimentally show that gene expression dynamics are less affected by loss of repressors when metabolism is reduced. Thus, layered repression provides robustness through error suppression and may provide an evolutionary route to a shorter reproductive cycle.

From Gregor Mendel to Eric Davidson: Mathematical Models and Basic Principles in Biology.

Mathematical models have been widespread in biology since its emergence as a modern experimental science in the 19th century. Focusing on models in developmental biology and heredity, this article (1) presents the properties and epistemological basis of pertinent mathematical models in biology from Mendel's model of heredity in the 19th century to Eric Davidson's model of developmental gene regulatory networks in the 21st; (2) shows that the models differ not only in their epistemologies but also in regard to explicitly or implicitly taking into account basic biological principles, in particular those of biological specificity (that became, in part, replaced by genetic information) and genetic causality. The article claims that models disregarding these principles did not impact the direction of biological research in a lasting way, although some of them, such as D'Arcy Thompson's models of biological form, were widely read and admired and others, such as Turing's models of development, stimulated research in other fields. Moreover, it suggests that successful models were not purely mathematical descriptions or simulations of biological phenomena but were based on inductive, as well as hypothetico-deductive, methodology. The recent availability of large amounts of sequencing data and new computational methodology tremendously facilitates system approaches and pattern recognition in many fields of research. Although these new technologies have given rise to claims that correlation is replacing experimentation and causal analysis, the article argues that the inductive and hypothetico-deductive experimental methodologies have remained fundamentally important as long as causal-mechanistic explanations of complex systems are pursued.

Quantitative Comparison of the Anterior-Posterior Patterning System in the Embryos of Five Drosophila Species.

Complex spatiotemporal gene expression patterns direct the development of the fertilized egg into an adult animal. Comparisons across species show that, in spite of changes in the underlying regulatory DNA sequence, developmental programs can be maintained across millions of years of evolution. Reciprocally, changes in gene expression can be used to generate morphological novelty. Distinguishing between changes in regulatory DNA that lead to changes in gene expression and those that do not is therefore a central goal of evolutionary developmental biology. Quantitative, spatially-resolved measurements of developmental gene expression patterns play a crucial role in this goal, enabling the detection of subtle phenotypic differences between species and the development of computations models that link the sequence of regulatory DNA to expression patterns. Here we report the generation of two atlases of cellular resolution gene expression measurements for the primary anterior-posterior patterning genes in Drosophila simulans and Drosophila virilis. By combining these data sets with existing atlases for three other Drosophila species, we detect subtle differences in the gene expression patterns and dynamics driving the highly conserved axis patterning system and delineate inter-species differences in the embryonic morphology. These data sets will be a resource for future modeling studies of the evolution of developmental gene regulatory networks.

Modularity, criticality, and evolvability of a developmental gene regulatory network.

The existence of discrete phenotypic traits suggests that the complex regulatory processes which produce them are functionally modular. These processes are usually represented by networks. Only modular networks can be partitioned into intelligible subcircuits able to evolve relatively independently. Traditionally, functional modularity is approximated by detection of modularity in network structure. However, the correlation between structure and function is loose. Many regulatory networks exhibit modular behaviour without structural modularity. Here we partition an experimentally tractable regulatory network-the gap gene system of dipteran insects-using an alternative approach. We show that this system, although not structurally modular, is composed of dynamical modules driving different aspects of whole-network behaviour. All these subcircuits share the same regulatory structure, but differ in components and sensitivity to regulatory interactions. Some subcircuits are in a state of criticality, while others are not, which explains the observed differential evolvability of the various expression features in the system.

Identification of evolutionarily conserved downstream core promoter elements required for the transcriptional regulation of Fushi tarazu target genes.

The regulation of transcription initiation is critical for developmental and cellular processes. RNA polymerase II (Pol II) is recruited by the basal transcription machinery to the core promoter where Pol II initiates transcription. The core promoter encompasses the region from -40 to +40 bp relative to the +1 transcription start site (TSS). Core promoters may contain one or more core promoter motifs that confer specific properties to the core promoter, such as the TATA box, initiator (Inr) and motifs that are located downstream of the TSS, namely, motif 10 element (MTE), the downstream core promoter element (DPE) and the Bridge, a bipartite core promoter element. We had previously shown that Caudal, an enhancer-binding homeodomain transcription factor and a key regulator of the Hox gene network, is a DPE-specific activator. Interestingly, pair-rule proteins have been implicated in enhancer-promoter communication at the engrailed locus. Fushi tarazu (Ftz) is an enhancer-binding homeodomain transcription factor encoded by the ftz pair-rule gene. Ftz works in concert with its co-factor, Ftz-F1, to activate transcription. Here, we examined whether Ftz and Ftz-F1 activate transcription with a preference for a specific core promoter motif. Our analysis revealed that similarly to Caudal, Ftz and Ftz-F1 activate the promoter containing a TATA box mutation to significantly higher levels than the promoter containing a DPE mutation, thus demonstrating a preference for the DPE motif. We further discovered that Ftz target genes are enriched for a combination of functional downstream core promoter elements that are conserved among Drosophila species. Thus, the unique combination (Inr, Bridge and DPE) of functional downstream core promoter elements within Ftz target genes highlights the complexity of transcriptional regulation via the core promoter in the transcription of different developmental gene regulatory networks.

Systematic analysis of the Myxococcus xanthus developmental gene regulatory network supports posttranslational regulation of FruA by C-signaling.

Upon starvation Myxococcus xanthus undergoes multicellular development. Rod-shaped cells move into mounds in which some cells differentiate into spores. Cells begin committing to sporulation at 24-30h poststarvation, but the mechanisms governing commitment are unknown. FruA and MrpC are transcription factors that are necessary for commitment. They bind cooperatively to promoter regions and activate developmental gene transcription, including that of the dev operon. Leading up to and during the commitment period, dev mRNA increased in wild type, but not in a mutant defective in C-signaling, a short-range signaling interaction between cells that is also necessary for commitment. The C-signaling mutant exhibited ~20-fold less dev mRNA than wild type at 30h poststarvation, despite a similar level of MrpC and only 2-fold less FruA. Boosting the FruA level twofold in the C-signaling mutant had little effect on the dev mRNA level, and dev mRNA was not less stable in the C-signaling mutant. Neither did high cooperativity of MrpC and FruA binding upstream of the dev promoter explain the data. Rather, our systematic experimental and computational analyses support a model in which C-signaling activates FruA at least ninefold posttranslationally in order to commit a cell to spore formation.

The evolution of gene regulatory networks controlling Arabidopsis thaliana L. trichome development

BackgroundThe variation in structure and function of gene regulatory networks (GRNs) participating in organisms development is a key for understanding species-specific evolutionary strategies. Even the tiniest modification of developmental GRN might result in a substantial change of a complex morphogenetic pattern. Great variety of trichomes and their accessibility makes them a useful model for studying the molecular processes of cell fate determination, cell cycle control and cellular morphogenesis. Nowadays, a large number of genes regulating the morphogenesis of A. thaliana trichomes are described. Here we aimed at a study the evolution of the GRN defining the trichome formation, and evaluation its importance in other developmental processes.ResultsIn study of the evolution of trichomes formation GRN we combined classical phylogenetic analysis with information on the GRN topology and composition in major plants taxa. This approach allowed us to estimate both times of evolutionary emergence of the GRN components which are mainly proteins, and the relative rate of their molecular evolution. Various simplifications of protein structure (based on the position of amino acid residues in protein globula, secondary structure type, and structural disorder) allowed us to demonstrate the evolutionary associations between changes in protein globules and speciations/duplications events. We discussed their potential involvement in protein-protein interactions and GRN function.ConclusionsWe hypothesize that the divergence and/or the specialization of the trichome-forming GRN is linked to the emergence of plant taxa. Information about the structural targets of the protein evolution in the GRN may predict switching points in gene networks functioning in course of evolution. We also propose a list of candidate genes responsible for the development of trichomes in a wide range of plant species.

Using ATAC-seq and RNA-seq to increase resolution in GRN connectivity.

Echinoderms have some of the most complete reconstructed developmental gene regulatory networks (GRN) of any embryo, accounting for the formation of most embryo tissues and organs. Yet, many nodes (genes and regulators) and their regulatory interactions are still to be uncovered. Traditionally, knockdown/knockout experiments are performed to determine regulator-gene interactions, which are individually validated by cis-regulatory analysis. Differential RNA-seq, combined with perturbation analysis, allows for genome-wide reconstruction of a GRN around given regulators; however, this level of resolution cannot determine direct interactions. ChiP-chip or ChIP-seq is better equipped for determining, genome-wide, whether binding of a given transcription factor (TF) to cis-regulatory elements occurs. Antibodies for the TFs of interest must be available, and if not, this presents a limiting factor. ATAC-seq identifies regions of open chromatin, that are typically trimethylated at H3K4, H3K36 and H3K79 (Kouzarides, 2007), for a given time point, condition, or tissue. This technology combined with RNA-seq and perturbation analysis provides high resolution of the possible functional interactions occurring during development. Additionally, ATAC-seq is less expensive than ChIP-seq, requires less starting material, and provides a global view of regulatory regions. This chapter provides detailed steps to identify potential regulatory relationships between the nodes of a GRN, given a well assembled genome, annotated with gene models, and ATAC-seq data combined with RNA-seq and knockdown experiments.

Conserved regulatory state expression controlled by divergent developmental gene regulatory networks in echinoids.

Evolution of the animal body plan is driven by changes in developmental gene regulatory networks (GRNs), but how networks change to control novel developmental phenotypes remains, in most cases, unresolved. Here, we address GRN evolution by comparing the endomesoderm GRN in two echinoid sea urchins, Strongylocentrotus purpuratus and Eucidaris tribuloides, with at least 268 million years of independent evolution. We first analyzed the expression of twelve transcription factors and signaling molecules of the S. purpuratus GRN in E. tribuloides embryos, showing that orthologous regulatory genes are expressed in corresponding endomesodermal cell fates in the two species. However, perturbation of regulatory genes revealed that important regulatory circuits of the S. purpuratus GRN are significantly different in E. tribuloides For example, mesodermal Delta/Notch signaling controls exclusion of alternative cell fates in E. tribuloides but controls mesoderm induction and activation of a positive feedback circuit in S. purpuratus These results indicate that the architecture of the sea urchin endomesoderm GRN evolved by extensive gain and loss of regulatory interactions between a conserved set of regulatory factors that control endomesodermal cell fate specification.

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.

Ancient mechanisms for the evolution of the bicoid homeodomain's function in fly development.

The ancient mechanisms that caused developmental gene regulatory networks to diversify among distantly related taxa are not well understood. Here we use ancestral protein reconstruction, biochemical experiments, and developmental assays of transgenic animals carrying reconstructed ancestral genes to investigate how the transcription factor Bicoid (Bcd) evolved its central role in anterior-posterior patterning in flies. We show that most of Bcd's derived functions are attributable to evolutionary changes within its homeodomain (HD) during a phylogenetic interval >140 million years ago. A single substitution from this period (Q50K) accounts almost entirely for the evolution of Bcd's derived DNA specificity in vitro. In transgenic embryos expressing the reconstructed ancestral HD, however, Q50K confers activation of only a few of Bcd's transcriptional targets and yields a very partial rescue of anterior development. Adding a second historical substitution (M54R) confers regulation of additional Bcd targets and further rescues anterior development. These results indicate that two epistatically interacting mutations played a major role in the evolution of Bcd's controlling regulatory role in early development. They also show how ancestral sequence reconstruction can be combined with in vivo characterization of transgenic animals to illuminate the historical mechanisms of developmental evolution.

Ankyrin domain encoding genes from an ancient horizontal transfer are functionally integrated into Nasonia developmental gene regulatory networks

BackgroundHow regulatory networks incorporate additional components and how novel genes are functionally integrated into well-established developmental processes are two important and intertwined questions whose answers have major implications for understanding the evolution of development. We recently discovered a set of lineage-restricted genes with strong and specific expression patterns along the dorsal-ventral (DV) axis of the embryo of the wasp Nasonia that may serve as a powerful system for addressing these questions. We sought to both understand the evolutionary history of these genes and to determine their functions in the Nasonia DV patterning system.ResultsWe have found that the novel DV genes are part of a large family of rapidly duplicating and diverging ankyrin domain-encoding genes that originated most likely by horizontal transfer from a prokaryote in a common ancestor of the wasp superfamily Chalcidoidea. We tested the function of those ankyrin-encoding genes expressed along the DV axis and found that they participate in early embryonic DV patterning. We also developed a new wasp model system (Melittobia) and found that some functional integration of ankyrin genes have been preserved for over 90 million years.ConclusionsOur results indicate that regulatory networks can incorporate novel genes that then become necessary for stable and repeatable outputs. Even a modest role in developmental networks may be enough to allow novel or duplicate genes to be maintained in the genome and become fully integrated network components.

Comparative Studies of Gene Expression Kinetics: Methodologies and Insights on Development and Evolution.

Across the animal kingdom, embryos of closely related species show high morphological similarity despite genetic and environmental distances. Deciphering the molecular mechanisms that underlie morphological conservation and those that support embryonic adaptation are keys to understand developmental robustness and evolution. Comparative studies of developmental gene regulatory networks can track the genetic changes that lead to evolutionary novelties. However, these studies are limited to a relatively small set of genes and demand extensive experimental efforts. An alternative approach enabled by next-generation sequencing, is to compare the expression kinetic of large sets of genes between different species. The advantages of these comparisons are that they can be done relatively easily, for any species and they provide information of all expressed genes. The challenge in these experiments is to compare the kinetic profiles of thousands of genes between species that develop in different rates. Here we review recent comparative studies that tackled the challenges of accurate staging and large-scale analyses using different computational approaches. These studies reveal how correct temporal scaling exposes the striking conservation of developmental gene expression between morphologically similar species. Different clustering approaches are used to address various comparative questions and identify the conservation and divergence of large gene sets. We discuss the unexpected contribution of housekeeping genes to the interspecies correlations and how this contribution distorts the hourglass pattern generated by developmental genes. Overall, we demonstrate how comparative studies of gene expression kinetics can provide novel insights into the developmental constraints and plasticity that shape animal body plans.

Inference of Developmental Gene Regulatory Networks Beyond Classical Model Systems: New Approaches in the Post-genomic Era

The advent of high-throughput sequencing (HTS) technologies has revolutionized the way we understand the transformation of genetic information into morphological traits. Elucidating the network of interactions between genes that govern cell differentiation through development is one of the core challenges in genome research. These networks are known as developmental gene regulatory networks (dGRNs) and consist largely of the functional linkage between developmental control genes, cis-regulatory modules, and differentiation genes, which generate spatially and temporally refined patterns of gene expression. Over the last 20 years, great advances have been made in determining these gene interactions mainly in classical model systems, including human, mouse, sea urchin, fruit fly, and worm. This has brought about a radical transformation in the fields of developmental biology and evolutionary biology, allowing the generation of high-resolution gene regulatory maps to analyze cell differentiation during animal development. Such maps have enabled the identification of gene regulatory circuits and have led to the development of network inference methods that can recapitulate the differentiation of specific cell-types or developmental stages. In contrast, dGRN research in non-classical model systems has been limited to the identification of developmental control genes via the candidate gene approach and the characterization of their spatiotemporal expression patterns, as well as to the discovery of cis-regulatory modules via patterns of sequence conservation and/or predicted transcription-factor binding sites. However, thanks to the continuous advances in HTS technologies, this scenario is rapidly changing. Here, we give a historical overview on the architecture and elucidation of the dGRNs. Subsequently, we summarize the approaches available to unravel these regulatory networks, highlighting the vast range of possibilities of integrating multiple technical advances and theoretical approaches to expand our understanding on the global gene regulation during animal development in non-classical model systems. Such new knowledge will not only lead to greater insights into the evolution of molecular mechanisms underlying cell identity and animal body plans, but also into the evolution of morphological key innovations in animals.

Convergent evolution of gene expression in two high-toothed stickleback populations.

Changes in developmental gene regulatory networks enable evolved changes in morphology. These changes can be in cis regulatory elements that act in an allele-specific manner, or changes to the overall trans regulatory environment that interacts with cis regulatory sequences. Here we address several questions about the evolution of gene expression accompanying a convergently evolved constructive morphological trait, increases in tooth number in two independently derived freshwater populations of threespine stickleback fish (Gasterosteus aculeatus). Are convergently evolved cis and/or trans changes in gene expression associated with convergently evolved morphological evolution? Do cis or trans regulatory changes contribute more to gene expression changes accompanying an evolved morphological gain trait? Transcriptome data from dental tissue of ancestral low-toothed and two independently derived high-toothed stickleback populations revealed significantly shared gene expression changes that have convergently evolved in the two high-toothed populations. Comparing cis and trans regulatory changes using phased gene expression data from F1 hybrids, we found that trans regulatory changes were predominant and more likely to be shared among both high-toothed populations. In contrast, while cis regulatory changes have evolved in both high-toothed populations, overall these changes were distinct and not shared among high-toothed populations. Together these data suggest that a convergently evolved trait can occur through genetically distinct regulatory changes that converge on similar trans regulatory environments.

Global analysis of primary mesenchyme cell cis-regulatory modules by chromatin accessibility profiling

BackgroundThe developmental gene regulatory network (GRN) that underlies skeletogenesis in sea urchins and other echinoderms is a paradigm of GRN structure, function, and evolution. This transcriptional network is deployed selectively in skeleton-forming primary mesenchyme cells (PMCs) of the early embryo. To advance our understanding of this model developmental GRN, we used genome-wide chromatin accessibility profiling to identify and characterize PMC cis-regulatory modules (CRMs).ResultsATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) analysis of purified PMCs provided a global picture of chromatin accessibility in these cells. We used both ATAC-seq and DNase-seq (DNase I hypersensitive site sequencing) to identify > 3000 sites that exhibited increased accessibility in PMCs relative to other embryonic cell lineages, and provide both computational and experimental evidence that a large fraction of these sites represent bona fide skeletogenic CRMs. Putative PMC CRMs were preferentially located near genes differentially expressed by PMCs and consensus binding sites for two key transcription factors in the PMC GRN, Alx1 and Ets1, were enriched in these CRMs. Moreover, a high proportion of candidate CRMs drove reporter gene expression specifically in PMCs in transgenic embryos. Surprisingly, we found that PMC CRMs were partially open in other embryonic lineages and exhibited hyperaccessibility as early as the 128-cell stage.ConclusionsOur work provides a comprehensive picture of chromatin accessibility in an early embryonic cell lineage. By identifying thousands of candidate PMC CRMs, we significantly enhance the utility of the sea urchin skeletogenic network as a general model of GRN architecture and evolution. Our work also shows that differential chromatin accessibility, which has been used for the high-throughput identification of enhancers in differentiated cell types, is a powerful approach for the identification of CRMs in early embryonic cells. Lastly, we conclude that in the sea urchin embryo, CRMs that control the cell type-specific expression of effector genes are hyperaccessible several hours in advance of gene activation.

The developmental transcriptome atlas of the spoon worm Urechis unicinctus (Echiurida: Annelida).

BackgroundEchiurida is one of the most intriguing major subgroups of annelida because, unlike most other annelids, echiurids lack metameric body segmentation as adults. For this reason, transcriptome analyses from various developmental stages of echiurid species can be of substantial value for understanding precise expression levels and the complex regulatory networks during early and larval development.ResultsA total of 914 million raw RNA-Seq reads were produced from 14 developmental stages of Urechis unicinctus and were de novo assembled into contigs spanning 63,928,225 bp with an N50 length of 2700 bp. The resulting comprehensive transcriptome database of the early developmental stages of U. unicinctus consists of 20,305 representative functional protein-coding transcripts. Approximately 66% of unigenes were assigned to superphylum-level taxa, including Lophotrochozoa (40%). The completeness of the transcriptome assembly was assessed using benchmarking universal single-copy orthologs; 75.7% of the single-copy orthologs were presented in our transcriptome database. We observed 3 distinct patterns of global transcriptome profiles from 14 developmental stages and identified 12,705 genes that showed dynamic regulation patterns during the differentiation and maturation of U. unicinctus cells.ConclusionsWe present the first large-scale developmental transcriptome dataset of U. unicinctus and provide a general overview of the dynamics of global gene expression changes during its early developmental stages. The analysis of time-course gene expression data is a first step toward understanding the complex developmental gene regulatory networks in U. unicinctus and will furnish a valuable resource for analyzing the functions of gene repertoires in various developmental phases.