A gene regulatory network controlling the embryonic specification of endoderm

How Does the Regulatory Genome Work?

The regulatory genome controls genome activity throughout the life of an organism. This requires that complex information processing functions are encoded in, and operated by, the regulatory genome. Although much remains to be learned about how the regulatory genome works, we here discuss two cases where regulatory functions have been experimentally dissected in great detail and at the systems level, and formalized by computational logic models. Both examples derive from the sea urchin embryo, but assess two distinct organizational levels of genomic information processing. The first example shows how the regulatory system of a single gene, endo16, executes logic operations through individual transcription factor binding sites and cis-regulatory modules that control the expression of this gene. The second example shows information processing at the gene regulatory network (GRN) level. The GRN controlling development of the sea urchin endomesoderm has been experimentally explored at an almost complete level. A Boolean logic model of this GRN suggests that the modular logic functions encoded at the single-gene level show compositionality and suffice to account for integrated function at the network level. We discuss these examples both from a biological-experimental point of view and from a computer science-informational point of view, as both illuminate principles of how the regulatory genome works.

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.

Encounters across networks: Windows into principles of genomic regulation

Gene regulatory networks account for the ability of the genome to program development in complex multi-cellular organisms. Such networks are based on principles of gene regulation by combinations of transcription factors that bind to specific cis-regulatory DNA sites to activate transcription. These cis-regulatory regions mediate logic processing at each network node, enabling progressive increases in organismal complexity with development. Gene regulatory network explanations of development have been shown to account for patterning and cell type diversification in fly and sea urchin embryonic systems, where networks are characterized by fast coupling between transcriptional inputs and changes in target gene transcription rates, and crucial cis-regulatory elements are concentrated relatively close to the protein coding sequences of the target genes, thus facilitating their identification. Stem cell-based development in post-embryonic mammalian systems also depends on gene networks, but differs from the fly and sea urchin systems. First, the number of regulatory elements per gene and the distances between regulatory elements and the genes they control are considerably larger, forcing searches via genome-wide transcription factor binding surveys rather than functional assays. Second, the intrinsic timing of network state transitions can be slowed considerably by the need to undo stem-cell chromatin configurations, which presumably add stability to stem-cell states but retard responses to transcription factor changes during differentiation. The dispersed, partially redundant cis-regulatory systems controlling gene expression and the slow state transition kinetics in these systems already reveal new insights and opportunities to extend understanding of the repertoire of gene networks and regulatory system logic.

Single cell transcriptomics: moving towards multi-omics

Single-cell multi-omics analysis helps characterize multiple layers of molecular features at a single-cell scale to provide insights into cellular processes and functions.

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.

Inhibition of microRNA suppression of Dishevelled results in Wnt pathway-associated developmental defects in sea urchin

MicroRNAs (miRNAs) are highly conserved, small non-coding RNAs that regulate gene expressions by binding to the 3' untranslated region of target mRNAs thereby silencing translation. Some miRNAs are key regulators of the Wnt signaling pathways, which impact developmental processes. This study investigates miRNA regulation of different isoforms of Dishevelled (Dvl/Dsh), which encode a key component in the Wnt signaling pathway. The sea urchin Dvl mRNA isoforms have similar spatial distribution in early development, but one isoform is distinctively expressed in the larval ciliary band. We demonstrated that Dvl isoforms are directly suppressed by miRNAs. By blocking miRNA suppression of Dvl isoforms, we observed dose-dependent defects in spicule length, patterning of the primary mesenchyme cells, gut morphology, and cilia. These defects likely result from increased Dvl protein levels, leading to perturbation of Wnt-dependent signaling pathways and additional Dvl-mediated processes. We further demonstrated that overexpression of Dvl isoforms recapitulated some of the Dvl miRNATP-induced phenotypes. Overall, our results indicate that miRNA suppression of Dvl isoforms plays an important role in ensuring proper development and function of primary mesenchyme cells and cilia.

Canonical Wnt/β-catenin and Notch signaling regulate animal/vegetal axial patterning in the cephalochordate amphioxus

In bilaterians, animal/vegetal axial (A/V) patterning is a fundamental early developmental event for establishment of animal/vegetal polarity and following specification of the germ layers (ectoderm, mesoderm, endoderm), of which the evolutionary origin is enigmatic. Understanding A/V axial patterning in a basal animal from each phylum would help to reconstruct the ancestral state of germ layer specification in bilaterians and thus, the evolution of mesoderm, the third intermediate cell layer. Herein, data show that the canonical Wnt/β-catenin (cWnt) and Notch signaling pathways control mesoderm specification from the early endomesoderm in the basal chordate amphioxus. Amphioxus belongs to the deuterostome, one of the main superphyla in Bilateria. In the present study, genes (tcf, dsh, axin, gsk3β) encoding cWnt components were expressed in the endomesoderm during the gastrula stages. Excess cWnt signaling by BIO, a GSK3 inhibitor, expanded the expression domains of outer endomesodermal genes that include future mesodermal ones and suppressed inner endomesodermal and ectodermal genes. Interfering Notch signaling by DAPT, a γ-secretase inhibitor, resulted in decreased expression of ectodermal and endomesodermal markers. These results suggest that cWnt and Notch have important roles in mesoderm specification in amphioxus embryos. The evolution of the mesoderm is also discussed.

Ectomesoderm and epithelial-mesenchymal transition-related genes in spiralian development

BACKGROUND

Spiralians (e.g., annelids, molluscs, and flatworms) possess two sources of mesoderm. One is from endodermal precursors (endomesoderm), which is considered to be the ancestral source in metazoans. The second is from ectoderm (ectomesoderm) and may represent a novel cell type in the Spiralia. In the mollusc Crepidula fornicata, ectomesoderm is derived from micromere daughters within the A and B cell quadrants. Their progeny lie along the anterolateral edges of the blastopore. There they undergo epithelial-mesenchymal transition (EMT), become rounded and undergo delamination/ingression. Subsequently, they assume the mesenchymal phenotype, and migrate beneath the surface ectoderm to differentiate various cell types, including muscles and pigment cells.

RESULTS

We examined expression of several genes whose homologs are known to regulate Type 1 EMT in other metazoans. Most of these genes were expressed within spiralian ectomesoderm during EMT.

CONCLUSIONS

We propose that spiralian ectomesoderm, which exhibits analogous cellular behaviors to other populations of mesenchymal cells, may be controlled by the same genes that drive EMT in other metazoans. Perhaps these genes comprise a conserved metazoan EMT gene regulatory network (GRN). This study represents the first step in elucidating the GRN controlling the development of a novel spiralian cell type (ectomesoderm). Developmental Dynamics 247:1097-1120, 2018. © 2018 Wiley Periodicals, Inc.

Notch signaling in the division of germ layers in bilaterian embryos

Bilaterian embryos are triploblastic organisms which develop three complete germ layers (ectoderm, mesoderm, and endoderm). While the ectoderm develops mainly from the animal hemisphere, there is diversity in the location from where the endoderm and the mesoderm arise in relation to the animal-vegetal axis, ranging from endoderm being specified between the ectoderm and mesoderm in echinoderms, and the mesoderm being specified between the ectoderm and the endoderm in vertebrates. A common feature is that part of the mesoderm segregates from an ancient bipotential endomesodermal domain. The process of segregation is noisy during the initial steps but it is gradually refined. In this review, we discuss the role of the Notch pathway in the establishment and refinement of boundaries between germ layers in bilaterians, with special focus on its interaction with the Wnt/β-catenin pathway.

Regulatory states in the developmental control of gene expression

A growing body of evidence shows that gene expression in multicellular organisms is controlled by the combinatorial function of multiple transcription factors. This indicates that not the individual transcription factors or signaling molecules, but the combination of expressed regulatory molecules, the regulatory state, should be viewed as the functional unit in gene regulation. Here, I discuss the concept of the regulatory state and its proposed role in the genome-wide control of gene expression. Recent analyses of regulatory gene expression in sea urchin embryos have been instrumental for solving the genomic control of cell fate specification in this system. Some of the approaches that were used to determine the expression of regulatory states during sea urchin embryogenesis are reviewed. Significant developmental changes in regulatory state expression leading to the distinct specification of cell fates are regulated by gene regulatory network circuits. How these regulatory state transitions are encoded in the genome is illuminated using the sea urchin endoderm-mesoderms cell fate decision circuit as an example. These observations highlight the importance of considering developmental gene regulation, and the function of individual transcription factors, in the context of regulatory states.

Gene regulatory networks and cell lineages that underlie the formation of skeletal muscle

Skeletal muscle in vertebrates is formed by two major routes, as illustrated by the mouse embryo. Somites give rise to myogenic progenitors that form all of the muscles of the trunk and limbs. The behavior of these cells and their entry into the myogenic program is controlled by gene regulatory networks, where paired box gene 3 (Pax3) plays a predominant role. Head and some neck muscles do not derive from somites, but mainly form from mesoderm in the pharyngeal region. Entry into the myogenic program also depends on the myogenic determination factor (MyoD) family of genes, but Pax3 is not expressed in these myogenic progenitors, where different gene regulatory networks function, with T-box factor 1 (Tbx1) and paired-like homeodomain factor 2 (Pitx2) as key upstream genes. The regulatory genes that underlie the formation of these muscles are also important players in cardiogenesis, expressed in the second heart field, which is a major source of myocardium and of the pharyngeal arch mesoderm that gives rise to skeletal muscles. The demonstration that both types of striated muscle derive from common progenitors comes from clonal analyses that have established a lineage tree for parts of the myocardium and different head and neck muscles. Evolutionary conservation of the two routes to skeletal muscle in vertebrates extends to chordates, to trunk muscles in the cephlochordate Amphioxus and to muscles derived from cardiopharyngeal mesoderm in the urochordate Ciona, where a related gene regulatory network determines cardiac or skeletal muscle cell fates. In conclusion, Eric Davidson's visionary contribution to our understanding of gene regulatory networks and their evolution is acknowledged.

Assessing regulatory information in developmental gene regulatory networks

Gene regulatory networks (GRNs) provide a transformation function between the static genomic sequence and the primary spatial specification processes operating development. The regulatory information encompassed in developmental GRNs thus goes far beyond the control of individual genes. We here address regulatory information at different levels of network organization, from single node to subcircuit to large-scale GRNs and discuss how regulatory design features such as network architecture, hierarchical organization, and cis-regulatory logic contribute to the developmental function of network circuits. Using specific subcircuits from the sea urchin endomesoderm GRN, for which both circuit design and biological function have been described, we evaluate by Boolean modeling and in silico perturbations the import of given circuit features on developmental function. The examples include subcircuits encoding positive feedback, mutual repression, and coherent feedforward, as well as signaling interaction circuitry. Within the hierarchy of the endomesoderm GRN, these subcircuits are organized in an intertwined and overlapping manner. Thus, we begin to see how regulatory information encoded at individual nodes is integrated at all levels of network organization to control developmental process.

EchinoBase: Tools for Echinoderm Genome Analyses

The echinoderms are a phylum of invertebrate deuterostome animals that constitute important research models for a number of biological disciplines. EchinoBase ( www.echinobase.org ) is an echinoderm-specific genome database and web information system that provides a platform for the interrogation and exploration of echinoderm genomic data. This chapter outlines the datasets available on EchinoBase; from assembled genomes and genome annotations, to spatial and quantitative expression data, as well as functional genomics datasets. We also highlight the bioinformatic tools available on the website to facilitate rapid inquiries using these data (genome browsers, precompiled BLAST databases, etc.), and suggest optimized strategies for performing these inquiries. We conclude with a perspective on how one could integrate various genomic resources to predict putative noncoding regulatory regions. The available datasets and analyses they permit provide the basic components required for developing an understanding of how echinoderm genomes are regulated, especially during early development, and provides a platform for comparative genomic inquiries among species in this phylum.

Coarse-graining as a downward causation mechanism

Downward causation is the controversial idea that 'higher' levels of organization can causally influence behaviour at 'lower' levels of organization. Here I propose that we can gain traction on downward causation by being operational and examining how adaptive systems identify regularities in evolutionary or learning time and use these regularities to guide behaviour. I suggest that in many adaptive systems components collectively compute their macroscopic worlds through coarse-graining. I further suggest we move from simple feedback to downward causation when components tune behaviour in response to estimates of collectively computed macroscopic properties. I introduce a weak and strong notion of downward causation and discuss the role the strong form plays in the origins of new organizational levels. I illustrate these points with examples from the study of biological and social systems and deep neural networks.This article is part of the themed issue 'Reconceptualizing the origins of life'.

Notch-mediated lateral inhibition is an evolutionarily conserved mechanism patterning the ectoderm in echinoids

Notch signaling is a crucial cog in early development of euechinoid sea urchins, specifying both non-skeletogenic mesodermal lineages and serotonergic neurons in the apical neuroectoderm. Here, the spatial distributions and function of delta, gcm, and hesc, three genes critical to these processes in euechinoids, are examined in the distantly related cidaroid sea urchin Eucidaris tribuloides. Spatial distribution and experimental perturbation of delta and hesc suggest that the function of Notch signaling in ectodermal patterning in early development of E. tr ibuloides is consistent with canonical lateral inhibition. Delta transcripts were observed in t he archenteron, apical ectoderm, and lateral ectoderm in gastrulating e mbryos of E. tribuloides. Perturbation of Notch signaling by either delta morpholino or treatment of DAPT downregulated hesc and upregulated delta and gcm, resulting in ectopic expression of delta and gcm. Similarly, hesc perturbation mirrored the effects of delta perturbation. Interestingly, perturbation of delta or hesc resulted in more cells expressing gcm and supernumerary pigment cells, suggesting that pigment cell proliferation is regulated by Notch in E. tribuloides. These results are consistent with an evolutionary scenario whereby, in the echinoid ancestor, Notch signaling was deployed in the ectoderm to specify neurogenic progenitors and controlled pigment cell proliferation in the dorsal ectoderm.

Programming Morphogenesis through Systems and Synthetic Biology

Mammalian tissue development is an intricate, spatiotemporal process of self-organization that emerges from gene regulatory networks of differentiating stem cells. A major goal in stem cell biology is to gain a sufficient understanding of gene regulatory networks and cell-cell interactions to enable the reliable and robust engineering of morphogenesis. Here, we review advances in synthetic biology, single cell genomics, and multiscale modeling, which, when synthesized, provide a framework to achieve the ambitious goal of programming morphogenesis in complex tissues and organoids.

New insights from a high-resolution look at gastrulation in the sea urchin, Lytechinus variegatus

BACKGROUND

Gastrulation is a complex orchestration of movements by cells that are specified early in development. Until now, classical convergent extension was considered to be the main contributor to sea urchin archenteron extension, and the relative contributions of cell divisions were unknown. Active migration of cells along the axis of extension was also not considered as a major factor in invagination.

RESULTS

Cell transplantations plus live imaging were used to examine endoderm cell morphogenesis during gastrulation at high-resolution in the optically clear sea urchin embryo. The invagination sequence was imaged throughout gastrulation. One of the eight macromeres was replaced by a fluorescently labeled macromere at the 32 cell stage. At gastrulation those patches of fluorescent endoderm cell progeny initially about 4 cells wide, released a column of cells about 2 cells wide early in gastrulation and then often this column narrowed to one cell wide by the end of archenteron lengthening. The primary movement of the column of cells was in the direction of elongation of the archenteron with the narrowing (convergence) occurring as one of the two cells moved ahead of its neighbor. As the column narrowed, the labeled endoderm cells generally remained as a contiguous population of cells, rarely separated by intrusion of a lateral unlabeled cell. This longitudinal cell migration mechanism was assessed quantitatively and accounted for almost 90% of the elongation process. Much of the extension was the contribution of Veg2 endoderm with a minor contribution late in gastrulation by Veg1 endoderm cells. We also analyzed the contribution of cell divisions to elongation. Endoderm cells in Lytechinus variagatus were determined to go through approximately one cell doubling during gastrulation. That doubling occurs without a net increase in cell mass, but the question remained as to whether oriented divisions might contribute to archenteron elongation. We learned that indeed there was a biased orientation of cell divisions along the plane of archenteron elongation, but when the impact of that bias was analyzed quantitatively, it contributed a maximum 15% to the total elongation of the gut.

CONCLUSIONS

The major driver of archenteron elongation in the sea urchin, Lytechinus variagatus, is directed movement of Veg2 endoderm cells as a narrowing column along the plane of elongation. The narrowing occurs as cells in the column converge as they migrate, so that the combination of migration and the angular convergence provide the major component of the lengthening. A minor contributor to elongation is oriented cell divisions that contribute to the lengthening but no more than about 15%.

Fitting structure to function in gene regulatory networks

Cascades of transcriptional regulation are the common source of the forward drive in all developmental systems. Increases in complexity and specificity of gene expression at successive stages are based on the collaboration of varied combinations of transcription factors already expressed in the cells to turn on new genes, and the logical relationships between the transcription factors acting and becoming newly expressed from stage to stage are best visualized as gene regulatory networks. However, gene regulatory networks used in different developmental contexts underlie processes that actually operate through different sets of rules, which affect the kinetics, synchronicity, and logical properties of individual network nodes. Contrasting early embryonic development in flies and sea urchins with adult mammalian hematopoietic development from stem cells, major differences are seen in transcription factor dosage dependence, the silencing or damping impacts of repression, and the impact of cellular regulatory history on the parts of the genome that are accessible to transcription factor action in a given cell type. These different features not only affect the kinds of models that can illuminate developmental mechanisms in the respective biological systems, but also reflect the evolutionary needs of these biological systems to optimize different aspects of development.

Parallel embryonic transcriptional programs evolve under distinct constraints and may enable morphological conservation amidst adaptation

Embryonic development evolves by balancing stringent morphological constraints with genetic and environmental variation. The design principle that allows developmental transcriptional programs to conserve embryonic morphology while adapting to environmental changes is still not fully understood. To address this fundamental challenge, we compare developmental transcriptomes of two sea urchin species, Paracentrotus lividus and Strongylocentrotus purpuratus, that shared a common ancestor about 40 million years ago and are geographically distant yet show similar morphology. We find that both developmental and housekeeping genes show highly dynamic and strongly conserved temporal expression patterns. The expression of other gene sets, including homeostasis and response genes, show divergent expression which could result from either evolutionary drift or adaptation to local environmental conditions. The interspecies correlations of developmental gene expressions are highest between morphologically similar developmental time points whereas the interspecies correlations of housekeeping gene expression are high between all the late zygotic time points. Relatedly, the position of the phylotypic stage varies between these two groups of genes: developmental gene expression shows highest conservation at mid-developmental stage, in agreement with the hourglass model while the conservation of housekeeping genes keeps increasing with developmental time. When all genes are combined, the relationship between conservation of gene expression and morphological similarity is partially masked by housekeeping genes and genes with diverged expression. Our study illustrates various transcriptional programs that coexist in the developing embryo and evolve under different constraints. Apparently, morphological constraints underlie the conservation of developmental gene expression while embryonic fitness requires the conservation of housekeeping gene expression and the species-specific adjustments of homeostasis gene expression. The distinct evolutionary forces acting on these transcriptional programs enable the conservation of similar body plans while allowing adaption.

Nuclear receptors connect progenitor transcription factors to cell cycle control

The specification and growth of organs is controlled simultaneously by networks of transcription factors. While the connection between these transcription factors with fate determinants is increasingly clear, how they establish the link with the cell cycle is far less understood. Here we investigate this link in the developing Drosophila eye, where two transcription factors, the MEIS1 homologue hth and the Zn-finger tsh, synergize to stimulate the proliferation of naïve eye progenitors. Experiments combining transcriptomics, open-chromatin profiling, motif analysis and functional assays indicate that these progenitor transcription factors exert a global regulation of the proliferation program. Rather than directly regulating cell cycle genes, they control proliferation through an intermediary layer of nuclear receptors of the ecdysone/estrogen-signaling pathway. This regulatory subnetwork between hth, tsh and nuclear receptors might be conserved from Drosophila to mammals, as we find a significant co-overexpression of their human homologues in specific cancer types.

Mouth development

A mouth is present in all animals, and comprises an opening from the outside into the oral cavity and the beginnings of the digestive tract to allow eating. This review focuses on the earliest steps in mouth formation. In the first half, we conclude that the mouth arose once during evolution. In all animals, the mouth forms from ectoderm and endoderm. A direct association of oral ectoderm and digestive endoderm is present even in triploblastic animals, and in chordates, this region is known as the extreme anterior domain (EAD). Further support for a single origin of the mouth is a conserved set of genes that form a 'mouth gene program' including foxA and otx2. In the second half of this review, we discuss steps involved in vertebrate mouth formation, using the frog Xenopus as a model. The vertebrate mouth derives from oral ectoderm from the anterior neural ridge, pharyngeal endoderm and cranial neural crest (NC). Vertebrates form a mouth by breaking through the body covering in a precise sequence including specification of EAD ectoderm and endoderm as well as NC, formation of a 'pre-mouth array,' basement membrane dissolution, stomodeum formation, and buccopharyngeal membrane perforation. In Xenopus, the EAD is also a craniofacial organizer that guides NC, while reciprocally, the NC signals to the EAD to elicit its morphogenesis into a pre-mouth array. Human mouth anomalies are prevalent and are affected by genetic and environmental factors, with understanding guided in part by use of animal models. WIREs Dev Biol 2017, 6:e275. doi: 10.1002/wdev.275 For further resources related to this article, please visit the WIREs website.

Stochastic modeling and numerical simulation of gene regulatory networks with protein bursting

Gene expression is inherently stochastic. Advanced single-cell microscopy techniques together with mathematical models for single gene expression led to important insights in elucidating the sources of intrinsic noise in prokaryotic and eukaryotic cells. In addition to the finite size effects due to low copy numbers, translational bursting is a dominant source of stochasticity in cell scenarios involving few short lived mRNA transcripts with high translational efficiency (as is typically the case for prokaryotes), causing protein synthesis to occur in random bursts. In the context of gene regulation cascades, the Chemical Master Equation (CME) governing gene expression has in general no closed form solution, and the accurate stochastic simulation of the dynamics of complex gene regulatory networks is a major computational challenge. The CME associated to a single gene self regulatory motif has been previously approximated by a one dimensional time dependent partial integral differential equation (PIDE). However, to the best of our knowledge, multidimensional versions for such PIDE have not been developed yet. Here we propose a multidimensional PIDE model for regulatory networks involving multiple genes with self and cross regulations (in which genes can be regulated by different transcription factors) derived as the continuous counterpart of a CME with jump process. The model offers a reliable description of systems with translational bursting. In order to provide an efficient numerical solution, we develop a semilagrangian method to discretize the differential part of the PIDE, combined with a composed trapezoidal quadrature formula to approximate the integral term. We apply the model and numerical method to study sustained stochastic oscillations and the development of competence, a particular case of transient differentiation attained by certain bacterial cells under stress conditions. We found that the resulting probability distributions are distinguishable from those characteristic of other transient differentiation processes. In this way, they can be employed as markers or signatures that identify such phenomena from bacterial population experimental data, for instance. The computational efficiency of the semilagrangian method makes it suitable for purposes like model identification and parameter estimation from experimental data or, in combination with optimization routines, the design of gene regulatory networks under molecular noise.

Nodal and BMP expression during the transition to pentamery in the sea urchin Heliocidaris erythrogramma: insights into patterning the enigmatic echinoderm body plan

BACKGROUND

The molecular mechanisms underlying the development of the unusual echinoderm pentameral body plan and their likeness to mechanisms underlying the development of the bilateral plans of other deuterostomes are of interest in tracing body plan evolution. In this first study of the spatial expression of genes associated with Nodal and BMP2/4 signalling during the transition to pentamery in sea urchins, we investigate Heliocidaris erythrogramma, a species that provides access to the developing adult rudiment within days of fertilization.

RESULTS

BMP2/4, and the putative downstream genes, Six1/2, Eya, Tbx2/3 and Msx were expressed in the earliest morphological manifestation of pentamery during development, the five hydrocoele lobes. The formation of the vestibular ectoderm, the specialized region overlying the left coelom that forms adult ectoderm, involved the expression of putative Nodal target genes Chordin, Gsc and BMP2/4 and putative BMP2/4 target genes Dlx, Msx and Tbx. The expression of Nodal, Lefty and Pitx2 in the right ectoderm, and Pitx2 in the right coelom, was as previously observed in other sea urchins.

CONCLUSION

That genes associated with Nodal and BMP2/4 signalling are expressed in the hydrocoele lobes, indicates that they have a role in the developmental transition to pentamery, contributing to our understanding of how the most unusual body plan in the Bilateria may have evolved. We suggest that the Nodal and BMP2/4 signalling cascades might have been duplicated or split during the evolution to pentamery.

Identification and in silico modeling of enhancers reveals new features of the cardiac differentiation network

Developmental patterning and tissue formation are regulated through complex gene regulatory networks (GRNs) driven through the action of transcription factors (TFs) converging on enhancer elements. Here, as a point of entry to dissect the poorly defined GRN underlying cardiomyocyte differentiation, we apply an integrated approach to identify active enhancers and TFs involved in Drosophila heart development. The Drosophila heart consists of 104 cardiomyocytes, representing less than 0.5% of all cells in the embryo. By modifying BiTS-ChIP for rare cells, we examined H3K4me3 and H3K27ac chromatin landscapes to identify active promoters and enhancers specifically in cardiomyocytes. These in vivo data were complemented by a machine learning approach and extensive in vivo validation in transgenic embryos, which identified many new heart enhancers and their associated TF motifs. Our results implicate many new TFs in late stages of heart development, including Bagpipe, an Nkx3.2 ortholog, which we show is essential for differentiated heart function.

Itinerary profiling to analyze a large number of protein-folding trajectories

Understanding how proteins fold through a vast number of unfolded states is a major subject in the study of protein folding. Herein, we present itinerary profiling as a simple method to analyze molecular dynamics trajectories, and apply this method to Trp-cage. In itinerary profiling, structural clusters included in a trajectory are represented by a bit sequence, and a number of trajectories, as well as the structural clusters, can be compared and classified. As a consequence, the structural clusters that characterize the foldability of trajectories were able to be identified. The connections between the clusters were then illustrated as a network and the structural features of the clusters were examined. We found that in the true folding funnel, Trp-cage formed a left-handed main-chain topology and the Trp6 side-chain was located at the front of the main-chain ring, even in the initial unfolded states. In contrast, in the false folding funnel of the pseudo-native states, in which the Trp6 side-chain is upside down in the protein core, Trp-cage had a right-handed main-chain topology and the Trp side-chain was at the back. The initial topological partition, as determined by the main-chain handedness and the location of the Trp residue, predetermines Trp-cage foldability and the destination of the trajectory to the native state or the pseudo-native states.

Divergence of ectodermal and mesodermal gene regulatory network linkages in early development of sea urchins

Developmental gene regulatory networks (GRNs) are assemblages of gene regulatory interactions that direct ontogeny of animal body plans. Studies of GRNs operating in the early development of euechinoid sea urchins have revealed that little appreciable change has occurred since their divergence ∼90 million years ago (mya). These observations suggest that strong conservation of GRN architecture was maintained in early development of the sea urchin lineage. Testing whether this holds for all sea urchins necessitates comparative analyses of echinoid taxa that diverged deeper in geological time. Recent studies highlighted extensive divergence of skeletogenic mesoderm specification in the sister clade of euechinoids, the cidaroids, suggesting that comparative analyses of cidaroid GRN architecture may confer a greater understanding of the evolutionary dynamics of developmental GRNs. Here I report spatiotemporal patterning of 55 regulatory genes and perturbation analyses of key regulatory genes involved in euechinoid oral-aboral patterning of nonskeletogenic mesodermal and ectodermal domains in early development of the cidaroid Eucidaris tribuloides These results indicate that developmental GRNs directing mesodermal and ectodermal specification have undergone marked alterations since the divergence of cidaroids and euechinoids. Notably, statistical and clustering analyses of echinoid temporal gene expression datasets indicate that regulation of mesodermal genes has diverged more markedly than regulation of ectodermal genes. Although research on indirect-developing euechinoid sea urchins suggests strong conservation of GRN circuitry during early embryogenesis, this study indicates that since the divergence of cidaroids and euechinoids, developmental GRNs have undergone significant, cell type-biased alterations.

Transcriptional Control of Developmental Cell Behaviors

Tissue-specific transcription regulators emerged as key developmental control genes, which operate in the context of complex gene regulatory networks (GRNs) to coordinate progressive cell fate specification and tissue morphogenesis. We discuss how GRNs control the individual cell behaviors underlying complex morphogenetic events. Cell behaviors classically range from mesenchymal cell motility to cell shape changes in epithelial sheets. These behaviors emerge from the tissue-specific, multiscale integration of the local activities of universal and pleiotropic effectors, which underlie modular subcellular processes including cytoskeletal dynamics, cell-cell and cell-matrix adhesion, signaling, polarity, and vesicle trafficking. Extrinsic cues and intrinsic cell competence determine the subcellular spatiotemporal patterns of effector activities. GRNs influence most subcellular activities by controlling only a fraction of the effector-coding genes, which we argue is enriched in effectors involved in reading and processing the extrinsic cues to contextualize intrinsic subcellular processes and canalize developmental cell behaviors. The properties of the transcription-cell behavior interface have profound implications for evolution and disease.

Qualitative Dynamical Modelling Can Formally Explain Mesoderm Specification and Predict Novel Developmental Phenotypes

Given the complexity of developmental networks, it is often difficult to predict the effect of genetic perturbations, even within coding genes. Regulatory factors generally have pleiotropic effects, exhibit partially redundant roles, and regulate highly interconnected pathways with ample cross-talk. Here, we delineate a logical model encompassing 48 components and 82 regulatory interactions involved in mesoderm specification during Drosophila development, thereby providing a formal integration of all available genetic information from the literature. The four main tissues derived from mesoderm correspond to alternative stable states. We demonstrate that the model can predict known mutant phenotypes and use it to systematically predict the effects of over 300 new, often non-intuitive, loss- and gain-of-function mutations, and combinations thereof. We further validated several novel predictions experimentally, thereby demonstrating the robustness of model. Logical modelling can thus contribute to formally explain and predict regulatory outcomes underlying cell fate decisions.

Molluscan models: Crepidula fornicata

Gastropod snails in the genus Crepidula have emerged as model systems for studying a metazoan super clade, the Spiralia. Recent work on one species in particular, Crepidula fornicata, has produced high-resolution cell lineage fate maps, details of morphogenetic events during gastrulation, key insights into the molecular underpinnings of early development, and the first demonstration of CRISPR/Cas9 genome editing in the Spiralia. Furthermore, invasive species of Crepidula are a significant ecological threat, while one of these, C. fornicata, is also being harvested for food. This review highlights progress towards developing these animals as models for evolutionary, developmental, and ecological studies. Such studies have contributed greatly to our understanding of biology in a major clade of bilaterians. This information may also help us to control and cultivate these snails.

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.

Mature maternal mRNAs are longer than zygotic ones and have complex degradation kinetics in sea urchin

Early in embryogenesis, maternally deposited transcripts are degraded and new zygotic transcripts are generated during the maternal to zygotic transition. Recent works have shown that early zygotic transcripts are short compared to maternal transcripts, in zebrafish and Drosophila species. The reduced zygotic transcript length was attributed to the short cell cycle in these organisms that prevents the transcription of long primary transcripts (intron delay). Here we study the length of maternal mRNAs and their degradation kinetics in two sea urchin species to further the understanding of maternal gene usage and processing. Early zygotic primary transcripts and mRNAs are shorter than maternal ones in the sea urchin, Strongylocentrotus purpuratus. Yet, while primary transcripts length increases when cell cycle lengthens, typical for intron delay, the relatively short length of zygotic mRNAs is consistent. The enhanced mRNA length is due to significantly longer maternal open reading frames and 3'UTRs compared to the zygotic lengths, a ratio that does not change with developmental time. This implies unique usage of both coding sequences and regulatory information in the maternal stage compared to the zygotic stages. We extracted the half-lifetimes due to maternal and zygotic degradation mechanisms from high-density time course of a set of maternal mRNAs in Paracentrotus lividus. The degradation rates due to maternal and zygotic degradation mechanisms are not correlated, indicating that these mechanisms are independent and relay on different regulatory information. Our studies illuminate specific structural and kinetic properties of sea urchin maternal mRNAs that might be broadly shared by other organisms.

Eric Davidson: Steps to a gene regulatory network for development

Eric Harris Davidson was a unique and creative intellectual force who grappled with the diversity of developmental processes used by animal embryos and wrestled them into an intelligible set of principles, then spent his life translating these process elements into molecularly definable terms through the architecture of gene regulatory networks. He took speculative risks in his theoretical writing but ran a highly organized, rigorous experimental program that yielded an unprecedentedly full characterization of a developing organism. His writings created logical order and a framework for mechanism from the complex phenomena at the heart of advanced multicellular organism development. This is a reminiscence of intellectual currents in his work as observed by the author through the last 30-35 years of Davidson's life.

Tunicates: exploring the sea shores and roaming the open ocean. A tribute to Thomas Huxley

This review is a tribute to the remarkable contributions of Thomas Huxley to the biology of tunicates, the likely sister group of vertebrates. In 1851, the great biologist and philosopher published two landmark papers on pelagic tunicates in the Philosophical Transactions of the Royal Society. They were dedicated to the description of the adult anatomy and life cycle of thaliaceans and appendicularians, the pelagic relatives of ascidians. In the first part of this review, we discuss the novel anatomical observations and evolutionary hypotheses made by Huxley, which would have a lasting influence on tunicate biology. We also briefly comment on the more philosophical reflections of Huxley on individuality. In the second part, we stress the originality and relevance of past and future studies of tunicates in the resolution of major biological issues. In particular, we focus on the complex relationship between genotype and phenotype and the phenomenon of developmental system drift. We propose that more than 150 years after Huxley's papers, tunicate embryos are still worth studying in their own right, independently of their evolutionary proximity to vertebrates, as they provide original and crucial insights into the process of animal evolution. Tunicates are still at the forefront of biological research.

Comparative Developmental Transcriptomics Reveals Rewiring of a Highly Conserved Gene Regulatory Network during a Major Life History Switch in the Sea Urchin Genus Heliocidaris

The ecologically significant shift in developmental strategy from planktotrophic (feeding) to lecithotrophic (nonfeeding) development in the sea urchin genus Heliocidaris is one of the most comprehensively studied life history transitions in any animal. Although the evolution of lecithotrophy involved substantial changes to larval development and morphology, it is not known to what extent changes in gene expression underlie the developmental differences between species, nor do we understand how these changes evolved within the context of the well-defined gene regulatory network (GRN) underlying sea urchin development. To address these questions, we used RNA-seq to measure expression dynamics across development in three species: the lecithotroph Heliocidaris erythrogramma, the closely related planktotroph H. tuberculata, and an outgroup planktotroph Lytechinus variegatus. Using well-established statistical methods, we developed a novel framework for identifying, quantifying, and polarizing evolutionary changes in gene expression profiles across the transcriptome and within the GRN. We found that major changes in gene expression profiles were more numerous during the evolution of lecithotrophy than during the persistence of planktotrophy, and that genes with derived expression profiles in the lecithotroph displayed specific characteristics as a group that are consistent with the dramatically altered developmental program in this species. Compared to the transcriptome, changes in gene expression profiles within the GRN were even more pronounced in the lecithotroph. We found evidence for conservation and likely divergence of particular GRN regulatory interactions in the lecithotroph, as well as significant changes in the expression of genes with known roles in larval skeletogenesis. We further use coexpression analysis to identify genes of unknown function that may contribute to both conserved and derived developmental traits between species. Collectively, our results indicate that distinct evolutionary processes operate on gene expression during periods of life history conservation and periods of life history divergence, and that this contrast is even more pronounced within the GRN than across the transcriptome as a whole.

Developmental gene regulatory networks in sea urchins and what we can learn from them

Sea urchin embryos begin zygotic transcription shortly after the egg is fertilized.  Throughout the cleavage stages a series of transcription factors are activated and, along with signaling through a number of pathways, at least 15 different cell types are specified by the beginning of gastrulation.  Experimentally, perturbation of contributing transcription factors, signals and receptors and their molecular consequences enabled the assembly of an extensive gene regulatory network model.  That effort, pioneered and led by Eric Davidson and his laboratory, with many additional insights provided by other laboratories, provided the sea urchin community with a valuable resource.  Here we describe the approaches used to enable the assembly of an advanced gene regulatory network model describing molecular diversification during early development.  We then provide examples to show how a relatively advanced authenticated network can be used as a tool for discovery of how diverse developmental mechanisms are controlled and work.

An experimentally validated network of nine haematopoietic transcription factors reveals mechanisms of cell state stability

Transcription factor (TF) networks determine cell-type identity by establishing and maintaining lineage-specific expression profiles, yet reconstruction of mammalian regulatory network models has been hampered by a lack of comprehensive functional validation of regulatory interactions. Here, we report comprehensive ChIP-Seq, transgenic and reporter gene experimental data that have allowed us to construct an experimentally validated regulatory network model for haematopoietic stem/progenitor cells (HSPCs). Model simulation coupled with subsequent experimental validation using single cell expression profiling revealed potential mechanisms for cell state stabilisation, and also how a leukaemogenic TF fusion protein perturbs key HSPC regulators. The approach presented here should help to improve our understanding of both normal physiological and disease processes.

Robustness and Accuracy in Sea Urchin Developmental Gene Regulatory Networks

Developmental gene regulatory networks robustly control the timely activation of regulatory and differentiation genes. The structure of these networks underlies their capacity to buffer intrinsic and extrinsic noise and maintain embryonic morphology. Here I illustrate how the use of specific architectures by the sea urchin developmental regulatory networks enables the robust control of cell fate decisions. The Wnt-βcatenin signaling pathway patterns the primary embryonic axis while the BMP signaling pathway patterns the secondary embryonic axis in the sea urchin embryo and across bilateria. Interestingly, in the sea urchin in both cases, the signaling pathway that defines the axis controls directly the expression of a set of downstream regulatory genes. I propose that this direct activation of a set of regulatory genes enables a uniform regulatory response and a clear cut cell fate decision in the endoderm and in the dorsal ectoderm. The specification of the mesodermal pigment cell lineage is activated by Delta signaling that initiates a triple positive feedback loop that locks down the pigment specification state. I propose that the use of compound positive feedback circuitry provides the endodermal cells enough time to turn off mesodermal genes and ensures correct mesoderm vs. endoderm fate decision. Thus, I argue that understanding the control properties of repeatedly used regulatory architectures illuminates their role in embryogenesis and provides possible explanations to their resistance to evolutionary change.

Quantifying collectivity

In biological function emerges from the interactions of components with only partially aligned interests. An example is the brain-a large aggregation of neurons capable of producing unitary, coherent output. A theory for how such aggregations produce coherent output remains elusive. A first question we might ask is how collective is the behavior of the components? Here we introduce two properties of collectivity and illustrate how these properties can be quantified using approaches from information theory and statistical physics. First, amplification quantifies the sensitivity of the large scale to information at the small scale and is related to the notion of criticality in statistical physics. Second, decomposability reveals the extent to which aggregate behavior is reducible to individual contributions or is the result of synergistic interactions among components forming larger subgroups. These measures facilitate identification of causally important components and subgroups that might be experimentally manipulated to study the evolution and controllability of biological circuits and their outputs.

Implications of Developmental Gene Regulatory Networks Inside and Outside Developmental Biology

The insight that the genomic control of developmental process is encoded in the form of gene regulatory networks has profound impacts on many areas of modern bioscience. Most importantly, it affects developmental biology itself, as it means that a causal understanding of development requires knowledge of the architecture of regulatory network interactions. Furthermore, it follows that functional changes in developmental gene regulatory networks have to be considered as a primary mechanism for evolutionary process. We here discuss some of the recent advances in gene regulatory network biology and how they have affected our current understanding of development, evolution, and regulatory genomics.

Ancestral state reconstruction by comparative analysis of a GRN kernel operating in echinoderms

Diverse sampling of organisms across the five major classes in the phylum Echinodermata is beginning to reveal much about the structure and function of gene regulatory networks (GRNs) in development and evolution. Sea urchins are the most studied clade within this phylum, and recent work suggests there has been dramatic rewiring at the top of the skeletogenic GRN along the lineage leading to extant members of the euechinoid sea urchins. Such rewiring likely accounts for some of the observed developmental differences between the two major subclasses of sea urchins-cidaroids and euechinoids. To address effects of topmost rewiring on downstream GRN events, we cloned four downstream regulatory genes within the skeletogenic GRN and surveyed their spatiotemporal expression patterns in the cidaroid Eucidaris tribuloides. We performed phylogenetic analyses with homologs from other non-vertebrate deuterostomes and characterized their spatiotemporal expression by quantitative polymerase chain reaction (qPCR) and whole-mount in situ hybridization (WMISH). Our data suggest the erg-hex-tgif subcircuit, a putative GRN kernel, exhibits a mesoderm-specific expression pattern early in Eucidaris development that is directly downstream of the initial mesodermal GRN circuitry. Comparative analysis of the expression of this subcircuit in four echinoderm taxa allowed robust ancestral state reconstruction, supporting hypotheses that its ancestral function was to stabilize the mesodermal regulatory state and that it has been co-opted and deployed as a unit in mesodermal subdomains in distantly diverged echinoderms. Importantly, our study supports the notion that GRN kernels exhibit structural and functional modularity, locking down and stabilizing clade-specific, embryonic regulatory states.

Sea Urchin Morphogenesis

In the sea urchin morphogenesis follows extensive molecular specification. The specification controls the many morphogenetic events and these, in turn, precede patterning steps that establish the larval body plan. To understand how the embryo is built it was necessary to understand those series of molecular steps. Here an example of the historical sequence of those discoveries is presented as it unfolded over the last 50 years, the years during which major progress in understanding development of many animals and plants was documented by CTDB. In sea urchin development a rich series of experimental studies first established many of the phenomenological components of skeletal morphogenesis and patterning without knowledge of the molecular components. The many discoveries of transcription factors, signals, and structural proteins that contribute to the shape of the endoskeleton of the sea urchin larva then followed as molecular tools became available. A number of transcription factors and signals were discovered that were necessary for specification, morphogenesis, and patterning. Perturbation of the transcription factors and signals provided the means for assembling models of the gene regulatory networks used for specification and controlled the subsequent morphogenetic events. The earlier experimental information informed perturbation experiments that asked how patterning worked. As a consequence it was learned that ectoderm provides a series of patterning signals to the skeletogenic cells and as a consequence the skeletogenic cells secrete a highly patterned skeleton based on their ability to genotypically decode the localized reception of several signals. We still do not understand the complexity of the signals received by the skeletogenic cells, nor do we understand in detail how the genotypic information shapes the secreted skeletal biomineral, but the current knowledge at least outlines the sequence of events and provides a useful template for future discoveries.

Quantitative developmental transcriptomes of the Mediterranean sea urchin Paracentrotus lividus

Embryonic development progresses through the timely activation of thousands of differentially activated genes. Quantitative developmental transcriptomes provide the means to relate global patterns of differentially expressed genes to the emerging body plans they generate. The sea urchin is one of the classic model systems for embryogenesis and the models of its developmental gene regulatory networks are of the most comprehensive of their kind. Thus, the sea urchin embryo is an excellent system for studies of its global developmental transcriptional profiles. Here we produced quantitative developmental transcriptomes of the sea urchin Paracentrotus lividus (P. lividus) at seven developmental stages from the fertilized egg to prism stage. We generated de-novo reference transcriptome and identified 29,817 genes that are expressed at this time period. We annotated and quantified gene expression at the different developmental stages and confirmed the reliability of the expression profiles by QPCR measurement of a subset of genes. The progression of embryo development is reflected in the observed global expression patterns and in our principle component analysis. Our study illuminates the rich patterns of gene expression that participate in sea urchin embryogenesis and provide an essential resource for further studies of the dynamic expression of P. lividus genes.

Sieve-based relation extraction of gene regulatory networks from biological literature

Background

Relation extraction is an essential procedure in literature mining. It focuses on extracting semantic relations between parts of text, called mentions. Biomedical literature includes an enormous amount of textual descriptions of biological entities, their interactions and results of related experiments. To extract them in an explicit, computer readable format, these relations were at first extracted manually from databases. Manual curation was later replaced with automatic or semi-automatic tools with natural language processing capabilities. The current challenge is the development of information extraction procedures that can directly infer more complex relational structures, such as gene regulatory networks.

Results

We develop a computational approach for extraction of gene regulatory networks from textual data. Our method is designed as a sieve-based system and uses linear-chain conditional random fields and rules for relation extraction. With this method we successfully extracted the sporulation gene regulation network in the bacterium Bacillus subtilis for the information extraction challenge at the BioNLP 2013 conference. To enable extraction of distant relations using first-order models, we transform the data into skip-mention sequences. We infer multiple models, each of which is able to extract different relationship types. Following the shared task, we conducted additional analysis using different system settings that resulted in reducing the reconstruction error of bacterial sporulation network from 0.73 to 0.68, measured as the slot error rate between the predicted and the reference network. We observe that all relation extraction sieves contribute to the predictive performance of the proposed approach. Also, features constructed by considering mention words and their prefixes and suffixes are the most important features for higher accuracy of extraction. Analysis of distances between different mention types in the text shows that our choice of transforming data into skip-mention sequences is appropriate for detecting relations between distant mentions.

Conclusions

Linear-chain conditional random fields, along with appropriate data transformations, can be efficiently used to extract relations. The sieve-based architecture simplifies the system as new sieves can be easily added or removed and each sieve can utilize the results of previous ones. Furthermore, sieves with conditional random fields can be trained on arbitrary text data and hence are applicable to broad range of relation extraction tasks and data domains.

De novo assembly and analysis of tissue-specific transcriptomes revealed the tissue-specific genes and profile of immunity from Strongylocentrotus intermedius

Strongylocentrotus intermedius is an important marine species in north China and Japan. Recent years, diseases are threating the sea urchin aquaculture industry seriously. To provide a genetic resource for S. intermedius as well as overview the immune-related genes of S. intermedius, we performed transcriptome sequencing of three cDNA libraries representing three tissues, coelomocytes, gut and peristomial membrane respectively. In total 138,421 contigs were assembled from all sequencing data. 96,764 contigs were annotated according to bioinformatics databases, including NT, nr, Swiss-Prot, KEGG, COG. 49,336 Contigs were annotated as CDS. In this study, we obtained 24,778 gene families from S. intermedius transcriptome. The gene expression analysis revealed that more genes were expressed in gut, more high expression level genes in coelomocytes when compared with other tissues. Specific expressed contigs in coelomocytes, gut, and peristomial membrane were 546, 1136, and 1012 respectively. Pathway analysis suggested 25, 17 and 36 potential specifically pathways may specific progressed in peristomial membrane, gut and coelomocytes respectively. Similarities and differences between S. intermedius and other echinoderms were analyzed. S. intermedius was more homology to Strongylocentrotus purpuratus than others sea urchin. Of 24,778 genes, 1074 genes are immune-related, immune genes were expressed with a higher level in coelomocytes than other tissues. Complement system may be the most important immune system in sea urchin. We also identified 2438 SSRs and 16,236 SNPs for S. intermedius. These results provide a transcriptome resource and foundation to study molecular mechanisms of sea urchin immune system.