Parallel embryonic transcriptional programs evolve under distinct constraints and may enable morphological conservation amidst adaptation

Single-cell phylotranscriptomics of developmental and cell type evolution

Single-cell phylotranscriptomics is an emerging tool to reveal the molecular and cellular mechanisms of evolution. We summarize its utility in studying the hourglass pattern of ontogenetic evolution and for understanding the evolutionary history of cell types. The developmental hourglass model suggests that the mid-embryonic stage is the most conserved period of development across species, which is supported by morphological and molecular studies. Single-cell phylotranscriptomic analysis has revealed previously underappreciated heterogeneity in transcriptome ages among lineages and cell types throughout development, and has identified the lineages and tissues that drive the whole-organism hourglass pattern. Single-cell transcriptome age analyses also provide important insights into the origin of germ layers, the different selective forces on tissues during adaptation, and the evolutionary relationships between cell types.

Developmental hourglass: Verification by numerical evolution and elucidation by dynamical-systems theory

Determining the general laws between evolution and development is a fundamental biological challenge. Developmental hourglasses have attracted increased attention as candidates for such laws, but the necessity of their emergence remains elusive. We conducted evolutionary simulations of developmental processes to confirm the emergence of the developmental hourglass and unveiled its establishment. We considered organisms consisting of cells containing identical gene networks that control morphogenesis and evolved them under selection pressure to induce more cell types. By computing the similarity between the spatial patterns of gene expression of two species that evolved from a common ancestor, a developmental hourglass was observed, that is, there was a correlation peak in the intermediate stage of development. The fraction of pleiotropic genes increased, whereas the variance in individuals decreased, consistent with previous experimental reports. Reduction of the unavoidable variance by initial or developmental noise, essential for survival, was achieved up to the hourglass bottleneck stage, followed by diversification in developmental processes, whose timing is controlled by the slow expression dynamics conserved among organisms sharing the hourglass. This study suggests why developmental hourglasses are observed within a certain phylogenetic range of species.

Integrating Complex Life Cycles in Comparative Developmental Biology

The goal of comparative developmental biology is identifying mechanistic differences in embryonic development between different taxa and how these evolutionary changes have led to morphological and organizational differences in adult body plans. Much of this work has focused on direct-developing species in which the adult forms straight from the embryo and embryonic modifications have direct effects on the adult. However, most animal lineages are defined by indirect development, in which the embryo gives rise to a larval body plan and the adult forms by transformation of the larva. Historically, much of our understanding of complex life cycles is viewed through the lenses of ecology and zoology. In this review, we discuss the importance of establishing developmental rather than morphological or ecological criteria for defining developmental mode and explicitly considering the evolutionary implications of incorporating complex life cycles into broad developmental comparisons of embryos across metazoans.

Widespread priming of transcriptional regulatory elements by incipient accessibility or RNA polymerase II pause in early embryos of the sea urchin Strongylocentrotus purpuratus

Transcriptional regulatory elements (TREs) are the primary nodes that control developmental gene regulatory networks. In embryo stages, larvae, and adult differentiated red spherule cells of the sea urchin Strongylocentrotus purpuratus, transcriptionally engaged TREs are detected by Precision Run-On Sequencing (PRO-seq), which maps genome-wide at base pair resolution the location of paused or elongating RNA polymerase II (Pol II). In parallel, TRE accessibility is estimated by the Assay for Transposase-Accessible Chromatin using Sequencing (ATAC-seq). Our analysis identifies surprisingly early and widespread TRE accessibility in 4-cell cleavage embryos that is not necessarily followed by concurrent or subsequent transcription. TRE transcriptional differences identified by PRO-seq provide more contrast among embryonic stages than ATAC-seq accessibility differences, in agreement with the apparent excess of accessible but inactive TREs during embryogenesis. Global TRE accessibility reaches a maximum around the 20-hour late blastula stage, which coincides with the consolidation of major embryo regionalizations and peak histone variant H2A.Z expression. A transcriptional potency model based on labile nucleosome TRE occupancy driven by DNA sequences and the prevalence of histone variants is proposed in order to explain the basal accessibility of transcriptionally inactive TREs during embryogenesis. However, our results would not reconcile well with labile nucleosome models based on simple A/T sequence enrichment. In addition, a large number of distal TREs become transcriptionally disengaged during developmental progression, in support of an early Pol II paused model for developmental gene regulation that eventually resolves in transcriptional activation or silencing. Thus, developmental potency in early embryos may be facilitated by incipient accessibility and transcriptional pause at TREs.

Nematode gene annotation by machine-learning-assisted proteotranscriptomics enables proteome-wide evolutionary analysis

Nematodes encompass more than 24,000 described species, which were discovered in almost every ecological habitat, and make up >80% of metazoan taxonomic diversity in soils. The last common ancestor of nematodes is believed to date back to ∼650-750 million years, generating a large and phylogenetically diverse group to be explored. However, for most species high-quality gene annotations are incomprehensive or missing. Combining short-read RNA sequencing with mass spectrometry-based proteomics and machine-learning quality control in an approach called proteotranscriptomics, we improve gene annotations for nine genome-sequenced nematode species and provide new gene annotations for three additional species without genome assemblies. Emphasizing the sensitivity of our methodology, we provide evidence for two hitherto undescribed genes in the model organism Caenorhabditis elegans Extensive phylogenetic systems analysis using this comprehensive proteome annotation provides new insights into evolutionary processes of this metazoan group.

Distinct regulatory states control the elongation of individual skeletal rods in the sea urchin embryo

Background

Understanding how gene regulatory networks (GRNs) control developmental progression is a key to the mechanistic understanding of morphogenesis. The sea urchin larval skeletogenesis provides an excellent platform to tackle this question. In the early stages of sea urchin skeletogenesis, skeletogenic genes are uniformly expressed in the skeletogenic lineage. Yet, during skeletal elongation, skeletogenic genes are expressed in distinct spatial sub‐domains. The regulation of differential gene expression during late skeletogenesis is not well understood.

Results

Here we reveal the dynamic expression of the skeletogenic regulatory genes that define a specific regulatory state for each pair of skeletal rods, in the sea urchin Paracentrotus lividus. The vascular endothelial growth factor (VEGF) signaling, essential for skeleton formation, specifically controls the migration of cells that form the postoral and distal anterolateral skeletogenic rods. VEGF signaling also controls the expression of regulatory genes in cells at the tips of the postoral rods, including the transcription factors Pitx1 and MyoD1. Pitx1 activity is required for normal skeletal elongation and for the expression of some of VEGF target genes.

Conclusions

Our study illuminates the fine‐tuning of the regulatory system during the transition from early to late skeletogenesis that gives rise to rod‐specific regulatory states.

The skeletogenic transcription factors form specific regulatory states in various skeletogenic sub‐populations.

Late VEGF signaling controls the regulatory states at the tips of the post‐oral and anterolateral skeletal rods.

VEGF signaling controls the expression of the transcription factors, MyoD1 and Pitx1.

Pitx1 activity is required for normal skeletal elongation and for the expression of some of VEGF target genes.

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.

Identification and prediction of developmental enhancers in sea urchin embryos

Background

The transcription of developmental regulatory genes is often controlled by multiple cis-regulatory elements. The identification and functional characterization of distal regulatory elements remains challenging, even in tractable model organisms like sea urchins.

Results

We evaluate the use of chromatin accessibility, transcription and RNA Polymerase II for their ability to predict enhancer activity of genomic regions in sea urchin embryos. ATAC-seq, PRO-seq, and Pol II ChIP-seq from early and late blastula embryos are manually contrasted with experimental cis-regulatory analyses available in sea urchin embryos, with particular attention to common developmental regulatory elements known to have enhancer and silencer functions differentially deployed among embryonic territories. Using the three functional genomic data types, machine learning models are trained and tested to classify and quantitatively predict the enhancer activity of several hundred genomic regions previously validated with reporter constructs in vivo.

Conclusions

Overall, chromatin accessibility and transcription have substantial power for predicting enhancer activity. For promoter-overlapping cis-regulatory elements in particular, the distribution of Pol II is the best predictor of enhancer activity in blastula embryos. Furthermore, ATAC- and PRO-seq predictive value is stage dependent for the promoter-overlapping subset. This suggests that the sequence of regulatory mechanisms leading to transcriptional activation have distinct relevance at different levels of the developmental gene regulatory hierarchy deployed during embryogenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07936-0.

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.

Something old, something new, something borrowed, something red: the origin of ecologically relevant novelties in Hemiptera

Comparative transcriptomics, applied in an evolutionary context, has transformed the possibilities for studying phenotypic evolution in non-model taxa. We review recent discoveries about the development of novel, ecologically relevant phenotypes in hemipteran insects. These discoveries highlight the diverse genomic substrates of novelty: 'something old', when novelty results from changes in the regulation of existing genes or gene duplication; 'something new', wherein lineage-restricted genes contribute to the evolution of new phenotypes; and 'something borrowed', showcasing contributions of horizontal gene transfer to the evolution of novelty, including carotenoid synthesis (resulting in 'something red'). These findings show the power and flexibility of comparative transcriptomic approaches for expanding beyond the 'toolkit' model for the evolution of development. We conclude by raising questions about the relationship between new genes and new traits and outlining a research framework for answering them in Hemiptera.

The developmental hourglass model and recapitulation: An attempt to integrate the two models

Recapitulation is a hypothetical concept that assumes embryogenesis of an animal parallels its own phylogenetic history, sequentially developing from more ancestral features to more derived ones. This concept predicts that the earliest developmental stage of various animals should represent the most evolutionarily conserved patterns. Recent transcriptome‐based studies, on the other hand, have reported that mid‐embryonic, organogenetic periods show the highest level of conservation (the developmental hourglass model). This, however, does not rule out the possibility that recapitulation would still be detected after the mid‐embryonic period. In accordance with this, recapitulation‐like morphological features are enriched in late developmental stages. Moreover, our recent chromatin accessibility‐based study provided molecular evidence for recapitulation in the mid‐to‐late embryogenesis of vertebrates, as newly evolved gene regulatory elements tended to be activated at late embryonic stages. In this review, we revisit the recapitulation hypothesis, together with recent molecular‐based studies that support the developmental hourglass model. We contend that the recapitulation hypothesis does not entirely contradict the developmental hourglass model and that these two may even coexist in later embryonic stages of vertebrates. Finally, we review possible mechanisms underlying the recapitulation pattern of chromatin accessibility together with the hourglass‐like evolutionary conservation in vertebrate embryogenesis.

Recapitulation pattern has been reported for chromatin accessibility during the mid‐to‐late embryogenesis.

The observed recapitulation pattern and the developmental hourglass model may coexist.

The possible evolutionary mechanisms underlying tendencies of embryonic evolution were discussed.

Genomic insights of body plan transitions from bilateral to pentameral symmetry in Echinoderms

Echinoderms are an exceptional group of bilaterians that develop pentameral adult symmetry from a bilaterally symmetric larva. However, the genetic basis in evolution and development of this unique transformation remains to be clarified. Here we report newly sequenced genomes, developmental transcriptomes, and proteomes of diverse echinoderms including the green sea urchin (L. variegatus), a sea cucumber (A. japonicus), and with particular emphasis on a sister group of the earliest-diverged echinoderms, the feather star (A. japonica). We learned that the last common ancestor of echinoderms retained a well-organized Hox cluster reminiscent of the hemichordate, and had gene sets involved in endoskeleton development. Further, unlike in other animal groups, the most conserved developmental stages were not at the body plan establishing phase, and genes normally involved in bilaterality appear to function in pentameric axis development. These results enhance our understanding of the divergence of protostomes and deuterostomes almost 500 Mya.

Li et al. investigate the evolution and genetic basis of the adult pentameral body plan in echinoderms using genomic, transcriptomic, and proteomic data. They determine that the last common ancestor of echinoderms contained an organized Hox cluster and endoskeleton genes, and suggest that cooption of bilateral development genes was involved in evolution of the pentameric body plan.

Developmental transcriptomes of the sea star, Patiria miniata, illuminate how gene expression changes with evolutionary distance

Understanding how changes in developmental gene expression alter morphogenesis is a fundamental problem in development and evolution. A promising approach to address this problem is to compare the developmental transcriptomes between related species. The echinoderm phylum consists of several model species that have significantly contributed to the understanding of gene regulation and evolution. Particularly, the regulatory networks of the sea star, Patiria miniata (P. miniata), have been extensively studied, however developmental transcriptomes for this species were lacking. Here we generated developmental transcriptomes of P. miniata and compared these with those of two sea urchins species. We demonstrate that the conservation of gene expression depends on gene function, cell type and evolutionary distance. With increasing evolutionary distance the interspecies correlations in gene expression decreases. The reduction is more severe in the correlations between morphologically equivalent stages (diagonal elements) than in the correlation between morphologically distinct stages (off-diagonal elements). This could reflect a decrease in the morphological constraints compared to other constraints that shape gene expression at large evolutionary divergence. Within this trend, the interspecies correlations of developmental control genes maintain their diagonality at large evolutionary distance, and peak at the onset of gastrulation, supporting the hourglass model of phylotypic stage conservation.

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.

Comparative Studies of Gene Expression Kinetics: Methodologies and Insights on Development and Evolution

Across the animal kingdom, embryos of closely related species show high morphological similarity despite genetic and environmental distances. Deciphering the molecular mechanisms that underlie morphological conservation and those that support embryonic adaptation are keys to understand developmental robustness and evolution. Comparative studies of developmental gene regulatory networks can track the genetic changes that lead to evolutionary novelties. However, these studies are limited to a relatively small set of genes and demand extensive experimental efforts. An alternative approach enabled by next-generation sequencing, is to compare the expression kinetic of large sets of genes between different species. The advantages of these comparisons are that they can be done relatively easily, for any species and they provide information of all expressed genes. The challenge in these experiments is to compare the kinetic profiles of thousands of genes between species that develop in different rates. Here we review recent comparative studies that tackled the challenges of accurate staging and large-scale analyses using different computational approaches. These studies reveal how correct temporal scaling exposes the striking conservation of developmental gene expression between morphologically similar species. Different clustering approaches are used to address various comparative questions and identify the conservation and divergence of large gene sets. We discuss the unexpected contribution of housekeeping genes to the interspecies correlations and how this contribution distorts the hourglass pattern generated by developmental genes. Overall, we demonstrate how comparative studies of gene expression kinetics can provide novel insights into the developmental constraints and plasticity that shape animal body plans.