Recent reconfiguration of an ancient developmental gene regulatory network in Heliocidaris sea urchins

Changes in developmental gene regulatory networks (dGRNs) underlie much of the diversity of life, but the evolutionary mechanisms that operate on regulatory interactions remain poorly understood. Closely related species with extreme phenotypic divergence provide a valuable window into the genetic and molecular basis for changes in dGRNs and their relationship to adaptive changes in organismal traits. Here we analyse genomes, epigenomes and transcriptomes during early development in two Heliocidaris sea urchin species that exhibit highly divergent life histories and in an outgroup species. Positive selection and chromatin accessibility modifications within putative regulatory elements are enriched on the branch leading to the derived life history, particularly near dGRN genes. Single-cell transcriptomes reveal a dramatic delay in cell fate specification in the derived state, which also has far fewer open chromatin regions, especially near conserved cell fate specification genes. Experimentally perturbing key transcription factors reveals profound evolutionary changes to early embryonic patterning events, disrupting regulatory interactions previously conserved for ~225 million years. These results demonstrate that natural selection can rapidly reshape developmental gene expression on a broad scale when selective regimes abruptly change. More broadly, even highly conserved dGRNs and patterning mechanisms in the early embryo remain evolvable under appropriate ecological circumstances.

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.

Extreme phenotypic divergence and the evolution of development

As analyses of developmental mechanisms extend to ever more species, it becomes important to understand not just what is conserved or altered during evolution, but why. Closely related species that exhibit extreme phenotypic divergence can be uniquely informative in this regard. A case in point is the sea urchin genus Heliocidaris, which contains species that recently evolved a life history involving nonfeeding larvae following nearly half a billion years of prior evolution with feeding larvae. The resulting shift in selective regimes produced rapid and surprisingly extensive changes in developmental mechanisms that are otherwise highly conserved among echinoderm species. The magnitude and extent of these changes challenges the notion that conservation of early development in echinoderms is largely due to internal constraints that prohibit modification and instead suggests that natural selection actively maintains stability of inherently malleable trait developmental mechanisms over immense time periods. Knowing how and why natural selection changed during the evolution of nonfeeding larvae can also reveal why developmental mechanisms do and do not change in particular ways.

Dorsal-ventral axis formation in sea urchin embryos

Most sea urchin species produce planktonic feeding larvae with distinct dorsal-ventral polarity. Such morphological indicators of polarity arise after gastrulation, when several morphogenesis and cell differentiation events occur differentially along the dorsal-ventral axis. For instance, the gut bends toward the ventral side where the mouth will form, skeletogenesis occurs initially near the ventral side with the forming skeleton extending dorsally, and pigment cells differentiate and embed in the dorsal ectoderm. The patterning mechanisms and gene regulatory networks underlying these events have been extensively studied. Two opposing TGF-β signaling pathways, Nodal and BMP, play key roles in all three germ layers to respectively pattern the sea urchin ventral and dorsal sides. In this chapter, I describe our current understanding of sea urchin dorsal-ventral patterning mechanisms. Additionally, differences in the patterning mechanisms observed in lecithotrophic sea urchins (nonfeeding larvae) and in cidaroid sea urchins are also discussed, along with evolutionary insights gained from comparative analyses.

Transcriptomic analysis of Nodal - and BMP- associated genes during development to the juvenile seastar in Parvulastra exigua (Asterinidae)

The molecular mechanisms underlying development of the pentameral body of adult echinoderms are poorly understood but are important to solve with respect to evolution of a unique body plan that contrasts with the bilateral body plan of other deuterostomes. As Nodal and BMP2/4 signalling is involved in axis formation in larvae and development of the echinoderm body plan, we used the developmental transcriptome generated for the asterinid seastar Parvulastra exigua to investigate the temporal expression patterns of Nodal and BMP2/4 genes from the embryo and across metamorphosis to the juvenile. For echinoderms, the Asteroidea represents the basal-type body architecture with a distinct (separated) ray structure. Parvulastra exigua has lecithotrophic development forming the juvenile soon after gastrulation providing ready access to the developing adult stage. We identified 39 genes associated with the Nodal and BMP2/4 network in the P. exigua developmental transcriptome. Clustering analysis of these genes resulted in 6 clusters with similar temporal expression patterns across development. A co-expression analysis revealed genes that have similar expression profiles as Nodal and BMP2/4. These results indicated genes that may have a regulatory relationship in patterning morphogenesis of the juvenile seastar. Developmental RNA-seq analyses of Parvulastra exigua show changes in Nodal and BMP2/4 signalling genes across the metamorphic transition. We provide the foundation for detailed analyses of this cascade in the evolution of the unusual pentameral echinoderm body and its deuterostome affinities.

Genetic basis for divergence in developmental gene expression in two closely related sea urchins

The genetic basis for divergence in developmental gene expression among species is poorly understood, despite growing evidence that such changes underlie many interesting traits. Here we quantify transcription in hybrids of Heliocidaris tuberculata and Heliocidaris erythrogramma, two closely related sea urchins with highly divergent developmental gene expression and life histories. We find that most expression differences between species result from genetic influences that affect one stage of development, indicating limited pleiotropic consequences for most mutations that contribute to divergence in gene expression. Activation of zygotic transcription is broadly delayed in H. erythrogramma, the species with the derived life history, despite its overall faster premetamorphic development. Altered expression of several terminal differentiation genes associated with the derived larval morphology of H. erythrogramma is based largely on differences in the expression or function of their upstream regulators, providing insights into the genetic basis for the evolution of key life history traits.

The development and metamorphosis of the indirect developing acorn worm Schizocardium californicum (Enteropneusta: Spengelidae)

Background

Enteropneusts are benthic marine invertebrates that belong to the deuterostome phylum Hemichordata. The two main clades of enteropneusts are defined by differences in early life history strategies. In the Spengelidae and Ptychoderidae, development is indirect via a planktotrophic tornaria larva. In contrast, development in the Harrimanidae is direct without an intervening larval life history stage. Most molecular studies in the development and evolution of the enteropneust adult body plan have been carried out in the harrimanid Saccoglossus kowalevskii. In order to compare these two developmental strategies, we have selected the spengelid enteropneust Schizocardium californicum as a suitable indirect developing species for molecular developmental studies. Here we describe the methods for adult collecting, spawning and larval rearing in Schizocardium californicum, and describe embryogenesis, larval development, and metamorphosis, using light microscopy, immunocytochemistry and confocal microscopy.

Results

Adult reproductive individuals can be collected intertidally and almost year-round. Spawning can be triggered by heat shock and large numbers of larvae can be reared through metamorphosis under laboratory conditions. Gastrulation begins at 17 h post-fertilization (hpf) and embryos hatch at 26 hpf as ciliated gastrulae. At 3 days post-fertilization (dpf), the tornaria has a circumoral ciliary band, mouth, tripartite digestive tract, protocoel, larval muscles and a simple serotonergic nervous system. The telotroch develops at 5 dpf. In the course of 60 days, the serotonergic nervous system becomes more elaborate, the posterior coeloms develop, and the length of the circumoral ciliary band increases. At the end of the larval stage, larval muscles disappear, gill slits form, and adult muscles develop. Metamorphosis occurs spontaneously when the larva reaches its maximal size (ca. 3 mm), and involves loss and reorganization of larval structures (muscles, nervous system, digestive tract), as well as development of adult structures (adult muscles, tripartite body organization).

Conclusions

This study will enable future research in S. californicum to address long standing questions related to the evolution of axial patterning mechanisms, germ layer induction, neurogenesis and neural patterning, the mechanisms of metamorphosis, the relationships between larval and adult body plans, and the evolution of metazoan larval forms.

Effects of Nodal inhibition on development of temnopleurid sea urchins

Adult rudiment formation in some temnopleurids begins with the formation of a cell mass that is pinched off the left ectoderm in early larval development. The cell mass forms the adult rudiment with the left coelomic pouch of the mesodermal region. However, details of the mechanisms to establish position of the cell mass are still unknown. We analyzed the inhibiting effect of Nodal, a factor for morphogenesis of the oral region and right side, for location of the cell mass, in four temnopleurids. Pulse inhibition, at least 5 min inhibition, during coelomic pouch formation allowed a cell mass to form on both sides, whereas treatments after that period did not. These results indicate that Nodal signaling controls the oral-aboral axis before gastrulation and then affects the position of the cell mass and adult rudiment up to coelomic pouch formation. They also indicate that the position of the adult rudiment under Nodal signaling pathways is conserved in temnopleurids, as adult rudiment formation is dependent on the cell mass.

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.

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

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

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.

Nodal signaling is required for mesodermal and ventral but not for dorsal fates in the indirect developing hemichordate, Ptychodera flava

Nodal signaling plays crucial roles in vertebrate developmental processes such as endoderm and mesoderm formation, and axial patterning events along the anteroposterior, dorsoventral and left-right axes. In echinoderms, Nodal plays an essential role in the establishment of the dorsoventral axis and left-right asymmetry, but not in endoderm or mesoderm induction. In protostomes, Nodal signaling appears to be involved only in establishing left-right asymmetry. Hence, it is hypothesized that Nodal signaling has been co-opted to pattern the dorsoventral axis of deuterostomes and for endoderm, mesoderm formation as well as anteroposterior patterning in chordates. Hemichordata, together with echinoderms, represent the sister taxon to chordates. In this study, we analyze the role of Nodal signaling in the indirect developing hemichordate Ptychodera flava. In particular, we show that during gastrulation nodal transcripts are detected in a ring of cells at the vegetal pole that gives rise to endomesoderm and in the ventral ectoderm at later stages of development. Inhibition of Nodal function disrupts dorsoventral fates and also blocks formation of the larval mesoderm. Interestingly, molecular analysis reveals that only mesodermal, apical and ventral gene expression is affected while the dorsal side appears to be patterned correctly. Taken together, this study suggests that the co-option of Nodal signaling in mesoderm formation and potentially in anteroposterior patterning has occurred prior to the emergence of chordates and that Nodal signaling on the ventral side is uncoupled from BMP signaling on the dorsal side, representing a major difference from the molecular mechanisms of dorsoventral patterning events in echinoderms.

Transcriptomic analysis of the highly derived radial body plan of a sea urchin

With their complex life cycle and highly derived body plan, echinoderms are unique among bilaterians. Although early development has been intensively studied, the molecular mechanisms underlying development of the adult echinoderm and its unusual radial body plan are largely unknown. To investigate the evolution of developmental changes in gene expression underlying radial body plan development and metamorphosis, we assembled a reference transcriptome de novo and used RNA-seq to measure gene expression profiles across larval, metamorphic, and postmetamorphic life cycle phases in the sea urchin Heliocidaris erythrogramma. Our results present a high-resolution view of gene expression dynamics during the complex transition from pre- to postmetamorphic development and suggest that distinct sets of regulatory and effector proteins are used during different life history phases. These analyses provide an important foundation for more detailed analyses of the evolution of the radial adult body of echinoderms.

Oral-aboral identity displayed in the expression of HpHox3 and HpHox11/13 in the adult rudiment of the sea urchin Holopneustes purpurescens

Hox genes are noted for their roles in specifying axial identity in bilateral forms. In the radial echinoderms, the axis whose identity Hox genes might specify remains unclear. From the expression of Hox genes in the development of the sea urchin Holopneustes purpurescens reported here and that reported previously, we clarify the axis that might be specified by Hox genes in echinoderms. The expression of HpHox11/13 here is described at three developmental stages. The expression is around the rim of the blastopore in gastrulae, in the archenteron wall and adjacent mesoderm in early vestibula larvae, and in a patch of mesoderm close to the archenteron wall in later vestibula larvae. The retained expression of HpHox11/13 in the patch of mesoderm in the later vestibula larvae is, we suggest, indicative of a posterior or an aboral growth zone. The expression of HpHox3 at the echinoid-rudiment stage, in contrast, is in oral mesoderm beneath the epineural folds, concentrated in sites where the first three adult spines form. With the expression of HpHox5 and HpHox11/13 reported previously, the expressions here support the role of Hox genes in specifying oral-aboral identity in echinoderms. How such specification and a posterior growth zone add support to a concept of the structural homology between echinoderms and chordates is discussed.

Coelomogenesis during the abbreviated development of the echinoid Heliocidaris erythrogramma and the developmental origin of the echinoderm pentameral body plan

The development of the coeloms is described in an echinoid with an abbreviated larval development and shows the early morphogenesis of the coeloms of the adult stage. The development is described from images obtained by laser scanning confocal microscopy. The development in Heliocidaris erythrogramma is asymmetric with a larger left coelom forming on the larval-left side and a smaller right coelom forming on the larval-right side. The right coelom forms after the development of the left coelom is well advanced. The hydrocoele forms from the anterior part of the left coelom. The five lobes of the hydrocoele from which the pentamery of the adult derives take shape on the outer, distal wall of the anterior part of the left coelom. The hydrocoele separates from the more posterior part of the left coelom, which becomes the left posterior coelom. The lobes of the hydrocoele are named, based on the site of the connexion of the stone canal to the hydrocoele. The mouth is assumed to form by penetration through only the outer, distal wall of the hydrocoele and the ectoderm. Both larval and adult polarities are evident in this larva. A comparison with coelomogenesis in the asteroid Parvulastra exigua, which also has an abbreviated development, leads to predictions of homology between the echinoderm and chordate phyla that do not require the hypothesis of a dorsoventral inversion event in chordates.

Ventralization of an indirect developing hemichordate by NiCl₂ suggests a conserved mechanism of dorso-ventral (D/V) patterning in Ambulacraria (hemichordates and echinoderms)

One of the earliest steps in embryonic development is the establishment of the future body axes. Morphological and molecular data place the Ambulacraria (echinoderms and hemichordates) within the Deuterostomia and as the sister taxon to chordates. Extensive work over the last decades in echinoid (sea urchins) echinoderms has led to the characterization of gene regulatory networks underlying germ layer specification and axis formation during embryogenesis. However, with the exception of recent studies from a direct developing hemichordate (Saccoglossus kowalevskii), very little is known about the molecular mechanism underlying early hemichordate development. Unlike echinoids, indirect developing hemichordates retain the larval body axes and major larval tissues after metamorphosis into the adult worm. In order to gain insight into dorso-ventral (D/V) patterning, we used nickel chloride (NiCl₂), a potent ventralizing agent on echinoderm embryos, on the indirect developing enteropneust hemichordate, Ptychodera flava. Our present study shows that NiCl₂ disrupts the D/V axis and induces formation of a circumferential mouth when treated before the onset of gastrulation. Molecular analysis, using newly isolated tissue-specific markers, shows that the ventral ectoderm is expanded at expense of dorsal ectoderm in treated embryos, but has little effect on germ layer or anterior-posterior markers. The resulting ventralized phenotype, the effective dose, and the NiCl₂ sensitive response period of Ptychodera flava, is very similar to the effects of nickel on embryonic development described in larval echinoderms. These strong similarities allow one to speculate that a NiCl₂ sensitive pathway involved in dorso-ventral patterning may be shared between echinoderms, hemichordates and a putative ambulacrarian ancestor. Furthermore, nickel treatments ventralize the direct developing hemichordate, S. kowalevskii indicating that a common pathway patterns both larval and adult body plans of the ambulacrarian ancestor and provides insight in to the origin of the chordate body plan.

Characterization and expression of a sea star otx ortholog (Protxβ1/2) in the larva of Patiriella regularis

A transcript of otx from the sea star Patiriella regularis (Protxβ1/2) was characterized and its expression in early bipinnaria larvae was documented by whole mount in situ hybridization (WMISH). The nucleotide sequence exhibited 94% identity with Amotxβ1/2 from the closely related species Patiria miniata. Protxβ1/2 was expressed strongly in the developing archenteron in the future fore and mid-gut regions. This was followed by expression of Protxβ1/2 in the developing enterocoels, mesodermal derivatives. This suggests a role for Protx in endomesoderm development. In coelom development, Protxβ1/2 was first expressed in the left coelom. Subsequently expression was evident in the right coelom, but localization was never as strong as in the left coelom. This asymmetry in Protxβ1/2 expression in the coeloms was evident up to the stage when they started to extend posteriorly. These data indicate that Protxβ1/2 may have a role in coelom development, particularly in the left coelom, a definitive adult structure.

Chapter 7. Axis formation and the rapid evolutionary transformation of larval form

Marine invertebrate embryos and larvae are diverse and can evolve rapidly, providing a link between early developmental and evolutionary mechanisms. We here discuss the role of evolutionary changes in axis formation, which is a crucial part of the patterning of marine embryos and larvae. We focus on sea urchin embryos, where axial features are well defined and subject to active current investigation. The genetic control of processes of formation of the three axial systems, animal-vegetal, dorsal-ventral, and left-right, is becoming established for species that undergo development via the feeding pluteus larva. These species represent the primitive condition among living sea urchins. We compare their developmental processes to the highly modified development of a species that has evolved a nonfeeding larva. This derived form has accelerated some elements of axis formation, and eliminated or modified others. Three features of embryonic/larval evolution stand out (1) evolution of developmental features occurs rapidly over geological time; (2) upstream gene regulatory systems of axis formation are conserved, whereas downstream features evolve rapidly; and (3) heterochronies play an important role.