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

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.


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, deep-sea hydrothermal vents. Furthermore, the isolates of Sulfurovum identified to date generally use H 2 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 S 0 as the sole electron donor (S 0 oxidation), as well as the terminal electron acceptor (S 0 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 S 0 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.

Bacterial Strains and Growth Conditions
Cultures of strain ST-419 were routinely grown in a modified MMJS medium at 37°C 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% N 2 /20% CO 2 (200 kPa). For sulfur reduction cultures, S 0 (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% H 2 /20% CO 2 (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% H 2 /20% CO 2 (200 kPa). S 0 granules were heattreated 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 S 0 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: µ = (lnN 2 − lnN 1 )/(t 2 − t 1 ), where N 2 and N 1 represent the number of cells mL −1 at time (hours of incubation) t 2 and t 1 , 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).

Chemical Analyses
For direct analysis, samples were taken from cultures with N 2 -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.

Transcriptome Analysis
Cells of strain ST-419 cultivated with different electron donors and acceptors under four types of catabolic conditions including S 0 oxidation, thiosulfate oxidation, S 0 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 Supple-

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.

Experiment with Dialysis Membranes to Test Cell Contact with S 0 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 S 0 oxidation and reduction processes. S 0 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.

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).

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).

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 S 0 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 S 0 had a significantly longer lag phase, which is necessary for cells to activate solid S 0 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: S 0 + 1.2 NO 3 -+ 0.4 H 2 O→SO 4 2-+ 0.6 N 2 + 0.8 H + [53]. Compared to S 0 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 S 0 (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.

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, S 0 , which initially floated on the medium, gradually sank to the bottom of the bottle as cell growth and sulfur reduction began. S 0 -reducing cultures consistently grew more slowly than nitrate-reducing cultures, but similar cell densities could be achieved after 24 h (>9 × 10 7 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.

General Features of Transcriptomes for S 0 Oxidation and Reduction Conditions
To examine changes in gene expression specifically associated with the oxidation and reduction of S 0 , we investigated the transcriptomes of strain ST-419 using cells harvested at four different conditions in two pairs: S 0 oxidation vs. thiosulfate oxidation (S 0 oxidation), and S 0 reduction vs. nitrate reduction (S 0 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 S 0 oxidation and reduction transcriptomes were evaluated using a heat map, which revealed that the expression profiles were separated among S 0 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 S 0 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 S 0 and nitrate, respectively (97.71-98.22% of protein coding genes), with 533 and 607 significantly increased and decreased, respectively, during reduction of S 0 compared with nitrate reduction ( Figure 2B). The list of the transcripts significantly expressed under both S 0 oxidation and reduction is given in Table S2. Among the significant upregulation transcripts, 144 genes were shared between S 0 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 S 0 oxidation and reduction, corresponding to 7.1% of the total CDS ( Figure 2C). The transcripts significantly expressed under both S 0 oxidation and reduction in common are shown in Table S2.

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, S 0 , which initially floated on the medium, gradually sank to the bottom of the bottle as cell growth and sulfur reduction began. S 0 -reducing cultures consistently grew more slowly of the transcripts significantly expressed under both S 0 oxidation and reduction is given in Table S2. Among the significant upregulation transcripts, 144 genes were shared between S 0 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 S 0 oxidation and reduction, corresponding to 7.1% of the total CDS ( Figure 2C). The transcripts significantly expressed under both S 0 oxidation and reduction in common are shown in Table S2.

Genes Potentially Involved in S 0 Activation and Their Differential Expression
Among the 144 genes significantly overexpressed in both S 0 oxidation and S 0 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 S 0 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 S 0 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 S 0 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 S 0 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 S 0 either oxidation or reduction ( Figure 3C,D). sponsible 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 S 0 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 S 0 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 S 0 either oxidation or reduction ( Figure 3C,D).  , TonB-dependent receptors (C) and signal transduction (D). S 0 _O represents S 0 oxidation; Na 2 S 2 O 3 _O represents thiosulfate oxidation; S 0 _N represents S 0 reduction and NaNO 3 _N represents nitrate reduction.

Genes Potentially Involved in S 0 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 S 0 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 S 0 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 S 0 oxidation and reduction, respectively ( Figure 4A; Table S2). Similarly, the other gene (IMZ28_RS01320) had the highest expression under S 0 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 S 0 [60,62]. Additionally, several genes coding for transporters were differentially expressed under both S 0 oxidation and S 0 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 S 0 by oxidative phosphorylation, when electrons released by catabolism (electrons derived from H2 in sulfur reduction, and from S 0 /H2S in sulfide oxidation) enter the respiratory chain.

Evaluation of the Necessity for Direct Cell-S 0 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 S 0 was not required for bacterial growth and sulfide production during sulfur reduction ( Figure 5). Compared to the Additionally, several genes coding for transporters were differentially expressed under both S 0 oxidation and S 0 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 ATPindependent 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 S 0 by oxidative phosphorylation, when electrons released by catabolism (electrons derived from H 2 in sulfur reduction, and from S 0 /H 2 S in sulfide oxidation) enter the respiratory chain.

Evaluation of the Necessity for Direct Cell-S 0 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 S 0 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 S 0 oxidation treatments when S 0 was sequestered in dialysis membranes. Thus, cells of the bacterium growing via S 0 reduction did not require direct contact to S 0 granules, whereas growing via S 0 oxidation did require direct contact.
Antioxidants 2023, 12, x FOR PEER REVIEW 11 of 2 condition without the dialysis membrane, both sulfur reduction and bacterial growth de creased 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). Cor respondingly, the final cell yield decreased by 41% and 22%, respectively ( Figure 5B). In contrast, no cell growth was observed in S 0 oxidation treatments when S 0 was sequestered in dialysis membranes. Thus, cells of the bacterium growing via S 0 reduction did not re quire direct contact to S 0 granules, whereas growing via S 0 oxidation did require direc contact.

Expression of the Genes Involved in S 0 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 operon showed obvious different transcription levels between thiosulfate and S 0 oxidation condi

Expression of the Genes Involved in S 0 Oxidation
High transcripts of the periplasmic Sox multienzyme complex, encoded by two gene operon, soxABXY 1 Z 1 and soxCDZ 2 Y 2 , 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 S 0 oxidation conditions ( Figure 6A). The gene cluster soxABXY 1 Z 1 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 soxH 1 (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 soxCDY 2 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 S 0 oxidation compared to thiosulfate oxidation, with the exception of the gene soxZ 2 , which showed no significant variation under both culture conditions ( Figure 6A). In addition, a soxH 2 -like gene (IMZ28_RS10560) adjacent to the soxCDY 2 Z 2 gene cluster was expressed in raised abundance under S 0 oxidation ( Figure 6A). This indicates the important role of soxCDY 2 Z 2 H 2 in the oxidation of S 0 .
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 S 0 oxidation. The gene hdrB (IMZ28_RS05750) was significantly upregulated by 1.52-fold under S 0 oxidation ( Figure 6A). Structural domain analysis showed that the hdrB-encoding protein has the cysteine functional motif, which is thought to be important for S 0 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 S 0 oxidation ( Figure 6A), which may be involved in the electron transfer during S 0 oxidation. In addition, a gene (IMZ28_RS07900) encoding MBLfold metallo-hydrolase showed a low expression under S 0 oxidation ( Figure 6A), although recently it was proposed to be involved in sulfur oxidation [28].

Expression of the Genes Involved in S 0 Reduction
Based on genome analysis, strain ST-419 contained four types of putative polysulfide reductase-encoding genes (psrA 1 B 1 C 1 , psrA 2 B 2 , psrA 3 B 3 C 3 and psrA 4 B 4 C 4 ), which may perform a sulfur reduction function. These four kinds of polysulfide reductase-encoding genes showed differential transcriptional expressions under S 0 reduction compared to nitrate reduction ( Figure 6B). Among them, the gene cluster psrA 3 B 3 C 3 was significantly upregulated by 4.03-7.12-fold with the high FPKM value of 1986-8312 under S 0 reduction ( Figure 6B). Three genes in the gene cluster psrA 1 B 1 C 1 also exhibited a high FPKM value of 2504-17,874, corresponding to an increase in gene expressions by 1.21-2.10-fold during S 0 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 psrA 2 B 2 and psrA 4 B 4 C 4 were quite low, with gene psrA 2 showing a significantly upregulated expression and gene psrA 4 showing a downregulated expression ( Figure 6B). These results indicated that the polysulfide reductase encoded by psrA 3 B 3 C 3 may be more important for sulfur reduction. lase, 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 S 0 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 S 0 oxidation ( Figure 6A). This indicates the important role of soxCDY2Z2H2 in the oxidation of S 0 .  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 S 0 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 S 0 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 S 0 reduction ( Figure 6B).

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 PsrA 1 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 psrA 1 B 1 C 1 shared the highest sequence similarity of 81-89% with the known sulfur respiration enzymes PsrABC of Sulfurovum sp. NBC37-1 [67]. The PsrA 2 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 psrA 2 B 2 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 PsrA 3 and PsrA 4 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 psrA 3 B 3 C 3 and psrA 4 B 4 C 4 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 S 0 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, 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 S 0 to H 2 S 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 S 0 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.

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 S 0 utilization in the genus Sulfurovum.
Considering the extremely low solubility and reactivity of S 0 [76,77], microorganisms most likely need a specific activation or solubilization mechanism to make S 0 available for their energy metabolism. Biofilm formation is an important way for cells to adsorb S 0 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 S 0 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 S 0 (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 S 0 ( Figure 3C,D). Therefore, biofilm formation by Sulfurovum may be a crucial step in activating S 0 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 S 0 (Figure 4) suggest their involvement in utilizing this substrate. Two oprD-like porin genes (IMZ28_RS01315 and IMZ28_RS01320) were highly expressed under S 0 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 S 0 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 S 0 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 transporta-tion [17,86,87]. However, the true membrane proteins involved in S 0 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 S 0 degradation/activation, or whether they facilitate S 0 transport into cells.
In addition to the direct contact of cells to insoluble S 0 , a non-contact or cooperative mechanism (coexistence of contact and non-contact activation) could be involved in S 0 activation. In this study, we found that strain ST-419 growing via S 0 reduction did not require direct contact to S 0 , whereas cells growing via S 0 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 S 0 without cell contact also likely exists. Furthermore, we found some genes associated with cysteine and glutathione synthesis which were differentially upregulated under S 0 oxidation and reduction (Table S5). Under S 0 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 S 0 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 S 0 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 S 0 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 soxABXY 1 Z 1 gene cluster was more expressed under thiosulfate oxidation condition, while the soxCDY 2 Z 2 gene cluster was more expressed under S 0 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 S 0 oxidation in strain ST-419 is performed solely by the gene cluster soxCDY 2 Z 2 . Indeed, Sulfurimonas sp. CVO from an oil field is able to oxidize S 0 to thiosulfate and sulfate with only soxCDY 2 Z 2 in the absence of soxABXY 1 Z 1 [62]. Thus, we propose that the S 0 oxidation pathway is performed by Sulfurovum species as follows ( Figure 8). First, circular S 8 is transformed into linear polysulfide (S n -, HS n -) 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 SoxY 2 , and generates a thiocysteine-S sulfate residue (SoxZ 2 Y 2 -S-S 7 -S -). The outer sulfone sulfur of the cystein-persulfide on SoxY 2 is then oxidized by SoxCD to form a cysteine-S sulfate residue (SoxZ 2 Y 2 -S-S 7 -SO 2 -), and subsequently oxidized to form SoxZ 2 Y 2 -S-S 7 -SO 3 2-. Finally, this complex is hydrolyzed to sulfate by SoxH 2 or by another way, and regenerates the SoxZ 2 Y 2 complex (SoxZ 2 Y 2 -SH) (Figure 8). At present, we are not able to resolve how the sulfonate group bound to SoxY 2 Z 2 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 soxCDY 2 Z 2 and its homologs can be retrieved in many other sulfuroxidizing Campylobacterota with conserved synteny as a thiol hydrolase. Another possibility might be that SoxCDY 2 Z 2 can catalyze the reaction by itself, as previously supposed in S. denitrificans DSM1251 [60], though there is no relevant experimental evidence. tinuously release the compounds containing a -SH group to catalyze the transition of the S 0 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 S 0 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.  During sulfur reduction, four polysulfide reductases of three groups including group I (psrA 1 B 1 C 1 ), group II (psrA 2 B 2 ) and group III (psrA 3 B 3 C 3 and psrA 4 B 4 C 4 ) were differentially expressed ( Figures 6B and 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 psrA 3 B 3 C 3 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 psrA 1 B 1 CDE and psrA 2 B 2 , 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 H 2 S by PsrA 3 B 3 C 3 (Figure 8). The end product of H 2 S 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.

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 S 0 activation, which would involve the reducing compounds such as Cys, GSH and H 2 S. Transcriptomic data indicated that the complex soxCDY 2 Z 2 H 2 plays a role in S 0 oxidization in the periplasm, while among the four polysulfide reductases, psrA 3 B 3 C 3 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 S 0 activation, and to validate the tentatively proposed models for sulfur oxidation and reduction.  Table S1: Primers for qRT-PCR used in this study; Table S2: he differentially up-regulated and down-regulated transcripts in S 0 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 S 0 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 S 0 oxidation and reduction.