Jungles: a new solution to the host/parasite phylogeny reconciliation problem

Strong phylogenetic congruence between Tulasnella fungi and their associated Drakaeinae orchids

The study of congruency between phylogenies of interacting species can provide a powerful approach for understanding the evolutionary history of symbiotic associations. Orchid mycorrhizal fungi can survive independently of orchids making cospeciation unlikely, leading us to predict that any congruence would arise from host-switches to closely related fungal species. The Australasian orchid subtribe Drakaeinae is an iconic group of sexually deceptive orchids that consists of approximately 66 species. In this study, we investigated the evolutionary relationships between representatives of all six Drakaeinae orchid genera (39 species) and their mycorrhizal fungi. We used an exome capture dataset to generate the first well-resolved phylogeny of the Drakaeinae genera. A total of 10 closely related Tulasnella Operational Taxonomic Units (OTUs) and previously described species were associated with the Drakaeinae orchids. Three of them were shared among orchid genera, with each genus associating with 1-6 Tulasnella lineages. Cophylogenetic analyses show Drakaeinae orchids and their Tulasnella associates exhibit significant congruence (p < 0.001) in the topology of their phylogenetic trees. An event-based method also revealed significant congruence in Drakaeinae-Tulasnella relationships, with duplications (35), losses (25), and failure to diverge (9) the most frequent events, with minimal evidence for cospeciation (1) and host-switches (2). The high number of duplications suggests that the orchids speciate independently from the fungi, and the fungal species association of the ancestral orchid species is typically maintained in the daughter species. For the Drakaeinae-Tulasnella interaction, a pattern of phylogenetic niche conservatism rather than coevolution likely explains the observed phylogenetic congruency in orchid and fungal phylogenies. Given that many orchid genera are characterized by sharing of fungal species between closely related orchid species, we predict that these findings may apply to a wide range of orchid lineages.

Divergent Viruses Discovered in Swine Alter the Understanding of Evolutionary History and Genetic Diversity of the Respirovirus Genus and Related Porcine Parainfluenza Viruses

Paramyxoviridae is a rapidly growing family of viruses, whose potential for cross-species transmission makes it difficult to predict the harm of newly emerging viruses to humans and animals. To better understand their diversity, evolutionary history, and co-evolution with their hosts, we analyzed a collection of porcine parainfluenza virus (PPIV) genomes to reconstruct the species classification basis and evolutionary history of the Respirovirus genus. We sequenced 17 complete genomes of porcine respirovirus 1 (also known as porcine parainfluenza virus 1; PPIV-1), thereby nearly tripling the number of currently available PPIV-1 genomes. We found that PPIV-1 was widely prevalent in China with two divergent lineages, PPIV-1a and PPIV-1b. We further provided evidence that a new species, porcine parainfluenza virus 2 (PPIV-2), had recently emerged in China. Our results pointed to a need for revising the current species demarcation criteria of the Respirovirus genus. In addition, we used PPIV-1 as an example to explore recombination and diversity of the Respirovirus genus. Interestingly, we only detected heterosubtypic recombination events between PPIV-1a and PPIV-1b with no intrasubtypic recombination events. The recombination hotspots highlighted a diverse geography-dependent genome structure of paramyxovirus infecting swine in China. Furthermore, we found no evidence of co-evolution between respirovirus and its host, indicating frequent cross-species transmission. In summary, our analyses showed that swine can be infected with a broad range of respiroviruses and recombination may serve as an important evolutionary mechanism for the Respirovirus genus’ greater diversity in genome structure than previously anticipated.

IMPORTANCE Livestock have emerged as critically underrecognized sources of paramyxovirus diversity, including pigs serving as the source of Nipah virus (NiV) and swine parainfluenza virus type 3, and goats and bovines harboring highly divergent viral lineages. Here, we identified a new species of Respirovirus genus named PPIV-2 in swine and proposed to revise the species demarcation criteria of the Respirovirus genus. We found heterosubtypic recombination events and high genetic diversity in PPIV-1. Further, we showed that genetic recombination may have occurred in the Respirovirus genus which may be associated with host range expansion. The continued expansion of Respirovirus genus diversity in livestock with relatively high human contact rates requires enhanced surveillance and ongoing evaluation of emerging cross-species transmission threats.

Linear Time Reconciliation With Bounded Transfers of Genes

Tree reconciliation is a general framework for investigating the mutual influence between gene and species trees according to the parsimony principle, that is, to each evolutionary event a cost is assigned and the goal is to find a reconciliation of minimum total cost. The resulting optimization problem is known as the reconciliation problem. Usually, the considered events are: co-divergence, gene Duplication, horizontal gene Transfer, and gene Loss (DTL model), while in a more conservative setting, gene transfers are not allowed (DL model). The reconciliation problem requires, in the DL model, time linear in the dimension of the two trees and at least quadratic time in the DTL model. Hence, it is reasonable to argue that the introduction of horizontal gene transfers increases the complexity of the problem. Instead, we introduce horizontal gene transfers with some constraints and prove that the problem is still linear in the dimension of the trees. Namely, we allow gene transfers of length bounded by k=2, on the basis of the observation that transfers are more likely to occur between closely related species than between distantly related ones. Then we extend the same reasonings to the case in which under additional constrains. In this paper we study also another problem related to the reconciliation one, that is optimally rooting one of the two trees when it is not, and also for it we prove similar results. The relevance of this contribution lies in showing that, in the transit from the DL to the DTL model, the computational time does not increase suddenly to quadratic but remains linear in the case when gene transfers are very short (i.e., happening between very close genes).

Efficiently sparse listing of classes of optimal cophylogeny reconciliations

Background

Cophylogeny reconciliation is a powerful method for analyzing host-parasite (or host-symbiont) co-evolution. It models co-evolution as an optimization problem where the set of all optimal solutions may represent different biological scenarios which thus need to be analyzed separately. Despite the significant research done in the area, few approaches have addressed the problem of helping the biologist deal with the often huge space of optimal solutions.

Results

In this paper, we propose a new approach to tackle this problem. We introduce three different criteria under which two solutions may be considered biologically equivalent, and then we propose polynomial-delay algorithms that enumerate only one representative per equivalence class (without listing all the solutions).

Conclusions

Our results are of both theoretical and practical importance. Indeed, as shown by the experiments, we are able to significantly reduce the space of optimal solutions while still maintaining important biological information about the whole space.

Primate malarias as a model for cross-species parasite transmission

Parasites regularly switch into new host species, representing a disease burden and conservation risk to the hosts. The distribution of these parasites also gives insight into characteristics of ecological networks and genetic mechanisms of host-parasite interactions. Some parasites are shared across many species, whereas others tend to be restricted to hosts from a single species. Understanding the mechanisms producing this distribution of host specificity can enable more effective interventions and potentially identify genetic targets for vaccines or therapies. As ecological connections between human and local animal populations increase, the risk to human and wildlife health from novel parasites also increases. Which of these parasites will fizzle out and which have the potential to become widespread in humans? We consider the case of primate malarias, caused by Plasmodium parasites, to investigate the interacting ecological and evolutionary mechanisms that put human and nonhuman primates at risk for infection. Plasmodium host switching from nonhuman primates to humans led to ancient introductions of the most common malaria-causing agents in humans today, and new parasite switching is a growing threat, especially in Asia and South America. Based on a wild host-Plasmodium occurrence database, we highlight geographic areas of concern and potential areas to target further sampling. We also discuss methodological developments that will facilitate clinical and field-based interventions to improve human and wildlife health based on this eco-evolutionary perspective.

Identifying duplications and lateral gene transfers simultaneously and rapidly

This paper deals with the problem of enumerating all minimum-cost LCA-reconciliations involving gene duplications and lateral gene transfers (LGTs) for a given species tree [Formula: see text] and a given gene tree [Formula: see text]. Previously, [Tofigh A, Hallett M, Lagergren J, Simultaneous identification of duplications and lateral gene transfers, IEEE/ACM Trans Comput Biol Bioinf 517-535, 2011.] gave a fixed-parameter algorithm for this problem that runs in [Formula: see text] time, where [Formula: see text] is the number of vertices in [Formula: see text], [Formula: see text] is the number of vertices in [Formula: see text], and [Formula: see text] is the minimum cost of an LCA-reconciliation between [Formula: see text] and [Formula: see text]. In this paper, by refining their algorithm, we obtain a new one for the same problem that finds and outputs the solutions in a compact form within [Formula: see text] time. In the most interesting case where [Formula: see text], our algorithm is [Formula: see text] times faster.

Indirect identification of horizontal gene transfer

Several implicit methods to infer horizontal gene transfer (HGT) focus on pairs of genes that have diverged only after the divergence of the two species in which the genes reside. This situation defines the edge set of a graph, the later-divergence-time (LDT) graph, whose vertices correspond to genes colored by their species. We investigate these graphs in the setting of relaxed scenarios, i.e., evolutionary scenarios that encompass all commonly used variants of duplication-transfer-loss scenarios in the literature. We characterize LDT graphs as a subclass of properly vertex-colored cographs, and provide a polynomial-time recognition algorithm as well as an algorithm to construct a relaxed scenario that explains a given LDT. An edge in an LDT graph implies that the two corresponding genes are separated by at least one HGT event. The converse is not true, however. We show that the complete xenology relation is described by an rs-Fitch graph, i.e., a complete multipartite graph satisfying constraints on the vertex coloring. This class of vertex-colored graphs is also recognizable in polynomial time. We finally address the question “how much information about all HGT events is contained in LDT graphs” with the help of simulations of evolutionary scenarios with a wide range of duplication, loss, and HGT events. In particular, we show that a simple greedy graph editing scheme can be used to efficiently detect HGT events that are implicitly contained in LDT graphs.

A broadscale analysis of host-symbiont cophylogeny reveals the drivers of phylogenetic congruence

Symbioses exert substantial biological influence, with great evolutionary and ecological relevance for disease, major evolutionary transitions, and the structure and function of ecological communities. Yet, much remains unknown about the patterns and processes that characterise symbioses. A major unanswered question is the extent to which symbiont phylogenies mirror those of their hosts and if patterns differ for parasites and mutualists. Addressing this question offers fundamental insights into evolutionary processes, such as whether symbionts typically codiverge with their hosts or if diversity is generated via host switches. Here, we perform a meta-analysis of host-symbiont phylogenetic congruence, encompassing 212 host-symbiont cophylogenetic studies that include ~10,000 species. Our analysis supersedes previous qualitative assessments by utilising a quantitative framework. We show that symbiont phylogeny broadly reflects host phylogeny across biodiversity and life-history, demonstrating a general pattern of phylogenetic congruence in host-symbiont interactions. We reveal two key aspects of symbiont life-history that promote closer ties between hosts and symbionts: vertical transmission and mutualism. Mode of symbiosis and mode of transmission are intimately interlinked, but vertical transmission is the dominant factor. Given the pervasiveness of symbioses, these findings provide important insights into the processes responsible for generating and maintaining the Earth's rich biodiversity.

Capybara: equivalence ClAss enumeration of coPhylogenY event-BAsed ReconciliAtions

MOTIVATION

Phylogenetic tree reconciliation is the method of choice in analyzing host-symbiont systems. Despite the many reconciliation tools that have been proposed in the literature, two main issues remain unresolved: (i) listing suboptimal solutions (i.e. whose score is 'close' to the optimal ones) and (ii) listing only solutions that are biologically different 'enough'. The first issue arises because the optimal solutions are not always the ones biologically most significant; providing many suboptimal solutions as alternatives for the optimal ones is thus very useful. The second one is related to the difficulty to analyze an often huge number of optimal solutions. In this article, we propose Capybara that addresses both of these problems in an efficient way. Furthermore, it includes a tool for visualizing the solutions that significantly helps the user in the process of analyzing the results.

AVAILABILITY AND IMPLEMENTATION

The source code, documentation and binaries for all platforms are freely available at https://capybara-doc.readthedocs.io/.

CONTACT

yishu.wang@univ-lyon1.fr or blerina.sinaimeri@inria.fr.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Wolbachia infection in wild mosquitoes (Diptera: Culicidae): implications for transmission modes and host-endosymbiont associations in Singapore

Background

Wolbachia are intracellular bacterial endosymbionts found in most insect lineages. In mosquitoes, the influence of these endosymbionts on host reproduction and arboviral transmission has spurred numerous studies aimed at using Wolbachia infection as a vector control technique. However, there are several knowledge gaps in the literature and little is known about natural Wolbachia infection across species, their transmission modes, or associations between various Wolbachia lineages and their hosts. This study aims to address these gaps by exploring mosquito-Wolbachia associations and their evolutionary implications.

Methods

We conducted tissue-specific polymerase chain reaction screening for Wolbachia infection in the leg, gut and reproductive tissues of wild mosquitoes from Singapore using the Wolbachia surface protein gene (wsp) molecular marker. Mosquito-Wolbachia associations were explored using three methods—tanglegram, distance-based, and event-based methods—and by inferred instances of vertical transmission and host shifts.

Results

Adult mosquitoes (271 specimens) representing 14 genera and 40 species were screened for Wolbachia. Overall, 21 species (51.2%) were found positive for Wolbachia, including five in the genus Aedes and five in the genus Culex. To our knowledge, Wolbachia infections have not been previously reported in seven of these 21 species: Aedes nr. fumidus, Aedes annandalei, Uranotaenia obscura, Uranotaenia trilineata, Verrallina butleri, Verrallina sp. and Zeugnomyia gracilis. Wolbachia were predominantly detected in the reproductive tissues, which is an indication of vertical transmission. However, Wolbachia infection rates varied widely within a mosquito host species. There was no clear signal of cophylogeny between the mosquito hosts and the 12 putative Wolbachia strains observed in this study. Host shift events were also observed.

Conclusions

Our results suggest that the mosquito-Wolbachia relationship is complex and that combinations of transmission modes and multiple evolutionary events likely explain the observed distribution of Wolbachia diversity across mosquito hosts. These findings have implications for a better understanding of the diversity and ecology of Wolbachia and for their utility as biocontrol agents.

A New Paradigm for Identifying Reconciliation-Scenario Altering Mutations Conferring Environmental Adaptation

An important goal in microbial computational genomics is to identify crucial events in the evolution of a gene that severely alter the duplication, loss, and mobilization patterns of the gene within the genomes in which it disseminates. In this article, we formalize this microbiological goal as a new pattern-matching problem in the domain of gene tree and species tree reconciliation, denoted "Reconciliation-Scenario Altering Mutation (RSAM) Discovery." We propose an [Formula: see text] time algorithm to solve this new problem, where m and n are the number of vertices of the input gene tree and species tree, respectively, and k is a user-specified parameter that bounds from above the number of optimal solutions of interest. The algorithm first constructs a hypergraph representing the k highest scoring reconciliation scenarios between the given gene tree and species tree, and then interrogates this hypergraph for subtrees matching a prespecified RSAM pattern. Our algorithm is optimal in the sense that the number of hypernodes in the hypergraph can be lower bounded by [Formula: see text]. We implement the new algorithm as a tool, called RSAM-finder, and demonstrate its application to the identification of RSAMs in toxins and drug resistance elements across a data set spanning hundreds of species.

Molecular fossils reveal ancient associations of dsDNA viruses with several phyla of fungi

Little is known about the infections of double-stranded DNA (dsDNA) viruses in fungi. Here, we use a paleovirological method to systematically identify the footprints of past dsDNA virus infections within the fungal genomes. We uncover two distinct groups of endogenous nucleocytoplasmic large DNA viruses (NCLDVs) in at least seven fungal phyla (accounting for about a third of known fungal phyla), revealing an unprecedented diversity of dsDNA viruses in fungi. Interestingly, one fungal dsDNA virus lineage infecting six fungal phyla is closely related to the giant virus Pithovirus, suggesting giant virus relatives might widely infect fungi. Co-speciation analyses indicate fungal NCLDVs mainly evolved through cross-species transmission. Taken together, our findings provide novel insights into the diversity and evolution of NCLDVs in fungi.

Inferring Pareto-optimal reconciliations across multiple event costs under the duplication-loss-coalescence model

BACKGROUND

Reconciliation methods are widely used to explain incongruence between a gene tree and species tree. However, the common approach of inferring maximum parsimony reconciliations (MPRs) relies on user-defined costs for each type of event, which can be difficult to estimate. Prior work has explored the relationship between event costs and maximum parsimony reconciliations in the duplication-loss and duplication-transfer-loss models, but no studies have addressed this relationship in the more complicated duplication-loss-coalescence model.

RESULTS

We provide a fixed-parameter tractable algorithm for computing Pareto-optimal reconciliations and recording all events that arise in those reconciliations, along with their frequencies. We apply this method to a case study of 16 fungi to systematically characterize the complexity of MPR space across event costs and identify events supported across this space.

CONCLUSION

This work provides a new framework for studying the relationship between event costs and reconciliations that incorporates both macro-evolutionary events and population effects and is thus broadly applicable across eukaryotic species.

Co-divergence and tree topology

In reconstructing the common evolutionary history of hosts and parasites, the current method of choice is the phylogenetic tree reconciliation. In this model, we are given a host tree H, a parasite tree P, and a function [Formula: see text] mapping the leaves of P to the leaves of H and the goal is to find, under some biologically motivated constraints, a reconciliation, that is a function from the vertices of P to the vertices of H that respects [Formula: see text] and allows the identification of biological events such as co-speciation, duplication and host switch. The maximum co-divergence problem consists in finding the maximum number of co-speciations in a reconciliation. This problem is NP-hard for arbitrary phylogenetic trees and no approximation algorithm is known. In this paper we consider the influence of tree topology on the maximum co-divergence problem. In particular we focus on a particular tree structure, namely caterpillar, and show that-in this case-the heuristics that are mostly used in the literature provide solutions that can be arbitrarily far from the optimal value. Then, we prove that finding the max co-divergence is equivalent to compute the maximum length of a subsequence with certain properties of a given permutation. This equivalence leads to two consequences: (1) it shows that we can compute efficiently in polynomial time the optimal time-feasible reconciliation and (2) it can be used to understand how much the tree topology influences the value of the maximum number of co-speciations.

Inferring Gene-Species Assignments in the Presence of Horizontal Gene Transfer

BACKGROUND

Microbial communities from environmental samples show great diversity as bacteria quickly responds to changes in their ecosystems. To assess the scenario of the actual changes, metagenomics experiments aimed at sequencing genomic DNA from such samples are performed. These new obtained sequences together with already known are used to infer phylogenetic trees assessing the taxonomic groups the species with these genes belong to. Here, we propose the first approach to the gene-species assignment problem by using reconciliation with horizontal gene transfer.

RESULTS

We propose efficient algorithms that search for optimal gene-species mappings taking into account gene duplication, loss and transfer events under two tractable models of HGT reconciliation.

CONCLUSIONS

We calculate both the optimal cost and all possible optimal scenarios. Furthermore as the number of optimal reconstructions can be large, we use a Monte-Carlo method for the inference of approximate distributions of gene-species assignments. We demonstrate the applicability on empirical and simulated datasets.

Insights into Circovirus Host Range from the Genomic Fossil Record

A diverse range of DNA sequences derived from circoviruses (family Circoviridae) has been identified in samples obtained from humans and domestic animals, often in association with pathological conditions. In the majority of cases, however, little is known about the natural biology of the viruses from which these sequences are derived. Endogenous circoviral elements (CVe) are DNA sequences derived from circoviruses that occur in animal genomes and provide a useful source of information about circovirus-host relationships. In this study, we screened genome assemblies of 675 animal species and identified numerous circovirus-related sequences, including the first examples of CVe derived from cycloviruses. We confirmed the presence of these CVe in the germ line of the elongate twig ant (Pseudomyrmex gracilis), thereby establishing that cycloviruses infect insects. We examined the evolutionary relationships between CVe and contemporary circoviruses, showing that CVe from ants and mites group relatively closely with cycloviruses in phylogenies. Furthermore, the relatively random interspersion of CVe from insect genomes with cyclovirus sequences recovered from vertebrate samples suggested that contamination might be an important consideration in studies reporting these viruses. Our study demonstrates how endogenous viral sequences can inform metagenomics-based virus discovery. In addition, it raises doubts about the role of cycloviruses as pathogens of humans and other vertebrates.IMPORTANCE Advances in DNA sequencing have dramatically increased the rate at which new viruses are being identified. However, the host species associations of most virus sequences identified in metagenomic samples are difficult to determine. Our analysis indicates that viruses proposed to infect vertebrates (in some cases being linked to human disease) may in fact be restricted to arthropod hosts. The detection of these sequences in vertebrate samples may reflect their widespread presence in the environment as viruses of parasitic arthropods.

Exploring the Robustness of the Parsimonious Reconciliation Method in Host-Symbiont Cophylogeny

The aim of this paper is to explore the robustness of the parsimonious host-symbiont tree reconciliation method under editing or small perturbations of the input. The editing involves making different choices of unique symbiont mapping to a host in the case where multiple associations exist. This is made necessary by the fact that the tree reconciliation model is currently unable to handle such associations. The analysis performed could however also address the problem of errors. The perturbations are re-rootings of the symbiont tree to deal with a possibly wrong placement of the root specially in the case of fast-evolving species. In order to do this robustness analysis, we introduce a simulation scheme specifically designed for the host-symbiont cophylogeny context, as well as a measure to compare sets of tree reconciliations, both of which are of interest by themselves.

Euphyllophyte Paleoviruses Illuminate Hidden Diversity and Macroevolutionary Mode of Caulimoviridae

Endogenous viral elements (paleoviruses) provide "molecular fossils" for studying the deep history and macroevolution of viruses. Endogenous plant pararetroviruses (EPRVs) are widespread in angiosperms, but little is known about EPRVs in earlier-branching plants. Here we use a large-scale phylogenomic approach to investigate the diversity and macroevolution of plant pararetroviruses (formally known as Caulimoviridae). We uncover an unprecedented and unappreciated diversity of EPRVs within the genomes of gymnosperms and ferns. The known angiosperm viruses constitute only a minor part of the Caulimoviridae diversity. By characterizing the distribution of EPRVs, we show that no major euphyllophyte lineages escape the activity of Caulimoviridae, raising the possibility that many exogenous Caulimoviridae remain to be discovered in euphyllophytes. We find that the copy numbers of EPRVs are generally high, suggesting that EPRVs might define a unique group of repetitive elements and represent important components of euphyllophyte genomes. Evolutionary analyses suggest an ancient origin of Caulimoviridae and at least three independent origins of Caulimoviridae in angiosperms. Our findings reveal the remarkable diversity of Caulimoviridae and have important implications for understanding the origin and macroevolution of plant pararetroviruses.IMPORTANCE Few viruses have been documented in plants outside angiosperms. Viruses can occasionally integrate into host genomes, forming endogenous viral elements (EVEs). Endogenous plant pararetroviruses (EPRVs) are widespread in angiosperms. In this study, we performed comprehensive comparative and phylogenetic analyses of EPRVs and found that EPRVs are present in the genomes of gymnosperms and ferns. We identified numerous EPRVs in gymnosperm and fern genomes, revealing an unprecedented depth in the diversity of plant pararetroviruses. Plant pararetroviruses mainly underwent cross-species transmission, and angiosperm pararetroviruses arose at least three times. Our study provides novel insights into the diversity and macroevolution of plant pararetroviruses.

A geography-aware reconciliation method to investigate diversification patterns in host/parasite interactions

Cospeciation studies aim at investigating whether hosts and symbionts speciate simultaneously or whether the associations diversify through host shifts. This problem is often tackled through reconciliation analyses that map the symbiont phylogeny onto the host phylogeny by mixing different types of diversification events. These reconciliations can be difficult to interpret and are not always biologically realistic. Researchers have underlined that the biogeographic histories of both hosts and symbionts influence the probability of cospeciation and host switches, but up to now no reconciliation software integrates geographic data. We present a new functionality in the Mowgli software that bridges this gap. The user can provide geographic information on both the host and symbiont extant and ancestral taxa. Constraints in the reconciliation algorithm have been implemented to generate biologically realistic codiversification scenarios. We apply our method to the fig/fig wasp association and infer diversification scenarios that differ from reconciliations ignoring geographic information. In addition, we updated the reconciliation viewer SylvX to visualize ancestral character states on the phylogenetic trees and highlight parts of reconciliations that are geographically inconsistent when not accounting for geographic constraints. We suggest that the comparison of reconciliations obtained with and without such constraints can help solving ambiguities in the biogeographic histories of the partners. With the development of robust methods in historical biogeography, and the advent of next-generation sequencing that leads to better-resolved trees, a geography-aware reconciliation method represents a substantial advance that is likely to be useful to researchers studying the evolution of biotic interactions and biogeography.

Host Specificity in Subarctic Aphids

Plants and herbivorous (or parasitic) insects form the majority of macroscopic life. The specificity of interaction between host plant and parasitic insect depends on the adaptations of both the host and the parasite. Over time, these interactions evolve and change as a result of an 'arms race' between host and parasite, and the resulting species-specific adaptations may be maintained, perpetuating these interactions across speciation events. This can lead to specialisation between species or clades. With speciation and species sorting over time, complex interactions evolve. Here, we elucidate a three-tier method to test these interactions using the aphids (Hemiptera: Aphididae) and plants of Churchill (Manitoba, Canada) as a model system. We analyzed these interactions by testing for three patterns in host specificity: monophagy, phylogenetic clustering, and cophylogeny. We defined monophagy strictly as one species feeding exclusively upon a single host plant species (an association likely driven by arms races in morphology, chemical resistance/tolerance, and visual appearance) and observed this in 7 of 22 aphid species. In all the remaining 'polyphagous' cases, there was a strong trend toward monophagy (80% of individuals were found on a single host plant species). Second, we observed two separate examples of phylogenetic clustering where groups of closely related aphid species fed upon individual plant species. Finally, we found no support for cophylogenetic relationships where both aphids and plants cospeciate to form congruent phylogenetic trees (evidence of coadaptation through an ongoing arms race). One explanation for uncovering species-specific interactions in a recently deglaciated, subarctic locality is that the species involved in the associations moved north together. Testing different levels of specificity in the most predominant species-species interactions on the planet will allow us to elucidate these patterns accurately and gives us insight into where to direct future research.

Time-consistent reconciliation maps and forbidden time travel

Background

In the absence of horizontal gene transfer it is possible to reconstruct the history of gene families from empirically determined orthology relations, which are equivalent to event-labeled gene trees. Knowledge of the event labels considerably simplifies the problem of reconciling a gene tree T with a species trees S, relative to the reconciliation problem without prior knowledge of the event types. It is well-known that optimal reconciliations in the unlabeled case may violate time-consistency and thus are not biologically feasible. Here we investigate the mathematical structure of the event labeled reconciliation problem with horizontal transfer.

Results

We investigate the issue of time-consistency for the event-labeled version of the reconciliation problem, provide a convenient axiomatic framework, and derive a complete characterization of time-consistent reconciliations. This characterization depends on certain weak conditions on the event-labeled gene trees that reflect conditions under which evolutionary events are observable at least in principle. We give an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {O}(|V(T)|\log (|V(S)|))$$\end{document}O(|V(T)|log(|V(S)|))-time algorithm to decide whether a time-consistent reconciliation map exists. It does not require the construction of explicit timing maps, but relies entirely on the comparably easy task of checking whether a small auxiliary graph is acyclic. The algorithms are implemented in C++ using the boost graph library and are freely available at https://github.com/Nojgaard/tc-recon.

Significance

The combinatorial characterization of time consistency and thus biologically feasible reconciliation is an important step towards the inference of gene family histories with horizontal transfer from orthology data, i.e., without presupposed gene and species trees. The fast algorithm to decide time consistency is useful in a broader context because it constitutes an attractive component for all tools that address tree reconciliation problems.

Tanglegrams: A Reduction Tool for Mathematical Phylogenetics

Many discrete mathematics problems in phylogenetics are defined in terms of the relative labeling of pairs of leaf-labeled trees. These relative labelings are naturally formalized as tanglegrams, which have previously been an object of study in coevolutionary analysis. Although there has been considerable work on planar drawings of tanglegrams, they have not been fully explored as combinatorial objects until recently. In this paper, we describe how many discrete mathematical questions on trees "factor" through a problem on tanglegrams, and how understanding that factoring can simplify analysis. Depending on the problem, it may be useful to consider a unordered version of tanglegrams, and/or their unrooted counterparts. For all of these definitions, we show how the isomorphism types of tanglegrams can be understood in terms of double cosets of the symmetric group, and we investigate their automorphisms. Understanding tanglegrams better will isolate the distinct problems on leaf-labeled pairs of trees and reveal natural symmetries of spaces associated with such problems.

Biologically feasible gene trees, reconciliation maps and informative triples

BACKGROUND

The history of gene families-which are equivalent to event-labeled gene trees-can be reconstructed from empirically estimated evolutionary event-relations containing pairs of orthologous, paralogous or xenologous genes. The question then arises as whether inferred event-labeled gene trees are biologically feasible, that is, if there is a possible true history that would explain a given gene tree. In practice, this problem is boiled down to finding a reconciliation map-also known as DTL-scenario-between the event-labeled gene trees and a (possibly unknown) species tree.

RESULTS

In this contribution, we first characterize whether there is a valid reconciliation map for binary event-labeled gene trees T that contain speciation, duplication and horizontal gene transfer events and some unknown species tree S in terms of "informative" triples that are displayed in T and provide information of the topology of S. These informative triples are used to infer the unknown species tree S for T. We obtain a similar result for non-binary gene trees. To this end, however, the reconciliation map needs to be further restricted. We provide a polynomial-time algorithm to decide whether there is a species tree for a given event-labeled gene tree, and in the positive case, to construct the species tree and the respective (restricted) reconciliation map. However, informative triples as well as DTL-scenarios have their limitations when they are used to explain the biological feasibility of gene trees. While reconciliation maps imply biological feasibility, we show that the converse is not true in general. Moreover, we show that informative triples neither provide enough information to characterize "relaxed" DTL-scenarios nor non-restricted reconciliation maps for non-binary biologically feasible gene trees.

Global patterns in coronavirus diversity

Since the emergence of Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV) and Middle East Respiratory Syndrom Coronavirus (MERS-CoV) it has become increasingly clear that bats are important reservoirs of CoVs. Despite this, only 6% of all CoV sequences in GenBank are from bats. The remaining 94% largely consist of known pathogens of public health or agricultural significance, indicating that current research effort is heavily biased towards describing known diseases rather than the ‘pre-emergent’ diversity in bats. Our study addresses this critical gap, and focuses on resource poor countries where the risk of zoonotic emergence is believed to be highest. We surveyed the diversity of CoVs in multiple host taxa from twenty countries to explore the factors driving viral diversity at a global scale. We identified sequences representing 100 discrete phylogenetic clusters, ninety-one of which were found in bats, and used ecological and epidemiologic analyses to show that patterns of CoV diversity correlate with those of bat diversity. This cements bats as the major evolutionary reservoirs and ecological drivers of CoV diversity. Co-phylogenetic reconciliation analysis was also used to show that host switching has contributed to CoV evolution, and a preliminary analysis suggests that regional variation exists in the dynamics of this process. Overall our study represents a model for exploring global viral diversity and advances our fundamental understanding of CoV biodiversity and the potential risk factors associated with zoonotic emergence.

Fast algorithm for the reconciliation of gene trees and LGT networks

In phylogenomics, reconciliations aim at explaining the discrepancies between the evolutionary histories of genes and species. Several reconciliation models are available when the evolution of the species of interest is modelled via phylogenetic trees; the most commonly used are the DL model, accounting for duplications and losses in gene evolution and yielding polynomially-solvable problems, and the DTL model, which also accounts for gene transfers and implies NP-hard problems. However, when dealing with non-tree-like evolutionary events such as hybridisations, phylogenetic networks - and not phylogenetic trees - should be used to model species evolution. Reconciliation models involving phylogenetic networks are still at their early days. In this paper, we propose a new reconciliation model in which the evolution of species is modelled by a special kind of phylogenetic networks - the LGT networks. Our model considers duplications, losses and transfers of genes, but restricts transfers to happen through some specific arcs of the network, called secondary arcs. Moreover, we provide a polynomial algorithm to compute the most parsimonious reconciliation between a gene tree and an LGT network under this model. Our method, when combined with quartet decomposition methods to detect putative "highways" of transfers, permits to refine their analyses by allowing to examine the two possible directions of a highway and even consider combinations of highways.

Whence river blindness? The domestication of mammals and host-parasite co-evolution in the nematode genus Onchocerca

The genus Onchocerca includes 34 described species and represents one of the largest genera of the filarial nematodes within the family Onchocercidae. Representative members of this genus are mainly parasites of ungulates, with some exceptions such as Onchocerca lupi and Onchocerca volvulus, infecting carnivores and/or humans. For a long time, the evolutionary relationships amongst onchocercids remained poorly studied, as the systematics of this genus was impaired by the high morphological variability of species included in the taxon. Although some molecular phylogenies were developed, these studies were mainly focused on bovine Onchocerca spp. and O. volvulus, including assessments of Wolbachia endosymbionts. In the present study, we analysed 13 Onchocerca spp. from a larger host spectrum using a panel of seven different genes. Analysis of the coxI marker supports its usefulness for the identification of species within the genus. The evolutionary history of the genus has been herein revised by multi-gene phylogenies, presenting three strongly supported clades of Onchocerca spp. Analyses of co-evolutionary scenarios between Onchocerca and their vertebrate hosts underline the effect of domestication on Onchocerca speciation. Our study indicates that a host switch event occurred between Bovidae, Canidae and humans. Cophylogenetic analyses between Onchocerca and the endosymbiotic bacterium Wolbachia indicate the strongest co-evolutionary pattern ever registered within the filarial nematodes. Finally, this dataset indicates that the clade composed by O. lupi, Onchocerca gutturosa, Onchocerca lienalis, Onchocerca ochengi and O. volvulus derived from recent speciation.

Testing cospeciation through large-scale cophylogenetic studies

Insects are involved in a multitude of interactions with other organisms, which make them ideal models for large-scale cophylogenetic studies. Once phylogenies of interacting lineages have been inferred, there are a number of questions we may wish to ask, such as what was the relationship between the partners in the past? Have they co-evolved for thousands or millions of years, or has one of the partners switched among different host species? To answer such questions, researchers may conduct cophylogenetic analysis, to explore the relationships between the phylogenies of interacting lineages and determine whether the match is significant or find explanations for observed differences. When combined with dating analyses, cophylogenetic analyses may support cospeciation of the partners or phylogenetic tracking. As they may reveal dynamics of host-pathogen coevolution, cophylogenetic studies may also help tackle global health issues (e.g. document the spread of disease causing pathogens). Cophylogenetic studies of parasitoids and their insect hosts may also help identify effective biocontrol agents. With the advent of next generation sequencing technologies and keeping in mind that systematic errors may occur, cophylogenetics will benefit from better-resolved trees, allowing more accurate reconciliation. However as trees become larger, current algorithms also become more computationally challenging. Nevertheless, both theoretical and methodological developments are leading to more accurate and powerful tests of cospeciation through cophylogenetic analysis.

Nematode endoparasites do not codiversify with their stick insect hosts

Host-parasite coevolution stems from reciprocal selection on host resistance and parasite infectivity, and can generate some of the strongest selective pressures known in nature. It is widely seen as a major driver of diversification, the most extreme case being parallel speciation in hosts and their associated parasites. Here, we report on endoparasitic nematodes, most likely members of the mermithid family, infecting different Timema stick insect species throughout California. The nematodes develop in the hemolymph of their insect host and kill it upon emergence, completely impeding host reproduction. Given the direct exposure of the endoparasites to the host's immune system in the hemolymph, and the consequences of infection on host fitness, we predicted that divergence among hosts may drive parallel divergence in the endoparasites. Our phylogenetic analyses suggested the presence of two differentiated endoparasite lineages. However, independently of whether the two lineages were considered separately or jointly, we found a complete lack of codivergence between the endoparasitic nematodes and their hosts in spite of extensive genetic variation among hosts and among parasites. Instead, there was strong isolation by distance among the endoparasitic nematodes, indicating that geography plays a more important role than host-related adaptations in driving parasite diversification in this system. The accumulating evidence for lack of codiversification between parasites and their hosts at macroevolutionary scales contrasts with the overwhelming evidence for coevolution within populations, and calls for studies linking micro- versus macroevolutionary dynamics in host-parasite interactions.

Phylogenetic framework for coevolutionary studies: a compass for exploring jungles of tangled trees

Phylogenetics is used to detect past evolutionary events, from how species originated to how their ecological interactions with other species arose, which can mirror cophylogenetic patterns. Cophylogenetic reconstructions uncover past ecological relationships between taxa through inferred coevolutionary events on trees, for example, codivergence, duplication, host-switching, and loss. These events can be detected by cophylogenetic analyses based on nodes and the length and branching pattern of the phylogenetic trees of symbiotic associations, for example, host-parasite. In the past 2 decades, algorithms have been developed for cophylogetenic analyses and implemented in different software, for example, statistical congruence index and event-based methods. Based on the combination of these approaches, it is possible to integrate temporal information into cophylogenetical inference, such as estimates of lineage divergence times between 2 taxa, for example, hosts and parasites. Additionally, the advances in phylogenetic biogeography applying methods based on parametric process models and combined Bayesian approaches, can be useful for interpreting coevolutionary histories in a scenario of biogeographical area connectivity through time. This article briefly reviews the basics of parasitology and provides an overview of software packages in cophylogenetic methods. Thus, the objective here is to present a phylogenetic framework for coevolutionary studies, with special emphasis on groups of parasitic organisms. Researchers wishing to undertake phylogeny-based coevolutionary studies can use this review as a "compass" when "walking" through jungles of tangled phylogenetic trees.

Inferring phylogeny and speciation of Gymnosporangium species, and their coevolution with host plants

Gymnosporangium species (Pucciniaceae, Pucciniales) cause serious diseases and significant economic losses to apple cultivars. Most of the reported species are heteroecious and complete their life cycles on two different plant hosts belonging to two unrelated genera, i.e. Juniperus and Malus. However, the phylogenetic relationships among Gymnosporangium species and the evolutionary history of Gymnosporangium on its aecial and telial hosts were still undetermined. In this study, we recognized species based on rDNA sequence data by using coalescent method of generalized mixed Yule-coalescent (GMYC) and Poisson Tree Processes (PTP) models. The evolutionary relationships of Gymnosporangium species and their hosts were investigated by comparing the cophylogenetic analyses of Gymnosporangium species with Malus species and Juniperus species, respectively. The concordant results of GMYC and PTP analyses recognized 14 species including 12 known species and two undescribed species. In addition, host alternations of 10 Gymnosporangium species were uncovered by linking the derived sequences between their aecial and telial stages. This study revealed the evolutionary process of Gymnosporangium species, and clarified that the aecial hosts played more important roles than telial hosts in the speciation of Gymnosporangium species. Host switch, losses, duplication and failure to divergence all contributed to the speciation of Gymnosporangium species.

Host and ecology both play a role in shaping distribution of digenean parasites of New Zealand whelks (Gastropoda: Buccinidae: Cominella)

Digenean parasites infecting four Cominella whelk species (C. glandiformis, C. adspersa, C. maculosa and C. virgata), which inhabit New Zealand's intertidal zone, were analysed using molecular techniques. Mitochondrial 16S and cytochrome oxidase 1 (COI) and nuclear rDNA ITS1 sequences were used to infer phylogenetic relationships amongst digenea. Host species were parasitized by a diverse range of digenea (Platyhelminthes, Trematoda), representing seven families: Echinostomatidae, Opecoelidae, Microphallidae, Strigeidae and three, as yet, undetermined families A, B and C. Each parasite family infected between one and three host whelk species, and infection levels were typically low (average infection rates ranged from 1·4 to 3·6%). Host specificity ranged from highly species-specific amongst the echinostomes, which were only ever observed infecting C. glandiformis, to the more generalist opecoelids and strigeids, which were capable of infecting three out of four of the Cominella species analysed. Digeneans displayed a highly variable geographic range; for example, echinostomes had a large geographic range stretching the length of New Zealand, from Northland to Otago, whereas Family B parasites were restricted to fairly small areas of the North Island. Our results add to a growing body of research identifying wide ranges in both host specificity and geographic range amongst intertidal, multi-host parasite systems.

Towards sub-quadratic time and space complexity solutions for the dated tree reconciliation problem

BACKGROUND

Recent coevolutionary analysis has considered tree topology as a means to reduce the asymptotic complexity associated with inferring the complex coevolutionary interrelationships that arise between phylogenetic trees. Targeted algorithmic design for specific tree topologies has to date been highly successful, with one recent formulation providing a logarithmic space complexity reduction for the dated tree reconciliation problem.

METHODS

In this work we build on this prior analysis providing a further asymptotic space reduction, by providing a new formulation for the dynamic programming table used by a number of popular coevolutionary analysis techniques. This model gives rise to a sub quadratic running time solution for the dated tree reconciliation problem for selected tree topologies, and is shown to be, in practice, the fastest method for solving the dated tree reconciliation problem for expected evolutionary trees. This result is achieved through the analysis of not only the topology of the trees considered for coevolutionary analysis, but also the underlying structure of the dynamic programming algorithms that are traditionally applied to such analysis.

CONCLUSION

The newly inferred theoretical complexity bounds introduced herein are then validated using a combination of synthetic and biological data sets, where the proposed model is shown to provide an [Formula: see text] space saving, while it is observed to run in half the time compared to the fastest known algorithm for solving the dated tree reconciliation problem. What is even more significant is that the algorithm derived herein is able to guarantee the optimality of its inferred solution, something that algorithms of comparable speed have to date been unable to achieve.

Metapopulation dominance and genomic-island acquisition of Bradyrhizobium with superior catabolic capabilities

Root nodule-forming rhizobia exhibit a bipartite lifestyle, replicating in soil and also within plant cells where they fix nitrogen for legume hosts. Host control models posit that legume hosts act as a predominant selective force on rhizobia, but few studies have examined rhizobial fitness in natural populations. Here, we genotyped and phenotyped Bradyrhizobium isolates across more than 800 km of the native Acmispon strigosus host range. We sequenced chromosomal genes expressed under free-living conditions and accessory symbiosis loci expressed in planta and encoded on an integrated 'symbiosis island' (SI). We uncovered a massive clonal expansion restricted to the Bradyrhizobium chromosome, with a single chromosomal haplotype dominating populations, ranging more than 700 km, and acquiring 42 divergent SI haplotypes, none of which were spatially widespread. For focal genotypes, we quantified utilization of 190 sole-carbon sources relevant to soil fitness. Chromosomal haplotypes that were both widespread and dominant exhibited superior growth on diverse carbon sources, whereas these patterns were not mirrored among SI haplotypes. Abundance, spatial range and catabolic superiority of chromosomal, but not symbiosis genotypes suggests that fitness in the soil environment, rather than symbiosis with hosts, might be the key driver of Bradyrhizobium dominance.

Breakdown of coevolution between symbiotic bacteria Wolbachia and their filarial hosts

Wolbachia is an alpha-proteobacterial symbiont widely distributed in arthropods. Since the identification of Wolbachia in certain animal-parasitic nematodes (the Onchocercidae or filariae), the relationship between arthropod and nematode Wolbachia has attracted great interest. The obligate symbiosis in filariae, which renders infected species susceptible to antibiotic chemotherapy, was held to be distinct from the Wolbachia-arthropod relationship, typified by reproductive parasitism. While co-evolutionary signatures in Wolbachia-arthropod symbioses are generally weak, reflecting horizontal transmission events, strict co-evolution between filariae and Wolbachia has been reported previously. However, the absence of close outgroups for phylogenetic studies prevented the determination of which host group originally acquired Wolbachia. Here, we present the largest co-phylogenetic analysis of Wolbachia in filariae performed to date including: (i) a screening and an updated phylogeny of Wolbachia; (ii) a co-phylogenetic analysis; and (iii) a hypothesis on the acquisition of Wolbachia infection. First, our results show a general overestimation of Wolbachia occurrence and support the hypothesis of an ancestral absence of infection in the nematode phylum. The accuracy of supergroup J is also underlined. Second, although a global pattern of coevolution remains, the signal is derived predominantly from filarial clades associated with Wolbachia in supergroups C and J. In other filarial clades, harbouring Wolbachia supergroups D and F, horizontal acquisitions and secondary losses are common. Finally, our results suggest that supergroup C is the basal Wolbachia clade within the Ecdysozoa. This hypothesis on the origin of Wolbachia would change drastically our understanding of Wolbachia evolution.

HGTree: database of horizontally transferred genes determined by tree reconciliation

The HGTree database provides putative genome-wide horizontal gene transfer (HGT) information for 2472 completely sequenced prokaryotic genomes. This task is accomplished by reconstructing approximate maximum likelihood phylogenetic trees for each orthologous gene and corresponding 16S rRNA reference species sets and then reconciling the two trees under parsimony framework. The tree reconciliation method is generally considered to be a reliable way to detect HGT events but its practical use has remained limited because the method is computationally intensive and conceptually challenging. In this regard, HGTree (http://hgtree.snu.ac.kr) represents a useful addition to the biological community and enables quick and easy retrieval of information for HGT-acquired genes to better understand microbial taxonomy and evolution. The database is freely available and can be easily scaled and updated to keep pace with the rapid rise in genomic information.

Cophylogenetic Reconciliation with ILP

In this paper, we present an integer linear programming (ILP) approach, called CoRe-ILP, for finding an optimal time consistent cophylogenetic host-parasite reconciliation under the cophylogenetic event model with the events cospeciation, duplication, sorting, host switch, and failure to diverge. Instead of assuming event costs, a simplified model is used, maximizing primarily for cospeciations and secondarily minimizing host switching events. Duplications, sortings, and failure to diverge events are not explicitly scored. Different from existing event based reconciliation methods, CoRe-ILP can use (approximate) phylogenetic branch lengths for filtering possible ancestral host-parasite interactions. Experimentally, it is shown that CoRe-ILP can successfully use branch length information and performs well for biological and simulated data sets. The results of CoRe-ILP are compared with the results of the reconciliation tools Jane 4, Treemap 3b, NOTUNG 2.8 Beta, and Ranger-DTL. Algorithm CoRe-ILP is implemented using IBM ILOG CPLEX Optimizer 12.6 and is freely available from http://pacosy.informatik.uni-leipzig.de/core-ilp.

Efficient algorithms for reconciling gene trees and species networks via duplication and loss events

Reconciliation methods explain topology differences between a species tree and a gene tree by evolutionary events other than speciations. However, not all phylogenies are trees: hybridization can occur and create new species and this results into reticulate phylogenies. Here, we consider the problem of reconciling a gene tree with a species network via duplication and loss events. Two variants are proposed and solved with effcient algorithms: the first one finds the best tree in the network with which to reconcile the gene tree, and the second one finds the best reconciliation between the gene tree and the whole network.

Inferring gene duplications, transfers and losses can be done in a discrete framework

In the field of phylogenetics, the evolutionary history of a set of organisms is commonly depicted by a species tree-whose internal nodes represent speciation events-while the evolutionary history of a gene family is depicted by a gene tree-whose internal nodes can also represent macro-evolutionary events such as gene duplications and transfers. As speciation events are only part of the events shaping a gene history, the topology of a gene tree can show incongruences with that of the corresponding species tree. These incongruences can be used to infer the macro-evolutionary events undergone by the gene family. This is done by embedding the gene tree inside the species tree and hence providing a reconciliation of those trees. In the past decade, several parsimony-based methods have been developed to infer such reconciliations, accounting for gene duplications ([Formula: see text]), transfers ([Formula: see text]) and losses ([Formula: see text]). The main contribution of this paper is to formally prove an important assumption implicitly made by previous works on these reconciliations, namely that solving the (maximum) parsimony [Formula: see text] reconciliation problem in the discrete framework is equivalent to finding a most parsimonious [Formula: see text] scenario in the continuous framework. In the process, we also prove several intermediate results that are useful on their own and constitute a theoretical toolbox that will likely facilitate future theoretical contributions in the field.