Dynamic Response of Pseudomonas putida S12 to Sudden Addition of Toluene and the Potential Role of the Solvent Tolerance Gene trgI

Pseudomonas putida S12 is exceptionally tolerant to various organic solvents. To obtain further insight into this bacterium’s primary defence mechanisms towards these potentially harmful substances, we studied its genome wide transcriptional response to sudden addition of toluene. Global gene expression profiles were monitored for 30 minutes after toluene addition. During toluene exposure, high oxygen-affinity cytochrome c oxidase is specifically expressed to provide for an adequate proton gradient supporting solvent efflux mechanisms. Concomitantly, the glyoxylate bypass route was up-regulated, to repair an apparent toluene stress-induced redox imbalance. A knock-out mutant of trgI, a recently identified toluene-repressed gene, was investigated in order to identify TrgI function. Remarkably, upon addition of toluene the number of differentially expressed genes initially was much lower in the trgI-mutant than in the wild-type strain. This suggested that after deletion of trgI cells were better prepared for sudden organic solvent stress. Before, as well as after, addition of toluene many genes of highly diverse functions were differentially expressed in trgI-mutant cells as compared to wild-type cells. This led to the hypothesis that TrgI may not only be involved in the modulation of solvent-elicited responses but in addition may affect basal expression levels of large groups of genes.


Introduction
Pseudomonas putida S12 is an exceptionally solvent-tolerant bacterium that thrives in the presence of organic solvents such as 1-octanol, toluene and benzene and as such is important for the production of industrially relevant chemicals [1][2][3][4].Several mechanisms of solvent tolerance have been identified in this organism, the most important of which is the solvent extrusion pump SrpABC that is located in the cytoplasmic membrane.This RND-type transporter actively removes solvent molecules from the membrane at the expense of the proton gradient [5,6].
Previously, we reported on the proteome and transcriptome responses of P. putida S12 to toluene [7,8].Several of the observed responses confirmed previous findings, e.g., the induction of the solvent extrusion pump SrpABC and differential expression of membrane-associated and stress-response genes.A novel finding was the differential expression of energy-management systems.This response may compensate for the dissipation of the proton gradient associated with solvent stress, since the accumulation of solvent molecules in the cytoplasmic membrane brings about increased permeability for protons (uncoupling effect).In addition, SrpABC-mediated solvent extrusion draws on the proton gradient.Moreover, a novel gene involved in solvent tolerance, trgI, was identified.The expression of trgI decreased immediately and very rapidly after sudden addition of toluene.In steady-state chemostats, i.e, under fully toluene-adapted conditions, expression of both the gene and its encoded protein was significantly down-regulated in the presence of 5 mM toluene [7,8].Deletion of trgI [7] resulted in a more rounded cell morphology and increased resistance to sudden toluene shock.Although the precise function of trgI remained obscure, we hypothesise that it is involved in the first line of defence against solvents [7].
The immediate down-regulation of trgI upon toluene exposure provoked the question how expression of other genes in P. putida S12 would respond to the sudden addition of toluene.This would provide valuable insights into the cellular functions involved in the early adaptation response to organic solvents.Although previously several-omics studies of solvent-exposed microorganisms have been published (for example [7][8][9][10][11][12]), none of these involved a genomewide monitoring of the early adaptation response.Instead, batch cultures were sampled at a single time-point, or steady-state chemostat cultures that were fully adapted to organic solvent were analysed.Therefore, we chose to study the global gene expression profiles of P. putida S12 during the first 30 minutes after the addition of toluene.In addition to wild type P. putida S12, the trgI knock-out mutant P. putida S12ΔtrgI was investigated to shed more light onto the role of this gene in the early solvent tolerance response.

Bacterial strains
The bacterial strains used in this study are Pseudomonas putida S12, which was originally isolated as a styrene utilising bacterium [13] and P. putida S12ΔtrgI.P. putida S12ΔtrgI is a trgI knock-out mutant that was constructed as described previously [7].

Standard culture conditions
Luria-Bertani broth (LB medium) [14] was used as the standard culturing medium.As a solid medium, LB with 1.5% (w/v) agar was used.Batch cultivation was routinely carried out in 100-ml Erlenmeyer flasks containing 25 ml of liquid medium, placed on a horizontally shaking incubator at 30°C.

Analysis of differential gene expression after a sudden addition of toluene
Culture conditions and sampling.Differential gene expression after a sudden addition of 5 mM toluene was analysed in early exponential phase cultures (optical density at 600 nm of 0.5-0.6) of P. putida strains S12 (wild-type) and S12ΔtrgI.The cells were grown in 100 ml of LB medium in 1-L bottles placed on a horizontally shaking incubator at 30°C.Samples of 1 ml were drawn with a syringe through the rubber-covered hole in the screw cap of the bottle immediately before (t = 0) and at set intervals (1, 2, 5, 10 and 30 minutes) after toluene addition.See [7] for details about sampling.
Transcriptome analysis.mRNA isolation, cDNA preparation and hybridisation for transcriptome analysis were performed as described previously [7].The custom made high-density microarrays used were based on the genome sequence of P. putida S12 (GenBank accession AYKV00000000.1).The microarray results, as well as the microarray itself, were made public in Gene Expression Omnibus, accession GSE64791.
Data analysis.Microarray data were imported into the GeneSpring GX 7.3.1 software package (Agilent Technologies, Santa Clara, CA, USA) using the GC RMA algorithm.After normalisation (signals below 0.01 were taken as 0.01; per chip: normalise to 50 th percentile; per gene: normalise to median) of the data, probe sets representing control genes were removed, as well as absolute non-changing loci (between 0.667-and 1.334-fold change).The resulting set of 6164 differentially expressed loci was used for further investigation.
The overall transcriptional activity change during toluene exposure was quantified by calculating the total number of differentially expressed genes for each time point after toluene addition, using the transcript levels at t = 0 as reference (Fig 1).
Overrepresentation and statistical analysis.Overrepresentation of specific groups of genes among the total response-groups was determined using the hyper-geometric test in the R statistical program environment.For example: in a group of genes selected for up-regulation following toluene exposure, the chance that an x-number of genes of COG group A are in that group of up-regulated genes is being assessed.Or, in other words, it can be assessed what the chance is that this group of up-regulated genes contains 50% COG group A while all genes contain 30% COG group A. The comparison between sets of two such groups (in this example 'upregulated' and 'COG group A') can be found in Table 1.

Results
The dynamic response of P. putida S12 to sudden toluene exposureglobal analysis of transcriptomics data The addition of toluene to exponentially growing P. putida S12 provoked an extensive change in the global transcriptional activity (Fig 1).Although this effect was observed in both the wildtype and the trgI-knockout strain, marked differences were observed in the timing and the total number of differentially expressed genes.In wild-type P. putida S12, the abundance of differentially expressed genes increased immediately after toluene addition, reaching a maximum after 5 minutes, after which it decreased and more or less stabilized 10 min.after toluene addition.In P. putida S12ΔtrgI, the expression response appeared to be delayed.The abundance of differentially expressed genes was initially much lower and only reached a maximum 10 minutes after toluene addition.Subsequently, it stabilized to a level above that in the wildtype strain.
Furthermore, in Fig 1 the abundance of the genes is presented that were differentially expressed in both the wild-type strain and P. putida S12ΔtrgI.This group likely represents genes responding to toluene irrespective of the presence of trgI.The abundance profile is remarkably similar to that of the differentially expressed genes in the trgI-knockout strain (Fig 1).The peak in the abundance profile of the wild-type strain 5 minutes after toluene addition is indicative of a generic rearrangement response that (partly) involves TrgI.
The fold changes of all genes whose expression altered after toluene addition were calculated for every possible combination of two time-points, i.e., between 0 and 1 minutes, 0 and 2 minutes, 5 and 30 minutes, etc.These combinations represent all possible time-frames within the experiment.The-apparently toluene-responsive-genes were furthermore classified by functional category to determine whether specific cellular functions were overrepresented within certain time frames.Functional classification was performed by the categories defined in the NCBI COG database (Clusters of Orthologous Groups of proteins; ftp://ftp.ncbi.nih.gov/pub/COG/COG/[15,16]).Fig 2 shows a summary of the results for both wild-type P. putida S12 and the trgI-knockout strain; a complete overview is shown in S1 Fig.
Overrepresentation of specific functional categories was most explicit in wild-type P. putida S12, for genes up-regulated immediately after addition of toluene.Among these genes, category C (Energy production and conversion), E (Amino acid transport and metabolism), I (Lipid metabolism) and O (Posttranslational modification, protein turnover and chaperones) were significantly overrepresented.In P. putida S12ΔtrgI only category O was overrepresented among the genes that responded immediately to toluene exposure, consistent with the observed delayed global transcriptional response of S12ΔtrgI.Category T (signal transduction) was overrepresented among the genes that were first upregulated and then down-regulated in response to toluene, in both wild-type P. putida S12 and S12ΔtrgI.In the wild-type strain, also category N (motility) was overrepresented among these genes.In S12ΔtrgI, category N was overrepresented among the genes that were down-regulated later after toluene addition.
Other categories overrepresented in both strains among this group of genes that were down-regulated later were J (translation) and H (coenzyme metabolism).Genes of category F (nucleotide transport and metabolism) were overrepresented in both strains early after toluene addition.In the wild-type strain, category K (transcription) genes were down-regulated a few minutes after addition of toluene, followed by up-regulation.In P. putida S12ΔtrgI, genes of this category were among the genes that were consistently up-regulated after toluene addition.Thus, clear differences were visible between wild-type S12 and S12ΔtrgI, both in the global transcriptional response to toluene and in the timing of expression of genes from different functional categories.

Effect of trgI deletion on the transcriptional response to toluene
The analysis of the global transcriptional response to toluene revealed clear differences between wild-type P. putida S12 and strain S12ΔtrgI.To establish which individual genes were affected by the trgI deletion, each gene was classified by expression behaviour (Table 1).The following criteria were applied: 1) expression response to toluene in wild-type P. putida S12; 2) expression level prior to toluene exposure; 3) expression response to toluene in P. putida S12ΔtrgI.Based on these criteria, four classes of genes can be discriminated by the effect of the trgI deletion on the intrinsic expression level and/or the response to toluene.Furthermore, a functional category (see above) was coupled to each individual gene to assess whether specific cellular functions were more affected than others by the trgI deletion (Table 1).
Only 1287 genes, i.e., around 20% of the total genome, were apparently unaffected by the trgI deletion with regard to both the intrinsic expression level and the response to toluene (group 1 in Table 1).Among these genes, particularly COGs O (posttranslational modification, protein turnover and chaperones), H (coenzyme metabolism) and D (cell division and chromosome partitioning) were overrepresented.This observation suggests that the trgI-deletion did not affect basic cell processes to a large extent.Also the number of genes that showed an unaltered response to toluene but an altered intrinsic expression level was relatively small (<500, group 2 in Table 1).Among these genes, COGs C (energy production and conversion) and O were overrepresented, indicating that the trgI deletion caused an intrinsic change in the energy metabolism, as well as in the chaperone machinery.
The large majority of the genes that were affected by the trgI-deletion showed an aberrant toluene response.These genes comprise a total of more than 4300 genes.More than 1700 genes appeared to have gained toluene responsiveness due to the trgI deletion (group 3 in Table 1).However, an even larger number (> 2600) appeared to have lost this property (group 4 in Table 1).The majority of this group (1594 genes) in addition showed an altered intrinsic expression level which may imply that these genes are expressed or repressed in a constitutive manner in P. putida S12ΔtrgI.Among the overrepresented functional groups in group 4, COG C can be linked to solvent tolerance which underlines the contribution of energy management systems to solvent tolerance.

Effect of toluene and deletion of trgI on genes and pathways associated with toluene tolerance
As indicated above, deletion of trgI affected the expression of a very large number of genes, suggesting that trgI exerted a high-level regulatory effect.This regulatory effect was furthermore strongly associated with the response to toluene.Although the trgI-affected genes covered many different functional groups, COGs C, E and I (energy, amino acid metabolism and lipid metabolism) appeared to be overrepresented (Table 1 and Fig 2).Since these functional categories can be rationally linked to solvent tolerance, the transcriptional response of genes and pathways within these COGs was studied in more detail, together with genes and pathways with an established role in solvent tolerance (also see S1 Table ).

Toluene-responsive genes involved in energy metabolism
The accumulation of organic solvents in the bacterial membrane puts a heavy burden on the cellular energy metabolism as a result of uncoupling and the energy requirement of the solvent extrusion pump.This was clearly reflected in the differential expression of many genes associated with cellular energy management systems, most notably cytochrome c oxidase-encoding genes (Fig 3).P. putida S12 harbours several different types of terminal cytochrome c oxidases: in addition to the high-oxygen affinity types cbb 3 -1 and cbb 3 -2 there is also the low affinity variant aa 3 [17].
The genes encoding the aa 3 -type were consistently expressed to a lower level than the genes encoding the cbb 3 -type, irrespective of the presence of toluene or trgI.The two cbb 3 variants were expressed to different levels depending on the presence of toluene.In S12ΔtrgI, the intrinsic expression level of both cbb 3 type-encoding genes was higher than in the wild-type strain.In the presence of toluene, the expression level of the cbb 3 -1 genes increased even further and subsequently decreased to the wild-type level.The slightly variable expression levels of the cbb 3 -2 genes consistently exceeded those of wild-type P. putida S12 upon toluene exposure.
Along with cytochrome c oxidase-encoding genes, several genes encoding subunits of the NADH dehydrogenase complex were expressed to a higher level in S12ΔtrgI within 10 minutes after toluene exposure (Fig 4).After 30 minutes, the expression level of the genes encoding the quinone oxidoreductase-subunit had returned to the wild-type level whereas the genes coding for the dehydrogenase subunit remained at a high level for at least 30 minutes after toluene exposure.

Toluene-responsive genes involved in lipid metabolism
The accumulation of organic solvent molecules in the bacterial membrane leads to cell lysis, releasing membrane fatty acids into the medium.In E. coli free fatty acids induce fatty acid degradation via fadR [18].P. putida S12 does not harbour a fadR homologue, but the gene psrA is responsible for the regulation of part of the fatty acid degradation pathway [18,19].In tolueneexposed P. putida S12 psrA as well as the genes of the fatty acid degradation pathway were upregulated after addition of toluene (S2 Fig, solid lines).In addition to fatty acid degradation, fatty acid synthesis is expected to occur in solvent-exposed P. putida S12 since newly synthesized fatty acids are required to counteract the toluene-induced degradation effects [20].Differential expression of genes involved in fatty acid biosynthesis, however, was variable (S3 Fig) .In wild-type P. putida S12, COG I genes were overrepresented after addition of toluene.In S12ΔtrgI, this was not the case (Fig 2 ) indicating less toluene-induced membrane damage in this strain, or a lower need for de-novo membrane fatty acid synthesis under toluene stress.It should be noted that the fatty acid biosynthetic genes are poorly annotated in the P. putida S12 genome that was used for the custom micro array, which prevented a more detailed analysis of the responses of these genes to toluene.

Toluene-responsive genes involved in acetyl-CoA metabolism
The genes encoding isocitrate lyase and malate synthase were found to be up-regulated in response to toluene.This effect was observed in both wild-type S12 and S12ΔtrgI (S4 Fig) although the level of up-regulation was higher in the latter.The encoded enzymes constitute the glyoxylate bypass, which converts acetyl-CoA into succinate and malate for biosynthesis of, a.o., amino acids of the aspartate family.This pathway is normally active during growth on C 2compounds and is controlled at the enzymatic level through NADH-mediated inhibition of isocitrate dehydrogenase [21].Generally, induction of the glyoxylate bypass genes and bypass activity are correlated [22], suggesting that NADH is present at surplus amounts during toluene stress.This may be unexpected in view of the high energy demand brought about by toluene stress, but it may also be an effect caused by malfunctioning respiratory chain components.It is well known that the functioning of membrane-embedded proteins is affected by accumulation of solvent molecules in the cytoplasmic membrane, but also due to compositional changes of the membrane in response to solvents [23].Possibly, the differential expression of respiratory chain components in response to toluene (as described above) is a compensatory response to compromised electron transport chain functioning.Alternatively, the glyoxylate pathway may be co-regulated with the degradation of fatty acids, or with valine, leucine and isoleucine degradation.In these pathways, the main degradation product is acetyl-CoA which effectively represents a C 2 substrate [18].It should be noted that multiple genes involved in the degradation of these amino acids were up-regulated in response to toluene and, again, the level of up-regulation was significantly higher in S12ΔtrgI (S5 Fig).

Toluene-responsive genes involved in solvent extrusion
Expression of srpABC, coding for the solvent extrusion pump SrpABC, and srpRS, the associated regulatory genes, increased immediately in wild type P. putida S12 after addition of toluene.This is in accordance with previous observations [7], although the expression level appeared to stabilize after 5 min. of toluene exposure rather than reaching a maximum after 30 min.as reported earlier [24] (Fig 5).Expression of other loci, encoding part of the putative toluene transport system Ttg2FEDC and the multidrug transporter ArpABC, increased immediately after addition of toluene, reaching maximum levels after 10 to 30 minutes of toluene exposure (Fig 5).All expression levels were similar in wild-type and S12ΔtrgI prior to toluene exposure, but increased much faster in the mutant, to much higher levels.It is unclear to which extent Ttg2FEDC is involved in toluene tolerance, and ArpABC has been established to play a role in antibiotic resistance rather than solvent tolerance [25].These observations indicate that trgI is involved in modulation of general responses to stress conditions including, but not restricted to, solvent-associated stress.These observations clearly indicate that, next to solventassociated stress, trgI is involved in modulation of other stress responses as well.Effect of toluene and deletion of trgI on genes and pathways that are not directly associated with toluene tolerance Differential expression of the arginine deiminase pathway in response to toluene.The ubiquitous arginine deiminase (ADI) pathway serves to generate energy from arginine fermentation under anoxic or carbon-limiting conditions [26].In the ADI pathway, arginine degradation via citrulline to ornithine yields 1 ATP per arginine molecule.Each ornithine generated is exported and exchanged for another arginine molecule via an arginine-ornithine specific antiporter.In wild-type P. putida S12, the genes encoding the enzymes for this pathway (arginine deiminase, ornithine carbamoyltransferase and carbamate kinase) and the arginine/ornithine antiporter were up-regulated immediately after addition of toluene (S6 Fig) .In contrast, another associated arginine/ornithine antiporter gene appeared to be down-regulated, together with aotM, aotP and aotQ, that constitute an arginine/ornithine importer.In P. putida S12ΔtrgI the expression of the ADI pathway genes was considerably higher; already in the absence of toluene, the expression was similar to that of the wild type after 30 minutes of toluene exposure.Upon addition of toluene, the expression levels increased even further.The expression of the arginine/ornithine antiporter gene decreased in strain S12ΔtrgI after an initial up-regulation, as did the expression of aotM, aotP and aotQ.Since the ADI pathway genes were up-regulated while the antiporter genes were down-regulated, it appears that the ADI pathway rather serves to accumulate intracellular ornithine than to generate energy under toluene-stressed conditions.Differential expression of protein folding genes in response to toluene.In the periplasm, proteins can be folded oxidatively via disulphide bond formation.The enzymes involved in this process are DsbAB, which oxidizes unfolded proteins and DsbCDG, which isomerizes disulfide bonds in misfolded proteins [27].In P. putida S12, dsbACDG have been annotated, as well as tlpA that likely also belongs to this system.In addition, chaperonins like DnaJ, DnaK, GroES, GroEL and HtpG play an important role in protein folding.All genes involved in this protein-damage repair system as well as a gene encoding a 33-kDa chaperonin were down-regulated in the wild-type strain within the first two minutes of toluene exposure, and subsequently up-regulated (S7 Fig and S8 Fig).In S12ΔtrgI, most if not all genes were up-regulated to an even higher extent.Hence, the protein damage-repair response is induced even more strongly in strain S12ΔtrgI during solvent exposure.

Deletion of trgI affects glucose and fructose metabolism
A remarkable physiological effect of the trgI deletion is loss of the ability to utilize glucose or fructose as the sole source of carbon and energy [6].Whereas fructose is imported via a PTStype transporter and further metabolized in the cytoplasm as fructose-1-phosphate, the initial glucose metabolism is more complex in Pseudomonads.Glucose enters the periplasm via porins OprB-1 or OprB-2 and is subsequently transported directly into the cytoplasm via an ABC transporter encoded by gtsABCD, or oxidized in the periplasm via gluconate (catalyzed by Gcd, glucose dehydrogenase) to 2-ketogluconate (catalysed by Gad) [28].Gluconate and 2-ketogluconate are transported to the cytoplasm via dedicated transporters encoded by gntP, respectively kguT and phosphorylated by GnuK, respectively KguK.
The expression of the fructose transporter genes was similar in wild-type P. putida S12 and in strain S12ΔtrgI whereas transcript levels of oprB-1, oprB-2 and gtsABCD were actually higher in S12ΔtrgI ( Fig 6).Gcd was expressed at a slightly lower level whereas gad, gntP, gnuK, kguT and kguK were not detectably transcribed.The expression of these genes is known to be induced by gluconate and ketogluconate [29,30].Therefore, the lack of detectable expression suggested that gluconate and ketogluconate could not be formed.Since both gcd and gtsABCD (as well as the fructose transporter genes) were transcribed, impaired transport of glucose (and fructose) into the periplasm appears to be the only viable explanation for the lack of growth on sugars.This phenomenon may be caused by structural modification of the outer cell structure, which may affect the functioning of the membrane-embedded transport proteins.A similar explanation has been reported for the impaired functioning of drugs efflux transporters of a mutant P. putida DOT-T1E strain with altered phospholipid head group composition [31].

Discussion
TrgI is a previously identified toluene-responsive gene of P. putida S12 with unknown function [7].To gain further insight into the solvent tolerance response mechanism of Pseudomonas and into a possible role of TrgI, the transcriptome responses of wild type S12 and of trgI-deletion mutant S12ΔtrgI [7] were monitored after addition of toluene.
In the first five minutes after toluene addition, the absolute number of differentially expressed genes was significantly smaller in a ΔtrgI strain than in a wild-type strain.Many genes in the mutant showed an altered intrinsic expression level.This suggested that TrgI is not only involved in toluene-elicited responses but also affects basal expression levels of large groups of genes.It appears that inactivation of trgI prepares the cells for a sudden addition of organic solvent by 'setting' basal expression of many genes to levels comparable to those in toluene-exposed wild type cells.However, the number of differentially expressed genes remained high in strain S12ΔtrgI.Thus, some feedback loop appears to be disturbed as well by the deletion of trgI.
NADH dehydrogenase-and cytochrome c oxidase genes were differentially expressed upon toluene exposure, as well as more structurally in the trgI-deletion mutant.This fully correlated with the proton gradient-dissipating effect of toluene [33,34].The relatively high expression of the cbb 3 -type cytochrome oxidase over the lower oxygen-affinity aa 3 -type in both wild-type and mutant strains suggested that in cultures exposed to toluene, high-affinity cytochrome c oxidase activity is specifically needed to provide for an adequate proton gradient.
Cells exposed to organic solvents are known to change their fatty acid composition to strengthen their membranes by increasing the degree of saturation by cis-trans isomerisation of unsaturated fatty acids or by newly synthesising long chain fatty acids [20,32,35].Eventually, accumulation of solvent molecules in the membranes results in cell lysis.Lysing cells release fatty acids that trigger expression of the genes encoding the fatty acid degradation pathway [33].However, we did not find evidence that fatty acid-synthesis was increased.This suggested that the capacity of the available fatty acid synthesis machinery was sufficient.
Genes encoding the Ttg2FEDC putative toluene transport system were rapidly up-regulated in response to toluene in P. putida S12.Previously, only the SrpABC solvent efflux pump was found to be specifically expressed in response to toluene [2,7].As SrpABC becomes fully expressed only after extended toluene exposure [25], Ttg2FEDC may cushion the first blow dealt by the sudden exposure to toluene.In P. putida S12ΔtrgI maximum expression levels of the Srp genes was reached much sooner, as with the Ttg2 genes and the genes coding for the antibiotic pump ArpABC.
Remarkably, the ADI (arginine deiminase) pathway appeared to be up-regulated in response to toluene, and to an even higher level in S12ΔtrgI than in the wildtype strain.Only one arginine-ornithine antiporter gene was down-regulated in the wild-type whereas both were down-regulated in S12ΔtrgI, suggesting that P. putida S12ΔtrgI accumulates the polyamine ornithine in response to toluene, providing additional protection against solvent stress [34].Moreover, ornithine may act as a precursor for proline, a compatible solute implicated in protection against water and solvent stress [36,37,38].Since the cultures were in mid-exponential phase, it is unlikely that either anoxic or carbon-limiting conditions occurred, which are the normal conditions for inducing the ADI pathway [26].
Genes involved in oxidative protein folding and several chaperonin-encoding genes were up-regulated in both wild-type P. putida S12 and S12ΔtrgI upon addition of toluene.Presumably, the up-regulation of those genes related to the occurrence of protein misfolding and random formation of disulfide bonds.These events may well occur in solvent-exposed cells as a result of oxidative stress caused by impaired respiration.Protein folding genes were up-regulated to a higher level in the trgI-knockout strain, suggesting this strain has a higher need of them.
Transcriptome analysis of P. putida S12 and S12ΔtrgI following sudden exposure to toluene has shed new light on solvent tolerance.The initial response to toluene in wild-type cells appeared to be fully geared towards survival before all solvent tolerance mechanisms are fully induced.The clear delay in timing of the overall transcriptional response in the trgI-deletion mutant was indicative for a milder response to sudden toluene exposure.This would be in agreement with an intrinsically improved toluene tolerance.The impact of the trgI deletion thus appeared to be relatively broad.
The diversity of responses associated with trgI, as well as its wide impact on gene expression, suggest an important regulatory role.Analyses using BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi) and SMART (http://smart.embl-heidelberg.de/)did not provide any clues to a potential function.A model of TrgI tertiary structure (http://zhanglab.ccmb.med.umich.edu/I-TASSER/) [39] showed highest structural similarity (TM-score of 0.646) with a molybdatedependent transcriptional regulator, modE, of E. coli [40], although TrgI is much smaller than ModE and does not exhibit clear DNA binding residues.Hence, one hypothesis may envision TrgI as a modulator of other transcription factors.S1 Table .Overview of all genes mentioned in main text.To make it easier to find all genes mentioned in the text in the GenBank and GEO databases, these genes are shown with their names, RPPX numbers, probe names, and additional remarks.(XLSX)

Fig 1 .
Fig 1. Summary of the transcriptomics results.Number of genes differentially expressed with an absolute log2(expression level at t = 0 / expression level at t = t) 0.5 in wild-type S12 (dashed line) and mutant S12ΔtrgI (solid line) at the indicated time points are shown.The dotted line shows the number of genes that is present in both comparisons.T = 0: no toluene present.doi:10.1371/journal.pone.0132416.g001

Table 1 .
Classification of genes by effect of trgI deletion on intrinsic expression level and toluene responsiveness.Intrinsic expression level is expression level at t = 0, immediately before addition of toluene.The expression level was defined as 'different' at a ratio of <-0.5 or >0.5.bToluene response is response immediately after addition of toluene until t = 30 min.The expression level was defined as 'different' when the correlation (of the expression pattern) between S12 and S12ΔtrgI was <0.8.c COGs were considered overrepresented if p<0.01For COG abbreviations, seeFig 2.