Original Research ARTICLE
Coupled Carbon, Sulfur, and Nitrogen Cycles Mediated by Microorganisms in the Water Column of a Shallow-Water Hydrothermal Ecosystem
- 1State Key Laboratory of Marine Environmental Science, Xiamen University, Xiamen, China
- 2College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
- 3Department of Limnology and Bio-Oceanography, University of Vienna, Vienna, Austria
- 4Department of Oceanography, National Sun Yat-sen University, Kaohsiung, Taiwan
Shallow-water hydrothermal vent ecosystems are distinctly different from deep-sea vents, as other than geothermal, sunlight is one of their primary sources of energy, so their resulting microbial communities differ to some extent. Yet compared with deep-sea systems, less is known about the active microbial community in shallow-water ecosystems. Thus, we studied the community compositions, their metabolic pathways, and possible coupling of microbially driven biogeochemical cycles in a shallow-water hydrothermal vent system off Kueishantao Islet, Taiwan, using high-throughput 16S rRNA sequences and metatranscriptome analyses. Gammaproteobacteria and Epsilonbacteraeota were the major active bacterial groups in the 16S rRNA libraries and the metatranscriptomes, and involved in the carbon, sulfur, and nitrogen metabolic pathways. As core players, Thiomicrospira, Thiomicrorhabdus, Thiothrix, Sulfurovum, and Arcobacter derived energy from the oxidation of reduced sulfur compounds and fixed dissolved inorganic carbon (DIC) by the Calvin-Benson-Bassham (CBB) or reverse tricarboxylic acid cycles. Sox-dependent and reverse sulfate reduction were the main pathways of energy generation, and probably coupled to denitrification by providing electrons to nitrate and nitrite. Sulfur-reducing Nautiliaceae members, accounting for a small proportion in the community, obtained energy by the oxidation of hydrogen, which also supplies metabolic energy for some sulfur-oxidizing bacteria. In addition, ammonia and nitrite oxidation is another type of energy generation in this hydrothermal system, with marker gene sequences belonging to Thaumarchaeota/Crenarchaeota and Nitrospina, respectively, and ammonia and nitrite oxidation was likely coupled to denitrification by providing substrate for nitrate and nitrite reduction to nitric oxide. Moreover, unlike the deep-sea systems, cyanobacteria may also actively participate in major metabolic pathways. This study helps us to better understand biogeochemical processes mediated by microorganisms and possible coupling of the carbon, sulfur, and nitrogen cycles in these unique ecosystems.
The discovery of marine hydrothermal vents greatly enhanced our understanding of microbial habitats and survival strategies as well as the origins of life. Microbial communities in deep-sea hydrothermal systems have been intensively studied (Brazelton and Baross, 2010; Xie et al., 2011; Grosche et al., 2015; Anantharaman et al., 2016) since the discovery of these vents in 1977. Most microbes in deep-sea hydrothermal vent ecosystems carry out chemosynthesis, which fixes carbon dioxide (CO2) into organic compounds using the energy released by chemical reactions; it does not require sunlight. However, in shallow-water hydrothermal vent ecosystems, generally at water depths less than 200 m, chemolithoautotrophy and photoautotrophy occur simultaneously (Maugeri et al., 2009; Zhang et al., 2012; Gomez-Saez et al., 2017; Tang et al., 2018). Previous surveys of bacterial 16S rRNA genes using tag pyrosequencing and clone libraries revealed a high abundance of chemoautotrophs within the classes Gammaproteobacteria and Epsilonproteobacteria (reclassified to a new phylum Epsilonbacteraeota; Waite et al., 2017) in shallow-water hydrothermal systems (Tang et al., 2013; Gomez-Saez et al., 2017). In addition, Cyanobacteria were also frequently found (Zhang et al., 2012; Tang et al., 2013; Gomez-Saez et al., 2017). Despite nearly 30 published studies on shallow-water hydrothermal systems, many open questions remain about the chemosynthetic and photosynthetic microbes, including the metabolic pathways they use, how the pathways are coupled with each other, and what factors control their ecology.
Shallow (water depth < 30 m) submarine hydrothermal activity has been observed within 1 km east of Kueishantao Islet, off Taiwan. This hydrothermal system has unique geochemical characteristics and is driven by both sunlight and geothermal energy; thus, it is an ideal ecosystem to study coupled metabolic pathways and microbially driven biogeochemical cycles in extreme environments. Gas emitted from the Kueishantao hydrothermal vents are composed of CO2, nitrogen (N2), methane (CH4) and small amounts of hydrogen sulfide (H2S) (Chen et al., 2005; Chen et al., 2016). The hydrothermal fluids originate with deep magmatic matter and meteoric water from the Kueishantao Islet (Liu et al., 2010), and mix with seawater to form the final hydrothermal fluids. Fractures are widely developed around the andesite-hosted hydrothermal vent and therefore relatively oxygen-rich seawater seeps through these fractures in the seafloor. A previous study indicated that H2S in the Kueishantao hydrothermal system mainly originates from thermal reductive reactions of seawater and sulfate radicals, suggesting that seawater is the initial source of H2S (Zhang, 2013). Thus, steep geochemical gradients form when reduced hydrothermal fluids meet the oxidized seawater. Electron donors in the gradients include sulfur (S0), thiosulfate (S2O32-), hydrogen (H2), organics, formate and fumarate, while nitrate, oxygen (O2), S0, and S2O32- are the major identified electron acceptors (Xie et al., 2011; Anantharaman et al., 2016; Tang et al., 2018). Thus, a series of redox reactions occur and drive the carbon, nitrogen, and sulfur cycles in this hydrothermal ecosystem, which consists of the vent fluids and the water surrounding the vent.
In this study, high-throughput 16S rRNA sequencing and metatranscriptome analyses were carried out to investigate the microbial community in the surface water immediately above a white hydrothermal vent and the bottom water next to the vent (Supplementary Figure S1). The potentially metabolically active bacterial compositions and metabolic pathways in the hydrothermal ecosystem were determined to improve our understanding of biogeochemical processes mediated by microorganisms and coupling of the carbon, sulfur, and nitrogen cycles in the water column of this unique ecosystem, driven by both sunlight and geothermal energy.
Materials and Methods
Study Sites and Sampling
A cluster of shallow hydrothermal vents is located within 1 km east of Kueishantao Islet (Supplementary Figure S1). A white vent was identified by scuba divers and its position was located using a Global Positioning System (24.83N, 121.96E). Two samples were collected in April 2014 from the surface water immediately above the vent (SW) and bottom water next to the vent (BW). All necessary permits were obtained for the described field studies, including the permits from the Coast Guard Administration of Taiwan and the Fisheries Management Office of the Yilan County.
On board, approximately 15 L of water per sample was pre-filtered through 3 μm pore-size polycarbonate membranes (EMD Millipore Corp., Darmstadt, Germany) and then collected in 0.22 μm Sterivex filter units (EMD Millipore Corp., Darmstadt, Germany) at a pressure of <0.03 MPa. The filtration was finished within 30 min to limit RNA degradation. As we could not assess whether the in situ expression profile had changed, a shorter filtering time (e.g., 10-15 min) may have been more appropriate. Each Sterivex was filled with RNAlater RNA stabilization solution (Ambion, United States), flash frozen and stored in liquid nitrogen until RNA extraction.
Salinity was obtained by the conversion of conductivity measurements from a Guildline salinometer (Autosal 8400B, Canada), and in-situ temperatures were determined by scuba divers using a thermocouple. pH values were measured with a pH meter (Radiometer PHM-85, Denmark) at 25°C and total alkalinity (TA) was measured with an alkalinity titrator. Nitrate, nitrite, and silicate were measured using a flow injection analyzer; the pink azo dye method was employed for nitrate and nitrite, and the silicomolybdenum blue method was used for silicate (Parsons, 1984; Pai et al., 1990). Dissolved inorganic carbon (DIC) was measured using a dissolved inorganic analyzer (AS-C3, Apollo SciTech, United States) with a precision of 0.1%.
RNA Extraction, PCR Amplification, and Sequencing
Total RNA was extracted using TRIzol reagent (Invitrogen, United States) according to Simms et al. (1993) with minor modifications, and treated with DNase (Qiagen, Valencia, CA, United States). DNA contamination was checked by amplifying the bacterial 16S rRNA genes with the universal primers 27F and 1492R. The purified RNA without DNA contamination was reverse transcribed to synthesize complementary DNA (cDNA) using the SuperScriptTM III First Strand Synthesis System kit with random primers (Invitrogen, United States) following the manufacturer’s specifications. The V3-V4 hypervariable region of bacterial 16S rRNA genes was amplified with primers 343F and 798R (Wilson et al., 1990; Rochelle et al., 1995). The forward primer 343F (bold letters) contained transposase sequences (underlined letters): TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGTACGGRAGGCAGCAG. The reverse primer 798R contained transposase sequences: GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGAGGGTATCTAATCCT. PCR reactions were conducted in triplicate. The 25 μL reactions contained 12.5 μL of 2 × KAPA HiFi Hotstart ReadyMix (Kapa Biosystems), 0.2 μM of each primer, and 2.5 μL of template cDNA (5 ng μL-1). The PCR thermal regime was initial denaturation at 95°C for 3 min, followed by 25 cycles of 30 s each of denaturation at 95°C, 30 s of annealing at 60°C, 30 s of extension at 72°C, and a final 10 min extension at 72°C. The products were purified with 0.8 volume of AMPure XP beads (Beckman Coulter, United States), and amplified with the P5 and P7 indexing primers. The products were further purified with 1 volume of AMPure XP beads, and then quantified with a Qubit dsDNA HS assay kit (Invitrogen, United States). Equimolar amounts of the PCR amplicons were pooled. Sequencing was conducted on an Illumina MiSeq sequencing system with MiSeq Reagent Kit v3 (Illumina, San Diego, CA, United States).
Raw data were filtered with MOTHUR v.1.32.1 (Schloss et al., 2009) for quality control. Paired-end DNA sequencing reads were chosen through the MOTHUR’s make.contigs routine and then assembled according to the overlap with no more than 1 base pair (bp) mismatch. Reads that either (1) contained more than one ambiguous nucleotide, (2) were shorter than 200 bp, or (3) had mononucleotide repeats greater than 8 bp were removed by screen.seqs routine. Sequences were clustered into operational taxonomic units with 97% similarity and the greengene database was used as a reference for taxonomic classification (Desantis et al., 2006). Based on the operational taxonomic unit assignment, diversity estimates (ACE, Chao 1, Shannon) were calculated using MOTHUR’s summary.single routine.
mRNA Purification, Sequencing, Assembly, and Annotation
The Ribo-Zero Magnetic Gold kit (Epicenter, Madison, WI, United States) was used to remove rRNA from the purified RNA without DNA contamination and obtain mRNA. Fragmentation buffer (Ambion, United States) was added to digest the mRNA to produce short fragments. Using these short fragments as templates, the first-strand cDNA was synthesized with random hexamer primers. The second-strand cDNA was synthesized using the first-strand cDNA as a template, then further purified with the QiaQuick PCR Purification kit (Qiagen, Valencia, CA, United States), and dissolved with elution buffer for end reparation and adding poly (A). After that, the short fragments were connected with sequencing adapters. Finally, the second-strand cDNA was degraded with an enzyme (UNG), and the product purified with a MiniElute PCR Purification kit (Qiagen, Valencia, CA, United States) and sequenced with an Illumina HiSeq2000 (Illumina, San Diego, CA, United States).
After removing low quality reads with 10% N, 20% low quality bases (<Q20) contaminated by adapter, or rRNA (Li et al., 2009), de novo assembly of the obtained clean reads was performed using Trinity (Grabherr et al., 2011). These sequences were further assembled using the sequence clustering software TGICL (Pertea et al., 2003) to obtain longer and non-redundant unigenes. For gene family clustering, the merged unigenes were classified into one cluster when the similarity between them was >70%. The merged unigenes were then blasted against public databases, including the National Center for Biotechnology Information non-redundant protein sequence database and Kyoto Encyclopedia of Genes and Genomes (blastx, E < 10-5). The protein with the highest sequence similarity was retrieved and annotated to each unigene. For annotated unigenes, protein coding sequence information was retrieved from the blast results. For unannotated unigenes, ESTscan was used to perform protein coding sequence predictions. The predicted protein coding sequences were then blasted (blastp, E < 10-5) against databases, including eggNOG (Huerta-Cepas et al., 2016), the Carbohydrate Active Enzyme database (Cantarel et al., 2009), and the Antibiotic Resistance Genes database (Liu and Pop, 2009), to obtain further functional information. To assess the quality of assembly and facilitate the subsequent relative abundance analysis, reads of each sample were mapped back to the merged unigenes using Bowtie2 (Langmead and Salzberg, 2012). Only unique reads (aligned concordantly exactly one time) mapped to each unigene were considered for relative abundance analysis, and the number of hits identified in the read-mapping step was normalized as a percentage to the total mapped reads of the individual transcriptomic dataset.
Results and Discussion
Physicochemical parameters were measured at three depths in the water column and from the hydrothermal vent fluids (Table 1 and Supplementary Figure S2). The water column pH values ranged from 6.24 to 7.05 because they were influenced by the acidic hydrothermal vent fluids (pH 5.34). Temperatures varied from 23.6°C to 24.1°C, and were much lower than the vent fluids (48°C). Salinity (34.19–34.52) and DIC concentrations (2028-2040 μmol L-1) were almost constant within the water column, but lower and higher, respectively, in the vent fluids (salinity: 33.76; DIC: 2352 μmol L-1). Total alkalinity ranged from 1637 to 2038 μmol L-1, with the minimum in the BW. The nitrate and nitrite concentrations varied from 0.14 to 2.86 μmol L-1 and 0.03 to 0.32 μmol L-1, respectively. The silicate concentrations decreased with depth, ranging from 5.3 to 15 μmol L-1. The phosphate concentrations were below detection limit. Overall, the shallow-water hydrothermal ecosystem was characterized by acidic vent fluids that were well mixed with seawater. Sulfide (S2-) concentrations were not measured in this study. Our previous studies suggested that compared with deep-sea vents, sulfide was relatively lower in this shallow vent system, varying from 0.01 to 0.85 mg L-1 (Zhang et al., 2012). Tang et al. (2013) also reported a similar range of sulfide concentrations from 0.11 mg L-1 in the water in the vent to 0.01 mg L-1 in the surface water.
Phylogenetic Identification of Active Bacteria Communities
A total of 65,925 and 67,591 qualified reads were rendered from the 16S rRNA libraries of the SW and BW, respectively (Supplementary Table S1). Estimates of bacterial community diversity indicated that there were no significant differences between SW and BW (Supplementary Table S2). Ribotypes of tags were identified phylogenetically and grouped by phylum, order, family, or genus. Gammaproteobacteria and Epsilonbacteraeota were the overwhelmingly dominant groups in the water column of the hydrothermal ecosystem (Figure 1A). The sulfide-oxidizing Thiomicrospira, which belongs to Thiotrichales, was the most abundant group within Gammaproteobacteria (Figure 1B) as well as the total libraries with 81.69% of the total tags in the BW library and 43.5% in the SW library. The second most abundant group in the total libraries was the family Helicobacteraceae (mostly unclassified) in Epsilonbacteraeota (Figure 1C), with 27.34% of the total tags in the SW library and 5.24% in the BW library. At the genus level, sulfur-oxidizing bacteria Sulfurimonas (Thiovulaceae), Arcobacter (Arcobacteraceae), and Sulfurovum (Sulfurovaceae) were the most abundant within Epsilonbacteraeota (Figure 1C). Some sulfur-oxidizing bacterial species of the genus Sulfurimonas were also nitrate-reducing bacteria that accept electrons from the oxidation of reduced inorganic sulfur compounds (Waite et al., 2017). These autotrophic denitrifiers have been frequently identified from diverse ecosystems, such as deep-sea hydrothermal vents and the central Baltic Sea (Takai et al., 2005; Brettar et al., 2006). Sulfur oxidation coupled with dissimilatory nitrate reduction is usually an important source of energy for DIC fixation in hydrothermal vents (Shao et al., 2010). All but one of the Sulfurimonas species can also use H2 as an energy source (Han and Perner, 2014, 2015). Autotrophic sulfur-oxidizing Gammaproteobacteria and Epsilonbacteraeota were the dominant microorganisms in this shallow-water hydrothermal ecosystem, and may significantly contribute to primary production utilizing reduced sulfur compounds as electron donors (Brazelton and Baross, 2010; Campbell et al., 2013).
FIGURE 1. Phylogenetic taxon distribution among the bacterial RNA-based libraries. (A) Relative abundance of bacterial phyla or classes in total tags of each library. (B) Relative abundance of bacterial orders or genera in total gammaproteobacterial tags. (C) Relative abundance of bacterial orders, families, or genera in total epsilonbacteraeota tags. (D) Relative abundance of bacterial orders or clades in total alphaproteobacterial tags.
The family Nautiliaceae was the second most abundant group within Epsilonbacteraeota with 3.55% of the total tags in the BW library and 2.05% in the SW library (Figure 1C), and contains the genera Lebetimonas and Nautilia as well as a large number of unclassified taxa. Members of this family are moderate thermophiles growing optimally between 40 and 60°C (Nakagawa and Takai, 2014). Under autotrophic conditions, Nautiliaceae members have an ability to grow anaerobically via respiratory S0 reduction with H2, utilizing H2 as an electron donor and S0 and other reduced compounds as electron acceptors (Campbell et al., 2006; Meyer and Huber, 2014; Nakagawa and Takai, 2014).
The Alphaproteobacteria were mainly composed of the SAR11 clade, which was more abundant in SW than BW (Figure 1D). The abundance of SAR11 in shallow hydrothermal vents may be the result of vent fluids mixing with seawater, since the SAR11 clade is widely distributed in surface seawater (Morris et al., 2002). It has been reported that the SAR11 clade does not possess genes mediating assimilatory sulfate reduction and thus they would require exogenous reduced sulfur for their metabolism (Tripp et al., 2008). So, high reduced sulfur conditions may also explain SAR11 abundance in the Kueishantao hydrothermal ecosystem. Rhodobacteraceae, which belongs to Rhodobacterales, was the second most abundant Alphaproteobacteria in BW (Figure 1D). Rhodobacteraceae comprises mainly aerobic photoheterotrophs and chemoheterotrophs, as well as purple non-sulfur bacteria that perform photosynthesis in anaerobic conditions; they are deeply involved in sulfur and carbon biogeochemical cycling (Pujalte et al., 2014).
Cyanobacteria, assigned to Synechococcus, had a relative abundance of 0.12%-1.7% in the SW and BW 16S rRNA libraries (Figure 1A), which could be the result of vent fluids mixing with seawater. Cyanobacteria carry out oxygenic photosynthesis in the surface water and may switch to anoxygenic photosynthesis at vents where H2S is high (Cohen et al., 1975, 1986; Padan, 1979). Further studies are needed to verify this switch. These results indicate the co-occurrence of chemoautotrophs and photoautotrophs/heterotrophs in a shallow-water hydrothermal vent, which is distinctly different from deep-sea hydrothermal ecosystems (Tarasov et al., 2005).
Overall, the dominant active microorganisms in the shallow hydrothermal ecosystem were chemoautotrophic bacteria that mostly have the potential to perform sulfur oxidation and reduction, followed by phototrophs. Therefore, the ecosystem is driven by energy derived mainly from sulfur redox reactions and light. Hydrogen might also be an important energy source for this ecosystem, as sulfur-oxidizing Sulfurimonas species and sulfur-reducing Nautiliaceae species utilize H2 as an electron donor and they were found in the active assemblages (Muyzer et al., 1995; Brinkhoff et al., 1999). These results suggest multiple types of energy generation are performed in this system. Biogeochemical conditions to some extent determine the community composition and type of energy generation. Our previous study indicated that CH4 concentration was the statistically significant variable that explains the community structure in this shallow-water hydrothermal vent (Zhang et al., 2012). It is because distinctly different CH4 concentrations shaped different redox environments between the surface and bottom waters and consequently influenced the distribution of the community composition (Zhang et al., 2012). This is consistent with the results of this study where we found higher relative abundance of sulfur-reducing Nautiliaceae species in the BW library than the SW library. This result is also consistent with the finding of one order of magnitude higher sulfide concentration in the vent water than in the surface water (Tang et al., 2013).
Major Metabolic Activities and Pathways
In metatranscriptomes, 16,714 and 34,995 unigenes were obtained from SW and BW, respectively (Supplementary Table S3). Previous studies suggested that there are six pathways for CO2 fixation (Hügler and Sievert, 2011). In this study, all enzymes included in the Calvin-Benson-Bassham (CBB) and the reductive tricarboxylic acid (rTCA) cycles were detected (Figure 2). Overall, the relative transcript abundance of genes encoding enzymes in the rTCA cycle did not show a significant difference between the two samples (same order of magnitude). The relative transcript abundance of the key gene encoding ATP-citrate lyase (EC:220.127.116.11) in the rTCA cycle was nearly two times higher in BW. In contrast, another key gene (EC:18.104.22.168, 2-oxoglutarate synthase) had a slightly higher relative abundance in SW. However, the relative transcript abundance of genes encoding enzymes involved in the CBB cycle was one order of magnitude higher in BW, except for the genes encoding ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) (EC:22.214.171.124) and transketolase (EC:126.96.36.199), which were one order of magnitude lower in BW (Figure 2). Moreover, the relative transcript abundance of genes encoding RuBisCO was distinctly higher (one to six orders of magnitude) than all other genes involved in the rTCA and CBB cycles. The rTCA cycle is universal in hydrothermal vent environments; it is performed by chemoautotrophic Epsilonbacteraeota (Hügler et al., 2005), and is the most economical pathway of carbon fixation in bacteria (Mangiapia and Scott, 2016). This pathway is particularly advantageous in an energy-limited environment (Oulas et al., 2016). The CBB cycle was the first carbon fixation pathway to be discovered about 60 years ago (Wilson and Calvin, 1955), and autotrophic Gammaproteobacteria mainly perform carbon fixation through this pathway (Hügler and Sievert, 2011).
FIGURE 2. Carbon fixation pathways identified in the metatranscriptomes based on Kyoto Encyclopedia of Genes and Genomes pathway maps. Enzyme classification (EC) numbers for each step are included in boxes. The box color indicates relative abundance of reads in total mapped reads for each enzyme; red indicates higher relative abundance and blue lower relative abundance. The surface water (SW) is always shown on the left and bottom water (BW) on the right. Key enzymes are marked with an asterisk. A bold border indicates dual functional proteins. Gray represents absence. rTCA, reductive tricarboxylic acid cycle; CBB, Calvin-Benson-Bassham cycle.
In hydrothermal ecosystems, sulfur redox reactions mediated by chemoautotrophic bacteria are the main processes for energy production and reliant upon hydrothermal vent fluids as an energy source (McCollom, 2000). In shallow-water hydrothermal vents, as metatranscriptome analyses revealed, oxidation of reduced sulfur compounds can play a key role in energy production. For example, gene encoding enzymes that catalyze sulfur compound oxidations from the most reduced sulfur species (–II) to the most oxidized inorganic form (+IV) were fully restored in the metatranscriptomes (Figure 3). The reverse sulfate reduction pathway includes the sulfide oxidation step involving flavocytochrome c sulfide dehydrogenase (FccAB), which initiates electron flow from sulfide to the transport chain. It also includes steps of sulfur oxidation to sulfite catalyzed by dissimilatory sulfite reductase (EC:188.8.131.52) and sulfite oxidation to sulfate catalyzed by sulfite dehydrogenase (EC:184.108.40.206). Additionally, sulfide/disulfide could be oxidized to sulfite catalyzed by assimilatory NADPH-sulfite reductase (EC:220.127.116.11) and ferredoxin-sulfite reductase (EC:18.104.22.168). The Sox pathway was another important sulfur oxidation pathway found in the Kueishantao white hydrothermal vent, as soxABCDXYZ genes encoding the thiosulfate oxidation multienzyme complex were retrieved from the metatranscriptomes (Figure 3). In addition, the sulfate adenylyltransferase (EC:22.214.171.124) and adenylylsulfate reductase (EC:126.96.36.199) were also identified in both BW and SW; these are involved in the energy-yielding dissimilatory sulfate reduction pathway (Figure 3). These results are consistent with the active bacterial community compositions, in which a large fraction of sulfur-oxidizing bacteria and a small number of sulfate-reducing bacteria were retrieved from the shallow hydrothermal vent. Overall, the relative transcript abundances of genes encoding the enzymes involved in these sulfur redox reactions were one to two orders of magnitude higher in BW than SW. The only exceptions were the soxB gene and the gene encoding sulfate adenylyltransferase, which were comparable between the two samples.
FIGURE 3. Sulfur metabolic pathways identified in the metatranscriptomes based on Kyoto Encyclopedia of Genes and Genomes pathway maps. Enzyme classification (EC) numbers for each step are included in boxes. The box color indicates relative abundance of reads in total mapped reads for each enzyme. The color scheme is the same as in Figure 2. Gray represents absence. SO42-, sulfate; SO32-, sulfite; S0, sulfur; S2O32-, thiosulfate; S2- sulfide; APS, adenylyl sulfate; FccAB, flavocytochrome c/sulfide dehydrogenase; SOX, thiosulfate oxidation multienzyme complex.
Denitrification, coupled to sulfur oxidation, usually plays an important role in most hydrothermal vent environments (Bourbonnais et al., 2012, 2014; Voss et al., 2013; Dang and Chen, 2017). Genes encoding dissimilatory nitrate reductase (EC:188.8.131.52, napAB and narGHI) and nitrite reductase (EC:184.108.40.206, nitric oxide-forming), involved in denitrification, were found in the metatranscriptomes from both SW and BW (Supplementary Table S4). The transcripts of genes encoding assimilatory nitrate reductase (EC:220.127.116.11, narB and nasA, and EC:18.104.22.168), NAD(P)H-nitrite reductase (EC:22.214.171.124), and ferredoxin-nitrite reductase (EC:126.96.36.199) were also retrieved from the metatranscriptomes (Figure 4). Nitrate and nitrite could be important electron acceptors when dissimilatory nitrate reduction (to ammonium, DNRA, or denitrification) is coupled to sulfur oxidation reactions (Nakagawa and Takai, 2008; Dang and Chen, 2017). For example, the sulfur and thiosulfate-oxidizing bacteria Sulfurovum sp., which were detected in the 16S rRNA libraries, have been reported to use nitrate or oxygen as electron acceptors (Inagaki et al., 2004; Giovannelli et al., 2016; Dang and Chen, 2017). In the present study, the transcripts of genes encoding nitrite reductase (ammonia-forming) related to DNRA (e.g., NirBD and NrfAH, according to the Kyoto Encyclopedia of Genes and Genomes) were not retrieved. Previous studies indicated that DNRA (nitrite ammonification) was mainly important in relatively reducing environments, as found in nutrient-rich coastal sediments (Rysgaard et al., 1996; Christensen et al., 2000; An and Gardner, 2002), while denitrification was more important in low to moderate organic-loading sediments (Fossing et al., 1995; Zopfi et al., 2001). For instance, a metagenomics analysis of sulfur-oxidizing Gammaproteobacteria from a coastal ecosystem in the eastern South Pacific suggested a coupling of sulfur oxidation and DNRA in oxygen-deficient waters (Murillo et al., 2014).
FIGURE 4. Nitrogen metabolic pathways identified in the metatranscriptomes based on Kyoto Encyclopedia of Genes and Genomes pathway maps. Enzyme classification (EC) numbers for each step are included in boxes. The box color indicates relative abundance of reads in total mapped reads for each enzyme. The color scheme is the same as in Figure 2. Gray represents absence.
As nitrification marker genes, the transcripts of ammonia monooxygenase (amoABC) genes (EC:188.8.131.52), responsible for ammonia oxidation to hydroxylamine, and nitrite oxidoreductase beta subunit (nxrB), responsible for nitrite oxidation to nitrate, were retrieved in the metatranscriptomes of SW and BW (Figure 4). Thus, there may be a coupling of nitrification and denitrification processes, resulting in nitrogen removal via the nitrite pathway.
In addition, the reduction of hydroxylamine and the hydrolysis of urea could produce ammonia in this shallow-water hydrothermal vent, as the genes encoding hydroxylamine oxidase (EC:184.108.40.206) and urease (EC:220.127.116.11) were identified in the metatranscriptomes (Figure 4). Ammonia can be stored in glutamine catalyzed by glutamine synthetase (EC:18.104.22.168) or glutamate catalyzed by glutamate dehydrogenase (EC:22.214.171.124 and EC:126.96.36.199). The former could be the main pathway, since the relative transcript abundance of the genes encoding glutamine synthetase was one to three orders of magnitude higher in SW and BW than glutamate dehydrogenase. Overall, the relative transcript abundances of genes encoding enzymes involved in the nitrogen cycle were up to one order of magnitude higher in BW. Detailed information on major carbon, sulfur, and nitrogen metabolic pathways are shown in (Supplementary Table S4).
Main Players in Major Metabolic Pathways
The sequences classified as Gammaproteobacteria, mainly Thiomicrospira, were the most abundant in the metatranscriptome dataset, which was consistent with the 16S rRNA libraries analysis that showed Thiomicrospira as having the highest abundance active population (Supplementary Table S5). Gene annotation and functional analysis indicated that the gammaproteobacterial genera Thiomicrospira, Thiomicrorhabdus, Hydrogenovibrio, and Thiothrix, Epsilonbacteraeota genera Lebetimonas and Caminibacter (Nautiliaceae), Sulfurovum and Nitratifractor (Sulfurovaceae), Sulfurimonas and Sulfuricurvum (Thiovulaceae), and Arcobacter (Arcobacteraceae), as well as Cyanobacteria (mainly Synechococcus and Prochlorococcus) and Archaea (Thaumarchaeota and Crenarchaeota) were the main autotrophic players in carbon fixation, and nitrogen and sulfur metabolism (Figure 5).
FIGURE 5. Distribution of relative transcript abundance of genes encoding enzymes included in carbon fixation, and nitrogen and sulfur metabolism among phylogenetic taxa. For the carbon fixation pathway, dark circles indicate reads of genes encoding key enzymes (marked with an asterisk in Figure 2); light circles indicate reads of genes encoding enzymes involved in each pathway. SW, surface water immediately above the vent; BW, bottom water next to the vent.
The most abundant transcript sequences belonged to Thiomicrospira (mainly T. crunogena XCL-2), most of which was involved in the CBB cycle (mainly RuBisCO sequences) (Figure 5). Thiomicrospira crunogena was originally isolated from the East Pacific Rise (Jannasch et al., 1985), and was subsequently detected in deep-sea hydrothermal vents (Wirsen et al., 1998) as well as shallow-water hydrothermal vents (Muyzer et al., 1995; Brinkhoff et al., 1999). It has been known as the representative ubiquitous chemolithoautotrophic sulfur-oxidizing bacteria (Scott et al., 2006) and has a remarkably high growth rate (Jannasch et al., 1985). In addition, the transcript sequences belonging to Thiothrix (mostly T. nivea DSM 5205), encoding RuBisCO, transketolase, and fructose-1,6-bisphosphatase in the CBB cycle, as well as FccAB and soxX, were retrieved from our metatranscriptome data. Thus, T. nivea could be chemolithoautotrophic via the CBB cycle or mixotrophic, since it has been reported to grow under heterotrophic conditions (Lapidus et al., 2011). A number of sequences classified as Thiomicrorhabdus within Gammaproteobacteria encoded assimilatory NAD(P)H-nitrite reductase and glutamine synthase, transketolase, and fructose-1,6-bisphosphate aldolase in the CBB cycle, as well as sulfate adenylyltransferase, FccAB, and soxABY.
The transcript sequences within Epsilonbacteraeota were involved in the rTCA pathway, dissimilatory/assimilatory nitrate reduction, and sulfur oxidizing/reducing processes (Figure 5). Nautiliaceae is a typical sulfur-reducing bacterial family living in hydrothermal vents, which can get energy from the oxidation of H2 or formate coupled with reduction of S0, to produce H2S (Hanson et al., 2013; Nakagawa and Takai, 2014). In our metatranscriptomes, they were abundant and involved in dissimilatory nitrate reduction, but very few were involved in sulfate reduction. For example, Lebetimonas had the most abundant transcript sequences in the total sequences involved in dissimilatory nitrate reduction (89.1% in SW and 54.1% in BW). All sequences belonging to Caminibacter were classified to C. mediatlanticus. Although they have been identified as sulfur or nitrate-reducing bacteria (Voordeckers et al., 2005), the sequences retrieved from our metatranscriptomes mainly related to the oxidation of sulfide to sulfur and the rTCA pathways. Sulfurovum (Sulfurovaceae) are typical (not strictly anaerobic) sulfur-oxidizing bacteria that can utilize nitrate or oxygen as electron acceptors (Inagaki et al., 2004). Nitratifractor, previously belonging to Nautiliaceae, were reclassified to Sulfurovaceae (Waite et al., 2017). In our metatranscriptomes, all sequences belonging to Nitratifractor were classified to N. salsuginis and mainly related to sulfur oxidation and rTCA pathways. Sulfurimonas denitrificans is a nitrate-reducing, sulfur-oxidizing species, but other Sulfurimonas species sequences retrieved from our metatranscriptomes were involved in assimilatory nitrate reduction, ammonia incorporation to glutamine, and glutamate pathways. The sequences assigned to Sulfuricurvum are associated with dissimilatory nitrate reduction, glutamine and cysteine synthesis, dissimilatory sulfate reduction, and rTCA pathways, although Sulfuricurvum species were reported to be sulfur-oxidizing bacteria (Kodama and Watanabe, 2004). Arcobacter had the most abundant transcript sequences in total sequences involved in sulfur metabolism, encoding sulfite dehydrogenase, FccAB, soxABCDXYZ, and cysteine synthase. Moreover, in our metatranscriptomes, hydrogenase sequences were identified in sulfur-reducing Lebetimonas and sulfur-oxidizing Sulfurimonas and Sulfuricurvum (Supplementary Table S6), suggesting that they might obtain energy through the oxidation of hydrogen.
Notably, Cyanobacteria (mainly Synechococcus and Prochlorococcus) participated in carbon (CBB cycle), nitrogen (ammonia assimilation), and sulfur (sulfide oxidation and sulfate assimilation) metabolism, although they were less abundant in the active population of the two 16S rRNA libraries (Figure 5 and Supplementary Table S5). Most Cyanobacteria are highly sensitive to sulfide toxicity (Oren et al., 1979), but some species are sulfide resistant (Cohen et al., 1986) or perform anoxygenic photosynthesis using sulfide rather than water as the terminal reductant (Frier et al., 1999; Kulp et al., 2008; Stal, 2012). A recent study indicated that under darkness and anoxygenic conditions, hydrogen in H2S accelerated the recovery of photosynthesis, and even enhanced photosynthetic rates at a given H2S concentration at low irradiance (Klatt et al., 2015). The geochemical properties of shallow-water hydrothermal environments possibly have retained many of the characteristics of the Earth’s early ocean (Baross and Hoffman, 1985) and thus Cyanobacteria in the habitat studied here may also preserve some of these characteristics. In addition, archaeal sequences were also retrieved, and belonged to Thaumarchaeota and Crenarchaeota, which were only involved in the nitrogen cycle, including amoABC and nitrite reductase (nirK). Sequences associated with the archaeal carbon fixation pathway, 3-hydroxypropionate/4-hydroxybutyrate (HP/HB) cycle, were not found (Figure 5).
Most of the reads within Alphaproteobacteria were associated with organic carbon metabolism and belonged to Rhodospirillales. However, the most abundant transcript sequences within Alphaproteobacteria involved in nitrogen and sulfur metabolism (sulfur oxidation, cysteine synthesis, and glutamine/glutamate synthesis) belonged to the order Rhodobacterales (Figure 5). Moreover, the transcript sequences encoding bacteriochlorophyll synthase were retrieved and assigned to Ahrensia, which is a typical aerobic photoheterotrophic genus in Rhodobacteraceae. The reads belonging to the SAR11 clade only accounted for ∼2% in SW and ∼7% in BW in our metatranscriptomes, despite the high proportions within Alphaproteobacteria in the two 16S rRNA libraries (∼73% in SW and ∼36% in BW); these transcript sequences are associated with glutamine and cysteine synthesis and dissimilatory sulfate reduction pathways. The most abundant transcript sequences within heterotrophic Gammaproteobacteria involved in nitrogen and sulfur metabolism belonged to Alteromonadales (Figure 5).
Most genes related to stress tolerance showed higher relative transcript abundance in BW, including genes encoding multiple molecular chaperones, anaerobic regulator proteins, and various antioxidant genes (Supplementary Figure S3). In addition, the relative transcript abundances of genes encoding DNA mismatch repair proteins were high in both SW and BW (Supplementary Figure S3). The DNA mismatch repair system plays a vital role in an organism’s response to DNA damage and maintains genomic stability. These results suggest that the microbial communities have evolved extensive DNA repair systems, such as heat shock stress response systems and regulators of anaerobiosis or antioxidant systems, to cope with the extreme conditions present at a hydrothermal vent (Kiley and Beinert, 1998; Xie et al., 2011).
High-throughput 16S rRNA sequences and metatranscriptome analyses revealed that Gammaproteobacteria and Epsilonbacteraeota were the most active bacterial populations involved in the major carbon, sulfur, and nitrogen metabolic pathways in a shallow-water hydrothermal ecosystem. The major sulfur oxidizers were Thiomicrospira, Thiomicrorhabdus, and Thiothrix from Gammaproteobacteria, and Arcobacter and Sulfurovum from Epsilonbacteraeota, which showed high transcript abundances in the genes involved in the SOX/reverse sulfate reduction pathway. The sulfur-reducing Nautiliaceae contributed very few transcripts to sulfate reduction, but showed a high level of transcription for genes involved in denitrification processes. In addition, Thiomicrorhabdus exhibited a range of genes related to assimilatory nitrate reduction.
We illustrated these major metabolic pathways and the possible coupling between microbially driven biogeochemical cycles in this ecosystem (Figure 6). Hydrogen sulfide contained in hydrothermal fluids from the Kueishantao vents is produced from thermal reduction of seawater sulfate radicals when seawater seeps through fractures in the seafloor (Figure 6). This suggested that seawater is the initial source of H2S and geothermal heat is the primary energy source. Consequently, chemolithoautotrophic microbes (mainly members within Gammaproteobacteria and Epsilonbacteraeota) derive energy from the oxidation of reduced sulfur compounds and fix DIC by the CBB and rTCA cycles. Sox-dependent and reverse sulfate reduction are the main pathways of energy generation, and are probably coupled to denitrification by the provision of electrons to nitrate and nitrite (Figure 6). Oxygen is also a possible electron acceptor for sulfur oxidation. In addition, hydrogen oxidation supplies metabolic energy for some sulfur-oxidizing (e.g., Sulfurimonas sp.) and sulfur-reducing (e.g., Lebetimonas sp.) bacteria, coupled to the reduction of nitrate and sulfur, respectively. Ammonia and nitrite oxidation are other types of energy generation carried out by Thaumarchaeota/Crenarchaeota and Nitrospina, respectively, in this hydrothermal system, and coupled to denitrification by providing nitrate and nitrite substrate (Figure 6). Furthermore, driven by light energy, Cyanobacteria and aerobic photoheterotrophs also actively participate in major metabolic pathways. We speculate that Cyanobacteria perform oxygenic photosynthesis, fixing CO2 through the CBB cycle and producing O2 in the surface water, and then switch to anoxygenic photosynthesis fixing CO2 and producing S0 or sulfite in the bottom water next to the vent (Figure 6). Our results indicate the co-occurrence of chemoautotrophs and photoautotrophs/heterotrophs in a shallow-water hydrothermal vent, which is distinctly different from deep-sea hydrothermal ecosystems. Overall, the oxidation of reduced sulfur compounds, using oxygen or nitrate as electron acceptors, provide significant energy for carbon fixation in this shallow-water hydrothermal vent ecosystem, which uses sunlight and geothermal as primary energy sources. This study helps us to better understand biogeochemical processes mediated by microorganisms and the possible coupling of the carbon, sulfur, and nitrogen cycles in this unique ecosystem.
FIGURE 6. Schematic diagram illustrating the coupling of the carbon, sulfur, and nitrogen cycles mediated by microorganisms in the shallow-water hydrothermal ecosystem. SO42-, sulfate; SO32-, sulfite; S/S0, sulfur; S2O32-, thiosulfate; S2-, sulfide; H2S, hydrogen sulfide; SO2, sulfur dioxide; NO3-, nitrate; NO2-, nitrite; NH4+, ammonium; NO, nitric oxide; H2, hydrogen; CH4, methane; CO2, carbon dioxide; O2, oxygen.
The 16S rRNA sequence and metatranscriptomic datasets were deposited in the Short Reads Archive (National Center for Biotechnology Information) under accession numbers SRP145248 and SRP145422.
YZ conceived and designed the research. YL, KT, LZ, ZZ, XX, C-TC, and DW conducted the experiments. YL, KT, YZ, LZ, ZZ, XX, C-TC, and NJ analyzed the data. YL, KT, LZ, and YZ wrote the paper. NJ, C-TC, and DW contributed to the interpretation of results and critical revision.
This work was funded by the National Key Research and Development Program (2016YFA0601400), NSFC projects (Joint Fund for the Promotion of Cross-Strait Science and Technology Cooperation “Coupling between microbially driven carbon, nitrogen, and sulfur cycle in andesite-hosted shallow-sea hydrothermal systems offshore NE Taiwan”, 41676125, 91751207, and 41721005), and the Fundamental Research Fund for the Central Universities (20720180123). This study was a contribution to the international IMBER project.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Bing-Jye Wang and Seawatch Co., Ltd. for their assistance in sampling. We also thank Kara Bogus, PhD, from Liwen Bianji, Edanz Editing China (www.liwenbianji.cn/ac), for editing the English text of a draft of this manuscript.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.02718/full#supplementary-material
An, S. M., and Gardner, W. S. (2002). Dissimilatory nitrate reduction to ammonium (DNRA) as a nitrogen link, versus denitrification as a sink in a shallow estuary (Laguna Madre/Baffin Bay, Texas). Mar. Ecol. Prog. Ser. 237, 41–50. doi: 10.3354/meps237041
Anantharaman, K., Breier, J. A., and Dick, G. J. (2016). Metagenomic resolution of microbial functions in deep-sea hydrothermal plumes across the Eastern Lau Spreading Center. ISME J. 10, 225–239. doi: 10.1038/ismej.2015.81
Baross, J. A., and Hoffman, S. E. (1985). Submarine hydrothermal vents and associated gradient environments as sites for the origin and evolution of life. Orig. Life 15, 327–345. doi: 10.1007/BF01808177
Bourbonnais, A., Juniper, S. K., Butterfield, D. A., Anderson, R. E., and Lehmann, M. F. (2014). Diversity and abundance of Bacteria and nirS-encoding denitrifiers associated with the Juan de Fuca Ridge hydrothermal system. Ann. Microbiol. 64, 1691–1705. doi: 10.1007/s13213-014-0813-3
Bourbonnais, A., Lehmann, M. F., Butterfield, D. A., and Juniper, S. K. (2012). Subseafloor nitrogen transformations in diffuse hydrothermal vent fluids of the Juan de Fuca Ridge evidenced by the isotopic composition of nitrate and ammonium. Geochem. Geophys. Geosyst. 13, 1–23. doi: 10.1029/2011gc003863
Brazelton, W. J., and Baross, J. A. (2010). Metagenomic comparison of two Thiomicrospira lineages inhabiting contrasting deep-sea hydrothermal environments. PLoS One 5:e13530. doi: 10.1371/journal.pone.0013530
Brettar, I., Labrenz, M., Flavier, S., Botel, J., Kuosa, H., Christen, R., et al. (2006). Identification of a Thiomicrospira denitrificans-like epsilonproteobacterium as a catalyst for autotrophic denitrification in the central Baltic Sea. Appl. Environ. Microbiol. 72, 1364–1372. doi: 10.1128/AEM.72.2.1364-1372.2006
Brinkhoff, T., Sievert, S. M., Kuever, J., and Muyzer, G. (1999). Distribution and diversity of sulfur-oxidizingThiomicrospira spp. at a shallow-water hydrothermal vent in the Aegean Sea (Milos, Greece). Appl. Environ. Microbiol. 65, 3843–3849.
Campbell, B. J., Polson, S. W., Lisa, Z. A., Williamson, S. J., Lee, C. K., Eric, W. K., et al. (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
Cantarel, B. L., Coutinho, P. M., Rancurel, C., Bernard, T., Lombard, V., and Henrissat, B. (2009). The Carbohydrate-active EnZymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res. 37, D233–D238. doi: 10.1093/nar/gkn663
Chen, C. T. A., Zeng, Z., Kuo, F. W., Yang, T. F., Wang, B. J., and Tu, Y. Y. (2005). Tide-influenced acidic hydrothermal system offshore NE Taiwan. Chem. Geol. 224, 69–81. doi: 10.1016/j.chemgeo.2005.07.022
Chen, X. G., Zhang, H. Y., Li, X., Chen, C. T. A., Yang, T. F., and Ye, Y. (2016). The chemical and isotopic compositions of gas discharge from shallow-water hydrothermal vents at Kueishantao, offshore northeast Taiwan. Geochem. J. 50, 341–355. doi: 10.2343/geochemj.2.0425
Christensen, P. B. S., Rysgaard, N. P., Sloth, T. D., and Schwærter, S. (2000). Sediment mineralization, nutrient fluxes, denitrification and dissimilatory nitrate reduction to ammonium in an estuarine fjord with sea cage trout farms. Aquat. Microb. Ecol. 21, 73–84. doi: 10.3354/ame021073
Cohen, Y., Jørgensen, B. B., Revsbech, N. P., and Poplawski, R. (1986). Adaptation to hydrogen sulfide of oxygenic and anoxygenic photosynthesis among Cyanobacteria. Appl. Environ. Microbiol. 51, 398–407.
Dang, H., and Chen, C. A. (2017). Ecological energetic perspectives on responses of nitrogen-transforming chemolithoautotrophic microbiota to changes in the marine environment. Front. Microbiol. 8:1246. doi: 10.3389/fmicb.2017.01246
Desantis, T. Z., Hugenholtz, P., Larsen, N., Rojas, M., Brodie, E. L., Keller, K., et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072. doi: 10.1128/AEM.03006-05
Fossing, H., Gallardo, V. A., Jørgensen, B. B., Hüttel, M., Nielsen, L. P., Schulz, H., et al. (1995). Concentration and transport of nitrate by the mat-forming sulphur bacterium Thioploca. Nature 374, 713–715. doi: 10.1038/374713a0
Frier, I., Rethmeier, J., and Fischer, U. (1999). “Molecular properties of soluble cytochrome c-552 and its participation in sulfur metabolism of Oscillatoria strain Bo32,” in The Phototrophic Prokaryotes, eds G. A. Peschek, W. Löffelhardt, and G. Schmetterer (Boston, MA: Springer).
Giovannelli, D., Chung, M., Staley, J., Starovoytov, V., Le Bris, N., and 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
Gomez-Saez, G. V., Pop Ristova, P., Sievert, S. M., Elvert, M., Hinrichs, K. U., and Buhring, S. I. (2017). Relative importance of chemoautotrophy for primary production in a light exposed marine shallow hydrothermal system. Front. Microbiol. 8:702. doi: 10.3389/fmicb.2017.00702
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883
Grosche, A., Sekaran, H., Perez-Rodriguez, I., Starovoytov, V., and 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
Hanson, T. E., Campbell, B. J., Kalis, K. M., Campbell, M. A., and Klotz, M. G. (2013). Nitrate ammonification by Nautilia profundicola AmH: experimental evidence consistent with a free hydroxylamine intermediate. Front. Microbiol. 4:180. doi: 10.3389/fmicb.2013.00180
Huerta-Cepas, J., Szklarczyk, D., Forslund, K., Cook, H., Heller, D., Walter, M. C., et al. (2016). eggNOG 4.5: a hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences. Nucleic Acids Res. 44, D286–D293. doi: 10.1093/nar/gkv1248
Hügler, M., Wirsen, C. O., Fuchs, G., Taylor, C. D., and Sievert, S. M. (2005). Evidence for autotrophic CO2 fixation via the reductive tricarboxylic acid cycle by members of the epsilon subdivision of proteobacteria. J. Bacteriol. 187, 3020–3027. doi: 10.1128/JB.187.9.3020-3027.2005
Inagaki, F., Takai, K., Nealson, K. H., and 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
Jannasch, H. W., Nelson, D. C., and Robertson, L. A. (1985). Thiomicrospira crunogena sp. nov. a colorless, sulfur-oxidizing bacterium from a deep-sea hydrothermal vent. J. Syst. Bacteriol. 35, 422–424. doi: 10.1099/00207713-35-4-422
Klatt, J. M., Haas, S., Yilmaz, P., Beer, D. D., and Polerecky, L. (2015). Hydrogen sulfide can inhibit and enhance oxygenic photosynthesis in a cyanobacterium from sulfidic springs. Environ. Microbiol. 17:3301. doi: 10.1111/1462-2920.12791
Kodama, Y., and Watanabe, K. (2004). Sulfuricurvum kujiense gen. nov., sp. nov., a facultatively anaerobic, chemolithoautotrophic, sulfur-oxidizing bacterium isolated from an underground crude-oil storage cavity. Int. J. Syst. Evol. Microbiol. 54, 2297–2300. doi: 10.1099/ijs.0.63243-0
Kulp, T. R., Hoeft, S. E., Asao, M., Madigan, M. T., Hollibaugh, J. T., Fisher, J. C., et al. (2008). Arsenic(iii) fuels anoxygenic photosynthesis in hot spring biofilms from mono lake, california. Science 321, 967–970. doi: 10.1126/science.1160799
Lapidus, A., Nolan, M., Lucas, S., Glavina Del Rio, T., Tice, H., Cheng, J. F., et al. (2011). Genome sequence of the filamentous, gliding Thiothrix nivea neotype strain (JP2(T)). Stand. Genomic Sci. 5, 398–406. doi: 10.4056/sigs.2344929
Li, R., Yu, C., Li, Y., Lam, T. W., Yiu, S. M., Kristiansen, K., et al. (2009). SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics 25, 1966–1967. doi: 10.1093/bioinformatics/btp336
Mangiapia, M., and Scott, K. (2016). From CO2 to cell: energetic expense of creating biomass using the Calvin–Benson–Bassham and reductive citric acid cycles based on genome data. FEMS Microbiol. Lett. 363;fnw054. doi: 10.1093/femsle/fnw054
Maugeri, T. L., Lentini, V., Gugliandolo, C., Italiano, F., Cousin, S., and Stackebrandt, E. (2009). Bacterial and archaeal populations at two shallow hydrothermal vents off Panarea Island (Eolian Islands. Italy). Extremophiles 13, 199–212. doi: 10.1007/s00792-008-0210-6
Morris, R. M., Rappé, M. S., Connon, S. A., Vergin, K. L., Siebold, W. A., and Carlson, C. A. (2002). SAR11 clade dominates ocean surface bacterioplankton communities. Nature 420, 806–810. doi: 10.1038/nature01240
Murillo, A. A., Ramirez-Flandes, S., DeLong, E. F., and Ulloa, O. (2014). Enhanced metabolic versatility of planktonic sulfur-oxidizing γ-proteobacteria in an oxygen-deficient coastal ecosystem. Front. Mar. Sci. 1:18. doi: 10.3389/fmars.2014.00018
Muyzer, G., Teske, A., Wirsen, C. O., and Jannasch, H. W. (1995). Phylogenetic relationships of Thiomicrospira species and their identification in deep-sea hydrothermal vent samples by denaturing gradient gel electrophoresis of 16S rDNA fragments. Arch. Microbiol. 164, 165–172. doi: 10.1007/BF02529967
Nakagawa, S., and Takai, K. (2014). “The family Nautiliaceae: the genera Caminibacter, Lebetimonas, and Nautilia,” in The Prokaryotes: Deltaproteobacteria and Epsilonproteobacteria, eds E. Rosenberg, E. F. DeLong, S. Lory, E. Stackebrandt, and F. Thompson (Berlin: Springer), 393–399.
Oren, A., Padan, E., and Malkin, S. (1979). Sulfide inhibition of Photosystem II in cyanobacteria (blue-green algae) and tobacco chloroplasts. Biochim. Biophys. Acta 546, 270–279. doi: 10.1016/0005-2728(79)90045-8
Oulas, A., Polymenakou, P. N., Seshadri, R., Tripp, H. J., Mandalakis, M., Paez-Espino, A. D., et al. (2016). Metagenomic investigation of the geologically unique Hellenic volcanic arc reveals a distinctive ecosystem with unexpected physiology. Environ. Microbiol. 18, 1122–1136. doi: 10.1111/1462-2920.13095
Pai, S. C., Yang, C. C., and Riley, J. P. (1990). Effects of acidity and molybdate concentration on the kinetics of the formation of the phosphoantimonylmolybdenum blue complex. Anal. Chim. Acta 299, 115–120. doi: 10.1016/S0003-2670(00)85116-8
Pertea, G., Huang, X., Liang, F., Antonescu, V., Sultana, R., Karamycheva, S., et al. (2003). TIGR Gene indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics 19, 651–652. doi: 10.1093/bioinformatics/btg034
Pujalte, M. J., Lucena, T., Ruvira, M. A., Arahal, D. R., and Macián, M. C. (2014). The Family Rhodobacteraceae. The Prokaryotes: Alphaproteobacteria and Betaproteobacteria. New York, NY: Springer. 439–512. doi: 10.1007/978-3-642-30197-1_377
Rochelle, P. A., Will, J. A. K., Fry, J. C., Jenkins, G. J. S., Parkes, R. J., Turley, C. M., et al. (1995). “Extraction and amplification of 16S rRNA genes from deep marine sediments and seawater to assess bacterial community diversity,” in Nucleic Acids in the Environment, eds J. D. van Elsas and J. T. Trevors (Berlin: Springer), 219–239.
Rysgaard, S., Risgaard-Petersen, N., and Sloth, N. P. (1996). Nitrification, denitrification, and nitrate ammonification in sediments of two coastal lagoons in southern France. Hydrobiologia 329, 133–141. doi: 10.1007/BF00034553
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Scott, K. M., Sievert, S. M., Abril, F. N., Ball, L. A., Barrett, C. J., Blake, R. A., et al. (2006). The genome of deep-sea vent chemolithoautotroph Thiomicrospira crunogena XCL-2. PLoS Biol. 4:e383. doi: 10.1371/journal.pbio.0040383
Shao, M. F., Zhang, T., and Fang, H. H. P. (2010). Sulfur-driven autotrophic denitrification: diversity, biochemistry, and engineering applications. Appl. Microbiol. Biotechnol. 88, 1027–1042. doi: 10.1007/s00253-010-2847-1
Takai, K., Campbell, B. J., Cary, S. C., Suzuki, M., Oida, H., Nunoura, T., et al. (2005). Enzymatic and genetic characterization of carbon and energy metabolisms by deep-sea hydrothermal chemolithoautotrophic isolates of Epsilonproteobacteria. Appl. Environ. Microbiol. 71, 7310–7320. doi: 10.1128/AEM.71.11.7310-7320.2005
Tang, K., Liu, K. S., Jiao, N. Z., Zhang, Y., and 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
Tang, K., Zhang, Y., Lin, D., Han, Y., Chen, C. A., Wang, D., et al. (2018). Cultivation-independent and cultivation-dependent analysis of microbes in the shallow-sea hydrothermal system off Kueishantao Island, Taiwan: unmasking heterotrophic bacterial diversity and functional capacity. Front. Microbiol. 9:279. doi: 10.3389/fmicb.2018.00279
Tarasov, V. G., Gebruk, A. V., Mironov, A. N., and Moskalev, L. I. (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
Tripp, H. J., Kitner, J. B., Schwalbach, M. S., Dacey, J. W., Wilhelm, L. J., and Giovannoni, S. J. (2008). SAR11 marine bacteria require exogenous reduced sulphur for growth. Nature 452, 741–744. doi: 10.1038//nature06776
Voordeckers, J. W., Starovoytov, V., and Vetriani, C. (2005). Caminibacter mediatlanticus sp. nov., a thermophilic, chemolithoautotrophic, nitrate-ammonifying bacterium isolated from a deep-sea hydrothermal vent on the mid-atlantic ridge. Int. J. Syst. Evol. Microbiol. 55, 773–779. doi: 10.1099/ijs.0.63430-0
Voss, M., Bange, H. W., Dippner, J. W., Middelburg, J. J., Montoya, J. P., and Ward, B. (2013). The marine nitrogen cycle: recent discoveries, uncertainties and the potential relevance of climate change. Philos. Trans. R. Soc. Lond. B Biol. Sci. 368:20130121 doi: 10.1098/rstb.2013.0121
Waite, D. W., Vanwonterghem, I., Rinke, C., Parks, D. H., Zhang, Y., Takai, K., et al. (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
Wirsen, C. O., Brinkhoff, T., Kuever, J., Muyzer, G., Molyneaux, S., and Jannasch, H. W. (1998). Comparison of a new Thiomicrospira strain from the mid-atlantic ridge with known hydrothermal vent isolates. Appl. Environ. Microbiol. 64, 4057–4059.
Xie, W., Wang, F., Guo, L., Chen, Z., Sievert, S. M., Meng, J., et al. (2011). Comparative metagenomics of microbial communities inhabiting deep-sea hydrothermal vent chimneys with contrasting chemistries. ISME J. 5, 414–426. doi: 10.1038/ismej.2010.144
Zhang, Y., Zhao, Z., Chen, C. T., Tang, K., Su, J., and 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
Zopfi, J., Kjaer, T., Nielsen, L. P., and Jorgensen, B. B. (2001). Ecology of Thioploca spp.: nitrate and sulfur storage in relation to chemical microgradients and influence of Thioploca spp. on the sedimentary nitrogen cycle. Appl. Environ. Microbiol. 67, 5530–5537. doi: 10.1128/AEM.67.12.5530-5537.2001
Keywords: shallow-water hydrothermal ecosystem, Kueishantao Islet, metatranscriptomics, 16S rRNA library, microbial community, metabolic pathway, biogeochemical cycle, coupling
Citation: Li Y, Tang K, Zhang L, Zhao Z, Xie X, Chen C-TA, Wang D, Jiao N and 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
Received: 08 June 2018; Accepted: 24 October 2018;
Published: 13 November 2018.
Edited by:Osvaldo Ulloa, Universidad de Concepción, Chile
Reviewed by:Luisa I. Falcon, Universidad Nacional Autónoma de México, Mexico
Alejandro A. Murillo, EMBL Heidelberg, Germany
Copyright © 2018 Li, Tang, Zhang, Zhao, Xie, Chen, Wang, Jiao and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yao Zhang, email@example.com
†These authors have contributed equally to this work