
==== Front
mSystems
mSystems
msystems
mSystems
2379-5077
American Society for Microbiology 1752 N St., N.W., Washington, DC

39166872
msystems00148-24
10.1128/msystems.00148-24
msystems.00148-24
Research Article
environmental-microbiologyEnvironmental MicrobiologyNovel isolates of hydrogen-oxidizing chemolithoautotrophic Sulfurospirillum provide insight to the functions and adaptation mechanisms of Campylobacteria in shallow-water hydrothermal vents
https://orcid.org/0000-0002-8340-9615
Wang Li 1 Conceptualization Data curation Formal analysis Funding acquisition Methodology Writing – original draft l-wang@shou.edu.cn

Cheng Xinyi 1 Methodology Project administration
Guo Yizhe 1 Data curation Software
Cao Junwei 1 Conceptualization Writing – review and editing
Sun Mingye 1 Project administration
Hwang Jiang-Shiou 2 Resources
https://orcid.org/0000-0001-5198-0422
Liu Rulong 1 Conceptualization Supervision Writing – review and editing rlliu@shou.edu.cn

Fang Jiasong 1 3 Funding acquisition Writing – review and editing
1 College of Oceanography and Ecological Science, Shanghai Ocean University , Shanghai, China
2 Institute of Marine Biology, National Taiwan Ocean University , Keelung, Taiwan
3 Laboratory for Marine Biology and Biotechnology, Qingdao National Laboratory for Marine Science and Technology , Qingdao, China
Editor Hernandez Marcela University of East Anglia , Norwich, United Kingdom

Address correspondence to Li Wang, l-wang@shou.edu.cn
Address correspondence to Rulong Liu, rlliu@shou.edu.cn
The authors declare no conflict of interest.

9 2024
21 8 2024
21 8 2024
9 9 e00148-2408 2 2024
24 7 2024
Copyright © 2024 Wang et al.
2024
Wang et al.
https://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

ABSTRACT

Enhancing the availability of representative isolates from hydrothermal vents (HTVs) is imperative for comprehending the microbial processes that propel the vent ecosystem. In recent years, Campylobacteria have emerged as the predominant and ubiquitous taxon across both shallow and deep-sea vent systems. Nevertheless, only a few isolates have been cultured, primarily originating from deep-sea HTVs. Presently, no cultivable isolates of Campylobacteria are accessible in shallow water vent systems (<200 m), which exhibit markedly distinct environmental conditions from their deep-sea counterparts. In this study, we enriched a novel isolate (genus Sulfurospirillum, Campylobacteria) from shallow-water HTVs of Kueishan Island. Genomic and physiological analysis revealed that this novel Campylobacteria species grows on a variety of substrate and carbon/energy sources. The pan-genome and phenotypic comparisons with 12 previously isolated Sulfurospirillum species from different environments supported the identification of functional features in Sulfurospirillum genomes crucial for adaptation to vent environments, such as sulfur oxidation, carbon fixation, biofilm formation, and benzoate/toluene degradation, as well as diverse genes related with signal transportation. To conclude, the metabolic characteristics of this novel Campylobacteria augment our understanding of Campylobacteria spanning from deep-sea to shallow-water vent systems.

IMPORTANCE

Campylobacteria emerge as the dominant and ubiquitous taxa within vent systems, playing important roles in the vent ecosystems. However, isolated representatives of Campylobacteria have been mainly from the deep-sea hydrothermal fields, leaving a significant knowledge gap regarding the functions, activities, and adaptation strategies of the vent microorganisms in shallow-water hydrothermal vents (HTVs). This study bridges this gap by providing insights into the phenomics and genomic diversity of genus Sulfurospirillum (order Campylobacterales, class Campylobacteria) based on data derived from a novel isolate obtained from shallow-water HTVs. Our mesophilic isolate of Sulfurospirillum not only augments the genus diversity of Campylobacteria pure cultures derived from vent systems but also serves as the inaugural reference isolate for Campylobacteria in shallow-water environments.

KEYWORDS

Sulfurospirillum
hydrogen-oxidizing
Campylobacteria
shallow-water hydrothermal vents
Chemoautolithotrophic
MOST | National Natural Science Foundation of China (NSFC) 41906134 Wang Li MOST | National Natural Science Foundation of China (NSFC) 42276149 Liu Rulong MOST | National Natural Science Foundation of China (NSFC) 92251303 Fang Jiasong cover-dateSeptember 2024
==== Body
pmcINTRODUCTION

Hydrothermal vents (HTVs), where geothermally heated water is expelled through fissures in Earth’s crust, constitute one of the most productive marine ecosystems. As hydrothermal fluid rises from the deep crust, a variety of chemicals and volatile gases enriched in the fluid escape to the surrounding seawater. This mixing of the expelled fluid and surrounding seawater occurs both below and above the seafloor, spanning vast distances. Chemolithoautotrophs within the vent system, including animal symbionts, microbial mats members, and free-living cells, utilize various reductive substrates (such as H2, H2S, CH4, and Fe2+) as energy sources to fix carbon (1, 2). Subsequently, this organic carbon sustains dense animal communities through symbiotic relationships with bacteria, as well as through grazing, suspension feeding, and subsequent trophic transfer. While vent ecosystems in the deep-sea (>200 m) primarily rely on microbial chemosynthesis (3), shallow-water hydrothermal vents (SW-HTVs) at depths <200 m exhibit distinct environmental conditions characterized by the presence of light, terrestrial inputs, tidal cycles, and meteoric water inputs. These factors contribute to the formation of highly diverse and complex microbial communities (4, 5).

Campylobacteria (formerly known as Epsilonproteobacteria) (6) are commonly detected chemoautotrophic taxa in vent systems. Despite the discovery of deep-sea HTVs in 1977, Campylobacteria have been identified as the dominant bacteria through molecular studies since the late 20th century (7–9). They have been found in various biomes of vent systems, including chimneys, subsurface, diffuse flow, fluids, or even as symbiont of metazoans (1, 2). Campylobacteria can comprise over 80% of microbial communities in the chimney structures (10, 11), diffuse flow vent emissions (12), and over 90% in subsurface sediment (13) and plume (14). Cultivation breakthrough have revealed their metabolic versatility, with sulfide and hydrogen gas as electron donors, and elemental sulfur and nitrate as electron acceptors for lithotrophic growth (15, 16). Phylogenetically diverse Campylobacteria representatives have been isolated from deep-sea HTVs, including genera Sulfurovum (17–21), Sulfurimonas (22–25), Caminibacter (26), Cetia (27), Lebetimonas (28), Nautilia (29–31), Nitratiruptor (32), Thioreductor (33), and Hydrogenimonas (34, 35). While Sulfurovum and Sulfurimonas belong to the order Campylobacterales, others are from the order Nautiliales, comprising moderately thermophilic bacteria. Unlike Nautiliales, most Campylobacterales representatives thrive at lower temperatures and tolerate a wider range of oxygen concentrations (15). In SW-HTVs, Campylobacteria have been identified as dominant bacteria in fluids and sediments with moderate to high temperatures (35–55°C) and high concentrations of hydrogen sulfide and other reduced sulfur compounds, but no Campylobacteria have been isolated from SW-HTVs (5). Considering the significant environmental differences between SW-HTVs and deep-sea HTVs (4, 5), bacteria in SW-HTVs may exhibit unique functions and adaptations compared with their deep-sea counterparts. Successful cultivation of Campylobacteria isolates from SW-HTVs may shed light on their ecological roles and adaptation mechanisms to shallow water vent systems.

Kueishan Island is a young volcanic island situated at the southernmost part of the Okinawa Trough of the western Pacific Ocean, housing a series of SW-HTVs off its coast. The hydrothermal fluid and volcanic gases emanating from the Kueishan Island SW-HTVs contain high concentrations of CO2 and H2S (36). Various sulfur formations, such as sulfur chimneys, sulfur sands, and sulfur balls, are widely distributed on the seabed (37). Previous studies have consistently identified Campylobacteria as the dominant taxa in the total bacterial community based on DNA or RNA analyses (38–43). In our study, we successfully isolated multiple strains of Campylobacteria from the Kueishan Island SW-HTVs using optimized isolation medium. These isolates were classified as novel species within the genus Sulfurospirillum. Focusing on one of the strains, we assessed its chemoautotrophic activities using H2 or simple organics like formate as the electron donor, and sulfur or nitrate as the electron acceptor. Through comprehensive physiological, genomic, and comparative genomic analysis, we uncovered unique metabolic characteristics and potential adaptation mechanisms of a novel isolate to SW-HTV environments.

MATERIALS AND METHODS

Enrichment and isolation

Sulfur-rich sediments of the Kueishan Island SW-HTV (121.96232°E, 24.83420°N), at depth of 21 m, were used as the inoculum in our study. Detailed sampling information and environmental parameters can be found in our previous studies (41, 42). The surface portion of each sediment core was discarded and the remaining portion was stored in sterile polypropylene bags and kept on ice during transportation to the laboratory.

In the laboratory, 25 g sediment sample was mixed with 20 mL MJ synthetic seawater (44) in a 120 mL serum bottle (Whatman) under a gas phase of pure N2 (200 kPa) at 28°C. After 48 h enrichment, 0.5 mL of the sediment slurry was transferred to a new 20 mL MMJHS medium for longer culture under a gas phase of N2/H2/CO2 (50:40:10, 200 kPa). Microbial pure cultures were obtained using a dilution-to-extinction approach. MMJHS medium was prepared by dissolving 1 g each of NaHCO3, Na2S2O3∙5H2O and NaNO3, 3 g S0, and 10 mL trace vitamin solution (45) in 1 L of MJ synthetic seawater (44). All the ingredients except sulfur were mixed, and the pH was adjusted to 6.5 with NaOH. Then, the medium was filter-sterilized with 0.2 µm filter (Whatman, United Kingdom). Next, 20 mL medium was distributed into 120 mL serum bottle, which was autoclaved together with sulfur power. Resazurin was added as a redox indicator to monitor bacterial growth. However, we did not add resazurin to the medium during subsequent substrate and growth factor tests to avoid potential disturbances to the optical density (OD) values.

Substrate test

The spectrum of electron donors and electron acceptors utilized was determined using MJ synthetic seawater containing 0.1% (wt/vol) NaHCO3 as the basal medium. Potential electron acceptors, including S0 (0.3%, wt/vol), Na2S2O3∙5H2O (0.1%, wt/vol), Na2SO3(0.1%, wt/vol), Na2SO4 (0.1%, wt/vol), NaNO3 (0.1%, wt/vol), and NaNO2 (0.1%, wt/vol), were tested individually using hydrogen as the electron donor in a mixed gas phase (N2/H2/CO2 = 50:40:10, 200 kPa). For experiments testing the capability of the isolates to utilize formate (0.1%, wt/vol), fumarate (0.1%, wt/vol), succinate (0.1%, wt/vol), malate (0.1%, wt/vol), pyruvate (0.1%, wt/vol), or lactate (0.1%, wt/vol) as an electron donor, S0 (0.3%, wt/vol) was added into the culture medium as electron acceptor, and N2 (200 kPa) was provided as the gas phase. Additionally, the ability of the isolates to grow by fermentation on glucose, mannitol, inositol, rhamnose, sucrose, melibiose, amygdalin, and arabinose was tested with API 20E test kits (BiomCrieux, Marcy-I'Etoile, France) under oxygen-free conditions.

To determine nitrogen utilization, 0.10% (wt/vol) NH4Cl or NaNO3 was added to MMJHS medium lacking of nitrogen sources, under a N2/H2/CO2 (50:40:10, 200 kPa) gas phase. Sulfur utilization was tested by providing S0 (0.3%, wt/vol) and/or Na2S2O3∙5H2O (0.1%, wt/vol) as the sulfur source to the medium, under a N2/H2/CO2 (50:40:10, 200 kPa) gas phase.

Morphology and growth factors

The morphology of the cells was examined using fluorescence microscope and transmission electron microscopy. To assess potential growth controlling factors such as pH, temperature, and NaCl concentration, changes in OD at 600 nm were monitored. Sensitivity to oxygen was determined by adding different concentrations of an oxygen scavenger to the medium under the gas phase composed of N2/H2/CO2 = 50:40:10, at 200 kPa. The oxygen scavenger utilized was a mixture of 2.5% Na2S·9H2O and 2.5% cysteine HCl.

Analysis of fatty acid, polar lipids, and respiratory lipoquinones

Fatty acid, respiratory quinones, and polar lipids were extracted from 100 mg, 400 mg and 1 g of freeze-dried cells, respectively. Fatty acids were saponified, extracted, and methylated using the standard protocol of MIDI (Sherlock Microbial Identification System, version 6.0) and identified by using the RTSBA6.0 database of the Microbial Identification System (46). Polar lipids of freeze-dried cells were extracted and separated on silica gel 60F254 aluminum-backed thin-layer plates (10 × 10 cm2; Merk 5554), which were dried for 30 min at 55°C and further analyzed according to Minnikin et al. (47). The first dimension of the solvent system was chloroform/methanol/water (65:24:4, by volume) and the second dimension was chloroform/glacial acetic acid/methanol/water (80:15:12:4, by volume). Other reagents such as α-naphthol, ninhydrin, and molybdenum blue (Sigma) were used to detect glycolipids, amino lipids, and phospholipids, according to Tindall (48). Phosphomolybdic acid (5%, wt/vol, dissolved in alcohol) was sprayed on the plates, which were then heated at 160°C for 10–15 min, to identify total lipids. Respiratory quinones were extracted using the method described by Minnikin et al. (47) and analyzed by HPLC as described by Tindall (48).

Molecular analysis

The 16S rRNA gene was amplified via PCR using primers Eubac 27F/1492R (49) and sequenced using Sanger dideoxynucleotide chain-termination method by Sangon Biotech Co. Ltd. (Shanghai, China). Unambiguously nucleotide sequences were used for phylogenetic analysis with MEGA package (50).

Genomic DNA was extracted using MiniBEST bacteria genomic DNA extraction kit Ver.3.0 (Takara BIO INC.) according to the manufacturer’s protocol. The purified genomic DNA was quantified using a TBS-380 fluorometer (Turner BioSystems Inc., Sunnyvale, CA). High-quality DNA (OD260/280 = 1.8–2.0, >20 µg) was sequenced using a combination of PacBio RS II Single MoleculeReal Time (SMRT) and Illumina sequencing platforms by Majorbio Bio-Pharm Technology Co. Ltd. (Shanghai, China). The Illumina data were employed to evaluate the genome complexity. After quality trimming, the reads were assembled into contigs using hierarchical genome assembly process (51). Subsequently, error correction of the PacBio assembly results was performed using the Illumina reads with Pilon. Glimmer (52) was utilized for CDS prediction, tRNA-scan-SE (53) for tRNA prediction, and Barrnap for rRNA prediction. Predicted CDSs were annotated from NR, Swiss-Prot, Pfam, GO (Gene Ontology), COG (Clusters of Orthologous Groups), and KEGG (Kyoto Encyclopedia of Genes and Genomes) databases using sequence alignment tools such as BLAST, Diamond, and HMMER. Each set of query proteins was aligned with the databases, and annotations of best-matched subjects (e-value < 10−5) were obtained for gene annotation. A phylogenomic tree of Campylobacteria genomes was reconstructed based on 400 universal marker genes using PhyloPhlAn (54).

Genome relatedness was assessed using the ANI (average nucleotide identity) (55), POCP (percentage of conserved proteins), and AAI (average amino acid identity) indexes. ANI value was calculated by pyANI program (56). POCP was calculated as described before (57). AAI was generated using the amino acid identity workflow of CompareM (http://github.com/dparks1134/CompareM).

[FeFe]- and [NiFe]-hydrogenase homologs were discretely distributed at the domain, phylum, and order levels of taxonomic classifications (58). Based on the phylogeny of the large subunits, [NiFe] hydrogenase was further classified into different groups according to a previous study (59). The annotated hydrogenases of our isolates’ genome were aligned with the reference hydrogenase sequences (59) by MUSCLE (60), and a phylogenetic tree was constructed in MEGA (50).

The software was executed with default settings unless other specified.

Computing the pangenome

The anvi’o pangenome workflow (https://merenlab.org/2016/11/08/pangenomics-v2) was employed to compute the pangenome with default settings (61), involving several steps: (i) generating an anvi’o genome database of Sulfurospirillum (anvi-gen-genomes-storage) to store DNA and amino acid sequences (ii), computing the pangenome (anvi-pan-genome) from a genome database by identifying “gene clusters” (iii), displaying the pangenome (anvi-display-pan) to visualize the distribution of gene clusters across genomes and interactively bin gene clusters into logical groups such as “core genes,” “unique genes of 1612,” and “accessory genes absenting in S. arcachonens and 1612” bins (iv), reporting summary output as a profile database including the names of gene clusters and their numbers in each bin (anvi-summarize) (v), obtaining the FASTA file with sequences for each bin (anvi-get-sequences-for-gene-clusters), followed by annotation by KEGG or COG. In this context, a “gene cluster” represents sequences of one or more predicted open reading frames grouped together based on their homology at the translated DNA sequence level. Gene clusters with more than one sequence may contain orthologous or paralogous sequences, or both, from one or more genomes analyzed in the pangenome. The figure of pan-genome profile was generated by BPGA with default parameters (62).

Data description

The assembled circular genomic data of strain 1612 was deposited in NCBI with the accession number CP140614, under an original name “Sulfurospirillum sp. 1612.” For genomic comparison, the following reference genomes were utilized: Sulfurospirillum deleyianum (DSM 6946T, GCA 000024885, 2.3 Mb), Sulfurospirillum multivorans (DSM 12446T, GCA 000568815, 3.2 Mb), Sulfurospirillum arcachonense (DSM 9755T, GCA 000597725, 2.7 Mb), Sulfurospirillum barnesii (DSM 10660T, GCA 000265295, 2.5 Mb), Sulfurospirillum arsenophilum (NBRC109478, GCF 000813345, 2.6 Mb), Sulfurospirillum halorespirans (DSM 13726T, GCA 001723605, 3.0 Mb), Sulfurospirillum cavolei (UCH003, GCA 001548055, 2.7 Mb), UCH001 (GCA 001548035, 2.6 Mb), JPD-1 (GCA 002309535, 2.8 Mb), SL2-2 (GCA 002205395, 2.9 Mb), ACSTCE (GCA 011769965, 2.7 Mb), and ACSDCE (GCA 011769985, 2.8 Mb).

RESULTS

Enrichment and isolation

A range of temperatures was applied to the cultures, encompassing 28°C, 35°C, 50°C, and 70°C. After several weeks, only the cultures incubated at 28°C showed slight turbidity, indicating probable microbial growth. Consequently, subsequent isolation efforts primarily focused on incubation at 28°C. Eventually, we successfully isolated seven strains, designated as strains 1612, 2216, 34, 4316-2, 3221-6, 1307, and 2301. As depicted in the phylogenetic tree of 16S rRNA gene (Supporting materials, Fig. S1), strains 1612, 34, 2216, 3221-6, and 4326-2 were similar, clustering into one clade, while strains 1307 and 2301 formed another clade. They were genetically related to the genus Sulfurospirillum, but the sequence similarity values were not high. For example, strain 1612 exhibited low similarity in its 16S rRNA gene sequence when compared to Sulfurospirillum carboxydovorans (91.5%), S. arcachonense (93.1%), and the type strain S. deleyianum (91.9%) (Supporting materials, Table S1). Strains 1612, 2216, 34, and 1307 were subsequently purified and preserved at the Marine Culture Collection of China (Xiamen, China) with the strain preservation numbers MCCC 1K08293, MCCC 1K08414, MCCC 1K07846, and MCCC 1K08950. Purity confirmation was achieved through microscopic examination and repeated partial sequencing of the 16S rRNA gene using both Sanger sequencing and high-throughput sequencing methods. Strain 1612 was chosen as the representative for subsequent phenomics and genomic analyses due to its faster growth rate.

Morphology, physiology, and growth characteristics of the isolate

Cells of strain 1612 were observed to be slightly curved rods, with lengths ranging from 1.0 to 2.5 µm, and each cell possessed a polar flagellum (Supporting materials, Fig. S2). The cell was stained Gram-negative. Strain 1612 exhibited growth at 4–50°C, with an optimum growth temperature of 28°C. No growth was observed below 4°C or above 50°C (Supporting materials, Fig. S3A). The growing pH of the bacteria ranged from 5.5 to 8.0, with an optimum growth at pH 6.5 (Supporting materials, Fig. S3B). NaCl requirement for growth was determined by using various concentrations of NaCl (1.0–6.0%, wt/vol) in the medium. The isolate grew in concentration ranging from 1.0% to 5.0% (wt/vol) NaCl, showing optimum growth at 3.0% (wt/vol) NaCl (Supporting materials, Fig. S3C).

The predominant cellular fatty acids identified in strain 1612 were C16:0 (25.19%), C18:1ω (24.67%), and C16:1ω (16.58%) (Supporting materials, Table S2). These profiles were consistent with those observed in most other Sulfurospirillum species, except for S. arcachonense. The polar lipids detected in strain 1612 included phosphatidylethanolamine (PE) and an unidentified lipids (AL). Two naphthoquinones compounds were identified: menaquinone with six isoprenoid units (MK-6; 4.89%) and monomethylmenaquinone with six isoprenoid units (MMK-6; 95.11%).

Classification based on phylogenomic analysis

A phylogenomic tree of existing Campylobacteria genomes was constructed to illustrate the taxonomic relatedness of the isolates (Fig. 1). Notably, Sulfurospirillum species are distinct from the other related genera like Campylobacter, Arcobacter, Helicobacter, and Wolinella, which are mostly host-associated. In contrast, reported pure cultures of Sulfurospirillum have been from diverse environments, including soil, freshwater sediment, sewage, and hydrocarbon-rich sites. In the phylogenetic tree, the newly isolated 1612 formed a clade with S. arcachonense, which was isolated from marine sediment.

Fig 1 Phylogenomic tree of Campylobacetial isolates. A total of 105 Campylobacetial reference genomes with the complete length were downloaded from NCBI. Combined with some other Sulfurospirillum isolates’ genomes, this phylogenetic tree was reconstructed using PhyloPhlAn software. Different colors correspond to the sources of isolates. The root was composed of three genomes of genus Hippea. The scale bar indicates 0.1 substitutions per amino acid position.

Phylogenetic tree displays relationships among various bacterial species, color-coded by habitat types including host-associated, marine environment, freshwater and sediment, hydrocarbon-rich, hydrothermal site, sewage, and soil.

We compared the ANI, AAI, and POCP values among the genomes of genus Sulfurospirillum. The ANI values between 1612 and other Sulfurospirillum genomes were around 71% (Supporting materials, Fig. S4; Table S3), far below the hypothesized species demarcation threshold of 95% (63), indicating that 1612 represents a distinct species. The ANI value (Fig. S4; Table S3) was also lower than the threshold for genus definition (73.98%) (64), but AAI and POCP values were around the threshold to assign 1612 into a novel genus (Supporting materials, Fig. S4; Table S4). The 63–65% of AAI value between 1612 and other strains (Supporting materials, Fig. S4; Table S4) was below the suggested 65% threshold to set a novel genus (65), but a little higher than the suggested 60% boundary mentioned in another study (66). The results of POCP values (55–63%, in Supporting materials, Table S5) were higher than 50%, a genus boundary for prokaryotic lineages (57). Considering the uncertainty of POCP value in delineating a new genus, we proposed strain 1612 to be a novel species within the genus Sulfurospirillum, named as Sulfurospirillum kueishanense sp. nov. However, the potential to establish a new genus out of genus Sulfurospirillum needs to be considered in the future.

Genomic features of the novel Campylobacteria isolate from SW-HTVs

The genome of strain 1612 (accession no. CP140614) comprised one circular chromosome with 2,377,931 bp containing 2365 open reading frames (Supporting materials, Table S6; Fig. S5). The total length of all the identified genes was 2,228,148, accounting for 93.7% of the total genomic length. The total GC-content was 37.84%. A total of 6 rRNAs (5S + 16S + 23S) and 39 tRNAs were detected. Approximately 90.19% of all identified genes, about 2,133 genes, could be annotated by the COG database. Excluding “function unknown,” the categories “energy production and conversion” and “amino acid transport and metabolism” were the most abundant, represented by about 227 and 200 genes, respectively (Supporting materials, Fig. S6). Using the KEGG database, approximately 58.01% of all identified genes, about 1372 gene, could be annotated. Most of the annotated genes were related to “metabolism” pathway (green histogram in Supporting materials, Fig. S7).

The genome of 1612 contained three genomic islands (Supporting materials, Table S6), termed GI01 (coordinates 1,478,760–1,492,033, total 13,273 bp), GI02 (coordinates 2,140,982–2,183,349, total 42,367 bp), and GI03 (coordinates 97,646–133,215, total 35,569 bp), possibly acquired via horizontal gene transfer. Most of the genes in the genomic islands encoded hypothetical proteins, but some appeared to encode the functional proteins related to the adaption mechanisms to the HTVs’ environment. For example, the GI02 included the formate transporter (gene2204, fdhC) and formate hydrogenlyase (FHL) complex (gene2205 to gene2214) (Supporting materials, Table S7), one of the hydrogenases identified in 1612.

Metabolisms and potential functions of the Campylobacteria isolate from SW-HTVs

By combining genome-scale analyses with known physiological traits, the metabolic network of strain 1612 was reconstructed (Fig. 2). The result showed that the bacteria may live a chemolithoautotrophic lifestyle with high plasticity on metabolisms, potentially contributing to several essential biogeochemical processes such as carbon fixation, hydrogen oxidation, sulfur reduction/oxidation and nitrate reduction.

Fig 2 Overview of the central metabolism pathways in strain 1612. The presence and absence (marked with question marks) of genes were predicted based on the annotations from KEGG. Created via BioRender.com. H2ase, hydrogenase; Nap, periplasmic nitrate reductase; Nrf, nitrite reductase; Nos, nitrous oxide reductase; TST, thiosulfate sulfurtransferase; Ttr, tetrathionate reductase; Sqr, sulfide-quinone oxidoreductase; Sox, sulfur-oxidation multienzyme; Psr, thiosulfate reductase; PFK, phosphofructokinase; FBP, Fructose 1,6-bisphosphate; pps, phosphoenolpyruvate synthase; PK, Pyruvate kinase.

Metabolic pathways in bacteria highlight glycolysis, gluconeogenesis, the pentose phosphate pathway, the Entner-Doudoroff pathway, TCA cycle, and respiratory chain components.

Carbon fixation in strain 1612 mainly occurs via the rTCA cycle (reductive citric acid cycle, or reductive tricarboxylic acid cycle), which reverses the reactions of the oxidative citric acid cycle (TCA cycle) and forms acetyl-CoA from two CO2 molecules (Fig. 2). Characteristic enzyme genes of the rTCA cycle, including fumarate reductase (frdAB), 2-oxoglutarate synthase (oorABCD), ATP-citrate lyase (aclAB), and pyruvate:ferredoxin oxidoreductase (porABCD), were all identified in the genome of strain 1612 (Fig. 2). The initial steps of gluconeogenesis pathway, forming fructose-6-phospahte from oxaloacetate or pyruvate (KEGG module M00003), were complete in 1612.

In addition to carbon fixation, strain 1612 also showed potential for using organic carbon. While glycolysis was incomplete due to the lack of some irreversible enzymes (such as phosphofructokinase in catalyzing phosphorylate fructose-6-phosphate to fructose-1,6-bisphosphate), the potential to degrade glucose was supplemented by the Entner-Doudoroff pathway (KEGG module M00008) in 1612.

Adenosine triphosphate (ATP) production in strain 1612 occurs through substrate-level phosphorylation and oxidative phosphorylation. Substrate level phosphorylation (KEGG module M00579) involves enzymes like phosphate acetyltransferase (Pta) producing acetyl-phosphate from acetyl-coenzyme A, and enzyme acetate kinase (AckA) producing acetate from acetyl-phosphate. Additionally, strain 1612 possesses an F-type ATPase to synthesize ATP through oxidative phosphorylation, utilizing a series of enzyme complexes in the respiratory chain and menaquinone to shuttle electrons along the electron transfer chain. The detected genes of respiratory complexes included the Complex I (NADH:quinone oxidoreductases I, nuoA ~N), Complex II (Succinate dehydrogenase/Fumarate reductase, sdhBC and frdABC), Complex III (Cytochrome bc1 reductase), Complex IV (Cytochrome caa3-type oxidase, coxABCD; and the cbb3-type, ccoNOPQ), and the Cyt bd oxidoreductase (cydAB) (Fig. 2). Both cbb3-type cytochrome c oxidase (67) and Cyt bd oxidoreductase (68) have a high affinity for oxygen and typically function under low-oxygen conditions. In contrast, cytochrome caa3-type oxidase has a low affinity for oxygen and is usually expressed and functional under high oxygen conditions in many bacterial species, such as Paracoccus denitrificans, Bradyrhizobium japonicum, and Rhodobacter sphaeroides (69–71).

Hydrogenase catalyzes reversible H2 oxidation, H2↔2H++2e− (72). In strain 1612, five [NiFe]-hydrogenase gene clusters and one hydrogenase accessory hyp gene cluster were detected (Fig. 3B). Based on the constructed phylogenetic tree, [NiFe]-hydrogenase detected in strain 1612 were identified as subgroup 1a, 1b, 2d, 4a, 4e, and 4f (Fig. 3A). The hydrogenase structural gene large subunit of subgroup 1a and 4f, corresponding to gene0156 and gene0407 in 1612, was not found in other Campylobacteria (Fig. 3A), indicating the presence of diverse hydrogenase in strain 1612.

Fig 3 Classification and phylogeny of hydrogenases identified from 1612 genome. (A) The neighbor-joining skeleton trees showed the phylogenetic relationships of hydrogenases structural gene large subunit identified from strain 1612 and previously reported bacteria (59). The central circle trees were colored by [NiFe]-hydrogenase subgroups. The black nodes indicated the well-supported nodes (bootstrap values >80). Detailed information was shown in Fig. S9; (B) The six hydrogenase clusters identified from the genome of 1612. The color of genes corresponded to the COG categories. The gene names of structural genes of hydrogenase were colored in red. The large subunit used to construct the phylogenetic tree was marked with a gene number. Each cluster included one or two groups of hydrogenases. The hydrogenase subgroups of each cluster were shown in the brackets below the cluster name. The gene hyfG of cluster 6 was annotated as hycE in KEGG, but we renamed it to “hyfG” as it belonged to group 4a hydrogenase in the phylogenetic tree.

Circular phylogenetic tree depicts different bacterial groups, color-coded and labeled with specific genes and clusters. Genetic clusters illustrate the arrangement and orientation of various genes within each cluster.

Strain 1612 also possesses genes encoding enzymes involved in oxidizing (Sqr, Sox, Ttr, and TST) or reducing (Psr) sulfur compound (Fig. 2). It also included genes for dissimilatory nitrate reduction to ammonium (DNRA) (such as nap and nrf genes involving in reducing nitrate and nitrite reduction to ammonium) and part of denitrification pathway (nos gene involving in N2O to nitrite and dinitrogen), suggesting a role in nitrogen cycling.

Physiological experiments validating the metabolic capabilities of strain 1612

The physiological experiments performed aimed to validate the annotated chemotrophic lifestyle of the isolated Campylobacteria species. The results demonstrated that strain 1612 exhibited robust growth when utilizing CO2 as the sole carbon source in conjunction with H2 as the electron donor, and either nitrate, thiosulfate or S0 as electron acceptors (Fig. 4A and B). However, no growth was observed when sulfate, sulfite or nitrite was employed as electron acceptors in conjunction with hydrogen as the sole electron donor and CO2 as the only carbon source (Fig. 4A and B). Strain 1612 also displayed the ability to utilize simple organic carbon sources, including formate, lactate, pyruvate, fumarate, succinate, and malate, as the electron donors (Fig. 4C, D and E). Each test was conducted with a single type of organic carbon added, while S0 served as the sole electron acceptor and CO2 as the only carbon source. Interestingly, when acetate was added, strain 1612 exhibited growth on formate as the sole electron donor in conjunction with nitrate, thiosulfate, or S0 as the sole electron acceptors (Supporting materials, Fig. S8A and B). Moreover, strain 1612 appeared growing better in the MMJHS medium containing both S0 and S2O32− (Supporting materials, Fig. S8D). However, no growth was observed when sucrose was used as electron donor and S0 as the electron acceptor (data not shown). The sensitivity to oxygen in strain 1612 was also tested. As shown in Fig. 4F, strain 1612 was not a strict anaerobe and can tolerate a certain level of oxygen, as evidenced by better growth in the absence of an oxygen scavenger.

Fig 4 The growth curves for the substrate testing. The spectrum of electron acceptor testing when using H2 as the sole electron donor (A&B). The electron donor testing was proceeded by adding lactate, pyruvate, fumarate, succinate, malate, or formate when S0 was the sole electron acceptor (C–E). The growth profile of 1612 under different concentrations of oxygen scavenger (F). All the “control” in the figure were referred to the negative control, meaning without bacterial inoculum. The electron donor and electron acceptor of the negative control for each testing were listed in pairs: H2/S0 (A), H2/NO3− (B), lactate/S0 (C), malate/S0 (D), H2/S0 (E), and H2/(S0, S2O32−, NO3−) (F).

Six line graphs plot bacterial growth over time under different conditions: sulfur compounds, nitrogen compounds, lactate and pyruvate, malate, succinate, and fumarate, hydrogen and formate, and varying scavenger concentrations, each with a control.

Pan-genome analysis and functional traits for bacterial adaptation in SW-HTVs

We analyzed 13 complete genomes of Sulfurospirillum strains, including the newly isolated 1612 in this study. These strains were sourced from diverse origins. All the genomes were utilized to construct the genomic phylogenetic tree (Fig. 1). The anvi’o pan-genomic pipeline distinguished 7,553 gene clusters (GCs) out of a total of 35,485 genes across the 13 genomes. The highest numbers of GCs were observed in the core genes shared among all genomes, followed by the unique genes of S. arcachonense, 1612, S. multivorans, and S. cavolei (Fig. 5A). Most of the core genes were associated with essential functions, such as energy production/conversion, translation, and amino acid transport/metabolism (Fig. 5C). The pan-genome profile illustrated that as the number of sequenced genomes increased, the pan-genome size expanded rather than reaching a plateau (Fig. 5B). Consequently, it can be inferred that the pan-genome of the genus Sulfurospirillum followed an open type, indicating a strong potential for horizontal gene transfer events. Considering the potential of 1612 to be a novel genus, we recalculated the pan-genome profile of genus Sulfurospirillum without including strain 1612. The result still indicated an open-type pan-genome (Supporting materials, Fig. S10).

Fig 5 Pan-genome analysis results. (A) Upset plot to highlight the number of genes sharing among genomes. The variable component of genomes contributing to a collection of gene clusters was marked by black dots. (B) The Pan-genome profile trends of the genus Sulfurospirillum. (C) The COG annotation results for the shared core genes.

Bar graph depicts the intersection size of gene families across different species. Line graph plots the number of total and core gene families versus the number of genomes. Pie chart plots the distribution of gene functions across various categories.

To discern the functional distinction between 1612 and other strains, we identified five bins including the core gene bin, two accessory gene bins, and two unique genes bins using the visualization platform of the Anvi’o pan-genomic workflow (Fig. 6A). The “core genes” comprised 12,889 genes belonging to 957 GCs (Fig. 6A). Most of its genes could be annotated by KEGG (79%) or COG (96%), while the ratio of annotated genes was much lower for the other bins (Supporting materials, Table S8). The GCs number of the “unique genes of 1612” bin and “unique genes of S. arcachonense” bin was comparable to the “core genes” (Fig. 6A). Protein families involved in “protein families: signaling and cellular processes,” “signal transduction,” “membrane transport,” and some “poorly characterized” genes had higher ratios in the unique genes of 1612 and S. arcachonense compared to the core genes (Fig. 6B). At a deeper functional annotation level, signal-related processes including “transporters,” “two-component system,” and “ABC transporters” were much more prevalent in the unique and accessory genes than the core genes (Supporting materials, Fig. S11). At the gene level, we also analyzed the KEGG orthology (KO) detected only in 1612. It was found that 1612 had specific ABC transporters of dipeptide, phosphate, and thiamine (Fig. 6B, purple block). Furthermore, several genes related to biofilm formation and benzoate/toluene degradation were exclusively detected in the unique genes of 1612 (Fig. 6B, blue and green blocks). Regarding energy metabolism, some genes associated with sulfur, nitrogen, methane, oxygen, and hydrogen were only detected in the “unique genes of 1612” (Fig. 6B, orange block, supporting materials, Table S9).

Fig 6 Pan-genome analysis for 13 Sulfurospirillum genomes using the visualization platform of Anvi’o. (A) The first 13 layers, in which dark color indicated the presence of a gene cluster and light color indicated its absence, represented individual genomes. ANI values among different genomes were represented on heatmap determined from the high similarity (red) and low similarity (white). Five specific groups of gene clusters were selected into bins such as the “core genes,” which were indicated in different colors in the outmost layer. The numbers of gene clusters and responding genes included in each bin were listed under the bin names, and separated with “/.” The dendrograms on the top represented the hierarchical clustering of genomes based on the frequencies of gene clusters. (B) The breakdown of gene annotation results shows the component of different KEGG level 1 categories in total gene numbers for each selected bin. Genes without annotated information were not shown in the figure but used in the calculation of relative abundance. Some annotation categories of “unique genes of 1612” bin were extended to show the information of its unique KOs.

Gene distribution among bacterial species includes unique genes, core genes, and accessory genes. Heatmap and table detail gene functions across different metabolic and cellular processes, including specific gene names and their roles.

The annotation results of COG were also carefully analyzed between the “unique genes of 1612” and “core gene” bins. The category “P” related to inorganic ion transport and metabolism was significantly higher in “unique genes of 1612” than in “core gene” bins (Supporting materials, Fig. S12B). A total of 72 COGs showed significant differences between these two bins (Supporting materials, Table S10). For “unique gene of 1612,” a significantly higher abundance of COGs was found in metabolism-related categories (such as P, C, and G). Specifically, rhodanese-related sulfurtransferase (COG0607 and COG2897), cytochrome c nitrite reductase (COG3303), Fe (COG4773 and COG1918), and phosphate (COG0226, COG0581, and COG0573) were abundant in the category “P.” In the category “C” related to energy production and conversion, hydrogenlyase (COG0650, COG0651, and COG5557), cytochrome c553/551/552 (COG2863 and COG4654), and napC (COG3005) were notably abundant. Additionally, the category “G” associated with carbohydrate transport and metabolism exhibited higher abundance in TRAP-type C4-dicarboxylate (e.g., succinate, fumarate, l-malate, and d-malate) transport system (COG1593 and COG3090), glycosyltransferase (COG1216), glycogen synthesis (COG0297 and COG0296), and alpha-amylase/alpha-mannosidase, GH57 family (COG 1449) (Supporting materials, Table S10; Fig. S13).

Comparison of key energy metabolic genes

Genes involved in energy metabolisms were compared among genomes, including 13 genomes of Sulfurospirillum and 6 genomes of mesophilic Campylobacterial isolates from deep-sea HTVs (Fig. 7). The genome of strain 1612 harbored characteristic enzyme genes of rTCA cycle for CO2 fixation (encoded by genes aclAB and oorABCD), sulfur oxidation system SOX, hydrogenase, and nitrate reductase, supporting the potential of chemoautotrophy. The genes of CO2 fixation (aclAB) and sulfur oxidation system (soxABCDYZ) were presented in the genomes of all collected isolates from deep-sea HTVs and 1612. Except for 1612, the potential of rTCA cycle for CO2 fixation (aclAB) was also detected in S. cavolei, but was absent in other Sulfurospirillum species.

Fig 7 Color-coded table indicating major functional genes and their abundance in 13 genomes of Sulfurospirillum and 6 genomes of mesophilic chemolithotrophic Campylobacteria. The y-axis indicated the genes likely involved in energy metabolism, and the x-axis indicated the strain designations. Red color corresponded to the complete genes identified in carbon, sulfur, and nitrogen nutrient cycles; pink indicated the partial identified in the gene cluster; white indicated an absence of relevant genes.

Heatmap depicts the presence of various genes in different bacterial species. Gene categories include sulfide oxidation, nitrogen fixation, acetate production, fermentation, glycolysis, cytochrome oxidase, and hydrogenase groups.

DISCUSSION

Our novel isolate 1612 from the Kueishan Island SW-HTVs enriched the species diversity of Campylobacteria. The two major orders, Nautiliales and Campylobacterales, within Campylobacteria, exhibited clear differences in their metabolic potential and core bioenergetics (6). Genera belonging to the order Nautiliales, such as Nautilia, Caminibacter, Lebetimonas, Cetia, Nitrosophilus, and Nitratiruptor, are commonly found in deep-sea HTVs, according to data from LPSN (List of Prokaryotic names with Standing in Nomenclature). In contrast, most cultured representatives of the order Campylobacterales thrive at lower temperatures and can typically tolerate a much greater range of oxygen concentrations (15). Genera Sulfurovum and Sulfuromonas are the frequently studied Campylobacterales in vent area, with few reports of genera Nitratifractor and Hydrogenimonas, based on data from LPSN. Although the initial enrichments of Campylobacteria isolated from deep-sea HTVs (including Am-H and Ex-18.2) (73) are recognized as Sulfurospirillum and have been tested to be moderately thermophilic sulfur-reducing heterotrophs or sulfur lithoautrotrophs using hydrogen as the electron donor, their draft genomes do not provide a basis for deeper comparison with other typical strains of Sulfurospirillum. In this study, our novel isolate (strain 1612) of genus Sulfurospirillum expands the strain diversity of Campylobacterales in hydrothermal vent areas. Through systematic genomic and phenotypic study, these results may contribute to a better understanding of mesophilic Campylobacteria in hydrothermal vent systems.

The ability to fix CO2 was only detected in strain 1612, but not in other Sulfurospirillum species. Genus Sulfurospirillum consists of versatile, often microaerophilic bacteria, capable of growing with various growth substrates (74), including hydrogen, formate, nitrate, sulfur compounds, and many toxic compounds (e.g., PCE, toxic arsenate, and selenate). Upon thorough comparison of phenotypic and genomic features within genus Sulfurospirillum, we found that the capacity for chemoautotrophy appears to be a unique metabolism for the novel Sulfurospirillum strain from the Kueishan Island SW-HTVs (strain 1612) (Table 1; Fig. 7). Although chemoautotrophic potential has been partially shown by detection of the aclAB gene in the genome of S. cavolei (75), our study is the first to clearly demonstrate the chemolithoautrophic capacity of the genus Sulfurospirillum via combined results of genomic features and experimental validations. Similar to the prevalent Campylobacterial isolated from deep-sea HTVs (76), such as the mesophilic Sulfurovum and Sulfurimonas, chemolithoautotrophy also appears to be a common feature of Campylobacteria in SW-HTVs. On the other side, considering the unique ability to fix CO2 in genus Sulfurospirillum and the monophyly formed by strain 1612 and S. arcachonense in the phylogenomic tree (Fig. 1), the potential to establish strain 1612 and our other isolates into a new genus out of genus Sulfurospirillum needs to be considered in future studies.

TABLE 1 Differential phenotypic characteristics of 1612 and other Sulfurospirillum speciesb,c

Characteristics	1	2	3	4	5	6	7	8	9	10	
Morphology	Slightly curved rods	Curved spiral rods	Curved/ helical rods	Vibrioid to spirillum	Vibrioid to spirillum	Vibrioid to spirillum	Slightly curved rods	Vibrioid to spirillum	Slightly curved rods	Vibrioid to spirillum	
Width (µm)	0.3–0.6	0.3–0.5	0.6	0.3	0.3	0.3	0.6	0.5–0.7	0.4–0.5	0.3	
Length (µm)	1.0–2.5	1.0–3.0	2.5–4.0	1.0–2.5	1–2	1–2	2.5–4.0	1.2–2.5	1.5–5.0	2.0–3.0	
Presence of flagella	+	+	+	+	+	+	+	+	+	+	
Colony color	Whitish	Yellow	ND	Whitish	ND	ND	ND	ND	ND	Greyish	
Chemoautotrophy	+	−	−	−	−	−	−	−	−	−	
Growth rage (optimum)	
 Temperature (°C)	20–37 (28)	20–36 (ND)	ND (30)	8–30 (26)	ND (33)	ND (20)	ND (25–30)	ND	20–40 (30)	9.9–31.9 (30.2)	
 NaCl (%)	2.0–4.0 (3.0)	<0.2 (ND)	ND (ND)	0.6–4.0 (1.2–2.0)	ND (0.05)	ND (0.1)	ND (ND)	1.3–3.5 (1.3)	<1.0 (0)	1.0–5.5 (2.0)	
 pH	5.5–8.0 (6.5)	ND (7.0–7.1)	ND (7.0–7.5)	6.1–8.2 (7.0–7.4)	ND (7.5)	ND (7.5)	ND (ND)	7.1–9.7 (8.5)	6.0–8.0 (7.0)	5.9–8.5 (7.3)	
Electron acceptors	
 Nitrate	+	+	+	−	+	+	+	+	+	+	
 Nitrite	−	+	−	−	+	+	+	+	−	+	
 Oxygen
 (microaerobic)	+a	+	ND	+	+	+	+	−	+	+	
 S0	+	+	ND	+	+	+	+	+	+	+	
 Sulfite	−	+	−	−	−	ND	−	ND	+	+	
 Thiosulfate	+	+	−	−	+	+	−	+	+	+	
 Sulfate	−	−	−	−	ND	ND	−	ND	−	ND	
 Arsenate	ND	−	+	−	+	+	+	+	+	+	
 PCE	−	ND	+	ND	−	ND	+	ND	−	+	
 Selenate	ND	−	+	ND	+	−	+	−	−	+	
 DMSO	ND	+	ND	−	−	ND	ND	−	+	+	
Electron donors	
 Hydrogen	+	+	+	+	+	+	+	+	+	+	
 Formate	+	+	+	ND	+	+	+	+	+	+	
 Fumarate	+	+	+	+	+	+	ND	+	+	+	
 Succinate,
 malate,	+	+	ND	+	ND	ND	ND	ND	+	+	
 Pyruvate	+	+	+	+	+	+	+	+	+	+	
 CO	ND	ND	ND	−	ND	ND	ND	ND	ND	+	
 Lactate	+	−	+	+	+	+	+	+	+	+	
Major quinone	MK-6 and MMK-6	MK-6 and MMK-6	MK	MK-6 and MMK-6	ND	ND	ND	ND	ND	ND	
Isolated from	Sulfur-rich sediment of shallow-water HTVs	Mud of forest pond	Activated sludge	Surface sediment of an intertidal mud flat	Freshwater marsh	Freshwater marsh	Soil	Bioreactor	Groundwater storaging crude oil	Marine sulphide-rich seep	
DNA G + C content (mol%)	37.8	40.6	41.5	32	40.8	40.9	41.8	47.6	42.7	30.5	
Gene bank accession no.	CP140614	GCA_000024885	GCA_000568815	GCA_000597725 (Draft)	GCA_000265295	GCF_000813345 (Draft)	GCA_001723605	ND	GCA_0015 48055	ND	
a Strain 1612 is a facultative aerobe. Its sensitivity to oxygen is evaluated by a simple test involving the addition of oxygen scavenger to the medium (results see Fig. 4F).

b Strains: 1, 1612 (this study); 2, S. deleyianum DSM 6946T (77, 78); 3, S. multivorans DSM 12446T (79, 80); 4, S. arcachonense DSM 9755T (81); 5 and 6, S. barnesii DSM 10660T and S. arsenophilum DSM 10659T, respectively (82); 7, S. halorespirans DSM 13726T (80); 8, S. alkalitolerans DSM 24537T (83); 9, S. cavolei DSM 18149 (75); and 10, ‘S. carboxydovorans’ DSM 16295 (84). +, positive; −, negative; ND, not data. All strains are able to utilize the following as electron donors: formate (+ acetate), hydrogen (+ acetate) and pyruvate. All strains are positive for the fermentation of fumarate.

c Abbreviations: MK, Menaquinon; MMK, monomethylmenaquinone.

Another important functional feature of strain 1612 is the presence of diverse subgroups of [NiFe] hydrogenases in its genome. Three major classes of hydrogenase, featured with a [NiFe], a [FeFe], or a [Fe] core, have been identified. Based on the phylogeny of the large subunits, [NiFe] hydrogenases are further classified into four groups: the membrane-bound and cytosolic H2-uptake [NiFe] hydrogenases (Group 1 and Group 2), bidirectional heteromultimeric cytosolic [NiFe]-hydrogenases (Group 3), and membrane-bound H2-evolving [NiFe] hydrogenases (Group4) (85). In strain 1612, the detected phylotypes of [NiFe]-hydrogenase belong to subgroups 1b, 2d, 4a, 4c, 4e, 1a, and 4f (Fig. 3A), suggesting that both H2-uptaking and evolving potentials present in the bacterium.

Hydrogen oxidation is widely distributed in Campylobacteria (58). As demonstrated in Wolinella succinogenes (86), most Campylobacteria only contain Group 1b hydrogenase genes (58). The newly isolated strain 1612 harbors both subgroup 1a, 1b and cytosolic H2-uptake subgroup 2d hydrogenase (Fig. 3). Subgroup 2d hydrogenase inferred to generate reductant for carbon fixation via rTCA (87). The observed diverse phylotypes of H2 utilization hydrogenase genes might be essential for strain 1612’s survival in hydrothermal vent ecosystem, as hydrogen oxidation is among the chemosynthetic reactions providing the greatest energy yields (88). Moreover, hydrogen oxidation seems to be a favorable energy source for autotrophic carbon fixation compared to the oxidation of sulfide/thiosulfate, despite considerably more energy being yielded through the oxidation of sulfide/thiosulfate than through hydrogen oxidation (free standard enthalpies are −797 kJ/mol H2S vs −237 kJ/mol H2 with O2 as electron acceptor) (89). This is because the redox potential of hydrogen is more negative than that of the reducing equivalent NAD(P)/H; in contrast to sulfide, a reverse electron transport is not required in conjunction with hydrogen oxidation. Thus, only a third of the energy is required for fixing 1 mol of carbon when oxidizing hydrogen compared to sulfide (1,060 kJ for hydrogen vs 3,500 kJ for sulfide) (90).

The Group 4a (Hyf) and its homologous Group 3 (Hyc) hydrogenase both function in the FHL complex during fermentation, which produces H2 in Escherichia coli (85). Similar with Group 4e (Ech) [NiFe]-hydrogenase, they are absent in genomes of typical Campylobacteria strains isolated from deep-sea HTVs, such as those belonging to genera Sulfurovum, Sulfuromonas, and Nitratifractor (Fig. 7). In contrast, the existence of hyf genes (Group 4a) have been reported in the most genomes of Sulfurospirillum, such as S. multivorans (91) and Sulfurospirillum diekertiae (92). The H2 production activities have also been demonstrated in Sulfurospirillum species, that is, S. multivorans, S. cavolei, S. deleyianum, and S. arsenophilum, which can fermentative grow on pyruvate and produce H2 after about 20 generations of transfer from the initial cultures (93). The similar structures of hyf gene clusters between strain 1612 and Sulfurospirillum species (93) support the possibilities of H2 production in strain 1612, yet the detailed mechanisms and in-situ activities of H2 production need to be further studied in the future.

Diverse respiratory oxygen reductases (terminal oxidases) are found in Campylobacteria isolated from SW-HTVs. In vent areas, Campylobacteria occupy high-sulfide, low-oxygen niche compared to chemoautotrophic Gammaproteobacteria (94, 95). Considering the sharp chemical and physical gradients formed during the mixing of reduced hydrothermal fluids with oxic seawater, the importance of oxygen for Campylobacteria is underscored. Recent short-term incubations have indicated the importance of oxygen in structuring natural communities and affecting carbon fixation efficiency for Campylobacteria (96). There are three enzyme superfamilies capable of acting as terminal respiratory oxygen reductases—heme-copper oxygen reductases, alternative oxidases, and Cyt bd oxidoreductases (97). Most terminal oxidases belong to the heme-copper oxygen reductases superfamily, such as cytochrome c oxidase (Cox), which uses a c-type cytochrome as its electron donor. Cox emzyme are universal in aerobic and facultative aerobic organisms, represented by caa3-type and cbb3-type oxidases (98). For our strain 1612 isolated from SW-HTVs, genes for respiratory oxygen reductases with both high (cbb3-type Cox and Cyt bd oxidoreductase) and low (caa3-type Cox) affinity for oxygen were identified in its genome (Fig. 2). This indicates that strain 1612 is a facultative aerobe, a conclusion supported by the physiological test with oxygen scavenger addition (Fig. 4F). In comparison, most isolates of genera Sulfurovum and Sulfurimonas from deep-sea HTVs only have the genes ccoNOPQ of cbb3-type Cox (Fig. 7). Higher concentrations of oxygen are reported in SW-HTVs than in their deep-sea counterparts (4, 99). This may be due to more effective mixing with oxygen-rich surface waters and the influence of tidal and wave actions. In Kueishan Island, high temporal variations in water temperature are attributable to diurnal tides (100). It is, therefore, reasonable to expect a similar variation in other parameters within SW-HTVs (99). The presence of more types of oxygen reductases in strain 1612 confirms the higher variable in oxygen concentration in SW-HTVs compared to the deep-sea environment.

Through our study, we identified unique genes in strain 1612 by comparing its genome with the ones of non-vent strains of Sulfurospirillum (Fig. 6B). These unique functional traits of strain 1612, possibly gained via genome extension events (e.g., gene horizontal transfer). Indeed, at least three genomic islands in genomes of 1612 (Supporting materials, Table S6) were found to correspond to the presence of horizontal gene transfer events. Such genomic plasticity increases micro-diversity and confers a competitive advantage enabling this lineage to thrive in changing physical–chemical gradients of hydrothermal systems (101). Reflecting the complex environmental conditions within SW-HTVs, the diverse signal genes identified in unique genes of strain 1612 may facilitate its adapting and dominance in vent area. This finding prompts us to focus on specific sources (e.g., dipeptide, phosphate, and thiamine) within the vent area to gain a deeper understanding of Campylobacteria. However, due to limitations in genome numbers and functional annotation, the results still lack a refined selection of target genes. Further physiological experiments involving transcriptome analysis and gene knockout techniques are necessary to elucidate the details of the adaptation mechanisms.

Description of Sulfurospirillum kueishanense sp. nov

Sulfurospirillum kueishanense (kuei.shan.en’ se. N.L. neut. adj. kueishanense, pertaining to the Kueishan Island, the western Pacific Ocean, the locality from which the strain was isolated). Slightly curved cells are 0.3–0.6 µm wide and 1.0–2.5 µm long. Gram-negative. Motile by polar flagellum. Colonies are whitish and round. The pH range is 5.5–8.0; optimum, 6.5. The temperature range is 20–37°C; growth optimum, 28°C. The salinity range is 2.0–4.0% NaC1; optimum growth occurs in the presence of 3.0% NaCl. Chemolithoautotrophy. Uses thiosulfate (incomplete reduction), sulfur, and nitrate as electron acceptors. Utilizes formate and H2, pyruvate, and fumarate as electron donors. MK-6 and MMK-6 are the sole quinones, with MMK-6 predominating. The predominant fatty acids in the membrane lipids include C18:1 ω7c, C16:0, and C18:0, with seldom C13:0 and C15:1. The GC content of the genomic DNA is 37.8 mol%. The type strain is 1612 (MCCC 1K08293). Isolated from the sulfur-rich sediment of SW-HTVs (Kueishan Island, western Pacific).

ACKNOWLEDGMENTS

This work was supported by the National Natural Science Foundation of China (grant nos. 41906134, 42276149, and 92251303). The authors thank students and technicians in J.-S.H. Laboratory of National Taiwan Ocean University for their assistance with sample collection. The authors thank all the reviewers for their suggestions and advice to improve the manuscript.

L.W. designed the sampling scheme, summarized the data, and wrote the original draft. R.L. guided the sampling scheme and revised the manuscript. L.W. and Y.G. performed bioinformatical analyses including sequence processing, genomic comparison, and phylogenetic analysis. X.C. assisted by M.S., conducted experimental procedures including DNA extraction, PCR amplification, and physiological testing. J.C. helped to design the physiological testing scheme. J.-S.H. managed the sampling. J.F. made comments and suggestions to the text. L.W., R.L., and J.F. also provided funding support for this work. All authors contributed to the article and approved the submitted version.

SUPPLEMENTAL MATERIAL

The following material is available online at https://doi.org/10.1128/msystems.00148-24.

10.1128/msystems.00148-24.SuF1 Supplemental figures msystems.00148-24-s0001.docx

Fig. S1 to S13.

10.1128/msystems.00148-24.SuF2 Supplemental tables msystems.00148-24-s0002.xlsx

Tables S1 to S10.

ASM does not own the copyrights to Supplemental Material that may be linked to, or accessed through, an article. The authors have granted ASM a non-exclusive, world-wide license to publish the Supplemental Material files. Please contact the corresponding author directly for reuse.
==== Refs
REFERENCES

1 Campbell BJ, Engel AS, Porter ML, Takai K. 2006. The versatile epsilon-proteobacteria: key players in sulphidic habitats. Nat Rev Microbiol 4 :458–468. doi:10.1038/nrmicro1414 16652138
2 Dick GJ. 2019. The microbiomes of deep-sea hydrothermal vents: distributed globally, shaped locally. Nat Rev Microbiol 17 :271–283. doi:10.1038/s41579-019-0160-2 30867583
3 Jannasch HW. 1984. Chemosynthesis: the nutritional basis for life at deep-sea vents. Oceanus 27 :73–78.
4 Tarasov VG, Gebruk AV, Mironov AN, Moskalev LI. 2005. Deep-sea and shallow-water hydrothermal vent communities: two different phenomena? Chem Geol 224 :5–39. doi:10.1016/j.chemgeo.2005.07.021
5 Price RE, Giovannelli D. 2017. A review of the geochemistry and microbiology of marine shallow-water hydrothermal vents, p 1–29. In Reference module in earth systems and environmental sciences. Elsevier.
6 Waite DW, Vanwonterghem I, Rinke C, Parks DH, Zhang Y, Takai K, Sievert SM, Simon J, Campbell BJ, Hanson TE, Woyke T, Klotz MG, Hugenholtz P. 2017. Comparative genomic analysis of the class Epsilonproteobacteria and proposed reclassification to Epsilonbacteraeota (phyl. nov.). Front Microbiol 8 :682. doi:10.3389/fmicb.2017.00682 28484436
7 Reysenbach AL, Longnecker K, Kirshtein J. 2000. Novel bacterial and archaeal lineages from an in situ growth chamber deployed at a Mid-Atlantic Ridge hydrothermal vent. Appl Environ Microbiol 66 :3798–3806. doi:10.1128/AEM.66.9.3798-3806.2000 10966393
8 Polz MF, Cavanaugh CM. 1995. Dominance of one bacterial phylotype at a Mid-Atlantic Ridge hydrothermal vent site. Proc Natl Acad Sci U S A 92 :7232–7236. doi:10.1073/pnas.92.16.7232 7543678
9 Moyer CL, Dobbs FC, Karl DM. 1995. Phylogenetic diversity of the bacterial community from a microbial mat at an active, hydrothermal vent system, Loihi Seamount, Hawaii. Appl Environ Microbiol 61 :1555–1562. doi:10.1128/aem.61.4.1555-1562.1995 7538279
10 Muto H, Takaki Y, Hirai M, Mino S, Sawayama S, Takai K, Nakagawa S. 2017. A simple and efficient RNA extraction method from deep-sea hydrothermal vent chimney structures. Microbes Environ 32 :330–335. doi:10.1264/jsme2.ME17048 29187693
11 Nakagawa T, Takai K, Suzuki Y, Hirayama H, Konno U, Tsunogai U, Horikoshi K. 2006. Geomicrobiological exploration and characterization of a novel deep-sea hydrothermal system at the TOTO caldera in the Mariana Volcanic Arc. Environ Microbiol 8 :37–49. doi:10.1111/j.1462-2920.2005.00884.x 16343320
12 Campbell BJ, Polson SW, Zeigler Allen L, Williamson SJ, Lee CK, Wommack KE, Cary SC. 2013. Diffuse flow environments within basalt- and sediment-based hydrothermal vent ecosystems harbor specialized microbial communities. Front Microbiol 4 :182. doi:10.3389/fmicb.2013.00182 23898323
13 Longnecker K, Reysenbach A-L. 2001. Expansion of the geographic distribution of a novel lineage of ɛ-Proteobacteria to a hydrothermal vent site on the Southern East Pacific Rise. FEMS Microbiol Ecol 35 :287–293. doi:10.1111/j.1574-6941.2001.tb00814.x 11311439
14 Nakagawa S, Takai K, Inagaki F, Hirayama H, Nunoura T, Horikoshi K, Sako Y. 2005. Distribution, phylogenetic diversity and physiological characteristics of ɛ-Proteobacteria in a deep-sea hydrothermal field. Environ Microbiol 7 :1619–1632. doi:10.1111/j.1462-2920.2005.00856.x 16156735
15 Sievert S, Vetriani C. 2012. Chemoautotrophy at deep-sea vents: past, present, and future. oceanog 25 :218–233. doi:10.5670/oceanog.2012.21
16 Yamamoto M, Takai K. 2011. Sulfur metabolisms in epsilon- and gamma-Proteobacteria in deep-sea hydrothermal fields. Front Microbiol 2 :192. doi:10.3389/fmicb.2011.00192 21960986
17 Xie S, Wang S, Li D, Shao Z, Lai Q, Wang Y, Wei M, Han X, Jiang L. 2019. 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 71 :004748. doi:10.1099/ijsem.0.004748
18 Mori K, Yamaguchi K, Hanada S. 2018. Sulfurovum denitrificans sp. nov., an obligately chemolithoautotrophic sulfur-oxidizing epsilonproteobacterium isolated from a hydrothermal field. Int J Syst Evol Microbiol 68 :2183–2187. doi:10.1099/ijsem.0.002803 29757127
19 Giovannelli D, Chung M, Staley J, Starovoytov V, Le Bris N, Vetriani C. 2016. Sulfurovum riftiae sp. nov., a mesophilic, thiosulfate-oxidizing, nitrate-reducing chemolithoautotrophic epsilonproteobacterium isolated from the tube of the deep-sea hydrothermal vent polychaete Riftia pachyptila. Int J Syst Evol Microbiol 66 :2697–2701. doi:10.1099/ijsem.0.001106 27116914
20 Mino S, Kudo H, Arai T, Sawabe T, Takai K, Nakagawa S. 2014. Sulfurovum aggregans sp. nov., a hydrogen-oxidizing, thiosulfate-reducing chemolithoautotroph within the Epsilonproteobacteria isolated from a deep-sea hydrothermal vent chimney, and an emended description of the genus Sulfurovum. Int J Syst Evol Microbiol 64 :3195–3201. doi:10.1099/ijs.0.065094-0 24966202
21 Inagaki F, Takai K, Nealson KH, Horikoshi K. 2004. Sulfurovum lithotrophicum gen. nov., sp. nov., a novel sulfur-oxidizing chemolithoautotroph within the epsilon-Proteobacteria isolated from Okinawa Trough hydrothermal sediments. Int J Syst Evol Microbiol 54 :1477–1482. doi:10.1099/ijs.0.03042-0 15388698
22 Wang S, Jiang L, Hu Q, Cui L, Zhu B, Fu X, Lai Q, Shao Z, Yang S. 2021. 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 12 :626705. doi:10.3389/fmicb.2021.626705 33717015
23 Hu Q, Wang S, Lai Q, Shao Z, Jiang L. 2021. Sulfurimonas indica sp. nov., a hydrogen- and sulfur-oxidizing chemolithoautotroph isolated from a hydrothermal sulfide chimney in the Northwest Indian Ocean. Int J Syst Evol Microbiol 71 :004575. doi:10.1099/ijsem.0.004575
24 Takai K, Suzuki M, Nakagawa S, Miyazaki M, Suzuki Y, Inagaki F, Horikoshi K. 2006. Sulfurimonas paralvinellae sp. nov., a novel mesophilic, hydrogen- and sulfur-oxidizing chemolithoautotroph within the Epsilonproteobacteria isolated from a deep-sea hydrothermal vent polychaete nest, reclassification of Thiomicrospira denitrificans as Sulfurimonas denitrificans comb. nov. and emended description of the genus Sulfurimonas. Int J Syst Evol Microbiol 56 :1725–1733. doi:10.1099/ijs.0.64255-0 16901999
25 Inagaki F, Takai K, Kobayashi H, Nealson KH, Horikoshi K. 2003. Sulfurimonas autotrophica gen. nov., sp. nov., a novel sulfur-oxidizing epsilon-proteobacterium isolated from hydrothermal sediments in the Mid-Okinawa Trough. Int J Syst Evol Microbiol 53 :1801–1805. doi:10.1099/ijs.0.02682-0 14657107
26 Alain K, Querellou J, Lesongeur F, Pignet P, Crassous P, Raguénès G, Cueff V, Cambon-Bonavita M-A. 2002. Caminibacter hydrogeniphilus gen. nov., sp. nov., a novel thermophilic, hydrogen-oxidizing bacterium isolated from an East Pacific Rise hydrothermal vent. Int J Syst Evol Microbiol 52 :1317–1323. doi:10.1099/00207713-52-4-1317 12148646
27 Grosche A, Sekaran H, Pérez-Rodríguez I, Starovoytov V, Vetriani C. 2015. Cetia pacifica gen. nov., sp. nov., a chemolithoautotrophic, thermophilic, nitrate-ammonifying bacterium from a deep-sea hydrothermal vent. Int J Syst Evol Microbiol 65 :1144–1150. doi:10.1099/ijs.0.000070 25604337
28 Takai K, Hirayama H, Nakagawa T, Suzuki Y, Nealson KH, Horikoshi K. 2005. Lebetimonas acidiphila gen. nov., sp. nov., a novel thermophilic, acidophilic, hydrogen-oxidizing chemolithoautotroph within the 'Epsilonproteobacteria', isolated from a deep-sea hydrothermal fumarole in the Mariana Arc. Int J Syst Evol Microbiol 55 :183–189. doi:10.1099/ijs.0.63330-0 15653874
29 Pérez-Rodríguez I, Ricci J, Voordeckers JW, Starovoytov V, Vetriani C. 2010. Nautilia nitratireducens sp. nov., a thermophilic, anaerobic, chemosynthetic, nitrate-ammonifying bacterium isolated from a deep-sea hydrothermal vent. Int J Syst Evol Microbiol 60 :1182–1186. doi:10.1099/ijs.0.013904-0 19667392
30 Alain K, Callac N, Guégan M, Lesongeur F, Crassous P, Cambon-Bonavita MA, Querellou J, Prieur D. 2009. Nautilia abyssi sp. nov., a thermophilic, chemolithoautotrophic, sulfur-reducing bacterium isolated from an East Pacific Rise hydrothermal vent. Int J Syst Evol Microbiol 59 :1310–1315. doi:10.1099/ijs.0.005454-0 19502307
31 Smith JL, Campbell BJ, Hanson TE, Zhang CL, Cary SC. 2008. Nautilia profundicola sp. nov., a thermophilic, sulfur-reducing epsilonproteobacterium from deep-sea hydrothermal vents. Int J Syst Evol Microbiol 58 :1598–1602. doi:10.1099/ijs.0.65435-0 18599701
32 Nakagawa S, Takai K, Inagaki F, Horikoshi K, Sako Y. 2005. Nitratiruptor tergarcus gen. nov., sp. nov. and Nitratifractor salsuginis gen. nov., sp. nov., nitrate-reducing chemolithoautotrophs of the epsilon-Proteobacteria isolated from a deep-sea hydrothermal system in the Mid-Okinawa Trough. Int J Syst Evol Microbiol 55 :925–933. doi:10.1099/ijs.0.63480-0 15774687
33 Nakagawa S, Inagaki F, Takai K, Horikoshi K, Sako Y. 2005. Thioreductor micantisoli gen. nov., sp. nov., a novel mesophilic, sulfur-reducing chemolithoautotroph within the epsilon-Proteobacteria isolated from hydrothermal sediments in the Mid-Okinawa Trough. Int J Syst Evol Microbiol 55 :599–605. doi:10.1099/ijs.0.63351-0 15774631
34 Mino S, Shiotani T, Nakagawa S, Takai K, Sawabe T. 2021. Hydrogenimonas urashimensis sp. nov., a hydrogen-oxidizing chemolithoautotroph isolated from a deep-sea hydrothermal vent in the Southern Mariana Trough. Syst Appl Microbiol 44 :126170. doi:10.1016/j.syapm.2020.126170 33340909
35 Takai K, Nealson KH, Horikoshi K. 2004. Hydrogenimonas thermophila gen. nov., sp. nov., a novel thermophilic, hydrogen-oxidizing chemolithoautotroph within the ε-Proteobacteria, isolated from a black smoker in a Central Indian Ridge hydrothermal field. Int J Syst Evol Microbiol 54 :25–32. doi:10.1099/ijs.0.02787-0 14742455
36 Yang TF, Lan TF, Lee H-F, Fu C-C, Chuang P-C, Lo C-H, Chen C-H, Chen C-TA, Lee C-S. 2005. Gas compositions and helium isotopic ratios of fluid samples around Kueishantao, NE offshore Taiwan and its tectonic implications. Geochem. J 39 :469–480. doi:10.2343/geochemj.39.469
37 Zeng Z, Chen C-TA, Yin X, Zhang X, Wang X, Zhang G, Wang X, Chen D. 2011. Origin of native sulfur ball from the Kueishantao hydrothermal field offshore northeast Taiwan: evidence from trace and rare earth element composition. J Asian Earth Sci 40 :661–671. doi:10.1016/j.jseaes.2010.10.019
38 Tang K, Liu K, Jiao N, Zhang Y, Chen C-T. 2013. Functional metagenomic investigations of microbial communities in a shallow-sea hydrothermal system. PLoS One 8 :e72958. doi:10.1371/journal.pone.0072958 23940820
39 Wang L, Shen Z, Cheng X, Hwang J-S, Guo Y, Sun M, Cao J, Liu R, Fang J. 2022. Metagenomic insights into the functions of microbial communities in sulfur-rich sediment of a shallow-water hydrothermal vent off Kueishan Island. Front Microbiol 13 :992034. doi:10.3389/fmicb.2022.992034 36532441
40 Li Y, Tang K, Zhang L, Zhao Z, Xie X, Chen C-TA, Wang D, Jiao N, Zhang Y. 2018. Coupled carbon, sulfur, and nitrogen cycles mediated by microorganisms in the water column of a shallow-water hydrothermal ecosystem. Front Microbiol 9 :2718. doi:10.3389/fmicb.2018.02718 30555427
41 Wang L, Cheung MK, Kwan HS, Hwang J-S, Wong CK. 2015. Microbial diversity in shallow-water hydrothermal sediments of Kueishan Island, Taiwan as revealed by pyrosequencing. J Basic Microbiol 55 :1308–1318. doi:10.1002/jobm.201400811 26132902
42 Wang L, Cheung MK, Liu R, Wong CK, Kwan HS, Hwang J-S. 2017. Diversity of total bacterial communities and chemoautotrophic populations in sulfur-rich sediments of shallow-water hydrothermal vents off Kueishan Island, Taiwan. Microb Ecol 73 :571–582. doi:10.1007/s00248-016-0898-2 27909749
43 Zhang Y, Zhao Z, Chen C-TA, Tang K, Su J, Jiao N. 2012. Sulfur metabolizing microbes dominate microbial communities in andesite-hosted shallow-sea hydrothermal systems. PLoS ONE 7 :e44593. doi:10.1371/journal.pone.0044593 22970260
44 Sako Y, Takai K, Ishida Y, Uchida A, Katayama Y. 1996. Rhodothermus obamensis sp. nov., a modern lineage of extremely thermophilic marine bacteria. Int J Syst Evol Microbiol 46 :1099–1104.
45 Balch WE, Fox GE, Magrum LJ, Woese CR, Wolfe RS. 1979. Methanogens: reevaluation of a unique biological group. Microbiol Rev 43 :260–296. doi:10.1128/mr.43.2.260-296.1979 390357
46 Sasser M. 1990. Identification of bacteria by gas chromatography of cellular fatty acids. USFCC Newsl 20 :1–6.
47 Minnikin DE, O’Donnell AG, Goodfellow M, Alderson G, Athalye M, Schaal A, Parlett JH. 1984. An integrated procedure for the extraction of bacterial isoprenoid quinones and polar lipids. J Microbiol Methods 2 :233–241. doi:10.1016/0167-7012(84)90018-6
48 Tindall BJ. 1990. Lipid composition of Halobacterium lacusprofundi. FEMS Microbiology Letters 66 :199–202. doi:10.1016/0378-1097(90)90282-U
49 Lane DJ. 1991. 16S/23S rRNA sequencing, p 115–175. In Stackebrandt E, Goodfellow M (ed), Nucleic acid techniques in bacterial systematics. Wiley, Chichester.
50 Tamura K, Stecher G, Kumar S. 2021. MEGA11: molecular evolutionary genetics analysis version 11. Mol Biol Evol 38 :3022–3027. doi:10.1093/molbev/msab120 33892491
51 Koren S, Walenz BP, Berlin K, Miller JR, Bergman NH, Phillippy AM. 2017. Canu: scalable and accurate long-read assembly via adaptivek-mer weighting and repeat separation. Genome Res 27 :722–736. doi:10.1101/gr.215087.116 28298431
52 Delcher AL, Bratke KA, Powers EC, Salzberg SL. 2007. Identifying bacterial genes and endosymbiont DNA with Glimmer. Bioinformatics 23 :673–679. doi:10.1093/bioinformatics/btm009 17237039
53 Borodovsky M, McIninch J. 1993. GENMARK: parallel gene recognition for both DNA strands. Comput Chem 17 :123–133. doi:10.1016/0097-8485(93)85004-V
54 Segata N, Börnigen D, Morgan XC, Huttenhower C. 2013. PhyloPhlAn is a new method for improved phylogenetic and taxonomic placement of microbes. Nat Commun 4 :2304. doi:10.1038/ncomms3304 23942190
55 Konstantinidis KT, Tiedje JM. 2005. Genomic insights that advance the species definition for prokaryotes. Proc Natl Acad Sci U S A 102 :2567–2572. doi:10.1073/pnas.0409727102 15701695
56 Pritchard L, Glover RH, Humphris S, Elphinstone JG, Toth IK. 2016. Genomics and taxonomy in diagnostics for food security: soft-rotting enterobacterial plant pathogens. Anal. Methods 8 :12–24. doi:10.1039/C5AY02550H
57 Qin QL, Xie BB, Zhang XY, Chen XL, Zhou BC, Zhou J, Oren A, Zhang YZ. 2014. A proposed genus boundary for the prokaryotes based on genomic insights. J Bacteriol 196 :2210–2215. doi:10.1128/JB.01688-14 24706738
58 Peters JW, Schut GJ, Boyd ES, Mulder DW, Shepard EM, Broderick JB, King PW, Adams MWW. 2015. [FeFe]- and [NiFe]-hydrogenase diversity, mechanism, and maturation. Biochim Biophys Acta 1853 :1350–1369. doi:10.1016/j.bbamcr.2014.11.021 25461840
59 Greening C, Biswas A, Carere CR, Jackson CJ, Taylor MC, Stott MB, Cook GM, Morales SE. 2016. Genomic and metagenomic surveys of hydrogenase distribution indicate H2 is a widely utilised energy source for microbial growth and survival. ISME J 10 :761–777. doi:10.1038/ismej.2015.153 26405831
60 Edgar RC. 2004. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32 :1792–1797. doi:10.1093/nar/gkh340 15034147
61 Eren AM, Kiefl E, Shaiber A, Veseli I, Miller SE, Schechter MS, Fink I, Pan JN, Yousef M, Fogarty EC, et al. . 2021. Community-led, integrated, reproducible multi-omics with anvi’o. Nat Microbiol 6 :3–6. doi:10.1038/s41564-020-00834-3 33349678
62 Chaudhari NM, Gupta VK, Dutta C. 2016. BPGA- an ultra-fast pan-genome analysis pipeline. Sci Rep 6 :24373. doi:10.1038/srep24373 27071527
63 Jain C, Rodriguez-R LM, Phillippy AM, Konstantinidis KT, Aluru S. 2018. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat Commun 9 :5114. doi:10.1038/s41467-018-07641-9 30504855
64 Barco RA, Garrity GM, Scott JJ, Amend JP, Nealson KH, Emerson D. 2020. A genus definition for Bacteria and Archaea based on a standard genome relatedness index. mBio 11 :e02475.31937639
65 Konstantinidis KT, Rosselló-Móra R, Amann R. 2017. Uncultivated microbes in need of their own taxonomy. ISME J 11 :2399–2406. doi:10.1038/ismej.2017.113 28731467
66 Luo C, Rodriguez-R LM, Konstantinidis KT. 2014. MyTaxa: an advanced taxonomic classifier for genomic and metagenomic sequences. Nucleic Acids Res 42 :e73–e73. doi:10.1093/nar/gku169 24589583
67 Pitcher RS, Brittain T, Watmough NJ. 2002. Cytochrome cbb(3) oxidase and bacterial microaerobic metabolism. Biochem Soc Trans 30 :653–658. doi:10.1042/bst0300653 12196157
68 D’mello R, Hill S, Poole RK. 1996. The cytochrome bd quinol oxidase in Escherichia coli has an extremely high oxygen affinity and two oxygen-binding haems: implications for regulation of activity in vivo by oxygen inhibition. Microbiology (Reading) 142 ( Pt 4) :755–763. doi:10.1099/00221287-142-4-755 8936304
69 Bosma G, Braster M, Stouthamer AH, van Verseveld HW. 1987. Isolation and characterization of ubiquinol oxidase complexes from Paracoccus denitrificans cells cultured under various limiting growth conditions in the chemostat. Eur J Biochem 165 :657–663. doi:10.1111/j.1432-1033.1987.tb11491.x 3036512
70 Gabel C, Maier RJ. 1993. Oxygen-dependent transcriptional regulation of cytochrome aa3 in Bradyrhizobium japonicum. J Bacteriol 175 :128–132. doi:10.1128/jb.175.1.128-132.1993 8380149
71 Arai H, Roh JH, Kaplan S. 2008. Transcriptome dynamics during the transition from anaerobic photosynthesis to aerobic respiration in Rhodobacter sphaeroides 2.4.1. J Bacteriol 190 :286–299. doi:10.1128/JB.01375-07 17965166
72 Jugder B-E, Welch J, Aguey-Zinsou K-F, Marquis CP. 2013. Fundamentals and electrochemical applications of [Ni–Fe]-uptake hydrogenases. RSC Adv 3 :8142. doi:10.1039/c3ra22668a
73 Campbell BJ, Jeanthon C, Kostka JE, Luther GW, Cary SC. 2001. Growth and phylogenetic properties of novel bacteria belonging to the epsilon subdivision of the Proteobacteria enriched from Alvinella pompejana and deep-sea hydrothermal vents. Appl Environ Microbiol 67 :4566–4572. doi:10.1128/AEM.67.10.4566-4572.2001 11571157
74 Goris T, Diekert G. 2016. The genus Sulfurospirillum, p 209–234. In Adrian L, Löffler FE (ed), Organohalide-respiring bacteria. Springer, Berlin, Heidelberg.
75 Kodama Y, Ha LT, Watanabe K. 2007. Sulfurospirillum cavolei sp nov, a facultatively anaerobic sulfur-reducing bacterium isolated from an underground crude oil storage cavity. Int J Syst Evol Microbiol 57 :827–831. doi:10.1099/ijs.0.64823-0 17392214
76 Nakagawa S, Takai K. 2008. Deep-sea vent chemoautotrophs: diversity, biochemistry and ecological significance. FEMS Microbiol Ecol 65 :1–14. doi:10.1111/j.1574-6941.2008.00502.x 18503548
77 Wolfe RS, Penning N. 1977. Reduction of sulfur by spirillum 5175 and syntrophism with Chlorobium. Appl Environ Microbiol 33 :427–433. doi:10.1128/aem.33.2.427-433.1977 848960
78 Schumacher W, Kroneck PMH. 1992. Anaerobic energy metabolism of the sulfur-reducing bacterium “Spirillum” 5175 during dissimilatory nitrate reduction to ammonia. Arch Microbiol 157 :464–470. doi:10.1007/BF00249106
79 Scholz-Muramatsu H, Neumann A, Messmer M, Moore E, Diekert G. 1995. Isolation and characterization of Dehalospirillum multivorans gen. nov., sp. Nov., a tetrachloroethene-utilizaing, strictly anaerobic bacterium. Arch Microbiol 163 :48–56. doi:10.1007/BF00262203
80 Luijten M, de Weert J, Smidt H, Boschker HTS, de Vos WM, Schraa G, Stams AJM. 2003. Description of Sulfurospirillum halorespirans sp. nov., an anaerobic, tetrachloroethene-respiring bacterium, and transfer of Dehalospirillum multivorans to the genus Sulfurospirillum as Sulfurospirillum multivorans comb. nov. Int J Syst Evol Microbiol 53 :787–793. doi:10.1099/ijs.0.02417-0 12807201
81 Finster K, Liesack W, Tindall BJ. 1997. Sulfurospirillum arcachonense sp. nov., a new-microaerophilic sulfur-reducing bacterium. Int J Syst Evol Microbiol 47 :1212–1217.
82 Stolz JF, Ellis DJ, Blum JS, Ahmann D, Lovley DR, Oremland RS. 1999. Note: Sulfurospirillum barnesii sp. nov. and Sulfurospirillum arsenophilum sp. nov., new members of the Sulfurospirillum clade of the ε-Proteobacteria. Int J Syst Evol Microbiol 49 :1177–1180. doi:10.1099/00207713-49-3-1177
83 Sorokin DY, Tourova TP, Muyzer G. 2013. Isolation and characterization of two novel alkalitolerant sulfidogens from a Thiopaq bioreactor, Desulfonatronum alkalitolerans sp. nov., and Sulfurospirillum alkalitolerans sp. nov. Extremophiles 17 :535–543. doi:10.1007/s00792-013-0538-4 23564266
84 Jensen A, Finster K. 2005. Isolation and characterization of Sulfurospirillum carboxydovorans sp nov., a new microaerophilic carbon monoxide oxidizing epsilon Proteobacterium. Antonie Van Leeuwenhoek 87 :339–353. doi:10.1007/s10482-004-6839-y 15928986
85 Vignais PM, Billoud B. 2007. Occurrence, classification, and biological function of hydrogenases: an overview. Chem Rev 107 :4206–4272. doi:10.1021/cr050196r 17927159
86 Dross F, Geisler V, Lenger R, Theis F, Krafft T, Fahrenholz F, Kojro E, Duchêne A, Tripier D, Juvenal K. 1992. The quinone-reactive Ni/Fe-hydrogenase of Wolinella succinogenes. Eur J Biochem 206 :93–102. doi:10.1111/j.1432-1033.1992.tb16905.x 1587288
87 Brugna-Guiral M, Tron P, Nitschke W, Stetter K-O, Burlat B, Guigliarelli B, Bruschi M, Giudici-Orticoni MT. 2003. [NiFe] hydrogenases from the hyperthermophilic bacterium Aquifex aeolicus: properties, function, and phylogenetics. Extremophiles 7 :145–157. doi:10.1007/s00792-002-0306-3 12664267
88 Amend JP, Shock EL. 2001. Energetics of overall metabolic reactions of thermophilic and hyperthermophilic Archaea and bacteria. FEMS Microbiol Rev 25 :175–243. doi:10.1111/j.1574-6976.2001.tb00576.x 11250035
89 Adam N, Perner M. 2018. Microbially Mediated hydrogen cycling in deep-sea hydrothermal vents. Front Microbiol 9 :2873. doi:10.3389/fmicb.2018.02873 30532749
90 Heijnen JJ, Van Dijken JP. 1992. In search of a thermodynamic description of biomass yields for the chemotrophic growth of microorganisms. Biotechnol Bioeng 39 :833–858. doi:10.1002/bit.260390806 18601018
91 Goris T, Schubert T, Gadkari J, Wubet T, Tarkka M, Buscot F, Adrian L, Diekert G. 2014. Insights into organohalide respiration and the versatile catabolism of Sulfurospirillum multivorans gained from comparative genomics and physiological studies. Environ Microbiol 16 :3562–3580. doi:10.1111/1462-2920.12589 25186071
92 Buttet GF, Murray AM, Goris T, Burion M, Jin B, Rolle M, Holliger C, Maillard J. 2018. Coexistence of two distinct Sulfurospirillum populations respiring tetrachloroethene-genomic and kinetic considerations. FEMS Microbiol Ecol 94 . doi:10.1093/femsec/fiy018
93 Kruse S, Goris T, Westermann M, Adrian L, Diekert G. 2018. Hydrogen production by Sulfurospirillum species enables syntrophic interactions of Epsilonproteobacteria. Nat Commun 9 :4872. doi:10.1038/s41467-018-07342-3 30451902
94 Macalady JL, Dattagupta S, Schaperdoth I, Jones DS, Druschel GK, Eastman D. 2008. Niche differentiation among sulfur-oxidizing bacterial populations in cave waters. ISME J 2 :590–601. doi:10.1038/ismej.2008.25 18356823
95 Meier DV, Pjevac P, Bach W, Hourdez S, Girguis PR, Vidoudez C, Amann R, Meyerdierks A. 2017. Niche partitioning of diverse sulfur-oxidizing bacteria at hydrothermal vents. ISME J 11 :1545–1558. doi:10.1038/ismej.2017.37 28375213
96 McNichol J, Dyksma S, Mußmann M, Seewald JS, Sylva SP, Sievert SM. 2022. Genus-specific carbon fixation activity measurements reveal distinct responses to oxygen among hydrothermal vent Campylobacteria. Appl Environ Microbiol 88 :e0208321. doi:10.1128/AEM.02083-21 34788061
97 Borisov VB, Gennis RB, Hemp J, Verkhovsky MI. 2011. The cytochrome bd respiratory oxygen reductases. Biochim Biophys Acta 1807 :1398–1413. doi:10.1016/j.bbabio.2011.06.016 21756872
98 Ekici S, Pawlik G, Lohmeyer E, Koch HG, Daldal F. 2012. Biogenesis of cbb3-type cytochrome c oxidase in Rhodobacter capsulatus. Biochim Biophys Acta 1817 :898–910. doi:10.1016/j.bbabio.2011.10.011 22079199
99 Dahms H-U, Schizas NV, James RA, Wang L, Hwang J-S. 2018. Marine hydrothermal vents as templates for global change scenarios. Hydrobiologia 818 :1–10. doi:10.1007/s10750-018-3598-8
100 Chen CTA, Zeng Z, Kuo FW, Yang TF, Wang BJ, Tu YY. 2005. Tide-influenced acidic hydrothermal system offshore NE Taiwan. Chem Geol 224 :69–81. doi:10.1016/j.chemgeo.2005.07.022
101 Nakagawa S, Takaki Y, Shimamura S, Reysenbach AL, Takai K, Horikoshi K. 2007. Deep-sea vent epsilon-proteobacterial genomes provide insights into emergence of pathogens. Proc Natl Acad Sci U S A 104 :12146–12150. doi:10.1073/pnas.0700687104 17615243
