Should we be worried about long-branch attraction in real data sets? Investigations using metazoan 18S rDNA

Maximum Likelihood Estimation for Unrooted 3-Leaf Trees: An Analytic Solution for the CFN Model

Maximum likelihood estimation is among the most widely-used methods for inferring phylogenetic trees from sequence data. This paper solves the problem of computing solutions to the maximum likelihood problem for 3-leaf trees under the 2-state symmetric mutation model (CFN model). Our main result is a closed-form solution to the maximum likelihood problem for unrooted 3-leaf trees, given generic data; this result characterizes all of the ways that a maximum likelihood estimate can fail to exist for generic data and provides theoretical validation for predictions made in Parks and Goldman (Syst Biol 63(5):798-811, 2014). Our proof makes use of both classical tools for studying group-based phylogenetic models such as Hadamard conjugation and reparameterization in terms of Fourier coordinates, as well as more recent results concerning the semi-algebraic constraints of the CFN model. To be able to put these into practice, we also give a complete characterization to test genericity.

A New Record of Pogonatum tahitense (Polytrichaceae) from Tibet, China: Taxonomic Description, Range Expansion, and Biogeographic History

The genus Pogonatum stands out as the most diverse within the family Polytrichaceae, encompassing over 50 species. Pogonatum tahitense has been recorded across various Pacific regions, including Hawaii in the United States and Tahiti in French Polynesia, as well as in Asia, such as in Taiwan in China, Java in Indonesia, and Sabah in Malaysia. In the current study, a specimen collected in Tibet, China, is described, confirming its taxonomic classification as P. tahitense through a comprehensive analysis integrating morphological evidence and molecular study based on sequences from the plastid (rbcL, rps4, trnL-F), mitochondrial (nad5), and nuclear (ITS2) regions. This documentation represents the first record of the species within mainland China. A time-calibrated, molecular-based phylogenetic analysis was conducted, employing various approaches for ancestral range inference. The findings suggest that P. tahitense originated during the Pleistocene epoch, approximately 1.8 mya, in Tibet, China.

Novel symmetry-preserving neural network model for phylogenetic inference

Motivation

Scientists world-wide are putting together massive efforts to understand how the biodiversity that we see on Earth evolved from single-cell organisms at the origin of life and this diversification process is represented through the Tree of Life. Low sampling rates and high heterogeneity in the rate of evolution across sites and lineages produce a phenomenon denoted "long branch attraction" (LBA) in which long nonsister lineages are estimated to be sisters regardless of their true evolutionary relationship. LBA has been a pervasive problem in phylogenetic inference affecting different types of methodologies from distance-based to likelihood-based.

Results

Here, we present a novel neural network model that outperforms standard phylogenetic methods and other neural network implementations under LBA settings. Furthermore, unlike existing neural network models in phylogenetics, our model naturally accounts for the tree isomorphisms via permutation invariant functions which ultimately result in lower memory and allows the seamless extension to larger trees.

Availability and implementation

We implement our novel theory on an open-source publicly available GitHub repository: https://github.com/crsl4/nn-phylogenetics.

Extensive nrDNA ITS polymorphism in Lycium: Non-concerted evolution and the identification of pseudogenes

The internal transcribed spacer (ITS) is one of the most extensively sequenced molecular markers in plant systematics due to its generally concerted evolution. While non-concerted evolution has been found in some plant taxa, such information is missing in Lycium. Molecular studies of six species and two variants of the genus Lycium revealed high levels of intra- and inter-individual polymorphism in the ITS, indicating non-concerted evolution. All genomic DNA ITS paralogues were identified as putative pseudogenes or functional paralogues through a series of comparisons of sequence features, including length and substitution variation, GC content, secondary structure stability, and the presence of conserved motifs in the 5.8S gene, and the rate of evolution. Approximately, 60% of ITS pseudogenes could be easily detected. Based on phylogenetic analysis, all pseudogenes were highly distinct from their corresponding functional copies, tended to evolve neutrally, and clustered randomly together in the evolutionary tree. The results probably suggest that this ITS non-concerted evolution is related to the recent divergence between tandem repeats within the Lycium genome and hybridization between species. Our study complements those of pseudogenes in plant taxa and provides a theoretical basis for the phylogeny and genetic origin of the genus Lycium while having important implications for the use of ITS molecular markers for phylogenetic reconstruction.

Evolutionary Rate Variation Among Lineages in Gene Trees has a Negative Impact on Species-Tree Inference

Phylogenetic analyses of genomic data provide a powerful means of reconstructing the evolutionary relationships among organisms, yet such analyses are often hindered by conflicting phylogenetic signals among loci. Identifying the signals that are most influential to species-tree estimation can help to inform the choice of data for phylogenomic analysis. We investigated this in an analysis of 30 phylogenomic data sets. For each data set, we examined the association between several branch-length characteristics of gene trees and the distance between these gene trees and the corresponding species trees. We found that the distance of each gene tree to the species tree inferred from the full data set was positively associated with variation in root-to-tip distances and negatively associated with mean branch support. However, no such associations were found for gene-tree length, a measure of the overall substitution rate at each locus. We further explored the usefulness of the best-performing branch-based characteristics for selecting loci for phylogenomic analyses. We found that loci that yield gene trees with high variation in root-to-tip distances have a disproportionately distant signal of tree topology compared with the complete data sets. These results suggest that rate variation across lineages should be taken into consideration when exploring and even selecting loci for phylogenomic analysis.[Branch support; data filtering; nucleotide substitution model; phylogenomics; substitution rate; summary coalescent methods.]

A combined molecular and morphological phylogeny of the Loricariinae (Siluriformes: Loricariidae), with emphasis on the Harttiini and Farlowellini

We present a combined molecular and morphological phylogenetic analysis of the Loricariinae, with emphasis on the Harttiini (Cteniloricaria, Harttia, and Harttiella) and Farlowellini (Aposturisoma, Farlowella, Lamontichthys, Pterosturisoma, Sturisoma, and Sturisomatichthys). Character sampling comprised seven molecular markers (the mitochondrial Cytb, nd2, 12S and 16S, and the nuclear MyH6, RAG1 and RAG2) and 196 morphological characters. A total of 1,059 specimens, and 159 tissue samples were analized, representing 100 species. A Bayesian Inference analysis was performed using the concatenated data matrix, which is comprised of 6,819 characters. The Loricariinae were found to comprise the tribes (Hartiini (Loricariini, Farlowellini)), the latter two elevated from subtribes. A Maximum Parsimony analysis was also performed using the same data matrix in order to reveal phenotypical synapomorphies to diagnose each clade. Two MP trees were found with a length of 14,704 steps, consistency index of 0.29 and retention index of 0.61, which were summarized in a strict consensus tree. Harttiini includes (Harttiella (Cteniloricaria, Harttia), and Farlowellini includes (Lamontichthys (Pterosturisoma (Sturisoma (Sturisomatichthys, Farlowella)))). Aposturisoma was recovered nested within Farlowella and is synonymyzed to the latter. Sturisoma was corroborated as strictly cis-Andean, while Sturisomatichthys encompasses, besides the valid species already included in the genus, the trans-Andean species once belonging to Sturisoma sensu lato. Identification keys and phylogenetic diagnoses of family-group taxa and genera of both the Harttiini and the Farlowellini are provided.

HOST-SWITCHING EVENTS IN LITOMOSOIDES CHANDLER, 1931 (FILARIOIDEA: ONCHOCERCIDAE) ARE NOT RAMPANT BUT CLADE DEPENDENT

The genus Litomosoides Chandler, 1931, includes species that as adults occur in the thoracic and abdominal cavity of mammalian hosts and are presumably vectored by mites. The vertebrate hosts include a variety of Neotropical mammals such as phyllostomid and mormoopid bats; cricetid, sciurid, and hystricognath rodents; and didelphid marsupials. It has been suggested that Litomosoides is not a monophyletic group and that rampant horizontal transfer explains their presence in disparate groups of mammals. Herein we present a phylogenetic reconstruction including mitochondrial genes of 13 vouchered species. This phylogeny is used to reconstruct the evolutionary history of these parasites and the ancestral states of key characters used in species classification, namely, the configuration of the spicules. The historical association of these filarioids with 6 groups of mammals, as well as their ancestral geographic distributions, were reconstructed using Bayesian statistical approaches comparing alternative models of biogeography and evolution and fossil states in selected nodes of the phylogeny. The optimal reconstruction suggests a model of dispersal, extinction, and cladogenesis (DEC) driving the evolution of Litomosoides; the results suggest an origin of Litomosoides in South America and association of ancestors with phyllostomids, and strong evidence of at least 2 host-switching events: 1 of these involving cricetid rodents and the other mormoopid bats. The latter event included a simultaneous geographic expansion of the parasite lineage across South and North America. The host-switching event from phyllostomid bats into cricetid rodents occurred once these rodents diversified across South America; subsequent diversification of the latter clade resulted in 2 branches, each showing expansion of the parasites back into North America. This result suggests that both parasites and cricetid rodents established an association in South America, underwent diversification, and then dispersed into North America. Further, this clade of cricetid-dwelling species includes parasites featuring the "sigmodontis" spicule type. The identification of a single host-switching event involving the disparate lineages of Chiroptera and Rodentia offers a framework to reconstruct the gene evolution and diversification of this lineage after the host-switching event. This will help in predicting the ability of these parasites to infect sympatric mammals.

Higher-level phylogenetic relationships of rove beetles (Coleoptera, Staphylinidae) inferred from mitochondrial genome sequences

Rove beetles (Staphylinidae) and allied families constitute a huge radiation of Coleoptera, but basal relationships in this group remain controversial. In this study, we newly sequenced eight mitogenomes of representatives of Staphylinidae by using next-generation sequencing method. Together with 99 existing mitogenomes of Staphyliniformia, (sub)family relationships were investigated with ML and Bayesian searches under various substitution models and data recoding schemes. The results consistently supported Scydmaenidae and Silphidae to be subordinate groups of Staphylinidae. Within the monophyletic Staphylinidae (including Scydmaenidae and Silphidae), the hypothesis of four major subfamily groups cannot be confirmed. Bayesian inferences under the site-heterogeneous mixture model generally supported the basal position of major clades corresponding to the Omaliine group. At the subfamily level, the monophyly of Pselaphinae, Oxytelinae, Scaphidiinae, Steninae and Staphylininae was supported. However, the subfamilies Omaliinae, Tachyporinae, Aleocharinae and Paederinae were each non-monophyletic.

A Novel Approach to Investigate the Effect of Tree Reconstruction Artifacts in Single-Gene Analysis Clarifies Opsin Evolution in Nonbilaterian Metazoans

Our ability to correctly reconstruct a phylogenetic tree is strongly affected by both systematic errors and the amount of phylogenetic signal in the data. Current approaches to tackle tree reconstruction artifacts, such as the use of parameter-rich models, do not translate readily to single-gene alignments. This, coupled with the limited amount of phylogenetic information contained in single-gene alignments, makes gene trees particularly difficult to reconstruct. Opsin phylogeny illustrates this problem clearly. Opsins are G-protein coupled receptors utilized in photoreceptive processes across Metazoa and their protein sequences are roughly 300 amino acids long. A number of incongruent opsin phylogenies have been published and opsin evolution remains poorly understood. Here, we present a novel approach, the canary sequence approach, to investigate and potentially circumvent errors in single-gene phylogenies. First, we demonstrate our approach using two well-understood cases of long-branch attraction in single-gene data sets, and simulations. After that, we apply our approach to a large collection of well-characterized opsins to clarify the relationships of the three main opsin subfamilies.

Testing adequacy for DNA substitution models

Background

Testing model adequacy is important before a DNA substitution model is chosen for phylogenetic inference. Using a mis-specified model can negatively impact phylogenetic inference, for example, the maximum likelihood method can be inconsistent when the DNA sequences are generated under a tree topology which is in the Felsentein Zone and analyzed with a mis-specified or inadequate model. However, model adequacy testing in phylogenetics is underdeveloped.

Results

Here we develop a simple, general, powerful and robust model test based on Pearson’s goodness-of-fit test and binning of site patterns. We demonstrate through simulation that this test is robust in its high power to reject the inadequate models for a large range of different ways of binning site patterns while the Type I error is controlled well. In the real data analysis we discovered many cases where models chosen by another method can be rejected by this new test, in particular, our proposed test rejects the most complex DNA model (GTR+I+ Γ) while the Goldman-Cox test fails to reject the commonly used simple models.

Conclusions

Model adequacy testing and bootstrap should be used together to assess reliability of conclusions after model selection and model fitting have already been applied to choose the model and fit it. The new goodness-of-fit test proposed in this paper is a simple and powerful model adequacy testing method serving such a regular model checking purpose. We caution against deriving strong conclusions from analyses based on inadequate models. At a minimum, those results derived from inadequate models can now be readly flagged using the new test, and reported as such.

A Critical Appraisal of the Placement of Xiphosura (Chelicerata) with Account of Known Sources of Phylogenetic Error

Horseshoe crabs (Xiphosura) are traditionally regarded as sister group to the clade of terrestrial chelicerates (Arachnida). This hypothesis has been challenged by recent phylogenomic analyses, but the non-monophyly of Arachnida has consistently been disregarded as artifactual. We re-evaluated the placement of Xiphosura among chelicerates using the most complete phylogenetic data set to date, expanding outgroup sampling, and including data from whole genome sequencing projects. In spite of uncertainty in the placement of some arachnid clades, all analyses show Xiphosura consistently nested within Arachnida as the sister group to Ricinulei (hooded tick spiders). It is apparent that the radiation of arachnids is an old one and occurred over a brief period of time, resulting in several consecutive short internodes, and thus is a potential case for the confounding effects of incomplete lineage sorting (ILS). We simulated coalescent gene trees to explore the effects of increasing levels of ILS on the placement of horseshoe crabs. In addition, common sources of systematic error were evaluated, as well as the effects of fast-evolving partitions and the dynamics of problematic long branch orders. Our results indicated that the placement of horseshoe crabs cannot be explained by missing data, compositional biases, saturation, or ILS. Interrogation of the phylogenetic signal showed that the majority of loci favor the derived placement of Xiphosura over a monophyletic Arachnida. Our analyses support the inference that horseshoe crabs represent a group of aquatic arachnids, comparable to aquatic mites, breaking a long-standing paradigm in chelicerate evolution and altering previous interpretations of the ancestral transition to the terrestrial habitat. Future studies testing chelicerate relationships should approach the task with a sampling strategy where the monophyly of Arachnida is not held as the premise.

Statistical Inconsistency of Maximum Parsimony for k-Tuple-Site Data

One of the main aims of phylogenetics is to reconstruct the "Tree of Life." In this respect, different methods and criteria are used to analyze DNA sequences of different species and to compare them in order to derive the evolutionary relationships of these species. Maximum parsimony is one such criterion for tree reconstruction, and it is the one which we will use in this paper. However, it is well known that tree reconstruction methods can lead to wrong relationship estimates. One typical problem of maximum parsimony is long branch attraction, which can lead to statistical inconsistency. In this work, we will consider a blockwise approach to alignment analysis, namely the so-called k-tuple analyses. For four taxa, it has already been shown that k-tuple-based analyses are statistically inconsistent if and only if the standard character-based (site-based) analyses are statistically inconsistent. So, in the four-taxon case, going from individual sites to k-tuples does not lead to any improvement. However, real biological analyses often consider more than only four taxa. Therefore, we analyze the case of five taxa for 2- and 3-tuple-site data and consider alphabets with two and four elements. We show that the equivalence of single-site data and k-tuple-site data then no longer holds. Even so, we can show that maximum parsimony is statistically inconsistent for k-tuple-site data and five taxa.

New patellogastropod mitogenomes help counteracting long-branch attraction in the deep phylogeny of gastropod mollusks

Long-branch attraction (LBA) is a well-known artifact in phylogenetic reconstruction. Sparse taxon sampling and extreme heterogeneity of evolutionary rates among lineages generate propitious situations for LBA, even defying probabilistic methods of phylogenetic inference. A clear example illustrating LBA challenges is the difficulty of reconstructing the deep gastropod phylogeny, particularly using mitochondrial (mt) genomes. Previous studies consistently obtained unorthodox phylogenetic relationships due to the LBA between the mitogenomes of patellogastropods (true limpets, represented only by Lottia digitalis), heterobranchs, and outgroup taxa. Here, we use the reconstruction of the gastropod mitogenomic phylogeny as a case exercise to test the effect of key methodological approaches proposed to counteract LBA, including the selection of slow-evolving representatives, the use of different outgroups, the application of site-heterogeneous evolutionary models, and the removal of fast-evolving sites. In this regard, we sequenced three new patellogastropod mt genomes, which displayed shorter branches than the one of Lottia as well as gene organizations more similar to that of the hypothetical gastropod ancestor. Phylogenetic analyses incorporating the mt genomes of Patella ferruginea, Patella vulgata, and Cellana radiata allowed eliminating the artificial clustering of Patellogastropoda and Heterobranchia that had prevailed in previous studies. Furthermore, the use of site-heterogeneous models with certain combinations of lineages within the outgroup allowed eliminating also the LBA between Heterobranchia and the outgroup, and recovering Apogastropoda (i.e., Caenogastropoda + Heterobranchia). Hence, for the first time, we were able to obtain a mitogenomic phylogeny of gastropods that is congruent with both morphological and nuclear datasets.

Embracing Uncertainty in Reconstructing Early Animal Evolution

The origin of animals, one of the major transitions in evolution, remains mysterious. Many key aspects of animal evolution can be reconstructed by comparing living species within a robust phylogenetic framework. However, uncertainty remains regarding the evolutionary relationships between two ancient animal lineages - sponges and ctenophores - and the remaining animal phyla. Comparative morphology and some phylogenomic analyses support the view that sponges represent the sister lineage to the rest of the animals, while other phylogenomic analyses support ctenophores, a phylum of carnivorous, gelatinous marine organisms, as the sister lineage. Here, we explore why different studies yield different answers and discuss the implications of the two alternative hypotheses for understanding the origin of animals. Reconstruction of ancient evolutionary radiations is devilishly difficult and will likely require broader sampling of sponge and ctenophore genomes, improved analytical strategies and critical analyses of the phylogenetic distribution and molecular mechanisms underlying apparently conserved traits. Rather than staking out positions in favor of the ctenophores-sister or the sponges-sister hypothesis, we submit that research programs aimed at understanding the biology of the first animals should instead embrace the uncertainty surrounding early animal evolution in their experimental designs.

Evolution of Eukaryal and Archaeal Pseudouridine Synthase Pus10

In archaea, pseudouridine (Ψ) synthase Pus10 modifies uridine (U) to Ψ at positions 54 and 55 of tRNA. In contrast, Pus10 is not found in bacteria, where modifications at those two positions are carried out by TrmA (U54 to m⁵U54) and TruB (U55 to Ψ55). Many eukaryotes have an apparent redundancy; their genomes contain orthologs of archaeal Pus10 and bacterial TrmA and TruB. Although eukaryal Pus10 genes share a conserved catalytic domain with archaeal Pus10 genes, their biological roles are not clear for the two reasons. First, experimental evidence suggests that human Pus10 participates in apoptosis induced by the tumor necrosis factor-related apoptosis-inducing ligand. Whether the function of human Pus10 is in place or in addition to of Ψ synthesis in tRNA is unknown. Second, Pus10 is found in earlier evolutionary branches of fungi (such as chytrid Batrachochytrium) but is absent in all dikaryon fungi surveyed (Ascomycetes and Basidiomycetes). We did a comprehensive analysis of sequenced genomes and found that orthologs of Pus10, TrmA, and TruB were present in all the animals, plants, and protozoa surveyed. This indicates that the common eukaryotic ancestor possesses all the three genes. Next, we examined 116 archaeal and eukaryotic Pus10 protein sequences to find that Pus10 existed as a single copy gene in all the surveyed genomes despite ancestral whole genome duplications had occurred. This indicates a possible deleterious gene dosage effect. Our results suggest that functional redundancy result in gene loss or neofunctionalization in different evolutionary lineages.

Can quartet analyses combining maximum likelihood estimation and Hennigian logic overcome long branch attraction in phylogenomic sequence data?

Systematic biases such as long branch attraction can mislead commonly relied upon model-based (i.e. maximum likelihood and Bayesian) phylogenetic methods when, as is usually the case with empirical data, there is model misspecification. We present PhyQuart, a new method for evaluating the three possible binary trees for any quartet of taxa. PhyQuart was developed through a process of reciprocal illumination between a priori considerations and the results of extensive simulations. It is based on identification of site-patterns that can be considered to support a particular quartet tree taking into account the Hennigian distinction between apomorphic and plesiomorphic similarity, and employing corrections to the raw observed frequencies of site-patterns that exploit expectations from maximum likelihood estimation. We demonstrate through extensive simulation experiments that, whereas maximum likeilihood estimation performs well in many cases, it can be outperformed by PhyQuart in cases where it fails due to extreme branch length asymmetries producing long-branch attraction artefacts where there is only very minor model misspecification.

Phylogenetic analysis of the genus Laparocerus, with comments on colonisation and diversification in Macaronesia (Coleoptera, Curculionidae, Entiminae)

The flightless Entiminae weevil genus Laparocerus is the species-richest genus, with 237 species and subspecies, inhabiting Macaronesia (Madeira archipelago, Selvagens, Canary Islands) and the continental 'Macaronesian enclave' in Morocco (one single polytypic species). This is the second contribution to gain insight of the genus and assist in its systematic revision with a mitochondrial phylogenetic analysis. It centres on the Canarian clade, adding the 12S rRNA gene to the combined set of COII and 16S rRNA used in our first contribution on the Madeiran clade (here re-analysed). The nuclear 28S rRNA was also used to produce an additional 4-gene tree to check coherency with the 3-gene tree. A total of 225 taxa (95%) has been sequenced, mostly one individual per taxa. Plausible explanations for incoherent data (mitochondrial introgressions, admixture, incomplete lineage sorting, etc.) are discussed for each of the monophyletic subclades that are coincident with established subgenera, or are restructured or newly described. The overall mean genetic divergence (p-distance) among species is 8.2%; the mean divergence within groups (subgenera) ranks from 2.9 to 7.0% (average 4.6%), and between groups, from 5.4% to 12.0% (average 9.2%). A trustful radiation event within a young island (1.72 Ma) was used to calibrate and produce a chronogram using the software RelTime. These results confirm the monophyly of both the Madeiran (36 species and subspecies) and the Canarian (196 species and subspecies) clades, which originated ca. 11.2 Ma ago, and started to radiate in their respective archipelagos ca. 8.5 and 7.7 Ma ago. The Madeiran clade seems to have begun in Porto Santo, and from there it jumped to the Desertas and to Madeira, with additional radiations. The Canarian clade shows a sequential star-shape radiation process generating subclades with a clear shift from East to West in coherence with the decreasing age of the islands. Laparocerus garretai from the Selvagens belongs to a Canarian subclade, and Laparocerus susicus from Morocco does not represent the ancestral continental lineage, but a back-colonisation from the Canaries to Africa. Dispersal processes, colonisation patterns, and ecological remarks are amply discussed. Diversification has been adaptive as well as non-adaptive, and the role of 'geological turbulence' is highlighted as one of the principal drivers of intra-island allopatric speciation. Based on the phylogenetic results, morphological features and distribution, five new monophyletic subgenera are described: Aridotroxsubg. n., Belicariussubg. n., Bencomiussubg. n., Canariotroxsubg. n., and Purpuraniussubg. n., totalling twenty subgenera in Laparocerus.

Motility, morphology and phylogeny of the plasmodial worm, Ceratomyxa vermiformis n. sp. (Cnidaria: Myxozoa: Myxosporea)

The Myxozoa demonstrate extensive morphological simplification and miniaturization relative to their free-living cnidarian ancestors. This is particularly pronounced in the highly derived myxosporeans, which develop as plasmodia and pseudoplasmodia. To date, motility in these stages has been linked with membrane deformation (e.g. as pseudopodia and mobile folds). Here we illustrate a motile, elongate plasmodium that undergoes coordinated undulatory locomotion, revealing remarkable convergence to a functional worm at the cellular level. Ultrastructural and confocal analyses of these plasmodia identify a highly differentiated external layer containing an actin-rich network, long tubular mitochondria, abundant microtubules, a secreted glycocalyx layer, and an internal region where sporogony occurs and which contains homogeneously distributed granular/fibrillar material. We consider how some of these features may support motility. We also describe the species based on spore morphology and SSU rDNA sequence data, undertake molecular phylogenetic analysis to place it within an early-diverging clade of the ceratomyxids, and evaluate the resultant implications for classification (validity of the genus Meglitschia) and for inferring early host environments (freshwater) of ceratomyxids.

Phylogenomic analysis of carangimorph fishes reveals flatfish asymmetry arose in a blink of the evolutionary eye

Background

Flatfish cranial asymmetry represents one of the most remarkable morphological innovations among vertebrates, and has fueled vigorous debate on the manner and rate at which strikingly divergent phenotypes evolve. A surprising result of many recent molecular phylogenetic studies is the lack of support for flatfish monophyly, where increasingly larger DNA datasets of up to 23 loci have either yielded a weakly supported flatfish clade or indicated the group is polyphyletic. Lack of resolution for flatfish relationships has been attributed to analytical limitations for dealing with processes such as nucleotide non-stationarity and incomplete lineage sorting (ILS). We tackle this phylogenetic problem using a sequence dataset comprising more than 1,000 ultraconserved DNA element (UCE) loci covering 45 carangimorphs, the broader clade containing flatfishes and several other specialized lineages such as remoras, billfishes, and archerfishes.

Results

We present a phylogeny based on UCE loci that unequivocally supports flatfish monophyly and a single origin of asymmetry. We document similar levels of discordance among UCE loci as in previous, smaller molecular datasets. However, relationships among flatfishes and carangimorphs recovered from multilocus concatenated and species tree analyses of our data are robust to the analytical framework applied and size of data matrix used. By integrating the UCE data with a rich fossil record, we find that the most distinctive carangimorph bodyplans arose rapidly during the Paleogene (66.0–23.03 Ma). Flatfish asymmetry, for example, likely evolved over an interval of no more than 2.97 million years.

Conclusions

The longstanding uncertainty in phylogenetic hypotheses for flatfishes and their carangimorph relatives highlights the limitations of smaller molecular datasets when applied to successive, rapid divergences. Here, we recovered significant support for flatfish monophyly and relationships among carangimorphs through analysis of over 1,000 UCE loci. The resulting time-calibrated phylogeny points to phenotypic divergence early within carangimorph history that broadly matches with the predictions of adaptive models of lineage diversification.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0786-x) contains supplementary material, which is available to authorized users.

The methanogenic redox cofactor F420 is widely synthesized by aerobic soil bacteria

F₄₂₀ is a low-potential redox cofactor that mediates the transformations of a wide range of complex organic compounds. Considered one of the rarest cofactors in biology, F₄₂₀ is best known for its role in methanogenesis and has only been chemically identified in two phyla to date, the Euryarchaeota and Actinobacteria. In this work, we show that this cofactor is more widely distributed than previously reported. We detected the genes encoding all five known F₄₂₀ biosynthesis enzymes (cofC, cofD, cofE, cofG and cofH) in at least 653 bacterial and 173 archaeal species, including members of the dominant soil phyla Proteobacteria, Chloroflexi and Firmicutes. Metagenome datamining validated that these genes were disproportionately abundant in aerated soils compared with other ecosystems. We confirmed through high-performance liquid chromatography analysis that aerobically grown stationary-phase cultures of three bacterial species, Paracoccus denitrificans, Oligotropha carboxidovorans and Thermomicrobium roseum, synthesized F₄₂₀, with oligoglutamate sidechains of different lengths. To understand the evolution of F₄₂₀ biosynthesis, we also analyzed the distribution, phylogeny and genetic organization of the cof genes. Our data suggest that although the F_o precursor to F₄₂₀ originated in methanogens, F₄₂₀ itself was first synthesized in an ancestral actinobacterium. F₄₂₀ biosynthesis genes were then disseminated horizontally to archaea and other bacteria. Together, our findings suggest that the cofactor is more significant in aerobic bacterial metabolism and soil ecosystem composition than previously thought. The cofactor may confer several competitive advantages for aerobic soil bacteria by mediating their central metabolic processes and broadening the range of organic compounds they can synthesize, detoxify and mineralize.

Species complexes and phylogenetic lineages of Hoferellus (Myxozoa, Cnidaria) including revision of the genus: A problematic case for taxonomy

Background

Myxozoans are metazoan parasites whose traditional spore morphology-based taxonomy conflicts DNA based phylogenies. Freshwater species of the genus Hoferellus are parasites of the excretory system, with several members infecting food and ornamental fish species, as well as amphibians. This study aims to increase our understanding of their molecular diversity and development, aspects about which little is known, and to generate a molecular diagnostic tool to discriminate between different pathogenic and non-pathogenic Hoferellus spp.

Methods

SSU and ITS rDNA phylogeny, along with morphological descriptions using light and electron microscopy were used to identify and characterize Hoferellus species collected from the urinary system of fishes and frogs. A PCR-based diagnostic assay was designed to differentiate between cryptic Hoferellus spp in cyprinid fishes commonly cultured in Central Europe.

Results

Our phylogenetic results separate the species of Hoferellus into two phylogenetic sublineages which are indistinguishable on the basis of generic morphological traits: 1) The Hoferellus sensu stricto sublineage, which is composed of the type species Hoferellus cyprini, Hoferellus carassii and a cryptic species, Hoferellus sp. detected only molecularly in common carp. 2) The Hoferellus sensu lato sublineage into which the new species we described in this study, Hoferellus gnathonemi sp. n. from the kidney of the elephantnose fish and Hoferellus anurae from reed frogs, are placed together with Hoferellus gilsoni previously sequenced from European eel. Apart from phylogenetic analyses, we also provide novel ultrastructural data on the phagocytotic nature of some Hoferellus plasmodia and on the elusive intracellular stages ascribed to the presporogonic development of this genus.

Conclusions

We provide molecular evidence of the polyphyly of the genus Hoferellus and provide novel morphological details of its members. Based on the presented data, we revise and propose emendation of the genus Hoferellus.

Phylogenetic relationships of Hemiptera inferred from mitochondrial and nuclear genes

Here, we reconstructed the Hemiptera phylogeny based on the expanded mitochondrial protein-coding genes and the nuclear 18S rRNA gene, separately. The differential rates of change across lineages may associate with long-branch attraction (LBA) effect and result in conflicting estimates of phylogeny from different types of data. To reduce the potential effects of systematic biases on inferences of topology, various data coding schemes, site removal method, and different algorithms were utilized in phylogenetic reconstruction. We show that the outgroups Phthiraptera, Thysanoptera, and the ingroup Sternorrhyncha share similar base composition, and exhibit "long branches" relative to other hemipterans. Thus, the long-branch attraction between these groups is suspected to cause the failure of recovering Hemiptera under the homogeneous model. In contrast, a monophyletic Hemiptera is supported when heterogeneous model is utilized in the analysis. Although higher level phylogenetic relationships within Hemiptera remain to be answered, consensus between analyses is beginning to converge on a stable phylogeny.

Utility of characters evolving at diverse rates of evolution to resolve quartet trees with unequal branch lengths: analytical predictions of long-branch effects

BACKGROUND

The detection and avoidance of "long-branch effects" in phylogenetic inference represents a longstanding challenge for molecular phylogenetic investigations. A consequence of parallelism and convergence, long-branch effects arise in phylogenetic inference when there is unequal molecular divergence among lineages, and they can positively mislead inference based on parsimony especially, but also inference based on maximum likelihood and Bayesian approaches. Long-branch effects have been exhaustively examined by simulation studies that have compared the performance of different inference methods in specific model trees and branch length spaces.

RESULTS

In this paper, by generalizing the phylogenetic signal and noise analysis to quartets with uneven subtending branches, we quantify the utility of molecular characters for resolution of quartet phylogenies via parsimony. Our quantification incorporates contributions toward the correct tree from either signal or homoplasy (i.e. "the right result for either the right reason or the wrong reason"). We also characterize a highly conservative lower bound of utility that incorporates contributions to the correct tree only when they correspond to true, unobscured parsimony-informative sites (i.e. "the right result for the right reason"). We apply the generalized signal and noise analysis to classic quartet phylogenies in which long-branch effects can arise due to unequal rates of evolution or an asymmetrical topology. Application of the analysis leads to identification of branch length conditions in which inference will be inconsistent and reveals insights regarding how to improve sampling of molecular loci and taxa in order to correctly resolve phylogenies in which long-branch effects are hypothesized to exist.

CONCLUSIONS

The generalized signal and noise analysis provides analytical prediction of utility of characters evolving at diverse rates of evolution to resolve quartet phylogenies with unequal branch lengths. The analysis can be applied to identifying characters evolving at appropriate rates to resolve phylogenies in which long-branch effects are hypothesized to occur.

Evolutionary origin of Ceratonova shasta and phylogeny of the marine myxosporean lineage

In order to clarify the phylogenetic relationships among the main marine myxosporean clades including newly established Ceratonova clade and scrutinizing their evolutionary origins, we performed large-scale phylogenetic analysis of all myxosporean species from the marine myxosporean lineage based on three gene analyses and statistical topology tests. Furthermore, we obtained new molecular data for Ceratonova shasta, C. gasterostea, eight Ceratomyxa species and one Myxodavisia species. We described five new species: Ceratomyxa ayami n. sp., C. leatherjacketi n. sp., C. synaphobranchi n. sp., C. verudaensis n. sp. and Myxodavisia bulani n. sp.; two of these formed a new, basal Ceratomyxa subclade. We identified that the Ceratomyxa clade is basal to all other marine myxosporean lineages, and Kudoa with Enteromyxum are the most recently branching clades. Topologies were least stable at the nodes connecting the marine urinary clade, the marine gall bladder clade and the Ceratonova clade. Bayesian inference analysis of SSU rDNA and the statistical tree topology tests suggested that Ceratonova is closely related to the Enteromyxum and Kudoa clades, which represent a large group of histozoic species. A close relationship between Ceratomyxa and Ceratonova was not supported, despite their similar myxospore morphologies. Overall, the site of sporulation in the vertebrate host is a more accurate predictor of phylogenetic relationships than the morphology of the myxospore.

Methods for analyzing the evolutionary relationship of NF-κB proteins using free, web-driven bioinformatics and phylogenetic tools

Phylogenetic analysis enables one to reconstruct the functional evolution of proteins. Current understanding of NF-κB signaling derives primarily from studies of a relatively small number of laboratory models-mainly vertebrates and insects-that represent a tiny fraction of animal evolution. As such, NF-κB has been the subject of limited phylogenetic analysis. The recent discovery of NF-κB proteins in "basal" marine animals (e.g., sponges, sea anemones, corals) and NF-κB-like proteins in non-metazoan lineages extends the origin of NF-κB signaling by several hundred million years and provides the opportunity to investigate the early evolution of this pathway using phylogenetic approaches. Here, we describe a combination of bioinformatic and phylogenetic analyses based on menu-driven, open-source computer programs that are readily accessible to molecular biologists without formal training in phylogenetic methods. These phylogenetically based comparisons of NF-κB proteins are powerful in that they reveal deep conservation and repeated instances of parallel evolution in the sequence and structure of NF-κB in distant animal groups, which suggest that important functional constraints limit the evolution of this protein.

Strepsiptera, phylogenomics and the long branch attraction problem

Insect phylogeny has recently been the focus of renewed interest as advances in sequencing techniques make it possible to rapidly generate large amounts of genomic or transcriptomic data for a species of interest. However, large numbers of markers are not sufficient to guarantee accurate phylogenetic reconstruction, and the choice of the model of sequence evolution as well as adequate taxonomic sampling are as important for phylogenomic studies as they are for single-gene phylogenies. Recently, the sequence of the genome of a strepsipteran has been published and used to place Strepsiptera as sister group to Coleoptera. However, this conclusion relied on a data set that did not include representatives of Neuropterida or of coleopteran lineages formerly proposed to be related to Strepsiptera. Furthermore, it did not use models that are robust against the long branch attraction artifact. Here we have sequenced the transcriptomes of seven key species to complete a data set comprising 36 species to study the higher level phylogeny of insects, with a particular focus on Neuropteroidea (Coleoptera, Strepsiptera, Neuropterida), especially on coleopteran taxa considered as potential close relatives of Strepsiptera. Using models robust against the long branch attraction artifact we find a highly resolved phylogeny that confirms the position of Strepsiptera as a sister group to Coleoptera, rather than as an internal clade of Coleoptera, and sheds new light onto the phylogeny of Neuropteroidea.

Evidence that ebolaviruses and cuevaviruses have been diverging from marburgviruses since the Miocene

An understanding of the timescale of evolution is critical for comparative virology but remains elusive for many RNA viruses. Age estimates based on mutation rates can severely underestimate divergences for ancient viral genes that are evolving under strong purifying selection. Paleoviral dating, however, can provide minimum age estimates for ancient divergence, but few orthologous paleoviruses are known within clades of extant viruses. For example, ebolaviruses and marburgviruses are well-studied mammalian pathogens, but their comparative biology is difficult to interpret because the existing estimates of divergence are controversial. Here we provide evidence that paleoviral elements of two genes (ebolavirus-like VP35 and NP) in cricetid rodent genomes originated after the divergence of ebolaviruses and cuevaviruses from marburgviruses. We provide evidence of orthology by identifying common paleoviral insertion sites among the rodent genomes. Our findings indicate that ebolaviruses and cuevaviruses have been diverging from marburgviruses since the early Miocene.

Maximum likelihood inference of small trees in the presence of long branches

The statistical basis of maximum likelihood (ML), its robustness, and the fact that it appears to suffer less from biases lead to it being one of the most popular methods for tree reconstruction. Despite its popularity, very few analytical solutions for ML exist, so biases suffered by ML are not well understood. One possible bias is long branch attraction (LBA), a regularly cited term generally used to describe a propensity for long branches to be joined together in estimated trees. Although initially mentioned in connection with inconsistency of parsimony, LBA has been claimed to affect all major phylogenetic reconstruction methods, including ML. Despite the widespread use of this term in the literature, exactly what LBA is and what may be causing it is poorly understood, even for simple evolutionary models and small model trees. Studies looking at LBA have focused on the effect of two long branches on tree reconstruction. However, to understand the effect of two long branches it is also important to understand the effect of just one long branch. If ML struggles to reconstruct one long branch, then this may have an impact on LBA. In this study, we look at the effect of one long branch on three-taxon tree reconstruction. We show that, counterintuitively, long branches are preferentially placed at the tips of the tree. This can be understood through the use of analytical solutions to the ML equation and distance matrix methods. We go on to look at the placement of two long branches on four-taxon trees, showing that there is no attraction between long branches, but that for extreme branch lengths long branches are joined together disproportionally often. These results illustrate that even small model trees are still interesting to help understand how ML phylogenetic reconstruction works, and that LBA is a complicated phenomenon that deserves further study.

EM for phylogenetic topology reconstruction on nonhomogeneous data

Background

The reconstruction of the phylogenetic tree topology of four taxa is, still nowadays, one of the main challenges in phylogenetics. Its difficulties lie in considering not too restrictive evolutionary models, and correctly dealing with the long-branch attraction problem. The correct reconstruction of 4-taxon trees is crucial for making quartet-based methods work and being able to recover large phylogenies.

Methods

We adapt the well known expectation-maximization algorithm to evolutionary Markov models on phylogenetic 4-taxon trees. We then use this algorithm to estimate the substitution parameters, compute the corresponding likelihood, and to infer the most likely quartet.

Results

In this paper we consider an expectation-maximization method for maximizing the likelihood of (time nonhomogeneous) evolutionary Markov models on trees. We study its success on reconstructing 4-taxon topologies and its performance as input method in quartet-based phylogenetic reconstruction methods such as QFIT and QuartetSuite. Our results show that the method proposed here outperforms neighbor-joining and the usual (time-homogeneous continuous-time) maximum likelihood methods on 4-leaved trees with among-lineage instantaneous rate heterogeneity, and perform similarly to usual continuous-time maximum-likelihood when data satisfies the assumptions of both methods.

Conclusions

The method presented in this paper is well suited for reconstructing the topology of any number of taxa via quartet-based methods and is highly accurate, specially regarding largely divergent trees and time nonhomogeneous data.

Sphaerospora sensu stricto: taxonomy, diversity and evolution of a unique lineage of myxosporeans (Myxozoa)

Myxosporeans (Myxozoa) are eukaryotic parasites, primarily of fish, whose classification is in a state of flux as taxonomists attempt to synthesize the traditional morphology-based system with emerging DNA sequence-based phylogenies. The genus Sphaerospora Thélohan, 1892, which includes pathogenic species that cause significant impacts on fisheries and aquaculture, is one of the most polyphyletic taxa and exemplifies the current challenges facing myxozoan taxonomists. The type species, S. elegans, clusters within the Sphaerospora sensu stricto clade, members of which share similar tissue tropism and long insertions in their variable rRNA gene regions. However, other morphologically similar sphaerosporids lie in different branches of myxozoan phylogenetic trees. Herein, we significantly extend taxonomic sampling of sphaerosporids with SSU+LSU rDNA and EF-2 sequence data for 12 taxa including three representatives of the morphologically similar genus Polysporoplasma Sitjà-Bobadilla et Álvarez-Pellitero, 1995. These taxa were sampled from different vertebrate host groups, biogeographic realms and environments. Our phylogenetic analyses and statistical tests of single and concatenated datasets revealed Sphaerospora s. s. as a strongly supported monophyletic lineage, that clustered sister to the whole myxosporean clade (freshwater+marine lineages). Generally, Sphaerospora s. s. rDNA sequences (up to 3.7 kb) are the longest of all myxozoans and indeed metazoans. The sphaerosporid clade has two lineages, which have specific morphological, biological and sequence traits. Lineage A taxa (marine Sphaerospora spp.) have a single binucleate sporoplasm and shorter AT-rich rDNA inserts. Lineage B taxa (freshwater/brackish Sphaerospora spp.+marine/brackish Polysporoplasma spp.) have 2-12 uninucleate sporoplasms and longer GC-rich rDNA inserts. Lineage B has four subclades that correlate with host group and habitat; all Polysporoplasma species, including the type species, cluster together in one of these subclades. We thus suppress the genus Polysporoplasma and the family Polysporoplasmidae and emend the generic diagnosis of the genus Sphaerospora. The combination of morphological, biological and DNA sequence data applied in this study helped to elucidate an important part of the taxonomic puzzle within the phylum Myxozoa.

Evaluating the Evolutionary Origins of Unexpected Character Distributions within the Bacterial Planctomycetes-Verrucomicrobia-Chlamydiae Superphylum

Recently, several characters that are absent from most bacteria, but which are found in many eukaryotes or archaea, have been identified within the bacterial Planctomycetes-Verrucomicrobia-Chlamydiae (PVC) superphylum. Hypotheses of the evolutionary history of such characters are commonly based on the inference of phylogenies of gene or protein families associated with the traits, estimated from multiple sequence alignments (MSAs). So far, studies of this kind have focused on the distribution of (i) two genes involved in the synthesis of sterol, (ii) tubulin genes, and (iii) c1 transfer genes. In many cases, these analyses have concluded that horizontal gene transfer (HGT) is likely to have played a role in shaping the taxonomic distribution of these gene families. In this article, we describe several issues with the inference of HGT from such analyses, in particular concerning the considerable uncertainty associated with our estimation of both gene family phylogenies (especially those containing ancient lineage divergences) and the Tree of Life (ToL), and the need for wider use and further development of explicit probabilistic models to compare hypotheses of vertical and horizontal genetic transmission. We suggest that data which is often taken as evidence for the occurrence of ancient HGT events may not be as convincing as is commonly described, and consideration of alternative theories is recommended. While focusing on analyses including PVCs, this discussion is also relevant for inferences of HGT involving other groups of organisms.

Long branch effects distort maximum likelihood phylogenies in simulations despite selection of the correct model

The aim of our study was to test the robustness and efficiency of maximum likelihood with respect to different long branch effects on multiple-taxon trees. We simulated data of different alignment lengths under two different 11-taxon trees and a broad range of different branch length conditions. The data were analyzed with the true model parameters as well as with estimated and incorrect assumptions about among-site rate variation. If length differences between connected branches strongly increase, tree inference with the correct likelihood model assumptions can fail. We found that incorporating invariant sites together with Γ distributed site rates in the tree reconstruction (Γ+I) increases the robustness of maximum likelihood in comparison with models using only Γ. The results show that for some topologies and branch lengths the reconstruction success of maximum likelihood under the correct model is still low for alignments with a length of 100,000 base positions. Altogether, the high confidence that is put in maximum likelihood trees is not always justified under certain tree shapes even if alignment lengths reach 100,000 base positions.