bacLIFE: a user-friendly computational workflow for genome analysis and prediction of lifestyle-associated genes in bacteria.

Bacteria have an extensive adaptive ability to live in close association with eukaryotic hosts, exhibiting detrimental, neutral or beneficial effects on host growth and health. However, the genes involved in niche adaptation are mostly unknown and their functions poorly characterized. Here, we present bacLIFE ( https://github.com/Carrion-lab/bacLIFE ) a streamlined computational workflow for genome annotation, large-scale comparative genomics, and prediction of lifestyle-associated genes (LAGs). As a proof of concept, we analyzed 16,846 genomes from the Burkholderia/Paraburkholderia and Pseudomonas genera, which led to the identification of hundreds of genes potentially associated with a plant pathogenic lifestyle. Site-directed mutagenesis of 14 of these predicted LAGs of unknown function, followed by plant bioassays, showed that 6 predicted LAGs are indeed involved in the phytopathogenic lifestyle of Burkholderia plantarii and Pseudomonas syringae pv. phaseolicola. These 6 LAGs encompassed a glycosyltransferase, extracellular binding proteins, homoserine dehydrogenases and hypothetical proteins. Collectively, our results highlight bacLIFE as an effective computational tool for prediction of LAGs and the generation of hypotheses for a better understanding of bacteria-host interactions.

Allelic variation in the indoleacetic acid-lysine synthase gene of the bacterial pathogen Pseudomonas savastanoi and its role in auxin production.

Indole-3-acetic acid (IAA) production is a pathogenicity/virulence factor in the Pseudomonas syringae complex, including Pseudomonas savastanoi. P. savastanoi pathovars (pvs.) genomes contain the iaaL gene, encoding an enzyme that catalyzes the biosynthesis of the less biologically active compound 3-indole-acetyl-ϵ-L-lysine (IAA-Lys). Previous studies have reported the identification of IAA-Lys in culture filtrates of P. savastanoi strains isolated from oleander (pv. nerii), but the conversion of IAA into a conjugate was not detectable in olive strains (pv. savastanoi). In this paper, we show the distribution of iaaL alleles in all available P. savastanoi genomes of strains isolated from woody hosts. Most strains encode two different paralogs, except for those isolated from broom (pv. retacarpa), which contain a single allele. In addition to the three previously reported iaaL alleles (iaaL Psv, iaaL Psn and iaaL Pto), we identified iaaL Psf, an exclusive allele of strains isolated from ash (pv. fraxini). We also found that the production of IAA-Lys in P. savastanoi pv. savastanoi and pv. nerii depends on a functional iaaL Psn allele, whereas in pv. fraxini depends on iaaL Psf. The production of IAA-Lys was detected in cultures of an olive strain heterologously expressing IaaLPsn-1, IaaLPsf-1 and IaaLPsf-3, but not when expressing IaaLPsv-1. In addition, Arabidopsis seedlings treated with the strains overproducing the conjugate, and thus reducing the free IAA content, alleviated the root elongation inhibitory effect of IAA. IAA-Lys synthase activity assays with purified allozymes confirmed the functionality and specificity of lysine as a substrate of IaaLPsn-1 and IaaLPsf-3, with IaaLPsf-3 showing the highest catalytic efficiency for both substrates. The IAA-Lys synthase activity of IaaLPsn-1 was abolished by the insertion of two additional tyrosine residues encoded in the inactive allozyme IaaLPsv-1. These results highlight the relevance of allelic variation in a phytohormone-related gene for the modulation of auxin production in a bacterial phytopathogen.

The Rhizobacterium Pseudomonas alcaligenes AVO110 Induces the Expression of Biofilm-Related Genes in Response to Rosellinia necatrix Exudates.

The rhizobacterium Pseudomonas alcaligenes AVO110 exhibits antagonism toward the phytopathogenic fungus Rosellinia necatrix. This strain efficiently colonizes R. necatrix hyphae and is able to feed on their exudates. Here, we report the complete genome sequence of P. alcaligenes AVO110. The phylogeny of all available P. alcaligenes genomes separates environmental isolates, including AVO110, from those obtained from infected human blood and oyster tissues, which cluster together with Pseudomonas otitidis. Core and pan-genome analyses showed that P. alcaligenes strains encode highly heterogenic gene pools, with the AVO110 genome encoding the largest and most exclusive variable region (~1.6 Mb, 1795 genes). The AVO110 singletons include a wide repertoire of genes related to biofilm formation, several of which are transcriptionally modulated by R. necatrix exudates. One of these genes (cmpA) encodes a GGDEF/EAL domain protein specific to Pseudomonas spp. strains isolated primarily from the rhizosphere of diverse plants, but also from soil and water samples. We also show that CmpA has a role in biofilm formation and that the integrity of its EAL domain is involved in this function. This study contributes to a better understanding of the niche-specific adaptations and lifestyles of P. alcaligenes, including the mycophagous behavior of strain AVO110.

Pseudomonas savastanoi pv. mandevillae pv. nov., a Clonal Pathogen Causing an Emerging, Devastating Disease of the Ornamental Plant Mandevilla spp.

Commercial production of the ornamental plant dipladenia (Mandevilla spp.) is threatened by dipladenia leaf and stem spot disease, caused by the bacterium Pseudomonas savastanoi. P. savastanoi includes four pathovars of woody hosts differentiated by a characteristic host range in olive, oleander, ash, and broom plants. However, isolates from dipladenia have not been ascribed to any particular lineage or P. savastanoi pathovar. Here we report that isolates from dipladenia represent a distinct, clonal lineage. First, dipladenia isolates display very similar plasmid profiles, including a plasmid encoding the iaaM gene for biosynthesis of indole-3-acetic acid. Second, multilocus sequence analysis and core genome single-nucleotide polymorphisms phylogenies showed a monophyletic origin for dipladenia isolates, which cluster with isolates from oleander (pathovar nerii) in a distinct clade well separated from other P. savastanoi strains. Metabolic profiling and cross-pathogenicity tests in olive, oleander, ash, broom, and dipladenia clearly distinguished dipladenia isolates from the four P. savastanoi pathovars. Comparative genomics of the draft genome sequence of the dipladenia strain Ph3 with the other four pathovars showed that Ph3 encodes very few strain-specific genes and a similar set of virulence genes to pv. nerii, including its repertoire of type III secretion system effectors. However, hierarchical clustering based on the catalog of effectors and their allelic variants clearly separated Ph3 from pv. nerii strains. Based on their distinctive pathogenicity profile, we propose a de novo pathovar for P. savastanoi isolates from dipladenia, P. savastanoi pv. mandevillae pv. nov., for which strain Ph3 (CFBP 8832PT) has been designated as the pathotype strain.

Genes ptz and idi, Coding for Cytokinin Biosynthesis Enzymes, Are Essential for Tumorigenesis and In Planta Growth by P. syringae pv. savastanoi NCPPB 3335.

The phytopathogenic bacterium Pseudomonas syringae pv. savastanoi elicits aerial tumors on olive plants and is also able to synthesize large amounts of auxins and cytokinins. The auxin indoleacetic acid was shown to be required for tumorigenesis, but there is only correlational evidence suggesting a role for cytokinins. The model strain NCPPB 3335 contains two plasmid-borne genes coding for cytokinin biosynthesis enzymes: ptz, for an isopentenyl transferase and idi, for an isopentenyl-diphosphate delta-isomerase. Phylogenetic analyses showed that carriage of ptz and idi is not strictly associated with tumorigenic bacteria, that both genes were linked when first acquired by P. syringae, and that a different allele of ptz has been independently acquired by P. syringae pv. savastanoi and closely related bacteria. We generated mutant derivatives of NCPPB 3335 cured of virulence plasmids or with site-specific deletions of genes ptz and/or idi and evaluated their virulence in lignified and micropropagated olive plants. Strains lacking ptz, idi, or both produced tumors with average volumes up to 29 times smaller and reached populations up to two orders of magnitude lower than those induced by strain NCPPB 3335; these phenotypes reverted by complementation with the cloned genes. Trans-zeatin was the most abundant cytokinin in culture filtrates of NCPPB 3335. Deletion of gene ptz abolished biosynthesis of trans-zeatin and dihydrozeatin, whereas a reduced but significant amount of isopentenyladenine was still detected in the medium, suggesting the existence of other genes contributing to cytokinin biosynthesis in P. syringae. Conversely, extracts from strains lacking gene idi contained significantly higher amounts of trans-zeatin than extracts from the wild-type strain but similar amounts of the other cytokinins. This suggests that Idi might promote tumorigenesis by ensuring the biosynthesis of the most active cytokinin forms, their correct balance in planta, or by regulating the expression of other virulence genes. Therefore, gene ptz, but not gene idi, is essential for the biosynthesis of high amounts of cytokinins in culture; however, both ptz and idi are individually essential for the adequate development of tumors on olive plants by Psv NCPPB 3335.

Host Range Determinants of Pseudomonas savastanoi Pathovars of Woody Hosts Revealed by Comparative Genomics and Cross-Pathogenicity Tests.

The study of host range determinants within the Pseudomonas syringae complex is gaining renewed attention due to its widespread distribution in non-agricultural environments, evidence of large variability in intra-pathovar host range, and the emergence of new epidemic diseases. This requires the establishment of appropriate model pathosystems facilitating integration of phenotypic, genomic and evolutionary data. Pseudomonas savastanoi pv. savastanoi is a model pathogen of the olive tree, and here we report a closed genome of strain NCPPB 3335, plus draft genome sequences of three strains isolated from oleander (pv. nerii), ash (pv. fraxini) and broom plants (pv. retacarpa). We then conducted a comparative genomic analysis of these four new genomes plus 16 publicly available genomes, representing 20 strains of these four P. savastanoi pathovars of woody hosts. Despite overlapping host ranges, cross-pathogenicity tests using four plant hosts clearly separated these pathovars and lead to pathovar reassignment of two strains. Critically, these functional assays were pivotal to reconcile phylogeny with host range and to define pathovar-specific genes repertoires. We report a pan-genome of 7,953 ortholog gene families and a total of 45 type III secretion system effector genes, including 24 core genes, four genes exclusive of pv. retacarpa and several genes encoding pathovar-specific truncations. Noticeably, the four pathovars corresponded with well-defined genetic lineages, with core genome phylogeny and hierarchical clustering of effector genes closely correlating with pathogenic specialization. Knot-inducing pathovars encode genes absent in the canker-inducing pv. fraxini, such as those related to indole acetic acid, cytokinins, rhizobitoxine, and a bacteriophytochrome. Other pathovar-exclusive genes encode type I, type II, type IV, and type VI secretion system proteins, the phytotoxine phevamine A, a siderophore, c-di-GMP-related proteins, methyl chemotaxis proteins, and a broad collection of transcriptional regulators and transporters of eight different superfamilies. Our combination of pathogenicity analyses and genomics tools allowed us to correctly assign strains to pathovars and to propose a repertoire of host range-related genes in the P. syringae complex.

Response of the Biocontrol Agent Pseudomonas pseudoalcaligenes AVO110 to Rosellinia necatrix Exudate.

The rhizobacterium Pseudomonas pseudoalcaligenes AVO110, isolated by the enrichment of competitive avocado root tip colonizers, controls avocado white root rot disease caused by Rosellinia necatrix Here, we applied signature-tagged mutagenesis (STM) during the growth and survival of AVO110 in fungal exudate-containing medium with the goal of identifying the molecular mechanisms linked to the interaction of this bacterium with R. necatrix A total of 26 STM mutants outcompeted by the parental strain in fungal exudate, but not in rich medium, were selected and named growth-attenuated mutants (GAMs). Twenty-one genes were identified as being required for this bacterial-fungal interaction, including membrane transporters, transcriptional regulators, and genes related to the metabolism of hydrocarbons, amino acids, fatty acids, and aromatic compounds. The bacterial traits identified here that are involved in the colonization of fungal hyphae include proteins involved in membrane maintenance (a dynamin-like protein and ColS) or cyclic-di-GMP signaling and chemotaxis. In addition, genes encoding a DNA helicase (recB) and a regulator of alginate production (algQ) were identified as being required for efficient colonization of the avocado rhizosphere.IMPORTANCE Diseases associated with fungal root invasion cause a significant loss of fruit tree production worldwide. The bacterium Pseudomonas pseudoalcaligenes AVO110 controls avocado white root rot disease caused by Rosellinia necatrix by using mechanisms involving competition for nutrients and niches. Here, a functional genomics approach was conducted to identify the bacterial traits involved in the interaction with this fungal pathogen. Our results contribute to a better understanding of the multitrophic interactions established among bacterial biocontrol agents, the plant rhizosphere, and the mycelia of soilborne pathogens.

Knots Untie: Molecular Determinants Involved in Knot Formation Induced by Pseudomonas savastanoi in Woody Hosts.

The study of the molecular basis of tree diseases is lately receiving a renewed attention, especially with the emerging perception that pathogens require specific pathogenicity and virulence factors to successfully colonize woody hosts. Pathosystems involving woody plants are notoriously difficult to study, although the use of model bacterial strains together with genetically homogeneous micropropagated plant material is providing a significant impetus to our understanding of the molecular determinants leading to disease. The gammaproteobacterium Pseudomonas savastanoi belongs to the intensively studied Pseudomonas syringae complex, and includes three pathogenic lineages causing tumorous overgrowths (knots) in diverse economically relevant trees and shrubs. As it occurs with many other bacteria, pathogenicity of P. savastanoi is dependent on a type III secretion system, which is accompanied by a core set of at least 20 effector genes shared among strains isolated from olive, oleander, and ash. The induction of knots of wild-type size requires that the pathogen maintains adequate levels of diverse metabolites, including the phytohormones indole-3-acetic acid and cytokinins, as well as cyclic-di-GMP, some of which can also regulate the expression of other pathogenicity and virulence genes and participate in bacterial competitiveness. In a remarkable example of social networking, quorum sensing molecules allow for the communication among P. savastanoi and other members of the knot microbiome, while at the same time are essential for tumor formation. Additionally, a distinguishing feature of bacteria from the P. syringae complex isolated from woody organs is the possession of a 15 kb genomic island (WHOP) carrying four operons and three other genes involved in degradation of phenolic compounds. Two of these operons mediate the catabolism of anthranilate and catechol and, together with another operon, are required for the induction of full-size tumors in woody hosts, but not in non-woody micropropagated plants. The use of transposon mutagenesis also uncovered a treasure trove of additional P. savastanoi genes affecting virulence and participating in diverse bacterial processes. Although there is still much to be learned on what makes a bacterium a successful pathogen of trees, we are already untying the knots.