Pyrosequencing as a tool for the detection ofPhytophthoraspecies: error rate and risk of false Molecular Operational Taxonomic Units

The Never-Ending Presence of Phytophthora Species in Italian Nurseries

Plant trade coupled with climate change has led to the increased spread of well-known and new Phytophthora species, a group of fungus-like organisms placed in the Kingdom Chromista. Their presence in plant nurseries is of particular concern because they are responsible for many plant diseases, with high environmental, economic and social impacts. This paper offers a brief overview of the current status of Phytophthora species in European plant nurseries. Focus was placed on Italian sites. Despite the increasing awareness of the risk of Phytophthora spread and the management strategies applied for controlling it, the complexity of the Phytophthora community in the horticulture industry is increasing over time. Since the survey carried out by Jung et al. (2016), new Phytophthora taxa and Phytophthora-host associations were identified. Phytophthorahydropathica, P. crassamura, P. pseudocryptogea and P. meadii were reported for the first time in European plant nurseries, while P. pistaciae, P. mediterranea and P. heterospora were isolated from Italian ornamental nurseries. Knowledge of Phytophthora diversity in plant nurseries and the potential damage caused by them will help to contribute to the development of early detection methods and sustainable management strategies to control Phytophthora spread in the future.

Diversity of Phytophthora Species Detected in Disturbed and Undisturbed British Soils Using High-Throughput Sequencing Targeting ITS rRNA and COI mtDNA Regions

Disease outbreaks caused by introduced Phytophthora species have been increasing in British forests and woodlands in recent years. A better knowledge of the Phytophthora communities already present in the UK is of great importance when developing management and mitigation strategies for these diseases. To do this, soils were sampled in “disturbed” sites, meaning sites frequently visited by the public, with recent and new plantings or soil disturbances versus more “natural” forest and woodland sites with little disturbance or management. Phytophthora diversity was assessed using high-throughput Illumina sequencing targeting the widely accepted barcoding Internal Transcribed Spacer 1 (ITS1) region of rRNA and comparing it with the mitochondrial cytochrome c oxidase I (COI) gene. Isolation of Phytophthora was run in parallel. Nothophytophthora spp. and Phytophthora spp. were detected in 79 and 41 of the 132 locations of the 14 studied sites when using ITS or COI, respectively. A total of 20 Phytophthora amplicon sequence variants (ASVs) were assigned to known Phytophthora species from eight clades (1a, 2, 2b, 3a, 5, 6b, 7a, 8b, 8c, 8d, 10a, and 10b) and 12 ASVs from six clades (1a, 2c, 3a, 3b, 6b, 7a, 8b, 8c, and 8d) when using ITS or COI, respectively. Only at two locations were the results in agreement for ITS, COI, and isolation. Additionally, 21 and 17 unknown Phytophthora phylotypes were detected using the ITS and COI, respectively. Several Phytophthora spp. within clades 7 and 8, including very important forest pathogens such as P. austrocedri and P. ramorum, were identified and found more frequently at “disturbed” sites. Additionally, eight ASVs identified as Nothophytophthora spp. were detected representing the first report of species within this new genus in Britain. Only three species not known to be present in Britain (P. castaneae, P. capsici, and P. fallax) were detected with the ITS primers and not with COI. To confirm the presence of these or any potential new Phytophthora species, sites should be re-sampled for confirmation. Additionally, there is a need to confirm if these species are a threat to British trees and try to establish any eradication measures required to mitigate Phytophthora spread in Britain.

Surprising low diversity of the plant pathogen Phytophthora in Amazonian forests

The genus Phytophthora represents a group of plant pathogens with broad global distribution. The majority of them cause the collar and root-rot of diverse plant species. Little is known about Phytophthora communities in forest ecosystems, especially in the Neotropical forests where natural enemies could maintain the huge plant diversity via negative density dependence. We characterized the diversity of soil-borne Phytophthora communities in the North French Guiana rainforest and investigated how they are structured by host identity and environmental factors. In this little-explored habitat, 250 soil cores were sampled from 10 plots hosting 10 different plant families across three forest environments (Terra Firme, Seasonally Flooded and White Sand). Phytophthora diversity was studied using a baiting approach and metabarcoding (High-Throughput Sequencing) on environmental DNA extracted from both soil samples and baiting-leaves. These three approaches revealed very similar communities, characterized by an unexpected low diversity of Phytophthora species, with the dominance of two cryptic species close to Phytophthora heveae. As expected, the Phytophthora community composition of the French Guiana rainforest was significantly impacted by the host plant family and environment. However, these plant pathogen communities are very small and are dominated by generalist species, questioning their potential roles as drivers of plant diversity in these Amazonian forests.

Diversity of Phytophthora Species Associated with Quercus ilex L. in Three Spanish Regions Evaluated by NGS

The diversity of Phytophthora species in declining Fagaceae forests in Europe is increasing in the last years. The genus Quercus is one of the most extended Fagaceae genera in Europe, and Q. ilex is the dominant tree in Spain. The introduction of soil-borne pathogens, such as Phytophthora in Fagaceae forests modifies the microbial community present in the rhizosphere, and has relevant environmental and economic consequences. A better understanding of the diversity of Phytophthora spp. associated with Q. ilex is proposed in this study by using Next Generation Sequencing (NGS) in six Q. ilex stands located in three regions in Spain. Thirty-seven Phytophthora phylotypes belonging to clades 1 to 12, except for clades 4, 5 and 11, are detected in this study, which represents a high diversity of Phytophthora species in holm oak Spanish forests. Phytophthora chlamydospora, P. citrophthora, P. gonapodyides, P. lacustris, P. meadii, P. plurivora, P. pseudocryptogea, P. psychrophila and P. quercina were present in the three regions. Seven phylotypes could not be associated with known Phytophthora species, so they were putatively named as Phytophthora sp. Most of the detected phylotypes corresponded to terrestrial Phytophthora species but aquatic species from clades 6 and 9 were also present in all regions.

Site-specific distribution of oak rhizosphere-associated oomycetes revealed by cytochrome c oxidase subunit II metabarcoding

The phylum Oomycota comprises important tree pathogens like Phytophthora quercina, involved in central European oak decline, and Phytophthora cinnamomi shown to affect holm oaks among many other hosts. Despite the importance to study the distribution, dispersal and niche partitioning of this phylum, metabarcoding surveys, and studies considering environmental factors that could explain oomycete community patterns are still rare. We investigated oomycetes in the rhizosphere of evergreen oaks in a Spanish oak woodland using metabarcoding based on Illumina sequencing of the taxonomic marker cytochrome c oxidase subunit II (cox2). We developed an approach amplifying a 333 bp long fragment using the forward primer Hud-F (Mycologia, 2000) and a reverse primer found using DegePrime (Applied and Environmental Microbiology, 2014). Factors reflecting topo-edaphic conditions and tree health were linked to oomycete community patterns. The majority of detected OTUs belonged to the Peronosporales. Most taxa were relatives of the Pythiaceae, but relatives of the Peronosporaceae and members of the Saprolegniales were also found. The most abundant OTUs were related to Globisporangium irregulare and P. cinnamomi, both displaying strong site-specific patterns. Oomycete communities were strongly correlated with the environmental factors: altitude, crown foliation, slope and soil skeleton and soil nitrogen. Our findings illustrate the significance of small scale variation in habitat conditions for the distribution of oomycetes and highlight the importance to study oomycete communities in relation to such ecological patterns.

High-throughput identification and diagnostics of pathogens and pests: Overview and practical recommendations

High-throughput identification technologies provide efficient tools for understanding the ecology and functioning of microorganisms. Yet, these methods have been only rarely used for monitoring and testing ecological hypotheses in plant pathogens and pests in spite of their immense importance in agriculture, forestry and plant community dynamics. The main objectives of this manuscript are the following: (a) to provide a comprehensive overview about the state-of-the-art high-throughput quantification and molecular identification methods used to address population dynamics, community ecology and host associations of microorganisms, with a specific focus on antagonists such as pathogens, viruses and pests; (b) to compile available information and provide recommendations about specific protocols and workable primers for bacteria, fungi, oomycetes and insect pests; and (c) to provide examples of novel methods used in other microbiological disciplines that are of great potential use for testing specific biological hypotheses related to pathology. Finally, we evaluate the overall perspectives of the state-of-the-art and still evolving methods for diagnostics and population- and community-level ecological research of pathogens and pests.

Canker and decline diseases caused by soil- and airborne Phytophthora species in forests and woodlands

Most members of the oomycete genus Phytophthora are primary plant pathogens. Both soil- and airborne Phytophthora species are able to survive adverse environmental conditions with enduring resting structures, mainly sexual oospores, vegetative chlamydospores and hyphal aggregations. Soilborne Phytophthora species infect fine roots and the bark of suberized roots and the collar region with motile biflagellate zoospores released from sporangia during wet soil conditions. Airborne Phytophthora species infect leaves, shoots, fruits and bark of branches and stems with caducous sporangia produced during humid conditions on infected plant tissues and dispersed by rain and wind splash. During the past six decades, the number of previously unknown Phytophthora declines and diebacks of natural and semi-natural forests and woodlands has increased exponentially, and the vast majority of them are driven by introduced invasive Phytophthora species. Nurseries in Europe, North America and Australia show high infestation rates with a wide range of mostly exotic Phytophthora species. Planting of infested nursery stock has proven to be the main pathway of Phytophthora species between and within continents. This review provides insights into the history, distribution, aetiology, symptomatology, dynamics and impact of the most important canker, decline and dieback diseases caused by soil- and airborne Phytophthora species in forests and natural ecosystems of Europe, Australia and the Americas.

Oomycete Communities Associated with Reed Die-Back Syndrome

Phragmites australis (Cav.) Trin. ex Steud. die-back is a widely-studied phenomenon that was first discovered in northern Europe and that, until recently, was almost unknown in the Mediterranean basin. It has been described as a complex syndrome affecting reed populations leading to their retreat and decline, with significant impacts on valuable ecosystem services. Among the factors that cause the decline, soil-living microorganisms can be crucial. The aims of this study were to analyze the diversity of oomycetes communities associated with reed stands, and to understand whether they could play a key role in the decline. Variations in the structure of oomycetes communities were studied by metabarcoding of the internal transcribed spacer (ITS) 1 region of ribosomal DNA, from the sediments of five Italian freshwater ecosystems. They were chosen to cover a large variability in terms of surface area, water depth, microclimate, and presence of documented reed retreat. From 96 samples collected from reed roots, rhizosphere, and bulk soil, we assembled 207661 ITS1 reads into 523 OTUs. We demonstrated that oomycete communities were structured by several factors, among which the most important was die-back occurrence. Our study also indicates that Pythiogeton spp. could be potentially involved in the development of die-back. The role of heavy metals in the soil was also explored, and cadmium concentration was shown to affect oomycetes distribution. This study represents a significant step forward for the characterization of microbial communities associated with reed die-back syndrome and helps to gain knowledge of the complexity of these important wet ecosystems.

Two previously unknown Phytophthora species associated with brown rot of Pomelo (Citrus grandis) fruits in Vietnam

Two distinct Phytophthora taxa were found to be associated with brown rot of pomelo (Citrus grandis), a new disease of this ancestral Citrus species, in the Vinh Long province, Mekong River Delta area, southern Vietnam. On the basis of morphological characters and using the ITS1-5.8S-ITS2 region of the rDNA and the cytochrome oxidase subunit 1 (COI) as barcode genes, one of the two taxa was provisionally named as Phytophthora sp. prodigiosa, being closely related to but distinct from P. insolita, a species in Phytophthora Clade 9, while the other one, was closely related to but distinct from the Clade 2 species P. meadii and was informally designated as Phytophthora sp. mekongensis. Isolates of P. sp. prodigiosa and P. sp. mekongensis were also obtained from necrotic fibrous roots of Volkamer lemon (C. volkameriana) rootstocks grafted with 'King' mandarin (Citrus nobilis) and from trees of pomelo, respectively, in other provinces of the Mekong River Delta, indicating a widespread occurrence of both Phytophthora species in this citrus-growing area. Koch's postulates were fulfilled via pathogenicity tests on fruits of various Citrus species, including pomelo, grapefruit (Citrus x paradisi), sweet orange (Citrus x sinensis) and bergamot (Citrus x bergamia) as well as on the rootstock of 2-year-old trees of pomelo and sweet orange on 'Carrizo' citrange (C. sinensis 'Washington Navel' x Poncirus trifoliata). This is the first report of a Phytophthora species from Clade 2 other than P. citricola and P. citrophthora as causal agent of fruit brown rot of Citrus worldwide and the first report of P. insolita complex in Vietnam. Results indicate that likely Vietnam is still an unexplored reservoir of Phytophthora diversity.

Characterization of Pythium spp. Associated with Asymptomatic Soybean in Southeastern Pennsylvania

Soybean production in Pennsylvania has increased substantially over the past 20 years and is a highly valued field crop, together with corn. Soilborne pathogens such as Pythium spp. can contribute to soybean stand establishment issues, particularly under the conservation tillage practices that are common in the state. In this study, we collected soil samples from eight asymptomatic soybean-corn rotation fields across six counties in southeastern Pennsylvania between May and June 2012. Pythium spp. were isolated via baiting, and tested for aggressiveness on both soybean and corn using laboratory assays. In addition to our culture-based survey, we also assessed the use of cytochrome oxidase subunit 1 pyrosequencing as a culture-independent method for measuring Pythium spp. diversity from environmental samples. Diversity estimates were consistent between the culture-based and pyrosequencing datasets; however, important methodological biases inherent to culture-independent methods may have led to some differences. Our results show that several Pythium spp. previously characterized as soybean or corn pathogens are present in southeastern Pennsylvania, including Pythium irregulare, P. sylvaticum, and P. ultimum var. sporangiiferum, with isolates showing aggressive phenotypes in lab assays.

Metabarcoding Analysis of Phytophthora Diversity Using Genus-Specific Primers and 454 Pyrosequencing

A metabarcoding method based on genus-specific primers and 454 pyrosequencing was utilized to investigate the genetic diversity of Phytophthora spp. in soil and root samples of potted plants, from eight nurseries. Pyrosequencing enabled the detection of 25 Phytophthora phylotypes distributed in seven different clades and provided a much higher resolution than a corresponding cloning/Sanger sequencing approach. Eleven of these phylotypes, including P. cactorum, P. citricola s.str., P. palmivora, P. palmivora-like, P. megasperma or P. gonapodyides, P. ramorum, and five putative new Phytophthora species phylogenetically related to clades 1, 2, 4, 6, and 7 were detected only with the 454 pyrosequencing approach. We also found an additional 18 novel records of a phylotype in a particular nursery that were not detected with cloning/Sanger sequencing. Several aspects confirmed the reliability of the method: (i) many identical sequence types were identified independently in different nurseries, (ii) most sequence types identified with 454 pyrosequencing were identical to those from the cloning/Sanger sequencing approach and/or perfectly matched GenBank deposited sequences, and (iii) the divergence noted between sequence types of putative new Phytophthora species and all other detected sequences was sufficient to rule out sequencing errors. The proposed method represents a powerful tool to study Phytophthora diversity providing that particular attention is paid to the analysis of 454 pyrosequencing raw read sequences and to the identification of sequence types.

The use of genus-specific amplicon pyrosequencing to assess phytophthora species diversity using eDNA from soil and water in Northern Spain

Phytophthora is one of the most important and aggressive plant pathogenic genera in agriculture and forestry. Early detection and identification of its pathways of infection and spread are of high importance to minimize the threat they pose to natural ecosystems. eDNA was extracted from soil and water from forests and plantations in the north of Spain. Phytophthora-specific primers were adapted for use in high-throughput Sequencing (HTS). Primers were tested in a control reaction containing eight Phytophthora species and applied to water and soil eDNA samples from northern Spain. Different score coverage threshold values were tested for optimal Phytophthora species separation in a custom-curated database and in the control reaction. Clustering at 99% was the optimal criteria to separate most of the Phytophthora species. Multiple Molecular Operational Taxonomic Units (MOTUs) corresponding to 36 distinct Phytophthora species were amplified in the environmental samples. Pyrosequencing of amplicons from soil samples revealed low Phytophthora diversity (13 species) in comparison with the 35 species detected in water samples. Thirteen of the MOTUs detected in rivers and streams showed no close match to sequences in international sequence databases, revealing that eDNA pyrosequencing is a useful strategy to assess Phytophthora species diversity in natural ecosystems.

An improved high throughput sequencing method for studying oomycete communities

Culture-independent studies using next generation sequencing have revolutionized microbial ecology, however, oomycete ecology in soils is severely lagging behind. The aim of this study was to improve and validate standard techniques for using high throughput sequencing as a tool for studying oomycete communities. The well-known primer sets ITS4, ITS6 and ITS7 were used in the study in a semi-nested PCR approach to target the internal transcribed spacer (ITS) 1 of ribosomal DNA in a next generation sequencing protocol. These primers have been used in similar studies before, but with limited success. We were able to increase the proportion of retrieved oomycete sequences dramatically mainly by increasing the annealing temperature during PCR. The optimized protocol was validated using three mock communities and the method was further evaluated using total DNA from 26 soil samples collected from different agricultural fields in Denmark, and 11 samples from carrot tissue with symptoms of Pythium infection. Sequence data from the Pythium and Phytophthora mock communities showed that our strategy successfully detected all included species. Taxonomic assignments of OTUs from 26 soil sample showed that 95% of the sequences could be assigned to oomycetes including Pythium, Aphanomyces, Peronospora, Saprolegnia and Phytophthora. A high proportion of oomycete reads was consistently present in all 26 soil samples showing the versatility of the strategy. A large diversity of Pythium species including pathogenic and saprophytic species were dominating in cultivated soil. Finally, we analyzed amplicons from carrots with symptoms of cavity spot. This resulted in 94% of the reads belonging to oomycetes with a dominance of species of Pythium that are known to be involved in causing cavity spot, thus demonstrating the usefulness of the method not only in soil DNA but also in a plant DNA background. In conclusion, we demonstrate a successful approach for pyrosequencing of oomycete communities using ITS1 as the barcode sequence with well-known primers for oomycete DNA amplification.

Optimization of Headspace Solid-Phase Microextraction Conditions for the Identification of Phytophthora cinnamomi Rands

A robust technique was developed to identify Phytophthora cinnamomi using headspace solid-phase microextraction (HS-SPME) combined with gas chromatography (GC) coupled to a flame ionization detector (FID) for analyzing volatile organic compounds (VOCs). Six fiber types were evaluated and results indicated that the three-phase fiber 50/30 μm divinylbenzene/carboxen/polydimethylsiloxane (DVB/CAR/PDMS) had the highest extraction efficiency for both polar and nonpolar GC columns. The maximum extraction efficiency (equilibrium absorption) was achieved 16 h after fiber exposure in the HS. Absorbed compounds on the fiber were completely desorbed in the GC injector after 5 min at 250°C. Compared with the nonpolar column, the polar column showed optimum separation of VOCs released from P. cinnamomi. Under the optimized HS-SPME and GC/FID conditions, lower detection limits for the four external standards was found to be between 1.57 to 27.36 ng/liter. Relative standard deviations <9.010% showed that the method is precise and reliable. The method also showed good linearity for the concentration range that was analyzed using four standards, with regression coefficients between 0.989 and 0.995, and the sensitivity of the method was 10⁴ times greater than that of the conventional HS method. In this study, the VOC profiles of six Phytophthora spp. and one Pythium sp. were characterized by the optimized HS-SPME-GC method. The combination of the VOCs creates a unique pattern for each pathogen; the chromatograms of different isolates of P. cinnamomi were the same and the specific VOC pattern of P. cinnamomi remained consistently independent of the growth medium used. The chromatograms and morphological studies showed that P. cinnamomi released specific VOCs at different stages of colony development. Using the optimized HS-SPME GC method, identification of P. cinnamomi from 15 in vivo diseased soil samples was as high as 100%. Results from this study demonstrate the feasibility of this method for identifying P. cinnamomi and the potential use of this method for physiological studies on P. cinnamomi.

Headspace Solid-Phase Microextraction and Gas Chromatography-Mass Spectrometry for Analysis of VOCs Produced by Phytophthora cinnamomi

Volatile organic compounds (VOCs) from Phytophthora cinnamomi-infected lupin seedlings were collected by headspace solid-phase microextraction (HS-SPME). The sampling was done 28 to 44, 52 to 68, and 76 to 92 h after inoculation (HAI). The HS-SPME samples were analyzed by gas chromatography-flame ionization detector (GC-FID) to assess the differences in volatile compounds between the P. cinnamomi-infected lupin seedlings and the control. Three specific peaks were identified after 52 to 68 h with the infected lupin seedlings, at which time there were no visible aboveground symptoms of infection. Subsequently, the VOCs of five different substrates (V8A, PDA, lupin seedlings, soil, and soil + lupin seedlings) infected with P. cinnamomi and the corresponding controls were analyzed by gas chromatography-mass spectrometry (GC/MS). A total of 87 VOCs were identified. Of these, the five most abundant that were unique to all five inoculated substrates included: 4-ethyl-2-methoxyphenol, 4-ethylphenol, butyrolactone, phenylethyl alcohol, and 3-hydroxy-2-butanone. Therefore, these metabolites can be used as markers for the identification of P. cinnamomi in different growing environments. Some VOCs were specific to a particular substrate; for example, 2,4,6-rrimethyl-heptanes, dl-6-methyl-5-hepten-2-ol, dimethyl trisulfide, 6,10-dimethyl- 5,9-undecadien-2-ol, and 2-methoxy-4-vinylphenol were specific to P. cinnamomi + V8A; heptanes and 5-methyl-3-heptaneone were specific to P. cinnamomi + PDA; 3-methyl-1-butanol, ethyl acetate, 2-methyl-propanoic acid, ethyl ester, and ethyl ester 2-methyl-butanoic acid were specific to P. cinnamomi-inoculated lupin seedlings; and benzyl alcohol and 4-ethyl-1, 2-dimethoxybenzene were only detected in the headspace of inoculated soil + lupin seedlings. Results from this investigation have multiple impacts as the volatile organic profiles produced by the pathogen can be utilized as an early warning system to detect the pathogen from contaminated field soil samples. Data from this investigation have also illuminated potential metabolic pathways utilized by the oomycete during infection which may serve as potential targets for the development of specific control strategies.

Impact of primer choice on characterization of orchid mycorrhizal communities using 454 pyrosequencing

Although the number of studies investigating mycorrhizal associations in orchids has increased in recent years, the fungal communities associating with orchids and how they differ between species and sites remain unclear. Recent research has indicated that individual orchid plants may associate with several fungi concurrently, implying that to study mycorrhizal associations in orchids the fungal community should be assessed, rather than the presence of individual dominant fungal species or strains. High-throughput sequencing methods, such as 454 pyrosequencing, are increasingly used as the primary tool for such analyses. However, many studies combine universal primers from previous phylogenetic or ecological studies to generate amplicons suitable for 454 pyrosequencing without first critically evaluating their performance, potentially resulting in biased fungal community descriptions. Here, following in silico primer analysis we evaluated the performance of different combinations of existing PCR primers to characterize orchid mycorrhizal communities using 454 pyrosequencing by analysis of both an artificially assembled community of mycorrhizal fungi isolated from diverse orchid species and root samples from three different orchid species (Anacamptis morio, Ophrys tenthredinifera and Serapias lingua). Our results indicate that primer pairs ITS3/ITS4OF and ITS86F/ITS4, targeting the internal transcribed spacer-2 (ITS-2) region, outperformed other tested primer pairs in terms of number of reads, number of operational taxonomic units recovered from the artificial community and number of different orchid mycorrhizal associating families detected in the orchid samples. Additionally, we show the complementary specificity of both primer pairs, making them highly suitable for tandem use when studying the diversity of orchid mycorrhizal communities.