From Nucleotides to Satellite Imagery: Approaches to Identify and Manage the Invasive Pathogen Xylella fastidiosa and Its Insect Vectors in Europe

Biological invasions represent some of the most severe threats to local communities and ecosystems. Among invasive species, the vector-borne pathogen Xylella fastidiosa is responsible for a wide variety of plant diseases and has profound environmental, social and economic impacts. Once restricted to the Americas, it has recently invaded Europe, where multiple dramatic outbreaks have highlighted critical challenges for its management. Here, we review the most recent advances on the identification, distribution and management of X. fastidiosa and its insect vectors in Europe through genetic and spatial ecology methodologies. We underline the most important theoretical and technological gaps that remain to be bridged. Challenges and future research directions are discussed in the light of improving our understanding of this invasive species, its vectors and host–pathogen interactions. We highlight the need of including different, complimentary outlooks in integrated frameworks to substantially improve our knowledge on invasive processes and optimize resources allocation. We provide an overview of genetic, spatial ecology and integrated approaches that will aid successful and sustainable management of one of the most dangerous threats to European agriculture and ecosystems.


Introduction
Biological invasions represent one of the most severe threats to local communities and ecosystems [1]. Invasions occur when species are intentionally or accidentally introduced from the native or historic range into a new area, spread in the novel environment successfully and adversely affect it [2]. The economic loss linked to this phenomenon is conspicuous: a potential estimate of the known losses due to alien species in Asia, Australia, Europe and North America is several billion dollars per year [3,4]. Even more prominent are the ecological and social impacts of alien organisms on native species and ecosystems, which include reduced biodiversity, decreased availability and quality of key natural resources, increased frequency of wildfires, flooding and pollution [5].
Invasive plants and animals have been a traditional focus of research and management efforts [2,6]; however, invasive microorganisms also represent a serious threat, as exemplified by the spread of alien pathogens inducing human, animals or plants diseases [7,8]. Microbial invasions in terrestrial ecosystems are spreading exceptionally rapidly due to unprecedented migrations of organisms following anthropogenic activities and climate change [7,8]. Several of them are pathogenic and impact communities and/or environments adversely [7][8][9], as recently highlighted by several epidemic coronavirus [10]. Pathogens linked to emerging infective diseases of plants harm primary economic resources [11]. One of the most dangerous plant pathogens is the bacterium Xylella fastidiosa (hereafter Xf) [12], which is responsible for a wide variety of diseases and has profound environmental, agricultural and economic impacts [13,14].
Xf is a Gram-negative bacterium that colonizes the plant xylem. Taxonomically, it is a single species but it is further categorized into subspecies and strains differing in genetic diversity, geographic distribution and host specialization ( [15], details in Section 2.2). Xf is the causal agent of many severe diseases in a wide range of wild plants and economically relevant agricultural crops; additionally, several other plants host this pathogen asymptomatically (ca. 350 plant host species [16]). Infected plants greatly vary in disease manifestation, timing of appearance and severity; a crucial role is played by host-pathogen interactions, whose mechanisms are still unclear [17]. In asymptomatic plant hosts, Xf lives in free-ranging small colonies in xylem vessels. Pathogenicity occurs when the bacterium rapidly multiplies locally forming biofilm clusters and compromises water and nutrients transport, leading to plant quick decline [17,18].
Once restricted to the Americas, Xf reached Europe in 2013, infecting the Apulian olive trees (southern Italy [19]). Later, the pest was identified in several plant species in France in 2015, Spain in 2016, Tuscany (central Italy) in 2018 and Portugal in 2019 [13]. Xf was also isolated in infected plants imported to Germany and the Netherlands from South America but it was quickly eradicated [15]. Following the European outbreak, Xf was initially assigned to the quarantine status in the European and Mediterranean Plant Protection Organization (EPPO) A1 list (pest absent in the EPPO region) and then transferred to the A2 List (locally present in the EPPO region) in September 2017 [20].
The long-range introduction of Xf in Europe, likely originating from Costa Rica, is linked to anthropogenic activities, particularly commercial trade of asymptomatic infected plant material (discussed in Section 2.2, Figure 1 [21]). However, Xf is a vector-borne bacterium that colonizes both plant xylem and insect gut: the small-scale spread of this bacterium among local plants is hence facilitated by xylem-sap feeding insects. Infection occurs when an insect feeds on an infected plant and subsequently on a healthy However, Xf is a vector-borne bacterium that colonizes both plant xylem and insect gut: the small-scale spread of this bacterium among local plants is hence facilitated by xylem-sap feeding insects. Infection occurs when an insect feeds on an infected plant and subsequently on a healthy plant [22]. Insects that ingest primarily xylem sap belong to four families within the plant [22]. Insects that ingest primarily xylem sap belong to four families within the Hemiptera: Cicadellidae (subfamily Cicadellinae, sharpshooter leafhoppers), Cercopidae (spittlebugs), Machaerotidae (tube-building spittlebugs) and Cicadidae (cicadas) [23]. In North and South America, primary insect vectors all belong to the Hemiptera order and include the Proconiini and Cicadellini insects and spittlebugs (Cercopidae) but other xylem-feeding bugs are also suspected to transmit the bacterium with the exception of cicadas [24][25][26]. In Europe, spittlebugs are the dominant group of potential Xf vectors [25,26]. They have been long known to transmit Xf, but only few studies have specifically addressed their active role in pathogen spread [27][28][29][30].
In the Apulian olive orchards, the Xf insect host is the meadow spittlebug Philaenus spumarius (Linnaeus, 1758) (Hemiptera: Aphrophoridae) [30][31][32], whose success in spreading Xf is based on its diet including a wide array of plant species, a long-life cycle and its abundance and invasiveness. Recently, other two species, Philaenus italosignus (Drosopoulos and Remane, 2000) (Hemiptera: Aphrophoridae) and Neophilaenus campestris (Fallén, 1805) (Hemiptera: Aphrophoridae), have also been confirmed to be vectors of Xf subsp. pauca ST53 to olive plants in Italy but only under experimental conditions [33]. Other co-occurring Hemiptera, such as Euscelis lineolatus (Brullé, 1832) (Hemiptera: Cicadellidae) and Cicada orni (Linnaeus, 1758) (Hemiptera: Cicadidae), are suspected to contribute to the spreading of this bacterium but genetic screenings (see Section 2.3) were negative and their impact is potentially more limited due to their restricted host range and seasonal or geographic constraints [22,23,25,34]. In the light of this unclear information, the current recommendation is to consider all spittlebugs as potential vectors until proven otherwise [22,23,34,35]. Clarifying the agents and routes of infection urgently needs further investigation [36].
No effective treatment to eradicate Xf is available thus far. Some chemical and biological control measures may temporarily reduce disease severity, as well as through insect vector control, but they may not be successful in eradicating Xf completely [37,38]. Therefore, prevention and containment are the most appropriate strategies to minimize the impact of Xf outbreaks in Europe. These policies are primarily based on accurate surveillance, followed by the eradication of potential Xf sources and the control of its insect vectors [39]. All these measures require a detailed knowledge of the biology, ecology and spatial distribution of Xf and its hosts.
Here, we aim to provide practical toolkit to scientists, stakeholders, policymakers and citizens to efficiently deal with such a threatening pathogen. We review the most recent advances on the identification, distribution and management of Xf and its insect vectors in Europe through genetic and spatial ecology methodologies ( Figure 2). Many studies have clarified the biology of Xf and its insect vectors but advances to successfully control them possibly without invasive procedures have been limited [16,23]. Latest developments Many studies have clarified the biology of Xf and its insect vectors but advances to successfully control them possibly without invasive procedures have been limited [16,23]. Latest developments in genetics and spatial ecology have proven to be markedly useful in the identification and management of invasive species, particularly through an integrated approach [40,41], which is currently largely missing Sustainability 2020, 12, 4508 4 of 38 in the Xf research. It is clear that a detailed and inclusive insight into ecological and evolutionary dynamics at multiple levels of biodiversity is crucial to identify, predict and prevent the spread of threatening organisms [42][43][44]. We discuss how genomic-scale methods would enhance ecological studies (and vice versa) of this and other invasive species and highlight promising directions worth exploring. To our knowledge, this is the first study to integrate genetic as well as ecological concepts and methodologies for this invasive bacterium. We conclude by providing indications on how to gain a comprehensive understanding of Xf, plants, insect vectors and host-pathogen interactions and aid a successful and sustainable management of one of the most dangerous threats to the European agriculture.

Genetic Approaches
Genetic applications are able to illuminate diversity, demography, ecology, adaptive potential and impact of invasive species; this essential information may not be accessible through other methods [45]. Despite the increase of theoretical and empirical research on invasive species in the last years, the genetic basis and evolutionary implications of biological invasions remain poorly understood [46][47][48]. The advent of the "omics" era is facilitating the application of genetic approaches for conservation and management purposes, for example by enabling a rapid discovery and taxonomic classification of alien species or tracking invasion routes [49][50][51]. Additionally, they provide powerful tools for the control of invasive species [52][53][54][55] and could help manage Xf and its insect vectors, as discussed in the following sections.

Sequencing Data
Due to the economic and ecological importance of Xf, draft and complete genomic resources have been available since the species' proper classification in 1987 [12]. As to March 2020, 31 assemblies (Box 1) for the European strains of this species have been released in the National Centre for Biotechnology Information (NCBI) genome database (Table 1 and Table S1). Additionally, a comparative genomic database with manually curated annotations of six strains is available ( Table 1): none of these six genomes is European but they provide a valuable source for comparative studies. The first and representative genome, Xf subspecies pauca 9a5c (citrus-specific phytopathogen from Brazil), has a total length of 2.73 Mb, includes a single circular genome and two plasmids, and has 2330 identified protein-coding genes but half of them have an unknown function ( [56], NCBI, Table 1). Subspecies, including the European strains, differ in genomic structure, content and diversity (Table 1, NCBI, Table S1). Only 900 genes are shared among the Xf subspecies pauca, fastidiosa and multiplex (core genomes, Box 1) while their accessory genomes (Box 1) are conspicuous (ca. 1500 genes, Table S1 [57,58]).
Xf was the first plant pathogen to enter in the genomic era [56]. The amount of existing genomic information has markedly increased since 2013 due to the rising concern following the outbreaks in Europe and the decreasing costs of sequencing technologies (NCBI , Table 1). However, while the quality and completeness of Xf assemblies has improved over the years, most of them are still Sustainability 2020, 12, 4508 5 of 38 fragmented, likely due to the presence of repeated regions [59]. In Europe, only three isolates (Box 1) of a single Xf strain have fully assembled and circularized genomes; the great majority of available genomic information are represented by fragmented assemblies (Table S1). Box 1. Definitions of some genetic terms used in the main text.

Term
Definition Accessory genome Part of genes that are not present in all strains of a species; typically includes strain specific adaptation such as antibiotic resistance.
Clade A group of organisms that include a common ancestor and all its descendants.

Competence
Ability of a microorganism to incorporate exogenous DNA in its genome.
Core genome Part of genes that are present in all strains of a species; typically includes housekeeping genes for cell envelope or regulatory functions.

Genetic vector
Vehicles (e.g., DNA plasmids, viruses, artificial chromosomes) for delivering foreign DNA into recipient cells.
Genome assembly A computational representation of a genome sequence. Genomes are not sequenced along the complete length of a chromosome but fragmented: these fragments are sequenced and then put back together (assembled). Xf whole genome sequences were initially obtained using random shotgun cloning strategies and Sanger sequencing [56,60,61]. Later, studies took advantage of the cheaper and more efficient massively parallel short reads sequencing approaches such as Illumina platforms [62][63][64][65][66][67]. In some studies, short (Illumina) reads strategies have been complemented by long reads (PacBio, Nanopore) and single molecule real time sequencing to reduce assembling complexity ( [59,64,66], Table S1, see [68,69] for detailed reviews of sequencing technologies). A combination of long (conserved genome organization, genome closure), short (high coverage and high quality polishing) and Sanger (gap bridging) sequencing, which optimizes the advantages vs. drawbacks of these methods is typically recommended [69].

Epigenetic variation
Annotation, the identification and localization of protein-coding genes and other functional regions in the genome, has been based on homology-based search tools and/or ab initio gene prediction algorithms (listed in Table S1). They all are automatic annotation pipelines, which can produce inaccurate results due to methodological inconsistencies, spelling mistakes, pseudogenes, hypothetical proteins and paralogs [70]. This likely contributed to the incomplete genomic representation of Xf ; additional studies and functional assays are required for a more accurate picture on gene content and its diversity within and among subspecies [71].
Vector genomes can illuminate the genetic basis of pathogen-vector interactions, their spread and invasion routes. The genome of the meadow spittlebug P. spumarius, the primary vector of Xf in Europe, has been recently sequenced through 10× genomics linked-reads [72], a sequencing technology based on short reads that preserves long-range information. The resulting genome assembly of 2.7 Gb is highly fragmented, likely due to the high heterozygosity; nonetheless, it is considered complete and representative of P. spumarius' gene content. Efforts to generate a chromosome-level assembly for this and additional ten species of xylem-feeding insects are underway (Table S2, Roberto Biello, unpublished work).

Genetic Diversity
A deep understanding of genetic relationships and diversity within Xf is of paramount importance for the management and prevention of outbreaks. In fact, while this pathogen is known to infect a wide range of plant species, its strains (Box 1) show specialization for a limited number of hosts. Some host-pathogen genotype associations remain asymptomatic and the emergence of new diseases and host ranges is associated with genomic recombination between different strains/subspecies [73].
Various approaches have been used to explore genetic structure in Xf : single-or multi-locus polymorphic genetic markers and, more recently, whole genome sequencing. Genomic analyses are being facilitated by the decreasing costs of high-throughput sequencing technologies and bioinformatic tools to quickly extract information associated to a group of organisms from available genomes and efficiently explore diversity in the large Xf pangenome (Box 1 [58,[74][75][76], Table S1, see [77] and Section 2.3 for further methodological discussion). Despite the copious efforts, a conclusive classification of Xf subspecies is still debated mainly due to methodological differences among studies, biased sampling limited to important crop plants and horizontal gene transfers that are favored by the natural competence of this organism [73]. According to the EPPO, three subspecies are formally recognized: Xf subspecies fastidiosa (main plant hosts: grapevine, almond), multiplex (stone fruits, share trees, olive, plums) and pauca (citrus, coffee, olive) [15]. Among them, only the first two are accepted by the International Society of Plant Pathology Committee on the Taxonomy of Plant Pathogenic Bacteria. Multi locus sequence typing (MLST, Box 1, discussed below), phylogenetic, genomic and ecological studies have indicated further genetic and biological diversity and the presence of other potential bacterial subspecies such as Xf subspecies sandyi (main plant hosts: oleander, magnolia), tashke (chitalpa) and morus (mulberry). Their relationships are not fully resolved and the same host species can be infected by various pathogen subspecies [15,31,74,76,[78][79][80][81][82][83][84][85]. In Europe, four Xf subspecies have been reported since the first outbreak in 2013: fastidiosa, multiplex, pauca and sandyi ( Figure 2 [15,31,39,80,[86][87][88]). A recent phylogenomic-calibrated tree estimated a divergence time Sustainability 2020, 12, 4508 7 of 38 among subspecies between 10,000 (Xf subspecies pauca vs. Xf subspecies fastidiosa plus Xf subspecies multiplex) and 3400 (Xf subspecies fastidiosa vs. Xf subspecies sandyi) years before present [58]. Each subspecies is under different selective pressures, contributing to their divergence [58,89].
Within subspecies, additional grouping of genotypes with distinct genetic and biological features is based on MLST, a largely accepted portable and robust genetic typing method to standardize genetic analyses [90]. As March 2020, 87 ST profiles have been reported worldwide, twelve from Europe (Xf MLST database, Table 1 and Figure 3 [63,76,78,82,87,[90][91][92][93][94][95]). Importantly, new STs have been described for Europe since 2013 representing hitherto undescribed genotypes, revealing that our current knowledge of this pathogen is still incomplete and/or Xf is rapidly evolving [96].
Although MLST revolutionized taxonomy in the face of genetic recombination, it is based on conserved genomic regions which may not have full discriminatory power. In fact, studies based on more polymorphic markers and genomic data showed a lack of ST monophyly and further genetic differentiation among subspecies and strains [66,[97][98][99]. Differentiation is sometimes reflected in genomic rearrangements, such as the inversion that distinguishes the Apulian isolates from the Brazilian ones of the same subspecies Xf pauca [66]. Divergence among subspecies/strains/isolates is in some cases pronounced enough to question their affiliation to the same group or even species [66,[97][98][99].
Sustainability 2020, 12, x FOR PEER REVIEW 7 of 38 revealing that our current knowledge of this pathogen is still incomplete and/or Xf is rapidly evolving [96].
Although MLST revolutionized taxonomy in the face of genetic recombination, it is based on conserved genomic regions which may not have full discriminatory power. In fact, studies based on more polymorphic markers and genomic data showed a lack of ST monophyly and further genetic differentiation among subspecies and strains [66,[97][98][99]. Differentiation is sometimes reflected in genomic rearrangements, such as the inversion that distinguishes the Apulian isolates from the Brazilian ones of the same subspecies Xf pauca [66]. Divergence among subspecies/strains/isolates is in some cases pronounced enough to question their affiliation to the same group or even species [66,[97][98][99]. and/or intercepted in Europe (data from Table S1 [31,95]). Xf was intercepted in Switzerland but information on the subspecies or ST is not available. NA, Not Available.
Whole genome sequencing has provided a significantly finer resolution than other approaches, which is necessary to discriminate between different microevolutionary scenarios and routes of dispersal. The four subspecies found in Europe have originally evolved in America in allopatry: Xf subspecies pauca in South America, Xf subspecies multiplex in North America, Xf subspecies fastidiosa in Southern Central America and Xf subspecies sandyi in Southern USA. Then, they have dispersed repeatedly in North, Central and South America, from which they reached Europe via Costa Rica (Xf  Europe (data from Table S1 [31,95]). Xf was intercepted in Switzerland but information on the subspecies or ST is not available. NA, Not Available.
Whole genome sequencing has provided a significantly finer resolution than other approaches, which is necessary to discriminate between different microevolutionary scenarios and routes of dispersal. The four subspecies found in Europe have originally evolved in America in allopatry: Xf Sustainability 2020, 12, 4508 8 of 38 subspecies pauca in South America, Xf subspecies multiplex in North America, Xf subspecies fastidiosa in Southern Central America and Xf subspecies sandyi in Southern USA. Then, they have dispersed repeatedly in North, Central and South America, from which they reached Europe via Costa Rica (Xf subspecies pauca and sandyi), California (Xf subspecies sandyi and fastidiosa) and the Southeastern USA (Xf subspecies multiplex; Figure 1 [58,62,64,66,76,82,90,94,[100][101][102][103][104][105]). The Xf subspecies multiplex has been introduced to Corsica multiple times but only once to Italy and Spain. The Italian Xf subspecies pauca infection is the result of a founder event by a single or few closely related lineages in 2008 [58,66,106]. Anthropogenic activities, particularly commercial trade of asymptomatic infected plant material, are the main vector of introduction of Xf in Europe ( Figure 1 [80]).
Xf exhibits high genomic variability and plasticity, allowing for a high adaptive potential [58,59,104,105,[107][108][109]. Host adaptation in this pathogen is mainly facilitated by multiple inter-and intra-subspecific homologous recombination events that created new genotypes able to invade new ecological niches [58,59]. Wild-type and experimentally generated hybrids show that transformation frequency is variable among strains/subspecies and increases with the concentration of environmental calcium, the number of sympatric strains and the lineage age [21,58,59,108,110]. Horizontal gene transfer is well described in the American strains. In Europe, the Xf subspecies pauca ST53 seems to have acquired genes from other subspecies or different species before its introduction in this continent, while evidence of post-colonization recombination is not clear in the Xf subspecies multiplex and fastidiosa. Highly hybrid strains intercepted in Europe from plants imported from Central and South America include Xf subspecies pauca CFBP8072, Xf subspecies sandyi XFC033 and Xf subspecies fastidiosa CFBP8073 [21,58,59,102]. Interestingly, all subspecies experienced recent intersubspecific recombination at common loci that are important to colonize the host environment [59] and some European strains show recombination with uncharacterized genotypes (e.g., CFPB8416 [21]).
Culture purification, genetic homology and circular contigs revealed the presence of numerous mobile elements such as plasmids and prophages (Box 1) in the Xf genome (18% in the strain 9a5c [111]). These genetic elements are differentially enriched among Xf strains, encode accessory modules conferring selective advantages in specific environments, may activate/inactivate host genes and facilitate horizontal gene transfers [57,112]. Although additional plasmids in Xf continue to be characterized, they are difficult to analyze due to their low resolution and repeated sequences (Box 1 [74,[113][114][115][116]). Additionally, their role during infections is still largely unknown; in fact, only few studies have moved beyond inference from homology with integrated in vivo and in vitro analyses [57,112].
Few studies have explored genetic diversity in the European Xf insect vectors so far and those were limited to few loci to clarify the insect's colonization history or the genetic basis of color polymorphism [117][118][119]. None of these studies investigated genetic structure of these insects in the light of their role as potential Xf vectors, yet it may underlie relevant biological and epidemiological differences. The upcoming availability of numerous genomic resources (see Section 2.1) will facilitate comparative studies and illuminate the genomic basis of Xf infection in insect hosts.

Diagnostic Protocols
No treatment or control measure is currently effective to eradicate Xf [95]. Therefore, fast and reliable diagnostic tools for an early and accurate discovery and identification of this pathogen are crucial for successful management of current and future epidemics in Europe. Phenotypic approaches or plant symptoms have proven to be inefficient monitoring and survey tools as early stages of infection can be asymptomatic and isolation or microscopy analyses is difficult and time-consuming [31,120]. A panel of standardized diagnostic methods to monitor the presence of Xf infections in Europe based on serological or genetic tests is available and annually revised [15].
Serological methods based on antibodies and immunofluorescence were the first to be established but are more expensive, time-consuming and less informative than nucleotides-based approaches [15]. In particular, they are not effective in determining subspecies [15,39,77]. Only a limited number Sustainability 2020, 12, 4508 9 of 38 of studies applied serological protocols to the European strains, often in combination with genetic strategies [93,[121][122][123][124][125].
Genetic methods have been the most widely used for a rapid detection of Xf so far. They include amplification of specific Xf genes through classic or more advanced polymerase chain reaction (PCR) approaches, single nucleotide polymorphisms (SNPs) and MLST [39,77,89,92,126]. These tests greatly vary in specificity, sensitivity, accuracy, flexibility, portability to the field and costs in terms of time, skilled labor and money; the choice of the most effective diagnostic technique largely depends on aims (e.g., quick pathogen detection vs. precise subspecies identification) and hosts [15,39,77,[126][127][128][129]. MLST is a classification and diagnostic scheme that is strongly recommended in the case of new outbreaks but is not effective for large-scale screenings and to discriminate among closely-related lineages [63,76,82,87,[90][91][92][93][94]98]. To this aim, SNPs provides deeper resolution, but their discovery is heavily based on an appropriate reference genome(s) and results from different studies might not be comparable [68]. Developing faster, more reliable and efficient genetic tools that allow the simultaneous detection of Xf and its identification at the subspecies/strain level as well as the detection of mixed infections of multiple subspecies in one sample has been the focus of the most recent assays [91,[129][130][131]. Other efforts have aimed to portability, low-cost and ease of use with the development of miniaturized devices such as lab-on-chip, which shows intermediate performance between ELISA and RT-PCR methods [132].
The genetic methods used to diagnose Xf are based on a limited set of loci, require prior knowledge of the analyzed sequences and are prone to false-positive and false-negative results [77,128]. Importantly, they are mainly restricted to regions of the genome that are present in all the analyzed isolates (most likely the core genome), potentially discarding useful information from the accessory genome [116]. Multiple approaches are often required to increase diagnostic power and resolution [77,128]. Genome approaches are not yet time and cost effective for the rapid identification of pathogens and for field applications, but they allow for the quick discovery of multiple, genome-wide makers (e.g., SNPs) for strain/isolate classification and for tracing and monitoring projects [77,129]. Recently, metagenomic approaches allowed detecting and characterizing Xf without the need for pathogen cultivation and amplification of specific bacterial genes, but subspecies or strain classification was limited to material with high bacterial concentrations [129][130][131][132].
A recurrent issue with the above-mentioned approaches is the low concentration or uneven distribution of bacteria in host tissues, as well as the presence of metabolites that impair amplification protocols [128,[131][132][133]. Several approaches were developed to reduce these challenges [78,91,130,[133][134][135][136], but they require expensive infrastructure and skilled operators. A fast and cheap alternative approach is targeting plant physiological markers involved in the early, pre-symptomatic host response to the pathogen in advanced PCRs; this method can be particularly useful in the field due to its portability [127,137,138].
Genetic diagnostic protocols to identify Xf within insect vectors largely overlap with the methods developed for plants host discussed above, including MLST for subspecies classification, with some modifications to improve isolation from insects [86]. Xf multiplies in the insect foregut and does not spread to other organs. This trait facilitates Xf isolation within vectors [86], although bacteria may exist at low concentrations and PCR inhibitors challenge Xf identification [13,139,140].
The classification of Xf infected insects is typically based on the morphological analysis of adult specimens, but genetic tools applicable to all life stages greatly facilitate the task. Standard species identification through cytochrome c oxidase I (COI) are routinely used ( [141], Consortium for the Barcode of Life CBOL, Table 1). Higher discriminatory power could be reached by sequencing the region at the 3 end of the COI, as well as the 12S and 16S genes, although there is still no evidence that these regions will allow the discrimination of cryptic insect species [142,143]. Recent advances in the mitochondrial genome sequencing [144][145][146], long-range PCR [147], nuclear whole genome sequencing and SNPs panels will dramatically improve pests screening in the near future.

Gene Expression Profiles
Gene expression is a crucial link between genotype and phenotype; its regulation plays a key role in adaptation and evolution. Transcriptomic studies in Xf and plants under different conditions contributed to shed light on this pathogen and the processes underlying host-bacteria interaction and disease emergence.
Gene expression profiling of Xf, together with genomic comparisons and mutagenesis (see Section 2.5), clarified the genetic basis and mechanisms of colonization and pathogenicity, although several processes remain elusive largely due to uncharacterized or nonexclusive genetic elements [17,56,106]. Hybridization-based microarray, quantitative PCR (qPCR) and RNAseq studies revealed that differences in genomic sequences between strains are reflected and amplified in gene expression, with a conserved core gene pool linked to survival in plant xylem and a flexible set including genes that favors adaptation in different hosts and putative pathogenicity factors: Xf plastically adapts its survival strategy to the environment [111,[148][149][150][151][152]. They also showed that plasmids and prophages are transcriptionally active and modulated independently from the rest of the genome [111,149].
Transcription control is environment-dependent and it should therefore be analyzed in vivo or in conditions that closely resemble natural settings to meaningfully represent the trait under analysis. As in planta studies are technically challenging in Xf, most analyses of gene expression in this bacterium have been performed in vitro in batch cultures, but this pathogen survives only in environments under continuous flow (xylem vessels and insects' feeding canal [17,110]). Whole transcriptome sequencing and qPCR in microfluidic chambers, which mimics the natural habitats of this bacterium, showed that gene expression and the regulatory role of environmental elements such as calcium are different between batch and flow settings [110]. This difference should be considered for planning new transcriptome studies and for correctly interpreting previous results.
Differently from other bacterial pathogens, Xf lacks virulence genes and effector proteins [60]; disease induction largely depends on host-specific factors [111,153]. Following infection, major plant pathways are activated including genes linked to immunity, bacteria and embolism sensing, calcium metabolism and lignification [125,[154][155][156]. The pathogen is perceived as an abiotic stress related to drought [110]; however, plant gene expression changes associated with Xf are clearly distinct from those associated to water stress [157,158]. Comparative transcriptome and physiological analyses in resistant and susceptible plants highlighted that resistance to this pathogen is linked to an advanced plant ability to quickly and effectively respond to water deficiency [138,[156][157][158][159]. Additionally, some resistant plants may control Xf without showing dramatic alterations in genes expression of the defense cascade. In such cases, further uncharacterized genetic elements or other mechanisms such as a lower anatomical cavitation susceptibility or an ionomic composition that is protective against disease might be involved [159,160].
Only a handful of the above-mentioned studies used high throughput sequencing approaches; among these, just few involved the European species and were restricted to host plants [161,162]. The application of RNAseq and its most recent developments [163] has been hampered by the cost of these sequencing technologies (now decreasing), the dilution issue when total RNA from both plant and pathogen is sequenced, the presence of repetitive regions and the low quality of most available genomes [137,152]. However, these high-throughput tools provide a more detailed and quantitative view of gene expression, alternative splicing and allele-specific expression compared to former techniques (such as microarrays), without the need of a priori knowledge on sequences [164]. It could also help uncovering gene regulatory networks, another unexplored avenue in Xf.
Gene expression profiles have been poorly analyzed in Xf insect vectors. Several studies have uncovered the genetic basis of phenotypic plasticity and adaptation of insect pests to their host plants [165][166][167][168] but none of them has focused on spittlebugs. In other herbivorous pests, transcriptomic data in caterpillars feeding on multiple plants showed two sets of highly expressed genes involved in the evolution of polyphagia [166]. In the green peach aphid Myzus persicae (Sulzer, 1776) (Hemiptera: Aphididae), the differential regulation of genes belonging to the aphid-expanded gene families underlies hosts plasticity [165]. Future studies on gene expression profiles in P. spumarius and other potential Xf insect vectors might help explain its wide ecological niche and its ability to transmit Xf to multiple plant species, with significant consequences on the development of management strategies.

Genetic Engineering
Exploratory genetic and genomic tools, such as the ones described in the previous sections, have significantly improved our understanding of Xf pathogenicity. However, complementary methods are needed to validate these findings, provide insights into the functional significance of observed variation, especially in the light of the fragmented annotation of this pathogen. Genetic toolkits such as target alteration of gene expression through mutagenesis and site-specific genetic engineering have rapidly expanded in the last years and have been widely used to study gene functions and for trait improvements [169].
Xf is naturally competent (Box 1 [108]) but this feature does not facilitate genetic editing of this organism [170][171][172]. The application of conventional techniques to modify expression patterns of specific bacterial genetic elements using transposon mutagenesis and homologous recombination (methods in [173]) has been limited by laborious protocols and the slow-growing nature of this bacterium [12,120,171,[174][175][176][177][178][179]. Recent developments have optimized robustness, efficiency and time effectiveness of genetic manipulation protocols [171]. Successful studies involving comparisons of wild-type and overexpressed or knocked-out Xf mutants and complementation analyses have clarified the mechanisms underlying virulence or host specificity and helped characterizing the function of unknown genetic elements [171,175,177,[180][181][182][183][184]. Some transgenic organisms were created for practical applications such as disease suppression, including a genetically modified Xf avirulent strain [180,185].
Host plants too have been selectively bred or engineered with genetic vectors (Box 1) containing target genes to increase their resistance to Xf infection [186][187][188][189][190]. Genetic engineering techniques such as transgenic approaches provide a high precision, fast and robust alternative to the laborious, lengthy and often unsuccessful breeding methods [191,192]. In Xf, ectopic expression of pathogen genes related to signal transduction or anti-apoptotic process in transgenic plants through indirect transformation techniques successfully reduced bacterial colonization and disease severity ("pathogen confusion", Box 1 [188][189][190]; methods reviewed in [191]). A significant decrease of disease symptoms has been obtained also by a transgene expressing an antimicrobial gene or a protein chimera with recognition and clearance domains that were specifically conveyed at the site of colonization [186,187].
Genetic engineering provides an opportunity to control Xf insect vectors efficiently. For example, paratransgenetic approaches make insects incapable of transmitting pathogens through engineered symbiotic bacteria delivering anti-pathogenic molecules [193][194][195][196]. This approach has been successfully applied in California (USA) to block Xf infection through the glassy-winged sharpshooter Homalodisca vitripennis (Germar, 1821) (Hemiptera: Cicadellidae) and its engineered bacterial symbiont that expresses two antimicrobial peptides [197,198].
All these studies focused on the non-European strains or species, which benefit from wider genetic resources and a more detailed knowledge of the mechanisms underlying infection and resistance. Thus far, only two engineered viral vectors to induce pathogen confusion have been planned for the Apulian olive tree but results have not been reported yet [199]. Genetic engineering has not been widely used in Europe also due to the stronger negative perception of genetically modified organisms and its stricter regulation [200,201]. Additionally, the success of genetic transformation techniques has been limited by illegitimate recombination, non-mendelian inheritance of transgenes, availability of few characterized genes and the quantitative nature of some trait of interest [191,192,200]. Caution must be taken when translating in vitro or in-house tests into expected patterns in the field, as recently shown by long-term studies in Xf [187]. The advent of precision breeding technology or more efficient cis-genetic tools such as engineered site-specific endonucleases (e.g., CRISP/Cas9), which edit genomes without introducing foreign genes and create non-transgenic plants, could help overcome these issues but potential implications (e.g., cascade effects on trophic networks and/or ecosystem functioning, legitimacy to engineer wild species, etc.) should be carefully considered [191,192,[200][201][202][203].

Spatial Ecology Approaches
Spatial ecology aims to uncover the causes and consequences of organism distribution and its changes over space and time [204,205]. Special attention has been devoted to biological invasions due to their significant adverse ecological and economic effects. However, there is no single tool or unified approach to understand and predict invasion impacts [206,207]. A deep understanding of ecological and evolutionary feedbacks between invasive species and native communities can identify targets for the eradication or minimization of alien organisms in addition or in support to field surveys [207].
Spatial ecological approaches such as Remote Sensing (RS), Geographical Information Systems (GIS) and Spatially Explicit Models (SEMs) have aided a successful identification and management of invasive species, especially when these methods were combined [6,[208][209][210][211][212][213][214][215][216]. These tools have provided crucial insights into the ecological factors influencing invasions, informed detection and surveillance of invasive species via direct or indirect indicators and predicted current or future invasion dynamics [209][210][211][212]. Most of their applications have involved alien pathogens of animals, plants and humans [208,209,[213][214][215][216]. However, these tools have proven to contribute valuably to effective management strategies in several plant diseases, for example by identifying suitable areas or ecological cues that are relevant for the alien species, detecting pathogens and their vectors, assessing crop quality for indications of biotic stressors and estimating future distributions ( [208,211,213,216], Sections 3.1-3.3). Among plant pathogens, researchers dealing with Xf have significantly benefited from spatial ecology tools first in the USA and more recently in Europe, leading to substantial improvements in the management of this concerning organism [95,[217][218][219][220][221]. However, previous studies that examined the available knowledge and management measures on Xf have generally not noticed the significant contribution of spatial ecology tools. In the following sections, we fill in this gap by introducing studies that took advantage of RS, GIS and SEM techniques for the identification and management of Xf and its insect vectors in Europe.

Remote Sensing
The early detection and updated mapping of the Xf distribution is critical to limit its spread and subsequent economic and ecological losses. RS approaches aimed to detect disease symptoms in plants include multispectral, hyperspectral, airborne digital color or video imagery (Box 2 and Table S3 [222,223]). Thanks to their high sensitivity, specificity, rapidity and affordability, these tools are useful to identify infected plants and delimit epidemic areas, particularly when scarce resources were available for extensive field surveys [224].

Box 2.
Definitions of some spatial ecology terms used in the main text.

Term Definition Area Under Curve (AUC)
A measurement of the discriminatory ability of classification models. The closer the AUC to 1, the better the predictive ability of the model.
Bayesian inference approach A method of statistical inference in which the Bayes' theorem is used to update the probability for a hypothesis as more evidence or information becomes available.
biomod2 An R package for ensemble forecasting of species distributions, enabling the treatment of a range of methodological uncertainties in models and the examination of species-environment associations.

Box 2. Cont.
Ensemble modeling It is a process where multiple diverse models are created to predict an outcome, either by using many different modeling algorithms or using different training datasets.
Hyperspectral image Differs from the multispectral image (Table S3) because it consists of hundreds of spectral bands of limited width.
Leaf Area Index (LAI) Total area of green elements (one side) in the canopy per unit horizontal ground area (m 2 /m 2 ).

Leaf chlorophyll content
The amount of chlorophyll per unit leaf area.

Leaf water content
The amount of water present in a leaf.
Maxent A software to model species niches and distributions by applying a machine-learning technique called maximum entropy modeling.
Multispectral image A stack of several digital images, each one corresponding to the radiation intensity measured in a specific wavelength interval (spectral band, Table S3), which has been reflected or emitted by ground surface elements, represented by the image pixels.

Normalized Difference Vegetation Index (NDVI)
Commonly used index calculated as (NIR-RED)/(NIR+RED), with NIR representing the reflectance measured at the near-infrared wavelength (750-950 nm, strongly reflected by vegetation, Table S3) and RED the corresponding at the red wavelength (620-700 nm, absorbed by vegetation, Table S3).

Phenocam
A digital camera set up at a fixed location to capture time-lapse images.

Solar induced chlorophyll fluorescence (SIF)
The reemitted solar light in the 650-850 nm range (Table S3) from the chlorophyll-a pigment, which is linked to the initial steps in photosynthesis.

Spectral reflectance
Generic term to express the reflecting ability of a surface at a given wavelength of the electromagnetic spectrum. The distribution of spectral reflectance along the spectrum is called "reflectance curve" which has a characteristic shape for different kinds of surfaces (i.e., water, soil, vegetation, etc.).
Thermal image Digital image representing the amount of radiation emitted by a surface in the thermal infrared range (7000-13,000 nm, Table S3), which is related with the temperature of the surface itself and other characteristics (i.e., emissivity).

True Skill Statistics (TSS)
It compares the number of correct forecasts minus those attributable to random guessing to that of a hypothetical set of ideal forecasts. It considers both omission and commission errors and success as a result of random guessing; its values range from −1 to +1, where +1 corresponds to perfect agreement and zero or less to a performance no better than random.

Vegetation index
It is a spectral transformation of two or more bands designed to enhance the contribution of vegetation properties and allow reliable spatial and temporal inter-comparisons of terrestrial photosynthetic activity and canopy structural variations.
In Xf, all RS studies have been limited to olive orchards in Southern Apulia thus far (Italy [225][226][227][228][229]). Here, plant external physiological alterations associated with Xf infections (e.g., canopy defoliation, leaf wilting and chlorosis) have been surveyed with a wide array of the state-of-the-art RS technologies based on passive and active optical imaging sensors (i.e., LIDAR, Laser Imaging Detection and Ranging) from different platforms, including satellites imagery, airplanes and drones [221,225,226]]. By analyzing differences in canopy density, damages due to biotic and abiotic stresses can be generally identified by means of indices such as the Normalized Difference Vegetation Index (Box 2 [230,231]. These methods are based on conspicuous alterations in physiological structure and functions of plants but they have not been used to survey insect vectors, which do not visibly affect plants health. Host plants presenting conspicuous stress indicators (leaf browning, wilting, chlorosis and desiccation) can be discriminated by Sentinel-2 satellites through spectral shifts [221,225,226,232,233]. This satellite imagery has been widely used to support analysis of land use and relevant geophysical variables (e.g., leaf area index, leaf chlorophyll content and leaf water content, Box 2). Compared to other satellite imaging methods (e.g., Landsat or MODIS), Sentinel-2 reports variations in vegetation characteristics with higher accuracy and narrower temporal intervals [234]. and allowed to monitor the spread of Xf in olive orchards at unprecedented resolution [225,226]. Such satellite-based monitoring has significantly supported Xf surveillance across large areas and informed management decisions [224]. However, satellites may detect Xf infections only in presence of pronounced symptoms [221,235]. Subtle variations in foliage composition and color are better identified by high-resolution airborne hyperspectral and thermal images (Box 2 and Table S3), which currently represent the best techniques to detect the early symptoms of Xf infection in plants [221].
Thermal anomalies occur earlier than foliage damages in density and chlorophyll content, but the required spatial and spectral resolution (wavelength range 7-12 µm, Table S3) is obtained only through airborne sensors [221,228]. A promising alternative could be the emerging Unmanned Aircraft Systems, which allows for cheaper and improved resolution mapping [236]. This method would enable a wider application of RS approaches, as recently shown in Italy [221,227]. This approach can be transferred to other regions dealing with Xf outbreaks.
It is worth mentioning proximal sensing by means of LIDAR sensor, which has been used to examine asymptomatic and symptomatic olive trees through field robots [228]. While this approach could not capture Xf presence successfully [228], the use of low-cost field robots coupled with more powerful sensors is a promising approach to inspect large areas continuously.
Advanced physiology-focused RS methods are critical in supporting large-scale Xf surveillance and inform management efforts [225,226,237]. For example, hyperspectral observations of the top of canopies at very fine spectral resolution (Table S3) by advanced airborne sensors (e.g., FLEX and "HyPlant" [238]) may detect variations in vegetation fluorescence (Box 2 [237][238][239]). This signal could be used to rapidly screen large areas and monitor Xf infections.
Finally, a detailed testing of RS applications to Xf in a wide array of environmental settings (e.g., managed or unmanaged orchards, soil type and morphology and environmental variables) and disease symptoms in different hosts is currently unreported. Then, advanced airborne sensors, integrated approaches and multivariate models including structural, spatial and spectral data are crucial for an early and effective identification of plants infected by Xf.

Geographic Information Systems
GIS tools store and analyze spatial information obtained from RS and field surveys in easily accessible databases, which have been widely used for the early and rapid assessment of Xf and its insect vector. Geo-referenced occurrences of Xf and/or its insect vectors combined with environmental variables have led to fine-scale maps that make it possible to identify the relevant abiotic and biotic factors underlying their presence and spread, predict invasion risks and estimate their epidemic dynamics to aid future surveillance campaigns. The first and currently most complete GIS database is provided by the Apulia region in a webGIS platform (SIT.Puglia, Table 2). This information is linked to a high-resolution digital map (orthophotos of 50 cm pixels) based on a soil usage, which, as to March 2020, is divided in three management areas: infected, containment and buffer zone (http://webapps.sit.puglia.it/freewebapps/DatiFasceXF/index.html).
Similarly, the French national database stores information collected during 2015-2019 in the south of France and Corsica ( [240,241], Agence nationale de sécurité sanitaire de l'alimentation, de l'environnement et du travail, ANSES; Table 2). Spain too has developed a national catalogue of Xf occurrences (Ministerio de Agricultura, Pesca y Alimentación, MAPA; Govern Illes Balears, GOIB; Table 2). Beyond the national borders, European and world databases of Xf georeferenced records are available from several authorities ( [242], European Food Safety Authority (EFSA), EPPO and Global Biodiversity Information Facility (GBIF), Table 2).
Numerous studies have benefited from this valuable GIS information to investigate the current and future distribution of Xf and its insect vectors [211,[243][244][245][246] (discussed in Section 3.3). However, not all national databases are freely accessible (e.g., the French database, Table 2). This barrier to data sharing strongly limits our ability to respond to rising biological invasions, particularly for highly threatening alien species that require immediate quarantine measures such as Xf [247]. A rapid dissemination of any invaders' information is key to deal with these extremely successful invasive pathogens. We highlight the need to move urgently beyond restricted data access to prevent efficiently further Xf outbreaks.
While RS and SEM approaches have been predominantly focused on Xf distribution, GIS studies mostly included Xf insect vectors, particularly P. spumarius [248,249]. Although methods to estimate abundance and the olive growth-related distribution of P. spumarius are based on a sweep net technique [250,251], new GIS approaches can be used to monitor Xf insect vectors at multiple time and spatial scales [248]. Integrated field surveys and GIS data indicated that P. spumarius prefers wild herbaceous plants and switches to olive trees only at an advanced developmental stage and when the summer season progresses and herbaceous plants dry out [248,249]. Agricultural landscape mosaics can drive the spatial distribution and activity of this spittlebug and influence local Xf outbreaks [249]. In particular, the presence of suitable habitats such as herbaceous and shrub plants nearby olive orchards can enhance the emergence and persistence of Xf infections [244].
Land-use/land-cover analysis and GIS applications analyzed the landscape transformation following Xf outbreaks which resulted in compromised agroecosystems, land abandonment and reduced ecosystem services ( [252], further details at https://land.copernicus.eu/pan-european/corineland-cover). To avoid these critical issues, GIS tools can assist the management of Xf infected areas not only in agricultural terms but also considering landscape and urban areas, for example to limit the spread of Xf vectors [253].
GIS approaches enable retrospective analyses to clarify the natural and anthropic factors affecting invasive species and guide the management of future outbreaks [252]. The integration of additional environmental variables (e.g., slope, altimetry and level and kind of management) would further aid the management of Xf and its insect vectors, for example by identifying habitat features facilitating insects' movements across agricultural landscapes. However, this kind of GIS analyses in Xf and its insect vectors have been limited to Italy [248,252] while this bacterium is threatening many other countries in Europe.

Spatially Explicit Models
SEMs investigate the impact of factors such as biotic and abiotic variables on spatial ecological processes (e.g., species distribution, dispersal and extinction) and assess the potential response of organisms to these factors [254]. Since the Xf emergence in southern Italy, several SEM studies have estimated the current and future distributions, abundance and population dynamics of Xf and its insect vectors in Europe, mostly taking advantage of correlative and mechanistic models [214].
Ecological niche models (ENMs) are statistical frameworks that link the geographical distribution of species (georeferenced presence records; e.g., individual locations and species' presence) to spatial variation in environmental conditions (geographic layers of environmental information; e.g., climate, land cover and soil attributes), identify the factors that enable species presence and predict potentially suitable areas [255]. ENMs estimate fundamental niches of species, and are applied when the aim is to know the potential distribution, as in the case of invasive species or projections in space and time. In European ENMs studies of Xf, presence records were obtained from online databases ( Table 2) and published data [98], sources of carefully verified information that is available thanks to collaborations between researchers and farmers. Environmental predictors were selected through statistical methods [13,14,86,219,220,240,243,256,257]. The most widely used predictors include elevation (WorldClim and European Environment Agency (EEA), Table 3), land cover/land use (Copernicus, ESA, Table 3) and climatic variables (WorldClim and Climatologies at High resolution for the Earth's Land Surface Areas (CHELSA), Table 3) used at an equatorial resolution of ca. 1 km [219,240,243,256,257], 4.5 km [14,86,220,258] or 10 km [13]. A high (ca. 1 km 2 ) or very high (<1 km 2 ) resolution should be preferred due to Xf 's small spatial range, particularly in the early stages of spread. Similar to other invasive species [9], the potential distribution of Xf and its insect vectors in Europe have been mainly estimated using the free software packages Maxent and biomod2 (Box 2 [13,14,129,258]). ENMs are validated through statistical analysis and field survey. For Xf, the area under curve (AUC) and the true skill statistics (TSS) have been widely used but no ground validation studies have been performed in Europe thus far (Box 2 [13,14,219,220,240,243,256]).
According to AUC and TSS (Box 2), the most accurate modeling algorithms used in Xf studies are Maxent, Artificial Neural Networks, Generalized Additive Models and Random Forests [13,220]. More recent developments aim to improve modeling efficiency by integrating multiple algorithms (ensemble models, Box 2 [259]), as recently used in Xf in Europe [13,220,258]. However, these ENMs are reliable only if they are accompanied by a careful parameterization of each algorithm/model underlying the ensemble, a step that is often overlooked. Researchers working with ENMs need to have a comprehensive understanding of how each algorithm works and how each algorithm responds with their particular data.
ENMs analyses indicated that the main descriptors of Xf current distribution are the minimum temperature of the coldest month, mean temperature of the warmest quarter and precipitation seasonality [95,219,220,240,243]. Xf had a high probability of occurring in the European areas characterized by warm temperatures and low rainfall during summer and mild winter temperatures [95,220]. In line with these predictions, the minimum threshold temperature for growth of Xf in plants is between 12 and 18 • C [217] and temperatures < 6 • C seem to kill this bacterium [260]. The mean temperature of driest and warmest quarter is likely related with the ability of the bacterium to support water stress or warm conditions [240]. Overall, the coldest period (month or quarter) is the most significant variable underlying Xf distribution in Europe [13,[218][219][220]240,258] and worldwide [218]. More specifically, a ENM study investigated the relationship between the potential distribution of each of the three Xf subspecies occurring in Europe and: (i) mean temperature of driest quarter and annual mean temperature (Xf subspecies fastidiosa); (ii) mean temperature of coldest quarter and mean temperature of driest quarter (Xf subspecies multiplex); and (iii) minimum temperature of coldest month and mean temperature of driest quarter (Xf subspecies pauca [13]). Overall, the main limiting variable is the minimum winter temperatures, with Xf subspecies pauca being particularly sensible to this predictor [13,220,243]. Xf subspecies fastidiosa showed a significant association with the mean temperature, consistently with its year-round susceptibility to warm temperatures [13,220]. Xf subspecies multiplex presented more suitable areas with warm temperatures during the dry season but is sensitive to cold temperatures in winter [13,220]. Finally, Xf subspecies pauca prefers warm climates with mild winters [13,220,243], low elevation (0-150 m a.s.l.), intensively cultivated areas with complex patterns, oak woodland and Mediterranean shrubland [219].
According to the most recent Xf ENMs studies, the currently established distribution of this plant pathogen in Europe (Apulia, Tuscany, French Riviera, Corsica, Alicante and Balearic Islands) is notably smaller than the extent of climatically suitable areas, which includes the Mediterranean coastal areas of Spain, Greece, Italy and France and the Atlantic coastal areas of France, Portugal and Spain [13,219,220,240,258]. Interestingly, the currently climatically suitable areas identified by ENMs diverge among Xf subspecies. The Xf subspecies fastidiosa presented suitable areas in the southern regions of France, Portugal and Spain, the Balearic Islands and large regions of Italy, Croatia, Greece, Cyprus, Crete and Turkey [13,220]. Xf subspecies multiplex showed the most extended suitable areas among the three subspecies [13,219,220]. In fact, with the exception of the northwestern areas, suitable climatic conditions for this species could be found in all Europe [13,220]. Finally, the Xf subspecies pauca disclosed the most restricted suitable distribution, including southern Italy, the Mediterranean coasts, Portugal, southern Spain and France [13,14,219,220]. These results indicate that the Xf subspecies multiplex could invade most of Europe while pauca and fastidiosa could be mostly restricted to the Mediterranean basin [13,219,220]. More recently, a ENMss analysis has focused on three Xf ST (Box 1): Xf subspecies fastidiosa ST1, Xf subspecies multiplex ST6 and Xf subspecies pauca ST53 [13]. While ST1 and ST53 results did not differ from the previous analyses significantly, ST6 showed a more limited potential distribution restricted to lower latitudes [13].
ENMs can also be used to predict Xf distributions assuming future scenarios of climate change [219,258]. The future distribution of the Xf subspecies pauca in the Mediterranean basin would not be affected or shift towards north ( [219,258], LB personal obs.) while the Xf subspecies multiplex will expand and/or shift northward by 2050, possibly threatening several of the most economically important wine, olive and fruit producing regions of Europe ( [258], LB personal obs.). A potential limit of these studies is the use of only one global circulation model while multiple ones should be included, the use of presence records from an introduced range rather than the native regions [219], a spatial resolution that is not appropriate for the study system (ca. 25 km 2 ) or a limited number of validation methods to assess the performance of their models [258].
Only few modelling studies are available for the insect vector P. spumarius. The most relevant climatic predictors of its presence in Europe included the maximum temperature of the warmest month, the temperature annual range and the precipitation of the coldest quarter [86]. However, a recent study in the UK showed that the optimal conditions for this spittlebug are warm areas [261]. Additionally, deciduous woodland and improved grassland, the spittlebug's preferred habitats, could support large populations of the insect vector and pose the greatest risk of Xf spread [261]. Another SEM approach to analyze the Xf spatiotemporal distribution are the mechanistic models, in which mathematical models describe the interaction between environmental factors and functional traits of the species that affect demography and distribution [207]. These models incorporate the mechanisms of spread and are typically used for small areas at high spatial resolution (<1 km 2 ). Such models are particularly useful to identify the most successful strategies to limit further pathogen spread. This approach successfully predicted Xf spread in several areas such as southern Apulia [246,262,263] and Corsica [264,265]. In Apulia, mechanistic models showed the importance of long-distance jumps in insect vectors movement, buffer zones width and the significant role of plants other than olive trees in increasing the spread rate [266]. Network analysis suggested that attempts to eradicate Xf from southern Italy may be pointless and that this region will provide a pathogen reservoir for further spread [246]. Therefore, the goal of management strategies in this region should be to keep Xf density as low as possible, primarily through control of P. spumarius and plant removal [13]. Another mechanistic model indicated that the current insect vector dispersal is less than 1 km/year but this is expected to increase substantially due to human activity and passive dispersal through strong winds [13]. These models can also be used to estimate the age of the infection. For example, a Susceptible-Infected-Removed (SIR) model suggested Xf was probably introduced in Corsica (France) much earlier than its discovery, and the pathogen may now be difficult to eradicate. Using a Bayesian inference approach (Box 2), it was in fact confirmed that Xf was introduced to South Corsica (near Ajaccio) around 1959, long before its first detection in 2015 [264]. Integrating genetic and demographic data from North Corsica and south-eastern mainland France would clarify the occurrence of multiple introductions and human-mediated long-distance dispersal. More recently, potential strategies to limit Xf infection in Apulia through insect vectors management across its entire life cycle were evaluated using SIR models [38]. A new stochastic simulation model incorporated epidemiological information and spread patterns of Xf subspecies pauca in olive trees in Apulia between 2013 and 2018 to explore various epidemiological, introduction routes and management scenarios [266]. This approach could be extended to other countries that are potentially suitable to this bacterium [13]. Other mechanistic studies focusing on this pathogen are still ongoing [106,267]. However, quantifiable data and measurements on its spread in Europe are scarce mainly due to differences in the bacterial strain, plant host, insect vector and environment among Xf infestations [263,265,268]. Therefore, predicting the extent of the spread as well as assessing the efficiency of management measures is extremely challenging.
All SEM studies reviewed here agree that the European climate can be particularly favorable to Xf, particularly for the subspecies fastidiosa and multiplex [13,220]. Caution is warranted when interpreting models projected into new zones or in future climate scenarios that are different than the calibration area. Additionally, Xf may adapt to a wide range of hosts [269] and this rapid evolution could affect future geographical distributions yet to date this information cannot be included in SEMs [13,220,243,258].

The Impact of Climate Change
Climate change can directly and/or indirectly affect plant pathogens and their vectors [270,271]. Likely, the modification of key abiotic factors will influence the distribution of insects, plants and pathogens in future environments. The direction of this change is less clear, as it may impact each species and their relationships differently, as seen in aphids and viruses linked to wheat [272].
Climate has a major role in the epidemiology of all the Xf -induced diseases [16,273]. Climate change will differentially affect Europe [274]. In Southern and Central Europe, higher temperatures and reduced mean summer precipitation are expected to increase the occurrence of heat waves and droughts, leading to progressive desertification [275]. Accordingly, Xf subspecies pauca is expected not to expand its distribution being intolerant to future climatic conditions [219,258]. Instead, in Northern and Eastern Europe, precipitations and humidity are expected to increase [274]. Thus, Xf subspecies multiplex will probably expand its distribution to those areas [258].
Climate change will likely expand the suitable geographic range of Xf insect vectors, with the risk of spreading the bacterium to previously unaffected areas, i.e., Xf subspecies multiplex [258]. Hemipterans are expected to benefit from future CO 2 levels but the effects on sap sucking insects are species-specific [276,277]. Some species may be positively influenced (e.g., Aphis gossypii (Glover, 1877) (Hemiptera: Aphididae); Myzus persicae; Rhopalosiphum maidis (Fitch, 1856) (Hemiptera: Aphididae); and Bemisia tabaci (Gennadius, 1889) (Hemiptera: Aleyrodidae) [277][278][279][280]), others might show a negative trend (e.g., Brevicoryne brassicae (Linnaeus, 1758) (Hemiptera: Aphididae); and Acyrthosiphon pisum (Harris, 1776) (Hemiptera: Aphididae) [277][278][279][280][281]) and further sap feeders are likely to be unaffected (e.g., Aphis nerii (Fonscolombe, 1841) (Hemiptera: Aphididae); Aphis oenotherae (Oestlund, 1887) (Hemiptera: Aphididae); and Aulacorthum solani (Kaltenbach, 1843) (Hemiptera: Aphididae) [277]). Overall, the spatial distribution of many insect and pest species is often limited by unsuitable climatic conditions and geographical barriers [282,283]. However, these species may be able to overcome these climatic constraints by taking advantage of artificial environments (e.g., greenhouses [284,285]), evolutionary adaptation including emergence of resistance [286] or human-mediated dispersal [287,288]. The effects of climate change on the distribution and behavior of Xf insect vectors in Europe is not known and it is therefore difficult to predict if their detrimental effects on global food security will increase or decrease in the future. Increased CO 2 and temperature will not only impact the distribution of Xf and its insect vectors, but also their interactions with plants. In fact, plant biochemistry and physiology will most probably be altered and the response against pathogens weakened under the predicted future climatic conditions. Considering this factor, therefore, direct and indirect dispersal of Xf will likely increase [289].

Genetic Approaches
Genetic approaches have significantly assisted the successful identification and management of Xf outbreaks (as reviewed in Section 2). A major limit to their application is the absence of complete genomes and robust annotations for most Xf strains and the large amount of uncharacterized genetic diversity [290]. Genetic protocols may not be equally effective in distantly related lineages but need to be adapted to different strains and hosts [110,168,187,291]. Automated annotators are benefiting from expanding datasets but this technology has not been updated in the last decade and its accuracy has decreased following the accumulation of fragmented genomes [292]. Recent technological developments such as the state-of-the-art whole genome/transcriptome sequencing, upgraded algorithms addressing several assembling issues (e.g., repetitive regions [293]) and CRISPR-based editing tools offer versatile and cost-effective procedures to obtain exploratory and diagnostic genetic tools that are specifically tailored for the analyzed sample. Gene regulatory networks can be explored through RT-qPCR [294,295] and Nanopore direct full-length RNA sequencing [296].
However, these cutting-edges tools and vast datasets are not easily accessible to non-bioinformatician researchers. A unified open-source workflow for automated whole bacterial genome studies [297] would enable a quick and effective identification and monitoring of outbreaks also from scientists without a computational background. At the same time, it would provide a wide amount of standardized genomic data for finer analyses.
The increasing global trade of plants could lead to further spread of this pathogen in new areas or to the introduction of novel genotypes in existing populations. Due to the high recombinant ability of this organism, sympatric populations of multiple bacterial strains could hybridize and generate more niche-adapted and/or virulent recombinants, as seen in South America [59,66,107]. Then, it is crucial to limit the import-export of susceptible or potentially infected material and to extensively screen plants for novel pathogen introductions and hybridization. The European areas that already host multiple strains of Xf (Corsica and Baleares, Figure 4) are at the highest risk of emergence of new variants and should be closely monitored. A more detailed knowledge of the hybrid strains and the recombination potential could possibly help to develop more efficient diagnostic protocols but also to create and/or spread strains that are avirulent or have limited ability to hybridize.
Another unexplored genetic mechanism in Xf is epigenetic variation (Box 1), which is known to regulate virulence in bacteria. We advocate future studies to illuminate the contribution of this kind of polymorphism to Xf evolution and pathogenicity, for example by methylation sequencing.
diversity [290]. Genetic protocols may not be equally effective in distantly related lineages but need to be adapted to different strains and hosts [110,168,187,291]. Automated annotators are benefiting from expanding datasets but this technology has not been updated in the last decade and its accuracy has decreased following the accumulation of fragmented genomes [292]. Recent technological developments such as the state-of-the-art whole genome/transcriptome sequencing, upgraded algorithms addressing several assembling issues (e.g., repetitive regions [293]) and CRISPR-based editing tools offer versatile and cost-effective procedures to obtain exploratory and diagnostic genetic tools that are specifically tailored for the analyzed sample. Gene regulatory networks can be explored through RT-qPCR [294,295] and Nanopore direct full-length RNA sequencing [296].
However, these cutting-edges tools and vast datasets are not easily accessible to nonbioinformatician researchers. A unified open-source workflow for automated whole bacterial genome studies [297] would enable a quick and effective identification and monitoring of outbreaks also from scientists without a computational background. At the same time, it would provide a wide amount of standardized genomic data for finer analyses.
The increasing global trade of plants could lead to further spread of this pathogen in new areas or to the introduction of novel genotypes in existing populations. Due to the high recombinant ability of this organism, sympatric populations of multiple bacterial strains could hybridize and generate more niche-adapted and/or virulent recombinants, as seen in South America [59,66,107]. Then, it is crucial to limit the import-export of susceptible or potentially infected material and to extensively screen plants for novel pathogen introductions and hybridization. The European areas that already host multiple strains of Xf (Corsica and Baleares, Figure 4) are at the highest risk of emergence of new variants and should be closely monitored. A more detailed knowledge of the hybrid strains and the recombination potential could possibly help to develop more efficient diagnostic protocols but also to create and/or spread strains that are avirulent or have limited ability to hybridize.
Another unexplored genetic mechanism in Xf is epigenetic variation (Box 1), which is known to regulate virulence in bacteria. We advocate future studies to illuminate the contribution of this kind of polymorphism to Xf evolution and pathogenicity, for example by methylation sequencing. . Spatial distribution of the European Xylella fastidiosa multi-locus sequence type (ST) that are at higher risk of hybridization due to geographical proximity or sympatry (details in the main text, data from Table S1 [32,98]).

Spatial Ecology
RS, GIS and SEMs methods have been rather scarcely used to investigate spatial patterns of Xf and its insect vectors so far, largely due to low resolution as well as difficult data accessibility and processing (as reviewed in Section 3). However, spatial ecology approaches are important tools for implementing effective management strategies. For example, they can help in obtaining a detailed knowledge of the biology and ecology of this pathogen and its insect vectors, identify early warning signals and optimize resources for screening [298][299][300].
Current hyperspectral satellite sensors (Box 2) lack the spatial resolution to distinguish individual tree foliage. Xf early detection would benefit from high (submeter) spatial resolution imaging spectroscopy, thermal data (Box 2) or aerial photos to assess subtle changes in spectral traits, a technology that can be potentially deployed at large scales with automated airborne platforms. Expected improvements include higher availability of multispectral optical imagery with increasing spatial, spectral and temporal resolutions, based on new satellites and sensors. Another area of development focuses on unmanned aerial vehicles and phenocams (Box 2) with increasing multispectral, hyperspectral and thermal imaging, as well as LIDAR. These tools would enable better estimates of the vegetation index (Box 2) thanks to their semiautomatic discrimination between physiological and pathological changes [301]. The resulting insights will refine our capacity to predict, detect and assess the distribution and impact of Xf and its insect vectors in Europe with unprecedented efficiency.
GIS developments are moving towards integrating additional territorial variables of interest (e.g., slope, altimetry, level and kind of management, insect vector distribution) and real-time geographic information [302] thanks to expanded fixed and mobile (e.g., on vehicles or pedestrians) sensor networks. More accessible data, as well as improved, free and user-friendly GIS software, are going to advance our ability to analyze Xf and its insect vectors, playing an important role in preventing the spread of Xf in Europe.
All SEM methods analyzed here could aid surveillance programs for Xf in Europe, for example by prioritizing areas to be inspected according to their climatic suitability. ENMs could inform the suitability of recipient areas before species are introduced and mechanistic model evaluate population changes over time according to landscapes. Future SEMs integrating insect vectors' distribution and/or population dynamics as predictors in Xf models as well as using several global circulation models to estimate their potential distributions in climate change scenarios (e.g., BCCCSM1-1, CCSM4, GISS-E2-R, MIROC5, HadGEM2-ES and MPI-ESMLR from WorldClim, Table 3) will likely help pursue these aims. Importantly, the early identification of invasive species is crucially based on obtaining the relevant information in a timely manner. As seen for Xf, a clear picture of ongoing invasions is too often acquired only several years after the first report of the responsible organism. This issue partly arises from researchers relying only on well-parameterized and validated results (e.g., ENMs), a time-consuming process that can negatively affect an effective early response. To pursue this aim, preliminary models too can inform initial risk assessments [6]. SEM tools are particularly useful in screening procedures and decision-making processes related to Xf and insect vectors management as they allow exploring, comparing and assessing the implications of different scenario in terms of likely success, cost and needed efforts before possibly implementing them.
Importantly, all the RS, GIS and SEMs methods heavily depend on data quality. Therefore, they need to be guided by a solid background on biology and ecology of Xf and insect vectors and by long-term research and fieldwork.

Integrating Genetic and Ecological Approaches
Multiple levels of biodiversity, from ecosystem to genetic diversity, need to be comprehensively included in studies that aim to accurately identify, understand, predict and prevent the spread of invasive species [42][43][44]. For example, genetic distances overlaid with environmental variables, as typical of landscape genetics studies, have illuminated current and historical movement of organisms, the underlying genetic and abiotic factors, potential colonization routes and predicted distribution shifts under future climate scenarios [43]. Additionally, maps are influential outreach instruments that are easily understood by non-experts such as citizens or conservation agencies and have greatly enhanced their understanding of the role of genetic and environmental factors on invasive species successes [43]. Then, an interdisciplinary approach combining diverse perspectives such as genetic and ecological data is one of the most powerful tools in conservation biology and epidemiology.
Accordingly, more recent studies focusing on Xf seek to genetically inform spatial ecology approaches ( Figure 2). The early detection of this pathogen was improved by the integration of photointerpretation of aerial/satellite images, accurate onsite data acquisition and rapid on-site pathogen detection in plant material and insect through genetic screenings [303].
Infected asymptomatic plants and insect vectors at earlier stages of infections were identified thanks to the combination of advanced PCRs, RS and SEMs leading to a highly precise spatial map of disease incidence ( [86,225], methods described in Sections 2 and 3). Additionally, an undisclosed higher genetic diversity and co-occurrences of Xf strains were identified within insect vectors, highlighting that our knowledge of this pathogen has been largely underestimated [86]. These insects hosting multiple, sympatric bacterial populations may provide other recombination opportunities, increasing the risk of emergence of new Xf variants through hybridization. These studies exemplify how combining different methodologies from diverse fields lead to a more accurate and effective detection and management of Xf.
In addition to these successful first combined approaches, incorporating interdisciplinary information more closely beyond simple reciprocal corroboration can provide even more effective and thorough tools for detecting and monitoring Xf. For example, genetic screenings alone allow identifying precisely Xf but are expensive in terms of both time and money (details in Section 2); spatially explicit ecological tools can delimit the area to be monitored. Inversely, environmental DNA can inform on multiple species (pathogen, plants hosts, insect vectors) distribution and abundance in a single sampling and reduce the costs associated to fieldwork and geographical imaging [42,43,304]. Transcriptomic profiling and remote sensing can supervise changes in crucial traits such as tree pigments levels in large areas at fine resolution, allowing the early identification of Xf infections. The distribution of genomic diversity and its association with climate maps and digital elevation models can identify adaptive variants and inform predictive modeling of host or pathogen responses to environmental factors such as climate changes [42,43,304]. These integrated models could shed light on the influence of external environment on epigenetic regulation, an information that would address a long-standing question in evolutionary biology and could help the management of this invasive species. A closer integration of Xf adaptive and hybridization potential as well as population size extrapolated from genomic and epigenetic data directly in the ENMs frameworks would greatly improve the prediction of outbreaks and climate responses (hybrid ENMs, e.g., [44]).
A critical limit of these integrated studies is the need of large amounts of genetic and ecological data as well as highly efficient statistical and computational resources to analyze these vast datasets and provide realistic projections. Most recent technical innovations and changes in data sharing policies are helping to affordably obtain large multidisciplinary datasets although some methodological and logistical challenges remain to be addressed [44]. Software packages that are more easily accessible to researchers from different fields would facilitate such integrations, which would be supported also by a more extensive data standardization (e.g., same resolution or format), an issue that is well known by bioinformaticians.
We advocate theoretical and empirical developments towards comprehensive interdisciplinary approaches to enable the integration of diverse perspectives that will improve our understanding of Xf and other invasive species and thus our ability to properly deal with this pathogen.

Conclusions
Xf is a highly threatening pathogen largely affecting the European agriculture. However, we are far from a comprehensive understanding of this invasive species, its hosts, ecological niches and their interactions. While previous management strategies addressing Xf have mainly relied on a monodisciplinary approach, we emphasize the importance of combining different perspectives such as genetic and spatial ecology and beyond. Including different, complimentary outlooks will substantially improve both our knowledge on invasive processes and resources allocation and thus significantly optimize diagnostic and management efforts. Even more imperative is an effective communication with stakeholders, especially those outside the academic community [305,306], to aid a successful and sustainable management of one of the most dangerous plant pathogens. This organism also offers an exceptional study system to explore the evolution of pathogenicity at different levels of biological complexity, from molecules to ecosystems, and to inform the development of integrated approaches for the identification and management of this and other invasive species.
Supplementary Materials: The following are available online at http://www.mdpi.com/2071-1050/12/11/4508/s1, Table S1. List of the available genomes for the European Xylella fastidiosa strains and associated features (data from NCBI). Strains without genomic resources are not included. Table S2. List of Xylella fastidiosa insect vector species (order Hemiptera) included in the BRIGIT project for which genomic resources will be available soon (details in the main text). Table S3. Spectral bands for the SENTINEL-2 sensors (S2A and S2B).