Draft Genome Sequence of Pseudomonas veronii Strain 1YdBTEX2

Fitness-Conditional Genes for Soil Adaptation in the Bioaugmentation Agent Pseudomonas veronii 1YdBTEX2

Strain inoculation (bioaugmentation) is a potentially useful technology to provide microbiomes with new functionalities. However, there is limited understanding of the genetic factors contributing to successful establishment of inoculants. This work aimed to characterize the genes implicated in proliferation of the monoaromatic compound-degrading Pseudomonas veronii 1YdBTEX2 in nonsterile polluted soils. We generated two independent mutant libraries by random minitransposon-delivered marker insertion followed by deep sequencing (Tn-seq) with a total of 5.0 × 10⁵ unique insertions. Libraries were grown in multiple successive cycles for up to 50 generations either in batch liquid medium or in two types of soil microcosms with different resident microbial content (sand or silt) in the presence of toluene. Analysis of gene insertion abundances at different time points (passed generations of metapopulation growth), in comparison to proportions at start and to in silico generated randomized insertion distributions, allowed to define ~800 essential genes common to both libraries and ~2,700 genes with conditional fitness effects in either liquid or soil (195 of which resulted in fitness gain). Conditional fitness genes largely overlapped among all growth conditions but affected approximately twice as many functions in liquid than in soil. This indicates soil to be a more promiscuous environment for mutant growth, probably because of additional nutrient availability. Commonly depleted genes covered a wide range of biological functions and metabolic pathways, such as inorganic ion transport, fatty acid metabolism, amino acid biosynthesis, or nucleotide and cofactor metabolism. Only sparse gene sets were uncovered whose insertion caused fitness decrease exclusive for soils, which were different between silt and sand. Despite detectable higher resident bacteria and potential protist predatory counts in silt, we were, therefore, unable to detect any immediately obvious candidate genes affecting P. veronii biological competitiveness. In contrast to liquid growth conditions, mutants inactivating flagella biosynthesis and motility consistently gained strong fitness advantage in soils and displayed higher growth rates than wild type. In conclusion, although many gene functions were found to be important for growth in soils, most of these are not specific as they affect growth in liquid minimal medium more in general. This indicates that P. veronii does not need major metabolic reprogramming for proliferation in soil with accessible carbon and generally favorable growth conditions. IMPORTANCE Restoring damaged microbiomes is still a formidable challenge. Classical widely adopted approaches consist of augmenting communities with pure or mixed cultures in the hope that these display their intended selected properties under in situ conditions. Ecological theory, however, dictates that introduction of a nonresident microbe is unlikely to lead to its successful proliferation in a foreign system such as a soil microbiome. In an effort to study this systematically, we used random transposon insertion scanning to identify genes and possibly, metabolic subsystems, that are crucial for growth and survival of a bacterial inoculant (Pseudomonas veronii) for targeted degradation of monoaromatic compounds in contaminated nonsterile soils. Our results indicate that although many gene functions are important for proliferation in soil, they are general factors for growth and not exclusive for soil. In other words, P. veronii is a generalist that is not a priori hindered by the soil for its proliferation and would make a good bioaugmentation candidate.

Genomic analysis of dibenzofuran-degrading Pseudomonas veronii strain Pvy reveals its biodegradative versatility

Certain industrial chemicals accumulate in the environment due to their recalcitrant properties. Bioremediation uses the capability of some environmental bacteria to break down these chemicals and attenuate the pollution. One such bacterial strain, designated Pvy, was isolated from sediment samples from a lagoon in Romania located near an oil refinery due to its capacity to degrade dibenzofuran (DF). The genome sequence of the Pvy strain was obtained using an Oxford Nanopore MiniION platform. According to the consensus 16S rRNA gene sequence that was compiled from six 16S rRNA gene copies contained in the genome and orthologous average nucleotide identity (OrthoANI) calculation, the Pvy strain was identified as Pseudomonas veronii, which confirmed the identification obtained with the aid of MALDI-TOF mass spectrometry and MALDI BioTyper. The genome was analyzed with respect to enzymes responsible for the overall biodegradative versatility of the strain. The Pvy strain was able to derive carbon from naphthalene (NP) and several aromatic compounds of natural origin, including salicylic, protocatechuic, p-hydroxybenzoic, trans-cinnamic, vanillic, and indoleacetic acids or vanillin, and was shown to degrade but not utilize DF. In total seven loci were found in the Pvy genome, which enables the strain to participate in the degradation of these aromatic compounds. Our experimental data also indicate that the transcription of the NP-dioxygenase α-subunit gene (ndoB), carried by the plasmid of the Pvy strain, is inducible by DF. These features make the Pvy strain a potential candidate for various bioremediation applications.

Mechanistic insights into bacterial metabolic reprogramming from omics-integrated genome-scale models

Understanding the adaptive responses of individual bacterial strains is crucial for microbiome engineering approaches that introduce new functionalities into complex microbiomes, such as xenobiotic compound metabolism for soil bioremediation. Adaptation requires metabolic reprogramming of the cell, which can be captured by multi-omics, but this data remains formidably challenging to interpret and predict. Here we present a new approach that combines genome-scale metabolic modeling with transcriptomics and exometabolomics, both of which are common tools for studying dynamic population behavior. As a realistic demonstration, we developed a genome-scale model of Pseudomonas veronii 1YdBTEX2, a candidate bioaugmentation agent for accelerated metabolism of mono-aromatic compounds in soil microbiomes, while simultaneously collecting experimental data of P. veronii metabolism during growth phase transitions. Predictions of the P. veronii growth rates and specific metabolic processes from the integrated model closely matched experimental observations. We conclude that integrative and network-based analysis can help build predictive models that accurately capture bacterial adaptation responses. Further development and testing of such models may considerably improve the successful establishment of bacterial inoculants in more complex systems.

Spatial organization of bacterial populations in response to oxygen and carbon counter-gradients in pore networks

Microbial activity in soil is spatially heterogeneous often forming spatial hotspots that contribute disproportionally to biogeochemical processes. Evidence suggests that bacterial spatial organization contributes to the persistence of anoxic hotspots even in unsaturated soils. Such processes are difficult to observe in situ at the microscale, hence mechanisms and time scales relevant for bacterial spatial organization remain largely qualitative. Here we develop an experimental platform based on glass-etched micrometric pore networks that mimics resource gradients postulated in soil aggregates to observe spatial organization of fluorescently tagged aerobic and facultative anaerobic bacteria. Two initially intermixed bacterial species, Pseudomonas putida and Pseudomonas veronii, segregate into preferential regions promoted by opposing gradients of carbon and oxygen (such persistent coexistence is not possible in well-mixed cultures). The study provides quantitative visualization and modeling of bacterial spatial organization within aggregate-like hotspots, a key step towards developing a mechanistic representation of bacterial community organization in soil pores.

Candidate nematicidal proteins in a new Pseudomonas veronii isolate identified by its antagonistic properties against Xiphinema index

The nematode Xiphinema index affects grape vines and transmits important viruses associated with fanleaf degeneration. Pseudomonas spp. are an extensive bacterial group in which important biodegradation and/or biocontrol properties can occur for several strains in the group. The aim of this study was to identify new Pseudomonas isolates with antagonist activity against X. index. Forty bacterial isolates were obtained from soil and root samples from Chilean vineyards. Thirteen new fluorescent pseudomonads were found and assessed for their antagonistic capability. The nematicide Pseudomonas protegens CHA0 was used as a control. Challenges of nematode individuals in King's B semi-solid agar Petri dishes facilitated the identification of the Pseudomonas veronii isolate R4, as determined by a 16S rRNA sequence comparison. This isolate was as effective as CHA0 as an antagonist of X. index, although it had a different lethality kinetic. Milk-induced R4 cultures exhibited protease and lipase activities in cell supernatants using both gelatin/tributyrin Petri dish assays and zymograms. Three proteins with these activities were isolated and subjected to mass spectrometry. Amino acid partial sequences enabled the identification of a 49-kDa protease similar to metalloprotease AprA and two lipases of 50 kDa and 69 kDa similar to LipA and ExoU, respectively. Electron microscopy analyses of challenged nematodes revealed degraded cuticle after R4 supernatant treatment. These results represent a new and unexplored property in this species associated with the presence of secretable lipases and protease, similar to characterized enzymes present in biocontrol pseudomonads.

The Genome of the Toluene-Degrading Pseudomonas veronii Strain 1YdBTEX2 and Its Differential Gene Expression in Contaminated Sand

The natural restoration of soils polluted by aromatic hydrocarbons such as benzene, toluene, ethylbenzene and m- and p-xylene (BTEX) may be accelerated by inoculation of specific biodegraders (bioaugmentation). Bioaugmentation mainly involves introducing bacteria that deploy their metabolic properties and adaptation potential to survive and propagate in the contaminated environment by degrading the pollutant. In order to better understand the adaptive response of cells during a transition to contaminated material, we analyzed here the genome and short-term (1 h) changes in genome-wide gene expression of the BTEX-degrading bacterium Pseudomonas veronii 1YdBTEX2 in non-sterile soil and liquid medium, both in presence or absence of toluene. We obtained a gapless genome sequence of P. veronii 1YdBTEX2 covering three individual replicons with a total size of 8 Mb, two of which are largely unrelated to current known bacterial replicons. One-hour exposure to toluene, both in soil and liquid, triggered massive transcription (up to 208-fold induction) of multiple gene clusters, such as toluene degradation pathway(s), chemotaxis and toluene efflux pumps. This clearly underlines their key role in the adaptive response to toluene. In comparison to liquid medium, cells in soil drastically changed expression of genes involved in membrane functioning (e.g., lipid composition, lipid metabolism, cell fatty acid synthesis), osmotic stress response (e.g., polyamine or trehalose synthesis, uptake of potassium) and putrescine metabolism, highlighting the immediate response mechanisms of P. veronii 1YdBTEX2 for successful establishment in polluted soil.

A draft genome sequence of Pseudomonas veronii R4: a grapevine (Vitis vinifera L.) root-associated strain with high biocontrol potential

A new plant commensal Pseudomonas veronii isolate (strain R4) was identified from a Xiphinema index biocontrol screen. Isolated from grapevine roots from vineyards in central Chile, the strain R4 exhibited a slower yet equivalently effective nematicide activity as the well-characterized P. protegens CHA0. Whole genome sequencing of strain R4 and comparative analysis among the available Pseudomonas spp. genomes allowed for the identification of gene clusters that encode putative extracellular proteases and lipase synthesis and secretion systems, which are proposed to mediate—at least in part—the observed nematicidal activity. In addition, R4 strain presented relevant gene clusters related to metal tolerance, which is typical in P. veronii. Bioinformatics analyses also showed gene clusters associated with plant growth promoting activity, such as indole-3-acetic acid synthesis. In addition, the strain R4 genome presented a metabolic gene clusters associated with phosphate and ammonia biotransformation from soil, which could improve their availability for plants.

N2 Gas Flushing Alleviates the Loss of Bacterial Diversity and Inhibits Psychrotrophic Pseudomonas during the Cold Storage of Bovine Raw Milk

The quality and safety of raw milk still remains a worldwide challenge. Culture-dependent methods indicated that the continuous N2 gas-flushing of raw milk reduced the bacterial growth during cold storage by up to four orders of magnitude, compared to cold storage alone. This study investigated the influence of N2 gas-flushing on bacterial diversity in bovine raw-milk samples, that were either cold stored at 6°C or additionally flushed with pure N2 for up to one week. Next-generation sequencing (NGS) of the V1-V2 hypervariable regions of 16S rRNA genes, derived from amplified cDNA, which was obtained from RNA directly isolated from raw-milk samples, was performed. The reads, which were clustered into 2448 operational taxonomic units (OTUs), were phylogenetically classified. Our data revealed a drastic reduction in the diversity of OTUs in raw milk during cold storage at 6°C at 97% similarity level; but, the N2-flushing treatment alleviated this reduction and substantially limited the loss of bacterial diversity during the same cold-storage period. Compared to cold-stored milk, the initial raw-milk samples contained less Proteobacteria (mainly Pseudomonadaceae, Moraxellaceae and Enterobacteriaceae) but more Firmicutes (mainly Ruminococcaceaea, Lachnospiraceae and Oscillospiraceaea) and Bacteroidetes (mainly Bacteroidales). Significant differences between cold-stored and additionally N2-flushed milk were mainly related to higher levels of Pseudomononadaceae (including the genera Pseudomonas and Acinetobacter) in cold-stored milk samples; furthermore, rare taxa were better preserved by the N2 gas flushing compared to the cold storage alone. No major changes in bacterial composition with time were found regarding the distribution of the major 9 OTUs, that dominated the Pseudomonas genus in N2-flushed or non-flushed milk samples, other than an intriguing predominance of bacteria related to P. veronii. Overall, this study established that neither bacteria causing milk spoilage nor any well-known human pathogen or anaerobe benefited from the N2 gas flushing even though the N2-flushed and non-flushed cold-stored milk differed in bacterial counts by up to 104-fold.

Degradation of Benzene by Pseudomonas veronii 1YdBTEX2 and 1YB2 Is Catalyzed by Enzymes Encoded in Distinct Catabolism Gene Clusters

Pseudomonas veronii 1YdBTEX2, a benzene and toluene degrader, and Pseudomonas veronii 1YB2, a benzene degrader, have previously been shown to be key players in a benzene-contaminated site. These strains harbor unique catabolic pathways for the degradation of benzene comprising a gene cluster encoding an isopropylbenzene dioxygenase where genes encoding downstream enzymes were interrupted by stop codons. Extradiol dioxygenases were recruited from gene clusters comprising genes encoding a 2-hydroxymuconic semialdehyde dehydrogenase necessary for benzene degradation but typically absent from isopropylbenzene dioxygenase-encoding gene clusters. The benzene dihydrodiol dehydrogenase-encoding gene was not clustered with any other aromatic degradation genes, and the encoded protein was only distantly related to dehydrogenases of aromatic degradation pathways. The involvement of the different gene clusters in the degradation pathways was suggested by real-time quantitative reverse transcription PCR.

Flow cytometry and real-time quantitative PCR as tools for assessing plasmid persistence

The maintenance of a plasmid in the absence of selection for plasmid-borne genes is not guaranteed. However, plasmid persistence can evolve under selective conditions. Studying the molecular mechanisms behind the evolution of plasmid persistence is key to understanding how plasmids are maintained under nonselective conditions. Given the current crisis of rapid antibiotic resistance spread by multidrug resistance plasmids, this insight is of high medical relevance. The conventional method for monitoring plasmid persistence (i.e., the fraction of plasmid-containing cells in a population over time) is based on cultivation and involves differentiating colonies of plasmid-containing and plasmid-free cells on agar plates. However, this technique is time-consuming and does not easily lend itself to high-throughput applications. Here, we present flow cytometry (FCM) and real-time quantitative PCR (qPCR) as alternative tools for monitoring plasmid persistence. For this, we measured the persistence of a model plasmid, pB10::gfp, in three Pseudomonas hosts and in known mixtures of plasmid-containing and -free cells. We also compared three performance criteria: dynamic range, resolution, and variance. Although not without exceptions, both techniques generated estimates of overall plasmid loss rates that were rather similar to those generated by the conventional plate count (PC) method. They also were able to resolve differences in loss rates between artificial plasmid persistence assays. Finally, we briefly discuss the advantages and disadvantages for each technique and conclude that, overall, both FCM and real-time qPCR are suitable alternatives to cultivation-based methods for routine measurement of plasmid persistence, thereby opening avenues for high-throughput analyses.