Genome-wide use of high- and low-affinity Tbrain transcription factor binding sites during echinoderm development

microRNA-1 regulates sea urchin skeletogenesis by directly targeting skeletogenic genes and modulating components of signaling pathways

microRNAs are evolutionarily conserved non-coding RNAs that direct post-transcriptional regulation of target transcripts. In vertebrates, microRNA-1 (miR-1) is expressed in muscle and has been found to play critical regulatory roles in vertebrate angiogenesis, a process that has been proposed to be analogous to sea urchin skeletogenesis. Results indicate that both miR-1 inhibitor and miR-1 mimic-injected larvae have significantly less F-actin enriched circumpharyngeal muscle fibers and fewer gut contractions. In addition, miR-1 regulates the positioning of skeletogenic primary mesenchyme cells (PMCs) and skeletogenesis of the sea urchin embryo. Interestingly, the gain-of-function of miR-1 leads to more severe PMC patterning and skeletal branching defects than its loss-of-function. The results suggest that miR-1 directly suppresses Ets1/2, Tbr, and VegfR7 of the skeletogenic gene regulatory network, and Nodal, and Wnt1 signaling components. This study identifies potential targets of miR-1 that impacts skeletogenesis and muscle formation and contributes to a deeper understanding of miR-1's function during development.

Low-affinity SPL binding sites contribute to subgenome expression divergence in allohexaploid wheat

Expression divergence caused by genetic variation and crosstalks among subgenomes of the allohexaploid bread wheat (Triticum aestivum. L., BBAADD) is hypothesized to increase its adaptability and/or plasticity. However, the molecular basis of expression divergence remains unclear. Squamosa promoter-binding protein-like (SPL) transcription factors are critical for a wide array of biological processes. In this study, we constructed expression regulatory networks by combining DAP-seq for 40 SPLs, ATAC-seq, and RNA-seq. Our findings indicate that a group of low-affinity SPL binding regions (SBRs) were targeted by diverse SPLs and caused different sequence preferences around the core GTAC motif. The SBRs including the low-affinity ones are evolutionarily conserved, enriched GWAS signals related to important agricultural traits. However, those SBRs are highly diversified among the cis-regulatory regions (CREs) of syntenic genes, with less than 8% SBRs coexisting in triad genes, suggesting that CRE variations are critical for subgenome differentiations. Knocking out of TaSPL7A/B/D and TaSPL15A/B/D subfamily further proved that both high- and low-affinity SBRs played critical roles in the differential expression of genes regulating tiller number and spike sizes. Our results have provided baseline data for downstream networks of SPLs and wheat improvements and revealed that CRE variations are critical sources for subgenome divergence in the allohexaploid wheat.

On the incongruence of genotype-phenotype and fitness landscapes

The mapping from genotype to phenotype to fitness typically involves multiple nonlinearities that can transform the effects of mutations. For example, mutations may contribute additively to a phenotype, but their effects on fitness may combine non-additively because selection favors a low or intermediate value of that phenotype. This can cause incongruence between the topographical properties of a fitness landscape and its underlying genotype-phenotype landscape. Yet, genotype-phenotype landscapes are often used as a proxy for fitness landscapes to study the dynamics and predictability of evolution. Here, we use theoretical models and empirical data on transcription factor-DNA interactions to systematically study the incongruence of genotype-phenotype and fitness landscapes when selection favors a low or intermediate phenotypic value. Using the theoretical models, we prove a number of fundamental results. For example, selection for low or intermediate phenotypic values does not change simple sign epistasis into reciprocal sign epistasis, implying that genotype-phenotype landscapes with only simple sign epistasis motifs will always give rise to single-peaked fitness landscapes under such selection. More broadly, we show that such selection tends to create fitness landscapes that are more rugged than the underlying genotype-phenotype landscape, but this increased ruggedness typically does not frustrate adaptive evolution because the local adaptive peaks in the fitness landscape tend to be nearly as tall as the global peak. Many of these results carry forward to the empirical genotype-phenotype landscapes, which may help to explain why low- and intermediate-affinity transcription factor-DNA interactions are so prevalent in eukaryotic gene regulation.

How do mutations change phenotypic traits and organismal fitness? This question is often addressed in the context of a classic metaphor of evolutionary theory—the fitness landscape. A fitness landscape is akin to a physical landscape, in which genotypes define spatial coordinates, and fitness defines the elevation of each coordinate. Evolution then acts like a hill-climbing process, in which populations ascend fitness peaks as a consequence of mutation and selection. It is becoming increasingly common to construct such landscapes using experimental data from high-throughput sequencing technologies and phenotypic assays, in systems such as macromolecules and gene regulatory circuits. Although these landscapes are typically defined by molecular phenotypes, and are therefore more appropriately referred to as genotype-phenotype landscapes, they are often used to study evolutionary dynamics. This requires the assumption that the molecular phenotype is a reasonable proxy for fitness, which need not be the case. For example, selection may favor a low or intermediate phenotypic value, causing incongruence between a fitness landscape and its underlying genotype-phenotype landscape. Here, we study such incongruence using a diversity of theoretical models and experimental data from gene regulatory systems. We regularly find incongruence, in that fitness landscapes tend to comprise more peaks than their underlying genotype-phenotype landscapes. However, using evolutionary simulations, we show that this increased ruggedness need not impede adaptation.

Evolutionary analyses of genes in Echinodermata offer insights towards the origin of metazoan phyla

Despite recent studies discussing the evolutionary impacts of gene duplications and losses among metazoans, the genomic basis for the evolution of phyla remains enigmatic. Here, we employ phylogenomic approaches to search for orthologous genes without known functions among echinoderms, and subsequently use them to guide the identification of their homologs across other metazoans. Our final set of 14 genes was obtained via a suite of homology prediction tools, gene expression data, gene ontology, and generating the Strongylocentrotus purpuratus phylome. The gene set was subjected to selection pressure analyses, which indicated that they are highly conserved and under negative selection. Their presence across broad taxonomic depths suggests that genes required to form a phylum are ancestral to that phylum. Therefore, rather than de novo gene genesis, we posit that evolutionary forces such as selection on existing genomic elements over large timescales may drive divergence and contribute to the emergence of phyla.

The arm of the starfish: The far-reaching applications of Patiria miniata as a model system in evolutionary, developmental, and regenerative biology

Many species of echinoderms have long been considered model research organisms in biology. Historically, much of this research has focused on the embryology of sea urchins and the use of their extensive gene regulatory networks as a tool to understand how the genome controls cell state specification and patterning. The establishment of Patiria miniata, the bat sea star, as a research organism has allowed us to expand on the concepts explored with sea urchins, viewing these genetic networks through a comparative lens, gaining great insight into the evolutionary mechanisms that shape developmental diversity. Extensive molecular tools have been developed in P. miniata, designed to explore gene expression dynamics and build gene regulatory networks. Echinoderms also have a robust set of bioinformatic and computational resources, centered around echinobase.org, an extensive database containing multiomic, developmental, and experimental resources for researchers. In addition to comparative evolutionary development, P. miniata is a promising system in its own right for studying whole body regeneration, metamorphosis and body plan development, as well as marine disease.

Regeneration of the larval sea star nervous system by wounding induced respecification to the sox2 lineage

The ability to restore lost body parts following traumatic injury is a fascinating area of biology that challenges current understanding of the ontogeny of differentiation. The origin of new cells needed to regenerate lost tissue, and whether they are pluripotent or have de- or trans-differentiated, remains one of the most important open questions . Additionally, it is not known whether developmental gene regulatory networks are reused or whether regeneration specific networks are deployed. Echinoderms, including sea stars, have extensive ability for regeneration, however, the technologies for obtaining transgenic echinoderms are limited and tracking cells involved in regeneration, and thus identifying the cellular sources and potencies has proven challenging. In this study, we develop new transgenic tools to follow the fate of populations of cells in the regenerating larva of the sea star Patiria miniata. We show that the larval serotonergic nervous system can regenerate following decapitation. Using a BAC-transgenesis approach we show that expression of the pan ectodermal marker, sox2, is induced in previously sox2 minus cells , even when cell division is inhibited. sox2+ cells give rise to new sox4+ neural precursors that then proceed along an embryonic neurogenesis pathway to reform the anterior nervous systems. sox2+ cells contribute to only neural and ectoderm lineages, indicating that these progenitors maintain their normal, embryonic lineage restriction. This indicates that sea star larval regeneration uses a combination of existing lineage restricted stem cells, as well as respecification of cells into neural lineages, and at least partial reuse of developmental GRNs to regenerate their nervous system.

Integration of 1:1 orthology maps and updated datasets into Echinobase

Echinobase (https://echinobase.org) is a central online platform that generates, manages and hosts genomic data relevant to echinoderm research. While the resource primarily serves the echinoderm research community, the recent release of an excellent quality genome for the frequently studied purple sea urchin (Strongylocentrotus purpuratus genome, v5.0) has provided an opportunity to adapt to the needs of a broader research community across other model systems. To this end, establishing pipelines to identify orthologous genes between echinoderms and other species has become a priority in many contexts including nomenclature, linking to data in other model organisms, and in internal functionality where data gathered in one hosted species can be associated with genes in other hosted echinoderms. This paper describes the orthology pipelines currently employed by Echinobase and how orthology data are processed to yield 1:1 ortholog mappings between a variety of echinoderms and other model taxa. We also describe functions of interest that have recently been included on the resource, including an updated developmental time course for S.purpuratus, and additional tracks for genome browsing. These data enhancements will increase the accessibility of the resource to non-echinoderm researchers and simultaneously expand the data quality and quantity available to core Echinobase users. Database URL: https://echinobase.org.

Genome-wide identification of binding sites and gene targets of Alx1, a pivotal regulator of echinoderm skeletogenesis

Alx1 is a conserved regulator of skeletogenesis in echinoderms and evolutionary changes in Alx1 sequence and expression have played a pivotal role in modifying programs of skeletogenesis within the phylum. Alx1 regulates a large suite of effector genes that control the morphogenetic behaviors and biomineral-forming activities of skeletogenic cells. To better understand the gene regulatory control of skeletogenesis by Alx1, we used genome-wide ChIP-seq to identify Alx1-binding sites and direct gene targets. Our analysis revealed that many terminal differentiation genes receive direct transcriptional inputs from Alx1. In addition, we found that intermediate transcription factors previously shown to be downstream of Alx1 all receive direct inputs from Alx1. Thus, Alx1 appears to regulate effector genes by indirect, as well as direct, mechanisms. We tested 23 high-confidence ChIP-seq peaks using GFP reporters and identified 18 active cis-regulatory modules (CRMs); this represents a high success rate for CRM discovery. Detailed analysis of a representative CRM confirmed that a conserved, palindromic Alx1-binding site was essential for expression. Our work significantly advances our understanding of the gene regulatory circuitry that controls skeletogenesis in sea urchins and provides a framework for evolutionary studies.

Cell type phylogenetics informs the evolutionary origin of echinoderm larval skeletogenic cell identity

The multiplicity of cell types comprising multicellular organisms begs the question as to how cell type identities evolve over time. Cell type phylogenetics informs this question by comparing gene expression of homologous cell types in distantly related taxa. We employ this approach to inform the identity of larval skeletogenic cells of echinoderms, a clade for which there are phylogenetically diverse datasets of spatial gene expression patterns. We determined ancestral spatial expression patterns of alx1, ets1, tbr, erg, and vegfr, key components of the skeletogenic gene regulatory network driving identity of the larval skeletogenic cell. Here we show ancestral state reconstructions of spatial gene expression of extant eleutherozoan echinoderms support homology and common ancestry of echinoderm larval skeletogenic cells. We propose larval skeletogenic cells arose in the stem lineage of eleutherozoans during a cell type duplication event that heterochronically activated adult skeletogenic cells in a topographically distinct tissue in early development.

Mechanism of cognate sequence discrimination by the ETS-family transcription factor ETS-1

Functional evidence increasingly implicates low-affinity DNA recognition by transcription factors as a general mechanism for the spatiotemporal control of developmental genes. Although the DNA sequence requirements for affinity are well-defined, the dynamic mechanisms that execute cognate recognition are much less resolved. To address this gap, here we examined ETS1, a paradigm developmental transcription factor, as a model for which cognate discrimination remains enigmatic. Using molecular dynamics simulations, we interrogated the DNA-binding domain of murine ETS1 alone and when bound to high-and low-affinity cognate sites or to nonspecific DNA. The results of our analyses revealed collective backbone and side-chain motions that distinguished cognate versus nonspecific as well as high- versus low-affinity cognate DNA binding. Combined with binding experiments with site-directed ETS1 mutants, the molecular dynamics data disclosed a triad of residues that respond specifically to low-affinity cognate DNA. We found that a DNA-contacting residue (Gln-336) specifically recognizes low-affinity DNA and triggers the loss of a distal salt bridge (Glu-343/Arg-378) via a large side-chain motion that compromises the hydrophobic packing of two core helices. As an intact Glu-343/Arg-378 bridge is the default state in unbound ETS1 and maintained in high-affinity and nonspecific complexes, the low-affinity complex represents a unique conformational adaptation to the suboptimization of developmental enhancers.

Genomic resources for the study of echinoderm development and evolution

Echinoderms are important research models for a wide range of biological questions. In particular, echinoderm embryos are exemplary models for dissecting the molecular and cellular processes that drive development and testing how these processes can be modified through evolution to produce the extensive morphological diversity observed in the phylum. Modern attempts to characterize these processes depend on some level of genomic analysis; from querying annotated gene sets to functional genomics experiments to identify candidate cis-regulatory sequences. Given how essential these data have become, it is important that researchers using available datasets or performing their own genome-scale experiments understand the nature and limitations of echinoderm genomic analyses. In this chapter we highlight the current state of echinoderm genomic data and provide methodological considerations for common approaches, including analysis of transcriptome and functional genomics datasets.

Methods for the experimental and computational analysis of gene regulatory networks in sea urchins

The discovery of gene regulatory networks (GRNs) has opened a gate to access the genomic mechanisms controlling development. GRNs are systems of transcriptional regulatory circuits that control the differential specification of cell fates during development by regulating gene expression. The experimental analysis of GRNs involves a collection of methods, each revealing aspects of the overall control process. This review provides an overview of experimental and computational methods that have been successfully applied for solving developmental GRNs in the sea urchin embryo. The key in this approach is to obtain experimental evidence for functional interactions between transcription factors and regulatory DNA. In the second part of this review, a more generally applicable strategy is discussed that shows a path from experimental evidence to annotation of regulatory linkages to the generation of GRN models.

Tuned polymerization of the transcription factor Yan limits off-DNA sequestration to confer context-specific repression

During development, transcriptional complexes at enhancers regulate gene expression in complex spatiotemporal patterns. To achieve robust expression without spurious activation, the affinity and specificity of transcription factor–DNA interactions must be precisely balanced. Protein–protein interactions among transcription factors are also critical, yet how their affinities impact enhancer output is not understood. The Drosophila transcription factor Yan provides a well-suited model to address this, as its function depends on the coordinated activities of two independent and essential domains: the DNA-binding ETS domain and the self-associating SAM domain. To explore how protein–protein affinity influences Yan function, we engineered mutants that increase SAM affinity over four orders of magnitude. This produced a dramatic subcellular redistribution of Yan into punctate structures, reduced repressive output and compromised survival. Cell-type specification and genetic interaction defects suggest distinct requirements for polymerization in different regulatory decisions. We conclude that tuned protein–protein interactions enable the dynamic spectrum of complexes that are required for proper regulation.

Conservation of DNA and ligand binding properties of retinoid X receptor from the placozoan Trichoplax adhaerens to human

Nuclear receptors are a superfamily of transcription factors restricted to animals. These transcription factors regulate a wide variety of genes with diverse roles in cellular homeostasis, development, and physiology. The origin and specificity of ligand binding within lineages of nuclear receptors (e.g., subfamilies) continues to be a focus of investigation geared toward understanding how the functions of these proteins were shaped over evolutionary history. Among early-diverging animal lineages, the retinoid X receptor (RXR) is first detected in the placozoan, Trichoplax adhaerens. To gain insight into RXR evolution, we characterized ligand- and DNA-binding activity of the RXR from T. adhaerens (TaRXR). Like bilaterian RXRs, TaRXR specifically bound 9-cis-retinoic acid, which is consistent with a recently published result and supports a conclusion that the ancestral RXR bound ligand. DNA binding site specificity of TaRXR was determined through protein binding microarrays (PBMs) and compared with human RXRɑ. The binding sites for these two RXR proteins were broadly conserved (∼85% shared high-affinity sequences within a targeted array), suggesting evolutionary constraint for the regulation of downstream genes. We searched for predicted binding motifs of the T. adhaerens genome within 1000 bases of annotated genes to identify potential regulatory targets. We identified 648 unique protein coding regions with predicted TaRXR binding sites that had diverse predicted functions, with enriched processes related to intracellular signal transduction and protein transport. Together, our data support hypotheses that the original RXR protein in animals bound a ligand with structural similarity to 9-cis-retinoic acid; the DNA motif recognized by RXR has changed little in more than 1 billion years of evolution; and the suite of processes regulated by this transcription factor diversified early in animal evolution.

From genome to anatomy: The architecture and evolution of the skeletogenic gene regulatory network of sea urchins and other echinoderms

The skeletogenic gene regulatory network (GRN) of sea urchins and other echinoderms is one of the most intensively studied transcriptional networks in any developing organism. As such, it serves as a preeminent model of GRN architecture and evolution. This review summarizes our current understanding of this developmental network. We describe in detail the most comprehensive model of the skeletogenic GRN, one developed for the euechinoid sea urchin Strongylocentrotus purpuratus, including its initial deployment by maternal inputs, its elaboration and stabilization through regulatory gene interactions, and its control of downstream effector genes that directly drive skeletal morphogenesis. We highlight recent comparative studies that have leveraged the euechinoid GRN model to examine the evolution of skeletogenic programs in diverse echinoderms, studies that have revealed both conserved and divergent features of skeletogenesis within the phylum. Last, we summarize the major insights that have emerged from analysis of the structure and evolution of the echinoderm skeletogenic GRN and identify key, unresolved questions as a guide for future work.

Towards a Dynamic Interaction Network of Life to unify and expand the evolutionary theory

The classic Darwinian theory and the Synthetic evolutionary theory and their linear models, while invaluable to study the origins and evolution of species, are not primarily designed to model the evolution of organisations, typically that of ecosystems, nor that of processes. How could evolutionary theory better explain the evolution of biological complexity and diversity? Inclusive network-based analyses of dynamic systems could retrace interactions between (related or unrelated) components. This theoretical shift from a Tree of Life to a Dynamic Interaction Network of Life, which is supported by diverse molecular, cellular, microbiological, organismal, ecological and evolutionary studies, would further unify evolutionary biology.

Omics approaches to study gene regulatory networks for development in echinoderms

Gene regulatory networks (GRNs) describe the interactions for a developmental process at a given time and space. Historically, perturbation experiments represent one of the key methods for analyzing and reconstructing a GRN, and the GRN governing early development in the sea urchin embryo stands as one of the more deeply dissected so far. As technology progresses, so do the methods used to address different biological questions. Next-generation sequencing (NGS) has become a standard experimental technique for genome and transcriptome sequencing and studies of protein-DNA interactions and DNA accessibility. While several efforts have been made toward the integration of different omics approaches for the study of the regulatory genome in many animals, in a few cases, these are applied with the purpose of reconstructing and experimentally testing developmental GRNs. Here, we review emerging approaches integrating multiple NGS technologies for the prediction and validation of gene interactions within echinoderm GRNs. These approaches can be applied to both 'model' and 'non-model' organisms. Although a number of issues still need to be addressed, advances in NGS applications, such as assay for transposase-accessible chromatin sequencing, combined with the availability of embryos belonging to different species, all separated by various evolutionary distances and accessible to experimental regulatory biology, place echinoderms in an unprecedented position for the reconstruction and evolutionary comparison of developmental GRNs. We conclude that sequencing technologies and integrated omics approaches allow the examination of GRNs on a genome-wide scale only if biological perturbation and cis-regulatory analyses are experimentally accessible, as in the case of echinoderm embryos.

ChIP-ping the branches of the tree: functional genomics and the evolution of eukaryotic gene regulation

Advances in the methods for detecting protein-DNA interactions have played a key role in determining the directions of research into the mechanisms of transcriptional regulation. The most recent major technological transformation happened a decade ago, with the move from using tiling arrays [chromatin immunoprecipitation (ChIP)-on-Chip] to high-throughput sequencing (ChIP-seq) as a readout for ChIP assays. In addition to the numerous other ways in which it is superior to arrays, by eliminating the need to design and manufacture them, sequencing also opened the door to carrying out comparative analyses of genome-wide transcription factor occupancy across species and studying chromatin biology in previously less accessible model and nonmodel organisms, thus allowing us to understand the evolution and diversity of regulatory mechanisms in unprecedented detail. Here, we review the biological insights obtained from such studies in recent years and discuss anticipated future developments in the field.

Echinoderm development and evolution in the post-genomic era

The highly recognizable animals within the phylum Echinodermata encompass an enormous disparity of adult and larval body plans. The extensive knowledge of sea urchin development has culminated in the description of the exquisitely detailed gene regulatory network (GRN) that governs the specification of various embryonic territories. This information provides a unique opportunity for comparative studies in other echinoderm taxa to understand the evolution and developmental mechanisms underlying body plan change. This review focuses on recent work that has utilized new genomic resources and systems-level experiments to address questions of evolutionary developmental biology. In particular, we synthesize the results of several recent studies from various echinoderm classes that have explored the development and evolution of the larval skeleton, which is a major feature that distinguishes the two predominant larval subtypes within the Phylum. We specifically examine the ways in which GRNs can evolve, either through cis regulatory and/or protein-level changes in transcription factors. We also examine recent work comparing evolution across shorter time scales that occur within and between species of sea urchin, and highlight the kinds of questions that can be addressed by these comparisons. The advent of new genomic and transcriptomic datasets in additional species from all classes of echinoderm will continue to empower the use of these taxa for evolutionary developmental studies.