Phylogenetic Inference Methods Research Articles

An automated convergence diagnostic for phylogenetic MCMC analyses.

Assessing convergence of Markov chain Monte Carlo (MCMC) based analyses is crucial but challenging, especially so in high dimensional and complex spaces such as the space of phylogenetic trees (treespace). In practice, it is assumed that the target distribution is the unique stationary distribution of the MCMC and convergence is achieved when samples appear to be stationary. Here we leverage recent advances in computational geometry of the treespace and introduce a method that combines classical statistical techniques and algorithms with geometric properties of the treespace to automatically evaluate and assess practical convergence of phylogenetic MCMC analyses. Our method monitors convergence across multiple MCMC chains and achieves high accuracy in detecting both practical convergence and convergence issues within treespace. Furthermore, our approach is developed to allow for real-time evaluation during the MCMC algorithm run, eliminating any of the chain post-processing steps that are currently required. Our tool therefore improves reliability and efficiency of MCMC based phylogenetic inference methods and makes analyses easier to reproduce and compare. We demonstrate the efficacy of our diagnostic via a well-calibrated simulation study and provide examples of its performance on real data sets. Although our method performs well in practice, a significant part of the underlying treespace probability theory is still missing, which creates an excellent opportunity for future mathematical research in this area. The open source package for the phylogenetic inference framework BEAST2, called ASM, that implements these methods, making them accessible through a user-friendly GUI, is available from https://github.com/rbouckaert/asm/. The open source Python package, called tetres, that provides an interface for these methods enabling their applications beyond BEAST2 can be accessed at https://github.com/bioDS/tetres/.

CMAPLE: efficient phylogenetic inference in the pandemic era.

We have recently introduced MAPLE (MAximum Parsimonious Likelihood Estimation), a new pandemic-scale phylogenetic inference method exclusively designed for genomic epidemiology. In response to the need for enhancing MAPLE's performance and scalability, here we present two key components: (1) CMAPLE software, a highly optimized C++ reimplementation of MAPLE with many new features and advancements; and (2) CMAPLE library, a suite of Application Programming Interfaces to facilitate the integration of the CMAPLE algorithm into existing phylogenetic inference packages. Notably, we have successfully integrated CMAPLE into the widely used IQ-TREE 2 software, enabling its rapid adoption in the scientific community. These advancements serve as a vital step towards better preparedness for future pandemics, offering researchers powerful tools for large-scale pathogen genomic analysis.

Testing ultraconserved elements (UCEs) for phylogenetic inference across bivalves (Mollusca: Bivalvia)

Bivalves constitute an important resource for fisheries and as cultural objects. Bivalve phylogenetics has had a long tradition using both morphological and molecular characters, and genomic resources are available for a good number of commercially important species. However, relationships among bivalve families have been unstable and major conflicting results exist between mitogenomics and results based on Sanger-based amplicon sequencing or phylotranscriptomics. Here we design and test an ultraconserved elements probe set for the class Bivalvia with the aim to use hundreds of loci without the need to sequence full genomes or transcriptomes, which are expensive and complex to analyze, and to open bivalve phylogenetics to museum specimens. Our probe set successfully captured 1,513 UCEs for a total of 263,800 bp with an average length of 174.59 ± 3.44 per UCE (ranging from 28 to 842 bp). Phylogenetic testing of this UCE probe set across Bivalvia and within the family Donacidae using different data matrices and methods for phylogenetic inference shows promising results at multiple taxonomic levels. In addition, our probe set was able to capture large numbers of UCEs for museum specimens collected before 1900 and from DNAs properly stored, of which many museums and laboratories are well stocked. Overall, this constitutes a novel and useful resource for bivalve phylogenetics.

Study of the Floristic, Morphological, and Genetic (atpF-atpH, Internal Transcribed Spacer (ITS), matK, psbK-psbI, rbcL, and trnH-psbA) Differences in Crataegus ambigua Populations in Mangistau (Kazakhstan).

This article studies the morphological parameters of vegetative and generative organs of different age groups of Crataegus ambigua from four populations in Western Karatau (Mangistau region, Kazakhstan). In this study, we examined four populations: Sultan Epe, Karakozaiym, Emdikorgan, and Samal, all located in various gorges of Western Karatau. Several phylogenetic inference methods were applied, using six genetic markers to reconstruct the evolutionary relationships between these populations: atpF-atpH, internal transcribed spacer (ITS), matK, psbK-psbI, rbcL, and trnH-psbA. We also used a statistical analysis of plants' vegetative and generative organs for three age groups (virgin, young, and adult generative). According to the age structure, Samal has a high concentration of young generative plants (42.3%) and adult generative plants (30.9%). Morphological analysis showed the significance of the parameters of the generative organs and separated the Samal population into a separate group according to the primary principal component analysis (PCoA) coordinates. The results of the floristic analysis showed that the Samal populations have a high concentration of species diversity. Comparative dendrograms using UPGMA (unweighted pair group method with arithmetic mean) showed that information gleaned from genetic markers and the psbK-psbI region can be used to determine the difference between the fourth Samal population and the other three.

Refining the phylogeny and taxonomy of the apple tribe Maleae (Rosaceae): insights from phylogenomic analyses of 563 plastomes and a taxonomic synopsis of Photinia and its allies in the Old World

This study addresses the longstanding absence of a comprehensive phylogenetic backbone for the apple tribe Maleae, a deficiency attributed to limited taxon and marker sampling. We conducted an extensive taxon sampling, incorporating 563 plastomes from a diverse range of 370 species encompassing 26 presently recognized genera. Employing a range of phylogenetic inference methods, including RAxML and IQ-TREE2 for Maximum Likelihood (ML) analyses, we established a robust phylogenetic framework for the Maleae tribe. Our phylogenomic investigations provided compelling support for three major clades within Maleae. By integrating nuclear phylogenetic data with morphological and chromosomal evidence, we propose an updated infra-tribal taxonomic system, comprising subtribe Malinae Reveal, subtribe Lindleyinae Reveal, and subtribe Vauqueliniinae B.B.Liu (subtr. nov.). Plastid phylogenetic analysis also confirmed the monophyly of most genera, except for Amelanchier, Malus, Sorbus sensu lato, and Stranvaesia. In addition, we present a comprehensive taxonomic synopsis of Photinia and its morphological allies in the Old World, recognizing 27 species and ten varieties within Photinia, three species and two varieties within Stranvaesia, and two species and three varieties within Weniomeles. Furthermore, we also lectotypified 12 names and made two new combinations, Photinia microphylla (J.E.Vidal) B.B.Liu and Weniomeles atropurpurea (P.L.Chiu ex Z.H.Chen & X.F.Jin) B.B.Liu.

Mitochondrial genomic characteristics and phylogenetic analysis of a brewing fungus, Rhizopus microsporus Tiegh. 1875 (Mucorales: Rhizopodaceae)

Rhizopus microsporus Tiegh. 1875 is widely used in a variety of industries, such as brewing, wine making, baking, and medicine production, as it has the capability to break down proteins and generate surface-active agents. To date, the mitochondrial genome features of early evolved fungi from the Rhizopus genus have not been extensively studied. Our research obtained a full mitochondrial genome of R. microsporus species, which was 43,837 bp in size and had a GC content of 24.93%. This genome contained 14 core protein-coding genes, 3 independent ORFs, 7 intronic ORFs, 24 tRNAs, and 2 rRNA genes. Through the use of the BI phylogenetic inference method, we were able to create phylogenetic trees for 25 early differentiation fungi which strongly supported the major clades; this indicated that R. microsporus is most closely related to Rhizopus oryzae.

Asymmetric Cluster-Based Measures for Comparative Phylogenetics.

Phylogenetic inference and reconstruction methods generate hypotheses on evolutionary history. Competing inference methods are frequently used, and the evaluation of the generated hypotheses is achieved using tree comparison costs. The Robinson-Foulds (RF) distance is a widely used cost to compare the topology of two trees, but this cost is sensitive to tree error and can overestimate tree differences. To overcome this limitation, a refined version of the RF distance called the Cluster Affinity (CA) distance was introduced. However, CA distances are symmetric and cannot compare different types of trees. These asymmetric comparisons occur when gene trees are compared with species trees, when disparate datasets are integrated into a supertree, or when tree comparison measures are used to infer a phylogenetic network. In this study, we introduce a relaxation of the original Affinity distance to compare heterogeneous trees called the asymmetric CA cost. We also develop a biologically interpretable cost, the Cluster Support cost that normalizes by cluster size across gene trees. The characteristics of these costs are similar to the symmetric CA cost. We describe efficient algorithms, derive the exact diameters, and use these to standardize the cost to be applicable in practice. These costs provide objective, fine-scale, and biologically interpretable values that can assess differences and similarities between phylogenetic trees.

Exploring the Distribution of Phylogenetic Networks Generated Under a Birth-Death-Hybridization Process

Gene-flow processes such as hybridization and introgression play important roles in shaping diversity across the tree of life. Recent studies extending birth-death models have made it possible to investigate patterns of reticulation in a macroevolutionary context. These models allow for different macroevolutionary patterns of gene flow events that can either add, maintain, or remove lineages—with the gene flow itself possibly being dependent on the relatedness between species—thus creating complex diversification scenarios. Further, many reticulate phylogenetic inference methods assume specific reticulation structures or phylogenies belonging to certain network classes. However, the distributions of phylogenetic networks under reticulate birth-death processes are poorly characterized, and it is unknown whether they violate common methodological assumptions. We use simulation techniques to explore phylogenetic network space under a birth-death-hybridization process where the hybridization rate can have a linear dependence on genetic distance. Specifically, we measured the number of lineages through time and role of hybridization in diversification along with the proportion of phylogenetic networks that belong to commonly used network classes (e.g., tree-child, tree-based, or level-1 networks). We find that the growth of phylogenetic networks and class membership are largely affected by assumptions about macroevolutionary patterns of gene flow. In accordance with previous studies, a lower proportion of networks belonged to these classes based on type and density of reticulate events. However, under a birth-death-hybridization process, these factors form an antagonistic relationship; the type of reticulation events that cause high membership proportions also lead to the highest reticulation density, consequently lowering the overall proportion of phylogenies in some classes. Further, we observed that genetic distance–dependent gene flow and incomplete sampling increase the proportion of class membership, primarily due to having fewer reticulate events. Our results can inform studies if their biological expectations of gene flow are associated with evolutionary histories that satisfy the assumptions of current methodology and aid in finding phylogenetic classes that are relevant for methods development.

Phylogenomics of Neogastropoda: The Backbone Hidden in the Bush.

The molluskan order Neogastropoda encompasses over 15,000 almost exclusively marine species playing important roles in benthic communities and in the economies of coastal countries. Neogastropoda underwent intensive cladogenesis in the early stages of diversification, generating a "bush" at the base of their evolutionary tree, which has been hard to resolve even with high throughput molecular data. In the present study to resolve the bush, we use a variety of phylogenetic inference methods and a comprehensive exon capture dataset of 1817 loci (79.6% data occupancy) comprising 112 taxa of 48 out of 60 Neogastropoda families. Our results show consistent topologies and high support in all analyses at (super)family level, supporting monophyly of Muricoidea, Mitroidea, Conoidea, and, with some reservations, Olivoidea and Buccinoidea. Volutoidea and Turbinelloidea as currently circumscribed are clearly paraphyletic. Despite our analyses consistently resolving most backbone nodes, 3 prove problematic: First, the uncertain placement of Cancellariidae, as the sister group to either a Ficoidea-Tonnoidea clade or to the rest of Neogastropoda, leaves monophyly of Neogastropoda unresolved. Second, relationships are contradictory at the base of the major "core Neogastropoda" grouping. Third, coalescence-based analyses reject monophyly of the Buccinoidea in relation to Vasidae. We analyzed phylogenetic signal of targeted loci in relation to potential biases, and we propose the most probable resolutions in the latter 2 recalcitrant nodes. The uncertain placement of Cancellariidae may be explained by orthology violations due to differential paralog loss shortly after the whole genome duplication, which should be resolved with a curated set of longer loci.

Morphological and phylogenetic resolution of Conoideocrella luteorostrata (Hypocreales: Clavicipitaceae), a potential biocontrol fungus for Fiorinia externa in United States Christmas tree production areas

ABSTRACT The entomopathogenic fungus Conoideocrella luteorostrata has recently been implicated in natural epizootics among exotic elongate hemlock scale (EHS) insects in Christmas tree farms in the eastern United States. Since 1913, C. luteorostrata has been reported from various plant feeding Hemiptera in the southeastern United States, but comprehensive morphological and phylogenetic studies of U.S. populations are lacking. The recovery of multiple strains of C. luteorostrata from mycosed EHS in North Carolina provided an opportunity to conduct pathogenicity assays and morphological and phylogenetic studies to investigate genus- and species-level boundaries among the Clavicipitaceae. Pathogenicity assays confirmed that C. luteorostrata causes mortality of EHS crawlers, an essential first step in developing this fungus as a biocontrol. Morphological studies revealed that conidia aligned with previous measurements of the Paecilomyces-like asexual state of C. luteorostrata, with conidiophore morphology consistent with historical observations. Additionally, a Hirsutella-like synanamorph was observed in select C. luteorostrata strains. In both a four-locus, 54-taxon Clavicipitaceae-wide phylogenetic analysis including D1–D2 domains of the nuclear 28S rRNA region (28S), elongation factor 1 alpha (EF1-α), DNA-directed RNA polymerase II subunit 1 (RPB1), and DNA-directed RNA polymerase II subunit 2 (RPB2) and a two-locus, 38-taxon (28S and EF1-α) phylogenetic analysis, all three Conoideocrella species were resolved as strongly supported monophyletic lineages across all loci and both methods (maximum likelihood and Bayesian inference) of phylogenetic inference except for 28S for C. tenuis. Despite the strong support for individual Conoideocrella species, none of the analyses supported the monophyly of Conoideocrella with the inclusion of Dussiella. Due to the paucity of RPB1 and RPB2 sequence data, EF1-α provided superior delimitation of intraspecies groupings for Conoideocrella and should be used in future studies. Further development of C. luteorostrata as a biocontrol against EHS will require additional surveys across diverse Hemiptera and expanded pathogenicity testing to clarify host range and efficacy of this fungus.

Is the Existence of Two Lineages for Hamadryas glauconome (Lepidoptera: Nymphalidae) True? Molecular and Ecological Evidence

The genus Hamadryas has a neotropical distribution. In 1983, the subspecies H. glauconome grisea from Mexico was recognized with subtle and subjective differences in color, size and distribution and limited to the northwest. Since then, there has been a debate about whether it is a different lineage from H. glauconome because adult-stage morphology studies have not found significant differences. This study aims to delimitate H. g. glauconome and H. g. grisea lineages with two sources of evidence: ecological and molecular—the former through ecological niche modeling using the accessible area for the species and estimating the minimum volume ellipsoid overlapping as a fundamental niche using occurrences databases. The molecular evidence is found through the methods of phylogenetic inference and the generalized mixed yule coalescent approach, using sequences of cytochrome oxidase I. Ecological and molecular evidence suggest that H. g. grisea is a different lineage from H. glauconome. Also, molecular evidence of a third lineage from the south of Texas needs further study. This study suggests that different evidence should be provided when morphology is not enough for delimiting species, especially in recently diverged species. Furthermore, the H. g. grisea cytochrome oxidase I sequence (658 bp) is published for the first time.

Reference Genome Choice and Filtering Thresholds Jointly Influence Phylogenomic Analyses.

Molecular phylogenies are a cornerstone of modern comparative biology and are commonly employed to investigate a range of biological phenomena, such as diversification rates, patterns in trait evolution, biogeography, and community assembly. Recent work has demonstrated that significant biases may be introduced into downstream phylogenetic analyses from processing genomic data; however, it remains unclear whether there are interactions among bioinformatic parameters or biases introduced through the choice of reference genome for sequence alignment and variant calling. We address these knowledge gaps by employing a combination of simulated and empirical data sets to investigate the extent to which the choice of reference genome in upstream bioinformatic processing of genomic data influences phylogenetic inference, as well as the way that reference genome choice interacts with bioinformatic filtering choices and phylogenetic inference method. We demonstrate that more stringent minor allele filters bias inferred trees away from the true species tree topology, and that these biased trees tend to be more imbalanced and have a higher center of gravity than the true trees. We find the greatest topological accuracy when filtering sites for minor allele count (MAC) >3-4 in our 51-taxa data sets, while tree center of gravity was closest to the true value when filtering for sites with MAC >1-2. In contrast, filtering for missing data increased accuracy in the inferred topologies; however, this effect was small in comparison to the effect of minor allele filters and may be undesirable due to a subsequent mutation spectrum distortion. The bias introduced by these filters differs based on the reference genome used in short read alignment, providing further support that choosing a reference genome for alignment is an important bioinformatic decision with implications for downstream analyses. These results demonstrate that attributes of the study system and dataset (and their interaction) add important nuance for how best to assemble and filter short-read genomic data for phylogenetic inference.

Phylogenomic conflict analyses of the plastid and mitochondrial genomes via deep genome skimming highlight their independent evolutionary histories: A case study in the cinquefoil genus Potentilla sensu lato (Potentilleae, Rosaceae)

Phylogenomic conflicts are widespread among genomic data, with most previous studies primarily focusing on nuclear datasets instead of organellar genomes. In this study, we investigate phylogenetic conflict analyses within and between plastid and mitochondrial genomes using Potentilla as a case study. We generated three plastid datasets (coding, noncoding, and all-region) and one mitochondrial dataset (coding regions) to infer phylogenies based on concatenated and multispecies coalescent (MSC) methods. Conflict analyses were then performed using PhyParts and Quartet Sampling (QS). Both plastid and mitochondrial genomes divided the Potentilla into eight highly supported clades, two of which were newly identified in this study. While most organellar loci were uninformative for the majority of nodes (bootstrap value < 70%), PhyParts and QS detected conflicting signals within the two organellar genomes. Regression analyses revealed that conflict signals mainly occurred among shorter loci, whereas longer loci tended to be more concordant with the species tree. In addition, two significant disagreements between the two organellar genomes were detected, likely attributed to hybridization and/or incomplete lineage sorting. Our results demonstrate that mitochondrial genes can fully resolve the phylogenetic relationships among eight major clades of Potentilla and are not always linked with plastome in evolutionary history. Stochastic inferences appear to be the primary source of observed conflicts among the gene trees. We recommend that the loci with short sequence length or containing limited informative sites should be used cautiously in MSC analysis, and suggest the joint application of concatenated and MSC methods for phylogenetic inference using organellar genomes.

AliSim-HPC: parallel sequence simulator for phylogenetics.

Sequence simulation plays a vital role in phylogenetics with many applications, such as evaluating phylogenetic methods, testing hypotheses, and generating training data for machine-learning applications. We recently introduced a new simulator for multiple sequence alignments called AliSim, which outperformed existing tools. However, with the increasing demands of simulating large data sets, AliSim is still slow due to its sequential implementation; for example, to simulate millions of sequence alignments, AliSim took several days or weeks. Parallelization has been used for many phylogenetic inference methods but not yet for sequence simulation. This paper introduces AliSim-HPC, which, for the first time, employs high-performance computing for phylogenetic simulations. AliSim-HPC parallelizes the simulation process at both multi-core and multi-CPU levels using the OpenMP and message passing interface (MPI) libraries, respectively. AliSim-HPC is highly efficient and scalable, which reduces the runtime to simulate 100 large gap-free alignments (30 000 sequences of one million sites) from over one day to 11 min using 256 CPU cores from a cluster with six computing nodes, a 153-fold speedup. While the OpenMP version can only simulate gap-free alignments, the MPI version supports insertion-deletion models like the sequential AliSim. AliSim-HPC is open-source and available as part of the new IQ-TREE version v2.2.3 at https://github.com/iqtree/iqtree2/releases with a user manual at http://www.iqtree.org/doc/AliSim.

DecentTree: Scalable Neighbour-Joining for the Genomic Era.

Neighbour-Joining is one of the most widely used distance-based phylogenetic inference methods. However, current implementations do not scale well for datasets with more than 10,000 sequences. Given the increasing pace of generating new sequence data, particularly in outbreaks of emerging diseases, and the already enormous existing databases of sequence data for which NJ is a useful approach, new implementations of existing methods are warranted. Here we present DecentTree, which provides highly optimised and parallel implementations of Neighbour-Joining and several of its variants. DecentTree is designed as a stand-alone application and a header-only library easily integrated with other phylogenetic software (e.g., it is integral in the popular IQ-TREE software). We show that DecentTree shows similar or improved performance over existing software (BIONJ, Quicktree, FastME, and RapidNJ), especially for handling very large alignments. For example, DecentTree is up to 6-fold faster than the fastest existing Neighbour-Joining software (e.g., RapidNJ) when generating a tree of 64,000 SARS-CoV-2 genomes. DecentTree is open source and freely available at https://github.com/iqtree/decenttree. All code and data used in this analysis are available on Github (https://github.com/asdcid/Comparison-of-neighbour-joining-software). Supplementary data are available at Bioinformatics online.

Molecular Characterization and Phylogenetic Analysis of Dengue Fever Viruses in Three Outbreaks in Tanzania Between 2017 and 2019.

Dengue is a disease of public health interest, and Tanzania experienced major outbreaks in 2014 and 2019. Here, we report our findings on the molecular characterization of dengue viruses (DENV) that circulated during two smaller outbreaks (2017 and 2018) and one major epidemic (2019) in Tanzania. We tested archived serum samples from 1,381 suspected dengue fever patients, with a median age of 29 (IQR:22-40) years, referred to the National Public Health Laboratory for confirmation of DENV infection. DENV serotypes were identified by reverse transcription polymerase chain reaction (RT-PCR), and specific genotypes were identified by sequencing the envelope glycoprotein gene and phylogenetic inference methods. DENV was confirmed in 823 (59.6%) cases. More than half (54.7%) of patients with dengue fever infection were males, and nearly three-quarters (73%) of the infected individuals were living in Kinondoni district, Dar es Salaam. DENV-3 Genotype III caused the two smaller outbreaks in 2017 and 2018, while DENV-1 Genotype V caused the 2019 epidemic. DENV-1 Genotype I was also detected in one patient in 2019. This study has demonstrated the molecular diversity of dengue viruses circulating in Tanzania. We found that contemporary circulating serotypes did not cause the major epidemic of 2019 but rather due to a serotype shift from DENV-3 (2017/2018) to DENV-1 in 2019. Such a change increases the risk for patients previously infected with a particular serotype to develop severe symptoms upon potential re-infection with a heterologous serotype due to antibody-dependent enhancement of infection. Therefore, the circulation of serotypes emphasizes the need to strengthen the country's dengue surveillance system for better management of patients, early detection of outbreaks, and vaccine development.

Parsimony and the rank of a flattening matrix

The standard models of sequence evolution on a tree determine probabilities for every character or site pattern. A flattening is an arrangement of these probabilities into a matrix, with rows corresponding to all possible site patterns for one set A of taxa and columns corresponding to all site patterns for another set B of taxa. Flattenings have been used to prove difficult results relating to phylogenetic invariants and consistency and also form the basis of several methods of phylogenetic inference. We prove that the rank of the flattening equals r^{nu _T(A|B)}, where r is the number of states and nu _T(A|B) is the minimum size of a vertex cut separating A from B. When T is binary the rank of the flattening equals r^{ell _T(A|B)} where ell _T(A|B) equals the parsimony length of the binary character separating A and B. We provide a direct proof that requires little more than undergraduate algebra, but note that the formula could also be derived from work by Casanellas and Fernández-Sánchez (2011) on phylogenetic invariants.

Revisiting chloroplast genomic landscape and annotation towards comparative chloroplast genomes of Rhamnaceae

BackgroundMassive parallel sequencing technologies have enabled the elucidation of plant phylogenetic relationships from chloroplast genomes at a high pace. These include members of the family Rhamnaceae. The current Rhamnaceae phylogenetic tree is from 13 out of 24 Rhamnaceae chloroplast genomes, and only one chloroplast genome of the genus Ventilago is available. Hence, the phylogenetic relationships in Rhamnaceae remain incomplete, and more representative species are needed.ResultsThe complete chloroplast genome of Ventilago harmandiana Pierre was outlined using a hybrid assembly of long- and short-read technologies. The accuracy and validity of the final genome were confirmed with PCR amplifications and investigation of coverage depth. Sanger sequencing was used to correct for differences in lengths and nucleotide bases between inverted repeats because of the homopolymers. The phylogenetic trees reconstructed using prevalent methods for phylogenetic inference were topologically similar. The clustering based on codon usage was congruent with the molecular phylogenetic tree. The groups of genera in each tribe were in accordance with tribal classification based on molecular markers. We resolved the phylogenetic relationships among six Hovenia species, three Rhamnus species, and two Ventilago species. Our reconstructed tree provides the most complete and reliable low-level taxonomy to date for the family Rhamnaceae. Similar to other higher plants, the RNA editing mostly resulted in converting serine to leucine. Besides, most genes were subjected to purifying selection. Annotation anomalies, including indel calling errors, unaligned open reading frames of the same gene, inconsistent prediction of intergenic regions, and misannotated genes, were identified in the published chloroplast genomes used in this study. These could be a result of the usual imperfections in computational tools, and/or existing errors in reference genomes. Importantly, these are points of concern with regards to utilizing published chloroplast genomes for comparative genomic analysis.ConclusionsIn summary, we successfully demonstrated the use of comprehensive genomic data, including DNA and amino acid sequences, to build a reliable and high-resolution phylogenetic tree for the family Rhamnaceae. Additionally, our study indicates that the revision of genome annotation before comparative genomic analyses is necessary to prevent the propagation of errors and complications in downstream analysis and interpretation.

Systematics of Miocene apes: State of the art of a neverending controversy

Hominoids diverged from cercopithecoids during the Oligocene in Afro-Arabia, initially radiating in that continent and subsequently dispersing into Eurasia. From the Late Miocene onward, the geographic range of hominoids progressively shrank, except for hominins, which dispersed out of Africa during the Pleistocene. Although the overall picture of hominoid evolution is clear based on available fossil evidence, many uncertainties persist regarding the phylogeny and paleobiogeography of Miocene apes (nonhominin hominoids), owing to their sparse record, pervasive homoplasy, and the decimated current diversity of this group. We review Miocene ape systematics and evolution by focusing on the most parsimonious cladograms published during the last decade. First, we provide a historical account of the progress made in Miocene ape phylogeny and paleobiogeography, report an updated classification of Miocene apes, and provide a list of Miocene ape species-locality occurrences together with an analysis of their paleobiodiversity dynamics. Second, we discuss various critical issues of Miocene ape phylogeny and paleobiogeography (hylobatid and crown hominid origins, plus the relationships of Oreopithecus) in the light of the highly divergent results obtained from cladistic analyses of craniodental and postcranial characters separately. We conclude that cladistic efforts to disentangle Miocene ape phylogeny are potentially biased by a long-branch attraction problem caused by the numerous postcranial similarities shared between hylobatids and hominids—despite the increasingly held view that they are likely homoplastic to a large extent, as illustrated by Sivapithecus and Pierolapithecus—and further aggravated by abundant missing data owing to incomplete preservation. Finally, we argue that—besides the recovery of additional fossils, the retrieval of paleoproteomic data, and a better integration between cladistics and geometric morphometrics—Miocene ape phylogenetics should take advantage of total-evidence (tip-dating) Bayesian methods of phylogenetic inference combining morphologic, molecular, and chronostratigraphic data. This would hopefully help ascertain whether hylobatid divergence was more basal than currently supported.

The Free Lunch is not over yet-systematic exploration of numerical thresholds in maximum likelihood phylogenetic inference.

Maximum likelihood (ML) is a widely used phylogenetic inference method. ML implementations heavily rely on numerical optimization routines that use internal numerical thresholds to determine convergence. We systematically analyze the impact of these threshold settings on the log-likelihood and runtimes for ML tree inferences with RAxML-NG, IQ-TREE, and FastTree on empirical datasets. We provide empirical evidence that we can substantially accelerate tree inferences with RAxML-NG and IQ-TREE by changing the default values of two such numerical thresholds. At the same time, altering these settings does not significantly impact the quality of the inferred trees. We further show that increasing both thresholds accelerates the RAxML-NG bootstrap without influencing the resulting support values. For RAxML-NG, increasing the likelihood thresholds and to 10 and 103, respectively, results in an average tree inference speedup of 1.9 ± 0.6 on Data collection 1, 1.8 ± 1.1 on Data collection 2, and 1.9 ± 0.8 on Data collection 2 for the RAxML-NG bootstrap compared to the runtime under the current default setting. Increasing the likelihood threshold to 10 in IQ-TREE results in an average tree inference speedup of 1.3 ± 0.4 on Data collection 1 and 1.3 ± 0.9 on Data collection 2. All MSAs we used for our analyses, as well as all results, are available for download at https://cme.h-its.org/exelixis/material/freeLunch_data.tar.gz. Our data generation scripts are available at https://github.com/tschuelia/ml-numerical-analysis.