Genomic data do not support comb jellies as the sister group to all other animals

Is Over-parameterization a Problem for Profile Mixture Models?

Biochemical constraints on the admissible amino acids at specific sites in proteins lead to heterogeneity of the amino acid substitution process over sites in alignments. It is well known that phylogenetic models of protein sequence evolution that do not account for site heterogeneity are prone to long-branch attraction (LBA) artifacts. Profile mixture models were developed to model heterogeneity of preferred amino acids at sites via a finite distribution of site classes each with a distinct set of equilibrium amino acid frequencies. However, it is unknown whether the large number of parameters in such models associated with the many amino acid frequency vectors can adversely affect tree topology estimates because of over-parameterization. Here, we demonstrate theoretically that for long sequences, over-parameterization does not create problems for estimation with profile mixture models. Under mild conditions, tree, amino acid frequencies, and other model parameters converge to true values as sequence length increases, even when there are large numbers of components in the frequency profile distributions. Because large sample theory does not necessarily imply good behavior for shorter alignments we explore the performance of these models with short alignments simulated with tree topologies that are prone to LBA artifacts. We find that over-parameterization is not a problem for complex profile mixture models even when there are many amino acid frequency vectors. In fact, simple models with few site classes behave poorly. Interestingly, we also found that misspecification of the amino acid frequency vectors does not lead to increased LBA artifacts as long as the estimated cumulative distribution function of the amino acid frequencies at sites adequately approximates the true one. In contrast, misspecification of the amino acid exchangeability rates can severely negatively affect parameter estimation. Finally, we explore the effects of including in the profile mixture model an additional "F-class" representing the overall frequencies of amino acids in the data set. Surprisingly, the F-class does not help parameter estimation significantly and can decrease the probability of correct tree estimation, depending on the scenario, even though it tends to improve likelihood scores.

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.

Highly conserved and extremely evolvable: BMP signalling in secondary axis patterning of Cnidaria and Bilateria

Bilateria encompass the vast majority of the animal phyla. As the name states, they are bilaterally symmetric, that is with a morphologically clear main body axis connecting their anterior and posterior ends, a second axis running between their dorsal and ventral surfaces, and with a left side being roughly a mirror image of their right side. Bone morphogenetic protein (BMP) signalling has widely conserved functions in the formation and patterning of the second, dorso-ventral (DV) body axis, albeit to different extents in different bilaterian species. Whilst initial findings in the fruit fly Drosophila and the frog Xenopus highlighted similarities amongst these evolutionarily very distant species, more recent analyses featuring other models revealed considerable diversity in the mechanisms underlying dorsoventral patterning. In fact, as phylogenetic sampling becomes broader, we find that this axis patterning system is so evolvable that even its core components can be deployed differently or lost in different model organisms. In this review, we will try to highlight the diversity of ways by which BMP signalling controls bilaterality in different animals, some of which do not belong to Bilateria. Future research combining functional analyses and modelling is bound to give us some understanding as to where the limits to the extent of the evolvability of BMP-dependent axial patterning may lie.

Decline and fall of the Ediacarans: late-Neoproterozoic extinctions and the rise of the modern biosphere

The end-Neoproterozoic transition marked a gradual but permanent shift between distinct configurations of Earth's biosphere. This interval witnessed the demise of the enigmatic Ediacaran Biota, ushering in the structured trophic webs and disparate animal body plans of Phanerozoic ecosystems. However, little consensus exists on the reality, drivers, and macroevolutionary implications of end-Neoproterozoic extinctions. Here we evaluate potential drivers of late-Neoproterozoic turnover by addressing recent findings on Ediacaran geochronology, the persistence of classical Ediacaran macrobionts into the Cambrian, and the existence of Ediacaran crown-group eumetazoans. Despite renewed interest in the possibility of Phanerozoic-style 'mass extinctions' in the latest Neoproterozoic, our synthesis of the available evidence does not support extinction models based on episodic geochemical triggers, nor does it validate simple ecological interpretations centred on direct competitive displacement. Instead, we argue that the protracted and indirect effects of early bilaterian innovations, including escalations in sediment engineering, predation, and the largely understudied impacts of reef-building, may best account for the temporal structure and possible selectivity of late-Neoproterozoic extinctions. We integrate these processes into a generalised model of early eumetazoan-dominated ecologies, charting the disruption of spatial and temporal isotropy on the Ediacaran benthos as a consequence of diversifying macrofaunal interactions. Given the nature of resource distribution in Ediacaran ecologies, the continuities among Ediacaran and Cambrian faunas, and the convergent origins of ecologically disruptive innovations among bilaterians we suggest that the rise of Phanerozoic-type biotas may have been unstoppable.

Divergence time estimates for the hypoxia-inducible factor-1 alpha (HIF1α) reveal an ancient emergence of animals in low-oxygen environments

Unveiling the tempo and mode of animal evolution is necessary to understand the links between environmental changes and biological innovation. Although the earliest unambiguous metazoan fossils date to the late Ediacaran period, molecular clock estimates agree that the last common ancestor (LCA) of all extant animals emerged ~850 Ma, in the Tonian period, before the oldest evidence for widespread ocean oxygenation at ~635-560 Ma in the Ediacaran period. Metazoans are aerobic organisms, that is, they are dependent on oxygen to survive. In low-oxygen conditions, most animals have an evolutionarily conserved pathway for maintaining oxygen homeostasis that triggers physiological changes in gene expression via the hypoxia-inducible factor (HIFa). However, here we confirm the absence of the characteristic HIFa protein domain responsible for the oxygen sensing of HIFa in sponges and ctenophores, indicating the LCA of metazoans lacked the functional protein domain as well, and so could have maintained their transcription levels unaltered under the very low-oxygen concentrations of their environments. Using Bayesian relaxed molecular clock dating, we inferred that the ancestral gene lineage responsible for HIFa arose in the Mesoproterozoic Era, ~1273 Ma (Credibility Interval 957-1621 Ma), consistent with the idea that important genetic machinery associated with animals evolved much earlier than the LCA of animals. Our data suggest at least two duplication events in the evolutionary history of HIFa, which generated three vertebrate paralogs, products of the two successive whole-genome duplications that occurred in the vertebrate LCA. Overall, our results support the hypothesis of a pre-Tonian emergence of metazoans under low-oxygen conditions, and an increase in oxygen response elements during animal evolution.

Hox, homology, and parsimony: An organismal perspective

Hox genes are important regulators in animal development. They often show a mosaic of conserved (e.g., longitudinal axis patterning) and lineage-specific novel functions (e.g., development of skeletal, sensory, or locomotory systems). Despite extensive research over the past decades, it remains controversial at which node in the animal tree of life the Hox cluster evolved. Its presence already in the last common metazoan ancestor has been proposed, although the genomes of both putative earliest extant metazoan offshoots, the ctenophores and the poriferans, are devoid of Hox sequences. The lack of Hox genes in the supposedly "simple"-built poriferans and their low number in cnidarians and the basally branching bilaterians, the xenacoelomorphs, seems to support the classical notion that the number of Hox genes is correlated with the degree of animal complexity. However, the 4-fold increase of the Hox cluster in xiphosurans, a basally branching chelicerate clade, as well as the situation in some teleost fishes that show a multitude of Hox genes compared to, e.g., human, demonstrates, that there is no per se direct correlation between organismal complexity and Hox number. Traditional approaches have tried to base homology on the morphological level on shared expression profiles of individual genes, but recent data have shown that, in particular with respect to Hox and other regulatory genes, complex gene-gene interactions rather than expression signatures of individual genes alone are responsible for shaping morphological traits during ontogeny. Accordingly, for sound homology assessments and reconstructions of character evolution on organ system level, additional independent datasets (e.g., morphological, developmental) need to be included in any such analyses. If supported by solid data, proposed structural homology should be regarded as valid and not be rejected solely on the grounds of non-parsimonious distribution of the character over a given phylogenetic topology.

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.

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.

Identifying the Best Approximating Model in Bayesian Phylogenetics: Bayes Factors, Cross-Validation or wAIC?

There is still no consensus as to how to select models in Bayesian phylogenetics, and more generally in applied Bayesian statistics. Bayes factors are often presented as the method of choice, yet other approaches have been proposed, such as cross-validation or information criteria. Each of these paradigms raises specific computational challenges, but they also differ in their statistical meaning, being motivated by different objectives: either testing hypotheses or finding the best-approximating model. These alternative goals entail different compromises, and as a result, Bayes factors, cross-validation, and information criteria may be valid for addressing different questions. Here, the question of Bayesian model selection is revisited, with a focus on the problem of finding the best-approximating model. Several model selection approaches were re-implemented, numerically assessed and compared: Bayes factors, cross-validation (CV), in its different forms (k-fold or leave-one-out), and the widely applicable information criterion (wAIC), which is asymptotically equivalent to leave-one-out cross-validation (LOO-CV). Using a combination of analytical results and empirical and simulation analyses, it is shown that Bayes factors are unduly conservative. In contrast, CV represents a more adequate formalism for selecting the model returning the best approximation of the data-generating process and the most accurate estimates of the parameters of interest. Among alternative CV schemes, LOO-CV and its asymptotic equivalent represented by the wAIC, stand out as the best choices, conceptually and computationally, given that both can be simultaneously computed based on standard Markov chain Monte Carlo runs under the posterior distribution. [Bayes factor; cross-validation; marginal likelihood; model comparison; wAIC.].

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.

Identifying and addressing methodological incongruence in phylogenomics: A review

The availability of phylogenetic data has greatly expanded in recent years. As a result, a new era in phylogenetic analysis is dawning-one in which the methods we use to analyse and assess our data are the bottleneck to producing valuable phylogenetic hypotheses, rather than the need to acquire more data. This makes the ability to accurately appraise and evaluate new methods of phylogenetic analysis and phylogenetic artefact identification more important than ever. Incongruence in phylogenetic reconstructions based on different datasets may be due to two major sources: biological and methodological. Biological sources comprise processes like horizontal gene transfer, hybridization and incomplete lineage sorting, while methodological ones contain falsely assigned data or violations of the assumptions of the underlying model. While the former provides interesting insights into the evolutionary history of the investigated groups, the latter should be avoided or minimized as best as possible. However, errors introduced by methodology must first be excluded or minimized to be able to conclude that biological sources are the cause. Fortunately, a variety of useful tools exist to help detect such misassignments and model violations and to apply ameliorating measurements. Still, the number of methods and their theoretical underpinning can be overwhelming and opaque. Here, we present a practical and comprehensive review of recent developments in techniques to detect artefacts arising from model violations and poorly assigned data. The advantages and disadvantages of the different methods to detect such misleading signals in phylogenetic reconstructions are also discussed. As there is no one-size-fits-all solution, this review can serve as a guide in choosing the most appropriate detection methods depending on both the actual dataset and the computational power available to the researcher. Ultimately, this informed selection will have a positive impact on the broader field, allowing us to better understand the evolutionary history of the group of interest.

The Hydra stem cell system - Revisited

Cnidarians are >600 million years old and are considered the sister group of Bilateria based on numerous molecular phylogenetic studies. Apart from Hydra, the genomes of all major clades of Cnidaria have been uncovered (e.g. Aurelia, Clytia, Nematostella and Acropora) and they reveal a remarkable completeness of the metazoan genomic toolbox. Of particular interest is Hydra, a model system of aging research, regenerative biology, and stem cell biology. With the knowledge gained from scRNA research, it is now possible to characterize the expression profiles of all cell types with great precision. In functional studies, our picture of the Hydra stem cell biology has changed, and we are in the process of obtaining a clear picture of the homeostasis and properties of the different stem cell populations. Even though Hydra is often compared to plant systems, the new data on germline and regeneration, but also on the dynamics and plasticity of the nervous system, show that Hydra with its simple body plan represents in a nutshell the prototype of an animal with stem cell lineages, whose properties correspond in many ways to Bilateria. This review provides an overview of the four stem cell lineages, the two epithelial lineages that constitute the ectoderm and the endoderm, as well as the multipotent somatic interstitial lineage (MPSC) and the germline stem cell lineage (GSC), also known as the interstitial cells of Hydra.

nRCFV: a new, dataset-size-independent metric to quantify compositional heterogeneity in nucleotide and amino acid datasets

MOTIVATION

Compositional heterogeneity-when the proportions of nucleotides and amino acids are not broadly similar across the dataset-is a cause of a great number of phylogenetic artefacts. Whilst a variety of methods can identify it post-hoc, few metrics exist to quantify compositional heterogeneity prior to the computationally intensive task of phylogenetic tree reconstruction. Here we assess the efficacy of one such existing, widely used, metric: Relative Composition Frequency Variability (RCFV), using both real and simulated data.

RESULTS

Our results show that RCFV can be biased by sequence length, the number of taxa, and the number of possible character states within the dataset. However, we also find that missing data does not appear to have an appreciable effect on RCFV. We discuss the theory behind this, the consequences of this for the future of the usage of the RCFV value and propose a new metric, nRCFV, which accounts for these biases. Alongside this, we present a new software that calculates both RCFV and nRCFV, called nRCFV_Reader.

AVAILABILITY AND IMPLEMENTATION

nRCFV has been implemented in RCFV_Reader, available at: https://github.com/JFFleming/RCFV_Reader . Both our simulation and real data are available at Datadryad: https://doi.org/10.5061/dryad.wpzgmsbpn .

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.

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.

From injury to patterning—MAPKs and Wnt signaling in Hydra

Hydra has a regenerative capacity that is not limited to individual organs but encompasses the entire body. Various global and integrative genome, transcriptome and proteome approaches have shown that many of the signaling pathways and transcription factors present in vertebrates are already present in Cnidaria, the sister group of Bilateria, and are also activated in regeneration. It is now possible to investigate one of the central questions of regeneration biology, i.e., how does the patterning system become activated by the injury signals that initiate regeneration. This review will present the current data obtained in Hydra and draw parallels with regeneration in Bilateria. Important findings of this global analysis are that the Wnt signaling pathway has a dual function in the regeneration process. In the early phase Wnt is activated generically and in a second phase of pattern formation it is activated in a position specific manner. Thus, Wnt signaling is part of the generic injury response, in which mitogen-activated protein kinases (MAPKs) are initially activated via calcium and reactive oxygen species (ROS). The MAPKs, p38, c-Jun N-terminal kinases (JNKs) and extracellular signal-regulated kinases (ERK) are essential for Wnt activation in Hydra head and foot regenerates. Furthermore, the antagonism between the ERK signaling pathway and stress-induced MAPKs results in a balanced induction of apoptosis and mitosis. However, the early Wnt genes are activated by MAPK signaling rather than apoptosis. Early Wnt gene activity is differentially integrated with a stable, β-Catenin-based gradient along the primary body axis maintaining axial polarity and activating further Wnts in the regenerating head. Because MAPKs and Wnts are highly evolutionarily conserved, we hypothesize that this mechanism is also present in vertebrates but may be activated to different degrees at the level of early Wnt gene integration.

Full-Length Transcriptome Maps of Reef-Building Coral Illuminate the Molecular Basis of Calcification, Symbiosis, and Circa-Dian Genes

Coral transcriptomic data largely rely on short-read sequencing, which severely limits the understanding of coral molecular mechanisms and leaves many important biological questions unresolved. Here, we sequence the full-length transcriptomes of four common and frequently dominant reef-building corals using the PacBio Sequel II platform. We obtain information on reported gene functions, structures, and expression profiles. Among them, a comparative analysis of biomineralization-related genes provides insights into the molecular basis of coral skeletal density. The gene expression profiles of the symbiont Symbiodiniaceae are also isolated and annotated from the holobiont sequence data. Finally, a phylogenetic analysis of key circadian clock genes among 40 evolutionarily representative species indicates that there are four key members in early metazoans, including cry genes; Clock or Npas2; cyc or Arntl; and tim, while per, as the fifth member, occurs in Bilateria. In summary, this work provides a foundation for further work on the manipulation of skeleton production or symbiosis to promote the survival of these important organisms.

A Practical Guide to Design and Assess a Phylogenomic Study

Over the last decade, molecular systematics has undergone a change of paradigm as high-throughput sequencing now makes it possible to reconstruct evolutionary relationships using genome-scale datasets. The advent of "big data" molecular phylogenetics provided a battery of new tools for biologists but simultaneously brought new methodological challenges. The increase in analytical complexity comes at the price of highly specific training in computational biology and molecular phylogenetics, resulting very often in a polarized accumulation of knowledge (technical on one side and biological on the other). Interpreting the robustness of genome-scale phylogenetic studies is not straightforward, particularly as new methodological developments have consistently shown that the general belief of "more genes, more robustness" often does not apply, and because there is a range of systematic errors that plague phylogenomic investigations. This is particularly problematic because phylogenomic studies are highly heterogeneous in their methodology, and best practices are often not clearly defined. The main aim of this article is to present what I consider as the ten most important points to take into consideration when planning a well-thought-out phylogenomic study and while evaluating the quality of published papers. The goal is to provide a practical step-by-step guide that can be easily followed by nonexperts and phylogenomic novices in order to assess the technical robustness of phylogenomic studies or improve the experimental design of a project.

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.

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.

Ontogeny, Phylotypic Periods, Paedomorphosis, and Ontogenetic Systematics

The key terms linking ontogeny and evolution are briefly reviewed. It is shown that their application and usage in the modern biology are often inconsistent and incorrectly understood even within the “evo-devo” field. For instance, the core modern reformulation that ontogeny not merely recapitulates, but produces phylogeny implies that ontogeny and phylogeny are closely interconnected. However, the vast modern phylogenetic and taxonomic fields largely omit ontogeny as a central concept. Instead, the common “clade-” and “tree-thinking” prevail, despite on the all achievements of the evo-devo. This is because the main conceptual basis of the modern biology is fundamentally ontogeny-free. In another words, in the Haeckel’s pair of “ontogeny and phylogeny,” ontogeny is still just a subsidiary for the evolutionary process (and hence, phylogeny), instead as in reality, its main driving force. The phylotypic periods is another important term of the evo-devo and represent a modern reformulation of Haeckel’s recapitulations and biogenetic law. However, surprisingly, this one of the most important biological evidence, based on the natural ontogenetic grounds, in the phylogenetic field that can be alleged as a “non-evolutionary concept.” All these observations clearly imply that a major revision of the main terms which are associated with the “ontogeny and phylogeny/evolution” field is urgently necessarily. Thus, “ontogenetic” is not just an endless addition to the term “systematics,” but instead a crucial term, without it neither systematics, nor biology have sense. To consistently employ the modern ontogenetic and epigenetic achievements, the concept of ontogenetic systematics is hereby refined. Ontogenetic systematics is not merely a “research program” but a key biological discipline which consistently links the enormous biological diversity with underlying fundamental process of ontogeny at both molecular and morphological levels. The paedomorphosis is another widespread ontogenetic-and-evolutionary process that is significantly underestimated or misinterpreted by the current phylogenetics and taxonomy. The term paedomorphosis is refined, as initially proposed to link ontogeny with evolution, whereas “neoteny” and “progenesis” are originally specific, narrow terms without evolutionary context, and should not be used as synonyms of paedomorphosis. Examples of application of the principles of ontogenetic systematics represented by such disparate animal groups as nudibranch molluscs and ophiuroid echinoderms clearly demonstrate that perseverance of the phylotypic periods is based not only on the classic examples in vertebrates, but it is a universal phenomenon in all organisms, including disparate animal phyla.

Ctenophores are direct developers that reproduce continuously beginning very early after hatching

A substantial body of literature reports that ctenophores exhibit an apparently unique life history characterized by biphasic sexual reproduction, the first phase of which is called larval reproduction or dissogeny. Whether this strategy is plastically deployed or a typical part of these species’ life history was unknown. In contrast to previous reports, we show that the ctenophore Mnemiopsis leidyi does not have separate phases of early and adult reproduction, regardless of the morphological transition to what has been considered the adult form. Rather, these ctenophores begin to reproduce at a small body size and spawn continuously from this point onward under adequate environmental conditions. They do not display a gap in productivity for metamorphosis or other physiological transition at a certain body size. Furthermore, nutritional and environmental constraints on fecundity are similar in both small and large animals. Our results provide critical parameters for understanding resource partitioning between growth and reproduction in this taxon, with implications for management of this species in its invaded range. Finally, we report an observation of similarly small-size spawning in a beroid ctenophore, which is morphologically, ecologically, and phylogenetically distinct from other ctenophores reported to spawn at small sizes. We conclude that spawning at small body size should be considered as the default, on-time developmental trajectory rather than as precocious, stress-induced, or otherwise unusual for ctenophores. The ancestral ctenophore was likely a direct developer, consistent with the hypothesis that multiphasic life cycles were introduced after the divergence of the ctenophore lineage.

Molecular dating of the blood pigment hemocyanin provides new insight into the origin of animals

The Neoproterozoic included changes in oceanic redox conditions, the configuration of continents and climate, extreme ice ages (Sturtian and Marinoan), and the rise of complex life forms. A much-debated topic in geobiology concerns the influence of atmospheric oxygenation on Earth and the origin and diversification of animal lineages, with the most widely popularized hypotheses relying on causal links between oxygen levels and the rise of animals. The vast majority of extant animals use aerobic metabolism for growth and homeostasis; hence, the binding and transportation of oxygen represent a vital physiological task. Considering the blood pigment hemocyanin (Hc) is present in sponges and ctenophores, and likely to be present in the common ancestor of animals, we investigated the evolution and date of Hc emergence using bioinformatics approaches on both transcriptomic and genomic data. Bayesian molecular dating suggested that the ancestral animal Hc gene arose approximately 881 Ma during the Tonian Period (1000-720 Ma), prior to the extreme glaciation events of the Cryogenian Period (720-635 Ma). This result is corroborated by a recently discovered fossil of a putative sponge ~890 Ma and modern molecular dating for the origin of metazoans of ~1,000-650 Ma (but does contradict previous inferences regarding the origin of Hc ~700-600 Ma). Our data reveal that crown-group animals already possessed hemocyanin-like blood pigments, which may have enhanced the oxygen-carrying capacity of these animals in hypoxic environments at that time or acted in the transport of hormones, detoxification of heavy metals, and immunity pathways.

Arrested in Glass: Actin within Sophisticated Architectures of Biosilica in Sponges

Actin is a fundamental member of an ancient superfamily of structural intracellular proteins and plays a crucial role in cytoskeleton dynamics, ciliogenesis, phagocytosis, and force generation in both prokaryotes and eukaryotes. It is shown that actin has another function in metazoans: patterning biosilica deposition, a role that has spanned over 500 million years. Species of glass sponges (Hexactinellida) and demosponges (Demospongiae), representatives of the first metazoans, with a broad diversity of skeletal structures with hierarchical architecture unchanged since the late Precambrian, are studied. By etching their skeletons, organic templates dominated by individual F-actin filaments, including branched fibers and the longest, thickest actin fiber bundles ever reported, are isolated. It is proposed that these actin-rich filaments are not the primary site of biosilicification, but this highly sophisticated and multi-scale form of biomineralization in metazoans is ptterned.

A Missing Link between Retrotransposons and Retroviruses

The origin and deep evolution of retroviruses remain largely unclear. It has been proposed that retroviruses might have originated from a Ty3/Gypsy retrotransposon, but all known Ty3/Gypsy retrotransposons are only distantly related to retroviruses. Retroviruses and some plant Athila/Tat elements (within Ty3/Gypsy retrotransposons) independently evolved a dual RNase H domain and an env/env-like gene. Here, we reported the discovery of a novel lineage of retrotransposons, designated Odin retrotransposons, in the genomes of eight sea anemones (order Actinaria) within the Cnidaria phylum. Odin retrotransposons exhibited unique genome features, encoding a dual RNase H domain (like retroviruses) but no env gene (like most Ty3/Gypsy retrotransposons). Phylogenetic analyses based on reverse transcriptase showed that Odin retrotransposons formed a sister group to lokiretroviruses, and lokiretroviruses and Odin retrotransposons together were sister to canonical retroviruses. Moreover, phylogenetic analyses based on RNase H and integrase also supported the hypothesis that Odin retrotransposons were sisters to lokiretroviruses. Lokiretroviruses and canonical retroviruses did not form a monophyletic group, indicating that lokiretroviruses and canonical retroviruses might represent two distinct virus families. Taken together, the discovery of Odin retrotransposons narrowed down the evolutionary gaps between retrotransposons and canonical retroviruses and lokiretroviruses.

Integrated phylogenomics and fossil data illuminate the evolution of beetles

Beetles constitute the most biodiverse animal order with over 380 000 described species and possibly several million more yet unnamed. Recent phylogenomic studies have arrived at considerably incongruent topologies and widely varying estimates of divergence dates for major beetle clades. Here, we use a dataset of 68 single-copy nuclear protein-coding (NPC) genes sampling 129 out of the 193 recognized extant families as well as the first comprehensive set of fully justified fossil calibrations to recover a refined timescale of beetle evolution. Using phylogenetic methods that counter the effects of compositional and rate heterogeneity, we recover a topology congruent with morphological studies, which we use, combined with other recent phylogenomic studies, to propose several formal changes in the classification of Coleoptera: Scirtiformia and Scirtoidea sensu nov., Clambiformia ser. nov. and Clamboidea sensu nov., Rhinorhipiformia ser. nov., Byrrhoidea sensu nov., Dryopoidea stat. res., Nosodendriformia ser. nov. and Staphyliniformia sensu nov., and Erotyloidea stat. nov., Nitiduloidea stat. nov. and Cucujoidea sensu nov., alongside changes below the superfamily level. Our divergence time analyses recovered a late Carboniferous origin of Coleoptera, a late Palaeozoic origin of all modern beetle suborders and a Triassic-Jurassic origin of most extant families, while fundamental divergences within beetle phylogeny did not coincide with the hypothesis of a Cretaceous Terrestrial Revolution.

Testing Phylogenetic Stability with Variable Taxon Sampling

Over the past three decades, computational capabilities have grown at such a rapid rate that they have given rise to many computationally heavy science fields such as phylogenomics. As increasingly more genomes are sequenced in the three domains of life, larger and more species-complete phylogenetic tree reconstructions are leading to a better understanding of the tree of life and the evolutionary histories in deep times. However, these large datasets pose unique challenges from a modeling and computational perspective: accurately describing the evolutionary process of thousands of species is still beyond the capability of current models, while the computational burden limits our ability to test multiple hypotheses. Thus, it is common practice to reduce the size of a dataset by selecting species to represent a clade (taxon sampling). Unfortunately, this process is subjective, and comparisons of large tree of life studies show that choice and number of species used in a dataset can alter the topology obtained. Thus, taxon sampling is, in itself, a process that needs to be fully investigated to determine its effect on phylogenetic stability. Here, we present the theory and practical application of an automated pipeline that can be easily implemented to explore the effect of taxon sampling on phylogenetic reconstructions. The application of this approach was recently discussed in a study of Terrabacteria and shows its power in investigating the accuracy of deep nodes of a phylogeny.

Rooting the Animal Tree of Life

Identifying our most distant animal relatives has emerged as one of the most challenging problems in phylogenetics. This debate has major implications for our understanding of the origin of multicellular animals and of the earliest events in animal evolution, including the origin of the nervous system. Some analyses identify sponges as our most distant animal relatives (Porifera-sister hypothesis), and others identify comb jellies (Ctenophora-sister hypothesis). These analyses vary in many respects, making it difficult to interpret previous tests of these hypotheses. To gain insight into why different studies yield different results, an important next step in the ongoing debate, we systematically test these hypotheses by synthesizing 15 previous phylogenomic studies and performing new standardized analyses under consistent conditions with additional models. We find that Ctenophora-sister is recovered across the full range of examined conditions, and Porifera-sister is recovered in some analyses under narrow conditions when most outgroups are excluded and site-heterogeneous CAT models are used. We additionally find that the number of categories in site-heterogeneous models is sufficient to explain the Porifera-sister results. Furthermore, our cross-validation analyses show CAT models that recover Porifera-sister have hundreds of additional categories and fail to fit significantly better than site-heterogenuous models with far fewer categories. Systematic and standardized testing of diverse phylogenetic models suggests that we should be skeptical of Porifera-sister results both because they are recovered under such narrow conditions and because the models in these conditions fit the data no better than other models that recover Ctenophora-sister.

Cambrian comb jellies from Utah illuminate the early evolution of nervous and sensory systems in ctenophores

Ctenophores are a group of predatory macroinvertebrates whose controversial phylogenetic position has prompted several competing hypotheses regarding the evolution of animal organ systems. Although ctenophores date back at least to the Cambrian, they have a poor fossil record due to their gelatinous bodies. Here, we describe two ctenophore species from the Cambrian of Utah, which illuminate the early evolution of nervous and sensory features in the phylum. Thalassostaphylos elegans has 16 comb rows, an oral skirt, and an apical organ with polar fields. Ctenorhabdotus campanelliformis has 24 comb rows, an oral skirt, an apical organ enclosed by a capsule and neurological tissues preserved as carbonaceous films. These are concentrated around the apical organ and ciliated furrows, which connect to a circumoral nerve ring via longitudinal axons. C. campanelliformis deviates from the neuroanatomy of living ctenophores and demonstrates a substantial complexity in the nervous system of Cambrian ctenophores.

•

Two species of rare fossil ctenophores are described from the Cambrian of Utah

•

Fossil ctenophores preserve remains of nervous tissue and sensory structures

•

Neurological structures include an oral nerve ring and giant longitudinal axons

•

Cambrian ctenophores had a more complex neuroanatomy than living species

A late-surviving stem-ctenophore from the Late Devonian of Miguasha (Canada)

Like other soft-bodied organisms, ctenophores (comb jellies) produce fossils only under exceptional taphonomic conditions. Here, we present the first record of a Late Devonian ctenophore from the Escuminac Formation from Miguasha in eastern Canada. Based on the 18-fold symmetry of this disc-shaped fossil, we assign it to the total-group Ctenophora. Our phylogenetic analyses suggest that the new taxon Daihuoides jakobvintheri gen. et sp. nov. falls near Cambrian stem ctenophores such as 'dinomischids' and 'scleroctenophorans'. Accordingly, Daihuoides is a Lazarus-taxon, which post-dates its older relatives by over 140 million years, and overlaps temporally with modern ctenophores, whose oldest representatives are known from the Early Devonian. Our analyses also indicate that the fossil record of ctenophores does not provide strong evidence for or against the phylogenomic hypothesis that ctenophores are sister to all other metazoans.

Evolution of a key enzyme of aerobic metabolism reveals Proterozoic functional subunit duplication events and an ancient origin of animals

The biological toolkits for aerobic respiration were critical for the rise and diversification of early animals. Aerobic life forms generate ATP through the oxidation of organic molecules in a process known as Krebs' Cycle, where the enzyme isocitrate dehydrogenase (IDH) regulates the cycle's turnover rate. Evolutionary reconstructions and molecular dating of proteins related to oxidative metabolism, such as IDH, can therefore provide an estimate of when the diversification of major taxa occurred, and their coevolution with the oxidative state of oceans and atmosphere. To establish the evolutionary history and divergence time of NAD-dependent IDH, we examined transcriptomic data from 195 eukaryotes (mostly animals). We demonstrate that two duplication events occurred in the evolutionary history of NAD-IDH, one in the ancestor of eukaryotes approximately at 1967 Ma, and another at 1629 Ma, both in the Paleoproterozoic Era. Moreover, NAD-IDH regulatory subunits β and γ are exclusive to metazoans, arising in the Mesoproterozoic. Our results therefore support the concept of an ''earlier-than-Tonian'' diversification of eukaryotes and the pre-Cryogenian emergence of a metazoan IDH enzyme.

Rapid radiation and rampant reticulation: Phylogenomics of South American Liolaemus lizards

Understanding the factors that cause heterogeneity among gene trees can increase the accuracy of species trees. Discordant signals across the genome are commonly produced by incomplete lineage sorting (ILS) and introgression, which in turn can result in reticulate evolution. Species tree inference using the multispecies coalescent is designed to deal with ILS and is robust to low levels of introgression, but extensive introgression violates the fundamental assumption that relationships are strictly bifurcating. In this study, we explore the phylogenomics of the iconic Liolaemus subgenus of South American lizards, a group of over 100 species mostly distributed in and around the Andes mountains. Using mitochondrial DNA (mtDNA) and genome-wide restriction-site associated DNA sequencing (RADseq; nDNA hereafter), we inferred a time-calibrated mtDNA gene tree, nDNA species trees, and phylogenetic networks. We found high levels of discordance between mtDNA and nDNA, which we attribute in part to extensive ILS resulting from rapid diversification. These data also reveal extensive and deep introgression, which combined with rapid diversification, explain the high level of phylogenetic discordance. We discuss these findings in the context of Andean orogeny and glacial cycles that fragmented, expanded, and contracted species distributions. Finally, we use the new phylogeny to resolve long-standing taxonomic issues in one of the most studied lizard groups in the New World.

Ionotropic Receptors as a Driving Force behind Human Synapse Establishment

The origin of nervous systems is a main theme in biology and its mechanisms are largely underlied by synaptic neurotransmission. One problem to explain synapse establishment is that synaptic orthologs are present in multiple aneural organisms. We questioned how the interactions among these elements evolved and to what extent it relates to our understanding of the nervous systems complexity. We identified the human neurotransmission gene network based on genes present in GABAergic, glutamatergic, serotonergic, dopaminergic, and cholinergic systems. The network comprises 321 human genes, 83 of which act exclusively in the nervous system. We reconstructed the evolutionary scenario of synapse emergence by looking for synaptic orthologs in 476 eukaryotes. The Human–Cnidaria common ancestor displayed a massive emergence of neuroexclusive genes, mainly ionotropic receptors, which might have been crucial to the evolution of synapses. Very few synaptic genes had their origin after the Human–Cnidaria common ancestor. We also identified a higher abundance of synaptic proteins in vertebrates, which suggests an increase in the synaptic network complexity of those organisms.

The IS630/Tc1/mariner transposons in three ctenophore genomes

Transposable elements (TEs) exert a significant effect on the structure and functioning of the genomes and also serve as a source of the new genes. The study of the TE diversity and evolution in different taxa is indispensable for the fundamental understanding of their roles in the genomes. IS630/Tc1/mariner (ITm) transposable elements represent the most prevalent and diverse group of DNA transposons. In this work, we studied the diversity, evolutionary dynamics and the phylogenetic relationships of the ITm transposons found in three ctenophore species: Mnemiopsis leidyi, Pleurobrachia bachei, Beroe ovata. We identified 29 ITm transposons, seven of which possess the terminal inverted repeats (TIRs) and an intact transposase, and, thus, are, presumably, active. Four other ITm transposons have the features of domesticated TEs. According to the results of the phylogenetic analysis, the ITm transposons of the ctenophores represent five groups - MLE/DD34D, TLE/DD34-38E, mosquito/DD37E, Visiror/DD41D and pogo/DDxD. Pogo/DDxD superfamily turnes out to be the most diverse and prevalent, since it accounts for more than 40% of the TEs identified. The data obtained in this research will fill the gap of knowledge of the diversity and evolution of the ITm transposons in the multicellular genomes and will lay the ground for the study of the TE effects on the evolution of the ctenophores.

Whole-genome microsynteny-based phylogeny of angiosperms

Plant genomes vary greatly in size, organization, and architecture. Such structural differences may be highly relevant for inference of genome evolution dynamics and phylogeny. Indeed, microsynteny-the conservation of local gene content and order-is recognized as a valuable source of phylogenetic information, but its use for the inference of large phylogenies has been limited. Here, by combining synteny network analysis, matrix representation, and maximum likelihood phylogenetic inference, we provide a way to reconstruct phylogenies based on microsynteny information. Both simulations and use of empirical data sets show our method to be accurate, consistent, and widely applicable. As an example, we focus on the analysis of a large-scale whole-genome data set for angiosperms, including more than 120 available high-quality genomes, representing more than 50 different plant families and 30 orders. Our 'microsynteny-based' tree is largely congruent with phylogenies proposed based on more traditional sequence alignment-based methods and current phylogenetic classifications but differs for some long-contested and controversial relationships. For instance, our synteny-based tree finds Vitales as early diverging eudicots, Saxifragales within superasterids, and magnoliids as sister to monocots. We discuss how synteny-based phylogenetic inference can complement traditional methods and could provide additional insights into some long-standing controversial phylogenetic relationships.

A single neuron subset governs a single coactive neuron circuit in Hydra vulgaris, representing a possible ancestral feature of neural evolution

The last common ancestor of Bilateria and Cnidaria is believed to be one of the first animals to develop a nervous system over 500 million years ago. Many of the genes involved in the neural function of the advanced nervous system in Bilateria are well conserved in Cnidaria. Thus, the cnidarian Hydra vulgaris is a good model organism for the study of the putative primitive nervous system in its last common ancestor. The diffuse nervous system of Hydra consists of several peptidergic neuron subsets. However, the specific functions of these subsets remain unclear. Using calcium imaging, here we show that the neuron subsets that express neuropeptide, Hym-176, function as motor circuits to evoke longitudinal contraction. We found that all neurons in a subset defined by the Hym-176 gene (Hym-176A) or its paralogs (Hym-176B) expression are excited simultaneously, followed by longitudinal contraction. This indicates not only that these neuron subsets have a motor function but also that a single molecularly defined neuron subset forms a single coactive circuit. This is in contrast with the bilaterian nervous system, where a single molecularly defined neuron subset harbors multiple coactive circuits, showing a mixture of neurons firing with different timings. Furthermore, we found that the two motor circuits, one expressing Hym-176B in the body column and the other expressing Hym-176A in the foot, are coordinately regulated to exert region-specific contraction. Our results demonstrate that one neuron subset is likely to form a monofunctional circuit as a minimum functional unit to build a more complex behavior in Hydra. This simple feature (one subset, one circuit, one function) found in Hydra may represent the simple ancestral condition of neural evolution.

Scalable Empirical Mixture Models That Account for Across-Site Compositional Heterogeneity

Biochemical demands constrain the range of amino acids acceptable at specific sites resulting in across-site compositional heterogeneity of the amino acid replacement process. Phylogenetic models that disregard this heterogeneity are prone to systematic errors, which can lead to severe long-branch attraction artifacts. State-of-the-art models accounting for across-site compositional heterogeneity include the CAT model, which is computationally expensive, and empirical distribution mixture models estimated via maximum likelihood (C10–C60 models). Here, we present a new, scalable method EDCluster for finding empirical distribution mixture models involving a simple cluster analysis. The cluster analysis utilizes specific coordinate transformations which allow the detection of specialized amino acid distributions either from curated databases or from the alignment at hand. We apply EDCluster to the HOGENOM and HSSP databases in order to provide universal distribution mixture (UDM) models comprising up to 4,096 components. Detailed analyses of the UDM models demonstrate the removal of various long-branch attraction artifacts and improved performance compared with the C10–C60 models. Ready-to-use implementations of the UDM models are provided for three established software packages (IQ-TREE, Phylobayes, and RevBayes).

Six-state Amino Acid Recoding is not an Effective Strategy to Offset Compositional Heterogeneity and Saturation in Phylogenetic Analyses

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.]

Relative Model Fit Does Not Predict Topological Accuracy in Single-Gene Protein Phylogenetics

It is regarded as best practice in phylogenetic reconstruction to perform relative model selection to determine an appropriate evolutionary model for the data. This procedure ranks a set of candidate models according to their goodness of fit to the data, commonly using an information theoretic criterion. Users then specify the best-ranking model for inference. Although it is often assumed that better-fitting models translate to increase accuracy, recent studies have shown that the specific model employed may not substantially affect inferences. We examine whether there is a systematic relationship between relative model fit and topological inference accuracy in protein phylogenetics, using simulations and real sequences. Simulations employed site-heterogeneous mechanistic codon models that are distinct from protein-level phylogenetic inference models, allowing us to investigate how protein models performs when they are misspecified to the data, as will be the case for any real sequence analysis. We broadly find that phylogenies inferred across models with vastly different fits to the data produce highly consistent topologies. We additionally find that all models infer similar proportions of false-positive splits, raising the possibility that all available models of protein evolution are similarly misspecified. Moreover, we find that the parameter-rich GTR (general time reversible) model, whose amino acid exchangeabilities are free parameters, performs similarly to models with fixed exchangeabilities, although the inference precision associated with GTR models was not examined. We conclude that, although relative model selection may not hinder phylogenetic analysis on protein data, it may not offer specific predictable improvements and is not a reliable proxy for accuracy.

Inferring the deep past from molecular data

There is an expectation that analyses of molecular sequences might be able to distinguish between alternative hypotheses for ancient relationships, but the phylogenetic methods used and types of data analyzed are of critical importance in any attempt to recover historical signal. Here, we discuss some common issues that can influence the topology of trees obtained when using overly simple models to analyze molecular data that often display complicated patterns of sequence heterogeneity. To illustrate our discussion, we have used three examples of inferred relationships which have changed radically as models and methods of analysis have improved. In two of these examples, the sister-group relationship between thermophilic Thermus and mesophilic Deinococcus, and the position of long-branch Microsporidia among eukaryotes, we show that recovering what is now generally considered to be the correct tree is critically dependent on the fit between model and data. In the third example, the position of eukaryotes in the tree of life, the hypothesis that is currently supported by the best available methods is fundamentally different from the classical view of relationships between major cellular domains. Since heterogeneity appears to be pervasive and varied among all molecular sequence data, and even the best available models can still struggle to deal with some problems, the issues we discuss are generally relevant to phylogenetic analyses. It remains essential to maintain a critical attitude to all trees as hypotheses of relationship that may change with more data and better methods.

The Roles of Protein Structure, Taxon Sampling, and Model Complexity in Phylogenomics: A Case Study Focused on Early Animal Divergences

Despite the long history of using protein sequences to infer the tree of life, the potential for different parts of protein structures to retain historical signal remains unclear. We propose that it might be possible to improve analyses of phylogenomic datasets by incorporating information about protein structure. We test this idea using the position of the root of Metazoa (animals) as a model system. We examined the distribution of “strongly decisive” sites (alignment positions that support a specific tree topology) in a dataset comprising >1500 proteins and almost 100 taxa. The proportion of each class of strongly decisive sites in different structural environments was very sensitive to the model used to analyze the data when a limited number of taxa were used but they were stable when taxa were added. As long as enough taxa were analyzed, sites in all structural environments supported the same topology regardless of whether standard tree searches or decisive sites were used to select the optimal tree. However, the use of decisive sites revealed a difference between the support for minority topologies for sites in different structural environments: buried sites and sites in sheet and coil environments exhibited equal support for the minority topologies, whereas solvent-exposed and helix sites had unequal numbers of sites, supporting the minority topologies. This suggests that the relatively slowly evolving buried, sheet, and coil sites are giving an accurate picture of the true species tree and the amount of conflict among gene trees. Taken as a whole, this study indicates that phylogenetic analyses using sites in different structural environments can yield different topologies for the deepest branches in the animal tree of life and that analyzing larger numbers of taxa eliminates this conflict. More broadly, our results highlight the desirability of incorporating information about protein structure into phylogenomic analyses.

The ABCs of the atypical Fam20 secretory pathway kinases

The study of extracellular phosphorylation was initiated in late 19th century when the secreted milk protein, casein, and egg-yolk protein, phosvitin, were shown to be phosphorylated. However, it took more than a century to identify Fam20C, which phosphorylates both casein and phosvitin under physiological conditions. This kinase, along with its family members Fam20A and Fam20B, defined a new family with altered amino acid sequences highly atypical from the canonical 540 kinases comprising the kinome. Fam20B is a glycan kinase that phosphorylates xylose residues and triggers peptidoglycan biosynthesis, a role conserved from sponges to human. The protein kinase, Fam20C, conserved from nematodes to humans, phosphorylates well over 100 substrates in the secretory pathway with overall functions postulated to encompass endoplasmic reticulum homeostasis, nutrition, cardiac function, coagulation, and biomineralization. The preferred phosphorylation motif of Fam20C is SxE/pS, and structural studies revealed that related member Fam20A allosterically activates Fam20C by forming a heterodimeric/tetrameric complex. Fam20A, a pseudokinase, is observed only in vertebrates. Loss-of-function genetic alterations in the Fam20 family lead to human diseases such as amelogenesis imperfecta, nephrocalcinosis, lethal and nonlethal forms of Raine syndrome with major skeletal defects, and altered phosphate homeostasis. Together, these three members of the Fam20 family modulate a diverse network of secretory pathway components playing crucial roles in health and disease. The overarching theme of this review is to highlight the progress that has been made in the emerging field of extracellular phosphorylation and the key roles secretory pathway kinases play in an ever-expanding number of cellular processes.

Evidence for sponges as sister to all other animals from partitioned phylogenomics with mixture models and recoding

Resolving the relationships between the major lineages in the animal tree of life is necessary to understand the origin and evolution of key animal traits. Sponges, characterized by their simple body plan, were traditionally considered the sister group of all other animal lineages, implying a gradual increase in animal complexity from unicellularity to complex multicellularity. However, the availability of genomic data has sparked tremendous controversy as some phylogenomic studies support comb jellies taking this position, requiring secondary loss or independent origins of complex traits. Here we show that incorporating site-heterogeneous mixture models and recoding into partitioned phylogenomics alleviates systematic errors that hamper commonly-applied phylogenetic models. Testing on real datasets, we show a great improvement in model-fit that attenuates branching artefacts induced by systematic error. We reanalyse key datasets and show that partitioned phylogenomics does not support comb jellies as sister to other animals at either the supermatrix or partition-specific level.

Resolving the Early Divergence Pattern of Teleost Fish Using Genome-Scale Data

Regarding the phylogenetic relationship of the three primary groups of teleost fishes, Osteoglossomorpha (bonytongues and others), Elopomorpha (eels and relatives), Clupeocephala (the remaining teleost fish), early morphological studies hypothesized the first divergence of Osteoglossomorpha, whereas the recent prevailing view is the first divergence of Elopomorpha. Molecular studies supported all the possible relationships of the three primary groups. This study analyzed genome-scale data from four previous studies: 1) 412 genes from 12 species, 2) 772 genes from 15 species, 3) 1,062 genes from 30 species, and 4) 491 UCE loci from 27 species. The effects of the species, loci, and models used on the constructed tree topologies were investigated. In the analyses of the data sets (1)–(3), although the first divergence of Clupeocephala that left the other two groups in a sister relationship was supported by concatenated sequences and gene trees of all the species and genes, the first divergence of Elopomorpha among the three groups was supported using species and/or genes with low divergence of sequence and amino-acid frequencies. This result corresponded to that of the UCE data set (4), whose sequence divergence was low, which supported the first divergence of Elopomorpha with high statistical significance. The increase in accuracy of the phylogenetic construction by using species and genes with low sequence divergence was predicted by a phylogenetic informativeness approach and confirmed by computer simulation. These results supported that Elopomorpha was the first basal group of teleost fish to have diverged, consistent with the prevailing view of recent morphological studies.

The origin of animals: an ancestral reconstruction of the unicellular-to-multicellular transition

How animals evolved from a single-celled ancestor, transitioning from a unicellular lifestyle to a coordinated multicellular entity, remains a fascinating question. Key events in this transition involved the emergence of processes related to cell adhesion, cell–cell communication and gene regulation. To understand how these capacities evolved, we need to reconstruct the features of both the last common multicellular ancestor of animals and the last unicellular ancestor of animals. In this review, we summarize recent advances in the characterization of these ancestors, inferred by comparative genomic analyses between the earliest branching animals and those radiating later, and between animals and their closest unicellular relatives. We also provide an updated hypothesis regarding the transition to animal multicellularity, which was likely gradual and involved the use of gene regulatory mechanisms in the emergence of early developmental and morphogenetic plans. Finally, we discuss some new avenues of research that will complement these studies in the coming years.