Improvement of molecular phylogenetic inference and the phylogeny of Bilateria

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

Conflicts in Mitochondrial Phylogenomics of Branchiopoda, with the First Complete Mitogenome of Laevicaudata (Crustacea: Branchiopoda)

Conflicting phylogenetic signals are pervasive across genomes. The potential impact of such systematic biases may be reduced by phylogenetic approaches accommodating for heterogeneity or by the exclusive use of homoplastic sites in the datasets. Here, we present the complete mitogenome of Lynceus grossipedia as the first representative of the suborder Laevicaudata. We employed a phylogenomic approach on the mitogenomic datasets representing all major branchiopod groups to identify the presence of conflicts and concordance across the phylogeny. We found pervasive phylogenetic conflicts at the base of Diplostraca. The homogeneity of the substitution pattern tests and posterior predictive tests revealed a high degree of compositional heterogeneity among branchiopod mitogenomes at both the nucleotide and amino acid levels, which biased the phylogenetic inference. Our results suggest that Laevicaudata as the basal clade of Phyllopoda was most likely an artifact caused by compositional heterogeneity and conflicting phylogenetic signal. We demonstrated that the exclusive use of homoplastic site methods combining the application of site-heterogeneous models produced correct phylogenetic estimates of the higher-level relationships among branchiopods.

Molecular phylogenies map to biogeography better than morphological ones

Phylogenetic relationships are inferred principally from two classes of data: morphological and molecular. Currently, most phylogenies of extant taxa are inferred from molecules and when morphological and molecular trees conflict the latter are often preferred. Although supported by simulations, the superiority of molecular trees has rarely been assessed empirically. Here we test phylogenetic accuracy using two independent data sources: biogeographic distributions and fossil first occurrences. For 48 pairs of morphological and molecular trees we show that, on average, molecular trees provide a better fit to biogeographic data than their morphological counterparts and that biogeographic congruence increases over research time. We find no significant differences in stratigraphic congruence between morphological and molecular trees. These results have implications for understanding the distribution of homoplasy in morphological data sets, the utility of morphology as a test of molecular hypotheses and the implications of analysing fossil groups for which molecular data are unavailable.

Using biogeographical and phylogenetic data, it is shown that molecular trees fit species geographical data better than trees inferred from morphology, and that these differences are not simply due to better tree resolution.

Evolution of methyltransferase like (METTL) proteins in Metazoa: A complex gene family involved in epitranscriptomic regulation and other epigenetic processes

The methyltransferase-like (METTL) proteins constitute a family of seven-beta-strand methyltransferases with S-adenosyl methionine-binding domains that modify DNA, RNA, and proteins. Methylation by METTL proteins contributes to the epigenetic, and in the case of RNA modifications, epitranscriptomic regulation of a variety of biological processes. Despite their functional importance, most investigations of the substrates and functions of METTLs within metazoans have been restricted to model vertebrate taxa. In the present work, we explore the evolutionary mechanisms driving the diversification and functional differentiation of 33 individual METTL proteins across Metazoa. Our results show that METTLs are nearly ubiquitous across the animal kingdom, with most having arisen early in metazoan evolution (i.e., occur in basal metazoan phyla). Individual METTL lineages each originated from single independent ancestors, constituting monophyletic clades, which suggests that each METTL was subject to strong selective constraints driving its structural and/or functional specialization. Interestingly, a similar process did not extend to the differentiation of nucleoside-modifying and protein-modifying METTLs (i.e., each METTL type did not form a unique monophyletic clade). The members of these two types of METTLs also exhibited differences in their rates of evolution. Overall, we provide evidence that the long-term evolution of METTL family members was driven by strong purifying selection, which in combination with adaptive selection episodes, led to the functional specialization of individual METTL lineages. This work contributes useful information regarding the evolution of a gene family that fulfills a variety of epigenetic functions, and can have profound influences on molecular processes and phenotypic traits.

Data curation and modeling of compositional heterogeneity in insect phylogenomics: A case study of the phylogeny of Dytiscoidea (Coleoptera: Adephaga)

Diving beetles and their allies are an almost ubiquitous group of freshwater predators. Knowledge of the phylogeny of the adephagan superfamily Dytiscoidea has significantly improved since the advent of molecular phylogenetics. However, despite recent comprehensive phylogenomic studies, some phylogenetic relationships among the constituent families remain elusive. In particular, the position of the family Hygrobiidae remains uncertain. We address these issues by re-analyzing recently published phylogenomic datasets for Dytiscoidea, using approaches to reduce compositional heterogeneity and adopting a site-heterogeneous mixture model. We obtained a consistent, well-resolved, and strongly supported tree. Consistent with previous studies, our analyses support Aspidytidae as the monophyletic sister group of Amphizoidae, and more importantly, Hygrobiidae as the sister of the diverse Dytiscidae, in agreement with morphology-based phylogenies. Our analyses provide a backbone phylogeny of Dytiscoidea, which lays the foundation for better understanding the evolution of morphological characters, life habits, and feeding behaviors of dytiscoid beetles.

New Phylogenomic Analysis of the Enigmatic Phylum Telonemia Further Resolves the Eukaryote Tree of Life

The resolution of the broad-scale tree of eukaryotes is constantly improving, but the evolutionary origin of several major groups remains unknown. Resolving the phylogenetic position of these “orphan” groups is important, especially those that originated early in evolution, because they represent missing evolutionary links between established groups. Telonemia is one such orphan taxon for which little is known. The group is composed of molecularly diverse biflagellated protists, often prevalent although not abundant in aquatic environments. Telonemia has been hypothesized to represent a deeply diverging eukaryotic phylum but no consensus exists as to where it is placed in the tree. Here, we established cultures and report the phylogenomic analyses of three new transcriptome data sets for divergent telonemid lineages. All our phylogenetic reconstructions, based on 248 genes and using site-heterogeneous mixture models, robustly resolve the evolutionary origin of Telonemia as sister to the Sar supergroup. This grouping remains well supported when as few as 60% of the genes are randomly subsampled, thus is not sensitive to the sets of genes used but requires a minimal alignment length to recover enough phylogenetic signal. Telonemia occupies a crucial position in the tree to examine the origin of Sar, one of the most lineage-rich eukaryote supergroups. We propose the moniker “TSAR” to accommodate this new mega-assemblage in the phylogeny of eukaryotes.

Animal Evolution: Last Word on Sponges-First?

A major problem in understanding animal evolution is where early branching phyla, especially sponges and comb jellies (sea gooseberries), sit in the tree of life. A new study seeks to overcome this problem by sampling more species and data cleansing.

Mitigating Anticipated Effects of Systematic Errors Supports Sister-Group Relationship between Xenacoelomorpha and Ambulacraria

Xenoturbella and the acoelomorph worms (Xenacoelomorpha) are simple marine animals with controversial affinities. They have been placed as the sister group of all other bilaterian animals (Nephrozoa hypothesis), implying their simplicity is an ancient characteristic [1, 2]; alternatively, they have been linked to the complex Ambulacraria (echinoderms and hemichordates) in a clade called the Xenambulacraria [3-5], suggesting their simplicity evolved by reduction from a complex ancestor. The difficulty resolving this problem implies the phylogenetic signal supporting the correct solution is weak and affected by inadequate modeling, creating a misleading non-phylogenetic signal. The idea that the Nephrozoa hypothesis might be an artifact is prompted by the faster molecular evolutionary rate observed within the Acoelomorpha. Unequal rates of evolution are known to result in the systematic artifact of long branch attraction, which would be predicted to result in an attraction between long-branch acoelomorphs and the outgroup, pulling them toward the root [6]. Other biases inadequately accommodated by the models used can also have strong effects, exacerbated in the context of short internal branches and long terminal branches [7]. We have assembled a large and informative dataset to address this problem. Analyses designed to reduce or to emphasize misleading signals show the Nephrozoa hypothesis is supported under conditions expected to exacerbate errors, and the Xenambulacraria hypothesis is preferred in conditions designed to reduce errors. Our reanalyses of two other recently published datasets [1, 2] produce the same result. We conclude that the Xenacoelomorpha are simplified relatives of the Ambulacraria.

Asgard archaea capable of anaerobic hydrocarbon cycling

Large reservoirs of natural gas in the oceanic subsurface sustain complex communities of anaerobic microbes, including archaeal lineages with potential to mediate oxidation of hydrocarbons such as methane and butane. Here we describe a previously unknown archaeal phylum, Helarchaeota, belonging to the Asgard superphylum and with the potential for hydrocarbon oxidation. We reconstruct Helarchaeota genomes from metagenomic data derived from hydrothermal deep-sea sediments in the hydrocarbon-rich Guaymas Basin. The genomes encode methyl-CoM reductase-like enzymes that are similar to those found in butane-oxidizing archaea, as well as several enzymes potentially involved in alkyl-CoA oxidation and the Wood-Ljungdahl pathway. We suggest that members of the Helarchaeota have the potential to activate and subsequently anaerobically oxidize hydrothermally generated short-chain hydrocarbons.

Evaluating Model Performance in Evolutionary Biology

Many fields of evolutionary biology now depend on stochastic mathematical models. These models are valuable for their ability to formalize predictions in the face of uncertainty and provide a quantitative framework for testing hypotheses. However, no mathematical model will fully capture biological complexity. Instead, these models attempt to capture the important features of biological systems using relatively simple mathematical principles. These simplifications can allow us to focus on differences that are meaningful, while ignoring those that are not. However, simplification also requires assumptions, and to the extent that these are wrong, so is our ability to predict or compare. Here, we discuss approaches for evaluating the performance of evolutionary models in light of their assumptions by comparing them against reality. We highlight general approaches, how they are applied, and remaining opportunities. Absolute tests of fit, even when not explicitly framed as such, are fundamental to progress in understanding evolution.

A Traveling Wave Ion Mobility Spectrometry (TWIMS) Study of the Robo1-Heparan Sulfate Interaction

Roundabout 1 (Robo1) interacts with its receptor Slit to regulate axon guidance, axon branching, and dendritic development in the nervous system and to regulate morphogenesis and many cell functions in the nonneuronal tissues. This interaction is known to be critically regulated by heparan sulfate (HS). Previous studies suggest that HS is required to promote the binding of Robo1 to Slit to form the minimal signaling complex, but the molecular details and the structural requirements of HS for this interaction are still unclear. Here, we describe the application of traveling wave ion mobility spectrometry (TWIMS) to study the conformational details of the Robo1-HS interaction. The results suggest that Robo1 exists in two conformations that differ by their compactness and capability to interact with HS. The results also suggest that the highly flexible interdomain hinge region connecting the Ig1 and Ig2 domains of Robo1 plays an important functional role in promoting the Robo1-Slit interaction. Moreover, variations in the sulfation pattern and size of HS were found to affect its binding affinity and selectivity to interact with different conformations of Robo1. Both MS measurements and CIU experiments show that the Robo1-HS interaction requires the presence of a specific size and pattern of modification of HS. Furthermore, the effect of N-glycosylation on the conformation of Robo1 and its binding modes with HS is reported. Graphical Abstract ᅟ.

Combined morphological and phylogenomic re-examination of malawimonads, a critical taxon for inferring the evolutionary history of eukaryotes

Modern syntheses of eukaryote diversity assign almost all taxa to one of three groups: Amorphea, Diaphoretickes and Excavata (comprising Discoba and Metamonada). The most glaring exception is Malawimonadidae, a group of small heterotrophic flagellates that resemble Excavata by morphology, but branch with Amorphea in most phylogenomic analyses. However, just one malawimonad, Malawimonas jakobiformis, has been studied with both morphological and molecular-phylogenetic approaches, raising the spectre of interpretation errors and phylogenetic artefacts from low taxon sampling. We report a morphological and phylogenomic study of a new deep-branching malawimonad, Gefionella okellyi n. gen. n. sp. Electron microscopy revealed all canonical features of 'typical excavates', including flagellar vanes (as an opposed pair, unlike M. jakobiformis but like many metamonads) and a composite fibre. Initial phylogenomic analyses grouped malawimonads with the Amorphea-related orphan lineage Collodictyon, separate from a Metamonada+Discoba clade. However, support for this topology weakened when more sophisticated evolutionary models were used, and/or fast-evolving sites and long-branching taxa (FS/LB) were excluded. Analyses of '-FS/LB' datasets instead suggested a relationship between malawimonads and metamonads. The 'malawimonad+metamonad signal' in morphological and molecular data argues against a strict Metamonada+Discoba clade (i.e. the predominant concept of Excavata). A Metamonad+Discoba clade should therefore not be assumed when inferring deep-level evolutionary history in eukaryotes.

Between a Pod and a Hard Test: The Deep Evolution of Amoebae

Amoebozoa is the eukaryotic supergroup sister to Obazoa, the lineage that contains the animals and Fungi, as well as their protistan relatives, and the breviate and apusomonad flagellates. Amoebozoa is extraordinarily diverse, encompassing important model organisms and significant pathogens. Although amoebozoans are integral to global nutrient cycles and present in nearly all environments, they remain vastly understudied. We present a robust phylogeny of Amoebozoa based on broad representative set of taxa in a phylogenomic framework (325 genes). By sampling 61 taxa using culture-based and single-cell transcriptomics, our analyses show two major clades of Amoebozoa, Discosea, and Tevosa. This phylogeny refutes previous studies in major respects. Our results support the hypothesis that the last common ancestor of Amoebozoa was sexual and flagellated, it also may have had the ability to disperse propagules from a sporocarp-type fruiting body. Overall, the main macroevolutionary patterns in Amoebozoa appear to result from the parallel losses of homologous characters of a multiphase life cycle that included flagella, sex, and sporocarps rather than independent acquisition of convergent features.

Segmentation in Tardigrada and diversification of segmental patterns in Panarthropoda

The origin and diversification of segmented metazoan body plans has fascinated biologists for over a century. The superphylum Panarthropoda includes three phyla of segmented animals-Euarthropoda, Onychophora, and Tardigrada. This superphylum includes representatives with relatively simple and representatives with relatively complex segmented body plans. At one extreme of this continuum, euarthropods exhibit an incredible diversity of serially homologous segments. Furthermore, distinct tagmosis patterns are exhibited by different classes of euarthropods. At the other extreme, all tardigrades share a simple segmented body plan that consists of a head and four leg-bearing segments. The modular body plans of panarthropods make them a tractable model for understanding diversification of animal body plans more generally. Here we review results of recent morphological and developmental studies of tardigrade segmentation. These results complement investigations of segmentation processes in other panarthropods and paleontological studies to illuminate the earliest steps in the evolution of panarthropod body plans.

Multiple measures could alleviate long-branch attraction in phylogenomic reconstruction of Cupressoideae (Cupressaceae)

Long-branch attraction (LBA) is a major obstacle in phylogenetic reconstruction. The phylogenetic relationships among Juniperus (J), Cupressus (C) and the Hesperocyparis-Callitropsis-Xanthocyparis (HCX) subclades of Cupressoideae are controversial. Our initial analyses of plastid protein-coding gene matrix revealed both J and C with much longer stem branches than those of HCX, so their sister relationships may be attributed to LBA. We used multiple measures including data filtering and modifying, evolutionary model selection and coalescent phylogenetic reconstruction to alleviate the LBA artifact. Data filtering by strictly removing unreliable aligned regions and removing substitution saturation genes and rapidly evolving sites could significantly reduce branch lengths of subclades J and C and recovered a relationship of J (C, HCX). In addition, using coalescent phylogenetic reconstruction could elucidate the LBA artifact and recovered J (C, HCX). However, some valid methods for other taxa were inefficient in alleviating the LBA artifact in J-C-HCX. Different strategies should be carefully considered and justified to reduce LBA in phylogenetic reconstruction of different groups. Three subclades of J-C-HCX were estimated to have experienced ancient rapid divergence within a short period, which could be another major obstacle in resolving relationships. Furthermore, our plastid phylogenomic analyses fully resolved the intergeneric relationships of Cupressoideae.

Exploring Phylogenetic Relationships within Myriapoda and the Effects of Matrix Composition and Occupancy on Phylogenomic Reconstruction

Myriapods, including the diverse and familiar centipedes and millipedes, are one of the dominant terrestrial arthropod groups. Although molecular evidence has shown that Myriapoda is monophyletic, its internal phylogeny remains contentious and understudied, especially when compared to those of Chelicerata and Hexapoda. Until now, efforts have focused on taxon sampling (e.g., by including a handful of genes from many species) or on maximizing matrix size (e.g., by including hundreds or thousands of genes in just a few species), but a phylogeny maximizing sampling at both levels remains elusive. In this study, we analyzed 40 Illumina transcriptomes representing 3 of the 4 myriapod classes (Diplopoda, Chilopoda, and Symphyla); 25 transcriptomes were newly sequenced to maximize representation at the ordinal level in Diplopoda and at the family level in Chilopoda. Ten supermatrices were constructed to explore the effect of several potential phylogenetic biases (e.g., rate of evolution, heterotachy) at 3 levels of gene occupancy per taxon (50%, 75%, and 90%). Analyses based on maximum likelihood and Bayesian mixture models retrieved monophyly of each myriapod class, and resulted in 2 alternative phylogenetic positions for Symphyla, as sister group to Diplopoda + Chilopoda, or closer to Diplopoda, the latter hypothesis having been traditionally supported by morphology. Within centipedes, all orders were well supported, but 2 deep nodes remained in conflict in the different analyses despite dense taxon sampling at the family level. Relationships among centipede orders in all analyses conducted with the most complete matrix (90% occupancy) are at odds not only with the sparser but more gene-rich supermatrices (75% and 50% supermatrices) and with the matrices optimizing phylogenetic informativeness or most conserved genes, but also with previous hypotheses based on morphology, development, or other molecular data sets. Our results indicate that a high percentage of ribosomal proteins in the most complete matrices, in conjunction with distance from the root, can act in concert to compromise the estimated relationships within the ingroup. We discuss the implications of these findings in the context of the ever more prevalent quest for completeness in phylogenomic studies.

Slit-Robo signaling

Slits are secreted proteins that bind to Roundabout (Robo) receptors. Slit-Robo signaling is best known for mediating axon repulsion in the developing nervous system. However, in recent years the functional repertoire of Slits and Robo has expanded tremendously and Slit-Robo signaling has been linked to roles in neurogenesis, angiogenesis and cancer progression among other processes. Likewise, our mechanistic understanding of Slit-Robo signaling has progressed enormously. Here, we summarize new insights into Slit-Robo evolutionary and system-dependent diversity, receptor-ligand interactions, signaling crosstalk and receptor activation.

Heterokont Predator Develorapax marinus gen. et sp. nov. - A Model of the Ochrophyte Ancestor

Heterotrophic lineages of Heterokonta (or stramenopiles), in contrast to a single monophyletic group of autotrophs, Ochrophyta, form several clades that independently branch off the heterokont stem lineage. The nearest neighbors of Ochrophyta in the phylogenetic tree appear to be almost exclusively bacterivorous, whereas the hypothesis of plastid acquisition by the ancestors of the ochrophyte lineage suggests an ability to engulf eukaryotic alga. In line with this hypothesis, the heterotrophic predator at the base of the ochrophyte lineage may be regarded as a model for the ochrophyte ancestor. Here, we present a new genus and species of marine free-living heterotrophic heterokont Develorapax marinus, which falls into an isolated heterokont cluster, along with the marine flagellate Developayella elegans, and is able to engulf eukaryotic cells. Together with environmental sequences D. marinus and D. elegans form a class-level clade Developea nom. nov. represented by species adapted to different environmental conditions and with a wide geographical distribution. The position of Developea among Heterokonta in large-scale phylogenetic tree is discussed. We propose that members of the Developea clade represent a model for transition from bacterivory to a predatory feeding mode by selection for larger prey. Presumably, such transition in the grazing strategy is possible in the presence of bacterial biofilms or aggregates expected in eutrophic environment, and has likely occurred in the ochrophyte ancestor.

Phylogenomic Insights into Animal Evolution

Animals make up only a small fraction of the eukaryotic tree of life, yet, from our vantage point as members of the animal kingdom, the evolution of the bewildering diversity of animal forms is endlessly fascinating. In the century following the publication of Darwin's Origin of Species, hypotheses regarding the evolution of the major branches of the animal kingdom - their relationships to each other and the evolution of their body plans - was based on a consideration of the morphological and developmental characteristics of the different animal groups. This morphology-based approach had many successes but important aspects of the evolutionary tree remained disputed. In the past three decades, molecular data, most obviously primary sequences of DNA and proteins, have provided an estimate of animal phylogeny largely independent of the morphological evolution we would ultimately like to understand. The molecular tree that has evolved over the past three decades has drastically altered our view of animal phylogeny and many aspects of the tree are no longer contentious. The focus of molecular studies on relationships between animal groups means, however, that the discipline has become somewhat divorced from the underlying biology and from the morphological characteristics whose evolution we aim to understand. Here, we consider what we currently know of animal phylogeny; what aspects we are still uncertain about and what our improved understanding of animal phylogeny can tell us about the evolution of the great diversity of animal life.

Rooting the tree of life: the phylogenetic jury is still out

This article aims to shed light on difficulties in rooting the tree of life (ToL) and to explore the (sociological) reasons underlying the limited interest in accurately addressing this fundamental issue. First, we briefly review the difficulties plaguing phylogenetic inference and the ways to improve the modelling of the substitution process, which is highly heterogeneous, both across sites and over time. We further observe that enriched taxon samplings, better gene samplings and clever data removal strategies have led to numerous revisions of the ToL, and that these improved shallow phylogenies nearly always relocate simple organisms higher in the ToL provided that long-branch attraction artefacts are kept at bay. Then, we note that, despite the flood of genomic data available since 2000, there has been a surprisingly low interest in inferring the root of the ToL. Furthermore, the rare studies dealing with this question were almost always based on methods dating from the 1990s that have been shown to be inaccurate for much more shallow issues! This leads us to argue that the current consensus about a bacterial root for the ToL can be traced back to the prejudice of Aristotle's Great Chain of Beings, in which simple organisms are ancestors of more complex life forms. Finally, we demonstrate that even the best models cannot yet handle the complexity of the evolutionary process encountered both at shallow depth, when the outgroup is too distant, and at the level of the inter-domain relationships. Altogether, we conclude that the commonly accepted bacterial root is still unproven and that the root of the ToL should be revisited using phylogenomic supermatrices to ensure that new evidence for eukaryogenesis, such as the recently described Lokiarcheota, is interpreted in a sound phylogenetic framework.

Chloroplast Phylogenomic Inference of Green Algae Relationships

The green algal phylum Chlorophyta has six diverse classes, but the phylogenetic relationship of the classes within Chlorophyta remains uncertain. In order to better understand the ancient Chlorophyta evolution, we have applied a site pattern sorting method to study compositional heterogeneity and the model fit in the green algal chloroplast genomic data. We show that the fastest-evolving sites are significantly correlated with among-site compositional heterogeneity, and these sites have a much poorer fit to the evolutionary model. Our phylogenomic analyses suggest that the class Chlorophyceae is a monophyletic group, and the classes Ulvophyceae, Trebouxiophyceae and Prasinophyceae are non-monophyletic groups. Our proposed phylogenetic tree of Chlorophyta will offer new insights to investigate ancient green algae evolution, and our analytical framework will provide a useful approach for evaluating and mitigating the potential errors of phylogenomic inferences.

Untangling the early diversification of eukaryotes: a phylogenomic study of the evolutionary origins of Centrohelida, Haptophyta and Cryptista

Assembling the global eukaryotic tree of life has long been a major effort of Biology. In recent years, pushed by the new availability of genome-scale data for microbial eukaryotes, it has become possible to revisit many evolutionary enigmas. However, some of the most ancient nodes, which are essential for inferring a stable tree, have remained highly controversial. Among other reasons, the lack of adequate genomic datasets for key taxa has prevented the robust reconstruction of early diversification events. In this context, the centrohelid heliozoans are particularly relevant for reconstructing the tree of eukaryotes because they represent one of the last substantial groups that was missing large and diverse genomic data. Here, we filled this gap by sequencing high-quality transcriptomes for four centrohelid lineages, each corresponding to a different family. Combining these new data with a broad eukaryotic sampling, we produced a gene-rich taxon-rich phylogenomic dataset that enabled us to refine the structure of the tree. Specifically, we show that (i) centrohelids relate to haptophytes, confirming Haptista; (ii) Haptista relates to SAR; (iii) Cryptista share strong affinity with Archaeplastida; and (iv) Haptista + SAR is sister to Cryptista + Archaeplastida. The implications of this topology are discussed in the broader context of plastid evolution.

Probabilistic models of eukaryotic evolution: time for integration

In spite of substantial work and recent progress, a global and fully resolved picture of the macroevolutionary history of eukaryotes is still under construction. This concerns not only the phylogenetic relations among major groups, but also the general characteristics of the underlying macroevolutionary processes, including the patterns of gene family evolution associated with endosymbioses, as well as their impact on the sequence evolutionary process. All these questions raise formidable methodological challenges, calling for a more powerful statistical paradigm. In this direction, model-based probabilistic approaches have played an increasingly important role. In particular, improved models of sequence evolution accounting for heterogeneities across sites and across lineages have led to significant, although insufficient, improvement in phylogenetic accuracy. More recently, one main trend has been to move away from simple parametric models and stepwise approaches, towards integrative models explicitly considering the intricate interplay between multiple levels of macroevolutionary processes. Such integrative models are in their infancy, and their application to the phylogeny of eukaryotes still requires substantial improvement of the underlying models, as well as additional computational developments.

Evolutionary relatedness does not predict competition and co-occurrence in natural or experimental communities of green algae

The competition-relatedness hypothesis (CRH) predicts that the strength of competition is the strongest among closely related species and decreases as species become less related. This hypothesis is based on the assumption that common ancestry causes close relatives to share biological traits that lead to greater ecological similarity. Although intuitively appealing, the extent to which phylogeny can predict competition and co-occurrence among species has only recently been rigorously tested, with mixed results. When studies have failed to support the CRH, critics have pointed out at least three limitations: (i) the use of data poor phylogenies that provide inaccurate estimates of species relatedness, (ii) the use of inappropriate statistical models that fail to detect relationships between relatedness and species interactions amidst nonlinearities and heteroskedastic variances, and (iii) overly simplified laboratory conditions that fail to allow eco-evolutionary relationships to emerge. Here, we address these limitations and find they do not explain why evolutionary relatedness fails to predict the strength of species interactions or probabilities of coexistence among freshwater green algae. First, we construct a new data-rich, transcriptome-based phylogeny of common freshwater green algae that are commonly cultured and used for laboratory experiments. Using this new phylogeny, we re-analyse ecological data from three previously published laboratory experiments. After accounting for the possibility of nonlinearities and heterogeneity of variances across levels of relatedness, we find no relationship between phylogenetic distance and ecological traits. In addition, we show that communities of North American green algae are randomly composed with respect to their evolutionary relationships in 99% of 1077 lakes spanning the continental United States. Together, these analyses result in one of the most comprehensive case studies of how evolutionary history influences species interactions and community assembly in both natural and experimental systems. Our results challenge the generality of the CRH and suggest it may be time to re-evaluate the validity and assumptions of this hypothesis.

Are 100 enough? Inferring acanthomorph teleost phylogeny using Anchored Hybrid Enrichment

BACKGROUND

The past decade has witnessed remarkable progress towards resolution of the Tree of Life. However, despite the increased use of genomic scale datasets, some phylogenetic relationships remain difficult to resolve. Here we employ anchored phylogenomics to capture 107 nuclear loci in 29 species of acanthomorph teleost fishes, with 25 of these species sampled from the recently delimited clade Ovalentaria. Previous studies employing multilocus nuclear exon datasets have not been able to resolve the nodes at the base of the Ovalentaria tree with confidence. Here we test whether a phylogenomic approach will provide better support for these nodes, and if not, why this may be.

RESULTS

After using a novel method to account for paralogous loci, we estimated phylogenies with maximum likelihood and species tree methods using DNA sequence alignments of over 80,000 base pairs. Several key relationships within Ovalentaria are well resolved, including 1) the sister taxon relationship between Cichlidae and Pholidichthys, 2) a clade containing blennies, grammas, clingfishes, and jawfishes, and 3) monophyly of Atherinomorpha (topminnows, flyingfishes, and silversides). However, many nodes in the phylogeny associated with the early diversification of Ovalentaria are poorly resolved in several analyses. Through the use of rarefaction curves we show that limited phylogenetic resolution among the earliest nodes in the Ovalentaria phylogeny does not appear to be due to a deficiency of data, as average global node support ceases to increase when only 1/3rd of the sampled loci are used in analyses. Instead this lack of resolution may be driven by model misspecification as a Bayesian mixed model analysis of the amino acid dataset provided good support for parts of the base of the Ovalentaria tree.

CONCLUSIONS

Although it does not appear that the limited phylogenetic resolution among the earliest nodes in the Ovalentaria phylogeny is due to a deficiency of data, it may be that both stochastic and systematic error resulting from model misspecification play a role in the poor resolution at the base of the Ovalentaria tree as a Bayesian approach was able to resolve some of the deeper nodes, where the other methods failed.

A Paleogene origin for crown passerines and the diversification of the Oscines in the New World

In this study, we present a detailed family-level phylogenetic hypothesis for the largest avian order (Aves: Passeriformes) and an unmatched multi-calibrated, relaxed clock inference for the diversification of crown passerines. Extended taxon sampling allowed the recovery of many challenging clades and elucidated their position in the tree. Acanthisittia appear to have diverged from all other passerines at the early Paleogene, which is considerably later than previously suggested. Thus, Passeriformes may be younger and represent an even more intense adaptive radiation compared to the remaining avian orders. Based on our divergence time estimates, a novel hypothesis for the diversification of modern Suboscines is proposed. According to this hypothesis, the first split between New and Old World lineages would be related to the severing of the Africa-South America biotic connection during the mid-late Eocene, implying an African origin for modern Eurylaimides. The monophyletic status of groups not recovered by any subsequent study since their circumscription, viz. Sylvioidea including Paridae, Remizidae, Hyliotidae, and Stenostiridae; and Muscicapoidea including the waxwing assemblage (Bombycilloidea) were notable topological findings. We also propose possible ecological interactions that may have shaped the distinct Oscine distribution patterns in the New World. The insectivorous endemic Oscines of the Americas, Vireonidae (Corvoidea), Mimidae, and Troglodytidae (Muscicapoidea), probably interfered with autochthonous Suboscines through direct competition. Thus, the Early Miocene arrival of these lineages before any other Oscines may have occupied the few available niches left by Tyrannides, constraining the diversification of insectivorous Oscines that arrived in the Americas later. The predominantly frugivorous-nectarivorous members of Passeroidea, which account for most of the diversity of New World-endemic Oscines, may not have been subjected to competition with Tyrannides. In fact, the vast availability of frugivory niches combined with weak competition with the autochthonous passerine fauna may have been crucial for passeroids to thrive in the New World.

A role for the Parkinson's disease protein DJ-1 as a chaperone and antioxidant in the anhydrobiotic nematode Panagrolaimus superbus

Mutations in the human DJ-1/PARK7 gene are associated with familial Parkinson's disease. DJ-1 belongs to a large, functionally diverse family with homologues in all biological kingdoms. Several activities have been demonstrated for DJ-1: an antioxidant protein, a redox-regulated molecular chaperone and a modulator of multiple cellular signalling pathways. The majority of functional studies have focussed on human DJ-1 (hDJ-1), but studies on DJ-1 homologues in Drosophila melanogaster, Caenorhabditis elegans, Dugesia japonica and Escherichia coli also provide evidence of a role for DJ-1 as an antioxidant. Here, we show that dehydration is a potent inducer of a dj-1 gene in the anhydrobiotic nematode Panagrolaimus superbus. Our secondary structure and homology modelling analyses shows that recombinant DJ-1 protein from P. superbus (PsuDJ-1.1) is a well-folded protein, which is similar in structure to the hDJ-1. PsuDJ-1.1 is a heat stable protein; with T1/2 unfolding transition values of 76 and 70 °C obtained from both circular dichroism (CD) and Fourier transform infrared spectroscopy (FTIR) measurements respectively. We found that PsuDJ-1.1 is an efficient antioxidant that also functions as a 'holdase' molecular chaperone that can maintain its chaperone function in a reducing environment. In addition to its chaperone activity, PsuDJ-1.1 may also be an important non-enzymatic antioxidant, capable of providing protection to P. superbus from oxidative damage when the nematodes are in a desiccated, anhydrobiotic state.

A transcriptome approach to ecdysozoan phylogeny

The monophyly of Ecdysozoa, which comprise molting phyla, has received strong support from several lines of evidence. However, the internal relationships of Ecdysozoa are still contended. We generated expressed sequence tags from a priapulid (penis worm), a kinorhynch (mud dragon), a tardigrade (water bear) and five chelicerate taxa by 454 transcriptome sequencing. A multigene alignment was assembled from 63 taxa, which comprised after matrix optimization 24,249 amino acid positions with high data density (2.6% gaps, 19.1% missing data). Phylogenetic analyses employing various models support the monophyly of Ecdysozoa. A clade combining Priapulida and Kinorhyncha (i.e. Scalidophora) was recovered as the earliest branch among Ecdysozoa. We conclude that Cycloneuralia, a taxon erected to combine Priapulida, Kinorhyncha and Nematoda (and others), are paraphyletic. Rather Arthropoda (including Onychophora) are allied with Nematoda and Tardigrada. Within Arthropoda, we found strong support for most clades, including monophyletic Mandibulata and Pancrustacea. The phylogeny within the Euchelicerata remained largely unresolved. There is conflicting evidence on the position of tardigrades: While Bayesian and maximum likelihood analyses of only slowly evolving genes recovered Tardigrada as a sister group to Arthropoda, analyses of the full data set, and of subsets containing genes evolving at fast and intermediate rates identified a clade of Tardigrada and Nematoda. Notably, the latter topology is also supported by the analyses of indel patterns.

Evaluating topological conflict in centipede phylogeny using transcriptomic data sets

Relationships between the five extant orders of centipedes have been considered solved based on morphology. Phylogenies based on samples of up to a few dozen genes have largely been congruent with the morphological tree apart from an alternative placement of one order, the relictual Craterostigmomorpha, consisting of two species in Tasmania and New Zealand. To address this incongruence, novel transcriptomic data were generated to sample all five orders of centipedes and also used as a test case for studying gene-tree incongruence. Maximum likelihood and Bayesian mixture model analyses of a data set composed of 1,934 orthologs with 45% missing data, as well as the 389 orthologs in the least saturated, stationary quartile, retrieve strong support for a sister-group relationship between Craterostigmomorpha and all other pleurostigmophoran centipedes, of which the latter group is newly named Amalpighiata. The Amalpighiata hypothesis, which shows little gene-tree incongruence and is robust to the influence of among-taxon compositional heterogeneity, implies convergent evolution in several morphological and behavioral characters traditionally used in centipede phylogenetics, such as maternal brood care, but accords with patterns of first appearances in the fossil record.

Selective neuronal staining in tardigrades and onychophorans provides insights into the evolution of segmental ganglia in panarthropods

BACKGROUND

Although molecular analyses have contributed to a better resolution of the animal tree of life, the phylogenetic position of tardigrades (water bears) is still controversial, as they have been united alternatively with nematodes, arthropods, onychophorans (velvet worms), or onychophorans plus arthropods. Depending on the hypothesis favoured, segmental ganglia in tardigrades and arthropods might either have evolved independently, or they might well be homologous, suggesting that they were either lost in onychophorans or are a synapomorphy of tardigrades and arthropods. To evaluate these alternatives, we analysed the organisation of the nervous system in three tardigrade species using antisera directed against tyrosinated and acetylated tubulin, the amine transmitter serotonin, and the invertebrate neuropeptides FMRFamide, allatostatin and perisulfakinin. In addition, we performed retrograde staining of nerves in the onychophoran Euperipatoides rowelli in order to compare the serial locations of motor neurons within the nervous system relative to the appendages they serve in arthropods, tardigrades and onychophorans.

RESULTS

Contrary to a previous report from a Macrobiotus species, our immunocytochemical and electron microscopic data revealed contralateral fibres and bundles of neurites in each trunk ganglion of three tardigrade species, including Macrobiotus cf. harmsworthi, Paramacrobiotus richtersi and Hypsibius dujardini. Moreover, we identified additional, extra-ganglionic commissures in the interpedal regions bridging the paired longitudinal connectives. Within the ganglia we found serially repeated sets of serotonin- and RFamid-like immunoreactive neurons. Furthermore, our data show that the trunk ganglia of tardigrades, which include the somata of motor neurons, are shifted anteriorly with respect to each corresponding leg pair, whereas no such shift is evident in the arrangement of motor neurons in the onychophoran nerve cords.

CONCLUSIONS

Taken together, these data reveal three major correspondences between the segmental ganglia of tardigrades and arthropods, including (i) contralateral projections and commissures in each ganglion, (ii) segmentally repeated sets of immunoreactive neurons, and (iii) an anteriorly shifted (parasegmental) position of ganglia. These correspondences support the homology of segmental ganglia in tardigrades and arthropods, suggesting that these structures were either lost in Onychophora or, alternatively, evolved in the tardigrade/arthropod lineage.

von Baer's law for the ages: lost and found principles of developmental evolution

In 1828, Karl Ernst von Baer formulated a series of empirically defined rules, which became widely known as the 'Law of Development' or 'von Baer's law of embryology'. This was one the most significant attempts to define the principles that connected morphological complexity and embryonic development. Understanding this relation is central to both evolutionary biology and developmental genetics. Von Baer's ideas have been both a source of inspiration to generations of biologists and a target of continuous criticism over many years. With advances in multiple fields, including paleontology, cladistics, phylogenetics, genomics, and cell and developmental biology, it is now possible to examine carefully the significance of von Baer's law and its predictions. In this review, I argue that, 185 years after von Baer's law was first formulated, its main concepts after proper refurbishing remain surprisingly relevant in revealing the fundamentals of the evolution-development connection, and suggest that their explanation should become the focus of renewed research.

Phylogenomics demonstrates that breviate flagellates are related to opisthokonts and apusomonads

Most eukaryotic lineages belong to one of a few major groups. However, several protistan lineages have not yet been robustly placed in any of these groups. Both the breviates and apusomonads are two such lineages that appear to be related to the Amoebozoa and Opisthokonta (i.e. the 'unikonts' or Amorphea); however, their precise phylogenetic positions remain unclear. Here, we describe a novel microaerophilic breviate, Pygsuia biforma gen. nov. sp. nov., isolated from a hypoxic estuarine sediment. Ultrastructurally, this species resembles the breviate genera Breviata and Subulatomonas but has two cell morphologies, adherent and swimming. Phylogenetic analyses of the small sub-unit rRNA gene show that Pygsuia is the sister to the other breviates. We constructed a 159-protein supermatrix, including orthologues identified in RNA-seq data from Pygsuia. Phylogenomic analyses of this dataset show that breviates, apusomonads and Opisthokonta form a strongly supported major eukaryotic grouping we name the Obazoa. Although some phylogenetic methods disagree, the balance of evidence suggests that the breviate lineage forms the deepest branch within Obazoa. We also found transcripts encoding a nearly complete integrin adhesome from Pygsuia, indicating that this protein complex involved in metazoan multicellularity may have evolved earlier in eukaryote evolution than previously thought.

Jumping and gliding rodents: mitogenomic affinities of Pedetidae and Anomaluridae deduced from an RNA-Seq approach

An RNA-Seq strategy was used to obtain the complete set of protein-coding mitochondrial genes from two rodent taxa. Thanks to the next generation sequencing (NGS) 454 approach, we determined the complete mitochondrial DNA genome from Graphiurus kelleni (Mammalia: Rodentia: Gliridae) and partial mitogenome from Pedetes capensis (Pedetidae), and compared them with published rodent and outgroup mitogenomes. We finished the mitogenome sequencing by a series of amplicons using conserved PCR primers to fill the gaps corresponding to tRNA, rRNA and control regions. Phylogenetic analyses of the mitogenomes suggest a well-supported rodent phylogeny in agreement with nuclear gene trees. Pedetes groups with Anomalurus into the clade Anomaluromorpha, while Graphiurus branches within the squirrel-related clade. Moreover, Pedetes+Anomalurus branch with Castor into the mouse-related clade. Our study demonstrates the utility of NGS for obtaining new mitochondrial genomes as well as the importance of choosing adequate models of sequence evolution to infer the phylogeny of rodents.

PhyloBayes MPI: phylogenetic reconstruction with infinite mixtures of profiles in a parallel environment

Modeling across site variation of the substitution process is increasingly recognized as important for obtaining more accurate phylogenetic reconstructions. Both finite and infinite mixture models have been proposed and have been shown to significantly improve on classical single-matrix models. Compared with their finite counterparts, infinite mixtures have a greater expressivity. However, they are computationally more challenging. This has resulted in practical compromises in the design of infinite mixture models. In particular, a fast but simplified version of a Dirichlet process model over equilibrium frequency profiles implemented in PhyloBayes has often been used in recent phylogenomics studies, while more refined model structures, more realistic and empirically more fit, have been practically out of reach. We introduce a message passing interface version of PhyloBayes, implementing the Dirichlet process mixture models as well as more classical empirical matrices and finite mixtures. The parallelization is made efficient thanks to the combination of two algorithmic strategies: a partial Gibbs sampling update of the tree topology and the use of a truncated stick-breaking representation for the Dirichlet process prior. The implementation shows close to linear gains in computational speed for up to 64 cores, thus allowing faster phylogenetic reconstruction under complex mixture models. PhyloBayes MPI is freely available from our website www.phylobayes.org.

TPMS: a set of utilities for querying collections of gene trees

Background

The information in large collections of phylogenetic trees is useful for many comparative genomic studies. Therefore, there is a need for flexible tools that allow exploration of such collections in order to retrieve relevant data as quickly as possible.

Results

In this paper, we present TPMS (Tree Pattern-Matching Suite), a set of programs for handling and retrieving gene trees according to different criteria. The programs from the suite include utilities for tree collection building, specific tree-pattern search strategies and tree rooting. Use of TPMS is illustrated through three examples: systematic search for incongruencies in a large tree collection, a short study on the Coelomata/Ecdysozoa controversy and an evaluation of the level of support for a recently published Mammal phylogeny.

Conclusion

TPMS is a powerful suite allowing to quickly retrieve sets of trees matching complex patterns in large collection or to root trees using more rigorous approaches than the classical midpoint method. As it is made of a set of command-line programs, it can be easily integrated in any sequence analysis pipeline for an automated use.

Neural markers reveal a one-segmented head in tardigrades (water bears)

BACKGROUND

While recent neuroanatomical and gene expression studies have clarified the alignment of cephalic segments in arthropods and onychophorans, the identity of head segments in tardigrades remains controversial. In particular, it is unclear whether the tardigrade head and its enclosed brain comprises one, or several segments, or a non-segmental structure. To clarify this, we applied a variety of histochemical and immunocytochemical markers to specimens of the tardigrade Macrobiotus cf. harmsworthi and the onychophoran Euperipatoides rowelli.

METHODOLOGY/PRINCIPAL FINDINGS

Our immunolabelling against serotonin, FMRFamide and α-tubulin reveals that the tardigrade brain is a dorsal, bilaterally symmetric structure that resembles the brain of onychophorans and arthropods rather than a circumoesophageal ring typical of cycloneuralians (nematodes and allies). A suboesophageal ganglion is clearly lacking. Our data further reveal a hitherto unknown, unpaired stomatogastric ganglion in Macrobiotus cf. harmsworthi, which innervates the ectodermal oesophagus and the endodermal midgut and is associated with the second leg-bearing segment. In contrast, the oesophagus of the onychophoran E. rowelli possesses no immunoreactive neurons, whereas scattered bipolar, serotonin-like immunoreactive cell bodies are found in the midgut wall. Furthermore, our results show that the onychophoran pharynx is innervated by a medullary loop nerve accompanied by monopolar, serotonin-like immunoreactive cell bodies.

CONCLUSIONS/SIGNIFICANCE

A comparison of the nervous system innervating the foregut and midgut structures in tardigrades and onychophorans to that of arthropods indicates that the stomatogastric ganglion is a potential synapomorphy of Tardigrada and Arthropoda. Its association with the second leg-bearing segment in tardigrades suggests that the second trunk ganglion is a homologue of the arthropod tritocerebrum, whereas the first ganglion corresponds to the deutocerebrum. We therefore conclude that the tardigrade brain consists of a single segmental region corresponding to the arthropod protocerebrum and, accordingly, that the tardigrade head is a non-composite, one-segmented structure.