# The First Genomic and Proteomic Characterization of a Deep-Sea Sulfate Reducer: Insights into the Piezophilic Lifestyle of Desulfovibrio piezophilus

Nathalie Pradel, Boyang Ji, G Gimenez, E Talla, P Lenoble, Marc Garel, Christian Tamburini, Patrick Fourquet, Régine Lebrun, Philippe Bertin, et al.

► To cite this version:

Nathalie Pradel, Boyang Ji, G Gimenez, E Talla, P Lenoble, et al.. The First Genomic and Proteomic Characterization of a Deep-Sea Sulfate Reducer: Insights into the Piezophilic Lifestyle of Desulfovibrio piezophilus. PLoS ONE, 2013, 8 (1), pp.e55130. 10.1371/journal.pone.0055130. [hal-00787201](https://hal.science/hal-00787201v1)

HAL Id: hal-00787201

https://hal.science/hal-00787201v1

Submitted on 9 Oct 2018

**HAL** is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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 CC BY 4.0 - Attribution - International License
# The First Genomic and Proteomic Characterization of a Deep-Sea Sulfate Reducer: Insights into the Piezophilic Lifestyle of *Desulfovibrio piezophilus*

Nathalie Pradel¹*, Boyang Ji², Grégory Gimenez³, Emmanuel Talla², Patricia Lenoble⁴, Marc Garel¹, Christian Tamburini¹, Patrick Fourquet⁵, Régine Lebrun⁶, Philippe Bertin⁷, Yann Denis⁶, Matthieu Pophillat⁵, Valérie Barbe⁴, Bernard Ollivier¹, Alain Dolla²*

¹ Aix-Marseille Université, Université du Sud Toulon-Var, CNRS/INSU, IRD, MIO, UM110, Marseille, France, ² Aix-Marseille Université, CNRS, LCB, UMR 7283, Marseille, France, ³ UMR 6236, Faculté de Médecine, Marseille, France, ⁴ Laboratoire de Finition C.E.A., Institut de Génomique - Genoscope, Evry, France, ⁵ Plate-forme Protéomique, CIML, Marseille, France, ⁶ Plate-formes Protéomique et Transcriptomique FR3479, IBiSA Marseille-Protéomique. IMM - CNRS, Marseille, France, ⁷ UMR 7156, CNRS, Université Louis Pasteur, Strasbourg, France

## Abstract

*Desulfovibrio piezophilus* strain C1TLV30T is a piezophilic anaerobe that was isolated from wood falls in the Mediterranean deep-sea. *D. piezophilus* represents a unique model for studying the adaptation of sulfate-reducing bacteria to hydrostatic pressure. Here, we report the 3.6 Mbp genome sequence of this piezophilic bacterium. An analysis of the genome revealed the presence of seven genomic islands as well as gene clusters that are most likely linked to life at a high hydrostatic pressure. Comparative genomics and differential proteomics identified the transport of solutes and amino acids as well as amino acid metabolism as major cellular processes for the adaptation of this bacterium to hydrostatic pressure. In addition, the proteome profiles showed that the abundance of key enzymes that are involved in sulfate reduction was dependent on hydrostatic pressure. A comparative analysis of orthologs from the non-piezophilic marine bacterium *D. salexigens* and *D. piezophilus* identified aspartic acid, glutamic acid, lysine, asparagine, serine and tyrosine as the amino acids preferentially replaced by arginine, histidine, alanine and threonine in the piezophilic strain. This work reveals the adaptation strategies developed by a sulfate reducer to a deep-sea lifestyle.

**Citation:** Pradel N, Ji B, Gimenez G, Talla E, Lenoble P, et al. (2013) The First Genomic and Proteomic Characterization of a Deep-Sea Sulfate Reducer: Insights into the Piezophilic Lifestyle of *Desulfovibrio piezophilus*. PLoS ONE 8(1): e55130. doi:10.1371/journal.pone.0055130

**Editor:** Valerie de Crécy-Lagard, University of Florida, United States of America

**Received:** November 23, 2012; **Accepted:** December 18, 2012; **Published:** January 30, 2013

**Copyright:** © 2013 Pradel et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** The work was supported by Institut de Recherche pour le Développement and National Centre for Scientific Research grants. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: nathalie.pradel@ird.fr (NP); dolla@imm.cnrs.fr (AD)

## Introduction

Sulfur cycling is an important biogeochemical process on Earth, especially in oceans where the majority of sulfur is present as sulfate (SO₄²⁻). Several data clearly indicate that microbial mediated sulfate reduction is a major terminal electron accepting process in the degradation of organic matter in marine sediments [1,2,3,4]. It thus emphasizes the importance of sulfate-reducing microorganisms in both the sulfur and carbon cycles. Because 60% of the Earth's surface is composed of deep-sea environments, we can postulate that a significant part, and even the majority of sulfur cycling occurs in the deep water ecosystem and, therefore, deep-sea sulfate-reducing microorganisms should be considered as key actors in Earth's sulfur cycling. Metagenomic studies have revealed that sulfate-reducing bacteria (SRB) are largely represented among the microorganisms that inhabit deep-sea environments [5,6,7]. However, a few thermophilic and mesophilic SRB have been isolated from deep-sea environments [8,9,10,11,12,13]. Among them, the *Desulfovibrio* genus is well-represented with *D. profundus*, *D. hydrothermalis*, and *D. piezophilus* [11,12,13].

Among the physico-chemical parameters that are characteristic of deep-sea environments, one of the most important is the

increase of the hydrostatic pressure with depth (1 MPa for every 100 m depth). Therefore, hydrostatic pressure could be an important factor that influences sulfate reduction and the global sulfur cycle, especially in the deep ocean. Within these ecosystems, the hydrostatic pressure can reach up to 110 MPa (Mariana Trench) with an average hydrostatic pressure of approximately 38 MPa, corresponding to depth of 3800 m. Organisms that experience optimal growth at hydrostatic pressure higher than atmospheric pressure (0.1 MPa) are called "piezophiles" [14]. Until now, factors involved in piezophilic lifestyle have mainly been analyzed in two facultative anaerobic heterotrophic genera, *Photobacterium* and *Shewanella* [15,16,17,18,19]. However, these two genera do not accurately represent the complexity and diversity of the deep-sea ecosystems; anaerobes and SRB, as discussed above, are also important members of these biotopes that have to be considered [1,20,21,22].

We recently isolated the sulfate-reducing anaerobe *Desulfovibrio piezophilus* C1TLV30T from Mediterranean deep-sea wood falls at a depth of 1700 m [13]. Because appreciation of the occurrence of deep-sea wood falls is increasing, it makes them important organic-rich biotopes in these depth environments [23]. The degradation of the wood falls is expected to be successive, beginning with
cellulose and other wood components degradation by heterotrophic bacteria. Diverse microbial communities would then develop through fermentation of other compounds; sulfate-reducing bacteria would, in turn, use the produced H₂ and led to an enrichment of hydrogen sulfide at the surface of the wood. The mesophilic piezophile *D. piezophilus*, which grows optimally at 30°C and 10 MPa, provides thus an unique opportunity to study the adaptation strategies of SRB to hydrostatic pressure.

In this paper, we report the genomic sequence of *D. piezophilus* CITLV30T and a comparative analysis of its sequence with the genomes of non-piezophilic *Desulfovibrio* species. In addition, we performed a differential proteome analysis after growth at atmospheric (0.1 MPa) and high (10 MPa) hydrostatic pressure. This analysis highlighted the importance of several gene clusters and metabolic pathways, such as amino acid metabolism and transport, in the adaptation to hydrostatic pressure and provided new insight into the piezophilic lifestyle and the evolution of the genus *Desulfovibrio*.

## Results

### General Features of the *D. piezophilus* C1TLV30T Genome

The *D. piezophilus* C1TLV30T genome consists of a single circular molecule of 3,644,098 bp with a GC content of approximately 50% (Fig. 1). Although the GC bases are not evenly distributed throughout the chromosome, the average GC content of *D. piezophilus* is quite low compared to the values found in other *Desulfovibrio* species, with the exception of *D. salexigens*. Seven genomic islands that range from 8 to 30 kbp were identified in the C1TLV30T genome (Fig. 1; Table S1). This number is in the same range as for other *Desulfovibrio* species (Table S8). The genes in these islands encode mainly prophage-like elements, proteins of unknown function, and hypothetical proteins. Very few insertion sequences (IS)/transposase elements were found (Fig. 1B), contrary to *Photobacterium profundum* SS9 that has many [24]. It indicates that *D. piezophilus* genome has not been subjected to many genomic rearrangements during its evolution. The low number of IS in *D. piezophilus* could be one characteristic of this phylogenetic group, as it is in the same range as in other *Desulfovibrio* species (Fig. 1B).

The *D. piezophilus* genome contains 3,424 coding DNA sequences (CDSs). Approximately 77.16% of the CDSs are classified in at least one Cluster of Orthologous Group of Proteins (COGs). Compared to representative non-piezophilic *Desulfovibrio* genomes (*D. salexigens*, *D. vulgaris*, and *D. aespoensis*), *D. piezophilus* has a larger proportion of genes that encode proteins involved in energy production and conversion (C), cell cycle control, cell division, chromosome partitioning (D), amino acid transport and metabolism (E), and post-translational modification, protein turnover, chaperones (O) inorganic ion transport and metabolism (P), secondary metabolites biosynthesis, transport and catabolism (Q) (Fig. S1). Those functional classes may be relevant for the cellular adaptation of *D. piezophilus* to hydrostatic pressure. Conversely, proteins that are involved in coenzyme transport and metabolism (H) are less represented in *D. piezophilus* than in other *Desulfovibrio* spp.

### Energy Metabolism and Sulfate Respiration in *D. piezophilus*

The main mode of energy generation in *D. piezophilus* is sulfate reduction. As expected, the *D. piezophilus* genome contains genes that encode proteins that are required for sulfate reduction, including sulfate transport, activation to adenosine-5'-phosphosulfate (APS) and APS reduction (Table S2). Energy conservation

in SRB involves a number of periplasmic and cytoplasmic enzymes, as well as transmembrane complexes [25,26] that are detected in the *D. piezophilus* genome (Table S2). *Desulfovibrio* species are characterized by a very high level of c-type cytochromes [26]. Although the *D. piezophilus* genome does not contain an extensive number of c-type cytochrome-encoding genes, it encodes several multiheme c-type cytochromes and one monoheme cytochrome (Table S2). One of these cytochromes (DESPIv2_11010) is closely related to a cytochrome from the piezophile *Photobacterium profundum*. Notably, this cytochrome was expressed at significantly higher levels under high (10 MPa) than low (0.1 MPa) hydrostatic pressure conditions (Table S3).

### Lipid Biosynthesis Pathways in *D. piezophilus*

A change in membrane fluidity is a well-established response to hydrostatic pressure exposure. Unsaturated fatty acids have been shown to be important for the adaptation to hydrostatic pressure [19]. We have previously shown increases in the chain length of membrane fatty acids (FAs) and in the abundance of monounsaturated FAs when *D. piezophilus* was cultured at a high hydrostatic pressure [13]. All of the genes necessary for the biosynthesis of FAs, phospholipids and peptidoglycans are present in the *D. piezophilus* genome. Notably, the *fabF* gene (DESPIv2_11609), which encodes a β-ketoacyl-ACP synthase II involved in the production of monounsaturated FAs, is present in the *D. piezophilus* genome. Its activity has been suggested to be responsive to pressure in piezophilic bacteria [27]. More interestingly, the *ole1* gene (DESPIv2_10162), which encodes a stearoyl-CoA desaturase, is also present. The activity of *ole1* could increase the proportion of unsaturated FAs. It has been suggested that the *ole1* gene product was involved in membrane fluidity at high hydrostatic pressure [28]. In addition, it has been shown that the *P. profundum* SS9 *ole1* ortholog increased its expression when pressure rose from atmospheric to high hydrostatic pressure [29]. Because the other genes that were identified as likely to be involved in lipid biosynthesis are also found in non-piezophilic bacteria, changes in the membrane lipid composition of *D. piezophilus* could also be attributed to the regulation of the classical lipid biosynthesis pathways.

### Comparative Genomics Reveal Specific Gene Content and Functions in *D. piezophilus*

A comparative genomic approach was carried out by searching for homologs within the non-redundant (nr) database. Among the 3,424 CDSs, 241 ORFans specific to *D. piezophilus* were found (no hits versus nr), and additional 106 genes that were not related to genes found in other sulfate-reducing microorganisms were identified (Table S4). Among the *D. piezophilus*-specific genes, 33 are found in predicted genomic islands, particularly GI5 and GI6, which contain 17 and 7 genes, respectively. The 241 ORFans encode proteins of unknown function that could play important roles in the piezophilic capacity of *D. piezophilus*. Among the latter 106 genes, 12% encoded proteins belonging to either the P or E COG categories (metabolite or ion metabolism and transport), 11% to the K or L COG categories (transcription, recombination and repair mechanisms), and 44% encoded proteins that did not belong to any COG category. Among these genes, some could be related to the hydrostatic pressure adaptation in this strain. Regarding the organization of the 347 genes, 19 gene clusters (at least 3 genes) were identified (Table S4). Interestingly, phylogenetic analyses revealed that three clusters (DESPIv2_11206-11237, DESPIv2_11831-11838, and DESPIv2_12041-12113) included genes that are closely related to genes from microorganisms able to growth under high hydrostatic pressure, such as *Colwellia*
A

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

B

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

<table><thead><tr><th></th><th>D. piezophilus</th><th>D. salexigens</th><th>D. aesponeensis</th><th>D. vulgaris<br>Hildenborough</th></tr></thead><tbody><tr><td>Size (base pairs)</td><td>3,646,098</td><td>4,289,847</td><td>3,629,109</td><td>3,570,858</td></tr><tr><td>Predicted CDSs</td><td>3,424</td><td>3,807</td><td>3,304</td><td>3,764</td></tr><tr><td>G+C content (%)</td><td>49.95</td><td>47.09</td><td>62.56</td><td>63.14</td></tr><tr><td>Coding density</td><td>89.39</td><td>88.67</td><td>89.04</td><td>87.36</td></tr><tr><td>Mean CDS Length (bp)</td><td>971.77</td><td>999.17</td><td>978.07</td><td>849.22</td></tr><tr><td>Mean CDS GC%</td><td>49.98</td><td>47.1</td><td>62.16</td><td>61.56</td></tr><tr><td>Number of tRNA</td><td>58</td><td>89</td><td>52</td><td>68</td></tr><tr><td>rRNA operons</td><td>3</td><td>5</td><td>2</td><td>5</td></tr><tr><td>Transposable elements</td><td>24</td><td>2</td><td>31</td><td>21</td></tr><tr><td>Nosferatu repeated sequences (%)</td><td>2.00</td><td>2.13</td><td>4.00</td><td>5.34</td></tr><tr><td>Proteins with predicted function</td><td>2,306 (68.8)</td><td>2,678 (73.0)</td><td>2,441 (73.9)</td><td>2,106 (62.3)</td></tr><tr><td>Inorganic ion transport and metabolism (%)</td><td>266 (7.9)</td><td>322 (8.1)</td><td>273 (7.4)</td><td>207 (5.5)</td></tr><tr><td>Proteins without predicted function:</td><td></td><td></td><td></td><td></td></tr><tr><td>Conserved hypothetical protein</td><td>672 (20.0)</td><td>784 (20.6)</td><td>748 (22.6)</td><td>890 (26.3)</td></tr><tr><td>Hypothetical proteins</td><td>376 (11.2)</td><td>243 (6.4)</td><td>115 (3.5)</td><td>384 (11.4)</td></tr></tbody></table>
**Figure 1.** Overview of the *D. piezophilus* genome. **A.** Representation of the circular chromosome of the *D. piezophilus* strain C1TLV30T in the form of concentric circles. The outer scale designates the coordinates in 500,000 base pairs, and the red ticks indicate physical gaps between two contigs in the same scaffold. Circles display (from outside to inside) in: 1) the predicted coding sequences on the plus strand; 2) the predicted coding sequences on the minus strand (the coding sequences are colored according to their putative roles as determined by COG classification); 3) rRNA; 4) tRNA; 5) transposable elements; 6) predicted genomic islands (GI, Table S1); 7) proteins that vary in abundance depending on the pressure (red, proteins more abundant at 10 MPa; green, proteins less abundant at 10 MPa); 8) GC content using a 5 kbp sliding window with step size of 1 kbp (red indicates that the GC content is higher than the average GC content and green indicates that the GC content is lower than the average GC content); 9) GC skew (calculated with the same parameters used for the GC content), with red indicating levels greater than 0 and green indicating levels less than 0.
**B.** General features of the *D. piezophilus* C1TLV30T genome compared with those of three other *Desulfovibrio* sp. *D. aespoensis* and *D. salexigens* were chosen as they belong to the same phylogenetic group (Fig. S3) and *D. salexigens* is also a mesophilic marine strain as *D. piezophilus*. *D. vulgaris* Hildenborough was introduced as model organism for *Desulfovibrio* species [47]. Twenty-one amino acid residues are represented by the tRNA genes in each organism. The numbers in parentheses are the percentages of the total number of proteins.
doi:10.1371/journal.pone.0055130.g001

*psychrerythraea* [30] or *Photobacterium profundum* [17] (Fig. S2 and Table S4). We then used qRT-PCR to determine whether the expression levels of these genes were dependent on hydrostatic pressure. The gene DESPIv2_11220, a representative of the gene cluster DESPIv2_11206-11237, appeared to be up-regulated when the cells were cultured under high (10 MPa) conditions compared with low (0.1 MPa, atmospheric pressure) hydrostatic pressure conditions, while the expression of the representative genes of the two other gene clusters, DESPIv2_12106 and DESPIv2_11838, were not affected by the pressure (Table S3).

To gain further insight into the piezophilic lifestyle in bacteria, we compared the distribution of the COG classes of *D. piezophilus* and *P. profundum*, the most well-studied piezophilic bacterium, to that of *D. salexigens*, the only non-piezophilic marine *Desulfovibrio* for which a genome sequence is available, and *D. vulgaris* Hildenborough, the best known *Desulfovibrio* strain. This comparison revealed 53 COGs that were present in both piezophilic strains but absent in the non-piezophilic strains (Table S5). Among them, the most highly represented functional categories are amino acid transport and metabolism (E), secondary metabolite biosynthesis, transport and catabolism (Q), general functional prediction only (R) and unknown function (S) categories. Interestingly, we found proteins involved in the tripartite TRAP-type mannitol/chloro compound transport system (COG4663, COG4664, and COG4665). TRAP transporters have been reported to be involved in several mechanisms used by bacteria to adapt to different habitats [31]. This suggests that the transport of solutes plays an important role in the piezophilic lifestyle of *D. piezophilus*.

**D. piezophilus** has a Different Amino Acid Composition than other *Desulfovibrio* spp

A tree-based amino acid (AA) composition analysis of 21 selected strains, including piezophilic and non-piezophilic archaea, Gammaproteobacteria and *Desulfovibrio* organisms, revealed clusters of three main groups (Fig. 2). These microorganisms were chosen according to their taxonomic affiliation and their habitat (Fig. S3, Table S6). The first group (Group I) was composed of archaeal organisms and included known piezophiles (*Pyrococcus abyssi* and *Thermococcus barophilus*), while the second group (Group II) was composed of all of the analyzed Gammaproteobacteria (with three known piezophilic strains : *Photobacterium profundum*, *Shewanella piezotolerans*, *Shewanella violacea*). The third group (Group III) was composed of 7 *Desulfovibrio* strains. Interestingly, *D. piezophilus* and *D. salexigens* appeared in isolated branches at the depth level of *Desulfovibrio*. While known piezophilic bacteria and archaea shared similar AA compositions with their closely related phylogenetic strains, *D. piezophilus* exhibited an amino acid composition that was distinct from that of the other *Desulfovibrio* species. It is known that the GC content and the GC3 influence the AA composition of an organism [32,33]. Because *D. piezophilus*, *D. salexigens*, some archaea (e.g., *P. abyssi*) and Gammaproteobac-

teria (e.g., *Vibrio vulnificus*) had similar average GC contents and GC3 (Table S6), these two factors could not be responsible for the observed clustering pattern. Taken together, these results strongly suggest that *D. piezophilus* has a different AA composition than other *Desulfovibrio* or piezophilic strains, while other known piezophilic strains have AA compositions similar to those of their most closely related non-piezophilic relatives (Fig. 2). This suggests that the AA composition of *D. piezophilus* is a specific adaptation for a piezophilic lifestyle.

To study the link between AA composition and piezophilic capabilities, we studied the site-specific AA replacements in orthologous proteins from *D. piezophilus* and *D. salexigens*. Because previous investigations have shown the importance of the amino acid composition for adaptations to temperature and salinity [34,35,36], we chose to compare these two *Desulfovibrio* strains (1,911 orthologous proteins containing a total of 632,496 AAs), as they are both mesophilic and marine strains (Table S6). We thus expect to avoid any bias that could had been introduced with respect to different optimal salinity and temperature of growth. Table 1 shows the distribution of the total AA substitutions regarding a single AA between the two compared organisms. The shaded boxes indicate the statistically significant AA substitutions (*P*<0.0005). These data revealed that aspartic acid (D), glutamic acid (E), lysine (K), asparagine (N), serine (S) and tyrosine (Y) were preferentially replaced (2.33%, 6.82%, 13.10%, 8.10%, 5.30% and 1.57%, respectively) by other AAs in *D. piezophilus* orthologs. In particular, E, K and N were substituted at the statistically highest binary exchange rate. All of these AA were replaced mainly by alanine (A), histidine (H), arginine (R) and threonine (T) in *D. piezophilus*. These replacements were less frequent when *D. salexigens* was compared with *D. vulgaris* Hildenborough, indicating that this analysis excludes bias due to random drift mutations.

Other studies have also reported that arginine and alanine were the “most piezophilic AAs” and that they played a major role in adaptation to the piezophilic lifestyle in archaea [37,38]. Other specific AA substitutions (F versus L, G versus H, I versus L, and L versus W) were also found to be significant and thus could also be factors in the adaptation of *Desulfovibrio* sp. to hydrostatic pressure.?The observed substitutions did not induce any significant variation of the average isoelectric point value of the *D. piezophilus* CDSs (approximately 5), which was similar to the value calculated for other *Desulfovibrio* species, suggesting that this parameter is not involved in the adaptation of *D. piezophilus* to hydrostatic pressure.

**Effect of Hydrostatic Pressure on the Proteome of D. piezophilus**

Pressure-dependent variations in the proteome of *D. piezophilus* were evaluated by 1DE- and 2DE-gels for the membrane and soluble fractions, respectively (Fig. 3). The abundance of 24 protein spots or bands increased, while that of 16 protein spots or bands decreased, as a specific consequence of cell growth at low
<figure><img src="image_4.png" /><figcaption>**Figure 2. Clustering of piezophilic and non-piezophilic prokaryotes based on amino-acid composition.** The neighbor-joining tree inferred from the amino acid composition of the selected genomes, including piezophilic prokaryotes. Based on the amino acid frequency, the distance matrix was calculated with the Pearson correlation distance, and the inferred tree was constructed with the Phylip package using the neighbor-joining (NJ) method. The list of organisms with their taxonomic affiliations and their GC and GC3 content is shown in Table S6. The organisms with similar amino acid compositions are circled with a dotted line.
doi:10.1371/journal.pone.0055130.g002</figcaption></figure>

(0.1 MPa) compared to high (10 MPa) hydrostatic pressure. The corresponding identifications are indicated in Tables 2 and 3. As shown in Fig. 1A, the genes that encoded these proteins were scattered across the genome, demonstrating that there was no specific hotspot for high or low hydrostatic pressure genes within the *D. piezophilus* chromosome. Most of these proteins (11 proteins) were involved in AA transport and metabolism. Among them, the ABC transporter glutamine-binding protein GlnH (DE-SPIv2_10610) was more abundant under low hydrostatic pressure (Table 2). The gene encoding this protein was found to be up-regulated by qRT-PCR when cells were cultured under low hydrostatic pressure (Table S3). Interestingly, two copies of the gene that encodes GlnH (DESPIv2_10600 and DESPIv2_10610) have been found in the genome of *D. piezophilus*. A similar change

in gene expression was also observed for DESPIv2_10600. The Gln transport thus could be an important mechanism for adaptation of *D. piezophilus* to hydrostatic pressure. It is also noteworthy that the abundance of HemK (DESPIv2_10174), a protein involved in glutamine transformation, was increased after growth at high hydrostatic pressure (Table 3). This observation emphasizes the possible importance of glutamate/glutamine metabolism and/or transport in the adaptation of the bacterium to hydrostatic pressure.

Three proteins involved in AA biosynthesis were found to be more abundant under high than low hydrostatic pressure: (i) alanine dehydrogenase, involved in alanine biosynthesis, (ii) the ATP phosphoribosyltransferase HisG (DESPIv2_12562), involved in histidine biosynthesis, and (iii) the ornithine carbamoyltransfer-
<table><caption>**Table 1.** Number of AA substitutions for each AA in the comparison of 1,911 orthologous proteins between *D. piezophilus* and *D. salexigens*.</caption><thead><tr><th></th><th>AA</th><th>A</th><th>C</th><th>D</th><th>E</th><th>F</th><th>G</th><th>H</th><th>I</th><th>K</th><th>L</th><th>M</th><th>N</th><th>P</th><th>Q</th><th>R</th><th>S</th><th>T</th><th>V</th><th>W</th><th>Y</th></tr></thead><tbody><tr><td>D. salexigens</td><td>A</td><td>20257</td><td>374</td><td>609</td><td>1124</td><td>212</td><td>1565</td><td>182</td><td>567</td><td>734</td><td>696</td><td>314</td><td>317</td><td>690</td><td>501</td><td>493</td><td>2167</td><td>1283</td><td>1389</td><td>52</td><td>103</td></tr><tr><td></td><td>C</td><td>417</td><td>3674</td><td>45</td><td>40</td><td>84</td><td>95</td><td>30</td><td>89</td><td>42</td><td>139</td><td>68</td><td>59</td><td>33</td><td>30</td><td>82</td><td>260</td><td>175</td><td>212</td><td>14</td><td>44</td></tr><tr><td></td><td>D</td><td><strong>708</strong></td><td>33</td><td>13857</td><td>2843</td><td>58</td><td>755</td><td><strong>227</strong></td><td>89</td><td>624</td><td>119</td><td>72</td><td>1054</td><td><strong>345</strong></td><td>446</td><td><strong>296</strong></td><td>759</td><td><strong>554</strong></td><td>140</td><td>9</td><td>42</td></tr><tr><td></td><td>E</td><td><strong>1543</strong></td><td>51</td><td><strong>3217</strong></td><td>16624</td><td>66</td><td><strong>713</strong></td><td><strong>296</strong></td><td>213</td><td>1362</td><td><strong>331</strong></td><td>160</td><td>563</td><td><strong>529</strong></td><td><strong>1244</strong></td><td><strong>738</strong></td><td><strong>983</strong></td><td><strong>896</strong></td><td><strong>416</strong></td><td>21</td><td>58</td></tr><tr><td></td><td>F</td><td>217</td><td>64</td><td>36</td><td>67</td><td>11404</td><td>86</td><td>186</td><td>487</td><td>52</td><td><strong>1400</strong></td><td>289</td><td>63</td><td>65</td><td>76</td><td>76</td><td>169</td><td>170</td><td>372</td><td>151</td><td>1264</td></tr><tr><td></td><td>G</td><td>1509</td><td>121</td><td>741</td><td>612</td><td>62</td><td>24896</td><td><strong>144</strong></td><td>105</td><td>425</td><td>164</td><td>75</td><td>498</td><td>221</td><td>233</td><td>333</td><td>897</td><td>352</td><td>173</td><td>19</td><td>40</td></tr><tr><td></td><td>H</td><td>170</td><td>26</td><td>142</td><td>209</td><td>155</td><td>88</td><td>4823</td><td>69</td><td>189</td><td>145</td><td>55</td><td>240</td><td>79</td><td>288</td><td>282</td><td>191</td><td>154</td><td>87</td><td>28</td><td>307</td></tr><tr><td></td><td>I</td><td>527</td><td>97</td><td>72</td><td>169</td><td>494</td><td>89</td><td>67</td><td>14289</td><td>166</td><td><strong>3471</strong></td><td>774</td><td>79</td><td>113</td><td>126</td><td>169</td><td>172</td><td>523</td><td>4167</td><td>47</td><td>130</td></tr><tr><td></td><td>K</td><td><strong>989</strong></td><td>57</td><td><strong>725</strong></td><td><strong>1452</strong></td><td><strong>96</strong></td><td><strong>556</strong></td><td><strong>353</strong></td><td>229</td><td>12658</td><td><strong>426</strong></td><td>233</td><td>713</td><td><strong>392</strong></td><td><strong>1253</strong></td><td><strong>2689</strong></td><td><strong>877</strong></td><td><strong>1001</strong></td><td><strong>401</strong></td><td>28</td><td>84</td></tr><tr><td></td><td>L</td><td>641</td><td>164</td><td>120</td><td>280</td><td>1225</td><td>161</td><td>174</td><td>2854</td><td>286</td><td>27777</td><td>1943</td><td>103</td><td>219</td><td>329</td><td>356</td><td>258</td><td>505</td><td>2010</td><td><strong>142</strong></td><td>297</td></tr><tr><td></td><td>M</td><td>293</td><td>56</td><td>79</td><td>141</td><td>275</td><td>72</td><td>67</td><td>740</td><td>191</td><td>2146</td><td>6418</td><td>79</td><td>57</td><td>185</td><td>152</td><td>145</td><td>243</td><td>547</td><td>46</td><td>89</td></tr><tr><td></td><td>N</td><td><strong>491</strong></td><td>54</td><td>1090</td><td>629</td><td>63</td><td><strong>596</strong></td><td><strong>366</strong></td><td>108</td><td>685</td><td><strong>158</strong></td><td>80</td><td>7465</td><td><strong>161</strong></td><td><strong>419</strong></td><td><strong>429</strong></td><td><strong>881</strong></td><td><strong>581</strong></td><td>133</td><td>13</td><td>107</td></tr><tr><td></td><td>P</td><td>717</td><td>32</td><td>264</td><td>406</td><td>55</td><td>212</td><td>89</td><td>115</td><td>283</td><td>254</td><td>56</td><td>109</td><td>13594</td><td>195</td><td>165</td><td>499</td><td>298</td><td>201</td><td>23</td><td>32</td></tr><tr><td></td><td>Q</td><td>436</td><td>29</td><td>419</td><td>1057</td><td>59</td><td>221</td><td>308</td><td>108</td><td>805</td><td>358</td><td>167</td><td>274</td><td>177</td><td>6142</td><td>668</td><td>400</td><td>367</td><td>180</td><td>25</td><td>58</td></tr><tr><td></td><td>R</td><td>437</td><td>67</td><td>257</td><td>495</td><td>73</td><td>293</td><td>333</td><td>137</td><td>1969</td><td>304</td><td>132</td><td>261</td><td>164</td><td>617</td><td>13408</td><td>433</td><td>373</td><td>172</td><td>65</td><td>89</td></tr><tr><td></td><td>S</td><td><strong>2729</strong></td><td>234</td><td>813</td><td>807</td><td>114</td><td><strong>1033</strong></td><td><strong>260</strong></td><td>197</td><td>703</td><td>318</td><td>157</td><td>743</td><td>515</td><td>450</td><td><strong>534</strong></td><td>11048</td><td><strong>1915</strong></td><td>369</td><td>29</td><td>90</td></tr><tr><td></td><td>T</td><td>1138</td><td>153</td><td>379</td><td>571</td><td>128</td><td>290</td><td>144</td><td>461</td><td>530</td><td>469</td><td>255</td><td>397</td><td>294</td><td>318</td><td>376</td><td>1630</td><td>11069</td><td>810</td><td>29</td><td>63</td></tr><tr><td></td><td>V</td><td>1541</td><td>225</td><td>115</td><td>320</td><td>364</td><td>187</td><td>82</td><td>4079</td><td>271</td><td>2081</td><td>638</td><td>107</td><td>209</td><td>182</td><td>211</td><td>365</td><td>946</td><td>17069</td><td>39</td><td>122</td></tr><tr><td></td><td>W</td><td>35</td><td>12</td><td>8</td><td>18</td><td>122</td><td>22</td><td>34</td><td>39</td><td>18</td><td>95</td><td>31</td><td>9</td><td>12</td><td>17</td><td>42</td><td>24</td><td>26</td><td>34</td><td>3310</td><td>117</td></tr><tr><td></td><td>Y</td><td><strong>140</strong></td><td>50</td><td>65</td><td>81</td><td>1381</td><td>46</td><td><strong>371</strong></td><td>132</td><td>71</td><td>311</td><td>97</td><td>98</td><td>50</td><td>83</td><td>94</td><td>139</td><td>100</td><td>129</td><td><strong>169</strong></td><td>7621</td></tr></tbody></table>

doi:10.1371/journal.pone.0055130.t001
<figure><img src="image_5.png" /><figcaption>**Figure 3. Differential proteome analyses depending on hydrostatic pressure.** Proteome profiles of *D. piezophilus*. **A.** Overview of the proteome profile obtained by 2DE-gels of soluble proteins recovered from *D. piezophilus* grown at 10 MPa. The proteins (100 µg) were separated by 2D gel electrophoresis on a pH gradient of 4–7 and a 10–18% polyacrylamide denaturing gradient, and the gel was silver-stained. A total of 727 protein spots were detected. **B.** Overlapping DiGE images of the *D. piezophilus* proteome. The proteins were separated on a pH gradient of 3–10 and a 12% polyacrylamide denaturing gel (Cy3-green = 10 MPa; Cy5-red = 0.1 Mpa). A total of 940 protein spots were detected. **C.** Sodium dodecyl sulfate-polyacrylamide gel electrophoresis of the membrane protein fraction of *D. piezophilus* recovered at low (0.1 MPa, Lane 1) or high (10 MPa, Lane 2) hydrostatic pressure. Molecular weight standard, MW. The arrows indicate the positions of the differentially abundant proteins according to the hydrostatic pressure growth conditions.

doi:10.1371/journal.pone.0055130.g003</figcaption></figure>
<table><caption>**Table 2.** *D. piezophilus* proteins with increased abundance under low (0.1 MPa) relative to high (10 MPa) hydrostatic pressure.</caption><thead><tr><th>Spot number</th><th>Gene Accession</th><th>Annotation</th><th>COGs category</th><th>DA* (±0.2)</th></tr></thead><tbody><tr><td>1</td><td>DESPIv2_11412<sup>#</sup></td><td>Periplasmic [NiFe] hydrogenase small subunit, HynA-1</td><td>C</td><td>2.6</td></tr><tr><td>2</td><td>DESPIv2_11809</td><td>Adenylylsulfate reductase, α subunit, AprA</td><td>C</td><td>2.4</td></tr><tr><td>3</td><td>DESPIv2_12525<sup>#</sup></td><td>ATP synthase subunit b, AtpF</td><td>C</td><td>2.2</td></tr><tr><td>4</td><td>DESPIv2_20239</td><td>Nitroreductase</td><td>C</td><td>1.5</td></tr><tr><td>5</td><td>DESPIv2_20445</td><td>Cobyrinic acid ac-diamide synthase</td><td>C</td><td>2.8</td></tr><tr><td>6</td><td>DESPIv2_20246</td><td>Extracellular solute-binding protein family 1, PotD</td><td>E</td><td>2.7</td></tr><tr><td>7</td><td>DESPIv2_10610</td><td>ABC transporter glutamine-binding protein GlnH</td><td>E</td><td>2.9</td></tr><tr><td>8</td><td>DESPIv2_10611</td><td>ATP-binding component of ABC superfamily, AapP</td><td>E</td><td>1.6</td></tr><tr><td>9</td><td>DESPIv2_10063</td><td>Polar amino acid ABC transporter subunit</td><td>E</td><td>4.5</td></tr><tr><td>10</td><td>DESPIv2_10036</td><td>Extracellular solute-binding protein, family 5, DppA</td><td>E</td><td>1.9</td></tr><tr><td>11</td><td>DESPIv2_12641<sup>#</sup></td><td>Extracellular ligand-binding receptor</td><td>E</td><td>1.9</td></tr><tr><td>12</td><td>DESPIv2_10448</td><td>Triosephosphate isomerase, TpiA</td><td>G</td><td>2.4</td></tr><tr><td>13</td><td>DESPIv2_11436</td><td>Uncharacterized aldolase aq_1554</td><td>G</td><td>4.3</td></tr><tr><td>14</td><td>DESPIv2_12350</td><td>Ribonuclease, Rne/Rng family</td><td>J</td><td>4.4</td></tr><tr><td>15</td><td>DESPIv2_10222<sup>#</sup></td><td>Peptidoglycan-associated lipoprotein Pal</td><td>M</td><td>2.0</td></tr><tr><td>16</td><td>DESPIv2_12413<sup>#</sup></td><td>Outer membrane chaperone Skp (OmpH)</td><td>M</td><td>6.5</td></tr><tr><td>17</td><td>DESPIv2_12605<sup>#</sup></td><td>Putative AsmA family protein</td><td>M</td><td>3.4</td></tr><tr><td>18</td><td>DESPIv2_11075<sup>#</sup></td><td>Serine protease do-like DegP</td><td>O</td><td>9.6</td></tr><tr><td>19</td><td>DESPIv2_10232</td><td>Basic membrane lipoprotein</td><td>R</td><td>1.8</td></tr><tr><td>20</td><td>DESPIv2_12411</td><td>Tetratricopeptide domain protein</td><td>R</td><td>—</td></tr><tr><td>21</td><td>DESPIv2_10558</td><td>Conserved protein of unknown function</td><td>S</td><td>3.6</td></tr><tr><td>22</td><td>DESPIv2_11823<sup>#</sup></td><td>Methyl-accepting chemotaxis sensory transducer</td><td>T</td><td>2.2</td></tr><tr><td>23</td><td>DESPIv2_10098<sup>#</sup></td><td>Outer membrane efflux protein</td><td>U</td><td>2.0</td></tr><tr><td>24</td><td>DESPIv2_10982</td><td>Conserved protein of unknown function</td><td></td><td>3.7</td></tr></tbody></table>

<table><caption>**Table 3.** *D. piezophilus* proteins with decreased abundance under low (0.1 MPa) relative to high (10 MPa) hydrostatic pressure.</caption><thead><tr><th>Spot number</th><th>Gene Accession</th><th>Annotation</th><th>COGs category</th><th>DA* (±0.2)</th></tr></thead><tbody><tr><td>25</td><td>DESPIv2_11808</td><td>Quinone-interacting membrane bound oxidoreductase, Flavin protein QmoA</td><td>C</td><td>3.9</td></tr><tr><td>26</td><td>DESPIv2_10343</td><td>Amino acid-binding ACT domain protein</td><td>E</td><td>2.4</td></tr><tr><td>27</td><td>DESPIv2_10492<sup>#</sup></td><td>Extracellular solute-binding protein family 3</td><td>E</td><td>2.8</td></tr><tr><td>28</td><td>DESPIv2_11485</td><td>Extracellular solute-binding protein family 3</td><td>E</td><td>6.4</td></tr><tr><td>29</td><td>DESPIv2_12435</td><td>Alanine dehydrogenase</td><td>E</td><td>2.6</td></tr><tr><td>30</td><td>DESPIv2_12562</td><td>ATP phosphoribosyltransferase HisG</td><td>E</td><td>5.6</td></tr><tr><td>31</td><td>DESPIv2_10845</td><td>Ornithine carbamoyltransferase ArgF</td><td>F</td><td>1.8</td></tr><tr><td>32</td><td>DESPIv2_11386</td><td>TRAP dicarboxylate transporter DctP subunit</td><td>G</td><td>2.3</td></tr><tr><td>33</td><td>DESPIv2_10174</td><td>Modification methylase, HemK family</td><td>J</td><td>2.8</td></tr><tr><td>34</td><td>DESPIv2_11695</td><td>50 S ribosomal protein L19</td><td>J</td><td>3.8</td></tr><tr><td>35</td><td>DESPIv2_12480</td><td>Heat shock protein Hsp20</td><td>O</td><td>1.7</td></tr><tr><td>36</td><td>DESPIv2_20109</td><td>Conserved exported protein of unknown function</td><td>P</td><td>3.7</td></tr><tr><td>37</td><td>DESPIv2_11432</td><td>Putative Methyltransferase type 12</td><td>R</td><td>2.3</td></tr><tr><td>38</td><td>DESPIv2_12341</td><td>CBS domain-containing membrane protein</td><td>R</td><td>2.5</td></tr><tr><td>39</td><td>DESPIv2_10688</td><td>Sensor protein</td><td>T</td><td>2.6</td></tr><tr><td>40</td><td>DESPIv2_10745</td><td>Conserved protein of unknown function</td><td></td><td>2.4</td></tr></tbody></table>

<footnote>*Differential abundance.</footnote>

<footnote>#Identified in the membrane fraction by 1DE-gels.</footnote>

<footnote>doi:10.1371/journal.pone.0055130.t002</footnote>

<footnote>*Differential abundance.</footnote>

<footnote>#Identified in the membrane fraction by 1DE-gels.</footnote>

<footnote>doi:10.1371/journal.pone.0055130.t003</footnote>
ase ArgF (DESPIv2_10845), involved in arginine biosynthesis. Therefore, as suggested by the genome analysis discussed above, AA transport and metabolism appear to play major roles in the adaptation of *D. piezophilus* to hydrostatic pressure.

The levels of three proteins that are involved in the sulfate reduction pathways and energy metabolism were found to vary in a pressure-dependent manner: the flavin protein QmoA (DESPIv2_11808), the [NiFe] hydrogenase small subunit HynA-1 (DESPIv2_11412), and the adenylylsulfate reductase alpha subunit AprA (DESPIv2_11809). Notably, none of the corresponding genes exhibited significant changes in expression (Table S3), suggesting that either translational regulation or protein stability changes in response to hydrostatic pressure.

The outer membrane chaperone OmpH (DESPIv2_12413) showed a large abundance increase (differential abundance (DA) = 6.53) under low hydrostatic pressure. Similarly, the putative serine protease do-like DegP (DESPIv2_11075) was also found to be more abundant under this condition (DA = 9.58). This protein was previously reported to be involved in the adaptation to atypical conditions [39,40]. Thus, it could act as a chaperone in the adaptation to low hydrostatic pressure. Increases in chaperone proteins in response to low hydrostatic pressure has been previously shown in the piezophilic bacterium *P. profundum* S9 [29]. The genes that encode OmpH and DegP did not exhibit significant pressure-dependent variations in expression levels (Table S3). The variations in the levels of these proteins could thus be linked to either the increased stability of the protein under low hydrostatic pressure or a translational regulatory mechanism.

DNA replication and translation processes have also been reported to play a role in the adaptation of microorganisms to hydrostatic pressure [18,41,42]. In the case of *D. piezophilus*, our differential proteomic analysis revealed few proteins involved in these processes, with the exception of a ribonuclease (DESPIv2_12350) involved in tRNA maturation and a 50S ribosomal protein L19 (DESPIv2_11695) that exhibited a significant variation in abundance that was dependent on the pressure. These variations could be linked to a translational regulation process that depends on hydrostatic pressure.

Five "unknown function proteins" could be good candidates for a specific adaptation to hydrostatic pressure. Among them, the gene encoding DESPIv2_20109 was up-regulated by a factor of 9.2 in high compared to low hydrostatic pressure conditions. Interestingly, DESPIv2_10558 was within GI3 and was located upstream of a gene-specific cluster (DESPIv2_10562 to 10565).

## Discussion

We present here the first sequenced genome of a piezophilic sulfate-reducing bacterium, *Desulfovibrio piezophilus* C1TLV30T. This genome sequence provides an important opportunity to better understand the influence of deep-sea marine microorganisms on the global sulfur cycle and to gain insights into the adaptation strategies utilized by these microorganisms to survive at different hydrostatic pressures.

Although the genomic characteristics of *D. piezophilus* are globally similar to those of other *Desulfovibrio* species, several differences may define or influence the piezophilic lifestyle of sulfate-reducing bacteria. Among the unique characteristics of the *D. piezophilus* genome are seven genomic islands that mainly include genes encoding proteins of unknown function. These genomic islands could be linked to the ability of the organism to adapt to its environment, and they may contain genes that are involved in piezophily. Additionally, the cluster regions DESPIv2_11206-11237, DESPIv2_11831-11838, and DESPIv2_12041-12113 include specific genes to *D. piezophilus* and genes closely related to genes found in the piezophilic *Photobacterium profundum* [29] and *Colwellia psychrerythraea*, this latter being able to grow under deep-sea pressure [30]. These genes are also expected to play important roles in the piezophilic lifestyle of *D. piezophilus*. Further studies, including genetic studies and analyses of protein structure/function relationships, will be required to determine the roles of these proteins in hydrostatic pressure adaptation.

In archaea, the piezophilic lifestyle has been directly associated with codon usage and tetranucleotide frequency [43]. However, this is not the case in *D. piezophilus*, in which these genomic patterns are similar to those of closely related non-piezophilic *Desulfovibrio* strains. However, we showed that the AA composition of *D. piezophilus* proteins is distinct from that of other *Desulfovibrio* species, suggesting that AAs are key factors in the adaptation to hydrostatic pressure. In addition, a comparative analysis with *D. salexigens* showed that AAs are replaced by A, R, H and T with significant frequency in the *D. piezophilus* orthologs. "Piezophilic AAs" like arginine have also been identified in the archaean *Pyrococcus* species [38,43]. These AA replacements could be important for maintaining protein stability and/or activity at high hydrostatic pressure.

Our work identified specific genes/proteins that may be important for the adaptation of sulfate reducers to hydrostatic pressure. Enzymes involved in sulfate respiration (QmoA, HynA-1, and AprA), key components of the anaerobic energy metabolism of these microorganisms, are differentially abundant depending on the hydrostatic pressure. This suggests that the hydrostatic pressure has an impact on the dissimilatory sulfate reduction pathway in *D. piezophilus*. *S. violacea* adapts its respiratory components in response to alterations in pressure to conduct aerobic respiration [16]. One can conclude that the modulation of respiration pathways is a general mechanism by which microorganisms adapt to variations in hydrostatic pressure. The transport of AAs and various solutes and AA metabolism were also identified as major biological processes affected by hydrostatic pressure. Previous studies on piezophilic *Photobacterium* spp. and *Shewanella* spp. have also suggested that the adaptation to hydrostatic pressure involved changes in the transport of solutes [15,19,42]. These data suggest that the modification of proteins that are involved in transport represents a major general evolutionary adaptation strategy utilized by microorganisms in the deep-sea environment. Notably, the enzymes involved in the biosynthesis of arginine, histidine, and alanine were found to be more abundant under high than low hydrostatic pressure in *D. piezophilus*. These AAs may function to maintain the pH of the cytosol, as previously suggested by Abe [15]. In *P. profundum*, it has been suggested that "piezolites" could compensate for high hydrostatic pressure by maintaining osmotic pressure in the cell [44]. As AA transport and metabolism appear to be important cellular processes for the piezophilic lifestyle, we can propose that amino acids could also act as "piezolites" in *D. piezophilus*. Additional studies are required to gain further insight into this hypothesis.

## Materials and Methods

### Sequencing, Gene Identification and Genome Analysis

The complete genome sequencing of *Desulfovibrio piezophilus* C1TLV30T was performed by Genoscope (CEA, Evry, France), and the final assembly of the C1TLV30T genome contained a single replicon. The annotated sequence data were deposited in the EMBL database (accession number FO203427). The computational prediction of coding sequences (CDS) and the identifica-
tion of other genome features were automatically performed by the MicroScope platform (https://www.genoscope.cns.fr/agc/microscope/home/index.php). Functional assignments were created using the MicroScope automated annotation pipeline [45] with extensive manual curation for biological processes and metabolic pathways supported by Interpro, Figfam, TMHMM and BLAST analyses. Information on additional comparative genomic analyses can be found in Text S1.

## Identification of Homologs and AA Substitution Analysis

Homologs were identified by BLAST alignment against the nr database (release of January 2012) with an E-value cutoff of 10⁻³. The AA composition of orthologous proteins was compared between *D. piezophilus* C1TLV30T and *D. salexigens* DSM 2638. Orthologous proteins were defined using a reciprocal best BLASTP hit strategy with *D. salexigens* DSM 2638 as a reference genome. Additional thresholds of 70% coverage and 30% identity were added for increased stringency. Based on the reciprocal best BLAST hit, each orthologous pair was aligned using MUSCLE [46]. A systematic analysis of ungapped positions using an in-house perl script reported the numbers of conserved and substituted AAs.?For each AA, significant deviation from the expected 50:50 ratio was determined by calculating a χ² with one degree of freedom.

## Cell Culture

*D. piezophilus* C1TLV30T was cultured at the optimal temperature growth of 30°C, under either 0.1 MPa (atmospheric pressure, low hydrostatic pressure) or 10 MPa (high hydrostatic pressure) in the medium described by Khelaifia et al. [13]. A first run of growth at high or low hydrostatic pressure, which constituted the adaptation phase, was used to inoculate the large-scale cultures (200 mL) that were used for further proteomic and transcriptomic experiments. These large-scale cultures were inoculated at a 1:30 ratio with the adapted cultures, and the cells were grown at high or low hydrostatic pressure until the OD₆₀₀ nm = 0.2. Cells from both conditions were then harvested by centrifugation at 6000 g for 20 min at 4°C and washed once with 20 mL of 0.1 M Tris-HCl 0.15 M NaCl buffer (pH 7.6). The pellets were subsequently used for protein and RNA extractions.

## Differential Proteomic Analysis

Proteins were extracted for analysis by either 1D gel electrophoresis for the membrane fraction or 2D gel electrophoresis for the soluble fraction as described in Text S1. After electrophoresis, the gels were silver-stained using the PlusOne Silver Staining Kit (GE Healthcare). Data were acquired from the stained gels using the Image Scanner II and Image Master 2D software from Amersham Biosciences. The intensity levels of the spots were determined by the relative spot volume of each protein compared to the normalized volume of proteins. The relative intensity of each protein spot was compared between the control and the experimental groups. An intensity ratio higher than 1.5 was considered significant, and the corresponding spot/band was selected for protein identification by MALDI-TOF mass spectrometry or LC-ESI-MS/MS. The soluble fractions were also analyzed by the DiGE methodology. Details are provided in Text S1.

## Transcript Analysis by qRT-PCR

Total RNA was isolated using the High Pure RNA Isolation Kit (Roche Diagnostics), and the cDNA was subsequently synthesized from four independent cultures. The cDNA was quantified with

the SsoFast Evagreen Supermix 2X Kit (BioRad). The real-time PCR runs were carried out on a CFX96 Real-Time System using a C1000 Thermal Cycler (BioRad). The primer sequences are described in Table S7. The Relative Expression Software Tool (REST) was used to calculate the relative expression of each gene for each condition, and the 16S rRNA gene was used as a reference for normalization. The quantification was performed in duplicate or triplicate for each cDNA preparation. Details are provided in Text S1.

## Supporting Information

**Figure S1** Frequency distribution of the COG categories between *D. piezophilus* and selected *Desulfovibrio* genomes. For each organism, COG category with a proportion away from the 2 SDs of the mean of other three organisms is considered as a significant difference and is marked by a star (red star, >2 SD; black star, <2 SD).

(TIF)

**Figure S2** Phylogenetic analysis of genes DE-SPIv2_11220 (A), DESPIv2_11837 (B). The 20 best BLAST hit proteins versus nr were retrieved from Genbank. The phylogenetic reconstruction was achieved on the trimmed alignment using PhyML with a WAG matrix. Branch support estimated using a bootstrap of 100.

(TIF)

**Figure S3** 16S RNA based Phylogenetic tree of *Desulfovibrio* strains.

(TIF)

**Table S1** The putative genomic islands in *D. piezophilus* genome.

(PDF)

**Table S2** List of genes encoding proteins involved in energy metabolism.

(PDF)

**Table S3** Gene expression levels by qRT-PCR.

(PDF)

**Table S4** The 347 genes with no homologs in sulfate reducers.

(PDF)

**Table S5** List of the 53 COGs present in *D. piezophilus* and *P. profundum* but absent in the non-piezophilic *D. vulgaris* Hildenborough and *D. salexigens* strains.

(PDF)

**Table S6** Average GC and GC3 content of genomes in this study and strains characteristics.

(PDF)

**Table S7** List of oligonucleotide primers used.

(PDF)

**Table S8** The genomic islands in *Desulfovibrio* species.

(PDF)

**Text S1** Supplemental information on Materials and Methods.

(DOC)

## Acknowledgments

We thank M. Belghazi for the MS and LC-MS/MS data analysis and collection using the Q-ToF Micromass mass spectrometer from the CAPM Nord, Marseille Protéomique, Institut Jean Roche.
## Author Contributions

Conceived and designed the experiments: NP AD. Performed the experiments: NP AD VB PL MG CT PF RL YD MP. Analyzed the

## References

1. Fauque G, Ollivier B (2004) in Microbial Diversity and Bioprospecting, ed Bull AT (ASM Press, Washington), pp. 169-176.

2. Jorgensen BB (1982) Ecology of the bacteria of the sulphur cycle with special reference to anoxic-oxic interface environments. Philos Trans R Soc Lond B Biol Sci 298: 543-561.

3. Vossmeyer A, Deusner C, Kato C, Inagaki F, Ferdelman TG (2012) Substrate-specific pressure-dependence of microbial sulfate reduction in deep-sea cold seep sediments of the Japan Trench. Front Microbiol 3: 253.

4. D'Hondt S, Jorgensen BB, Miller DJ, Batzke A, Blake R, et al. (2004) Distributions of microbial activities in deep subseafloor sediments. Science 306: 2216-2221.

5. Hu Y, Fu C, Yin Y, Cheng G, Lei F, et al. (2010) Construction and preliminary analysis of a deep-sea sediment metagenomic fosmid library from Qiongdongnan basin, south China sea. Mar Biotechnol 12: 719-727.

6. Martin-Cuadrado AB, López-García P, Alba JC, Moreira D, Monticelli L, et al. (2007) Metagenomics of the deep Mediterranean, a warm bathypelagic habitat. PLoS One 2: e914.

7. Stokke R, Roalkvam I, Lanzen A, Hafidason H, Steen IH. (2012) Integrated metagenomic and metaproteomic analyses of an ANME-1-dominated community in marine cold seep sediments. Environ Microbiol 14: 1333-1346.

8. Alain K, Postec A, Grinsard E, Lesongeur F, Prieur D, et al. (2010) *Thermodesulfatator atlanticus* sp. nov., a thermophilic, chemolithoautotrophic, sulfate-reducing bacterium isolated from a mid-atlantic ridge hydrothermal vent. Int J Syst Evol Microbiol 60: 33-38.

9. Audifrin C, Cayol JL, Joulian C, Casalot L, Thomas P, et al. (2003) *Desulfonauticus submarinus* gen. nov., sp. nov., a novel sulfate-reducing bacterium isolated from a deep-sea hydrothermal vent. Int J Syst Evol Microbiol 53: 1585-1590.

10. Jeanthon C, L'Haridon S, Cueff V, Banta A, Reysenbach AL, et al. (2002) *Thermodesulfobacterium hydrogeniphilum* sp. nov., a thermophilic, chemolithoautotrophic, sulfate-reducing bacterium isolated from a deep-sea hydrothermal vent at Guaymas Basin, and emendation of the genus *Thermodesulfobacterium*. Int J Syst Evol Microbiol 52: 765-772.

11. Alazard D, Dukan S, Urios A, Verfé F, Bouabida N, et al. (2003) *Desulfocibrio hydrothermalis* sp. nov., a novel sulfate-reducing bacterium isolated from hydrothermal vents. Int J Syst Evol Microbiol 53: 173-178.

12. Bale SJ, Goodman K, Rochelle PA, Marchesi JR, Fry JC, et al. (1997) *Desulfovibrio profundus* sp. nov., a novel barophilic sulphate-reducing bacterium from deep sediment layers in the Japan Sea. Int J Syst Bacteriol 47: 515-521.

13. Khelaifia S, Fardeau ML, Pradel N, Aussignargues C, Garel M, et al. (2011) *Desulfovibrio piezophilus*, sp. nov., a novel piezophilic sulfate-reducing bacterium isolated from wood falls in Mediterranean Sea. Int J Syst Evol Microbiol 61: 2706-2711.

14. Yayanos AA (1995) Microbiology to 10,500 meters in the deep sea. Annu Rev Microbiol 49: 777-805.

15. Abe F (2007) Exploration of the effects of high hydrostatic pressure on microbial growth, physiology and survival: perspectives from piezophysiology. Biosci Biotechnol Biochem 71: 2347-2357.

16. Chikuma S, Kasahara R, Kato C, Tamegai H (2007) Bacterial adaptation to high pressure: a respiratory system in the deep-sea bacterium *Shewanella violacea* DSS12. FEMS Microbiol Lett 267: 108-112.

17. El-Hajj ZW, Tryfona T, Alcock DJ, Hasan F, Lauro FM, et al. (2009) Importance of proteins controlling initiation of DNA replication in the growth of the high-pressure-loving bacterium *Photobacterium profundum* S89. J Bacteriol 191: 6383-6393.

18. Lauro FM, Tran K, Vezzi A, Vito N, Valle G, et al. (2008) Large-scale transposon mutagenesis of *Photobacterium profundum* S89 reveals new genetic loci important for growth at low temperature and high pressure. J Bacteriol 190: 1699-1709.

19. Simonato F, Campanaro S, Lauro FM, Vezzi A, D'Angelo M, et al. (2006) Piezophilic adaptation: a genomic point of view. J Biotechnol 126: 11-25.

20. Elsgaard L, Guezennec J, Benbouzid-Rollet N, Prieur D (1995) Mesophilic sulfate-reducing bacteria from three deep-sea hydrothermal vent sites. Oceanol Acta 18: 95-104.

21. Kallmeyer J, Boetius A (2004) Effects of temperature and pressure on sulfate reduction and anaerobic oxidation of methane in hydrothermal sediments of Guaymas basin. Appl Environ Microbiol 70: 1231-1233.

22. Kaneko R, Hayashi T, Tanahashi M, Naganuma T (2007) Phylogenetic diversity and distribution of dissimilatory sulfite reductase genes from deep-sea sediment cores. Marine Biotechnol 9: 429-436.

data: NP AD BJ ET GG VB MP BO. Contributed reagents/materials/ analysis tools: PB BO. Wrote the paper: NP AD ET GG.

23. Fagervold SK, Galand PE, Zbinden M, Gaill F, Lebaron P, et al. (2012) Sunken woods on the ocean floor provide diverse specialized habitats for microorganisms. FEMS Microbiol Ecol, in press.

24. Campanaro S, Vezzi A, Vito N, Lauro FM, D'Angelo M, et al. (2005) Laterally transferred elements and high pressure adaptation in *Photobacterium profundum* strains. BMC Genomics 6: 122.

25. Keller KL, Wall JD (2011) Genetics and molecular biology of the electron flow for sulfate respiration in *Desulfovibrio*. Front Microbiol 2: 135.

26. Pereira IA, Ramos AR, Grein F, Marques MC, da Silva SM, et al. (2011) A comparative genomic analysis of energy metabolism in sulfate reducing bacteria and archaea. Front Microbiol 2: 69.

27. Bartlett DH (2002) Pressure effects on *in vivo* microbial processes. Biochim Biophys Acta 1595: 367-381.

28. Fernandes PM, Domitrovic T, Kao CM, Kurtenbach E (2004) Genomic expression pattern in *Saccharomyces cerevisiae* cells in response to high hydrostatic pressure. FEBS Lett 556: 153-160.

29. Vezzi A, Campanaro S, D'Angelo M, Simonato F, Vito N, et al. (2005) Life at depth: *Photobacterium profundum* genome sequence and expression analysis. Science 307: 1459-1461.

30. Methé BA, Nelson KE, Deming JW, Momen B, Melamud E, et al. (2005) The psychrophilic lifestyle as revealed by the genome sequence of *Colwellia psychrethrytrae* 34H through genomic and proteomic analyses. Proc Natl Acad Sci USA 102: 10913-10918.

31. Bakermans C, Sloup RE, Zarka DG, Tiedje JM, Thomashow MF (2009) Development and use of genetic system to identify genes required for efficient low-temperature growth of *Psychrobacter arcticus* 273-4. Extremophiles 13: 21-30.

32. Basak S, Banerjee T, Gupta SK, Ghosh TC (2004) Investigation on the causes of codon and amino acid usages variation between thermophilic *Aquifex aeolicus* and mesophilic *Bacillus subtilis*. J Biomol Struct Dyn 22: 205-214.

33. Wilquet V, Van de Castelee M (1999) The role of the codon first letter in the relationship between genomic GC content and protein amino acid composition. Res Microbiol 150: 21-32.

34. Casanueva A, Tuffin M, Cary C, Cowan DA (2010) Molecular adaptations to psychrophily: the impact of 'omic' technologies. Trends Microbiol 18: 374-381.

35. Mavromatis K, Ivanova N, Anderson I, Lykidis A, Hooper SD, et al. (2009) Genome analysis of the anaerobic thermohalophilic bacterium *Halothermothrix oreni*. PLoS One 4: e4192.

36. Rhodes ME, Fitz-Gibbon ST, Oren A, House CH (2010) Amino acid signatures of salinity on an environmental scale with a focus on the Dead Sea. Environ Microbiol 12: 2613-2623.

37. DiGiulio M (2005) A comparison of proteins from *Pyrococcus furiosus* and *Pyrococcus abyssi*: barophily in the physicochemical properties of amino acids and in the genetic code. Gene 346: 1-6.

38. Jun X, Lupeng L, Minjuan X, Oger P, Fengping W, et al. (2011) Complete genome sequence of the obligate piezophilic hyperthermophilic Archaeon *Pyrococcus yayanosi* CH1. J Bacteriol 193: 4297-4298.

39. Jomaa A, Iwanczyk J, Tran J, Ortega J (2009) Characterization of the autoclavage process of the *Escherichia coli* HtrA protein: implications for its physiological role. J Bacteriol 191: 1924-1932.

40. Skorko-Glonek J, Sobiecka-Szkatula A, Narkiewicz J, Lipinska B (2008) The proteolytic activity of the HtrA (DegP) protein from *Escherichia coli* at low temperatures. Microbiology 154: 3649-3658.

41. El-Hajj ZW, Alcock D, Tryfona T, Lauro FM, Sawyer L, et al. (2010) Insights into piezophily from genetic studies on the deep-sea bacterium, *Photobacterium profundum* S89. Ann N Y Acad Sci 1189: 143-148.

42. Oger PM, Jebbar M (2010) The many ways of coping with pressure. Res Microbiol 161: 799-809.

43. DiGiulio M (2005) The ocean abysses witnessed the origin of the genetic code. Gene 346: 7-12.

44. Martin DD, Bartlett DH, Roberts MF (2002) Solute accumulation in the deep-sea bacterium *Photobacterium profundum*. Extremophiles 6: 507-514.

45. Vallenet D, Engelen S, Mornico D, Cruveiller S, Fleury L, et al. (2009) MicroScope: a plateform for microbial genome annotation and comparative genomics. Database (Oxford) 2009: bap021.

46. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32: 1792-1797.

47. Heidelberg J, Seshadri R, Haveman S, Hemme C, Paulsen I, et al. (2004) The genome sequence of the anaerobic sulfate-reducing bacterium *Desulfovibrio vulgaris* Hildenborough. Nat Biotechnol 22: 554-559.