The MRP Method

Multiomics analysis reveals B. MO1 as a distinct Babesia species and provides insights into its evolution and virulence

Babesiosis, caused by protozoan parasites of the genus Babesia , is an emerging tick-borne disease of significance for both human and animal health. Babesia parasites infect erythrocytes of vertebrate hosts where they develop and multiply rapidly to cause the pathological symptoms associated with the disease. The identification of various Babesia species underscores the ongoing risk of new zoonotic pathogens capable of infecting humans, a concern amplified by anthropogenic activities and environmental shifts impacting the distribution and transmission dynamics of parasites, their vectors, and reservoir hosts. One such species, Babesia MO1, previously implicated in severe cases of human babesiosis in the midwestern United States, was initially considered closely related to B. divergens , the predominant agent of human babesiosis in Europe. Yet, uncertainties persist regarding whether these pathogens represent distinct variants of the same species or are entirely separate species. We show that although both B. MO1 and B. divergens share similar genome sizes, comprising three nuclear chromosomes, one linear mitochondrial chromosome, and one circular apicoplast chromosome, major differences exist in terms of genomic sequence divergence, gene functions, transcription profiles, replication rates and susceptibility to antiparasitic drugs. Furthermore, both pathogens have evolved distinct classes of multigene families, crucial for their pathogenicity and adaptation to specific mammalian hosts. Leveraging genomic information for B. MO1, B. divergens , and other members of the Babesiidae family within Apicomplexa provides valuable insights into the evolution, diversity, and virulence of these parasites. This knowledge serves as a critical tool in preemptively addressing the emergence and rapid transmission of more virulent strains.

Decoupling speciation and extinction reveals both abiotic and biotic drivers shaped 250 million years of diversity in crocodile-line archosaurs

Whereas living representatives of Pseudosuchia, crocodylians, number fewer than 30 species, more than 700 pseudosuchian species are known from their 250-million-year fossil record, displaying far greater ecomorphological diversity than their extant counterparts. With a new time-calibrated tree of >500 species, we use a phylogenetic framework to reveal that pseudosuchian evolutionary history and diversification dynamics were directly shaped by the interplay of abiotic and biotic processes over hundreds of millions of years, supported by information theory analyses. Speciation, but not extinction, is correlated with higher temperatures in terrestrial and marine lineages, with high sea level associated with heightened extinction in non-marine taxa. Low lineage diversity and increased speciation in non-marine species is consistent with opportunities for niche-filling, whereas increased competition may have led to elevated extinction rates. In marine lineages, competition via increased lineage diversity appears to have driven both speciation and extinction. Decoupling speciation and extinction, in combination with ecological partitioning, reveals a more complex picture of pseudosuchian evolution than previously understood. As the number of species threatened with extinction by anthropogenic climate change continues to rise, the fossil record provides a unique window into the drivers that led to clade success and those that may ultimately lead to extinction.

Ecological Transitions and the Shape of the Decapod Tree of Life

Understanding the processes that shaped the distribution of species richness across the Tree of Life is a central macroevolutionary research agenda. Major ecological innovations, including transitions between habitats, may help to explain the striking asymmetries of diversity that are often observed between sister clades. Here, we test the impact of such transitions on speciation rates across decapod crustaceans, modelling diversification dynamics within a phylogenetic framework. Our results show that, while terrestrial lineages have higher speciation rates than either marine or freshwater lineages, there is no difference between mean speciation rates in marine and freshwater lineages across Decapoda. Partitioning our data by infraorder reveals that those clades with habitat heterogeneity have higher speciation rates in freshwater and terrestrial lineages, with freshwater rates up to 1.5 times faster than marine rates, and terrestrial rates approximately four times faster. This averaging out of marine and freshwater speciation rates results from the varying contributions of different clades to average speciation rates. However, with the exception of Caridea, we find no evidence for any causal relationship between habitat and speciation rate. Our results demonstrate that while statistical generalisations about ecological traits and evolutionary rates are valuable, there are many exceptions. Hence, while freshwater and terrestrial lineages typically speciate faster than their marine relatives, there are many atypically slow freshwater lineages and fast marine lineages across Decapoda. Future work on diversification patterns will benefit from the inclusion of fossil data, as well as additional ecological factors.

Comparing Phylogenetic Trees Side by Side Through iPhyloC, a New Interactive Web-Based Framework

Current frameworks of side-by-side phylogenetic trees comparison face two issues: (1) accepting mainly binary trees as input and (2) assuming input trees having identical or highly overlapping taxa. However, cladistic comparative studies often lead with multiple nontotally resolved trees with nonidentical sets of taxa. We tackle these issues in this study, presenting the iPhyloC, an interactive web-based framework for comparing phylogenetic trees side by side. iPhyloC supports automatic identification of the common taxa in the input trees, comparison options between them, intuitive design, high usability, scalability to large trees, and cross-platform support. iPhyloC was tested using different trees and a supertree depicting the phylogenetic relationships within the insect order Diptera as examples.

BuM 2.0: Software for online generation of matrices for Brooks Parsimony Analysis

Brooks Parsimony Analysis is a cladistic biogeographic method that aims to extract biogeographic information from phylogenetic trees, depicting from a group of cladograms a general pattern of relationships among the areas the taxa inhabit. We present here BuM 2.0, an online framework to automatically create matrices for Brooks Parsimony Analysis (BPA) using Baum Ragans algorithm for Matrix Representation with Parsimony.

Freshwater transitions and symbioses shaped the evolution and extant diversity of caridean shrimps

Understanding the processes that shaped the strikingly irregular distribution of species richness across the Tree of Life is a major research agenda. Changes in ecology may go some way to explain the often strongly asymmetrical fates of sister clades, and we test this in the caridean shrimps. First appearing in the Lower Jurassic, there are now ~3500 species worldwide. Carideans experienced several independent transitions to freshwater from marine habitats, while many of the marine species have also evolved a symbiotic lifestyle. Here we use diversification rate analyses to test whether these ecological traits promote or inhibit diversity within a phylogenetic framework. We demonstrate that speciation rates are more than twice as high in freshwater clades, whilst symbiotic ecologies are associated with lower speciation rates. These lower rates amongst symbiotic species are of concern given that symbioses often occur in some of the most diverse, delicately balanced and threatened marine ecosystems.

Katie Davis et al. test the hypothesis that ecological traits are linked to diversification in caridean shrimps. They find that transitions from marine to freshwater habitats contributed to higher diversification rates, whereas symbiosis is associated with a slight decrease in diversification rates.

Improved Heuristics for Minimum-Flip Supertree Construction

The utility of the matrix representation with flipping (MRF) supertree method has been limited by the speed of its heuristic algorithms. We describe a new heuristic algorithm for MRF supertree construction that improves upon the speed of the previous heuristic by a factor of n (the number of taxa in the supertree). This new heuristic makes MRF tractable for large-scale supertree analyses and allows the first comparisons of MRF with other supertree methods using large empirical data sets. Analyses of three published supertree data sets with between 267 to 571 taxa indicate that MRF supertrees are equally or more similar to the input trees on average than matrix representation with parsimony (MRP) and modified mincut supertrees. The results also show that large differences may exist between MRF and MRP supertrees and demonstrate that the MRF supertree method is a practical and potentially more accurate alternative to the nearly ubiquitous MRP super-tree method.

Fin modules: an evolutionary perspective on appendage disparity in basal vertebrates

BACKGROUND

Fishes are extremely speciose and also highly disparate in their fin configurations, more specifically in the number of fins present as well as their structure, shape, and size. How they achieved this remarkable disparity is difficult to explain in the absence of any comprehensive overview of the evolutionary history of fish appendages. Fin modularity could provide an explanation for both the observed disparity in fin configurations and the sequential appearance of new fins. Modularity is considered as an important prerequisite for the evolvability of living systems, enabling individual modules to be optimized without interfering with others. Similarities in developmental patterns between some of the fins already suggest that they form developmental modules during ontogeny. At a macroevolutionary scale, these developmental modules could act as evolutionary units of change and contribute to the disparity in fin configurations. This study addresses fin disparity in a phylogenetic perspective, while focusing on the presence/absence and number of each of the median and paired fins.

RESULTS

Patterns of fin morphological disparity were assessed by mapping fin characters on a new phylogenetic supertree of fish orders. Among agnathans, disparity in fin configurations results from the sequential appearance of novel fins forming various combinations. Both median and paired fins would have appeared first as elongated ribbon-like structures, which were the precursors for more constricted appendages. Among chondrichthyans, disparity in fin configurations relates mostly to median fin losses. Among actinopterygians, fin disparity involves fin losses, the addition of novel fins (e.g., the adipose fin), and coordinated duplications of the dorsal and anal fins. Furthermore, some pairs of fins, notably the dorsal/anal and pectoral/pelvic fins, show non-independence in their character distribution, supporting expectations based on developmental and morphological evidence that these fin pairs form evolutionary modules.

CONCLUSIONS

Our results suggest that the pectoral/pelvic fins and the dorsal/anal fins form two distinct evolutionary modules, and that the latter is nested within a more inclusive median fins module. Because the modularity hypotheses that we are testing are also supported by developmental and variational data, this constitutes a striking example linking developmental, variational, and evolutionary modules.

Global cooling as a driver of diversification in a major marine clade

Climate is a strong driver of global diversity and will become increasingly important as human influences drive temperature changes at unprecedented rates. Here we investigate diversification and speciation trends within a diverse group of aquatic crustaceans, the Anomura. We use a phylogenetic framework to demonstrate that speciation rate is correlated with global cooling across the entire tree, in contrast to previous studies. Additionally, we find that marine clades continue to show evidence of increased speciation rates with cooler global temperatures, while the single freshwater clade shows the opposite trend with speciation rates positively correlated to global warming. Our findings suggest that both global cooling and warming lead to diversification and that habitat plays a role in the responses of species to climate change. These results have important implications for our understanding of how extant biota respond to ongoing climate change and are of particular importance for conservation planning of marine ecosystems.

In many groups of organisms, speciation rates are higher when global temperatures are warmer. Here, Davis et al. find that marine crustaceans in the Anomura clade have higher speciation rates during cooler periods, whereas their freshwater relatives have the more typical relationship of higher speciation rates in warmer periods.

Phylogenomic species tree estimation in the presence of incomplete lineage sorting and horizontal gene transfer

BACKGROUND

Species tree estimation is challenged by gene tree heterogeneity resulting from biological processes such as duplication and loss, hybridization, incomplete lineage sorting (ILS), and horizontal gene transfer (HGT). Mathematical theory about reconstructing species trees in the presence of HGT alone or ILS alone suggests that quartet-based species tree methods (known to be statistically consistent under ILS, or under bounded amounts of HGT) might be effective techniques for estimating species trees when both HGT and ILS are present.

RESULTS

We evaluated several publicly available coalescent-based methods and concatenation under maximum likelihood on simulated datasets with moderate ILS and varying levels of HGT. Our study shows that two quartet-based species tree estimation methods (ASTRAL-2 and weighted Quartets MaxCut) are both highly accurate, even on datasets with high rates of HGT. In contrast, although NJst and concatenation using maximum likelihood are highly accurate under low HGT, they are less robust to high HGT rates.

CONCLUSION

Our study shows that quartet-based species-tree estimation methods can be highly accurate under the presence of both HGT and ILS. The study suggests the possibility that some quartet-based methods might be statistically consistent under phylogenomic models of gene tree heterogeneity with both HGT and ILS.

Consensus and confusion in molluscan trees: evaluating morphological and molecular phylogenies

Mollusks are the most morphologically disparate living animal phylum, they have diversified into all habitats, and have a deep fossil record. Monophyly and identity of their eight living classes is undisputed, but relationships between these groups and patterns of their early radiation have remained elusive. Arguments about traditional morphological phylogeny focus on a small number of topological concepts but often without regard to proximity of the individual classes. In contrast, molecular studies have proposed a number of radically different, inherently contradictory, and controversial sister relationships. Here, we assembled a data set of 42 unique published trees describing molluscan interrelationships. We used these data to ask several questions about the state of resolution of molluscan phylogeny compared with a null model of the variation possible in random trees constructed from a monophyletic assemblage of eight terminals. Although 27 different unique trees have been proposed from morphological inference, the majority of these are not statistically different from each other. Within the available molecular topologies, only four studies to date have included the deep sea class Monoplacophora; but 36.4% of all trees are not significantly different. We also present supertrees derived from two data partitions and three methods, including all available molecular molluscan phylogenies, which will form the basis for future hypothesis testing. The supertrees presented here were not constructed to provide yet another hypothesis of molluscan relationships, but rather to algorithmically evaluate the relationships present in the disparate published topologies. Based on the totality of available evidence, certain patterns of relatedness among constituent taxa become clear. The internodal distance is consistently short between a few taxon pairs, particularly supporting the relatedness of Monoplacophora and the chitons, Polyplacophora. Other taxon pairs are rarely or never found in close proximity, such as the vermiform Caudofoveata and Bivalvia. Our results have specific utility for guiding constructive research planning to better test relationships in Mollusca as well as other problematic groups. Taxa with consistently proximate relationships should be the focus of a combined approach in a concerted assessment of potential genetic and anatomical homology, whereas unequivocally distant taxa will make the most constructive choices for exemplar selection in higher level phylogenomic analyses.

Accurate phylogenetic tree reconstruction from quartets: a heuristic approach

Supertree methods construct trees on a set of taxa (species) combining many smaller trees on the overlapping subsets of the entire set of taxa. A ‘quartet’ is an unrooted tree over taxa, hence the quartet-based supertree methods combine many -taxon unrooted trees into a single and coherent tree over the complete set of taxa. Quartet-based phylogeny reconstruction methods have been receiving considerable attentions in the recent years. An accurate and efficient quartet-based method might be competitive with the current best phylogenetic tree reconstruction methods (such as maximum likelihood or Bayesian MCMC analyses), without being as computationally intensive. In this paper, we present a novel and highly accurate quartet-based phylogenetic tree reconstruction method. We performed an extensive experimental study to evaluate the accuracy and scalability of our approach on both simulated and biological datasets.

Sequencing of the smallest Apicomplexan genome from the human pathogen Babesia microti

We have sequenced the genome of the emerging human pathogen Babesia microti and compared it with that of other protozoa. B. microti has the smallest nuclear genome among all Apicomplexan parasites sequenced to date with three chromosomes encoding ∼3500 polypeptides, several of which are species specific. Genome-wide phylogenetic analyses indicate that B. microti is significantly distant from all species of Babesidae and Theileridae and defines a new clade in the phylum Apicomplexa. Furthermore, unlike all other Apicomplexa, its mitochondrial genome is circular. Genome-scale reconstruction of functional networks revealed that B. microti has the minimal metabolic requirement for intraerythrocytic protozoan parasitism. B. microti multigene families differ from those of other protozoa in both the copy number and organization. Two lateral transfer events with significant metabolic implications occurred during the evolution of this parasite. The genomic sequencing of B. microti identified several targets suitable for the development of diagnostic assays and novel therapies for human babesiosis.

MRL and SuperFine+MRL: new supertree methods

Background

Supertree methods combine trees on subsets of the full taxon set together to produce a tree on the entire set of taxa. Of the many supertree methods, the most popular is MRP (Matrix Representation with Parsimony), a method that operates by first encoding the input set of source trees by a large matrix (the "MRP matrix") over {0,1, ?}, and then running maximum parsimony heuristics on the MRP matrix. Experimental studies evaluating MRP in comparison to other supertree methods have established that for large datasets, MRP generally produces trees of equal or greater accuracy than other methods, and can run on larger datasets. A recent development in supertree methods is SuperFine+MRP, a method that combines MRP with a divide-and-conquer approach, and produces more accurate trees in less time than MRP. In this paper we consider a new approach for supertree estimation, called MRL (Matrix Representation with Likelihood). MRL begins with the same MRP matrix, but then analyzes the MRP matrix using heuristics (such as RAxML) for 2-state Maximum Likelihood.

Results

We compared MRP and SuperFine+MRP with MRL and SuperFine+MRL on simulated and biological datasets. We examined the MRP and MRL scores of each method on a wide range of datasets, as well as the resulting topological accuracy of the trees. Our experimental results show that MRL, coupled with a very good ML heuristic such as RAxML, produced more accurate trees than MRP, and MRL scores were more strongly correlated with topological accuracy than MRP scores.

Conclusions

SuperFine+MRP, when based upon a good MP heuristic, such as TNT, produces among the best scores for both MRP and MRL, and is generally faster and more topologically accurate than other supertree methods we tested.

Vertebral evolution and the diversification of squamate reptiles

Taxonomic, morphological, and functional diversity are often discordant and independent components of diversity. A fundamental and largely unanswered question in evolutionary biology is why some clades diversify primarily in some of these components and not others. Dramatic variation in trunk vertebral numbers (14 to >300) among squamate reptiles coincides with different body shapes, and snake-like body shapes have evolved numerous times. However, whether increased evolutionary rates or numbers of vertebrae underlie body shape and taxonomic diversification is unknown. Using a supertree of squamates including 1375 species, and corresponding vertebral and body shape data, we show that increased rates of evolution in vertebral numbers have coincided with increased rates and disparity in body shape evolution, but not changes in rates of taxonomic diversification. We also show that the evolution of many vertebrae has not spurred or inhibited body shape or taxonomic diversification, suggesting that increased vertebral number is not a key innovation. Our findings demonstrate that lineage attributes such as the relaxation of constraints on vertebral number can facilitate the evolution of novel body shapes, but that different factors are responsible for body shape and taxonomic diversification.

Accuracy of phylogeny reconstruction methods combining overlapping gene data sets

Background

The availability of many gene alignments with overlapping taxon sets raises the question of which strategy is the best to infer species phylogenies from multiple gene information. Methods and programs abound that use the gene alignment in different ways to reconstruct the species tree. In particular, different methods combine the original data at different points along the way from the underlying sequences to the final tree. Accordingly, they are classified into superalignment, supertree and medium-level approaches. Here, we present a simulation study to compare different methods from each of these three approaches.

Results

We observe that superalignment methods usually outperform the other approaches over a wide range of parameters including sparse data and gene-specific evolutionary parameters. In the presence of high incongruency among gene trees, however, other combination methods show better performance than the superalignment approach. Surprisingly, some supertree and medium-level methods exhibit, on average, worse results than a single gene phylogeny with complete taxon information.

Conclusions

For some methods, using the reconstructed gene tree as an estimation of the species tree is superior to the combination of incomplete information. Superalignment usually performs best since it is less susceptible to stochastic error. Supertree methods can outperform superalignment in the presence of gene-tree conflict.

Phylogenetic analysis of sexual systems in Inuleae (Asteraceae)

From an ancestor with bisexual flowers, plants with unisexual flowers, or even unisexual individuals have evolved in different lineages of angiosperms. The Asteraceae tribe Inuleae includes hermaphroditic, monoecious, dioecious, and gynomonoecious species. Gynomonoecy, the sexual system in which female and bisexual flowers occur on the same plant, is prevalent in the Asteraceae. We inferred one large gene phylogeny (ndhF) and two supertrees to investigate whether gynomonoecy was a stage in the evolution from hermaphroditism to monoecy. We identified transitions in sexual system evolution using the stochastic character mapping method. From gynomonoecious ancestors, both hermaphroditic and monoecious descendants have evolved. Gynomonoecy was not restricted to a stage in the evolution toward monoecy because the number of transitions and the rate of change from monoecy to gynomonoecy were much higher than the opposite. We also investigated one hypothesized association among female flowers and the development of a petaloid ray as an explanation of gynomonoecy maintenance in Asteraceae. We found that peripheral female flowers and petaloid rays were phylogenetically correlated. However, empirical evidence shows that a causal relationship between these traits is not clear.

PhySIC_IST: cleaning source trees to infer more informative supertrees

Background

Supertree methods combine phylogenies with overlapping sets of taxa into a larger one. Topological conflicts frequently arise among source trees for methodological or biological reasons, such as long branch attraction, lateral gene transfers, gene duplication/loss or deep gene coalescence. When topological conflicts occur among source trees, liberal methods infer supertrees containing the most frequent alternative, while veto methods infer supertrees not contradicting any source tree, i.e. discard all conflicting resolutions. When the source trees host a significant number of topological conflicts or have a small taxon overlap, supertree methods of both kinds can propose poorly resolved, hence uninformative, supertrees.

Results

To overcome this problem, we propose to infer non-plenary supertrees, i.e. supertrees that do not necessarily contain all the taxa present in the source trees, discarding those whose position greatly differs among source trees or for which insufficient information is provided. We detail a variant of the PhySIC veto method called PhySIC_IST that can infer non-plenary supertrees. PhySIC_IST aims at inferring supertrees that satisfy the same appealing theoretical properties as with PhySIC, while being as informative as possible under this constraint. The informativeness of a supertree is estimated using a variation of the CIC (Cladistic Information Content) criterion, that takes into account both the presence of multifurcations and the absence of some taxa. Additionally, we propose a statistical preprocessing step called STC (Source Trees Correction) to correct the source trees prior to the supertree inference. STC is a liberal step that removes the parts of each source tree that significantly conflict with other source trees. Combining STC with a veto method allows an explicit trade-off between veto and liberal approaches, tuned by a single parameter.

Performing large-scale simulations, we observe that STC+PhySIC_IST infers much more informative supertrees than PhySIC, while preserving low type I error compared to the well-known MRP method. Two biological case studies on animals confirm that the STC preprocess successfully detects anomalies in the source trees while STC+PhySIC_IST provides well-resolved supertrees agreeing with current knowledge in systematics.

Conclusion

The paper introduces and tests two new methodologies, PhySIC_IST and STC, that demonstrate the interest in inferring non-plenary supertrees as well as preprocessing the source trees. An implementation of the methods is available at: .

Conditioned genome reconstruction: how to avoid choosing the conditioning genome

Genome phylogenies can be inferred from data on the presence and absence of genes across taxa. Logdet distances may be a good method, because they allow expected genome size to vary across the tree. Recently, Lake and Rivera proposed conditioned genome reconstruction (calculation of logdet distances using only those genes present in a conditioning genome) to deal with unobservable genes that are absent from every taxon of interest. We prove that their method can consistently estimate the topology for almost any choice of conditioning genome. Nevertheless, the choice of conditioning genome is important for small samples. For real bacterial genome data, different choices of conditioning genome can result in strong bootstrap support for different tree topologies. To overcome this problem, we developed supertree methods that combine information from all choices of conditioning genome. One of these methods, based on the BIONJ algorithm, performs well on simulated data and may have applications to other supertree problems. However, an analysis of 40 bacterial genomes using this method supports an incorrect clade of parasites. This is a common feature of model-based gene content methods and is due to parallel gene loss.

Supertree bootstrapping methods for assessing phylogenetic variation among genes in genome-scale data sets

Nonparamtric bootstrapping methods may be useful for assessing confidence in a supertree inference. We examined the performance of two supertree bootstrapping methods on four published data sets that each include sequence data from more than 100 genes. In "input tree bootstrapping," input gene trees are sampled with replacement and then combined in replicate supertree analyses; in "stratified bootstrapping," trees from each gene's separate (conventional) bootstrap tree set are sampled randomly with replacement and then combined. Generally, support values from both supertree bootstrap methods were similar or slightly lower than corresponding bootstrap values from a total evidence, or supermatrix, analysis. Yet, supertree bootstrap support also exceeded supermatrix bootstrap support for a number of clades. There was little overall difference in support scores between the input tree and stratified bootstrapping methods. Results from supertree bootstrapping methods, when compared to results from corresponding supermatrix bootstrapping, may provide insights into patterns of variation among genes in genome-scale data sets.