# The PEG-responding desiccome of the alder microsymbiont Frankia alni

Kais Ghedira, Emna Harigua-Souiai, Cherif Ben Hamda, Pascale Fournier, Petar Pujic, Sihem Guesmi, Ikram Guizani, Guylaine Miotello, J. Armengaud, Philippe Normand, et al.

► To cite this version:

Kais Ghedira, Emna Harigua-Souiai, Cherif Ben Hamda, Pascale Fournier, Petar Pujic, et al. The PEG-responding desiccome of the alder microsymbiont Frankia alni. Scientific Reports, 2018, 8 (1), pp.759 - 759. 10.1038/s41598-017-18839-0. pasteur-01729024

HAL Id: pasteur-01729024

https://riip.hal.science/pasteur-01729024

Submitted on 12 Mar 2018

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
SCIENTIFIC REPORTS

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

OPEN

# The PEG-responding desiccome of the alder microsymbiont *Frankia alni*

Received: 18 September 2017

Accepted: 5 December 2017

Published online: 15 January 2018

Kais Ghedira<fnref n="1" />, Emna Harigua-Souiai ID², Cherif Ben Hamda<fnref n="1" />,³, Pascale Fournier⁴, Petar Pujic⁴, Sihem Guesmi⁵,⁶, Ikram Guizani², Guylaine Miotello⁷, Jean Armengaud⁷, Philippe Normand ID⁴ & Haïtham Sghaier ID⁵,⁸

Actinorhizal plants are ecologically and economically important. Symbiosis with nitrogen-fixing bacteria allows these woody dicotyledonous plants to colonise soils under nitrogen deficiency, water-stress or other extreme conditions. However, proteins involved in xerotolerance of symbiotic microorganisms have yet to be identified. Here we characterise the polyethylene glycol (PEG)-responding desiccome from the most geographically widespread Gram-positive nitrogen-fixing plant symbiont, *Frankia alni*, by next-generation proteomics, taking advantage of a Q-Exactive HF tandem mass spectrometer equipped with an ultra-high-field Orbitrap analyser. A total of 2,052 proteins were detected and quantified. Under osmotic stress, PEG-grown *F. alni* cells increased the abundance of envelope-associated proteins like ABC transporters, mechano-sensitive ion channels and Clustered Regularly Interspaced Short Palindromic Repeats CRISPR-associated (cas) components. Conjointly, dispensable pathways, like nitrogen fixation, aerobic respiration and homologous recombination, were markedly down-regulated. Molecular modelling and docking simulations suggested that the PEG is acting on *Frankia* partly by filling the inner part of an up-regulated osmotic-stress large conductance mechanosensitive channel.

Actinobacteria belonging to the genus *Frankia* do establish nitrogen-fixing nodular symbiosis with the roots of 23 angiosperm genera that are collectively called “actinorhizals”<fnref n="1" />. These plants form root nodules in which *Frankia* fixes nitrogen, thus permitting them to thrive in pioneer soils poor in nitrogen and organic matter, such as glacial moraines, lava fields, forest burnouts or anthropogenic sites such as mine spoils or hydrodam dykes². *Frankia* establishes a symbiotic association with the roots of several dicotyledonous plants. The different *Frankia* lineages form a coherent cluster at the root of the aerobic actinobacteria phylum³, and *F. alni* in particular establishes symbiosis with alder (*Alnus*) and bayberry (*Morella*) species⁴.

The interaction has evolved over several million years with a sophisticated dialogue that does not imply acylated N-acetyl-glucosamine oligomeric Nod factors⁵. The *Frankia* determinants of symbiosis are still poorly known, which is for the most part due to the lack of a genetic transformation system. Transcriptomics has shown genes coding for nitrogenase (*nif*), hydrogenase uptake (*hup*), hopanoids (*shc*, *hpn*), iron-sulfur (*suf*) clusters as among the most up-regulated⁶. Proteomics has also been used to analyze the symbiosis⁷ to identify up-regulated peptides and it has shown the presence of Nif, Hup, Suf, Hop proteins as expected but also several transporters, regulators and various proteins involved in stress responses.

<footnote n="1">¹Laboratory of Bioinformatics, Biomathematics and Biostatistics - LR16IPT09, Institut Pasteur de Tunis, Université de Tunis el Manar, Tunis, 1002, Tunisia. ²Laboratory of Molecular Epidemiology and Experimental Pathology - LR11IPT04, Institut Pasteur de Tunis, Université de Tunis el Manar, Tunis, 1002, Tunisia. ³Université de Lyon, Université Lyon 1, Lyon; CNRS, UMR 5557, Ecologie Microbienne, UMR1418, INRA, 69622 Cedex, Villeurbanne, France. ⁴Laboratory "Energy and Matter for Development of Nuclear Sciences" (LR16CNSTN02), National Center for Nuclear Sciences and Technology (CNSTN), Sidi Thabet Technopark, 2020, Tunisia. ⁵National Agronomy Institute (INAT), Avenue Charles Nicolle, 1082, Tunis, Mahrajène, Tunisia. ⁶Laboratoire Innovations Technologiques pour la Détection et le Diagnostic (Li2D), Service de Pharmacologie et Immunoanalyse (SPI), CEA, INRA, F-30207, Bagnols sur Cèze, France. ⁷Associated with Laboratory "Biotechnology and Nuclear Technology" (LR16CNSTN01) & Laboratory "Biotechnology and Bio-Geo Resources Valorization" (LR11ES31), Sidi Thabet Technopark, 2020, Tunisia. Correspondence and requests for materials should be addressed to P.N. (email: philippe.normand@univ-lyon1.fr)</footnote>
Osmotic stress is a constant challenge for bacteria living in a range of soils⁸,⁹, from those having low salinity in rainy cold latitudes to highly saline ones in warm and dry parts of the world. It also affects symbionts¹⁰ that must alternate between two main biotopes, the soil and the root tissues. Soils, especially poor soils, have low osmotic potential while plant tissues have a much higher osmotic potential. Most studies on the effect of salinity on actinorhizal symbionts have focused on the plant, especially *Casuarina* to illustrate various adaptations such as sodium partitioning or proline accumulation¹¹. Only a few recent studies on *Frankia* have shown the effects of salinity on ammonium (NH₄⁺) assimilation¹², on cell wall/membrane biogenesis functions and on some transport proteins¹³. *Alnus*-, *Casuarina*- and *Elaeagnus*- infective isolates grow best at 50 mM but are nevertheless able to grow well in medium containing up to 200 mM NaCl, but not at 500 mM¹⁴. The effect of NaCl on one strain, Cc16, was apparently less severe than that of another osmolyte, sucrose¹⁵. Nevertheless, little is known about the molecular adaptations and the present study was undertaken to better understand how *Frankia* coped with osmolytes using proteogenomics, an approach recently used to better decipher how *Geodermatophilaceae* coped with desiccation¹⁶.

NaCl causes various impairments to cells, some related to osmotic stress and some related to ionic stress and these two effects are clearly distinguishable¹⁷. PEG is a non-ionic polymer of various molecular weights (MW) that has been used often for osmotic studies when it is sought to have only the osmotic effect.

Mechanisms involved in desiccation tolerance have been well studied in Cyanobacteria¹⁸, Actinobacteria¹⁹,²⁰ as well as xerotolerant members of the *Bosea*, *Chelatococcus*, *Deinococcus* and *Methylobacterium* genera²¹. Briefly, the different studies have shown that desiccation damages cell membranes, proteins and DNA due to oxidative stress and production of reactive oxygen species (ROS)²⁰,²¹. Strategies adopted by these microorganisms for desiccation tolerance include prevention of ROS damage, osmoprotection through the accumulation of sucrose and/or trehalose¹⁰,²², the uptake of exogenous glycine betaine²³, the production of extracellular polysaccharides (EPS) that reduce water loss²⁴, or the ability to limit protein oxidation during dehydration²¹. Previous studies⁸,²⁰,²¹,²⁵⁻²⁷ highlighted changes in expression of several genes following water stress including genes for synthesis of trehalose/sucrose, sugar transporters, chaperone genes (*groES/EL, dnaK/J*), oxidative stress protection genes (*dps*, thioneredoxin), ABC transporters, dehydrogenases, esterases, proteases, hydrolases and lyases. A *Casuarina*-infective *Frankia* strain was recently shown to have few proteins upregulated under salt stress, among which cell wall/membrane biogenesis functions and some transport proteins¹³.

The desiccome, a term coined by Potts and colleagues²⁸, can be defined as the set of genes, proteins and metabolites that are necessary for desiccation tolerance²⁹. Here we report the first PEG-responding desiccome of a nitrogen-fixing symbiotic bacterium, *F. alni*. We highlight the role of ABC transporters, mechanosensitive ion channels and Cas components in response to this stress.

## Results

**Growth in polyethylene glycol (PEG).** Upon inoculation in BAP- medium at a 0.1 OD₆₀₀, growth of *F. alni* strain ACN14a was inhibited at all concentrations of PEG tested. Nitrogen fixation was active without PEG but even the lowest concentration of PEG completely inhibited it (data not shown). We chose an inoculation density of 0.1 to have sufficient biomass and after 7 days, growth without PEG had reached an OD₆₀₀ of 0.677, i.e. a 0.74 fold in density increase. In the presence of PEG 0.9%, OD₆₀₀ reached 0.42 or a 0.47 fold increase. Growth in the presence of PEG was thus 63% that without PEG. Hyphae, vesicles and sporangia looked similar under both conditions.

**Comprehensive proteome coverage of *F. alni* cells.** Thanks to the high speed and high resolution of this analyser, a total of 295,788 high quality MS/MS spectra were recorded. Amongst these, a set of 208,018 MS/MS spectra (70.3%) was assigned to peptide sequences. This dataset allowed detecting a total of 20,825 peptide sequences, which pointed at 2,454 polypeptides, with 2,052 proteins certified with at least two distinct peptide sequences. These proteins and their characteristics are listed in Supplementary Data, Table S1. On the average, a total of 20,801 spectral counts were measured per sample with a remarkably low standard deviation (5.3%). With an average of 10 peptides and 100 spectral counts per polypeptide, an extensive coverage of the proteome was reached that compares favourably to most proteomic studies of bacteria¹⁶,³⁰,³¹.

**Differentially expressed proteins.** When comparing both conditions and compiling the five biological replicates of each condition, a total of 1,951 proteins were detected whatever the conditions. A set of 29 proteins was specifically detected only in the reference while 85 were detected only in the PEG-treated cells. Differentially detected protein abundances were identified using a Fold Change (FC) threshold of 1.5 and *p*-value below 0.05. Four classes were delineated as follows: i) Blue class proteins for which identifications satisfied both, the fold change (>1.5) and statistical criteria (*p*-value < 0.05); ii) Orange class proteins for which identifications did not meet the fold criterion but have low *p*-values; iii) Green class proteins for which identifications satisfied the fold criterion but not the statistical criterion; iv) and finally Red class for which identifications did not meet the fold and *p*-value criteria. The Blue class comprised 294 proteins. A total of 211 proteins were more abundant in PEG-treated cells compared to the reference while 83 were found less abundant. Figure 1 shows the distribution of the proteins affected by the PEG treatment compared to the control within the *F. alni* genome and their expression levels. Genes coding for more abundant proteins upon PEG-induction are highlighted in red whereas those more abundant in the control were highlighted in cyan in Fig. 1. Tables 1 and 2 list the most up- and down-regulated PEG-modulated proteins of *F. alni*, respectively. Twelve ribosome-associated translation proteins (FRAAL1099-50S ribosomal subunit protein L30 (2.10), *etc.* ) were also up-regulated by PEG but are not shown in Table 1.

**Re-annotation of differentially expressed proteins following polyethylene glycol (PEG) treatment.** Among the blue class differentially expressed proteins (294 proteins) identified under PEG treatment,
<figure><img src="image_3.png" /><figcaption>**Figure 1.** *F. alni* proteome expression following polyethylene glycol (PEG) treatment. From the outside in are 1-the coordinates in bp; 2-the coordinates in 100 genes; 3-the symbiotic genes (in green the hopanoid genes, in blue the uptake hydrogenase genes, in yellow the iron-sulfur genes, in gray the cellulase genes, in red the nitrogenase genes); 4-selected genes that play a role in stress response with oxidative stress genes (orange), heat-resistance genes (red), osmotic-resistance genes (cyan), UV-resistance (violet), metals-resistance (grey), Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) genes (black), poorly characterised or universal stress genes (yellow); 5-Genes coding for proteins more abundant upon PEG-induction (red), or more abundant in the control (cyan); 6-Symbiosis up-regulated genes (in green) from Alloisio *et al.*<sup>6</sup>; 7-Core genome at a level of 70% of conserved amino acids over 80% of the length of the shortest sequence in selected symbiotic *Frankia* strains of cluster1a (ACN14a, CpI1, QA3), cluster1c (CcI3, CcI6), cluster2 (Dg1, BMG5.1) and cluster3 (EAN1pec, BMG5.12, BCU110501), in yellow.</figcaption></figure>

53 were characterised as hypothetical proteins and/or conserved hypothetical proteins. We therefore conducted a new analysis using PSI-BLAST tool to re annotate them. Seventeen ORFs remained as proteins with unknown function because of their low similarity to characterised proteins and the function of 36 ORFs could be deduced from closely-related known proteins (Supplementary Data, Table S2. Re-annotation of differentially expressed proteins following polyethylene glycol (PEG) treatment based on a Psi-Blast approach.). Among re-annotated ORFs, several envelope-associated proteins were identified including the most up-regulated protein (FRAAL2567 (5.67)), a putative auto-transporter adhesin. These results consolidate previous data (see Table 1) that have evidenced the extreme importance of envelope-associated proteins in response to water stress.

**Functional analysis.** The functional enrichment analysis applied on the lists of differentially expressed proteins following PEG treatment, covering Gene Ontology and pathway annotations, was performed through the STRING database³². Among the enriched KEGG pathways within the blue class differentially expressed proteins, we denote the general metabolic pathways (map01100, FDR = 0.000147), ribosome (map03010, FDR = 0.000147) and sesquiterpenoid and triterpenoid biosynthesis (map00909, FDR = 0.0032). Among the most significantly enriched biological processes associated to the blue class of differentially expressed proteins, we found: nitrogen compound metabolic process (GO.0006807, FDR = 7.79e⁻⁰⁵), metabolic process (GO.0008152, FDR = 0.00022), protein metabolic process (GO.0019538, FDR = 0.00022), cellular nitrogen compound metabolic process (GO.0034641, FDR = 0.00022), cellular metabolic process (GO.0044237, FDR = 0.00022) and translation
<table><caption>**Table 1.** Summary of the most up-regulated polyethylene glycol (PEG)-induced proteins of *F. alni*.</caption><thead><tr><th>Response mechanism</th><th>Protein</th><th>Product</th><th>Tfold</th></tr></thead><tbody><tr><td rowspan="4">Membrane transport proteins</td><td>FRAAL2567</td><td>Putative autotransporter adhesion</td><td>5.67</td></tr><tr><td>FRAAL2549</td><td>Secondary metabolite secretion ATP-binding transport</td><td>4.80</td></tr><tr><td>FRAAL1420</td><td>Nitrate/sulfonate ATP-binding (ABC) transporter</td><td>2.89</td></tr><tr><td>FRAAL2085</td><td>Trk system potassium uptake protein TrkA</td><td>1.65</td></tr><tr><td rowspan="9">Intracellular osmolytes concentration</td><td>FRAAL0095</td><td>Osmotic-stress mechanosensitive channel MscL</td><td>4.30</td></tr><tr><td>FRAAL1888</td><td>Glutamine amidotransferase</td><td>2.20</td></tr><tr><td>FRAAL1197</td><td>2-deoxyribose-5-phosphate aldolase, NAD(P)-linked</td><td>2.15</td></tr><tr><td>FRAAL5155</td><td>Peptidase S51, dipeptidase E</td><td>1.91</td></tr><tr><td>FRAAL2547</td><td>Non-ribosomal peptide synthase</td><td>2.20</td></tr><tr><td>FRAAL2542</td><td>Non-ribosomal peptide synthase</td><td>1.80</td></tr><tr><td>FRAAL2545</td><td>Non-ribosomal peptide synthase</td><td>1.60</td></tr><tr><td>FRAAL2500</td><td>Polyketide cyclase</td><td>1.79</td></tr><tr><td rowspan="3">Molecular chaperones</td><td>FRAAL4431</td><td>Clp-family ATP-binding protease</td><td>4.00</td></tr><tr><td>FRAAL4781</td><td>Acyl-coenzyme A thioesterase</td><td>3.00</td></tr><tr><td>FRAAL6258</td><td>Acyl carrier protein (ACP)</td><td>2.33</td></tr><tr><td rowspan="5">Membrane lipid biosynthesis</td><td>FRAAL1744</td><td>Phospholipase D</td><td>1.82</td></tr><tr><td>FRAAL3585</td><td>Esterase</td><td>1.75</td></tr><tr><td>FRAAL5087</td><td>1-acylglycerol-3-phosphate O-acyltransferase</td><td>1.50</td></tr><tr><td>FRAAL2171</td><td>Acyl-CoA N-acyltransferase</td><td>1.50</td></tr><tr><td>FRAAL1222</td><td>Carbonic anhydrase</td><td>2.80</td></tr><tr><td rowspan="3">Detoxification of Reactive-Oxygen Species (ROS) and DNA repair</td><td>FRAAL6022</td><td>Gamma-aminobutyraldehyde dehydrogenase</td><td>2.17</td></tr><tr><td>FRAAL1417</td><td>Methylmalonate-semialdehyde dehydrogenase (acylating)</td><td>1.50</td></tr><tr><td>FRAAL5703</td><td>Global regulator (repressor) for SOS regulon LexA</td><td>1.56</td></tr></tbody></table>

<table><caption>**Table 2.** Summary of the most down-regulated polyethylene glycol (PEG)-modulated proteins of *F. alni*.</caption><thead><tr><th>Response mechanism</th><th>Protein</th><th>Product</th><th>Tfold</th></tr></thead><tbody><tr><td rowspan="10">Nitrogen fixation</td><td>FRAAL6814</td><td>Nitrogenase-associated homocitrate synthase NifV</td><td>-6.50</td></tr><tr><td>FRAAL6803</td><td>FeMo cofactor biosynthesis protein NifB</td><td>-5.24</td></tr><tr><td>FRAAL6810</td><td>Nitrogenase iron-molybdenum cofactor biosynthesis NifE</td><td>-3.80</td></tr><tr><td>FRAAL6809</td><td>Nitrogenase iron-molybdenum cofactor biosynthesis NifN</td><td>-3.50</td></tr><tr><td>FRAAL2491</td><td>squalene-hopenecyclase Shc2</td><td>-2.17</td></tr><tr><td>FRAAL6797</td><td>Ferredoxin</td><td>-2.00</td></tr><tr><td>FRAAL6799</td><td>2-oxoglutarate ferredoxinoxidoreductase beta-subunit KorA</td><td>-1.74</td></tr><tr><td>FRAAL1432</td><td>Squalene-hopenecyclase Shc1</td><td>-1.62</td></tr><tr><td>FRAAL4563</td><td>Transport protein associated with Fe-S cluster assembly SufB</td><td>-1.60</td></tr><tr><td rowspan="8">Respiration and energy production and conversion</td><td>FRAAL4147</td><td>Cytochrome c oxidase polypeptide I (AA3 subunit 1) CtaD</td><td>-3.70</td></tr><tr><td>FRAAL5625</td><td>Poly(3-hydroxybutyrate) depolymerasePha</td><td>-2.83</td></tr><tr><td>FRAAL6005</td><td>Biotin carboxyl carrier protein</td><td>-2.00</td></tr><tr><td>FRAAL6082</td><td>Succinate-semialdehyde dehydrogenase I, NADP-dependent</td><td>-2.00</td></tr><tr><td>FRAAL5481</td><td>Malate synthase G AceB</td><td>-1.75</td></tr><tr><td>FRAAL0315</td><td>Aldehyde dehydrogenase</td><td>-1.71</td></tr><tr><td>FRAAL1829</td><td>[NiFe] uptake hydrogenase, large subunit HupL</td><td>-1.70</td></tr><tr><td>FRAAL1039</td><td>NADH-quinoneoxidoreductase chain H NuoH</td><td>-1.50</td></tr><tr><td rowspan="2">DNA repair system</td><td>FRAAL6528</td><td>ATP-dependent protease, DNA repair protein</td><td>-2.36</td></tr><tr><td>FRAAL5811</td><td>ATP-dependent DNA helicase RecG</td><td>-1.50</td></tr></tbody></table>

(GO.0006412, FDR = 0.000706). Outputs of biological processes enrichment data were analysed with GOPlot³³ to generate insightful plots. Figure 2 and Supplementary Figures 1 and 2 highlight the most enriched biological processes associated to the blue differentially expressed, orange and red class proteins, respectively. There are no significantly enriched biological processes in the green class proteins.

Among the top blue class enriched biological processes (Fig. 2), we denote the nitrogen compound metabolic process, which includes proteins belonging to Nif family (e.g., FRAAL6811/NifK), Rpm family (e.g., FRAAL1926/RpmA), Rpl family (e.g., FRAAL5214/RplT) and other proteins like FRAAL1410/MtnD, FRAAL5703/LexA, FRAAL5734/RimO, FRAAL1085/RpsS, FRAAL1076/RpsL, FRAAL1447/UreC, FRAAL5238/PyrB, FRAAL5828/LeuD, FRAAL4973/HisL, FRAAL2146/Apt, FRAAL4574/Tal and FRAAL1197/DeoC.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Chord diagram showing the most enriched biological processes (GO terms) with their differentially expressed blue class proteins following PEG treatment. In each chord, enriched GO biological processes are shown on the right, and *F. alni* blue class differentially expressed proteins contributing to this enrichment are shown on the left. On the left side of the circle, each differentially expressed protein (locus tag (FRAAL)) is represented by a rectangle which colour is correlated to the value of the Tfold (logFC). PEG up-regulated proteins are displayed in red whereas down-regulated proteins are displayed in blue. Chords connect protein names with biological process GO term groups. Each GO term is represented by one coloured line.</figcaption></figure>

**F. alni** protein-protein interactions. Supplementary Figures 3, 4, 5 and 6 represent *F. alni* protein-protein interactions within the blue, orange, green and red class proteins, respectively. Green nodes colour in Supplementary Figure 3 denotes up-regulated proteins whereas red nodes colour denotes down-regulated proteins. For instance, the network corresponding to the blue class proteins (Supplementary Figure 3) highlights an interaction between an up-regulated putative Clp-family ATP-binding protease (FRAAL4431) and a down-regulated putative signal peptide (FRAAL4430).

**Cluster of Orthologous Genes (COGs).** Figure 3 and Supplementary Figures 7, 8 and 9 show the distribution of the blue, orange, green and red classes proteins into the COG functional categories, respectively. Green colour within the Fig. 3 denotes up-regulated proteins whereas red colour denotes down-regulated proteins. We were not able to assign several hypothetical proteins to any COG category. For the blue class (Fig. 3), five COG categories/classes contain only up-regulated proteins—D (cell cycle control, cell division, chromosome partitioning), F (nucleotide transport and metabolism), S (function unknown), U (intracellular trafficking, secretion and vesicular transport) and W (extracellular structures). These results provide parallel pieces of evidence (see also Table 1) in support of the importance of protein clusters related to cellular processes and signalling and metabolism as well as some poorly characterised proteins during growth of *F. alni* under water stress.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Predicted COG classes of up- and down-regulated *F. alni* blue class proteins. Up- and down-regulated COG classes in green and red colours, respectively. A: RNA processing and modification, B: Chromatin structure and dynamics, C: Energy production and conversion, D: Cell cycle control and mitosis, E: Amino Acid metabolism and transport, F: Nucleotide metabolism and transport, G: Carbohydrate metabolism and transport, H: Coenzyme metabolism, I: Lipid metabolism, J: Translation, K: Transcription, L: Replication and repair, M: Cell wall/membrane/envelope biogenesis, N: Cell motility, O: Post-translational modification, protein turnover, chaperone functions, P: Inorganic ion transport and metabolism, Q: Secondary Structure, R: General Functional Prediction only, S: Function Unknown, T: Signal Transduction, U: Intracellular trafficking and secretion, V: Defence mechanisms, W: Extracellular structures, Y: Nuclear structure, Z: Cytoskeleton, (-): No COG assigned.</figcaption></figure>

**Comparison between differentially expressed proteins following polyethylene glycol (PEG)-induced water stress and symbiosis.** We compared up- and down-regulated proteins related to tolerance to PEG-induced water stress and symbiosis in *F. alni*. Common up- (128) and down-regulated (37) proteins were identified and highlighted in Fig. 4, yielding evidence of common pathways between xerotolerance and symbiosis in *F. alni*.

**FRAAL0095-polyethylene glycol (PEG) interaction investigation.** Among the most up-regulated proteins, we chose the FRAAL0095 (Q0RUG3) as a case study. The amino acid sequence of this putative mech-anosensitive channel of large conductance (MscL) was submitted to the Swiss Model server³⁴. Templates search identified the structure of the MscL of *Mycobacterium tuberculosis* (PDBid: 2OAR) as the template with the high-est sequence identity (47.6%) with the target. It is crystallised as a homo-pentamer with three gold ions. This template was selected for model building with respect to its oligo-state. The resulting structure model shown in Supplementary Figure 10 exhibited satisfactory quality scores with a Global Model Quality Estimation (GMQE) score of 0.72. This metric, comprised between zero and one, reflects the expected accuracy of the model with respect to the target-template alignment³⁴. The higher the GMQE the more accurate is the model. Another quality assessment was performed through the Ramachandran diagrams (see Supplementary Figure 11) of the homo-pentamer model of FRAAL0095. It presented 11 residues (1.7%) within the outliers region, distributed on all five chains. They mainly occur at the top of the channel (extra-cellular part), or in the trans-membrane region. The global quality of the model was thus considered as satisfactory.

Docking of the PEGbasic and related probes (Fig. 5) was performed on this structure model. The resulting docking poses were merged and used by the SOM-BSfinder tuned version to define a 3D map of the preferential interaction spots of the PEG and its related probes with the protein surface. The resulting map and hot spots are shown in Fig. 6(A–F). Three populated regions, called Consensual Clusters (CCs), could be identified at equiv-alent positions between equivalent helices from pairs of adjacent chains of the homo-pentamer (Fig. 6(C)). A fourth CC, ranked as least important and coloured in red on Fig. 6(C), was detected in a higher position than the other CCs.

Moderate and long chain PEG molecules are polymers with the repeating unit -CH₂CH₂O-. Their length determines their molecular weights which are computed as follows: MW = 18.02 + 44.05*n*, where n in the num-ber of the repeating units. For instance, PEG400 is composed of 8–9 repeating units and PEG8000 is composed of 181–182 repeating units. PEG400 was herein used as a basic model for moderate to long chain PEGs to conduct molecular docking simulations. In fact, the moderate length of PEG400 can be handled by docking programs for accurate calculations while larger sizes (e.g., PEG8000) would be time consuming with poor quality (data not shown). We hypothesised that data collected with PEG400 along with the different probes will provide insightful results on the PEG-protein interactions for our case study. PEG400 (see Fig. 5) was thus docked on the structure model of the FRAAL0095. The SOM-BSfinder analysis returned a 3D map with a different shape as compared to the one obtained with the PEGbasic and related probes. It presented higher global density, which can be explained by the length and the flexibility of the polyethylene chain of PEG400. Still, a CC was detected between each two equivalent helices from adjacent chains (Fig. 6(F)), which is consistent with the results obtained with the PEG-basic and the probes. Moreover, an additional CC (coloured in green on Fig. 6(F)) was identified in a higher position than the other CCs. The docking results suggested a folding of the PEG400 inside the channel formed by the five chains. This is consistent with the literature³⁵,³⁶. In fact, it was proposed that moderate (PEG800) to long chain (PEG4700 and higher) PEGs would fold themselves and form a corona that mainly establishes hydrophobic
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Differentially expressed proteins following polyethylene glycol (PEG) treatment versus differentially expressed proteins involved in symbiosis. Proteomic data of *in vitro* PEG-treated cells on the x-axis with the symbiotic transcriptome⁶ (21 dpi) on the y-axis (scale log10). The symbiotic transcriptome data presented here are those with a FC above 2 and below 0.5 are represented together with those proteomic with FC values above 1 and below -1.</figcaption></figure>

interactions with proteins³⁵⁻³⁸. We further identified the residues that may interact with the PEG400 as those located at a maximal distance of 4.0 Å from any of its atoms; and obtained residues alanine 25 (A25), alanine 28 (A28), valine 29 (V29), threonine 31 (T31), valine 34 (V34), lysine 35 (K35), arginine 65 (R65), glutamate 67 (E67), asparagine 70 (N70), serine 71 (S71), alanine 74 (A74) and phenylalanine 75 (F75). These interacting residues were mostly hydrophobic (A, V, F) or polar with neutral side chains (T, S, N). On the protein surface, they constituted a series of symmetric hydrophobic pockets (see Supplementary Figure 10) that define what have previously been described as the hydrophobic core of MscL³⁹. Hydrophobic and polar residues appeared to be central to the PEG-FRAAL0095 interactions.

## Discussion

Little is known about the means by which *Frankia* cells tolerate desiccation in the soil¹⁰. In view of the fact that increasing concentrations of greenhouse gases in the atmosphere and climate change are expected to modify the global water cycle and to increase the odds of worsening drought in the next decades⁴⁰, an understanding of the responses of *F. alni* to desiccation and thus to osmotic stress becomes an important issue in the attempt to adapt symbiosis to desiccated niches. In the present work we have analysed the expression of proteins in the alder symbiont *F. alni* in response to PEG-induced water stress to characterise the impact of desiccation on nitrogen-fixing symbiotic bacteria. We used liquid medium since it easier to have reproducible conditions with healthy cells devoid of agar.

The results of the experiments summarised in Tables 1 and 2 indicated six most up-regulated and three most down-regulated PEG-modulated response mechanisms, respectively. The most remarkable differences in the analysed proteomes of *F. alni* cells grown with and without PEG could be assigned to envelope-associated proteins (e.g., FRAAL2567 (5.67) and FRAAL2549 (4.8); cf. Table 1). FRAAL2567 is a putative autotransporter adhesin, which is concordant with previous studies indicating the importance of proteins associated with the modulation of the structure and function of the three-dimensional extracellular matrix under stress conditions—desiccation, UV irradiation and oxidation. In addition, cell envelope remodelling, transcriptional and translational regulation and salt-stress responsive hypothetical proteins were shown to be important for salt-tolerant *Frankia* strains¹³. A previous study²⁰ that concords with present results (Table 1), showed that *Rhodococcus jostii* RHA1 desiccation-specific transcriptome had several genes associated with cellular processes including lipid metabolism and cell envelope modification. ABC transporters, like FRAAL2549, are integral membrane proteins that typically use cellular energy to translocate solutes across cellular membranes in all phyla⁴¹. Through the transport of molecules, ABC transporters are involved in diverse cellular processes including osmotic homeostasis and resistance to xenotoxins⁴². Moreover, up-regulation of a potassium transporter (FRAAL2085 (1.65)) was observed. Orthologs of this protein were reported to scavenge K⁺ ions from the environment and maintain cell hydrostatic pressure in response to desiccation. In addition, MscL, one of the best studied highly conserved mechanosensitive channels throughout the bacterial kingdom⁴³, was found to be highly up-regulated in *F. alni*
<figure><img src="image_7.png" /><figcaption>(i) PEGbasic</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>(ii) PEG400</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>(iii) Probe 0</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>(iv) Probe 1</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>(v) Probe 2</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>(vi) Probe 3</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>(vii) Probe 4</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Figure 6.** Surface mapping of the FRAAL0095 and identified hot-spots. Panels in the first line present results obtained with the docking of PEG basic and derived probes. Panels in the second line present results obtained with the docking of PEG 400. (A) A face-view of the protein structure along with the u-Matrix obtained with SOM-BSfinder analysis of the PEGbasic docking results, with a zoom on the u-Matrix. (B) A top-view of the protein structure along the u-Matrix of the PEGbasic docking results. (C) A top-view of the protein structure along with the CCs (hot-spots) identified using the SOM-BSfinder analysis of the PEGbasic docking results. (D) A face-view of the protein structure along with the u-Matrix obtained with SOM-BSfinder analysis of the PEG400 docking results, with a zoom on the u-Matrix. (E) A top-view of the protein structure along the u-Matrix of the PEG400 docking results. (F) A top-view of the protein structure along with the CCs (hot-spots) identified using the SOM-BSfinder analysis of the PEG400 docking results.</figcaption></figure>

following exposure to PEG (FRAAL0095 (4.30)). This transmembrane channel acts through switching from closed to open conformation under severe osmotic challenges by opening a water-filled pore⁴⁴,⁴⁵. Molecular modelling revealed potential interacting sites of the PEG molecule inside the channel formed by five monomers of MscL of *F. alni*. These sites were located between each pair of adjacent chains in a symmetric frame. Interactions with the PEG mainly included hydrophobic and polar residues that define the hydrophobic core that is *sine qua none* to the channel activity⁴⁶. Recently, the co-crystal structure of the C-terminal SH3 domain of myosin IB from *Entamoeba histolytica* with PEG was reported⁴⁷. Gautam and colleagues⁴⁷ brought insightful elements for a better understanding of the interactions of PEG with protein surfaces and reported in particular the corona-like binding mode previously described⁴⁸⁻⁵¹. Our findings are consistent with the literature, which adds confidence to our computational approach. In fact, long-chain PEGs are flexible polymers with multiple rotatable bonds and high degrees of freedom, which makes them a challenging case for molecular docking software. In the present work, we established a protocol to circumvent this issue by using the repeating unit of the PEG and its chemical features
as probes in order to infer a simplistic model for docking. Then, the SOM-based method permitted to generate an accurate map of the binding hot spots on the MscL surface, along with a prediction of the overall binding mode. Our results suggest that long-chain PEGs may act as an inner coat of the channel<sup>39,52</sup>. This mode of binding may either induce an increase in water molecules conductance through the channel or induce a mechanistic blocking of its entrance. Our findings constitute a response element to the role of MscL (FRAAL0095) in the *F. alni* desic- come, and a starting point to conduct further investigations on the PEG-FRAAL0095 interactions.

As expected, envelope stress would induce a cascade of other defence pathways. Most notable is the induction of one CRISPR-associated proteins (FRAAL0457 (1.67) being part of the Green class) supporting published obser- vations that expression of CRISPR-Cas components in several bacterial species can be induced following envelope stress<sup>53</sup>. Supplementary Table S3 lists the CRISPR-associated PEG-induced proteins of *F. alni*. Also, the present analysis revealed the up-regulation of proteins associated with cell osmoprotection. For instance, FRAAL1888 (2.20) and FRAAL2171 (1.5) are involved in the synthesis of the dipeptide N-acetylglutaminylglutamine amide (NAGGN) that plays a key role in osmoprotection<sup>54</sup>. In addition to its role in NAGGN biosynthesis, FRAAL1888 is involved in the conversion of glutamine into glutamate, one of the most abundant and up-regulated amino acids in *Alnus* and *Casuarina* nodules, for which a role as osmolyte was postulated by Brooks and Benson<sup>55</sup>. Also, the FRAAL5155 peptidase (1.91), by the cleavage of dipeptides into single peptides, would provide needed intra- cellular osmolytes to offset the high extracellular osmotic pressure.

In addition, consistent with previous studies, we observed a significant up-regulation of FRAAL4431 (4) encoding a Clp protease. Indeed, ClpP proteases play an essential role in removing damaged proteins from low GC Gram-positive bacteria under stress conditions such as osmotic shock<sup>56</sup>. Furthermore, proteolysis of such proteins would release amino acids that would in turn provide osmolytes to help compensate the osmotic shock. Similarly, a Clp-protease was induced during desiccation and dehydration in *Pyropia columbina*<sup>57</sup>. FRAAL1222 encoding a carbonic anhydrase and presumably involved in stress response was quantified in larger amount (2.8) in PEG-grown *F. alni* cells. It is expected that this enzyme, by hydrating CO₂ to carbonate (H₂CO₃), would increase the concentration of protons within the cell and decrease the pH. A previous study<sup>18</sup> showed an increase in the expression of an orphan carbonic anhydrase in the cyanobacterium *Microcoleus vaginatus* using a whole genome transcriptional time course assay in response to hydration. The observed increase in several ROS-scavenging proteins and translation-associated ribosome proteins provided evidence in support of a higher demand of detoxification and protein synthesis during growth under PEG-induced water stress conditions. For example, FRAAL6022 (2.17) encodes a gamma-aminobutyraldehyde dehydrogenase (GABA synthesis). GABA can act as a hydroxyl radical scavenger<sup>58</sup>.

*Deinococcus radiodurans*' xerotolerance was explained by the phenomenon of gene sharing⁵⁹: during land col- onization by *Terrabacteria*, gene products (e.g., DNA repair genes, etc.) contributing to this bacterium's ionizing radiation were recruited to serve an additional function that is desiccation tolerance⁶⁰. In *F. alni*, we observed an upregulation of the SOS response genes transcriptional repressor LexA (1.56) with an expected downreg- ulation of some DNA repair proteins (Table 2). Furthermore, the investigation of the relationship between desiccation tolerance and nitrogen metabolism in *Anabaena* sp. PCC 7120 suggested that terrestrial cyanobac- teria may acclimate to desiccation stress *via* nitrogen (N₂) fixation by using desiccation inducible genes⁶¹. In our study, TFold analysis pinpointed three classes of differentially down-regulated proteins (e.g., FRAAL6814 (−6.5), FRAAL4147 (−3.7) and FRAAL6528 (−2.36); cf. Table 2). The remarkable decrease in protein expres- sion related to nitrogen fixation (more than sixfold), aerobic respiration (more than threefold) and homologous recombination (more than twofold) quantitatively reflects the magnitude of the impact of nitrogen-, oxygen- and homologous-recombination-associated energy demands in PEG-grown *F. alni* cells.

Identified common pathways between symbiosis⁶¹ and desiccation tolerance (Fig. 4) are in accordance with the fact that plant-growth-promoting (PGP) rhizobacteria augment plant tolerance to drought⁶². In line with this, it was recently demonstrated that bacterial endophytic communities promote date palm (*Phoenix dactylifera* L.) growth under drought conditions⁶³. Also, it was shown that plant growth promotion ability exerted by bacteria is a drought-induced effect⁶⁴. In addition, a previous work²⁵ showed an induction of the threonine dehydratase (DR_0567)—converting threonine into NH₃/NH₄⁺—and the nodulation efficiency protein D (NfeD, DR_2142 (WP_010888773.1)) following desiccation of *D. radiodurans* R₁.

Penetration into host tissues implies many stresses, one of which is osmotic since *Alnus glutinosa*, as most alder species, thrives on river banks where water has low salt levels, while root tissues have a higher, isotonic osmotic potential. Furthermore, alder synthesises many peptides upon entry of *Frankia*, many of which are defen- sins that bind to the cell membrane and modify its porosity⁶⁵ and thus the ability to cope with osmolytes. Rhizobia have been shown to have many determinants associated with desiccation tolerance, many of which have similar- ities to those seen upregulated in *Frankia* such as transporters⁶⁶, the mechanosensitive channel⁶⁷ or disaccharide accumulation⁶⁸. However the link between osmotic response and symbiosis establishment was not very strong, implying the responses seen here would have more usefulness for saprotrophic soil existence.

## Conclusions

The present paper is the first high-throughput proteomic study of *F. alni* subjected to water-related stress. This research, analysing the most representative proteome of a *Frankia* strain, contributes to a better comprehension of environmental stress adaptation particularly in desiccated soils. In the future, our data might be used in further comparative proteogenomic studies of nitrogen-fixing plant symbionts.

## Materials and Methods

**Bacterial growth and treatment with polyethylene glycol (PEG).** *F. alni* strain ACN14a⁶⁹ was grown in liquid BAP- medium without ammonium as described earlier⁷⁰. The cells were then syringed thoroughly with needles of decreasing width (21G–27G), the cells OD₆₀₀ measured and the suspension diluted prior to inoculation.
An OD₆₀₀ inoculation of 0.1 was used to make a growth dynamic and to monitor ARA activity at 0.3, 0.9 and 2.8% PEG8000 (Sigma). Inoculations were then made at the same OD₆₀₀ in fresh BAP- medium without or with PEG8000 (Sigma) added to a final concentration of 0.9% w.vol⁻¹. The cells were grown in 125 mL Erlenmeyers with 40 mL medium at 28°C on an orbital shaker at 200 rpm for 7 days. The cells were harvested, their OD₆₀₀ measured, an aliquot inoculated onto LB agar plates to detect contaminations and the cells observed under the microscope. The cells were then sedimented (1,500 × g for 5 minutes) and the pellets frozen until proteomics analysis.

**Proteomic sample preparation and nanoscale liquid chromatography coupled to tandem mass spectrometry (nano LC-MS/MS) analysis of tryptic peptides.** The protein content of cells grown without PEG and in the presence of PEG was established by a shotgun procedure. Following trypsin proteolysis, peptides were analysed by nanoLC-MS/MS with a Q-Exactive HF tandem mass spectrometer incorporating an ultra-high field Orbitrap analyser. Five independent biological replicates were analysed for both conditions: PEG-treated bacteria and untreated bacteria. A volume of 5 µL of LDS (1X) was added per mg of bacterial pellet prior a 5 min heat treatment at 99°C followed by a treatment for 5 min in an ultrasonic bath. Cells and debris were transferred into a 2 mL Precellys (Bertin Technologies, F-78180 Montigny le Bretonneux) tube containing 200 mg of glass beads and the samples were subjected to 3 cycles of grinding at 6,500 rpm for 20 sec by means of a Precellys grinder (Bertin technologies). Samples were then centrifuged for 40 sec at 16,000 g. The resulting supernatants were incubated for 10 min at 99 °C and were subjected to a short SDS-PAGE migration as previously described⁷¹. The polyacrylamide bands containing the whole solubilised protein content of each sample were processed as previously described⁷². Briefly, they were subjected to DTT reduction, treated with iodoacetamide and then, proteolysed with Sequencing Grade Trypsin (Roche, F-38240 Meylan) using 0.01% of proteaseMAX detergent (Promega, F-69260 Charbonnières-les-Bains). The resulting peptides (10 µL) were analysed in data-dependent mode with an ESI-Q Exactive HF mass spectrometer (ThermoFisher Scientific, F-91963, Courtaboeuf) equipped with an ultra-high field Orbitrap analyser and coupled to an Ultimate 3000 176 RSL Nano LC System (ThermoFisher). Peptides were injected onto a reverse phase Acclaim PepMap 100 C18 column (3 µm, 100 Å, 75 µm id × 500 mm) and resolved at a flow rate of 0.2 µL/min with a 60 min gradient of CH3CN in presence of 0.1% formic acid. A dataset of 10 nanoLC-MS/MS runs were recorded. The Q-Exactive HF instrument was operated with Top20 standard parameters and a dynamic exclusion of 10 sec as previously described⁷³. MS/MS spectra were searched using MASCOT 2.2.04 software (Matrix Science, London, W1U 7GB, UK) against the *F. alni* database (5,804 protein sequences) with the following parameters: full-trypsin specificity, maximum of two missed cleavages, mass tolerances of 5 ppm on the parent ion and 0.02 Da on the higher energy collisional dissociation-induced peptidic fragments, fixed modification of carboxyamidomethylated cysteine (+57.0215), and oxidised methionine (+15.9949) as dynamic modifications. Peptide-to-MS/MS spectrum matching with a MASCOT score below a *p*-value of 0.05 were selected and assigned to unique peptide sequence following the parsimony principle. A protein was considered valid when at least two different peptides were detected. The false-positive rate for protein identification was estimated by a search with a reverse decoy database to be below 0.1% using the same parameters. Proteins were quantified based on their spectral counts. The normalised spectral abundance factor (NSAF) was calculated by dividing the spectral count for each observed protein by the polypeptide theoretical mass, as described previously³⁰ and is presented as percentage of the NSAF sum considering all proteins. Proteome comparison between both conditions was done taking into account the five biological replicates with the TFold module from the PatternLab software and standard normalisation, as previously described⁷¹.

**Homology-based functional analysis and Clusters of Orthologous Groups of proteins (COGs) prediction.** Protein sequences of the PEG-responding differentially expressed proteins were retrieved from the National Center for Biotechnology Information (NCBI) database⁷⁴ and uploaded into the STRING database³² analysis tool for Gene Ontologies and metabolic pathways functional enrichments. Biological processes and metabolic pathways with False Discovery Rate (FDR) ≤ 0.005 were considered as significant. *F. alni* Clusters of Orthologous Groups of proteins (COGs) were inferred for an important set of proteins based on BLASTP⁷⁵ best reciprocal hits with *Frankia* sp. EAN1pec for which these clusters are known and available in the database of COGs⁷⁶.

**Protein-protein interaction network.** The sequences of the 294 proteins of *F. alni* belonging to the PatternLab blue class (for which identifications satisfied both the fold (>1.5) and statistical criteria (*p*-value < 0.05)) were downloaded and saved in a multi-FASTA file. This file was then uploaded in the STRING database³² and *F. alni* was chosen as a query microorganism. A total of 294 hits with 100% identity were detected and their protein-protein interactions were predicted. In order to generate protein-protein interactions integrating the expression data information of each protein, data were integrated into the Cytoscape tool⁷⁷. Expression values were imported as node attributes. The same steps were followed to analyse protein-protein interaction networks of the other detected proteins (belonging to the PatternLab orange, green and red classes (see Results section for details)). All single nodes were not represented with Cytoscape which considers only interacting proteins.

**Re-annotation of differentially expressed hypothetical proteins.** Differentially expressed proteins that were annotated as hypothetical (or conserved hypothetical) were re-annotated. For this, we conducted new analyses using the Position-Specific Iterated BLAST (PSI-BLAST) program⁷⁸ at http://blast.ncbi.nlm.nih.gov/ with default parameters against a database of non-redundant protein sequences (nr).

**Molecular modelling and docking.** *Comparative modelling.* Protein sequence of the FRAAL0095, one of the most up-regulated *Frankia* genes (4.30), identified as coding for an osmotic-stress large conductance mechanosensitive channel, was submitted to the Swiss Model Server³⁴. Search for templates of experimentally
determined protein structures within the Swiss-Model Template Library (SMTL) was performed. Identified templates, ranked according to their sequence identity rate (IR) and their coverage rate (CR) with the target sequence, were examined. The template(s) with the highest IR and CR were selected for model building. The model with the most satisfactory quality assessment scores was considered for the subsequent analysis. Ramachandran diagrams of the selected model were generated using the RAMPAGE server⁷⁹, and used to assess the 3D model quality.

**Molecular docking and surface mapping.** Molecular docking of the PEG and related molecules (herein called ligands) on the structure model of the protein target (herein called receptor) was performed using AutoDock Vina 1.1.2⁸⁰. The latter program requires input files of the ligands and the receptor in PDBQT (Protein Data Bank⁸¹, Partial Charge ('Q'), and Atom Type ('T')) format. Default parameters were used and a maximum of 20 lowest-energy poses were kept for each ligand. Two forms of the PEG were considered: (i) a short form having the formula $C_2H_{10}O_3$ and is referred to as PEGbasic and (ii) a longer form having the formula $C_{18}H_{38}O_{10}$ to which we refer as PEG400. Additionally, five chemical substructures of the PEG were generated and used as probes to map hot spots on the protein surface. Simplified Molecular Input Line Entry System (SMILES) conformations of PEGbasic, PEG400 and the probes were generated manually, then converted to PDBQT format using the OpenBabel package (http://openbabel.org)⁸². The PDB file of the target was used to generate the corresponding PDBQT format by adding hydrogen atoms and atomic partial charges using the OpenBabel package⁸². Each ligand was docked 500 times on the protein. Docking results were analysed using a customised version of the SOM-BSfinder method⁸³. Herein, input data consisted in atomic coordinates of the docked PEGbasic, PEG400, and the five probes. Otherwise, default parameters were used as described previously⁸³. PyMol 1.8.2.1 (Schrodinger, LLC)⁸⁴ was used to visualise the results and generate the figures.

**Proteomics data repository.** The mass spectrometry proteomic data were deposited at the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository with the data set identifiers PXD007226 and DOI 10.6019/PXD007226.

## References

1. Normand, P., Benson, D. R., Berry, A. M. & Tisa, L. S. In *The Prokaryotes: Actinobacteria* (eds Eugene Rosenberg et al.) 339–356 (Springer Berlin Heidelberg, 2014).

2. Benson, D. R. & Silvester, W. B. Biology of *Frankia* strains, actinomycete symbionts of actinorhizal plants. *Microbiol Rev* 57, 293–319 (1993).

3. Sen, A. et al. Phylogeny of the class Actinobacteria revisited in the light of complete genomes. The orders 'Frankiales' and Micrococcales should be split into coherent entities: proposal of *Frankiales* ord. nov., *Geodermatophilales* ord. nov., *Acidothermales* ord. nov. and *Nakamurellales* ord. nov. *Int J Syst Evol Microbiol* 64, 3821–3832, https://doi.org/10.1099/ijss.0.063966-0 (2014).

4. Nouioui, I. et al. Proposal of a type strain for *Frankia alni* (Woronin 1866) Von Tubeuf 1895, emended description of *Frankia alni*, and recognition of *Frankia casuarinae* sp. nov. and *Frankia elaeagni* sp. nov. *Int J Syst Evol Microbiol* 66, 5201–5210, https://doi.org/10.1099/ijsem.0.001496 (2016).

5. Normand, P. et al. Genome characteristics of facultatively symbiotic *Frankia* sp. strains reflect host range and host plant biogeography. *Genome Res* 17, 7–15, https://doi.org/10.1101/gr.5798407 (2007).

6. Alloisio, N., Kucho, K.-i., Pujic, P. & Normand, P. In *Biological Nitrogen Fixation* 757–768 (John Wiley & Sons, Inc, 2015).

7. Mastronunzio, J. E., Huang, Y. & Benson, D. R. Diminished exoproteome of *Frankia* spp. in culture and symbiosis. *Appl Environ Microbiol* 75, 6721–6728, https://doi.org/10.1128/aem.01559-09 (2009).

8. Lebre, P. H., De Maayer, P. & Cowan, D. A. Xerotolerant bacteria: surviving through a dry spell. *Nat Rev Micro* 15, 285–296, https://doi.org/10.1038/nrmicro.2017.16 (2017).

9. Kakumanu, M. L., Cantrell, C. L. & Williams, M. A. Microbial community response to varying magnitudes of desiccation in soil: a test of the osmolyte accumulation hypothesis. *Soil Biol Biochem* 57, 644–653 (2013).

10. Burleigh, S. H. & Dawson, J. O. Desiccation tolerance and trehalose production in *Frankia* hyphae. *Soil Biol Biochem* 26, 593–598 (1994).

11. Selvakesavan, R. K. et al. Intraspecies variation in sodium partitioning, potassium and proline accumulation under salt stress in *Casuarina equisetifolia* Forst. *Symbiosis* 70, 117–127, https://doi.org/10.1007/s13199-016-0424-9 (2016).

12. Srivastava, Å., Singh, A., Singh, S. S. & Mishra, A. K. Salt stress-induced changes in antioxidative defense system and proteome profiles of salt-tolerant and sensitive *Frankia* strains. *J Environ Sci Health. Part A, Toxic/hazardous substances & environmental engineering* 52, 420–428, https://doi.org/10.1080/10934529.2016.1270672 (2017).

13. Oshone, R. et al. Genomic, transcriptomic, and proteomic approaches towards understanding the molecular mechanisms of salt tolerance in *Frankia* strains isolated from *Casuarina* trees. *BMC Genom* 18, 633, https://doi.org/10.1186/s12864-017-4056-0 (2017).

14. Dawson, J. O. & Gibson, A. H. Sensitivity of selected *Frankia* isolates from *Casuarina*, *Allocasuarina* and North American host plants to sodium chloride. *Physiol Plant* 70, 272–278, https://doi.org/10.1111/j.1399-3054.1987.tb06144.x (1987).

15. Oshone, R., Mansour, S. R. & Tisa, L. S. Effect of salt stress on the physiology of *Frankia* sp strain Cc16. *J Biosci* 38, 699–702 (2013).

16. Sghaier, H. et al. Stone-dwelling actinobacteria *Blastococcus saxobidens*, *Modestobacter marinus* and *Geodermatophilus obscurus* proteogenomes. *The ISME J* 10, 21–29, https://doi.org/10.1038/ismej.2015.108 (2016).

17. Umezawa, T., Mizuno, K. & Fujimura, T. Discrimination of genes expressed in response to the ionic or osmotic effect of salt stress in soybean with cDNA-AFLP. *Plant, Cell Environ* 25, 1617–1625 (2002).

18. Rajeev, L. et al. Dynamic cyanobacterial response to hydration and dehydration in a desert biological soil crust. *ISME J* 7, 2178–2191, https://doi.org/10.1038/ismej.2013.83 (2013).

19. Mohammadipanah, F. & Wink, J. *Actinobacteria* from arid and desert habitats: Diversity and biological activity. *Front Microbiol* 6, 1541, https://doi.org/10.3389/fmicb.2015.01541 (2015).

20. LeBlanc, J. C., Gonçalves, E. R. & Mohn, W. W. Global response to desiccation stress in the soil actinomycete *Rhodococcus jostii* RHA1. *Appl Environ Microbiol* 74, 2627–2636, https://doi.org/10.1128/aem.02711-07 (2008).

21. Fredrickson, J. K. et al. Protein oxidation: key to bacterial desiccation resistance? *ISME J* 2, 393–403 (2008).

22. McIntyre, H. J. et al. Trehalose biosynthesis in *Rhizobium leguminosarum* bv. trifolii and its role in desiccation tolerance. *Appl Environ Microbiol* 73, 3984–3992, https://doi.org/10.1128/aem.00412-07 (2007).

23. Boncompagni, E., Osteras, M., Poggi, M. C. & le Rudulier, D. Occurrence of choline and glycine betaine uptake and metabolism in the family *rhizobiaceae* and their roles in osmoprotection. *Appl Environ Microbiol* 65, 2072–2077 (1999).

24. De Philippis, R. & Vincenzini, M. Exocellular polysaccharides from cyanobacteria and their possible applications. *FEMS Microbiol Rev* 22, 151–175, https://doi.org/10.1111/j.1574-6976.1998.tb00365.x (1998).
25. Tanaka, M. et al. Analysis of *Deinococcus radiodurans*'s transcriptional response to ionizing radiation and desiccation reveals novel proteins that contribute to extreme radioresistance. *Genetics* **168**, 21–33, https://doi.org/10.1534/genetics.104.029249 (2004).

26. Cytryn, E. J. *et al.* Transcriptional and physiological responses of *Bradyrhizobium japonicum* to desiccation-induced stress. *J Bacteriol* **189**, 6751–6762, https://doi.org/10.1128/jb.00533-07 (2007).

27. Jeon, J. M., Lee, H. I., Sadowsky, M. J., Sugawara, M. & Chang, W. S. Characterization of a functional role of the *Bradyrhizobium japonicum* Isocitrate Lyase in desiccation tolerance. *Int J Molec Sci* **16**, 16695–16709, https://doi.org/10.3390/ijms160716695 (2015).

28. Potts, M., Slaughter, S. M., Hunneke, F. U., Garst, J. F. & Helm, R. F. Desiccation tolerance of prokaryotes: application of principles to human cells. *Integ Comp Biol* **45**, 800–809, https://doi.org/10.1093/icb/45.5.800 (2005).

29. Leprince, O. & Buitink, J. Desiccation tolerance: from genomics to the field. *Plant Sci* **179**, 554–564, https://doi.org/10.1016/j.plantsci.2010.02.011 (2010).

30. Christie-Oleza, J. A., Fernandez, B., Nogales, B., Bosch, R. & Armengaud, J. Proteomic insights into the lifestyle of an environmentally relevant marine bacterium. *ISME J* **6**, 124–135, https://doi.org/10.1038/ismej.2011.86 (2012).

31. Rubiano-Labrador, C. et al. Proteogenomic insights into salt tolerance by a halotolerant alpha-proteobacterium isolated from an Andean saline spring. *J Proteom* **97**, 36–47, https://doi.org/10.1016/j.jprot.2013.05.020 (2014).

32. Szklarczyk, D. et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. *Nucl Acids Res* **45**, D362–D368, https://doi.org/10.1093/nar/gkw937 (2017).

33. Walter, W., Sanchez-Cabo, F. & Ricote, M. GOplot: an R package for visually combining expression data with functional analysis. *Bioinformatics* **31**, 2912–2914, https://doi.org/10.1093/bioinformatics/btv300 (2015).

34. Bienert, S. et al. The SWISS-MODEL Repository-new features and functionality. *Nucl Acids Res* **45**, D313–D319, https://doi.org/10.1093/nar/gkw1132 (2017).

35. Hašek, J. Poly(ethylene glycol) interactions with proteins. *Zeitschrift für Kristall Suppl* **23**, 613–618 (2006).

36. Wu, J. et al. Binding characteristics between polyethylene glycol (PEG) and proteins in aqueous solution. *J Mat Chem B* **2**, 2983–2992, https://doi.org/10.1039/C4TB00253A (2014).

37. Knowles, D. B. et al. Chemical interactions of polyethylene glycols (PEGs) and glycerol with protein functional groups: applications to effects of PEG and glycerol on protein processes. *Biochemistry* **54**, 3528–3542, https://doi.org/10.1021/acs.biochem.5b00246 (2015).

38. Lim, K. & Herron, J. N. In *Poly(Ethylene Glycol) chemistry: biotechnical and biomedical applications* (ed. Milton Harris, J.) 29–56 (Springer US, 1992).

39. Perozo, E., Cortes, D., Sompornpisut, P., Kloda, A. & Martinac, B. Open channel structure of MscL and the gating mechanism of mechanosensitive channels. *Nature* **418**, 942–948 (2002).

40. Prudhomme, C. et al. Hydrological droughts in the 21st century, hotspots and uncertainties from a global multimodel ensemble experiment. *Proc Nat Acad Sci USA* **111**, 3262–3267, https://doi.org/10.1073/pnas.1222473110 (2014).

41. Jones, P. M. & George, A. M. The ABC transporter structure and mechanism: perspectives on recent research. *Cell Molec Life Sci*: CMLS **61**, 682–699, https://doi.org/10.1007/s00018-003-3336-9 (2008).

42. Martinoia, E. et al. Multifunctionality of plant ABC transporters—more than just detoxifiers. *Planta* **214**, 345–355 (2002).

43. Zhong, D. & Blount, P. Phosphatidylinositol is crucial for the mechanosensitivity of *Mycobacterium tuberculosis* MscL. *Biochemistry* **52**, 5415–5420, https://doi.org/10.1021/bi400790J (2013).

44. Sukharev, S., Blount, P., Martinac, B. & Kung, C. Mechanosensitive channels of *Escherichia coli*: the MscL gene, protein, and activities. *Annu Rev Physiol* **59**, 633–657 (1997).

45. Hamill, O. & Martinac, B. Molecular basis of mechanotransduction in living cells. *Physiol Rev* **81**, 685–740 (2001).

46. Häse, C., Le Dain, A. & Martinac, B. Molecular dissection of the large mechanosensitive ion channel (MscL) of *E. coli*: mutants with altered channel gating and pressure sensitivity. *J Membr Biol* **157**, 17–25 (1997).

47. Gautam, G., Rehman, S., Pandey, P. & Gourinath, S. Crystal structure of the PEG-bound SH3 domain of myosin IB from *Entamoeba histolytica* reveals its mode of ligand recognition. *Acta Crystallogr D Struct Biol* **73**, 672–682, https://doi.org/10.1107/S2059798317009639 (2017).

48. Lim, S. T., Jane, J. L., Rajagopalan, S. & Seib, P. A. Effect of Starch Granule Size on Physical Properties of Starch-Filled Polyethylene Film. *Biotechnol Prog* **8**, 51–57 (1992).

49. Hašek, J. Poly (ethylene glycol) interactions with proteins. *Z. Kristallogr. Suppl* **23**, 613–618 (2006).

50. Wu, J. et al. Binding characteristics between polyethylene glycol (PEG) and proteins in aqueous solution. *J Mat Chem B* **2**, 2983–2992 (2014).

51. Knowles, D. B. et al. Chemical Interactions of Polyethylene Glycols (PEGs) and Glycerol with Protein Functional Groups: Applications to Effects of PEG and Glycerol on Protein Processes. *Biochemistry* **54**, 3528–3542 (2015).

52. Haswell, E., Phillips, R. & Rees, D. Mechanosensitive channels: what can they do and how do they do it? *Structure* **19**, 1356–1369, https://doi.org/10.1016/j.str.2011.09.005. (2011).

53. Sampson, T. R. & Weiss, D. S. CRISPR-Cas systems: new players in gene regulation and bacterial physiology. *Front Cell Inf Microbiol* **4**, 37, https://doi.org/10.3389/fcimb.2014.00037 (2014).

54. Sagot, B. et al. Osmotically induced synthesis of the dipeptide N-acetylglutaminylglutamine amide is mediated by a new pathway conserved among bacteria. *Proc Nat Acad Sci USA* **107**, 12652–12657, https://doi.org/10.1073/pnas.1003063107 (2010).

55. Brooks, J. M. & Benson, D. R. Comparative metabolomics of root nodules infected with *Frankia* sp. strains and uninfected roots from *Alnus glutinosa* and *Casuarina cunninghamiana* reflects physiological integration. *Symbiosis* **70**, 87–96, https://doi.org/10.1007/S13199-016-0379-x (2016).

56. Frees, D., Savijoki, K., Varmanen, P. & Ingmer, H. Clp ATPases and ClpP proteolytic complexes regulate vital biological processes in low GC, Gram-positive bacteria. *Molec Microbiol* **63**, 1285–1295, https://doi.org/10.1111/j.1365-2958.2007.05598.x (2007).

57. Contreras-Porcia, L., Thomas, D., Flores, V. & Correa, J. A. Tolerance to oxidative stress induced by desiccation in *Porphyra columbia* (Bangiales, Rhodophyta). *J Exp Bot* **62**, 1815–1829, https://doi.org/10.1093/jxb/erq364 (2011).

58. Smirnoff, N. & Cumbes, Q. J. Hydroxyl radical scavenging activity of compatible solutes. *Phytochemistry* **28**, 1057–1060, https://doi.org/10.1016/0031-9422(89)80182-7 (1989).

59. Sghaier, H., Mitomo, H. & Narumi, I. Genomic confession of *Deinococcus radiodurans*: it started out as a radiation-resistant organism. *Viva Origino* **33**, 243–257 (2005).

60. Sghaier, H., Narumi, I., Satoh, K., Ohba, H. & Mitomo, H. Problems with the current deinococcal hypothesis: an alternative theory. *Theory Biosci* **126**, 43–45, https://doi.org/10.1007/s12064-007-0004-x (2007).

61. Katoh, H. Desiccation-inducible genes are related to N(2)-fixing system under desiccation in a terrestrial cyanobacterium. *Biochim Biophys Acta* **1817**, 1263–1269, https://doi.org/10.1016/j.bbabio.2012.03.029 (2012).

62. Yang, J., Kloepper, J. W. & Ryu, C. M. Rhizosphere bacteria help plants tolerate abiotic stress. *Trends Plant Sci* **14**, 1–4, https://doi.org/10.1016/j.tplants.2008.10.004 (2009).

63. Cherif, H. et al. Oasis desert farming selects environment-specific date palm root endophytic communities and cultivable bacteria that promote resistance to drought. *Environ Microbiol Reports* **7**, 668–678, https://doi.org/10.1111/1758-2229.12304 (2015).

64. Rolli, E. et al. Improved plant resistance to drought is promoted by the root-associated microbiome as a water stress-dependent trait. *Environ Microbiol* **17**, 316–331, https://doi.org/10.1111/1462-2920.12439 (2015).

65. Carro, L. et al. *Alnus* peptides modify membrane porosity and induce the release of nitrogen-rich metabolites from nitrogen-fixing *Frankia*. *ISME J* **9**, 1723–1733, https://doi.org/10.1038/ismej.2014.257 (2015).
66. Vanderlinde, E. M. et al. Identification of a novel ABC transporter required for desiccation tolerance, and biofilm formation in *Rhizobium leguminosarum* bv. *vioiae* 3841. *FEMS Microbiol Ecol* **71**, 327–340, https://doi.org/10.1111/j.1574-6941.2009.00824.x (2010).

67. Balleza, D., Gomez-Lagunas, F. & Quinto, C. Cloning and functional expression of an MscL ortholog from *Rhizobium etli*: characterization of a mechanosensitive channel. *J Membrane Biol* **234**, 13–27, https://doi.org/10.1007/s00232-010-9235-8 (2010).

68. Gouffi, K., Pica, N., Pichereau, V. & Blanco, C. Disaccharides as a new class of nonaccumulated osmoprotectants for *Sinorhizobium meliloti*. *Appl Environ Microbiol* **65**, 1491–1500 (1999).

69. Normand, P. & Lalonde, M. Evaluation of *Frankia* strains isolated from provenances of two *Alnus* species. *Can J Microbiol* **28**, 1133–1142, https://doi.org/10.1139/m82-168 (1982).

70. Carro, L. et al. Physiological effects of major up-regulated *Alnus glutinosa* peptides on *Frankia* sp. ACN14a. *Microbiology* **162**, 1173–1184, https://doi.org/10.1099/mic.0.000291 (2016).

71. Hartmann, E. M., Allain, F., Gaillard, J. C., Pible, O. & Armengaud, J. Taking the shortcut for high-throughput shotgun proteomic analysis of bacteria. *Meth Molec Biol* **1197**, 275–285, https://doi.org/10.1007/978-1-4939-1261-2_16 (2014).

72. Clair, G., Armengaud, J. & Duport, C. Restricting fermentative potential by proteome remodeling: an adaptive strategy evidenced in *Bacillus cereus*. *Molec Cell Proteom: MCP* **11**, M111 013102, https://doi.org/10.1074/mcp.M111.013102 (2012).

73. Klein, G. et al. RNA-binding proteins are a major target of silica nanoparticles in cell extracts. *Nanotoxicology* **10**, 1555–1564, https://doi.org/10.1080/17435390.2016.1244299 (2016).

74. Database Resources of the National Center for Biotechnology Information. *Nucl Acids Res* **45**, D12–D17, https://doi.org/10.1093/nar/gkw1071 (2017).

75. Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. Basic local alignment search tool. *J Molec Biol* **215**, 403–410, https://doi.org/10.1016/s0022-2836(05)80360-2 (1990).

76. Tatusov, R. L., Galperin, M. Y., Natale, D. A. & Koonin, E. V. The COG database: a tool for genome-scale analysis of protein functions and evolution. *Nucl Acids Res* **28**, 33–36 (2000).

77. Shannon, P. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. *Genome Res* **13**, 2498–2504, https://doi.org/10.1101/gr.1239303 (2003).

78. Altschul, S. F. et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. *Nucl Acids Res* **25**, 3389–3402 (1997).

79. Lovell, S. C. et al. Structure validation by Calpha geometry: phi, psi and Cbeta deviation. *Proteins* **50**, 437–450, https://doi.org/10.1002/prot.10286 (2003).

80. Trott, O. & Olson, A. J. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. *J Comput Chem* **31**, 455–461, https://doi.org/10.1002/jcc.21334 (2010).

81. Rose, P. W. et al. The RCSB protein data bank: integrative view of protein, gene and 3D structural information. *Nucl Acids Res* **45**, D271–D281, https://doi.org/10.1093/nar/gkw1000 (2017).

82. O'Boyle, N. M. et al. Open Babel: an open chemical toolbox. *J Cheminform* **3**, 33, https://doi.org/10.1186/1758-2946-3-33 (2011).

83. Harigua-Souiai, E. et al. Identification of binding sites and favorable ligand binding moieties by virtual screening and self-organizing map analysis. *BMC Bioinform* **16**, 93, https://doi.org/10.1186/s12859-015-0518-z (2015).

84. Schrodinger, LLC. *The PyMOL Molecular Graphics System, Version 1.8.2.1* (2015).

## Acknowledgements

This work was supported by the Tunisian National Center for Nuclear Sciences and Technology (CNSTN), the Institut Pasteur de Tunis (Tunis) and the Ministry of Higher Education and Scientific Research in Tunisia.

## Author Contributions

K.G. carried out bioinformatics analysis of data and drafted the manuscript. E.H.-S. performed the structural bioinformatics study and took part in writing the manuscript. C.B.H. helped to analyse the results and generated some of the figures. P.F. did the PEG tests. P.P. planned the PEG tests and analysed the results. S.G. collected data and researched literature. I.G. helped in the structural bioinformatics study. M.G. performed the protein extractions, the tandem mass spectrometry analysis and analysed the proteomics results. J.A. contributed to the proteomics results, performed the statistical analysis, and took part in writing the manuscript. P.N. planned the study, made genome annotation and took part in writing the manuscript. H.S. participated in the design and coordination of the study and drafted the manuscript. All authors read and approved the final manuscript.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-017-18839-0.

**Competing Interests:** The authors declare that they have 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_15.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) 2018