Abstract: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.

Conserved and specific gene expression patterns in the embryonic development of tardigrades

The Emergence of the Tardigrade Hypsibius exemplaris as a Model System

Weighted gene co-expression network analysis of embryos and first instar larvae of the horseshoe crab Tachypleus tridentatus uncovers development gene networks

Expression patterns of FGF and BMP pathway genes in the tardigrade

Comparative analyses of the transcriptome among three development stages of Zeugodacus tau larvae (Diptera: Tephritidae)

Weighted Gene Co-Expression Network Analysis Reveals Potential Genes Involved in Early Metamorphosis Process in Sea Cucumber Apostichopus Japonicus.

RNA-seq analysis reveals changes in mRNA expression during development in Daphnia mitsukuri

Ecdysteroid-dependent molting in tardigrades

Recent Reconfiguration of an Ancient Developmental Gene Regulatory Network in Heliocidaris Sea Urchins

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

A comparative analysis of transcription factor expression during metazoan embryonic development

Temporal dynamics of gene expression during metamorphosis in two distant Drosophila species

The hourglass and the early conservation models - co-existing evolutionary patterns in vertebrate development

Expression of distal limb patterning genes in Hypsibius exemplaris indicate regionalization and suggest distal identity of tardigrade legs

Gene regulatory dynamics during the development of a paleopteran insect, the mayfly Cloeon dipterum

Evolution of maternal and early zygotic transcript regulation across Drosophila

Deep Transcriptomic Analysis Reveals the Dynamic Developmental Progression During Early Development of Channel Catfish (ictalurus Punctatus)

Comparative transcriptomic evidence of physiological changes and potential relationships in vertebrates under different dormancy states

Transcriptome Sequencing Analysis Reveals Dynamic Changes in Major Biological Functions during the Early Development of Clearhead Icefish, Protosalanx chinensis

Characterization And Functional Analysis Of Hsp18.3 Gene In The Red Flour Beetle, Tribolium Castaneum

Single-animal, single-tube RNA extraction for comparison of relative transcript levels via qRT-PCR in the tardigrade Hypsibius exemplaris.