RECONSTRUCTING THE HISTORY OF HOST-PARASITE ASSOCIATIONS USING GENERALISED PARSIMONY

Abstract — In reconstructing the history of host-parasite associations, it is necessary to consider several different processes, such as cospeciation and host switching, that may affect an association. A simple reconstruction method is to maximise the number of host-parasite cospeciations. However, maximum cospeciation reconstruction may require the postulation of a large number of other kinds of events, such as parasite extinction or exclusion from certain hosts. A more sophisticated method associates each kind of event with a cost or weight which is inversely related to the likelihood of that kind of event occurring. I present a method of the latter type that distinguishes between two different processes: host tracking, of which cospeciation is a special case, and host switching. Given a relative weight for these two types of events, it is possible to convert the host phytogeny into a cost matrix, allowing for host switching, and use generalised-parsimony algorithms to find minimum-cost reconstructions of the history of the host-parasite association. Different relative switch weights give different minimum-cost reconstructions; the optimal switch weight can be found by maximising the fit between the tracking events and the parasite phytogeny, controlling for the number of postulated switches. As an empirical application of the method, data on an association between pocket gophers and their parasitic chewing lice were re-examined. Although these data have been extensively analysed previously, the generalised parsimony approach throws new light on the history of the association.     

Biogeography explains cophylogenetic patterns in toucan chewing lice

Historically, comparisons of host and parasite phylogenies have concentrated on cospeciation. However, many of these comparisons have demonstrated that the phylogenies of hosts and parasites are seldom completely congruent, suggesting that phenomena other than cospeciation play an important role in the evolution of host-parasite assemblages. Other coevolutionary phenomena, such as host switching, parasite duplication (speciation on the host), sorting (extinction), and failure to speciate can also influence host-parasite assemblages. Using mitochondrial and nuclear protein-coding DNA sequences, I reconstructed the phylogeny of ectoparasitic toucan chewing lice in the Austrophilopterus cancellosus subspecies complex and compared this phylogeny with the phylogeny of the hosts, the Ramphastos toucans, to reconstruct the history of coevolutionary events in this host-parasite assemblage. Three salient findings emerged. First, reconstructions of host and louse phylogenies indicate that they do not branch in parallel, and their cophylogenetic history shows little or no significant cospeciation. Second, members of monophyletic Austrophilopterus toucan louse lineages are not necessarily restricted to monophyletic host lineages. Often, closely related lice are found on more distantly related but sympatric toucan hosts. Third, the geographic distribution of the hosts apparently plays a role in the speciation of these lice. These results suggest that for some louse lineages biogeography may be more important than host associations in structuring louse populations and species, particularly when host life history (e.g., hole nesting) or parasite life history (e.g., phoresis) might promote frequent host switching events between syntopic host species. These findings highlight the importance of integrating biogeographic information into cophylogenetic studies.     

A Bayesian framework for the analysis of cospeciation

Abstract.— Information on the history of cospeciation and host switching for a group of host and parasite species is contained in the DNA sequences sampled from each. Here, we develop a Bayesian framework for the analysis of cospeciation. We suggest a simple model of host switching by a parasite on a host phylogeny in which host switching events are assumed to occur at a constant rate over the entire evolutionary history of associated hosts and parasites. The posterior probability density of the parameters of the model of host switching are evaluated numerically using Markov chain Monte Carlo. In particular, the method generates the probability density of the number of host switches and of the host switching rate. Moreover, the method provides information on the probability that an event of host switching is associated with a particular pair of branches. A Bayesian approach has several advantages over other methods for the analysis of cospeciation. In particular, it does not assume that the host or parasite phylogenies are known without error; many alternative phylogenies are sampled in proportion to their probability of being correct.     

Evolutionary relationships, cospeciation, and host switching in avian malaria parasites

We used phylogenetic analyses of cytochrome b sequences of malaria parasites and their avian hosts to assess the coevolutionary relationships between host and parasite lineages. Many lineages of avian malaria parasites have broad host distributions, which tend to obscure cospeciation events. The hosts of a single parasite or of closely related parasites were nonetheless most frequently recovered from members of the same host taxonomic family, more so than expected by chance. However, global assessments of the relationship between parasite and host phylogenetic trees, using Component and ParaFit, failed to detect significant cospeciation. The event-based approach employed by TreeFitter revealed significant cospeciation and duplication with certain cost assignments for these events, but host switching was consistently more prominent in matching the parasite tree to the host tree. The absence of a global cospeciation signal despite conservative host distribution most likely reflects relatively frequent acquisition of new hosts by individual parasite lineages. Understanding these processes will require a more refined species concept for malaria parasites and more extensive sampling of parasite distributions across hosts. If parasites can disperse between allopatric host populations through alternative hosts, cospeciation may not have a strong influence on the architecture of host-parasite relationships. Rather, parasite speciation may happen more often in conjunction with the acquisition of new hosts followed by divergent selection between host lineages in sympatry. Detailed studies of the phylogeographic distributions of hosts and parasites are needed to characterize these events.     

STATISTICAL TESTS OF HOST-PARASITE COSPECIATION

A history of cospeciation (synchronous speciation) among ecologically associated, but otherwise distantly related, species is often revealed by a strong correspondence of their phylogenies. In this paper, we present several tests of cospeciation that use maximum-likelihood and Bayesian methods of phylogenetic estimation. The hypotheses tested include: (1) topological agreement of phylogenies for coevolving groups; (2) identical speciation times of associated species; and (3) identical evolutionary rates in genes of associated species. These tests are applied to examine a possible instance of host-parasite coevolution among pocket gophers and lice using mitochondrial COI DNA sequences. The observed differences between gopher and louse trees cannot be explained by sampling error and are consistent with a rate of host switching about one-third the host speciation rate. A subset of the gopher-louse data is consistent with a common history of evolution (i.e., the topologies and speciation times are identical). However, the relative rate of nucleotide substitution is two to four times higher in the lice than in the gophers.     

Monogeneans of West African cichlid fish: evolution and cophylogenetic interactions

The goals of this paper were to investigate phylogenetic and evolutionary patterns of cichlid fish from West Africa and their Cichlidogyrus and Scutogyrus monogenean parasites, to uncover the presence of host-parasite cospeciation and to assess the level of morphological adaptation in parasites. This required the following steps, each one representing specific objectives of this paper: (1) to build phylogenetic trees for Cichlidogyrus and Scutogyrus species based on ribosomal DNA sequences, (2) to investigate phylogenetic relationships within West African cichlid fish based on the analysis of mitochondrial cytochrome b DNA sequences, (3) to investigate host-parasite cophylogenetic history to gain clues on parasite speciation process, and (4) to investigate the link between the morphology of the attachment apparatus and parasite phylogeny. Phylogenetic analyses supported the monophyletic origin of the Cichlidogyrus/Scutogyrus group, and suggested that Cichlidogyrus is polyphyletic and that Scutogyrus is monophyletic. The phylogeny of Cichlidae supported the separation of mouthbrooders and substrate-brooders and is consistent with the hypothesis that the mouthbrooding behavior of Oreochromis and Sarotherodon evolved from substrate-brooding behavior. The mapping of morphological characters of the haptor onto the parasite phylogenetic tree suggests that the attachment organ has evolved from a very simple form to a more complex one. The cophylogenetic analyses indicated a significant fit between trees using distance-based tests, but no significant cospeciation signal using tree-based tests, suggesting the presence of parasite duplications and host switches on related host species. This shed some light on the diversification process of Cichlidogyrus species parasitizing West African cichlids.     

Clade-limited colonization in brood parasitic finches (Vidua spp.)

The African brood parasitic finches (Vidua spp.) are host specialists that mimic the songs and nestling mouth markings of their finch hosts (family Estrildidae). Although recent molecular analyses suggest rapid speciation associated with host switches in some members of this group, the association of different Vidua lineages with particular host genera suggests the possibility of cospeciation at higher levels in the host and parasite phylogenies. We compared a phylogeny of all Vidua species with a phylogeny of their estrildid finch hosts and compared divergence time estimates for the two groups. Basal divergences among extant members of the Vidulidae and among Vidua species are more recent than those among host genera and species, respectively, allowing a model of cospeciation to be rejected at most or all levels of the Vidua phylogeny. Nonetheless, some tests for cospeciation indicated significant congruence between host and parasite tree topologies. This result may be an artifact of clade-limited colonization. Host switches in parasitic finches have most often involved new hosts in the same or a closely related genus, an effect that increases the apparent congruence of host and parasites trees.     

When do parasites fail to speciate in response to host speciation?

Cospeciation generally increases the similarity between host and parasite phylogenies. Incongruence between host and parasite phylogenies has previously been explained in terms of host switching, sorting, and duplication events. Here, we describe an additional process, failure of the parasite to speciate in response to host speciation, that may be important in some host-parasite systems. Failure to speciate is likely to occur when gene flow among parasite populations is much higher than that of their hosts. We reconstructed trees from mitochondrial and nuclear DNA sequences for pigeons and doves (Aves: Columbiformes) and their feather lice in the genus Columbicola (Insecta: Phthiraptera). Although comparisons of the trees from each group revealed a significant amount of cospeciation, there was also a significant degree of incongruence. Cophylogenetic analyses generally indicated that host switching may be an important process in the history of this host-parasite association. Using terminal sister taxon comparisons, we also identified three apparent cases where the host has speciated but the associated parasite has not. In two of these cases of failure to speciate, these comparisons involve allopatric sister taxa of hosts whose lice also occur on hosts sympatric with both of the allopatric sisters. These additional hosts for generalist lice may promote gene flow with lice on the allopatric sister species. Relative rate comparisons for the mitochondrial cytochrome oxidase I gene indicate that molecular substitution occurs about 11 times faster in lice than in their avian hosts.     

Three-Dimensional Cost-Matrix Optimization and Maximum Cospeciation

In recent years, event-based approaches have been gaining ground in coevolutionary and biogeographical inference. Unlike pattern-based methods, event-based protocols deal directly with evolutionary events, such as dispersals and host switches. Three protocols have been proposed to date: (1) a coevolutionary method based on optimization of a standard two-dimensional cost matrix; (2) dispersal–vicariance analysis, based on optimization of a three-dimensional cost matrix; and (3) the maximum cospeciation method, thus far not considered a cost matrix method. I describe here general three-dimensional cost matrix optimization algorithms and how they can be applied to the maximum cospeciation problem. The new algorithms demonstrate that all existing event-based protocols, as well as possible future methods based on more complicated process models, can be incorporated into the three-dimensional cost matrix optimization framework.     

Coevolution between Lamellodiscus (Monogenea: Diplectanidae) and Sparidae (Teleostei): the study of a complex host-parasite system

Abstract.— Host-parasite coevolution was studied between Sparidae (Teleostei) fishes and their parasites of the genus Lamellodiscus (Monogenea, Diplectanidae) in the northwestern Mediterranean Sea. Molecular phylogenies were reconstructed for both groups. The phylogenetic tree of the Sparidae was obtained from previously published 16S mitochondrial DNA (mtDNA) sequences associated with new cytochrome-b mtDNA sequences via a "total evidence" procedure. The phylogeny of Lamellodiscus species was reconstructed from 18S rDNA sequences that we obtained. Host-parasite coevolution was studied through different methods: TreeFitter, TreeMap, and a new method, ParaFit. If the cost of a host switch is not assumed to be high for parasites, all methods agree on the absence of widespread cospeciation processes in this host-parasite system. Host-parasite associations were interpreted to be due more to ecological factors than to coevolutionary processes. Host specificity appeared not to be related to host-parasite cospeciation.     

Inferring processes of coevolutionary diversification in a community of Panamanian strangler figs and associated pollinating wasps*

The fig and pollinator wasp obligate mutualism is diverse (～750 described species), ecologically important, and ancient (～80 Ma). Once thought to be an example of strict one‐to‐one cospeciation, current thinking suggests genera of pollinator wasps codiversify with corresponding sections of figs, but the degree to which cospeciation or other processes contribute to the association at finer scales is unclear. Here, we use genome‐wide sequence data from a community of Panamanian strangler figs and associated wasp pollinators to estimate the relative contributions of four evolutionary processes generating cophylogenetic patterns in this mutualism: cospeciation, host switching, pollinator speciation, and pollinator extinction. Using a model‐based approach adapted from the study of gene family evolution, our results demonstrate the importance of host switching of pollinator wasps at this fine phylogenetic and regional scale. Although we estimate a modest amount of cospeciation, simulations reveal the number of putative cospeciation events to be consistent with what would be expected by chance. Additionally, model selection tests identify host switching as a critical parameter for explaining cophylogenetic patterns in this system. Our study demonstrates a promising approach through which the history of evolutionary association between interacting lineages can be rigorously modeled and tested in a probabilistic phylogenetic framework.     

Pseudocospeciation of the mycoparasite Cosmospora with their fungal hosts

Species of Cosmospora are parasites of other fungi (mycoparasites), including species belonging to the Xylariales. Based on prior taxonomic work, these fungi were determined to be highly host specific. We suspected that the association of Cosmospora and their hosts could not be a result of random chance, and tested the cospeciation of Cosmospora and the their hosts with contemporary methods (e.g., ParaFit, PACo, and Jane). The cophylogeny of Cosmospora and their hosts was found to be congruent, but only host‐parasite links in more recent evolutionary lineages of the host were determined as coevolutionary. Reconciliation reconstructions determined at least five host‐switch events early in the evolution of Cosmospora. Additionally, the rates of evolution between Cosmospora and their hosts were unequal. This pattern is more likely to be explained by pseudocospeciation (i.e., host switches followed by cospeciation), which also produces congruent cophylogenies.     

When can host shifts produce congruent host and parasite phylogenies? A simulation approach

Congruence between host and parasite phylogenies is often taken as evidence for cospeciation. However, 'pseudocospeciation', resulting from host-switches followed by parasite speciation, may also generate congruent trees. To investigate this process and the conditions favouring its appearance, we here simulated the adaptive radiation of a parasite onto a new range of hosts. A very high congruence between the host tree and the resulting parasite trees was obtained when parasites switched between closely related hosts. Setting a shorter time lag for speciation after switches between distantly related hosts further increased the degree of congruence. The shape of the host tree, however, had a strong impact, as no congruence could be obtained when starting with highly unbalanced host trees. The strong congruences obtained were erroneously interpreted as the result of cospeciations by commonly used phylogenetic software packages despite the fact that all speciations resulted from host-switches in our model. These results highlight the importance of estimating the age of nodes in host and parasite phylogenies when testing for cospeciation and also demonstrate that the results obtained with software packages simulating evolutionary events must be interpreted with caution.     

Host shift and cospeciation rate estimation from co‐phylogenies

Host shifts can cause novel infectious diseases, and is a key process in diversification. Disentangling the effects of host shift vs. those of cospeciation is non‐trivial as both can result in phylogenic congruence. We develop a new framework based on network analysis and Approximate Bayesian Computation to quantify host shift and cospeciation rates in host‐parasite systems. Our method enables estimation of the expected time to the next host shift or cospeciation event. We then apply it to avian haemosporidian parasite systems and to the pocket gophers‐chewing lice system, and demonstrate that both host shift and cospeciation can be reliably estimated by our method. We confirm that host shifts have shaped the evolutionary history of avian haemosporidian parasites and have played a minor role in the gopher–chewing lice system. Our method is promising for predicting the rate of potential host shifts and thus the emergence of novel infectious diseases.     

Evolutionary relationships between digeneans of the family Brachycladiidae Odhner, 1905 and their marine mammal hosts: A cophylogenetic study

Cophylogenetic studies examine the congruence between host and parasite phylogenies. There are few studies that quantify the relative contribution of coevolutionary events, i.e. duplication, loss, failure-to-diverge, host-switching and spreading in trophically-transmitted parasites at the marine realm. We addressed this issue in the Brachycladiidae, a cosmopolitan digenean family specific to marine mammals. We used, for the first time, distance-based and event-based methods to explicitly test the coevolutionary events that have shaped the current brachycladiid-marine mammal associations. Parasite phylogeny was constructed using mtDNA ND3 sequences of nine brachycladiid species, and host phylogeny using cytochrome b sequences of 104 mammalian species. A total of 50 host-parasite links were identified. Distance-based methods supported the hypothesis of a global non-random association of host and parasite phylogenies. Significant individual links (i.e., 24 out of 50) were those related to Campula oblonga, Nasitrema delphini, N. globicephalae and Brachycladium atlanticum and their associated taxa from the Delphinoidea. Regarding event-based methods, we explored 54 schemes using different combinations of costs for each potential coevolutionary event. Three coevolutionary scenarios were identified across all schemes and in all cases the number of loss events (87–156) was the most numerous, followed by failure-to-diverge (40), duplication (3–6), host-switching (0–3) and cospeciation (0–2). We developed a framework to interpret the evolution of this host-parasite system and confirmed that failure-to-diverge and colonization with or without subsequent diversification could have been decisive in the establishment of the associations between brachycladiids and marine mammals.     

A priori and a posteriori methods in comparative evolutionary studies of host–parasite associations

Brooks parsimony analysis (BPA) and reconciliation methods in studies of host–parasite associations differ fundamentally, despite using the same null hypothesis. Reconciliation methods may eliminate or modify input data to maximize fit of single parasite clades to a null hypothesis of cospeciation, by invoking different a priori assumptions, including a known host phylogeny. By examining the degree of phylogenetic congruence among multiple parasite clades, using hosts as analogs of taxa but not presuming a host phylogeny or any degree of cospeciation a priori, BPA modifies the null hypothesis of cospeciation if necessary to maintain the integrity of the input data. Two exemplars illustrate critical empirical differences between reconciliation methods and BPA: (1) reconciliation methods rather than BPA may select the incorrect general host cladogram for a set of data from different clades of parasites, (2) BPA rather than reconciliation methods provides the most parsimonious interpretation of all available data, and (3) secondary BPA, proposed in 1990, when applied to data sets in which host‐switching produces hosts with reticulate histories, provides the most parsimonious and biologically realistic interpretations of general host cladograms. The extent to which these general host cladograms, based on cospeciation among different parasite clades inhabiting the same hosts, correspond to host phylogeny can be tested, a posteriori, by comparison with a host phylogeny generated from nonparasite data. These observations lead to the conclusion that BPA and reconciliation methods are designed to implement different research programs based on different epistemologies. BPA is an a posteriori method that is designed to assess the host context of parasite speciation events, whereas reconciliation methods are a priori methods that are designed to fit parasite phylogenies to a host phylogeny. Host‐switching events are essential for explaining complex histories of host–parasite associations. BPA assumes coevolutionary complexity (historical contingency), relying on parsimony as an a posteriori explanatory tool to summarize complex results, whereas reconciliation methods, which embody formalized assumptions of maximum cospeciation, are based on a priori conceptual parsimony. Modifications of basic reconciliation methods, embodied in TreeMap 1.0 and TreeMap 2.02, represent the addition of weighting schemes in which the researcher specifies allowed departures from cospeciation a priori, with the result that TreeMap results more closely agree with BPA results than do reconciled tree analysis results.     

Evolutionary history of nematodes associated with sweat bees

Organisms that live in close association with other organisms make up a large part of the world’s diversity. One driver of this diversity is the evolution of host-species specificity, which can occur via reproductive isolation following a host-switch or, given the correct circumstances, via cospeciation. In this study, we explored the diversity and evolutionary history of Acrostichus nematodes that are associated with halictid bees in North America. First, we conducted surveys of bees in Virginia, and found six halictid species that host Acrostichus. To test the hypothesis of cospeciation, we constructed phylogenetic hypotheses of Acrostichus based on three genes. We found Acrostichus puri and Acrostichus halicti to be species complexes comprising cryptic, host-specific species. Although several nodes in the host and symbiont phylogenies were congruent and tests for cospeciation were significant, the host’s biogeography, the apparent patchiness of the association across the host’s phylogeny, and the amount of evolution in the nematode sequence suggested a mixture of cospeciation, host switching, and extinction events instead of strict cospeciation. Cospeciation can explain the relationships between Ac. puri and its augochlorine hosts, but colonization of Halictus hosts is more likely than cospeciation. The nematodes are vertically transmitted, but sexual transmission is also likely. Both of these transmission modes may explain host-species specificity and congruent bee and nematode phylogenies. Additionally, all halictid hosts come from eusocial or socially polymorphic lineages, suggesting that sociality may be a factor in the suitability of hosts for Acrostichus.     

Plant-insect interactions: double-dating associated insect and plant lineages reveals asynchronous radiations

An increasing number of plant-insect studies using phylogenetic analysis suggest that cospeciation events are rare in plant-insect systems. Instead, nonrandom patterns of phylogenetic congruence are produced by phylogenetically conserved host switching (to related plants) or tracking of particular resources or traits (e.g., chemical). The dominance of host switching in many phytophagous insect groups may make the detection of genuine cospeciation events difficult. One important test of putative cospeciation events is to verify whether reciprocal speciation is temporally plausible. We explored techniques for double-dating of both plant and insect phylogenies. We use dated molecular phylogenies of a psyllid (Hemiptera)-Genisteae (Fabaceae) system, a predominantly monophagous insect-plant association widespread on the Atlantic Macaronesian islands. Phylogenetic reconciliation analysis suggests high levels of parallel cladogenesis between legumes and psyllids. However, dating using molecular clocks calibrated on known geological ages of the Macaronesian islands revealed that the legume and psyllid radiations were not contemporaneous but sequential. Whereas the main plant radiation occurred some 8 million years ago, the insect radiation occurred about 3 million years ago. We estimated that >60% of the psyllid speciation has resulted from host switching between related hosts. The only evidence for true cospeciation is in the much more recent and localized radiation of genistoid legumes in the Canary Islands, where the psyllid and legume radiations have been partially contemporaneous. The identification of specific cospeciation events over this time period, however, is hindered by the phylogenetic uncertainty in both legume and psyllid phylogenies due to the apparent rapidity of the species radiations.     

Cophylogeny: insights from fish-parasite systems

Host-parasite cophylogeny is a topic that has grasped the attention of scientists since the end of the 19th century, but the development of dedicated analytical methods only arose in the last 30 years. Research on host-parasite systems and on the development of more and more sophisticated numerical methods to estimate the degree of cospeciation has thus progressed, permitting the elaboration of evolutionary scenarios. The main outcome of these studies is that the expected clear pattern of cospeciation between many hosts and parasites is often obscure. In practice, much attention has been devoted to few host-parasite systems. Particularly, aquatic host-parasite associations have not been so extensively studied, and, after briefly reviewing the main analytical methods, this paper focuses on host-monogenean systems, because this kind of interaction is expected to be an ideal model for cophylogeny studies. But is it? And what does it tell us about the evolutionary and ecological forces driving cospeciation in the open sea? Biogeography should also be considered when possible, and it has been useful for explaining some patterns of cospeciation. It should thus be more deeply exploited in the future. We need new methods and new biological models that better, if not fully, depict patterns and thereby permit deeper understanding of processes within cophylogenetic patterns.     

Linking coevolutionary history to ecological process: doves and lice

Abstract Many host-specific parasites are restricted to a limited range of host species by ecological barriers that impede dispersal and successful establishment. In some cases, microevolutionary differentiation is apparent on top of host specificity, as evidenced by significant parasite population genetic structure among host populations. Ecological barriers responsible for specificity and genetic structure can, in principle, reinforce macroevolutionary processes that generate congruent host-parasite phylogenies. However, few studies have explored both the micro- and macroevolutionary ramifications of close association in a single host-parasite system. Here we compare the macroevolutionary histories of two genera of feather lice (Phthiraptera: Ischnocera) that both parasitize New World pigeons and doves (Aves: Columbiformes). Earlier work has shown that dove body lice (genus Physconelloides ) are more host specific and have greater population genetic structure than dove wing lice ( Columbicola ). We reconstructed phylogenies for representatives of the two genera of lice and their hosts, using nuclear and mitochondrial DNA sequences. The phylogenies were well resolved and generally well supported. We compared the phylogenies of body lice and wing lice to the host phylogeny using reconciliation analyses. We found that dove body lice show strong evidence of cospeciation whereas dove wing lice do not. Although the ecology of body and wing lice is very similar, differences in their dispersal ability may underlie these joint differences in host specificity, population genetic structure, and coevolutionary history.