Phylogenetic approaches in coevolution and biogeography

I review phylogenetic approaches to problems in coevolution and biogeography, illustrating with case studies. In coevolution, genealogical trees are essential in differentiating between ancient and recent associations, in identifying cospeciation events, and in studying host-switching patterns. Cospeciating associations are of particular interest because they allow powerful tests of molecular clocks and accurate comparison of evolutionary rates across groups of organisms. In biogeography, phylogenies can help reconstruct the distribution history of individual groups and identify past geological events that have affected the evolution of entire communities. Parsimony analysis in coevolution and biogeography should be based on identification of different types of events, each of which is associated with a specific cost. Similar event-based methods are applicable to coevolutionary and biogeographic inference, as well as in the mapping of gene trees onto organism trees. The discussed examples span a variety of organisms and spatiotemporal scales: primate pin worms, HIV, pocket gophers and their lice, aphids and their bacterial symbionts, gall wasps and their host plants, the root of the tree of life, the historical biogeography of the Holarctic, and the geographical origin of our own species.     

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.     

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.     

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.     

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.     

Analytical approaches to measuring cospeciation of host and parasites: through a glass, darkly

Studies of cophylogenetic associations between hosts and parasites have become increasingly common. Historically, congruence between host and parasite phylogenies has been seen as evidence for cospeciation. Analyses of such coevolutionary relationships, however, are made extremely difficult by the complex interplay of cospeciation, host switching, sorting (extinction), duplication (intrahost speciation) and inertia (lack of parasite speciation) events, all of which may produce incongruence between host and parasite phylogenies. Here we review several methods of analysing cospeciation. We illustrate these methods with an example from a Procellariiformes (seabird) and chewing louse (Halipeurus) association.     

Unravelling evolutionary relationships between epifoliar Meliolaceae and angiosperms

Epifoliar fungi are a group of poorly studied fungal symbionts that coinhabit the surface of living plants. Meliolaceae is the largest group of epifoliar fungi and has been considered as obligate parasites. We investigated the taxonomy of Meliolaceae and the coevolutionary events with their host plants using time-calibrated cophylogeny based on large subunit, small subunit, and internal transcribed spacer (ITS) sequence data obtained from 17 different fungal taxa and rbcL, ITS, and trnH-psbA sequence data from their corresponding hosts. Nine new fungal species are introduced in this paper and Appendiculella is synonymized under Asteridiella. The dominant coevolutionary events during the Cretaceous and Cenozoic are cospeciation and host shift, respectively. We hypothesize that the evolutionary history of epifoliar fungi can be divided into three major periods: origins of families, formations of genera, and diversification of species. The rise of angiosperms prompted the evolution of modern epifoliar fungi and the diversification of orders of Angiospermae fostered the formation of epifoliar fungal genera. Phylogenetically, epifoliar fungal genera can be delimited according to their coevolutionary patterns and divergent periods.     

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.     

Structure optimization in a three-dimensional off-lattice protein model

We studied a three-dimensional off-lattice AB model with two species of monomers, hydrophobic (A) and hydrophilic (B), and present two optimization algorithms: face-centered-cubic (FCC)-lattice pruned-enriched-Rosenbluth method (PERM) and subsequent conjugate gradient (PERM++) minimization and heuristic conjugate gradient (HCG) simulation based on "off-trap" strategy. In PERM++, we apply the PERM to the FCC-lattice to produce the initial conformation, and conjugate gradient minimization is then used to reach the minimum energy state. Both algorithms have been tested in the three-dimensional AB model for all sequences with lengths 13 < or = n < or = 55. The numerical results show that the proposed methods are very promising for finding the ground states of proteins. In several cases, we renew the putative ground states energy values.     

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.     

Histories of host shifts and cospeciation among free‐living parasitoids of Rhagoletis flies

Host shifts by specialist insects can lead to reproductive isolation between insect populations that use different hosts, promoting diversification. When both a phytophagous insect and its ancestrally associated parasitoid shift to the same novel host plant, they may cospeciate. However, because adult parasitoids are free living, they can also colonize novel host insects and diversify independent of their ancestral host insect. Although shifts of parasitoids to new insect hosts have been documented in ecological time, the long‐term importance of such shifts to parasitoid diversity has not been evaluated. We used a genus of flies with a history of speciation via host shifting (Rhagoletis [Diptera: Tephritidae]) and three associated hymenopteran parasitoid genera (Diachasma, Coptera and Utetes) to examine cophylogenetic relationships between parasitoids and their host insects. We inferred phylogenies of Rhagoletis, Diachasma, Coptera and Utetes and used distance‐based cophylogenetic methods (ParaFit and PACo) to assess congruence between fly and parasitoid trees. We used an event‐based method with a free‐living parasitoid cost model to reconstruct cophylogenetic histories of each parasitoid genus and Rhagoletis. We found that the current species diversity and host–parasitoid associations between the Rhagoletis flies and parasitoids are the primary result of ancient cospeciation events. Parasitoid shifts to ancestrally unrelated hosts primarily occur near the branch tips, suggesting that host shifts contribute to recent parasitoid species diversity but that these lineages may not persist over longer time periods. Our analyses also stress the importance of biologically informed cost models when investigating the coevolutionary histories of hosts and free‐living parasitoids.     

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.     

Estimators of repeatability

Summary Analysis of variance and principal components methods have been suggested for estimating repeatability. In this study, six estimation procedures are compared: ANOVA, principal components based on the sample covariance matrix and also on the sample correlation matrix, a related multivariate method (structural analysis) based on the sample covariance matrix and also on the sample correlation matrix, and maximum likelihood estimation. A simulation study indicates that when the standard linear model assumptions are met, the estimators are quite similar except when the repeatability is small. Overall, maximum likelihood appears the preferred method. If the assumption of equal variance is relaxed, the methods based on the sample correlation matrix perform better although others are surprisingly robust. The structural analysis method (with sample correlation matrix) appears to be best.Paper number 776 from the Department of Meat and Animal Science, University of Wisconsin-Madison.     

Cophylogenetic analysis of the relationship between anemonefish Amphiprion (Perciformes: Pomacentridae) and their symbiotic host anemones (Anthozoa: Actiniaria)

ABSTRACT

The association between anemonefish and anemone is a classical example of mutualism in coral reefs. Although mutualism is probably the key innovation that triggered the adaptive radiation of anemonefish into a wide range of habitats, the coevolutionary history between the groups has not been thoroughly tested in a phylogenetic framework. We examined the evolutionary history of the association via distance-based (Parafit and PACo) and event-based methods (Core-PA, Jane). Mitochondrial DNA sequences (COI mtDNA, Cytb, 16S rDNA and 12S rDNA) were used to reconstruct the phylogenies of tropical Amphiprion species and their host anemones by using maximum likelihood with best-fit models selected. Neither distance-based analyses nor event-based analyses revealed global significant congruence between the phylogenies of the hosts and the symbionts, and thus no evidence for coevolution between anemone-anemonefish. However, at the individual pair level, the fish showed some dependence on anemone hosts. Even though living in close association and benefiting from each other, the change of genetic composition of one species (anemonefish) does not always evolve in response to changes in the other (anemones). These findings expand our understanding of the pattern and the role of evolutionary events to allow a better prediction of the future of the anemonefish-anemone relationship.     

A Class-Information-Based Penalized Matrix Decomposition for Identifying Plants Core Genes Responding to Abiotic Stresses

In terms of making genes expression data more interpretable and comprehensible, there exists a significant superiority on sparse methods. Many sparse methods, such as penalized matrix decomposition (PMD) and sparse principal component analysis (SPCA), have been applied to extract plants core genes. Supervised algorithms, especially the support vector machine-recursive feature elimination (SVM-RFE) method, always have good performance in gene selection. In this paper, we draw into class information via the total scatter matrix and put forward a class-information-based penalized matrix decomposition (CIPMD) method to improve the gene identification performance of PMD-based method. Firstly, the total scatter matrix is obtained based on different samples of the gene expression data. Secondly, a new data matrix is constructed by decomposing the total scatter matrix. Thirdly, the new data matrix is decomposed by PMD to obtain the sparse eigensamples. Finally, the core genes are identified according to the nonzero entries in eigensamples. The results on simulation data show that CIPMD method can reach higher identification accuracies than the conventional gene identification methods. Moreover, the results on real gene expression data demonstrate that CIPMD method can identify more core genes closely related to the abiotic stresses than the other methods.     

Diversification and coevolution in brood pollination mutualisms: Windows into the role of biotic interactions in generating biological diversity

Brood pollination mutualisms—interactions in which specialized insects are both the pollinators (as adults) and seed predators (as larvae) of their host plants—have been influential study systems for coevolutionary biology. These mutualisms include those between figs and fig wasps, yuccas and yucca moths, leafflowers and leafflower moths, globeflowers and globeflower flies, Silene plants and Hadena and Perizoma moths, saxifrages and Greya moths, and senita cacti and senita moths. The high reciprocal diversity and species‐specificity of some of these mutualisms have been cited as evidence that coevolution between plants and pollinators drives their mutual diversification. However, the mechanisms by which these mutualisms diversify have received less attention. In this paper, we review key hypotheses about how these mutualisms diversify and what role coevolution between plants and pollinators may play in this process. We find that most species‐rich brood pollination mutualisms show significant phylogenetic congruence at high taxonomic scales, but there is limited evidence for the processes of both cospeciation and duplication, and there are no unambiguous examples known of strict‐sense contemporaneous cospeciation. Allopatric speciation appears important across multiple systems, particularly in the insects. Host‐shifts appear to be common, and widespread host‐shifts by pollinators may displace other pollinator lineages. There is relatively little evidence for a “coevolution through cospeciation” model or that coevolution promotes speciation in these systems. Although we have made great progress in understanding the mechanisms by which brood pollination mutualisms diversify, many opportunities remain to use these intriguing symbioses to understand the role of biotic interactions in generating biological diversity.     

Developing Subdomain Allocation Algorithms Based on Spatial and Communicational Constraints to Accelerate Dust Storm Simulation

Dust storm has serious disastrous impacts on environment, human health, and assets. The developments and applications of dust storm models have contributed significantly to better understand and predict the distribution, intensity and structure of dust storms. However, dust storm simulation is a data and computing intensive process. To improve the computing performance, high performance computing has been widely adopted by dividing the entire study area into multiple subdomains and allocating each subdomain on different computing nodes in a parallel fashion. Inappropriate allocation may introduce imbalanced task loads and unnecessary communications among computing nodes. Therefore, allocation is a key factor that may impact the efficiency of parallel process. An allocation algorithm is expected to consider the computing cost and communication cost for each computing node to minimize total execution time and reduce overall communication cost for the entire simulation. This research introduces three algorithms to optimize the allocation by considering the spatial and communicational constraints: 1) an Integer Linear Programming (ILP) based algorithm from combinational optimization perspective; 2) a K-Means and Kernighan-Lin combined heuristic algorithm (K&K) integrating geometric and coordinate-free methods by merging local and global partitioning; 3) an automatic seeded region growing based geometric and local partitioning algorithm (ASRG). The performance and effectiveness of the three algorithms are compared based on different factors. Further, we adopt the K&K algorithm as the demonstrated algorithm for the experiment of dust model simulation with the non-hydrostatic mesoscale model (NMM-dust) and compared the performance with the MPI default sequential allocation. The results demonstrate that K&K method significantly improves the simulation performance with better subdomain allocation. This method can also be adopted for other relevant atmospheric and numerical modeling.     

利用改进的3DMax算法重构染色体3D结构

近年来,随着高通量染色体构象捕获(Hi-C)等技术的发展和高通量测序成本的降低,全基因组交互作用的数据量快速增长,交互作用图谱分辨率不断提高,促使染色体和基因组三维结构建模的研究取得了很大进展,已经提出了几种从染色体构象捕捉数据中构建单个染色体或整个基因组结构的方法。文中通过对在 Hi-C 数据基础上对染色体三维结构重建的相关文献进行分析,总结了重建染色体三维空间结构的经典算法3DMax的原理,并且提出了一种新的随机梯度上升算法:XNadam,是Nadam优化方法的一个变体,将其应用于3DMax算法中,以便提高3DMax算法的性能,从而用于预测染色体三维结构。