Search

Innexins facilitate cell-cell communication by forming gap junctions or non-junctional hemichannels, which play important roles in metabolic, chemical, ionic, and electrical coupling. The lack of knowledge regarding the evolution and role of these channels in ctenophores (comb jellies), the likely sister group to the rest of animals, represents a substantial gap in our understanding of the evolution of intercellular communication in animals. Here we identify and phylogenetically characterize the complete set of innexins of four ctenophores:Mnemiopsis leidyi, Hormiphora californensis, Pleurobrachia bachei, andBeroe ovata. Our phylogenetic analyses suggest that ctenophore innexins diversified independently from those of other animals and were established early in the emergence of ctenophores. We identified a four-innexin genomic cluster, which was present in the last common ancestor of these four species and has been largely maintained in these lineages. Evidence from correlated spatial and temporal gene expression of theM. leidyiinnexin cluster suggest that this cluster has been maintained due to constraints related to gene regulation. We describe basic electrophysiological properties of putative ctenophore hemichannels from muscle cells using intracellular recording techniques, showing substantial overlap with the properties of bilaterian innexin channels. Together, our results suggest that the last common ancestor of animals had gap junctional channels also capable of forming functional innexin hemichannels, and that innexin genes have independently evolved in major lineages throughout Metazoa.

Hox and ParaHox transcription factors are important for specifying cell fates along the primary body axes during the development of most animals. Within Cnidaria, much of the research on Hox/ParaHox genes has focused on Anthozoa (anemones and corals) and Hydrozoa (hydroids) and has concentrated on the evolution and function of cnidarian Hox genes in relation to their bilaterian counterparts. Here we analyze together the full complement of Hox and ParaHox genes from species representing all four medusozoan classes (Staurozoa, Cubozoa, Hydrozoa, and Scyphozoa) and both anthozoan classes (Octocorallia and Hexacorallia). Our results show that Hox genes involved in patterning the directive axes of anthozoan polyps are absent in the stem leading to Medusozoa. For the first time, we show spatial and temporal expression patterns of Hox and ParaHox genes in the upside-down jellyfish Cassiopea xamachana (Scyphozoa), which are consistent with diversification of medusozoan Hox genes both from anthozoans and within medusozoa. Despite unprecedented taxon sampling, our phylogenetic analyses, like previous studies, are characterized by a lack of clear homology between most cnidarian and bilaterian Hox and Hox-related genes. Unlike previous studies, we propose the hypothesis that the cnidarian–bilaterian ancestor possessed a remarkably large Hox complement and that extensive loss of Hox genes was experienced by both cnidarian and bilaterian lineages.

Cnidarians are critical members of aquatic communities and have been an experimental system for a diversity of research areas ranging from development to biomechanics to global change biology. Yet we still lack a well-resolved, taxonomically balanced, cnidarian tree of life to place this research in appropriate phylogenetic context. To move towards this goal, we combined data from 26 new anthozoan transcriptomes with 86 previously published cnidarian and outgroup datasets to generate two 748-locus alignments containing 123,051 (trimmed) and 449,935 (untrimmed) amino acids. We estimated maximum likelihood phylogenies for both matrices under partitioned and unpartitioned site-homogeneous and site-heterogenous models of substitution. We used the resulting topology to constrain a phylogenetic analysis of 1,814 small subunit ribosomal (18S) gene sequences from GenBank. Our results confirm the position of Ceriantharia (tube-dwelling anemones), a historically recalcitrant group, as sister to the rest of Hexacorallia across all phylogenies regardless of data matrix or model choice. We also find unanimous support for the sister relationship of Endocnidozoa and Medusozoa and propose the name Operculozoa for the clade uniting these taxa. Our 18S hybrid phylogeny provides insight into relationships of 15% of extant taxa. Together these data are an invaluable resource for comparative cnidarian research and provide perspective to guide future refinement of cnidarian systematics.

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.

Six-state amino acid recoding strategies are commonly applied to combat the effects of compositional heterogeneity and substitution saturation in phylogenetic analyses. While these methods have been endorsed from a theoretical perspective, their performance has never been extensively tested. Here, we test the effectiveness of six-state recoding approaches by comparing the performance of analyses on recoded and non-recoded data sets that have been simulated under gradients of compositional heterogeneity or saturation. In our simulation analyses, non-recoding approaches consistently outperform six-state recoding approaches. Our results suggest that six-state recoding strategies are not effective in the face of high saturation. Furthermore, while recoding strategies do buffer the effects of compositional heterogeneity, the loss of information that accompanies six-state recoding outweighs its benefits. In addition, we evaluate recoding schemes with 9, 12, 15, and 18 states and show that these consistently outperform six-state recoding. Our analyses of other recoding schemes suggest that under conditions of very high compositional heterogeneity, it may be advantageous to apply recoding using more than six states, but we caution that applying any recoding should include sufficient justification. Our results have important implications for the more than 90 published papers that have incorporated six-state recoding, many of which have significant bearing on relationships across the tree of life. [Compositional heterogeneity; Dayhoff 6-state recoding; S&R 6-state recoding; six-state amino acid recoding; substitution saturation.]

Extreme environmental gradients represent excellent study systems to better understand the variables that mediate patterns of genomic variation between populations. They also allow for more accurate predictions of how future environmental change might affect marine species. The Persian/Arabian Gulf is extreme in both temperature and salinity, whereas the adjacent Gulf of Oman has conditions more typical of tropical oceans. The sea urchin Echinometra sp. EZ inhabits both of these seas and plays a critical role in coral reef health as a grazer and bioeroder, but, to date, there have been no population genomic studies on this or any urchin species in this unique region. E sp. EZ’s life history traits (e.g., large population sizes, large reproductive clutches, and long life spans), in theory, should homogenize populations unless nonneutral processes are occurring. Here, we generated a draft genome and a restriction site-associated DNA sequencing data set from seven populations along an environmental gradient across the Persian/Arabian Gulf and the Gulf of Oman. The estimated genome size of E. sp. EZ was 609 Mb and the heterozygosity was among the highest recorded for an echinoderm at 4.5%. We recovered 918 high-quality SNPs from 85 individuals which we then used in downstream analyses. Population structure analyses revealed a high degree of admixture between all sites, although there was population differentiation and significant pairwise FST values between the two seas. Preliminary results suggest migration is bidirectional between the seas and nine candidate loci were identified as being under putative natural selection, including one collagen gene. This study is the first to investigate the population genomics of a sea urchin from this extreme environmental gradient and is an important contribution to our understanding of the complex spatial patterns that drive genomic divergence.

Medusozoa is a widely distributed ancient lineage that harbors one-third of Cnidaria diversity divided into 4 classes. This clade is characterized by the succession of stages and modes of reproduction during metagenic lifecycles, and includes some of the most plastic body plans and life cycles among animals. The characterization of traditional genomic features, such as chromosome numbers and genome sizes, was rather overlooked in Medusozoa and many evolutionary questions still remain unanswered. Modern genomic DNA sequencing in this group started in 2010 with the publication of the Hydra vulgaris genome and has experienced an exponential increase in the past 3 years. Therefore, an update of the state of Medusozoa genomics is warranted. We reviewed different sources of evidence, including cytogenetic records and high-throughput sequencing projects. We focused on 4 main topics that would be relevant for the broad Cnidaria research community: (i) taxonomic coverage of genomic information; (ii) continuity, quality, and completeness of high-throughput sequencing datasets; (iii) overview of the Medusozoa specific research questions approached with genomics; and (iv) the accessibility of data and metadata. We highlight a lack of standardization in genomic projects and their reports, and reinforce a series of recommendations to enhance future collaborative research.

To date, genomic analyses in amoebozoans have been mostly limited to model organisms or medically important lineages. Consequently, the vast diversity of Amoebozoa genomes remain unexplored. A draft genome of Cochliopodium minus, an amoeba characterized by extensive cellular and nuclear fusions, is presented. C. minus has been a subject of recent investigation for its unusual sexual behavior. Cochliopodium’s sexual activity occurs during vegetative stage making it an ideal model for studying sexual development, which is sorely lacking in the group. Here we generate a C. minus draft genome assembly. From this genome, we detect a substantial number of lateral gene transfer (LGT) instances from bacteria (15%), archaea (0.9%) and viruses (0.7%) the majority of which are detected in our transcriptome data. We identify the complete meiosis toolkit genes in the C. minus genome, as well as the absence of several key genes involved in plasmogamy and karyogamy. Comparative genomics of amoebozoans reveals variation in sexual mechanism exist in the group. Similar to complex eukaryotes, C. minus (some amoebae) possesses Tyrosine kinases and duplicate copies of SPO11. We report a first example of alternative splicing in a key meiosis gene and draw important insights on molecular mechanism of sex in C. minus using genomic and transcriptomic data.

Cnidocytes are the explosive stinging cells found only in cnidarians (corals, jellyfish, etc). Specialized for prey capture and defense, cnidocytes are morphologically complex and vary widely in form and function across taxa; how such diversity evolved is unknown. Using CRISPR/Cas9-mediated genome editing in the burrowing sea anemoneNematostella vectensis, we show that a single transcription factor (NvSox2) acts as a binary switch between two alternative cnidocyte fates. Knockout ofNvSox2caused a complete transformation of nematocytes (piercing cells) into spirocytes (ensnaring cells). The type of spirocyte induced byNvSox2knockout (robust spirocyte) is not normally found inN. vectensisbut is common in sea anemones from other habitats. Homeotic control of cell fate provides a mechanistic explanation for the discontinuous distribution of cnidocyte types across cnidarians and demonstrates how simple counts of cell types can underestimate biodiversity.

Cnidocytes (“stinging cells”) are an unequivocally novel cell type used by cnidarians (corals, jellyfish, and their kin) to immobilize prey. Although they are known to share a common evolutionary origin with neurons, the developmental program that promoted the emergence of cnidocyte fate is not known. Using functional genomics in the sea anemone, Nematostella vectensis, we show that cnidocytes evolved by suppression of neural fate in a subset of neurons expressing RFamide. We further show that a single regulatory gene, a C2H2-type zinc finger transcription factor (ZNF845), coordinates both the gain of novel (cnidocyte-specific) traits and the inhibition of ancestral (neural) traits during cnidocyte development and that this gene arose by domain shuffling in the stem cnidarian. Thus, we uncover a mechanism by which a truly novel regulatory gene (ZNF845) promoted the origin of a truly novel cell type (cnidocyte) through duplication of an ancestral cell lineage (neuron) and inhibition of its ancestral identity (RFamide).SignificanceIn this study, we demonstrate how new cell types can arise in animals through duplication of an ancestral (old) cell type followed by functional divergence of the new daughter cell. Specifically, we show that stinging cells in cnidarians (jellyfish and corals) evolved by duplication of an ancestral neuron followed by inhibition of the RFamide neuropeptide it once secreted. This is the first evidence that stinging cells evolved from a specific subtype of neurons and suggests some neurons may be easier to co-opt for novel functions than others.

Placozoa are a morphologically simplistic group of marine animals found globally in tropical and subtropical environments. They consist of two named species, Trichoplax adhaerens and more recently Hoilungia hongkongensis, both with roughly six morphologically distinct cell types. With a sequenced genome, a limited number of cell types, and a simple flattened morphology, Trichoplax is an ideal model organism from which to explore the biology of an animal with a cellular complexity analagous to that of the earliest animals. Using a new approach for identification of gene expression patterns, this research looks at the relationship of Chordin/TgfΒ signaling and the axial patterning system of Placozoa. Our results suggest that placozoans have an oral–aboral axis similar to cnidarians and that the parahoxozoan ancestor (common ancestor of Placozoa and Cnidaria) was likely radially symmetric.

Graphics workloads make highly dynamic use of resources such as execution units (EUs), and thus can benefit from fast, fine-grain dynamic voltage and frequency scaling (DVFS) and retentive sleep. This paper presents a 14-nm graphics processing unit (GPU) prototype with modified EUs which include an integrated voltage regulator (IVR). The IVR enables energy-efficient EU turbo operation, data retention, and $V_{\text {MIN}}$ optimization per EU. Silicon measurements show that IVR-enabled EU turbo operation offers up to 32% (average 29%) energy reduction at constant performance.

Regeneration is one of the most fascinating and yet least understood processes of animals. Echinoderms, one of the closest related invertebrate groups to humans, can contribute to our understanding of the genetic basis of regenerative processes. Amongst echinoderms, sea cucumbers have the ability to grow back most of their body parts following injury, including the intestine and nervous tissue. The cellular and molecular events underlying these abilities in sea cucumbers have been most extensively studied in the species Holothuria glaberrima. However, research into the regenerative abilities of this species have been impeded due to the lack of adequate genomic resources. Here, we report the first draft genome assembly of H. glaberrima and demonstrate its value for future genetic studies. Using only short sequencing reads, we assembled the genome into 2,960,762 scaffolds totaling 1.5 gigabases with an N50 of 15 kilobases. Our BUSCO assessment of the genome resulted in 882 (90.2%) complete and partial genes from 978 genes queried. We incorporated transcriptomic data from several different life history stages to annotate 41,076 genes in our final assembly. To demonstrate the usefulness of the genome, we fully annotated the melanotransferrin (Mtf) gene family, which have a potential role in regeneration of the sea cucumber intestine. Using these same data, we extracted the mitochondrial genome, which showed high conservation to that of other holothuroids. Thus, these data will be a critical resource for ongoing studies of regeneration and other studies in sea cucumbers.

Following the completion of the genome sequencing and gene prediction of Mnemiopsis leidyi, a lobate ctenophore that is native to the coastal waters of the western Atlantic Ocean, we developed and implemented the Mnemiopsis Genome Project Portal (MGP Portal), a comprehensive Web-based data portal for navigating the genome sequence and gene annotations. In the years following the first release of the MGP Portal, it has become evident that the inclusion of data from significant published studies on Mnemiopsis has been critical to its adoption as the centralized resource for this emerging model organism. With this most recent update, the Portal has significantly expanded to include in situ images, temporal developmental expression profiles and single-cell expression data. Recent enhancements also include implementations of an updated BLAST interface, new graphical visualization tools and updates to gene pages that integrate all new data types. Database URL: https://research.nhgri.nih.gov/mnemiopsis/

The genomes of non-bilaterian metazoans are key to understanding the molecular basis of early animal evolution. However, a full comprehension of how animal-specific traits, such as nervous systems, arose is hindered by the scarcity and fragmented nature of genomes from key taxa, such as Porifera. Ephydatia muelleri is a freshwater sponge found across the northern hemisphere. Here, we present its 326 Mb genome, assembled to high contiguity (N50: 9.88 Mb) with 23 chromosomes on 24 scaffolds. Our analyses reveal a metazoan-typical genome architecture, with highly shared synteny across Metazoa, and suggest that adaptation to the extreme temperatures and conditions found in freshwater often involves gene duplication. The pancontinental distribution and ready laboratory culture of E. muelleri make this a highly practical model system which, with RNAseq, DNA methylation and bacterial amplicon data spanning its development and range, allows exploration of genomic changes both within sponges and in early animal evolution., Reconstructing the early molecular evolution of animals requires genomic resources for non-bilaterian animals. Here, the authors present the chromosome-level genome of a freshwater sponge together with analyses of its genome architecture, methylation, developmental gene expression, and microbiome.

The complete mitogenome of Echinometra sp. EZ has been described and fully annotated in this study. Phylogenetic analysis of cytochrome c oxidase subunit I (COI) from six Echinometra species confirms that our sample is E. sp. EZ. The mitogenome is 15,698 bp in length and contains 13 protein-coding genes, 22 tRNAs, 2 rRNAs, and a non-coding region with an identical organization to other Echinoidea. The E. sp. EZ mitogenome shared ∼99.1% identity to the published Echinometra mathaei mitogenome, differing by 147 SNPs. The E. sp. EZ mitogenome will serve as a resource that can be applied to disentangling the Echinometra species complex and to future population genetic studies of this ecologically important sea urchin species.

ANISEED (https://www.aniseed.cnrs.fr) is the main model organism database for the worldwide community of scientists working on tunicates, the vertebrate sister-group. Information provided for each species includes functionally-annotated gene and transcript models with orthology relationships within tunicates, and with echinoderms, cephalochordates and vertebrates. Beyond genes the system describes other genetic elements, including repeated elements and cis-regulatory modules. Gene expression profiles for several thousand genes are formalized in both wild-type and experimentally-manipulated conditions, using formal anatomical ontologies. These data can be explored through three complementary types of browsers, each offering a different view-point. A developmental browser summarizes the information in a gene- or territory-centric manner. Advanced genomic browsers integrate the genetic features surrounding genes or gene sets within a species. A Genomicus synteny browser explores the conservation of local gene order across deuterostome. This new release covers an extended taxonomic range of 14 species, including for the first time a non-ascidian species, the appendicularian Oikopleura dioica. Functional annotations, provided for each species, were enhanced through a combination of manual curation of gene models and the development of an improved orthology detection pipeline. Finally, gene expression profiles and anatomical territories can be explored in 4D online through the newly developed Morphonet morphogenetic browser.

BackgroundMutations in gene regulatory networks often lead to genetic divergence without impacting gene expression or developmental patterning. The rules governing this process of developmental systems drift, including the variable impact of selective constraints on different nodes in a gene regulatory network, remain poorly delineated.ResultsHere we examine developmental systems drift within the cardiopharyngeal gene regulatory networks of two tunicate species,Corella inflataandCiona robusta.Cross-species analysis of regulatory elements suggests thattrans-regulatory architecture is largely conserved between these highly divergent species. In contrast,cis-regulatory elements within this network exhibit distinct levels of conservation. In particular, while most of the regulatory elements we analyzed showed extensive rearrangements of functional binding sites, the enhancer for the cardiopharyngeal transcription factorFoxFis remarkably well-conserved. Even minor alterations in spacing between binding sites lead to loss ofFoxFenhancer function, suggesting that boundtrans-factors form position-dependent complexes.ConclusionsOur findings reveal heterogeneous levels of divergence across cardiopharyngealcis-regulatory elements. These distinct levels of divergence presumably reflect constraints that are not clearly associated with gene function or position within the regulatory network. Thus, levels ofcis-regulatory divergence or drift appear to be governed by distinct structural constraints that will be difficult to predict based on network architecture.

Six-state amino acid recoding strategies are commonly applied to combat the effects of compositional heterogeneity and substitution saturation in phylogenetic analyses. While these methods have been endorsed from a theoretical perspective, their performance has never been extensively tested. Here, we test the effectiveness of 6-state recoding approaches by comparing the performance of analyses on recoded and non-recoded datasets that have been simulated under gradients of compositional heterogeneity or saturation. In all of our simulation analyses, non-recoding approaches greatly outperformed 6-state recoding approaches. Our results suggest that 6-state recoding strategies are not effective in the face of high saturation. Further, while recoding strategies do buffer the effects of compositional heterogeneity, the loss of information that accompanies 6-state recoding outweighs its benefits, even in the most compositionally heterogeneous datasets. In addition, we evaluate recoding schemes with 9, 12, 15, and 18 states and show that these all outperform 6-state recoding. Our results have important implications for the more than 70 published papers that have incorporated 6-state recoding, many of which have significant bearing on relationships across the tree of life.

Addressing the origin of axial-patterning machinery is essential for understanding the evolution of animal form. Historically, sponges, a lineage that branched off early in animal evolution, were thought to lack Hox and ParaHox genes, suggesting that these critical axial-patterning genes arose after sponges diverged. However, a recent study has challenged this long-held doctrine by claiming to identify ParaHox genes (Cdx family) in two calcareous sponge species, Sycon ciliatum and Leucosolenia complicata. We reanalyzed the main data sets in this paper and analyzed an additional data set that expanded the number of bilaterians represented and removed outgroup homeodomains. As in the previous study, our Neighbor-Joining analyses of the original data sets recovered a clade that included sponge and Cdx genes, whereas Bayesian analyses placed these sponge genes within the NKL subclass of homeodomains. Unlike the original study, only one of our two maximum-likelihood analyses was congruent with Cdx genes in sponges. Our analyses of our additional data set led to the sponge genes consistently being placed within the NKL subclass of homeodomains regardless of method or model. Our results show more support for these sponge genes belonging to the NKL subclass, and therefore imply that Hox and ParaHox genes arose after Porifera diverged from the rest of animals.

Horseshoe crabs are xiphosuran chelicerates, the sister group to arachnids. As such, they are important for understanding the most recent common ancestor of Euchelicerata and the evolution and diversification of Arthropoda. Limulus polyphemus is the most investigated of the four extant species of horseshoe crabs, and the structure and function of its visual system have long been a major focus of studies critical for understanding the evolution of visual systems in arthropods. Likewise, studies of genes encoding Limulus opsins, the protein component of the visual pigments, are critical for understanding opsin evolution and diversification among chelicerates, where knowledge of opsins is limited, and more broadly among arthropods. In the present study, we sequenced and assembled a high quality nuclear genomic sequence of L. polyphemus and used these data to annotate the full repertoire of Limulus opsins. We conducted a detailed phylogenetic analysis of Limulus opsins, including using gene structure and synteny information to identify relationships among different opsin classes. We used our phylogeny to identify significant genomic events that shaped opsin evolution and therefore the visual system of Limulus. We also describe the tissue expression patterns of the 18 opsins identified and show that transcripts encoding a number, including a peropsin, are present throughout the central nervous system. In addition to significantly extending our understanding of photosensitivity in Limulus and providing critical insight into the genomic evolution of horseshoe crab opsins, this work provides a valuable genomic resource for addressing myriad questions related to xiphosuran physiology and arthropod evolution.