Next Article in Journal
The Role of Fucoxanthin as a Potent Nrf2 Activator via Akt/GSK-3β/Fyn Axis against Amyloid-β Peptide-Induced Oxidative Damage
Next Article in Special Issue
Thioredoxin-2 Regulates SqrR-Mediated Polysulfide-Responsive Transcription via Reduction of a Polysulfide Link in SqrR
Previous Article in Journal
Noninvasive Methods to Detect Reactive Oxygen Species as a Proxy of Seed Quality
Previous Article in Special Issue
The Transcriptional Repressor PerR Senses Sulfane Sulfur by Cysteine Persulfidation at the Structural Zn2+ Site in Synechococcus sp. PCC7002
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analysis of Cyclooctasulfur Oxidation and Reduction by the Neutrophilic Chemolithoautotrophic Sulfurovum indicum from Deep-Sea Hydrothermal Ecosystems

1
Key Laboratory of Marine Genetic Resources, Third Institute of Oceanography, Ministry of Natural Resources, Xiamen 361005, China
2
State Key Laboratory Breeding Base of Marine Genetic Resources, Key Laboratory of Marine Genetic Resources of Fujian Province, Xiamen 361005, China
3
Sino-French Laboratory of Deep-Sea Microbiology (MicrobSea), Xiamen 361005, China
4
Department of Bioengineering and Biotechnology, Huaqiao University, Xiamen 361021, China
5
CNRS, Université Brest, Ifremer, Unité Biologie et Ecologie des Ecosystèmes Marins Profonds BEEP, UMR 6197, IRP 1211 MicrobSea, IUEM, Rue Dumont d’Urville, F-29280 Plouzané, France
6
Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519000, China
*
Authors to whom correspondence should be addressed.
Antioxidants 2023, 12(3), 627; https://doi.org/10.3390/antiox12030627
Submission received: 20 January 2023 / Revised: 28 February 2023 / Accepted: 1 March 2023 / Published: 3 March 2023
(This article belongs to the Special Issue Reactive Sulfur Species in Microorganisms)

Abstract

:
Chemolithoautotrophic Campylobacterota are widespread and predominant in worldwide hydrothermal vents, and they are key players in the turnover of zero-valence sulfur. However, at present, the mechanism of cyclooctasulfur activation and catabolism in Campylobacterota bacteria is not clearly understood. Here, we investigated these processes in a hydrothermal vent isolate named Sulfurovum indicum ST-419. A transcriptome analysis revealed that multiple genes related to biofilm formation were highly expressed during both sulfur oxidation and reduction. Additionally, biofilms containing cells and EPS coated on sulfur particles were observed by SEM, suggesting that biofilm formation may be involved in S0 activation in Sulfurovum species. Meanwhile, several genes encoding the outer membrane proteins of OprD family were also highly expressed, and among them, gene IMZ28_RS00565 exhibited significantly high expressions by 2.53- and 7.63-fold changes under both conditions, respectively, which may play a role in sulfur uptake. However, other mechanisms could be involved in sulfur activation and uptake, as experiments with dialysis bags showed that direct contact between cells and sulfur particles was not mandatory for sulfur reduction activity, whereas cell growth via sulfur oxidation did require direct contact. This indirect reaction could be ascribed to the role of H2S and/or other thiol-containing compounds, such as cysteine and GSH, which could be produced in the culture medium during sulfur reduction. In the periplasm, the sulfur-oxidation-multienzyme complexes soxABXY1Z1 and soxCDY2Z2 are likely responsible for thiosulfate oxidation and S0 oxidation, respectively. In addition, among the four psr gene clusters encoding polysulfide reductases, only psrA3B3C3 was significantly upregulated under the sulfur reduction condition, implying its essential role in sulfur reduction. These results expand our understanding of the interactions of Campylobacterota with the zero-valence sulfur and their adaptability to deep-sea hydrothermal environments.

1. Introduction

Sulfur occurs in many oxidation states and is of central importance in biogeochemical cycles [1,2]. As an essential component of biomass, sulfur is assimilated into organic compounds and also participates in energy-yielding processes as an electron acceptor or donor in chemolithautotrophic, photolithoautotrophic and heterotrophic microorganisms [3,4,5,6]. Elemental sulfur is present in a variety of environments, including marine sediments, microbial mats, cold/hot springs, oxygen-minimum zones, glacial shields, volcanic soils and hydrothermal vents, where it always accumulates in amounts visible to the naked eye [7,8,9,10,11,12]. In these cases, sulfur occurs mostly in the form of zero-valence sulfur (S0), with cyclooctasulfur (S8) as the most stable and common form [13].
As elemental sulfur is extremely poorly soluble in water [14], an initial activation reaction, by direct or indirect interaction between sulfur and cells, is required prior to its uptake [15]. In acidophilic sulfur-oxidizing bacteria, a direct contact between cells and sulfur is necessary [16,17]. For example, in members of the genera Acidithiobacillus and Acidiphilium, extracellular S0 is activated by reacting with reactive thiol groups (-SH) of special outer membrane proteins, forming -SnH (n ≥ 2) complexes and then transporting it to the periplasm, where it is oxidized [18,19]. Extracellular polymers (e.g., exopolysaccharides EPS), cell surface proteins, extracellular vesicles and flagella play a key role in the attachment of bacterial cells to the S0 surface during this process [20,21,22]. In some sulfur-reducing microorganisms such as Acidilobus sulfurireducens 18D70 and Desulfurella amilsii TR1, the physical contact of cells and S0 is not mandatory [23,24]. These indirect interactions could occur through the formation of polysulfides generated by the nucleophilic attack on sulfur by sulfide, or through extracellular electron transport via pili or redox cofactors present in biofilms [25].
When S0 is activated and transported into the cell, it is catalyzed by different enzymes for oxidation or reduction to supply energy for growth. The bacterial oxidation of S0 has been extensively studied in various sulfur-oxidizing bacteria [26,27,28]. In the periplasmic space, various enzymes, such as sulfur-oxidizing-multienzyme complex (Sox), tetrathionate hydrolase (TetH) and thiosulfate dehydrogenase (Tsd), participate in the further oxidation of thiol-bound sulfane sulfur atoms (R-S-SnH) [17,29]. The Sox pathway, which is governed by the conserved operon SoxXYZABCD, operates in photo- and chemo-lithotrophic Alphaproteobacteria and Campylobacterota [7,30]. This enzymatic complex includes the c-type cytochromes SoxAX, the sulfur-compound-binding module SoxYZ, the sulfate thiol esterase SoxB, and a sulfur dehydrogenase Sox(CD)2 [31]. In the cytoplasmic space, sulfur can be oxidized to sulfite by a series of enzymes, such as sulfur oxygenase reductase (Sor), thiosulfate sulfurtransferase (Tst), heterodisulfide reductase (Hdr) and sulfur dioxygenase (Sdo) [19,27]. For the reduction of sulfur, three main enzymes including polysulfide reductase (Psr), sulfur reductase (Sre) and sulfide dehydrogenase (Sudh) have been purified and characterized in a limited number of S0 reducers, such as Wolinella succinogenes, Pyrococcus furiosus, Acidianus ambivalens and members of the Desulfurellaceae family [24,32,33]. The membrane-bound polysulfide reductase consists of three subunits (PsrABC). PsrA is the catalytic subunit of the Psr protein for polysulfide reduction to sulfide. PsrB on the periplasmic side of the membrane contains several [Fe-S] clusters, which likely mediate the electron transfer from the membrane anchor (PsrC) to the catalytic subunit (PsrA) of the Psr [15]. Psr contains menaquinone, which could serve as an electron acceptor for the hydrogenase, and could also be an electron donor for polysulfide/S0 reduction [32]. Our recent work revealed that two sulfur-reduction pathways performed by periplasmic and cytoplasmic Psr coexisted in chemolithoautotrophic Sulfurimonas from deep-sea hydrothermal vents, coupled with different energy conservation pathways [34].
The phylum Campylobacterota (formerly Epsilonproteobacteria) includes representatives that are primary producers of hydrothermal ecosystems and accounts for 66–98% of the microorganisms associated with the substrates of some hydrothermal vents [35,36,37]. The genus Sulfurovum, which belongs to Campylobacterota, was first isolated and identified in 2004 from deep-sea hydrothermal vent sediments [38]. It is widely distributed in sulfur-rich marine environments and is very abundant in the fluids of deep-sea hydrothermal vents worldwide [39,40]. In particular, members of the genus Sulfurovum are frequently found as endo- and/or epi-symbionts of various deep-sea hydrothermal vent animals [41]. It also has been suggested that Sulfurovum and Sulfurimonas species especially occupy the niche of S0/S8 oxidation in marine sediments and in hydrothermal vents [6,42]. All these observations tend to indicate that Sulfurovum taxa play an important ecological role in deep-sea hydrothermal vents. Furthermore, the isolates of Sulfurovum identified to date generally use H2 or reduced-sulfur compounds as electron donors, and oxygen, nitrate or elemental sulfur as terminal electron acceptors [38,43].
In our previous study, a deep-sea bacterium designated Sulfurovum indicum ST-419 was isolated from a hydrothermal plume on the Wocan-1 hydrothermal site of the northwestern Indian Ocean [43]. Wocan-1 on the Carlsberg Ridge (60°68′ E, 6°56′ N) was discovered in March 2017 during the COMRA DY 38 oceanic scientific cruise [43]. Interestingly, this strain can grow vigorously with S0 as the sole electron donor (S0 oxidation), as well as the terminal electron acceptor (S0 reduction). Sulfur-related energy metabolism in Sulfurovum is particularly important for issues related to the ecology and biogeochemistry of deep-sea hydrothermal ecosystems. In this study, we aimed to elucidate the mechanisms involved in the activation and transport of cylcooctasulfur, in addition to the S0 oxidation and reduction pathways. Understanding the interaction of Sulfurovum with cylcooctasulfur will help us gain insight into the role of Campylobacterota in various ecological niches.

2. Materials and Methods

2.1. Bacterial Strains and Growth Conditions

Cultures of strain ST-419 were routinely grown in a modified MMJS medium at 37 ℃ as previously described [35]. For sulfur oxidation cultures, either elemental sulfur (5%, w/v) or thiosulfate (10 mM) was added as the sole electron donor, and nitrate (10 mM) as the sole electron acceptor. In both cases, the headspace of the anaerobic bottles was filled with 80% N2/20% CO2 (200 kPa). For sulfur reduction cultures, S0 (5%, w/v) was added as the sole electron acceptor with molecular hydrogen as the sole electron donor in anaerobic bottles filled with a gas phase of 80% H2/20% CO2 (200 kPa). As a control of sulfur reduction, nitrate reduction was set in parallel, with nitrate (10 mM) as the sole electron acceptor under a gas phase of 80% H2/20% CO2 (200 kPa). S0 granules were heat-treated prior to use by incubation at 95 °C for 2 h, then stored at room temperature [34]. Bacterial growth was monitored by cell counting over the incubation using a phase-contrast microscope (Eclipse 80i, Nikon, Tokyo, Japan). The anaerobic bottles under S0 oxidation and reduction were submerged in an ultrapure water solution and sonicated (WVR Ultrasonic cleaner, 80 W) for 1 min. All cultivation tests were performed in triplicate. Growth rates (µ, h−1) were calculated at multiple time points on replicate cultures and averaged using: µ = (lnN2 − lnN1)/(t2 − t1), where N2 and N1 represent the number of cells mL−1 at time (hours of incubation) t2 and t1, respectively. Generation times (tg; measured in hours) were calculated using: tg = (ln2)/µ [44]. The morphology and size of the cells were observed with a transmission electron microscopy (Model JEM-1230; JEOL, Tokyo, Japan).

2.2. Chemical Analyses

For direct analysis, samples were taken from cultures with N2-flushed syringes. Concentrations of thiosulfate, sulfate and nitrate were measured after filtration (0.2 µm, cellulose-acetate) using ion chromatography (ICS-2100, Thermo Scientific, Waltham, MA, USA) [45]. The concentrations of hydrogen and nitrogen in the headspace of the bottles were determined by gas chromatography (CP-2002, GL-Science, Tokyo, Japan) [46]. Dissolved sulfide concentrations were determined using the methylene blue assay as previously described [47]. All cultivation tests were performed in triplicate.

2.3. Transcriptome Analysis

Cells of strain ST-419 cultivated with different electron donors and acceptors under four types of catabolic conditions including S0 oxidation, thiosulfate oxidation, S0 reduction and nitrate reduction were harvested during the exponential growth phase, and centrifuged at 8000 rpm for 10 min at 4 °C. Resulting bacterial pellets obtained were rapidly frozen in liquid nitrogen and used for total RNA extraction. Transcriptome analysis was performed by Novogene (Tianjin, China). Detailed protocols of these procedures, including library preparation, clustering and sequencing, and the data analyses, are described in the Supplementary Materials. RNA-Seq reads were deposited in GenBank under accession numbers SRR22829317 to SRR22829328.

2.4. Quantitative Real-Time PCR Analyses

Fifteen genes were chosen for verification of RNA-Seq data by qRT-PCR. Total RNAs from the four culture conditions described above were extracted using TRIzol® reagent (Invitrogen, Waltham, MA, USA) and were then reverse-transcribed into cDNA using PremixScript™ RT reagent Kit (Takara, Otsu, Shiga, Japan). RNA degradation and contamination were monitored on 1% agarose gels. RNA purity and concentration were checked using the NanoPhotometer® spectrophotometer (IMPLEN, Calabasas, CA, USA). The expression levels of the 15 selected genes were determined using SYBR® Premix Ex Taq™ (RR420, Takara, Japan) with a Light cycler®480 (Roche, Switzerland) according to the manufacturer’s recommendations. After the cycling protocol, a melting curve was generated in order to detect single-gene-specific peaks and to check for the presence of primer–dimer peaks. The amplification efficiency was analyzed as follows: E = 10(−1/Slope) – 1, and amplification efficiencies were approximately ranged from 98.86 to 99.99%. The 16S rRNA gene was used as an internal reference, and the relative gene expression was calculated using the 2−ΔΔCt method [48]. Specific primers for the 15 genes and 16S rRNA were designed using Primer 6.0 as shown in Table S1. Three independent biological replicates were performed for each condition, and three technical replicates were performed for each reaction.

2.5. Experiment with Dialysis Membranes to Test Cell Contact with S0 Granules

Elemental sulfur was enclosed in dialysis membranes in batch cultures to examine the need for physical contact between cells and bulk-solid sulfur during the S0 oxidation and reduction processes. S0 was added to dialysis membranes (Spectrum Laboratories, USA) with pore sizes of 6 to 8 kDa and 12 to 14 kDa, and then secured with dialysis clips. Before use, all dialysis membranes were incubated at 80 °C in sterile deionized water for 24 h to remove preservatives, and this process was repeated three times, with the deionized water being replaced each time [49]. Cultures grown with elemental sulfur exposed fully to the medium (no dialysis membrane) were used as positive controls. Uninoculated media containing sulfur that was sequestered in dialysis membranes were used as negative controls. The sulfide and sulfate production, as well as the cell density, were monitored as described above. All of the treatments were performed in triplicate.

2.6. Sequence and Phylogenetic Analyses

The gene clusters encoding proteins involved in polysulfide reductases were identified by a local BLASTP protein similarity search using previously characterized gene sequences (e.g., polysulfide reductase from W. Succinogenes and Sulfurovum sp. NBC37-1) or by ontology using functional search terms (e.g., sulfur, polysulfide). For the BLASTP analysis, we used an amino acid similarity cutoff of >30%, alignment coverage > 80% and an e-value cutoff of 1E-5. Sequences were then aligned using ClustalW [50] and viewed in an on-site program (http://www.bio-soft.net/sms/index.html, accessed on 17 August 2022). A phylogenetic tree was reconstructed by the maximum-likelihood method using RAxML version 8.2.11 with the GTR + CAT model [51]. PSI-BLAST and DELTA-BLAST [52] were used for the domain analysis of outer membrane proteins. BOCTOPUS (boctopus.bioinfo.se) and PRED-TMMB (http://bioinformatics.biol.uoa.gr////PRED-TMBB, accessed on 12 October 2022) were used to search for beta-barrel topology. Signal peptides were predicted using SignalP 4.0 (http://www.cbs.dtu.dk/services/SignalP/, accessed on 26 October 2022).

2.7. Statistical Analysis

The significant differences among groups were subjected to one-way analysis of variance (one-way ANOVA) and multiple comparisons using the SPSS 19.0 program. A statistical significance was defined in our study by p < 0.05 (indicated by * in all figures) or p < 0.01 (indicated by ** in all figures).

3. Results

3.1. Growth Kinetics under Two Different Sulfur Oxidation Conditions

Strain ST-419 can grow chemolithoautotrophically with cyclooctasulfur or thiosulfate as the sole electron donor and nitrate as the terminal electron acceptor in an MMJS medium, as we reported previously [43]. The growth curves of strain ST-419 as well as the substrate and product concentrations over time are shown in Figure 1 (Figure 1A,B). The S0 granules initially floated on the medium due to the hydrophobic surface, and then gradually dipped into the medium and sank to the bottom of the bottles with the onset of cell growth and sulfur oxidation. As shown in Figure 1A,B, cultures grown with S0 had a significantly longer lag phase, which is necessary for cells to activate solid S0 prior to uptake and metabolizing. Over 24 h of incubation, elemental sulfur was gradually converted to sulfate, the final product of sulfur oxidation, reaching a maximum concentration of 8.97 mM, and accompanied by active cell growth and nitrate reduction from 8 h to 22 h until nitrate completely converted into nitrogen (Figure 1A). This reaction is congruent with the following equation: S0 + 1.2 NO3 + 0.4 H2O→SO42– + 0.6 N2 + 0.8 H+ [53]. Compared to S0 oxidation, the overall growth was better in the thiosulfate-oxidizing culture. The growth rate of thiosulfate-grown cultures was 1.55-times higher than that of cultures grown on S0 (0.11 h−1 vs. 0.06 h−1), corresponding to doubling times of ~6.93 h and ~10.75 h, respectively (Figure 1A,B). The short-rod morphology and size of the cells was similar under both experimental conditions.

3.2. Growth Kinetics by Sulfur and Nitrate Reduction Coupled with Hydrogen Oxidization

Strain ST-419 can also grow chemolithoautotrophically with hydrogen as the sole energy source and elemental sulfur or nitrate as the sole electron acceptor [43]. Similarly, S0, which initially floated on the medium, gradually sank to the bottom of the bottle as cell growth and sulfur reduction began. S0-reducing cultures consistently grew more slowly than nitrate-reducing cultures, but similar cell densities could be achieved after 24 h (>9 × 107 cells mL−1) (Figure 1E,F). Growth kinetics revealed that the doubling times in sulfur-reducing culture and in nitrate-reducing culture were approximately 7.70 h and 3.87 h, respectively (Figure 1E,F). The dissolved sulfide concentration reached up to 2.46 mM in cultures grown by sulfur reduction after a 24 h incubation (Figure 1C). Cell density, hydrogen consumption and sulfide production increased in parallel (Figure 1C–F). The cell morphology was also similar under both culture conditions as a short rod without obvious change in cell size.

3.3. General Features of Transcriptomes for S0 Oxidation and Reduction Conditions

To examine changes in gene expression specifically associated with the oxidation and reduction of S0, we investigated the transcriptomes of strain ST-419 using cells harvested at four different conditions in two pairs: S0 oxidation vs. thiosulfate oxidation (S0 oxidation), and S0 reduction vs. nitrate reduction (S0 reduction). A principal component analysis (PCA) showed that sample variability among the experimental treatments was higher than in the biological replicates, with replicates clustering well-separated along the main axis (54.2% of total variance) (Figure S1A). The expression patterns of S0 oxidation and reduction transcriptomes were evaluated using a heat map, which revealed that the expression profiles were separated among S0 oxidation and reduction against controls samples, respectively (Figure S1B). A total of 2137 and 2140 gene transcripts were detected in the transcriptomes of strain ST-419 during oxidation of S0 and thiosulfate, respectively, corresponding to 97.62–97.76% of all protein coding genes. Among those transcripts, 430 differentially transcribed genes (DEGs) (padj < 0.05) were identified, of which 212 were significantly increased and 218 significantly decreased (Figure 2A). Similarly, 2139 and 2250 transcripts were detected during reduction of S0 and nitrate, respectively (97.71–98.22% of protein coding genes), with 533 and 607 significantly increased and decreased, respectively, during reduction of S0 compared with nitrate reduction (Figure 2B). The list of the transcripts significantly expressed under both S0 oxidation and reduction is given in Table S2. Among the significant upregulation transcripts, 144 genes were shared between S0 oxidation and reduction, accounting for 6.6% of the total predicted coding sequences (CDS) of the genome (Figure 2C). Similarly, 155 significant downregulated transcripts were shared between S0 oxidation and reduction, corresponding to 7.1% of the total CDS (Figure 2C). The transcripts significantly expressed under both S0 oxidation and reduction in common are shown in Table S2.

3.4. Genes Potentially Involved in S0 Activation and Their Differential Expression

Among the 144 genes significantly overexpressed in both S0 oxidation and S0 reduction, 29 gene transcripts associated with biofilm formation including exopolysaccharide synthesis, bacterial secretion, signal transduction and TonB-dependent transfer system were found (Figure 3). Among them, most of these genes were related to exopolysaccharide synthesis, including polysaccharide synthetase and glycosyltransferase, which were significantly upregulated by 1.42–2.68- and 1.53–21.79-fold in both conditions, respectively (Figure 3A). For bacterial secretion, we identified several genes encoding efflux transporters (IMZ28_RS00120, IMZ28_RS06520 and IMZ28_RS06525) belonging to the HlyD family, and the outer membrane efflux protein, TolC (IMZ28_RS06515), which were significantly upregulated by 1.65–2.83-fold and 3.85–9.64-fold in S0 oxidation and reduction, respectively (Figure 3B). Genes associated with type II secretion systems including GspH (IMZ28_RS01610), GspF (IMZ28_RS03550) and GspG (IMZ28_RS10950) were also significantly upregulated in S0 oxidation and reduction conditions (Figure 3B). The HlyD and TolC in type I secretion systems, as well as Gsp in type II secretion systems, are responsible for exporting polysaccharides to the cell surface [54,55]. In addition, a gene (IMZ28_RS08888) related to type I secretion systems, containing a cadherin tandem repeat domain, was upregulated by 1.51- and 1.77-fold under S0 oxidation and reduction (this difference was statistically different), respectively (Figure 3B). Two type IV pilin-related genes including a pilus assembly protein MshL (IMZ28_RS83530) and a pilus assembly ATPase CpaE (IMZ28_RS00215) were also significantly upregulated in S0 oxidation and reduction conditions. Cadherin-like domains and type IV pili are frequently involved in surface adhesion [56,57]. Furthermore, some of these common genes were implicated in the signal transduction, such as the two-component system and TonB-dependent receptors in the presence of S0 either oxidation or reduction (Figure 3C,D).

3.5. Genes Potentially Involved in S0 Uptake and Their Differential Expression

A total of seven genes encoding the outer membrane OprD-like porins were detected in the transcriptomic data (Figure 4A). Porins of the OprD family are diverse and allow the facilitated uptake of a variety of specific substrates [58,59], so they may be involved in the uptake and transport of sulfur into cells as detailed below. Three genes encoding porin (IMZ28_RS00565, IMZ28_RS07310 and IMZ28_RS07390) were significantly upregulated by 1.77–2.53- and 2.66–7.62-fold in S0 oxidation and reduction compared with the control samples, respectively (Figure 4A). Among them, the gene IMZ28_RS00565 had the highest number of transcripts with FPKM values of 9110 and 12,412 for the two conditions, respectively, except two genes expressed at much lower levels (Table S2). Furthermore, the OprD protein encoded by gene IMZ28_RS00565 shared 19% sequence identity with its homolog in Sulfurimonas denitrificans DSM1251 [60]. Protein domain analysis showed that the N-terminal amino acids 42-403 of IMZ28_RS00565 possessed high similarity to the major outer membrane protein OprD porin of Campylobacterota (pfam05538). The modeling of the structural topology revealed several potential transmembrane domains that could form a beta-barrel structure and a signal peptide, pointing to a role as an outer membrane protein [61]. However, based on the low overall similarity, a functional classification cannot be inferred at present. Additionally, two oprD genes (IMZ28_RS01315 and IMZ28_RS01320) were also highly expressed in S0 oxidation and reduction conditions. Among them, gene IMZ28_RS01315 had a higher expression level with FPKM values of 8743 and 9068 under both conditions, and was upregulated by 1.90- and 1.57-fold under S0 oxidation and reduction, respectively (Figure 4A; Table S2). Similarly, the other gene (IMZ28_RS01320) had the highest expression under S0 oxidation and reduction with FPKM values of 40,107 and 46,447, respectively, and was upregulated by 1.15- and 2.02-fold under both conditions (Figure 4A; Table S2). Furthermore, two oprD genes (IMZ28_RS01315 and IMZ28_RS01320) encoding proteins shared the highest amino acid sequence identity of 34% and 45% with those of S. denitrificans DSM1251, respectively. A phylogenetic analysis showed that they clustered together with the homologous genes of S. denitrificans DSM1251 and Sulfurimonas sp. CVO (Figure 4B), which have recently been proposed to be involved in the uptake or transport of S0 [60,62].
Additionally, several genes coding for transporters were differentially expressed under both S0 oxidation and S0 reduction conditions (Table S3). Among them, four genes encoding substrate binding proteins (IMZ28_RS00125, IMZ28_RS02590, IMZ28_RS09850 and IMZ28_RS09855) were detected (Table S3), which belonged to ABC-type (ATP binding cassette) transporters burning ATP to fuel substrate transport [63]. Three tripartite ATP-independent periplasmic (TRAP)-type transporters (IMZ28_RS00555 IMZ28_RS00560 and IMZ28_RS06185) were also identified (Table S3). These transporters could obtain energy to actively channel substrates from the extracellular environment to the cytoplasm by combining them with the thermodynamically favorable transport of a solute such as Na+ [64]. These different types of transporters are involved in the transport of various compounds, and their substrate spectrum might cover sulfur compounds. In addition, three DEGs encoding ATPase (IMZ28_RS01810, IMZ28_RS07155 and IMZ28_RS10785) were also discovered (Table S3). They could be involved in the synthesis of ATP in the process of oxidation and reduction of S0 by oxidative phosphorylation, when electrons released by catabolism (electrons derived from H2 in sulfur reduction, and from S0/H2S in sulfide oxidation) enter the respiratory chain.

3.6. Evaluation of the Necessity for Direct Cell-S0 Contact Using Dialysis Bag Incubation

To investigate whether direct physical contact between cells and sulfur granules is essential for sulfur oxidation and respiration by strain ST-419, dialysis membranes with two kinds of pore sizes (6–8 and 12–14 kDa) were used to separate cells from elemental sulfur. The results showed that direct cell contact with solid S0 was not required for bacterial growth and sulfide production during sulfur reduction (Figure 5). Compared to the condition without the dialysis membrane, both sulfur reduction and bacterial growth decreased to varying degrees in the presence of the dialysis bags (Figure 5). In detail, with the membranes of pore sizes 6–8 kDa and 12–14 kDa, sulfide production decreased by 49% and 34% compared to the control without membrane, respectively (Figure 5A). Correspondingly, the final cell yield decreased by 41% and 22%, respectively (Figure 5B). In contrast, no cell growth was observed in S0 oxidation treatments when S0 was sequestered in dialysis membranes. Thus, cells of the bacterium growing via S0 reduction did not require direct contact to S0 granules, whereas growing via S0 oxidation did require direct contact.

3.7. Expression of the Genes Involved in S0 Oxidation

High transcripts of the periplasmic Sox multienzyme complex, encoded by two gene operon, soxABXY1Z1 and soxCDZ2Y2, were detected, implying that strain ST-419 uses the Sox system for the oxidation of reduced sulfur compounds. However, both sox operons showed obvious different transcription levels between thiosulfate and S0 oxidation conditions (Figure 6A). The gene cluster soxABXY1Z1 was transcribed under both growth conditions overall, but transcripts were more abundant with thiosulfate (Figure 6A). In its upstream in the genome, there are genes dsrE (IMZ28_RS08580), soxW (IMZ28_RS08585) and soxH1 (IMZ28_RS08590), which were significantly upregulated by 2.17–3.83-fold under thiosulfate oxidation (Figure 6A), suggesting their involvement in this process. DsrE and SoxW are, respectively, the periplasmic and cytoplasmic thioredoxins, which take part in the electron transport process [65]. SoxH is annotated as a putative metallo-hydrolase, which is responsible for releasing sulfate from SoxY [31]. Notably, the soxCDY2 genes in the second gene cluster were quite highly expressed with a FPKM value of >15,000, and significantly upregulated by 1.47–1.69-fold under S0 oxidation compared to thiosulfate oxidation, with the exception of the gene soxZ2, which showed no significant variation under both culture conditions (Figure 6A). In addition, a soxH2-like gene (IMZ28_RS10560) adjacent to the soxCDY2Z2 gene cluster was expressed in raised abundance under S0 oxidation (Figure 6A). This indicates the important role of soxCDY2Z2H2 in the oxidation of S0.
Among the significantly upregulated genes, three genes including hdrB (IMZ28_RS05750) and two genes encoding rhodanese-like proteins (IMZ28_RS02435 and IMZ28_RS02440) might also be involved in S0 oxidation. The gene hdrB (IMZ28_RS05750) was significantly upregulated by 1.52-fold under S0 oxidation (Figure 6A). Structural domain analysis showed that the hdrB-encoding protein has the cysteine functional motif, which is thought to be important for S0 oxidation in acidophilic sulfur oxidizers [66]. Two genes encoding rhodanese-like protein (IMZ28_RS02435 and IMZ28_RS02440) were also significantly upregulated by 2.77–2.82-fold in S0 oxidation (Figure 6A), which may be involved in the electron transfer during S0 oxidation. In addition, a gene (IMZ28_RS07900) encoding MBL-fold metallo-hydrolase showed a low expression under S0 oxidation (Figure 6A), although recently it was proposed to be involved in sulfur oxidation [28].

3.8. Expression of the Genes Involved in S0 Reduction

Based on genome analysis, strain ST-419 contained four types of putative polysulfide reductase-encoding genes (psrA1B1C1, psrA2B2, psrA3B3C3 and psrA4B4C4), which may perform a sulfur reduction function. These four kinds of polysulfide reductase-encoding genes showed differential transcriptional expressions under S0 reduction compared to nitrate reduction (Figure 6B). Among them, the gene cluster psrA3B3C3 was significantly upregulated by 4.03–7.12-fold with the high FPKM value of 1986–8312 under S0 reduction (Figure 6B). Three genes in the gene cluster psrA1B1C1 also exhibited a high FPKM value of 2504–17,874, corresponding to an increase in gene expressions by 1.21–2.10-fold during S0 reduction compared to nitrate reduction (Figure 6B), although there were no significant differences between the two conditions. In contrast, the transcriptional levels of genes in psrA2B2 and psrA4B4C4 were quite low, with gene psrA2 showing a significantly upregulated expression and gene psrA4 showing a downregulated expression (Figure 6B). These results indicated that the polysulfide reductase encoded by psrA3B3C3 may be more important for sulfur reduction.
In addition, three genes encoding for sulfide:quinone oxidoreductase (Sqr) belonging to type II, type IV and type VI, which are well known to be responsible for the oxidation of sulfide [67], showed different expression profiles during S0 reduction (Figure 6B). Among them, the expression of type IV sqr (IMZ28_RS09870) and type VI sqr (IMZ28_RS10500) was higher and significantly upregulated by 2.07–2.32-fold during S0 reduction (Figure 6B), implying their certain role in the process of sulfur reduction. The involvement of Sqr towards the direction of the reduction of sulfur compounds has been shown in Sulfurovum sp. NBC37-1 and in the thermophilic bacterium Thermovibrio ammonificans [68,69], which oxidizes the quinone pool and contributes to the reduction of elemental sulfur. In addition, the gene encoding NADH-dependent sulfur reductase (Nsr, IMZ28_RS07000), which has been proposed to be involved in sulfur reduction in T. ammonificans [69], showed a low expression level and was downregulated during S0 reduction (Figure 6B).

3.9. qRT-PCR Analyses of Sulfur Oxidation and Reduction Gene Transcripts

In order to verify the differential transcription of the genes putatively involved in sulfur metabolism, we performed qRT-PCR on 15 selected genes encoding: the beta lactamase precursor soxH1 (IMZ28_RS08590), the thiosulfohydrolase soxB (IMZ28_RS08600), the sulfur oxidation c-type cytochrome soxA (IMZ28_RS08605), the thiosulfate oxidation carrier complex protein soxZ1 (IMZ28_RS08610), the thiosulfate oxidation carrier protein soxY1 (IMZ28_RS08615), the sulfur oxidation c-type cytochrome soxX (IMZ28_RS08620), the beta-lactamase domain protein soxH2 (IMZ28_RS10560), the thiosulfate oxidation carrier protein soxZ2 (IMZ28_RS10570), the c-type cytochrome soxY2 (IMZ28_RS10575), the sulfite dehydrogenase soxD (IMZ28_RS10580), the sulfur oxidation molybdopterin C protein soxC (IMZ28_RS10585), the molybdopterin oxidoreductase psrA1 (IMZ28_RS08440), the molybdopterin oxidoreductase psrA2 (IMZ28_RS03620), the molybdopterin oxidoreductase psrA3 (IMZ28_RS06785) and the molybdopterin oxidoreductase psrA4 (IMZ28_RS07170) (Table S1). qRT-PCR results indicated that the transcripts for soxH1, soxB, soxA, soxZ1, soxY1 and soxX were higher in thiosulfate oxidation than S0 oxidation, and among them, the expressions of soxH1 and soxB were significantly upregulated by 3.31- and 1.47-fold under thiosulfate oxidation (Figure S3A). The expression levels of soxH2, soxZ2, soxY2, soxD and soxC were greater upon S0 oxidation than thiosulfate oxidation. Among them, genes of soxY2, soxD and soxC were upregulated by 1.49-, 1.58- and 1.67-fold, respectively, under S0 oxidation as compared to thiosulfate oxidation (Figure S3A). Furthermore, qRT-PCR quantification of the Psr large-subunit-encoding genes (psrA1, psrA2, psrA3 and psrA4) revealed that only psrA3 was significantly upregulated by 3.94-fold in S0 reduction versus nitrate reduction, whereas psrA1 and psrA2 showed a weak upregulation and psrA4 showed a downregulation (Figure S3B). These results indicate that the expression patterns from RNA-seq analysis are reliable.

3.10. Phylogenetic and Sequence Analyses of Diverse Polysulfide Reductases

The phylogenetic analysis based on deduced amino acid sequences encoded by four psrA-like genes showed that they formed three separate clades (group I, group II and group III) (Figure 7). Furthermore, they shared a common ancestor with the PsrA of W. succinogenes, the thiosulfate reductase PhsA of Salmonella enterica and the sulfur reductase SreA of A. ambivalens (Figure 7), which have been experimentally demonstrated to reduce sulfur or polysulfide [70,71,72]. The PsrA1 in group I is present in all Sulfurovum species and formed a monophyletic cluster, suggesting this group might be conserved within genus Sulfurovum (Figure 7). The genes in psrA1B1C1 shared the highest sequence similarity of 81–89% with the known sulfur respiration enzymes PsrABC of Sulfurovum sp. NBC37-1 [67]. The PsrA2 from genus Sulfurovum clustered into one group, and clustered with its homologs from the genus Sulfurimonas and the sulfur reductase SreA from Aquifex aeolicus (Figure 7). The gene cluster psrA2B2 shared 57–68% sequence similarity with the homologs from Sulfurimonas hydrogeniphila NW10, which was recently proposed to be involved in sulfur reduction [34]. As for group III, phylogenetic analysis indicated that PsrA3 and PsrA4 of strain ST-419 clustered with the homologous gene from Sulfurovum sp. HSL1-3, and furthermore formed one branch with other polysulfide reductase from the genera Sulfurimonas and Nautilia (Figure 7). The gene cluster of psrA3B3C3 and psrA4B4C4 shared 20–65% sequence similarity, and both had the highest similarity of 15–46% and 13–40%, respectively, with the PsrABC from W. succinogenes [70]. Furthermore, we detected the distribution of the psr gene cluster in genus Sulfurovum from 5fivepure isolated genomes and 122 metagenomes currently available on NCBI database. The result revealed that the four psr genotypes were present in Sulfurovum spp. from various environmental niches, including hydrothermal vents, terrestrial biofilm, hydrocarbon-rich habitats and marine sediments (Table S4). This provides evidence that the genotype of strain ST-419 is shared by other Sulfurovum spp., and suggests that this mechanism of sulfur reduction may be environmentally significant.
The primary amino acid sequences of the four PsrAs from S. indicum ST-419 were aligned with those of homologous enzymes that were biochemically characterized and shown to reduce S0 to H2S in vitro (Figure S4). Sequence comparison showed that several important cysteine residues are conserved at different positions, as indicated by one or two asterisks on the Figure S4. Among them, a conserved Cys176 residue was found to be necessary for the sulfur reductase activity in A. aeolicus [73]; this residue is assumed to be redox-active and involved in a Cys-S intermediate during the catalytic cycle. This conserved cysteine was also found in the PsrA of W. succinogenes (Cys-173) and the PhsA of S. enterica (Cys-178), which is responsible for sulfur reduction [34,73]. In the genome of strain ST-419, the Cys residue is conserved in the four PsrA proteins at positions 171, 247, 183 and 195, respectively (two asterisks, Figure S4). This conserved cysteine is also detected in the PsrA of Sulfurovum sp. NBC37-1 (Cys-257) (Figure S4). Overall, these observations strongly suggest that the Psr-like enzymes might be directly involved in S0 reduction in strain ST-419. Furthermore, none of the four PsrA proteins contained a typical twin arginine motif or a signal peptide, suggesting that they probably work in the cytoplasm.

4. Discussion

Elemental sulfur is abundant in hydrothermal vents, and its associated catabolism by the dominant chemolithoautotrophic Campylobacterota remains poorly understood. The genus Sulfurovum has been confirmed as one of the most predominant members of deep-sea hydrothermal prokaryotic residents that support the unique chemoautotrophic ecosystem [37], and is strongly involved in the oxidation of reduced inorganic sulfur compounds [26,38,74,75]. In the present study, we investigated the mechanisms of cyclooctasulfur activation and metabolisms when it is used as an electron donor or acceptor during sulfur oxidation and reduction in the deep-sea hydrothermal vent bacterium S. indicum. To our knowledge, this is the first report about the activation and metabolism mechanisms associated with S0 utilization in the genus Sulfurovum.
Considering the extremely low solubility and reactivity of S0 [76,77], microorganisms most likely need a specific activation or solubilization mechanism to make S0 available for their energy metabolism. Biofilm formation is an important way for cells to adsorb S0 in some bacteria and archaea, as previously observed [78,79,80,81]. During this process, flagella play an important role in initial biofilm development, which have been suggested to mediate attachment to sulfur particles in acidophilic and neutrophilic sulfur-oxidizing bacteria, including Sulfurimonas of the phylum Campylobacterota [18,60,82,83]. In this study, obvious biofilms were also observed on the surface of S0 particles by SEM (Figure S2). However, unlike the previously studied bacteria, strain ST-419 and other members of Sulfurovum entirely lacked the genes for surface-associated flagellar proteins and the bacterial chemotactic system (Figure S5). In contrast, genes related to biofilm formation including those involved in exopolysaccharide synthesis, bacteria secretion, signal transduction and the TonB-dependent transfer system were significantly upregulated in the presence of S0 (Figure 3). It is well-known that exopolysaccharides are one of the main components of biofilms and play a key role in the initial bacterial attachment [84,85]. The bacterial secretion systems of type I (such as HlyD and TolC family proteins) and type II were also markedly overexpressed, with the possible effect of facilitating surface adhesion, as suggested by another study [78]. Additionally, the genes involved in signal transduction and TonB-dependent receptors, which are essential for biofilm formation, were also clearly induced in the presence of S0 (Figure 3C,D). Therefore, biofilm formation by Sulfurovum may be a crucial step in activating S0 for both sulfur oxidation and reduction as supported by morphological and transcriptomic data.
The high expressions of frequently observed outer membrane proteins (OMP) in cells oxidizing and reducing solid S0 (Figure 4) suggest their involvement in utilizing this substrate. Two oprD-like porin genes (IMZ28_RS01315 and IMZ28_RS01320) were highly expressed under S0 oxidation and reduction (Figure 4), and they shared more than 30% sequence identity with the homologous genes in S. denitrificans DSM 1251 and Sulfurimonas sp. CVO, which were shown to be involved in the utilization of solid S0 by a transcriptomic and/or proteomic analyses [60,62]. Furthermore, in contrast to what was found in the genus Sulfurimonas, we found another highly expressed gene oprD (IMZ28_RS00565) in strain ST-419, which showed a 2.53- and 7.63-fold increase under both conditions, respectively (Figure 4; Table S2). All three proteins possessed the conserved domain of the OprD family, which is likely to play an important role in Campylobacterota as the homologs are frequently found in the genomes of strains in this phylum. Thus, we proposed that these porins may play a key role in S0 uptake in Campylobacterota. In addition, in acidophilic sulfur-oxidizing bacteria, some outer membrane proteins, such as Omp40 and Omp44, are also considered to play a role in sulfur attachment and transportation [17,86,87]. However, the true membrane proteins involved in S0 transport have not yet been experimentally confirmed in either acidophilic or neutrophilic sulfur-oxidizing bacteria [17,60,62]. Thus, further research is needed to determine whether these outer membrane proteins are directly involved in S0 degradation/activation, or whether they facilitate S0 transport into cells.
In addition to the direct contact of cells to insoluble S0, a non-contact or cooperative mechanism (coexistence of contact and non-contact activation) could be involved in S0 activation. In this study, we found that strain ST-419 growing via S0 reduction did not require direct contact to S0, whereas cells growing via S0 oxidation did require direct access (Figure 5), which is consistent with our previous reports in other species of Sulfurovum and Sulfurimonas [34,75]. Thus, a mechanism of activation of S0 without cell contact also likely exists. Furthermore, we found some genes associated with cysteine and glutathione synthesis which were differentially upregulated under S0 oxidation and reduction (Table S5). Under S0 oxidation conditions, sulfur activation requires a direct cell contact, as cellular access to bulk sulfur contributes to the efficiency of the overall process by keeping cells and their substrate in close proximity, and avoiding the oxidation of the -SH-group-containing compounds (Figure 8). Under S0 reduction conditions, the sulfide produced by the metabolism, and/or the release of compounds containing -SH groups such as cysteine, GSH [88,89], could be involved in the activation of sulfur through a nucleophilic attack (Figure 8). However, in natural environments, contact and non-contact activation mechanisms can always coexist, in a way we propose naming cooperative activation, which is to some extent similar to the metal oxidation mechanism observed in acidophilic sulfur-oxidizing bacteria [20,21,90]. In the cooperative activation, the bacterial cells would continuously release the compounds containing a -SH group to catalyze the transition of the S0 ring to the -S-S group. Continuous release would lead the bacteria to consume substantial amounts of the high-energy substrate. To minimize substrate consumption, the bacteria would get as close as possible to the elemental sulfur surface. Thus, we propose that the cells should be in contact with S0 before the release of the reducing reagents. However, it is still not known which reducing substances the bacterium releases, and this question is expected to be solved by transcriptome analysis of dialysis bags or by metabolomics.
Transcriptome analysis of cultures growing by sulfur oxidation showed that the soxABXY1Z1 gene cluster was more expressed under thiosulfate oxidation condition, while the soxCDY2Z2 gene cluster was more expressed under S0 oxidation conditions (Figure 6A). Consistently, these two sox gene clusters were also differentially regulated by different sulfur compounds in Allochromatium vinosum and S. denitrificans [6,60,91,92]. It is possible that S0 oxidation in strain ST-419 is performed solely by the gene cluster soxCDY2Z2. Indeed, Sulfurimonas sp. CVO from an oil field is able to oxidize S0 to thiosulfate and sulfate with only soxCDY2Z2 in the absence of soxABXY1Z1 [62]. Thus, we propose that the S0 oxidation pathway is performed by Sulfurovum species as follows (Figure 8). First, circular S8 is transformed into linear polysulfide (Sn-, HSn) in a currently unknown way, and then is transported into the cellular periplasm by a transporter, possibly via an OprD-like porin. In the periplasmic space, polysulfide is covalently bound to a cysteine residue of SoxY2, and generates a thiocysteine-S sulfate residue (SoxZ2Y2-S-S7-S-). The outer sulfone sulfur of the cystein-persulfide on SoxY2 is then oxidized by SoxCD to form a cysteine-S sulfate residue (SoxZ2Y2-S-S7-SO2-), and subsequently oxidized to form SoxZ2Y2-S-S7-SO32-. Finally, this complex is hydrolyzed to sulfate by SoxH2 or by another way, and regenerates the SoxZ2Y2 complex (SoxZ2Y2-SH) (Figure 8). At present, we are not able to resolve how the sulfonate group bound to SoxY2Z2 is hydrolyzed to form sulfate; possibly by the putative periplasmic metallo-hydrolase encoded by IMZ28_RS10560. This scenario seems plausible as this gene is located next to soxCDY2Z2 and its homologs can be retrieved in many other sulfur-oxidizing Campylobacterota with conserved synteny as a thiol hydrolase. Another possibility might be that SoxCDY2Z2 can catalyze the reaction by itself, as previously supposed in S. denitrificans DSM1251 [60], though there is no relevant experimental evidence.
During sulfur reduction, four polysulfide reductases of three groups including group I (psrA1B1C1), group II (psrA2B2) and group III (psrA3B3C3 and psrA4B4C4) were differentially expressed (Figure 6B and Figure 7). Protein domain analysis showed that all these proteins were located in the cytoplasm, implying that S. indicum performed a cytoplasmic sulfur reduction. Among them, the transcript of psrA3B3C3 was highly abundant and significantly upregulated, implying its essential role in sulfur reduction. This is significantly different from our recent report of a Sulfurimonas isolate, which used both periplasmic and cytoplasmic polysulfide reductases, encoded by genes psrA1B1CDE and psrA2B2, respectively, to perform cyclooctasulfur reduction [34]. Together with the results in Sulfurimonas spp., cytoplasmic sulfur reduction seems to be a crucial catabolic pathway in the phylum Campylobacterota. The overall mechanism proposed for sulfur respiration in strain ST-419 is summarized in Figure 8. The reduction of elemental sulfur could be conducted by an electron transport chain from molecular hydrogen, via the upregulated membrane-bound hydrogenase, with menaquinone (MK) as electron carriers in the membrane, to reduce polysulfide from producing H2S by PsrA3B3C3 (Figure 8). The end product of H2S then diffuses outside the cell and helps convert bulk sulfur to dissolved polysulfide. The reduction of sulfur and HS- diffusion out of the cell allows the formation of a proton gradient. However, at this point, it is difficult to explain the coexistence of four different cytoplasmic polysulfide reductases. We speculate that these polysulfide reductases may facilitate the host adaptation to variations in the dynamic environments of hydrothermal vents, which needs further investigations.

5. Conclusions

In this report, we investigated the processes of activation, uptake and subsequent oxidation and reduction of cyclooctasulfur in Sulfurovum indicum, a neutrophilic chemolithoautotrophic bacterium of the phylum Campylobacterota that is predominant in deep-sea hydrothermal ecosystems. We described the key genes and metabolic pathways involved in biofilm formation, sulfur uptake, periplasmic sulfur oxidation and cytoplasmic sulfur reduction, coupled with nitrate reduction and hydrogen oxidation, respectively. Our isolate and other Sulfurovum genomes entirely lack the genes for surface-associated flagellar proteins and bacterial chemotactic systems. We propose that a cooperative mechanism may exist for S0 activation, which would involve the reducing compounds such as Cys, GSH and H2S. Transcriptomic data indicated that the complex soxCDY2Z2H2 plays a role in S0 oxidization in the periplasm, while among the four polysulfide reductases, psrA3B3C3 plays a more important role in sulfur reduction in cytoplasm. These mechanisms may be applicable to other Campylobacterota. The results of this study provide a better understanding of how cells derive energy from elemental sulfur, which is abundant in the current and past marine ecosystem. In the future, further genetic and biochemical investigations will be needed to confirm the genes involved in S0 activation, and to validate the tentatively proposed models for sulfur oxidation and reduction.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/antiox12030627/s1, Figure S1: Principal component analysis (PCA) and heatmap analysis of RNA-seq transcriptome data demonstrated S0-utilized gene expression patterns. A PCA plot of all RNA-seq samples showing how samples subjected to different treatments clustered separately (A). A heatmap including the gene cluster analysis done with log10 (FPKM+1) values (B); Figure S2: Scanning electron micrographs of S0 particles collected from cultures of strain ST-419 grown by S0 oxidation (B) and by S0 reduction (D). Smooth sulfur granules from uninoculated media was the control (A,C); Figure S3: qRT-PCR analyses of the expression of genes encoding S0 oxidation (A) and S0 reduction (B) of strain ST-419; Figure S4: Multiple sequence alignment of the region containing the proposed molybdenum ligand from the catalytic subunit of enzymes similar to polysulfide/sulfur/thiosulfate reductases; Figure S5: The clusters of orthologous genes (COGs) in the genome of strain ST-419. Table S1: Primers for qRT-PCR used in this study; Table S2: he differentially up-regulated and down-regulated transcripts in S0 oxidation and reduction compared to the controls, as well as the common genes including up and down regulated in both conditions; Table S3: Differential expression of genes involved in transport and ATPase under S0 oxidation and reduction; Table S4: Distribution of genes encoding polysulfide reductase in the genomes of cultured isolates and metagenome-assembled genomes (MAGs) of the genus Sulfurovum; Table S5: Differential expression of genes involved in cysteine and glutathione synthesis under S0 oxidation and reduction.

Author Contributions

Conceptualization, S.W. and S.X.; investigation, S.W., L.J. and Z.S.; supervision, L.J., K.A. and Z.S.; formal analysis, S.W., L.C., S.X. and Z.S.; visualization, S.W.; writing—original draft, S.W.; writing—review and editing, S.W., L.J., K.A. and Z.S.; funding acquisition, S.W., L.J., K.A. and Z.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Financial support was provided by the China Ocean Mineral Resources R&D Association (COMRA) program (No. DY135-B2-01 to Z.S.), the National Natural Science Foundation of China (No. 42176134 to L.J.), China Postdoctoral Science Foundation (No. 2021TQ0397 to S.W.), National Key Research and Development Program of China (No. 2018YFC0310701 to Z.S.), the Scientific Research Foundation of Third Institute of Oceanography, MNR (No. 2019021 to Z.S.), the IRP 1211 MicrobSea (No. IRP1211 to K.A. and Z.S.) the National Natural Science Foundation of China (No. 91951201 to X.Z.) and the National Key R&D Program of China (No. 2021YFF0501304 to X.Z.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article and Supplementary Material.

Acknowledgments

We thank Li Gu (Third Institute of Oceanography) for the support with the scanning electron microscopy analysis.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Sievert, S.M.; Kiene, R.P.; Schulz-Vogt, H. The sulfur cycle. Oceanography 2007, 20, 117–123. [Google Scholar] [CrossRef] [Green Version]
  2. Wasmund, K.; Mußmann, M.; Loy, A. The life sulfuric: Microbial ecology of sulfur cycling in marine sediments. Environ. Microbiol. Rep. 2017, 9, 323–344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Friedrich, C.G.; Bardischewsky, F.; Rother, D.; Quentmeier, A.; Fischer, J. Prokaryotic sulfur oxidation. Curr. Opin. Microbiol. 2005, 8, 253–259. [Google Scholar] [CrossRef] [PubMed]
  4. Gregersen, L.H.; Bryant, D.A.; Frigaard, N.U. Mechanisms and evolution of oxidative sulfur metabolism in green sulfur bacteria. Front. Microbiol. 2011, 2, 116. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Sievert, S.M.; Vetriani, C. Chemoautotrophy at deep-sea vents: Past, present, and future. Oceanography 2012, 25, 218–233. [Google Scholar] [CrossRef]
  6. Meier, D.V.; Pjevac, P.; Bach, W.; Hourdez, S.; Girguis, P.R.; Vidoudez, C.; Amann, R.; Meyerdierks, A. Niche partitioning of diverse sulfur-oxidizing bacteria at hydrothermal vents. ISME J. 2017, 11, 1545–1558. [Google Scholar] [CrossRef] [Green Version]
  7. Taylor, C.D.; Wirsen, C.O. Microbiology and ecology of filamentous sulfur formation. Science 1997, 277, 1483–1485. [Google Scholar] [CrossRef]
  8. Engel, A.S.; Porter, M.L.; Stern, L.A.; Quinlan, S.; Bennett, P.C. Bacterial diversity and ecosystem function of filamentous microbial mats from aphotic (cave) sulfidic springs dominated by chemolithoautotrophic-Epsilonproteobacteria. FEMS Microbiol. Ecol. 2004, 51, 31–53. [Google Scholar] [CrossRef] [Green Version]
  9. Lavik, G.; Stührmann, T.; Brüchert, V.; Plas, A.V.; der Mohrholz, V.; Lam, P.; Mußmann, M.; Fuchs, B.M.; Amann, R.; Lass, U.; et al. Detoxification of sulphidic African shelf waters by blooming chemolithotrophs. Nature 2009, 457, 581. [Google Scholar] [CrossRef]
  10. Findlay, A.J.; Gartman, A.; MacDonald, D.J.; Hanson, T.E.; Shaw, T.J.; Luther, G.W. Distribution and size fractionation of elemental sulfur in aqueous environments: The Chesapeake Bay and mid-Atlantic ridge. Geochim. Cosmochim. Acta 2014, 142, 334–348. [Google Scholar] [CrossRef]
  11. Lau, G.E.; Cosmidis, J.; Grasby, S.E.; Trivedi, C.B.; Spear, J.R.; Templeton, A.S. Low-temperature formation and stabilization of rare allotropes of cyclooctasulfur (β-S8 and γ-S8) in the presence of organic carbon at a sulfur-rich glacial site in the Canadian High Arctic. Geochim. Cosmochim. Acta 2017, 200, 218–231. [Google Scholar] [CrossRef] [Green Version]
  12. Cai, R.; He, W.; Liu, R.; Zhang, J.; Zhang, X.; Sun, C. Deep-sea in situ insights into the formation of zero-valent sulfur driven by a bacterial thiosulfate oxidation pathway. mBio 2022, 13, e0014322. [Google Scholar] [CrossRef]
  13. Kafantaris, F.C.; Druschel, K.G. Kinetics of the nucleophilic dissolution of hydrophobic and hydrophilic elemental sulfur sols by sulfide. Geochim. Cosmochim. Acta 2019, 269, 554–565. [Google Scholar] [CrossRef]
  14. Steudel, R. Mechanism for the formation of elemental sulfur from aqueous sulfide in chemical and microbiological desulfurization processes. Ind. Eng. Chem. Res. 1996, 35, 1417–1423. [Google Scholar] [CrossRef]
  15. Zhang, L.; Qiu, Y.Y.; Zhou, Y.; Chen, G.H.; van Loosdrecht, M.C.M.; Jiang, F. Elemental sulfur as electron donor and/or acceptor: Mechanisms, applications and perspectives for biological water and wastewater treatment. Water. Res. 2021, 202, 117373. [Google Scholar] [CrossRef]
  16. Arredondo, R.; García, A.; Jerez, C. Partial removal of lipopolysaccharide from thiobacillus ferrooxidans affects its adhesion to solids. Appl. Environ. Microbiol. 1994, 6, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Rohwerder, T.; Sand, W. The sulfane sulfur of persulfides is the actual substrate of the sulfur-oxidizing enzymes from Acidithiobacillus and Acidiphilium spp. Microbiology 2003, 149, 1699–1710. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Mangold, S.; Valdés, J.; Holmes, D.; Dopson, M. Sulfur metabolism in the extreme acidophile Acidithiobacillus caldus. Front. Microbiol. 2011, 2, 17. [Google Scholar] [CrossRef] [Green Version]
  19. Chen, L.; Ren, Y.; Lin, J.; Liu, X.; Pang, X.; Lin, J. Acidithiobacillus caldus sulfur oxidation model based on transcriptome analysis between the wild type and sulfur oxygenase reductase defective mutant. PloS ONE 2012, 7, e39470. [Google Scholar] [CrossRef]
  20. Bouchal, P.; Zdrahal, Z.; Helanova, S.; Janiczek, O.; Hallberg, K.B.; Mandl, M. Proteomic and bioinformatic analysis of iron- and sulfur-oxidizing Acidithiobacillus ferrooxidans using immobilized pH gradients mass spectrometry. Proteomics 2010, 6, 4278–4285. [Google Scholar] [CrossRef]
  21. He, H.; Xia, J.; Huang, G.; Jiang, H.C.; Tao, X.X.; Zhao, Y.D.; He, W. Analysis of the elemental sulfur bio-oxidation by Acidithiobacillus ferrooxidans with sulfur K-edge XANES. World J. Microbiol. Biotechnol. 2010, 27, 1927–1931. [Google Scholar] [CrossRef]
  22. Zhan, Y.; Yang, M.; Zhang, S.; Zhao, D.; Duan, J.; Wang, W.; Yan, L. Iron and sulfur oxidation pathways of Acidithiobacillus ferrooxidans. World J. Microbiol. Biotechnol. 2019, 35, 60. [Google Scholar] [CrossRef] [PubMed]
  23. Boyd, E.S.; Druschel, G.K. Involvement of intermediate sulfur species in biological reduction of elemental sulfur under acidic, hydrothermal conditions. Appl. Environ. Microbiol. 2013, 79, 2061–2068. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Florentino, A.P.; Stams, A.J.; Sánchez-Andrea, I. Genome sequence of Desulfurella amilsii strain TR1 and comparative genomics of Desulfurellaceae family. Front. Microbiol. 2017, 8, 222. [Google Scholar] [CrossRef] [Green Version]
  25. Florentino, A.P.; Pereira, I.A.; Boeren, S.; van den Born, M.; Stams, A.J.; SánchezAndrea, I. Insight into the sulfur metabolism of Desulfurella amilsii by differential proteomics. Environ. Microbiol. 2019, 21, 209–225. [Google Scholar] [CrossRef] [Green Version]
  26. Yamamoto, M.; Takai, K. Sulfur metabolisms in epsilon- and gamma-proteobacteria in deep-sea hydrothermal fields. Front. Microbiol. 2011, 2, 192. [Google Scholar] [CrossRef] [Green Version]
  27. Wang, R.; Lin, J.Q.; Liu, X.M.; Pang, X.; Zhang, C.J.; Yang, C.L.; Gao, X.Y.; Lin, C.M.; Li, Y.Q.; Li, Y. Sulfur oxidation in the acidophilic autotrophic Acidithiobacillus spp. Front. Microbiol. 2019, 9, 3290. [Google Scholar] [CrossRef]
  28. Zhang, J.; Liu, R.; Xi, S.; Cai, R.; Sun, C. A novel bacterial thiosulfate oxidation pathway provides a new clue about the formation of zero-valent sulfur in deep sea. ISME J. 2020, 14, 2261–2274. [Google Scholar] [CrossRef]
  29. Yin, H.; Zhang, X.; Li, X.; He, Z.; Liang, Y.; Guo, X.; Hu, Q.; Xiao, Y.; Cong, J.; Ma, L. Whole-genome sequencing reveals novel insights into sulfur oxidation in the extremophile Acidithiobacillus thiooxidans. BMC Microbiol. 2014, 14, 179. [Google Scholar] [CrossRef] [Green Version]
  30. Friedrich, C.G.; Rother, D.; Bardischewsky, F.; Quentmeier, A.; Fischer, J. Oxidation of reduced inorganic sulfur compounds by bacteria: Emergence of a common mechanism? Appl. Environ. Microbiol. 2001, 67, 2873–2882. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Grabarczyk, D.B.; Berks, B.C. Intermediates in the Sox sulfur oxidation pathway are bound to a sulfane conjugate of the carrier protein SoxYZ. PloS ONE 2017, 12, e0173395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Ma, K.; Adams, M. Sulfide dehydrogenase from the hyperthermophilic archaeon Pyrococcus furiosus: A new multifunctional enzyme involved in the reduction of elemental sulfur. J. Bacteriol. 1994, 176, 6509–6517. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Hedderich, R.; Klimmek, O.; Kröger, A.; Dirmeier, R.; Keller, M.; Stetter, K.O. Anaerobic respiration with elemental sulfur and with disulfides. FEMS Microbiol. Rev. 1998, 22, 353–381. [Google Scholar] [CrossRef]
  34. Wang, S.; Jiang, L.; Hu, Q.; Liu, X.; Shao, Z. Elemental sulfur reduction by a deep-sea hydrothermal vent campylobacterium Sulfurimonas sp. NW10. Environ. Microbiol. 2021, 23, 965–979. [Google Scholar] [CrossRef]
  35. Takai, K.; Inagaki, F.; Nakagawa, S.; Hirayama, H.; Nunoura, T.; Sako, Y.; Nealson, K.H.; Horikoshi, K. Isolation and phylogenetic diversity of members of previously uncultivated ε-Proteobacteria in deep-sea hydrothermal fields. FEMS Microbiol. Lett. 2003, 218, 167–174. [Google Scholar] [PubMed] [Green Version]
  36. Vetriani, C.; Voordeckers, J.W.; Crespo-Medina, M.; O’Brien, C.E.; Giovannelli, D.; Lutz, R.A. Deep-sea hydrothermal vent Epsilonproteobacteria encode a conserved and widespread nitrate reduction pathway (Nap). ISME J. 2014, 8, 1510–1521. [Google Scholar] [CrossRef]
  37. Dick, G.J. The microbiomes of deep-sea hydrothermal vents: Distributed globally, shaped locally. Nat. Rev. Microbiol. 2019, 17, 271–283. [Google Scholar] [CrossRef]
  38. Inagaki, F.; Takai, K.; Nealson, K.H.; Horikoshi, K. Sulfurovum lithotrophicum gen. nov., sp. nov., a novel sulfur-oxidizing chemolithoautotroph within the ε-Proteobacteria isolated from Okinawa trough hydrothermal sediments. Int. J. Syst. Evol. Microbiol. 2004, 54, 1477–1482. [Google Scholar] [CrossRef] [Green Version]
  39. Reveillaud, J.; Reddington, E.; McDermott, J.; Algar, C.; Meyer, J.L.; Sylva, S.; Seewald, J.; German, C.R.; Huber, J.A. Subseafloor microbial communities in hydrogen-rich vent fluids from hydrothermal systems along the Mid-Cayman Rise. Environ. Microbiol. 2016, 18, 1970–1987. [Google Scholar] [CrossRef] [Green Version]
  40. Anderson, R.E.; Reveillaud, J.; Reddington, E.; Delmont, T.O.; Eren, A.M.; McDermott, J.M.; Seewald, J.S.; Huber, J.A. Genomic variation in microbial populations inhabiting the marine subseafloor at deep-sea hydrothermal vents. Nat. Commun. 2017, 8, 1114. [Google Scholar] [CrossRef] [Green Version]
  41. Watsuji, T.O.; Nakagawa, S.; Tsuchida, S.; Toki, T.; Hirota, A.; Tsunogai, U.; Takai, K. Diversity and function of epibiotic microbial communities on the galatheid crab, Shinkaia crosnieri. Microbes. Environ. 2010, 25, 288–294. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Pjevac, P.; Kamyshny, A.; Dyksma, S.; Mußmann, M. Microbial consumption of zero-valence sulfur in marine benthic habitats. Environ. Microbiol. 2014, 16, 3416–3430. [Google Scholar] [CrossRef] [PubMed]
  43. Xie, S.; Wang, S.; Li, D.; Shao, Z.; Lai, Q.; Wang, Y.; Wei, M.; Han, X.; Jiang, L. Sulfurovum indicum sp. nov. a novel hydrogen- and sulfur-oxidizing chemolithoautotroph isolated from a deep-sea hydrothermal plume in the northwestern indian ocean. Int. J. Syst. Evol. Microbiol. 2021, 71, 004748. [Google Scholar] [CrossRef] [PubMed]
  44. Pfaffl, M.W. A new mathematical model for relative quantification in realtime RT-PCR. Nucleic Acids Res. 2001, 29, e45. [Google Scholar] [CrossRef]
  45. Lahme, S.; Enning, D.; Callbeck, C.M.; Vega, D.M.; Curtis, T.P.; Head, I.M.; Hubert, C.R.J. Metabolites of an oil field sulfide-oxidizing, nitrate-reducing Sulfurimonas sp. cause severe corrosion. Appl. Environ. Microbiol. 2019, 85, e01891-18. [Google Scholar] [CrossRef] [Green Version]
  46. Callbeck, C.M.; Lavik, G.; Ferdelman, T.G.; Fuchs, B.; Gruber-Vodicka, H.R.; Hach, P.F.; Littmann, S.; Schoffelen, N.J.; Kalvelage, T.; Thomsen, S.; et al. Oxygen minimum zone cryptic sulfur cycling sustained by offshore transport of key sulfur oxidizing bacteria. Nat. Commun. 2018, 9, 1729. [Google Scholar] [CrossRef] [Green Version]
  47. Boyd, E.S.; Jackson, R.A.; Encarnacion, G.; Zahn, J.A.; Beard, T.; Leavitt, W.D.; Pi, Y.; Zhang, C.L.; Pearson, A.; Geesey, G.G. Isolation, characterization, and ecology of sulfur-respiring crenarchaea inhabiting acid-sulfate-chloride-containing geothermal springs in Yellowstone National Park. Appl. Environ. Microbiol. 2007, 73, 6669–6677. [Google Scholar] [CrossRef] [Green Version]
  48. Ma, N.; Sha, Z.L.; Sun, C.M. Formation of cadmium sulfide nanoparticles mediates cadmium resistance and light utilization of the deep-sea bacterium Idiomarina sp. OT37-5b. Environ. Microbiol. 2021, 23, 934–948. [Google Scholar] [CrossRef]
  49. Amenabar, M.J.; Boyd, E.S. Mechanisms of mineral substrate acquisition in a thermoacidophile. Appl. Environ. Microbiol. 2018, 84, e00334-18. [Google Scholar] [CrossRef] [Green Version]
  50. Thompson, J.D.; Gibson, T.J.; Higgins, D.G. Multiple sequence alignment using ClustalW and ClustalX. Curr. Protoc. Bioinform. 2002, 2, 3. [Google Scholar] [CrossRef]
  51. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Boratyn, G.M.; Schäffer, A.A.; Agarwala, R.; Altschul, S.F.; Lipman, D.J.; Madden, T.L. Domain enhanced lookup time accelerated BLAST. Biol. Direct 2012, 7, 12. [Google Scholar] [CrossRef] [Green Version]
  53. Slobodkina, G.B.; Mardanov, A.V.; Ravin, N.V.; Frolova, A.A.; Chernyh, N.A.; Bonch-Osmolovskaya, E.A.; Slobodkin, A.I. Respiratory ammonification of nitrate coupled to anaerobic oxidation of elemental sulfur in deep-sea autotrophic thermophilic bacteria. Front. Microbiol. 2017, 8, 87. [Google Scholar] [CrossRef] [PubMed]
  54. Flemming, H.C.; Wingender, J. The biofilm matrix. Nat. Rev. Microbiol. 2010, 8, 623–633. [Google Scholar] [CrossRef] [PubMed]
  55. Bogino, P.C.; Oliva, M.d.L.; Sorroche, F.G.; Giordano, W. The role of bacterial biofilms and surface components in plant-cacterial associations. Int. J. Mol. Sci. 2013, 14, 15838–15859. [Google Scholar] [CrossRef] [Green Version]
  56. Vozza, N.F.; Abdian, P.L.; Russo, D.M.; Mongiardini, E.J.; Lodeiro, A.R.; Molin, S.; Zorreguieta, A. A rhizobium leguminosarum CHDL- (cadherin-like-) lectin participates in assembly and remodeling of the biofilm matrix. Front. Microbiol. 2016, 7, 1608. [Google Scholar] [CrossRef] [Green Version]
  57. Craig, L.; Forest, K.T.; Maier, B. Type IV Pili: Dynamics, biophysics and functional consequences. Nat. Rev. Microbiol. 2019, 17, 429–440. [Google Scholar] [CrossRef]
  58. Sandeep, T.; Hancock, R.E.W. Involvement of two related porins, OprD and OpdP, in the uptake of arginine by Pseudomonas aeruginosa. FEMS Microbiol. Lett. 2006, 260, 23–29. [Google Scholar]
  59. Tommassen, J. Assembly of outer-membrane proteins in bacteria and mitochondria. Microbiology. 2010, 156, 2587–2596. [Google Scholar] [CrossRef] [Green Version]
  60. Götz, F.; Pjevac, P.; Markert, S.; McNichol, J.; Becher, D.; Schweder, T.; Mussmann, M.; Sievert, S.M. Transcriptomic and proteomic insight into the mechanism of cyclooctasulfur- versus thiosulfate-oxidation by the chemolithoautotroph Sulfurimonas denitrificans. Environ. Microbiol. 2019, 21, 244–258. [Google Scholar] [CrossRef] [Green Version]
  61. Buonfiglio, V.; Polidoro, M.; Flora, L.; Citro, G.; Valenti, P.; Orsi, N. Identification of two outer membrane proteins involved in the oxidation of sulphur compounds in Thiobacillus ferrooxidans. FEMS Microbiol. Rev. 1993, 11, 43–50. [Google Scholar] [CrossRef]
  62. Lahme, S.; Callbeck, C.M.; Eland, L.E.; Wipat, A.; Enning, D.; Head, I.M.; Hubert, C.R.J. Comparison of sulfide-oxidizing Sulfurimonas strains reveals a new mode of thiosulfate formation in subsurface environments. Environ. Microbiol. 2020, 22, 1784–1800. [Google Scholar] [CrossRef] [PubMed]
  63. Heide, V.D.T.; Poolman, B. ABC transporters: One, two or four extracytoplasmic substrate-binding sites? EMBO Rep. 2002, 3, 938–943. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Mulligan, C.; Fischer, M.; Thomas, G.H. Tripartite ATP-independent periplasmic (TRAP) transporters in bacteria and archaea. FEMS Microbiol. Rev. 2011, 35, 68–86. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Bagchi, A.; Ghosh, T.C. Structural insight into the interactions of SoxV, SoxW and SoxS in the process of transport of reductants during sulfur oxidation by the novel global sulfur oxidation reaction cycle. Biophys. Chem. 2006, 119, 7–13. [Google Scholar] [CrossRef] [PubMed]
  66. Qi, H.; Chen, H.; Ao, J.; Zhou, H.; Chen, X. Identification of differentially expressed genes in Sulfobacillus sp. TPY grown on either elemental sulphur or Fe2+. J. Gen. Appl. Microbiol. 2010, 56, 389–397. [Google Scholar] [CrossRef] [Green Version]
  67. Wang, S.; Jiang, L.; Hu, Q.; Cui, L.; Zhu, B.; Fu, X.; Lai, Q.; Shao, Z.; Yang, S. Characterization of Sulfurimonas hydrogeniphila sp. nov., a novel bacterium predominant in deep-sea hydrothermal vents and comparative genomic analyses of the genus Sulfurimonas. Front. Microbiol. 2021, 12, 626705. [Google Scholar] [CrossRef]
  68. Yamamoto, M.; Nakagawa, S.; Shimamura, S.; Takai, K.; Horikoshi, K. Molecular characterization of inorganic sulfur-compound metabolism in the deep-sea epsilonproteobacterium Sulfurovum sp. NBC37-1. Environ. Microbiol. 2010, 12, 1144–1153. [Google Scholar] [CrossRef]
  69. Jelen, B.; Giovannelli, D.; Falkowski, P.G.; Vetriani, C. Elemental sulfur reduction in the deep-sea vent thermophile, Thermovibrio ammonificans. Environ. Microbiol. 2018, 20, 2301–2316. [Google Scholar] [CrossRef]
  70. Krafft, T.; Bokranz, M.; Klimmek, O.; Schröder, I.; Fahrenholz, F.; Kojro, E.; Kröger, A. Cloning and nucleotide sequence of the psrA gene of Wolinella succinogenes polysulphide reductase. Eur. J. Biochem. 1992, 206, 503–510. [Google Scholar] [CrossRef]
  71. Hinsley, A.P.; Berks, B.C. Specificity of respiratory pathways involved in the reduction of sulfur compounds by Salmonella enterica. Microbiology 2002, 148, 3631–3638. [Google Scholar] [CrossRef] [Green Version]
  72. Laska, S.; Lottspeich, F.; Kletzin, A. Membrane-bound hydrogenase and sulfur reductase of the hyperthermophilic and acidophilic archaeon Acidianus ambivalens. Microbiology 2003, 149, 2357–2371. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Guiral, M.; Tron, P.; Aubert, C.; Gloter, A.; Iobbi-Nivol, C.; Giudici-Orticoni, M.T. A membrane-bound multienzyme, hydrogen-oxidizing, and sulfur-reducing complex from the hyperthermophilic bacterium Aquifex aeolicus. J. Biol. Chem. 2005, 280, 42004–42015. [Google Scholar] [CrossRef] [Green Version]
  74. Takai, K.; Campbell, B.J.; Cary, S.C.; Suzuki, M.; Oida, H.; Nunoura, T.; Hirayama, H.; Nakagawa, S.; Suzuki, Y.; Inagaki, F.; et al. Enzymatic and genetic characterization of carbon and energy metabolisms by deep-sea hydrothermal chemolithoautotrophic isolates of Epsilonproteobacteria. Appl. Environ. Microbiol. 2005, 71, 7310–7320. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Wang, S.; Jiang, L.; Xie, S.; Alain, K.; Wang, Z.; Wang, J.; Liu, D.; Shao, Z. Disproportionation of inorganic sulfur compounds by mesophilic chemolithoautotrophic Campylobacterota. Msystems 2022, 21, e0095422. [Google Scholar] [CrossRef] [PubMed]
  76. Franz, B.; Lichtenberg, H.; Hormes, J.; Modrow, H.; Dahl, C.; Prange, A. Utilization of solid ‘elemental’ sulfur by the phototrophic purple sulfur bacterium Allochromatium vinosum: A sulfur K-edge X-ray absorption spectroscopy study. Microbiology 2007, 153, 1268–1274. [Google Scholar] [CrossRef] [Green Version]
  77. Blumentals, I.I.; Itoh, M.; Olson, G.J.; Kelly, R.M. Role of polysulfides in reduction of elemental sulfur by the hyperthermophilic archaebacterium Pyrococcus furiosus. Appl. Environ. Microbiol. 1990, 56, 1255–1262. [Google Scholar] [CrossRef] [Green Version]
  78. Wu, X.; Pedersen, K.; Edlund, J.; Eriksson, L.; Åström, M.; Andersson, A.F.; Bertilsson, S.; Dopson, M. Potential for hydrogen-oxidizing chemolithoautotrophic and diazotrophic populations to initiate biofilm formation in oligotrophic, deep terrestrial subsurface waters. Microbiome 2017, 5, 37. [Google Scholar] [CrossRef] [Green Version]
  79. Kerou, M.; Offre, P.; Valledor, L.; Abby, S.S.; Melcher, M.; Nagler, M.; Weckwerth, W.; Schleper, C. Proteomics and comparative genomics of Nitrososphaera viennensis reveal the core genome and adaptations of archaeal ammonia oxidizers. Proc. Natl. Acad. Sci. USA 2016, 113, E7937–E7946. [Google Scholar] [CrossRef] [Green Version]
  80. Zhang, R.; Bellenberg, S.; Neu, T.R.; Sand, W.; Vera, M. The Biofilm Lifestyle of Acidophilic Metal/Sulfur-Oxidizing Microorganisms. In Biotechnology of Extremophiles: Grand Challenges in Biology and Biotechnology; Rampelotto, P., Ed.; Springer: Cham, Switzerland, 2016; Volume 1. [Google Scholar]
  81. Zhang, R.; Neu, T.R.; Li, Q.; Blanchard, V.; Zhang, Y.; Schippers, A.; Sand, W. Insight into interactions of thermoacidophilic archaea with elemental sulfur: Biofilm dynamics and EPS analysis. Front. Microbiol. 2019, 10, 896. [Google Scholar] [CrossRef] [Green Version]
  82. Svensson, S.L.; Pryjma, M.; Gaynor, E.C. Flagella-mediated adhesion and extracellular DNA release contribute to biofilm formation and stress tolerance of Campylobacter jejuni. PLoS ONE 2014, 9, e106063. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Jarrell, K.F.; Stark, M.; Nair, D.B.; Chong, J.P.J. Flagella and pili are both necessary for efficient attachment of Methanococcus maripaludis to surfaces. FEMS Microbiol. Lett. 2011, 319, 44–50. [Google Scholar] [CrossRef] [PubMed]
  84. Harneit, K.G.; Ksel, A.; Kock, D.; Klock, J.H.; Gehrke, T.; Sand, W. Adhesion to metal sulfide surfaces by cells of Acidithiobacillus ferrooxidans, Acidithiobacillus thiooxidans and Leptospirillum ferrooxidans. Hydrometallurgy 2006, 83, 245–254. [Google Scholar] [CrossRef]
  85. Zhu, H.; Tian, B.; Liu, W.; Zhang, S.; Cao, C.; Zhang, Y.; Zou, W. A three-stage culture process for improved exopolysaccharide production by tremella fuciformis. Bioresour. Technol. 2012, 116, 526–528. [Google Scholar] [CrossRef]
  86. Manchur, M.A.; Kikumoto, M.; Kanao, T.; Takada, J.; Kamimura, K. Characterization of an OmpA-like outer membrane protein of the acidophilic iron-oxidizing bacterium, Acidithiobacillus ferrooxidans. Extremophiles 2011, 15, 403–410. [Google Scholar] [CrossRef] [Green Version]
  87. Ramírez, P.; Guiliani, N.; Valenzuela, L.; Beard, S.; Jerez, C.A. Differential protein expression during growth of Acidithiobacillus ferrooxidans on ferrous iron, sulfur compounds, or metal sulfides. Appl. Environ. Microbiol. 2004, 70, 4491–4498. [Google Scholar] [CrossRef] [Green Version]
  88. Ma, N.; Sun, C. Cadmium sulfide nanoparticle biomineralization and biofilm formation mediate cadmium resistance of the deep-sea bacterium Pseudoalteromonas sp. MT33b. Environ. Microbiol. Rep. 2021, 13, 325–336. [Google Scholar] [CrossRef] [PubMed]
  89. Ares, A.; Sakai, S.; Sasaki, T.; Mitarai, S.; Nunoura, T. Sequestration and efflux largely account for Cd and Cu resistance in the deep sea Epsilonproteobacterium, Nitratiruptor sp. SB155-2. Environ. Microbiol. 2021, 16255, 1462–2920. [Google Scholar]
  90. Quatrini, R.; Jedlicki, E.; Holmes, D.S. Genomic insights into the iron uptake mechanisms of the biomining microorganism Acidithiobacillus ferrooxidans. J. Ind. Microbiol. Biotechnol. 2005, 32, 606–614. [Google Scholar] [CrossRef] [PubMed]
  91. Hensen, D.; Sperling, D.; Trüper, H.G.; Brune, D.C.; Dahl, C. Thiosulphate oxidation in the phototrophic sulphur bacterium Allochromatium vinosum. Mol. Microbiol. 2006, 62, 794–810. [Google Scholar] [CrossRef]
  92. Pjevac, P.; Meier, D.V.; Markert, S.; Hentschker, C.; Schweder, T.; Becher, D.; Gruber-Vodicka, H.R.; Richter, M.; Bach, W.; Amann, R.; et al. Metaproteogenomic profiling of microbial communities colonizing actively venting hydrothermal chimneys. Front. Microbiol. 2018, 9, 680. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Chemolithoautotrophic growth of S. indicum ST-419 in the presence of S0 as the sole electron donor (A) or the sole electron acceptor (C). Cultures of strain ST-419 by thiosulfate oxidation (B) and nitrate reduction (D) served as controls, respectively. Cell density in S0 reduction and nitrate reduction is shown in (E,F). Error bars indicate the standard deviation between triplicate cultures. Closed squares, diamonds, circles, triangles and inverted triangles indicate the concentrations of NO3, SO42−, S2O32−, N2 and H2. Open squares and circles indicate H2S concentration and cell numbers, respectively.
Figure 1. Chemolithoautotrophic growth of S. indicum ST-419 in the presence of S0 as the sole electron donor (A) or the sole electron acceptor (C). Cultures of strain ST-419 by thiosulfate oxidation (B) and nitrate reduction (D) served as controls, respectively. Cell density in S0 reduction and nitrate reduction is shown in (E,F). Error bars indicate the standard deviation between triplicate cultures. Closed squares, diamonds, circles, triangles and inverted triangles indicate the concentrations of NO3, SO42−, S2O32−, N2 and H2. Open squares and circles indicate H2S concentration and cell numbers, respectively.
Antioxidants 12 00627 g001
Figure 2. Volcano plots showing differentially expressed genes of S. indicum ST-419 under S0 oxidation (A) and S0 reduction (B) compared to the control. (C) Venn diagrams showing the differential expression gene occurring in both conditions including upregulated (Left) and downregulated (Right). % represents the proportion of differential expression genes shared in the two conditions among the CDS of the genome.
Figure 2. Volcano plots showing differentially expressed genes of S. indicum ST-419 under S0 oxidation (A) and S0 reduction (B) compared to the control. (C) Venn diagrams showing the differential expression gene occurring in both conditions including upregulated (Left) and downregulated (Right). % represents the proportion of differential expression genes shared in the two conditions among the CDS of the genome.
Antioxidants 12 00627 g002
Figure 3. Transcriptomic analysis of genes associated with biofilm formation of strain ST-419 between S0 oxidation and reduction. Heat map analysis of differentially expressed genes related to exopolysaccharide synthesis (A), bacteria secretion (B), TonB-dependent receptors (C) and signal transduction (D). S0_O represents S0 oxidation; Na2S2O3_O represents thiosulfate oxidation; S0_N represents S0 reduction and NaNO3_N represents nitrate reduction.
Figure 3. Transcriptomic analysis of genes associated with biofilm formation of strain ST-419 between S0 oxidation and reduction. Heat map analysis of differentially expressed genes related to exopolysaccharide synthesis (A), bacteria secretion (B), TonB-dependent receptors (C) and signal transduction (D). S0_O represents S0 oxidation; Na2S2O3_O represents thiosulfate oxidation; S0_N represents S0 reduction and NaNO3_N represents nitrate reduction.
Antioxidants 12 00627 g003
Figure 4. Transcriptomic analysis of genes associated with outer membrane proteins of strain ST-419 under S0 oxidation and reduction (A). Phylogenetic relationships of outer membrane protein from strain ST-419 with homologous genes in the genera of Sulfurimonas, Acinetobacter and Pseudomonas (B). The three OprD proteins in red fonts have the highest expressions under S0 oxidation and reduction. The OprD from Pseudomonas aeruginosa PAO1 in bold is the known characterized enzyme.
Figure 4. Transcriptomic analysis of genes associated with outer membrane proteins of strain ST-419 under S0 oxidation and reduction (A). Phylogenetic relationships of outer membrane protein from strain ST-419 with homologous genes in the genera of Sulfurimonas, Acinetobacter and Pseudomonas (B). The three OprD proteins in red fonts have the highest expressions under S0 oxidation and reduction. The OprD from Pseudomonas aeruginosa PAO1 in bold is the known characterized enzyme.
Antioxidants 12 00627 g004
Figure 5. Sulfide concentrations (A) and cell densities (B) in cultures of strain ST-419 grown chemolithoautotrophically with H2 as the sole electron donor and S0 as the sole electron acceptor. S0 was provided in the bulk medium (control) or was sequestered in dialysis membranes (pore sizes of 6 to 8 kDa or 12 to 14 kDa) to prevent physical contact with the bulk S0.
Figure 5. Sulfide concentrations (A) and cell densities (B) in cultures of strain ST-419 grown chemolithoautotrophically with H2 as the sole electron donor and S0 as the sole electron acceptor. S0 was provided in the bulk medium (control) or was sequestered in dialysis membranes (pore sizes of 6 to 8 kDa or 12 to 14 kDa) to prevent physical contact with the bulk S0.
Antioxidants 12 00627 g005
Figure 6. The relative abundance and expression of genes related to S0 oxidation (A) and S0 reduction (B) in strain ST-419. Light green bars show the average values of strain ST-419 grown on S0 oxidation and reduction, and light orange bars show the average of cultures grown on thiosulfate oxidation and nitrate reduction. Statistically significant differences are denoted by one asterisk (padj < 0.05) or two asterisks (padj < 0.01). The symbols of “I”, “Ⅱ” and “Ⅲ” represente the different gene cluster of polysulfide reductases. Normalized gene expression is plotted as fragments per kilobase and million reads (FPKM). Error bars indicate the standard deviation from triplicate cultures.
Figure 6. The relative abundance and expression of genes related to S0 oxidation (A) and S0 reduction (B) in strain ST-419. Light green bars show the average values of strain ST-419 grown on S0 oxidation and reduction, and light orange bars show the average of cultures grown on thiosulfate oxidation and nitrate reduction. Statistically significant differences are denoted by one asterisk (padj < 0.05) or two asterisks (padj < 0.01). The symbols of “I”, “Ⅱ” and “Ⅲ” represente the different gene cluster of polysulfide reductases. Normalized gene expression is plotted as fragments per kilobase and million reads (FPKM). Error bars indicate the standard deviation from triplicate cultures.
Antioxidants 12 00627 g006
Figure 7. Phylogenetic relationships of putative PsrAs from strain ST-419 with various sequences of the sulfur reductase/polysulfide reductase/thiosulfate reductase family. The polysulfide reductase present in cultured strains of the genus Sulfurovum are shown in red and the homologous genes in strain ST-419 are marked with an asterisk. Biochemically characterized enzymes are shown in bold. The red stars represented the polysulfide reductases of strain ST-419. Amino acid sequences were derived from the non-redundant protein database of NCBI and accession numbers are shown in the rear. Bootstrap values based on 1000 replicates are shown at branch nodes. Bar = 1.0 substitutions per protein position.
Figure 7. Phylogenetic relationships of putative PsrAs from strain ST-419 with various sequences of the sulfur reductase/polysulfide reductase/thiosulfate reductase family. The polysulfide reductase present in cultured strains of the genus Sulfurovum are shown in red and the homologous genes in strain ST-419 are marked with an asterisk. Biochemically characterized enzymes are shown in bold. The red stars represented the polysulfide reductases of strain ST-419. Amino acid sequences were derived from the non-redundant protein database of NCBI and accession numbers are shown in the rear. Bootstrap values based on 1000 replicates are shown at branch nodes. Bar = 1.0 substitutions per protein position.
Antioxidants 12 00627 g007
Figure 8. Proposed model for extracellular S0 activation, transportation and intracellular oxidation and reduction in strain ST-419. OM, outer membrane; PS, periplasmic space; IM, inner membrane; CS, cytoplasmic space; Sox, sulfur-oxidation-multienzyme complex; Q, quinone molecules; Psr, polysulfide reductase; TC, two-component system; TBDT, TonB-dependent transporters. Hypothetical reactions concerning energy conservation are indicated by a question mark.
Figure 8. Proposed model for extracellular S0 activation, transportation and intracellular oxidation and reduction in strain ST-419. OM, outer membrane; PS, periplasmic space; IM, inner membrane; CS, cytoplasmic space; Sox, sulfur-oxidation-multienzyme complex; Q, quinone molecules; Psr, polysulfide reductase; TC, two-component system; TBDT, TonB-dependent transporters. Hypothetical reactions concerning energy conservation are indicated by a question mark.
Antioxidants 12 00627 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, S.; Jiang, L.; Cui, L.; Alain, K.; Xie, S.; Shao, Z. Transcriptome Analysis of Cyclooctasulfur Oxidation and Reduction by the Neutrophilic Chemolithoautotrophic Sulfurovum indicum from Deep-Sea Hydrothermal Ecosystems. Antioxidants 2023, 12, 627. https://doi.org/10.3390/antiox12030627

AMA Style

Wang S, Jiang L, Cui L, Alain K, Xie S, Shao Z. Transcriptome Analysis of Cyclooctasulfur Oxidation and Reduction by the Neutrophilic Chemolithoautotrophic Sulfurovum indicum from Deep-Sea Hydrothermal Ecosystems. Antioxidants. 2023; 12(3):627. https://doi.org/10.3390/antiox12030627

Chicago/Turabian Style

Wang, Shasha, Lijing Jiang, Liang Cui, Karine Alain, Shaobin Xie, and Zongze Shao. 2023. "Transcriptome Analysis of Cyclooctasulfur Oxidation and Reduction by the Neutrophilic Chemolithoautotrophic Sulfurovum indicum from Deep-Sea Hydrothermal Ecosystems" Antioxidants 12, no. 3: 627. https://doi.org/10.3390/antiox12030627

APA Style

Wang, S., Jiang, L., Cui, L., Alain, K., Xie, S., & Shao, Z. (2023). Transcriptome Analysis of Cyclooctasulfur Oxidation and Reduction by the Neutrophilic Chemolithoautotrophic Sulfurovum indicum from Deep-Sea Hydrothermal Ecosystems. Antioxidants, 12(3), 627. https://doi.org/10.3390/antiox12030627

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop