Complete genome sequence of the fire blight pathogen Erwinia pyrifoliae DSM 12163T and comparative genomic insights into plant pathogenicity

New Insights and Approach Toward the Genetic Diversity and Strain Typing of Erwinia pyrifoliae Based on rsxC, an Electron Transport Gene

Erwinia pyrifoliae, a causal agent of black shoot blight in apple and pear trees, is a plant pathogenic bacterium first reported in South Korea. The symptoms of black shoot blight are very similar to those of the fire blight disease in apple and pear trees caused by E. amylovora, as E. pyrifoliae has a genetically very close relationship with E. amylovora. Recently, there have been reports that E. pyrifoliae causes disease in European strawberries, resulting in severe fruit loss that aroused great concern about its spread, distribution, and host range. Therefore, it is essential to establish a trustworthy approach to understanding the distribution patterns of E. pyrifoliae based on the genetic background to strengthen the barrier of potential spreading risks, although advanced methods have been provided to accurately detect E. pyrifoliae and E. amylovora. Consequently, this study discovered a noble and noteworthy gene, rsxC, capable of providing the pathogen genotype by comparing E. pyrifoliae genomic sequences in the international representative genome archive. Different numbers of 40-unit amino acid repeats in this gene among the strains induced intraspecific traits in RsxC. By comparing their repeat pattern, E. pyrifoliae isolates were divided into two main groups, branching into several clades via sequence alignment of 35 E. pyrifoliae isolates from various apple orchards from 2020 to 2021 in South Korea. The newly discovered quadraginta amino acid repeat within this gene would be a valuable genetic touchstone for determining the genotype and distribution pattern of E. pyrifoliae strains, ultimately leading to exploring their evolution. The function of amino acid repeats and the biological significance of strains need to be elucidated further.

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.

Evidence of Greater Competitive Fitness of Erwinia amylovora over E. pyrifoliae in Korean Isolates

Erwinia amylovora and E. pyrifoliae are the causative agents of destructive diseases in both apple and pear trees viz. fire blight and black shoot blight, respectively. Since the introduction of fire blight in Korea in 2015, the occurrence of both pathogens has been independently reported. The co-incidence of these diseases is highly probable given the co-existence of their pathogenic bacteria in the same trees or orchards in a city/ district. Hence, this study evaluated whether both diseases occurred in neighboring orchards and whether they occurred together in a single orchard. The competition and virulence of the two pathogens was compared using growth rates in vitro and in planta. Importantly, E amylovora showed significantly higher colony numbers than E. pyrifoliae when they were co-cultured in liquid media and co-inoculated into immature apple fruits and seedlings. In a comparison of the usage of major carbon sources, which are abundant in immature apple fruits and seedlings, E. amylovora also showed better growth rates than E. pyrifoliae. In virulence assays, including motility and a hypersensitive response (HR), E. amylovora demonstrated a larger diameter of travel from the inoculation site than E. pyrifoliae in both swarming and swimming motilities. E. amylovora elicited a HR in tobacco leaves when diluted from 1:1 to 1:16 but E. pyrifoliae does not elicit a HR when diluted at 1:16. Therefore, E. amylovora was concluded to have a greater competitive fitness than E. pyrifoliae.

A Novel Signaling Pathway Connects Thiamine Biosynthesis, Bacterial Respiration, and Production of the Exopolysaccharide Amylovoran in Erwinia amylovora

Erwinia amylovora is a plant pathogen causing necrotrophic fire blight disease of apple, pear, and other rosaceous plants. This bacterium colonizes host vascular tissues via the production of exopolysaccharides (EPSs) including amylovoran. It is well-established that the nearly ubiquitous plasmid pEA29 of E. amylovora is an essential virulence factor, but the underlying mechanism remains uncharacterized. Here, we demonstrated that pEA29 was required for E. amylovora to produce amylovoran and to form a biofilm, and this regulation was dependent on the thiamine biosynthesis operon thiOSGF. We then conducted carbohydrate and genetic analyses demonstrating that the thiamine-mediated effect on amylovoran production was indirect, as cells lacking thiOSGF produced an EPS that did not contain glucuronic acid, one of the key components of amylovoran, whereas the transcriptional activity and RNA levels of the amylovoran biosynthesis genes were not altered. Alternatively, addition of exogenous thiamine restored amylovoran production in the pEA29-cured strain of E. amylovora and positively impacted amylovoran production in a dose-dependent manner. Individual deletion of several chromosomal thiamine biosynthesis genes also affected amylovoran production, implying that a complete thiamine biosynthesis pathway is required for the thiamine-mediated effect on amylovoran production in E. amylovora. Finally, we determined that an imbalanced tricarboxylic acid cycle negatively affected amylovoran production, which was restored by addition of exogenous thiamine or overexpression of the thiOSGF operon. In summary, our report revealed a novel signaling pathway that impacts E. amylovora virulence in which thiamine biosynthesis enhances bacterial respiration that provides energetic requirements for the biosynthesis of EPS amylovoran.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Virulence Genetics of an Erwinia amylovora Putative Polysaccharide Transporter Family Member

The Gram-negative enterobacterium Erwinia amylovora causes fire blight disease in apple and pear trees. Lipopolysaccharides and the exopolysaccharide amylovoran are essential E. amylovora virulence factors. We found that mutations in rfbX disrupted amylovoran production and virulence in apple fruits and tree shoots and that the deletion of yibD suppressed the rfbX mutant phenotype. The level of expression of yibD was about 10-fold higher in the ΔrfbX mutant than the wild type. A forward genetic suppressor screen in the ΔrfbX mutant uncovered multiple mutations in yibD and supported the conclusion that the virulence defect of rfbX mutants is due to reduced amylovoran production. The yibD and rfbX genes are expressed as a two-gene operon, yibD rfbX The rfbX gene encodes a previously uncharacterized putative polysaccharide subunit transporter, while yibD encodes a predicted glycosyltransferase. Mutation of rfbX did not have a detectable effect on lipopolysaccharide patterns; however, the overexpression of yibD in both the wild-type and ΔyibD ΔrfbX genetic backgrounds disrupted both amylovoran and lipopolysaccharide production. Additionally, the overexpression of yibD in the ΔyibD ΔrfbX mutant inhibited bacterial growth in amylovoran-inducing medium. This growth inhibition phenotype was used in a forward genetic suppressor screen and reverse-genetics tests to identify several genes involved in lipopolysaccharide production, which, when mutated, restored the ability of the ΔyibD ΔrfbX mutant overexpressing yibD to grow in amylovoran-inducing medium. Remarkably, all the lipopolysaccharide gene mutants tested were defective in lipopolysaccharide and amylovoran production. These results reveal a genetic connection between amylovoran and lipopolysaccharide production in E. amylovora IMPORTANCE This study discovered previously unknown genetic connections between exopolysaccharide and lipopolysaccharide production in the fire blight pathogen Erwinia amylovora This represents a step forward in our understanding of the biology underlying the production of these two macromolecules. Fire blight is an economically important disease that impacts the production of apples and pears worldwide. Few fire blight control measures are available, and growers rely heavily on antibiotic applications at bloom time. Both exopolysaccharide and lipopolysaccharide are E. amylovora virulence factors. Our results indicate that the overexpression of the yibD gene in E. amylovora disrupts both lipopolysaccharide production and exopolysaccharide production. This effect could potentially be used as the basis for the development of an antivirulence treatment for the prevention of fire blight disease.

Genomic divergence between Dickeya zeae strain EC2 isolated from rice and previously identified strains, suggests a different rice foot rot strain

Rice foot rot caused by Dickeya zeae is an important bacterial disease of rice worldwide. In this study, we identified a new strain EC2 from rice in Guangdong province, China. This strain differed from the previously identified strain from rice in its biochemical characteristics, pathogenicity, and genomic constituents. To explore genomic discrepancies between EC2 and previously identified strains from rice, a complete genome sequence of EC2 was obtained and used for comparative genomic analyses. The complete genome sequence of EC2 is 4,575,125 bp in length. EC2 was phylogenetically closest to previously identified Dickeya strains from rice, but not within their subgroup. In terms of secretion systems, genomic comparisons revealed that EC2 harbored only type I (T1SS), typeⅡ (T2SS), and type VI (T6SS) secretion systems. The flagella cluster of this strain possessed specific genomic characteristics like other D. zeae strains from Guangdong and from rice; within this locus, the genetic diversity among strains from rice was much lower than that of within strains from non-rice hosts. Unlike other strains from rice, EC2 lost the zeamine cluster, but retained the clustered regularly interspaced short palindromic repeats-1 (CRISPR-1) array. Compared to the other D. zeae strains containing both exopolysaccharide (EPS) and capsular polysaccharide (CPS) clusters, EC2 harbored only the CPS cluster, while the other strains from rice carried only the EPS cluster. Furthermore, we found strain MS1 from banana, carrying both EPS and CPS clusters, produced significantly more EPS than the strains from rice, and exhibited different biofilm-associated phenotypes. Comparative genomics analyses suggest EC2 likely evolved through a pathway different from the other D. zeae strains from rice, producing a new type of rice foot rot pathogen. These findings emphasize the emergence of a new type of D. zeae strain causing rice foot rot, an essential step in the early prevention of this rice bacterial disease.

Comparative Genome Analysis Reveals Natural Variations in the Genomes of Erwinia pyrifoliae, a Black Shoot Blight Pathogen in Apple and Pear

Erwinia pyrifoliae is a Gram-negative bacterial plant pathogen that causes black shoot blight in apple and pear. Although earlier studies reported the genome comparison of Erwinia species, E. pyrifoliae strains for such analysis were isolated in 1996. In 2014, the strain E. pyrifoliae EpK1/15 was newly isolated in the apple tree showing black shoot blight in South Korea. This study aimed to better understand the similarities and differences caused by natural variations at the genomic level between newly isolated E. pyrifoliae EpK1/15 and the strain Ep1/96, which were isolated almost 20 years apart. Several comparative genomic analyses were conducted, and Clusters of Orthologous Groups of proteins (COG) database was used to classify functional annotation for each strain. E. pyrifoliae EpK1/15 had similarities with the Ep1/96 strain in stress-related genes, Tn3 transposase of insertion sequences, type III secretion systems, and small RNAs. The most remarkable difference to emerge from this comparison was that although the draft genome of E. pyrifoliae EpK1/15 was almost conserved, Epk1/15 strain had at least three sorts of structural variations in functional annotation according to COG database; chromosome inversion, translocation, and duplication. These results indicate that E. pyrifoliae species has gone natural variations within almost 20 years at the genomic level, and we can trace their similarities and differences with comparative genomic analysis.

ClpXP-Dependent RpoS Degradation Enables Full Activation of Type III Secretion System, Amylovoran Production, and Motility in Erwinia amylovora

Erwinia amylovora, the causal agent of fire blight disease of apple and pear, employs intracellular proteases, including Lon and ClpXP, for posttranslational regulation of various cellular proteins. It has been shown that Lon plays a critical role in E. amylovora virulence by directly targeting type III secretion system (T3SS) proteins and the Rcs phosphorelay system. In this study, we genetically examined the role of ClpXP and its potential interaction with Lon in E. amylovora. Mutation in clpXP diminished the expression of the T3SS, reduced exopolysaccharide amylovoran production and motility, and resulted in delayed disease progress. Western blot analyses showed highly accumulated RpoS proteins in the clpXP mutant. Moreover, mutation of rpoS in the clpXP mutant background rescued the expression of the T3SS and amylovoran production, suggesting that ClpXP-dependent RpoS degradation positively affects virulence traits. Interestingly, lack of both ClpXP and Lon resulted in significantly reduced virulence but increased expression of the T3SS and amylovoran production. However, this phenomenon was independent of RpoS accumulation, suggesting that ClpXP and Lon are indispensable for full virulence in E. amylovora.

Draft Genome Sequence of Pantoea agglomerans JM1, a Strain Isolated from Soil Polluted by Industrial Production of Beta-Lactam Antibiotics That Exhibits Valacyclovir-Like Hydrolase Activity

Strain Pantoea agglomerans JM1 was isolated from the soil of a microbiome that had been exposed to polluting pharmaceuticals. The bacterium exhibited enzymatic activities useful for the biotransformation of beta-lactams. The genome of the strain was assembled and described, and the gene encoding valacyclovir-like hydrolase was identified.

Role of the type VI secretion systems during disease interactions of Erwinia amylovora with its plant host

Background

Type VI secretion systems (T6SS) are widespread among Gram-negative bacteria and have a potential role as essential virulence factors or to maintain symbiotic interactions. Three T6SS gene clusters were identified in the genome of E. amylovora CFBP 1430, of which T6SS-1 and T6SS-3 represent complete T6SS machineries, while T6SS-2 is reduced in its gene content.

Results

To assess the contribution of T6SSs to virulence and potential transcriptomic changes of E. amylovora CFBP 1430, single and double mutants in two structural genes were generated for T6SS-1 and T6SS-3. Plant assays showed that mutants in T6SS-3 were slightly more virulent in apple shoots while inducing less disease symptoms on apple flowers, indicating that T6SSs have only a minor effect on virulence of E. amylovora CFBP 1430. The mutations led under in vitro conditions to the differential expression of type III secretion systems, iron acquisition, chemotaxis, flagellar, and fimbrial genes. Comparison of the in planta and in vitro transcriptome data sets revealed a common differential expression of three processes and a set of chemotaxis and motility genes. Additional experiments proved that T6SS mutants are impaired in their motility.

Conclusion

These results suggest that the deletion of T6SSs alters metabolic and motility processes. Nevertheless, the difference in lesion development in apple shoots and flower necrosis of T6SS mutants was indicative that T6SSs influences the disease progression and the establishment of the pathogen on host plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4010-1) contains supplementary material, which is available to authorized users.

Recurrent evolution of gut symbiotic bacteria in pentatomid stinkbugs

BACKGROUND

Diverse animals are intimately associated with microbial symbionts. How such host-symbiont associations have evolved is a fundamental biological issue. Recent studies have revealed a variety of evolutionary relationships, such as obligatory, facultative, and free-living, of gut bacterial symbiosis within the stinkbug family Pentatomidae, although the whole evolutionary picture remains elusive.

RESULTS

Here we investigated a comprehensive assembly of Japanese pentatomid stinkbugs representing 28 genera, 35 species, and 143 populations. Polymerase chain reaction (PCR), cloning, and sequencing of bacterial 16S rRNA gene from their midgut symbiotic organ consistently detected a single bacterial species from each of the insect samples, indicating a general tendency toward monosymbiotic gut association. Bacterial sequences detected from different populations of the same species were completely or nearly identical, indicating that the majority of the gut symbiotic associations are stably maintained at the species level. Furthermore, bacterial sequences detected from different species in the same genus tended to form well-supported clades, suggesting that host-symbiont associations are often stable even at the genus level. Meanwhile, when we compared such sequences with published sequences available in DNA databases, we found a number of counter-examples to such stable host-symbiont relationships; i.e., symbionts from different host species in the same genus may be phylogenetically distant, and symbionts from the same host species may be phylogenetically diverse. Likewise, symbionts of diverse pentatomid species may be closely related to symbionts of other stinkbug families, and symbionts of diverse pentatomid species may even be allied to free-living bacteria. Molecular evolutionary analyses revealed that higher molecular evolutionary rates, higher AT nucleotide compositions, and smaller genome sizes tended to be associated with the pentatomid symbionts constituting the stable lineages, whereas these traits were rarely observed in the pentatomid symbionts of promiscuous type.

CONCLUSIONS

These results indicate that gut symbiotic bacteria have evolved repeatedly and dynamically in the stinkbug family Pentatomidae, which have plausibly entailed frequent symbiont acquisitions, losses, replacements and transfers, while establishing a number of relatively stable host-symbiont associations. The diverse host-symbiont relationships observed in the Pentatomidae will provide an ideal arena for investigating the evolution of symbiosis experimentally and theoretically.

Horizontal Gene Acquisitions, Mobile Element Proliferation, and Genome Decay in the Host-Restricted Plant Pathogen Erwinia Tracheiphila

Modern industrial agriculture depends on high-density cultivation of genetically similar crop plants, creating favorable conditions for the emergence of novel pathogens with increased fitness in managed compared with ecologically intact settings. Here, we present the genome sequence of six strains of the cucurbit bacterial wilt pathogen Erwinia tracheiphila (Enterobacteriaceae) isolated from infected squash plants in New York, Pennsylvania, Kentucky, and Michigan. These genomes exhibit a high proportion of recent horizontal gene acquisitions, invasion and remarkable amplification of mobile genetic elements, and pseudogenization of approximately 20% of the coding sequences. These genome attributes indicate that E. tracheiphila recently emerged as a host-restricted pathogen. Furthermore, chromosomal rearrangements associated with phage and transposable element proliferation contribute to substantial differences in gene content and genetic architecture between the six E. tracheiphila strains and other Erwinia species. Together, these data lead us to hypothesize that E. tracheiphila has undergone recent evolution through both genome decay (pseudogenization) and genome expansion (horizontal gene transfer and mobile element amplification). Despite evidence of dramatic genomic changes, the six strains are genetically monomorphic, suggesting a recent population bottleneck and emergence into E. tracheiphila's current ecological niche.

Comparative genome sequencing to assess the genetic diversity and virulence attributes of 15 Vibrio anguillarum isolates

Vibrio anguillarum is the causative agent of vibriosis, a deadly haemorrhagic septicaemic disease affecting various marine and fresh/brackish water fish, bivalves and crustaceans. However, the diversity and virulence mechanisms of this pathogen are still insufficiently known. In this study, we aimed to increase our understanding of V. anguillarum diversity and virulence through comparative genome analysis of 15 V. anguillarum strains, obtained from different hosts or non-host niches and geographical regions, among which 10 and 5 strains were found to be virulent and avirulent, respectively, against sea bass larvae. First, the 15 draft genomes were annotated and screened for putative virulence factors, including genes encoding iron uptake systems, transport systems and non-ribosomal peptide synthetases. Second, comparative genome analysis was performed, focusing on single nucleotide polymorphisms (SNPs) and small insertions and deletions (InDels). Five V. anguillarum strains showed a remarkably high nucleotide identity. However, these strains comprise both virulent and avirulent strains towards sea bass larvae, suggesting that differences in virulence may be caused by subtle nucleotide variations. Clearly, the draft genome sequence of these 15 strains represents a starting point for further genetic research of this economically important fish pathogen.

Genetic islands in pome fruit pathogenic and non-pathogenic Erwinia species and related plasmids

New pathogenic bacteria belonging to the genus Erwinia associated with pome fruit trees (Erwinia, E. piriflorinigrans, E. uzenensis) have been increasingly described in the last years, and comparative analyses have found that all these species share several genetic characteristics. Studies at different level (whole genome comparison, virulence genes, plasmid content, etc.) show a high intraspecies homogeneity (i.e., among E. amylovora strains) and also abundant similarities appear between the different Erwinia species: presence of plasmids of similar size in the pathogenic species; high similarity in several genes associated with exopolysaccharide production and hence, with virulence, as well as in some other genes, in the chromosomes. Many genetic similarities have been observed also among some of the plasmids (and genomes) from the pathogenic species and E. tasmaniensis or E. billingiae, two epiphytic species on the same hosts. The amount of genetic material shared in this genus varies from individual genes to clusters, genomic islands and genetic material that even may constitute a whole plasmid. Recent research on evolution of erwinias point out the horizontal transfer acquisition of some genomic islands that were subsequently lost in some species and several pathogenic traits that are still present. How this common material has been obtained and is efficiently maintained in different species belonging to the same genus sharing a common ecological niche provides an idea of the origin and evolution of the pathogenic Erwinia and the interaction with non-pathogenic species present in the same niche, and the role of the genes that are conserved in all of them.

Draft Genomes, Phylogenetic Reconstruction, and Comparative Genomics of Two Novel Cohabiting Bacterial Symbionts Isolated from Frankliniella occidentalis

Obligate bacterial symbionts are widespread in many invertebrates, where they are often confined to specialized host cells and are transmitted directly from mother to progeny. Increasing numbers of these bacteria are being characterized but questions remain about their population structure and evolution. Here we take a comparative genomics approach to investigate two prominent bacterial symbionts (BFo1 and BFo2) isolated from geographically separated populations of western flower thrips, Frankliniella occidentalis. Our multifaceted approach to classifying these symbionts includes concatenated multilocus sequence analysis (MLSA) phylogenies, ribosomal multilocus sequence typing (rMLST), construction of whole-genome phylogenies, and in-depth genomic comparisons. We showed that the BFo1 genome clusters more closely to species in the genus Erwinia, and is a putative close relative to Erwinia aphidicola. BFo1 is also likely to have shared a common ancestor with Erwinia pyrifoliae/Erwinia amylovora and the nonpathogenic Erwinia tasmaniensis and genetic traits similar to Erwinia billingiae. The BFo1 genome contained virulence factors found in the genus Erwinia but represented a divergent lineage. In contrast, we showed that BFo2 belongs within the Enterobacteriales but does not group closely with any currently known bacterial species. Concatenated MLSA phylogenies indicate that it may have shared a common ancestor to the Erwinia and Pantoea genera, and based on the clustering of rMLST genes, it was most closely related to Pantoea ananatis but represented a divergent lineage. We reconstructed a core genome of a putative common ancestor of Erwinia and Pantoea and compared this with the genomes of BFo bacteria. BFo2 possessed none of the virulence determinants that were omnipresent in the Erwinia and Pantoea genera. Taken together, these data are consistent with BFo2 representing a highly novel species that maybe related to known Pantoea.

Commonalities and differences of T3SSs in rhizobia and plant pathogenic bacteria

Plant pathogenic bacteria and rhizobia infect higher plants albeit the interactions with their hosts are principally distinct and lead to completely different phenotypic outcomes, either pathogenic or mutualistic, respectively. Bacterial protein delivery to plant host plays an essential role in determining the phenotypic outcome of plant-bacteria interactions. The involvement of type III secretion systems (T3SSs) in mediating animal- and plant-pathogen interactions was discovered in the mid-80's and is now recognized as a multiprotein nanomachine dedicated to trans-kingdom movement of effector proteins. The discovery of T3SS in bacteria with symbiotic lifestyles broadened its role beyond virulence. In most T3SS-positive bacterial pathogens, virulence is largely dependent on functional T3SSs, while in rhizobia the system is dispensable for nodulation and can affect positively or negatively the mutualistic associations with their hosts. This review focuses on recent comparative genome analyses in plant pathogens and rhizobia that uncovered similarities and variations among T3SSs in their genetic organization, regulatory networks and type III secreted proteins and discusses the evolutionary adaptations of T3SSs and type III secreted proteins that might account for the distinguishable phenotypes and host range characteristics of plant pathogens and symbionts.

PathogenFinder--distinguishing friend from foe using bacterial whole genome sequence data

Although the majority of bacteria are harmless or even beneficial to their host, others are highly virulent and can cause serious diseases, and even death. Due to the constantly decreasing cost of high-throughput sequencing there are now many completely sequenced genomes available from both human pathogenic and innocuous strains. The data can be used to identify gene families that correlate with pathogenicity and to develop tools to predict the pathogenicity of newly sequenced strains, investigations that previously were mainly done by means of more expensive and time consuming experimental approaches. We describe PathogenFinder (http://cge.cbs.dtu.dk/services/PathogenFinder/), a web-server for the prediction of bacterial pathogenicity by analysing the input proteome, genome, or raw reads provided by the user. The method relies on groups of proteins, created without regard to their annotated function or known involvement in pathogenicity. The method has been built to work with all taxonomic groups of bacteria and using the entire training-set, achieved an accuracy of 88.6% on an independent test-set, by correctly classifying 398 out of 449 completely sequenced bacteria. The approach here proposed is not biased on sets of genes known to be associated with pathogenicity, thus the approach could aid the discovery of novel pathogenicity factors. Furthermore the pathogenicity prediction web-server could be used to isolate the potential pathogenic features of both known and unknown strains.

Comparative genome analysis of Enterobacter cloacae

The Enterobacter cloacae species includes an extremely diverse group of bacteria that are associated with plants, soil and humans. Publication of the complete genome sequence of the plant growth-promoting endophytic E. cloacae subsp. cloacae ENHKU01 provided an opportunity to perform the first comparative genome analysis between strains of this dynamic species. Examination of the pan-genome of E. cloacae showed that the conserved core genome retains the general physiological and survival genes of the species, while genomic factors in plasmids and variable regions determine the virulence of the human pathogenic E. cloacae strain; additionally, the diversity of fimbriae contributes to variation in colonization and host determination of different E. cloacae strains. Comparative genome analysis further illustrated that E. cloacae strains possess multiple mechanisms for antagonistic action against other microorganisms, which involve the production of siderophores and various antimicrobial compounds, such as bacteriocins, chitinases and antibiotic resistance proteins. The presence of Type VI secretion systems is expected to provide further fitness advantages for E. cloacae in microbial competition, thus allowing it to survive in different environments. Competition assays were performed to support our observations in genomic analysis, where E. cloacae subsp. cloacae ENHKU01 demonstrated antagonistic activities against a wide range of plant pathogenic fungal and bacterial species.

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.

Genome sequence of Pectobacterium carotovorum subsp. carotovorum strain PCC21, a pathogen causing soft rot in Chinese cabbage

Pectobacterium carotovorum is a plant-pathogenic enterobacterium responsible for soft rot in various commercially important plants. Here we report the complete genome sequence and automatic annotation of strain PCC21.

Erwinia amylovora loop-mediated isothermal amplification (LAMP) assay for rapid pathogen detection and on-site diagnosis of fire blight

Several molecular methods have been developed for the detection of Erwinia amylovora, the causal agent of fire blight in pear and apple, but none are truly applicable for on-site use in the field. We developed a fast, reliable and field applicable detection method using a novel target on the E. amylovora chromosome that we identified by applying a comparative genomic pipeline. The target coding sequences (CDSs) are both uniquely specific for and all-inclusive of E. amylovora genotypes. This avoids potential false negatives that can occur with most commonly used methods based on amplification of plasmid gene targets, which can vary among strains. Loop-mediated isothermal AMPlification (LAMP) with OptiGene Genie II chemistry and instrumentation proved to be an exceptionally rapid (under 15 min) and robust method for detecting E. amylovora in orchards, as well as simple to use in the plant diagnostic laboratory. Comparative validation results using plant samples from inoculated greenhouse trials and from natural field infections (of regional and temporal diverse origin) showed that our LAMP had an equivalent or greater performance regarding sensitivity, specificity, speed and simplicity than real-time PCR (TaqMan), other LAMP assays, immunoassays and plating, demonstrating its utility for routine testing.

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.

Erwinia amylovora CRISPR elements provide new tools for evaluating strain diversity and for microbial source tracking

Clustered regularly interspaced short palindromic repeats (CRISPRs) comprise a family of short DNA repeat sequences that are separated by non repetitive spacer sequences and, in combination with a suite of Cas proteins, are thought to function as an adaptive immune system against invading DNA. The number of CRISPR arrays in a bacterial chromosome is variable, and the content of each array can differ in both repeat number and in the presence or absence of specific spacers. We utilized a comparative sequence analysis of CRISPR arrays of the plant pathogen Erwinia amylovora to uncover previously unknown genetic diversity in this species. A total of 85 E. amylovora strains varying in geographic isolation (North America, Europe, New Zealand, and the Middle East), host range, plasmid content, and streptomycin sensitivity/resistance were evaluated for CRISPR array number and spacer variability. From these strains, 588 unique spacers were identified in the three CRISPR arrays present in E. amylovora, and these arrays could be categorized into 20, 17, and 2 patterns types, respectively. Analysis of the relatedness of spacer content differentiated most apple and pear strains isolated in the eastern U.S. from western U.S. strains. In addition, we identified North American strains that shared CRISPR genotypes with strains isolated on other continents. E. amylovora strains from Rubus and Indian hawthorn contained mostly unique spacers compared to apple and pear strains, while strains from loquat shared 79% of spacers with apple and pear strains. Approximately 23% of the spacers matched known sequences, with 16% targeting plasmids and 5% targeting bacteriophage. The plasmid pEU30, isolated in E. amylovora strains from the western U.S., was targeted by 55 spacers. Lastly, we used spacer patterns and content to determine that streptomycin-resistant strains of E. amylovora from Michigan were low in diversity and matched corresponding streptomycin-sensitive strains from the background population.

Detection of AI-2 receptors in genomes of Enterobacteriaceae suggests a role of type-2 quorum sensing in closed ecosystems

The LuxS enzyme, an S-ribosyl-homocysteine lyase, catalyzes the production of the signal precursor for autoinducer-2 mediated quorum sensing (QS-2) in Vibrio. Its widespread occurrence among bacteria is often considered the evidence for a universal language for interspecies communication. Presence of the luxS gene and production of the autoinducer-2 (AI-2) signal have repeatedly been the only evidences presented to assign a functional QS-2 to the most diverse species. In fact, LuxS has a primary metabolic role as part of the activated methyl cycle. In this review we have analyzed the distribution of QS-2 related genes in Enterobacteriaceae by moving the focus of the investigation from AI-2 production to the detection of potential AI-2 receptors. The latter are common in pathogens or endosymbionts of animals, but were also found in a limited number of Enterobacteriaceae of the genera Enterobacter, Klebsiella, and Pantoea that live in close association with plants or fungi. Although a precise function of QS-2 in these species has not been identified, they all show an endophytic or endosymbiontic lifestyle that suggests a role of type-2 quorum sensing in the adaptation to closed ecosystems.

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.

Fire blight: applied genomic insights of the pathogen and host

The enterobacterial phytopathogen Erwinia amylovora causes fire blight, an invasive disease that threatens a wide range of commercial and ornamental Rosaceae host plants. The response elicited by E. amylovora in its host during disease development is similar to the hypersensitive reaction that typically leads to resistance in an incompatible host-pathogen interaction, yet no gene-for-gene resistance has been described for this host-pathogen system. Comparative genomic analysis has found an unprecedented degree of genetic uniformity among strains of E. amylovora, suggesting that the pathogen has undergone a recent genetic bottleneck. The genome of apple, an important host of E. amylovora, has been sequenced, creating new opportunities for the study of interactions between host and pathogen during fire blight development and for the identification of resistance genes. This review includes recent advances in the genomics of both host and pathogen.

Comparative genomics of the Type VI secretion systems of Pantoea and Erwinia species reveals the presence of putative effector islands that may be translocated by the VgrG and Hcp proteins

Background

The Type VI secretion apparatus is assembled by a conserved set of proteins encoded within a distinct locus. The putative effector proteins Hcp and VgrG are also encoded within these loci. We have identified numerous distinct Type VI secretion system (T6SS) loci in the genomes of several ecologically diverse Pantoea and Erwinia species and detected the presence of putative effector islands associated with the hcp and vgrG genes.

Results

Between two and four T6SS loci occur among the Pantoea and Erwinia species. While two of the loci (T6SS-1 and T6SS-2) are well conserved among the various strains, the third (T6SS-3) locus is not universally distributed. Additional orthologous loci are present in Pantoea sp. aB-valens and Erwinia billingiae Eb661. Comparative analysis of the T6SS-1 and T6SS-3 loci showed non-conserved islands associated with the vgrG and hcp, and vgrG genes, respectively. These regions had a G+C content far lower than the conserved portions of the loci. Many of the proteins encoded within the hcp and vgrG islands carry conserved domains, which suggests they may serve as effector proteins for the T6SS. A number of the proteins also show homology to the C-terminal extensions of evolved VgrG proteins.

Conclusions

Extensive diversity was observed in the number and content of the T6SS loci among the Pantoea and Erwinia species. Genomic islands could be observed within some of T6SS loci, which are associated with the hcp and vgrG proteins and carry putative effector domain proteins. We propose new hypotheses concerning a role for these islands in the acquisition of T6SS effectors and the development of novel evolved VgrG and Hcp proteins.

Comparative Genomics of Erwinia amylovora and Related Erwinia Species-What do We Learn?

Erwinia amylovora, the causal agent of fire blight disease of apples and pears, is one of the most important plant bacterial pathogens with worldwide economic significance. Recent reports on the complete or draft genome sequences of four species in the genus Erwinia, including E. amylovora, E. pyrifoliae, E. tasmaniensis, and E. billingiae, have provided us near complete genetic information about this pathogen and its closely-related species. This review describes in silico subtractive hybridization-based comparative genomic analyses of eight genomes currently available, and highlights what we have learned from these comparative analyses, as well as genetic and functional genomic studies. Sequence analyses reinforce the assumption that E. amylovora is a relatively homogeneous species and support the current classification scheme of E. amylovora and its related species. The potential evolutionary origin of these Erwinia species is also proposed. The current understanding of the pathogen, its virulence mechanism and host specificity from genome sequencing data is summarized. Future research directions are also suggested.

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.

Metabolic versatility and antibacterial metabolite biosynthesis are distinguishing genomic features of the fire blight antagonist Pantoea vagans C9-1

BACKGROUND

Pantoea vagans is a commercialized biological control agent used against the pome fruit bacterial disease fire blight, caused by Erwinia amylovora. Compared to other biocontrol agents, relatively little is currently known regarding Pantoea genetics. Better understanding of antagonist mechanisms of action and ecological fitness is critical to improving efficacy.

PRINCIPAL FINDINGS

Genome analysis indicated two major factors Contribute to biocontrol activity: competition for limiting substrates and antibacterial metabolite production. Pathways for utilization of a broad diversity of sugars and acquisition of iron were identified. Metabolism of sorbitol by P. vagans C9-1 may be a major metabolic feature in biocontrol of fire blight. Biosynthetic genes for the antibacterial peptide pantocin A were found on a chromosomal 28-kb genomic island, and for dapdiamide E on the plasmid pPag2. There was no evidence of potential virulence factors that could enable an animal or phytopathogenic lifestyle and no indication of any genetic-based biosafety risk in the antagonist.

CONCLUSIONS

Identifying key determinants contributing to disease suppression allows the development of procedures to follow their expression in planta and the genome sequence contributes to rationale risk assessment regarding the use of the biocontrol strain in agricultural systems.

Diversity, evolution, and functionality of clustered regularly interspaced short palindromic repeat (CRISPR) regions in the fire blight pathogen Erwinia amylovora

The clustered regularly interspaced short palindromic repeat (CRISPR)/Cas system confers acquired heritable immunity against mobile nucleic acid elements in prokaryotes, limiting phage infection and horizontal gene transfer of plasmids. In CRISPR arrays, characteristic repeats are interspersed with similarly sized nonrepetitive spacers derived from transmissible genetic elements and acquired when the cell is challenged with foreign DNA. New spacers are added sequentially and the number and type of CRISPR units can differ among strains, providing a record of phage/plasmid exposure within a species and giving a valuable typing tool. The aim of this work was to investigate CRISPR diversity in the highly homogeneous species Erwinia amylovora, the causal agent of fire blight. A total of 18 CRISPR genotypes were defined within a collection of 37 cosmopolitan strains. Strains from Spiraeoideae plants clustered in three major groups: groups II and III were composed exclusively of bacteria originating from the United States, whereas group I generally contained strains of more recent dissemination obtained in Europe, New Zealand, and the Middle East. Strains from Rosoideae and Indian hawthorn (Rhaphiolepis indica) clustered separately and displayed a higher intrinsic diversity than that of isolates from Spiraeoideae plants. Reciprocal exclusion was generally observed between plasmid content and cognate spacer sequences, supporting the role of the CRISPR/Cas system in protecting against foreign DNA elements. However, in several group III strains, retention of plasmid pEU30 is inconsistent with a functional CRISPR/Cas system.

Gnom(Cmp): a quantitative approach for comparative analysis of closely related genomes of bacterial pathogens

Comparative genome analysis is a powerful approach to understanding the biology of infectious bacterial pathogens. In this study, a quantitative approach, referred to as Gnom(Cmp), was developed to study the microevolution of bacterial pathogens. Although much more time-consuming than existing tools, this procedure provides a much higher resolution. Gnom(Cmp) accomplishes this by establishing genome-wide heterogeneity genotypes, which are then quantified and comparatively analyzed. The heterogeneity genotypes are defined as chromosomal base positions that have multiple variants within particular genomes, resulted from DNA duplications and subsequent mutations. To prove the concept, the procedure was applied on the genomes of 15 Staphylococcus aureus strains, focusing extensively on two pairs of hVISA/VISA strains. hVISA refers to heteroresistant vancomycin-intermediate S. aureus strains and VISA is their VISA mutants. hVISA/VISA displays some remarkable properties. hVISA is susceptible to vancomycin, but VISA mutants emerge soon after a short period of vancomycin therapy, therefore making the pathogen a great model organism for fast-evolving bacterial pathogens. The analysis indicated that Gnom(Cmp) could reveal variants within the genomes, which can be analyzed within the global genome context. Gnom(Cmp) discovered evolutionary hotspots and their dynamics among many closely related, even isogenic genomes. The analysis thus allows the exploration of the molecular mechanisms behind hVISA/VISA evolution, providing a working hypotheses for experimental testing and validation.

Development of Prevotella nigrescens-specific PCR primers based on the nucleotide sequence of a Pn23 DNA probe

A previous study reported the cloning of a putative Prevotella nigrescens-specific DNA probe, Pn23, using random shotgun method. The present study evaluated the species-specificity of Pn23 for P. nigrescens using the clinical strains of Prevotella intermedia and P. nigrescens to develop P. nigrescens-specific polymerase chain reaction (PCR) primers. Southern blot analysis showed that the DNA probe, Pn23, detected only the genomic DNA of P. nigrescens strains. PCR showed that the two sets of PCR primers, Pn23-F1/Pn23-R1 and Pn23-F2/Pn23-R2, had species-specificity for P. nigrescens. Interestingly, the two sets of PCR primers, Pn23-F6/Pn23-R6 and Pn23-F7/Pn23-R7, had strain-specificity for P. nigrescens ATCC 33563. The detection limits of the four primer sets were 40 or 4 pg of the purified genomic DNA of P. nigrescens ATCC 33563. These results suggest that the DNA probe, Pn23, and the two sets of PCR primers, Pn23-F1/Pn23-R1 and Pn23-F2/Pn23-R2, can be useful for the detection of P. nigrescens in the molecular epidemiological studies of oral infectious diseases.

Genome sequence of an Erwinia amylovora strain with pathogenicity restricted to Rubus plants

Here, we present the genome of a strain of Erwinia amylovora, the fire blight pathogen, with pathogenicity restricted to Rubus spp. Comparative genomics of ATCC BAA-2158 with E. amylovora strains from non-Rubus hosts identified significant genetic differences but support the inclusion of this strain within the species E. amylovora.

Complete genome sequence of Japanese erwinia strain ejp617, a bacterial shoot blight pathogen of pear

The Japanese Erwinia strain Ejp617 is a plant pathogen that causes bacterial shoot blight of pear in Japan. Here, we report the complete genome sequence of strain Ejp617 isolated from Nashi pears in Japan to provide further valuable insight among related Erwinia species.

Evolutionary insights from Erwinia amylovora genomics

Evolutionary genomics is coming into focus with the recent availability of complete sequences for many bacterial species. A hypothesis on the evolution of virulence factors in the plant pathogen Erwinia amylovora, the causative agent of fire blight, was generated using comparative genomics with the genomes E. amylovora, Erwinia pyrifoliae and Erwinia tasmaniensis. Putative virulence factors were mapped to the proposed genealogy of the genus Erwinia that is based on phylogenetic and genomic data. Ancestral origin of several virulence factors was identified, including levan biosynthesis, sorbitol metabolism, three T3SS and two T6SS. Other factors appeared to have been acquired after divergence of pathogenic species, including a second flagellar gene and two glycosyltransferases involved in amylovoran biosynthesis. E. amylovora singletons include 3 unique T3SS effectors that may explain differential virulence/host ranges. E. amylovora also has a unique T1SS export system, and a unique third T6SS gene cluster. Genetic analysis revealed signatures of foreign DNA suggesting that horizontal gene transfer is responsible for some of these differential features between the three species.

Genome sequence of the biocontrol agent Pantoea vagans strain C9-1

Pantoea vagans is a Gram-negative enterobacterial plant epiphyte of a broad range of plants. Here we report the 4.89-Mb genome sequence of P. vagans strain C9-1 (formerly Pantoea agglomerans), which is commercially registered for biological control of fire blight, a disease of pear and apple trees caused by Erwinia amylovora.