To obtain CSOB from a newly discovered hydrothermal vent on the Carlsberg Ridge, active vent chimney samples were ground on board and inoculated into sealed bottles filled with MMJHS medium with thiosulfate and hydrogen as the energy sources. After about one and half month incubation on board, enriched bacterial cultures were transferred into 10 ml MMJHS medium, and incubated at 28°C in the laboratory. After 2 days of incubation, bacterial growth was obvious with cells in form of short rods. The bacterium was subsequently purified with the dilution-to-extinction method. The culture in the serum bottle showing growth at the highest dilution was designated as strain NW10^T. The purity of the culture was further confirmed by microscopic examination and 16S rRNA gene sequencing. Interestingly, the bacterium can produce large amount of elemental sulfur in the form of extracellular granules as determined below.

Morphological observations by phase-contrast light microscopy and transmission electron microscope showed that cells of strain NW10^T were Gram-negative, motile and short rod-shaped with a size of 0.4–0.8 μm wide and 0.8–3.5 μm long (Supplementary Figure S1). Spore formation was not found during the culture incubation. Cells in older cultures tended to form aggregates. These morphological features are shared with other vent species S. autotrophica OK10^T and S. paralvinellae GO25^T (Table 1). In addition, on solid MMJHS medium agar plates, strain NW10^T formed small, white, round-shaped colonies with smooth boundaries.

Sequence comparison of the 16S rRNA gene sequences obtained from PCR amplification (1,447 bp) showed that strain NW10^T is closest to the only two deep-sea vent bacterial species S. paralvinellae GO25^T and S. autotrophica OK10^T, with 95.8 and 95.2% sequence identities, respectively. The maximum-likelihood phylogenetic tree based on 16S rRNA gene sequences showed that strain NW10^T clusters with S. paralvinellae GO25^T within the genus Sulfurimonas (Figure 1). This was further confirmed through the phylogenomic trees constructed with the neighbor-joining and minimum evolution methods (Supplementary Figures S2, S3). The phylogenomic tree based on the up-to-date bacterial core gene sequences showed that strain NW10^T forms a branch with S. autotrophica OK10^T and S. paralvinellae GO25^T (Figure 2), supporting further that strain NW10^T should belong to the Sulfurimonas genus and likely represents a novel species.

Growth experiment showed that strain NW10^T could grow in the range of temperatures (4–45°C), salinities (2–4% (w/v) NaCl), pH (5.0–9.0) and oxygen concentrations (1–20%). Strain NW10^T also could use nitrate as sole electron acceptor in the absence of oxygen. The optimal growth occurred at 33°C, 6% O₂, pH 6.0–6.5 and 3% (w/v) NaCl (Table 1). Under optimal growth conditions, the maximum cell concentration was approximately 5.0 × 10⁸ cells/ml and the doubling time was approximately 6 h. Chemoautotrophic growth showed that strain NW10^T could grow with thiosulfate, sulfide, elemental sulfur, and hydrogen as energy sources, but not with sulfite and tetrathionate. Strain NW10^T grew better with hydrogen as the sole energy source, which brought about the highest cell concentration far more than other energy sources. Thus, hydrogen is possibly the preferred energy source for this bacterium. Similarly, hydrogen also supported the best growth of S. paralvinellae GO25^T (Takai et al., 2006). When hydrogen was used as the energy source, strain NW10^T could respire element sulfur. The product of sulfur reduction was sulfide, which reached up to 42 μM at the late exponential phase in the medium (Wang et al., 2020a). As a chemolithoautotroph, strain NW10^T could not grow using any of the tested organic compounds as the carbon source. In addition, none of these organic compounds could be used as an energy source.

The predominant cellular fatty acids of strain NW10^T were C_16:1 ω7c (31.6%), C_16:0 (28.2%), C_18:1 ω7c (18.4%) and 3-OH C_14:0 (6.4%), which are similar to those of S. paralvinellae GO25^T, S. lithotrophicum GYSZ_1^T and S. gotlandica GD1^T, and distinctly different from those of other Sulfurimonas species. The significant difference between strain NW10^T and its closest relative S. paralvinellae GO25^T was that its prominent fatty acid in strain NW10^T was C_16:1 ω7c, representing 31.6%, whereas in S. paralvinellae GO25^T, C_18:1 ω7c was prominent, accounting for 37.0%. In addition, C_16:1 ω5c + t (1.3%) was observed in strain NW10^T, but not in S. paralvinellae GO25^T. The fatty acid profile of strain NW10^T and the related type strains are detailed in Table 1.

When strain NW10^T was incubated with hydrogen and thiosulfate as mixed electron donors and oxygen as the sole electron acceptor, elemental sulfur occurred in the culture at the mid-exponential growth phase, and accumulated in the late exponential phase and during the stationary phase (Figure 3A). When grown in MMJHS medium with neutral pH such as 7.0 (unbuffered) or weakly acidic pH such as 5.5 (in buffered medium), strain NW10^T produced large amount of naked-eye elemental sulfur in the culture. However, there was no accumulation of elemental sulfur observed in alkaline conditions, such as at pH 8.0. Elemental sulfur production was also affected by the oxygen concentration in the gaseous phase and the result showed the maximum accumulation of elemental sulfur occurred at 6% oxygen concentration. Morphological observations under SEM showed that extracellular sulfur was mainly shaped as crystalline bars (Figures 3B,C), with length varying from 10 to 100 μm. EDS analysis revealed that these crystals mainly contained elemental sulfur (Figure 3D). Raman spectroscopy further showed the presence of α-S₈ when compared against a standard chemical product of S₈⁰ (Figure 3E). This confirmed that the extracellular aggregates produced by strain NW10^T were biogenic sulfur and composed of crystalline α-S₈. These findings suggested that strain NW10^T performs incomplete oxidation of thiosulfate with elemental sulfur as an oxidation intermediate accumulated outside of the cells.

Strain NW10^T has a single circular chromosome (Supplementary Figure S4) with complete genome size of 2,432,011 bp with GC content of 37.3%, which is similar to that of S. paralvinellae GO25^T (Table 1). No plasmid was detected in the genome of this bacterium. Total 2,472 genes were predicted, which contained 2,367 protein coding genes and 57 RNA genes. The RNA genes cover 45 tRNAs and 12 rRNAs. ANI and DDH were calculated to identify the genomic similarities of strain NW10^T with other species of the genus Sulfurimonas. Pairwise ANI values between strain NW10^T and its closest relatives, S. paralvinellae GO25^T and S. autotrophica OK10^T, were 74.50 and 81.15%, respectively. The predicted DDH value between strain NW10^T and S. paralvinellae GO25^T was 19.20% and the value between strain NW10^T and S. autotrophica OK10^T was 24.70%. These findings supported the classification of strain NW10^T as one distinct species according to the cut-off thresholds of ANI (95–96%) and DDH (70%) for delineation of prokaryotic species (Chun et al., 2018).

Pan-genomic studies were carried out to investigate the genotypic features of 11 Sulfurimonas species, including five vent bacteria-S. hydrogeniphila NW10^T, Sulfurimonas sp. NW8N, Sulfurimonas sp. S2-6, S. paralvinellae GO25^T and S. autotrophica OK10^T, and six non-vent marine bacteria-Sulfurimonas sp. B2, S. xiamenensis 1-1N^T, S. lithotrophica GYSZ_1^T, S. denitrificans DSM1251^T, S. hongkongensis AST-10 and S. gotlandica GD1^T (Table 2). Comparative analyses based on orthologous groups of proteins revealed 886 core genes shared by all the eleven Sulfurimonas genomes (Figure 4A). The percentages of core genes in each genome ranged from 32.7 to 47.6%, indicating that the 11 strains of Sulfurimonas shared a low percentage of common functional proteins. The percentages of unique genes in each genome varied from 10.8 to 23.0%, however, a relative high proportion of distributed genes in each genome was detected, ranging from 36.2 to 51.0% (Figure 4B). The distributed genes usually offer species diversity, environmental adaptation and other characteristics for bacteria (Innamorati et al., 2020). In addition, considering that the definition of unique genes is subject to change, a singleton today may be reclassified into a multi-gene cluster when new genome data is included in the future (Zhang and Sievert, 2014), we decided to focus on the distributed genes to investigate common adaptation characteristics of the Sulfurimonas genus to deep-sea hydrothermal vent environments.

The niche-specific genes were retrieved based on comparative genome analysis. We compared the genes common in five hydrothermal vent strains with those shared by non-hydrothermal vent strains. Overall, 57 unique genes specific of the hydrothermal vent isolates were found (Supplementary Table S2). The hydrothermal vent-specific genes encoded four major functions, including signal transduction, energy production and conversion, inorganic ion transport and cell wall/membrane/envelope biogenesis (Supplementary Table S2).

As for signal transduction, some transcriptional regulators belonging to the Per-ARNT-Sim (PAS) family were shared by the hydrothermal vent strains (Supplementary Table S2). Two histidine kinases (clusters 1,373 and 1,406), as the sensing component of the two-component signal transduction system, and one response regulator (cluster 1,293) were also found across the five vent strains (Supplementary Table S2), which may be responsible for sensing certain hydrothermal vent environmental conditions. In addition, the genomes of vent strains encoded relatively high numbers of signaling proteins (Supplementary Table S2), and particularly the genes encoding proteins with EAL and GGDEF domains, likely involved in the synthesis and hydrolysis of the intracellular signaling compound cyclic diguanylate. Regarding energy conservation, some genes involved in energy metabolism such as pyruvate dehydrogenase complex, which converts pyruvate to acetyl-CoA, NADH and CO₂, were also shared by vent strains (Supplementary Table S2). We also found three rhodanese-related sulfur transferase (clusters 1,452, 1,501, and 1,507), involved in sulfur metabolism, uniquely present in the genomes of all vent strains. In addition, the vent strains shared multiple transporters, including the ABC transporter systems (clusters 1,349, 1,362, and 1,399) and some metal ion transporters such as Zn²⁺, Mg²⁺, Cu²⁺ and K⁺ (Supplementary Table S2). Finally, we found three outer membrane proteins TolC (clusters 1,471, 1,475, and 1,511), previously associated with different efflux systems and type I protein secretion (Pérez-Llarena and Bou, 2016), shared by all vent strains.

Strain NW10^T shared the highest 16S rRNA gene sequence similarity (95.8%) with S. paralvinellae GO25^T and formed a phylogenetic subcluster within the genus Sulfurimonas, indicating that this strain should belong to the Sulfurimonas genus. Furthermore, ANI and DDH analyses clearly showed that strain NW10^T could be differentiated genetically from other previously described Sulfurimonas species. Phenotypically, strain NW10^T was different from its closest relative S. paralvinellae GO25^T in many characteristics (Table 1), such as the growth conditions and maximum growth rate. The utilization patterns of electron donor were also different, in that strain NW10^T could oxidize sulfide and elemental sulfur, whereas S. paralvinellae GO25^T could not (Table 1). The combined evidence in phenotypic, chemotaxonomic and phylogenetic features conclusively demonstrates that strain NW10^T represents a novel species in the genus Sulfurimonas, for which the name Sulfurimonas hydrogeniphila sp. nov. is proposed.

In deep-sea hydrothermal environments, bacteria of this species are widely spread, as supported by ecological distribution search in different environments through the Integrated Microbial Next Generation Sequencing (IMNGS) (Lagkouvardos et al., 2016). Across the 274,621 deep-sea hydrothermal vent samples gathered in the IMNGS platform, the relative abundance of Sulfurimonas sp. NW10-like sequences (>97% similarity of 16S rRNA gene in length 1,447 bps) was more than 1% total sequence reads in 49% of the vent samples and more than 0.1–1% in 39% of the vent samples (Wang et al., 2020a). The highest relative abundance was found in a deep-sea hydrothermal vent fluid sample from the Pacific Ocean, accounting for 22.4% of the total 16S sequences (Wang et al., 2020a). The relative high abundance and wide distribution of the new species highlight its importance in the biogeochemical cycles of sulfur in in situ hydrothermal environments.

In sulfur-oxidizing bacteria, thiosulfate is usually oxidized by a Sox multi-enzyme system (SoxABCDXYZ) located in the periplasm. Two kinds of Sox pathways have been described according to the presence or absence of SoxCD. When SoxCD is present, it acts as a sulfur dehydrogenase and completely oxidizes thiosulfate to sulfate without formation of sulfur globule. Otherwise, sulfur is formed as an intermediate without SoxCD (Frigaard and Dahl, 2008). In this study, we found the soxCD genes were present in all Sulfurimonas genomes (Table 3). In addition, previous study showed that Sulfurimonas species did perform a complete thiosulfate oxidation without elemental sulfur as an intermediate (Inagaki et al., 2003; Sievert et al., 2008a; Labrenz et al., 2013). However, strain NW10^T seems to be an exception in genus Sulfurimonas. Despite the presence of SoxCD, it performed incomplete oxidation of thiosulfate, with elemental sulfur accumulation outside of the cells (Figure 3).

The accumulation of extracellular S⁰ was significantly influenced by culture conditions, such as pH. At alkaline pH, strain NW10^T completely oxidized thiosulfate to sulfate, while under neutral and acidic conditions, thiosulfate was oxidized to sulfate, with formation of elemental sulfur as an intermediate. The phenomenon was also observed in different strains of Hydrogenovibrio genus (Javor et al., 1990; Houghton et al., 2016; Jiang et al., 2017). Considering that the vent fluids of black chimneys are typically acidic, it is likely that the Sulfurimonas species inhabiting on vent chimneys can generate extracellular S⁰ in situ. In addition, oxygen concentrations can also significantly influence the elemental sulfur generation, with a maximum accumulation of extracellular S⁰ occurred at 6% oxygen.

Furthermore, the structures of extracellular S⁰ produced by NW10^T were mainly in the form of crystalline bars with clear edges composed of α-S₈, which was significantly different from the sulfur globules usually resulting from most chemotrophic and phototrophic bacteria activity (Dahl and Prange, 2006; Jiang et al., 2017; Cron et al., 2019). α-S₈ is the thermodynamically most stable form of elemental sulfur at ambient temperature and pressure, and has been found in very diverse environments, such as marine sediments, water columns, euxinic lakes, sulfidic caves, hydrothermal vents, as well as cold or hot springs (Roy and Trudinger, 1970; Taylor et al., 1999; Steudel, 2000; Gleeson et al., 2012; Findlay et al., 2014; Hamilton et al., 2015). This crystalline sulfur could serve as an important intermediate in the biogeochemical sulfur cycle, and be further consumed by a wide diversity of microorganisms, such as sulfur oxidizer, sulfur reducer, or sulfur disproportionator. Our results highlight the potential ecological significance of extracellular S⁰ produced by Sulfurimonas in deep-sea hydrothermal ecosystems.

The mechanisms involved in the formation of extracellular S⁰ by SOB are still enigmatic compared with those of intracellular S⁰, although a variety of chemotrophic and phototrophic SOB and archaea can produce extracellular S⁰ (Dahl and Friedrich, 2008). So far, a large number of studies focused on the analysis of chemical form and structure of microbial extra-and intra-cellular S⁰ (Pickering et al., 2001; George et al., 2008; Berg et al., 2014). To our knowledge, the enzymes catalyzing the oxidation of sulfide or thiosulfate to S⁰ are usually located in the bacterial periplasm, but sulfur globules are observed outside of the cells, and sometimes even keeping the cells at a distance (Gregersen et al., 2011; Marnocha et al., 2016; Cron et al., 2019). Therefore, it has been proposed that reduced sulfur compounds could be initially oxidized to soluble polysulfide intermediates in the periplasm, and then be transported outside the cells to form sulfur globules. It is still not clear how extracellular S⁰ can accumulate outside of the cells (Cron et al., 2019). Recently, increasing number of studies indicate that SOB could excrete soluble organics to help form and stabilize S⁰ in the environment (Cosmidis et al., 2019; Cron et al., 2019; Marnocha et al., 2019). This could explain many of the puzzling fact observed in microbial extracellular S⁰, such as coated by organics on the surface, and growing extracellularly at a distance from the cells (Hanson et al., 2016; Marnocha et al., 2016, 2019). Here, the extracellular S⁰ produced by strain NW10^T mainly contained elemental sulfur, with extremely low amounts of carbon and oxygen. In addition, there were no known homologs of sulfur globule proteins (SGP) genes in its genome, and comparative genome analysis showed no significant difference in Sox pathways between strain NW10^T and other Sulfurimonas species and the only difference is that SoxC protein of strain NW10^T lacks 18 bases at the N-terminal. The mechanism of extracellular S⁰ production requires further investigations by means of multiple omics in future.

Bacteria of the Sulfurimonas genus have been found to colonize a broad range of natural habitats from terrestrial, coastal sediment, shallow waters to deep-sea hydrothermal vents globally (Grote et al., 2008; Dahle et al., 2013; Meier et al., 2017). In this report to understand their environmental adaptation, we carried out comparative genomic analyses between vent and non-vent marine Sulfurimonas strains. Regarding energy conservation, all Sulfurimonas genomes contained the broad suite of genes encoding the enzymes capable of oxidizing thiosulfate, sulfite, sulfide and hydrogen (Table 3). In relation to sulfide oxidation, Sulfurimonas species contained genes encoding diverse types of Sqrs, allocated into Types II, III, IV, V, and VI. The variation and diversification of these Sqrs are presumed to play important roles in sulfide oxidation, sulfide assimilation, energy generation, heavy metal tolerance, detoxification and sulfide signaling (Marcia et al., 2010). Type IV and VI Sqrs were relatively conserved in genus Sulfurimonas (Table 3), likely to cope with different concentrations of sulfide, such as free sulfide in hydrothermal vents and metal sulfide in marine sediments (Han and Perner, 2015). Additionally, Sulfurimonas bacteria had one or more copies of Type II Sqrs (Table 3), which may compensate for the function of other Sqrs under specific environmental conditions (Han and Perner, 2015). Distinct roles of Type II Sqr have been proposed in different microorganisms. For example, it may be involved in heavy metal tolerance in the yeast Saccharomyces pombe (Vande-Weghe and Ow, 1999), sulfur assimilation in the non-pathogenic bacterium Pseudomonas putida KT2440 (Shibata and Kobayashi, 2006), and sulfide signaling in mammalian cells (Shahak and Hauska, 2008). In addition, some Sulfurimonas species harbored Type III and V Sqrs besides Type II, IV and VI Sqrs (Table 3), which may function in sulfide oxidation to enhance energy generation or detoxification and sulfide signaling (Han and Perner, 2015). It is interesting to find that Type V Sqr is only present in vent strains, such as NW10^T, S2-6 and S. autotrophica OK10^T. Type V Sqrs are a group of archaeal proteins, predominantly found in Creanarcheaota and Euryarchaeota (Sousa et al., 2018). It was only characterized in the archaeon Acidianus ambivalens and has been confirmed to possess Sqr activity with high substrate affinity (Brito et al., 2009). In addition, this type of Sqr in A. ambivalens exhibited the highest activity at 70°C, and only 3% activity remained at 25°C (Brito et al., 2009). This indicates that Type V Sqr may tolerate elevated temperatures. Phylogenetic relationship showed that Sulfurimonas has likely acquired its Type V Sqr gene from archaea through horizontal gene transfer (Figure 5; Han and Perner, 2015). Further experiments will be required to confirm the temperature tolerance attributed to Type V Sqr. In addition, genes encoding rhodanese-related sulfurtransferase were significantly enriched in vent strains (Supplementary Table S1). Previous study has indicated that the rhodanese-related sulfurtransferase serves as a polysulfide-sulfur transferase at lower polysulfide concentration in Wolinella succinogenes (Klimmek et al., 1991) and was possibly involved in polysulfide reduction in Desulfurella amilsii (Florentino et al., 2019). Here, we hypothesize that these sulfur transferases from vent strains may play a key role in the sulfur/polysulfide respiration process, where elemental sulfur or polysulfides are abundant.

Hydrogen oxidation is popular in bacteria of the Sulfurimonas genus. Members of this genus harbor a relatively broad suite of hydrogenases (Table 3), which can catalyze the hydrogen oxidation coupled to various electron acceptors reduction, including oxygen, nitrate and elemental sulfur. Almost all Sulfurimonas members contained multiple copies of Group I hydrogenases for hydrogen uptake, likely exhibiting different hydrogen affinities (Grote et al., 2012; Labrenz et al., 2013). Group II hydrogenases likely participate in hydrogen sensing or energy conversion at a low concentration of hydrogen (Campbell et al., 2009; Sievert et al., 2008b; Han and Perner, 2015). Previous study showed that Group II hydrogenases were not found in the hydrothermal vent Sulfurimonas isolates, but present in the marine water and sediments, thus speculated that in these habitats hydrogen concentration is relatively low than in hydrothermal vents, and that this group of hydrogenases may be more important under low hydrogen concentrations (Han and Perner, 2015). Yet, in this study, Group II hydrogenase was found in vent Sulfurimonas species including strains NW10^T, NW8N, and S2-6. Hence, it is unlikely that Group II hydrogenase in vent Sulfurimonas is specialized to function at low hydrogen concentrations, and the elucidation of the role of Group II hydrogenase in Sulfurimonas species will require examination of its activity under different hydrogen concentrations.

Additionally, four Sulfurimonas species harbored Group IV hydrogenases, possibly involved in hydrogen evolution or energy conversion (Vignais and Billoud, 2007). In addition to Hyc and Ech, Coo subtype of group IV hydrogenase was first found in vent strain Sulfurimonas sp. S2-6, containing the cluster CooLUHF (Table 3). This energy-converting hydrogenases can couple CO and H₂ metabolism with energy conservation. The Coo hydrogenase cluster has been identified in 30 bacterial representatives, at particularly high frequency in Deltaproteobacteria (many sulfate reducers, e.g., Desulfovibrio vulgaris), Alpha- (e.g., Rhodospirillum rubrum) and Campylobacteria (Nautilia profundicola), Firmicutes (e.g., Carboxydothermus hydrogenoformans), Betaproteobacteria and Gammaproteobacteria (Schoelmerich and Müller, 2019). Phylogenetic analysis based on the large subunit CooH indicated that this hydrogenase was clustered together with those of R. rubrum, D. vulgaris and C. hydrogenoformans (Figure 6), which have been experimentally demonstrated to oxidize CO, producing carbon dioxide and hydrogen as products, through CooMKLXUHF complexes (Bonam and Ludden, 1987; Wu et al., 2005; Caffrey et al., 2007). Further, protein domain analyses by DELTA-BLAST and PSI-BLAST in NCBI showed that CooH contains three functional domains: two respiratory-chain NADH dehydrogenases and one nickel-dependent hydrogenase. The Coo hydrogenases could be involved in H₂ production from CO according to the overall equation CO + H₂O→CO₂ + H₂ (Heidelberg et al., 2004). In any case, this is the first report of Coo hydrogenase found in the genus Sulfurimonas and the function of H₂ production from Coo hydrogenase in Sulfurimonas species needs further experimental confirmation.

In a word, the evolution and functions of different types of Sqr and hydrogenases within one host remain enigmatic and need further investigations, especially to elucidate their relevance in relation to host adaptation to hydrothermal environments. In addition to energy metabolism, comparative genome analysis revealed other vent-specific gene signatures related to signal transduction mechanisms and inorganic ion transporter mechanisms, including unique two-component signal transduction system and a relative abundance of signaling proteins, the ABC transporter system and metal ion transporter (such as Zn²⁺, Mg²⁺, Cu²⁺ and K⁺ transporter), and outer membrane protein TolC. Overall, this versatile energy metabolism, environmental sensing systems, and multiple transporter mechanisms could contribute to the wide spreading and high adaptability of these organisms to different hydrothermal vent fields globally.

Sulfurimonas hydrogeniphila (hy.dro.ge.ni.phi’la. N.L. neut. n. hydrogenum hydrogen; Gr. adj. philos loving; N.L. fem. adj. hydrogeniphila hydrogen-loving, because growth prefers hydrogen).

Cells are Gram-negative short rods (0.8–3.5 × 0.4–0.8 μm) that are motile by means of a polar flagellum. Anaerobic to aerobic, even with the air in the headspace (optimum 6%). Growth occurs at 4-45°C (optimum 33°C), pH 5.0–9.0 (optimum pH 6.0–6.5), and with 2.0–4.0% (w/v) NaCl concentration (optimum 3.0%). Obligate chemolithoautotrophic growth occurs with H₂, thiosulfate, sulfide, elemental sulfur as electron donors, and oxygen, nitrate, and S⁰ can be utilized as an electron acceptor. Organic substrates are not utilized as carbon sources and energy sources. Major cellular fatty acids are C_16:1 ω7c (31.6%), C_16:0 (28.2%) and C_18:1 ω7c (18.4%).

The type strain, NW10^T (=MCCC 1A13987^T = KTCC 15781^T) was isolated from the deep-sea hydrothermal sulfide chimneys in the Carlsberg Ridge of Northwest Indian Ocean. The G + C content of its genomic DNA is 37.3 mol%.

The chimney samples were collected near an active hydrothermal vent on the Carlsberg Ridge (60°31′E, 6°21′N), at a depth of 2936 m, by a human operated vehicle “Jiao-long” during COMRA DY 38 oceanic scientific cruise in March 2017. Aboard the research vessel Xiang-Yang-Hong No. 9, chimney samples were immediately transferred into MMJHS medium under a gas phase mixture of 80% H₂/18% CO₂/2% O₂ (200 kPa) and then incubated at 28°C according to the previous description (Inagaki et al., 2003). After successful enrichment with MMJHS medium, the well-grown culture was further purified using the dilution-to-extinction technique with the same medium. MMJS medium consisted of NaCl (30 g l^–1), KCl (0.33 g l^–1), NH₄Cl (0.25 g l^–1), MgCl₂⋅6H₂O (4.18 g l^–1), CaCl₂⋅2H₂O (0.14 g l^–1), K₂HPO₄ (0.14 g l^–1), NaHCO₃ (1 g l^–1), Na₂S₂O₃⋅5H₂O (10 mM), Wolfe’s vitamins (1 ml l^–1) and trace element solution (10 ml l^–1).

The genomic DNA was prepared according to the method described previously (Jiang et al., 2010) and the 16S rRNA gene was amplified by PCR primers described previously (Lane, 1991). The sequence was compared with those of other type strains using the EzTaxon-e server (Yoon et al., 2017). The 16S rRNA gene sequences of the related taxa were obtained from the GenBank database. Phylogenetic trees were constructed by using MEGA 6.0 (Tamura et al., 2013) using the neighbor-joining (Saitou and Nei, 1987), maximum-likelihood (Felsenstein, 1981) and minimum evolution methods (Rzhetsky and Nei, 1992) after multiple alignments of the data by CLUSTAL_W. Evolutionary distances were calculated using Kimura’s two-parameter model and bootstrap values were determined based on 1,000 replications. The phylogenomic tree was constructed based on an up-to-date 92 bacterial core gene sets by UBCG version 3.0 (Na et al., 2018). Genome sequences of reference taxa were retrieved from the NCBI database and the 92 concatenated core genes were extracted, aligned and concatenated using default parameters. The tree topology was supported by the maximum-likelihood method for 100 bootstrap replications using RAxML version 8.2.11 (Stamatakis, 2014) with the GTR + CAT model.

Cell morphology was observed under a transmission electron microscopy (Model JEM-1230, JEOL, Japan) with cultures grown in MMJHS medium at 28°C for 1 day. The physiological characterization of the isolate was tested on MMJHS medium (Inagaki et al., 2003). After autoclaving, the medium (10 ml) was dispensed into 50 ml serum bottles, then sealed with a butyl-rubber stopper under a gas phase of 80% H₂/18% CO₂/2% O₂ (200 kPa). All cultivation experiments were performed in triplicate, unless otherwise specified. The growth was measured by direct cell counting using a phase contrast microscope (Eclipse 80i, Nikon, Japan). Growth at different temperatures was examined in MMJHS at 4, 10, 15, 20, 25, 28, 30, 33, 35, 37, 45, 50, and 60°C. The growth salinity range was examined by adjusting the concentrations of NaCl between 0 and 9.0% (w/v), at 0.5 (w/v) intervals. To determine the effect of pH on growth, the pH of MMJHS medium was adjusted from 4.5 to 9.0 with a 0.5 pH unit interval by using 30 mM acetate/acetic acid buffer (pH 4.0–5.0), MES (pH 5.0–6.0), PIPES (pH 6.0–7.0), HEPES (pH 7.0–7.5), Tris and CAPSO (pH 8.0 and above). The effect of O₂ on growth was examined by adjusting the oxygen concentration (0, 1, 2, 4, 6, 8, 10% at 200 kPa and 20% at 100 kPa) in the headspace gas. To test the anaerobic growth, 10 mM nitrate was added as an alternative electron acceptor.

The ability of sulfur oxidation of strain NW10^T was tested using the following sulfur compounds as the sole energy source in the MMJS medium, including thiosulfate (10 mM), sulfite (5 mM), thiocyanate (5 mM), tetrathionate (5 mM), elemental sulfur (1% w/v), or sodium sulfide (50, 100, 500 μM, and 1 mM) with a gas phase of 74% N₂/20% CO₂/6% O₂ (200 kPa). Molecular hydrogen was also tested in MMJH medium in the absence of thiosulfate under a gas phase of 74% H₂/20% CO₂/6% O₂ (200 kPa). To determine the utilization of other electron acceptors, each of the potential electron acceptors, such as thiosulfate (10 mM), tetrathionate (10 mM), sulfite (2 mM and 10 mM), elemental sulfur (1%, w/v), nitrate (10 mM), nitrite (1 mM and 5 mM), selenate (5 mM), arsenate (5 mM), fumarate (10 mM), and ferric citrate (20 mM) was examined with MMJHS medium under 80% H₂/20% CO₂ (200 kPa). Heterotrophic growth was examined in MMJHS medium in the absence of NaHCO₃ under a gas phase of 94% N₂/6% O₂ (200 kPa). Each of the following potential organic carbon compounds was tested: 0.1% (w/v) peptone, starch, tryptone, yeast extract, casein, and casamino acids, 5 mM formate, acetate, propionate, citrate, tartrate, fumarate, succinate, malate, and pyruvate, 5 mM of each of 20 amino acids, 0.02% (w/v) glucose, galactose, sucrose, fructose, lactose, maltose, and trehalose. In an attempt to determine the alternative energy source, these organic compounds were used as an energy source in MMJ medium to replace thiosulfate under a gas phase of 74% N₂/20% CO₂/6% O₂ (200 kPa).

For analyses of fatty acids, cells grown on MMJHS medium at 33°C for 24 h were saponified, methylated, and extracted following the standard MIDI protocol (Sherlock Microbial Identification System, version 6.0B). The fatty acids were analyzed by gas chromatography (Agilent Technologies 6850) and then the result was identified using the TSBA6.0 database of the Microbial Identification System.

Scanning electron microscope (SEM) (S-3400N; Hitachi, Japan) and Raman spectra (XploRA; Horiba JY, France) were used to identify the shape, components and structure of extracellular sulfur produced by strain NW10^T. For SEM analysis, a milky white suspension was collected using polycarbonate filters (Merck Millipore, pore size 3.0 μm), rinsed three times with deionized water and observed by SEM at 5 kV. Energy-Dispersive Spectrum (EDS) (model 550i, IXRF systems, United States) equipment with SEM was employed at an accelerating voltage of 5 keV for 30 s. For Raman analysis, about 5 ml samples were concentrated by centrifugation and rinsed three times with deionized water. The supernatant was decanted and the pellet was lyophilized overnight. A small amount of powder was put on the glass slide and observed by microconfocal method. Raman spectra were collected using a Horiba XploRA Raman spectrometer coupled with a 532 nm laser source. For α-S₈ standard, commercial precipitated sulfur (purity > 99%, Thermo Fisher Scientific) was used.

The complete genome of strain NW10^T was sequenced by Tianjin Biochip Corporation (Tianjin, PR China), using the single molecule real-time (SMRT) technology on the Pacific Biosciences (PacBio) sequencing platform. The sequenced reads were filtered, and high quality paired-end reads were assembled to construct a circular genome with SOAPdenovo (version 2.04)¹. The G + C content of the chromosomal DNA was determined according to the genome sequence. Gene prediction was carried out by Glimmer program (Delcher et al., 2007). rRNA identification was performed with the RNAmmer 1.2 software (Lagesen et al., 2007), and tRNAscan-SE (version 1.21) was used to identify the tRNA genes (Schattner et al., 2005). Gene prediction and annotation were carried out using NCBI Prokaryotic Genomes Annotation Pipeline (PGAP) and the Rapid Annotation using Subsystem Technology (RAST) pipeline² (Aziz et al., 2008). The functional annotation and metabolic pathways were analyzed by searching against KEGG and COG databases. To further clarify the genetic relatedness between strain NW10^T and related species of the genus Sulfurimonas, the average nucleotide identity (ANI) value between two genomes was calculated using the web service of EZGenome³ (Richter and Rosselló-Móra, 2009). The predicted in silico DNA-DNA hybridization (DDH) values were determined online⁴ using the Genome-to-Genome Distance Calculator (GGDC) (Auch et al., 2010).

To avoid the possible deviations due to different annotation methods, we used RAST server for re-annotation. A pan-genome for the eleven genomes was determined by BPGA (Chaudhari et al., 2016) pipeline to identify orthologous groups among Sulfurimonas strains and to extrapolate the pan-genome models of applying default parameters. Orthologous clusters were assigned by grouping all protein sequences in the 11 genomes using USEARCH based on their sequence similarity (E-value < 10^–5, >50% coverage) and each protein was assigned to one protein family. The pan genome analysis complied the set of core genes shared among all strains, a set of distributed genes shared with more than two but not all strains, and unique genes only found in a single strain. COG and KEGG distributions of the core, distributed and unique gene families were calculated based on representative sequences.