Skip to main content


Front. Microbiol., 22 February 2018
Sec. Aquatic Microbiology
Volume 9 - 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

Kai Tang1*†, Yao Zhang1†, Dan Lin1, Yu Han1, Chen-Tung A. Chen2, Deli Wang1, Yu-Shih Lin2, Jia Sun1, Qiang Zheng1 and Nianzhi Jiao1*
  • 1State Key Laboratory of Marine Environmental Science, Institute of Marine Microbes and Ecospheres, Xiamen University, Xiamen, China
  • 2Department of Oceanography, National Sun Yat-sen University, Kaohsiung, Taiwan

Shallow-sea hydrothermal systems experience continuous fluctuations of physicochemical conditions due to seawater influx which generates variable habitats, affecting the phylogenetic composition and metabolic potential of microbial communities. Until recently, studies of submarine hydrothermal communities have focused primarily on chemolithoautotrophic organisms, however, there have been limited studies on heterotrophic bacteria. Here, fluorescence in situ hybridization, high throughput 16S rRNA gene amplicon sequencing, and functional metagenomes were used to assess microbial communities from the shallow-sea hydrothermal system off Kueishantao Island, Taiwan. The results showed that the shallow-sea hydrothermal system harbored not only autotrophic bacteria but abundant heterotrophic bacteria. The potential for marker genes sulfur oxidation and carbon fixation were detected in the metagenome datasets, suggesting a role for sulfur and carbon cycling in the shallow-sea hydrothermal system. Furthermore, the presence of diverse genes that encode transporters, glycoside hydrolases, and peptidase indicates the genetic potential for heterotrophic utilization of organic substrates. A total of 408 cultivable heterotrophic bacteria were isolated, in which the taxonomic families typically associated with oligotrophy, copiotrophy, and phototrophy were frequently found. The cultivation-independent and -dependent analyses performed herein show that Alphaproteobacteria and Gammaproteobacteria represent the dominant heterotrophs in the investigated shallow-sea hydrothermal system. Genomic and physiological characterization of a novel strain P5 obtained in this study, belonging to the genus Rhodovulum within Alphaproteobacteria, provides an example of heterotrophic bacteria with major functional capacity presented in the metagenome datasets. Collectively, in addition to autotrophic bacteria, the shallow-sea hydrothermal system also harbors many heterotrophic bacteria with versatile genetic potential to adapt to the unique environmental conditions.


Approximately 50–60 shallow-sea hydrothermal systems at depths of less than 200 m below sea level are currently known, occurring near active coastal or submarine volcanoes, with systems located along arcs, mid-ocean ridges, and in island arc-related environments and even in continental margins (Tarasov et al., 2005). While extensive microbial surveys of the deep-sea hydrothermal systems have been conducted since the first discovery of deep-sea hydrothermal vents nearly 40 years ago (Corliss et al., 1979), attention has also been paid to their shallow-sea counterparts, which are much easier to access and can often be explored via SCUBA diving (Tarasov et al., 2005). The biological data for nearly 30 shallow-sea hydrothermal vent ecosystems have been published (Tang, 2014). Previous surveys of 16S rRNA genes using tag pyrosequencing and clone libraries have revealed the composition of microbial communities in shallow-sea hydrothermal environments, including the shallow-sea hydrothermal systems located at Kueishantao Island off NE Taiwan (Zhang et al., 2012; Tang et al., 2013), Eolian Islands (Vulcano Island and Panarea Island) in Italy (Manini et al., 2008; Maugeri et al., 2009, 2010, 2013b), Ambitle Island (Meyer-Dombard et al., 2012), D. João de Castro Bank, Azores (Mohandass et al., 2012), Milos Island in Greece (Brinkhoff et al., 1999; Sievert et al., 1999, 2000; Bayraktarov et al., 2013; Giovannelli et al., 2013; Price et al., 2013), Taketomi Island in Japan (Hirayama et al., 2007), and Eyjafjordur in Iceland (Marteinsson et al., 2001). These investigations showed that there were frequently a high abundance of autotrophs within the classes Gammaproteobacteria and Campylobacteria (previously termed Epsilonproteobacteria) (Waite et al., 2017) in the shallow-sea systems.

The submarine hydrothermal systems harbor chemolithoautotrophic bacteria and archaea communities typically associated with locally reduced gasses sulfide (H2S), methane (CH4), and hydrogen (H2) (Grzymski et al., 2008; McDermott et al., 2015; Brazelton and Baross, 2010), which are considered to support primary productivity through carbon dioxide (CO2) fixation. The sulfur-reducing chemolithoautotrophs Nautiliales-like organisms within Campylobacteria and sulfide-oxidizing chemolithoautotrophs Thiomicrospira-like organisms within Gammaproteobacteria dominated and exhibited distinct zonation within the water columns of the shallow hydrothermal system off Kueishantao Island, Taiwan (Zhang et al., 2012; Tang et al., 2013). Chemolithoautotrophs Nautiliales-like and Thiomicrospira-like organisms are frequently found in the other shallow-sea systems (Hirayama et al., 2007), and were found in some deep-sea hydrothermal system as well (Scott et al., 2006; Campbell et al., 2009; Brazelton and Baross, 2010; Yamamoto and Takai, 2011). They might possess the reductive tricarboxylic acid cycle (rTCA) and the Calvin-Benson-Bassham (CBB) cycle for carbon fixation in response to the available energy source in the form of the oxidation of reduced sulfur compounds and hydrogen in the environment, which is possibly fueled by geochemical energy with hydrogen and reduced sulfur, respectively (Tang et al., 2013). Chemolithoautotrophic organisms were also active within the steep geochemical gradients of the shallow-sea hydrothermal sediments, which were possibly involved in sulfide oxidation and sulfate reduction (Bayraktarov et al., 2013; Giovannelli et al., 2013). The sediments often harbored abundant Campylobacteria (such as genus Sulfurovum that encompasses sulfur- and thiosulfate-oxidizers bacteria) similar to those seen in deep-sea vents, but also other proteobacterial lineages that are distinct from those of deep-sea vents (Giovannelli et al., 2013). These studies provided some indication of a potential biogeochemical function for chemolithoautotrophic organisms in the shallow-sea hydrothermal system with significant ramifications for sulfur and carbon cycles (Zhang et al., 2012; Bayraktarov et al., 2013; Giovannelli et al., 2013; Tang et al., 2013).

Inorganic carbon is the primary carbon source assimilated by autotrophic bacteria in submarine hydrothermal systems, but hydrothermal fluids can also carry elevated concentrations of dissolved and particulate organic matter (Bennett et al., 2011; Yang et al., 2012). In addition, fluids may cool down and change physicochemical fluctuating conditions by mixing rapidly with seawater (Yang et al., 2012). The enrichment of nutrients and the temperature and dissolved oxygen gradients in the mixed hydrothermal fluids and seawater might support active heterotrophic microbes. Recent studies have shown that deep-sea hydrothermal systems are inhabited by versatile heterotrophic Alphaproteobacteria and Gammaproteobacteria, which are significantly distinct from heterotrophic lineages common in the deep-sea environment, some of which have the potential for alkane degradation (Meier et al., 2016). A bacterium closely related to a human pathogenic Vibrio species was isolated from surrounding sulfide chimneys of a hydrothermal vent along the East Pacific Rise, and its genomic information provided new insights on how species adapt to the deep-sea environment (Hasan et al., 2015). Several studies have shown the presence of physiologically, metabolically, and phylogenetically diverse heterotrophic communities in shallow-sea hydrothermal systems (Giovannelli et al., 2013; Meyer-Dombard et al., 2013; Price et al., 2013), in which some seemed mostly to be involved in arsenic and iron redox cycling (Meyer-Dombard et al., 2013; Price et al., 2013). However, in marine hydrothermal systems, compared to chemolithoautotrophic organisms, heterotrophic bacteria distribution, diversity, and metabolic capacity have still been poorly investigated, and isolates with reference genomes from the shallow-sea hydrothermal system are scarce.

In the present work, we revisited the shallow-sea hydrothermal system near Kueishantao Island and collected samples from the water column above the yellow and white vents, in addition to samples of sandy sediments, rocks, and dead vent deposits near the vent sites. Microbial compositions were characterized using culture-independent methods including fluorescence in situ hybridization (FISH) and 16S rRNA gene amplicon sequencing. Metagenomic analysis was used to determine the functional and genetic potential of bacteria. The various heterotrophic bacteria with potentially different trophic strategies were isolated from the investigated shallow-sea hydrothermal systems. Some of these isolates may be previously undiscovered bacterial species. Further studies aimed to characterize the physiological and genomic features of bacterial strain P5, a bacterioplankton species with a considerable capacity for adaptation to the shallow-sea hydrothermal environment.

Materials and Methods

Sample Collection

Sampling was performed by scuba divers in May 2015 in two different regions of the shallow-sea hydrothermal system located near Kueishantao Island (N 24.834, E 121.962), Taiwan, China: the white vent area and the yellow vent area (Figure 1). All necessary permits were obtained for the described field studies. The vents were observed at depths of approximately 10 and 7 m, respectively. Four samples were collected from the following locations (water samples): W_0m, W_5m, W_surface, and W_outside, which are 0 m above the white vent, 5 m above the white vent, the surface water directly above the vent, and surface water 6 m laterally away from the white vent, respectively. Four more samples were collected from the following locations (water samples): Y_0m, Y_5m, Y_surface, and Y_outside, which are 0 m above the yellow vent, 5 m above the yellow vent, the surface water directly above the yellow vent, and the surface water 6 m laterally away from the yellow vent. The sandy and rocky sediments nearby the yellow vent were named YS_S and YS_R, respectively. The sandy and rocky sediments on the seafloor were named S_S and S_R. The rocky sediment collected nearby the dead vent was named DS_R (an expired vent without the plume and gas discharging). The 2-l water samples were filtered through 3 and 0.2 μm pore-size polycarbonate filters (EMD Millipore Corp., Darmstadt, Germany) for further analysis. Sediment samples were frozen on site and kept frozen during transportation and storage.


FIGURE 1. Pictures of the shallow-sea hydrothermal system off Kueishantao Island, Taiwan. (A) Map of the sample sites and magnified view of the sample map of the locations indicated by YV (yellow vent area) and WV (white vent area); (B) the surface water in the shallow-sea hydrothermal system.

Abundance of Bacteria Determined by CARD-FISH

Fifty-milliliter water samples were immediately fixed with freshly prepared paraformaldehyde (4% final concentration) and stored at 4°C overnight prior to filtration through 0.2 μm pore size polycarbonate filters (Millipore, United States). Filtered samples were stored at -20°C for later analysis by FISH with horseradish peroxidase-labeled oligonucleotide probes (CARD-FISH). Picoplankton abundance was determined by DAPI staining, and bacteria were enumerated by CARD-FISH. Filters were embedded in low-gelling-point agarose and incubated with either lysozyme for the bacteria probe mix (Eub338: 5′-GCT GCC TCC CGT AGG AGT-3′, Eub338II: 5′-GCA GCC ACC CGT AGG TGT-3′ or Eub338III: 5′-GCT GCC ACC CGT AGG TGT-3′) or for the negative control probe (Non338: 5′-ACTCCT ACG GGA GGC AGC-3′) (Teira et al., 2004). Filters were cut into four pieces and hybridized with HRP-labeled oligonucleotide probes and tyramide Alexa488 for signal amplification following the previously described protocol (Teira et al., 2004).

DNA Extraction

Total DNA of water samples was extracted from the samples with the FastDNA SPIN Kit for Soil (MP Biomedicals, Illkirch, France) according to the manufacturer’s instructions. DNA from the sediment samples was isolated using PowerMax Soil DNA Isolation Kits (MoBio Laboratories Inc., United States) according the manufacturer’s instructions. The concentration and purity of DNA were evaluated using a NanoDrop spectrophotometer (ND-1000, Thermo Fisher Scientific, Waltham, MA, United States). DNA extracts were stored at -80°C until further analysis. The following sequencing was performed at the Chinese National Human Genome Center in Shanghai.

Amplification and Sequencing of the Bacterial 16S rRNA Gene

Amplicons were generated using fusion degenerate primers 343F (Wilson et al., 1990) and 798R (Rochelle et al., 1995) with ligated overhang Illumina adapter consensus sequences. The libraries were prepared in accordance with the instructions included with the Illumina Nextera XT Index kit (Illumina, United States). Pooled amplicons were purified with the Agencourt AMPure XP purification system (Beckman, United States) and analyzed with an Agilent bioanalyzer 2100 (Agilent Technologies, United States) to confirm appropriate amplicon size. Finally, amplicons were diluted to 10 nM, quantified and sequenced on the Illumina MiSeq platform (reagent kit v.3; Illumina, United States).

Metagenomes Sequencing

For metagenome sequencing, 1 μg of sample DNA was sheared to 500 bp by Covaris M220 (Covaris, United States). The library was constructed by NEBNext UltraTM DNA Library Prep Kit (NEB, United States). Finally, 10 nM sequencing library was used generate cluster in cBot using TruSeq PE Cluster Kit (Illumina, United States), and sequenced by Illumina Hiseq 2500 for 2 × 250 bp data.

16S rRNA Gene Sequence Data Analysis

Acquired Illumina reads were filtered by Meta Genome Rapid Annotation using Subsystem Technology (MG-RAST1) QC pipelines to remove the replicated reads (Meyer et al., 2008). The filtered reads were used for the following bioinformatic analysis. For taxonomic analysis, the SILVA small subunit (SSU) database implemented in MG-RAST was used as annotation source for 16S rRNA reads to analyze the bacterial populations in samples using an E-value cutoff of 1e-05, minimum identity cutoff of 60%, and minimum alignment length cutoff of 150 bp.

Shotgun Metagenomic Sequence Data Analysis

Functional profiles were identified using the SEED subsystems annotation source of the MG-RAST, with 1e-05 as maximum e-value, a minimum identity of 60%, and a minimum alignment length of 15 amino acids. To remove the bias of average genome size on the sampling of genes from a given metagenomic community, the raw functional gene hits were normalized to the number of recA gene encoding recombinase A hits in the respective database. The PRIMER-6 package was used to calculate the Bray–Curtis similarity matrices of metagenomes and generate non-metric multidimensional scaling plots. PERMANOVA implanted in PAST v3.05 was carried out to compare samples from each environment (Hammer et al., 2001). To determine whether the relative abundances of functional genes differed significantly between sample categories, we conducted multiple t-tests with P-values calculated using a Holm–Sidak correction (Shaffer, 1995) for multiple comparisons implemented in Prism v6.

Strain Isolation and Culture

All reagents used in bacterial cultures were obtained from Sigma–Aldrich (United States) unless otherwise specified. For the cultivation of bacteria, 200 μL of the water sample was spread onto at least one of the agar plates. The widely used medium for routine culture of marine bacteria were selected (Joint et al., 2010), including low-nutrient R2A (Difco, United States), and nutrient-rich 2216E (Becton-Dickinson, United States), RO (1 g peptone, 1 g yeast extract, 1 g natrium aceticum, 1 g sodium acetate per liter artificial seawater (ASW) with vitamins and trace elements, pH 7.8–8.0), SYPG (30 g NaCl, 0.5 g yeast extract, 0.25 g tryptic peptone, 0.1 g glutamic acid monosodium salt per liter, pH 7.5), SC (1 g D-glucose, 1 g alkapolpeg-600, 1 g L-malic acid, 1 g D-aspartic acid, 1 g yeast extract per liter ASW with vitamins and trace elements, pH7.2–7.5), NS (0.5 g sodium nitrate, 0.5 g sodium sulfate, 0.1 g yeast extract per liter artificial seawater with vitamins and trace elements, pH 7.2–7.5), YTSS (2.0 g yeast extract, 1.25 g tryptone, 20 g sea salt (Sigma) per liter distilled water and autoclave, pH 7.0). All agar plates were incubated at 28°C or under anaerobic conditions at 24°C. Bacterial colonies were picked from the plates and purified further on 2216E plates. The bacterial 16S rRNA gene sequencing and sequence analysis were performed using a previous method (Embley, 1991). The sequence data have been submitted to the NCBI database under accession numbers AM988866 through AM989325. Genomic DNA was extracted from 200 μL of culture using a TIANamp Bacteria DNA Kit (Tiangen, China). The primer pair pufLF (5′-CTKTTCGACTTCTGGGTSGG-3′) and pufMR (5′-CCATSGTCCAGCGCCAGAA-3′) were designed to amplify pufL and pufM through PCR to allow detection of the photosynthetic reaction center genes in the strains (Cho et al., 2007).

Genome Sequencing and Analysis

Whole genome sequencing of strain P5 was accomplished using a hybrid approach, combining Illumina short read data with PacBio long read data (Koren et al., 2012). For PacBio sequencing, 5 μg of sample DNA was sheared to 10 Kb by a Covaris® g-TUBE® (Covaris, United States). A PacBio® SMRTbellTM Template Prep Kit (PacBio, United States) was used to construct the library. The sequencing primers were annealed using a PacBio DNA/Polymerase Kit (PacBio, United States) and polymerase combined with the SMRTbell templates. We obtained long read data from PacBio RS II PacBio RS II (PacBio, United States). The raw Illumina data were filtered by the FASTX-Toolkit to remove the adapters, N bases, and low-quality reads. The clean data were assembled using Velvet v1.2.03 with default parameters. The PacBio long reads were assembled by RS HGAP assembly 3. The complete genome was finally gap closed by Sanger sequencing. The final assembled genomes of P5 were automatically annotated and analyzed through the IMG/ER2. The comparison and visualization of multiple genomes was conducted with BRIG (Alikhan et al., 2011).

Physiological and Biochemical Analysis of Strain P5

Strain P5 was isolated from surface seawater off Kueishantao Island, northeast Taiwan by 2216E medium. The strain was adaptively grown for 3 days (inoculated into a new bottle everyday) in ASW containing a 1 mM concentration of dissolved organic carbon (DOC). The DOC concentration was adjusted by ASW supplemented with full strength medium (5 g peptone and 1 g yeast extract per liter ASW and the final concentration was determined using a total organic carbon analyzer (Shimadzu, Japan). Bacteria (105 cells/mL) were then washed three times with autoclaved ASW and used to inoculate growth medium consisting of artificial seawater (ASW) base combined with substrates including 0.45 μM DOC, 2.5 mM NaHCO3, and 1 mM Na2S2O3 in phosphate-buffered saline (pH 7.4) (autotrophic culture: NaHCO3 and Na2S2O3; heterotrophic culture: DOC). Cultures were incubated at 22°C at 160 rpm/min in the darkness or light (12 h)/dark (12 h) cycle for 4 days. To determine bacterial cell density, cultures were stained with SYBR Green I (1:100 dilution; Molecular Probes, United States) for 15 min, and measured via flow cytometry (BD Accuri C6, United States). All culture experiments were performed in triplicate.

The strain P5 was cultured in 2216E and freeze-dried by a Freeze Dry System (Labconco Corp., Czechia) for 48 h, a moderate amount of chloroform was added, and cells were broken under a Ultrasonicator SM-650D (Shunma, China). The Ultrasonicator was operated for a total of 10 min containing 9 s for running, and 5 s for stopping. The power was 22% (total power 650 W). The supernatant was measured using a UV-Vis spectrophotometer (Agilent 8453, United States) after being centrifuged at 7006 g/min for 10 min. Bacteriochlorophyll a was measured with a FIRe Fluorometer Induction and Relaxation System (Satlantic, Canada) at 800 nm.

Nucleotide Sequence Accession Numbers

The metagenomic datasets are publicly available in the MG-RAST system under project identifiers W_0m (4671689.3 and 4668298.3), W_5m (4671690.3 and 4668304.3), W_surface (4671687.3 and 4668303.3), W_outside (4671688.3 and 4668305.3), Y_0m (4671691.3 and 4668302.3 for 16S rDNA and functional metagenome, respectively), Y_5m (4671692.3 and 4668301.3), Y_surface (4671693.3 and 4668300.3), Y_outside (4671694.3 and 4668299.3), YS_S (4647123.3 and 4644370.3), YS_R (4647124.3 and 4644369.3), S_S (4647121.3 and 4644371.3), S_R (4647122.3 and 4644368.3), and DS_R (4647120.3 and 4644372.3). They have also been deposited in the NCBI Short Read Archive: 16S rDNA reads, SRR5229862-SRR5229869 for W_0m, W_5m, W_surface, W_outside, Y_0m, Y_5m, Y_surface, and Y_outside, respectively, YS_S (SRR5149598), YS_R (SRR5149595), S_S (SRR5149599), S_R (SRR5149606), DS_R (SRR5149594), metagenomic DNA reads, SRR5229878-SRR5229885 for W_0m, W_5m, W_surface, W_outside, Y_0m, Y_5m, Y_surface, and Y_outside, respectively, YS_S (SRR5149596), YS_R (SRR5149597), S_S (SRR5149607), S_R (SRR5149602), DS_R (SRR5149604). The complete genome sequence of strain P5 has been deposited into GenBank under the accession numbers CP015039 (chromosome) and CP015040-CP015043 (plasmids).

Results and Discussion

Environmental Parameters

The fluids in the shallow-sea hydrothermal system were slightly acidic (approximately pH 6) and oxic (Supplementary Table S1). Salinity, dissolved inorganic carbon, nitrate, nitrite, and phosphate concentrations were nearly oceanic (Supplementary Table S1). Elemental sulfur (S0) is naturally enriched in the shallow-sea hydrothermal fluids near Kueishantao Island. Compared to deep-sea vents, the sulfide concentration was lower in this shallow vent system (Chen et al., 2005; Chan et al., 2016). Compared to the common marine environment, the concentartion of DIC was higher, while the DOC centration was similar in the shallow-sea hydrothermal system (Supplementary Table S1). These physico-chemical properties of the hydrothermal system might support not only autotrophic bacteria, but also heterotrophic bacteria.

Cell Densities

The DAPI-based total cell counts were on average on the order of 106 cells mL-1 (Supplementary Table S2). The general bacterial probe mixture EUB338 hybridized 90–94% of all DAPI-stained cells in the white vent area, while the percentage of bacterial probe in the yellow vent area was 57–87% (Supplementary Table S2). Compared to the white vent (temperature up to 58°C inside the vent), the relatively high temperature (up to 116°C) near the yellow vent might result in lower bacterial abundance and their proportions of total microbial cells at the Y_0m site (0 m above the yellow vent) (only 57%).

Bacterial Community Composition

The investigation of the microbial communities using 16S rRNA gene amplicon sequencing revealed that a significant proportion of the sequences in the water samples (14–89%) and sediment samples (21–47%) could either not be unassigned or were assigned to unclassified groups, indicating the presence of so far uncharacterized bacteria within this geothermal ecosystem. Campylobacteria were apparently present in relatively high abundance in the sediment samples, and represented the main bacterial groups at the Y_0m site and the W_0m 318 site (water samples, 0 m above the vents) (Figure 2). However, the relative abundance of 16S rRNA sequences of Campylobacteria decreased in the water column above the vents (Figure 2). Few sequences of Campylobacteria were classified into those from sulfur-reducing chemolithoautotrophs Nautiliales-like organisms that were previously known to be a dominant group in the white vent (Zhang et al., 2012; Tang et al., 2013). The shallow-sea hydrothermal vents varied greatly in chemical composition of gasses discharging such as CO2 and H2S (Chen et al., 2016). In the previous one, the DIC concentrations was about 6 mM in 0 m above the white vent (Zhang et al., 2012), however, there was only a little higher (∼2.7 mM) than that at nearby seawater in this study. On the other hand, the H2S concentrations have declined since 2014 (Chen et al., 2016). In addition, physical conditions such as tide could influence the hydrothermal ecosystem (Chen et al., 2005). These might influence on the microbial community composition. Gammaproteobacteria were found abundantly at the surface waters and outside the white vent (W_surface site and W_outside site) (Figure 2). However, the very few sequences of sulfur-oxidizing Gammaproteobacteria Thiomicrospira-like organisms were detected. Members of the family Alteromonadaceae and Oceanospirillaceae within Gammaproteobacteria are the most abundant organisms in W_surface and W_outside, respectively. Y_5m (5 m above the yellow vent) and Y_surface (the surface water immediately above the yellow vent) have abundant members of Alphaproteobacteria (Figure 2). Flavobacteriia were present at all sites (Figure 2). There are some relatively rare taxa found in the hydrothermal system including Deltaproteobacteria, Beatproteobacteria, Actinobacteria, Cyanobacteria, Aquificae, Chorobi, Spirochaetes, and Tenericutes. Deltaproteobacteria, and Firmicutes are mainly found in sediments samples (Figure 2).


FIGURE 2. The main class-level bacterial composition among the samples. The map plot depicts the relative percentage of each bacteria within each sample. Water samples: W_0m, W_5m, W_surface, W_outside, Y_0m, Y_5m, Y_surface, and Y_outside; sediment samples: YS_S, S_R, S_S, S_R, and DS_R.

Overall, homology-based taxonomic assignments of shotgun metagenomics data predicted by the MG-RAST showed that Gammaproteobacteria and Alphaproteobacteria were abundant classes in the white and yellow vent water samples, making up approximately 60% of the microbial assemblage, whereas Campylobacteria were abundant class in the vent sites and sediment samples. The most dominant proteobacterial orders overall were Rhodobacterales and Rhizobiales within the class Alphaproteobacteria; Alteromonadales, Pseudomonadales, and Oceanospirillales within the class Gammaproteobacteria. However, many shotgun metagenomic sequences (mostly > 50% of the total sequences) still remained functionally unannotated and taxonomically unassigned due to the limitations of the available tools and the paucity of reference genomes. These could possibly lead to a discrepancy between shotgun metagenomics data taxonomic assignment and 16S rRNA-based classification.

Metabolic Potentials and Functional Biomarkers of Microbial Communities

Based on the relative abundance of SEED subsystems, multidimensional scaling plots show that the hydrothermal water samples in this study clustered together at the functional level, and were separate from the hydrothermal sediments samples (Figure 3 and Supplementary Table S3); permutational multivariate analysis of variance (PERMANOVA) indicates that this difference is significant (P = 0.009). We observed a set of SEED functions “Carbohydrates” and “Amino Acids and Derivatives” that contributed strongly to the difference between the communities in the sediments and water samples; that set of SEED functions is more abundant in the hydrothermal water samples. Compared to a metagenome from the previous investigation in the white vent area, “Carbohydrates” were more abundant in the metagenomic datasets in this study, but there were fewer “Protein metabolism” results. Compared to the bacterial community compositions from the previous investigation, Flavobacteriia were more abundant in this study. The Flavobacteriia genomes contain diverse and abundant glycoside hydrolases (GHs) that are involved in carbohydrate metabolisms (Tang et al., 2017). Two categories, “Carbohydrates” and “Protein metabolism,” together contributed to approximately 26% of the differences between communities in the two metagenomes. However, “Sulfur metabolism” and “Photosynthesis” contribute less to their differences. The hydrothermal water samples were functionally distinct from those in the costal and open seawater samples (PERMANOVA, P < 0.01 in all cases). The genes associated with the cell wall and capsule category, and virulence, disease, and defense category were more abundant in the hydrothermal water samples than those in the costal and open seawater samples metagenomes, which are main contributors to the dissimilarity between metagenomes (Supplementary Table S3).


FIGURE 3. Multidimensional scaling plots of samples using Bray–Curtis similarity according to the SEED subsystem. Color represents different sample areas and each habitat label type (n) is indicated [data highlighted by two circles from this study; W_0m (2010) and W_surface (2010) data from our previous study (4) (Tang et al., 2013), H1 and H2 data from references (Grzymski et al., 2008; McDermott et al., 2015), other data from a reference (Rusch et al., 2007), Details in Supplementary Table S3. Samples from each of the respective environments clustered together based on their functional profile. The stress values are reported in the top right corner and represent the goodness-of-fit.

The functional metagenomes suggested that oxidation of reduced sulfur compounds in the hydrothermal system could occur through the Sox multienzyme system (Figure 4), which catalyzes the complete oxidation of reduced sulfur compounds to sulfate (Yamamoto and Takai, 2011). Several enzymes have been proposed to as possible oxidizers of inorganic sulfur compounds, including sulfide quinone oxidoreductase and sulfite oxidase, which oxidize sulfide to elemental sulfur and sulfite to sulfate, respectively (Figure 4) (Yamamoto and Takai, 2011). Genes encoding polysulfide reductase (Psr) are present in the sediments metagenome, resulting in the reduction of polysulfide derived from elemental sulfur to sulfide (Yamamoto and Takai, 2011). Previous metagenome analyses have suggested that Campylobacteria in the vent can gain energy from sulfur-reduction catalyzed by Psr to fix CO2 by the rTCA cycle (Tang et al., 2013). However, the relative gene abundance of Psr in the water samples was lower than those in the sediments (Figure 4). The Deltaproteobacteria contributed to Psr gene sequences in the sediment metagenomes. In addition, few key genes encoding ATP-dependent citrate lyase for the rTCA cycle (Hügler and Sievert, 2011) were found in any of the metagenome datasets in this study, whereas genes encoding ribulose-1,5-bisphosphate carboxylase (RuBisCO) were found; those genes mediated the CBB cycle for carbon fixation. The relative abundances of RuBisCO (normalized to recA genes) in the hydrothermal water samples, on average, were 0.12 and 0.14 in the yellow vent and white vent areas, respectively, which were both lower than found (0.98) in a previous investigation (Tang et al., 2013). In contrast, heterotrophic metabolism might be predominant both in the water column and the sediment, as indicated by a set of transporters, peptidases, and GHs genes in all the metagenomic datasets (Figure 4). Although the relative abundances of genes encoding GH and peptidase between the sediment samples and water samples display no statistical difference (t-test, P > 0.01), they are more diverse in the waters than in the sediment of hydrothermal system, indicating a wider spectrum of substrate utilization for bacteria in the waters of this shallow-sea system. Similarly, more ATP-binding cassette systems (ABC transporters), symporters, and TonB-dependent receptors (TBDTs) were found in the hydrothermal waters, which allowed them to import organic matter efficiently (Figure 4). All metagenomes were found to possess exporters involved in the efflux of heavy metals and metabolites (Supplementary Table S4). Predicted substrates for the transport systems in the metagenomes include a variety of carbohydrates, carboxylic acids, amino acids, peptides, metals, and other nutrients (Supplementary Table S4).


FIGURE 4. Heat map of functional composition among the samples. The value of functional gene relative to recA (single copy control gene) is assigned with a color relative to the maximum value among all comparisons of each gene. The colors represent the minimum (blue), middle (white, 0.5), and maximum (red) values listed on the right. The number in the box represented the components of each functional category.

All metagenomes were found to contain genes encoding sulfate adenylyltransferase, adenylsulfate kinase, and adenylsulfate reductase, which are required for assimilatory sulfate reduction to supply sulfur for biosynthesis in aerobic marine bacteria (Figure 4). Functional metagenomic analyses indicated that the relative abundance of genes encoding Psr, Ni-Fe hydrogenase, and periplasmic nitrate reductase exhibited statistically significant differences between hydrothermal sediments and water samples (t-test, P < 0.05) and it was more abundant in sediment samples. The genes encoding Ni–Fe hydrogenase in the metagenome enabled bacteria to use H2 as an energy source (Petersen et al., 2011). There is some possibility that nitrate could be used as alternative electron acceptor in the presence of periplasmic nitrate reductase (Canfield et al., 2010). Thus, chemolithotrophs might contribute to the chemical transformations of elements in the sediment.

Overall, the Rubisco gene sequences were mainly homologous to those in the order Rhodobacterales, Rhizobiales, Methylococcales, and Thiotrichales, in contrast to previous metagenomes where most of the Rubisco gene sequences homologs were affiliated with Thiomicrospira-like (Tang et al., 2013). Cyanobacteria also contribute to Rubisco gene pool at the sea surface of the hydrothermal system. The order Desulfovibrionales within Deltaproteobacteria in the vent sites and sediment samples contributed to genes encoding NiFe hydrogenase. The majority of genes encoding TBDTs were closely related to Flavobacteriia and Gammaproteobacteria (primarily members of the Alteromonadales and the Pseudomonadales) in the water samples, while TBDTs sequences from the Campylobacteria and Deltaproteobacteria class were found in sediment samples. The members of the order Rhodobacterale contributed greatly to the diversity of ABC transporters in the metagenomes, and the members of Flavobacteriia contribute greatly to the diversity of peptidases and glycoside hydrolases in the metagenomes. The sox sequences in the metagenomic data were mainly homologous to those found in the orders Rhodobacterale and the class Campylobacteria. Genes encoding periplasmic nitrate reductase and dissimilatory nitrite reductase were more prevalent in the sediment dataset than in the water sample dataset, and they were mainly homologous to sequences from the members of Campylobacteria and Deltaproteobacteria.

Taxonomic Assignment of the Isolated Strains

A total of 408 isolates were identified by near full-length 16S rRNA gene sequence analysis (Supplementary Table S5). Of these isolates, 187 were from the yellow-vent area, 111 were from the white-vent area, and 65 were from the surface water above the vents. It was noted that 42 isolates were obtained from the waters above the dead vent. A previous study showed that heterotrophic bacteria dominated inactive deep-sea hydrothermal system (Sylvan et al., 2012). The isolates were distributed into seven different bacterial classes from four different phyla, namely Proteobacteria, Bacteroidetes, Actinobacteria, and Firmicutes. However, the failure to isolate Epsilonprobacteria strains may be because the optimal conditions for their culture have yet to be established. The majority of the characterized strains belonged to the Alphaproteobacteria (174 strains), Gammaproteobacteria (145 strains), and Actinobacteria (46 strains). The rest of the strains were assigned to Betaproteobacteria, Flavobacteriia, Cytophaga, and Bacillus. The spectrum of different genera was greatest within the classes Alphaproteobacteria (31 genera), followed by Gammaproteobacteria (15 genera), Flavobacteriia (14 genera), and Actinobacteria (12 genera) (Figure 5). For the classes Alphaproteobacteria, Erythrobacter (81 strains) and Paracoccus (31 strains) were the predominant genus (Figure 5). For the classes Gammaproteobacteria, Vibrio (41 strains) and Pseudoalteromonas (36 strains) constituted a high proportion (Figure 5). For the classes Actinobacteria, Microbacterium (16 strains) was dominant. These taxonomic groups are commonly detected in ocean environments (Haggerty and Dinsdale, 2017).


FIGURE 5. The diversity of the bacterial isolates. (A) Gammaproteobacteria, (B) Alphaproteobacteria, (C) Flavobacteriia, and (D) Actinobacteria. Each color represents the percentage of the taxon in the total isolates.

These isolates may possess novel functions and exhibit a broad range of ecological attributes and life-history strategies. For example, some strains belonged to oligotrophic taxa (including members of Erythrobacteraceae and Sphingomonadaceae groups); some strains belonged to the copiotrophic taxa (including members of Vibrio and Alteromonadaceae phyla) and may be able to rapidly consume labile DOC (Church, 2009; Kirchman, 2015). The members of Flavobacteriaceae have the ability to degrade complex high-weight molecule organic compounds (Buchan et al., 2014). Those in the genera Erythrobacter and Novosphingobium can to metabolize nutrient-poor and recalcitrant carbon substrates (Church, 2009; Kirchman, 2015). Furthermore, pufL and pufM genes, which encode photoreaction center L and M polypeptides, respectively, in Erythrobacter were identified as aerobic anoxygenic phototrophic bacterial gene biomarkers. Previous studies have suggested that they have the capacity to undergo photoheterotrophy in marine environments (Jiao et al., 2007).

On the basis of the taxonomic assignment by the EzTaxon classifier using annotated 16S rRNA gene sequences (Chun et al., 2007), a total of 58 strains from the classes Alphaproteobacteria, Gammaproteobacteria, Flavobacteriia and Actinobacteria might represent new bacterial species not yet validly described with less than 98% maximum identity with their closest BLAST hits (Supplementary Table S5). Among them, the information on the first complete genome of a cultivated actinomycete strain JLT9 isolated from the shallow-sea hydrothermal system was reported, which contained various sulfur oxidation genes (Figure 6A) (Han et al., 2017). The complete genome of the Maribacter sp. T28 harbored the xylanolytic, alginolytic and pectinolytic enzymes responsible for polysaccharide degradation (Figure 6B) (Genbank accession, CP018760) (Zhan et al., 2017). The genome data suggested that Marivivens sp. JLT3646 within Alphaproteobacteria has the potential to degrade aromatic monomers (Figure 6C) (Genbank accession, CP018572 and CP018573 for chromosome and plasmid, respectively) (Chen et al., 2017).


FIGURE 6. Transmission electron micrograph images of representative strains isolated from the shallow-sea hydrothermal system. (A) Strain JLT9 belonging to the genus of Serinicoccus within Actinobacteria, (B) Strain T28 belonging to the genus of Maribacter within Flavobacteriia, (C) Strain JL3646 belonging to the genus of Marivivens within Alphaproteobacteria, (D) Strain P5 belonging to the genus of Rhodovulum within Alphaproteobacteria.

The gene sequences homologs affiliated with Rhodobacterales were abundant, accounting up to 20.57% of the total sequences in the water sample. Strain P5 (Figure 6D) exhibited a 97.48% 16S rRNA sequence similarity with Rhodovulum adriaticum DSM 2781 within the order Rhodobacterales (Figure 7). Total of eight strains of Rhodovulum were isolated form the white and yellow vents area (Supplementary Table S5). The members of Rhodovulum were also frequently found near the shallow-sea submarine vents of Panarea Island (Maugeri et al., 2013a). The newly discovered strain P5 of great metabolic versatility could be considered to be representative of heterotrophic bacteria isolate from the shallow-sea hydrothermal system.


FIGURE 7. Phylogenetic tree based on alignment of 16S rRNA gene sequences showing the relationship of strain P5 (shown in bold print) within the genus Rhodovulum. The trees were constructed by the Neighbor-Joining method using MEGA version 6.0 software (closed circles represent the same tree branch using the algorithm of Maximum Likelihood and Maximum Parsimony) and rooted by using Rhodospirillum rubrum ATCC 11170T as the outgroup. Numbers at nodes represent bootstrap values (based on 100 resamplings). Bootstrap values more than 50% are noted. The GenBank accession numbers for 16S rRNA gene sequences are shown in parentheses. Bar, 1 nt substitution per 100 nt.

Biochemical and Physiological Characteristics of Strain P5

Strain P5 is a Gram-negative, spindle-shaped, purple bacterium (1.5–2.0 μm in length and 0.9–1.0 μm in width) (Figure 6D) that can grow at the temperature range of 20–40°C (optimum, 26–34°C), in the pH range of 5–9 (optimum, 6–8), in the presence of 0–6.0% (w/v) NaCl (optimum, 3.0%), and in aerobic or microanaerobic conditions (Supplementary Table S6). Electron microscopy of ultrathin sections revealed the presence of a vesicular type internal photosynthetic membrane that is a common feature in other Rhodovulum species (Figure 8A) (Masuda et al., 1999; Srinivas et al., 2007; Kompantseva et al., 2010; Maugeri et al., 2013a). The in vivo absorption spectrum of intact P5 cells exhibited four major peaks at 508, 588, and 804 and 850 nm, thus confirming the presence of bacteriochlorophyll a and carotenoids (Figure 8B). The functionality of the photosynthetic apparatus was investigated by infrared kinetic fluorescence measurements. P5 displayed clear induction of bacteriochlorophyll a with FV/FM 0.749 ± 0.003 (mean ± SD; n = 3), which confirmed that its fully functional photosynthetic reaction centers are connected to an efficient electron transfer chain (Figure 8C). These results suggested that P5 had phototrophy capability similar to other Rhodovulum species (Masuda et al., 1999; Srinivas et al., 2007; Kompantseva et al., 2010; Maugeri et al., 2013a). Strain P5 is capable of heterotrophic, autotrophic and photoheterotrophic growth (Figure 8D).


FIGURE 8. Morphological, biochemical and physiological characterization of bacterial strain P5. (A) Transmission electron micrograph images of strain P5. The external surface of the cytoplasmic membrane of strain P5 is of vesicle type. (B) In vivo absorption spectrum of strains P5 and reference strain. (C) Fluorescence induction and relaxation kinetics recorded by infrared fluorometry. (D) Growth rates of bacteria.

Genomic Features of Strain P5

Strain P5 exhibited a genome of 4,137,334 bp (one chromosome and four plasmids) with a G+C content of 64.64% that coded 4,243 protein-coding and 68 RNA genes (including nine rRNA operons). Bidirectional BLAST analyses showed that approximately 3,000 genes exhibited > 50% sequence identity between strain P5 and two sequenced genomes of Rhodovulum species [Rhodovulum sulfidophilum DSM 1374 (Nagao et al., 2015b) and Rhodovulum sulfidophilum DSM 2351 (Nagao et al., 2015a)] (Figure 9), indicating identical or equivalent function.


FIGURE 9. Multiple genome comparison of strain P5 (chromosome) and Rhodovulum sulfidophilum strains. The composite genome comparison Figure was generated using BRIG after performing a BLASTn analysis of strain P5 as the reference genome. Each genome mapping to the reference is represented as a colored ring, with a solid color representing greater than 50% sequence identity. The innermost rings show GC content (black) and GC skew (purple/green) in strain P5. The remaining rings indicate genome regions of R. sulfidophilum DSM 2351 (green), R. sulfidophilum DSM 1374 (blue), and strain P5 (purple). The predicted major functional genes in P5 are labeled on the outermost ring with black arcs.

Strain P5 was predicted to possess complete central carbon metabolic pathways, including glycolysis, the pentose phosphate pathway, the Entner–Doudoroff pathway, and the tricarboxylic acid cycle. The P5 genome harbored genes encoding 16 GHs and 30 peptidases as well as a urease gene cluster, indicating the potential for the degradation of carbohydrates, proteins, peptides, and urea (Supplementary Table S7). The most abundant transporter systems in the P5 genome were ABC transporter components, followed by TRAP-type transport components, and the exporters and antiporters, accounting for 8.1% of the total protein-coding genes of the chromosome in P5. Predicted uptaking substrates for the complete transporter systems comprised a variety of carbohydrates, carboxylic acids, amino acids, peptides, alkanesulfonate, metals, and other nutrients. The genome harbored two genes sets encoding type IV secretion machinery systems that are used to deliver proteins and DNA into the extracellular environment, whereas the other two Rhodovulum species lack these genes (Figure 9). This organism has genes for the complete CBB required for autotrophic carbon fixation. Well-known electron donors utilized for microbial autotrophic growth are hydrogen and reduced forms of sulfur (sulfide, S0, and thiosulfate). Strain P5 possessed the complete repertoire of genes for the oxidation of reduced sulfur compounds, which increased niche-specialization competitiveness in the shallow-sea environments. These genes encoded enzymes for the oxidation of reduced sulfur compounds including the Sox enzyme complex for oxidation of reduced sulfur to sulfate (SO42-), sulfide quinone oxidoreductase, mediating the oxidation of sulfide to elemental sulfur, rhodanese sulfurtranferase for oxidation of thiosulfate (S2O32-) to sulfite (SO32-), and reverse dissimilatory sulfite reductase for oxidation of elemental sulfur to sulfite, adenosine 5-phosphosulfate reductase and sulfate adenylyltransferase for oxidation of sulfite to sulfate (SO42-). Strain P5 has two gene clusters encoding different Ni–Fe hydrogenases, which is predicted to catalyze the reversible oxidation of hydrogen gas and enables bacteria to use molecular hydrogen as a source of energy. The heterodisulfide reductase-related proteins (Hdr) are only identified in the P5 genome (Figure 9), and are likely candidates to be involved in energy coupling through electron bifurcation from diverse electron donors such as formate or H2 via formate dehydrogenase or Hdr-associated hydrogenase (Kaster et al., 2011). Similar to other Rhodovulum species, the genome of strain P5 also contains all of the genes required for photosynthesis gene clusters (PGCs), including those involved in biosynthesis of bacteriochlorophyll a, carotenoids, light-harvesting systems, and reaction center components (Supplementary Table S7).

The genomes of P5 contain one large gene set that encodes a flagellum system and genes that encode methyl-accepting chemotaxis proteins, indicating that P5 may use them to facilitate their movement toward nutrient rich zones. Genes encoding gas vesicle proteins have been identified, providing buoyancy to cells as flotation devices in order to obtain the optimum amount of light and nutrients at a suitable depth in the environment. Several chemotaxis proteins were located near the gas vesicle proteins (gvp) gene cluster. In contrast, two Rhodovulum species lack the gvp gene cluster (Figure 9).

Strain P5 carries several prophage-like elements, in which a Mu-like head group phage has been induced successfully (Lin et al., 2016). The presence of clustered regularly interspaced palindromic repeat (CRISPR) arrays and their associated Cas genes in P5 and two other Rhodovulum species form a system that is possibly involved in bacterial defense against phages or plasmids. However, CRISPR sequences revealed no similarity between strains with regards to the numbers of repeats and spacer sequences, indicating that the histories of phage infection are different. Both bacteria have a predicted gene transfer agent (Supplementary Table S7).

With respect to morphological and biochemical traits as well as phylogenetic relationships and genomic analysis, isolate P5 is clearly a strain belonging to the genus Rhodovulum.


In this study, few 16S rRNA gene sequences were related to Thiomicrospira- and Nautiliales-like chemolithoautotrophic bacteria, both of which were previously reported to be abundant in the investigated Kueishantao shallow-sea hydrothermal system (Zhang et al., 2012; Tang et al., 2013). In contrast, heterotrophic bacteria in the investigated hydrothermal system were abundant. This hints at substantial spatial and/or temporal variability in the composition of the microbial communities of the shallow-sea hydrothermal system could be enormous. The genetic potential of the microbial community was analyzed using marker genes for the carbon, nitrogen, and sulfur metabolism. In the present investigation, we detected the potential for chemotrophic CO2 fixation mainly through the CBB cycle. Sulfur oxidation and reduction marker genes were present. Genes encoding enzymes involved in the denitrification pathway and H2 utilization were also detected. The metagenomes contained abundant genes responsible for heterotrophic utilization of organic substrates including transporters, glycoside hydrolases, and peptidases genes. Cultivation attempts targeting heterotrophs resulted in the isolation of 408 heterotrophic strains that are typically considered as oligotrophic, copiotrophic, or phototrophic bacteria, in which a novel species Rhodovulum sp. P5 was isolated. Physiological analysis indicated that P5 is capable of heterotrophy, autotrophy, and phototrophy. The P5 genome possessed the complete CBB cycle, PGCs, Ni-Fe hydrogenases, and the complete repertoire of genes involved in the oxidation of reduced sulfur compounds that are putative metabolic potentials of heterotrophs associated with the adaptions to the shallow-sea hydrothermal system. The genome harbored two gene sets encoding type IV secretion machinery systems, the gvp gene cluster, and the genes encoding heterodisulfide reductase that were not present in two Rhodovulum isolates from other habitats, which increased niche-specialization competitiveness in the shallow-sea environments. This study unveiled that the shallow-sea hydrothermal system harbored diverse microbial communities and their potential functions, enabling us to conduct more focused studies on heterotrophic activity in situ and ecophysiological features of isolated bacteria to ultimately get a complete picture on the tight coupling between microbes and biogeochemical cycling in the shallow-sea hydrothermal ecosystems in the future.

Author Contributions

KT and NJ conceived and designed the experiments. YZ, DL, YH, C-TC, DW, Y-SL, JS, and QZ conducted the experiments. KT and YZ analyzed the data. All of the authors assisted in writing the manuscript, discussed the results, and commented on the manuscript.


This study was supported by the National Key Research and Development Program of China (2016YFA0601100), the National Programme on Global Change and Air-Sea Interaction (GASI-03-01-02-05), the National Natural Science Foundation of China project (41276131, 41676070, 41776167, and 91751207), and the Fundamental Research Funds for the Central Universities (20720150078). 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 the Marine Research Station, Institute of Cellular and Organismic Biology, Academia Sinica for providing lab space for sample processing.

Supplementary Material

The Supplementary Material for this article can be found online at:


  1. ^
  2. ^


Alikhan, N. F., Petty, N. K., Ben Zakour, N. L., and Beatson, S. A. (2011). BLAST ring image generator (BRIG): simple prokaryote genome comparisons. BMC Genomics 12:402. doi: 10.1186/1471-2164-12-402

PubMed Abstract | CrossRef Full Text | Google Scholar

Bayraktarov, E., Price, R. E., Ferdelman, T. G., and Finster, K. (2013). The pH and pCO2 dependence of sulfate reduction in shallow-sea hydrothermal CO2-venting sediments (Milos Island, Greece). Front. Microbiol. 4:111. doi: 10.3389/fmicb.2013.00111

PubMed Abstract | CrossRef Full Text | Google Scholar

Bennett, S. A., Statham, P. J., Green, D. R. H., Le Bris, N., Mcdermott, J. M., Prado, F., et al. (2011). Dissolved and particulate organic carbon in hydrothermal plumes from the East Pacific Rise, 9°50′N. Deep Sea Res. Part I Oceanogr. 58, 922–931. doi: 10.1016/j.dsr.2011.06.010

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Brinkhoff, T., Sievert, S. M., Kuever, J., and Muyzer, G. (1999). Distribution and diversity of sulfur-oxidizing Thiomicrospira spp. at a shallow-water hydrothermal vent in the Aegean Sea (Milos, Greece). Appl. Environ. Microbiol. 65, 3843–3849.

PubMed Abstract | Google Scholar

Buchan, A., LeCleir, G. R., Gulvik, C. A., and González, J. M. (2014). Master recyclers: features and functions of bacteria associated with phytoplankton blooms. Nat. Rev. Microbiol. 12, 686–698. doi: 10.1038/nrmicro3326

PubMed Abstract | CrossRef Full Text | Google Scholar

Campbell, B. J., Smith, J. L., Hanson, T. E., Klotz, M. G., Stein, L. Y., Lee, C. K., et al. (2009). Adaptations to submarine hydrothermal environments exemplified by the genome of Nautilia profundicola. PLoS Genet. 5:e1000362. doi: 10.1371/journal.pgen.1000362

PubMed Abstract | CrossRef Full Text | Google Scholar

Canfield, D. E., Glazer, A. N., and Falkowski, P. G. (2010). The evolution and future of Earth’s nitrogen cycle. Science 330, 192–196. doi: 10.1126/science.1186120

PubMed Abstract | CrossRef Full Text | Google Scholar

Chan, B. K., Wang, T. W., Chen, P. C., Lin, C. W., Chan, T. Y., and Tsang, L. M. (2016). Community structure of macrobiota and environmental parameters in shallow water hydrothermal vents off Kueishan Island, Taiwan. PLoS One 11:e0148675. doi: 10.1371/journal.pone.0148675

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Chen, X., Zhan, P., Han, Y., Lin, D., Sun, J., and Tang, K. (2017). Complete genome sequence of Marivivens sp. JLT3646, a potential aromatic compound degrader. Mar. Genomics (in press). doi: 10.1016/j.margen.2017.06.002

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Cho, J. C., Stapels, M. D., Morris, R. M., Vergin, K. L., Schwalbach, M. S., Givan, S. A., et al. (2007). Polyphyletic photosynthetic reaction centre genes in oligotrophic marine Gammaproteobacteria. Environ. Microbiol. 9, 1456–1463. doi: 10.1111/j.1462-2920.2007.01264.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Chun, J., Lee, J. H., Jung, Y., Kim, M., Kim, S., Kim, B. K., et al. (2007). EzTaxon: a web-based tool for the identification of prokaryotes based on 16S ribosomal RNA gene sequences. Int. J. Syst. Evol. Microbiol. 57, 2259–2261. doi: 10.1099/ijs.0.64915-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Church, M. J. (2009). The trophic tapestry of the sea. Proc. Natl. Acad. Sci. U.S.A. 106, 15519–15520. doi: 10.1073/pnas.0908881106

PubMed Abstract | CrossRef Full Text | Google Scholar

Corliss, J. B., Dymond, J., Gordon, L. I., Edmond, J. M., Von Herzen, R. I., Ballard, R. D., et al. (1979). Submarine thermal sprirngs on the galapagos rift. Science 203, 1073–1083. doi: 10.1126/science.203.4385.1073

PubMed Abstract | CrossRef Full Text | Google Scholar

Embley, T. M. (1991). The linear PCR reaction: a simple and robust method for sequencing amplified rRNA genes. Lett. Appl. Microbiol. 13, 171–174. doi: 10.1111/j.1472-765X.1991.tb00600.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Giovannelli, D., d’Errico, G., Manini, E., Yakimov, M., and Vetriani, C. (2013). Diversity and phylogenetic analyses of bacteria from a shallow-water hydrothermal vent in Milos island (Greece). Front. Microbiol. 4:184. doi: 10.3389/fmicb.2013.00184

PubMed Abstract | CrossRef Full Text | Google Scholar

Grzymski, J. J., Murray, A. E., Campbell, B. J., Kaplarevic, M., Gao, G. R., Lee, C., et al. (2008). Metagenome analysis of an extreme microbial symbiosis reveals eurythermal adaptation and metabolic flexibility. Proc. Natl. Acad. Sci. U.S.A. 105, 17516–17521. doi: 10.1073/pnas.0802782105

PubMed Abstract | CrossRef Full Text | Google Scholar

Haggerty, J. M., and Dinsdale, E. A. (2017). Distinct biogeographical patterns of marine bacterial taxonomy and functional genes. Glob. Ecol. Biogeogr. 26, 177–190. doi: 10.1111/geb.12528

CrossRef Full Text | Google Scholar

Hammer,Ø., Harper, D. A. T., and Ryan, P. D. (2001). PAST: paleontological statistics software package for education and data analysis. Palaeontol. Electronica 4, 1–9.

Google Scholar

Han, Y., Lin, D., Yu, L., Chen, X., Sun, J., and Tang, K. (2017). Complete genome sequence of Serinicoccus sp. JLT9, an actinomycete isolated from the shallow-sea hydrothermal system. Mar. Genomics 32, 19–21. doi: 10.1016/j.margen.2016.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasan, N. A., Grim, C. J., Grim, C. J., Lipp, E. K., Lipp, E. K., Rivera, I. N. G., et al. (2015). Deep-sea hydrothermal vent bacteria related to human pathogenic Vibrio species. Proc. Natl. Acad. Sci. U.S.A. 112, E2813–E2819. doi: 10.1073/pnas.1503928112

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirayama, H., Sunamura, M., Takai, K., Nunoura, T., Noguchi, T., Oida, H., et al. (2007). Culture-dependent and -independent characterization of microbial communities associated with a shallow submarine hydrothermal system occurring within a coral reef off Taketomi Island, Japan. Appl. Environ. Microbiol. 73, 7642–7656. doi: 10.1128/AEM.01258-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Hügler, M., and Sievert, S. M. (2011). Beyond the calvin cycle: autotrophic carbon fixation in the ocean. Annu. Rev. Mar. Sci. 3, 261–289. doi: 10.1146/annurev-marine-120709-142712

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiao, N., Zhang, Y., Zeng, Y., Hong, N., Liu, R., Chen, F., et al. (2007). Distinct distribution pattern of abundance and diversity of aerobic anoxygenic phototrophic bacteria in the global ocean. Environ. Microbiol. 9, 3091–3099. doi: 10.1111/j.1462-2920.2007.01419.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Joint, I., Mühling, M., and Querellou, J. (2010). Culturing marine bacteria – an essential prerequisite for biodiscovery. Microb. Biotechnol. 3, 564–575. doi: 10.1111/j.1751-7915.2010.00188.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaster, A. K., Moll, J., Parey, K., and Thauer, R. K. (2011). Coupling of ferredoxin and heterodisulfide reduction via electron bifurcation in hydrogenotrophic methanogenic archaea. Proc. Natl. Acad. Sci. U.S.A. 108, 2981–2986. doi: 10.1073/pnas.1016761108

PubMed Abstract | CrossRef Full Text | Google Scholar

Kirchman, D. L. (2015). Growth rates of microbes in the oceans. Ann. Rev. Mar. Sci. 8, 285–309. doi: 10.1146/annurev-marine-122414-033938

PubMed Abstract | CrossRef Full Text | Google Scholar

Kompantseva, E. I., Komova, A. V., and Kostrikina, N. A. (2010). Rhodovulum steppense sp. nov., an obligately haloalkaliphilic purple nonsulfur bacterium widespread in saline soda lakes of Central Asia. Int. J. Syst. Evol. Microbiol. 60, 1210–1214. doi: 10.1099/ijs.0.014639-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Koren, S., Schatz, M. C., Walenz, B. P., Martin, J., Howard, J. T., Ganapathy, G., et al. (2012). Hybrid error correction and de novo assembly of single-molecule sequencing reads. Nat. Biotechnol. 30, 693–700. doi: 10.1038/nbt.2280

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, D., Tang, K., Han, Y., Li, C., and Chen, X. (2016). Genome sequence of an inducible phage in Rhodovulum sp. P5 isolated from the shallow-sea hydrothermal system. Mar. Genomics 30, 93–95. doi: 10.1016/j.margen.2016.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Manini, E., Luna, G. M., Corinaldesi, C., Zeppilli, D., Bortoluzzi, G., Caramanna, G., et al. (2008). Prokaryote diversity and virus abundance in shallow hydrothermal vents of the Mediterranean Sea (Panarea Island) and the Pacific Ocean (north Sulawesi-Indonesia). Microb. Ecol. 55, 626–639. doi: 10.1007/s00248-007-9306-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Marteinsson, V. T., Hauksdóttir, S., Hobel, C. F., Kristmannsdóttir, H., Hreggvidsson, G. O., and Kristjánsson, J. K. (2001). Phylogenetic diversity analysis of subterranean hot springs in Iceland. Appl. Environ. Microbiol. 67, 4242–4248. doi: 10.1128/AEM.67.9.4242-4248.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Masuda, S., Yoshida, M., Nagashima, K. V., Shimada, K., and Matsuura, K. (1999). A new cytochrome subunit bound to the photosynthetic reaction center in the purple bacterium, Rhodovulum sulfidophilum. J. Biol. Chem. 274, 10795–10801. doi: 10.1074/jbc.274.16.10795

PubMed Abstract | CrossRef Full Text | Google Scholar

Maugeri, T. L., Gugliandolo, C., and Lentini, V. (2013a). Diversity of prokaryotes at a shallow submarine vent of Panarea Island (Italy) by high-throughput sequencing. Atti. Accad. Pelorit. Pericol. Cl. Sci. Fis. 91:A1. doi: 10.1478/AAPP.912A1

CrossRef Full Text | Google Scholar

Maugeri, T. L., Lentini, V., Spanò, A., and Gugliandolo, C. (2013b). Abundance and diversity of picocyanobacteria in shallow hydrothermal vents of Panarea Island (Italy). Geomicrobiol. J. 30, 93–99. doi: 10.1080/01490451.2011.653088

CrossRef Full Text | Google Scholar

Maugeri, T. L., Lentini, V., Gugliandolo, C., Cousin, S., and Stackebrandt, E. (2010). Microbial diversity at a hot, shallow-sea hydrothermal vent in the Southern Tyrrhenian Sea (Italy). Geomicrobiol. J. 27, 380–390. doi: 10.1080/01490450903451518

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

McDermott, J. M., Seewald, J. S., German, C. R., and Sylva, S. P. (2015). Pathways for abiotic organic synthesis at submarine hydrothermal fields. Proc. Natl. Acad. Sci. U.S.A. 112, 7668–7672. doi: 10.1073/pnas.1506295112

PubMed Abstract | CrossRef Full Text | Google Scholar

Meier, D. V., Bach, W., Girguis, P. R., Gruber-Vodicka, H. R., Reeves, E. P., Richter, M., et al. (2016). Heterotrophic Proteobacteria in the vicinity of diffuse hydrothermal venting. Environ. Microbiol. 18, 4348–4368. doi: 10.1111/1462-2920.13304

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyer, F., Paarmann, D., D’Souza, M., Olson, R., Glass, E., Kubal, M., et al. (2008). The metagenomics RAST server– a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics 9:386. doi: 10.1186/1471-2105-9-386

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyer-Dombard, D. R., Amend, J. P., and Osburn, M. R. (2013). Microbial diversity and potential for arsenic and iron biogeochemical cycling at an arsenic rich, shallow-sea hydrothermal vent (Tutum Bay, Papua New Guinea. Chem. Geol. 348, 37–47. doi: 10.1016/j.chemgeo.2012.02.024

CrossRef Full Text | Google Scholar

Meyer-Dombard, D. R., Price, R. E., Pichler, T., and Amend, J. P. (2012). Prokaryotic populations in arsenic-rich shallow-sea hydrothermal sediments of Ambitle Island, Papua New Guinea. Geomicrobiol. J. 29, 1–17. doi: 10.1080/01490451.2010.520073

CrossRef Full Text | Google Scholar

Mohandass, C., Rajasabapathy, R., Ravindran, C., Colaco, A., Santos, R. S., and Meena, R. M. (2012). Bacterial diversity and their adaptations in the shallow water hydrothermal vent at D. Joao de Castro Seamount (DJCS), Azores, Portugal. Cah. Biol. Mar. 53, 65–76.

Google Scholar

Nagao, N., Hirose, Y., Misawa, N., Ohtsubo, Y., and Umekage, S. (2015a). Complete genome sequence of Rhodovulum sulfidophilum DSM 2351, an extracellular nucleic acid-producing bacterium. Genome Announc. 3:e00388-15. doi: 10.1128/genomeA.00388-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagao, N., Yamamoto, J., Komatsu, H., Suzuki, H., Hirose, Y., Umekage, S., et al. (2015b). The gene transfer agent-like particle of the marine phototrophic bacterium Rhodovulum sulfidophilum. Biochem. Biophys. Rep. 4, 369–374. doi: 10.1016/j.bbrep.2015.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Petersen, J. M., Zielinski, F. U., Pape, T., Seifert, R., Moraru, C., Amann, R., et al. (2011). Hydrogen is an energy source for hydrothermal vent symbioses. Nature 476, 176–180. doi: 10.1038/nature10325

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, R. E., Lesniewski, R., Nitzsche, K. S., Meyerdierks, A., Saltikov, C., Pichler, T., et al. (2013). Archaeal and bacterial diversity in an arsenic-rich shallow-Sea hydrothermal system undergoing phase separation. Front. Microbiol. 4:158. doi: 10.3389/fmicb.2013.00158

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

Google Scholar

Rusch, D. B., Halpern, A. L., Sutton, G., Heidelberg, K. B., Williamson, S., Yooseph, S., et al. (2007). The sorcerer II global ocean sampling expedition: northwest Atlantic through eastern tropical pacific. PLoS Biol. 5:e77. doi: 10.1371/journal.pbio.0050077

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Shaffer, J. P. (1995). Multiple hypothesis testing. Annu. Rev. Inc. 46, 561–584.

Google Scholar

Sievert, S. M., Brinkhoff, T., Muyzer, G., Ziebis, W., and Kuever, J. (1999). Spatial heterogeneity of bacterial populations along an environmental gradient at a shallow submarine hydrothermal vent near Milos Island (Greece). Appl. Environ. Microbiol. 65, 3834–3842. doi: 10.1128/AEM.66.7.3102-3109.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Sievert, S. M., Kuever, J., and Muyzer, G. (2000). Identification of 16S ribosomal DNA-defined bacterial populations at a shallow submarine hydrothermal vent near Milos Island (Greece). Appl. Environ. Microbiol. 66, 3102–3109.

PubMed Abstract | Google Scholar

Srinivas, T. N., Kumar, P. A., Sasikala, C. H., and Ramana, C. H. V. (2007). Rhodovulum imhoffii sp. nov. Int. J. Syst. Evol. Microbiol. 57, 228–232. doi: 10.1099/ijs.0.64470-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Sylvan, J. B., Toner, B. M., and Edwards, K. J. (2012). Life and death of deep-sea vents: bacterial diversity and ecosystem succession on inactive hydrothermal sulfides. mBio 3:e00279-11. doi: 10.1128/mBio.00279-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, K. (2014). “Microbial communities in a shallow-sea hydrothermal system,” in Encyclopedia of Metagenomics, ed. E. N. Karen (Berlin: Springer), 1–8.

Google Scholar

Tang, K., Lin, Y., Han, Y., and Jiao, N. (2017). Characterization of potential polysaccharide utilization systems in the marine Bacteroidetes Gramella Flava JLT2011 using a multi-omics approach. Front. Microbiol. 8:220. doi: 10.3389/fmicb.2017.00220

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, K., Liu, K., Jiao, N., Zhang, Y., and Chen, C. T. A. (2013). Functional metagenomic investigations of microbial communities in a shallow-sea hydrothermal system. PLoS One 8:e72958. doi: 10.1371/journal.pone.0072958

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Teira, E., Reinthaler, T., Pernthaler, A., Pernthaler, J., and Herndl, G. J. (2004). Combining catalyzed reporter deposition-fluorescence in situ hybridization and microautoradiography to detect substrate utilization by bacteria and archaea in the Deep Ocean. Appl. Environ. Microbiol. 70, 4411–4414. doi: 10.1128/AEM.70.7.4411

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilson, K. H., Blitchington, R. B., and Greene, R. C. (1990). Amplification of bacterial 16S ribosomal DNA with polymerase chain reaction. J. Clin. Microbiol. 28, 1942–1946.

Google Scholar

Yamamoto, M., and 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, L., Hong, H., Guo, W., Chen, C. T. A., Pan, P. I., and Feng, C. C. (2012). Absorption and fluorescence of dissolved organic matter in submarine hydrothermal vents off NE Taiwan. Mar. Chem. 12, 64–71. doi: 10.1016/j.marchem.2011.10.003

CrossRef Full Text | Google Scholar

Zhan, P., Tang, K., Chen, X. F., and Yu, L. M. (2017). Complete genome sequence of Maribacter sp. T28, a polysaccharide-degrading marine flavobacteria. J. Biotechnol. 259, 1–5. doi: 10.1016/j.jbiotec.2017.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Zhao, Z., Chen, C. T. A., 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: shallow-sea hydrothermal system, metagenomics, genomics, heterotrophic bacteria, microbial community, Rhodovulum

Citation: Tang K, Zhang Y, Lin D, Han Y, Chen C-TA, Wang D, Lin Y-S, Sun J, Zheng Q and Jiao N (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

Received: 26 April 2017; Accepted: 07 February 2018;
Published: 22 February 2018.

Edited by:

Stefan M. Sievert, Woods Hole Oceanographic Institution, United States

Reviewed by:

Jinjun Kan, Stroud Water Research Center, United States
Qi-Long Qin, Shandong University, China

Copyright © 2018 Tang, Zhang, Lin, Han, Chen, Wang, Lin, Sun, Zheng and Jiao. 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 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: Kai Tang, Nianzhi Jiao,

These authors have contributed equally to this work.