RecPhyloXML: a format for reconciled gene trees

[genus]_[species]; Presenting phylogenies to facilitate synthesis

Each published phylogeny is a potential contribution to the synthesis of the Tree of Life and countless downstream projects. Steps are needed for fully synthesizable science, but only a minority of studies achieve these. We here review the range of phylogenetic presentation and note aspects that hinder further analysis. We provide simple suggestions on publication that would greatly enhance utilizability, and propose a formal grammar for phylogeny terminal format. We suggest that each published phylogeny should be accompanied by at minimum the single preferred result in machine readable tree (e.g. Newick) form in the supplement, a simple task fulfilled by fewer than half of studies. Further, the tree should be clear from the file name and extension; the orientation (rooted or unrooted) should match the figures; terminals labels should include genus and species IDs; underscores should separate strings within-field (instead of white spaces); and if other informational fields are added these should be separated by a unique delimiting character (we suggest multiple underscores or the vertical pipe character, |) and ordered consistently. These requirements are largely independent of phylogenetic study aims, while we note other requirements for synthesis (e.g. removal of species repeats and uninformative terminals) that are not necessarily the responsibility of authors. Machine readable trees show greater variation in terminal formatting than typical phylogeny images (owing presumably to greater scrutiny of the latter), and thus are complex and laborious to parse. Since the majority of existing studies have provided only images, we additionally review typical variation in plotting style, information that will be necessary for developing the automated phylogeny transcription tools needed for their eventual inclusion in the Tree of Life.

AleRax: a tool for gene and species tree co-estimation and reconciliation under a probabilistic model of gene duplication, transfer, and loss

MOTIVATION

Genomes are a rich source of information on the pattern and process of evolution across biological scales. How best to make use of that information is an active area of research in phylogenetics. Ideally, phylogenetic methods should not only model substitutions along gene trees, which explain differences between homologous gene sequences, but also the processes that generate the gene trees themselves along a shared species tree. To conduct accurate inferences, one needs to account for uncertainty at both levels, that is, in gene trees estimated from inherently short sequences and in their diverse evolutionary histories along a shared species tree.

RESULTS

We present AleRax, a software that can infer reconciled gene trees together with a shared species tree using a simple, yet powerful, probabilistic model of gene duplication, transfer, and loss. A key feature of AleRax is its ability to account for uncertainty in the gene tree and its reconciliation by using an efficient approximation to calculate the joint phylogenetic-reconciliation likelihood and sample reconciled gene trees accordingly. Simulations and analyses of empirical data show that AleRax is one order of magnitude faster than competing gene tree inference tools while attaining the same accuracy. It is consistently more robust than species tree inference methods such as SpeciesRax and ASTRAL-Pro 2 under gene tree uncertainty. Finally, AleRax can process multiple gene families in parallel thereby allowing users to compare competing phylogenetic hypotheses and estimate model parameters, such as duplication, transfer, and loss probabilities for genome-scale datasets with hundreds of taxa.

AVAILABILITY AND IMPLEMENTATION

GNU GPL at https://github.com/BenoitMorel/AleRax and data are made available at https://cme.h-its.org/exelixis/material/alerax_data.tar.gz.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a crucial window into early vertebrate evolution1-3. Here we investigate the complex history, timing and functional role of genome-wide duplications4-7 and programmed DNA elimination8,9 in vertebrates in the light of a chromosome-scale genome sequence for the brown hagfish Eptatretus atami. Combining evidence from syntenic and phylogenetic analyses, we establish a comprehensive picture of vertebrate genome evolution, including an auto-tetraploidization (1RV) that predates the early Cambrian cyclostome-gnathostome split, followed by a mid-late Cambrian allo-tetraploidization (2RJV) in gnathostomes and a prolonged Cambrian-Ordovician hexaploidization (2RCY) in cyclostomes. Subsequently, hagfishes underwent extensive genomic changes, with chromosomal fusions accompanied by the loss of genes that are essential for organ systems (for example, genes involved in the development of eyes and in the proliferation of osteoclasts); these changes account, in part, for the simplification of the hagfish body plan1,2. Finally, we characterize programmed DNA elimination in hagfish, identifying protein-coding genes and repetitive elements that are deleted from somatic cell lineages during early development. The elimination of these germline-specific genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline and pluripotency functions, paralleling findings in lampreys10,11. Reconstruction of the early genomic history of vertebrates provides a framework for further investigations of the evolution of cyclostomes and jawed vertebrates.

AGO, a Framework for the Reconstruction of Ancestral Syntenies and Gene Orders

Reconstructing ancestral gene orders from the genome data of extant species is an important problem in comparative and evolutionary genomics. In a phylogenomics setting that accounts for gene family evolution through gene duplication and gene loss, the reconstruction of ancestral gene orders involves several steps, including multiple sequence alignment, the inference of reconciled gene trees, and the inference of ancestral syntenies and gene adjacencies. For each of the steps of such a process, several methods can be used and implemented using a growing corpus of, often parameterized, tools; in practice, interfacing such tools into an ancestral gene order reconstruction pipeline is far from trivial. This chapter introduces AGO, a Python-based framework aimed at creating ancestral gene order reconstruction pipelines allowing to interface and parameterize different bioinformatics tools. The authors illustrate the features of AGO by reconstructing ancestral gene orders for the X chromosome of three ancestral Anopheles species using three different pipelines. AGO is freely available at https://github.com/cchauve/AGO-pipeline .

The evolution and spread of sulfur cycling enzymes reflect the redox state of the early Earth

The biogeochemical sulfur cycle plays a central role in fueling microbial metabolisms, regulating the Earth's redox state, and affecting climate. However, geochemical reconstructions of the ancient sulfur cycle are confounded by ambiguous isotopic signals. We use phylogenetic reconciliation to ascertain the timing of ancient sulfur cycling gene events across the tree of life. Our results suggest that metabolisms using sulfide oxidation emerged in the Archean, but those involving thiosulfate emerged only after the Great Oxidation Event. Our data reveal that observed geochemical signatures resulted not from the expansion of a single type of organism but were instead associated with genomic innovation across the biosphere. Moreover, our results provide the first indication of organic sulfur cycling from the Mid-Proterozoic onwards, with implications for climate regulation and atmospheric biosignatures. Overall, our results provide insights into how the biological sulfur cycle evolved in tandem with the redox state of the early Earth.

The monoaminergic system is a bilaterian innovation

Monoamines like serotonin, dopamine, and adrenaline/noradrenaline (epinephrine/norepinephrine) act as neuromodulators in the nervous system. They play a role in complex behaviours, cognitive functions such as learning and memory formation, as well as fundamental homeostatic processes such as sleep and feeding. However, the evolutionary origin of the genes required for monoaminergic modulation is uncertain. Using a phylogenomic approach, in this study, we show that most of the genes involved in monoamine production, modulation, and reception originated in the bilaterian stem group. This suggests that the monoaminergic system is a bilaterian novelty and that its evolution may have contributed to the Cambrian diversification.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a critical window into early vertebrate evolution. Here, we investigate the complex history, timing, and functional role of genome-wide duplications in vertebrates in the light of a chromosome-scale genome of the brown hagfish Eptatretus atami. Using robust chromosome-scale (paralogon-based) phylogenetic methods, we confirm the monophyly of cyclostomes, document an auto-tetraploidization (1RV) that predated the origin of crown group vertebrates ~517 Mya, and establish the timing of subsequent independent duplications in the gnathostome and cyclostome lineages. Some 1RV gene duplications can be linked to key vertebrate innovations, suggesting that this early genomewide event contributed to the emergence of pan-vertebrate features such as neural crest. The hagfish karyotype is derived by numerous fusions relative to the ancestral cyclostome arrangement preserved by lampreys. These genomic changes were accompanied by the loss of genes essential for organ systems (eyes, osteoclast) that are absent in hagfish, accounting in part for the simplification of the hagfish body plan; other gene family expansions account for hagfishes' capacity to produce slime. Finally, we characterise programmed DNA elimination in somatic cells of hagfish, identifying protein-coding and repetitive elements that are deleted during development. As in lampreys, the elimination of these genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline/pluripotency functions. Reconstruction of the early genomic history of vertebrates provides a framework for further exploration of vertebrate novelties.

Detection of interphylum transfers of the magnetosome gene cluster in magnetotactic bacteria

Magnetosome synthesis in magnetotactic bacteria (MTB) is regarded as a very ancient evolutionary process that dates back to deep-branching phyla. Magnetotactic bacteria belonging to one of such phyla, Nitrospirota, contain the classical genes for the magnetosome synthesis (e.g., mam, mms) and man genes, which were considered to be specific for this group. However, the recent discovery of man genes in MTB from the Thermodesulfobacteriota phylum has raised several questions about the inheritance of these genes in MTB. In this work, three new man genes containing MTB genomes affiliated with Nitrospirota and Thermodesulfobacteriota, were obtained. By applying reconciliation with these and the previously published MTB genomes, we demonstrate that the last common ancestor of all Nitrospirota was most likely not magnetotactic as assumed previously. Instead, our findings suggest that the genes for magnetosome synthesis were transmitted to the phylum Nitrospirota by horizontal gene transfer (HGT), which is the first case of the interphylum transfer of magnetosome genes detected to date. Furthermore, we provide evidence for the HGT of magnetosome genes from the Magnetobacteriaceae to the Dissulfurispiraceae family within Nitrospirota. Thus, our results imply a more significant role of HGT in the MTB evolution than deemed before and challenge the hypothesis of the ancient origin of magnetosome synthesis.

Thirdkind: displaying phylogenetic encounters beyond 2-level reconciliation

MOTIVATION

Reconciliation between a host and its symbiont phylogenies or between a species and a gene phylogenies is a prevalent approach in evolution, however no simple generic tool (i.e. virtually usable by all reconciliation software, from host/symbiont to species/gene comparisons) is available to visualize reconciliation results. Moreover there is no tool to visualize 3-levels reconciliations, i.e. to visualize 2 nested reconciliations as for example in a host/symbiont/gene complex.

RESULTS

Thirdkind is a light and easy to install command line software producing svg files displaying reconciliations, including 3-levels reconciliations. It takes a standard format recPhyloXML as input, and is thus usable with most reconciliation software.

AVAILABILITY AND IMPLEMENTATION

https://github.com/simonpenel/thirdkind/wiki.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Deciphering Microbial Gene Family Evolution Using Duplication-Transfer-Loss Reconciliation and RANGER-DTL

Phylogenetic reconciliation has emerged as a principled, highly effective technique for investigating the origin, spread, and evolutionary history of microbial gene families. Proper application of phylogenetic reconciliation requires a clear understanding of potential pitfalls and sources of error, and knowledge of the most effective reconciliation-based tools and protocols to use to maximize accuracy. In this book chapter, we provide a brief overview of Duplication-Transfer-Loss (DTL) reconciliation, the standard reconciliation model used to study microbial gene families and provide a step-by-step computational protocol to maximize the accuracy of DTL reconciliation and minimize false-positive evolutionary inferences.

DoubleRecViz: a web-based tool for visualizing transcript-gene-species tree reconciliation

MOTIVATION

A phylogenetic tree reconciliation is a mapping of one phylogenetic tree onto another which represents the co-evolution of two sets of taxa (e.g. parasite-host co-evolution, gene-species co-evolution). The reconciliation framework was extended to allow modeling the co-evolution of three sets of taxa such as transcript-gene-species co-evolutions. Several web-based tools have been developed for the display and manipulation of phylogenetic trees and co-phylogenetic trees involving two trees, but there currently exists no tool for visualizing the joint reconciliation between three phylogenetic trees.

RESULTS

Here, we present DoubleRecViz, a web-based tool for visualizing double reconciliations between phylogenetic trees at three levels: transcript, gene and species. DoubleRecViz extends the RecPhyloXML model-developed for gene-species tree reconciliation-to represent joint transcript-gene and gene-species tree reconciliations. It is implemented using the Dash library, which is a toolbox that provides dynamic visualization functionalities for web data visualization in Python.

AVAILABILITY AND IMPLEMENTATION

DoubleRecViz is available through a web server at https://doublerecviz.cobius.usherbrooke.ca. The source code and information about installation procedures are also available at https://github.com/UdeS-CoBIUS/DoubleRecViz.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

GeneRax: A Tool for Species-Tree-Aware Maximum Likelihood-Based Gene Family Tree Inference under Gene Duplication, Transfer, and Loss

Inferring phylogenetic trees for individual homologous gene families is difficult because alignments are often too short, and thus contain insufficient signal, while substitution models inevitably fail to capture the complexity of the evolutionary processes. To overcome these challenges, species-tree-aware methods also leverage information from a putative species tree. However, only few methods are available that implement a full likelihood framework or account for horizontal gene transfers. Furthermore, these methods often require expensive data preprocessing (e.g., computing bootstrap trees) and rely on approximations and heuristics that limit the degree of tree space exploration. Here, we present GeneRax, the first maximum likelihood species-tree-aware phylogenetic inference software. It simultaneously accounts for substitutions at the sequence level as well as gene level events, such as duplication, transfer, and loss relying on established maximum likelihood optimization algorithms. GeneRax can infer rooted phylogenetic trees for multiple gene families, directly from the per-gene sequence alignments and a rooted, yet undated, species tree. We show that compared with competing tools, on simulated data GeneRax infers trees that are the closest to the true tree in 90% of the simulations in terms of relative Robinson–Foulds distance. On empirical data sets, GeneRax is the fastest among all tested methods when starting from aligned sequences, and it infers trees with the highest likelihood score, based on our model. GeneRax completed tree inferences and reconciliations for 1,099 Cyanobacteria families in 8 min on 512 CPU cores. Thus, its parallelization scheme enables large-scale analyses. GeneRax is available under GNU GPL at https://github.com/BenoitMorel/GeneRax (last accessed June 17, 2020).  

Treerecs: an integrated phylogenetic tool, from sequences to reconciliations

MOTIVATION

Gene and species tree reconciliation methods are used to interpret gene trees, root them and correct uncertainties that are due to scarcity of signal in multiple sequence alignments. So far, reconciliation tools have not been integrated in standard phylogenetic software and they either lack performance on certain functions, or usability for biologists.

RESULTS

We present Treerecs, a phylogenetic software based on duplication-loss reconciliation. Treerecs is simple to install and to use. It is fast and versatile, has a graphic output, and can be used along with methods for phylogenetic inference on multiple alignments like PLL and Seaview.

AVAILABILITY AND IMPLEMENTATION

Treerecs is open-source. Its source code (C++, AGPLv3) and manuals are available from https://project.inria.fr/treerecs/.

Zombi: a phylogenetic simulator of trees, genomes and sequences that accounts for dead linages

Summary

Here we present Zombi, a tool to simulate the evolution of species, genomes and sequences in silico, that considers for the first time the evolution of genomes in extinct lineages. It also incorporates various features that have not to date been combined in a single simulator, such as the possibility of generating species trees with a pre-defined variation of speciation and extinction rates through time, simulating explicitly intergenic sequences of variable length and outputting gene tree—species tree reconciliations.

Availability and implementation

Source code and manual are freely available in https://github.com/AADavin/ZOMBI/.

Supplementary information

Supplementary data are available at Bioinformatics online.