# Escherichia coli under Ionic Silver Stress: An Integrative Approach to Explore Transcriptional, Physiological and Biochemical Responses

Claire Saulou-Berion, Ignacio González, Brice Enjalbert, Jean-Nicolas

Audinot, Isabelle Fourquaux, Frédéric Jamme, Muriel Cocaign-Bousquet,

Muriel Mercier-Bonin, Laurence Girbal

► To cite this version:

Claire Saulou-Berion, Ignacio González, Brice Enjalbert, Jean-Nicolas Audinot, Isabelle Fourquaux, et al.. Escherichia coli under Ionic Silver Stress: An Integrative Approach to Explore Transcriptional, Physiological and Biochemical Responses. PLoS ONE, 2015, 12 (12), pp.e0145748. 10.1371/journal.pone.0145748. hal-01510596

HAL Id: hal-01510596

https://hal.science/hal-01510596

Submitted on 19 Apr 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# *Escherichia coli* under Ionic Silver Stress: An Integrative Approach to Explore Transcriptional, Physiological and Biochemical Responses

Claire Saulou-Bérion¹,²,³ na, Ignacio Gonzalez¹,²,³, Brice Enjalbert¹,²,³, Jean-Nicolas Audinot⁴, Isabelle Fourquaux⁵, Frédéric Jamme⁶,⁷, Muriel Cocaign-Bousquet¹,²,³, Muriel Mercier-Bonin¹,²,³nb, Laurence Girbal¹,²,³*

1 Université de Toulouse, INSA, UPS, INPT, LISBP, Toulouse, France, 2 INRA, UMR792 Ingénierie des Systèmes Biologiques et des Procédés, Toulouse, France, 3 CNRS, UMR5504, Toulouse, France,

4 Luxembourg Institute of Science and Technology (LIST), Material Research & Technology Department (MRT), Belvaux, Luxembourg, 5 Faculté de Médecine Rangueil, Centre de Microscopie Electronique Appliquée à la Biologie (CMEAB), Toulouse Cedex, France, 6 INRA, UAR1008, CEPIA, Nantes, France,

7 Synchrotron SOLEIL, Gif-sur-Yvette, France

<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

<figure><img src="image_3.png" /></figure>

## OPEN ACCESS

**Citation:** Saulou-Bérion C, Gonzalez I, Enjalbert B, Audinot J-N, Fourquaux I, Jamme F, et al. (2015) *Escherichia coli* under Ionic Silver Stress: An Integrative Approach to Explore Transcriptional, Physiological and Biochemical Responses. PLoS ONE 10(12): e0145748. doi:10.1371/journal.pone.0145748

**Editor:** Ying-Jan Wang, National Cheng Kung University, TAIWAN

**Received:** July 10, 2015

**Accepted:** December 8, 2015

**Published:** December 22, 2015

**Copyright:** © 2015 Saulou-Bérion et al. This is an open access article distributed under the terms of the [Creative Commons Attribution License](http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All raw and processed transcriptomic data were deposited on Gene Expression Omnibus data repository (GSE67735).

**Funding:** Funding was received from the French National Research Agency (http://www.agence-nationale-recherche.fr/), grant number: ANR-07-BLAN-0196-01 to CSB. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## Abstract

For a better understanding of the systemic effect of sub-lethal micromolar concentrations of ionic silver on *Escherichia coli*, we performed a multi-level characterization of cells under Ag⁺-mediated stress using an integrative biology approach combining physiological, biochemical and transcriptomic data. Physiological parameters, namely bacterial growth and survival after Ag⁺ exposure, were first quantified and related to the accumulation of intracellular silver, probed for the first time by nano secondary ion mass spectroscopy at sub-micrometer lateral resolution. Modifications in *E. coli* biochemical composition were evaluated under Ag⁺-mediated stress by *in situ* synchrotron Fourier-transform infrared micro-spectroscopy and a comprehensive transcriptome response was also determined. Using multivariate statistics, correlations between the physiological parameters, the extracellular concentration of AgNO₃ and the intracellular silver content, gene expression profiles and micro-spectroscopic data were investigated. We identified Ag⁺-dependent regulation of gene expression required for growth (e.g. transporter genes, transcriptional regulators, ribosomal proteins), for ionic silver transport and detoxification (e.g. *copA*, *cueO*, *mgtA*, *nhaR*) and for coping with various types of stress (*dnaK*, *pspA*, *metA*, *R*, oxidoreductase genes). The silver-induced shortening of the acyl chain of fatty acids, mostly encountered in cell membrane, was highlighted by microspectroscopy and correlated with the down-regulated expression of genes involved in fatty acid transport (*fadL*) and synthesis/modification of lipid A (*lpxA* and *arnA*). The increase in the disordered secondary structure of proteins in the presence of Ag⁺ was assessed through the conformational shift shown for amides I and II,

<footnote>na Current address: INRA, UMR782 Génie et Microbiologie des Procédés Alimentaires, AgroParisTech, Thiverval-Grignon, France</footnote>

<footnote>nb Current address: Toxalim UMR INRA/INPT/UPS 1331, Toulouse, France</footnote>

<footnote>* Laurence.girbal@insa-toulouse.fr</footnote>
**Competing Interests:** The authors have declared that no competing interests exist.

and further correlated with the up-regulated expression of peptidase (*hfq*) and chaperone (*dnaJ*), and regulation of transpeptidase expression (*ycfS* and *ycbB*). Interestingly, as these transpeptidases act on the structural integrity of the cell wall, regulation of their expression may explain the morphological damage reported under Ag⁺-mediated stress. This result clearly demonstrates that the cell membrane is a key target of ionic silver.

Introduction

Ionic silver (Ag⁺) has been reported to exhibit strong biocidal activity in a wide range of bacteria and fungi, even at low micromolar concentrations [1–6]. Since ancient times, silver and its derived compounds have been empirically exploited for their therapeutic effect on clinical diseases (burn wounds, skin infections, newborn eye gonorrhea [7, 8]). In recent decades, silver derivatives have been used for new medical applications (e.g., silver sulfadiazine cream, silver-coated nylon fibers, dressings or catheters), notably against antibiotic-resistant bacteria [4].

The precise mechanisms of ionic silver antibacterial action are not fully understood. Several modes of action have been proposed, the most likely being the interaction of the Ag⁺ ions with electron-donating groups such as thiols, carboxylates, amides, imidazoles, indoles and hydroxyls [1, 9–11]. This type of chemical interaction could target membrane-bound and cytoplasmic proteins and enzymes, as well as bases in DNA. These interactions have been reported to cause (i) structural changes in the bacterial membranes [1], (ii) inhibition of the enzymes of the respiratory chain, and consequently, decoupling of respiration from ATP synthesis [12], (iii) formation of reactive oxygen species [13, 14], and (iv) lesions in DNA, affecting chromosome replication ability [1]. A recent study by our group was performed on the model bacterium *Escherichia coli* by synchrotron Fourier-transform infrared (sFTIR) microspectroscopy to elucidate the biochemical effects of ionic silver at single-cell scale [15]. These results showed that the observed biochemical changes were closely dependent on the concentration of Ag⁺ (lethal vs. sub-lethal concentration), growth conditions (stress applied from cell inoculation on, or as a pulse in the mid-exponential phase), and in some conditions, were associated with morphological damage. Indeed, the inhibitory action of silver ions was shown to be targeted against lipids and, at high concentrations, against proteins through a loss in α-helix conformation. The heterogeneity of Ag⁺-mediated effects within the bacterial population was also highlighted. In another study performed at the level of gene expression [16], the effect of silver ions on *E. coli* cells was compared to that of copper ions, highlighting shared up-regulation of genes coding for efflux pumps (i.e., up-regulation of the genes coding for CopA, a P-type ATPase efflux pump, and for CusCFBA, a Cu⁺ efflux pump) and a detoxifying oxidase (i.e., up-regulation of the gene coding for the CueO oxidase). More recently, the response of *E. coli* to nano-sized silver was investigated using the same type of transcriptomic approach [17]. This study provided evidence that exposure to this form of silver stimulates anaerobic respiration (i.e., up-regulation of anaerobic respiration-related reductases) and induces up-regulation of copper/silver resistance genes, including *cusBCF*, *copA* and *cueO*. However, a comprehensive understanding of the effect of ionic silver on *E. coli* bacteria is still lacking.

In this context, following our previously described results [15], the present work was devoted to multi-level characterization of *E. coli* exposed to Ag⁺-mediated stress in which an integrative biology approach, based on the combination of physiological, biochemical and transcriptomic data sets, was used for the first time. Bacterial growth and survival after Ag⁺ exposure were first quantified and related to the accumulation of intracellular silver, detected by
secondary ion mass spectroscopy at high lateral resolution (NanoSIMS) [18–20]. The whole transcriptomic response of *E. coli* cells under ionic silver-mediated stress was then characterized. Based on the modifications in *E. coli* biochemical composition observed by *in situ* synchrotron Fourier-transform infrared microspectroscopy and depicted in detail in [15], clear correlations were established between (i) variations in the cellular physiological parameters, (ii) variations in the biochemical characteristics of cell fatty acids and proteins, and (iii) regulation of gene expression. This challenging approach is shown in the present study to be powerful for determining key genetic markers of the *E. coli* response to ionic silver. The expression of genes involved in general growth-related functions was globally down-regulated whereas that of genes related to stress responses (metal ion stress, general stress and oxidative stress) was up-regulated. Regulation of gene expression after Ag⁺ exposure was also correlated with modifications in *E. coli*, including cell morphology and biochemical composition, such as acyl chain length and protein secondary structure.

## Materials and Methods

### Bacterial strain and growth conditions

The laboratory strain *E. coli* K12 MG1655 (genotype: F- lambda- ilvG- rfb-50 rph-1 [21] was used in this work. As previously described in [15], a first inoculum was obtained by adding 1 mL of stock culture to 5 mL of Luria Bertani broth (Biokar Diagnostics, PRS Panreac) in a 15-mL half-opened BD Falcon™ tube, followed by 8 h of incubation at 37°C under continuous shaking at 120 rpm. This preculture was then used to inoculate 40 mL of M9 medium (Merck Prolabo) complemented with 2.5 g/L of anhydrous glucose (GPR Rectapur) with approximately 1.5x10⁷ CFU/mL (corresponding to an optical density at 595 nm [OD₅₀₀ nm] of 0.02). This second preculture was incubated overnight under the same conditions (37°C, 120 rpm), then used to inoculate cultures to test the effect of ionic silver. The AgNO₃ concentrations used in all the following experiments described here were 0, 5.0, 6.5 and 8.5 µM. These cultures were performed at 37°C and 120 rpm in 200 mL of M9 broth complemented with 2.5 g/L of anhydrous glucose, inoculated at 10⁸ CFU/mL (corresponding to an OD₅₀₀ nm of 0.12).

### Ionic silver treatment of *E. coli* cultures

Ionic silver stock solutions were prepared from AgNO₃ powder (Fisher Scientific Bioblock) as previously described in [15] and stored in the dark. The final tested concentrations were obtained by diluting the corresponding AgNO₃ stock solution directly into the M9 growth medium. A control experiment was simultaneously performed by adding deionised water to the M9 broth. In both cases, a stabilization period of 24 h at 37°C and 120 rpm was included before cell inoculation. Bacterial growth kinetics were then monitored by measuring optical density at 595 nm (OD₅₀₀ nm). This made it possible to calculate the specific growth rate µ (h⁻¹) throughout the course of the experiment, notably the rate reached at sampling time (i.e., after 3 h of culture; µt = 3h) and that of the entire exponential phase of growth (µexpo). Experiments for each condition were performed in triplicate.

Samples were taken after 3 h of culture (corresponding to the mid-exponential phase for the growth without silver) for each condition. Cells were harvested by centrifugation with a Hettich Universal 320R centrifuge at 4,000 rpm for 5 min at 4°C, and washed twice with sterile saline solution (0.15 M NaCl). After incubation of cells for 24 h at 37°C, the cultivable cells were enumerated by plating on LB agar medium and expressed as colony-forming units per millilitre (CFU/mL). To account for variations in the number of cells recovered from each culture, the results were normalised by the total cell concentration and expressed in CFU/mL/ODU. Reported data correspond to the mean ± SD values of duplicates from three independent experiments.
## Biochemical composition of *E. coli* cells subjected to treatments with different concentrations of ionic silver

Modifications in the biochemical composition of *E. coli* cells grown in the presence of the different concentrations of ionic silver under study were evaluated at single-cell scale by *in situ* synchrotron FTIR micro-spectroscopy using the SMIS beamline at the SOLEIL Synchrotron (Gif-sur-Yvette, France), as previously described in detail in [15]. Briefly, the spatial resolution reached the micrometre order (4.1 µm), thanks to the use of a zinc-selenide attenuated total reflectance hemispherical internal reflection element, and optimized signal detection allowed by the use of the synchrotron source. Twenty spectra were acquired for each condition tested. The main informative regions observed were (i) the region 3,100–2,700 cm⁻¹ representative of fatty acids mostly encountered in cell membranes; (ii) the region 1,800–1,480 cm⁻¹ assigned to C = O stretching vibration, and amide I and amide II groups characteristic of the peptide bonds of cell proteins; and (iii) the region 1,450–1,150 cm⁻¹ attributed to PO₂⁻ groups mostly found in nucleic acids and, in some cases, to C-O stretching modes. To improve the detection of shifts in the characteristic absorption bands, the sFTIR spectra were pre-processed by calculating the second derivative (Savitsky–Golay, order 3, 9 points) on the whole infrared range (4,000–800 cm⁻¹).

## Elemental composition mapmaking of *E. coli* cells subjected to ionic silver treatments

Elemental maps of silver were obtained by NanoSIMS50 analysis from semi-thin sections of resin-embedded bacterial cells, by adapting an experimental protocol previously used for yeast cells [22]. Washed cells sampled after 3 h of culture were fixed for 1 h at 4°C using a 2.5% glutaraldehyde solution in cacodylate buffer (0.1 M, pH = 7.4), then harvested by centrifugation at 1,000 rpm for 10 min and embedded in a 1.5% agar solution. A second fixation step was performed for 16 h at 4°C with the same solution, followed by three washings in a cacodylate buffer (0.2 M, pH = 7.4). The post-fixation step was performed in a solution of osmium tetroxide (2%) in cacodylate buffer (0.2 M, pH = 7.4) for 2 h at room temperature, followed by three washings as previously described. Cell dehydration was achieved through a series of incubations in graded ethanol solutions (10 min each in 50%, 70% and 80%, then 2×15 min each in 90% and 95% ethanol). Infiltration and embedding were performed by incubating cells in several changes of LRWhite resin/ethanol mixtures, in which the concentration of resin was gradually increased to 100% (resin/ethanol 1:3, 1:2, 2:3 and pure resin, each bath for 12 h at 4°C). Samples were left to polymerize for 24 h at 60°C and finally sectioned into 300-nm thick slices (Ultracut Reichert) and mounted on silicon plots (Siltronix, Archamps, France).

Nano-secondary ion mass spectrometry (Nano-SIMS) analyses were then performed with a NanoSIMS50 (Cameca, Gennevilliers, France). The Cs⁺ primary beam was accelerated at 8 keV and used for rasterizing a 10 x 10 µm² area on the polarized sample surface (– 8 keV). The primary ion intensity was 1.0–0.8 pA corresponding to a probe working diameter in the range of 50 to 80 nm (lateral resolution; [19]) and leading to the erosion of 1–2 nm of matter. The masses simultaneously recorded using a multi-detection system were ¹²C¹⁴N⁻ (m = 26.00307), ³¹P⁻ (m = 30.97376), ³²S⁻ (m = 31.97207) and ¹⁰⁷Ag⁻ (m = 106.905095), whose intensities were expressed in counts per second (c/s). The distribution maps were recorded in a pixel format of 256 x 256 image points, with an acquisition time of 40 ms per pixel. Mass resolution (M/ΔM) was above 5,000 and mass calibration of silver was achieved using silver foil (Goodfellow, Huntingdon, UK). The isotopic ratio of ¹⁰⁷Ag⁻ to ¹⁰⁹Ag⁻ was checked to ensure that no isobaric interference had caused inaccurate measurements.
For each sample, at least 30 bacterial cells were selected on different images and analyzed individually. Data were treated with the module ImageJ plugin to process and analyze the captured images [23]. After selection of regions of interest (ROI), which here corresponded to a single bacterium, the total number of counts of each element detected by SIMS was calculated pixel by pixel. Digital masks of bacterial cells were defined using the ¹²C¹⁴N⁻ image, which provided both the highest contrast and widest ranges of grey values. Phosphorus and sulphur signals were indicative of cell components [22, 23]. The ¹⁰⁷Ag⁻ signal enabled the detection of silver atoms present in the cells. To obtain comparable Ag intensity values independently from variations in the intensity of the NanoSIMS50 signal occurring during measurements, the ratio of ¹⁰⁷Ag⁻ to ¹²C¹⁴N⁻ intensity counts (Ag/CN, dimensionless) was calculated for each selected bacterium, based on the assumption of a relatively constant ¹²C¹⁴N⁻ content in the selected area (i.e., the surface of the cell section) [18]. Although the real Ag concentration within cells cannot be assessed, this method made it possible to compare silver uptake between the different samples under study. The mean ± SD of these measurements was then determined, resulting in a “mean Ag/CN intensity ratio” for each sample.

## Microarray procedure and gene expression data analysis

At 3 h of culture, 60 mL of culture (at OD₅₀₅ nm = 0.73 ± 0.07, 0.63 ± 0.02, 0.12 ± 0.01 and 0.11 ± 0.02 for cultures subjected to 0, 5.0, 6.5 and 8.5 µM AgNO₃, respectively) were sampled, placed in 50-mL BD Falcon™ tubes and centrifuged at 4,000 rpm for 5 min at 4°C, after which the pellet was frozen in liquid nitrogen. The pellet was resuspended in 1 mL of TE buffer (Tris-HCl 10 mM, pH 8, EDTA 1 mM, 1 mg lysozyme) and incubated for 5 min at room temperature. Total RNA was extracted with an RNeasy midi kit (Qiagen), including the DNase I treatment described in the manufacturer's instructions. The total quantity and integrity of RNA were checked with a Nanodrop® and an Agilent BioAnalyzer, respectively. A double strand cDNA synthesis kit (InvitroGen) was used to produce cDNA from 1.5 µg of total RNA. Following the manufacturer's instructions, 1 µg of cDNA was labeled using the one color DNA labeling kit and 2 µg of labeled cDNA were hybridized onto E. coli K-12 gene expression arrays (Nimblegen, Roche) for 17 h at 42°C. Arrays were later washed and scanned using a MS200 Microarray Scanner (Nimblegen, Roche). The pictures were analyzed with DEVA 1.2.1 software. All array procedures were performed by the GeT-Biopuces platform (http://get.genotoul.fr). Raw and processed data were deposited in the Gene Expression Omnibus data repository (GSE67735).

Raw probe intensities (three replicates for each concentration of ionic silver) were analyzed in an R computing environment [24] using the affy [25] and limma [26] package of BioConductor. Raw data were subjected to an RMA-base background correction [27]. After background correction, intra-replicate quantile normalization was performed for each concentration of ionic silver. A set of probes in whose background the ranks were roughly the same across all the 12 arrays was selected. The median value of the intensities of the probe set for each array was used as a scaling factor for normalisation between growth conditions. After normalisation, the expression of a gene was calculated by RMA summarisation [27] for each concentration of ionic silver. Differential gene expression was evaluated using a modified t-test combined with an empirical Bayes method [28]. The p-values for the contrast of interest were adjusted for multiple testing by the Benjamini-Hochberg False Discovery Rate (FDR) [29]. A p-value threshold of 1% and a fold change < 0.5 or > 2 were used for differential expression. Agglomerative hierarchical clustering of the biological samples on the extracted genes for heat map display was derived using the Euclidean distance as the similarity measure and the Ward agglomeration method using the stats R package [24]. Gene ontology enrichment analysis of the extracted genes was performed using the topGO R package [30].
## Data integration

To identify and interpret correlation patterns in several large data matrices from the same biological system under investigation, sparse partial least squares regression (sPLSR) [31] was applied using the mixOmics R package [32]. To this end, three data sets were used:

1. Variations in internal and external silver contents and physiological parameters from the designed study (in replicates) corresponding to Table 1.

2. The gene expression data of size (12x3486) containing expression values (in triplicate for each concentration of ionic silver) for the 3,486 genes differentially expressed in at least one of the three concentration of ionic silver under ionic-silver mediated stress. Two gene expression subsets were also used, restricted to genes either involved in fatty acid metabolism (size 12x106) or protein folding (size 12x104).

3. The sFTIR data of size (75x271) were provided by the second derivative spectra from the sFTIR microscopy data, as indicated above and detailed in [15]. The three wavenumber regions of interest were selected for integration: (i) 3,100–2,700 cm⁻¹ (fatty acids), (ii) 1,800–1,480 cm⁻¹ (representative characteristic of C = O stretching vibration, amide I and amide II bands of peptide bonds) and (iii) 1,450–1,150 cm⁻¹ (PO₂⁻ groups and, for a part, C-O stretching modes).

In sPLSR, the two-block data matrices to be integrated were denoted X(n x p) and Y(n x q), where p and q are the total number of variables measured on the same n subjects. When two data sets to be integrated displayed different sample numbers n, the dataset with the lowest sample number was completed by a bootstrap method to reach the higher n value. We analysed the data sets using sPLSR in an attempt to explain/predict the Y variables with respect to X and to perform simultaneous variable selection in the two datasets [31]. The optimum number of PLS components in the sPLSR models and the number of variables to be selected on each dimension were found by cross-validation by calculating the Q2 criterion and the minimum root mean square error (RMSE), respectively. Correlation circle plots and relevance networks resulting from the sPLSR approach [33] were used to help interpret the results. To highlight the strongest variable associations only, variables with a correlation greater than 0.5 in absolute value were chosen to infer the relevance networks (see figure legends). This threshold was arbitrarily chosen to obtain biologically interpretable networks that were neither too sparse nor too dense. The networks obtained were then used as inputs for Cytoscape [34] for visualisation.

To correlate the intracellular silver content and physiological parameters with biochemical changes within bacterial cells, sPLSR analysis was performed with variations in AgNO₃ concentrations and physiological parameters data as X and the selected sFTIR bands as Y. When sPLSR was used to integrate the gene expression data and the sFTIR spectra, the sPLSR was run with the gene expression as X and the sFTIR spectra as Y. Variations in AgNO₃ concentrations and physiological parameters and gene expression were integrated by running sPLSR using the AgNO₃ concentrations and physiological parameters data as X and the gene expression as Y.

## Results

### Physiological and biochemical responses of *E. coli* to ionic silver

Physiological parameters of *E. coli* cells grown in the presence of different concentrations of ionic silver ([AgNO₃] = 5.0, 6.5 or 8.5 µM in the growth medium from inoculation) were assessed and compared to those of the control experiment (i.e., bacteria grown without any
<table><caption>**Table 1.** Physiological parameters of *E. coli* cultures (in M9 medium supplemented with 2.5 g/L of glucose) as a function of the concentration of extracellular AgNO₃. The value μt = 3h corresponds to the specific growth rate calculated between 2.5 and 3.5 h of culture. The μexpo value is the specific growth rate over the entire exponential growth phase (start and end specified in brackets) after the beginning of growth. The mean ± SD values of three independent experiments are given, except at a concentration of 6.5 μM AgNO₃ for which the values of all three experiments are detailed.</caption><thead><tr><th>[AgNO<sub>3</sub>] in the medium (μM)</th><th>Lag time (h)</th><th>μ<sub>t = 3h</sub> (h<sup>-1</sup>)</th><th>Number of cultivable cells at 3 h (CFU/mL/ODU)</th><th>μ<sub>expo</sub> (h<sup>-1</sup>)</th><th>Maximum absorbance (at 595 nm)</th></tr></thead><tbody><tr><td>0 (control)</td><td>0.5 ± 0.0</td><td>0.62 ± 0.01 (from 2.5 to 3.5 h)</td><td>5.2 ± 1.3 x 10<sup>8</sup></td><td>0.63 ± 0.01 (from 0.5 to 5.5 h)</td><td>3.31 ± 0.12</td></tr><tr><td>5.0</td><td>0.5 ± 0.0</td><td>0.63 ± 0.03 (from 2.5 to 3.5 h)</td><td>4.3 ± 1.2 x 10<sup>8</sup></td><td>0.59 ± 0.03 (from 0.5 to 5.5 h)</td><td>3.19 ± 0.11</td></tr><tr><td>6.5</td><td>9</td><td>0</td><td>2.4 x 10<sup>8</sup></td><td>0.54 (from 10.0 to 14.5 h)</td><td>3.13</td></tr><tr><td>6.5</td><td>17</td><td>0</td><td>2.3 x 10<sup>8</sup></td><td>0.44 (from 18.0 to 24.5 h)</td><td>3.36</td></tr><tr><td>6.5</td><td>> 24</td><td>0</td><td>4.3 x 10<sup>7</sup></td><td>0</td><td>0.11</td></tr><tr><td>8.5</td><td>> 24</td><td>0.0 ± 0.0</td><td>1.1 ± 0.2 x 10<sup>7</sup></td><td>0.0 ± 0.0</td><td>0.10 ± 0.01</td></tr><tr><td>100<sup>a</sup></td><td>—</td><td>0.0 ± 0.0</td><td>0.0 ± 0.0</td><td>0.0 ± 0.0</td><td>—</td></tr></tbody></table>

doi:10.1371/journal.pone.0145748.t001

added silver) (Table 1). A concentration of 5.0 μM of AgNO₃ had a weak effect on *E. coli* growth, as the cells reached a μexpo of 0.59 ± 0.03 h⁻¹ (0.63 ± 0.01 h⁻¹ for the control). At a concentration of 6.5 μM of AgNO₃, the three replicates exhibited quite different physiological behaviours. In one replicate, no growth was observed before the end of the experiment (i.e. 24 h), although bacterial cells were still able to grow on LB agar medium after washing (Table 1). In the two other replicates, exponential growth occurred at a lower μexpo than the control (0.54 and 0.44 h⁻¹) after a lag phase of 9 h and 17 h, respectively. At the end of these two cultures, the maximum concentration of biomass reached a similar value to that obtained at 5.0 and 0 μM. At a concentration of 8.5 μM of AgNO₃, no bacterial growth was observed over a period of 24 h. Although these cells did not grow in the presence of ionic silver, after washing (i.e. removal of the silver), a significant number of cells was still cultivable (10⁷ CFU/mL/ODU). For the purpose of comparison, at a concentration of 100 μM of AgNO₃, 100% of the cells were no longer cultivable, showing that 100 μM of AgNO₃ was a lethal concentration (Table 1). Taken together, these results show that 6.5 μM is a pivotal value for the physiological effect of ionic silver, and that concentrations of up to 8.5 μM of AgNO₃ can be considered as sub-lethal.

To visualise the presence of intracellular silver after 3 h of culture, the silver element was mapped by NanoSIMS on at least 30 individual cells in each sample (Fig 1). The distribution of the ¹²C¹⁴N⁻ cluster was simultaneously recorded as a “fingerprint” of bacterial cells, while the signals of phosphorus and sulphur made it possible to highlight interactions between the ionic silver and cell components (Fig 1). As shown in the images of bacterial cells grown in the presence of AgNO₃, silver mainly overlapped sulphur. The distribution of phosphorus was quite homogeneous, which was consistent with its ubiquitous presence in the cell (nucleic acids, phosphorylated metabolites, etc.).

To obtain a semi-quantitative value, the intensity ratio between the two signals (Ag/CN) was then calculated for each cell [35] and averaged for each group (Fig 2). A significant increase in the intracellular silver content was observed in cells treated with 5.0 μM AgNO₃, as the Ag/CN ratio was multiplied by a factor of 15 compared to untreated cells (Fig 2). The Ag/CN ratio measured in bacteria treated with 6.5 μM of AgNO₃ increased by a factor 1.5 compared to the cells treated with 5.0 μM of AgNO₃, while it was of the same order of magnitude as the latter in cells treated with 8.5-μM of AgNO₃ (Fig 2).

<footnote><sup>a</sup>: values from [15].</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Distribution maps of CN (in grey), P (in blue), S (in green) and Ag (in red) elements detected by NanoSIMS analysis in *E. coli* cells after 3 h of culture without ionic silver (control) or in the presence of 5.0, 6.5 and 8.5 µM of AgNO₃. The fourth column corresponds to overlapping of the three-coloured images (Composite Image mode [23]). Scale bar: 2 µm.</figcaption></figure>

doi:10.1371/journal.pone.0145748.g001

Integrative analysis of the physiological parameters and the biochemical
response of *E. coli* after Ag⁺ exposure

The first approach used to integrate the data associated variations in the concentration of extra-
cellular AgNO3, intracellular silver content (reflected by the Ag/CN ratio), physiological
parameters (specific growth rate µt = 3h and the number of cultivable cells expressed in CFU/
mL/ODU) and modifications in E. coli biochemical composition, evaluated by *in situ* sFTIR
micro-spectroscopy at the single-cell scale (data and peak assignments from [15]; see the exper-
imental section for further details). Raw sFTIR spectra recorded on bacteria grown in the
absence or in the presence of ionic silver at sub-lethal concentrations, are presented in S1 Fig.
The peaks were identified based on our previous study [15].

As expected, the intracellular silver content was correlated with the concentration of extra-
cellular AgNO3, and both were anti-correlated with the specific growth rate µt = 3h and with the
number of cultivable cells (Fig 3A). Furthermore, the intracellular silver content was correlated
with wavenumber characteristics of fatty acids and peptide bonds (Fig 3A and 3B). In particu-
lar, the Ag/CN ratio was positively correlated with the wavenumbers 2,873, 2,869 and 2,865
cm⁻¹, assigned to the symmetric C-H stretching of—CH3 groups in lipids, but negatively corre-
lated with the wavenumbers 2,850 and 2,846 cm⁻¹, representative of the C-H symmetric stretch
of >CH2 groups. This showed that the presence of silver in E. coli cells led to an increase in
CH3 over CH2 groups, corresponding to a shortening of the aliphatic chains of bacterial lipids.
Positive and negative correlations of the Ag/CN ratio were also observed with respectively, the
up-shift (2,939 cm⁻¹) and the downshift (2,915 cm⁻¹) of the asymmetric C-H stretching of
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Variation in the Ag/CN ratio with the extracellular AgNO₃ concentration deduced from NanoSIMS analysis. For each condition, the number of bacterial cells used for calculation (n) is indicated.</figcaption></figure>

doi:10.1371/journal.pone.0145748.g002

>CH₂ groups (2,930 cm⁻¹), suggesting an increased proportion of “cis” over “trans” conformation in unsaturated fatty acids of cells subjected to ionic silver-mediated stress.

In the spectral range from 1,800 to 1,700 cm⁻¹, it will be noted that the wavenumbers 1,743 and 1,739 cm⁻¹, which are characteristic of C = O stretching vibrations in lipid esters [36], were positively correlated with the specific growth rate and the number of cultivable cells but negatively correlated with the intracellular silver content and the extracellular concentration of AgNO₃ (Fig 3B). Concerning the peptide bond, the intracellular silver content and the extracellular concentration of AgNO₃ were correlated with wavenumber 1,670 cm⁻¹, the AgNO₃ concentration with 1,681 cm⁻¹ while 1,697, 1,677 and 1,673 cm⁻¹ were anti-correlated with the
<figure><img src="image_6.png" /><figcaption>Fig 3. Results of the sPLSR integrative approach between the Y matrix, corresponding to AgNO₃ concentrations and physiological parameters (i.e., specific growth rate μt = 3h, number of cultivable cells (CFU/mL/ODU) and intracellular silver content illustrated by the Ag/CN ratio at 3 h culture) and the X matrix composed of the selected sFTIR data (i.e. the regions 3,100–2,800 cm⁻¹, 1,700–1,600 cm⁻¹, 1,600–1,480 cm⁻¹ and 1,300–1,150 cm⁻¹, representative of fatty acids, C = O stretching vibration, amide I and amide II bands and PO2- groups and C-O stretching mode, respectively). The regions 2,800–2,700 cm⁻¹, 1,800–1,700 cm⁻¹ and 1,480–1,300 cm⁻¹ are grouped under “Others”. (A) Correlation circle plots for</figcaption></figure>
dimensions 1 and 2. The subsets of correlated wavenumbers are represented by thick points coloured according to the legend, while the physiological parameters are represented by their name in black. (B) Relevance networks resulting from the sPLSR approach. A threshold of 0.5 was used for relevant correlations. Green and red edges indicate positive and negative correlations, respectively. Physiological parameters and wavenumbers are represented as circles and ellipses, respectively.

doi:10.1371/journal.pone.0145748.g003

number of cultivable cells (Fig 3B). These wavenumbers, assigned to anti-parallel pleated β-sheets of amide I, unordered random coil and turn forms of amide I [36], suggest a disordered secondary structure of proteins after Ag⁺ exposure. Some nucleic acid damage is possible as wavenumber 1,592 cm⁻¹ [36] was positively correlated with the Ag/CN ratio and with the extracellular AgNO₃ concentration (Fig 3B). The extracellular concentration of AgNO₃ and intracellular silver content were positively correlated with wavenumbers 2,750 and 2,753 cm⁻¹ and negatively correlated with 2,734 cm⁻¹, while the peak at 2,780 cm⁻¹ was positively correlated with the Ag/CN ratio and negatively correlated with the specific growth rate and the number of cultivable cells (Fig 3B). These four wavenumbers are characteristic of N-H stretching [36]. The wavenumbers 1,400 and 1,396 cm⁻¹, attributed to symmetric CH₃ bending of the methyl groups of proteins, and the wavenumber 1,164 cm⁻¹ attributed to C-O stretching in serine, threonine and tyrosine of proteins [36], were positively correlated with the specific growth rate and the number of cultivable cells but negatively correlated with the Ag/CN ratio and the extracellular concentration of AgNO₃ (Fig 3B). Finally, wavenumber 1,781 cm⁻¹ could not be attributed to any biochemical component of the cell.

## Transcriptional response of *E. coli* to ionic silver

To analyse the impact of ionic silver on global gene expression, transcriptomic analyses were performed after a 3-h exposure to concentrations of 5.0, 6.5 and 8.5 µM of ionic silver, and compared to control conditions, i.e. without the addition of any silver. At a concentration of 5.0 µM of silver, the yield of total extracted RNA was not affected whereas at concentrations of 6.5 and 8.5 µM, 2-fold less yield was measured. In addition, at these two AgNO₃ concentrations, the proportion of mRNA to total RNA (estimated by summing all the spot intensities in the array) was reduced more than 3-fold. The reduced proportion of mRNA may be linked to a decrease in transcription and/or an increase in mRNA degradation at the highest concentrations of ionic silver.

Differentially expressed genes were determined by comparing expression under the three above mentioned conditions. A total of 3,486 genes (corresponding to 81% of the total of 4,254 *E. coli* genes) were significantly differentially expressed (p-value < 1% with a fold change < 0.5 or > 2) at least in one of the three comparisons (S1 Table). A comparative transcriptomic analysis was first performed to determine the extent of overlap between our transcriptomic response to ionic silver (this study) and transcriptomes of *E. coli* cells responding to other stressful conditions depicted in the literature (i.e., subjected to from sub-lethal to lethal concentrations of other metal ions, under oxidative conditions or after a heat shock) (Table 2). High correlations were observed between 6.5 µM AgNO₃, 8.5 µM AgNO₃, 6 mM H₂O₂ and heat shock at 45°C (correlation coefficients > 0.35). Under all these conditions, the imposed stress led to growth arrest. No strong similarity was observed between our silver ion response and the other metal ion responses reported in the literature except between 8.5 µM AgNO₃ and 2 mM copper ion (correlation coefficient of 0.15) (Table 2). It is interesting to note that the 5.0 µM AgNO₃ treatment did not correlate with the 8.5 µM treatment and only slightly with the 6.5 µM treatment. These results suggest that the transcriptomic response really only took place at a concentration of more than 5 µM.
<table><caption>**Table 2. Pearson correlation between transcriptomic matrix for a selection of stress and ion response conditions.** The similarity between the log2 matrix of the gene expression ratio vs. the untreated control was determined under the following conditions: oxidative stress (6 mM H₂O₂, time 10 minutes from [53]), heat shock (45°C, time 10 minutes from [54]), silver ion response (this study, in bold), copper ion response (2 mM, from [55]), cadmium ion response (0.1 mM, from [56]), zinc ion response (0.2 mM, from [48]) and cobalt ion response (0.25 mM, from [57]). Correlations coefficients ≥ 0.35 are in bold.</caption><thead><tr><th></th><th>H₂O₂, 6 mM</th><th>HS, 45°C</th><th>Ag, 8.5 µM</th><th>Ag, 6.5 µM</th><th>Cu, 2 mM</th><th>Cd, 0.1 mM</th><th>Zn, 0.2 mM</th><th>Ag, 5 µM</th><th>Co, .25 mM</th></tr></thead><tbody><tr><td>H₂O₂, 6 mM</td><td><strong>1.00</strong></td><td><strong>0.54</strong></td><td><strong>0.47</strong></td><td><strong>0.35</strong></td><td>0.19</td><td>0.04</td><td>-0.01</td><td>-0.16</td><td>-0.05</td></tr><tr><td>HS, 45°C</td><td><strong>0.54</strong></td><td><strong>1.00</strong></td><td><strong>0.52</strong></td><td><strong>0.45</strong></td><td>0.24</td><td>0.10</td><td>-0.05</td><td>-0.14</td><td>-0.02</td></tr><tr><td>Ag, 8.5 µM</td><td><strong>0.47</strong></td><td><strong>0.52</strong></td><td><strong>1.00</strong></td><td><strong>0.78</strong></td><td>0.15</td><td>0.09</td><td>0.07</td><td>0.02</td><td>-0.02</td></tr><tr><td>Ag, 6.5 µM</td><td><strong>0.35</strong></td><td><strong>0.45</strong></td><td><strong>0.78</strong></td><td><strong>1.00</strong></td><td>0.25</td><td>0.23</td><td>0.09</td><td>0.14</td><td>0.04</td></tr><tr><td>Cu, 2 mM</td><td>0.19</td><td>0.24</td><td>0.15</td><td>0.25</td><td><strong>1.00</strong></td><td>0.11</td><td>0.13</td><td>-0.05</td><td>0.07</td></tr><tr><td>Cd, 0.1 mM</td><td>0.04</td><td>0.10</td><td>0.09</td><td>0.23</td><td>0.11</td><td><strong>1.00</strong></td><td>-0.05</td><td>0.06</td><td>0.05</td></tr><tr><td>Zn, 0.2 mM</td><td>-0.01</td><td>-0.05</td><td>0.07</td><td>0.09</td><td>0.13</td><td>-0.05</td><td><strong>1.00</strong></td><td>0.09</td><td>0.02</td></tr><tr><td>Ag, 5 µM</td><td>-0.16</td><td>-0.14</td><td>0.02</td><td>0.14</td><td>-0.05</td><td>0.06</td><td>0.09</td><td><strong>1.00</strong></td><td>0.06</td></tr><tr><td>Co, 0.25 mM</td><td>-0.05</td><td>-0.02</td><td>-0.02</td><td>0.04</td><td>0.07</td><td>0.05</td><td>0.02</td><td>0.06</td><td><strong>1.00</strong></td></tr></tbody></table>

doi:10.1371/journal.pone.0145748.t002

To get further insight into the progressive regulation of gene expression with an increase in the concentration of AgNO₃, a clustering procedure was implemented to distinguish different gene expression regulation profiles (Fig 4A). Clustering of gene expression ratios between conditions (5.0 vs. 0, 6.5 vs. 0 and 8.5 vs. 0 µM AgNO₃) differentiated six clusters (Fig 4B, S1 Table). All clusters exhibited a weak transcriptomic response in the first comparison “5.0 vs. 0 µM AgNO₃” confirming the lack of strong transcriptomic response at 5.0 µM. Above 5.0 µM, the profiles differed in the six clusters. The first five clusters (1 to 5) generally corresponded to down-regulated gene expression over the whole range of AgNO₃ concentrations but differed in the magnitude of the decrease between 8.5 µM AgNO₃ and 0 µM: clusters 1 and 2 corresponded to moderate decreases (average 15-fold, 53% of differentially expressed genes), cluster 5 included small decreases (average 4-fold, 25% of differentially expressed genes) while the very large decreases (average 200-fold) corresponding to 17% of the differentially expressed genes were grouped in clusters 3 and 4. Cluster 6 (136 genes, 4% of differentially expressed genes) contained genes up-regulated at 6.5 µM and/or 8.5 µM of AgNO₃. Therefore, except for 4% of the differentially expressed genes, gene expression was down-regulated in the presence of ionic silver at a concentration > 5.0 µM. This is in agreement with the smaller proportion of mRNA of total RNA at high AgNO₃ concentrations.

We then examined the function of genes whose regulation differed from the general trend of moderate down-regulation. Enrichment in functional category was thus determined (p-value < 0.05) in clusters 3 and 4 of strongly down-regulated expression, in cluster 5 of less than the average down-regulated expression, and in cluster 6 of up-regulated gene expression (Table 3). Pathways related to growth, i.e. metabolic processes, transports, membrane biogenesis, cell shape regulation and the translational process, were over-represented in the strongly down-regulated subgroup of genes.

In addition to the general down-regulation of growth-associated functions, specific responses to ionic-silver mediated stress were also observed. Cell motility, chemotaxis and flagellum organization were strongly affected. The most striking feature was the response to metal ions: genes involved in response to copper (*copA*, a copper efflux transport and *cueO*, a copper oxidase) and cadmium (*zinT*, a zinc and cadmium binding protein and *ygiW*, unknown function) were enriched in the group of up-regulated genes, whereas the Cu(I) efflux pump (*cusCFBA*) and the zinc transporter (*zupT*) were less down-regulated than average. We found up-regulated expression of transport of magnesium (*mgtA*), manganese (*mntH*) and calcium (*chaA* and *nhaA*). A general response to stress was observed with the enrichment of categories
A

<figure><img src="image_7.png" /></figure>

B

<figure><img src="image_8.png" /><figcaption>**Fig 4.** Ionic silver-mediated expression regulation. (A) Hierarchical clustering of expression ratios between growth conditions with 5.0, 6.5 or 8.5 µM of AgNO₃ in the medium and the control (i.e., no ionic silver added) for the 3,486 differentially-expressed genes in one of the three ionic-silver mediated stress conditions. The resulting heat map shows the samples in rows and the genes in columns, red indicates up-regulation, green down-regulation, and black no change. (B) The expression of genes in each cluster is plotted as a function of the concentration of AgNO₃ in the medium to show the general trend of expression regulation. The red line in the individual clusters represents the median calculated from the expression levels of the genes. Cluster 1: 1,131 genes, cluster 2: 746 genes, cluster 3: 516 genes, cluster 4: 76 genes, cluster 5: 881 genes and cluster 6: 136 genes.</figcaption></figure>

doi:10.1371/journal.pone.0145748.g004
<table><caption>**Table 3. Functional analysis of the genetic response of E. coli cells to ionic silver.** Selected functional categories enriched (p-value < 0.05) in groups of genes with generally strongly down-regulated expression (clusters 3 and 4), weakly down-regulated expression (cluster 5) and up-regulated expression (cluster 6) are included. For each category, associated genes are also listed.</caption><thead><tr><th>Functional category</th><th>Gene</th></tr></thead><tbody><tr><td colspan="2"><strong>Strong down-regulation (clusters 3 and 4)</strong></td></tr><tr><td colspan="2"><strong>Metabolic process</strong></td></tr><tr><td>Glycolysis</td><td>sucA, lpd, eno, fbaAB, gpmAM, pykF, pgi</td></tr><tr><td>Glyoxylate cycle</td><td>aceA,B,K, iclR</td></tr><tr><td>Galactitol</td><td>gatA,B,C,D</td></tr><tr><td>Amino acids</td><td>folD, ilvH,I,D,B, hisA,B,C,D,F,G,H,I, trpC,D,E, aroA,C,G,K, pheA</td></tr><tr><td>Bases</td><td>adk, fold, guaA, gmk, purB,C,D,F,K,H,L,M,N,T,U, carA,B, pyrCDGH, relA, spoT, thyA</td></tr><tr><td>Glycogen energy reserve</td><td>glgB,C,P,X</td></tr><tr><td colspan="2"><strong>Transport</strong></td></tr><tr><td>Maltose</td><td>malE,F,K</td></tr><tr><td>Proteins and peptides</td><td>dppA,B,C,D,F, oppA,B,C,D,F, secA,B, lolB,C,E, tatA,B</td></tr><tr><td>Iron</td><td>fecIR, fep, fes, fhuC,E, feoAB, tolC, tonB, entACDE, fiu, efeO, fre</td></tr><tr><td>Lipopolysaccharides</td><td>lptA,B,C,D</td></tr><tr><td>Electrons</td><td>frdA, napA, narG, nirB, nuoA,B,C,E,F,G,H,I,J,K,L,M,N, cyoA,B, C, sucB,C,D, sdhA,C,D, aceA,B,K, lpd, dld, ppc</td></tr><tr><td>Protons</td><td>atpB,D,F,H,I</td></tr><tr><td colspan="2"><strong>Membrane biogenesis</strong></td></tr><tr><td>Lipopolysaccharides</td><td>rfbAC,D,X, waaF,G,Q, wbbI,J,K, rffEG, rfaD, rmlB, rfc, rfe, glf, cld, lpxD,H,P, kdsD, glmU, gmhB, yijP</td></tr><tr><td>Cell outer membrane</td><td>bamA,D, lptD, hlpA, surA</td></tr><tr><td><strong>Cell shape regulation</strong></td><td>murB,C,E, mreB,C,D, rodZ, dacA, ftsW, mpl, mraY, yafK, glmU, ispU</td></tr><tr><td colspan="2"><strong>Translation</strong></td></tr><tr><td>Ribosome generation</td><td>rsuA, rluB,D,E</td></tr><tr><td>tRNA synthetase</td><td>aspS, gltX, glyS,Q, hisS, lysS, metG, pheS, trpS</td></tr><tr><td>Translation termination</td><td>frr, prfA,B,C</td></tr><tr><td><strong>Cell motility, chemotaxis and flagellum organization</strong></td><td>cheA,B,R,W,Y,Z, dppA, flgA,B,C,D,E,F,G,H,I,J,K,L,M,N, flhA,C, fliC,D,F,G,H,I,J,K,L,M,N,O,S,T, frdA, motA,B, rbsB, tap, tar, tsr, yhjH</td></tr><tr><td colspan="2"><strong>Small down-regulation (cluster 5)</strong></td></tr><tr><td colspan="2"><strong>Transport</strong></td></tr><tr><td>Copper ion transmembrane</td><td>zupT, cusA,B,C,F</td></tr><tr><td>Organic anion</td><td>abgT, glcA, phnC,J,L,O,P, gltS, ychM</td></tr><tr><td>Protein secretion (type II system)</td><td>hofC, yghF, gspC,E,F,G,J</td></tr><tr><td>Gluconate</td><td>ygbN, dsdX, gntP, idnT, yjhF</td></tr><tr><td colspan="2"><strong>Metabolic process</strong></td></tr><tr><td>Aromatic compound</td><td>paaF,G,H,J, hcaC,D,F, mhpC,D,E,F, tyrR, feaB, hcaB, tnaA</td></tr><tr><td>Organic acids</td><td>dgoA,D,K, yjjM, cynT,S,X</td></tr><tr><td>Arabinose</td><td>araA,C, sgbE, fucl,K</td></tr><tr><td>Purines</td><td>yagR,S, allA,C, xdhA,B,C,D</td></tr><tr><td>Arginine</td><td>speA, astB,C,D,E</td></tr><tr><td colspan="2"><strong>Up-regulation (cluster 6)</strong></td></tr><tr><td colspan="2"><strong>Response to stimulus</strong></td></tr></tbody></table>

(Continued)
<table><caption>**Table 3. (Continued)**</caption><thead><tr><th>Functional category</th><th>Gene</th></tr></thead><tbody><tr><td>Stress</td><td>clpB, dnaJ,K, htpG,X, lon, pspA,B,C,D,G, raiA, psiE, dps, ibpA, B, ydel, rmf, bhsA, ariR, yodD, zinT, ygaW, ygiW, yhbO, ytfE, tisB, mntH</td></tr><tr><td>Heat</td><td>dnaJ,K, htpG, lon, pspA, ibpA,B, yhbO</td></tr><tr><td>Hydrogen peroxide</td><td>ydel, ariR, yodD, ygiW</td></tr><tr><td>Copper ion</td><td>cueO, copA</td></tr><tr><td>Cadmium ion</td><td>zintT, ygiW</td></tr><tr><td>Acidity</td><td>asr, yodD</td></tr><tr><td><strong>Transport</strong></td><td></td></tr><tr><td>Divalent metal ion</td><td>chaA, mgtA, mntH</td></tr><tr><td>Inorganic anion</td><td>yfeO, narU, sbp</td></tr><tr><td><strong>Biofilm formation</strong></td><td>ydeH,I, bhsA, ariR, yodD, msqR, ygiW</td></tr><tr><td><strong>Homeostasis</strong></td><td></td></tr><tr><td>Cellular sodium ion</td><td>chaA, nhaA</td></tr><tr><td><strong>Biosynthetic process</strong></td><td></td></tr><tr><td>Methionine</td><td>metA,B,F,R, mmuM</td></tr><tr><td>Arginine</td><td>argA,B,C</td></tr></tbody></table>

doi:10.1371/journal.pone.0145748.t003

associated with stress in general and in particular heat shock proteins in the group of up-regulated genes. At the metabolic level, higher expression of the methionine- and arginine-related pathways was observed. Other pathways were also down-regulated by ionic silver, i.e. the catabolism of aromatic compounds, organic acids and sugars and the metabolism of purine. Altogether, these results revealed a marked impact of ionic silver on the gene expression profile, reducing general growth-related functions and favouring stress responses.

## Integrative analysis of the physiological parameters and the transcriptional response of *E. coli* after Ag⁺ exposure

To investigate the physiological significance of the Ag⁺-mediated modifications in gene expression, we analysed correlation patterns between the gene expression profiles, the extracellular silver concentration in the medium (AgNO₃ concentration) and in the bacterial cells (Ag/CN ratio), and the physiological parameters (specific growth rate µt = 3h and the number of cultivable cells after 3 h of culture), using sparse partial least squares regression (sPLSR) analysis [31] (Fig 5). Among the correlated variables, most gene expression regulations were positively correlated with the specific growth rate and the number of cultivable cells and, as expected, negatively correlated with the external and intracellular silver contents. Expression was mostly down-regulated, and concerned (i) many transporters [32 genes including transport of copper (cusAB), nickel (nikB,C,D,E)], (ii) 10 transcriptional regulators (dctR, dsdC, feaR, hcaR, metJ, mhpR, sgrR, yfeR, yfhH and yqhC), (iii) genes involved in oxidation-reduction processes [27 genes including nitrite/nitrate reductase (napD, narI,J, nirB and nrfA), fumarate reductase (frdB,C,D), dimethyl sulfoxide reductase (dmsB,C), xanthine oxidoreductase (xdhB,D), hydrogenase (hypA,E) and cupric reductase (ndh)], (iv) genes coding for membrane/periplasmic proteins (nikA, ompW, yihG,H, yjjB, yohJ and zraP) and ribosomal proteins (rplB,J and rpsC,J,N). The AgNO₃ concentration and the Ag/CN ratio were rarely positively correlated with the regulation of gene expression. All the regulations positively correlated with the AgNO₃ concentration and the Ag/CN ratio were up-regulated in the presence of ionic silver. The genes were shown to be involved in responses to general stress (dnaK and pspA) and copper (copA and
<figure><img src="image_9.png" /><figcaption>**Fig 5.** Relevance networks resulting from the sPLSR approach between the Y matrix corresponding to AgNO₃ concentrations and physiological parameters (i.e., specific growth rate µt = 3h, number of cultivable cells (CFU/mL/ODU) and intracellular silver content illustrated by the Ag/CN ratio at 3 h culture) and the X matrix representing the expressions of the 3,486 differentially expressed genes in one of the three AgNO₃-mediated stress growth conditions. A threshold of 0.5 was used for relevant correlations. Green and red edges indicate positive and negative correlations, respectively. Physiological parameters are represented as circles and genes as rectangles.</figcaption></figure>

doi:10.1371/journal.pone.0145748.g005

*cueO*, ion transport (magnesium *mgtA*, sodium *nhaR*), and methionine biosynthesis (*metA*, *R*) and putatively in iron homeostasis (*yqji*) and in adhesion (*ydeS*, *R*). The combination of the different regulation of gene expression constitutes the genetic pattern of the physiological response of *E. coli* to Ag⁺-mediated stress.

## Integrative analysis of the biochemical and transcriptional responses of *E. coli* after Ag⁺ exposure

As shown above, ionic silver stress mainly affected fatty acids (chain length, saturation and conformation) and peptide bonds (secondary structure). We next focused on regulation of the expression of genes specifically involved in fatty acid metabolism and protein folding in the presence of Ag⁺ ions to provide additional insights into Ag⁺-related functional genomics.

Regarding fatty acid metabolism, we observed down-regulation of genes involved not only in the initiation of fatty acid biosynthesis (*accA*, *B*, *C*, *D*, *fabD* and *fabH*) and activation of fatty acid synthase (*acp H*, *S*, *T*) but also in fatty acid catabolism (aerobic β-oxidation via *fadA*, *B*, *D*) and transport (*fadL*, a long-chain fatty acid outer membrane transporter) (S1 Table). Interestingly, fatty acid elongating genes (*fabG*, *fabZ*, *fabI*, *fabF*, *fabB* but not *fabA*) were also down-regulated. FabB is the first enzyme in the unsaturated fatty acid elongation pathway while FabF initiates elongation of the unsaturated fatty acid palmitoleate in cis-vaccenate. This general down-regulation of fatty acid metabolism could be directly linked to the general shortening of the acyl chain, shown in the presence of ionic silver. Moreover, although FabA and FabZ are isoform enzymes, only FabA can cause unsaturation in fatty acids. The relative level of FabA compared to that of FabB and FabZ controls the ratio of saturated to unsaturated fatty acids in
E. coli [37, 38]. Fewer unsaturated bonds would be obtained at increasing concentrations of ionic silver if (i) FabA activity was lower than FabZ and/or if (ii) FabB activity was low. At the expression level, we found no such regulatory profiles, *fabA* was expressed at similar or higher levels than *fabZ* at 6.5 and 8.5 µM AgNO₃ and *fabB* expression, despite its down-regulation, was not significantly lower than *fabA* (S2 Fig).

A search for correlations between selected sFTIR data and gene expression regulation revealed negative correlations between the fatty acid wavenumbers 2,873, 2,869 and/or 2,865 cm⁻¹ and the expression of eight genes (Fig 6): *acpH* (of the *acpTSH* cluster involved in activation of fatty acid synthase), *acnB* (aconitase), *ydiF* (predicted acetyl-CoA transferase), *paaF* (of the phenylacetate metabolism) and *ychM* (predicted inner membrane transporter) and more interestingly, *fadL*, a long-chain fatty acid outer membrane transporter, *lpxA* and *arnA*, catalyzing biosynthesis and phosphate modification of lipid A, respectively. These results show that a higher peak at 2,873–2,865 cm⁻¹ corresponding to shortening of an acyl chain was correlated with a lower concentration of fatty acid transporter at the outer membrane and (modified) lipid A.

Concerning protein folding, higher expression of chaperones *dnaJ* and *dnaK* was observed in Ag⁺-treated cells, whereas *groS* expression was less affected than the average (note that *groL* and *grpE* expression exhibited no specific regulatory pattern) (S1 Table). A marked decrease was observed in the expression of peptidylprolyl isomerase genes (*fklB*, *ppiB*, *surA* and *tig*), which function as protein folding chaperones for newly synthesised proteins. The bacterial response to low level of ionic silver was thus correlated with specific higher expression of chaperones involved in the rescue of misfolded proteins but not with those involved in folding newly formed proteins. Up-regulated expression of proteases (*lon*, *ftsH* and *htpX*) in the presence of ionic silver could be responsible for the degradation of unfolded/misfolded proteins. In the same way, expression of *clpA*, *P*, *X*, and *hslU*, *V* proteases did not follow the general trend of decreased expression. However, expression of *degP* involved in abnormal protein degradation in the periplasm nevertheless decreased.

As in the case of fatty acids, correlation analysis was performed to link variations in gene expression to protein structure-related sFTIR data (Fig 7). The up-regulated expression of protease *htpX* and chaperone *dnaJ* was identified as markers of the unordered random coil form of amide I (i.e., a positive correlation between expression regulation of *htpX* and *dnaJ* and the 1,681 cm⁻¹ wavenumber). Concerning amide II, the β-sheet characteristic wavenumbers (1,519 and 1,523 cm⁻¹) were not correlated with the expression of the previously mentioned protease and chaperone, but were strongly correlated with the L,D transpeptidase *ycfS*. The α-helix characteristic wavenumber 1,658 cm⁻¹ was also positively correlated with *ycfS*, while wavenumbers 1,654 and 1,650 cm⁻¹ were both positively correlated with the protease *degP*, and the L,D transpeptidase *ycbB*.

## Discussion

In this work, a multi-level characterization of E. coli response upon exposure to micromolar concentration of ionic silver was performed using an integrative biology approach to identify and connect different systemic effects at the transcriptional, physiological and biochemical levels.

First, it is important to note that the concentrations of AgNO₃ used in this study were sub-lethal. Although cell growth was affected to different extents at 6.5 and 8.5 µM of AgNO₃, after 3 hours of culture and contrary to a lethal concentration (e.g. a concentration of 100 µM of AgNO₃ [15]), a proportion of the E. coli population (at least 10⁷ cells/mL/ODU at 6.5 and 8.5 µM of AgNO₃) remained cultivable after removal of the silver. The effect of ionic silver on
<figure><img src="image_10.png" /><figcaption>**Fig 6.** Relevance networks resulting from the sPLSR approach between the Y matrix composed of selected sFTIR data in the region of fatty acids (3,100–2,800 cm⁻¹) and the X matrix corresponding to the expression of genes involved in fatty acid metabolism. The gene list was restricted to the 106 differentially-expressed genes annotated as involved in fatty acid and phospholipid metabolism and lipid transport (GO:0006631, GO:0006644 and GO:0006869). A threshold of 0.5 was used for relevant correlations. Green and red edges indicate positive and negative correlations respectively. Wavenumbers are represented as ellipses and genes as rectangles.</figcaption></figure>

doi:10.1371/journal.pone.0145748.g006

*E. coli* growth was shown to be, at least to some extent, linked to the accumulation of silver inside bacterial cells, as shown by NanoSIMS and elemental composition mapmaking. The overlap of silver with sulphur supported the hypothesis that Ag⁺ ions interact with intracellular proteins through the formation of silver-associated protein aggregates. The effect of ionic silver on proteins was further demonstrated (see below). To the best of our knowledge, this is the first time NanoSIMS has been used to probe silver intracellular uptake in prokaryotes, because previous studies were restricted to yeasts [22], daphnids [39, 40] and human cell lines [18, 40]. In fact, thanks to the excellent sensitivity of SIMS at nanoscale, it was possible to detect, localize, and quantify silver in small bacteria (through the Ag/CN ratio) at very low concentrations, i.e. when 5.0 µM AgNO₃ was added to the medium. Notably, it was shown that a tiny increase in intracellular silver content between cells treated with 5.0 or 6.5 µM of AgNO₃ could affect cell growth. This may be due to the phenomenon of silver "bioaccumulation", which explains the efficacy of the Ag⁺ ion even at low micromolar concentrations [2]. As reported previously [41], the amount of intracellular ionic silver was a better marker of the toxicity of silver nanoparticles than the quantity of Ag⁺ ions dissolved from these nanoparticles in the bacterial
<figure><img src="image_11.png" /><figcaption>**Fig 7.** Relevance networks resulting from the sPLSR approach between the Y matrix composed of selected sFTIR data in the region of amide I and amide II (1,700–1,480 cm⁻¹), and the X matrix corresponding to the expression of genes involved in protein folding. The gene list was restricted to the 104 differentially-expressed genes annotated as chaperone (GO:0051087), peptidase (GO:0008233), heat shock protein (GO:0031072), peptidyl-prolyl isomerase (GO:0003755) and protein disulfide isomerase (GO:0003756). A threshold of 0.5 was used for relevant correlations. Green and red edges indicate positive and negative correlations respectively. Wavenumbers are represented as ellipses and genes as rectangles.</figcaption></figure>

doi:10.1371/journal.pone.0145748.g007

environment. This corroborates the interest of NanoSIMS for detecting intracellular bioavailable silver. For bacterial cells treated with 8.5 µM AgNO₃, the Ag/CN ratio was lower, probably indicating intracellular silver release due to damaged and more “leaky” *E. coli* cell walls. This type of morphological alterations have already been described [1] in *E. coli* cells exposed to 10 µg/mL AgNO₃ (i.e., 92.7 µM) as well as in our earlier study [15]. It can also be hypothesised that an active system of transport was set up by cells to excrete silver [12].

A search for correlations between sFTIR data and intracellular silver content confirmed our previous conclusions obtained with the extracellular AgNO₃ concentration, in particular silver-mediated shortening of acyl chains in fatty acids of *E. coli* after exposure to sub-lethal concentrations of ionic silver (i.e., < 10 µM) [15]. However, in the present study, our insight into intracellular silver composition and the use of a multivariate statistical data analysis for identification of covariance (sPLSR) allowed more subtle conditions to be correctly probed and distinguished, thereby enhancing our initial approach, which only took the extracellular concentration of AgNO₃ into account and was based on the combination of principal component
analysis (PCA) and partial least squares (PLS) regression. In particular, variations in the con-
formation of unsaturated fatty acids were highlighted. In fact, in the presence of low concentra-
tions of intracellular silver (i.e., sub-lethal AgNO₃ concentrations) and in agreement with the
results of our earlier study [15], a higher “cis” conformation content was observed, indicating a
slight increase in bacterial membrane fluidity. It can be hypothesised that these changes in
membrane fluidity facilitate silver excretion and correspond to an active response by bacterial
cells. Many studies focused on the membrane unsaturated fatty acids of bacterial cells subjected
to various adverse environments have reported “cis-to-trans” isomerisation (for example, [42]
in *E. coli*; [43], [44] and [45] *Pseudomonas* or *Vibrio*), further illustrated by our group in *E. coli*
at a lethal concentration of silver (100 µM AgNO₃) [15]. However, an increase in “cis” confor-
mation was also reported by [46], in a study of the effects of sub-lethal concentrations of euge-
nol on the fatty acid profile of *E. coli* O157:H7 ATCC 43888 membranes. In a previous study
[47], the same authors also observed this type of modification by evaluating the stress induced
by a sub-lethal mix of ethanol and hexanal on the membrane fatty acid composition of *E. coli*
555. Alterations of the protein secondary structure (i.e., loss of α-helix structures) were clearly
demonstrated in the present study, while not detected at similar AgNO₃ concentrations in our
earlier work [15]. Lastly, new findings, such as the disruption of lipid esters, were also obtained
thanks to our optimised and more sensitive method.

The transcriptomic response was weak at 5.0 µM but massive and strong at 6.5 µM and
8.5 µM AgNO3 with a majority of down-regulated genes mainly involved in growth-related
functions. Down-regulated expression of transporter genes, transcriptional regulators and
ribosomal proteins was positively correlated with the decrease in the growth rate and in the
number of cultivable cells. Although in our study, the global transcriptomic response generally
differed from the copper-mediated one, we showed that the AgNO3 concentration and Ag/CN
ratio were positively correlated with up-regulation of genes specifically involved in the response
to copper stress. This confirms that mechanisms to cope with silver and copper share a primary
implication of the efflux transport *copA* in the resistance process and oxidation via the cuproux
oxidase CueO [16]. But we also showed that the ionic silver concentration was positively corre-
lated with up-regulated markers of the general stress response, cation transport and methio-
nine biosynthesis. The general response to stress contributes to the response to silver as shown
in *E. coli* cells exposed to Zn [48], ethanol and sodium hydroxide [42]. In addition, our integra-
tive approach revealed that within the general stress response, up-regulated expression of chap-
erones and proteases in the presence of silver ions was correlated with the increase in
disordered secondary structure of proteins.

Only global transcriptomic responses associated with growth arrest conditions were
strongly correlated with one another, indicating that the global variation in gene expression
reflects a growth rate effect rather than a direct effect of silver. Therefore, at concentrations of
6.5 µM and 8.5 µM of AgNO3 corresponding to non-growing cells, the observed genetic
response probably corresponded to the overlap between the direct specific transcriptional
responses to coping with ionic silver stress and the indirect general transcriptional response
linked to changes in the growth rate. In the future, it would thus be interesting to conduct a
study of *E. coli* response to ionic silver-mediated stress in chemostat culture at imposed growth
rates in order to remove the growth rate effect and specifically focus on the silver-mediated
effect [49].

Of note, under our experimental conditions, we did not find any clear specific up-regulated
expression of the DNA repair processes (excision repair, nucleases, recombination and error-
prone DNA polymerase) or of anaerobic respiration after exposure to silver ions, as previously
reported for silver nanoparticles [11, 17]. On the contrary, we found strong down-regulation of
genes involved in anaerobic respiration (nitrate/nitrite reductase, fumarate reductase enriched
in clusters 3 and 4, Fig 4). Consequently, the antibacterial effect of ionic silver appears to differ, at least to some extent, from the effect exerted by nano-sized silver, as recently reported [41]. Furthermore, the down-regulated expression of many oxidase-reductase genes was shown to be correlated with the presence of silver ions. This finding, added to the correlation found between the AgNO₃ concentration and up-regulation of methionine biosynthesis, strongly suggests the occurrence of oxidative stress in bacterial cells exposed to ionic silver [41, 50, 51]. Taken together, these results show that the response to ionic silver stress at the transcriptomic level corresponds to a particular combination of components involved in a range of different responses (growth rate, metal ion stress, general stress and oxidative stress).

The effect of ionic silver, i.e. reducing the amount of unsaturated fatty acids [15], was shown here to be not related to the transcriptomic regulation of the three key enzymes FabA, FabB and FabZ. However, we cannot exclude the possibility that ionic silver plays a role in regulating the synthesis (at the post-transcriptional level) and the activity (at the post-translational level) of these enzymes. Furthermore, the silver-induced shortening of the acyl chain length was shown to be correlated with the down-regulated expression of genes involved in the unsaturated fatty acid elongation pathway, fatty acid synthesis and transport, and also catalysing bio-synthesis and phosphate modification of lipid A. This result clearly shows that cell membrane is a key target of ionic silver. In addition, in an earlier study by our group, the cell ultrastructure of *E. coli* was shown to be damaged under Ag⁺-mediated stress, with possible detachment of the outer membrane at lethal and sub-lethal concentrations of ionic silver [15]. Using correlation analysis of changes in gene expression and protein structure-related sFTIR data, we are now able to provide a plausible explanation for these morphological changes. Indeed, the conformational shift (from α-helix to β-sheet) of amides I and II in the presence of ionic silver was correlated with gene expression regulation of L,D transpeptidases YcfS and YcbB. Interestingly, these L,D transpeptidases act on the structural integrity of the cell wall. In particular, YcfS contributes to the covalent attachment of peptidoglycan to the outer membrane whereas YcbB is involved in the synthesis of peptidoglycan cross-links [52]. ycfS expression was up-regulated in our experiments, while ycbB expression was less affected than the average of all other genes. Therefore, the specific regulation of the expression of YcfS and YcbB in the presence of ionic silver could be at play in the morphological changes observed. Therefore, a molecular approach aimed at modulating the expression of YcfS and YcbB can be recommended to control the effects of silver on *E. coli* when ionic silver is used in antibacterial coatings, for instance.

In conclusion, the multi-level characterization of *E. coli* exposed to Ag⁺-mediated stress was achieved using an integrative biology approach, based on the combination of physiological, biochemical and transcriptomic data sets. From an analytical point of view, the elemental composition mapmaking provided by NanoSIMS was used for the first time to unambiguously detect and localize intracellular silver within individual bacterial cells in interaction with cell components like proteins, due to NanoSIMS high sensitivity and sub-micrometer lateral resolution. The semi-quantitative data obtained with this powerful tool combined with the evaluation of biochemical composition of the cell at the single-cell scale thanks to sFTIR microspectroscopy was extremely useful in further explaining the physiological results. Finally, it was possible to define a genetic pattern of the ionic-silver induced response specifically dedicated to coping with intracellular silver but also associated with changes in growth and survival, cell membrane and protein composition and structure. As Gram-positive bacteria are generally said to be less sensitive to ionic silver than Gram-negative bacteria, notably due to differences in the structure of the cell wall, it would be interesting to use the same approach to explore the biological response of Gram-positive bacteria under silver stress.
**Supporting Information**

**S1 Fig.** (A) **FTIR raw spectra** (means of 20 spectra) in the 4,000–800 cm⁻¹ region recorded at a single-cell scale in *E. coli* grown without (control) or in the presence of ionic silver at sub-lethal concentrations. The regions of interest are identified in the figure. (B) Focus on the fatty acid region (3,100–2,800 cm⁻¹) with labels on the peaks attributed to the symmetric C-H stretching of—CH₃ and >CH₂ groups and to the up- and downshift of the asymmetric C-H stretching of >CH₂ bands. (C) Focus on the region including the bands attributed to C = O of esters, amides I and II and PO₂⁻ of nucleic acids and C-O stretching modes (1,800–1,200 cm⁻¹). The bands characteristic of C = O in lipid esters, amide I of peptide bonds and NH₂ of adenine are labeled.

(TIF)

**S2 Fig. Expression of fabA, fabZ and fabB in response to ionic-silver mediated stress.** Intensity (in log₂) reflects the level of gene expression in each growth condition, i.e. with 0, 5.0, 6.5 and 8.5 µM of AgNO₃ in the medium.

(TIF)

**S1 Table.** Fold changes in expression and p-values between different concentrations of AgNO₃ and the associated cluster number for all the genes.

(XLSX)

**Acknowledgments**

The research described in this paper built on results obtained at the French national synchrotron facility SOLEIL (Gif-sur-Yvette, France), using the SMIS beamline (proposal no 20090556). We are very grateful to Dr. Paul Dumas for the development of the SMIS beamline and for helpful discussions. The authors gratefully acknowledge the “BioPleasure” project of the French National Research Agency (ANR-07-BLAN-0196-01) for funding.

**Author Contributions**

Conceived and designed the experiments: CSB JNA IF FJ MCB MMB LG. Performed the experiments: CSB JNA IF FJ MMB LG. Analyzed the data: CSB IG BE JNA IF FJ MCB MMB LG. Wrote the paper: CSB IG BE JNA IF FJ MCB MMB LG.

**References**

1. Feng QL, Wu J, Chen GQ, Cui FZ, Kim TN, Kim JO. A mechanistic study of the antibacterial effect of silver ions on Escherichia coli and Staphylococcus aureus. Journal of biomedical materials research. 2000; 52(4):662-8. PMID: 11033548.

2. Lansdown ABG. Silver I: its antibacterial properties and mechanisms of action. J Wound Care 2002; 11 (4):125-30. PMID: 11998592

3. Matsumura Y, Yoshikata K, Kunisaki S, Tsuchido T. Mode of bactericidal action of silver zeolite and its comparison with that of silver nitrate. Appl Environ Microb. 2003; 69(7):4278-81. doi: 10.1128/Aem.69.7.4278-4281.2003 WOS:000184082100082.

4. Silver S, Phung LT, Silver G. Silver as biocides in burn and wound dressings and bacterial resistance to silver compounds. J Ind Microbiol Biot. 2006; 33(7):627-34. doi: 10.1007/s10295-006-0139-7 WOS:000239443800020.

5. Kierans M, Staines AM, Bennett H, Gadd GM. Silver Tolerance and Accumulation in Yeasts. Biol Met. 1991; 4(2):100-6. doi: 10.1007/Bf01135386 WOS:A1991FK55400005. PMID: 1878280

6. Yang HC, Pon LA. Toxicity of metal ions used in dental alloys: A study in the yeast Saccharomyces cerevisiae. Drug Chem Toxicol. 2003; 26(2):75-85. doi: 10.1081/Dct-120020403 WOS:000183020500001. PMID: 12816393
7. Lansdown ABG. Silver II: toxicity in mammals and how its products aid wound repair. J Wound Care. 2002; 11(5):173-7. PMID: 12055941

8. Klasen HJ. Historical review of the use of silver in the treatment of burns. I. Early uses. Burns: journal of the International Society for Burn Injuries. 2000; 26(2):117-30. PMID: 10716354.

9. Liau SY, Read DC, Pugh WJ, Furr JR, Russell AD. Interaction of silver nitrate with readily identifiable groups: relationship to the antibacterial action of silver ions. Letters in applied microbiology. 1997; 25 (4):279-83. PMID: 9351278.

10. Gordon O, Slenters TV, Brunetto PS, Villaruz AE, Sturdevant DE, Otto M, et al. Silver Coordination Polymers for Prevention of Implant Infection: Thiol Interaction, Impact on Respiratory Chain Enzymes, and Hydroxyl Radical Induction. Antimicrob Agents Ch. 2010; 54(10):4208-18. doi: 10.1128/Aac. 01830-09 WOS:000281907200023.

11. Radzig MA, Nadtochenko VA, Koksharova OA, Kiwi J, Lipasova VA, Khmel IA. Antibacterial effects of silver nanoparticles on gram-negative bacteria: influence on the growth and biofilms formation, mechanisms of action. Colloids and surfaces B, Biointerfaces. 2013; 102:300-6. doi: 10.1016/j.colSURFB.2012. 07.039 PMID: 23006569.

12. Holt KB, Bard AJ. Interaction of silver(I) ions with the respiratory chain of Escherichia coli: an electro-chemical and scanning electrochemistry microscopy study of the antimicrobial mechanism of micromolar Ag+. Biochemistry. 2005; 44(39):13214-23. doi: 10.1021/bi0508542 PMID: 16185089.

13. Inoue Y, Hoshino M, Takahashi H, Noguchi T, Murata T, Kanzaki Y, et al. Bactericidal activity of Ag-zeolite mediated by reactive oxygen species under aerated conditions. Journal of inorganic biochemistry. 2002; 92(1):37-42. PMID: 12230986.

14. Park HJ, Kim JY, Kim J, Lee JH, Hahn JS, Gu MB, et al. Silver-ion-mediated reactive oxygen species generation affecting bactericidal activity. Water Res. 2009; 43(4):1027-32. doi: 10.1016/j.watres.2008. 12.002 WOS:000264328600017. PMID: 19073336

15. Saulou C, Jamme F, Girbal L, Maranges C, Fourquaux I, Cocaign-Bousquet M, et al. Synchrotron FTIR microspectroscopy of Escherichia coli at single-cell scale under silver-induced stress conditions. Analytical and bioanalytical chemistry. 2013; 405(8):2685-97. doi: 10.1007/s00216-013-6725-4 PMID: 23354575.

16. Wu MY, Suryanarayanan K, van Ooij WJ, Oerther DB. Using microbial genomics to evaluate the effectiveness of silver to prevent biofilm formation. Water science and technology: a journal of the International Association on Water Pollution Research. 2007; 55(8-9):413-9. PMID: 17547012.

17. Du H, Lo TM, Sitompul J, Chang MW. Systems-level analysis of Escherichia coli response to silver nanoparticles: the roles of anaerobic respiration in microbial resistance. Biochemical and biophysical research communications. 2012; 424(4):657-62. doi: 10.1016/j.bbrc.2012.06.134 PMID: 22771582.

18. Gutleb AC, Freitas J, Murk AJ, Verhaegen S, Ropstad E, Udelhoven T, et al. NanoSIMS50-a powerful tool to elucidate cellular localization of halogenated organic compounds. Analytical and bioanalytical chemistry. 2012; 404(9):2693-8. doi: 10.1007/s00216-012-6066-8 WOS:000310339400012. PMID: 22566200

19. Hoppe P, Cohen S, Meibom A. NanoSIMS: Technical Aspects and Applications in Cosmochemistry and Biological Geochemistry. Geostand Geoanal Res. 2013; 37(2):111-54. doi: 10.1111/j.1751-908X. 2013.00239.x WOS:000319514900002.

20. Mueller CW, Weber PK, Kilburn MR, Hoeschen C, Kleber M, Pett-Ridge J. Advances in the Analysis of Biogeochemical Interfaces: NanoSIMS to Investigate Soil Microenvironments. Adv Agron. 2013; 121:1-46. doi: 10.1016/B978-0-12-407685-3.00001-3 WOS:000322846500001.

21. Blattner FR. The Complete Genome Sequence of Escherichia coli K-12. Science. 1997; 277 (5331):1453-62. doi: 10.1126/science.277.5331.1453 PMID: 9278503

22. Saulou C, Despax B, Raynaud P, Zanna S, Syeux A, Marcus P, et al. Plasma-Mediated Nanosilver-Organosilicon Composite Films Deposited on Stainless Steel: Synthesis, Surface Characterization, and Evaluation of Anti-Adhesive and Anti-Microbial Properties on the Model Yeast Saccharomyces cerevisiae. Plasma Process Polym. 2012; 9(3):324-38. doi: 10.1002/ppap.201100033 WOS:000301225300011.

23. Lechene C, Hillion F, McMahon G, Benson D, Kleinfeld AM, Kampf JP, et al. High-resolution quantitative imaging of mammalian and bacterial cells using stable isotope mass spectrometry. J Biol 2006; 5 (6):20-50. PMID: 17010211

24. R Core Team (2012). R: A language and environment for statistical computing. R Foundation for Statistical Computing V, Austria. ISBN 3-900051-07-0, URL http://www.R-project.org/.

25. Gautier L, Cope L, Bolstad BM, Irizarry RA. affy—analysis of Affymetrix GeneChip data at the probe level. Bioinformatics. 2004; 20(3):307-15. doi: 10.1093/bioinformatics/btg405 PMID: 14960456.
26. Smyth GK. Limma: linear models for microarray data. In Bioinformatics and Computational Biology Solutions using R and Bioconductor. Gentleman R., Carey V., Huber W., Irizarry R., and Dudoit S., eds. (Springer, New York). 2005:397-420.

27. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003; 4(2):249-64. doi: 10.1093/biostatistics/4.2.249 PMID: 12925520

28. Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Statistical applications in genetics and molecular biology. 2004; 3:Article3. doi: 10.2202/1544-6115.1027 PMID: 16646809.

29. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B. 1995; 57:289-300.

30. Alexa A, Rahnenfuhrer J. topGO: Enrichment analysis for Gene Ontology. In: 2.2.0 Rpv, editor. 2010.

31. Le Cao KA, Rossouw D, Robert-Granie C, Besse P. A Sparse PLS for Variable Selection when Integrating Omics Data. Statistical applications in genetics and molecular biology. 2008; 7(1). WOS:000261334800001.

32. Le Cao KA, Gonzalez I, Dejean S. integrOmics: an R package to unravel relationships between two omics datasets. Bioinformatics. 2009; 25(21):2855-6. doi: 10.1093/bioinformatics/btp515 PMID: 19706745; PubMed Central PMCID: PMC2781751.

33. Gonzalez I, Le Cao KA, Davis MD, Dejean S. Insightful graphical outputs to explore relationships between two 'omics' data sets. BioData Mining. 2013; 5:19.

34. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003; 13(11):2498-504. doi: 10.1101/Gr.1239303 WOS:000186357000016. PMID: 14597658

35. Eybe T, Bohn T, Audinot JN, Udelhoven T, Cauchie HM, Migeon HN, et al. Uptake visualization of deltamethrin by NanoSIMS and acute toxicity to the water flea Daphnia magna. Chemosphere. 2009; 76(1):134-40. doi: 10.1016/j.chemosphere.2009.02.005 WOS:000266897200020. PMID: 19278714

36. Movasaghi Z, Rehman S, Rehman IU. Fourier transform infrared (FTIR) spectroscopy of biological tissues. Appl Spectrosc Rev. 2008; 43(2):134-79. doi: 10.1080/05704920701829043 WOS:000252893900004.

37. Heath RJ, Rock CO. Roles of the FabA and FabZ beta-hydroxyacyl-acyl carrier protein dehydratases in Escherichia coli fatty acid biosynthesis. The Journal of biological chemistry. 1996; 271(44):27795-801. PMID: 8910376.

38. Magnuson K, Jackowski S, Rock CO, Cronan JE Jr. Regulation of fatty acid biosynthesis in Escherichia coli. Microbiological reviews. 1993; 57(3):522-42. PMID: 8246839; PubMed Central PMCID: PMC372925.

39. Georgantzopoulou A, Balachandran YL, Rosenkranz P, Dusinska M, Lankoff A, Wojewodzka M, et al. Ag nanoparticles: size- and surface-dependent effects on model aquatic organisms and uptake evaluation with NanoSIMS. Nanotoxicology. 2013; 7(7):1168-78. doi: 10.3109/17435390.2012.715312 WOS:000325223500001. PMID: 22834480

40. Audinot JN, Georgantzopoulou A, Piret JP, Gutleb AC, Dowsett D, Migeon HN, et al. Identification and localization of nanoparticles in tissues by mass spectrometry. Surf Interface Anal. 2013; 45(1):230-3. doi: 10.1002/Sia.5099 WOS:000319048700059.

41. Ivask A, ElBadawy A, Kaweeteerawat C, Boren D, Fischer H, Ji ZX, et al. Toxicity Mechanisms in Escherichia coli Vary for Silver Nanoparticles and Differ from Ionic Silver. Acs Nano. 2014; 8(1):374-86. doi: 10.1021/Nn4044047 WOS:000330542900037. PMID: 24341736

42. Moen B, Janbu AO, Langsrud S, Langsrud O, Hobman JL, Constantinidou C, et al. Global responses of Escherichia coli to adverse conditions determined by microarrays and FT-IR spectroscopy. Canadian journal of microbiology. 2009; 55(6):714-28. doi: 10.1139/w09-016 PMID: 19767843.

43. Loffeld B, Keweloh H. cis/trans Isomerization of unsaturated fatty acids as possible control mechanism of membrane fluidity in Pseudomonas putida P8. Lipids. 1996; 31(8):811-5. doi: 10.1007/Bf02522976 WOS:A1996VD35200003. PMID: 8869883

44. Heipieper HJ, Meinhardt F, Segura A. The cis-trans isomerase of unsaturated fatty acids in Pseudomonas and Vibrio: biochemistry, molecular biology and physiological function of a unique stress adaptive mechanism. FEMS microbiology letters. 2003; 229(1):1-7. PMID: 14659535.

45. Hartig C, Loffhagen N, Harms H. Formation of trans fatty acids is not involved in growth-linked membrane adaptation of Pseudomonas putida. Appl Environ Microbiol. 2005; 71(4):1915-22. doi: 10.1128/AEM.71.4.1915-1922.2005 PMID: 15812020; PubMed Central PMCID: PMC1082513.

46. Di Pasqua R, Hoskins N, Betts G, Mauriello G. Changes in membrane fatty acids composition of microbial cells induced by addiction of thymol, carvacrol, limonene, cinnamaldehyde, and eugenol in the
growing media. J Agr Food Chem. 2006; 54(7):2745-9. doi: 10.1021/Jf052722I
WOS:000236628500042.

47. Patrignani F, Lucci L, Belletti N, Gardini F, Guerzoni AE, Lanciotti R. Effects of sub-lethal concentrations of hexanal and 2-(E)-hexenal on membrane fatty acid composition and volatile compounds of Listeria monocytogenes, Staphylococcus aureus, Salmonella enteritidis and Escherichia coli. Int J Food Microbiol. 2008; 123(1-2):1-8. doi: 10.1016/j.ijfoodmicro.2007.09.009 WOS:000254982100001. PMID: 18055050

48. Lee LJ, Barrett JA, Poole RK. Genome-wide transcriptional response of chemostat-cultured Escherichia coli to zinc. Journal of bacteriology. 2005; 187(3):1124-34. doi: 10.1128/JB.187.3.1124-1134. 2005 PMID: 15659689; PubMed Central PMCID: PMC545701.

49. Ishii N, Nakahigashi K, Baba T, Robert M, Soga T, Kanai A, et al. Multiple high-throughput analyses monitor the response of *Escherichia coli* to perturbations. Science. 2007; 316(5824):593-7. doi: 10.1126/science.1132067 PMID: 17379776.

50. Flatley J, Barrett J, Pullan ST, Hughes MN, Green J, Poole RK. Transcriptional responses of Escherichia coli to S-nitrosoglutathione under defined chemostat conditions reveal major changes in methionine biosynthesis. Journal of Biological Chemistry. 2005; 280(11):10065-72. doi: 10.1074/jbc.M410393200 WOS:000227559600041. PMID: 15647275

51. Partridge JD, Scott C, Tang Y, Poole RK, Green J. Escherichia coli transcriptome dynamics during the transition from anaerobic to aerobic conditions. Journal of Biological Chemistry. 2006; 281(38):27806-15. doi: 10.1074/jbc.M603450200 WOS:000240534400016. PMID: 16857675

52. Magnet S, Dubost L, Marie A, Arthur M, Gutmann L. Identification of the L,D-transpeptidases for peptidoglycan cross-linking in Escherichia coli. Journal of bacteriology. 2008; 190(13):4782-5. doi: 10.1128/JB.00025-08 PMID: 18456808; PubMed Central PMCID: PMC2446776.

53. McLean S, Bowman LAH, Sanguinetti G, Read RC, Poole RK. Peroxynitrite Toxicity in Escherichia coli K12 Elicits Expression of Oxidative Stress Responses and Protein Nitration and Nitrosylation. Journal of Biological Chemistry. 2010; 285(27):20724-31. doi: 10.1074/jbc.M109.085506 WOS:000279228600036. PMID: 20427277

54. Jozefczuk S, Klie S, Catchpole G, Szymanski J, Cuadros-Inostroza A, Steinhauser D, et al. Metabolomic and transcriptomic stress response of Escherichia coli. Mol Syst Biol. 2010; 6. Artn 364 doi: 10.1038/Msb.2010.18 WOS:000278575700001.

55. Kershaw CJ, Brown NL, Constantinidou C, Patel MD, Hobman JL. The expression profile of Escherichia coli K-12 in response to minimal, optimal and excess copper concentrations. Microbiol-Sgm. 2005; 151:1187-98. doi: 10.1099/mic.0.27650-0 WOS:000228509100017.

56. Helbig K, Grosse C, Nies DH. Cadmium toxicity in glutathione mutants of Escherichia coli. Journal of bacteriology. 2008; 190(15):5439-54. doi: 10.1128/JB.00272-08 PMID: 18539742; PubMed Central PMCID: PMC2493267.

57. Fantino JR, Py B, Fontecave M, Barras F. A genetic analysis of the response of Escherichia coli to cobalt stress. Environ Microbiol. 2010; 12(10):2846-57. doi: 10.1111/j.1462-2920.2010.02265.x WOS:000282375400017. PMID: 20545747