1. Introduction
Phytophthora, a genus prominent among oomycetes, encompasses over 150 species, with the majority acting as plant pathogens causing collar and root diseases in various plant species. These diseases often lead to the decline and eventual demise of the afflicted plants [
1].
Phytophthora cinnamomi is the causal agent of several epidemics across Europe and globally and it is considered one of the most devastating plant pathogens in the world [
2]. Its impact spans forestry, horticulture, and nursery industries [
3,
4]. With a host range encompassing approximately 5,000 woody plant species across 70 countries,
P. cinnamomi poses a substantial threat [
5,
6,
7,
8].
Although initially prevalent in tropical and subtropical regions,
P. cinnamomi has demonstrated adaptability to cooler and drier environments [
6,
7]. Its ability to thrive saprophytically in soil or persist asymptomatically in non-host plants significantly contributes to its long-term survival.
P. cinnamomi is reported particularly damaging to oaks [
9,
10], to chestnuts where it causes ink disease [
11] and common walnuts (
Juglans regia L.). This pathogen has led to substantial economic losses in all walnut growing regions like southern Europe, the USA, and Chile, [
5,
7,
12,
13,
14,
15]. Italy alone has witnessed the uprooting of over 150 hectares of common walnut orchards due to
P. cinnamomi attacks [
16].
Efforts to combat
P. cinnamomi include research on resistant walnut rootstocks [
17], the exploration of control strategies such as phosphite application [
12,
18] and containment and/or eradication of spot infections [
19]. Early detection remains pivotal in managing P
. cinnamomi diseases, especially in tree cultivation where plant diagnosis often occurs late. A timely diagnosis, in these cases, could be performed identifying the pathogen directly in soil surrounding plants, thus, allowing prevention and control of primary sources of infection.
The attention surrounding this pathogen is demonstrated by the number of identification methods developed over the years, starting from the first PCRs in 2003 up to sniffer dogs in 2023 [
20]. Molecular genetic assays designed for
P. cinnamomi detection, include polymerase chain reaction (PCR) [
21], nested [
22] and real-time PCR [
23,
24,
25], loop-mediated isothermal amplification (LAMP) and Recombinase Polymerase Amplification Assay (RPA) [
26,
27]. However, many lack specificity, [
28] and only a few based on real-time PCR capable of quantifying pathogens [
24,
25].
Naturally infested soils were not quantified in these previous studies; one soil type was usually studied following inoculation with the pathogens. Thus, there are no reports on how P. cinnamomi population in soil fluctuate in actual production fields.
The selection of diagnostic methods is typically driven by application needs, it depends on specificity and/or sensibility required, on necessity to have field application and/or pathogen quantification, moreover a suitable DNA extraction should be implemented, especially for environmental samples. For the correct management of tree cultivation, like walnut orchards in Northern Italy, we need a new P. cinnamomi early detection method especially efficient for soil samples analysis, specific enough to distinguish P. cinnamomi from other Phytophthora present in soil, sensible enough to detect the pathogen in very low amount typical of soil samples and with the possibility to quantify the pathogen biomass to know the level of infection.
In this context, the choice of genetic locus and in silico sequence analysis are crucial. While the internal transcribed spacer (ITS) of nuclear DNA is commonly used, it may not effectively discriminate between closely related Phytophthora species. For diploid organisms like Phytophthora, nuclear single-copy genes, like the ras-related Ypt1 gene [
29], offer advantages for phylogenetic analysis and species discrimination. Moreover, this kind of genes, devoid of intergenomic concerted evolution, could provide accurate quantification of pathogen biomass, when analyzed in quantitative real-time PCR.
Successful molecular diagnostics hinge on obtaining sufficient, high-quality DNA from samples. Soil DNA extraction is often challenged by the low amount of DNA of single organism present in soil and by the presence of polymerase inhibitors such as polyphenols, polysaccharides, and humic acids [
30].
Overcoming these obstacles, a highly sensitive and specific qPCR methodology based on the ras-related Ypt1 gene was developed. The objectives of this study were to (i) develop a qPCR assay for the sensitive and specific detection and quantification of P. cinnamomi, (ii) adapt the assay to the analysis of soil samples by improving DNA extraction method and sampling procedure for early detection of P. cinnamomi in field, (iii) validate the qPCR assays in walnut orchards infected by P. cinnamomi.
2. Materials and Methods
2.1. Phytophthora Isolates and Isolation
The isolates used in this study are listed in
Table 1. A total of 50
Phytophthora isolates belonging to clades 1, 2 ,4, 6, 7, 8, 9, 10 plus two
Pythium species (
Table 1) were used to test the specificity and sensitivity of the primers set developed for
P. cinnamomi diagnosis. Several
P. cinnamomi isolates obtained from woody plants (walnut, oak, or chestnut) were used (
Table 1). DNA of the
P. cinnamomi isolate CREADC-Om274 from walnut was used as reference to set up real-time PCR conditions.
P. cinnamomi isolates were isolated from symptomatic plants from several Italian regions. Most of them were obtained from naturally infected common walnut (Juglans regia L.) trees grown in commercially fruit orchards of northern Italy. Stock cultures are preserved in the dark both in sterile distilled water at room temperature and at 10 °C, as well as in oatmeal agar (OA - Sigma-Aldrich, Saint Louis, MI, USA) slant tubes with mineral oil at 15 ± 2°Cat CREA-DC culture collection in Rome. Pure cultures were obtained by transferring single hyphal tip from the edge of the colonies onto potato dextrose agar (PDA) (Oxoid, Basingstoke, UK).
Tissue fragments were obtained from collars or stems of infected woody material cut from the margins of necrotic lesions. Tissue fragments obtained from healthy plants were used as control. Small tissue fragments of about 3-5 mm × 3-5 mm previously surface disinfested for 1 min in a 1 % NaOCl solution, rinsed for 5 min in sterile distilled water, were either plated onto P5ARPH selective medium [
31] or placed in a 1.5 ml Eppendorf ® tube for DNA extraction and stored at -20 °C until use.
2.2. DNA Extraction from Pure Colony and Plant Tissue
For all Phytophthora and Pythium isolates used in this study, mycelial DNA was extracted from pure cultures grown on PDA at 25 °C for 5 days in the dark. Mycelium was scraped and ground to a fine powder under liquid nitrogen, placed in a 1.5 ml sterile Eppendorf® tube, and stored at -20 °C until use. Total DNA was extracted using Wizard genomic DNA purification kit (Promega, Madison, WI, USA) following manufacturer’s instructions.
DNA extraction from plant tissue was performed on approximately 100 mg of tissue fragments. Samples were homogenized by grinding in liquid nitrogen and total DNA was extracted using the DNeasy Plant Mini kit (QIAGEN GmbH, Hilden, Germany) following the manufacturer’s instructions.
2.3. Real-Time PCR Primer and Probe Design for P. cinnamomi
The single copy ras-related protein gene
Ypt1 was chosen to design primers and probe. For this purpose, 51
Ypt1 gene sequences of
Phytophthora spp. representing all ten clades of Phytophthora philogeny and one Ypt
1 sequence of
Pythium aphanidermatum, were retrieved from NCBI GenBank Database (
Table S1) and aligned using the multiple sequence comparison by log-expectation (MUSCLE) method [
32], to find the unique polymorphic regions of
P. cinnamomi. Based on these regions, primers and probes were designed by PRIMER3 0.4.0 (
http://frodo.wi.mit.edu/primer3/) [
33].The probe was labelled at the 5′ end with 6-carboxyfluorescein (FAM) as a reporter dye and modified at the 3′ end with the quencher Black Hole Quencher1 (BHQ1). Primers and probe melting temperatures (TM), were calculated with PRIMER 3 software [
34].
2.4. qPCR Conditions
All Real-time PCR reactions were run in MultiplateTM PCR Plates 96-well clear (Bio-Rad, UK), using a CFX96 C1000 Thermal Cycler Real-Time System (Bio-Rad- Hercules, California, USA). Data acquisition and analysis were obtained by the supplied Bio-Rad CFX Manager software version 3.0 (3.0.1224.1015) according to the manufacturer’s instructions.
Each 15-μl reaction contained 1μl of genomic DNA,1x GoTaq® G2 Hot Start Buffer (Promega), 5 mM MgCl2, 0.2 mM each dNTP, 0.33 μM of each primer, 0.13 μM of the probe, and GoTaq® G2 Hot Start DNA Polymerase (Promega). Negative control reactions contained 1µl of sterile distilled water. Reactions were performed under the following conditions: 10 min at 95 °C, followed by 40 cycles at 95 °C for 20 s, and 62 °C for 20 s. Fluorescence was monitored in each PCR cycle during the annealing–extension phase at 62 °C. The cycle threshold (Ct) value was calculated automatically using software version 3.0 (3.0.1224.1015) by determining the PCR cycle number at which the reporter fluorescence exceeded the background. Triplicate reactions were performed in each assay, and each assay was repeated at least twice.
Nested real-time PCR conditions were as follow: first-round PCR was performed with primers YPh1-fwd and YPh1-rev (Table 3) for
Phytophthora spp. amplification in conventional PCR following conditions described by Schena [
21]. The second round was carried out with 1μl of the first-round PCR product as template in real-time PCR using the primer pair and the probe developed in this study at the conditions described above.
2.5. Validation of the Real-TimePCR Method
To assess the analytical sensitivity a log-linear standard curve was generated with the following concentrations of P. cinnamomi isolate CREADC-Om274 genomic DNA: 5, 2, 1 ng/μl, 500, 200, 100, 50, 20, 10, 5, 2, 1 pg/μl, 500, 200, 100, 50 fg/μl, by plotting logarithms of known concentrations of target DNA against the Ct values, considering three replicates for each concentration level. The resulting regression equations were used to calculate P. cinnamomi DNA amount in unknown samples.
The limit of detection (LOD, expressed in ng) was determined as the lowest amount of target genomic DNA that is amplified in 100% of the replicates. Linearity of the method was evaluated on three different P. cinnamomi isolates: CREADC-Om274 from walnut, CREADC-Om139 from chestnut, and CREADC-Om144 from oak.
To determine the Ct Cut-off value, i.e. the Ct above which signals are considered negative, we have analyzed serial dilutions, i.e. 10, 1 ng/μl, 100, 10, 1 pg/μl, 100, 50, 10 fg/μl, of P. cinnamomi isolate CREADC-Om274 genomic DNA, each concentration (group) with five replicates. The cut-off cycle was obtained from the mean Ct of the last group of samples with at least 3 replicates positive for DNA concentration out of five, plus 0.5 (to consider the difference in threshold chosen between runs).
The analytical specificity of the qPCR was tested using 1ng/μl DNA gDNA from the 50
Phytophthora spp. and the two
Pythium spp. isolates listed in
Table 1. Prior to specificity test, all DNA samples were subjected to a conventional PCR with primers ITS6 and ITS4 according to Cooke [
35] to check their ability to be amplified. For pathogen quantification, a standard curve with tenfold dilutions from 1 ng/μl to 100 fg/μl of
P. cinnamomi isolate CREADC-Om274 genomic DNA is performed.
2.6. Soil DNA Extraction and Sampling
For soil DNA extraction from pot plants, Quick-DNATM Fecal/Soil Microbe Midi Prep Kit (Zymo Research, Irvine, CA, USA) on 3-5 g of soil was used following manufacturer’s instructions. For field soil samples the following modifications were made: 15g of soil for each sample were placed in a 50 ml Falcon® containing Bashing BeadsTM, 100 µl of Proteinase K (20 mg/ml; Sigma, Saint Louis, MO, USA), and 27 ml of lysis buffer (0.1 M Tris-HCl pH 8, 0.1 M EDTA, 0.1 M Na2HPO4, 1.5 M NaCl, 1% CTAB (hexadecyltrimethylammonium bromide; pH 8). Homogenization was performed by Fast-Prep 24 5G (MP Biomedicals, USA) at the speed of 6.0 m/s, using Adapter BigPrep for 50 ml Falcon®, for 40 s. Then, 6 ml of 10% SDS were added, and samples were incubated for 2 hours at 65 °C followed by centrifugation at 6,000 × g for 10 min. The supernatant was taken and filtered with the Zymo Spin V-E columns, provided by the Quick-DNATM Fecal/Soil Microbe Midi Prep Kit (Zymo Research, Irvine, CA, USA). Subsequent DNA purification steps were performed following Zymo Research kit manufacturer’s instructions. Total DNA was suspended in 150 µl DNA Elution Buffer. This volume was concentrated to 45 µl by ethanol precipitation. Soil DNA extractions was performed in 3 replicates per sample. After the extraction procedures, the concentration and quality of DNA was checked by Qubit with the dsDNA HS Assay Kit (Invitrogen, Carlsbad, CA, USA) and by Nanodrop ND-1000 (Thermo Fisher scientific Inc., USA). DNA was stored at -20 °C until use.
To test the effect of the matrix on the amplification efficiency, DNA extraction of the different soils analysed were spiked with known concentrations of P. cinnamomi isolate CREADC-Om274 genomic DNA, namely from 1 ng/μl to 100 fg/μl, as standard curves.
Sampling: Soil samples (250 g) were collected as shown in
Figure 1. In the orchard, walnut trees are arranged in rows, at 4-5 m from each other, with approximately 7-8 m between one row and the other. Soil samples were collected at 0,5 m, 1,5 m, 2,5 m, and 3,5 m from the tree along the row at the two different depths: 20 cm and 40 cm (
Figure 1). Soil samples were also collected between one row and the other (inter-row) at 1,5m from the tree at the two depths of 20 cm and 40 cm.
Non-dried soil samples were extracted upon arrival to the laboratory or after storage at -80 °C. Diseased plant material as well as soil samples were collected in spring and fall.
4. Discussion
Soilborne diseases are a major limiting factor for the cultivation of most crops and are costly and difficult to manage. Practical and economic methods of disease control are limited once a crop has been established. Detecting plant pathogens in soil is a critical aspect of plant disease management and agricultural sustainability. Plant pathogens, including bacteria, fungi, oomycetes, nematodes, and viruses, can reside in the soil and pose significant threats to crops by causing diseases that reduce yield and quality. Early detection and accurate identification of soil-borne pathogens result thus essential for implementing effective control measures and minimizing economic losses.
P. cinnamomi, is a very invasive soil-borne plant pathogen affecting thousands of known hosts, including ornamental plants, horticultural and tree crops, and natural ecosystems. In tree crops like walnut or chestnut, the decline caused by
P. cinnamomi starts from roots often remaining unnoticed for years. When symptoms like wilting, yellowing, and retention of dried foliage become visible on the tree crown, rot of feeder roots, larger roots, crowns, and stems [
5] are so severe to be fatal. Plant detection analysis at the crown and stems level is destructive, often tardive and can exacerbate disease spread, resulting in significant economic losses.
Early detection of P.
cinnamomi in soil is a promising approach to managing infected tree orchards. However, detecting plant pathogens in soil presents several challenges
: i) Pathogen Diversity: The soil is a complex environment with diverse microbial communities. Phytophthora populations in soils exhibit significant variability in terms of diversity and species richness due to the polyphagy of these species and their ability to adapt to different environments, [
36]. Distinguishing
P. cinnamomi from other Phytophthora species requires highly specific methods;
ii) Low Abundance: Pathogens are often present in low numbers, making detection difficult. Sensitive techniques like real-time PCR or Next Generation Sequencing are essential for identifying these low-abundance pathogens;
iii) Soil Interference: Soil components can inhibit detection methods. Sample preparation and DNA extraction techniques must minimize these interferences to ensure accurate results;
iv) Environmental Variability: Soil properties such as pH, moisture, and organic matter content can affect pathogen survival and detection. Standardizing sampling and analysis protocols is crucial for reliable detection.
To address these challenges, we developed a detection procedure, including a new qPCR assay, for early detection of P. cinnamomi in field soil samples from tree crop orchards.
To assure specificity, the qPCR developed here has been based on the ras-related protein gene
Ypt1 [
37]. which exhibits high interspecific genetic variability and low intraspecific variation, facilitating diagnostic assays [
37].
In addition, this gene is well represented in GenBank Database providing reference material for
Phytophthora species identification. In the present study, the alignment of 51
Ypt1 gene sequences belonging to 48 different
Phytophthora species and two
Pythium species (
Table 2) has shown that the intron 3 region of
Ypt1 gene is the most suitable one for designing primers and probes that specifically identify
P. cinnamomi. Specificity was tested
in vitro against 50 isolates belonging to 21
Phytophthora species of clades 1, 2 ,4, 6, 7, 8, 9, and 10, including 10 isolates of
P. cinnamomi (
Table 1).[
28] underlined the importance to directly test the efficiency and specificity of the primers against the most closely related species. Particular attention was given to
Phytophthora species of clade 7, with emphasis on the same subclade (7c) of
P. cinnamomi such as
P. parvispora, to avoid false positives.
The qPCR method’s lower limit of detection (LOD) is 200 fg of total genomic DNA, comparable to other real-time diagnostic tools for
Phytophthora [
24,
37,
38].
We perform a nested real-time PCR approach to increase sensitivity. However, the increase in sensitivity was not so significant likely because single copies of a target DNA were amplified by a single round of PCR, and the nested PCR only improved the signal strength without increasing sensitivity, as already reported by Schena [
37]. Further lowering the sensitivity limit was deemed impractical; instead, improving DNA extraction methods was more effective.
The presence of interference components is a problem of DNA extractions from all environmental samples. DNA extraction from plant material can be challenging, especially from woody plants like walnut or chestnut, containing high amount of tannins that inhibit DNA polymerase in PCR reactions. Our extraction method and the qPCR assay effectively detect
P. cinnamomi in woody host plants obtaining values around 10-20 pg of pathogen DNA /mg host tissue (
Table 3). Sexual oospores, asexual chlamydospores, intracellular hyphal aggregates and lignituber formations are thought to enable
P. cinnamomi survival for long periods under adverse conditions [
39]. Expression of pathogen amount as DNA amount would encompass all these different forms of the pathogen including mycelium.
Soil DNA extraction is even more challenging, due to variability in soil physical and chemical composition and microbiome profile. Main problems are the occurrence of Taq polymerase inhibitors and the low concentration of the target pathogen DNA. We improved DNA extraction by analyzing larger soil samples (10-15 grams). Actually, a further improvement in soil DNA extraction could be obtained by increasing the amount of starting soil material even further, like 20-100 gram of soil. The problem of Taq polymerase inhibitors that were not eliminated with the DNA extraction procedure has been minimized by using Taq polymerase specifically resistant to inhibitors in the qPCR assay. Despite these improvements, still a reduction in sensitivity of the method can occur, especially in agricultural exploited field samples, characterized by low biomass (low DNA) and high inhibitors (
Figure 4). The solution has been to normalize the sensitivity of each specific soil by building standard curve in presence of soil extracts.
Moreover, standardizing sampling protocols for soil sample collection is essential for consistent results. We studied the pathogen distribution around infected walnut trees in Northern Italy orchards, to check if there were typical patterns of
P. cinnamomi presence associated with infected trees and to optimize sampling procedures. Results show that infected trees were consistently associated with
P.
cinnamomi DNA while no pathogen DNA was detected around healthy plants. Infected trees were always associated with
P. cinnamomi at 50 cm from the tree both at 20 cm and 40 cm in depth. At furthest distances from the infected tree, the presence of
P. cinnamomi seem more casual, changing at 15 days apart or in spring/autumn. Phytophthora species seem to stay quite above ground, for example
P. infestans did not percolate through the soil but instead remained at the surface [
40].
Quantifying pathogens in open field soils is challenging due to soil complexity and environmental factors. Precisely because of this complexity, simulation in controlled environment with artificial inoculation in pots couldn't have reflected situation in open field (already just analysis of P. cinnamomi DNA in pot with Rhododendron show us that in pot the pathogen remains more confined and concentrate and detection is easier). Thus, we worked directly in field, necessarily with natural infections.
Results show that
P. cinnamomi infective soil could contain different amount of the pathogen ranging from the qPCR detection limit, set at 0,25 pg of pathogen DNA/g of soil, to 267 pg of pathogen DNA/g of soil (
Table 5). Walnut Tas 11/18 is an exception with so high value of
P. cinnamomi DNA amount (
Table 5), with much lower (<10pg) amount of pathogen DNA/g of soil being the rule. With nested qPCR we show that even lower amount of
P. cinnamomi in soil, <25 pg of pathogen DNA/g of soil (
Table 5), is associated with
P. cinnamomi infected tree.
Correlation between pathogen populations amount and development of symptoms as number of lesions on the plant has already been showed for a Phytophthora species [
40]. More studies would be needed to establish the correlation between amount of
P. cinnamomi and rate of disease on walnut, even if in case of tree this correlation would be difficult to study since the disease develop within years. Anyhow the availability of a quantification method represents the starting point for such studies. Moreover, the use of a quantification method may help the management of samples with DNA values straddling the detection limit, very common for soil samples.
This work shows that the presence of P. cinnamomi DNA is associated with walnut infected trees, and this already helps growers to make informed walnut management decisions.
In conclusion, we developed a qPCR method to detect and quantify P. cinnamomi DNA in mycelium, plant tissue, and especially in field soil samples. Pathogen quantification may contribute to estimate disease potential risk and to set up adequate control strategies avoiding pathogen dissemination This molecular approach is a valuable tool for managing P. cinnamomi in agricultural commercial activities including walnut production. In addition, this method allows the detection of P. cinnamomi in soil prior to plantation/cultivation to prevent future damage.