The Transcriptome of the Sea Urchin Embryo

Expression Regulation Mechanisms of Sea Urchin (Strongylocentrotus intermedius) Under the High Temperature: New Evidence for the miRNA-mRNA Interaction Involvement

In the context of global warming and continuous high temperatures in the northern part of China during summer, the mortality rate of our main breeding species, Strongylocentrotus intermedius, reached 80% in 2020. How sea urchins respond to high temperatures is of great concern to academia and industry. In this study, we examined the antioxidant enzyme activities of different color tube-footed sea urchins under heat stress and compared their transcriptome and microRNA (miRNA) profiles using RNA-Seq. The results showed that the antioxidant enzyme activities of sea urchins were altered by thermal stress, and the changes in peroxidase activities of red tube-footed sea urchins were particularly significant. Investigations revealed that 1,079 differentially expressed genes (DEGs), 11 DE miRNAs, and 104 “DE miRNA-DEG” pairs in total were detected in sea urchins under high temperature stress. Several mRNA and miRNAs were significantly changed (e.g. HSP70, DnaJ11, HYAL, CALR, miR-184-p5, miR-92a, miR-92c, and miR-124-p5), suggesting these genes and miRNAs exerted important functions in response to high temperature. At the transcriptional level, red tube-footed sea urchins were found to be more sensitive to high temperature and could respond to high temperature rapidly. DE miRNA-mRNA network showed that miR-92b-3p and PC-5p-7420 were the most corresponding miRNAs. Five mRNAs (DnaJ11, SAR1B, CALR, HYOU1, TUBA) may be potential markers of sea urchin response to high temperature. Possible interaction between miRNA-mRNA could be linked to protein folding in the endoplasmic reticulum, Phagosomes, and calcium transport. This study provides a theoretical basis for the molecular mechanism of sea urchin heat tolerance and information that will aid in the selection and breeding of sea urchins with high temperature tolerance.

The developmental transcriptome for Lytechinus variegatus exhibits temporally punctuated gene expression changes

Embryonic development is arguably the most complex process an organism undergoes during its lifetime, and understanding this complexity is best approached with a systems-level perspective. The sea urchin has become a highly valuable model organism for understanding developmental specification, morphogenesis, and evolution. As a non-chordate deuterostome, the sea urchin occupies an important evolutionary niche between protostomes and vertebrates. Lytechinus variegatus (Lv) is an Atlantic species that has been well studied, and which has provided important insights into signal transduction, patterning, and morphogenetic changes during embryonic and larval development. The Pacific species, Strongylocentrotus purpuratus (Sp), is another well-studied sea urchin, particularly for gene regulatory networks (GRNs) and cis-regulatory analyses. A well-annotated genome and transcriptome for Sp are available, but similar resources have not been developed for Lv. Here, we provide an analysis of the Lv transcriptome at 11 timepoints during embryonic and larval development. Temporal analysis suggests that the gene regulatory networks that underlie specification are well-conserved among sea urchin species. We show that the major transitions in variation of embryonic transcription divide the developmental time series into four distinct, temporally sequential phases. Our work shows that sea urchin development occurs via sequential intervals of relatively stable gene expression states that are punctuated by abrupt transitions.

Transcriptome analysis to characterize the genes related to gonad growth and fatty acid metabolism in the sea urchin Strongylocentrotus intermedius

BACKGROUND

Sea urchin gonads of both sexes, commonly termed "roe", are highly valued seafood delicacies, and Strongylocentrotus intermedius is considered one of the tastiest sea urchins. In order to produce high-quality gonads for consumption and clarify the mechanism of gonad growth and development of the sea urchin, more genetic information, especially at the transcriptome level, is needed.

OBJECTIVE

A more thorough understanding of sea urchin gonad growth and development in both sexes could enable regulation of these processes at several stages with the aim of suppressing gametogenesis in order to produce high-quality gonads for consumption.

METHODS

The adult sea urchins S. intermedius were cultured for 3 months, and were sampled for the gonadal transcriptome analysis which has been performed on the RNAs of three male and female adults of S. intermedius in each gonad development stage.

RESULTS

Illumina sequencing raw sequence data was deposited in the NCBI Sequence Read Archive (SRA) database (PRJNA532998). It generated 560,196,356 raw reads and 548,956,944 clean reads were acquired, which were assembled into 107,850 transcripts with 44,124 genes. Comparative analysis showed the differentially expressed genes (DEGs) from 114 to 2566. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were used to determine the functional significance of these DEGs. We have selected 9 genes related to growth and 12 genes related to fatty acid biosynthesis and metabolism in sea urchin gonads.

CONCLUSION

These data for sea urchins were intended to provide markers for gonad growth and development that can be accumulated for use in aquaculture applications.

Glycosylation Fosters Interactions between Model Sea Urchin Spicule Matrix Proteins. Implications for Embryonic Spiculogenesis and Biomineralization

The formation of embryonic mineralized skeletal elements (spicules) in the sea urchin requires the participation of proteins, many of which may interact with one another and assist in the creation of an extracellular matrix wherein mineral formation takes place. To probe this, we created a sea urchin spicule recombinant model protein pair system wherein we tested the interactions between two major spicule proteins, SpSM50 and the glycoprotein, SpSM30B/C. Both proteins are strong hydrogelators that manipulate early and later events in mineral formation. We discovered that the anionic glycan moieties of SpSM30B/C are required for interaction with the SpSM50 protein and that these interactions are Ca(II)-independent. In addition, when these proteins form a complex, they create hybrid hydrogel particles that are physically distinct from their individual counterparts. Thus, glycan-mediated interactions play an important role in in vitro spicule protein assembly and most likely within the spicule itself.

The developmental transcriptomes of two sea biscuit species with differing larval types

BACKGROUND

Larval developmental patterns are extremely varied both between and within phyla, however the genetic mechanisms leading to this diversification are poorly understood. We assembled and compared the developmental transcriptomes for two sea biscuit species (Echinodermata: Echinoidea) with differing patterns of larval development, to provide a resource for investigating the evolution of alternate life cycles. One species (Clypeaster subdepressus) develops via an obligately feeding larva which metamorphoses 3-4 weeks after fertilization; the other (Clypeaster rosaceus) develops via a rare, intermediate larval type-facultative feeding- and can develop through metamorphosis entirely based on egg provisioning in under one week.

RESULTS

Overall, the two transcriptomes are highly similar, containing largely orthologous contigs with similar functional annotation. However, we found distinct differences in gene expression patterns between the two species. Larvae from C. rosaceus, the facultative planktotroph, turned genes on at earlier stages and had less differentiation in gene expression between larval stages, whereas, C. subdepressus showed a higher degree of stage-specific gene expression.

CONCLUSION

This study is the first genetic analysis of a species with facultatively feeding larvae. Our results are consistent with known developmental differences between the larval types and raise the question of whether earlier onset of developmental genes is a key step in the evolution of a reduced larval period. By publishing a transcriptome for this rare, intermediate, larval type, this study adds developmental breadth to the current genetic resources, which will provide a valuable tool for future research on echinoderm development as well as studies on the evolution of development in general.

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.

Transcriptomic identification of starfish neuropeptide precursors yields new insights into neuropeptide evolution

Neuropeptides are evolutionarily ancient mediators of neuronal signalling in nervous systems. With recent advances in genomics/transcriptomics, an increasingly wide range of species has become accessible for molecular analysis. The deuterostomian invertebrates are of particular interest in this regard because they occupy an ‘intermediate' position in animal phylogeny, bridging the gap between the well-studied model protostomian invertebrates (e.g. Drosophila melanogaster, Caenorhabditis elegans) and the vertebrates. Here we have identified 40 neuropeptide precursors in the starfish Asterias rubens, a deuterostomian invertebrate from the phylum Echinodermata. Importantly, these include kisspeptin-type and melanin-concentrating hormone-type precursors, which are the first to be discovered in a non-chordate species. Starfish tachykinin-type, somatostatin-type, pigment-dispersing factor-type and corticotropin-releasing hormone-type precursors are the first to be discovered in the echinoderm/ambulacrarian clade of the animal kingdom. Other precursors identified include vasopressin/oxytocin-type, gonadotropin-releasing hormone-type, thyrotropin-releasing hormone-type, calcitonin-type, cholecystokinin/gastrin-type, orexin-type, luqin-type, pedal peptide/orcokinin-type, glycoprotein hormone-type, bursicon-type, relaxin-type and insulin-like growth factor-type precursors. This is the most comprehensive identification of neuropeptide precursor proteins in an echinoderm to date, yielding new insights into the evolution of neuropeptide signalling systems. Furthermore, these data provide a basis for experimental analysis of neuropeptide function in the unique context of the decentralized, pentaradial echinoderm bauplan.

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.

Contribution of transcription to animal early development

In mature gametes and during the oocyte-to-embryo transition, transcription is generally silenced and gene expression is post-transcriptionally regulated. However, we recently discovered that major transcription can occur immediately after fertilization, prior to pronuclear fusion, and in the first cell division of the oocyte-to-embryo transition in the nematode Ascaris suum. We postulate that the balance between transcriptional and post-transcriptional regulation during the oocyte-to-embryo transition may largely be determined by cell cycle length and thus the time available for the genome to be transcribed.

microRNAs regulate β-catenin of the Wnt signaling pathway in early sea urchin development

Development of complex multicellular organisms requires careful regulation at both transcriptional and post-transcriptional levels. Post-transcriptional gene regulation is in part mediated by a class of non-coding RNAs of 21-25 nucleotides in length known as microRNAs (miRNAs). β-catenin, regulated by the canonical Wnt signaling pathway, has a highly evolutionarily conserved function in patterning early metazoan embryos, in forming the Anterior-Posterior axis, and in establishing the endomesoderm. Using reporter constructs and site-directed mutagenesis, we identified at least three miRNA binding sites within the 3' untranslated region (3'UTR) of the sea urchin β-catenin. Further, blocking these three miRNA binding sites within the β-catenin 3'UTR to prevent regulation of endogenous β-catenin by miRNAs resulted in a minor increase in β-catenin protein accumulation that is sufficient to induce aberrant gut morphology and circumesophageal musculature. These phenotypes are likely the result of increased transcript levels of Wnt responsive endomesodermal regulatory genes. This study demonstrates the importance of miRNA regulation of β-catenin in early development.

Isolation and assessment of signaling proteins from synchronized cultures during egg activation and through the egg-to-embryo transition in sea urchins

Sea urchins are an excellent model system for investigating fertilization mechanisms and fundamental cell biological phenomenon such as release from quiescence, cell division, secretion, and basic signal transduction. The ease of gamete collection, fertilization, and culture is complemented by exquisite developmental synchronicity and the ability to carry out both large-scale biochemical studies and single-cell experiments. In particular, fertilization in echinoderms serves as a paradigm for a digital signaling event-a one-time only switch that launches the egg into the developmental pathway. Sperm-induced egg activation is dependent on the release of calcium from internal stores and subsequent effects on a myriad of cellular events such as exocytosis, cytoskeletal remodeling, and cell cycle reentry. Here we describe methods to investigate individual signaling proteins as well as global proteomic and phosphoproteomic changes involved in the initial steps of egg activation through the egg-to-embryo transition.

Analysis of the skin transcriptome in two oujiang color varieties of common carp

Background

Body color and coloration patterns are important phenotypic traits to maintain survival and reproduction activities. The Oujiang color varieties of common carp (Cyprinus carpio var. color), with a narrow distribution in Zhejiang Province of China and a history of aquaculture for over 1,200 years, consistently exhibit a variety of body color patterns. The molecular mechanism underlying diverse color patterns in these variants is unknown. To the practical end, it is essential to develop molecular markers that can distinguish different phenotypes and assist selective breeding.

Methodology/Principal Findings

In this exploratory study, we conducted Roche 454 transcriptome sequencing of two pooled skin tissue samples of Oujiang common carp, which correspond to distinct color patterns, red with big black spots (RB) and whole white (WW), and a total of 737,525 sequence reads were generated. The reads obtained in this study were co-assembled jointly with common carp Roche 454 sequencing reads downloaded from NCBI SRA database, resulting in 43,923 isotigs and 546,676 singletons. Over 31 thousand (31,445; 71.6%) isotigs were found with significant BLAST matches (E<1e-10) to the nr protein database, which corresponds to 12,597 annotated zebrafish genes. A total of 70,947 isotigs and singletons (transcripts) were annotated with Gene Ontology, and 60,221 transcripts were found with corresponding EC numbers. Out of 145 zebrafish pigmentation genes, orthologs for 117 were recovered in Oujiang color carp transcriptome, including 18 found only among singletons. Our transcriptome analysis revealed over 52,902 SNPs in Oujiang common carp, and identified 63 SNP markers that are putatively unique either for RB or WW.

Conclusions

The transcriptome of Oujiang color varieties of common carp obtained through this study, along with the pigmentation genes recovered and the color pattern-specific molecular markers developed, will facilitate future research on the molecular mechanism of color patterns and promote aquaculture of Oujiang color varieties of common carp through molecular marker assisted-selective breeding.

Tools for sea urchin genomic analysis

The Sea Urchin Genome Project Web site, SpBase ( http://SpBase.org ), in association with a suite of publicly available sequence comparison tools provides a platform from which to analyze genes and genomic sequences of sea urchin. This information system is specifically designed to support laboratory bench studies in cell and molecular biology. In particular these tools and datasets have supported the description of the gene regulatory networks of the purple sea urchin S. purpuratus. This chapter details methods to undertake in the first steps to find genes and noncoding regulatory sequences for further analysis.

ATP-binding cassette transporters protect sea urchin gametes and embryonic cells against the harmful effects of ultraviolet light

Embryos of marine organisms whose development occurs externally are particularly sensitive to ultraviolet (UV) light (bands A and B, respectively, UVA and UVB). ATP-binding cassette (ABC) transporters are the first line of cellular defense against chemical or physical stress. The present work investigated the involvement of ABC transporters on UVA or UVB effects on eggs, spermatozoa, and embryonic cells of the sea urchin Echinometra lucunter. Gametes or embryos were exposed to UVA (3.6-14.4 kJ m(-2)) or UVB (0.112-14.4 kJ m(-2)), and embryonic development was monitored by optical microscopy at different developmental stages in the presence or absence of the ABC-transporter blockers reversin205 (ABCB1 blocker) or MK571 (ABCC1 blocker). E. lucunter eggs, spermatozoa and embryos were resistant to UVA exposure. Resistance to the harmful effects of UVB was strongly associated to ABC transporter activity (embryos > eggs > spermatozoa). ABCB1 or ABCC1 blockage promoted the injurious effects of UVA on spermatozoa. ABCC1 transporter blockage increased UVB-dependent damage in eggs while ABCB1 transporter inhibition increased harmful effects of UVB in embryonic cells. ABC-transporter activity was not, however, affected by UVB exposure. In conclusion, the present study is the first report on the protective role of ABC transporters against harmful effects of UVA and UVB on sea urchin eggs and embryonic cells.

Proteomic responses of sea urchin embryos to stressful ultraviolet radiation

Solar ultraviolet radiation (UVR, 290-400 nm) penetrates into seawater and can harm shallow-dwelling and planktonic marine organisms. Studies dating back to the 1930s revealed that echinoids, especially sea urchin embryos, are powerful models for deciphering the effects of UVR on embryonic development and how embryos defend themselves against UV-induced damage. In addition to providing a large number of synchronously developing embryos amenable to cellular, biochemical, molecular, and single-cell analyses, the purple sea urchin, Strongylocentrotus purpuratus, also offers an annotated genome. Together, these aspects allow for the in-depth study of molecular and biochemical signatures of UVR stress. Here, we review the effects of UVR on embryonic development, focusing on the early-cleavage stages, and begin to integrate data regarding single-protein responses with comprehensive proteomic assessments. Proteomic studies reveal changes in levels of post-translational modifications to proteins that respond to UVR, and identify proteins that can then be interrogated as putative targets or components of stress-response pathways. These responsive proteins are distributed among systems upon which targeted studies can now begin to be mapped. Post-transcriptional and translational controls may provide early embryos with a rapid, fine-tuned response to stress during early stages, especially during pre-blastula stages that rely primarily on maternally derived defenses rather than on responses through zygotic gene transcription.

ATP-binding cassette (ABC) transporter expression and localization in sea urchin development

BACKGROUND

ATP-binding cassette (ABC) transporters are membrane proteins that regulate intracellular concentrations of myriad compounds and ions. There are >100 ABC transporter predictions in the Strongylocentrotus purpuratus genome, including 40 annotated ABCB, ABCC, and ABCG "multidrug efflux" transporters. Despite the importance of multidrug transporters for protection and signaling, their expression patterns have not been characterized in deuterostome embryos.

RESULTS

Sea urchin embryos expressed 20 ABCB, ABCC, and ABCG transporter genes in the first 58 hr of development, from unfertilized egg to early prism. We quantified transcripts of ABCB1a, ABCB4a, ABCC1, ABCC5a, ABCC9a, and ABCG2b, and found that ABCB1a mRNA was 10-100 times more abundant than other transporter mRNAs. In situ hybridization showed ABCB1a was expressed ubiquitously in embryos, while ABCC5a was restricted to secondary mesenchyme cells and their precursors. Fluorescent protein fusions showed localization of ABCB1a on apical cell surfaces, and ABCC5a on basolateral surfaces.

CONCLUSIONS

Embryos use many ABC transporters with predicted functions in cell signaling, lysosomal and mitochondrial homeostasis, potassium channel regulation, pigmentation, and xenobiotic efflux. Detailed characterization of ABCB1a and ABCC5a revealed that they have different temporal and spatial gene expression profiles and protein localization patterns that correlate to their predicted functions in protection and development, respectively.

Pantropic retroviruses as a transduction tool for sea urchin embryos

Sea urchins are an important model for experiments at the intersection of development and systems biology, and technical innovations that enhance the utility of this model are of great value. This study explores pantropic retroviruses as a transduction tool for sea urchin embryos, and demonstrates that pantropic retroviruses infect sea urchin embryos with high efficiency and genomically integrate at a copy number of one per cell. We successfully used a self-inactivation strategy to both insert a sea urchin-specific enhancer and disrupt the endogenous viral enhancer. The resulting self-inactivating viruses drive global and persistent gene expression, consistent with genomic integration during the first cell cycle. Together, these data provide substantial proof of principle for transduction technology in sea urchin embryos.

Molecular profiling of marine fauna: integration of omics with environmental assessment of the world's oceans

Many species that contribute to the commercial and ecological richness of our marine ecosystems are harbingers of environmental change. The ability of organisms to rapidly detect and respond to changes in the surrounding environment represents the foundation for application of molecular profiling technologies towards marine sentinel species in an attempt to identify signature profiles that may reside within the transcriptome, proteome, or metabolome and that are indicative of a particular environmental exposure event. The current review highlights recent examples of the biological information obtained for marine sentinel teleosts, mammals, and invertebrates. While in its infancy, such basal information can provide a systems biology framework in the detection and evaluation of environmental chemical contaminant effects on marine fauna. Repeated evaluation across different seasons and local marine environs will lead to discrimination between signature profiles representing normal variation within the complex milieu of environmental factors that trigger biological response in a given sentinel species and permit a greater understanding of normal versus anthropogenic-associated modulation of biological pathways, which prove detrimental to marine fauna. It is anticipated that incorporation of contaminant-specific molecular signatures into current risk assessment paradigms will lead to enhanced wildlife management strategies that minimize the impacts of our industrialized society on marine ecosystems.

The genomic regulatory control of skeletal morphogenesis in the sea urchin

A central challenge of developmental and evolutionary biology is to understand how anatomy is encoded in the genome. Elucidating the genetic mechanisms that control the development of specific anatomical features will require the analysis of model morphogenetic processes and an integration of biological information at genomic, cellular and tissue levels. The formation of the endoskeleton of the sea urchin embryo is a powerful experimental system for developing such an integrated view of the genomic regulatory control of morphogenesis. The dynamic cellular behaviors that underlie skeletogenesis are well understood and a complex transcriptional gene regulatory network (GRN) that underlies the specification of embryonic skeletogenic cells (primary mesenchyme cells, PMCs) has recently been elucidated. Here, we link the PMC specification GRN to genes that directly control skeletal morphogenesis. We identify new gene products that play a proximate role in skeletal morphogenesis and uncover transcriptional regulatory inputs into many of these genes. Our work extends the importance of the PMC GRN as a model developmental GRN and establishes a unique picture of the genomic regulatory control of a major morphogenetic process. Furthermore, because echinoderms exhibit diverse programs of skeletal development, the newly expanded sea urchin skeletogenic GRN will provide a foundation for comparative studies that explore the relationship between GRN evolution and morphological evolution.

Select microRNAs are essential for early development in the sea urchin

microRNAs (miRNAs) are small noncoding RNAs that mediate post-transcriptional gene regulation and have emerged as essential regulators of many developmental events. The transcriptional network during early embryogenesis of the purple sea urchin, Strongylocentrotus purpuratus, is well described and can serve as an excellent model to test functional contributions of miRNAs in embryogenesis. We examined the loss of function phenotypes of major components of the miRNA biogenesis pathway. Inhibition of de novo synthesis of Drosha and Dicer in the embryo led to consistent developmental defects, a failure to gastrulate, and embryonic lethality, including changes in the steady state levels of transcription factors and signaling molecules involved in germ layer specification. We annotated and profiled small RNA expression from the ovary and several early embryonic stages by deep sequencing followed by computational analysis. miRNAs as well as a large population of putative piRNAs (piwi-interacting RNAs) had dynamic accumulation profiles through early development. Defects in morphogenesis caused by loss of Drosha could be rescued with four miRNAs. Taken together our results indicate that post-transcriptional gene regulation directed by miRNAs is functionally important for early embryogenesis and is an integral part of the early embryonic gene regulatory network in S. purpuratus.

Expression pattern of polyketide synthase-2 during sea urchin development

Polyketide synthases (PKSs) are a large group of proteins responsible for the biosynthesis of polyketide compounds, which are mainly found in bacteria, fungi, and plants. Polyketides have a wide array of biological functions, including antibiotic, antifungal, predator defense, and light responses. In this study, we describe the developmental expression pattern of pks2, one of two pks found in the sea urchin genome. Throughout development, pks2 expression was restricted to skeletogenic cells and their precursors. Pks2 was first detected during the blastula stage. The transcript level peaked at hatched blastula, when all skeletogenic cell precursors expressed pks2. This was followed by a steady decline in expression in the skeletogenic cells on the aboral side of the embryo. By the prism stage, pks2 expression was limited to only 3-4 skeletogenic cells localized on the oral side.

Transcriptional increase and misexpression of 14-3-3 epsilon in sea urchin embryos exposed to UV-B

Members of the 14-3-3 protein family are involved in many important cellular events, including stress response, survival and apoptosis. Genes of the 14-3-3 family are conserved from plants to humans, and some members are responsive to UV radiation. Here, we report the isolation of the complete cDNA encoding the 14-3-3 epsilon isoform from Paracentrotus lividus sea urchin embryos, referred to as Pl14-3-3ε, and the phylogenetic relationship with other homologues described in different phyla. Pl14-3-3ε mRNA levels were measured by QPCR during development and found to increase from the mesenchyme blastula to the prism stage. In response to UV-B (312 nm) exposure, early stage embryos collected 2 h later showed a 2.3-fold (at 400 J/m(2)) and a 2.7-fold (at 800 J/m(2)) increase in Pl14-3-3ε transcript levels compared with controls. The spatial expression of Pl14-3-3ε mRNA, detected by whole mount in situ hybridization in both control and UV-B exposed embryos, harvested at late developmental stages, showed transcripts to be located in the archenteron of gastrula stage and widely distributed in all germ layers, respectively. The Pl14-3-3ε mRNA delocalization parallels the failure in archenteron elongation observed morphologically, as well as the lack of specific endoderm markers, investigated by indirect immuno-fluorescence on whole mount embryos. Results confirm the involvement of 14-3-3ε in the stress response elicited by UV-B and demonstrate, for the first time, its contribution at the transcriptional level in the sea urchin embryo.

Whole-genome positive selection and habitat-driven evolution in a shallow and a deep-sea urchin

Comparisons of genomic sequence between divergent species can provide insight into the action of natural selection across many distinct classes of proteins. Here, we examine the extent of positive selection as a function of tissue-specific and stage-specific gene expression in two closely-related sea urchins, the shallow-water Strongylocentrotus purpuratus and the deep-sea Allocentrotus fragilis, which have diverged greatly in their adult but not larval habitats. Genes that are expressed specifically in adult somatic tissue have significantly higher dN/dS ratios than the genome-wide average, whereas those in larvae are indistinguishable from the genome-wide average. Testis-specific genes have the highest dN/dS values, whereas ovary-specific have the lowest. Branch-site models involving the outgroup S. franciscanus indicate greater selection (ω_FG) along the A. fragilis branch than along the S. purpuratus branch. The A. fragilis branch also shows a higher proportion of genes under positive selection, including those involved in skeletal development, endocytosis, and sulfur metabolism. Both lineages are approximately equal in enrichment for positive selection of genes involved in immunity, development, and cell–cell communication. The branch-site models further suggest that adult-specific genes have experienced greater positive selection than those expressed in larvae and that ovary-specific genes are more conserved (i.e., experienced greater negative selection) than those expressed specifically in adult somatic tissues and testis. Our results chart the patterns of protein change that have occurred after habitat divergence in these two species and show that the developmental or functional context in which a gene acts can play an important role in how divergent species adapt to new environments.

Practical computational methods for regulatory genomics: a cisGRN-Lexicon and cisGRN-browser for gene regulatory networks

The CYRENE Project focuses on the study of cis-regulatory genomics and gene regulatory networks (GRN) and has three components: a cisGRN-Lexicon, a cisGRN-Browser, and the Virtual Sea Urchin software system. The project has been done in collaboration with Eric Davidson and is deeply inspired by his experimental work in genomic regulatory systems and gene regulatory networks. The current CYRENE cisGRN-Lexicon contains the regulatory architecture of 200 transcription factors encoding genes and 100 other regulatory genes in eight species: human, mouse, fruit fly, sea urchin, nematode, rat, chicken, and zebrafish, with higher priority on the first five species. The only regulatory genes included in the cisGRN-Lexicon (CYRENE genes) are those whose regulatory architecture is validated by what we call the Davidson Criterion: they contain functionally authenticated sites by site-specific mutagenesis, conducted in vivo, and followed by gene transfer and functional test. This is recognized as the most stringent experimental validation criterion to date for such a genomic regulatory architecture. The CYRENE cisGRN-Browser is a full genome browser tailored for cis-regulatory annotation and investigation. It began as a branch of the Celera Genome Browser (available as open source at http://sourceforge.net/projects/celeragb /) and has been transformed to a genome browser fully devoted to regulatory genomics. Its access paradigm for genomic data is zoom-to-the-DNA-base in real time. A more recent component of the CYRENE project is the Virtual Sea Urchin system (VSU), an interactive visualization tool that provides a four-dimensional (spatial and temporal) map of the gene regulatory networks of the sea urchin embryo.

The maternal-to-zygotic transition: a play in two acts

All animal embryos pass through a stage during which developmental control is handed from maternally provided gene products to those synthesized from the zygotic genome. This maternal-to-zygotic transition (MZT) has been extensively studied in model organisms, including echinoderms, nematodes, insects, fish,amphibians and mammals. In all cases, the MZT can be subdivided into two interrelated processes: first, a subset of maternal mRNAs and proteins is eliminated; second, zygotic transcription is initiated. The timing and scale of these two events differ across species, as do the cellular and morphogenetic processes that sculpt their embryos. In this article, we discuss conserved and distinct features within the two component processes of the MZT.

Building developmental gene regulatory networks

Animal development is an elaborate process programmed by genomic regulatory instructions. Regulatory genes encode transcription factors and signal molecules, and their expression is under the control of cis-regulatory modules that define the logic of transcriptional responses to the inputs of other regulatory genes. The functional linkages among regulatory genes constitute the gene regulatory networks (GRNs) that govern cell specification and patterning in development. Constructing such networks requires identification of the regulatory genes involved and characterization of their temporal and spatial expression patterns. Interactions (activation/repression) among transcription factors or signals can be investigated by large-scale perturbation analysis, in which the function of each gene is specifically blocked. Resultant expression changes are then integrated to identify direct linkages, and to reveal the structure of the GRN. Predicted GRN linkages can be tested and verified by cis-regulatory analysis. The explanatory power of the GRN was shown in the lineage specification of sea urchin endomesoderm. Acquiring such networks is essential for a systematic and mechanistic understanding of the developmental process.

Learning global models of transcriptional regulatory networks from data

Organisms must continually adapt to changing cellular and environmental factors (e.g., oxygen levels) by altering their gene expression patterns. At the same time, all organisms must have stable gene expression patterns that are robust to small fluctuations in environmental factors and genetic variation. Learning and characterizing the structure and dynamics of Regulatory Networks (RNs), on a whole-genome scale, is a key problem in systems biology. Here, we review the challenges associated with inferring RNs in a solely data-driven manner, concisely discuss the implications and contingencies of possible procedures that can be used, specifically focusing on one such procedure, the Inferelator. Importantly, the Inferelator explicitly models the temporal component of regulation, can learn the interactions between transcription factors and environmental factors, and attaches a statistically meaningful weight to every edge. The result of the Inferelator is a dynamical model of the RN that can be used to model the time-evolution of cell state.

Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx

BACKGROUND

New methods are needed for genomic-scale analysis of emerging model organisms that exemplify important biological questions but lack fully sequenced genomes. For example, there is an urgent need to understand the potential for corals to adapt to climate change, but few molecular resources are available for studying these processes in reef-building corals. To facilitate genomics studies in corals and other non-model systems, we describe methods for transcriptome sequencing using 454, as well as strategies for assembling a useful catalog of genes from the output. We have applied these methods to sequence the transcriptome of planulae larvae from the coral Acropora millepora.

RESULTS

More than 600,000 reads produced in a single 454 sequencing run were assembled into approximately 40,000 contigs with five-fold average sequencing coverage. Based on sequence similarity with known proteins, these analyses identified approximately 11,000 different genes expressed in a range of conditions including thermal stress and settlement induction. Assembled sequences were annotated with gene names, conserved domains, and Gene Ontology terms. Targeted searches using these annotations identified the majority of genes associated with essential metabolic pathways and conserved signaling pathways, as well as novel candidate genes for stress-related processes. Comparisons with the genome of the anemone Nematostella vectensis revealed approximately 8,500 pairs of orthologs and approximately 100 candidate coral-specific genes. More than 30,000 SNPs were detected in the coral sequences, and a subset of these validated by re-sequencing.

CONCLUSION

The methods described here for deep sequencing of the transcriptome should be widely applicable to generate catalogs of genes and genetic markers in emerging model organisms. Our data provide the most comprehensive sequence resource currently available for reef-building corals, and include an extensive collection of potential genetic markers for association and population connectivity studies. The characterization of the larval transcriptome for this widely-studied coral will enable research into the biological processes underlying stress responses in corals and evolutionary adaptation to global climate change.

Transcriptionally active gene fragments derived from potentially fast-evolving donor genes in the rice genome

The unprecedented complexity of the transcriptomic data obtained in recent years creates opportunities for new genomic studies aimed at interpolating regulatory code of gene expression and tracing genome evolution. We report here the identification and characterization of a set of 851 intergenic loci that represent transcribed gene fragments (TGFs) ectopically duplicated from 1030 non-transposable element (non-TE) donor genes in the rice genome. We analyzed the genomic context of the TGFs and donor genes. We show that the TGFs have adopted transcriptional orientation and pattern independent of the donor genes. We further show that TGFs have undergone relaxed purifying selection, consistent with their being pseudogenized. We found that the donor genes, which are biased toward certain molecular functions, exhibit an accelerated evolution rate comparing to the genome average. Our results demonstrated a large number of actively TGFs in the rice genome and shed light on the origin, mode of action and function of the TGFs.

Chordin is required for neural but not axial development in sea urchin embryos

The oral-aboral (OA) axis in the sea urchin is specified by the TGFbeta family members Nodal and BMP2/4. Nodal promotes oral specification, whereas BMP2/4, despite being expressed in the oral territory, is required for aboral specification. This study explores the role of Chordin (Chd) during sea urchin embryogenesis. Chd is a secreted BMP inhibitor that plays an important role in axial and neural specification and patterning in Drosophila and vertebrate embryos. In Lytechinus variegatus embryos, Chd and BMP2/4 are functionally antagonistic. Both are expressed in overlapping domains in the oral territory prior to and during gastrulation. Perturbation shows that, surprisingly, Chd is not involved in OA axis specification. Instead, Chd is required both for normal patterning of the ciliary band at the OA boundary and for development of synaptotagmin B-positive (synB) neurons in a manner that is reciprocal with BMP2/4. Chd expression and synB-positive neural development are both downstream from p38 MAPK and Nodal, but not Goosecoid. These data are summarized in a model for synB neural development.

Expression of multiple Src family kinases in sea urchin eggs and their function in Ca2+ release at fertilization

Egg activation at fertilization in deuterostomes requires a rise in intracellular Ca(2+), which is released from the egg's endoplasmic reticulum. In sea urchins, a Src Family Kinase (SpSFK1) is necessary for the PLCgamma-mediated signaling event that initiates this Ca(2+) release (Giusti, A.F., O'Neill, F.J., Yamasu, K., Foltz, K.R. and Jaffe, L.A., 2003. Function of a sea urchin egg Src family kinase in initiating Ca2+ release at fertilization. Dev. Biol. 256, 367-378.). Annotation of the Strongylocentrotus purpuratus genome sequence led to the identification of additional, predicted SFKs (Bradham, C.A., Foltz, D.R., Beane, W.S., Amone, M.I., Rizzo, F., Coffman, J.A., Mushegian, A., Goel, M., Morales, J., Geneviere, A.M., Lapraz, F., Robertson, A.J., Kelkar, H., Loza-Coll, M., Townley, I.K., Raisch, M., Roux, M.M., Lepage, T., Gache, C., McClay, D.R., Manning, G., 2006. The sea urchin kinome: a first look. Dev. Biol. 300, 180-193.; Roux, M.M., Townley, I.K., Raisch, M., Reade, A., Bradham, C., Humphreys, G., Gunaratne, H.J., Killian, C.E., Moy, G., Su, Y.H., Ettensohn, C.A., Wilt, F., Vacquier, V.D., Burke, R.D., Wessel, G. and Foltz, K.R., 2006. A functional genomic and proteomic perspective of sea urchin calcium signaling and egg activation. Dev. Biol. 300, 416-433.). Here, we describe the cloning and characterization of these 4 additional SFKs and test their function during the initial Ca(2+) release at fertilization using the dominant-interfering microinjection method coupled with Ca(2+) recording. While two of the new SFKs (SpFrk and SpSFK3) are necessary for Ca(2+) release, SpSFK5 appears dispensable for early egg to embryo transition events. Interestingly, SpSFK7 may be involved in preventing precocious release of Ca(2+). Binding studies indicate that only SpSFK1 is capable of direct interaction with PLCgamma. Immunolocalization studies suggest that one or more SpSFK and PLCgamma are localized to the egg cortex and at the site of sperm-egg interaction. Collectively, these data indicate that more than one SFK is involved in the Ca(2+) release pathway at fertilization.

SpBase: the sea urchin genome database and web site

SpBase is a system of databases focused on the genomic information from sea urchins and related echinoderms. It is exposed to the public through a web site served with open source software (http://spbase.org/). The enterprise was undertaken to provide an easily used collection of information to directly support experimental work on these useful research models in cell and developmental biology. The information served from the databases emerges from the draft genomic sequence of the purple sea urchin, Strongylocentrotus purpuratus and includes sequence data and genomic resource descriptions for other members of the echinoderm clade which in total span 540 million years of evolutionary time. This version of the system contains two assemblies of the purple sea urchin genome, associated expressed sequences, gene annotations and accessory resources. Search mechanisms for the sequences and the gene annotations are provided. Because the system is maintained along with the Sea Urchin Genome resource, a database of sequenced clones is also provided.

cDNA sequences for transcription factors and signaling proteins of the hemichordate Saccoglossus kowalevskii: efficacy of the expressed sequence tag (EST) approach for evolutionary and developmental studies of a new organism

We describe a collection of expressed sequence tags (ESTs) for Saccoglossus kowalevskii, a direct-developing hemichordate valuable for evolutionary comparisons with chordates. The 202,175 ESTs represent 163,633 arrayed clones carrying cDNAs prepared from embryonic libraries, and they assemble into 13,677 continuous sequences (contigs), leaving 10,896 singletons (excluding mitochondrial sequences). Of the contigs, 53% had significant matches when BLAST was used to query the NCBI databases (< or = 10(-10)), as did 51% of the singletons. Contigs most frequently matched sequences from amphioxus (29%), chordates (67%), and deuterostomes (87%). From the clone array, we isolated 400 full-length sequences for transcription factors and signaling proteins of use for evolutionary and developmental studies. The set includes sequences for fox, pax, tbx, hox, and other homeobox-containing factors, and for ligands and receptors of the TGFbeta, Wnt, Hh, Delta/Notch, and RTK pathways. At least 80% of key sequences have been obtained, when judged against gene lists of model organisms. The median length of these cDNAs is 2.3 kb, including 1.05 kb of 3' untranslated region (UTR). Only 30% are entirely matched by single contigs assembled from ESTs. We conclude that an EST collection based on 150,000 clones is a rich source of sequences for molecular developmental work, and that the EST approach is an efficient way to initiate comparative studies of a new organism.

Remodeling the model organism: matrix metalloproteinase functions in invertebrates

The matrix metalloproteinase (MMP) family of extracellular proteases is conserved throughout the animal kingdom. Studies of invertebrate MMPs have demonstrated they are involved in tissue remodeling. In Drosophila, MMPs are required for histolysis, tracheal growth, tissue invasion, axon guidance, and dendritic remodeling. Recent work demonstrates that MMPs also participate in Drosophila tumor invasion. In Caenorhabditis elegans an MMP is involved in anchor cell invasion; a Hydra MMP is important for regeneration and maintaining cell identity; and a sea urchin MMP degrades matrix to allow hatching. In worms and in flies, MMPs are regulated by the JNK pathway.

Noncoding but nonexpendable: transcriptional regulation by large noncoding RNA in eukaryotes

Genome sequencing and annotation has advanced our understanding of genome organization and gene structure but initially only allowed predictions of how many genes might be present. Mechanisms such as alternative splicing reveal that these predictions only scratch the surface of the true nature of the transcriptome. Several thousand expressed partial gene fragments have been cloned but were considered transcriptional noise or cloning artifacts. We now know that genomes are indeed expressed at much higher levels than was previously predicted, and much of the additional transcription maps to intergenic regions, intron sequences, and untranslated regions of mRNAs. These transcripts are expressed from either the sense or the antisense strand and can be confirmed by conventional techniques. In addition to the already established roles for small RNAs in gene regulation, large noncoding RNAs (ncRNAs) are also emerging as potent regulators of gene expression. In this review, we summarize several illustrative examples of gene regulatory mechanisms that involve large ncRNAs. We describe several distinct regulatory mechanisms that involve large ncRNAs, such as transcriptional interference and promoter inactivation, as well as indirect effects on transcription regulatory proteins and in genomic imprinting. These diverse functions for large ncRNAs are likely to be only the first of many novel regulatory mechanisms emerging from this growing field.

The grateful dead: damage-associated molecular pattern molecules and reduction/oxidation regulate immunity

The response to pathogens and damage in plants and animals involves a series of carefully orchestrated, highly evolved, molecular mechanisms resulting in pathogen resistance and wound healing. In metazoans, damage- or pathogen-associated molecular pattern molecules (DAMPs, PAMPs) execute precise intracellular tasks and are also able to exert disparate functions when released into the extracellular space. The emergent consequence for both inflammation and wound healing of the abnormal extracellular persistence of these factors may underlie many clinical disorders. DAMPs/PAMPs are recognized by hereditable receptors including the Toll-like receptors, the NOD1-like receptors and retinoic-acid-inducible gene I-like receptors, as well as the receptor for advanced glycation end products. These host molecules 'sense' not only pathogens but also misfolded/glycated proteins or exposed hydrophobic portions of molecules, activating intracellular cascades that lead to an inflammatory response. Equally important are means to not only respond to these molecules but also to eradicate them. We have speculated that their destruction through oxidative mechanisms normally exerted by myeloid cells, such as neutrophils and eosinophils, or their persistence in the setting of pathologic extracellular reducing environments, maintained by exuberant necrotic cell death and/or oxidoreductases, represent important molecular means enabling chronic inflammatory states.

Systems perspectives on mRNA processing

The application of genomic technologies to the study of mRNA processing is increasingly conducted in metazoan organisms in order to understand the complex events that occur during and after transcription. Large-scale systems analyses of mRNA-protein interactions and mRNA dynamics have revealed specificity in mRNA transcription, splicing, transport, translation, and turnover, and have begun to make connections between the different layers of mRNA processing. Here, we review global studies of post-transcriptional processes and discuss the challenges facing our understanding of mRNA regulation in metazoan organisms. In parallel, we examine genome-scale investigations that have expanded our knowledge of RNA-binding proteins and the networks of mRNAs that they regulate.

A global view of gene expression in lithium and zinc treated sea urchin embryos: new components of gene regulatory networks

BACKGROUND

The genome of the sea urchin Strongylocentrotus purpuratus has recently been sequenced because it is a major model system for the study of gene regulatory networks. Embryonic expression patterns for most genes are unknown, however.

RESULTS

Using large-scale screens on arrays carrying 50% to 70% of all genes, we identified novel territory-specific markers. Our strategy was based on computational selection of genes that are differentially expressed in lithium-treated embryos, which form excess endomesoderm, and in zinc-treated embryos, in which endomesoderm specification is blocked. Whole-mount in situ hybridization (WISH) analysis of 700 genes indicates that the apical organ region is eliminated in lithium-treated embryos. Conversely, apical and specifically neural markers are expressed more broadly in zinc-treated embryos, whereas endomesoderm signaling is severely reduced. Strikingly, the number of serotonergic neurons is amplified by at least tenfold in zinc-treated embryos. WISH analysis further indicates that there is crosstalk between the Wnt (wingless int), Notch, and fibroblast growth factor signaling pathways in secondary mesoderm cell specification and differentiation, similar to signaling cascades that function during development of presomitic mesoderm in mouse embryogenesis. We provide differential expression data for more than 4,000 genes and WISH patterns of more than 250 genes, and more than 2,400 annotated WISH images.

CONCLUSION

Our work provides tissue-specific expression patterns for a large fraction of the sea urchin genes that have not yet been included in existing regulatory networks and await functional integration. Furthermore, we noted neuron-inducing activity of zinc on embryonic development; this is the first observation of such activity in any organism.

Dissection of the Bradyrhizobium japonicum NifA+sigma54 regulon, and identification of a ferredoxin gene (fdxN) for symbiotic nitrogen fixation

Hierarchically organized regulatory proteins form a complex network for expression control of symbiotic and accessory genes in the nitrogen-fixing soybean symbiont Bradyrhizobium japonicum. A genome-wide survey of regulatory interactions was made possible with the design of a custom-made gene chip. Here, we report the first use of the microarray in a comprehensive and complete characterization of the B. japonicum NifA+sigma(54) regulon which forms an important node in the entire network. Comparative transcript profiles of anaerobically grown wild-type, nifA, and rpoN (1/2) mutant cells were complemented with a position-specific frequency matrix-based search for NifA- and sigma(54)-binding sites plus a simple operon definition. One of the newly identified NifA+sigma(54)-dependent genes, fdxN, encodes a ferredoxin required for efficient symbiotic nitrogen fixation, which makes it a candidate for being a direct electron donor to nitrogenase. The fdxN gene has an unconventional, albeit functional sigma(54 )promoter with the dinucleotide GA instead of the consensus GC motif at position -12. A GC-containing mutant promoter and the atypical GA-containing promoter of the wild type were disparately activated. Expression analyses were also carried out with two other NifA+sigma(54) targets (ectC; ahpC). Incidentally, the tiling-like design of the microarray has helped to arrive at completely revised annotations of the ectC- and ahpC-upstream DNA regions, which are now compatible with promoter locations. Taken together, the approaches used here led to a substantial expansion of the NifA+sigma(54 )regulon size, culminating in a total of 65 genes for nitrogen fixation and diverse other processes.

Genome-wide transcription and the implications for genomic organization

Recent evidence of genome-wide transcription in several species indicates that the amount of transcription that occurs cannot be entirely accounted for by current sets of genome-wide annotations. Evidence indicates that most of both strands of the human genome might be transcribed, implying extensive overlap of transcriptional units and regulatory elements. These observations suggest that genomic architecture is not colinear, but is instead interleaved and modular, and that the same genomic sequences are multifunctional: that is, used for multiple independently regulated transcripts and as regulatory regions. What are the implications and consequences of such an interleaved genomic architecture in terms of increased information content, transcriptional complexity, evolution and disease states?

A basal deuterostome genome viewed as a natural experiment

With the determination of its genome sequence the utility of the sea urchin model system increases. The phylogenetic position of the sea urchin among the deuterostomes allows for informative comparisons to vertebrate research models. A combined whole genome shotgun and bacterial artificial chromosome based strategy yielded a high quality draft genome sequence of 814 Mb. The predicted gene set estimated to include 23,300 genes was annotated and compared to those of other metazoan animals. Gene family expansions in the innate immune system are large and offer a first glimpse of how the long-lived sea urchin defends itself. The gene sets of the sea urchin place it firmly among the deuterostomes and indicate that various gene family-specific expansions and contractions characterize the evolution of animal genomes rather than the invention of new genes.