Search

Wintersweet (Chimonanthus praecox) is one of the most important ornamental plants, whose color is mainly determined by the middle tepals. However, the molecular mechanisms underlying the intriguing flower color development among different wintersweet groups are still largely unknown. In addition, wintersweet belongs to magnoliids, and the issue remains as to how to determine the phylogenetic position of magnoliids conclusively. Here, the whole genome of red flower wintersweet, a new wintersweet type, was sequenced and assembled with high-quality. The genome comprised 11 super-scaffolds (chromosomes) with a total size of 737.03 Mb. Based on the analyses of the long branch attraction (LBA), incomplete lineage sorting (ILS), sparse taxon sampling and other factors, we suggest that a bifurcating tree may not fully represent the complex early diversification of the angiosperms, and that magnoliids are most likely to be the sister of the eudicots. The wintersweet genome appears to have undergone two rounds of whole-genome duplication (WGD) events: a younger WGD event representing an independent event specific to the Calycanthaceae species; and an ancient WGD event shared by Laurales. By integrating genomic, transcriptomic, and metabolomic data, CpANS1 and transcription factor CpMYB1 were considered to play key roles in regulating tepal color development, whereas CpMYB1 needs to form a complex with bHLH and WD40 to fully perform its regulatory function. The present study not only provided novel insights into the evolution of magnoliids and the molecular mechanism for flower color development, but also laid the foundation for subsequent functional genomics study and molecular breeding of wintersweet.

Phylogenetic analyses may suffer from multiple sources of error leading to conflict between genes and methods of inference. The evolutionary history of the mollusc clade Vetigastropoda makes them susceptible to these conflicts, their higher level phylogeny remaining largely unresolved. Originating over 350 Ma, vetigastropods were the dominant marine snails in the Paleozoic. Multiple extinction events and new radiations have resulted in both very long and very short branches and a large extant diversity of over 4000 species. This is the perfect setting of a hard phylogenetic question in which sources of conflict can be explored. We present 41 new transcriptomes across the diversity of vetigastropods (62 terminals total), and provide the first genomic-scale phylogeny for the group. We find that deep divergences differ from previous studies in which long branch attraction was likely pervasive. Robust results leading to changes in taxonomy include the paraphyly of the order Lepetellida and the family Tegulidae. Tectinae subfam. nov. is designated for the clade comprising Tectus, Cittarium, and Rochia. For two early divergences, topologies disagreed between concatenated analyses using site heterogeneous models versus concatenated partitioned analyses and summary coalescent methods. We investigated rate and composition heterogeneity among genes, as well as missing data by locus and by taxon, none of which had an impact on the inferred topologies. We also found no evidence for ancient introgression throughout the phylogeny. We further tested whether uninformative genes and over-partitioning were responsible for this discordance by evaluating the phylogenetic signal of individual genes using likelihood mapping, and by analyzing the most informative genes with a full multispecies coalescent (MSC) model. We find that most genes are not informative at the two conflicting nodes, but neither this nor gene-wise partitioning are the cause of discordant results. New method implementations that simultaneously integrate amino acid profile mixture models and the MSC might be necessary to resolve these and other recalcitrant nodes in the Tree of Life. [Fissurellidae; Haliotidae; likelihood mapping; multispecies coalescent; phylogenetic signal; phylogenomic conflict; site heterogeneity; Trochoidea.]

Long branch attraction (LBA) is a prevalent form of bias in phylogenetic estimation but the reasons for it are only partially understood. We argue here that it is largely due to differences in the sizes of the model spaces corresponding to different trees. Trees with long branches together allow much more flexible internal branch length parameter estimation. Consequently, although each tree has the same number of parameters, trees with long branches together have larger effective model spaces. The problem of LBA becomes particularly pronounced with partitioned data. Formulation of tree estimation as model selection leads us to propose bootstrap bias corrections as cross-checks on estimation when long branches end up being estimated together. [Bootstrap; long branch attraction; maximum likelihood; model selection; partitioned model; phylogenetics.]

Long-branch attraction is a systematic artifact that results in erroneous groupings of fast-evolving taxa. The combination of short, deep internodes in tandem with long-branch attraction artifacts has produced empirically intractable parts of the Tree of Life. One such group is the arthropod subphylum Chelicerata, whose backbone phylogeny has remained unstable despite improvements in phylogenetic methods and genome-scale data sets. Pseudoscorpion placement is particularly variable across data sets and analytical frameworks, with this group either clustering with other long-branch orders or with Arachnopulmonata (scorpions and tetrapulmonates). To surmount long-branch attraction, we investigated the effect of taxonomic sampling via sequential deletion of basally branching pseudoscorpion superfamilies, as well as varying gene occupancy thresholds in supermatrices. We show that concatenated supermatrices and coalescent-based summary species tree approaches support a sister group relationship of pseudoscorpions and scorpions, when more of the basally branching taxa are sampled. Matrix completeness had demonstrably less influence on tree topology. As an external arbiter of phylogenetic placement, we leveraged the recent discovery of an ancient genome duplication in the common ancestor of Arachnopulmonata as a litmus test for competing hypotheses of pseudoscorpion relationships. We generated a high-quality developmental transcriptome and the first genome for pseudoscorpions to assess the incidence of arachnopulmonate-specific duplications (e.g., homeobox genes and miRNAs). Our results support the inclusion of pseudoscorpions in Arachnopulmonata (new definition), as the sister group of scorpions. Panscorpiones (new name) is proposed for the clade uniting Scorpiones and Pseudoscorpiones.

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

Fabales is a cosmopolitan angiosperm order that consists of four families: Leguminosae (Fabaceae), Polygalaceae, Surianaceae, and Quillajaceae. Despite the great interest in this group, a convincing phylogeny of the order is still not available. Therefore, the aim of this study was to explicitly test for possible long branch attraction (LBA) problems within Fabales for the first time, and determine whether low tree stemminess and unequal branch lengths could worsen this problem. Supermatrix analysis of Fabales was carried out using previously published plastid matK, trnL, rbcL, and newly sequenced nuclear sqd1 regions for 678 taxa in total, including 43 outgroup taxa from families of Fabidae. We employed additional analyses, such as simulations, network analyses, sampling different outgroup taxa (random or real), removing fast evolving sites and fast evolving taxa, and molecular clock rooting, to identify both LBA and (or) rooting problems. These analyses clearly show that the Fabales phylogeny has been influenced by the sampling of outgroup taxa, but not LBA. However, network analyses show that even though it is weak, there is a consistent phylogenetic signal among the rapidly radiated Fabales families, which can be traced by further analyses. While, molecular clock rooting analysis yielded a (Leguminosae(Polygalaceae(Surianaceae+Quillajaceae))) topology with strong support for the first time here, supermatrix analyses yielded a ((Leguminosae+Polygalaceae)(Surianaceae+Quillajaceae)) with low-moderate support. © 2020, Canadian Science Publishing. All rights reserved.

Due to morphological reduction and absence of amplifiable plastid genes, the identification of photosynthetic relatives of heterotrophic plants is problematic. Although nuclear and mitochondrial gene sequences may offer a welcome alternative source of phylogenetic markers, the presence of rate heterogeneity in these genes may introduce bias/systematic error in phylogenetic analyses. We examine the phylogenetic position of Thismiaceae based on nuclear 18S rDNA and mitochondrial atpA DNA sequence data, as well as using parsimony, likelihood and Bayesian inference methods. Significant differences in evolutionary rates of these genes between closely related taxa lead to conflicting results: while parsimony analyses of 18S rDNA and combined data strongly support the monophyly of Thismiaceae, Bayesian inference, with and without a relaxed molecular clock, as well as the Swofford-Olsen-Waddell-Hillis (SOWH) test confidently reject this hypothesis. We show that rate heterogeneity in our data leads to long-branch attraction artifacts in parsimony analysis. However, using model-based inference methods the question of whether Thismiaceae are monophyletic remains elusive. On the one hand maximum likelihood nonparametric bootstrapping and parametric hypothesis tests fail to support a paraphyletic Thismiaceae, on the other hand Bayesian inference methods (both without and with a relaxed clock) significantly reject a monophyletic Thismiaceae. These results show that an adequate sampling, the use of rate homogeneous data, and the application of different inference methods are important factors for developing phylogenetic hypotheses of myco-heterotrophic plants. (C) The Willi Hennig Society 2009. ispartof: Cladistics vol:25 issue:1 pages:64-77 ispartof: location:United States status: published

Background Seisonidea (also Seisonacea or Seisonidae) is a group of small animals living on marine crustaceans (Nebalia spec.) with only four species described so far. Its monophyletic origin with mostly free-living wheel animals (Monogononta, Bdelloidea) and endoparasitic thorny-headed worms (Acanthocephala) is widely accepted. However, the phylogenetic relationships inside the Rotifera-Acanthocephala clade (Rotifera sensu lato or Syndermata) are subject to ongoing debate, with consequences for our understanding of how genomes and lifestyles might have evolved. To gain new insights, we analyzed first drafts of the genome and transcriptome of the key taxon Seisonidea. Results Analyses of gDNA-Seq and mRNA-Seq data uncovered two genetically distinct lineages in Seison nebaliae Grube, 1861 off the French Channel coast. Their mitochondrial haplotypes shared only 82% sequence identity despite identical gene order. In the nuclear genome, distinct linages were reflected in different gene compactness, GC content and codon usage. The haploid nuclear genome spans ca. 46 Mb, of which 96% were reconstructed. According to ~ 23,000 SuperTranscripts, gene number in S. nebaliae should be within the range published for other members of Rotifera-Acanthocephala. Consistent with this, numbers of metazoan core orthologues and ANTP-type transcriptional regulatory genes in the S. nebaliae genome assembly were between the corresponding numbers in the other assemblies analyzed. We additionally provide evidence that a basal branching of Seisonidea within Rotifera-Acanthocephala could reflect attraction to the outgroup. Accordingly, rooting via a reconstructed ancestral sequence led to monophyletic Pararotatoria (Seisonidea+Acanthocephala) within Hemirotifera (Bdelloidea+Pararotatoria). Conclusion Matching genome/transcriptome metrics with the above phylogenetic hypothesis suggests that a haploid nuclear genome of about 50 Mb represents the plesiomorphic state for Rotifera-Acanthocephala. Smaller genome size in S. nebaliae probably results from subsequent reduction. In contrast, genome size should have increased independently in monogononts as well as bdelloid and acanthocephalan stem lines. The present data additionally indicate a decrease in gene repertoire from free-living to epizoic and endoparasitic lifestyles. Potentially, this reflects corresponding steps from the root of Rotifera-Acanthocephala via the last common ancestors of Hemirotifera and Pararotatoria to the one of Acanthocephala. Lastly, rooting via a reconstructed ancestral sequence may prove useful in phylogenetic analyses of other deep splits. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07857-y.

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.

Summary Species of the family Tylenchidae are encountered in large numbers in soils. The genus Labrys was recently described characterised by a remarkable lip pattern that differs from all other known Tylenchidae genera. Here we describe a curious new species, Labrys fujianensis sp. n., that morphologically fits the genus Labrys but which is genetically divergent. The phylogeny was inferred based on 18S and 28S rDNA and light and scanning electron microscopy were used to extract detailed morphologies. The phylogenetic position of this species and its phenotypic convergence are discussed. The possibility of a long-branch attraction artefact was inspected both by removal of variable nucleotide sites and monophyletic testing of topologies. The results confirmed the divergent positioning of the presented species and it is demonstrated that the genetic diversity in Tylenchidae may be much higher than expected due to morphological homoplasy.

It is commonly assumed that increasing the number of characters has the potential to resolve evolutionary radiations. Here, we studied photosynthetic stramenopiles (Ochrophyta) using alignments of heterogeneous origin mitochondrion, plastid, and nucleus. Surprisingly while statistical support for the relationships between the six major Ochrophyta lineages increases when comparing the mitochondrion (6,762 sites) and plastid (21,692 sites) trees, it decreases in the nuclear (209,105 sites) tree. Statistical support is not simply related to the data set size but also to the quantity of phylogenetic signal available at each position and our ability to extract it. Here, we show that this ability for current phylogenetic methods is limited, because conflicting results were obtained when varying taxon sampling. Even though the use of a better fitting model improved signal extraction and reduced the observed conflicts, the plastid data set provided higher statistical support for the ochrophyte radiation than the larger nucleus data set. We propose that the higher support observed in the plastid tree is due to an acceleration of the evolutionary rate in one short deep internal branch, implying that more phylogenetic signal per position is available to resolve the Ochrophyta radiation in the plastid than in the nuclear data set. Our work therefore suggests that, in order to resolve radiations, beyond the obvious use of data sets with more positions, we need to continue developing models of sequence evolution that better extract the phylogenetic signal and design methods to search for genes/characters that contain more signal specifically for short internal branches.

International audience; Phylogenomics seeks to use next‐generation data to robustly infer an organism's evolutionary history. Yet, the practical caveats of phylogenomics motivate investigation of improved efficiency, particularly when quality of phylogenies are questionable. To achieve improvements, one goal is to maintain or enhance the quality of phylogenetic inference while severely reducing dataset size. We approach this by assessing which kinds of loci in phylogenomic alignments provide the majority of support for a phylogenetic inference of cockroaches in Blaberoidea. We examine locus substitution rate, saturation, evolutionary divergence, rate heterogeneity, stabilizing selection, and a priori information content as traits that may determine optimality. Our controlled experimental design is based on 265 loci for 102 blaberoidean taxa and 22 outgroup species. Loci with high substitution rate, low saturation, low sequence distance, low rate heterogeneity, and strong stabilizing selection derive more support for phylogenetic relationships. We found that some phylogenetic information content estimators may not be meaningful for assessing information content a priori. We use these findings to design concatenated datasets with an optimized subsample of 100 loci. The tree inferred from the optimized subsample alignment was largely identical to that inferred from all 265 loci but with less evidence of long branch attraction, improved statistical support, and potential 4‐6x improvements to computation time. Supported by phylogenetic and morphological evidence, we erect three newly named clades (Anallactinae Evangelista & Wipfler subfam. nov., Orkrasomeria tax. nov. Evangelista, Wipfler, & Béthoux and Hemithyrsocerini Evangelista tribe nov.) and propose other taxonomic modifications. The diagnosis of Pseudophyllodromiidae Grandcolas, 1996 is modified to accommodate Anallactinae and Pseudophyllodromiinae Vickery & Kevan, 1983. The diagnosis of Ectobiidae Brunner von Wattenwyl, 1865 is modified to add novel morphological characters.

Chaetonotidae is the most diverse and widely distributed family of the order Chaetonotida (Gastrotricha) and includes both marine and freshwater species. Although the family is regarded as a sister taxon to the exclusively marine Xenotrichulidae, the type of environment, marine or freshwater, where Chaetonotidae originated is still not known. Here, we reconstructed the phylogeny of the family based on molecular sequence data and mapped both morphological and ecological characters to determine the ancestral environment of the first members of the family. Our results revealed that the freshwater genus Bifidochaetus is the earliest branching lineage in the paraphyletic Chaetonotidae (encompassing Dasydytidae and Neogosseidae). Moreover, we reconstructed Lepidochaetus-Cephalionotus clade as a monophyletic sister group to the remaining chaetonotids, which supports Kisielewski's morphological based hypothesis concerning undifferentiated type of body scales as a most primary character in Chaetonotidae. We also found that reversals to marine habitats occurred independently in different Chaetonotidae lineages, thus marine species in the genera Heterolepidoderma, Halichaetonotus, Aspidiophorus and subgenera Chaetonotus (Schizochaetonotus) or Chaetonotus (Marinochaetus) should be assumed as having secondarily invaded the marine environment. Character mapping revealed a series of synapomorphies that define the clade that includes Chaetonotidae (with Dasydytidae and Neogosseidae), the most important of which may be those linked to reproduction.

Despite their ecological and evolutionary importance as key components of tropical ecosystems, the phylogeny of fig trees is still unresolved. We use restriction-site-associated DNA (RAD) sequencing (ca420kb) and 102 morphological characters to elucidate the relationships between 70 species ofFicusrepresenting all known subgenera and sections and five outgroups. We compare morphological and molecular results to highlight discrepancies and reveal possible inference bias. We analyse marker and taxon properties that may bias molecular inferences, with existing softwares and a new approach based on iterative principal component analysis to reduce variance between clusters of samples. For the first time, with both molecular and morphological data, we recover a monophyletic subgenusUrostigmaand a clade with all gynodioecious fig trees. However, our analyses show that it is not possible to homogenize evolutionary rates and GC content for all taxa prior to phylogenetic inference and that four competing positions for the root of the molecular tree are possible. The placement of the long-branched sectionPharmacosyceaas sister to all other fig trees is not supported by morphological data and considered as a result of a long branch attraction artefact to the outgroups. Regarding morphological features and indirect evidence from the pollinator tree of life, the topology that divides the genusFicusinto monoeciousversusgynodioecious species appears most likely. Active pollination is inferred as the ancestral state for all topologies, ambiguity remains for ancestral breeding system including for the favored topology, and it appears most likely that the ancestor of fig trees was a freestanding tree. Increasing sampling may improve results and would be at least as relevant as maximizing the number of sequenced regions given the strong heterogeneity in evolutionary rates, and to a lesser extent, base composition among species. Despite morphological plasticity and frequent homoplasy of multiple characters, we advocate giving a central role to morphology in our understanding of the evolution ofFicus, especially as it can help detect insidious systematic errors that tend to become more pronounced with larger molecular data sets.

The majority strand of mitochondrial genomes of crustaceans usually exhibits negative GC skews. Most isopods exhibit an inversed strand asymmetry, believed to be a consequence of an inversion of the replication origin (ROI). Recently, we proposed that an additional ROI event in the common ancestor of Cymothoidae and Corallanidae families resulted in a double-inverted skew (negative GC), and that taxa with homoplastic skews cluster together in phylogenetic analyses (long-branch attraction, LBA). Herein, we further explore these hypotheses, for which we sequenced the mitogenome of Asotana magnifica (Cymothoidae), and tested whether our conclusions were biased by poor taxon sampling and inclusion of outgroups. (1) The new mitogenome also exhibits a double-inverted skew, which supports the hypothesis of an additional ROI event in the common ancestor of Cymothoidae and Corallanidae families. (2) It exhibits a unique gene order, which corroborates that isopods possess exceptionally destabilized mitogenomic architecture. (3) Improved taxonomic sampling failed to resolve skew-driven phylogenetic artefacts. (4) The use of a single outgroup exacerbated the LBA, whereas both the use of a large number of outgroups and complete exclusion of outgroups ameliorated it.

Resolving deep divergences in the tree of life is challenging even for analyses of genome-scale phylogenetic data sets. Relationships between Basidiomycota subphyla, the rusts and allies (Pucciniomycotina), smuts and allies (Ustilaginomycotina), and mushroom-forming fungi and allies (Agaricomycotina) were found particularly recalcitrant both to traditional multigene and genome-scale phylogenetics. Here, we address basal Basidiomycota relationships using concatenated and gene tree-based analyses of various phylogenomic data sets to examine the contribution of several potential sources of bias. We evaluate the contribution of biological causes (hard polytomy, incomplete lineage sorting) versus unmodeled evolutionary processes and factors that exacerbate their effects (e.g., fast-evolving sites and long-branch taxa) to inferences of basal Basidiomycota relationships. Bayesian Markov Chain Monte Carlo and likelihood mapping analyses reject the hard polytomy with confidence. In concatenated analyses, fast-evolving sites and oversimplified models of amino acid substitution favored the grouping of smuts with mushroom-forming fungi, often leading to maximal bootstrap support in both concatenation and coalescent analyses. On the contrary, the most conserved data subsets grouped rusts and allies with mushroom-forming fungi, although this relationship proved labile, sensitive to model choice, to different data subsets and to missing data. Excluding putative long-branch taxa, genes with high proportions of missing data and/or with strong signal failed to reveal a consistent trend toward one or the other topology, suggesting that additional sources of conflict are at play. While concatenated analyses yielded strong but conflicting support, individual gene trees mostly provided poor support for any resolution of rusts, smuts, and mushroom-forming fungi, suggesting that the true Basidiomycota tree might be in a part of tree space that is difficult to access using both concatenation and gene tree-based approaches. Inference-based assessments of absolute model fit strongly reject best-fit models for the vast majority of genes, indicating a poor fit of even the most commonly used models. While this is consistent with previous assessments of site-homogenous models of amino acid evolution, this does not appear to be the sole source of confounding signal. Our analyses suggest that topologies uniting smuts with mushroom-forming fungi can arise as a result of inappropriate modeling of amino acid sites that might be prone to systematic bias. We speculate that improved models of sequence evolution could shed more light on basal splits in the Basidiomycota, which, for now, remain unresolved despite the use of whole genome data.

Maximum likelihood and maximum parsimony are two key methods for phylogenetic tree reconstruction. Under certain conditions, each of these two methods can perform more or less efficiently, resulting in unresolved or disputed phylogenies. We show that a neural network can distinguish between four-taxon alignments that were evolved under conditions susceptible to either long-branch attraction or long-branch repulsion. When likelihood and parsimony methods are discordant, the neural network can provide insight as to which tree reconstruction method is best suited to the alignment. When applied to the contentious case of Strepsiptera evolution, our method shows robust support for the current scientific view, that is, it places Strepsiptera with beetles, distant from flies.

Assessing support for molecular phylogenies is difficult because the data is heterogeneous in quality and overwhelming in quantity. Traditionally, node support values (bootstrap frequency, Bayesian posterior probability) are used to assess confidence in tree topologies. Other analyses to assess the quality of phylogenetic data (e.g. Lento plots, saturation plots, trait consistency) and the resulting phylogenetic trees (e.g. internode certainty, parameter permutation tests, topological tests) exist but are rarely applied. Here we argue that a single qualitative analysis is insufficient to assess support of a phylogenetic hypothesis and relate data quality to tree quality. We use six molecular markers to infer the phylogeny of Blattodea and apply various tests to assess relationship support, locus quality, and the relationship between the two. We use internode-certainty calculations in conjunction with bootstrap scores, alignment permutations, and an approximately unbiased (AU) test to assess if the molecular data unambiguously support the phylogenetic relationships found. Our results show higher support for the position of Lamproblattidae, high support for the termite phylogeny, and low support for the position of Anaplectidae, Corydioidea and phylogeny of Blaberoidea. We use Lento plots in conjunction with mutation-saturation plots, calculations of locus homoplasy to assess locus quality, identify long branch attraction, and decide if the tree's relationships are the result of data biases. We conclude that multiple tests and metrics need to be taken into account to assess tree support and data robustness.

Catfishes (Teleostei: Siluriformes) are a highly diverse order within Ostariophysi that is distributed worldwide. At the base of this clade emerge three lineages with well-defined monophylies: Diplomystidae, Loricarioidei, and Siluroidei. Morphological phylogeny studies place the Diplomystidae as the earliest branching of these three lineages, but studies based on molecular phylogenetics consistently find the fast-evolving Loricarioidei instead. The high lineage evolutionary rate heterogeneity in this order and the fact that the lineage placed closest to the root in the molecular phylogenies is fast evolving, including many long branches, raises the possibility that the discrepancy between morphological and molecular phylogenies may be the result of a long branch attraction inference artifact. We test this hypothesis by using a 10-gene dataset to evaluate the arrangement of the three main siluriform lineages, and apply the LS3 and LS4 taxon sequence subsampling methods to reduce evolutionary rate heterogeneity among lineages. The initial and complete dataset supports the basal branching of Loricarioidei as in all previous molecular phylogenies, but once lineage rate heterogeneity is reduced with LS3 or LS4 through the removal of sequences disrupting homogeneity, the phylogeny shows Diplomystidae as the earliest branching group, with high supports, as proposed by morphological phylogeny. The result obtained with LS3, however, introduces the misplacement of one of the species with the highest amount of missing data, Scoloplax sp. Because the sequence sub-selection criterion of LS4 has been optimized to reduce data removal, the phylogeny resulting from the LS4-processed data is in agreement with the known intra-lineage relationships in addition to supporting the morphologically-based rooting hypothesis. Our results are the first instance in which a consensus between molecular and morphological phylogeny is reached concerning the root of this order.

Background Many algorithms and programs are available for phylogenetic reconstruction of families of proteins. Methods used widely at present use either a number of distance-based principles or character-based principles of maximum parsimony or maximum likelihood. Results We developed a novel program, named PQ, for reconstructing protein and nucleic acid phylogenies following a new character-based principle. Being tested on natural sequences PQ improves upon the results of maximum parsimony and maximum likelihood. Working with alignments of 10 and 15 sequences, it also outperforms the FastME program, which is based on one of the distance-based principles. Among all tested programs PQ is proved to be the least susceptible to long branch attraction. FastME outperforms PQ when processing alignments of 45 sequences, however. We confirm a recent result that on natural sequences FastME outperforms maximum parsimony and maximum likelihood. At the same time, both PQ and FastME are inferior to maximum parsimony and maximum likelihood on simulated sequences. PQ is open source and available to the public via an online interface. Conclusions The software we developed offers an open-source alternative for phylogenetic reconstruction for relatively small sets of proteins and nucleic acids, with up to a few tens of sequences. Electronic supplementary material The online version of this article (10.1186/s12859-018-2399-4) contains supplementary material, which is available to authorized users.

The Guianas are one of the most diverse regions of the Neotropics, hosting a particularly high rate of freshwater fish endemism. The present distributional patterns of freshwater fish species in the major catchments of the Guianas (comprising Guyana, Suriname and French Guiana) were analysed to reveal the faunal relationships between rivers, evaluate different hypotheses concerning biogeographical units, and redefine the boundaries of the Guianese freshwater ecoregions. A parsimony analysis of endemicity was performed using a data partitioning strategy to alleviate some drawbacks inherent to the method (e.g. long branch attraction artefact, heterotachy), and take into account alternative parsimony models assigning different constraints on state changes for the different species. A strong spatial element was present in the data with a structuring of species along a west-east gradient. Two main biogeographical units were highlighted: one to the west, ranging from the Essequibo to the Commewijne rivers and including the Proto Berbice and Surinamese regions, and one to the east ranging from the Maroni to the Oyapock rivers and including the Western, Central and Eastern French Guiana regions. Each ecoregion possessed distinctive fish assemblages, and three to four potential zones of faunal exchanges between Amazonian and Guianese rivers have been confirmed.

Premise of the study Phylogenomic studies employing large numbers of genes, including those based on plastid genomes (plastomes), are becoming common. Nonphotosynthetic plants such as mycoheterotrophs (which rely on root-associated fungi for essential nutrients, including carbon) tend to have highly elevated rates of plastome evolution, substantial genome reduction, or both. Mycoheterotroph plastomes therefore provide excellent test cases for investigating how extreme conditions impact phylogenomic inference. Methods We used parsimony and likelihood analysis of protein-coding gene sets from published and newly completed plastomes to infer the phylogenetic placement of taxa from the 10 angiosperm families in which mycoheterotrophy evolved. Key results Despite multiple very long branches that reflect elevated substitution rates, and frequently patchy gene recovery due to genome reduction, inferred phylogenetic placements of most mycoheterotrophic lineages in DNA-based likelihood analyses are both well supported and congruent with other studies. Amino-acid-based likelihood placements are broadly consistent with DNA-based inferences, but extremely rate-elevated taxa can have unexpected placements-albeit with weak support. In contrast, parsimony analysis is strongly misled by long-branch attraction among many distantly related mycoheterotrophic monocots. Conclusions Mycoheterotrophic plastomes provide challenging cases for phylogenomic inference, as substitutional rates can be elevated and genome reduction can lead to sparse gene recovery. Nonetheless, diverse likelihood frameworks provide generally well-supported and mutually concordant phylogenetic placements of mycoheterotrophs, consistent with recent phylogenetic studies and angiosperm-wide classifications. Previous predictions of parallel photosynthesis loss within families are supported for Burmanniaceae, Ericaceae, Gentianaceae, and Orchidaceae. Burmanniaceae and Thismiaceae should not be combined as a single family in Dioscoreales.

Phytophagous stink bugs are globally distributed and many harbor vertically inherited bacterial symbionts that are extracellular, yet little is known about how the symbiont’s genomes have evolved under this transmission strategy. Genome reduction is common in insect intracellular symbionts but limited genome sampling of the extracellular symbionts of distantly related stink bugs has precluded inferring patterns of extracellular symbiont genome evolution. To address this knowledge gap, we completely sequenced the genomes of the uncultivable bacterial symbionts of four neotropical stink bugs of the Edessa genus. Phylogenetic and comparative analyses indicated that the symbionts form a clade within the Pantoea genus and their genomes are highly reduced (∼0.8 Mb). Furthermore, genome synteny analysis and a jackknife approach for phylogenetic reconstruction, which corrected for long branch attraction artifacts, indicated that the Edessa symbionts were the result of a single symbiotic event that was distinct from the symbiosis event giving rise to Candidatus “Pantoea carbekii,” the extracellular symbiont of the invasive pentatomid stink bug, Halyomorpha halys. Metabolic functions inferred from the Edessa symbiont genomes suggests a shift in genomic composition characteristic of its lifestyle in that they retained many host-supportive functions while undergoing dramatic gene loss and establishing a stable relationship with their host insects. Given the undersampled nature of extracellular insect symbionts, this study is the first comparative analysis of these symbiont genomes from four distinct Edessa stink bug species. Finally, we propose the candidate name “Candidatus Pantoea edessiphila” for the species of these symbionts with strain designations according to their host species.

Inversions of the origin of replication (ORI) in mitochondrial genomes produce asymmetrical mutational pressures that can cause strong base composition skews. Due to skews often being overlooked, the total number of crustacean lineages that underwent ORI events remains unknown. We analysed skews, cumulative skew plots, conserved sequence motifs, and mitochondrial architecture of all 965 available crustacean mitogenomes (699 unique species). We found indications of an ORI in 159 (22.7%) species, and mapped these to 23 ORI events: 16 identified with confidence and 7 putative (13 newly proposed, and for 5 we improved the resolution). Two ORIs occurred at or above the order level: Isopoda and Copepoda. Shifts in skew plots are not a precise tool for identifying the replication mechanism. We discuss how ORIs can produce mutational bursts in mitogenomes and show how these can interfere with various types of evolutionary studies. Phylogenetic analyses were plagued by artefactual clustering, and ORI lineages exhibited longer branches, a higher number of synonymous substitutions, higher mutational saturation, and higher compositional heterogeneity. ORI events also affected codon usage and protein properties. We discuss how this may have caused erroneous interpretation of data in previous studies that did not account for skew patterns.