Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Highly Accessed Research article

Parasitic plants have increased rates of molecular evolution across all three genomes

Lindell Bromham*, Peter F Cowman and Robert Lanfear

Author Affiliations

Centre for Macroevolution and Macroecology, Research School of Biology, Australian National University, Canberra, A.C.T, 0200, Australia

For all author emails, please log on.

BMC Evolutionary Biology 2013, 13:126  doi:10.1186/1471-2148-13-126


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/13/126


Received:21 March 2013
Accepted:5 June 2013
Published:19 June 2013

© 2013 Bromham et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Theoretical models and experimental evidence suggest that rates of molecular evolution could be raised in parasitic organisms compared to non-parasitic taxa. Parasitic plants provide an ideal test for these predictions, as there are at least a dozen independent origins of the parasitic lifestyle in angiosperms. Studies of a number of parasitic plant lineages have suggested faster rates of molecular evolution, but the results of some studies have been mixed. Comparative analysis of all parasitic plant lineages, including sequences from all three genomes, is needed to examine the generality of the relationship between rates of molecular evolution and parasitism in plants.

Results

We analysed DNA sequence data from the mitochondrial, nuclear and chloroplast genomes for 12 independent evolutionary origins of parasitism in angiosperms. We demonstrated that parasitic lineages have a faster rate of molecular evolution than their non-parasitic relatives in sequences for all three genomes, for both synonymous and nonsynonymous substitutions.

Conclusions

Our results prove that raised rates of molecular evolution are a general feature of parasitic plants, not confined to a few taxa or specific genes. We discuss possible causes for this relationship, including increased positive selection associated with host-parasite arms races, relaxed selection, reduced population size or repeated bottlenecks, increased mutation rates, and indirect causal links with generation time and body size. We find no evidence that faster rates are due to smaller effective populations sizes or changes in selection pressure. Instead, our results suggest that parasitic plants have a higher mutation rate than their close non-parasitic relatives. This may be due to a direct connection, where some aspect of the parasitic lifestyle drives the evolution of raised mutation rates. Alternatively, this pattern may be driven by an indirect connection between rates and parasitism: for example, parasitic plants tend to be smaller than their non-parasitic relatives, which may result in more cell generations per year, thus a higher rate of mutations arising from DNA copy errors per unit time. Demonstration that adoption of a parasitic lifestyle influences the rate of genomic evolution is relevant to attempts to infer molecular phylogenies of parasitic plants and to estimate their evolutionary divergence times using sequence data.

Keywords:
Mutation; Substitution; Comparative; Population size; Molecular rates

Background

Theoretical models have led to the prediction that parasites should evolve higher mutation rates in order to out-evolve their hosts [1]. This prediction has received support from experimental populations of bacteria [2,3]. But does it also apply to multicellular eukaryotes? A number of molecular phylogenetic studies have noted that parasitic taxa have much longer branch lengths, suggesting elevated rates of molecular evolution [4-7]. While it would not be surprising to find that specific genes involved in host-parasite interaction experience elevated substitution rates because of specific selection pressures [8], these studies report higher substitution rates in “housekeeping” genes that are not specifically connected to host-parasite interaction. It is therefore pertinent to ask whether being a parasite causes an increase in the genome-wide rate of molecular evolution in multicellular organisms.

Parasitic plants represent an ideal case for testing the hypothesis that parasites have faster rates of molecular evolution. They show a range of adaptations to parasitism and vary in the type and degree of interaction with their hosts [9,10]. Non-photosynthetic holoparasites are completely dependent on their host for energy, whereas hemiparasites retain photosynthetic ability. Hemiparasites can be either obligate (host-dependent) or facultative (able to be free-living). Parasitic plants also differ in their physical connection to their host (e.g., through the stems or roots), and the degree of reduction of their own anatomy (e.g., loss of stems, leaves or roots).

Faster substitution rates have been reported for a number of genes in several parasitic plant lineages [6,7,11-14]. Furthermore, the mode of parasitism has been suggested to influence rates of molecular evolution, for example holoparasites have been shown to have faster rates than confamilial hemiparasites in three plastid genes of the diverse broomrape family Orobanchaceae [4]. However, not all studies have demonstrated a general increase in substitution rate in parasitic plants. For example, a study of substitution rates in the 18S rRNA gene found that less than half of the parasitic taxa had elevated rates of molecular evolution (3 out of 16 hemiparasitic taxa, and 8 out of 17 holoparasites [7]), and only some of the nuclear genes in a study of substitution rates in the holoparasitic root parasite Balanophora showed significantly increased rates [13]. So the generality and possible causes of the increase in rates in parasite genomes remain unknown.

There are several reasons why a wide-ranging comparative study of rates of molecular evolution in parasitic plants is needed. Most importantly, a comparative study allows us to move from anecdotal observation to general principle. Although faster rates of molecular evolution have been noted in several different species of parasitic plants, we cannot tell from those studies alone whether accelerated rates are a general consequence of their parasitic lifestyle, or whether the faster rates in these particular species and genes are incidental to parasitism. By examining rate variation for all available parasitic lineages in the same analysis we can test the generality of the pattern.

A phylogenetic comparative study also provides a statistically rigorous way to test the hypothesis that parasitic lifestyle increases the rate of molecular evolution. Previously published studies have used a range of data, different methods, and different statistical approaches. In particular, many have used overlapping comparisons between parasites and non-parasites [6,13,15], or treated each parasitic lineage as an independent data point [7], violating assumptions of statistical independence inherent in the tests [16]. It is therefore important to analyse all of the available data within a single statistical framework to describe general patterns in molecular evolution in all independent lineages of parasitic plants.

A comparative study also allows us to compare all three plant genomes, to check for consistency of patterns of molecular evolution. For example, we might expect that the chloroplast genome would evolve rapidly if released from selective constraint, but it is not clear that the same explanation would apply to genes from the mitochondrial or nuclear genomes [7]. Studies that have reported higher rates in parasitic plants have not tested different possible explanations of this pattern, such as positive selection, relaxed purifying selection, or changes in population size. By comparing different substitutions rates within the same analytical framework, we can examine the reasons for the proposed acceleration in rates in parasitic plant taxa, by teasing apart the patterns due to mutation rate, selection and population processes.

There are many independent origins of parasitic lifestyle in angiosperms [17], so we are able to analyse data from 12 separate lineages of parasites. These lineages vary widely in growth form, life cycle, distribution and level of specialisation to parasitism. We are able to combine sequence data from chloroplast, nuclear, and mitochondrial sequences for representative taxa from all 12 lineages with data from their nonparasitic relatives. This provides sufficient data to allow the relationship between parasitism and faster rates to be examined within a single phylogenetic comparative framework.

Results

Parasitic plant clades had significantly higher substitution rates in all measures of branch length for all three genomes (Figure 1). Ten out of twelve nuclear comparisons showed a faster substitution rate in the parasitic lineage compared to its nonparasitic relative (Wilcoxon signed ranks test: Z = 2.786, p = 0.005). For the mitochondrial genes, all but one comparison had longer branch lengths in the parasitic clade (Z = 2.628, p = 0.009) and both nonsynonymous (dN: Z = 2.942, p = 0.003) and synonymous rates (Z = 2.157, p = 0.031) were significantly higher in the parasitic lineages.

thumbnailFigure 1. Scatter plot of comparisons between phylogenetically independent pairs of parasitic plants and their nonparasitic relatives. Points were calculated as the value of the nonparasitic clade subtracted from the parasitic clade, so that if the parasite has a faster rate then the difference is greater than 0. Each line represents a single substitution class in one of the three plant genomes, and each comparison is represented either by a red dot if the parasitic clade has the greater substitution rate or a blue dot if the nonparasitic clade has the higher substitution rate. The points have been slightly jittered to allow overlapping points to be seen clearly.

For the chloroplast sequences, there were fewer sequences available, and neither the rRNA nor the protein-coding sequences had significantly different rates in parasites under the Wilcoxon signed ranks test. However, there was a clear trend in the data. For those comparisons with 16S rRNA (rrn16), seven out of eight comparisons has faster rates in the parasitic lineage. For the chloroplast protein-coding alignment, five out of six comparisons had faster rates in the parasitic lineage. We tested for a significant pattern across all comparisons by asking whether the parasite had the faster rate in whichever chloroplast genes were available, considering only the sign but not the magnitude of the branch length difference. Ten out of twelve comparisons had faster rates in the parasitic clade in the chloroplast sequences available for that comparison (two tailed sign test, p = 0.039). Because there are only 6 comparisons with chloroplast protein coding genes, it is difficult to draw any conclusions regarding synonymous and nonsynonymous substitution rates, but they do follow the same pattern as the mitochondrial and nuclear genomes, with both raised dS and dN in most parasite lineages, but not increased dN/dS.

We did not find a consistent difference in GC content between sequences from parasitic and non-parasitic taxa for nuclear (signs test: p = 0.77), mitochondrial (p = 0.15) or chloroplast (p = 0.15) alignments. This suggests that GC bias is unlikely to explain the consistent differences in rate between parasites and their non-parasitic relatives.

Discussion

Our analysis provides strong support for the hypothesis that parasitic plants have consistently faster rates of molecular evolution than their non-parasitic relatives in mitochondrial, chloroplast, and nuclear sequences. The pattern is detected in all classes of substitutions tested, including ribosomal rRNA genes, and both synonymous and nonsynonymous substitutions in protein coding genes.

What could increase the rate of substitution so consistently in such varied lineages of parasitic plants? The parasitic plants included in this study range from fungus-like holoparasites, such as the achlorophyllous members of the Apodanthaceae that exist entirely within host tissues except when flowering, to perennial shrubs and trees, such as the sandalwoods (e.g. Nestronia) [18]. Here, we will consider several possible explanations for faster substitution rates in parasites.

Relaxed selection

By drawing on metabolites from another individual, a parasite may be under less pressure to maintain their own means to produce those resources. Indeed, there may be selection pressure to lose unnecessary features, in order to save resources or accelerate the life cycle. Many holoparasites have lost or greatly reduced their leaves and stems, and in extreme cases the main vegetative body of the parasitic plant consists of roots growing as a mycelium-like mass inside the host tissues (e.g. Rafflesiaceae, Cytinaceae, Apodanthaceae). In such circumstances we would expect the genes that make products and structures needed for photosynthesis, such as leaves and chlorophyll, to decay into pseudogenes or be lost from the genome [19-22].

A gene that has no selective pressure to maintain function should acquire substitutions at a rate determined by the mutation rate. In this case, we would expect the nonsynonymous substitution rate to approximate the synonymous rate, and their ratio to be approximately 1. We do not find any evidence to support this effect in our study, where values of the ratio of dN to dS are not consistently closer to one in the parasitic lineages, even in the available protein-coding genes of the chloroplast [14]. It is possible that dN/dS values less than one in the parasitic lineage reflect past (negative) selection. Since our substitution rates are measured from the common ancestor of the sister clades, the genes in the parasitic clade will have been in autotrophic plants for some part of their history. In this case, we might expect dN/dS ratios to be higher (closer to one) in most parasitic plants when compared to their non-parasitic relatives, as they would have accumulated relatively more nonsynonymous substitutions during the parasitic part of their evolutionary history. But there is no clear trend toward higher dN/dS in parasitic plant lineages in any of the sequences: parasitic lineages have higher dN/dS in 9 out of 12 comparisons for mitochondrial genes (Z = 1.294, p = 0.196), and only 1 out of 6 comparisons for chloroplast genes (Z = −1.26, p = 0.21).

Loss of function in photosynthesis genes might not follow a simple decay from the adoption of parasitic lifestyle onwards. Chloroplast genes in different species of Orobanche (broomrape), a non-photosynthetic root parasite, show variable rates of decay into pseudogenes [23], and parasitic plants show different degrees of loss of genes from the plastid genome [20]. Some parasitic plants retain low levels of photosynthetic activity, insufficient to live independently, but enough to tide them over while waiting to find a suitable host [24]. Adaptation to parasitism may come about through changes to gene expression levels, leaving less of a signature in the coding sequences of genes [25].

Chloroplast genes in parasites are not necessarily divorced from the effects of selection. Chloroplasts genomes are not solely directed toward photosynthesis, they also contain genes that play a variety of key metabolic roles [19,26]. In any case, it is difficult to imagine that all the genes in this study are under relaxed purifying selection in parasites. The nuclear genes analysed code for essential parts of the protein synthesis machinery, and the mitochondrial genes are associated with gene expression and metabolism. Furthermore, we would not expect relaxation of selection to have a marked effect on synonymous substitution rates. So it does not seem very likely that the acceleration in both synonymous and nonsynonymous substitution rates across all three genomes in a wide range of parasitic plants can be explained through general relaxed selection on gene function.

Positive selection

Hosts may use a number of different strategies to resist parasitic plants, operating at all stages of parasite establishment, ranging from immune-like responses, expression of parasite-specific toxins, and silencing host genes that make products essential for parasite growth [27]. So parasites will probably have to battle their hosts on multiple fronts, suggesting that many parts of parasite physiology and development, and consequently many different genes, will be under the strong selective pressure of a host-parasite arms-race. Indeed, transcriptomics of tissues at the host-parasite interface suggest that many pathways upregulated in the parasites are “general purpose” genes involved in core processes in both parasite and host species [28]. So selection could raise the substitution rate in a wide variety of genes involved in successful exploitation of hosts, including genes that might be considered “house keeping genes” [5]. The arms race between host and parasite may account for the divergence of parasitic lineages onto different host plant populations, such that parasites can be more successfully grown on their local host plants than those from another population [23].

Positive selection on protein-coding genes could raise the non-synonymous substitution rate, as it could drive the fixation of mutations that altered the functional properties of the protein. It is possible, therefore, that positive selection could explain part of the elevation of non-synonymous substitution rates we observed in parasitic plant lineages. However, we would not expect the synonymous substitution rate to be elevated by a host-parasite arms race of this kind. While some synonymous mutations could experience a degree of selective pressure [29], for example for translational efficiency [30,31], it is difficult to see how this could differ sufficiently between parasitic and autotrophic plants to explain the acceleration in synonymous substitution rates in parasites seen in this study. So if positive selection had raised the non-synonymous substitution rate but left the synonymous substitution rate relatively unchanged in parasitic lineages, then we would expect these lineages to show increased dN/dS relative to non-parasitic lineages. We observe no consistent changes in dN/dS associated with parasitic lineages in our data, suggesting that positive selection is unlikely to explain any significant proportion of the elevated parasitic substitution rates that we have observed.

Reduction in effective population size

Effective population size (Ne) affects patterns and rates of substitution, so should always be considered when seeking an explanation for consistent patterns in molecular evolution [32-34]. Parasitic taxa might undergo more frequent and severe founder effects, which might decrease their long-term effective population size and so increase the rate of substitution [5,35]. Endosymbiotic bacteria and fungi have higher substitution rates than their free-living relatives, which has been interpreted as a result of consistently lower population sizes and frequent population bottlenecks [36]. Could small population size or frequent population bottlenecks explain the faster rates in parasitic plants?

It is not clear that parasitic plants will typically have significantly smaller effective population sizes than their autotrophic relatives. Many are animal pollinated, indeed the flowers may be the most conspicuous part of the life cycle, like the gigantic fly-pollinated Rafflesia flowers, and many species have edible seeds that are likely to be dispersed by mammals or birds [18]. Some parasitic plants are widespread generalists that can tolerate a range of conditions and grow on a variety of host plants, to the extent that they are economically important weeds of many different crops [37]. So there seems no a priori reason to suspect that parasitic plants are less likely to be outcrossed or widely dispersed compared to their autotrophic relatives. In any case, lower average population sizes in parasitic lineages should mimic the effects of relaxed purifying selection, on a genome-wide scale, because selection becomes less efficient at removing slightly deleterious mutations in smaller populations [38]. We don’t detect any evidence of this in our study, because there is no significant increase in the ratio of nonsynonymous to synonymous substitutions in the parasitic taxa.

Increased mutation rate

Parasites will often have much faster rates of molecular evolution than their hosts, potentially due to their much faster rate of generation turnover [39], which may help them get ahead in the arms race [40]. But in this study, we have not compared the rate of molecular evolution of parasites to their hosts, but of parasites to their free-living relatives. Synonymous substitution rate (dS) is commonly interpreted as reflecting the mutation rate, because changes to the gene sequence that do not cause a change in the protein-product are likely to be invisible to selection [41]. However, other factors could influence the synonymous substitution rate, such as GC biased gene conversion [42] and codon-usage bias [29], and synonynomous substitution rate can vary across the genome [29,43]. But while dS is not a perfect predictor of mutation rate, changes in mutation rate should be reflected in changes in the rate of both synonymous and nonsynonymous substitution. So elevated mutation rates in parasites could produce the pattern we observe in this study, raising both the synonymous and non-synonymous substitution rates, but not influencing the dN/dS ratio. If elevated mutation rates can explain the patterns we have observed, we must explain how adopting a parasitic lifestyle could result in an increase in the mutation rate.

Simple models that set up an arms race between host and parasite loci predict that increased mutation rate will allow parasites to “chase” the host’s changing defences [1]. Mutator alleles can speed evolution in a clonally reproducing population even when they are not fixed in the population [44]. However, the increase in mutation predicted by chase models is typically transient [45], as the positively selected mutations recombine away from mutator alleles. So alleles that cause a global increase in mutation rate are unlikely to be maintained by selection for the novel traits they generate [46]. The link between mutation rate and adaptive evolution is complex even in asexual lineages [47,48], and in sexual populations, it is not clear the extent to which adaptation will be mutation limited [49]. It seems unlikely that a sexually reproducing parasitic plant would gain enough advantage from increased mutation rate, in terms of evading host defence, for selection on mutator alleles to outweigh the costs of producing more deleterious mutations.

Alternatively, increased mutation rate in parasites might results from an erosion of the DNA repair systems that maintain replication fidelity and fix incidental damage. For example, it has been suggested that the lower bound of the mutation rate is set by the efficiency of selection for maintenance or improvement to DNA repair, which will be primarily limited by the effective population size [50]. In this case, we might expect to see that dS correlates with Ne, so if parasitic plants had consistently lower population sizes it may increase dS and dN, potentially accounting for the lack of correlation between parasitism and dN/dS. However, it is difficult to assess this explanation without an independent way of estimating Ne for these species.

The two explanations just discussed both concern the role of selection in shaping mutation rates, whether increased positive selection for novelty or less efficient negative selection on DNA copy fidelity and repair. In addition, there are many other species traits that can influence mutation rates, so if these traits differ consistently between parasites and their nonparasitic relatives, it could potentially explain the pattern in rates of molecular evolution. The most obvious candidate is number of DNA replications per unit time.

Species with shorter generation times tend to have higher rates of molecular evolution, presumably because they accumulate more DNA replication errors per unit time [51]. In plants, it has been widely reported that annual plants and herbs have faster substitution rates than perennials and woody plants, respectively [52], though it is not clear exactly how this relates to number of genome copies per unit time [29]. Plant height also correlates with rate of molecular evolution in plants, which may be because taller plants tend to have longer generation times (and so fewer meiosis per unit of time) and slower average growth rates (and so fewer mitoses per unit of time) [53]. Highly reduced holoparasites that do not produce stems, such as Pilostyles and Apodanthes, are clearly much shorter than their freeliving relatives, such as Hibiscus and Theobroma (Additional file 1: Table S1). Similarly, the fungus-like Cynomorium coccineum is much smaller than its nonparasitic relative Peridiscus lucidus, a South American tree species. So if parasitic plants tend to have smaller bodies, shorter lives or more rapid generation turnover than their nonparasitic relatives, then they might also have higher mutation rates because their genomes will undergo more replications per unit time.

Additional file 1: Table S1. Nuclear comparisons. Details of comparisons between parasitic plants and their nonparasitic relatives. For each comparison (Comp.), we list the parasitic species and non-parasitic species included, with the GenBank accession numbers for sequences used in this analysis. We also list the parasitic mode for the parasitic species, whether Holoparasitic (Holo, for nonautotrophic species) or Hemiparasitic (Hemi, for species capable of photosynthesis). Note that some taxa include sequences from congeneric species.

Format: XLSX Size: 41KB Download fileOpen Data

Although we do not have reliable measures of average height of all species in our comparisons, we can make a broad assessment based on growth form. In the majority of comparisons, the parasitic species are clearly smaller than the nonparasitic species – specifically in the Apodanthaceae (comparison 1), Cytinaceae (2), Rafflesiaceae (3), Cynomoriaceae (4), Mitrastemonaceae (6), Lauraceae (10), Hydnoraceae (11), and the clade consisting of the parasitic relatives of the Olacaceae (12). These comparisons all show higher rates in the parasitic species for the nuclear genes and mitochondrial genes except for comparison 6 (Mitrastemonaceae) which has slower rates for the parasitic species for mitochondrial dS and chloroplast 16S rRNA.

While it provides a plausible explanation of the patterns we observe, the plant height hypothesis does not provide a perfect fit to all of the substitution rate differences in this study. For example, the perennial hemiparasite Krameria lanceolata has lower nuclear rates and mitochondrial dS than the annual herb Kallstroemia parviflora, which is consistent with larger size in the parasitic taxon, but it also has faster mitochondrial dN. For the chloroplast comparison, Krameria has lower substitution rates than its larger relatives, the tree Guaiacum sanctum and the woody shrub Larrea cuneifolia. For the remaining two comparisons (9 and 10), there is not a clear difference in body size, but the parasites have the higher substitution rates in all alignments. So while the patterns we observe are broadly consistent with an effect of plant size on rates of mutation rates, the match to the observations is not exact, and there may be other factors at play.

Another potential source of an indirect link between molecular evolution and parasitic lifestyle is through environmental effects. It has been suggested that plant lineages at low latitudes have increased rates of molecular evolution, potentially due to the mutagenic effect of UV radiation, or to higher growth rates and faster generation turnover [53-55]. But latitude or environmental energy is unlikely to provide an explanation for higher substitution rates in parasitic plants. Some parasitic plant groups are primarily tropical, such as the Rafflesiaceae, found in South East Asia. But most of the parasitic lineages included in this study have a broader distribution, including both temperate and tropical species, for example the Orobanchaceae (broom rapes) and Cuscuta (dodder). Similarly, net diversification rate in plants has been linked to both synonymous and nonsynonymous rates [56,57], but since parasitic plant clades tend to be less diverse than their nonparasitic sister clades [58], this seems unlikely to provide an explanation for our results.

Conclusions

The study of parasitic plants has many important outcomes. By reducing the water and nutrients available to host plant growth, parasitic plants have a serious impact on agriculture, which is estimated to impact on the food supply of over 100 million people [27]. Parasitic plants also provide an excellent case study for examining many evolutionary processes, including diversification and specialization, genome evolution, and host-parasite interactions. Demonstrating consistently higher substitution rates in parasitic plants not only provides a window on molecular evolution, it also has practical implications for the use of DNA sequence data in evolutionary and ecological research. Molecular markers have played a key role not only in establishing the phylogeny and systematics of parasitic plants, but in understanding broader evolutionary patterns, such specialization to different host plants. Consistent differences in rate of molecular evolution in parasitic lineages should be considered when interpreting molecular phylogenies and date estimates.

Methods

Data

We gathered DNA sequence data from GenBank for each of the parasitic plant clades identified by Barkman et al. [17]. We focus on angiosperms, so we do not include the gymnosperm Parasitaxus, nor the liverwort Aneura mirabilis, both of which appear to form parasitic attachments to host plants via a fungal partner [20,59]. Nor do we include mycoheterotrophic angiosperm species, such as those in the Ericaceae [60], which derive part of their nutrition through a relationship with fungi, because the line between fully parasitic mycoheterotrophs and plants that have mycorrhizal partners is not always clear [61,62].

We use a sister-clades approach, identifying clades of parasitic plants and their autotrophic (non-parasitic) relatives. Each lineage of the sister clade has had the same amount of time since their last common ancestor to accumulate substitutions, therefore any difference in branch length (estimated number of substitutions since the last common ancestor) is likely to reflect a difference in the net substitution rate [16,56]. If parasitic plants have a higher rate of molecular evolution, then we would expect to see a longer total branch length in the parasitic clade compared to its non-parasitic sister clade [51].

The position of parasitic clades within angiosperm phylogeny has been controversial, and even the higher-level relationships are debated [63-65]. As a consequence, there have been many changes to higher-level taxonomy. For simplicity, we report the GenBank taxonomy specified for each sequence used.

We identified close non-parasitic relatives of each of the parasitic lineages, using published phylogenies and taxonomies as a guide [17,64-72]. For the purposes of this study, we do not need to know the exact phylogenetic position of parasitic taxa, as long as we can be sure that we are comparing phylogenetically independent origins of parasitic lineages to their autotrophic relatives. For each independently evolved parasitic lineage we aimed to select a related autotrophic clade that could be placed outside the parasitic clade with certainty, so we did not always choose the closest relative where phylogeny was uncertain. This gave us 12 phylogenetically independent comparisons between parasitic and non-parasitic lineages (Tables 1 and 2). Each parasitic and non-parasitic sister clade was represented by sequences from up to nine different species (see Additional file 1: Table S1, Additional file 2: Table S2 and Additional file 3: Table S3). We always balanced the number of species in the sister clades, to minimize the effects of node density on substitution rate estimation [16].

Table 1. Estimates of branch length for sister clades of parasitic (P) and autotrophic (non-parasitic, NP) plants for nuclear and mitochondrial sequences

Table 2. Estimates of branch length for sister clades of parasitic (P) and autotrophic (non-parasitic, NP) plants for chloroplast sequences

Additional file 2: Table S2. Mitochondrial comparisons. Details of comparisons between parasitic plants and their nonparasitic relatives. For each comparison (Comp.), we list the parasitic species and non-parasitic species included, with the GenBank accession numbers for sequences used in this analysis. We also list the parasitic mode for the parasitic species, whether Holoparasitic (Holo, for nonautotrophic species) or Hemiparasitic (Hemi, for species capable of photosynthesis). Note that some taxa include sequences from congeneric species.

Format: XLSX Size: 46KB Download fileOpen Data

Additional file 3: Table S3. Chloroplast comparisons.Details of comparisons between parasitic plants and their nonparasitic relatives. For each comparison (Comp.), we list the parasitic species and non-parasitic species included, with the GenBank accession numbers for sequences used in this analysis. We also list the parasitic mode for the parasitic species, whether Holoparasitic (Holo, for nonautotrophic species) or Hemiparasitic (Hemi, for species capable of photosynthesis). Note that some taxa include sequences from congeneric species.

Format: XLSX Size: 40KB Download fileOpen Data

We targeted genes for which sequences were available for the largest number of the parasitic taxa and their non-parasitic relatives: maturase R (matR), cytochrome oxidase subunit I (coxI), NADH dehydrogenase subunit 1 (nad1) and ATPase F1 alpha subunit (atp1) from the mitochondrial genome; 18S small subunit ribosomal RNA (rrn18) and 26S ribosomal RNA (rrn26) from the nuclear genome; ribulose-1, 5-bisphosphate carboxylase/oxygenase large subunit (rbcL), maturase K (matK) and 16S ribosomal RNA (rrn16) from the chloroplast genome. We used the GeneFinder script [73] to find the longest available sequences for each gene on GenBank. Some gene sequences were unavailable for some taxa (in particular, chloroplast sequences were often not available for holoparasites). A complete list of the sequences, along with their GenBank accession numbers, is available in the Supplementary Material (Additional file 1: Table S1, Additional file 2: Table S2 and Additional file 3: Table S3).

Analysis

We aligned sequences in Geneious version 6.0 [74]. We aligned protein-coding exons in frame, using the amino acid translations as a guide, with the MUSCLE translation alignment plugin [75], and we aligned non-coding sequences using the standard MUSCLE alignment plugin in Geneious. We then adjusted all alignments by eye. When we identified pseudogenes in chloroplast protein sequences we removed them from the alignments. We produced three separate genome alignments by concatenating the genes from each genome (mitochondrial, chloroplast, and nuclear). For each genome alignment, substitution models and partitioning schemes were chosen using AICc in Partition Finder [76]. The best substitution models and partition schemes were used for phylogenetic estimation in RAxML v7.0.4 [77]. Separate maximum likelihood analyses were run for each alignment.

Next we created a single maximum likelihood (ML) tree topology and corresponding alignment for each of the three genomes (nuclear, plastid, and mitochondrial) for each of the 12 sister pairs (36 topologies and alignments in total). For each comparison, we then extracted the parasitic and non-parasitic taxa plus a closely related non-parasite outgroup from a different order from each of the genome trees using the APE package [78], and produced corresponding alignment files. For each alignment, we only included gene regions that had coverage for both the parasitic and non-parasitic clades. We removed the small introns of the coxI and nad1 genes in the mitochondrial alignment.

We estimated total branch length for all alignments using maximum likelihood and a GTR + G model in the baseml program of the PAML package [79]. For protein-coding genes (mitochondrial matR, coxI, nad1 and atp1; chloroplast rbcL and matK), we estimated synonymous, and nonsynonymous branch lengths using the GY94 codon substitution model [80] in the codeml program of the PAML package [79], with dN/dS values free to vary across the tree. In addition, we estimated clade-specific dN/dS ratios (ω) for the parasitic and non-parasitic clade in each comparison.

To compare the rate of accumulation of substitutions in parasitic lineages to their nonparasitic sister lineage, we calculated the total branch length of the parasitic and non-parasitic clades in each sister pair using a python script. Total branch lengths for ribosomal RNA genes and protein-coding genes, and total dS and dN branch lengths for protein-coding genes were calculated by summing all edge lengths for each clade from the shared node with its sister clade. Since estimates of substitution rates can be influenced by base composition bias, we also calculated GC content for each sister clade. We tested for a significant association between parasitism and rate of molecular evolution using a Wilcoxon signed ranks test in the SPlus package v8.2.

Availability of supporting data

Data has been deposited in the Data Dryad Repository: http://dx.doi.org/10.5061/dryad.fc74k webcite.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed to data collection and analysis and manuscript preparation and all authors read and approved the final manuscript.

Acknowledgements

We thank Katherine Byron, Simon Ho and Art Wuster for their contributions to earlier versions of this project.

References

  1. Haraguchi Y, Sasaki A: Host-parasite arms-race in mutation modifications - indefinate escalation despite a heavy load?

    J Theor Biol 1996, 183(2):121-137. PubMed Abstract | Publisher Full Text OpenURL

  2. Pal C, Maciá MD, Oliver A, Schachar I, Buckling A: Coevolution with viruses drives the evolution of bacterial mutation rates.

    Nature 2007, 450(7172):1079-1081. PubMed Abstract | Publisher Full Text OpenURL

  3. Paterson S, Vogwill T, Buckling A, Benmayor R, Spiers AJ, Thomson NR, Quail M, Smith F, Walker D, Libberton B: Antagonistic coevolution accelerates molecular evolution.

    Nature 2010, 464(7286):275-278. PubMed Abstract | Publisher Full Text OpenURL

  4. Young ND, De Pamphilis CW: Rate variation in parasitic plants: correlated and uncorrelated patterns among plastid genes of different function.

    BMC Evol Biol 2005, 5:16. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  5. Dowton M, Austin AD: Increased genetic diversity in mitochondrial genes is correlated with the evolution of parasitism in the Hymenoptera.

    J Mol Evol 1995, 41:958-965. PubMed Abstract OpenURL

  6. Duff RJ, Nickrent DL: Characterization of mitochondrial small-subunit riibosomal RNAs from holoparasitic plants.

    J Mol Evol 1997, 45(6):631. PubMed Abstract | Publisher Full Text OpenURL

  7. Lemaire B, Huysmans S, Smets E, Merckx V: Rate accelerations in nuclear 18S rDNA of mycoheterotrophic and parasitic angiosperms.

    J Plant Res 2011, 124(5):561-576. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Jiggins FM, Hurst GDD, Yang Z: Host-symbiont conflicts: positive selection on an outer membrane protein of parasitic but not mutualistic Rickettsiaceae.

    Mol Biol Evol 2002, 19:1341-1349. PubMed Abstract | Publisher Full Text OpenURL

  9. Kuijt J: The biology of flowering plants, vol. Berkeley: University of California Press; 1969. OpenURL

  10. Westwood JH, Yoder JI, Timko MP, De Pamphilis CW: The evolution of parasitism in plants.

    Trends Plant Sci 2010, 15(4):227-235. PubMed Abstract | Publisher Full Text OpenURL

  11. Nickrent DL, Starr EM: High rates of nucleotide substitution in nuclear small-subunit (18S) rDNA from holoparasitic flowering plants.

    J Mol Evol 1994, 39(1):62. PubMed Abstract OpenURL

  12. Nickrent DL, Garcia MA: On the brink of holoparasitism: plastome evolution in dwarf mistletoes (Arceuthobium, Viscaceae).

    J Mol Evol 2009, 68(6):603-615. PubMed Abstract | Publisher Full Text OpenURL

  13. Su H-J, Hu J-M: Rate heterogeneity in six protein-coding genes from the holoparasite Balanophora (Balanophoraceae) and other taxa of Santalales.

    Ann Bot 2012, 110:1137-1147. PubMed Abstract | Publisher Full Text OpenURL

  14. McNeal JR, Kuehl JV, Boore JL, De Pamphilis CW: Complete plastid genome sequences suggest strong selection for retention of photosynthetic genes in the parasitic plant genus Cuscuta.

    BMC Plant Biol 2007, 7(1):57. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  15. Wolfe KH, Katz-Downie DS, Morden CW, Palmer JD: Evolution of the plastid ribosomal RNA operon in a nongreen parasitic plant: accelerated sequence evolution, altered promoter structure, and tRNA pseudogenes.

    Plant Mol Biol 1992, 18(6):1037-1048. PubMed Abstract | Publisher Full Text OpenURL

  16. Lanfear R, Welch JJ, Bromham L: Watching the clock: Studying variation in rates of molecular evolution.

    Trends Ecol Evol 2010, 25(9):495-503. PubMed Abstract | Publisher Full Text OpenURL

  17. Barkman TJ, McNeal JR, Lim SH, Coat G, Croom HB, Young ND, De Pamphilis CW: Mitochondrial DNA suggests at least 11 origins of parasitism in angiosperms and reveals genomic chimerism in parasitic plants.

    BMC Evol Biol 2007, 7:248. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  18. Nickrent DL: Parasitic plants connection. Southern Illinois: Southern Illinois University; 1997. OpenURL

  19. Krause K: From chloroplasts to cryptic plastids: evolution of plastid genomes in parasitic plants.

    Curr Genet 2008, 54(3):111-121. PubMed Abstract | Publisher Full Text OpenURL

  20. Wickett NJ, Zhang Y, Hansen SK, Roper JM, Kuehl JV, Plock SA, Wolf PG, DePamphilis CW, Boore JL, Goffinet B: Functional gene losses occur with minimal size reduction in the plastid genome of the parasitic liverwort Aneura mirabilis.

    Mol Biol Evol 2008, 25(2):393-401. PubMed Abstract | Publisher Full Text OpenURL

  21. Randle CP, Wolfe AD: The evolution and expression of rbcL in holoparasitic sister-genera Harveya and Hyobanche (Orobanchaceae).

    Am J Bot 2005, 92(9):1575-1585. PubMed Abstract | Publisher Full Text OpenURL

  22. Young ND, DePamphilis CW: Purifying selection detected in the plastid gene matK and flanking ribozyme regions within a group II intron of nonphotosynthetic plants.

    Mol Biol Evol 2000, 17(12):1933-1941. PubMed Abstract | Publisher Full Text OpenURL

  23. Thorogood CJ, Rumsey FJ, Harris SA, Hiscock SJ: Host-driven divergence in the parasitic plant Orobanche minor Sm. (Orobanchaceae).

    Mol Ecol 2008, 17(19):4289-4303. PubMed Abstract | Publisher Full Text OpenURL

  24. Van der Kooij T, Krause K, Dörr I, Krupinska K: Molecular, functional and ultrastructural characterisation of plastids from six species of the parasitic flowering plant genus Cuscuta.

    Planta 2000, 210(5):701-707. PubMed Abstract | Publisher Full Text OpenURL

  25. Funk H, Berg S, Krupinska K, Maier U, Krause K: Complete DNA sequences of the plastid genomes of two parasitic flowering plant species, Cuscuta reflexa and Cuscuta gronovii.

    BMC Plant Biol 2007, 7(1):45. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  26. Barbrook AC, Howe CJ, Purton S: Why are plastid genomes retained in non-photosynthetic organisms?

    Trends Plant Sci 2006, 11(2):101-108. PubMed Abstract | Publisher Full Text OpenURL

  27. Yoder JI, Scholes JD: Host plant resistance to parasitic weeds; recent progress and bottlenecks.

    Curr Opp Plant Biol 2010, 13(4):478-484. Publisher Full Text OpenURL

  28. Honaas L, Wafula E, Yang Z, Der J, Wickett N, Altman N, Taylor C, Yoder J, Timko M, Westwood J, et al.: Functional genomics of a generalist parasitic plant: Laser microdissection of host-parasite interface reveals host-specific patterns of parasite gene expression.

    BMC Plant Biol 2013, 13(1):9. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  29. Gaut B, Yang L, Takuno S, Eguiarte LE: The patterns and causes of variation in plant nucleotide substitution rates.

    Annu Rev Ecol Evol Syst 2011, 42:245-266. Publisher Full Text OpenURL

  30. Gu WJ, Wang XF, Zhai CY, Xie XY, Zhou T: Selection on synonymous sites for increased accessibility around miRNA binding sites in plants.

    Mol Biol Evol 2012, 29(10):3037-3044. PubMed Abstract | Publisher Full Text OpenURL

  31. Qiu S, Bergero R, Zeng K, Charlesworth D: Patterns of codon usage bias in Silene latifolia.

    Mol Biol Evol 2011, 28(1):771-780. PubMed Abstract | Publisher Full Text OpenURL

  32. Lynch M: The origins of genome architecture. Sunderland, Mass: Sinauer Assoc; 2007. OpenURL

  33. Woolfit M: Effective population size and the rate and pattern of nucleotide substitutions.

    Biol Lett 2009, 5(3):417-420. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Charlesworth B: Effective population size and patterns of molecular evolution and variation.

    Nature Rev Genet 2009, 10:195-205. PubMed Abstract | Publisher Full Text OpenURL

  35. Ebert D: Experimental evolution of parasites.

    Science 1998, 282(5393):1432-1436. PubMed Abstract | Publisher Full Text OpenURL

  36. Woolfit M, Bromham L: Increased rates of sequence evolution in endosymbiotic bacteria and fungi with small effective population sizes.

    Mol Biol Evol 2003, 20(9):1545-1555. PubMed Abstract | Publisher Full Text OpenURL

  37. Musselman LJ: The biology of Striga, Orobanche, and other root-parasitic weeds.

    Annu Rev Phytopathol 1980, 18(1):463-489. Publisher Full Text OpenURL

  38. Lahti DC, Johnson NA, Ajie BC, Otto SP, Hendry AP, Blumstein DT, Coss RG, Donohue K, Foster SA: Relaxed selection in the wild.

    Trends Ecol Evol 2009, 24(9):487-496. PubMed Abstract | Publisher Full Text OpenURL

  39. Hafner M, Sudman P, Villablanca F, Spradling T, Demastes J, Nadler S: Disparate rates if molecular evolution in cospeciating hosts and parasites.

    Science 1994, 265:1087-1090. PubMed Abstract | Publisher Full Text OpenURL

  40. Kaltz O, Shykoff JA: Local adaptation in host–parasite systems.

    Heredity 1998, 81(4):361-370. Publisher Full Text OpenURL

  41. Lanfear R, Ho SYW, Love D, Bromham L: Mutation rate influences diversification rate in birds.

    Proc natl Acad Sci USA 2010, 107(47):20423-20428. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Lartillot N: Interaction between Selection and Biased Gene Conversion in Mammalian Protein-Coding Sequence Evolution Revealed by a Phylogenetic Covariance Analysis.

    Mol Biol Evol 2013, 30(2):356-368. PubMed Abstract | Publisher Full Text OpenURL

  43. DeRose-Wilson L, Gaut B: Transcription-related mutations and GC content drive variation in nucleotide substitution rates across the genomes of Arabidopsis thaliana and Arabidopsis lyrata.

    BMC Evol Biol 2007, 7(1):66. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  44. Taddei F, Radman M, Maynard-Smith J, Toupance B, Gouyon P, Godelle B: Role of mutator alleles in adaptive evolution.

    Nature 1997, 387(6634):700-702. PubMed Abstract | Publisher Full Text OpenURL

  45. Giraud A, Matic I, Tenaillon O, Clara A, Radman M, Fons M, Taddei F: Costs and benefits of high mutation rates: adaptive evolution of bacteria in the mouse gut.

    Science 2001, 291(5513):2606-2608. PubMed Abstract | Publisher Full Text OpenURL

  46. Metzgar D, Wills C: Evidence for the adaptive evolution of mutation rates.

    Cell 2000, 101(6):581. PubMed Abstract | Publisher Full Text OpenURL

  47. Rainey PB: Evolutionary genetics: the economics of mutation.

    Curr Biol 1999, 9(10):R371-R373. PubMed Abstract | Publisher Full Text OpenURL

  48. Denamur E, Tenaillon O, Deschamps C, Skurnik D, Ronco E, Gaillard JL, Picard B, Branger C, Matic I: Intermediate mutation frequencies favor evolution of multidrug resistance in Escherichia coli.

    Genetics 2005, 171(2):825-827. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  49. Weissman DB, Barton NH: Limits to the rate of adaptive substitution in sexual populations.

    PLoS Genet 2012, 8(6):e1002740. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  50. Lynch M: The Lower Bound to the Evolution of Mutation Rates.

    Genome Biol Evol 2011, 3:1107-1118. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  51. Bromham L: Why do species vary in their rate of molecular evolution?

    Biol Lett 2009, 5:401-404. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Smith SA, Donoghue MJ: Rates of molecular evolution are linked to life history in flowering plants.

    Science 2008, 322:86-89. PubMed Abstract | Publisher Full Text OpenURL

  53. Lanfear R, Ho SYW, Davies TJ, Moles ATA L, Swenson NG, Warman L, Zanne AE, Allen AP: Taller plants have lower rates of molecular evolution: the rate of mitosis hypothesis.

    Nature Communications 2013, 4:1879. PubMed Abstract | Publisher Full Text OpenURL

  54. Davies TJ, Savolainen V, Chase MW, Moat J, Barraclough TG: Environmental energy and evolutionary rates in flowering plants.

    Proc R Soc Lond B 2004, 271:2195-2220. Publisher Full Text OpenURL

  55. Wright SD, Keeling J, Gillman L: The road from Santa Rosalia: a faster tempo of evolution in tropical climates.

    Proc natl Acad Sci USA 2006, 103:7718-7722. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  56. Barraclough TG, Savolainen V: Evolutionary rates and species diversity in flowering plants.

    Evolution 2001, 55:677-683. PubMed Abstract | Publisher Full Text OpenURL

  57. Duchene D, Bromham L: Rates of molecular evolution and diversification in plants: chloroplast substitution rates correlate with species richness in the Proteaceae.

    BMC Evol Biol 2013, 13:65. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  58. Hardy NB, Cook LG: Testing for Ecological Limitation of Diversification: A Case Study Using Parasitic Plants.

    Am Nat 2012, 180(4):438-449. PubMed Abstract | Publisher Full Text OpenURL

  59. Feild TS, Brodribb TJ: A unique mode of parasitism in the conifer coral tree Parasitaxus ustus (Podocarpaceae).

    Plant Cell Environ 2005, 28(10):1316-1325. Publisher Full Text OpenURL

  60. Braukmann T, Stefanovicá S: Plastid genome evolution in mycoheterotrophic Ericaceae.

    Plant Mol Biol 2012, 1:16. OpenURL

  61. Leake JR, Cameron DD: Physiological ecology of mycoheterotrophy.

    New Phytol 2010, 185(3):601-605. PubMed Abstract | Publisher Full Text OpenURL

  62. Merckx V, Freudenstein JV: Evolution of mycoheterotrophy in plants: a phylogenetic perspective.

    New Phytol 2010, 185(3):605-609. PubMed Abstract | Publisher Full Text OpenURL

  63. Filipowicz N, Renner SS: The worldwide holoparasitic Apodanthaceae confidently placed in the Cucurbitales by nuclear and mitochondrial gene trees.

    BMC Evol Biol 2010, 10(1):219. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  64. Nickrent DL, Der JP, Anderson FE: Discovery of the photosynthetic relatives of the “Maltese mushroom” Cynomorium.

    BMC Evol Biol 2005, 5(1):38. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  65. Nickrent D, Blarer A, Qiu Y-L, Vidal-Russell R, Anderson F: Phylogenetic inference in Rafflesiales: the influence of rate heterogeneity and horizontal gene transfer.

    BMC Evol Biol 2004, 4:1471-2148. OpenURL

  66. Barkman TJ, Lim S-H, Salleh KM, Nais J: Mitochondrial DNA sequences reveal the photosynthetic relatives of Rafflesia, the world‘s largest flower.

    Proc natl Acad Sci USA 2004, 101:787-792. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  67. Vr M, Nickrent DL: Molecular phylogenetic relationships of Olacaceae and related Santalales.

    Syst Bot 2008, 33(1):97-106. Publisher Full Text OpenURL

  68. Olmstead RG, Wolfe AD, Young ND, Elisons WJ, Reeves PA: Disintegration of the Scrophulariaceae.

    Am J Bot 2001, 88(2):348-361. PubMed Abstract | Publisher Full Text OpenURL

  69. Rohwer JG: Toward a phylogenetic classification of the Lauraceae: evidence from matK sequences.

    Syst Bot 2000, 25(1):60-71. Publisher Full Text OpenURL

  70. Sheahan M, Chase M: A phylogenetic analysis of Zygophyllaceae R. Br. based on morphological, anatomical and rbcL DNA sequence data.

    Bot J Linn Soc 1996, 122:279-300. OpenURL

  71. Soltis D, Soltis P, Chase M, Mort M, Albach D, Zanis M, Savolainen V, Hahn W, Hoot S, Fay M: Angiosperm phylogeny inferred from 18s rDNA, rbcL, and atpB sequences.

    Bot J Linn Soc 2000, 133:381-461. OpenURL

  72. Stefanovic S, Olmstead RG: Testing the phylogenetic position of a parasitic plant (Cuscuta, Convolvulaceae, Asteridae): Bayesian inference and the parametric bootstrap on data drawn from three genomes.

    Syst Biol 2004, 53(3):384-399. PubMed Abstract | Publisher Full Text OpenURL

  73. Lanfear R, Bromham L: Estimating phylogenies for species assemblages: a complete phylogeny for the past and present native birds of New Zealand.

    Mol Phylog Evol 2011, 61:958-963. Publisher Full Text OpenURL

  74. Drummond AJ, Ashton B, Buxton S, Cheung M, Cooper A, Duran C, Field M, Heled J, Kearse M, Markowitz S, Moir R, et al.:

    Geneious. 2011.

    http://www.geneious.com webcite

    OpenURL

  75. Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput.

    Nucleic Acids Res 2004, 32:1792-1797. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  76. Lanfear R, Calcott B, Ho SYW, Guindon S: PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses.

    Molecular Biology and Evolution 2012, 29(6):1695-1701. PubMed Abstract | Publisher Full Text OpenURL

  77. Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models.

    Bioinformatics 2006, 22(21):2688-2690. PubMed Abstract | Publisher Full Text OpenURL

  78. Paradis E, Claude J: K. S: APE: analyses of phylogenetics and evolution in R language.

    Bioinformatics 2004, 20:289-290. PubMed Abstract | Publisher Full Text OpenURL

  79. Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood.

    Mol Biol Evol 2007, 24:1586-1591.

    http://abacus.gene.ucl.ac.uk/software/paml.html webcite

    PubMed Abstract | Publisher Full Text OpenURL

  80. Goldman N, Yang Z: A codon-based model of nucleotide substitution for protein-coding DNA sequences.

    Mol Biol Evol 1994, 11:725-736. PubMed Abstract | Publisher Full Text OpenURL