Cross-kingdom host shifts of phytomyxid parasites
BMC Evolutionary Biology volume 14, Article number: 33 (2014)
Phytomyxids (plasmodiophorids and phagomyxids) are cosmopolitan, obligate biotrophic protist parasites of plants, diatoms, oomycetes and brown algae. Plasmodiophorids are best known as pathogens or vectors for viruses of arable crops (e.g. clubroot in brassicas, powdery potato scab, and rhizomania in sugar beet). Some phytomyxid parasites are of considerable economic and ecologic importance globally, and their hosts include important species in marine and terrestrial environments. However most phytomyxid diversity remains uncharacterised and knowledge of their relationships with host taxa is very fragmentary.
Our molecular and morphological analyses of phytomyxid isolates–including for the first time oomycete and sea-grass parasites–demonstrate two cross-kingdom host shifts between closely related parasite species: between angiosperms and oomycetes, and from diatoms/brown algae to angiosperms. Switching between such phylogenetically distant hosts is generally unknown in host-dependent eukaryote parasites. We reveal novel plasmodiophorid lineages in soils, suggesting a much higher diversity than previously known, and also present the most comprehensive phytomyxid phylogeny to date.
Such large-scale host shifts between closely related obligate biotrophic eukaryote parasites is to our knowledge unique to phytomyxids. Phytomyxids may readily adapt to a wide diversity of new hosts because they have retained the ability to covertly infect alternative hosts. A high cryptic diversity and ubiquitous distribution in agricultural and natural habitats implies that in a changing environment phytomyxids could threaten the productivity of key species in marine and terrestrial environments alike via host shift speciation.
Phytomyxea (phytomyxids) are a poorly known group of obligate biotrophic, endobiotic parasites of plants, diatoms, brown algae, and oomycetes [1, 2]. Their phylogenetic position was long debated but molecular phylogenies now robustly place them within the eukaryote supergroup Rhizaria, as sister group to the omnivorous vampyrellid amoebae [3, 4]. The group is subdivided into two orders: Phagomyxida (phagomyxids) and Plasmodiophorida (plasmodiophorids; [5, 6]). Phagomyxids are very poorly known with only four characterised lineages, comprising two pairs of closely related species in the genera Phagomyxa (parasites of marine diatoms; ) and Maullinia (parasites of brown algae; [6, 8]). In addition to these described species, several morphologically uncharacterised DNA-lineages have been detected by environmental sequencing [9, 10].
Plasmodiophorids are better known because they include plant parasites causing significant diseases of crops including brassicas, potatoes, and grain crops (e.g. maize, rice, wheat, sorghum). The best studied species is the clubroot-causing Plasmodiophora brassicae, a parasite of crucifers which accounts for up to 10% loss of the worldwide production of Brassica crops . Other well-studied species include Spongospora subterranea, which causes powdery scab of potato and can serve as a vector for Potato Mop Top Virus . The otherwise symptomless Polymyxa graminis transmits economically important viruses to a number of grain crops  while Polymyxa betae is the vector for beet necrotic yellow vein virus, the cause of sugar beet “rhizomania” . Given the impact of plasmodiophorids on staple crops it is not surprising that research has been focussed on the species causing economic damage, while research on species from outside of agricultural environments has been rare.
Phytomyxids are biotrophic parasites, i.e. they infect and multiply in a host without killing it, as opposed to hemi-biotrophic or nectrotrophic parasites, which ultimately kill the host . Upon infection phytomyxids use an arsenal of methods to manipulate host physiology and to escape defence responses [16, 17]. A dependence on living hosts to survive and to complete their life cycle is invariably associated with high host specificity because of the highly specialised interactions involved [16, 17].
Phytomyxids have a complex life cycle consisting of two functionally different zoosporic stages (primary and secondary), which upon infection of a suitable host develop into multinucleate plasmodia inside the infected host cell (for a detailed review of the life cycle see e.g. ). These plasmodia either develop into a zoosporangium which rapidly releases large numbers of secondary zoospores or into resting spores which are the survival stage of the parasite and from which primary zoospores are released [1, 18, 19]. The resting spores can be aggregated into so-called sporosori and this arrangement is characteristic for each genus. Because of their complex six-stage life cycle [1, 18], their small size (3-6 μm; ), and the fact that some life cycle stages are difficult to find and identify morphologically, there have been no recent studies of plasmodiophorid biodiversity; previous studies date back to the pre-molecular era [1, 20, 21]. Phytomyxids are not well represented in DNA databases: sequences are only available for nine of the 41 described species. An ecologically and evolutionary important group missing from molecular phylogenies are the phytomyxean oomycete parasites. These were originally placed in the now superseded taxon Plasmodiophoromycetes, then still included within Fungi . The taxonomy of one genus of oomycete parasites–Woronina spp.–was studied two decades ago using ultrastructural studies, but Woronina has not previously been included in molecular phylogenies including more recently discovered parasites of other stramenopile hosts [7, 8].
Molecular biodiversity studies are important for understanding the evolution, abundance and the possible economic and ecological impact of parasites, especially as global parasite numbers are thought to be hugely underestimated . This unknown parasite biodiversity may have a large impact in a changing environment, where the emergence of existing and ‘novel’ diseases can have devastating effects [23–25]. Therefore, a targeted study of phytomyxid biodiversity is overdue, because of their known potential to damage crops directly by causing diseases or indirectly by transmitting viruses [2, 15].
To further investigate phytomyxid biodiversity we used a parallel approach to screen selected soil ecosystems: microscopically examining potential hosts, coupled with targeted molecular screening of rhizosphere-associated and root-free soil from the same sites. Soil samples were obtained from sites ranging from intensively farmed to those with restricted anthropogenic input; we hypothesised that sites harbouring different host communities will also differ in parasite biodiversity. We also aimed to test whether phytomyxids that parasitise oomycetes group with phagomyxids or plasmodiophorids using a combination of isolation, microscopic identification, and 18S rDNA clone libraries. Through analysing 18S rDNA sequences we identified cross kingdom host shifts within Plasmodiophorida and Phagomyxida.
New phytomyxid isolates and phylogenetic analyses
18S rDNA sequences were generated for 20 phytomyxid isolates for which resting spores could be characterised by microscopy (Figure 1, Additional file 1: Table S1). These isolates were identified following morphotaxonomic concepts (Figure 2; ). Resting spore size, shape, and arrangement are the primary defining taxonomic characters of phytomyxids, along with the identity of their host species [1, 26]. The new sequences were included in Bayesian and Maximum Likelihood phylogenetic analyses with phytomyxid sequences from previously characterised specimens, phylogenetically informative environmental sequences [9, 10, 27], and the closely related environmental lineage ‘Novel Clade 9’ . Our new isolates represent six previously described species as well as two undescribed lineages (Maullinia sp., Figure 2e; Polymyxa cf. graminis 2; Additional file 1: Table S1). We provide the first 18S rDNA sequences for Sorosphaerula viticola (four isolates from the vine Vitis spp., Figure 2a), Ligniera junci (two isolates from the rush Juncus spp., Figure 2d), Woronina pythii from the oomycete Pythium sp. (Figure 2b, c), and Plasmodiophora diplantherae from the angiosperm sea grass Halodule wrightii (Figure 2f). We also sequenced five isolates of Polymyxa graminis from the grass Poa sp., two of Plasmodiophora brassicae from the crucifers Sinapis sp. and Brassica sp., and a plasmodiophorid we identified as P. graminis from Poa sp. (referred to as Polymyxa cf. graminis 2) which formed a distinct clade branching basally to P. graminis, S. veronicae, P. betae, and S. viticola in the 18S phylogeny. The two orders Phagomyxida and Plasmodiophorida, as previously conceived, were robustly supported. Seven distinct lineages could be resolved within the Phagomyxida and eleven within the Plasmodiophorida. According to this analysis, all plasmodiophorid genera (as currently circumscribed) that are represented in our phylogeny by more than one species are polyphyletic (Figure 1).
Cross-kingdom host shifts of closely related Phytomyxids
We show that plasmodiophorids infecting angiosperms are paraphyletic with respect to those infecting oomycetes, as are stramenopile-associated Phagomyxida to the angiosperm parasite Plasmodiophora diplantherae. P. diplantherae branches robustly within phagomyxids, sister to two environmental sequences which were found in clone libraries of deep-sea/anoxic saline sediments [9, 10]. The oomycete parasite Woronina pythii grouped within plasmodiophorids, closely related to Spongospora nasturtii and Plasmodiophora brassicae, both parasites of cruciferous plants (Figure 3). However, the branching order of these lineages within Plasmodiophorida is not resolved; Woronina may in fact be the most basally branching plasmodiophorid lineage. Light microscopy of W. pythii showed that size, arrangement and development of the resting spores from sporogenic plasmodia were very similar to that seen in other plasmodiophorids (e.g. P. brassicae, S. veronicae, S. viticola) (Figure 2B, C). The resting spores of P. diplantherae (Figure 2F) on the other hand strongly resembled those of the phagomyxid Maullinia sp. (Figure 2E) in ornamentation and a loose aggregation at maturity.
Furthermore, plasmodiophorids are known to rarely infect a range of so-called alternative hosts in which very limited (i.e. to single host cells) infections can be observed. Therefore, we conducted a comprehensive literature study to identify and summarise possible alternative hosts (Figure 4, Additional file 1: Table S2). Interactions with alternative hosts remain highly understudied; however one characteristic common to those known is that resting spores are rarely or never formed. According to these data some plasmodiophorid species can infect host plants belonging to two or more plant families; Ligniera junci, Polymyxa graminis, Plasmodiophora brassicae, and Spongospora subterranea are known to parasitize both monocot and eudicot host species (Figure 4). Plasmodiophorid parasites of oomycetes and the phagomyxid M. ectocarpii (Figure 2G) can infect a range of distantly related oomycetes or brown algae respectively (Additional file 1: Table S2).
Despite an extensive sampling of living plants in this study, no alternative host associations in nature could be confirmed by microscopic observation. However, direct and clone sequencing of rhizosphere samples revealed phytomyxid sequences associated with possible alternative host plants (Additional file 1: Table S3). Sequences closely related to Polymyxa graminis 1 (primary hosts Poaceae) were found associated with Taraxacum sp., Veronica sp., and Vitis sp.; sequences with a 96.4-97.2% similarity to Polymyxa betae (primary hosts Amaranthaceae) with Vitis sp.; a Spongospora nasturtii-like sequence (96.1% similarity to AF245217; primary host Nasturtium spp., Brassicacae) from Carex sp.; and S. subterranea-like sequences (98.4% similarity, primary hosts Solanaceae) from Trifolium sp. The nature of these associations remains unclear; microscopy of the associated root samples was not conclusive and sequence data alone does not provide sufficient information about organismal associations.
The division of Phytomyxea into Phagomyxida and Plasmodiophorida was confirmed by our phylogenetic (Figure 1) and morphological analyses, and also corresponded to the environment in which the host species were found rather than host taxonomy: plasmodiophorids are a monophyletic lineage of parasites of land plants and soil-dwelling oomycetes, while phagomyxids are a monophyletic lineage infecting marine plants and stramenopiles. We tested our plasmodiophorid-specific primer combination on phagomyxid isolates (Maullinia sp., Maullina ectocarpii, P. diplantherae), and on a selection of marine algae (brown, red, green) and marine (littoral) sediment from the South coast of the UK (n?=?38). All results were negative: no plasmodiophorids were detected in marine environments.
Environmental sequence analysis of plasmodiophorids
Plasmodiophorid-specific PCRs were carried out on 385 individual rhizosphere and soil samples (Additional file 1: Table S4), of which 153 gave positive results. After pooling of amplicons for cloning, a total of 303 18S rDNA reads were obtained from 21 clone libraries and direct sequencing of rhizosphere samples (n?=?19). After quality and chimera checks 295 of these sequences were analysed further. A total of 81 different plasmodiophorid sequence-types (18S-types, defined as described in the Methods) were detected (Additional file 1: Table S1, Additional file 1: Table S3). These grouped into 16 phylogenetically well supported clades within the Plasmodiophorida (Figure 5). Ten clades included sequences from at least one microscopically confirmed isolate, whereas six comprised only environmental 18S-types. Clade 4 (S. veronicae) and Clade 8 (environmental) are monotypic clades. The environmental lineages Clade 2, Clade 11, Clade 12, Clade 14 each consist of two 18S-types. Clade 3 (incl. S. viticola), Clade 7 (incl. L. junci) and Clade 16 (environmental) comprise three 18S-types each, while Clade 6 (incl. P. cf. graminis 2) consists of four 18S-types. Clade 10 (all P. brassicae) contains five, Clade 9 (incl. S. subterranea) six, Clade 15 (incl. S. nasturtii) seven and Clade 5 (incl. P. betae) eight, respectively. Clade 13 which contains W. pythii was shown to comprise eleven distinct 18S-types from five independent samples, indicating a substantial diversity of probable oomycete parasites. Clade 1 (incl. P. graminis 1) includes thirty 18S-types of which five belong to isolates.
All six unidentified environmental lineages (Clades 2, 8, 11, 12, 14, 16) were each found in one geographically/ecologically distinct sample set only: Clade 2 was found associated with plant roots of a glacial succession site, Clade 8 (monotypic) and Clade 12 were found only in South Africa fynbos soil samples, and Clade 11, Clade 14 and Clade 16 in samples from Volga flood plains.
Positive PCRs from rhizosphere samples generally could be directly sequenced without cloning, while soil samples needed to be cloned even if the diversity within the individual samples was generally low (Additional file 1: Table S4). Every soil sample contained at least two 18S-types (Additional file 1: Table S4). Nevertheless the sequenced plasmodiophorid 18S-types detected in any individual clone library were never assigned to more than four of the 16 clades defined in Figure 5.
From a German vineyard sampled repeatedly over a period of two years only five 18S-types (Figure 5, Additional file 1: Table S4) could be identified from 78 sequence reads obtained from four independent soil clone libraries. At the same time, three different species of plasmodiophorids were morphologically identified in plant samples from the same site: Sorosphaerula viticola was abundant in Vitis sp. roots, P. betae was identified from the roots of a Chenopodium sp., and Polymyxa graminis was frequently seen in the roots of Poa species. The latter was the most abundant sequence type from the soil clones and rhizosphere samples, while S. viticola and P. betae were not detected in the soil clone libraries, but were found when rhizosphere samples from primary host plants were sequenced. Therefore the 18S-types detected morphologically were a subset of those detected by sequencing: all species identified by microscopy were detected in the respective soil/rhizosphere/plant samples at least once, but not all DNA types could be assigned to actual isolates. Across all samples, plasmodiophorids were generally more easily detected in the rhizosphere of their assumed host than in neighbouring bulk soil.
Two host shifts between plant and stramenopile hosts were found in the phytomyxids (Figures 1 and 3). Prior to this study it was thought that Phagomyxida were exclusively parasites of stramenopiles and Plasmodiophorida exclusively plant parasites . However, we have shown that Woronina pythii, a parasite of the oomycete Pythium spp., branches within the Plasmodiophorida (Figure 1). Conversely, Plasmodiophora diplantherae, a parasite of an angiosperm sea grass, branches robustly amongst parasites of diatoms and brown algae. There is to our knowledge no other group of obligate biotrophic eukaryote parasites showing such pronounced cross-kingdom host shifts. In the eukaryotes cross-kingdom host shifts are known from less host-dependent opportunistic or necrotrophic parasites belonging to fungi and oomycetes [28, 30], while host shifts within other biotrophic parasites seem to occur only between relatively closely related hosts [31–34]. Biotrophic oomycetes include parasites of both nematodes (Haptoglossa spp.) and brown algae (Eurychasma spp.), but these parasites show distinct morphologies and adaptations to their hosts  suggesting strongly divergent adaptations to different host conditions. Phytomyxid parasites on the other hand have a very similar lifecycle, morphology and development inside the host (Figure 2) without any obvious specialised adaptation to a particular group of hosts. Thus, in phytomyxids, closely related lineages can be parasites of a wide range of fundamentally different hosts (autotrophic or heterotrophic) in different environments (marine, soil, freshwater).
Our morphological analyses are concordant with the cross-kingdom host jumps as revealed by molecular phylogenies. Resting spores and development of W. pythii strongly resemble that of plasmodiophorid parasites of vascular plants, especially Spongospora spp. and Sorosphaerula spp. (Figure 2a-d). The resting spores of P. diplantherae are morphologically very similar to those recently identified in Maullinia sp. from the kelp Durviellea antarctica in their ornamentation and loose aggregation (Figure 2e-g). Such loose arrangement of the resting spores is also seen in the plasmodiophorid P. brassicae, but ultrastructurally P. diplantherae is significantly different from P. brassicae (e.g. larger nucleus in P. diplantherae and differing numbers and arrangement of the synaptomal complexes) and all other plasmodiophorids . Ultrastructural studies have shown that the timing and morphology of meiosis in Woronina phythii are very similar to those in the closely related P. brassicae (Figure 1, ). Taken together, our new sequence data, phylogenetic analyses, and complementary morphological investigations question existing concepts of phytomyxid parasitism. Based on our phylogenetic analysis (Figure 1) we infer that in the marine environment parasitism of stramenopiles developed first, followed by a host switch towards angiosperm hosts. The direction of the host switch within plasmodiophorids is equivocal, because of the unresolved branching order of the deep branching taxa, including Woronina. However, we raise the hypothesis that the ancestral plasmodiophorid might have been a parasite of stramenopiles, which then bifurcated into lineages with plant hosts and those with oomycete hosts. Further sampling of the lineage diversity at the base of the phytomyxid radiation and determination of their parasitic status would enable this hypothesis to be tested further, as would characterisation of Novel Clade 9, the immediate sister group to phytomyxids , and multi-gene phylogenomic analyses of characterised members of the whole vampyrellid-phytomyxid radiation.
During the last decade it has become clear that cases of co-speciation between parasites and their hosts are rare and that the diversification of obligate parasites is more often shaped by host shifts . However, how parasites and their hosts co-evolve alongside each other remains a key issue in understanding both partners [36, 37], and can also influence policy decisions on quarantine and control measures by authorities . Phytomyxids require their primary host plant to complete their life cycle and to successfully reproduce and persist in the environment, but they can also develop infections at a low level in many alternative, non-primary host plants Figure 4, Additional file 1: Table S1; [39–41]. However, the nature of such interactions remains enigmatic. We detected phytomyxid sequences in the rhizosphere of potential alternative host plants, but could not microscopically confirm any infections of alternative hosts. The ability to infect diverse hosts without causing a distinct, visible response or by causing restricted infections over a number of generations could be a highly advantageous trait when the primary host is transient or rare. Restricted infections of phytomyxids in alternative hosts breach the first barrier to extending host range and may open the way to integrated adaptation to new hosts [42, 43]. Based on our findings, we hypothesise that the main evolutionary driver of phytomyxids is related to their ability to infect a wide range of hosts, while the main constraint is ecological. We also conclude that speciation in phytomyxids is based on host shifts rather than other co-evolutionary processes [44, 45] because of the fundamental (taxonomic) differences between the respective host species.
Phagomyxida are apparently exclusively marine and Plasmodiophorida are apparently exclusively non-marine, but we show that potentially the host range of both groups is larger than previously thought. A division in marine and freshwater/soil lineages is common in protists; very few free-living protist 18S-types are found in both marine and non-marine habitats , and in some protist groups evolutionary transitions from one environment to the other appear to be infrequent [47–49]. Ecological structuring also operates at taxonomically lower levels within plasmodiophorids: we revealed distinct, (novel) lineages and radiations from geographically and ecologically separated samples (South Africa, Volga, alpine glacial succession site, Figure 5). This supports the hypothesis that ecological barriers phylogenetically compartmentalise the diversification of phytomyxids.
Overall, we show a high 18S-type diversity of plasmodiophorids in soil ecosystems (Figure 5). Isolates from P. graminis, for example, seem to be highly variable at this level, e.g. average similarities 96.9?±?1.5% between the 30 individual 18S-types of Clade 1 (Figure 1). However, within this clade we could identify no pattern in 18S-type–host-plant species association, although the branching order is not always well resolved. It is possible that this cryptic diversity displays the natural variation within P. graminis, as result of evolutionary failure to speciate along with their hosts . But on the other hand it is possible that sequence clusters deriving from within similar environments in this clade represent P. graminis strains (or cryptic species) with different, locally advantageous host preferences, a phenomenon reported previously for P. gramins[50, 51], but otherwise are currently unknown. The nature of species complexes or cryptic species which are specialised on specific hosts/environments is increasingly recognised for fungal and oomycete plant pathogens [52, 53] and is likely to be valid for phytomyxids in general and P. graminis in particular [50, 51]. Cryptic speciation and infection of novel hosts via host shift speciation is considered to be the most frequent reason for the development of emerging diseases in fungi and oomycetes [42, 53, 54]. Host shifts between closely related parasites can also lead to more severe disease symptoms . Our findings raise the possibility that host shifts and subsequent adaptation of phytomyxids to new hosts could become more frequent or intensified in changing environments.
The overall diversity of phytomyxids in individual samples was relatively low, especially when compared to the diversity found within many groups of free-living protists, where often tens of relatively closely related but biologically distinct 18S-types can be found in a few cm3 of soil [55–58]. In soil samples of the intensively sampled vineyard (a total of 300 g soil DNA was extracted and analysed during the two year experiment) only five 18S-types belonging to two clades (clade 1, clade 13) could be identified via soil clone libraries while at the same site three species were identified microscopically and from targeted rhizosphere clones (P. graminis–clade1, P. betae–clade 5, S. viticola–clade 3). This contrasts with the detection of eleven 18S-types belonging to five species found in a smaller volume of targeted root and rhizosphere samples from the same site, highlighting that the soil clone libraries on their own underestimate the abundance and biodiversity of plasmodiophorids.
Within each of the two major phytomyxid clades there have been cross-kingdom host shifts between closely related parasite taxa: from stramenopile to angiosperm hosts in phagomyxids and between the same two host groups by plasmodiophorids. Such large-scale host shifts between closely related obligate biotrophic eukaryote parasites is to our knowledge unique to phytomyxids. The fundamental constraint in lineage radiation seems to be the environment: plasmodiophorids and phagomyxids are monophyletic groups with, respectively, non-marine and marine hosts. Environment type also seems to compartmentalise lineage radiation within plasmodiophorids at lower taxonomic levels. Our results indicate that the current research focus on a few economically significant phytomyxid lineages may lead to a limited understanding of this highly varied and largely enigmatic group of parasites of plants and stramenopiles. Some phytomyxid plant pathogens also have the ability to develop at least part of their life cycle in alternative hosts. Further work into phytomyxid diversity and their interactions with hosts in marine and non-marine (e.g. freshwater) habitats will determine (i) if high level taxonomic host shifts are more frequent, (ii) whether other high level eukaryote taxa (phyla, supergroups) are used as hosts, and (iii) whether the diversity and spatial occurrence of more closely related (i.e. around species-level) phytomyxid lineages is most strongly determined by host specificity or a wider range of environmental factors.
Samples were collected from a variety of environments (Additional file 1: Table S5), including 40 individual soil samples and 10 plant samples collected repeatedly in 2009/2010 from a German vineyard (Bernkastel, Moselle valley). The soil samples were analysed individually, whereas the plant samples (blocks of 30?×?30 cm) were pooled into plant families prior further processing to keep the sample size manageable to allow a parallel DNA-based and microscopic study. Potential host plants for phytomyxids were collected independently since 2004 to increase the taxon sampling for phylogenies. Plant samples were inspected microscopically, and if a plasmodiophorid parasite was found, DNA was extracted and sequenced (Additional file 1: Table S1).
DNA extraction and PCR
Soil samples were sieved through a 2 mm mesh size metal sieve (sterilised by flaming with Ethanol between samples) to remove stones and other large debris before DNA extraction. Root samples were picked directly from sample blocks using sterile tweezers, carefully rinsed with sterile distilled water, and used for DNA extraction. This approach of root sampling allowed tracing the selected root to the respective plant also reduced the risk of washing away fine roots which would result in a considerable loss of root biomass available for analyses.
Each soil sample was extracted in three replicates (0.5 g each). A fourth replicate of each soil sample was spiked with +/- 0.5 cm2 fungal mycelium as positive control (Roesleria subterranea strain IB 2005/0506), ensuring that negative results are not false negatives caused by PCR inhibitors. A negative control and a positive control without soil (R. subterranea strain IB 2005/0506) were included in every extraction.
Soil DNA and DNA from plant roots was extracted as described  or using the MoBio UltraClean Soil DNA isolation kit (South Africa samples). Each DNA extract was tested individually using fungal primers (V9G, LS266)  which served as yet another quality control for the DNA extraction to exclude false negative results because fungi are present in all samples, hence all samples should amplify with fungal primers.
To detect plasmodiophorids a nested PCR using two sets of 18S rDNA primers was used. First round PCR Primers were Plas1f (5′-TCAgTgAATCTgCggATggC-3′) and Plas 1r (5′-ggTgCSKCKAgRTVCAAgAggC-3′); second round primers were Plas2f (5′-TggATgTACgAGAgTACTACATgg-3′) and Plas2r (5′-CgTTgAACCTAgCATTgTAgCg-3′).
PCR-reactions (20 μL) contained final concentrations of 0.2 mM dNTPs, 0.1 μM of each primer, 1?×?PCR-reaction-buffer, 2.5 mM MgCl2, 2 μg mL-1 Bovine Serum Albumin, 0.5 U of Taq-Polymerase (Dream Taq, Fermentas) and 3 μl DNA. A touchdown PCR- protocol was used: 96°C for 4 min initial denaturation, followed by two cycles of 96°C for 25 s, 65°C for 25 s and 72°C for 1.5 min succeeded by two cycles each with an annealing temperature of 60°C and 58°C and finally 30 cycles with an annealing temperature of 54°C and completed by a final elongation step at 72°C for 10 min.
PCR products were cloned using the Fermentas insT/A cloning Kit (Thermo scientific, Germany). Positive clones were identified using the PCR primers described above. PCR products were PEG-purified (http://www.mcdb.lsa.umich.edu/labs/olsen/files/PCR.pdf) and sequenced using the specific primers (Macrogen, Korea). Sequences used in this study were deposited in Genbank accession numbers KF111152-KF111238.
Sequence manipulation, chimera check, and sequence selection
Sequences were initially blasted against GenBank (blastn) to identify phytomyxid sequences. A primary alignment of all sequences was created using MAFFT  implemented in Geneious pro (version 5.6.4; ) and this alignment was improved manually using BioEdit (version 22.214.171.124 ). To exclude chimeras Bellerophon  was used. We also compared trees computed from the alignment divided into half, and compared the respective positions of all sequences on these trees. Sequences which showed significantly different branching orders on the two trees were excluded. To group similar and/or identical sequences PhyML  and MrBayes  trees were calculated and highly similar sequences grouped and subsequently removed: sequences from individual clone libraries with a similarity of 98% or higher were treated as belonging to a single 18S-type; only one representative sequence of each 18S-type was used for subsequent analyses. The cut off-value of 98% sequence similarity was not applied to sequences from isolates or for sequences from independent clone libraries, as these sequences were considered to be ecologically and phylogenetically informative.
For phylogenetic analysis sequences were re-aligned using the l-ins-i algorithm in MAFFT and refined by eye in MacGDE. The refined alignment was analysed in RAxML BlackBox (GTR model with CAT approximation (all parameters estimated from the data); bootstrap values were mapped onto the tree with the highest likelihood value. Bayesian consensus trees were also constructed using MrBayes v 3.1.2 in parallel mode. Two separate MC3 runs with randomly generated starting trees were carried out for 6 M generations each with one cold and three heated chains. The evolutionary model included a GTR substitution matrix, a four-category autocorrelated gamma correction and the covarion model. All parameters were estimated from the data. Trees were sampled every 100 generations. 2 M generations were discarded as “burn-in” (trees sampled before the likelihood plots reached a plateau) and consensus trees constructed from the returning sample.
Plant roots were carefully washed with tap water to remove soil and organic matter attached to the roots. Samples were transferred into petri dishes with sterile, distilled water and roots that were still part of the plant of interest were selected using forceps. These roots where then either used directly for microscopy (in water) or stained with cotton blue. The roots were screened using 250 × magnification and DIC on a Nikon Optiphot light microscope, and phytomyxean structures confirmed using a 1250 × magnification.
Bait tests for the isolation of oomycete parasitic phytomyxids
A few grains of soil were added to a 9 cm petri dish filled with a mixture of sterile tap water and distilled water (1:1). Hemp seed (3 seeds per 9 cm petri dish), niger seed (Guizotia abyssinica, 3-5 seeds per 9 cm petri dish) and sesame seeds (3-5 seeds per 9 cm petri dish) were used as bait respectively. From Lohbach, Rotmoos and Bernkastel soils 20 parallel samples were tested with each bait. Samples were incubated at 10°C in the dark and checked for oomycete growth once a week for one month. The low incubation temperature slowed the growth of zygomycetes (esp. Mucor spp. and Mortierella spp.), hyphomycetes (esp. Trichoderma spp., Fusarium spp.) and bacteria which out-competed the oomycetes quickly when incubated at 20°C. Baiting cultures were first screened for oomycete growth using a reflected light microscope with 250?×?magnification. Samples which contained oomycetes were screened in an inverted microscope for the presence of cytosori. If resting spores could be identified, a microscope slide was prepared and the presence of the parasite confirmed using a Nikon Optiphot light microscope (DIC, 100?×?oil). Samples with confirmed infections were frozen for subsequent DNA extraction and analyses.
Host range meta-study
Potential primary and secondary hosts of phytomyxids were determined  supplemented by more recent studies on the host range of S. subterranea[41, 67], Polymyxa graminis[39, 50, 68], Polymyxa betae[69, 70] and Plasmodiophora brassicae. Host plants were pooled into plant families (Additional file 1: Table S2) and host–phytomyxid associations collated in an edges table (using directed, but un-weighted edges) suitable for import into Gephi 0.8.1 beta (https://gephi.org/). Nodes were defined by default from the edges table in the data laboratory module, and the initial network graph was created using the “Force Atlas” layout (settings differing from default: repulsion strength?=?10000; auto stabilisation strength?=?100; gravity?=?100). The resulting network was manually modified in Gephi until the final layout shown in Figure 4 was obtained. Fine adjustments to the layout (e.g. font colours, background colours) were made using Adobe Illustrator.
Availability of supporting data
Sequences used in this study were deposited in Genbank accession numbers KF111152-KF111238. Any associated data are summarised in Additional file 1.
Braselton JP: Plasmodiophoromycota. The Mycota VII Part A Systematics and Evolution. Edited by: McLaughlin DJ, McLaughlin EG, Lemke PA. 2001, Berlin: Springer-Verlag, 81-91.
Neuhauser S, Kirchmair M, Gleason FH: Ecological roles of the parasitic phytomyxids (plasmodiophorids) in marine ecosystems - a review. Mar Freshw Res. 2011, 62: 365-371. 10.1071/MF10282.
Bass D, Chao EEY, Nikolaev S, Yabuki A, Ishida K-i, Berney C, Pakzad U, Wylezich C, Cavalier-Smith T: Phylogeny of novel naked filose and reticulose cercozoa: granofilosea cl. n. and proteomyxidea revised. Protist. 2009, 160: 75-109. 10.1016/j.protis.2008.07.002.
Burki F, Kudryavtsev A, Matz MV, Aglyamova GV, Bulman S, Fiers M, Keeling PJ, Pawlowski J: Evolution of Rhizaria: new insights from phylogenomic analysis of uncultivated protists. BMC Evol Biol. 2010, 10: 377-10.1186/1471-2148-10-377. doi:10.1186/1471-2148-10-377
Bulman SR, Kuhn SF, Marshall JW, Schnepf E: A phylogenetic analysis of the SSU rRNA from members of the plasmodiophorida and phagomyxida. Protist. 2001, 152: 43-51. 10.1078/1434-4610-00042.
Goecke F, Wiese J, Nunez A, Labes A, Imhoff JF, Neuhauser S: A Novel phytomyxean parasite associated with galls on the bull-kelp Durvillaea antarctica (Chamisso) Hariot. PLoS One. 2012, 7 (9): e45358-10.1371/journal.pone.0045358. doi:10.1371/journal.pone.0045358
Schnepf E, Kühn SF, Bulman S: Phagomyxa bellerocheae sp. nov. and Phagomyxa odontellae sp. nov., Plasmodiophoromycetes feeding on marine diatoms. Helgol Mar Res. 2000, 54: 237-241. 10.1007/s101520000056.
Maier I, Parodi E, Westermeier R, Muller DG: Maullinia ectocarpii gen. et sp nov (Plasmodiophorea), an intracellular parasite in Ectocarpus siliculosus (Ectocarpales, Phaeophyceae) and other filamentous brown algae. Protist. 2000, 151: 225-238. 10.1078/1434-4610-00021.
Takishita K, Miyake H, Kawato M, Maruyama T: Genetic diversity of microbial eukaryotes in anoxic sediment around fumaroles on a submarine caldera floor based on the small-subunit rDNA phylogeny. Extremophiles. 2005, 9: 185-196. 10.1007/s00792-005-0432-9.
Takishita K, Tsuchiya M, Kawato M, Oguri K, Kitazato H, Maruyama T: Genetic diversity of microbial eukaryotes in anoxic sediment of the saline meromictic lake Namako-ike (Japan): on the detection of anaerobic or anoxic-tolerant lineages of eukaryotes. Protist. 2007, 158: 51-64. 10.1016/j.protis.2006.07.003.
Dixon GR: The occurrence and economic impact of Plasmodiophora brassicae and clubroot disease. J Plant Growth Regul. 2009, 28: 194-202. 10.1007/s00344-009-9090-y.
Merz U, Falloon RE: Review: powdery scab of potato-increased knowledge of pathogen biology and disease epidemiology for effective disease management. Potato Res. 2009, 52: 17-37. 10.1007/s11540-008-9105-2.
Kanyuka K, Ward E, Adams MJ: Polymyxa graminis and the cereal viruses it transmits: a research challenge. Mol Plant Pathol. 2003, 4: 393-406. 10.1046/j.1364-3703.2003.00177.x.
Varrelmann M: Occurrence, spread and pathogenicity of different forms of the Rhizomania virus (Beet necrotic yellow vein virus, BNYVV) - review on biology and variability of Rhizomania and on detection of isolates possibly overcoming resistance. Zuckerindustrie. 2007, 132: 113-120.
Neuhauser S, Bulman S, Kirchmair M: Plasmodiophorids: the challenge to understand soil-borne, obligate biotrophs with a multiphasic life cycle. Molecular Identification of Fungi. Edited by: Gherbawy Y, Voigt K. 2010, Berlin Heidelberg: Springer, 51-78.
Ludwig-Mueller J, Prinsen E, Rolfe SA, Scholes JD: Metabolism and plant hormone action during clubroot disease. J Plant Growth Regul. 2009, 28: 229-244. 10.1007/s00344-009-9089-4.
Siemens J, Graf H, Bulman S, In O, Ludwig-Mueller J: Monitoring expression of selected Plasmodiophora brassicae genes during clubroot development in Arabidopsis thaliana. Plant Pathol. 2009, 58: 130-136. 10.1111/j.1365-3059.2008.01943.x.
Kageyama K, Asano T: Life cycle of Plasmodiophora brassicae. J Plant Growth Regul. 2009, 28: 203-211. 10.1007/s00344-009-9101-z.
Neuhauser S, Kirchmair M, Gleason FH: The ecological potentials of Phytomyxea (“plasmodiophorids”) in aquatic food webs. Hydrobiologia. 2011, 659: 23-35. 10.1007/s10750-010-0508-0.
Karling JS: The Plasmodiophorales. 1968, New York: Hafner Publishing Company
Cook W: A monograph of the Plasmodiophorales. Arch Protistenkunde. 1933, 80: 179-254.
Dylewski DP, Miller CE: The ultrastructure of mitosis during sporangiogenesis in Woronina pythii (Plasmodiophoromycetes). Protoplasma. 1984, 121: 42-53. 10.1007/BF01279751.
Lafferty KD, Porter JW, Ford SE: Are diseases increasing in the ocean?. Annu Rev Ecol Evol Syst. 2004, 31-54.
Poulin R, Paterson R, Townsend C, Tompkins D, Kelly D: Biological invasions and the dynamics of endemic diseases in freshwater ecosystems. Freshwat Biol. 2010, 56: 676-688.
Gachon CMM, Sime-Ngando T, Strittmatter M, Chambouvet A, Kim GH: Algal diseases: spotlight on a black box. Trends Plant Sci. 2010, 15: 633-640. 10.1016/j.tplants.2010.08.005.
Karling JS: A synopsis of Rozella and Rozellopsis. Mycologia. 1942, 34: 193-208. 10.2307/3754811.
Lesaulnier C, Papamichail D, McCorkle S, Ollivier B, Skiena S, Taghavi S, Zak D, Van Der Lelie D: Elevated atmospheric CO2 affects soil microbial diversity associated with trembling aspen. Environ Microbiol. 2008, 10: 926-941. 10.1111/j.1462-2920.2007.01512.x.
Beakes GW, Glockling SL, Sekimoto S: The evolutionary phylogeny of the oomycete “fungi”. Protoplasma. 2012, 249: 3-19.
Cronk QCB: Plant evolution and development in a post-genomic context. Nat Rev Genet. 2001, 2: 607-619. 10.1038/35084556.
Van Baarlen P, Van Belkum A, Summerbell RC, Crous PW, Thomma BP: Molecular mechanisms of pathogenicity: how do pathogenic microorganisms develop cross‐kingdom host jumps?. FEMS Microbiol Rev. 2007, 31: 239-277. 10.1111/j.1574-6976.2007.00065.x.
Goeker M, Voglmayr H, Riethmueller A, Oberwinkler F: How do obligate parasites evolve? A multi-gene phylogenetic analysis of downy mildews. Fungal Genet Biol. 2007, 44: 105-122. 10.1016/j.fgb.2006.07.005.
Lara E, Belbahri L: SSU rRNA reveals major trends in oomycete evolution. Fungal Divers. 2011, 49: 93-100. 10.1007/s13225-011-0098-9.
Takamatsu S: Phylogeny and evolution of the powdery mildew fungi (Erysiphales, Ascomycota) inferred from nuclear ribosomal DNA sequences. Mycoscience. 2004, 45: 147-157. 10.1007/S10267-003-0159-3.
Van der Merwe M, Ericson L, Walker J, Thrall PH, Burdon JJ: Evolutionary relationships among species of Puccinia and Uromyces (Pucciniaceae, Uredinales) inferred from partial protein coding gene phylogenies. Mycol Res. 2007, 111: 163-175. 10.1016/j.mycres.2006.09.015.
Braselton JP, Short FT: Karyotypic analysis of Plasmodiophora diplantherae. Mycologia. 1985, 77: 940-945. 10.2307/3793306.
de Vienne DM, Refrégier G, López-Villavicencio M, Tellier A, Hood ME, Giraud T: Cospeciation vs host-shift speciation: methods for testing, evidence from natural associations and relation to coevolution. New Phytol. 2013, 198: 347-385. 10.1111/nph.12150.
Yoder JB, Nuismer SL: When does coevolution promote diversification?. Am Nat. 2010, 176: 802-817. 10.1086/657048.
Gilbert GS, Webb CO: Phylogenetic signal in plant pathogen–host range. Proc Natl Acad Sci U S A. 2007, 104: 4979-4983. 10.1073/pnas.0607968104.
Desoignies N, Stocco C, Bragard C, Legreve A: A new phenotype of Polymyxa betae in Arabidopsis thaliana. Eur J Plant Pathol. 2011, 131: 27-38. 10.1007/s10658-011-9783-5.
Ludwig-Muller J, Bennett RN, Kiddle G, Ihmig S, Ruppel M, Hilgenberg W: The host range of Plasmodiophora brassicae and its relationship to endogenous glucosinolate content. New Phytol. 1999, 141: 443-458. 10.1046/j.1469-8137.1999.00368.x.
Qu X, Christ BJ: The host range of Spongospora subterranea f. sp subterranea in the United States. Am J Potato Res. 2006, 83: 343-347. 10.1007/BF02871595.
Giraud T, Gladieux P, Gavrilets S: Linking the emergence of fungal plant diseases with ecological speciation. Trends Ecol Evol. 2010, 25: 387-395. 10.1016/j.tree.2010.03.006.
Antonovics J, Boots M, Ebert D, Koskella B, Poss M, Sadd BM: The origin of specificity by means of natural selection: evolved and nonhost resistance in host–pathogen interactions. Evolution. 2013, 67: 1-9. 10.1111/j.1558-5646.2012.01793.x.
Johnson KP, Adams RJ, Page RD, Clayton DH: When do parasites fail to speciate in response to host speciation?. Syst Biol. 2003, 52: 37-47. 10.1080/10635150390132704.
Libeskind-Hadas R, Charleston MA: On the computational complexity of the reticulate cophylogeny reconstruction problem. J Comput Biol. 2009, 16: 105-117. 10.1089/cmb.2008.0084.
Richards TA, Bass D: Molecular screening of free-living microbial eukaryotes: diversity and distribution using a meta-analysis. Curr Opin Microbiol. 2005, 8: 240-252. 10.1016/j.mib.2005.04.010.
Logares R, Brate J, Bertilsson S, Clasen JL, Shalchian-Tabrizi K, Rengefors K: Infrequent marine-freshwater transitions in the microbial world. Trends Microbiol. 2009, 17: 414-422. 10.1016/j.tim.2009.05.010.
Bråte J, Logares R, Berney C, Ree DK, Klaveness D, Jakobsen KS, Shalchian-Tabrizi K: Freshwater Perkinsea and marine-freshwater colonizations revealed by pyrosequencing and phylogeny of environmental rDNA. ISME J. 2010, 4: 1144-1153. 10.1038/ismej.2010.39.
Shalchian‐Tabrizi K, Bråte J, Logares R, Klaveness D, Berney C, Jakobsen KS: Diversification of unicellular eukaryotes: cryptomonad colonizations of marine and fresh waters inferred from revised 18S rRNA phylogeny. Environ Microbiol. 2008, 10: 2635-2644. 10.1111/j.1462-2920.2008.01685.x.
Legreve A, Vanpee B, Delfosse P, Maraite H: Host range of tropical and sub-tropical isolates of Polymyxa graminis. Eur J Plant Pathol. 2000, 106: 379-389. 10.1023/A:1008784823899.
Smith MJ, Adams MJ, Ward E: Ribosomal DNA analyses reveal greater sequence variation in Polymyxa species than previously thought and indicate the possibility of new ribotype-host-virus associations. Environ Microbiol Rep. 2013, 5: 143-150. 10.1111/1758-2229.12026.
Crous P, Groenewald J: Hosts, species and genotypes: opinions versus data. Australas Plant Pathol. 2005, 34: 463-470. 10.1071/AP05082.
Rouxel M, Mestre P, Comont G, Lehman BL, Schilder A, Delmotte F: Phylogenetic and experimental evidence for host-specialized cryptic species in a biotrophic oomycete. New Phytol. 2013, 197: 251-263. 10.1111/nph.12016.
Jiang RH, Tyler BM: Mechanisms and evolution of virulence in oomycetes. Annu Rev Phytopathol. 2012, 50: 295-318. 10.1146/annurev-phyto-081211-172912.
Bass D, Richards TA, Matthai L, Marsh V, Cavalier-Smith T: DNA evidence for global dispersal and probable endemicity of protozoa. BMC Evol Biol. 2007, 7: 162-10.1186/1471-2148-7-162.
Howe AT, Bass D, Vickerman K, Chao EE, Cavalier-Smith T: Phylogeny, taxonomy, and astounding genetic diversity of glissomonadida ord. nov., the dominant gliding zooflagellates in soil (Protozoa: Cercozoa). Protist. 2009, 160: 159-189. 10.1016/j.protis.2008.11.007.
Lara E, Mitchell EA, Moreira D, López García P: Highly diverse and seasonally dynamic protist community in a pristine peat bog. Protist. 2011, 162: 14-32. 10.1016/j.protis.2010.05.003.
von der Heyden S, Cavalier-Smith T: Culturing and environmental DNA sequencing uncover hidden kinetoplastid biodiversity and a major marine clade within ancestrally freshwater Neobodo designis. Int J Syst Evol Microbiol. 2005, 55: 2605-2621. 10.1099/ijs.0.63606-0.
Neuhauser S, Huber L, Kirchmair M: A DNA based method to detect the grapevine root-rotting fungus Roesleria subterranea in soil and root samples. Phytopathol Mediterr. 2009, 48: 59-72.
Gerrits van Den Ende AHGG, De Hoog GS: Variability and molecular diagnostics of the neurotropic species Cladophialophora bantiana. Stud Mycol. 1999, 43: 151-162.
Katoh K, Kuma K-I, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005, 33: 511-518. 10.1093/nar/gki198.
Drummond A, Ashton B, Buxton S, Cheung M, Cooper A, Duran C, Heled J, Kearse M, Markowitz S, Moir R, et al: Geneious, version 5.6.4. 2012
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999, 41: 95-98.
Huber T, Faulkner G, Hugenholtz P: Bellerophon: a program to detect chimeric sequences in multiple sequence alignments. Bioinformatics. 2004, 20: 2317-2319. 10.1093/bioinformatics/bth226.
Guindon S, Delsuc F, Dufayard J-F, Gascuel O: Estimating maximum likelihood phylogenies with PhyML. Methods Mol Biol (Clifton, NJ). 2009, 537: 113-137. 10.1007/978-1-59745-251-9_6.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Andersen BAB, Nicolaisen M, Nielsen SL: Alternative hosts for potato mop-top virus, genus Pomovirus and its vector Spongospora subterranea f.sp subterranea. Potato Res. 2002, 45: 37-43. 10.1007/BF02732217.
Legreve A, Schmit J, Bragard C, Maraite H: The role of climate and alternative hosts in the epidemiology of rhizomanina. Symposium of the International Working Group on Plant Viruses with Fungal Vectors. Edited by: Rush C. 2005, Bologna, Italy: , 129-132.
Barr KJ, Asher MJC: The host range of Polymyxa betae in Britain. Plant Pathol. 1992, 41: 64-68. 10.1111/j.1365-3059.1992.tb02317.x.
Hugo SA, Henry CM, Harju V: The role of alternative hosts of Polymyxa betae in transmission of beet necrotic yellow vein virus (BNYVV) in England. Plant Pathol. 1996, 45: 662-666. 10.1046/j.1365-3059.1996.d01-182.x.
The authors acknowledge Lars Huber, Mark Hoffmann, Fred Unmann and Martin Schletterer for sample collection; Allison Walker for providing P. diplantherae material, and F. Goecke for providing material of Maullinia sp. Sigrid Neuhauser gratefully acknowledges funding by the Austrian Science Fund (FWF) through a Hertha-Firnberg research grant (T379-B16) and an Erwin Schrödinger grant (J3175-B20), the University Innsbruck (718001 Nachwuchsförderung) and the Forschungszentrum Berglandwirtschaft (University Innsbruck). DB was supported by NERC New Investigator (NE/H000887/1) and Standard Research (NE/H009426/1) grants. DB and SN were also supported by the SynTax project ‚Phytomyxids–a hidden world of parasites and symbionts’.
The authors declare that they have no competing interests.
SN and MK designed research, SN performed research, SN, MK, SB and DB analyzed data and wrote the paper. All authors read and approved the final manuscript.
Sigrid Neuhauser, Martin Kirchmair contributed equally to this work.
Electronic supplementary material
Additional file 1: Table S1: Microscopically observed isolates. Table S2. Data and references used for meta-analysis of phytomyxid host range, including data used to compile Figure 4. Table S3. Dataset of environmental sequences shown in Figure 5, including information on samples and abundances of the individual 18-S types in the clone libraries. Table S4. Sample numbers, positive detections using the plasmodiophorid specific primer strategy, and 18S rDNA types within each sample pool. Table S5. Samples, sampling dates and sampling sites. (PDF 186 KB)
About this article
Cite this article
Neuhauser, S., Kirchmair, M., Bulman, S. et al. Cross-kingdom host shifts of phytomyxid parasites. BMC Evol Biol 14, 33 (2014). https://doi.org/10.1186/1471-2148-14-33
- Plant pathogen
- Host range