Occurrence of Pseudomonas syringae pvs. actinidiae, actinidifoliorum and Other P. syringae Strains on Kiwifruit in Northern Spain

Pseudomonas syringae pv. actinidiae (Psa), the agent causing bacterial canker of kiwifruit, has been present in the Principality of Asturias (PA), Northern Spain, since 2013, although with restricted distribution. In this study, 53 strains collected in kiwifruit orchards in PA during the period 2014–2020 were characterized by a polyphasic approach including biochemical and phylogenetic analysis. Thirty-three strains, previously identified by PCR as Psa, have been found to be a homogeneous group in phylogenetic analysis, which seems to indicate that there have been few introductions of the pathogen into the region. Two strains were confirmed as P. syringae pv. actinidifoliorum (Pfm), so this is the first report of Pfm in the PA. The remaining 18 strains were found to be close to P. avellanae and P. syringae pv. antirrhini or to strains described as Pfm look-alikes. Pathogenicity tests carried out on peppers with a selection of strains have shown that both Psa and Pfm caused clear damage, while the 18 atypical strains caused variable lesions. It would be necessary to carry out pathogenicity testing of atypical strains on kiwifruit plants to study the role of these strains in the kiwifruit pathosystem to evaluate their pathogenic potential in this crop.


Introduction
Kiwifruit (Actinidia deliciosa (Chev.)Liang and Ferguson), although originally from China, is now cultivated in many different zones of the world, including in New Zealand, Chile, the USA, Korea, Japan and several European countries (Italy, Greece, France, Spain and Portugal).In Spain, kiwifruit cultivation was first introduced in the northwest region at the end of the 1960s.In the Principality of Asturias (PA), which is situated in the center of Spain's northern coast, cultivation of the fruit was not established until 1972-1973, and it was not until the 1980s that kiwifruit orchards became relatively widespread.The latest data highlight the PA as the second largest kiwifruit producer in Spain, with a total of 280 ha producing up to 3823 t.These data make kiwifruit the second most important commercial fruit crop in this region, preceded only by apple (cider and dessert apple) with a total production of 13,086 t, despite variations due to the pronounced biannual alternation of this species [1].
Diseases have been widely observed in this crop.Due to the recent domestication of kiwifruit, in clear contrast with other, more established crops (e.g., vineyards or apples), the emergence of fungal and bacterial agents affecting kiwifruit plants was documented once commercial kiwifruit orchards became significant [2].
Among its bacterial pathogens, Pseudomonas viridiflava and P. syringae pv.syringae cause leaf spots and bud and flower rot that decrease production [3][4][5][6].P. syringae includes more than 60 well-characterized pathovars including two that have been described as pathogenic agents specific to this host: P. syringae pv.actinidiae (Psa) and P. syringae pv.actinidifoliorum (Pfm) cause bacterial canker and bacterial spots, respectively [7][8][9].Bacterial canker is the most damaging disease at present; it is considered the greatest risk to kiwifruit culture and is responsible for both production and plant losses.
Different biovars of Psa are dispersed around the world and were initially classified on the basis of the multilocus sequence analysis (MLSA) method, the detection of type III secretion system effector genes and the presence of phaseolotoxin and/or coronatine phytotoxins [26][27][28][29][30]. Biovar 1 (Psa 1) includes strains detected in Japan that produce phaseolotoxin; biovar 2 (Psa 2) includes Korean strains that produce coronatine; biovar 3 (Psa 3) includes the Italian aggressive strains that do not produce phytotoxins and were responsible for the 2008 epidemic outbreak; biovar 4 (Psa 4), currently considered by Cunty et al. [9] to be a new pathovar, P. syringae pv.actinidifoliorum (Pfm) which includes low virulence strains; biovar 5 (Psa 5) includes Japanese endemic strains isolated in 2012 without producing phytotoxins; and biovar 6 (Psa 6), found in Japan in 2015, could also be an endemic lineage similar to biovar 5 but differs in that it produces phytotoxins, phaseolotoxin and coronatine [29,31].
Among the different Psa biovars, the most virulent is biovar 3, also named Psa-V, causing severe damage in kiwifruit orchards worldwide [9,10,27,32].In 2013, McCann et al. [2] sequenced the genome of more than 30 Psa strains, concluding that Psa-V isolates from different sources are very similar in phylogenetic terms, with the exception of a Chinese isolate which is also virulent.Currently, three different clades are considered in Psa with inter-and intra-clade pathogenic variability [33].
This group of virulent bacteria is of worldwide interest as a model of study, as it is a disease that has been described recently, which allows its diversification and evolution patterns to be studied [2].Furthermore, its great ease of dissemination, the difficulty of controlling it and the considerable crop damage it causes are other reasons for subjecting it to detailed investigation.
In Europe, Psa has been included in the EPPO (European Plant Protection Organization)'s A2 list of pests recommended for regulation as quarantine pests since 2012 [34].In Spain, bacterial canker affecting Actinidia spp. was reported for the first time in Galicia [18,19].Official information shows that Psa was detected in PA in 2013; later, an outbreak detected in the Basque Country was eradicated in 2015.Psa was also detected in Navarra in 2018, and in 2019, it was again detected in the Basque Country and Asturias, with the infected plants being destroyed in the Basque Country.In 2020, two outbreaks were detected, one in Asturias and the other in Navarra [35][36][37][38][39].In addition, Pfm was identified in Galicia [40], and a recent study [41] cites the presence of a bacteria actinidifoliorum look-alike, pending confirmation, in Spain, although without specifying the region.
In this study, we present the characterization and identification by means of biochemical and genetic approaches of Psa and other related P. syringae strains isolated from kiwifruit orchards in PA over the years 2014-2020.The results obtained provide information about the taxonomy and phenotypic variation of the strains involved in kiwifruit canker, which will be useful for further epidemiological studies and for strengthening biosecurity strategies in the management of bacterial pathogens in kiwifruit orchards.

Bacterial Isolates
A total of 53 bacterial isolates obtained from samples of kiwifruit analyzed between 2014 and 2020 in the PA, specifically on A. deliciosa cv.'Hayward' from commercial kiwifruit orchards and one domestic orchard, as well as one on cv.'Tomuri' male vines, were investigated in the present study.Of these, thirty-three had previously been identified as Psa according to EPPO Diagnostic Protocols PM7/120 (1) [42], eighteen of which were isolated in 2014, six in 2015, three in 2018 and the remaining six in 2020.In addition, 20 strains isolated in 2014 were included in this study, since they amplified one of the two bands expected in the PCR test developed by Gallelli et al. [43] and the single band expected in the PCR developed by Rees-George et al. [44].
Bacterial isolates were coded and preserved by two different methods, at −80 • C with 50% DMSO (Dimethyl Sulfoxide), and freeze-dried, in the collection of the Phytopathology Laboratory of SERIDA (LPPA).

Detection of Genes Involved in Toxins and Levan Production 2.2.1. Detection of Toxin Genes
Fragments of the tox-argK gene cluster involved in the biosynthesis of the phaseolotoxin and the cfl gene encoding coronatine were amplified [49][50][51].As controls, P. syringae pv.tomato CECT (Spanish Type Culture Collection) 126, LPPA 696, P. syringae pv.phaseolicola LPPA 597 and LPPA 1407 were used.

Detection of Levansucrase Genes
The analysis of the distribution of lsc genes in the isolates under study was carried out by PCR amplification with specific primers.lscA/B/C genes were analyzed using LPPA 28 and LPPA 275 as controls, respectively [55].

Amplification of Housekeeping Genes
Three genes were amplified: the gltA gene, also known as cts encoding citrate synthase, rpoD, encoding RNA polymerase sigma factor and gyrB, encoding the B subunit of DNA gyrase [56,57].The amplified fragments were sequenced by Secugen (Spain) or Eurofins (Germany).Sequences were submitted to GenBank and the accession numbers are shown in Supplementary Table S1.The sequences were compared with those deposited in databases through the BLAST algorithm [58].

Phylogenetic Analysis
Sequences of concatenated and individual amplified genes were aligned using Clustal W [59], and phylogenetic trees were constructed based on each individual locus and the three concatenated genes by the maximum likelihood method with the Tamura-Nei model.Their topological robustness was evaluated by bootstrap analysis based on 1000 replicates using MEGA 11 software [60].Pseudomonas asturiensis LPPA 221 T was used as the outgroup.
Sequences of Psa, Pfm and other pvs. of P. syringae obtained from the GenBank databases were included for comparison in the phylogenetic analysis (Supplementary Table S1).
The number of segregating sites (S) and the mean of the nucleotide diversity (π), defined as the average number of nucleotide differences by site between sequences of the whole population, were calculated both for the individual genes and the concatenated sequences, also using the Kimura two-parameters model in MEGA 11.Data on the evolutionary divergence of Psa-Pfm and strains of P. syringae are shown in Supplementary Tables S2 and S3 [60].

Pathogenicity Assays
Pathogenicity was tested by inoculation of the isolates on marketable fruits of Capsicum annuum cv.'California' according to Abelleira et al. [40].Twelve LPPA strains (2727, 2983, 3697, 3700, 2668, 2669, 2723, 2725, 2757, 2758, 2759, 2760) were selected for the assay.One pepper per strain, previously cleaned with ethanol, was inoculated by ten 2 mm punctures made in the epidermis and then filled with 15 µL of a bacterial suspension of 10 8 cfu/mL.Strains NZ 10627 and CFBP 8039 as positive controls and 15 µL of Luria broth as negative control were used.Inoculated peppers were kept at 25 • C for 7 days.Koch's postulates were fulfilled by identification of the re-isolated bacteria from the symptomatic lesions on inoculated peppers.

Bacterial Isolates
Bacterial strains recovered over 2014-2020 in several kiwifruit orchards in PA were first grouped into Psa (33 strains) and atypical strains (20 strains) according to the results of the two PCR tests performed.Atypical strains amplified the expected band (280 bp) in the simple PCR [44] but only the upper one (492 bp) in the duplex PCR test [43].This result is consistent with that described for Pfm [9]; however, our study indicates that only two out of twenty strains can be assigned to this pathovar.Data for the total of 53 strains are shown in Table 1.
Strains LPPA 2983-2985 were isolated from A. deliciosa cv.Tomuri, while all others were isolated from A. deliciosa cv.Hayward.
The presence of Psa in Asturias was only notified by the Plant Health Service of the Government of PA in 2013 in the municipalities of Pravia and Langreo [35].Some other locations have been added since then: the area of Salas in 2018 [37], the area of Carreño in 2019 [38] and Piloña and Villaviciosa in 2020 [39].However, there is no indexed bibliography publishing these results except for the recent reference by Moran et al. [41].These authors describe new genetic lineages of a P. syringae pv.actinidifoliorum look-alike isolated from the main kiwifruit-producing areas in the north and east of Spain, although without providing details of the locality of origin of the strains under study and without revealing the presence of either Psa or Pfm in the series they collected.Therefore, this is the first article in which the strains collected in the PA have been studied and in which isolates of both Psa and Pfm have been found.

Phenotypic Characterization of the Isolates
All the strains analyzed were negative for Gram stain, oxidase and arginine dihydrolase and positive for hypersensitivity on tobacco leaves, oxidation of glucose and also the use of sucrose and xylose as sole carbon sources.The use of sucrose as a source of carbon rules out that any of the atypical strains correspond to P. viridiflava, which, as was mentioned in the introduction, is one of the bacteria responsible for damage in kiwifruit.Typical P. viridiflava strains are levan negative, but strains of these bacteria isolated from kiwifruit and other crops have been described in Asturias as levan positive.However, these levan-positive strains remain negative when utilizing sucrose as the sole carbon source [6].
Concerning the Psa group, phenotypic characterization evidenced that the series under study shows differences with those described by other authors in several characteristics, such as the use of xylose, trigonelline, mannitol, sorbitol and inositol as the only source of carbon and casein, tween 80, gelatine and esculin hydrolysis [7,11,13].However, Abelleira et al. [40] have already found a virulent strain of Psa that was esculin positive.
The LPPA 2768 and LPPA 2769 strains have been identified in the present study as Pfm on the basis of the biochemical and physiological tests, as both showed fluorescence and were esculin positive.
The Ps group showed greater heterogenicity.However, all of them were esculin and gelatin positive, and only one strain was not fluorescent.Biochemical features of the isolates under study are compiled in Table 2.These results would emphasize the hypothesis that at least most of the atypical strains collected in 2014 would not correspond to Psa.Furthermore, fluorescence appeared to be a key factor for selecting Psa strains, since 100% of those identified as such were nonfluorescent compared to 5% of the atypical strains.

Phytotoxin Gene Detection
None of the 53 strains studied amplified the gene fragments corresponding to phaseolotoxin or coronatine toxins, whereas the reference strains P. syringae pvs.tomato and phaseolicola were found to be positive.These results are consistent with what is expected for virulent Psa strains corresponding to biovar 3.
Phytotoxins were considered the major virulence determinants in P. syringae [54] and frequently their production increases plant damage.Coronatine and phaseolotoxin generally induce chlorosis, whereas syringomicin and syringopeptin cause necrosis [61].
None of the 53 strains amplified syrB, syrD or sypA.On the basis of these results, we would rule out that any of the atypical strains are P. syringae pv.syringae, which is another of the bacteria described as pathogenic in kiwifruit [4,5].

Levansucrase Gene Detection
Levan production is a significant taxonomic feature of the fluorescent Pseudomonas group to which P. syringae belongs and is considered a virulence factor in diseases caused by phytopathogenic bacteria [62].Levan is a high-molecular-weight polysaccharide synthesized by several levansucrase isoenzymes (Lsc), so the presence of genes involved in the production of these enzymes was studied.The lscB gene is located in a plasmid and produces an extracellular enzyme, and the lscC gene is positioned on the chromosome and its product is a periplasmic levansucrase, while lscA is considered a cryptic gene with chromosomal location and is not related to levansucrase activities [55].It has been suggested that lscA may be an ancestral levansucrase gene that lost the ability to express the enzyme in P. syringae [63].
All Psa strains amplified at least one of the lsc genes, specifically lscC (18 strains, 54%) or both lscB and C (15 strains, 45%), which is consistent with levan production.By contrast, none of the Psa strains amplified the lscA gene.
Neither of the two strains of Pfm amplified any of the lsc genes even though they produced levan, so the existence of additional levansucrase isoenzymes cannot be ruled out, as has already been mentioned in previous work with other P. syringae strains [64].
Moreover, in the atypical strain series, we have found all the possible combinations, which means bacteria with and without lsc genes and producing or not producing levan.Six out of eighteen atypical strains (33%) amplified lscC, but only one of them, LPPA 2759, is also levan positive, so it is possible that the gene was not functional in the remaining five.Two strains (11%) amplified both lscB and C genes and were levan positive.A total of 50% of the strains in this group did not amplify the three screened genes, but LPPA 2758 was positive for levan production.In the latter case, and as previously said for Pfm, it is possible that isoenzymes other than those tested are involved [65].

Multilocus Sequencing Analysis (MLSA)
The majority of the strains identified as Psa had at least 99% similarity to the deposited Psa sequences when compared to the sequences of each gene by BLAST.
This was true for the gltA gene, except for strain LPPA 2980, which had a similarity percentage of 98.77%.With reference to the gyrB gene, the similarity percentage was >99%, except in the case of LPPA 2717 and LPPA 2718, for which the values were 98.40% and 98.45%, respectively.In contrast, the rpoD region of all the strains identified as Psa shared more than 99% sequence similarity with the deposited Psa sequences.
In the heterogeneous group of eighteen atypical strains, in the case of the gltA gene, for ten strains, the relatively closest matches were Psa and P. avellanae; in three strains, the closest matches were cf.Pfm and P. avellanae; in another three, it was cf.Pfm and in two others, it was P. syringae pv.tomato.With gyrB, eleven strains had P. avellanae as their closest relative, five had cf.Pfm and P. avellanae and two had P. syringae pv.tomato.Finally, in the case of rpoD, ten strains had P. avellanae as the closest match, two had cf.Pfm and P. avellanae, four had cf.Pfm and P. syringae pv.tomato, while two were closest to P. syringae pvs.antirrhini and tomato.

Phylogenetic Trees
Phylogenetic trees generated with the concatenated sequences of the gyrB, gltA and rpoD genes of the Psa and Pfm strains (Figure 1) involved 2263 bp (879 gltA, 603 gyrB and 781 rpoD).The analysis based on each of the genes separately (Supplementary Figure S1A-C) allows the differentiation of Psa and Pfm from the atypical strains.Estimates of the evolutionary divergence between the Psa and Pfm sequences are shown in Supplementary Table S2.
As shown in Figure 1, all isolates previously identified as Psa clustered with Psa strain CFBP 7181 with a 99% bootstrap.On the other hand, two strains of those considered atypical Psa, as they did not amplify the 230 bp band in the duplex PCR, clustered with CFBP 8161 corresponding to Pfm lineage 1.This pathovar was previously referred to as Psa biovar 4 and was subsequently considered a new pathovar [9].The fact that this pathovar only amplified the 492 bp band of the two expected for Psa in duplex PCR had been previously described by Cunty et al. [9], but in our case, 18 more strains showed this result, although they were not clustered with Pfm.This implies that both the simple and duplex PCR allow the identification of Psa but not of Pfm.
The numbers of segregating sites (S) and the nucleotide diversity (π) of the concatenated sequence of the three housekeeping genes (gltA, gyrB and rpoD) with respect to the three individual genes are shown in Table 3. GyrB was the gene that contributed most to nucleotide diversity, a result that concurs with that already obtained by Yin et al. [65].However, gltA was the gene that contributed the least to nucleotide diversity, both at the individual level and with the three concatenated genes, rpoD being the one that is closest to the value of π calculated for the three concatenated genes.The Psa strains analyzed seem to be a homogeneous group, as we had previously assumed since the plants analyzed in 2014 had only two origins: Galicia, which is the first Spanish geographic area where the disease appeared [18], and Italy, which is where the outbreak produced by Psa biovar 3 was described [13].
This result contrasts with the findings of other authors in other geographical areas, showing high genetic variability in Psa, whether related to the origin of the strains [66] or not [67].In Europe, it seems that the clonal expansion of Psa was followed by a broad genomic diversification, as reported by Figueira et al. [68].In Portugal, the presence of two genetically distinct subpopulations of Psa biovar 3 has been recently described [69,70].In the phylogenetic trees made with the 18 atypical strains, sequences of the pathovars that were found to be most closely related by BLAST were included.Figure 2 shows the phylogenetic tree based on concatenated partial sequences of the gltA, gyrB and rpoD genes, and the trees with the individual genes are included in Supplementary Figure S2A-C.Estimates of evolutionary divergence between atypical strain sequences are shown in Supplementary Table S3.
As can be seen in Figure 2, nine strains were grouped into a cluster with a bootstrap of 99%.Three strains isolated from kiwifruit samples collected in the same orchard were grouped with the IVIA 4447 strain that had been deposited as unconfirmed (cf) Pfm.We cannot know if the geographical origin of these strains is the same since Morán et al. [41] did not specify this data in their article, in which they only refer to "isolates from Asturias and other areas of Spain".In addition, we have not confirmed them as Pfm since they did not group with the Pfm control strain CFBP 8161.Two strains isolated from the same orchard and year, LPPA 2722 and LPPA 2759, were grouped with P. avellanae as the closest.
LPPA 2756 and LPPA 2757, also from the same orchard and year, were grouped with P. syringae pv.antirrhini as the closest.None of these strains have fully matched any of the nearby pathovars.S1.S1.S1.S1.
Concerning nucleotide diversity (Table 4), it can be observed that the gene that contributed the least to nucleotide diversity was gltA, while that which contributed the most was gyrB.However, rpoD was the gene that showed a value of nucleotide diversity closest to that obtained with the set of the three concatenated genes.This result matches with that obtained for Psa in the studied series (Table 3) and confirms that rpoD is a good phylogenetic marker, as described by other authors [71].

Pathogenicity Assays
Necrotic areas were observed in the pathogenicity tests carried out on yellow peppers (Capsicum annum cv 'California') inoculated with representative strains of Psa (LPPA 2727, LPPA 3700, NZ 10627E) and Pfm (LPPA 2769, CFBP 8039), without finding major differences between the Asturian strains; the Psa control strain presented the least browning.Control fruits showed no lesions (Figure 3).The inoculated bacteria were recovered and reidentified from the lesions produced.
Life 2024, 14, x FOR PEER REVIEW 12 o

Pathogenicity Assays
Necrotic areas were observed in the pathogenicity tests carried out on yellow pepp (Capsicum annum cv 'California') inoculated with representative strains of Psa (LPPA 27 LPPA 3700, NZ 10627E) and Pfm (LPPA 2769, CFBP 8039), without finding major diff ences between the Asturian strains; the Psa control strain presented the least browni Control fruits showed no lesions (Figure 3).The inoculated bacteria were recovered a reidentified from the lesions produced.When comparing lesions on peppers from Figures 3 and 4, slight differences in inoculation results are observed, with those in Figure 3 being more intense.This res suggests that it is possible that these atypical bacteria have low virulence and theref may cause mild damage to kiwifruit, so the economic impact would be milder on kiwifruit crop than that caused by Psa.
The important contribution that Abelleira et al. [40] made by testing different ho in their inoculations has allowed us to work with pepper, which has several advantag its lower acquisition and maintenance cost, the availability of the material to be inocula at any time of the year and the greater speed in obtaining the results.
Further work is needed to identify atypical strains by analyzing the whole genom On the other hand, it is necessary to clarify the role played by these atypical bacte in the symptoms observed in the orchards.Pseudomonas syringae is found in differ environments; it has been isolated from cultivated and wild plants, as well as in we [72][73][74][75][76], but also in the environment, e.g., in snow, clouds, soil, etc. [77][78][79].So, i necessary to clarify the role that P. syringae strains isolated from kiwifruit may have in kiwifruit pathosystem.From what has been indicated in this study, it seems that th strains could be of low virulence, at least in the cv.Hayward, but it would be advisable carry out inoculations on kiwifruit plants in order to follow the development of symptoms and determine if any of them may be of importance to kiwifruit culture a therefore, require some type of treatment.
However, it should be noted that the treatment of bacteriosis occurring in agricult is difficult.Nowadays, many studies are being carried out in order to find a treatment t can mitigate the damage caused by Psa, bearing in mind that the use of antibiotic restricted for health reasons as part of the EU's "One Health" strategy.Therefore, effo are directed towards the development of biological control agents (BCAs).In this conte phages that infect the bacterium have been described to control this disease [80][81][82][83][84][85].When comparing lesions on peppers from Figures 3 and 4, slight differences in the inoculation results are observed, with those in Figure 3 being more intense.This result suggests that it is possible that these atypical bacteria have low virulence and therefore may cause mild damage to kiwifruit, so the economic impact would be milder on the kiwifruit crop than that caused by Psa.
The important contribution that Abelleira et al. [40] made by testing different hosts in their inoculations has allowed us to work with pepper, which has several advantages: its lower acquisition and maintenance cost, the availability of the material to be inoculated at any time of the year and the greater speed in obtaining the results.
Further work is needed to identify atypical strains by analyzing the whole genome.On the other hand, it is necessary to clarify the role played by these atypical bacteria in the symptoms observed in the orchards.Pseudomonas syringae is found in different environments; it has been isolated from cultivated and wild plants, as well as in weeds [72][73][74][75][76], but also in the environment, e.g., in snow, clouds, soil, etc. [77][78][79].So, it is necessary to clarify the role that P. syringae strains isolated from kiwifruit may have in the kiwifruit pathosystem.From what has been indicated in this study, it seems that these strains could be of low virulence, at least in the cv.Hayward, but it would be advisable to carry out inoculations on kiwifruit plants in order to follow the development of the symptoms and determine if any of them may be of importance to kiwifruit culture and, therefore, require some type of treatment.
However, it should be noted that the treatment of bacteriosis occurring in agriculture is difficult.Nowadays, many studies are being carried out in order to find a treatment that can mitigate the damage caused by Psa, bearing in mind that the use of antibiotics is restricted for health reasons as part of the EU's "One Health" strategy.Therefore, efforts are directed towards the development of biological control agents (BCAs).In this context, phages that infect the bacterium have been described to control this disease [80][81][82][83][84][85].

1.
The presence in the Principality of Asturias of P. syringae pv.actinidiae has been confirmed, and P. s. pv.actinidifoliorum has been detected for the first time.

2.
The Psa strains under study have been shown to be a homogeneous group, which seems to indicate that there have been few introductions of the pathogen into the region.

3.
The simple + duplex PCR method used for the detection of Psa has allowed the correct identification of 33 Psa strains.However, it is not specific to Pfm, since of the twenty strains that initially gave the result described for this pathovar, only two were identified as Pfm. 4.
A total of 18 atypical strains were grouped in a phylogenetic tree with P. avellanae, P. syringae pv.antirrhini and a group of strains described as close to cf.Pfm.

5.
The pathogenicity tests carried out on pepper gave similar results for the atypical bacteria tested, so it will be necessary to carry the same tests out on kiwifruit plants to clarify their role in the kiwifruit pathosystem.

Life 2024 , 18 Figure 1 .
Figure 1.Phylogenetic tree of Psa and Pfm strains based on concatenated partial sequences of the gltA, gyrB and rpoD genes by using the maximum likelihood method and Tamura-Nei model.Bootstrap values ≥ 50% (1000 replicates) are indicated at branch points.Bar scale, substitutions per site.This analysis involved 38 nucleotide sequences, with a total of 2263 positions.Analyses were conducted in MEGA11 [60].P. asturiensis LPPA 221 T was used as outgroup and P. syringae pv.actinidiae CFBP 7811 and P. syringae pv.actinidifoliorum CFBP 8161 as controls.Accession numbers of the sequences are shown in Supplementary TableS1.

Figure 1 .
Figure 1.Phylogenetic tree of Psa and Pfm strains based on concatenated partial sequences of the gltA, gyrB and rpoD genes by using the maximum likelihood method and Tamura-Nei model.Bootstrap values ≥ 50% (1000 replicates) are indicated at branch points.Bar scale, substitutions per site.This analysis involved 38 nucleotide sequences, with a total of 2263 positions.Analyses were conducted in MEGA11 [60].P. asturiensis LPPA 221 T was used as outgroup and P. syringae pv.actinidiae CFBP 7811 and P. syringae pv.actinidifoliorum CFBP 8161 as controls.Accession numbers of the sequences are shown in Supplementary TableS1.

Figure 2 .
Figure 2. Phylogenetic tree of Ps strains based on concatenated partial sequences of the gltA, gyrB and rpoD genes by using the maximum likelihood method and Tamura-Nei model.Bootstrap values ≥ 50% (1000 replicates) are indicated at branch points.Bar scale, substitutions per site.This analysis involved 25 nucleotide sequences, with a total of 2092 positions.Analyses were conducted in MEGA11 [60].P. syringae pv.actinidiae CFBP 7811, P. syringae pv.actinidifoliorum CFBP 8161, P. syringae pv.tomato CFBP 2212, P. avellanae NCPPB 4222 and P. syringae pv.antirrhini strain 126 were included as controls and P. asturiensis 221 T as outgroup.Accession numbers of the sequences are shown in Supplementary TableS1.

Figure 2 .
Figure 2. Phylogenetic tree of Ps strains based on concatenated partial sequences of the gltA, gyrB and rpoD genes by using the maximum likelihood method and Tamura-Nei model.Bootstrap values ≥ 50% (1000 replicates) are indicated at branch points.Bar scale, substitutions per site.This analysis involved 25 nucleotide sequences, with a total of 2092 positions.Analyses were conducted in MEGA11 [60].P. syringae pv.actinidiae CFBP 7811, P. syringae pv.actinidifoliorum CFBP 8161, P. syringae pv.tomato CFBP 2212, P. avellanae NCPPB 4222 and P. syringae pv.antirrhini strain 126 were included as controls and P. asturiensis 221 T as outgroup.Accession numbers of the sequences are shown in Supplementary TableS1.

Table 1 .
Origin and features of Pseudomonas syringae isolates from the Principality of Asturias, from 2014 to 2020, that were included in this study.

Table 2 .
Biochemical features of the isolates under study.

Table 3 .
Segregating sites and nucleotide diversity of Psa + Pfm strains under study.
m = number of sequences, n = number of positions, S = number of segregating sites, π = nucleotide diversity.

Table 4 .
Segregating sites and nucleotide diversity of atypical P. syringae strains under study.
Funding:Institutional Review Board Statement: Not applicable.Informed Consent Statement: Not applicable.