Genome sequencing provides new insights on the distribution of Erwinia amylovora lineages in northern Italy

Intraspecific diversity of Erwinia amylovora strains from northern Algeria

BACKGROUND

Fire blight, caused by Erwinia amylovora, is the most destructive bacterial disease affecting plants in the Rosaceae family, leading to significant economic losses. In Algeria, this disease has been reported since 2010. This study aimed to investigate the origin of fire blight in Algeria, in order to increase knowledge of the epidemiology of this serious disease and contribute to its management. A comprehensive characterization of 18 E. amylovora isolates recovered from northern Algeria between 2016 and 2021 to evaluate their phenotypical and genotypical diversity was conducted.

RESULTS

Phenotypic differences, particularly in growth kinetics, virulence, and fatty acid profiles, allowed differentiation of strains into five groups, possibly indicating distinct introduction events. Genetic characterization revealed that only one strain lacked the ubiquitous plasmid pEA29, which is correlated with reduced virulence, while none harbored the pEI70 plasmid. Phylogenetic analysis using concatenated sequences of the recA, groEL, rpoS, ams, and hrpN genes grouped Algerian strains with those from a broadly prevalent clade. CRISPR genotyping identified a novel CR1 pattern and three genotypes, two of them previously unreported.

CONCLUSIONS

This study represents the first phenotypic, genetic, and phylogenetic investigation of E. amylovora strains in the region, and provides valuable information on the possible pathways of the introduction of this fire blight pathogen in northern Africa. The findings suggest one or more introduction events from a common ancestor, likely originating in northern Italy, followed by dispersal in various regions of Algeria.

Examination of Large Chromosomal Inversions in the Genome of Erwinia amylovora Strains Reveals Worldwide Distribution and North America-Specific Types

Erwinia amylovora is a relatively homogeneous species with low genetic diversity at the nucleotide level. However, phenotypic differences and genomic structural variations among E. amylovora strains have been documented. In this study, we identified 10 large chromosomal inversion (LCI) types in the Spiraeoideae-infecting (SI) E. amylovora strains by combining whole genome sequencing and PCR-based molecular markers. It was found that LCIs were mainly caused by homologous recombination events among seven rRNA operons (rrns) in SI E. amylovora strains. Although ribotyping results identified inter- and intra-variations in the internal transcribed spacer (ITS1 and ITS2) regions among rrns, LCIs tend to occur between rrns transcribed in the opposite directions and with the same tRNA content (tRNA-Glu or tRNA-Ile/Ala) in ITS1. Based on the LCI types, physical/estimated replichore imbalance (PRI/ERI) was examined and calculated. Among the 117 SI strains evaluated, the LCI types of Ea1189, CFBP1430, and Ea273 were the most common, with ERI values at 1.31, 7.87, and 4.47°, respectively. These three LCI types had worldwide distribution, whereas the remaining seven LCI types were restricted to North America (or certain regions of the United States). Our results indicated ongoing chromosomal recombination events in the SI E. amylovora population and showed that LCI events are mostly symmetrical, keeping the ERI less than 15°. These findings provide initial evidence about the prevalence of certain LCI types in E. amylovora strains, how LCI occurs, and its potential evolutionary advantage and history, which might help track the movement of the pathogen.

Novel approach toward the understanding of genetic diversity based on the two types of amino acid repeats in Erwinia amylovora

Erwinia amylovora is a notorious plant pathogenic bacterium of global concern that has devastated the apple and pear production industry worldwide. Nevertheless, the approaches available currently to understand the genetic diversity of E. amylovora remain unsatisfactory because of the lack of a trustworthy index and data covering the globally occurring E. amylovora strains; thus, their origin and distribution pattern remains ambiguous. Therefore, there is a growing need for robust approaches for obtaining this information via the comparison of the genomic structure of Amygdaloideae-infecting strains to understand their genetic diversity and distribution. Here, the whole-genome sequences of 245 E. amylovora strains available from the NCBI database were compared to identify intraspecific genes for use as an improved index for the simple classification of E. amylovora strains regarding their distribution. Finally, we discovered two kinds of strain-typing protein-encoding genes, i.e., the SAM-dependent methyltransferase and electron transport complex subunit RsxC. Interestingly, both of these proteins carried an amino acid repeat in these strains: SAM-dependent methyltransferase comprised a single-amino-acid repeat (asparagine), whereas RsxC carried a 40-amino-acid repeat, which was differentially distributed among the strains. These noteworthy findings and approaches may enable the exploration of the genetic diversity of E. amylovora from a global perspective.

Spontaneous Resistance of Erwinia amylovora Against Bacteriophage Y2 Affects Infectivity of Multiple Phages

Broad application of antibiotics gave rise to increasing numbers of antibiotic resistant bacteria. Therefore, effective alternatives are currently investigated. Bacteriophages, natural predators of bacteria, could work as such an alternative. Although phages can be highly effective at eliminating specific bacteria, phage resistance can be observed after application. The nature of this resistance, however, can differ depending on the phage. Exposing Erwinia amylovora CFBP 1430, the causative agent of fire blight, to the different phages Bue1, L1, S2, S6, or M7 led to transient resistance. The bacteria reversed to a phage sensitive state after the phage was eliminated. When wild type bacteria were incubated with Y2, permanently resistant colonies (1430^Y2R) formed spontaneously. In addition, 1430^Y2R revealed cross-resistance against other phages (Bue1) or lowered the efficiency of plating (L1, S2, and S6). Pull down experiments revealed that Y2 is no longer able to bind to the mutant suggesting mutation or masking of the Y2 receptor. Other phages tested were still able to bind to 1430^Y2R. Bue1 was observed to still adsorb to the mutant, but no host lysis was found. These findings indicated that, in addition to the alterations of the Y2 receptor, the 1430^Y2R mutant might block phage attack at different stage of infection. Whole genome sequencing of 1430^Y2R revealed a deletion in the gene with the locus tag EAMY_2231. The gene, which encodes a putative galactosyltransferase, was truncated due to the resulting frameshift. The mutant 1430^Y2R was monitored for potential defects or fitness loss. Weaker growth was observed in LB medium compared to the wild type but not in minimal medium. Strain 1430^Y2R was still highly virulent in blossoms even though amylovoran production was observed to be reduced. Additionally, LPS structures were analyzed and were clearly shown to be altered in the mutant. Complementation of the truncated EAMY_2231 in trans restored the wild type phenotype. The truncation of EAMY_2231 can therefore be associated with manifold modifications in 1430^Y2R, which can affect different phages simultaneously.