Genome evolution and divergence in cis-regulatory architecture is associated with condition-responsive development in horned dung beetles.

Phenotypic plasticity is thought to be an important driver of diversification and adaptation to environmental variation, yet the genomic mechanisms mediating plastic trait development and evolution remain poorly understood. The Scarabaeinae, or true dung beetles, are a species-rich clade of insects recognized for their highly diversified nutrition-responsive development including that of cephalic horns-evolutionarily novel, secondary sexual weapons that exhibit remarkable intra- and interspecific variation. Here, we investigate the evolutionary basis for horns as well as other key dung beetle traits via comparative genomic and developmental assays. We begin by presenting chromosome-level genome assemblies of three dung beetle species in the tribe Onthophagini (> 2500 extant species) including Onthophagus taurus, O. sagittarius, and Digitonthophagus gazella. Comparing these assemblies to those of seven other species across the order Coleoptera identifies evolutionary changes in coding sequence associated with metabolic regulation of plasticity and metamorphosis. We then contrast chromatin accessibility in developing head horn tissues of high- and low-nutrition O. taurus males and females and identify distinct cis-regulatory architectures underlying nutrition- compared to sex-responsive development, including a large proportion of recently evolved regulatory elements sensitive to horn morph determination. Binding motifs of known and new candidate transcription factors are enriched in these nutrition-responsive open chromatin regions. Our work highlights the importance of chromatin state regulation in mediating the development and evolution of plastic traits, demonstrates gene networks are highly evolvable transducers of environmental and genetic signals, and provides new reference-quality genomes for three species that will bolster future developmental, ecological, and evolutionary studies of this insect group.

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.

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.

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.

The maternal-zygotic transition and zygotic activation of the Mnemiopsis leidyi genome occurs within the first three cleavage cycles.

The maternal-zygotic transition (MZT) describes the developmental reprogramming of gene expression marked by the degradation of maternally supplied gene products and activation of the zygotic genome. While the timing and duration of the MZT vary among taxa, little is known about early-stage transcriptional dynamics in the non-bilaterian phylum Ctenophora. We sought to better understand the extent of maternal mRNA loading and subsequent differential transcript abundance during the earliest stages of development by performing comprehensive RNA-sequencing-based analyses of mRNA abundance in single- and eight-cell stage embryos in the lobate ctenophore Mnemiopsis leidyi. We found 1,908 contigs with significant differential abundance between single- and eight-cell stages, of which 1,208 contigs were more abundant at the single-cell stage and 700 contigs were more abundant at the eight-cell stage. Of the differentially abundant contigs, 267 were exclusively present in the eight-cell samples, providing strong evidence that both the MZT and zygotic genome activation (ZGA) have commenced by the eight-cell stage. Many highly abundant transcripts encode genes involved in molecular mechanisms critical to the MZT, such as maternal transcript degradation, serine/threonine kinase activity, and chromatin remodeling. Our results suggest that chromosomal restructuring, which is critical to ZGA and the initiation of transcriptional regulation necessary for normal development, begins by the third cleavage within 1.5 hr post-fertilization in M. leidyi.

Gene regulatory networks underlying the development and evolution of plasticity in horned beetles.

Horned beetles have emerged as a powerful study system with which to investigate the developmental mechanisms underlying environment-responsive development and its evolution. We begin by reviewing key advances in our understanding of the diverse roles played by transcription factors, endocrine regulators, and signal transduction pathways in the regulation of horned beetle plasticity. We then explore recent efforts aimed at understanding how such condition-specific expression may be regulated in the first place, as well as how the differential expression of master regulators may instruct conditional expression of downstream target genes. Here, we focus on the significance of chromatin remodeling as a powerful but thus far understudied mechanism able to facilitate trait-, sex-, and species-specific responses to environmental conditions.

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 patterns across development and evolution; however, little is known about the genetic and molecular mechanisms that influence chromatin configuration itself. Because cis and trans influences can both theoretically influence the accessibility of the epigenome, we sought to better characterize the role that both of these mechanisms play in altering chromatin accessibility in two closely related sea urchin species. Using hybrids of the two species, 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 three important developmental stages, and compared our results to similar patterns in the transcriptome. We found extensive cis- and trans-based influences on evolutionary changes in chromatin, with cis effects slightly more numerous and larger in effect. Genetic mechanisms influencing gene expression and chromatin configuration are correlated, but differ in several important ways. Maternal influences also appear to have more of an effect on chromatin accessibility than on gene expression, persisting well past the maternal-to-zygotic transition. Furthermore, chromatin accessibility near GRN genes appears to be regulated differently than the rest of the epigenome, and indicates that trans factors may 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 the regulation of chromatin is more genetically complex than was previously appreciated.

Near-Chromosomal-Level Genome Assembly of the Sea Urchin Echinometra lucunter, a Model for Speciation in the Sea.

Echinometra lucunter, the rock-boring sea urchin, is a widely distributed echinoid and a model for ecological studies of reproduction, responses to climate change, and speciation. We present a near chromosome-level genome assembly of E. lucunter, including 21 scaffolds larger than 10 Mb predicted to represent each of the chromosomes of the species. The 760.4 Mb assembly includes a scaffold N50 of 30.0 Mb and BUSCO (benchmarking universal single-copy orthologue) single copy and a duplicated score of 95.8% and 1.4%, respectively. Ab-initio gene model prediction and annotation with transcriptomic data constructed 33,989 gene models composing 50.4% of the assembly, including 37,036 transcripts. Repetitive elements make up approximately 39.6% of the assembly, and unresolved gap sequences are estimated to be 0.65%. Whole genome alignment with Echinometra sp. EZ revealed high synteny and conservation between the two species, further bolstering Echinometra as an emerging genus for comparative genomics studies. This genome assembly represents a high-quality genomic resource for future evolutionary and developmental studies of this species and more broadly of echinoderms.

A Chromosome-level Genome Assembly of the Highly Heterozygous Sea Urchin Echinometra sp. EZ Reveals Adaptation in the Regulatory Regions of Stress Response Genes.

Echinometra is the most widespread genus of sea urchin and has been the focus of a wide range of studies in ecology, speciation, and reproduction. However, available genetic data for this genus are generally limited to a few select loci. Here, we present a chromosome-level genome assembly based on 10x Genomics, PacBio, and Hi-C sequencing for Echinometra sp. EZ from the Persian/Arabian Gulf. The genome is assembled into 210 scaffolds totaling 817.8 Mb with an N50 of 39.5 Mb. From this assembly, we determined that the E. sp. EZ genome consists of 2n = 42 chromosomes. BUSCO analysis showed that 95.3% of BUSCO genes were complete. Ab initio and transcript-informed gene modeling and annotation identified 29,405 genes, including a conserved Hox cluster. E. sp. EZ can be found in high-temperature and high-salinity environments, and we therefore compared E. sp. EZ gene families and transcription factors associated with environmental stress response ("defensome") with other echinoid species with similar high-quality genomic resources. While the number of defensome genes was broadly similar for all species, we identified strong signatures of positive selection in E. sp. EZ noncoding elements near genes involved in environmental response pathways as well as losses of transcription factors important for environmental response. These data provide key insights into the biology of E. sp. EZ as well as the diversification of Echinometra more widely and will serve as a useful tool for the community to explore questions in this taxonomic group and beyond.

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.

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.

Transcriptomic analysis of sea star development through metamorphosis to the highly derived pentameral body plan with a focus on neural transcription factors.

The Echinodermata is characterized by a secondarily evolved pentameral body plan. While the evolutionary origin of this body plan has been the subject of debate, the molecular mechanisms underlying its development are poorly understood. We assembled a de novo developmental transcriptome from the embryo through metamorphosis in the sea star Parvulastra exigua. We use the asteroid model as it represents the basal-type echinoderm body architecture. Global variation in gene expression distinguished the gastrula profile and showed that metamorphic and juvenile stages were more similar to each other than to the pre-metamorphic stages, pointing to the marked changes that occur during metamorphosis. Differential expression and gene ontology (GO) analyses revealed dynamic changes in gene expression throughout development and the transition to pentamery. Many GO terms enriched during late metamorphosis were related to neurogenesis and signalling. Neural transcription factor genes exhibited clusters with distinct expression patterns. A suite of these genes was up-regulated during metamorphosis (e.g. Pax6, Eya, Hey, NeuroD, FoxD, Mbx, and Otp). In situ hybridization showed expression of neural genes in the CNS and sensory structures. Our results provide a foundation to understand the metamorphic transition in echinoderms and the genes involved in development and evolution of pentamery.

Chromosomal-Level Genome Assembly of the Sea Urchin Lytechinus variegatus Substantially Improves Functional Genomic Analyses.

Lytechinus variegatus is a camarodont sea urchin found widely throughout the western Atlantic Ocean in a variety of shallow-water marine habitats. Its distribution, abundance, and amenability to developmental perturbation make it a popular model for ecologists and developmental biologists. Here, we present a chromosomal-level genome assembly of L. variegatus generated from a combination of PacBio long reads, 10× Genomics sequencing, and HiC chromatin interaction sequencing. We show L. variegatus has 19 chromosomes with an assembly size of 870.4 Mb. The contiguity and completeness of this assembly are reflected by a scaffold length N50 of 45.5 Mb and BUSCO completeness score of 95.5%. Ab initio and transcript-informed gene modeling and annotation identified 27,232 genes with an average gene length of 12.6 kb, comprising an estimated 39.5% of the genome. Repetitive regions, on the other hand, make up 45.4% of the genome. Physical mapping of well-studied developmental genes onto each chromosome reveals nonrandom spatial distribution of distinct genes and gene families, which provides insight into how certain gene families may have evolved and are transcriptionally regulated in this species. Lastly, aligning RNA-seq and ATAC-seq data onto this assembly demonstrates the value of highly contiguous, complete genome assemblies for functional genomics analyses that is unattainable with fragmented, incomplete assemblies. This genome will be of great value to the scientific community as a resource for genome evolution, developmental, and ecological studies of this species and the Echinodermata.