The importance of multilocus sequence typing: cautionary tales from the bacterium Xylella fastidiosa

Metagenomic Sequencing for Identification of Xylella fastidiosa from Leaf Samples

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 [C_q], >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.

Genetic Diversity of Xylella fastidiosa in Mexican Vineyards

Xylella fastidiosa is a xylem-inhabiting phytopathogenic bacterium that affects diverse agriculturally relevant crops. In Mexico, X. fastidiosa has been reported in the states of Baja California, Coahuila, and Querétaro. In order to determine the genetic diversity of this bacterium in Mexico, 408 grapevine samples were collected from the main producing states in México. For X. fastidiosa identification, real-time PCR and three-loci end-point PCR were employed. The genotyping of the subspecies was carried out using multilocus sequence typing and analysis, based on seven housekeeping genes: leuA, petC, malF, cysG, holC, nuoL, and gltT. The resulting sequences were compared with those present in extant databases. The presence of X. fastidiosa subsp. fastidiosa in the states of Baja California (sequence type 1), Coahuila (sequence type 1), and Querétaro was confirmed. The isolates from northern Mexico bear high similarity to grapevine isolates from the United States. However, the isolates from Querétaro showed significant differences with currently known sequences, showing that there is genetic variability among the X. fastidiosa subsp. fastidiosa populations from grapevines in northern and central Mexico.

Xylella fastidiosa: climate suitability of European continent

The bacterium Xylella fastidiosa (Xf) is a plant endophyte native to the Americas that causes diseases in many crops of economic importance (grapevine, Citrus, Olive trees etc). Xf has been recently detected in several regions outside of its native range including Europe where little is known about its potential geographical expansion. We collected data documenting the native and invaded ranges of the Xf subspecies fastidiosa, pauca and multiplex and fitted bioclimatic species distribution models (SDMs) to assess the potential climate suitability of European continent for those pathogens. According to model predictions, the currently reported distribution of Xf in Europe is small compared to the large extent of climatically suitable areas. The regions at high risk encompass the Mediterranean coastal areas of Spain, Greece, Italy and France, the Atlantic coastal areas of France, Portugal and Spain as well as the southwestern regions of Spain and lowlands in southern Italy. The extent of predicted climatically suitable conditions for the different subspecies are contrasted. The subspecies multiplex, and to a certain extent the subspecies fastidiosa, represent a threat to most of Europe while the climatically suitable areas for the subspecies pauca are mostly limited to the Mediterranean basin. These results provide crucial information for the design of a spatially informed European-scale integrated management strategy, including early detection surveys in plants and insect vectors and quarantine measures.

Identification of genetic relationships and subspecies signatures in Xylella fastidiosa

BACKGROUND

The phytopathogenic bacterium Xylella fastidiosa was thought to be restricted to the Americas where it infects and kills numerous hosts. Its detection worldwide has been blooming since 2013 in Europe and Asia. Genetically diverse, this species is divided into six subspecies but genetic traits governing this classification are poorly understood.

RESULTS

SkIf (Specific k-mers Identification) was designed and exploited for comparative genomics on a dataset of 46 X. fastidiosa genomes, including seven newly sequenced individuals. It was helpful to quickly check the synonymy between strains from different collections. SkIf identified specific SNPs within 16S rRNA sequences that can be employed for predicting the distribution of Xylella through data mining. Applied to inter- and intra-subspecies analyses, it identified specific k-mers in genes affiliated to differential gene ontologies. Chemotaxis-related genes more prevalently possess specific k-mers in genomes from subspecies fastidiosa, morus and sandyi taken as a whole group. In the subspecies pauca increased abundance of specific k-mers was found in genes associated with the bacterial cell wall/envelope/plasma membrane. Most often, the k-mer specificity occurred in core genes with non-synonymous SNPs in their sequences in genomes of the other subspecies, suggesting putative impact in the protein functions. The presence of two integrative and conjugative elements (ICEs) was identified, one chromosomic and an entire plasmid in a single strain of X. fastidiosa subsp. pauca. Finally, a revised taxonomy of X. fastidiosa into three major clades defined by the subspecies pauca (clade I), multiplex (clade II) and the combination of fastidiosa, morus and sandyi (clade III) was strongly supported by k-mers specifically associated with these subspecies.

CONCLUSIONS

SkIf is a robust and rapid software, freely available, that can be dedicated to the comparison of sequence datasets and is applicable to any field of research. Applied to X. fastidiosa, an emerging pathogen in Europe, it provided an important resource to mine for identifying genetic markers of subspecies to optimize the strategies attempted to limit the pathogen dissemination in novel areas.

An Experimental Test of the Host-Plant Range of Nonrecombinant Strains of North American Xylella fastidiosa subsp. multiplex

Nonrecombinant strains of Xylella fastidiosa subsp. multiplex (those lacking evidence of significant intersubspecific homologous recombination) infect the xylem of a wide range of native and nonnative trees in North America. However, the degree to which different strains have a specialized host range remains poorly understood. We tested eight strains isolated from five different tree species (almond, olive, sweetgum, and plum in California and oak in Washington, DC). Experiments were conducted in greenhouses in Riverside, CA, and each strain was tested on 11 to 15 of the 17 plant species tested. Hosts infected by the most strains were plum (5 of 8 strains) and almond (4 of 8), while their congener peach was only infected by 1 of 8. No strains infected oleander or mulberry. All strains successfully infected their original host, with peach, olive (1 of 7), and sweetgum (2 of 6) only infected by such strains. Of the 90 total strain-novel-host combinations tested, 11 resulted in unambiguous infection, 2 gave ambiguous results, and the remaining 77 failed to result in symptoms or bacterial spread. All eight strains had a unique host range, including two pairs of strains with the same multilocus sequence typing sequence type, providing strong evidence of extensive plant-host specialization. There was little evidence that host relatedness was driving host specificity.

Novel amplification targets for rapid detection and differentiation of Xylella fastidiosa subspecies fastidiosa and multiplex in plant and insect tissues

Xylella fastidiosa is an insect-transmitted bacterial plant pathogen which causes a variety of economically important diseases worldwide. Molecular identification of X. fastidiosa is used for quarantine screening, surveillance, and research applications; many of which require subspecies level differentiation of pathogen isolates. This study describes quantitative PCR (qPCR) and isothermal amplification assays which can rapidly identify X. fastidiosa isolates belonging to the fastidiosa and multiplex subspecies. The TaqMan qPCR primers described here are used to differentiate X. fastidiosa strains by subspecies in plant and insect tissue in a single reaction, with the inclusion of a general amplification control probe to identify potential false negative samples. This TaqMan qPCR protocol can identify between 10³ and 10⁴ cfu/ml concentrations of X. fastidiosa at the subspecies level in a variety of sample types. Additionally, loop-mediated isothermal amplification (LAMP) targets were designed for faster detection of X. fastidiosa subspecies fastidiosa and multiplex, applicable to a field setting. These assays are effective for strain differentiation in artificially and naturally inoculated plant material, and in field collected insect vectors.

Using insects to detect, monitor and predict the distribution of Xylella fastidiosa: a case study in Corsica

We sampled ca 2500 specimens of Philaenus spumarius (Hemiptera: Aphrophoridae) throughout Corsica without a priori knowledge on the presence of symptoms on plants. We screened 448 specimens for the presence of Xylella fastidiosa (Xf) using qPCR and a custom nested PCR. qPCR appeared versatile and under-estimated the prevalence of Xf. Nested PCR showed that Xf was present in all populations. Molecular results were validated by prediction on the distribution of Xf made from tests conducted on plants, which shows the pertinence of using vectors in risk assessment studies. Xf was detected in tenerals and adults. Thus, P. spumarius could acquire Xf from its host plant, mostly Cistus monspeliensis in Corsica, which may act as reservoir for the next season. This contrasts with other observations and suggests that management strategies may have to be adapted on a case-by-case basis. At least two genetic entities and several variants of Xf not yet identified on plants were present in the insects, which suggests ancient introductions of Xf and a probable underestimation of the current diversity of the strains present in Corsica. Interestingly 6% of the specimens carried two subspecies of Xf. Studies are required to better characterize the strains present in Corsica and to determine how the disease was introduced, spread and why no sign of a potential epidemic was detected earlier. This study shows that, when sensitive enough methods are implemented, spittlebugs (and more specifically P. spumarius for which species distribution modelling shows it could be a good sentinel for Europe) can be used to predict and better assess the exact distribution of Xf. Furthermore, Xf multiply only in their foregut and does not become circulative, which facilitates its detection.

Updated pest categorisation of Xylella fastidiosa

Following a request from the European Commission, the EFSA Plant Health Panel updated its pest categorisation of Xylella fastidiosa, previously delivered as part of the pest risk assessment published in 2015. X. fastidiosa is a Gram‐negative bacterium, responsible for various plant diseases, including Pierce's disease, phony peach disease, citrus variegated chlorosis, olive quick decline syndrome, almond leaf scorch and various other leaf scorch diseases. The pathogen is endemic in the Americas and is present in Iran. In the EU, it is reported in southern Apulia in Italy, on the island of Corsica and in the Provence‐Alpes‐Côte d'Azur region in France, as well as in the Autonomous region of Madrid, the province of Alicante and the Balearic Islands in Spain. The reported status is ‘transient, under eradication’, except for the Balearic Islands, Corsica and southern of Apulia, where the status is ‘present with a restricted distribution, under containment’. The pathogen is regulated under Council Directive 2000/29/EC and through emergency measures under http://eur-lex.europa.eu/legal-content/EN/TXT/?uri=CELEX:32015D0789 (as amended http://eur-lex.europa.eu/legal-content/EN/TXT/?uri=CELEX:32017D2352). The pest could enter the EU via host plants for planting and via infectious insect vectors. The host range includes hundreds of host species listed in the EFSA host plant database. In the EU, host plants are widely distributed and climatic conditions are favourable for its establishment. X. fastidiosa can spread by movement of host plants for planting and infectious insect vectors. X. fastidiosa is known to cause severe direct damage to major crops including almonds, citrus, grapevines, olives, stone fruits and also forest trees, landscape and ornamental trees, with high impacts. The criteria assessed by the Panel for consideration as a potential Union quarantine pest are met (the pathogen is present in the EU, but it has a restricted distribution and is under official control). X. fastidiosa is not considered as a regulated non‐quarantine pest (RNQP) as the pathogen may spread also via insect vector transmission.

Ionome changes in Xylella fastidiosa-infected Nicotiana tabacum correlate with virulence and discriminate between subspecies of bacterial isolates

Characterization of ionomes has been used to uncover the basis of nutrient utilization and environmental adaptation of plants. Here, ionomic profiles were used to understand the phenotypic response of a plant to infection by genetically diverse isolates of Xylella fastidiosa, a gram-negative, xylem-limited bacterial plant pathogen. In this study, X. fastidiosa isolates were used to infect a common model host (Nicotiana tabacum 'SR1'), and leaf and sap concentrations of eleven elements together with plant colonization and symptoms were assessed. Multivariate statistical analysis revealed that changes in the ionome were significantly correlated with symptom severity and bacterial populations in host petioles. Moreover, plant ionome modification by infection could be used to differentiate the X. fastidiosa subspecies with which the plant was infected. This report establishes host ionome modification as a phenotypic response to infection.

Intersubspecific recombination in Xylella fastidiosa Strains native to the United States: infection of novel hosts associated with an unsuccessful invasion

The bacterial pathogen Xylella fastidiosa infects xylem and causes disease in many plant species in the Americas. Different subspecies of this bacterium and different genotypes within subspecies infect different plant hosts, but the genetics of host adaptation are unknown. Here we examined the hypothesis that the introduction of novel genetic variation via intersubspecific homologous recombination (IHR) facilitates host shifts. We investigated IHR in 33 X. fastidiosa subsp. multiplex isolates previously identified as recombinant based on 8 loci (7 multilocus sequence typing [MLST] loci plus 1 locus). We found significant evidence of introgression from X. fastidiosa subsp. fastidiosa in 4 of the loci and, using published data, evidence of IHR in 6 of 9 additional loci. Our data showed that IHR regions in 2 of the 4 loci were inconsistent (12 mismatches) with X. fastidiosa subsp. fastidiosa alleles found in the United States but consistent with alleles from Central America. The other two loci were consistent with alleles from both regions. We propose that the recombinant forms all originated via genomewide recombination of one X. fastidiosa subsp. multiplex ancestor with one X. fastidiosa subsp. fastidiosa donor from Central America that was introduced into the United States but subsequently disappeared. Using all of the available data, 5 plant hosts of the recombinant types were identified, 3 of which also supported non-IHR X. fastidiosa subsp. multiplex, but 2 were unique to recombinant types from blueberry (7 isolates from Georgia, 3 from Florida); and blackberry (1 each from Florida and North Carolina), strongly supporting the hypothesis that IHR facilitated a host shift to blueberry and possibly blackberry.

Recent evolutionary radiation and host plant specialization in the Xylella fastidiosa subspecies native to the United States

The bacterial pathogen, Xylella fastidiosa, infects many plant species in the Americas, making it a good model for investigating the genetics of host adaptation. We used multilocus sequence typing (MLST) to identify isolates of the native U.S. subsp. multiplex that were largely unaffected by intersubspecific homologous recombination (IHR) and to investigate how their evolutionary history influences plant host specialization. We identified 110 "non-IHR" isolates, 2 minimally recombinant "intermediate" ones (including the subspecific type), and 31 with extensive IHR. The non-IHR and intermediate isolates defined 23 sequence types (STs) which we used to identify 22 plant hosts (73% trees) characteristic of the subspecies. Except for almond, subsp. multiplex showed no host overlap with the introduced subspecies (subspecies fastidiosa and sandyi). MLST sequences revealed that subsp. multiplex underwent recent radiation (<25% of subspecies age) which included only limited intrasubspecific recombination (ρ/θ = 0.02); only one isolated lineage (ST50 from ash) was older. A total of 20 of the STs grouped into three loose phylogenetic clusters distinguished by nonoverlapping hosts (excepting purple leaf plum): "almond," "peach," and "oak" types. These host differences were not geographical, since all three types also occurred in California. ST designation was a good indicator of host specialization. ST09, widespread in the southeastern United States, only infected oak species, and all peach isolates were ST10 (from California, Florida, and Georgia). Only ST23 had a broad host range. Hosts of related genotypes were sometimes related, but often host groupings crossed plant family or even order, suggesting that phylogenetically plastic features of hosts affect bacterial pathogenicity.

A conjugative 38 kB plasmid is present in multiple subspecies of Xylella fastidiosa

A ≈ 38kB plasmid (pXF-RIV5) was present in the Riv5 strain of Xylella fastidiosa subsp. multiplex isolated from ornamental plum in southern California. The complete nucleotide sequence of pXF-RIV5 is almost identical to that of pXFAS01 from X. fastidiosa subsp. fastidiosa strain M23; the two plasmids vary at only 6 nucleotide positions. BLAST searches and phylogenetic analyses indicate pXF-RIV5 and pXFAS01 share some similarity to chromosomal and plasmid (pXF51) sequences of X. fastidiosa subsp. pauca strain 9a5c and more distant similarity to plasmids from a wide variety of bacteria. Both pXF-RIV5 and pXFAS01 encode homologues of a complete Type IV secretion system involved in conjugation and DNA transfer among bacteria. Mating pair formation proteins (Trb) from Yersinia pseudotuberculosis IP31758 are the mostly closely related non-X. fastidiosa proteins to most of the Trb proteins encoded by pXF-RIV5 and pXFAS01. Unlike many bacterial conjugative plasmids, pXF-RIV5 and pXFAS01 do not carry homologues of known accessory modules that confer selective advantage on host bacteria. However, both plasmids encode seven hypothetical proteins of unknown function and possess a small transposon-associated region encoding a putative transposase and associated factor. Vegetative replication of pXF-RIV5 and pXFAS01 appears to be under control of RepA protein and both plasmids have an origin of DNA replication (oriV) similar to that of pRP4 and pR751 from Escherichia coli. In contrast, conjugative plasmids commonly encode TrfA and have an oriV similar to those found in IncP-1 incompatibility group plasmids. The presence of nearly identical plasmids in single strains from two distinct subspecies of X. fastidiosa is indicative of recent horizontal transfer, probably subsequent to the introduction of subspecies fastidiosa to the United States in the late 19(th) century.

Detecting genetic introgression: high levels of intersubspecific recombination found in Xylella fastidiosa in Brazil

Documenting the role of novel mutation versus homologous recombination in bacterial evolution, and especially in the invasion of new hosts, is central to understanding the long-term dynamics of pathogenic bacteria. We used multilocus sequence typing (MLST) to study this issue in Xylella fastidiosa subsp. pauca from Brazil, a bacterium causing citrus variegated chlorosis (CVC) and coffee leaf scorch (CLS). All 55 citrus isolates typed (plus one coffee isolate) defined three similar sequence types (STs) dominated by ST11 (85%), while the remaining 22 coffee isolates defined two STs, mainly ST16 (74%). This low level of variation masked unusually large allelic differences (>1% divergence with no intermediates) at five loci (leuA, petC, malF, cysG, and holC). We developed an introgression test to detect whether these large differences were due to introgression via homologous recombination from another X. fastidiosa subspecies. Using additional sequencing around these loci, we established that the seven randomly chosen MLST targets contained seven regions of introgression totaling 2,172 bp of 4,161 bp (52%), only 409 bp (10%) of which were detected by other recombination tests. This high level of introgression suggests the hypothesis that X. fastidiosa subsp. pauca became pathogenic on citrus and coffee (crops cultivated in Brazil for several hundred years) only recently after it gained genetic variation via intersubspecific recombination, facilitating a switch from native hosts. A candidate donor is the subspecies infecting plum in the region since 1935 (possibly X. fastidiosa subsp. multiplex). This hypothesis predicts that nonrecombinant native X. fastidiosa subsp. pauca (not yet isolated) does not cause disease in citrus or coffee.