An early cell shape transition drives evolutionary expansion of the human forebrain.

The human brain has undergone rapid expansion since humans diverged from other great apes, but the mechanism of this human-specific enlargement is still unknown. Here, we use cerebral organoids derived from human, gorilla, and chimpanzee cells to study developmental mechanisms driving evolutionary brain expansion. We find that neuroepithelial differentiation is a protracted process in apes, involving a previously unrecognized transition state characterized by a change in cell shape. Furthermore, we show that human organoids are larger due to a delay in this transition, associated with differences in interkinetic nuclear migration and cell cycle length. Comparative RNA sequencing (RNA-seq) reveals differences in expression dynamics of cell morphogenesis factors, including ZEB2, a known epithelial-mesenchymal transition regulator. We show that ZEB2 promotes neuroepithelial transition, and its manipulation and downstream signaling leads to acquisition of nonhuman ape architecture in the human context and vice versa, establishing an important role for neuroepithelial cell shape in human brain expansion.

Contrasting the development of larval and adult body plans during the evolution of biphasic lifecycles in sea urchins.

Biphasic lifecycles are widespread among animals, but little is known about how the developmental transition between larvae and adults is regulated. Sea urchins are a unique system for studying this phenomenon because of the stark differences between their bilateral larval and pentaradial adult body plans. Here, we use single-cell RNA sequencing to analyze the development of Heliocidaris erythrogramma (He), a sea urchin species with an accelerated, non-feeding mode of larval development. The sequencing time course extends from embryogenesis to roughly a day before the onset of metamorphosis in He larvae, which is a period that has not been covered by previous datasets. We find that the non-feeding developmental strategy of He is associated with several changes in the specification of larval cell types compared to sea urchins with feeding larvae, such as the loss of a larva-specific skeletal cell population. Furthermore, the development of the larval and adult body plans in sea urchins may utilize largely different sets of regulatory genes. These findings lay the groundwork for extending existing developmental gene regulatory networks to cover additional stages of biphasic lifecycles.

An RNA interference approach for functional studies in the sea urchin and its use in analysis of nodal signaling gradients.

Dicer substrate interfering RNAs (DsiRNAs) destroy targeted transcripts using the RNA-Induced Silencing Complex (RISC) through a process called RNA interference (RNAi). This process is ubiquitous among eukaryotes. Here we report the utility of DsiRNA in embryos of the sea urchin Lytechinus variegatus (Lv). Specific knockdowns phenocopy known morpholino and inhibitor knockdowns, and DsiRNA offers a useful alternative to morpholinos. Methods are described for the design of specific DsiRNAs that lead to destruction of targeted mRNA. DsiRNAs directed against pks1, an enzyme necessary for pigment production, show how successful DsiRNA perturbations are monitored by RNA in situ analysis and by qPCR to determine relative destruction of targeted mRNA. DsiRNA-based knockdowns phenocopy morpholino- and drug-based inhibition of nodal and lefty. Other knockdowns demonstrate that the RISC operates early in development as well as on genes that are first transcribed hours after gastrulation is completed. Thus, DsiRNAs effectively mediate destruction of targeted mRNA in the sea urchin embryo. The approach offers significant advantages over other widely used methods in the urchin in terms of cost, and ease of procurement, and offers sizeable experimental advantages in terms of ease of handling, injection, and knockdown validation.

Hybrid Epigenomes Reveal Extensive Local Genetic Changes to Chromatin Accessibility Contribute to Divergence in Embryonic Gene Expression Between Species.

Chromatin accessibility plays an important role in shaping gene expression, yet little is known about the genetic and molecular mechanisms that influence the evolution of chromatin configuration. Both local (cis) and distant (trans) genetic influences can in principle influence chromatin accessibility and are based on distinct molecular mechanisms. We, therefore, sought to characterize the role that each of these plays in altering chromatin accessibility in 2 closely related sea urchin species. Using hybrids of Heliocidaris erythrogramma and Heliocidaris tuberculata, and adapting a statistical framework previously developed for the analysis of cis and trans influences on the transcriptome, we examined how these mechanisms shape the regulatory landscape at 3 important developmental stages, and compared our results to similar analyses of the transcriptome. We found extensive cis- and trans-based influences on evolutionary changes in chromatin, with cis effects generally larger in effect. Evolutionary changes in accessibility and gene expression are correlated, especially when expression has a local genetic basis. Maternal influences appear to have more of an effect on chromatin accessibility than on gene expression, persisting well past the maternal-to-zygotic transition. Chromatin accessibility near gene regulatory network genes appears to be distinctly regulated, with trans factors appearing to play an outsized role in the configuration of chromatin near these genes. Together, our results represent the first attempt to quantify cis and trans influences on evolutionary divergence in chromatin configuration in an outbred natural study system and suggest that chromatin regulation is more genetically complex than was previously appreciated.

Developmental single-cell transcriptomics in the Lytechinus variegatus sea urchin embryo.

Using scRNA-seq coupled with computational approaches, we studied transcriptional changes in cell states of sea urchin embryos during development to the larval stage. Eighteen closely spaced time points were taken during the first 24 h of development of Lytechinus variegatus (Lv). Developmental trajectories were constructed using Waddington-OT, a computational approach to 'stitch' together developmental time points. Skeletogenic and primordial germ cell trajectories diverged early in cleavage. Ectodermal progenitors were distinct from other lineages by the 6th cleavage, although a small percentage of ectoderm cells briefly co-expressed endoderm markers that indicated an early ecto-endoderm cell state, likely in cells originating from the equatorial region of the egg. Endomesoderm cells also originated at the 6th cleavage and this state persisted for more than two cleavages, then diverged into distinct endoderm and mesoderm fates asynchronously, with some cells retaining an intermediate specification status until gastrulation. Seventy-nine out of 80 genes (99%) examined, and included in published developmental gene regulatory networks (dGRNs), are present in the Lv-scRNA-seq dataset and are expressed in the correct lineages in which the dGRN circuits operate.

Genomic insights into Plasmodium vivax population structure and diversity in central Africa.

BACKGROUND: Though Plasmodium vivax is the second most common malaria species to infect humans, it has not traditionally been considered a major human health concern in central Africa given the high prevalence of the human Duffy-negative phenotype that is believed to prevent infection. Increasing reports of asymptomatic and symptomatic infections in Duffy-negative individuals throughout Africa raise the possibility that P. vivax is evolving to evade host resistance, but there are few parasite samples with genomic data available from this part of the world. METHODS: Whole genome sequencing of one new P. vivax isolate from the Democratic Republic of the Congo (DRC) was performed and used in population genomics analyses to assess how this central African isolate fits into the global context of this species. RESULTS: Plasmodium vivax from DRC is similar to other African populations and is not closely related to the non-human primate parasite P. vivax-like. Evidence is found for a duplication of the gene PvDBP and a single copy of PvDBP2. CONCLUSION: These results suggest an endemic P. vivax population is present in central Africa. Intentional sampling of P. vivax across Africa would further contextualize this sample within African P. vivax diversity and shed light on the mechanisms of infection in Duffy negative individuals. These results are limited by the uncertainty of how representative this single sample is of the larger population of P. vivax in central Africa.

Microbiome reduction and endosymbiont gain from a switch in sea urchin life history.

Animal gastrointestinal tracts harbor a microbiome that is integral to host function, yet species from diverse phyla have evolved a reduced digestive system or lost it completely. Whether such changes are associated with alterations in the diversity and/or abundance of the microbiome remains an untested hypothesis in evolutionary symbiosis. Here, using the life history transition from planktotrophy (feeding) to lecithotrophy (nonfeeding) in the sea urchin Heliocidaris, we demonstrate that the lack of a functional gut corresponds with a reduction in microbial community diversity and abundance as well as the association with a diet-specific microbiome. We also determine that the lecithotroph vertically transmits a Rickettsiales that may complement host nutrition through amino acid biosynthesis and influence host reproduction. Our results indicate that the evolutionary loss of a functional gut correlates with a reduction in the microbiome and the association with an endosymbiont. Symbiotic transitions can therefore accompany life history transitions in the evolution of developmental strategies.

Evolutionary Changes in the Chromatin Landscape Contribute to Reorganization of a Developmental Gene Network During Rapid Life History Evolution in Sea Urchins.

Chromatin configuration is highly dynamic during embryonic development in animals, exerting an important point of control in transcriptional regulation. Yet there exists remarkably little information about the role of evolutionary changes in chromatin configuration to the evolution of gene expression and organismal traits. Genome-wide assays of chromatin configuration, coupled with whole-genome alignments, can help address this gap in knowledge in several ways. In this study we present a comparative analysis of regulatory element sequences and accessibility throughout embryogenesis in three sea urchin species with divergent life histories: a lecithotroph Heliocidaris erythrogramma, a closely related planktotroph H. tuberculata, and a distantly related planktotroph Lytechinus variegatus. We identified distinct epigenetic and mutational signatures of evolutionary modifications to the function of putative cis-regulatory elements in H. erythrogramma that have accumulated nonuniformly throughout the genome, suggesting selection, rather than drift, underlies many modifications associated with the derived life history. Specifically, regulatory elements composing the sea urchin developmental gene regulatory network are enriched for signatures of positive selection and accessibility changes which may function to alter binding affinity and access of developmental transcription factors to these sites. Furthermore, regulatory element changes often correlate with divergent expression patterns of genes involved in cell type specification, morphogenesis, and development of other derived traits, suggesting these evolutionary modifications have been consequential for phenotypic evolution in H. erythrogramma. Collectively, our results demonstrate that selective pressures imposed by changes in developmental life history rapidly reshape the cis-regulatory landscape of core developmental genes to generate novel traits and embryonic programs.

The hard clam genome reveals massive expansion and diversification of inhibitors of apoptosis in Bivalvia.

BACKGROUND: Inhibitors of apoptosis (IAPs) are critical regulators of programmed cell death that are essential for development, oncogenesis, and immune and stress responses. However, available knowledge regarding IAP is largely biased toward humans and model species, while the distribution, function, and evolutionary novelties of this gene family remain poorly understood in many taxa, including Mollusca, the second most speciose phylum of Metazoa. RESULTS: Here, we present a chromosome-level genome assembly of an economically significant bivalve, the hard clam Mercenaria mercenaria, which reveals an unexpected and dramatic expansion of the IAP gene family to 159 members, the largest IAP gene repertoire observed in any metazoan. Comparative genome analysis reveals that this massive expansion is characteristic of bivalves more generally. Reconstruction of the evolutionary history of molluscan IAP genes indicates that most originated in early metazoans and greatly expanded in Bivalvia through both lineage-specific tandem duplication and retroposition, with 37.1% of hard clam IAPs located on a single chromosome. The expanded IAPs have been subjected to frequent domain shuffling, which has in turn shaped their architectural diversity. Further, we observed that extant IAPs exhibit dynamic and orchestrated expression patterns among tissues and in response to different environmental stressors. CONCLUSIONS: Our results suggest that sophisticated regulation of apoptosis enabled by the massive expansion and diversification of IAPs has been crucial for the evolutionary success of hard clam and other molluscan lineages, allowing them to cope with local environmental stresses. This study broadens our understanding of IAP proteins and expression diversity and provides novel resources for studying molluscan biology and IAP function and evolution.

TBX5-encoded T-box transcription factor 5 variant T223M is associated with long QT syndrome and pediatric sudden cardiac death.

Long QT syndrome (LQTS) is a genetic disease resulting in a prolonged QT interval on a resting electrocardiogram, predisposing affected individuals to polymorphic ventricular tachycardia and sudden death. Although a number of genes have been implicated in this disease, nearly one in four individuals exhibiting the LQTS phenotype are genotype-negative. Whole-exome sequencing identified a missense T223M variant in TBX5 that cosegregates with prolonged QT interval in a family with otherwise genotype-negative LQTS and sudden death. The TBX5-T223M variant was absent among large ostensibly healthy populations (gnomAD) and predicted to be pathogenic by in silico modeling based on Panther, PolyPhen-2, Provean, SIFT, SNAP2, and PredictSNP prediction tools. The variant was located in a highly conserved region of TBX5 predicted to be part of the DNA-binding interface. A luciferase assay identified a 57.5% reduction in the ability of TBX5-T223M to drive expression at the atrial natriuretic factor promotor compared to wildtype TBX5 in vitro. We conclude that the variant is pathogenic in this family, and we put TBX5 forward as a disease susceptibility allele for genotype-negative LQTS. The identification of this familial variant may serve as a basis for the identification of previously unknown mechanisms of LQTS with broader implications for cardiac electrophysiology.

Identifying branch-specific positive selection throughout the regulatory genome using an appropriate proxy neutral.

BACKGROUND: Adaptive changes in cis-regulatory elements are an essential component of evolution by natural selection. Identifying adaptive and functional noncoding DNA elements throughout the genome is therefore crucial for understanding the relationship between phenotype and genotype. RESULTS: We used ENCODE annotations to identify appropriate proxy neutral sequences and demonstrate that the conservativeness of the test can be modulated during the filtration of reference alignments. We applied the method to noncoding Human Accelerated Elements as well as open chromatin elements previously identified in 125 human tissues and cell lines to demonstrate its utility. Then, we evaluated the impact of query region length, proxy neutral sequence length, and branch count on test sensitivity and specificity. We found that the length of the query alignment can vary between 150 bp and 1 kb without affecting the estimation of selection, while for the reference alignment, we found that a length of 3 kb is adequate for proper testing. We also simulated sequence alignments under different classes of evolution and validated our ability to distinguish positive selection from relaxation of constraint and neutral evolution. Finally, we re-confirmed that a quarter of all non-coding Human Accelerated Elements are evolving by positive selection. CONCLUSION: Here, we introduce a method we called adaptiPhy, which adds significant improvements to our earlier method that tests for branch-specific directional selection in noncoding sequences. The motivation for these improvements is to provide a more sensitive and better targeted characterization of directional selection and neutral evolution across the genome.

Ocean acidification induces distinct transcriptomic responses across life history stages of the sea urchin Heliocidaris erythrogramma.

Ocean acidification (OA) from seawater uptake of rising carbon dioxide emissions impairs development in marine invertebrates, particularly in calcifying species. Plasticity in gene expression is thought to mediate many of these physiological effects, but how these responses change across life history stages remains unclear. The abbreviated lecithotrophic development of the sea urchin Heliocidaris erythrogramma provides a valuable opportunity to analyse gene expression responses across a wide range of life history stages, including the benthic, post-metamorphic juvenile. We measured the transcriptional response to OA in H. erythrogramma at three stages of the life cycle (embryo, larva, and juvenile) in a controlled breeding design. The results reveal a broad range of strikingly stage-specific impacts of OA on transcription, including changes in the number and identity of affected genes; the magnitude, sign, and variance of their expression response; and the developmental trajectory of expression. The impact of OA on transcription was notably modest in relation to gene expression changes during unperturbed development and much smaller than genetic contributions from parentage. The latter result suggests that natural populations may provide an extensive genetic reservoir of resilience to OA. Taken together, these results highlight the complexity of the molecular response to OA, its substantial life history stage specificity, and the importance of contextualizing the transcriptional response to pH stress in light of normal development and standing genetic variation to better understand the capacity for marine invertebrates to adapt to OA.

A comparative analysis of egg provisioning using mass spectrometry during rapid life history evolution in sea urchins.

A dramatic life history switch that has evolved numerous times in marine invertebrates is the transition from planktotrophic (feeding) to lecithotrophic (nonfeeding) larval development-an evolutionary tradeoff with many important developmental and ecological consequences. To attain a more comprehensive understanding of the molecular basis for this switch, we performed untargeted lipidomic and proteomic liquid chromatography-tandem mass spectrometry on eggs and larvae from three sea urchin species: the lecithotroph Heliocidaris erythrogramma, the closely related planktotroph Heliocidaris tuberculata, and the distantly related planktotroph Lytechinus variegatus. We identify numerous molecular-level changes possibly associated with the evolution of lecithotrophy in H. erythrogramma. We find the massive lipid stores of H. erythrogramma eggs are largely composed of low-density, diacylglycerol ether lipids that, contrary to expectations, appear to support postmetamorphic development and survivorship. Rapid premetamorphic development in this species may instead be powered by upregulated carbohydrate metabolism or triacylglycerol metabolism. We also find proteins involved in oxidative stress regulation are upregulated in H. erythrogramma eggs, and apoB-like lipid transfer proteins may be important for echinoid oogenic nutrient provisioning. These results demonstrate how mass spectrometry can enrich our understanding of life history evolution and organismal diversity by identifying specific molecules associated with distinct life history strategies and prompt new hypotheses about how and why these adaptations evolve.

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

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

Expression of genes and proteins of the pax-six-eya-dach network in the metamorphic sea urchin: Insights into development of the enigmatic echinoderm body plan and sensory structures.

BACKGROUND: Photoreception-associated genes of the Pax-Six-Eya-Dach network (PSEDN) are deployed for many roles in addition to photoreception development. In this first study of PSEDN genes during development of the pentameral body in sea urchins, we investigated their spatial expression in Heliocidaris erythrogramma. RESULTS: Expression of PSEDN genes in the hydrocoele of early (Dach, Eya, Six1/2) and/or late (Pax6, Six3/6) larvae, and the five hydrocoele lobes, the first morphological expression of pentamery, supports a role in body plan development. Pax6, Six1/2, and Six3/6 were localized to the primary and/or secondary podia and putative sensory/neuronal cells. Six1/2 and Six3/6 were expressed in the neuropil region in the terminal disc of the podia. Dach was localized to spines. Sequential up-regulation of gene expression as new podia and spines formed was evident. Rhabdomeric opsin and pax6 protein were localized to cells in the primary podia and spines. CONCLUSIONS: Our results support roles for PSEDN genes in development of the pentameral body plan, contributing to our understanding of how the most unusual body plan in the Bilateria may have evolved. Development of sensory cells within the Pax-Six expression field is consistent with the role of these genes in sensory cell development in diverse species. Developmental Dynamics 247:239-249, 2018. © 2017 Wiley Periodicals, Inc.

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

BACKGROUND: The molecular mechanisms underlying the development of the unusual echinoderm pentameral body plan and their likeness to mechanisms underlying the development of the bilateral plans of other deuterostomes are of interest in tracing body plan evolution. In this first study of the spatial expression of genes associated with Nodal and BMP2/4 signalling during the transition to pentamery in sea urchins, we investigate Heliocidaris erythrogramma, a species that provides access to the developing adult rudiment within days of fertilization. RESULTS: BMP2/4, and the putative downstream genes, Six1/2, Eya, Tbx2/3 and Msx were expressed in the earliest morphological manifestation of pentamery during development, the five hydrocoele lobes. The formation of the vestibular ectoderm, the specialized region overlying the left coelom that forms adult ectoderm, involved the expression of putative Nodal target genes Chordin, Gsc and BMP2/4 and putative BMP2/4 target genes Dlx, Msx and Tbx. The expression of Nodal, Lefty and Pitx2 in the right ectoderm, and Pitx2 in the right coelom, was as previously observed in other sea urchins. CONCLUSION: That genes associated with Nodal and BMP2/4 signalling are expressed in the hydrocoele lobes, indicates that they have a role in the developmental transition to pentamery, contributing to our understanding of how the most unusual body plan in the Bilateria may have evolved. We suggest that the Nodal and BMP2/4 signalling cascades might have been duplicated or split during the evolution to pentamery.

Gene expression and adaptive noncoding changes during human evolution.

BACKGROUND: Despite evidence for adaptive changes in both gene expression and non-protein-coding, putatively regulatory regions of the genome during human evolution, the relationship between gene expression and adaptive changes in cis-regulatory regions remains unclear. RESULTS: Here we present new measurements of gene expression in five tissues of humans and chimpanzees, and use them to assess this relationship. We then compare our results with previous studies of adaptive noncoding changes, analyzing correlations at the level of gene ontology groups, in order to gain statistical power to detect correlations. CONCLUSIONS: Consistent with previous studies, we find little correlation between gene expression and adaptive noncoding changes at the level of individual genes; however, we do find significant correlations at the level of biological function ontology groups. The types of function include processes regulated by specific transcription factors, responses to genetic or chemical perturbations, and differentiation of cell types within the immune system. Among functional categories co-enriched with both differential expression and noncoding adaptation, prominent themes include cancer, particularly epithelial cancers, and neural development and function.

The phylogeny of extant starfish (Asteroidea: Echinodermata) including Xyloplax, based on comparative transcriptomics.

Multi-locus phylogenetic studies of echinoderms based on Sanger and RNA-seq technologies and the fossil record have provided evidence for the Asterozoa-Echinozoa hypothesis. This hypothesis posits a sister relationship between asterozoan classes (Asteroidea and Ophiuroidea) and a similar relationship between echinozoan classes (Echinoidea and Holothuroidea). Despite this consensus around Asterozoa-Echinozoa, phylogenetic relationships within the class Asteroidea (sea stars or starfish) have been controversial for over a century. Open questions include relationships within asteroids and the status of the enigmatic taxon Xyloplax. Xyloplax is thought by some to represent a newly discovered sixth class of echinoderms - and by others to be an asteroid. To address these questions, we applied a novel workflow to a large RNA-seq dataset that encompassed a broad taxonomic and genomic sample. This study included 15 species sampled from all extant orders and 13 families, plus four ophiuroid species as an outgroup. To expand the taxonomic coverage, the study also incorporated five previously published transcriptomes and one previously published expressed sequence tags (EST) dataset. We developed and applied methods that used a range of alignment parameters with increasing permissiveness in terms of gap characters present within an alignment. This procedure facilitated the selection of phylogenomic data subsets from large amounts of transcriptome data. The results included 19 nested data subsets that ranged from 37 to 4,281loci. Tree searches on all data subsets reconstructed Xyloplax as a velatid asteroid rather than a new class. This result implies that asteroid morphology remains labile well beyond the establishment of the body plan of the group. In the phylogenetic tree with the highest average asteroid nodal support several monophyletic groups were recovered. In this tree, Forcipulatida and Velatida are monophyletic and form a clade that includes Brisingida as sister to Forcipulatida. Xyloplax is consistently recovered as sister to Pteraster. Paxillosida and Spinulosida are each monophyletic, with Notomyotida as sister to the Paxillosida. Valvatida is recovered as paraphyletic. The results from other data subsets are largely consistent with these results. Our results support the hypothesis that the earliest divergence event among extant asteroids separated Velatida and Forcipulatacea from Valvatacea and Spinulosida.

EchinoDB, an application for comparative transcriptomics of deeply-sampled clades of echinoderms.

BACKGROUND: One of our goals for the echinoderm tree of life project (http://echinotol.org) is to identify orthologs suitable for phylogenetic analysis from next-generation transcriptome data. The current dataset is the largest assembled for echinoderm phylogeny and transcriptomics. We used RNA-Seq to profile adult tissues from 42 echinoderm specimens from 24 orders and 37 families. In order to achieve sampling members of clades that span key evolutionary divergence, many of our exemplars were collected from deep and polar seas. DESCRIPTION: A small fraction of the transcriptome data we produced is being used for phylogenetic reconstruction. Thus to make a larger dataset available to researchers with a wide variety of interests, we made a web-based application, EchinoDB (http://echinodb.uncc.edu). EchinoDB is a repository of orthologous transcripts from echinoderms that is searchable via keywords and sequence similarity. CONCLUSIONS: From transcripts we identified 749,397 clusters of orthologous loci. We have developed the information technology to manage and search the loci their annotations with respect to the Sea Urchin (Strongylocentrotus purpuratus) genome. Several users have already taken advantage of these data for spin-off projects in developmental biology, gene family studies, and neuroscience. We hope others will search EchinoDB to discover datasets relevant to a variety of additional questions in comparative biology.

The impact of gene expression variation on the robustness and evolvability of a developmental gene regulatory network.

Regulatory interactions buffer development against genetic and environmental perturbations, but adaptation requires phenotypes to change. We investigated the relationship between robustness and evolvability within the gene regulatory network underlying development of the larval skeleton in the sea urchin Strongylocentrotus purpuratus. We find extensive variation in gene expression in this network throughout development in a natural population, some of which has a heritable genetic basis. Switch-like regulatory interactions predominate during early development, buffer expression variation, and may promote the accumulation of cryptic genetic variation affecting early stages. Regulatory interactions during later development are typically more sensitive (linear), allowing variation in expression to affect downstream target genes. Variation in skeletal morphology is associated primarily with expression variation of a few, primarily structural, genes at terminal positions within the network. These results indicate that the position and properties of gene interactions within a network can have important evolutionary consequences independent of their immediate regulatory role.