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.     

Comparing host and parasite phylogenies: gyrodactylus flatworms jumping from goby to goby

The combination of exceptionally high species diversity, high host specificity, and a complex reproduction system raises many questions about the underlying mechanisms triggering speciation in the flatworm genus Gyrodactylus. The coevolutionary history with their goby hosts was investigated using both topology- and distance-based approaches; phylogenies were constructed of the V4 region of the 18S rRNA and the complete ITS rDNA region for the parasites, and 12S and 16S mtDNA fragments for the hosts. The overall fit between both trees was significant according to the topology-based programs (TreeMap 1.0, 2.0 beta and TreeFitter), but not according to the timed analysis in TreeMap 2.0 beta and the distance-based method (ParaFit). An absolute timing of speciation events in host and parasite ruled out the possibility of synchronous speciation for the gill parasites, favouring the distance-based result. Based on this information together with the biological background of host and parasite, the following TreeMap solution was selected. The group of gill parasites evolved from a host switch from G. arcuatus, parasitizing the three-spined stickleback onto the gobies, followed by several host-switching events among the respective goby hosts. The timing of these events is estimated to date back to the Late Pleistocene, suggesting a role for refugia-mediated mixing of parasite species. In contrast, it is suggested that co-speciation in the fin-parasites resulted in several host-associated species complexes. This illustrates that phylogenetically conserved host-switching mimics the phylogenetic signature of co-speciation, confounding topology-based programs.     

Ending a decade of deception: a valiant failure, a not-so-valiant failure, and a success story

Prior studies involving two methods, Brooks Parsimony Analysis (BPA) and TreeMap, have found BPA to be the more reliable method. Recent criticisms leveled at these studies argue that the tests were unfairly created and biased in favor of BPA. The authors of a recent critique offered new exemplars to demonstrate flaws in BPA, plus a simple fix to correct the flaws found in TreeMap. A re‐evaluation of their exemplars clearly shows that the authors' calculations are incorrect, their understanding of the methods is lacking, and that their simple fix does not work. Additional analyses using TreeMap 2.02 are run to show that TreeMap 2.02, like TreeMap 1.0, cannot adequately deal with widespread parasites, contrary to the claims of its supporters. Furthermore, the exemplars corroborate previous findings that BPA, when calculated correctly, is more reliable than TreeMap1.0 and TreeMap 2.02 and therefore the method of choice in coevolutionary and biogeographic studies.     

Brooks Parsimony Analysis: a valiant failure

A recent comparison of two methods for examining correlated host and parasite phylogenies, namely TreeMap 1.0 and Brooks Parsimony Analysis concluded that the latter method performed better and is to be preferred. Reevaluation of the examples contrived for that study demonstrates that the two methods were only compared on one kind of problem (widespread parasite) for which there is an easy fix for TreeMap 1.0. Other kinds of problems like host‐switching among sister taxa or host‐switching over great distances across a host tree befuddle BPA even as they are readily resolved parsimoniously by TreeMap 1.0. These difficulties, compounded with inaccurate counting of ad hoc hypotheses required by its solutions render BPA unsuitable for comparison of host and parasite phylogenies.     

Extending phylogenetic studies of coevolution: secondary Brooks parsimony analysis, parasites, and the Great Apes

Dowling recently compared the empirical properties of Brooks parsimony analysis (BPA) and the leading method for studying phylogenetic aspects of coevolution, reconciled tree analysis (using the computer program TreeMap), based on a series of simulations. Like the majority of authors who have compared BPA with other methods, however, Dowling considered only the form of BPA proposed in 1981 and did not take into account various modifications of the method proposed from 1986 to 2002. This leaves some doubt as to the robustness of his assessments of both the superiority of BPA and its shortcomings. We provide a précis of the principles of contemporary BPA, including ways to implement it algorithmically, using either Wagner algorithm‐based or Hennigian argumentation‐based approaches, followed by an empirical example. Our study supports Dowling's fundamental conclusions about the superiority of primary BPA relative to TreeMap. However, his conclusions about the shortcomings of BPA due to inclusive ORing (i.e., the production of ghost taxa) are incorrect, as secondary BPA eliminates inclusive ORing from the method. Secondary BPA provides a more complete account of the evolutionary associations between the parasite groups and their hosts than does primary BPA, without sacrificing any indirectly generated information about host phylogeny. Secondary BPA of two groups of nematodes inhabiting Great Apes shows that TreeMap analysis underestimated the amount of cospeciation in the evolution of the nematode genus Enterobius.     

PARALLEL PHYLOGENIES: RECONSTRUCTING THE HISTORY OF HOST-PARASITE ASSEMBLAGES

Abstract— A method for reconstructing the history of a host-parasite assemblage is described. This method has the advantage of making explicit the relationship between the host and parasite trees, and it allows a visually intuitive representation of that history. It also enables host switches to be incorporated as an explanation of the observed pattern of host-parasite associations, without the spurious overestimates of the number of host switches that can be obtained using Brooks parsimony analysis (BPA). Reconstructions that maximize the number of cospeciation events have the greatest explanatory power and are hence preferred over reconstructions with fewer cospeciation events. A heuristic algorithm to find a single maximal reconstruction, and an exact algorithm to find all such reconstructions are presented. Two empirical applications of the method are given.     

寄生虫及其宿主协同进化的研究进展

本文对寄生虫及其宿主协同进化的研究进行了回顾,将其发展分为三个阶段：1．寄生虫-宿主协同进化的初步认识;2．协同进化模式及其内在机制的探索;3．协同进化机制研究方法的发展。目前的研究主要集中于协同进化生物学意义的进一步深入探讨。同时,对协同进化的有关概念、方法和本学科的发展进行了简单阐述和讨论。     

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.     

SELECTION AND STRAIN SPECIFICITY OF COMPATIBILITY BETWEEN SNAIL INTERMEDIATE HOSTS AND THEIR PARASITIC SCHISTOSOMES

Many theoretical models of host-parasite coevolution assume that variation in host resistance to parasite infection is, at least partially, genetically determined and specific to the strain of infecting parasite. However, very few experimental studies have been conducted to test this assumption in animal-parasite systems. Biomphalaria glabrata snails serve as the intermediate hosts of Schistosoma mansoni. Although some snails are resistant to infection, there is no evidence of fixation of resistance in field populations. Two possible explanations for this are high fitness costs associated with resistance and a dynamic coevolution between parasite and host, perhaps involving matching alleles or gene-for-gene interactions. Two strains of B. glabrata were artificially selected for either resistance or susceptibility to each of two strains of S. mansoni parasite for three generations. Third-generation snails were then were exposed to either the parasite strain to which they had been selected or to a different parasite strain. In both host strains, resistance and susceptibility (compatibility) were found to be heritable. Moreover, compatibility to one parasite strain was not associated with compatibility to another strain, implying no genetic trade-off. Our results are discussed in terms of potential mechanisms of resistance in this host-parasite system and their implications to general coevolutionary theory.     

A protocol for temporal calibration of general area cladograms

Aim To describe a protocol for incorporating a temporal dimension into historical biogeographical analysis, while maintaining the essential independence of all datasets, involving the generation of general area cladograms. Location Global. Methods General area cladograms (GACs) are a reconstruction of the evolutionary history of a set of areas and unrelated clades within those areas. Nodes on a GAC correspond to speciation events in a group of taxa; general nodes are those at which multiple unrelated clades speciate. We undertake temporal calibration of GACs using molecular clock estimates of splitting events between extant taxa as well as first appearance data from the fossil record. We present two examples based on re‐analysis of previously published data: first, a temporally calibrated GAC generated from secondary Brooks parsimony analysis (BPA) of six extant bird clades from the south‐west of North America using molecular clock estimates of divergence times; and second, an analysis of African Neogene mammals based on a phylogenetic analysis for comparing trees (PACT) analysis. Results A hypothetical example demonstrates how temporal calibration reveals potentially critical information about the timing of both unique and general events, while also illustrating instances of incongruence between dates generated from molecular clock estimates and fossils. For the African Neogene mammal dataset, our analysis reveals that most mammal clades underwent geodispersal associated with the Neogene climatic optimum (c. 16 Ma) and vicariant speciation in central Africa correlated with increased aridity and cooler temperatures around 2.5 Ma. Main conclusions Temporally calibrated GACs are valuable tools for assessing whether coordinated patterns of speciation are associated with large‐scale climatic or tectonic phenomena.     

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.     

Evolution of spatially structured host–parasite interactions

Spatial structure has dramatic effects on the demography and the evolution of species. A large variety of theoretical models have attempted to understand how local dispersal may shape the coevolution of interacting species such as host–parasite interactions. The lack of a unifying framework is a serious impediment for anyone willing to understand current theory. Here, we review previous theoretical studies in the light of a single epidemiological model that allows us to explore the effects of both host and parasite migration rates on the evolution and coevolution of various life‐history traits. We discuss the impact of local dispersal on parasite virulence, various host defence strategies and local adaptation. Our analysis shows that evolutionary and coevolutionary outcomes crucially depend on the details of the host–parasite life cycle and on which life‐history trait is involved in the interaction. We also discuss experimental studies that support the effects of spatial structure on the evolution of host–parasite interactions. This review highlights major similarities between some theoretical results, but it also reveals an important gap between evolutionary and coevolutionary models. We discuss possible ways to bridge this gap within a more unified framework that would reconcile spatial epidemiology, evolution and coevolution.     

Phylogenetic analysis of European species of Proteocephalus (Cestoda: Proteocephalidea): compatibility of molecular and morphological data, and parasite-host coevolution.

The phylogeny of European species of the tapeworm genus Proteocephalus was studied, based on partial 18S rDNA and morphological data. The group was found to be monophyletic. The analysis showed a low informative value of available morphological characters in comparison with molecular data. The morphological matrix resulted in a poorly resolved tree which is, however, compatible with the topology (Proteocephalus osculatus (Proteocephalus torulosus (Proteocephalus macrocephalus, Proteocephalus filicollis) (Proteocephalus tetrastomus, Proteocephalus percae, Proteocephalus longicollis))) based on the 18S rDNA data. A comparison performed by the program TreeMap showed a lack of significant congruency between parasite and host phylogenies. Therefore, the distribution of species in their hosts appears to be independent of the phylogeny and it is likely to be a result of host-switching, rather than co-speciation events.     

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.     

Marine Area Relationships from Twenty Sponge Phylogenies. A Comparison of Methods and Coding Strategies

Published phylogenies of 20 marine sponge groups are used to build general area cladograms of marine areas of endemism under three different methods for which algorithms adapted for personal computers are available, viz. COMPONENT, BPA and TAS, and two different coding strategies, Assumption 0 (A0) and "no assumption" (NA). The latter is a recently proposed procedure for handling the distributions of widespread taxa by treating these as separate areas of endemism, rather than as suites of smaller constituent areas. The 20 phylogenies contained large numbers of problem data which prevented an exhaustive search for all possible equally "best" general area cladograms. The Nelson consensus trees and their equivalents in parsimony analysis for all six attempts (viz. three different methodologies under two different coding strategies) were compared using their fit with the 20 sponge phylogenies as a measure of quality. Fit was determined using the number of "cospeciations" between a general area cladogram and a taxon area cladogram computed with TreeMap 1.0. No single method or coding strategy yielded a clearly better fit, each cladogram fitting variously better or worse with various phylogenies. In general, fit with NA coding was higher than with A0 coding, but random tree tests failed to generate statistically significant support for the conclusion that NA coding improves fit. Assuming that available sponge phylogenies are representative of marine benthic groups, software and hardware limitations are serious obstacles to a successful development of marine general area cladograms under any method or coding strategy.     

Have chondracanthid copepods co-speciated with their teleost hosts?

Chondracanthid copepods parasitise many teleost species and have a mobile larval stage. It has been suggested that copepod parasites, with free-living infective stages that infect hosts by attaching to their external surfaces, will have co-evolved with their hosts. We examined copepods from the genus Chondracanthus and their teleost hosts for evidence of a close co-evolutionary association by comparing host and parasite phylogenies using TreeMap analysis. In general, significant co-speciation was observed and instances of host switching were rare. The prevalence of intra-host speciation events was high relative to other such studies and may relate to the large geographical distances over which hosts are spread.     

Coevolutionary interactions in a host–parasite system

Interactions between parasitic cuckoos and their hosts represent a classic example of coevolution, where adaptations in the parasite to exploit the host select for defences, which in turn select for new parasite adaptations. Current interactions between the two parties may be at an evolutionary equilibrium or, alternatively, a coevolutionary arms race may be taking place. By taking into account the effect of gene flow in 15 European magpie ( Pica pica ) populations, we studied the coevolutionary interactions with its brood parasite, the great spotted cuckoo ( Clamator glandarius ). Our results suggest that, in Europe, magpies and cuckoos are engaged in an ongoing coevolutionary process because, despite controlling for the large amounts of gene flow among different magpie populations, we still found a positive relationship between host defence (i.e. foreign egg recognition and rejection) and parasite selection pressure.     

Weighted ancestral area analysis and a solution of the redundant distribution problem

A new cladistic method for the estimation of ancestral areas is based on reversible parsimony in combination with a weighting scheme that weights steps in positionally plesiomorphic branches more highly than steps in positionally apomorphic branches. By applying this method to cladograms of human mitochondrial DNA, the method is superior to previously proposed algorithms. The method is also an appropriate tool for the solution of the redundant distribution problem in area cladograms. Under the assumption of allopatric speciation, redundant distributions, i.e., sympatry of sister groups, show that dispersal has occurred; thus, the ancestral area of at least one sister group was smaller than the combined distribution of its descendants. With the weighted ancestral area analysis, the ancestral areas can be confined and at least some dispersal events can be distinguished from possible vicariance events. As applied to a cladogram of the Polypteridae, weighted ancestral area analysis is superior to Brooks parsimony analysis (assumption 0) and component analysis under assumptions 1 and 2 (Nelson and Platnick, 1981, Systematics and biogeography: Cladistics and vicariance. Columbia Univ. Press, New York.) in resolving redundancies. The results of the weighted ancestral area analysis may differ from the results of dispersal-vicariance analysis, because the rules of dispersal-vicariance analysis indirectly favor the questionable assumption that the ancestral species occupied only one unit area.     

Seabird and louse coevolution: complex histories revealed by 12S rRNA sequences and reconciliation analyses

We investigated the coevolutionary history of seabirds (orders Procellariiformes and Sphenisciformes) and their lice (order Phthiraptera). Independent trees were produced for the seabirds (tree derived from 12S ribosomal RNA, isoenzyme, and behavioral data) and their lice (trees derived from 12S rRNA data). Brook's parsimony analysis (BPA) supported a general history of cospeciation (consistency index = 0.84, retention index = 0.81). We inferred that the homoplasy in the BPA was caused by one intrahost speciation, one potential host-switching, and eight or nine sorting events. Using reconciliation analysis, we quantified the cost of fitting the louse tree onto the seabird tree. The reconciled trees postulated one host-switching, nine cospeciation, three or four intrahost speciation, and 11 to 14 sorting events. The number of cospeciation events was significantly more than would be expected from chance alone (P < 0.01). The sequence data were used to test for rate heterogeneity for both seabirds and lice. Neither data set displayed significant rate heterogeneity. An examination of the codivergent nodes revealed that seabirds and lice have cospeciated synchronously and that lice have evolved at approximately 5.5 times the rate of seabirds. The degree of sequence divergence supported some of the postulated intrahost speciation events (e.g., Halipeurus predated the evolution of their present hosts). The sequence data also supported some of the postulated host-switching events. These results demonstrate the value of sequence data and reconciliation analyses in unraveling complex histories between hosts and their parasites.     

Genetic correlations and the coevolutionary dynamics of three-species systems

The majority of species interact with at least several others. We develop simple genetic models of coevolution between three species where interactions are mediated by quantitative traits. We assume that one of the species has two quantitative traits, each of which governs its interaction with one of the other two species. We use this model to explore how genetic correlations between the two traits in the multivariate species shape the evolutionary dynamics and outcomes of three species interactions. Our results suggest that genetic correlations are most important when at least one of the interactions is between a predator and prey or parasite and host. In these cases, genetic correlations between traits lead to a wide variety of novel coevolutionary outcomes and dynamics. In particular, genetic correlations can affect the existence and stability of coevolutionary equilibrium points, and they can lead to recurrent or permanent maladaptation. When the three species interact only as competitors or mutualists, however, genetic correlations have no effect on the outcome of coevolution. In all cases, our results reveal the surprising conclusion that both positive and negative genetic correlations between traits have qualitatively identical effects on coevolutionary dynamics.