Molecular Detection of 10 of the Most Unwanted Alien Forest Pathogens in Canada Using Real-Time PCR

Invasive alien tree pathogens can cause significant economic losses as well as large-scale damage to natural ecosystems. Early detection to prevent their establishment and spread is an important approach used by several national plant protection organizations (NPPOs). Molecular detection tools targeting 10 of the most unwanted alien forest pathogens in Canada were developed as part of the TAIGA project (http://taigaforesthealth.com/). Forest pathogens were selected following an independent prioritization. Specific TaqMan real-time PCR detection assays were designed to function under homogeneous conditions so that they may be used in 96- or 384-well plate format arrays for high-throughput testing of large numbers of samples against multiple targets. Assays were validated for 1) specificity, 2) sensitivity, 3) precision, and 4) robustness on environmental samples. All assays were highly specific when evaluated against a panel of pure cultures of target and phylogenetically closely-related species. Sensitivity, evaluated by assessing the limit of detection (with a threshold of 95% of positive samples), was found to be between one and ten target gene region copies. Precision or repeatability of each assay revealed a mean coefficient of variation of 3.4%. All assays successfully allowed detection of target pathogen on positive environmental samples, without any non-specific amplification. These molecular detection tools will allow for rapid and reliable detection of 10 of the most unwanted alien forest pathogens in Canada.


Introduction
Invasive alien tree pathogens can cause significant economic losses as well as large-scale damage to natural ecosystems. Over the last century, Canada has experienced the dramatic consequences of introductions of alien forest pathogens. The pathogens responsible for white pine blister rust (Cronartium ribicola J.C. Fisch), beech bark disease (Cryptococcus fagisuga Lindinger), and Dutch elm disease (Ophiostoma ulmi (Buisman) Melin & Nannf. and O. novoulmi (Brasier)) were accidentally introduced into Canada and resulted in the death of millions of Pinus strobus, P. monticola, P. albicaulis, Fagus grandifolia Ehrh. and Ulmus americana L. trees throughout their distribution range. Despite public and institutional awareness of alien forest species, it is expected that their number and impact will keep increasing in the future [1,2]. In order to implement quarantine and enforce mitigation measures following the introduction of exotic pathogens, national plant protection organizations (NPPOs) such as the Canadian Food and Inspection Agency (CFIA) need rapid, reliable, sensitive and accurate detection methods. The challenge for NPPOs is to be able to detect pathogens at their different life stages, including those that have the capacity to remain latent on asymptomatic tissues. Molecular detection using real-time PCR approaches allows for rapid, reliable and sensitive detection while simultaneously processing of large numbers of samples.
Real-time PCR has become the gold standard in pathogen detection in many fields, e.g. medicine, animal health, agriculture as well as forestry. It is sensitive enough to detect minute amounts of DNA from the target organism mixed with environmental material or host DNA, and the use of hydrolysis probes (TaqMan probes) offers an additional level of specificity, thereby enabling discrimination between closely related species with few polymorphic sites [3]. Real-time PCR also provides accurate quantification of target DNA in processed samples, which is directly proportional to the biomass of the targeted organism. For all these reasons, real-time PCR has been increasingly used to prevent and mitigate the introduction and dispersal of exotic and invasive plant pathogens.
So far, real-time PCR assays of forest pathogens have been mostly developed for single specific pathogens. However, real-time PCR using TaqMan probes offers the opportunity for multiplexing (multiple reactions in one tube [4]) and arraying (multiple reactions in separate tubes but on a single support), allowing for the simultaneous detection of a range of different pathogens in a large number of samples by performing a single real-time PCR run (e.g. [5][6][7]). Multiplexing usually requires extensive fine-tuning to avoid cross-reactivity and/or loss of sensitivity [8,9]. This critical step can be circumvented by using arrays of assays operating under the same real-time PCR conditions, but running in as many tubes as there are assays included in the arrays.
The objective of the present study therefore was to develop and validate a set of sensitive, specific and precise real-time PCR assays for the rapid detection of 10 of the most unwanted alien forest fungal pathogens in Canada selected from a list of over 100 tree pathogens regulated by international, continental, and national phytosanitary organizations. Selection was based on the pathogen's i) history of invasiveness, type and degree of damage/symptoms/pathogenicity, ii) host species and estimated economic impact, iii) dispersal pathways, establishment and adaptability, and iv) likelihood of establishment in Canada. Our goal was to develop assays that could be used either in simplex or in plate-based arrays.

Isolates selection
For each target pathogen, we built a panel of isolates encompassing multiple isolates of the target species and isolates from closely related species (sister species). Our selection of sister species was based on phylogenies found in recent scientific peer-reviewed studies [10][11][12][13][14][15] as well as on the advice of particular taxonomic group specialists. Cultures were obtained from collections (CBS, ATCC, as well as private ones) and stored in replicates at FPInnovations in Vancouver and the CFIA in Ottawa. When available, type cultures from referenced culture collections along with isolates provided by taxonomic authorities were preferred. To capture intraspecific genetic diversity, isolates from different hosts and different geographic origins were used when available. The list of isolates used to design the assays is presented in Table 1.

DNA extraction
For all isolates (except for Ceratocystis species), DNA was extracted using the Qiagen's DNeasy Plant Mini Kit (Qiagen, Valencia, CA, USA) according to the manufacturer's instructions. Ceratocystis DNA was extracted using a modified version of Zolan and Pukkila's phenol/chloroform extraction protocol [16]. A small piece of mycelium was homogenized in 400 μl of extraction buffer (100mM Tris-HCl pH 9.5, 1.4 M NaCl, 20 mM EDTA, 2% CTAB, 1% PEG 8000, and 0.25% β-mercaptoethanol). Samples were then incubated at 65°C for 1h (vortexing every 15 minutes). Next, 400 μl of phenol:chloroform:isoamyl alcohol (25:24:1) were added to the homogenate, vortexed for 10 s and centrifuged at 13,000 x g for 10 min. Supernatant was mixed by inversion with 70 μl of 7.5M ammonium acetate and 600 μl of ice-cold isopropanol, incubated at -20°C for a minimum of 1h, and then centrifuged at 13,000 x g for 20 min (4°C). DNA was rinsed with 800 μl ice-cold 70% ethanol and centrifuged at 13,000 x g for 5 min (4°C). DNA was then incubated at 55°C to evaporate any remaining ethanol and re-suspended in 50 μl 10 mM Tris-HCl, pH 8.0. DNA was visualized on agarose gel stained with ethidium bromide, and DNA concentration was measured using the Qubit 2.0 Fluorometer with the dsDNA BR Assay Kit (Life Technologies, Carlsbad, CA, USA) according to the manufacturer's instructions.

DNA sequencing and phylogenetic analyses
The internal transcribed spacer (ITS) gene, recognized as the universal DNA barcode for fungi [17], was systematically amplified and sequenced for all isolates. NCBI nucleotide blast of the ITS sequences was performed to detect misidentification and potential contamination of isolates. A list of the different gene regions sequenced along with the primers used is presented in Table 2. PCR reactions were performed in a final volume of 25 μl and contained 1X PCR buffer, 1.5 mM MgCl 2 , 200 μM of each dNTP (Invitrogen), 0.4 μM of each primer (Integrated DNA Technologies Inc., Coralville, IA, USA), 1 U of Platinum Taq DNA polymerase (Invitrogen), and 1 μl of template DNA. Sequencing of both DNA strands was performed by the Centre de recherche du Centre Hospitalier Universitaire de Québec (CHUQ) sequencing platform on an ABI 3730xl (Applied Biosystems, Foster City, CA, USA) using the specific forward and reverse primers.
Alignments were used to guide the development of assays and were performed with the ClustalW algorithm implemented in BioEdit v7.1.3.0 [18]. Evolutionary relationships between targets and their sister species were inferred from DNA sequences of ITS. Phylogenetic trees were reconstructed by using the maximum likelihood method with the Tamura-Nei model implemented in MEGA5 [19]. Statistical support of nodes was assessed by performing 500 bootstrap replicates.
For the Phytophthora ramorum and Phytophthora kernoviae targets, we did not perform the gene region sequencing and phylogenetic analyses described above. Instead, the detection assays for these two species were designed in genes unique to these target species that were identified by using a comparative genomics approach developed in the TAIGA project (http:// taigaforesthealth.com/Home.aspx) (S1 File).
SYBRGreen-based real-time PCR quantification for standardization of isolates' DNA concentration DNA concentration of all isolates was standardized following a qPCR quantification using genus general primers. To do so, we quantified the number of target gene copies that were initially present (before any PCR amplification) in the sample, which directly relates to the abundance of the pathogen prior to DNA extraction. This quantification allowed us to work with samples having a standardized DNA concentration for specificity validation. It assured us we had DNA in high enough concentration in all samples to confirm assay discrimination against closely related species. Genus general primers were designed using Oligo Explorer v1.2 and Oligo Analyzer v1.2 (Gene Link, NY, USA) in a conserved gene region for all closely related species. The following criteria were also used to guide primer design: 1) length between 18 and 25 bp; 2) melting temperature (T m ) close to 60°C (using the nearest neighbor algorithm); 3) absence of polymorphism within targeted species; and 4) minimal secondary structure (especially dimer formation at the 3' end). Primer pairs were designed such that PCR products were shorter than 200 bp ( Table 2). Real-time PCR was performed with an Applied Biosystems 7500 Fast Real-Time PCR System (Life Technologies, Carlsbad, CA, USA). All reactions were performed in a final volume of 10 μl and contained 1X QuantiTect SYBR Green PCR Master Mix (Qiagen, Valencia, CA, USA), 0.5 μM of each of the genus general primers (Table 2), and 1 μl of template DNA. Real-time PCR thermocycling conditions were set at 95°C for 15 min, followed by 50 cycles at 95°C for 15 s, 58°C (primer Tm-2°C) for 30 s, and 65°C for 90 s. Fluorescence was read at the end of the extension step.
Gene copy number quantification was then performed using a Java program based on linear regression of efficiency [20] and sample DNA concentration was adjusted to 5,000 gene copies per μl, whenever possible.

Target-specific TaqMan-based real-time PCR assays
All the molecular detection assays targeting prioritized tree pathogens are based on the Taq-Man technology. The following strategies were used to design all of our detection assays. Based on the sequences recovered, we targeted the gene that allowed for the best discrimination at the species level, i.e. the gene that maximized the number of single nucleotide polymorphisms (SNPs) between species while keeping a low level of intraspecific variability. Primer and probe designs were performed using Oligo Explorer v1.2 and Oligo Analyzer v1.2. Each set of primer pair and probe was designed so that there was minimal secondary structure (especially dimer formation at the 3' end) and amplicon length did not exceed 350 base pairs (Table 3). Primers and probes were also designed to ascertain that interspecific SNPs were preferentially localized at the 3' end of the primers for maximum discrimination effect of the primer-template annealing [21]. The real-time PCR master mix used, QuantiTect Multiplex PCR NoROX Master Mix (Qiagen), possesses features that allow for the use of short oligonucleotides when necessary. By allowing the design of shorter primers and probes, these elements increase the SNP specificity of the primer and probe. All probes were labelled with fluorescein (6-FAM) at the 5' end and with the quencher Iowa Black FQ (ZEN-IBFQ). All primers and TaqMan probes were manufactured by Integrated DNA Technologies Inc. All assays were designed to work under the same thermocycling conditions, offering the opportunity to array them into 96-or 384-well plates machine formats, based on the user's needs.
Validating the specificity of the tree pathogen TaqMan assays Specificity validation of all the assays was performed using the panels of isolates presented in Table 1 and Fig 1. For target species belonging to same genera (Ceratocystis and Phytophthora), we used the whole genera panel to evaluate each of the assay's specificity. Real-time PCR amplification was conducted using 1X QuantiTect Multiplex PCR NoROX Master Mix, with 0.6 μM of each primer, 0.1 μM of TaqMan probe, and 5,000 gene copies of template DNA, whenever possible, in a final reaction volume of 10 μl. Two technical replicates were performed for all reactions using an Applied Biosystems 7500 Fast Real-Time PCR System. Real-time PCR thermocycling conditions were set at 95°C for 15 min, followed by 50 cycles at 95°C for 15 s and 60°C for 90 s. Fluorescence was read at each cycle, at the end of the extension step.
Validating the sensitivity of the tree pathogen TaqMan assays Sensitivity of the TaqMan assays was evaluated in terms of both efficiency and limit of detection (LOD). For each target assay, experiments were conducted to 1) determine if Ct values were proportional to the amount of target template DNA (efficiency) and 2) evaluate the LOD, which is the smallest amount of target DNA that can be detected for each of the assays. One isolate for each of the target species was selected, and TaqMan assay sensitivity was assessed on parallel sets of serial dilutions from the DNA stock.
To assess efficiency of the amplification reaction, TaqMan assays were run with serial dilutions of template DNA from the target species, ranging from 1 to 15,000 copies of the target gene region, as quantified using the species-specific primers. Standard curves were obtained by plotting the values of Ct against the log value of the target gene region copy number. Amplification reaction efficiency was calculated using the following formula: where E represents the amplification reaction efficiency and slope is the slope value of the line derived from the standard curve plot. Estimation of the LOD was done by performing 20 replicates of the TaqMan real-time PCR reactions for each of the following DNA concentrations: 1, 3, 5, and 10 copies per μL. The lowest DNA concentration with a level of 95% successful amplification was identified as the LOD. Validating the precision of the tree pathogen TaqMan assays Precision (or repeatability) of the assays refers to the robustness of the assay with the same samples repeatedly analyzed in the same manner [24]. Ct values from different real-time PCR runs on different isolates of target species, assessed with a standardized concentration of 5,000 gene copies, were compiled and used to determine the precision of the assays. For each assay, mean Ct value, standard deviation and coefficient of variation were calculated.

Validating the tree pathogen TaqMan assays on environmental samples
The complete list of all environmental samples, including the source, is presented in Table 4. Because of the phytosanitary risks of infected material, environmental samples were supplied by collaborators as purified DNA samples. Since the objective was to test the assays' performance in a variety of different conditions, collaborators were free to use the routine DNA extraction protocols implemented in their respective laboratories instead of a unique standardized DNA extraction protocol. The efficiency of the DNA extraction was assessed for each sample by performing a control real-time TaqMan PCR reaction that targeted either the plant cytochrome oxidase gene (for primers and probes sequences, see Table 3). All reactions were performed in a final volume of 10 μL and contained 1X QuantiTect Multiplex PCR NoROX Master Mix, with 0.6 μM of each primer, 0.1 μM of TaqMan probe, and 1 μL of template DNA. Real-time PCR thermocycling conditions were set at 95°C for 15 min, followed by 50 cycles at 95°C for 15 s and 60°C for 90 s. Fluorescence was read at the end of each cycle. TaqMan assays were performed with three technical replicates for each related environmental sample. Reactions were performed as described earlier, using 1 μL of environmental DNA.        Positive (using target species' DNA from pure culture) and negative (no template DNA) controls were included in all qPCR runs. Target gene region copy numbers were calculated by translating Ct values, using standard curve equations. Positive results were all confirmed by Sanger sequencing of the real-time TaqMan PCR product.

Results and Discussion
Assay design and development Development of the detection assays was based on two strategies targeting 1) unique SNPs or 2) unique genes. The SNP-based approach uses alignment of genes present in all species, but exploits the presence of SNPs between the target species and the close relatives. It was used for all assays reported in this paper except for P. ramorum and P. kernoviae. For both of these species, comparative genomics was used to identify genes uniquely found in the target species to design the detection assay. Detail about this TAIGA strategy, the related genomic resources and bioinformatics pipeline are available on the TAIGA project website (http:// taigaforesthealth.com/). A crucial step in the development of the SNP-based detection assays is the identification of appropriate target DNA regions. Readily amplified genes across taxa of a group were sequenced, and genes showing interspecific variability were selected for assay development ( Table 2). As a result, different genes were selected for each target species, some of them being single-or low-copy genes (e.g. β-tubulin, EF1, RPB2 and Hsp60), and others being multi-copy genes (e.g. IGS, ITS, Cluster62 and Cluster97).
In order to standardize the DNA concentration of all isolates, a genus general real-time PCR SYBRGreen I assay targeting the selected DNA region was designed and validated for each target group. Using the linear regression of efficiency (LRE) quantification approach [20], DNA concentration of isolates was determined and standardized to 5,000 target gene region copies, which usually translates into a Ct value ranging between 20 and 25.
To ensure the repeatability of the qPCR experiments described by other teams and to ease the interpretation of results, most of the Minimum Information for Quantitative Real-Time PCR Experiments (MIQE), as described by Bustin et al. (2009) [24], is presented in this paper.

Specificity of the tree pathogen TaqMan assays
Probes' and primers' specificity was first tested in silico using BLAST on the NCBI nucleotide collection (nr/nt) database. A wet lab was then performed to assess candidate sets' panel specificity on DNA samples from target and sister species isolates. During the first round of specificity validation, whenever an unexpected amplification was observed, two hypotheses were explored. First, it could be due to trace contamination of the DNA sample with target species' DNA, which sometimes happens during sample manipulation. This was suspected when cycle threshold (Ct) values were much higher (around 35-37) than those of the target species isolates. When contamination was suspected, a SYBRGreen real-time PCR reaction along with a melting curve gradient was performed. When melting curves of positive and suspected false positive samples were identical, the SYBRGreen real-time PCR reaction product of the false positive was sequenced and aligned with reference sequences to confirm the contamination. In such cases, the contaminated DNA sample was discarded and fresh DNA was re-extracted from a pure culture of the isolate. If the first hypothesis was confirmed, we concluded we had a real false positive reaction due to a lack of specificity. In such cases, further screening of primers and probes was conducted.
All our final assays were 100% specific successfully discriminating the target species of tree pathogen from the closely related species. In cases where we were unable to obtain culture or DNA of some closely related species, we still performed in silico specificity validation using sequences obtained from the public domain. Despite the current results, we cannot rule out potential cross reactivity of the present assays with evolutionarily related species that have not been described yet. Some of the target species belong to what can be considered as orphan and poorly resolved taxonomic groups (such as Mycosphaerella, with newly described species M. musivoides P.E. Busby & G. Newc and M. wasatchii P.E. Busby & G. Newc [25]). Instead of being a drawback, molecular cross reactivity with cryptic species can represent an opportunity to isolate and describe novel fungal species with similar or different pathogenic behaviors [26].

Sensitivity of the tree pathogen TaqMan assay
Overall, the assays we developed have high efficiency and sensitivity, with limits of detection varying between 1 and 10 target gene region copies (Table 5). In the present study, no difference in sensitivity values was observed between assays targeting single-or low-copy genes and those targeting multi-copy genes. For all tested species, Ct values were proportional to the amount of template DNA used for the real-time PCR reaction. The standard curves generated by plotting the log of DNA (copies) against the Ct value determined by qPCR display linearity across the whole range of dilutions assessed, with a correlation coefficient (r 2 ) ranging from 0.950 for the P. kernoviae assay to 0.999 for the G. morbida assay (Fig 2). Moreover, PCR amplification efficiencies ranged between 83 and 97%, which is considered to be an acceptable range [27], except for P. kernoviae. The low amplification efficiency (73%) of the P. kernoviae assay could be due to a number of factors, such as the presence of inhibitors in the DNA samples, suboptimal primer and probe design (e.g. presence of non-specific products and primer dimers), and pipetting errors. Experimental investigations ruled out the presence of non-specific products and dimers, sample contamination, inappropriate dilution series and pipetting errors. Another possible explanation for this reduced efficiency is the presence of a secondary structure in the region targeted by this assay. Amplification efficiency can vary across a genome [28,29]. Genomic regions resistant to amplification by PCR correlate with high GC contents [30,31] that do not denature efficiently under routine amplification conditions. However, the GC content of that region was around 55%, which is not considered a problem.
The limit of detection (LOD) was defined as the lowest concentration of target DNA at which 95% of the positive samples were detected [24]. According to Bustin et al. (2009), and assuming an even Poisson distribution, the lowest theoretically possible detection limit is three DNA copies per PCR reaction to provide a positive signal in 95% of the PCR reactions performed. Our assays revealed LOD values varying between 1 and 10 target gene region copies (Table 5). Those values are comparable to what has been reported by others working on trace detection for regulatory or public health applications, looking either for the presence of genetically modified DNA [32][33][34][35], virus DNA in blood samples [36,37], or human DNA [38,39], where LOD values varying between 1 and 25 copies were obtained.
As an additional validation step, we were able to compare results from our F. circinatum assay with those published for a similar qPCR test developed by Ioos et al. (2009) [40]. The published assay was not compliant with our set of real-time PCR conditions and had to be redesigned to suit qPCR standardized parameters. Our assay targets a different segment of the intergenic spacer region than the one used elsewhere [40]. To conduct a fair comparison, a subset of the environmental samples used by Ioos et al. (2009) [40] was obtained and tested with our assay. The results we obtained were actually very similar. Although our assay had a slightly delayed detection threshold of approximately 3 Ct values compared with that of Ioos et al.

Precision of the tree pathogen TaqMan assays
The precision of our ten tree pathogen TaqMan assays is shown in Fig 3. All assays have a mean Ct value ranging between 23 and 26 for 5,000 target gene region copies. This value depends on amplicon size and primers and probe properties. Using the mean Ct value and the standard deviation, we also calculated a coefficient of variation for each assay, which varied between 0.7% for the S. pseudotuberosa assay and 8.2% for the F. circinatum assay. Those values clearly demonstrate that our assays have a high degree of repeatability, an important advantage when dealing with possible regulatory issues. Validation of the tree pathogen TaqMan assays on environmental samples All tree pathogen assays successfully detected target pathogens from the positive environmental samples provided by collaborators (Tables 4 and 6). Negative environmental DNA samples were available for 6 out of the 10 assays; no false positive results were obtained with any of these.
One of the most critical steps when dealing with environmental samples is the quality of the DNA sample, which may vary with the extraction protocol used. In spite of this, detection limits as low as one target gene region copy were obtained (Table 6). For some samples, we obtained a positive result, i.e. a detectable Ct value using the TaqMan specific assays, but the calculated target gene region copy number was less than one. These results might be explained by the fact that we used standard curve equations to extrapolate those copy number values. Therefore, there exists a certain level of imprecision that has a more important effect on samples with a low level of target pathogen DNA. This imprecision caused by the extrapolation of target gene region copy values can also explain why, in some other cases (e.g. C. fagacearum, G. morbida, F. circinatum and S. pseudotuberosa), we obtained a slightly higher value for the target gene region than the one obtained with the genus assay ( Table 6). The opposite result was also seen; for some environmental samples (e.g. C. laricicola, C. polonica, G. abietina (EU race), R. necatrix, S. pseudotuberosa and P. ramorum), the genus gene copy number was much higher    than the target gene region copy number. This is most probably due to the presence of more than one species from the targeted genus in the environmental samples. For example, G. abietina (EU race) positive environmental samples were obtained from infected Pinus resinosa samples from the province of Québec, Canada. We know that the North American race of G. abietina, which is detected and counted with the genus assay but not with the EU race specific assay, might also be present on red pines [41]. Other examples are Rhododendon sp., Griselinia sp., and Magnolia x soulangeana plant tissues infected with P. ramorum. Those plants are well known to be hosts for other Phytophthora species [42,43], which might explain the high level of quantification at the genus level, concomitant with a low level of P. ramorum itself. Inhibitors from plant material or insect specimens co-extracted with DNA are a source of contamination that can impact PCR amplification accuracy [44][45][46], which is the variation between observed and expected data [24]. This property was evaluated for some of those assays (C. polonica and C. laricola) in a previous study [47]. In those specific cases, the presence of environmental DNA or any other co-extracted compound had no effect on the performance of the TaqMan assays. However, we are aware that assay performance could vary slightly depending on the material it is tested against. In fact, inhibition may be caused by the presence of different compounds such as acidic plant polysaccharides [48,49], plant phenolics [50], the contamination of DNA samples with co-extracted polyphenol-bound proteins from the insect cuticle [51], or with phenolics and tannins found in the digestive tracts of xylophagous insects [45]. Accuracy evaluation should therefore be one of the initial steps for any user of those assays dealing with new environmental material.
We developed sensitive and specific molecular assays for ten alien tree pathogens identified as high priority potential threats for Canadian forests: Ceratocystis fagacearum, Ceratocystis laricicola, Ceratocystis polonica, Fusarium circinatum, Gremmeniella abietina (EU race), Geosmithia morbida, Phytophthora kernoviae, Phytophthora ramorum, Rosellinia necatrix and Sclerotinia pseudotuberosa. All of these assays are specific, i.e. they have the ability to amplify a unique DNA fragment of interest without amplifying or detecting non-target sequences. Detection assays were already available for some of the target tree pathogens selected. In some cases, they were included in our tree pathogen TaqMan assay panel (e.g. C. polonica and C. laricicola [47]). However, in most cases, existing assays had to be redesigned either because 1) they were not compliant with our real-time PCR conditions (C. fagacearum [52], F. circinatum [40,53], P. kernoviae [54], P. ramorum [9,[54][55][56][57][58][59], S. pseudotuberosa [60]), 2) they were not tested against all closely related species (F. circinatum [61], P. kernoviae [7,62], P. ramorum [7,55,[63][64][65], R. necatrix [66][67][68]), or 3) they did not target the specific race of interest (G. abietina EU race [69]). All assays were designed to be used under the same real-time PCR conditions, using the same chemistry and the same thermocycling parameters. Therefore, they can be performed in micro-well plates arrayed in any machine format to suit individual users' needs and to increase throughput. Reactions for multiple samples, targeting multiple pathogens, can be performed in a single real-time PCR run, which is an important advantage under operational conditions where testing a large number of samples against of large number of targets is required. Molecular detection of these pathogenic species directly from plant material or insect vectors represents a powerful tool to prevent their introduction and establishment as potential invasive species.
Supporting Information S1 File. Identification of unique gene models to Phytophthora kernoviae and Phytophthora ramorum. (DOCX)