Comparative genomic analysis of Erwinia amylovora reveals novel insights in phylogenetic arrangement, plasmid diversity, and streptomycin resistance.

Erwinia amylovora is a destructive pathogen of Rosaceous plants and an economic concern worldwide. Herein, we report 93 new E. amylovora genomes from North America, Europe, the Mediterranean, and New Zealand. This new genomic information demonstrates the existence of three primary clades of Amygdaloideae (apple and pear) infecting E. amylovora and suggests all three independently originate from North America. The comprehensive sequencing also identified and confirmed the presence of 7 novel plasmids ranging in size from 2.9 to 34.7 kbp. While the function of the novel plasmids is unknown, the plasmids pEAR27, pEAR28, and pEAR35 encoded for type IV secretion systems. The strA-strB gene pair and the K43R point mutation at codon 43 of the rpsL gene have been previously documented to confer streptomycin resistance. Of the sequenced isolates, rpsL-based streptomycin resistance was more common and was found with the highest frequency in the Western North American clade.

Comparative genomics of Spiraeoideae-infecting Erwinia amylovora strains provides novel insight to genetic diversity and identifies the genetic basis of a low-virulence strain.

Erwinia amylovora is the causal agent of fire blight, one of the most devastating diseases of apple and pear. Erwinia amylovora is thought to have originated in North America and has now spread to at least 50 countries worldwide. An understanding of the diversity of the pathogen population and the transmission to different geographical regions is important for the future mitigation of this disease. In this research, we performed an expanded comparative genomic study of the Spiraeoideae-infecting (SI) E. amylovora population in North America and Europe. We discovered that, although still highly homogeneous, the genetic diversity of 30 E. amylovora genomes examined was about 30 times higher than previously determined. These isolates belong to four distinct clades, three of which display geographical clustering and one of which contains strains from various geographical locations ('Widely Prevalent' clade). Furthermore, we revealed that strains from the Widely Prevalent clade displayed a higher level of recombination with strains from a clade strictly from the eastern USA, which suggests that the Widely Prevalent clade probably originated from the eastern USA before it spread to other locations. Finally, we detected variations in virulence in the SI E. amylovora strains on immature pear, and identified the genetic basis of one of the low-virulence strains as being caused by a single nucleotide polymorphism in hfq, a gene encoding an important virulence regulator. Our results provide insights into the population structure, distribution and evolution of SI E. amylovora in North America and Europe.

Conservation of Erwinia amylovora pathogenicity-relevant genes among Erwinia genomes.

The Erwinia genus comprises species that are plant pathogens, non-pathogen, epiphytes, and opportunistic human pathogens. Within the genus, Erwinia amylovora ranks among the top 10 plant pathogenic bacteria. It causes the fire blight disease and is a global threat to commercial apple and pear production. We analyzed the presence/absence of the E. amylovora genes reported to be important for pathogenicity towards Rosaceae within various Erwinia strains genomes. This simple bottom-up approach, allowed us to correlate the analyzed genes to pathogenicity, host specificity, and make useful considerations to drive targeted studies.

Phylogeography and population structure of the biologically invasive phytopathogen Erwinia amylovora inferred using minisatellites.

Erwinia amylovora causes a major disease of pome fruit trees worldwide, and is regulated as a quarantine organism in many countries. While some diversity of isolates has been observed, molecular epidemiology of this bacterium is hindered by a lack of simple molecular typing techniques with sufficiently high resolution. We report a molecular typing system of E. amylovora based on variable number of tandem repeats (VNTR) analysis. Repeats in the E. amylovora genome were identified with comparative genomic tools, and VNTR markers were developed and validated. A Multiple-Locus VNTR Analysis (MLVA) was applied to E. amylovora isolates from bacterial collections representing global and regional distribution of the pathogen. Based on six repeats, MLVA allowed the distinction of 227 haplotypes among a collection of 833 isolates of worldwide origin. Three geographically separated groups were recognized among global isolates using Bayesian clustering methods. Analysis of regional outbreaks confirmed presence of diverse haplotypes but also high representation of certain haplotypes during outbreaks. MLVA analysis is a practical method for epidemiological studies of E. amylovora, identifying previously unresolved population structure within outbreaks. Knowledge of such structure can increase our understanding on how plant diseases emerge and spread over a given geographical region.

Molecular analyses of Erwinia amylovora strains isolated in Russia, Poland, Slovenia and Austria describing further spread of fire blight in Europe.

Fire blight, a bacteriosis of apple and pear, was assayed with molecular tools to associate its origin in Russia, Slovenia and south-eastern Austria with neighboring countries. The identification of all investigated strains was confirmed by MALDI-TOF mass spectroscopy except one. Independent isolation was verified by the level of amylovoran synthesis and by the number of short sequence DNA repeats in plasmid pEA29. DNA of gently lysed E. amylovora strains from Russia, Slovenia, Austria, Hungary, Italy, Spain, Croatia, Poland, Central Europe and Iran was treated with restriction enzymes XbaI and SpeI to create typical banding patterns for PFGE analysis. The pattern Pt2 indicated that most Russian E. amylovora strains were related to strains from Turkey and Iran. Strains from Slovenia exhibited patterns Pt3 and Pt2, both present in the neighboring countries. Strains were also probed for the recently described plasmid pEI70 detected in Pt1 strains from Poland and in Pt3 strains from other countries. The distribution of pattern Pt3 suggests distribution of fire blight from Belgium and the Netherlands to Central Spain and Northern Italy and then north to Carinthia. The PFGE patterns indicate that trade of plants may have introduced fire blight into southern parts of Europe proceeded by sequential spread.

Comparative genomics of 12 strains of Erwinia amylovora identifies a pan-genome with a large conserved core.

The plant pathogen Erwinia amylovora can be divided into two host-specific groupings; strains infecting a broad range of hosts within the Rosaceae subfamily Spiraeoideae (e.g., Malus, Pyrus, Crataegus, Sorbus) and strains infecting Rubus (raspberries and blackberries). Comparative genomic analysis of 12 strains representing distinct populations (e.g., geographic, temporal, host origin) of E. amylovora was used to describe the pan-genome of this major pathogen. The pan-genome contains 5751 coding sequences and is highly conserved relative to other phytopathogenic bacteria comprising on average 89% conserved, core genes. The chromosomes of Spiraeoideae-infecting strains were highly homogeneous, while greater genetic diversity was observed between Spiraeoideae- and Rubus-infecting strains (and among individual Rubus-infecting strains), the majority of which was attributed to variable genomic islands. Based on genomic distance scores and phylogenetic analysis, the Rubus-infecting strain ATCC BAA-2158 was genetically more closely related to the Spiraeoideae-infecting strains of E. amylovora than it was to the other Rubus-infecting strains. Analysis of the accessory genomes of Spiraeoideae- and Rubus-infecting strains has identified putative host-specific determinants including variation in the effector protein HopX1(Ea) and a putative secondary metabolite pathway only present in Rubus-infecting strains.

Comparative analysis of the Hrp pathogenicity island of Rubus- and Spiraeoideae-infecting Erwinia amylovora strains identifies the IT region as a remnant of an integrative conjugative element.

The Hrp pathogenicity island (hrpPAI) of Erwinia amylovora not only encodes a type III secretion system (T3SS) and other genes required for pathogenesis on host plants, but also includes the so-called island transfer (IT) region, a region that originates from an integrative conjugative element (ICE). Comparative genomic analysis of the IT regions of two Spiraeoideae- and three Rubus-infecting strains revealed that the regions in Spiraeoideae-infecting strains were syntenic and highly conserved in length and genetic information, but that the IT regions of the Rubus-infecting strains varied in gene content and length, showing a mosaic structure. None of the ICEs in E. amylovora strains were complete, as conserved ICE genes and the left border were missing, probably due to reductive genome evolution. Comparison of the hrpPAI region of E. amylovora strains to syntenic regions from other Erwinia spp. indicates that the hrpPAI and the IT regions are the result of several insertion and deletion events that have occurred within the ICE. It also suggests that the T3SS was present in a common ancestor of the pathoadapted Erwinia spp. and that insertion and deletion events in the IT region occurred during speciation.

Lipopolysaccharide biosynthesis genes discriminate between Rubus- and Spiraeoideae-infective genotypes of Erwinia amylovora.

Comparative genomic analysis revealed differences in the lipopolysaccharide (LPS) biosynthesis gene cluster between the Rubus-infecting strain ATCC BAA-2158 and the Spiraeoideae-infecting strain CFBP 1430 of Erwinia amylovora. These differences corroborate rpoB-based phylogenetic clustering of E. amylovora into four different groups and enable the discrimination of Spiraeoideae- and Rubus-infecting strains. The structure of the differences between the two groups supports the hypothesis that adaptation to Rubus spp. took place after species separation of E. amylovora and E. pyrifoliae that contrasts with a recently proposed scenario, based on CRISPR data, in which the shift to domesticated apple would have caused an evolutionary bottleneck in the Spiraeoideae-infecting strains of E. amylovora which would be a much earlier event. In the core region of the LPS biosynthetic gene cluster, Spiraeoideae-infecting strains encode three glycosyltransferases and an LPS ligase (Spiraeoideae-type waaL), whereas Rubus-infecting strains encode two glycosyltransferases and a different LPS ligase (Rubus-type waaL). These coding domains share little to no homology at the amino acid level between Rubus- and Spiraeoideae-infecting strains, and this genotypic difference was confirmed by polymerase chain reaction analysis of the associated DNA region in 31 Rubus- and Spiraeoideae-infecting strains. The LPS biosynthesis gene cluster may thus be used as a molecular marker to distinguish between Rubus- and Spiraeoideae-infecting strains of E. amylovora using primers designed in this study.

Eop1 from a Rubus strain of Erwinia amylovora functions as a host-range limiting factor.

Strains of Erwinia amylovora, the bacterium causing the disease fire blight of rosaceous plants, are separated into two groups based on host range: Spiraeoideae and Rubus strains. Spiraeoideae strains have wide host ranges, infecting plants in many rosaceous genera, including apple and pear. In the field, Rubus strains infect the genus Rubus exclusively, which includes raspberry and blackberry. Based on comparisons of limited sequence data from a Rubus and a Spiraeoideae strain, the gene eop1 was identified as unusually divergent, and it was selected as a possible host specificity factor. To test this, eop1 genes from a Rubus strain and a Spiraeoideae strain were cloned and mutated. Expression of the Rubus-strain eop1 reduced the virulence of E. amylovora in immature pear fruit and in apple shoots. Sequencing the orfA-eop1 regions of several strains of E. amylovora confirmed that forms of eop1 are conserved among strains with similar host ranges. This work provides evidence that eop1 from a Rubus-specific strain can function as a determinant of host specificity in E. amylovora.

Complete genome sequence of the fire blight pathogen Erwinia amylovora CFBP 1430 and comparison to other Erwinia spp.

Fire blight, caused by the enterobacterium Erwinia amylovora, is a devastating disease of rosaceous plants that has global economic importance for apple and pear production and trade. The complete genome of E. amylovora CFBP 1430 was sequenced, annotated, and compared with the genomes of other Erwinia spp. Several singleton and shared features of the E. amylovora CFBP 1430 genome were identified that offer a first view into evolutionary aspects within the genus Erwinia. Comparative genomics identified or clarified virulence and fitness determinants and secretion systems. Novel insights revealed in the genome of E. amylovora CFBP 1430 hold potential for exploitation to improve the design of more effective fire blight control strategies.

Erwinia amylovora: modern methods for detection and differentiation.

Erwinia amylovora is the causative agent of fire blight, a very destructive disease of numerous members of the rosaceae. The primary route of infection for host species, including commercially grown apple and pear, is the newly opened blossom. Susceptibility of flowers to infection for only a few days creates narrow window for infection. Not surprisingly, the risk of disease is related to E. amylovora population size. As a result, methods that supply quick, accurate and sensitive quantification of the pathogen population are important tools for determining the need for and the efficacy of disease control intervention. Plating samples and assessing colony-forming units constitutes an accurate and sensitive but slow method. Endpoint PCR is quick and sensitive but is not particularly amenable to quantification. We describe a real-time PCR procedure that provides all requirements. This method is based on chromosomal genes rather than on the pEa29 plasmid and so can be used to measure isolates that have been cured of the plasmid. The method has been used very successfully in directly quantify whole E. amylovora cells, in a variety of tissues from the orchard environment.

Polymerase chain reaction fingerprinting of Erwinia amylovora has a limited phylogenetic value but allows the design of highly specific molecular markers.

Erwinia amylovora, the causal agent of fire blight, is genetically very homogeneous, and current methodologies provide insufficient or contradictory information about the probable dispersal routes of the pathogen. With the final aim to obtain specific and reliable molecular markers for different lineages of the pathogen, we studied the molecular basis of rep-polymerase chain reaction (PCR) polymorphism using seven different arbitrary primers to fingerprint 93 E. amylovora strains from different countries, including Spain. Polymorphism was very low, and was displayed by only 11 E. amylovora strains, which produced 22 polymorphic bands. Five of 11 polymorphic bands cloned contained DNA that was present in more than 85% of the strains, whereas six bands were due to DNA present exclusively in the strains producing the rep-PCR polymorphism. Also, five of the polymorphic bands were due to the possession of either the ubiquitous plasmid pEA29, of plasmid pEU30, which was exclusively found in strains from North America, or of a 35-kb cryptic plasmid, present only in 28 strains from Northern Spain. We designed primer pairs from several cloned polymorphic bands that allowed the specific identification of the strains producing the polymorphism. Our results indicate that rep-PCR is not adequate for constructing genealogies of E. amylovora, although the strategy illustrated here, as well as the designed primers, can be used effectively in epidemiological studies with this pathogen.

Phenotypic diversity of Erwinia amylovora in Bulgaria.

Fifty-one strains of Erwinia amylovora isolated from nine host plants in Bulgaria were characterized phenotypically and identified by the API 20E and BIOLOG system. The identification was confirmed by PCR amplification of a specific region of the plasmid pEA29 and the genome ams region. The phenotypic diversity of the strains was studied on the basis of their API 20E and BIOLOG metabolic profiles, as well as of their SDS-PAGE protein profile. Metabolic diversity among the strains was established, but no connection with the origin of the strains was revealed. The Bulgarian strains showed API 20E metabolic profiles not found in previous studies of E. amylovora. The strains formed a homogenous group on the basis of their protein profiles. All the strains were sensitive to the antibiotics streptomycin, tetracycline and oxytetracycline. This study was an initial step towards an investigation of the diversity and evolution in the Bulgarian population of E. amylovora, and it was the first characterization of E. amylovora strains isolated from different host plants in the period 1995-2005 in Bulgaria.

Genetic differences between blight-causing Erwinia species with differing host specificities, identified by suppression subtractive hybridization.

PCR-based subtractive hybridization was used to isolate sequences from Erwinia amylovora strain Ea110, which is pathogenic on apples and pears, that were not present in three closely related strains with differing host specificities: E. amylovora MR1, which is pathogenic only on Rubus spp.; Erwinia pyrifoliae Ep1/96, the causal agent of shoot blight of Asian pears; and Erwinia sp. strain Ejp556, the causal agent of bacterial shoot blight of pear in Japan. In total, six subtractive libraries were constructed and analyzed. Recovered sequences included type III secretion components, hypothetical membrane proteins, and ATP-binding proteins. In addition, we identified an Ea110-specific sequence with homology to a type III secretion apparatus component of the insect endosymbiont Sodalis glossinidius, as well as an Ep1/96-specific sequence with homology to the Yersinia pestis effector protein tyrosine phosphatase YopH.

Characterization of Erwinia amylovora strains from different host plants using repetitive-sequences PCR analysis, and restriction fragment length polymorphism and short-sequence DNA repeats of plasmid pEA29.

AIMS: The three main aims of the study were the assessment of the genetic relationship between a deviating Erwinia amylovora strain isolated from Amelanchier sp. (Maloideae) grown in Canada and other strains from Maloideae and Rosoideae, the investigation of the variability of the PstI fragment of the pEA29 plasmid using restriction fragment length polymorphism (RFLP) analysis and the determination of the number of short-sequence DNA repeats (SSR) by DNA sequence analysis in representative strains. METHODS AND RESULTS: Ninety-three strains obtained from 12 plant genera and different geographical locations were examined by repetitive-sequences PCR using Enterobacterial Repetitive Intergenic Consensus, BOX and Repetitive Extragenic Palindromic primer sets. Upon the unweighted pair group method with arithmetic mean analysis, a deviating strain from Amelanchier sp. was analysed using amplified ribosomal DNA restriction analysis (ARDRA) analysis and the sequencing of the 16S rDNA gene. This strain showed 99% similarity to other E. amylovora strains in the 16S gene and the same banding pattern with ARDRA. The RFLP analysis of pEA29 plasmid using MspI and Sau3A restriction enzymes showed a higher variability than that previously observed and no clear-cut grouping of the strains was possible. The number of SSR units reiterated two to 12 times. The strains obtained from pear orchards showing for the first time symptoms of fire blight had a low number of SSR units. CONCLUSIONS: The strains from Maloideae exhibit a wider genetic variability than previously thought. The RFLP analysis of a fragment of the pEA29 plasmid would not seem a reliable method for typing E. amylovora strains. A low number of SSR units was observed with first epidemics of fire blight. SIGNIFICANCE AND IMPACT OF THE STUDY: The current detection techniques are mainly based on the genetic similarities observed within the strains from the cultivated tree-fruit crops. For a more reliable detection of the fire blight pathogen also in wild and ornamentals Rosaceous plants the genetic features of deviating E. amylovora strains have to be studied in detail.

New host plants of Erwinia amylovora in Bulgaria.

Nine strains of Erwinia amylovora were isolated from new host plants in Bulgaria--chokeberry and strawberry. The strains were characterized morphologically and biochemically using the API 20E and BIOLOG system. It was established that they showed three different API 20E metabolic profiles, not found by previous studies of E. amylovora. All strains were identified as E. amylovora due to their metabolic fingerprint patterns obtained by the BIOLOG system. The identification was confirmed by PCR amplification of a specific region of plasmid pEA29 and genome ams-region. This study is the first characterization and identification of E. amylovora strains isolated from chokeberry and strawberry by the API 20E and BIOLOG system and by polymerase chain reaction.

Stability of short sequence repeats and their application for the characterization of Erwinia amylovora strains.

A motif of eight nucleotides (GAATTACA) reiterated 3 to 15 times within the PstI fragment of the pEa29 plasmid was found in Erwinia amylovora strains representing a valuable typing method for this pathogen. The stability of short sequence DNA repeat (SSR) numbers was investigated to determine the suitability of this marker for strain differentiation. The number of SSR units was found to be stable under laboratory and certain stress conditions. This meets the requirements for a suitable genetic marker that should be stable upon cultivation of strains. Therefore, this SSR marker was used for strain differentiation from SSR-3 to SSR-15 and a large number of E. amylovora strains from Austria was screened for their SSR numbers for epidemiological identification purposes. Traceability was possible if strains had very high or very low SSR numbers.

Molecular differentiation of Erwinia amylovora strains from North America and of two Asian pear pathogens by analyses of PFGE patterns and hrpN genes.

In order to determine a possible genomic divergence of Erwinia amylovora'fruit tree' and raspberry strains from North America, several isolates were differentiated by pulsed-field gel electrophoresis (PFGE) analysis, the size of short DNA sequence repeats (SSRs) and the nucleotide and deduced amino acid sequences of their hrpN genes. By PFGE analysis European strains are highly related, whereas strains from North America were diverse and were further distinguished by the SSR numbers from plasmid pEA29. The E. amylovora strains from Europe showed identical HrpN sequences in contrast to the American isolates from fruit trees and raspberry. Those were related to each other, but distinguishable by their HrpN patterns. The Asian pear pathogens differed in HrpN among each other and from E. amylovora. Erwinia pyrifoliae isolates and the Erwinia strains from Japan were ordered via their HrpN sequences in agreement with the PFGE patterns. For all three pathogens, dendrograms from PFGE and sequence data indicate an evolutionary diversity within the species in spite of a genetic conservation for parts of the hrpN genes suggesting a long persistence of the Asian pear pathogens in Korea and Japan as well as of fire blight in North America. Some of the divergent American E. amylovora isolates share PFGE patterns with the relatively uniform European strains.