The Pseudomonas fluorescens complex includes Pseudomonas strains that have been taxonomically assigned to more than fifty different species, many of which have been described as plant growth-promoting rhizobacteria (PGPR) with potential applications in biocontrol and biofertilization. So far the phylogeny of this complex has been analyzed according to phenotypic traits, 16S rDNA, MLSA and inferred by whole-genome analysis. However, since most of the type strains have not been fully sequenced and new species are frequently described, correlation between taxonomy and phylogenomic analysis is missing. In recent years, the genomes of a large number of strains have been sequenced, showing important genomic heterogeneity and providing information suitable for genomic studies that are important to understand the genomic and genetic diversity shown by strains of this complex. Based on MLSA and several whole-genome sequence-based analyses of 93 sequenced strains, we have divided the P. fluorescens complex into eight phylogenomic groups that agree with previous works based on type strains. Digital DDH (dDDH) identified 69 species and 75 subspecies within the 93 genomes. The eight groups corresponded to clustering with a threshold of 31.8% dDDH, in full agreement with our MLSA. The Average Nucleotide Identity (ANI) approach showed inconsistencies regarding the assignment to species and to the eight groups. The small core genome of 1,334 CDSs and the large pan-genome of 30,848 CDSs, show the large diversity and genetic heterogeneity of the P. fluorescens complex. However, a low number of strains were enough to explain most of the CDSs diversity at core and strain-specific genomic fractions. Finally, the identification and analysis of group-specific genome and the screening for distinctive characters revealed a phylogenomic distribution of traits among the groups that provided insights into biocontrol and bioremediation applications as well as their role as PGPR.

Microbial communities play a crucial role in supporting plant health and productivity. Reproducible, natural plant-associated microbiomes can help disentangle microbial dynamics across time and space. Here, using a sequential propagation strategy, we generated a complex and reproducible wheat rhizosphere microbiome (RhizCom) to study successional dynamics and interactions between the soil and heritable seed-borne rhizosphere microbiomes (SbRB) in a microcosm. Using 16S rRNA sequencing and genome-resolved shotgun metagenomics, we find that SbRB surpassed native soil microbes as the dominant rhizosphere-associated microbiome source. SbRB genomes were enriched in host-associated traits including degradation of key saccharide (niche partitioning) and cross-feeding interactions that supported partner strains (niche facilitation). In vitro co-culture experiments confirmed that helper SbRB strains facilitated the growth of partner bacteria on disaccharides as sole carbon source. These results reveal the importance of seed microbiota dynamics in microbial succession and community assembly, which could inform strategies for crop microbiome manipulation.

Diesel is a complex pollutant composed of a mixture of aliphatic and aromatic hydrocarbons. Because of this complexity, diesel bioremediation requires multiple microorganisms, which harbor the catabolic pathways to degrade the mixture. By enrichment cultivation of rhizospheric soil from a diesel-polluted site, we have isolated a bacterial consortium that can grow aerobically with diesel and different alkanes and polycyclic aromatic hydrocarbons (PAHs) as the sole carbon and energy source. Microbiome diversity analyses based on 16S rRNA gene showed that the diesel-degrading consortium consists of 76 amplicon sequence variants (ASVs) and it is dominated by Pseudomonas, Aquabacterium, Chryseobacterium, and Sphingomonadaceae. Changes in microbiome composition were observed when growing on specific hydrocarbons, reflecting that different populations degrade different hydrocarbons. Shotgun metagenome sequence analysis of the consortium growing on diesel has identified redundant genes encoding enzymes implicated in the initial oxidation of alkanes (AlkB, LadA, CYP450) and a variety of hydroxylating and ring-cleavage dioxygenases involved in aromatic and polyaromatic hydrocarbon degradation. The phylogenetic assignment of these enzymes to specific genera allowed us to model the role of specific populations in the diesel-degrading consortium. Rhizoremediation of diesel-polluted soil microcosms using the consortium, resulted in an important enhancement in the reduction of total petroleum hydrocarbons (TPHs), making it suited for rhizoremediation applications.

Background: Plant-beneficial bacterial inoculants are of great interest in agriculture as they have the potential to promote plant growth and health. However, the inoculation of the rhizosphere microbiome often results in a suboptimal or transient colonization, which is due to a variety of factors that influence the fate of the inoculant. To better understand the fate of plant-beneficial inoculants in complex rhizosphere microbiomes, composed by hundreds of genotypes and multifactorial selection mechanisms, controlled studies with high-complexity soil microbiomes are needed. Results: We analysed early compositional changes in a taxa-rich natural soil bacterial community under both exponential nutrient-rich and stationary nutrient-limited growth conditions (i.e. growing and stable communities, respectively) following inoculation with the plant-beneficial bacterium Pseudomonas protegens in a bulk soil or a wheat rhizosphere environment. P. protegens successfully established under all conditions tested and was more abundant in the rhizosphere of the stable community. Nutrient availability was a major factor driving microbiome composition and structure as well as the underlying assembly processes. While access to nutrients resulted in communities assembled mainly by homogeneous selection, stochastic processes dominated under the nutrient-deprived conditions. We also observed an increased rhizosphere selection effect under nutrient-limited conditions, resulting in a higher number of amplicon sequence variants (ASVs) whose relative abundance was enriched. The inoculation with P. protegens produced discrete changes, some of which involved other Pseudomonas. Direct competition between Pseudomonas strains partially failed to replicate the observed differences in the microbiome and pointed to a more complex interaction network. Conclusions: The results of this study show that nutrient availability is a major driving force of microbiome composition, structure and diversity in both the bulk soil and the wheat rhizosphere and determines the assembly processes that govern early microbiome development. The successful establishment of the inoculant was facilitated by the wheat rhizosphere and produced discrete changes among other members of the microbiome. Direct competition between Pseudomonas strains only partially explained the microbiome changes, indicating that indirect interactions or spatial distribution in the rhizosphere or soil interface may be crucial for the survival of certain bacteria.

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.

Polychlorinated biphenyls (PCBs) are widespread persistent pollutants that cause several adverse health effects. Aerobic bioremediation of PCBs involves the activity of either one bacterial species or a microbial consortium. Using multiple species will enhance the range of PCB congeners co-metabolized since different PCB-degrading microorganisms exhibit different substrate specificity. We have isolated a bacterial consortium by successive enrichment culture using biphenyl (analog of PCBs) as the sole carbon and energy source. This consortium is able to grow on biphenyl, benzoate, and protocatechuate. Whole-community DNA extracted from the consortium was used to analyze biodiversity by Illumina sequencing of a 16S rRNA gene amplicon library and to determine the metagenome by whole-genome shotgun Illumina sequencing. Biodiversity analysis shows that the consortium consists of 24 operational taxonomic units (≥97% identity). The consortium is dominated by strains belonging to the genus Pseudomonas, but also contains betaproteobacteria and Rhodococcus strains. whole-genome shotgun (WGS) analysis resulted in contigs containing 78.3 Mbp of sequenced DNA, representing around 65% of the expected DNA in the consortium. Bioinformatic analysis of this metagenome has identified the genes encoding the enzymes implicated in three pathways for the conversion of biphenyl to benzoate and five pathways from benzoate to tricarboxylic acid (TCA) cycle intermediates, allowing us to model the whole biodegradation network. By genus assignment of coding sequences, we have also been able to determine that the three biphenyl to benzoate pathways are carried out by Rhodococcus strains. In turn, strains belonging to Pseudomonas and Bordetella are the main responsible of three of the benzoate to TCA pathways while the benzoate conversion into TCA cycle intermediates via benzoyl-CoA and the catechol meta-cleavage pathways are carried out by beta proteobacteria belonging to genera such as Achromobacter and Variovorax. We have isolated a Rhodococcus strain WAY2 from the consortium which contains the genes encoding the three biphenyl to benzoate pathways indicating that this strain is responsible for all the biphenyl to benzoate transformations. The presented results show that metagenomic analysis of consortia allows the identification of bacteria active in biodegradation processes and the assignment of specific reactions and pathways to specific bacterial groups.