Intra-Laboratory Evaluation of DNA Extraction Methods and Assessment of a Droplet Digital PCR for the Detection of Xanthomonas citri pv. citri on Different Citrus Species

Xanthomonas citri pv. citri (Xcc) and X. citri pv. aurantifolii (Xca), causal agents of citrus bacterial canker, are both regulated by the European Union to prevent their introduction. Xcc is responsible for severe outbreaks of citrus production worldwide, therefore, a prompt and reliable detection is advisable for the early detection of this bacterium either in symptomatic or asymptomatic plant material. The current EPPO (European and Mediterranean Plant Protection Organization) diagnostic protocol, PM 7/44(1), includes several diagnostic tests even if new assays have been developed in the latter years for which validation data are needed. Recently, a test performance study was organized within the Valitest EU Project to validate Xcc diagnostic methods and provide evidence on the most reliable assays; however, the influence of DNA extraction methods (DEM) on the reliability of the detection has never been assessed. In this study we evaluate four different DEM, by following two different approaches: (i) a comparison by real-time PCR standard curves of bacterial DNA versus bacterial DNA added to plant DNA (lemon, leaves and fruit; orange fruit); and (ii) the evaluation of performance criteria of spiked samples (plant extract added with ten-fold diluted bacterial suspensions at known concentrations). Droplet digital PCR is developed and compared with real-time PCR, as the detection method.


Introduction
Citrus bacterial canker (CBC), caused by Xanthomonas citri pv. citri (Xcc) and X. citri pv. aurantifolii (Xca), is one of the most serious diseases of several host species of the Rutaceae family, among which is included Citrus, Fortunella, Poncirus and their hybrids [1]. Xcc and Xca are the causal agents of the "Asiatic citrus canker" and the "South American citrus canker" (or sometimes false canker), respectively, and the symptoms they cause are morphologically indistinguishable [2]. CBC is mostly a leaf-and fruit-spotting disease, but when conditions are highly favorable, infections cause defoliation, shoot dieback and fruit drop. The two pathovars are distinguished in pathotypes, Xcc in pathotypes A, A w , and A* and Xca in pathotypes B and C. Pathotype A has a wide host range of Citrus and other related genera and a worldwide distribution causing severe outbreaks. Pathotypes A w and A* have a narrow host range and are limited, respectively, in Florida (USA) and Asia. Pathotypes B and C both infect C. aurantifolia, have a restricted host-range and are reported in South America (Argentina, Brazil, Paraguay and Uruguay), and no significant outbreaks have been reported in recent years, suggesting that it has been outcompeted [3].

Samples Homogeneity
The evaluation of homogeneity was performed for samples SET1 and SET2 (see material methods for details). The lower DNA (1 fg/reaction, SET1) and bacterial (10 2-10 cfu/mL, SET2) concentrations showed a high SD (Standard Deviation) (>1), consequently, these samples were not further considered for the standard curve evaluation (SET1) or for the evaluation of performance criteria (SET2), respectively. In particular, all samples at 1 fg/reaction and 10 cfu/mL, and some samples at 10 2 cfu/mL, showed a high incidence of poorly homogeneous samples and/or Ct values very close to those recorded for the previous dilution (∆Ct was less than three cycles as expected).

Development of Droplet Digital PCR
The optimal annealing/elongation temperature was identified at 58 • C, either with 5 or 8 µL of DNA per reaction. At this temperature, the three following parameters were satisfied: (i) positive droplets with the highest fluorescence amplitude; (ii) a better separation between positive/negative droplets; and (iii) less rain (i.e., droplets ranging between the positive and negative ones) (data not shown). The optimal selected volume was 8 µL, as the increase in DNA volume resulted in a higher number of positive droplets with respect to 5 µL. The detection was successful testing the SET1 and SET2 samples, for all DEM/plant matrices combinations (Tables 1 and 2) and samples of bacterial DNA from strain Xcc NCPPB 3234 (Table 3). Figure 1 shows the results obtained on samples of SET2 from 10 3 to 10 cfu/mL. The optimal annealing/elongation temperature was identified at 58 °C, either with 5 or 8 µL of DNA per reaction. At this temperature, the three following parameters were satisfied: (i) positive droplets with the highest fluorescence amplitude; (ii) a better separation between positive/ negative droplets; and (iii) less rain (i.e., droplets ranging between the positive and negative ones) (data not shown). The optimal selected volume was 8 µL, as the increase in DNA volume resulted in a higher number of positive droplets with respect to 5 µL. The detection was successful testing the SET1 and SET2 samples, for all DEM/ plant matrices combinations (Tables 1 and 2). Table 3 shows the results of Figure 1 shows the results obtained on samples of SET2 from 10 3 to 10 cfu/mL.   . For real-time PCR are reported the average of Ct values and the number of positive wells on the total assessed; for ddPCR are reported copies/mL and positive droplets. Grey boxes indicate inconsistent results between the replicates (∆Ct less than 3 cycles between ten-decimal dilutions) and SD > 1. The efficiency, r 2 , and the slope of each real-time PCR are shown. The number of total droplets in ddPCR was >10.000 for all data reported; "-" = not tested samples. DEM = DNA extraction method.  Table 2. Results obtained by real-time PCR and ddPCR in spiked samples of lemon fruit and leaves and orange fruit of SET2 extracted with the four different DEM (Plant = DNeasy Plant Mini kit, Mericon = DNeasy Mericon Food Kit, CTAB = CTAB-based method, Quick = QuickPick SML Plant DNA kit). For real-time PCR are reported the Ct values mean ± SD and the number of positive replicates/number of replicates analyzed; for ddPCR are reported the positive droplets mean, the copies/µL mean ± SD and the number of positive replicates/number of replicates analyzed. Grey boxes indicate inconsistent results between the replicates (∆Ct less than 3 cycles between ten-decimal dilutions) and SD > 1. The number of total droplets in dd-PCR was >10.000 for all data considered for the analyses. Not tested samples are indicate as follows (-). DEM = DNA extraction method.   Table 3. Results obtained by real-time PCR and droplet digital PCR in samples of ten-fold dilution of bacterial DNA from strain Xcc NCPPB 3234, pathotype A were prepared as a standard reference curve. For real-time PCR are reported the average of Ct values and the number of positive wells on the total assessed; for ddPCR are reported copies/mL and positive droplets. Grey boxes indicate inconsistent results between the replicates (∆Ct less than 3 cycles between ten-decimal dilutions) and SD > 1. The efficiency, r 2 , and the slope of each real-time PCR are shown. The number of total droplets in ddPCR was >10.000 for all data reported; "-" = not tested samples.

Evaluation of SET1 Samples
The raw data of real-time PCR and ddPCR results for SET1 are reported in Table 1, according to the DEM and the different plant matrices. The Table 3 shows the result of the standard curve obtained by real-time PCR on bacterial DNA from a pure culture. The amplification efficiency of DNA from the Xcc pure culture showed an adequate value (E = 98.17%); this was comparable to the efficiency shown by most of the combinations of DEM/plant matrices of SET1, which give acceptable values between 90-105%. The exceptions were orange fruit/Mericon (E = 88%), orange fruit/Quick (E = 89%), and lemon leaf/Ctab (E = 111%), whose E values slightly deviate from the acceptable values.
The comparison between the standard curves of bacterial DNA and SET1 (Figure 2), highlight an overlapping for all plant species/matrices, that was particularly evident at the higher bacterial concentration. On the other hand, at lower bacterial concentration some differences were observed due to a higher difference among the Ct (in particular at 10 fg/real-time PCR reaction) ( Figure 2). This was predictable at low target concentrations due to the greater variability between technical replicates.
The range of Ct values for SET1 samples and for bacterial DNA are graphically represented in Figure 3. The range of Ct were not particularly variable among the different combinations of DEM and plant matrices, showing comparable results with bacterial DNA.
The Bland-Altman plot reported in Figure 4 shows the percentage variation of the average of sample concentrations of SET1 (represented by the black line of Figure 4) with respect to a null variation (represented by the red line in Figure 4). The percentage variation indicates whether the method underestimates or overestimates the overall Ct values. A good performance of the method implies that the black line approaches the red one, indicating a low variability. All DEM/plant matrices combinations were within the 95% of confidence interval included in the graph by black dashed lines. In particular, the Plant DEM gave the best performance with all plant matrices, except for the orange fruit.    In order to assess the effect of inhibitors related to the different DEM/plant matrices combination, the SET1 was analyzed by ddPCR, a technique known to be less affected by inhibitors with respect to real-time PCR [11]. The results showed that at the lowest DNA target concentration (10 fg/ddPCR), the orange fruit matrix resulted in the highest number of negative results; however, among the DEM, Plant showed reliable results with all plant matrices as shown in Table 1. Quick showed accurate results with the exception for orange fruit, Ctab was reliable only in samples of lemon leaf and Mericon in none (Table 1).
with respect to a null variation (represented by the red line in Figure 4). The percentage variation indicates whether the method underestimates or overestimates the overall Ct values. A good performance of the method implies that the black line approaches the red one, indicating a low variability. All DEM/plant matrices combinations were within the 95% of confidence interval included in the graph by black dashed lines. In particular, the Plant DEM gave the best performance with all plant matrices, except for the orange fruit.

Evaluation of Healthy Samples and Performance Criteria of SET2
The analysis by real-time PCR [10] of 20 healthy randomly selected samples, sporadically yielded high Ct values and/or inconsistencies between the two technical replicates (i.e., high Ct value/NA) ( Table S1). As reported in Table 4 high Ct/inconsistencies were in different percentages depending on the plant matrices or DEM. In particular, the higher percentages occurred by using the Ctab method (43%), with respect to the Mericon and Plant (respectively, 28% and 25%), with lower values for the Quick (11 %). The healthy status of these 20 samples was assessed by droplet digital PCR, giving a negative result (data not shown). The results of the real-time PCR and ddPCR on SET2 samples are reported in Table 2. The homogeneity showed an acceptable range of variation within bacterial concentrations from 10 7 to 10 3 cfu/mL of SET2 samples ( Table 2). The results of performance criteria ( Table 5) showed 100% of diagnostic specificity, diagnostic sensitivity, and accuracy for all DEM/ plant matrices combinations.
The analytical sensitivity (ASE) was evaluated on SET2 considering different DEM/ plant matrices. In particular, ASE was estimated for real-time PCR [10] and the ddPCR in a range of 10 2 -10 3 cfu/mL and 10-10 2 cfu/mL, respectively, depending on the method used for the DNA extraction (Table 6). By real-time PCR, Quick and Plant showed better results, than Ctab and Mericon. The ASE of 10 2 cfu/mL was obtained by Quick and Plant except for one plant matrix (lemon fruit and orange fruit, respectively) which gave 10 3 cfu/mL. Using Ctab, the ASE was 10 2 cfu/mL only for one out of three plant matrices (lemon leaf) and 10 3 cfu/mL for the other matrices, whereas for Mericon the ASE was at 10 3 cfu/mL for all plant matrices. By ddPCR, using Plant, Mericon and Ctab, the ASE was 10 cfu/mL for two out of three matrices (10 2 cfu/mL was obtained for lemon fruit and orange fruit depending on the DEM) whereas by Quick the ASE was 10 cfu/mL for all plant matrices (Table 6). Table 5. Values of performance criteria obtained by real-time PCR diagnostic sensitivity, diagnostic specificity and accuracy. Positive agreement (PA), negative agreement (NA), positive deviation (PD) and negative deviation (ND) obtained by testing spiked samples contaminated from 10 7 -10 3 cfu/mL of Xcc SET2 and non-spiked healthy samples.

Discussion
The diagnosis of plant pathogenic bacteria represents one of the more important challenges aimed to prevent the dissemination of world-wide diseases. The adoption of a reliable detection method is crucial to correctly determine the presence of a pest before its establishment in a pest-free area and/or its spread. To standardize the use of diagnostics, in plant pathology it is suggested that their validation follow UNI EN ISO 16140 [12], ISO/IEC 17025 [13] and EPPO protocols that allow to evaluate the performance of tests and/or laboratories [14,15]. Intra-or inter-laboratory studies permit the evaluation of several parameters for each diagnostic test in terms of sensitivity, specificity, repeatability, and reproducibility; however, a crucial step is the preparation of DNA extracts that could interfere with the ability of a test to detect the target, regardless of its performance (e.g., sensitivity, specificity).
Furthermore, the EPPO protocols for diagnostics generally do not require DNA quality control as a mandatory step, therefore the DNA is not subjected to control before performing a test. Neither of the majority of diagnostic methods for phytopathogenic bacteria foresee the use of internal DNA amplification controls in each sample (i.e., the mitochondrial cytochrome oxidase, cox gene; the region coding for the 5.8S rRNA, 5.8S rDNA gene), rather, additional samples (PIC, positive isolation) are provided to guarantee the success of the DNA extraction. Moreover, many commercial kits based on standardized procedures have been developed which can ensure the successful outcome of the analysis, and today many kits based on different principles or systems (i.e., columns, beads, manual, automatized) are available. However, the choice of the most appropriate extraction kit is related to various factors ranging from technical aspects (e.g., type of plant matrix) to the economic availability of the laboratory. Taking this last aspect into account, still today many designated laboratories that carry out official analyses use home-made DEM based on the use of Ctab, as for the detection of Xylella fastidiosa [16], one of the most alerted phytopathogenic bacteria in Europe and worldwide.
Xcc, along with Xca, are listed by the EU as a priority pest. Both pathogens were included in the EURL-BAC work program in 2019-2020 with the aim of standardizing, inventorying test protocols, developing and validating detection tests, and providing reference material to standardize detection/identification and facilitate disclosure to national reference laboratories (NRLs). In the case of Xcc, a TPS validation was recently organized within the Valitest Project (XCC-1) (GA n • 773139), aimed at validating several diagnostic methods for Xcc detection. The obtained results provided useful indications on the validated diagnostic methods, highlighting their performance. An adequate availability of diagnostic tools for the detection of Xcc emerged from the report by the Valitest Project (XCC-1) [17], highlighting that the real-time PCR techniques are more sensitive than respective conventional PCR, except for the real-time PCR of Mavrodieva et al. [18]. This validation activity did not consider the DNA extraction and the influence that this step may have for the reliability of the detection both in conditions of high and low levels of infection of the plant sample, or in relation to different plant matrices. To avoid the occurrence of false negative results in the molecular detection test, the DNA extraction step must ensure the effective removal of DNA polymerase inhibitors [8,9]. In this study, the possible influence of several combinations of DEM/plant matrices were assessed with the real-time PCR detection test of Cubero and Graham [11]. The real-time PCR developed by Ròbene et al. [19] was not taken into consideration as the primers/probe had not yet been published at the time of the TPS execution. The Valitest TPS has been considered reliable for the detection of Xcc for both tests [17]. For this study, Xcc pathotype A was selected, considering its broader plant host-range and wider spread in the infected areas, compared to the A W and A* pathotypes [5]. The overall results showed that at high-medium bacterial (10 7 -10 3 cfu/mL) or DNA (1 ng-100 fg/PCR reaction) concentrations all DEM allow a reliable detection of Xcc from the tested plant matrices (lemon and orange fruits and lemon leaves). This evidence was highlighted: (i) by the overlapping of the standard curves of bacterial DNA and the samples of SET1; (ii) by the acceptable efficiency values for most samples of SET1 showing a scarce influence of the DNA quality on the amplification efficiency; and (iii) by the value of 100% of diagnostic sensitivity, diagnostic specificity and of the accuracy for samples of SET2.
At the lowest concentrations (10 fg/PCR reaction-SET1 and 10 cfu/mL-SET2), most of the samples did not show an acceptable homogeneity or a reliable result by real-time PCR, which resulted in variable results, frequently non-consistent, depending on the plant matrices /methods combination. This evidence shows that at this condition the DNA extraction methodology can affect the reliability of the detection tests. The evaluation of the ASE of SET2 (spiked samples) by the real-time PCR of Cubero and Graham [10] showed a range of 10 2 -10 3 cfu/mL, in accordance with the results of the XCC-1 TPS Valitest [17]; however, this value ranged depending on the different DEM/plant matrices combination, showing the best ASE values for Quick and Plant, intermediate for Ctab and the worst for Mericon.
The ddPCR is notoriously less affected by inhibitors [11]. The presence of inhibitors has the potential to increase errors, reducing the method resolution and producing false negative results [20]. Our ddPCR outcome showed negative results with most of the samples of the Xcc SET1 at 1fg/ddPCR due to the limit of detection of the method. At 10fg/ddPCR, negative results were obtained by Mericon with all matrices, by Ctab with orange and lemon fruit, by Quick with orange fruit, whereas Plant detected the target by all matrices. By assessing SET2, the ddPCR showed a ten-fold higher ASE with respect to real-time PCR, with the best ASE values obtained when Quick was used for the DEM.
The overall results indicate that the use of commercial kits is adequate for a reliable DNA extraction for all the tested plant matrices at a high concentration of the target. Particular attention must be paid to the choice of the most appropriate extraction method when processing samples with an expected low bacterial load, as in the case of asymptomatic samples. In this case, the use of Quick and Plant seems the most reliable DEM. An adequate choice of the DNA extraction procedure allows to improve the sensitivity of the method by one decimal factor. In this regard, the adaptation of a ddPCR system according to the Cubero and Graham test [10] has shown that the latter is more sensitive than the respective real-time PCR method. Dupas et al. (2019) [21] reported that the sensitivity of ddPCR compared to real-time PCR is controversial, because some authors have reported an increasing in sensitivity but others showed a 10-100-fold lower sensitivity. The benefits of ddPCR seem to be related to the studied pathosystem [22], e.g., the detection of Ralstonia solanacearum was improved using ddPCR compared to real-time PCR but was comparable for Erwinia amylovora. The ddPCR adapted from Harper et al. (2010) [23] for the detection of Xylella fastidiosa showed the same sensitivity of the corresponding real-time PCR for several plant species (Olea europaea, Polygala myrtifolia and Rosmarinus officinalis) and an improved sensitivity for others (Quercus ilex and Lavandula angustifolia) [21], that are known to be rich in plant inhibitors. An improved sensitivity of the ddPCR may therefore depend both on the higher tolerance of the ddPCR to inhibitors and on the lower coefficient of variation of the ddPCR compared to real-time PCR, especially at low target concentrations [12]. Zhao et al. [11] developed a ddPCR for the detection of Xcc reporting its potential for quantitative detection with a high precision and accuracy with respect to a real-time PCR assay, and its suitability for the diagnosis of pathogens in field samples with a complex matrix of inhibitors.
Finally, in our own laboratory conditions late Ct values occurred by real-time PCR on healthy plant matrices, with a high incidence for Ctab and low for Quick, supporting the argument that Quickis the most suitable method for DNA extraction. A ddPCR was performed to evaluate if the late Ct were due to background contaminations and the obtained negative results confirmed that these late Ct values were false signals. We suggest establishing a cut-off for real-time PCR [10] to avoid the occurrence of false-positive, due to late Ct values, that interfere with the detection of the bacterium. This aspect must be considered particularly for the analysis of asymptomatic material, where the choice of an adequate DEM could be crucial for the success of the diagnosis. In conclusion, this intra-laboratory study highlighted that all DEM can be used for the analysis of symptomatic material; however, the choice of Quickand Plant seem the most appropriate for the successful detection of Xcc in samples with a low concentration of the pathogen, such as in asymptomatic material. It is worth noting that the surveillance protocol should consider possible asymptomatic but infected material (EFSA) [6]. Indeed, it must be taken into account that the major risk factors for the long-distance spread capacity of Xcc are represented by the commercial shipment of diseased or contaminated fruit [24] or plant propagating material (e.g., budwood, rootstock, seedlings) [2,25].
Although molecular diagnostic tests allow for reaching the levels of sensitivity adequate for the detection of Xcc in low concentrations, the validation data produced in this study show that the wrong choice of the extraction method can affect the outcome of official analyses, especially in critical samples (i.e., low concentration of the pathogen, pest-free areas). Conversely, applying reliable diagnostics will help protect pest-free areas from accidental introduction of this quarantine pathogen.

Bacterial Strain
The strain NCPPB 3234 (Xcc) was used for the samples' preparation. This strain, belonging to the National Collection of Plant Pathogenic Bacteria (NCPPB Fera Science Ltd -UK), was isolated in 1982 in Japan from Citrus spp. and belong to pathotype A. The lyophilized strain was grown in NGA (nutrient agar added with 0.25% d-glucose) at 28 • C. Bacterial suspensions were prepared in phosphate saline buffer (PBS 10 mM, pH = 7.2) and spectrophotometrically (DS-11 Fx+, Spectrophotometer-Fluorometer Denovix Inc., Wilmington, DE, USA) measured at a concentration corresponding to about 10 8 colony forming units (CFU) mL −1 (OD 660 = 0.1). The ten-fold dilution of bacterial DNA was prepared as a standard reference curve.

Plant Material
Plant material of the host plants was collected from plants in open fields (leaves) or bought in trade markets (fruits) located in the Latium and Campania Regions (Italy) ( Table 7). The plant extract was prepared as reported in [6] with some modifications as indicated below. Thirty grams of leaf tissue and/or fruit surface layer was collected for the orange and lemon (Table 7) and cut into smaller pieces. Each sample was subdivided into 3 subsamples of 10 g, washed with PBS in sterile bags (Bioreba, Reinach, Switzerland), and shaken for 20 min at room temperature. After incubation, the supernatant was filtered to remove the plant material and centrifuged for 20 min at 10,000× g. The pellet of each subsample was resuspended in 10 mL of PBS. Finally, the three subsamples were combined into one sample to proceed with the steps reported below. The healthy status of the plant extract was assessed for each plant matrix by real-time PCR [10] before proceeding with samples preparation and used only in the case of negative results. The plant extracts were used for the plant DNA extraction performed with four different DEM as reported below.

Sample Preparation
Samples prepared for the intra-laboratory study were grouped in SET1 and SET2. SET1 was prepared by amending the DNA of the plant extract (orange fruit, lemon leaves and lemon fruit) obtained with the four different DEM, with a ten-fold dilution of Xcc DNA (from 1 ng to 1 fg/real-time PCR or ddPCR reaction). SET2 consisted of spiked samples prepared by adding bacterial suspensions at known concentrations (CFU/mL) to healthy plant extracts (orange fruit, lemon leaves and lemon fruit). The obtained spiked samples (bacterial suspension + plant extract) were extracted by DEM. The bacterial suspension was added in the different plant matrices at a final concentration of Xcc 10 7 , 10 6 , 10 5 , 10 4 , 10 3 , 10 2 , and 10 cfu/mL. Moreover, a plant extract without the addition of bacteria was prepared (NIC). Three independent biological replicates (for each plant matrix, extracted with the different DEM) were prepared for each bacterial concentration. The number of colonies forming units (CFUs) was determined by plating 100 µL of bacterial suspensions on an NGA medium and incubating at about 27 • C and determined after 4 days. Twenty samples for each plant matrix randomly selected were extracted with the four DEM and evaluated by real-time PCR and droplet digital PCR as non-target samples for the evaluation of the performance criteria.

DNA Extraction
Total bacterial DNA was extracted from 1 mL of bacterial cultures of Xcc NCPPB 3234 using Gentra Puregene Yeast/Bact. Kit (Qiagen, Venlo, The Netherlands). The DNA concentration was evaluated by Qubit (dsDNA HS Assay kit, Invitrogen, Waltham, MA, USA). Plant DNA extraction was performed following the manufacturers' instructions of the commercial kits, i.e., DNeasy Plant Mini kit (Qiagen, Venlo, The Netherlands), DNeasy Mericon TM Food Kit (Qiagen, Venlo, The Netherlands), and QuickPick™ SML Plant DNA kit (QRET Technologies Ltd., Turku, Finland), manual version, from 500 mL of plant extract. The Ctab DEM was performed as described for Xylella fastidiosa [16]. The DNA was stored ≤−15 • C until analysis.

Real-Time (Real-Time PCR)
The Cubero and Graham [11] real-time PCR reaction was performed according to the TPS instructions for the detection of Xcc (TPS code: XCC1) [17] received in the frame of the EU project Valitest (GA n • 773139). TaqMan reactions were carried out in a 25 µL reaction mixture containing 0.4 mM of each primer (J-pth3/4), (Eurofins Genomics Germany GmbH, Ebersberg, Germany) 1× GoTaq®Probe qPCR Master Mix (Promega Corporation, Madison, WI, USA), 0.2 µM TaqMan probe J-Taqpth2 (Eurofins), 5 µL of DNA extracts; the thermal profile consisted of 45 cycles at 95 • C-10 min, 95 • C-15 s, and 60 • C 1 min. All real-time PCR runs were performed including at least a control sample of PAC (positive amplification control), NIC (negative internal control) and NAC (negative amplification control). The samples of SET1 were run in one biological replicate, each in three technical replicates (n = 3). The samples of SET 2 were run in three biological replicates (for each plant matrix), each in technical triplicate (n = 9).

Droplet Digital PCR (ddPCR)
The ddPCR was developed adapting the method of Cubero and Graham [10]. For this purpose, the optimal annealing/extension temperature was identified by observing three different parameters: (i) the highest fluorescence amplitude; (ii) the positive/negative droplets disengagement; and (iii) less rain (i.e., droplets ranging between the positive and negative ones).
Following Biorad ddPCR TM Supermix datasheet information, two thermal gradients were tested to determine the optimal hybridization temperature ranging from 53 • C to 60 • C on 5 µL or 8 µL of DNA Xcc, strain NCPPB 3234. The optimized ddPCR reaction mix retained consisted of 1xddPCRTM Supermix for Probes no dUTP ( Bio-Rad Laboratories Inc., Hercules, CA, USA), 900 nM Jpht3, Jpth4 primers, 250 nM P-Jpth2 probe and 8 µL of the DNA sample. The optimal thermocycler conditions retained consisted of 95 • C 10 min, then 40 cycles of two-steps of 94 • C for 30 s and 58 • C for 1 min, followed by 98 • C for 10 min. A temperature ramp of 2.5 • C/s was fixed in all PCR steps. The ddPCR reactions with fewer than 10,000 droplets generated were excluded from the analyses; a result was considered positive if at least two positive droplets were individuated.
The samples of SET1 were run in one biological replicate. The samples of SET 2 were run in three biological replicates (for each plant matrix) (n = 3). In particular, there were tested samples from 100 pg to 10 fg/ddPCR for SET1 and samples from 10 3 cfu/mL to 10 cfu/mL for SET2, respectively. Ten healthy samples from each tested plant species previously assessed by real-time PCR were also tested by ddPCR, QX 200 Droplet Reader, QX 200 Droplet Generator, PX1 PCR Plate Sealer, QuantaSoft 1.7.4.0917 (Biorad, Bio-Rad Laboratories Inc., Hercules, CA, USA).

Samples Homogeneity and Performance Criteria Evaluation
The homogeneity was evaluated for samples of SET1 and SET2 by considering the (SD) of Ct values obtained for each sample by real-time PCR. The replicates were considered sufficiently homogeneous with SD < 1.
The following performance criteria, diagnostic specificity, diagnostic sensitivity, accuracy, repeatability, and analytical sensitivity were evaluated on samples of SET2 according to the EPPO PM 7/98 (4) [13]. In particular, the following target samples were assessed: 3 biological replicates (in 3 repetitions, n = 9) from 10 7 to 10 3 cfu/mL excluding the samples that did not respect the homogeneity parameters (SD > 1). A total of 45 target samples were assessed. As the non-targets, 20 healthy samples of each of the plant matrices were tested. Some interference was observed on healthy matrices by real-time PCR (high Ct value in duplicates or inconsistencies between the two replicates (i.e., Ct/NA value) whose results are reported in Table S1.

Statistical Methods
The sample mean ± the standard deviation (SD) was calculated for all data of the real-time PCR and ddPCR. Excel Microsoft 365®was employed as the statistic package. The analysis of qualitative data was performed by the calculation of performance criteria as reported in the PM 7/122 (1). The Bland-Altman test was performed through R version 4.0.2, using the packages "base" and "stats" [26]. Funding: This research was funded by the European Union Commission (Grant agreement SI2.809508) and was performed in the frame of "EURL-BAC-European Union Reference Laboratory for Bacteriology" activity (2019-2020). Part of these data were uploaded to the EPPO database of the European Plant Protection Organization (EPPO) on diagnostic experience as validation data.
Institutional Review Board Statement: Not applicable in this study.