- Research article
- Open Access
A multilocus timescale for oomycete evolution estimated under three distinct molecular clock models
BMC Evolutionary Biology volume 14, Article number: 101 (2014)
Molecular clock methodologies allow for the estimation of divergence times across a variety of organisms; this can be particularly useful for groups lacking robust fossil histories, such as microbial eukaryotes with few distinguishing morphological traits. Here we have used a Bayesian molecular clock method under three distinct clock models to estimate divergence times within oomycetes, a group of fungal-like eukaryotes that are ubiquitous in the environment and include a number of devastating pathogenic species. The earliest fossil evidence for oomycetes comes from the Lower Devonian (~400 Ma), however the taxonomic affinities of these fossils are unclear.
Complete genome sequences were used to identify orthologous proteins among oomycetes, diatoms, and a brown alga, with a focus on conserved regulators of gene expression such as DNA and histone modifiers and transcription factors. Our molecular clock estimates place the origin of oomycetes by at least the mid-Paleozoic (~430-400 Ma), with the divergence between two major lineages, the peronosporaleans and saprolegnialeans, in the early Mesozoic (~225-190 Ma). Divergence times estimated under the three clock models were similar, although only the strict and random local clock models produced reliable estimates for most parameters.
Our molecular timescale suggests that modern pathogenic oomycetes diverged well after the origin of their respective hosts, indicating that environmental conditions or perhaps horizontal gene transfer events, rather than host availability, may have driven lineage diversification. Our findings also suggest that the last common ancestor of oomycetes possessed a full complement of eukaryotic regulatory proteins, including those involved in histone modification, RNA interference, and tRNA and rRNA methylation; interestingly no match to canonical DNA methyltransferases could be identified in the oomycete genomes studied here.
Eukaryotic diversity is primarily microbial, with multicellularity restricted to a few distinct lineages (plants, animals, fungi, and some algae). While the Proterozoic fossil record contains an abundance of organic-walled, often ornamented, microfossils interpreted as eukaryotes, evidence for the origins and diversification of specific lineages of microbial eukaryotes is rare, especially for those groups with few diagnostic morphological characters . Molecular clock methods therefore provide the only avenue for elucidating the evolutionary history of some lineages. With the recognition that a single rate (“strict”) molecular clock as originally proposed by Zuckerkandl and Pauling [2, 3] was often inadequate in light of rate variation among organisms, early studies suggested the use of local clocks or the removal of lineages that violated the assumption of rate homogeneity (reviewed in ). The continued development of molecular clock methodologies over the past two decades has allowed for the estimation of divergence times under more complex models of rate variation. Initial “relaxed clock” methods, such as non-parametric rate smoothing  and penalized likelihood , allowed rates to vary but sought to minimize large differences between parent and descendent branches. Additionally, Bayesian relaxed clock methods allow rates to vary among lineages but assume autocorrelation by drawing the rate of a descendent branch from a distribution whose mean is determined by the rate of the parent branch [7, 8]; other Bayesian methods relax this assumption of autocorrelation for the co-estimation of phylogeny and divergence times . Most recently, a random local clock model approach has been proposed which allows rate changes to occur along any branch in a phylogeny; this method allows users to directly test various local clock scenarios against a strict clock model of no rate changes .
In addition to improved modeling of rate variation, newer molecular clock methods are also able to better incorporate calibration uncertainty into the estimation of divergence times. Early methods treated fossil calibrations as fixed points (from which rates were derived); newer methods utilize probability distributions to better reflect the paleontological uncertainty of a fossil’s phylogenetic position in relation to modern organisms [11, 12], as well as variance around the numerical age of geologic formations. However, some authors have already shown that modeling fossil probability distributions under different assumptions can have significant impacts on divergence time estimation , illustrating that rate calibration is still an important source of potential error in molecular clock studies.
In this study we have focused on the fungal-like oomycetes (Peronosporomycetes sensu), a group of heterotrophic eukaryotes closely related to diatoms, brown algae, and other stramenopiles . A close relationship among stramenopiles, alveolates, and several photosynthetic eukaryotes with red algal-derived plastids was previously suggested as the supergroup Chromalveolata . However molecular studies have supported a grouping of stramenopiles and alveolates with the non-photosynthetic rhizarians (“SAR” sensu), excluding other photosynthetic lineages; the recently revised eukaryote classification has now formalized the Sar supergroup . Many oomycetes are saprotrophic in aquatic and terrestrial ecosystems, however several devastating pathogens are known, such as Phytophthora infestans, the causal agent of late blight in solanaceous plant hosts . Some orders are primarily pathogenic, such as the Peronosporales and Albuginales, while others are composed of both pathogenic and saprotrophic members, such as the Pythiales, Saprolegniales, Leptomitales, and Rhipidiales . Several basal lineages, such as the Eurychasmales and Haliphthorales, are known primarily as pathogens of marine algae and crustaceans, leading some to suggest that the oomycetes may be “hard-wired” for pathogenic lifestyles .
The earliest robust fossil evidence of oomycetes comes from the Lower Devonian (Pragian, ~408 Ma) Rhynie Chert . Thick-walled, ornamented structures interpreted as oogonium-antheridium complexes , as well as thin-walled polyoosporous oogonia , are well preserved in association with degraded plant debris and cyanobacteria-dominated microbial mats. More recent oomycete fossils occur in the Carboniferous, where evidence for endophytic  and perhaps parasitic [25, 26] interactions with plant hosts is more compelling. Additionally, the fossil species Combresomyces cornifer originally described from Lower Carboniferous chert in central France  has also been identified in Middle Triassic silicified peat from Antarctica , providing an intriguing example of geographic range and morphological stasis over roughly 90 million years of oomycete evolution .
This is the first study to estimate divergence times within the oomycetes using molecular clock methods. Previous studies have typically included a single representative within a larger study of eukaryotic evolution [30–32], or have used oomycetes to root the analysis [33, 34]. As there is little a priori information on the tempo of evolution within oomycetes, here we estimate divergence times under three distinct molecular clock models: a single-rate strict clock, a relaxed clock with uncorrelated rates modeled under a lognormal distribution (UCLD), and a random local clock model. The availability of several complete genome sequences for oomycetes, diatoms, and a brown alga allowed us to carefully curate a dataset of 40 orthologs for divergence time estimation; we chose to focus on known regulators of eukaryotic gene expression to investigate their presence and level of conservation within pathogenic oomycetes. While the performance of the three models differed, the estimated divergence times suggested that oomycetes diverged from other stramenopiles by at least the mid-Paleozoic, and that two major lineages, the peronosporaleans and saprolegnialeans, diverged in the early Mesozoic, approximately 200 Ma after the first appearance of oomycetes in the fossil record.
Regulators of gene expression in oomycetes
Complete genome sequences from eighteen species were examined (Table 1). A total of 70 genes involved in the regulation of gene expression were examined for homology in Phytophthora infestans (Table 2); homologs of two genes (Drosha-like; TFIIH, Ssl1 subunit) could not be identified in P. infestans but were present in other oomycetes. In general, oomycetes possess a full complement of canonical transcription factors and genes involved in chromatin modification, including multiple histone acetyltransferases, deacetylases, and methyltransferases (Table 2). Proteins known to be involved in post-transcriptional gene silencing  were identified in our search, including homologs of Argonaut, Dicer, RNA-dependent RNA polymerase, double-stranded RNA binding proteins, and an RNaseIII-domain containing protein (Table 2). A recent study has shown that these genes are expressed and functional in P. infestans. However, unlike the previous study, we were able to identify a second Dicer-like homolog in the genomes of other oomycetes that is absent in P. infestans; these sequences showed more similarity to human and Drosophila Drosha proteins than to other Dicer homologs (data not shown). Two distinct groups of Argonaut proteins were identified in the oomycetes, as well as two types of double-stranded RNA binding proteins (Table 2). While no homologs to canonical eukaryotic DNA methyltransferases could be identified, a homolog of DNA methyltransferase 1-associated protein was present in all the genomes analyzed here. Several genes involved in RNA methylation were also found (Table 2).
Divergence time analyses
Robust orthology relationships could be determined for 52 out of the initial 70 datasets; 40 of these datasets contained minimal missing data and were used to estimate divergence times (see Additional file 1 for a list of genes included in the analysis). Calibration priors were modeled with a gamma distribution in order to assign higher probabilities to divergence times somewhat older than the hard bound (offset value); initial tests with lognormal priors produced very similar divergence times (data not shown). Five independent analyses of 50 million generations each were run under each of the three models, with the random local clock model being the most computationally intensive. Strict clock and UCLD analyses run on an iMac (10.8.5) desktop with a 2.7 GHz Intel core i5 processor took approximately seven days. Random local clock analyses run on a Linux (Mint14) desktop with a 3.3 GHz Xeon quad core processor took approximately 30 days. Posterior distributions on parameters were identical across all five runs under the strict clock model. Parameter distributions were consistent and overlapping for all five runs under the UCLD model with only one run deviating for the estimate of the root height (700 Ma versus approximately 500 Ma in the other four runs), however all runs showed weak evidence of convergence even after 50 million generations. One run under the random local clock model failed to converge; of the four successful runs, parameter distributions were consistent and overlapping with only one run deviating for the rate estimate (1.76 × 10−3 versus 1.88 × 10−3 for the other three runs). Log and tree files for two of the five runs with the highest effective sample size (ESS) for the likelihood parameter were then combined; under the strict clock model, all five runs performed equally, so the first two runs were combined. Analyses run without data (Prior Only) resulted in time estimates that were markedly different from those obtained with the full dataset for the majority of nodes (Table 3), suggesting that our divergence time estimates were driven by the data themselves and not by settings on the calibration priors. Divergence times among oomycete lineages were consistent among all three models (Table 3), however estimates under the UCLD model may have been influenced by poor mixing as several parameters showed ESS values less than 200 (Tables 3 and 4). The resulting timetree suggests an origin for oomycetes in the mid-Paleozoic, with a divergence between two major lineages, the peronosporaleans and saprolegnialeans, in the early Mesozoic (Figure 1). A complete list of divergence times with 95% confidence intervals for each node under each model is presented in Additional file 2.
Models for estimating divergence times under a molecular clock have become more complex over the past two decades. In this study we have used three distinct models, a single-rate strict clock, a UCLD relaxed clock, and a random local clock, to estimate divergence times among the fungal-like oomycetes. Analyses run under the strict clock model performed robustly, with all parameters showing evidence of thorough sampling (ESS > > 1000) and chain convergence. Because we had no a priori expectation of rate homogeneity among oomycetes or between oomycetes and ochrophytes, we also estimated divergence times under “relaxed” clock models. Both the UCLD and random local clock models indicated moderate to high levels of rate variation among lineages (as shown by the coefficient of variation parameter, Table 4), suggesting that a strict clock model was not appropriate for our dataset regardless of performance of the MCMC. In addition, an analysis of Bayes factors suggested that the two relaxed clock methods were a better fit for the data (ln Bayes factor in favor of relaxed clock models over strict clock >100). Rates estimated under the UCLD model appeared to be strongly influenced by the calibration priors, leading to rates 1.5 to 3.5 times higher in the ochrophyte lineages than in the oomycetes (data not shown). However, UCLD analyses failed to converge even after 50 million generations, thus limiting our ability to interpret parameter and divergence time estimates. Only a few parameters showed signs of poor mixing in the random local clock analyses (ESS < 200), but in general there was good evidence of chain convergence under this model, with the trade-off of long computational times.
Despite differences in performance among the three clock models, divergence time estimates among oomycetes were strikingly consistent (Table 3 and Additional file 2), and all models estimated a mid-Paleozoic origin for oomycetes (Figure 1). Our estimate for the divergence of oomycetes from other stramenopiles is somewhat consistent with results from a study of ochrophyte evolution using small subunit ribosomal DNA data , but is considerably younger than estimates generated from broader studies of eukaryote evolution [31, 32]. However, it seems likely that the times recovered here for the divergence between oomycetes and ochrophytes, as well as the root node, may be underestimated, for several reasons. A recent simulation study of relaxed clock models showed that the deepest nodes in a tree tend to be underestimated when shallow calibrations are used , which reflects our reliance on diatom calibrations to estimate divergences throughout the tree. Also, the posterior distributions recovered for the ingroup (node g in Figure 1) and root (node h) time estimates overlapped with their respective prior distributions, and were tightly constrained by the lower limit of 408 Ma imposed by the priors (data not shown). In addition, the long branch connecting the origin of oomycetes (node g) to the divergence between the peronosporaleans and saprolegnialeans (node d), as well as the long branch in the calibration taxa (between nodes e and f), may have influenced rate estimates under the UCLD and random local clock models. As a result, divergence times estimated for these nodes were sensitive to the model, particularly the ochrophyte estimates under the UCLD clock (Table 3 and Additional file 2); however, given the poor performance of the UCLD analysis, it is difficult to assess the reliability of these estimates. Additional sequence data from basal oomycetes such as Eurychasma dicksonii and Haliphthoros sp. , as well as from more ochrophyte calibration taxa, will help break up these long branches and led to more reliable rate estimates. The oldest accepted oomycete fossils come from the Lower Devonian Rhynie Chert, which is thought to have been a non-marine hot spring environment [21, 22]. Phylogenetic evidence suggests that the earliest diverging oomycetes were likely marine [15, 20], therefore the origin of this group may have occurred some time prior to the appearance of fossils in non-marine environments.
Fossil evidence of oomycetes also occurs throughout the Carboniferous, particularly in association with lycophytes (reviewed in ). While previous authors have suggested affinities with certain taxonomic groups (e.g., [25, 26]), the divergence times estimated here indicate that modern peronosporalean and saprolegnialean lineages originated much later, in the mid to late Mesozoic (Figure 1). Modern saprolegnialeans, such as Saprolegnia parasitica, are commonly associated with freshwater environments, and can be devastating pathogens of fish, amphibians, crustaceans, and insects ; saprotrophic species, such as Thraustotheca clavata, are also known from this group. In contrast, modern peronosporaleans are predominately terrestrial and many are significant plant pathogens. Two species included in our analysis, Hyaloperonospora arabidopsidis and Albugo laibachii, are obligate biotrophs who are fully dependent on their host (Arabidopsis). Phytophthora species cause disease on a wide variety of plants, and significant effort has been undertaken to understand their mechanisms of virulence and host specificity (reviewed in ). While it is undesirable to extrapolate as to the likely hosts for early diverging lineages, it does seem reasonable to suggest that host availability was not a constraining factor in oomycete diversification. Particularly for the modern plant pathogenic oomycetes, both fossil and molecular clock evidence suggests that the major lineages of angiosperms had diversified by the mid-Cretaceous , prior to our estimates for divergences among the peronosporaleans. The evolution of pathogenic lifestyles, therefore, may have been in response to certain environmental changes, or may have been facilitated by the horizontal transfer of pathogenicity-related genes from true Fungi [53–55] or from bacteria [45, 56], as has been suggested previously.
In this study, we chose to focus on conserved regulators of eukaryotic gene expression to examine their presence and level of conservation in pathogenic oomycetes. Mechanisms of gene expression regulation are highly conserved across eukaryotes and were most likely present in the last common ancestor, including epigenetic and RNA-based processes for transcriptional and post-transcriptional gene silencing [57–59]. Although we have not conducted an exhaustive survey here, our results suggest that the common ancestor of oomycetes possessed a full complement of regulatory proteins, including those involved in histone modification, RNA interference, and tRNA and rRNA methylation. Surprisingly, no orthologs of canonical DNA methyltransferases could be identified in the genomes of oomycetes. A single putative DNA methylase is present in the genome of Pythium ultimum (T014901), but no orthologs could be detected in the other oomycete genomes. Gene silencing studies in Phytophthora infestans have failed to detect evidence of cytosine methylation [60, 61], however recent work in P. sojae does suggest the presence of methylated DNA . DNA methyltransferases also appear to be absent from the Ectocarpus genome , as well as from the model eukaryotes Saccharomyces cerevisiae and Caenorhabditis elegans, however several are known from diatoms [64, 65]. Further study is therefore needed to confirm the presence and mechanism of DNA methylation in oomycetes.
This is the first study to estimate divergence times among the fungal-like oomycetes. The consistency of our time estimates under three distinct molecular clock models suggests that the resulting timetree likely recovers the main divergences among lineages, which occurred in the mid to late Mesozoic. Our estimates for the origin of oomycetes and the divergence of stramenopiles from other eukaryotes may have been underestimated due to the limited fossil information available for the taxa included in this study. Additional information from the oomycete fossil record, especially from the diverse Cretaceous assemblages, as well as new sequence data from basal oomycete lineages and other under-sampled eukaryotes , may help future molecular clock studies better estimate evolutionary rates.
Reference sequences for canonical eukaryotic transcription factors and proteins involved in post-transcriptional gene silencing, DNA and RNA methylation, and chromatin modification were obtained from the Gene Database at NCBI (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/gene) for human, Drosophila, Saccharomyces, and/or Arabidopsis. The reference protein sequences were then used to search for homologs in the genome of Phytophthora infestans T30-4 . Additional reference sequences were also obtained from a study of gene silencing in P. infestans. Both the eukaryotic reference sequences and the putative P. infestans homologs were used to search the available genomes of oomycetes, diatoms, and a brown alga (Table 1); outgroup sequences were obtained from Tetrahymena thermophila when available. All potential homologs of equivalent BLAST e-values within each genome were included for orthology assessment.
Protein domains were determined for all potential homologs using Pfam . Sequences that did not contain the appropriate domains for proper protein function were removed from each dataset except in cases where the protein sequence appeared truncated due to genome misannotation, particularly for Hyaloperonospora arabidopsidis. Each dataset was aligned under default settings in ClustalX v2 , and preliminary neighbor-joining phylogenies were generated under a Poisson correction with pairwise deletion of alignment gaps in MEGA v5 . Sequences within each dataset were considered orthologous if they shared protein domains and their phylogeny reflected known species relationships. In datasets with species-specific paralogs, one sequence was arbitrarily chosen to represent the ortholog for divergence time estimation. In cases where orthology was ambiguous or no homolog could be identified, the sequence was coded as missing data. A complete list of protein accession numbers per gene for each genome is available in Additional file 1.
Divergence time analysis
Protein datasets with robust orthology were used to co-estimate phylogeny and divergence times using Bayesian inference in BEAST v1.7.5 . Initial runs of 10 million generations were used under each clock model to evaluate settings on priors and to generate a user tree for subsequent analyses. For the final analyses, each protein dataset was treated as a separate partition under a WAG substitution model; a Yule speciation process was assumed with a uniform distribution on the birthrate (0–100; initial value 0.01). For the strict clock analyses, the rate parameter (clock.rate) was modeled with an exponential prior distribution (mean 1.0, initial value 0.01). For the UCLD relaxed clock model, an exponential prior distribution (mean 0.1, initial value 0.01) was used for the mean rate (ucld.mean) and standard deviation (ucld.stdev). Several parameters control the rate and number of rate changes under the random local clock model; a Poisson distribution (mean 0.693) was used as the prior for the number of local clocks (rateChanges), an exponential prior distribution (mean 1.0, initial value 0.001) was used for the relative rates among local clocks (localclocks.relativerates), and an exponential prior distribution (mean 1.0, initital value 0.01) was used for the rate (clock.rate). Five independent analyses were run for 50 million generations each, under all three clock models; log and tree files from the two runs with the highest parameter ESS values per model were combined (after removing burn-in from each run) using LogCombiner v1.7.5. Tracer v1.5  was used to evaluate convergence, estimate the appropriate burn-in for each run, and calculate Bayes factors for model comparisons. Analyses were also repeated without data (priors only) to determine the impact of calibration settings on the resulting divergence time estimates; three independent runs of 50 million generations each were performed under each clock model. Trees were visualized in FigTree v1.4 .
Fossil evidence from diatoms and oomycetes was used to calibrate the molecular clock analyses; all calibrations were modeled with a gamma prior distribution (shape 2.0) with the offset value set as the uppermost boundary of the time interval (stage) containing the relevant fossil. The value for the scale parameter was set so that the age at the 95% quantile was roughly equivalent to the lowermost boundary on the geologic epoch containing the relevant fossil. Appropriate geological times were obtained from the International Commission on Stratigraphy chronostratigraphic chart, January 2013 version (http://stratigraphy.org). Fossil evidence from the Late Cretaceous (Campanian) pennate diatoms  provided a minimum age of 72.1 Ma on the divergence between Thalassiosira and Phaeodactylum (5-95% quantiles = 74–100 Ma). Early Jurassic (Toarcian) diatom fossils  provided a minimum age of 174 Ma on the divergence between diatoms and Ectocarpus (5-95% quantiles = 176–202 Ma). The Rhynie chert oomycete fossils  were used to define a minimum divergence time of 408 Ma between oomycetes and ochrophytes (5-95% quantiles = 418–550 Ma). A wide uniform prior distribution (408–1750 Ma; initial value 635 Ma) was used for the root age as there are few robust estimates on the divergence between alveolates and stramenopiles. Beast XML-formatted data files have been deposited in Dryad .
Availability of supporting data
The data sets supporting the results of this article are available in the Dryad Digital Repository, http://0-dx.doi.org.brum.beds.ac.uk/10.5061/dryad.39mc5.
Knoll AH: The fossil record of microbial life. Fundamentals of Geobiology. Edited by: Knoll AH, Canfield DE, Konhauser KO. 2012, Oxford: Wiley-Blackwell
Zuckerkandl E, Pauling L: Molecular disease, evolution, and genic heterogeneity. Horizons in Biochemistry. Edited by: Kasha M, Pullman B. 1962, New York: Academic Press, 189-225.
Zuckerkandl E, Pauling L: Evolutionary divergence and convergence in proteins. Evolving Genes and Proteins. Edited by: Bryson V, Vogel HJ. 1965, New York: Academic Press, 97-166.
Kumar S: Molecular clocks: four decades of evolution. Nat Rev Genet. 2005, 6 (8): 654-662. 10.1038/nrg1659.
Sanderson MJ: A nonparametric approach to estimating divergence times in the absence of rate constancy. Mol Biol Evol. 1997, 14 (12): 1218-1231. 10.1093/oxfordjournals.molbev.a025731.
Sanderson MJ: Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach. Mol Biol Evol. 2002, 19 (1): 101-109. 10.1093/oxfordjournals.molbev.a003974.
Thorne JL, Kishino H, Painter IS: Estimating the rate of evolution of the rate of molecular evolution. Mol Biol Evol. 1998, 15 (12): 1647-1657. 10.1093/oxfordjournals.molbev.a025892.
Aris-Brosou S, Yang Z: Effects of models of rate evolution on estimation of divergence dates with special reference to the metazoan 18S ribosomal RNA phylogeny. Syst Biol. 2002, 51 (5): 703-714. 10.1080/10635150290102375.
Drummond AJ, Ho SYW, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4 (5): e88-10.1371/journal.pbio.0040088.
Drummond AJ, Suchard MA: Bayesian random local clocks, or one rate to rule them all. BMC Biol. 2010, 8: 114-10.1186/1741-7007-8-114.
Ho SYW, Phillips MJ: Accounting for calibration uncertainty in phylogenetic estimation of evolutionary divergence times. Syst Biol. 2009, 58 (3): 367-380. 10.1093/sysbio/syp035.
Parham JF, Donoghue PCJ, Bell CJ, Calway TD, Head JJ, Holroyd PA, Inoue JG, Irmis RB, Joyce WG, Ksepka DT, Patané JSL, Smith ND, Tarver JE, van Tuinen M, Yang Z, Angielczyk KD, Greenwood JM, Hipsley CA, Jacobs L, Makovicky PJ, Müller J, Smith KT, Theodor JM, Warnock RCM, Benton MJ: Best practices for justifying fossil calibrations. Syst Biol. 2012, 61 (2): 346-359. 10.1093/sysbio/syr107.
Warnock RCM, Yang Z, Donoghue PCJ: Exploring uncertainty in the calibration of the molecular clock. Biol Lett. 2012, 8 (1): 156-159. 10.1098/rsbl.2011.0710.
Dick MW: Straminipilous Fungi. 2001, Dordrecht: Kluwer Academic Publishers
Beakes GW, Glocklin SL, Sekimoto S: The evolutionary phylogeny of the oomycete "fungi". Protoplasma. 2012, 249 (1): 3-19. 10.1007/s00709-011-0269-2.
Cavalier-Smith T: Principles of protein and lipid targeting in secondary symbiogenesis: Euglenoid, Dinoflagellate, and Sporozoan plastid origins and the eukaryote family tree. J Eukaryot Microbiol. 1999, 46 (4): 347-366. 10.1111/j.1550-7408.1999.tb04614.x.
Burki F, Shalchian-Tabrizi K, Minge M, Skjæveland Å, Nikolaev SI, Jakobsen KS, Pawlowski J: Phylogenomics reshuffles the eukaryotic supergroups. PLoS ONE. 2007, 2 (8): e790-10.1371/journal.pone.0000790.
Adl SM, Simpson AGB, Lane CE, Lukeš J, Bass D, Bowser SS, Brown MW, Burki F, Dunthorn M, Hampl V, Heiss A, Hoppenrath M, Lara E, le Gall L, Lynn DH, McManus H, Mitchell EAD, Mozley-Stanridge SE, Parfrey LW, Pawlowski J, Rueckert S, Shadwick L, Schoch CL, Smirnov A, Spiegel FW: The Revised Classification of Eukaryotes. J Eukaryot Microbiol. 2012, 59 (5): 429-514. 10.1111/j.1550-7408.2012.00644.x.
Fry W: Phytophthora infestans: the plant (and R gene) destroyer. Mol Plant Pathol. 2008, 9 (3): 385-402. 10.1111/j.1364-3703.2007.00465.x.
Beakes GW, Sekimoto S: The evolutionary phylogeny of oomycetes - insights gained from studies of holocarpic parasites of algae and invertebrates. Oomycete Genetics and Genomics: Diversity, Interactions, and Research Tools. Edited by: Lamour K, Kamoun S. 2009, Hoboken, NJ: John Wiley & Sons, Inc, 1-24.
Krings M, Taylor TN, Dotzler N: The fossil record of the Peronosporomycetes (Oomycota). Mycologia. 2011, 103 (3): 445-457. 10.3852/10-278.
Taylor TN, Krings M, Kerp H: Hassiella monospora gen. et sp. nov., a microfungus from the 400 million year old Rhynie chert. Mycol Res. 2006, 110 (6): 628-632. 10.1016/j.mycres.2006.02.009.
Krings M, Taylor TN, Taylor EL, Kerp H, Hass H, Dotzler N, Harper CJ: Microfossils from the Lower Devonian Rhynie Chert with Suggested Affinities to the Peronosporomycetes. J Paleontol. 2012, 86 (2): 358-367. 10.1666/11-087.1.
Krings M, Taylor TN, Dotzler N, Decombeix A-L: Galtierella biscalithecae nov. gen. et sp., a Late Pennsylvanian endophytic water mold (Peronosporomycetes) from France. Comptes Rendus Palevol. 2010, 9 (1–2): 5-11.
Strullu-Derrien C, Kenrick P, Rioult JP, Strullu DG: Evidence of parasitic Oomycetes (Peronosporomycetes) infecting the stem cortex of the Carboniferous seed fern Lyginopteris oldhamia. Proc R Soc B Biol Sci. 2011, 278 (1706): 675-680. 10.1098/rspb.2010.1603.
Stidd BM, Cosentino K: Albugo-like oogonia from the American Carboniferous. Science. 1975, 190 (4219): 1092-1093. 10.1126/science.190.4219.1092.
Dotzler N, Krings M, Agerer R, Galtier J, Taylor TN: Combresomyces cornifer gen. sp. nov., an endophytic peronosporomycete in Lepidodendron from the Carboniferous of central France. Mycol Res. 2008, 112 (9): 1107-1114. 10.1016/j.mycres.2008.03.003.
Schwendemann AB, Taylor TN, Taylor EL, Krings M, Dotzler N: Combresomyces cornifer from the Triassic of Antarctica: Evolutionary stasis in the Peronosporomycetes. Rev Palaeobot Palynol. 2009, 154 (1–4): 1-5.
Krings M, Taylor TN, Dotzler N: Fungal endophytes as a driving force in lan plant evolution: evidence from the fossil record. Biocomplexity of Plant-Fungal Interactions. Edited by: Southworth D. 2012, Ames, Iowa: John Wiley & Sons, Inc
Berney C, Pawlowski J: A molecular time-scale for eukaryote evolution recalibrated with the continuous microfossil record. Proc R Soc B Biol Sci. 2006, 273 (1596): 1867-1872. 10.1098/rspb.2006.3537.
Parfrey LW, Lahr DJG, Knoll AH, Katz LA: Estimating the timing of early eukaryotic diversification with multigene molecular clocks. Proc Natl Acad Sci. 2011, 108 (33): 13624-13629. 10.1073/pnas.1110633108.
Hackett JD, Yoon HS, Butterfield NJ, Sanderson MJ, Bhattacharya D, Falkowski PG, Knoll AH: Plastid endosymbiosis: sources and timing of the major events. Evolution of Primary Producers in the Sea. 2007, Burlington, MA: Elsevier Academic Press, 109-132.
Phillips N, Calhoun S, Moustafa A, Bhattacharya D, Braun EL: Genomic insights into evolutionary relationships among heterokont lineages emphasizing the Phaeophyceae. J Phycol. 2008, 44 (1): 15-18. 10.1111/j.1529-8817.2007.00435.x.
Brown JW, Sorhannus U: A molecular genetic timescale for the diversification of autotrophic stramenopiles (Ochrophyta): substantive underestimation of putative fossil ages. PLoS ONE. 2010, 5 (9): e12759-10.1371/journal.pone.0012759.
Siomi H, Siomi MC: On the road to reading the RNA-interference code. Nature. 2009, 457 (7228): 396-404. 10.1038/nature07754.
Vetukuri RR, Avrova AO, Grenville-Briggs LJ, Van West P, Soderbom F, Savenkov EI, Whisson SC, Dixelius C: Evidence for involvement of Dicer-like, Argonaute and histone deacetylase proteins in gene silencing in Phytophthora infestans. Mol Plant Pathol. 2011, 12 (8): 772-785. 10.1111/j.1364-3703.2011.00710.x.
Kemen E, Gardiner A, Schultz-Larsen T, Kemen AC, Balmuth AL, Robert-Seilaniantz A, Bailey K, Holub E, Studholme DJ, MacLean D, Jones JDG: Gene gain and loss during evolution of obligate parasitism in the white rust pathogen of Arabidopsis thaliana. PLoS Biol. 2011, 9 (7): e1001094-10.1371/journal.pbio.1001094.
Cock JM, Sterck L, Rouze P, Scornet D, Allen AE, Amoutzias G, Anthouard V, Artiguenave F, Aury J-M, Badger JH, Beszteri B, Billiau K, Bonnet E, Bothwell JH, Bowler C, Boyen C, Brownlee C, Carrano CJ, Charrier B, Cho GY, Coelho SM, Collen J, Corre E, Da Silva C, Delage L, Delaroque N, Dittami SM, Doulbeau S, Elias M, Farnham G, et al: The Ectocarpus genome and the independent evolution of multicellularity in brown algae. Nature. 2010, 465 (7298): 617-621. 10.1038/nature09016.
Baxter L, Tripathy S, Ishaque N, Boot N, Cabral A, Kemen E, Thines M, Ah-Fong A, Anderson R, Badejoko W, Bittner-Eddy P, Boore JL, Chibucos MC, Coates M, Dehal P, Delehaunty K, Dong S, Downton P, Dumas B, Fabro G, Fronick C, Fuerstenberg SI, Fulton L, Gaulin E, Govers F, Hughes L, Humphray S, Jiang RHY, Judelson H, Kamoun S: Signatures of adaptation to obligate biotrophy in the Hyaloperonospora arabidopsidis genome. Science. 2010, 330 (6010): 1549-1551. 10.1126/science.1195203.
Bowler C, Allen AE, Badger JH, Grimwood J, Jabbari K, Kuo A, Maheswari U, Martens C, Maumus F, Otillar RP, Rayko E, Salamov A, Vandepoele K, Beszteri B, Gruber A, Heijde M, Katinka M, Mock T, Valentin K, Verret F, Berges JA, Brownlee C, Cadoret J-P, Chiovitti A, Choi CJ, Coesel S, De Martino A, Detter JC, Durkin C, Falciatore A, et al: The Phaeodactylum genome reveals the evolutionary history of diatom genomes. Nature. 2008, 456 (7219): 239-244. 10.1038/nature07410.
Lamour KH, Mudge J, Gobena D, Hurtado-Gonzales OP, Schmutz J, Kuo A, Miller NA, Rice BJ, Raffaele S, Cano LM, Bharti AK, Donahoo RS, Finley S, Huitema E, Hulvey J, Platt D, Salamov A, Savidor A, Sharma R, Stam R, Storey D, Thines M, Win J, Haas BJ, Dinwiddie DL, Jenkins J, Knight JR, Affourtit JP, Han CS, Chertkov O: Genome sequencing and mapping reveal loss of heterozygosity as a mechanism for rapid adaptation in the vegetable pathogen Phytophthora capsici. Mol Plant Microbe Interact. 2012, 25 (10): 1350-1360. 10.1094/MPMI-02-12-0028-R.
Haas BJ, Kamoun S, Zody MC, Jiang RHY, Handsaker RE, Cano LM, Grabherr M, Kodira CD, Raffaele S, Torto-Alalibo T, Bozkurt TO, Ah-Fong AMV, Alvarado L, Anderson VL, Armstrong MR, Avrova A, Baxter L, Beynon J, Boevink PC, Bollmann SR, Bos JIB, Bulone V, Cai G, Cakir C, Carrington JC, Chawner M, Conti L, Costanzo S, Ewan R, Fahlgren N, et al: Genome sequence and analysis of the Irish potato famine pathogen Phytophthora infestans. Nature. 2009, 461 (7262): 393-398. 10.1038/nature08358.
Tyler BM, Tripathy S, Zhang X, Dehal P, Jiang RHY, Aerts A, Arredondo FD, Baxter L, Bensasson D, Beynon JL, Tyler BM, Tripathy S, Zhang X, Dehal P, Jiang RHY, Aerts A, Arredondo FD, Baxter L, Bensasson D, Beynon JL, Chapman J, Damasceno CMB, Dorrance AE, Dou D, Dickerman AW, Dubchak IL, Garbelotto M, Gijzen M, Gordon SG, Govers F, et al: Phytophthora genome sequences uncover evolutionary origins and mechanisms of pathogenesis. Science. 2006, 313 (5791): 1261-1266. 10.1126/science.1128796.
Levesque CA, Brouwer H, Cano L, Hamilton J, Holt C, Huitema E, Raffaele S, Robideau G, Thines M, Win J, Levesque CA, Brouwer H, Cano L, Hamilton J, Holt C, Huitema E, Raffaele S, Robideau G, Thines M, Win J, Zerillo M, Beakes G, Boore J, Busam D, Dumas B, Ferriera S, Fuerstenberg S, Gachon C, Gaulin E, Govers F, et al: Genome sequence of the necrotrophic plant pathogen Pythium ultimum reveals original pathogenicity mechanisms and effector repertoire. Genome Biol. 2010, 11 (7): R73-10.1186/gb-2010-11-7-r73.
Jiang RHY, de Bruijn I, Haas BJ, Belmonte R, Löbach L, Christie J, van den Ackerveken G, Bottin A, Bulone V, Díaz-Moreno SM, Dumas B, Fan L, Gaulin E, Govers F, Grenville-Briggs LJ, Horner NR, Levin JZ, Mammella M, Meijer HJG, Morris P, Nusbaum C, Oome S, Phillips AJ, van Rooyen D, Rzeszutek E, Saraiva M, Secombes CJ, Seidl MF, Snel B, Stassen JHM, et al: Distinctive expansion of potential virulence genes in the genome of the oomycete fish pathogen Saprolegnia parasitica. PLoS Genet. 2013, 9 (6): e1003272-10.1371/journal.pgen.1003272.
Stover NA, Krieger CJ, Binkley G, Dong Q, Fisk DG, Nash R, Sethuraman A, Weng S, Cherry JM: Tetrahymena Genome Database (TGD): a new genomic resource for Tetrahymena thermophila research. Nucleic Acids Res. 2006, 34 (S1): D500-D503.
Armbrust EV, Berges JA, Bowler C, Green BR, Martinez D, Putnam NH, Zhou S, Allen AE, Apt KE, Bechner M, Brzezinski MA, Chaal BK, Chiovitti A, Davis AK, Demarest MS, Detter JC, Glavina T, Goodstein D, Hadi MZ, Hellsten U, Hildebrand M, Jenkins BD, Jurka J, Kapitonov VV, Kröger N, Lau WWY, Lane TW, Larimer FW, Lippmeier JC, Lucas S, et al: The genome of the diatom Thalassiosira pseudonana: ecology, evolution, and metabolism. Science. 2004, 306 (5693): 79-86. 10.1126/science.1101156.
Battistuzzi FU, Filipski A, Hedges SB, Kumar S: Performance of relaxed-clock methods in estimating evolutionary divergence times and their credibility intervals. Mol Biol Evol. 2010, 27 (6): 1289-1300. 10.1093/molbev/msq014.
Grenville-Briggs L, Gachon CMM, Strittmatter M, Sterck L, Küpper FC, van West P: A Molecular Insight into Algal-Oomycete Warfare: cDNA Analysis of Ectocarpus siliculosus Infected with the Basal Oomycete Eurychasma dicksonii. PLoS ONE. 2011, 6 (9): e24500-10.1371/journal.pone.0024500.
Sekimoto S, Hatai K, Honda D: Molecular phylogeny of an unidentified Haliphthoros-like marine oomycete and Haliphthoros milfordensis inferred from nuclear-encoded small- and large-subunit rRNA genes and mitochondrial-encoded cox2 gene. Mycoscience. 2007, 48 (4): 212-221. 10.1007/S10267-007-0357-5.
Jiang RHY, Tyler BM: Mechanisms and evolution of virulence in oomycetes. Annu Rev Phytopathol. 2012, 50 (1): 295-318. 10.1146/annurev-phyto-081211-172912.
Bell CD, Soltis DE, Soltis PS: The age and diversification of the angiosperms re-revisited. Am J Bot. 2010, 97 (8): 1296-1303. 10.3732/ajb.0900346.
Luis P, Gauthier A, Trouvelot S, Poinssot B, Frettinger P: Identification of Plasmopara viticola genes potentially involved in pathogenesis on grapevine suggests new similarities between oomycetes and true Fungi. Phytopathology. 2013, 103 (10): 1035-1044. 10.1094/PHYTO-06-12-0121-R.
Morris PF, Schlosser LR, Onasch KD, Wittenschlaeger T, Austin R, Provart N: Multiple horizontal gene transfer events and domain fusions have created novel regulatory and metabolic networks in the oomycete genome. PLoS ONE. 2009, 4 (7): e6133-10.1371/journal.pone.0006133.
Richards TA, Soanes DM, Jones MDM, Vasieva O, Leonard G, Paszkiewicz K, Foster PG, Hall N, Talbot NJ: Horizontal gene transfer facilitated the evolution of plant parasitic mechanisms in the oomycetes. Proc Natl Acad Sci. 2011, 108 (37): 15258-15263. 10.1073/pnas.1105100108.
Belbahri L, Calmin G, Mauch F, Andersson JO: Evolution of the cutinase gene family: evidence for lateral gene transfer of a candidate Phytophthora virulence factor. Gene. 2008, 408 (1–2): 1-8.
Collins L, Chen XS: Ancestral RNA: The RNA biology of the eukaryotic ancestor. RNA Biol. 2009, 6 (5): 495-502. 10.4161/rna.6.5.9551.
Katz LA: Origin and diversification of eukaryotes. Annu Rev Microbiol. 2012, 66 (1): 411-427. 10.1146/annurev-micro-090110-102808.
Shabalina SA, Koonin EV: Origins and evolution of eukaryotic RNA interference. Trends Ecol Evol. 2008, 23 (10): 578-587. 10.1016/j.tree.2008.06.005.
van West P, Kamoun S, van’t Klooster JW, Govers F: Internuclear gene silencing in Phytophthora infestans. Mol Cell. 1999, 3 (3): 339-348. 10.1016/S1097-2765(00)80461-X.
van West P, Shepherd SJ, Walker CA, Li S, Appiah AA, Grenville-Briggs LJ, Govers F, Gow NAR: Internuclear gene silencing in Phytophthora infestans is established through chromatin remodelling. Microbiology. 2008, 154 (5): 1482-1490. 10.1099/mic.0.2007/015545-0.
Spangler M: Cytosine methylation of Phytophthora sojae by methylated DNA immunoprecipitation. MS thesis. 2012, Biological Sciences Department: Bowling Green State Universit
Goll MG, Bestor TH: Eukaryotic cytosine methyltransferases. Annu Rev Biochem. 2005, 74 (1): 481-514. 10.1146/annurev.biochem.74.010904.153721.
De Riso V, Raniello R, Maumus F, Rogato A, Bowler C, Falciatore A: Gene silencing in the marine diatom Phaeodactylum tricornutum. Nucleic Acids Res. 2009, 37 (14): e96-10.1093/nar/gkp448.
Montsant A, Allen AE, Coesel S, Martino AD, Falciatore A, Mangogna M, Siaut M, Heijde M, Jabbari K, Maheswari U, Rayko E, Vardi A, Apt KE, Berges JA, Chiovitti A, Davis AK, Thamatrakoln K, Hadi MZ, Lane TW, Lippmeier JC, Martinez D, Parker MS, Pazour GJ, Saito MA, Rokhsar DS, Armbrust EV, Bowler C: Identification and comparative genomic analysis of signaling and regulatory components in the diatom Thalassiosira pseudonana. J Phycol. 2007, 43 (3): 585-604. 10.1111/j.1529-8817.2007.00342.x.
del Campo J, Sieracki ME, Molestina R, Keeling P, Massana R, Ruiz-Trillo I: The others: our biased perspectives of eukaryotic genomes. Trends Ecol Evol. 2014, 29 (5): 252-259. 10.1016/j.tree.2014.03.006.
Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, Heger A, Holm L, Sonnhammer ELL, Eddy SR, Bateman A, Finn RD: The Pfam protein families database. Nucleic Acids Res. 2012, 40 (D1): D290-D301. 10.1093/nar/gkr1065.
Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23 (21): 2947-2948. 10.1093/bioinformatics/btm404.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.
Drummond AJ, Suchard MA, Xie D, Rambaut A: Bayesian Phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012, 29 (8): 1969-1973. 10.1093/molbev/mss075.
Rambaut A, Drummond AJ: Tracer version 1.5, available at. http://tree.bio.ed.ac.uk/software/tracer/,
Rambaut A, Drummond AJ: FigTree version 1.4, available at. http://tree.bio.ed.ac.uk/software/figtree,
Kooistra W, Gersonde R, Medlin LK, Mann DG: The origin and evolution of the Diatoms: their adaptation to a planktonic existence. Evolution of Primary Producers in the Sea. Edited by: Falkowski PG, Knoll AH. 2007, Burlington, MA: Elsevier Academic Press, 207-249.
Sims PA, Mann DG, Medlin LK: Evolution of the diatoms: insights from fossil, biological and molecular data. Phycologia. 2006, 45 (4): 361-402. 10.2216/05-22.1.
Matari N, Blair JE: Data from: a multilocus timescale for oomycete evolution estimated under three distinct molecular clock models. Dryad Digital Repository. 2014, http://0-dx.doi.org.brum.beds.ac.uk/10.5061/dryad.39mc5,
The authors would like to thank Dr. Chris Lane (University of Rhode Island) for access to unpublished genomic data from Achlya hypogyna and Thraustotheca clavata, and Dr. Vipaporn Phuntumart (Bowling Green State University) for discussion of unpublished methylation data from Phytophthora sojae. We also thank Jason Brooks and Anthony Weaver for assistance with the random local clock analyses, and Dr. Jorge Mena-Ali (F&M) for helpful conversations regarding parameter settings. This work was supported by the US Department of Agriculture National Institute of Food and Agriculture (2011-68004-30104 and 2010-65110-20488 to J.E.B) and by a grant to Franklin & Marshall College from the Howard Hughes Medical Institute Science Education Program.
The authors declare that they have no competing interests.
JEB and NHM designed the experiment and performed the data mining. JEB performed the divergence time analyses and drafted the manuscript. Both authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Matari, N.H., Blair, J.E. A multilocus timescale for oomycete evolution estimated under three distinct molecular clock models. BMC Evol Biol 14, 101 (2014) doi:10.1186/1471-2148-14-101
- Divergence times
- Bayesian inference
- Molecular clock
- Gene expression regulation