Rooting the Animal Tree of Life

Unveiling the hidden viromes across the animal tree of life: insights from a taxonomic classification pipeline applied to invertebrates of 31 metazoan phyla

Invertebrates constitute the majority of animal species on Earth, including most disease-causing agents or vectors, with more diverse viromes when compared to vertebrates. Recent advancements in high-throughput sequencing have significantly expanded our understanding of invertebrate viruses, yet this knowledge remains biased toward a few well-studied animal lineages. In this study, we analyze invertebrate DNA and RNA viromes for 31 phyla using 417 publicly available RNA-Seq data sets from diverse environments in the marine-terrestrial and marine-freshwater gradients. This study aims to (i) estimate virome compositions at the family level for the first time across the animal tree of life, including the first exploration of the virome in several phyla, (ii) quantify the diversity of invertebrate viromes and characterize the structure of invertebrate-virus infection networks, and (iii) investigate host phylum and habitat influence on virome differences. Results showed that a set of few viral families of eukaryotes, comprising Retroviridae, Flaviviridae, and several families of giant DNA viruses, were ubiquitous and highly abundant. Nevertheless, some differences emerged between phyla, revealing for instance a less diverse virome in Ctenophora compared to the other animal phyla. Compositional analysis of the viromes showed that the host phylum explained over five times more variance in composition than its habitat. Moreover, significant similarities were observed between the viromes of some phylogenetically related phyla, which could highlight the influence of co-evolution in shaping invertebrate viromes.IMPORTANCEThis study significantly enhances our understanding of the global animal virome by characterizing the viromes of previously unexamined invertebrate lineages from a large number of animal phyla. It showcases the great diversity of viromes within each phylum and investigates the role of habitat shaping animal viral communities. Furthermore, our research identifies dominant virus families in invertebrates and distinguishes phyla with analogous viromes. This study sets the road toward a deeper understanding of the virome across the animal tree of life.

Early Animal Origin of BACE1 APP/Aβ Proteolytic Function

Alzheimer's disease is characterized, in part, by the accumulation of β-amyloid (Aβ) in the brain. Aβ is produced via the proteolysis of APP by BACE1 and γ-secretase. Since BACE1 is the rate-limiting enzyme in the production of Aβ, and a target for therapeutics, it is of interest to know when its proteolytic function evolved and for what purpose. Here, we take a functional evolutionary approach to show that BACE1 likely evolved from a gene duplication event near the base of the animal clade and that BACE1 APP/Aβ proteolytic function evolved during early animal diversification, hundreds of millions of years before the evolution of the APP/Aβ substrate. Our examination of BACE1 APP/Aβ proteolytic function includes cnidarians, ctenophores, and choanoflagellates. The most basal BACE1 ortholog is found in cnidarians, while ctenophores, placozoa, and choanoflagellates have genes equally orthologous to BACE1 and BACE2. BACE1 from a cnidarian (Hydra) can cleave APP to release Aβ, pushing back the date of the origin of its function to near the origin of animals. We tested more divergent BACE1/2 genes from a ctenophore (Mnemiopsis) and a choanoflagellate (Monosiga), and neither has this activity. These findings indicate that the specific proteolytic function of BACE1 evolved during the very earliest diversification of animals, most likely after a gene-duplication event.

Phylogenomic reconstruction of Solenogastres (Mollusca, Aplacophora) informs hypotheses on body size evolution

Body size is a fundamental characteristic of animals that impacts every aspect of their biology from anatomical complexity to ecology. In Mollusca, Solenogastres has been considered important to understanding the group's early evolution as most morphology-based phylogenetic reconstructions placed it as an early branching molluscan lineage. Under this scenario, molluscs were thought to have evolved from a small, turbellarian-like ancestor and small (i.e., macrofaunal) body size was inferred to be plesiomorphic for Solenogastres. More recently, phylogenomic studies have shown that aplacophorans (Solenogastres + Caudofoveata) form a clade with chitons (Polyplacophora), which is sister to all other molluscs, suggesting a relatively large-bodied (i.e., megafaunal) ancestor for Mollusca. Meanwhile, recent investigations into aplacophoran phylogeny have called the assumption that the last common ancestor of Solenogastres was small-bodied into question, but sampling of meiofaunal species was limited, biasing these studies towards large-bodied taxa and leaving fundamental questions about solenogaster body size evolution unanswered. Here, we supplemented available data with transcriptomes from eight diverse meiofaunal species of Solenogastres and conducted phylogenomic analyses on datasets of up to 949 genes. Maximum likelihood analyses support the meiofaunal family Meiomeniidae as the sister group to all other solenogasters, congruent with earlier ideas of a small-bodied ancestor of Solenogastres. In contrast, Bayesian Inference analyses support the large-bodied family Amphimeniidae as the sister group to all other solenogasters. Investigation of phylogenetic signal by comparing site-wise likelihood scores for the two competing hypotheses support the Meiomeniidae-first topology. In light of these results, we performed ancestral character state reconstruction to explore the implications of both hypotheses on understanding of Solenogaster evolution and review previous hypotheses about body size evolution and its potential consequences for solenogaster biology. Both hypotheses imply that body size evolution has been highly dynamic over the course of solenogaster evolution and that their relatively static body plan has successfully allowed for evolutionary transitions between meio-, macro- and megafaunal size ranges.

The promise and pitfalls of synteny in phylogenomics

Reconstructing the tree of life remains a central goal in biology. Early methods, which relied on small numbers of morphological or genetic characters, often yielded conflicting evolutionary histories, undermining confidence in the results. Investigations based on phylogenomics, which use hundreds to thousands of loci for phylogenetic inquiry, have provided a clearer picture of life's history, but certain branches remain problematic. To resolve difficult nodes on the tree of life, 2 recent studies tested the utility of synteny, the conserved collinearity of orthologous genetic loci in 2 or more organisms, for phylogenetics. Synteny exhibits compelling phylogenomic potential while also raising new challenges. This Essay identifies and discusses specific opportunities and challenges that bear on the value of synteny data and other rare genomic changes for phylogenomic studies. Synteny-based analyses of highly contiguous genome assemblies mark a new chapter in the phylogenomic era and the quest to reconstruct the tree of life.

Extracellular DNA traps in a ctenophore demonstrate immune cell behaviors in a non-bilaterian

The formation of extracellular DNA traps (ETosis) is a first response mechanism by specific immune cells following exposure to microbes. Initially characterized in vertebrate neutrophils, cells capable of ETosis have been discovered recently in diverse non-vertebrate taxa. To assess the conservation of ETosis between evolutionarily distant non-vertebrate phyla, we observed and quantified ETosis using the model ctenophore Mnemiopsis leidyi and the oyster Crassostrea gigas. Here we report that ctenophores - thought to have diverged very early from the metazoan stem lineage - possess immune-like cells capable of phagocytosis and ETosis. We demonstrate that both Mnemiopsis and Crassostrea immune cells undergo ETosis after exposure to diverse microbes and chemical agents that stimulate ion flux. We thus propose that ETosis is an evolutionarily conserved metazoan defense against pathogens.

Molecular evolution of the Thrombospondin superfamily

Thrombospondins (TSPs) are multidomain, calcium-binding glycoproteins that have wide-ranging roles in vertebrates in cell interactions, extracellular matrix (ECM) organisation, angiogenesis, tissue remodelling, synaptogenesis, and also in musculoskeletal and cardiovascular functions. Land animals encode five TSPs, which assembly co-translationally either as trimers (subgroup A) or pentamers (subgroup B). The vast majority of research has focused on this canonical TSP family, which evolved through the whole-genome duplications that took place early in the vertebrate lineage. With benefit of the growth in genome- and transcriptome-predicted proteomes of a much wider range of animal species, examination of TSPs throughout metazoan phyla has revealed extensive conservation of subgroup B-type TSPs in invertebrates. In addition, these searches established that canonical TSPs are, in fact, one branch within a TSP superfamily that includes other clades designated mega-TSPs, sushi-TSPs and poriferan-TSPs. Despite the apparent simplicity of poriferans and cnidarians as organisms, these phyla encode a greater diversity of TSP superfamily members than vertebrates. We discuss here the molecular characteristics of the TSP superfamily members, current knowledge of their expression profiles and functions in invertebrates, and models for the evolution of this complex ECM superfamily.

The origin and evolution of Wnt signalling

The Wnt signal transduction pathway has essential roles in the formation of the primary body axis during development, cellular differentiation and tissue homeostasis. This animal-specific pathway has been studied extensively in contexts ranging from developmental biology to medicine for more than 40 years. Despite its physiological importance, an understanding of the evolutionary origin and primary function of Wnt signalling has begun to emerge only recently. Recent studies on very basal metazoan species have shown high levels of conservation of components of both canonical and non-canonical Wnt signalling pathways. Furthermore, some pathway proteins have been described also in non-animal species, suggesting that recruitment and functional adaptation of these factors has occurred in metazoans. In this Review, we summarize the current state of research regarding the evolutionary origin of Wnt signalling, its ancestral function and the characteristics of the primal Wnt ligand, with emphasis on the importance of genomic studies in various pre-metazoan and basal metazoan species.

The ctenophore Mnemiopsis leidyi deploys a rapid injury response dating back to the last common animal ancestor

Regenerative potential is widespread but unevenly distributed across animals. However, our understanding of the molecular mechanisms underlying regenerative processes is limited to a handful of model organisms, restricting robust comparative analyses. Here, we conduct a time course of RNA-seq during whole body regeneration in Mnemiopsis leidyi (Ctenophora) to uncover gene expression changes that correspond with key events during the regenerative timeline of this species. We identified several genes highly enriched in this dataset beginning as early as 10 minutes after surgical bisection including transcription factors in the early timepoints, peptidases in the middle timepoints, and cytoskeletal genes in the later timepoints. We validated the expression of early response transcription factors by whole mount in situ hybridization, showing that these genes exhibited high expression in tissues surrounding the wound site. These genes exhibit a pattern of transient upregulation as seen in a variety of other organisms, suggesting that they may be initiators of an ancient gene regulatory network linking wound healing to the initiation of a regenerative response.

Functional analysis of ctenophore Shaker K+ channels: N-type inactivation in the animal roots

Here we explore the evolutionary origins of fast N-type ball-and-chain inactivation in Shaker (Kv1) K+ channels by functionally characterizing Shaker channels from the ctenophore (comb jelly) Mnemiopsis leidyi. Ctenophores are the sister lineage to other animals and Mnemiopsis has >40 Shaker-like K+ channels, but they have not been functionally characterized. We identified three Mnemiopsis channels (MlShak3-5) with N-type inactivation ball-like sequences at their N termini and functionally expressed them in Xenopus oocytes. Two of the channels, MlShak4 and MlShak5, showed rapid inactivation similar to cnidarian and bilaterian Shakers with rapid N-type inactivation, whereas MlShak3 inactivated ∼100-fold more slowly. Fast inactivation in MlShak4 and MlShak5 required the putative N-terminal inactivation ball sequences. Furthermore, the rate of fast inactivation in these channels depended on the number of inactivation balls/channel, but the rate of recovery from inactivation did not. These findings closely match the mechanism of N-type inactivation first described for Drosophila Shaker in which 1) inactivation balls on the N termini of each subunit can independently block the pore, and 2) only one inactivation ball occupies the pore binding site at a time. These findings suggest classical N-type activation evolved in Shaker channels at the very base of the animal phylogeny in a common ancestor of ctenophores, cnidarians, and bilaterians and that fast-inactivating Shakers are therefore a fundamental type of animal K+ channel. Interestingly, we find evidence from functional co-expression experiments and molecular dynamics that MlShak4 and MlShak5 do not co-assemble, suggesting that Mnemiopsis has at least two functionally independent N-type Shaker channels.

Brief History of Placozoa

Placozoans are morphologically the simplest free-living animals. They represent a unique window of opportunities to understand both the origin of the animal organization and the rules of life for the system and synthetic biology of the future. However, despite more than 100 years of their investigations, we know little about their organization, natural habitats, and life strategies. Here, we introduce this unique animal phylum and highlight some directions vital to broadening the frontiers of the biomedical sciences. In particular, understanding the genomic bases of placozoan biodiversity, cell identity, connectivity, reproduction, and cellular bases of behavior are critical hot spots for future studies.

Brief History of Ctenophora

Ctenophores are the descendants of the earliest surviving lineage of ancestral metazoans, predating the branch leading to sponges (Ctenophore-first phylogeny). Emerging genomic, ultrastructural, cellular, and systemic data indicate that virtually every aspect of ctenophore biology as well as ctenophore development are remarkably different from what is described in representatives of other 32 animal phyla. The outcome of this reconstruction is that most system-level components associated with the ctenophore organization result from convergent evolution. In other words, the ctenophore lineage independently evolved as high animal complexities with the astonishing diversity of cell types and structures as bilaterians and cnidarians. Specifically, neurons, synapses, muscles, mesoderm, through gut, sensory, and integrative systems evolved independently in Ctenophora. Rapid parallel evolution of complex traits is associated with a broad spectrum of unique ctenophore-specific molecular innovations, including alternative toolkits for making an animal. However, the systematic studies of ctenophores are in their infancy, and deciphering their remarkable morphological and functional diversity is one of the hot topics in biological research, with many anticipated surprises.

Ocean to Tree: Leveraging Single-Molecule RNA-Seq to Repair Genome Gene Models and Improve Phylogenomic Analysis of Gene and Species Evolution

Understanding gene evolution across genomes and organisms, including ctenophores, can provide unexpected biological insights. It enables powerful integrative approaches that leverage sequence diversity to advance biomedicine. Sequencing and bioinformatic tools can be inexpensive and user-friendly, but numerous options and coding can intimidate new users. Distinct challenges exist in working with data from diverse species but may go unrecognized by researchers accustomed to gold-standard genomes. Here, we provide a high-level workflow and detailed pipeline to enable animal collection, single-molecule sequencing, and phylogenomic analysis of gene and species evolution. As a demonstration, we focus on (1) PacBio RNA-seq of the genome-sequenced ctenophore Mnemiopsis leidyi, (2) diversity and evolution of the mechanosensitive ion channel Piezo in genetic models and basal-branching animals, and (3) associated challenges and solutions to working with diverse species and genomes, including gene model updating and repair using single-molecule RNA-seq. We provide a Python Jupyter Notebook version of our pipeline (GitHub Repository: Ctenophore-Ocean-To-Tree-2023 https://github.com/000generic/Ctenophore-Ocean-To-Tree-2023 ) that can be run for free in the Google Colab cloud to replicate our findings or modified for specific or greater use. Our protocol enables users to design new sequencing projects in ctenophores, marine invertebrates, or other novel organisms. It provides a simple, comprehensive platform that can ease new user entry into running their evolutionary sequence analyses.

Parallel Evolution of Transcription Factors in Basal Metazoans

Transcription factors (TFs) play a pivotal role as regulators of gene expression, orchestrating the formation and maintenance of diverse animal body plans and innovations. However, the precise contributions of TFs and the underlying mechanisms driving the origin of basal metazoan body plans, particularly in ctenophores, remain elusive. Here, we present a comprehensive catalog of TFs in 2 ctenophore species, Pleurobrachia bachei and Mnemiopsis leidyi, revealing 428 and 418 TFs in their respective genomes. In contrast, morphologically simpler metazoans have a reduced TF representation compared to ctenophores, cnidarians, and bilaterians: the sponge Amphimedon encodes 277 TFs, and the placozoan Trichoplax adhaerens encodes 274 TFs. The emergence of complex ctenophore tissues and organs coincides with significant lineage-specific diversification of the zinc finger C2H2 (ZF-C2H2) and homeobox superfamilies of TFs. Notable, the lineages leading to Amphimedon and Trichoplax exhibit independent expansions of leucine zipper (BZIP) TFs. Some lineage-specific TFs may have evolved through the domestication of mobile elements, thereby supporting alternative mechanisms of parallel TF evolution and body plan diversification across the Metazoa.

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

Phylogenetic and Protein Structure Analyses Provide Insight into the Evolution and Diversification of the CD36 Domain "Apex" among Scavenger Receptor Class B Proteins across Eukarya

The cluster of differentiation 36 (CD36) domain defines the characteristic ectodomain associated with class B scavenger receptor (SR-B) proteins. In bilaterians, SR-Bs play critical roles in diverse biological processes including innate immunity functions such as pathogen recognition and apoptotic cell clearance, as well as metabolic sensing associated with fatty acid uptake and cholesterol transport. Although previous studies suggest this protein family is ancient, SR-B diversity across Eukarya has not been robustly characterized. We analyzed SR-B homologs identified from the genomes and transcriptomes of 165 diverse eukaryotic species. The presence of highly conserved amino acid motifs across major eukaryotic supergroups supports the presence of a SR-B homolog in the last eukaryotic common ancestor. Our comparative analyses of SR-B protein structure identify the retention of a canonical asymmetric beta barrel tertiary structure within the CD36 ectodomain across Eukarya. We also identify multiple instances of independent lineage-specific sequence expansions in the apex region of the CD36 ectodomain-a region functionally associated with ligand-sensing. We hypothesize that a combination of both sequence expansion and structural variation in the CD36 apex region may reflect the evolution of SR-B ligand-sensing specificity between diverse eukaryotic clades.

Chemical cognition: chemoconnectomics and convergent evolution of integrative systems in animals

Neurons underpin cognition in animals. However, the roots of animal cognition are elusive from both mechanistic and evolutionary standpoints. Two conceptual frameworks both highlight and promise to address these challenges. First, we discuss evidence that animal neural and other integrative systems evolved more than once (convergent evolution) within basal metazoan lineages, giving us unique experiments by Nature for future studies. The most remarkable examples are neural systems in ctenophores and neuroid-like systems in placozoans and sponges. Second, in addition to classical synaptic wiring, a chemical connectome mediated by hundreds of signal molecules operates in tandem with neurons and is the most information-rich source of emerging properties and adaptability. The major gap-dynamic, multifunctional chemical micro-environments in nervous systems-is not understood well. Thus, novel tools and information are needed to establish mechanistic links between orchestrated, yet cell-specific, volume transmission and behaviors. Uniting what we call chemoconnectomics and analyses of the cellular bases of behavior in basal metazoan lineages arguably would form the foundation for deciphering the origins and early evolution of elementary cognition and intelligence.

Tissue Biology: In Search of a New Paradigm

Animal tissues are made up of multiple cell types that are increasingly well-characterized, yet our understanding of the core principles that govern tissue organization is still incomplete. This is in part because many observable tissue characteristics, such as cellular composition and spatial patterns, are emergent properties, and as such, they cannot be explained through the knowledge of individual cells alone. Here we propose a complex systems theory perspective to address this fundamental gap in our understanding of tissue biology. We introduce the concept of cell categories, which is based on cell relations rather than cell identity. Based on these notions we then discuss common principles of tissue modularity, introducing compositional, structural, and functional tissue modules. Cell diversity and cell relations provide a basis for a new perspective on the underlying principles of tissue organization in health and disease.

Compositional shifts associated with major evolutionary transitions in plants

Heterogeneity in gene trees, morphological characters, and composition has been associated with several major plant clades. Here, we examine heterogeneity in composition across a large transcriptomic dataset of plants to better understand whether locations of shifts in composition are shared across gene regions and whether directions of shifts within clades are shared across gene regions. We estimate mixed models of composition for both nucleotide and amino acids across a recent large-scale transcriptomic dataset for plants. We find shifts in composition across both nucleotide and amino acid datasets, with more shifts detected in nucleotides. We find that Chlorophytes and lineages within experience the most shifts. However, many shifts occur at the origins of land, vascular, and seed plants. While genes in these clades do not typically share the same composition, they tend to shift in the same direction. We discuss potential causes of these patterns. Compositional heterogeneity has been highlighted as a potential problem for phylogenetic analysis, but the variation presented here highlights the need to further investigate these patterns for the signal of biological processes.

Compositionally Constrained Sites Drive Long-Branch Attraction

Accurate phylogenies are fundamental to our understanding of the pattern and process of evolution. Yet, phylogenies at deep evolutionary timescales, with correspondingly long branches, have been fraught with controversy resulting from conflicting estimates from models with varying complexity and goodness of fit. Analyses of historical as well as current empirical datasets, such as alignments including Microsporidia, Nematoda, or Platyhelminthes, have demonstrated that inadequate modeling of across-site compositional heterogeneity, which is the result of biochemical constraints that lead to varying patterns of accepted amino acids along sequences, can lead to erroneous topologies that are strongly supported. Unfortunately, models that adequately account for across-site compositional heterogeneity remain computationally challenging or intractable for an increasing fraction of contemporary datasets. Here, we introduce "compositional constraint analysis," a method to investigate the effect of site-specific constraints on amino acid composition on phylogenetic inference. We show that more constrained sites with lower diversity and less constrained sites with higher diversity exhibit ostensibly conflicting signals under models ignoring across-site compositional heterogeneity that lead to long-branch attraction artifacts and demonstrate that more complex models accounting for across-site compositional heterogeneity can ameliorate this bias. We present CAT-posterior mean site frequencies (PMSF), a pipeline for diagnosing and resolving phylogenetic bias resulting from inadequate modeling of across-site compositional heterogeneity based on the CAT model. CAT-PMSF is robust against long-branch attraction in all alignments we have examined. We suggest using CAT-PMSF when convergence of the CAT model cannot be assured. We find evidence that compositionally constrained sites are driving long-branch attraction in two metazoan datasets and recover evidence for Porifera as the sister group to all other animals. [Animal phylogeny; cross-site heterogeneity; long-branch attraction; phylogenomics.].

Major Revisions in Pancrustacean Phylogeny and Evidence of Sensitivity to Taxon Sampling

The clade Pancrustacea, comprising crustaceans and hexapods, is the most diverse group of animals on earth, containing over 80% of animal species and half of animal biomass. It has been the subject of several recent phylogenomic analyses, yet relationships within Pancrustacea show a notable lack of stability. Here, the phylogeny is estimated with expanded taxon sampling, particularly of malacostracans. We show small changes in taxon sampling have large impacts on phylogenetic estimation. By analyzing identical orthologs between two slightly different taxon sets, we show that the differences in the resulting topologies are due primarily to the effects of taxon sampling on the phylogenetic reconstruction method. We compare trees resulting from our phylogenomic analyses with those from the literature to explore the large tree space of pancrustacean phylogenetic hypotheses and find that statistical topology tests reject the previously published trees in favor of the maximum likelihood trees produced here. Our results reject several clades including Caridoida, Eucarida, Multicrustacea, Vericrustacea, and Syncarida. Notably, we find Copepoda nested within Allotriocarida with high support and recover a novel relationship between decapods, euphausiids, and syncarids that we refer to as the Syneucarida. With denser taxon sampling, we find Stomatopoda sister to this latter clade, which we collectively name Stomatocarida, dividing Malacostraca into three clades: Leptostraca, Peracarida, and Stomatocarida. A new Bayesian divergence time estimation is conducted using 13 vetted fossils. We review our results in the context of other pancrustacean phylogenetic hypotheses and highlight 15 key taxa to sample in future studies.

Origin of the OAS-RNase L innate immune pathway before the rise of jawed vertebrates via molecular tinkering

Discriminating self from nonself is fundamental to immunity. Yet, it remains largely elusive how the mechanisms of self and nonself discrimination originated. Sensing double-stranded RNA as nonself, the 2',5'-oligoadenylate synthetase (OAS)-ribonuclease L (RNase L) pathway represents a crucial component of innate immunity. Here, we combine phylogenomic and functional analyses to show that the functional OAS-RNase L pathway likely originated through tinkering with preexisting proteins before the rise of jawed vertebrates during or before the Silurian period (444 to 419 Mya). Multiple concerted losses of OAS and RNase L occurred during the evolution of jawed vertebrates, further supporting the ancient coupling between OAS and RNase L. Moreover, both OAS and RNase L genes evolved under episodic positive selection across jawed vertebrates, suggesting a long-running evolutionary arms race between the OAS-RNase L pathway and microbes. Our findings illuminate how an innate immune pathway originated via molecular tinkering.

Confusion will be my epitaph: genome-scale discordance stifles phylogenetic resolution of Holothuroidea

Sea cucumbers (Holothuroidea) are a diverse clade of echinoderms found from intertidal waters to the bottom of the deepest oceanic trenches. Their reduced skeletons and limited number of phylogenetically informative traits have long obfuscated morphological classifications. Sanger-sequenced molecular datasets have also failed to constrain the position of major lineages. Noteworthy, topological uncertainty has hindered a resolution for Neoholothuriida, a highly diverse clade of Permo-Triassic age. We perform the first phylogenomic analysis of Holothuroidea, combining existing datasets with 13 novel transcriptomes. Using a highly curated dataset of 1100 orthologues, our efforts recapitulate previous results, struggling to resolve interrelationships among neoholothuriid clades. Three approaches to phylogenetic reconstruction (concatenation under both site-homogeneous and site-heterogeneous models, and coalescent-aware inference) result in alternative resolutions, all of which are recovered with strong support and across a range of datasets filtered for phylogenetic usefulness. We explore this intriguing result using gene-wise log-likelihood scores and attempt to correlate these with a large set of gene properties. While presenting novel ways of exploring and visualizing support for alternative trees, we are unable to discover significant predictors of topological preference, and our efforts fail to favour one topology. Neoholothuriid genomes seem to retain an amalgam of signals derived from multiple phylogenetic histories.

Ancient Rapid Radiation Explains Most Conflicts Among Gene Trees and Well-Supported Phylogenomic Trees of Nostocalean Cyanobacteria

Prokaryotic genomes are often considered to be mosaics of genes that do not necessarily share the same evolutionary history due to widespread horizontal gene transfers (HGTs). Consequently, representing evolutionary relationships of prokaryotes as bifurcating trees has long been controversial. However, studies reporting conflicts among gene trees derived from phylogenomic data sets have shown that these conflicts can be the result of artifacts or evolutionary processes other than HGT, such as incomplete lineage sorting, low phylogenetic signal, and systematic errors due to substitution model misspecification. Here, we present the results of an extensive exploration of phylogenetic conflicts in the cyanobacterial order Nostocales, for which previous studies have inferred strongly supported conflicting relationships when using different concatenated phylogenomic data sets. We found that most of these conflicts are concentrated in deep clusters of short internodes of the Nostocales phylogeny, where the great majority of individual genes have low resolving power. We then inferred phylogenetic networks to detect HGT events while also accounting for incomplete lineage sorting. Our results indicate that most conflicts among gene trees are likely due to incomplete lineage sorting linked to an ancient rapid radiation, rather than to HGTs. Moreover, the short internodes of this radiation fit the expectations of the anomaly zone, i.e., a region of the tree parameter space where a species tree is discordant with its most likely gene tree. We demonstrated that concatenation of different sets of loci can recover up to 17 distinct and well-supported relationships within the putative anomaly zone of Nostocales, corresponding to the observed conflicts among well-supported trees based on concatenated data sets from previous studies. Our findings highlight the important role of rapid radiations as a potential cause of strongly conflicting phylogenetic relationships when using phylogenomic data sets of bacteria. We propose that polytomies may be the most appropriate phylogenetic representation of these rapid radiations that are part of anomaly zones, especially when all possible genomic markers have been considered to infer these phylogenies. [Anomaly zone; bacteria; horizontal gene transfer; incomplete lineage sorting; Nostocales; phylogenomic conflict; rapid radiation; Rhizonema.].

The monoaminergic system is a bilaterian innovation

Monoamines like serotonin, dopamine, and adrenaline/noradrenaline (epinephrine/norepinephrine) act as neuromodulators in the nervous system. They play a role in complex behaviours, cognitive functions such as learning and memory formation, as well as fundamental homeostatic processes such as sleep and feeding. However, the evolutionary origin of the genes required for monoaminergic modulation is uncertain. Using a phylogenomic approach, in this study, we show that most of the genes involved in monoamine production, modulation, and reception originated in the bilaterian stem group. This suggests that the monoaminergic system is a bilaterian novelty and that its evolution may have contributed to the Cambrian diversification.

Recording cilia activity in ctenophores: effects of nitric oxide and low molecular weight transmitters

Cilia are the major effectors in Ctenophores, but very little is known about their transmitter control and integration. Here, we present a simple protocol to monitor and quantify cilia activity and provide evidence for polysynaptic control of cilia coordination in ctenophores. We also screened the effects of several classical bilaterian neurotransmitters (acetylcholine, dopamine, L-DOPA, serotonin, octopamine, histamine, gamma-aminobutyric acid (GABA), L-aspartate, L-glutamate, glycine), neuropeptide (FMRFamide), and nitric oxide (NO) on cilia beating in Pleurobrachia bachei and Bolinopsis infundibulum. NO and FMRFamide produced noticeable inhibitory effects on cilia activity, whereas other tested transmitters were ineffective. These findings further suggest that ctenophore-specific neuropeptides could be major candidates for signal molecules controlling cilia activity in representatives of this early-branching metazoan lineage.

Ancient gene linkages support ctenophores as sister to other animals

A central question in evolutionary biology is whether sponges or ctenophores (comb jellies) are the sister group to all other animals. These alternative phylogenetic hypotheses imply different scenarios for the evolution of complex neural systems and other animal-specific traits1-6. Conventional phylogenetic approaches based on morphological characters and increasingly extensive gene sequence collections have not been able to definitively answer this question7-11. Here we develop chromosome-scale gene linkage, also known as synteny, as a phylogenetic character for resolving this question12. We report new chromosome-scale genomes for a ctenophore and two marine sponges, and for three unicellular relatives of animals (a choanoflagellate, a filasterean amoeba and an ichthyosporean) that serve as outgroups for phylogenetic analysis. We find ancient syntenies that are conserved between animals and their close unicellular relatives. Ctenophores and unicellular eukaryotes share ancestral metazoan patterns, whereas sponges, bilaterians, and cnidarians share derived chromosomal rearrangements. Conserved syntenic characters unite sponges with bilaterians, cnidarians, and placozoans in a monophyletic clade to the exclusion of ctenophores, placing ctenophores as the sister group to all other animals. The patterns of synteny shared by sponges, bilaterians, and cnidarians are the result of rare and irreversible chromosome fusion-and-mixing events that provide robust and unambiguous phylogenetic support for the ctenophore-sister hypothesis. These findings provide a new framework for resolving deep, recalcitrant phylogenetic problems and have implications for our understanding of animal evolution.

Nitric oxide suppresses cilia activity in ctenophores

Cilia are the major effectors in Ctenophores, but very little is known about their transmitter control and integration. Here, we present a simple protocol to monitor and quantify cilia activity in semi-intact preparations and provide evidence for polysynaptic control of cilia coordination in ctenophores. Next, we screen the effects of several classical bilaterian neurotransmitters (acetylcholine, dopamine, L-DOPA, serotonin, octopamine, histamine, gamma-aminobutyric acid (GABA), L-aspartate, L-glutamate, glycine), neuropeptides (FMRFamide), and nitric oxide (NO) on cilia beating in Pleurobrachia bachei and Bolinopsis infundibulum . Only NO inhibited cilia beating, whereas other tested transmitters were ineffective. These findings further suggest that ctenophore-specific neuropeptides could be major candidate signaling molecules controlling cilia activity in representatives of this early-branching metazoan lineage.

Neurons that connect without synapses

The ctenophore nerve net suggests a complex evolutionary history of the animal nervous system.

Highly Dynamic Gene Family Evolution Suggests Changing Roles for PON Genes Within Metazoa

Change in gene family size has been shown to facilitate adaptation to different selective pressures. This includes gene duplication to increase dosage or diversification of enzymatic substrates and gene deletion due to relaxed selection. We recently found that the PON1 gene, an enzyme with arylesterase and lactonase activity, was lost repeatedly in different aquatic mammalian lineages, suggesting that the PON gene family is responsive to environmental change. We further investigated if these fluctuations in gene family size were restricted to mammals and approximately when this gene family was expanded within mammals. Using 112 metazoan protein models, we explored the evolutionary history of the PON family to characterize the dynamic evolution of this gene family. We found that there have been multiple, independent expansion events in tardigrades, cephalochordates, and echinoderms. In addition, there have been partial gene loss events in monotremes and sea cucumbers and what appears to be complete loss in arthropods, urochordates, platyhelminths, ctenophores, and placozoans. In addition, we show the mammalian expansion to three PON paralogs occurred in the ancestor of all mammals after the divergence of sauropsida but before the divergence of monotremes from therians. We also provide evidence of a novel PON expansion within the brushtail possum. In the face of repeated expansions and deletions in the context of changing environments, we suggest a range of selective pressures, including pathogen infection and mitigation of oxidative damage, are likely influencing the diversification of this dynamic gene family across metazoa.

Hydrodynamics in early animal evolution

Choanoflagellates and sponges feed by filtering microscopic particles from water currents created by the flagella of microvillar collar complexes situated on the cell bodies of the solitary or colonial choanoflagellates and on the choanocytes in sponges. The filtering mechanism has been known for more than a century, but only recently has the filtering process been studied in detail and also modelled, so that a detailed picture of the water currents has been obtained. In the solitary and most of the colonial choanoflagellates, the water flows freely around the cells, but in some forms, the cells are arranged in an open meshwork through which the water can be pumped. In the sponges, the choanocytes are located in choanocyte chambers (or choanocyte areas) with separate incurrent and excurrent canals/pores located in a larger body, which enables a fixed pattern of water currents through the collar complexes. Previous theories for the origin of sponges show evolutionary stages with choanocyte chambers without any opening or with only one opening, which makes separation of incurrent and excurrent impossible, and such stages must have been unable to feed. Therefore a new theory is proposed, which shows a continuous evolutionary lineage in which all stages are able to feed by means of the collar complexes.

The origin of human epithelial tissue

The histological structure of human epithelial tissue is complex, but all epithelia share three major features: cohesion, polarity and attachment. These functions are mainly achieved by the presence of specialized structures such as intercellular junctions, polarity protein complexes and basement membranes. In the present review, we have analyzed the presence of each of these structures in several groups of animals that are considered to be at the base of the animal evolution tree. Interestingly, these characters seem to have evolved independently, and a careful histological and structural analysis of the phylogenetic tree shows different groups of animals in which epithelia are absent and groups in which cells show only some of the specialized structures found in differentiated epithelia. These findings could contribute to understand how epithelial tissues evolved and determine their current protective functions.

Gradistics: An underappreciated dimension in evolutionary space

The growth of complexity is an unsolved and underappreciated problem. We consider possible causes of this growth, hypotheses testing, molecular mechanisms, complexity measures, cases of simplification, and significance for biomedicine. We focus on a general ability of regulation, which is based on the growing information storage and processing capacities, as the main proxy of complexity. Natural selection is indifferent to complexity. However, complexification can be inferred from the same first principle, on which natural selection is founded. Natural selection depends on potentially unlimited reproduction under limited environmental conditions. Because of the demographic pressure, the simple ecological niches become fulfilled and diversified (due to species splitting and divergence). Diversification increases complexity of biocenoses. After the filling and diversification of simple niches, the more complex niches can arise. This is the 'atomic orbitals' (AO) model. Complexity has many shortcomings but it has an advantage. This advantage is ability to regulatory adaptation, including behavioral, formed in the evolution by means of genetic adaptation. Regulatory adaptation is much faster than genetic one because it is based on the information previously accumulated via genetic adaptation and learning. Regulatory adaptation further increases complexity of biocenoses. This is the 'regulatory advantage' (RA) model. The comparison of both models allows testable predictions. We focus on the animal evolution because of the appearance of higher regulatory level (nervous system), which is absent in other lineages, and relevance to humans (including biomedical aspects).

Improving Orthologous Signal and Model Fit in Datasets Addressing the Root of the Animal Phylogeny

There is conflicting evidence as to whether Porifera (sponges) or Ctenophora (comb jellies) comprise the root of the animal phylogeny. Support for either a Porifera-sister or Ctenophore-sister tree has been extensively examined in the context of model selection, taxon sampling, and outgroup selection. The influence of dataset construction is comparatively understudied. We re-examine five animal phylogeny datasets that have supported either root hypothesis using an approach designed to enrich orthologous signal in phylogenomic datasets. We find that many component orthogroups in animal datasets fail to recover major lineages as monophyletic with the exception of Ctenophora, regardless of the supported root. Enriching these datasets to retain orthogroups recovering ≥3 major lineages reduces dataset size by up to 50% while retaining underlying phylogenetic information and taxon sampling. Site-heterogeneous phylogenomic analysis of these enriched datasets recovers both Porifera-sister and Ctenophora-sister positions, even with additional constraints on outgroup sampling. Two datasets which previously supported Ctenophora-sister support Porifera-sister upon enrichment. All enriched datasets display improved model fitness under posterior predictive analysis. While not conclusively rooting animals at either Porifera or Ctenophora, we do see an increase in signal for Porifera-sister and a decrease in signal for Ctenophore-sister when data are filtered for orthologous signal. Our results indicate that dataset size and construction as well as model fit influence animal root inference.

Functional analysis in a model sea anemone reveals phylogenetic complexity and a role in cnidocyte discharge of DEG/ENaC ion channels

Ion channels of the DEG/ENaC family share a similar structure but serve strikingly diverse biological functions, such as Na⁺ reabsorption, mechanosensing, proton-sensing, chemosensing and cell-cell communication via neuropeptides. This functional diversity raises the question of the ancient function of DEG/ENaCs. Using an extensive phylogenetic analysis across many different animal groups, we found a surprising diversity of DEG/ENaCs already in Cnidaria (corals, sea anemones, hydroids and jellyfish). Using a combination of gene expression analysis, electrophysiological and functional studies combined with pharmacological inhibition as well as genetic knockout in the model cnidarian Nematostella vectensis, we reveal an unanticipated role for a proton-sensitive DEG/ENaC in discharge of N. vectensis cnidocytes, the stinging cells typifying all cnidarians. Our study supports the view that DEG/ENaCs are versatile channels that have been co-opted for diverse functions since their early occurrence in animals and that respond to simple and ancient stimuli, such as omnipresent protons.

Exploring genome gene content and morphological analysis to test recalcitrant nodes in the animal phylogeny

An accurate phylogeny of animals is needed to clarify their evolution, ecology, and impact on shaping the biosphere. Although datasets of several hundred thousand amino acids are nowadays routinely used to test phylogenetic hypotheses, key deep nodes in the metazoan tree remain unresolved: the root of animals, the root of Bilateria, and the monophyly of Deuterostomia. Instead of using the standard approach of amino acid datasets, we performed analyses of newly assembled genome gene content and morphological datasets to investigate these recalcitrant nodes in the phylogeny of animals. We explored extensively the choices for assembling the genome gene content dataset and model choices of morphological analyses. Our results are robust to these choices and provide additional insights into the early evolution of animals, they are consistent with sponges as the sister group of all the other animals, the worm-like bilaterian lineage Xenacoelomorpha as the sister group of the other Bilateria, and tentatively support monophyletic Deuterostomia.

Phylotranscriptomics interrogation uncovers a complex evolutionary history for the planarian genus Dugesia (Platyhelminthes, Tricladida) in the Western Mediterranean

The Mediterranean is one of the most biodiverse areas of the Paleartic region. Here, basing on large data sets of single copy orthologs obtained from transcriptomic data, we investigated the evolutionary history of the genus Dugesia in the Western Mediterranean area. The results corroborated that the complex paleogeological history of the region was an important driver of diversification for the genus, speciating as microplates and islands were forming. These processes led to the differentiation of three main biogeographic clades: Iberia-Apennines-Alps, Corsica-Sardinia, and Iberia-Africa. The internal relationships of these major clades were analysed with several representative samples per species. The use of large data sets regarding the number of loci and samples, as well as state-of-the-art phylogenomic inference methods allowed us to answer different unresolved questions about the evolution of particular groups, such as the diversification path of D. subtentaculata in the Iberian Peninsula and its colonization of Africa. Additionally, our results support the differentiation of D. benazzii in two lineages which could represent two species. Finally, we analysed here for the first time a comprehensive number of samples from several asexual Iberian populations whose assignment at the species level has been an enigma through the years. The phylogenies obtained with different inference methods showed a branching topology of asexual individuals at the base of sexual clades. We hypothesize that this unexpected topology is related to long-term asexuality. This work represents the first phylotranscriptomic analysis of Tricladida, laying the first stone of the genomic era in phylogenetic studies on this taxonomic group.

Nitric oxide signaling in ctenophores

Nitric oxide (NO) is one of the most ancient and versatile signal molecules across all domains of life. NO signaling might also play an essential role in the origin of animal organization. Yet, practically nothing is known about the distribution and functions of NO-dependent signaling pathways in representatives of early branching metazoans such as Ctenophora. Here, we explore the presence and organization of NO signaling components using Mnemiopsis and kin as essential reference species. We show that NO synthase (NOS) is present in at least eight ctenophore species, including Euplokamis and Coeloplana, representing the most basal ctenophore lineages. However, NOS could be secondarily lost in many other ctenophores, including Pleurobrachia and Beroe. In Mnemiopsis leidyi, NOS is present both in adult tissues and differentially expressed in later embryonic stages suggesting the involvement of NO in developmental mechanisms. Ctenophores also possess soluble guanylyl cyclases as potential NO receptors with weak but differential expression across tissues. Combined, these data indicate that the canonical NO-cGMP signaling pathways existed in the common ancestor of animals and could be involved in the control of morphogenesis, cilia activities, feeding and different behaviors.

Transposon-derived transcription factors across metazoans

Transposable elements (TE) could serve as sources of new transcription factors (TFs) in plants and some other model species, but such evidence is lacking for most animal lineages. Here, we discovered multiple independent co-options of TEs to generate 788 TFs across Metazoa, including all early-branching animal lineages. Six of ten superfamilies of DNA transposon-derived conserved TF families (ZBED, CENPB, FHY3, HTH-Psq, THAP, and FLYWCH) were identified across nine phyla encompassing the entire metazoan phylogeny. The most extensive convergent domestication of potentially TE-derived TFs occurred in the hydroid polyps, polychaete worms, cephalopods, oysters, and sea slugs. Phylogenetic reconstructions showed species-specific clustering and lineage-specific expansion; none of the identified TE-derived TFs revealed homologs in their closest neighbors. Together, our study established a framework for categorizing TE-derived TFs and informing the origins of novel genes across phyla.

Step in Time: Conservation of Circadian Clock Genes in Animal Evolution

Over the past few decades, the molecular mechanisms responsible for circadian phenotypes of animals have been studied in increasing detail in mammals, some insects, and other invertebrates. Particular circadian proteins and their interactions are shared across evolutionary distant animals, resulting in a hypothesis for the canonical circadian clock of animals. As the number of species for which the circadian clockwork has been described increases, the circadian clock in animals driving cyclical phenotypes becomes less similar. Our focus in this review is to develop and synthesize the current literature to better understand the antiquity and evolution of the animal circadian clockwork. Here, we provide an updated understanding of circadian clock evolution in animals, largely through the lens of conserved genes characterized in the circadian clock identified in bilaterian species. These comparisons reveal extensive variation within the likely composition of the core clock mechanism, including losses of many genes, and that the ancestral clock of animals does not equate to the bilaterian clock. Despite the loss of these core genes, these species retain circadian behaviors and physiology, suggesting novel clocks have evolved repeatedly. Additionally, we highlight highly conserved cellular processes (e.g., cell division, nutrition) that intersect with the circadian clock of some animals. The conservation of these processes throughout the animal tree remains essentially unknown, but understanding their role in the evolution and maintenance of the circadian clock will provide important areas for future study.

Alternative neural systems: What is a neuron? (Ctenophores, sponges and placozoans)

How to make a neuron, a synapse, and a neural circuit? Is there only one 'design' for a neural architecture with a universally shared genomic blueprint across species? The brief answer is "No." Four early divergent lineages from the nerveless common ancestor of all animals independently evolved distinct neuroid-type integrative systems. One of these is a subset of neural nets in comb jellies with unique synapses; the second lineage is the well-known Cnidaria + Bilateria; the two others are non-synaptic neuroid systems in sponges and placozoans. By integrating scRNA-seq and microscopy data, we revise the definition of neurons as synaptically-coupled polarized and highly heterogenous secretory cells at the top of behavioral hierarchies with learning capabilities. This physiological (not phylogenetic) definition separates 'true' neurons from non-synaptically and gap junction-coupled integrative systems executing more stereotyped behaviors. Growing evidence supports the hypothesis of multiple origins of neurons and synapses. Thus, many non-bilaterian and bilaterian neuronal classes, circuits or systems are considered functional rather than genetic categories, composed of non-homologous cell types. In summary, little-explored examples of convergent neuronal evolution in representatives of early branching metazoans provide conceptually novel microanatomical and physiological architectures of behavioral controls in animals with prospects of neuro-engineering and synthetic biology.

Phylogenomics: Is less more when using large-scale datasets?

Phylogenetic studies have traditionally placed the simple Xenoacoelomorph worms as the sister group of all other animals with bilateral body symmetry. A new study shows that misidentification of orthologous genes might have been the source of at least some support for this placement.

A metazoan-specific C-terminal motif in EXC-4 and Gα-Rho/Rac signaling regulate cell outgrowth during tubulogenesis in C. elegans

Chloride intracellular channels (CLICs) are conserved proteins for which the cellular and molecular functions remain mysterious. An important insight into CLIC function came from the discovery that Caenorhabditis elegans EXC-4/CLIC regulates morphogenesis of the excretory canal (ExCa) cell, a single-cell tube. Subsequent work showed that mammalian CLICs regulate vascular development and angiogenesis, and human CLIC1 can rescue exc-4 mutants, suggesting conserved function in biological tube formation (tubulogenesis) and maintenance. However, the cell behaviors and signaling pathways regulated by EXC-4/CLICs during tubulogenesis in vivo remain largely unknown. We report a new exc-4 mutation, affecting a C-terminal residue conserved in virtually all metazoan CLICs, that reveals a specific role for EXC-4 in ExCa outgrowth. Cell culture studies suggest a function for CLICs in heterotrimeric G protein (Gα/β/γ)-Rho/Rac signaling, and Rho-family GTPases are common regulators of cell outgrowth. Using our new exc-4 mutant, we describe a previously unknown function for Gα-encoding genes (gpa-12/Gα12/13, gpa-7/Gαi, egl-30/Gαq and gsa-1/Gαs), ced-10/Rac and mig-2/RhoG in EXC-4-mediated ExCa outgrowth. Our results demonstrate that EXC-4/CLICs are primordial players in Gα-Rho/Rac-signaling, a pathway that is crucial for tubulogenesis in C. elegans and in vascular development.

Mass spectrometry of short peptides reveals common features of metazoan peptidergic neurons

The evolutionary origins of neurons remain unknown. Although recent genome data of extant early-branching animals have shown that neural genes existed in the common ancestor of animals, the physiological and genetic properties of neurons in the early evolutionary phase are still unclear. Here, we performed a mass spectrometry-based comprehensive survey of short peptides from early-branching lineages Cnidaria, Porifera and Ctenophora. We identified a number of mature ctenophore neuropeptides that are expressed in neurons associated with sensory, muscular and digestive systems. The ctenophore peptides are stored in vesicles in cell bodies and neurites, suggesting volume transmission similar to that of cnidarian and bilaterian peptidergic systems. A comparison of genetic characteristics revealed that the peptide-expressing cells of Cnidaria and Ctenophora express the vast majority of genes that have pivotal roles in maturation, secretion and degradation of neuropeptides in Bilateria. Functional analysis of neuropeptides and prediction of receptors with machine learning demonstrated peptide regulation of a wide range of target effector cells, including cells of muscular systems. The striking parallels between the peptidergic neuronal properties of Cnidaria and Bilateria and those of Ctenophora, the most basal neuron-bearing animals, suggest a common evolutionary origin of metazoan peptidergic nervous systems.

Renewed perspectives on the sedentary-pelagic last common bilaterian ancestor

Various evaluations of the last common bilaterian ancestor (lcba) currently suggest that it resembled either a microscopic, non-segmented motile adult; or, on the contrary, a complex segmented adult motile urbilaterian. These fundamental inconsistencies remain largely unexplained. A majority of multidisciplinary data regarding sedentary adult ancestral bilaterian organization is overlooked. The sedentary-pelagic model is supported now by a number of novel developmental, paleontological and molecular phylogenetic data: (1) data in support of sedentary sponges, in the adult stage, as sister to all other Metazoa; (2) a similarity of molecular developmental pathways in both adults and larvae across sedentary sponges, cnidarians, and bilaterians; (3) a cnidarian-bilaterian relationship, including a unique sharing of a bona fide Hox-gene cluster, of which the evolutionary appearance does not connect directly to a bilaterian motile organization; (4) the presence of sedentary and tube-dwelling representatives of the main bilaterian clades in the early Cambrian; (5) an absence of definite taxonomic attribution of Ediacaran taxa reconstructed as motile to any true bilaterian phyla; (6) a similarity of tube morphology (and the clear presence of a protoconch-like apical structure of the Ediacaran sedentary Cloudinidae) among shells of the early Cambrian, and later true bilaterians, such as semi-sedentary hyoliths and motile molluscs; (7) recent data that provide growing evidence for a complex urbilaterian, despite a continuous molecular phylogenetic controversy. The present review compares the main existing models and reconciles the sedentary model of an urbilaterian and the model of a larva-like lcba with a unified sedentary(adult)-pelagic(larva) model of the lcba.

Recoding amino acids to a reduced alphabet may increase or decrease phylogenetic accuracy

Common molecular phylogenetic characteristics such as long branches and compositional heterogeneity can be problematic for phylogenetic reconstruction when using amino acid data. Recoding alignments to reduced alphabets before phylogenetic analysis has often been used both to explore and potentially decrease the effect of such problems. We tested the effectiveness of this strategy on topological accuracy using simulated data on four-taxon trees. We simulated alignments in phylogenetically challenging ways to test the phylogenetic accuracy of analyses using various recoding strategies together with commonly-used homogeneous models. We tested three recoding methods based on amino acid exchangeability, and another recoding method based on lowering the compositional heterogeneity among alignment sequences as measured by the Chi-squared statistic. Our simulation results show that on trees with long branches where sequences approach saturation, accuracy was not greatly affected by exchangeability-based recodings, but Chi-squared-based recoding decreased accuracy. We then simulated sequences with different kinds of compositional heterogeneity over the tree. Recoding often increased accuracy on such alignments. Exchangeability-based recoding was rarely worse than not recoding, and often considerably better. Recoding based on lowering the Chi-squared value improved accuracy in some cases but not in others, suggesting that low compositional heterogeneity by itself is not sufficient to increase accuracy in the analysis of these alignments. We also simulated alignments using site-specific amino acid profiles, making sequences that had compositional heterogeneity over alignment sites. Exchangeability-based recoding coupled with site-homogeneous models had poor accuracy for these datasets but Chi-squared-based recoding on these alignments increased accuracy. We then simulated datasets that were compositionally both site- and tree-heterogeneous, like many real datasets. The effect on accuracy of recoding such doubly problematic datasets varied widely, depending on the type of compositional tree-heterogeneity and on the recoding scheme. Interestingly, analysis of unrecoded compositionally heterogeneous alignments with the NDCH or CAT models was generally more accurate than homogeneous analysis, whether recoded or not. Overall, our results suggest that making trees for recoded amino acid datasets can be useful, but they need to be interpreted cautiously as part of a more comprehensive analysis. The use of better fitting models like NDCH and CAT, which directly account for the patterns in the data, may offer a more promising long-term solution for analysing empirical data.

Multigenerational laboratory culture of pelagic ctenophores and CRISPR-Cas9 genome editing in the lobate Mnemiopsis leidyi

Despite long-standing experimental interest in ctenophores due to their unique biology, ecological influence and evolutionary status, previous work has largely been constrained by the periodic seasonal availability of wild-caught animals and difficulty in reliably closing the life cycle. To address this problem, we have developed straightforward protocols that can be easily implemented to establish long-term multigenerational cultures for biological experimentation in the laboratory. In this protocol, we describe the continuous culture of the Atlantic lobate ctenophore Mnemiopsis leidyi. A rapid 3-week egg-to-egg generation time makes Mnemiopsis suitable for a wide range of experimental genetic, cellular, embryological, physiological, developmental, ecological and evolutionary studies. We provide recommendations for general husbandry to close the life cycle of Mnemiopsis in the laboratory, including feeding requirements, light-induced spawning, collection of embryos and rearing of juveniles to adults. These protocols have been successfully applied to maintain long-term multigenerational cultures of several species of pelagic ctenophores, and can be utilized by laboratories lacking easy access to the ocean. We also provide protocols for targeted genome editing via microinjection with CRISPR-Cas9 that can be completed within ~2 weeks, including single-guide RNA synthesis, early embryo microinjection, phenotype assessment and sequence validation of genome edits. These protocols provide a foundation for using Mnemiopsis as a model organism for functional genomic analyses in ctenophores.

Comparative Single-Cell Transcriptomics Reveals Novel Genes Involved in Bivalve Embryonic Shell Formation and Questions Ontogenetic Homology of Molluscan Shell Types

Mollusks are known for their highly diverse repertoire of body plans that often includes external armor in form of mineralized hardparts. Representatives of the Conchifera, one of the two major lineages that comprises taxa which originated from a uni-shelled ancestor (Monoplacophora, Gastropoda, Cephalopoda, Scaphopoda, Bivalvia), are particularly relevant regarding the evolution of mollusk shells. Previous studies have found that the shell matrix of the adult shell (teleoconch) is rapidly evolving and that the gene set involved in shell formation is highly taxon-specific. However, detailed annotation of genes expressed in tissues involved in the formation of the embryonic shell (protoconch I) or the larval shell (protoconch II) are currently lacking. Here, we analyzed the genetic toolbox involved in embryonic and larval shell formation in the quagga mussel Dreissena rostriformis using single cell RNA sequencing. We found significant differences in genes expressed during embryonic and larval shell secretion, calling into question ontogenetic homology of these transitory bivalve shell types. Further ortholog comparisons throughout Metazoa indicates that a common genetic biomineralization toolbox, that was secondarily co-opted into molluscan shell formation, was already present in the last common metazoan ancestor. Genes included are engrailed, carbonic anhydrase, and tyrosinase homologs. However, we found that 25% of the genes expressed in the embryonic shell field of D. rostriformis lack an ortholog match with any other metazoan. This indicates that not only adult but also embryonic mollusk shells may be fast-evolving structures. We raise the question as to what degree, and on which taxonomic level, the gene complement involved in conchiferan protoconch formation may be lineage-specific or conserved across taxa.