Analysis of the P. lividus sea urchin genome highlights contrasting trends of genomic and regulatory evolution in deuterostomes

Sea urchins are emblematic models in developmental biology and display several characteristics that set them apart from other deuterostomes. To uncover the genomic cues that may underlie these specificities, we generated a chromosome-scale genome assembly for the sea urchin Paracentrotus lividus and an extensive gene expression and epigenetic profiles of its embryonic development. We found that, unlike vertebrates, sea urchins retained ancestral chromosomal linkages but underwent very fast intrachromosomal gene order mixing. We identified a burst of gene duplication in the echinoid lineage and showed that some of these expanded genes have been recruited in novel structures (water vascular system, Aristotle's lantern, and skeletogenic micromere lineage). Finally, we identified gene-regulatory modules conserved between sea urchins and chordates. Our results suggest that gene-regulatory networks controlling development can be conserved despite extensive gene order rearrangement.

The Degenerate Tale of Ascidian Tails

Ascidians are invertebrate chordates, with swimming chordate tadpole larvae that have distinct heads and tails. The head contains the small brain, sensory organs, including the ocellus (light) and otolith (gravity) and the presumptive endoderm, while the tail has a notochord surrounded by muscle cells and a dorsal nerve cord. One of the chordate features is a post-anal tail. Ascidian tadpoles are nonfeeding, and their tails are critical for larval locomotion. After hatching the larvae swim up toward light and are carried by the tide and ocean currents. When competent to settle, ascidian tadpole larvae swim down, away from light, to settle and metamorphose into a sessile adult. Tunicates are classified as chordates because of their chordate tadpole larvae; in contrast, the sessile adult has a U-shaped gut and very derived body plan, looking nothing like a chordate. There is one group of ascidians, the Molgulidae, where many species are known to have tailless larvae. The Swalla Lab has been studying the evolution of tailless ascidian larvae in this clade for over 30 years and has shown that tailless larvae have evolved independently several times in this clade. Comparison of the genomes of two closely related species, the tailed Molgula oculata and tailless Molgula occulta reveals much synteny, but there have been multiple insertions and deletions that have disrupted larval genes in the tailless species. Genomics and transcriptomics have previously shown that there are pseudogenes expressed in the tailless embryos, suggesting that the partial rescue of tailed features in their hybrid larvae is due to the expression of intact genes from the tailed parent. Yet surprisingly, we find that the notochord gene regulatory network is mostly intact in the tailless M. occulta, although the notochord does not converge and extend and remains as an aggregate of cells we call the "notoball." We expect that eventually many of the larval gene networks will become evolutionarily lost in tailless ascidians and the larval body plan abandoned, with eggs developing directly into an adult. Here we review the current evolutionary and developmental evidence on how the molgulids lost their tails.

A cis-regulatory change underlying the motor neuron-specific loss of Ebf expression in immotile tunicate larvae

Many species in the tunicate family Molgulidae have independently lost their swimming larval form and instead develop as tailless, immotile larvae. These larvae do not develop structures that are essential for swimming such as the notochord, otolith, and tail muscles. However, little is known about neural development in these nonswimming larvae. Here, we studied the patterning of the Motor Ganglion (MG) of Molgula occulta, a nonswimming species. We found that spatial patterns of MG neuron regulators in this species are conserved, compared with species with swimming larvae, suggesting that the gene networks regulating their expression are intact despite the loss of swimming. However, expression of the key motor neuron regulatory gene Ebf (Collier/Olf/EBF) was reduced in the developing MG of M. occulta when compared with molgulid species with swimming larvae. This was corroborated by measuring allele-specific expression of Ebf in hybrid embryos from crosses of M. occulta with the swimming species M. oculata. Heterologous reporter construct assays in the model tunicate species Ciona robusta revealed a specific cis-regulatory sequence change that reduces expression of Ebf in the MG, but not in other cells. Taken together, these data suggest that MG neurons are still specified in M. occulta larvae, but their differentiation might be impaired due to reduction of Ebf expression levels.

Regulation of Neurogenesis by FGF Signaling and Neurogenin in the Invertebrate Chordate Ciona

Neurogenesis is a complex sequence of cellular processes and behaviors driven by the coordinated expression of conserved effectors. The bipolar tail neurons (BTNs) of Ciona develop according to a highly dynamic, yet highly stereotyped developmental program and thus could serve as an accessible model system for neurogenesis, including underlying cell behaviors like neuronal delamination, migration, and polarized axon outgrowth. Here we investigate both the upstream events that shape BTN neurogenesis through spatiotemporal regulation of the conserved proneural factor Neurog, spatiotemporal, and the gene expression profile of differentiating BTNs downstream of Neurog activity. We show that, although early FGF signaling is required for Neurog expression and BTN specification, Fgf8/17/18 is expressed in tail tip cells at later stages and suppresses sustained Neurog expression in the anterior BTN (aBTN) lineage, such that only one cell (the one furthest from the source of Fgf8/17/18) maintains Neurog expression and becomes a neuron. Curiously, Fgf8/17/18 might not affect neurogenesis of the posterior BTNs (pBTNs), which are in direct contact with the Fgf8/17/18-expressing cells. Finally, to profile gene expression associated with BTN neurogenesis we performed RNAseq of isolated BTN lineage cells in which BTN neurogenesis was enhanced or suppressed by perturbing Neurog function. This allowed us to identify several candidate genes that might play conserved roles in neurogenesis and neuronal migration in other animals, including mammals.

Effector gene expression underlying neuron subtype-specific traits in the Motor Ganglion of Ciona

The central nervous system of the Ciona larva contains only 177 neurons. The precise regulation of neuron subtype-specific morphogenesis and differentiation observed during the formation of this minimal connectome offers a unique opportunity to dissect gene regulatory networks underlying chordate neurodevelopment. Here we compare the transcriptomes of two very distinct neuron types in the hindbrain/spinal cord homolog of Ciona, the Motor Ganglion (MG): the Descending decussating neuron (ddN, proposed homolog of Mauthner Cells in vertebrates) and the MG Interneuron 2 (MGIN2). Both types are invariantly represented by a single bilaterally symmetric left/right pair of cells in every larva. Supernumerary ddNs and MGIN2s were generated in synchronized embryos and isolated by fluorescence-activated cell sorting for transcriptome profiling. Differential gene expression analysis revealed ddN- and MGIN2-specific enrichment of a wide range of genes, including many encoding potential "effectors" of subtype-specific morphological and functional traits. More specifically, we identified the upregulation of centrosome-associated, microtubule-stabilizing/bundling proteins and extracellular guidance cues part of a single intrinsic regulatory program that might underlie the unique polarization of the ddNs, the only descending MG neurons that cross the midline. Consistent with our predictions, CRISPR/Cas9-mediated, tissue-specific elimination of two such candidate effectors, Efcab6-related and Netrin1, impaired ddN polarized axon outgrowth across the midline.

Development of a goosegrass (Eleusine indica) draft genome and application to weed science research

BACKGROUND

Genomes are vital to the study of genomics, population genetics, and evolution of species. To date, only one genome (Echinochloa crus-galli) for C4 annual weedy grass species has been sequenced. Research was conducted to develop a draft genome of goosegrass (Eleusine indica; 2n = 2x = 18), one of the most common and troublesome weeds in the world.

RESULTS

A draft assembly of an approximately 492 Mb whole-genome sequence of goosegrass was obtained by de novo assembly of paired-end and mate-paired reads generated by Illumina sequencing of total genomic DNA. The genome was assembled into 24,072 scaffolds with N50 = 233,459 bp. More than 99% of transcriptome sequences were mapped to the goosegrass draft genome, and 95% of the commonly conserved plant genes were present. The assembled genome contains 25,467 unique protein-coding genes. Genes associated with herbicide resistance were obtained and variant calling allowed the detection of 754,409 single nucleotide polymorphisms. In addition, we also report 115,417 simple sequence repeats which can be deployed in population genetics and phylogenetic analysis.

CONCLUSION

This is the first report of genome sequence of goosegrass. Our assembly was able to identify all major herbicide-resistance related genes and develop a useful tool for other genomic and evolutionary analysis. © 2019 Society of Chemical Industry.

ANISEED 2017: extending the integrated ascidian database to the exploration and evolutionary comparison of genome-scale datasets

ANISEED (www.aniseed.cnrs.fr) is the main model organism database for tunicates, the sister-group of vertebrates. This release gives access to annotated genomes, gene expression patterns, and anatomical descriptions for nine ascidian species. It provides increased integration with external molecular and taxonomy databases, better support for epigenomics datasets, in particular RNA-seq, ChIP-seq and SELEX-seq, and features novel interactive interfaces for existing and novel datatypes. In particular, the cross-species navigation and comparison is enhanced through a novel taxonomy section describing each represented species and through the implementation of interactive phylogenetic gene trees for 60% of tunicate genes. The gene expression section displays the results of RNA-seq experiments for the three major model species of solitary ascidians. Gene expression is controlled by the binding of transcription factors to cis-regulatory sequences. A high-resolution description of the DNA-binding specificity for 131 Ciona robusta (formerly C. intestinalis type A) transcription factors by SELEX-seq is provided and used to map candidate binding sites across the Ciona robusta and Phallusia mammillata genomes. Finally, use of a WashU Epigenome browser enhances genome navigation, while a Genomicus server was set up to explore microsynteny relationships within tunicates and with vertebrates, Amphioxus, echinoderms and hemichordates.

Using ATAC-seq and RNA-seq to increase resolution in GRN connectivity

Echinoderms have some of the most complete reconstructed developmental gene regulatory networks (GRN) of any embryo, accounting for the formation of most embryo tissues and organs. Yet, many nodes (genes and regulators) and their regulatory interactions are still to be uncovered. Traditionally, knockdown/knockout experiments are performed to determine regulator-gene interactions, which are individually validated by cis-regulatory analysis. Differential RNA-seq, combined with perturbation analysis, allows for genome-wide reconstruction of a GRN around given regulators; however, this level of resolution cannot determine direct interactions. ChiP-chip or ChIP-seq is better equipped for determining, genome-wide, whether binding of a given transcription factor (TF) to cis-regulatory elements occurs. Antibodies for the TFs of interest must be available, and if not, this presents a limiting factor. ATAC-seq identifies regions of open chromatin, that are typically trimethylated at H3K4, H3K36 and H3K79 (Kouzarides, 2007), for a given time point, condition, or tissue. This technology combined with RNA-seq and perturbation analysis provides high resolution of the possible functional interactions occurring during development. Additionally, ATAC-seq is less expensive than ChIP-seq, requires less starting material, and provides a global view of regulatory regions. This chapter provides detailed steps to identify potential regulatory relationships between the nodes of a GRN, given a well assembled genome, annotated with gene models, and ATAC-seq data combined with RNA-seq and knockdown experiments.

The crowns have eyes: multiple opsins found in the eyes of the crown-of-thorns starfish Acanthaster planci

Background

Opsins are G protein-coupled receptors used for both visual and non-visual photoreception, and these proteins evolutionarily date back to the base of the bilaterians. In the current sequencing age, phylogenomic analysis has proven to be a powerful tool, facilitating the increase in knowledge about diversity within the opsin subclasses and, so far, at least nine types of opsins have been identified. Within echinoderms, opsins have been studied in Echinoidea and Ophiuroidea, which do not possess proper image forming eyes, but rather widely dispersed dermal photoreceptors. However, most species of Asteroidea, the starfish, possess true eyes and studying them will shed light on the diversity of opsin usage within echinoderms and help resolve the evolutionary history of opsins.

Results

Using high-throughput RNA sequencing, we have sequenced and analyzed the transcriptomes of different Acanthaster planci tissue samples: eyes, radial nerve, tube feet and a mixture of tissues from other organs. At least ten opsins were identified, and eight of them were found significantly differentially expressed in both eyes and radial nerve, with R-opsin being the most highly expressed in the eye.

Conclusion

This study provides new important insight into the involvement of opsins in visual and nonvisual photoreception. Of relevance, we found the first indication of an r-opsin photopigment expressed in a well-developed visual eye in a deuterostome animal. Additionally, we provided tissue specific A. planci transcriptomes that will aid in future Evo Devo studies.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1276-0) contains supplementary material, which is available to authorized users.

Developmental system drift in motor ganglion patterning between distantly related tunicates

BACKGROUND

The larval nervous system of the solitary tunicate Ciona is a simple model for the study of chordate neurodevelopment. The development and connectivity of the Ciona motor ganglion have been studied in fine detail, but how this important structure develops in other tunicates is not well known.

METHODS AND RESULTS

By comparing gene expression patterns in the developing MG of the distantly related tunicate Molgula occidentalis, we found that its patterning is highly conserved compared to the Ciona MG. MG neuronal subtypes in Molgula were specified in the exact same positions as in Ciona, though the timing of subtype-specific gene expression onset was slightly shifted to begin earlier, relative to mitotic exit and differentiation. In transgenic Molgula embryos electroporated with Dmbx reporter plasmids, we were also able to characterize the morphology of the lone pair of descending decussating neurons (ddNs) in Molgula, revealing the same unique contralateral projection seen in Ciona ddNs and their putative vertebrate homologs the Mauthner cells. Although Dmbx expression labels the ddNs in both species, cross-species transgenic assays revealed significant changes to the regulatory logic underlying Dmbx transcription. We found that Dmbx cis-regulatory DNAs from Ciona can drive highly specific reporter gene expression in Molgula ddNs, but Molgula sequences are not active in Ciona ddNs.

CONCLUSIONS

This acute divergence in the molecular mechanisms that underlie otherwise functionally conserved cis-regulatory DNAs supports the recently proposed idea that the extreme genetic plasticity observed in tunicates may be attributed to the extreme rigidity of the spatial organization of their embryonic cell lineages.

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.

The complete mitochondrial genome of the common pandora Pagellus erythrinus (Perciformes: Sparidae)

The common pandora (Pagellus erythrinus, Linnaeus 1758), one of the most popular sea bream species in the Mediterranean Sea, has high potential for aquaculture development. In this investigation, we analyzed the complete mitochondrial genome of P. erythrinus. The sequence has 16,828 bp in length and consists of 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a two non-coding regions (D-loop and L-origin). The overall nucleotide composition is: 27.5% of A, 28.2% of C, 27.5% of T, and 16.8% of G.

Evolutionary loss of melanogenesis in the tunicate Molgula occulta

BACKGROUND

Analyzing close species with diverse developmental modes is instrumental for investigating the evolutionary significance of physiological, anatomical and behavioral features at a molecular level. Many examples of trait loss are known in metazoan populations living in dark environments. Tunicates are the closest living relatives of vertebrates and typically present a lifecycle with distinct motile larval and sessile adult stages. The nervous system of the motile larva contains melanized cells associated with geotactic and light-sensing organs. It has been suggested that these are homologous to vertebrate neural crest-derived melanocytes. Probably due to ecological adaptation to distinct habitats, several species of tunicates in the Molgulidae family have tailless (anural) larvae that fail to develop sensory organ-associated melanocytes. Here we studied the evolution of Tyrosinase family genes, indispensible for melanogenesis, in the anural, unpigmented Molgula occulta and in the tailed, pigmented Molgula oculata by using phylogenetic, developmental and molecular approaches.

RESULTS

We performed an evolutionary reconstruction of the tunicate Tyrosinase gene family: in particular, we found that M. oculata possesses genes predicted to encode one Tyrosinase (Tyr) and three Tyrosinase-related proteins (Tyrps) while M. occulta has only Tyr and Tyrp.a pseudogenes that are not likely to encode functional proteins. Analysis of Tyr sequences from various M. occulta individuals indicates that different alleles independently acquired frameshifting short indels and/or larger mobile genetic element insertions, resulting in pseudogenization of the Tyr locus. In M. oculata, Tyr is expressed in presumptive pigment cell precursors as in the model tunicate Ciona robusta. Furthermore, a M. oculata Tyr reporter gene construct was active in the pigment cell precursors of C. robusta embryos, hinting at conservation of the regulatory network underlying Tyr expression in tunicates. In contrast, we did not observe any expression of the Tyr pseudogene in M. occulta embryos. Similarly, M. occulta Tyr allele expression was not rescued in pigmented interspecific M. occulta × M. oculata hybrid embryos, suggesting deleterious mutations also to its cis-regulatory sequences. However, in situ hybridization for transcripts from the M. occulta Tyrp.a pseudogene revealed its expression in vestigial pigment cell precursors in this species.

CONCLUSIONS

We reveal a complex evolutionary history of the melanogenesis pathway in tunicates, characterized by distinct gene duplication and loss events. Our expression and molecular data support a tight correlation between pseudogenization of Tyrosinase family members and the absence of pigmentation in the immotile larvae of M. occulta. These results suggest that relaxation of purifying selection has resulted in the loss of sensory organ-associated melanocytes and core genes in the melanogenesis biosynthetic pathway in M. occulta.

Complete plastid genome sequence of goosegrass (Eleusine indica) and comparison with other Poaceae

Eleusine indica, also known as goosegrass, is a serious weed in at least 42 countries. In this paper we report the complete plastid genome sequence of goosegrass obtained by de novo assembly of paired-end and mate-paired reads generated by Illumina sequencing of total genomic DNA. The goosegrass plastome is a circular molecule of 135,151bp in length, consisting of two single-copy regions separated by a pair of inverted repeats (IRs) of 20,919 bases. The large (LSC) and the small (SSC) single-copy regions span 80,667 bases and 12,646 bases, respectively. The plastome of goosegrass has 38.19% GC content and includes 108 unique genes, of which 76 are protein-coding, 28 are transfer RNA, and 4 are ribosomal RNA. The goosegrass plastome sequence was compared to eight other species of Poaceae. Although generally conserved with respect to Poaceae, this genomic resource will be useful for evolutionary studies within this weed species and the genus Eleusine.

ANISEED 2015: a digital framework for the comparative developmental biology of ascidians

Ascidians belong to the tunicates, the sister group of vertebrates and are recognized model organisms in the field of embryonic development, regeneration and stem cells. ANISEED is the main information system in the field of ascidian developmental biology. This article reports the development of the system since its initial publication in 2010. Over the past five years, we refactored the system from an initial custom schema to an extended version of the Chado schema and redesigned all user and back end interfaces. This new architecture was used to improve and enrich the description of Ciona intestinalis embryonic development, based on an improved genome assembly and gene model set, refined functional gene annotation, and anatomical ontologies, and a new collection of full ORF cDNAs. The genomes of nine ascidian species have been sequenced since the release of the C. intestinalis genome. In ANISEED 2015, all nine new ascidian species can be explored via dedicated genome browsers, and searched by Blast. In addition, ANISEED provides full functional gene annotation, anatomical ontologies and some gene expression data for the six species with highest quality genomes. ANISEED is publicly available at: http://www.aniseed.cnrs.fr.

Opsin evolution in the Ambulacraria

Opsins--G-protein coupled receptors involved in photoreception--have been extensively studied in the animal kingdom. The present work provides new insights into opsin-based photoreception and photoreceptor cell evolution with a first analysis of opsin sequence data for a major deuterostome clade, the Ambulacraria. Systematic data analysis, including for the first time hemichordate opsin sequences and an expanded echinoderm dataset, led to a robust opsin phylogeny for this cornerstone superphylum. Multiple genomic and transcriptomic resources were surveyed to cover each class of Hemichordata and Echinodermata. In total, 119 ambulacrarian opsin sequences were found, 22 new sequences in hemichordates and 97 in echinoderms (including 67 new sequences). We framed the ambulacrarian opsin repertoire within eumetazoan diversity by including selected reference opsins from non-ambulacrarians. Our findings corroborate the presence of all major ancestral bilaterian opsin groups in Ambulacraria. Furthermore, we identified two opsin groups specific to echinoderms. In conclusion, a molecular phylogenetic framework for investigating light-perception and photobiological behaviors in marine deuterostomes has been obtained.

Divergent mechanisms regulate conserved cardiopharyngeal development and gene expression in distantly related ascidians

Ascidians present a striking dichotomy between conserved phenotypes and divergent genomes: embryonic cell lineages and gene expression patterns are conserved between distantly related species. Much research has focused on Ciona or Halocynthia spp. but development in other ascidians remains poorly characterized. In this study, we surveyed the multipotent myogenic B7.5 lineage in Molgula spp. Comparisons to the homologous lineage in Ciona revealed identical cell division and fate specification events that result in segregation of larval, cardiac, and pharyngeal muscle progenitors. Moreover, the expression patterns of key regulators are conserved, but cross-species transgenic assays uncovered incompatibility, or ‘unintelligibility’, of orthologous cis-regulatory sequences between Molgula and Ciona. These sequences drive identical expression patterns that are not recapitulated in cross-species assays. We show that this unintelligibility is likely due to changes in both cis- and trans-acting elements, hinting at widespread and frequent turnover of regulatory mechanisms underlying otherwise conserved aspects of ascidian embryogenesis.

DOI:http://dx.doi.org/10.7554/eLife.03728.001

When two species have features that look similar, this may be because the features arise by the same processes during development. Other features may look similar yet develop by different mechanisms. ‘Developmental system drift’ refers to the process where a physical feature remains unaltered during evolution, but the underlying pathway that controls its development is changed. However, to date, there have been only a few experimental studies that support this idea.

Ascidians—also commonly known as sea squirts—are vase-like marine creatures, which start off as tadpole-like larvae that swim around until they find a place to settle down and attach themselves. Once attached, the sea squirts lose the ability to swim and start feeding, typically by filtering material out of the seawater. Sea squirts and their close relatives are the invertebrates (animals without backbones) that are most closely related to all vertebrates (animals with backbones), including humans. Furthermore, although different species of sea squirt have almost identical embryos, their genomes are very different.

Stolfi et al. have now studied whether developmental system drift may have occurred during the evolution of ascidians, by analyzing different species of sea squirt named Molgula and Ciona. Stolfi et al. compared the genomes of Molgula and Ciona and studied the expression of genes in the cells that give rise to the heart and the muscles of the head. As an embryo develops, specific genes are switched on or off, and these patterns of gene activation were broadly identical in the two species of sea squirt examined.

Enhancers are sequences of DNA that control when and how a gene is switched on. Given the similarities between the development of heart and head muscle cells in the different sea squirts, Stolfi et al. looked to see if the mechanisms of gene expression, and therefore the enhancers, were also conserved. Unexpectedly, this was not the case. When enhancers from Molgula were introduced into Ciona (and vice versa), these sequences were unable to switch on gene expression—thus enhancers from one sea squirt species could not function in the other.

Stolfi et al. conclude that the developmental systems may have drifted considerably during evolution of the sea squirts, in spite of their nearly identical embryos. This reinforces the view that different paths can lead to the formation of similar physical features.

DOI:http://dx.doi.org/10.7554/eLife.03728.002