Bacterial motility plays a crucial role in competitiveness and colonization in the rhizosphere. In this work, Chromatin ImmunoPrecipitation Sequencing (ChIP-seq) analysis has been used to identify genes putatively regulated by the transcriptional regulatory protein FleQ in Pseudomonas fluorescens F113 and Pseudomonas putida KT2440. This protein was previously identified as a master regulator of flagella and biofilm formation in both strains. This work has demonstrated that FleQ from both bacteria are conserved and functionally equivalent for motility regulation. Furthermore, the ChIP-seq analysis has shown that FleQ is a global regulator with the identification of 121 and 103 FleQ putative binding sites in P. fluorescens F113 and P. putida KT2440 respectively. Putative genes regulated by FleQ included, as expected, flagellar and motility-related genes and others involved in adhesion and exopolysaccharide production. Surprisingly, the ChIP-seq analysis also identified iron homeostasis-related genes for which positive regulation was shown by RT-qPCR. The results also showed that FleQ from P. fluorescens F113 shares an important part of its direct regulon with AmrZ, a global regulator also implicated in environmental adaption. Although AmrZ also regulates motility and iron uptake, the overlap occurred mostly with the iron-related genes, since both regulators control a different set of motility-related genes.

The transcriptional regulator AmrZ is a global regulatory protein conserved within the pseudomonads. AmrZ can act both as a positive and a negative regulator of gene expression, controlling many genes implicated in environmental adaption. Regulated traits include motility, iron homeostasis, exopolysaccharides production and the ability to form biofilms. In Pseudomonas fluorescens F113, an amrZ mutant presents a pleiotropic phenotype, showing increased swimming motility, decreased biofilm formation and very limited ability for competitive colonization of rhizosphere, its natural habitat. It also shows different colony morphology and binding of the dye Congo Red. The amrZ mutant presents severely reduced levels of the messenger molecule cyclic-di-GMP (c-di-GMP), which is consistent with the motility and biofilm formation phenotypes. Most of the genes encoding proteins with diguanylate cyclase (DGCs) or phosphodiesterase (PDEs) domains, implicated in c-di-GMP turnover in this bacterium, appear to be regulated by AmrZ. Phenotypic analysis of eight mutants in genes shown to be directly regulated by AmrZ and encoding c-di-GMP related enzymes, showed that seven of them were altered in motility and/or biofilm formation. The results presented here show that in P. fluorescens, AmrZ determines c-di-GMP levels through the regulation of a complex network of genes encoding DGCs and PDEs.

The Pseudomonas fluorescens complex of species includes plant-associated bacteria with potential biotechnological applications in agriculture and environmental protection. Many of these bacteria can promote plant growth by different means, including modification of plant hormonal balance and biocontrol. The P. fluorescens group is currently divided into eight major subgroups in which these properties and many other ecophysiological traits are phylogenetically distributed. Therefore, a rapid phylogroup assignment for a particular isolate could be useful to simplify the screening of putative inoculants. By using comparative genomics on 71 P. fluorescens genomes, we have identified nine markers which allow classification of any isolate into these eight subgroups, by a presence/absence PCR test. Nine primer pairs were developed for the amplification of these markers. The specificity and sensitivity of these primer pairs were assessed on 28 field isolates, environmental samples from soil and rhizosphere and tested by in silico PCR on 421 genomes. Phylogenomic analysis validated the results: the PCR-based system for classification of P. fluorescens isolates has a 98.34% of accuracy and it could be used as a rapid and simple assay to evaluate the potential of any P. fluorescens complex strain.

Biofilms are complex structures that are crucial during host–bacteria interaction and colonization. Bacteria within biofilms are surrounded by an extracellular matrix (ECM) typically composed of proteins, polysaccharides, lipids, and DNA. Pseudomonads contain a variety of ECM components, some of which have been extensively characterized. However, neither the ECM composition of plant-associated pseudomonads nor their phylogenetic distribution within the genus has been so thoroughly studied. In this work, we use in silico methods to describe the ECM composition of Pseudomonas fluorescens F113, a plant growth-promoting rhizobacteria and model for rhizosphere colonization. These components include the polysaccharides alginate, poly-N-acetyl-glucosamine (PNAG) and levan; the adhesins LapA, MapA and PsmE; and the functional amyloids in Pseudomonas. Interestingly, we identified novel components: the Pseudomonas acidic polysaccharide (Pap), whose presence is limited within the genus; and a novel type of Flp/Tad pilus, partially different from the one described in P. aeruginosa. Furthermore, we explored the phylogenetic distribution of the most relevant ECM components in nearly 600 complete Pseudomonas genomes. Our analyses show that Pseudomonas populations contain a diverse set of gene/gene clusters potentially involved in the formation of their ECMs, showing certain commensal versus pathogen lifestyle specialization.

The genome of Pseudomonas fluorescens F113, a model rhizobacterium and a plant growth-promoting agent, encodes three putative type VI secretion systems (T6SSs); F1-, F2- and F3-T6SS. Bioinformatic analysis of the F113 T6SSs has revealed that they belong to group 3, group 1.1, and group 4a, respectively, similar to those previously described in Pseudomonas aeruginosa. In addition, in silico analyses allowed us to identify genes encoding a total of five orphan VgrG proteins and eight putative effectors (Tfe), some with their cognate immunity protein (Tfi) pairs. Genes encoding Tfe and Tfi are found in the proximity of P. fluorescens F113 vgrG, hcp, eagR and tap genes. RNA-Seq analyses in liquid culture and rhizosphere have revealed that F1- and F3-T6SS are expressed under all conditions, indicating that they are active systems, while F2-T6SS did not show any relevant expression under the tested conditions. The analysis of structural mutants in the three T6SSs has shown that the active F1- and F3-T6SSs are involved in interbacterial killing while F2 is not active in these conditions and its role is still unknown.. A rhizosphere colonization analysis of the double mutant affected in the F1- and F3-T6SS clusters showed that the double mutant was severely impaired in persistence in the rhizosphere microbiome, revealing the importance of these two systems for rhizosphere adaption.

The complete genome sequence of Rhodococcus sp. WAY2 (WAY2) consists of a circular chromosome, three linear replicons and a small circular plasmid. The linear replicons contain typical actinobacterial invertron-type telomeres with the central CGTXCGC motif. Comparative phylogenetic analysis of the 16S rRNA gene along with phylogenomic analysis based on the genome-to-genome blast distance phylogeny (GBDP) algorithm and digital DNA–DNA hybridization (dDDH) with other Rhodococcus type strains resulted in a clear differentiation of WAY2, which is likely a new species. The genome of WAY2 contains five distinct clusters of bph, etb and nah genes, putatively involved in the degradation of several aromatic compounds. These clusters are distributed throughout the linear plasmids. The high sequence homology of the ring-hydroxylating subunits of these systems with other known enzymes has allowed us to model the range of aromatic substrates they could degrade. Further functional characterization revealed that WAY2 was able to grow with biphenyl, naphthalene and xylene as sole carbon and energy sources, and could oxidize multiple aromatic compounds, including ethylbenzene, phenanthrene, dibenzofuran and toluene. In addition, WAY2 was able to co-metabolize 23 polychlorinated biphenyl congeners, consistent with the five different ring-hydroxylating systems encoded by its genome. WAY2 could also use n-alkanes of various chain-lengths as a sole carbon source, probably due to the presence of alkB and ladA gene copies, which are only found in its chromosome. These results show that WAY2 has a potential to be used for the biodegradation of multiple organic compounds.

Dual flagellar systems have been described in several bacterial genera, but the extent of their prevalence has not been fully explored. Bradyrhizobium diazoefficiens USDA 110T possesses two flagellar systems, the subpolar and the lateral flagella. The lateral flagellum of Bradyrhizobium displays no obvious role, since its performance is explained by cooperation with the subpolar flagellum. In contrast, the lateral flagellum is the only type of flagella present in the related Rhizobiaceae family. In this work, we have analyzed the phylogeny of the Bradyrhizobium genus by means of Genome-to-Genome Blast Distance Phylogeny (GBDP) and Average Nucleotide Identity (ANI) comparisons of 128 genomes and divided it into 13 phylogenomic groups. While all the Bradyrhizobium genomes encode the subpolar flagellum, none of them encodes only the lateral flagellum. The simultaneous presence of both flagella is exclusive of the B. japonicum phylogenomic group. Additionally, 292 Rhizobiales order genomes were analyzed and both flagellar systems are present together in only nine genera. Phylogenetic analysis of 150 representative Rhizobiales genomes revealed an uneven distribution of these flagellar systems. While genomes within and close to the Rhizobiaceae family only possess the lateral flagellum, the subpolar flagellum is exclusive of more early-diverging families, where certain genera also present both flagella.

Rhizosphere colonization by bacteria involves molecular and cellular mechanisms, such as motility and chemotaxis, biofilm formation, metabolic versatility, or biosynthesis of secondary metabolites, among others. Nonetheless, there is limited knowledge concerning the main regulatory factors that drive the rhizosphere colonization process. Here we show the importance of the AmrZ and FleQ transcription factors for adaption in the plant growth-promoting rhizobacterium (PGPR) and rhizosphere colonization model Pseudomonas ogarae F113. RNA-Seq analyses of P. ogarae F113 grown in liquid cultures either in exponential and stationary growth phase, and rhizosphere conditions, revealed that rhizosphere is a key driver of global changes in gene expression in this bacterium. Regarding the genetic background, this work has revealed that a mutation in fleQ causes considerably more alterations in the gene expression profile of this bacterium than a mutation in amrZ under rhizosphere conditions. The functional analysis has revealed that in P. ogarae F113, the transcription factors AmrZ and FleQ regulate genes involved in diverse bacterial functions. Notably, in the rhizosphere, these transcription factors antagonistically regulate genes related to motility, biofilm formation, nitrogen, sulfur, and amino acid metabolism, transport, signalling, and secretion, especially the type VI secretion systems. These results define the regulon of two important bifunctional transcriptional regulators in pseudomonads during the process of rhizosphere colonization.