Genomic and metabolic differences between Pseudomonas putida populations inhabiting sugarcane rhizosphere or bulk soil

Pseudomonas putida is one of 13 major groups of Pseudomonas spp. and contains numerous species occupying diverse niches and performing many functions such as plant growth promotion and bioremediation. Here we compared a set of 19 P. putida isolates obtained from sugarcane rhizosphere or bulk soil using a population genomics approach aiming to assess genomic and metabolic differences between populations from these habitats. Phylogenomics placed rhizosphere versus bulk soil strains in separate clades clustering with different type strains of the P. putida group. Multivariate analyses indicated that the rhizosphere and bulk soil isolates form distinct populations. Comparative genomics identified several genetic functions (GO-terms) significantly different between populations, including some exclusively present in the rhizosphere or bulk soil strains, such as D-galactonic acid catabolism and cellulose biosynthesis, respectively. The metabolic profiles of rhizosphere and bulk soil populations analyzed by Biolog Ecoplates also differ significantly, most notably by the higher oxidation of D-galactonic/D-galacturonic acid by the rhizosphere population. Accordingly, D-galactonate catabolism operon (dgo) was present in all rhizosphere isolates and absent in the bulk soil population. This study showed that sugarcane rhizosphere and bulk soil harbor different populations of P. putida and identified genes and functions potentially associated with their soil niches.