Topology-dependent asymmetry in systematic errors affects phylogenetic placement of Ctenophora and Xenacoelomorpha

The evolutionary relationships of two animal phyla, Ctenophora and Xenacoelomorpha, have proved highly contentious. Ctenophora have been proposed as the most distant relatives of all other animals (Ctenophora-first rather than the traditional Porifera-first). Xenacoelomorpha may be primitively simple relatives of all other bilaterally symmetrical animals (Nephrozoa) or simplified relatives of echinoderms and hemichordates (Xenambulacraria). In both cases, one of the alternative topologies must be a result of errors in tree reconstruction. Here, using empirical data and simulations, we show that the Ctenophora-first and Nephrozoa topologies (but not Porifera-first and Ambulacraria topologies) are strongly supported by analyses affected by systematic errors. Accommodating this finding suggests that empirical studies supporting Ctenophora-first and Nephrozoa trees are likely to be explained by systematic error. This would imply that the alternative Porifera-first and Xenambulacraria topologies, which are supported by analyses designed to minimize systematic error, are the most credible current alternatives

Systematic errors in phylogenetic trees

The effort to reconstruct the tree of life was revolutionized by the use of sequences of proteins and nucleic acids. Phylogenetic trees are now routinely inferred using hundreds of thousands of amino acid or nucleotide characters. It thus seems surprising that many aspects of the tree of life are still controversial; conflicting results between large scale phylogenomic studies show that errors remain common despite large datasets. These errors often result from systematic biases in the way sequences evolve. While the resulting systematic errors are well understood, it requires careful efforts to reduce their effects

Systematic errors in orthology inference and their effects on evolutionary analyses

Summary The availability of complete sets of genes from many organisms makes it possible to identify genes unique to (or lost from) certain clades. This information is used to reconstruct phylogenetic trees; identify genes involved in the evolution of clade specific novelties; and for phylostratigraphy—identifying ages of genes in a given species. These investigations rely on accurately predicted orthologs. Here we use simulation to produce sets of orthologs that experience no gains or losses. We show that errors in identifying orthologs increase with higher rates of evolution. We use the predicted sets of orthologs, with errors, to reconstruct phylogenetic trees; to count gains and losses; and for phylostratigraphy. Our simulated data, containing information only from errors in orthology prediction, closely recapitulate findings from empirical data. We suggest published downstream analyses must be informed to a large extent by errors in orthology prediction that mimic expected patterns of gene evolution

Automated, phylogeny-based genotype delimitation of the Hepatitis Viruses HBV and HCV

Background The classification of hepatitis viruses still predominantly relies on ad hoc criteria, i.e., phenotypic traits and arbitrary genetic distance thresholds. Given the subjectivity of such practices coupled with the constant sequencing of samples and discovery of new strains, this manual approach to virus classification becomes cumbersome and impossible to generalize. Methods Using two well-studied hepatitis virus datasets, HBV and HCV, we assess if computational methods for molecular species delimitation that are typically applied to barcoding biodiversity studies can also be successfully deployed for hepatitis virus classification. For comparison, we also used ABGD, a tool that in contrast to other distance methods attempts to automatically identify the barcoding gap using pairwise genetic distances for a set of aligned input sequences. Results—Discussion We found that the mPTP species delimitation tool identified even without adapting its default parameters taxonomic clusters that either correspond to the currently acknowledged genotypes or to known subdivision of genotypes (subtypes or subgenotypes). In the cases where the delimited cluster corresponded to subtype or subgenotype, there were previous concerns that their status may be underestimated. The clusters obtained from the ABGD analysis differed depending on the parameters used. However, under certain values the results were very similar to the taxonomy and mPTP which indicates the usefulness of distance based methods in virus taxonomy under appropriate parameter settings. The overlap of predicted clusters with taxonomically acknowledged genotypes implies that virus classification can be successfully automated

Codon-based analysis of selection pressure and genetic structure in the Psammobates tentorius (Bell, 1828) species complex, and phylogeny inferred from both codons and amino acid sequences

This study used codon analysis (dN/dS and Tv/Ti) to investigate selection pressure and genetic structure in the highly polymorphic Psammobates tentorius species complex, and amino acid sequences to construct a phylogeny tree for it. Our results revealed a strong selection signal at node ‘C2 + C3’, possibly driven by aridity intensification resulting from the development of the Benguela Current. A similar signal was noticed at C3, possibly due to the same driving force. These findings suggest that environmental selection pressure favoured those groups and that further cladogenic events were possible. Selection pressure was also found to be high at C1, C4 and C7, which may indicate that they are also favoured by the current selection pressure. The codon-based phylogeny did not retrieve any potentially undescribed species, but nonetheless provided support for the validity of the seven distinct clades retrieved with the DNA sequence data. The amino acid sequence-based phylogeny generally supported the seven lineages as valid putative species. Investigation at the genomic scale could, however, help to solve the issue. In general, we found the codon, dN, dS, Tv, Ti and amino acid sequence-based phylogenetic inferences useful in species delimitation and recommend their use in species delimitation studies

Is reproductive strategy a key factor in understanding the evolutionary history of Southern Ocean Asteroidea (Echinodermata)?

Life traits such as reproductive strategy can be determining factors of species evolutionary history and explain the resulting diversity patterns. This can be investigated using phylogeographic analyses of genetic units. In this work, the genetic structure of five asteroid genera with contrasting reproductive strategies (brooding: Diplasterias, Notasterias and Lysasterias versus broadcasting: Psilaster and Bathybiaster) was investigated in the Southern Ocean. Over 1,400 mtDNA cytochrome C oxidase subunit I (COI) sequences were analysed using five species delineation methods (ABGD, ASAP, mPTP, sGMYC and mGMYC), two phylogenetic reconstructions (ML and BA), and molecular clock calibrations, in order to examine the weight of reproductive strategy in the observed differences among phylogeographic patterns. We hypothesised that brooding species would show higher levels of genetic diversity and species richness along with a clearer geographic structuring than broadcasting species. In contrast, genetic diversity and species richness were not found to be significantly different between brooders and broadcasters, but broadcasters are less spatially structured than brooders supporting our initial hypothesis and suggesting more complex evolutionary histories associated to this reproductive strategy. Broadcasters' phylogeography can be explained by different scenarios including deep‐sea colonisation routes, bipolarity or cosmopolitanism, and sub‐Antarctic emergence for the genus Bathybiaster; Antarctic‐ New Zealand faunal exchanges across the Polar Front for the genus Psilaster. Brooders' phylogeography could support the previously formulated hypothesis of a past trans‐Antarctic seaway established between the Ross and the Weddell seas during the Plio‐Pleistocene. Our results also show, for the first time, that the Weddell Sea is populated by a mixed asteroid fauna originating from both the East and West Antarctic

Evolution and networks in ancient and widespread symbioses between Mucoromycotina and liverworts

Like the majority of land plants, liverworts regularly form intimate symbioses with arbuscular mycorrhizal fungi (Glomeromycotina). Recent phylogenetic and physiological studies report that they also form intimate symbioses with Mucoromycotina fungi and that some of these, like those involving Glomeromycotina, represent nutritional mutualisms. To compare these symbioses, we carried out a global analysis of Mucoromycotina fungi in liverworts and other plants using species delimitation, ancestral reconstruction, and network analyses. We found that Mucoromycotina are more common and diverse symbionts of liverworts than previously thought, globally distributed, ancestral, and often co-occur with Glomeromycotina within plants. However, our results also suggest that the associations formed by Mucoromycotina fungi are fundamentally different because, unlike Glomeromycotina, they may have evolved multiple times and their symbiotic networks are un-nested (i.e., not forming nested subsets of species). We infer that the global Mucoromycotina symbiosis is evolutionarily and ecologically distinctive

One thousand DNA barcodes of piranhas and pacus reveal geographic structure and unrecognised diversity in the Amazon

Piranhas and pacus (Characiformes: Serrasalmidae) are a charismatic but understudied family of Neotropical fishes. Here, we analyse a DNA barcode dataset comprising 1,122 specimens, 69 species, 16 genera, 208 localities, and 34 major river drainages in order to make an inventory of diversity and to highlight taxa and biogeographic areas worthy of further sampling effort and conservation protection. Using four methods of species discovery - incorporating both tree and distance based techniques - we report between 76 and 99 species-like clusters, i.e. between 20% and 33% of a priori identified taxonomic species were represented by more than one mtDNA lineage. There was a high degree of congruence between clusters, with 60% supported by three or four methods. Pacus of the genus Myloplus exhibited the most intraspecific variation, with six of the 13 species sampled found to have multiple lineages. Conversely, piranhas of the Serrasalmus rhombeus group proved difficult to delimit with these methods due to genetic similarity and polyphyly. Overall, our results recognise substantially underestimated diversity in the serrasalmids, and emphasise the Guiana and Brazilian Shield rivers as biogeographically important areas with multiple cases of across-shield and within-shield diversifications. We additionally highlight the distinctiveness and complex phylogeographic history of rheophilic taxa in particular, and suggest multiple colonisations of these habitats by different serrasalmid lineages. © 2018 The Author(s)

Peeking through the trapdoor: Historical biogeography of the Aegean endemic spider Cyrtocarenum Ausserer, 1871 with an estimation of mtDNA substitution rates for Mygalomorphae

The Aegean region, located in the Eastern Mediterranean, is an area of rich biodiversity and endemism. Its position, geographical configuration and complex geological history have shaped the diversification history of many animal taxa. Mygalomorph spiders have drawn the attention of researchers, as excellent model systems for phylogeographical investigations. However, phylogeographic studies of spiders in the Aegean region are scarce. In this study, we focused on the phylogeography of the endemic ctenizid trap-door spider Cyrtocarenum Ausserer, 1871. The genus includes two morphologically described species: C. grajum (C.L. Koch, 1836) and C. cunicularium (Olivier, 1811). We sampled 60 specimens from the distributions of both species and analyzed four mitochondrial and two nuclear markers. Cyrtocarenum served as an example to demonstrate the importance of natural history traits in the inference of phylogeographic scenarios. The mtDNA substitution rates inferred for the genus are profoundly higher compared to araneomorph spiders and other arthropods, which seems tightly associated with their biology. We evaluate published mtDNA substitution rates followed in the literature for mygalomorph spiders and discuss potential pitfalls. Following gene tree (maximum likelihood, Bayesian inference) and species tree approaches (*BEAST), we reconstructed a time-calibrated phylogeny of the genus. These results, combined with a biogeographical ancestral-area analysis, helped build a biogeographic scenario that describes how the major palaeogeographic and palaeoclimatic events of the Aegean may have affected the distribution of Cyrtocarenum lineages. The diversification of the genus seems to have begun in the Middle Miocene in the present west Aegean area, while major phylogenetic events occurred at the Miocene-Pliocene boundary for C. cunicularium, probably related to the Messinian Salinity Crisis. Our results also demonstrate the clear molecular distinction of the two morphologically described species, but possible cryptic lineages may exist within C. cunicularium. © 2016 Elsevier Inc

Molecular phylogeny of three Mesalina (Reptilia: Lacertidae) species (M. guttulata, M. brevirostris and M. bahaeldini) from North Africa and the Middle East: Another case of paraphyly?

Mesalina is a widespread lacertid genus occurring throughout the Saharo-Sindian region from North Africa to Pakistan. It has been through a series of taxonomic revisions, but the phylogenetic relationships among the species remain unclear. In this study we estimate the phylogeographic structure of M. guttulata across most of its distributional range and we evaluate the relationships between M. guttulata and the sympatric species M. brevirostris and M. bahaeldini using partial mitochondrial DNA (mtDNA) sequences (cyt b and 16S). M. guttulata and M. brevirostris represent species complexes, whereas M. bahaeldini considered before as M. guttulata is a recently described species with very restricted distribution. Here we present the first evidence that M. guttulata is a paraphyletic taxon with respect to M. bahaeldini, while M. brevirostris proves to be a polytypic species or even a species complex, confirming previous studies. Although mtDNA markers have several properties that make them suitable for phylogeographic studies, they are not free of difficulties. Phylogeographic inferences within and between closely related species can be mislead by introgression and retention of ancestral polymorphism (incomplete lineage sorting). However, the present distribution pattern, the estimated times of divergence and the significant variation in morphology within M. guttulata led us to accept that the paraphyletic pattern observed, is most likely due to inaccurate taxonomy. Our hypothesis is that what has hitherto been considered as intraspecific variation, actually reflects species-level variation. Furthermore, our biogeographic analyses and the estimated time of divergences suggest that the present distribution of M. guttulata was the result of several dispersal and vicariant events, which are associated with historical changes (climatic oscillations and paleogeographic barriers) of late Miocene and Pliocene period. © 2008 Elsevier Inc. All rights reserved