Influence of rehydration on transcriptome during resuscitation of desiccated Pseudomonas putida KT2440

Pseudomonas putida KT2440 is a desiccation-sensitive bacterium that loses culturability after 15 days of air desiccation. We have previously shown that P. putida KT2440 can develop a viable but nonculturable (VBNC) state after being exposed to desiccation stress and eventually recover when desiccated cells are rehydrated for at least 24 h. To determine which genes of transport, oxidation-reduction, and transcription processes could be involved in the return of P. putida KT2440 to the culturable state, a transcriptome analysis was carried out comparing the gene expression of non-desiccated samples with samples subjected to desiccation followed by 20 min of rehydration or desiccation followed by 24 h of rehydration. Desiccation stress triggered a VBNC state of P. putida. The major response was detected after 24 h of rehydration with 148 upregulated and 42 downregulated genes. During the VBNC state, P. putida activated transmembrane transport processes like that of siderophores through a TonB-dependent transporter and putative polyhydric alcohol transport systems. Prolonged rehydration with distilled water resuscitated P. putida KT2440 cells activating the catabolism of phenylalanine/tyrosine to provide energy and carbon for ubiquinone biosynthesis while maintaining a reduced protein synthesis. On the other hand, the interruption of the TonB-dependent receptor gene (PP_1446) increased desiccation survival of the mutant strain. The activation of the iron transport system (TonB-dependent siderophore receptor) and alcohol transport can be helping the VBNC state of P. putida. Activation of catabolism of phenylalanine/tyrosine and reduced protein synthesis was needed for resuscitation from the VBNC state.


Background
Beneficial interactions between plants and bacteria greatly favor plant growth and its protection against diseases and pests (Kloepper et al. 1989;Hayat et al. 2010). A successful interaction not only relies on the capabilities of the beneficial bacteria but also on abiotic factors such as temperature, pH, salinity, and water availability that can disrupt the bacterial-plant interaction (Molina-Romero et al. 2017;Enebe and Babalola 2018). Water limiting conditions affect the viability of the cells and cause loss of membrane integrity, protein denaturation, accumulation of reactive oxygen species (ROS), and DNA damage leading to bacterial death (Potts 1994;Potts et al. 2005;Lebre et al. 2017). Pseudomonas putida KT2440 is a metabolically versatile Gram-negative, gamma-proteobacterium. P. putida KT2440 is a plant growth-promoting rhizobacterium (PGPR) able to colonize the rhizosphere of several plants, degrade organic pollutants, and fight plant diseases (Molina et al. 2000;Espinosa-Urgel et al. 2002;Nelson et al. 2002;Belda et al. 2016).
Despite the metabolic adaptability of P. putida KT2440, its survival decreases drastically under reduced water availability (Manzanera et al. 2002;Pazos-Rojas et al. 2019). P. putida KT2440 loses culturability after 15 days of desiccation in a controlled chamber (50% relative humidity and 30°C) entering into the viable but nonculturable (VBNC) state but can return to the cultivable state after interaction with maize root exudates or 24 h of rehydration (Molina-Romero et al. 2017;Pazos-Rojas et al. 2019). The VBNC state can be considered a long-term survival strategy of bacterial cells or as a deteriorated state of the cells that lead to death (Ducret et al. 2014). Desiccated P. putida cells have impaired culturability and membrane integrity but maintain active gene expression as shown by the expression and synthesis of an introduced GFP reporter protein (Pazos-Rojas et al. 2019).
VBNC bacteria might recover culturability under suitable conditions. For example, Legionella pneumophila is resuscitated by contact with amoebae or by amending agar plates with compounds that reduce the formation of reactive oxygen species (ROS) (Ducret et al. 2014). The culturability of Lactobacillus harbinensis is recovered by the addition of catalase rather than by a temperature upshift (Liu et al. 2018). Escherichia coli and P. aeruginosa resuscitate by incubation in LB broth at 37°C after UV (Zhang et al. 2015) or chlorine treatment (Chen et al. 2018). Acidovorax citrulli and P. aeruginosa regain culturability immediately after copper ion chelation (Bédard et al. 2014;Kan et al. 2019).
Regarding the mechanisms to cope with desiccation stress, some bacteria maintain the integrity of the membrane during drought stress by increasing de novo synthesis of trehalose that acts as an osmoprotectant (Purvis et al. 2005;Zhang and Yan 2012). P. putida and other soil bacteria induce the synthesis of exopolysaccharides and biofilm formation during water limitation to retain water and to ensure a moist microenvironment to handle drying conditions (Roberson and Firestone 1992;Chang et al. 2007). The mechanism through which P. putida KT2440 enter the VBNC state upon desiccation stress and resuscitate after rehydration remains poorly investigated. Here, we studied the transcriptome of P. putida KT2440 subjected to desiccation stress followed by rehydrating to identify genes potentially involved in the resuscitation of VBNC P. putida cells. Results suggest that activation of the iron transport system through a TonB-dependent receptor could be helping the VBNC state of P. putida KT2440.

Results
Culturable cell counts of P. putida KT2440 under air desiccation (30°C and 50% relative humidity) and transcriptional analysis Bacterial samples that were subjected to air desiccation dry completely becoming powder after 6 days. After 15 days of air desiccation, the culturable cell counts of P. putida KT2440 decreased from 10 10 CFU/mL to 0 CFU/ mL (Fig. 1). However, when desiccated bacterial samples were rehydrated for at least 24 h, the culturable cell counts were as high as 10 4 CFU/mL. We previously described that P. putida KT2440 can enter the viable but nonculturable (VBNC) state during desiccation stress and return to the culturable state after 24 h of rehydration (Pazos-Rojas et al. 2019). To determine which genes could be involved in the return of P. putida KT2440 to the culturable state, a transcriptome analysis was carried out comparing the gene expression of non-desiccated samples with samples subjected to desiccation followed by either 20-min rehydration (VBNC state) or 24-h rehydration (resuscitation from the VBNC state). The major response of air-desiccated cells was detected after 24 h of rehydration with 148 upregulated genes and 42 downregulated genes (Fig. 2). The response after 20 min of rehydration was lower in terms of differentially expressed genes, with only 6 upregulated genes. These genes were related to transmembrane transport and oxidation-reduction processes ( Table 1). The complete list of differentially expressed genes after 24 h of rehydration and their fold changes are presented in Table S1 in the supplemental material. A total of 56% of differentially expressed genes after 24 h of rehydration are associated with biological process GO terms. The function of the upregulated genes is mainly involved in transport, oxidationreduction processes, and regulation of transcription (Fig. 3). The functional annotation and overrepresentation of GO terms (P < 0.05) of the upregulated functional group of genes belonging to the transport process included the peptide ABC transporter substrate-binding protein (PP_4146), microcin C ABC transporter substrate-binding protein (PP_4147), membrane protein (PP_3561), EmrB/QacA family drug Fig. 1 Culturable cell counts of P. putida KT2440 subjected to 18 days of desiccation. Means and standard deviation were calculated for triplicates resistance transporter (PP_2067), MFS transporter (PP_ 2651), 4-hydroxybenzoate transporter (PP_1376), ribose ABC transporter ATP-binding subunit (PP_2455), C4dicarboxylate transport protein (PP_1188), pyoverdine ABC transporter ATP-binding protein/permease (PP_4216), and nickel ABC transporter ATP-binding protein (PP_3346) ( Table S1). The functional annotation clustering also revealed the enrichment of gene products involved in regulation of transcription such as RNA polymerase sigma-70 factor (PP_3086), AraC family transcriptional regulator (PP_ 2072), LuxR family transcriptional regulator (PP_0574), transcriptional regulator (PP_0175), 2-ketogluconate utilization repressor(ptxS) (PP_3380), GntR family transcriptional regulators (PP_2253 and PP_0204), LysR family transcriptional regulator (PP_2250), and TetR family transcriptional regulators (PP_4754, PP_2475, PP_2951). Another overrepresented functional group was the biosynthesis of amino acids, particularly phenylalanine and tyrosine biosynthesis which involved the upregulation of shikimate 5-dehydrogenase (PP_ 3768), 3-dehydroquinate synthase (aroB) (PP_5078), and fructose-bisphosphate aldolase (fba) (PP_4960) (Supplemental Figure S3). The upregulation of 2-Cys peroxiredoxin (tpx) (PP_3587), thioredoxin (trxC) (PP_5069), and glutaredoxin (grxC) (PP_5054) suggested activation of the antioxidant defense or cell stress response since these enzymes can also perform the function of molecular chaperones. The downregulated genes after 24 h of rehydration included genes related to translation, oxidation-reduction processes, regulation of transcription, transport, and biosynthetic processes (Fig. 3). The downregulation of PP_2467, PP_0454, PP_0446, PP_0459 encoding for 50S ribosomal proteins, and the downregulation of PP_4709 (rpsO gene) encoding 30S ribosomal protein suggested a limited translation capability of these cells. Furthermore, the PP_4000 gene encoding serine-tRNA ligase was also downregulated. There was downregulation of the LysR family transcriptional regulator (PP_2695), integration host factor subunit beta (gene ihfB, PP_1773), and RNA polymerase sigma 70 factor (gene rpoD, PP_0387). Since RNA polymerase sigma factor RpoD is the primary sigma factor during exponential growth, its repression agrees with the non-growing state of desiccated-24-h rehydrated cells caused by the lack of nutrients in the sterile distilled water. The repression of rpoD and genes encoding 50S and 30S proteins was in line with the fact that sigma 70-RNA polymerase transcribes genes encoding ribosomal proteins. The PP_1914, PP_3839, PP_3433, and PP_4258 genes coding for oxidoreductases were downregulated in desiccated-24-h rehydrated cells (Supplementary Table S1). The downregulation of the PP_1914 gene which codes for the 3-oxoacyl-ACP reductase subunit probably affected fatty acid biosynthesis while the downregulation of PP_3433 (4hydroxyphenylpyruvate dioxygenase) might be impairing the phenylalanine/tyrosine catabolism at the phenylpyruvate node diverting carbon flow to phenylacetaldehyde and phenyllactate nodes (Supplemental Figure S4).
Desiccation tolerance of P. putida KT2440 and its mutants Considering that the six upregulated genes after 20 min of rehydration could be involved in triggering the cascade leading to the resuscitation from the viable but nonculturable state, four for which P. putida KT2440 mutant strains (ΔPP_3451, ΔPP_0880, ΔPP_1446, ΔPP_0256) were available in the Pseudomonas Reference Culture Collection Fig. 2 Venn diagrams showing statistically significant (p < 0.05) up-and downregulated genes in the desiccated-rehydrated cells relative to nondesiccated control cells from P. putida KT2440 in microarray analyses. Values shown represent more than 2-fold differential expression. Twenty minutes and 24 h indicate the bacterial cells which were desiccated for 18 days and then rehydrated for 20 min or 24 h, respectively, before RNA extraction. Microarray analyses were performed in triplicate for each condition of the EEZ-CSIC (http://artemisa.eez.csic.es/prcc/) were studied in greater detail to evaluate their involvement in desiccation-rehydration tolerance (Fig. 4). The ΔPP_0880 and ΔPP_1446 mutants are interrupted in a dppCdipeptide ABC transporter permease and a TonBdependent receptor, respectively. The PP_0256 gene encodes a molybdopterin oxidoreductase subunit and the PP_ 3451 gene encodes a hypothetical protein (Table 1). There was no difference in the desiccation tolerance between the mutant strains and the wild type as almost all dropped to 0 CFU/mL after 15 days of desiccation, except for the ΔPP_ 1446 mutant which survived 18 days of desiccation with 10 3 CFU/mL (Fig. 4). When 18 days-desiccated samples were rehydrated for 24 h, the Log 10 culturable cells (CFU/ mL) were 3.17 ± 0.03, 2.81 ± 0.30, 3.32 ± 0.02, 3.66 ± 0.07, 3.50 ± 0.19 for P. putida KT2440, ΔPP_0256, ΔPP_3451, ΔPP_0880, and ΔPP_1446, respectively. This result suggests that none of the mutations were essential for the resuscitation of the viable but non-culturable state. Exposure of ΔPP_3451, ΔPP_0880, and ΔPP_0256 mutants to 18 days of desiccation resulted in similar cell viability relative to the wild-type cells subjected to identical conditions (Fig. 4); therefore, we decided to explore if these selected genes that were found to be upregulated upon 20min rehydration were also overexpressed in desiccated cells prior rehydration. Thus, the qRT-PCR analysis was performed for the parental strain subjected to 18 days of desiccation (Fig. 5). Compared with non-desiccated samples (BD), the transcription levels of PP_1446 and PP_0256 genes increased significantly in the desiccated cells without rehydration (WoR condition). This overexpression of PP_ 1446 (TonB-dependent receptor) and PP_0256 (molybdopterin oxidoreductase subunit) was maintained regardless of posterior 20-min or 24-h rehydration (Fig. 5). PP_3451 gene expression increased significantly only after 20 min and 24 h of rehydration while the expression of PP_0880 did not increase under any of the conditions evaluated. A qRT-PCR analysis of PP_2676 and PP_2662 was also performed (Supplemental Figure S1) and the results show that both are expressed in desiccated cells and after rehydration periods, thereby both displayed a similar expression trend to that observed by microarray (Table 1). In summary, five of the six genes displayed a similar expression profile obtained by both microarray and qRT-PCR at 20-min rehydration.

Discussion
In other research works where P. putida is subjected to mild water limitation, bacteria upregulate gene expression of processes related to transcription, biosynthesis, and degradation of surface polysaccharides, and cell envelope; while differentially regulate the transport and binding proteins (Gülez et al. 2012;Gulez et al. 2014). In this work, desiccated P. putida KT2440 cells upon 24 h of rehydration upregulated biological processes related to transport, oxidation-reduction, regulation of transcription, and biosynthetic process while downregulated biological processes related with translation, oxidationreduction, and the regulation of transcription. This upand downregulation of genes within similar role categories is common in the water stress response of P. putida KT2440 (Gülez et al. 2012). Genes encoding for products classified within the same biological process can act on different pathways. Thus, PP_3768 and PP_3839 belonging to the same biological process GO term, oxidation-reduction process, were up-and downregulated, respectively (Fig. 3, Supplemental Table S1) probably playing roles within specific pathways. The PP_3768 gene encoding shikimate 5-dehydrogenase was upregulated because it is required for phenylalanine/tyrosine biosynthesis while the PP_3839 gene encoding alcohol dehydrogenase was downregulated to decrease, probably, the reduction and oxidation of aldehydes and alcohols.
During water stress, Pseudomonas cells produce alginate to create a hydrated microenvironment while accumulating intracellular organic osmolytes like trehalose, betaine, ectoine, N-acetylglutaminylglutamine amide, Ala, Gly, Glu, Ser, Thr, Asp, glucosylglycerol, and mannitol (Pocard et al. 1994;Kets et al. 1996;Paul and Nair 2008;Kurz et al. 2010). After 18 days of desiccation, desiccated P. putida cells lost their cultivability and probably accumulated a high number of compatible solutes. Twenty-four hours of rehydration with sterile distilled water resuscitated dormant cells with bacterial counts of 10 3 -10 4 CFU/mL but it could also cause a hypoosmotic shock driving water into the cells. To cope with the resulting increase of turgor pressure, bacteria probably release the stockpile osmolytes to reduce the driving force for water entry. The sudden addition of water to the environment, such as during a rainstorm, is a typical hypoosmotic stress faced by soil bacteria and it was exactly what we were looking for with this experiment. It is important to mention that distilled water does not lyse P. putida cells as observed by transmission electron microscopy analysis (Pazos-Rojas et al. 2019). Thus, after 24 h of starvation and osmotic down-shock, P. putida activated the RrbsA-I ribose ABC transporter (PP_2455), GlnL/GlnG twocomponent system sensor histidine kinase/phosphatase (PP_5047) which senses N 2 availability, and the SpuB Glutamylpolyamine synthetase (PP_5183), and DctA C4 dicarboxylate (succinate, malate, L-aspartate, oxaloacetate) transport protein (PP_1188) for carbon and nitrogen uptake. Although distilled water contains no nutrients, trace amounts of ribose, amino acids, and organic compounds released by dead cells or during the hypoosmotic shock could be available in this medium. Starvation induces the stringent response producing the alarmones guanosine tetraphosphate (ppGpp) and guanosine pentaphosphate (pppGpp) through the enzymes RelA and SpoT homologs (Liu et al. 2017). Although the transcription levels of relA and spoT were not affected after 24 h of rehydration, the inhibition of ribosome protein biosynthesis and repression of RpoD, two targets of (p)ppGpp, suggest an active stringent response in Fig. 5 qRT-PCR analysis of genes PP_3451, PP_0880, PP_1446, and PP_0256 in desiccated P. putida KT2440 wild-type cells followed by rehydration. "BD" means before desiccation and refers to non-desiccated cells. WoR indicates that RNA was extracted from desiccated bacteria prior to rehydration. Twenty minutes and 24 h indicate that 18 days-desiccated bacteria were rehydrated for 20 min or 24 h before RNA extraction. Error bars represent the standard deviation of three replicates. *Indicates statistical significance (p < 0.005) these rehydrated cells (Zhang et al. 2019). Besides, the activation of pyoverdine transporter, a siderophore upregulated by (p)ppGpp, agrees with the stringent control (Vogt et al. 2011;Zhang et al. 2019). Other potential starvation-survival strategy used by desiccated-24-h rehydrated cells was the activation of phenylalanine/tyrosine biosynthesis (Supplemental Figure S3) and phenylalanine degradation (Supplemental Figure S4). This result suggests that desiccatedrehydrated P. putida cells obtain nitrogen from the catabolism of this amino acid (Herrera et al. 2012). Furthermore, the repression of pde encoding 3',5'-cyclic-AMP phosphodiesterase (PP_4917) can be increasing the cAMP pool mediating the utilization of phenylalanine and tyrosine as nitrogen sources by cAMP-CRP (Herrera et al. 2012;Green et al. 2014). The carbon flow coming from phenylalanine can be deflecting, at the phenylpyruvate node, to ubiquinone biosynthesis, which agrees with the activation of UbiA (PP_ 5318), the regulatory enzyme of the overall ubiquinone pathway (Supplemental Figure S4). Ubiquinone is used as a cofactor of quinone-dependent alcohol dehydrogenases (PP_2674, PP_2679) and electron carriers in the oxidative phosphorylation. The complex I (NADH dehydrogenase) of 24-h rehydrated cells was activated as suggested by the upregulation of the nuoG gene (PP_4124) while the CcoP-II gene of cytochrome c oxidase complex was repressed indicating that the electron transport chain was working partially and that the energy state of the cells was low. Finally, the reduced protein synthesis in desiccated-24-h rehydrated cells, as suggested by the repression of genes encoding 50S and 30S ribosomal proteins, is a hallmark of "persister" or VBNC cells (Song and Wood 2020).
P. putida KT2440 lost cultivability after 15 days of air desiccation ( Fig. 1), entering the VBNC state while maintaining the upregulation of five of six genes detected after 20 min of rehydration (Table 1). The function of these genes could be of importance to maintain the viable but nonculturable (VBNC) state or to resuscitate from the VBNC state since five of these genes were highly upregulated in desiccated (VBNC cells) and desiccated-shortly rehydrated cells (Fig. 5, supplemental Figure S1). It is interesting that PP_2676, a gene encoding for a substrate-binding protein of an ABC importer, was upregulated in all conditions tested (desiccated, desiccated-20-min rehydrated, and desiccated-24-h rehydrated cells) by both microarray and qRT-PCR analysis. Unfortunately, the ΔPP_2776 mutant strain was not available to evaluate its behavior under desiccation stress. The TonB-dependent receptor (PP_1446), substrate-binding protein (PP_2676), ethylene glycol porin (PP_2662), hypothetical protein (PP_3451), and molybdopterin oxidoreductase subunit (PP_0256) genes were found to be upregulated by the qRT-PCR analysis in both desiccated (WoR) and desiccated-20-min rehydrated cells (Fig. 5, Supplemental Figure S1). Molybdopterin oxidoreductases play essential roles in bacteria, some of them are required for persistence during pathogenesis, protection of periplasmic proteins from oxidative damage, and degradation of nicotine (Williams et al. 2015;Gennaris et al. 2015;Xia et al. 2018). The role of molybdopterin oxidoreductase in desiccated P. putida cells is unknown but it will be interesting to explore any protecting effect on proteins of the periplasm against oxidation stress. The porin, PP_2662, and the substratebinding protein, PP_2676, are involved in the exchange of ethylene glycol and its oxidation products across the outer membrane during the degradation of ethylene glycol and other alcohols via quinoprotein alcohol dehydrogenases (Mückschel et al. 2012;Vallon et al. 2015). The porin, PP_ 2662, and quinoprotein alcohol dehydrogenases PdeE and PdeH (PP_2674, PP_2679) of P. putida KT2440 are strongly induced by ethylene glycol but not by its oxidized byproduct glyoxylic acid (Mückschel et al. 2012). The substratebinding protein, PP_2676, under control of the PedR1 regulator, is highly induced in response to chloramphenicol and naphthalene (Fernandez et al. 2012). We presume that PP_ 2662 and PP_2676 could be involved in the transport and degradation of polyhydric alcohols, e.g., glucosylglycerol, which are accumulated during desiccation stress as compatible solutes. The nutrient starvation experienced by P. putida cells during the 18 days of desiccation probably increased the (p)ppGpp pool, a positive regulator of TonBdependent transporters, increasing the uptake of siderophore, vitamin B12, nickel complexes, or carbohydrates (Schauer et al. 2008;Zhang et al. 2019). Since water removal triggers iron oxidation of proteins (Labrude et al. 1987;Kamnev et al. 2020), desiccated cells, in turn, activated the TonB-dependent receptor (PP_1446) probably to scavenge iron. During dormant states, other bacteria also maintain active transcription of genes encoding antioxidant defenses and ion transport. For example, VBNC Vibrio cholerae cells upregulate genes responsible for the transport of iron (iron (III) ABC transporter) and preservation of DNA metabolism (Asakura et al. 2007). Novosphingobium sp. LH128 dormant cells overexpress genes encoding for catalases and 1-Cys peroxiredoxin and genes involved in the transport of ions across the membrane, such as homologs of cation/multidrug efflux pump proteins, Na+/H+ 327 antiporter, copper resistance protein, and bacterioferritin (Fida et al. 2017). In line with this need of activating an antioxidant response, supplementation of culture medium with the antioxidant oxyrase or ROS scavengers resuscitates enterohemorrhagic Escherichia coli and Legionella pneumophila from the VBNC state (Reissbrodt et al. 2002;Ducret et al. 2014) suggesting that bacterial cells could be encountering oxidative stress during the formation of the VBNC state. Thus, relieving the oxidative stress with catalase or non-enzyme peroxide-degrading compounds restores cell culturability (Mizunoe et al. 1999;Liu et al. 2018). The TonB-dependent siderophore receptor is involved in the uptake of essential nutrients, including iron, which is essential for biofilm formation (Gulez et al. 2014;Luscher et al. 2018). Iron might be involved in the resuscitation of non-culturable bacteria since medium supplementation with ferrioxamine E resuscitates the VBNC Salmonella enterica cells (Reissbrodt et al. 2002). The mechanism of this effect is not clear. The authors proposed that ferrioxamine E can modulate the release and use of ferric iron preventing the generation oxygen free radicals (Reissbrodt et al. 2000;Reissbrodt et al. 2002). We would propose for desiccated P. putida cells (VBNC cells) that the accumulation of elevated ROS concentrations oxidize the Fe(II) catalytic cofactor of a wide variety of non-redox enzymes and the [4Fe-4S] clusters of dehydratases triggering the mismetallation of these enzymes thereby disrupting metabolism (Gu and Imlay 2013). The bacteria in response activate iron-import, normally with ferric iron-chelating substances under aerobic conditions but delivering ferrous iron into the cytoplasm, to restore the activity of mismetallated enzymes thus preserving the fundamental cellular functions of cells (Beauchene et al. 2017). Bacteria do not activate ferrous iron transport systems under aerobic conditions but reduce a significant proportion of ferric-iron assimilated from the media to ferrous-iron during aerobic iron assimilation (Kamnev et al. 2014;Beauchene et al. 2017). Therefore, the activation of the ferriciron transport system observed in this work by the overexpression of the TonB-dependent receptor could be a potential strategy used by P. putida KT2440 to provide the ferrous-iron required for vital functions during the VBNC state. In fact, activation of ferric iron-import genes is one of the strategies used by E. coli for the re-metallation of superoxide damaged mononuclear iron enzymes. The resultant apoproteins can either quickly re-bind iron or inappropriately bind zinc (Gu and Imlay 2013). Curiously, a mutant lacking the TonB-dependent receptor gene (PP_1446) was more resistant to desiccation so that cells remained viable longer (Fig. 4). Since intracellular iron catalyzes the Fenton chemistry, the scarce amount of intracellular iron accumulated into the TonB-dependent receptor mutant cells could be limiting the ROS and Fenton chemistry damage during desiccation stress, enabling the PP_1446 mutant strain to survive longer than the parental strain (Frawley and Fang 2014). Although this explanation could appear paradoxical, it is not. We speculate that the iron import during the desiccated state does not provide the substrate for Fenton reactions but cofactors for un-metallated apoproteins through an unknown mechanism that directly delivers the Fe 2+ into the active sites of damaged iron enzymes. Further experiments should be performed to explain the observed behavior of the PP_1446 mutant strain. Finally, we cannot discharge transcriptional changes associated with the change from LB to water. Even stationary phase cells of P. putida were washed, resuspended in distilled water, and maintained in this milieu for 1 h before RNA extraction of "day 0" samples. We cannot discharge the putative transcriptional changes associated with nutritional stress, low-salt stress, and iron starvation induced by transfer cells from LB medium to distilled water. So, the findings of this study can carry out a potential response associated with the mentioned stresses. Further experiments with stationary phase cells resuspended in distilled water and maintained in this setting (without desiccation) for 18 days as non-desiccated control cells will be required to discard the transcriptional changes associated with the change from LB to water.

Conclusions
Desiccation causes loss of culturability of P. putida KT2440 making cells enter the VBNC state as a longterm survival strategy. During this dormant state, P. putida maintains active transmembrane transport processes like siderophore and putative polyhydric alcohol transport systems. Prolonged rehydration with distilled water resuscitates P. putida KT2440 cells activating the catabolism of phenylalanine/tyrosine to provide energy and carbon for ubiquinone biosynthesis while maintaining reduced protein synthesis.

Bacterial strains and growth conditions
Pseudomonas putida KT2440, a plasmid-free derivative of P. putida mt-2 was used in the present study (Nelson et al. 2002). Mutant strains were donated by the Consejo Superior de Investigaciones Científicas, Estación Experimental del Zaidín (CSIC-EEZ). The collection of mutants of P. putida KT2440 was obtained by random mutagenesis using the miniTn5 (Km) transposon. The transposon insertion was identified by arbitrary PCR . The strains used in this work are shown in Table 2. P. putida KT2440 was grown at 30°C in LB (Luria-Bertani) liquid medium supplemented with 100 μg/mL chloramphenicol and constant agitation. Mutant strains were grown in LB medium supplemented with 100 μg/mL chloramphenicol and 50 μg/mL kanamycin.
Bacterial desiccation assays P. putida KT2440 strain and mutant strains were grown until the stationary phase in LB liquid medium supplemented with antibiotics because at this growth phase cells exhibit higher survival to desiccation than cells in the exponential phase (Muñoz-Rojas et al. 2006). Thus, cultures were incubated for 24 h to ensure cells had reached the stationary phase as shown in the growth curve in the supplementary Figure S2. The growth curve of wild type and mutant strains was similar, all of them reached Log (CFU/mL) = 9.77 ± 0.37 in 24 h. Three independent experiments of desiccation assays were performed as described in an established protocol (Molina-Romero et al. 2017). The cells were centrifuged, washed twice with sterile distilled water, and resuspended in distilled water with the initial volume. Subsequently, the cell suspension was aliquoted in 1.5 mL tubes with 500 μL each tube and covered with a sterile cotton cap. The tubes with cell suspension were placed in a desiccation chamber for 18 days at 30°C and 50% relative humidity. Every 3 days after the beginning of desiccation, bacterial survival was determined using five sample tubes. Desiccated samples were rehydrated for 20 min with sterile distilled water to the initial volume. After letting the biomass-powder suspensions sit at room temperature for 20 min, samples were mixed by pipetting up and down for 10 s, then serially diluted and streaked onto LB-agar plates. Twenty minutes of rehydration was the minimal time necessary to disaggregate the agglomeration of dried cells. After 18 days of desiccation, samples were rehydrated for 20 min or 24 h. Twenty-four hours of rehydration was the minimal time required for dormant P. putida cells to return to the cultivable state according to the time profile of bacterial resuscitation (Pazos-Rojas et al. 2019). The number of CFU/ml was monitored by the drop plate method. Three replicates per sample were used.
RNA extraction, reverse transcription (RT), and labeling for microarrays RNA extraction: Total RNA was extracted from P. putida KT2440 at day 0 (these are non-desiccated control cells that were grown until stationary phase, washed, aliquoted into 1.5 mL tubes, and maintained at this milieu 1 h) and from cells after 18 days of desiccation and which had been rehydrated for 20 min or 24 h before RNA extraction. A pool of 30 sample tubes was used for RNA extraction of each condition. The cell suspensions were extracted using TRI Reagent (Ambion) and BCP (1-Bromo-3-chloropropane); RNA was precipitated with isopropanol and washed with 70% ethanol. Finally, RNA was resuspended in DEPC water. To eliminate DNA contamination, the samples were treated with TURBO DNA-free kit (Invitrogen). RNA concentration and purity were determined using a NanoDrop 2000/ 2000c spectrophotometer (Thermo Scientific).
cDNA Synthesis and labeling cDNA probes were generated by reverse transcription using SuperScript II (Invitrogen). The reaction mixture contained 20 μg of total RNA and 7.5 μg of pd(N)6 random hexamer primer (Amersham). Alkaline hydrolysis was used to eliminate non-transcript RNA in cDNA. Primers and dNTPs not incorporated were eliminated using the PCR purification kit (Qiagen). Then, cDNA was labeled: Cy3 fluorophore (CyTM3 monofunctional reactive dye; Amersham Pharmacia) was used to label samples before desiccation and CyTM5 fluorophore (CyTM5 monofunctional reactive dye; Amersham Pharmacia) was used to label samples after desiccation.

Hybridization conditions and microarrays data analysis
Hybridization: Microarray slides (Progenika, Biopharma) were cleaned and immersed in a stop solution (5× SSC; 0.1% SDS; 1% BSA) and washed twice with milli-Q water (Rodríguez-Herva et al. 2010). Labeled cDNA was rehydrated with 35 μL of hybridization solution (5× SSC; 25% (v/v) formamide; 0.5% (p/v) SDS; 5%(p/v) dextran sulfate; 5× solution Denhardt's: Ficoll; polyvinylpyrrolidone; BSA). The mixture was placed on microarray slides and put in a hybridization chamber (ArrayIt Hybridization Cassette) for 20 h at 45°C immersed in a water bath. For each condition evaluated, three slides were used to eliminate irregular data. Microarray data acquisition and analysis: The slide was scanned in a 4100A Biochip Scanner (Axon Corporation, CA) and image data was analyzed in GenePix Pro 5.1 software. Differential gene expression and statistical analysis were performed with AlmaZen Software (Bioalma). Gene expression data were normalized with the following parameters: foldchange <-2 AND exp_level >64 AND p value >0.05 OR foldchange >2 AND exp_level >64 AND p value <0.05. Genes that showed differential expression were classified by Gene Ontology (The Pseudomonas Genome Database, http://www.pseudomonas.com) using biological process terms. Enrichment analysis was performed using the free online platform DAVID (Huang et al. 2018).

Real-time quantitative reverse transcription-PCR (qRT-PCR)
To validate microarray results upon 20 min rehydration, quantitative reverse transcription-PCR (qRT-PCR) was performed using the primers listed in Table 3. RNA extraction was performed using the hot acidic phenolchloroform method and genomic DNA was removed using the Invitrogen TURBO DNA-free Kit (Thermo Fisher Scientific) (Baez et al. 2009). Total RNA was extracted from P. putida KT2440 at day 0 (these are nondesiccated control cells that were grown until stationary phase, washed, aliquoted into 1.5 ml tubes, and maintained at this setting for 1 h), and from cells subjected to 18 days of desiccation. Thus, RNA was isolated from desiccated cells without rehydration (WoR), desiccated cells with 20-min rehydration, and desiccated cells with 24 h rehydration. A pool of 15 sample tubes was used for RNA extraction of each condition. cDNA was generated by using 2 μg RNA in a total volume of 20 μL with 10 μM of specific DNA primers (antisense primers Table 3) according to the protocol of the High Capacity cDNA Reverse Transcription Kit (Applied Biosystems). Real-time qRT-PCR was performed using an Applied Biosystems 7500 Fast Real-Time PCR System. SYBR Green PCR Master Mix was used as a signal reporter. Each reaction was composed of 20 ng of cDNA and 6 μM of each pair of primers in a total volume of 20 μL. The following amplification conditions were used: 1 cycle of 5 min at 95°C followed by 30 cycles of denaturation at 95°C for 30 s and annealing and extension at 55°C for 30 s. Each reaction was performed in triplicate. The expression of the rpoN gene (encoding the sigma factor 54) was used as an endogenous control to normalize the amount of mRNA obtained from a target gene (Yuste et al. 2006;Follonier et al. 2013). Data were analyzed using the 2 −ΔΔC T method (Livak and Schmittgen 2001).

Statistical analysis
The statistical analysis of qRT-PCR results was performed using Sigma Plot (Handel Scientific Software). Data were analyzed globally by ANOVA single factor, and significant treatment effects were determined by the Tukey-Kramer post hoc test. the data and designed the figures. LILL and LAPR conducted the experiments. LELC performed qRT-PCR of desiccated cells. YEMG and VQH supervised the work and administrated the project. LAPR, JDLT, and PVD manufactured the samples and performed the microarray experiments. AB took the lead in writing the manuscript. The author(s) read and approved the final manuscript.

Funding
Funding was provided by the intramural program of the Vicerrectoria de Investigación y Estudios de Posgrado, BUAP through grants PV16-ID00450, PV17-ID00337. Lilia Isela Lopez Lara received a scholarship (number 273274) from the National Council of Science and Technology (CONACyT).

Competing interests
The authors declare that they have no conflicts of interest.

Availability of data and materials
The complete list of differentially expressed genes at 24 h rehydration and their fold change are presented in Table S1.