A combined morphological and molecular phylogeny for sea urchins (Echinoidea: Echinodermata)

Taxonomic Diversity, Predicted Metabolic Pathway, and Interaction Pattern of Bacterial Community in Sea Urchin Anthocidaris crassispina

Bacterial assemblages associated with sea urchin are critical to their physiology and ecology within marine ecosystems. In this study, we characterized the bacterial communities in wild sea urchin Anthocidaris crassispina captured in Daya Bay, South China Sea. A total of 363 amplicon sequence variants belonging to nine phyla and 141 genera were classified from intestine, body surface, and surrounding seawater samples. Proteobacteria, Firmicutes, and Bacteroidetes were the dominant bacteria phyla found in this study. A network analysis of bacterial interspecies interactions revealed varying complexity, stability, connectivity, and relationship patterns across the samples, with the most intricate network observed in the surrounding seawater. Metagenomic predictions highlighted the distinct bacterial metabolic pathways, with significant differences between intestine and seawater samples. Notably, pathways associated with polysaccharide degradation, including chitin derivatives, starch, and CoM biosynthesis, were markedly abundant, underscoring the gut microbiota's key role in digesting algae. In addition, other metabolic pathways in intestine samples were linked to immune response regulation of sea urchins. Overall, this study provides a comprehensive overview of the bacterial community structure and potential functional roles in A. crassispina.

Characterization of thyrotropin-releasing hormone producing neurons in sea urchin, from larva to juvenile

Most sea urchin species are indirect developers, going through a larval stage called pluteus. The pluteus possesses its own nervous system, consisting mainly of the apical organ neurons (controlling metamorphosis and settlement) and ciliary band neurons (controlling swimming behavior and food collection). Additional neurons are located in various areas of the gut. In recent years, the molecular complexity of this apparently "simple" nervous system has become apparent, with at least 12 neuronal populations identified through scRNA-sequencing in the species Strongylocentrotus purpuratus. Among these, there is a cluster of neurosecretory cells that produce a thyrotropin-releasing hormone-type neuropeptide (TRHergic) and that are also photosensory (expressing a Go-Opsin). However, much less is known about the organization of the nervous system in other sea urchin species. The aim of this work was to thoroughly characterize the localization of the TRHergic cells from early pluteus to juvenile stages in the Mediterranean sea urchin species Paracentrotus lividus combining immunostaining and whole mount in situ hybridization. We also compared the localization of TRHergic cells in early plutei of two other sea urchin species, Arbacia lixula and Heliocidaris tuberculata. This work provides new information on the anatomy and development of the nervous system in sea urchins. Moreover, by comparing the molecular signature of the TRHergic cells in P. lividus and S. purpuratus, we have obtained new insights how TRH-type neuropeptide signaling evolved in relatively closely related species.

Phylogeny, ancestral ranges and reclassification of sand dollars

Classification of the Class Echinoidea is under significant revision in light of emerging molecular phylogenetic evidence. In particular, the sister-group relationships within the superorder Luminacea (Echinoidea: Irregularia) have been considerably updated. However, the placement of many families remains largely unresolved due to a series of incongruent evidence obtained from morphological, paleontological, and genetic data for the majority of extant representatives. In this study, we investigated the phylogenetic relationships of 25 taxa, belonging to eleven luminacean families. We proposed three new superfamilies: Astriclypeoidea, Mellitoidea, and Taiwanasteroidea (including Dendrasteridae, Taiwanasteridae, Scutellidae, and Echinarachniidae), instead of the currently recognized superfamily Scutelloidea Gray, 1825. In light of the new data obtained from ten additional species, the historical biogeography reconstructed shows that the tropical western Pacific and eastern Indian Oceans are the cradle for early sand dollar diversification. Hothouse conditions during the late Cretaceous and early Paleogene were coupled with diversification events of major clades of sand dollars. We also demonstrate that Taiwan fauna can play a key role in terms of understanding the major Cenozoic migration and dispersal events in the evolutionary history of Luminacea.

The evolution of neurosensation provides opportunities and constraints for phenotypic plasticity

Phenotypic plasticity is widely regarded as important for enabling species resilience to environmental change and for species evolution. However, insight into the complex mechanisms by which phenotypic plasticity evolves in nature is limited by our ability to reconstruct evolutionary histories of plasticity. By using part of the molecular mechanism, we were able to trace the evolution of pre-feeding phenotypic plasticity across the class Echinoidea and identify the origin of plasticity at the base of the regular urchins. The neurosensory foundation for plasticity was ancestral within the echinoids. However, coincident development of the plastic trait and the neurosensory system was not achieved until the regular urchins, likely due to pleiotropic effects and linkages between the two colocalized systems. Plasticity continues to evolve within the urchins with numerous instances of losses associated with loss of sensory abilities and neurons, consistent with a cost of maintaining these capabilities. Thus, evidence was found for the neurosensory system providing opportunities and constraints to the evolution of phenotypic plasticity.

Mitogenome-wise codon usage pattern from comparative analysis of the first mitogenome of Blepharipa sp. (Muga uzifly) with other Oestroid flies

Uziflies (Family: Tachinidae) are dipteran endoparasites of sericigenous insects which cause major economic loss in the silk industry globally. Here, we are presenting the first full mitogenome of Blepharipa sp. (Acc: KY644698, 15,080 bp, A + T = 78.41%), a dipteran parasitoid of Muga silkworm (Antheraea assamensis) found in the Indian states of Assam and Meghalaya. This study has confirmed that Blepharipa sp. mitogenome gene content and arrangement is similar to other Tachinidae and Sarcophagidae flies of Oestroidea superfamily, typical of ancestral Diptera. Although, Calliphoridae and Oestridae flies have undergone tRNA translocation and insertion, forming unique intergenic spacers (IGS) and overlapping regions (OL) and a few of them (IGS, OL) have been conserved across Oestroidea flies. The Tachinidae mitogenomes exhibit more AT content and AT biased codons in their protein-coding genes (PCGs) than the Oestroidea counterpart. About 92.07% of all (3722) codons in PCGs of this new species have A/T in their 3rd codon position. The high proportion of AT and repeats in the control region (CR) affects sequence coverage, resulting in a short CR (Blepharipa sp.: 168 bp) and a smaller tachinid mitogenome. Our research unveils those genes with a high AT content had a reduced effective number of codons, leading to high codon usage bias. The neutrality test shows that natural selection has a stronger influence on codon usage bias than directed mutational pressure. This study also reveals that longer PCGs (e.g., nad5, cox1) have a higher codon usage bias than shorter PCGs (e.g., atp8, nad4l). The divergence rates increase nonlinearly as AT content at the 3rd codon position increases and higher rate of synonymous divergence than nonsynonymous divergence causes strong purifying selection. The phylogenetic analysis explains that Blepharipa sp. is well suited in the family of insectivorous tachinid maggots. It's possible that biased codon usage in the Tachinidae family reduces the effective number of codons, and purifying selection retains the core functions in their mitogenome, which could help with efficient metabolism in their endo-parasitic life style and survival strategy.

Phylogenomic analyses of echinoid diversification prompt a re-evaluation of their fossil record

Echinoids are key components of modern marine ecosystems. Despite a remarkable fossil record, the emergence of their crown group is documented by few specimens of unclear affinities, rendering their early history uncertain. The origin of sand dollars, one of its most distinctive clades, is also unclear due to an unstable phylogenetic context. We employ 18 novel genomes and transcriptomes to build a phylogenomic dataset with a near-complete sampling of major lineages. With it, we revise the phylogeny and divergence times of echinoids, and place their history within the broader context of echinoderm evolution. We also introduce the concept of a chronospace - a multidimensional representation of node ages - and use it to explore methodological decisions involved in time calibrating phylogenies. We find the choice of clock model to have the strongest impact on divergence times, while the use of site-heterogeneous models and alternative node prior distributions show minimal effects. The choice of loci has an intermediate impact, affecting mostly deep Paleozoic nodes, for which clock-like genes recover dates more congruent with fossil evidence. Our results reveal that crown group echinoids originated in the Permian and diversified rapidly in the Triassic, despite the relative lack of fossil evidence for this early diversification. We also clarify the relationships between sand dollars and their close relatives and confidently date their origins to the Cretaceous, implying ghost ranges spanning approximately 50 million years, a remarkable discrepancy with their rich fossil record.

Toxicological Impact of Rare Earth Elements (REEs) on the Reproduction and Development of Aquatic Organisms Using Sea Urchins as Biological Models

The growing presence of lanthanides in the environment has drawn the attention of the scientific community on their safety and toxicity. The sources of lanthanides in the environment include diagnostic medicine, electronic devices, permanent magnets, etc. Their exponential use and the poor management of waste disposal raise serious concerns about the quality and safety of the ecosystems at a global level. This review focused on the impact of lanthanides in marine organisms on reproductive fitness, fertilization and embryonic development, using the sea urchin as a biological model system. Scientific evidence shows that exposure to lanthanides triggers a wide variety of toxic insults, including reproductive performance, fertilization, redox metabolism, embryogenesis, and regulation of embryonic gene expression. This was thoroughly demonstrated for gadolinium, the most widely used lanthanide in diagnostic medicine, whose uptake in sea urchin embryos occurs in a time- and concentration-dependent manner, correlates with decreased calcium absorption and primarily affects skeletal growth, with incorrect regulation of the skeletal gene regulatory network. The results collected on sea urchin embryos demonstrate a variable sensitivity of the early life stages of different species, highlighting the importance of testing the effects of pollution in different species. The accumulation of lanthanides and their emerging negative effects make risk assessment and consequent legislative intervention on their disposal mandatory.

Phylogeography, colouration, and cryptic speciation across the Indo-Pacific in the sea urchin genus Echinothrix

The sea urchins Echinothrix calamaris and Echinothrix diadema have sympatric distributions throughout the Indo-Pacific. Diverse colour variation is reported in both species. To reconstruct the phylogeny of the genus and assess gene flow across the Indo-Pacific we sequenced mitochondrial 16S rDNA, ATPase-6, and ATPase-8, and nuclear 28S rDNA and the Calpain-7 intron. Our analyses revealed that E. diadema formed a single trans-Indo-Pacific clade, but E. calamaris contained three discrete clades. One clade was endemic to the Red Sea and the Gulf of Oman. A second clade occurred from Malaysia in the West to Moorea in the East. A third clade of E. calamaris was distributed across the entire Indo-Pacific biogeographic region. A fossil calibrated phylogeny revealed that the ancestor of E. diadema diverged from the ancestor of E. calamaris ~ 16.8 million years ago (Ma), and that the ancestor of the trans-Indo-Pacific clade and Red Sea and Gulf of Oman clade split from the western and central Pacific clade ~ 9.8 Ma. Time since divergence and genetic distances suggested species level differentiation among clades of E. calamaris. Colour variation was extensive in E. calamaris, but not clade or locality specific. There was little colour polymorphism in E. diadema.

A Total-Evidence Dated Phylogeny of Echinoidea Combining Phylogenomic and Paleontological Data

Phylogenomic and paleontological data constitute complementary resources for unraveling the phylogenetic relationships and divergence times of lineages, yet few studies have attempted to fully integrate them. Several unique properties of echinoids (sea urchins) make them especially useful for such synthesizing approaches, including a remarkable fossil record that can be incorporated into explicit phylogenetic hypotheses. We revisit the phylogeny of crown group Echinoidea using a total-evidence dating approach that combines the largest phylogenomic data set for the clade, a large-scale morphological matrix with a dense fossil sampling, and a novel compendium of tip and node age constraints. To this end, we develop a novel method for subsampling phylogenomic data sets that selects loci with high phylogenetic signal, low systematic biases, and enhanced clock-like behavior. Our results demonstrate that combining different data sources increases topological accuracy and helps resolve conflicts between molecular and morphological data. Notably, we present a new hypothesis for the origin of sand dollars, and restructure the relationships between stem and crown echinoids in a way that implies a long stretch of undiscovered evolutionary history of the crown group in the late Paleozoic. Our efforts help bridge the gap between phylogenomics and phylogenetic paleontology, providing a model example of the benefits of combining the two. [Echinoidea; fossils; paleontology; phylogenomics; time calibration; total evidence.].

Characterization of the Gut Microbiota of the Antarctic Heart Urchin (Spatangoida) Abatus agassizii

Abatus agassizii is an irregular sea urchin species that inhabits shallow waters of South Georgia and South Shetlands Islands. As a deposit-feeder, A. agassizii nutrition relies on the ingestion of the surrounding sediment in which it lives barely burrowed. Despite the low complexity of its feeding habit, it harbors a long and twice-looped digestive tract suggesting that it may host a complex bacterial community. Here, we characterized the gut microbiota of specimens from two A. agassizii populations at the south of the King George Island in the West Antarctic Peninsula. Using a metabarcoding approach targeting the 16S rRNA gene, we characterized the Abatus microbiota composition and putative functional capacity, evaluating its differentiation among the gut content and the gut tissue in comparison with the external sediment. Additionally, we aimed to define a core gut microbiota between A. agassizii populations to identify potential keystone bacterial taxa. Our results show that the diversity and the composition of the microbiota, at both genetic and predicted functional levels, were mostly driven by the sample type, and to a lesser extent by the population location. Specific bacterial taxa, belonging mostly to Planctomycetacia and Spirochaetia, were differently enriched in the gut content and the gut tissue, respectively. Predictive functional profiles revealed higher abundance of specific pathways, as the sulfur cycle in the gut content and the amino acid metabolism, in the gut tissue. Further, the definition of a core microbiota allowed to obtain evidence of specific localization of bacterial taxa and the identification of potential keystone taxa assigned to the Desulfobacula and Spirochaeta genera as potentially host selected. The ecological relevance of these keystone taxa in the host metabolism is discussed.

Characterization and phylogenetic analysis of the complete mitochondrial genome of Stephnometra indica (Pelmatozoa: Crinoidea)

The crinoid Stephnometra indica inhabits coral reefs and surrounding waters in the western Pacific Ocean. In this study, the complete mitochondrial genome of S. indica (15,878 bp) was confirmed to contain 13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes. The gene order of S. indica was identical to that of previously published Crinoidea. The nucleotide composition of the S. indica mitogenome was biased toward A + T nucleotide 73.68%. There were four unassigned regions (UASs) found in S. indica and all four UAS regions were AT-biased. The phylogenetic relationship illustrated that S. indica is more closely related to Antedon mediterranea within the family Crinoidea.

A phylogenomic resolution of the sea urchin tree of life

BACKGROUND

Echinoidea is a clade of marine animals including sea urchins, heart urchins, sand dollars and sea biscuits. Found in benthic habitats across all latitudes, echinoids are key components of marine communities such as coral reefs and kelp forests. A little over 1000 species inhabit the oceans today, a diversity that traces its roots back at least to the Permian. Although much effort has been devoted to elucidating the echinoid tree of life using a variety of morphological data, molecular attempts have relied on only a handful of genes. Both of these approaches have had limited success at resolving the deepest nodes of the tree, and their disagreement over the positions of a number of clades remains unresolved.

RESULTS

We performed de novo sequencing and assembly of 17 transcriptomes to complement available genomic resources of sea urchins and produce the first phylogenomic analysis of the clade. Multiple methods of probabilistic inference recovered identical topologies, with virtually all nodes showing maximum support. In contrast, the coalescent-based method ASTRAL-II resolved one node differently, a result apparently driven by gene tree error induced by evolutionary rate heterogeneity. Regardless of the method employed, our phylogenetic structure deviates from the currently accepted classification of echinoids, with neither Acroechinoidea (all euechinoids except echinothurioids), nor Clypeasteroida (sand dollars and sea biscuits) being monophyletic as currently defined. We show that phylogenetic signal for novel resolutions of these lineages is strong and distributed throughout the genome, and fail to recover systematic biases as drivers of our results.

CONCLUSIONS

Our investigation substantially augments the molecular resources available for sea urchins, providing the first transcriptomes for many of its main lineages. Using this expanded genomic dataset, we resolve the position of several clades in agreement with early molecular analyses but in disagreement with morphological data. Our efforts settle multiple phylogenetic uncertainties, including the position of the enigmatic deep-sea echinothurioids and the identity of the sister clade to sand dollars. We offer a detailed assessment of evolutionary scenarios that could reconcile our findings with morphological evidence, opening up new lines of research into the development and evolutionary history of this ancient clade.

Mind the gap! The mitochondrial control region and its power as a phylogenetic marker in echinoids

BACKGROUND

In Metazoa, mitochondrial markers are the most commonly used targets for inferring species-level molecular phylogenies due to their extremely low rate of recombination, maternal inheritance, ease of use and fast substitution rate in comparison to nuclear DNA. The mitochondrial control region (CR) is the main non-coding area of the mitochondrial genome and contains the mitochondrial origin of replication and transcription. While sequences of the cytochrome oxidase subunit 1 (COI) and 16S rRNA genes are the prime mitochondrial markers in phylogenetic studies, the highly variable CR is typically ignored and not targeted in such analyses. However, the higher substitution rate of the CR can be harnessed to infer the phylogeny of closely related species, and the use of a non-coding region alleviates biases resulting from both directional and purifying selection. Additionally, complete mitochondrial genome assemblies utilizing next generation sequencing (NGS) data often show exceptionally low coverage at specific regions, including the CR. This can only be resolved by targeted sequencing of this region.

RESULTS

Here we provide novel sequence data for the echinoid mitochondrial control region in over 40 species across the echinoid phylogenetic tree. We demonstrate the advantages of directly targeting the CR and adjacent tRNAs to facilitate complementing low coverage NGS data from complete mitochondrial genome assemblies. Finally, we test the performance of this region as a phylogenetic marker both in the lab and in phylogenetic analyses, and demonstrate its superior performance over the other available mitochondrial markers in echinoids.

CONCLUSIONS

Our target region of the mitochondrial CR (1) facilitates the first thorough investigation of this region across a wide range of echinoid taxa, (2) provides a tool for complementing missing data in NGS experiments, and (3) identifies the CR as a powerful, novel marker for phylogenetic inference in echinoids due to its high variability, lack of selection, and high compatibility across the entire class, outperforming conventional mitochondrial markers.

Paleogenomics of echinoids reveals an ancient origin for the double-negative specification of micromeres in sea urchins

Establishing a timeline for the evolution of novelties is a common, unifying goal at the intersection of evolutionary and developmental biology. Analyses of gene regulatory networks (GRNs) provide the ability to understand the underlying genetic and developmental mechanisms responsible for the origin of morphological structures both in the development of an individual and across entire evolutionary lineages. Accurately dating GRN novelties, thereby establishing a timeline for GRN evolution, is necessary to answer questions about the rate at which GRNs and their subcircuits evolve, and to tie their evolution to paleoenvironmental and paleoecological changes. Paleogenomics unites the fossil record and all aspects of deep time, with modern genomics and developmental biology to understand the evolution of genomes in evolutionary time. Recent work on the regulatory genomic basis of development in cidaroid echinoids, sand dollars, heart urchins, and other nonmodel echinoderms provides an ideal dataset with which to explore GRN evolution in a comparative framework. Using divergence time estimation and ancestral state reconstructions, we have determined the age of the double-negative gate (DNG), the subcircuit which specifies micromeres and skeletogenic cells in Strongylocentrotus purpuratus We have determined that the DNG has likely been used for euechinoid echinoid micromere specification since at least the Late Triassic. The innovation of the DNG thus predates the burst of post-Paleozoic echinoid morphological diversification that began in the Early Jurassic. Paleogenomics has wide applicability for the integration of deep time and molecular developmental data, and has wide utility in rigorously establishing timelines for GRN evolution.

Echinoderm phylogeny: how congruent are morphological and molecular estimates?

Single data sets, whether derived from morphological or molecular evidence, provide one-off estimates of the correct phylogeny. Their reliability can only be gauged by statistical approaches such as bootstrapping or clade decay, but these test only whether there are sufficient characters in the data matrix to justify the groupings identified. They do not test whether the characters themselves are reliable. Consequently, confidence in the correctness of phylogenetic interpretations comes primarily from discovering the same (or statistically indistinguishable) patterns from independent data sets.

Congruence studies are most advanced for echinoids, where four independent data sets (two morphological and two molecular) provide strong corroboration for a single phylogenetic scheme. Analysis of all four data sets combined generates a highly robust hypothesis of relationships. The situation is very different for asteroids. Two analyses based on morphological data have reached very different conclusions. Three independent molecular data sets also have been compiled, but none has a statistically reliable signal concerning higher taxon relationships. Even combining all three molecular data sets fails to generate a statistically robust solution, implying that the major lines of asteroids diverged rapidly from one another. For ophiuroids, both morphological and molecular data generate topologies that for the most part lack statistical robustness. There is currently no cladistic analysis of holothurian relationships based on morphological data, and only a few taxa have been sequenced. The molecular data is, however, congruent and does permit an initial assessment of relationships. Nothing definite can be deduced about crinoid relationships since even fewer molecular sequences are known and morphological analysis remains sketchy.

Class-level relationships derived from two morphological and two molecular data sets also show considerable congruence, though a single definitive solution has yet to emerge.

Rates and relations of mitochondrial genome evolution across the Echinoidea, with special focus on the superfamily Odontophora

In order to better characterize the placement of genus Tripneustes, as a representative of the Toxopneustidae family within the broader sea urchin mitochondrial (MT) phylogeny, the complete MT genome of Tripneustes gratilla was generated and compared with all published echinoid MT genomes currently available on NCBI GenBank. The MT genome phylogeny supports the existence of the superfamily Odontophora (consisting of the families Strongylocentrotidae, Echinometridae, and Toxopneustidae). A relaxed molecular‐clock time calibration suggests a split between the three key Odontophore MT lineages occurred during the late Eocene/Oligocene. Major global oceanographic changes have been inferred during this time frame, potentially driving species diversification through environmental selection pressures. To test for signatures of selection acting on the mitochondria, the historical rate of gene evolution of individual MT genes was assessed through a branch‐site comparison of nonsynonymous to synonymous substitution ratios (ω). Models of positive selection and neutral evolution, as compared via a likelihood ratio test, show no evidence of strong historical positive selection on mitochondrial genes at the genesis of the Odontophora. However, while pairwise ω comparison revealed signatures of strong negative selection, relatively elevated ω values were observed within the Strongylocentrotus genus.

Diversification of spatiotemporal expression and copy number variation of the echinoid hbox12/pmar1/micro1 multigene family

Changes occurring during evolution in the cis-regulatory landscapes of individual members of multigene families might impart diversification in their spatiotemporal expression and function. The archetypal member of the echinoid hbox12/pmar1/micro1 family is hbox12-a, a homeobox-containing gene expressed exclusively by dorsal blastomeres, where it governs the dorsal/ventral gene regulatory network during embryogenesis of the sea urchin Paracentrotus lividus. Here we describe the inventory of the hbox12/pmar1/micro1 genes in P. lividus, highlighting that gene copy number variation occurs across individual sea urchins of the same species. We show that the various hbox12/pmar1/micro1 genes group into three subfamilies according to their spatiotemporal expression, which ranges from broad transcription throughout development to transient expression in either the animal hemisphere or micromeres of the early embryo. Interestingly, the promoter regions of those genes showing comparable expression patterns are highly similar, while differing from those of the other subfamilies. Strikingly, phylogenetic analysis suggests that the hbox12/pmar1/micro1 genes are species-specific, exhibiting extensive divergence in their noncoding, but not in their coding, sequences across three distinct sea urchin species. In spite of this, two micromere-specific genes of P. lividus possess a TCF/LEF-binding motif in a similar position, and their transcription relies on Wnt/β-catenin signaling, similar to the pmar1 and micro1 genes, which in other sea urchin species are involved in micromere specification. Altogether, our findings suggest that the hbox12/pmar1/micro1 gene family evolved rather rapidly, generating paralogs whose cis-regulatory sequences diverged following multiple rounds of duplication from a common ancestor.

Best Dressed Test: A Study of the Covering Behavior of the Collector Urchin Tripneustes gratilla

Many sea urchin genera exhibit cryptic covering behaviors. One such behavior has been documented in the sea urchin Tripneustes gratilla, and previous studies have theorized that this behavior serves as protection from UV radiation. However, other hypotheses have been presented such as protection from predators or added weight to help T. gratilla resist strong currents. A field study was conducted in October-November 2015 in Moorea, French Polynesia to assess urchin covering behavior in natural habitats. The study found that urchins partially underneath rocks covered more, and with more algae, than urchins totally underneath rocks. To test if this behavior was driven by light intensity, a series of 30-minute experimental trials were run on 10 individuals in bright and dim conditions. Individuals were given red and clear plastic, and percent cover of each was recorded. These tests were repeated once fifty percent of spines had been removed from the urchin, in order to determine whether spine loss affects T. gratilla covering behavior. The study found that urchins had a distinct preference for cover that best protects them from UV radiation. Spine loss did not significantly affect urchin ability to cover, and urchins with removed spines still preferred opaque cover. Additionally, covering behavior was mapped onto a phylogeny of echinoderms to determine how it might have evolved. Understanding urchin covering behavior more fully is a step towards an understanding of the evolution of cryptic behavior across species.

Phylogeography of the sand dollar genus Mellita: cryptic speciation along the coasts of the Americas

Sand dollars of the genus Mellita are members of the sandy shallow-water fauna. The genus ranges in tropical and subtropical regions on the two coasts of the Americas. To reconstruct the phylogeography of the genus we sequenced parts of the mitochondrial cytochrome oxidase I and of 16S rRNA as well as part of the nuclear 28S rRNA gene from a total of 185 specimens of all ten described morphospecies from 31 localities. Our analyses revealed the presence of eleven species, including six cryptic species. Sequences of five morphospecies do not constitute monophyletic molecular units and thus probably represent ecophenotypic variants. The fossil-calibrated phylogeny showed that the ancestor of Mellita diverged into a Pacific lineage and an Atlantic+Pacific lineage close to the Miocene/Pliocene boundary. Atlantic M. tenuis, M. quinquiesperforata and two undescribed species of Mellita have non-overlapping distributions. Pacific Mellita consist of two highly divergent lineages that became established at different times, resulting in sympatric M. longifissa and M. notabilis. Judged by modern day ranges, not all divergence in this genus conforms to an allopatric speciation model. Only the separation of M. quinquiesperforata from M. notabilis is clearly due to vicariance as the result of the completion of the Isthmus of Panama. The molecular phylogeny calibrated on fossil evidence estimated this event as having occurred ~3 Ma, thus providing evidence that, contrary to a recent proposal, the central American Isthmus was not completed until this date.

Phylogenomics of strongylocentrotid sea urchins

Background

Strongylocentrotid sea urchins have a long tradition as model organisms for studying many fundamental processes in biology including fertilization, embryology, development and genome regulation but the phylogenetic relationships of the group remain largely unresolved. Although the differing isolating mechanisms of vicariance and rapidly evolving gamete recognition proteins have been proposed, a stable and robust phylogeny is unavailable.

Results

We used a phylogenomic approach with mitochondrial and nuclear genes taking advantage of the whole-genome sequencing of nine species in the group to establish a stable (i.e. concordance in tree topology among multiple lies of evidence) and robust (i.e. high nodal support) phylogenetic hypothesis for the family Strongylocentrotidae. We generated eight draft mitochondrial genome assemblies and obtained 13 complete mitochondrial genes for each species. Consistent with previous studies, mitochondrial sequences failed to provide a reliable phylogeny. In contrast, we obtained a very well-supported phylogeny from 2301 nuclear genes without evidence of positive Darwinian selection both from the majority of most-likely gene trees and the concatenated fourfold degenerate sites: ((P. depressus, (M. nudus, M. franciscanus), (H. pulcherrimus, (S. purpuratus, (S. fragilis, (S. pallidus, (S. droebachiensis, S. intermedius)). This phylogeny was consistent with a single invasion of deep-water environments followed by a holarctic expansion by Strongylocentrotus. Divergence times for each species estimated with reference to the divergence times between the two major clades of the group suggest a correspondence in the timing with the opening of the Bering Strait and the invasion of the holarctic regions.

Conclusions

Nuclear genome data contains phylogenetic signal informative for understanding the evolutionary history of this group. However, mitochondrial genome data does not. Vicariance can explain major patterns observed in the phylogeny. Other isolating mechanisms are appropriate to explore in this system to help explain divergence patterns not well supported by vicariance, such as the effects of rapidly evolving gamete recognition proteins on isolating populations. Our findings of a stable and robust phylogeny, with the increase in mitochondrial and nuclear comparative genomic data, provide a system in which we can enhance our understanding of molecular evolution and adaptation in this group of sea urchins.

Population genetics of cis-regulatory sequences that operate during embryonic development in the sea urchin Strongylocentrotus purpuratus

Despite the fact that noncoding sequences comprise a substantial fraction of functional sites within all genomes, the evolutionary mechanisms that operate on genetic variation within regulatory elements remain poorly understood. In this study, we examine the population genetics of the core, upstream cis-regulatory regions of eight genes (AN, CyIIa, CyIIIa, Endo16, FoxB, HE, SM30 a, and SM50) that function during the early development of the purple sea urchin, Strongylocentrotus purpuratus. Quantitative and qualitative measures of segregating variation are not conspicuously different between cis-regulatory and closely linked "proxy neutral" noncoding regions containing no known functional sites. Length and compound mutations are common in noncoding sequences; conventional descriptive statistics ignore such mutations, under-representing true genetic variation by approximately 28% for these loci in this population. Patterns of variation in the cis-regulatory regions of six of the genes examined (CyIIa, CyIIIa, Endo16, FoxB, AN, and HE) are consistent with directional selection. Genetic variation within annotated transcription factor binding sites is comparable to, and frequently greater than, that of surrounding sequences. Comparisons of two paralog pairs (CyIIa/CyIIIa and AN/HE) suggest that distinct evolutionary processes have operated on their cis-regulatory regions following gene duplication. Together, these analyses provide a detailed view of the evolutionary mechanisms operating on noncoding sequences within a natural population, and underscore how little is known about how these processes operate on cis-regulatory sequences.

Fixed, free, and fixed: the fickle phylogeny of extant Crinoidea (Echinodermata) and their Permian-Triassic origin

Although the status of Crinoidea (sea lilies and featherstars) as sister group to all other living echinoderms is well-established, relationships among crinoids, particularly extant forms, are debated. All living species are currently placed in Articulata, which is generally accepted as the only crinoid group to survive the Permian-Triassic extinction event. Recent classifications have recognized five major extant taxa: Isocrinida, Hyocrinida, Bourgueticrinina, Comatulidina and Cyrtocrinida, plus several smaller groups with uncertain taxonomic status, e.g., Guillecrinus, Proisocrinus and Caledonicrinus. Here we infer the phylogeny of extant Crinoidea using three mitochondrial genes and two nuclear genes from 59 crinoid terminals that span the majority of extant crinoid diversity. Although there is poor support for some of the more basal nodes, and some tree topologies varied with the data used and mode of analysis, we obtain several robust results. Cyrtocrinida, Hyocrinida, Isocrinida are all recovered as clades, but two stalked crinoid groups, Bourgueticrinina and Guillecrinina, nest among the featherstars, lending support to an argument that they are paedomorphic forms. Hence, they are reduced to families within Comatulida. Proisocrinus is clearly shown to be part of Isocrinida, and Caledonicrinus may not be a bourgueticrinid. Among comatulids, tree topologies show little congruence with current taxonomy, indicating that much systematic revision is required. Relaxed molecular clock analyses with eight fossil calibration points recover Articulata with a median date to the most recent common ancestor at 231-252mya in the Middle to Upper Triassic. These analyses tend to support the hypothesis that the group is a radiation from a small clade that passed through the Permian-Triassic extinction event rather than several lineages that survived. Our tree topologies show various scenarios for the evolution of stalks and cirri in Articulata, so it is clear that further data and taxon sampling are needed to recover a more robust phylogeny of the group.

Evolution of a novel muscle design in sea urchins (Echinodermata: Echinoidea)

The sea urchin (Echinodermata: Echinoidea) masticatory apparatus, or Aristotle's lantern, is a complex structure composed of numerous hard and soft components. The lantern is powered by various paired and unpaired muscle groups. We describe how one set of these muscles, the lantern protractor muscles, has evolved a specialized morphology. This morphology is characterized by the formation of adaxially-facing lobes perpendicular to the main orientation of the muscle, giving the protractor a frilled aspect in horizontal section. Histological and ultrastructural analyses show that the microstructure of frilled muscles is largely identical to that of conventional, flat muscles. Measurements of muscle dimensions in equally-sized specimens demonstrate that the frilled muscle design, in comparison to that of the flat muscle type, considerably increases muscle volume as well as the muscle's surface directed towards the interradial cavity, a compartment of the peripharyngeal coelom. Scanning electron microscopical observations reveal that the insertions of frilled and flat protractor muscles result in characteristic muscle scars on the stereom, reflecting the shapes of individual muscles. Our comparative study of 49 derived "regular" echinoid species using magnetic resonance imaging (MRI) shows that frilled protractor muscles are found only in taxa belonging to the families Toxopneustidae, Echinometridae, and Strongylocentrotidae. The onset of lobe formation during ontogenesis varies between species of these three families. Because frilled protractor muscles are best observed in situ, the application of a non-invasive imaging technique was crucial for the unequivocal identification of this morphological character on a large scale. Although it is currently possible only to speculate on the functional advantages which the frilled muscle morphology might confer, our study forms the anatomical and evolutionary framework for future analyses of this unusual muscle design among sea urchins.

Specific expression of a TRIM-containing factor in ectoderm cells affects the skeletal morphogenetic program of the sea urchin embryo

In the indirect developing sea urchin embryo, the primary mesenchyme cells (PMCs) acquire most of the positional and temporal information from the overlying ectoderm for skeletal initiation and growth. In this study, we characterize the function of the novel gene strim1, which encodes a tripartite motif-containing (TRIM) protein, that adds to the list of genes constituting the epithelial-mesenchymal signaling network. We report that strim1 is expressed in ectoderm regions adjacent to the bilateral clusters of PMCs and that its misexpression leads to severe skeletal abnormalities. Reciprocally, knock down of strim1 function abrogates PMC positioning and blocks skeletogenesis. Blastomere transplantation experiments establish that the defects in PMC patterning, number and skeletal growth depend upon strim1 misexpression in ectoderm cells. Furthermore, clonal expression of strim1 into knocked down embryos locally restores skeletogenesis. We also provide evidence that the Otp and Pax2/5/8 regulators, as well as FGFA, but not VEGF, ligand act downstream to strim1 in ectoderm cells, and that strim1 triggers the expression of the PMC marker sm30, an ectoderm-signaling dependent gene. We conclude that the strim1 function elicits specific gene expression both in ectoderm cells and PMCs to guide the skeletal biomineralization during morphogenesis.

Origin and evolutionary plasticity of the gastric caecum in sea urchins (Echinodermata: Echinoidea)

Background

The digestive tract of many metazoan invertebrates is characterized by the presence of caeca or diverticula that serve secretory and/or absorptive functions. With the development of various feeding habits, distinctive digestive organs may be present in certain taxa. This also holds true for sea urchins (Echinodermata: Echinoidea), in which a highly specialized gastric caecum can be found in members of a derived subgroup, the Irregularia (cake urchins, sea biscuits, sand dollars, heart urchins, and related forms). As such a specialized caecum has not been reported from "regular" sea urchin taxa, the aim of this study was to elucidate its evolutionary origin.

Results

Using morphological data derived from dissection, magnetic resonance imaging, and extensive literature studies, we compare the digestive tract of 168 echinoid species belonging to 51 extant families. Based on a number of characters such as topography, general morphology, mesenterial suspension, and integration into the haemal system, we homologize the gastric caecum with the more or less pronounced dilation of the anterior stomach that is observed in most "regular" sea urchin taxa. In the Irregularia, a gastric caecum can be found in all taxa except in the Laganina and Scutellina. It is also undeveloped in certain spatangoid species.

Conclusions

According to our findings, the sea urchin gastric caecum most likely constitutes a synapomorphy of the Euechinoidea. Its occurrence in "regular" euechinoids is linked to the presence of an additional festoon of the anterior stomach in ambulacrum III. Both structures, the additional festoon and the gastric caecum, are absent in the sister taxon to the Euechinoidea, the Cidaroida. Since the degree of specialization of the gastric caecum is most pronounced in the predominantly sediment-burrowing irregular taxa, we hypothesize that its evolution is closely linked to the development of more elaborate infaunal lifestyles. We provide a comprehensive study of the origin and evolutionary plasticity of a conspicuous digestive tract structure, the gastric caecum, in a major taxon of the extant invertebrate macrozoobenthos.

Sniffing out new data and hypotheses on the form, function, and evolution of the echinopluteus post-oral vibratile lobe

The performance requirements of ciliary band feeding explain the convoluted forms of many marine invertebrate larvae. Convolutions increase surface area and therefore feeding rates per unit body volume. We review recent advances in morphology, neural development, and behavior at settlement of the echinoid Lytechinus pictus and provide new ultrastructural and expression data on larvae of its congener, L. variegatus. Larvae of the echinometrid Colobocentrotus atratus contain neurons identified by their expression of nitric oxide synthase (NOS), indicating that this character is not unique to Lytechinus. We hypothesize that in some echinoids the convoluted shape of the post-oral vibratile lobe (POVL) covaries with the distribution of identified sensory neurons to enable olfaction during settlement. An analysis of variation in structural elaboration of the post-oral transverse ciliary band (PTB) within Echinoida and in feeding larvae of other echinoderm classes indicates that only echinoids, but not all echinoids, possess this novel character; larvae that do are distributed heterogeneously within the class. In recognition of this specialized function for the POVL and surrounding ectoderm, and because it is lobate and grows toward the mouth, we propose naming this structure the adoral lobe.

The major yolk protein is synthesized in the digestive tract and secreted into the body cavities in sea urchin larvae

Major yolk protein (MYP), a transferrin superfamily protein contained in yolk granules of sea urchin eggs, also occurs in the coelomic fluid of male and female adult sea urchins regardless of their reproductive cycle. MYP in the coelomic fluid (CFMYP; 180 kDa) has a zinc-binding capacity and has a higher molecular mass than MYP in eggs (EGMYP; 170 kDa). CFMYP is thought to be synthesized in the digestive tract and secreted into the coelomic fluid where it is involved in the transport of zinc derived from food. To clarify when and where MYP synthesis starts, we investigated the expression of MYP during larval development and growth in Pseudocentrotus depressus. MYP mRNA was detected using RT-PCR in the early 8-arm pluteus stage and its expression persisted until after metamorphosis. Real-time RT-PCR revealed that MYP mRNA increased exponentially from the early 8-arm stage to metamorphosis. Western blotting showed that maternal EGMYP disappeared by the 4-arm stage and that newly synthesized CFMYP was present at and after the mid 8-arm stage. In the late 8-arm larvae, MYP mRNA was detected in the digestive tract using in situ hybridization, and the protein was found in the somatocoel and the blastocoel-derived space between the somatocoel and epidermis using immunohistochemistry. These results suggest that CFMYP is synthesized in the digestive tract and secreted into the body cavities at and after the early 8-arm stage. We assume that in larvae, CFMYP transports zinc derived from food via the body cavities to various tissues, as suggested for adults.

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.

Systematic comparison and reconstruction of sea urchin (Echinoidea) internal anatomy: a novel approach using magnetic resonance imaging

Background

Traditional comparative morphological analyses and subsequent three-dimensional reconstructions suffer from a number of drawbacks. This is particularly evident in the case of soft tissue studies that are technically demanding, time-consuming, and often prone to produce artefacts. These problems can partly be overcome by employing non-invasive, destruction-free imaging techniques, in particular micro-computed tomography or magnetic resonance imaging.

Results

Here, we employed high-field magnetic resonance imaging techniques to gather numerous data from members of a major marine invertebrate taxon, the sea urchins (Echinoidea). For this model study, 13 of the 14 currently recognized high-ranking subtaxa (orders) of this group of animals were analyzed. Based on the acquired datasets, interactive three-dimensional models were assembled. Our analyses reveal that selected soft tissue characters can even be used for phylogenetic inferences in sea urchins, as exemplified by differences in the size and shape of the gastric caecum found in the Irregularia.

Conclusion

The main focus of our investigation was to explore the possibility to systematically visualize the internal anatomy of echinoids obtained from various museum collections. We show that, in contrast to classical preparative procedures, magnetic resonance imaging can give rapid, destruction-free access to morphological data from numerous specimens, thus extending the range of techniques available for comparative studies of invertebrate morphology.

A method for constraining the age of origination of derived characters

Fossils are the physical records of the history of morphological character evolution on Earth and can provide valuable information concerning the sequence and timing of origination of derived characters. Knowledge of the timing of origination of synapomorphies makes it possible to estimate when unobserved character changes occurred in the geological past. Here we present a method for estimating the temporal interval during which synapomorphies evolved. The method requires either direct inclusion of fossil taxa (with or without extant taxa) in cladistic analyses based on morphological or combined data, or indirectly using the "molecular scaffold approach." Second, characters of interest are mapped on a most parsimonious tree and "minimum age node mapping" is used to place minimum ages on the nodes of the tree. Finally, characters of interest are evaluated for younger and/or older temporal constraints on the time of their origination; application of the older bound assumes ancestry of fossil terminals included in the tree. A key is provided herein describing the method. Among other applications, this approach has the potential to provide a powerful test of purported evolutionary cause-effect relationships. For example, the method has the ability to discover that derived characters of suggested adaptational significance may considerably pre-date the cause(s) that are hypothesized to have favored their establishment.

Gene expression patterns in a novel animal appendage: the sea urchin pluteus arm

The larval arms of echinoid plutei are used for locomotion and feeding. They are composed of internal calcite skeletal rods covered by an ectoderm layer bearing a ciliary band. Skeletogenesis includes an autonomous molecular differentiation program in primary mesenchyme cells (PMCs), initiated when PMCs leave the vegetal plate for the blastocoel, and a patterning of the differentiated skeletal units that requires molecular cues from the overlaying ectoderm. The arms represent a larval feature that arose in the echinoid lineage during the Paleozoic and offers a subject for the study of gene co-option in the evolution of novel larval features. We isolated new molecular markers in two closely related but differently developing species, Heliocidaris tuberculata and Heliocidaris erythrogramma. We report the expression of a larval arm-associated ectoderm gene tetraspanin, as well as two new PMC markers, advillin and carbonic anhydrase. Tetraspanin localizes to the animal half of blastula stage H. tuberculata and then undergoes a restriction into the putative oral ectoderm and future location of the postoral arms, where it continues to be expressed at the leading edge of both the postoral and anterolateral arms. In H. erythrogramma, its expression initiates in the animal half of blastulae and expands over the entire ectoderm from gastrulation onward. Advillin and carbonic anhydrase are upregulated in the PMCs postgastrulation and localized to the leading edge of the growing larval arms of H. tuberculata but do not exhibit coordinated expression in H. erythrogramma larvae. The tight spatiotemporal regulation of these genes in H. tuberculata along with other ontogenetic and phylogenetic evidence suggest that pluteus arms are novel larval organs, distinguishable from the processes of skeletogenesis per se. The dissociation of expression control in H. erythrogramma suggest that coordinate gene expression in H. tuberculata evolved as part of the evolution of pluteus arms, and is not required for larval or adult development.

Testing the molecular clock: molecular and paleontological estimates of divergence times in the Echinoidea (Echinodermata)

The phylogenetic relationships of 46 echinoids, with representatives from 13 of the 14 ordinal-level clades and about 70% of extant families commonly recognized, have been established from 3 genes (3,226 alignable bases) and 119 morphological characters. Morphological and molecular estimates are similar enough to be considered suboptimal estimates of one another, and the combined data provide a tree that, when calibrated against the fossil record, provides paleontological estimates of divergence times and completeness of their fossil record. The order of branching on the cladogram largely agrees with the stratigraphic order of first occurrences and implies that their fossil record is more than 85% complete at family level and at a resolution of 5-Myr time intervals. Molecular estimates of divergence times derived from applying both molecular clock and relaxed molecular clock models are concordant with estimates based on the fossil record in up to 70% of cases, with most concordant results obtained using Sanderson's semiparametric penalized likelihood method and a logarithmic-penalty function. There are 3 regions of the tree where molecular and fossil estimates of divergence time consistently disagree. Comparison with results obtained when molecular divergence dates are estimated from the combined (morphology + gene) tree suggests that errors in phylogenetic reconstruction explain only one of these. In another region the error most likely lies with the paleontological estimates because taxa in this region are demonstrated to have a very poor fossil record. In the third case, morphological and paleontological evidence is much stronger, and the topology for this part of the molecular tree differs from that derived from the combined data. Here the cause of the mismatch is unclear but could be methodological, arising from marked inequality of molecular rates. Overall, the level of agreement reached between these different data and methodological approaches leads us to believe that careful application of likelihood and Bayesian methods to molecular data provides realistic divergence time estimates in the majority of cases (almost 80% in this specific example), thus providing a remarkably well-calibrated phylogeny of a character-rich clade of ubiquitous marine benthic invertebrates.

Multiple origins of the ascidian-Prochloron symbiosis: molecular phylogeny of photosymbiotic and non-symbiotic colonial ascidians inferred from 18S rDNA sequences

In the tropics, certain didemnid ascidians harbor the prokaryotic photosymbiont Prochloron. To date, this photosymbiosis has been found in four didemnid genera that include non-symbiotic species. Here, we report the molecular phylogeny of symbiotic and non-symbiotic didemnids based on their 18S rDNA sequences. The data cover all four genera containing symbiotic species and one other genus comprised of only non-symbiotic species. Near-complete nucleotide sequences of 18S rDNAs were determined for four non-didemnid species and 52 didemnid samples (five genera), including 48 photosymbiotic samples collected from the Ryukyu Archipelago, the Great Barrier Reef, Hawaii, and Bali. Our phylogenetic trees indicated a monophyletic origin of the family Didemnidae, as well as each of the didemnid genera. The results strongly support the hypothesis that establishment of the ascidian-Prochloron symbiosis occurred independently in the Didemnidae lineage at least once in each of the genera that possess symbiotic species.

Systematics of Chaetognatha under the light of molecular data, using duplicated ribosomal 18S DNA sequences

While the phylogenetic position of Chaetognatha has became central to the question of early bilaterian evolution, the internal systematics of the phylum are still not clear. The phylogenetic relationships of the chaetognaths were investigated using newly obtained small subunit ribosomal RNA nuclear 18S (SSU rRNA) sequences from 16 species together with 3 sequences available in GenBank. As previously shown with the large subunit ribosomal RNA 28S gene, two classes of Chaetognatha SSU rRNA gene can be identified, suggesting a duplication of the whole ribosomal cluster; allowing the rooting of one class of genes by another in phylogenetic analyses. Maximum Parsimony, Maximum Likelihood and Bayesian analyses of the molecular data, and statistical tests showed (1) that there are three main monophyletic groups: Sagittidae/Krohnittidae, Spadellidae/Pterosagittidae, and Eukrohniidae/Heterokrohniidae, (2) that the group of Aphragmophora without Pterosagittidae (Sagittidae/Krohnittidae) is monophyletic, (3) the Spadellidae/Pterosagittidae and Eukrohniidae/Heterokrohniidae families are very likely clustered, (4) the Krohnittidae and Pterosagittidae groups should no longer be considered as families as they are included in other groups designated as families, (5) suborder Ctenodontina is not monophyletic and the Flabellodontina should no longer be considered as a suborder, and (6) the Syngonata/Chorismogonata and the Monophragmophora/Biphragmophora hypotheses are rejected. Such conclusions are considered in the light of morphological characters, several of which are shown to be prone to homoplasy.

Larval ectoderm, organizational homology, and the origins of evolutionary novelty

Comprehending the origin of marine invertebrate larvae remains a key domain of research for evolutionary biologists, including the repeated origin of direct developmental modes in echinoids. In order to address the latter question, we surveyed existing evidence on relationships of homology between the ectoderm territories of two closely related sea urchin species in the genus Heliocidaris that differ in their developmental mode. Additionally, we explored a recently articulated idea about homology called 'organizational homology' (Müller 2003. In: Müller GB, Newman SA, editors. Origination of organismal form: beyond the gene in developmental and evolutionary biology. Cambridge, MA: A Bradford Book, The MIT Press. p 51-69. ) in the context of this specific empirical case study. Applying the perspective of organizational homology to our experimental system of congeneric echinoids has led us to a new hypothesis concerning the ectoderm evolution in these species. The extravestibular ectoderm of the direct developer Heliocidaris erythrogramma is a novel developmental territory that arose as a fusion of the oral and aboral ectoderm territories found in indirect developing echinoids such as Heliocidaris tuberculata. This hypothesis instantiates a theoretical principle concerning the origin of developmental modules, 'integration', which has been neglected because the opposite theoretical principle, 'parcellation', is more readily observable in events such as gene duplication and divergence (Wagner 1996. Am Zool 36:36-43).

An evolutionary constraint: strongly disfavored class of change in DNA sequence during divergence of cis-regulatory modules

The DNA of functional cis-regulatory modules displays extensive sequence conservation in comparisons of genomes from modestly distant species. Patches of sequence that are several hundred base pairs in length within these modules are often seen to be 80-95% identical, although the flanking sequence cannot even be aligned. However, it is unlikely that base pairs located between the transcription factor target sites of cis-regulatory modules have sequence-dependent function, and the mechanism that constrains evolutionary change within cis-regulatory modules is incompletely understood. We chose five functionally characterized cis-regulatory modules from the Strongylocentrotus purpuratus (sea urchin) genome and obtained orthologous regulatory and flanking sequences from a bacterial artificial chromosome genome library of a congener, Strongylocentrotus franciscanus. As expected, single-nucleotide substitutions and small indels occur freely at many positions within the regulatory modules of these two species, as they do outside the regulatory modules. However, large indels (>20 bp) are statistically almost absent within the regulatory modules, although they are common in flanking intergenic or intronic sequence. The result helps to explain the patterns of evolutionary sequence divergence characteristic of cis-regulatory DNA.

Structure, expression, and transcriptional regulation of the Strongylocentrotus franciscanus spec gene family encoding intracellular calcium-binding proteins

The mechanisms by which gene expression patterns emerge during evolution are poorly understood. The sea urchin spec genes offer a useful means to investigate evolutionary mechanisms. Genes of the spec family from Strongylocentrotus purpuratus and Lytechinus pictus have identical patterns of aboral ectoderm-specific expression but exhibit species-specific differences in copy number, genomic structure, temporal expression, and cis-regulatory architecture. Here, we identify spec genes from a phylogenetic intermediate, Strongylocentrotus franciscanus, to gain insight into the evolution of the spec gene family and its transcriptional regulation. We identified two spec genes in the S. franciscanus genome, sfspec1a and sfspec1b, that were orthologous to spec1 from S. purpuratus. sfspec1b transcripts began to accumulate at the blastula stage and became progressively more abundant; this was reminiscent of spec expression in L. pictus but different from that in S. purpuratus. As expected, sfspec1b expression was restricted to aboral ectoderm cells. The six-exon structure of the sfspec1b genomic locus was identical to that of the S. purpuratus spec genes and was bounded by two repeat-spacer-repeat (RSR) repetitive sequence elements, which are conserved features of S. purpuratus spec genes and function as transcriptional enhancers. The enhancer activity of the sfspec1b RSRs was comparable to that of their S. purpuratus counterparts, although the placement and orientation of crucial cis-regulatory elements within the RSRs differed. We discovered a spec gene in S. franciscanus that was only distantly related to other spec genes but was highly conserved in S. purpuratus. Unexpectedly, this gene was expressed exclusively in endoderm lineages. Our results show that the evolution of spec cis-regulatory elements is highly dynamic and that substantial alterations can occur when maintaining or grossly modifying gene expression patterns.

Creation of cis-regulatory elements during sea urchin evolution by co-option and optimization of a repetitive sequence adjacent to the spec2a gene

The creation, preservation, and degeneration of cis-regulatory elements controlling developmental gene expression are fundamental genome-level evolutionary processes about which little is known. Here, we identify critical differences in cis-regulatory elements controlling the expression of the sea urchin aboral ectoderm-specific spec genes. We found multiple copies of a repetitive sequence element termed RSR in genomes of species within the Strongylocentrotidae family, but RSRs were not detected in genomes of species outside Strongylocentrotidae. spec genes in Strongylocentrotus purpuratus are invariably associated with RSRs, and the spec2a RSR functioned as a transcriptional enhancer and displayed greater activity than did spec1 or spec2c RSRs. Single-base pair differences at two cis-regulatory elements within the spec2a RSR increased the binding affinities of four transcription factors, SpCCAAT-binding factor at one element and SpOtx, SpGoosecoid, and SpGATA-E at another. The cis-regulatory elements to which these four factors bound were recent evolutionary acquisitions that acted to either activate or repress transcription, depending on the cell type. These elements were found in the spec2a RSR ortholog in Strongylocentrotus pallidus but not in RSR orthologs of Strongylocentrotus droebachiensis or Hemicentrotus pulcherrimus. Our results indicated that a dynamic pattern of cis-regulatory element evolution exists for spec genes despite their conserved aboral ectoderm expression.