EchinoBase: Tools for Echinoderm Genome Analyses

Echinobase: a resource to support the echinoderm research community

Echinobase (www.echinobase.org) is a model organism knowledgebase serving as a resource for the community that studies echinoderms, a phylum of marine invertebrates that includes sea urchins and sea stars. Echinoderms have been important experimental models for over 100 years and continue to make important contributions to environmental, evolutionary, and developmental studies, including research on developmental gene regulatory networks. As a centralized resource, Echinobase hosts genomes and collects functional genomic data, reagents, literature, and other information for the community. This third-generation site is based on the Xenbase knowledgebase design and utilizes gene-centric pages to minimize the time and effort required to access genomic information. Summary gene pages display gene symbols and names, functional data, links to the JBrowse genome browser, and orthology to other organisms and reagents, and tabs from the Summary gene page contain more detailed information concerning mRNAs, proteins, diseases, and protein-protein interactions. The gene pages also display 1:1 orthologs between the fully supported species Strongylocentrotus purpuratus (purple sea urchin), Lytechinus variegatus (green sea urchin), Patiria miniata (bat star), and Acanthaster planci (crown-of-thorns sea star). JBrowse tracks are available for visualization of functional genomic data from both fully supported species and the partially supported species Anneissia japonica (feather star), Asterias rubens (sugar star), and L. pictus (painted sea urchin). Echinobase serves a vital role by providing researchers with annotated genomes including orthology, functional genomic data aligned to the genomes, and curated reagents and data. The Echinoderm Anatomical Ontology provides a framework for standardizing developmental data across the phylum, and knowledgebase content is formatted to be findable, accessible, interoperable, and reusable by the research community.

Thyroid hormone membrane receptor binding and transcriptional regulation in the sea urchin Strongylocentrotus purpuratus

Thyroid hormones (THs) are small amino acid derived signaling molecules with broad physiological and developmental functions in animals. Specifically, their function in metamorphic development, ion regulation, angiogenesis and many others have been studied in detail in mammals and some other vertebrates. Despite extensive reports showing pharmacological responses of invertebrate species to THs, little is known about TH signaling mechanisms outside of vertebrates. Previous work in sea urchins suggests that non-genomic mechanisms are activated by TH ligands. Here we show that several THs bind to sea urchin (Strongylocentrotus purpuratus) cell membrane extracts and are displaced by ligands of RGD-binding integrins. A transcriptional analysis across sea urchin developmental stages shows activation of genomic and non-genomic pathways in response to TH exposure, suggesting that both pathways are activated by THs in sea urchin embryos and larvae. We also provide evidence associating TH regulation of gene expression with TH response elements in the genome. In ontogeny, we found more differentially expressed genes in older larvae compared to gastrula stages. In contrast to gastrula stages, the acceleration of skeletogenesis by thyroxine in older larvae is not fully inhibited by competitive ligands or inhibitors of the integrin membrane receptor pathway, suggesting that THs likely activate multiple pathways. Our data confirms a signaling function of THs in sea urchin development and suggests that both genomic and non-genomic mechanisms play a role, with genomic signaling being more prominent during later stages of larval development.

Twinkle twinkle brittle star: the draft genome of Ophioderma brevispinum (Echinodermata: Ophiuroidea) as a resource for regeneration research

Background

Echinoderms are established models in experimental and developmental biology, however genomic resources are still lacking for many species. Here, we present the draft genome of Ophioderma brevispinum, an emerging model organism in the field of regenerative biology. This new genomic resource provides a reference for experimental studies of regenerative mechanisms.

Results

We report a de novo nuclear genome assembly for the brittle star O. brevispinum and annotation facilitated by the transcriptome assembly. The final assembly is 2.68 Gb in length and contains 146,703 predicted protein-coding gene models. We also report a mitochondrial genome for this species, which is 15,831 bp in length, and contains 13 protein-coding, 22 tRNAs, and 2 rRNAs genes, respectively. In addition, 29 genes of the Notch signaling pathway are identified to illustrate the practical utility of the assembly for studies of regeneration.

Conclusions

The sequenced and annotated genome of O. brevispinum presented here provides the first such resource for an ophiuroid model species. Considering the remarkable regenerative capacity of this species, this genome will be an essential resource in future research efforts on molecular mechanisms regulating regeneration.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08750-y).

The arm of the starfish: The far-reaching applications of Patiria miniata as a model system in evolutionary, developmental, and regenerative biology

Many species of echinoderms have long been considered model research organisms in biology. Historically, much of this research has focused on the embryology of sea urchins and the use of their extensive gene regulatory networks as a tool to understand how the genome controls cell state specification and patterning. The establishment of Patiria miniata, the bat sea star, as a research organism has allowed us to expand on the concepts explored with sea urchins, viewing these genetic networks through a comparative lens, gaining great insight into the evolutionary mechanisms that shape developmental diversity. Extensive molecular tools have been developed in P. miniata, designed to explore gene expression dynamics and build gene regulatory networks. Echinoderms also have a robust set of bioinformatic and computational resources, centered around echinobase.org, an extensive database containing multiomic, developmental, and experimental resources for researchers. In addition to comparative evolutionary development, P. miniata is a promising system in its own right for studying whole body regeneration, metamorphosis and body plan development, as well as marine disease.

Regeneration of the larval sea star nervous system by wounding induced respecification to the sox2 lineage

The ability to restore lost body parts following traumatic injury is a fascinating area of biology that challenges current understanding of the ontogeny of differentiation. The origin of new cells needed to regenerate lost tissue, and whether they are pluripotent or have de- or trans-differentiated, remains one of the most important open questions . Additionally, it is not known whether developmental gene regulatory networks are reused or whether regeneration specific networks are deployed. Echinoderms, including sea stars, have extensive ability for regeneration, however, the technologies for obtaining transgenic echinoderms are limited and tracking cells involved in regeneration, and thus identifying the cellular sources and potencies has proven challenging. In this study, we develop new transgenic tools to follow the fate of populations of cells in the regenerating larva of the sea star Patiria miniata. We show that the larval serotonergic nervous system can regenerate following decapitation. Using a BAC-transgenesis approach we show that expression of the pan ectodermal marker, sox2, is induced in previously sox2 minus cells , even when cell division is inhibited. sox2+ cells give rise to new sox4+ neural precursors that then proceed along an embryonic neurogenesis pathway to reform the anterior nervous systems. sox2+ cells contribute to only neural and ectoderm lineages, indicating that these progenitors maintain their normal, embryonic lineage restriction. This indicates that sea star larval regeneration uses a combination of existing lineage restricted stem cells, as well as respecification of cells into neural lineages, and at least partial reuse of developmental GRNs to regenerate their nervous system.

Echinobase: leveraging an extant model organism database to build a knowledgebase supporting research on the genomics and biology of echinoderms

Echinobase (www.echinobase.org) is a third generation web resource supporting genomic research on echinoderms. The new version was built by cloning the mature Xenopus model organism knowledgebase, Xenbase, refactoring data ingestion pipelines and modifying the user interface to adapt to multispecies echinoderm content. This approach leveraged over 15 years of previous database and web application development to generate a new fully featured informatics resource in a single year. In addition to the software stack, Echinobase uses the private cloud and physical hosts that support Xenbase. Echinobase currently supports six echinoderm species, focused on those used for genomics, developmental biology and gene regulatory network analyses. Over 38 000 gene pages, 18 000 publications, new improved genome assemblies, JBrowse genome browser and BLAST + services are available and supported by the development of a new echinoderm anatomical ontology, uniformly applied formal gene nomenclature, and consistent orthology predictions. A novel feature of Echinobase is integrating support for multiple, disparate species. New genomes from the diverse echinoderm phylum will be added and supported as data becomes available. The common code development design of the integrated knowledgebases ensures parallel improvements as each resource evolves. This approach is widely applicable for developing new model organism informatics resources.

Developmental gene regulatory network connections predicted by machine learning from gene expression data alone

Gene regulatory network (GRN) inference can now take advantage of powerful machine learning algorithms to complement traditional experimental methods in building gene networks. However, the dynamical nature of embryonic development-representing the time-dependent interactions between thousands of transcription factors, signaling molecules, and effector genes-is one of the most challenging arenas for GRN prediction. In this work, we show that successful GRN predictions for a developmental network from gene expression data alone can be obtained with the Priors Enriched Absent Knowledge (PEAK) network inference algorithm. PEAK is a noise-robust method that models gene expression dynamics via ordinary differential equations and selects the best network based on information-theoretic criteria coupled with the machine learning algorithm Elastic Net. We test our GRN prediction methodology using two gene expression datasets for the purple sea urchin, Stronglyocentrotus purpuratus, and cross-check our results against existing GRN models that have been constructed and validated by over 30 years of experimental results. Our results find a remarkably high degree of sensitivity in identifying known gene interactions in the network (maximum 81.58%). We also generate novel predictions for interactions that have not yet been described, which provide a resource for researchers to use to further complete the sea urchin GRN. Published ChIPseq data and spatial co-expression analysis further support a subset of the top novel predictions. We conclude that GRN predictions that match known gene interactions can be produced using gene expression data alone from developmental time series experiments.

Polarized Dishevelled dissolution and reassembly drives embryonic axis specification in sea star oocytes

The organismal body axes that are formed during embryogenesis are intimately linked to intrinsic asymmetries established at the cellular scale in oocytes.1 However, the mechanisms that generate cellular asymmetries within the oocyte and then transduce that polarity to organismal scale body axes are poorly understood outside of select model organisms. Here, we report an axis-defining event in meiotic oocytes of the sea star Patiria miniata. Dishevelled (Dvl) is a cytoplasmic Wnt pathway effector required for axis development in diverse species,2-4 but the mechanisms governing its function and distribution remain poorly defined. Using time-lapse imaging, we find that Dvl localizes uniformly to puncta throughout the cell cortex in Prophase I-arrested oocytes but becomes enriched at the vegetal pole following meiotic resumption through a dissolution-reassembly mechanism. This process is driven by an initial disassembly phase of Dvl puncta, followed by selective reformation of Dvl assemblies at the vegetal pole. Rather than being driven by Wnt signaling, this localization behavior is coupled to meiotic cell cycle progression and influenced by Lamp1+ endosome association and Frizzled receptors pre-localized within the oocyte cortex. Our results reveal a cell cycle-linked mechanism by which maternal cellular polarity is transduced to the embryo through spatially regulated Dvl dynamics.

Regeneration in Echinoderms: Molecular Advancements

Which genes and gene signaling pathways mediate regenerative processes? In recent years, multiple studies, using a variety of animal models, have aimed to answer this question. Some answers have been obtained from transcriptomic and genomic studies where possible gene and gene pathway candidates thought to be involved in tissue and organ regeneration have been identified. Several of these studies have been done in echinoderms, an animal group that forms part of the deuterostomes along with vertebrates. Echinoderms, with their outstanding regenerative abilities, can provide important insights into the molecular basis of regeneration. Here we review the available data to determine the genes and signaling pathways that have been proposed to be involved in regenerative processes. Our analyses provide a curated list of genes and gene signaling pathways and match them with the different cellular processes of the regenerative response. In this way, the molecular basis of echinoderm regenerative potential is revealed, and is available for comparisons with other animal taxa.

Dynamically expressed single ELAV/Hu orthologue elavl2 of bees is required for learning and memory

Changes in gene expression are a hallmark of learning and memory consolidation. Little is known about how alternative mRNA processing, particularly abundant in neuron-specific genes, contributes to these processes. Prototype RNA binding proteins of the neuronally expressed ELAV/Hu family are candidates for roles in learning and memory, but their capacity to cross-regulate and take over each other's functions complicate substantiation of such links. Honey bees Apis mellifera have only one elav/Hu family gene elavl2, that has functionally diversified by increasing alternative splicing including an evolutionary conserved microexon. RNAi knockdown demonstrates that ELAVL2 is required for learning and memory in bees. ELAVL2 is dynamically expressed with altered alternative splicing and subcellular localization in mushroom bodies, but not in other brain regions. Expression and alternative splicing of elavl2 change during memory consolidation illustrating an alternative mRNA processing program as part of a local gene expression response underlying memory consolidation.

Classifying domain-specific text documents containing ambiguous keywords

A keyword-based search of comprehensive databases such as PubMed may return irrelevant papers, especially if the keywords are used in multiple fields of study. In such cases, domain experts (curators) need to verify the results and remove the irrelevant articles. Automating this filtering process will save time, but it has to be done well enough to ensure few relevant papers are rejected and few irrelevant papers are accepted. A good solution would be fast, work with the limited amount of data freely available (full paper body may be missing), handle ambiguous keywords and be as domain-neutral as possible. In this paper, we evaluate a number of classification algorithms for identifying a domain-specific set of papers about echinoderm species and show that the resulting tool satisfies most of the abovementioned requirements. Echinoderms consist of a number of very different organisms, including brittle stars, sea stars (starfish), sea urchins and sea cucumbers. While their taxonomic identifiers are specific, the common names are used in many other contexts, creating ambiguity and making a keyword search prone to error. We try classifiers using Linear, Naïve Bayes, Nearest Neighbor, Tree, SVM, Bagging, AdaBoost and Neural Network learning models and compare their performance. We show how effective the resulting classifiers are in filtering irrelevant articles returned from PubMed. The methodology used is more dependent on the good selection of training data and is a practical solution that can be applied to other fields of study facing similar challenges.

Database URL The code and date reported in this paper are freely available at http://xenbaseturbofrog.org/pub/Text-Topic-Classifier/

The Development and Neuronal Complexity of Bipinnaria Larvae of the Sea Star Asterias rubens

Free-swimming planktonic larvae are a key stage in the development of many marine phyla, and studies of these organisms have contributed to our understanding of major genetic and evolutionary processes. Although transitory, these larvae often attain a remarkable degree of tissue complexity, with well-defined musculature and nervous systems. Among the best studied are larvae belonging to the phylum Echinodermata, but with work largely focused on the pluteus larvae of sea urchins (class Echinoidea). The greatest diversity of larval strategies among echinoderms is found in the class Asteroidea (sea stars), organisms that are rapidly emerging as experimental systems for genetic and developmental studies. However, the bipinnaria larvae of sea stars have only been studied in detail in a small number of species and although they have been relatively well described neuro-anatomically, they are poorly understood neurochemically. Here, we have analyzed embryonic development and bipinnaria larval anatomy in the common North Atlantic sea star Asterias rubens, using a variety of staining methods in combination with confocal microscopy. Importantly, the chemical complexity of the nervous system of bipinnaria larvae was revealed through use of a diverse set of antibodies, with identification of at least three centers of differing neurochemical signature within the previously described nervous system: the anterior apical organ, oral region, and ciliary bands. Furthermore, the anatomy of the musculature and sites of cell division in bipinnaria larvae was analyzed. Comparisons of developmental progression and molecular anatomy across the Echinodermata provided a basis for hypotheses on the shared evolutionary and developmental processes that have shaped this group of animals. We conclude that bipinnaria larvae appear to be remarkably conserved across ∼200 million years of evolutionary time and may represent a strong evolutionary and/or developmental constraint on species utilizing this larval strategy.

A nomenclature for echinoderm genes

Echinoderm embryos and larvae are prominent experimental model systems for studying developmental mechanisms. High-quality, assembled, annotated genome sequences are now available for several echinoderm species, including representatives from most classes. The increased availability of these data necessitates the development of a nomenclature that assigns universally interpretable gene symbols to echinoderm genes to facilitate cross-species comparisons of gene functions, both within echinoderms and across other phyla. This paper describes the implementation of an improved set of echinoderm gene nomenclature guidelines that both communicates meaningful orthology information in protein-coding gene symbols and names and establishes continuity with nomenclatures developed for major vertebrate model organisms, including humans. Differences between the echinoderm gene nomenclature guidelines and vertebrate guidelines are examined and explained. This nomenclature incorporates novel solutions to allow for several types of orthologous relationships, including the single echinoderm genes with multiple vertebrate co-orthologs that result from whole-genome-duplication events. The current version of the Echinoderm Gene Nomenclature Guidelines can be found at https://www.echinobase.org/gene/static/geneNomenclature.jsp

Database URL https://www.echinobase.org/

The tweety Gene Family: From Embryo to Disease

The tweety genes encode gated chloride channels that are found in animals, plants, and even simple eukaryotes, signifying their deep evolutionary origin. In vertebrates, the tweety gene family is highly conserved and consists of three members—ttyh1, ttyh2, and ttyh3—that are important for the regulation of cell volume. While research has elucidated potential physiological functions of ttyh1 in neural stem cell maintenance, proliferation, and filopodia formation during neural development, the roles of ttyh2 and ttyh3 are less characterized, though their expression patterns during embryonic and fetal development suggest potential roles in the development of a wide range of tissues including a role in the immune system in response to pathogen-associated molecules. Additionally, members of the tweety gene family have been implicated in various pathologies including cancers, particularly pediatric brain tumors, and neurodegenerative diseases such as Alzheimer’s and Parkinson’s disease. Here, we review the current state of research using information from published articles and open-source databases on the tweety gene family with regard to its structure, evolution, expression during development and adulthood, biochemical and cellular functions, and role in human disease. We also identify promising areas for further research to advance our understanding of this important, yet still understudied, family of genes.

Integration of 1:1 orthology maps and updated datasets into Echinobase

Echinobase (https://echinobase.org) is a central online platform that generates, manages and hosts genomic data relevant to echinoderm research. While the resource primarily serves the echinoderm research community, the recent release of an excellent quality genome for the frequently studied purple sea urchin (Strongylocentrotus purpuratus genome, v5.0) has provided an opportunity to adapt to the needs of a broader research community across other model systems. To this end, establishing pipelines to identify orthologous genes between echinoderms and other species has become a priority in many contexts including nomenclature, linking to data in other model organisms, and in internal functionality where data gathered in one hosted species can be associated with genes in other hosted echinoderms. This paper describes the orthology pipelines currently employed by Echinobase and how orthology data are processed to yield 1:1 ortholog mappings between a variety of echinoderms and other model taxa. We also describe functions of interest that have recently been included on the resource, including an updated developmental time course for S.purpuratus, and additional tracks for genome browsing. These data enhancements will increase the accessibility of the resource to non-echinoderm researchers and simultaneously expand the data quality and quantity available to core Echinobase users. Database URL: https://echinobase.org.

CrustyBase: an interactive online database for crustacean transcriptomes

Transcriptome sequencing has opened the field of genomics to a wide variety of researchers, owing to its efficiency, applicability across species and ability to quantify gene expression. The resulting datasets are a rich source of information that can be mined for many years into the future, with each dataset providing a unique angle on a specific context in biology. Maintaining accessibility to this accumulation of data presents quite a challenge for researchers.The primary focus of conventional genomics databases is the storage, navigation and interpretation of sequence data, which is typically classified down to the level of a species or individual. The addition of expression data adds a new dimension to this paradigm - the sampling context. Does gene expression describe different tissues, a temporal distribution or an experimental treatment? These data not only describe an individual, but the biological context surrounding that individual. The structure and utility of a transcriptome database must therefore reflect these attributes. We present an online database which has been designed to maximise the accessibility of crustacean transcriptome data by providing intuitive navigation within and between datasets and instant visualization of gene expression and protein structure.The site is accessible at https://crustybase.org and currently holds 10 datasets from a range of crustacean species. It also allows for upload of novel transcriptome datasets through a simple web interface, allowing the research community to contribute their own data to a pool of shared knowledge.

Maternal control of early patterning in sea urchin embryos

Sea urchin development has been studied extensively for more than a century and considered regulative since the first experimental evidence. Further investigations have repeatedly supported this standpoint by revealing the presence of inductive mechanisms that alter cell fate decisions at early cleavage stages and flexibility of development in response to environmental conditions. Some features indicate that sea urchin development is not completely regulative, but actually includes determinative events. In 16-cell embryos, mesomeres and macromeres represent multipotency, while the cell fate of most vegetal micromeres is restricted. It is known that the mature sea urchin eggs are polarized by the asymmetrical distribution of some maternal mRNAs and proteins. Spatially-distributed maternal factors are necessary for the orientation of the primary animal-vegetal axis, which is established by both maternal and zygotic mechanisms later in development. The secondary dorsal-ventral axis is conditionally specified later in development. Dorsal-ventral polarity is very liable during the early cleavages, though more recent data argue that its direction may be oriented by maternal asymmetry. In this review, we focus on the role of maternal factors in initial embryonic patterning during the first cleavages of sea urchin embryos before activation of the embryonic genome.

Drosophila SLC22 Orthologs Related to OATs, OCTs, and OCTNs Regulate Development and Responsiveness to Oxidative Stress

The SLC22 family of transporters is widely expressed, evolutionarily conserved, and plays a major role in regulating homeostasis by transporting small organic molecules such as metabolites, signaling molecules, and antioxidants. Analysis of transporters in fruit flies provides a simple yet orthologous platform to study the endogenous function of drug transporters in vivo. Evolutionary analysis of Drosophila melanogaster putative SLC22 orthologs reveals that, while many of the 25 SLC22 fruit fly orthologs do not fall within previously established SLC22 subclades, at least four members appear orthologous to mammalian SLC22 members (SLC22A16:CG6356, SLC22A15:CG7458, CG7442 and SLC22A18:CG3168). We functionally evaluated the role of SLC22 transporters in Drosophila melanogaster by knocking down 14 of these genes. Three putative SLC22 ortholog knockdowns-CG3168, CG6356, and CG7442/SLC22A-did not undergo eclosion and were lethal at the pupa stage, indicating the developmental importance of these genes. Additionally, knocking down four SLC22 members increased resistance to oxidative stress via paraquat testing (CG4630: p < 0.05, CG6006: p < 0.05, CG6126: p < 0.01 and CG16727: p < 0.05). Consistent with recent evidence that SLC22 is central to a Remote Sensing and Signaling Network (RSSN) involved in signaling and metabolism, these phenotypes support a key role for SLC22 in handling reactive oxygen species.

Ion channels and signaling pathways used in the fast polyspermy block

Fertilization of an egg by multiple sperms, polyspermy, is lethal to most sexually reproducing species. To combat the entry of additional sperm into already fertilized eggs, organisms have developed various polyspermy blocks. One such barrier, the fast polyspermy block, uses a fertilization-activated depolarization of the egg membrane to electrically inhibit supernumerary sperm from entering the egg. The fast block is commonly used by eggs of oviparous animals with external fertilization. In this review, we discuss the history of the fast block discovery, as well as general features shared by all organisms that use this polyspermy block. Given the diversity of habitats of external fertilizers, the fine details of the fast block-signaling pathways differ drastically between species, including the identity of the depolarizing ions. We highlight the known molecular mediators of these signaling pathways in amphibians and echinoderms, with a fine focus on ion channels that signal these fertilization-evoked depolarizations. We also discuss the investigation for a fast polyspermy block in mammals and teleost fish, and we outline potential fast block triggers. Since the first electrical recordings made on eggs in the 1950s, the fields of developmental biology and electrophysiology have substantially matured, and yet we are only now beginning to discern the intricate molecular mechanisms regulating the fast block to polyspermy.

The role of thyroglobulin in thyroid hormonogenesis

In humans, the thyroid hormones T₃ and T₄ are synthesized in the thyroid gland in a process that crucially involves the iodoglycoprotein thyroglobulin. The overall structure of thyroglobulin is conserved in all vertebrates. Upon thyroglobulin delivery from thyrocytes to the follicular lumen of the thyroid gland via the secretory pathway, multiple tyrosine residues can become iodinated to form mono-iodotyrosine (MIT) and/or di-iodotyrosine (DIT); however, selective tyrosine residues lead to preferential formation of T₄ and T₃ at distinct sites. T₄ formation involves oxidative coupling between two DIT side chains, and de novo T₃ formation involves coupling between an MIT donor and a DIT acceptor. Thyroid hormone synthesis is stimulated by TSH activating its receptor (TSHR), which upregulates the activity of many thyroid gene products involved in hormonogenesis. Additionally, TSH regulates post-translational changes in thyroglobulin that selectively enhance its capacity for T₃ formation - this process is important in iodide deficiency and in Graves disease. 167 different mutations, many of which are newly discovered, are now known to exist in TG (encoding human thyroglobulin) that can lead to defective thyroid hormone synthesis, resulting in congenital hypothyroidism.

Techniques for analyzing gene expression using BAC-based reporter constructs

To characterize the complex regulatory control of gene expression using fluorescent protein reporters, it is often necessary to analyze large genomic regions. Bacteria artificial chromosome (BAC) vectors, which are able to support DNA fragments of up to 300kb, provide stable platforms for experimental manipulation. Using phage-based systems of homologous recombination, BACs can be efficiently engineered for a variety of aims. These include expressing fluorescent proteins to delineate gene expression boundaries using high-resolution, in vivo microscopy, tracing cell lineages using stable fluorescent proteins, perturbing endogenous protein function by expressing dominant negative forms, interfering with development by mis-expressing transcription factors, and identifying regulatory regions through deletion analysis. Here, we present a series of protocols for identifying BAC clones that contain genes of interest, modifying BACs for use as reporter constructs, and preparing BAC DNA for microinjection into fertilized eggs. Although the protocols here are tailored for use in echinoderm embryonic and larval stages, these methods are easily adaptable for use in other transgenic systems. As fluorescent protein technology continues to expand, so do the potential applications for recombinant BACs.

Spatial and temporal patterns of gene expression during neurogenesis in the sea urchin Lytechinus variegatus

BACKGROUND

The sea urchin is a basal deuterostome that is more closely related to vertebrates than many organisms traditionally used to study neurogenesis. This phylogenetic position means that the sea urchin can provide insights into the evolution of the nervous system by helping resolve which developmental processes are deuterostome innovations, which are innovations in other clades, and which are ancestral. However, the nervous system of echinoderms is one of the least understood of all major metazoan phyla. To gain insights into echinoderm neurogenesis, spatial and temporal gene expression data are essential. Then, functional data will enable the building of a detailed gene regulatory network for neurogenesis in the sea urchin that can be compared across metazoans to resolve questions about how nervous systems evolved.

RESULTS

Here, we analyze spatiotemporal gene expression during sea urchin neurogenesis for genes that have been shown to be neurogenic in one or more species. We report the expression of 21 genes expressed in areas of neurogenesis in the sea urchin embryo from blastula stage (just before neural progenitors begin their specification sequence) through pluteus larval stage (when much of the nervous system has been patterned). Among those 21 gene expression patterns, we report expression of 11 transcription factors and 2 axon guidance genes, each expressed in discrete domains in the neuroectoderm or in the endoderm. Most of these genes are expressed in and around the ciliary band. Some including the transcription factors Lv-mbx, Lv-dmrt, Lv-islet, and Lv-atbf1, the nuclear protein Lv-prohibitin, and the guidance molecule Lv-semaa are expressed in the endoderm where they are presumably involved in neurogenesis in the gut.

CONCLUSIONS

This study builds a foundation to study how neurons are specified and evolved by analyzing spatial and temporal gene expression during neurogenesis in a basal deuterostome. With these expression patterns, we will be able to understand what genes are required for neural development in the sea urchin. These data can be used as a starting point to (1) build a spatial gene regulatory network for sea urchin neurogenesis, (2) identify how subtypes of neurons are specified, (3) perform comparative studies with the sea urchin, protostome, and vertebrate organisms.

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.