# Metagenomic Sequencing for Identification of Xylella fastidiosa from Leaf Samples

Verónica Román-Reyna, Enora Dupas, Sophie Cesbron, Guido Marchi, Sara Campigli, Mary Ann Hansen, Elizabeth Bush, Melanie Prarat, Katherine Shiplett, Melanie Ivey, et al.

► To cite this version:

Verónica Román-Reyna, Enora Dupas, Sophie Cesbron, Guido Marchi, Sara Campigli, et al. Metagenomic Sequencing for Identification of Xylella fastidiosa from Leaf Samples. mSystems, 2021, 6 (5), pp.e0059121. 10.1128/mSystems.00591-21. hal-03541606

HAL Id: hal-03541606

https://hal.inrae.fr/hal-03541606

Submitted on 16 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Metagenomic Sequencing for Identification of *Xylella fastidiosa* from Leaf Samples

Verónica Román-Reyna,ª,ᵇ Enora Dupas,ᶜ,ᵈ Sophie Cesbron,ᶜ Guido Marchi,ᵉ Sara Campigli,ᵉ Mary Ann Hansen,ᵉ Elizabeth Bush,ᵉ Melanie Prarat,ᵉ Katherine Shiplett,ᵉ Melanie L. Lewis Ivey,ᵇ Joy Pierzynski,<fnref n="1" /> Sally A. Miller,ᵇ,ᵉ Francesca Peduto Hand,ª Marie-Agnès Jacques,ᶜ Jonathan M. Jacobsª,ᵇ

**ABSTRACT** *Xylella fastidiosa* (*Xf*) is a globally distributed plant-pathogenic bacterium. The primary control strategy for *Xf* diseases is eradicating infected plants; therefore, timely and accurate detection is necessary to prevent crop losses and further pathogen dispersal. Conventional *Xf* diagnostics primarily relies on quantitative PCR (qPCR) assays. However, these methods do not consider new or emerging variants due to pathogen genetic recombination and sensitivity limitations. We developed and tested a metagenomics pipeline using in-house short-read sequencing as a complementary approach for affordable, fast, and highly accurate *Xf* detection. We used metagenomics to identify *Xf* to the strain level in single- and mixed-infected plant samples at concentrations as low as 1 pg of bacterial DNA per gram of tissue. We also tested naturally infected samples from various plant species originating from Europe and the United States. We identified *Xf* subspecies in samples previously considered inconclusive with real-time PCR (quantification cycle [Cq], >35). Overall, we showed the versatility of the pipeline by using different plant hosts and DNA extraction methods. Our pipeline provides taxonomic and functional information for *Xf* diagnostics without extensive knowledge of the disease. This pipeline demonstrates that metagenomics can be used for early detection of *Xf* and incorporated as a tool to inform disease management strategies.

**IMPORTANCE** Destructive *Xylella fastidiosa* (*Xf*) outbreaks in Europe highlight this pathogen's capacity to expand its host range and geographical distribution. The current disease diagnostic approaches are limited by a multiple-step process, biases to known sequences, and detection limits. We developed a low-cost, user-friendly metagenomic sequencing tool for *Xf* detection. In less than 3 days, we were able to identify *Xf* subspecies and strains in field-collected samples. Overall, our pipeline is a diagnostics tool that could be easily extended to other plant-pathogen interactions and implemented for emerging plant threat surveillance.

**KEYWORDS** *Xylella fastidiosa*, metagenomics, diagnostics, short-read sequencing

*Xylella fastidiosa* (*Xf*) is a globally distributed insect-transmitted plant-pathogenic bacterium causing diseases on a large host range. To date, 595 plant species belonging to 85 botanical families have been reported as *Xf* hosts (1), some of which are of major socioeconomic interest, such as grapevine, olive, citrus, coffee, and almond (2). *Xf* colonizes the

**Citation** Román-Reyna V, Dupas E, Cesbron S, Marchi G, Campigli S, Hansen MA, Bush E, Prarat N, Shiplett K, Ivey MLL, Pierzynski J, Miller SA, Peduto Hand F, Jacques M-A, Jacobs JM. 2021. Metagenomic sequencing for identification of *Xylella fastidiosa* from leaf samples. *mSystems* 6:e00591-21. https://doi.org/10.1128/mSystems.00591-21.

**Editor** Davide Bulgarelli, University of Dundee

**Ad Hoc Peer Reviewer** Alessio Aprile, Università del Salento

**Copyright** © 2021 Román-Reyna et al. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Address correspondence to Jonathan M. Jacobs, jacobs.1080@osu.edu.

**Received** 12 May 2021

**Accepted** 8 September 2021

**Published** 26 October 2021

<footnote>ªDepartment of Plant Pathology, The Ohio State University, Columbus, Ohio, USA</footnote>

<footnote>ᵇInfectious Disease Institute, The Ohio State University, Columbus, Ohio, USA</footnote>

<footnote>ᶜUniversity of Angers, Institut Agro, INRAE, IRHS, SFR QUASAV, Angers, France</footnote>

<footnote>ᵈFrench Agency for Food, Environmental, and Occupational Health & Safety, Plant Health Laboratory, Angers, France</footnote>

<footnote>ᵉDepartment of Agriculture, Food, Environment, and Forestry, University of Florence, Florence, Italy</footnote>

<footnote>>fSchool of Plant and Environmental Sciences, Virginia Tech, Blacksburg, Virginia, USA</footnote>

<footnote n="9">⁹Animal Disease Diagnostic Laboratory, Ohio Department of Agriculture, Reynoldsburg, Ohio, USA</footnote>

<footnote>ᵇDepartment of Plant Pathology, The Ohio State University, Wooster, Ohio, USA</footnote>

<footnote n="1">¹C. Wayne Ellett Plant and Pest Diagnostic Clinic, Department of Plant Pathology, The Ohio State University, Reynoldsburg, Ohio, USA</footnote>
xylem vessels of plants, where it forms biofilms (3) that, together with tyloses and gums pro-
duced by the plant in response to the infection (4), limit water translocation. Infected hosts
display symptoms of leaf scorches and plant dieback finally followed by plant death (3).

Xf was first described in and limited to the Americas but recently emerged in Europe, highlighting the pathogen's capacity to expand its host range and geographical distribution (2, 5). The pathogen was reported in Italy in 2013, where is currently devastating Apulian olive production, then detected in France in 2015, Spain in 2016 and Portugal in 2018, on both cultivated and spontaneous Mediterranean plant species (2). The primary control strategy for diseases caused by Xf includes eradication of infected hosts via early detection. Therefore, developing methods with the lowest limit of detection are critical to prevent major losses for growers and future pathogen dispersal.

The diagnostics of diseases caused by fastidious pathogens such as Xf are difficult.
This difficulty is increased, as infected plants may remain asymptomatic for very long
periods of time, which is associated with low bacterial concentrations and an irregular
distribution of the pathogens in the plants (6). It is of major interest to develop reliable
and highly sensitive tools for detection and detailed identification that can be used
directly on plant extracts. Current standards for Xf diagnostics primarily rely on quanti-
tative real-time PCR (qPCR) assays to detect and sometimes identify the bacterium (7-12),
followed by the amplification and sequencing of two, for subspecies identification, to seven
housekeeping genes (cysG, gltT, holC, leuA, malF, nuoL, and petC) for sequence type (ST)
determination and phylogeny reconstruction (2) (Fig. 1A). Five subspecies are proposed in X.
fastidiosa, i.e., fastidiosa, multiplex, pauca, morus, and sandyi (13-15). However, whole-ge-
nome analyses revealed similarities of the subspecies fastidiosa, morus, and sandyi, which
cluster into one clade. Moreover, genome analysis indicated a high frequency of horizontal
gene transfer and recombination among Xf subspecies (14-16).

Plant samples infected by more than one Xf strain belonging to several subspecies are not uncommon and are not easy to detect (17, 18). Nevertheless, current methods do not consider new or emerging variants resulting from pathogen genetic recombination (14). For example, qPCR with high quantification cycle (Cq) values (>35) are considered inconclusive (2), making decisions about disease control difficult. A complementary tool for diagnostics is the use of next-generation sequencing (19) (Fig. 1). Because this approach can be directly used on plant extracts, it is not biased toward known sequences and provides more information about the pathogen genome, such as virulence traits. Metagenomics, the study of genetic material from environmental samples, beyond whole-genome sequencing, allows for the detection of strains from several subspecies and ST at the same time from the host (20). Recently, the use of long-read sequencing as diagnostic tool identified Xf subspecies and ST from infected samples (12, 21).

In this study, we developed and tested a metagenomics pipeline using in-house
short-read sequencing as a complementary approach for affordable and accurate Xf
detection. We were able to use metagenomics to identify Xf to the strain level in single- and
mixed-infected plant samples, at concentrations as low as 1 pg of bacterial DNA per gram of
tissue. In addition, we tested naturally infected field samples from Europe and the United
States. We identified Xf subspecies in samples with Cq values equal to and greater than 37,
which is beyond the threshold of detection for the standard and certified qPCR methods (2).
Overall, we developed a robust diagnostics pipeline that could be easily extended to other
pathogens and implemented for surveillance of emerging agricultural threats.

RESULTS

**Metagenomics for the diagnostics pipeline.** We developed and tested a metagenomics pipeline for *Xf* detection and subspecies identification (Fig. 1). We tested this pipeline based on three types of DNA samples—from bacterial colonies in culture, spiked plant samples, and naturally infected plant samples (Fig. 1B). To recover and identify *Xf* subspecies and compare them to the already sequenced genomes, we developed a pipeline that uses six different tools and custom-made databases (22) (Fig. 1C). The pipeline recovers *Xf* reads with the software Kraken 2 and a custom-made database (22). The database has
A

Sequenced-based detection

Conventional detection

<figure><img src="image_3.png" /></figure>

iSeq 100 sequencing

<figure><img src="image_4.png" /></figure>

B

Sample and DNA preparation

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

Sample preparation

C

Analysis for pathogen identification

<figure><img src="image_8.png" /><figcaption>**FIG 1** Metagenomics for the diagnostic pipeline. (A) Sequenced-based detection. Two approaches were used for *Xf* detection—conventional detection and iSeq 100 sequencing. For conventional (Continued on next page)</figcaption></figure>
user-specified genomes for *Xf* read identification. The user-specified genomes belonged to *Xylella* (n = 81), *Xanthomonas* sp. (n = 10), *Escherichia coli* (n = 1), and several plant sequences from NCBI (Table S1). The custom-made database had plant sequences to avoid false-positive results because we found that some NCBI *Xf* genomes contained plant genomic DNA sequences. The plant DNA sequence hits had 100% identity to plant 18S or chloroplast reads. We could not remove all plant reads from the 81 NCBI *Xf* genomes. Therefore, the plant reads in the database serve as a filter to ensure that plant reads were not misidentified as *Xf* reads.

After Kraken 2, the pipeline *de novo* assembled the recovered *Xf* reads into contigs with the program SPAdes (23). The pipeline used the *Xf* contigs for four different analyses: (i) subspecies identification, (ii) phylogeny reconstruction, (iii) identification of the already sequenced genetically closest strains, and (iv) alleles for multilocus sequence type (MLST) profile and virulence-related genes determination. The pipeline used *Xf* contigs and the tool SendSketch from the BBMap software to identify subspecies. Then it used Pyani and LINbase software to reconstruct phylogeny by average nucleotide identity (ANI) (24, 25). Next, it assigned *Xf* strains to each *Xf* contig to identify the closest strains with the tool BBSplit from the BBMap software. Finally, to identify specific genes or alleles, the pipeline used local BLAST search with two types of subjects; (i) one subject was the reported MLST allele genes, and (ii) the other subject was the protein sequences from genes associated with virulence.

**Xf identification from in silico prepared samples.** To test the pipeline sensitivity, we used *in silico* samples with target (e.g., *Xf*) and nontarget DNAs (e.g., nonhost plant). The samples included variable amounts of nonhost barley (*Hordeum vulgare*) sequenced reads *in silico* spiked with *Xylella fastidiosa* subsp. *fastidiosa* (*Xff*) CFBP 7970 reads. We obtained a strong linear correlation between the Kraken 2 results and the proportion of spiked *Xf* sequence reads (percentage of *Xylella* reads spiked in sample, y = 103.21x – 0.0127; R² = 1) (Fig. S1A). We recovered *Xf* reads and assembled them as contigs using SPAdes. With the *Xf* contigs, we performed BLAST analysis to identify MLST alleles and virulence genes. We were able to identify one to four MLST-related genes for samples spiked with 0.5 to 2.4% *Xf* reads (Table S3). This result indicated that we cannot capture the full MLST gene set for ST identification with less than 2.4% *Xf* reads (Table S3). We calculated, for all samples, the percentage of gene similarity to the virulence-related genes (Table S4). The percentage of gene similarity increased with the higher number of spiked *Xf* reads. Samples with a lower number of *Xf* reads had a low genome coverage to recover and analyze complete gene sequences (Tables S4 and S6).

We then identified *Xf* subspecies using the *Xf* contigs. Since the *in silico* samples only identified *Xff* reads, we expected that SendSketch assigned all contigs to *Xff*. However, we found that 9 to 15% of *Xff* contigs were instead assigned to *Xfm*. Based on these results, we did two additional analyses to determine the best approach to analyze the *Xf* subspecies composition. For the first analysis, we hypothesized that complete assembled genomes

**FIG 1** Legend (Continued)

detection, samples were analyzed using qPCR assays, Harper's test or tetraplex Dupas's test, and MLST involving Sanger sequencing of seven housekeeping genes. iSeq 100 libraries were prepared according to the manufacturer's instructions. After 17 h of sequencing, demultiplexed samples were recovered from the machine and used for subsequent analysis. (B) Sample preparation. The samples used for the pipeline were DNA-extracted from bacterial strains in culture, spiked plant material, and naturally infected samples. (C) Pathogen identification via metagenomic analysis. Demultiplexed fastq reads from all samples were then used for metagenomic analysis. We created a database to recover *Xf* reads using Kraken 2. The database contained *Xylella*, *Xanthomonas*, and *Escherichia coli* genomes. We also added plant plastid genomes to remove false-positive results. *Xf* reads were recovered from the fastq files. The *Xf* recovered reads were *de novo* assembled to obtain *Xf* contigs, using SPAdes. The *Xf* contigs were used in four different analyses, subspecies identification, phylogeny reconstruction, identification of the genetically closest strain with a sequenced genome, and alleles from specific genes. To determine subspecies, we used the tool SendSketch. To reconstruct phylogeny, we calculated ANI using Pyani and the website tool LINbase (https://linbase.org/). To determine the genetically closest known *Xf* strain, we detected the number of hits to each *Xf* strain using BBSplit. To identify specific gene alleles, we calculated the percentage of identity to the seven MLST genes (*cysG*, *gltT*, *holC*, *leuA*, *malF*, *nuoL*, *petC*) and the percentage of similarity to 17 virulence-related proteins using local BLAST+. Graphics were created with BioRender.
would reduce the percentage of reads assigned to other subspecies. To test this, we created a smaller Kraken2 database with 30 genomes instead of 81. These 30 *Xf* genomes had complete assemblies. We recovered 1% to 2% fewer *Xf* reads with the new database, and the subspecies distribution remained the same (data not shown). The results indicated that *Xf* subspecies classification is not related to the level of genome assembly. Therefore, the original Kraken2 database with 81 NCBI *Xf* genomes was retained for all further analyses.

For the second analysis, we manually assessed the SendSketch sensitivity to mixed infections with new *in-silico* samples. The new samples included a set amount of barley reads *in silico* spiked with variable amounts of *Xff* CFBP 7970 and *X. fastidiosa* subsp. *multiplex* (*Xfm*) CFBP8418 reads (Table S2). For these new *in silico* samples, we recovered *Xf* reads with Kraken 2, assembled the reads as contigs, and ran SendSketch to identify subspecies. When using BLAST, a certain number of *Xf* contigs mapped equally (100% identity) to *Xff* and *Xfm* (*Xf* core contigs). We observed that *Xf* core contigs are directly proportional to the total of *Xf* recovered reads and samples with a higher *Xff* to *Xfm* spiked read ratio (Table S2, Fig. S1B and C). Moreover, the tool SendSketch randomly assigned the subspecies to *Xf* contigs with 100% identity to *Xff* and *Xfm*. Consequently, we developed a manual correction to separate single from mixed infections. We only used samples with either *Xfm* or *Xff*, *X. fastidiosa* subsp. *pauca* (*Xfp*) was not part of the analysis. The correction consists of calculating the logarithm of the *Xfm:Xff* contig ratio. Corrected log ratios from 0.081 to 0.4 are considered a mixed infection. Log ratios below -0.012 are identified as single *Xff*, and those higher than 0.43 are *Xfm* single infection (Fig. S1B and C).

To evaluate the pipeline with samples free of *Xf*, we used extracted DNAs of two healthy plant samples and a non-*Xylella* control (i.e., barley leaves infiltrated with *Xanthomonas*). For the artificially inoculated barley samples, Kraken 2 software recovered 20 to 30% of the total reads as *Xanthomonas*. For all four *Xf*-free samples, Kraken 2 recovered 6 to 19 *Xf* reads (Table S2). All these *Xf* reads corresponded to plant reads based on the BLAST Web tool from NCBI. Based on these results, the pipeline considers a sample *Xf*-free when it cannot recover more than 19 *Xf* reads (Fig. S2).

**Xf identification from isolated bacteria.** To test the capacity of iSeq 100 sequencing, we used six genomic DNAs (gDNAs) from isolated bacteria and two known *Xf* gDNAs as controls. The six *Xf* gDNAs were isolated from Italian field samples (See Materials and Methods). The two control *Xf* gDNA samples were *Xff* CFBP 7970 (CFBP 7970 iSeq100) and *Xfm* CFBP 8418 (CFBP 8418 iSeq100). All eight gDNA samples were sequenced with the iSeq 100 system. For all eight samples, Kraken 2 recovered 99% of the total reads as *Xf* reads (Table S2). We assembled the *Xf* reads as contigs and classified them into subspecies. For CFBP 7970 and CFBP 8418, for which a genome was already available, 54 to 78% of the contigs corresponded to *Xf* core contigs, 44 to 65% to their respective subspecies, and 2% to the closest subspecies. On average, within the six Italian samples, 20% of the contigs corresponded to *Xf* core contigs, 25% to *Xfm*, and 2% to *Xff*.

The ANI values were consistent with the *Xf* contig abundance (Fig. 2, Table S2). All six Italian samples and CFBP 8418 iSeq100 had 99 to 100% identity to *Xfm* and less than 97% identity to *Xff* and *Xfp*. The control sample, CFBP 7970 iseq100, had 100% identity to *Xff* and less than 98% identity to *Xfm* and *Xfp* (Fig. 2, Table S2).

For strain identification, we used the program BBSplit and the Harvest suite. For each sample, we selected the top three closest strains based on the program BBSplit output. Then, we used these closest strains and the sample to compare the number of single nucleotide polymorphisms (SNPs) with the Harvest suite (26). For the sample CFBP 8418 iSeq100, the closest strain with 30 SNPs was *Xfm* CFBP 8418 (Tables S3 and S5). For the sample CFBP 7970 iSeq100, the closest strain with fewer SNPs was *Xff* CFBP 7970. The six Italian samples had the same three closest *Xfm* strains, TOS5, TOS4, and TOS14. All six samples had fewer SNPs than the strain *Xfm* TOS4. The three TOS strains and the Italian samples were isolated from the outbreak area of Monte Argentario, Tuscany, Italy (27).

We performed BLAST analysis to identify MLST alleles and virulence genes for all eight isolated bacteria with the assembled *Xf* contigs. For virulence genes, the sample CFBP 7970 iSeq100 had 100% similarity to all *Xff* virulence genes except for rpfE (96.4%)
<figure><img src="image_9.png" /><figcaption>**FIG 2** Phylogenetic reconstruction of isolated bacteria used in this study. The cluster analysis is based on average nucleotide identity values from Pyani. Branch colors indicate different *Xf* subspecies—*Xff* (green), *Xfm* (blue), *Xfp* (gray), *Xf* subspecies *morus* and *sandyi* (black). The sequenced gDNA from isolated bacteria are indicated in blue or green. The *Xf* genomes obtained from NCBI are indicated in black. The cluster was plotted using ComplexHeatmap R package.</figcaption></figure>

(Table S4). The sample CFBP 8418 iSeq100 had 100% similarity to all *Xfm* virulence genes except for pilB (99.8%). The six Italian samples had the same similarity percentages for all *Xfm* virulence genes except for hemagglutinin (95.7 to 100%). We were able to identify all virulence genes and to complete the allelic profiles for ST identification (Table S3). As expected, the ST identified for the sample CFBP 7970 iSeq100 was ST2, and the one for CFBP 8418 iSeq100 was ST6. The six Italian samples had the same ST87 number.

**Xf identification from spiked plant samples.** We tested the pipeline with DNA extracted from grapevine petioles and midribs artificially inoculated with known bacterial concentrations of the strain *Xff* CFBP 7970, *Xfm* CFBP 8418, or an equal mix of both strains. Kraken 2 output recovered 0.01 to 74.8% of the total sequences as *Xf* reads (Table S2). The percentage of recovered *Xf* reads had a positive correlation with log₁₀ CFU values (R² = 0.9876) (Fig. 3A) and C_q values (R² = 0.9141) (data not shown). The pipeline detected *Xf* with the lowest bacterial concentration tested in this study (1 × 10⁴ CFU/ml), equivalent to a C_q value of 28.85 and 1.62 pg·µl⁻¹.
A

<figure><img src="image_10.png" /></figure>

B

Percentage of Xf recovered reads

Kraken 2 to recover Xf reads

<figure><img src="image_11.png" /></figure>

C

SendSketch to determine Xf subspecies distribution

<figure><img src="image_12.png" /></figure>

**FIG 3** Spiked samples with different dilutions and mixed samples. (A) Scatterplot comparing the log₁₀ CFU with the percentage of Xff CFBP7970 recovered reads from the total read number. The dotted

(Continued on next page)
After *Xf* contig assembly, we were able to identify the subspecies for all samples, and the log ratio separated single from mixed infections (Fig. 3B and C). The log ratio for *Xff* single infections varied between −1.56 and −0.23. The log-ratio for *Xfm* single infection was 1.15, while that for mixed-strain infection was 0.08. The ANI values confirmed the *Xff* and *Xfm* subspecies for single-infected samples (Fig. S3, Table S2). The mixed sample (*Xff* CFBP 7970 + *Xfm* CFBP 8418) had a higher ANI value for *Xfm*. This result was consistent with a higher number of *Xfm* contigs for the mixed-strain sample (Fig. 3C, Table S2).

Based on BBSplit results, the genetically closest strains sequenced for most of the *Xff* single infections was CFBP 7970, followed by ATCC 35879, GV230, and TPD4 (Table S5). For *Xfm* single infection, the genetically closest *Xfm* strains were Dixon and CFBP 8418. For the mixed infection, 90% of *Xf* contigs were assigned to *Xfm* Dixon and 6% to *Xff* strains.

With the assembled *Xf* contigs, we performed BLAST analysis to identify MLST alleles and virulence genes. We identified the ST number for three of the seven artificially inoculated samples. The mixed-infected sample of grapevine inoculated with the strain CFBP 8418 was identified as ST6 (Table S4). The grapevine sample inoculated with CFBP 7970 (10⁸ CFU/ml) was ST2. We could not assign an ST to the grapevine sample inoculated with the strain CFBP 7970 (10⁷ CFU/ml) because we only identified six of the seven MLST alleles. In contrast to MLST analysis, we detected at least two virulence genes per sample (Table S4). The sample with the lowest CFU values, CFBP 7970 (10⁴ CFU/ml), had 42% similarity to *Xff* hemagglutinin and 41% similarity to *Xfm* pilQ. For the remaining *Xff* single-infected samples, the percentage of similarity to a single subspecies increased with the higher CFU number, which is also associated with higher genome coverage (Table S6). For the *Xfm* single infection, all the virulence genes had 100% similarity to *Xfm*. For the mixed infection, all the virulence genes had 100% similarity to *Xfm* and, on average, 98% to *Xff*.

**Xf identification from field-collected samples.** Finally, we tested the iSeq 100 sequencing capacity with European and American field samples (Table S2). We used 24 samples with $C_q$ values ranging from 21 to 40 based on Harper's qPCR assay. We used three samples that were negative based on the same qPCR assay. The DNAs from the 27 samples were extracted from six different hosts—*Olea europaea*, *Polygala myrtifolia* (France and Italy), *Quercus ilex*, *Spartium junceum*, *Rhamnus alaternus*, and *Vitis vinifera*.

Kraken 2 recovered 0.004 to 1.43% of the total reads as *Xf* (Table S2). We assembled the *Xf* reads into 1 to 2,896 contigs. We found all samples had at least one contig with at least 400 bp. The limit of *Xf* detection with Harper and tetraplex qPCR corresponded to 30 to 37 $C_q$ values (17). Therefore, we evaluated 16 samples that either had fewer than 30 *Xf* contigs, were classified as *Xf*-negative, or had $C_q$ higher than 30 (2). We used each *Xf* contig from the 16 samples as the query for a nucleotide BLAST search using the Web tool from NCBI. Of the 16 samples, 11 gave 100% identity to *Xf* genomes. Hence, these 11 samples were considered *Xf* positive. All contigs from the other five samples, FR1-Pm, FR1-Oe, IT6-Sj, IT11-Sj, and US1-Vv, had 100% identity to chloroplast and 18S plant sequences but none to *Xf*. Therefore, these five samples were considered *Xf* negative (Fig. 4). With our pipeline, we were able to detect *Xf* in samples considered inconclusive by qPCR according to Harper's (2).

We then used the contigs from the 22 *Xf* positive samples for subspecies classification. Overall, the samples had 50% to 100% of contigs classified as *Xf* core contigs (Fig. 4). Three French samples (FR2-Qi, FR2-Pm, FR4-Oe) and five Italian samples (IT5-Sj, IT7-Sj, IT8-Sj, IT9-Sj, IT6-Ra) had 1 to 6 contigs assigned as *Xfm*. The French sample FR3-Oe, seven Italian samples (IT2-Sj, IT3-Sj, IT4-Sj, IT12-Sj, IT4-Ra, IT2-Pm, IT3-Pm), and the U.S. sample, US2-Vv, were *Xfm* single-infected based on the manual log correction (log ratio > 0.43) (Table S2). The sample

**FIG 3** Legend (Continued)

line indicates a logarithmic trendline, $y = 0.4222\ln(x) + 7.9511$; $R^2 = 9.876$. The green numbers indicate the $C_q$ values for each sample. (B) Percentage of *Xf* recovered reads by Kraken 2 from single-infected (*Xff*, 1e8; *Xfm*, 1e7 CFU) and mixed-infected (1e7 CFU) samples. Teal bars indicate *Xf* reads, and orange bars indicate unclassified reads. Unclassified reads show no similarity to *Xf* reads, such as plant or other microorganisms reads. (C) Proportion of *Xf* subspecies from total *Xf* contigs in single and mixed infections. The black dots indicated the log ratio as a manual correction to detect single and mixed infection. *Xfm* and *Xff* are indicated in blue and green, respectively.
<figure><img src="image_13.png" /><figcaption>**FIG 4** *Xf* from field-collected samples and read mapping to subspecies from database. Stacked columns indicate the number of *Xf* contigs with 100% identify to each *Xf* subspecies. Samples indicated in red were *Xf*-negative with our pipeline. *Xff* represents the sum of *Xf* subsp. *fastidiosa*, *Xf* subsp. *morus*, and *Xf* subsp. *sandyi*; *Xfm*, *Xf* subsp. *multiplex*; *Xfp*, *Xf* subsp. *pauca*. "Only core" indicates samples that only have *Xf* core contigs. The *C*<sub>a</sub> values are indicated on the top of each bar. ND, not determined; N, negative for *Xf* based on qPCR. Sample code and hosts are indicated on the x axis. Each country of origin is indicated in the sample ID—France (FR), Italy (IT), and United States (US)—along with the host from which they were isolated.</figcaption></figure>

FR2-Oe had 17% of contigs assigned to *Xfp*, and the sample US3-Vv was *Xff* single-infected (log ratio < 0.08).

Then, we determined the ANI values and *Xf* strain compositions for the 22 *Xf*-positive samples. Both results were consistent with the subspecies identification. The Italian samples had 99 to 100% ANI to *Xfm*. Five French samples (FR1-Qi, FR2-Pm, FR3-Oe, FR4-Oe, FR2-Qi) had 99 to 100% ANI to *Xfm*, and FR2-Oe had 98% ANI to *Xfp* (Fig. 5A, Table S2). The sample US2-Vv had 99% ANI to *Xfm*, and US3-Vv had 100% to *Xff* (Fig. 5A). For strain distribution, all the French and U.S. and three Italian samples (IT2-Sj, IT4-Sj, IT5-Sj) had more than 40% *Xf* contigs had 100% identity to one strain (Table S5, Fig. 5B). Except for IT2-Sj, the Italian samples had most of the contigs assigned to the three strains TOS4, TOS5, and TOS14. The sample IT2-Sj had more reads assigned to *Xfm* RAAR14. Overall, the tools SendSketch, Pyani, and Bbsplit validated the qPCR subspecies results for field samples.

We performed BLAST analysis to identify MLST alleles and virulence genes for the 22 infected samples. For 16 of 22 *Xf*-positive samples, we found the percentage of gene similarity to be 26 to 100% for at least one virulence gene (Table S4). Distinct from single-gene analysis, we only identified some MLST-related alleles for four samples (US3-Vv, FR2-Oe, IT2-Pm, IT3-Pm); consequently, we could not identify the ST number (Table S3).

To compare some of our results with a high-performance, deep-sequencing Illumina platform as a control, we selected nine samples for resequencing with the MiSeq platform
<figure><img src="image_14.png" /><figcaption>**FIG 5** Metagenomics analyses of field sample identifies of bacterial subspecies. (A) The dendrogram indicates the distance and cluster analysis based on ANI values using NCBI whole genomes and assembled *Xf* contigs. Branch colors represent each *Xf* subspecies. Blue, gray, and green names indicate iSeq 100-sequenced samples. (B) The heatmap shows the percentage of unique contigs assigned to each *Xf* strain. The cluster and strain distribution were plotted using ComplexHeatmap R package.</figcaption></figure>

using the same iSeq 100 libraries from this study. The nine samples were IT3-Pm, IT5-Sj, FR2-Oe, FR3-Oe, US1-Vv, US2-Vv, FR1-Pm, FR2-Pm, and IT9-Sj (Table S2, Fig. S4). We analyzed the MiSeq sequences with our pipeline and recovered 0.005 to 0.792% of total reads as *Xf* reads with Kraken 2. We assembled the *Xf* reads into contigs and manually assessed all samples with less than 30 *Xf* contigs. The NCBI blastn analysis indicated that the samples US1-Vv and FR1-Pm had contigs with 100% identity to plant reads; therefore, we confirmed they were *Xf*-negative samples. The other seven samples were considered *Xf*-positive. We followed the pipeline to identify and determine subspecies, phylogeny, genetically closest sequenced genome strains, MLST profile, and virulence-related genes. The results for subspecies and phylogeny identification were the same between MiSeq and iSeq 100 sequencers (Fig. S4), but there were some differences for the other three analysis results (Table S2). For the genetically closest sequenced genome strain analysis, all samples gave the same strain distribution as iSeq 100 results, except FR2-Oe, which showed *Xfp* OLS0478 instead of *Xfp* COF0407 as the most abundant strain. These two *Xfp* strains are phylogenetically close. For MLST analysis, we identified four more alleles for the sample IT3-Pm with the MiSeq platform than with iSeq 100 (Table S3), while we only detected two alleles in the sample US2-Vv sequenced with Miseq. We were only able to detect MLST alleles for the sample FR2-Oe with the iSeq 100 platform. We were able to calculate the percentage of gene similarity for more virulence genes with the MiSeq platform than with the iSeq 100 platform. The variation between Illumina platforms was not consistent.

DISCUSSION

In this study, we developed a user-friendly metagenomic pipeline to identify and determine *Xylella fastidiosa* subspecies from field-collected samples without the need for pathogen isolation. We demonstrated the flexibility of the pipeline by using seven different plant hosts and three DNA extraction methods. We recovered and assembled *Xf* reads into contigs from total DNA samples. We used percentage of similarity to a single subspecies to identify *Xf* subspecies and validated the results through phylogeny and
strain proximity. Finally, we examined potential virulence-related genes among all sequenced samples.

To recover *Xf* reads from field samples with the tool Kraken 2, we used *Xf* genomes available on NCBI. We found plant plastid reads in all genomes obtained from pure *Xf* cultures, except for *Xff* CFBP7970. We decided to add plant reads in the database to filter out potential plant contamination. We still recovered *Xf* reads for some plant samples reported as *Xf*-negative. Therefore, we manually examined the contigs of samples with fewer than 30 *Xf* contigs to determine if they have a low *Xf* concentration or are negative samples. More than the *Cq* values, the contig evaluation will be necessary to determine if a sample is truly *Xf*-negative.

We observed that the Kraken 2 tool not only recovered *Xf* reads but also identified *Xf* subspecies. We decided not to use Kraken2 to identify subspecies because we found that the recovery is affected by incomplete NCBI *Xf* genome subspecies information. Kraken 2 uses by default the NCBI taxonomy to classify reads; if the genomes used to build the database do not have subspecies information, it will keep most of the reads at the species level. To improve the subspecies resolution, we decided to use contigs instead of reads and the tool SendSketch (BBMap tool). Contigs or assembled reads increase the coverage and reduce false-positive reads (28). We used SendSketch because it uses the MinHash algorithm to be fast and it takes into account whole genomes and does not use taxonomy.

After we identified samples as *Xf*-positive, we defined *Xf* subspecies. We observed that some samples had mapped contigs to both *Xff* and *Xfm*. These are contigs most likely associated with core sequences, as only 3% of *Xff* and *Xfm* genomes are different. However, it is also possible to have a percentage of annotation error due to sequencing contamination (29, 30). To determine if the samples were single- or mixed-infected, we corrected the results by calculating the log *Xfm:Xff* contig ratio (see Materials and Methods). Once we defined the log values for *Xff*, *Xfm*, and *Xfm+Xff* (mixed) infections, we validated the presence of single *Xfm* infections in all the tested European samples. Identifying *Xf* subspecies in ornamental and crop plants is essential for the correct application of eradication measures or for plant movements within the European Union territory according to regulation (EU) 2020/1201 (31).

The number of reads generated by the iSeq 100 sequencer highlighted some limitations of the pipeline. For example, the ST were not determined for any field sample because we could not recover the complete sequences of the seven genes. This is probably associated with low genome coverage. The low number of reads also hampers deep SNP diversity, inter-subspecific homologous recombination analyses, and microbiome studies. Other sequencing systems, with higher read output than the iSeq 100, can also be used with this pipeline, as the input is fastq files. With a high number of reads, the pipeline will provide better resolution to recover the MSLT genes and have enough sequencing depth to describe the microbial diversity in the sample. Complete diagnostic analysis should include microbiome diversity analysis to provide the context of microbial interactions during disease.

Some of the diagnostic tools for *Xf* diagnostics are qPCR and Sanger sequencing. These tools require amplification of known *Xf* genome regions but, as is the case of MLST, do not consider new or emerging variants. Moreover, these tools introduce bias due to primer design, have unresolved results with high *Cq* values, and may take longer since they follow a multistep process. Our pipeline complements these conventional tools by obtaining metagenomic data directly from symptomatic or asymptomatic samples and increasing the detection power. We found some discrepancies between the number of recovered *Xf* reads and *Cq* values. These differences could be caused by PCR inhibitors or the genomic target region for qPCR that underestimates the bacterial concentration (32, 33). Metagenomics sequencing is becoming a more affordable and faster approach for diagnostics. For example, the whole detection/identification with qPCR and MLST scheme could have an estimated cost of $52 to $54 per sample and takes three to 4 days to detect one to seven genes. With the iSeq 100, it could cost $50 to $70 (when having 12 samples in the same run) and take 2 days but while also allowing a complete genomic analysis of the plant and its pathogenic
and commensal microbiota. Metagenomics for pathogen diagnostics will serve as a model for future diagnostics programs, and we can eventually expect this to become readily accessible to teams across Europe as sequencing becomes more affordable. Our work sets the stage in preparedness for this event. We know that many European laboratories may not have access to metagenomics for pathogen identification. We feel that this and other approaches can serve as a platform for epidemic preparedness.

We are aware that long-read sequencing technologies are becoming affordable, accessible technologies for doing disease diagnostics. One of these platforms is from Oxford Nanopore Technologies (ONT). ONT allows for long-read sequencing of lengths over 5 kb, in contrast to short-read read sequencing, such as Illumina, which has a maximum of 300 to 500 bp. Long-read sequencing, compared to short-read, provides longer contiguous sequences, which is critical for important repetitive virulence factors or insertion elements. Some disadvantages of long-read compared to short-read sequencing include (i) fewer total reads generated and (ii) often higher inherent error rates due to the sequencing chemistry (34). We are concerned that fewer reads would not provide sufficient coverage to multiplex samples and would capture low concentrations of pathogen reads, while Illumina has more reads to reach a lower limit of detection. Higher error rates could affect classification methods for *Xf* subspecies, as their genomes display less than 3% difference for ANI. However, the technologies are swiftly advancing with new ONT chemistries (R10) that reduce the sequencing errors competitively with Illumina (34). Overall, a complete diagnostic approach should include long- and short-read sequencing to provide at the same time information about the genomes, pathogen abundance, and nucleotide changes.

In conclusion, our pipeline provides *Xf* taxonomy and functional information for diagnostics without extensive knowledge of the host or pathogen. The pipeline databases used for the analysis can be public repositories or privately collected gDNA and could be adapted by the user and tailored to different plant pathogens. The sequencing can be adapted to be an in-house system, as the library preparation and sequencing are user-friendly and not limited by the DNA quality or quantity. The analysis can be adjusted to detect several pathogens simultaneously. Our pipeline can be used for early detection of *Xf* or other crop pathogens and be incorporated as part of management strategies.

## MATERIALS AND METHODS

**Xylella fastidiosa strains.** The strain *Xf* subsp. *fastidiosa* (*Xff*) CFBP 7970, isolated in the United States (Florida) in 2013 from *Vitis vinifera*, and *Xf* subsp. *multiplex* (*Xfm*) CFBP 8418, isolated in France in 2015 from *Spartium junceum*, were provided by the French Collection of Plant-Associated Bacteria (CIRM-CFBP; https://www6.inra.fr/cirm_eng/CFBP-Plant-Associated-Bacteria) and used as controls for whole-genome sequencing. Both strains were cultivated on modified PWG medium (Gelrite, 12 g · liter⁻¹; soytone, 4 g · liter⁻¹; Bacto tryptone, 1 g · liter⁻¹; MgSO₄·7H₂O, 0.4 g · liter⁻¹; K₂HPO₄·1.2 g · liter⁻¹; KH₂PO₄·1 g · liter⁻¹; hemin chloride [0.1% in NaOH, 0.05 M], 10 ml · liter⁻¹; bovine serum albumin [BSA] [7.5%], 24 ml · liter⁻¹; L-glutamine, 4 g · liter⁻¹) at 28°C for 1 week.

**Artificially inoculated plant samples.** For artificial inoculations, 10 ml of a calibrated CFBP 7970 or CFBP 8418 strain suspension was spiked in 2 g of detached *V. vinifera* leaves. Sterile water was used for negative controls. The DNA extraction was performed using a CTAB-based extraction protocol (2) with slight modifications in order to concentrate bacterial DNA. After a 20-min centrifugation of the sample at 20,000 × g, the pellet was resuspended in 1 ml of CTAB buffer. At the end of the extraction protocol, the pellet was resuspended in 50 µl of sterile demineralized water. *Xf* presence in the infected samples was checked using Harper's qPCR assay (8).

**Plant material and bacterial gDNA.** Healthy plant material of *Vitis vinifera* (2 g of leaf petioles) was spiked with 10 ml of a calibrated CFBP 7970 or CFBP 8418 strain suspension. Sterile water was used as the negative control. The DNA extraction was performed using a CTAB-based extraction protocol (2) with slight modifications in order to concentrate bacterial DNA. After a 20-min of centrifugation at 20,000 × g, the plant macerate sample was resuspended in 1 ml of CTAB buffer. At the end of the extraction protocol, the pellet was resuspended in 50 µl of sterile demineralized water. *Xf* presence in the infected samples was checked using Harper's qPCR assay (8).

Naturally infected samples were collected in Europe and the United States. Symptomatic samples of *Olea europaea*, *Polygala myrtifolia*, and *Quercus ilex* were collected in October 2018 in Corsica (France) and in September 2019 in the French Riviera. *Xf* detection and DNA extraction of whole infected plant tissue were performed as mentioned above. *Xf* subspecies were identified using the tetraplex qPCR (17).

Twig tissues, leaf petioles, or green shoots of *Rhamnus alaternus*, *Spartium junceum*, and *P. myrtifolia* growing in the *Xf* outbreak zone of Monte Argentario (Grosseto, Tuscany, Italy) were collected during 2019 and 2020 (35). *Xf* was detected using Harper's qPCR assay (8). For samples with a $C_q$ value lower than 30, *Xf* isolation was attempted on buffered charcoal yeast extract (BCYE) agar according to PM 7/24-4 (2). Bacterial isolates
that became visible to the unaided eye within 3 days of incubation at 28°C were discarded; those that became visible thereafter were streaked twice for purity on BCYE agar and identified as *Xf* based on qPCR results (8). Reactions were carried out after boiling the bacterial suspension for 10 min. The DNA of one isolate among those that tested positive by qPCR from each plant was extracted using the CTAB-based protocol and further characterized to the subspecies and ST levels following the multilocus sequence typing (MLST) approach (36). The GoTaq probe qPCR master mix (Promega; A6102) and GoTaq G2 (Promega; M784B) polymerase were used for qPCR and conventional PCR experiments, respectively. The bacterial DNA of three isolates from *R. alaternus* (IT1-Ra to IT3-Ra), one from *S. junceum* (IT1-Sj), one from *P. myrtifolia* (IT1-Pm), and one from *Prunus dulcis* (IT1-Pd) were sequenced in this study.

**Vitis vinifera** DNA samples were received from the Virginia Tech Plant Disease Clinic (Virginia, USA). All samples were collected from Virginia vineyards in 2019. The US2-Vv sample was collected from a vineyard in Greene County, and the US3-Vv sample was from a vineyard in Isle of Wight County. The DNA extraction and *Xf* detection protocols are based on instructions from the Virginia Tech Plant Disease Clinic (VTPDC). Approximately 50 to 100 mg of grape leaf or petiole tissue was excised from each sample using a razor blade. The excised tissue was transferred to lysing matrix A tubes (MP Biomedicals; 6910-500) and ground using a FastPrep 24 device (MP Biomedicals; 116004500). DNA was extracted using an Isolate II plant DNA extraction kit (Bioline; 52070) the following manufacturer's recommendations and CTAB lysis buffer. For *Xf* detection, Harper's qPCR was performed on the StepOnePlus system (Life Technologies; 4376600) with a Sensi-FAST probe Hi-ROX qPCR kit (Bioline; 82005) (8).

**Pipeline controls.** To test the pipeline, seven samples were used as negative controls. The controls were two DNA samples from healthy barley and wild grass leaves that were grown in a greenhouse; two DNA samples from barley leaves were infiltrated with a bacterial suspension (10⁸ CFU/ml) of *Xanthomonas translucens* pv. translucens UPB886. Three healthy samples were collected in France in 2020 and in the United States in 2019. Petioles and midribs were collected from healthy *Olea europaea* plants in a non-*Xf*-infected area (Angers, France) and from *Polygala myrtifolia* plants that were purchased from a local nursery. *V. vinifera* leaves were collected from the vineyard in Greene County, VA, USA. The DNA from these samples was extracted using the CTAB method as mentioned above. The absence of *Xf* in the healthy plants was confirmed using Harper's qPCR assay (8).

For *in silico* pipeline controls, two types of positive controls were used. To validate the detection of different concentrations of *Xf*, barley fasta sequence files were *in silico* mixed with reads of the *Xff* CFBP 7970 sequenced in this study. The final proportion of *Xff* reads in the sample ranged from 0.2 to 2.4% of total reads. To validate the detection limits for single and mixed infections, the barley fasta sequence files were *in silico* mixed with different proportions of *Xff* CFBP 7970 and *Xfm* CFBP 8418 reads, to get from percentage CFBP7970:99% CFBP8418 to 99% CFBP7970:1% CFBP8418.

**iSeq 100 sequencing.** iSeq 100 sequencing libraries were prepared according to the Illumina reference guide for the Nextera DNA Flex library prep and Nextera DNA CD indexes. In brief, 200 to 500 ng of DNA was quantified by spectrophotometry and used for library preparation. Then, the libraries were diluted to have the same starting concentration prior to sample pooling. Eight to twelve libraries were mixed together, and 1 nM pooled library was used per run. The sequencing settings were paired-ended (PE) read type, 151 read cycles, and 8 index cycles. In the iSeq 100 system, the Illumina Generate FASTQ analysis module for base calling and demultiplexing was selected. After sequencing for 17 h, fastq paired-end read files were extracted from the machine for subsequent analysis.

**Pipeline for *Xylella* sp. detection, classification, and quantification via metagenomic analysis.** Fastq files and the program Kraken 2 were used to recover *Xylella* reads (22). The Kraken 2 command options were -paired, -minimum-hit-groups 5, -report, and -db. The database was created with 92 NCBI genomes, 79 from *Xf*, 2 from *Xylella taiwanensis*, 10 from *Xanthomonas* spp., and 1 from *Escherichia coli* (Table S1). The 81 *Xylella* genomes were used to recover *Xylella* reads. The last 11 genomes were added to remove reads common to *Proteobacteria* that might give a false-positive match. The tool SendSketch with the nucleotide (nt) server was used to make sure the 81 NCBI *Xf* genomes were not contaminated with plant reads (BBMap; B. Bushnell; [sourceforge.net/projects/bbmap/](http://sourceforge.net/projects/bbmap/); 30 October 2019). A total of 49 NCBI sequences from plant 18S and chloroplast were added to the Kraken 2 customized database to avoid extracting reads annotated as plant reads (Table S1).

The reads classified as *Xylella* were extracted with the script *extract_kraken_reads* from the KrakenTools suite ([https://github.com/jennifertu717/KrakenTools](https://github.com/jennifertu717/KrakenTools)). The extracted *Xf* reads were used for downstream analysis. First, *Xf* reads were *de novo* assembled with the software SPAdes (23) using default settings and the option -only-assembler. Second, the *Xf* contigs were the query sequences in the Basic Local Alignment Search Tool website (NCBI) to confirm if they were *Xf* reads or misclassified plant reads. The blastn parameters were nucleotide collection (nonredundant/nucleotide [nr/nt]) as database and megaBLAST program selection.

The *Xf*-positive contigs were used in four different analyses: (i) to identify subspecies, (ii) to reconstruct phylogeny, (iii) to identify the genetically closest strains already sequenced, and (iv) to identify alleles from specific genes and the MLST profile.

To identify subspecies, the tool SendSketch was run with the parameters, mode=sequence, records=2, format=3, minani=100, minhit=1, address=ref, and level=0 (BBMap; B. Bushnell B.; [sourceforge.net/projects/bbmap/](http://sourceforge.net/projects/bbmap/)). Only contigs with 100% average nucleotide identity (ANI) were used to identify *Xf* subspecies. The *Xf* contigs with no hits were considered core sequences. For visualization, results were plotted using stacked bars.

To reconstruct phylogeny, ANI was calculated using the software Pyani (0.2.10) and LINbase (24, 25). For Pyani, the option -ANIm was set, and for LINbase, "Identify using a gene sequence" was set as the identification method. The R package ComplexHeatmap was used to visualize the ANI cluster analysis with the parameters clustering_distance_rows = robust_dist and clustering_method_rows = "average." Robust_dist was
a function suggested by the ComplexHeatmap Complete Reference (37). We confirmed the scale using hclust from the R library stats.

To identify the genetically closest already sequenced *Xf* strains, the tool BBSplit was run with *Xf*-contigs (BBMap; B. Bushnell; sourceforge.net/projects/bbmap/). The tool used 81 *Xf* genomes from NCBI and the options minratio=1 and ambig=best. For comparisons, each sample was normalized to its total *Xf* contigs and plotted using the R package ComplexHeatmap. For the isolated bacterial genomes, the most abundant strains were used as references to identify genomics variants using Harvest suite tools (26).

To determine specific gene alleles, the percentage of identity was calculated using *Xf* contigs as the query and the blastn algorithm (nucleotide-nucleotide BLAST 2.8.1+). The database contained complete nucleotide sequences for all alleles for the seven genes used for ST identification (*cysG*, *gltT*, *leuA*, *malF*, *nuoL*, *holC*, and *petC*) (36). All 147 alleles were downloaded from the website PubMLST (38) (last updated, 2019-03-06). To determine the presence of reported *Xf* virulence-related genes or those common to several plant-pathogenic bacteria (39, 40), the percentage of similarity was calculated using *Xf* contigs as the query using the blastx algorithm (nucleotide-nucleotide BLAST 2.8.1+). The databases contained complete amino acid sequences for *gumBCDE*, *pilBMQTVW*, *rpfCEFG*, *tolC*, 6-phosphogluconolactonase (*pgl*), and hemagglutinin from the *Xfm* M12 and *Xff* M23 NCBI genomes. For both blastn and blastx scripts the default settings were used, and two were modified (evaluation 0.1; outfmt "6 qseqid sseqid pident evalue").

**MiSeq sequencing.** To validate the iSeq 100 results, the same iSeq 100 libraries were used for MiSeq deep sequencing. Nine samples were selected, at least one from each iSeq 100 run, making sure not to use the same i5 and i7 tags. These nine libraries were sent to the Animal Disease Diagnostic Laboratory (Ohio Department of Agriculture, Reynoldsburg, Ohio) for sequencing. Library preparation was performed using an Illumina DNA Flex kit, and 2 × 250 sequencing was performed on the MiSeq platform using V3 chemistry. The pipeline described above was used to analyze the MiSeq fastq files.

**Data availability.** The assembled genomes of the six *Xf* isolates were deposited in NCBI under the BioProject number PRJNA728043 and BioSample numbers SAMN19067399, SAMN19067400, SAMN19067398, SAMN19067401, SAMN19067403, and SAMN19067402. The scripts for Kraken 2 can be found at https://github.com/DerrickWood/kraken2/wiki/Manual; the scripts to analyze Kraken 2 outputs can be found at https://github.com/jenniferlu717/KrakenTools; the scripts for SendSketch and BBSplit are from the user guide at https://jgi.doe.gov/data-and-tools/bbtools/bb-tools-user-guide/. To calculate the SNPs for each genome, we used the Harvest manual at https://harvest.readthedocs.io/en/latest/. For building the trees, we used the scripts from the ComplexHeatmap reference manual at https://jokergoo.github.io/ComplexHeatmap-reference/book/a-single-heatmap.html. All iSeq 100 and MiSeq raw reads can be made available upon request.

## SUPPLEMENTAL MATERIAL

Supplemental material is available online only.

**FIG S1**, EPS file, 2.4 MB.

**FIG S2**, EPS file, 0.9 MB.

**FIG S3**, EPS file, 0.7 MB.

**FIG S4**, EPS file, 1.1 MB.

**TABLE S1**, XLSX file, 0.02 MB.

**TABLE S2**, XLSX file, 0.03 MB.

**TABLE S3**, XLSX file, 0.01 MB.

**TABLE S4**, XLSX file, 0.02 MB.

**TABLE S5**, XLSX file, 0.02 MB.

**TABLE S6**, XLSX file, 0.01 MB.

## ACKNOWLEDGMENTS

We are grateful for funding support from the Ohio Department of Agriculture Specialty Crops block grant (AGR-SCG-19-03) and USDA NIFA FACT (2021-67021-34343) to J.M.J. Sara Campigli was financed by a grant from the Phytosanitary Service of the Tuscany Region (Italy).

We thank The Ohio Supercomputer Center for providing high performance computing resources.

We also thank Stephen Cohen, Jeff Chang, and Alexandra Weinsberg for their valuable input for the experiment development and editing.

## REFERENCES

1. EFSA. 2020. Update of the Xylella spp. host plant database: systematic literature search up to 30 June 2019. EFSA J 18:e06114. https://doi.org/10.2903/j.efsa.2020.6114.

2. EPPO. 2019. PM 7/24 (4) Xylella fastidiosa. EPPO Bull 49:175-227. https://doi.org/10.1111/epp.12575.

3. Chatterjee S, Wistrom C, Lindow SE. 2008. A cell-cell signaling sensor is required for virulence and insect transmission of Xylella fastidiosa. Proc Natl Acad Sci U S A 105:2670-2675. https://doi.org/10.1073/pnas.0712236105.

4. Sabella E, Aprile A, Genga A, Siciliano T, Nutricati E, Nicoli F, Vergine M, Negro C, De Bellis L, Luvisi A. 2019. Xylem cavitation susceptibility and
refilling mechanisms in olive trees infected by Xylella fastidiosa. *Sci Rep* 9: 9602. https://doi.org/10.1038/s41598-019-46092-0.

5. Pierce NB, Newton B. 1892. The California vine disease: a preliminary report of investigations. G.P.O: Washington, DC.

6. EFSA. 2018. Updated pest categorisation of Xylella fastidiosa. EFSA J 16: 5357. https://doi.org/10.2903/j.efsa.2018.5357.

7. Francis M, Lin H, Rosa JC-L, Doddapaneni H, Civerolo EL. 2006. Genome-based PCR primers for specific and sensitive detection and quantification of Xylella fastidiosa. Eur J Plant Pathol 115:203-213. https://doi.org/10.1007/s10658-006-9009-4.

8. Harper SJ, Ward LI, Clover GRG. 2010. Development of LAMP and real-time PCR methods for the rapid detection of Xylella fastidiosa for quarantine and field applications. Phytopathology 100:1282-1288. https://doi.org/10.1094/PHYTO-06-10-0168.

9. Li W, Teixeira DC, Hartung JS, Huang Q, Duan Y, Zhou L, Chen J, Lin H Lopes S, Ayres AJ, Levy L. 2013. Development and systematic validation of qPCR assays for rapid and reliable differentiation of Xylella fastidiosa strains causing citrus variegated chlorosis. J Microbiol Methods 92:79-89. https://doi.org/10.1016/j.mimet.2012.10.008.

10. Ouyang P, Arif M, Fletcher J, Melcher U, Corona FMO. 2013. Enhanced reliability and accuracy for field deployable bioforensic detection and discrimination of Xylella fastidiosa subsp. paucu, causal agent of citrus variegated chlorosis using Razor Ex technology and TaqMan quantitative PCR. PLoS One 8:e81647. https://doi.org/10.1371/journal.pone.0081647.

11. Ito T, Suzaki K. 2017. Universal detection of mycoplasmas and Xylella spp. by TaqMan singleplex and multiplex real-time PCR with dual priming oligonucleotides. PLoS One 12:e0185427. https://doi.org/10.1371/journal.pone.0185427.

12. Bonants P, Griekspoor Y, Houwers I, Krijger M, van der Zouwen P, van der Lee TAJ, van der Wolf J. 2019. Development and evaluation of a triplex TaqMan assay and next-generation sequence analysis for improved detection of Xylella in plant material. Plant Dis 103:645-655. https://doi.org/10.1094/PDIS-08-18-1433-RE.

13. Schaad NW, Postnikova E, Lacy G, Fatmi M, Chang C-J. 2004. Xylella fastidiosa subspecies: X. fastidiosa subsp. piercii, subsp. nov., X. fastidiosa subsp. multiplex subsp. nov., and X. fastidiosa subsp. paucu subsp. nov. Syst Appl Microbiol 27:290-300. https://doi.org/10.1078/0723-2020-00263.

14. Nunney L, Schuenzel EL, Scally M, Bromley RE, Stouthamer R. 2014. Large-scale intersubspecific recombination in the plant-pathogenic bacterium Xylella fastidiosa Is associated with the host shift to mulberry. Appl Environmental Microbiol 80:3025-3033. https://doi.org/10.1128/AEM.04112-13.

15. Denancé N, Briand M, Gaborieau R, Gaillard S, Jacques M-A. 2019. Identification of genetic relationships and subspecies signatures in Xylella fastidiosa. BMC Genomics 20:239. https://doi.org/10.1186/s12864-019-5565-9.

16. Marcelletti S, Scortichini M. 2016. Genome-wide comparison and taxonomic relatedness of multiple Xylella fastidiosa strains reveal the occurrence of three subspecies and a new Xylella species. Arch Microbiol 198: 803-812. https://doi.org/10.1007/s00203-016-1245-1.

17. Dupas E, Briand M, Jacques M-A, Cesbrun S. 2019. Novel tetraplex quantitative PCR assays for simultaneous detection and identification of Xylella fastidiosa subspecies in plant tissues. Front Plant Sci 10:1732. https://doi.org/10.3389/fpls.2019.01732.

18. Denancé N, Legendre B, Briand M, Olivier V, de Boisseson C, Poliakoff F, Jacques M-A. 2017. Several subspecies and sequence types are associated with the emergence of Xylella fastidiosa in natural settings in France. Plant Pathol 66:1054-1064. https://doi.org/10.1111/ppa.12695.

19. Jones S, Baizan-Edge A, MacFarlane S, Torrance L. 2017. Viral diagnostics in plants using next generation sequencing: computational analysis in practice. Front Plant Sci 8:1770. https://doi.org/10.3389/fpls.2017.01770.

20. Mechan Llontop ME, Sharma P, Aguilera Flores M, Yang S, Pollok J, Tian L, Huang C, Rideout S, Heath LS, Li S, Vinatzer BA. 2020. Strain-level identification of bacterial tomato pathogens directly from metagenomic sequences. Phytopathology 110:768-779. https://doi.org/10.1094/PHYTO-09-19-0351-R.

21. Faino L, Scala V, Albanese A, Modesti V, Grottoli A, Pucci N, L'Aurora A, Reverberi M, Loreti S. 2019. Nanopore sequencing for the detection and the identification of Xylella fastidiosa subspecies and sequence types from naturally infected plant material. bioRxiv https://doi.org/10.1111/ppa.13416.

22. Wood DE, Salzberg SL. 2014. Kraken: ultrafast metagenomic sequence classification using exact alignments. Genome Biol 15:R46. https://doi.org/10.1186/gb-2014-15-3-r46.

23. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, Pyshkin AV, Sirotkin AV, Vyahhi N, Tesler G, Alekseyev MA, Pevzner PA. 2012. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol 19:455-477. https://doi.org/10.1089/cmb.2012.0021.

24. Tian L, Huang C, Mazloom R, Heath LS, Vinatzer BA. 2020. LINbase: a Web server for genome-based identification of prokaryotes as members of crowdsourced taxa. Nucleic Acids Res 48:W529-W537. https://doi.org/10.1093/nar/gkaa190.

25. Pritchard L, Glover RH, Humphris S, Elphinstone JG, Toth IK. 2016. Genomics and taxonomy in diagnostics for food security: soft-rotting enterobacterial plant pathogens. Anal Methods 8:12-24. https://doi.org/10.1039/C5AY02550H.

26. Treangen TJ, Ondov BD, Koren S, Phillippy AM. 2014. The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes. Genome Biol 15:524. https://doi.org/10.1186/s13059-014-0524-x.

27. Giampetruzzi A, D'Attoma G, Zicca S, Abou Kubaa R, Rizzo D, Boscia D, Saldarelli P, Saponari M. 2019. Draft genome sequence resources of three strains (TOS4, TOS5, and TOS14) of Xylella fastidiosa infecting different host plants in the newly discovered outbreak in Tuscany, Italy. Phytopathology 109:1516-1518. https://doi.org/10.1094/PHYTO-04-19-0108-A.

28. Ayling L, Clark MD, Leggett RM. 2020. New approaches for metagenome assembly with short reads. Brief Bioinform 21:584-594. https://doi.org/10.1093/bib/bbz020.

29. Nunney L, Elfekih S, Stouthamer R. 2012. The importance of multilocus sequence typing: cautionary tales from the bacterium Xylella fastidiosa. Phytopathology 102:456-460. https://doi.org/10.1094/PHYTO-10-11-0298.

30. Kircher M, Sawyer S, Meyer M. 2012. Double indexing overcomes inaccuracies in multiplex sequencing on the Illumina platform. Nucleic Acids Res 40:e3. https://doi.org/10.1093/nar/gkr771.

31. Nunney L, Vickerman DB, Bromley RE, Russell SA, Hartman JR, Morano LD, Stouthamer R. 2013. Recent evolutionary radiation and host plant specialization in the Xylella fastidiosa subspecies native to the United States. Appl Environmental Microbiol 79:2189-2200. https://doi.org/10.1128/AEM.03208-12.

32. Modesti V, Pucci N, Lucchesi S, Campus L, Loreti S. 2017. Experience of the Latium region (Central Italy) as a pest-free area for monitoring of Xylella fastidiosa: distinctive features of molecular diagnostic methods. Eur J Plant Pathol 148:557-566. https://doi.org/10.1007/s10658-016-1111-7.

33. Schrader C, Schielke A, Ellerbroek L, Johne R. 2012. PCR inhibitors: occurrence, properties and removal. J Appl Microbiol 113:1014-1026. https://doi.org/10.1111/j.1365-2672.2012.05384.x.

34. Karst SM, Ziels RM, Kirkegaard RH, Sørensen EA, McDonald D, Zhu Q, Knight R, Albertsen M. 2021. High-accuracy long-read amplicon sequences using unique molecular identifiers with Nanopore or PacBio sequencing. Nat Methods 18:165-169. https://doi.org/10.1038/s41592-020-01041-y.

35. Marchi G, Rizzo D, Ranaldi F, Ghelardini L, Ricciolini M, Scarpelli I, Drosera L, Goti E, Capretti P, Surico G. 2018. First detection of Xylella fastidiosa subsp. DNA in Tuscany (Italy). Phytopathol Mediterr 57: 363-364. https://doi.org/10.14601/Phytopathol_Mediterr-24454.

36. Yuan X, Morano L, Bromley R, Spring-Pearson S, Stouthamer R, Nunney L. 2010. Multilocus sequence typing of Xylella fastidiosa causing Pierce's disease and oleander leaf scorch in the United States. Phytopathology 100:601-611. https://doi.org/10.1094/PHYTO-100-6-0601.

37. Gu Z, Eils R, Schlesner M. 2016. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics 32:2847-2849. https://doi.org/10.1093/bioinformatics/btw313.

38. Jolley KA, Bray JE, Maiden MCJ. 2018. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res 3:124. https://doi.org/10.12688/wellcomeopenres.14826.1.

39. da Silva FR, Vettore AL, Kemper EL, Leite A, Arruda P. 2001. Fastidian gum: the Xylella fastidiosa exopolysaccharide possibly involved in bacterial pathogenicity. FEMS Microbiol Lett 203:165-171. https://doi.org/10.1111/j.1574-6968.2001.tb10836.x.

40. Chatterjee S, Almeida RPP, Lindow S. 2008. Living in two worlds: the plant and insect lifestyles of Xylella fastidiosa. Annu Rev Phytopathol 46:243-271. https://doi.org/10.1146/annurev.phyto.45.062806.094342.