Improved Modeling of Compositional Heterogeneity Supports Sponges as Sister to All Other Animals

Piezo proteins, mechano reception and behaviour in Hydra

Cells are constantly exposed to mechanical forces. These stimuli must be detected and transduced to maintain homeostasis. Due that reason, mechano-sensory systems (MS) are present in all the organisms to generate appropriate responses. Piezo proteins are a recently discovered family of mechano-gated ion channels that responds to mechanical changes of the plasma membrane, allowing the influx of cations, mainly Ca2+. Piezo MS channels are widely represented in Metazoa, acting in several physiological systems. Hydra sp. is a freshwater member of the phylum Cnidaria which is assumed as the sister group of Bilateria. Despite the existence of Piezo channels in Hydra is known, their physiological roles remain unknown. In this work we delve into the physiological relevance of MS Piezo in responses associated to mechanical stimuli in Hydra sp. We analysed the effects of Jedi1 (a specific agonist of Piezo1) on the contractile behaviours, and cnidocyst discharge, and compared them with responses caused by natural stimuli. The results show that the activation of Piezo channels increases the contractile behaviour, stimulating the occurrence of contraction burst, a sudden kind of retraction observed in response to touch and osmotic stress. Cnidocyst discharge was also induced by Jedi1, resembling the response caused by the contact of the prey. The effects of both Jedi1 and natural stimuli were avoided in the presence of the inhibitor of MS channels GdCl3. The bioinformatic analysis shows that the protein predicted by Hydra genome has the characteristic motifs of Piezo proteins. These results are consistent with the existence of Piezo channels in Hydra, unveiling their physiological roles. Because of the relevance of Piezo channels in several pathological conditions and the high level of conservation in metazoans, Hydra could provide a significative experimental model to assay biological and pharmacological issues.

The Chromosome-level Genome of the Ctenophore Mnemiopsis leidyi A. Agassiz, 1865 Reveals a Unique Immune Gene Repertoire

Ctenophora are basal marine metazoans, the sister group of all other animals. Mnemiopsis leidyi is one of the most successful invasive species worldwide with intense ecological and evolutionary research interest. Here, we generated a chromosome-level genome assembly of M. leidyi with a focus on its immune gene repertoire. The genome was 247.97 Mb, with N50 16.84 Mb, and 84.7% completeness. Its karyotype was 13 chromosomes. In this genome and that of two other ctenophores, Bolinopsis microptera and Hormiphora californensis, we detected a high number of protein domains related to potential immune receptors. Among those, proteins containing Toll/interleukin-1 (TIR2) domain, NACHT domain, Scavenger Receptor Cystein-Rich (SRCR) domain, or C-type Lectin domain (CTLD) were abundant and presented unique domain architectures in M. leidyi. M. leidyi seems to lack bona fide Toll-like Receptors, but it does possess a repertoire of 15 TIR2 domain-containing genes. Besides, we detected a bona fide NOD-like receptor and 38 NACHT domain-containing genes. In order to verify the function of those domain-containing genes, we exposed M. leidyi to the pathogen Vibrio coralliilyticus. Among the differentially expressed genes, we identified potential immune receptors, including four TIR2 domain-containing genes, all of which were upregulated in response to pathogen exposure. To conclude, many common immune receptor domains, highly conserved across metazoans, are already present in Ctenophora. These domains have large expansions and unique architectures in M. leidyi, findings consistent with the basal evolutionary position of this group, but still might have conserved functions in immunity and host-microbe interaction.

CAT-Posterior Mean Site Frequencies Improves Phylogenetic Modeling Under Maximum Likelihood and Resolves Tardigrada as the Sister of Arthropoda Plus Onychophora

Tardigrada, the water bears, are microscopic animals with walking appendages that are members of Ecdysozoa, the clade of molting animals that also includes Nematoda (round worms), Nematomorpha (horsehair worms), Priapulida (penis worms), Kinorhyncha (mud dragons), Loricifera (loricated animals), Arthropoda (insects, spiders, centipedes, crustaceans, and their allies), and Onychophora (velvet worms). The phylogenetic relationships within Ecdysozoa are still unclear, with analyses of molecular and morphological data yielding incongruent results. Accounting for across-site compositional heterogeneity using mixture models that partition sites in frequency categories, CATegories (CAT)-based models, has been shown to improve fit in Bayesian analyses. However, CAT-based models such as CAT-Poisson or CAT-GTR (where CAT is combined with a General Time Reversible matrix to account for replacement rate heterogeneity) have proven difficult to implement in maximum likelihood. Here, we use CAT-posterior mean site frequencies (CAT-PMSF), a new method to export dataset-specific mixture models (CAT-Poisson and CAT-GTR) parameterized using Bayesian methods to maximum likelihood software. We developed new maximum likelihood-based model adequacy tests using parametric bootstrap and show that CAT-PMSF describes across-site compositional heterogeneity better than other across-site compositionally heterogeneous models currently implemented in maximum likelihood software. CAT-PMSF suggests that tardigrades are members of Panarthropoda, a lineage also including Arthropoda and Onychophora. Within Panarthropoda, our results favor Tardigrada as sister to Onychophora plus Arthropoda (the Lobopodia hypothesis). Our results illustrate the power of CAT-PMSF to model across-site compositionally heterogeneous datasets in the maximum likelihood framework and clarify the relationships between the Tardigrada and the Ecdysozoa.

Phylogenomics resolves a 100-year-old debate regarding the evolutionary history of caddisflies (Insecta: Trichoptera)

Trichoptera (caddisfly) phylogeny provides an interesting example of aquatic insect evolution, with rich ecological diversification, especially for underwater architecture. Trichoptera provide numerous critical ecosystem services and are also one of the most important groups of aquatic insects for assessing water quality. The phylogenetic relationships of Trichoptera have been debated for nearly a century. In particular, the phylogenetic position of the "cocoon-makers" within Trichoptera has long been contested. Here, we designed a universal single-copy orthologue and sets of ultraconserved element markers specific for Trichoptera for the first time. Simultaneously, we reconstructed the phylogenetic relationship of Trichoptera based on genome data from 111 species, representing 29 families and 71 genera. Our phylogenetic inference clarifies the probable phylogenetic relationships of "cocoon-makers" within Integripalpia. Hydroptilidae is considered as the basal lineage within Integripalpia, and the families Glossosomatidae, Hydrobiosidae, and Rhyacophilidae formed a monophyletic clade, sister to the integripalpian subterorder Phryganides. The resulting divergence time and ancestral state reconstruction suggest that the most recent common ancestor of Trichoptera appeared in the early Permian and that diversification was strongly correlated with habitat change.

BAD2matrix: Phylogenomic matrix concatenation, indel coding, and more

Premise

Common steps in phylogenomic matrix production include biological sequence concatenation, morphological data concatenation, insertion/deletion (indel) coding, gene content (presence/absence) coding, removing uninformative characters for parsimony analysis, recording with reduced amino acid alphabets, and occupancy filtering. Existing software does not accomplish these tasks on a phylogenomic scale using a single program.

Methods and Results

BAD2matrix is a Python script that performs the above-mentioned steps in phylogenomic matrix construction for DNA or amino acid sequences as well as morphological data. The script works in UNIX-like environments (e.g., LINUX, MacOS, Windows Subsystem for LINUX).

Conclusions

BAD2matrix helps simplify phylogenomic pipelines and can be downloaded from https://github.com/dpl10/BAD2matrix/tree/master under a GNU General Public License v2.

Sequential Bayesian Phylogenetic Inference

The ideal approach to Bayesian phylogenetic inference is to estimate all parameters of interest jointly in a single hierarchical model. However, this is often not feasible in practice due to the high computational cost. Instead, phylogenetic pipelines generally consist of sequential analyses, whereby a single point estimate from a given analysis is used as input for the next analysis (e.g., a single multiple sequence alignment is used to estimate a gene tree). In this framework, uncertainty is not propagated from step to step, which can lead to inaccurate or spuriously confident results. Here, we formally develop and test a sequential inference approach for Bayesian phylogenetic inference, which uses importance sampling to generate observations for the next step of an analysis pipeline from the posterior distribution produced in the previous step. Our sequential inference approach presented here not only accounts for uncertainty between analysis steps but also allows for greater flexibility in software choice (and hence model availability) and can be computationally more efficient than the traditional joint inference approach when multiple models are being tested. We show that our sequential inference approach is identical in practice to the joint inference approach only if sufficient information in the data is present (a narrow posterior distribution) and/or sufficiently many important samples are used. Conversely, we show that the common practice of using a single point estimate can be biased, for example, a single phylogeny estimate can transform an unrooted phylogeny into a time-calibrated phylogeny. We demonstrate the theory of sequential Bayesian inference using both a toy example and an empirical case study of divergence-time estimation in insects using a relaxed clock model from transcriptome data. In the empirical example, we estimate 3 posterior distributions of branch lengths from the same data (DNA character matrix with a GTR+Γ+I substitution model, an amino acid data matrix with empirical substitution models, and an amino acid data matrix with the PhyloBayes CAT-GTR model). Finally, we apply 3 different node-calibration strategies and show that divergence time estimates are affected by both the data source and underlying substitution process to estimate branch lengths as well as the node-calibration strategies. Thus, our new sequential Bayesian phylogenetic inference provides the opportunity to efficiently test different approaches for divergence time estimation, including branch-length estimation from other software.

Acoelomorph flatworm monophyly is a long-branch attraction artefact obscuring a clade of Acoela and Xenoturbellida

Acoelomorpha is a broadly accepted clade of bilaterian animals made up of the fast-evolving, morphologically simple, mainly marine flatworm lineages Acoela and Nemertodermatida. Phylogenomic studies support Acoelomorpha's close relationship with the slowly evolving and similarly simplistic Xenoturbella, together forming the phylum Xenacoelomorpha. The phylogenetic placement of Xenacoelomorpha amongst bilaterians is controversial, with some studies supporting Xenacoelomorpha as the sister group to all other bilaterians, implying that their simplicity may be representative of early bilaterians. Others propose that this placement is an error resulting from the fast-evolving Acoelomorpha, and instead suggest that they are the degenerate sister group to Ambulacraria. Perhaps as a result of this debate, internal xenacoelomorph relationships have been somewhat overlooked at a phylogenomic scale. Here, I employ a highly targeted approach to detect and overcome possible phylogenomic error in the relationship between Xenoturbella and the fast-evolving acoelomorph flatworms. The results indicate that the subphylum Acoelomorpha is a long-branch attraction artefact obscuring a previously undiscovered clade comprising Xenoturbella and Acoela, which I name Xenacoela. The findings also suggest that Xenacoelomorpha is not the sister group to all other bilaterians. This study provides a template for future efforts aimed at discovering and correcting unrecognized long-branch attraction artefacts throughout the tree of life.

A taxon-rich and genome-scale phylogeny of Opisthokonta

Ancient divergences within Opisthokonta-a major lineage that includes organisms in the kingdoms Animalia, Fungi, and their unicellular relatives-remain contentious. To assess progress toward a genome-scale Opisthokonta phylogeny, we conducted the most taxon rich phylogenomic analysis using sets of genes inferred with different orthology inference methods and established the geological timeline of Opisthokonta diversification. We also conducted sensitivity analysis by subsampling genes or taxa from the full data matrix based on filtering criteria previously shown to improve phylogenomic inference. We found that approximately 85% of internal branches were congruent across data matrices and the approaches used. Notably, the use of different orthology inference methods was a substantial contributor to the observed incongruence: analyses using the same set of orthologs showed high congruence of 97% to 98%, whereas different sets of orthologs resulted in somewhat lower congruence (87% to 91%). Examination of unicellular Holozoa relationships suggests that the instability observed across varying gene sets may stem from weak phylogenetic signals. Our results provide a comprehensive Opisthokonta phylogenomic framework that will be useful for illuminating ancient evolutionary episodes concerning the origin and diversification of the 2 major eukaryotic kingdoms and emphasize the importance of investigating effects of orthology inference on phylogenetic analyses to resolve ancient divergences.

Discovering Fragile Clades and Causal Sequences in Phylogenomics by Evolutionary Sparse Learning

Phylogenomic analyses of long sequences, consisting of many genes and genomic segments, reconstruct organismal relationships with high statistical confidence. But, inferred relationships can be sensitive to excluding just a few sequences. Currently, there is no direct way to identify fragile relationships and the associated individual gene sequences in species. Here, we introduce novel metrics for gene-species sequence concordance and clade probability derived from evolutionary sparse learning models. We validated these metrics using fungi, plant, and animal phylogenomic datasets, highlighting the ability of the new metrics to pinpoint fragile clades and the sequences responsible. The new approach does not necessitate the investigation of alternative phylogenetic hypotheses, substitution models, or repeated data subset analyses. Our methodology offers a streamlined approach to evaluating major inferred clades and identifying sequences that may distort reconstructed phylogenies using large datasets.

Evolutionary origin of the nervous system from Ctenophora prospective

Nervous system is one of the key adaptations underlying the evolutionary success of the majority of animal groups. Ctenophores (or comb jellies) are gelatinous marine invertebrates that were probably the first lineage to diverge from the rest of animals. Due to the key phylogenetic position and multiple unique adaptations, the noncentralized nervous system of comb jellies has been in the center of the debate around the origin of the nervous system in the animal kingdom and whether it happened only once or twice. Here, we discuss the latest findings in ctenophore neuroscience and multiple challenges on the way to build a clear evolutionary picture of the origin of the nervous system.

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.

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.

Evolution of the ribbon-like organization of the Golgi apparatus in animal cells

The "ribbon," a structural arrangement in which Golgi stacks connect to each other, is considered to be restricted to vertebrate cells. Although ribbon disruption is linked to various human pathologies, its functional role in cellular processes remains unclear. In this study, we investigate the evolutionary origin of the Golgi ribbon. We observe a ribbon-like architecture in the cells of several metazoan taxa suggesting its early emergence in animal evolution predating the appearance of vertebrates. Supported by AlphaFold2 modeling, we propose that the evolution of Golgi reassembly and stacking protein (GRASP) binding by golgin tethers may have driven the joining of Golgi stacks resulting in the ribbon-like configuration. Additionally, we find that Golgi ribbon assembly is a shared developmental feature of deuterostomes, implying a role in embryogenesis. Overall, our study points to the functional significance of the Golgi ribbon beyond vertebrates and underscores the need for further investigations to unravel its elusive biological roles.

Rediscovering the unusual, solitary bryozoan Monobryozoon ambulans Remane, 1936: first molecular and new morphological data clarify its phylogenetic position

BACKGROUND

One of the most peculiar groups of the mostly colonial phylum Bryozoa is the taxon Monobryozoon, whose name already implies non-colonial members of the phylum. Its peculiarity and highly unusual lifestyle as a meiobenthic clade living on sand grains has fascinated many biologists. In particular its systematic relationship to other bryozoans remains a mystery. Despite numerous searches for M. ambulans in its type locality Helgoland, a locality with a long-lasting marine station and tradition of numerous courses and workshops, it has never been reencountered until today. Here we report the first observations of this almost mythical species, Monobryozoon ambulans.

RESULTS

For the first time since 1938, we present new modern, morphological analyses of this species as well as the first ever molecular data. Our detailed morphological analysis confirms most previous descriptions, but also ascertains the presence of special ambulatory polymorphic zooids. We consider these as bud anlagen that ultimately consecutively separate from the animal rendering it pseudo-colonial. The remaining morphological data show strong ties to alcyonidioidean ctenostome bryozoans. Our morphological data is in accordance with the phylogenomic analysis, which clusters it with species of Alcyonidium as a sister group to multiporate ctenostomes. Divergence time estimation and ancestral state reconstruction recover the solitary state of M. ambulans as a derived character that probably evolved in the Late Cretaceous. In this study, we also provide the entire mitogenome of M. ambulans, which-despite the momentary lack of comparable data-provides important data of a unique and rare species for comparative aspects in the future.

CONCLUSIONS

We were able to provide first sequence data and modern morphological data for the unique bryozoan, M. ambulans, which are both supporting an alcyonidioidean relationship within ctenostome bryozoans.

Ant backbone phylogeny resolved by modelling compositional heterogeneity among sites in genomic data

Ants are the most ubiquitous and ecologically dominant arthropods on Earth, and understanding their phylogeny is crucial for deciphering their character evolution, species diversification, and biogeography. Although recent genomic data have shown promise in clarifying intrafamilial relationships across the tree of ants, inconsistencies between molecular datasets have also emerged. Here I re-examine the most comprehensive published Sanger-sequencing and genome-scale datasets of ants using model comparison methods that model among-site compositional heterogeneity to understand the sources of conflict in phylogenetic studies. My results under the best-fitting model, selected on the basis of Bayesian cross-validation and posterior predictive model checking, identify contentious nodes in ant phylogeny whose resolution is modelling-dependent. I show that the Bayesian infinite mixture CAT model outperforms empirical finite mixture models (C20, C40 and C60) and that, under the best-fitting CAT-GTR + G4 model, the enigmatic Martialis heureka is sister to all ants except Leptanillinae, rejecting the more popular hypothesis supported under worse-fitting models, that place it as sister to Leptanillinae. These analyses resolve a lasting controversy in ant phylogeny and highlight the significance of model comparison and adequate modelling of among-site compositional heterogeneity in reconstructing the deep phylogeny of insects.

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.

DNA Sequences Are as Useful as Protein Sequences for Inferring Deep Phylogenies

Inference of deep phylogenies has almost exclusively used protein rather than DNA sequences based on the perception that protein sequences are less prone to homoplasy and saturation or to issues of compositional heterogeneity than DNA sequences. Here, we analyze a model of codon evolution under an idealized genetic code and demonstrate that those perceptions may be misconceptions. We conduct a simulation study to assess the utility of protein versus DNA sequences for inferring deep phylogenies, with protein-coding data generated under models of heterogeneous substitution processes across sites in the sequence and among lineages on the tree, and then analyzed using nucleotide, amino acid, and codon models. Analysis of DNA sequences under nucleotide-substitution models (possibly with the third codon positions excluded) recovered the correct tree at least as often as analysis of the corresponding protein sequences under modern amino acid models. We also applied the different data-analysis strategies to an empirical dataset to infer the metazoan phylogeny. Our results from both simulated and real data suggest that DNA sequences may be as useful as proteins for inferring deep phylogenies and should not be excluded from such analyses. Analysis of DNA data under nucleotide models has a major computational advantage over protein-data analysis, potentially making it feasible to use advanced models that account for among-site and among-lineage heterogeneity in the nucleotide-substitution process in inference of deep phylogenies.

Stepwise emergence of the neuronal gene expression program in early animal evolution

The assembly of the neuronal and other major cell type programs occurred early in animal evolution. We can reconstruct this process by studying non-bilaterians like placozoans. These small disc-shaped animals not only have nine morphologically described cell types and no neurons but also show coordinated behaviors triggered by peptide-secreting cells. We investigated possible neuronal affinities of these peptidergic cells using phylogenetics, chromatin profiling, and comparative single-cell genomics in four placozoans. We found conserved cell type expression programs across placozoans, including populations of transdifferentiating and cycling cells, suggestive of active cell type homeostasis. We also uncovered fourteen peptidergic cell types expressing neuronal-associated components like the pre-synaptic scaffold that derive from progenitor cells with neurogenesis signatures. In contrast, earlier-branching animals like sponges and ctenophores lacked this conserved expression. Our findings indicate that key neuronal developmental and effector gene modules evolved before the advent of cnidarian/bilaterian neurons in the context of paracrine cell signaling.

Data-specific substitution models improve protein-based phylogenetics

Calculating amino-acid substitution models that are specific for individual protein data sets is often difficult due to the computational burden of estimating large numbers of rate parameters. In this study, we tested the computational efficiency and accuracy of five methods used to estimate substitution models, namely Codeml, FastMG, IQ-TREE, P4 (maximum likelihood), and P4 (Bayesian inference). Data-specific substitution models were estimated from simulated alignments (with different lengths) that were generated from a known simulation model and simulation tree. Each of the resulting data-specific substitution models was used to calculate the maximum likelihood score of the simulation tree and simulated data that was used to calculate the model, and compared with the maximum likelihood scores of the known simulation model and simulation tree on the same simulated data. Additionally, the commonly-used empirical models, cpREV and WAG, were assessed similarly. Data-specific models performed better than the empirical models, which under-fitted the simulated alignments, had the highest difference to the simulation model maximum-likelihood score, clustered further from the simulation model in principal component analysis ordination, and inferred less accurate trees. Data-specific models and the simulation model shared statistically indistinguishable maximum-likelihood scores, indicating that the five methods were reasonably accurate at estimating substitution models by this measure. Nevertheless, tree statistics showed differences between optimal maximum likelihood trees. Unlike other model estimating methods, trees inferred using data-specific models generated with IQ-TREE and P4 (maximum likelihood) were not significantly different from the trees derived from the simulation model in each analysis, indicating that these two methods alone were the most accurate at estimating data-specific models. To show the benefits of using data-specific protein models several published data sets were reanalysed using IQ-TREE-estimated models. These newly estimated models were a better fit to the data than the empirical models that were used by the original authors, often inferred longer trees, and resulted in different tree topologies in more than half of the re-analysed data sets. The results of this study show that software availability and high computation burden are not limitations to generating better-fitting data-specific amino-acid substitution models for phylogenetic analyses.

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.

Morphological and transcriptional effects of crude oil and dispersant exposure on the marine sponge Cinachyrella alloclada

Marine sponges play important roles in benthic ecosystems. More than providing shelter and food to other species, they help maintain water quality by regulating nitrogen and ammonium levels in the water, and bioaccumulate heavy metals. This system, however, is particularly sensitive to sudden environmental changes including catastrophic pollution event such as oil spills. Hundreds of oil platforms are currently actively extracting oil and gas in the Gulf of Mexico. To test the vulnerability of the benthic ecosystems to oil spills, we utilized the Caribbean reef sponge, Cinachyrella alloclada, as a novel experimental indicator. We have exposed organisms to crude oil and oil dispersant for up to 24 h and measured resultant gene expression changes. Our findings indicate that 1-hour exposure to water accommodated fractions (WAF) was enough to elicit massive shifts in gene expression in sponges and host bacterial communities (8052 differentially expressed transcripts) with the up-regulation of stress related pathways, cancer related pathways, and cell integrity pathways. Genes that were upregulated included heat shock proteins, apoptosis, oncogenes (Rab/Ras, Src, CMYC), and several E3 ubiquitin ligases. 24-hour exposure of chemically enhanced WAF (CE-WAF) had the greatest impact to benthic communities, resulting in mostly downregulation of gene expression (4248 differentially expressed transcripts). Gene deregulation from 1-hour treatments follow this decreasing trend of toxicity: WAF > CE-WAF > Dispersant, while the 24-hour treatment showed a shift to CE-WAF > Dispersant > WAF in our experiments. Thus, this study supports the development of Cinachyrella alloclada as a research model organism and bioindicator species for Florida reefs and underscores the importance of developing more efficient and safer ways to remove oil in the event of a spill catastrophe.

The compact genome of the sponge Oopsacas minuta (Hexactinellida) is lacking key metazoan core genes

BACKGROUND

Explaining the emergence of the hallmarks of bilaterians is a central focus of evolutionary developmental biology-evodevo-and evolutionary genomics. For this purpose, we must both expand and also refine our knowledge of non-bilaterian genomes, especially by studying early branching animals, in particular those in the metazoan phylum Porifera.

RESULTS

We present a comprehensive analysis of the first whole genome of a glass sponge, Oopsacas minuta, a member of the Hexactinellida. Studying this class of sponge is evolutionary relevant because it differs from the three other Porifera classes in terms of development, tissue organization, ecology, and physiology. Although O. minuta does not exhibit drastic body simplifications, its genome is among the smallest of animal genomes sequenced so far, and surprisingly lacks several metazoan core genes (including Wnt and several key transcription factors). Our study also provides the complete genome of a symbiotic Archaea dominating the associated microbial community: a new Thaumarchaeota species.

CONCLUSIONS

The genome of the glass sponge O. minuta differs from all other available sponge genomes by its compactness and smaller number of encoded proteins. The unexpected loss of numerous genes previously considered ancestral and pivotal for metazoan morphogenetic processes most likely reflects the peculiar syncytial tissue organization in this group. Our work further documents the importance of convergence during animal evolution, with multiple convergent evolution of septate-like junctions, electrical-signaling and multiciliated cells in metazoans.

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.

Neurons that connect without synapses

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

Bacterial aerobic methane cycling by the marine sponge-associated microbiome

BACKGROUND

Methanotrophy by the sponge-hosted microbiome has been mainly reported in the ecological context of deep-sea hydrocarbon seep niches where methane is either produced geothermically or via anaerobic methanogenic archaea inhabiting the sulfate-depleted sediments. However, methane-oxidizing bacteria from the candidate phylum Binatota have recently been described and shown to be present in oxic shallow-water marine sponges, where sources of methane remain undescribed.

RESULTS

Here, using an integrative -omics approach, we provide evidence for sponge-hosted bacterial methane synthesis occurring in fully oxygenated shallow-water habitats. Specifically, we suggest methane generation occurs via at least two independent pathways involving methylamine and methylphosphonate transformations that, concomitantly to aerobic methane production, generate bioavailable nitrogen and phosphate, respectively. Methylphosphonate may be sourced from seawater continuously filtered by the sponge host. Methylamines may also be externally sourced or, alternatively, generated by a multi-step metabolic process where carnitine, derived from sponge cell debris, is transformed to methylamine by different sponge-hosted microbial lineages. Finally, methanotrophs specialized in pigment production, affiliated to the phylum Binatota, may provide a photoprotective function, closing a previously undescribed C1-metabolic loop that involves both the sponge host and specific members of the associated microbial community.

CONCLUSION

Given the global distribution of this ancient animal lineage and their remarkable water filtration activity, sponge-hosted methane cycling may affect methane supersaturation in oxic coastal environments. Depending on the net balance between methane production and consumption, sponges may serve as marine sources or sinks of this potent greenhouse gas. Video Abstract.

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.

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.

The Aging Human Liver: The Weal and Woe of Evolutionary Legacy

Aging is characterized by the progressive decline of biological integrity and its compensatory mechanisms as well as immunological dysregulation. This goes along with an increasing risk of frailty and disease. Against this background, we here specifically focus on the aging of the human liver. For the first time, we shed light on the intertwining evolutionary underpinnings of the liver's declining regenerative capacity, the phenomenon of inflammaging, and the biotransformation capacity in the process of aging. In addition, we discuss how aging influences the risk for developing nonalcoholic fatty liver disease, hepatocellular carcinoma, and/or autoimmune hepatitis, and we describe chronic diseases as accelerators of biological aging.

On the origin of the functional versatility of macrophages

Macrophages represent the most functionally versatile cells in the animal body. In addition to recognizing and destroying pathogens, macrophages remove senescent and exhausted cells, promote wound healing, and govern tissue and metabolic homeostasis. In addition, many specialized populations of tissue-resident macrophages exhibit highly specialized functions essential for the function of specific organs. Sometimes, however, macrophages cease to perform their protective function and their seemingly incomprehensible response to certain stimuli leads to pathology. In this study, we address the question of the origin of the functional versatility of macrophages. To this end, we have searched for the evolutionary origin of macrophages themselves and for the emergence of their characteristic properties. We hypothesize that many of the characteristic features of proinflammatory macrophages evolved in the unicellular ancestors of animals, and that the functional repertoire of macrophage-like amoebocytes further expanded with the evolution of multicellularity and the increasing complexity of tissues and organ systems. We suggest that the entire repertoire of macrophage functions evolved by repurposing and diversification of basic functions that evolved early in the evolution of metazoans under conditions barely comparable to that in tissues of multicellular organisms. We believe that by applying this perspective, we may find an explanation for the otherwise counterintuitive behavior of macrophages in many human pathologies.

Near-Infrared 810 nm Light Affects Porifera Chondrosia reniformis (Nardo, 1847) Regeneration: Molecular Implications and Evolutionary Considerations of Photobiomodulation-Animal Cell Interaction

Chemotrophic choice as a metabolic source of energy has characterised animal cell evolution. However, light interactions with animal cell photoacceptors that are able to increase energetic metabolism (photo-biomodulation (PBM)) have been previously described. In the present study, we cut three specimens of Chondrosia reniformis into four equal parts (12 fragments), and we irradiated the regenerating edge of six fragments with the previously characterised 810 nm near-infrared light, delivered at 1 W, 60 J/cm², 1 W/cm², and 60 J in a continuous-wave mode for 60 s through a flat-top hand-piece with a rounded spot-size area of 1 cm². Six fragments were irradiated with 0 W for 60 s as the controls. We performed irradiation at the time 0 h and every 24 h for a total of five administrations. We monitored the regeneration process for five days (120 h) in aquaria by examining the macroscopic and histological changes. We analysed the gene expression profile of the inflammatory processes, apoptosis, heat stress, growth factors, and collagen production and determined oxidative stress enzyme activity and the total prokaryotic symbiont content. PBM sped up C. reniformis regeneration when compared to the controls. Particularly, transforming growth factor TGF3 and TGF6 upregulation during the early phase of regeneration and TGF5 upregulation 120 h postinjury in the irradiated samples supports the positive effect of PBM in sponge tissue recovery. Conversely, the expression of TGF4, a sponge fibroblast growth factor homologue, was not affected by irradiation, indicating that multiple, independent pathways regulate the TGF genes. The results are consistent with our previous data on a wide range of organisms and humans, suggesting that PBM interaction with primary and secondary cell targets has been conserved through the evolution of life forms.

Resolving tricky nodes in the tree of life through amino acid recoding

Genomic data allowed a detailed resolution of the Tree of Life, but "tricky nodes" such as the root of the animals remain unresolved. Genome-scale datasets are heterogeneous as genes and species are exposed to different pressures, and this can negatively impacts phylogenetic accuracy. We use simulated genomic-scale datasets and show that recoding amino acid data improves accuracy when the model does not account for the compositional heterogeneity of the amino acid alignment. We apply our findings to three datasets addressing the root of the animal tree, where the debate centers on whether sponges (Porifera) or comb jellies (Ctenophora) represent the sister of all other animals. We show that results from empirical data follow predictions from simulations and suggest that, at the least in phylogenies inferred from amino acid sequences, a placement of the ctenophores as sister to all the other animals is best explained as a tree reconstruction artifact.

Phylogenomics and the first higher taxonomy of Placozoa, an ancient and enigmatic animal phylum

Placozoa is an ancient phylum of extraordinarily unusual animals: miniscule, ameboid creatures that lack most fundamental animal features. Despite high genetic diversity, only recently have the second and third species been named. While prior genomic studies suffer from incomplete placozoan taxon sampling, we more than double the count with protein sequences from seven key genomes and produce the first nuclear phylogenomic reconstruction of all major placozoan lineages. This leads us to the first complete Linnaean taxonomic classification of Placozoa, over a century after its discovery: This may be the only time in the 21st century when an entire higher taxonomy for a whole animal phylum is formalized. Our classification establishes 2 new classes, 4 new orders, 3 new families, 1 new genus, and 1 new species, namely classes Polyplacotomia and Uniplacotomia; orders Polyplacotomea, Trichoplacea, Cladhexea, and Hoilungea; families Polyplacotomidae, Cladtertiidae, and Hoilungidae; and genus Cladtertia with species Cladtertia collaboinventa, nov. Our likelihood and gene content tree topologies refine the relationships determined in previous studies. Adding morphological data into our phylogenomic matrices suggests sponges (Porifera) as the sister to other animals, indicating that modest data addition shifts this node away from comb jellies (Ctenophora). Furthermore, by adding the first genomic protein data of the exceptionally distinct and branching Polyplacotoma mediterranea, we solidify its position as sister to all other placozoans; a divergence we estimate to be over 400 million years old. Yet even this deep split sits on a long branch to other animals, suggesting a bottleneck event followed by diversification. Ancestral state reconstructions indicate large shifts in gene content within Placozoa, with Hoilungia hongkongensis and its closest relatives having the most unique genetics.

Filtering artifactual signal increases support for Xenacoelomorpha and Ambulacraria sister relationship in the animal tree of life

Conflicting studies place a group of bilaterian invertebrates containing xenoturbellids and acoelomorphs, the Xenacoelomorpha, as either the primary emerging bilaterian phylum^1,2,3,4,5,6 or within Deuterostomia, sister to Ambulacraria.^7,8,9,10,11 Although their placement as sister to the rest of Bilateria supports relatively simple morphology in the ancestral bilaterian, their alternative placement within Deuterostomia suggests a morphologically complex ancestral bilaterian along with extensive loss of major phenotypic traits in the Xenacoelomorpha. Recent studies have questioned whether Deuterostomia should be considered monophyletic at all.^10,12,13 Hidden paralogy and poor phylogenetic signal present a major challenge for reconstructing species phylogenies.^{14,15,16,17,18} Here, we assess whether these issues have contributed to the conflict over the placement of Xenacoelomorpha. We reanalyzed published datasets, enriching for orthogroups whose gene trees support well-resolved clans elsewhere in the animal tree.¹⁶ We find that most genes in previously published datasets violate incontestable clans, suggesting that hidden paralogy and low phylogenetic signal affect the ability to reconstruct branching patterns at deep nodes in the animal tree. We demonstrate that removing orthogroups that cannot recapitulate incontestable relationships alters the final topology that is inferred, while simultaneously improving the fit of the model to the data. We discover increased, but ultimately not conclusive, support for the existence of Xenambulacraria in our set of filtered orthogroups. At a time when we are progressing toward sequencing all life on the planet, we argue that long-standing contentious issues in the tree of life will be resolved using smaller amounts of better quality data that can be modeled adequately.¹⁹.

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.

A Novel Iridovirus Discovered in Deep-Sea Carnivorous Sponges

Carnivorous sponges (family Cladorhizidae) use small invertebrates as their main source of nutrients. We discovered a novel iridovirus (carnivorous sponge-associated iridovirus, CaSpA-IV) in Chondrocladia grandis and Cladorhiza oxeata specimens collected in the Arctic and Atlantic oceans at depths of 537–852 m. The sequenced viral genome (~190,000 bp) comprised 185 predicted ORFs, including those encoding 26 iridoviral core proteins, and phylogenetic analyses showed that CaSpA-IV is a close relative to members of the genus Decapodiridovirus and highly identical to a partially sequenced virus pathogenic to decapod shrimps. CaSpA-IV was found in various anatomical regions of six C. grandis (sphere, stem, root) from the Gulf of Maine and Baffin Bay and of two C. oxeata (sphere, secondary axis) from Baffin Bay. Partial MCP sequencing revealed a divergent virus (CaSpA-IV-2) in one C. oxeata. The analysis of a 10 nt long tandem repeat showed a number of repeats consistent across sub-sections of the same sponges but different between animals, suggesting the presence of different strains. As the genetic material of crustaceans, particularly from the zooplanktonic copepod order Calanoida, was identified in the investigated samples, further studies are required to elucidate whether CaSpA-IV infects the carnivorous sponges, their crustacean prey, or both.

Polyzoa is back: The effect of complete gene sets on the placement of Ectoprocta and Entoprocta

The phylogenomic approach has largely resolved metazoan phylogeny and improved our knowledge of animal evolution based on morphology, paleontology, and embryology. Nevertheless, the placement of two major lophotrochozoan phyla, Entoprocta (Kamptozoa) and Ectoprocta (Bryozoa), remains highly controversial: Originally considered as a single group named Polyzoa (Bryozoa), they were separated on the basis of morphology. So far, each new study of lophotrochozoan evolution has still consistently proposed different phylogenetic positions for these groups. Here, we reinvestigated the placement of Entoprocta and Ectoprocta using highly complete datasets with rigorous contamination removal. Our results from maximum likelihood, Bayesian, and coalescent analyses strongly support the topology in which Entoprocta and Bryozoa form a distinct clade, placed as a sister group to all other lophotrochozoan clades: Annelida, Mollusca, Brachiopoda, Phoronida, and Nemertea. Our study favors the evolutionary scenario where Entoprocta, Cycliophora, and Bryozoa constitute one of the earliest branches among Lophotrochozoa and thus supports the Polyzoa hypothesis.

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.

Phototransduction in a marine sponge provides insights into the origin of animal vision

Most organisms respond to light. Here, we investigate the origin of metazoan phototransduction by comparing well-characterized opsin-based photosystems in neural animals with those in the sponge Amphimedon queenslandica. Although sponges lack neurons and opsins, they can respond rapidly to light. In Amphimedon larvae, this is guided by the light-sensing posterior pigment ring. We first use cell-type-specific transcriptomes to reveal that genes that characterize eumetazoan Gt- and Go-mediated photosystems are enriched in the pigment ring. We then apply a suite of signaling pathway agonists and antagonists to swimming larvae exposed to directional light. These experiments implicate metabotropic glutamate receptors, phospholipase-C, protein kinase C, and voltage-gated calcium channels in larval phototaxis; the inhibition of phospholipase-C, a key transducer of the Gq-mediated pathway, completely reverses phototactic behavior. Together, these results are consistent with aneural sponges sharing with neural metazoans an ancestral set of photosignaling pathways.

•

Amphimedon larvae are negatively phototactic but lack neurons and opsins

•

Sponge larval photosensory cells are enriched in conserved phototransduction genes

•

Conserved photosignaling pathways appear to be controlling larval phototaxis

•

Phototactic behavior is reversed by the inhibition of phospholipase-C

Damaged Dickinsonia specimens provide clues to Ediacaran vendobiont biology

Recently reported specimens of the enigmatic Ediacaran fossil Dickinsonia from Russia show damage and repair that provides evidence of how they grew, and of their biological affinities. Marginal and terminal areas of wilting deformation are necrotic zones separating regenerated growth, sometimes on two divergent axes, rather than a single axis. Necrotic zones of damage to Dickinsonia are not a thick scar or callus, like a wound or amputation. Nor are they smooth transitions to a regenerated tail or arm. The wilted necrotic zone is most like damage by freezing, salt, or sunburn of leaves and lichens, compatible with evidence of terrestrial habitat from associated frigid and gypsic paleosols. Dickinsonia did not regrow by postembryonic addition of modules from a subterminal or patterned growth zone as in earthworms, myriapods, trilobites, crustaceans, and lizards. Rather Dickinsonia postembryonic regrowth from sublethal damage was from microscopic apical and lateral meristems, as in plants and lichens. Considered as fungal, Dickinsonia, and perhaps others of Class Vendobionta, were more likely Glomeromycota or Mucoromycotina, rather than Ascomycota or Basidiomycota.

Neural Cell Type Diversity in Cnidaria

Neurons are the fundamental building blocks of nervous systems. It appears intuitive that the human brain is made up of hundreds, if not thousands different types of neurons. Conversely, the seemingly diffuse nerve net of Cnidaria is often assumed to be simple. However, evidence that the Cnidaria nervous system is indeed simple is sparse. Recent technical advances make it possible to assess the diversity and function of neurons with unprecedented resolution. Transgenic animals expressing genetically encoded Calcium sensors allow direct physiological assessments of neural responses within the nerve net and provide insight into the spatial organization of the nervous system. Moreover, response and activity patterns allow the characterization of cell types on a functional level. Molecular and genetic identities on the other hand can be assessed combining single-cell transcriptomic analysis with correlations of gene expression in defined neurons. Here I review recent advances on these two experimental strategies focusing on Hydra, Nematostella, and Clytia.

Marine Sponge Endosymbionts: Structural and Functional Specificity of the Microbiome within Euryspongia arenaria Cells

Sponge microbiomes are typically profiled by analyzing the community DNA of whole tissues, which does not distinguish the taxa residing within sponge cells from extracellular microbes. To uncover the endosymbiotic microbiome, we separated the sponge cells to enrich the intracellular microbes. The intracellular bacterial community of sponge Euryspongia arenaria was initially assessed by amplicon sequencing, which indicated that it hosts three unique phyla not found in the extracellular and bulk tissue microbiomes. These three phyla account for 66% of the taxonomically known genera in the intracellular microbiome. The shotgun metagenomic analysis extended the taxonomic coverage to viruses and eukaryotes, revealing the most abundant signature taxa specific to the intracellular microbiome. Functional KEGG pathway annotation demonstrated that the endosymbiotic microbiome hosted the greatest number of unique gene orthologs. The pathway profiles distinguished the intra- and extracellular microbiomes from the tissue and seawater microbiomes. Carbohydrate-active enzyme analysis further discriminated each microbiome based on their representative and dominant enzyme families. One pathway involved in digestion system and family esterase had a consistently higher level in intracellular microbiome and could statistically differentiate the intracellular microbiome from the others, suggesting that triacylglycerol lipases could be the key functional component peculiar to the endosymbionts. The identified higher abundance of lipase-related eggNOG categories further supported the lipid-hydrolyzing metabolism of endosymbiotic microbiota. Pseudomonas members, reported as lipase-producing bacteria, were only in the endosymbiotic microbiome, meanwhile Pseudomonas also showed a greater abundance intracellularly. Our study aided a comprehensive sponge microbiome that demonstrated the taxonomic and functional specificity of endosymbiotic microbiota.

IMPORTANCE Sponges host abundant microbial symbionts that can produce an impressive number of novel bioactive metabolites. However, knowledge on intracellular (endosymbiotic) microbiota is scarce. We characterize the composition and function of the endosymbiotic microbiome by separation of sponge cells and enrichment of intracellular microbes. We uncover a noteworthy number of taxa exclusively in the endosymbiotic microbiome. We unlock the unique pathways and enzymes of endosymbiotic taxa. This study achieves a more comprehensive sponge microbial community profile, which demonstrates the structural and functional specificity of the endosymbiotic microbiome. Our findings not only open the possibility to reveal the low abundant and the likely missed microbiota when directly sequencing the sponge bulk tissues, but also warrant future in-depth exploration within single sponge cells.

Phylogenomic conflict analyses in the apple genus Malus s.l. reveal widespread hybridization and allopolyploidy driving diversification, with insights into the complex biogeographic history in the Northern Hemisphere

Phylogenomic evidence from an increasing number of studies has demonstrated that different data sets and analytical approaches often reconstruct strongly supported but conflicting relationships. In this study, 785 single-copy nuclear genes and 75 complete plastomes were used to infer the phylogenetic relationships and estimate the historical biogeography of the apple genus Malus sensu lato, an economically important lineage disjunctly distributed in the Northern Hemisphere and involved in known and suspected hybridization and allopolyploidy events. The nuclear phylogeny recovered the monophyly of Malus s.l. (including Docynia); however, the genus was supported to be biphyletic in the plastid phylogeny. An ancient chloroplast capture event in the Eocene in western North America best explains the cytonuclear discordance. Our conflict analysis demonstrated that ILS, hybridization, and allopolyploidy could explain the widespread nuclear gene tree discordance. One deep hybridization event (Malus doumeri) and one recent event (Malus coronaria) were detected in Malus s.l. Furthermore, our historical biogeographic analysis integrating living and fossil data supported a widespread East Asian-western North American origin of Malus s.l. in the Eocene, followed by several extinction and dispersal events in the Northern Hemisphere. We also propose a general workflow for assessing phylogenomic discordance and biogeographic analysis using deep genome skimming data sets.