# Metabolic modulations of Pseudomonas graminis in response to H2O2 in cloud water

Nolwenn Wirgot, Marie Lagree, Mounir Traïkia, Ludovic Besaury, Pierre Amato, Isabelle Canet, Martine Sancelme, Cyril C. Jousse, Binta Diémé, Bernard Lyan, et al.

► To cite this version:

Nolwenn Wirgot, Marie Lagree, Mounir Traïkia, Ludovic Besaury, Pierre Amato, et al. Metabolic modulations of Pseudomonas graminis in response to H2O2 in cloud water. Scientific Reports, 2019, 9, pp.12799. 10.1038/s41598-019-49319-2. hal-02279722

HAL Id: hal-02279722

https://hal.science/hal-02279722

Submitted on 24 Sep 2019

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
OPEN

# Metabolic modulations of *Pseudomonas graminis* in response to H₂O₂ in cloud water

Received: 3 April 2019

Accepted: 26 July 2019

Published online: 05 September 2019

Nolwenn Wirgot<fnref n="1" />, Marie Lagrée<fnref n="1" />,², Mounir Traïkia<fnref n="1" />,², Ludovic Besaury<fnref n="1" />, Pierre Amato<fnref n="1" />,³, Isabelle Canet<fnref n="1" />, Martine Sancelme<fnref n="1" />, Cyril Jousse<fnref n="1" />,², Binta Diémé<fnref n="1" />,², Bernard Lyan² & Anne-Marie Delort<fnref n="1" />,²

In cloud water, microorganisms are exposed to very strong stresses especially related to the presence of reactive oxygen species including H₂O₂ and radicals, which are the driving force of cloud chemistry. In order to understand how the bacterium *Pseudomonas graminis* isolated from cloud water respond to this oxidative stress, it was incubated in microcosms containing a synthetic solution of cloud water in the presence or in the absence of H₂O₂. P. graminis metabolome was examined by LC-MS and NMR after 50 min and after 24 hours of incubation. After 50 min, the cells were metabolizing H₂O₂ while this compound was still present in the medium, and it was completely biodegraded after 24 hours. Cells exposed to H₂O₂ had a distinct metabolome as compared to unexposed cells, revealing modulations of certain metabolic pathways in response to oxidative stress. These data indicated that the regulations observed mainly involved carbohydrate, glutathione, energy, lipid, peptides and amino-acids metabolisms. When cells had detoxified H₂O₂ from the medium, their metabolome was not distinguishable anymore from unexposed cells, highlighting the capacity of resilience of this bacterium. This work illustrates the interactions existing between the cloud microbial metabolome and cloud chemistry.

Clouds are oxidative environments with a redox potential of up to more than 200 mV<fnref n="1" />. The oxidant capacity of the atmosphere is often evaluated based on the concentration of H₂O₂ in the gas and aqueous phases of clouds: this is highly variable depending on the season, day-night cycles, and the level of air pollution². According to Deguillaume et al.<fnref n="1" />, the concentration of H₂O₂ can vary from 0 to 58 µM in cloud waters collected at the puy de Dôme Mountain station. H₂O₂ is the main source of hydroxyl radicals in clouds³; this oxidizes the organic matter (OM) into molecules such as organic acids or aldehydes and, ultimately, to CO₂. Short chain organic acids and aldehydes account for most of the identified carbon compounds in cloud waters (about 10% of the OM). The most abundant are acetate, formate, succinate, malonate, oxalate and formaldehyde<fnref n="1" />,⁴. More recently, amino acids have been quantified to represent about 9% of the OM⁵ in clouds. Consistently, carboxylic acids, aldehydes and amino acids are also present at high concentrations in fog and aerosols⁴,⁶,⁷. The precise composition of the rest of the OM in the atmosphere is still not well defined. Recent investigation using FT-ICR-MS approaches provided a global description based on exact mass formulae and H/C/O/N molecular ratios⁸,⁹.

Although cloud water is probably a stressful environment for living organisms, in particular due to the presence of reactive oxygen species (ROS) such as H₂O₂, it was shown that cloud microorganisms largely survive exposure to H₂O₂ at its natural atmospheric concentration¹⁰,¹¹ and that some remain alive and metabolically active in clouds¹²⁻¹⁵. ATP (Adenosine TriPhosphate) is currently used as a proxy for metabolic activity and it was quantified in cloud waters¹⁶⁻¹⁹. Different studies have shown that airborne microorganisms are able to bio-transform mono and dicarboxylic acids (oxalic, malonic, succinic, glutaric, adipic, pimelic, pinic acids, acetic, formic acids) and also formaldehyde and methanol¹⁶,¹⁹⁻²⁵. The measured biotransformation rates were shown to be competitive with rates involving radical reactions. It was thus proposed that microbial activity could be an alternative route to radical chemistry³,⁴,¹⁵,²⁶,²⁷. In addition, cloud microorganisms impact H₂O₂ concentration by

<footnote n="1">¹Université Clermont Auvergne, CNRS, Sigma-Clermont, Institut de Chimie de Clermont-Ferrand, F-63000, Clermont-Ferrand, France. ²Plateforme d'Exploration du Métabolisme, Université Clermont Auvergne & I.N.R.A site de Theix, Clermont-Ferrand, France. Correspondence and requests for materials should be addressed to A.-M.D. (email: A-marie.delort@uca.fr)</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** (a) Temporal evolution of H₂O₂ concentration in the presence (circles) or absence of *P. graminis* 13b-3 (triangles). (b) Evolution of ATP concentration in *P. graminis* 13b-3 incubated with H₂O₂ (triangles) or without H₂O₂ (squares). Error bars are standard errors of the means of three replicates. Where error bars are not visible, they are smaller than the symbol.</figcaption></figure>

detoxifying it into H₂O and O₂ through catalases activity, and consequently they can indirectly limit the efficiency of radical chemistry¹¹,¹⁹.

The ability to manage ROS is a permanent challenge for aerobic organisms that use molecular oxygen as the terminal electron acceptor. Reactive by-products of oxygen, such as superoxide anion radical, hydroxyl radicals and hydrogen peroxide, are generated in aerobically grown cells. Oxidative stress is defined as an imbalance in favor of the prooxidants and disfavoring the antioxidants²⁸; resulting in cellular damages damage (DNA, lipids, proteins and carbohydrates). To avoid such damages, microorganisms deploy various strategies involving specific enzyme like catalases, superoxide dismutases and the production of oxidant scavengers (glutathione, pigments, etc.). In addition, bacteria exposed to ROS are suspected to reorganize their entire functioning under oxidative stress. Modulation of the global cell metabolism can be observed using metabolomics. This powerful tool has been used for investigating oxidative stress response in plants²⁹ and bacteria³⁰,³¹.

In a recent study, we observed that H₂O₂ modulated the global metabolic functioning of microorganisms in cloud water, as shown by strong correlations between H₂O₂ and cellular ATP concentrations¹¹. The objective of this work was to go further in elucidating precisely the metabolic modifications involved, with the aim to understand how cloud microorganisms may interact with cloud chemistry under very oxidative conditions. Since microorganisms use or produce organic compounds modulations of their metabolism can potentially influence the reactivity of organic compounds in clouds.

We used a metabolomics approach to investigate the response to H₂O₂ exposure of *Pseudomonas graminis* strain PDD-13b-3, a bacterium that we previously isolated from cloud water and whose genome is sequenced³². *Pseudomonas* species are among the most frequent bacteria observed in clouds, and they account for some of the most active cells in these environments¹³,¹⁸.

The combination of NMR and LC-MS analyses allowed identifying and discriminating the metabolites and the metabolic pathways impacted by this stress. Because our main objective was to highlight the modifications of the metabolism when the impact of H₂O₂ was maximum, we monitored the ATP concentration over time, and we sampled cells for metabolomics study at the time corresponding to the lowest ATP level (t = 50 min), when H₂O₂ was not totally degraded yet. We also investigated the metabolic status of the cells after 24 hours, when cells had recovered initial ATP concentration and no H₂O₂ was left. The interactions between the bacterial metabolome and cloud chemistry implied under these conditions are further discussed.

## Results

**Monitoring of hydrogen peroxide and ATP concentrations.** H₂O₂ degradation by *Pseudomonas graminis* PDD-13b-3 was monitored during the incubation by measuring its concentration until it was exhausted (Fig. 1a); H₂O₂ was completely biodegraded after two hours, and the corresponding degradation rate was on the same order of magnitude as compared to Wirgot *et al.*¹¹ (in the range of 10⁻¹⁷ mol.min⁻¹ cell⁻¹).

The ATP concentration was closely monitored during the first 13 hours of incubation, and after 24 hours (Fig. 1b). The addition of H₂O₂ (time “0”) caused a rapid decrease of ATP concentration (170 pM to 25 pM in approximately 30 min), and this lasted until H₂O₂ was totally degraded (represented on the graph with a vertical bar around 120 min). After 400 min, *P. graminis* 13b-3 recovered an ATP concentration that was close to the initial value (approximately 150 pM). After 24 hours, the ATP concentration reached a concentration of 125 pM in all samples, indistinctly between cells exposed or unexposed to hydrogen peroxide.

To assess these two different response phases, the cell extract contents (i.e. the metabolomes) after 50 min and 24 h of incubation in the presence and in the absence of added H₂O₂ were compared.

**Impact of H₂O₂ exposure on *P. graminis* metabolome.** *Effect of exposure time on metabolic profile.* First, statistical analyses (PCA, Figs S1 and S2; PLS-DA, Figs 2 and 3) were done without any preconceptions
<figure><img src="image_3.png" /><figcaption>**Figure 2.** PLS-DA analyses (UV scaling). (a) NMR and (b) LC-MS example of the positive ionization mode. Blue circles correspond to the samples incubated in the presence of H₂O₂ and extracted at 50 min; red circles correspond to the samples incubated in the absence of H₂O₂ and extracted at 50 min; green circles correspond to the samples extracted at 24 hours with and without H₂O₂.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 3.** PLS-DA analyses. Score plots (UV scaling) obtained for bacteria metabolomes extracted at 50 min of incubation in the presence (blue circles) or in the absence (red circles) of H₂O₂: (a) NMR and (b) LC-MS (example of the positive ionization mode).</figcaption></figure>

on the dataset. Therefore, all samples i.e., extracted at 50 min and 24 h incubated in the presence/absence of H₂O₂, respectively, were considered, in order to highlight potential discriminations between the samples and to identify which compounds caused these significant differences. PLS-DA analyses were performed on LC-MS and NMR data from (i) cells incubated with H₂O₂ at 50 min, (ii) cells incubated without H₂O₂ at 50 min, (iii) cells incubated with H₂O₂ after 24 h of incubation and (iv) cells incubated without H₂O₂ at 24 h. Three classes are observed where samples at 24 h are grouped in a single class (Fig. 2); the model is valid and robust. In addition, a univariate statistical test (Wilcoxon-Mann-Whitney) based on relative intensities of ions and buckets for LC-MS and NMR analyses was carried out to confirm this grouping of samples at 24 h. At 50 min, most of buckets or ions were statistically different reflecting a majority of p-values strictly inferior to 0.05. The situation at 24 h was completely different, and the majority of data were not discriminating with p-value strictly greater than 0.05 (Fig. S3).

*Identification of key metabolites, markers of H₂O₂ exposure.* Based on statistical analyses, it is clear that the metabolome of cells incubated with or without H₂O₂ were different at time 50 min, while there is no significant difference at time 24 h between the two groups of cells (exposed or not to the oxidant). To highlight the key metabolites involved in the response to the presence of hydrogen peroxide, we focused our study on the sample extracted at 50 min. In order to identify key markers of H₂O₂ exposure, PLS-DA was done on data from samples
collected at 50 min, incubated or not in the presence of H₂O₂. A clear discrimination was observed between the two groups confirming that H₂O₂ had a significant effect on bacteria (Figs 3 and S4). Finally, using UV scaling lists of VIP scores 338 discriminating buckets (NMR), 119 and 79 discriminating ions (LC-MS, positive and negative mode) (VIP value >1) were obtained. In addition to the selection based on the value of VIP score, a second selection was done based on the results of the Wilcoxon Mann Whitney test (p-value < 0.05). Finally, 229 buckets, 93 (positive mode) and 76 (negative mode) ions were kept for the identification step. Table 1 summarizes the 60 VIP identified metabolites. Ratios S (Stressed: in the presence of H₂O₂)/R (Reference: in the absence of H₂O₂) were calculated for each compound based on relative intensities (ions and/or buckets) to evaluate if the compound was more/less present in the samples with H₂O₂.

The use of two different and complementary analytical techniques (NMR and LC-MS, see Tables S1, S2, S3) allowed detecting various metabolites, which are key markers of the response of *P. graminis* 13b-3 to hydrogen peroxide stress (Table 1): amino-acids, saccharides, amines, lipids, peptides and miscellaneous compounds. NMR analysis allowed detecting most of “miscellaneous” metabolites such as AMP, ATP, NADH, betaine, acetate, lactate, succinate, citrate, β-hydroxybutyrate, pipecolic acid, inosinic acid, nicotinamide ribotide, most sugars and also glutathione. LC-MS allowed the detection of lipids, some amino acids, peptides and others compounds such as glycerol-3P, spermidine and xylose. Some metabolites were detected with both NMR and LC-MS analyses. It is the case for amino-acids (leucine, valine, alanine, methionine, tyrosine, and glutamate), saccharide (trehalose) and miscellaneous compounds (pantothenate and UMP).

Discussion

Our results clearly show the metabolic response of *Pseudomonas graminis* 13b-3 to oxidative conditions such as existing in clouds.

First, metabolic pathways of *Pseudomonas graminis* are highly impacted under hydrogen peroxide stress. These impacted pathways (Figs 4–6) were mapped using Kegg data base along with the strain's genome annotation³²,³³. Carbohydrates pathway is displayed in Fig. 4. Carbohydrates are important sources of energy in living cells; their metabolism was highly impacted by the presence of H₂O₂. This metabolism is complex showing interconversion of various poly-, oligo- and di- and mono-saccharides. Glycogen and trehalose (as well as its derivative trehalose-6 phosphate) concentrations were decreased in the presence of H₂O₂, while those of maltose, glucose, cellobiose and fructose were increased. These metabolic changes showed that this bacterium accelerates the bio-degradation of its storage compounds (glycogen and trehalose) into maltose and glucose. Glucose and maltose entered the glycolysis and formed fructose from fructose 6-phosphate.

Soluble sugars, especially glucose and fructose, were involved in the responses to a number of stresses including oxidative stress due to the presence of ROS. It has been suggested earlier that sugar-signalling and sugar-modulated gene expressions are involved in the control of oxidative stress expression³⁴.

Carboxylic acid concentrations were modified in the presence of H₂O₂: lactate, citrate and succinate were found at higher amounts while acetate concentration is decreased (Fig. 5). This result is consistent from that observed in *E. coli* exposed to H₂O₂ where fumaric acid, succinic and malic acids are also more abundant³⁰.

Three classes of lipids were detected by LC-MS analysis (palmitoleic acid, three lyso PE (lyso-phosphoryl-ethanolamine) and two lyso PC (lyso-phosphoryl-choline)). Palmitoleic acid is part of fatty acid metabolism, while lyso PC and lyso PE are part of glycerophospholipid metabolism (Fig. 5). Except for one type of lyso PC, all the VIP scores were high (average around 1.80), so lipids were strongly discriminating metabolites in the expression of oxidative stress. In addition, the low values of S/R ratios measured at 50 min (average around 0.40) for all lipids indicated that these metabolite concentrations were lower in the presence of hydrogen peroxide (Table 1). In parallel, the concentration of glycerol and glycerol-3-phosphate, two intermediates of the phospholipid metabolism, were decreased under oxidative conditions. This observation is consistent with the fact that lipids are major targets of oxidative stress. Free radicals can attack directly polyunsaturated fatty acids in membranes and initiate lipid peroxidation. A primary effect of lipid peroxidation is a decrease in membrane fluidity, which alters membrane properties and can disrupt membrane-bound proteins significantly³⁵. Membrane damage disrupts cellular integrity, causing loss of the proton motive force that drives ATP formation and the uptake of many nutrients³⁶. However, further studies should be conducted in the future to explore in more detail lipid metabolism using different extraction and HPLC conditions. These will give more information about the oxidative state of lipids under oxidative conditions.

Energy metabolism was strongly impacted by the presence of hydrogen peroxide; at 50 min, when bacteria were exposed to the oxidant, the concentrations of adenosine, monophosphate (AMP) and adenosine triphosphate (ATP) were decreased (Table 1, Fig. 5) as well as that of inosinic acid, an intermediate of ATP and AMP synthesis. This observation is consistent with the results shown in Fig. 1 and also in Wirgot *et al*¹¹. Indeed, at 50 min the ATP concentration was about five times lower in the sample incubated with hydrogen peroxide than in the reference. The present study showed a decrease of 2.5-fold (S/R = 0.4). ATP is listed among the highest VIP metabolites (2.1 score). A few studies also reported a decrease of ATP concentrations in cells exposed to H₂O₂; these include *Streptomyces pneumoniae*³⁷ and cultured mammalian and plant cells³⁸⁻⁴⁰.

NADH concentration was also decreased in the presence of H₂O₂, which could result from a tuning of the cell redox potential (NAD⁺/NADH, NADP⁺/NADPH ratios). As H₂O₂ impacts glutathione metabolism, glycolysis, TCA cycle and DNA repair system, it also indirectly impacted these NAD⁺/NADH, NADP⁺/NADPH co-enzymes¹¹,⁴¹,⁴². NADH decrease could also have resulted from modulations for rather favoring the production of NADPH, involved in many anti-oxidant enzymatic reactions, at the expense of NADH as reported in *Pseudomonas fluorescens*³¹,⁴².

A large number of amino-acids were identified as VIP metabolites in this study (Table 1, Fig. 5). The concentrations of 75% of the amino-acids detected (isoleucine, leucine, valine, lysine, arginine, methionine, aspartate and cyclic compounds like tryptophan, phenylalanine, and tyrosine) were increased at 50 min under oxidative
<table><thead><tr><th rowspan="2">Identified VIP metabolite</th><th colspan="2">NMR</th><th colspan="3">LC-MS</th><th rowspan="2">ID metrics (class)</th></tr><tr><th>VIP Score</th><th>Ratio S/R</th><th>Ionization mode</th><th>VIP Score</th><th>Ratio S/R</th></tr></thead><tbody><tr><td colspan="7"><strong>Amino acids</strong></td></tr><tr><td>Isoleucine</td><td>1.82-1.89</td><td>1.61-1.67</td><td>Pos</td><td>2.40</td><td>2.47</td><td>II (9)</td></tr><tr><td>Leucine</td><td>1.11-1.43</td><td>1.23-1.32</td><td></td><td></td><td></td><td></td></tr><tr><td>Valine</td><td>1.93-2.10</td><td>4.35-4.55</td><td>- Neg</td><td>-2.15</td><td>4.99-5.26</td><td>I (9)</td></tr><tr><td>Lysine</td><td>1.12-1.03</td><td>1.12-1.13</td><td></td><td></td><td></td><td>(9)</td></tr><tr><td>Arginine</td><td>1.16-1.49</td><td>1.16-1.47</td><td></td><td></td><td></td><td>(7)</td></tr><tr><td>Proline</td><td>1.23-1.49</td><td>0.75-0.63</td><td></td><td></td><td></td><td>(3)</td></tr><tr><td>Glutamate</td><td>1.43-1.39</td><td>0.68-0.71</td><td>Pos - Neg</td><td>1.12-1.48</td><td>0.72-0.70</td><td>I (5)</td></tr><tr><td>Methionine</td><td>1.33</td><td>1.22</td><td>Pos - Neg</td><td>1.83-1.60</td><td>1.41-1.52</td><td>I (9)</td></tr><tr><td>Alanine</td><td>1.81</td><td>0.88</td><td>Pos</td><td>1.33</td><td>0.72</td><td>II (9)</td></tr><tr><td>Beta-Alanine</td><td>1.10</td><td>1.41</td><td></td><td></td><td></td><td>(5)</td></tr><tr><td>Aspartate</td><td>1.03-1.13</td><td>1.11-1.16</td><td></td><td></td><td></td><td>(9)</td></tr><tr><td>Asparagine</td><td>1.21</td><td>0.82</td><td></td><td></td><td></td><td>(5)</td></tr><tr><td>Phenylalanine</td><td></td><td></td><td>Pos - Neg</td><td>2.57-2.13</td><td>2.29-2.08</td><td>I</td></tr><tr><td>Tyrosine</td><td>1.40-1.53</td><td>1.44-1.52</td><td>Neg</td><td>1.58</td><td>2.10</td><td>I (9)</td></tr><tr><td>Glycine</td><td></td><td>Detected but not VIP</td><td>Pos</td><td>1.73</td><td>1.38</td><td>I (9)</td></tr><tr><td>Citrulline</td><td>1.01-2.00</td><td>0.84-0.81</td><td></td><td></td><td></td><td>(5)</td></tr><tr><td>Tryptophan</td><td>2.17</td><td>6.98</td><td></td><td></td><td></td><td>(5)</td></tr><tr><td colspan="7"><strong>Saccharides</strong></td></tr><tr><td>Trehalose*</td><td>1.74-1.69</td><td>0.85-0.85</td><td>Pos</td><td>1.96</td><td>0.74</td><td>II (9)</td></tr><tr><td>Trehalose 6-P*</td><td></td><td></td><td>Pos - Neg</td><td>2.27-2.04</td><td>0.52-0.54</td><td>II (9)</td></tr><tr><td>Cellobiose*</td><td>1.82-1.01</td><td>1.26-1.21</td><td>Pos</td><td>1.96</td><td>0.74</td><td>II (9)</td></tr><tr><td>Cellobiose 6-P*</td><td></td><td></td><td>Pos - Neg</td><td>2.27-2.04</td><td>0.52-0.54</td><td>II (9)</td></tr><tr><td>Glycogen</td><td>1.80-1.76</td><td>0.65-0.84</td><td></td><td></td><td></td><td>(9)</td></tr><tr><td>Alpha-glucose</td><td>1.71-1.69</td><td>1.52-2.19</td><td></td><td></td><td></td><td>(9)</td></tr><tr><td>Fructose</td><td>1.15-1.22</td><td>1.59-1.79</td><td></td><td></td><td></td><td>(2)</td></tr><tr><td>Maltose</td><td>1.40</td><td>1.65</td><td></td><td></td><td></td><td>(2)</td></tr><tr><td>Xylose</td><td></td><td></td><td>Neg</td><td>1.54</td><td>0.45</td><td>II (9)</td></tr><tr><td colspan="7"><strong>Amines</strong></td></tr><tr><td>Putrescine</td><td>1.16-1.23</td><td>1.32-1.16</td><td></td><td></td><td></td><td>(7)</td></tr><tr><td>Carnitine</td><td>1.16</td><td>1.06</td><td>Pos</td><td>1.48</td><td>1.76</td><td>I (9)</td></tr><tr><td>Spermidine</td><td></td><td></td><td>Pos</td><td>1.04</td><td>1.46</td><td>I</td></tr><tr><td colspan="7"><strong>Lipids</strong></td></tr><tr><td>Palmitoleic acid</td><td></td><td></td><td>Pos - Neg</td><td>1.77-1.52</td><td>0.36-0.41</td><td>I</td></tr><tr><td>Lyso PE (16:1)</td><td></td><td></td><td>Pos - Neg</td><td>2.07-1.71</td><td>0.44-0.50</td><td>II (3)</td></tr><tr><td>Lyso PC (16:1)</td><td></td><td></td><td>Pos</td><td>1.59 -</td><td>0.65</td><td>II (3)</td></tr><tr><td>Lyso PE (16:0)</td><td></td><td></td><td>Pos - Neg</td><td>1.97-1.72</td><td>0.41-0.36</td><td>I</td></tr><tr><td>Lyso PC (16:0)</td><td></td><td></td><td>Pos</td><td>1.73</td><td>0.41</td><td>II (3)</td></tr><tr><td>Lyso PE (18:1)</td><td></td><td></td><td>Pos - Neg</td><td>2.07-1.55</td><td>0.40-0.39</td><td>I</td></tr><tr><td colspan="7"><strong>Miscellaneous</strong></td></tr><tr><td>Acetate</td><td>1.01</td><td>0.84</td><td></td><td></td><td></td><td>(6)</td></tr><tr><td>Lactate</td><td>1.09</td><td>1.41</td><td></td><td></td><td></td><td>(9)</td></tr><tr><td>Succinate</td><td>1.05</td><td>1.40</td><td></td><td></td><td></td><td>(6)</td></tr><tr><td>Citrate</td><td>1.33-1.28</td><td>1.15-1.22</td><td></td><td></td><td></td><td>(3)</td></tr><tr><td>Beta-hydroxybutyrate</td><td>1.59-1.44</td><td>1.25-1.26</td><td></td><td></td><td></td><td>(5)</td></tr><tr><td>Pipecolic acid</td><td>1.96-1.59</td><td>1.46-1.25</td><td></td><td></td><td></td><td>(4)</td></tr><tr><td>Pantothenic acid</td><td>1.05-1.02</td><td>1.40-1.21</td><td>Pos - Neg</td><td>2.31-1.76</td><td>1.90-1.92</td><td>I (9)</td></tr><tr><td>Inosinic acid</td><td>1.73-1.41</td><td>0.71-0.72</td><td></td><td></td><td></td><td>(9)</td></tr><tr><td>Nicotinamide ribotide</td><td>1.19-1.03</td><td>0.86</td><td></td><td></td><td></td><td>(3)</td></tr><tr><td>AMP</td><td>1.73-1.38</td><td>0.72-0.76</td><td></td><td></td><td></td><td>(9)</td></tr><tr><td>ATP</td><td>2.09-2.17</td><td>0.43-0.40</td><td></td><td></td><td></td><td>(9)</td></tr><tr><td>NADH</td><td>1.37-1.38</td><td>0.75-0.76</td><td></td><td></td><td></td><td>(2)</td></tr><tr><td>UMP</td><td>2.04-1.92</td><td>0.43-0.45</td><td>Neg</td><td>2.03</td><td>0.47</td><td>I (9)</td></tr><tr><td>Glycerol</td><td>1.03</td><td>0.90</td><td></td><td></td><td></td><td>(9)</td></tr><tr><td>Glycerol 3-P</td><td></td><td></td><td>Neg</td><td>1.65</td><td>0.49</td><td>I</td></tr><tr><td colspan="7">Continued</td></tr></tbody></table>
<table><caption>**Table 1.** Metabolites identified according to LC-MS and NMR analyses. VIP scores indicate the relative influence of the corresponding metabolites on the discrimination between the two conditions (with and without H₂O₂) for both analytical techniques. Ratios S (in the presence of H₂O₂)/R (in the absence of H₂O₂) were determined from relative intensities in profiles at 50 min. Classes and scores are according to Sumner *et al.*<sup>59,60</sup>. *Under our analytical conditions, Trehalose/Cellobiose and Trehalose 6-P/Cellobiose 6-P were not differentiated.</caption><thead><tr><th rowspan="2">Identified VIP metabolite</th><th colspan="2">NMR</th><th colspan="3">LC-MS</th><th rowspan="2">ID metrics (class)</th></tr><tr><th>VIP Score</th><th>Ratio S/R</th><th>Ionization mode</th><th>VIP Score</th><th>Ratio S/R</th></tr></thead><tbody><tr><td>Betaine</td><td>2.13</td><td>4.46</td><td>Pos</td><td>2.70</td><td>5.14</td><td>I (2)</td></tr><tr><td>Glutathione</td><td>1.08</td><td>1.28</td><td></td><td></td><td></td><td>(9)</td></tr><tr><td colspan="7"><strong>Peptides</strong></td></tr><tr><td>Alanyl-alanine</td><td></td><td></td><td>Pos</td><td>2.01</td><td>0.65</td><td>I</td></tr><tr><td>Isoleucyl/leucyl-histidine</td><td></td><td></td><td>Pos</td><td>1.99</td><td>2.35</td><td>II (3)</td></tr><tr><td>Isoleucyl/leucyl-aspartic acid</td><td></td><td></td><td>Pos</td><td>2.58</td><td>6.14</td><td>II (3)</td></tr><tr><td>Valyl-isoleucine/leucine</td><td></td><td></td><td>Pos</td><td>1.74</td><td>4.89</td><td>II (3)</td></tr><tr><td>Valyl-serine</td><td></td><td></td><td>Pos</td><td>1.17</td><td>0.71</td><td>I</td></tr><tr><td>Valyl-aspartic acid</td><td></td><td></td><td>Pos</td><td>2.00</td><td>1.95</td><td>I</td></tr><tr><td>Gln-Asp-Thr-Pro</td><td></td><td></td><td>Pos</td><td>2.32</td><td>4.31</td><td>II (3)</td></tr><tr><td>Arg-Cys-Ser-Trp</td><td></td><td></td><td>Pos</td><td>1.85</td><td>0.57</td><td>II (3)</td></tr></tbody></table>

<figure><img src="image_5.png" /><figcaption>**Figure 4.** Carbohydrate metabolism showing metabolites over-produced in the presence of hydrogen peroxide (red boxes) and those under-produced (purple boxes). Likely metabolic intermediates are indicated in white boxes with dotted frame.</figcaption></figure>

stress. S/R ratios of valine and tryptophan were particularly high (4.5 and 6.98 respectively) compared to the others with a high VIP score (around 2.0) suggesting that these could be strong key markers of the response of *P. graminis* to H₂O₂ stress. Some compounds, which are connected to amino acid metabolism, were also over-produced, for instance carnitine and pipecolate linked to lysine metabolism or betaine linked to glycine metabolism. Betaine concentration was particularly increased (S/R = 4.46) and presented a high VIP value (2.13). This metabolite is also found in the case of cold shock<sup>43</sup>. However, some amino acid concentrations were decreased in the presence of H₂O₂, including alanine, glutamate, citrulline, proline and asparagine. Jozefczuk *et al.*<sup>30</sup> also observed an overproduction of most of the amino acids in *E. coli* exposed to H₂O₂, although some differences can
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Central metabolism: Citrate cycle, lipids, amino acids and peptides. The metabolites, which are over-produced in the presence of hydrogen peroxide, are indicated in red boxes, while those which are under-produced are indicated in purple boxes. Putative metabolic intermediates are indicated in white boxes with dotted frame.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Glutathione metabolism. The metabolites, which are over-produced in the presence of hydrogen peroxide are indicated in red boxes, while those, which are under-produced, are indicated in purple boxes. Putative metabolic intermediates are indicated in white boxes with dotted frame.</figcaption></figure>
be noticed concerning the exact nature of these amino acids (for instance alanine is overproduced in *E. coli* while it is underproduced in our case).

Some amino acids (glutamate, glycine) or amines linked to amino-acid pathways (carnitine, betaine) are well known osmo- or cryo-protectants. These compatible solutes accumulate in the cells facing osmotic or cold shocks. Carnitine is known as an antioxidant compound<sup>44</sup>. Except for glutamate whose concentration was decreasing, glycine, carnitine and betaine were accumulating in *P. graminis* 13b-3 in the presence of H<sub>2</sub>O<sub>2</sub>. This is consistent with studies that report a concomitant response to oxidative stress and cold shock<sup>45</sup>.

Interestingly, the modifications of amino acid concentrations were accompanied by the modulation of dipep-
tide and tetrapeptide concentrations (Table 1, Fig. 5). Ala-Ala, Val-Ser and [Arg-Cys-ser-tyrp] concentrations
decreased under stress conditions while the concentration of Ile(leu)-His, Ile(leu)-Asp, Val-Ile(leu), Val-Asp
and [Gln-Asp-Thr-Pro] increased. It is worth noting that some of these peptides had the highest VIP and S/R
values of all the metabolites of interest. It is particularly true for ILe(leu)-Asp (VIP score = 2.58, S/R = 6.14) and
[Gln-Asp-Thr-Pro] (VIP score = 2.32, S/R = 4.31). These peptides seem, thus, of major interest and might be the
more important key markers of H<sub>2</sub>O<sub>2</sub> stress for *Pseudomonas graminis* 13b-3. The discovery of dipeptides and
tetrapeptides was unexpected as they have never been described as a response to oxidative stress. Only a very few
previous studies report the presence of short peptides in bacterial cells. Kol et al.<sup>46</sup> reported an overproduction
of one dipeptide (Gly-Pro) and four tripeptides (Arg-Gly-Pro, Lys-Gly-Pro, Glu-Gly-Pro and Ala-Gly-Pro) in
response to a salt shock in *Steptomyces coelicolor* based on a metabolomics approach. These authors suggested
that these peptides all containing proline could be hydrolyzed by the cell to produce proline, which is a known
osmoprotectant. The overproduction of these peptides was concomitant with the increase of proline concentra-
tions. The presence of proline-containing short peptides was also described earlier in *Listeria monocytogenes*<sup>47</sup>,
*Lactobacillus casei*<sup>48</sup> *Oenococcus oeni*<sup>49</sup> and *Bacillus subtilis*<sup>50</sup> facing osmotic shocks. In our case, only one tetra-
peptide [Gln-Asp-Thr-Pro] contained proline and - although it was overproduced in the presence of H<sub>2</sub>O<sub>2</sub> - pro-
line concentration was not increased but decreased. The mechanism involved there seems thus quite different,
and the exact role of these dipeptides and tetrapeptides remains unknown. Using a metabolomics approach, Jousse
et al.<sup>43</sup> observed the presence of short peptides in *Pseudomonas syringae* 32b-74 under cold shock. In this latter
case, these peptides were not specifically containing proline and they were not all overproduced. In parallel,
the free amino-acid concentrations were mainly decreased under cold conditions. It was suggested that short
peptides could play unknown roles in regulating pathways. In the present study, the modulation of dipeptide
and tetrapeptide concentrations could be also related to regulating pathways. In addition, when peptides were
accumulating, we cannot exclude their potential role in the protection against oxidative damages. For instance, it
was shown earlier that linear glycine-containing dipeptides could scavenge radicals in human erythrocytes when
exposed to oxidative stress<sup>51</sup>.

As expected during oxidative stress, glutathione GSH concentration increased in *P. graminis* (Fig. 6, Table 1). However, we can note that glutathione was not among the highest VIP metabolites (1.08 score) and the S/R ratio was not very high (1.28). Glutathione is widely found in many living organisms including microorganisms; it has various functions involved in the defense against oxidative stress<sup>11,44</sup>. Consistently with the increase of glutathione concentration, other metabolites directly linked to the same pathway, such as glycine, methionine, putrescine, arginine, spermidine, beta-alanine, UMP and pantothenate, were also found increased in the presence of oxidants in our study (Fig. 6).

Putrescine and spermidine are among the most widely distributed cellular polyamines and these are essential for normal cellular growth and multiplication of both prokaryotic and eukaryotic cells. Spermine and spermidine can act as free radical scavengers and reduce cellular damages induced by oxygen radicals<sup>52</sup>. Polyamines therefore play an important adjunctive role in protecting cells from the toxic effects of ROS.

This work shows of the plasticity of *P. graminis* 13b-3 metabolome in clouds. *Pseudomonas graminis* responded to oxidative stress by modifying many metabolic pathways (carbohydrates, lipids, amino-acids, amines, glutathione, peptides). However, *P. graminis* was short after able to recover, as attested by increasing ATP concentration (Fig. 1) and using similar metabolic pathways as the reference cells (Fig. 2) when the medium had been completely detoxified. These results clearly demonstrate the flexibility of *P. graminis* after H₂O₂ exposure. This is consistent with the fact that this bacterium and others isolated from clouds survive when exposed to H₂O₂ at high concentrations and were even able to grow after H₂O₂ exposure¹⁰. Previous experiments also showed a very strong correlation between ATP and cloud H₂O₂ concentrations measured directly in cloud water¹¹. This suggests a great plasticity of the cloud microorganism metabolome which constantly adapts to the H₂O₂ concentrations. This aspect is particularly important in the atmospheric context as H₂O₂ concentration is highly variable in clouds depending on various parameters (e.g. photolytic activity, temperature and concentrations of trace gases, such as VOCs, CO, O₃, and NOₓ). Sunlight is particularly important; as a result, diurnal and seasonal variations of H₂O₂ concentration are observed. In the long term, the atmospheric concentration of H₂O₂ might increase due to the increasing concentration of pollutants such as NOₓ².

The very high oxidative condition present in clouds (up to 200 mV) is probably one of the major constraints experienced by cloud microorganisms. Consistently, a recent metatrancriptomics analysis of the cloud microbiome highlighted in particular the expression of genes involved in oxidative stress response, antioxidant activity, as well as glutathione biosynthesis in bacteria⁵³.

Finally our work points out the interactions between the bacterial metabolome and cloud chemistry.

Hydrogen peroxide is a key component in the chemical reactivity in the atmosphere, driving most of the radical reactions as a major source of hydroxyl radicals (°OH)² (Fig. 7). In the water phase of clouds, H₂O₂ can interact with iron and/or UV light to produce °OH radicals by Fenton, photo-Fenton or direct photolysis reactions, even though the major part of °OH radical is likely transferred from the gas phase to the water phase of clouds. OH is the main oxidant for organics in cloud water³. This OM is very complex, and its molecular composition remains still unknown. Two major classes of small organic compounds have been identified and quantified, namely
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Simplified scheme of the interactions of H₂O₂ with the bacterium metabolome and the organic matter (OM) present in the cloud droplet. The modulation of the metabolome by H₂O₂, in particular, amino acid and carboxylic acid pathways, could have an impact on cloud chemistry by changing the amino acid and carboxylic acid concentrations in the cloud droplet. Carbox acid: carboxylic acids. Dotted lines indicate unknown pathways (transport, reactions...).</figcaption></figure>

carboxylic acids and aldehydes (10%) and amino acids (9%)¹,⁴,⁵. These compounds are also major components of fog and aerosols⁴,⁷ and contribute to atmospheric chemistry and the global carbon and nitrogen cycles. The origins of carboxylic acids are direct emission, transformation in the gas phase and transfer to the cloud droplet, or they can also result from the transformation of the OM by radical chemistry directly in cloud water⁴. The origin of amino acids is not clearly established as they have been discovered only recently in clouds⁵. Free amino acids, peptides and proteins have been found in aerosols and are likely issued from biological aerosols particles (BAPs)⁷. The presence of "peptide-like compounds" has been shown in cloud water by global FT-ICR-MS analyses although their concentration has not been quantified⁸,⁹. Both amino acids and carboxylic acids can interact with °OH radicals in the water phase of clouds for further reactions³⁻⁵.

It was shown that microbial metabolism could be an alternative route to radical chemistry in the transformation pathways of carbon compounds and oxidants in clouds¹⁵,²⁷. Ariya *et al.*²⁰ showed that dicarboxylic acids (C₂–C₉), namely oxalic, malonic, succinic, glutaric, adipic, pimelic, and pinic acids could be degraded efficiently by airborne-microorganisms. Using microcosms mimicking the cloud environment or real cloud samples, Vaitilingom *et al.*¹⁹,²²,²³ and Husarova *et al.*²⁴ showed that the biodegradation rates of organic acids (acetate, succinate, malonate) and of C1 compounds (formate, formaldehyde, methanol) were within the same concentration range as photodegradation rates. In addition, microorganisms isolated from cloud water have been shown to produce organic compounds, such as pyruvate from lactic acid²¹. The biotransformation of amino acids present in clouds has not been published yet, but it is very likely to occur as it was shown that microorganisms can grow using the carbon and nitrogen sources present in cloud water, thus building proteins from amino acids²¹. In addition, the presence of transcripts of genes involved in amino acid pathways has been shown recently in clouds⁵³. We have also confirmed here that H₂O₂ can be transformed by catalases, consequently microorganism activity could decrease the efficiency of radical chemistry because it reduces the source of radicals¹¹,¹⁹.

In this work, we have clearly shown that H₂O₂ directly impacts the metabolome of a cloud bacterium (Fig. 7). Many pathways were modified including glutathione, saccharide, lipids, energy (ATP, NADH), peptides, amino acids and carboxylic acids pathways. The metabolites of these pathways could be taken up or excreted by the bacterium in the cloud droplet, and thus potentially change the cloud OM composition. This is particularly true for amino acids and carboxylic acids, which are major components of the cloud OM (Fig. 7). Lactate, succinate and citrate were overproduced while acetate was underproduced in the presence of H₂O₂. Among the amino acids whose concentrations are modulated under oxidative stress, valine and tryptophan concentrations were drastically increased (around 4.5- to 7-fold). Peptides are also key markers of H₂O₂ stress and could change the "peptide-like" concentrations in cloud water. We have shown that *P. graminis* has a great metabolome plasticity and the capacity of acclimate to a range of H₂O₂ concentrations. This is quite important regarding atmospheric
context as H₂O₂ concentrations are highly variable. This means that not only cloud microorganisms can compete with radical processes but this competition could be modulated by H₂O₂ concentration, i.e. by atmospheric scenarios. We are not able to quantify this modulation at the moment; future work is needed, notably looking at amino acids and peptide transformations. Although this work was based on a single strain, we can expect that this phenomenon is true for all cloud metabolomes as we previously showed a very strong correlation between ATP and H₂O₂ concentration measured in real clouds, ATP being a proxy of the whole cell metabolism¹¹.

## Material and Methods

The general workflow describing the experiment protocol used for the metabolomics experiments is presented in Fig. S5.

Three independent experiments were performed (three batches of twelve samples each). Furthermore, the biodegradation of hydrogen peroxide as well as the evolution of ATP concentration were followed for each experiment over time. These experiments were conducted in marine artificial cloud solution as the majority of clouds sampled at the puy de Dôme station are from marine origin²³. The chosen incubation temperature was 17 °C, because it is the average temperature in summer at the puy de Dôme station (1465 m).

In nature, microorganisms are aerosolized and then condense water to become part of the cloud water system, therefore they are directly transferred from their original rich environment (soil, vegetation, surface waters) to the cloud medium containing H₂O₂ and have no time to grow first in this new medium before being exposed to H₂O₂. This is the reason why we grew the cells in R2A, washed them twice with the cloud medium to eliminate the nutrients and salts (we checked by NMR by there was no nutrient left originating from R2A such as glucose), and transferred them directly in cloud water with or without H₂O₂ to really see the impact of H₂O₂ on their metabolism.

**Chemical reagents.** Hydrogen peroxide solution was freshly prepared just before each experiment from a commercial solution (H₂O₂; 30%; not stabilized; Fluka Analytical). Acetonitrile (CH₃CN; 99.99%; Optima LC-MS; Fisher Scientific), methanol (CH₃OH; 99.9%; Fluka Analytical) and water (H₂O; Optima LC-MS; Fisher Scientific) were used in the proportions 2:2:1, respectively, for metabolites' extraction. Water (H₂O; Biosolve) and acetonitrile (CH₃CN; 99.97%) in the proportions 1:1 (v/v), acidified at 0.1% with a formic acid solution (HCOOH; 99%; Biosolve), were used for LC mobile phase. Deuterium oxide (D₂O; 99.96%D; Eurisotop), phosphate buffer (KH₂PO₄; 99.50% & K₂HPO₄; 99.50%, at a final concentration of 30 mM; Sigma Aldrich) and sodium azide (NaN₃; 99.99%; Sigma Aldrich) were used for NMR sample preparation. Carbon source, nitrogen compounds and salts used for artificial cloud solution are commercially available (Table S4).

Standard solutions used for LC/MS validation which were colchicine, tryptophan, phenylalanine, amino anthracene and creatinine were prepared from stock solution (0.5 g/L for each compound in water/acetonitrile 1:1 (v/v) solution containing 0.1% formic acid and stored at -20 °C) diluted in the same solvent at a final concentration of 5.0 10⁻³ g/L before use.

**Basic physiological characteristics of *Pseudomonas graminis* PDD-13b-3.** The bacterial strain PDD-13b-3 (GenBank DQ512786) was isolated from cloud water collected on the 7th of August 2004 by impaction from puy de Dôme summit (1465 m asl; France), by culture on R2A medium at 17 °C; it was identified as a species of *Pseudomonas graminis*¹⁷. Its draft genome was sequenced recently and is available under the Bioproject accession number PRJNA362581³². In the laboratory, this bacterium is able to grow at temperatures at 5 °C, 17 °C, 27 °C and 37 °C (lower and higher temperatures were not investigated) in R2 medium, with generation times of 8.6, 4.45 and 3.7 hours at 5 °C, 17 °C and 37 °C respectively, and an optimum for growth at 27 °C with 1.2 h/generation²⁸. Its cellular ATP content measured during exponential growth in R2 at 17 °C is 5.5 (± 0.7) 10⁻¹⁸ mol cell⁻¹. Previous work showed that this strain is able to uptake and use as unique C sources acetate, formate, L-lactate, succinate, formaldehyde, and to a lesser extent also methanol²¹.

**Microbial cultivation and incubations in microcosms.** For this work, bacteria were grown aerobically in 100 mL of R2A medium under agitation (200 rpm) at 17 °C for 17 hours. Then, cells in the exponential growth phase were collected by centrifugation for 3 min at 10,500 g. The supernatant was removed, and the bacterial pellet was suspended and washed twice in artificial cloud solution²³. The concentration of cells was estimated by optical density measurement at 575 nm and diluted to obtain a concentration of ~10⁷ cell.mL⁻¹; biomass was further precisely quantified by flow cytometry (BD Facscalibur Becton-Dickinson; λexc = 488 nm; λem = 530 nm) on SYBR green stained samples. The artificial cloud solution was used to mimic cloud conditions; it contains the major carbon and nitrogen compounds and salts present in clouds originating from the ocean when samples are collected at the sampling site (Table S4). H₂O₂ concentration was fixed at 200 µM for our experiments. The cell concentration/chemical concentration ratio (Carbon compounds, ions and oxidants) was kept identical to that of a real cloud. We showed that biodegradation rates are kept constant when the cell/substrate ratios are constant²².

Amber Erlenmeyer's flasks containing 150 mL of cloud water solution were inoculated with 10⁷ cells.mL⁻¹ and incubated at 17 °C, 110 rpm, and half were added with 200 µM of hydrogen peroxide. Two Erlenmeyer flasks (with and without H₂O₂) were used to monitor the evolution of hydrogen peroxide and ATP concentrations over time. The experiment was repeated three times, resulting in three biological replicates. Hydrogen peroxide concentration was measured by fluorescence spectroscopy (Safire II TECAN®; λexc = 320 nm, λem = 390 nm) as in Wirgot et al.¹¹. This method is based on a reaction between hydrogen peroxide, horse radish peroxidase and p-hydroxyphenylacetic acid. Bioluminescence was used to determine ATP concentration (Glomax® 20/20 single tube luminometer from Promega); this is based on the endergonic oxidation of D-luciferin by the enzyme luciferase, and producing an amount of photons directly proportional to that of ATP in the sample¹¹. The protocol used was adapted from ATP Biomass Kit HS by Biothema (Sweden).
Twelve Erlenmeyer flasks were used to study the modulation of the bacterial metabolome in the presence or absence of H₂O₂. After 50 minutes and 24 hours of incubation (i.e., time after H₂O₂ addition), respectively, sub-samples were collected, cell pellets were obtained by centrifugation (12000 g, 4 min, 4 °C) and rinsed twice in 0.8% sodium chloride solution (w/v).

Centrifugation followed by a quick quenching with cold solvent was chosen instead of filtration because it was more adapted to the high volumes of samples (12 samples of 150 mL). The advantage of centrifugation is that all the samples are centrifuged simultaneously. A simultaneous treatment is important because the incubation is stopped at 50 min, therefore any additional minutes between the samples due to extraction/quenching procedure will lead to very heterogeneous results.

**Samples preparation for metabolomics.** The metabolomes were extracted from pellets using 1.2 mL of water/methanol/acetonitrile [1:2:2] solution kept in an ice bath<sup>43</sup>. After a last step of centrifugation intended to remove cell debris (5 min, 4 °C, 12000 g), 750 µL of the supernatants were dedicated to NMR and four aliquots of 50 µL were dedicated to LC-MS analyses. Moreover, five aliquots of 10 µL of each sample were mixed together for serving as quality controls through LC-MS analyses.

For NMR analyses, the supernatants were evaporated at 10 °C for approximately 48 hours using a refrigerated centrivap concentrator equipped with a centrivap cold trap at -50 °C (LabConco; ThermoScientific). The extracts were suspended into 300 µL of deuterium oxide (D₂O, 99.96%D) and dried two times in order to eliminate H₂O traces. Finally, the extracts were suspended in 90% of phosphate buffer (30 mM) prepared in D₂O (99.96%D).

**Metabolic profiling.** *LC-MS analyses for metabolite profiling.* LC-MS analyses were performed on the Metabolic profiler® platform (Bruker) using a fast LC system (Agilent 1200 series) coupled to a MicroToF mass spectrometer equipped with electrospray ionization (ESI) (Bruker). Chromatography separation was performed on an Acquity HSST3 C18 column (150\*2.1 mm 1.8 µm; Waters). The flow rate was fixed at 0.2 mL/min with 1 L of water and acetonitrile containing 0.1% of formic acid for mobile phases A and B, respectively. The gradient elution was carried out as follows: 0–2 min, 100% A; 2–15 min linear gradient to 0% A; 15–22 min, 0% A; 22–22.10 min linear gradient back to 100% A; 4 min equilibration wash with 100% A. The injection volumes for all samples were 6 µL and the column temperature was set to 30 °C. The mass spectrometer operated with full scan (50–1000 m/z) in positive and negative ion modes. Nitrogen was used as the drying gas. The gas nebulizer (nitrogen) pressure was 40 psi, the desolvation gas flow rate was 9 L/min and temperature was maintained at 200 °C. Capillary tension was 4500 V. Finally, softwares used for data acquisition were OtofControl 3.2 and Hystar 3.2 (Bruker Daltonics).

Matrix of relative intensities for all samples (Stressed and Reference) corresponding to each identified metabolite obtained for LC-MS analyses is shown in Table S5. Means, standard deviations and ratios S/R were calculated.

*1D NMR analyses for metabolite profiling.* One-dimensional ¹H NMR experiments were done with an Avance III 500 MHz NMR spectrometer using a 5 mm Prodigy TCI z-gradient (¹H/¹⁹F/¹³C/¹⁵N) probe (Bruker Biospin Wissenbourg, France). Samples were analyzed on 5 mm tubes refrigerated at 6 °C on an auto-sampler (Sample Jet). NMR spectra were recorded at 25 °C K. A standard one-dimensional *noe* spectroscopy sequence (*noesy*gppr1d with water presaturation and gradients) was used with low power irradiation of the water resonance during the recycle delay of 2 s and the mixing time of 10 ms. A total of 1024 scans was collected with a 90° impulsion time of 8.70 µs, acquisition time of 3.3 s, spectral window of 64 K data points zero-filled to 128 K before Fourier transformation with 0.3 Hz line broadening. All spectra were finally processed with Topspin version 3.5pl5 (Bruker Daltonics).

**Data treatments.** *LC/MS preprocessing of chromatograms.* MS raw data were processed using the web-based platform Galaxy framework including preprocessing, normalization and quality control steps<sup>54</sup>. XCMS software implanted in the Galaxy platform was used to process raw data for feature detection, alignment and framing<sup>55</sup>. This step produced a table of features characterized by sample retention time, mass to charge (m/z) ratio and intensity (i.e., peak area). XCMS parameters used for data processing are shown in Table S6. After we observed an analytical deviation of signals, intensities of peaks were corrected with a regression model in order to remove signal drift and batch-effects<sup>56</sup>. The stability of signal intensities across batches was then evaluated by relative standard deviation (RSD) of each feature in QCs (quality control standard). Features with RSD in QCs higher than in samples were excluded. We kept only features with RSD in QCs below 30% for further multivariate analysis.

*NMR preprocessing of spectra.* 1D ¹H spectra were divided into regions (buckets) of 0.01 ppm width. Over the chemical range of 0.5 (right border) to 10 (left border) ppm, values were not used in the bucketing. Furthermore, the area around residual water signal (4.6 to 5.1 ppm) was excluded. Bucketing was performed using the AMIX software (Bruker GmbH). The signal intensity in each bucket was normalized by the total intensities of spectra. Finally, a data matrix of relative intensities characterized by chemical shifts was obtained and used for multivariate analysis.

Matrix of relative intensities for all samples (Stressed and Reference) corresponding to each identified metabolite obtained for NMR analyses is shown in Table S7. Means, standard deviations and ratios S/R were calculated.

**Statistical analysis.** The data matrices obtained for LC-MS (positive and negative mode) and NMR (¹H) were subjected to multivariate and univariate analyses. First, multivariate analyses were performed using SIMCA-P software (v12, Umetrics) to identify a potential discrimination between the samples incubated in the presence or absence of H₂O₂ and for the two times of incubation. Principal Component Analysis (PCA) (Figs S1 and S2) and Partial Least Square Discriminant Analysis (PLS-DA, Figs 2 and 3) were done. A selection
of discriminant ions and buckets was done based on variable importance in projection (VIP) values. The values strictly greater than unity were kept for the identification step. In addition to multivariate analyses, univariate analysis was done (Wilcoxon-Mann-Whitney test); ions and buckets with a p-value lower than 0.05 were kept for the identification step.

**Identification of discriminants variables.** UPLC-MS/MS analyses for structural identification of the metabolites. CAMERA R package was used to annotate adduct peaks, isotopes and in-source fragments of compound spectra in positive and negative mode⁵⁷. QCs were reinjected in order to identify the most discriminant features in multivariate analysis. Each feature was identified from in-source fragments and adducts. Additional structural information was obtained by performing MS/MS fragmentation using UPLC-qTOF-MS/MS. MS/MS fragmentation patterns were compared to simulations through Massfrontier (v.7; Thermo Scientific) software. When standards were available, these were analysed to confirm putative identifications.

1D and 2D NMR analyses for structural identification of the metabolites. For better sensitivity and resolution of spectra, a NMR spectrometer Bruker Avance III 950 MHz equipped with a 5 mm TCI (¹H/¹³C/¹⁵N/²H) cryoprobe with z-gradient coil probe (Bruker Biospin Wissenbourg, France) was used, in order to identify discriminant buckets. For 1D ¹H-spectra, a standard one-dimensional spectroscopy sequence (noesygpr1d) was used with low power irradiation of the water resonance during the recycle delay of 4 s and the mixing time of 10 ms. For each spectrum, eight scans were collected with an 90° impulsion time of 7.7 µs, an acquisition time of 3.3 s, a spectral window of 10 000 Hz and 64 K data points zero-filled to 128 K before Fourier transformation with 0.3 Hz line broadening. All 2D homonuclear (¹H-¹H COSY, ¹H-¹H TOCSY, ¹H-¹H JRES) and heteronuclear (¹H-¹³C HSQC and HMBC) experiments were performed with quadrature phase detection in dimensions, using state-TPPI or QF detection mode in the indirect one. For each of the 512 increments in the indirect dimension, 2 K data points were collected and 32 or 64 transients were accumulated in the direct dimension. A ¹³C decoupling (GARP) was performed during acquisition time for heteronuclear experiments. A π/2 shifted square sine-bell function was applied in both dimensions before Fourier transformation. All NMR spectra were recorded at 300 K. Spectra were processed (phase and baseline correction) with Topspin version 3.5pl5.

Putative identifications for LC/MS and NMR discriminant variables were performed using web data base searches such as KEGG³³ and HMDB⁵⁸.

**Data Availability**

NMR and MS data are available via https://doi.org/10.15454/Q7VKT5 and https://doi.org/10.15454/ZAWDKQ.

**References**

1. Deguillaume, L. et al. Classification of clouds sampled at the puy de Dôme (France) based on 10 yr of monitoring of their physicochemical properties. *Atmos. Chem. Phys.* **14**, 1485–1506, https://doi.org/10.5194/acp-14-1485-2014 (2014).

2. Vione, D., Maurino, V., Minero, C. & Pelizzetti, E. The atmospheric chemistry of hydrogen peroxide: a review. *Ann. Chim-Rome* **93**, 477–488 (2003).

3. Herrmann, H. et al. Tropospheric aqueous-phase chemistry: kinetics, mechanisms, and its coupling to a changing gas phase. *Chem. Rev.* **115**, 4259–334, https://doi.org/10.1021/cr500447k (2015).

4. Herckes, P., Valsaraj, K. T. & Collett, J. L. A review of observations of organic matter in fogs and clouds: Origin, processing and fate. *Atmos. Res.* **132–133**, 434–449, https://doi.org/10.1016/j.atmosres.2013.06.005 (2013).

5. Bianco, A., Voyard, G., Deguillaume, L., Mailhot, G. & Brigante, M. Improving the characterization of dissolved organic carbon in cloud water: Amino acids and their impact on the oxidant capacity. *Sci. Rep.* **6**, 37420, https://doi.org/10.1038/srep37420 (2016).

6. Zhang, Q. & Anastasio, C. Free and combined amino compounds in atmospheric fine particles (PM2.5) and fog waters from Northern California. *Atmos. Environ.* **37**, 2247–2258, https://doi.org/10.1016/S1352-2310(03)00127-4 (2003).

7. Matos, J. T. V., Duarte, R. M. B. O. & Duarte, A. C. challenges in the identification and characterization of free amino acids and proteinaceous compounds in atmospheric aerosols: A critical review. *Trends in Anal. Chem.* **75**, 97–107, https://doi.org/10.1016/j.trac.2015.08.0042016 (2016).

8. Zhao, Y., Hallar, A. G. & Mazzoleni, L. R. Atmospheric organic matter in clouds: exact masses and molecular formula identification using ultrahigh-resolution FT-ICR mass spectrometry. *Atmos. Chem. Phys.* **13**, 12343–12362, https://doi.org/10.5194/acp-13-12343-2013 (2013).

9. Bianco, A. et al. Molecular Characterization of Cloud Water Samples Collected at the puy de Dôme (France) by Fourier Transform Ion Cyclotron Resonance Mass Spectrometry. *Environ. Sci. Technol.*, https://doi.org/10.1021/acs.est.8b01964 (2018).

10. Joly, M. et al. Survival of microbial isolates from clouds toward simulated atmospheric stress factors. *Atmos. Environ.*, https://doi.org/10.1016/j.atmosenv.2015.07.009 (2015).

11. Wirgot, N., Vinatier, V., Deguillaume, L., Sancelme, M. & Delort, A.-M. H₂O₂ modulates the energetic metabolism of the cloud microbiome. *Atmos. Chem. Phys.* **17**, 14841–14851, https://doi.org/10.5194/acp-17-14841-2017 (2017).

12. Sattler, B., Puxbaum, H. & Psenner, R. Bacterial growth in supercooled cloud droplets. *Geophys. Res. Lett.* **28**, 239–242, https://doi.org/10.1029/2000GL011684 (2001).

13. Amato, P. et al. Active microorganisms thrive among extremely diverse communities in cloud water. *Plos One* **12**(8), e0182869 (2017).

14. Hill, K. A. et al. Processing of atmospheric nitrogen by clouds above a forest environment. *J. Geophys. Res.* **112**, https://doi.org/10.1029/2006JD008002 (2007).

15. Hu, W. et al. Bacteria in atmospheric waters: Detection, characteristics and implications. *Atmos. Environ.* **179**, 201–221, https://doi.org/10.1016/j.atmosenv.2018.02.026 (2018).

16. Amato, P. et al. Microbial population in cloud water at the Puy de Dôme: implications for the chemistry of clouds. *Atmos. Environ.* **39**, 4143–4153, https://doi.org/10.1016/j.atmosenv.2005.04.002, (2005).

17. Amato, P. et al. Microorganisms isolated from the water phase of tropospheric clouds at the Puy de Dôme: major groups and growth abilities at low temperatures. *FEMS Microbiol. Ecol.* **59**, 242–254, https://doi.org/10.1111/j.1574-6941.2006.00199 (2007).

18. Vaitilingom, M. et al. Long-term features of cloud microbiology at the puy de Dôme (France). *Atmos. Environ.* **56**, 88–100, https://doi.org/10.1016/j.atmosenv.2012.03.072 (2012).

19. Vaitilingom, M. et al. Potential impact of microbial activity on the oxidant capacity and organic carbon budget in clouds. *Proc. Natl. Acad. Sci USA* **110**, 559–564, https://doi.org/10.1073/pnas.1205743110 (2013).
20. Ariya, P. A., Nepotchatkykh, O., Ignatova, O. & Amyot, M. Microbiological degradation of atmospheric organic compounds. *Geophys. Res. Lett.* **29**, 34-1-34-4, https://doi.org/10.1029/2002GL015637 (2002).

21. Amato, P. et al. A fate for organic acids, formaldehyde and methanol in cloud water: their biotransformation by micro-organisms. *Atmos. Chem. Phys.* **7**, 5253-5276, https://doi.org/10.5194/acp-7-4159-2007 (2007).

22. Vaïtilingom, M. et al. Contribution of microbial activity to carbon chemistry in clouds. *Appl. Environ. Microb.* **76**, 23-29, https://doi.org/10.1128/AEM.01127-09 (2010).

23. Vaïtilingom, M. et al. Atmospheric chemistry of carboxylic acids: microbial implication versus photochemistry. *Atmos. Chem. Phys.* **11**, 8721-8733, https://doi.org/10.5194/acp-11-8721-2011 (2011).

24. Husárová, S. et al. Biotransformation of methanol and formaldehyde by bacteria isolated from clouds. Comparison with radical chemistry. *Atmos. Environ.* **45**, 6093-6102, https://doi.org/10.1016/j.atmosenv.2011.06.035 (2011).

25. Krumins, V., Mainelis, G., Kerkhof, L. J. & Fennell, D. E. Substrate-Dependent rRNA Production in an Airborne Bacterium. *Environ. Sci. Technol.* **1**, 376-381, https://doi.org/10.1021/ez500245y (2014).

26. Ervens, B. Modeling the Processing of Aerosol and Trace Gases in Clouds and Fogs. *Chem. Rev.* **115**(10), 4157-4198, https://doi.org/10.1021/cr5005887 (2015).

27. Delort, A.-M et al. Impacts on cloud chemistry. In "Microbiology of aerosols", Wiley, Editors: A.-M. Delort, P. Amato, Chapter 3.2, 221-248 (2017).

28. Sies, H. Oxidative stress: introduction. In *Oxidative Stress: Oxidants and antioxidants*. ed. SIES. H. Academic Press. London, 15-22 (1991).

29. Noctor, G., Lelarge-Trouverie, C. & Mhamdi, A. The metabolomics of oxidative stress. *Phytochemistry.* **112**, 33-53, https://doi.org/10.1016/j.phytochem.2014.09.002 (2015).

30. Jozefczuk, S. et al. Metabolomic and transcriptomic stress response of *Escherichia coli*. *Mol. Syst. Biol.* **6**, 364, https://doi.org/10.1038/msb.2010.18 (2010).

31. Mailloux, R. J., Lemaire, J. & Appanna, V. D. Metabolic networks to combat oxidative stress in *Pseudomonas fluorescens*. *Antonie Van Leeuwenhoek.* **99**(3), 433-442, https://doi.org/10.1007/s10482-010-9538-x (2011).

32. Besaury, L., Amato, P., Wirgot, N., Sancelme, M. & Delort, A.-M. Draft genome sequence of *Pseudomonas graminis* PDD-13b-3, a model strain isolated from cloud water. *Genome Announcement.* **5**, e00464-17 (2017).

33. Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M. & Tanabe, M. KEGG as a reference resource for gene and protein annotation. *Nucleic Acids Res.* **44**(D1), D457-462, https://doi.org/10.1093/nar/gkv1070 (2016).

34. Couée, I., Sulmon, C., Gouesbet, G. & El Amrani, A. Involvement of soluble sugars in reactive oxygen species balance and responses to oxidative stress in plants. *J. Exp. Bot.* **57**(3), 449-459, https://doi.org/10.1093/jxb/erj027 (2006).

35. Demple, B. Regulation of Bacterial Oxidative Stress Genes. *Annu. Rev. Genet.* **25**, 315-337, https://doi.org/10.1146/annurev.ge.25.120191.001531 (1991).

36. Farr, S. B. & Kogoma, T. Oxidative stress responses in *Escherichia coli* and *Salmonella typhimurium*. *Microbiol. Rev.* **55**, 561-585 (1991).

37. Perricone, C. D., Park, S., Imlay, J. A. & Weiser, J. N. Factors contributing to hydrogen peroxide resistance in *Streptococcus pneumoniae* include pyruvate oxidase (SpxB) and avoidance of the toxic effects of the Fenton reaction. *J. Bacteriol.* **185**, 6815-6825, https://doi.org/10.1128/JB.185.23.6815-6825.2003 (2003).

38. Spragg, R. G., Hinshaw, D. B., Hyslop, P. A., Schraufstatter, I. U. & Cochrane, C. G. Alterations in adenosine triphosphate and energy charge in cultured endothelial and P388D1 cells after oxidant injury. *J. Clin. Invest.* **76**, 1471-1476, https://doi.org/10.1172/JCI112126 (1985).

39. Hyslop, P. A. et al. Mechanisms of oxidant-mediated cell injury the glycolytic and mitochondrial pathways of ADP phosphorylation are major intracellular targets inactivated by hydrogen peroxide. *J. Biol. Chem.* **263**, 1665-1675 (1988).

40. Tiwari, B. S., Belenghi, B. & Levine, A. Oxidative stress increased respiration and generation of reactive oxygen species, resulting in ATP depletion, opening of mitochondrial permeability transition, and programmed cell death. *Plant Physiol.* **128**, 1271-1281, https://doi.org/10.1104/pp.010999 (2002).

41. Oka, S. I., Hsu, C. P. & Sadoshima, J. Regulation of cell survival and death by pyridine nucleotides. *Circ. Res.* **111**, 611-627, https://doi.org/10.1161/CIRCRESAHA.111.247932 (2012).

42. Singh, R., Mailloux, R. J., Puiseux-Dao, S. & Appanna, V. D. Oxidative stress evokes a metabolic adaptation that favors increased NADPH synthesis and decreased NADH production in *Pseudomonas fluorescens*. *J. Bacteriol.* **189**, 6665-6675, https://doi.org/10.1128/JB.00555-07 (2007).

43. Jousse, C. et al. Metabolomics study of the response to cold shock in a strain of *Pseudomonas syringae* isolated from cloud water. *Metabolomics*, https://doi.org/10.1007/s11306-017-1295-7 (2017).

44. Davies, K. J. A. Oxidative stress antioxidant defenses and damage removal repair and replacement systems. *IUBMB Life.* **50**, 279-289, https://doi.org/10.1080/713803728 (2000).

45. Zhang, J., Li, Y., Chen, W., Du, G. C. & Chen, J. Glutathione improves the cold resistance of *Lacotococcus sanfranciscensis* by physiological regulation. *Food microbiol.* **31**, 285-292, https://doi.org/10.1016/j.fm.2012.04.006 (2012).

46. Kol, S. et al. Metabolomic Characterization of the Salt Stress Response in *Streptomyces coelicolor*. *Appl. Environ. Microbiol.* **76**, 2574-2581, https://doi.org/10.1128/AEM.01992-09 (2010).

47. Amezaga, M.-R. et al. The role of peptide metabolism in the growth of *Listeria monocytogenes* ATCC 23074 at high osmolarity. *Microbiology.* **141**, 41-49, https://doi.org/10.1099/00221287-141-1-41 (1995).

48. Piuri, M., Sanchez-Rivas, C. & Ruzal, S. M. Adaptation to high salt in *Lactobacillus*: role of peptides and proteolytic enzymes. *J. Appl. Microbiol.* **95**, 372-379, https://doi.org/10.1046/j.1365-2672.2003.01971.x (2003).

49. Le Marrec, C., Bon, E. & Lonvaud-Funel, A. Tolerance to high osmolality of the lactic acid bacterium *Oenococcus oeni* and identification of potential osmoprotectants. *Int. J. Food Microbiol.* **115**, 335-342, https://doi.org/10.1016/j.ijfoodmicro.2006.12.039 (2007).

50. Zaprasis, A. et al. Osmoprotection of *Bacillus subtilis* through import and proteolysis of pro-line-containing peptides. *Appl. Environ. Microbiol.* https://doi.org/10.1128/AEM.01934-12 (2013).

51. Zheng, L., Dong, H., Su, G., Zhao, Q. & Zhao, M. Radical scavenging activities of Tyr-Trp-Cys and Met-Gly and their protective effects against AAPH-induced oxidative damage in human erythrocytes. *Food Chem.* **197**, 807-813, https://doi.org/10.1016/j.foodchem.2015.11.012 (2016).

52. Ha, H. C. et al. The natural polyamine spermine functions directly as a free radical scavenger. *Proc. Natl. Acad. Sci.* **95**(19), 11140-11145, https://doi.org/10.1073/pnas.95.19.11140 (1998).

53. Amato, P. et al. Metatranscriptomics exploration of microbial functioning in clouds. *Sci. Rep.* **9**, 4383, https://doi.org/10.1038/s41598-019-41032-4.2019 (2019).

54. Giacomoni, F. et al. Workflow4Metabolomics: a collaborative research infrastructure for computational metabolomics. *Bioinformatics.* **31**, 1493-1495, https://doi.org/10.1093/bioinformatics/btu813 (2015).

55. Smith, C. A., Want, E. J., O'Maille, G., Abagyan, R. & Siuzdak, G. XCMS: Processing Mass Spectrometry Data for Metabolite Profiling Using Nonlinear Peak Alignment, Matching, and Identification. *Anal. Chem.* **78**, 779-787, https://doi.org/10.1021/ac051437y (2006).
56. Dunn, W. B. et al. The Human Serum Metabolome (HUSERMET) Consortium: Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry. *Nature Protocols*. **6**, 1060–1083, https://doi.org/10.1038/nprot.2011.335 (2011).

57. Kuhl, C., Tautenhahn, R., Böttcher, C., Larson, T. R. & Neumann, S. CAMERA: An Integrated Strategy for Compound Spectra Extraction and Annotation of Liquid Chromatography/Mass Spectrometry Data Sets. *Anal. Chem.* **84**(1), 283–289, https://doi.org/10.1021/ac202450g (2012).

58. Wishart, D. S. et al. HMDB 3.0- The Human Metabolome Database in 2013. *Nucleic Acids Res.* **41**(1), 801–807, https://doi.org/10.1093/nar/gks1065 (2013).

59. Sumner, L. W. et al. Proposed minimum reporting standards for chemical analysis. *Metabolomics*. **3**, 211–221, doi.org/1.1007/s11306-007-0082-2 (2007).

60. Sumner, L. W. et al. Proposed quantitative and alphanumeric metabolite identification metrics. *Metabolomics*. **10**, 1047–1049, https://doi.org/10.1007/s11306-014-0739-6 (2014).

## Acknowledgements

Financial support from the TGIR-RMN-THC Fr3050 CNRS for conducting the research is gratefully acknowledged. N. Wirgot thanks the French Ministry of Research for her PhD scholarship. Part of the work was funded in the frame of the ANR BIOCAP French program. Acknowledgment to people developing the web-based platform Galaxy. We used this platform for LC-MS data treatment and also for the statistical analyses. Barbara Ervens is acknowledged for fruitful discussions and corrections of the manuscript.

## Author Contributions

Anne-Marie Delort designed the experiments and coordinated the work. Nolwenn Wirgot, Marie Lagrée, Mounir Traikia, Ludovic Besaury, Isabelle Canet, Martine Sancelme and Bernard Lyan performed the experiments and analysed the data. Nolwenn Wirgot, Cyril Jousse, Binta Diémé, Pierre Amato and Anne-Marie Delort wrote the manuscript.

## Additional Information

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-019-49319-2.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019